Genetic and transcriptomic analyses provide new insights on the early antiviral response to VHSV in resistant and susceptible rainbow trout

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.

Immunological insights into the resistance of Nile tilapia strains to an infection with tilapia lake virus

The emergence of viral diseases affecting fish and causing very high mortality can lead to the disruption of aquaculture production. Recently, this occurred in Nile tilapia aquaculture where a disease caused by a systemic infection with a novel virus named tilapia lake virus (TiLV) caused havoc in cultured populations. With mortality surpassing 90% in young tilapia, the disease caused by TiLV has become a serious challenge for global tilapia aquaculture. In order to partly mitigate the losses, we explored the natural resistance to TiLV-induced disease in three genetic strains of tilapia which were kept at the University of Göttingen, Germany. We used two strains originating from Nilotic regions (Lake Mansala (MAN) and Lake Turkana (ELM)) and one from an unknown location (DRE). We were able to show that the virus is capable of overcoming the natural resistance of tilapia when injected, providing inaccurate mortality results that might complicate finding the resistant strains. Using the cohabitation infection model, we found an ELM strain that did not develop any clinical signs of the infection, which resulted in nearly 100% survival rate. The other two strains (DRE and MAN) showed severe clinical signs and much lower survival rates of 29.3% in the DRE strain and 6.7% in the MAN strain. The disease resistance of tilapia from the ELM strain was correlated with lower viral loads both at the mucosa and internal tissues. Our results suggest that the lower viral load could be caused by a higher magnitude of a mx1-based antiviral response in the initial phase of infection. The lower pro-inflammatory responses also found in the resistant strain might additionally contribute to its protection from developing pathological changes related to the disease. In conclusion, our results suggest the possibility of using TiLV-resistant strains as an ad hoc, cost-effective solution to the TiLV challenge. However, as the fish from the disease-resistant strain still retained significant virus loads in liver and brain and thus could become persistent virus carriers, they should be used within an integrative approach also combining biosecurity, diagnostics and vaccination measures.\.

Genome-Wide Identification and Expression Analysis of Potential Antiviral Tripartite Motif Proteins (TRIMs) in Grass Carp (Ctenopharyngodon idella)

Tripartite motif proteins (TRIMs), especially B30.2 domain-containing TRIMs (TRIMs-B30.2), are increasingly well known for their antiviral immune functions in mammals, while antiviral TRIMs are far from being identified in teleosts. In the present study, we identified a total of 42 CiTRIMs from the genome of grass carp, Ctenopharyngodon idella, an important cultured teleost in China, based on hmmsearch and SMART analysis. Among these CiTRIMs, the gene loci of 37 CiTRIMs were located on different chromosomes and shared gene collinearities with homologous counterparts from human and zebrafish genomes. They possessed intact conserved RBCC or RB domain assemblies at their N-termini and eight different domains, including the B30.2 domain, at their C-termini. A total of 19 TRIMs-B30.2 were identified, and most of them were clustered into a large branch of CiTRIMs in the dendrogram. Tissue expression analysis showed that 42 CiTRIMs were universally expressed in various grass carp tissues. A total of 11 significantly differentially expressed CiTRIMs were found in two sets of grass carp transcriptomes during grass carp reovirus (GCRV) infection. Three of them, including Cibtr40, CiTRIM103 and CiTRIM109, which all belonged to TRIMs-B30.2, were associated with the type I interferon response during GCRV infection by weighted network co-expression and gene expression trend analyses, suggesting their involvement in antiviral immunity. These findings may offer useful information for understanding the structure, evolution, and function of TRIMs in teleosts and provide potential antiviral immune molecule markers for grass carp.

Evidence of Transcriptional Shutoff by Pathogenic Viral Haemorrhagic Septicaemia Virus in Rainbow Trout

The basis of pathogenicity of viral haemorrhagic septicaemia virus (VHSV) was analysed in the transcriptome of a rainbow trout cell line inoculated with pathogenic and non-pathogenic VHSV isolates. Although both VHSV isolates showed similar viral replication patterns, the number of differentially expressed genes was 42-fold higher in cells inoculated with the non-pathogenic VHSV at 3 h post inoculation (hpi). Infection with the non-pathogenic isolate resulted in Gene Ontologies (GO) enrichment of terms such as immune response, cytokine-mediated signalling pathway, regulation of translational initiation, unfolded protein binding, and protein folding, and induced an over-representation of the p53, PPAR, and TGF-β signalling pathways. Inoculation with the pathogenic isolate resulted in the GO enrichment of terms related to lipid metabolism and the salmonella infection KEGG pathway involved in the rearrangement of the cytoskeleton. Antiviral response was evident at 12hpi in cells infected with the pathogenic isolate. Overall, the data showed a delay in the response of genes involved in immune responses and viral sensing in cells inoculated with the pathogenic isolate and suggest transcriptional shutoff and immune avoidance as a critical mechanism of pathogenicity in VHSV. These pathways offer opportunities to further understand and manage VHSV pathogenicity in rainbow trout.

