A 44K microarray dataset of the changing transcriptome in developing Atlantic salmon (Salmo salar L.)

Transcriptome Profiling of Gene Expression in Atlantic Salmon (Salmo salar) at Early Stage of Development

For Atlantic salmon development, the most critical phase is the early development stage from egg to fry through alevin. However, the studies investigating the early development of Atlantic salmon based on RNA-seq are scarce and focus only on one stage of development. Therefore, using the RNA-seq technology, the assessment of different gene expressions of various early development stages (egg, alevin, and fry) was performed on a global scale. Over 22 GB of clean data was generated from 9 libraries with three replicates for each stage with over 90% mapping efficiency. A total of 5534 genes were differentially expressed, among which 19, 606, and 826 genes were specifically expressed in each stage, respectively. The transcriptome analysis showed that the number of differentially expressed genes (DEGs) increased as the Atlantic salmon progressed in development from egg to fry stage. In addition, gene ontology enrichment demonstrated that egg and alevin stages are characterized by upregulation of genes involved in spinal cord development, neuron projection morphogenesis, axonogenesis, and cytoplasmic translation. At the fry stage, upregulated genes were enriched in the muscle development process (muscle cell development, striated muscle cell differentiation, and muscle tissue development), immune system (defense response and canonical NF-kappaB signal transduction), as well as epidermis development. These results suggest that the early development of Atlantic salmon is characterized by a dynamic shift in gene expression and DEGs between different stages, which provided a solid foundation for the investigation of Atlantic salmon development.

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.

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.

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 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.

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

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

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

Simple Summary

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

Abstract

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

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Transcriptomic Profiling of the Adaptive and Innate Immune Responses of Atlantic Salmon to Renibacterium salmoninarum Infection

Bacterial Kidney Disease (BKD), which is caused by a Gram-positive, intracellular bacterial pathogen (Renibacterium salmoninarum), affects salmonids including Atlantic salmon (Salmo salar). However, the transcriptome response of Atlantic salmon to BKD remained unknown before the current study. We used a 44K salmonid microarray platform to characterise the global gene expression response of Atlantic salmon to BKD. Fish (~54 g) were injected with a dose of R. salmoninarum (H-2 strain, 2 × 10⁸ CFU per fish) or sterile medium (control), and then head kidney samples were collected at 13 days post-infection/injection (dpi). Firstly, infection levels of individuals were determined through quantifying the R. salmoninarum level by RNA-based TaqMan qPCR assays. Thereafter, based on the qPCR results for infection level, fish (n = 5) that showed no (control), higher (H-BKD), or lower (L-BKD) infection level at 13 dpi were subjected to microarray analyses. We identified 6,766 and 7,729 differentially expressed probes in the H-BKD and L-BKD groups, respectively. There were 357 probes responsive to the infection level (H-BKD vs. L-BKD). Several adaptive and innate immune processes were dysregulated in R. salmoninarum-infected Atlantic salmon. Adaptive immune pathways associated with lymphocyte differentiation and activation (e.g., lymphocyte chemotaxis, T-cell activation, and immunoglobulin secretion), as well as antigen-presenting cell functions, were shown to be differentially regulated in response to BKD. The infection level-responsive transcripts were related to several mechanisms such as the JAK-STAT signalling pathway, B-cell differentiation and interleukin-1 responses. Sixty-five microarray-identified transcripts were subjected to qPCR validation, and they showed the same fold-change direction as microarray results. The qPCR-validated transcripts studied herein play putative roles in various immune processes including pathogen recognition (e.g., tlr5), antibacterial activity (e.g., hamp and camp), regulation of immune responses (e.g., tnfrsf11b and socs1), T-/B-cell differentiation (e.g., ccl4, irf1 and ccr5), T-cell functions (e.g., rnf144a, il13ra1b and tnfrsf6b), and antigen-presenting cell functions (e.g., fcgr1). The present study revealed diverse immune mechanisms dysregulated by R. salmoninarum in Atlantic salmon, and enhanced the current understanding of Atlantic salmon response to BKD. The identified biomarker genes can be used for future studies on improving the resistance of Atlantic salmon to BKD.

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

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

Transcriptomic Profiling in Fins of Atlantic Salmon Parasitized with Sea Lice: Evidence for an Early Imbalance Between Chalimus-Induced Immunomodulation and the Host's Defense Response

Parasitic sea lice (e.g., Lepeophtheirus salmonis) cause costly outbreaks in salmon farming. Molecular insights into parasite-induced host responses will provide the basis for improved management strategies. We investigated the early transcriptomic responses in pelvic fins of Atlantic salmon parasitized with chalimus I stage sea lice. Fin samples collected from non-infected (i.e. pre-infected) control (PRE) and at chalimus-attachment sites (ATT) and adjacent to chalimus-attachment sites (ADJ) from infected fish were used in profiling global gene expression using 44 K microarrays. We identified 6568 differentially expressed probes (DEPs, FDR < 5%) that included 1928 shared DEPs between ATT and ADJ compared to PRE. The ATT versus ADJ comparison revealed 90 DEPs, all of which were upregulated in ATT samples. Gene ontology/pathway term network analyses revealed profound changes in physiological processes, including extracellular matrix (ECM) degradation, tissue repair/remodeling and wound healing, immunity and defense, chemotaxis and signaling, antiviral response, and redox homeostasis in infected fins. The QPCR analysis of 37 microarray-identified transcripts representing these functional themes served to confirm the microarray results with a significant positive correlation (p < 0.0001). Most immune/defense-relevant transcripts were downregulated in both ATT and ADJ sites compared to PRE, suggesting that chalimus exerts immunosuppressive effects in the salmon's fins. The comparison between ATT and ADJ sites demonstrated the upregulation of a suite of immune-relevant transcripts, evidencing the salmon's attempt to mount an anti-lice response. We hypothesize that an imbalance between immunomodulation caused by chalimus during the early phase of infection and weak defense response manifested by Atlantic salmon makes it a susceptible host for L. salmonis.

