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.     

The influence of breeding phenology on the genetic structure of four pond‐breeding salamanders

Understanding metapopulation dynamics requires knowledge about local population dynamics and movement in both space and time. Most genetic metapopulation studies use one or two study species across the same landscape to infer population dynamics; however, using multiple co‐occurring species allows for testing of hypotheses related to different life history strategies. We used genetic data to study dispersal, as measured by gene flow, in three ambystomatid salamanders (Ambystoma annulatum , A. maculatum , and A. opacum ) and the Central Newt (Notophthalmus viridescens louisianensis ) on the same landscape in Missouri, USA . While all four salamander species are forest dependent organisms that require fishless ponds to reproduce, they differ in breeding phenology and spatial distribution on the landscape. We use these differences in life history and distribution to address the following questions: (1) Are there species‐level differences in the observed patterns of genetic diversity and genetic structure? and (2) Is dispersal influenced by landscape resistance? We detected two genetic clusters in A. annulatum and A. opacum on our landscape; both species breed in the fall and larvae overwinter in ponds. In contrast, no structure was evident in A. maculatum and N. v. louisianensis , species that breed during the spring. Tests for isolation by distance were significant for the three ambystomatids but not for N. v. louisianensis . Landscape resistance also contributed to genetic differentiation for all four species. Our results suggest species‐level differences in dispersal ability and breeding phenology are driving observed patterns of genetic differentiation. From an evolutionary standpoint, the observed differences in dispersal distances and genetic structure between fall breeding and spring breeding species may be a result of the trade‐off between larval period length and size at metamorphosis which in turn may influence the long‐term viability of the metapopulation. Thus, it is important to consider life history differences among closely related and ecologically similar species when making management decisions.     

Landscape resistance and habitat combine to provide an optimal model of genetic structure and connectivity at the range margin of a small mammal

We evaluated the effect of habitat and landscape characteristics on the population genetic structure of the white‐footed mouse. We develop a new approach that uses numerical optimization to define a model that combines site differences and landscape resistance to explain the genetic differentiation between mouse populations inhabiting forest patches in southern Québec. We used ecological distance computed from resistance surfaces with Circuitscape to infer the effect of the landscape matrix on gene flow. We calculated site differences using a site index of habitat characteristics. A model that combined site differences and resistance distances explained a high proportion of the variance in genetic differentiation and outperformed models that used geographical distance alone. Urban and agriculture‐related land uses were, respectively, the most and the least resistant landscape features influencing gene flow. Our method detected the effect of rivers and highways as highly resistant linear barriers. The density of grass and shrubs on the ground best explained the variation in the site index of habitat characteristics. Our model indicates that movement of white‐footed mouse in this region is constrained along routes of low resistance. Our approach can generate models that may improve predictions of future northward range expansion of this small mammal.     

Urban land use limits regional bumble bee gene flow

Potential declines in native pollinator communities and increased reliance on pollinator‐dependent crops have raised concerns about native pollinator conservation and dispersal across human‐altered landscapes. Bumble bees are one of the most effective native pollinators and are often the first to be extirpated in human‐altered habitats, yet little is known about how bumble bees move across fine spatial scales and what landscapes promote or limit their gene flow. In this study, we examine regional genetic differentiation and fine‐scale relatedness patterns of the yellow‐faced bumble bee, Bombus vosnesenskii, to investigate how current and historic habitat composition impact gene flow. We conducted our study across a landscape mosaic of natural, agricultural and urban/suburban habitats, and we show that B. vosnesenskii exhibits low but significant levels of differentiation across the study system (F_ST = 0.019, D_est = 0.049). Most importantly, we reveal significant relationships between pairwise F_ST and resistance models created from contemporary land use maps. Specifically, B. vosnesenskii gene flow is most limited by commercial, industrial and transportation‐related impervious cover. Finally, our fine‐scale analysis reveals significant but declining relatedness between individuals at the 1–9 km spatial scale, most likely due to local queen dispersal. Overall, our results indicate that B. vosnesenskii exhibits considerable local dispersal and that regional gene flow is significantly limited by impervious cover associated with urbanization.     

