Transcriptome profiling of posterior kidney of brown trout, Salmo trutta, during proliferative kidney disease

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.

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.

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.

Potential biomarkers of endocrine and habitat disruption identified via RNA-Seq in Salvelinus fontinalis with proximity to fracking operations in Pennsylvania headwater stream ecosystems

Unconventional natural gas development (fracking) has been a rapidly expanding technique used for the extraction of natural gas from the Marcellus Shale formation in Pennsylvania. There remains a knowledge gap regarding the ecological impacts of fracking, especially regarding the long-term health of native Brook trout (Salvelinus fontinalis) populations. During the summer of 2015, Brook trout were sampled from twelve streams located in forested, northwestern Pennsylvania in order to evaluate the impacts of fracking on Brook trout. Four stream sites were undisturbed (no fracking activity), three had a developed well pad without fracking activity, and five had active fracking with natural gas production. Liver tissue was isolated from two to five fish per stream and underwent RNA-Seq analysis to identify differentially expressed genes between ecosystems with differing fracking status. Data were analyzed individually and with samples pooled within-stream to account for hierarchical data structure and variation in sample coverage within streams. Differentially expressed and differentially alternatively spliced genes had functions related to lipid and steroid metabolism, mRNA processing, RNA polymerase and protein regulation. Unique to our study, genes related to xenobiotic and stress responses were found as well as potential markers for endocrine disruption and saline adaptation that were identified in watersheds with active fracking activity. These results support the utility of RNA-Seq to assess trout health and suggest detrimental impacts of fracking on sensitive trout populations.

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.

Kinetics of Parasite-Specific Antibody and B-Cell-Associated Gene Expression in Brown Trout, Salmo trutta during Proliferative Kidney Disease

Tetracapsuloides bryosalmonae, a myxozoan endoparasite often causes chronic infection in brown trout. Antiparasite immunity mediated by antibodies and B cells is known as an important determinant of host survival and parasite proliferation during chronic infections. Accordingly, studying their time course during proliferative kidney disease (PKD) might be helpful in improving our understanding of its chronic nature. Therefore, we conducted this study to examine parasite specific serum antibody and B-cell-mediated response in laboratory-infected brown trout at different time points. Brown trout were exposed to the spores of T. bryosalmonae, derived from infected bryozoans. Samples were collected at different time points and processed for indirect ELISA, histopathology, and qRT-PCR. T. bryosalmonae specific antibody was detected at 4 weeks post exposure (wpe) and it persisted until 17 wpe. Additionally, the expressions of C4A, CD34, CD79A, BLNK, CD74, BCL7, and CD22 were differentially regulated in the important immune organs, kidney and spleen. To our knowledge, this is the first study addressing anti-T. bryosalmonae antibody response in brown trout at different time points. The results from this study provide valuable insights into the processes leading to changes in B cell development, inflammation and antibody production during the course of PKD in brown trout.

RNA-Seq-Based Analysis of Cortisol-Induced Differential Gene Expression Associated with Piscirickettsia salmonis Infection in Rainbow Trout (Oncorhynchus mykiss) Myotubes

Salmonid rickettsial septicemia (SRS) is the major infectious disease of the Chilean salmonid aquaculture industry caused by Piscirickettsia salmonis. Intensive farming conditions generate stress and increased susceptibility to diseases, being skeletal muscle mainly affected. However, the interplay between pathogen infection and stress in muscle is poorly understood. In this study, we perform an RNA-seq analysis on rainbow trout myotubes that are pretreated for 3 h with cortisol (100 ng/mL) and then infected with P. salmonis strain LF-89 for 8 h (MOI 50). Twelve libraries are constructed from RNA samples (n = 3 per group) and sequenced on Illumina HiSeq 4000. A total of 704,979,454 high-quality reads are obtained, with 70.25% mapped against the reference genome. In silico DETs include 175 total genes-124 are upregulated and 51 are downregulated. GO enrichment analysis reveals highly impacted biological processes related to apoptosis, negative regulation of cell proliferation, and innate immune response. These results are validated by RT-qPCR of nine candidate transcripts. Furthermore, cortisol pretreatment significantly stimulated bacterial gene expression of ahpC and 23s compared to infection. In conclusion, for the first time, we describe a transcriptomic response of trout myotubes infected with P. salmonis by inducing apoptosis, downregulating cell proliferation, and intrinsic immune-like response that is differentially regulated by cortisol.

It's a hard knock life for some: Heterogeneity in infection life history of salmonids influences parasite disease outcomes

