Transcriptome profiling of the antiviral immune response in Atlantic cod macrophages

Water-in-oil adjuvant challenges in fish vaccination: An experimental inactivated adjuvanted vaccine against betanodavirus infection in Senegalese sole

The extensive growth of intensive fish farming has led to a massive spread of infectious diseases. Nervous necrosis virus (NNV) is the causative agent of the viral encephalo- and retinopathy disease which has become a major threat for fish farming all over the globe. The devastating mortality rates recorded in disease outbreaks, especially when infected specimens are at early stages of development, have a high economic impact on the sector. Currently, vaccines are the most cost-effective preventing tool in the fight against viruses. Inactivated vaccines have the advantage of simplicity in their development at the same time as present the antigen in a similar manner than the natural infection in the host. Nevertheless, they usually trigger weaker immune responses needing adjuvants to boost their effectiveness. In this work, we have intraperitoneally vaccinated Senegalese sole juveniles (Solea senegalensis) with a previously designed inactivated vaccine against NNV based on binary ethylenimine (BEI), mixed or not with an oil-adjuvant. Our results demonstrated the potential activation of different immune pathways when the vaccine was administered alone compared to the oil-adjuvanted vaccine, both resulting in an equivalent partial improvement in survival following a NNV challenge. However, whilst the vaccine alone led to a significant increase in specific antibodies, in the adjuvanted version those antibodies were kept basal although with a slight improvement in their neutralization capacity. At transcriptional level, neither vaccine (adjuvanted or not) triggered the immune system activation during the vaccination period. However, after NNV infection, the BEI-inactivated vaccines alone and oil-adjuvanted both elicited the stimulation of antiviral responsive genes (rtp3, herc4), antigen presentation molecules (mhcii) and T-cell markers (cd8a) in the head-kidney. Additionally, the oil-adjuvanted vaccine appears to stimulate mediator cytokines (il6) and B-cell markers (ight and ighm). Surprisingly, when the adjuvant was administered alone, fish showed the highest survival rates concomitantly with a lack of NNV-IgM production, pointing to the possible induction of different immune pathways than the B-cell responses via antibodies by the adjuvant. Since this combined vaccine did not succeed in the full extension of protection against the pathogen, further studies should be performed focusing on unravelling the molecular mechanisms through which adjuvants trigger the immune response, both independently and when added to a vaccine antigen.

Lost and Found: The Family of NF-κB Inhibitors Is Larger than Assumed in Salmonid Fish

NF-κB signalling is largely controlled by the family of 'inhibitors of NF-κB' (IκB). The relevant databases indicate that the genome of rainbow trout contains multiple gene copies coding for iκbα (nfkbia), iκbε (nfkbie), iκbδ (nkfbid), iκbζ (nfkbiz), and bcl3, but it lacks iκbβ (nfkbib) and iκbη (ankrd42). Strikingly, three nfkbia paralogs are apparently present in salmonid fish, two of which share a high sequence identity, while the third putative nfkbia gene is significantly less like its two paralogs. This particular nfkbia gene product, iκbα, clusters with the human IκBβ in a phylogenetic analysis, while the other two iκbα proteins from trout associate with their human IκBα counterpart. The transcript concentrations were significantly higher for the structurally more closely related nfkbia paralogs than for the structurally less similar paralog, suggesting that iκbβ probably has not been lost from the salmonid genomes but has been incorrectly designated as iκbα. In the present study, two gene variants coding for iκbα (nfkbia) and iκbε (nfkbie) were prominently expressed in the immune tissues and, particularly, in a cell fraction enriched with granulocytes, monocytes/macrophages, and dendritic cells from the head kidney of rainbow trout. Stimulation of salmonid CHSE-214 cells with zymosan significantly upregulated the iκbα-encoding gene while elevating the copy numbers of the inflammatory markers interleukin-1-beta and interleukin-8. Overexpression of iκbα and iκbε in CHSE-214 cells dose-dependently quenched both the basal and stimulated activity of an NF-κB promoter suggesting their involvement in immune-regulatory processes. This study provides the first functional data on iκbε-versus the well-researched iκbα factor-in a non-mammalian model species.

Ancestral reconstruction reveals catalytic inactivation of activation-induced cytidine deaminase concomitant with cold water adaption in the Gadiformes bony fish

BACKGROUND

Antibody affinity maturation in vertebrates requires the enzyme activation-induced cytidine deaminase (AID) which initiates secondary antibody diversification by mutating the immunoglobulin loci. AID-driven antibody diversification is conserved across jawed vertebrates since bony and cartilaginous fish. Two exceptions have recently been reported, the Pipefish and Anglerfish, in which the AID-encoding aicda gene has been lost. Both cases are associated with unusual reproductive behavior, including male pregnancy and sexual parasitism. Several cold water fish in the Atlantic cod (Gadinae) family carry an aicda gene that encodes for a full-length enzyme but lack affinity-matured antibodies and rely on antibodies of broad antigenic specificity. Hence, we examined the functionality of their AID.

RESULTS

By combining genomics, transcriptomics, immune responsiveness, and functional enzymology of AID from 36 extant species, we demonstrate that AID of that Atlantic cod and related fish have extremely lethargic or no catalytic activity. Through ancestral reconstruction and functional enzymology of 71 AID enzymes, we show that this enzymatic inactivation likely took place relatively recently at the emergence of the true cod family (Gadidae) from their ancestral Gadiformes order. We show that this AID inactivation is not only concordant with the previously shown loss of key adaptive immune genes and expansion of innate and cell-based immune genes in the Gadiformes but is further reflected in the genomes of these fish in the form of loss of AID-favored sequence motifs in their immunoglobulin variable region genes.

CONCLUSIONS

Recent demonstrations of the loss of the aicda gene in two fish species challenge the paradigm that AID-driven secondary antibody diversification is absolutely conserved in jawed vertebrates. These species have unusual reproductive behaviors forming an evolutionary pressure for a certain loss of immunity to avoid tissue rejection. We report here an instance of catalytic inactivation and functional loss of AID rather than gene loss in a conventionally reproducing vertebrate. Our data suggest that an expanded innate immunity, in addition to lower pathogenic pressures in a cold environment relieved the pressure to maintain robust secondary antibody diversification. We suggest that in this unique scenario, the AID-mediated collateral genome-wide damage would form an evolutionary pressure to lose AID function.

