Evaluation of the impact of camelina oil-containing diets on the expression of genes involved in the innate anti-viral immune response in Atlantic cod (Gadus morhua)

Genetic variation and structural diversity in major seed proteins among and within Camelina species

Genetic variation in seed protein composition, seed protein gene expression and predictions of seed protein physiochemical properties were documented in C. sativa and other Camelina species. Seed protein diversity was examined in six Camelina species (C. hispida, C. laxa, C. microcarpa, C. neglecta, C. rumelica and C. sativa). Differences were observed in seed protein electrophoretic profiles, total seed protein content and amino acid composition between the species. Genes encoding major seed proteins (cruciferins, napins, oleosins and vicilins) were catalogued for C. sativa and RNA-Seq analysis established the expression patterns of these and other genes in developing seed from anthesis through to maturation. Examination of 187 C. sativa accessions revealed limited variation in seed protein electrophoretic profiles, though sufficient to group the majority into classes based on high MW protein profiles corresponding to the cruciferin region. C. sativa possessed four distinct types of cruciferins, named CsCRA, CsCRB, CsCRC and CsCRD, which corresponded to orthologues in Arabidopsis thaliana with members of each type encoded by homeologous genes on the three C. sativa sub-genomes. Total protein content and amino acid composition varied only slightly; however, RNA-Seq analysis revealed that CsCRA and CsCRB genes contributed > 95% of the cruciferin transcripts in most lines, whereas CsCRC genes were the most highly expressed cruciferin genes in others, including the type cultivar DH55. This was confirmed by proteomics analyses. Cruciferin is the most abundant seed protein and contributes the most to functionality. Modelling of the C. sativa cruciferins indicated that each type possesses different physiochemical attributes that were predicted to impart unique functional properties. As such, opportunities exist to create C. sativa cultivars with seed protein profiles tailored to specific technical applications.

Interacting Effects of Sea Louse (Lepeophtheirus salmonis) Infection and Formalin-Killed Aeromonas salmonicida on Atlantic Salmon Skin Transcriptome

Lepeophtheirus salmonis (sea lice) and bacterial co-infection threatens wild and farmed Atlantic salmon performance and welfare. In the present study, pre-adult L. salmonis-infected and non-infected salmon were intraperitoneally injected with either formalin-killed Aeromonas salmonicida bacterin (ASAL) or phosphate-buffered saline (PBS). Dorsal skin samples from each injection/infection group (PBS/no lice, PBS/lice, ASAL/no lice, and ASAL/lice) were collected at 24 h post-injection and used for transcriptome profiling using a 44K salmonid microarray platform. Microarray results showed no clear inflammation gene expression signatures and revealed extensive gene repression effects by pre-adult lice (2,189 down and 345 up-regulated probes) in the PBS-injected salmon (PBS/lice vs. PBS/no lice), which involved basic cellular (e.g., RNA and protein metabolism) processes. Lice repressive effects were not observed within the group of ASAL-injected salmon (ASAL/lice vs. ASAL/no lice); on the contrary, the observed skin transcriptome changes –albeit of lesser magnitude (82 up and 1 down-regulated probes)– suggested the activation in key immune and wound healing processes (e.g., neutrophil degranulation, keratinocyte differentiation). The molecular skin response to ASAL was more intense in the lice-infected (ASAL/lice vs. PBS/lice; 272 up and 11 down-regulated probes) than in the non-infected fish (ASAL/no lice vs. PBS/no lice; 27 up-regulated probes). Regardless of lice infection, the skin’s response to ASAL was characterized by the putative activation of both antibacterial and wound healing pathways. The transcriptomic changes prompted by ASAL+lice co-stimulation (ASAL/lice vs. PBS/no lice; 1878 up and 3120 down-regulated probes) confirmed partial mitigation of lice repressive effects on fundamental cellular processes and the activation of pathways involved in innate (e.g., neutrophil degranulation) and adaptive immunity (e.g., antibody formation), as well as endothelial cell migration. The qPCR analyses evidenced immune-relevant genes co-stimulated by ASAL and lice in an additive (e.g., mbl2b, bcl6) and synergistic (e.g., hampa, il4r) manner. These results provided insight on the physiological response of the skin of L. salmonis-infected salmon 24 h after ASAL stimulation, which revealed immunostimulatory properties by the bacterin with potential applications in anti-lice treatments for aquaculture. As a simulated co-infection model, the present study also serves as a source of candidate gene biomarkers for sea lice and bacterial co-infection.

