The sockeye salmon genome, transcriptome, and analyses identifying population defining regions of the genome

Conserved islands of divergence associated with adaptive variation in sockeye salmon are maintained by multiple mechanisms

Local adaptation is facilitated by loci clustered in relatively few regions of the genome, termed genomic islands of divergence. The mechanisms that create and maintain these islands and how they contribute to adaptive divergence is an active research topic. Here, we use sockeye salmon as a model to investigate both the mechanisms responsible for creating islands of divergence and the patterns of differentiation at these islands. Previous research suggested that multiple islands contributed to adaptive radiation of sockeye salmon. However, the low-density genomic methods used by these studies made it difficult to fully elucidate the mechanisms responsible for islands and connect genotypes to adaptive variation. We used whole genome resequencing to genotype millions of loci to investigate patterns of genetic variation at islands and the mechanisms that potentially created them. We discovered 64 islands, including 16 clustered in four genomic regions shared between two isolated populations. Characterisation of these four regions suggested that three were likely created by structural variation, while one was created by processes not involving structural variation. All four regions were small (< 600 kb), suggesting low recombination regions do not have to span megabases to be important for adaptive divergence. Differentiation at islands was not consistently associated with established population attributes. In sum, the landscape of adaptive divergence and the mechanisms that create it are complex; this complexity likely helps to facilitate fine-scale local adaptation unique to each population.

SNP Polymorphisms Are Associated with Environmental Factors in Sockeye Salmon Populations Across the Northwest Pacific: Insights from Redundancy Analysis

The SNP variation in sockeye salmon across the Asian part of its range was studied in 23 samples from 16 lake-river systems of the West Pacific Coast to improve understanding of genetic adaptation in response to spawning watersheds conditions. Identification of candidate SNPs and environmental factors that can contribute to local adaptations in sockeye salmon populations was carried out using redundancy analysis (RDA), a powerful tool for landscape genetics proven to be effective in genotype-environment association studies. Climatic and hydrographic indices (7 indices in total), reflecting abiotic conditions in freshwater habitats of sockeye salmon and characterizing the temperature regime in the river basin, its variability during the year, the amount of precipitation, as well as the height of the maximum tide in the estuary, were used as predictor factors. Among the 45 analyzed SNPs, several loci (ALDOB-135, HGFA, and RAG3-93) correlated with predictors gradients along the northwest Pacific coast were identified. The putative candidate loci localized in genes involved in the immune and inflammatory responses, as well as genes encoding temperature-sensitive enzymes and some hormones regulating ion homeostasis in fish during the anadromous migration and smoltification, were potentially associated with environmental conditions in natal rivers. The findings could have implications for aquaculture, conservation, and resource management in the context of global climate change.

Revealing the evolutionary history and contemporary population structure of Pacific salmon in the Fraser River through genome resequencing

The Fraser River once supported massive salmon returns. However, over the last century, the largest returns have consistently been less than half of the recorded historical maximum. There is substantial interest from surrounding communities and governments to increase salmon returns for both human use and functional ecosystems. To generate resources for this endeavor, we resequenced genomes of Chinook (Oncorhynchus tshawytscha), coho (Oncorhynchus kisutch), and sockeye salmon (Oncorhynchus nerka) from the Fraser River at moderate coverage (∼16×). A total of 954 resequenced genomes were analyzed, with 681 collected specifically for this study from tissues sampled between 1997 and 2021. An additional 273 were collected from previous studies. At the species level, Chinook salmon appeared to have 1.6-2.1× more SNPs than coho or sockeye salmon, respectively. This difference may be attributable to large historical declines of coho and sockeye salmon. At the population level, 3 Fraser River genetic groups were identified for each species using principal component and admixture analyses. These were consistent with previous research and supports the continued use of these groups in conservation and management efforts. Environmental factors and a migration barrier were identified as major factors influencing the boundaries of these genetic groups. Additionally, 20 potentially adaptive loci were identified among the genetic groups. This information may be valuable in new management and conservation efforts. Furthermore, the resequenced genomes are an important resource for contemporary genomics research on Fraser River salmon and have been made publicly available.

Genomic vulnerability of a freshwater salmonid under climate change