Heterogeneity in immunity occurs across numerous disease systems with individuals from the same population having diverse disease outcomes. Proliferative kidney disease (PKD) caused by Tetracapsuloides bryosalmonae, is a persistent parasitic disease negatively impacting both wild and farmed salmonids. Little is known of how PKD is spread or maintained within wild susceptible populations. We investigated an aspect of fish disease that has been largely overlooked, that is, the role of the host phenotypic heterogeneity in disease outcome. We examined how host susceptibility to T. bryosalmonae infection, and the disease PKD, varied across different infection life-history stages and how it differs between naïve, re-infected and persistently infected hosts. We investigated the response to parasite exposure in host phenotypes with (a) different ages and (b) heterogeneous infection life histories. Among (a) the age phenotypes were young-of-the-year (YOY) fish and juvenile 1+ fish (fish older than one) and, for (b) juvenile 1+ infection survivors were either re-exposed or not re- exposed to the parasite and response phenotypes were assigned post-hoc dependant on infection status. In fish not re-exposed this included fish that cleared infection (CI) or had a persistent infection (PI). In fish re-exposed these included fish that were re-infected (RI), or re-exposed and uninfected (RCI). We assessed both parasite-centric (infection prevalence, parasite burden, malacospore transmission) and host-centric parameters (growth rates, disease severity, infection tolerance and the immune response). In (a), YOY fish, parasite success and disease severity were greater and differences in the immune response occurred, demonstrating an ontogenetic decline of susceptibility in older fish. In (b), in PI and RI fish, parasite success and disease severity were comparable. However, expression of several adaptive immunity markers was greater in RI fish, indicating concomitant immunity, as re-exposure did not intensify infection. We demonstrate the relevance of heterogeneity in infection life history on disease outcome and describe several distinctive features of immune ontogeny and protective immunity in this model not previously reported. The relevance of such themes on a population level requires greater research in many aquatic disease systems to generate clearer framework for understanding the spread and maintenance of aquatic pathogens.

Development of Fish Parasite Vaccines in the OMICs Era: Progress and Opportunities

Globally, parasites are increasingly being recognized as catastrophic agents in both aquaculture sector and in the wild aquatic habitats leading to an estimated annual loss between 1.05 billion and 9.58 billion USD. The currently available therapeutic and control measures are accompanied by many limitations. Hence, vaccines are recommended as the "only green and effective solution" to address these concerns and protect fish from pathogens. However, vaccine development warrants a better understanding of host-parasite interaction and parasite biology. Currently, only one commercial parasite vaccine is available against the ectoparasite sea lice. Additionally, only a few trials have reported potential vaccine candidates against endoparasites. Transcriptome, genome, and proteomic data at present are available only for a limited number of aquatic parasites. Omics-based interventions can be significant in the identification of suitable vaccine candidates, finally leading to the development of multivalent vaccines for significant protection against parasitic infections in fish. The present review highlights the progress in the immunobiology of pathogenic parasites and the prospects of vaccine development. Finally, an approach for developing a multivalent vaccine for parasitic diseases is presented. Data sources to prepare this review included Pubmed, google scholar, official reports, and websites.

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.

Modulation of local and systemic immune responses in brown trout (Salmo trutta) following exposure to Myxobolus cerebralis

Myxobolus cerebralis, the etiological agent of Whirling Disease (WD), is a freshwater myxozoan parasite with considerable economic and ecological relevance for salmonids. There are differences in disease susceptibility between species and strains of salmonids. Recently, we have reported that the suppressor of cytokine signaling SOCS1 and SOCS3 are key in modulating rainbow trout (Oncorhynchus mykiss) immune responses and that resistant fish apparently exhibit effective Th17 cell response after exposure to M. cerebralis. It is unclear whether such molecules and pathways are also involved in the immune response of M. cerebralis infected brown trout (Salmo trutta). Hence, this study aimed to explore their role during immune modulation in infected brown trout, which is considered resistant to this parasite. Fish were exposed to the triactinomyxon (TAM) stages of M. cerebralis and quantitative real-time PCR (RT-qPCR) was carried out to examine local (caudal fin) and systemic (head kidney, spleen) immune transcriptional changes associated with WD over time in infected and control fish. All of the immune genes in the three tissues studied were differentially expressed in infected fish at multiple time points. Brown trout reduced the parasite load and demonstrated effective immune responses, likely by keeping pro-inflammatory and anti-inflammatory cytokines in balance whilst stimulating efficient Th17-mediated immunity. This study increases knowledge on the brown trout immune response to M. cerebralis and helps us to understand the underlying mechanisms of WD resistance.

Transcriptome analysis of goldfish (Carassius auratus) in response to Gyrodactylus kobayashii infection

