The prokaryotic and eukaryotic microbiome of Pacific oyster spat is shaped by ocean warming but not acidification

Pacific oysters (Magallana gigas, a.k.a. Crassostrea gigas), the most widely farmed oysters, are under threat from climate change and emerging pathogens. In part, their resilience may be affected by their microbiome, which, in turn, may be influenced by ocean warming and acidification. To understand these impacts, we exposed early-development Pacific oyster spat to different temperatures (18°C and 24°C) and pCO2 levels (800, 1,600, and 2,800 µatm) in a fully crossed design for 3 weeks. Under all conditions, the microbiome changed over time, with a large decrease in the relative abundance of potentially pathogenic ciliates (Uronema marinum) in all treatments with time. The microbiome composition differed significantly with temperature, but not acidification, indicating that Pacific oyster spat microbiomes can be altered by ocean warming but is resilient to ocean acidification in our experiments. Microbial taxa differed in relative abundance with temperature, implying different adaptive strategies and ecological specializations among microorganisms. Additionally, a small proportion (~0.2% of the total taxa) of the relatively abundant microbial taxa were core constituents (>50% occurrence among samples) across different temperatures, pCO2 levels, or time. Some taxa, including A4b bacteria and members of the family Saprospiraceae in the phyla Chloroflexi (syn. Chloroflexota) and Bacteroidetes (syn. Bacteroidota), respectively, as well as protists in the genera Labyrinthula and Aplanochytrium in the class Labyrinthulomycetes, and Pseudoperkinsus tapetis in the class Ichthyosporea were core constituents across temperatures, pCO2 levels, and time, suggesting that they play an important, albeit unknown, role in maintaining the structural and functional stability of the Pacific oyster spat microbiome in response to ocean warming and acidification. These findings highlight the flexibility of the spat microbiome to environmental changes.IMPORTANCEPacific oysters are the most economically important and widely farmed species of oyster, and their production depends on healthy oyster spat. In turn, spat health and productivity are affected by the associated microbiota; yet, studies have not scrutinized the effects of temperature and pCO2 on the prokaryotic and eukaryotic microbiomes of spat. Here, we show that both the prokaryotic and, for the first time, eukaryotic microbiome of Pacific oyster spat are surprisingly resilient to changes in acidification, but sensitive to ocean warming. The findings have potential implications for oyster survival amid climate change and underscore the need to understand temperature and pCO2 effects on the microbiome and the cascading effects on oyster health and productivity.

Population genomics of harbour seal Phoca vitulina from northern British Columbia through California and comparison to the Atlantic subspecies

The harbour seal Phoca vitulina is a ubiquitous pinniped species found throughout coastal waters of the Northern Hemisphere. Harbour seal impacts on ecosystem dynamics may be significant due to their high abundance and food web position. Two subspecies exist in North America, P. v. richardii in the Pacific Ocean and P. v. vitulina in the Atlantic. Strong natal philopatry of harbour seals can result in fine-scale genetic structure and isolation by distance. Management of harbour seals is expected to benefit from improved resolution of seal population structure and dynamics. Here, we use genotyping-by-sequencing to genotype 146 harbour seals from the eastern Pacific Ocean (i.e. British Columbia (BC), Oregon and California) and the western Atlantic Ocean (i.e. Québec, Newfoundland and Labrador). Using 12,742 identified variants, we confirm the recently identified elevated genetic diversity in the eastern Pacific relative to the western Atlantic and greatest differentiation between the subspecies. Further, we demonstrate that this is independent of reference genome bias or other potential technical artefacts. Coast-specific analyses with 8933 and 3828 variants in Pacific and Atlantic subspecies, respectively, identify divergence between BC and Oregon-California, and between Québec and Newfoundland-Labrador. Unexpected PCA outlier clusters were observed in two populations due to cryptic relatedness of individuals; subsequently, closely related samples were removed. Admixture analysis indicates an isolation-by-distance signature where Oregon seals contained some of the BC signature, whereas California did not. Additional sampling is needed in the central and north coast of BC to determine whether a discrete separation of populations exists within the region.

Genomic diversity of wild and cultured Yesso scallop Mizuhopecten yessoensis from Japan and Canada