Understanding the adaptive potential of populations and species is pivotal for minimizing the loss of biodiversity in this era of rapid climate change. Adaptive potential has been estimated in various ways, including based on levels of standing genetic variation, presence of potentially beneficial alleles, and/or the severity of environmental change. Kokanee salmon, the non-migratory ecotype of sockeye salmon (Oncorhynchus nerka), is culturally and economically important and has already been impacted by the effects of climate change. To assess its climate vulnerability moving forward, we integrated analyses of standing genetic variation, genotype-environment associations, and climate modeling based on sequence and structural genomic variation from 224 whole genomes sampled from 22 lakes in British Columbia and Yukon (Canada). We found that variables for extreme temperatures, particularly warmer temperatures, had the most pervasive signature of selection in the genome and were the strongest predictors of levels of standing variation and of putatively adaptive genomic variation, both sequence and structural. Genomic offset estimates, a measure of climate vulnerability, were significantly correlated with higher increases in extreme warm temperatures, further highlighting the risk of summer heat waves that are predicted to increase in frequency in the future. Levels of standing genetic variation, an important metric for population viability and resilience, were not correlated with genomic offset. Nonetheless, our combined approach highlights the importance of integrating different sources of information and genomic data to formulate more comprehensive and accurate predictions on the vulnerability of populations and species to future climate change.

The late-evolving salmon and trout join the GnRH1 club

Although it is known that the whitefish, an ancient salmonid, expresses three distinct gonadotropin-releasing hormone (GnRH) forms in the brain, it has been thought that the later-evolving salmonids (salmon and trout) had only two types of GnRH: GnRH2 and GnRH3. We now provide evidence for the expression of GnRH1 in the gonads of Atlantic salmon by rapid amplification of cDNA ends, real-time quantitative PCR and immunohistochemistry. We examined six different salmonid genomes and found that each assembly has one gene that likely encodes a viable GnRH1 prepropeptide. In contrast to both functional GnRH2 and GnRH3 paralogs, the GnRH1 homeolog can no longer express the hormone. Furthermore, the viable salmonid GnRH1 mRNA is composed of only three exons, rather than the four exons that build the GnRH2 and GnRH3 mRNAs. Transcribed gnrh1 is broadly expressed (in 17/18 tissues examined), with relative abundance highest in the ovaries. Expression of the gnrh2 and gnrh3 mRNAs is more restricted, primarily to the brain, and not in the gonads. The GnRH1 proximal promoter presents composite binding elements that predict interactions with complexes that contain diverse cell fate and differentiation transcription factors. We provide immunological evidence for GnRH1 peptide in the nucleus of 1-year-old type A spermatogonia and cortical alveoli oocytes. GnRH1 peptide was not detected during other germ cell or reproductive stages. GnRH1 activity in the salmonid gonad may occur only during early stages of development and play a key role in a regulatory network that controls mitotic and/or meiotic processes within the germ cell.

Allele surfing causes maladaptation in a Pacific salmon of conservation concern

How various factors, including demography, recombination or genome duplication, may impact the efficacy of natural selection and the burden of deleterious mutations, is a central question in evolutionary biology and genetics. In this study, we show that key evolutionary processes, including variations in i) effective population size (Ne) ii) recombination rates and iii) chromosome inheritance, have influenced the genetic load and efficacy of selection in Coho salmon (Oncorhynchus kisutch), a widely distributed salmonid species on the west coast of North America. Using whole genome resequencing data from 14 populations at different migratory distances from their southern glacial refugium, we found evidence supporting gene surfing, wherein reduced Ne at the postglacial recolonization front, leads to a decrease in the efficacy of selection and a surf of deleterious alleles in the northernmost populations. Furthermore, our results indicate that recombination rates play a prime role in shaping the load along the genome. Additionally, we identified variation in polyploidy as a contributing factor to within-genome variation of the load. Overall, our results align remarkably well with expectations under the nearly neutral theory of molecular evolution. We discuss the fundamental and applied implications of these findings for evolutionary and conservation genomics.

Insights from a chum salmon (Oncorhynchus keta) genome assembly regarding whole-genome duplication and nucleotide variation influencing gene function

