Effect of triploidy on liver gene expression in coho salmon (Oncorhynchus kisutch) under different metabolic states

BACKGROUND

Triploid coho salmon are excellent models for studying gene dosage and the effects of increased cell volume on gene expression. Triploids have an additional haploid genome in each cell and have fewer but larger cells than diploid coho salmon to accommodate the increased genome size. Studying gene expression in triploid coho salmon provides insight into how gene expression may have been affected after the salmonid-specific genome duplication which occurred some 90 MYA. Triploid coho salmon are sterile and consequently can live longer and grow larger than diploid congeners in many semelparous species (spawning only once) because they never reach maturity and post-spawning mortality is averted. Triploid fishes are also of interest to the commercial sector (larger fish are more valuable) and to fisheries management since sterile fish can potentially minimize negative impacts of escaped fish in the wild.

RESULTS

The vast majority of genes in liver tissue had similar expression levels between diploid and triploid coho salmon, indicating that the same amount of mRNA transcripts were being produced per gene copy (positive gene dosage effects) within a larger volume cell. Several genes related to nutrition and compensatory growth were differentially expressed between diploid and triploid salmon, indicating that some loci are sensitive to cell size and/or DNA content per cell. To examine how robust expression between ploidies is under different conditions, a genetic/metabolic modifier in the form of different doses of a growth hormone transgene was used to assess gene expression under conditions that the genome has not naturally experienced or adapted to. While many (up to 1400) genes were differentially expressed between non-transgenic and transgenic fish, relatively few genes were differentially expressed between diploids and triploids with similar doses of the transgene. These observations indicate that the small effect of ploidy on gene expression is robust to large changes in physiological state.

CONCLUSIONS

These findings are of interest from a gene regulatory perspective, but also valuable for understanding phenotypic effects in triploids, transgenics, and triploid transgenics that could affect their utility in culture conditions and their fitness and potential consequences of release into nature.

Design and characterization of an 87k SNP genotyping array for Arctic charr (Salvelinus alpinus)

We have generated a high-density, high-throughput genotyping array for characterizing genome-wide variation in Arctic charr (Salvelinus alpinus). Novel single nucleotide polymorphisms (SNPs) were identified in charr from the Fraser, Nauyuk and Tree River aquaculture strains, which originated from northern Canada and fish from Iceland using high coverage sequencing, reduced representation sequencing and RNA-seq datasets. The array was designed to capture genome-wide variation from a diverse suite of Arctic charr populations. Cross validation of SNPs from various sources and comparison with previously published Arctic charr SNP data provided a set of candidate SNPs that generalize across populations. Further candidate SNPs were identified based on minor allele frequency, association with RNA transcripts, even spacing across intergenic regions and association with the sex determining (sdY) gene. The performance of the 86,503 SNP array was assessed by genotyping Fraser, Nauyuk and Tree River strain individuals, as well as wild Icelandic Arctic charr. Overall, 63,060 of the SNPs were polymorphic within at least one group and 36.8% were unique to one of the four groups, suggesting that the array design allows for characterization of both within and across population genetic diversity. The concordance between sdY markers and known phenotypic sex indicated that the array can accurately determine the sex of individuals based on genotype alone. The Salp87k genotyping array provides researchers and breeders the opportunity to analyze genetic variation in Arctic charr at a more detailed level than previously possible.

A genetic linkage map for the salmon louse (Lepeophtheirus salmonis): evidence for high male:female and inter-familial recombination rate differences

