Search

The use of high‐throughput, low‐density sequencing approaches has dramatically increased in recent years in studies of eco‐evolutionary processes in wild populations and domestication in commercial aquaculture. Most of these studies focus on identifying panels of SNP loci for a single downstream application, whereas there have been few studies examining the trade‐offs for selecting panels of markers for use in multiple applications. Here, we detail the use of a bioinformatic workflow for the development of a dual‐purpose SNP panel for parentage and population assignment, which included identifying putative SNP loci, filtering for the most informative loci for the two tasks, designing effective multiplex PCR primers, optimizing the SNP panel for performance, and performing quality control steps for downstream applications. We applied this workflow to two adjacent Alaskan Sockeye Salmon populations and identified a GTseq panel of 142 SNP loci for parentage and 35 SNP loci for population assignment. Only 50–75 panel loci were necessary for >95% accurate parentage, whereas population assignment success, with all 172 panel loci, ranged from 93.9% to 96.2%. Finally, we discuss the trade‐offs and complexities of the decision‐making process that drives SNP panel development, optimization, and testing., The use of high‐throughput, low‐density sequencing approaches has dramatically increased in recent years in studies of eco‐evolutionary processes in wild populations and domestication in commercial aquaculture. Here, we detail the use of a bioinformatic workflow for the development of a dual‐purpose SNP panel for parentage and population assignment in Alaskan Sockeye Salmon. We discuss the trade‐offs and complexities of the decision‐making process that drives SNP panel development, optimization, and testing.

Targeted amplicon sequencing methods, such as genotyping-in-thousands by sequencing (GT-seq), facilitate rapid, accurate, and cost-effective analysis of hundreds of genetic loci in thousands of individuals. Development of GT-seq panels is nontrivial, but studies describing trade-offs associated with different steps of GT-seq panel development are rare. Here, we construct a dual-purpose GT-seq panel for walleye (Sander vitreus), discuss trade-offs associated with different development and genotyping approaches, and provide suggestions for researchers constructing their own GT-seq panels. Our GT-seq panel was developed using an ascertainment set consisting of restriction site-associated DNA data from 954 individuals sampled from 23 populations in Minnesota and Wisconsin, USA. We conducted simulations to test the utility of all loci for parentage analysis and genetic stock identification and designed 600 primer pairs to maximize joint accuracy for these analyses. We then performed three rounds of primer optimization to remove loci that overamplified and our final panel consisted of 436 loci. We also explored different approaches for DNA extraction, multiplexed polymerase chain reaction (PCR) amplification, and cleanup steps during the GT-seq process and discovered the following: (i) inexpensive Chelex extractions performed well for genotyping; (ii) the exonuclease I and shrimp alkaline phosphatase (ExoSAP) procedure included in some current protocols did not improve results substantially and was probably unnecessary; and (iii) it was possible to PCR amplify panels separately and combine them prior to adapter ligation. Well-optimized GT-seq panels are valuable resources for conservation genetics and our findings and suggestions should aid in their construction in myriad taxa.

Whole-genome duplication (WGD) is hypothesized to be an important evolutionary mechanism that can facilitate adaptation and speciation. Genomes that exist in states of both diploidy and residual tetraploidy are of particular interest, as mechanisms that maintain the ploidy mosaic after WGD may provide important insights into evolutionary processes. The Salmonidae family exhibits residual tetraploidy, and this, combined with the evolutionary diversity formed after an ancestral autotetraploidization event, makes this group a useful study system. In this study, we generate a novel linkage map for cisco (Coregonus artedi), an economically and culturally important fish in North America and a member of the subfamily Coregoninae, which previously lacked a high-density haploid linkage map. We also conduct comparative genomic analyses to refine our understanding of chromosomal fusion/fission history across salmonids. To facilitate this comparative approach, we use the naming strategy of protokaryotype identifiers (PKs) to associate duplicated chromosomes to their putative ancestral state. The female linkage map for cisco contains 20,292 loci, 3,225 of which are likely within residually tetraploid regions. Comparative genomic analyses revealed that patterns of residual tetrasomy are generally conserved across species, although interspecific variation persists. To determine the broad-scale retention of residual tetrasomy across the salmonids, we analyze sequence similarity of currently available genomes and find evidence of residual tetrasomy in seven of the eight chromosomes that have been previously hypothesized to show this pattern. This interspecific variation in extent of rediploidization may have important implications for understanding salmonid evolutionary histories and informing future conservation efforts.

Chinook salmon (Oncorhynchus tshawytscha) are migratory fish that are highly valued for subsistence, sport, and commercial fisheries throughout their native range. Populations of Chinook salmon in western Alaska have exhibited long-term declines, leading to restrictions on harvests. Management priorities require greater resolution for genetic stock identification (GSI) than is available with current methods. We leveraged RADseq, TaqMan, and GT-seq data originating from multiple sources, collected through time, to develop a set of GT-seq panels containing 1092 single nucleotide polymorphisms (SNPs) that improved GSI resolution in western Alaska for at-sea and in-river sampling. We generated a dense linkage map to ensure that markers selected for panels spanned the entire genome. In addition, we identified multiple RADseq markers that were associated with sex; these aligned to a 5-centimorgan (cM) region on the sex chromosome. Finally, we developed a bioinformatic pipeline to streamline analysis of GT-seq data that is capable of genotyping microhaplotypes and paralogs, both of which can improve GSI resolution over traditional single-SNP data. Our panels and pipeline provide tools for management agencies to rapidly and easily analyze large-scale genotyping projects.

