Effect of plant-based diets with varying ratios of ω6 to ω3 fatty acids on growth performance, tissue composition, fatty acid biosynthesis and lipid-related gene expression in Atlantic salmon (Salmo salar)

Little is known about how variation in omega-6 to omega-3 (ω6:ω3) fatty acid (FA) ratios affects lipid metabolism and eicosanoid synthesis in salmon, and the potential underlying molecular mechanisms. The current study examined the impact of five plant-based diets (12-week exposure) with varying ω6:ω3 (0.3-2.7) on the growth, tissue lipid composition (muscle and liver), and hepatic transcript expression of lipid metabolism and eicosanoid synthesis-related genes in Atlantic salmon. Growth performance and organ indices were not affected by dietary ω6:ω3. The liver and muscle FA composition was highly reflective of the diet (ω6:ω3 of 0.2-0.8 and 0.3-1.9, respectively) and suggested elongation and desaturation of the ω3 and ω6 precursors 18:3ω3 and 18:2ω6. Furthermore, proportions of ω6 and ω3 PUFA in both tissues showed significant positive correlations with dietary inclusion (% of diet) of soy and linseed oils, respectively. Compound-specific stable isotope analysis (CSIA) further demonstrated that liver long-chain polyunsaturated fatty acid (LC-PUFA) synthesis (specifically 20:5ω3 and 20:4ω6) was largely driven by dietary 18:3ω3 and 18:2ω6, even when 20:5ω3 and 22:6ω3 were supplied at levels above minimum requirements. In addition, significant positive and negative correlations were identified between the transcript expression of LC-PUFA synthesis-related genes and liver ω6 and ω3 LC-PUFA, respectively, further supporting FA biosynthesis. Liver ω3 LC-PUFA also correlated negatively with the eicosanoid synthesis-related transcripts pgds and cox1. This is the first study to use CSIA, hepatic transcriptome, and tissue lipid composition analyses concurrently to demonstrate the impact of plant-based diets with varying ω6:ω3 on farmed Atlantic salmon.

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.

Minimizing marine ingredients in diets of farmed Atlantic salmon (Salmo salar): Effects on growth performance and muscle lipid and fatty acid composition

Due to limited fish meal and fish oil resources and their high costs for the aquaculture industry, it is necessary to find alternative sustainable sources of protein and lipids. Therefore, seven different diets were formulated with different levels of animal by-products, vegetable proteins, fish oil and rapeseed oil, to feed farmed Atlantic salmon, and their effects on growth performance, muscle lipid class, and fatty acid composition were examined. Protein sources included anchovy, poultry, feather, blood, corn, soy and wheat. Growth performance indicated that the diet with the lowest fish meal and fish oil content resulted in the lowest weight gain and final weight, followed by the diet containing the highest level of animal by-products. The lipid class analysis showed no statistical difference in the muscle total lipid content using different diets. However, significant statistical differences were observed among the main lipid classes; triacylglycerols, phospholipids, and sterols. The diet containing 1.4% omega-3 long-chain fatty acids resulted in the highest content of triacylglycerols and phospholipids. Diets containing medium and low levels of fish oil and fish meal, respectively, led to as high a level of ω3 fatty acids in muscle as when fish were fed diets with high levels of fish meal and fish oil. The results of this study suggest that feeding a diet containing low levels of fish meal and moderate levels of fish oil does not significantly affect ω3 fatty acid composition in muscle. Fish meal could be reduced to 5% without affecting growth as long as there was a minimum of 5% fish oil, and animal by-products did not exceed 26% of the diet.

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.

Functional Annotation of All Salmonid Genomes (FAASG): an international initiative supporting future salmonid research, conservation and aquaculture

We describe an emerging initiative - the 'Functional Annotation of All Salmonid Genomes' (FAASG), which will leverage the extensive trait diversity that has evolved since a whole genome duplication event in the salmonid ancestor, to develop an integrative understanding of the functional genomic basis of phenotypic variation. The outcomes of FAASG will have diverse applications, ranging from improved understanding of genome evolution, to improving the efficiency and sustainability of aquaculture production, supporting the future of fundamental and applied research in an iconic fish lineage of major societal importance.

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.

Successive Losses of Central Immune Genes Characterize the Gadiformes' Alternate Immunity

Great genetic variability among teleost immunomes, with gene losses and expansions of central adaptive and innate components, has been discovered through genome sequencing over the last few years. Here, we demonstrate that the innate Myxovirus resistance gene (Mx) is lost from the ancestor of Gadiformes and the closely related Stylephorus chordatus, thus predating the loss of Major Histocompatibility Complex class II (MHCII) in Gadiformes. Although the functional implication of Mx loss is still unknown, we demonstrate that this loss is one of several ancient events appearing in successive order throughout the evolution of teleost immunity. In particular, we find that the loss of Toll-like receptor 5 predates the loss of Mx involving the entire Paracanthopterygii lineage. Using a time-calibrated phylogeny, we show that loss of MHCII and Mx overlap with major paleoclimatic and geological events indicating that these genetic changes were adaptive responses to the changing environment at the time.

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.

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.

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

The inclusion of plant meals in diets of farmed Atlantic salmon can elicit inflammatory responses in the distal intestine (DI). For the present work, fish were fed a standard fish meal (FM) diet or a diet with partial replacement of FM with solvent-extracted camelina meal (CM) (8, 16, or 24 % CM inclusion) during a 16-week feeding trial. A significant decrease in growth performance was seen in fish fed all CM inclusion diets (Hixson et al. in Aquacult Nutr 22:615-630, 2016). A 4x44K oligonucleotide microarray experiment was carried out and significance analysis of microarrays (SAM) and rank products (RP) methods were used to identify differentially expressed genes between the DIs of fish fed the 24 % CM diet and those fed the FM diet. Twelve features representing six known transcripts and two unknowns were identified as CM responsive by both SAM and RP. The six known transcripts (including thioredoxin and ependymin), in addition to tgfb, mmp13, and GILT, were studied using qPCR with RNA templates from all four experimental diet groups. All six microarray-identified genes were confirmed to be CM responsive, as was tgfb and mmp13. Histopathological analyses identified signs of inflammation in the DI of salmon fed CM-containing diets, including lamina propria and sub-epithelial mucosa thickening, infiltration of eosinophilic granule cells, increased goblet cells and decreased enterocyte vacuolization. All of these were significantly altered in 24 % CM compared to all other diets, with the latter two also altered in 16 % CM compared with 8 % CM and control diet groups. Significant correlation was seen between histological parameters as well as between five of the qPCR analyzed genes and histological parameters. These molecular biomarkers of inflammation arising from long-term dietary CM exposure will be useful in the development of CM-containing diets that do not have deleterious effects on salmon growth or physiology.

Cold-induced changes in stress hormone and steroidogenic transcript levels in cunner (Tautogolabrus adspersus), a fish capable of metabolic depression

The cunner (Tautogolabrus adspersus) is a fish with a wide latitudinal distribution that is capable of going into metabolic depression during the winter months, and thus, represents a unique model to investigate the impacts of cold temperatures on the stress response. In this study, we measured resting (pre-stress) plasma cortisol levels in 10 °C and 0 °C acclimated cunner from Newfoundland, and both catecholamine and cortisol levels after they were given a standardized handling stress (i.e. 1 min air exposure). In addition, we cloned and characterized cDNAs for several key genes of the cortisol-axis [cytochrome P450scc, steroidogenic acute regulatory protein (StAR) and a glucocorticoid receptor (GR) most likely to be an ortholog of the teleost GR2], determined the tissue distribution of their transcripts, and measured their constitutive (i.e. pre-stress) transcript levels in individuals acclimated to both temperatures. In cunner acclimated to 0 °C, post-stress epinephrine and norepinephrine levels were much lower (by approximately 9- and 5-fold, respectively) compared to 10 °C acclimated fish, and these fish had relatively low resting cortisol levels (~15 ngml(-1)) and showed a typical post-stress response. In contrast, those acclimated to 10 °C had quite high resting cortisol levels (~75 ngml(-1)) that actually decreased (to ~20 ngml(-1)) post-stress before returning to pre-stress levels. Finally, fish acclimated to 10 °C had higher P450scc transcript levels in the head kidney and lower levels of GR transcript in both the head kidney and liver. Taken together, these results suggest that: (1) temperature has a profound effect on the stress response of this species; and (2) although the ancestors of this species inhabited warm waters (i.e. they are members of the family Labridae), populations of cunner from colder regions may show signs of stress at temperatures as low as 10 °C.

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.