A salmon louse (Lepeophtheirus salmonis salmonis) genetic linkage map was constructed to serve as a genomic resource for future investigations into the biology of this important marine parasitic copepod species, and to provide insights into the inheritance patterns of genetic markers in this species. SNP genotyping of 8 families confirmed the presence of 15 linkage groups based upon the assignment of 93,773 markers. Progeny sample size weight adjusted map sizes in males (with the exception of SL12 and SL15) ranged in size from 96.50 cM (SL11) to 134.61 cM (SL06), and total combined map steps or bins ranged from 143 (SL09) to 203 (SL13). The SL12 male map was the smallest linkage group with a weight-averaged size of 3.05 cM with 6 recombination bins. Male:female specific recombination rate differences are 10.49:1 and represent one of the largest reported sex-specific differences for any animal species. Recombination ratio differences (M:F) ranged from 1.0 (SL12) to 29:1 (SL15). The number of markers exhibiting normal Mendelian segregation within the sex linkage group SL15 was extremely low (N = 80) in comparison to other linkage groups genotyped [range: 1459 (SL12)-10206 markers (SL05)]. Re-evaluation of Mendelian inheritance patterns of markers unassigned to any mapping parent according to hemizygous segregation patterns (models presented) identified matches for many of these markers to hemizygous patterns. The greatest proportion of these markers assigned to SL15 (N increased to 574). Inclusion of the hemizygous markers revised SL15 sex-specific recombination rate differences to 28:1. Recombination hot- and coldspots were identified across all linkage groups with all linkage groups possessing multiple peaks. Nine of 13 linkage groups evaluated possessed adjacent domains with hot-coldspot transitional zones. The most common pattern was for one end of the linkage to show elevated recombination in addition to internal regions. For SL01 and SL06, however, a terminal region with high recombination was not evident while a central domain possessing extremely high-recombination levels was present. High levels of recombination were weakly coupled to higher levels of SNP variation within domains, but this association was very strong for the central domains of SL01 and SL06. From the pooled paternal half-sib lots (several virgin females placed with 1 male), only 1 or two surviving family lots were obtained. Surviving families possessed parents where both the male and female possessed either inherently low or high recombination rates. This study provides insight into the organization of the sea louse genome, and describes large differences in recombination rate that exist among individuals of the same sex, and between the sexes. These differences in recombination rate may be coupled to the capabilities of this species to adapt to environmental and pharmaceutical treatments, given that family survivorship appears to be enhanced when parents have similar recombination levels.

Caligus rogercresseyi acetylcholinesterase types and variants: a potential marker for organophosphate resistance

Background

Control of the sea louse Caligus rogercresseyi in the Chilean salmonid industry is reliant on chemical treatments. Azamethiphos was introduced in 2013, although other organophosphates were previously used. In 2014, reduced sensitivity to azamethiphos was detected in the Los Lagos Region using bioassays. The main target of organophosphates is the enzyme acetylcholinesterase (AChE). Mutations in the AChE gene are the main cause of organophosphate resistance in arthropods, including other sea lice. In the present study, we aimed to characterize C. rogercresseyi AChE(s) gene(s) and to study the association between AChE variants and azamethiphos resistance in this sea louse species.

Methods

Samples of adult male and female C. rogercresseyi were collected in the Los Lagos Region in 2014. Twenty-four hour exposure bioassays with azamethiphos were performed to select sensitive and resistant lice. The full-length cDNA coding sequences encoding for two AChEs in C. rogercresseyi were molecularly characterized. One of the AChE genes was screened by direct sequencing in the azamethiphos-selected lice to search for variants. An additional louse sampling was performed before and after an azamethiphos treatment in the field in 2017 to validate the findings.

Results

The molecular analysis revealed two putative AChEs in C. rogercresseyi. In silico analysis and 3D modelling of the protein sequences identified both of them as invertebrate AChE type 1; they were named C. rogercresseyi AChE1a and 1b. AChE1a had the characteristics of the main synaptic AChE, while AChE1b lacked some of the important amino acids of a typical AChE. A missense change found in the main synaptic AChE (1a), F318F/V (F290 in Torpedo californica), was associated with survival of C. rogercresseyi at high azamethiphos concentrations (bioassays and field treatment). The amino acid change was located in the acyl pocket of the active-site gorge of the protein.

Conclusions

