Habitat distribution influences dispersal and fine-scale genetic population structure of eastern foxsnakes (Mintonius gloydi) across a fragmented landscape

Genetic Variability and Admixture Zones in the Italian Populations of Turkey Oak (Quercus cerris L.)

The Turkey oak (Quercus cerris L.) is widely distributed in Italy, where it is the ecologically dominant oak on sandy and acidic soil. In this work, we analysed 23 natural populations by means of eight SSR (microsatellite) markers, to obtain the first synthetic map of genetic variability for this species and to study its dispersion during the Holocene, due to the possibility that at least one refugium during the Last Glacial Maximum was in Italy. The analyses showed a good amount of genetic variability together with fair differentiation between populations, as indicated by F_ST = 0.059. A Bayesian analysis of the amount of admixture among populations revealed the presence of four putative gene pools of origin and a rough subdivision of the populations according to their geographic location, as confirmed by the spatial analysis. No evidence for the existence of putative refugial populations was found; however, this study paves the way for the planning of conservation strategies also with regard to the relationship between Turkey oak and other oak species in Italy.

New strategies for characterizing genetic structure in wide-ranging, continuously distributed species: A Greater Sage-grouse case study

Characterizing genetic structure across a species’ range is relevant for management and conservation as it can be used to define population boundaries and quantify connectivity. Wide-ranging species residing in continuously distributed habitat pose substantial challenges for the characterization of genetic structure as many analytical methods used are less effective when isolation by distance is an underlying biological pattern. Here, we illustrate strategies for overcoming these challenges using a species of significant conservation concern, the Greater Sage-grouse (Centrocercus urophasianus), providing a new method to identify centers of genetic differentiation and combining multiple methods to help inform management and conservation strategies for this and other such species. Our objectives were to (1) describe large-scale patterns of population genetic structure and gene flow and (2) to characterize genetic subpopulation centers across the range of Greater Sage-grouse. Samples from 2,134 individuals were genotyped at 15 microsatellite loci. Using standard STRUCTURE and spatial principal components analyses, we found evidence for four or six areas of large-scale genetic differentiation and, following our novel method, 12 subpopulation centers of differentiation. Gene flow was greater, and differentiation reduced in areas of contiguous habitat (eastern Montana, most of Wyoming, much of Oregon, Nevada, and parts of Idaho). As expected, areas of fragmented habitat such as in Utah (with 6 subpopulation centers) exhibited the greatest genetic differentiation and lowest effective migration. The subpopulation centers defined here could be monitored to maintain genetic diversity and connectivity with other subpopulation centers. Many areas outside subpopulation centers are contact zones where different genetic groups converge and could be priorities for maintaining overall connectivity. Our novel method and process of leveraging multiple different analyses to find common genetic patterns provides a path forward to characterizing genetic structure in wide-ranging, continuously distributed species.

Why did the chicken not cross the road? Anthropogenic development influences the movement of a grassland bird

Movement and selection are inherently linked behaviors that form the foundation of a species' space-use patterns. Anthropogenic development in natural ecosystems can result in a variety of behavioral responses that can involve changes in either movement (speed or direction of travel) or selection (resources used), which in turn may cause population-level consequences including loss of landscape connectivity. Understanding how a species alters these different behaviors in response to human activity is essential for effective conservation. In this study, we investigated the effects of anthropogenic development such as roads, power lines and oil wells on the greater prairie-chicken (Tympanuchus cupido) movement and selection behaviors in the post-nesting and non-breeding season. Our first objective was to assess using integrated step selection analysis (iSSA) if greater prairie-chickens altered their movement behaviors or their selection patterns when encountering oil wells, power lines, or roads. Our second objective was to determine whether prairie-chickens avoided crossing linear features such as roads or power lines by comparing the number of crossing events in greater prairie-chicken movement tracks to the number of movements that crossed these features in simulated movement tracks. Based on the iSSA analysis, we found that greater prairie-chickens avoided oil wells, power lines, and roads in both seasons, and altered their rate of movement when near anthropogenic structures. However, changes in speed varied by season, with prairie-chickens increasing their movement rates in the post-nesting season when near to development and decreasing movement rates in the non-breeding season. Furthermore, prairie-chickens crossed roads and power lines at much lower rates than expected. These changes in behavior can result in habitat loss for greater prairie-chickens, as well as the potential loss of landscape connectivity. By considering both movement and selection, we were able to develop an ecological understanding of how increasing human activity may influence the space use of this species of conservation concern. Furthermore, this research provides insight into the decision-making processes by animals when they encounter anthropogenic development.

Comparing the population structure of the specialist Butler’s Gartersnake (Thamnophis butleri) and the generalist Eastern Gartersnake (Thamnophis sirtalis sirtalis) in Ontario (Canada) and Michigan (USA)

