Functional Genomic Analysis of the Impact of Camelina (Camelina sativa) Meal on Atlantic Salmon (Salmo salar) Distal Intestine Gene Expression and Physiology

Gill and Liver Transcript Expression Changes Associated With Gill Damage in Atlantic Salmon (Salmo salar)

Gill damage represents a significant challenge in the teleost fish aquaculture industry globally, due to the gill’s involvement in several vital functions and direct contact with the surrounding environment. To examine the local and systemic effects accompanying gill damage (which is likely to negatively affect gill function) of Atlantic salmon, we performed a field sampling to collect gill and liver tissue after several environmental insults (e.g., harmful algal blooms). Before sampling, gills were visually inspected and gill damage was scored; gill scores were assigned from pristine [gill score 0 (GS0)] to severely damaged gills (GS3). Using a 44K salmonid microarray platform, we aimed to compare the transcriptomes of pristine and moderately damaged (i.e., GS2) gill tissue. Rank Products analysis (5% percentage of false-positives) identified 254 and 34 upregulated and downregulated probes, respectively, in GS2 compared with GS0. Differentially expressed probes represented genes associated with functions including gill remodeling, wound healing, and stress and immune responses. We performed gill and liver qPCR for all four gill damage scores using microarray-identified and other damage-associated biomarker genes. Transcripts related to wound healing (e.g., neb and klhl41b) were significantly upregulated in GS2 compared with GS0 in the gills. Also, transcripts associated with immune and stress-relevant pathways were dysregulated (e.g., downregulation of snaclec 1-like and upregulation of igkv3) in GS2 compared with GS0 gills. The livers of salmon with moderate gill damage (i.e., GS2) showed significant upregulation of transcripts related to wound healing (i.e., chtop), apoptosis (e.g., bnip3l), blood coagulation (e.g., f2 and serpind1b), transcription regulation (i.e., pparg), and stress-responses (e.g., cyp3a27) compared with livers of GS0 fish. We performed principal component analysis (PCA) using transcript levels for gill and liver separately. The gill PCA showed that PC1 significantly separated GS2 from all other gill scores. The genes contributing most to this separation were pgam2, des, neb, tnnt2, and myom1. The liver PCA showed that PC1 significantly separated GS2 from GS0; levels of hsp70, cyp3a27, pparg, chtop, and serpind1b were the highest contributors to this separation. Also, hepatic acute phase biomarkers (e.g., serpind1b and f2) were positively correlated to each other and to gill damage. Gill damage-responsive biomarker genes and associated qPCR assays arising from this study will be valuable in future research aimed at developing therapeutic diets to improve farmed salmon welfare.

Growth Performance, Growth-Related Genes, Digestibility, Digestive Enzyme Activity, Immune and Stress Responses of de novo Camelina Meal in Diets of Red Seabream (Pagrus major)

Simple Summary

Fish meal (FM) is the major protein source in aquafeed to achieve sustainable aquaculture production. However, the supply of FM is low due to high cost and low availability. There is ongoing conceited research to identify alternative viable protein sources to replace the finite FM. Novel camelina meal (CM) is a plant protein with high amino acids levels and has been tested as an alternative protein source in livestock feeding. However, there is limited information on supplementation of the CM in aquaculture diets. In this study, four diets were formulated to contain 0% plant protein, 205 g/kg soybean meal, and two levels of CM (205 g/kg and 330 g/kg) in diets fed to red seabream. Results indicated that adding CM 205 g/kg in diets-maintained growth performance, digestive enzyme activities, and nutrient digestibility regulated the immunity and stress resistance and modulated the growth-related genes in red seabream. These findings provide the first step in using novel CM and are essential for future practical formulations of feed for red seabream and other marine fish species.

Abstract