A Multidisciplinary Approach Evaluating Soybean Meal-Induced Enteritis in Rainbow Trout Oncorhynchus mykiss

This study evaluated a diverse range of markers of feeding stress to obtain a more precise assessment of the welfare of rainbow trout in relation to inadequate husbandry conditions. A feeding stress model based on dietary soybean meal was employed to identify suitable minimally invasive “classical” stress markers, together with molecular signatures. In a 56-day feeding experiment, rainbow trout were fed diets containing different levels of soybean meal. The impact of these different soybean meal diets on rainbow trout was assessed by water quality analyses, clinical health observations, classic growth and performance parameters, gut histopathology, blood-parameter measurements and multigene-expression profiling in RNA from whole blood. Soybean meal-induced enteritis was manifested phenotypically by an inflammatory reaction in the posterior section of the intestine and by diarrhoea in some trout. These inflammatory changes were associated with decreased supranuclear vacuolation. The haematocrit values and the levels of plasma cortisol and circulating lymphocytes in the blood were increased in trout that had consumed high amounts of SBM. Notably, the increased haematocrit depended significantly on the bodyweight of the individual trout. The transcript levels of certain genes (e.g., MAP3K1, LYG, NOD1, STAT1 and HSP90AB) emerged as potentially useful indicators in the blood of rainbow trout providing valuable information about inadequate nutrition. The expression-profiling findings provide a basis for improved, minimally invasive monitoring of feeding regimens in trout farming and may stimulate the development of practical detection devices for innovative aquaculture operations.

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.

Transcriptome Profiling of Atlantic Salmon Adherent Head Kidney Leukocytes Reveals That Macrophages Are Selectively Enriched During Culture

The Atlantic salmon (Salmo salar) is an economically important fish, both in aquaculture and in the wild. In vertebrates, macrophages are some of the first cell types to respond to pathogen infection and disease. While macrophage biology has been characterized in mammals, less is known in fish. Our previous work identified changes in the morphology, phagocytic ability, and miRNA profile of Atlantic salmon adherent head kidney leukocytes (HKLs) from predominantly “monocyte-like” at Day 1 of in vitro culture to predominantly “macrophage-like” at Day 5 of culture. Therefore, to further characterize these two cell populations, we examined the mRNA transcriptome profile in Day 1 and Day 5 HKLs using a 44K oligonucleotide microarray. Large changes in the transcriptome were revealed, including changes in the expression of macrophage and immune-related transcripts (e.g. csf1r, arg1, tnfa, mx2), lipid-related transcripts (e.g. fasn, dhcr7, fabp6), and transcription factors involved in macrophage differentiation and function (e.g. klf2, klf9, irf7, irf8, stat1). The in silico target prediction analysis of differentially expressed genes (DEGs) using miRNAs known to change expression in Day 5 HKLs, followed by gene pathway enrichment analysis, supported that these miRNAs may be involved in macrophage maturation by targeting specific DEGs. Elucidating how immune cells, such as macrophages, develop and function is a key step in understanding the Atlantic salmon immune system. Overall, the results indicate that, without the addition of exogenous factors, the adherent HKL cell population differentiates in vitro to become macrophage-like.

Influence of Varying Dietary ω6 to ω3 Fatty Acid Ratios on the Hepatic Transcriptome, and Association with Phenotypic Traits (Growth, Somatic Indices, and Tissue Lipid Composition), in Atlantic Salmon (Salmo salar)

Simple Summary

