Comparative transcriptomics and host-specific parasite gene expression profiles inform on drivers of proliferative kidney disease

The Tilapia Cyst Tissue Enclosing the Proliferating Myxobolus bejeranoi Parasite Exhibits Cornified Structure and Immune Barrier Function

Myxozoa, a unique group of obligate endoparasites within the phylum Cnidaria, can cause emerging diseases in wild and cultured fish populations. Recently, the myxozoan Myxobolus bejeranoi has been identified as a prevalent pathogen infecting the gills of cultured hybrid tilapia, leading to systemic immune suppression and considerable mortality. Here, we employed a proteomic approach to examine the impact of M. bejeranoi infection on fish gills, focusing on the structure of the granulomata, or cyst, formed around the proliferating parasite to prevent its spread to surrounding tissue. Enrichment analysis showed increased immune response and oxidative stress in infected gill tissue, most markedly in the cyst's wall. The intense immune reaction included a consortium of endopeptidase inhibitors, potentially combating the myxozoan arsenal of secreted proteases. Analysis of the cyst's proteome and histology staining indicated that keratin intermediate filaments contribute to its structural rigidity. Moreover, we uncovered skin-specific proteins, including a grainyhead-like transcription factor and a teleost-specific S100 calcium-binding protein that may play a role in epithelial morphogenesis and cysts formation. These findings deepen our understanding of the proteomic elements that grant the cyst its distinctive nature at the critical interface between the fish host and myxozoan parasite.

The myxozoans Myxobolus cerebralis and Tetracapsuloides bryosalmonae modulate rainbow trout immune responses: quantitative shotgun proteomics at the portals of entry after single and co-infections

Introduction

Little is known about the proteomic changes at the portals of entry in rainbow trout after infection with the myxozoan parasites, Myxobolus cerebralis, and Tetracapsuloides bryosalmonae. Whirling disease (WD) is a severe disease of salmonids, caused by the myxosporean M. cerebralis, while, proliferative kidney disease (PKD) is caused by T. bryosalmonae, which instead belongs to the class Malacosporea. Climate change is providing more suitable conditions for myxozoan parasites lifecycle, posing a high risk to salmonid aquaculture and contributing to the decline of wild trout populations in North America and Europe. Therefore, the aim of this study was to provide the first proteomic profiles of the host in the search for evasion strategies during single and coinfection with M. cerebralis and T. bryosalmonae.

Methods

One group of fish was initially infected with M. cerebralis and another group with T. bryosalmonae. After 30 days, half of the fish in each group were co-infected with the other parasite. Using a quantitative proteomic approach, we investigated proteomic changes in the caudal fins and gills of rainbow trout before and after co-infection.

Results

In the caudal fins, 16 proteins were differentially regulated post exposure to M. cerebralis, whereas 27 proteins were differentially modulated in the gills of the infected rainbow trout post exposure to T. bryosalmonae. After co-infection, 4 proteins involved in parasite recognition and the regulation of host immune responses were differentially modulated between the groups in the caudal fin. In the gills, 11 proteins involved in parasite recognition and host immunity, including 4 myxozoan proteins predicted to be virulence factors, were differentially modulated.

Discussion

The results of this study increase our knowledge on rainbow trout co-infections by myxozoan parasites and rainbow trout immune responses against myxozoans at the portals of entry, supporting a better understanding of these host-parasite interactions.

The Molecular Mechanisms Employed by the Parasite Myxobolus bejeranoi (Cnidaria: Myxozoa) from Invasion through Sporulation for Successful Proliferation in Its Fish Host

Myxozoa is a unique group of obligate endoparasites in the phylum Cnidaria that can cause emerging diseases in wild and cultured fish populations. Recently, we identified a new myxozoan species, Myxobolus bejeranoi, which infects the gills of cultured tilapia while suppressing host immunity. To uncover the molecular mechanisms underlying this successful parasitic strategy, we conducted transcriptomics analysis of M. bejeranoi throughout the infection. Our results show that histones, which are essential for accelerated cell division, are highly expressed even one day after invasion. As the infection progressed, conserved parasitic genes that are known to modulate the host immune reaction in different parasitic taxa were upregulated. These genes included energy-related glycolytic enzymes, as well as calreticulin, proteases, and miRNA biogenesis proteins. Interestingly, myxozoan calreticulin formed a distinct phylogenetic clade apart from other cnidarians, suggesting a possible function in parasite pathogenesis. Sporogenesis was in its final stages 20 days post-exposure, as spore-specific markers were highly expressed. Lastly, we provide the first catalog of transcription factors in a Myxozoa species, which is minimized compared to free-living cnidarians and is dominated by homeodomain types. Overall, these molecular insights into myxozoan infection support the concept that parasitic strategies are a result of convergent evolution.

