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.

BackgroundWhole 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.ResultsWe 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. Ohnolog expression levels have evolved mostly asymmetrically, by diverting one ohnolog copy down a path towards lower expression and possible pseudogenization. Loss of expression in one ohnolog is significantly associated with transposable element insertions in promoters and likely driven by selection on gene dosage including selection on stoichiometric balance. We also find symmetric expression shifts, and these are associated with genes under strong evolutionary constraints such as ribosome subunit genes. This possibly reflects selection operating to achieve a 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 likely source of novel binding sites driving tissue-specific gains in expression.ConclusionsOur 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.

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.

Growth hormone (GH) transgenic fish often exhibit remarkable transformations in growth rate and other phenotypes relative to wild-type. The 5750A transgenic coho salmon strain exhibits strong sexually dimorphic growth, with females possessing growth stimulation at a level typical of that seen for both sexes in other strains harbouring the same gene construct (e.g. M77), while males display a modest level of growth stimulation. GH mRNA levels were significantly higher in females than in males of the 5750A strain but equivalent in the M77 strain, indicating sex and transgene insertion locus altered transgene expression. We found that acute estradiol treatments did not influence GH expression in either strain (5750A and M77) or the transgene promoter (metallothionein-B), suggesting that estradiol level was not a significant factor influencing transgene activity. The feminization of XX and XY fish of the 5750A and M77 strains generated all-female groups and resulted in equalized growth of the two genetic sexes, suggesting that the presence of the Y chromosome was not directly capable of influencing the GH transgene–mediated growth in a physiological female conditions. These data suggest that the difference in growth rate seen between the sexes in the 5750A strain arises from non-estradiol-mediated sex influences on gene regulation at the transgene locus. This study shows how genetic factors and transgene insertion sites can influence transgene expression with significant consequent effects on phenotype. Supplementary Information The online version contains supplementary material available at 10.1007/s10126-020-10012-5.

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.

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.

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

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.

Background Triploid coho salmon are excellent models for studying gene dosage and the effects of increased cell volume on gene expression. Triploids have an additional haploid genome in each cell and have fewer but larger cells than diploid coho salmon to accommodate the increased genome size. Studying gene expression in triploid coho salmon provides insight into how gene expression may have been affected after the salmonid-specific genome duplication which occurred some 90 MYA. Triploid coho salmon are sterile and consequently can live longer and grow larger than diploid congeners in many semelparous species (spawning only once) because they never reach maturity and post-spawning mortality is averted. Triploid fishes are also of interest to the commercial sector (larger fish are more valuable) and to fisheries management since sterile fish can potentially minimize negative impacts of escaped fish in the wild. Results The vast majority of genes in liver tissue had similar expression levels between diploid and triploid coho salmon, indicating that the same amount of mRNA transcripts were being produced per gene copy (positive gene dosage effects) within a larger volume cell. Several genes related to nutrition and compensatory growth were differentially expressed between diploid and triploid salmon, indicating that some loci are sensitive to cell size and/or DNA content per cell. To examine how robust expression between ploidies is under different conditions, a genetic/metabolic modifier in the form of different doses of a growth hormone transgene was used to assess gene expression under conditions that the genome has not naturally experienced or adapted to. While many (up to 1400) genes were differentially expressed between non-transgenic and transgenic fish, relatively few genes were differentially expressed between diploids and triploids with similar doses of the transgene. These observations indicate that the small effect of ploidy on gene expression is robust to large changes in physiological state. Conclusions These findings are of interest from a gene regulatory perspective, but also valuable for understanding phenotypic effects in triploids, transgenics, and triploid transgenics that could affect their utility in culture conditions and their fitness and potential consequences of release into nature. Electronic supplementary material The online version of this article (10.1186/s12864-019-5655-8) contains supplementary material, which is available to authorized users.

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.

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.

A thorough reconstruction of historical processes is essential for a comprehensive understanding of 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 recently, 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.

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.

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 two-year life history of 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 a centromere drive or a combination of 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.

Animal domestication and artificial selection give rise to gradual changes at the genomic level in populations. Subsequent footprints of selection known as selection signatures or selective sweeps have been traced in the genomes of many animal livestock species by exploiting variations in linkage disequilibrium patterns and/or reduction of genetic diversity.Domestication of most aquatic species is recent in comparison with land animals, and salmonids are one of the most important fish species in aquaculture. Coho salmon (Oncorhynchus kisutch), cultivated primarily in Chile, has been subject to breeding programs to improve growth, disease resistance traits, and flesh color. This study aimed to identify selection signatures that may be involved in adaptation to culture conditions and traits of productive interest. To do so, individuals of two domestic populations cultured in Chile were genotyped with 200 thousand SNPs, and analyses were conducted using iHS, XP-EHH and CLR. Several signatures of selection on different chromosomal regions were detected across both populations. Some of the identified regions under selection contained genes such anapc2, alad, chp2 and myn that have been previously associated with body weight in Atlantic salmon or sec24d and robo1 that have been associated with disease resistance to Piscirickettsia salmonis in coho salmon. Findings in our study can contribute to an integrated genome-wide map of selection signatures, to help identify the genetic mechanisms of phenotypic diversity in coho salmon.