The Yesso scallop Mizuhopecten yessoensis is an important aquaculture species that was introduced to Western Canada from Japan to establish an economically viable scallop farming industry. This highly fecund species has been propagated in Canadian aquaculture hatcheries for the past 40 years, raising questions about genetic diversity and genetic differences among hatchery stocks. In this study, we compare cultured Canadian and wild Japanese populations of Yesso scallop using double-digest restriction site-associated DNA (ddRAD) sequencing to genotype 21,048 variants in 71 wild-caught scallops from Japan, 65 scallops from the Vancouver Island University breeding population, and 37 scallops obtained from a commercial farm off Vancouver Island, British Columbia. The wild scallops are largely comprised of equally unrelated individuals, whereas cultured scallops are comprised of multiple families of related individuals. The polymorphism rate estimated in wild scallops was 1.7%, whereas in the cultured strains, it ranged between 1.35 and 1.07%. Interestingly, heterozygosity rates were highest in the cultured populations, which is likely due to shellfish hatchery practices of crossing divergent strains to gain benefits of heterosis and to avoid inbreeding. Evidence of founder effects and drift was observed in the cultured strains, including high genetic differentiation between cultured populations and between cultured populations and the wild population. Cultured populations had effective population sizes ranging from 9 to 26 individuals whereas the wild population was estimated at 25,048-56,291 individuals. Further, a depletion of low-frequency variants was observed in the cultured populations. These results indicate significant genetic diversity losses in cultured scallops in Canadian breeding programs.

Chinook salmon diversity contributes to fishery stability and trade-offs with mixed-stock harvest

Variation among populations in life history and intrinsic population characteristics (i.e., population diversity) helps maintain resilience to environmental change and dampen interannual variability in ecosystem services. As a result, ecological variation, and the processes that generate it, is considered central to strategies for managing risks to ecosystems in an increasingly variable and uncertain world. However, characterizing population diversity is difficult, particularly in large and remote regions, which often prevents its formal consideration in management advice. We combined genetic stock identification of archived scale and tissue samples with state-space run-reconstruction models to estimate migration timing and annual return abundance for eight geographically and genetically distinct Chinook salmon populations within the Canadian portion of the Yukon River. We found that among-population variation in migration timing and return abundances resulted in aggregate return migrations that were 2.1 times longer and 1.4 times more stable than if they had composed a single homogeneous population. We then fit state-space spawner-recruitment models to the annual return abundances to characterize among-population diversity in intrinsic productivity and population size and their consequences for the fisheries they support. Productivity and carrying capacity varied among populations by approximately 2.4-fold (2.9 to 6.9 recruits per spawner) and three-fold (8800 to 27,000 spawners), respectively. This diversity implies an equilibrium trade-off between harvesting of the population aggregate and the conservation of individual populations whereby the harvest rate predicted to maximize aggregate harvests comes at the cost of overfishing ~40% of the populations but with a relatively low risk of extirpating the weakest ones. Our findings illustrate how population diversity in one of the largest salmon-producing river basins in the world contributes to fishery stability and food security in a region where salmon have high cultural and subsistence value. More generally, our work demonstrates the utility of molecular analyses of archived biological material for characterizing diversity in biological systems and its benefits and consequences for trade-offs in decision-making.

In-field genetic stock identification of overwintering coho salmon in the Gulf of Alaska: Evaluation of Nanopore sequencing for remote real-time deployment

Genetic stock identification (GSI) from genotyping‐by‐sequencing of single nucleotide polymorphism (SNP) loci has become the gold standard for stock of origin identification in Pacific salmon. The sequencing platforms currently applied require large batch sizes and multiday processing in specialized facilities to perform genotyping by the thousands. However, recent advances in third‐generation single‐molecule sequencing platforms, such as the Oxford Nanopore minION, provide base calling on portable, pocket‐sized sequencers and promise real‐time, in‐field stock identification of variable batch sizes. Here we evaluate utility and comparability to established GSI platforms of at‐sea stock identification of coho salmon (Oncorhynchus kisutch) using targeted SNP amplicon sequencing on the minION platform during a high‐sea winter expedition to the Gulf of Alaska. As long read sequencers are not optimized for short amplicons, we concatenate amplicons to increase coverage and throughput. Nanopore sequencing at‐sea yielded data sufficient for stock assignment for 50 out of 80 individuals. Nanopore‐based SNP calls agreed with Ion Torrent‐based genotypes in 83.25%, but assignment of individuals to stock of origin only agreed in 61.5% of individuals, highlighting inherent challenges of Nanopore sequencing, such as resolution of homopolymer tracts and indels. However, poor representation of assayed salmon in the queried baseline data set contributed to poor assignment confidence on both platforms. Future improvements will focus on lowering turnaround time and cost, increasing accuracy and throughput, as well as augmentation of the existing baselines. If successfully implemented, Nanopore sequencing will provide an alternative method to the large‐scale laboratory approach by providing mobile small batch genotyping to diverse stakeholders.