Chum salmon are ecologically important to Pacific Ocean ecosystems and commercially important to fisheries. To improve the genetic resources available for this species, we sequenced and assembled the genome of a male chum salmon using Oxford Nanopore read technology and the Flye genome assembly software (contig N50: ∼2 Mbp, complete BUSCOs: ∼98.1%). We also resequenced the genomes of 59 chum salmon from hatchery sources to better characterize the genome assembly and the diversity of nucleotide variants impacting phenotype variation. With genomic sequences from a doubled haploid individual, we were able to identify regions of the genome assembly that have been collapsed due to high sequence similarity between homeologous (duplicated) chromosomes. The homeologous chromosomes are relics of an ancient salmonid-specific genome duplication. These regions were enriched with genes whose functions are related to the immune system and responses to toxins. From analyzing nucleotide variant annotations of the resequenced genomes, we were also able to identify genes that have increased levels of variants thought to moderately impact gene function. Genes related to the immune system and the detection of chemical stimuli (olfaction) had increased levels of these variants based on a gene ontology enrichment analysis. The tandem organization of many of the enriched genes raises the question of why they have this organization.

A homology guide for Pacific salmon genus Oncorhynchus resolves patterns of ohnolog retention, resolution and local adaptation following the salmonid-specific whole-genome duplication event

Salmonid fishes have emerged as a tractable model to study whole-genome duplications (WGDs) as this group has undergone four rounds of WGDs. While most of the salmonid genome has returned to a diploid state, a significant proportion of genes are maintained as duplicates and are referred to as ohnologs. The fact that much of the modern salmonid gene repertoire is comprised of ohnologs, while other genes have returned to their singleton state creates complications for genetic studies by obscuring homology relationships. The difficulty this creates is particularly prominent in Pacific salmonids belonging to genus Oncorhynchus who are the focus of intense genetics-based conservation and management efforts owing to the important ecological and cultural roles these fish play. To address this gap, we generated a homology guide for six species of Oncorhynchus with available genomes and used this guide to describe patterns of ohnolog retention and resolution. Overall, we find that ohnologs comprise approximately half of each species modern gene repertoires, which are functionally enriched for genes involved in DNA binding, while the less numerous singleton genes are heavily enriched in dosage-sensitive processes such as mitochondrial metabolism. Additionally, by reanalyzing published expression data from locally adapted strains of O. mykiss, we show that numerous ohnologs exhibit adaptive expression profiles; however, ohnologs are not more likely to display adaptive signatures than either paralogs or singletons. Finally, we demonstrate the utility of our homology guide by investigating the evolutionary relationship among genes highlighted as playing a role in salmonid life-history traits or gene editing targets.

The genomic basis of reproductive and migratory behaviour in a polymorphic salmonid

Recent ecotypic differentiation provides unique opportunities to investigate the genomic basis and architecture of local adaptation, while offering insights into how species form and persist. Sockeye salmon (Oncorhynchus nerka) exhibit migratory and resident ("kokanee") ecotypes, which are further distinguished into shore-spawning and stream-spawning reproductive ecotypes. Here, we analysed 36 sockeye (stream-spawning) and kokanee (stream- and shore-spawning) genomes from a system where they co-occur and have recent common ancestry (Okanagan Lake/River in British Columbia, Canada) to investigate the genomic basis of reproductive and migratory behaviour. Examination of the genomic landscape of differentiation, differences in allele frequencies and genotype-phenotype associations revealed three main blocks of sequence differentiation on chromosomes 7, 12 and 20, associated with migratory behaviour, spawning location and spawning timing. Structural variants identified in these same areas suggest they could contribute to ecotypic differentiation directly as causal variants or via maintenance of their genomic architecture through recombination suppression mechanisms. Genes in these regions were related to spatial memory and swimming endurance (SYNGAP, TPM3), as well as eye and brain development (including SIX6), potentially associated with differences in migratory behaviour and visual habitats across spawning locations, respectively. Additional genes (GREB1L, ROCK1) identified here have been associated with timing of migration in other salmonids and could explain variation in timing of O. nerka spawning. Together, these results based on the joint analysis of sequence and structural variation represent a significant advance in our understanding of the genomic landscape of ecotypic differentiation at different stages in the speciation continuum.

A chromosome-level genome of Brachymystax tsinlingensis provides resources and insights into salmonids evolution

