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

BACKGROUND: Antibody affinity maturation in vertebrates requires the enzyme activation-induced cytidine deaminase (AID) which initiates secondary antibody diversification by mutating the immunoglobulin loci. AID-driven antibody diversification is conserved across jawed vertebrates since bony and cartilaginous fish. Two exceptions have recently been reported, the Pipefish and Anglerfish, in which the AID-encoding aicda gene has been lost. Both cases are associated with unusual reproductive behavior, including male pregnancy and sexual parasitism. Several cold water fish in the Atlantic cod (Gadinae) family carry an aicda gene that encodes for a full-length enzyme but lack affinity-matured antibodies and rely on antibodies of broad antigenic specificity. Hence, we examined the functionality of their AID. RESULTS: By combining genomics, transcriptomics, immune responsiveness, and functional enzymology of AID from 36 extant species, we demonstrate that AID of that Atlantic cod and related fish have extremely lethargic or no catalytic activity. Through ancestral reconstruction and functional enzymology of 71 AID enzymes, we show that this enzymatic inactivation likely took place relatively recently at the emergence of the true cod family (Gadidae) from their ancestral Gadiformes order. We show that this AID inactivation is not only concordant with the previously shown loss of key adaptive immune genes and expansion of innate and cell-based immune genes in the Gadiformes but is further reflected in the genomes of these fish in the form of loss of AID-favored sequence motifs in their immunoglobulin variable region genes. CONCLUSIONS: Recent demonstrations of the loss of the aicda gene in two fish species challenge the paradigm that AID-driven secondary antibody diversification is absolutely conserved in jawed vertebrates. These species have unusual reproductive behaviors forming an evolutionary pressure for a certain loss of immunity to avoid tissue rejection. We report here an instance of catalytic inactivation and functional loss of AID rather than gene loss in a conventionally reproducing vertebrate. Our data suggest that an expanded innate immunity, in addition to lower pathogenic pressures in a cold environment relieved the pressure to maintain robust secondary antibody diversification. We suggest that in this unique scenario, the AID-mediated collateral genome-wide damage would form an evolutionary pressure to lose AID function.

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

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

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.

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.

A chromosome-scale Mytilus edulis genome assembly for aquaculture, marine ecology, and evolution.

The smooth-shelled blue mussel, Mytilus edulis is part of the Mytilus species complex, encompassing at least three putative species: M. edulis, Mytilus galloprovincialis, and Mytilus trossulus. These three species occur on both sides of the Atlantic and hybridize in nature, and both M. edulis and M. galloprovincialis are important aquaculture species. They are also invasive species in many parts of the world. Here, we present a chromosome-level assembly of M. edulis. We used a combination of PacBio sequencing and Dovetail's Omni-C technology to generate an assembly with 14 long scaffolds containing 94% of the predicted length of the M. edulis genome (1.6 out of 1.7 Gb). Assembly statistics were as follows: total length = 1.65 Gb, N50 = 116 Mb, L50 = 7, and L90 = 13. BUSCO analysis showed 92.55% eukaryote BUSCOs identified. AB-Initio annotation using RNA-seq from mantle, gills, muscle, and foot predicted 47,128 genes. These gene models were combined with IsoSeq validation resulting in 45,379 full CDS protein sequences and 129,708 isoforms. Using GBS and shotgun sequencing, we also sequenced several eastern Canadian populations of Mytilus to characterize single-nucleotide as well as structural variance. This high-quality genome for M. edulis provides a platform to develop tools that can be used in breeding, molecular ecology and evolution to address questions of both commercial and environmental perspectives.

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.

Rejection of Lepeophtheirus salmonis driven in part by chitin sensing is not impacted by seawater acclimitization in Coho salmon (Oncorhynchus kisutch).

