Zooplankton on the move: first results on the quantification of dispersal of zooplankton in a set of interconnected ponds

In systems of interconnected ponds or lakes, the dispersal of zooplankton may be mediated by the active population component, with rivulets and overflows functioning as dispersal pathways and the dispersal being unidirectional. Such systems offer the possibility to study the impact of dispersal rates on local population dynamics and community structure, and provide opportunities to quantify dispersal in the field in a straightforward manner. In this study, dispersal of active zooplankton populations among interconnected ponds was quantified directly in the field by sampling the small waterways connecting the ponds. The number of dispersing zooplankton sampled in connecting elements was on average high (almost 7000 ind h^–1). However, the contribution of dispersing individuals to total population size in the target ponds was very limited (< 1% 24 h^–1.). Only a weak diel pattern in dispersal rates was observed.     

Biological Invasions Across Spatial Scales: Intercontinental, Regional, and Local Dispersal of Cladoceran Zooplankton

The frequency of dispersal of invertebrates among lakes depends upon perspective and spatial scale. Effective passive dispersal requires both the transport of propagules and the establishment of populations large enough to be detected. At a global scale, biogeographic patterns of cladoceran zooplankton species suggest that effective dispersal among continents was originally rare, but greatly increased in the past century with expanded commerce. Genetic analysis allows some reconstruction of past dispersal events. Allozyme and mitochondrial DNA comparisons among New World and Old-World populations of several exotic cladocerans have provided estimates for likely source populations of colonists, their dispersal corridors, and timing of earlier dispersal events. Detecting the Old-World tropical exotic Daphnia lumholtzi early in its invasion of North America has allowed detailed analysis of its spatial spread. Twelve years of collection records indicate a rapid invasion of reservoirs in the United States, by both regional spread and long-distance jumps to new regions. Combining landscape features with zooplankton surveys from south-central US reservoirs revealed higher colonization rates of D. lumholtzi at lower landscape positions, a result which can be explained by either greater propagule load or by higher susceptibility of these downstream reservoirs. Because invaded reservoirs provide a source of propagules for nearby floodplain ponds, the rarity of this species in ponds suggests limitation by local environments. Such analyses of invading species over multiple spatial scales allow a better understanding of ecological processes governing invasion dynamics.     

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.     

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.     

Influence of landscape on the population genetic structure of the alpine butterfly parnassius smintheus (Papilionidae)

Four microsatellite DNA markers were developed which were used to examine the relationship between landscape and population genetic structure among a set of populations of the butterfly Parnassius smintheus located in the foothills of the Canadian Rockies. Detailed information on the dispersal of adult butterflies among this same set of populations was available. Simple and partial Mantel tests were used to examine the relationships between genetic distances, predicted rates of dispersal, and a number of landscape variables, all measured pairwise for 17 sample sites. Nei's standard genetic distance was negatively correlated with predicted dispersal. We observed a significant pattern of isolation by distance at a very small spatial scale. The distance between sites that was through forest was a stronger predictor of genetic distance than the distance through open meadow, indicating a significant effect of landscape on population genetic structure beyond that of simple isolation by distance. Our results suggest that rises in the tree-line in alpine areas, caused by global warming, will lead to reduced gene flow among populations of P. smintheus.     

A meta‐analysis of isolation by distance: relic or reference standard for landscape genetics?

Isolation by distance (IBD) has been a common measure of genetic structure among populations and is based on Euclidean distances among populations. Whereas IBD does not incorporate geographic complexity (e.g. dispersal barriers, corridors) that may better predict genetic structure, a new approach (landscape genetics) joins landscape ecology with population genetics to better model genetic structure. Should IBD be set aside or should it persist as the most simple model in landscape genetics? We evaluated the status of IBD by collecting and analyzing results of 240 IBD data sets among diverse taxa and study systems. IBD typically represented a low proportion of variance in genetic structure (mean r²=0.22) in part because many studies included relatively few populations (mean=11). The number of populations studied (N) was asymptotically related to IBD significance; a study with 9 populations has only 50% probability of significance, while one with >23 populations will have 90% probability of significance. Surprisingly, ectothermic animals were significantly (p=0.0018) more likely to have significant IBD than endotherms, which suggests a metabolic basis underlying gene flow rates. We also observed marginally significant effects on IBD significance for a) taxa in general and b) dispersal modes within actively‐dispersing endotherms. Other factors analyzed (genetic markers, genetic distances, habitats, active or passive dispersal, plant growth form) did not significantly affect IBD, likely related to typical N. For multiple reasons we conclude that IBD should continue as the simplest reference standard against which all other, more complex models should be compared in landscape genetics research.     