Understanding how gene flow influences adaptive divergence is important for predicting adaptive responses. Theoretical studies suggest that when gene flow is high, clustering of adaptive genes in fewer genomic regions would protect adaptive alleles from recombination and thus be selected for, but few studies have tested it with empirical data. Here, we used restriction site-associated sequencing to generate genomic data for six fish species with contrasting life histories from six reaches of the Upper Mississippi River System, USA. We used four differentiation-based outlier tests and three genotype-environment association analyses to define neutral single nucleotide polymorphisms (SNPs) and outlier SNPs that were putatively under selection. We then examined the distribution of outlier SNPs along the genome and investigated whether these SNPs were found in genomic islands of differentiation and inversions. We found that gene flow varied among species, and outlier SNPs were clustered more tightly in species with higher gene flow. The two species with the highest overall F

Many salmonids have a male heterogametic (XX/XY) sex determination system, and they are supposed to have a conserved master sex determining gene (sdY), that interacts at the protein level with Foxl2 leading to the blockage of the synergistic induction of Foxl2 and Nr5a1 of the cyp19a1a promoter. However, this hypothesis of a conserved master sex determining role of sdY in salmonids is still challenged by a few exceptions, one of them being the presence of some naturally occurring “apparent” XY Chinook salmon females. Here we show that XY Chinook salmon females have a sdY gene (sdY-N183), which has one missense mutation leading to a substitution of a conserved isoleucine to an asparagine (SdY I183N). In contrast, Chinook salmon males have both a non-mutated sdY-I183 gene and the missense mutation sdY-N183 gene. The 3D model of SdY-N183 predicts that the I183N hydrophobic to hydrophilic amino acid change leads to a local modification of the β-sandwich structure of SdY. Using in vitro cell transfection assays we found that SdY-N183, like SdY-I183, is preferentially localized in the cytoplasm. However, compared to SdY-I183, SdY-N183 is more prone to degradation, its nuclear translocation by Foxl2 is reduced and SdY-N183 is unable to significantly repress the synergistic Foxl2/Nr5a1 induction of the cyp19a1a promoter. Altogether our results suggest that the sdY-N183 gene of XY Chinook females is a non-functional gene and that SdY-N183 is no longer able to promote testicular differentiation by impairing the synthesis of estrogens in the early differentiating gonads of wild Chinook salmon XY females.

Whole-genome duplications (WGDs) have occurred repeatedly and broadly throughout the evolutionary history of eukaryotes. However, the effects of WGD on genome function and evolution remain unclear. The salmonid WGD that occurred approximately 88 million years ago presents an excellent opportunity for studying the effects of WGD as ∼10–15% of each salmonid genome still exhibits tetrasomic inheritance. Herein, we utilized the rainbow trout (Oncorhynchus mykiss) genome assembly and brain transcriptome data to examine the fate of gene pairs (ohnologs) following the salmonid whole-genome duplication. We find higher sequence identity between ohnologs located within known tetrasomic regions than between ohnologs found in disomic regions, and that tetrasomically inherited ohnologs showed greater similarity in patterns of gene expression and per ohnolog were lower expressed, than disomically inherited ohnologs. Enrichment testing for Gene Ontology terms identified 49 over-represented terms in tetrasomically inherited ohnologs compared to disomic ohnologs. However, why these ohnologs are retained as tetrasomic is difficult to answer. It could be that we have identified salmonid specific “dangerous duplicates”, that is, genes that cannot take on new roles following WGD. Alternatively, there may be adaptive advantages for retaining genes as functional duplicates in tetrasomic regions, as presumably, movement of these genes into disomic regions would affect both their sequence identity and their gene expression patterns.

Migratory behaviour patterns in animals are controlled by a complex genetic architecture. Rainbow trout (Oncorhynchus mykiss) is a salmonid fish that spawns in streams but exhibits three primary life history pathways: stream-resident (fluvial), lake-migrant (adfluvial) and ocean-migrant (anadromous). Previous studies examining fluvial and anadromous O. mykiss have identified several genes associated with life history divergence including the presence of an inversion complex within chromosome 5 (Omy05) that appears to maintain a suite of linked genes controlling migratory behaviour. However, adfluvial trout are migratory without being anadromous, and the genetic basis for this life history has not been investigated from evolutionary perspectives. We sampled wild, native nonanadromous rainbow trout occupying connected stream and lake habitats in a southwest Alaskan watershed to determine whether these fish exhibit genetic divergence between fluvial and adfluvial ecotypes, and whether that divergence parallels that documented in fluvial and anadromous O. mykiss. Data from restriction site-associated DNA (RAD) sequencing revealed an association between frequencies of both the Omy05 inversion complex and other single nucleotide polymorphisms (SNPs) with habitat type (stream or lake), supporting the genetic divergence of fluvial and adfluvial individuals in sympatry. The presence of a genetic basis for migration into lakes, analogous to that documented for anadromy, indicates that the adfluvial ecotype must be recognized separately from the fluvial form of O. mykiss even though neither is anadromous. These results highlight the genetic architecture underlying migration and the importance of chromosomal inversions in promoting and sustaining intraspecific diversity.

