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.

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

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.

Parallelism in eco-morphology and gene expression despite variable evolutionary and genomic backgrounds in a Holarctic fish.

Understanding the extent to which ecological divergence is repeatable is essential for predicting responses of biodiversity to environmental change. Here we test the predictability of evolution, from genotype to phenotype, by studying parallel evolution in a salmonid fish, Arctic charr (Salvelinus alpinus), across eleven replicate sympatric ecotype pairs (benthivorous-planktivorous and planktivorous-piscivorous) and two evolutionary lineages. We found considerable variability in eco-morphological divergence, with several traits related to foraging (eye diameter, pectoral fin length) being highly parallel even across lineages. This suggests repeated and predictable adaptation to environment. Consistent with ancestral genetic variation, hundreds of loci were associated with ecotype divergence within lineages of which eight were shared across lineages. This shared genetic variation was maintained despite variation in evolutionary histories, ranging from postglacial divergence in sympatry (ca. 10-15kya) to pre-glacial divergence (ca. 20-40kya) with postglacial secondary contact. Transcriptome-wide gene expression (44,102 genes) was highly parallel across replicates, involved biological processes characteristic of ecotype morphology and physiology, and revealed parallelism at the level of regulatory networks. This expression divergence was not only plastic but in part genetically controlled by parallel cis-eQTL. Lastly, we found that the magnitude of phenotypic divergence was largely correlated with the genetic differentiation and gene expression divergence. In contrast, the direction of phenotypic change was mostly determined by the interplay of adaptive genetic variation, gene expression, and ecosystem size. Ecosystem size further explained variation in putatively adaptive, ecotype-associated genomic patterns within and across lineages, highlighting the role of environmental variation and stochasticity in parallel evolution. Together, our findings demonstrate the parallel evolution of eco-morphology and gene expression within and across evolutionary lineages, which is controlled by the interplay of environmental stochasticity and evolutionary contingencies, largely overcoming variable evolutionary histories and genomic backgrounds.

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.

The Atlantic salmon genome provides insights into rediploidization.

The whole-genome duplication 80 million years ago of the common ancestor of salmonids (salmonid-specific fourth vertebrate whole-genome duplication, Ss4R) provides unique opportunities to learn about the evolutionary fate of a duplicated vertebrate genome in 70 extant lineages. Here we present a high-quality genome assembly for Atlantic salmon (Salmo salar), and show that large genomic reorganizations, coinciding with bursts of transposon-mediated repeat expansions, were crucial for the post-Ss4R rediploidization process. Comparisons of duplicate gene expression patterns across a wide range of tissues with orthologous genes from a pre-Ss4R outgroup unexpectedly demonstrate far more instances of neofunctionalization than subfunctionalization. Surprisingly, we find that genes that were retained as duplicates after the teleost-specific whole-genome duplication 320 million years ago were not more likely to be retained after the Ss4R, and that the duplicate retention was not influenced to a great extent by the nature of the predicted protein interactions of the gene products. Finally, we demonstrate that the Atlantic salmon assembly can serve as a reference sequence for the study of other salmonids for a range of purposes.

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.

Effect of triploidy on liver gene expression in coho salmon (Oncorhynchus kisutch) under different metabolic states.

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.

A genetic linkage map for the salmon louse (Lepeophtheirus salmonis): evidence for high male:female and inter-familial recombination rate differences.