Nutritional immunomodulation of Atlantic salmon response to Renibacterium salmoninarum bacterin

We investigated the immunomodulatory effect of varying levels of dietary ω6/ω3 fatty acids (FA) on Atlantic salmon (Salmo salar) antibacterial response. Two groups were fed either high-18:3ω3 or high-18:2ω6 FA diets for 8 weeks, and a third group was fed for 4 weeks on the high-18:2ω6 diet followed by 4 weeks on the high-18:3ω3 diet and termed “switched-diet”. Following the second 4 weeks of feeding (i.e., at 8 weeks), head kidney tissues from all groups were sampled for FA analysis. Fish were then intraperitoneally injected with either a formalin-killed Renibacterium salmoninarum bacterin (5 × 10⁷ cells mL⁻¹) or phosphate-buffered saline (PBS control), and head kidney tissues for gene expression analysis were sampled at 24 h post-injection. FA analysis showed that the head kidney profile reflected the dietary FA, especially for C₁₈ FAs. The qPCR analyses of twenty-three genes showed that both the high-ω6 and high-ω3 groups had significant bacterin-dependent induction of some transcripts involved in lipid metabolism (ch25ha and lipe), pathogen recognition (clec12b and tlr5), and immune effectors (znrf1 and cish). In contrast, these transcripts did not significantly respond to the bacterin in the “switched-diet” group. Concurrently, biomarkers encoding proteins with putative roles in biotic inflammatory response (tnfrsf6b) and dendritic cell maturation (ccl13) were upregulated, and a chemokine receptor (cxcr1) was downregulated with the bacterin injection regardless of the experimental diets. On the other hand, an inflammatory regulator biomarker, bcl3, was only significantly upregulated in the high-ω3 fed group, and a C-type lectin family member (clec3a) was only significantly downregulated in the switched-diet group with the bacterin injection (compared with diet-matched PBS-injected controls). Transcript fold-change (FC: bacterin/PBS) showed that tlr5 was significantly over 2-fold higher in the high-18:2ω6 diet group compared with other diet groups. FC and FA associations highlighted the role of DGLA (20:3ω6; anti-inflammatory) and/or EPA (20:5ω3; anti-inflammatory) vs. ARA (20:4ω6; pro-inflammatory) as representative of the anti-inflammatory/pro-inflammatory balance between eicosanoid precursors. Also, the correlations revealed associations of FA proportions (% total FA) and FA ratios with several eicosanoid and immune receptor biomarkers (e.g., DGLA/ARA significant positive correlation with pgds, 5loxa, 5loxb, tlr5, and cxcr1). In summary, dietary FA profiles and/or regimens modulated the expression of some immune-relevant genes in Atlantic salmon injected with R. salmoninarum bacterin. The modulation of Atlantic salmon responses to bacterial pathogens and their associated antigens using high-ω6/high-ω3 diets warrants further investigation.

Integrated analysis of immune parameters, miRNA-mRNA interaction, and immune genes expression in the liver of rainbow trout following infectious hematopoietic necrosis virus infection

Rainbow trout (Oncorhynchus mykiss) is an important economical cold-water fish worldwide. However, infection with infectious hematopoietic necrosis virus (IHNV) has severely restricted the development of aquaculture and caused huge economic losses. Currently, little is known about the immune defense mechanisms of rainbow trout against IHNV. In this study, we detected the changes of immune parameters over different post-infection periods (6-, 12-, 24-, 48-, 72-, 96-, 120-, and 144 hours post-infection (hpi)), mRNA and miRNA expression profiles under 48 hpi (T48L) compared to control (C48L), and key immune-related genes expression patterns in rainbow trout liver following IHNV challenge through biochemical methods, RNA-seq, and qRT-PCR, and the function of miR-330-y was verified by overexpression and silencing in vitro and in vivo. The results revealed that alkaline phosphatase (AKP), alanine aminotransferase (ALT), catalase (CAT), and total superoxide dismutase (T-SOD) activities, and lysozyme (LZM) content showed significant peaks at 48 hpi, whereas malondialdehyde (MDA) content and aspartate aminotransferase (AST) activity decreased continuously during infection, and acid phosphatase (ACP) activity varied slightly. From RNA-seq, a total of 6844 genes and 86 miRNAs were differentially expressed, and numerous immune-related differentially expressed genes (DEGs) involved in RIG-I-like receptor signaling pathway, Toll-like receptor signaling pathway, NOD-like receptor signaling pathway, cytokine-cytokine receptor interaction, and antigen processing and presentation were significantly upregulated in T48Lm group, including IFIH1, DHX58, MAVS, TRAF3, IRF3, IRF7, MX1, TLR3, TLR8, MYD88, NOD1, NOD2, IL-8, CXCR1, CD209, CD83, and TAP1. Integrated analysis identified seven miRNAs (miR-425-x, miR-185-x, miR-338-x, miR-330-y, miR-361-x, miR-505-y, and miR-191-x) that target at least three key immune-related DEGs. Expression analysis showed that IFIH1, DHX58, IRF3, IRF7, MX1, TLR3, TLR8, and MYD88 showed a marked increase after 24 hpi during infection. Further research confirmed TAP1 as one of the targets of miR-330-y, overexpression of miR-330-y with mimics or agomir significantly reduced the expression levels of TAP1, IRF3, and IFN, and the opposite effects were obtained by inhibitor. These results facilitate in-depth understanding of the immune mechanisms in rainbow trout against IHNV.

The Skeleton and Biomineralization Mechanism as Part of the Innate Immune System of Stony Corals

Stony corals are among the most important calcifiers in the marine ecosystem as they form the coral reefs. Coral reefs have huge ecological importance as they constitute the most diverse marine ecosystem, providing a home to roughly a quarter of all marine species. In recent years, many studies have shed light on the mechanisms underlying the biomineralization processes in corals, as characterizing the calicoblast cell layer and genes involved in the formation of the calcium carbonate skeleton. In addition, considerable advancements have been made in the research field of coral immunity as characterizing genes involved in the immune response to pathogens and stressors, and the revealing of specialized immune cells, including their gene expression profile and phagocytosis capabilities. Yet, these two fields of corals research have never been integrated. Here, we discuss how the coral skeleton plays a role as the first line of defense. We integrate the knowledge from both fields and highlight genes and proteins that are related to biomineralization and might be involved in the innate immune response and help the coral deal with pathogens that penetrate its skeleton. In many organisms, the immune system has been tied to calcification. In humans, immune factors enhance ectopic calcification which causes severe diseases. Further investigation of coral immune genes which are involved in skeleton defense as well as in biomineralization might shed light on our understanding of the correlation and the interaction of both processes as well as reveal novel comprehension of how immune factors enhance calcification.

