Population genetics of Mioscirtus wagneri,a grasshopper showing a highly fragmented distribution

The genetic consequences of population fragmentation and isolation are major issues in conservation biology. In this study we analyse the genetic variability and structure of the Iberian populations of Mioscirtus wagneri, a specialized grasshopper exclusively inhabiting highly fragmented hypersaline low grounds. For this purpose we have used seven species‐specific microsatellite markers to type 478 individuals from 24 localities and obtain accurate estimates of their genetic variability. Genetic diversity was relatively low and we detected genetic signatures suggesting that certain populations of M. wagneri have probably passed through severe demographic bottlenecks. We have found that the populations of this grasshopper show a strong genetic structure even at small geographical scales, indicating that they mostly behave as isolated populations with low levels of gene flow among them. Thus, several populations can be regarded as independent and genetically differentiated units which require adequate conservation strategies to avoid eventual extinctions that in highly isolated localities are not likely to be compensated for with the arrival of immigrants from neighbouring populations. Overall, our results show that these populations probably represent the ‘fragments’ of a formerly more widespread population and highlight the importance of protecting Iberian hypersaline environments due to the high number of rare and endangered species they sustain.     

Testing the role of ancient and contemporary landscapes on structuring genetic variation in a specialist grasshopper

Understanding the processes underlying spatial patterns of genetic diversity and structure of natural populations is a central topic in evolutionary biogeography. In this study, we combine data on ancient and contemporary landscape composition to get a comprehensive view of the factors shaping genetic variation across the populations of the scrub‐legume grasshopper (Chorthippus binotatus binotatus) from the biogeographically complex region of southeast Iberia. First, we examined geographical patterns of genetic structure and employed an approximate Bayesian computation (ABC) approach to compare different plausible scenarios of population divergence. Second, we used a landscape genetic framework to test for the effects of (1) Late Miocene paleogeography, (2) Pleistocene climate fluctuations, and (3) contemporary topographic complexity on the spatial patterns of population genetic differentiation. Genetic structure and ABC analyses supported the presence of three genetic clusters and a sequential west‐to‐east splitting model that predated the last glacial maximum (LGM, c. 21 Kya). Landscape genetic analyses revealed that population genetic differentiation was primarily shaped by contemporary topographic complexity, but was not explained by any paleogeographic scenario or resistance distances based on climate suitability in the present or during the LGM. Overall, this study emphasizes the need of integrating information on ancient and contemporary landscape composition to get a comprehensive view of their relative importance to explain spatial patterns of genetic variation in organisms inhabiting regions with complex biogeographical histories.     

Population genomic analysis suggests strong influence of river network on spatial distribution of genetic variation in invasive saltcedar across the southwestern United States

Understanding the complex influences of landscape and anthropogenic elements that shape the population genetic structure of invasive species provides insight into patterns of colonization and spread. The application of landscape genomics techniques to these questions may offer detailed, previously undocumented insights into factors influencing species invasions. We investigated the spatial pattern of genetic variation and the influences of landscape factors on population similarity in an invasive riparian shrub, saltcedar (Tamarix L.) by analysing 1,997 genomewide SNP markers for 259 individuals from 25 populations collected throughout the southwestern United States. Our results revealed a broad‐scale spatial genetic differentiation of saltcedar populations between the Colorado and Rio Grande river basins and identified potential barriers to population similarity along both river systems. River pathways most strongly contributed to population similarity. In contrast, low temperature and dams likely served as barriers to population similarity. We hypothesize that large‐scale geographic patterns in genetic diversity resulted from a combination of early introductions from distinct populations, the subsequent influence of natural selection, dispersal barriers and founder effects during range expansion.     

Spatial scale affects landscape genetic analysis of a wetland grasshopper