A salmon louse (Lepeophtheirus salmonis salmonis) genetic linkage map was constructed to serve as a genomic resource for future investigations into the biology of this important marine parasitic copepod species, and to provide insights into the inheritance patterns of genetic markers in this species. SNP genotyping of 8 families confirmed the presence of 15 linkage groups based upon the assignment of 93,773 markers. Progeny sample size weight adjusted map sizes in males (with the exception of SL12 and SL15) ranged in size from 96.50 cM (SL11) to 134.61 cM (SL06), and total combined map steps or bins ranged from 143 (SL09) to 203 (SL13). The SL12 male map was the smallest linkage group with a weight-averaged size of 3.05 cM with 6 recombination bins. Male:female specific recombination rate differences are 10.49:1 and represent one of the largest reported sex-specific differences for any animal species. Recombination ratio differences (M:F) ranged from 1.0 (SL12) to 29:1 (SL15). The number of markers exhibiting normal Mendelian segregation within the sex linkage group SL15 was extremely low (N = 80) in comparison to other linkage groups genotyped [range: 1459 (SL12)-10206 markers (SL05)]. Re-evaluation of Mendelian inheritance patterns of markers unassigned to any mapping parent according to hemizygous segregation patterns (models presented) identified matches for many of these markers to hemizygous patterns. The greatest proportion of these markers assigned to SL15 (N increased to 574). Inclusion of the hemizygous markers revised SL15 sex-specific recombination rate differences to 28:1. Recombination hot- and coldspots were identified across all linkage groups with all linkage groups possessing multiple peaks. Nine of 13 linkage groups evaluated possessed adjacent domains with hot-coldspot transitional zones. The most common pattern was for one end of the linkage to show elevated recombination in addition to internal regions. For SL01 and SL06, however, a terminal region with high recombination was not evident while a central domain possessing extremely high-recombination levels was present. High levels of recombination were weakly coupled to higher levels of SNP variation within domains, but this association was very strong for the central domains of SL01 and SL06. From the pooled paternal half-sib lots (several virgin females placed with 1 male), only 1 or two surviving family lots were obtained. Surviving families possessed parents where both the male and female possessed either inherently low or high recombination rates. This study provides insight into the organization of the sea louse genome, and describes large differences in recombination rate that exist among individuals of the same sex, and between the sexes. These differences in recombination rate may be coupled to the capabilities of this species to adapt to environmental and pharmaceutical treatments, given that family survivorship appears to be enhanced when parents have similar recombination levels.

Subcellular localization and characterization of estrogenic pathway regulators and mediators in Atlantic salmon spermatozoal cells.

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

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 (O. kisutch), and sockeye salmon (O. nerka) from the Fraser River at moderate coverage (∼16x). 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.1x 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, three Fraser River genetic groups were identified for each species using principal component and admixture analyses, which is 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.

Genomics of sablefish (Anoplopoma fimbria): expressed genes, mitochondrial phylogeny, linkage map and identification of a putative sex gene.

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

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.

In-field genetic stock identification of overwintering coho salmon in the Gulf of Alaska: Evaluation of Nanopore sequencing for remote real-time deployment.

Genetic stock identification (GSI) from genotyping-by-sequencing of single nucleotide polymorphism (SNP) loci has become the gold standard for stock of origin identification in Pacific salmon. The sequencing platforms currently applied require large batch sizes and multiday processing in specialized facilities to perform genotyping by the thousands. However, recent advances in third-generation single-molecule sequencing platforms, such as the Oxford Nanopore minION, provide base calling on portable, pocket-sized sequencers and promise real-time, in-field stock identification of variable batch sizes. Here we evaluate utility and comparability to established GSI platforms of at-sea stock identification of coho salmon (Oncorhynchus kisutch) using targeted SNP amplicon sequencing on the minION platform during a high-sea winter expedition to the Gulf of Alaska. As long read sequencers are not optimized for short amplicons, we concatenate amplicons to increase coverage and throughput. Nanopore sequencing at-sea yielded data sufficient for stock assignment for 50 out of 80 individuals. Nanopore-based SNP calls agreed with Ion Torrent-based genotypes in 83.25%, but assignment of individuals to stock of origin only agreed in 61.5% of individuals, highlighting inherent challenges of Nanopore sequencing, such as resolution of homopolymer tracts and indels. However, poor representation of assayed salmon in the queried baseline data set contributed to poor assignment confidence on both platforms. Future improvements will focus on lowering turnaround time and cost, increasing accuracy and throughput, as well as augmentation of the existing baselines. If successfully implemented, Nanopore sequencing will provide an alternative method to the large-scale laboratory approach by providing mobile small batch genotyping to diverse stakeholders.

Regulation and expression of sexual differentiation factors in embryonic and extragonadal tissues of Atlantic salmon.