Cloning and characterization of aquaglyceroporin genes from rainbow smelt (Osmerus mordax) and transcript expression in response to cold temperature

Aquaglyceroporins (GLPs) are integral membrane proteins that facilitate passive movement of water, glycerol and urea across cellular membranes. In this study, GLP-encoding genes were characterized in rainbow smelt (Osmerus mordax mordax), an anadromous teleost that accumulates high glycerol and modest urea levels in plasma and tissues as an adaptive cryoprotectant mechanism in sub-zero temperatures. We report the gene and promoter sequences for two aqp10b paralogs (aqp10ba, aqp10bb) that are 82% identical at the predicted amino acid level, and aqp9b. Aqp10bb and aqp9b have the 6 exon structure common to vertebrate GLPs. Aqp10ba has 8 exons; there are two additional exons at the 5' end, and the promoter sequence is different from aqp10bb. Molecular phylogenetic analysis suggests that the aqp10b paralogs arose from a gene duplication event specific to the smelt lineage. Smelt GLP transcripts are ubiquitously expressed; however, aqp10ba transcripts were highest in kidney, aqp10bb transcripts were highest in kidney, intestine, pyloric caeca and brain, and aqp9b transcripts were highest in spleen, liver, red blood cells and kidney. In cold-temperature challenge experiments, plasma glycerol and urea levels were significantly higher in cold- compared to warm-acclimated smelt; however, GLP transcript levels were generally either significantly lower or remained constant. The exception was significantly higher aqp10ba transcript levels in kidney. High aqp10ba transcripts in smelt kidney that increase significantly in response to cold temperature in congruence with plasma urea suggest that this gene duplicate may have evolved to allow the re-absorption of urea to concomitantly conserve nitrogen and prevent freezing.

Characterization and expression analyses of five interferon regulatory factor transcripts (Irf4a, Irf4b, Irf7, Irf8, Irf10) in Atlantic cod (Gadus morhua)

The interferon regulatory factor (IRF) family of genes encodes a group of transcription factors that have important roles not only in regulating the expression of Type I interferons (IFNs) and other genes in the IFN pathway, but also in growth, development and the regulation of oncogenesis. In this study, several IRF family members (Irf4a, Irf4b, Irf7, Irf8, Irf10) in Atlantic cod (Gadus morhua) were characterized at the cDNA and putative amino acid levels, allowing for phylogenetic analysis of these proteins in teleost fish, as well as the development of gene-specific primers used in RT-PCR and quantitative PCR (QPCR) analyses. Two Atlantic cod Irf10 splice variants were identified and their presence confirmed by sequencing of the Irf10 genomic region. RT-PCR showed that Irf7, Irf8 and both Irf10 transcripts were expressed in all 15 cod tissues tested, while Irf4a and Irf4b were absent in some tissues. QPCR analysis of spleen expression expanded upon this, and upon previous work. All IRF transcripts in the study were responsive to stimulation by the viral mimic poly(I:C), and all except Irf4a were responsive to exposure to formalin-killed Aeromonas salmonicida (ASAL). These IRF genes, thus, are likely important in the cod immune response to both viral and bacterial infections. Increased temperature (10 °C to 16 °C) was also observed to modulate the antibacterial responses of all IRF transcripts, and the antiviral responses of Irf4b and Irf10-v2. This research supports earlier studies which reported that elevated temperature modulates the expression of many immune genes in Atlantic cod.

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)

To improve sustainability of aquaculture, especially for carnivorous species like Atlantic cod, replacement of fish oil-based diets with vegetable oil-based diets has been studied. The use of vegetable oil in fish feeds can significantly change the fatty acid composition of fish tissues, and given the importance of fatty acids in inflammation and immunity, this change could potentially impact the immune response and health of the fish. The oilseed Camelina sativa is a promising source for this vegetable oil, because of the high oil content of its seeds (40%), a higher n-3 fatty acid content than most other oilseeds, and a high amount of γ-tocopherol. This study aims to investigate the effect of the replacement of dietary fish oil with oil from Camelina sativa on the immune response of Atlantic cod, as measured by the gene expression in spleen. Juvenile cod were fed on a fish oil-based diet (FO) or one of two diets in which camelina oil replaced 40% or 80% of fish oil (40CO and 80CO respectively) for 67 days, after which they were injected with either the viral mimic polyriboinosinic polyribocytidylic acid (pIC), or phosphate-buffered saline (PBS) as a control. Microarray analysis was used to determine the effect of the diet on the basal spleen transcriptome (pre-injection), and on the response to pIC (24 h post-injection). No marked differences in the spleen transcriptome were found between the three diets, either before or after injection with pIC. All fish, regardless of diet, showed a strong anti-viral response 24 h after pIC injection, with more than 500 genes having a significant difference of expression of 2-fold or higher compared to the PBS-injected fish for the FO, 40CO and 80CO diets. Gene Ontology annotation analysis of the three pIC-responsive gene lists indicated they were highly similar, and that the term 'immune system process' was significantly enriched in the pIC-responsive gene lists for all three diets. QPCR analysis for 5 genes with a known function in the anti-viral innate immune response (LGP2, STAT1, IRF1, ISG15 and viperin) showed modestly (smaller than 2-fold) up-regulated basal expression of LGP2, IRF1 and STAT1 in fish fed 40CO compared to the other diets. After pIC injection, all 5 genes were significantly and strongly up-regulated in pIC-injected fish compared to PBS-injected fish, but no significant differences were found between any of the diets. In conclusion, replacement of up to 80% of fish oil with camelina oil in Atlantic cod diets does not have a strong effect on basal spleen gene expression. Atlantic cod fed on camelina oil-containing diets are capable of mounting a strong anti-viral immune response, which is comparable to that in cod fed with a fish oil diet.

Toxicological effect of single contaminants and contaminant mixtures associated with plant ingredients in novel salmon feeds

Increasing use of plant feed ingredients may introduce contaminants not previously associated with farming of salmonids, such as pesticides and PAHs from environmental sources or from thermal processing of oil seeds. To screen for interaction effects of contaminants newly introduced in salmon feeds, Atlantic salmon primary hepatocytes were used. The xCELLigence cytotoxicity system was used to select optimal dosages of the PAHs benzo(a)pyrene and phenanthrene, the pesticides chlorpyrifos and endosulfan, and combinations of these. NMR and MS metabolic profiling and microarray transcriptomic profiling was used to identify novel biomarkers. Lipidomic and transcriptomic profiling suggested perturbation of lipid metabolism, as well as endocrine disruption. The pesticides gave the strongest responses, despite having less effect on cell viability than the PAHs. Only weak molecular responses were detected in PAH-exposed hepatocytes. Chlorpyrifos suppressed the synthesis of unsaturated fatty acids. Endosulfan affected steroid hormone synthesis, while benzo(a)pyrene disturbed vitamin D3 metabolism. The primary mixture effect was additive, although at high concentrations the pesticides acted in a synergistic fashion to decrease cell viability and down-regulate CYP3A and FABP4 transcription. This work highlights the usefulness of 'omics techniques and multivariate data analysis to investigate interactions within mixtures of contaminants with different modes of action.

An investigation of appetite-related peptide transcript expression in Atlantic cod (Gadus morhua) brain following a Camelina sativa meal-supplemented feeding trial

