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.     

Isolation by landscape in populations of a prized edible mushroom Tricholoma matsutake

Tricholoma matsutake, a wild edible ectomycorrhizal mushroom, is revered for its distinguished flavor and iconic significance. Here, we test for landscape effects on T. matsutake gene flow and population structure in the Eastern Himalayas. Using single-nucleotide polymorphic (SNP) DNA markers, isolation by distance patterns were tested on eight populations within and between watersheds. We find that high, treeless ridgelines are effective barriers to gene flow, even at distances less than 65 km, whereas populations located within watersheds are structured at greater distances. Mantel tests demonstrated a significant positive correlation between F _st and a “landscape distance” measured as the shortest distance between population pairs below treeline r = 0.574, P = 0.002, whereas strict euclidian distances do not correlate. AMOVA analysis revealed significant partitioning with 91% of the genetic variance found within populations and 7% found between watersheds, indicative of sexually recombining populations with limited gene flow between watersheds. We show that landscape is an important determinant of air-dispersed ectomycorrhizal species population structure in heterogeneous landscapes.     

Landscape genetic analyses reveal fine‐scale effects of forest fragmentation in an insular tropical bird

Within the framework of landscape genetics, resistance surface modelling is particularly relevant to explicitly test competing hypotheses about landscape effects on gene flow. To investigate how fragmentation of tropical forest affects population connectivity in a forest specialist bird species, we optimized resistance surfaces without a priori specification, using least‐cost (LCP) or resistance (IBR) distances. We implemented a two‐step procedure in order (i) to objectively define the landscape thematic resolution (level of detail in classification scheme to describe landscape variables) and spatial extent (area within the landscape boundaries) and then (ii) to test the relative role of several landscape features (elevation, roads, land cover) in genetic differentiation in the Plumbeous Warbler (Setophaga plumbea). We detected a small‐scale reduction of gene flow mainly driven by land cover, with a negative impact of the nonforest matrix on landscape functional connectivity. However, matrix components did not equally constrain gene flow, as their conductivity increased with increasing structural similarity with forest habitat: urban areas and meadows had the highest resistance values whereas agricultural areas had intermediate resistance values. Our results revealed a higher performance of IBR compared to LCP in explaining gene flow, reflecting suboptimal movements across this human‐modified landscape, challenging the common use of LCP to design habitat corridors and advocating for a broader use of circuit theory modelling. Finally, our results emphasize the need for an objective definition of landscape scales (landscape extent and thematic resolution) and highlight potential pitfalls associated with parameterization of resistance surfaces.     

Connectivity and gene flow among Eastern Tiger Salamander (Ambystoma tigrinum) populations in highly modified anthropogenic landscapes

Fragmented landscapes resulting from anthropogenic habitat modification can have significant impacts on dispersal, gene flow, and persistence of wildlife populations. Therefore, quantifying population connectivity across a mosaic of habitats in highly modified landscapes is critical for the development of conservation management plans for threatened populations. Endangered populations of the eastern tiger salamander (Ambystoma tigrinum) in New York and New Jersey are at the northern edge of the species’ range and remaining populations persist in highly developed landscapes in both states. We used landscape genetic approaches to examine regional genetic population structure and potential barriers to migration among remaining populations. Despite the post-glacial demographic processes that have shaped genetic diversity in tiger salamander populations at the northern extent of their range, we found that populations in each state belong to distinct genetic clusters, consistent with the large geographic distance that separates them. We detected overall low genetic diversity and high relatedness within populations, likely due to recent range expansion, isolation, and relatively small population sizes. Nonetheless, landscape connectivity analyses reveal habitat corridors among remaining breeding ponds. Furthermore, molecular estimates of population connectivity among ponds indicate that gene flow still occurs at regional scales. Further fragmentation of remaining habitat will potentially restrict dispersal among breeding ponds, cause the erosion of genetic diversity, and exacerbate already high levels of inbreeding. We recommend the continued management and maintenance of habitat corridors to ensure long-term viability of these endangered populations.     

Effects of slope and riparian habitat connectivity on gene flow in an endangered Panamanian frog, Atelopus varius