Population differences in Chinook salmon (Oncorhynchus tshawytscha) DNA methylation: Genetic drift and environmental factors

Local adaptation and phenotypic differences among populations have been reported in many species, though most studies focus on either neutral or adaptive genetic differentiation. With the discovery of DNA methylation, questions have arisen about its contribution to individual variation in and among natural populations. Previous studies have identified differences in methylation among populations of organisms, although most to date have been in plants and model animal species. Here we obtained eyed eggs from eight populations of Chinook salmon (Oncorhynchus tshawytscha) and assayed DNA methylation at 23 genes involved in development, immune function, stress response, and metabolism using a gene-targeted PCR-based assay for next-generation sequencing. Evidence for population differences in methylation was found at eight out of 23 gene loci after controlling for developmental timing in each individual. However, we found no correlation between freshwater environmental parameters and methylation variation among populations at those eight genes. A weak correlation was identified between pairwise DNA methylation dissimilarity among populations and pairwise F ST based on 15 microsatellite loci, indicating weak effects of genetic drift or geographic distance on methylation. The weak correlation was primarily driven by two genes, GTIIBS and Nkef. However, single-gene Mantel tests comparing methylation and pairwise F ST were not significant after Bonferroni correction. Thus, population differences in DNA methylation are more likely related to unmeasured oceanic environmental conditions, local adaptation, and/or genetic drift. DNA methylation is an additional mechanism that contributes to among population variation, with potential influences on organism phenotype, adaptive potential, and population resilience.

Relative genomic impacts of translocation history, hatchery practices, and farm selection in Pacific oyster Crassostrea gigas throughout the Northern Hemisphere

Pacific oyster Crassostrea gigas, endemic to coastal Asia, has been translocated globally throughout the past century, resulting in self-sustaining introduced populations (naturalized). Oyster aquaculture industries in many parts of the world depend on commercially available seed (hatchery-farmed) or naturalized/wild oysters to move onto a farm (naturalized-farmed). It is therefore important to understand genetic variation among populations and farm types. Here, we genotype naturalized/wild populations from France, Japan, China, and most extensively in coastal British Columbia, Canada. We also genotype cultured populations from throughout the Northern Hemisphere to compare with naturalized populations. In total, 16,942 markers were identified using double-digest RAD-sequencing in 182 naturalized, 112 hatchery-farmed, and 72 naturalized-farmed oysters (n = 366). Consistent with previous studies, very low genetic differentiation was observed around Vancouver Island (mean F _ST = 0.0019) and low differentiation between countries in the Japan-Canada-France historical translocation lineage (France-Canada F _ST = 0.0024; Japan-Canada F _ST = 0.0060). Chinese populations were more differentiated (China-Japan F _ST = 0.0241). Hatchery-propagated populations had higher interindividual relatedness suggesting family structure. Within-population inbreeding was not detected on farms, but nucleotide diversity and polymorphism rate were lower in one farm population. Moving oysters from nature onto farms did not result in strong within-generation selection. Private alleles at substantial frequency were identified in several hatchery populations grown in BC, suggesting nonlocal origins. Tests of selection identified outlier loci consistent with selective differences associated with domestication, in some cases consistently identified in multiple farms. Top outlier candidates were nearby genes involved in calcium signaling and calmodulin activity. Implications of potential introgression from hatchery-farmed oysters depend on whether naturalized populations are valued as a locally adapted resource or as an introduced, invasive species. Given the value of the industry in BC and the challenges the industry faces (e.g., climate change, crop losses, biotic stressors), this remains an important question.

Sex-Specific Co-expression Networks and Sex-Biased Gene Expression in the Salmonid Brook Charr Salvelinus fontinalis