The present study demonstrates the presence of two types of AChE1 genes in C. rogercresseyi. Although enzymatic assays are needed, AChE1a is most probably the main synaptic AChE. The function of AChE1b is unknown, but evidence points to a scavenger role. The AChE1a F/V318 variant is most probably involved in organophosphate resistance, and can be a good marker for resistance monitoring.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3151-7) contains supplementary material, which is available to authorized users.

The Arctic charr (Salvelinus alpinus) genome and transcriptome assembly

Arctic charr have a circumpolar distribution, persevere under extreme environmental conditions, and reach ages unknown to most other salmonids. The Salvelinus genus is primarily composed of species with genomes that are structured more like the ancestral salmonid genome than most Oncorhynchus and Salmo species of sister genera. It is thought that this aspect of the genome may be important for local adaptation (due to increased recombination) and anadromy (the migration of fish from saltwater to freshwater). In this study, we describe the generation of a new genetic map, the sequencing and assembly of the Arctic charr genome (GenBank accession: GCF_002910315.2) using the newly created genetic map and a previous genetic map, and present several analyses of the Arctic charr genes and genome assembly. The newly generated genetic map consists of 8,574 unique genetic markers and is similar to previous genetic maps with the exception of three major structural differences. The N50, identified BUSCOs, repetitive DNA content, and total size of the Arctic charr assembled genome are all comparable to other assembled salmonid genomes. An analysis to identify orthologous genes revealed that a large number of orthologs could be identified between salmonids and many appear to have highly conserved gene expression profiles between species. Comparing orthologous gene expression profiles may give us a better insight into which genes are more likely to influence species specific phenotypes.

Regulatory processes that control haploid expression of salmon sperm mRNAs

Objective

Various stages of mRNA processing are necessary for functionally important genes required during late-stage sperm differentiation. Protein–RNA complexes form that edit, stabilize, store, deliver, localize and regulate translation of sperm mRNAs. These regulatory processes are often directed by recognition sequence elements and the particular composition of the proteins associated with the mRNAs. Previous work has shown that the cAMP response element modulator (CREM), estrogen receptor-alpha (ERα) and forkhead box L2A (FOXL2A) proteins are present in late-stage salmon sperm. Here we investigate whether these and other regulatory proteins might control processing of mRNAs not expressed until the haploid stage of development. We also examine regulatory processes that prepare and present mRNAs that generate unique products essential for differentiating sperm (i.e. for flagellar assembly and function).

Results

We provide evidence for potential sperm-specific recognition elements in 5′-untranslated regions (utrs) that may bind CREM, ERα, FOXL2A, Y-box and other proteins. We show that changes within the 5′-utrs and open reading frames of some sperm genes lead to distinct protein termini that may provide specific interfaces necessary for localization and function within the paternal gamete.

Electronic supplementary material

The online version of this article (10.1186/s13104-018-3749-z) contains supplementary material, which is available to authorized users.

High level efficacy of lufenuron against sea lice (Lepeophtheirus salmonis) linked to rapid impact on moulting processes

Drug resistance in the salmon louse Lepeophtheirus salmonis is a global issue for Atlantic salmon aquaculture. Multiple resistance has been described across most available compound classes with the exception of the benzoylureas. To target this gap in effective management of L. salmonis and other species of sea lice (e.g. Caligus spp.), Elanco Animal Health is developing an in-feed treatment containing lufenuron (a benzoylurea) to be administered prior to seawater transfer of salmon smolts and to provide long-term protection of salmon against sea lice infestations. Benzoylureas disrupt chitin synthesis, formation, and deposition during all moulting events. However, the mechanism(s) of action are not yet fully understood and most research completed to date has focused on insects. We exposed the first parasitic stage of L. salmonis to 700 ppb lufenuron for three hours and observed over 90% reduction in survival to the chalimus II life stage on the host, as compared to vehicle controls. This agrees with a follow up in vivo administration study on the host, which showed >95% reduction by the chalimus I stage. Transcriptomic responses of salmon lice exposed to lufenuron included genes related to moulting, epithelial differentiation, solute transport, and general developmental processes. Global metabolite profiles also suggest that membrane stability and fluidity is impacted in treated lice. These molecular signals are likely the underpinnings of an abnormal moulting process and cuticle formation observed ultrastructurally using transmission electron microscopy. Treated nauplii-staged lice exhibited multiple abnormalities in the integument, suggesting that the coordinated assembly of the epi- and procuticle is impaired. In all cases, treatment with lufenuron had rapid impacts on L. salmonis development. We describe multiple experiments to characterize the efficacy of lufenuron on eggs, larvae, and parasitic stages of L. salmonis, and provide the most comprehensive assessment of the physiological responses of a marine arthropod to a benzoylurea chemical.