Species differing in life-history attributes vary in their responses to features within a shared landscape. We evaluated genetic structure of sympatric gartersnake species in southwestern Ontario (Canada) and southeastern Michigan (USA), where habitat fragmentation is high due to agriculture and urbanization. We surveyed genetic structure of a habitat specialist, Butler’s Gartersnake (Thamnophis butleri (Cope, 1889)), and a habitat generalist, Eastern Gartersnake (Thamnophis sirtalis sirtalis (Linnaeus, 1758)), using DNA microsatellites. Bayesian clustering, discriminant analysis of principal components, and pairwise population comparisons revealed genetic differentiation among three major regional clusters of Butler’s Gartersnake with evidence of further division within one. Genetic clustering of Butler’s Gartersnake suggest that inhospitable habitat limits dispersal. Eastern Gartersnakes showed less structure, with assignment tests implying a single genetic cluster. We found positive significant Mantel’s r for both species in the smallest distance class (<15 km), but significant isolation by distance for Butler’s Gartersnake only. These findings together imply that connectivity for Eastern Gartersnakes is less impacted by habitat loss and fragmentation or that we were less able to detect their effects. Our study shows the value of multispecies comparisons in studies seeking to understand the underlying causes of genetic structure in natural populations.

Intact landscape promotes gene flow and low genetic structuring in the threatened Eastern Massasauga Rattlesnake

Genetic structuring of wild populations is dependent on environmental, ecological, and life-history factors. The specific role environmental context plays in genetic structuring is important to conservation practitioners working with rare species across areas with varying degrees of fragmentation. We investigated fine-scale genetic patterns of the federally threatened Eastern Massasauga Rattlesnake (Sistrurus catenatus) on a relatively undisturbed island in northern Michigan, USA. This species often persists in habitat islands throughout much of its distribution due to extensive habitat loss and distance-limited dispersal. We found that the entire island population exhibited weak genetic structuring with spatially segregated variation in effective migration and genetic diversity. The low level of genetic structuring contrasts with previous studies in the southern part of the species' range at comparable fine scales (~7 km), in which much higher levels of structuring were documented. The island population's genetic structuring more closely resembles that of populations from Ontario, Canada, that occupy similarly intact habitats. Intrapopulation variation in effective migration and genetic diversity likely corresponds to the presence of large inland lakes acting as barriers and more human activity in the southern portion of the island. The observed genetic structuring in this intact landscape suggests that the Eastern Massasauga is capable of sufficient interpatch movements to reduce overall genetic structuring and colonize new habitats. Landscape mosaics with multiple habitat patches and localized barriers (e.g., large water bodies or roads) will promote gene flow and natural colonization for this declining species.

Contribution of Connectivity Assessments to Green Infrastructure (GI)

A major goal of green infrastructure (GI) is to provide functional networks of habitats and ecosystems to maintain biodiversity long-term, while at the same time optimizing landscape and ecosystem functions and services to meet human needs. Traditionally, connectivity studies are informed by movement ecology with species-specific attributes of the type and timing of movement (e.g., dispersal, foraging, mating) and movement distances, while spatial environmental data help delineate movement pathways across landscapes. To date, a range of methods and approaches are available that (a) are relevant across any organism and movement type independent of time and space scales, (b) are ready-to-use as standalone freeware or custom GIS implementation, and (c) produce appealing visual outputs that facilitate communication with land managers. However, to enhance the robustness of connectivity assessments and ensure that current trends in connectivity modeling contribute to GI with their full potential, common denominators on which to ground planning and design strategies are required. Likewise, comparable, repeatable connectivity assessments will be needed to put results of these scientific tools into practice for multi-functional GI plans and implementation. In this paper, we discuss use and limitations of state-of-the-art connectivity methods in contributing to GI implementation.

Road avoidance and its energetic consequences for reptiles

Roads are one of the most widespread human-caused habitat modifications that can increase wildlife mortality rates and alter behavior. Roads can act as barriers with variable permeability to movement and can increase distances wildlife travel to access habitats. Movement is energetically costly, and avoidance of roads could therefore impact an animal's energy budget. We tested whether reptiles avoid roads or road crossings and explored whether the energetic consequences of road avoidance decreased individual fitness. Using telemetry data from Blanding's turtles (Emydoidea blandingii; 11,658 locations of 286 turtles from 15 sites) and eastern massasaugas (Sistrurus catenatus; 1,868 locations of 49 snakes from 3 sites), we compared frequency of observed road crossings and use of road-adjacent habitat by reptiles to expected frequencies based on simulated correlated random walks. Turtles and snakes did not avoid habitats near roads, but both species avoided road crossings. Compared with simulations, turtles made fewer crossings of paved roads with low speed limits and more crossings of paved roads with high speed limits. Snakes made fewer crossings of all road types than expected based on simulated paths. Turtles traveled longer daily distances when their home range contained roads, but the predicted energetic cost was negligible: substantially less than the cost of producing one egg. Snakes with roads in their home range did not travel further per day than snakes without roads in their home range. We found that turtles and snakes avoided crossing roads, but road avoidance is unlikely to impact fitness through energetic expenditures. Therefore, mortality from vehicle strikes remains the most significant impact of roads on reptile populations.

Quantifying functional connectivity: The role of breeding habitat, abundance, and landscape features on range-wide gene flow in sage-grouse