Most landscape genetic studies assess the impact of landscape elements on species' dispersal and gene flow. Many of these studies perform their analysis on all possible population pairs in a study area and do not explicitly consider the effects of spatial scale and population network topology on their results. Here, we examined the effects of spatial scale and population network topology on the outcome of a landscape genetic analysis. Additionally, we tested whether the relevant spatial scale of landscape genetic analysis could be defined by population network topology or by isolation‐by‐distance (IBD) patterns. A data set of the wetland grasshopper Stethophyma grossum, collected in a fragmented agricultural landscape, was used to analyse population network topology, IBD patterns and dispersal habitats, using least‐cost transect analysis. Landscape genetic analyses neglecting spatial scale and population network topology resulted in models with low fits, with which a most likely dispersal habitat could not be identified. In contrast, analyses considering spatial scale and population network topology resulted in high model fits by restricting landscape genetic analysis to smaller scales (0–3 km) and neighbouring populations, as represented by a Gabriel graph. These models also successfully identified a likely dispersal habitat of S. grossum. The above results suggest that spatial scale and potentially population network topology should be more explicitly considered in future landscape genetic analyses.     

Fine-scale spatial genetic structure and within population male-biased gene-flow in the grasshopper <Emphasis Type="Italic">Mioscirtus wagneri</Emphasis>

Dispersal is a life history trait that plays a key role in population dynamics, determining gene flow and influencing the size, structure and persistence of populations. For these reasons, the study of the genetic consequences of dispersal can be considered a central topic in both conservation and population genetics. In this study we examine the patterns of fine-scale genetic structure within two populations of the grasshopper Mioscirtus wagneri (Orhoptera: Acrididae). For this purpose, we have used seven species-specific microsatellite markers to type 266 individuals from two populations (Peña Hueca and El Salobral) located in Central Spain. We have found subtle genetic differentiation between some sampling patches and significant kinship structures up to 25 m distance which were particularly patent for females. In Peña Hueca locality, patterns of isolation-by-distance at both the patch scale and the individual level have also revealed an association between genetic differentiation/similarity and geographical distance in females but not in males. Overall, these data suggest a fine-scale spatial genetic substructure in the studied populations which seems to be mainly driven by female philopatry. Such pattern of within population genetic structure together with the inferred restricted dispersal distances is likely to contribute to reduce effective population sizes and inter-population gene flow. This can erode genetic variability and limit the colonization ability of this orthoptera, factors which can ultimately compromise the long-term persistence of their small size and isolated populations.     

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

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

Dendritic connectivity shapes spatial patterns of genetic diversity: a simulation‐based study

Landscape features notoriously affect spatial patterns of biodiversity. For instance, in dendritic ecological networks (such as river basins), dendritic connectivity has been proposed to create unique spatial patterns of biodiversity. Here, we compared genetic datasets simulated under a lattice‐like, a dendritic and a circular landscape to test the influence of dendritic connectivity on neutral genetic diversity. The circular landscape had a level of connectivity similar to that of the dendritic landscape, so as to isolate the influence of dendricity on genetic diversity. We found that genetic diversity and differentiation varied strikingly among the three landscapes. For instance, the dendritic landscape generated higher total number of alleles and higher global F_st than the lattice‐like landscape, and these indices also varied between the dendritic and the circular landscapes, suggesting an effect of dendricity. Furthermore, in the dendritic landscape, allelic richness was higher in highly connected demes (e.g. confluences in rivers) than in low‐connected demes (e.g. upstream and downstream populations), which was not the case in the circular landscape, hence confirming the major role of dendricity. This led to bell‐shaped distributions of allelic richness along an upstream–downstream gradient. Conversely, genetic differentiation (F_st) was lower in highly than in low‐connected demes (which was not observed in circular landscape), and significant patterns of isolation by distance (IBD) were also observed in the dendritic landscape. We conclude that in dendritic networks, the combined influence of dendricity and connectivity generates unique spatial patterns of neutral genetic diversity, which has implications for population geneticists and conservationists.     

Geographical and Ecological Drivers of Mitonuclear Genetic Divergence in a Mediterranean Grasshopper