A 200K SNP chip reveals a novel Pacific salmon louse genotype linked to differential efficacy of emamectin benzoate

Antiparasitic drugs such as emamectin benzoate (EMB) are relied upon to reduce the parasite load, particularly of the sea louse Lepeophtheirus salmonis, on farmed salmon. The decline in EMB treatment efficacy for this purpose is an important issue for salmon producers around the world, and particularly for those in the Atlantic Ocean where widespread EMB tolerance in sea lice is recognized as a significant problem. Salmon farms in the Northeast Pacific Ocean have not historically experienced the same issues with treatment efficacy, possibly due to the relatively large population of endemic salmonid hosts that serve to both redistribute surviving lice and dilute populations potentially under selection by introducing naïve lice to farms. Frequent migration of lice among farmed and wild hosts should limit the effect of farm-specific selection pressures on changes to the overall allele frequencies of sea lice in the Pacific Ocean. A previous study using microsatellites examined L. salmonis oncorhynchi from 10 Pacific locations from wild and farmed hosts and found no population structure. Recently however, a farm population of sea lice was detected where EMB bioassay exposure tolerance was abnormally elevated. In response, we have developed a Pacific louse draft genome that complements the previously-released Atlantic louse sequence. These genomes were combined with whole-genome re-sequencing data to design a highly sensitive 201,279 marker SNP array applicable for both subspecies (90,827 validated Pacific loci; 153,569 validated Atlantic loci). Notably, kmer spectrum analysis of the re-sequenced samples indicated that Pacific lice exhibit a large within-individual heterozygosity rate (average of 1 in every 72 bases) that is markedly higher than that of Atlantic individuals (1 in every 173 bases). The SNP chip was used to produce a high-density map for Atlantic sea louse linkage group 5 that was previously shown to be associated with EMB tolerance in Atlantic lice. Additionally, 478 Pacific louse samples from farmed and wild hosts obtained between 2005 and 2014 were also genotyped on the array. Clustering analysis allowed us to detect the apparent emergence of an otherwise rare genotype at a high frequency among the lice collected from two farms in 2013 that had reported elevated EMB tolerance. This genotype was not observed in louse samples collected from the same farm in 2010, nor in any lice sampled from other locations prior to 2013. However, this genotype was detected at low frequencies in louse samples from farms in two locations reporting elevated EMB tolerance in 2014. These results suggest that a rare genotype present in Pacific lice may be locally expanded in farms after EMB treatment. Supporting this hypothesis, 437 SNPs associated with this genotype were found to be in a region of linkage group 5 that overlaps the region associated with EMB resistance in Atlantic lice. Finally, five of the top diagnostic SNPs within this region were used to screen lice that had been subjected to an EMB survival assay, revealing a significant association between these SNPs and EMB treatment outcome. To our knowledge this work is the first report to identify a genetic link to altered EMB efficacy in L. salmonis in the Pacific Ocean.

Chinook salmon (Oncorhynchus tshawytscha) genome and transcriptome

