Climate change can impair bacterial pathogen defences in sablefish via hypoxia-mediated effects on adaptive immunity

Low-oxygen levels (hypoxia) in aquatic habitats are becoming more common because of global warming and eutrophication. However, the effects on the health/disease status of fishes, the world's largest group of vertebrates, are unclear. Therefore, we assessed how long-term hypoxia affected the immune function of sablefish, an ecologically and economically important North Pacific species, including the response to a formalin-killed Aeromonas salmonicida bacterin. Sablefish were held at normoxia or hypoxia (100% or 40% air saturated seawater, respectively) for 6-16 weeks, while we measured a diverse array of immunological traits. Given that the sablefish is a non-model organism, this involved the development of a species-specific methodological toolbox comprised of qPCR primers for 16 key immune genes, assays for blood antibacterial defences, the assessment of blood immunoglobulin (IgM) levels with ELISA, and flow cytometry and confocal microscopy techniques. We show that innate immune parameters were typically elevated in response to the bacterial antigens, but were not substantially affected by hypoxia. In contrast, hypoxia completely prevented the ∼1.5-fold increase in blood IgM level that was observed under normoxic conditions following bacterin exposure, implying a serious impairment of adaptive immunity. Since the sablefish is naturally hypoxia tolerant, our results demonstrate that climate change-related deoxygenation may be a serious threat to the immune competency of fishes.

Vegetable omega-3 and omega-6 fatty acids differentially modulate the antiviral and antibacterial immune responses of Atlantic salmon

The immunomodulatory effects of omega-3 and omega-6 fatty acids are a crucial subject of investigation for sustainable fish aquaculture, as fish oil is increasingly replaced by terrestrial vegetable oils in aquafeeds. Unlike previous research focusing on fish oil replacement with vegetable alternatives, our study explored how the omega-6 to omega-3 polyunsaturated fatty acid (PUFA) ratio in low-fish oil aquafeeds influences Atlantic salmon's antiviral and antibacterial immune responses. Atlantic salmon were fed aquafeeds rich in soy oil (high in omega-6) or linseed oil (high in omega-3) for 12 weeks and then challenged with bacterial (formalin-killed Aeromonas salmonicida) or viral-like (polyriboinosinic polyribocytidylic acid) antigens. The head kidneys of salmon fed high dietary omega-3 levels exhibited a more anti-inflammatory fatty acid profile and a restrained induction of pro-inflammatory and neutrophil-related genes during the immune challenges. The high-omega-3 diet also promoted a higher expression of genes associated with the interferon-mediated signaling pathway, potentially enhancing antiviral immunity. This research highlights the capacity of vegetable oils with different omega-6 to omega-3 PUFA ratios to modulate specific components of fish immune responses, offering insights for future research on the intricate lipid nutrition-immunity interplay and the development of novel sustainable low-fish oil clinical aquaculture feeds.

Characterization and transcript expression analyses of four Atlantic salmon (Salmo salar) serpinh1 paralogues provide evidence of evolutionary divergence

Atlantic salmon (Salmo salar) are not only the world's most economically important farmed fish in terms of total value, but also a salmonid, which means that they are invaluable for studies of the evolutionary fate of genes following multiple whole-genome duplication (WGD) events. In this study, four paralogues of the molecular chaperone serpinh1 were characterized in Atlantic salmon, as while this gene is considered to be a sensitive biomarker of heat stress in salmonids, mammalian studies have also identified it as being essential for collagen structural assembly and integrity. The four salmon paralogues were cloned and sequenced so that in silico analyses at the nucleotide and deduced amino acid levels could be performed. In addition, qPCR was used to measure: paralogue- and sex-specific constitutive serpinh1 expression across 17 adult tissues; and their expression in the liver and head kidney of male Atlantic salmon as affected by stress phenotype (high vs. low responder), increased temperature, and injection with a multi-valent vaccine. Compared to the other three paralogues, serpinh1a-2 had a unique constitutive expression profile across the 17 tissues. Although stress phenotype had minimal impact on the transcript expression of the four paralogues, injection with a commercial vaccine containing several formalin inactivated bacterins increased the expression of most paralogues (by 1.1 to 4.5-fold) across both tissues. At 20 °C, the expression levels of serpinh1a-1 and serpinh1a-2 were generally lower (by -1.1- to -1.6-fold), and serpinh1b-1 and serpinh1b-2 were 10.2- to 19.0-fold greater, in comparison to salmon held at 12 °C. With recent studies suggesting a putative link between serpinh1 and upper thermal tolerance in salmonids, the current research is a valuable first step in elucidating the potential mechanisms involved. This research: supports the use of serpinh1b-1 and serpinh1b-2 as a biomarkers of heat stress in salmon; and provides evidence of neo- and/or subfunctionalization between the paralogues, and important insights into how multiple genome duplication events can potentially lead to evolutionary divergence.

