Current and Historical Drivers of Landscape Genetic Structure Differ in Core and Peripheral Salamander Populations

With predicted decreases in genetic diversity and greater genetic differentiation at range peripheries relative to their cores, it can be difficult to distinguish between the roles of current disturbance versus historic processes in shaping contemporary genetic patterns. To address this problem, we test for differences in historic demography and landscape genetic structure of coastal giant salamanders (Dicamptodon tenebrosus) in two core regions (Washington State, United States) versus the species' northern peripheral region (British Columbia, Canada) where the species is listed as threatened. Coalescent-based demographic simulations were consistent with a pattern of post-glacial range expansion, with both ancestral and current estimates of effective population size being much larger within the core region relative to the periphery. However, contrary to predictions of recent human-induced population decline in the less genetically diverse peripheral region, there was no genetic signature of population size change. Effects of current demographic processes on genetic structure were evident using a resistance-based landscape genetics approach. Among core populations, genetic structure was best explained by length of the growing season and isolation by resistance (i.e. a ‘flat’ landscape), but at the periphery, topography (slope and elevation) had the greatest influence on genetic structure. Although reduced genetic variation at the range periphery of D. tenebrosus appears to be largely the result of biogeographical history rather than recent impacts, our analyses suggest that inherent landscape features act to alter dispersal pathways uniquely in different parts of the species' geographic range, with implications for habitat management

Environmental and anthropogenic drivers of connectivity patterns: A basis for prioritizing conservation efforts for threatened populations

Ecosystem fragmentation and habitat loss have been the focus of landscape management due to restrictions on contemporary connectivity and dispersal of populations. Here, we used an individual approach to determine the drivers of genetic differentiation in caribou of the Canadian Rockies. We modelled the effects of isolation by distance, landscape resistance and predation risk and evaluated the consequences of individual migratory behaviour (seasonally migratory vs. sedentary) on gene flow in this threatened species. We applied distance‐based and reciprocal causal modelling approaches, testing alternative hypotheses on the effects of geographic, topographic, environmental and local population‐specific variables on genetic differentiation and relatedness among individuals. Overall, gene flow was restricted to neighbouring local populations, with spatial coordinates, local population size, groups and elevation explaining connectivity among individuals. Landscape resistance, geographic distances and predation risk were correlated with genetic distances, with correlations threefold higher for sedentary than for migratory caribou. As local caribou populations are increasingly isolated, our results indicate the need to address genetic connectivity, especially for populations with individuals displaying different migratory behaviours, whilst maintaining quality habitat both within and across the ranges of threatened populations

Data from: Rangewide landscape genetics of an endemic Pacific northwestern salamander

A species' genetic structure often varies in response to ecological and landscape processes that differ throughout the species' geographic range, yet landscape genetics studies are rarely spatially replicated. The Cope's giant salamander (Dicamptodon copei) is a neotenic, dispersal-limited amphibian with a restricted geographic range in the Pacific northwestern USA. We investigated which landscape factors affect D. copei gene flow in three regions spanning the species' range, which vary in climate, landcover and degree of anthropogenic disturbance. Least cost paths and Circuitscape resistance analyses revealed that gene flow patterns vary across the species' range, with unique combinations of landscape variables affecting gene flow in different regions. Populations in the northern coastal portions of the range had relatively high gene flow, largely facilitated by stream and river networks. Near the southeastern edge of the species' range, gene flow was more restricted overall, with relatively less facilitation by streams and more limitation by heat load index and fragmented forest cover. These results suggested that the landscape is more difficult for individuals to disperse through at the southeastern edge of the species' range, with terrestrial habitat desiccation factors becoming more limiting to gene flow. We suggest that caution be used when attempting to extrapolate landscape genetic models and conservation measures from one portion of a species' range to another

Molecular evidence for historical and recent population size reductions of tiger salamanders (Ambystoma tigrinum) in Yellowstone National Park

Abstract Population declines caused by natural and anthropogenic factors can quickly erode genetic diversity in natural populations. In this study, we examined genetic variation within 10 tiger salamander populations across northern Yellowstone National Park in Wyoming and Montana, USA using eight microsatellite loci. We tested for the genetic signature of population decline using heterozygosity excess, shifts in allele frequencies, and low ratios of allelic number to allelic size range (M-ratios). We found different results among the three tests. All 10 populations had low M-ratios, five had shifts in allele frequencies and only two had significant heterozygosity excesses. These results support theoretical expectations of different temporal signatures among bottleneck tests and suggest that both historical fish stocking, recent, sustained drought, and possibly an emerging amphibian disease have contributed to declines in effective population size

New microsatellite markers for examining genetic variation in peripheral and core populations of the Coastal Giant Salamander (Dicamptodon tenebrosus)

The Coastal Giant Salamander (Dicamptodon tenebrosus) is classified as threatened at the northern periphery of its range in British Columbia (BC), Canada, primarily due to forestry practices and habitat fragmentation. Characterising dispersal behaviour and population connectivity is therefore a priority for this region, while genetic differentiation in core versus peripheral locations remains unstudied in this wide-ranging species. We present seven new polymorphic microsatellite markers for use in population genetic analyses of D. tenebrosus. We examine locus characteristics and genetic variation in 12 streams at the species' northern range limit in BC, and within two regions representing sub-peripheral (North Cascades) and core localities (South Cascades) in Washington State, United States. In BC, the number of alleles per locus ranged from 2-5 and observed heterozygosity ranged from 0.044-0.825. Genetic differentiation was highest between BC and the South Cascades, and intermediate between BC and the North Cascades. Across loci, mean allelic richness was similar across regions, while private allelic richness was highest in the core locality (corrected for sample size). These new microsatellite loci will be a valuable addition to existing markers for detailed landscape and population genetic analyses of D. tenebrosus across its range

Appendix C. All model results for landscape correlations with G'ST.

All model results for landscape correlations with G'ST

Appendix A. Description of landscape variables used for least-cost path or regression analyses.

Description of landscape variables used for least-cost path or regression analyses

Dicamptodon copei microsatellite genotype data

Dicamptodon copei tissue samples were collected from the field in 2006-2008. DNA was extracted using DNeasy tissue kits (Qiagen) according to the manufacturer’s protocols. Eleven microsatellite loci were amplified using PCR with fluorescently labeled primers, using PCR conditions from Steele et al. (2008) (Table S1, Supporting Information). Microsatellite products were run on an ABI 3730XL automated sequencer (Applied Biosystems) at the Washington State University LBB1 core facility and genotyped using ABI GENEMAPPER 3.7 software. The data include microsatellite fragment lengths for each study site (rows) and (diploid) locus (columns)