Landscape effects on gene flow for a climate‐sensitive montane species,the American pika

Climate change is arguably the greatest challenge to conservation of our time. Most vulnerability assessments rely on past and current species distributions to predict future persistence but ignore species' abilities to disperse through landscapes, which may be particularly important in fragmented habitats and crucial for long‐term persistence in changing environments. Landscape genetic approaches explore the interactions between landscape features and gene flow and can clarify how organisms move among suitable habitats, but have suffered from methodological uncertainties. We used a landscape genetic approach to determine how landscape and climate‐related features influence gene flow for American pikas (Ochotona princeps) in Crater Lake National Park. Pikas are heat intolerant and restricted to cool microclimates; thus, range contractions have been predicted as climate changes. We evaluated the correlation between landscape variables and genetic distance using partial Mantel tests in a causal modelling framework, and used spatially explicit simulations to evaluate methods of model optimization including a novel approach based on relative support and reciprocal causal modelling. We found that gene flow was primarily restricted by topographic relief, water and west‐facing aspects, suggesting that physical restrictions related to small body size and mode of locomotion, as well as exposure to relatively high temperatures, limit pika dispersal in this alpine habitat. Our model optimization successfully identified landscape features influencing resistance in the simulated data for this landscape, but underestimated the magnitude of resistance. This is the first landscape genetic study to address the fundamental question of what limits dispersal and gene flow in the American pika.     

A multiscale analysis of gene flow for the New England cottontail,an imperiled habitat specialist in a fragmented landscape

Landscape features of anthropogenic or natural origin can influence organisms' dispersal patterns and the connectivity of populations. Understanding these relationships is of broad interest in ecology and evolutionary biology and provides key insights for habitat conservation planning at the landscape scale. This knowledge is germane to restoration efforts for the New England cottontail (Sylvilagus transitionalis), an early successional habitat specialist of conservation concern. We evaluated local population structure and measures of genetic diversity of a geographically isolated population of cottontails in the northeastern United States. We also conducted a multiscale landscape genetic analysis, in which we assessed genetic discontinuities relative to the landscape and developed several resistance models to test hypotheses about landscape features that promote or inhibit cottontail dispersal within and across the local populations. Bayesian clustering identified four genetically distinct populations, with very little migration among them, and additional substructure within one of those populations. These populations had private alleles, low genetic diversity, critically low effective population sizes (3.2–36.7), and evidence of recent genetic bottlenecks. Major highways and a river were found to limit cottontail dispersal and to separate populations. The habitat along roadsides, railroad beds, and utility corridors, on the other hand, was found to facilitate cottontail movement among patches. The relative importance of dispersal barriers and facilitators on gene flow varied among populations in relation to landscape composition, demonstrating the complexity and context dependency of factors influencing gene flow and highlighting the importance of replication and scale in landscape genetic studies. Our findings provide information for the design of restoration landscapes for the New England cottontail and also highlight the dual influence of roads, as both barriers and facilitators of dispersal for an early successional habitat specialist in a fragmented landscape.     

Ecological resistance surfaces predict fine‐scale genetic differentiation in a terrestrial woodland salamander