Landscape genetic analyses reveal cryptic population structure and putative selection gradients in a large-scale estuarine environment

Disentangling the relative contributions of selective and neutral processes underlying phenotypic and genetic variation under natural, environmental conditions remains a central challenge in evolutionary ecology. However, much of the variation that could be informative in this area of research is likely to be cryptic in nature; thus, the identification of wild populations suitable for study may be problematic. We use a landscape genetics approach to identify such populations of three-spined stickleback inhabiting the Saint Lawrence River estuary. We sampled 1865 adult fish over multiple years. Individuals were genotyped for nine microsatellite loci, and georeferenced multilocus data were used to infer population groupings, as well as locations of genetic discontinuities, under a Bayesian model framework ( geneland ). We modelled environmental data using nonparametric multiple regression to explain genetic differentiation as a function of spatio-ecological effects. Additionally, we used genotype data to estimate dispersal and gene flow to parameterize a simple model predicting adaptive vs. plastic divergence between demes. We demonstrate a bipartite division of the genetic landscape into freshwater and maritime zones, independent of geographical distance. Moreover, we show that the greatest proportion of genetic variation (31.5%) is explained by environmental differences. However, the potential for either adaptive or plastic divergence between demes is highly dependent upon the strength of migration and selection. Consequently, we highlight the utility of landscape genetics as a tool for hypothesis generation and experimental design, to identify focal populations and putative selection gradients, in order to distinguish between phenotypic plasticity and local adaptation.     

Population structure of Columbia spotted frogs (Rana luteiventris) is strongly affected by the landscape

Landscape features such as mountains, rivers, and ecological gradients may strongly affect patterns of dispersal and gene flow among populations and thereby shape population dynamics and evolutionary trajectories. The landscape may have a particularly strong effect on patterns of dispersal and gene flow in amphibians because amphibians are thought to have poor dispersal abilities. We examined genetic variation at six microsatellite loci in Columbia spotted frogs (Rana luteiventris) from 28 breeding ponds in western Montana and Idaho, USA, in order to investigate the effects of landscape structure on patterns of gene flow. We were particularly interested in addressing three questions: (i) do ridges act as barriers to gene flow? (ii) is gene flow restricted between low and high elevation ponds? (iii) does a pond equal a ‘randomly mating population’ (a deme)? We found that mountain ridges and elevational differences were associated with increased genetic differentiation among sites, suggesting that gene flow is restricted by ridges and elevation in this species. We also found that populations of Columbia spotted frogs generally include more than a single pond except for very isolated ponds. There was also evidence for surprisingly high levels of gene flow among low elevation sites separated by large distances. Moreover, genetic variation within populations was strongly negatively correlated with elevation, suggesting effective population sizes are much smaller at high elevation than at low elevation. Our results show that landscape features have a profound effect on patterns of genetic variation in Columbia spotted frogs.     

Dispersal and genetic structure in the American marten, Martes americana

Natal dispersal in a vagile carnivore, the American marten (Martes americana), was studied by comparing radio-tracking data and microsatellite genetic structure in two populations occupying contrasting habitats. The genetic differentiation determined among groups of individuals using F(ST) indices appeared to be weak in both landscapes, and showed no increase with geographical distance. Genetic structure investigated using pairwise genetic distances between individuals conversely showed a pattern of isolation by distance (IBD), but only in the population occurring in a homogeneous high-quality habitat, therefore showing the advantage of individual-based analyses in detecting within-population processes and local landscape effects. The telemetry study of juveniles revealed a leptokurtic distribution of dispersal distances in both populations, and estimates of the mean squared parent-offspring axial distance (sigma2) inferred both from the genetic pattern of IBD and from the radio-tracking survey showed that most juveniles make little contribution to gene flow.     

Scale-dependence and mechanisms of dispersal in freshwater zooplankton

