Isolation-by-distance in landscapes: considerations for landscape genetics

In landscape genetics, isolation-by-distance (IBD) is regarded as a baseline pattern that is obtained without additional effects of landscape elements on gene flow. However, the configuration of suitable habitat patches determines deme topology, which in turn should affect rates of gene flow. IBD patterns can be characterized either by monotonically increasing pairwise genetic differentiation (for example, F_ST) with increasing interdeme geographic distance (case-I pattern) or by monotonically increasing pairwise genetic differentiation up to a certain geographical distance beyond which no correlation is detectable anymore (case-IV pattern). We investigated if landscape configuration influenced the rate at which a case-IV pattern changed to a case-I pattern. We also determined at what interdeme distance the highest correlation was measured between genetic differentiation and geographic distance and whether this distance corresponded to the maximum migration distance. We set up a population genetic simulation study and assessed the development of IBD patterns for several habitat configurations and maximum migration distances. We show that the rate and likelihood of the transition of case-IV to case-I F_ST–distance relationships was strongly influenced by habitat configuration and maximum migration distance. We also found that the maximum correlation between genetic differentiation and geographic distance was not related to the maximum migration distance and was measured across all deme pairs in a case-I pattern and, for a case-IV pattern, at the distance where the F_ST–distance curve flattens out. We argue that in landscape genetics, separate analyses should be performed to either assess IBD or the landscape effects on gene flow.     

Historical influences dominate the population genetic structure of a sedentary marine fish,Atlantic wolffish (Anarhichas lupus), across the North Atlantic Ocean

Genetic variation was assessed in Atlantic wolffish, Anarhichas lupus, across the North Atlantic Ocean using microsatellite and amplified fragment length polymorphism (AFLP) markers. Despite unusual life history attributes such as large benthic eggs, large larvae, a limited pelagic stage and relatively sedentary adults, which suggest potential for strong population structure, range‐wide F_ST values were comparable to other marine fishes (≤0.035). Nevertheless, both significant genetic differentiation among regions and isolation by distance were observed, suggesting limited dispersal in this species. AFLP loci, evaluated on a subset of samples, revealed slightly higher F_ST values, but similar patterns of differentiation and isolation‐by‐distance estimates, compared to microsatellites. The genetic structure of Atlantic wolffish has likely been shaped by its post‐glacial history of recolonization, North Atlantic current patterns and continuity of habitat on continental shelves.     

On the Apportionment of Population Structure

Measures of population differentiation, such as F_ST, are traditionally derived from the partition of diversity within and between populations. However, the emergence of population clusters from multilocus analysis is a function of genetic structure (departures from panmixia) rather than of diversity. If the populations are close to panmixia, slight differences between the mean pairwise distance within and between populations (low F_ST) can manifest as strong separation between the populations, thus population clusters are often evident even when the vast majority of diversity is partitioned within populations rather than between them. For any given F_ST value, clusters can be tighter (more panmictic) or looser (more stratified), and in this respect higher F_ST does not always imply stronger differentiation. In this study we propose a measure for the partition of structure, denoted E_ST, which is more consistent with results from clustering schemes. Crucially, our measure is based on a statistic of the data that is a good measure of internal structure, mimicking the information extracted by unsupervised clustering or dimensionality reduction schemes. To assess the utility of our metric, we ranked various human (HGDP) population pairs based on F_ST and E_ST and found substantial differences in ranking order. E_ST ranking seems more consistent with population clustering and classification and possibly with geographic distance between populations. Thus, E_ST may at times outperform F_ST in identifying evolutionary significant differentiation.     

A meta-analysis of the effects of plant traits and geographical scale on the magnitude of adaptive differentiation as measured by the difference between QST and FST

