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.

Population-size history inferences from the coho salmon (Oncorhynchus kisutch) genome.

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.

Long-distance migration is a major factor driving local adaptation at continental scale in Coho salmon.

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.

Re-evaluating Coho salmon (Oncorhynchus kisutch) conservation units in Canada using genomic data.

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 (Oncorhynchus kisutch) from 149 sampling locations and a total of 24,542 high-quality SNPs obtained using genotyping-by-sequencing and mapped to the Coho salmon reference genome to (1) evaluate the current delineation of CUs for Coho in Canada and (2) compare patterns of population structure observed using neutral and outlier loci from genotype-environment association analyses to determine whether separate CUs that capture adaptive diversity are needed. Our results reflected CU boundaries on the whole, with the majority of sampling locations managed in the same CU clustering together within genetic groups. However, additional groups that are not currently represented by CUs were also uncovered. We observed considerable overlap in the genetic clusters identified using neutral or candidate loci, indicating a general congruence in patterns of genetic variation driven by local adaptation and gene flow in this species. Consequently, we suggest that the current CU boundaries for Coho salmon are largely well-suited for meeting the Canadian Wild Salmon Policy's objective of defining biologically distinct groups, but we highlight specific areas where CU boundaries may be refined.

Strong Parallel Differential Gene Expression Induced by Hatchery Rearing Weakly Associated with Methylation Signals in Adult Coho Salmon (O. kisutch).

Human activities and resource exploitation led to a massive decline of wild salmonid populations, consequently, numerous conservation programs have been developed to supplement wild populations. However, many studies documented reduced fitness of hatchery-born relative to wild fish. Here, by using both RNA sequencing and Whole Genome Bisulfite Sequencing of hatchery and wild-born adult Coho salmon (Oncorhynchus kisutch) originating from two previously studied river systems, we show that early-life hatchery-rearing environment-induced significant and parallel gene expression differentiation is maintained until Coho come back to their natal river for reproduction. A total of 3,643 genes differentially expressed and 859 coexpressed genes were downregulated in parallel in hatchery-born fish from both rivers relative to their wild congeners. Among those genes, 26 displayed a significant relationship between gene expression and the median gene body methylation and 669 single CpGs displayed a significant correlation between methylation level and the associated gene expression. The link between methylation and gene expression was weak suggesting that DNA methylation is not the only player in mediating hatchery-related expression differences. Yet, significant gene expression differentiation was observed despite 18 months spent in a common environment (i.e., the sea). Finally, the differentiation is observed in parallel in two different river systems, highlighting the fact that early-life environment may account for at least some of the reduced fitness of the hatchery salmon in the wild. These results illustrate the relevance and importance of considering both epigenome and transcriptome to evaluate the costs and benefits of large-scale supplementation programs.

Chinook and Coho salmon hybrids linked to habitat and climatic changes on Vancouver Island, British Columbia.

Between 2013 and 2019, 63 presumed Chinook salmon Oncorhynchus tshawytscha sampled primarily in the Strait of Georgia (0.63% of total sample) were identified as potential Chinook-Coho (Oncorhynchus kisutch) hybrids by the presence of anomalous microsatellite genotypes. Their hybrid origin was confirmed by single nucleotide polymorphism amplification of two species-specific amplicons. Mitochondrial DNA indicated that most of these fish resulted from the hybridization of Coho salmon females and Chinook salmon males. Although no diagnostic external features were identified, several individuals displayed an abnormal scale arrangement on the caudal peduncle. One hybrid juvenile examined for meristics exhibited a pyloric caeca count intermediate between published values for Chinook and Coho salmon. Most hybrids originated in the Cowichan River during the 2014 brood year. Their prevalence in the watershed is a naturally occurring event, likely exacerbated by prolonged low water levels which limit habitat and delay Chinook salmon spawning, in addition to the differential abundance of the parental species. This research is the first to document ongoing natural hybridization (Chinook-Coho salmon crosses) and link it to habitat and climatic changes, and includes the identification of eight F1 adults and two juvenile backcross or F2 hybrids. The potential negative impacts of hybridization, particularly in Coho salmon through potential introgression, warrant hybrid identification as an ecosystem monitoring tool within a survey program.

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