The study of the neutral and/or selective processes driving genetic variation in natural populations is central to determine the evolutionary history of species and lineages and understand how they interact with different historical and contemporary components of landscape heterogeneity. Here, we combine nuclear and mitochondrial data to study the processes shaping genetic divergence in the Mediterranean esparto grasshopper (Ramburiella hispanica). Our analyses revealed the presence of three main lineages, two in Europe that split in the Early-Middle Pleistocene and one in North Africa that diverged from the two European ones after the Messinian. Lineage-specific potential distribution models and tests of environmental niche differentiation suggest that the phylogeographic structure of the species was driven by allopatric divergence due to the re-opening of the Gibraltar strait at the end of the Messinian (Europe–Africa split) and population fragmentation in geographically isolated Pleistocene climatic refugia (European split). Although we found no evidence for environment as an important driver of genetic divergence at the onset of lineage formation, our analyses considering the spatial distribution of populations and different aspects of landscape composition suggest that genetic differentiation at mitochondrial loci was largely explained by environmental dissimilarity, whereas resistance-based estimates of geographical distance were the only predictors of genetic differentiation at nuclear markers. Overall, our study shows that although historical factors have largely shaped concordant range-wide patterns of mitonuclear genetic structure in the esparto grasshopper, different contemporary processes (neutral gene flow vs. environmental-based selection) seem to be governing the spatial distribution of genetic variation in the two genomes.     

Spatiotemporally explicit demographic modelling supports a joint effect of historical barriers to dispersal and contemporary landscape composition on structuring genomic variation in a red‐listed grasshopper

Inferring the processes underlying spatial patterns of genomic variation is fundamental to understand how organisms interact with landscape heterogeneity and to identify the factors determining species distributional shifts. Here, we use genomic data (restriction site‐associated DNA sequencing) to test biologically informed models representing historical and contemporary demographic scenarios of population connectivity for the Iberian cross‐backed grasshopper Dociostaurus hispanicus, a species with a narrow distribution that currently forms highly fragmented populations. All models incorporated biological aspects of the focal taxon that could hypothetically impact its geographical patterns of genomic variation, including (a) spatial configuration of impassable barriers to dispersal defined by topographic landscapes not occupied by the species; (b) distributional shifts resulting from the interaction between the species bioclimatic envelope and Pleistocene glacial cycles; and (c) contemporary distribution of suitable habitats after extensive land clearing for agriculture. Spatiotemporally explicit simulations under different scenarios considering these aspects and statistical evaluation of competing models within an Approximate Bayesian Computation framework supported spatial configuration of topographic barriers to dispersal and human‐driven habitat fragmentation as the main factors explaining the geographical distribution of genomic variation in the species, with no apparent impact of hypothetical distributional shifts linked to Pleistocene climatic oscillations. Collectively, this study supports that both historical (i.e., topographic barriers) and contemporary (i.e., anthropogenic habitat fragmentation) aspects of landscape composition have shaped major axes of genomic variation in the studied species and emphasizes the potential of model‐based approaches to gain insights into the temporal scale at which different processes impact the demography of natural populations.     

Population structure and landscape genetics in the endangered subterranean rodent Ctenomys porteousi

In order to devise adequate conservation and management strategies for endangered species, it is important to incorporate a reliable understanding of its spatial population structure, detecting the existence of demographic partitions throughout its geographical range and characterizing the distribution of its genetic diversity. Moreover, in species that occupy fragmented habitats it is essential to know how landscape characteristics may affect the genetic connectivity among populations. In this study we use eight microsatellite markers to analyze population structure and gene flow patterns in the complete geographic range of the endangered rodent Ctenomys porteousi. Also, we use landscape genetics approaches to evaluate the effects of landscape configuration on the genetic connectivity among populations. In spite of geographical proximity of the sampling sites (8–27 km between the nearest sites) and the absence of marked barriers to individual movement, strong population structure and low values of gene flow were observed. Genetic differentiation among sampling sites was consistent with a simple model of isolation by distance, where peripheral areas showed higher population differentiation than those sites located in the central area of the species’ distribution. Landscape genetics analysis suggested that habitat fragmentation at regional level has affected the distribution of genetic variation among populations. The distance of sampling sites to areas of the landscape having higher habitat connectivity was the environmental factor most strongly related to population genetic structure. In general, our results indicate strong genetic structure in C. porteousi, even at a small spatial scale, and suggest that habitat fragmentation could increase the population differentiation.     

Genetic structure of brown trout, salmo trutta l., at the southern limit of the distribution range of the anadromous form