When unifying genomic resources among studies and comparing data between species, there is often no better resource than a genome sequence. Having a reference genome for the Chinook salmon (Oncorhynchus tshawytscha) will enable the extensive genomic resources available for Pacific salmon, Atlantic salmon, and rainbow trout to be leveraged when asking questions related to the Chinook salmon. The Chinook salmon's wide distribution, long cultural impact, evolutionary history, substantial hatchery production, and recent wild-population decline make it an important research species. In this study, we sequenced and assembled the genome of a Chilliwack River Hatchery female Chinook salmon (gynogenetic and homozygous at all loci). With a reference genome sequence, new questions can be asked about the nature of this species, and its role in a rapidly changing world.

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.

Functional Annotation of All Salmonid Genomes (FAASG): an international initiative supporting future salmonid research, conservation and aquaculture

We describe an emerging initiative - the 'Functional Annotation of All Salmonid Genomes' (FAASG), which will leverage the extensive trait diversity that has evolved since a whole genome duplication event in the salmonid ancestor, to develop an integrative understanding of the functional genomic basis of phenotypic variation. The outcomes of FAASG will have diverse applications, ranging from improved understanding of genome evolution, to improving the efficiency and sustainability of aquaculture production, supporting the future of fundamental and applied research in an iconic fish lineage of major societal importance.

Host–parasite transcriptomics during immunostimulant-enhanced rejection of salmon lice (Lepeophtheirus salmonis) by Atlantic salmon (Salmo salar)

Salmon lice ( Lepeophtheirus salmonis) are important ectoparasites of wild and farmed salmonids and cause major losses to the salmon farming industry throughout the Northern Hemisphere. With the emergence of resistance to several commonly used parasiticides, novel control strategies and integration of multiple treatment options are needed, including host immunostimulation. Here, we investigate the effects of a functional feed containing a peptidoglycan and nucleotide formulation on L. salmonis infection of Atlantic salmon ( Salmo salar) by characterizing lice infection levels, the expression of several host immune genes, and the parasite transcriptomic response to the immunostimulated host. Although initial infection intensities were low, the low dose (LD) immunostimulant diet reduced the total lice burden by 50% relative to controls. Immunostimulant fed hosts up-regulated interleukin-1β in the skin and spleen. This gene has been implicated in successful responses of several salmonid species to salmon lice but is typically not observed in Atlantic salmon, suggesting a favorable influence on the immune response. Lice infecting LD immunostimulated salmon overexpressed genes putatively involved in parasite immunity, including carboxylesterases, and underexpressed genes putatively involved in feeding (e.g., proteases). These lice response genes further improve the characterization of the transcriptome of the non-model parasite by identifying genes potentially involved in evading host immunity.

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.

Cypermethrin exposure induces metabolic and stress-related gene expression in copepodid salmon lice (Lepeophtheirus salmonis)

Cypermethrin has been administered for decades to control salmon lice (Lepeophtheirus salmonis) infestations in Atlantic salmon farming regions globally. However, resistance to cypermethrin and other available therapeutants has threatened the sustainability of this growing industry. To better understand the effects of cypermethrin on L. salmonis, a 38K oligonucleotide microarray and RT-qPCR analyses were applied to pools of copepodid larvae exposed to 1.0ppb cypermethrin or seawater controls for 24h. Phenotypic assessments and global gene expression profiles showed a significant disruption of homeostasis in copepodid L. salmonis exposed to cypermethrin. Multiple degradative enzymes were overexpressed in cypermethrin-treated lice including five trypsin-like serine proteases and three cytochrome p450s CYP3a24 (p=0.03, fold change (FC)=3.8; GenBank accession no. JP326960.1), CYP6w1 (p=0.008, FC=5.3; GenBank accession no. JP317875.1), and CYP6d4 (p=0.01; FC=7.9; GenBank accession no. JP334550.1). These enzymes represent preliminary markers for understanding the physiological response of L. salmonis to cypermethrin exposure. A general stress response was also observed in cypermethrin-treated lice which included differential expression of cell signaling genes involved in the induction of cell growth, solute transport, and metabolism. Lastly, a consensus-based analysis was completed with two previously published L. salmonis transcriptome studies revealing genes that respond to cypermethrin, emamectin benzoate (another delousing agent) and hyposalinity. This included concordant differential expression of heat shock beta-1, ammonium transporter Rh types B, and 72kDa type IV collagenase across different L. salmonis studies. This is currently the most comprehensive transcriptome assessment of chemical exposure on the first infectious stage of L. salmonis, providing novel markers for studying drug resistance and general stress in this important parasite.

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.

