Search

Artificial propagation and wild release may influence the genetic integrity of wild populations. This practice has been prevalent in fisheries for centuries and is often termed 'stocking'. In the Laurentian Great Lakes (Great Lakes here-on), walleye populations faced declines from the 1950s to the 1970s, prompting extensive stocking efforts for restoration. By the mid-2010s, walleye populations showed signs of recovery, but the genetic legacy of stocking on population structure at the genomic level remains unclear. Using a dataset of 45,600 genome-aligned SNP loci genotyped in 1075 walleye individuals, we investigated the genetic impacts of over 50 years of stocking across the Great Lakes. Population structure was associated with both natural geographic barriers and stocking from non-native sources. Admixture between Lake Erie walleye and walleye from the re-populated Tittabawassee River indicate that stocking may have re-distributed putatively adaptive alleles around the Great Lakes. Genome scans identified F ST outliers and evidence of selective sweeps, indicating local adaptation of spawning populations is likely. Notably, one genomic region showed strong differentiation between Muskegon River and walleye from the Tittabawassee River, which was re-populated by Muskegon strain walleye, suggesting admixture and selection both impact the observed genetic diversity. Overall, our study underscores how artificial propagation and translocations can significantly alter the evolutionary trajectory of populations. The findings highlight the complex interplay between stocking practices and population genetic diversity, emphasising the need for careful management strategies to preserve the genetic integrity of wild populations amidst conservation efforts., (© 2024 The Author(s). Molecular Ecology published by John Wiley & Sons Ltd.)

Blooms of Alexandrium catenella threaten to disrupt subsistence, recreational, and commercial shellfish harvest in Alaska, as the paralytic shellfish toxins (PSTs) produced pose a serious public health risk and can lead to costly shutdowns for shellfish farmers. Current methods of PST detection in the region range from monitoring programs utilizing net tows to detect A. catenella to direct shellfish tissue testing via mouse bioassay (MBA) for commercial aquaculture harvest, as well as various optional testing methods for subsistence and recreational harvesters. The efficacy and feasibility of these methods vary, and they have not been directly compared in Southeast Alaska. In this study, we sought to assess and compare A. catenella and PST early detection methods to determine which can provide the most effective and accurate warning of A. catenella blooms or PST events. We found microscope counts to be variable and prone to missing lower numbers of A. catenella, which may be indicative of bloom formation. However, quantitative polymerase chain reaction (qPCR) significantly correlated with microscope counts and was able to effectively detect even low numbers of A. catenella on all sampling days. Paralytic shellfish toxin concentrations measured by enzyme-linked immunosorbent assay and MBA significantly correlated with each other, qPCR, and some microscope counts. These results show that qPCR is an effective tool for both monitoring A. catenella and serving as a proxy for PSTs. Further work is needed to refine qPCR protocols in this system to provide bloom warnings on an actionable timescale for the aquaculture industry and other shellfish harvesters. Integr Environ Assess Manag 2024;20:2189-2202. © 2024 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC). This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA., (© 2024 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC). This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.)

Summary: We developed loco-pipe, a Snakemake pipeline that seamlessly streamlines a set of essential population genomic analyses for low-coverage whole genome sequencing (lcWGS) data. loco-pipe is highly automated, easily customizable, massively parallelized, and thus is a valuable tool for both new and experienced users of lcWGS., Availability and Implementation: loco-pipe is published under the GPLv3. It is freely available on GitHub (github.com/sudmantlab/loco-pipe) and archived on Zenodo (doi.org/10.5281/zenodo.10425920)., Competing Interests: None declared., (© The Author(s) 2024. Published by Oxford University Press.)

Background: The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems., Results: We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary., Conclusions: Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species., (© 2024. The Author(s).)

