Gene flow simulations demonstrate resistance of long-lived species to genetic erosion from habitat fragmentation

Habitat fragmentation restricts the movement of individuals across a landscape. In terrestrial and aquatic systems, barriers to movement can modify population and community dynamics at local or regional scales. This study contrasted life history traits related to lifespan with habitat fragmentation to determine impacts on species population genetic structure in the Neuse River Basin, USA. For this, we simulated gene flow among evenly-spaced populations in a river network and tracked individual and population genetics for 200 years. The modeled scenarios represent a full cross between five life history strategies and four riverscapes representing varying degrees of fragmentation. The five life history strategies include species (based on freshwater mussels) with average lifespans ranging from 10 to 50 years and age at maturity from 2 to 6 years. The movement landscapes included a (1) panmictic, (2) stepping-stone landscape allowing movement to only neighboring populations during each dispersal event, (3) partially-fragmented landscape divided by dams currently in the network, and (4) fully-fragmented landscape. Results suggest species with shorter lifespans have higher population genetic structure in fragmented landscapes than species with longer lifespans. Furthermore, species with shorter lifespans in highly fragmented landscapes may be harboring genetic degradation or decline as allele fixation and loss. Although anthropogenic fragmentation of many river systems is only 100–200 years old, the simulation indicates that species can respond genetically in that period of time. Additionally, the time frame of the simulation suggests that genetic impacts of habitat fragmentation in some species present in the Neuse River Basin may not yet be manifesting and restoration activities could be successful.     

Modelling single nucleotide effects in phosphoglucose isomerase on dispersal in the Glanville fritillary butterfly: coupling of ecological and evolutionary dynamics

Dispersal comprises a complex life-history syndrome that influences the demographic dynamics of especially those species that live in fragmented landscapes, the structure of which may in turn be expected to impose selection on dispersal. We have constructed an individual-based evolutionary sexual model of dispersal for species occurring as metapopulations in habitat patch networks. The model assumes correlated random walk dispersal with edge-mediated behaviour (habitat selection) and spatially correlated stochastic local dynamics. The model is parametrized with extensive data for the Glanville fritillary butterfly. Based on empirical results for a single nucleotide polymorphism (SNP) in the phosphoglucose isomerase (Pgi) gene, we assume that dispersal rate in the landscape matrix, fecundity and survival are affected by a locus with two alleles, A and C, individuals with the C allele being more mobile. The model was successfully tested with two independent empirical datasets on spatial variation in Pgi allele frequency. First, at the level of local populations, the frequency of the C allele is the highest in newly established isolated populations and the lowest in old isolated populations. Second, at the level of sub-networks with dissimilar numbers and connectivities of patches, the frequency of C increases with decreasing network size and hence with decreasing average metapopulation size. The frequency of C is the highest in landscapes where local extinction risk is high and where there are abundant opportunities to establish new populations. Our results indicate that the strength of the coupling of the ecological and evolutionary dynamics depends on the spatial scale and is asymmetric, demographic dynamics having a greater immediate impact on genetic dynamics than vice versa.     

Both population size and patch quality affect local extinctions and colonizations

Currently, the habitat of many species is fragmented, resulting in small local populations with individuals occasionally dispersing between the remaining habitat patches. In a solitary bee metapopulation, extinction probability was related to both local bee population sizes and pollen resources measured as host plant population size. Patch size, on the other hand, had no additional predictive power. The turnover rate of local bee populations in 63 habitat patches over 4 years was high, with 72 extinction events and 31 colonization events, but the pollen plant population was stable with no extinctions or colonizations. Both pollen resources and bee populations had strong and independent effects on extinction probability, but connectivity was not of importance. Colonizations occurred more frequently within larger host plant populations. For metapopulation survival of the bee, large pollen plant populations are essential, independent of current bee population size.     

Genetic (RAPD) diversity in Peromyscus maniculatus populations in a naturally fragmented landscape

