Dual-Organ Transcriptomic Analysis of Rainbow Trout Infected With Ichthyophthirius multifiliis Through Co-Expression and Machine Learning

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.

Transcriptome analysis revealed immune responses in the kidney of Atlantic salmon (Salmo salar) co-infected with sea lice (Lepeophtheirus salmonis) and infectious salmon anemia virus

Sea lice (Lepeophtheirus salmonis) and infectious salmon anemia virus (ISAv) are two of the most important pathogens in Atlantic salmon (Salmo salar) farming and typically cause substantial economic losses to the industry. However, the immune interactions between hosts and these pathogens are still unclear, especially in the scenario of co-infection. In this study, we artificially infected Atlantic salmon with sea lice and ISAv, and investigated the gene expression patterns of Atlantic salmon head kidneys in response to both lice only and co-infection with lice and ISAv by transcriptomic analysis. The challenge experiment indicated that co-infection resulted in a cumulative mortality rate of 47.8 %, while no mortality was observed in the lice alone infection. We identified 240 differentially expressed genes (DEGs) under the lice alone infection, of which 185 were down-regulated and 55 were up-regulated, while a total of 994 DEGs were identified in the co-infection, of which 206 were down-regulated and 788 were significantly up-regulated. The pathway enrichment analysis revealed that single-infection significantly suppressed the innate immune system (e.g., the complement system), whereas co-infection induced a strong immune response, leading to the activation of immune-related signaling pathways such as Toll-like receptors and NOD-like receptors pathways, as well as significant upregulation of genes related to the activation of interferon and MH class I protein complex. Our results provide the first global transcriptomic study of gene expression in the Atlantic salmon head kidney in response to co-infection with sea lice and ISAv, and provided the baseline knowledge for understanding the immune responses during co-infection.

Identification of Genes Related to Resistance to Ichthyophthirius multifiliis Based on Co-expression Network Analysis in Grass Carp

The ciliate protozoan Ichthyophthirius multifiliis is an essential parasite causing white spot disease in grass carp, leading to significant economic losses. Understanding the molecular basis of grass carp's response to I. multifiliis has important scientific and environmental values. The transcriptional network analysis offers a valuable strategy to decipher the changes in gene expression in grass carp infected with I. multifiliis. Our goal was to screen the genes and pathways involved in resistance to I. multifiliis in grass carp. The different traits exhibited by grass carp infected with I. multifiliis may be caused by the differences in gene expression among grass carp individuals. Herein, to reveal those resistance-associated genes against I. multifiliis infection, we performed RNA sequencing using weighted gene co-expression network analysis (WGCNA). The biological function analysis and hub gene annotation for highly relevant modules revealed that different pathogen recognition and clearance responses resulted in different resistance to I. multifiliis infection. Furthermore, gene enrichment analysis revealed that I. multifiliis invasion in the disease-resistant group mainly activated immune pathways, including scavenger receptor activity and kappa B kinase/NF-kappa B signaling. By the annotation of the highly correlated module of the hub gene, we revealed that the apoptosis and ribosome biogenesis-related genes were enriched in the disease-resistant grass carp. The results of the dark grey module showed that several genes were mainly enriched in the two-component system (ko02020) and steroid biosynthesis (ko00100), suggesting that they are resistance-associated and energy metabolism-associated genes. In the disease resistance group, hub genes mainly included Nlrc3, fos, AAP8, HAP2, HAX, cho2, and zgc:113,036. This study revealed the gene network associated with disease resistance after I. multifiliis infection. The disease resistance-related pathways and central genes identified in this study are candidate references for breeders breeding disease-resistant. The results of this study may also provide some references for the development of drugs to antagonize I. multifiliis infection.

Diapause-Linked Gene Expression Pattern and Related Candidate Duplicated Genes of the Mountain Butterfly Parnassius glacialis (Lepidoptera: Papilionidae) Revealed by Comprehensive Transcriptome Profiling