Profiling the transcriptome response of Atlantic salmon head kidney to formalin-killed Renibacterium salmoninarum

Renibacterium salmoninarum is a Gram-positive, intracellular bacterial pathogen that causes Bacterial Kidney Disease (BKD) in Atlantic salmon (Salmo salar). The host transcriptomic response to this immune-suppressive pathogen remains poorly understood. To identify R. salmoninarum-responsive genes, Atlantic salmon were intraperitoneally injected with a low (5 × 10⁵ cells/kg, Low-Rs) or high (5 × 10⁷ cells/kg; High-Rs) dose of formalin-killed R. salmoninarum bacterin or phosphate-buffered saline (PBS control); head kidney samples were collected before and 24 h after injection. Using 44K microarray analysis, we identified 107 and 345 differentially expressed probes in response to R. salmoninarum bacterin (i.e. High-Rs vs. PBS control) by Significance Analysis of Microarrays (SAM) and Rank Products (RP), respectively. Twenty-two microarray-identified genes were subjected to qPCR assays, and 17 genes were confirmed as being significantly responsive to the bacterin. There was an up-regulation in expression of genes playing putative roles as immune receptors and antimicrobial effectors. Genes with putative roles as pathogen recognition (e.g. clec12b and tlr5) or immunoregulatory (e.g. tnfrsf6b and tnfrsf11b) receptors were up-regulated in response to R.salmoninarum bacterin. Also, chemokines and a chemokine receptor showed opposite regulation [up-regulation of effectors (i.e. ccl13 and ccl) and down-regulation of cxcr1] in response to the bacterin. The present study identified and validated novel biomarker genes (e.g. ctsl1, lipe, cldn4, ccny) that can be used to assess Atlantic salmon response to R. salmoninarum, and will be valuable in the development of tools to combat BKD.

The initiation of puberty in Atlantic salmon brings about large changes in testicular gene expression that are modulated by the energy status

Background

When puberty starts before males reach harvest size, animal welfare and sustainability issues occur in Atlantic salmon (Salmo salar) aquaculture. Hallmarks of male puberty are an increased proliferation activity in the testis and elevated androgen production. Examining transcriptional changes in salmon testis during the transition from immature to maturing testes may help understanding the regulation of puberty, potentially leading to procedures to modulate its start. Since differences in body weight influence, via unknown mechanisms, the chances for entering puberty, we used two feed rations to create body weight differences.

Results

Maturing testes were characterized by an elevated proliferation activity of Sertoli cells and of single undifferentiated spermatogonia. Pituitary gene expression data suggest increased Gnrh receptor and gonadotropin gene expression, potentially responsible for the elevated circulating androgen levels in maturing fish. Transcriptional changes in maturing testes included a broad variety of signaling systems (e.g. Tgfβ, Wnt, insulin/Igf, nuclear receptors), but also, activation of metabolic pathways such as anaerobic metabolism and protection against ROS. Feed restriction lowered the incidence of puberty. In males maturing despite feed restriction, plasma androgen levels were higher than in maturing fish receiving the full ration. A group of 449 genes that were up-regulated in maturing fully fed fish, was up-regulated more prominently in testis from fish maturing under caloric restriction. Moreover, 421 genes were specifically up-regulated in testes from fish maturing under caloric restriction, including carbon metabolism genes, a pathway relevant for nucleotide biosynthesis and for placing epigenetic marks.

Conclusions

Undifferentiated spermatogonia and Sertoli cell populations increased at the beginning of puberty, which was associated with the up-regulation of metabolic pathways (e.g. anaerobic and ROS pathways) known from other stem cell systems. The higher androgen levels in males maturing under caloric restriction may be responsible for the stronger up-regulation of a common set of (449) maturation-associated genes, and the specific up-regulation of another set of (421) genes. The latter opened regulatory and/or metabolic options for initiating puberty despite feed restriction. As a means to reduce the incidence of male puberty in salmon, however, caloric restriction seems unsuitable.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5869-9) contains supplementary material, which is available to authorized users.

A pathway-focused RT-qPCR array study on immune relevant genes in rainbow trout (Oncorhynchus mykiss) harboring cecropin P1 transgene

Recently, our laboratory had produced five families of transgenic rainbow trout harboring cecropin P1 transgene, and via repeated challenge studies these fish exhibited a significant elevation of resistance to infection by microbial pathogens. By cDNA microarray and mRNA deep sequencing (mRNA-seq) analyses on two of the five families of cecropin P1 transgenic fish, differentially expressed genes (DEGs) relevant to the innate and adaptive immune pathways in three different immune-related tissues, (i.e. spleen, kidney and liver) were profiled. These results supported our hypothesis that in addition to its direct microbicidal activity, the transgene product of cecropin P1 induces immunomodulatory activity in the transgenic host. Here, we have adapted the technique of quantitative reverse transcription real time PCR (RT-qPCR) array to analyze the expression of genes relevant to the innate and adaptive immune pathways in the rest three families. A RT-qPCR array was constructed with oligonucleotide primers of fifty-two innate/adaptive immune relevant DEGs shown to be the most perturbed by cecropin P1 transgene product in previous studies. Messenger RNA isolated from the spleen, kidney and liver of transgenic fish and non-transgenic fish control were studied on this array. Results of RT-qPCR array revealed that statistically significant perturbations of gene expression were detected in pathways of cytokine/chemokine signaling, Toll-like receptor signaling, complement cascade, antigen processing/presentation, lysosomal phagocytosis and leukocyte trans-endothelial migration in the transgenic spleen; extracellular matrix (ECM) organization and leukocyte trans-endothelial migration pathways in the transgenic kidney; lysosomal activity pathway in the transgenic liver. Furthermore, genes related to the pathways of the peroxisome proliferator-activated receptors (PPAR) signaling, lipid metabolism process and arachidonic acid metabolism were also impacted in the transgenic liver. Findings of the current study are in good agreement with those discoveries in previous two transgenic families by cDNA microarray and mRNA-seq analyses.