A thorough reconstruction of historical processes is essential for a comprehensive understanding of 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 recently, 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., Author summary Reconstruction of a species’ past demographic history from genetic data can highlight historical factors that have shaped the distribution of genetic diversity along its genome and its geographic range. Here, we combine genotyping-by-sequencing with demographic modelling to address these issues in the Coho salmon, a Pacific salmon of conservation concern in some parts of its range, notably in the south. Our demographic reconstructions reveal a linear decrease in genetic diversity toward the north of the species range, supporting the hypothesis of a northern route of postglacial recolonization from a single major southern refugium. As predicted by theory, we also observed a higher proportion of deleterious mutations in the most distant populations from this refugium. Beyond this general pattern, among-site variation in the proportion of deleterious mutations is consistent with different local trends in effective population sizes. Our results highlight the potential importance of understanding historical factors that have shaped geographic patterns of the distribution of deleterious mutations in order to implement effective management programs for the conservation of wild populations. Such fundamental knowledge of human historical demography is now having major impacts on health sciences, and we argue it is time to integrate such approaches in conservation science as well.

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.

The characterization of runs of homozygosity (ROH), using high-density single nucleotide polymorphisms (SNPs) allows inferences to be made about the past demographic history of animal populations and the genomic ROH has become a common approach to characterize the inbreeding. We aimed to analyze and characterize ROH patterns and compare different genomic and pedigree-based methods to estimate the inbreeding coefficient in two pure lines (POP A and B) and one recently admixed line (POP C) of coho salmon (Oncorhynchus kisutch) breeding nuclei, genotyped using a 200 K Affymetrix Axiom&reg, myDesign Custom SNP Array. A large number and greater mean length of ROH were found for the two &ldquo, pure&rdquo, lines and the recently admixed line (POP C) showed the lowest number and smaller mean length of ROH. The ROH analysis for different length classes suggests that all three coho salmon lines the genome is largely composed of a high number of short segments (&lt, 4 Mb), and for POP C no segment &gt, 16 Mb was found. A high variable number of ROH, mean length and inbreeding values across chromosomes, positively the consequence of artificial selection. Pedigree-based inbreeding values tended to underestimate genomic-based inbreeding levels, which in turn varied depending on the method used for estimation. The high positive correlations between different genomic-based inbreeding coefficients suggest that they are consistent and may be more accurate than pedigree-based methods, given that they capture information from past and more recent demographic events, even when there are no pedigree records available.

Sockeye salmon (Oncorhynchus nerka) is a commercially and culturally important species to the people that live along the northern Pacific Ocean coast. There are two main sockeye salmon ecotypes-the ocean-going (anadromous) ecotype and the fresh-water ecotype known as kokanee. The goal of this study was to better understand the population structure of sockeye salmon and identify possible genomic differences among populations and between the two ecotypes. In pursuit of this goal, we generated the first reference sockeye salmon genome assembly and an RNA-seq transcriptome data set to better annotate features of the assembly. Resequenced whole-genomes of 140 sockeye salmon and kokanee were analyzed to understand population structure and identify genomic differences between ecotypes. Three distinct geographic and genetic groups were identified from analyses of the resequencing data. Nucleotide variants in an immunoglobulin heavy chain variable gene cluster on chromosome 26 were found to differentiate the northwestern group from the southern and upper Columbia River groups. Several candidate genes were found to be associated with the kokanee ecotype. Many of these genes were related to ammonia tolerance or vision. Finally, the sex chromosomes of this species were better characterized, and an alternative sex-determination mechanism was identified in a subset of upper Columbia River kokanee.

Sexual reproduction is a ubiquitous basic feature of life and genetic sex determination is thus widespread, at least among eukaryotes. Understanding the remarkable diversity of sex determination mechanisms, however, is limited by the paucity of empirical studies. Here, we traced back the evolution of sex determination in an entire clade of vertebrates and uncovered that the northern pike (Esox lucius) master sex-determining gene initiated from a 65 to 90 million-year-old gene duplication and remained sex-linked on undifferentiated sex chromosomes for at least 56 million years. Contrasting with its ancient origin, we identified several independent species- or population-specific transitions of sex determination mechanisms in this lineage, including an unexpected complete and recent Y-chromosome loss in some North American northern pike populations. These findings highlight the diversity of the evolutionary fates of master sex-determining genes and raise the importance of careful considerations of population demographic history in sex determination studies. Our study also puts forward the hypothesis that occasional sex reversals and genetic bottlenecks provide a non-adaptive explanation for sex determination transitions.

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.