Transcriptome profiling in head kidney of rainbow trout (Oncorhynchus mykiss) after infection with the low-virulent Nagano genotype of infectious hematopoietic necrosis virus

Infectious hematopoietic necrosis virus (IHNV) causes clinical diseases and mortality in a wide variety of salmonid species. Here, we studied transcriptional responses in rainbow trout infected by the IHNV-Nagano strain isolated in Korea. RNA-seq-based transcriptome analysis of head kidney tissues cataloged differentially expressed genes. Enrichment analysis of gene ontology annotations was performed, and a total of fifteen biological process terms were commonly identified at all time points. In the Kyoto Encyclopedia of Genes and Genomes pathway analysis, pathogen recognition receptor (PRR) signaling pathways such as the retinoic-acid-inducible gene-I-like receptor signaling pathway and the Toll-like receptor signaling pathway were identified at all time points. The nucleotide-binding oligomerization-domain-like receptor signaling pathway and cytosolic DNA-sensing pathway were identified at days 1 and 3. Protein-protein interaction network and centrality analyses revealed that the immune system, signaling molecules, and interaction pathways were upregulated at days 1 and 3, with the highest centrality of tumor necrosis factor. Cancer, cellular community, and endocrine system pathways were downregulated, with the highest centrality of fibronectin 1 at day 5. STAT1 was upregulated from days 1 to 5 with a high centrality. The reproducibility and repeatability of the transcriptome analysis were validated by RT-qPCR. IHNV-Nagano infection dynamically changed the transcriptome profiles in the head kidney of rainbow trout and induced a defense mechanism by regulating the immune and inflammatory pathways through PRR signaling at an early stage. Downregulated pathways involved in extracellular matrix formation and focal adhesion at day 5 indicated the possible failure of wound healing, which is important in the pathogenesis of IHNV infection.

Combining Multiple Approaches and Models to Dissect the Genetic Architecture of Resistance to Infections in Fish

Infectious diseases represent a major threat for the sustainable development of fish farming. Efficient vaccines are not available against all diseases, and growing antibiotics resistance limits the use of antimicrobial drugs in aquaculture. It is therefore important to understand the basis of fish natural resistance to infections to help genetic selection and to develop new approaches against infectious diseases. However, the identification of the main mechanisms determining the resistance or susceptibility of a host to a pathogenic microbe is challenging, integrating the complexity of the variation of host genetics, the variability of pathogens, and their capacity of fast evolution and adaptation. Multiple approaches have been used for this purpose: (i) genetic approaches, QTL (quantitative trait loci) mapping or GWAS (genome-wide association study) analysis, to dissect the genetic architecture of disease resistance, and (ii) transcriptomics and functional assays to link the genetic constitution of a fish to the molecular mechanisms involved in its interactions with pathogens. To date, many studies in a wide range of fish species have investigated the genetic determinism of resistance to many diseases using QTL mapping or GWAS analyses. A few of these studies pointed mainly toward adaptive mechanisms of resistance/susceptibility to infections; others pointed toward innate or intrinsic mechanisms. However, in the majority of studies, underlying mechanisms remain unknown. By comparing gene expression profiles between resistant and susceptible genetic backgrounds, transcriptomics studies have contributed to build a framework of gene pathways determining fish responsiveness to a number of pathogens. Adding functional assays to expression and genetic approaches has led to a better understanding of resistance mechanisms in some cases. The development of knock-out approaches will complement these analyses and help to validate putative candidate genes critical for resistance to infections. In this review, we highlight fish isogenic lines as a unique biological material to unravel the complexity of host response to different pathogens. In the future, combining multiple approaches will lead to a better understanding of the dynamics of interaction between the pathogen and the host immune response, and contribute to the identification of potential targets of selection for improved resistance.

Different transcriptional response between susceptible and resistant common carp (Cyprinus carpio) fish hints on the mechanism of CyHV-3 disease resistance

Background

Infectious disease outbreaks form major setbacks to aquaculture production and to further development of this important sector. Cyprinid herpes virus-3 (CyHV-3) is a dsDNA virus widely hampering production of common carp (Cyprinus carpio), one of the most farmed fish species worldwide. Genetically disease resistant strains are highly sought after as a sustainable solution to this problem. To study the genetic basis and cellular pathways underlying disease resistance, RNA-Seq was used to characterize transcriptional responses of susceptible and resistant fish at day 4 after CyHV-3 infection.

Results

In susceptible fish, over four times more differentially expressed genes were up-regulated between day 0 and 4 compared to resistant fish. Susceptible and resistant fish responded distinctively to infection as only 55 (9%) of the up-regulated genes were shared by these two fish types. Susceptible fish elicited a typical anti-viral response, involving interferon and interferon responsive genes, earlier than resistant fish did. Furthermore, chemokine profiles indicated that the two fish types elicited different cellular immunity responses. A comparative phylogenetic approach assisted in chemokine copies annotation pointing to different orthologous copies common to bony-fishes and even carp-specific paralogs that were differentially regulated and contributed to the different response of these two fish types. Susceptible fish up-regulated more ccl19 chemokines, which attract T-cells and macrophages, the anti-viral role of which is established, whereas resistant fish up-regulated more cxcl8/il8 chemokines, which attract neutrophils, the antiviral role of which is unfamiliar.