Characterization of the European Sea Bass (Dicentrarchus labrax) Gonadal Transcriptome During Sexual Development

The European sea bass is one of the most important cultured fish in Europe and has a marked sexual growth dimorphism in favor of females. It is a gonochoristic species with polygenic sex determination, where a combination between still undifferentiated genetic factors and environmental temperature determines sex ratios. The molecular mechanisms responsible for gonadal sex differentiation are still unknown. Here, we sampled fish during the gonadal developmental period (110 to 350 days post fertilization, dpf), and performed a comprehensive transcriptomic study by using a species-specific microarray. This analysis uncovered sex-specific gonadal transcriptomic profiles at each stage of development, identifying larger number of differentially expressed genes in ovaries when compared to testis. The expression patterns of 54 reproduction-related genes were analyzed. We found that hsd17β10 is a reliable marker of early ovarian differentiation. Further, three genes, pdgfb, snx1, and nfy, not previously related to fish sex differentiation, were tightly associated with testis development in the sea bass. Regarding signaling pathways, lysine degradation, bladder cancer, and NOD-like receptor signaling were enriched for ovarian development while eight pathways including basal transcription factors and steroid biosynthesis were enriched for testis development. Analysis of the transcription factor abundance showed an earlier increase in females than in males. Our results show that, although many players in the sex differentiation pathways are conserved among species, there are peculiarities in gene expression worth exploring. The genes identified in this study illustrate the diversity of players involved in fish sex differentiation and can become potential biomarkers for the management of sex ratios in the European sea bass and perhaps other cultured species.

Changes in the liver transcriptome of farmed Atlantic salmon (Salmo salar) fed experimental diets based on terrestrial alternatives to fish meal and fish oil

Background

Dependence on marine natural resources threatens the sustainability of Atlantic salmon aquaculture. In the present study, Atlantic salmon fed for 14 weeks with an experimental diet based on animal by-products and vegetable oil (ABP) exhibited reduced growth performance compared with others fed a fish meal/fish oil based experimental diet (MAR) and a plant protein/vegetable oil-based experimental diet (VEG). To characterize the molecular changes underlying the differences in growth performance, we conducted a 44 K microarray study of the liver transcriptome of the three dietary groups.

Results

The microarray experiment identified 122 differentially expressed features (Rank Products, PFP < 10%). Based on their associated Gene Ontology terms, 46 probes were classified as metabolic and growth-relevant genes, 25 as immune-related, and 12 as related to oxidation-reduction processes. The microarray results were validated by qPCR analysis of 29 microarray-identified transcripts. Diets significantly modulated the transcription of genes involved in carbohydrate metabolism (gck and pfkfb4), cell growth and proliferation (sgk2 and htra1), apoptosis (gadd45b), lipid metabolism (fabp3, idi1, sqs), and immunity (igd, mx, ifit5, and mhcI). Hierarchical clustering and linear correlation analyses were performed to find gene expression patterns among the qPCR-analyzed transcripts, and connections between them and muscle and liver lipid composition. Overall, our results indicate that changes in the liver transcriptome and tissue lipid composition were driven by cholesterol synthesis up-regulation by ABP and VEG diets, and the lower carbohydrate intake in the ABP group. Two of the microarray-identified genes (sgk2 and htra1) might be key to explaining glucose metabolism regulation and the dietary-modulation of the immune system in fish. To evaluate the potential of these genes as predictive biomarkers, we subjected the qPCR data to a stepwise discriminant analysis. Three sets of no more than four genes were found to be able to predict, with high accuracy (67–94%), salmon growth and fatty acid composition.

Conclusions

This study provides new findings on the impact of terrestrial animal and plant products on the nutrition and health of farmed Atlantic salmon, and a new method based on gene biomarkers for potentially predicting desired phenotypes, which could help formulate superior feeds for the Atlantic salmon aquaculture industry.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5188-6) contains supplementary material, which is available to authorized users.

RNA-Seq analysis of differentially expressed genes relevant to innate and adaptive immunity in cecropin P1 transgenic rainbow trout (Oncorhynchus mykiss)

Background

In the past years, our laboratory successfully generated transgenic rainbow trout bearing cecropin P1 transgene. These fish exhibited resistant characteristic to infection by Aeromonas salmonicida, Infectious Hematopoietic Necrosis Virus (IHNV) and Ceratomyxa shasta (a parasitic pathogen). Previously, treating rainbow trout macrophage cells (RTS-11) with cecropin B, pleurocidin and CF17, respectively, resulted in elevated expression of two pro-inflammatory genes, e.g. cyclooxygenase-2 (cox-2) and interleukin-1β (il-1β). In addition, a profiling of global gene expression by 44 k salmonid microarray analysis was conducted, and the results showed that immune relevant processes have been perturbed in cecopin P1 transgenic rainbow trout. Therefore, we hypothesized that cecropin P1 may not only eliminate pathogens directly, but also modulate the host immune systems, leading to increased resistance against pathogen infections. To confirm this hypothesis, we performed de novo mRNA deep sequencing (RNA-Seq) to analyze the transcriptomic expression profiles in three immune competent tissues of cecropin P1 transgenic rainbow trout.

Results

