Search

Clam gardens are habitat modifications established by coastal Indigenous Peoples of northwest North America to enhance intertidal clam habitat productivity to provide secure and reliable local food resources. These gardens were established long ago and are mostly unmaintained in present times. To determine whether unmaintained clam gardens still provide a more productive and beneficial habitat than unmodified clam beaches, Pacific littleneck clams (Leukoma staminea) were transplanted to unmaintained clam gardens or unmodified reference beaches and then evaluated for growth, survival, and transcriptomic signatures of the gill and digestive gland after 16 weeks. All beaches in the study were characterized for sediment grain-size distribution, carbonate concentration, and organic content to identify sediment characteristics that may differ between clam garden and reference sites, as well as potentially contribute to differences in clam productivity. Clam growth and survival, sediment grain-size distribution, and carbonate and organic content were not significantly affected by unmaintained clam gardens. Wide variations in survival, growth, and sediment characteristics were observed among beaches. Examining across beach location resulted in the identification of significant negative effects of small rocks (2.00 – 4.75 mm), and silt (< 63 µm) on both growth and survival. Sand (250 – 500 µm) had a significant positive influence on both growth and survival while fine sand (125 – 250 µm) had a significant positive effect only on growth. Coarse sand (0.50 – 1.00 mm) had a significant positive effect on survival, while very fine sand (63 – 125 µm), carbonate, and organic content all had significant negative effects. To evaluate transcriptomic effects, ade novotranscriptome forL. stamineawas assembled. The final assembly contained 52,000 putative transcripts and those specific to each of the two tissues were identified. This revealed that similar functional categories were enriched in tissue-specific genes, but each tissue had its own transcripts comprising the categories. Transcriptomic analysis revealed differential expression in individuals from clam gardens, and although this effect was small in terms of the numbers of genes, specific response genes were identified consistently in both tissues. In summary, while the unmaintained clam gardens did not impact clam growth and survival over the 16 weeks of the study, it did have an effect at the level of the transcriptome. Furthermore, correlations of transcripts associated with either high or low survival provide new insights into ecological associations of these genes in this non-model organism. In summary, localized environmental factors are likely to have a greater influence on Pacific littleneck clam physiology, growth, and survival than the presence/absence of unmaintained clam gardens.

Despite widespread biodiversity losses, an understanding of how most taxa will respond to future climate change is lacking. Here we integrate genomics and environmental modelling to assess climate change responses in an ecologically and economically important Arctic species. Environmentally associated genomic diversity and machine learning are used to identify highly vulnerable populations of anadromous (migratory) Arctic charr, and we reconstruct estimates of effective population size spanning the twentieth century to identify past climate-associated declines. We uncover past region-wide declines in effective population size that correspond to decreases in temperature and community biomass in the Northwest Atlantic. We find vulnerable populations near the southern range limit, indicating northward shifts and a possible loss of commercially important life-history variation in response to climate change. The genomic approach used here to investigate climate change response identifies past and future declines that impact species persistence, ecosystem stability and food security in the Arctic. Genomics and environmental modelling are integrated to assess past and future changes in Arctic charr populations in response to changing climate. Southern population vulnerability suggests climate change may lead to northward shifts and the loss of important life-history variation.

Coho salmon (Oncorhynchus kisutch) are a culturally and economically important species that return from multiyear ocean migrations to spawn in rivers that flow to the Northern Pacific Ocean. Southern stocks of coho salmon have significantly declined over the past quarter century, and unfortunately, conservation efforts have not reversed this trend. To assist in stock management and conservation efforts, we generated two chromosome-level genome assemblies and sequenced 24 RNA-seq libraries to better annotate the coho salmon genome assemblies. We also resequenced the genomes of 83 coho salmon across their North American range to identify nucleotide variants, characterize the broad effects of isolation-by-distance using a genome-wide association analysis approach, and understand the demographic histories of these salmon by modeling population size from genome-wide data. We observed that more than 13% of all SNPs were associated with latitude (before multiple test correction), likely an affect of isolation-by-distance. From demographic history modeling, we estimated that the SNP latitudinal gradient likely developed as recently as 8,000 years ago. In addition, we identified four genes each harboring multiple SNPs associated with latitude; all of these SNPs were also predicted to modify the function of the gene. Three of these genes have roles in cell junction maintenance and may be involved in osmoregulation. This signifies that ocean salinity may have been a factor influencing coho salmon recolonization after the last glaciation period – generating the current pattern of variation in these three genes.