Functional connectivity, quantified using landscape genetics, can inform conservation through the identification of factors linking genetic structure to landscape mechanisms. We used breeding habitat metrics, landscape attributes, and indices of grouse abundance, to compare fit between structural connectivity and genetic differentiation within five long-established Sage-Grouse Management Zones (MZ) I-V using microsatellite genotypes from 6,844 greater sage-grouse (Centrocercus urophasianus) collected across their 10.7 million-km2 range. We estimated structural connectivity using a circuit theory-based approach where we built resistance surfaces using thresholds dividing the landscape into "habitat" and "nonhabitat" and nodes were clusters of sage-grouse leks (where feather samples were collected using noninvasive techniques). As hypothesized, MZ-specific habitat metrics were the best predictors of differentiation. To our surprise, inclusion of grouse abundance-corrected indices did not greatly improve model fit in most MZs. Functional connectivity of breeding habitat was reduced when probability of lek occurrence dropped below 0.25 (MZs I, IV) and 0.5 (II), thresholds lower than those previously identified as required for the formation of breeding leks, which suggests that individuals are willing to travel through undesirable habitat. The individual MZ landscape results suggested terrain roughness and steepness shaped functional connectivity across all MZs. Across respective MZs, sagebrush availability (<10%-30%; II, IV, V), tree canopy cover (>10%; I, II, IV), and cultivation (>25%; I, II, IV, V) each reduced movement beyond their respective thresholds. Model validations confirmed variation in predictive ability across MZs with top resistance surfaces better predicting gene flow than geographic distance alone, especially in cases of low and high differentiation among lek groups. The resultant resistance maps we produced spatially depict the strength and redundancy of range-wide gene flow and can help direct conservation actions to maintain and restore functional connectivity for sage-grouse.

A landscape genetic analysis of swamp rabbits (Sylvilagus aquaticus) suggests forest canopy cover enhances gene flow in an agricultural matrix

Habitat loss and fragmentation pose a continued and immediate threat to wildlife and create a persistent need for ecological information at the landscape scale to guide conservation efforts. Landscape features influence population connectivity for many species and genetic analyses can be employed to determine which of these features are most important. Because population connectivity through dispersal is important to the persistence of swamp rabbits (Sylvilagus aquaticus (Bachman, 1837)) at the northern edge of their range, we used a landscape genetic approach to relate gene flow to landscape features that may impact dispersal success. We tested resistance values for attributes of land cover, watercourse corridors, canopy cover, and roads and used causal modeling and redundancy analysis to relate these representations of landscapes to genetic distance for swamp rabbits in southern Illinois, USA. Models that included canopy cover had the strongest correlations with genetic distance and were supported by our methods whereas other models were not. We concluded that high tree canopy cover enhances gene flow and landscape connectivity for swamp rabbits in southern Illinois. Our study provides important empirical evidence that landscape variables may impact the habitat connectivity of swamp rabbits. Preserving dispersal routes for swamp rabbits should focus on improving canopy cover, in both bottomland and upland, to connect suitable habitat.

Population genetic analysis reveals barriers and corridors for gene flow within and among riparian populations of a rare plant

Landscape features and life-history traits affect gene flow, migration and drift to impact on spatial genetic structure of species. Understanding this is important for managing genetic diversity of threatened species. This study assessed the spatial genetic structure of the rare riparian Grevillea sp. Cooljarloo (Proteaceae), which is restricted to a 20 km² region impacted by mining in the northern sandplains of the Southwest Australian Floristic Region, an international biodiversity hotspot. Within creek lines and floodplains, the distribution is largely continuous. Models of dispersal within riparian systems were assessed by spatial genetic analyses including population level partitioning of genetic variation and individual Bayesian clustering. High levels of genetic variation and weak isolation by distance within creek line and floodplain populations suggest large effective population sizes and strong connectivity, with little evidence for unidirectional gene flow as might be expected from hydrochory. Regional clustering of creek line populations and strong divergence among creek line populations suggest substantially lower levels of gene flow among creek lines than within creek lines. There was however a surprising amount of genetic admixture in floodplain populations, which could be explained by irregular flooding and/or movements by highly mobile nectar-feeding bird pollinators. Our results highlight that for conservation of rare riparian species, avoiding an impact to hydrodynamic processes, such as water tables and flooding dynamics, may be just as critical as avoiding direct impacts on the number of plants.

Consequences of population topology for studying gene flow using link-based landscape genetic methods

Many landscape genetic studies aim to determine the effect of landscape on gene flow between populations. These studies frequently employ link‐based methods that relate pairwise measures of historical gene flow to measures of the landscape and the geographical distance between populations. However, apart from landscape and distance, there is a third important factor that can influence historical gene flow, that is, population topology (i.e., the arrangement of populations throughout a landscape). As the population topology is determined in part by the landscape configuration, I argue that it should play a more prominent role in landscape genetics. Making use of existing literature and theoretical examples, I discuss how population topology can influence results in landscape genetic studies and how it can be taken into account to improve the accuracy of these results. In support of my arguments, I have performed a literature review of landscape genetic studies published during the first half of 2015 as well as several computer simulations of gene flow between populations. First, I argue why one should carefully consider which population pairs should be included in link‐based analyses. Second, I discuss several ways in which the population topology can be incorporated in response and explanatory variables. Third, I outline why it is important to sample populations in such a way that a good representation of the population topology is obtained. Fourth, I discuss how statistical testing for link‐based approaches could be influenced by the population topology. I conclude the article with six recommendations geared toward better incorporating population topology in link‐based landscape genetic studies.