De novo sequencing of mRNA of the rainbow trout spleen, liver and kidney tissues were conducted by second-generation Illumina system, followed by Trinity assembly. Tissue specific unigenes were obtained, and annotated according to the Gene Ontology (GO) and the Nucleotide Basic Local Alignment Search Tool (BLAST). Over 2000 differentially expressed genes (DEGs) were determined by normalized ratio of Reads Per Kilobase of transcript per million mapped reads (RPKM) among the transgenic and non-transgenic fish in a tissue specific manner, and there were 82 DEGs in common among the three tissues. In addition, the enrichment analysis according to Gene Ontology Biological Process (GO:BP), and Kyoto Encyclopedia of Genes and Genomes (KEGG) based pathway analysis associated with innate/adaptive immunity of fish were also performed to illustrate the altered immune-related functions in each tissue.

Conclusions

According to the RNA-Seq data, the correlations between alteration of gene expression profiles and the functional perturbations of the host immune processes were revealed. In comparison with the results of cDNA microarray analysis conducted by Lo et al., the overall results supported our hypothesis that the gene product of cecropin P1 transgene may not only directly eliminate pathogens, but also modulate the host immune system. Results of this study present valuable genetic information for Oncorhynchus mykiss, and will benefit future studies on the immunology of this fish species.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5141-8) contains supplementary material, which is available to authorized users.

Transcriptome profiling of antiviral immune and dietary fatty acid dependent responses of Atlantic salmon macrophage-like cells

BACKGROUND

Due to the limited availability and high cost of fish oil in the face of increasing aquaculture production, there is a need to reduce usage of fish oil in aquafeeds without compromising farm fish health. Therefore, the present study was conducted to determine if different levels of vegetable and fish oils can alter antiviral responses of salmon macrophage-like cells (MLCs). Atlantic salmon (Salmo salar) were fed diets containing 7.4% (FO7) or 5.1% (FO5) fish oil. These diets were designed to be relatively low in EPA + DHA (i.e. FO7: 1.41% and FO5: 1%), but near the requirement level, and resulting in comparable growth. Vegetable oil (i.e. rapeseed oil) was used to balance fish oil in experimental diets. After a 16-week feeding trial, MLCs isolated from fish in these dietary groups were stimulated by a viral mimic (dsRNA: pIC) for 6 h (qPCR assay) and 24 h (microarray and qPCR assays).

RESULTS

The fatty acid composition of head kidney leukocytes varied between the two dietary groups (e.g. higher 20:5n-3 in the FO7 group). Following microarray assays using a 44K salmonid platform, Rank Products (RP) analysis showed 14 and 54 differentially expressed probes (DEP) (PFP < 0.05) between the two diets in control and pIC groups (FO5 vs. FO7), respectively. Nonetheless, Significance Analysis of Microarrays (SAM, FDR < 0.05) identified only one DEP between pIC groups of the two diets. Moreover, we identified a large number (i.e. 890 DEP in FO7 and 1128 DEP in FO5 overlapping between SAM and RP) of pIC-responsive transcripts, and several of them were involved in TLR-/RLR-dependent and cytokine-mediated pathways. The microarray results were validated as significantly differentially expressed by qPCR assays for 2 out of 9 diet-responsive transcripts and for all of the 35 selected pIC-responsive transcripts.

CONCLUSION

Fatty acid-binding protein adipocyte (fabp4) and proteasome subunit beta type-8 (psmb8) were significantly up- and down-regulated, respectively, in the MLCs of fish fed the diet with a lower level of fish oil, suggesting that they are important diet-responsive, immune-related biomarkers for future studies. Although the different levels of dietary fish and vegetable oils involved in this study affected the expression of some transcripts, the immune-related pathways and functions activated by the antiviral response of salmon MLCs in both groups were comparable overall. Moreover, the qPCR revealed transcripts responding early to pIC (e.g. lgp2, map3k8, socs1, dusp5 and cflar) and time-responsive transcripts (e.g. scarb1-a, csf1r, traf5a, cd80 and ctsf) in salmon MLCs. The present study provides a comprehensive picture of the putative molecular pathways (e.g. RLR-, TLR-, MAPK- and IFN-associated pathways) activated by the antiviral response of salmon MLCs.

Functional 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.

Comparing the transcriptomes of embryos from domesticated and wild Atlantic salmon (Salmo salar L.) stocks and examining factors that influence heritability of gene expression

BACKGROUND

Due to selective breeding, domesticated and wild Atlantic salmon are genetically diverged, which raises concerns about farmed escapees having the potential to alter the genetic composition of wild populations and thereby disrupting local adaptation. Documenting transcriptional differences between wild and domesticated stocks under controlled conditions is one way to explore the consequences of domestication and selection. We compared the transcriptomes of wild and domesticated Atlantic salmon embryos, by using a custom 44k oligonucleotide microarray to identify perturbed gene pathways between the two stocks, and to document the inheritance patterns of differentially-expressed genes by examining gene expression in their reciprocal hybrids.

RESULTS

Data from 24 array interrogations were analysed: four reciprocal cross types (W♀ × W♂, D♀ × W♂; W♀ × D♂, D♀ × D♂) × six biological replicates. A common set of 31,491 features on the microarrays passed quality control, of which about 62 % were assigned a KEGG Orthology number. A total of 6037 distinct genes were identified for gene-set enrichment/pathway analysis. The most highly enriched functional groups that were perturbed between the two stocks were cellular signalling and immune system, ribosome and RNA transport, and focal adhesion and gap junction pathways, relating to cell communication and cell adhesion molecules. Most transcripts that were differentially expressed between the stocks were governed by additive gene interaction (33 to 42 %). Maternal dominance and over-dominance were also prevalent modes of inheritance, with no convincing evidence for a stock effect.

CONCLUSIONS