Camelina sativa is a hardy oilseed crop with seeds that contain high levels of ω3 polyunsaturated fatty acids and protein, which are critical components of fish feed. Camelina might thus be used as a cheaper and more sustainable supplement to fish-based products in aquaculture. Atlantic cod, Gadus morhua, is a species of interest in the aquaculture industry due to a decrease in wild populations and subsequent collapse of some cod fisheries. As cod are carnivorous fish, it is necessary to determine how this species physiologically tolerates plant-based diets. In this study, juvenile Atlantic cod were subjected to 13 weeks of either 15 or 30% camelina meal (CM)-supplemented diets or a control fish meal feed. Growth and food intake were evaluated and the mRNA expression of appetite-related hormones [pro-melanin-concentrating hormone (pmch), hypocretin (synonym: orexin, hcrt), neuropeptide Y (npy) and cocaine- and amphetamine-regulated transcript (cart)] was assessed using quantitative real-time PCR in brain regions related to food intake regulation (telencephalon/preoptic area, optic tectum/thalamus and hypothalamus). CM inclusion diets caused decreases in both growth and food intake in Atlantic cod. Optic tectum pmch transcript expression was significantly higher in fish fed the 30% CM diet compared to fish fed the 15% CM diet. In the hypothalamus, compared to fish fed the control diet, hcrt expression was significantly higher in fish fed the 30% CM diet, while npy transcript expression was significantly higher in fish fed the 15% CM diet. cart mRNA expression was not affected by diet in any brain region. Further studies are needed to determine which factors (e.g. anti-nutritional factors, palatability and nutritional deficits) contribute to reduced feed intake and growth, as well as the maximum CM inclusion level that does not negatively influence feed intake, growth rate and the transcript expression of appetite-related factors in Atlantic cod.

Identification of a gene set to evaluate the potential effects of loud sounds from seismic surveys on the ears of fishes: a study with Salmo salar

Functional genomic studies were carried out on the inner ear of Atlantic salmon Salmo salar following exposure to a seismic airgun. Microarray analyses revealed 79 unique transcripts (passing background threshold), with 42 reproducibly up-regulated and 37 reproducibly down-regulated in exposed v. control fish. Regarding the potential effects on cellular energetics and cellular respiration, altered transcripts included those with roles in oxygen transport, the glycolytic pathway, the Krebs cycle and the electron transport chain. Of these, a number of transcripts encoding haemoglobins that are important in oxygen transport were up-regulated and among the most highly expressed. Up-regulation of transcripts encoding nicotinamide riboside kinase 2, which is also important in energy production and linked to nerve cell damage, points to evidence of neuronal damage in the ear following noise exposure. Transcripts related to protein modification or degradation also indicated potential damaging effects of sound on ear tissues. Notable in this regard were transcripts associated with the proteasome-ubiquitin pathway, which is involved in protein degradation, with the transcript encoding ubiquitin family domain-containing protein 1 displaying the highest response to exposure. The differential expression of transcripts observed for some immune responses could potentially be linked to the rupture of cell membranes. Meanwhile, the altered expression of transcripts for cytoskeletal proteins that contribute to the structural integrity of the inner ear could point to repair or regeneration of ear tissues including auditory hair cells. Regarding potential effects on hormones and vitamins, the protein carrier for thyroxine and retinol (vitamin A), namely transthyretin, was altered at the transcript expression level and it has been suggested from studies in mammalian systems that retinoic acid may play a role in the regeneration of damaged hair cells. The microarray experiment identified the transcript encoding growth hormone I as up-regulated by loud sound, supporting previous evidence linking growth hormone to hair cell regeneration in fishes. Quantitative (q) reverse transcription (RT) polymerase chain reaction (qRT-PCR) analyses confirmed dysregulation of some microarray-identified transcripts and in some cases revealed a high level of biological variability in the exposed group. These results support the potential utility of molecular biomarkers to evaluate the effect of seismic surveys on fishes with studies on the ears being placed in a priority category for development of exposure-response relationships. Knowledge of such relationships is necessary for addressing the question of potential size of injury zones.

Variation in embryonic mortality and maternal transcript expression among Atlantic cod (Gadus morhua) broodstock: a functional genomics study

Early life stage mortality is an important issue for Atlantic cod aquaculture, yet the impact of the cod maternal (egg) transcriptome on egg quality and mortality during embryonic development is poorly understood. In the present work, we studied embryonic mortality and maternal transcript expression using eggs from 15 females. Total mortality at 7days post-fertilization (7 dpf, segmentation stage) was used as an indice of egg quality. A 20,000 probe (20K) microarray experiment compared the 7hours post-fertilization (7 hpf, ~2-cell stage) egg transcriptome of the two lowest quality females (>90% mortality at 7 dpf) to that of the highest quality female (~16% mortality at 7 dpf). Forty-three microarray probes were consistently differentially expressed in both low versus high quality egg comparisons (25 higher expressed in low quality eggs, and 18 higher expressed in high quality eggs). The microarray experiment also identified many immune-relevant genes [e.g. interferon (IFN) pathway genes ifngr1 and ifrd1)] that were highly expressed in eggs of all 3 females regardless of quality. Twelve of the 43 candidate egg quality-associated genes, and ifngr1, ifrd1 and irf7, were included in a qPCR study with 7 hpf eggs from all 15 females. Then, the genes that were confirmed by qPCR to be greater than 2-fold differentially expressed between 7 hpf eggs from the lowest and highest quality females (dcbld1, ddc, and acy3 more highly expressed in the 2 lowest quality females; kpna7 and hacd1 more highly expressed in the highest quality female), and the 3 IFN pathway genes, were included in a second qPCR study with unfertilized eggs. While some maternal transcripts included in these qPCR studies were associated with extremes in egg quality, there was little correlation between egg quality and gene expression when all females were considered. Both dcbld1 and ddc showed greater than 100-fold differences in transcript expression between females and were potentially influenced by family. The Atlantic cod ddc (dopa decarboxylase) complete cDNA was characterized, and has a 1461bp open reading frame encoding a 486 amino acid protein that contains all eight residues of the conserved pyridoxal 5'-phosphate binding site including the catalytic lysine. This study provides valuable new information and resources related to the Atlantic cod egg transcriptome. Some of these microarray-identified, qPCR-confirmed, Atlantic cod egg transcripts (e.g. ddc, kpna7) play important roles during embryonic development of other vertebrate species, and may have similar functions in Atlantic cod.

Family-specific differences in growth rate and hepatic gene expression in juvenile triploid growth hormone (GH) transgenic Atlantic salmon (Salmo salar)

Growth hormone transgenic (GHTg) Atlantic salmon (Salmo salar) have enhanced growth when compared to their non-transgenic counterparts, and this trait can be beneficial for aquaculture production. Biological confinement of GHTg Atlantic salmon may be achieved through the induction of triploidy (3N). The growth rates of triploid GH transgenic (3NGHTg) Atlantic salmon juveniles were found to significantly vary between families in the AquaBounty breeding program. In order to characterize gene expression associated with enhanced growth in juvenile 3NGHTg Atlantic salmon, a functional genomics approach (32K cDNA microarray hybridizations followed by QPCR) was used to identify and validate liver transcripts that were differentially expressed between two fast-growing 3NGHTg Atlantic salmon families (AS11, AS26) and a slow-growing 3NGHTg Atlantic salmon family (AS25); juvenile growth rate was evaluated over a 45-day period. Of 687 microarray-identified differentially expressed features, 143 (116 more highly expressed in fast-growing and 27 more highly expressed in slow-growing juveniles) were identified in the AS11 vs. AS25 microarray study, while 544 (442 more highly expressed in fast-growing and 102 more highly expressed in slow-growing juveniles) were identified in the AS26 vs. AS25 microarray study. Forty microarray features (39 putatively associated with fast growth and 1 putatively associated with slow growth) were present in both microarray experiment gene lists. The expression levels of 15 microarray-identified transcripts were studied using QPCR with individual RNA samples to validate microarray results and to study biological variability of transcript expression. The QPCR results agreed with the microarray results for 12 of 13 putative fast-growth associated transcripts, but QPCR did not validate the microarray results for 2 putative slow-growth associated transcripts. Many of the 39 microarray-identified genes putatively associated at the transcript expression level with fast-growing 3NGHTg salmon juveniles (including APOA1, APOA4, B2M, FADSD6, FTM, and GAPDH) are involved in metabolism, iron homeostasis and oxygen transport, and immune- or stress-related responses. The results of this study increase our knowledge of family-specific impacts on growth rate and hepatic gene expression in juvenile 3NGHTg Atlantic salmon. In addition, this study provides a suite of putative rapid growth rate-associated transcripts that may contribute to the development of molecular markers [e.g. intronic, exonic or regulatory region single nucleotide polymorphisms (SNPs)] for the selection of GHTg Atlantic salmon broodstock that can be utilized to produce sterile triploids of desired growth performance for future commercial applications.