Developing approaches for linear mixed modeling in landscape genetics through landscape-directed dispersal simulations

Dispersal can impact population dynamics and geographic variation, and thus, genetic approaches that can establish which landscape factors influence population connectivity have ecological and evolutionary importance. Mixed models that account for the error structure of pairwise datasets are increasingly used to compare models relating genetic differentiation to pairwise measures of landscape resistance. A model selection framework based on information criteria metrics or explained variance may help disentangle the ecological and landscape factors influencing genetic structure, yet there are currently no consensus for the best protocols. Here, we develop landscape-directed simulations and test a series of replicates that emulate independent empirical datasets of two species with different life history characteristics (greater sage-grouse; eastern foxsnake). We determined that in our simulated scenarios, AIC and BIC were the best model selection indices and that marginal R2 values were biased toward more complex models. The model coefficients for landscape variables generally reflected the underlying dispersal model with confidence intervals that did not overlap with zero across the entire model set. When we controlled for geographic distance, variables not in the underlying dispersal models (i.e., nontrue) typically overlapped zero. Our study helps establish methods for using linear mixed models to identify the features underlying patterns of dispersal across a variety of landscapes.

The relationship between least-cost and resistance distance

Least-cost modelling and circuit theory are common analogs used in ecology and evolution to model gene flow or animal movement across landscapes. Least-cost modelling estimates the least-cost distance, whereas circuit theory estimates resistance distance. The bias added in choosing one method over the other has not been well documented. We designed an experiment to test whether both methods were linearly related. We also tested the sensitivity of these metrics to variation in Euclidean distance, spatial autocorrelation, the number of pixels representing the landscape, and data aggregation. We found that least-cost and resistance distance were not linearly related unless a transformation was applied. Resistance distance was less sensitive to the number of pixels representing a landscape and was also less sensitive than least-cost distance to the Euclidean distance between nodes. Spatial autocorrelation did not affect either method or the relationship between methods. Resistance distance was more sensitive to aggregation in any form compared to least-cost distance. Therefore, the metric used to infer movement or gene flow and the manipulations applied to the data used to calculate these metrics may govern findings.

Differential influences of local subpopulations on regional diversity and differentiation for greater sage-grouse (Centrocercus urophasianus)

The distribution of spatial genetic variation across a region can shape evolutionary dynamics and impact population persistence. Local population dynamics and among-population dispersal rates are strong drivers of this spatial genetic variation, yet for many species we lack a clear understanding of how these population processes interact in space to shape within-species genetic variation. Here, we used extensive genetic and demographic data from 10 subpopulations of greater sage-grouse to parameterize a simulated approximate Bayesian computation (ABC) model and (i) test for regional differences in population density and dispersal rates for greater sage-grouse subpopulations in Wyoming, and (ii) quantify how these differences impact subpopulation regional influence on genetic variation. We found a close match between observed and simulated data under our parameterized model and strong variation in density and dispersal rates across Wyoming. Sensitivity analyses suggested that changes in dispersal (via landscape resistance) had a greater influence on regional differentiation, whereas changes in density had a greater influence on mean diversity across all subpopulations. Local subpopulations, however, varied in their regional influence on genetic variation. Decreases in the size and dispersal rates of central populations with low overall and net immigration (i.e. population sources) had the greatest negative impact on genetic variation. Overall, our results provide insight into the interactions among demography, dispersal and genetic variation and highlight the potential of ABC to disentangle the complexity of regional population dynamics and project the genetic impact of changing conditions.

Anthropogenic Habitats Facilitate Dispersal of an Early Successional Obligate: Implications for Restoration of an Endangered Ecosystem

Landscape modification and habitat fragmentation disrupt the connectivity of natural landscapes, with major consequences for biodiversity. Species that require patchily distributed habitats, such as those that specialize on early successional ecosystems, must disperse through a landscape matrix with unsuitable habitat types. We evaluated landscape effects on dispersal of an early successional obligate, the New England cottontail (Sylvilagus transitionalis). Using a landscape genetics approach, we identified barriers and facilitators of gene flow and connectivity corridors for a population of cottontails in the northeastern United States. We modeled dispersal in relation to landscape structure and composition and tested hypotheses about the influence of habitat fragmentation on gene flow. Anthropogenic and natural shrubland habitats facilitated gene flow, while the remainder of the matrix, particularly development and forest, impeded gene flow. The relative influence of matrix habitats differed between study areas in relation to a fragmentation gradient. Barrier features had higher explanatory power in the more fragmented site, while facilitating features were important in the less fragmented site. Landscape models that included a simultaneous barrier and facilitating effect of roads had higher explanatory power than models that considered either effect separately, supporting the hypothesis that roads act as both barriers and facilitators at all spatial scales. The inclusion of LiDAR-identified shrubland habitat improved the fit of our facilitator models. Corridor analyses using circuit and least cost path approaches revealed the importance of anthropogenic, linear features for restoring connectivity between the study areas. In fragmented landscapes, human-modified habitats may enhance functional connectivity by providing suitable dispersal conduits for early successional specialists.

Assessing Genetic Structure in Common but Ecologically Distinct Carnivores: The Stone Marten and Red Fox