Genetic variation at 33 protein loci was investigated in 41 wild brown trout populations from four river basins in Galicia (northwest Spain) to analyse the amount and distribution of genetic diversity in a marginal area, located in the distribution limit of the anadromous form of this species. The genetic diversity detected within populations (H between 0 and 6%) lies within the range quoted for this species in previous reports. The Mino, the most southern river basin analysed, showed a significantly lower genetic diversity and the highest genetic differentiation among the river basins studied. The hierarchical gene diversity analysis showed high population differentiation in a restricted area (GST = 27%), mostly due to differences among populations within basins (GSC = 22%). The reduction of GST observed when the isolated samples were excluded from the analysis (GST = 17%) showed the importance of habitat fragmentation on the heterogeneity detected. Gene flow among populations was comparatively evaluated by three indirect methods, which in general revealed low figures of absolute number of migrants per generation, slightly higher than 1. The gene flow among basins reflected a positive relationship with geographical distance. This trend was confirmed by the significant correlation observed between geographical and genetic distances, including all population pairs, which suggests a component of isolation by distance in brown trout genetic structure. Nevertheless, the nonsignificant intrabasin correlation demonstrates the complexity of genetic relationships among populations in this species. The model of genetic structure in brown trout is discussed in the light of the results obtained.     

The interplay of dispersal limitation,rivers, and historical events shapes the genetic structure of an Amazonian frog

Disentangling the impact of landscape features such as rivers and historical events on dispersal is a challenging but necessary task to gain a comprehensive picture of the evolution of diverse biota such as that found in Amazonia. Adenomera andreae, a small, territorial, terrestrial frog species of the Amazonian forest represents a good model for such studies. We combined cytochrome b sequences with 12 microsatellites to investigate the genetic structure at two contrasted spatial scales in French Guiana: along a ～6‐km transect, to evaluate dispersal ability, and between paired bank populations along a ～65‐km stretch of the Approuague river, to test the effect of rivers as barriers to dispersal. We observed significant spatial genetic structure between individuals at a remarkably small geographical scale, and conclude that the species has a restricted dispersal ability that is probably tied to its life‐history traits. Mitochondrial and microsatellite data also indicate a high level of differentiation among populations on opposite banks of the river, and, in some cases, among populations on the same riverbank. These results suggest that the observed population structure in A. andreae is the result of restricted dispersal abilities combined with the action of rivers and Quaternary population isolation. Given that Amazonia hosts a great portion of anurans, as well as other small vertebrates, that display life‐history traits comparable with A. andreae, we argue that our analyses provide new insights into the complex interactions among evolutionary processes shaping Amazonian biodiversity. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106 , 356–373.     

Population genetic structure of two ecologically distinct multimammate rats: the commensal Mastomys natalensis and the wild Mastomys erythroleucus in southeastern Senegal

Using the same set of microsatellite markers, we compared the population genetic structure of two Mastomys species, one being exclusively commensal in southeastern Senegal, and the other being continuously distributed outside villages in this region. Both species were sampled in the same landscape context and at the same spatial scale. According to the expectations based on the degree of habitat patchiness (which is higher for commensal populations in this rural area), genetic diversity was lower and genetic differentiation was higher in commensal populations of Mastomys natalensis than in wild populations of Mastomys erythroleucus. Contrasting estimates of effective dispersal and current migration rates corroborates previous data on differences in social structure between the two species. Isolation-by-distance analyses showed that human-mediated dispersal is not a major factor explaining the pattern of genetic differentiation for M. natalensis, and that gene flow is high and random between M. erythroleucus populations at the spatial scale considered.     

Genetic and morphological divergence at different spatiotemporal scales in the grasshopper Mioscirtus wagneri (Orthoptera: Acrididae)

The study of the association between morphological and genetic divergence can provide important information on the factors determining population differentiation and gene flow at different spatiotemporal scales. In this study we analyze the congruence between morphological and genetic divergence in the Iberian populations of Mioscirtus wagneri, a specialized grasshopper exclusively inhabiting highly fragmented hypersaline low grounds. We have found strong morphological variation among the studied localities and among mtDNA- and microsatellite-based genetic clusters. However, we have detected some cases of morphological convergence between highly differentiated populations. By contrast, certain genetically homogeneous populations at both mtDNA and microsatellite markers showed significant morphological differentiation which may be explained by phenotypic plasticity or divergent selection pressures acting at different spatiotemporal scales. Mantel tests also revealed that morphological divergence was associated with microsatellite- but not with mtDNA-based genetic distances. Overall, this study suggests that morphological traits can provide additional information on the underlying population genetic structure when only data on scarcely variable mtDNA markers is available. Thus, morphology can retain useful information on genetic structure and has the benefit over molecular methods of being inexpensive, offering a preliminary/complementary useful criterion for the establishment of management units necessary to guide conservation policies.     