The genetic underpinnings of incipient speciation, including the genomic mechanisms which contribute to morphological and ecological differentiation and reproductive isolation, remain poorly understood. The repeated evolution of consistently, phenotypically distinct morphs of Arctic Charr (Salvelinus alpinus) within the Quaternary period offer an ideal model to study the repeatability of evolution at the genomic level. Sympatric morphs of Arctic Charr are found across this species' circumpolar distribution. However, the specific genetic mechanisms driving this morph differentiation are largely unknown despite the cultural and economic importance of the anadromous morph. We used a newly designed 87k SNP chip to investigate the character and consistency of the genomic differences among sympatric morphs within three recently deglaciated and geographically proximate lakes in Labrador, Canada. We found genetically distinct small and large morph Arctic Charr in all three lakes consistent with resident and anadromous morphs, respectively. A degree of reproductive isolation among sympatric morphs is likely given genome-wide distributions of outlier SNPs and high genome-wide FST s. Across all lakes, outlier SNPs were largely nonoverlapping suggesting a lack of genetic parallelism driving morph differentiation. Alternatively, several genes and paralogous copies of the same gene consistently differentiated morphs across multiple lakes suggesting their importance to the manifestation of morphs. Our results confirm the utility of Arctic Charr as a model for investigating the predictability of evolution and support the importance of both genetic parallelism and nonparallelism to the incipient speciation of Arctic Charr morphs.

The resiliency of populations and species to environmental change is dependent on the maintenance of genetic diversity, and as such, quantifying diversity is central to combating ongoing widespread reductions in biodiversity. With the advent of next‐generation sequencing, several methods now exist for resolving fine‐scale population structure, but the comparative performance of these methods for genetic assignment has rarely been tested. Here, we evaluate the performance of sequenced microsatellites and a single nucleotide polymorphism (SNP) array to resolve fine‐scale population structure in a critically important salmonid in north eastern Canada, Arctic Charr (Salvelinus alpinus). We also assess the utility of sequenced microsatellites for fisheries applications by quantifying the spatial scales of movement and exploitation through genetic assignment of fishery samples to rivers of origin and comparing these results with a 29‐year tagging dataset. Self‐assignment and simulation‐based analyses of 111 genome‐wide microsatellite loci and 500 informative SNPs from 28 populations of Arctic Charr in north‐eastern Canada identified largely river‐specific genetic structure. Despite large differences (~4X) in the number of loci surveyed between panels, mean self‐assignment accuracy was similar with the microsatellite loci and the SNP panel (>90%). Subsequent analysis of 996 fishery‐collected samples using the microsatellite panel revealed that larger rivers contribute greater numbers of individuals to the fishery and that coastal fisheries largely exploit individuals originating from nearby rivers, corroborating results from traditional tagging experiments. Our results demonstrate the efficacy of sequence‐based microsatellite genotyping to advance understanding of fine‐scale population structure and harvest composition in northern and understudied species.

Copepods encompass numerous ecological roles including parasites, detrivores and phytoplankton grazers. Nonetheless, copepod genome assemblies remain scarce. Lepeophtheirus salmonis is an economically and ecologically important ectoparasitic copepod found on salmonid fish. We present the 695.4 Mbp L. salmonis genome assembly containing ≈60% repetitive regions and 13,081 annotated protein-coding genes. The genome comprises 14 autosomes and a ZZ-ZW sex chromosome system. Assembly assessment identified 92.4% of the expected arthropod genes. Transcriptomics supported annotation and indicated a marked shift in gene expression after host attachment, including apparent downregulation of genes related to circadian rhythm coinciding with abandoning diurnal migration. The genome shows evolutionary signatures including loss of genes needed for peroxisome biogenesis, presence of numerous FNII domains, and an incomplete heme homeostasis pathway suggesting heme proteins to be obtained from the host. Despite repeated development of resistance against chemical treatments L. salmonis exhibits low numbers of many genes involved in detoxification. publishedVersion

The post-glacial colonization of Gander Lake in Newfoundland, Canada, by Arctic Charr (Salvelinus alpinus) provides the opportunity to study the genomic basis of adaptation to extreme deep-water environments. Colonization of deep-water (>50 m) habitats often requires extensive adaptation to cope with novel environmental challenges from high hydrostatic pressure, low temperature, and low light, but the genomic mechanisms underlying evolution in these environments are rarely known. Here, we compare genomic divergence between a deep-water morph adapted to depths of up to 288 m and a larger, piscivorous pelagic morph occupying shallower depths. Using both a SNP array and resequencing of whole nuclear and mitochondrial genomes, we find clear genetic divergence (FST = 0.11-0.15) between deep and shallow water morphs, despite an absence of morph divergence across the mitochondrial genome. Outlier analyses identified many diverged genomic regions containing genes enriched for processes such as gene expression and DNA repair, cardiac function, and membrane transport. Detection of putative copy number variants (CNVs) uncovered 385 genes with CNVs distinct to piscivorous morphs, and 275 genes with CNVs distinct to deep-water morphs, enriched for processes associated with synapse assembly. Demographic analyses identified evidence for recent and local morph divergence, and ongoing reductions in diversity consistent with postglacial colonization. Together, these results show that Arctic Charr morph divergence has occurred through genome-wide differentiation and elevated divergence of genes underlying multiple cellular and physiological processes, providing insight into the genomic basis of adaptation in a deep-water habitat following postglacial recolonization.

