Search

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 in Canada and the United States 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 a chromosome-level genome assembly. We also resequenced the genomes of 83 coho salmon across the North American range to identify nucleotide variants and understand the demographic histories of these salmon by modeling effective population size from genome-wide data. From demographic history modeling, we observed reductions in effective population sizes between 3,750 and 8,000 years ago for several northern sampling sites, which may correspond to bottleneck events during recolonization after glacial retreat.

Many natural and managed organisms will require substantial functional genetic variation to respond to selection in the face of rapid environmental change. Pacific salmon have experienced strong fluctuations in critical fitness traits over the past five decades. We examined genetic and phenotypic variability over three generations in a pedigreed hatchery population of coho salmon (Oncorhynchus kisutch) by monitoring seven fitness-related traits. Three-year-old adult return numbers varied more than fivefold, and jack (2-year-old males) numbers varied 13-fold. Body sizes of Inch Creek coho salmon decreased consistently such that fish were only 40.7% as heavy in 2015 as in 2006, and female reproductive traits also decreased. During the study period, the majority of families produced returning adult progeny, and effective population size was relatively constant. Heritability estimates for phenotypic traits were significantly greater than zero except for condition factor, and the estimated heritability for jacking was 0.42. The Inch Creek coho salmon population harbours substantial heritability for fitness and reproductive traits and thus likely retains substantial capacity for adaptation despite many years of hatchery propagation.

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.