Plant oils are routinely used in fish feeds as a fish oil replacement. However, these terrestrial alternatives typically contain high levels of ω6 fatty acids (FA) and, thus, high ω6 to ω3 (ω6:ω3) FA ratios, which influence farmed fish and their consumers. The ω6:ω3 ratio is known to affect many biological processes (e.g., inflammation, FA metabolism) and human diseases; however, its impacts on fish physiology and the underlying molecular mechanisms are less well understood. In this study, we used 44 K microarrays to examine which genes and molecular pathways are altered by variation in dietary ω6:ω3 in Atlantic salmon. Our microarray study showed that several genes related to immune response, lipid metabolism, cell proliferation, and translation were differentially expressed between the two extreme ω6:ω3 dietary treatments. We also revealed that the PPARα activation-related transcript helz2 is a potential novel molecular biomarker of tissue variation in ω6:ω3. Further, correlation analyses illustrated the relationships between liver transcript expression and tissue (liver, muscle) lipid composition, and other phenotypic traits in salmon fed low levels of fish oil. This nutrigenomic study enhanced the current understanding of Atlantic salmon gene expression response to varying dietary ω6:ω3.

Abstract

The importance of dietary omega-6 to omega-3 (ω6:ω3) fatty acid (FA) ratios for human health has been extensively examined. However, its impact on fish physiology, and the underlying molecular mechanisms, are less well understood. This study investigated the influence of plant-based diets (12-week exposure) with varying ω6:ω3 (0.4–2.7) on the hepatic transcriptome of Atlantic salmon. Using 44 K microarray analysis, genes involved in immune and inflammatory response (lect2a, itgb5, helz2a, p43), lipid metabolism (helz2a), cell proliferation (htra1b), control of muscle and neuronal development (mef2d) and translation (eif2a, eif4b1, p43) were identified; these were differentially expressed between the two extreme ω6:ω3 dietary treatments (high ω6 vs. high ω3) at week 12. Eight out of 10 microarray-identified transcripts showed an agreement in the direction of expression fold-change between the microarray and qPCR studies. The PPARα activation-related transcript helz2a was confirmed by qPCR to be down-regulated by high ω6 diet compared with high ω3 diet. The transcript expression of two helz2 paralogues was positively correlated with ω3, and negatively with ω6 FA in both liver and muscle, thus indicating their potential as biomarkers of tissue ω6:ω3 variation. Mef2d expression in liver was suppressed in the high ω6 compared to the balanced diet (ω6:ω3 of 2.7 and 0.9, respectively) fed fish, and showed negative correlations with ω6:ω3 in both tissues. The hepatic expression of two lect2 paralogues was negatively correlated with viscerosomatic index, while htra1b correlated negatively with salmon weight gain and condition factor. Finally, p43 and eif2a were positively correlated with liver Σω3, while these transcripts and eif4b2 showed negative correlations with 18:2ω6 in the liver. This suggested that some aspects of protein synthesis were influenced by dietary ω6:ω3. In summary, this nutrigenomic study identified hepatic transcripts responsive to dietary variation in ω6:ω3, and relationships of transcript expression with tissue (liver, muscle) lipid composition and other phenotypic traits.

The transcriptomic responses of Atlantic salmon (Salmo salar) to high temperature stress alone, and in combination with moderate hypoxia

BACKGROUND

Increases in ocean temperatures and in the frequency and severity of hypoxic events are expected with climate change, and may become a challenge for cultured Atlantic salmon and negatively affect their growth, immunology and welfare. Thus, we examined how an incremental temperature increase alone (Warm & Normoxic-WN: 12 → 20 °C; 1 °C week- 1), and in combination with moderate hypoxia (Warm & Hypoxic-WH: ~ 70% air saturation), impacted the salmon's hepatic transcriptome expr\ession compared to control fish (CT: 12 °C, normoxic) using 44 K microarrays and qPCR.

RESULTS

Overall, we identified 2894 differentially expressed probes (DEPs, FDR < 5%), that included 1111 shared DEPs, while 789 and 994 DEPs were specific to WN and WH fish, respectively. Pathway analysis indicated that the cellular mechanisms affected by the two experimental conditions were quite similar, with up-regulated genes functionally associated with the heat shock response, ER-stress, apoptosis and immune defence, while genes connected with general metabolic processes, proteolysis and oxidation-reduction were largely suppressed. The qPCR assessment of 41 microarray-identified genes validated that the heat shock response (hsp90aa1, serpinh1), apoptosis (casp8, jund, jak2) and immune responses (apod, c1ql2, epx) were up-regulated in WN and WH fish, while oxidative stress and hypoxia sensitive genes were down-regulated (cirbp, cyp1a1, egln2, gstt1, hif1α, prdx6, rraga, ucp2). However, the additional challenge of hypoxia resulted in more pronounced effects on heat shock and immune-related processes, including a stronger influence on the expression of 14 immune-related genes. Finally, robust correlations between the transcription of 19 genes and several phenotypic traits in WH fish suggest that changes in gene expression were related to impaired physiological and growth performance.