Landscape genetics and spatial pattern of phenotypic variation of Eristalis tenax across Europe

A study of population connectivity of the migratory insect species, such as dronefly Eristalis tenax (Diptera, Syrphidae), has an essential importance in understanding the relative influence of the evolutionary forces and environmental features that interact in the spatial distribution of molecular and morphological diversity. However, specific study aiming to understand spatial genetic structure of dronefly populations and its migratory potential is lacking. Hence, we studied a spatial pattern of genetic and phenotypic variation of seven European populations of E. tenax incorporating landscape genetic methods using allozyme data, wing size and shape and abdominal colour pattern. Based on the observed lack of genotypic structuring, we suggested that there has been sufficient long‐distance gene flow to effectively homogenize population structuring at a broader geographical scale. Wing shape similarity among populations and an overlap of abdominal colour variation showed no clear clustering related to geography, which is in congruence with genetic data. However, genetic (F_ST values) and phenotypic (wing size) data and landscape genetics indicated subdivision between the Balkan populations (four Serbian samples) and populations from Central (Germany and Switzerland) and Northern (Finland) Europe. These findings indicated a potential connection between the Central and Northern Europe supporting the Central European origin of the flies caught in Finland. Thus, by performing spatial analysis and combining genetic–morphological approach, we shed light on the movement pattern in complex landscapes and thus provided the necessary guidelines to a broad‐scale analysis of this widespread generalist pollinator.     

Linking genetic and ecological differentiation in an ungulate with a circumpolar distribution

Genetic differentiation among populations may arise from the disruption of gene flow due to local adaptation to distinct environments and/or neutral accumulation of mutations and genetic drift resulted from geographical isolation. Quantifying the role of these processes in determining the genetic structure of natural populations remains challenging. Here, we analyze the relative contribution of isolation‐by‐resistance (IBR), isolation‐by‐environment (IBE), genetic drift and historical isolation in allopatry during Pleistocene glacial cycles on shaping patterns of genetic differentiation in caribou/reindeer populations Rangifer tarandus across the entire distribution range of the species. Our study integrates analyses at range‐wide and regional scales to partial out the effects of historical and contemporary isolation mechanisms. At the circumpolar scale, our results indicate that genetic differentiation is predominantly explained by IBR and historical isolation. At a regional scale, we found that IBR, IBE and population size significantly explained the spatial distribution of genetic variation among populations belonging to the Euro‐Beringian lineage within North America. In contrast, genetic differentiation among populations within the North American lineage was predominantly explained by IBR and population size, but not IBE. We also found discrepancies between genetic and ecotype designation across the Holarctic species distribution range. Overall, these results indicate that multiple isolating mechanisms have played roles in shaping the spatial distribution of genetic variation across the distribution range of a large mammal with high potential for gene flow. Considering multiple spatial scales and simultaneously testing a comprehensive suite of potential isolating mechanisms, our study contributes to understand the ecological and evolutionary processes underlying organism–landscape interactions.     

A cryptic genetic boundary in remnant populations of a long‐lived,bird‐pollinated shrub Banksia sphaerocarpa var. caesia (Proteaceae)

Plant species often exhibit genetic structure at multiple spatial scales. Detection of this structure may depend on the sampling strategy used. We intensively sampled a common, naturally patchy Banksia species within a 200 km² region, in order to assess patterns of genetic diversity and structure at multiple spatial scales. In total, 1321 adult shrubs from 37 geographical populations were genotyped using eight highly polymorphic microsatellite markers developed for the species. Genetic structure was detected at three spatial scales. First, we identified a stark and unexpected division of the landscape into two genetic subregions, one to the north‐east and one to the south‐west of the sampling grid. This differentiation was based on sudden, highly structured differences in common allele frequencies, the cause of which is unknown but that may relate to physical and reproductive barriers to gene flow, localised selection, and/or historical processes. Second, we observed genetic differentiation of populations within these subregions, reflecting previously described patterns of restricted pollen flow in this species. Finally, fine‐scale genetic structure, although weak, was observed within some of the populations (mean S_P = 0.01837). When feasible, intensive sampling may uncover cryptic patterns of genetic structure that would otherwise be overlooked when sampling at broader spatial scales. Further studies using a similar sampling strategy may reveal this type of discontinuity to be a feature of other south‐western Australian taxa and has implications for our understanding of evolution in this landscape as well as conservation into the future. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 115 , 241–255.     