The identification of populations and spatial genetic patterns is important for ecological and conservation research, and spatially explicit individual-based methods have been recognised as powerful tools in this context. Mammalian carnivores are intrinsically vulnerable to habitat fragmentation but not much is known about the genetic consequences of fragmentation in common species. Stone martens (Martes foina) and red foxes (Vulpes vulpes) share a widespread Palearctic distribution and are considered habitat generalists, but in the Iberian Peninsula stone martens tend to occur in higher quality habitats. We compared their genetic structure in Portugal to see if they are consistent with their differences in ecological plasticity, and also to illustrate an approach to explicitly delineate the spatial boundaries of consistently identified genetic units. We analysed microsatellite data using spatial Bayesian clustering methods (implemented in the software BAPS, GENELAND and TESS), a progressive partitioning approach and a multivariate technique (Spatial Principal Components Analysis-sPCA). Three consensus Bayesian clusters were identified for the stone marten. No consensus was achieved for the red fox, but one cluster was the most probable clustering solution. Progressive partitioning and sPCA suggested additional clusters in the stone marten but they were not consistent among methods and were geographically incoherent. The contrasting results between the two species are consistent with the literature reporting stricter ecological requirements of the stone marten in the Iberian Peninsula. The observed genetic structure in the stone marten may have been influenced by landscape features, particularly rivers, and fragmentation. We suggest that an approach based on a consensus clustering solution of multiple different algorithms may provide an objective and effective means to delineate potential boundaries of inferred subpopulations. sPCA and progressive partitioning offer further verification of possible population structure and may be useful for revealing cryptic spatial genetic patterns worth further investigation.

Polygamy and an absence of fine-scale structure in Dendroctonus ponderosae (Hopk.) (Coleoptera: Curcilionidae) confirmed using molecular markers

An understanding of mating systems and fine-scale spatial genetic structure is required to effectively manage forest pest species such as Dendroctonus ponderosae (mountain pine beetle). Here we used genome-wide single-nucleotide polymorphisms to assess the fine-scale genetic structure and mating system of D. ponderosae collected from a single stand in Alberta, Canada. Fine-scale spatial genetic structure was absent within the stand and the majority of genetic variation was best explained at the individual level. Relatedness estimates support previous reports of pre-emergence mating. Parentage assignment tests indicate that a polygamous mating system better explains the relationships among individuals within a gallery than the previously reported female monogamous/male polygynous system. Furthermore, there is some evidence to suggest that females may exploit the galleries of other females, at least under epidemic conditions. Our results suggest that current management models are likely to be effective across large geographic areas based on the absence of fine-scale genetic structure.

Landscape characteristics influencing the genetic structure of greater sage-grouse within the stronghold of their range: a holistic modeling approach

Given the significance of animal dispersal to population dynamics and geographic variability, understanding how dispersal is impacted by landscape patterns has major ecological and conservation importance. Speaking to the importance of dispersal, the use of linear mixed models to compare genetic differentiation with pairwise resistance derived from landscape resistance surfaces has presented new opportunities to disentangle the menagerie of factors behind effective dispersal across a given landscape. Here, we combine these approaches with novel resistance surface parameterization to determine how the distribution of high- and low-quality seasonal habitat and individual landscape components shape patterns of gene flow for the greater sage-grouse (Centrocercus urophasianus) across Wyoming. We found that pairwise resistance derived from the distribution of low-quality nesting and winter, but not summer, seasonal habitat had the strongest correlation with genetic differentiation. Although the patterns were not as strong as with habitat distribution, multivariate models with sagebrush cover and landscape ruggedness or forest cover and ruggedness similarly had a much stronger fit with genetic differentiation than an undifferentiated landscape. In most cases, landscape resistance surfaces transformed with 17.33-km-diameter moving windows were preferred, suggesting small-scale differences in habitat were unimportant at this large spatial extent. Despite the emergence of these overall patterns, there were differences in the selection of top models depending on the model selection criteria, suggesting research into the most appropriate criteria for landscape genetics is required. Overall, our results highlight the importance of differences in seasonal habitat preferences to patterns of gene flow and suggest the combination of habitat suitability modeling and linear mixed models with our resistance parameterization is a powerful approach to discerning the effects of landscape on gene flow.

Assessing the permeability of landscape features to animal movement: using genetic structure to infer functional connectivity

Human-altered environments often challenge native species with a complex spatial distribution of resources. Hostile landscape features can inhibit animal movement (i.e., genetic exchange), while other landscape attributes facilitate gene flow. The genetic attributes of organisms inhabiting such complex environments can reveal the legacy of their movements through the landscape. Thus, by evaluating landscape attributes within the context of genetic connectivity of organisms within the landscape, we can elucidate how a species has coped with the enhanced complexity of human altered environments. In this research, we utilized genetic data from eastern chipmunks (Tamias striatus) in conjunction with spatially explicit habitat attribute data to evaluate the realized permeability of various landscape elements in a fragmented agricultural ecosystem. To accomplish this we 1) used logistic regression to evaluate whether land cover attributes were most often associated with the matrix between or habitat within genetically identified populations across the landscape, and 2) utilized spatially explicit habitat attribute data to predict genetically-derived Bayesian probabilities of population membership of individual chipmunks in an agricultural ecosystem. Consistency between the results of the two approaches with regard to facilitators and inhibitors of gene flow in the landscape indicate that this is a promising new way to utilize both landscape and genetic data to gain a deeper understanding of human-altered ecosystems.