Our data indicate that even at this relatively early developmental stage, transcriptional differences exist between the two stocks and affect pathways that are relevant to wild versus domesticated environments. Many of the identified differentially perturbed pathways are involved in organogenesis, which is expected to be an active process at the eyed egg stage. The dominant effects are more largely due to the maternal line than to the origin of the stock. This finding is particularly relevant in the context of potential introgression between farmed and wild fish, since female escapees tend to have a higher spawning success rate compared to males.

Infectious hematopoietic necrosis virus (IHNV) persistence in Sockeye Salmon: influence on brain transcriptome and subsequent response to the viral mimic poly(I:C)

Background

Sockeye Salmon are an iconic species widely distributed throughout the North Pacific. A devastating pathogen of Sockeye Salmon is infectious hematopoietic necrosis virus (IHNV, genus Novirhabdovirus, family Rhabdoviridae). It has been postulated that IHNV is maintained in salmon populations by persisting over the life of its host and/or by residing in natural reservoirs other than its susceptible hosts. Herein we demonstrate the presence of IHNV in the brain of Sockeye Salmon that survived an experimentally-induced outbreak, suggesting the presence of viral persistence in this susceptible species. To understand the viral persistent state in Sockeye Salmon we profiled the transcriptome to evaluate the host response in asymptomatic carriers and to determine what effects (if any) IHNV exposure may have on subsequent virus challenges.

Results

A laboratory disease model to simulate a natural IHNV outbreak in Sockeye Salmon resulted in over a third of the population incurring acute IHN disease and mortality during the first four months after initial exposure. Nine months post IHNV exposure, despite the absence of disease and mortality, a small percentage (<4 %) of the surviving population contained IHNV in brain. Transcriptome analysis in brain of asymptomatic virus carriers and survivors without virus exhibited distinct transcriptional profiles in comparison to naïve fish. Characteristic for carriers was the up-regulation of genes involved in antibody production and antigen presentation. In both carriers and survivors a down-regulation of genes related to cholesterol biosynthesis, resembling an antiviral mechanism observed in higher vertebrates was revealed along with differences in nervous system development. Moreover, following challenge with poly(I:C), survivors and carriers displayed an elevated antiviral immune response in comparison to naïve fish.

Conclusions

IHN virus persistence was identified in Sockeye Salmon where it elicited a unique brain transcriptome profile suggesting an ongoing adaptive immune response. IHNV carriers remained uncompromised in mounting efficient innate antiviral responses when exposed to a viral mimic. The capacity of IHNV to reside in asymptomatic hosts supports a virus carrier hypothesis and if proven infectious, could have significant epidemiological consequences towards maintaining and spreading IHNV among susceptible host populations.

Electronic supplementary material

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

The developmental transcriptome of contrasting Arctic charr (Salvelinus alpinus) morphs

Species and populations with parallel evolution of specific traits can help illuminate how predictable adaptations and divergence are at the molecular and developmental level. Following the last glacial period, dwarfism and specialized bottom feeding morphology evolved rapidly in several landlocked Arctic charrSalvelinus alpinuspopulations in Iceland.  To study the genetic divergence between small benthic morphs and limnetic morphs, we conducted RNA-sequencing charr embryos at four stages in early development. We studied two stocks with contrasting morphologies: the small benthic (SB) charr from Lake Thingvallavatn and Holar aquaculture (AC) charr.The data reveal significant differences in expression of several biological pathways during charr development. There was also an expression difference between SB- and AC-charr in genes involved in energy metabolism and blood coagulation genes. We confirmed differing expression of five genes in whole embryos with qPCR, includinglysozymeandnatterin-likewhich was previously identified as a fish-toxin of a lectin family that may be a putative immunopeptide. We also verified differential expression of 7 genes in the developing head that associated consistently with benthic v.s.limnetic morphology (studied in 4 morphs). Comparison of single nucleotide polymorphism (SNP) frequencies reveals extensive genetic differentiation between the SB and AC-charr (~1300 with more than 50% frequency difference). Curiously, three derived alleles in the otherwise conserved 12s and 16s mitochondrial ribosomal RNA genes are found in benthic charr.The data implicate multiple genes and molecular pathways in divergence of small benthic charr and/or the response of aquaculture charr to domestication. Functional, genetic and population genetic studies on more freshwater and anadromous populations are needed to confirm the specific loci and mutations relating to specific ecological traits in Arctic charr.

The developmental transcriptome of contrasting Arctic charr ( Salvelinus alpinus) morphs

Species and populations with parallel evolution of specific traits can help illuminate how predictable adaptations and divergence are at the molecular and developmental level. Following the last glacial period, dwarfism and specialized bottom feeding morphology evolved rapidly in several landlocked Arctic charr Salvelinus alpinus populations in Iceland.  

To study the genetic divergence between small benthic morphs and limnetic morphs, we conducted RNA-sequencing charr embryos at four stages in early development. We studied two stocks with contrasting morphologies: the small benthic (SB) charr from Lake Thingvallavatn and Holar aquaculture (AC) charr.

The data reveal significant differences in expression of several biological pathways during charr development. There was also an expression difference between SB- and AC-charr in genes involved in energy metabolism and blood coagulation genes. We confirmed differing expression of five genes in whole embryos with qPCR, including lysozyme and natterin-like which was previously identified as a fish-toxin of a lectin family that may be a putative immunopeptide. We also verified differential expression of 7 genes in the developing head that associated consistently with benthic v.s.limnetic morphology (studied in 4 morphs). Comparison of single nucleotide polymorphism (SNP) frequencies reveals extensive genetic differentiation between the SB and AC-charr (~1300 with more than 50% frequency difference). Curiously, three derived alleles in the otherwise conserved 12s and 16s mitochondrial ribosomal RNA genes are found in benthic charr.

The data implicate multiple genes and molecular pathways in divergence of small benthic charr and/or the response of aquaculture charr to domestication. Functional, genetic and population genetic studies on more freshwater and anadromous populations are needed to confirm the specific loci and mutations relating to specific ecological traits in Arctic charr.