Networks of co-expressed genes produce complex phenotypes associated with functional novelty. Sex differences in gene expression levels or in the structure of gene co-expression networks can cause sexual dimorphism and may resolve sexually antagonistic selection. Here we used RNA-sequencing in the salmonid Brook Charr Salvelinus fontinalis to characterize sex-specific co-expression networks in the liver of 47 female and 53 male offspring. In both networks, modules were characterized for functional enrichment, hub gene identification, and associations with 15 growth, reproduction, and stress-related phenotypes. Modules were then evaluated for preservation in the opposite sex, and in the congener Arctic Charr Salvelinus alpinus Overall, more transcripts were assigned to a module in the female network than in the male network, which coincided with higher inter-individual gene expression and phenotype variation in the females. Most modules were preserved between sexes and species, including those involved in conserved cellular processes (e.g., translation, immune pathways). However, two sex-specific male modules were identified, and these may contribute to sexual dimorphism. To compare with the network analysis, differentially expressed transcripts were identified between the sexes, revealing a total of 16% of expressed transcripts as sex-biased. For both sexes, there was no overrepresentation of sex-biased genes or sex-specific modules on the putative sex chromosome. Sex-biased transcripts were also not overrepresented in sex-specific modules, and in fact highly male-biased transcripts were enriched in preserved modules. Comparative network analysis and differential expression analyses identified different aspects of sex differences in gene expression, and both provided new insights on the genes underlying sexual dimorphism in the salmonid Brook Charr.

Effects of the vertically transmitted microsporidian Facilispora margolisi and the parasiticide emamectin benzoate on salmon lice (Lepeophtheirus salmonis)

Background

Microsporidia are highly specialized, parasitic fungi that infect a wide range of eukaryotic hosts from all major taxa. Infections cause a variety of damaging effects on host physiology from increased stress to death. The microsporidian Facilispora margolisi infects the Pacific salmon louse (Lepeophtheirus salmonis oncorhynchi), an economically and ecologically important ectoparasitic copepod that can impact wild and cultured salmonids.

Results

Vertical transmission of F. margolisi was demonstrated by using PCR and in situ hybridization to identify and localize microsporidia in female L. salmonis and their offspring. Spores and developmental structures of F. margolisi were identified in 77% of F₁ generation copepods derived from infected females while offspring from uninfected females all tested negative for the microsporidia. The transcriptomic response of the salmon louse to F. margolisi was profiled at both the copepodid larval stage and the pre-adult stage using microarray technology. Infected copepodids differentially expressed 577 transcripts related to stress, ATP generation and structural components of muscle and cuticle. The infection also impacted the response of the copepodid to the parasiticide emamectin benzoate (EMB) at a low dose of 1.0 ppb for 24 h. A set of 48 transcripts putatively involved in feeding and host immunomodulation were up to 8-fold underexpressed in the F. margolisi infected copepodids treated with EMB compared with controls or either stressor alone. Additionally, these infected lice treated with EMB also overexpressed 101 transcripts involved in stress resistance and signalling compared to the other groups. In contrast, infected pre-adult lice did not display a stress response, suggesting a decrease in microsporidian virulence associated with lice maturity. Furthermore, copepodid infectivity and moulting was not affected by the microsporidian infection.

Conclusions

This study demonstrated that F. margolisi is transmitted vertically between salmon louse generations and that biological impacts of infection differ depending on the stage of the copepod host. The infection caused significant perturbations of larval transcriptomes and therefore must be considered in future studies in which impacts to host development and environmental factors are assessed. Fitness impacts are probably minor, although the interaction between pesticide exposure and microsporidian infection merits further study.

Electronic supplementary material

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

Sex Chromosome Evolution, Heterochiasmy, and Physiological QTL in the Salmonid Brook Charr Salvelinus fontinalis