[This corrects the article DOI: 10.1371/journal.pone.0240935.].

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.

Epigenomic modifications induced by hatchery rearing persist in germ line cells of adult salmon after their oceanic migration.

Human activities induce direct or indirect selection pressure on natural population and may ultimately affect population's integrity. While numerous conservation programs aimed to minimize human-induced genomic variation, human-induced environmental variation may generate epigenomic variation potentially affecting fitness through phenotypic modifications. Major questions remain pertaining to how much epigenomic variation arises from environmental heterogeneity, whether this variation can persist throughout life, and whether it can be transmitted across generations. We performed whole genome bisulfite sequencing (WGBS) on the sperm of genetically indistinguishable hatchery and wild-born migrating adults of Coho salmon (Oncorhynchus kisutch) from two geographically distant rivers at different epigenome scales. Our results showed that coupling WGBS with fine-scale analyses (local and chromosomal) allows the detection of parallel early-life hatchery-induced epimarks that differentiate wild from hatchery-reared salmon. Four chromosomes and 183 differentially methylated regions (DMRs) displayed a significant signal of methylation differentiation between hatchery and wild-born Coho salmon. Moreover, those early-life epimarks persisted in germ line cells despite about 1.5 year spent in the ocean following release from hatchery, opening the possibility for transgenerational inheritance. Our results strengthen the hypothesis that epigenomic modifications environmentally induced during early-life development persist in germ cells of adults until reproduction, which could potentially impact their fitness.

Correction: Demographic history shaped geographical patterns of deleterious mutation load in a broadly distributed Pacific Salmon.

[This corrects the article DOI: 10.1371/journal.pgen.1008348.].

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

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.

Demographic history shaped geographical patterns of deleterious mutation load in a broadly distributed Pacific 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.

Insights on the concept of indicator populations derived from parentage-based tagging in a large-scale coho salmon application in British Columbia, Canada.

For Pacific salmon, the key fisheries management goal in British Columbia (BC) is to maintain and restore healthy and diverse Pacific salmon populations, making conservation of salmon biodiversity the highest priority for resource management decision-making. Salmon status assessments are often conducted on coded-wire-tagged subsets of indicator populations based on assumptions of little differentiation within or among proximal populations. In the current study of southern BC coho salmon (Oncorhynchus kisutch) populations, parentage-based tagging (PBT) analysis provided novel information on migration and life-history patterns to test the assumptions of biological homogeneity over limited (generally < 100 km) geographic distances and, potentially, to inform management of fisheries and hatchery broodstocks. Heterogeneity for location and timing of fishery captures, family productivity, and exploitation rate was observed over small geographic scales, within regions that are, or might be expected to be, within the area encompassed by a single-tagged indicator population. These results provide little support for the suggestion that information gained from tagged indicator populations is representative of marine distribution, productivity, and exploitation patterns of proximal populations.

Variation in migration pattern, broodstock origin, and family productivity of coho salmon hatchery populations in British Columbia, Canada, derived from parentage-based tagging.

In salmonid parentage-based tagging (PBT) applications, entire hatchery broodstocks are genotyped, and subsequently, progeny can be nonlethally sampled and assigned back to their parents using parentage analysis, thus identifying their hatchery of origin and brood year (i.e., age). Inter- and intrapopulation variability in migration patterns, life history traits, and fishery contributions can be determined from PBT analysis of samples derived from both fisheries and escapements (portion of a salmon population that does not get caught in fisheries and returns to its natal river to spawn). In the current study of southern British Columbia coho salmon (Oncorhynchus kisutch) populations, PBT analysis provided novel information on intrapopulation heterogeneity among males in the total number of progeny identified in fisheries and escapements, the proportion of progeny sampled from fisheries versus escapement, the proportion of two-year-old progeny (jacks) produced, and the within-season return time of progeny. Fishery recoveries of coho salmon revealed heterogeneity in migration patterns among and within populations, with recoveries from north and central coast fisheries distinguishing "northern migrating" from "resident" populations. In northern migrating populations, the mean distance between fishery captures of sibs (brothers and sisters) was significantly less than the mean distance between nonsibs, indicating the possible presence of intrapopulation genetic heterogeneity for migration pattern. Variation among populations in productivity and within populations in fish catchability indicated that population selection and broodstock management can be implemented to optimize harvest benefits from hatcheries. Application of PBT provided valuable information for assessment and management of hatchery-origin coho salmon in British Columbia.