The late-evolving salmon and trout join the GnRH1 club

Although it is known that the whitefish, an ancient salmonid, expresses three distinct gonadotropin-releasing hormone (GnRH) forms in the brain, it has been thought that the later-evolving salmonids (salmon and trout) had only two types of GnRH: GnRH2 and GnRH3. We now provide evidence for the expression of GnRH1 in the gonads of Atlantic salmon by rapid amplification of cDNA ends, real-time quantitative PCR and immunohistochemistry. We examined six different salmonid genomes and found that each assembly has one gene that likely encodes a viable GnRH1 prepropeptide. In contrast to both functional GnRH2 and GnRH3 paralogs, the GnRH1 homeolog can no longer express the hormone. Furthermore, the viable salmonid GnRH1 mRNA is composed of only three exons, rather than the four exons that build the GnRH2 and GnRH3 mRNAs. Transcribed gnrh1 is broadly expressed (in 17/18 tissues examined), with relative abundance highest in the ovaries. Expression of the gnrh2 and gnrh3 mRNAs is more restricted, primarily to the brain, and not in the gonads. The GnRH1 proximal promoter presents composite binding elements that predict interactions with complexes that contain diverse cell fate and differentiation transcription factors. We provide immunological evidence for GnRH1 peptide in the nucleus of 1-year-old type A spermatogonia and cortical alveoli oocytes. GnRH1 peptide was not detected during other germ cell or reproductive stages. GnRH1 activity in the salmonid gonad may occur only during early stages of development and play a key role in a regulatory network that controls mitotic and/or meiotic processes within the germ cell.

Impact of stress phenotype, elevated temperature, and bacterin exposure on male Atlantic salmon (Salmo salar) growth, stress, and immune biomarker gene expression

In this study, postsmolt male Atlantic salmon, previously identified as low responders (LRs) or high responders (HRs) based on poststress cortisol levels, had their head kidney and liver sampled at 12°C and 20°C before injection (time 0) and after injection (i.e., at 12- and 24-h postinjection, respectively) with either Forte Micro (a multivalent vaccine containing bacterin, to capture peak antibacterial responses) or an equal volume of PBS. Quantitative real-time PCR (qPCR) was then used to measure the expression of 15 biomarker genes in the head kidney and 12 genes in the liver at each temperature/sampling point. Target transcripts were chosen that were related to growth, stress, and innate antibacterial immune responses. Many temperature, phenotype, and injection effects were found for individual genes within these three broad categories, and multivariate statistical analyses (i.e., principal component analysis and permutational multivariate analysis of variance) were used to look for overall patterns in transcript expression. These analyses revealed that HR salmon at 20°C mounted a more robust response (P < 0.05) for the 10 head kidney immune-related transcripts when injected with Forte Micro than LR salmon. In contrast, the seven liver stress-related transcripts displayed a greater response (P = 0.057) in LR versus HR fish with Forte Micro at 12°C. Overall, although this research did find some differences between LR and HR fish, it does not provide strong (conclusive) evidence that the selection of a particular phenotype would have major implications for the health of salmon over the temperature range examined.NEW & NOTEWORTHY This is the first paper to describe the impact of both temperature and bacterial stimulation on head kidney and liver transcript expression in Atlantic salmon characterized as LRs versus HRs. Notably, we found that HR salmon at 20°C mounted a more robust innate antibacterial immune response than LR salmon. In addition, LR fish at 12°C may (P = 0.057) exhibit higher expression of stress-related transcripts in response to vaccine injection relative to HR fish.