Arthropods comprise the most populous groups of animals and show an astonishing diversity. Aquatic arthropods belonging to the sub-class Copepoda encompass a range of ecological roles from endo- and ectoparasites to grazers of phytoplankton that link primary producers to higher trophic levels. Despite the important role of copepods in central ecosystem services and their impact as parasites, representative genomic data and genome assemblies are scarce. This limits our opportunities to understand both the specific biology of individual species and unifying copepod genomic features that may govern their capacities to, for example, adapt to a changing environment. Among the copepod parasites we find Lepeophtheirus salmonis, an important ectoparasite that represents a threat to wild salmonid stocks and causes large annual losses to the salmon farming industry. Here we present the salmon louse genome - the first genome of a parasitic copepod that is fully sequenced and annotated. The 695.4 Mbp assembly was validated by a genetic linkage map and comprises 13 autosomes that recombine almost exclusively in males, one autosome that is shielded from recombination in both sexes and a ZZ-ZW style sex chromosome system. The genome assembly contains approximate 60% repetitive regions and comprise 13081 annotated predicted protein-coding genes. The predicted gene set appears to be quite complete as 92.4% of the expected Arthropod genes were found by a BUSCO. The gene annotations were validated by transcriptome sequencing that corresponds to the expected function of selected tissues. Transcriptome sequencing further revealed a marked shift in the gene expression pattern at the transition from the planktonic dispersal phase to the parasitic lifestyle after host attachment. Among other features, genes related to circadian rhythm are downregulated upon attaching to a host - probably reflecting abandoning a planktonic life with diurnal migration. The genome shows several evolutionary signatures including a large expansion of FNII domains, commonly considered vertebrate specific, and an incomplete heme homeostasis pathway suggesting that heme proteins are obtained from the host. The salmon louse has repeatedly demonstrated a large capacity to develop resistance against chemical treatments. Nonetheless, it exhibits low reduced numbers of several genes commonly involved in detoxification; cytochrome P450, ATP-Binding Cassette type transporters and Glutathione S-transferases. Interestingly, only one gene family with a putative detoxification role was expanded: the major vault proteins. Finally, the salmon louse has lost the ability to sustain peroxisomes, a loss apparently shared in the Caligid family but not among copepods in general.

We present a chromosome-level, long-read genome assembly as a reference for northern pike (Esox lucius) where 97.5% of the genome is chromosome-anchored and N50 falls at 37.5 Mb. Whole-genome resequencing was genotyped using this assembly for 47 northern pike representing six North American populations from Alaska to New Jersey. We discovered that a disproportionate frequency of genetic polymorphism exists among populations east and west of the North American Continental Divide (NACD), indicating reproductive isolation across this barrier. Genome-wide analysis of heterozygous SNP density revealed a remarkable lack of genetic variation with 1 polymorphic site every 6.3kb in the Yukon River drainage and one every 16.5kb east of the NACD. Observed heterozygosity (Ho), nucleotide diversity (π), and Tajima’s D are depressed in populations east of the NACD (east vs. west:Ho: 0.092 vs 0.31; π: 0.092 vs 0.28; Tajima’s D: -1.61 vs -0.47). We confirm the presence of the master sex determining (MSD) gene,amhby, in the Yukon River drainage and in an invasive population in British Columbia and confirm its absence in populations east of the NACD. We also describe an Alaskan population whereamhbyis present but not associated with male gender determination. Our results support that northern pike originally colonized North America through Beringia, that Alaska provided an unglaciated refugium for northern pike during the last ice age, and southeast of the NACD was colonized by a small founding population(s) that lostamhby.

The estimation of linkage disequilibrium between molecular markers within a population is critical when establishing the minimum number of markers required for association studies, genomic selection, and inferring historical events influencing different populations. This work aimed to evaluate the extent and decay of linkage disequilibrium in a coho salmon breeding population using a high-density SNP array. Linkage disequilibrium was estimated between a total of 93,502 SNPs found in 64 individuals (33 dams and 31 sires) from the breeding population. The markers encompass all 30 coho salmon chromosomes and comprise 1,684.62 Mb of the genome. The average density of markers per chromosome ranged from 48.31 to 66 per 1 Mb. The minor allele frequency averaged 0.26 (with a range from 0.22 to 0.27). The overall average linkage disequilibrium among SNPs pairs measured as r 2 was 0.10. The Average r 2 value decreased with increasing physical distance, with values ranging from 0.21 to 0.07 at a distance lower than 1 kb and up to 10 Mb, respectively. An r 2 threshold of 0.2 was reached at distance of approximately 40 Kb. Chromosomes Okis05, Okis15 and Okis28 showed high levels of linkage disequilibrium (>0.20 at distances lower than 1 Mb). Average r 2 values were lower than 0.15 for all chromosomes at distances greater than 4 Mb. An effective population size of 43 was estimated for the population 10 generations ago, and 325, for 139 generations ago. Based on the effective number of chromosome segments, we suggest that at least 74,000 SNPs would be necessary for an association mapping study and genomic predictions. Therefore, the SNP panel used allowed us to capture high-resolution information in the farmed coho salmon population. Furthermore, based on the contemporary N e, a new mate allocation strategy is suggested to increase the effective population size.

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.

