Skip to main content

Transcriptome of Sphaerospora molnari (Cnidaria, Myxosporea) blood stages provides proteolytic arsenal as potential therapeutic targets against sphaerosporosis in common carp



Parasites employ proteases to evade host immune systems, feed and replicate and are often the target of anti-parasite strategies to disrupt these interactions. Myxozoans are obligate cnidarian parasites, alternating between invertebrate and fish hosts. Their genes are highly divergent from other metazoans, and available genomic and transcriptomic datasets are limited. Some myxozoans are important aquaculture pathogens such as Sphaerospora molnari replicating in the blood of farmed carp before reaching the gills for sporogenesis and transmission. Proliferative stages cause a massive systemic lymphocyte response and the disruption of the gill epithelia by spore-forming stages leads to respiratory problems and mortalities. In the absence of a S. molnari genome, we utilized a de novo approach to assemble the first transcriptome of proliferative myxozoan stages to identify S. molnari proteases that are upregulated during the first stages of infection when the parasite multiplies massively, rather than in late spore-forming plasmodia. Furthermore, a subset of orthologs was used to characterize 3D structures and putative druggable targets.


An assembled and host filtered transcriptome containing 9436 proteins, mapping to 29,560 contigs was mined for protease virulence factors and revealed that cysteine proteases were most common (38%), at a higher percentage than other myxozoans or cnidarians (25–30%). Two cathepsin Ls that were found upregulated in spore-forming stages with a presenilin like aspartic protease and a dipeptidyl peptidase. We also identified downregulated proteases in the spore-forming development when compared with proliferative stages including an astacin metallopeptidase and lipases (qPCR). In total, 235 transcripts were identified as putative proteases using a MEROPS database. In silico analysis of highly transcribed cathepsins revealed potential drug targets within this data set that should be prioritised for development.


In silico surveys for proteins are essential in drug discovery and understanding host-parasite interactions in non-model systems. The present study of S. molnari’s protease arsenal reveals previously unknown proteases potentially used for host exploitation and immune evasion. The pioneering dataset serves as a model for myxozoan virulence research, which is of particular importance as myxozoan diseases have recently been shown to emerge and expand geographically, due to climate change.


The relationship between parasites and their hosts is under constant pressure; parasites must invade, replicate and feed whilst avoiding the host immune system. Proteases are the weapon of choice for parasites to overcome these challenges within the host, and can be specifically adapted for cleaving host proteins or modifying their own proteins for immune system avoidance [1,2,3]. Proteases are often high priority proteins for investigation as they have essential roles in development, invasion or feeding [4]. However, proteases are involved in other cellular functions e.g. transport and activation of other peptidases, and it can often be unclear which peptidases are essential to parasite survival or success [5, 6]. Drug or interference targets can be difficult to identify in a wash of uncharacterized proteins, however proteases linked to an essential cellular pathway or localised to a particular organelle e.g. lysosome can be considered useful targets for life cycle or development disruption [5, 7].

Anti-parasite drugs currently available have been identified by screening sets of compounds in vitro culture systems and by borrowing compounds that have worked in other pathogens and applying those to a new parasite model [8]. Firstly, this limits progress to organisms and life stages that can be isolated and cultured; secondly it relies on applicable compounds having been found in related organisms; and thirdly it limits discovery as it looks at one target at a time for feasibility. In silico drug target discovery in contrast has the attractive attributes of speed, low cost and no requirement for living parasites. In the case of non-model organisms this is likely the first step before prioritising any protein for further experimentation with the aim of anti-parasite treatment development.

Myxozoans are parasitic cnidarians that are important pathogens to both wild and cultured fish populations and yet there are no drug targets specified for this group and limited proteolytic studies to examine activity or function of selected proteins [9, 10]. Myxozoans are suggested to have reduced genomes compared to their free living cnidarian relatives [11, 12] which could have an impact on the range and diversity of the peptidases expressed. Many aspects of myxozoan biology are still unknown or inferred by comparison with other parasites to infer biology such as their metabolism (Thelohanellus kitauei - [12]), their replication [13] or proteins interacting with the host immune system (reviewed in [14]). Myxozoans are entirely parasitic in their life cycle, they alternate between a vertebrate and an invertebrate host with two entirely different types of transmissible spores in each developmental phase [15,16,17,18]. Myxospores are often hardy stages that are capable of being exposed to the environment for long periods of time waiting for uptake by their invertebrate hosts. The actinospores are generally more fragile and only viable for a limited period of time as they are released into the water column to encounter a suitable vertebrate host [19]. There are two main sources of material for genomic and transcriptomic analysis, plasmodia or cysts of developing myxospores from the vertebrate [11, 12] or actinospores released from their invertebrate host [11]. Spore development represents the final step prior to transmission with the genetic arsenal related to their production of durable spores often expressed in cysts, separated from the host immune response by connective tissue, while actinospores are collected from the environment, prior to infecting their vertebrate hosts. Therefore, they do not provide many insights into what proteins are helping the parasite feed or replicate or evade immune detection.

Sphaerosporids are a major clade of the Myxosporea, with a large proportion found in bony and cartilaginous fish, and amphibians [20,21,22,23,24]. A specific trait that has only been identified in this clade is the presence of large, extracellular stages circulating in the blood stream of their fish hosts [25,26,27]. The parasites not only use the blood for transport to their target organ but proliferate within it and are present almost all year round (Fig. 1, [26, 28, 30]). Sphaerospora molnari is a parasite of the common carp in Central Europe with motile blood stages that provoke a strong immune response [29] and are a likely co-factor for developing Swim Bladder Inflammation [30]. S. molnari blood stages (SMBS) are prime targets for parasite intervention therapy, as they are 1) responsible for massive proliferation in the earliest stages of infection of fish, 2) freely circulating in the blood and any drug targeting the SMBS would not need to be applied to host tissue or taken up by host cells; 3) they are circulating in the blood for an extended period and therefore there is a longer window for application of anti-parasite therapies. In addition, preliminary protein studies on SMBS show a high level of sequence divergence even in highly conserved proteins such as actin [28] and therefore SMBS could potentially have proteases that are highly divergent from their hosts as well as other cnidarians which would aid protein target assay development. This study examines protease families and groups present in the transcriptome of SMBS to investigate their diversity and divergence. We compare key protease groups with examples known from other parasites that have been successfully flagged as drug or anti-parasite targets. In addition, we provide gene expression data for selected candidates with the goal of identifying stage-specific proteases of interest for future functional studies.

Fig. 1

