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.

Estimating contemporary effective population size in non-model species using linkage disequilibrium across thousands of loci.

Contemporary effective population size (Ne) can be estimated using linkage disequilibrium (LD) observed across pairs of loci presumed to be selectively neutral and unlinked. This method has been commonly applied to data sets containing 10-100 loci to inform conservation and study population demography. Performance of these Ne estimates could be improved by incorporating data from thousands of loci. However, these thousands of loci exist on a limited number of chromosomes, ensuring that some fraction will be physically linked. Linked loci have elevated LD due to limited recombination, which if not accounted for can cause Ne estimates to be downwardly biased. Here, we present results from coalescent and forward simulations designed to evaluate the bias of LD-based Ne estimates ([Ncirc ]e). Contrary to common perceptions, increasing the number of loci does not increase the magnitude of linkage. Although we show it is possible to identify some pairs of loci that produce unusually large r(2) values, simply removing large r(2) values is not a reliable way to eliminate bias. Fortunately, the magnitude of bias in [Ncirc ]e is strongly and negatively correlated with the process of recombination, including the number of chromosomes and their length, and this relationship provides a general way to adjust for bias. Additionally, we show that with thousands of loci, precision of [Ncirc ]e is much lower than expected based on the assumption that each pair of loci provides completely independent information.

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.

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.