Chitosan nanoparticle formulation attenuates poly (I:C) induced innate immune responses against inactivated virus vaccine in Atlantic salmon (Salmo salar)

Many vaccine formulations, in particular vaccines based on inactivated virus, needs adjuvants to boost immunogenicity. In aquaculture, mineral and plant oil are used as adjuvant in commercial vaccines, and the advent of oil-adjuvanted vaccines was crucial to aquaculture development. Nevertheless, some of these approved vaccines display suboptimal performance in the field compared to experimental conditions. Therefore, there is a need to improve adjuvants and delivery methods for fish vaccines against viruses. We used RNA sequencing of Atlantic salmon head kidney to analyse the difference in gene expression 24 h after injection of different experimental vaccine formulations. We compared five different formulations in addition to a PBS control: inactivated virus alone (group V), soluble poly (I:C) (group P), nanoparticles containing poly (I:C) (group N), soluble poly (I:C) + inactivated virus (group PV) and finally nanoparticles containing poly (I:C) + inactivated virus (group NV). Our results showed poly (I:C)'s ability as adjuvant and its capacity influence innate immune genes expression in Atlantic salmon. Soluble poly (I:C) upregulated multiple immune related genes and was more effective compared to poly (I:C) formulated into chitosan nanoparticles (more than 10 fold increase in differentially expressed genes, DEGs). However, inclusion of inactivated ISA virus in the nanoparticle vaccine, increased the number of DEGs fivefold suggesting a synergistic effect of adjuvant and antigen. Our results indicate that the way poly (I:C) is formulated and the presence of antigen is important for the magnitude of the innate immune response in Atlantic salmon.

BEI Inactivated Vaccine Induces Innate and Adaptive Responses and Elicits Partial Protection upon Reassortant Betanodavirus Infection in Senegalese Sole

Nervous necrosis virus (NNV), the causative agent of viral encephalopathy and retinopathy (VER), is one of the most threatening viruses affecting marine and freshwater fish species worldwide. Senegalese sole is a promising fish species in Mediterranean aquaculture but also highly susceptible to NNV and VER outbreaks, that puts its farming at risk. The development of vaccines for aquaculture is one of best tools to prevent viral spread and sudden outbreaks, and virus inactivation is the simplest and most cost-effective method available. In this work, we have designed two inactivated vaccines based on the use of formalin or binary ethylenimine (BEI) to inactivate a reassortant NNV strain. After vaccination, the BEI-inactivated vaccine triggered the production of specific IgM-NNV antibodies and stimulated innate and adaptive immune responses at transcriptional level (rtp3, mx, mhcii and tcrb coding genes). Moreover, it partially improved survival after an NNV in vivo challenge, reducing the mid-term viral load and avoiding the down-regulation of immune response post-challenge. On the other hand, the formalin-inactivated vaccine improved the survival of fish upon infection without inducing the production of IgM-NNV antibodies and only stimulating the expression of herc4 and mhcii genes (in head-kidney and brain, respectively) during the vaccination period; this suggests that other immune-related pathways may be involved in the partial protection provoked. Although these vaccines against NNV showed encouraging results, further studies are needed to improve sole protection and to fully understand the underlying immune mechanism.

Immunogene expression analysis in betanodavirus infected-Senegalese sole using an OpenArray® platform

The transcriptomic response of Senegalese sole (Solea senegalensis) triggered by two betanodaviruses with different virulence to that fish species has been assessed using an OpenArray® platform based on TaqMan™ quantitative PCR. The transcription of 112 genes per sample has been evaluated at two sampling times in two organs (head kidney and eye/brain-pooled samples). Those genes were involved in several roles or pathways, such as viral recognition, regulation of type I (IFN-1)-dependent immune responses, JAK-STAT cascade, interferon stimulated genes, protein ubiquitination, virus responsive genes, complement system, inflammatory response, other immune system effectors, regulation of T-cell proliferation, and proteolysis and apoptosis. The highly virulent isolate, wSs160.3, a wild type reassortant containing a RGNNV-type RNA1 and a SJNNV-type RNA2 segments, induced the expression of a higher number of genes in both tested organs than the moderately virulent strain, a recombinant harbouring mutations in the protruding domain of the capsid protein. The number of differentially expressed genes was higher 2 days after the infection with the wild type isolate than at 3 days post-inoculation. The wild type isolate also elicited an exacerbated interferon 1 response, which, instead of protecting sole against the infection, increases the disease severity by the induction of apoptosis and inflammation-derived immunopathology, although inflammation seems to be modulated by the complement system. Furthermore, results derived from this study suggest a potential important role for some genes with high expression after infection with the highly virulent virus, such as rtp3, sacs and isg15. On the other hand, the infection with the mutant does not induce immune response, probably due to an altered recognition by the host, which is supported by a different viral recognition pathway, involving myd88 and tbkbp1.

Innovation in Nucleotide-Binding Oligomerization-Like Receptor and Toll-Like Receptor Sensing Drives the Major Histocompatibility Complex-II Free Atlantic Cod Immune System

The absence of MHC class II antigen presentation and multiple pathogen recognition receptors in the Atlantic cod has not impaired its immune response however how underlying mechanisms have adapted remains largely unknown. In this study, ex vivo cod macrophages were challenged with various bacterial and viral microbe-associated molecular patterns (MAMP) to identify major response pathways. Cytosolic MAMP-PRR pathways based upon the NOD-like receptors (NLRs) and RIG-I-like receptors (RLRs) were identified as the critical response pathways. Our analyses suggest that internalization of exogenous ligands through scavenger receptors drives both pathways activating transcription factors like NF-kB (Nuclear factor-kappa B) and interferon regulatory factors (IRFs). Further, ligand-dependent differential expression of a unique TLR25 isoform and multiple NLR paralogues suggests (sub)neofunctionalization toward specific immune defensive strategies. Our results further demonstrate that the unique immune system of the Atlantic cod provides an unprecedented opportunity to explore the evolutionary history of PRR-based signaling in vertebrate immunity.