Ontogeny of the Digestive System of Atlantic Salmon (Salmo salar L.) and Effects of Soybean Meal from Start-Feeding

Despite a long history of rearing Atlantic salmon in hatcheries in Norway, knowledge of molecular and physiological aspects of juvenile development is still limited. To facilitate introduction of alternative feed ingredients and feed additives during early phases, increased knowledge regarding the ontogeny of the digestive apparatus in salmon is needed. In this study, we characterized the development of the gastrointestinal tract and accessory digestive organs for five months following hatch by using histological, biochemical and molecular methods. Furthermore, the effects of a diet containing 16.7% soybean meal (SBM) introduced at start-feeding were investigated, as compared to a fishmeal based control diet. Salmon yolk sac alevins and fry were sampled at 18 time points from hatch until 144 days post hatch (dph). Histomorphological development was investigated at 7, 27, 46, 54 and 144 dph. Ontogenetic expression patterns of genes encoding key digestive enzymes, nutrient transporters, gastrointestinal peptide hormones and T-cell markers were analyzed from 13 time points by qPCR. At 7 dph, the digestive system of Atlantic salmon alevins was morphologically distinct with an early stomach, liver, pancreas, anterior and posterior intestine. About one week before the yolk sac was internalized and exogenous feeding was started, gastric glands and developing pyloric caeca were observed, which coincided with an increase in gene expression of gastric and pancreatic enzymes and nutrient transporters. Thus, the observed organs seemed ready to digest external feed well before the yolk sac was absorbed into the abdominal cavity. In contrast to post-smolt Atlantic salmon, inclusion of SBM did not induce intestinal inflammation in the juveniles. This indicates that SBM can be used in compound feeds for salmon fry from start-feeding to at least 144 dph and/or 4-5 g body weight.

A comparison between egg trancriptomes of cod and salmon reveals species-specific traits in eggs for each species

Fish in use in aquaculture display large variation in gamete biology. To reach better understanding around this issue, this study aims at identifying if species specific "egg life history traits" can be hidden in the unfertilized egg. This was done by investigating egg transcriptome differences between Atlantic salmon and Atlantic cod. Salmon and cod eggs were selected due to their largely differencing phenotypes. An oligo microarray analysis was performed on ovulated eggs from cod (n = 8) and salmon (n = 7). The arrays were normalized to a similar spectrum for both arrays. Both arrays were re-annotated with SWISS-Prot and KEGG genes to retrieve an official gene symbol and an orthologous KEGG annotation, in salmon and cod arrays this represented 14,009 and 7,437 genes respectively. The probe linked to the highest gene expression for that particular KEGG annotation was used to compare expression between species. Differential expression was calculated for genes that had an annotation with score >300, resulting in a total of 2,457 KEGG annotations (genes) being differently expressed between the species (FD > 2). This analysis revealed that immune, signal transduction and excretory related pathways were overrepresented in salmon compared to cod. The most overrepresented pathways in cod were related to regulation of genetic information processing and metabolism. To conclude this analysis clearly point at some distinct transcriptome repertoires for cod and salmon and that these differences may explain some of the species-specific biological features for salmon and cod eggs.

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.

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

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

Altered gene expression patterns of innate and adaptive immunity pathways in transgenic rainbow trout harboring Cecropin P1 transgene

Background

We have recently developed several homozygous families of transgenic rainbow trout harbouring cecropin P1 transgene. These fish exhibit resistance characteristic to infection by Aeromonas salmonicida and infectious hematopoietic necrosis virus (IHNV). In our earlier studies we have reported that treatment of a rainbow trout macrophage cell line (RTS11) with a linear cationic α-helical antimicrobial peptide (e.g., cecropin B) resulted in elevated levels of expression of two pro-inflammatory relevant genes (e.g., IL-1β and COX-2). Therefore, we hypothesized that in addition to the direct antimicrobial activity of cecropin P1 in the disease resistant transgenic rainbow trout, this antimicrobial peptide may also affect the expression of immune relevant genes in the host. To confirm this hypothesis, we launched a study to determine the global gene expression profiles in three immune competent organs of cecropin P1 transgenic rainbow trout by using a 44k salmonid microarray.

Results

From the microarray data, a total of 2480 genes in the spleen, 3022 in the kidney, and 2102 in the liver were determined as differentially expressed genes (DEGs) in the cecropin P1 transgenic rainbow trout when compared to the non-transgenics. There were 478 DEGs in common among three tissues. Enrichment analyses conducted by two different bioinformatics tools revealed a tissue specific profile of functional pathway perturbation. Many of them were directly related to innate immune system such as phagocytosis, lysosomal processing, complement activation, antigen processing/presentation, and leukocyte migration. Perturbation of other biological functions that might contribute indirectly to host immunity was also observed.

Conclusions

The gene product of cecropin P1 transgene produced in the disease resistant transgenic rainbow trout not only can kill the pathogens directly but also exert multifaceted immunomodulatory properties to boost host immunity. The identified genes involved in different pathways related to immune function are valuable indicators associated with enhanced host immunity. These genes may serve as markers for selective breeding of rainbow trout or other aquaculture important fish species bearing traits of disease resistance.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-887) contains supplementary material, which is available to authorized users.

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

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

Infectious disease, shifting climates, and opportunistic predators: cumulative factors potentially impacting wild salmon declines