In the original publication [...]

Understanding how gene flow influences adaptive divergence is important for predicting adaptive responses. Theoretical studies suggest that when gene flow is high, clustering of adaptive genes in fewer genomic regions would protect adaptive alleles from among-population recombination and thus be selected for, but few studies have tested this hypothesis with empirical data. Here, we used RADseq to generate genomic data for six fish species with contrasting life histories from six reaches of the Upper Mississippi River System, USA. We then conducted genome scans for genomic islands of divergence to examine the distribution of adaptive loci and investigated whether these loci were found in inversions. We found that gene flow varied among species, and adaptive loci were clustered more tightly in species with higher gene flow. For example, the two species with the highest overall FST (0.03 - 0.07) and therefore lowest gene flow showed little evidence of clusters of adaptive loci, with adaptive loci spread uniformly across the genome. In contrast, nearly all adaptive loci in the species with the lowest FST (0.0004) were found in a single large putative inversion. Two other species with intermediate gene flow (FST ~ 0.004) also showed clustered genomic architectures, with most islands of divergence clustered on a few chromosomes. These results provide important empirical evidence to support the hypothesis that increasingly clustered architectures of local adaptation are associated with high gene flow. Our study utilized a unique system with species spanning a large gradient of life histories to highlight the importance of gene flow in shaping adaptive divergence.

Dam construction and longitudinal river habitat fragmentation disrupt important life histories and movement of aquatic species. This is especially true for Oncorhynchus mykiss that exhibits both migratory (steelhead) and non-migratory (resident rainbow) forms. While the negative effects of dams on salmonids have been extensively documented, few studies have had the opportunity to compare population genetic diversity and structure prior to and following dam removal. Here we examine the impacts of the removal of two dams on the Elwha River on the population genetics of O. mykiss. Genetic data were produced from &gt, 1200 samples collected prior to dam removal from both life history forms, and post-dam removal from steelhead. We identified three genetic clusters prior to dam removal primarily explained by isolation due to dams and natural barriers. Following dam removal, genetic structure decreased and admixture increased. Despite large O. mykiss population declines after dam construction, we did not detect shifts in population genetic diversity or allele frequencies of loci putatively involved in migratory phenotypic variation. Steelhead descendants from formerly below and above dammed populations recolonized the river rapidly after dam removal, suggesting that dam construction did not significantly reduce genetic diversity underlying O. mykiss life history strategies. These results have significant evolutionary implications for the conservation of migratory adaptive potential in O. mykiss populations above current anthropogenic barriers.

Variation in age at maturity is an important contributor to life history and demographic variation within and among species. The optimal age at maturity can vary by sex, and the ability of each sex to evolve towards its fitness optimum depends on the genetic architecture of maturation. Using GWAS of RAD sequencing data, we show that age at maturity in Chinook salmon exhibits sex-specific genetic architecture, with age at maturity in males influenced by large (up to 20 Mb) male-specific haplotypes. These regions showed no such effect in females. We also provide evidence for translocation of the sex-determining gene between two different chromosomes. This has important implications for sexually antagonistic selection, particularly that sex linkage of adaptive genes may differ within and among populations based on chromosomal location of the sex-determining gene. Our findings will facilitate research into the genetic causes of shifting demography in Chinook salmon as well as a better understanding of sex determination in this species and Pacific salmon in general.

Many studies exclude loci that exhibit linkage disequilibrium (LD); however, high LD can signal reduced recombination around genomic features such as chromosome inversions or sex-determining regions. Chromosome inversions and sex-determining regions are often involved in adaptation, allowing for the inheritance of co-adapted gene complexes and for the resolution of sexually antagonistic selection through sex-specific partitioning of genetic variants. Genomic features such as these can escape detection when loci with LD are removed; in addition, failing to account for these features can introduce bias to analyses. We examined patterns of LD using network analysis to identify an overlapping chromosome inversion and sex-determining region in chum salmon. The signal of the inversion was strong enough to show up as false population substructure when the entire dataset was analyzed, while the effect of the sex-determining region on population structure was only obvious after restricting analysis to the sex chromosome. Understanding the extent and geographic distribution of inversions is now a critically important part of genetic analyses of natural populations. Our results highlight the importance of analyzing and understanding patterns of LD in genomic dataset and the perils of excluding or ignoring loci exhibiting LD. Blindly excluding loci in LD would have prevented detection of the sex-determining region and chromosome inversion while failing to understand the genomic features leading to high-LD could have resulted in false interpretations of population structure.