A PCR assay detects a male-specific duplicated copy of Anti-Müllerian hormone (amh) in the lingcod (Ophiodon elongatus)

BACKGROUND

Anti-Müllerian hormone (amh) or Müllerian-inhibiting substance (mis) is a member of the transforming growth factor-β family of hormones. This gene plays a key role in vertebrate male sex-determination by inhibiting the development of the Müllerian ducts, and has been shown to be the master sex-determinant in the Patagonian pejerrey.

RESULTS

In the lingcod, Ophiodon elongatus, both males and females share one copy of amh, however we have identified a second duplicate copy that appears solely in the male individuals. We have developed a PCR-based assay targeting the TGF-β domain of amh that provides a simple method with which to sex lingcod from a small amount of tissue. An analysis across 57 individuals gave a 100% success rate in identifying the phenotypic sex.

CONCLUSIONS

We present a simple method to sex lingcod through non-lethal tissue sampling. A third, independent, male-specific duplication of amh in a teleost fish has been identified in the lingcod.

Multi-tissue transcriptome profiles for coho salmon (Oncorhynchus kisutch), a species undergoing rediploidization following whole-genome duplication

Salmonids are an important family of fish both from economic and basic research perspectives, and have been subjected to extensive research at whole-animal and molecular levels. Most research to date has been conducted on Atlantic salmon (Salmo salar) and rainbow trout (Oncorhynchus mykiss), but more recently other salmonids have become a focus of study due to their interesting life histories and because of their potential for use in commercial aquaculture. However, molecular biology and genetic analyses for these emerging species are currently hampered due to the lack of extensive genomic resources. To overcome some of these limitations, we have constructed a 43,228 sequence transcriptome from 13 tissues from coho salmon, Oncorhynchus kisutch using de novo transcriptome assembly methods. The transcriptome profiling analysis has provided data distinguishing allelic variation from paralogues that arose during the recent whole-genome duplication event in this family, thus allowing simplified analysis of gene-specific expression. Additionally, 1599 novel coho sequences have been identified through comparison with transcriptomes from two other salmonids species (Atlantic salmon and rainbow trout), and with northern pike. The transcriptome presented here will be useful for genomic analysis of coho salmon and other closely related salmonid species.

Differential modulation of resistance biomarkers in skin of juvenile and mature pink salmon, Oncorhynchus gorbuscha by the salmon louse, Lepeophtheirus salmonis

Juvenile pink salmon larger than 0.7 g reject the sea louse, Lepeophtheirus salmonis, and are considered resistant to the infection. Robust innate defense responses in the skin contribute to the observed resistance. In contrast adult pink salmon captured at sea or shortly before spawning carry large numbers of the parasite, suggesting inability to control the infection. The purpose of this research is to better understand these apparently contradictory conclusions by comparing a suite of genetic and cellular markers of resistance to L. salmonis in the skin of juvenile and mature pink salmon. The expression of major histocompatibility factor II, C-reactive protein, interleukin-1β, interleukin-8 and cyclooxygenase-2 was down-regulated in mature but not juvenile pink salmon. Similarly, skin at the site of parasite attachment in juvenile salmon was highly populated with MHIIβ(+) and IL-1β(+) cells that were either absent, or at reduced levels at similar sites in mature salmon. In addition, mucocyte density was relatively low in the skin of mature salmon, irrespective of louse infection. In juveniles, the higher mucocyte density decreased following louse attachment. We show that in mature pink salmon, genetic and histological responses in skin are depressed and speculate that salmonid defense against L. salmonis is modulated by maturation.

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.