Identification of an IRF10 gene in common carp (Cyprinus carpio L.) and analysis of its function in the antiviral and antibacterial immune response

Background

Interferon (IFN) regulatory factors (IRFs), as transcriptional regulatory factors, play important roles in regulating the expression of type I IFN and IFN- stimulated genes (ISGs) in innate immune responses. In addition, they participate in cell growth and development and regulate oncogenesis.

Results

In the present study, the cDNA sequence of IRF10 in common carp (Cyprinus carpio L.) was characterized (abbreviation, CcIRF10). The predicted protein sequence of CcIRF10 shared 52.7–89.2% identity with other teleost IRF10s and contained a DNA-binding domain (DBD), a nuclear localization signal (NLS) and an IRF-associated domain (IAD). Phylogenetic analysis showed that CcIRF10 had the closest relationship with IRF10 of Ctenopharyngodon idella. CcIRF10 transcripts were detectable in all examined tissues, with the highest expression in the gonad and the lowest expression in the head kidney. CcIRF10 expression was upregulated in the spleen, head kidney, foregut and hindgut upon polyinosinic:polycytidylic acid (poly I:C) and Aeromonas hydrophila stimulation and induced by poly I:C, lipopolysaccharide (LPS) and peptidoglycan (PGN) in peripheral blood leucocytes (PBLs) and head kidney leukocytes (HKLs) of C. carpio. In addition, overexpression of CcIRF10 was able to decrease the expression of the IFN and IFN-stimulated genes PKR and ISG15.

Conclusions

These results indicate that CcIRF10 participates in antiviral and antibacterial immunity and negatively regulates the IFN response, which provides new insights into the IFN system of C. carpio.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-020-02674-z.

Integration of RNA-seq and RNAi provides a novel insight into the immune responses of Epinephelus coioides to the impB gene of Pseudomonas plecoglossicida

Pseudomonas plecoglossicida is a Gram-negative bacterium that causes visceral white spot disease in Epinephelus coioides and leads to severe aquatic economic losses. The RNA-seq results of a previous study showed that the expression of the impB gene in P. plecoglossicida was significantly upregulated during infection. Four shRNAs were designed and synthesized to silence the impB gene in P. plecoglossicida, and the maximum silencing efficiency was 95.2%. Intraperitoneal injection of the impB-RNAi strain of P. plecoglossicida did not cause E. coioides death, and the spleens of infected fish did not show significant clinical symptoms. Although the injection of the mutant strain increased the antibody titer in E. coioides serum, it could not effectively protect E. coioides against wild strain infection. Compared with E. coioides infected with the wild type strain, the RNA-seq results for E. coioides infected with the impB-RNAi strain differed greatly. The KEGG enrichment analysis showed that key genes of the chemokine signalling pathway of E. coioides were downregulated by the silencing of impB in P. plecoglossicida. Infection with the impB-RNAi strain of P. plecoglossicida through injection did not produce good immune protection against E. coioides. The present study provides a novel insight into the immune responses of E. coioides to the impB gene of P. plecoglossicida.

Profiling the transcriptome response of Atlantic salmon head kidney to formalin-killed Renibacterium salmoninarum

Renibacterium salmoninarum is a Gram-positive, intracellular bacterial pathogen that causes Bacterial Kidney Disease (BKD) in Atlantic salmon (Salmo salar). The host transcriptomic response to this immune-suppressive pathogen remains poorly understood. To identify R. salmoninarum-responsive genes, Atlantic salmon were intraperitoneally injected with a low (5 × 10⁵ cells/kg, Low-Rs) or high (5 × 10⁷ cells/kg; High-Rs) dose of formalin-killed R. salmoninarum bacterin or phosphate-buffered saline (PBS control); head kidney samples were collected before and 24 h after injection. Using 44K microarray analysis, we identified 107 and 345 differentially expressed probes in response to R. salmoninarum bacterin (i.e. High-Rs vs. PBS control) by Significance Analysis of Microarrays (SAM) and Rank Products (RP), respectively. Twenty-two microarray-identified genes were subjected to qPCR assays, and 17 genes were confirmed as being significantly responsive to the bacterin. There was an up-regulation in expression of genes playing putative roles as immune receptors and antimicrobial effectors. Genes with putative roles as pathogen recognition (e.g. clec12b and tlr5) or immunoregulatory (e.g. tnfrsf6b and tnfrsf11b) receptors were up-regulated in response to R.salmoninarum bacterin. Also, chemokines and a chemokine receptor showed opposite regulation [up-regulation of effectors (i.e. ccl13 and ccl) and down-regulation of cxcr1] in response to the bacterin. The present study identified and validated novel biomarker genes (e.g. ctsl1, lipe, cldn4, ccny) that can be used to assess Atlantic salmon response to R. salmoninarum, and will be valuable in the development of tools to combat BKD.

Impact of rearing temperature on the innate antiviral immune response of growth hormone transgenic female triploid Atlantic salmon (Salmo salar)

AquAdvantage Salmon (growth hormone transgenic female triploid Atlantic salmon) are a faster-growing alternative to conventional farmed diploid Atlantic salmon. To investigate optimal rearing conditions for their commercial production, a laboratory study was conducted in a freshwater recirculating aquaculture system (RAS) to examine the effect of rearing temperature (10.5 °C, 13.5 °C, 16.5 °C) on their antiviral immune and stress responses. When each temperature treatment group reached an average weight of 800 g, a subset of fish were intraperitoneally injected with either polyriboinosinic polyribocytidylic acid (pIC, a viral mimic) or an equal volume of sterile phosphate-buffered saline (PBS). Blood and head kidney samples were collected before injection and 6, 24 and 48 h post-injection (hpi). Transcript abundance of 7 antiviral biomarker genes (tlr3, lgp2, stat1b, isg15a, rsad2, mxb, ifng) was measured by real-time quantitative polymerase chain reaction (qPCR) on head kidney RNA samples. Plasma cortisol levels from blood samples collected pre-injection and from pIC and PBS groups at 24 hpi were quantified by ELISA. While rearing temperature and treatment did not significantly affect circulating cortisol, all genes tested were significantly upregulated by pIC at all three temperatures (except for tlr3, which was only upregulated in the 10.5 °C treatment). Target gene activation was generally observed at 24 hpi, with most transcript levels decreasing by 48 hpi in pIC-injected fish. Although a high amount of biological variability in response to pIC was evident across all treatments, rearing temperature significantly influenced transcript abundance and/or fold-changes comparing time- and temperature-matched pIC- and PBS-injected fish for several genes (tlr3, lgp2, stat1b, isg15a, rsad2 and ifng) at 24 hpi. As an example, significantly higher fold-changes of rsad2, isg15a and ifng were found in fish reared at 10.5 °C when compared to 16.5 °C. Multivariate analysis confirmed that rearing temperature modulated antiviral immune response. The present experiment provides novel insight into the relationship between rearing temperature and innate antiviral immune response in AquAdvantage Salmon.

