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.

Reduced representation (RRL) sequencing approaches (e.g., RADSeq, genotyping by sequencing) require decisions about how much to invest in genome coverage and sequencing depth, as well as choices of values for adjustable bioinformatics parameters. To empirically explore the importance of these "simple" methodological decisions, we generated two independent sequencing libraries for the same 142 individual lake whitefish (Coregonus clupeaformis) using a nextRAD RRL approach: (1) a larger number of loci at low sequencing depth based on a 9mer (library A); and (2) fewer loci at higher sequencing depth based on a 10mer (library B). The fish were selected from populations with different levels of expected genetic subdivision. Each library was analyzed using the STACKS pipeline followed by three types of population structure assessment (F ST , DAPC and ADMIXTURE) with iterative increases in the stringency of sequencing depth and missing data requirements, as well as more specific a priori population maps. Library B was always able to resolve strong population differentiation in all three types of assessment regardless of the selected parameters, largely due to retention of more loci in analyses. In contrast, library A produced more variable results; increasing the minimum sequencing depth threshold (-m) resulted in a reduced number of retained loci, and therefore lost resolution at high-m values for F ST and ADMIX-TURE, but not DAPC. When detecting fine population differentiation, the population map influenced the number of loci and missing data, which generated artefacts in all downstream analyses tested. Similarly, when examining fine scale population subdivision, library B was robust to changing parameters but library A lost resolution depending on the parameter set. We used library B to examine actual subdivision in our study populations. All three types of analysis found complete subdivision among populations in Lake Huron, ON and Dore Lake, SK, Canada using 10,640 SNP loci. Weak population subdivision was detected in Lake Huron with fish from sites in the northwest , Search Bay, North Point and Hammond Bay,

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.

Round whitefish (Prosopium cylindraceum) have a broad, disjunct range across northern North America and Eurasia, and little is known about their genetic population structure. We performed genetic analyses of round whitefish from 17 sites across its range using nine microsatellites, two mitochondrial DNA (mtDNA) loci, and 4918 to 8835 single-nucleotide polymorphism (SNP) loci. Our analyses identified deep phylogenetic division between eastern and western portions of the range, likely indicative of origins from at least two separate Pleistocene glacial refugia. Regionally, microsatellites and SNPs identified congruent patterns in subdivision, and population structure was consistent with expectations based on hydrologic connectivity. Within the Laurentian Great Lakes, Lake Huron and Lake Ontario were identified as key areas of interest. Lake Huron appears to be a contemporary source population for several other Great Lakes, and Lake Ontario contains a genetically discrete group of round whitefish. In all cases, multiple genetic markers yielded similar patterns, but SNPs offered substantially enhanced resolution. We conclude that round whitefish have population subdivision on several scales important for understanding their evolutionary history and conservation planning.

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.