Congruent population structure across paralogous and nonparalogous loci in Salish Sea chum salmon (Oncorhynchus keta).

Whole-genome duplications are major evolutionary events with a lasting impact on genome structure. Duplication events complicate genetic analyses as paralogous sequences are difficult to distinguish; consequently, paralogs are often excluded from studies. The effects of an ancient whole-genome duplication (approximately 88 MYA) are still evident in salmonids through the persistence of numerous paralogous gene sequences and partial tetrasomic inheritance. We use restriction site-associated DNA sequencing on 10 collections of chum salmon from the Salish Sea in the USA and Canada to investigate genetic diversity and population structure in both tetrasomic and rediploidized regions of the genome. We use a pedigree and high-density linkage map to identify paralogous loci and to investigate genetic variation across the genome. By applying multivariate statistical methods, we show that it is possible to characterize paralogous loci and that they display similar patterns of population structure as the diploidized portion of the genome. We find genetic associations with the adaptively important trait of run-timing in both sets of loci. By including paralogous loci in genome scans, we can observe evolutionary signals in genomic regions that have routinely been excluded from population genetic studies in other polyploid-derived species.

Major histocompatibility complex diversity is positively associated with stream water temperatures in proximate populations of sockeye salmon.

Local adaptation to heterogeneous environments generates population diversity within species, significantly increasing ecosystem stability and flows of ecosystem services. However, few studies have isolated the specific mechanisms that create and maintain this diversity. Here, we examined the relationship between water temperature in streams used for spawning and genetic diversity at a gene involved in immune function [the major histocompatibility complex (MHC)] in 14 populations of sockeye salmon (Oncorhynchus nerka) sampled across the Wood River basin in south-western Alaska. The largest influence on MHC diversity was lake basin, but we also found a significant positive correlation between average water temperature and MHC diversity. This positive relationship between temperature and MHC diversity appears to have been produced by natural selection at very local scales rather than neutral processes, as no correlation was observed between temperature and genetic diversity at 90 neutral markers. Additionally, no significant relationship was observed between temperature variability and MHC diversity. Although lake basin was the largest driver of differences in MHC diversity, our results also demonstrate that fine-scale differences in water temperature may generate variable selection regimes in populations that spawn in habitats separated by as little as 1 km. Additionally, our results indicated that some populations may be reaching a maximum level of MHC diversity. We postulate that salmon from populations in warm streams may delay spawning until late summer to avoid thermal stress as well as the elevated levels of pathogen prevalence and virulence associated with warm temperatures earlier in the summer.

[Genetic Differentiation of Sockeye Salmon Oncorhynchus nerka from Kamchatka River Basin and the Lake-River Systems of the West Coast of the Bering Sea as Inferred from Data on Single Nucleotide Polymorphism].

The variability of 45 single nucleotide polymorphism loci (SNP) was studied in sockeye salmon from the Kamchatka River basin and four lake-river systems of the west coast of the Bering Sea. Based on the genetic differentiation estimates for the largest sockeye salmon populations of Eastern Kamchatka and Chukotka, the examined samples were combined into two regional groups represented by the population of the Kamchatka River drainage, which included numerous local subpopulations and seasonal races, and the northern population grouping from the rivers of Olutorsko-Navarinsky raion, wherein the sockeye salmon from Maynypilginskaya Lake-River system was relatively isolated. Considerable divergence was observed between the island (Sarannoe Lake, Bering Island) and continental populations. Genetic heterogeneity was revealed and groups of early- and late-maturing individuals were isolated in the sample of late-run sockeye salmon from Kamchatka River. In Apuka River, subdivision of the spawning run into two genetically distinct spatial and temporal groupings was also observed. The results suggest that the differentiation of sockeye salmon samples by single nucleotide substitution frequencies was largely due to differences in the direction and strength of local selection at some loci in the population complexes and intrapopulation groupings from the examined river basins of Eastern Kamchatka, Chukotka, and Commander Islands.

Parallel signatures of selection in temporally isolated lineages of pink salmon.