Targeted amplicon sequencing methods, such as genotyping-in-thousands by sequencing (GT-seq), facilitate rapid, accurate, and cost-effective analysis of hundreds of genetic loci in thousands of individuals, but studies describing detailed workflows of GTseq panel development are rare. Here, we develop a dual-purpose GT-seq panel for walleye (Sander vitreus) and discuss trade-offs associated with different development and genotyping approaches. Our GT-seq panel was developed using restriction site-associated DNA data from 954 individuals sampled from 23 populations in Minnesota and Wisconsin, USA. We then conducted simulations to test the utility of loci for parentage analysis and genetic stock identification and designed 600 primer pairs to maximize joint accuracy for these analyses. We conducted three rounds of primer optimization to remove loci that overamplified and our final panel consisted of 436 loci. Optimization focused on reducing variation in amplification rate among loci and minimizing the proportion of off-target sequence, both of which are important considerations for developing large GT-seq panels. We also explored different approaches for DNA extraction, multiplexed polymerase chain reaction (PCR) amplification, and cleanup steps during the GT-seq process and discovered the following: (1) inexpensive Chelex extractions performed well for genotyping, (2) the exonuclease I and shrimp alkaline phosphatase (ExoSAP) procedure included in some current protocols did not improve results substantially and was likely unnecessary, and (3) it was possible to PCR amplify panels separately and combine them prior to adapter ligation. Well-optimized GT-seq panels are valuable resources for conservation genetics and our findings should aid in their construction in myriad taxa.

Pink salmon, the most abundant Pacific salmon, have an obligate 2-year life cycle that leads to reproductively isolated even- and odd-year lineages. Using new and existing data, we examined the genetic structure of both lineages across their distributional range by genotyping 16 681 single-nucleotide polymorphisms for 383 individuals originating from seven pairs of even- and odd-year populations. Distinct differences in standing pools of genetic variation were identified between the lineages; we observed higher levels of heterozygosity, allelic richness, and significantly more private alleles in the odd-year lineage. However, the patterns of population structure were concordant between lineages: the Asian and northern Alaska populations displayed little differentiation but differed significantly from populations in southcentral Alaska and the Pacific Northwest. Our population structure results, in the context of known paleoecological information, suggest that both lineages occupied a northern Beringial refugium as well as a Cascadian refugium in North America during the Last Glacial Maximum. These results highlight the influence of historical patterns of habitat availability on contemporary population structure and support the hypothesis of a pre-glacial origin of the lineages.

Next generation sequencing techniques have revolutionized the collection of genome and transcriptome data from nonmodel organisms. This manuscript details the application of restriction site-associated DNA sequencing (RADseq) to generate a marker-dense genetic map for Brook Trout (Salvelinus fontinalis). The consensus map was constructed from three full-sib families totaling 176 F1 individuals. The map consisted of 42 linkage groups with a total female map size of 2502.5 cM, and a total male map size of 1863.8 cM. Synteny was confirmed with Atlantic Salmon for 38 linkage groups, with Rainbow Trout for 37 linkage groups, Arctic Char for 36 linkage groups, and with a previously published Brook Trout linkage map for 39 linkage groups. Comparative mapping confirmed the presence of 8 metacentric and 34 acrocentric chromosomes in Brook Trout. Six metacentric chromosomes seem to be conserved with Arctic Char suggesting there have been at least two species-specific fusion and fission events within the genus Salvelinus. In addition, the sex marker (sdY; sexually dimorphic on the Y chromosome) was mapped to Brook Trout BC35, which is homologous with Atlantic Salmon Ssa09qa, Rainbow Trout Omy25, and Arctic Char AC04q. Ultimately, this linkage map will be a useful resource for studies on the genome organization of Salvelinus, and facilitates comparisons of the Salvelinus genome with Salmo and Oncorhynchus.

A common challenge for fisheries management is resolving the relative contribution of closely related populations where accuracy of genetic assignment may be limited. An overlooked method for increasing assignment accuracy is the use of multi-SNP (single nucleotide polymorphism) haplotypes rather than single-SNP genotypes. Haplotypes increase power for detecting population structure, and loci derived from next-generation sequencing methods often contain multiple SNPs. We evaluated the utility of multi-SNP haplotyping for mixture analysis in western Alaska Chinook salmon (Oncorhynchus tshawytscha). Multi-SNP haplotype data increased the accuracy of mixture analysis for closely related populations by up to seven percentage points relative to single-SNP genotype data for a set of 500 loci; 90% accuracy was achievable with as few as 150 loci with multi-SNP haplotypes but required at least 300 loci with single-SNP genotypes. Individual assignment to reporting groups showed an even greater increase in accuracy of up to 17 percentage points when multi-SNP haplotypes were used. Haplotyping multiple SNPs shows promise to improve the accuracy of assigning unknown fish to population of origin whenever haplotype data are available.

Pinpointing genetic alleles that influence the timing of salmon migration might aid restoration efforts