Whole-genome duplication (WGD) can have large impacts on genome evolution, and much remains unknown about these impacts. This includes the mechanisms of coping with a duplicated sex determination system and whether this has an impact on increasing the diversity of sex determination mechanisms. Other impacts include sexual conflict, where alleles having different optimums in each sex can result in sequestration of genes into nonrecombining sex chromosomes. Sex chromosome development itself may involve sex-specific recombination rate (i.e., heterochiasmy), which is also poorly understood. The family Salmonidae is a model system for these phenomena, having undergone autotetraploidization and subsequent rediploidization in most of the genome at the base of the lineage. The salmonid master sex determining gene is known, and many species have nonhomologous sex chromosomes, putatively due to transposition of this gene. In this study, we identify the sex chromosome of Brook Charr Salvelinus fontinalis and compare sex chromosome identities across the lineage (eight species and four genera). Although nonhomology is frequent, homologous sex chromosomes and other consistencies are present in distantly related species, indicating probable convergence on specific sex and neo-sex chromosomes. We also characterize strong heterochiasmy with 2.7-fold more crossovers in maternal than paternal haplotypes with paternal crossovers biased to chromosome ends. When considering only rediploidized chromosomes, the overall heterochiasmy trend remains, although with only 1.9-fold more recombination in the female than the male. Y chromosome crossovers are restricted to a single end of the chromosome, and this chromosome contains a large interspecific inversion, although its status between males and females remains unknown. Finally, we identify quantitative trait loci (QTL) for 21 unique growth, reproductive, and stress-related phenotypes to improve knowledge of the genetic architecture of these traits important to aquaculture and evolution.

Genomic patterns of diversity and divergence of two introduced salmonid species in Patagonia, South America

Invasive species have become widespread in aquatic environments throughout the world, yet there are few studies that have examined genomic variation of multiple introduced species in newly colonized environments. In this study, we contrast genomic variation in two salmonid species (anadromous Chinook Salmon, Oncorhynchus tshawytscha, 11,579 SNPs and resident Brook Charr Salvelinus fontinalis, 13,522 SNPs) with differing invasion success after introduction to new environments in South America relative to populations from their native range in North America. Estimates of genetic diversity were not significantly different between introduced and source populations for either species, indicative of propagule pressure that has been shown to maintain diversity in founding populations relative to their native range. Introduced populations also demonstrated higher connectivity and gene flow than those in their native range. Evidence for candidate loci under divergent selection was observed, but was limited to specific introduced populations and was not widely evident. Patterns of genomic variation were consistent with general dispersal potential of each species and therefore also the notion that life history variation may contribute to both invasion success and subsequent genetic structure of these two salmonids in Patagonia.

Enhanced transcriptomic responses in the Pacific salmon louse Lepeophtheirus salmonis oncorhynchi to the non-native Atlantic Salmon Salmo salar suggests increased parasite fitness

Background

Outcomes of infections with the salmon louse Lepeophtheirus salmonis vary considerably among its natural hosts (Salmo, Oncorhynchus spp.). Host-parasite interactions range from weak to strong host responses accompanied by high to low parasite abundances, respectively. Parasite behavioral studies indicate that the louse prefers the host Atlantic Salmon (Salmo salar), which is characterized by a weak immune response, and that this results in enhanced parasite reproduction and growth rates. Furthermore, parasite-derived immunosuppressive molecules (e.g., proteases) have been detected at higher amounts in response to the mucus of Atlantic Salmon relative to Coho Salmon (Oncorhynchus kisutch). However, the host-specific responses of the salmon louse have not been well characterized in either of the genetically distinct sub-species that occur in the Atlantic and Pacific Oceans.

Results

We assessed and compared the transcriptomic feeding response of the Pacific salmon louse (L. salmonis oncorhynchi,) while parasitizing the highly susceptible Atlantic Salmon and Sockeye Salmon (Oncorhynchus nerka) or the more resistant Coho Salmon (Oncorhynchus kisutch) using a 38 K oligonucleotide microarray. The response of the louse was enhanced both in the number of overexpressed genes and in the magnitude of expression while feeding on the non-native Atlantic Salmon, compared to either Coho or Sockeye Salmon. For example, putative virulence factors (e.g., cathepsin L, trypsin, carboxypeptidase B), metabolic enzymes (e.g., cytochrome B, cytochrome C), protein synthesis enzymes (e.g., ribosomal protein P2, 60S ribosomal protein L7), and reproduction-related genes (e.g., estrogen sulfotransferase) were overexpressed in Atlantic-fed lice, indicating heightened parasite fitness with this host species. In contrast, responses in Coho- or Sockeye-fed lice were more similar to those of parasites deprived of a host. To test for host acclimation by the parasite, we performed a reciprocal host transfer experiment and determined that the exaggerated response to Atlantic Salmon was independent of the initial host species, confirming our conclusion that the Pacific salmon louse exhibits an enhanced response to Atlantic Salmon.

Conclusions