Dietary Immunostimulant CpG Modulates MicroRNA Biomarkers Associated with Immune Responses in Atlantic Salmon (Salmo salar)

MicroRNAs (miRNAs) are key regulators in fish immune responses. However, no study has previously characterized the impact of polyriboinosinic polyribocytidylic acid (pIC) and formalin-killed typical Aeromonas salmonicida (ASAL) on miRNA expression in Atlantic salmon fed a commercial diet with and without immunostimulant CpG. To this end, first, we performed small RNA deep sequencing and qPCR analyses to identify and confirm pIC- and/or ASAL-responsive miRNAs in the head kidney of salmon fed a control diet. DESeq2 analyses identified 12 and 18 miRNAs differentially expressed in pIC and ASAL groups, respectively, compared to the controls. Fifteen of these miRNAs were studied by qPCR; nine remained significant by qPCR. Five miRNAs (miR-27d-1-2-5p, miR-29b-2-5p, miR-146a-5p, miR-146a-1-2-3p, miR-221-5p) were shown by qPCR to be significantly induced by both pIC and ASAL. Second, the effect of CpG-containing functional feed on miRNA expression was investigated by qPCR. In pre-injection samples, 6 of 15 miRNAs (e.g., miR-181a-5-3p, miR-462a-3p, miR-722-3p) had significantly lower expression in fish fed CpG diet than control diet. In contrast, several miRNAs (e.g., miR-146a-1-2-3p, miR-192a-5p, miR-194a-5p) in the PBS- and ASAL-injected groups had significantly higher expression in CpG-fed fish. Multivariate statistical analyses confirmed that the CpG diet had a greater impact on miRNA expression in ASAL-injected compared with pIC-injected fish. This study identified immune-relevant miRNA biomarkers that will be valuable in the development of diets to combat infectious diseases of salmon.

Molecular Regulation of Host Defense Responses Mediated by Biological Anti-TMV Agent Ningnanmycin

Ningnanmycin (NNM) belongs to microbial pesticides that display comprehensive antiviral activity against plant viruses. NNM treatment has been shown to efficiently delay or suppress the disease symptoms caused by tobacco mosaic virus (TMV) infection in local-inoculated or systemic-uninoculated tobacco leaves, respectively. However, the underlying molecular mechanism of NNM-mediated antiviral activity remains to be further elucidated. In this study, 414 differentially expressed genes (DEGs), including 383 which were up-regulated and 31 down-regulated, caused by NNM treatment in TMV-infected BY-2 protoplasts, were discovered by RNA-seq. In addition, KEGG analysis indicated significant enrichment of DEGs in the plant-pathogen interaction and MAPK signaling pathway. The up-regulated expression of crucial DEGs, including defense-responsive genes, such as the receptor-like kinase FLS2, RLK1, and the mitogen-activated protein kinase kinase kinase MAPKKK, calcium signaling genes, such as the calcium-binding protein CML19, as well as phytohormone responsive genes, such as the WRKY transcription factors WRKY40 and WRKY70, were confirmed by RT-qPCR. These findings provided valuable insights into the antiviral mechanisms of NNM, which indicated that the agent induces tobacco systemic resistance against TMV via activating multiple plant defense signaling pathways.

Full characterization and transcript expression profiling of the interferon regulatory factor (IRF) gene family in Atlantic cod (Gadus morhua)

Atlantic cod (Gadus morhua) represents a unique immune system among teleost fish, making it a species of interest for immunological studies, and especially for investigating the evolutionary history of immune gene families. The interferon regulatory factor (IRF) gene family encodes transcription factors which function in the interferon pathway, but also in areas including leukocyte differentiation, cell growth, autoimmunity, and development. We previously characterized several IRF family members in Atlantic cod (Irf4a, Irf4b, Irf7, Irf8, and two Irf10 splice variants) at the cDNA and putative amino acid levels, and in the current study we took advantage of the new and improved Atlantic cod genome assembly in combination with rapid amplification of cDNA ends (RACE) to characterize the remaining family members (i.e. Irf3, Irf5, Irf6, Irf9, and two Irf2 splice variants). Real-time quantitative PCR (QPCR) was used to investigate constitutive expression of all IRF transcripts during embryonic development, suggesting several putative maternal transcripts, and potential stage-specific roles. QPCR studies also showed 11 of 13 transcripts were responsive to stimulation with poly(I:C), while 6 of 13 transcripts were responsive to lipopolysaccharide (LPS) in Atlantic cod head kidney macrophages, indicating roles for cod IRF family members in both antiviral and antibacterial responses. This study is the first to investigate expression of the complete IRF family in Atlantic cod, and suggests potential novel roles for several of these transcription factors within immunity as well as in early development of this species.

Aeromonas salmonicida subsp. salmonicida Early Infection and Immune Response of Atlantic Cod (Gadus morhua L.) Primary Macrophages