Developmental cycle of Sphaerospora molnari within its host Cyprinus carpio. Sphaerospora molnari blood stages and infected gill images by A.S. Holzer, (blood stages modified from [28, 29]). Common carp (host) royalty free stock image (


This transcriptome is the first next generation sequencing dataset from any sphaerosporid species, and also the first dataset from a highly proliferative, extrasporogonic developmental myxozoan stage. Pooled host (Cyprinus carpio) blood cells and Sphaerospora molnari blood stages (SMBS) from 8 infected fish were used for this transcriptome. Illumina HiSeq sequencing yielded a total of 52,040,447 clean paired reads, mapping these to the gene models of C. carpio removed 14,849,448 reads (BioProject PRJNA522909). A trinity assembly of these remaining reads gave 157,506 transcripts, (mean length 766 nt, 39.75% GC content). 127,741 of these transcripts (81.1%) were found in the carp genome with an e-value of 1e− 05 or more and a percentage identity of > 75%. The remaining 29,765 transcripts were therefore assumed to be S. molnari, and were translated into 29,588 proteins (Table 1). To examine the presence of potential chimeric sequences, we amplified a substantial number (n = 15) of our flagged proteases and ribosomal DNA to verify their sequence according to our assembly. Sanger sequencing of the complete ribosomal DNA yielded 13,486 bp (Genbank Acc. Nr MK533682), and large fragments of all flagged proteases were also verified from cDNA (Supplementary Data).

Table 1 Summary of S. molnari blood stages (SMBS) transcriptome dataset

Screening with BUSCO, identified that S. molnari retained and is expressing at least half of the 978 benchmark metazoan genes [31] (53% of the single copy metazoan genes). We completed the same analysis for myxozoans Myxobolus cerebralis and Kudoa iwatai, and non-myxozoans Polypodium hydriforme, Edwardsiella lineata, Nematostella vectensis (Table 2). S. molnari had the highest number of complete BUSCOs of all the myxozoan datasets, whereas N. vectensis had the highest overall (905/978, 93%), similar results were found for the single copy BUSCOs. S. molnari had the lowest number of missing genes (374/978, 38.2%) within the myxozoans, in comparison N. vectensis was only missing 30 genes (30/978, 3%) (Table 2).

Table 2 Benchmarking Universal Single Copy Orthologs (BUSCO) identified in datasets

We queried the transcriptomic dataset of SMBS for proteases using representative protease sequences downloaded from the MEROPS database. Less than 1% of all transcripts had a strong sequence match. There were 235 homologs identified in SMBS representing 45 peptidase families, the majority of the proteases were cysteine (38%), followed by metalloproteases (31%), serine (15%), threonine (14%) and aspartic groups (2%) (Fig. 2). Families that were highly represented in SMBS (Table 3) were C01 Papain-like proteases, C12 and C19 – ubiquitinyl hydrolases, M24 – aminopeptidases and dipeptidyl peptidases, M67 – ubiquitin releasing proteases often associated with proteasomal degradation and T01 – proteasome proteases. There were many families that were absent in all examined cnidarians, and even more, missing from only the myxozoans compared to the free living species e.g. S28 lysosomal carboxypeptidase.

Fig. 2

Comparative summary of expressed proteases identified in three myxozoan datasets, two non-myxozoan parasitic cnidarian datasets and a free living cnidarian species. Pie charts showing diversity and abundance of protease clans and families in transcriptomes of Sphaerospora molnari and other myxozoans (Myxobolus cerebralis and Kudoa iwatai), two non-myxozoan parasitic cnidarians (Edwardsiella lineata and Polypodium hydriforme) and free living species Nematostella vectensis

Table 3 Protease families identified in S. molnari bood stages (SMBS) and other cnidarian species

To more closely examine the proteases present in SMBS, we looked at enzymes that were in highly represented families, and were transcripts with a high number of reads mapping to them (TPM) or had high similarity to proteases known from other parasite species. In particular, we examined cysteine proteases in the MEROPS family C01 – cathepsin L, aspartic proteases in the family A22 – Signal peptide and Presenilin-like proteases, metallopeptidases in the M12 – the metzincins, and S09 the Prolyl oligopeptidase family.

Cathepsins: S. molnari’s transcriptome revealed eight cathepsin-L like sequences by sequence homology, however, five were excluded from our further analysis due to 1) incomplete transcript and 2) missing or uncharacterised active sites at either substrate S1 or S2 sites, or 3) a sequence homology that appeared to be closest to cathepsin L but in fact had chymotrypsin-like folds (Ser-His-Asn) or Gly-His-Asn catalytic triads. Three cathepsin Ls analysed (CathL1–3) were all propeptides with signal peptides which may indicate later activation within the cell (e.g. lysosome) or extracellularly [32]. All had conserved a glutamine for the oxyanion hole known for cysteine peptidases, however, Sm_CL3 did not retain a stabilising asparagine close to the His active site, this was replaced by a negatively charged aspartic acid. We aligned the predicted tertiary structures of S. molnari cathepsin Ls to a procathepsin L from a metazoan parasite with x-ray crystallography evidence (PDB: 2O6X, [33]) (Fig. 3a-c). Two of S. molnari’s cathepsin Ls, Sm_CL 1 and 2 were able to be aligned to the crystal structure with high confidence, however, there was a marked difference in the distribution of hydrophobic resides (Fig. 3b-c). In particular, there were higher numbers of hydrophilic residues at the active site compared to Fasciola hepatica (Fig. 3d-f). The number of charged residues were similar overall, however the distribution of positive and negatively charged amino acids was different between all three proteins (Fig. 3g-i).

Fig. 3

Predicted tertiary structure of SMBS Cathepsin L 1 and 2 and comparison with Fasciola hepatica cathepsin L. aFasciola hepatica cathepsin L crystalised structure 2O6X showing hydrophobic residues (orange). b Sm_CL 1 predicted structure based on Phyre2 model aligned to F. hepatica pdb. c Sm_CL2 predicted structure based on Phyre2 model aligned to F. hepatica pdb. d-f Closer view of substrate binding site 2, D - F. hepatica cathepsin, E – Sm_CL1, F – Sm_CL2. g-i Charged residues white = neutral, red = negative, blue = positive of G - fasciola hepatica, H - Sm_CL1, I - Sm_CL2

Aspartic: We identified two aspartic proteases in the transcriptome of SMBS, one presenilin-like and a signal peptide peptidase-like sequence (Sm_SPPL). The presenilin-like protease had a lower TPM value than Sm_SPPL (12.66272 compared to 146.829). We therefore focused our analysis on Sm_SPPL, which clustered with other invertebrates (Fig. 4a) and conformed to the structure of a signal peptide peptidase including active residues in transmembrane domains 6, 7 and 9 in the predicted model (Fig. 4b, c, d). Sm_SPPL has the appropriate active site residues and was able to align the predicted tertiary structure of a SPPL. Sm_SPPL was distinguished from a presenilin by the presence of QPALLY motif in its last transmembrane domain, like others in this group [34].

Fig. 4

Aspartic protease of S. molnari – signal peptide peptidase. a Phylogenetic relationship of S. molnari SPPL (Sm_SPPL), (ML, 1000 bootstrap support values shown at nodes). b Predicted tertiary structure of Sm_SPPL according to Phyre2, showing nine transmembrane domains with active site highlighted in pink. c Closer view of active site and interacting residues. d Schematic of Sm_SPPL structure showing location of active sites and motifs within sequence

Metallopeptidases: Seventy two metallopeptidases in 16 of 29 families were flagged in the S. molnari proliferative stage peptidases. All of the S. molnari metallopeptidase families were shared by one or both of the other myxozoan species (M. cerebralis or K. iwatai) examined as well as two or three of the free living species (Table 3). The highest proportion of SMBS metallopeptidases were in the M12 family (adamalysins and reprolysins). We examine five of these metallopeptidases in more detail (Sm_MP1 TPM = 14.5; Sm_ MP2 TPM = 16.17; Sm_MP3 TPM = 22.07; Sm_MP4 TPM = 2.76; Sm_MP5 TPM = 0.79). All except one of the five have the methionine turn and zinc binding motif HEXXHXXGXXH with either a serine or an aspargine binding residue. Two of the target metallopeptidases had signal peptides, and almost all had reprolysin and disintegrin domains. Sm_MP4 has a potential C-terminal transmembrane tail that may anchor it similarly to other well-known ADAM proteins.

Serine: Multiple serine proteases were identified within SMBS transcriptome, the largest was MEROPS family S09. The transcriptome of SMBS yielded a serine protease with sequence and structural homology to dipeptidyl peptidase 4 (Fig. 5b, c). SMBS_DPP contained the catalytic triad (Ser-633, Asp-711, His-743), a predicted transmembrane domain at the N-terminus. Phyre2 was able to model the sequence with 100% confidence according to the models of other dipeptidyl proteases (DPPIV and VIIII). The sequence clustered to Thelohanellus kitauei DPPIV and Hydra vulgaris POP isoforms (Fig. 5a) however, both of the H. vulgaris isoforms contained a transmembrane domain and a high homology to DPPIVs rather than POPs.

Fig. 5

Phylogenetic placement and predictive structure of S. molnari dipeptidyl peptidase. a maximum likelihood phylogenetic analysis of S. molnari serine protease (MEROPS family S09) showing grouping with dipeptidyl peptidases rather than prolyl oligopeptidases with Thelohanellus kitauei. b Predictive structure of Sm_DPP as a monomer. c Overlayed predictive structures of human dipeptidyl peptidase IV dimer (pdb id: 2G5T, yellow) and Sm_DPP (blue)

We then examined the expression of eight key proteases in blood stages compared to spore forming gill stages by qPCR and also in silico expression (TPM). Three cathepsins Sm_CL1, 2 and 3; a presenilin like aspartic protease (Sm_SP1); a dipeptidyl peptidase (Sm_DDPIV); a metallopeptidase (Sm_MP1) and two lipases (Sm_Lip1 and 2) were used as target proteases. Half of our candidates were upregulated in blood stages rather than gill stages (Fig. 6). Expression of Sm_CL3 was over 500x higher in sporogonic gill stages than presporogonic SMBS (Fig. 6a). Sm_CL3 was predicted to contain a transmembrane domain and had a relatively low TPM value according to our assembly (2.3503) compared to the other two cathepsins Sm_CL1 and 2 (TPM = 782.87 and 77.28, respectively), both of which had signal peptides and were transcribed as procathepsins. Both Sm_CL1 and 2 showed similar expression in the blood and the gills. An astacin metallopeptidase and two lipases were also upregulated in blood compared to gill samples (Fig. 6b). In contrast, both the aspartic and serine proteases, Sm_SP1 and Sm_DPPIV respectively were significantly upregulated in spore forming stages in the gills. Sm_DPPIV did not have a high relative expression within our transcriptome (TPM = 3.99), however when we compared its expression in blood stages vs. gill stages, we found it was more highly expressed in sporogonic gill stages by almost 100 fold. Sm_SP1 had a higher TPM value than Sm_DPPIV (146.829) and was 42 times more highly expressed in the gills than blood.

