A tale of two fish: Comparative transcriptomics of resistant and susceptible steelhead following exposure to Ceratonova shasta highlights differences in parasite recognition

Subtle transcriptomic response of Eurasian perch (Perca fluviatilis) associated with Triaenophorus nodulosus plerocercoid infection

Determining the physiological effects of parasites and characterizing genes involved in host responses to infections are essential to improving our understanding of host-parasite interactions and their ecological and evolutionary consequences. This task, however, is complicated by high diversity and complex life histories of many parasite species. The use of transcriptomics in the context of wild-caught specimens can help ameliorate this by providing both qualitative and quantitative information on gene expression patterns in response to parasites in specific host organs and tissues. Here, we evaluated the physiological impact of the widespread parasite, the pike tapeworm (Triaenophorus nodulosus), on its second intermediate host, the Eurasian perch (Perca fluviatilis). We used an RNAseq approach to analyse gene expression in the liver, the target organ of T. nodulosus plerocercoids, and spleen which is one of the main immune organs in teleost fishes. We compared perch collected from multiple lakes consisting of individuals with (n = 8) and without (n = 6) T. nodulosus plerocercoids in the liver. Results revealed a small number of differentially expressed genes (DEGs, adjusted p-value ≤0.05) in both spleen (n = 22) and liver (n = 10). DEGs in spleen consisted of mostly upregulated immune related genes (e.g., JUN, SIK1, THSB1), while those in the liver were often linked to metabolic functions (e.g., FABP1, CADM4, CDAB). However, Gene Ontology (GO) analysis showed lack of functional enrichment among DEGs. This study demonstrates that Eurasian perch displays a subtle response at a gene expression level to T. nodulosus plerocercoid infection. Given that plerocercoids are low-metabolic activity transmission stages, our results suggest that moderate T. nodulosus plerocercoid infection most likely does not provoke an extensive host immune response and have relatively low physiological costs for the host. Our findings illustrate that not all conspicuous infections have severe effects on host gene regulation.

Plasma proteomic profiling of bacterial cold water disease-resistant and -susceptible rainbow trout lines and biomarker discovery

Genetic variation for disease resistance is present in salmonid fish; however, the molecular basis is poorly understood, and biomarkers of disease susceptibility/resistance are unavailable. Previously, we selected a line of rainbow trout for high survival following standardized challenge with Flavobacterium psychrophilum (Fp), the causative agent of bacterial cold water disease. The resistant line (ARS-Fp-R) exhibits over 60 percentage points higher survival compared to a reference susceptible line (ARS-Fp-S). To gain insight into the differential host response between genetic lines, we compared the plasma proteomes from day 6 after intramuscular challenge. Pooled plasma from unhandled, PBS-injected, and Fp-injected groups were simultaneously analyzed using a TMT 6-plex label, and the relative abundance of 513 proteins was determined. Data are available via ProteomeXchange, with identifier PXD041308, and the relative protein abundance values were compared to mRNA measured from a prior, whole-body RNA-seq dataset. Our results identified a subset of differentially abundant intracellular proteins was identified, including troponin and myosin, which were not transcriptionally regulated, suggesting that these proteins were released into plasma following pathogen-induced tissue damage. A separate subset of high-abundance, secreted proteins were transcriptionally regulated in infected fish. The highest differentially expressed protein was a C1q family member (designated complement C1q-like protein 3; C1q-LP3) that was upregulated over 20-fold in the infected susceptible line while only modestly upregulated, 1.8-fold, in the infected resistant line. Validation of biomarkers was performed using immunoassays and C1q-LP3, skeletal muscle troponin C, cathelcidin 2, haptoglobin, leptin, and growth and differentiation factor 15 exhibited elevated concentration in susceptible line plasma. Complement factor H-like 1 exhibited higher abundance in the resistant line compared to the susceptible line in both control and challenged fish and thus was a baseline differentiator between lines. C1q-LP3 and STNC were elevated in Atlantic salmon plasma following experimental challenge with Fp. In summary, these findings further the understanding of the differential host response to Fp and identifies salmonid biomarkers that may have use for genetic line evaluation and on-farm health monitoring.

Infectious agents and their physiological correlates in early marine Chinook salmon (Oncorhynchus tshawytscha)