Brachymystax tsinlingensis Li, 1966 is an endangered freshwater fish with economic, ecological, and scientific values. Study of the genome of B. tsinlingensis might be particularly insightful given that this is the only Brachymystax species with genome. We present a high-quality chromosome-level genome assembly and protein-coding gene annotation for B. tsinlingensis with Illumina short reads, Nanopore long reads, Hi-C sequencing reads, and RNA-seq reads from 5 tissues/organs. The final chromosome-level genome size is 2,031,709,341 bp with 40 chromosomes. We found that the salmonids have a unique GC content and codon usage, have a slower evolutionary rate, and possess specific positively selected genes. We also confirmed the salmonids have undergone a whole-genome duplication event and a burst of transposon-mediated repeat expansion, and lost HoxAbβ Hox cluster, highly expressed genes in muscle may partially explain the migratory habits of B. tsinlingensis. The high-quality B. tsinlingensis assembled genome could provide a valuable reference for the study of other salmonids as well as aid the conservation of this endangered species.

An update of the salmon louse (Lepeophtheirus salmonis) reference genome assembly

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.

Genotyping-in-Thousands by sequencing panel development and application to inform kokanee salmon (Oncorhynchus nerka) fisheries management at multiple scales

The ability to differentiate life history variants is vital for estimating fisheries management parameters, yet traditional survey methods can be inaccurate in mixed-stock fisheries. Such is the case for kokanee, the freshwater resident form of sockeye salmon (Oncorhynchus nerka), which exhibits various reproductive ecotypes (stream-, shore-, deep-spawning) that co-occur with each other and/or anadromous O. nerka in some systems across their pan-Pacific distribution. Here, we developed a multi-purpose Genotyping-in-Thousands by sequencing (GT-seq) panel of 288 targeted single nucleotide polymorphisms (SNPs) to enable accurate kokanee stock identification by geographic basin, migratory form, and reproductive ecotype across British Columbia, Canada. The GT-seq panel exhibited high self-assignment accuracy (93.3%) and perfect assignment of individuals not included in the baseline to their geographic basin, migratory form, and reproductive ecotype of origin. The GT-seq panel was subsequently applied to Wood Lake, a valuable mixed-stock fishery, revealing high concordance (>98%) with previous assignments to ecotype using microsatellites and TaqMan® SNP genotyping assays, while improving resolution, extending a long-term time-series, and demonstrating the scalability of this approach for this system and others.

The pink salmon genome: Uncovering the genomic consequences of a two-year life cycle

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.

Genome-Wide Investigation of the Multiple Origins Hypothesis for Deep-Spawning Kokanee Salmon (Oncorhynchus nerka) across its Pan-Pacific Distribution

Salmonids have emerged as important study systems for investigating molecular processes underlying parallel evolution given their tremendous life history variation. Kokanee, the resident form of anadromous sockeye salmon (Oncorhynchus nerka), have evolved multiple times across the species' pan-Pacific distribution, exhibiting multiple reproductive ecotypes including those that spawn in streams, on lake-shores, and at lake depths >50 m. The latter has only been detected in 5 locations in Japan and British Columbia, Canada. Here, we investigated the multiple origins hypothesis for deep-spawning kokanee, using 9721 single nucleotide polymorphisms distributed across the genome analyzed for the vast majority of known populations in Japan (Saiko Lake) and Canada (Anderson, Seton, East Barrière Lakes) relative to stream-spawning populations in both regions. We detected 397 outlier loci, none of which were robustly identified in paired-ecotype comparisons in Japan and Canada independently. Bayesian clustering and principal components analyses based on neutral loci revealed 6 distinct clusters, largely associated with geography or translocation history, rather than ecotype. Moreover, a high level of divergence between Canadian and Japanese populations, and between deep- and stream-spawning populations regionally, suggests the deep-spawning ecotype independently evolved on the 2 continents. On a finer level, Japanese kokanee populations exhibited low estimates of heterozygosity, significant levels of inbreeding, and reduced effective population sizes relative to Canadian populations, likely associated with transplantation history. Along with preliminary evidence for hybridization between deep- and stream-spawning ecotypes in Saiko Lake, these findings should be considered within the context of on-going kokanee fisheries management in Japan.

Genome-wide analysis reveals demographic and life-history patterns associated with habitat modification in landlocked, deep-spawning sockeye salmon (Oncorhynchus nerka)