Gene expression and pathologic alterations in juvenile rainbow trout due to chronic dietary TCDD exposure

The goal of this project was to use functional genomic methods to identify molecular biomarkers as indicators of the impact of TCDD exposure in rainbow trout. Specifically, we investigated the effects of chronic dietary TCDD exposure on whole juvenile rainbow trout global gene expression associated with histopathological analysis. Juvenile rainbow trout were fed Biodiet starter with TCDD added at 0, 0.1, 1, 10 and 100 ppb (ngTCDD/g food), and fish were sampled from each group at 7, 14, 28 and 42 days after initiation of feeding. 100 ppb TCDD caused 100% mortality at 39 days. Fish fed with 100 ppb TCDD food had TCDD accumulation of 47.37 ppb (ngTCDD/g fish) in whole fish at 28 days. Histological analysis from TCDD-treated trout sampled from 28 and 42 days revealed that obvious lesions were found in skin, oropharynx, liver, gas bladder, intestine, pancreas, nose and kidney. In addition, TCDD caused anemia in peripheral blood, decreases in abdominal fat, increases of remodeling of fin rays, edema in pericardium and retrobulbar hemorrhage in the 100 ppb TCDD-treated rainbow trout compared to the control group at 28 days. Dose- and time-dependent global gene expression analyses were performed using the cGRASP 16,000 (16K) cDNA microarray. TCDD-responsive whole body transcripts identified in the microarray experiments have putative functions involved in various biological processes including growth, cell proliferation, metabolic process, and immune system processes. Nine microarray-identified genes were selected for QPCR validation. CYP1A3 and CYP1A1 were common up-regulated genes and HBB1 was a common down-regulated gene among each group based on microarray data, and their QPCR validations are consistent with microarray data for the 10 and 100 ppb TCDD treatment groups after 28 days exposure (p<0.05). In addition, in the 100 ppb group at 28 days, expression of complement component C3-1 and trypsin-1 precursor have a more than 10-fold induction from the microarray experiments, and their QPCR validations are consistent and showed significant induction in the 100 ppb group at 28 days (p<0.05). Overall, lesion in nasal epithelium is a novel and significant result in this study, and TCDD-responsive rainbow trout transcripts identified in the present study may lead to the development of new molecular biomarkers for assessing the potential impacts of environmental TCDD on rainbow trout populations.

The impact of a moderate chronic temperature increase on spleen immune-relevant gene transcription depends on whether Atlantic cod (Gadus morhua) are stimulated with bacterial versus viral antigens

Exposure to elevated temperature is an inherent feature of Atlantic cod (Gadus morhua) sea-cage culture in some regions (e.g., Newfoundland) and may also become an increasingly prevalent challenge for wild fish populations because of accelerated climate change. Therefore, understanding how elevated temperatures impacts the immune response of this commercially important species may help to reduce the potential negative impacts of such challenges. Previously, we investigated the impacts of moderately elevated temperature on the antiviral responses of Atlantic cod (Hori et al. 2012) and reported that elevated temperature modulated the spleen transcriptome response to polyriboinosinic polyribocytidylic acid (pIC, a viral mimic). Herein, we report a complementary microarray study that investigated the impact of the same elevated temperature regime on the Atlantic cod spleen transcriptome response to intraperitoneal (IP) injection of formalin-killed Aeromonas salmonicida (ASAL). Fish were held at two different temperatures (10 °C and 16 °C) prior to immune stimulation and sampled 6 and 24 h post-injection (HPI). In this experiment, we identified 711 and 666 nonredundant ASAL-responsive genes at 6HPI and 24HPI, respectively. These included several known antibacterial genes, including hepcidin, cathelicidin, ferritin heavy subunit, and interleukin 8. However, we only identified 15 differentially expressed genes at 6HPI and 2 at 24HPI (FDR 1%) when comparing ASAL-injected fish held at 10 °C versus 16 °C. In contrast, the same comparisons with pIC-injected fish yielded 290 and 339 differentially expressed genes (FDR 1%) at 6HPI and 24HPI, respectively. These results suggest that moderately elevated temperature has a lesser effect on the Atlantic cod spleen transcriptome response to ASAL (i.e., the antibacterial response) than to pIC (i.e., antiviral response). Thus, the impacts of high temperatures on the cod's immune response may be pathogen dependent.

A moderate increase in ambient temperature modulates the Atlantic cod (Gadus morhua) spleen transcriptome response to intraperitoneal viral mimic injection

BACKGROUND

Atlantic cod (Gadus morhua) reared in sea-cages can experience large variations in temperature, and these have been shown to affect their immune function. We used the new 20K Atlantic cod microarray to investigate how a water temperature change which, simulates that seen in Newfoundland during the spring-summer (i.e. from 10°C to 16°C, 1°C increase every 5 days) impacted the cod spleen transcriptome response to the intraperitoneal injection of a viral mimic (polyriboinosinic polyribocytidylic acid, pIC).

RESULTS

The temperature regime alone did not cause any significant increases in plasma cortisol levels and only minor changes in spleen gene transcription. However, it had a considerable impact on the fish spleen transcriptome response to pIC [290 and 339 significantly differentially expressed genes between 16°C and 10°C at 6 and 24 hours post-injection (HPI), respectively]. Seventeen microarray-identified transcripts were selected for QPCR validation based on immune-relevant functional annotations. Fifteen of these transcripts (i.e. 88%), including DHX58, STAT1, IRF7, ISG15, RSAD2 and IκBα, were shown by QPCR to be significantly induced by pIC.

CONCLUSIONS

The temperature increase appeared to accelerate the spleen immune transcriptome response to pIC. We found 41 and 999 genes differentially expressed between fish injected with PBS vs. pIC at 10°C and sampled at 6HPI and 24HPI, respectively. In contrast, there were 656 and 246 genes differentially expressed between fish injected with PBS vs. pIC at 16°C and sampled at 6HPI and 24HPI, respectively. Our results indicate that the modulation of mRNA expression of genes belonging to the NF-κB and type I interferon signal transduction pathways may play a role in controlling temperature-induced changes in the spleen's transcript expression response to pIC. Moreover, interferon effector genes such as ISG15 and RSAD2 were differentially expressed between fish injected with pIC at 10°C vs. 16°C at 6HPI. These results substantially increase our understanding of the genes and molecular pathways involved in the negative impacts of elevated ambient temperature on fish health, and may also be valuable to our understanding of how accelerated global climate change could impact cold-water marine finfish species.

The mRNA expression of cortisol axis related genes differs in Atlantic cod (Gadus morhua) categorized as high or low responders

