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 stocking on wild fish populations has long been of interest as negative consequences such as reduced fitness and loss of genetic diversity are commonly reported outcomes. In an attempt to sustain a fishery, managers implemented nearly five decades of extensive stocking of over a million Muskellunge (Esox masquinongy), a native species in the Lower St. Lawrence River (Québec, Canada). We investigated the effect of this stocking on population genetic structure and allelic diversity in the St. Lawrence River in addition to tributaries and several stocked inland lakes. Using genotype by sequencing, we genotyped 643 individuals representing 22 locations and combined this information with forward simulations to investigate the genetic consequences of long-term stocking. Individuals native to the St. Lawrence watershed were genetically differentiated from stocking sources and tributaries, and inland lakes were naturally differentiated from the main river. Empirical data and simulations within the St. Lawrence River revealed weak stocking effects on admixture patterns. Our data suggest that the genetic structure associated with stocked fish was diluted into its relatively large effective population size. This interpretation is also consistent with a hypothesis that selection against introgression was in operation and relatively efficient within the large St. Lawrence River system. In contrast, smaller populations from adjacent tributaries and lakes displayed greater stocking-related admixture that resulted in comparatively higher heterozygosity than the St. Lawrence. Finally, individuals from inland lakes that were established by stocking maintained a close affinity with their source populations. This study illustrated a benefit of combining extensive genomic data with forward simulations for improved inference regarding population-level genetic effects of long-term stocking, and its relevance for fishery management decision making.