Fig. 6

Real-time PCR of selected proteases in S. molnari blood stages (non-sporogonic) and gill stages (sporogonic). Relative abundance to two housekeeping genes (Elongation Factor 2 and Glyceraldehyde-3-phosphate dehydrogenase) in cDNA samples of circulating blood stages (n = 3) and spore forming gill stages in common carp (n = 3) including 95% confidence intervals and average for each marker. a Cathepsins (Sm_CL1, 2 and 3) in blood stages (red) and sporogonic gill stages (blue); b Other proteases (A22 – Signal peptide and Presenilin like peptidases (Sm_SP1); M12 metallopeptidase (Sm_MP1); S09 – Prolyl oligopeptidase (Sm_DPPIV 1) and Lipases (Sm_Lipase 1 and 2) in blood stages (red) and sporogonic gill stages (blue)


This transcriptome analysis of S. molnari is the first from this entire clade of sphaerosporid myxozoans, this is the only group for which an invertebrate host is unknown and is the only dataset from an extrasporogonic stage of development. Therefore, it offers a unique insight into the mechanisms of myxozoan development and host interactions. Our focus on proteases, was primarily in the pursuit of identifying proteins that would be worthy of further investigation to understand their role in the host-parasite relationship.

The transcriptome of S. molnari blood stages (SMBS) appears to follow a similar trend as other parasitic and free living cnidarians in its heavy AT richness for many of its genes [35] and its overall reduction in genes [11, 12]. Gene divergence and AT richness of the SMBS dataset aids distinction of host and parasite genes in many cases and could also aid targeted gene therapy as seen in other anti-parasite drug design e.g. Plasmodium falciparum [36]. The number of proteins identified could have been limited due to the extensive host filtering we conducted on the dataset, divergent genes and/or the nature of the particular life stage of SMBS. Separation of the SMBS from the host immune cells would significantly improve further transcriptomic analysis to combine physical and bioinformatic filtering. More than half of the known BUSCO metazoan genes were identified in this dataset, similarly to M. cerebralis triactinomyxon stages and sporogonic K. iwatai transcriptomes. Although using a different comparative dataset (CEGMA, previous benchmarking metazoan dataset) Chang et al. [11] also noted a reduced number of core eukaryotic genes in the transcriptomes and genomes of K. iwatai (70% of CEGS, compared to our 56% of BUSCOs) and M. cerebralis (39% of CEGS compared to 55% of BUSCOs) whilst P. hydriforme consistently had more than its myxozoan relatives and was closer in numbers to the free living species (90% of CEGs and 89% of BUSCOs). BUSCO genes are also used to estimate the completeness of a genomic or transcriptomic dataset, whilst it is unlikely that we sequenced the extent of the expressed transcripts of S. molnari due to overwhelming host contamination, it is interesting to note an overall reduction in this gene set for myxozoans in general. For comparison, other early branching metazoans and cnidarians have retained a high percentage of these conserved proteins, e.g. Placozoa species Hoilungia hongkongensis had 90–95.3% [37]; sponge Halisarca caerulea had 92.6% [38]. Free living cnidarians have variable accounts of their retained complete BUSCO genes, e.g. Acropora digitifera 69% (675/978) [39]; Exaiptasia pallida 84.5% (826/978) [39]; Scleractinian coral Porites lutea 56.4% [40]. Other metazoan parasites have also retained a higher level of genes e.g. the trematode Microphallus sp. maintained close to 70% of genes from the same metazoan dataset [41]. The transcriptomes of parasitic nematode Teladorsagia circumcincta were found to have only 23–38% of the conserved gene set yet appeared to have good matches to the representative genome [42]. Furthermore, cestode parasites have a range of conserved eukaryotic genes from Taenia multiceps (81.8%) to Echinococcus multilocularis (91.6%) [43]. Additionally, despite removing protein sequences that were more than 90% similar to each other in the S. molnari dataset, there were a high number of duplicated single copy BUSCO genes (Table 2). This may be caused by the pooling of wild-type individuals sampled from different locations and populations. De novo assemblers such as Trinity, designed to generate alternative transcripts, are more sensitive to sequencing errors and/or highly heterozygous datasets, and thus multiple loci may have been assembled into ‘isoforms’ of the same parent transcript increasing the amount of duplication [44, 45].

We examined some specific examples of proteases of interest within the S. molnari transcriptome that were homologous to other known parasitic proteases involved in important host-parasite functions. There are examples of parasitic proteases that have been exploited as drug targets in all different groups and protease families. Two of the aspartic peptidases we identified in the SMBS transcriptome were also in groups that had been identified as potential drug targets. A presenilin-like and a signal peptide peptidase-like protease (Sm_SPPL) were found in SMBS both are in the MEROPS protease family A22. The proteases are closely related, SPPLs are transmembrane proteases with their catalytic sites buried in bordering transmembrane domains, and have varied functions in eukaryotes depending on their localization within the cell [46, 47]. In humans, they are involved in the processing of peptides for the MHC I epitopes and their central role in processing signal peptides indicates an important role in signaling and protein modification as part of the endoplasmic reticulum [34]. In contrast to SPPLs, presenilins require a co-factor and multiple subunits to form a catalytic γ -secretase complex [46]. Both aspartic proteases have been identified as potential drug targets in Plasmodium falciparum [48], and can potentially play a role in the cleavage of peptides for antigen presentation on a cell surface [49]. SPPL blockers have been shown experimentally in parasite-derived homologs e.g. LY411,575 (reduced Plasmodium berghei in mice and humans [50]), and bioinformatically as suggested in protozoan species Toxoplasma gondii, Leishmania infantum, and Trypanosoma cruzi [51]. Both S. molnari proteases are divergent from their hosts and in the case of Sm_SPPL have been shown to be highly expressed during its sporogenesis development in the host’s gills. This localization may be useful for targeting future therapeutic assays at such an aspartic protease in this system.

Another group of proteases that are commonly targeted for anti-parasite therapy are the metallopeptidases. These have been associated with many functions of parasitic development, they are linked to digestion e.g. falcilysin in Plasmodium falciparum [52], host immune system disruption e.g. matrix and secreted metalloproteases [1, 53]; or host immune system evasion e.g. Leishmania GP63 [2]. SMBS had a large variety of metallopeptidases, it is also important to note families that are missing in the myxozoan species, for example we did not identify any M08 homologues which have been originally characterised in Leishmania to aid immune system avoidance in promastigotes [54]. Nor did we identify any matrix metallopeptidases (M10 family) which have also been previously identified as protozoan parasite drug targets [55, 56]. One of the largest groups we identified in SMBS was M12, the metzincins. These have been found in other parasitic organisms feeding in their host’s blood and are considered important for digestion or nutrient uptake which may be their role in SMBS as they are circulating in the host’s blood. The diversity of SMBS metzincins also indicates varied roles for these metallopeptidases in these developmental stages. Although the function of each specific protease are not known, and they were not considered to be highly expressed (TPM values range from 0.7–22.07), they share some valuable attributes that make them strong candidates for future research. Primarily, they share a key domain and binding site and none of the potential target sequences (Sm_MP1–5) had “cysteine rich” domains identified in their sequence which has been discussed as a possible hindrance to recombinant protein production [57]. Recombinant protein application or other interference with parasite metzincins has been shown to be a successful anti-tick treatment [58] and may potentially be employed with any of these SMBS metallopeptidases. A recombinant protein with epitopes of multiple SMBS metzincins, considering they are so divergent from each other as well as any recognisable fish host homolog, could hold potential for a therapeutic experiment.

In the serine proteases, the largest family found in SMBS was S09. This family has multiple groups of which only some have been identified as putative drug targets. One such group is the dipeptidyl peptidases (DPPs) and prolyl oligopeptidases (POPs), which specialize in the cleavage of proline residues (DPPs cleave at N-terminus, POPs cleave at C-terminal). The S09 protease we looked at as a potential drug target in SMBS had similarities to both groups however we determined it to be a DPP according to its sequence, phylogenetic position and its structural homology to DPPIV (Fig. 5). DPPIV has also been identified in bee and snake venoms and possibly influences vasodilation and constriction, something that may be useful to SMBS in host blood vessels. DPPIV exists as a homodimer (potentially tetramer) targets proline residues and has suggested roles in antigen presentation on the cell surface, cell adhesion and collagen binding activity [59]. In protozoans, DPPIVs have been shown to play a role in the encystment of Giardia, a process that can be blocked by the application of inhibitors [60]. DDPIVs have been associated and inhibited in other blood-feeding parasites including Haemonchus contortus where it is suspected to play a role in fibrogen breakdown and coagulation in the blood [61]. Due to the choice and variety of DPPIV inhibitors that are available for both human and veterinary medical applications, we propose that SMBS DPPIV could be a good target for assessing the impact of its inhibition on the proliferation and metabolism of S. molnari in its host. This is based partially on its unique sequence as evidenced from its relation to other homologs, its predicted structure which could aid the application of known inhibitors, and its expression in a developmental stage that could be feeding on blood.

One of the most well researched protease groups in parasites are cysteine proteases, commonly investigated for their potential as therapeutic targets in parasites as they are known to be involved in a variety of pathways as well as being associated with parasitic development and proliferation [62]. In myxozoans, cysteine proteases have been suggested to be involved in the proteolytic destruction of host tissue based on activity and substrate assays [10, 63, 64]. They are also associated with host hemoglobin degradation in blood feeding parasites, e.g. many parasitic helminths [65] and maybe the case for S. molnari which proliferates in host blood and interacts with erythrocytes [29]. The replacement of a stabilising asparagine close to the His active site of Sm_CL3 was also seen in a cathepsin L of the myxozoan Kudoa thyrsites. K. thrysites is suggested to use its cathepsins for the degradation of host tissue, i.e. “Milky Flesh Disease” [9]. S. molnari cathepsins appear to have all of their disulfide bridges in the three isoforms we examined in contrast to the cathepsin Ls reported from K. thyrsites [9]. Sm_CL1, 2 and 3 had close homology to Fasciola hepatica cathepsins, these proteases are expressed throughout F. hepatica’s development with distinct roles in feeding and invasion based on their binding sites [62, 66]. All Sm_CLs had higher numbers of hydrophilic residues at the active site compared to F. hepatica (Fig. 3a-c). The number of charged residues were similar to each other overall, however the distribution of positive and negatively charged amino acid was different between all three proteins (Fig. 3). These changes could potentially impact the activity and substrate affinity for the proteases. Cysteine proteases are particularly studied in parasites due to their roles in moulting, encysting and digestion across parasite taxa, and in particular they are often targeted for anti-parasite therapies or inhibition [62, 65, 67,68,69]. Here, we present the first expression analysis of cathepsins in SMBS that are potentially important in the feeding, immune evasion or tissue penetration of this parasite in its host. Characterisation of these proteases (substrate, inhibitors, abundance in the proteome) will advance our knowledge of the roles these cathepsins play in the development of this parasite and further inform our prioritization of protease targets for intervention and control assays.

Cathepsin Ls were the most highly expressed (using TPM as indication) proteases identified within our transcriptome, that combined with the large number of in vitro and in vivo experiments on other parasite cathepsins and their effect on parasite survival make them prime candidates for further development. Some examples use peptides mimicking compounds such as aziridine-2,3-dicarboxylate-based inhibitors to inhibit cathepsin L activity e.g. Trypansoma brucei rhodesiense [4] or non-peptidic inhibitors e.g. chalcones in F. hepatica [70]. S. molnari is unusual in the number of cathepsin L’s expressed in its blood stages, not just isoforms but distinct proteases. Until in vitro assays are able to better pinpoint which cathepsins are linked to SMBS survival, we rely on sequence similarity and predictive modeling to infer homology to other cathepsins including those that have been inhibited successfully. Sm_CL1 and 2 as described above, had signal peptides for potential secretion and they differed in amino acid composition (37.0% identity). Sm_CL2 had more hydrophobic residues at the S1 active site and Sm_CL1 was more highly expressed in the blood stage than Sm_CL2. Both had predicted structures that aligned to FhCL1 F. hepatica’s cathepsin’s structure which has been inhibited as well as another of its cathepsins FhCL3 with flavonoid compounds both experimentally and with computational docking [70]. Both F. hepatica cathepsins exhibited hydrophobic interactions with compound C34, and were successfully inhibited although at different rates depending on the concentration of the compound. Experimental evidence would be needed, however we propose Sm_CL1 and 2 as good candidates for inhibition with flavonoids such as C34 in future in vitro assays. The divergence in sequence identity from their host, their expression in both blood and gill stages and their similarity to proteases that have been successfully blocked in such assays are all good evidence that these could be the first drug targets for S. molnari.

Comparison of the expression of some of the proteases in non-spore forming blood stages and gill sporogonic stages gave further insights into the function of some groups. The cathepsin Ls that we examined each had a different expression profile during SMBS development, with Sm_CL3 being the main cathepsin expressed in the gills. Localisation and expression would influence which cathepsin to target for future work, and may indicate roles in sporogenesis for Sm_CL3 compared to CL1 and CL2 in the blood. CL1 and 2 may act as other cathepsins in blood feeding parasites and break down hemoglobin for sustenance. Lipases and metallopeptidases are also associated with feeding in parasites and were comparatively higher in SMBS than in the gills. Astacin metallopeptidases have also been shown to break down complex proteins such as hemoglobin or migrate through tissues [4, 71, 72], and its upregulation in the proliferative developmental stage rather than spore forming could indicate a role in nutrient acquisition or penetration of target tissues (Fig. 1) in SMBS. Lipases have been suggested to play role in parasite cell surface protein presentation as well as digestion of host protein [73, 74], the two candidates we show are upregulated in SMBS should be further investigated for their substrate and localisation within the SMBS to further characterise if they are integral to parasite survival. Disruption of such enzymes could be of interest to uncover the proteomic arsenal SMBS use to digest and navigate its host for survival. Conversely, the aspartic (Sm_SP1) and serine protease (Sm_DPPIV) would be more informative to investigate where they appear to be highly expressed – the gills. These two proteases have been associated with antigen presentation and host immune system avoidance, it could be significant that they are relatively abundant when the parasite is bound within host target tissue, rather than its extracellular blood stage, where it would be directly in conflict with the host immune system. Localisation and further characterisation will be vital to learning more about the role these proteases play in the developmental stages we identified them in and determine if they are interesting targets with anti-parasitic potential.


We produced and explored the first transcriptomic dataset of early proliferative myxozoan stages to date and identified family expansions in cysteine, metallo and threonine clans. We did not identify any myxozoan-specific radiations in particular groups, however, all of the myxozoan proteases we examined were highly divergent from each other as well as from other cnidarians. Vaccine development against a number of metazoan parasites is based on proteases as antigens of interest [75,76,77]. However, with regard to myxozoans, the function and involvement of these enzymes in host-parasite interaction first need to be elucidated as a major lack of knowledge exists with regard to metabolomics and the molecular means of host interaction. The general strategy for therapeutically targeting proteases is to identify a specific inhibitor — generally a small molecule — that blocks the active site. Discovery efforts for new inhibitors have typically been based on the structure of known protease substrates, presenting a substantial challenge for the development of peptidomimetic compounds that have the pharmacokinetic characteristics needed to be suitable as a drug.

This study advances our knowledge of myxozoan protease sequence, predicted structure and in some cases hydrophobicity and amino acid changes. This information furthers investigation into the potential role these proteases play in the development, sustenance and host immune evasion of these important parasite stages. Vaccination plays an important role in commercial large-scale fish farming and is a key reason for the success of salmon aquaculture, however, available vaccines are aimed at bacterial and viral pathogens, while parasite vaccines for fish are still inexistent [78], likely due to an insufficient knowledge of potential parasite target molecules of fish parasites, when compared with parasites of human or veterinary importance. It is now more than timely to explore new genomic data for such targets as epidemical models predict major emerging disease outbreaks and an increased geographic range of myxozoan species such as T. bryosalmonae, in relation to temperature increase [79, 80], with new records from northern European territories [81,82,83] and recent reports of massive fish killings from the Yellowstone river [84, 85], Hutchins et al. 2017). Exploring the enzymes expressed during early establishment and proliferation of myxozoan infections is essential to finding putatively relevant vaccine targets that can inhibit rapid multiplication of cryptic parasite stages in fish, long before the onset of disease. The proteases we discuss here are putative targets for further research, confirmation of their expression in different stages of S. molnari’s life cycle (SMBS vs. gill sporogonic stages vs. extracellular secretion) would be an invaluable method of testing their activity and function and therefore their use in anti-parasitic development.