Landscape genetics has seen tremendous advances since its introduction, but parameterization and optimization of resistance surfaces still poses significant challenges. Despite increased availability and resolution of spatial data, few studies have integrated empirical data to directly represent ecological processes as genetic resistance surfaces. In our study, we determine the landscape and ecological factors affecting gene flow in the western slimy salamander (Plethodon albagula). We used field data to derive resistance surfaces representing salamander abundance and rate of water loss through combinations of canopy cover, topographic wetness, topographic position, solar exposure and distance from ravine. These ecologically explicit composite surfaces directly represent an ecological process or physiological limitation of our organism. Using generalized linear mixed‐effects models, we optimized resistance surfaces using a nonlinear optimization algorithm to minimize model AIC. We found clear support for the resistance surface representing the rate of water loss experienced by adult salamanders in the summer. Resistance was lowest at intermediate levels of water loss and higher when the rate of water loss was predicted to be low or high. This pattern may arise from the compensatory movement behaviour of salamanders through suboptimal habitat, but also reflects the physiological limitations of salamanders and their sensitivity to extreme environmental conditions. Our study demonstrates that composite representations of ecologically explicit processes can provide novel insight and can better explain genetic differentiation than ecologically implicit landscape resistance surfaces. Additionally, our study underscores the fact that spatial estimates of habitat suitability or abundance may not serve as adequate proxies for describing gene flow, as predicted abundance was a poor predictor of genetic differentiation.     

Divergent landscape effects on population connectivity in two co-occurring amphibian species

The physical and environmental attributes of landscapes often shape patterns of population connectivity by influencing dispersal and gene flow. Landscape effects on movement are typically evaluated for single species. However, inferences from multiple species are required for multi‐species management strategies increasingly being applied in conservation. In this study, I compared the spatial genetic patterns of two amphibian species across the northeastern United States and estimated the influence of specific landscape features on the observed genetic structure. The spotted salamander (Ambystoma maculatum) and wood frog (Rana sylvatica) share many ecological attributes related to habitat use, phenology and site fidelity. However, I hypothesized that important differences in their movement patterns and life history would create distinct genetic patterns for each species. Using 14 microsatellite loci, I tested for differences in the level of genetic differentiation between the two species across 22 breeding ponds. The effects of eight landscape features were also estimated by evaluating 32 landscape resistance models. Spotted salamanders exhibited significantly higher genetic differentiation than wood frogs. Different landscape features were also identified as potential drivers of the genetic patterns in each species, with little overlap in model support between species. Collectively, these results provide strong evidence that these two amphibian species interact with the landscape in measurably different ways. The distinct genetic patterns observed are consistent with key differences in movement ability and life history between A. maculatum and R. sylvatica. These results highlight the importance of considering more than one species when assessing the impacts of the landscape matrix on population connectivity, even for ecologically similar species within the same habitats.     

NONSTATIONARY PATTERNS OF ISOLATION‐BY‐DISTANCE: INFERRING MEASURES OF LOCAL GENETIC DIFFERENTIATION WITH BAYESIAN KRIGING

Patterns of isolation‐by‐distance (IBD) arise when population differentiation increases with increasing geographic distances. Patterns of IBD are usually caused by local spatial dispersal, which explains why differences of allele frequencies between populations accumulate with distance. However, spatial variations of demographic parameters such as migration rate or population density can generate nonstationary patterns of IBD where the rate at which genetic differentiation accumulates varies across space. To characterize nonstationary patterns of IBD, we infer local genetic differentiation based on Bayesian kriging. Local genetic differentiation for a sampled population is defined as the average genetic differentiation between the sampled population and fictive neighboring populations. To avoid defining populations in advance, the method can also be applied at the scale of individuals making it relevant for landscape genetics. Inference of local genetic differentiation relies on a matrix of pairwise similarity or dissimilarity between populations or individuals such as matrices of between pairs of populations. Simulation studies show that maps of local genetic differentiation can reveal barriers to gene flow but also other patterns such as continuous variations of gene flow across habitat. The potential of the method is illustrated with two datasets: single nucleotide polymorphisms from human Swedish populations and dominant markers for alpine plant species.     

Contemporary human‐altered landscapes and oceanic barriers reduce bumble bee gene flow