In contrast to other teleosts, Atlantic cod (Gadus morhua) has an expanded repertoire of MHC-I and TLR components, but lacks the MHC-II, the invariant chain/CD74, and CD4⁺ T cell response, essential for production of antibodies and prevention of bacterial infectious diseases. The mechanisms by which G. morhua fight bacterial infections are not well understood. Aeromonas salmonicida subsp. salmonicida is a recurrent pathogen in cultured and wild fish, and has been reported in Atlantic cod. Macrophages are some of the first responders to bacterial infection and the link between innate and adaptive immune response. Here, we evaluated the viability, reactive oxygen species (ROS) production, cell morphology, and gene expression of cod primary macrophages in response to A. salmonicida infection. We found that A. salmonicida infects cod primary macrophages without killing the cod cells. Likewise, infected Atlantic cod macrophages up-regulated key genes involved in the inflammatory response (e.g., IL-1β and IL-8) and bacterial recognition (e.g., BPI/LBP). Nevertheless, our results showed a down-regulation of genes related to antimicrobial peptide and ROS production, suggesting that A. salmonicida utilizes its virulence mechanisms to control and prevent macrophage anti-bacterial activity. Our results also indicate that Atlantic cod has a basal ROS production in non-infected cells, and this was not increased after contact with A. salmonicida. Transmission electron microscopy results showed that A. salmonicida was able to infect the macrophages in a high number, and release outer membrane vesicles (OMV) during intracellular infection. These results suggest that Atlantic cod macrophage innate immunity is able to detect A. salmonicida and trigger an anti-inflammatory response, however A. salmonicida controls the cell immune response to prevent bacterial clearance, during early infection.

Whole transcriptome analysis of the Atlantic cod vaccine response reveals subtle changes in adaptive immunity

Atlantic cod has lost the Major Histocompatibility complex class II pathway - central to pathogen presentation, humoral response and immunity. Here, we investigate the immunological response of Atlantic cod subsequent to dip vaccination with Vibrioanguillarum bacterin using transcriptome sequencing. The experiment was conducted on siblings from an Atlantic cod family found to be highly susceptible towards vibriosis where vaccination has demonstrated improved pathogen resistance. Gene expression analyses at 2, 4, 21 and 42 days post vaccination revealed GO-term enrichment for muscle, neuron and metabolism-related pathways. In-depth characterization of immune-related GO terms demonstrated down-regulation of MHCI antigen presentation, C-type lectin receptor signaling and granulocyte activation over time. Phagocytosis, interferon-gamma signaling and negative regulation of innate immunity were increasingly up-regulated over time. Individual differentially expressed immune genes implies weak initiation of acute phase proteins with little or no inflammation. Furthermore, gene expression indicates presence of T-cells, NK-like cells, B-cells and monocytes/macrophages. Three MHCI transcripts were up-regulated with B2M and SEC61. Overall, we find no clear immune-related transcriptomic response which could be attributed to Atlantic cod's alternative immune system. However, we cannot rule out that this response is related to vaccination protocol/sampling strategy. Earlier functional studies demonstrate significant memory in Atlantic cod post dip vaccination and combined with our results indicate the presence of other adaptive immunity mechanisms. In particular, we suggest that further investigations should look into CD8+ memory T-cells, γδ T-cells, T-cell independent memory or memory induced through NK-like/other lymphoid cells locally in the mucosal lining for this particular vaccination strategy.

Characterization and Transcript Expression Analyses of Atlantic Cod Viperin

Viperin is a key antiviral effector in immune responses of vertebrates including the Atlantic cod (Gadus morhua). Using cloning, sequencing and gene expression analyses, we characterized the Atlantic cod viperin at the nucleotide and hypothetical amino acid levels, and its regulating factors were investigated. Atlantic cod viperin cDNA is 1,342 bp long, and its predicted protein contains 347 amino acids. Using in silico analyses, we showed that Atlantic cod viperin is composed of 5 exons, as in other vertebrate orthologs. In addition, the radical SAM domain and C-terminal sequences of the predicted Viperin protein are highly conserved among various species. As expected, Atlantic cod Viperin was most closely related to other teleost orthologs. Using computational modeling, we show that the Atlantic cod Viperin forms similar overall protein architecture compared to mammalian Viperins. qPCR revealed that viperin is a weakly expressed transcript during embryonic development of Atlantic cod. In adults, the highest constitutive expression of viperin transcript was found in blood compared with 18 other tissues. Using isolated macrophages and synthetic dsRNA (pIC) stimulation, we tested various immune inhibitors to determine the possible regulating pathways of Atlantic cod viperin. Atlantic cod viperin showed a comparable pIC induction to other well-known antiviral genes (e.g., interferon gamma and interferon-stimulated gene 15-1) in response to various immune inhibitors. The pIC induction of Atlantic cod viperin was significantly inhibited with 2-Aminopurine, Chloroquine, SB202190, and Ruxolitinib. Therefore, endosomal-TLR-mediated pIC recognition and signal transducers (i.e., PKR and p38 MAPK) downstream of the TLR-dependent pathway may activate the gene expression response of Atlantic cod viperin. Also, these results suggest that antiviral responses of Atlantic cod viperin may be transcriptionally regulated through the interferon-activated pathway.

Changes in the liver transcriptome of farmed Atlantic salmon (Salmo salar) fed experimental diets based on terrestrial alternatives to fish meal and fish oil

Background

Dependence on marine natural resources threatens the sustainability of Atlantic salmon aquaculture. In the present study, Atlantic salmon fed for 14 weeks with an experimental diet based on animal by-products and vegetable oil (ABP) exhibited reduced growth performance compared with others fed a fish meal/fish oil based experimental diet (MAR) and a plant protein/vegetable oil-based experimental diet (VEG). To characterize the molecular changes underlying the differences in growth performance, we conducted a 44 K microarray study of the liver transcriptome of the three dietary groups.

Results

The microarray experiment identified 122 differentially expressed features (Rank Products, PFP < 10%). Based on their associated Gene Ontology terms, 46 probes were classified as metabolic and growth-relevant genes, 25 as immune-related, and 12 as related to oxidation-reduction processes. The microarray results were validated by qPCR analysis of 29 microarray-identified transcripts. Diets significantly modulated the transcription of genes involved in carbohydrate metabolism (gck and pfkfb4), cell growth and proliferation (sgk2 and htra1), apoptosis (gadd45b), lipid metabolism (fabp3, idi1, sqs), and immunity (igd, mx, ifit5, and mhcI). Hierarchical clustering and linear correlation analyses were performed to find gene expression patterns among the qPCR-analyzed transcripts, and connections between them and muscle and liver lipid composition. Overall, our results indicate that changes in the liver transcriptome and tissue lipid composition were driven by cholesterol synthesis up-regulation by ABP and VEG diets, and the lower carbohydrate intake in the ABP group. Two of the microarray-identified genes (sgk2 and htra1) might be key to explaining glucose metabolism regulation and the dietary-modulation of the immune system in fish. To evaluate the potential of these genes as predictive biomarkers, we subjected the qPCR data to a stepwise discriminant analysis. Three sets of no more than four genes were found to be able to predict, with high accuracy (67–94%), salmon growth and fatty acid composition.