Aim Understanding how heterogeneous landscapes shape genetic structure not only sheds light on processes involved in population divergence and speciation, but can also guide management strategies to promote and maintain genetic connectivity of populations of endangered species. This study aimed to (1) identify barriers and corridors for gene flow among populations of the endangered frog, Atelopus varius and (2) assess the relative contributions of alternative landscape factors to patterns of genetic variation among these populations in a hypothesis testing framework. Location This study took place in western Panama and included all nine of the remaining known populations of A. varius at the time of study. Methods The influence of landscape variables on gene flow among populations was examined by testing for correlations between alternative landscape‐resistance scenarios and genetic distance. Fifteen alternative hypotheses about the influence of (1) riparian habitat corridors, (2) steep slopes, and (3) climatic suitability on patterns of genetic structure were tested in a causal modelling framework, using Mantel and partial‐Mantel tests, along with an analysis of molecular variation. Results Only the hypothesis attributing resistance to dispersal across steep slopes (genetic isolation by slope distance) was fully supported by the causal modelling approach. However, the analysis of molecular variance and the paths of least‐slope among populations suggest that riparian habitat connectivity may influence genetic structure as well. Main conclusions These results suggest that patterns of genetic variation among A. varius populations are affected by the slope of the landscape such that areas with steep slopes act as barriers to gene flow. In contrast, areas of low slope, such as streams and mountain ridges, appear to be important corridors for gene flow, especially among high elevation populations. These results engender important considerations for the management of this critically endangered species.     

Population structure of spotted salamanders (Ambystoma maculatum) in a fragmented landscape

Understanding the impacts of landscape-level processes on the population biology of amphibians is critical, especially for species inhabiting anthropogenically modified landscapes. Many pond-breeding amphibians are presumed to exist as metapopulations, but few studies demonstrate the extent and consequences of this metapopulation structure. Gene flow measures may facilitate the construction of more realistic models of population structure than direct measures of migration. This is especially true for species that are cryptic, such as many amphibians. We used eight polymorphic microsatellite loci to determine the genetic population structure of spotted salamanders ( Ambystoma maculatum ) breeding at 17 ponds in northeastern Ohio, a landscape fragmented by roads, agriculture, urban areas and the Cuyahoga River. Using a variety of analyses (Bayesian clustering, F -statistics, AMOVA) we generated a model of salamander population genetic structure. Our data revealed patterns of genetic connectivity that were not predicted by geographical distances between ponds (no isolation by distance). We also tested for a relationship between population structure and several indices of landscape resistance, but found no effect of potential barriers to dispersal on genetic connectivity. Strong overall connectivity among ponds, despite the hostile habitat matrix, may be facilitated by a network of riparian corridors associated with the Cuyahoga River; however, high gene flow in this system may indicate a general ability to disperse and colonize beyond particular corridors.     

Gene flow in complex landscapes: testing multiple hypotheses with causal modeling

Predicting population-level effects of landscape change depends on identifying factors that influence population connectivity in complex landscapes. However, most putative movement corridors and barriers have not been based on empirical data. In this study, we identify factors that influence connectivity by comparing patterns of genetic similarity among 146 black bears (Ursus americanus), sampled across a 3,000-km(2) study area in northern Idaho, with 110 landscape-resistance hypotheses. Genetic similarities were based on the pairwise percentage dissimilarity among all individuals based on nine microsatellite loci (average expected heterozygosity=0.79). Landscape-resistance hypotheses describe a range of potential relationships between movement cost and land cover, slope, elevation, roads, Euclidean distance, and a putative movement barrier. These hypotheses were divided into seven organizational models in which the influences of barriers, distance, and landscape features were statistically separated using partial Mantel tests. Only one of the competing organizational models was fully supported: patterns of genetic structure are primarily related to landscape gradients of land cover and elevation. The alternative landscape models, isolation by barriers and isolation by distance, are not supported. In this black bear population, gene flow is facilitated by contiguous forest cover at middle elevations.     

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.     

Population connectivity and genetic diversity of American marten (Martes americana) in the United States northern Rocky Mountains in a climate change context

Climate change is likely to alter population connectivity, particularly for species associated with higher elevation environments. The goal of this study is to predict the potential effects of future climate change on population connectivity and genetic diversity of American marten populations across a 30.2 million hectare region of the in the US northern Rocky Mountains. We use a landscape resistance model validated from empirical landscape genetics modeling to predict the current and expected future extent and fragmentation of American marten dispersal habitat under five climate change scenarios, corresponding to climatic warming of between 0.7 and 3.3 °C, consistent with expected climate change by year 2080. We predict the regions of the current and future landscapes where gene flow is expected to be governed by isolation by distance and the regions where population fragmentation is expected to limit gene flow. Finally, we predict changes in the strength and location of predicted movement corridors, fracture zones and the location of dispersal barriers across the study area in each scenario. We found that under the current climate, gene flow is predicted to be limited primarily by distance (isolation), and landscape structure does not significantly limit gene flow, resulting in very high genetic diversity over most of the study area. Projected climatic warming substantially reduces the extent and increases the fragmentation of marten populations in the western and northwestern parts of the study area. In contrast, climate change is not predicted to fragment the extensive higher elevation mountain massifs in central Idaho, the northern U.S. continental divide, and Greater Yellowstone Ecosystem. In addition, we show locations in the study area that are important corridors in the current landscape that remain intact across the climate change scenarios.     