Lepeophtheirus salmonis is a naturally occurring marine parasite of salmonid fishes in the Northern hemisphere, and a major problem in salmonid aquaculture. In addition to the direct effects on host fish, L. salmonis may act as a vector for diseases. Here, the microbial community of L. salmonis recovered from whole genome shotgun sequencing was compared between lice sampled from both the Atlantic and the Pacific, laboratory-reared and wild lice, in addition to lice displaying resistance towards chemical treatments. The analysis shows clear differences in the metagenomic composition between the Atlantic and the Pacific Ocean, whereas the resistance status of the L. salmonis or the cultivation did not have significant impact.

Understanding the extent to which ecological divergence is repeatable is essential for predicting responses of biodiversity to environmental change. Here we test the predictability of evolution, from genotype to phenotype, by studying parallel evolution in a salmonid fish, Arctic charr (Salvelinus alpinus), across eleven replicate sympatric ecotype pairs (benthivorous-planktivorous and planktivorous-piscivorous) and two evolutionary lineages. We found considerable variability in eco-morphological divergence, with several traits related to foraging (eye diameter, pectoral fin length) being highly parallel even across lineages. This suggests repeated and predictable adaptation to environment. Consistent with ancestral genetic variation, hundreds of loci were associated with ecotype divergence within lineages of which eight were shared across lineages. This shared genetic variation was maintained despite variation in evolutionary histories, ranging from postglacial divergence in sympatry (ca. 10-15kya) to pre-glacial divergence (ca. 20-40kya) with postglacial secondary contact. Transcriptome-wide gene expression (44,102 genes) was highly parallel across replicates, involved biological processes characteristic of ecotype morphology and physiology, and revealed parallelism at the level of regulatory networks. This expression divergence was not only plastic but in part genetically controlled by parallel cis-eQTL. Lastly, we found that the magnitude of phenotypic divergence was largely correlated with the genetic differentiation and gene expression divergence. In contrast, the direction of phenotypic change was mostly determined by the interplay of adaptive genetic variation, gene expression, and ecosystem size. Ecosystem size further explained variation in putatively adaptive, ecotype-associated genomic patterns within and across lineages, highlighting the role of environmental variation and stochasticity in parallel evolution. Together, our findings demonstrate the parallel evolution of eco-morphology and gene expression within and across evolutionary lineages, which is controlled by the interplay of environmental stochasticity and evolutionary contingencies, largely overcoming variable evolutionary histories and genomic backgrounds., Author summary A renowned natural model system for adaptive evolution is the repeated and rapid divergence of fishes into different sympatric trophic and morphological specialists, known as ecotypes. The drivers and constraints of these repeated divergences are complex and not well understood but it is often observed that postglacial fishes diverge in predictable patterns. Here we use a framework of parallel (or convergent) evolution to test the predictability of divergence in the most variable northern freshwater fish, the Arctic charr. Using a hierarchy of replication—from individuals to divergent phylogeographic lineages—we detect parallel evolution of foraging-related traits despite variation in genomic backgrounds and evolutionary histories. The level of phenotypic parallelism can be explained by the complex interplay of environment, shared genetic variation, and variability in gene expression. While phenotypic divergence may be determined and/or constraint by genetic and molecular divergence, the direction of change is largely determined by molecular parallelism and shared adaptive genetic variation. These are in turn associated with environmental similarity, reflected as ecosystem size. We suggest that gene expression facilitates parallel ecotype evolution, but that the extent of parallelism is further influenced by the level of shared genetic variation and ecological opportunity.

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.

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.

The estimation of linkage disequilibrium between molecular markers within a population is critical when establishing the minimum number of markers required for association studies, genomic selection, and inferring historical events influencing different populations. This work aimed to evaluate the extent and decay of linkage disequilibrium in a coho salmon breeding population using a high-density SNP array. Linkage disequilibrium was estimated between a total of 93,502 SNPs found in 64 individuals (33 dams and 31 sires) from the breeding population. The markers encompass all 30 coho salmon chromosomes and comprise 1,684.62 Mb of the genome. The average density of markers per chromosome ranged from 48.31 to 66 per 1 Mb. The minor allele frequency averaged 0.26 (with a range from 0.22 to 0.27). The overall average linkage disequilibrium among SNPs pairs measured as r2 was 0.10. The Average r2 value decreased with increasing physical distance, with values ranging from 0.21 to 0.07 at a distance lower than 1 kb and up to 10 Mb, respectively. An r2 threshold of 0.2 was reached at distance of approximately 40 Kb. Chromosomes Okis05, Okis15 and Okis28 showed high levels of linkage disequilibrium (>0.20 at distances lower than 1 Mb). Average r2 values were lower than 0.15 for all chromosomes at distances greater than 4 Mb. An effective population size of 43 was estimated for the population 10 generations ago, and 325, for 139 generations ago. Based on the effective number of chromosome segments, we suggest that at least 74,000 SNPs would be necessary for an association mapping study and genomic predictions. Therefore, the SNP panel used allowed us to capture high-resolution information in the farmed coho salmon population. Furthermore, based on the contemporary Ne, a new mate allocation strategy is suggested to increase the effective population size.

