Temporal effective size estimates of a managed walleye Sander vitreus population and implications for genetic-based management

A population genetic window into the past and future of the walleye Sander vitreus: relation to historic walleye and the extinct "blue pike" S. v. "glaucus"

Background

Conserving genetic diversity and local adaptations are management priorities for wild populations of exploited species, which increasingly are subject to climate change, habitat loss, and pollution. These constitute growing concerns for the walleye Sander vitreus, an ecologically and economically valuable North American temperate fish with large Laurentian Great Lakes' fisheries. This study compares genetic diversity and divergence patterns across its widespread native range using mitochondrial (mt) DNA control region sequences and nine nuclear DNA microsatellite (μsat) loci, examining historic and contemporary influences. We analyze the genetic and morphological characters of a putative endemic variant– “blue pike” S. v. “glaucus” –described from Lakes Erie and Ontario, which became extinct. Walleye with turquoise-colored mucus also are evaluated, since some have questioned whether these are related to the “blue pike”.

Results

Walleye populations are distinguished by considerable genetic divergence (mean F_ST mtDNA = 0.32 ± 0.01, μsat = 0.13 ± 0.00) and substantial diversity across their range (mean heterozygosity mtDNA = 0.53 ± 0.02, μsat = 0.68 ± 0.03). Southern populations markedly differ, possessing unique haplotypes and alleles, especially the Ohio/New River population that houses the oldest haplotype and has the most pronounced divergence. Northern formerly glaciated populations have greatest diversity in Lake Erie (mean heterozygosity mtDNA = 0.79 ± 0.00, μsat = 0.72 ± 0.01). Genetic diversity was much less in the historic Lake Erie samples from 1923–1949 (mean heterozygosity mtDNA = 0.05 ± 0.01, μsat = 0.47 ± 0.06) than today. The historic “blue pike” had no unique haplotypes/alleles and there is no evidence that it comprised a separate taxon from walleye. Turquoise mucus walleye also show no genetic differentiation from other sympatric walleye and no correspondence to the “blue pike”.

Conclusions

Contemporary walleye populations possess high levels of genetic diversity and divergence, despite habitat degradation and exploitation. Genetic and previously published tagging data indicate that natal homing and spawning site philopatry led to population structure. Population patterns were shaped by climate change and drainage connections, with northern ones tracing to post-glacial recolonization. Southerly populations possess unique alleles and may provide an important genetic reservoir. Allelic frequencies of Lake Erie walleye from ~70–90 years ago significantly differed from those today, suggesting population recovery after extensive habitat loss, pollution, and exploitation. The historic “blue pike” is indistinguishable from walleye, indicating that taxonomic designation is not warranted.

Estimation of effective population size in continuously distributed populations: there goes the neighborhood

Use of genetic methods to estimate effective population size (Ne) is rapidly increasing, but all approaches make simplifying assumptions unlikely to be met in real populations. In particular, all assume a single, unstructured population, and none has been evaluated for use with continuously distributed species. We simulated continuous populations with local mating structure, as envisioned by Wright's concept of neighborhood size (NS), and evaluated performance of a single-sample estimator based on linkage disequilibrium (LD), which provides an estimate of the effective number of parents that produced the sample (Nb). Results illustrate the interacting effects of two phenomena, drift and mixture, that contribute to LD. Samples from areas equal to or smaller than a breeding window produced estimates close to the NS. As the sampling window increased in size to encompass multiple genetic neighborhoods, mixture LD from a two-locus Wahlund effect overwhelmed the reduction in drift LD from incorporating offspring from more parents. As a consequence, never approached the global Ne, even when the geographic scale of sampling was large. Results indicate that caution is needed in applying standard methods for estimating effective size to continuously distributed populations.

Census (N C) and genetically effective (N e) population size in a lake-resident population of brown trout Salmo trutta

Census (N(C)) and effective population size (N(e)) were estimated for a lake-resident population of brown trout Salmo trutta as 576 and 63, respectively. The point estimate of the ratio of effective to census population size (N(e):N(C)) for this population is 0.11 with a range of 0.06-0.26, suggesting that N(e):N(C) ratio for lake-resident populations agree more with estimates for fishes with anadromous life histories than the small ratios observed in many marine fishes.