Gyrodactylid monogeneans are widespread parasites of teleost fishes, and infection with these parasites results in high host morbidity and mortality in aquaculture. To comprehensively elucidate the immune mechanisms against Gyrodactylus kobayashii, the transcriptome profiles of goldfish (Carassius auratus) skin after challenge with G. kobayashii were first investigated using next-generation sequencing. Approximately 21 million clean reads per library were obtained, and the average percentage of these clean reads mapped to the reference genome was 82.25%. A total of 556 differentially expressed genes (DEGs), including 344 upregulated and 212 downregulated genes, were identified, and 380 DEGs were successfully annotated and assigned to 95 signaling pathways in Kyoto Encyclopedia of Genes and Genomes (KEGG). In addition, 14 pathways associated with immune response were identified mainly including mTOR signaling pathway, cytokine-cytokine receptor interaction, intestinal immune network for IgA production, toll-like receptor signaling pathway, and phagosome. Twelve genes were selected and validated using qRT-PCR. A similar trend of these genes between RNA-Seq and qRT-PCR was observed, indicating that RNA-Seq data was reliable. Besides, the ALP activity and NO content in serum were significantly higher in the infected goldfish compared with the non-infected goldfish. In summary, this study provides better understandings of immune defense mechanisms of goldfish against G. kobayashii, which will support future molecular research on gyrodactylids and facilitate the prevention and treatment of gyrodactylosis in aquaculture.

Transcriptome Analysis Elucidates the Key Responses of Bryozoan Fredericella sultana during the Development of Tetracapsuloides bryosalmonae (Myxozoa)

Bryozoans are sessile, filter-feeding, and colony-building invertebrate organisms. Fredericella sultana is a well known primary host of the myxozoan parasite Tetracapsuloides bryosalmonae. There have been no attempts to identify the cellular responses induced in F. sultana during the T. bryosalmonae development. We therefore performed transcriptome analysis with the aim of identifying candidate genes and biological pathways of F. sultana involved in the response to T. bryosalmonae. A total of 1166 differentially up- and downregulated genes were identified in the infected F. sultana. Gene ontology of biological processes of upregulated genes pointed to the involvement of the innate immune response, establishment of protein localization, and ribosome biogenesis, while the downregulated genes were involved in mitotic spindle assembly, viral entry into the host cell, and response to nitric oxide. Eukaryotic Initiation Factor 2 signaling was identified as a top canonical pathway and MYCN as a top upstream regulator in the differentially expressed genes. Our study provides the first transcriptional profiling data on the F. sultana zooid's response to T. bryosalmonae. Pathways and upstream regulators help us to understand the complex interplay in the infected F. sultana. The results will facilitate the elucidation of innate immune mechanisms of bryozoan and will lay a foundation for further analyses on bryozoan-responsive candidate genes, which will be an important resource for the comparative analysis of gene expression in bryozoans.

Identification and Expression Profiling of Toll-Like Receptors of Brown Trout (Salmo trutta) during Proliferative Kidney Disease

Proliferative kidney disease is an emerging disease among salmonids in Europe and North America caused by the myxozoan parasite Tetracapsuloides bryosalmonae. The decline of endemic brown trout (Salmo trutta) in the Alpine streams of Europe is fostered by T. bryosalmonae infection. Toll-like receptors (TLRs) are a family of pattern recognition receptors that acts as sentinels of the immune system against the invading pathogens. However, little is known about the TLRs' response in salmonids against the myxozoan infection. In the present study, we identified and evaluated TLR1, TLR19, and TLR13-like genes of brown trout using data-mining and phylogenetic analysis. The expression pattern of TLRs was examined in the posterior kidney of brown trout infected with T. bryosalmonae at various time points. Typical Toll/interleukin-1 receptor protein domain was found in all tested TLRs. However, TLR13-like chr2 had a short amino acid sequence with no LRR domain. Phylogenetic analysis illustrated that TLR orthologs are conserved across vertebrates. Similarly, a conserved synteny gene block arrangement was observed in the case of TLR1 and TLR19 across fish species. Interestingly, all tested TLRs showed their maximal relative expression from 6 to 10 weeks post-exposure to the parasite. Our results suggest that these TLRs may play an important role in the innate defense mechanism of brown trout against the invading T. bryosalmonae.

The Malacosporean Myxozoan Parasite Tetracapsuloides bryosalmonae: A Threat to Wild Salmonids

Tetracapsuloides bryosalmonae is a myxozoan parasite responsible for proliferative kidney disease (PKD) in a wide range of salmonids. PKD, characterized by high mortality and morbidity, is well known for affecting aquaculture operations and wild salmonid populations across Europe and North America. The life cycle of T. bryosalmonae revolves around freshwater bryozoan and salmonid fish hosts. In recent years, T. bryosalmonae has been reported among wild salmonids from the European countries where it has not been reported previously. T. bryosalmonae is believed to be a possible reason for the diminishing wild salmonid populations in the natural water bodies of many European countries. Climate crisis driven rising water temperature can further accelerate the distribution of T. bryosalmonae. Expansion of the geographical distribution of T. bryosalmonae may further advocate the decline of wild salmonid populations, especially brown trout (Salmo trutta) in their habitats. Mathematical models are used to understand the pattern and distribution of T. bryosalmonae among the host in the natural water bodies. The present manuscript not only summarizes the incidences of T. bryosalmonae among the wild salmonid populations, but also discusses the contemporary understanding about the development of T. bryosalmonae in its hosts and the influences of various factors in the spread of the disease in the wild.