Much of the world's terrestrial landscapes are being altered by humans in the form of agriculture, urbanization and pastoral systems, with major implications for biodiversity. Bumble bees are one of the most effective pollinators in both natural and cultivated landscapes, but are often the first to be extirpated in human‐altered habitats. Yet, little is known about the role of natural and human‐altered habitats in promoting or limiting bumble bee gene flow. In this study, I closely examine the genetic structure of the yellow‐faced bumble bee, Bombus vosnesenskii, across the southwestern US coast and find strong evidence that natural oceanic barriers, as well as contemporary human‐altered habitats, limit bee gene flow. Heterozygosity and allelic richness were lower in island populations, while private allelic richness was higher in island populations compared to mainland populations. Genetic differentiation, measured for three indices across the 1000 km study region, was significantly greater than the null expectation (F_ST = 0.041, F’_ST = 0.044 and D_est = 0.155) and correlated with geographic distance. Furthermore, genetic differentiation patterns were most strongly correlated with contemporary (2011) not past (2006, 2001) resistance maps calibrated for high dispersal limitation over oceans, impervious habitat and croplands. Despite the incorporation of dramatic elevation gradients, the analyses reveal that oceans and contemporary human land use, not mountains, are the primary dispersal barriers for B. vosnesenskii gene flow. These findings reinforce the importance of maintaining corridors of suitable habitat across the distribution range of native pollinators to promote their persistence and safeguard their ability to provide essential pollination services.     

Genetic evidence for landscape effects on dispersal in the army ant Eciton burchellii

Inhibited dispersal, leading to reduced gene flow, threatens populations with inbreeding depression and local extinction. Fragmentation may be especially detrimental to social insects because inhibited gene flow has important consequences for cooperation and competition within and among colonies. Army ants have winged males and permanently wingless queens; these traits imply male‐biased dispersal. However, army ant colonies are obligately nomadic and have the potential to traverse landscapes. Eciton burchellii, the most regularly nomadic army ant, is a forest interior species: colony raiding activities are limited in the absence of forest cover. To examine whether nomadism and landscape (forest clearing and elevation) affect population genetic structure in a montane E. burchellii population, we reconstructed queen and male genotypes from 25 colonies at seven polymorphic microsatellite loci. Pairwise genetic distances among individuals were compared to pairwise geographical and resistance distances using regressions with permutations, partial Mantel tests and random forests analyses. Although there was no significant spatial genetic structure in queens or males in montane forest, dispersal may be male‐biased. We found significant isolation by landscape resistance for queens based on land cover (forest clearing), but not on elevation. Summed colony emigrations over the lifetime of the queen may contribute to gene flow in this species and forest clearing impedes these movements and subsequent gene dispersal. Further forest cover removal may increasingly inhibit Eciton burchellii colony dispersal. We recommend maintaining habitat connectivity in tropical forests to promote population persistence for this keystone species.     

Dispersal ability and habitat requirements determine landscape‐level genetic patterns in desert aquatic insects

Species occupying the same geographic range can exhibit remarkably different population structures across the landscape, ranging from highly diversified to panmictic. Given limitations on collecting population‐level data for large numbers of species, ecologists seek to identify proximate organismal traits—such as dispersal ability, habitat preference and life history—that are strong predictors of realized population structure. We examined how dispersal ability and habitat structure affect the regional balance of gene flow and genetic drift within three aquatic insects that represent the range of dispersal abilities and habitat requirements observed in desert stream insect communities. For each species, we tested for linear relationships between genetic distances and geographic distances using Euclidean and landscape‐based metrics of resistance. We found that the moderate‐disperser Mesocapnia arizonensis (Plecoptera: Capniidae) has a strong isolation‐by‐distance pattern, suggesting migration–drift equilibrium. By contrast, population structure in the flightless Abedus herberti (Hemiptera: Belostomatidae) is influenced by genetic drift, while gene flow is the dominant force in the strong‐flying Boreonectes aequinoctialis (Coleoptera: Dytiscidae). The best‐fitting landscape model for M. arizonensis was based on Euclidean distance. Analyses also identified a strong spatial scale‐dependence, where landscape genetic methods only performed well for species that were intermediate in dispersal ability. Our results highlight the fact that when either gene flow or genetic drift dominates in shaping population structure, no detectable relationship between genetic and geographic distances is expected at certain spatial scales. This study provides insight into how gene flow and drift interact at the regional scale for these insects as well as the organisms that share similar habitats and dispersal abilities.     