Comparative landscape genetics of two river frog species occurring at different elevations on Mount Kilimanjaro

Estimating population connectivity and species' abilities to disperse across the landscape is crucial for understanding the long‐term persistence of species in changing environments. Surprisingly, few landscape genetic studies focused on tropical regions despite the alarming extinction rates within these ecosystems. Here, we compared the influence of landscape features on the distribution of genetic variation of an Afromontane frog, Amietia wittei, with that of its more broadly distributed lowland congener, Amietia angolensis, on Mt. Kilimanjaro, Tanzania. We predicted high gene flow in the montane species with movements enhanced through terrestrial habitats of the continuous rainforest. In contrast, dispersal might be restricted to aquatic corridors and reduced by anthropogenic disturbance in the lowland species. We found high gene flow in A. wittei relative to other montane amphibians. Nonetheless, gene flow was lower than in the lowland species which showed little population structure. Least‐cost path analysis suggested that dispersal is facilitated by stream networks in both species, but different landscape features were identified to influence connectivity among populations. Contrary to a previous study, gene flow in the lowland species was negatively correlated with the presence of human settlements. Also, genetic subdivision in A. wittei did not coincide with specific physical barriers as in other landscape genetic studies, suggesting that factors other than topography may contribute to population divergence. Overall, these results highlight the importance of a comparative landscape genetic approach for assessing the influence of the landscape matrix on population connectivity, particularly because nonintuitive results can alter the course of conservation and management.     

Landscape connectivity influences gene flow in a roe deer population inhabiting a fragmented landscape: an individual-based approach

Changes in agricultural practices and forest fragmentation can have a dramatic effect on landscape connectivity and the dispersal of animals, potentially reducing gene flow within populations. In this study, we assessed the influence of woodland connectivity on gene flow in a traditionally forest-dwelling species--the European roe deer--in a fragmented landscape. From a sample of 648 roe deer spatially referenced within a study area of 55 x 40 km, interindividual genetic distances were calculated from genotypes at 12 polymorphic microsatellite loci. We calculated two geographical distances between each pair of individuals: the Euclidean distance (straight line) and the 'least cost distance' (the trajectory that maximizes the use of wooded corridors). We tested the correlation between genetic pairwise distances and the two types of geographical pairwise distance using Mantel tests. The correlation was better using the least cost distance, which takes into account the distribution of wooded patches, especially for females (the correlation was stronger but not significant for males). These results suggest that in a fragmented woodland area roe deer dispersal is strongly linked to wooded structures and hence that gene flow within the roe deer population is influenced by the connectivity of the landscape.     

整合电路理论的生态廊道及其重要性识别

景观连接度被认为是影响诸多生态过程的一个重要因素.基于最小累积阻力模型的最小成本路径识别方法可以有效识别异质性景观中的功能连接,已被广泛应用到景观的功能连接评价与生态廊道模拟的研究中.基于电路理论的连接度模型用电阻代替了图论中的边、用电阻距离代替成本距离,来衡量异质性景观的功能连接.本文以SIMMAP 2.0软件生成的模拟景观为对象,借助于Linkage Mapper工具和Circuitscape软件,探讨如何将最小累积阻力模型与基于电路理论的连接度模型相结合来识别生态廊道及其景观要素的相对重要性.结果表明: 两种模型在应用中各有优势,互为补充.最小成本路径方法可以有效识别栖息地之间的最小成本廊道,基于电路理论的连接度模型通过电流密度的计算可以有效识别对景观连接性有重要影响的景观要素和“夹点”地区,并且“夹点”的位置不受廊道宽度的影响,在廊道重要性识别研究中具有明显优势.该方法可为区域生态保护规划和生态廊道设计提供科学依据.     

Urbanisation versus agriculture: a comparison of local genetic diversity and gene flow between wood mouse Apodemus sylvaticus populations in human‐modified landscapes