Landscape genetics for the empirical assessment of resistance surfaces: the European pine marten (Martes martes) as a target-species of a regional ecological network

Coherent ecological networks (EN) composed of core areas linked by ecological corridors are being developed worldwide with the goal of promoting landscape connectivity and biodiversity conservation. However, empirical assessment of the performance of EN designs is critical to evaluate the utility of these networks to mitigate effects of habitat loss and fragmentation. Landscape genetics provides a particularly valuable framework to address the question of functional connectivity by providing a direct means to investigate the effects of landscape structure on gene flow. The goals of this study are (1) to evaluate the landscape features that drive gene flow of an EN target species (European pine marten), and (2) evaluate the optimality of a regional EN design in providing connectivity for this species within the Basque Country (North Spain). Using partial Mantel tests in a reciprocal causal modeling framework we competed 59 alternative models, including isolation by distance and the regional EN. Our analysis indicated that the regional EN was among the most supported resistance models for the pine marten, but was not the best supported model. Gene flow of pine marten in northern Spain is facilitated by natural vegetation, and is resisted by anthropogenic landcover types and roads. Our results suggest that the regional EN design being implemented in the Basque Country will effectively facilitate gene flow of forest dwelling species at regional scale.

Inferred vs realized patterns of gene flow: an analysis of population structure in the Andros Island Rock Iguana

Ecological data, the primary source of information on patterns and rates of migration, can be integrated with genetic data to more accurately describe the realized connectivity between geographically isolated demes. In this paper we implement this approach and discuss its implications for managing populations of the endangered Andros Island Rock Iguana, Cyclura cychlura cychlura. This iguana is endemic to Andros, a highly fragmented landmass of large islands and smaller cays. Field observations suggest that geographically isolated demes were panmictic due to high, inferred rates of gene flow. We expand on these observations using 16 polymorphic microsatellites to investigate the genetic structure and rates of gene flow from 188 Andros Iguanas collected across 23 island sites. Bayesian clustering of specimens assigned individuals to three distinct genotypic clusters. An analysis of molecular variance (AMOVA) indicates that allele frequency differences are responsible for a significant portion of the genetic variance across the three defined clusters (F_st =  0.117, p0.01). These clusters are associated with larger islands and satellite cays isolated by broad water channels with strong currents. These findings imply that broad water channels present greater obstacles to gene flow than was inferred from field observation alone. Additionally, rates of gene flow were indirectly estimated using BAYESASS 3.0. The proportion of individuals originating from within each identified cluster varied from 94.5 to 98.7%, providing further support for local isolation. Our assessment reveals a major disparity between inferred and realized gene flow. We discuss our results in a conservation perspective for species inhabiting highly fragmented landscapes.

Physiological vagility and its relationship to dispersal and neutral genetic heterogeneity in vertebrates

Vagility is the inherent power of movement by individuals. Vagility and the available duration of movement determine the dispersal distance individuals can move to interbreed, which affects the fine-scale genetic structure of vertebrate populations. Vagility and variation in population genetic structure are normally explained by geographic variation and not by the inherent power of movement by individuals. We present a new, quantitative definition for physiological vagility that incorporates aerobic capacity, body size, body temperature and the metabolic cost of transport, variables that are independent of the physical environment. Physiological vagility is the speed at which an animal can move sustainably based on these parameters. This meta-analysis tests whether this definition of physiological vagility correlates with empirical data for maximal dispersal distances and measured microsatellite genetic differentiation with distance {[F(ST)/[1-F(ST))]/ln distance} for amphibians, reptiles, birds and mammals utilizing three locomotor modes (running, flying, swimming). Maximal dispersal distance and physiological vagility increased with body mass for amphibians, reptiles and mammals utilizing terrestrial movement. The relative slopes of these relationships indicate that larger individuals require longer movement durations to achieve maximal dispersal distances. Both physiological vagility and maximal dispersal distance were independent of body mass for flying vertebrates. Genetic differentiation with distance was greatest for terrestrial locomotion, with amphibians showing the greatest mean and variance in differentiation. Flying birds, flying mammals and swimming marine mammals showed the least differentiation. Mean physiological vagility of different groups (class and locomotor mode) accounted for 98% of the mean variation in genetic differentiation with distance in each group. Genetic differentiation with distance was not related to body mass. The physiological capacity for movement (physiological vagility) quantitatively predicts genetic isolation by distance in the vertebrates examined.

Landscape influences on dispersal behaviour: a theoretical model and empirical test using the fire salamander, Salamandra infraimmaculata