Harvest in walleye Sander vitreus fisheries is size-selective and could influence phenotypic traits of spawners; however, contributions of individual spawners to recruitment are unknown. We used parentage analyses using single nucleotide polymorphisms to test whether parental traits were related to the probability of offspring survival in Escanaba Lake, Wisconsin. From 2017 to 2020, 1339 adults and 1138 juveniles were genotyped and 66% of the offspring were assigned to at least one parent. Logistic regression indicated the probability of reproductive success (survival of age-0 to first fall) was positively (but weakly) related to total length and growth rate in females, but not age. No traits analyzed were related to reproductive success for males. Our analysis identified the model with the predictors' growth rate and year for females and the models with year and age and year for males as the most likely models to explain variation in reproductive success. Our findings indicate that interannual variation (i.e., environmental conditions) likely plays a key role in determining the probability of reproductive success in this population and provide limited support that female age, length, and growth rate influence recruitment., Competing Interests: The authors have no conflicts of interest to disclose., (© 2024 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Fish hatcheries are widely used to enhance fisheries and supplement declining wild populations. However, substantial evidence suggests that hatchery fish are subject to differential selection pressures compared to their wild counterparts. Domestication selection, or adaptation to the hatchery environment, poses a risk to wild populations if traits specific to success in the hatchery environment have a genetic component and there is subsequent introgression between hatchery and wild fish. Few studies have investigated domestication selection in hatcheries on a genomic level, and even fewer have done so in parallel across multiple hatchery-wild population pairs. In this study, we used low-coverage whole-genome sequencing to investigate signals of domestication selection in three separate hatchery populations of Chinook salmon, Oncorhynchus tshawytscha , after approximately seven generations of divergence from their corresponding wild progenitor populations. We sequenced 192 individuals from populations across Southeast Alaska and estimated genotype likelihoods at over six million loci. We discovered a total of 14 outlier peaks displaying high genetic differentiation ( F ST ) between hatchery-wild pairs, although no peaks were shared across the three comparisons. Peaks were small (53 kb on average) and often displayed elevated absolute genetic divergence ( D xy ) and linkage disequilibrium, suggesting some level of domestication selection has occurred. Our study provides evidence that domestication selection can lead to genetic differences between hatchery and wild populations in only a few generations. Additionally, our data suggest that population-specific adaptation to hatchery environments likely occurs through different genetic pathways, even for populations with similar standing genetic variation. These results highlight the need to collect paired genotype-phenotype data to understand how domestication may be affecting fitness and to identify potential management practices that may mitigate genetic risks despite multiple pathways of domestication., Competing Interests: The authors declare there are no competing interests., (© 2024 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Local adaptation is facilitated by loci clustered in relatively few regions of the genome, termed genomic islands of divergence. The mechanisms that create and maintain these islands and how they contribute to adaptive divergence is an active research topic. Here, we use sockeye salmon as a model to investigate both the mechanisms responsible for creating islands of divergence and the patterns of differentiation at these islands. Previous research suggested that multiple islands contributed to adaptive radiation of sockeye salmon. However, the low-density genomic methods used by these studies made it difficult to fully elucidate the mechanisms responsible for islands and connect genotypes to adaptive variation. We used whole genome resequencing to genotype millions of loci to investigate patterns of genetic variation at islands and the mechanisms that potentially created them. We discovered 64 islands, including 16 clustered in four genomic regions shared between two isolated populations. Characterisation of these four regions suggested that three were likely created by structural variation, while one was created by processes not involving structural variation. All four regions were small (< 600 kb), suggesting low recombination regions do not have to span megabases to be important for adaptive divergence. Differentiation at islands was not consistently associated with established population attributes. In sum, the landscape of adaptive divergence and the mechanisms that create it are complex; this complexity likely helps to facilitate fine-scale local adaptation unique to each population., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Diet analysis is a vital tool for understanding trophic interactions and is frequently used to inform conservation and management. Molecular approaches can identify diet items that are impossible to distinguish using more traditional visual-based methods. Yet, our understanding of how different variables, such as predator species or prey ration size, influence molecular diet analysis is still incomplete. Here, we conducted a large feeding trial to assess the impact that ration size, predator species, and temperature had on digestion rates estimated with visual identification, qPCR, and metabarcoding. Our trial was conducted by feeding two rations of Chinook salmon (Oncorhynchus tshawytscha) to two piscivorous fish species (largemouth bass [Micropterus salmoides] and channel catfish [Ictalurus punctatus]) held at two different temperatures (15.5 and 18.5°C) and sacrificed at regular intervals up to 120 h from the time of ingestion to quantify the prey contents remaining in the digestive tract. We found that ration size, temperature, and predator species all influenced digestion rate, with some indication that ration size had the largest influence. DNA-based analyses were able to identify salmon smolt prey in predator gut samples for much longer than visual analysis (~12 h for visual analysis vs. ~72 h for molecular analyses). Our study provides evidence that modelling the persistence of prey DNA in predator guts for molecular diet analyses may be feasible using a small set of controlling variables for many fish systems., (© 2023 John Wiley & Sons Ltd. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.)

How to identify the drivers of population connectivity remains a fundamental question in ecology and evolution. Answering this question can be challenging in aquatic environments where dynamic lake and ocean currents coupled with high levels of dispersal and gene flow can decrease the utility of modern population genetic tools. To address this challenge, we used RAD-Seq to genotype 959 yellow perch ( Perca flavescens ), a species with an ~40-day pelagic larval duration (PLD), collected from 20 sites circumscribing Lake Michigan. We also developed a novel, integrative approach that couples detailed biophysical models with eco-genetic agent-based models to generate "predictive" values of genetic differentiation. By comparing predictive and empirical values of genetic differentiation, we estimated the relative contributions for known drivers of population connectivity (e.g., currents, behavior, PLD). For the main basin populations (i.e., the largest contiguous portion of the lake), we found that high gene flow led to low overall levels of genetic differentiation among populations ( F ST  = 0.003). By far the best predictors of genetic differentiation were connectivity matrices that were derived from periods of time when there were strong and highly dispersive currents. Thus, these highly dispersive currents are driving the patterns of population connectivity in the main basin. We also found that populations from the northern and southern main basin are slightly divergent from one another, while those from Green Bay and the main basin are highly divergent ( F ST  = 0.11). By integrating biophysical and eco-genetic models with genome-wide data, we illustrate that the drivers of population connectivity can be identified in high gene flow systems., Competing Interests: The authors declare no conflicts of interest., (© 2023 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.)