Cortisol is a major stress hormone in fish and is known, under normal or stressful conditions, to affect several physiological processes including growth and immunity. Thus, efforts have been made for several cultured finfish species, including the Atlantic cod, to determine whether fish with a high or low cortisol response to stress can be identified and selected. However, we have a limited understanding of the mechanisms that determine these two phenotypes. Thus, we measured total and free plasma cortisol levels in high and low responding cod when subjected to a 30 s handling stress, and the mRNA expression of four key genes in the glucocorticoid (i.e. cortisol) stress axis both pre- and post-stress. The cortisol data is consistent with our previous findings for cod, with high responding (HR) fish having ∼3-fold higher total and free plasma cortisol levels when compared to low responding (LR) fish. Three of the transcripts studied encode key proteins involved in steroidogenesis (StAR, P450scc and 3βHSD), and the constitutive mRNA expression of all three genes was significantly higher (∼2-fold) in the head kidney of HR fish when compared to LR cod. The other gene of interest was the glucocorticoid receptor (GR). We partly cloned and characterized a cDNA from Atlantic cod likely to be this fish's ortholog of the teleost GR1, and showed that while there was no difference in hepatic constitutive GR mRNA expression between groups, HR fish had liver GR mRNA levels that were significantly (1.8-fold) higher at 3 h post-stress as compared to LR fish. Our results suggest that the different magnitude of cortisol response between LR and HR fish is at least partially determined by the capacity of the interrenal tissue to produce steroids.

Expression analysis of glycerol synthesis-related liver transcripts in rainbow smelt (Osmerus mordax) exposed to a controlled decrease in temperature

Rainbow smelt (Osmerus mordax) accumulate high glycerol levels to avoid freezing at subzero temperatures. Glyceroneogenesis is activated by low temperature and occurs in liver via a branch in glycolysis and gluconeogenesis. In this study, carbohydrate and liver transcript levels of 21 genes potentially associated with glycerol production were assessed during a controlled warm to cold transition. Smelt were held at 8°C (warm smelt; non-glycerol accumulating) or subjected to a controlled decrease in water temperature from 8° to 0°C (cold smelt; glycerol accumulating) and sampled at the end of the temperature decrease and 1 mo later. In cold smelt compared with warm smelt, liver glycogen levels were lower and phosphoglucomutase transcript levels were higher. Plasma glycerol levels were higher and increased over time in cold smelt; in cold smelt, liver phosphofructokinase and pyruvate dehydrogenase kinase transcript levels increased over time. These findings imply that glycerol production is being fueled by glycogen degradation and inhibition of pyruvate oxidation serves to channel metabolic flux toward glycerol as opposed to complete glycolysis. Plasma glucose and liver glucose-6-phosphatase transcript levels were higher. Lipoprotein lipase transcript levels were higher, suggesting enhanced lipid breakdown to fuel energy metabolism. Glutamine synthetase transcript levels were higher, perhaps to store nitrogen for biosynthesis in spring.

Identification and validation of differentially expressed transcripts in a hepatocyte model of cold-induced glycerol production in rainbow smelt (Osmerus mordax)

Rainbow smelt ( Osmerus mordax ) avoid freezing by producing antifreeze protein (AFP) and accumulating glycerol. Glyceroneogenesis occurs in liver via a branch in glycolysis and gluconeogenesis and is activated by low temperature. Hepatocytes were isolated from the livers of fish acclimated to 8°C. Cells were incubated at warm (8°C; nonglycerol accumulating) or cold (0.4°C; glycerol accumulating) temperature over a 72-h time course. Reciprocal suppression subtractive hybridization libraries enriched for cold-responsive transcripts were constructed at 72 h. Microarray analyses using a 16K salmonid cDNA array were performed at 24, 48, and 72 h. Expression of type II AFP and 21 carbohydrate, amino acid, or lipid metabolism-related transcripts were validated using quantitative RT-PCR. Type II AFP transcript levels were not directly temperature related. In cold cells, levels of the glucose synthesis transcript were transiently higher. Increased glycerol production was not associated with increased phosphofructokinase or cytosolic glycerol-3-phosphate dehydrogenase transcript levels. Levels of transcripts (phosphoenolpyruvate carboxykinase, mitochondrial malate dehydrogenase, alanine aminotransferase, glutamate dehydrogenase, and aquaglyceroporin 9) associated with mobilization of amino acids to fuel glycerol accumulation were all transiently higher, suggesting a common regulatory mechanism. In cold compared with warm cells, pyruvate dehydrogenase kinase [an inhibitor of pyruvate dehydrogenase (PDH)] transcript levels were 20-fold higher. Potent inhibition of PDH would direct pyruvate and oxaloacetate derived from amino acids to glycerol, as opposed to oxidation via the citric acid cycle. Levels of a transcript potentially encoding glycerol-3-phosphatase, an enzyme not yet characterized in any vertebrate species, were higher following cold incubation. Finally, this study also presents the novel finding of increased glutamine synthetase transcript levels in response to low temperature.

Identification and molecular cloning of Atlantic cod (Gadus morhua) activating transcription factor 3 (ATF3) transcript and its induction in spleen following intraperitoneal polyriboinosinic polyribocytidylic acid injection

Activating transcription factor 3 (ATF3) participates in cellular processes to adapt to various extra- and intra-cellular changes including the modulation of immunity to prevent uncontrolled immune responses to pathogens. In teleost fishes, the involvement of ATF3 in immune response has not been documented. In this study, the putative Atlantic cod (Gadus morhua) ATF3 transcript was identified by performing rapid amplification of cDNA ends (RACE) based on unknown expressed sequence tags (ESTs) that are potentially inducible by polyriboinosinic polyribocytidylic acid (pIC, a synthetic double-stranded RNA viral mimic) in Atlantic cod. ATF3-like ESTs were the most abundant unknown transcript (i.e. lacking significant BLAST hits) generated from a previously constructed cDNA library enriched for pIC inducible transcripts in Atlantic cod spleen. The full-length cDNA of cod ATF3 consists of 2329 nucleotides with an open reading frame (ORF) of 735 bp encoding 244 amino acids. The deduced amino acid sequence of Atlantic cod ATF3 shares over 45% identity with its putative orthologs from other vertebrates. In addition, the presence of a conserved basic region leucine zipper (bZIP) domain in the deduced Atlantic cod ATF3-like protein further supports its identity as an ATF3 homolog. In the spleen of Atlantic cod challenged with intraperitoneal (IP) injections of pIC, the time-course transcript expression of ATF3 was studied using quantitative reverse transcription-polymerase chain reaction (QPCR). At 6 h following the pIC injection, the relative expression level of ATF3 mRNA was significantly up-regulated in comparison to a pre-injected control (61.9-fold) and its time-matched saline-injected control (97.3-fold). At 24 h following the pIC injection, the mRNA expression level of cod ATF3 had subsided and was no longer significantly different from its pre-injected control, but significantly higher (1.88-fold) than its time-matched saline-injected control. Collectively, these results suggest that ATF3 may be involved in the modulation of innate anti-viral response in Atlantic cod.

Development and experimental validation of a 20K Atlantic cod (Gadus morhua) oligonucleotide microarray based on a collection of over 150,000 ESTs

The collapse of Atlantic cod (Gadus morhua) wild populations strongly impacted the Atlantic cod fishery and led to the development of cod aquaculture. In order to improve aquaculture and broodstock quality, we need to gain knowledge of genes and pathways involved in Atlantic cod responses to pathogens and other stressors. The Atlantic Cod Genomics and Broodstock Development Project has generated over 150,000 expressed sequence tags from 42 cDNA libraries representing various tissues, developmental stages, and stimuli. We used this resource to develop an Atlantic cod oligonucleotide microarray containing 20,000 unique probes. Selection of sequences from the full range of cDNA libraries enables application of the microarray for a broad spectrum of Atlantic cod functional genomics studies. We included sequences that were highly abundant in suppression subtractive hybridization (SSH) libraries, which were enriched for transcripts responsive to pathogens or other stressors. These sequences represent genes that potentially play an important role in stress and/or immune responses, making the microarray particularly useful for studies of Atlantic cod gene expression responses to immune stimuli and other stressors. To demonstrate its value, we used the microarray to analyze the Atlantic cod spleen response to stimulation with formalin-killed, atypical Aeromonas salmonicida, resulting in a gene expression profile that indicates a strong innate immune response. These results were further validated by quantitative PCR analysis and comparison to results from previous analysis of an SSH library. This study shows that the Atlantic cod 20K oligonucleotide microarray is a valuable new tool for Atlantic cod functional genomics research.