When populations reside within a heterogeneous landscape, isolation by distance may not be a good predictor of genetic divergence if dispersal behaviour and therefore gene flow depend on landscape features. Commonly used approaches linking landscape features to gene flow include the least cost path (LCP), random walk (RW), and isolation by resistance (IBR) models. However, none of these models is likely to be the most appropriate for all species and in all environments. We compared the performance of LCP, RW and IBR models of dispersal with the aid of simulations conducted on artificially generated landscapes. We also applied each model to empirical data on the landscape genetics of the endangered fire salamander, Salamandra infraimmaculata, in northern Israel, where conservation planning requires an understanding of the dispersal corridors. Our simulations demonstrate that wide dispersal corridors of the low-cost environment facilitate dispersal in the IBR model, but inhibit dispersal in the RW model. In our empirical study, IBR explained the genetic divergence better than the LCP and RW models (partial Mantel correlation 0.413 for IBR, compared to 0.212 for LCP, and 0.340 for RW). Overall dispersal cost in salamanders was also well predicted by landscape feature slope steepness (76%), and elevation (24%). We conclude that fire salamander dispersal is well characterised by IBR predictions. Together with our simulation findings, these results indicate that wide dispersal corridors facilitate, rather than hinder, salamander dispersal. Comparison of genetic data to dispersal model outputs can be a useful technique in inferring dispersal behaviour from population genetic data.

Genetic diversity and population demography of the Chinese crocodile lizard (Shinisaurus crocodilurus) in China

The Chinese crocodile lizard Shinisaurus crocodilurus is a critically endangered species, listed in Appendix II of CITES. Its populations and habitat in China have undergone significant changes in recent years. Understanding the genetic variability and phylogeography of this species is very important for successful conservation. In this study, samples were taken from 11 wild ponds and two captive populations in China. We sequenced mitochondrial CYTB, partial ND6, and partial tRNA-Glu and genotyped 10 microsatellite loci. Our analyses of these data showed low genetic variability, no strong isolation caused by distance, and a lack of a phylogeographic structure in this species. Based on our results, the basal divergence between two clades of S. crocodilurus in China may have been caused by the formation of the Pearl River system. We found a population expansion in one of these clades. Microsatellite analysis indicated the presence of three clusters, separated by significant genetic differences. We found that most individuals in the two captive populations were from the Luokeng (Guangdong) and Guangxi wild source populations, respectively.

Conservation genetics of the endangered Spotted Turtle (Clemmys guttata) illustrate the risks of “bottleneck tests”

Studies of population genetics in turtles have suggested that turtles do not experience genetic impacts of bottlenecks as strongly as expected. However, recent studies cast doubt on two commonly used tests implemented in the program BOTTLENECK, suggesting that these findings should be re-evaluated. The Spotted Turtle (Clemmys guttata (Schneider, 1792)) is endangered both globally and within Canada, but genetic data required to develop effective recovery strategies are unavailable. Here, we conducted the first study of population genetic structure in C. guttata. We then used multiple small populations of C. guttata as replicates to test whether the commonly used program BOTTLENECK could detect the genetic signature of bottlenecks in our study populations, which are all thought to have experienced significant declines in the past 2–3 generations (75 years). Turtles (n = 256) were genotyped at 11 microsatellite loci. A suite of Bayesian population genetics analyses and a principal coordinates analysis identified a minimum of 6 distinct genetic populations and a maximum of 10 differentiated subpopulations across the sampled Canadian range of C. guttata, which corresponded to demographically independent units. BOTTLENECK failed to detect population declines. A literature review found that bottleneck tests in 17 of 18 previous genetic studies of tortoises and freshwater turtles were based on suboptimal sampling, potentially confounding their results. High retention of genetic diversity (allelic richness and heterozygosity) in isolated populations of C. guttata and other turtle species is encouraging for species recovery, but conclusions about the prevalence of genetic bottlenecks in such populations should be re-examined.

Population genetic structure of the endangered Butler’s Gartersnake (Thamnophis butleri): does the Short-headed Gartersnake (Thamnophis brachystoma) exist in Canada?

Understanding population genetic structure is fundamental to conservation of endangered species. It is particularly important when working with species that are morphologically conserved because strong genetic divisions could represent cryptic species. Butler’s Gartersnake (Thamnophis butleri (Cope, 1889)) is an endangered species in Canada, having a fragmented distribution and being restricted to southwestern Ontario. Furthermore, it is difficult to distinguish morphologically from a closely related species, the Short-headed Gartersnake (Thamnophis brachystoma (Cope, 1892)). We use mitochondrial DNA (mtDNA) and seven microsatellite DNA loci to evaluate the genetic structure of Canadian T. butleri populations and to test for the presence of T. brachystoma in one of these populations. All individuals had the same mtDNA haplotype, and there was no evidence of multiple, syntopic genetic clusters, thereby rejecting the hypothesis that T. butleri and T. brachystoma co-exist in Canada. Two different model-based assignment tests using microsatellite DNA data suggest that there are four to five genetically distinct clusters of T. butleri (FST from 0.12 to 0.20). We provide the first population genetic study of T. butleri in Canada and refute the presence of T. brachystoma. Our results may provide guidance on recovery strategies for this species and identify areas to target fine-scale genetic analyses.

The influence of landscape on gene flow in the eastern massasauga rattlesnake (Sistrurus c. catenatus): insight from computer simulations