Molecular methods including metabarcoding and quantitative polymerase chain reaction have shown promise for estimating species abundance by quantifying the concentration of genetic material in field samples. However, the relationship between specimen abundance and detectable concentrations of genetic material is often variable in practice. DNA mixture analysis represents an alternative approach to quantify specimen abundance based on the presence of unique alleles in a sample. The DNA mixture approach provides novel opportunities to inform ecology and conservation by estimating the absolute abundance of target taxa through molecular methods; yet, the challenges associated with genotyping many highly variable markers in mixed-DNA samples have prevented its widespread use. To advance molecular approaches for abundance estimation, we explored the utility of microhaplotypes for DNA mixture analysis by applying a 125-marker panel to 1179 Chinook salmon (Oncorhynchus tshawytscha) smolts from the Sacramento-San Joaquin Delta, California, USA. We assessed the accuracy of DNA mixture analysis through a combination of mock mixtures containing DNA from up to 20 smolts and a trophic ecological application enumerating smolts in predator diets. Mock DNA mixtures of up to 10 smolts could reliably be resolved using microhaplotypes, and increasing the panel size would likely facilitate the identification of more individuals. However, while analysis of predator gastrointestinal tract contents indicated DNA mixture analysis could discern the presence of multiple prey items, poor and variable DNA quality prevented accurate genotyping and abundance estimation. Our results indicate that DNA mixture analysis can perform well with high-quality DNA, but methodological improvements in genotyping degraded DNA are necessary before this approach can be used on marginal-quality samples., (© 2023 John Wiley & Sons Ltd.)

Inferences made from molecular data support regional stock assessment goals by providing insights into the genetic population dynamics of enigmatic species. Population genomics metrics, such as genetic diversity and population connectivity, serve as useful proxies for species health and stability. Sleeper sharks (genus Somniosus) are ecologically important deep-sea predators, estimated to reach ages of 250 to 300 yr and taking decades to reach sexual maturity. The subgenus Somniosus (Somniosus) is comprised of 3 species: S. pacificus, S. microcephalus, and S. antarcticus. Given the life history strategy of somniosids, they are vulnerable to overfishing and population declines. Further, data to assess the stocks of these species are limited. To address this deficiency, we used the reduced representation library method Restriction-site Associated DNA sequencing (RADseq) to conduct phylogenomic and population genomics analyses, providing novel information for use in stock assessments. Our results strongly support the species status of S. microcephalus (N = 79), but recover S. antarcticus (N = 2) intermixed within the S. pacificus (N = 170) clade. Population genomics analyses reveal genetic homogeneity within S. pacificus and S. microcephalus, and estimates of effective population size were in the hundreds for both species. Kinship analysis identified 2 first-degree relative pairs within our dataset (1 within each species). Our results contribute new information for stock assessments of these uniquely long-lived species by providing the strongest molecular evidence to date for the synonymization of S. antarcticus and S. pacificus, as well as estimating population genomic metrics for each supported species within the Somniosus (Somniosus) subgenus., (Published by Oxford University Press on behalf of The American Genetic Association 2022.)

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 ST (0.0303-0.0720) and therefore lowest gene flow showed little evidence of clusters of outlier SNPs, with outlier SNPs in these species spreading uniformly across the genome. In contrast, nearly all outlier SNPs in the species with the lowest F ST (0.0003) were found in a single large putative inversion. Two other species with intermediate gene flow (F ST  ~ 0.0025-0.0050) also showed clustered genomic architectures, with most islands of differentiation clustered on a few chromosomes. Our results provide important empirical evidence to support the hypothesis that increasingly clustered architecture of local adaptation is associated with high gene flow., (© 2021 John Wiley & Sons Ltd.)

Conservation and management professionals often work across jurisdictional boundaries to identify broad ecological patterns. These collaborations help to protect populations whose distributions span political borders. One common limitation to multijurisdictional collaboration is consistency in data recording and reporting. This limitation can impact genetic research, which relies on data about specific markers in an organism's genome. Incomplete overlap of markers between separate studies can prevent direct comparisons of results. Standardized marker panels can reduce the impact of this issue and provide a common starting place for new research. Genotyping-in-thousands (GTSeq) is one approach used to create standardized marker panels for nonmodel organisms. Here, we describe the development, optimization, and early assessments of a new GTSeq panel for use with walleye ( Sander vitreus ) from the Great Lakes region of North America. High genome-coverage sequencing conducted using RAD capture provided genotypes for thousands of single nucleotide polymorphisms (SNPs). From these markers, SNP and microhaplotype markers were chosen, which were informative for genetic stock identification (GSI) and kinship analysis. The final GTSeq panel contained 500 markers, including 197 microhaplotypes and 303 SNPs. Leave-one-out GSI simulations indicated that GSI accuracy should be greater than 80% in most jurisdictions. The false-positive rates of parent-offspring and full-sibling kinship identification were found to be low. Finally, genotypes could be consistently scored among separate sequencing runs >94% of the time. Results indicate that the GTSeq panel that we developed should perform well for multijurisdictional walleye research throughout the Great Lakes region., (© 2022 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