Humans by their very nature alter the distribution of species. Be it introduction of exotic species, habitat alterations or construction of barriers, anthropogenic changes provide novel experimental systems for the molecular ecologist to study evolutionary change. These events often provide a contradiction. Effective population sizes are generally low, and introduced populations are typically characterized by reduced diversity consistent with theoretical predictions of population bottlenecks and founder effects. However, despite reduced diversity, rapid change sometimes occurs. Identification of genomic regions associated with these rapid adaptive responses to novel selection pressures provides a window into genomic regions important in adaptive diversity, both in the novel and native ranges. These studies also provide an important means to estimate the pace of evolutionary change. In this issue, Willoughby et al. () compared the heterozygosity of steelhead (the anadromous form of rainbow trout Oncorhynchus mykiss) introduced into Lake Michigan in the late 1880s to the putative source population from the ancestral California range. After 25 generations of isolation in Lake Michigan, Willoughby et al. () found consistent genomewide reductions in genetic diversity as estimated by a measure of pooled heterozygosity. Despite this overall reduction in heterozygosity, three chromosomal regions showed signals of rapid adaptation and contained genes associated with osmoregulatory and wound-healing functions.

Many studies exclude loci exhibiting linkage disequilibrium (LD); however, high LD can signal reduced recombination around genomic features such as chromosome inversions or sex-determining regions. Chromosome inversions and sex-determining regions are often involved in adaptation, allowing for the inheritance of co-adapted gene complexes and for the resolution of sexually antagonistic selection through sex-specific partitioning of genetic variants. Genomic features such as these can escape detection when loci with LD are removed; in addition, failing to account for these features can introduce bias to analyses. We examined patterns of LD using network analysis to identify an overlapping chromosome inversion and sex-determining region in chum salmon. The signal of the inversion was strong enough to show up as false population substructure when the entire dataset was analyzed, while the signal of the sex-determining region was only obvious after restricting genetic analysis to the sex chromosome. Understanding the extent and geographic distribution of inversions is now a critically important part of genetic analyses of natural populations. The results of this study highlight the importance of analyzing and understanding patterns of LD in genomic dataset and the perils of ignoring or excluding loci exhibiting LD.

Interpretation of high-throughput sequence data requires an understanding of how decisions made during bioinformatic data processing can influence results. One source of bias that is often cited is PCR clones (or PCR duplicates). PCR clones are common in restriction site-associated sequencing (RAD-seq) data sets, which are increasingly being used for molecular ecology. To determine the influence PCR clones and the bioinformatic handling of clones have on genotyping, we evaluate four RAD-seq data sets. Data sets were compared before and after clones were removed to estimate the number of clones present in RAD-seq data, quantify how often the presence of clones in a data set causes genotype calls to change compared to when clones were removed, investigate the mechanisms that lead to genotype call changes and test whether clones bias heterozygosity estimates. Our RAD-seq data sets contained 30%-60% PCR clones, but 95% of RAD-tags had five or fewer clones. Relatively few genotypes changed once clones were removed (5%-10%), and the vast majority of these changes (98%) were associated with genotypes switching from a called to no-call state or vice versa. PCR clones had a larger influence on genotype calls in individuals with low read depth but appeared to influence genotype calls at all loci similarly. Removal of PCR clones reduced the number of called genotypes by 2% but had almost no influence on estimates of heterozygosity. As such, while steps should be taken to limit PCR clones during library preparation, PCR clones are likely not a substantial source of bias for most RAD-seq studies.

Whole genome duplications have occurred in the recent ancestors of many plants, fish, and amphibians, resulting in a pervasiveness of paralogous loci and the potential for both disomic and tetrasomic inheritance in the same genome. Paralogs can be difficult to reliably genotype and are often excluded from genotyping-by-sequencing (GBS) analyses; however, removal requires paralogs to be identified which is difficult without a reference genome. We present a method for identifying paralogs in natural populations by combining two properties of duplicated loci: 1) the expected frequency of heterozygotes exceeds that for singleton loci, and 2) within heterozygotes, observed read ratios for each allele in GBS data will deviate from the 1:1 expected for singleton (diploid) loci. These deviations are often not apparent within individuals, particularly when sequence coverage is low; but, we postulated that summing allele reads for each locus over all heterozygous individuals in a population would provide sufficient power to detect deviations at those loci. We identified paralogous loci in three species: Chinook salmon (Oncorhynchus tshawytscha) which retains regions with ongoing residual tetrasomy on eight chromosome arms following a recent whole genome duplication, mountain barberry (Berberis alpina) which has a large proportion of paralogs that arose through an unknown mechanism, and dusky parrotfish (Scarus niger) which has largely re-diploidized following an ancient whole genome duplication. Importantly, this approach only requires the genotype and allele-specific read counts for each individual, information which is readily obtained from most GBS analysis pipelines. This article is protected by copyright. All rights reserved.