A 60-day experiment was designed to assess the effect of different ratios of fish meal (FM): camelina meal plant protein (CM) on growth response and relative gene expression of growth-promoting factors, feed utilization potency, digestive enzymes activities, apparent digestibility (ADC), stress response, non-specific immunity of Pagrus major. Four isonitrogenous (490.7 g/kg of crude protein) and isolipidic (91.5 g/kg total lipid) experimental diets were formulated and designated as camelina meal (CM0), soyabean meal (SBM20.5), CM20.5, and CM33 based on protein contents. At the end of the feed trial, significantly higher (p < 0.05) weight gain, specific growth rate, and feed intake but lower feed conversion ratio were recorded in fish fed CM0, SBM20.5, and CM20.5 than fish fed CM33. The lowest growth, feed utilization, enzyme activity, and digestibility were recorded in fish fed CM33. Significantly higher pepsin, amylase, and protease activities were observed in fish fed CM0, SBM20.5, and CM20.5 diets than fish fed CM33. The highest ADC of protein was recorded in fish fed CM0, SBM20.5, and CM20.5 diets. Hematocrit levels were depressed CM33 while total serum protein, total cholesterol, triglyceride, blood urea nitrogen, total bilirubin, aspartate aminotransferase, and alanine aminotransferase were not significantly changed by the inclusion of CM. Non-specific immune variables (lysozyme activity, peroxidase activity in serum and nitro blue tetrazolium) in fish fed CM0, SBM20.5, and CM20.5 were significantly higher than in fish fed CM33 diet. The superoxide dismutase of fish fed CM20.5 was not significantly different from CM0 and SBM20.5 (p > 0.05). Catalase and low salinity stress test show that CM0, SBM20.5, and CM20.5 were not significantly (p > 0.05) different, while CM33 was significantly lower than the rest of the diets. TBARs show that CM20.5 and CM33 diets were significantly different (p < 0.05), but CM20.5 was not significantly different from SBM20.5. Significantly higher hepatic IGF-1 and IGF-2 mRNA expression was found in fish-fed diet groups CM0, SBM20.5, and CM20.5 than fish fed CM33. The present study indicated that the addition of CM up 205 kg/kg to diet maintains growth, digestive enzymes, nutrient digestibility, immunity, stress resistance, and feed utilization efficiency of red sea bream.

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.

Phenotypic Examination of Camelina sativa (L.) Crantz Accessions from the USDA-ARS National Genetics Resource Program

Camelina sativa (L.) Crntz. is a hardy self-pollinated oilseed plant that belongs to the Brassicaceae family; widely grown throughout the northern hemisphere until the 1940s for production of vegetable oil but was later displaced by higher-yielding rapeseed and sunflower crops. However, interest in camelina as an alternative oil source has been renewed due to its high oil content that is rich in polyunsaturated fatty acids, antioxidants as well as its ability to grow on marginal lands with minimal requirements. For this reason, our group decided to screen the existing (2011) National Genetic Resources Program (NGRP) center collection of camelina for its genetic diversity and provide a phenotypic evaluation of the cultivars available. Properties evaluated include seed and oil traits, developmental and mature morphologies, as well as chromosome content. Selectable marker genes were also evaluated for potential use in biotech manipulation. Data is provided in a raw uncompiled format to allow other researchers to analyze the unbiased information for their own studies. Our evaluation has determined that the NGRP collection has a wide range of genetic potential for both breeding and biotechnological manipulation purposes. Accessions were identified within the NGRP collection that appear to have desirable seed harvest weight (5.06 g/plant) and oil content (44.1%). Other cultivars were identified as having fatty acid characteristics that may be suitable for meal and/or food use, such as low (<2%) erucic acid content, which is often considered for healthy consumption and ranged from a high of 4.79% to a low of 1.83%. Descriptive statistics are provided for a breadth of traits from 41 accessions, as well as raw data, and key seed traits are further explored. Data presented is available for public use.

Assessing Diversity in the Camelina Genus Provides Insights into the Genome Structure of Camelina sativa

Camelina sativa (L.) Crantz an oilseed crop of the Brassicaceae family is gaining attention due to its potential as a source of high value oil for food, feed or fuel. The hexaploid domesticated C. sativa has limited genetic diversity, encouraging the exploration of related species for novel allelic variation for traits of interest. The current study utilized genotyping by sequencing to characterize 193 Camelina accessions belonging to seven different species collected primarily from the Ukrainian-Russian region and Eastern Europe. Population analyses among Camelina accessions with a 2n = 40 karyotype identified three subpopulations, two composed of domesticated C. sativa and one of C. microcarpa species. Winter type Camelina lines were identified as admixtures of C. sativa and C. microcarpa. Eighteen genotypes of related C. microcarpa unexpectedly shared only two subgenomes with C. sativa, suggesting a novel or cryptic sub-species of C. microcarpa with 19 haploid chromosomes. One C. microcarpa accession (2n = 26) was found to comprise the first two subgenomes of C. sativa suggesting a tetraploid structure. The defined chromosome series among C. microcarpa germplasm, including the newly designated C. neglecta diploid née C. microcarpa, suggested an evolutionary trajectory for the formation of the C. sativa hexaploid genome and re-defined the underlying subgenome structure of the reference genome.

Liver Transcriptome Profiling Reveals That Dietary DHA and EPA Levels Influence Suites of Genes Involved in Metabolism, Redox Homeostasis, and Immune Function in Atlantic Salmon (Salmo salar)