Studying the effect of similar environments on diverse genetic backgrounds has long been a goal of evolutionary biologists with studies typically relying on experimental approaches. Pink salmon, a highly abundant and widely ranging salmonid, provide a naturally occurring opportunity to study the effects of similar environments on divergent genetic backgrounds due to a strict two-year semelparous life history. The species is composed of two reproductively isolated lineages with overlapping ranges that share the same spawning and rearing environments in alternate years. We used restriction-site-associated DNA (RAD) sequencing to discover and genotype approximately 8000 SNP loci in three population pairs of even- and odd-year pink salmon along a latitudinal gradient in North America. We found greater differentiation within the odd-year than within the even-year lineage and greater differentiation in the southern pair from Puget Sound than in the northern Alaskan population pairs. We identified 15 SNPs reflecting signatures of parallel selection using both a differentiation-based method (BAYESCAN) and an environmental correlation method (BAYENV). These SNPs represent genomic regions that may be particularly informative in understanding adaptive evolution in pink salmon and exploring how differing genetic backgrounds within a species respond to selection from the same natural environment.

[Intra- and interpopulation variability of southwestern Kamchatka sockeye salmon Oncorhynchus nerka inferred from the data on single nucleotide polymorphism].

The variability of 45 single nucleotide polymorphism (SNP) loci was studied in nine samples of the sockeye salmon Oncorhynchus nerka from the rivers of southwestern Kamchatka. The Wahlund effect, gametic disequilibrium at some loci, and a decrease in interpopulation genetic diversity estimates observed in samples from the Bolshaya River outlet are explained in terms of the samples' heterogeneity. Partitioning of mixed samples using some biological characteristics of the individuals led to a noticeable decrease in the frequency of these phenomena. It was demonstrated that the allelic diversity between the populations within the river Plotnikovs accounted for the larger part of genetic variation, as compared to the differentiation between the basins. The SNP loci responsible for intra- and interpopulation differentiation of sockeye salmon from the rivers of southwestern Kamchatka were identified. Some recommendations for field population genetic studies of Asian sockeye salmon were formulated.

Secondary contact and changes in coastal habitat availability influence the nonequilibrium population structure of a salmonid (Oncorhynchus keta).

Numerous empirical studies have reported lack of migration-drift equilibrium in wild populations. Determining the causes of nonequilibrium population structure is challenging because different evolutionary processes acting at a variety of spatiotemporal scales can produce similar patterns. Studies of contemporary populations in northern latitudes suggest that nonequilibrium population structure is probably caused by recent colonization of the region after the last Pleistocene ice age ended ~13,000 years ago. The chum salmon's (Oncorhynchus keta) range was fragmented by dramatic environmental changes during the Pleistocene. We investigated the population structure of chum salmon on the North Alaska Peninsula (NAP) and, using both empirical data and simulations, evaluated the effects of colonization timing and founder population heterogeneity on patterns of genetic differentiation. We screened 161 single nucleotide polymorphisms and found evidence of nonequilibrium population structure when the slope of the isolation-by-distance relationship was examined at incremental spatial scales. In addition, simulations suggested that this pattern closely matched models of recent colonization of the NAP by secondary contact. Our results agree with geological and archaeological data indicating that the NAP was a dynamic landscape that may have been more recently colonized than during the last deglaciation because of dramatic changes in coastal hydrology over the last several thousand years.

Linkage mapping with paralogs exposes regions of residual tetrasomic inheritance in chum salmon (Oncorhynchus keta).

Gene sequence similarity due to shared ancestry after a duplication event, that is paralogy, complicates the assessment of genetic variation, as sequences originating from paralogs can be difficult to distinguish. These confounded sequences are often removed prior to further analyses, leaving the underlying loci uncharacterized. Salmonids have only partially rediploidized subsequent to a whole-genome duplication; residual tetrasomic inheritance has been observed in males. We present a maximum-likelihood-based method to resolve confounded paralogous loci by observing the segregation of alleles in gynogenetic haploid offspring and demonstrate its effectiveness by constructing two linkage maps for chum salmon (Oncorhynchus keta), with and without these newly resolved loci. We find that the resolved paralogous loci are not randomly distributed across the genome. A majority are clustered in expanded subtelomeric regions of 14 linkage groups, suggesting a significant fraction of the chum salmon genome may be missed by the exclusion of paralogous loci. Transposable elements have been proposed as drivers of genome evolution and, in salmonids, may have an important role in the rediploidization process by driving differentiation between homeologous chromosomes. Consistent with that hypothesis, we find a reduced fraction of transposable element annotations among paralogous loci, and these loci predominately occur in the genomic regions that lag in the rediploidization process.