Understanding how gene flow shapes contemporary population structure requires the explicit consideration of landscape composition and configuration. New landscape genetic approaches allow us to link such heterogeneity to gene flow within and among populations. However, the attribution of cause is difficult when landscape features are spatially correlated, or when genetic patterns reflect past events. We use spatial Bayesian clustering and landscape resistance analysis to identify the landscape features that influence gene flow across two regional populations of the eastern massasauga rattlesnake, Sistrurus c. catenatus. Based on spatially explicit simulations, we inferred how habitat distribution modulates gene flow and attempted to disentangle the effects of spatially confounded landscape features. We found genetic clustering across one regional landscape but not the other, and also local differences in the effect of landscape on gene flow. Beyond the effects of isolation-by-distance, water bodies appear to underlie genetic differentiation among individuals in one regional population. Significant effects of roads were additionally detected locally, but these effects are possibly confounded with the signal of water bodies. In contrast, we found no signal of isolation-by-distance or landscape effects on genetic structure in the other regional population. Our simulations imply that these local differences have arisen as a result of differences in population density or tendencies for juvenile rather than adult dispersal. Importantly, our simulations also demonstrate that the ability to detect the consequences of contemporary anthropogenic landscape features (e.g. roads) on gene flow may be compromised when long-standing natural features (e.g. water bodies) co-exist on the landscape.

Tracking climate change in a dispersal-limited species: reduced spatial and genetic connectivity in a montane salamander

Tropical montane taxa are often locally adapted to very specific climatic conditions, contributing to their lower dispersal potential across complex landscapes. Climate and landscape features in montane regions affect population genetic structure in predictable ways, yet few empirical studies quantify the effects of both factors in shaping genetic structure of montane-adapted taxa. Here, we considered temporal and spatial variability in climate to explain contemporary genetic differentiation between populations of the montane salamander, Pseudoeurycea leprosa. Specifically, we used ecological niche modelling (ENM) and measured spatial connectivity and gene flow (using both mtDNA and microsatellite markers) across extant populations of P. leprosa in the Trans-Mexican Volcanic Belt (TVB). Our results indicate significant spatial and genetic isolation among populations, but we cannot distinguish between isolation by distance over time or current landscape barriers as mechanisms shaping population genetic divergences. Combining ecological niche modelling, spatial connectivity analyses, and historical and contemporary genetic signatures from different classes of genetic markers allows for inference of historical evolutionary processes and predictions of the impacts future climate change will have on the genetic diversity of montane taxa with low dispersal rates. Pseudoeurycea leprosa is one montane species among many endemic to this region and thus is a case study for the continued persistence of spatially and genetically isolated populations in the highly biodiverse TVB of central Mexico.

Landscape genetics of a top neotropical predator

Habitat loss and fragmentation as a consequence of human activities is a worldwide phenomenon and one of the major threats to global biodiversity. Habitat loss and fragmentation is particularly a concern in the biodiverse tropics, where deforestation is occurring at unprecedented rates. Although insects are one of the most diverse and functionally important groups in tropical ecosystems, the quantitative effect of landscape features on their gene flow remains unknown. Here, we used a robust landscape genetics approach to quantify the effect of ten landscape features (deforestation, mature forests, other forest types, the River Chagres, streams, stream banks, roads, sea, lakes and swamps) and interactions between them, on the gene flow of a neotropical forest keystone species, the army ant Eciton burchellii. The influence of landscape on E. burchellii's gene flow reflected the different dispersal capability of its sexes; aerial for males and pedestrian for females, and the different depths of population history inferred from microsatellites and mitochondrial DNA. In contrast to the gene flow-facilitating effect of mature forests, deforested areas were found to be strong barriers for E. burchellii's gene flow. Other forest types were found to be gene flow facilitators but only when interacting with mature secondary forests, therefore indicating the importance of mature forests for the survival of E. burchelii and its associate species. The River Chagres was identified as a major historical gene flow barrier for E. burchellii, suggesting that an important loss of connectivity may occur because of large artificial waterways such as the Panama Canal.

Phylogeography and population genetic structure of the Ornate Dragon Lizard, Ctenophorus ornatus

Species inhabiting ancient, geologically stable landscapes that have been impacted by agriculture and urbanisation are expected to have complex patterns of genetic subdivision due to the influence of both historical and contemporary gene flow. Here, we investigate genetic differences among populations of the granite outcrop-dwelling lizard Ctenophorus ornatus, a phenotypically variable species with a wide geographical distribution across the south-west of Western Australia. Phylogenetic analysis of mitochondrial DNA sequence data revealed two distinct evolutionary lineages that have been isolated for more than four million years within the C. ornatus complex. This evolutionary split is associated with a change in dorsal colouration of the lizards from deep brown or black to reddish-pink. In addition, analysis of microsatellite data revealed high levels of genetic structuring within each lineage, as well as strong isolation by distance at multiple spatial scales. Among the 50 outcrop populations’ analysed, non-hierarchical Bayesian clustering analysis revealed the presence of 23 distinct genetic groups, with outcrop populations less than 4 km apart usually forming a single genetic group. When a hierarchical analysis was carried out, almost every outcrop was assigned to a different genetic group. Our results show there are multiple levels of genetic structuring in C. ornatus, reflecting the influence of both historical and contemporary evolutionary processes. They also highlight the need to recognise the presence of two evolutionarily distinct lineages when making conservation management decisions on this species.