Comparison of coded-wire tagging with parentage-based tagging and genetic stock identification in a large-scale coho salmon fisheries application in British Columbia, Canada.

Wild Pacific salmon, including Coho salmon Onchorynchus kisutch, have been supplemented with hatchery propagation for over 50 years in support of increased ocean harvest and conservation of threatened populations. In Canada, the Wild Salmon Policy for Pacific salmon was established with the goal of maintaining and restoring healthy and diverse Pacific salmon populations, making conservation of wild salmon and their habitats the highest priority for resource management decision-making. A new approach to the assessment and management of wild coho salmon, and the associated hatchery production and fishery management is needed. Implementation of parentage-based tagging (PBT) may overcome problems associated with coded-wire tag-based (CWT) assessment and management of coho salmon fisheries, providing at a minimum information equivalent to that derived from the CWT program. PBT and genetic stock identification (GSI) were used to identify coho salmon sampled in fisheries (8,006 individuals) and escapements (1,692 individuals) in British Columbia to specific conservation units (CU), populations, and broodyears. Individuals were genotyped at 304 single nucleotide polymorphisms (SNPs) via direct sequencing of amplicons. Very high accuracy of assignment to population (100%) via PBT for 543 jack (age 2) assigned to correct age and collection location and 265 coded-wire tag (CWT, age 3) coho salmon assigned to correct age and release location was observed, with a 40,774-individual, 267-population baseline available for assignment. Coho salmon from un-CWTed enhanced populations contributed 65% of the catch in southern recreational fisheries in 2017. Application of a PBT-GSI system of identification to individuals in 2017 fisheries and escapements provided high-resolution estimates of stock composition, catch, and exploitation rate by CU or population, providing an alternate and more effective method in the assessment and management of Canadian-origin coho salmon relative to CWTs, and an opportunity for a genetic-based system to replace the current CWT system for coho salmon assessment.

Chinook salmon (Oncorhynchus tshawytscha) genome and transcriptome.

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

Parallel epigenetic modifications induced by hatchery rearing in a Pacific salmon.

Wild stocks of Pacific salmonids have experienced sharp declines in abundance over the past century. Consequently, billions of fish are released each year for enhancing abundance and sustaining fisheries. However, the beneficial role of this widely used management practice is highly debated since fitness decrease of hatchery-origin fish in the wild has been documented. Artificial selection in hatcheries has often been invoked as the most likely explanation for reduced fitness, and most studies to date have focused on finding signatures of hatchery-induced selection at the DNA level. We tested an alternative hypothesis, that captive rearing induces epigenetic reprogramming, by comparing genome-wide patterns of methylation and variation at the DNA level in hatchery-reared coho salmon (Oncorhynchus kisutch) with those of their wild counterparts in two geographically distant rivers. We found a highly significant proportion of epigenetic variation explained by the rearing environment that was as high as the one explained by the river of origin. The differentially methylated regions show enrichment for biological functions that may affect the capacity of hatchery-born smolts to migrate successfully in the ocean. Shared epigenetic variation between hatchery-reared salmon provides evidence for parallel epigenetic modifications induced by hatchery rearing in the absence of genetic differentiation between hatchery and natural-origin fish for each river. This study highlights epigenetic modifications induced by captive rearing as a potential explanatory mechanism for reduced fitness in hatchery-reared salmon.