There is tremendous variation in life-history strategies among anadromous salmonids. Species that enter the ocean environment at small sizes (< 20 g) are likely under more physiological pressure from pathogens; however, little data is available on responses at these early stages. With this in mind, we performed salmon louse challenges with Coho salmon either immediately after seawater entry (SW; ca. 10 g) or after 30 days in SW (ca. 20 g). Irrespective of size or time in SW, parasites were rapidly rejected by the host, with > 90% of all parasites lost by 16 days post-infection (dpi). Rejection was concomitant with host epithelial granulomatous infiltrations that initially targeted the embedded frontal filament (4 dpi) and the entire parasite by 10 dpi. Illumina sequencing, followed by functional enrichment analysis, revealed a concerted defense response in the fin within 1 dpi that included multiple innate and adaptive immunity components. Strikingly, early indications of an allergic-type inflammatory response were associated with chitin sensing pathways orchestrated by early overexpression of the IgE-receptor, fcer1g. Additionally, there was profound overexpression of several classes of c-type lectin receptors, including dectin-2, mincle, and dc-sign at 1 dpi onward. These profiles and upregulation of cellular effector markers were corroborated by histopathological evaluation, revealing the simultaneous presence of mast cell/eosinophilic granular cells, sacciform cells, macrophages/histiocytes, and granulocytes in fin. At 10 dpi and concurrent with parasite expulsion, there was evidence of immunoregulation in addition to tissue remodelling pathways. At 16 dpi, the response was effectively abrogated. Simultaneous profiling of the parasite transcriptome revealed early induction of chitin metabolism and immunomodulation, toxin production and ECM degradation; however, after 7 dpi, these were replaced with overexpression of stress and immune defense genes. These data present the first evidence for Coho salmon demonstrating chitin- and sugar moiety-sensing as key drivers of salmon louse rejection.

Unveiling the Impact of Gene Presence/Absence Variation in Driving Inter-Individual Sequence Diversity within the CRP-I Gene Family in Mytilus spp.

Mussels (Mytilus spp.) tolerate infections much better than other species living in the same marine coastal environment thanks to a highly efficient innate immune system, which exploits a remarkable diversification of effector molecules involved in mucosal and humoral responses. Among these, antimicrobial peptides (AMPs) are subjected to massive gene presence/absence variation (PAV), endowing each individual with a potentially unique repertoire of defense molecules. The unavailability of a chromosome-scale assembly has so far prevented a comprehensive evaluation of the genomic arrangement of AMP-encoding loci, preventing an accurate ascertainment of the orthology/paralogy relationships among sequence variants. Here, we characterized the CRP-I gene cluster in the blue mussel Mytilus edulis, which includes about 50 paralogous genes and pseudogenes, mostly packed in a small genomic region within chromosome 5. We further reported the occurrence of widespread PAV within this family in the Mytilus species complex and provided evidence that CRP-I peptides likely adopt a knottin fold. We functionally characterized the synthetic peptide sCRP-I H1, assessing the presence of biological activities consistent with other knottins, revealing that mussel CRP-I peptides are unlikely to act as antimicrobial agents or protease inhibitors, even though they may be used as defense molecules against infections from eukaryotic parasites.

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.

RNA-Seq Analysis of the Growth Hormone Transgenic Female Triploid Atlantic Salmon (Salmo salar) Hepatic Transcriptome Reveals Broad Temperature-Mediated Effects on Metabolism and Other Biological Processes.

This study examined the impact of rearing temperature (10.5, 13.5 or 16.5°C) on the hepatic transcriptome of AquAdvantage Salmon (growth hormone transgenic female triploid Atlantic salmon) at an average weight of 800 g. Six stranded PE libraries were Illumina-sequenced from each temperature group, resulting in an average of over 100 M raw reads per individual fish. RNA-sequencing (RNA-seq) results showed the greatest difference in the number of differentially expressed transcripts (1750 DETs), as revealed by both DESeq2 and edgeR (q < 0.05; fold-change > |1.5|), was between the 10.5 and 16.5°C temperature groups. In contrast, 172 and 52 DETs were found in the 10.5 vs. 13.5°C and the 13.5 vs. 16.5°C comparisons, respectively. Considering the DETs between the 10.5 and 16.5°C groups, 282 enriched gene ontology (GO) terms were identified (q < 0.05), including "response to stress", "immune system process", "lipid metabolic process", "oxidation-reduction process", and "cholesterol metabolic process", suggesting elevated temperature elicited broad effects on multiple biological systems. Pathway analysis using ClueGO showed additional impacts on amino acid and lipid metabolism. There was a significant positive correlation between RNA-seq and real-time quantitative polymerase chain reaction (RT-qPCR) results for 8 of 9 metabolic-related transcripts tested. RT-qPCR results also correlated to changes in fillet tissue composition previously reported in these salmon (e.g., methionine and lysine concentrations positively correlated with hsp90ab1 transcript expression), suggesting that rearing temperature played a significant role in mediating metabolic/biosynthetic pathways of AquAdvantage Salmon. Many transcripts related to lipid/fatty acid metabolism (e.g., elovl2, fabpi, hacd2, mgll, s27a2, thrsp) were downregulated at 16.5°C compared to both other temperature groups. Additionally, enrichment of stress-, apoptosis- and catabolism-relevant GO terms at 16.5°C suggests that this temperature may not be ideal for commercial production when using freshwater recirculating aquaculture systems (RAS). This study relates phenotypic responses to transcript-specific findings and therefore aids in the determination of an optimal rearing temperature for AquAdvantage Salmon. With approval to grow and sell AquAdvantage Salmon in the United States and Canada, the novel insights provided by this research can help industry expansion by promoting optimal physiological performance and health.

