Coalescent methods reconstruct contributions of natural colonization and stocking to origins of Michigan inland Cisco (Coregonus artedi)

Fish population structure in previously glaciated regions is often influenced by natural colonization processes and human-mediated dispersal, including fish stocking. Endemic populations are of conservation interest because they may contain rare and unique genetic variation. While coregonines are native to certain Michigan inland lakes, some were stocked with fish from Great Lakes sources, calling into question the origin of extant populations. While most stocking targeted lake whitefish (Coregonus clupeaformis), cisco (C. artedi) were also stocked from the Great Lakes to inland waterbodies. We used population genetic data (microsatellite genotypes and mitochondrial (mt)DNA sequences), coalescent modeling, and approximate Bayesian computation to investigate the origins of 12 inland Michigan cisco populations. The spatial distribution of mtDNA haplotypes suggests Michigan is an introgression zone for two ancestral cisco lineages associated with separate glacial refugia. Low levels of genetic diversity and high levels of genetic divergence were observed for populations located well inland of the Great Lakes relative to populations occupying waterbodies near the Great Lakes. Estimates of recent Great Lakes gene flow ranged from 27 to 48% for populations near the Great Lakes shoreline but were substantially lower (under 8%) for populations further inland. Inland lakes with elevated recent gene flow estimates may have been recipients of stocked coregonine fry, including cisco. Low levels of genetic diversity paired with a high likelihood of endemism as indicated by strong genetic divergence and low Great Lakes population inputs suggest the analyzed cisco populations occupying southern Michigan kettle lakes are of elevated conservation interest.