The impact of directional selection on specific trait types in plant species, and how a species’ life history mediates this response to selection remains understudied. Discovering such interactions is however crucial for understanding the interplay between ecological and genetic processes underlying local adaptation in plants, and to evaluate a species’ evolutionary potential with respect to changing environments. Furthermore, it remains unclear whether the degree of adaptive differentiation generally increases with the geographical distance between plant populations. Here, we present a weighted mixed model based meta-analysis aimed at unraveling the potential interactions between plant trait types, life history characteristics and Q_ST–F_ST comparisons, and assessing the effect of geographical scale on population differentiation. Based on 51 studies we found that Q_ST values exceeded their corresponding F_ST values in 71.74 % out of 401 cases. Furthermore, different trait types were found to be differently susceptible to natural selection and the magnitude of Q_ST–F_ST comparisons was mediated by a plant species’ life span. These findings may be closely related to the genetic architectures of trait types and life histories, with the proportion of large-effect genes likely shaping the response to natural selection. Q_ST–F_ST values also increased with increasing distance between populations, pinpointing the combined effects of environmental differentiation and isolation by distance on the magnitude of population divergence. Finally, our model showed an inverse relationship between F_ST and Q_ST–F_ST values, presumably resulting from isolation by distance, the exchange of advantageous alleles, or genetic correlations among traits.     

Genetically distinct populations of northern shrimp,Pandalus borealis,in the North Atlantic: adaptation to different temperatures as an isolation factor

The large‐scale population genetic structure of northern shrimp, Pandalus borealis, was investigated over the species’ range in the North Atlantic, identifying multiple genetically distinct groups. Genetic divergence among sample localities varied among 10 microsatellite loci (range: F_ST = ?0.0002 to 0.0475) with a highly significant average (F_ST = 0.0149; P < 0.0001). In contrast, little or no genetic differences were observed among temporal replicates from the same localities (F_ST = 0.0004; P = 0.33). Spatial genetic patterns were compared to geographic distances, patterns of larval drift obtained through oceanographic modelling, and temperature differences, within a multiple linear regression framework. The best‐fit model included all three factors and explained approximately 29% of all spatial genetic divergence. However, geographic distance and larval drift alone had only minor effects (2.5–4.7%) on large‐scale genetic differentiation patterns, whereas bottom temperature differences explained most (26%). Larval drift was found to promote genetic homogeneity in parts of the study area with strong currents, but appeared ineffective across large temperature gradients. These findings highlight the breakdown of gene flow in a species with a long pelagic larval phase (up to 3 months) and indicate a role for local adaptation to temperature conditions in promoting evolutionary diversification and speciation in the marine environment.     

Genetic and phenotypic variation in a colourful treefrog across five geographic barriers

Aim Understanding the patterns and processes underlying phenotype in a polytypic species provides key insights into microevolutionary mechanisms of diversification. The red‐eyed treefrog, Agalychnis callidryas, exhibits strong regional differentiation in colour pattern, corresponding to five admixed mitochondrial DNA clades. We evaluated spatial diversity patterns across multiple, putative barriers to examine the fine‐scale processes that mediate phenotypic divergence between some regions while maintaining homogeneity between others. Location We examined patterns of phenotypic diversification among 17 sites that span five putative biogeographic barriers in lower Central America (Costa Rica and Panama). Methods We tested the extent to which genetic distance (F_ST) derived from six multilocus nuclear genotypes covaried with measures of phenotypic distance (leg coloration) within and between biogeographic regions. We used linear regression analyses to determine the role of geographic and genetic factors in structuring spatial patterns of phenotypic diversity. Results The factors that best explained patterns of phenotypic diversity varied among biogeographic regions. We identified one geographic barrier that impeded gene exchange and resulted in concordant phenotypic divergence across the Continental Divide, isolating Caribbean and Pacific populations. Across Caribbean Costa Rican populations, one barrier structured phenotypic but not genetic diversity patterns, indicating a role for selection. In other regions, the putative barriers had no determining effect on either genetic or leg colour structure. Main conclusions The processes mediating the distribution and diversification of colour pattern in this polytypic, wide‐ranging treefrog varied among biogeographic regions. Spatially varying selection combined with the isolating effects of geographic factors probably resulted in the patchy distribution of colour diversity across Costa Rican and Panamanian populations.     

Is life history a barrier to dispersal? Contrasting patterns of genetic differentiation along an oceanographically complex coast

Extreme variation in early life‐history strategies is considered a moderately good predictor of genetic subdivision and hence dispersal for a range of marine species. In reality, however, a good deal of population differentiation must reflect historical effects, more subtle variation in life histories, and, particularly, the interaction of larvae with oceanographic processes. Using a combination of allozyme and microsatellite markers, we show that the large‐scale genetic structure of populations of three species (direct and planktonically developing cushion stars and a planktonic developing sea anemone that is also asexually viviparous) varies consistently, in line with the predicted capacity for dispersal within three geographic regions. We detected high levels of genetic subdivision for the direct developing cushion star (F_ST = 0.6), low levels for the planktonically developing cushion star (F_ST = 0.009), and intermediate levels for the sexual/asexual sea anmone (F_ST = 0.19). These patterns are exhibited despite the highly variable patterns of current movement and the presence of biogeographic barriers. Our results suggest that, although there is large scale genetic differentiation for two species, patterns of population connectivity are remarkably consistent within major regions and do not reflect variation in major oceanographic processes or genetic discontinuity coincident with biogeographic boundaries. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 106–116.     