A comprehensive analysis of teleost MHC class I sequences

BACKGROUND

MHC class I (MHCI) molecules are the key presenters of peptides generated through the intracellular pathway to CD8-positive T-cells. In fish, MHCI genes were first identified in the early 1990's, but we still know little about their functional relevance. The expansion and presumed sub-functionalization of cod MHCI and access to many published fish genome sequences provide us with the incentive to undertake a comprehensive study of deduced teleost fish MHCI molecules.

RESULTS

We expand the known MHCI lineages in teleosts to five with identification of a new lineage defined as P. The two lineages U and Z, which both include presumed peptide binding classical/typical molecules besides more derived molecules, are present in all teleosts analyzed. The U lineage displays two modes of evolution, most pronouncedly observed in classical-type alpha 1 domains; cod and stickleback have expanded on one of at least eight ancient alpha 1 domain lineages as opposed to many other teleosts that preserved a number of these ancient lineages. The Z lineage comes in a typical format present in all analyzed ray-finned fish species as well as lungfish. The typical Z format displays an unprecedented conservation of almost all 37 residues predicted to make up the peptide binding groove. However, also co-existing atypical Z sub-lineage molecules, which lost the presumed peptide binding motif, are found in some fish like carps and cavefish. The remaining three lineages, L, S and P, are not predicted to bind peptides and are lost in some species.

CONCLUSIONS

Much like tetrapods, teleosts have polymorphic classical peptide binding MHCI molecules, a number of classical-similar non-classical MHCI molecules, and some members of more diverged MHCI lineages. Different from tetrapods, however, is that in some teleosts the classical MHCI polymorphism incorporates multiple ancient MHCI domain lineages. Also different from tetrapods is that teleosts have typical Z molecules, in which the residues that presumably form the peptide binding groove have been almost completely conserved for over 400 million years. The reasons for the uniquely teleost evolution modes of peptide binding MHCI molecules remain an enigma.

Chemokine receptors in Atlantic salmon

Teleost sequence data have revealed that many immune genes have evolved differently when compared to other vertebrates. Thus, each gene family needs functional studies to define the biological role of individual members within major species groups. Chemokine receptors, being excellent markers for various leukocyte subpopulations, are one such example where studies are needed to decipher individual gene function. The unique salmonid whole genome duplication that occurred approximately 95 million years ago has provided salmonids with many additional duplicates further adding to the complexity and diversity. Here we have performed a systematic study of these receptors in Atlantic salmon with particular focus on potential inflammatory receptors. Using the preliminary salmon genome data we identified 48 chemokine or chemokine-like receptors including orthologues to the ten receptors previously published in trout. We found expressed support for 40 of the bona fide salmon receptors. Eighteen of the chemokine receptors are duplicated, and when tested against a diploid sister group the majority were shown to be remnants of the 4R whole genome duplication with subsequent high sequence identity. The salmon chemokine receptor repertoire of 40 expressed bona fide genes is comparably larger than that found in humans with 23 receptors. Diversification has been a major driving force for these duplicate genes with the main variability residing in ligand binding and signalling domains.

Signatures of resistance to Lepeophtheirus salmonis include a TH2-type response at the louse-salmon interface