BACKGROUND: The products of cyp19, dax, foxl2, mis, sf1 and sox9 have each been associated with sex-determining processes among vertebrates. We provide evidence for expression of these regulators very early in salmonid development and in tissues outside of the hypothalamic-pituitary-adrenal/gonadal (HPAG) axis. Although the function of these factors in sexual differentiation have been defined, their roles in early development before sexual fate decisions and in tissues beyond the brain or gonad are essentially unknown. RESULTS: Bacterial artificial chromosomes containing salmon dax1 and dax2, foxl2b and mis were isolated and the regulatory regions that control their expression were characterized. Transposon integrations are implicated in the shaping of the dax and foxl2 loci. Splice variants for cyp19b1 and mis in both embryonic and adult tissues were detected and characterized. We found that cyp19b1 transcripts are generated that contain 5'-untranslated regions of different lengths due to cryptic splicing of the 3'-end of intron 1. We also demonstrate that salmon mis transcripts can encode prodomain products that present different C-termini and terminate before translation of the MIS hormone. Regulatory differences in the expression of two distinct aromatases cyp19a and cyp19b1 are exerted, despite transcription of their transactivators (ie; dax1, foxl2, sf1) occurring much earlier during embryonic development. CONCLUSIONS: We report the embryonic and extragonadal expression of dax, foxl2, mis and other differentiation factors that indicate that they have functions that are more general and not restricted to steroidogenesis and gonadogenesis. Spliced cyp19b1 and mis transcripts are generated that may provide regulatory controls for tissue- or development-specific activities. Selection of cyp19b1 transcripts may be regulated by DAX-1, FOXL2 and SF-1 complexes that bind motifs in intron 1, or by signals within exon 2 that recruit splicing factors, or both. The potential translation of proteins bearing only the N-terminal MIS prodomain may modulate the functions of other TGF ß family members in different tissues. The expression patterns of dax1 early in salmon embryogenesis implicate its role as a lineage determination factor. Other roles for these factors during embryogenesis and outside the HPAG axis are discussed.

Assessing the effects of genotype-by-environment interaction on epigenetic, transcriptomic, and phenotypic response in a Pacific salmon.

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.

Parentage-based tagging combined with genetic stock identification is a cost-effective and viable replacement for coded-wire tagging in large-scale assessments of marine Chinook salmon fisheries in British Columbia, Canada.

Wild Pacific salmon, including Chinook salmon Oncorhynchus tshawytscha, have been supplemented with hatchery propagation for over 50 years in support of increased ocean harvest, mitigation for hydroelectric development, 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. For policy implementation, a new approach to the assessment and management of Chinook salmon and the associated hatchery production and fisheries management are needed. Implementation of genetic stock identification (GSI) and parentage-based tagging (PBT) for marine fisheries assessment may overcome problems associated with coded-wire tag-based (CWT) assessment and management of Chinook salmon fisheries, providing at a minimum information equivalent to that derived from the CWT program. GSI and PBT were used to identify Chinook salmon sampled in 2018 and 2019 marine fisheries (18,819 individuals genotyped) in British Columbia to specific conservation units (CU), populations, and broodyears. Individuals were genotyped at 391 single nucleotide polymorphisms via direct sequencing of amplicons. Very high accuracy of assignment to population and age (>99.5%) via PBT was observed for 1994 Chinook salmon of ages 2-4 years, with a 105,722-individual, 380-population baseline available for assignment. Application of a GSI-PBT system of identification to individuals in 2019 fisheries provided high-resolution estimates of stock composition, catch, and exploitation rate by CU or population, with fishery exploitation rates directly comparable to those provided by CWTs for 13 populations. GSI and PBT provide an alternate, cheaper, and more effective method in the assessment and management of Canadian-origin Chinook salmon relative to CWTs, and an opportunity for a genetics-based system to replace the current CWT system for salmon assessment.

Comparative regulomics supports pervasive selection on gene dosage following whole genome duplication.

BACKGROUND: 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. RESULTS: 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. 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. CONCLUSIONS: 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.

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