Understanding the extent to which evolution is predictable under multifarious selection is a longstanding question in evolutionary biology. However, the interplay of stochastic and contingent factors influencing the extent of parallelism in nature is not well understood. To test the predictability of evolution, we studied a ‘natural experiment’ on different organismal levels across lakes and evolutionary lineages of a freshwater salmonid fish, Arctic charr (Salvelinus alpinus). We identified significant phenotypic parallelism between Arctic charr ecotype pairs within a continuum of parallel evolution and highly parallel adaptive morphological traits. Variability in phenotypic predictability was explained by complex demographic histories, differing genomic backgrounds and genomic responses to selection, variable genetic associations with ecotype, and environmental variation. Remarkably, gene expression was highly similar across ecotype replicates, and explained the observed parallelism continuum. Our findings suggest that parallel evolution by non-parallel evolutionary routes is possible when the regulatory molecular phenotype compensates for divergent histories.

Highlights- First study to highlight parallel epigenetic modifications induced by hatchery rearing as a potential explanatory mechanism for rapid change in fitnessSummaryA puzzling question in conservation biology is how to maintain overall fitness of individuals bred in captive environment upon release into the wild, especially for rehabilitating declining or threatened species [1,2]. For salmonid species, a heritable change in fitness related traits and gene expression has been reported to occur in a single generation of captivity in hatchery environment [3–5]. Such rapid changes are congruent with models of inadvertent domestication selection which may lead to maladaptation in the natural environment [4]. Arguably, the underlying mechanism by which captivity may induce fitness difference between wild and captive congeners is still poorly understood. Short-term selection on complex phenotypic traits is expected to induce subtle changes in allele frequency over multiple loci [7–9]. Yet, most studies investigating the molecular basis for rapid change in fitness related traits occurring in hatchery have concentrated their effort on finding evidence for selection at the genome level by identifying loci with large effect.Numerous wild stocks of Pacific anadromous salmon and trout (genusOncorhynchusandSalmo) have experienced fluctuating abundance over the past century, with a series of sharp declines [6–8]. With the objectives of preserving ecosystem integrity, enhancing declining populations and sustaining fisheries, conservation hatcheries have been flourishing. This is particularly true along the North American Pacific coast where billions of salmonids, all species included, are released each year. Despite substantial improvement of production management, the beneficial ecological role of hatcheries in enhancing and restoring wild stocks is still debated, mainly because of the reduced fitness and maladaptation of hatchery-fish when released in the wild [3,5,9]. Although previous studies showed that domestication selection was involved in such fitness impairment, they also observed that different environmental conditions (e.g. reduced fish density) significantly modulated the physiological acclimation to hatchery environment [4].Environmental stimuli are especially relevant during early embryonic development, which also correspond to a sensitive methylation reprogramming window in vertebrates [10,11]. It is therefore plausible that differences in rearing environment during early development may result in epigenetic modifications that could in turn impact on fitness. However, the only epigenetic study to date pertaining to captive rearing in salmonids and performed using methylation-sensitive amplified fragments (MSAP) failed to identify significant changes in methylation profile associated with hatchery rearing [12]Here, we used a higher resolution approach to compare the genome-wide pattern of methylation in hatchery-reared juvenile (smolt) Coho Salmon with that of their wild counterparts in two geographically distant rivers in British Columbia, Canada. Using a reduced representation bisulfite sequencing (RRBS) approach covering an average per individual of about 70 million cytosines in CpG context, we identified 100 methylated regions (DMRs) that differed in parallel between hatchery and natural origin salmon in both rivers. The total variance of epigenetic variation among individuals explained by river or origin and rearing environment in a RDA model was 16% (adj.R2=0.16), and both variables equally explained about 8% of the variance after controlling for each other. The gene ontology analysis revealed that regions with different methylation levels between hatchery and natural origin salmon showed enrichment for ion homeostasis, synaptic and neuromuscular regulation, immune and stress response, and control of locomotion functions. We further identified 15,044 SNPs that allowed detection of significant differences between either rivers or sexes. However, no effect of rearing environment was observed, confirming that hatchery and natural origin fish of a given river belong to the same panmictic population, as expected based on the hatchery programs applied in these rivers (see Supplementary experimental procedures). Moreover, neither a standard genome-scan approach nor a polygenic statistical framework allowed detection of selective effects within a single generation between hatchery and natural origin salmon. Therefore, this is the first study to demonstrate that parallel epigenetic modifications induced by hatchery rearing during early development may represent a potential explanatory mechanism for rapid change in fitness-related traits previously reported in salmonids.