Communities of organisms form as a result of both interspecific and abiotic interactions within local habitat patches and dispersal among patches in a region. Local processes are expected to play a dominant role when dispersal occurs much more often than extinction. We performed two field experiments to examine rates and mechanisms of dispersal in freshwater pond zooplankton communities. First, we tested the effect of distance from a source on the rate of colonization of artificial habitat by placing wading pools at 5, 10, 30 and 60 m from two natural fishless ponds and observing the succession of zooplankton. Seventy-eight percent of the species in the source ponds that were capable of living in the pools colonized at least once during the experiment. A new species was found in the pools on average once every four days, suggesting that colonization events occur on the order of days to weeks for many species. Colonization rates declined further from the source at one pond but not the other, and the effect of distance was relatively weak at both ponds. This suggests that many species disperse broadly over short distances. The second experiment tested the role of animal vectors for zooplankton dispersal by restricting access to the pools. Eight treatments were imposed that excluded potential animal vectors along a body size gradient from large mammals to small insects. While the treatments affected zooplankton colonization, many species invaded even when all animals larger than 1 mm were excluded. Animal vectors may therefore be less important for dispersal than wind. Our results suggest that zooplankton are highly effective dispersers over short distances, and can disperse via several mechanisms. Local interactions should therefore play a dominant role in structuring these communities at small regional scales.     

Spatial genetic structure in a metapopulation of the land snail Cepaea nemoralis (Gastropoda: Helicidae)