Differentially expressed transcripts of Tetracapsuloides bryosalmonae (Cnidaria) between carrier and dead-end hosts involved in key biological processes: novel insights from a coupled approach of FACS and RNA sequencing

Tetracapsuloides bryosalmonae is a malacosporean endoparasite that infects a wide range of salmonids and causes proliferative kidney disease (PKD). Brown trout serves as a carrier host whereas rainbow trout represents a dead-end host. We thus asked if the parasite adapts to the different hosts by changing molecular mechanisms. We used fluorescent activated cell sorting (FACS) to isolate parasites from the kidney of brown trout and rainbow trout following experimental infection with T. bryosalmonae. The sorted parasite cells were then subjected to RNA sequencing. By this approach, we identified 1120 parasite transcripts that were expressed differentially in parasites derived from brown trout and rainbow trout. We found elevated levels of transcripts related to cytoskeleton organisation, cell polarity, peptidyl-serine phosphorylation in parasites sorted from brown trout. In contrast, transcripts related to translation, ribonucleoprotein complex biogenesis and subunit organisation, non-membrane bounded organelle assembly, regulation of protein catabolic process and protein refolding were upregulated in rainbow trout-derived parasites. These findings show distinct molecular adaptations of parasites, which may underlie their distinct outcomes in the two hosts. Moreover, the identification of these differentially expressed transcripts may enable the identification of novel drug targets that may be exploited as treatment against T. bryosalmonae. We here also describe for the first time how FACS based isolation of T. bryosalmonae cells from infected kidney of fish fosters research and allows to define differentially expressed parasite transcripts in carrier and dead-end fish hosts.

Apoptotic gene loss in Cnidaria is associated with transition to parasitism

The phylum Cnidaria consists of several morphologically diverse classes including Anthozoa, Cubozoa, Hydrozoa, Polypodiozoa, Scyphozoa, Staurozoa, and Myxozoa. Myxozoa comprises two subclasses of obligate parasites-Myxosporea and Malacosporea, which demonstrate various degrees of simplification. Myxosporea were previously reported to lack the majority of core protein domains of apoptotic proteins including caspases, Bcl-2, and APAF-1 homologs. Other sequenced Cnidaria, including the parasite Polypodium hydriforme from Polypodiozoa do not share this genetic feature. Whether this loss of core apoptotic proteins is unique to Myxosporea or also present in its sister subclass Malacosporea was not previously investigated. We show that the presence of core apoptotic proteins gradually diminishes from free-living Cnidaria to Polypodium to Malacosporea to Myxosporea. This observation does not favor the hypothesis of catastrophic simplification of Myxosporea at the genetic level, but rather supports a stepwise adaptation to parasitism that likely started from early parasitic ancestors that gave rise to Myxozoa.

Red blood cells in proliferative kidney disease-rainbow trout (Oncorhynchus mykiss) infected by Tetracapsuloides bryosalmonae harbor IgM+ red blood cells