Population structure of sea-type and lake-type sockeye salmon and kokanee in the Fraser River and Columbia River drainages.

Population structure of three ecotypes of Oncorhynchus nerka (sea-type Sockeye Salmon, lake-type Sockeye Salmon, and Kokanee) in the Fraser River and Columbia River drainages was examined with microsatellite variation, with the main focus as to whether Kokanee population structure within the Fraser River drainage suggested either a monophyletic or polyphyletic origin of the ecotype within the drainage. Variation at 14 microsatellite loci was surveyed for sea-type and lake-type Sockeye Salmon and Kokanee sampled from 121 populations in the two river drainages. An index of genetic differentiation, FST, over all populations and loci was 0.087, with individual locus values ranging from 0.031 to 0.172. Standardized to an ecotype sample size of 275 individuals, the least genetically diverse ecotype was sea-type Sockeye Salmon with 203 alleles, whereas Kokanee displayed the greatest number of alleles (260 alleles), with lake-type Sockeye Salmon intermediate (241 alleles). Kokanee populations from the Columbia River drainage (Okanagan Lake, Kootenay Lake), the South Thompson River (a major Fraser River tributary) drainage populations, and the mid-Fraser River populations all clustered together in a neighbor-joining analysis, indicative of a monophyletic origin of the Kokanee ecotype in these regions, likely reflecting the origin of salmon radiating from a refuge after the last glaciation period. However, upstream of the mid-Fraser River populations, there were closer relationships between the lake-type Sockeye Salmon ecotype and the Kokanee ecotype, indicative of the Kokanee ecotype evolving independently from the lake-type Sockeye Salmon ecotype in parallel radiation. Kokanee population structure within the entire Fraser River drainage suggested a polyphyletic origin of the ecotype within the drainage. Studies employing geographically restricted population sampling may not outline accurately the phylogenetic history of salmonid ecotypes.

Harnessing the Power of Genomics to Secure the Future of Seafood.

Best use of scientific knowledge is required to maintain the fundamental role of seafood in human nutrition. While it is acknowledged that genomic-based methods allow the collection of powerful data, their value to inform fisheries management, aquaculture, and biosecurity applications remains underestimated. We review genomic applications of relevance to the sustainable management of seafood resources, illustrate the benefits of, and identify barriers to their integration. We conclude that the value of genomic information towards securing the future of seafood does not need to be further demonstrated. Instead, we need immediate efforts to remove structural roadblocks and focus on ways that support integration of genomic-informed methods into management and production practices. We propose solutions to pave the way forward.

Detection of quantitative trait loci (QTL) related to grilsing and late sexual maturation in Atlantic salmon (Salmo salar).

In Atlantic salmon aquaculture, early sexual maturation represents a major problem for producers. This is especially true for grilse, which mature after one sea winter before reaching a desirable harvest weight, rather than after two sea winters. Salmon maturing as grilse have a much lower market value than later maturing individuals. For this reason, most companies desire fish that grow fast and mature late. Marker-assisted selection has the potential to improve the efficiency of selection against early maturation and for late sexual maturation; however, studies identifying age of sexual maturation-related genetic markers are lacking for Atlantic salmon. Therefore, we used a 6.5K single-nucleotide polymorphism (SNP) array to genotype five families from the Mainstream Canada broodstock program and search for SNPs associated with early (grilsing) or late sexual maturation. There were 529 SNP loci that were variable across all five families, and this was the set that was used for quantitative trait loci (QTL) analysis. GridQTL identified two chromosomes, Ssa10 and Ssa21, containing QTL related to grilsing. In contrast, only one QTL, on Ssa18, was found linked to late maturation in Atlantic salmon. Our previous work on these five families did not identify genome-wide significant growth-related QTL on Ssa10, Ssa21, or Ssa18. Therefore, taken together, these results suggest that both grilsing and late sexual maturation are controlled independently of one another and also from growth-related traits. The identification of genomic regions associated with grilsing or late sexual maturation provide an opportunity to incorporate this information into selective breeding programs that will enhance Atlantic salmon farming.