Conclusions

Taken together, by pointing out transcriptional differences between susceptible and resistant fish in response to CyHV-3 infection, this study unraveled possible genes and pathways that take part in disease resistance mechanisms in fish and thus, enhances our understanding of fish immunogenetics and supports the development of sustainable and safe aquaculture.

Fish antiviral tripartite motif (TRIM) proteins

Tripartite motif (TRIM) family or RBCC proteins comprises characteristic zinc-binding domains (a RING (R), a B-box type 1 (B1) and a B-box type 2 (B2)) and coiled-coil (CC) domain followed by a C-terminus variable domain. There are about 80 different TRIM proteins in human, but more than 200 in zebrafish with several large gene expansions (ftr >70 genes; btr >30 genes; trim35 > 30 genes). Repertoires of trim genes in fish are variable across fishes, but they have been remarkably diversified independently in a number of species. In mammals, TRIM proteins are involved in antiviral immunity through an astonishing diversity of mechanisms, from direct viral restriction to modulation of immune signaling and more recently autophagy. In fish, the antiviral role of TRIM proteins remains poorly understood. In zebrafish, fish specific TRIMs so called fintrims show a signature of positive selection in the C terminus SPRY domain, reminding features of mammalian antiviral trims such as TRIM5. Expression studies show that a number of trim genes, including many fintrims, can be induced during viral infections, and may play a role in antiviral defence. Some of them trigger antiviral activity in vitro against DNA and RNA viruses, such as FTR83 that also up-regulates the expression of type I IFN in zebrafish larvae. The tissue distribution of TRIM expression suggests that they may be involved in the regionalization of antiviral immunity, providing a particular protection to sensitive areas exposed to invading pathogens.

Resistance of pearlspot larvae, Etroplus suratensis, to redspotted grouper nervous necrosis virus by immersion challenge

Viral nervous necrosis (VNN) affects more than 120 species mostly belonging to the order Perciformes. However, none of the brackishwater species belonging to the family Cichlidae under the order Perciformes are reported to be susceptible. Hence, the present experiment was undertaken to study the susceptibility of the brackishwater cichlid, pearlspot, Etroplus suratensis to NNV. Thirty-day-old pearlspot larvae were infected with NNV by immersion. Mortality was recorded till 14 days post-infection, and the infected larvae were subjected to nested RT-PCR and histology. The virus was isolated from infected larvae using SSN-1 cells. To study the replication of the virus in vitro, primary cultured brain cells of E. suratensis and IEK cells were infected with NNV. No mortality was observed in any of the control or experimentally infected larvae. However, the experimentally infected larvae were positive for NNV by nested RT-PCR and the virus was isolated using SSN-1 cells. Further, the infected pearlspot brain cells and IEK cells showed cytopathic effect at second and third passage of the virus and they were positive for NNV by nested RT-PCR. Pearlspot is relatively resistant to VNN although the virus could replicate in the larvae and in cell culture.

Quantitative trait loci for resistance to Flavobacterium psychrophilum in rainbow trout: effect of the mode of infection and evidence of epistatic interactions

Background

Bacterial cold-water disease, which is caused by Flavobacterium psychrophilum, is one of the major diseases that affect rainbow trout (Oncorhynchus mykiss) and a primary concern for trout farming. Better knowledge of the genetic basis of resistance to F. psychrophilum would help to implement this trait in selection schemes and to investigate the immune mechanisms associated with resistance. Various studies have revealed that skin and mucus may contribute to response to infection. However, previous quantitative trait loci (QTL) studies were conducted by using injection as the route of infection. Immersion challenge, which is assumed to mimic natural infection by F. psychrophilum more closely, may reveal different defence mechanisms.

Results

Two isogenic lines of rainbow trout with contrasting susceptibilities to F. psychrophilum were crossed to produce doubled haploid F2 progeny. Fish were infected with F. psychrophilum either by intramuscular injection (115 individuals) or by immersion (195 individuals), and genotyped for 9654 markers using RAD-sequencing. Fifteen QTL associated with resistance traits were detected and only three QTL were common between the injection and immersion. Using a model that accounted for epistatic interactions between QTL, two main types of interactions were revealed. A “compensation-like” effect was detected between several pairs of QTL for the two modes of infection. An “enhancing-like” interaction effect was detected between four pairs of QTL. Integration of the QTL results with results of a previous transcriptomic analysis of response to F. psychrophilum infection resulted in a list of potential candidate immune genes that belong to four relevant functional categories (bacterial sensors, effectors of antibacterial immunity, inflammatory factors and interferon-stimulated genes).

Conclusions

These results provide new insights into the genetic determinism of rainbow trout resistance to F. psychrophilum and confirm that some QTL with large effects are involved in this trait. For the first time, the role of epistatic interactions between resistance-associated QTL was evidenced. We found that the infection protocol used had an effect on the modulation of defence mechanisms and also identified relevant immune functional candidate genes.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0431-9) contains supplementary material, which is available to authorized users.