Much progress has been made regarding our understanding of aromatase regulation, estrogen synthesis partitioning and communication between the germinal and somatic compartments of the differentiating gonad. We now know that most of the enzymatic and signaling apparatus required for steroidogenesis is endogenously expressed within germ cells. However, less is known about the expression and localization of steroidogenic components within mature spermatozoa. We have assembled a sperm library presenting 197,015 putative transcripts. Co-expression clustering analysis revealed that 6687 genes were present at higher levels in sperm in comparison to fifteen other salmon tissue libraries. The sperm transcriptome is highly complex containing the highest proportion of unannotated genes (45%) of the tissues analyzed. Our analysis of highly expressed genes in late-stage sperm revealed dedication to tasks involving chromatin remodeling, flagellogenesis and proteolysis. In addition, using various different embedding and microscopic techniques, we examined the morphology of salmon spermatozoa and characterized expression and localization of several estrogenic regulatory and signaling proteins by immunohistochemistry. We provide evidence for the endogenous synthesis and localization of aromatase (CYP19A and CYP19B1) and potential mediators of estrogen [i.e., ER-alpha and soluble adenylyl cyclase (sAC)] or phosphate (i.e., CREB and FOXL2A) signaling. Partitioning of select transcripts that encode AR-beta, FSH and the LH receptor, but not AR-alpha, LH or the FSH receptor, further points to localized specificity of function in the steroidogenic circuitry of the sperm cell. These results open new avenues of investigation to further our understanding of the intra- and intercellular regulatory processes that guide sperm development and biology.

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

The whole-genome duplication 80 million years ago of the common ancestor of salmonids (salmonid-specific fourth vertebrate whole-genome duplication, Ss4R) provides unique opportunities to learn about the evolutionary fate of a duplicated vertebrate genome in 70 extant lineages. Here we present a high-quality genome assembly for Atlantic salmon (Salmo salar), and show that large genomic reorganizations, coinciding with bursts of transposon-mediated repeat expansions, were crucial for the post-Ss4R rediploidization process. Comparisons of duplicate gene expression patterns across a wide range of tissues with orthologous genes from a pre-Ss4R outgroup unexpectedly demonstrate far more instances of neofunctionalization than subfunctionalization. Surprisingly, we find that genes that were retained as duplicates after the teleost-specific whole-genome duplication 320 million years ago were not more likely to be retained after the Ss4R, and that the duplicate retention was not influenced to a great extent by the nature of the predicted protein interactions of the gene products. Finally, we demonstrate that the Atlantic salmon assembly can serve as a reference sequence for the study of other salmonids for a range of purposes. Supplementary information The online version of this article (doi:10.1038/nature17164) contains supplementary material, which is available to authorized users., The genome sequence is presented for the Atlantic salmon (Salmo salar), providing information about a rediploidization following a salmonid-specific whole-genome duplication event that resulted in an autotetraploidization. Supplementary information The online version of this article (doi:10.1038/nature17164) contains supplementary material, which is available to authorized users., A window on salmonid genome evolution William Davidson and colleagues report sequencing and assembly of the Atlantic salmon genome, which they demonstrate as a useful reference to also improve the genome assembly of other salmanoids. Their analyses provide insights into duplicate retention patterns across two rounds of whole-genome duplication that have occurred in this lineage. Supplementary information The online version of this article (doi:10.1038/nature17164) contains supplementary material, which is available to authorized users.

We report sets of 19 and 16 microsatellite loci for the examination of the population genetics of Boreogadus saida and Arctogadus glacialis, respectively. Six of these loci were developed from a collection of 9,497 expressed sequences from B. saida while the remaining loci were found in the literature and optimized for use in B. saida and A. glacialis. The numbers of alleles observed for each locus ranged from 3 to 33 in B. saida and 1–22 in A. glacialis. Observed heterozygosities ranged from 0.02 to 0.93 in B. saida and 0.17–1.0 in A. glacialis. Species specific differences were observed for the loci providing new tools for the identification of these two morphologically similar arctic gadids. The loci presented here can be used to distinguish between the two species and fill fundamental biological knowledge gaps, thus promoting conservation of these important fishes.

Sea lice are common parasites of both farmed and wild salmon. Salmon farming constitutes an important economic market in North America, South America, and Northern Europe. Infections with sea lice can result in significant production losses. A compilation of genomic information on different genera of sea lice is an important resource for understanding their biology as well as for the study of population genetics and control strategies. We report on over 150,000 expressed sequence tags (ESTs) from five different species (Pacific Lepeophtheirus salmonis (49,672 new ESTs in addition to 14,994 previously reported ESTs), Atlantic L. salmonis (57,349 ESTs), Caligus clemensi (14,821 ESTs), Caligus rogercresseyi (32,135 ESTs), and Lernaeocera branchialis (16,441 ESTs)). For each species, ESTs were assembled into complete or partial genes and annotated by comparisons to known proteins in public databases. In addition, whole mitochondrial (mt) genome sequences of C. clemensi (13,440 bp) and C. rogercresseyi (13,468 bp) were determined and compared to L. salmonis. Both nuclear and mtDNA genes show very high levels of sequence divergence between these ectoparastic copepods suggesting that the different species of sea lice have been in existence for 37–113 million years and that parasitic association with salmonids is also quite ancient. Our ESTs and mtDNA data provide a novel resource for the study of sea louse biology, population genetics, and control strategies. This genomic information provides the material basis for the development of a 38K sea louse microarray that can be used in conjunction with our existing 44K salmon microarray to study host–parasite interactions at the molecular level. This report represents the largest genomic resource for any copepod species to date. Electronic supplementary material The online version of this article (doi:10.1007/s10126-011-9398-z) contains supplementary material, which is available to authorized users.

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.