Conclusions

This study provides new findings on the impact of terrestrial animal and plant products on the nutrition and health of farmed Atlantic salmon, and a new method based on gene biomarkers for potentially predicting desired phenotypes, which could help formulate superior feeds for the Atlantic salmon aquaculture industry.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5188-6) contains supplementary material, which is available to authorized users.

Applications of transcriptomics and proteomics in understanding fish immunity

With the development of intensive aquaculture, economic losses increasingly result from fish mortality due to pathogen infection. In recent years, a growing number of researchers have used transcriptomic and proteomic analyses to study fish immune responses to exogenous pathogen infection. Integrating transcriptomic and proteomic analyses provides a better understanding of the fish immune system including gene expression, regulation, and the intricate biological processes underlying immune responses against infection. This review focuses on the recent advances in the fields of transcriptomics and proteomics, which have contributed to our understanding of fish immunity to exogenous pathogens.

Discovery of microRNAs associated with the antiviral immune response of Atlantic cod macrophages

MicroRNAs (miRNAs) are known to play important immunoregulatory roles in teleosts, although miRNAs involved in the antiviral immune response of Atlantic cod (Gadus morhua) were previously uncharacterised. Using deep sequencing and qPCR, the present study was conducted to identify miRNAs responsive to the viral mimic, polyriboinosinic polyribocytidylic acid (pIC) in Atlantic cod macrophages. Macrophage samples isolated from Atlantic cod (n=3) and treated with pIC or phosphate buffered saline (PBS control) for 24 and 72h were used for miRNA profiling. Following deep sequencing, DESeq2 analyses identified four (miR-731-3p, miR-125b-3-3p, miR-150-3p and miR-462-3p) and two (miR-2188-3p and miR-462-3p) significantly differentially expressed miRNAs at 24 and 72h post-stimulation (HPS), respectively. Sequencing-identified miRNAs were subjected to qPCR validation using a larger number of biological replicates (n=6) exposed to pIC or PBS over time (i.e. 12, 24, 48 and 72 HPS). As in sequencing, miR-731-3p, miR-462-3p and miR-2188-3p showed significant up-regulation by pIC. The sequencing results were not qPCR-validated for miR-125b-3-3p and miR-150-3p as up- and down-regulated miRNAs at 24 HPS, respectively; however, qPCR results showed significant up-regulation in response to pIC stimulation at later time points (i.e. 48 and/or 72 HPS). We also used qPCR to assess the expression of other miRNAs that were previously shown as immune responsive in other vertebrates. qPCR results at 48 and/or 72 HPS revealed that miR-128-3-5p, miR-214-1-5p and miR-451-3p were induced by pIC, whereas miR-30b-3p and miR-199-1-3p expression were repressed in response to pIC. The present study identified ten pIC-stimulated miRNAs, suggesting them as important in antiviral immune responses of Atlantic cod macrophages. Some pIC-responsive miRNAs identified in this study were predicted to target putative immune-related genes of Atlantic cod (e.g. miR-30b-3p targeting herc4), although the regulatory functions of these miRNAs need to be validated by future studies.

Omics and cytokine discovery in fish: Presenting the Yellowtail kingfish (Seriola lalandi) as a case study

A continued programme of research is essential to overcome production bottlenecks in any aquacultured fish species. Since the introduction of genetic and molecular techniques, the quality of immune research undertaken in fish has greatly improved. Thousands of species specific cytokine genes have been discovered, which can be used to conduct more sensitive studies to understand how fish physiology is affected by aquaculture environments or disease. Newly available transcriptomic technologies, make it increasingly easier to study the immunogenetics of farmed species for which little data exists. This paper reviews how the application of transcriptomic procedures such as RNA Sequencing (RNA-Seq) can advance fish research. As a case study, we present some preliminary findings using RNA-Seq to identify cytokine related genes in Seriola lalandi. These will allow in-depth investigations to understand the immune responses of these fish in response to environmental change or disease and help in the development of therapeutic approaches.

Transcriptome profiling of antiviral immune and dietary fatty acid dependent responses of Atlantic salmon macrophage-like cells

BACKGROUND

Due to the limited availability and high cost of fish oil in the face of increasing aquaculture production, there is a need to reduce usage of fish oil in aquafeeds without compromising farm fish health. Therefore, the present study was conducted to determine if different levels of vegetable and fish oils can alter antiviral responses of salmon macrophage-like cells (MLCs). Atlantic salmon (Salmo salar) were fed diets containing 7.4% (FO7) or 5.1% (FO5) fish oil. These diets were designed to be relatively low in EPA + DHA (i.e. FO7: 1.41% and FO5: 1%), but near the requirement level, and resulting in comparable growth. Vegetable oil (i.e. rapeseed oil) was used to balance fish oil in experimental diets. After a 16-week feeding trial, MLCs isolated from fish in these dietary groups were stimulated by a viral mimic (dsRNA: pIC) for 6 h (qPCR assay) and 24 h (microarray and qPCR assays).

RESULTS

The fatty acid composition of head kidney leukocytes varied between the two dietary groups (e.g. higher 20:5n-3 in the FO7 group). Following microarray assays using a 44K salmonid platform, Rank Products (RP) analysis showed 14 and 54 differentially expressed probes (DEP) (PFP < 0.05) between the two diets in control and pIC groups (FO5 vs. FO7), respectively. Nonetheless, Significance Analysis of Microarrays (SAM, FDR < 0.05) identified only one DEP between pIC groups of the two diets. Moreover, we identified a large number (i.e. 890 DEP in FO7 and 1128 DEP in FO5 overlapping between SAM and RP) of pIC-responsive transcripts, and several of them were involved in TLR-/RLR-dependent and cytokine-mediated pathways. The microarray results were validated as significantly differentially expressed by qPCR assays for 2 out of 9 diet-responsive transcripts and for all of the 35 selected pIC-responsive transcripts.

CONCLUSION