Transcriptome Analysis of Atlantic Salmon (Salmo salar) Skin in Response to Sea Lice and Infectious Salmon Anemia Virus Co-Infection Under Different Experimental Functional Diets.

Sea lice (Lepeophtheirus salmonis) are ectoparasitic copepods that cause significant economic loss in marine salmoniculture. In commercial salmon farms, infestation with sea lice can enhance susceptibility to other significant pathogens, such as the highly contagious infectious salmon anemia virus (ISAv). In this study, transcriptomic analysis was used to evaluate the impact of four experimental functional feeds (i.e. 0.3% EPA/DHA+high-ω6, 0.3% EPA/DHA+high-ω6+immunostimulant (IS), 1% EPA/DHA+high-ω6, and 1% EPA/DHA+high-ω3) on Atlantic salmon (Salmo salar) during a single infection with sea lice (L. salmonis) and a co-infection with sea lice and ISAv. The overall objectives were to compare the transcriptomic profiles of skin between lice infection alone with co-infection groups and assess differences in gene expression response among animals with different experimental diets. Atlantic salmon smolts were challenged with L. salmonis following a 28-day feeding trial. Fish were then challenged with ISAv at 18 days post-sea lice infection (dpi), and maintained on individual diets, to establish a co-infection model. Skin tissues sampled at 33 dpi were subjected to RNA-seq analysis. The co-infection's overall survival rates were between 37%-50%, while no mortality was observed in the single infection with lice. With regard to the infection status, 756 and 1303 consensus differentially expressed genes (DEGs) among the four diets were identified in "lice infection vs. pre-infection" and "co-infection vs. pre-infection" groups, respectively, that were shared between the four experimental diets. The co-infection groups (co-infection vs. pre-infection) included up-regulated genes associated with glycolysis, the interferon pathway, complement cascade activity, and heat shock protein family, while the down-regulated genes were related to antigen presentation and processing, T-cell activation, collagen formation, and extracellular matrix. Pathway enrichment analysis conducted between infected groups (lice infection vs. co-infection) resulted in several immune-related significant GO terms and pathways unique to this group, such as "autophagosome", "cytosolic DNA-sensing pathway" and "response to type I interferons". Understanding how experimental functional feeds can impact the host response and the trajectory of co-infections will be an essential step in identifying efficacious intervention strategies that account for the complexities of disease in open cage culture.

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 degrees C) and exposed to heat stress (1 h at 25 degrees C followed by recovery at 11 degrees 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.

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.

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.

Vaccine-Induced Protection Against Furunculosis Involves Pre-emptive Priming of Humoral Immunity in Arctic Charr.

With respect to salmonid aquaculture, one of the most important bacterial pathogens due to high mortality and antibiotic usage is the causative agent of typical furunculosis, Aeromonas salmonicida spp. salmonicida (Asal). In Atlantic salmon, Salmo salar, the host response during infections with Asal is well-documented, with furunculosis outbreaks resulting in significant mortality in commercial settings. However, less is known about the host-pathogen interactions in the emerging aquaculture species, Arctic charr Salvelinus alpinus. Furthermore, there is no data on the efficacy or response of this species after vaccination with commonly administered vaccines against furunculosis. To this end, we examined the immunological response of S. alpinus during infection with Asal, with or without administration of vaccines (Forte Micro®, Forte Micro® + Renogen®, Elanco Animal Health). Artic charr (vaccinated or unvaccinated) were i.p.-injected with a virulent strain of Asal (106 CFUs/mL) and tissues were collected pre-infection/post-vaccination, 8, and 29 days post-infection. Unvaccinated Arctic charr were susceptible to Asal with 72% mortalities observed after 31 days. However, there was 72-82% protection in fish vaccinated with either the single or dual-vaccine, respectively. Protection in vaccinated fish was concordant with significantly higher serum IgM concentrations, and following RNA sequencing and transcriptome assembly, differential expression analysis revealed several patterns and pathways associated with the improved survival of vaccinated fish. Most striking was the dramatically higher basal expression of complement/coagulation factors, acute phase-proteins, and iron hemostasis proteins in pre-challenged, vaccinated fish. Remarkably, following Asal infection, this response was abrogated and instead the transcriptome was characterized by a lack of immune-stimulation compared to that of unvaccinated fish. Furthermore, where pathways of actin assembly and FcγR-mediated phagocytosis were significantly differentially regulated in unvaccinated fish, vaccinated fish showed either the opposite regulation (ForteMicro®), or no impact at all (ForteMicro®Renogen®). The present data indicates that vaccine-induced protection against Asal relies on the pre-activation and immediate control of humoral immune parameters that is coincident with reduced activation of apoptotic (e.g., NF-κB) and actin-associated pathways.

Functional genomic analysis of the response of Atlantic cod (Gadus morhua) spleen to the viral mimic polyriboinosinic polyribocytidylic acid (pIC).

In order to improve our understanding of how Atlantic cod (Gadus morhua) respond to viruses, we characterized immune-related gene expression in spleen tissues following stimulation with a synthetic double-stranded RNA polyriboinosinic polyribocytidylic acid (pIC). We used reciprocal suppression subtractive hybridization (SSH) cDNA libraries and quantitative RTPCR (QPCR) to identify and quantify pIC-responsive transcripts. A total of 3874 expressed sequence tags (ESTs) were generated from SSH libraries enriched for genes responsive to pIC. Thirteen immune-relevant genes from the libraries were subjected to QPCR. Genes confirmed as up-regulated by pIC included interferon stimulated gene 15, a small inducible cytokine, interferon regulatory factors (1, 7, and 10), MHC class I, viperin, and ATP-dependent helicase LGP2. Alpha-1-microglobulin (bikunin) was down-regulated, suggesting that pIC may suppress the acute phase response. Since the SSH libraries built for this study identified genes involved in the antiviral response, they are important resources for studying the responses of Atlantic cod to viruses. Evidence is provided for the existence of a RIG-I-like RNA helicase viral recognition pathway in Atlantic cod. Taken together, our data show that Atlantic cod can recognize double-stranded RNA and mount a rapid and potent interferon pathway response that is similar to that observed in other fish species and higher vertebrates.

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.

Transcriptomic responses of Atlantic salmon (Salmo salar) to environmental enrichment during juvenile rearing.

Captive rearing programs (hatcheries) are often used in conservation and management efforts for at-risk salmonid fish populations. However, hatcheries typically rear juveniles in environments that contrast starkly with natural conditions, which may lead to phenotypic and/or genetic changes that adversely affect the performance of juveniles upon their release to the wild. Environmental enrichment has been proposed as a mechanism to improve the efficacy of population restoration efforts from captive-rearing programs; in this study, we examine the influence of environmental enrichment during embryo and yolk-sac larval rearing on the transcriptome of Atlantic salmon (Salmo salar). Full siblings were reared in either a hatchery environment devoid of structure or an environment enriched with gravel substrate. At the end of endogenous feeding by juveniles, we examined patterns of gene transcript abundance in head tissues using the cGRASP-designed Agilent 4×44K microarray. Significance analysis of microarrays (SAM) indicated that 808 genes were differentially transcribed between the rearing environments and a total of 184 gene ontological (GO) terms were over- or under-represented in this gene list, several associated with mitosis/cell cycle and muscle and heart development. There were also pronounced differences among families in the degree of transcriptional response to rearing environment enrichment, suggesting that gene-by-environment effects, possibly related to parental origin, could influence the efficacy of enrichment interventions.