Patagonia is an understudied area, especially when it comes to population genomic studies with relevance to fishery management. However, the dynamic and heterogeneous landscape in this area can harbor an important but cryptic genetic population structure. Once such information is revealed, it can be integrated into the management of infrequently investigated species. Eleginops maclovinus is a protandrous hermaphrodite species with economic importance for local communities that are currently managed as a single genetic unit. In this study, we sampled five locations distributed across a salinity cline from Northern Patagonia to investigate the genetic population structure of E . maclovinus . We used restriction site-associated DNA (RAD) sequencing and outlier tests to obtain neutral and adaptive loci, using F ST and GEA approaches. We identified a spatial pattern of structuration with gene flow and spatial selection by environmental association. Neutral and adaptive loci showed two and three genetic groups, respectively. The effective population sizes estimated ranged from 572 (Chepu) to 14,454 (Chaitén) and were influenced more by locality than by salinity cline. We found loci putatively associated with salinity suggesting that salinity may act as a selective driver in E . maclovinus populations. These results suggest a complex interaction between genetic drift, gene flow, and natural selection in this area. Our findings also suggest several evolutionary significant units in this area, and the information should be integrated into the management of this species. We discussed the significance of these results for fishery management and suggest future directions to improve our understanding of how E . maclovinus has adapted to the dynamic waters of Northern Patagonia., Competing Interests: All authors claim no conflict of interest., (© 2022 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

Understanding patterns of genetic structure and adaptive variation in natural populations is crucial for informing conservation and management. Past genetic research using 11 microsatellite loci identified six genetic stocks of lake whitefish ( Coregonus clupeaformis ) within Lake Michigan, USA. However, ambiguity in genetic stock assignments suggested those neutral microsatellite markers did not provide adequate power for delineating lake whitefish stocks in this system, prompting calls for a genomics approach to investigate stock structure. Here, we generated a dense genomic dataset to characterize population structure and investigate patterns of neutral and adaptive genetic diversity among lake whitefish populations in Lake Michigan. Using Rapture sequencing, we genotyped 829 individuals collected from 17 baseline populations at 197,588 SNP markers after quality filtering. Although the overall pattern of genetic structure was similar to the previous microsatellite study, our genomic data provided several novel insights. Our results indicated a large genetic break between the northwestern and eastern sides of Lake Michigan, and we found a much greater level of population structure on the eastern side compared to the northwestern side. Collectively, we observed five genomic islands of adaptive divergence on five different chromosomes. Each island displayed a different pattern of population structure, suggesting that combinations of genotypes at these adaptive regions are facilitating local adaptation to spatially heterogenous selection pressures. Additionally, we identified a large linkage disequilibrium block of ~8.5 Mb on chromosome 20 that is suggestive of a putative inversion but with a low frequency of the minor haplotype. Our study provides a comprehensive assessment of population structure and adaptive variation that can help inform the management of Lake Michigan's lake whitefish fishery and highlights the utility of incorporating adaptive loci into fisheries management., Competing Interests: The authors declare no conflict of interest., (© 2022 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Differences in genomic architecture between populations, such as chromosomal inversions, may play an important role in facilitating adaptation despite opportunities for gene flow. One system where chromosomal inversions may be important for eco-evolutionary dynamics is in freshwater fishes, which often live in heterogenous environments characterized by varying levels of connectivity and varying opportunities for gene flow. In the present study, reduced representation sequencing was used to study possible adaptation in n = 345 walleye ( Sander vitreus ) from three North American waterbodies: Cedar Bluff Reservoir (Kansas, USA), Lake Manitoba (Manitoba, Canada), and Lake Winnipeg (Manitoba, Canada). Haplotype and outlier-based tests revealed a putative chromosomal inversion that contained three expressed genes and was nearly fixed in walleye assigned to Lake Winnipeg. These patterns exist despite the potential for high gene flow between these proximate Canadian lakes, suggesting that the inversion may be important for facilitating adaptive divergence between the two lakes despite gene flow. However, a specific adaptive role for the putative inversion could not be tested with the present data. Our study illuminates the importance of genomic architecture consistent with local adaptation in freshwater fishes. Furthermore, our results provide additional evidence that inversions may facilitate local adaptation in many organisms that inhabit connected but heterogenous environments., Competing Interests: The authors declare no conflicts of interest., (© 2022 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

Chronic Wasting Disease (CWD), a well-described transmissible spongiform encephalopathy of the Cervidae family, is associated with the aggregation of an abnormal isoform (PrP CWD ) of the naturally occurring host prion protein (PrP C ). Variations in the PrP gene ( PRNP ) have been associated with CWD rate of infection and disease progression. We analysed 568 free-ranging white-tailed deer ( Odocoileus virginianus ) from 9 CWD-positive Michigan counties for PRNP polymorphisms. Sampling included 185 CWD-positive, 332 CWD non-detected, and an additional 51 CWD non-detected paired to CWD-positives by sex, age, and harvest location. We found 12 polymorphic sites of which 5 were non-synonymous and resulted in a change in amino acid composition. Thirteen haplotypes were predicted, of which 11 have previously been described. Using logistic regression, consistent with other studies, we found haplotypes C (OR = 0.488, 95% CI = 0.321-0.730, P < 0.001) and F (OR = 0.122, 95% CI = 0.007-0.612, P < 0.05) and diplotype BC (OR = 0.340, 95% CI = 0.154-0.709, P < 0.01) were less likely to be found in deer infected with CWD. As has also been documented in other studies, the presence of a serine at amino acid 96 was less likely to be found in deer infected with CWD (P < 0.001, OR = 0.360 and 95% CI = 0.227-0.556). Identification of PRNP polymorphisms associated with reduced vulnerability to CWD in Michigan deer and their spatial distribution can help managers design surveillance programmesand identify and prioritize areas for CWD management.

Small, isolated populations present a challenge for conservation. The dueling effects of selection and drift in a limited pool of genetic diversity make the responses of small populations to environmental perturbations erratic and difficult to predict. This is particularly true at the edge of a species range, where populations often persist at the limits of their environmental tolerances. Populations of cisco, Coregonus artedi , in inland lakes have experienced numerous extirpations along the southern edge of their range in recent decades, which are thought to result from environmental degradation and loss of cold, well-oxygenated habitat as lakes warm. Yet, cisco extirpations do not show a clear latitudinal pattern, suggesting that local environmental factors and potentially local adaptation may influence resilience. Here, we used genomic tools to investigate the nature of this pattern of resilience. We used restriction site-associated DNA capture (Rapture) sequencing to survey genomic diversity and differentiation in southern inland lake cisco populations and compared the frequency of deleterious mutations that potentially influence fitness across lakes. We also examined haplotype diversity in a region of the major histocompatibility complex involved in stress and immune system response. We correlated these metrics to spatial and environmental factors including latitude, lake size, and measures of oxythermal habitat and found significant relationships between genetic metrics and broad and local factors. High levels of genetic differentiation among populations were punctuated by a phylogeographic break and residual patterns of isolation-by-distance. Although the prevalence of deleterious mutations and inbreeding coefficients was significantly correlated with latitude, neutral and non-neutral genetic diversity were most strongly correlated with lake surface area. Notably, differences among lakes in the availability of estimated oxythermal habitat left no clear population genomic signature. Our results shed light on the complex dynamics influencing these isolated populations and provide valuable information for their conservation., Competing Interests: The authors declare no conflict of interest., (© 2021 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Fragmentation of river systems by dams can have substantial genetic impacts on fish populations. However, genetic structure can exist naturally at small scales through processes other than isolation by physical barriers. We sampled individuals from five native fish species with varying life histories above and below a dam in the lower Boardman River, Michigan, USA, and used RADseq to investigate processes influencing genetic structure in this system. Species assessed were white sucker Catostomus commersonii , yellow perch Perca flavescens , walleye Sander vitreus , smallmouth bass Micropterus dolomieu , and rock bass Ambloplites rupestris . We detected significant differentiation within each species, but patterns of population structure varied substantially. Interestingly, genetic structure did not appear to be solely the result of fragmentation by the dam. While genetic structure in yellow perch and walleye generally coincided with "above dam" and "below dam" sampling locations, samples from our other three species did not. Specifically, samples from rock bass, smallmouth bass, and, to a much lesser extent, white sucker, aligned with a putative Great Lakes (GL) group that contained mostly individuals sampled below the dam and a putative Boardman River (BR) group that contained individuals sampled both above and below the dam, with some evidence of admixture among groups. We hypothesize that the GL and BR groups formed prior to dam construction and our samples largely represent a mixed stock that was sampled sympatrically outside of the spawning season. Support for this hypothesis is especially strong in smallmouth bass, where GL fish were 151 mm smaller than BR fish on average, suggesting a potential ontogenetic habitat shift of young GL fish into the lower river for feeding and/or refuge. Our study illuminates the complex dynamics shaping genetic structure in fragmented river systems and indicates that conclusions drawn for a single species cannot be generalized., Competing Interests: 5The authors declare no conflict of interest., (© 2021 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Stocking of fish is an important tool for maintaining fisheries but can also significantly alter population genetic structure and erode the portfolio of within-species diversity that is important for promoting resilience and adaptability. Walleye ( Sander vitreus ) are a highly valued sportfish in the midwestern United States, a region characterized by postglacial recolonization from multiple lineages and an extensive history of stocking. We leveraged genomic data and recently developed analytical approaches to explore the population structure of walleye from two midwestern states, Minnesota and Wisconsin. We genotyped 954 walleye from 23 populations at ~20,000 loci using genotyping by sequencing and tested for patterns of population structure with single-SNP and microhaplotype data. Populations from Minnesota and Wisconsin were highly differentiated from each other, with additional substructure found in each state. Population structure did not consistently adhere to drainage boundaries, as cases of high intra-drainage and low inter-drainage differentiation were observed. Low genetic structure was observed between populations from the upper Wisconsin and upper Chippewa river watersheds, which are found as few as 50 km apart and were likely homogenized through historical stocking. Nevertheless, we were able to differentiate these populations using microhaplotype-based co-ancestry analysis, providing increased resolution over previous microsatellite studies and our other single SNP-based analyses. Although our results illustrate that walleye population structure has been influenced by past stocking practices, native ancestry still exists in most populations and walleye populations may be able to purge non-native alleles and haplotypes in the absence of stocking. Our study is one of the first to use genomic tools to investigate the influence of stocking on population structure in a nonsalmonid fish and outlines a workflow leveraging recently developed analytical methods to improve resolution of complex population structure that will be highly applicable in many species and systems., Competing Interests: The authors declare no conflicts of interest., (© 2020 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

The understanding of the evolution of variable sex determination mechanisms across taxa requires comparative studies among closely related species. Following the fate of a known master sex-determining gene, we traced the evolution of sex determination in an entire teleost order (Esociformes). We discovered that the northern pike ( Esox lucius ) master sex-determining gene originated from a 65 to 90 million-year-old gene duplication event and that it remained sex linked on undifferentiated sex chromosomes for at least 56 million years in multiple species. We identified several independent species- or population-specific sex determination transitions, including a recent loss of a Y chromosome. These findings highlight the diversity of evolutionary fates of master sex-determining genes and the importance of population demographic history in sex determination studies. We hypothesize that occasional sex reversals and genetic bottlenecks provide a non-adaptive explanation for sex determination transitions., Competing Interests: QP, RF, EJ, HD, AH, BK, ER, FG, WL, LB, MT, SC, ES, GD, Fv, SC, JL, HV, KO, RG, CE, JL, CR, HH, RT, PD, KN, LJ, HP, CK, EI, VT, MS, JP, YG No competing interests declared, (© 2021, Pan et al.)

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., (© 2020 John Wiley & Sons Ltd.)

Understanding the evolutionary impacts of harvest on fish populations is important for informing fisheries management and conservation and has become a growing research topic over the last decade. However, the dynamics of fish populations are highly complex, and phenotypes can be influenced by many biotic and abiotic factors. Therefore, it is vital to collect robust data and explore multiple alternative hypotheses before concluding that fish populations are influenced by harvest. In their recently published manuscript, Bowles et al, Evolutionary Applications, 13(6):1128 conducted age/growth and genomic analysis of walleye ( Sander vitreus ) populations sampled 13-15 years (1-2.5 generations) apart and hypothesized that observed phenotypic and genomic changes in this time period were likely due to harvest. Specifically, Bowles et al. (2020) documented differential declines in size-at-age in three exploited walleye populations compared to a separate, but presumably less-exploited, reference population. Additionally, they documented population genetic differentiation in one population pair, homogenization in another, and outlier loci putatively under selection across time points. Based on their phenotypic and genetic results, they hypothesized that selective harvest had led to fisheries-induced evolution (referred to as nascent changes) in the exploited populations in as little as 1-2.5 generations. We re-analyzed their data and found that (a) sizes declined across both exploited and reference populations during the time period studied and (b) observed genomic differentiation in their study was the result of inadequate data filtering, including retaining individuals with high amounts of missing data and retaining potentially undersplit and oversplit loci that created false signals of differentiation between time points. This re-analysis did not provide evidence for phenotypic or genetic changes attributable to harvest in any of the study populations, contrasting the hypotheses presented by Bowles et al. (2020). Our comment highlights the potential pitfalls associated with conducting age/growth analyses with low sample sizes and inadequately filtering genomic datasets., (© 2020 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd. This article has been contributed to by US Government employees and their work is in the public domain in the USA.)

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., (Copyright © 2020 Blumstein et al.)

Yellow perch, Perca flavescens, is an ecologically and economically important species native to a large portion of the northern United States and southern Canada and is also a promising candidate species for aquaculture. However, no yellow perch reference genome has been available to facilitate improvements in both fisheries and aquaculture management practices. By combining Oxford Nanopore Technologies long-reads, 10X Genomics Illumina short linked reads and a chromosome contact map produced with Hi-C, we generated a high-continuity chromosome-scale yellow perch genome assembly of 877.4 Mb. It contains, in agreement with the known diploid chromosome yellow perch count, 24 chromosome-size scaffolds covering 98.8% of the complete assembly (N50 = 37.4 Mb, L50 = 11). We also provide a first characterization of the yellow perch sex determination locus that contains a male-specific duplicate of the anti-Mullerian hormone type II receptor gene (amhr2by) inserted at the proximal end of the Y chromosome (chromosome 9). Using this sex-specific information, we developed a simple PCR genotyping assay which accurately differentiates XY genetic males (amhr2by + ) from XX genetic females (amhr2by - ). Our high-quality genome assembly is an important genomic resource for future studies on yellow perch ecology, toxicology, fisheries and aquaculture research. In addition, characterization of the amhr2by gene as a candidate sex-determining gene in yellow perch provides a new example of the recurrent implication of the transforming growth factor beta pathway in fish sex determination, and highlights gene duplication as an important genomic mechanism for the emergence of new master sex determination genes., (© 2020 John Wiley & Sons Ltd.)

Effective resource management depends on our ability to partition diversity into biologically meaningful units. Recent evolutionary divergence, however, can often lead to ambiguity in morphological and genetic differentiation, complicating the delineation of valid conservation units. Such is the case with the "coregonine problem," where recent postglacial radiations of coregonines into lacustrine habitats resulted in the evolution of numerous species flocks, often with ambiguous taxonomy. The application of genomics methods is beginning to shed light on this problem and the evolutionary mechanisms underlying divergence in these ecologically and economically important fishes. Here, we used restriction site-associated DNA (RAD) sequencing to examine genetic diversity and differentiation among sympatric forms in the Coregonus artedi complex in the Apostle Islands of Lake Superior, the largest lake in the Laurentian Great Lakes. Using 29,068 SNPs, we were able to clearly distinguish among the three most common forms for the first time, as well as identify putative hybrids and potentially misidentified specimens. Population assignment rates for these forms using our RAD data were 93%-100% with the only mis-assignments arising from putative hybrids, an improvement from 62% to 77% using microsatellites. Estimates of pairwise differentiation ( F ST : 0.045-0.056) were large given the detection of hybrids, suggesting that reduced fitness of hybrid individuals may be a potential mechanism for the maintenance of differentiation. We also used a newly built C. artedi linkage map to look for islands of genetic divergence among forms and found widespread differentiation across the genome, a pattern indicative of long-term drift, suggesting that these forms have been reproductively isolated for a substantial amount of time. The results of this study provide valuable information that can be applied to develop well-informed management strategies and stress the importance of re-evaluating conservation units with genomic tools to ensure they accurately reflect species diversity., Competing Interests: None declared., (© 2020 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

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., (Published 2019. This article is a U.S. Government work and is in the public domain in the USA.)

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., (© 2019 John Wiley & Sons Ltd.)

In many species, males may exhibit alternative life histories to circumvent the costs of intrasexual competition and female courtship. While the evolution and underlying genetic and physiological mechanisms behind alternative reproductive tactics are well studied, there has been less consideration of the ecological factors that regulate their prevalence. Here, we examine six decades of age composition records from thirty-six populations of sockeye salmon (Oncorhynchus nerka) to quantify associations between spawning habitat characteristics and the prevalence of precocious sneakers known as 'jacks'. Jack prevalence was independent of neutral genetic structure among stream populations, but varied among habitat types and as a function of continuous geomorphic characteristics. Jacks were more common in streams relative to beaches and rivers, and their prevalence was negatively associated with stream width, depth, elevation, slope and area, but positively related to bank cover. Behavioural observations showed that jacks made greater use of banks, wood and shallows than guard males, indicating that their reproductive success depends on the availability of such refuges. Our results emphasize the role of the physical habitat in shaping reproductive tactic frequencies among populations, likely through local adaptation in response to variable fitness expectations under different geomorphic conditions., (© 2017 European Society For Evolutionary Biology. Journal of Evolutionary Biology © 2017 European Society For Evolutionary Biology.)

A whole-genome duplication (WGD) doubles the entire genomic content of a species and is thought to have catalysed adaptive radiation in some polyploid-origin lineages. However, little is known about general consequences of a WGD because gene duplicates (i.e., paralogs) are commonly filtered in genomic studies; such filtering may remove substantial portions of the genome in data sets from polyploid-origin species. We demonstrate a new method that enables genome-wide scans for signatures of selection at both nonduplicated and duplicated loci by taking locus-specific copy number into account. We apply this method to RAD sequence data from different ecotypes of a polyploid-origin salmonid (Oncorhynchus nerka) and reveal signatures of divergent selection that would have been missed if duplicated loci were filtered. We also find conserved signatures of elevated divergence at pairs of homeologous chromosomes with residual tetrasomic inheritance, suggesting that joint evolution of some nondiverged gene duplicates may affect the adaptive potential of these genes. These findings illustrate that including duplicated loci in genomic analyses enables novel insights into the evolutionary consequences of WGDs and local segmental gene duplications., (© 2017 John Wiley & Sons Ltd.)

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. Conversely, neutral structure observed throughout the rest of the genome was low and not partitioned by ecotype. One island on linkage group So13 was particularly large and contained six SNPs with F ST  > 0.14 (average F ST of neutral SNPs = 0.01). Functional annotation revealed that the peak of this island contained a nonsynonymous 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., (© 2016 John Wiley & Sons Ltd.)

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., (© The American Genetic Association 2016. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

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.

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., (© 2016 European Society For Evolutionary Biology. Journal of Evolutionary Biology © 2016 European Society For Evolutionary Biology.)

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., (© 2015 John Wiley & Sons Ltd.)

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., (© The American Genetic Association 2015. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Twelve eulachon (Thaleichthys pacificus, Osmeridae) populations ranging from Cook Inlet, Alaska and along the west coast of North America to the Columbia River were examined by restriction-site-associated DNA (RAD) sequencing to elucidate patterns of neutral and adaptive variation in this high geneflow species. A total of 4104 single-nucleotide polymorphisms (SNPs) were discovered across the genome, with 193 putatively adaptive SNPs as determined by F(ST) outlier tests. Estimates of population structure in eulachon with the putatively adaptive SNPs were similar, but provided greater resolution of stocks compared with a putatively neutral panel of 3911 SNPs or previous estimates with 14 microsatellites. A cline of increasing measures of genetic diversity from south to north was found in the adaptive panel, but not in the neutral markers (SNPs or microsatellites). This may indicate divergent selective pressures in differing freshwater and marine environments between regional eulachon populations and that these adaptive diversity patterns not seen with neutral markers could be a consideration when determining genetic boundaries for conservation purposes. Estimates of effective population size (N(e)) were similar with the neutral SNP panel and microsatellites and may be utilized to monitor population status for eulachon where census sizes are difficult to obtain. Greater differentiation with the panel of putatively adaptive SNPs provided higher individual assignment accuracy compared to the neutral panel or microsatellites for stock identification purposes. This study presents the first SNPs that have been developed for eulachon, and analyses with these markers highlighted the importance of integrating genome-wide neutral and adaptive genetic variation for the applications of conservation and management., (© 2015 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

The genes of the major histocompatibility complex (MHC) are an important component of the vertebrate immune system and can provide insights into the role of pathogen-mediated selection in wild populations. Here, we examined variation at the MHC class II peptide-binding region in 27 populations of sockeye salmon (Oncorhynchus nerka), distributed among three distinct spawning ecotypes, from a complex of interconnected rivers and lakes in south-western Alaska. We also obtained genotypes from 90 putatively neutral single nucleotide polymorphisms for each population to compare the relative roles of demography and selection in shaping the observed MHC variation. We found that MHC divergence was generally partitioned by spawning ecotype (lake beaches, rivers and streams) and was 30 times greater than variation at neutral markers. Additionally, we observed substantial differences in modes of selection and diversity among ecotypes, with beach populations displaying higher levels of directional selection and lower MHC diversity than the other two ecotypes. Finally, the level of MHC differentiation in our study system was comparable to that observed over much larger geographic ranges, suggesting that MHC variation does not necessarily increase with increasing spatial scale and may instead be driven by fine-scale differences in pathogen communities or pathogen virulence. The low levels of neutral structure and spatial proximity of populations in our study system indicate that MHC differentiation can be maintained through strong selective pressure even when ample opportunities for gene flow exist., (© 2014 John Wiley & Sons Ltd.)

Numerous studies of population genetics in salmonids and other anadromous fishes have revealed that population structure is generally organized into geographic hierarchies (isolation by distance), but significant structure can exist in proximate populations due to varying selective pressures (isolation by adaptation). In Chignik Lakes, Alaska, anadromous Dolly Varden char (Salvelinus malma) spawn in nearly all accessible streams throughout the watershed, including those draining directly to an estuary, Chignik Lagoon, into larger rivers, and into lakes. Collections of Dolly Varden fry from 13 streams throughout the system revealed low levels of population structure among streams emptying into freshwater. However, much stronger genetic differentiation was detected between streams emptying into freshwater and streams flowing directly into estuarine environments. This fine-scale reproductive isolation without any physical barriers to migration is likely driven by differences in selection pressures across freshwater and estuarine environments. Estuary tributaries had fewer larger, older juveniles, suggesting an alternative life history of smolting and migration to the marine environment at a much smaller size than occurs in the other populations. Therefore, genetic data were consistent with a scenario where isolation by adaptation occurs between populations of Dolly Varden in the study system, and ecological data suggest that this isolation may partially be a result of a novel Dolly Varden life history of seawater tolerance at a smaller size than previously recognized.

Recent advances in population genomics have made it possible to detect previously unidentified structure, obtain more accurate estimates of demographic parameters, and explore adaptive divergence, potentially revolutionizing the way genetic data are used to manage wild populations. Here, we identified 10 944 single-nucleotide polymorphisms using restriction-site-associated DNA (RAD) sequencing to explore population structure, demography, and adaptive divergence in five populations of Chinook salmon (Oncorhynchus tshawytscha) from western Alaska. Patterns of population structure were similar to those of past studies, but our ability to assign individuals back to their region of origin was greatly improved (>90% accuracy for all populations). We also calculated effective size with and without removing physically linked loci identified from a linkage map, a novel method for nonmodel organisms. Estimates of effective size were generally above 1000 and were biased downward when physically linked loci were not removed. Outlier tests based on genetic differentiation identified 733 loci and three genomic regions under putative selection. These markers and genomic regions are excellent candidates for future research and can be used to create high-resolution panels for genetic monitoring and population assignment. This work demonstrates the utility of genomic data to inform conservation in highly exploited species with shallow population structure.

Differences in intermale aggression have been repeatedly reported for DBA/1 and C57BL/10 mice. The results of rreciprocal crosses combined with cross-fostering procedures suggest an involvement of the Y chromosome. In the present study, the length of the Y chromosome relative to that of chromosome 19 was ascertained in five sublines of DBA/1Bg, three sublines of C57BL/10Bg, and C57BL/10.DBA/1-Y congenic stock of mice, which carries the DBA/1Bg Y chromosome. With respect to the length of the Y chromosome relative to that of chromosome of 19, two of the DBA/1 sublines are shorter than the other three DBA/1 sublines, and all DBA/1 sublines are shorter than the three C57BL/10 sublines. This is attributable primarily to the length of the Y chromosome. The C57BL/10 sublines and the BL10.D1-Y congenic stock tested exhibit the same relative lengths of the Y chromosome, suggesting that its length has changed on the C57BL/10 genetic background. There is a parallel dependence on autosomal background of the effect of the Y chromosome on intermale aggression.