An integrated approach to gene discovery and marker development in Atlantic cod (Gadus morhua)

Atlantic cod is a species that has been overexploited by the capture fishery. Programs to domesticate this species are underway in several countries, including Canada, to provide an alternative route for production. Selective breeding programs have been successfully applied in the domestication of other species, with genomics-based approaches used to augment conventional methods of animal production in recent years. Genomics tools, such as gene sequences and sets of variable markers, also have the potential to enhance and accelerate selective breeding programs in aquaculture, and to provide better monitoring tools to ensure that wild cod populations are well managed. We describe the generation of significant genomics resources for Atlantic cod through an integrated genomics/selective breeding approach. These include 158,877 expressed sequence tags (ESTs), a set of annotated putative transcripts and several thousand single nucleotide polymorphism markers that were developed from, and have been shown to be highly variable in, fish enrolled in two selective breeding programs. Our EST collection was generated from various tissues and life cycle stages. In some cases, tissues from which libraries were generated were isolated from fish exposed to stressors, including elevated temperature, or antigen stimulation (bacterial and viral) to enrich for transcripts that are involved in these response pathways. The genomics resources described here support the developing aquaculture industry, enabling the application of molecular markers within selective breeding programs. Marker sets should also find widespread application in fisheries management.

Characterization and expression studies of Gaduscidin-1 and Gaduscidin-2; paralogous antimicrobial peptide-like transcripts from Atlantic cod (Gadus morhua)

Piscidins are a family of antimicrobial peptides (AMPs) from fish that constitute an important component of their innate immune system. Based on previously generated Atlantic cod (Gadus morhua) expressed sequence tags (ESTs), we identified sequences representing two paralogous AMP-like transcripts. These Atlantic cod paralogues were termed gaduscidins (GAD-1 and GAD-2), derived from the genus name Gadus. We obtained full-length cDNA sequences for these putative AMP-encoding transcripts using bi-directional rapid amplification of cDNA ends (RACE). GAD-1 and GAD-2 putative peptides exhibit sequence similarity with members of the piscidin family from teleost fish. Quantitative reverse transcription - polymerase chain reaction (QPCR) was utilized in transcript expression studies of GAD-1 and GAD-2. We examined the constitutive expression of these transcripts in six tissues (head kidney, blood, brain, gill, pyloric caecum, and spleen) of non-stressed juvenile cod; GAD-1 and GAD-2 transcripts were detected in all six tissues, with the highest expression of both transcripts being in spleen, head kidney, and gill. Transcript expression of GAD-1 and GAD-2 was also examined in immune tissues (spleen and head kidney) following intraperitoneal (IP) injection of formalin-killed, atypical Aeromonas salmonicida (ASAL) or phosphate-buffered saline (PBS control). Both transcripts were weakly (less than 4-fold) up-regulated by ASAL in spleen but non-responsive to ASAL in head kidney. Since GAD-1 and GAD-2 transcripts are highly expressed constitutively in immune-relevant tissues (e.g. spleen and head kidney), and are weakly induced in spleen following IP injection with bacterial antigens, they may represent important components of innate immunity in Atlantic cod.

The fight between the teleost fish immune response and aquatic viruses

Teleost fish represent a transition point on the phylogenetic spectrum between invertebrates that depend only on innate immunity and mammals that heavily depend on adaptive immunity. The major mechanisms of the teleost fish innate immune response are suggested to be similar to mammals, although fine details of the process require further studies. Within the innate immune response the type I interferon (IFN) system is an essential innate antiviral component that protects fish from some virus infections. The current progress of cloning and functional characterization of fish antiviral genes is promising in further elucidation of the fish antiviral response. The adaptive immune system of fish utilizes cellular components more or less similar to mammals. Teleost fish produce IgM as a primary antibody response and lack isotype switching to mount virus-specific antibodies during the infection process. Despite this, the development of successful fish rhabdoviral vaccines suggest that vaccination may prove to be an effective way of promoting fish adaptive immune responses to viruses. This paper reviews the bony fish antiviral response with specific discussion on the evolutionary mechanisms that allow aquatic viruses to co-exist with their host. Detailed aspects of the teleost type I IFN system are also addressed.

Transcriptome responses to heat stress in the nucleated red blood cells of the rainbow trout (Oncorhynchus mykiss)

The retention of a nucleus in the mature state of fish red blood cells (RBCs) and the ability to easily collect and manipulate blood in nonterminal experiments make blood an ideal tissue on which to study the cellular stress response in fish. Through the use of the cGRASP 16K salmonid microarray, we investigated differences in RBC global gene transcription in fish held under control conditions (11°C) and exposed to heat stress (1 h at 25°C followed by recovery at 11°C). Repeated blood sampling (via a dorsal aorta cannula) enables us to examine the individual stress response over time. Samples were taken preheat stress (representing individual control) and at 4 and 24 h postheat stress (representing early and late transcriptional regulation). Approximately 3,000 microarray features had signal above threshold when hybridized with RBC RNA-derived targets, and cannulation did not have a detectable effect on RBC mRNA expression at the investigated time points. Genes involved in the stress response, immune response, and apoptosis were among those showing the highest dysregulation during both early and late transcriptional regulation. Additionally, genes related to the differentiation and development of blood cells were transcriptionally upregulated at the 24 h time point. This study provides a broader understanding of the mechanisms underpinning the stress response in fish and the discovery of novel genes that are regulated in a stress specific manner. Moreover, salmonid transcripts that are consistently dysregulated in blood in response to heat stress are potential candidates of nonlethal biomarkers of exposure to this particular stressor.

Atlantic cod (Gadus morhua) CC chemokines: Diversity and expression analysis

Chemokines are a large, diverse group of small cytokines that can be classified into several families, including the CC chemokines that are characterized by two adjacent cysteines near their amino terminus. CC chemokines play a pivotal role in host defense mechanisms by inducing leukocyte chemotaxis under physiological and inflammatory conditions. Analysis of CC chemokines from teleost fishes indicates that the number of CC chemokine genes and their tissue expression patterns vary largely in this group of vertebrates. Here we describe 32 distinct CC chemokine sequences from Atlantic cod (Gadus morhua) identified by analysis of approximately 206,000 ESTs. Phylogenetic analysis of Atlantic cod CC chemokines placed these sequences in seven clusters, most likely resulting from species-specific gene duplications, and two unique sequences; 12 of these CC chemokines, including at least one member of each cluster, were analyzed by QPCR using four immune-related tissues (head kidney, liver, spleen and blood) obtained from unstimulated, polyriboinosinic polyribocytidylic acid (pIC)-stimulated and formalin-killed atypical Aeromonas salmonicida-stimulated individuals. EST abundance and QPCR analysis indicate that the expression of closely related CC chemokines GmSCYA101 and GmSCYA102, GmSCYA108 and GmSCYA109 or GmSCYA122 and GmSCYA124 can be highly tissue-specific despite substantial sequence identity. Stimulation with the viral mimic pIC or formalin-killed atypical A. salmonicida resulted in increased expression of most of the CC chemokines, indicating that they can be regarded as either inducible (inflammatory) or dual-function rather than constitutive (homeostatic). Tissue specificity, and the level of induction, varied broadly; for example, GmSCYA123 was at least 4-fold up-regulated by both inducers in all tissues analyzed, whereas pIC increased the expression of GmSCYA124 in liver over 1500 times.

Impact of asymptomatic nodavirus carrier state and intraperitoneal viral mimic injection on brain transcript expression in Atlantic cod (Gadus morhua)