Emerging diseases are impacting animals under high-density culture, yet few studies assess their importance to wild populations. Microparasites selected for enhanced virulence in culture settings should be less successful maintaining infectivity in wild populations, as once the host dies, there are limited opportunities to infect new individuals. Instead, moderately virulent microparasites persisting for long periods across multiple environments are of greatest concern. Evolved resistance to endemic microparasites may reduce susceptibilities, but as barriers to microparasite distributions are weakened, and environments become more stressful, unexposed populations may be impacted and pathogenicity enhanced. We provide an overview of the evolutionary and ecological impacts of infectious diseases in wild salmon and suggest ways in which modern technologies can elucidate the microparasites of greatest potential import. We present four case studies that resolve microparasite impacts on adult salmon migration success, impact of river warming on microparasite replication, and infection status on susceptibility to predation. Future health of wild salmon must be considered in a holistic context that includes the cumulative or synergistic impacts of multiple stressors. These approaches will identify populations at greatest risk, critically needed to manage and potentially ameliorate the shifts in current or future trajectories of wild populations.

Divergent immunity and energetic programs in the gills of migratory and resident Oncorhynchus mykiss

Divergent life history strategies occur in steelhead or rainbow trout Oncorhynchus mykiss, and many populations produce both migrant (anadromous fish that move to the ocean after rearing) and resident (do not migrate and remain in fresh water) individuals. Mechanisms leading to each type are only partially understood; while the general tendency of a population is heritable, individual tendency may be plastic, influenced by local environment. Steelhead hatchery programmes aim to mitigate losses in wild stocks by producing trout that will migrate to the ocean and not compete with wild trout for limited freshwater resources. To increase our understanding of gill function in these migratory or resident phenotypes, here we compare gill transcriptome profiles of hatchery-released fish either at the release site (residents) or five river kilometres downstream while still in full fresh water (migrants). To test whether any of these genes can be used as predictive markers for smoltification, we compared these genes between migrant-like and undifferentiated trout while still in the hatchery in a common environment (prerelease). Results confirmed the gradual process of smoltification, and the importance of energetics, gill remodelling and ion transport capacity for migrants. Additionally, residents overexpressed transcripts involved in antiviral defences, potentially for immune surveillance via dendritic cells in the gills. The best smoltification marker candidate was protein s100a4, expression of which was highly correlated with Na(+) , K(+) ATPase (NKA) activity and smolt-like morphology in pre- and postrelease trout gills.

Comparative transcriptomics of Atlantic Salmo salar, chum Oncorhynchus keta and pink salmon O. gorbuscha during infections with salmon lice Lepeophtheirus salmonis

Background

Salmon species vary in susceptibility to infections with the salmon louse (Lepeophtheirus salmonis). Comparing mechanisms underlying responses in susceptible and resistant species is important for estimating impacts of infections on wild salmon, selective breeding of farmed salmon, and expanding our knowledge of fish immune responses to ectoparasites. Herein we report three L. salmonis experimental infection trials of co-habited Atlantic Salmo salar, chum Oncorhynchus keta and pink salmon O. gorbuscha, profiling hematocrit, blood cortisol concentrations, and transcriptomic responses of the anterior kidney and skin to the infection.

Results

In all trials, infection densities (lice per host weight (g)) were consistently highest on chum salmon, followed by Atlantic salmon, and lowest in pink salmon. At 43 days post-exposure, all lice had developed to motile stages, and infection density was uniformly low among species. Hematocrit was reduced in infected Atlantic and chum salmon, and cortisol was elevated in infected chum salmon. Systemic transcriptomic responses were profiled in all species and large differences in response functions were identified between Atlantic and Pacific (chum and pink) salmon. Pink and chum salmon up-regulated acute phase response genes, including complement and coagulation components, and down-regulated antiviral immune genes. The pink salmon response involved the largest and most diverse iron sequestration and homeostasis mechanisms. Pattern recognition receptors were up-regulated in all species but the active components were often species-specific. C-type lectin domain family 4 member M and acidic mammalian chitinase were specifically up-regulated in the resistant pink salmon.

Conclusions

Experimental exposures consistently indicated increased susceptibility in chum and Atlantic salmon, and resistance in pink salmon, with differences in infection density occurring within the first three days of infection. Transcriptomic analysis suggested candidate resistance functions including local inflammation with cytokines, specific innate pattern recognition receptors, and iron homeostasis. Suppressed antiviral immunity in both susceptible and resistant species indicates the importance of future work investigating co-infections of viral pathogens and lice.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-200) contains supplementary material, which is available to authorized users.

Influence of developmental stage and genotype on liver mRNA levels among wild, domesticated, and hybrid rainbow trout (Oncorhynchus mykiss)

BACKGROUND

Release of domesticated strains of fish into nature may pose a threat to wild populations with respect to their evolved genetic structure and fitness. Understanding alterations that have occurred in both physiology and genetics as a consequence of domestication can assist in evaluating the risks posed by introgression of domesticated genomes into wild genetic backgrounds, however the molecular causes of these consequences are currently poorly defined. The present study has examined levels of mRNA in fast-growing pure domesticated (D), slow-growing age-matched pure wild (Wa), slow-growing size-matched pure wild (Ws), and first generation hybrid cross (W/D) rainbow trout (Oncorhynchus mykiss) to investigate the influence of genotype (domesticated vs. wild, and their interactions in hybrids) and developmental stage (age- or size-matched animals) on genetic responses (i.e. dominant vs. recessive) and specific physiological pathways.

RESULTS

Highly significant differences in mRNA levels were found between domesticated and wild-type rainbow trout genotypes (321 mRNAs), with many mRNAs in the wild-domesticated hybrid progeny showing intermediate levels. Differences were also found between age-matched and size-matched wild-type trout groups (64 mRNAs), with unique mRNA differences for each of the wild-type groups when compared to domesticated trout (Wa: 114 mRNAs, Ws: 88 mRNAs), illustrating an influence of fish developmental stage affecting findings when used as comparator groups to other genotypes. Analysis of differentially expressed mRNAs (found for both wild-type trout to domesticated comparisons) among the genotypes indicates that 34.8% are regulated consistent with an additive genetic model, whereas 39.1% and 26.1% show a recessive or dominant mode of regulation, respectively. These molecular data are largely consistent with phenotypic data (growth and behavioural assessments) assessed in domesticated and wild trout strains.