Landscape genetics of an early successional specialist in a disturbance‐prone environment

Species that specialize in disturbed habitats may have considerably different dispersal strategies than those adapted to more stable environments. However, little is known of the dispersal patterns and population structure of such species. This information is important for conservation because many postfire specialists are at risk from anthropogenic changes to natural disturbance regimes. We used microsatellite markers to assess the effect of landscape variation and recent disturbance history on dispersal by a small mammal species that occupies the early seral stage of vegetation regeneration in burnt environments. We predicted that a postfire specialist would be able to disperse over multiple habitat types (generalist) and not exhibit sex‐biased dispersal, as such strategies should enable effective colonization of spatially and temporally variable habitat. We found significant differentiation between sites that fitted an isolation‐by‐distance pattern and spatial autocorrelation of multilocus genotypes to a distance of 2–3 km. There was no consistent genetic evidence for sex‐biased dispersal. We tested the influence of different habitat‐ and fire‐specific landscape resistance scenarios on genetic distance between individuals and found a significant effect of fire. Our genetic data supported recently burned vegetation having greater conductance for gene flow than unburnt habitat, but variation in habitat quality between vegetation types and occupied patches had no effect on gene flow. Postfire specialists must evolve an effective dispersal ability to move over distances that would ensure access to early successional stage vegetation. Natural disturbance and natural heterogeneity may therefore not influence population genetic structure as negatively as expected.     

Habitat use,but not gene flow,is influenced by human activities in two ecotypes of Egyptian fruit bat (Rousettus aegyptiacus)

Understanding the ecological, behavioural and evolutionary response of organisms to changing environments is of primary importance in a human‐altered world. It is crucial to elucidate how human activities alter gene flow and what are the consequences for the genetic structure of a species. We studied two lineages of the Egyptian fruit bat (Rousettus aegyptiacus) throughout the contact zone between mesic and arid Ecozones in the Middle East to evaluate the species' response to the growing proportion of human‐altered habitats in the desert. We integrated population genetics, morphometrics and movement ecology to analyse population structure, morphological variation and habitat use from GPS‐ or radio‐tagged individuals from both desert and Mediterranean areas. We classified the spatial distribution and environmental stratification by describing physical–geographical conditions and land cover. We analysed this information to estimate patch occupancy and used an isolation‐by‐resistance approach to model gene flow patterns. Our results suggest that lineages from desert and Mediterranean habitats, despite their admixture, are isolated by environment and by adaptation supporting their classification as ecotypes. We found a positive effect of human‐altered habitats on patch occupancy and habitat use of fruit bats by increasing the availability of roosting and foraging areas. While this commensalism promotes the distribution of fruit bats throughout the Middle East, gene flow between colonies has not been altered by human activities. This discrepancy between habitat use and gene flow patterns may, therefore, be explained by the breeding system of the species and modifications of natal dispersal patterns.     

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

Dispersal is a fundamental attribute of species in nature and shapes population dynamics, evolutionary trajectories and genetic variation across spatial and temporal scales. It is increasingly clear that landscape features have large impacts on dispersal patterns. Thus, understanding how individuals and species move through landscapes is essential for predicting impacts of landscape alterations. Information on dispersal patterns, however, is lacking for many taxa, particularly reptiles. Eastern foxsnakes (Mintoinus gloydi) are marsh and prairie specialists that avoid agricultural fields, but they have persisted across a fragmented region in southwestern Ontario and northern Ohio. Here, we combined habitat suitability modelling with population genetic analyses to infer how foxsnakes disperse through a habitat mosaic of natural and altered landscape features. Boundary regions between the eight genetic clusters, identified through assignment tests, were comprised of low suitability habitat (e.g. agricultural fields). Island populations were grouped into a single genetic cluster, and comparatively low F(ST) values between island and mainland populations suggest open water presents less of a barrier than nonsuitable terrestrial habitat. Isolation by resistance and least-cost path analysis produced similar results with matrices of pairwise individual genetic distance significantly more correlated to matrices of resistance values derived from habitat suitability than models with an undifferentiated landscape. Spatial autocorrelation results matched better with assignment results when incorporating resistance values rather than straight-line distances. All analyses used in our study produced similar results suggesting that habitat degradation limits dispersal for foxsnakes, which has had a strong effect on the genetic population structure across this region.     

