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

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

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

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

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

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

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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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.

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.

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

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

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

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

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

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

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

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

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

BACKGROUND: Atlantic cod (Gadus morhua) reared in sea-cages can experience large variations in temperature, and these have been shown to affect their immune function. We used the new 20K Atlantic cod microarray to investigate how a water temperature change which, simulates that seen in Newfoundland during the spring-summer (i.e. from 10°C to 16°C, 1°C increase every 5 days) impacted the cod spleen transcriptome response to the intraperitoneal injection of a viral mimic (polyriboinosinic polyribocytidylic acid, pIC). RESULTS: The temperature regime alone did not cause any significant increases in plasma cortisol levels and only minor changes in spleen gene transcription. However, it had a considerable impact on the fish spleen transcriptome response to pIC [290 and 339 significantly differentially expressed genes between 16°C and 10°C at 6 and 24 hours post-injection (HPI), respectively]. Seventeen microarray-identified transcripts were selected for QPCR validation based on immune-relevant functional annotations. Fifteen of these transcripts (i.e. 88%), including DHX58, STAT1, IRF7, ISG15, RSAD2 and IκBα, were shown by QPCR to be significantly induced by pIC. CONCLUSIONS: The temperature increase appeared to accelerate the spleen immune transcriptome response to pIC. We found 41 and 999 genes differentially expressed between fish injected with PBS vs. pIC at 10°C and sampled at 6HPI and 24HPI, respectively. In contrast, there were 656 and 246 genes differentially expressed between fish injected with PBS vs. pIC at 16°C and sampled at 6HPI and 24HPI, respectively. Our results indicate that the modulation of mRNA expression of genes belonging to the NF-κB and type I interferon signal transduction pathways may play a role in controlling temperature-induced changes in the spleen's transcript expression response to pIC. Moreover, interferon effector genes such as ISG15 and RSAD2 were differentially expressed between fish injected with pIC at 10°C vs. 16°C at 6HPI. These results substantially increase our understanding of the genes and molecular pathways involved in the negative impacts of elevated ambient temperature on fish health, and may also be valuable to our understanding of how accelerated global climate change could impact cold-water marine finfish species.

Dynamic expression profiles of virus-responsive and putative antimicrobial peptide-encoding transcripts during Atlantic cod (Gadus morhua) embryonic and early larval development.

Early life stage mortality is one of the problems faced by Atlantic cod aquaculture. However, our understanding of immunity in early life stage fish is still incomplete, and the information available is restricted to a few species. In the present work we investigated the expression of immune-relevant transcripts in Atlantic cod during early development. The transcripts subjected to QPCR analysis in the present study were previously identified as putative anti-viral or anti-bacterial genes in Atlantic cod using suppression subtractive hybridization (SSH) libraries, QPCR, and/or microarrays. Of the 11 genes involved in this study, only atf3, cxc chemokine and gaduscidin-1 were not detected at the transcript level in all developmental stages investigated from unfertilized egg to early larval stage. Adam22, hamp, il8, irf1, irf7, lgp2, sacsin, and stat1 transcripts were detected in unfertilized egg and 7h post-fertilization (~2-cell stage) embryos, showing maternal contribution of these immune-relevant transcripts to the early embryonic transcriptome. The Atlantic cod genes included in this study presented diverse transcript expression profiles throughout embryonic and early larval development. For example, adam22 and sacsin transcripts rose abruptly during blastula/gastrula stage and were then expressed at relatively high levels through subsequent embryonic and early larval developmental stages. A peak in irf1 and irf7 transcript expression during early segmentation suggests that these interferon pathway genes play developmental stage-specific roles during cod embryogenesis. Stat1 had increasing transcript expression throughout blastula/gastrula, segmentation, and early larval developmental stages. Atf3, cxc chemokine, gaduscidin-1, and il8 transcripts rose approximately 2-3 fold during hatching, supporting the hypothesis that there is preparation at the immune-relevant transcript expression level to deal with environmental pathogens that may be encountered during early larval development. The specific roles that interferon pathway and other immune-relevant genes play in early life stage cod, and the potential impact of their dynamic transcript expression on immune competence of Atlantic cod embryos and larvae, remain unclear and warrant further study.