Basal and immune-responsive transcript expression of two Atlantic salmon interferon regulatory factor 2 (irf2) paralogues

Atlantic salmon (Salmo salar) is one of the most economically important aquaculture species globally. However, disease has become a prevalent threat to this industry. A thorough understanding of the genes and molecular pathways involved in the immune responses of Atlantic salmon is imperative for selective breeding of disease-resistant broodstock, as well as developing new diets and vaccines to mitigate the impact of disease. Members of the interferon regulatory factor (IRF) family of transcription factors play roles in the induction of interferons and other cytokines involved in host immune responses to intracellular and parasitic pathogens. IRF family members also play diverse roles in other biological processes, such as stress response, reproduction and development. The current study focused on one member of the IRF family: interferon regulatory factor 2 (irf2). As previously shown, due to the genome duplication that occurred ∼80 million years ago in the salmonid lineage, there are two irf2 paralogues in the Atlantic salmon genome. In silico analyses at the cDNA and deduced amino acid levels were conducted followed by phylogenetic tree construction with IRF2 amino acid sequences from various ray-finned fishes, cartilaginous fish and tetrapods. qPCR was then used to analyze paralogue-specific irf2 constitutive expression across 17 adult tissues, as well as responses to the viral mimic pIC (i.e., synthetic double-stranded RNA analog) in cultured macrophage-like cells (in vitro) and to infection with the Gram-negative bacterium Moritella viscosa in skin samples (in vivo). The qPCR studies showed sex- and paralogue-specific differences in expression across tissues. For example, expression of both paralogues was higher in ovary than in testes; expression (considering both sexes together) was highest for irf2-1 in gonad and for irf2-2 in hindgut. Both irf2 paralogues were responsive to pIC stimulation, but varied in their induction level, with irf2-1 having an overall stronger response than irf2-2. Only one paralogue, irf2-2, was significantly responsive to M. viscosa infection. Differences in irf2-1 and irf2-2 transcript expression levels constitutively across tissues, and in response to pIC and M. viscosa, may suggest neo- or subfunctionalization of the duplicated genes. This novel information expands current knowledge and provides insight into how genome duplication events may impact host regulation of important immune markers.

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.

Influence of Supplemental Dietary Cholesterol on Growth Performance, Indices of Stress, Fillet Pigmentation, and Upper Thermal Tolerance of Female Triploid Atlantic Salmon (Salmo salar)