Urbanisation and agriculture dramatically modify the landscapes available for use by wildlife, affecting key aspects of their ecology such as survival, foraging, predation, competition and reproductive success. Relatively little is known about the effects of urbanisation and agriculture on the genetic structure, gene flow and genetic diversity of wild species. Here, landscape genetic techniques were applied to compare local genetic diversity and gene flow between wood mouse populations in urban and arable landscapes. Using nine microsatellite markers, individuals were genotyped from six arable and seven urban sample sites. Inter‐population genetic differentiation was significantly greater in urban than arable habitat, while allele richness, private allele richness and heterozygosity were higher for arable sample sites, with varying degrees of significance. These suggest that urban habitat was sufficiently fragmented to limit gene flow. To test the effect of landscape features on gene flow, several cost‐distance measures were generated. Overland distance and Euclidean distance correlated best with inter‐population genetic differentiation in arable habitat, whereas distances that accommodated differences in habitat quality better explained differentiation in urban habitat. There was no evidence that margins adjacent to roads, rivers or railways facilitated gene flow. Together, the results indicate that urban landscapes expose wood mice to greater fragmentation in habitat quality than arable areas, leading to greater population isolation that is not mitigated by the presence of dispersal corridors.     

Landscape Genetics of American Beaver in Coastal Oregon

American beaver (Castor canadensis) have been translocated for population restoration, reduction of human-wildlife conflict, and enhancement of ecosystem function. Yet few studies have assessed dispersal of beaver, making it difficult to determine at what scale translocations are appropriate. Genetic studies can provide inferences about gene flow, and thus dispersal. We used a landscape genetic approach to evaluate whether landscape features influenced gene flow among beaver in the Coast Range of western Oregon, USA, using samples collected April–September 2014. We collected genetic samples from live-captured (n = 232), road-killed (n = 2) and trapper-provided (n = 58) tissue samples and genotyped them at 10 microsatellite loci. We mapped records of beaver translocations into or within the study area during the twentieth century to consider the effect of those movements on genetic structure. We used population assignment tests to delimit genetic clusters, evaluated correspondence of those clusters with watershed boundaries and translocation history, and then estimated differentiation between clusters and between watersheds using model-based and model-free approaches. We evaluated how individual genetic differences varied with geographic distance, and investigated related pairs within clusters. We developed landscape resistance models incorporating slope, distance to water, and watershed boundaries at 2 scales, and estimated effective distances between sample locations with least cost path and circuit theoretic analyses. We evaluated the correlation of individual genetic distances with effective distances using a pseudo-bootstrapping approach. Landscape genetic models did not explain spatial variation in genetic structure better than geographic distance, but hierarchical genetic structure corresponded with watershed boundaries and suggested influences from historical translocations. Pairwise individual genetic distances were positively correlated with geographic distances to 61 km; highly-related pairs mostly were detected <1 km apart (median = 1.0 km, = 14.6 ± 2.3 [SE] km, n = 77). We concluded that slope and distance to water did not strongly limit dispersal and gene flow by beaver in this system, but concordance of genetic structure with watershed boundaries suggests that dispersal is more common within than between watersheds. Genetic differentiation of beaver within this topographically complex system was much greater than reported in a study at similar spatial scales in relatively flat topography. We recommend that translocation efforts of American beaver in topographically complex landscapes occur within watersheds when possible but conclude that dispersal can occur across watersheds. © 2021 The Wildlife Society. This article is a U.S. Government work and is in the public domain in the USA.     

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.     

Surprisingly little population genetic structure in a fungus‐associated beetle despite its exploitation of multiple hosts

In heterogeneous environments, landscape features directly affect the structure of genetic variation among populations by functioning as barriers to gene flow. Resource‐associated population genetic structure, in which populations that use different resources (e.g., host plants) are genetically distinct, is a well‐studied example of how environmental heterogeneity structures populations. However, the pattern that emerges in a given landscape should depend on its particular combination of resources. If resources constitute barriers to gene flow, population differentiation should be lowest in homogeneous landscapes, and highest where resources exist in equal proportions. In this study, we tested whether host community diversity affects population genetic structure in a beetle (Bolitotherus cornutus) that exploits three sympatric host fungi. We collected B. cornutus from plots containing the three host fungi in different proportions and quantified population genetic structure in each plot using a panel of microsatellite loci. We found no relationship between host community diversity and population differentiation in this species; however, we also found no evidence of resource‐associated differentiation, suggesting that host fungi are not substantial barriers to gene flow. Moreover, we detected no genetic differentiation among B. cornutus populations separated by several kilometers, even though a previous study demonstrated moderate genetic structure on the scale of a few hundred meters. Although we found no effect of community diversity on population genetic structure in this study, the role of host communities in the structuring of genetic variation in heterogeneous landscapes should be further explored in a species that exhibits resource‐associated population genetic structure.     