An integrated linkage map reveals candidate genes underlying adaptive variation in Chinook salmon (Oncorhynchus tshawytscha).

Salmonids are an important cultural and ecological resource exhibiting near worldwide distribution between their native and introduced range. Previous research has generated linkage maps and genomic resources for several species as well as genome assemblies for two species. We first leveraged improvements in mapping and genotyping methods to create a dense linkage map for Chinook salmon Oncorhynchus tshawytscha by assembling family data from different sources. We successfully mapped 14 620 SNP loci including 2336 paralogs in subtelomeric regions. This improved map was then used as a foundation to integrate genomic resources for gene annotation and population genomic analyses. We anchored a total of 286 scaffolds from the Atlantic salmon genome to the linkage map to provide a framework for the placement 11 728 Chinook salmon ESTs. Previously identified thermotolerance QTL were found to colocalize with several candidate genes including HSP70, a gene known to be involved in thermal response, as well as its inhibitor. Multiple regions of the genome with elevated divergence between populations were also identified, and annotation of ESTs in these regions identified candidate genes for fitness related traits such as stress response, growth and behaviour. Collectively, these results demonstrate the utility of combining genomic resources with linkage maps to enhance evolutionary inferences.

Short reads and nonmodel species: exploring the complexities of next-generation sequence assembly and SNP discovery in the absence of a reference genome.

How practical is gene and SNP discovery in a nonmodel species using short read sequences? Next-generation sequencing technologies are being applied to an increasing number of species with no reference genome. For nonmodel species, the cost, availability of existing genetic resources, genome complexity and the planned method of assembly must all be considered when selecting a sequencing platform. Our goal was to examine the feasibility and optimal methodology for SNP and gene discovery in the sockeye salmon (Oncorhynchus nerka) using short read sequences. SOLiD short reads (up to 50 bp) were generated from single- and pooled-tissue transcriptome libraries from ten sockeye salmon. The individuals were from five distinct populations from the Wood River Lakes and Mendeltna Creek, Alaska. As no reference genome was available for sockeye salmon, the SOLiD sequence reads were assembled to publicly available EST reference sequences from sockeye salmon and two closely related species, rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar). Additionally, de novo assembly of the SOLiD data was carried out, and the SOLiD reads were remapped to the de novo contigs. The results from each reference assembly were compared across all references. The number and size of contigs assembled varied with the size reference sequences. In silico SNP discovery was carried out on contigs from all four EST references; however, discovery of valid SNPs was most successful using one of the two conspecific references.

Transcriptome sequencing and high-resolution melt analysis advance single nucleotide polymorphism discovery in duplicated salmonids.

Until recently, single nucleotide polymorphism (SNP) discovery in nonmodel organisms faced many challenges, often depending upon a targeted-gene approach and Sanger sequencing of many individuals. The advent of next-generation sequencing technologies has dramatically improved discovery, but validating and testing SNPs for use in population studies remain labour intensive. Here, we detail a SNP discovery and validation pipeline that incorporates 454 pyrosequencing, high-resolution melt analysis (HRMA) and 5' nuclease genotyping. We generated 4.59×10(8) bp of redundant sequence from transcriptomes of two individual chum salmon, a highly valued species across the Pacific Rim. Nearly 26000 putative SNPs were identified--some as heterozygotes and some as homozygous for different nucleotides in the two individuals. For validation, we selected 202 templates containing single putative SNPs and conducted HRMA on 10 individuals from each of 19 populations from across the species range. Finally, 5' nuclease genotyping validated 37 SNPs that conformed to Hardy-Weinberg equilibrium expectations. Putative SNPs expressed as heterozygotes in an ascertainment individual had more than twice the validation rate of those homozygous for different alleles in the two fish, suggesting that many of the latter may have been paralogous sequence variants. Overall, this validation rate of 37/202 suggests that we have found more than 4500 templates containing SNPs for use in this population set. We anticipate using this pipeline to significantly expand the number of SNPs available for the studies of population structure and mixture analyses as well as for the studies of adaptive genetic variation in nonmodel organisms.