CONCLUSION

Increasing temperature to 20 °C alone, and in combination with hypoxia, resulted in the differential expression of genes involved in similar pathways in Atlantic salmon. However, the expression responses of heat shock and immune-relevant genes in fish exposed to 20 °C and hypoxia were more affected, and strongly related to phenotypic characteristics (e.g., growth). This study provides valuable information on how these two environmental challenges affect the expression of stress-, metabolic- and immune-related genes and pathways, and identifies potential biomarker genes for improving our understanding of fish health and welfare.

Influence of Dietary Long-Chain Polyunsaturated Fatty Acids and ω6 to ω3 Ratios on Head Kidney Lipid Composition and Expression of Fatty Acid and Eicosanoid Metabolism Genes in Atlantic Salmon (Salmo salar)

The interaction of dietary eicosapentaenoic acid and docosahexaenoic acid (EPA+DHA) levels with omega-6 to omega-3 ratios (ω6:ω3), and their impact on head kidney lipid metabolism in farmed fish, are not fully elucidated. We investigated the influence of five plant-based diets (12-week exposure) with varying EPA+DHA levels (0.3, 1.0, or 1.4%) and ω6:ω3 (high ω6, high ω3, or balanced) on tissue lipid composition, and transcript expression of genes involved in fatty acid and eicosanoid metabolism in Atlantic salmon head kidney. Tissue fatty acid composition was reflective of the diet with respect to C₁₈ PUFA and MUFA levels (% of total FA), and ω6:ω3 (0.5–1.5). Fish fed 0.3% EPA+DHA with high ω6 (0.3% EPA+DHA↑ω6) had the highest increase in proportions (1.7–2.3-fold) and in concentrations (1.4-1.8-fold) of arachidonic acid (ARA). EPA showed the greatest decrease in proportion and in concentration (by ~½) in the 0.3% EPA+DHA↑ω6 fed fish compared to the other treatments. However, no differences were observed in EPA proportions among salmon fed the high ω3 (0.3 and 1.0% EPA+DHA) and balanced (1.4% EPA+DHA) diets, and DHA proportions were similar among all treatments. Further, the transcript expression of elovl5a was lowest in the 0.3% EPA+DHA↑ω6 fed fish, and correlated positively with 20:3ω3, 20:4ω3 and EPA:ARA in the head kidney. This indicates that high dietary 18:3ω3 promoted the synthesis of ω3 LC-PUFA. Dietary EPA+DHA levels had a positive impact on elovl5a, fadsd5 and srebp1 expression, and these transcripts positively correlated with tissue ΣMUFA. This supported the hypothesis that LC-PUFA synthesis is positively influenced by tissue MUFA levels in Atlantic salmon. The expression of pparaa was higher in the 0.3% EPA+DHA↑ω6 compared to the 0.3% EPA+DHA↑ω3 fed fish. Finally, significant correlations between head kidney fatty acid composition and the expression of eicosanoid synthesis-related transcripts (i.e., 5loxa, 5loxb, cox1, cox2, ptges2, ptges3, and pgds) illustrated the constitutive relationships among fatty acids and eicosanoid metabolism in 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.

Interaction between dietary fatty acids and genotype on immune response in Atlantic salmon (Salmo salar) after vaccination: A transcriptome study