The myxozoan parasite Tetracapsuloides bryosalmonae is the causative agent of proliferative kidney disease (PKD)-a disease of salmonid fishes, notably of the commercially farmed rainbow trout Oncorhynchus mykiss. Both wild and farmed salmonids are threatened by this virulent/deadly disease, a chronic immunopathology characterized by massive lymphocyte proliferation and hyperplasia, which manifests as swollen kidneys in susceptible hosts. Studying the immune response towards the parasite helps us understand the causes and consequences of PKD. While examining the B cell population during a seasonal outbreak of PKD, we unexpectedly detected the B cell marker immunoglobulin M (IgM) on red blood cells (RBCs) of infected farmed rainbow trout. Here, we studied the nature of this IgM and this IgM⁺ cell population. We verified the presence of surface IgM via parallel approaches: flow cytometry, microscopy, and mass spectrometry. The levels of surface IgM (allowing complete resolution of IgM^- RBCs from IgM⁺ RBCs) and frequency of IgM⁺ RBCs (with up to 99% of RBCs being positive) have not been described before in healthy fishes nor those suffering from disease. To assess the influence of the disease on these cells, we profiled the transcriptomes of teleost RBCs in health and disease. Compared to RBCs originating from healthy fish, PKD fundamentally altered RBCs in their metabolism, adhesion, and innate immune response to inflammation. In summary, RBCs play a larger role in host immunity than previously appreciated. Specifically, our findings indicate that the nucleated RBCs of rainbow trout interact with host IgM and contribute to the immune response in PKD.

Identification of in vivo induced antigens of the malacosporean parasite Tetracapsuloides bryosalmonae (Cnidaria) using in vivo induced antigen technology

Tetracapsuloides bryosalmonae is a malacosporean endoparasite that causes proliferative kidney disease (PKD) in wild and farmed salmonids in Europe and North America. The life cycle of T. bryosalmonae completes between invertebrate bryozoan and vertebrate fish hosts. Inside the fish, virulence factors of T. bryosalmonae are induced during infection or interactions with host cells. T. bryosalmonae genes expressed in vivo are likely to be important in fish pathogenesis. Herein, we identify in vivo induced antigens of T. bryosalmonae during infection in brown trout (Salmo trutta) using in vivo induced antigen technology (IVIAT). Brown trout were exposed to the spores of T. bryosalmonae and were sampled at different time points. The pooled sera were first pre-adsorbed with antigens to remove false positive results. Subsequently, adsorbed sera were used to screen a T. bryosalmonae cDNA phage expression library. Immunoscreening analysis revealed 136 immunogenic T. bryosalmonae proteins induced in brown trout during parasite development. They are involved in signal transduction, transport, metabolism, ion-protein binding, protein folding, and also include hypothetical proteins, of so far unknown functions. The identified in vivo induced antigens will be useful in the understanding of T. bryosalmonae pathogenesis during infection in susceptible hosts. Some of the antigens found may have significant implications for the discovery of candidate molecules for the development of potential therapies and preventive measures against T. bryosalmonae in salmonids.

Genome-wide alternative splicing profile in the posterior kidney of brown trout (Salmo trutta) during proliferative kidney disease

BACKGROUND

The cnidarian myxozoan parasite Tetracapsuloides bryosalmonae causes chronic proliferative kidney disease (PKD) in salmonids. This parasite is a serious threat to wild and cultured salmonids. T. bryosalmonae undergoes intra-luminal sporogonic development in the kidney of brown trout (Salmo trutta) and the viable spores are released via urine. We investigated the alternative splicing pattern in the posterior kidney of brown trout during PKD.

RESULTS

RNA-seq data were generated from the posterior kidney of brown trout collected at 12 weeks post-exposure to T. bryosalmonae. Subsequently, this data was mapped to the brown trout genome. About 153 significant differently expressed alternatively spliced (DEAS) genes, (delta PSI = 5%, FDR P-value < 0.05) were identified from 19,722 alternatively spliced events. Among the DEAS genes, the least and most abundant alternative splicing types were alternative 5' splice site (5.23%) and exon skipping (70.59%), respectively. The DEAS genes were significantly enriched for sodium-potassium transporter activity and ion homeostasis (ahcyl1, atp1a3a, atp1a1a.1, and atp1a1a.5). The protein-protein interaction network analysis enriched two local network clusters namely cation transporting ATPase C-terminus and Sodium/potassium ATPase beta chain cluster, and mixed inclusion of Ion homeostasis and EF-hand domain cluster. Furthermore, the human disease-related salmonella infection pathway was significantly enriched in the protein-protein interaction network.

CONCLUSION

This study provides the first baseline information about alternative splicing in brown trout during PKD. The generated data lay a foundation for further functional molecular studies in PKD - brown trout infection model. The information generated from the present study can help to develop therapeutic strategies for PKD in the future.