Disease outbreaks with the salmon louse Lepeophtheirus salmonis cause significant economic losses in mariculture operations worldwide. Variable innate immune responses at the louse-attachment site contribute to differences in susceptibility among species such that members of Salmo spp. are more susceptible to infection than those of some Oncorhynchus spp. Relatively little is known about the mechanisms that contribute to disease resistance or susceptibility to L. salmonis in salmon. Here, we utilize histochemistry and transcriptomics in a comparative infection model with susceptible (Atlantic, sockeye) and resistant (coho) salmon. At least three cell populations (MHIIβ+, IL1β+, TNFα+) were activated in coho salmon skin during L. salmonis infection. Locally elevated expression of several pro-inflammatory mediators (e.g. IL1β, IL8, TNFα, COX2, C/EBPβ), and tissue repair enzymes (MMP9, MMP13) were detected in susceptible and resistant species. However, responses specific to coho salmon (e.g. IL4, IL6, TGFβ) or responses shared among susceptible salmon (e.g. SAP, TRF, Cath in Atlantic and sockeye salmon) provide evidence for species-specific pathways contributing to resistance or susceptibility, respectively. Our results confirm the importance of an early pro-inflammatory TH1-type pathway as an initial host response during infection with Pacific sea lice, and demonstrate subsequent regulatory TH2-type processes as candidate defense mechanisms in the skin of resistant coho salmon.

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.

Atlantic salmon possesses two clusters of type I interferon receptor genes on different chromosomes, which allows for a larger repertoire of interferon receptors than in zebrafish and mammals

Mammalian type I interferons (IFNs) signal through a receptor composed of the IFNAR1 and IFNAR2 chains. In zebrafish two-cysteine IFNs utilize a receptor composed of CRFB1 and CRFB5, while four-cysteine IFNs signal through a receptor formed by CRFB2 and CRFB5. In the present work two CRFB clusters were identified in different chromosomes of Atlantic salmon. Genes of three CRFB5s, one CRFB1, one CRFB2 and the novel CRFB5x were identified, cloned and studied functionally. All CRFBs were expressed in 10 different organs, but the relative expression of CRFBs varied. Mx-reporter assay was used to study which CRFBs might be involved in receptors for salmon IFNa, IFNb and IFNc. The results of Mx-reporter assays suggest that IFNa signals through a receptor composed of CRFB1a as the long chain and either CRFB5a, CRFB5b or CRFB5c as the short chain; IFNc signals through a receptor with CRFB5a or CRFB5c as the short chain while IFNb may signal through a receptor with CRFB5x as a short chain. Taken together, the present work demonstrates that Atlantic salmon has a more diverse repertoire of type I IFN receptors compared to zebrafish or mammals.

The genome and linkage map of the northern pike (Esox lucius): conserved synteny revealed between the salmonid sister group and the Neoteleostei

The northern pike is the most frequently studied member of the Esociformes, the closest order to the diverse and economically important Salmoniformes. The ancestor of all salmonids purportedly experienced a whole-genome duplication (WGD) event, making salmonid species ideal for studying the early impacts of genome duplication while complicating their use in wider analyses of teleost evolution. Studies suggest that the Esociformes diverged from the salmonid lineage prior to the WGD, supporting the use of northern pike as a pre-duplication outgroup. Here we present the first genome assembly, reference transcriptome and linkage map for northern pike, and evaluate the suitability of this species to provide a representative pre-duplication genome for future studies of salmonid and teleost evolution. The northern pike genome sequence is composed of 94,267 contigs (N50 = 16,909 bp) contained in 5,688 scaffolds (N50 = 700,535 bp); the total scaffolded genome size is 878 million bases. Multiple lines of evidence suggest that over 96% of the protein-coding genome is present in the genome assembly. The reference transcriptome was constructed from 13 tissues and contains 38,696 transcripts, which are accompanied by normalized expression data in all tissues. Gene-prediction analysis produced a total of 19,601 northern pike-specific gene models. The first-generation linkage map identifies 25 linkage groups, in agreement with northern pike's diploid karyotype of 2N = 50, and facilitates the placement of 46% of assembled bases onto linkage groups. Analyses reveal a high degree of conserved synteny between northern pike and other model teleost genomes. While conservation of gene order is limited to smaller syntenic blocks, the wider conservation of genome organization implies the northern pike exhibits a suitable approximation of a non-duplicated Protacanthopterygiian genome. This dataset will facilitate future studies of esocid biology and empower ongoing examinations of the Atlantic salmon and rainbow trout genomes by facilitating their comparison with other major teleost groups.