A pivotal matter to aquaculture is the sourcing of sustainable resources as ingredients to aquafeeds. Levels of plant delivered oils as source of fatty acids (FA) in aquafeeds have reached around 70% resulting in reduced levels of long-chain omega-3 polyunsaturated fatty acids (LC n-3 PUFA), such as eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA), in salmon fillet composition. EPA and DHA can modulate inflammation and immune response, so it is crucial to understand how fish immune response is affected by low LC n-3 PUFA diet and if this diet can have a detrimental effect on vaccine response. Atlantic salmon (Salmo salar) can produce EPA/DHA from α-linolenic acid (ALA) and this endogenous capacity can be explored to develop families with higher tolerance to low LC n-3 PUFA diets. Here we analyze innate and adaptive immune response in Atlantic salmon to a commercial vaccine after being fed low levels of EPA and DHA, and we also compare three strains of salmon selected by their endogenous capacity of synthesizing LC- n-3 PUFA. A total of 2,890 differentially expressed genes (DEGs) were identified (p-value adjusted < 0.1) when comparing vaccinated fish against control non-vaccinated. Gene ontology (GO) and KEGG analysis with 442 up/downregulated genes revealed that most DEGs were both related to immune response as well as part of important immune related pathways, as "Toll-like receptor" and "Cytokine-Cytokine interaction". Adaptive response was also addressed by measuring antigen specific IgM, and titers were significantly higher than in the pre-immune fish at 62 days post-immunization. However, diet and strain had no/little effect on vaccine-specific IgM or innate immune responses. Atlantic salmon therefore display robustness in its response to vaccination even when feed low levels of LC n-3 PUFA.

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.

Nutrigenomics and immune function in fish: new insights from omics technologies

The interplay between nutrition and immune system is well recognised, however the true integration of research between nutrition, animal energy status and immune function is still far from clear. In fish nutrition, especially for species maintained in aquaculture, formulated feeds are significantly different from the natural diet with recent changes in nutrient sources, especially with protein and oil sources now being predominated by terrestrial derived ingredients. Additionally, many feeds are now incorporated to health management and termed functional feeds, which are believed to improve fish health, reduce disease outbreaks and/or improve post-infection recovery. Using new omics technologies, including transcriptomics (microarray and RNA-seq) and proteomics, the impacts of nutrition on the immune system is becoming clearer. By using molecular pathway enrichment analysis, modules of genes can indicate how both local (intestinal) and systemic immune function are being altered. Although great progress has been made to define the changes in host immune function, understanding the interplay between fish nutrition, intestinal microbiome and immune system is only just beginning to emerge.

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.

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 the antiviral immune response in Atlantic cod macrophages

A study was conducted to determine the transcriptome response of Atlantic cod (Gadus morhua) macrophages to the viral mimic, polyriboinosinic polyribocytidylic acid (pIC), using a 20K Atlantic cod microarray platform and qPCR. We identified 285 significantly up-regulated and 161 significantly down-regulated probes in cod macrophages 24 h after pIC stimulation. A subset of 26 microarray-identified transcripts was subjected to qPCR validation using samples treated with pIC or phosphate-buffered saline (control) over time (3, 6, 12, 24, 48 h), and 77% of them showed a significant response to pIC. The microarray and qPCR analyses in this study showed that pIC induced the expression of cod macrophage transcripts involved in RLR- and TLR-dependent pathogen recognition (e.g. tlr3, tlr7, mda5 and lgp2), as well as signal transducers (e.g. stat1 and nfkbia) and transcription activators (e.g. irf7 and irf10) in the MyD88-independent and dependent signalling pathways. Several immune effectors (e.g. isg15s, viperin, herc4, mip2 and ccl13) were significantly up-regulated in pIC-stimulated cod macrophages. The expression of some transcripts (e.g. irf7, irf10, viperin) was significantly up-regulated by pIC as early as 12 h. All pIC-induced transcripts had peak expression at either 24 h (e.g. tlr7, irf7, mip2) or 48 h (e.g. tlr3, lgp2, stat1). This study suggests possible roles of both vertebrate-conserved (e.g. tlr3 as an up-regulated gene) and fish-specific (tlr22g as a down-regulated gene) receptors in dsRNA recognition, and the importance of conserved and potentially fish-specific interferon stimulated genes in cod macrophages.

Transcriptome profiling reveals that feeding wild zooplankton to larval Atlantic cod (Gadus morhua) influences suites of genes involved in oxidation-reduction, mitosis, and selenium homeostasis