CONCLUSIONS

The present molecular data are concordant with domestication having clearly altered rainbow trout genomes and consequent phenotype from that of native wild populations. Although mainly additive responses were noted in hybrid progeny, the prevalence of dominant and non-additive responses reveals that introgression of domesticated and wild genotypes alters the type of genetic control of mRNA levels from that of wild-type, which may lead to disruption of gene regulation systems important for developing phenotypes for optimal fitness in nature. A clear influence of both fish age and size (developmental stage) on mRNA levels was also noted in this study, which highlights the importance of examining multiple control samples to provide a comprehensive understanding of changes observed between strains possessing differences in growth rate.

Early response of gene expression in the distal intestine of Atlantic salmon (Salmo salar L.) during the development of soybean meal induced enteritis

Plant products in general and soybeans in particular can challenge the function and health of the intestinal tract. Salmonids develop an intestinal inflammation when fed diets containing soybean meal (SBM) and certain other legume ingredients. In the present study a 44K oligonucleotide salmonid microarray, qPCR and histology were used to investigate early response mechanisms in the distal intestine of Atlantic salmon (Salmo salar L.) during the first week of oral exposure to a diet containing 20% extracted SBM. The distal intestine transcriptome was profiled on days 1, 2, 3, 5 and 7 and compared to a control group fed fishmeal as the sole protein source. Histological evaluation of the distal intestine revealed the first signs of inflammation on day 5. The most prominent gene expression changes were seen on days 3 and 5. Up-regulation in immune-related genes was observed during the first 5 days, including GTPase IMAP family members, NF-kB-related genes and regulators of T cell and B cell function. Many functional genes involved in lipid metabolism, proteolysis, transport, metabolism and detoxification were initially up-regulated on days 1-3, possibly as an attempt by the tissue to compensate for the initiating immune response. Cell repair and extracellular matrix remodeling genes were up-regulated (heparanase, collagenase) on days 3 and 5. Down regulation of genes related to endocytosis, exocytosis, detoxification, transporters and metabolic processes from day 5 indicated initiation of dysfunction of digestive and metabolic functions that may occur as a result of inflammation or as a response to the introduction of soybean meal in the diet. This is the first study conducting transcriptomic profiling to characterize early responses during the development of SBMIE. Switching Atlantic salmon from a fishmeal to a 20% SBM diet resulted in rapid changes to the intestinal transcriptome, indicating an immune reaction with subsequent impaired epithelial barrier function and other vital intestinal functions.

Iron-sulfur cluster scaffold (ISCU) gene is duplicated in salmonid fish and tissue and temperature dependent expressed in rainbow trout

The iron-sulfur cluster protein ISCU is a scaffold protein tasked with the building and mediation of iron-sulfur [Fe-S]-clusters. These are crucial for [Fe-S]-enzymes, which are involved in essential biological cell processes like metabolism or ion transport. Analysis of ISCU in rainbow trout (Oncorhynchus mykiss) and maraena whitefish (Coregonus maraena) revealed the existence of two gene variants in each of the two salmonids. This study presents the characterization of the duplicated ISCU cDNA sequences in both species as well as the comparative functional analysis of the genes in healthy and affected fish of two rainbow trout strains differing in trait robustness under regional aquaculture conditions. Coding sequences of trout ISCUA and ISCUB genes are spanning over five exons. Open reading frames (ORF) of trout (ISCUA: 495bp, ISCUB: 498bp) and whitefish (ISCUA and ISCUB: 495bp) genes encode for evolutionary highly conserved proteins and share 72% sequence similarity with human ISCU. Transcriptome analyses comparing healthy fish of the local rainbow trout strain BORN and the import strain TCO revealed strain-specific expression patterns for ISCU. Expression analyses by quantitative RT-PCR indicated remarkable differences between the transcript level of the gene variants ISCUA and ISCUB. Moderate temperature challenge (8°C and 23°C) suggests a generally higher transcript level of the two gene variants at 8°C in the liver, spleen, and gill of both strains. However, no remarkable differences between the strains occurred in the temperature-dependent ISCU gene expression profiles. The experimental infection with Aeromonas salmonicida resulted in a different ISCU gene expression in the gill and trunk kidney of both strains after two weeks, suggesting a specific role of the scaffold gene in rainbow trout strain BORN, regarding the recovery after infection. Although results partially reflect the expected strain- and tissue-specific ISCUA and ISCUB regulation in rainbow trout, the data do not support the assumed association of ISCU with the trait robustness.

A genomic and ecotoxicological perspective of DNA array studies in aquatic environmental risk assessment

Ecotoxicogenomics is developing into a key tool for the assessment of environmental impacts and environmental risk assessment for aquatic ecosystems. This review aims to report achievements and drawbacks of this technique and to explore potential conceptual and experimental procedures to improve future investigations. Ecotoxicogenomic literature evidences the ability of genomic technologies to characterize toxicant specific gene transcriptome patterns that can be used to identify major toxicants affecting aquatic species. They also contribute decisively to the development of new molecular biomarkers and, in many cases, to the determination of new possible gene targets. Primary transcriptomic responses obtained after short exposures provided more information of putative gene targets than secondary responses obtained after long, chronic exposures, which in turn are usually more accurate to describe actual environmental impacts in natural populations. Several problems need to be addressed in future investigations: the lack of studies (and genomic information) on key ecological species and taxa, the need to better understand the different transcriptomic responses to high and low doses and, especially, short and long exposures, and the need to improve experimental designs to minimize false transcriptome interpretations of target genes.