Conservation units (CUs) are important tools for supporting the implementation of standardized management practices for exploited species. Following the adoption of the Wild Salmon Policy in Canada, CUs were defined for Pacific salmon based on characteristics related to ecotype, life history and genetic variation using microsatellite markers as indirect measures of local adaptation. Genomic data sets have the potential to improve the definition of CUs by reducing variance around estimates of population genetic parameters, thereby increasing the power to detect more subtle patterns of population genetic structure and by providing an opportunity to incorporate adaptive information more directly with the identification of variants putatively under selection. We used one of the largest genomic data sets recently published for a nonmodel species, comprising 5662 individual Coho salmon (

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 potentially significant genetic consequences associated with the loss of migratory capacity of diadromous fishes that have become landlocked in freshwater are poorly understood. Consistent selective pressures associated with freshwater residency may drive repeated differentiation both between allopatric landlocked and anadromous populations and within landlocked populations (resulting in sympatric morphs). Alternatively, the strong genetic drift anticipated in isolated landlocked populations could hinder consistent adaptation, limiting genetic parallelism. Understanding the degree of genetic parallelism underlying differentiation has implications for both the predictability of evolution and management practices. We employed an 87k single-nucleotide polymorphism (SNP) array to examine the genetic characteristics of landlocked and anadromous Arctic char (

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.

Salmon lice have plagued the salmon farming industry and have negatively impacted salmon populations in the wild. In response, researchers have generated high density genetic maps, genome assemblies, transcriptomes, and whole-genome resequencing data to better understand this parasite. In this study, we used long-read sequencing technology to update the previous genome assemblies of Atlantic Ocean salmon lice with a more contiguous assembly and a more comprehensive gene catalog of Pacific Ocean salmon lice. We were also able to further characterize genomic features previously identified from other studies by using published resequenced genomes of 25 Atlantic and 15 Pacific salmon lice. One example was further characterizing the ZW sex chromosomes. For both the Atlantic and Pacific Ocean salmon lice subspecies, we found that the female W-chromosome is only a small fraction of the Z-chromosome and that the vast majority of the W and Z-chromosome do not contain conserved regions (i.e. pseudoautosomal regions). However, conserved orthologous protein sequences can still be identified between the W- and Z-chromosomes.

Chum salmon (Oncorhynchus keta) is the species with the widest geographic range of the anadromous Pacific salmonids,. Chum salmon is the second largest of the Pacific salmon, behind Chinook salmon, and considered the most plentiful Pacific salmon by overall biomass. This species is of significant commercial and economic importance: on average the commercial chum salmon fishery has the second highest processed value of the Pacific salmon within British Columbia. The aim of this work was to establish genomic baseline resources for this species. Our first step to accomplish this goal was to generate a chum salmon reference genome assembly from a doubled-haploid chum salmon. Gene annotation of this genome was facilitated by an extensive RNA-seq database we were able to create from multiple tissues. Range-wide resequencing of chum salmon genomes allowed us to categorize genome-wide geographic variation, which in turn reinforced the idea that genetic differentiation was best described on a regional, rather than at a stock-specific, level. Within British Columbia, chum salmon regional groupings were described at the conservation unit (CU) level, and there may be substructure within particular CUs. Genome wide associations of phenotypic sex to SNP genetic markers identified two clear peaks, a very strong peak on Linkage Group 15, and another on Linkage Group 3. With these new resources, we were better able to characterize the sex-determining region and gain further insights into sex determination in chum salmon and the general biology of this species.

Inferring the genomic basis of local adaptation is a long-standing goal of evolutionary biology. Beyond its fundamental evolutionary implications, such knowledge can guide conservation decisions for populations of conservation and management concern. Here, we investigated the genomic basis of local adaptation in the Coho salmon (Oncorhynchus kisutch) across its entire North American range. We hypothesized that extensive spatial variation in environmental conditions and the species homing behavior may promote the establishment of local adaptation. We genotyped 7,829 individuals representing 217 sampling locations at more than 100,000 high-quality RADseq loci to investigate how recombination might affect the detection of loci putatively under selection and took advantage of the precise description of the demographic history of the species from our previous work to draw accurate population genomic inferences about local adaptation. Results indicated that genetic differentiation scans and genetic-environment association analyses were both significantly affected by variation in recombination rate as low recombination regions displayed an increased number of outliers. By taking these confounding factors into consideration, we revealed that migration distance was the primary selective factor driving local adaptation and partial parallel divergence among distant populations. Moreover, we identified several candidates SNP associated with long distance migration and altitude including a gene known to be involved in adaptation to altitude in other species. The evolutionary implications of our findings are discussed along with conservation applications.

Inferring the genomic basis of local adaptation is a long-standing goal of evolutionary biology. Beyond its fundamental evolutionary implications, such knowledge can guide conservation decisions for populations of conservation and management concern. Here, we investigated the genomic basis of local adaptation in the Coho salmon (Oncorhynchus kisutch) across its entire North American range. We hypothesized that extensive spatial variation in environmental conditions and the species' homing behaviour may promote the establishment of local adaptation. We genotyped 7829 individuals representing 217 sampling locations at more than 100,000 high-quality RADseq loci to investigate how recombination might affect the detection of loci putatively under selection and took advantage of the precise description of the demographic history of the species from our previous work to draw accurate population genomic inferences about local adaptation. The results indicated that genetic differentiation scans and genetic-environment association analyses were both significantly affected by variation in recombination rate as low recombination regions displayed an increased number of outliers. By taking these confounding factors into consideration, we revealed that migration distance was the primary selective factor driving local adaptation and partial parallel divergence among distant populations. Moreover, we identified several candidate single nucleotide polymorphisms associated with long-distance migration and altitude including a gene known to be involved in adaptation to altitude in other species. The evolutionary implications of our findings are discussed along with conservation applications.

Pink salmon (Oncorhynchus gorbuscha) adults are the smallest of the five Pacific salmon native to the western Pacific Ocean. Pink salmon are also the most abundant of these species and account for a large proportion of the commercial value of the salmon fishery worldwide. A strict two-year life-history of most pink salmon generates temporally isolated populations that spawn either in even-years or odd-years. To uncover the influence of this genetic isolation, reference genome assemblies were generated for each year-class and whole genome re-sequencing data was collected from salmon of both year-classes. The salmon were sampled from six Canadian rivers and one Japanese river. At multiple centromeres we identified peaks of Fst between year-classes that were millions of base-pairs long. The largest Fst peak was also associated with a million base-pair chromosomal polymorphism found in the odd-year genome near a centromere. These Fst peaks may be the result of centromere drive or a combination or reduced recombination and genetic drift, and they could influence speciation. Other regions of the genome influenced by odd-year and even-year temporal isolation and tentatively under selection were mostly associated with genes related to immune function, organ development/maintenance, and behaviour.

Males and females often differ in their fitness optima for shared traits that have a shared genetic basis, leading to sexual conflict. Morphologically differentiated sex chromosomes can resolve this conflict and protect sexually antagonistic variation, but they accumulate deleterious mutations. However, how sexual conflict is resolved in species that lack differentiated sex chromosomes is largely unknown. Here we present a chromosome-anchored genome assembly for rainbow trout (Oncorhynchus mykiss) and characterize a 55-Mb double-inversion supergene that mediates sex-specific migratory tendency through sex-dependent dominance reversal, an alternative mechanism for resolving sexual conflict. The double inversion contains key photosensory, circadian rhythm, adiposity and sex-related genes and displays a latitudinal frequency cline, indicating environmentally dependent selection. Our results show sex-dependent dominance reversal across a large autosomal supergene, a mechanism for sexual conflict resolution capable of protecting sexually antagonistic variation while avoiding the homozygous lethality and deleterious mutations associated with typical heteromorphic sex chromosomes.

The avermectins, emamectin benzoate (EMB) and ivermectin (IVM) have been commonly used in North America over the last two decades to control the salmon louse, Lepeophtheirus salmonis, infections in farmed Atlantic salmon. Emamectin benzoate, trade name SLICE ™, was used heavily in the Eastern Canadian industry between the years 2000–2008, due to its long lasting protection and efficacy against all parasitic life stages. However, over reliance on this drug soon resulted in reduced sensitivity in many L. salmonis populations, resulting more recently in uses of higher treatment dosages and switching to the use of IVM. For these reasons, we investigated the effects of different dosages of EMB and multiple IVM treatments on baseline immunophysiological indicators, anti-viral responses and protection against subsequent salmon lice exposure in salmon smolts. Different doses of EMB or repeated treatment with IVM did not affect feeding behaviour in salmon, however by the end of the second IVM treatment, some neurotoxicity was observed. A single (1×) EMB dose (50 μg/Kg) administered for 7 consecutive days had no significant effect on L. salmonis abundance and development, whereas triple the dosage (150 μg/Kg) significantly reduced lice development, thereby eliminating subsequent stress responses in salmon associated with lice development to pre-adult stages. Emamectin benzoate and IVM treatment did not significantly impact expression of resting antigen presentation molecules in salmon (MH class I or II), however they did inhibit short-term (6 h) induced responses to the ISA virus. The impact of gender, as previously shown, had the greatest effect on louse transcriptomic regulation, but avermectin treatment also caused perturbations in gene expression. Transcriptome differences between lice on control and 1× EMB treated fish were larger than those observed for IVM or 3× EMB. Nearly half of the transcripts differentially expressed by IVM were also affected by one of the EMB treatments. Transcriptomic results from the louse suggest a high degree of similarity and concordance within and across studies in avermectin treatment, with gender of louse and dosage of drug significantly impacting the outcomes.

A current challenge for environmental DNA (eDNA) applications is how to account for an environmental (or false-positive) background in surveys. We performed two controlled experiments in the Goldstream Hatchery in British Columbia using a validated coho salmon (Oncorhynchus kisutch) eDNA assay (eONKI4). In the density experiment at high copy number, eDNA in 2 L water samples was measured from four 10 kL tanks containing 1 to 65 juvenile coho salmon. At these densities, we obtained a strong positive 1:1 relationship between predicted copy number/L and coho salmon biomass (g/L). The dilution experiment simulated a situation where fish leave a pool environment, and water from upstream continues to flow through at rates of 141–159 L/min. Here, three coho salmon were placed in four 10 kL tanks, removed after nine days, and the amount of remaining eDNA was measured at times coinciding with dilutions of 20, 40, 80, 160, and 1000 kL. The dilution experiment demonstrates a novel method using Binomial–Poisson distributions to detect target species eDNA at low copy number in the presence of an environmental background. This includes determination of the limit of blank with background (LOB-B) with a controlled false positive rate, and limit of detection with background (LOD-B) with a controlled false negative rate, which provides a statistically robust “Detect” or “No Detect” assessment for eDNA surveys.

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.

Genotype-by-environment (GxE) interactions are non-parallel reaction norms among individuals with different genotypes in response to different environmental conditions. GxE interactions are an extension of phenotypic plasticity and consequently studying such interactions improves our ability to predict effects of different environments on phenotype as well as the fitness of genetically distinct organisms and their capacity to interact with ecosystems. Growth hormone transgenic coho salmon grow much faster than non-transgenics when raised in tank environments, but show little difference in growth when reared in nature-like streams. We used this model system to evaluate potential mechanisms underlying this growth rate GxE interaction, performing RNA-seq to measure gene transcription and whole-genome bisulfite sequencing to measure gene methylation in liver tissue. Gene ontology (GO) term analysis revealed stress as an important biological process potentially influencing growth rate GxE interactions. While few genes with transcription differences also had methylation differences, in promoter or gene regions, many genes were differentially methylated between tank and stream environments. A GO term analysis of differentially methylated genes between tank and stream environments revealed increased methylation in the stream environment of more than 95% of the differentially methylated genes, many with biological processes unrelated to liver function. The lower nutritional condition of the stream environment may cause increased negative regulation of genes less vital for liver tissue function than when fish are reared in tanks with unlimited food availability. These data show a large effect of rearing environment both on gene expression and methylation, but it is less clear that the detected epigenetic marks are responsible for the observed altered growth and physiological responses.

Landscape genetic analyses of wildlife populations can exclude variation in a broad suite of potential spatiotemporal correlates, including consideration of how such variation might have similarly influenced people over time. Grizzly bear (Ursus arctos) populations in what is now known as coastal British Columbia, Canada, provide an opportunity to examine the possible effects of a complex set of landscape and human influences on genetic structure. In this collaboration among the Nuxalk, HaíÉ«zaqv, Kitasoo/Xai'xais, Gitga'at, and Wuikinuxv First Nations and conservation scientists, we characterized patterns of genetic differentiation in the grizzly bear, a species of high cultural value, by genotyping 22 microsatellite loci from noninvasively collected hair samples over a 23,500 km² area. We identified three well-differentiated groups. Resistance surfaces, which incorporated past and present human use, settlement, and landscape resistant features, could not explain this pattern of genetic variation. Notably, however, we detected spatial alignment between Indigenous language families and grizzly bear genetic groups. Grizzly bears sampled within an area represented by a given language family were significantly similar to those sampled within that language family (P = 0.001) and significantly divergent to those sampled outside the language family (P = 0.001). This spatial co-occurrence suggests that grizzly bear and human groups have been shaped by the landscape in similar ways, creating a convergence of grizzly bear genetic and human linguistic diversity. Additionally, grizzly bear management units designated by the provincial government currently divide an otherwise continuous group and exclude recently colonized island populations that are genetically continuous with adjacent mainland groups. This work provides not only insight into how ecological and geographic conditions can similarly shape the distribution of people and wildlife but also new genetic evidence to support renewed, locally led management of grizzly bears into the future.

Whole genome duplication (WGD) events have played a major role in eukaryotic genome evolution, but the consequence of these extreme events in adaptive genome evolution is still not well understood. To address this knowledge gap we used a comparative phylogenetic model and transcriptomic data from seven species to infer selection on gene expression in duplicated genes (ohnologs) following the salmonid WGD 80-100 million years ago. We find rare cases of tissue-specific expression evolution but pervasive expression evolution affecting many tissues, reflecting strong selection on maintenance of genome stability following genome doubling. Although ohnolog expression levels have evolved mostly asymmetrically, by diverting one ohnolog copy down a path towards pseudogenization, strong evolutionary constraints have frequently also favoured symmetric shifts in gene dosage of both copies, likely to achieve gene dose reduction while avoiding accumulation of ‘toxic mutations’. Mechanistically, ohnolog regulatory divergence is dictated by the number of bound transcription factors in promoters, with transposable elements being one source of novel binding sites driving tissue-specific gains in expression. Our results imply pervasive adaptive expression evolution following WGD to overcome the immediate challenges posed by genome doubling and to exploit the long-term genetic opportunities for novel phenotype evolution.

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.

New computational methods and next‐generation sequencing (NGS) approaches have enabled the use of thousands or hundreds of thousands of genetic markers to address previously intractable questions. The methods and massive marker sets present both new data analysis challenges and opportunities to visualize, understand, and apply population and conservation genomic data in novel ways. The large scale and complexity of NGS data also increases the expertise and effort required to thoroughly and thoughtfully analyze and interpret data. To aid in this endeavor, a recent workshop entitled “Population Genomic Data Analysis,” also known as “ConGen 2017,” was held at the University of Montana. The ConGen workshop brought 15 instructors together with knowledge in a wide range of topics including NGS data filtering, genome assembly, genomic monitoring of effective population size, migration modeling, detecting adaptive genomic variation, genomewide association analysis, inbreeding depression, and landscape genomics. Here, we summarize the major themes of the workshop and the important take‐home points that were offered to students throughout. We emphasize increasing participation by women in population and conservation genomics as a vital step for the advancement of science. Some important themes that emerged during the workshop included the need for data visualization and its importance in finding problematic data, the effects of data filtering choices on downstream population genomic analyses, the increasing availability of whole‐genome sequencing, and the new challenges it presents. Our goal here is to help motivate and educate a worldwide audience to improve population genomic data analysis and interpretation, and thereby advance the contribution of genomics to molecular ecology, evolutionary biology, and especially to the conservation of biodiversity.

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.

In teleost fish as in mammals, humoral adaptive immunity is based on B lymphocytes expressing highly diverse immunoglobulins (IG). During B cell differentiation, IG loci are subjected to genomic rearrangements of V, D, and J genes, producing a unique antigen receptor expressed on the surface of each lymphocyte. During the course of an immune response to infections or immunizations, B cell clones specific of epitopes from the immunogen are expanded and activated, leading to production of specific antibodies. Among teleost fish, salmonids comprise key species for aquaculture. Rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar) are especially important from a commercial point of view and have emerged as critical models for fish immunology. The growing interest to capture accurate and comprehensive antibody responses against common pathogens and vaccines has resulted in recent efforts to sequence the IG repertoire in these species. In this context, a unified and standardized nomenclature of salmonid IG heavy chain (IGH) genes is urgently required, to improve accuracy of annotation of adaptive immune receptor repertoire dataset generated by high-throughput sequencing (AIRRseq) and facilitate comparisons between studies and species. Interestingly, the assembly of salmonids IGH genomic sequences is challenging due to the presence of two large size duplicated IGH loci and high numbers of IG genes and pseudogenes. We used data available for Atlantic salmon to establish an IMGT standardized nomenclature of IGH genes in this species and then applied the IMGT rules to the rainbow trout IGH loci to set up a nomenclature, which takes into account the specificities of Salmonid loci. This unique, consistent nomenclature for Salmonid IGH genes was then used to construct IMGT sequence reference directories allowing accurate annotation of AIRRseq data. The complex issues raised by the genetic diversity of salmon and trout strains are discussed in the context of IG repertoire annotation.

A thorough reconstruction of historical processes is essential for a comprehensive understanding the mechanisms shaping patterns of genetic diversity. Indeed, past and current conditions influencing effective population size have important evolutionary implications for the efficacy of selection, increased accumulation of deleterious mutations, and loss of adaptive potential. Here, we gather extensive genome-wide data that represent the extant diversity of the Coho salmon (Oncorhynchus kisutch) to address two objectives. We demonstrate that a single glacial refugium is the source of most of the present-day genetic diversity, with detectable inputs from a putative secondary micro-refugium. We found statistical support for a scenario whereby ancestral populations located south of the ice sheets expanded in postglacial time, swamping out most of the diversity from other putative micro-refugia. Demographic inferences revealed that genetic diversity was also affected by linked selection in large parts of the genome. Moreover, we demonstrate that the recent demographic history of this species generated regional differences in the load of deleterious mutations among populations, a finding that mirrors recent results from human populations and provides increased support for models of expansion load. We propose that insights from these historical inferences should be better integrated in conservation planning of wild organisms, which currently focuses largely on neutral genetic diversity and local adaptation, with the role of potentially maladaptive variation being generally ignored.

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.

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.

Traits with different fitness optima in males and females cause sexual conflict when they have a shared genetic basis. Heteromorphic sex chromosomes can resolve this conflict and protect sexually antagonistic polymorphisms but accumulate deleterious mutations. However, many taxa lack differentiated sex chromosomes, and how sexual conflict is resolved in these species is largely unknown. Here we present a chromosome-anchored genome assembly for rainbow trout (Oncorhynchus mykiss) and characterize a 56 Mb double-inversion supergene that mediates sex-specific migration through sex-dependent dominance, a mechanism that reduces sexual conflict. The double-inversion contains key photosensory, circadian rhythm, adiposity, and sexual differentiation genes and displays frequency clines associated with latitude and temperature, revealing environmental dependence. Our results constitute the first example of sex-dependent dominance across a large autosomal supergene, a novel mechanism for sexual conflict resolution capable of protecting polygenic sexually antagonistic variation while avoiding the homozygous lethality and deleterious mutation load of heteromorphic sex chromosomes.

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

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.

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.

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.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Carotenoids are primarily responsible for the characteristic red flesh coloration of salmon. Flesh coloration is an economically and evolutionarily significant trait that varies inter- and intra-specifically, yet the underlying genetic mechanism is unknown. Chinook salmon ( Oncorhynchus tshawytscha ) represents an ideal system to study carotenoid variation as, unlike other salmonids, they exhibit extreme differences in carotenoid utilization due to genetic polymorphisms. Here, we crossed populations of Chinook salmon with fixed differences in flesh coloration (red versus white) for a genome-wide association study to identify loci associated with pigmentation. Here, the beta-carotene oxygenase 2-like ( BCO2-l ) gene was significantly associated with flesh colour, with the most significant single nucleotide polymorphism explaining 66% of the variation in colour. BCO2 gene disruption is linked to carotenoid accumulation in other taxa, therefore we hypothesize that an ancestral mutation partially disrupting BCO2-l activity (i.e. hypomorphic mutation) allowed the deposition and accumulation of carotenoids within Salmonidae. Indeed, we found elevated transcript levels of BCO2-l in white Chinook salmon relative to red. The long-standing mystery of why salmon are red, while no other fishes are, is thus probably explained by a hypomorphic mutation in the proto-salmonid at the time of divergence of red-fleshed salmonid genera (approx. 30 Ma).

Molecular genetic tools have become standard in biological studies of both model and non-model species. This has created a growing need for sequence information, a resource hitherto limited for many species. With new sequencing technologies this is rapidly changing, and whole genome shotgun sequencing has become a realistic goal for many species. However, present sequencing protocols require more DNA than can be extracted from single individuals of many small metazoans, potentially forcing sequencing projects to perform sequencing on samples derived from several individuals. A pertinent question thus arises: can wild samples be used or is inbreeding necessary? In the present study we compare assemblies generated using sequence data from inbred and wild Lepeophtheirus salmonis. The results indicate not only that measures to reduce the genetic variability may significantly improve the final assemblies but also that deeper coverage to some extent can compensate for the detrimental effects of natural sequence...

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

We describe an emerging initiative - the ‘Functional Analysis 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, through 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.

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.

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. Electronic supplementary material The online version of this article (doi:10.1186/s13104-016-2030-6) 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.

Whole genome duplications (WGDs) are considered to have been a driving force in the generation of evolutionary diversity that is characteristic of higher eukaryotes. The ancestor of salmonids underwent two additional WGDs compared to mammals, one (3R) at the base of the teleost radiation and another (4R) in the common ancestor of extant salmonids. We have chosen the fatty acid binding protein (fabp) gene family as a model to study the fate of duplicated genes in teleosts following WGDs. As previously described for zebrafish, we identified two copies (fabp7a and fabp7b) of the brain-type fabp gene in several fish including rainbow smelt, but there was only a single transcript in northern pike, the closest relative of the salmonids, and two rather than the expected four fabp7 genes in Atlantic salmon, rainbow trout and grayling. A phylogenetic analysis revealed that a loss of the fabp7a gene occurred in the common ancestor of the northern pike and salmonids after it had diverged from the rainbow smelt, and that the 4R WGD then gave rise to the fabp7bI and fabp7bII observed in salmonids. This is supported by genetic mapping that placed the Atlantic salmon duplicated fabp7b genes on homeologous chromosomes. There was no evidence of neo-functionalization in the salmonid fabp7bI and fabp7bII genes based on dN/dS ratios and an examination of amino acid substitutions. Atlantic salmon fabp7bI and fabp7bII genes are both expressed broadly like fabp7b expression in northern pike. However, only Atlantic salmon fabp7bII, like its counterpart in northern pike and zebrafish, was expressed in the liver. A comparison of ~2000bp upstream of Atlantic salmon fabp7b gene duplicates revealed an insertion of 62bp in fabp7bI relative to fabp7bII. The presence of predicted transcription factor binding sites in this insertion sequence may explain the differential expression of the fabp7b gene duplicates in Atlantic salmon liver.