Y-chromosome haplotypes are associated with variation in size and age at maturity in male Chinook salmon.

Variation in size and age at maturity is an important component of life history that is influenced by both environmental and genetic factors. In salmonids, large size confers a direct reproductive advantage through increased fecundity and egg quality in females, while larger males gain a reproductive advantage by monopolizing access to females. In addition, variation in size and age at maturity in males can be associated with different reproductive strategies; younger smaller males may gain reproductive success by sneaking among mating pairs. In both sexes, there is a trade-off between older age and increased reproductive success and increased risk of mortality by delaying reproduction. We identified four Y-chromosome haplogroups that showed regional- and population-specific variation in frequency using RADseq data for 21 populations of Alaska Chinook salmon. We then characterized the range-wide distribution of these haplogroups using GT-seq assays. These haplogroups exhibited associations with size at maturity in multiple populations, suggesting that lack of recombination between X and Y-chromosomes has allowed Y-chromosome haplogroups to capture different alleles that influence size at maturity. Ultimately, conservation of life history diversity in Chinook salmon may require conservation of Y-chromosome haplotype diversity.

Resolving allele dosage in duplicated loci using genotyping-by-sequencing data: A path forward for population genetic analysis.

Whole-genome duplications have occurred in the recent ancestors of many plants, fish and amphibians. Signals of these whole-genome duplications still exist in the form of paralogous loci. Recent advances have allowed reliable identification of paralogs in genotyping-by-sequencing (GBS) data such as that generated from restriction-site-associated DNA sequencing (RADSeq); however, excluding paralogs from analyses is still routine due to difficulties in genotyping. This exclusion of paralogs may filter a large fraction of loci, including loci that may be adaptively important or informative for population genetic analyses. We present a maximum-likelihood method for inferring allele dosage in paralogs and assess its accuracy using simulated GBS, empirical RADSeq and amplicon sequencing data from Chinook salmon. We accurately infer allele dosage for some paralogs from a RADSeq data set and show how accuracy is dependent upon both read depth and allele frequency. The amplicon sequencing data set, using RADSeq-derived markers, achieved sufficient depth to infer allele dosage for all paralogs. This study demonstrates that RADSeq locus discovery combined with amplicon sequencing of targeted loci is an effective method for incorporating paralogs into population genetic analyses.

Rank and order: evaluating the performance of SNPs for individual assignment in a non-model organism.

Single nucleotide polymorphisms (SNPs) are valuable tools for ecological and evolutionary studies. In non-model species, the use of SNPs has been limited by the number of markers available. However, new technologies and decreasing technology costs have facilitated the discovery of a constantly increasing number of SNPs. With hundreds or thousands of SNPs potentially available, there is interest in comparing and developing methods for evaluating SNPs to create panels of high-throughput assays that are customized for performance, research questions, and resources. Here we use five different methods to rank 43 new SNPs and 71 previously published SNPs for sockeye salmon: F(ST), informativeness (I(n)), average contribution to principal components (LC), and the locus-ranking programs BELS and WHICHLOCI. We then tested the performance of these different ranking methods by creating 48- and 96-SNP panels of the top-ranked loci for each method and used empirical and simulated data to obtain the probability of assigning individuals to the correct population using each panel. All 96-SNP panels performed similarly and better than the 48-SNP panels except for the 96-SNP BELS panel. Among the 48-SNP panels, panels created from F(ST), I(n), and LC ranks performed better than panels formed using the top-ranked loci from the programs BELS and WHICHLOCI. The application of ranking methods to optimize panel performance will become more important as more high-throughput assays become available.

Multiplex preamplification PCR and microsatellite validation enables accurate single nucleotide polymorphism genotyping of historical fish scales.

Incorporating historical tissues into the study of ecological, conservation and management questions can broaden the scope of population genetic research by enhancing our understanding of evolutionary processes and anthropogenic influences on natural populations. Genotyping historical and low-quality samples has been plagued by challenges associated with low amounts of template DNA and the potential for pre-existing DNA contamination among samples. We describe a two-step process designed to (i) accurately genotype large numbers of historical low-quality scale samples in a high-throughput format and (ii) screen samples for pre-existing DNA contamination. First, we describe how an efficient multiplex preamplification PCR of 45 single nucleotide polymorphisms (SNPs) can generate highly accurate genotypes with low failure and error rates in subsequent SNP genotyping reactions of individual historical scales from sockeye salmon (Oncorhynchus nerka). Second, we demonstrate how the method can be modified for the amplification of microsatellite loci to detect pre-existing DNA contamination. A total of 760 individual historical scale and 182 contemporary fin clip samples were genotyped and screened for contamination. Genotyping failure and error rates were exceedingly low and similar for both historical and contemporary samples. Pre-existing contamination in 21% of the historical samples was successfully identified by screening the amplified microsatellite loci. The advantages of automation, low failure and error rates, and ability to multiplex both the preamplification and subsequent genotyping reactions combine to make the protocol ideally suited for efficiently genotyping large numbers of potentially contaminated low-quality sources of DNA.

SNP genotyping by the 5'-nuclease reaction: advances in high-throughput genotyping with nonmodel organisms.

Population genetics studies play an increasingly important role in the management and conservation of nonmodel organisms. Unlike studies with model organisms, a typical population genetics study of a nonmodel organism may be conducted by analyzing thousands or hundreds of thousands of individuals for several dozen single nucleotide polymorphisms (SNPs). The use of robust, robotically mediated TaqMan reactions provides substantial advantages in these types of studies. We describe the methods and laboratory setup for analyzing a sustained high throughput of SNP assays in routine university or natural resource agency laboratories with a handful of thermal cyclers. Agencies sustain rates of nearly 150,000 assays per week using uniplex reactions with the Applied Biosystems 7900HT Fast Real-Time PCR System (AB 7900HT). We further describe the medium-density array run on the BioMark from Fluidigm, which increases this rate to over 500,000 assays per week by multiplexing 96 samples for 96 SNPs.