Combining population genomics and forward simulations to investigate stocking impacts: A case study of Muskellunge (Esox masquinongy) from the St. Lawrence River basin

Understanding the genetic and evolutionary impacts of fish stocking on wild populations has long been of interest as negative consequences such as reduced fitness and loss of genetic resources are commonly reported. Nearly five decades of extensive stocking of over a million Muskellunge (Esox masquinongy) in the Lower St. Lawrence River (Quebec, Canada) was implemented by managers in an attempt to sustain a fishery. We investigated the effect of stocking on this native species genetic structure and allelic diversity in the St. Lawrence river and its tributaries as well as several stocked inland lakes. Using Genotype-By-Sequencing (GBS), we genotyped 643 individuals representing 22 sampling sites and combined this information with forward simulations to investigate the genetic consequences of stocking. Individuals native to the St. Lawrence watershed were genetically divergent from the sources used for stocking and both the St. Lawrence tributaries and inland lakes were also naturally divergent from the main stem. Empirical data and simulations revealed weak effects of stocking on admixture patterns within the St. Lawrence despite intense stocking in the past, whereas footprints of admixture were higher in the smaller stocked populations from tributaries and lakes. Altogether, our data suggests that selection against introgression has been relatively efficient within the large St. Lawrence River. In contrast, the smaller populations from adjacent tributaries and lakes still displayed stocking related admixture which apparently resulted in higher genetic diversity, suggesting that, while stocking stopped at the same time, its impact remained higher in these populations. Finally, the origin of populations from inland lakes that were established by stocking confirmed their close affinity with these source populations. This study illustrates the benefit of combining extensive genomic data with forward simulations for improved inferences regarding outcomes of population enhancement by stocking, as well as its relevance for fishery management decision making.