Nodaviruses and other RNA viruses have a profoundly negative impact on the global aquaculture industry. Nodaviruses target nervous tissue causing viral nervous necrosis, a disease characterized by neurological damage, swimming abnormalities, and morbidity. This study used functional genomic techniques to study the Atlantic cod (Gadus morhua) brain transcript expression responses to asymptomatic high nodavirus carrier state and intraperitoneal injection of polyriboinosinic polyribocytidylic acid (pIC). Reciprocal suppression subtractive hybridization (SSH) cDNA libraries enriched for virus-responsive brain transcripts were constructed and characterized. We generated 1,938 expressed sequence tags (ESTs) from a forward brain SSH library (enriched for transcripts upregulated by nodavirus and/or pIC) and 1,980 ESTs from a reverse brain SSH library (enriched for transcripts downregulated by nodavirus and/or pIC). To examine the effect of nodavirus carrier state on individual brain gene expression in asymptomatic cod, 27 transcripts of interest were selected for quantitative reverse transcription-polymerase chain reaction (QPCR) studies. Transcripts found to be >10-fold upregulated in individuals with a high nodavirus carrier state relative to those in a no/low nodavirus carrier state were identified as ISG15, IL8, DHX58 (alias LGP2), ZNFX1, RSAD2 (alias viperin), and SACS (sacsin, alias spastic ataxia of Charlevoix-Saguenay). These and other SSH-identified transcripts were also found by QPCR to be significantly (P < 0.05) upregulated by pIC compared with saline-injected controls within 72 h of injection. Several transcripts identified in the reverse SSH library, including two putative ubiquitination pathway members (HERC4 and SUMO2), were found to be significantly (P < 0.05) downregulated in individuals with a high nodavirus carrier state. Our data shows that Atlantic cod brains have a strong interferon pathway response to asymptomatic high nodavirus carrier state and that many interferon pathway and other immune relevant transcripts are significantly induced in brain by both nodavirus and pIC.

Heat-shock responsive genes identified and validated in Atlantic cod (Gadus morhua) liver, head kidney and skeletal muscle using genomic techniques

BACKGROUND

Daily and seasonal changes in temperature are challenges that fish within aquaculture settings cannot completely avoid, and are known to elicit complex organismal and cellular stress responses. We conducted a large-scale gene discovery and transcript expression study in order to better understand the genes that are potentially involved in the physiological and cellular aspects of stress caused by heat-shock. We used suppression subtractive hybridization (SSH) cDNA library construction and characterization to identify transcripts that were dysregulated by heat-shock in liver, skeletal muscle and head kidney of Atlantic cod. These tissues were selected due to their roles in metabolic regulation, locomotion and growth, and immune function, respectively. Fish were exposed for 3 hours to an 8 degrees C elevation in temperature, and then allowed to recover for 24 hours at the original temperature (i.e. 10 degrees C). Tissue samples obtained before heat-shock (BHS), at the cessation of heat-shock (CS), and 3, 12, and 24 hours after the cessation of heat-shock (ACS), were used for reciprocal SSH library construction and quantitative reverse transcription - polymerase chain reaction (QPCR) analysis of gene expression using samples from a group that was transferred but not heat-shocked (CT) as controls.

RESULTS

We sequenced and characterized 4394 ESTs (1524 from liver, 1451 from head kidney and 1419 from skeletal muscle) from three "forward subtracted" libraries (enriched for genes up-regulated by heat-shock) and 1586 from the liver "reverse subtracted" library (enriched for genes down-regulated by heat-shock), for a total of 5980 ESTs. Several cDNAs encoding putative chaperones belonging to the heat-shock protein (HSP) family were found in these libraries, and "protein folding" was among the gene ontology (GO) terms with the highest proportion in the libraries. QPCR analysis of HSP90alpha and HSP70-1 (synonym: HSPA1A) mRNA expression showed significant up-regulation in all three tissues studied. These transcripts were more than 100-fold up-regulated in liver following heat-shock. We also identified HSP47, GRP78 and GRP94-like transcripts, which were significantly up-regulated in all 3 tissues studied. Toll-like receptor 22 (TLR22) transcript, found in the liver reverse SSH library, was shown by QPCR to be significantly down-regulated in the head kidney after heat-shock.

CONCLUSION

Chaperones are an important part of the cellular response to stress, and genes identified in this work may play important roles in resistance to thermal-stress. Moreover, the transcript for one key immune response gene (TLR22) was down-regulated by heat-shock, and this down-regulation may be a component of heat-induced immunosuppression.

Stimulation of growth and changes in the hepatic transcriptome by 17β-estradiol in the yellow perch (Perca flavescens)

The effects of dietary 17β-estradiol (E2) on growth and liver transcriptomics were investigated in the yellow perch ( Perca flavescens). After a 3-mo treatment, E2 significantly stimulated an increase in length and weight of juvenile male and female perch relative to control animals. The increase was significantly greater in females compared with males. Separate, unnormalized cDNA libraries were constructed from equal quantities of RNA from 6 male and 6 female livers of E2-treated and control perch, and 3,546 and 3,719 expressed sequence tags (ESTs) were obtained, respectively. To characterize E2-regulated transcripts, EST frequencies between libraries were calculated within contiguous sequences that were assembled from the combined ESTs of both libraries. Frequencies were also determined in EST transcript groupings produced by aligning all of the ESTs from both libraries at the nucleotide level. From these analyses, there were 28 annotated transcripts that were regulated by 75% between libraries and for which there were at least 5 ESTs of the same transcript between libraries. Regulation of a subset ( 14 ) of these transcripts was confirmed by quantitative reverse transcription-polymerase chain reaction (QPCR). Transcripts that were upregulated by E2 included reproduction-related proteins, binding proteins, and proteases and protease inhibitors. While not part of the transcript frequency analysis, QPCR showed significant upregulation of estrogen receptor esr1 and of insulin-like growth factor I (IGF-I) in E2 livers. E2-downregulated transcripts represented a variety of functional categories including components of the respiratory chain, lipid transport and metabolism, glycolysis, amino acid and nitrogen metabolism, binding proteins, a hydrolytic enzyme, and a transcriptional regulator. In perch it appears that exogenous estrogen drastically shifts liver metabolism toward the production of lipoproteins and carbohydrate binding proteins, and that the growth-promoting action may involve an increase in hepatic IGF-I production.

Infectious salmon anaemia virus (ISAV) isolates induce distinct gene expression responses in the Atlantic salmon (Salmo salar) macrophage/dendritic-like cell line TO, assessed using genomic techniques