Salmonids are the most widely studied group of fish, and in the last few years, genomics technologies have begun to contribute to this rich biology. The first salmonid microarrays appeared in 2004 and since then several dozen studies have demonstrated the utility of genomic approaches. The widespread use of the genomics research on Atlantic salmon project 16 k array and greatly expanded genome resources have led to the development of an experimental 5 k oligo (70-mer) array and a 32 k cDNA microarray in the near future. In this paper, the authors examined some of the procedures used in the development of past arrays and reexamined them in light of new genomic data available. Some preliminary control experiments of the new 5 k array were investigated that examine oligo designs based on distance from the polyA tail, the effects of mismatches and cross-species hybridization specificity. Beneficial approaches are also identified in the development of the new 32 k cDNA array.

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.

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. Electronic supplementary material The online version of this article (doi:10.1186/s12862-015-0309-1) contains supplementary material, which is available to authorized users.

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.

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.

Inherited symbionts are ubiquitous in insects and can have important consequences for the fitness of their hosts. Many inherited symbionts defend their hosts against parasites or other natural enemies; however, the means by which most symbionts confer protection is virtually unknown. We examine the mechanisms of defence in a recently discovered case of symbiont-mediated protection, where the bacterial symbiont Spiroplasma defends the fly Drosophila neotestacea from a virulent nematode parasite, Howardula aoronymphium. Using quantitative PCR of Spiroplasma infection intensities and whole transcriptome sequencing, we attempt to distinguish between the following modes of defence: symbiont-parasite competition, host immune priming and the production of toxic factors by Spiroplasma. Our findings do not support a model of exploitative competition between Howardula and Spiroplasma to mediate defence, nor do we find strong support for host immune priming during Spiroplasma infection. Interestingly, we recovered sequence for putative toxins encoded by Spiroplasma, including a novel putative ribosome-inactivating protein, transcripts of which are up-regulated in response to nematode exposure. Protection via the production of toxins may be a widely used and important mechanism in heritable defensive symbioses in insects.

Background The sablefish (order: Scorpaeniformes) is an economically important species in commercial fisheries of the North Pacific and an emerging species in aquaculture. Aside from a handful of sequences in NCBI and a few published microsatellite markers, little is known about the genetics of this species. The development of genetic tools, including polymorphic markers and a linkage map will allow for the successful development of future broodstock and mapping of phenotypes of interest. The significant sexual dimorphism between females and males makes a genetic test for early identification of sex desirable. Results A full mitochondrial genome is presented and the resulting phylogenetic analysis verifies the placement of the sablefish within the Scorpaeniformes. Nearly 35,000 assembled transcript sequences are used to identify genes and obtain polymorphic SNP and microsatellite markers. 360 transcribed polymorphic loci from two sablefish families produce a map of 24 linkage groups. The sex phenotype maps to sablefish LG14 of the male map. We show significant conserved synteny and conservation of gene-order between the threespine stickleback Gasterosteus aculeatus and sablefish. An additional 1843 polymorphic SNP markers are identified through next-generation sequencing techniques. Sex-specific markers and sequence insertions are identified immediately upstream of the gene gonadal-soma derived factor (gsdf), the master sex determinant locus in the medaka species Oryzias luzonensis. Conclusions The first genomic resources for sablefish provide a foundation for further studies. Over 35,000 transcripts are presented, and the genetic map represents, as far as we can determine, the first linkage map for a member of the Scorpaeniformes. The observed level of conserved synteny and comparative mapping will allow the use of the stickleback genome in future genetic studies on sablefish and other related fish, particularly as a guide to whole-genome assembly. The identification of sex-specific insertions immediately upstream of a known master sex determinant implicates gsdf as an excellent candidate for the master sex determinant for sablefish.

Background Classical major histocompatibility complex (MHC) class II molecules play an essential role in presenting peptide antigens to CD4+ T lymphocytes in the acquired immune system. The non-classical class II DM molecule, HLA-DM in the case of humans, possesses critical function in assisting the classical MHC class II molecules for proper peptide loading and is highly conserved in tetrapod species. Although the absence of DM-like genes in teleost fish has been speculated based on the results of homology searches, it has not been definitively clear whether the DM system is truly specific for tetrapods or not. To obtain a clear answer, we comprehensively searched class II genes in representative teleost fish genomes and analyzed those genes regarding the critical functional features required for the DM system. Results We discovered a novel ancient class II group (DE) in teleost fish and classified teleost fish class II genes into three major groups (DA, DB and DE). Based on several criteria, we investigated the classical/non-classical nature of various class II genes and showed that only one of three groups (DA) exhibits classical-type characteristics. Analyses of predicted class II molecules revealed that the critical tryptophan residue required for a classical class II molecule in the DM system could be found only in some non-classical but not in classical-type class II molecules of teleost fish. Conclusions Teleost fish, a major group of vertebrates, do not possess the DM system for the classical class II peptide-loading and this sophisticated system has specially evolved in the tetrapod lineage.