Mitigating human impacts including climate change on proliferative kidney disease in salmonids of running waters

Over the last two decades, an increasing number of reports have identified a decline in salmonid populations, possibly linked to infection with the parasite Tetracapsuloides bryosalmonae and the corresponding disease, that is, proliferative kidney disease (PKD). The life cycle of this myxozoan parasite includes sessile bryozoan species as invertebrate host, which facilitates the distribution of the parasite in running waters. As the disease outcome is temperature dependent, the impact of the disease on salmonid populations is increasing with global warming due to climate change. The goal of this review is to provide a detailed overview of measures to mitigate the effects of PKD on salmonid populations. It first summarizes the parasite life cycle, temperature-driven disease dynamics and new immunological and molecular research into disease resistance and, based on this, discusses management possibilities. Sophisticated management actions focusing on local adaptation of salmonid populations, restoration of the riverine ecosystem and keeping water temperatures cool are necessary to reduce the negative effects of PKD. Such actions include temporary stocking with PKD-resistant salmonids, as this may assist in conserving current populations that fail to reproduce.

A myxozoan genome reveals mosaic evolution in a parasitic cnidarian

BACKGROUND

Parasite evolution has been conceptualized as a process of genetic loss and simplification. Contrary to this model, there is evidence of expansion and conservation of gene families related to essential functions of parasitism in some parasite genomes, reminiscent of widespread mosaic evolution-where subregions of a genome have different rates of evolutionary change. We found evidence of mosaic genome evolution in the cnidarian Myxobolus honghuensis, a myxozoan parasite of fish, with extremely simple morphology.

RESULTS

We compared M. honghuensis with other myxozoans and free-living cnidarians, and determined that it has a relatively larger myxozoan genome (206 Mb), which is less reduced and less compact due to gene retention, large introns, transposon insertion, but not polyploidy. Relative to other metazoans, the M. honghuensis genome is depleted of neural genes and has only the simplest animal immune components. Conversely, it has relatively more genes involved in stress resistance, tissue invasion, energy metabolism, and cellular processes compared to other myxozoans and free-living cnidarians. We postulate that the expansion of these gene families is the result of evolutionary adaptations to endoparasitism. M. honghuensis retains genes found in free-living Cnidaria, including a reduced nervous system, myogenic components, ANTP class Homeobox genes, and components of the Wnt and Hedgehog pathways.

CONCLUSIONS

Our analyses suggest that the M. honghuensis genome evolved as a mosaic of conservative, divergent, depleted, and enhanced genes and pathways. These findings illustrate that myxozoans are not as genetically simple as previously regarded, and the evolution of some myxozoans is driven by both genomic streamlining and expansion.

Proteases as Therapeutic Targets Against the Parasitic Cnidarian Ceratonova shasta: Characterization of Molecules Key to Parasite Virulence In Salmonid Hosts

Proteases and their inhibitors play critical roles in host-parasite interactions and in the outcomes of infections. Ceratonova shasta is a myxozoan pathogen that causes enteronecrosis in economically important salmonids from the Pacific Northwest of North America. This cnidarian parasite has host-specific genotypes with varying virulence, making it a powerful system to decipher virulence mechanisms in myxozoans. Using C. shasta genome and transcriptome, we identified four proteases of different catalytic types: cathepsin D (aspartic), cathepsin L and Z-like (cysteine) and aminopeptidase-N (metallo); and a stefin (cysteine protease inhibitor), which implied involvement in virulence and hence represent target molecules for the development of therapeutic strategies. We characterized, annotated and modelled their 3D protein structure using bioinformatics and computational tools. We quantified their expression in C. shasta genotype 0 (low virulence, no mortality) and IIR (high virulence and mortality) in rainbow trout Oncorhynchus mykiss, to demonstrate that there are major differences between the genotypes during infection and parasite development. High proliferation of genotype IIR was associated with high expression of the cathepsin D and the stefin, likely correlated with high nutrient demands and to regulate cell metabolism, with upregulation preceding massive proliferation and systemic dispersion. In contrast, upregulation of the cathepsin L and Z-like cysteine proteases may have roles in host immune evasion in genotype 0 infections, which are associated with low proliferation, low inflammation and non-destructive development. In contrast to the other proteases, C. shasta aminopeptidase-N appears to have a prominent role in nematocyst formation in both genotypes, but only during sporogenesis. Homology searches of C. shasta proteases against other myxozoan transcriptomes revealed a high abundance of cathepsin L and aminopeptidase homologs suggesting common gene requirements across species. Our study identified molecules of potential therapeutic significance for aquaculture and serves as a baseline for future research aimed at functional characterisation of these targets.