Loci that can be used to screen for sex in salmon can provide important information for study of both wild and cultured populations. Here, we tested for associations between sex and genotypes at thousands of loci available from a genotyping-by-sequencing (GBS) dataset to discover sex-associated loci in sockeye salmon (Oncorhynchus nerka). We discovered 7 sex-associated loci, developed high-throughput assays for 2 loci, and tested the utility of these 2 assays in 8 collections of sockeye salmon sampled throughout North America. We also screened an existing assay based on the master sex-determining gene in salmon (sdY) in these collections. The ability of GBS-derived loci to assign fish to their phenotypic sex varied substantially among collections suggesting that recombination between the loci that we discovered and the sex-determining gene has occurred. Assignment accuracy to phenotypic sex was much higher with the sdY assay but was still less than 100%. Alignment of sequences from GBS-derived loci to draft genomes for 2 salmonids provided strong evidence that many of these loci are found on chromosomes orthologous to the known sex chromosome in sockeye salmon. Our study is the first to describe the approximate location of the sex-determining region in sockeye salmon and indicates that sdY is also the master sex-determining gene in this species. However, discordances between sdY genotypes and phenotypic sex and the variable performance of GBS-derived loci warrant more research.

Understanding the genetic mechanisms that facilitate adaptive radiation is an important component of evolutionary biology. Here, we genotyped 82 neutral SNPs, seven SNPs in islands of divergence identified in a previous study (island SNPs), and a region of the major histocompatibility complex (MHC) in 32 populations of sockeye salmon to investigate whether conserved genes and genomic regions are involved in adaptive radiation. Populations representing three ecotypes were sampled from seven drainages with differing habitats and colonization histories spanning a range of 2,000 km. We found strong signatures of parallel selection across drainages at the island SNPs and MHC, suggesting that the same loci undergo divergent selection during adaptive radiation. However, patterns of differentiation at most island SNPs and the MHC were not associated with ecotypes, suggesting that these loci are responding differently to a mosaic of selective pressures. Our study provides some of the first evidence that conserved genomic islands may be involved in adaptive divergence of salmon populations. Additionally, our data provide further support for the hypothesis that sockeye salmon inhabiting rivers unconnected to lakes harbour similar genetic diversity across large distances, are likely the ancestral form of the species, and have repeatedly recolonized lake systems as they have become available after glacial recession. Finally, our results highlight the value and importance of validating outlier loci by screening additional populations and regions, a practice that will hopefully become more common in the future.

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.

In their recently corrected manuscript, "Breaking RAD: An evaluation of the utility of restriction site associated DNA sequencing for genome scans of adaptation", Lowry et al. argue that genome scans using RADseq will miss many loci under selection due to a combination of sparse marker density and low levels of linkage disequilibrium in most species. We agree that marker density and levels of LD are important considerations when designing a RADseq study; however, we dispute that RAD-based genome scans are as prone to failure as Lowry et al. suggest. Their arguments ignore the flexible nature of RADseq; the availability of different restriction enzymes and capacity for combining restriction enzymes ensures that a well-designed study should be able to generate enough markers for efficient genome coverage. We further believe that simplifying assumptions about linkage disequilibrium in their simulations are invalid in many species. Finally, it is important to note that the alternative methods proposed by Lowry et al. have limitations equal to or greater than RADseq. The wealth of studies with positive impactful findings that have used RAD genome scans instead supports the argument that properly conducted RAD genome scans are an effective method for gaining insight into ecology and evolution, particularly for non-model organisms and those with large or complex genomes.

Linkage mapping is often used to identify genes associated with phenotypic traits and for aiding genome assemblies. Still, many emerging maps do not locate centromeres - an essential component of the genomic landscape. Here, we demonstrate that for genomes with strong chiasma interference, approximate centromere placement is possible by phasing the same data used to generate linkage maps. Assuming one obligate crossover per chromosome arm, information about centromere location can be revealed by tracking the accumulated recombination frequency along linkage groups, similar to half-tetrad analyses. We validate the method on a linkage map for sockeye salmon (Oncorhynchus nerka) with known centromeric regions. Further tests suggest that the method will work well in other salmonids and other eukaryotes. However, the method performed weakly when applied to a male linkage map (rainbow trout; O. mykiss) characterized by low and unevenly distributed recombination - a general feature of male meiosis in many species. Further, a high frequency of double crossovers along chromosome arms in barley reduced resolution for locating centromeric regions on most linkage groups. Despite these limitations, our method should work well for high-density maps in species with strong recombination interference and will enrich many existing and future mapping resources.

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.

Sex-bias in gene expression is a mechanism that can generate phenotypic variance between the sexes, however, relatively little is known about how patterns of sex-bias vary during development, and how variable sex-bias is between different populations. To that end, we measured sex-bias in gene expression in the brain transcriptome of rainbow trout (Oncorhynchus mykiss) during the first two years of development. Our sampling included from the fry stage through to when O. mykiss either migrate to the ocean or remain resident and undergo sexual maturation. Samples came from two F1 lines: One from migratory steelhead trout and one from resident rainbow trout. All samples were reared in a common garden environment and RNA sequencing (RNA-seq) was used to estimate patterns of gene expression. A total of 1,716 (4.6% of total) genes showed evidence of sex-bias in gene expression in at least one time point. The majority (96.7%) of sex-biased genes were differentially expressed during the second year of development, indicating that patterns of sex-bias in expression are tied to key developmental events, such as migration and sexual maturation. Mapping of differentially expressed genes to the O. mykiss genome revealed that the X chromosome is enriched for female upregulated genes, and this may indicate a lack of dosage compensation in rainbow trout. There were many more sex-biased genes in the migratory line than the resident line suggesting differences in patterns of gene expression in the brain between populations subjected to different forces of selection. Overall, our results suggest that there is considerable variation in the extent and identity of genes exhibiting sex-bias during the first two years of life. These differentially expressed genes may be connected to developmental differences between the sexes, and/or between adopting a resident or migratory life history.