Evolutionary and demographic processes shaping geographic patterns of genetic diversity in a keystone species,the African forest elephant (Loxodonta cyclotis)

The past processes that have shaped geographic patterns of genetic diversity may be difficult to infer from current patterns. However, in species with sex differences in dispersal, differing phylogeographic patterns between mitochondrial (mt) and nuclear (nu) DNA may provide contrasting insights into past events. Forest elephants (Loxodonta cyclotis) were impacted by climate and habitat change during the Pleistocene, which likely shaped phylogeographic patterns in mitochondrial (mt) DNA that have persisted due to limited female dispersal. By contrast, the nuclear (nu) DNA phylogeography of forest elephants in Central Africa has not been determined. We therefore examined the population structure of Central African forest elephants by genotyping 94 individuals from six localities at 21 microsatellite loci. Between forest elephants in western and eastern Congolian forests, there was only modest genetic differentiation, a pattern highly discordant with that of mtDNA. Nuclear genetic patterns are consistent with isolation by distance. Alternatively, male‐mediated gene flow may have reduced the previous regional differentiation in Central Africa suggested by mtDNA patterns, which likely reflect forest fragmentation during the Pleistocene. In species like elephants, male‐mediated gene flow erases the nuclear genetic signatures of past climate and habitat changes, but these continue to persist as patterns in mtDNA because females do not disperse. Conservation implications of these results are discussed.     

Effects of landscape and history on diversification of a montane, stream-breeding amphibian

Aim The aim of this study was to understand the roles of landscape features in shaping patterns of contemporary and historical genetic diversification among populations of the Andean tree frog (Hypsiboas andinus) across spatial scales. Location Andes mountains, north‐western Argentina, South America. Methods Mitochondrial DNA control region sequences were utilized to assess genetic differentiation among populations and calculate population pair‐wise genetic distances. Three models of movement, namely traditional straight‐line distance and two effective distances based on habitat classification, were examined to determine which of these explained the most variation in pair‐wise population genetic differentiation. The two habitat classifications were based on digital vegetation and hydrology layers that were generated from a 90‐m resolution digital elevation model (DEM) and known relationships between elevation and habitat. Mantel tests were conducted to test for correlations between geographic and genetic distance matrices and to estimate the percentage variation explained by each type of geographic distance. To investigate the location of possible barriers to gene flow, we used Monmonier’s maximum difference algorithm as implemented in barrier 2.2. Results At both geographic scales, effective distances explained more variation in genetic differentiation than did straight‐line distance. The least‐cost distances based on the simple classification performed better than the more detailed habitat classification. We controlled for the effects of historical range fragmentation determined from previous nested clade analyses, and therefore evaluated the effect of different distances on the genetic variation attributable to more recent factors. Effective distances identified populations that were highly divergent as a result of isolation in unsuitable habitats. The proposed locations of barriers to gene flow identified using Monmonier’s maximum difference algorithm corresponded well with earlier analyses and supported findings from our partial Mantel tests. Main conclusions Our results indicate that landscape features have been important in both historical and contemporary genetic structuring of populations of H. andinus at both large and small spatial scales. A landscape genetic perspective offers novel insights not provided by traditional phylogeographic studies: (1) effective distances can better explain patterns of differentiation in populations, especially in heterogeneous landscapes where barriers to dispersal may be common; and (2) least‐cost path analysis can help to identify corridors of movement between populations that are biologically more realistic.     

Functional connectivity in replicated urban landscapes in the land snail (Cornu aspersum)