Landscape features affect gene flow of Scottish Highland red deer (Cervus elaphus)

Landscape features have been shown to strongly influence dispersal and, consequently, the genetic population structure of organisms. Studies quantifying the effect of landscape features on gene flow of large mammals with high dispersal capabilities are rare and have mainly been focused at large geographical scales. In this study, we assessed the influence of several natural and human-made landscape features on red deer gene flow in the Scottish Highlands by analysing 695 individuals for 21 microsatellite markers. Despite the relatively small scale of the study area (115 × 87 km), significant population structure was found using F -statistics ( F _ST = 0.019) and the program structure , with major differentiation found between populations sampled on either side of the main geographical barrier (the Great Glen). To assess the effect of landscape features on red deer population structure, the ArcMap GIS was used to create cost-distance matrices for moving between populations, using a range of cost values for each of the landscape features under consideration. Landscape features were shown to significantly affect red deer gene flow as they explained a greater proportion of the genetic variation than the geographical distance between populations. Sea lochs were found to be the most important red deer gene flow barriers in our study area, followed by mountain slopes, roads and forests. Inland lochs and rivers were identified as landscape features that might facilitate gene flow of red deer. Additionally, we explored the effect of choosing arbitrary cell cost values to construct least cost-distance matrices and described a method for improving the selection of cell cost values for a particular landscape feature.     

A comparison of regression methods for model selection in individual‐based landscape genetic analysis

Anthropogenic migration barriers fragment many populations and limit the ability of species to respond to climate‐induced biome shifts. Conservation actions designed to conserve habitat connectivity and mitigate barriers are needed to unite fragmented populations into larger, more viable metapopulations, and to allow species to track their climate envelope over time. Landscape genetic analysis provides an empirical means to infer landscape factors influencing gene flow and thereby inform such conservation actions. However, there are currently many methods available for model selection in landscape genetics, and considerable uncertainty as to which provide the greatest accuracy in identifying the true landscape model influencing gene flow among competing alternative hypotheses. In this study, we used population genetic simulations to evaluate the performance of seven regression‐based model selection methods on a broad array of landscapes that varied by the number and type of variables contributing to resistance, the magnitude and cohesion of resistance, as well as the functional relationship between variables and resistance. We also assessed the effect of transformations designed to linearize the relationship between genetic and landscape distances. We found that linear mixed effects models had the highest accuracy in every way we evaluated model performance; however, other methods also performed well in many circumstances, particularly when landscape resistance was high and the correlation among competing hypotheses was limited. Our results provide guidance for which regression‐based model selection methods provide the most accurate inferences in landscape genetic analysis and thereby best inform connectivity conservation actions.     

Ancient trade routes shaped the genetic structure of horses in eastern Eurasia

Animal exchange networks have been shown to play an important role in determining gene flow among domestic animal populations. The Silk Road is one of the oldest continuous exchange networks in human history, yet its effectiveness in facilitating animal exchange across large geographical distances and topographically challenging landscapes has never been explicitly studied. Horses are known to have been traded along the Silk Roads; however, extensive movement of horses in connection with other human activities may have obscured the genetic signature of the Silk Roads. To investigate the role of the Silk Roads in shaping the genetic structure of horses in eastern Eurasia, we analysed microsatellite genotyping data from 455 village horses sampled from 17 locations. Using least‐cost path methods, we compared the performance of models containing the Silk Roads as corridors for gene flow with models containing single landscape features. We also determined whether the recent isolation of former Soviet Union countries from the rest of Eurasia has affected the genetic structure of our samples. The overall level of genetic differentiation was low, consistent with historically high levels of gene flow across the study region. The spatial genetic structure was characterized by a significant, albeit weak, pattern of isolation by distance across the continent with no evidence for the presence of distinct genetic clusters. Incorporating landscape features considerably improved the fit of the data; however, when we controlled for geographical distance, only the correlation between genetic differentiation and the Silk Roads remained significant, supporting the effectiveness of this ancient trade network in facilitating gene flow across large geographical distances in a topographically complex landscape.