Understanding the genetic architecture of phenotypic traits can provide important information about the mechanisms and genomic regions involved in local adaptation and speciation. Here, we used genotyping-by-sequencing and a combination of previously published and newly generated data to construct sex-specific linkage maps for sockeye salmon (Oncorhynchus nerka). We then used the denser female linkage map to conduct quantitative trait locus (QTL) analysis for 4 phenotypic traits in 3 families. The female linkage map consisted of 6322 loci distributed across 29 linkage groups and was 4082 cM long, and the male map contained 2179 loci found on 28 linkage groups and was 2291 cM long. We found 26 QTL: 6 for thermotolerance, 5 for length, 9 for weight, and 6 for condition factor. QTL were distributed nonrandomly across the genome and were often found in hotspots containing multiple QTL for a variety of phenotypic traits. These hotspots may represent adaptively important regions and are excellent candidates for future research. Comparing our results with studies in other salmonids revealed several regions with overlapping QTL for the same phenotypic trait, indicating these regions may be adaptively important across multiple species. Altogether, our study demonstrates the utility of genomic data for investigating the genetic basis of important phenotypic traits. Additionally, the linkage map created here will enable future research on the genetic basis of phenotypic traits in salmon.

Regions of the genome displaying elevated differentiation (genomic islands of divergence) are thought to play an important role in local adaptation, especially in populations experiencing high gene flow. However, the characteristics of these islands as well as the functional significance of genes located within them remain largely unknown. Here, we used data from thousands of SNPs aligned to a linkage map to investigate genomic islands of divergence in three ecotypes of sockeye salmon (Oncorhynchus nerka) from a single drainage in southwestern Alaska. We found ten islands displaying high differentiation among ecotypes, contrasting neutral structure throughout the rest of the genome which was low and not partitioned by ecotype. One island on linkage group 13 was particularly large and contained six SNPs with FST > 0.14 (average FST of neutral SNPs = 0.01). Functional annotation revealed that the peak of this island contained a non-synonymous mutation in a gene involved in growth in other species (TULP4). The islands that we discovered were relatively small (80 - 402 Kb), loci found in islands did not show reduced levels of diversity, and loci in islands displayed slightly elevated linkage disequilibrium. These attributes suggest that the islands discovered here were likely generated by divergence hitchhiking, however, we cannot rule out the possibility that other mechanisms may have produced them. Our results suggest that islands of divergence serve an important role in local adaptation with gene flow and represent a significant advance towards understanding the genetic basis of ecotypic differentiation. This article is protected by copyright. All rights reserved.

Daphnia pulex is a widely used toxicological model and is known for its sensitivity to cadmium (Cd). Recent research suggests that microRNAs (miRNAs) play a critical role in animal responses to heavy metals. To investigate the functions of D. pulex miRNAs under Cd exposure, we analyzed the miRNA profiles of D. pulex after 48 h using miRNA microarrays and validated our findings by q-PCR. miRNA dpu-let-7 was identified as a stably expressed gene and used as a reference. We identified 22 and 21 differentially expressed miRNAs under low (20 μg/L CdCl2) and high-exposure (40 μg/L CdCl2) concentrations compared to controls, respectively. Cellular functions of predicted miRNA target Cd-responsive genes included oxidative stress, ion transport, mitochondrial damage, and DNA repair. An insulin-related network was also identified in relation to several Cd-responsive miRNAs. The expression of three predicted target genes for miR-71 and miR-210 were evaluated, and expression of two of them (SCN2A and SLC31A1) was negatively correlated with the expression of their regulator miRNAs. We show miR-210 is hypoxia-responsive in D. pulex and propose Cd and hypoxia induce miR-210 via a same HIF1α modulated pathway. Collectively, this research advances our understanding on the role of miRNAs in response to heavy-metal exposure.