The salmon aquaculture industry must be proactive at developing mitigation tools/strategies to offset the potential negative impacts of climate change. Therefore, this study examined if additional dietary cholesterol could enhance salmon production at elevated temperatures. We hypothesized that supplemental cholesterol could aid in maintaining cell rigidity, reducing stress and the need to mobilize astaxanthin muscle stores, and improving salmon growth and survival at high rearing temperatures. Accordingly, postsmolt female triploid salmon were exposed to an incremental temperature challenge (+0.2°C day^-1) to mimic conditions that they experience in sea cages in the summer, with temperature held at both 16 and 18°C for several weeks [i.e., 3 weeks at 16°C, followed by an increase at 0.2°C day^-1 to 18°C (10 days), then 5 weeks at 18°C] to prolong their exposure to elevated temperatures. From 16°C onwards, the fish were fed either a control diet, or one of two nutritionally equivalent experimental diets containing supplemental cholesterol [+1.30%, experimental diet ^#1 (ED1); or +1.76%, experimental diet ^#2 (ED2)]. Adding cholesterol to the diet did not affect the salmon's incremental thermal maximum (IT_Max), growth, plasma cortisol, or liver stress-related transcript expression. However, ED2 appeared to have a small negative impact on survival, and both ED1 and ED2 reduced fillet "bleaching" above 18°C as measured using SalmoFan™ scores. Although the current results suggest that supplementing salmon diets with cholesterol would have few/minimal benefits for the industry, ≤ 5% of the female triploid Atlantic salmon used in this study irrespective of diet died before temperature reached 22°C. These latter data suggest that it is possible to produce all female populations of reproductively sterile salmon that can withstand summer temperatures in Atlantic Canada.

Nutritional immunomodulation of Atlantic salmon response to Renibacterium salmoninarum bacterin

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

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.

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

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

Atlantic Salmon (Salmo salar) bacterial and viral innate immune responses are not impaired by florfenicol or tetracycline administration

Antibiotics are used to treat bacterial infections in fish aquaculture, and these drugs can interact with immune cells/the immune system and potentially leave fish vulnerable to viral, fungal, parasitic, or other bacterial infections. However, the effects of antibiotics on fish immunity have largely been overlooked by the aquaculture industry. We tested, at 12 and 20 °C, whether tetracycline and florfenicol (the most commonly used antibiotics in commercial aquaculture), affected the Atlantic salmon's capacity to respond to bacterial or viral stimulation. Atlantic salmon were acclimated to 12 or 20 °C and fed with tetracycline or florfenicol (100 and 10 mg kg of body weight^-1 day^-1, respectively) medicated feed for 15 or 10 days, respectively. Thereafter, we evaluated their immune function prior to, and after, an intraperitoneal injection of Forte Micro (containing inactivated cultures of Aeromonas salmonicida, Vibrio anguillarum, Vibrio ordalii and Vibrio salmonicida) or the viral mimic polyriboinosinic polyribocytidylic acid (pIC). We measured the transcript expression levels of 8 anti-bacterial and 8 anti-viral putative biomarker genes, and the innate (leukocyte respiratory burst, plasma lysozyme activity and hemolytic activity of the alternative complement pathway) and cellular (relative number of erythrocytes, lymphocytes and thrombocytes, and granulocytes such as monocytes and neutrophils) responses to these challenges. Overall, we only found a few minor effects of either tetracycline or florfenicol on immune gene expression or function at either temperature. Although several studies have reported that antibiotics may negatively affect fish immune responses, our results show that industry-relevant dietary tetracycline and florfenicol treatments do not substantially impact the salmon's innate immune responses. Currently, this is the most comprehensive study on the effects of antibiotics administrated according to industry protocols on immune function in Atlantic salmon.

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

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

Global gene expression responses of Atlantic salmon skin to Moritella viscosa

Moritella viscosa is a Gram-negative pathogen that causes large, chronic ulcers, known as winter-ulcer disease, in the skin of several fish species including Atlantic salmon. We used a bath challenge approach to profile the transcriptome responses of M. viscosa-infected Atlantic salmon skin at the lesion (Mv-At) and away from the lesion (Mv-Aw) sites. M. viscosa infection was confirmed through RNA-based qPCR assays. RNA-Seq identified 5212 and 2911 transcripts differentially expressed in the Mv-At compared to no-infection control and Mv-Aw groups, respectively. Also, there were 563 differentially expressed transcripts when comparing the Mv-Aw to control samples. Our results suggest that M. viscosa caused massive and strong, but largely infection site-focused, transcriptome dysregulations in Atlantic salmon skin, and its effects beyond the skin lesion site were comparably subtle. The M. viscosa-induced transcripts of Atlantic salmon were mainly involved in innate and adaptive immune response-related pathways, whereas the suppressed transcripts by this pathogen were largely connected to developmental and cellular processes. As validated by qPCR, M. viscosa dysregulated transcripts encoding receptors, signal transducers, transcription factors and immune effectors playing roles in TLR- and IFN-dependent pathways as well as immunoregulation, antigen presentation and T-cell development. This study broadened the current understanding of molecular pathways underlying M. viscosa-triggered responses of Atlantic salmon, and identified biomarkers that may assist to diagnose and combat this pathogen.

Distinct early life stage gene expression effects of hybridization among European and North American farmed and wild Atlantic salmon populations

Due to multi-generation domestication selection, farmed and wild Atlantic salmon diverge genetically, which raises concerns about potential genetic interactions among escaped farmed and wild populations and disruption of local adaptation through introgression. When farmed strains of distant geographic origin are used, it is unknown whether the genetic consequences posed by escaped farmed fish will be greater than if more locally derived strains are used. Quantifying gene transcript expression differences among divergent farmed, wild and F₁ hybrids under controlled conditions is one of the ways to explore the consequences of hybridization. We compared the transcriptomes of fry at the end of yolk sac absorption of a European (EO) farmed ("StofnFiskur", Norwegian strain), a North American (NA) farmed (Saint John River, NB strain), a Newfoundland (NF) wild population with EO ancestry, and related F₁ hybrids using 44K microarrays. Our findings indicate that the wild population showed greater transcriptome differences from the EO farmed strain than that of the NA farmed strain. We also found the largest differences in global gene expression between the two farmed strains. We detected the fewest differentially expressed transcripts between F₁ hybrids and domesticated/wild maternal strains. We also found that the differentially expressed genes between cross types over-represented GO terms associated with metabolism, development, growth, immune response, and redox homeostasis processes. These findings suggest that the interbreeding of escaped EO/NA farmed and NF wild population would alter gene transcription, and the consequences of hybridization would be greater from escaped EO farmed than NA farmed salmon, resulting in potential effects on the wild populations.

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 Profiling of Atlantic Salmon (Salmo salar) Parr With Higher and Lower Pathogen Loads Following Piscirickettsia salmonis Infection

Salmonid rickettsial septicemia (SRS), caused by Piscirickettsia salmonis, is one of the most devastating diseases of salmonids. However, the transcriptomic responses of Atlantic salmon (Salmon salar) in freshwater to an EM-90-like isolate have not been explored. Here, we infected Atlantic salmon parr with an EM-90-like isolate and conducted time-course qPCR analyses of pathogen load and four biomarkers (campb, hampa, il8a, tlr5a) of innate immunity on the head kidney samples. Transcript expression of three of these genes (except hampa), as well as pathogen level, peaked at 21 days post-injection (DPI). Multivariate analyses of infected individuals at 21 DPI revealed two infection phenotypes [lower (L-SRS) and higher (H-SRS) infection level]. Five fish from each group (Control, L-SRS, and H-SRS) were selected for transcriptome profiling using a 44K salmonid microarray platform. We identified 1,636 and 3,076 differentially expressed probes (DEPs) in the L-SRS and H-SRS groups compared with the control group, respectively (FDR = 1%). Gene ontology term enrichment analyses of SRS-responsive genes revealed the activation of a large number of innate (e.g. “phagocytosis”, “defense response to bacterium”, “inflammatory response”) and adaptive (e.g. “regulation of T cell activation”, “antigen processing and presentation of exogenous antigen”) immune processes, while a small number of general physiological processes (e.g. “apoptotic process”, development and metabolism relevant) was enriched. Transcriptome results were confirmed by qPCR analyses of 42 microarray-identified transcripts. Furthermore, the comparison of individuals with differing levels of infection (H-SRS vs. L-SRS) generated insights into the biological processes possibly involved in disease resistance or susceptibility. This study demonstrated a low mortality (~30%) EM-90-like infection model and broadened the current understanding of molecular pathways underlying P. salmonis-triggered responses of Atlantic salmon, identifying biomarkers that may assist to diagnose and combat this pathogen.

Characterization of miRNAs in Embryonic, Larval, and Adult Lumpfish Provides a Reference miRNAome for Cyclopterus lumpus

Simple Summary

Lumpfish (Cyclopterus lumpus) is an emergent aquaculture species, and its miRNA repertoire is still unknown. miRNAs are critical post-transcriptional modulators of teleost gene expression. Therefore, a lumpfish reference miRNAome was characterized by small RNA sequencing and miRDeep analysis of samples from different organs and developmental stages. The resulting miRNAome, an essential reference for future expression analyses, consists of 443 unique mature miRNAs from 391 conserved and eight novel miRNA genes. Enrichment of specific miRNAs in particular organs and developmental stages indicates that some conserved lumpfish miRNAs regulate organ and developmental stage-specific functions reported in other teleosts.

Abstract

MicroRNAs (miRNAs) are endogenous small RNA molecules involved in the post-transcriptional regulation of protein expression by binding to the mRNA of target genes. They are key regulators in teleost development, maintenance of tissue-specific functions, and immune responses. Lumpfish (Cyclopterus lumpus) is becoming an emergent aquaculture species as it has been utilized as a cleaner fish to biocontrol sea lice (e.g., Lepeophtheirus salmonis) infestation in the Atlantic Salmon (Salmo salar) aquaculture. The lumpfish miRNAs repertoire is unknown. This study identified and characterized miRNA encoding genes in lumpfish from three developmental stages (adult, embryos, and larvae). A total of 16 samples from six different adult lumpfish organs (spleen, liver, head kidney, brain, muscle, and gill), embryos, and larvae were individually small RNA sequenced. Altogether, 391 conserved miRNA precursor sequences (discovered in the majority of teleost fish species reported in miRbase), eight novel miRNA precursor sequences (so far only discovered in lumpfish), and 443 unique mature miRNAs were identified. Transcriptomics analysis suggested organ-specific and age-specific expression of miRNAs (e.g., miR-122-1-5p specific of the liver). Most of the miRNAs found in lumpfish are conserved in teleost and higher vertebrates, suggesting an essential and common role across teleost and higher vertebrates. This study is the first miRNA characterization of lumpfish that provides the reference miRNAome for future functional studies.

Reverse Transcription-Quantitative Real-Time Polymerase Chain Reaction (RT-qPCR) for Gene Expression Analyses

The reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) is considered to be the gold standard for gene expression research. However, for this claim to be valid, RT-qPCR studies must test and optimize the quality of its RNA templates and assays. This chapter describes the experimental procedures required to generate reliable and reproducible gene expression results using RT-qPCR.