The mountain butterfly Parnassius glacialis is a representative species of the genus Parnassius, which probably originated in the high-altitude Qinhai-Tibet Plateau in the Miocene and later dispersed eastward into relatively low-altitude regions of central to eastern China. However, little is known about the molecular mechanisms underlying the long-term evolutionary adaptation to heterogeneous environmental conditions of this butterfly species. In this study, we obtained the high-throughput RNA-Seq data from twenty-four adult individuals in eight localities, covering nearly all known distributional areas in China, and firstly identified the diapause-linked gene expression pattern that is likely to correlate with local adaptation in adult P. glacialis populations. Secondly, we found a series of pathways responsible for hormone biosynthesis, energy metabolism and immune defense that also exhibited unique enrichment patterns in each group that are probably related to habitat-specific adaptability. Furthermore, we also identified a suite of duplicated genes (including two transposable elements) that are mostly co-expressed to promote the plastic responses to different environmental conditions. Together, these findings can help us to better understand this species' successful colonization to distinct geographic areas from the western to eastern areas of China, and also provide us with some insights into the evolution of diapause in mountain Parnassius butterfly species.

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.

Social complexity, life-history and lineage influence the molecular basis of castes in vespid wasps

A key mechanistic hypothesis for the evolution of division of labour in social insects is that a shared set of genes co-opted from a common solitary ancestral ground plan (a genetic toolkit for sociality) regulates caste differentiation across levels of social complexity. Using brain transcriptome data from nine species of vespid wasps, we test for overlap in differentially expressed caste genes and use machine learning models to predict castes using different gene sets. We find evidence of a shared genetic toolkit across species representing different levels of social complexity. We also find evidence of additional fine-scale differences in predictive gene sets, functional enrichment and rates of gene evolution that are related to level of social complexity, lineage and of colony founding. These results suggest that the concept of a shared genetic toolkit for sociality may be too simplistic to fully describe the process of the major transition to sociality.

Differentially expressed genes in head kidney of Pelteobagrus fulvidraco following Vibrio cholerae challenge

The yellow catfish (Pelteobagrus fulvidraco) is a freshwater fish with high economic value in eastern China. Nevertheless, pathogens causing bacterial diseases in P. fulvidraco have brought about huge economic loss and high mortality in artificial aquaculture. For disease control, it is critical to further understand the immune system of yellow catfish and immune-related genes with which they respond to pathogenic infections. In this study, high-throughput sequencing methods were used to analyze the transcriptomic spectrum of the head kidney from P. fulvidraco challenged by Vibrio cholera. A total of 45,544 unique transcript fragments (unigenes) were acquired after assembly and annotation, with an average length of 1,373 bp. Additionally, 674 differentially expressed genes (DEGs) were identified after stimulation with V. cholerae, 353 and 321 genes were identified as remarkably up- or downregulated, respectively. To further study the immune-related DEGs, we performed KEGG enrichment and GO enrichment. The results showed gene regulation of response to stimulus, immune response, immune system progress, response to external stimuli and cellular response to stimuli. Analysis of KEGG enrichment is important to identify chief immune related pathways. Real-time quantitative reverse transcription-PCR (qRT-PCR) results indicated 10 immune response genes that were found to be upregulated compared to a control group after 6 h of V. cholerae challenging. In summary, the results of our study are helpful to determine the defense mechanisms and immune system responses of yellow catfish in reaction to bacterial challenges.

Multi-Organ Transcriptome Response of Lumpfish (Cyclopterus lumpus) to Aeromonas salmonicida Subspecies salmonicida Systemic Infection

Lumpfish is utilized as a cleaner fish to biocontrol sealice infestations in Atlantic salmon farms. Aeromonas salmonicida, a Gram-negative facultative intracellular pathogen, is the causative agent of furunculosis in several fish species, including lumpfish. In this study, lumpfish were intraperitoneally injected with different doses of A. salmonicida to calculate the LD50. Samples of blood, head-kidney, spleen, and liver were collected at different time points to determine the infection kinetics. We determined that A. salmonicida LD50 is 102 CFU per dose. We found that the lumpfish head-kidney is the primary target organ of A. salmonicida. Triplicate biological samples were collected from head-kidney, spleen, and liver pre-infection and at 3- and 10-days post-infection for RNA-sequencing. The reference genome-guided transcriptome assembly resulted in 6246 differentially expressed genes. The de novo assembly resulted in 403,204 transcripts, which added 1307 novel genes not identified by the reference genome-guided transcriptome. Differential gene expression and gene ontology enrichment analyses suggested that A. salmonicida induces lethal infection in lumpfish by uncontrolled and detrimental blood coagulation, complement activation, inflammation, DNA damage, suppression of the adaptive immune system, and prevention of cytoskeleton formation.