This study characterized global transcriptomic responses of Pacific salmon lice during infection of susceptible and resistant hosts. Similar parasite responses during infection of Coho or Sockeye Salmon, despite differences in natural immunity to infection between these host species, indicate that host susceptibility status alone does not drive the parasite response. We identified an enhanced louse response after feeding on Atlantic Salmon, characterized by up-regulation of virulence factors, energy metabolism and reproductive-associated transcripts. In contrast, the responses of lice infecting Coho or Sockeye Salmon were weaker, with reduced expression of virulence factors. These observations indicate that the response of the louse is independent of host susceptibility and suggest that co-evolutionary host-parasite relationships may influence contemporary host-parasite interactions. This research improves our understanding of the susceptibility of Atlantic Salmon and may assist in the development of novel control measures against the salmon louse.

Electronic supplementary material

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

Salmonid Chromosome Evolution as Revealed by a Novel Method for Comparing RADseq Linkage Maps

Whole genome duplication (WGD) can provide material for evolutionary innovation. Family Salmonidae is ideal for studying the effects of WGD as the ancestral salmonid underwent WGD relatively recently, ∼65 Ma, then rediploidized and diversified. Extensive synteny between homologous chromosome arms occurs in extant salmonids, but each species has both conserved and unique chromosome arm fusions and fissions. Assembly of large, outbred eukaryotic genomes can be difficult, but structural rearrangements within such taxa can be investigated using linkage maps. RAD sequencing provides unprecedented ability to generate high-density linkage maps for nonmodel species, but can result in low numbers of homologous markers between species due to phylogenetic distance or differences in library preparation. Here, we generate a high-density linkage map (3,826 markers) for the Salvelinus genera (Brook Charr S. fontinalis), and then identify corresponding chromosome arms among the other available salmonid high-density linkage maps, including six species of Oncorhynchus, and one species for each of Salmo, Coregonus, and the nonduplicated sister group for the salmonids, Northern Pike Esox lucius for identifying post-duplicated homeologs. To facilitate this process, we developed MapComp to identify identical and proximate (i.e. nearby) markers between linkage maps using a reference genome of a related species as an intermediate, increasing the number of comparable markers between linkage maps by 5-fold. This enabled a characterization of the most likely history of retained chromosomal rearrangements post-WGD, and several conserved chromosomal inversions. Analyses of RADseq-based linkage maps from other taxa will also benefit from MapComp, available at: https://github.com/enormandeau/mapcomp/

Sex-biased gene expression and sequence conservation in Atlantic and Pacific salmon lice (Lepeophtheirus salmonis)

BACKGROUND

Salmon lice, Lepeophtheirus salmonis (Copepoda: Caligidae), are highly important ectoparasites of farmed and wild salmonids, and cause multi-million dollar losses to the salmon aquaculture industry annually. Salmon lice display extensive sexual dimorphism in ontogeny, morphology, physiology, behavior, and more. Therefore, the identification of transcripts with differential expression between males and females (sex-biased transcripts) may help elucidate the relationship between sexual selection and sexually dimorphic characteristics.

RESULTS

Sex-biased transcripts were identified from transcriptome analyses of three L. salmonis populations, including both Atlantic and Pacific subspecies. A total of 35-43 % of all quality-filtered transcripts were sex-biased in L. salmonis, with male-biased transcripts exhibiting higher fold change than female-biased transcripts. For Gene Ontology and functional analyses, a consensus-based approach was used to identify concordantly differentially expressed sex-biased transcripts across the three populations. A total of 127 male-specific transcripts (i.e. those without detectable expression in any female) were identified, and were enriched with reproductive functions (e.g. seminal fluid and male accessory gland proteins). Other sex-biased transcripts involved in morphogenesis, feeding, energy generation, and sensory and immune system development and function were also identified. Interestingly, as observed in model systems, male-biased L. salmonis transcripts were more frequently without annotation compared to female-biased or unbiased transcripts, suggesting higher rates of sequence divergence in male-biased transcripts.

CONCLUSIONS

Transcriptome differences between male and female L. salmonis described here provide key insights into the molecular mechanisms controlling sexual dimorphism in L. salmonis. This analysis offers targets for parasite control and provides a foundation for further analyses exploring critical topics such as the interaction between sex and drug resistance, sex-specific factors in host-parasite relationships, and reproductive roles within L. salmonis.

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.