Many migratory traits are heritable, but there is a paucity of evidence identifying the molecular mechanisms underlying differentiation in alternative migratory tactics, or in linking variation in gene expression to migratory behaviors. To that end, we examined differential gene expression in the brain transcriptome between young steelhead trout that had undergone the smoltification process, and resident rainbow trout ( Oncorhynchus mykiss ) from Sashin Creek, Alaska. Samples were sequenced from two time points: immediately before (at 20 months of age) and during (2 years of age) the presumed peak of smoltification. Smolt and resident individuals came from two genetic crosses, one where both parents were migratory, and another where both parents were residents. A total of 533 (1.9%) genes were differentially expressed between crosses, or between smolt and resident samples. These genes include some candidate migratory genes (such as POMC ), as well as genes with no previous known involvement in the migratory process. Progeny from resident parents showed more upregulated genes than progeny from migrant parents at both time points. Pathway analysis showed enrichment in 227 biological pathways between cross type, and 171 biological pathways were enriched between residents and smolts. Enriched pathways had connections to many biofunctions, and most were only enriched in one contrast. However, pathways connected to phototransduction were enriched between both cross type and migratory tactics in 11 out of 12 contrasts, suggesting there are fundamental differences in how smolts and residents process light in the brain. The genes and pathways described herein constitute an a priori candidate list for future studies of migration in other populations of O. mykiss , and other migratory species.

Little information has been gathered regarding the ontogenetic changes that contribute to differentiation between resident and migrant individuals, particularly before the onset of gross morphological and physiological changes in migratory individuals. The aim of this study was to evaluate gene expression during early development in Oncorhynchus mykiss populations with different life histories, in a tissue known to integrate environmental cues to regulate complex developmental processes and behaviours. We sampled offspring produced from migrant and resident parents, collecting whole embryos prior to the beginning of first feeding, and brain tissue at three additional time points over the first year of development. RNA sequencing for 32 individuals generated a reference transcriptome of 30 177 genes that passed count thresholds. Differential gene expression between migrant and resident offspring was observed for 1982 genes. The greatest number of differentially expressed genes occurred at 8 months of age, in the spring a full year before the obvious physiological transformation from stream-dwelling parr to sea water-adaptable smolts begins for migrant individuals. Sex and age exhibited considerable effects on differential gene expression between migrants and resident offspring. Differential gene expression was observed in genes previously associated with migration, but also in genes previously unassociated with early life history divergence. Pathway analysis revealed coordinated differential expression in genes related to phototransduction, which could modulate photoperiod responsiveness and variation in circadian rhythms. The role for early differentiation in light sensitivity and biological rhythms is particularly intriguing in understanding early brain processes involved in differentiation of migratory and resident life history types.

Fitness related traits often show spatial variation across populations of widely distributed species. Comparisons of genetic variation among populations in putatively neutral DNA markers and in phenotypic traits susceptible to selection (QST FST analysis) can be used to determine to what degree differentiation among populations can be attributed to selection or genetic drift. Traditionally, QST FST analyses require a large number of populations to achieve sufficient statistical power; however, new methods have been developed that allow QST FST comparisons to be conducted on as few as two populations if their pedigrees are informative. This study compared genetic and morphological divergence in three strains of brook trout Salvelinus fontinalis that were historically or currently used for stocking in the Lake Superior Basin. Herein we examined if morphological divergence among populations showed temporal variation, and if divergence could be attributed to selection or was indistinguishable from genetic drift. Multivariate QST FST analysis showed evidence for divergent selection between populations. Univariate analyses suggests that the pattern observed in the multivariate analyses was largely driven by divergent selection for length and weight, and moreover by divergence between the Assinica strain and each of the Iron River and Siskiwit strains rather than divergent selection between each population pair. While it could not be determined if divergence was due to natural selection or inadvertent artificial selection in hatcheries, selected differences were consistent with patterns of domestication commonly found in salmonids.

Daphnia pulex, the crustacean with the first sequenced genome, is an important organism that has been widely used in ecological and toxicological research. MicroRNAs (miRNAs) are 21–25 nucleotide small non-coding RNAs that are involved in a myriad of physiological processes. In this research, we predicted 75 D. pulex miRNAs by sequence homology and secondary structure identification from the full genome sequence. Fourteen predicted miRNAs were selected for quantitative real time polymerase chain reaction (RT-PCR) validation. Out of these, eight (mir-8, mir-9, mir-12, mir-92, mir-100, mir-133, mir-153 and mir-283) were successfully amplified and validated. Next, expression levels were quantified at three different life stages (days 4, 8 and 12 of age) using U6 spliceosomal RNA as a reference gene. The expression of mir-8, mir-9, mir-12, mir-92 and mir-100 significantly differed across time suggesting these microRNAs might play a critical role during D. pulex development. This is the first study to identify and validate miRNAs in D. pulex, which is an important first step in further studies that evaluate their roles in development and response to environmental and ecological stimuli.

We deep-sequenced the transcriptome of Chinook salmon (Oncorhynchus tshawytscha) that yielded 2.5 million high-quality reads (combined for four fish) with an average length of 378 bp. De novo assembly resulted in 44,264 contigs with an average length of 567 bp and an average depth of 29 reads. Nearly half (42 %) of the contigs were annotated through alignment against protein, gene ontology (GO) and taxonomic databases using BLASTX. Overrepresented GO categories included metabolism (32 %), biosynthesis (11 %), transport (7 %), transcription (5 %) and other important pathways (response to stress, lipid metabolism and reproduction: 3 %). We identified 3,793 putative single nucleotide polymorphisms (SNPs) in silico, of which 718 were annotated. We characterized a sample of 54 annotated SNPs within contigs with transition-to-transversion ratios