Urban areas are highly fragmented and thereby exert strong constraints on individual dispersal. Despite this, some species manage to persist in urban areas, such as the garden snail, Cornu aspersum, which is common in cityscapes despite its low mobility. Using landscape genetic approaches, we combined study area replication and multiscale analysis to determine how landscape composition, configuration and connectivity influence snail dispersal across urban areas. At the overall landscape scale, areas with a high percentage of roads decreased genetic differentiation between populations. At the population scale, genetic differentiation was positively linked with building surface, the proportion of borders where wooded patches and roads appeared side by side and the proportion of borders combining wooded patches and other impervious areas. Analyses based on pairwise genetic distances validated the isolation‐by‐distance and isolation‐by‐resistance models for this land snail, with an equal fit to least‐cost paths and circuit‐theory‐based models. Each of the 12 landscapes analysed separately yielded specific relations to environmental features, whereas analyses integrating all replicates highlighted general common effects. Our results suggest that urban transport infrastructures facilitate passive snail dispersal. At a local scale, corresponding to active dispersal, unfavourable habitats (wooded and impervious areas) isolate populations. This work upholds the use of replicated landscapes to increase the generalizability of landscape genetics results and shows how multiscale analyses provide insight into scale‐dependent processes.     

The use and misuse of regression models in landscape genetic analyses

The field of landscape genetics has been rapidly evolving, adopting and adapting analytical frameworks to address research questions. Current studies are increasingly using regression‐based frameworks to infer the individual contributions of landscape and habitat variables on genetic differentiation. This paper outlines appropriate and inappropriate uses of multiple regression for these purposes, and demonstrates through simulation the limitations of different analytical frameworks for making correct inference. Of particular concern are recent studies seeking to explain genetic differences by fitting regression models with effective distance variables calculated independently on separate landscape resistance surfaces. When moving across the landscape, organisms cannot respond independently and uniquely to habitat and landscape features. Analyses seeking to understand how landscape features affect gene flow should model a single conductance or resistance surface as a parameterized function of relevant spatial covariates, and estimate the values of these parameters by linking a single set of resistance distances to observed genetic dissimilarity via a loss function. While this loss function may involve a regression‐like step, the associated nuisance parameters are not interpretable in terms of organismal movement and should not be conflated with what is actually of interest: the mapping between spatial covariates and conductance/resistance. The growth and evolution of landscape genetics as a field has been rapid and exciting. It is the goal of this paper to highlight past missteps and demonstrate limitations of current approaches to ensure that future use of regression models will appropriately consider the process being modeled, which will provide clarity to model interpretation.     

Limited influence of local and landscape factors on finescale gene flow in two pond‐breeding amphibians

Dispersal and gene flow within animal populations are influenced by the composition and configuration of the landscape. In this study, we evaluated hypotheses about the impact of natural and anthropogenic factors on genetic differentiation in two amphibian species, the spotted salamander (Ambystoma maculatum) and the wood frog (Lithobates sylvaticus) in a commercial forest in central Maine. We conducted this analysis at two scales: a local level, focused on factors measured at each breeding pond, and a landscape level, focused on factors measured between ponds. We investigated the effects of a number of environmental factors in six categories including Productivity, Physical, Land Composition, Land Configuration, Isolation and Location. Embryos were sampled from 56 spotted salamander breeding ponds and 39 wood frog breeding ponds. We used a hierarchical Bayesian approach in the program GESTE at each breeding pond and a random forest algorithm in conjunction with a network analysis between the ponds. We found overall high genetic connectivity across distances up to 17 km for both species and a limited effect of natural and anthropogenic factors on gene flow. We found the null models best explained patterns of genetic differentiation at a local level and found several factors at the landscape level that weakly influenced gene flow. This research indicates multiscale investigations that incorporate local and landscape factors are valuable for understanding patterns of gene flow. Our findings suggest that dispersal rates in this system are high enough to minimize genetic structuring and that current forestry practices do not significantly impede dispersal.