Identification and characterization of warm temperature acclimation proteins (Wap65s) in rainbow trout (Oncorhynchus mykiss)

Hemopexin is a vital glycoprotein for processing excessive iron in blood and functions as an iron scavenger in mammals. Teleosts however, unlike mammals, have two known hemopexin paralogs called warm temperature acclimation-related 65 kDa protein (Wap65-1 and Wap65-2, collectively termed Wap65s). Although Wap65s in rainbow trout have been considered notable biomarkers with significantly higher and/or lower expression under conditions of stress or disease, the individual roles, similarities and differences between the two paralogs are not well known. The aim of this study was to gain an understanding of the characteristics and functions of trout Wap65s from the perspective of iron-metabolism, physiological roles, and relevant immunological responses. The expression of Wap65-1 and -2 in this study was determined in the face of challenges by Aeromonas salmonicida, infectious hematopoietic necrosis virus (IHNV), and iron-dextran. Immuno-histochemistry (IHC) was employed to localize the major cell types for Wap65-2 expression, and trout leukocytes were isolated and incubated with LPS and OxLDL for comprehending the immunological characteristics of Wap65-2. We demonstrate that Wap65-1 is expressed only in the liver but Wap65-2 is systemically expressed in most organs and tissues. Interestingly, Wap65-1 expression was not significantly changed under A. salmonicida and iron-dextran administration, but was significantly decreased under IHNV. In contrast, Wap65-2 was up-regulated in all challenged groups, however with different expression patterns in the blood and liver. These results suggested that the two paralogs may participate in different biological roles. IHC showed that Wap65-2 antibody had high affinity for leukocyte-like cells, and macrophages but not lymphocytes significantly increased expression under LPS and OxLDL stimulation. These results support the conclusion that trout Wap65-2, not Wap65-1 may have conventional hemopexin functions such as reported in mammals including effects on iron metabolism, inflammation, and acute-phase protein.

Identification, classification and functional characterization of HSP70s in rainbow trout (Oncorhynchus mykiss) through multi-omics approaches

Heat shock protein 70s (HSP70s) are known to play vital biological processes in rainbow trout. However, information on the numerous roles and classification of many different HSP70s is insufficient. The purpose of this study was to investigate the characteristics of all HSP70s in rainbow trout using multi-dimensional genomic and transcriptomic analyses for inspecting HSP70 homologs, phylogenetic characteristics, DNA motifs, and transcription factor binding sites (TFBSs). Also, the transcriptomic results in conditions of acute thermal stress and Ichthyophthirius multifiliis infection were used to characterize the expression of all HSP70 homologs, and the isoforms of the most sensitive HSP70 were predicted in silico. A total of 23 HSP70s were identified, and they were divided into seven evolutionary groups (groups 1-7). Groups 1 and 2 had relatively longer phylogenetic distances compared to the other groups, which can speculate origin of groups 1 and 2 HSP70s would be different compared to others. With transcriptomic profiling, most HSPs belonging to group 3 showed highly sensitive responses to I. multifiliis infection, not thermal stress, but the group 6 HSP70s had the opposite expression tendencies. Likewise, the composition of the TFBS in each HSP70 was consistent with its group classification. Since TFBSs are widely known to influence transcriptomic expression, they could be one of the major reasons for the different patterns of expression within the HSP70 groups. Moreover, this study demonstrated several isoforms of HSP70a, by far the most sensitive HSP70s, under several stress environments such as hypoxia, thermal, and overcrowding stress. This is an important fundamental study to expand the understanding of HSP70s in rainbow trout as well as for selecting the most sensitive biomarkers for types of stress.