The early marine life of Pacific salmon is believed to be a critical period limiting population-level survival. Recent evidence suggests that some infectious agents are associated with survival but linkages with underlying physiological mechanisms are lacking. While challenge studies can demonstrate cause and effect relationships between infection and pathological change or mortality, in some cases pathological change may only manifest in the presence of environmental stressors; thus, it is important to gain context from field observations. Herein, we examined physiological correlates with infectious agent loads in Chinook salmon during their first ocean year. We measured physiology at the molecular (gene expression), metabolic (plasma chemistry) and cellular (histopathology) levels. Of 46 assayed infectious agents, 27 were detected, including viruses, bacteria and parasites. This exploratory study identified.

a strong molecular response to viral disease and pathological change consistent with jaundice/anemia associated with Piscine orthoreovirus,strong molecular signals of gill inflammation and immune response associated with gill agents `Candidatus Branchiomonas cysticola' and Parvicapsula pseudobranchicola,a general downregulation of gill immune response associated with Parvicapsula minibicornis complementary to that of P. pseudobranchicola.Importantly, our study provides the first evidence that the molecular activation of viral disease response and the lesions observed during the development of the PRV-related disease jaundice/anemia in farmed Chinook salmon are also observed in wild juvenile Chinook salmon.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Ceratonova shasta: a cnidarian parasite of annelids and salmonids

The myxozoan Ceratonova shasta was described from hatchery rainbow trout over 70 years ago. The parasite continues to cause severe disease in salmon and trout, and is recognized as a barrier to salmon recovery in some rivers. This review incorporates changes in our knowledge of the parasite's life cycle, taxonomy and biology and examines how this information has expanded our understanding of the interactions between C. shasta and its salmonid and annelid hosts, and how overarching environmental factors affect this host–parasite system. Development of molecular diagnostic techniques has allowed discrimination of differences in parasite genotypes, which have differing host affinities, and enabled the measurement of the spatio-temporal abundance of these different genotypes. Establishment of the C. shasta life cycle in the laboratory has enabled studies on host–parasite interactions and the availability of transcriptomic data has informed our understanding of parasite virulence factors and host defences. Together, these advances have informed the development of models and management actions to mitigate disease.

Infection by the Parasite Myxobolus bejeranoi (Cnidaria: Myxozoa) Suppresses the Immune System of Hybrid Tilapia

Myxozoa (Cnidaria) is a large group of microscopic obligate endoparasites that can cause emerging diseases, affecting wild fish populations and fisheries. Recently, the myxozoan Myxobolus bejeranoi was found to infect the gills of hybrid tilapia (Nile tilapia (Oreochromis niloticus) × Jordan/blue tilapia (O. aureus)), causing high morbidity and mortality. Here, we used comparative transcriptomics to elucidate the molecular processes occurring in the fish host following infection by M. bejeranoi. Fish were exposed to pond water containing actinospores for 24 h and the effects of minor, intermediate, and severe infections on the sporulation site, the gills, and on the hematopoietic organs, head kidney and spleen, were compared. Enrichment analysis for GO and KEGG pathways indicated immune system activation in gills at severe infection, whereas in the head kidney a broad immune suppression included deactivation of cytokines and GATA3 transcription factor responsible for T helper cell differentiation. In the spleen, the cytotoxic effector proteins perforin and granzyme B were downregulated and insulin, which may function as an immunomodulatory hormone inducing systemic immune suppression, was upregulated. These findings suggest that M. bejeranoi is a highly efficient parasite that disables the defense mechanisms of its fish host hybrid tilapia.

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.

Severe Natural Outbreak of Cryptocaryon irritans in Gilthead Seabream Produces Leukocyte Mobilization and Innate Immunity at the Gill Tissue

The protozoan parasite Cryptocaryon irritans causes marine white spot disease in a wide range of fish hosts, including gilthead seabream, a very sensitive species with great economic importance in the Mediterranean area. Thus, we aimed to evaluate the immunity of gilthead seabream after a severe natural outbreak of C. irritans. Morphological alterations and immune cell appearance in the gills were studied by light microscopy and immunohistochemical staining. The expression of several immune-related genes in the gills and head kidney were studied by qPCR, including inflammatory and immune cell markers, antimicrobial peptides (AMP), and cell-mediated cytotoxicity (CMC) molecules. Serum humoral innate immune activities were also assayed. Fish mortality reached 100% 8 days after the appearance of the C. irritans episode. Gill filaments were engrossed and packed without any space between filaments and included parasites and large numbers of undifferentiated and immune cells, namely acidophilic granulocytes. Our data suggest leukocyte mobilization from the head kidney, while the gills show the up-regulated transcription of inflammatory, AMPs, and CMC-related molecules. Meanwhile, only serum bactericidal activity was increased upon infection. A potent local innate immune response in the gills, probably orchestrated by AMPs and CMC, is triggered by a severe natural outbreak of C. irritans.

The cnidarian parasite Ceratonova shasta utilizes inherited and recruited venom-like compounds during infection

Background

Cnidarians are the most ancient venomous organisms. They store a cocktail of venom proteins inside unique stinging organelles called nematocysts. When a cnidarian encounters chemical and physical cues from a potential threat or prey animal, the nematocyst is triggered and fires a harpoon-like tubule to penetrate and inject venom into the prey. Nematocysts are present in all Cnidaria, including the morphologically simple Myxozoa, which are a speciose group of microscopic, spore-forming, obligate parasites of fish and invertebrates. Rather than predation or defense, myxozoans use nematocysts for adhesion to hosts, but the involvement of venom in this process is poorly understood. Recent work shows some myxozoans have a reduced repertoire of venom-like compounds (VLCs) relative to free-living cnidarians, however the function of these proteins is not known.

Methods

We searched for VLCs in the nematocyst proteome and a time-series infection transcriptome of Ceratonova shasta, a myxozoan parasite of salmonid fish. We used four parallel approaches to detect VLCs: BLAST and HMMER searches to preexisting cnidarian venom datasets, the machine learning tool ToxClassifier, and structural modeling of nematocyst proteomes. Sequences that scored positive by at least three methods were considered VLCs. We then mapped their time-series expressions in the fish host and analyzed their phylogenetic relatedness to sequences from other venomous animals.

Results

We identified eight VLCs, all of which have closely related sequences in other myxozoan datasets, suggesting a conserved venom profile across Myxozoa, and an overall reduction in venom diversity relative to free-living cnidarians. Expression of the VLCs over the 3-week fish infection varied considerably: three sequences were most expressed at one day post-exposure in the fish’s gills; whereas expression of the other five VLCs peaked at 21 days post-exposure in the intestines, coinciding with the formation of mature parasite spores with nematocysts. Expression of VLC genes early in infection, prior to the development of nematocysts, suggests venoms in C. shasta have been repurposed to facilitate parasite invasion and proliferation within the host. Molecular phylogenetics suggested some VLCs were inherited from a cnidarian ancestor, whereas others were more closely related to sequences from venomous non-Cnidarian organisms and thus may have gained qualities of venom components via convergent evolution. The presence of VLCs and their differential expression during parasite infection enrich the concept of what functions a “venom” can have and represent targets for designing therapeutics against myxozoan infections.

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.

Intestinal Transcriptomic and Histologic Profiling Reveals Tissue Repair Mechanisms Underlying Resistance to the Parasite Ceratonova shasta

BACKGROUND

Myxozoan parasites infect fish worldwide causing significant disease or death in many economically important fish species, including rainbow trout and steelhead trout (Oncorhynchus mykiss). The myxozoan Ceratonova shasta is a parasite of salmon and trout that causes ceratomyxosis, a disease characterized by severe inflammation in the intestine resulting in hemorrhaging and necrosis. Populations of O. mykiss that are genetically fixed for resistance or susceptibility to ceratomyxosis exist naturally, offering a tractable system for studying the immune response to myxozoans. The aim of this study was to understand how steelhead trout that are resistant to the disease respond to C. shasta once it has become established in the intestine and identify potential mechanisms of resistance.

RESULTS

Sequencing of intestinal mRNA from resistant steelhead trout with severe C. shasta infections identified 417 genes differentially expressed during the initial stage of the infection compared to uninfected control fish. A strong induction of interferon-gamma and interferon-stimulated genes was evident, along with genes involved in cell adhesion and migration. A total of 11,984 genes were differentially expressed during the late stage of the infection, most notably interferon-gamma, interleukin-6, and immunoglobulin transcripts. A distinct hardening of the intestinal tissue and a strong inflammatory reaction in the intestinal submucosa including severe hyperplasia and inflammatory cell infiltrates were observed in response to the infection. The massive upregulation of caspase-14 early in the infection, a protein involved in keratinocyte differentiation might reflect the rapid onset of epithelial repair mechanisms, and the collagenous stratum compactum seemed to limit the spread of C. shasta within the intestinal layers. These observations could explain the ability of resistant fish to eventually recover from the infection.

CONCLUSIONS

Our results suggest that resistance to ceratomyxosis involves both a rapid induction of key immune factors and a tissue response that limits the spread of the parasite and the subsequent tissue damage. These results improve our understanding of the myxozoan-host dialogue and provide a framework for future studies investigating the infection dynamics of C. shasta and other myxozoans.