We assessed the effects of long-term habitat fragmentation on genetic (random amplified polymorphic DNA) diversity in 11 Peromyscus maniculatus populations in the Lake Superior watershed. We analysed genetic structure at two spatial scales and the effect of island size and isolation on genetic diversity. At the regional scale, island populations differed from mainland populations (FST = 0.36), but mainland populations did not differ from each other (FST = 0.01). At the local scale, populations of the main island of Isle Royale differed from adjacent islet populations (P < 0.001; Monte Carlo approximation of Fisher's exact test), but not from each other (combined P = 0.63). Although geographical distance and genetic distance were positively correlated (P < 0.01; Mantel test), cluster analysis revealed some inconsistencies. Finally, genetic diversity was inversely related to isolation (P = 0.01), but had an unexpectedly negative relationship with island area (P = 0.03). The genetic structure of P. maniculatus populations in portions of the Lake Superior watershed appears to have been affected by long-term habitat fragmentation.     

Effects of geographic distance,sea barriers and habitat on the genetic structure and diversity of all-hybrid water frog populations

The history of population size and migration patterns leaves its mark in the genetics of populations. We investigate the genetic structure of the edible frog, Pelophylax esculentus in the Danish archipelago and adjacent countries. This frog is of particular interest because it is a hybrid that, in this area, forms all-hybrid populations of diploid (LR) and triploid (LLR and LRR) genomotypes with no (or very few) adults of the parental species (LL and RR). This study is the first to cover the entire geographic range of Danish, Swedish and German all-hybrid populations, documenting their extent and providing a broad picture of their diversity of neutral genetic markers and genomotype proportions. With 18 microsatellite markers, we found that genetic diversity declines northwards in agreement with the glacial refuge and central-marginal hypotheses; however, populations on small and medium-sized islands are no less diverse than those on large islands and continental peninsulas. Isolation by distance exists across the archipelago with limited influence of fragmentation by brackish seawater. The extremely low genetic diversity in all-hybrid populations, compared with adjacent populations, may be responsible for the maintenance of their special breeding system. We also show large variation among ponds in proportions of LLR, LR and LRR genomotypes, but little geographic pattern in their distribution. Instead, we found relationships between the genomotype proportions and some of 15 habitat parameters monitored. Body size differences among LLR, LR and LRR further suggest ecological differences.     

Habitat continuity and stepping‐stone oceanographic distances explain population genetic connectivity of the brown alga Cystoseira amentacea

Effective predictive and management approaches for species occurring in a metapopulation structure require good understanding of interpopulation connectivity. In this study, we ask whether population genetic structure of marine species with fragmented distributions can be predicted by stepping‐stone oceanographic transport and habitat continuity, using as model an ecosystem‐structuring brown alga, Cystoseira amentacea var. stricta. To answer this question, we analysed the genetic structure and estimated the connectivity of populations along discontinuous rocky habitat patches in southern Italy, using microsatellite markers at multiple scales. In addition, we modelled the effect of rocky habitat continuity and ocean circulation on gene flow by simulating Lagrangian particle dispersal based on ocean surface currents allowing multigenerational stepping‐stone dynamics. Populations were highly differentiated, at scales from few metres up to thousands of kilometres. The best possible model fit to explain the genetic results combined current direction, rocky habitat extension and distance along the coast among rocky sites. We conclude that a combination of variable suitable habitat and oceanographic transport is a useful predictor of genetic structure. This relationship provides insight into the mechanisms of dispersal and the role of life‐history traits. Our results highlight the importance of spatially explicit modelling of stepping‐stone dynamics and oceanographic directional transport coupled with habitat suitability, to better describe and predict marine population structure and differentiation. This study also suggests the appropriate spatial scales for the conservation, restoration and management of species that are increasingly affected by habitat modifications.     

Gap-crossing movements predict species occupancy in Amazonian forest fragments

In fragmented landscapes, species persistence within isolated habitat patches is governed by a myriad of species life‐history, habitat patch and landscape characteristics. We investigated the inter‐specific variation in non‐forest gap‐crossing abilities of an entire tropical forest‐dependent avifauna. We then related this measure of dispersal ability to species life‐history characteristics and occupancy data from 31 variable‐sized forest patches sampled within the same fragmented forest landscape. A total of 5436 gap‐crossing movements of 231 forest‐dependent bird species were observed across ten linear forest gaps of varying widths, adjacent to large areas of undisturbed forest. Species persistence in isolated fragments was strongly linked to gap‐crossing ability. The most capable gap‐crossers were medium to large‐bodied species in the large insectivore, frugivore and granivore guilds, matching the most prevalent subset of species in small forest patches. However, some competent gap‐crossing species failed to occur in small patches, and minimum forest‐patch area requirements were more important in determining patch occupancy for these species. Narrow forest gaps (4–70 m) created by roads and power‐lines may become territory boundaries, thereby eliminating home‐range gap‐crossing movements for many forest species, but permit rarer dispersal events. Wider gaps (>70 m) may inhibit gap‐crossing behaviour for all but the most vagile species. Although patch size and quality may be the most important factors in structuring species assemblages in forest fragments, our results show that the degree of patch isolation and permeability of the surrounding matrix also explain which species can persist in forest isolates. Reducing the number and width of forest‐dividing gaps; maintaining and/or creating forest corridors and increasing matrix permeability through the creation and maintenance of ‘stepping‐stone’ structures will maximise the species retention in fragmented tropical forest landscapes.     

Patch occupancy and abundance of local populations in landscapes differing in degree of habitat fragmentation: a case study of the colonial black‐headed gull,Chroicocephalus ridibundus

Aim This study investigated whether habitat fragmentation at the landscape level influences patch occupancy and abundance of the black‐headed gull, Chroicocephalus ridibundus, and whether the response of the species to environmental factors is consistent across replicated landscape plots. Location Water bodies (habitat patches) in southern Poland. Methods Surveys were conducted in two landscape types (four plots in each): (1) more‐fragmented landscape, in which habitat patches were small (mean size 2.2–6.2 ha) and far apart (mean distance 2.5–3.1 km); and (2) less‐fragmented landscape, in which habitat patches were large (mean size 9.2–16.5 ha) and separated by short distances (mean 0.9–1.4 km). Observations were performed twice in 284 potential habitat patches during the 2007 breeding season. Results Colonies were significantly more frequent and larger in the less‐fragmented landscapes than in the more‐fragmented ones. Probability of patch occupancy and number of breeding birds were positively related with patch size and these relationships were especially strong in the more‐fragmented landscapes. In the less‐fragmented landscapes, the occurrence of black‐headed gulls was negatively related to the distance to the nearest local population, but in the more‐fragmented landscapes such a relationship was not detected. As distance to the nearest habitat patch increased, the probability of the patch occupancy decreased in the more‐fragmented landscapes. Moreover, abundance was negatively influenced by distance to the nearest habitat patch, especially strongly in more‐fragmented landscapes. Proximity of corridors (rivers) positively influenced the occupation of patches regardless of landscape type. The number of islets positively influenced occupancy and abundance of local populations, and this relationship was stronger in the more‐fragmented landscapes. Main conclusions Our results are in agreement with predictions from metapopulation theory and are the first evidence that populations of black‐headed gulls may have a metapopulation structure. However, patch occupancy and abundance were differentially affected by explanatory variables in the more‐fragmented landscapes than in the less‐fragmented ones. This implies that it is impossible to derive, a priori, predictions about presence/abundance patterns based on only a single landscape.     

Increased inbreeding and strong kinship structure in Taxus baccata estimated from both AFLP and SSR data

Habitat fragmentation can have severe genetic consequences for trees, such as increased inbreeding and decreased effective population size. In effect, local populations suffer from reduction of genetic variation, and thus loss of adaptive capacity, which consequently increases their risk of extinction. In Europe, Taxus baccata is among a number of tree species experiencing strong habitat fragmentation. However, there is little empirical data on the population genetic consequences of fragmentation for this species. This study aimed to characterize local genetic structure in two natural remnants of English yew in Poland based on both amplified fragment length polymorphism (AFLP) and microsatellite (SSR) markers. We introduced a Bayesian approach that estimates the average inbreeding coefficient using AFLP (dominant) markers. Results showed that, in spite of high dispersal potential (bird-mediated seed dispersal and wind-mediated pollen dispersal), English yew populations show strong kinship structure, with a spatial extent of 50–100 m, depending on the population. The estimated inbreeding levels ranged from 0.016 to 0.063, depending on the population and marker used. Several patterns were evident: (1) AFLP markers showed stronger kinship structure than SSRs; (2) AFLP markers provided higher inbreeding estimates than SSRs; and (3) kinship structure and inbreeding were more pronounced in denser populations regardless of the marker used. Our results suggest that, because both kinship structure and (bi-parental) inbreeding exist in populations of English yew, gene dispersal can be fairly limited in this species. Furthermore, at a local scale, gene dispersal intensity can be more limited in a dense population.     

The fragmented population structure of the Sardinian chalk hill blue butterfly (Lepidoptera, Lycaenidae)

The endemic Sardinian chalk hill blue butterfly, Polyommatus coridon gennargenti, is considered vulnerable to extinction because of its low genetic variation and restricted distribution. The species also has a fragmented distribution, which follows the patchy distribution pattern of its larval host-plant. A preliminary investigation of the population structure of P. coridon gennargenti was carried out on a small network of four local populations by means of capture–recapture methods. Estimated population sizes and movement rates among the four adjacent local populations suggest that this taxon has a metapopulation structure composed of loosely connected small local populations. Natural fragmentation, isolation, and traditional land use contribute to the vulnerability of P. coridon gennargenti to extinction. Low effective population sizes and restricted movement between habitat patches lead to inbreeding and an increased vulnerability to extinction of this island population.     

Intertidal estimates of sea urchin abundance reveal congruence in spatial structure for a guild of consumers

We hypothesized congruence in the spatial structure of abundance data sampled across multiple scales for an ecological guild of consumers that exploit similar nutritional and habitat resources. We tested this hypothesis on the spatial organization of abundance of an herbivorous guild of sea urchins. We also examined whether the amount of local along‐shore rocky habitat can explain the observed spatial patterns of abundance. Standardized estimates of abundance of four intertidal sea urchins—Diadema cf. savignyi, Echinometra mathaei, Parechinus angulosus, and Stomopneustes variolaris—were determined by six observers at 105 sites across 2,850 km of coast of South Africa. For each species and observer, wavelet analysis was used on abundance estimates, after controlling for potential biases, to examine their spatial structure. The relationship between local sea urchin abundance and the amount of upstream and downstream rocky habitat, as defined by the prevailing ocean current, was also investigated. All species exhibited robust structure at scales of 75–220 km, despite variability among observers. Less robust structure in the abundances of three species was detected at larger scales of 430–898 km. Abundance estimates of sympatric populations of two species (D. cf. savignyi and E. mathaei) were positively correlated with the amount of rocky habitat upstream of the site, suggesting that upstream populations act as larval sources across a wide range of scales. No relationship between abundance and habitat size was found for P. angulosus or S. variolaris. Within the range of scales examined, we found robust congruence in spatial structure in abundance at the lower, but not the larger, range of scales for all four species. The relationship between abundance and upstream habitat availability in two species suggests that larval supply from upstream populations was probably the mechanism linking habitat size and abundance.     

Rampant drift in artificially fragmented populations of the endangered tidewater goby (Eucyclogobius newberryi)

Habitat fragmentation and its genetic consequences are a critically important issue in evaluating the evolutionary penalties of human habitat modification. Here, we examine the genetic structure and diversity in naturally subdivided and artificially fragmented populations of the endangered tidewater goby (Eucyclogobius newberryi), a small fish restricted to discrete coastal lagoons and estuaries in California, USA. We use five naturally fragmented coastal populations from a 300‐ km spatial scale as a standard to assess migration and drift relative to eight artificially fragmented bay populations from a 30‐ km spatial scale. Using nine microsatellite loci in 621 individuals, and a 522‐base fragment of mitochondrial DNA control region from 103 individuals, we found striking differences in the relative influences of migration and drift on genetic variation at these two scales. Overall, the artificially fragmented populations exhibited a consistent pattern of higher genetic differentiation and significantly lower genetic diversity relative to the naturally fragmented populations. Thus, even in a species characterized by habitat isolation and subdivision, further artificial fragmentation appears to result in substantial population genetic consequences and may not be sustainable.     

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.     

Genetic structure of a native cyprinid in a reservoir‐altered stream network

1. Reservoirs modify riverine ecosystems worldwide, and often with deleterious impacts on native biota. The immediate effects of reservoirs on native fish species below dams and in impounded reaches have received considerable attention, but it is unclear how reservoirs may affect fish species at larger spatial and temporal scales. Documented declines of stream fish populations in direct tributaries of reservoirs suggest reservoir pools may reduce gene flow among historically connected populations. 2. Because of increased predator densities in reservoirs and the extent of habitat alteration in impounded reaches, I predicted reservoir habitats would reduce gene flow among small‐bodied fish populations separated by reservoir habitat. I used microsatellite markers to assess the spatial genetic structure of populations of the red shiner (Cyprinella lutrensis), in a reservoir‐fragmented stream network (Lake Texoma, U.S.A.). I also tested the prediction that populations in two direct tributaries that have experienced population declines would have low genetic diversity. Individuals were collected from six sites upstream of the reservoir, three sites in the reservoir and three sites in direct tributaries of the reservoir during 2008 and 2009. 3. Results indicate that most populations were isolated by distance with little divergence among populations. In one direct tributary population, however, there was substantial genetic divergence, and genetic diversity was significantly lower than in other populations. Gene flow also seemed to be lower in reservoir habitats than in intact stream habitats, suggesting reservoir habitats may be reducing gene flow among the reservoir‐separated populations. These results indicate that reservoirs may reduce gene flow among reservoir‐fragmented stream fish populations, altering the evolutionary trajectories of fragmented populations.     

Spatial genetic structuring in a vagile species, the European wood mouse

We examined the genetic structure of natural populations of the European wood mouse Apodemus sylvaticus at the microgeographic (<3 km) and macrogeographic (>30 km) scales. Ecological and behavioural studies indicate that this species exhibits considerable dispersal relative to its home-range size. Thus, there is potential for high gene flow over larger geographic areas. As levels of population genetic structure are related to gene flow, we hypothesized that population genetic structuring at the microgeographic level should be negligible, increasing only with geographic distance. To test this, four sites were sampled within a microgeographic scale with two additional samples at the macrogeographic level. Individuals ( n =415) were screened and analysed for seven polymorphic microsatellite loci. Contrary to our hypothesis, significant levels of population structuring were detected at both scales. Comparing genetic differentiation with geographic distance suggests increasing genetic isolation with distance. However, this distance effect was non-significant being confounded by surprisingly high levels of differentiation among microgeographic samples. We attribute this pattern of genetic differentiation to the effect of habitat fragmentation, splitting large populations into components with small effective population sizes resulting in enhanced genetic drift. Our results indicate that it is incorrect to assume genetic homogeneity among populations even where there is no evidence of physical barriers and dispersal can occur freely. In the case of A. sylvaticus , it is not clear whether dispersal does not occur across habitat barriers or behavioural dispersal occurs without consequent gene flow.     

Population genetics and fitness in fragmented populations of the dioecious and endangered <Emphasis Type="Italic">Silene otites</Emphasis> (Caryophyllaceae)

Population fragmentation is often correlated with loss of genetic diversity and reduced fitness. Obligate out-crossing (dioecy) is expected to enhance genetic diversity, reduce genetic differentiation, and avoid inbreeding depression through frequent gene flow. However, in highly fragmented populations dioecy has only diminishing effects upon genetic structure as pollination limitations (e.g. flight distance of pollinators) most often restrict inter-population gene flow in insect pollinated species. In fragmented dry grasslands in northeastern Germany, we analysed genetic structure, fitness, and habitat quality of the endangered dioecious Silene otites (Caryophyllaceae). Using AFLP markers, a high level of differentiation among ten populations was found (F _st = 0.36), while the intra-population genetic diversities (H _E = 0.165–0.240) were similar as compared to hermaphroditic species. There was neither a correlation between geographic and genetic distance nor between genetic diversity and population size, which indicates reduced gene flow among populations and random genetic drift. Plant size was positively correlated with genetic diversity. Seed set and number of juveniles were positively related to population size. Higher total coverage resulted in reduced plant fitness, and the number of juveniles was negatively correlated to cryptogam cover. Additionally, we found a sex ratio bias towards more male plants in larger populations. Overall, our results indicate that on a regional geographic scale dioecy does not necessarily prevent genetic erosion in the case of habitat fragmentation, especially in the absence of long distance seed and pollen dispersal capacity.     

Population genetic structure of male black grouse (Tetrao tetrix L.) in fragmented vs. continuous landscapes

We investigated the association of habitat fragmentation with genetic structure of male black grouse Tetrao tetrix. Using 14 microsatellites, we compared the genetic differentiation of males among nine localities in continuous lowland habitats in Finland to the genetic differentiation among 14 localities in fragmented habitats in the Alps (France, Switzerland and Italy). In both areas, we found significant genetic differentiation. However, the average differentiation, measured as theta, was more than three times higher in the Alps than in Finland. The greater differentiation found in the Alps is probably due to the presence of mountain ridges rising above natural habitats of the species, which form barriers to gene flow, and to a higher influence of genetic drift resulting from lower effective sizes in highly fragmented habitats. The detection of isolation by distance in the Alps suggests that gene flow among populations does occur. The genetic variability measured as gene diversity HE and allelic richness A was lower in the Alps than in Finland. This could result from the higher fragmentation and/or from the fact that populations in the Alps are isolated from the main species range and have a lower effective size than in Finland. This study suggests that habitat fragmentation can affect genetic structure of avian species with relatively high dispersal propensities.     

Genetic differentiation and gene flow among populations of the alpine butterfly, Parnassius smintheus, vary with landscape connectivity

Levels of gene flow among populations vary both inter- and intraspecifically, and understanding the ecological bases of variation in levels of gene flow represents an important link between the ecological and evolutionary dynamics of populations. The effects of habitat spatial structure on gene flow have received considerable attention; however, most studies have been conducted at a single spatial scale and without background data on how individual movement is affected by landscape features. We examined the influence of habitat connectivity on inferred levels of gene flow in a high-altitude, meadow-dwelling butterfly, Parnassius smintheus. For this species, we had background data on the effects of landscape structure on both individual movement and on small-scale population genetic differentiation. We compared genetic differentiation and patterns of isolation by distance, based on variation at seven microsatellite loci, among three regions representing two levels of connectivity of high-altitude, nonforested habitats. We found that reduced connectivity of habitats, resulting from more forest cover at high altitudes, was associated with greater genetic differentiation among populations (higher estimated FST), a breakdown of isolation by distance, and overall lower levels of inferred gene flow. These observed differences were consistent with expectations based on our knowledge of the movement behaviour of this species and on previous population genetic analyses conducted at the smaller spatial scale. Our results indicate that the role of gene flow may vary among groups of populations depending on the interplay between individual movement and the structure of the surrounding landscape.     

Genetic variation and structure in the Mediterranean shrubs Myrtus communis and Pistacia lentiscus in different landscape contexts

Studies concerning different habitat configurations can provide insights into the complex interactions between species’ life‐history traits and the environment and can help to predict patterns in population genetics. In this study, we compared patterns of genetic variation in two Mediterranean shrub species (Myrtus communis and Pistacia lentiscus) that co‐occur in populations within three contrasting landscape contexts: continuous, fragmented‐connected and fragmented‐isolated populations. Analysing variation at microsatellites loci, our results revealed weak responses to the landscape contexts. We rather found a population‐specific response in both study species. However, despite both study species sharing similar levels of genetic diversity, Myrtus displayed higher levels of homozygosity and genetic differentiation among populations, stronger patterns of within‐population spatial genetic structure, lower values of mutation‐scaled effective population size and stronger evidence for recent genetic bottlenecks than Pistacia. This result highlights the influence of past events (e.g. historical connectivity, fluctuations in population size) and local factors (e.g. microhabitat availability for recruitment, habitat quality, plant density, native fauna) and that the landscape configuration per se (i.e. fragment size and/or isolation) might not completely determine the species’ genetic patterns.