Habitat fragmentation is a major force affecting demography and genetic structure of wild populations, especially in agricultural landscapes. The land snail Cepaea nemoralis (L.) was selected to investigate the impact of habitat fragmentation on the spatial genetic structure of an organism with limited dispersal ability. Genetic and morphological patterns were investigated at a local scale of a 500 m transect and a mesoscale of 4 x 4 km in a fragmented agricultural landscape while accounting for variation in the landscape using least-cost models. Analysis of microsatellite loci using expected heterozygosity (HE), pairwise genetic distance (FST/1-FST) and spatial autocorrelograms (Moran's I) as well as shell characteristics revealed spatial structuring at both scales and provided evidence for a metapopulation structure. Genetic diversity was related to morphological diversity regardless of landscape properties. This pointed to bottlenecks caused by founder effects after (re)colonization. Our study suggests that metapopulation structure depended on both landscape features and the shape of the dispersal function. A range of genetic spatial autocorrelation up to 80 m at the local scale and up to 800 m at the mesoscale indicated leptokurtic dispersal patterns. The metapopulation dynamics of C. nemoralis resulted in a patchwork of interconnected, spatially structured subpopulations. They were shaped by gene flow which was affected by landscape features, the dispersal function and an increasing role of genetic drift with distance.     

Zooplankton community structure and environmental conditions in a set of interconnected ponds

We studied the zooplankton community structure in a set of 33 interconnected shallow ponds that are restricted to a relatively small area (`De Maten', Genk, Belgium, 200 ha). As the ponds share the same water source, geology and history, and as the ponds are interconnected (reducing chance effects of dispersal with colonisation), differences in zooplankton community structure can be attributed to local biotic and abiotic interactions. We studied zooplankton community, biotic (phytoplankton, macrophyte cover, fish densities, macroinvertebrate densities), abiotic (turbidity, nutrient concentrations, pH, conductivity, iron concentration) and morphometric (depth, area, perimeter) characteristics of the different ponds. Our results indicate that the ponds differ substantially in their zooplankton community structure, and that these differences are strongly related to differences in trophic structure and biotic interactions, in concordance with the theory of alternative equilibria. Ponds in the clear-water state are characterised by large Daphnia species and species associated with the littoral zone, low chlorophyll-a concentrations, low fish densities and high macroinvertebrate densities. Ponds in the turbid-water state are characterised by high abundances of rotifers, cyclopoid copepods and the opposite environmental conditions. Some ponds show an intermediate pattern, with a dominance of small Daphnia species. Our results show that interconnected ponds may differ strongly in zooplankton community composition, and that these differences are related to differences in predation intensity (top-down) and habitat diversity (macrophyte cover).     

Inference of hazel grouse population structure using multilocus data: a landscape genetic approach

In conservation and management of species it is important to make inferences about gene flow, dispersal and population structure. In this study, we used 613 georeferenced tissue samples from hazel grouse (Bonasa bonasia) where each individual was genotyped at 12 microsatellite loci to make inference on population genetic structure, gene flow and dispersal in northern Sweden. Observed levels of genetic diversity suggest that Swedish hazel grouse do not suffer loss of genetic diversity compared with other grouse species. We found significant F(IS) (deviation from Hardy-Weinberg expectations) over the entire sample using jack-knifed estimators over loci, which is most likely explained by a Wahlund effect. With the use of spatial autocorrelation methods, we detected significant isolation by distance among individuals. Neighbourhood size was estimated in the order of 62-158 individuals corresponding to a dispersal distance of 950-1500 m. Using a spatial statistical model for landscape genetics to infer the number of populations and the spatial location of genetic discontinuities between these populations we found indications that Swedish hazel grouse are divided into a northern and a southern population. We could not find a sharp border between these two populations and none of the observed borders appeared to coincide with any potential geographical barriers.These results imply that gene flow appears somewhat unrestricted in the boreal taiga forests of northern Sweden and that the two populations of hazel grouse in Sweden may be explained by the post-glacial reinvasion history of the Scandinavian Peninsula.     

Landscape genetics of a specialized grasshopper inhabiting highly fragmented habitats: a role for spatial scale

Aim The study of geographical discontinuities in the distribution of genetic variability in natural populations is a central topic in both evolutionary and conservation research. In this study, we aimed to analyse (1) the factors associated with genetic diversity at the landscape spatial scale in the highly specialized grasshopper Mioscirtus wagneri and (2) to identify the relative contribution of alternative factors to the observed patterns of genetic structure in this species. Location La Mancha region, Central Spain. Methods We sampled 28 populations of the grasshopper M. wagneri and genotyped 648 individuals at seven microsatellite loci. We employed a causal modelling approach to identify the most influential variables associated with genetic differentiation within a multiple hypothesis‐testing framework. Results We found that genetic diversity differs among populations located in different river basins and decreases with population isolation. Causal modelling analyses showed variability in the relative influence of the studied landscape features across different spatial scales. When a highly isolated population is considered, the analyses suggested that geographical distance is the only factor explaining the genetic differentiation between populations. When that population is excluded, the causal modelling analysis revealed that elevation and river basins are also relevant factors contributing to explaining genetic differentiation between the studied populations. Main conclusions These results indicate that the spatial scale considered and the inclusion of outlier populations may have important consequences on the inferred contribution of alternative landscape factors on the patterns of genetic differentiation even when all populations are expected to similarly respond to landscape structure. Thus, a multiscale perspective should also be incorporated into the landscape genetics framework to avoid biased conclusions derived from the spatial scale analysed and/or the geographical distribution of the studied populations.     

Streams over mountains: influence of riparian connectivity on gene flow in the Pacific jumping mouse (Zapus trinotatus)

In species affiliated with heterogeneous habitat, we expect gene flow to be restricted due to constraints placed on individual movement by habitat boundaries. This is likely to impact both individual dispersal and connectivity between populations. In this study, a GIS-based landscape genetics approach was used, in combination with fine-scale spatial autocorrelation analysis and the estimation of recent intersubpopulation migration rates, to infer patterns of dispersal and migration in the riparian-affiliated Pacific jumping mouse (Zapus trinotatus). A total of 228 individuals were sampled from nine subpopulations across a system of three rivers and genotyped at eight microsatellite loci. Significant spatial autocorrelation among individuals revealed a pattern of fine-scale spatial genetic structure indicative of limited dispersal. Geographical distances between pairwise subpopulations were defined following four criteria: (i) Euclidean distance, and three landscape-specific distances, (ii) river distance (distance travelled along the river only), (iii) overland distance (similar to Euclidean, but includes elevation), and (iv) habitat-path distance (a least-cost path distance that models movement along habitat pathways). Pairwise Mantel tests were used to test for a correlation between genetic distance and each of the geographical distances. Significant correlations were found between genetic distance and both the overland and habitat-path distances; however, the correlation with habitat-path distance was stronger. Lastly, estimates of recent migration rates revealed that migration occurs not only within drainages but also across large topographic barriers. These results suggest that patterns of dispersal and migration in Pacific jumping mice are largely determined by habitat connectivity.     

Drivers of plant species’ potential to spread: the importance of demography versus seed dispersal

Understanding the ability of plants to spread is important for assessing conservation strategies, landscape dynamics, invasiveness and ability to cope with climate change. While long‐distance seed dispersal is often viewed as a key process in population spread, the importance of inter‐specific variation in demography is less explored. Indeed, the relative importance of demography vs seed dispersal in determining population spread is still little understood. We modelled species’ potential for population spread in terms of annual migration rates for a set of species inhabiting dry grasslands of central Europe. Simultaneously, we estimated the importance of demographic (population growth rate) versus long‐distance dispersal (99th percentile dispersal distance) characteristics for among‐species differences in modelled population spread. In addition, we assessed how well simple proxy measures related to demography (the number and survival of seedlings, the survival of flowering individuals) and dispersal (plant height, terminal velocity and wind speed during dispersal) predicted modelled spread rates. We found that species’ demographic rates were the more powerful predictors of species’ modelled potential to spread than dispersal. Furthermore, our simple proxies were correlated with modelled species spread rates and together their predictive power was high. Our findings highlight that for understanding variation among species in their potential for population spread, detailed information on local demography and dispersal might not always be necessary. Simple proxies or assumptions that are based primarily on species demography could be sufficient.     

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.     

Landscape genetics of an endangered lemur (Propithecus tattersalli) within its entire fragmented range

Habitat fragmentation may strongly reduce individuals’ dispersal among resource patches and hence influence population distribution and persistence. We studied the impact of landscape heterogeneity on the dispersal of the golden‐crowned sifaka (Propithecus tattersalli), an endangered social lemur species living in a restricted and highly fragmented landscape. We combined spatial analysis and population genetics methods to describe population units and identify the environmental factors which best predict the rates and patterns of genetic differentiation within and between populations. We used non‐invasive methods to genotype 230 individuals at 13 microsatellites in all the main forest fragments of its entire distribution area. Our analyses suggest that the Manankolana River and geographical distance are the primary structuring factors, while a national road crossing the region does not seem to impede gene flow. Altogether, our results are in agreement with a limited influence of forest habitat connectivity on gene flow patterns (except for North of the species’ range), suggesting that dispersal is still possible today among most forest patches for this species. Within forest patches, we find that dispersal is mainly among neighbouring social groups, hence confirming previous behavioural observations.     

Fine-scale spatial genetic structure and dispersal among spotted salamander (Ambystoma maculatum) breeding populations

We examined fine-scale genetic variation among breeding aggregations of the spotted salamander (Ambystoma maculatum) to quantify dispersal, interpopulation connectivity and population genetic structure. Spotted salamanders rely on temporary ponds or wetlands for aggregate breeding. Adequate breeding sites are relatively isolated from one another and field studies suggest considerable adult site fidelity; therefore, we expected to find population structure and differentiation at small spatial scales. We used microsatellites to estimate population structure and dispersal among 29 breeding aggregations in Tompkins County, New York, USA, an area encompassing 1272 km(2). Bayesian and frequency-based analyses revealed fine-scale genetic structure with two genetically defined demes: the North deme included seven breeding ponds, and the South deme included 13 ponds. Nine ponds showed evidence of admixture between these two genetic pools. Bayesian assignment tests for detection of interpopulation dispersal indicate that immigration among ponds is common within demes, and that certain populations serve as sources of immigrants to neighbouring ponds. Likewise, spatial genetic correlation analyses showed that populations < or = 4.8 km distant from each other show significant genetic correlation that is not evident at higher scales. Within-population levels of relatedness are consistently larger than expected if mating were completely random across ponds, and in the case of a few ponds, within-population processes such as inbreeding or reproductive skew contribute significantly to differentiation from neighbouring ponds. Our data underscore the importance of these within-population processes as a source of genetic diversity across the landscape, despite considerable population connectivity. Our data further suggest that spotted salamander breeding groups behave as metapopulations, with population clusters as functional units, but sufficient migration among demes to allow for potential rescue and recolonization. Amphibian habitats are becoming increasingly fragmented and a clear understanding of dispersal and patterns of population connectivity for taxa with different ecologies and life histories is crucial for their conservation.     

Directed dispersal by rotational shepherding supports landscape genetic connectivity in a calcareous grassland plant

Directed dispersal by animal vectors has been found to have large effects on the structure and dynamics of plant populations adapted to frugivory. Yet, empirical data are lacking on the potential of directed dispersal by rotational grazing of domestic animals to mediate gene flow across the landscape. Here, we investigated the potential effect of large‐flock shepherding on landscape‐scale genetic structure in the calcareous grassland plant Dianthus carthusianorum, whose seeds lack morphological adaptations to dispersal to animals or wind. We found a significant pattern of genetic structure differentiating population within grazed patches of three nonoverlapping shepherding systems and populations of ungrazed patches. Among ungrazed patches, we found a strong and significant effect of isolation by distance (r = 0.56). In contrast, genetic distance between grazed patches within the same herding system was unrelated to geographical distance but significantly related to distance along shepherding routes (r = 0.44). This latter effect of connectivity along shepherding routes suggests that gene flow is spatially restricted occurring mostly between adjacent populations. While this study used nuclear markers that integrate gene flow by pollen and seed, the significant difference in the genetic structure between ungrazed patches and patches connected by large‐flock shepherding indicates the potential of directed seed dispersal by sheep across the landscape.