Natural selection contributes to geographic patterns of thermal plasticity in Plantago lanceolata

A long‐standing debate in evolutionary biology concerns the relative importance of different evolutionary forces in explaining phenotypic diversification at large geographic scales. For example, natural selection is typically assumed to underlie divergence along environmental gradients. However, neutral evolutionary processes can produce similar patterns. We collected molecular genetic data from 14 European populations of Plantago lanceolata to test the contributions of natural selection versus neutral evolution to population divergence in temperature‐sensitive phenotypic plasticity of floral reflectance. In P. lanceolata, reflectance plasticity is positively correlated with latitude/altitude. We used population pairwise comparisons between neutral genetic differentiation (F_ST and Jost's D) and phenotypic differentiation (P_ST) to assess the contributions of geographic distance and environmental parameters of the reproductive season in driving population divergence. Data are consistent with selection having shaped large‐scale geographic patterns in thermal plasticity. The aggregate pattern of P_ST versus F_ST was consistent with divergent selection. F_ST explained thermal plasticity differences only when geographic distance was not included in the model. Differences in the extent of cool reproductive season temperatures, and not overall temperature variation, explained plasticity differences independent of distance. Results are consistent with the hypothesis that thermal plasticity is adaptive where growing seasons are shorter and cooler, that is, at high latitude/altitude.     

Genetic variability of the stable fly assessed on a global scale using amplified fragment length polymorphism

The stable fly, Stomoxys calcitrans (L.) (Diptera: Muscidae), is a blood‐feeding, economically important pest of animals and humans worldwide. Improved management strategies are essential and their development would benefit from studies on genetic diversity of stable flies. Especially if done on a global scale, such research could generate information necessary for the development and application of more efficient control methods. Herein we report on a genetic study of stable flies using amplified fragment length polymorphism, with samples of 10–40 individuals acquired from a total of 25 locations in the Nearctic, Neotropic, Palearctic, Afrotropic and Australasian biogeographical regions. We hypothesized that genetic differentiation would exist across geographical barriers. Although F_ST (0.33) was moderately high, the G_ST (0.05; representing genetic diversity between individuals) was very low; N_m values (representing gene flow) were high (9.36). The mismatch distribution and tests of neutrality suggested population expansion, with no genetic differentiation between locations. The analysis of molecular variance (AMOVA) results showed the majority of genetic diversity was within groups. The mantel test showed no correlation between geographic and genetic distance; this strongly supports the AMOVA results. These results suggest that stable flies did not show genetic differentiation but are panmictic, with no evidence of isolation by distance or across geographical barriers.     

Genetic variability and structure of the water vole Arvicola amphibius across four metapopulations in northern Norway

Water vole Arvicola amphibius populations have recently experienced severe decline in several European countries as a consequence of both reduction in suitable habitat and the establishment of the alien predator American mink Neovison vison. We used DNA microsatellite markers to describe the genetic structure of 14 island populations of water vole off the coast of northern Norway. We looked at intra‐ and inter‐population levels of genetic variation and examined the effect of distance among pairs of populations on genetic differentiation (isolation by distance). We found a high level of genetic differentiation (measured by F_ST) among populations overall as well as between all pairs of populations. The genetic differentiation between populations was positively correlated with geographic distance between them. A clustering analysis grouped individuals into 7 distinct clusters and showed the presence of 3 immigrants among them. Our results suggest a small geographic scale for evolutionary and population dynamic processes in our water vole populations.     

Geographic Influence on Genetic Structure in the Widespread Neotropical Tree Simarouba amara (Simaroubaceae)

This study aimed to assess the population genetic structure of a widespread Neotropical tree species, Simarouba amara, at local, regional and continental spatial scales. We used five microsatellite loci to examine genetic variation in 14 natural populations (N?=?478 individuals) of this vertebrate-dispersed rain forest tree species in Panama, Ecuador, and French Guiana. Estimates of genetic differentiation (F _st and R _st) were significant among all but one population pair and global differentiation was moderate (F _st?=?0.25, R _st?=?0.33) with 94% of genetic variation ascribed to differences among three main geographic regions (Central America, Western Ecuador, Amazon basin). There was no evidence of isolation by distance within regions. Allele-size mutations contributed significantly (R _ST > F _ST) to the divergences between cis- and trans-Andean populations, highlighting the role of the northern Andean cordilleras as an important geographic barrier for this species.     

Positive correlation between vegetation dissimilarity and genetic differentiation of Carex sempervirens

Synthesizing genetic data at population level and vegetation data at community level may give insight into how ecological and evolutionary processes associated with different vegetation influence genetic diversity and differentiation of plant populations. We correlated population genetic patterns of Carex sempervirens with community vegetation patterns in abandoned subalpine grassland in the Swiss Alps. Within-population genetic diversity (percentage of polymorphic bands and Nei's gene diversity) of C. sempervirens was not significantly correlated with plant richness, evenness or Shannon's diversity index (Pearson correlation coefficient |r|<0.32, P>0.10). However, the genetic distance (F_ST) between C. sempervirens populations was significantly positively correlated with the vegetation dissimilarity between communities (Mantel's r=0.23, P<0.01). The correlation between the population genetic differentiation and the vegetation dissimilarity was not due to the parallel effects of geographic isolation or site conditions, because F_ST was not correlated with the geographic distance or the pairwise differences in any of the measured site condition parameters. One likely mechanism is that different plant communities were associated with different selective forces, which, in turn, influenced the genetic differentiation between C. sempervirens populations. Another possibility is that the vegetation heterogeneity (dissimilarity) generated ecological barriers against gene flow and thus enhanced the genetic differentiation between C. sempervirens populations.     

Environmental and climatic determinants of molecular diversity and genetic population structure in a coenagrionid damselfly

Identifying environmental factors that structure intraspecific genetic diversity is of interest for both habitat preservation and biodiversity conservation. Recent advances in statistical and geographical genetics make it possible to investigate how environmental factors affect geographic organisation and population structure of molecular genetic diversity within species. Here we present a study on a common and wide ranging insect, the blue tailed damselfly Ischnuraelegans, which has been the target of many ecological and evolutionary studies. We addressed the following questions: (i) Is the population structure affected by longitudinal or latitudinal gradients?; (ii) Do geographic boundaries limit gene flow?; (iii) Does geographic distance affect connectivity and is there a signature of past bottlenecks?; (iv) Is there evidence of a recent range expansion and (vi) what is the effect of geography and climatic factors on population structure? We found low to moderate genetic sub-structuring between populations (mean F_ST = 0.06, D_est = 0.12), and an effect of longitude, but not latitude, on genetic diversity. No significant effects of geographic boundaries (e.g. water bodies) were found. F_ST-and D_est-values increased with geographic distance; however, there was no evidence for recent bottlenecks. Finally, we did not detect any molecular signatures of range expansions or an effect of geographic suitability, although local precipitation had a strong effect on genetic differentiation. The population structure of this small insect has probably been shaped by ecological factors that are correlated with longitudinal gradients, geographic distances, and local precipitation. The relatively weak global population structure and high degree of genetic variation within populations suggest that I. elegans has high dispersal ability, which is consistent with this species being an effective and early coloniser of new habitats.     

Life history,larval dispersal,and connectivity in coral reef fish among the Scattered Islands of the Mozambique Channel

The Western Indian Ocean harbors one of the world’s most diverse marine biota yet is threatened by exploitation with few conservation measures in place. Primary candidates for conservation in the region are the Scattered Islands (Îles Éparses), a group of relatively pristine and uninhabited islands in the Mozambique Channel. However, while optimal conservation strategies depend on the degree of population connectivity among spatially isolated habitats, very few studies have been conducted in the area. Here, we use highly variable microsatellite markers from two damselfishes (Amphiprion akallopisos and Dascyllus trimaculatus) with differing life history traits [pelagic larval duration (PLD), adult habitat] to compare genetic structure and connectivity among these islands using classic population structure indices as well as Bayesian clustering methods. All classical fixation indexes F _ST, R _ST, G′_ST, and Jost’s D show stronger genetic differentiation among islands for A. akallopisos compared to D. trimaculatus, consistent with the former species’ shorter PLD and stronger adult site attachment, which may restrict larval dispersal potential. In agreement with these results, the Bayesian analysis revealed clear genetic differentiation among the islands in A. akallopisos, separating the southern group (Bassas da India and Europa) from the center (Juan de Nova) and northern (Îles Glorieuses) islands, but not for D. trimaculatus. Local oceanographic patterns such as eddies that occur along the Mozambique Channel appear to parallel the results reported for A. akallopisos, but such features seem to have little effect on the genetic differentiation of D. trimaculatus. The contrasting patterns of genetic differentiation between species within the same family highlight the importance of accounting for diverse life history traits when assessing community-wide connectivity, an increasingly common consideration in conservation planning.

     

The Distribution and Host Shifts of Cotton-Melon Aphids in Northern China

Aphis gossypii Glover (Hemiptera: Aphididae) is a serious pest of cotton in northern China. A microsatellite analysis was used to characterize the genetic structure of A. gossypii populations from different geographic, host plant, and seasonal populations in 2014. Among 906 individuals, 507 multilocus genotypes were identified, with genotypic richness values of 0.07–1.00 for the populations. We observed moderate levels of genetic differentiation among geographic populations (F_ST = 0.103; 95% confidence interval: 0.065–0.145) and host plant populations (F_ST = 0.237; 95% confidence interval: 0.187–0.296). A Mantel test of isolation by distance revealed no significant correlations between Slatkin’s linearized F_ST and the natural logarithm of geographic distance. A Bayesian analysis of population genetic structures identified three clusters. An analysis of molecular variance revealed significant differences among the three clusters (F = 0.26596, P < 0.0001), among seasons (F = 0.04244, P = 0.00381), and among host populations (F = 0.12975, P = 0.0029). Thus, the A. gossypii populations in northern China exhibit considerable genotypic diversity. Additionally, our findings indicated that the 31 analyzed populations could be classified as one of three host biotypes (i.e., cotton, cucumber, and pomegranate biotypes). There were also clear seasonal effects on population genetic structure diversity among aphids collected from Anyang.     

Gene flow in malaria vector Anopheles stephensi (Diptera: Culicidae) across the Aravalli Hills in North-west India

The population structure of An. stephensi in North-west India was studied to assess the impact of the Aravalli Hills, as a barrier to gene flow using microsatellite markers. Large and significant genetic differentiation was found along the sides of, as well as across, the Aravalli Hills as the mean F_ST and R_ST on west vs. east of the Aravalli Hills were 0.213, 0.112 and 0.179, 0.056, respectively. Similarly, across the hills, mean values of F_ST and R_ST were 0.100 and 0.094, respectively. Genetic diversity on both sides did not vary significantly. The F_ST values were more sensitive than R_ST values, indicating that genetic drift might have caused genetic differentiation between populations. A positive correlation (r = 0.0149 and 0.157, respective to F_ST and R_ST) was found between genetic differentiations and geographic distances irrespective of the hills. Low level of gene flow was found along both sides (Nm = 0.92 and 0.14; west vs. east of Aravalli Hills, respectively) as compared to across the Aravalli Hills (Nm = 2.25). It was found that the Aravalli Hills are not working as an effective barrier to gene flow for An. Stephensi, maybe because of the low average height and discontinuous hills, however, the distance is playing a major role for differentiation between populations due to active mode of dispersal of An. stephensi mosquitoes which have a short flight range. All this information should help draw the strategies for genetic control of mosquitoes using transgenic mosquitoes.     

Population structure of the malaria vector Anopheles sinensis (Diptera: Culicidae) in China: two gene pools inferred by microsatellites

Background

Anopheles sinensis is a competent malaria vector in China. An understanding of vector population structure is important to the vector-based malaria control programs. However, there is no adequate data of A. sinensis population genetics available yet.

Methodology/Principal Findings

This study used 5 microsatellite loci to estimate population genetic diversity, genetic differentiation and demographic history of A. sinensis from 14 representative localities in China. All 5 microsatellite loci were highly polymorphic across populations, with high allelic richness and heterozygosity. Hardy–Weinberg disequilibrium was found in 12 populations associated with heterozygote deficits, which was likely caused by the presence of null allele and the Wahlund effect. Bayesian clustering analysis revealed two gene pools, grouping samples into two population clusters; one includes six and the other includes eight populations. Out of 14 samples, six samples were mixed with individuals from both gene pools, indicating the coexistence of two genetic units in the areas sampled. The overall differentiation between two genetic pools was moderate (F _ST = 0.156). Pairwise differentiation between populations were lower within clusters (F _ST = 0.008–0.028 in cluster I and F _ST = 0.004–0.048 in cluster II) than between clusters (F _ST = 0.120–0.201). A reduced gene flow (Nm = 1–1.7) was detected between clusters. No evidence of isolation by distance was detected among populations neither within nor between the two clusters. There are differences in effective population size (Ne = 14.3-infinite) across sampled populations.

Conclusions/Significance

Two genetic pools with moderate genetic differentiation were identified in the A. sinensis populations in China. The population divergence was not correlated with geographic distance or barrier in the range. Variable effective population size and other demographic effects of historical population perturbations could be the factors affecting the population differentiation. The structured populations may limit the migration of genes under pressures/selections, such as insecticides and immune genes against malaria.     

High gene flow between populations of Macrotermes michaelseni (Isoptera, Termitidae)

Termite alates are thought to be poor active flyers, and this should lead to considerable genetic differentiation on small spatial scales. However, using four microsatellite loci for the termite Macrotermes michaelseni we found low values of genetic differentiation (F_ST) across a spatial scale of even more than 50 km. Genetic differentiation between populations increased with spatial distance up to 50 km. Furthermore, up to this distance, the scatter around the linear regression of genetic differentiation versus spatial distance increased with spatial distance. This suggests that across such spatial distances gene flow and genetic drift are of about equal importance, and near equilibrium. Using a regional F_ST as well as the distance between populations with non-significant F_ST-values (up to 25 km), gene flow is sufficiently high so that populations may be regarded as panmictic on spatial scales of 25 to 50 km. The apparent contradiction between dispersal distances observed in the field and estimates of gene flow from genetic markers may be due to the masses of swarming alates. Assuming a leptokurtic distribution of dispersal distances, atleast some alates are expected to travel considerable distances, most likely by passive drift. Received 25 January 2005; revised 11 April 2005; accepted 26 April 2005.     

Analysis of F ST outliers at allozyme loci in Pacific salmon: implications for natural selection

Natural selection has been invoked to explain the observed geographic distribution of allozyme allele frequencies for a number of teleost species. The effects of selection on allozyme loci in three species of Pacific salmon were tested. A simulation-based approach to estimate the null distribution of population differentiation (F _ST) and test for F _ST outliers was used. This approach showed that a majority of allozyme loci conform to neutral expectations predicted by the simulation model, with relatively few F _ST outliers found. No consistent F _ST outlier loci were found across species. Analysis of population sub-groups based on geography and genetic identity reduced the number of outlier loci for some species, indicating that large geographic groups may include genetically divergent populations and/or that there is geographic heterogeneity in selection pressure upon allozyme loci. Two outlier allozyme loci found in this analysis, lactate dehydrogenase-B and malic enzyme, have been shown to be influenced by selection in other teleost species. This approach is also useful in identifying allozyme loci (or other genetic markers) that meet assumptions for population genetic study.     

Little to no genetic structure in the ectomycorrhizal basidiomycete Suillus spraguei (Syn. S. pictus) across parts of the northeastern USA

Population genetic studies of ectomycorrhizal fungi at scales larger than 100 km are still relatively rare with highly variable results. In this study, we determined the population genetic structure of the ectomycorrhizal basidiomycete Suillus spraguei over distances up to 600 km in northeastern USA forests. S. spraguei associates exclusively with five-needled pines and only with white pine (Pinus strobus) in the eastern USA. We used six microsatellite loci to assess the genetic structure between eight sites sampled in the Adirondack Park of New York and seven sites sampled in other forests of New York, Pennsylvania, and Massachusetts. Except for one site, little to no genetic differentiation was detected in pairwise comparisons of the sites (F _ST?=?0 to 0.05). Only one site was moderately differentiated from most other sites (F _ST?=?0.02 to 0.15). The Mantel test showed no significant correlation between genetic and geographic distances (isolation by distance; R ²?=?0.003, P?=?0.3). The STRUCTURE analysis also supported the presence of a single cluster (K?=?1).