Single nucleotide polymorphisms across a species' range: implications for conservation studies of Pacific salmon.

Studies of the oceanic and near-shore distributions of Pacific salmon, whose migrations typically span thousands of kilometres, have become increasingly valuable in the presence of climate change, increasing hatchery production and potentially high rates of bycatch in offshore fisheries. Genetics data offer considerable insights into both the migratory routes as well as the evolutionary histories of the species. However, these types of studies require extensive data sets from spawning populations originating from across the species' range. Single nucleotide polymorphisms (SNPs) have been particularly amenable for multinational applications because they are easily shared, require little interlaboratory standardization and can be assayed through increasingly efficient technologies. Here, we discuss the development of a data set for 114 populations of chum salmon through a collaboration among North American and Asian researchers, termed PacSNP. PacSNP is focused on developing the database and applying it to problems of international interest. A data set spanning the entire range of species provides a unique opportunity to examine patterns of variability, and we review issues associated with SNP development. We found evidence of ascertainment bias within the data set, variable linkage relationships between SNPs associated with ancestral groupings and outlier loci with alleles associated with latitude.

Gene mapping of isozyme loci in chum salmon.

Recombination values were used to calculate the gene-centromere map distances for four electrophoretically detected loci, Aat3, Idh1, Idh4, and Mpi, in chum salmon (Oncorhynchus keta). We also report the results from 39 pairwise examinations for joint segregation for 10 loci in nine testcross families. Only two loci assorted nonrandomly--either Aat1 or Aat2 with Gpt. Gene-centromere distances for Aat3 and Mpi differed significantly from those reported previously for rainbow trout (Salmo gairdneri), a closely related species. This difference indicates either the presence of chromosome rearrangements or a different rate of recombination between the species. These results contrast with the conservation of linkage distances previously reported within and between other salmonid genera.

Gene-centromere mapping of 312 loci in pink salmon by half-tetrad analysis.

We estimated recombination rates between 312 loci and their centromeres in gynogenetic diploid pink salmon (Oncorhynchus gorbuscha) that we produced by initiating development with irradiated sperm and blocking the maternal second meiotic division. Amplified fragment length polymorphisms (AFLPs) were significantly more centromeric than loci identified by three other techniques (allozymes, microsatellites, and PCR using primer sequences from interspersed nuclear elements). The near absence of AFLPs in distal regions could limit their utility in constructing linkage maps. A large proportion of loci had frequency of second division segregation (y) values approaching 1.0, indicating near complete crossover interference on many chromosome arms. As predicted from models of chromosomal evolution in salmonids based upon results with allozyme loci, all duplicated microsatellite loci that shared alleles (isoloci) had y values of nearly 1.0.

Inheritance of nuclear DNA markers in gynogenetic haploid pink salmon.

We describe the inheritance of 460 PCR-based loci in the polyploid-derived pink salmon (Oncorhynchus gorbuscha) genome using gynogenetic haploid embryos. We detected a length polymorphism in a growth hormone gene (GH-2) intron that is caused by an 81 bp insertion homologous to the 3' end of the salmonid short interspersed repetitive element (SINE) SmaI. Such insertion polymorphisms within species bring into question the use of SINEs as phylogenetic markers. We confirmed that a microsatellite locus encodes a PCR-null allele that is responsible for an apparent deficit of heterozygotes in a population sample from Prince William Sound. Another set of microsatellite primers amplified alleles of the same molecular weight from both loci of a duplicated pair. In our analysis of several PCR-based multilocus techniques, we failed to detect evidence of comigrating fragments produced by duplicated loci. Segregation analysis of PCR-based markers using gynogenetic haploid embryos ensures that the interpretation of molecular variation is not complicated by heterozygosity, diploidy, or gene duplication. We urge investigators to test the inheritance of polymorphisms in salmonids prior to using them to measure genetic variation.