Drivers of genetic differentiation in a generalist insect‐pollinated herb across spatial scales

The isolation‐by‐distance model (IBD) predicts that genetic differentiation among populations increases with geographic distance. Yet, empirical studies show that a variety of ecological, topographic and historical factors may override the effect of geographic distance on genetic variation. This may particularly apply to species with narrow but highly heterogeneous distribution ranges, such as those occurring along elevational gradients. Using nine SSR markers, we study the genetic differentiation of the montane pollination‐generalist herb, Erysimum mediohispanicum. Because the effects of any given factor may depend on the geographic scale considered, we investigate the contribution of different environmental and historical factors at three different spatial scales. We evaluate five competing models that put forward the role of geographic distance, local environmental factors [biotic interactions (IBEb) and climatic variables (IBEa)], landscape resistance (IBR) and phylogeographic patterns (IBP), respectively. We find significant IBD regardless of the spatial scale and the genetic distance estimator considered. However, IBEa and IBP also play a prominent role in shaping genetic differentiation patterns at the larger spatial scales, and IBR is significant at the fine spatial scale. Overall, our results highlight the importance of combining different estimators, statistical approaches and spatial scales to disentangle the relative importance of the various ecological factors contributing to the shaping of genetic divergence patterns in natural populations.     

Influence of geography and environment on patterns of genetic differentiation in a widespread submerged macrophyte,Eurasian watermilfoil (Myriophyllum spicatum L., Haloragaceae)

The effects of geographic and environmental variables on the pattern of genetic differentiation have been thoroughly studied, whereas empirical studies on aquatic plants are rare. We examined the spatial genetic differentiation of 58 Myriophyllum spicatum populations distributed throughout China with 12 microsatellite loci, and we analyzed its association with geographic distance, geographic barriers, and environmental dissimilarity using causal modeling and multiple matrix regression with randomization (MMRR) analysis. Two genetic clusters were identified, and their geographic distribution suggested mountain ranges as a barrier to gene flow. The causal modeling revealed that both climate and geographic barriers significantly influenced genetic divergence among M. spicatum populations and that climate had the highest regression coefficient according to the MMRR analysis. This study showed that geography and environment together played roles in shaping the genetic structure of M. spicatum and that the influence of environment was greater. Our findings emphasized the potential importance of the environment in producing population genetic differentiation in aquatic plants at a large geographic scale.     

Landscape features impact connectivity between soil populations: a comparative study of gene flow in earthworms

Landscape features are known to alter the spatial genetic variation of aboveground organisms. Here, we tested the hypothesis that the genetic structure of belowground organisms also responds to landscape structure. Microsatellite markers were used to carry out a landscape genetic study of two endogeic earthworm species, Allolobophora chlorotica (N = 440, eight microsatellites) and Aporrectodea icterica (N = 519, seven microsatellites), in an agricultural landscape in the North of France, where landscape features were characterized with high accuracy. We found that habitat fragmentation impacted genetic variation of earthworm populations at the local scale. A significant relationship was observed between genetic diversity (H_e, A_r) and several landscape features in A. icterica populations and A. chlorotica. Moreover, a strong genetic differentiation between sites was observed in both species, with a low degree of genetic admixture and high F_st values. The landscape connectivity analysis at the regional scale, including isolation by distance, least‐cost path and cost‐weighted distance approaches, showed that genetic distances were linked to landscape connectivity in A. chlorotica. This indicates that the fragmentation of natural habitats has shaped their dispersal patterns and local effective population sizes. Landscape connectivity analysis confirmed that a priori favourable habitats such as grasslands may constitute dispersal corridors for these species.