Animal and sample collection for transcriptome analysis

Common Carp (Cyprinus carpio) were obtained from two localities, Štrmilov in Czech Republic (49.1644° N, 15.2031° E) and Hortobágy in Hungary (47.3542° N, 21.0000° E) during 2013–2015. All fish were obtained commercially, were less than 2 years old and were transported live to the laboratory. Fish were anaesthetized with clove oil and blood was taken, all animal procedures were performed in accordance with Czech legislation (section 29 of Protection of Animals Against Cruelty Act No. 246/1992) and approved by the Czech Ministry of Agriculture. 8 fish were found to be infected with Sphaerospora molnari blood stages (SMBS), and whole blood was centrifuged for 5 min at 3500 rpm in heparinized hematocrit tubes to isolate host white blood cells mixed with SMBS.

For qPCR analysis

Presporogonic S. molnari blood stages and sporogonic gill stages were isolated from experimentally infected fish (n = 3) in the Czech Republic or recirculation system in Hungary (n = 3). Specific parasite free fish tissues were selected from laboratory cultures within the Czech Republic (n = 3). Fish were euthanized according to the ethics license and methods above.

RNA isolation and Transcriptome assembly

Total RNA was isolated from blood stages/host white blood cell mixtures with Macherey-Nagel NucleoSpin Kit II (Biotech, Czech Republic) for transcriptome sequencing at Beijing Genomics Institute (Bejing Genomics Institute, Hong Kong) with an Illumina HiSeq 2000 (75 bp paired end reads). Reads were filtered for bacterial contaminants and then aligned to the Common carp genome ( with bowtie2, using -very-sensitive parameter. All remaining reads were then assembled with Trinity v2.4.0 and the transcripts were compared with the carp genome again for further host transcript identification (BLAST parameters: tblastx, 1e− 05), those with a percentage identity of > 75% were removed to create a “non-host” dataset which was translated to protein using OrfPredictor [86]. Redundancy of the translated non-host dataset was removed with CD-Hit using 0.9 cutoff [87] and blastp (1e− 05) annotated with NCBI nr database.


The non-redundant dataset was screened for assembly quality and completeness by identifying BUSCO genes ( using the metazoan 09 dataset [31]. Proteases and inhibitors were identified using the “meropsscan” dataset ( using blastp, 1e 05 [88]. Signal peptides and transmembrane domains were predicted using Signal P and TmHmm (

Comparative assemblies

Transcriptomic sequences for other cnidarian species were downloaded from NCBI: Myxozoans – Myxobolus cerebralis triactinomyxon stages from Tubifex host (PRJNA258474), Kudoa iwatai myxospores from cysts of Sparus aurata tissue (PRJNA248713). Non-myxozoan cnidarians with endoparasitic life stages - Polypodium hydriforme non-parasitic stolons (PRJNA251648), and Edwardsiella lineata (mix of parasitic and free living life stages) (downloaded from EdwardBase site when was active in 2014 [89]. Finally, to compare with a completely free living relative, genome derived RefSeq mRNA sequences of the anthozoan Nematostella vectensis PRJNA19965 were downloaded from NCBI. Transcripts were translated into peptide sequences by OrfFinder [86] and searched for proteases as above, the redundancy was removed (0.9 cutoff) and the non redundant dataset was screened for BUSCO genes as above.

Target protease groups: Four groups of proteases were more closely examined (A22 – Signal peptide peptidases and presenilins, C01-Cathepsin Ls, M12 – metalloendopeptidases, S09 – Prolyl oligopeptidases and dipeptidyl peptidases). Representative S. molnari sequences with other sequences from Genbank and Uniprot (including fish sequences) were analysed phylogenetically with fish and other cnidarian and parasitic-derived sequences using RaxML (L + G + I). Tertiary structures were predicted for key proteases using the Phyre2 server [90] and models were compared and manipulated in PyMol ver. 1.4.1 (The PyMOL Molecular Graphic System).

Sanger sequencing of key predicted proteases and rDNA

RNA was extracted from 3 biological replicates of S. molnari proliferative blood stages (pooled samples from several individuals) and spore-forming stages (individual fish). Total host+parasite RNA was isolated using the Nucleospin RNA Kit (Machery-Nagel) including a DNase treatment step. RNA concentration and purity was checked using a Nanodrop (ND-1000) Spectrophotometer (NanoDrop Technologies) and cDNA was synthesised using the Transcriptor High Fidelity cDNA synthesis Kit (Roche). Primers were designed to amplify full length sequences of selected proteases and ribosomal DNA as single amplicons or with long overlaps between individual sections, to confirm the assembled transcriptome sequences.

Quantification of stage-specific expression of candidate proteases

Gene-specific primers were designed to amplify short, 70–150 bp regions suitable for qPCR (Supplementary Data). All primers were tested for functionality and specificity using conventional PCR prior to performing qPCR. qPCR was performed using the FastStart Universal Sybr Green Master Mix (Rox) on LightCycler® 480 Real-Time PCR System (Roche). Reactions contained 12.5 ul of FastStart Universal SYBR Green PCR Master Mix, Roche, Germany (2X conc.), 1 μl of each forward and reverse primer (10 μM conc.), 5.5 μl of PCR grade water, and approx. 500 ng of cDNA, resulting in a final volume of 25 ul. Cycling conditions were as follows: Denaturation at 95 °C for 5 min, followed by 50 cycles of 95 °C for 10 s, and 58 °C for 10 s and 72 °C for 10 s. Melting curve analysis were performed after each qRT-PCR to ensure primer specificity. The relative expression ratio of each sample was calculated according to Pfaffl [91], based on the take-off deviation of sample versus controls at each time point and normalized relative to Elongation Factor 2 and Glyceraldehyde-3-phosphate dehydrogenase (housekeeping genes, [29]). Confidence intervals, and box plots made in R.

Availability of data and materials

Raw unfiltered sequence reads are deposited to Genbank under Bioproject PRJNA522909. Primers used for qPCR are provided as supplementary material. Host and parasite transcripts are available through Dryad depository link Sanger sequenced ribosomal DNA was submitted to Genbank under MK533682.



Sphaerospora molnari Blood Stages


Open Reading Frame


Benchmarking Universal Single Copy Orthologs


Core Eukaryotic Genes Mapping Approach




Transcripts Per Kilobase Million


Dipeptidyl peptidases


Prolyl oligopeptidases


  1. 1.

    Culley FJ, Brown A, Conroy DM, Sabroe I, Pritchard DL, Williams TJ. Eotaxin is specifically cleaved by hookworm metalloproteases preventing its action in vitro and in vivo. J Immunol. 2000;165(11):6447–53.

    CAS  PubMed  Google Scholar 

  2. 2.

    Isnard A, Shio MT, Olivier M. Impact of Leishmania metalloprotease GP63 on macrophage signaling. Front Cell Infect Microbiol. 2012;2:72.

    PubMed  PubMed Central  Google Scholar 

  3. 3.

    Schmid-Hempel P. Immune defence, parasite evasion strategies and their relevance for 'macroscopic phenomena' such as virulence. Philos Trans R Soc Lond Ser B Biol Sci. 2008;364(1513):85–98.

    Google Scholar 

  4. 4.

    McKerrow JH, Caffrey C, Kelly B, Loke P, Sajid M. Proteases in parasitic diseases. Annu Rev Pathol. 2006;1:497–536.

    CAS  PubMed  Google Scholar 

  5. 5.

    Doyle MA, Gasser RB, Woodcraft BJ, Hall RS, Ralph SA. Drug target prediction and prioritization: using orthology to predict essentiality in parasite genomes. BMC Genomics. 2010;11:222.

    PubMed  PubMed Central  Google Scholar 

  6. 6.

    Steverding D, Sexton DW, Wang X, Gehrke SS, Wagner GK, Caffrey CR. Trypanosoma brucei: Chemical evidence that cathepsin L is essential for survival and a relevant drug target. Int J Parasitol. 2012;42(5):481–8.

    CAS  PubMed  Google Scholar 

  7. 7.

    Gluzman IY, Francis SE, Oksman A, Smith CE, Duffin KL, Goldberg DE. Order and specificity of the Plasmodium falciparum hemoglobin degradation pathway. J Clin Invest. 1994;93:1602–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Müller J, Hemphill A. Drug target identification in protozoan parasites. Expert Opin Drug Discovery. 2016;11(8):815–24.

    Google Scholar 

  9. 9.

    Funk VA, Olafson RW, Raap M, Smith D, Aitken L, Haddow JD, Wang D, Dawson-Coates JA, Burke RD, Miller KM. Identification, characterization and deduced amino acid sequence of the dominant protease from Kudoa paniformis and K. thyrsites: A unique cytoplasmic cysteine protease. Comp Biochem Physiol B Biochem Mol Biol. 2008;149(3):477–89.

    PubMed  Google Scholar 

  10. 10.

    Shin SP, Zenke K, Yokoyama H. Characterization of proteases isolated from Kudoa septempunctata. J Vet Res. 2015;55(3):175–9.

    Google Scholar 

  11. 11.

    Chang ES, Neuhof M, Rubinstein ND, Diamant A, Philippe H, Huchon D, Cartwright P. Genomic insights into the evolutionary origin of Myxozoa within Cnidaria. Proc Natl Acad Sci U S A. 2015;112(48):14912–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Yang Y, Xiong J, Zhou Z, Huo F, Miao W, Ran C, Liu Y, Zhang J, Feng J, Wang M, Wang M, Wang L, Yao B. The genome of the myxosporean Thelohanellus kitauei shows adaptations to nutrient acquisition within its fish host. Genome Biol Evo. 2014;12:3182–98.

    Google Scholar 

  13. 13.

    Feist SW, Morris DJ, Alama-Bermejo G, Holzer AS. In: Okamura B, Gruhl A, Bartholomew JL, editors. Cellular process in myxozoans in Myxozoa Evolution. Ecology and Development. Cham: Springer International Publishing; 2015.

    Google Scholar 

  14. 14.

    Schmidt-Posthaus H, Wahli T. Host and environmental influences on development of disease. In: Okamura B, Gruhl A, Bartholomew J, editors. Myxozoan evolution, ecology and development. Cham: Springer; 2015.

    Google Scholar 

  15. 15.

    Lom J, Dyková I. Protozoan parasites of fishes. Amsterdam: Elsevier Science Publishers; 1992.

    Google Scholar 

  16. 16.

    Okamura B, Gruhl A, Reft AJ. Cnidarian origins of the Myxozoa. In: Myxozoan evolution, ecology and development. Cham: Springer; 2015. p. 45–68.

    Google Scholar 

  17. 17.

    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 Ecol. 2018;27:1651–66.

    PubMed  Google Scholar 

  18. 18.

    Wolf K, Markiw ME. Biology contravenes taxonomy in the Myxozoa: New discoveries show alternation of invertebrate and vertebrate hosts. Science. 1984;225(4669):1449–52.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Eszterbauer E, Atkinson S, Diamant A, Morris D, El-Matbouli M, Hartikainen H. Myxozoan life cycles: practical approaches and insights. In: Myxozoan evolution, ecology and development. Cham: Springer; 2015. p. 175–98.

    Google Scholar 

  20. 20.

    Arthur JR, Lom J. Sphaerospora araii n. sp. (Myxosporea: Sphaerosporidae) from the kidney of a longnose skate (Raja rhina Jordan and Gilbert) from the Pacific Ocean off Canada. Can J Zool. 1985;63:2902–6

    Google Scholar 

  21. 21.

    Bartošová P, Fiala I, Jirků M, Cinková M, Caffara M, Fioravanti ML, Atkinson SD, Bartholomew JL, Holzer AS. Sphaerospora sensu stricto: Taxonomy, diversity and evolution of a unique lineage of myxosporeans (Myxozoa). Mol Phylogenet Evol. 2013;68(1):93–105.

    PubMed  Google Scholar 

  22. 22.

    Desser SS, Lom J, Dyková I. Developmental stages of Sphaerospora ohlmacheri (Whinery, 1893) n.comb. (Myxozoa: Myxosporea) in the renal tubules of bullfrog tadpoles, Rana catesbeiana, from Lake of Two Rivers, Algonquin Park, Ontario. Can J Zool. 2011;64:2213–7.

    Google Scholar 

  23. 23.

    Jirků M, Fiala I, Modrý D. Tracing the genus Sphaerospora: rediscovery, redescription and phylogeny of the Sphaerospora ranae (Morelle,) n. comb. (Myxosporea, Sphaerosporidae), with emendation of the genus Sphaerospora. Parasitol. 2007;134(12):1727–39.

    Google Scholar 

  24. 24.

    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(1):347.

    PubMed  PubMed Central  Google Scholar 

  25. 25.

    Baska F, Molnár K. Blood stages of Sphaerospora spp. (Myxosporea) in cyprinid fishes. Dis Aq Org. 1988;5:23–8.

    Google Scholar 

  26. 26.

    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. Parasitol. 1983;86:529–35.

    Google Scholar 

  27. 27.

    Lom J, Pavlásková M, Dyková I. Notes on kidney infecting species of the genus Sphaerospora Thelohan (Myxosporea), including a new species Sphaerospora gobionis sp. nov., and on myxosporean life cycle stages in stages in the blood of some freshwater fish. J Fish Dis. 1985;8:221–32.

    Google Scholar 

  28. 28.

    Hartigan A, Estensoro I, Vancová M, et al. New cell motility model observed in parasitic cnidarian Sphaerospora molnari (Myxozoa:Myxosporea) blood stages in fish. Sci Rep. 2016;6:39093.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Korytář T, Wiegertjes G, Zusková E, Tomanová A, Lisnerová M, Patra S, Sieranski V, Šíma R, Born-Torrijos A, Wentzel A, Blasco-Monleon S, Yanes-Roca C, Policar T, Holzer AS. The kinetics of cellular and humoral immune responses of common carp to presporogonic development of the myxozoan Sphaerospora molnari. Parasit Vectors. 2019;2:208.

    Google Scholar 

  30. 30.

    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.

    PubMed  PubMed Central  Google Scholar 

  31. 31.

    Simão F, Waterhouse R, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31(19):3210–2.

    PubMed  Google Scholar 

  32. 32.

    Verma S, Dixit R, Pandey KC. Cysteine proteases: modes of activation and future prospects as pharmacological targets. Front Pharmacol. 2016;7:107.

    PubMed  PubMed Central  Google Scholar 

  33. 33.

    Stack CM, Caffrey CR, Donnelly SM, Seshaadri A, Lowther J, Tort JF, Collins PR, Robinson MW, Xu W, McKerrow JH, Craik CS, Geiger SR, Marion R, Brinen LS, Dalton JP. Structural and functional relationships in the virulence-associated cathepsin L proteases of the parasitic liver fluke, Fasciola hepatica. J Biol Chem. 2007;283(15):9896–908.

    PubMed  PubMed Central  Google Scholar 

  34. 34.

    Weihofen A, Binn K, Lemberg MK, Ashman K, Martoglio B. Identification of signal peptide peptidase, a presenilin-type aspartic protease. Science. 2002;296:2215–8.

    CAS  PubMed  Google Scholar 

  35. 35.

    Gul IS, Staal J, Hulpiau P, De Keuckelaere E, Kamm K, Deroo T, Sanders E, Staes K, Driege Y, Saeys Y, Beyaert R, Technau U, Schierwater B, van Roy F. GC content of early metazoan genes and its impact on gene expression levels in mammalian cell lines. Genome Biol Evol. 2018;10(3):909–17.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Yanow SK, Purcell LA, Lee M, Spithill TW. Genomics-based drug design targets the AT-rich malaria parasite: implications for antiparasite chemotherapy. Pharmacogenomics. 2007;8(9):1267–72.

    CAS  PubMed  Google Scholar 

  37. 37.

    Eitel M, Francis WR, Varoqueaux F, Daraspe J, Osigus HJ, Krebs S, Vargas S, Blum H, Williams GA, Schierwater B, Wörheide G. Comparative genomics and the nature of placozoan species. PLoS Biol. 2018;16(7):e2005359.

    PubMed  PubMed Central  Google Scholar 

  38. 38.

    Kenny NJ, de Goeij JM, de Bakker DM, Whalen CG, Berezikov E, Riesgo A. Towards the identification of ancestrally shared regenerative mechanisms across the Metazoa: a Transcriptomic case study in the Demosponge Halisarca caerulea. Mar Genomics. 2018;37:135–47.

    PubMed  Google Scholar 

  39. 39.

    Celis JS, Wibberg D, Ramírez-Portilla C, Rupp O, Sczyrba A, Winkler A, Kalinowski J, Wilke T. Binning enables efficient host genome reconstruction in cnidarian holobionts. GigaScience. 2018;7:7.

    Google Scholar 

  40. 40.

    Pootakham W, Naktang C, Sonthirod C, Yoocha T, Sangsrakru D, Jomchai N, Putchim L, Tangphatsornruang S. Development of a novel reference transcriptome for scleractinian coral Porites lutea using single-molecule long-read isoform sequencing (Iso-Seq). Front Mar Sci. 2018;5

  41. 41.

    Bankers L, Neiman M, et al. G3 (Bethseda). 2017;7(3):871–80.

    CAS  Google Scholar 

  42. 42.

    McNeilly TN, Few D, Burgess STG, Wright H, Bartley DJ, Bartley Y, Nisbet AJ. Niche-specific gene expression in a parasitic nematode; increased expression of immunomodulators in Teladorsagia circumcincta larvae derived from host mucosa. Sci Rep. 2017;7:7214.

    PubMed  PubMed Central  Google Scholar 

  43. 43.

    Li W, Liu B, Yang Y, Ren Y, Wang S, Liu C, Zhang N, Qu Z, Yang W, Zhang Y, Yan H, Jiang F, Li L, Li S, Jia W, Yin H, Cai X, Liu T, DP MM, Fan W, Fu B. The genome of tapeworm Taenia multiceps sheds light on understanding parasitic mechanism and control of coenurosis disease. DNA Res. 2018;25(5):499–510.

    PubMed  PubMed Central  Google Scholar 

  44. 44.

    Gayral P, Melo-Ferreira J, Glémin S, Bienre N, Carneiro M, Nabholz B, Lourenco JM, Alves PC, Ballenghein M, Faivre N. Reference-free population genomics from next-generation transcriptome data and the vertebrate–invertebrate gap. PLoS Genet. 2013;9:e1003457.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Martin JA, Wang Z. Next-generation transcriptome assembly. Nat Rev Genet. 2011;12:671.

    CAS  PubMed  Google Scholar 

  46. 46.

    Schrul B, Kapp K, Sinning I, Dobberstein B. Signal peptide peptidase (SPP) assembles with substrates and misfolded membrane proteins into distinct oligomeric complexes. Biochem J. 2010;427(3):523–34.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Voss M, Schröder B, Fluhrer R. Mechanism, specificity, and physiology of signal peptide peptidase (SPP) and SPP-like proteases. Biochim Biophys Acta. 2013;1828(12):2828–39.

    CAS  PubMed  Google Scholar 

  48. 48.

    Li X, Chen H, Bahamontes-Rosa N, Kun JFJ, Traore B, Crompton PD, Chishti AH. Plasmodium falciparum signal peptide peptidase is a promising drug target against blood stage malaria. Biochem Biophys Res Commun. 2009;380(3):454–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  49. 49.

    Lázaro S, Gamarra D, Del Val M. Proteolytic enzymes involved in MHC class I antigen processing: A guerrilla army that partners with the proteasome. Mol Immunol. 2015;68(2):72–6.

    PubMed  Google Scholar 

  50. 50.

    Parvanova I, Epiphanio S, Fauq A, Golde TE, Prudêncio M, Mota MM. A small molecule inhibitor of signal peptide peptidase inhibits Plasmodium development in the liver and decreases malaria severity. PLoS One. 2009;4(4):e5078.

    PubMed  PubMed Central  Google Scholar 

  51. 51.

    Harbut MB, Patel BA, Yeung BKS, McNamara CW, Bright AT, Ballard J, Greenbaum DC. Targeting the ERAD pathway via inhibition of signal peptide peptidase for antiparasitic therapeutic design. Proc Natl Acad Sci U S A. 2012;109(52):21486–91.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Eggleson KK, Duffin KL, Goldberg DE. Identification and characterization of falcilysin, a metallopeptidase involved in hemoglobin catabolism within the malaria parasite Plasmodium falciparum. J Biol Chem. 1999;274(45):32411–7.

    CAS  PubMed  Google Scholar 

  53. 53.

    Bruschi F, Pinto B. The significance of matrix metalloproteinases in parasitic infections involving the central nervous system. Pathogens. 2013;2(1):105–29.

    PubMed  PubMed Central  Google Scholar 

  54. 54.

    McGwire BS, Chang KP, Engman DM. Migration through the extracellular matrix by the parasitic protozoan Leishmania is enhanced by surface metalloprotease gp63. Infect Immun. 2003;71(2):1008–10.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Geurts N, Opdenakker G, Van den Steen PE. Matrix metalloproteinases as therapeutic targets in protozoan parasitic infections. Pharmacol Ther. 2012;133(3):257–79.

    CAS  PubMed  Google Scholar 

  56. 56.

    Piña-Vázquez C, Reyes-López M, Ortíz-Estrada G, de la Garza M, Serrano-Luna J. Host-parasite interaction: parasite-derived and -induced proteases that degrade human extracellular matrix. J Parasitol Res. 2012;2012:748206.

    PubMed  PubMed Central  Google Scholar 

  57. 57.

    Ramos OH, Selistre-de-Araujo HS. Snake venom metalloproteases--structure and function of catalytic and disintegrin domains. Comp Biochem Physiol C Toxicol Pharmacol. 2006;142(3–4):328–46.

    CAS  PubMed  Google Scholar 

  58. 58.

    Decrem Y, Beaufays J, Blasioli V, Lahaye K, Brossard M, Vanhamme L, Godfroid E. A family of putative metalloproteases in the salivary glands of the tick Ixodes ricinus. FEBS J. 2008;275(7):1485–99.

    CAS  PubMed  Google Scholar 

  59. 59.

    Wagner L, Klemann C, Stephan M, von Hörsten S. Unravelling the immunological roles of dipeptidyl peptidase 4 (DPP4) activity and/or structure homologue (DASH) proteins. Clin Exp Immunol. 2016;184:265–83.

    CAS  PubMed  PubMed Central  Google Scholar 

  60. 60.

    Touz MC, Nores MJ, Slavin I, Piacenza L, Acosta D, Carmona C, Lujan HD. Membrane-associated dipeptidyl peptidase IV is involved in encystation-specific gene expression during Giardia differentiation. Biochem J. 2002;364(3):703–10.

    CAS  PubMed  PubMed Central  Google Scholar 

  61. 61.

    Geldhof P, Knox D. The intestinal contortin structure in Haemonchus contortus: an immobilised anticoagulant? Int J Parasitol. 2008;38:1579–88.

    CAS  PubMed  Google Scholar 

  62. 62.

    Robinson MW, Dalton JP. Cysteine proteases of pathogenic organisms. Landes Bioscience. 2011.

  63. 63.

    Konagaya S. Studies on the jellied meat of fish, with special reference to that of yellowfin tuna. Bull Tokai Reg Fish Res Lab. 1984;114:1–101.

    Google Scholar 

  64. 64.

    Martone CB, Spivak E, Busconi L, Folco EJE, Sánchez JJ. A cysteine protease from myxosporean degrades host myofibrils in vitro. Comp Biochem Physiol B Biochem Mol Biol. 1999;123:267–72.

    Google Scholar 

  65. 65.

    Caffrey CR, Goupil L, Rebello KM, Dalton JP, Smith D. Cysteine proteases as digestive enzymes in parasitic helminths. PLoS Neg Trop Dis. 2018;12(8):e0005840.

    Google Scholar 

  66. 66.

    Dalton J, Caffrey CR, Sajid M, Stack C, Donnelly S, Loukas A, Don T, Mckerrow J, Halton DW, Brindley P. Proteases in trematode biology. In: Parasitic Flatworms: Molecular biology, Biochemistry, Immunology and Physiology; 2006. p. 348–68.

    Google Scholar 

  67. 67.

    McKerrow JH. The diverse roles of cysteine proteases in parasites and their suitability as drug targets. PLoS Negl Trop Dis. 2018;12(8):e0005639.

    PubMed  PubMed Central  Google Scholar 

  68. 68.

    Selzer PM, Pingel S, Hseih I, Ugele B, Chan VJ, Engel JC, Bogyo M, Russell DG, Sakanari JA, McKerrow JH. Cysteine protease inhibitors as chemotherapy: lessons from a parasite target. Proc Natl Acad Sci U S A. 1999;96(20):11015–22.

    CAS  PubMed  PubMed Central  Google Scholar 

  69. 69.

    Siqueira-Neto JL, Debnath A, McCall LI, Bernatchez JA, Ndao M, et al. Cysteine proteases in protozoan parasites. PLoS Negl Trop Dis. 2018;12(8):e0006512.

    PubMed  PubMed Central  Google Scholar 

  70. 70.

    Ferraro F, Merlino A, dell’ Oca N, Gil J, Tort JF, Gonzalez M, Cerecetto H, Cabrera M, Corvo I. Identification of chalcones as Fasciola hepatica cathepsin L inhibitors using a comprehensive experimental and computational approach. PLoS Negl Trop Dis. 2016;10(7):e0004834.

    PubMed  PubMed Central  Google Scholar 

  71. 71.

    Williamson AL, Lustigman S, Oksov Y, Deumic V, Plieskatt J, Mendez S, Zhan B, Bottazzi ME, Hotez PJ, Loukas A. Ancylostoma caninum MTP-1, an astacin-like metalloprotease secreted by infective hookworm larvae, is involved in tissue migration. Infect Immun. 2006;74(2):961–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  72. 72.

    International Helminth Genomes Consortium. Comparative genomics of the major parasitic worms. Nat Genet. 2018;51(1):163–74.

    PubMed Central  Google Scholar 

  73. 73.

    Hahn J, Seeber F, Kolodziej H, Ignatius R, Laue M, Aebischer T, Klotz C. High sensitivity of Giardia duodenalis to tetrahydrolipstatin (orlistat) in vitro. PLoS One. 2013;8(8):e71597.

    CAS  PubMed  PubMed Central  Google Scholar 

  74. 74.

    Shakarian AM, McGugan GC, Joshi MB, Bowers L, Ganim C, Barowski J, Dwyer DM. Identification, characterization, and expression of a unique secretory lipase from the human pathogen Leishmania donovani. Mol Cell Biochem. 2010;341(1–2):17–31.

    CAS  PubMed  PubMed Central  Google Scholar 

  75. 75.

    Chaves SP, Gomes DCO, De-Simone SG, Rossi-Bergmann B, de Matos HL. Serine proteases and vaccines against Leishmaniasis: a dual role. J Vaccine Vaccin. 2015;6:1.

    Google Scholar 

  76. 76.

    Stutzer C, Richards SA, Ferreira M, Baron S, Maritz-Olivier C. Metazoan parasite vaccines: present status and future prospects. Front Cell Infect Microbiol. 2018;8:67.

    PubMed  PubMed Central  Google Scholar 

  77. 77.

    de Vries E, Bakker N, Krijgsveld J, Knox DP, Heck AJ, Yatsuda AP. An AC-5 cathepsin B-like protease purified from Haemonchus contortus excretory secretory products shows protective antigen potential for lambs. Vet Res. 2009;40(4):1–11.

    Google Scholar 

  78. 78.

    Sommerset I, Krossøy B, Biering E, Frost P. Vaccines for fish aquaculture. Expert Rev Vaccines. 2005;4:89–101.

    CAS  Article  PubMed  Google Scholar 

  79. 79.

    Carraro L, Bertuzzo E, Mari L, Fontes I, Hartikainen H, Strepparava N, Schmidt-Posthaus H, Wahli T, Jokela J, Gatto M, Rinaldo A. Integrated field, laboratory, and theoretical study of PKD spread in a Swiss prealpine river. Proc Natl Acad Sci U S A. 2017;114(45):11992–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  80. 80.

    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(4):735–53.

    Google Scholar 

  81. 81.

    Bruneaux M, Visse M, Gross R, Pukk L, Saks L, Vasemägi A. Parasite infection and decreased thermal tolerance: impact of proliferative kidney disease on a wild salmonid fish in the context of climate change. Funct Ecol. 2017;31(1):216–26.

    Google Scholar 

  82. 82.

    Debes PV, Gross R, Vasemägi A. Quantitative genetic variation in, and environmental effects on, pathogen resistance and temperature-dependent disease severity in a wild trout. Am Nat. 2017;190(2):244–65.

    PubMed  Google Scholar 

  83. 83.

    Skovgaard A, Buchmann K. Tetracapsuloides bryosalmonae and PKD in juvenile wild salmonids in Denmark. Dis Aq Org. 2012;101(1):33–42.

    CAS  Google Scholar 

  84. 84.

    Robbins J. Tiny invader, deadly to fish, shuts down a river in Montana. In: The New York Times; 2016. Accessed: 10 Oct 2016.

    Google Scholar 

  85. 85.

    Hutchins PR, Sepulveda AJ, Martin RM, Hopper LR. A probe-based quantitative PCR assay for detecting Tetracapsuloides bryosalmonae in fish tissue and environmental DNA water samples. Cons Gen Resource. 2018;10(3):317–9.

    Google Scholar 

  86. 86.

    Min XJ, Butler G, Storms R, Tsang A. OrfPredictor: predicting protein-coding regions in EST-derived sequences. Nucleic Acids Res. 2005;33:W677–80.

    CAS  PubMed  PubMed Central  Google Scholar 

  87. 87.

    Li W. Godzik Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22(13):1658–9.

    CAS  PubMed  Google Scholar 

  88. 88.

    Rawlings ND, Waller M, Barrett AJ. Bateman MEROPS: the database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 2014;42:D503–9.

    CAS  PubMed  Google Scholar 

  89. 89.

    Stefanik DJ, Lubinski TJ, Granger BR, Byrd AL, Reitzel AM, DeFilippo L, Lorenc A, Finnerty JR. Production of a reference transcriptome and transcriptomic database (EdwardsiellaBase) for the lined sea anemone, Edwardsiella lineata, a parasitic cnidarian. BMC Genomics. 2014;15:71.

    PubMed  PubMed Central  Google Scholar 

  90. 90.

    Kelley LA, Mezulis S, Yates CM, Wass MN, Sternberg MJE. The Phyre2 web portal for protein modeling, prediction and analysis. Nat Protoc. 2015;10:845.

    CAS  PubMed  PubMed Central  Google Scholar 

  91. 91.

    Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:e45.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references


The authors acknowledge BGI Genomics for their sequencing services and are grateful for the support of our fish farmers in Czech Republic and Hungary.


The study was funded by the European Commission Horizon 2020 Research and Innovation Action (project no. 634429, ParaFishControl) for travel, reagents, sequencing and salary support, the Czech Science Foundation (project no. 19-28399X; AQUAPARA-OMICS) for travel and lab consumables and the Hungarian National Research, Development and Innovation Office (project no. NN124220) funded travel and lab consumables. No funding bodies played a role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information




ASH conceived and planned the project; AH, EE and ASH collected RNA and DNA samples; AH performed transcriptome analyses and structural modelling, AK and HP performed control sequencing of proteases and rDNA, AK designed and performed qPCR assays; AH and ASH wrote the manuscript. All authors agreed to the final version of it.

Corresponding author

Correspondence to Ashlie Hartigan.

Ethics declarations

Ethics approval and consent to participate

Fish were obtained commercially, with the owners of the carp ponds providing us with the fish directly on site. Fish were transported live to the laboratory and thereafter, fish were anaesthetized with clove oil and blood was taken. All animal procedures were performed in accordance with Czech legislation (section 29 of Protection of Animals Against Cruelty Act No. 246/1992) and approved by the Czech Ministry of Agriculture. We declare that animal handling complied with the relevant European and international guidelines on animal welfare, namely Directive 2010/63/EU on the protection of animals used for scientific purposes and the guidelines and recommendations of the Federation of Laboratory Animal Science Associations.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file1.

Supplementary Data 1: Table of primers used in qPCR assay including house keeping genes.

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Hartigan, A., Kosakyan, A., Pecková, H. et al. Transcriptome of Sphaerospora molnari (Cnidaria, Myxosporea) blood stages provides proteolytic arsenal as potential therapeutic targets against sphaerosporosis in common carp. BMC Genomics 21, 404 (2020).

Download citation


  • Myxozoa
  • In silico screening
  • Proteases
  • Aquaculture
  • Parasite
  • Drug targets