Phenotypic stress response does not influence the upper thermal tolerance of male Atlantic salmon (Salmo salar)

Fish can be identified as either low responders (LR) or high responders (HR) based on post-stress cortisol levels and whether they exhibit a proactive or reactive stress coping style, respectively. In this study, male Atlantic salmon (Salmo salar) from 17 families reared at 9 °C were repeatedly exposed to an acute handling stress over a period of four months, with plasma cortisol levels measured at 1 h post-stress. Fish were identified as either LR or HR if the total Z-score calculated from their cortisol responses fell into the lower or upper quartile ranges, respectively; with intermediate responders (IR) classified as the remainder. Salmon characterized as LR, IR or HR were then subjected to an incremental thermal challenge, where temperature was raised at 0.2 °C day^-1 from their acclimation temperature (12 °C) to mimic natural sea-cage farming conditions during the summer in Newfoundland. Interestingly, feed intake remained high up to 22 °C, while previous studies have shown a decrease in salmon appetite after ∼16-18 °C. After the first three mortalities were recorded at elevated temperature, a subset of LR and HR salmon were exposed to another acute handling stress event at 23.6 °C. Basal and post-stress measurements of plasma cortisol, glucose and lactate did not differ between stress response phenotypes at this temperature. In the end, the average incremental thermal maximum (IT_Max) of LR and HR fish was not different (25.1 °C). In comparison, the critical thermal maximum (CT_Max; temperature increased at 2 °C h^-1) of the remaining IR fish that had been held at 12 °C was 28.5 °C. Collectively, these results: 1) show that this population of Atlantic salmon is very thermally tolerant, and further question the relevance of CT_Max in assessing responses to real-world temperature changes; and 2) indicate that characterization of stress phenotype at 9 °C is not predictive of their stress response or survival at high temperatures. Therefore, selection of fish based on phenotypic stress response at low temperatures may not be beneficial to incorporate into Atlantic salmon breeding programs, especially if the goal is to improve growth performance and survival at high temperatures in sea-cages.

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

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

