Microsatellite variation, population structure, and bottlenecks in the threatened copperbelly water snake

Habitat loss and isolation is pervasive in the Midwest U.S. Wetlands are experiencing particularly dramatic declines, yet there is a paucity of information on the genetic impacts of these losses to obligate wetland vertebrates. We quantified the genetic variation of extant populations of a shallow wetland specialist and evaluated potential reductions in population size (i.e. bottlenecks) using seven polymorphic microsatellite DNA markers. We analyzed 228 copperbelly water snakes (Nerodia erythrogaster neglecta), representing populations from three states. Moderate genetic differentiation exists among all three regions (F (ST) = 0.12, P < 0.001), with evidence for low levels of differentiation within the federally protected Ohio region (F (ST) = 0.025, P = 0.007), and moderate to strong differentiation within the Indiana region (F (ST) = 0.23, P < 0.001). Furthermore, Bayesian clustering (i.e. STRUCTURE) supports the separation of the Indiana sites, both from each other and from all other sampling sites. However, it does not support the separation of the Ohio sites from the Kentucky sites. Differentiation among sampling sites did not appear to be related to geographic distance, but rather depended on the quality of terrestrial corridors used for dispersal. Mode shifts in allele frequencies and excess heterozygosity tests were negative, while M-ratio tests were nearly all positive, indicating the likelihood of historical rather than contemporary population bottlenecks. However, potential subspecific intergradation in the Kentucky region may have artificially lowered the M-ratio, and we suggest caution when using the M-ratio approach if intergradation is suspected. Our results have conservation implications for wetland management and management of the copperbelly populations, and emphasizes the importance of protecting wetland complexes.