BACKGROUND

Larval nutrition and growth are key issues for wild and cultured cod. While it was shown previously that larval cod fed wild zooplankton grow faster than those fed only rotifers, the mechanisms involved in this enhanced growth are not completely understood. We used microarrays to identify larval cod transcripts that respond to feeding with small amounts of wild zooplankton (5-10 % of live prey items). The larval transcriptome was compared between 3 treatment groups [fed rotifers (RA), rotifers with protein hydrolysate (RA-PH), or rotifers with zooplankton (RA-Zoo)] at 9-10 mm length [26-30 days post-hatch (dph)] to identify a robust suite of zooplankton-responsive genes (i.e. differentially expressed between RA-Zoo and both other groups).

RESULTS

The microarray experiment identified 147 significantly up-regulated and 156 significantly down-regulated features in RA-Zoo compared with both RA and RA-PH. Gene ontology terms overrepresented in the RA-Zoo responsive gene set included "response to selenium ion" and several related to cell division and oxidation-reduction. Ten selenoprotein-encoding genes, and 2 genes involved in thyroid hormone generation, were up-regulated in RA-Zoo compared with both other groups. Hierarchical clustering of RA-Zoo responsive genes involved in oxidation-reduction and selenium homeostasis demonstrated that only the zooplankton treatment had a considerable and consistent impact on the expression of these genes. Fourteen microarray-identified genes were selected for QPCR involving 9-13 mm larvae, and 13 of these were validated as differentially expressed between RA-Zoo and both other groups at ~9 mm. In contrast, in age-matched (34-35 dph; ~11 mm RA and RA-PH, ~13 mm RA-Zoo) and size-matched (~13 mm) older larvae, only 2 and 3 genes, respectively, showed the same direction of RA-Zoo-responsive change as in ~9 mm larvae.

CONCLUSIONS

The modulation of genes involved in selenium binding, redox homeostasis, and thyroid hormone generation in ~9 mm RA-Zoo larvae in this study may be in response to the relatively high levels of selenium, iodine, and LC-PUFA (potentially causing oxidative stress) in zooplankton. Nonetheless, only a subset of zooplankton-responsive genes in ~9 mm larvae remained so in older larvae, suggesting that the observed transcriptome changes are largely involved in initiating the period of growth enhancement.

Atlantic salmon (Salmo salar) liver transcriptome response to diets containing Camelina sativa products

Due to increasing demand for fish oil (FO) and fish meal (FM) in aquafeeds, more sustainable alternatives such as plant-derived oils and proteins are needed. Camelina sativa products are viable feed ingredients given the high oil and crude protein content in the seed. Atlantic salmon were fed diets with complete or partial replacement of FO and/or FM with camelina oil (CO) and/or camelina meal (CM) in a 16-week trial [Control diet: FO; Test diets: 100% CO replacement of FO (100CO), or 100CO with solvent-extracted FM (100COSEFM), 10% CM (100CO10CM), or SEFM+10% CM (100COSEFM10CM)]. Diet composition, growth, and fatty acid analyses for this feeding trial were published previously. A 44K microarray experiment identified liver transcripts that responded to 100COSEFM10CM (associated with reduced growth) compared to controls, yielding 67 differentially expressed features (FDR<5%). Ten microarray-identified genes [cpt1, pcb, bar, igfbp-5b (2 paralogues), btg1, dnph1, lect-2, clra, klf9, and fadsd6a], and three additional genes involved in lipid metabolism [elovl2, elovl5 (2 paralogues), and fadsd5], were subjected to QPCR with liver templates from all 5 dietary treatments. Of the microarray-identified genes, only bar was not QPCR validated. Both igfbp-5b paralogues were significantly down-regulated, and fadsd6a was significantly up-regulated, in all 4 camelina-containing diet groups compared with controls. Multivariate statistics were used to correlate hepatic desaturase and elongase gene expression data with tissue fatty acid profiles, indicating the involvement of these genes in LC-PUFA biosynthesis. This nutrigenomic study provides molecular biomarkers for use in developing novel aquafeeds using camelina products.