Infectious salmon anaemia virus (ISAV) is a marine orthomyxovirus of significant interest not only as a cause of a fatal disease of farmed Atlantic salmon resulting in severe economic losses to the aquaculture industry, but also as the only poikilothermic orthomyxovirus. ISAV targets vascular endothelial cells and macrophages, and is known to influence the expression of both innate and adaptive immune response relevant genes. ISAV isolates from different geographic regions have been shown to vary considerably in their pathogenicity for Atlantic salmon. This study aimed to characterize the Atlantic salmon TO macrophage/dendritic-like cell responses to infection with a selection of ISAV isolates of different genotypes and pathogenicity phenotypes. The first TO infection trial used ISAV isolates NBISA01 and RPC/NB-04-085-1 of high and low pathogenicity, respectively, and global gene expression analyses were carried out using approximately 16,000 gene (16K) salmonid cDNA microarrays to compare RNA samples extracted from TO cells harvested 24 and 72h post-infection versus time-matched uninfected controls. Overall, the microarray experiment showed that RPC/NB-04-085-1-infected cells had a higher total number of reproducibly dysregulated genes (88 genes: the sum of genes greater than 2-fold up- or down-regulated in all four replicate microarrays of a given comparison) than the NBISA01-infected cells (10 genes) for the combined sampling points (i.e. 24 and 72h). This microarray experiment identified several salmon genes that were differentially regulated by NBISA01 and RPC/NB-04-085-1, and which may be useful as molecular biomarkers of ISAV infection. An initial quantitative reverse transcription-polymerase chain reaction (QRT-PCR) study involving 25 microarray-identified genes confirmed the differences in the level of dysregulation of host transcripts between the two ISAV isolates (i.e. NBISA01 and RPC/NB-04-085-1). A second TO infection trial was run using a selection of four clinical ISAV isolates (Norway-810/9/99, a high pathogenicity isolate of European genotype; RPC/NB-04-085-1, a low pathogenicity isolate of European genotype; NBISA01, a high pathogenicity isolate of North American genotype; and RPC/NB-01-0593-1, an intermediate pathogenicity isolate of North American genotype), and UV-inactivated RPC/NB-04-085-1, with sampling at 24, 36, 48, 72, 96, and 120h post-infection. The microarray-identified, QRT-PCR validated suite of 24 molecular biomarkers of response to ISAV were used in a second QRT-PCR experiment to assess the TO cell gene expression responses to the four ISAV isolates at all six time points in the infection. The QRT-PCR data showed that RPC/NB-04-085-1 caused the highest fold changes of most immune-relevant genes [such as interferon-inducible protein Gig1, Mx1 protein, interferon-induced protein with tetratricopeptide repeats 5, Radical S-adenosyl methionine domain-containing protein (viperin), and several genes involved in the ISGylation pathway], followed by Norway-810/9/99. NBISA01 and RPC/NB-01-0593-01 (both of North American genotype) showed low fold up-regulation of transcripts that were highly induced by RPC/NB-04-085-1 isolate. These findings show that ISAV isolates have strain-specific variations in their ability to induce immune response genes.

RNA viruses in the sea

Viruses are ubiquitous in the sea and appear to outnumber all other forms of marine life by at least an order of magnitude. Through selective infection, viruses influence nutrient cycling, community structure, and evolution in the ocean. Over the past 20 years we have learned a great deal about the diversity and ecology of the viruses that constitute the marine virioplankton, but until recently the emphasis has been on DNA viruses. Along with expanding knowledge about RNA viruses that infect important marine animals, recent isolations of RNA viruses that infect single-celled eukaryotes and molecular analyses of the RNA virioplankton have revealed that marine RNA viruses are novel, widespread, and genetically diverse. Discoveries in marine RNA virology are broadening our understanding of the biology, ecology, and evolution of viruses, and the epidemiology of viral diseases, but there is still much that we need to learn about the ecology and diversity of RNA viruses before we can fully appreciate their contributions to the dynamics of marine ecosystems. As a step toward making sense of how RNA viruses contribute to the extraordinary viral diversity in the sea, we summarize in this review what is currently known about RNA viruses that infect marine organisms.

Domestication and growth hormone transgenesis cause similar changes in gene expression in coho salmon (Oncorhynchus kisutch)

Domestication has been extensively used in agricultural animals to modify phenotypes such as growth rate. More recently, transgenesis of growth factor genes [primarily growth hormone (GH)] has also been explored as a rapid approach to accelerating performance of agricultural species. Growth rates of many fishes respond dramatically to GH gene transgenesis, whereas genetic engineering of domestic mammalian livestock has resulted in relatively modest gains. The most dramatic effects of GH transgenesis in fish have been seen in relatively wild strains that have undergone little or no selection for enhanced growth, whereas genetic modification of livestock necessarily has been performed in highly domesticated strains that already possess very rapid growth. Such fast-growing domesticates may be refractory to further stimulation if the same regulatory pathways are being exploited by both genetic approaches. By directly comparing gene expression in wild-type, domestic, and GH transgenic strains of coho salmon, we have found that domestication and GH transgenesis are modifying similar genetic pathways. Genes in many different physiological pathways show modified expression in domestic and GH transgenic strains relative to wild-type, but effects are strongly correlated. Genes specifically involved in growth regulation (IGF1, GHR, IGF-II, THR) are also concordantly regulated in domestic and transgenic fish, and both strains show elevated levels of circulating IGF1. Muscle expression of GH in nontransgenic strains was found to be elevated in domesticated fish relative to wild type, providing a possible mechanism for growth enhancement. These data have implications for genetic improvement of existing domesticated species and risk assessment and regulation of emerging transgenic strains.

Identification and analysis of differentially expressed genes in immune tissues of Atlantic cod stimulated with formalin-killed, atypical Aeromonas salmonicida

Physiological changes, elicited in animal immune tissues by exposure to pathogens, may be studied using functional genomics approaches. We created and characterized reciprocal suppression subtractive hybridization (SSH) cDNA libraries to identify differentially expressed genes in spleen and head kidney tissues of Atlantic cod (Gadus morhua) challenged with intraperitoneal injections of formalin-killed, atypical Aeromonas salmonicida. Of 4,154 ESTs from four cDNA libraries, 10 genes with immune-relevant functional annotations were selected for QPCR studies using individual fish templates to assess biological variability. Genes confirmed by QPCR as upregulated by A. salmonicida included interleukin-1 beta, interleukin-8, a small inducible cytokine, interferon regulatory factor 1 (IRF1), ferritin heavy subunit, cathelicidin, and hepcidin. This study is the first large-scale discovery of bacteria-responsive genes in cod and the first to demonstrate upregulation of IRF1 in fish immune tissues as a result of bacterial antigen stimulation. Given the importance of IRF1 in vertebrate immune responses to viral and bacterial pathogens, the full-length cDNA sequence of Atlantic cod IRF1 was obtained and compared with putative orthologous sequences from other organisms. Functional annotations of assembled SSH library ESTs showed that bacterial antigen stimulation caused changes in many biological processes including chemotaxis, regulation of apoptosis, antimicrobial peptide production, and iron homeostasis. Moreover, differences in spleen and head kidney gene expression responses to the bacterial antigens pointed to a potential role for the cod spleen in blood-borne pathogen clearance. Our data show that Atlantic cod immune tissue responses to bacterial antigens are similar to those seen in other fish species and higher vertebrates.

The immune and stress responses of Atlantic cod to long-term increases in water temperature

Sea-caged cod are limited in their movements in the water column, and thus can be exposed to large seasonal ( approximately 0-20 degrees C) temperature fluctuations. To investigate the physiological response of Atlantic cod to summer-like increases in temperature, we exposed 10 degrees C acclimated juvenile cod to a graded thermal challenge (1 degrees C increase every 5 days) and measured: (1) plasma cortisol and glucose levels; (2) the respiratory burst activity of blood leukocytes; and (3) the expression of specific immune-related genes [MHC Class I, Interleukin-1beta (IL-1beta), beta2-microglobulin (beta2-M), Immunoglobulin M (IgM)-light (L) and -heavy (H) chains] in the blood using quantitative reverse transcription-polymerase chain reaction (QRT-PCR). The experiment was stopped at 19.1 degrees C, with 26.7% of the fish surviving to this point. Plasma glucose levels increased slightly at 16 and 18 degrees C (by 1.39- and 1.74-fold, respectively), in contrast, cortisol levels were elevated significantly (by 2.9-fold) at 16 degrees C but returned to control levels thereafter. The effect of increasing temperature on the expression of immune related genes in blood cells (leukocytes) was variable and depended on the gene of interest. The expression of IgM-H remained stable for the duration of the experiment. In contrast, IL-1beta expression was increased significantly (by approximately 25-fold) at 19 degrees C as compared to time-matched control fish, and changes in the expression of beta2-M, MHC Class I and IgM-L followed a pattern similar to that seen for cortisol: increasing at 16 degrees C (by 4.2-, 5.3- and 17-fold, respectively), but returning to pre-stress levels by 19 degrees C. Interestingly, increasing temperatures had no effect on respiratory burst activity. This study is the first to examine the effects of a chronic regimen of increasing temperature on the stress physiology and immunology of a marine teleost, and suggests that immune function is influenced by complex interactions between thermal effects and temperature-induced stress (elevated circulating cortisol levels).