Fatty acid-binding protein adipocyte (fabp4) and proteasome subunit beta type-8 (psmb8) were significantly up- and down-regulated, respectively, in the MLCs of fish fed the diet with a lower level of fish oil, suggesting that they are important diet-responsive, immune-related biomarkers for future studies. Although the different levels of dietary fish and vegetable oils involved in this study affected the expression of some transcripts, the immune-related pathways and functions activated by the antiviral response of salmon MLCs in both groups were comparable overall. Moreover, the qPCR revealed transcripts responding early to pIC (e.g. lgp2, map3k8, socs1, dusp5 and cflar) and time-responsive transcripts (e.g. scarb1-a, csf1r, traf5a, cd80 and ctsf) in salmon MLCs. The present study provides a comprehensive picture of the putative molecular pathways (e.g. RLR-, TLR-, MAPK- and IFN-associated pathways) activated by the antiviral response of salmon MLCs.

CXCR3.1 and CXCR3.2 Differentially Contribute to Macrophage Polarization in Teleost Fish

The study of multiple copies of chemokine receptor genes in various teleosts has long appealed to investigators seeking to understand the evolution of the immune system. The CXCR CXCR3 gene has two isoforms, CXCR3.1 and CXCR3.2, which are both expressed in macrophages. The distinct roles of teleost CXCR3s have not been identified previously. In this article, we found that CXCR3.1 and CXCR3.2 differentially contributed to macrophage polarization in the teleosts: ayu (Plecoglossus altivelis), grass carp (Ctenopharyngodon idella), and spotted green pufferfish (Tetraodon nigroviridis). In ayu macrophages, the P. altivelis CXCR3.1 (PaCXCR3.1) gene was constitutively expressed, whereas the P. altivelis CXCR3.2 (PaCXCR3.2) gene was induced postinfection with Escherichia coli Upon E. coli infection, PaCXCR3.1⁺ and PaCXCR3.2⁺ macrophages showed an M1 and an M2 phenotype, respectively. CXCL9-11-like proteins mediated M1 and M2 polarization by interacting with the PaCXCR3.1 and PaCXCR3.2 proteins on macrophages, respectively. The transcription factors P. altivelis STAT1 and P. altivelis STAT3 were activated in PaCXCR3.1⁺ and PaCXCR3.2⁺ macrophages, respectively. Furthermore, the prognosis of septic ayu adoptively transferred with PaCXCR3.2⁺ macrophages was improved. Our data reveal a previously unknown mechanism for macrophage polarization, suggesting that redundant genes may regulate crucial functions in the teleost immune system.

The dietary replacement of marine ingredients by terrestrial animal and plant alternatives modulates the antiviral immune response of Atlantic salmon (Salmo salar)

The effects of replacing marine ingredients by terrestrial ingredients on the health of Atlantic salmon (Salmo salar) are poorly understood. During a 14-week trial, Atlantic salmon fed a fish meal-fish oil based diet (MAR) showed similar growth performance to others fed a plant protein/vegetable oil based diet (VEG), whereas poorer performance was observed in those fed an animal by-product meal/vegetable oil based diet (ABP). At the end of the trial, salmon were injected with either phosphate-buffered saline (PBS) or the viral mimic polyriboinosinic polyribocytidylic acid (pIC) and sampled for head kidney RNA after 24 h. The levels of 27 immune-related transcripts, and of 5 others involved in eicosanoid synthesis (including paralogues in both cases) were measured in the head kidney of the salmon using qPCR. All of the assayed immune-related genes and cox2 were pIC-induced, while the other eicosanoid synthesis-related genes were pIC-repressed. Linear regression was used to establish correlations between different immune transcripts, elucidating the cascade of responses to pIC and specialization among paralogues. Regarding the effect of diet on the antiviral immune response, pIC-treated fish fed diets ABP and VEG showed higher transcript levels of tlr3, irf1b, stat1a, isg15b, and gig1 compared to those fed diet MAR. We infer that the observed dietary immunomodulation could be due to the lower proportion of arachidonic acid (ARA), eicosapentaenoic acid (EPA), and docosahexaenoic acid (DHA) in diets ABP and VEG. Furthermore, our results suggest a major role of dietary ARA in Atlantic salmon immunity, as low ARA proportion in diet VEG coincided with the highest pIC-induction of some immune transcripts (tlr7, stat1c, mxb, and gig1) and the lowest levels of transcripts encoding eicosanoid-synthesizing enzymes (5loxa, 5loxb, and pgds). In contrast, the high ARA/EPA ratio of diet ABP appeared to favor increased expression of transcripts involved in the synthesis of pro-inflammatory eicosanoids (5loxa and 5loxb) and chemotaxis (ccl19b). In conclusion, our findings show that nutritionally balanced plant-based diets may enhance the immune response of Atlantic salmon. Future studies should explore the possible advantages of plant-based diets in Atlantic salmon exposed to a viral infection.

Transcriptome profiling of spleen provides insights into the antiviral mechanism in Schizothorax prenanti after poly (I: C) challenge

Schizothorax prenanti (S. prenanti) is an important economical cold-water fish species in southwestern China, but it is susceptible to various pathogens infection. In order to clearly elucidate the antiviral mechanism, in this study, we have analyzed the transcriptome of S. prenanti spleen after challenge with the virus mimic, poly (I:C) (pIC), using next generation sequencing technology (RNA-seq). A total of 313 differential expressed genes (DEGs) in spleen at 12 h were obtained after pIC treatment, including 268 significantly up-regulated unigenes (fold change > 2) and 45 significantly down-regulated unigenes (fold change > 2). Through the immune-related DEGs (IRDs) screening, 47 IRDs were used to establish heat map, which intuitively showed a significantly difference after pIC treatment. To validate the RNA-seq data and observe gene expression, the expression levels of 14 IRDs were detected by qPCR after pIC treatment at 0, 4, 8, 12, and 24 h. The results indicated that the qPCR data presented a positive line correlation with RNA-seq data, and the 14 IRDs were responsive to pIC stimulation except IL-1β. Thus, based on the RNA-seq and qPCR data, we inferred that MDA5- and Jak-mediated signaling pathways may involve in the antiviral signaling transduction, and induce type I IFNs and ISGs to block virus invasion, respectively. Unfortunately, TLR3 and TLR22, as receptors of virus dsRNA, were no significantly expressed in this study. Nonetheless, our study still provides useful mRNA sequences of antiviral immunity for further immunological research, and facilitates improving disease restriction in S. prenanti.