Human-mediated habitat fragmentation in freshwater ecosystems can negatively impact genetic diversity, demography, and life history of native biota, while disrupting the behavior of species that are dependent on spatial connectivity to complete their life cycles. In the Alouette River system (British Columbia, Canada), dam construction in 1928 impacted passage of anadromous sockeye salmon (Oncorhynchus nerka), with the last records of migrants occurring in the 1930s. Since that time, O. nerka persisted as a resident population in Alouette Reservoir until experimental water releases beginning in 2005 created conditions for migration; two years later, returning migrants were observed for the first time in ~70 years, raising important basic and applied questions regarding life-history variation and population structure in this system. Here, we investigated the genetic distinctiveness and population history of Alouette Reservoir O. nerka using genome-wide SNP data (n = 7,709 loci) collected for resident and migrant individuals, as well as for neighboring anadromous sockeye salmon and resident kokanee populations within the Fraser River drainage (n = 312 individuals). Bayesian clustering and principal components analyses based on neutral loci revealed five distinct clusters, largely associated with geography, and clearly demonstrated that Alouette Reservoir resident and migrant individuals are genetically distinct from other O. nerka populations in the Fraser River drainage. At a finer level, there was no clear evidence for differentiation between Alouette Reservoir residents and migrants; although we detected eight high-confidence outlier loci, they all mapped to sex chromosomes suggesting that differences were likely due to uneven sex ratios rather than life history. Taken together, these data suggest that contemporary Alouette Reservoir O. nerka represents a landlocked sockeye salmon population, constituting the first reported instance of deep-water spawning behavior associated with this life-history form. This finding punctuates the need for reassessment of conservation status and supports ongoing fisheries management activities in Alouette Reservoir.

The use of non-lethal sampling for transcriptomics to assess the physiological status of wild fishes

Fishes respond to different abiotic and biotic stressors through changes in gene expression as a part of an integrated physiological response. Transcriptomics approaches have been used to quantify gene expression patterns as a reductionist approach to understand responses to environmental stressors in animal physiology and have become more commonly used to study wild fishes. We argue that non-lethal sampling for transcriptomics should become the norm for assessing the physiological status of wild fishes, especially when there are conservation implications. Processes at the level of the transcriptome provide a "snapshot" of the cellular conditions at a given time; however, by using a non-lethal sampling protocol, researchers can connect the transcriptome profile with fitness-relevant ecological endpoints such as reproduction, movement patterns and survival. Furthermore, telemetry is a widely used approach in fisheries to understand movement patterns in the wild, and when combined with transcriptional profiling, provides arguably the most powerful use of non-lethal sampling for transcriptomics in wild fishes. In this review, we discuss the different tissues that can be successfully incorporated into non-lethal sampling strategies, which is particularly useful in the context of the emerging field of conservation transcriptomics. We briefly describe different methods for transcriptional profiling in fishes from high-throughput qPCR to whole transcriptome approaches. Further, we discuss strategies and the limitations of using transcriptomics for non-lethally studying fishes. Lastly, as 'omics' technology continues to advance, transcriptomics paired with different omics approaches to study wild fishes will provide insight into the factors that regulate phenotypic variation and the physiological responses to changing environmental conditions in the future.

EigenGWAS: An online visualizing and interactive application for detecting genomic signatures of natural selection

Detecting genetic regions under selection in structured populations is of great importance in ecology, evolutionary biology and breeding programmes. We recently proposed EigenGWAS, an unsupervised genomic scanning approach that is similar to F ST but does not require grouping information of the population, for detection of genomic regions under selection. The original EigenGWAS is designed for the random mating population, and here we extend its use to inbred populations. We also show in theory and simulation that eigenvalues, the previous corrector for genetic drift in EigenGWAS, are overcorrected for genetic drift, and the genomic inflation factor is a better option for this adjustment. Applying the updated algorithm, we introduce the new EigenGWAS online platform with highly efficient core implementation. Our online computational tool accepts plink data in a standard binary format that can be easily converted from the original sequencing data, provides the users with graphical results via the R-Shiny user-friendly interface. We applied the proposed method and tool to various data sets, and biologically interpretable results as well as caveats that may lead to an unsatisfactory outcome are given. The EigenGWAS online platform is available at www.eigengwas.com, and can be localized and scaled up via R (recommended) or docker.