Transcripts for dax1, foxl2, mis and sf1 are co-expressed in the somatic companion cells of teleost germ cells. These regulatory factors function, in part, to modulate the transcription of aromatase, particularly cyp19a, the terminal enzyme of estrogen biosynthesis. At least two separate aromatase loci exist in teleost fish that encode distinct isoforms. The activity of two forms, cyp19a and cyp19b1, is predominantly associated with the ovary and the brain, respectively. We isolated sequences that compose the proximal promoters of cyp19a, cyp19b1 and foxl2a, to identify potential transcription factor binding motifs to define sex-specific regulatory profiles for each gene. We also provide evidence for the translation and immunological localization of DAX-1, FOXL2 and MIS to the endoplasmic reticulum and accumulation within secretory vesicles of the salmon oocyte. We found no evidence for the expression of CYP19A or CYP19B1 in the oocyte at the one-year-old stage. However, synthesis of both aromatases was localized to testicular germ and soma cells at this early stage of development. Production of these regulatory factors in the germ cells may serve to modulate the transcription and activity of endogenous aromatase and/or contribute to the differentiation of the neighbouring companion cells through secretory signaling.

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.

Background Salmonids are one of the most intensely studied fish, in part due to their economic and environmental importance, and in part due to a recent whole genome duplication in the common ancestor of salmonids. This duplication greatly impacts species diversification, functional specialization, and adaptation. Extensive new genomic resources have recently become available for Atlantic salmon (Salmo salar), but documentation of allelic versus duplicate reference genes remains a major uncertainty in the complete characterization of its genome and its evolution. Results From existing expressed sequence tag (EST) resources and three new full-length cDNA libraries, 9,057 reference quality full-length gene insert clones were identified for Atlantic salmon. A further 1,365 reference full-length clones were annotated from 29,221 northern pike (Esox lucius) ESTs. Pairwise dN/dS comparisons within each of 408 sets of duplicated salmon genes using northern pike as a diploid out-group show asymmetric relaxation of selection on salmon duplicates. Conclusions 9,057 full-length reference genes were characterized in S. salar and can be used to identify alleles and gene family members. Comparisons of duplicated genes show that while purifying selection is the predominant force acting on both duplicates, consistent with retention of functionality in both copies, some relaxation of pressure on gene duplicates can be identified. In addition, there is evidence that evolution has acted asymmetrically on paralogs, allowing one of the pair to diverge at a faster rate.

Genomic resources in rainbow smelt (Osmerus mordax) enable us to examine the genome duplication process in salmonids and test hypotheses relating to the fate of duplicated genes. They further enable us to pursue physiological and ecological studies in smelt. A bacterial artificial chromosome library containing 52,410 clones with an average insert size of 146 kb was constructed. This library represents an 11-fold average coverage of the rainbow smelt (O. mordax) genome. In addition, several complementary deoxyribonucleic acid libraries were constructed, and 36,758 sequences were obtained and combined into 12,159 transcripts. Over half of these transcripts have been identified, several of which have been associated with cold adaptation. These basic resources show high levels of similarity (86%) to salmonid genes and provide initial support for genome duplication in the salmonid ancestor. They also facilitate identification of genes important to fish and direct us toward new technologies for other studies in fish biology.

Nuclear deoxyribonucleic acid sequences from approximately 15,000 salmon louse expressed sequence tags (ESTs), the complete mitochondrial genome (16,148bp) of salmon louse, and 16S ribosomal ribonucleic acid (rRNA) and cytochrome oxidase subunit I (COI) genes from 68 salmon lice collected from Japan, Alaska, and western Canada support a Pacific lineage of Lepeophtheirus salmonis that is distinct from that occurring in the Atlantic Ocean. On average, nuclear genes are 3.2% different, the complete mitochondrial genome is 7.1% different, and 16S rRNA and COI genes are 4.2% and 6.1% different, respectively. Reduced genetic diversity within the Pacific form of L. salmonis is consistent with an introduction into the Pacific from the Atlantic Ocean. The level of divergence is consistent with the hypothesis that the Pacific form of L. salmonis coevolved with Pacific salmon (Onchorhynchus spp.) and the Atlantic form coevolved with Atlantic salmonids (Salmo spp.) independently for the last 2.5-11 million years. The level of genetic divergence coincides with the opportunity for migration of fish between the Atlantic and Pacific Ocean basins via the Arctic Ocean with the opening of the Bering Strait, approximately 5 million years ago. The genetic differences may help explain apparent differences in pathogenicity and environmental sensitivity documented for the Atlantic and Pacific forms of L. salmonis.