The Atlantic salmon's stress- and immune-related transcriptional responses to moderate hypoxia, an incremental temperature increase, and these challenges combined

The marine environment is predicted to become warmer, and more hypoxic, and these conditions may negatively impact the health and survival of coastal fish species, including wild and farmed Atlantic salmon (Salmo salar). Thus, we examined how: (1) moderate hypoxia (∼70% air saturation) at 12°C for 3 weeks; (2) an incremental temperature increase from 12°C to 20°C (at 1°C week-1) followed by 4 weeks at 20°C; and (3) treatment "2" combined with moderate hypoxia affected transcript expression in the liver of post-smolts as compared to control conditions (normoxia, 12°C). Specifically, we assessed the expression of 45 genes related to the heat shock response, oxidative stress, apoptosis, metabolism and immunity using a high-throughput qPCR approach (Fluidigm Biomark™ HD). The expression profiles of 27 "stress"-related genes indicated that: (i) moderate hypoxia affected the expression of several stress genes at 12°C; (ii) their expression was impacted by 16°C under normoxic conditions, and this effect increased until 20°C; (iii) the effects of moderate hypoxia were not additive to those at temperatures above 16°C; and (iv) long-term (4 weeks) exposure to 20°C, with or without hypoxia, resulted in a limited acclimatory response. In contrast, the expression of 15 immune-related genes was not greatly affected until temperatures reached 20°C, and this effect was particularly evident in fish exposed to the added challenge of hypoxia. These results provide valuable information on how these two important environmental factors affect the "stress" physiology and immunology of Atlantic salmon, and we identify genes that may be useful as hypoxia and/or temperature biomarkers in salmonids and other fishes.

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

