Demographic patterns of walleye (Sander vitreus) reproductive success in a Wisconsin population.

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.

High-density genomic data reveal fine-scale population structure and pronounced islands of adaptive divergence in lake whitefish (Coregonus clupeaformis) from Lake Michigan.

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.

Incomplete bioinformatic filtering and inadequate age and growth analysis lead to an incorrect inference of harvested-induced changes.

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.

Absence of PCB Hot Spot Effect in Walleye Sander vitreus from Lower Green Bay of Lake Michigan.

Under certain conditions, polychlorinated biphenyl (PCB) concentration in individuals of one sex of an adult fish population may exceed that of the other sex by more than a factor of two. This phenomenon, known as the PCB hot spot effect, has been postulated to be contingent upon the following two conditions: (1) presence of a PCB hot spot in the bottom sediments of the aquatic ecosystem, such that prey PCB concentrations in the hot spot region are substantially higher than prey PCB concentrations in locations distant from the hot spot, and (2) habitat use varying between the sexes, such that individuals of one sex inhabit the hot spot region to a considerably greater degree than individuals of the other sex. To test whether PCB concentrations in walleye Sander vitreus from lower Green Bay of Lake Michigan displayed a PCB hot spot effect, whole-fish PCB concentrations were determined in ten female and ten male adult walleye from the population spawning in the Fox River, the main tributary to lower Green Bay. In addition, mark-recapture data for the Fox River walleye population were analyzed to determine differences in spatial distributions between the sexes. Results revealed that the ratio of mean PCB concentration in males to mean PCB concentration in females was only 1.13, indicating the absence of a PCB hot spot effect. This result was corroborated by the mark-recapture data analysis, which showed no significant difference in habitat use between the sexes. Thus, although condition 1 was met, condition 2 was not met. Consequently, the PCB hot spot effect was not observed in the Fox River walleye population. Lack of a significant difference in PCB congener distributions between the sexes further corroborated our conclusions.