The optimal dietary requirement of omega-3 long-chain polyunsaturated fatty acids (ω3 LC-PUFA), namely docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA), for Atlantic salmon that promotes growth and health warrants careful investigation. We used 44K microarrays to study the influence of increasing levels of dietary DHA + EPA (0, 1.0, and 1.4% of the diet, as formulated) in the presence of high linoleic acid (LA) on Atlantic salmon growth and liver transcriptome. After a 14-week feeding trial, Atlantic salmon fed diet ω3LC0 (i.e. 0% of DHA + EPA) showed significantly lower final weight and weight gain, and higher feed conversion ratio compared with ω3LC1.0 and ω3LC1.4 diet groups. The microarray experiment identified 55 and 77 differentially expressed probes (Rank Products analyses; PFP < 10%) in salmon fed diets ω3LC1.4 and ω3LC1.0 compared with those fed diet ω3LC0, respectively. The comparison between ω3LC1.4 and ω3LC1.0 revealed 134 differentially expressed probes. The microarray results were confirmed by qPCR analyses of 22 microarray-identified transcripts. Several key genes involved in fatty acid metabolism including LC-PUFA synthesis were upregulated in fish fed ω3LC0 compared with both other groups. Hierarchical clustering and linear regression analyses of liver qPCR and fatty acid composition data demonstrated significant correlations. In the current study, 1.0% ω3 LC-PUFA seemed to be the minimum requirement for Atlantic salmon based on growth performance; however, multivariate statistical analyses (PERMANOVA and SIMPER) showed that fish fed ω3LC1.0 and ω3LC1.4 diets had similar hepatic fatty acid profiles but marked differences in the transcript expression of biomarker genes involved in redox homeostasis (mgst1), immune responses (mxb, igmb, irf3, lect2a, srk2, and lyz2), and LC-PUFA synthesis (srebp1, fadsd5, and elovl2). This research has provided new insights into dietary requirement of DHA and EPA and their impact on physiologically important pathways in addition to lipid metabolism in Atlantic salmon.

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.

The interbranchial lymphoid tissue likely contributes to immune tolerance and defense in the gills of Atlantic salmon

Central and peripheral immune tolerance is together with defense mechanisms a hallmark of all lymphoid tissues. In fish, such tolerance is especially important in the gills, where the intimate contact between gill tissue and the aqueous environment would otherwise lead to continual immune stimulation by innocuous antigens. In this paper, we focus on the expression of genes associated with immune regulation by the interbranchial lymphoid tissue (ILT) in an attempt to understand its role in maintaining immune homeostasis. Both healthy and virus-challenged fish were investigated, and transcript levels were examined from laser-dissected ILT, gills, head kidney and intestine. Lack of Aire expression in the ILT excluded its involvement in central tolerance and any possibility of its being an analogue to the thymus. On the other hand, the ILT appears to participate in peripheral immune tolerance due to its relatively high expression of forkhead box protein 3 (Foxp3) and other genes associated with regulatory T cells (Tregs) and immune suppression.

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.

De novo assembly of Sockeye salmon kidney transcriptomes reveal a limited early response to piscine reovirus with or without infectious hematopoietic necrosis virus superinfection

Background

Piscine reovirus (PRV) has been associated with the serious disease known as Heart and Skeletal Muscle Inflammation (HSMI) in cultured Atlantic salmon Salmo salar in Norway. PRV is also prevalent in wild and farmed salmon without overt disease manifestations, suggesting multifactorial triggers or PRV variant-specific factors are required to initiate disease. In this study, we explore the head kidney transcriptome of Sockeye salmon Oncorhynchus nerka during early PRV infection to identify host responses in the absence of disease in hopes of elucidating mechanisms by which PRV may directly alter host functions and contribute to the development of a disease state. We further investigate the role of PRV as a coinfecting agent following superinfection with infectious hematopoietic necrosis virus (IHNV) – a highly pathogenic rhabdovirus endemic to the west coast of North America.

Results

Challenge of Sockeye salmon with PRV resulted in high quantities of viral transcripts to become present in the blood and kidney of infected fish without manifestations of disease. De novo transcriptome assembly of over 2.3 billion paired RNA-seq reads from the head kidneys of 36 fish identified more than 320,000 putative unigenes, of which less than 20 were suggested to be differentially expressed in response to PRV at either 2 or 3 weeks post challenge by DESeq2 and edgeR analysis. Of these, only one, Ependymin, was confirmed to be differentially expressed by qPCR in an expanded sample set. In contrast, IHNV induced substantial transcriptional changes (differential expression of > 20,000 unigenes) which included transcripts involved in antiviral and inflammatory response pathways. Prior infection with PRV had no significant effect on host responses to superinfecting IHNV, nor did host responses initiated by IHNV exposure influence increasing PRV loads.

Conclusions

PRV does not substantially alter the head kidney transcriptome of Sockeye salmon during early (2 to 3 week) infection and dissemination in a period of significant increasing viral load, nor does the presence of PRV change the host transcriptional response to an IHNV superinfection. Further, concurrent infections of PRV and IHNV do not appear to significantly influence the infectivity or severity of IHNV associated disease, or conversely, PRV load.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3196-y) contains supplementary material, which is available to authorized users.