To React or Not to React: The Dilemma of Fish Immune Systems Facing Myxozoan Infections

Myxozoans are microscopic, metazoan, obligate parasites, belonging to the phylum Cnidaria. In contrast to the free-living lifestyle of most members of this taxon, myxozoans have complex life cycles alternating between vertebrate and invertebrate hosts. Vertebrate hosts are primarily fish, although they are also reported from amphibians, reptiles, trematodes, mollusks, birds and mammals. Invertebrate hosts include annelids and bryozoans. Most myxozoans are not overtly pathogenic to fish hosts, but some are responsible for severe economic losses in fisheries and aquaculture. In both scenarios, the interaction between the parasite and the host immune system is key to explain such different outcomes of this relationship. Innate immune responses contribute to the resistance of certain fish strains and species, and the absence or low levels of some innate and regulatory factors explain the high pathogenicity of some infections. In many cases, immune evasion explains the absence of a host response and allows the parasite to proliferate covertly during the first stages of the infection. In some infections, the lack of an appropriate regulatory response results in an excessive inflammatory response, causing immunopathological consequences that are worse than inflicted by the parasite itself. This review will update the available information about the immune responses against Myxozoa, with special focus on T and B lymphocyte and immunoglobulin responses, how these immune effectors are modulated by different biotic and abiotic factors, and on the mechanisms of immune evasion targeting specific immune effectors. The current and future design of control strategies for myxozoan diseases is based on understanding this myxozoan-fish interaction, and immune-based strategies such as improvement of innate and specific factors through diets and additives, host genetic selection, passive immunization and vaccination, are starting to be considered.

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

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

Advances and Discoveries in Myxozoan Genomics

Myxozoans are highly diverse and globally distributed cnidarian endoparasites in freshwater and marine habitats. They have adopted a heteroxenous life cycle, including invertebrate and fish hosts, and have been associated with diseases in aquaculture and wild fish stocks. Despite their importance, genomic resources of myxozoans have proven difficult to obtain due to their miniaturized and derived genome character and close associations with fish tissues. The first 'omic' datasets have now become the main resource for a better understanding of host-parasite interactions, virulence, and diversity, but also the evolutionary history of myxozoans. In this review, we discuss recent genomic advances in the field and outline outstanding questions to be answered with continuous and improved efforts of generating myxozoan genomic data.

Know your enemy - transcriptome of myxozoan Tetracapsuloides bryosalmonae reveals potential drug targets against proliferative kidney disease in salmonids

The myxozoan Tetracapsuloides bryosalmonae is a widely spread endoparasite that causes proliferative kidney disease (PKD) in salmonid fish. We developed an in silico pipeline to separate transcripts of T. bryosalmonae from the kidney tissue of its natural vertebrate host, brown trout (Salmo trutta). After stringent filtering, we constructed a partial transcriptome assembly T. bryosalmonae, comprising 3427 transcripts. Based on homology-restricted searches of the assembled parasite transcriptome and Atlantic salmon (Salmo salar) proteome, we identified four protein targets (Endoglycoceramidase, Legumain-like protease, Carbonic anhydrase 2, Pancreatic lipase-related protein 2) for the development of anti-parasitic drugs against T. bryosalmonae. Earlier work of these proteins on parasitic protists and helminths suggests that the identified anti-parasitic drug targets represent promising chemotherapeutic candidates also against T. bryosalmonae, and strengthen the view that the known inhibitors can be effective in evolutionarily distant organisms. In addition, we identified differentially expressed T. bryosalmonae genes between moderately and severely infected fish, indicating an increased abundance of T. bryosalmonae sporogonic stages in fish with low parasite load. In conclusion, this study paves the way for future genomic research in T. bryosalmonae and represents an important step towards the development of effective drugs against PKD.