Simple Summary

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

Abstract

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

RNA-seq analysis of the mantle transcriptome from Mytilus edulis during a seasonal spawning event in deep and shallow water culture sites on the northeast coast of Newfoundland, Canada

The blue mussel (Mytilus edulis) has global commercial and ecological importance both in wild and cultured conditions. However there is a qualitative and quantitative lack of knowledge of the molecular mechanisms associated with its reproductive physiology, especially with reference to environmental interactions. Here we initiated a transcriptomic analysis (RNA-sequencing (RNA-seq)) of the mantle from both sexes sampled during a seasonal spawning event and from two culture depths (shallow-5 m; deep- 15 m). Mantle libraries were produced from 3 males and 3 females sampled from each of two shallow sites and two deep sites for a total of 12 replicate male and 12 replicate female libraries (24 total libraries). Overall a total of 2.3 billion raw 100 base reads with an average of 96.5 million reads/library were obtained and assembled into 296,118 transcripts with an average length of 568 bp. Overall, 315 transcripts from male libraries and 25 from female libraries were found to be upregulated in deep water as compared to shallow (edgeR adjusted p value ≤ 0.05). Conversely, 126 transcripts from male libraries and 135 from female libraries were found to be significantly downregulated at the same depth. Thirteen transcripts were selected for qPCR validation based on importance in reproduction, antimicrobial defense and metabolism. Of these, 9 RNA-seq identified transcripts were shown by qPCR to be differentially expressed between groups: 2 were upregulated in deep compared with shallow water (dhx38, mt-co1), 2 were upregulated for female compared with male mantle (pias2, mapkap1) and 6 genes (fndc3a, acbd3, klhl10, ccnb3, armc4, mt-co1) showed to be upregulated in males compared to females. The majority of qPCR studied transcripts were identified as involved in gamete development based on the UniProt database. This study further characterizes the importance of the mantle transcriptome during reproductive activities of M. edulis.