Conservation genetics of the eastern yellow-bellied racer (Coluber constrictor flaviventris) and bullsnake (Pituophis catenifer sayi): River valleys are critical features for snakes at northern range limits.

On the North American Great Plains, several snake species reach their northern range limit where they rely on sparsely distributed hibernacula located in major river valleys. Independent colonization histories for the river valleys and barriers to gene flow caused by the lack of suitable habitat between them may have produced genetically differentiated snake populations. To test this hypothesis, we used 10 microsatellite loci to examine the population structure of two species of conservation concern in Canada: the eastern yellow-bellied racer (Coluber constrictor flaviventris) and bullsnake (Pituophis catenifer sayi) in 3 major river valleys in southern Saskatchewan. Fixation indices (FST) showed that populations in river valleys were significantly differentiated for both species (racers, FST = 0.096, P = 0.001; bullsnakes FST = 0.045-0.157, P = 0.001). Bayesian assignment (STRUCTURE) and ordination (DAPC) strongly supported genetically differentiated groups in the geographically distinct river valleys. Finer-scale subdivision of populations within river valleys was not apparent based on our data, but is a topic that should be investigated further. Our findings highlight the importance of major river valleys for snakes at the northern extent of their ranges, and raise the possibility that populations in each river valley may warrant separate management strategies.

Fine-Scale Ecological and Genetic Population Structure of Two Whitefish (Coregoninae) Species in the Vicinity of Industrial Thermal Emissions.

Thermal pollution from industrial processes can have negative impacts on the spawning and development of cold-water fish. Point sources of thermal effluent may need to be managed to avoid affecting discrete populations. Correspondingly, we examined fine-scale ecological and genetic population structure of two whitefish species (Coregonus clupeaformis and Prosopium cylindraceum) on Lake Huron, Canada, in the immediate vicinity of thermal effluent from nuclear power generation. Niche metrics using δ13C and δ15N stable isotopes showed high levels of overlap (48.6 to 94.5%) in resource use by adult fish captured in areas affected by thermal effluent compared to nearby reference locations. Isotopic niche size, a metric of resource use diversity, was 1.3- to 2.8-fold higher than reference values in some thermally affected areas, indicative of fish mixing. Microsatellite analyses of genetic population structure (Fst, STRUCTURE and DAPC) indicated that fish captured at all locations in the vicinity of the power plant were part of a larger population extending beyond the study area. In concert, ecological and genetic markers do not support the presence of an evolutionarily significant unit in the vicinity of the power plant. Thus, future research should focus on the potential impacts of thermal emissions on development and recruitment.

Impacts of degraded DNA on restriction enzyme associated DNA sequencing (RADSeq).

Degraded DNA from suboptimal field sampling is common in molecular ecology. However, its impact on techniques that use restriction site associated next-generation DNA sequencing (RADSeq, GBS) is unknown. We experimentally examined the effects of in situDNA degradation on data generation for a modified double-digest RADSeq approach (3RAD). We generated libraries using genomic DNA serially extracted from the muscle tissue of 8 individual lake whitefish (Coregonus clupeaformis) following 0-, 12-, 48- and 96-h incubation at room temperature posteuthanasia. This treatment of the tissue resulted in input DNA that ranged in quality from nearly intact to highly sheared. All samples were sequenced as a multiplexed pool on an Illumina MiSeq. Libraries created from low to moderately degraded DNA (12-48 h) performed well. In contrast, the number of RADtags per individual, number of variable sites, and percentage of identical RADtags retained were all dramatically reduced when libraries were made using highly degraded DNA (96-h group). This reduction in performance was largely due to a significant and unexpected loss of raw reads as a result of poor quality scores. Our findings remained consistent after changes in restriction enzymes, modified fold coverage values (2- to 16-fold), and additional read-length trimming. We conclude that starting DNA quality is an important consideration for RADSeq; however, the approach remains robust until genomic DNA is extensively degraded.