Transcriptomic responses to emamectin benzoate in Pacific and Atlantic Canada salmon lice Lepeophtheirus salmonis with differing levels of drug resistance

Salmon lice Lepeophtheirus salmonis are an ecologically and economically important parasite of wild and farmed salmon. In Scotland, Norway, and Eastern Canada, L. salmonis have developed resistance to emamectin benzoate (EMB), one of the few parasiticides available for salmon lice. Drug resistance mechanisms can be complex, potentially differing among populations and involving multiple genes with additive effects (i.e., polygenic resistance). Indicators of resistance development may enable early detection and countermeasures to avoid the spread of resistance. Here, we collect sensitive Pacific L. salmonis and sensitive and resistant Atlantic L. salmonis from salmon farms, propagate in laboratory (F1), expose to EMB in bioassays, and evaluate either baseline (Atlantic only) or induced transcriptomic differences between populations. In all populations, induced responses were minor and a cellular stress response was not identified. Pacific lice did not upregulate any genes in response to EMB, but downregulated degradative enzymes and transport proteins at 50 ppb EMB. Baseline differences between sensitive and now resistant Atlantic lice were much greater than responses to exposures. All resistant lice overexpressed degradative enzymes, and resistant males, the most resistant group, overexpressed collagenases to the greatest extent. These results indicate an accumulation of baseline expression differences related to resistance.

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.

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.

Transcriptomics of coping strategies in free-swimming Lepeophtheirus salmonis (Copepoda) larvae responding to abiotic stress

The salmon louse Lepeophtheirus salmonis is a marine ectoparasite of wild and farmed salmon in the Northern Hemisphere. Infections of farmed salmon are of economic and ecological concern. Nauplius and copepodid salmon lice larvae are free-swimming and disperse in the water column until they encounter a host. In this study, we characterized the sublethal stress responses of L. salmonis copepodid larvae by applying a 38K oligonucleotide microarray to profile transcriptomes following 24 h exposures to suboptimal salinity (30–10 parts per thousand (‰)) or temperature (16–4 °C) environments. Hyposalinity exposure resulted in large-scale gene expression changes relative to those elicited by a thermal gradient. Subsequently, transcriptome responses to a more finely resolved salinity gradient between 30 ‰ and 25 ‰ were profiled. Minimal changes occurred at 29 ‰ or 28 ‰, a threshold of response was identified at 27 ‰, and the largest response was at 25 ‰. Differentially expressed genes were clustered by pattern of expression, and clusters were characterized by functional enrichment analysis. Results indicate larval copepods adopt two distinct coping strategies in response to short-term hyposaline stress: a primary response using molecular chaperones and catabolic processes at 27 ‰; and a secondary response up-regulating ion pumps, transporters, a different suite of chaperones and apoptosis-related transcripts at 26 ‰ and 25 ‰. The results further our understanding of the tolerances of L. salmonis copepodids to salinity and temperature gradients and may assist in the development of salmon louse management strategies.

Ecological transcriptomics of lake-type and riverine sockeye salmon (Oncorhynchus nerka)

BACKGROUND

There are a growing number of genomes sequenced with tentative functions assigned to a large proportion of the individual genes. Model organisms in laboratory settings form the basis for the assignment of gene function, and the ecological context of gene function is lacking. This work addresses this shortcoming by investigating expressed genes of sockeye salmon (Oncorhynchus nerka) muscle tissue. We compared morphology and gene expression in natural juvenile sockeye populations related to river and lake habitats. Based on previously documented divergent morphology, feeding strategy, and predation in association with these distinct environments, we expect that burst swimming is favored in riverine population and continuous swimming is favored in lake-type population. In turn we predict that morphology and expressed genes promote burst swimming in riverine sockeye and continuous swimming in lake-type sockeye.

RESULTS

We found the riverine sockeye population had deep, robust bodies and lake-type had shallow, streamlined bodies. Gene expression patterns were measured using a 16 k microarray, discovering 141 genes with significant differential expression. Overall, the identity and function of these genes was consistent with our hypothesis. In addition, Gene Ontology (GO) enrichment analyses with a larger set of differentially expressed genes found the "biosynthesis" category enriched for the riverine population and the "metabolism" category enriched for the lake-type population.

CONCLUSIONS

This study provides a framework for understanding sockeye life history from a transcriptomic perspective and a starting point for more extensive, targeted studies determining the ecological context of genes.