Genetic diversity of wild populations of Rheum tanguticum endemic to China as revealed by ISSR analysis

Rheum tanguticum is an important but endangered traditional Chinese medicine endemic to China. The wild resources have been declining. Establishing the genetic diversity of the species would assist in its conservation and breeding program. Inter-simple sequence repeats (ISSR) markers were used to assess the genetic diversity and population genetic structure in 13 wild populations of R. tanguticum from Qinghai Province. Thirteen selected primers produced 329 discernible bands, with 326 (92.94%) being polymorphic, indicating high genetic diversity at the species level. The Nei's gene diversity (He) was estimated to be 0.1724 within populations (range 0.1026–0.2104), and 0.2689 at the species level. Analysis of molecular variance (AMOVA) showed that the genetic variation was found mainly within populations (71.02%), but variance among populations was only 28.98%. In addition, Nei's differentiation coefficients (G_st) was found to be high (0.3585), confirming the relatively high level of genetic differentiation among populations. Mantel test revealed a significant correlation between genetic and geographic distances (r = 0.573, P = 0.002), and the unweighted pair-group method using arithmetic average (UPGMA) clustering and Principal coordinates analysis (PCoA) demonstrated similar results. Meanwhile, the genetic diversity of R. tanguticum positively correlated with altitude and annual mean precipitation, but negatively correlated with latitude and annual mean temperature. This result might be an explanation that the natural distribution of R. tanguticum is limited to alpine cold areas. We propose conservation strategy and breeding program for this plant.     

Comparison of the extent of genetic variation of Vallisneria natans and its sympatric congener V. spinulosa in lakes of the middle-lower reaches of the Yangtze River

Vallisneria natans and Vallisneria spinulosa are two morphologically very similar and sympatrically dominant submerged macrophytes in lakes of the middle-lower reaches of the Yangtze River. Genetic variation was compared based on a total of 196 individuals from six V. natans populations and 201 individuals from seven V. spinulosa populations. Using eight ISSR primers, a total of 139 and 129 DNA fragments were generated with 121 being polymorphic in V. natans and 99 in V. spinulosa. The two species maintained higher genetic variation both at the species and population levels in comparison with other aquatic macrophytes. A higher level of genetic diversity among populations was found in V. natans than in V. spinulosa: the percentage of polymorphic loci (PPL) in V. natans was 52-62% vs. 38-47% in V. spinulosa; gene diversity (H) was 0.21 in V. natans vs. 0.17 in V. spinulosa.Both an analysis of molecular variance (AMOVA) and F-estimation (F_ST) showed that most of the total genetic variation resided within populations of both species (AMOVA: 85% and 80%; F_ST: 0.132 and 0.202), indicating low genetic differentiation between populations. Principal coordinates analysis (PCA) indicated evident gene flow between populations of both species. The outcrossing reproductive mode and pervasive gene flow might have played important roles in maintaining high genetic diversity and in shaping low population differentiation of the two Vallisneria species, while the extent of clonal growth might account for the different levels of population divergence between them.     

Genetic diversity of Astragalus nitidiflorus, a critically endangered endemic of SE Spain,and implications for its conservation

Inter-simple sequence repeats (ISSR) markers were used to assess the genetic diversity and population structure in five populations of Astragalus nitidiflorus, a critically endangered species endemic to southeast Spain. Eight primers amplified 78 bands with 40 (51.3%) being polymorphic. Statistical results indicated a low genetic diversity at the population and species level, with percentages of polymorphic bands (PPB) ranging from 28.2 to 37.2% (an average of 31.8%), and means of gene diversity (H_E) of 0.129 and 0.171 respectively. The Shannon’s index (SI) ranged from 0.160 to 0.214 at the population level and was 0.260 at the species level. A low level of genetic differentiation among populations was detected, based on the Shannon’s information index (0.297), the coefficient of genetic differentiation between populations (G_ST = 0.2418) and AMOVA analysis (Φ_ST = 0.255). The estimated gene flow (Nm) was 0.789. The high genetic connectivity found among populations of A. nitidiflorus is an evidence of a recent habitat fragmentation. In addition, a bottleneck event in the past has been revealed, with a subsequent reduction of population size and a loss of genetic variation. Based on these results, the conservation strategy of A. nitidiflorus was proposed.     

Genetic diversity of Sinojackia dolichocarpa (Styracaceae), a species endangered and endemic to China,detected by inter-simple sequence repeat (ISSR)

Sinojackia dolichocarpa, a species endangered and endemic to China, is distributed only in the regional area of Shimen and Sangzhi in Hunan Province. Inter-simple sequence repeat (ISSR) markers were used to investigate the genetic diversity within and among the four natural populations of S. dolichocarpa. Leaf samples were collected from 84 individuals. Thirteen ISSR primers selected from 80 primers gave rise to 137 discernible DNA bands of which 100 (72.99%) were polymorphic. On average each primer gave rise to 10.5 bands including 7.7 bands with polymorphic profile. At the species level, high genetic diversity was detected (PPB: 72.99%; H_E: 0.2255; Ho: 0.3453). However, relatively low genetic diversity existed within populations. Population in Maozhuhe (MZH) exhibits the greatest level of variability (PPB: 40.38%, H_E: 0.1566, Ho: 0.2330), whereas the population in Jingguanmen (JGM) finds its own variability at the lowest level (PPB: 30.66%; H_E: 0.1078; Ho: 0.1601). A high level of genetic differentiation among populations was revealed by Nei's gene diversity statistics (45.30%), Shannon's information measure (45.24%) and analysis of molecular variance (AMOVA) (52.88%). The main factors responsible for the high level of differentiation among populations are probably related to the selfing reproductive system and the isolation of populations. The strong genetic differentiation among populations indicates that the management for the conservation of genetic variability in S. dolichocarpa should aim to preserve every population.     

Genetic diversity within and among populations of the endangered and endemic species Primula merrilliana in China

Primula merrilliana Schltr. is an endangered and narrowly-distributed endemic species of southern Anhui Province in China. In this study, the level of genetic variation and the pattern of genetic structure in six natural populations of P. merrilliana were assessed by using ISSR (inter-simple sequence repeats) markers. Based on ten primers, 137 clear and reproducible DNA fragments were generated, of which 109 were polymorphic. The statistical results indicated that there was a relatively low genetic diversity within populations, and a high genetic differentiation among populations (G_ST = 0.53, Φ_ST = 0.49). The level of population genetic diversity was correlated to habitat type and the gene flow (N_m) was low with only 0.45. The unexpected genetic structure of P. merrilliana may be explained by limited gene flow that was caused by habitat fragmentation and limited seeds and pollen dispersal ability, self-compatible breeding system and biennial life form.     

Spatial Genetic Structure of Campanula sabatia, a Threatened Narrow Endemic Species of the Mediterranean Basin

Campanula sabatia is endemic to NW Italy (Liguria) and it was included in the European Red list of endangered species due to the heavy human pressure on its habitat. AFLP markers were used to detect the genetic diversity within and among ten populations (totalling 83 individuals) representing the range of the species. In spite of the limited distribution of this endemic taxon, high levels in percentage of polymorphic bands (PPB), gene diversity (H _S and H _T) and Shannon??s information index (I) were detected both at population (PPB?=?60?%, H _S?=?0.1853, I?=?0.2836) and at species level (PPB?=?100?%, H _T?=?0.2415, I?=?0.3871). The coefficient of genetic differentiation among population (G??_ST) was 0.1935, while the level of gene flow (N?? _m) was estimated to be 2.0832. AMOVA analysis identified a genetic variation within populations of 83?% of the total. Bayesian clustering methods assigned individuals to two geographical groups partly found within the same population, probably due to a high rate of genetic exchange among its populations. Conservation measures are suggested on the basis of the genetic diversity detected to ensure an effective protection for this endemic species.     

Mating system and population genetic structure of Bruguiera gymnorrhiza (Rhizophoraceae), a viviparous mangrove species in China

In order to explain the diversity patterns and develop the conservation strategies, the population genetic structures and the mating systems of Bruguiera gymnorrhiza from the coastlines of south China were investigated in this study. The mating system parameters were analyzed using progeny arrays for allozyme markers. The multilocus outcrossing rates (tm) ranged from 0.845 (Fugong) to 0.267 (Dongzhai harbor). High allozyme variations within the five collected populations were determined and compared with the published data of other plant species with the mixed mating systems. At species level, the percentage of polymorphic loci (P) was 80%, the average number of alleles per locus (A) was 2.440, and the heterozygosity (He) was 0.293. The total gene diversity within each population (H_S = 0.2782) and the coefficient of genetic differentiation (G_ST = 0.0579) among the populations were estimated. On the basis of this population genetic structure, it is suggested that the gene flow (Nm = 3.85) is quite high, which is possibly related to its water-dispersed hypocotyls. It is also suggested that the mating system of this species is of mixed mating.     

Genetic structure in peripheral Western European populations of the endangered species Cochlearia pyrenaica (Brassicaceae)

Cochlearia pyrenaica is one of the most endangered plant species in Europe, listed in many European and regional conservation policy documents (e.g. Spain, France, Belgium, Switzerland). To study its genetic structure, define its conservation units and propose a management strategy for this species, amplified fragment length polymorphism markers were used to analyse the genetic diversity within and between five representative populations of the species distribution in Western Europe (Cantabrian Range, North of Spain; Pyrenees, France; Wallonia, Belgium). Low levels of genetic diversity were revealed by the population percentage of polymorphic bands (PPB?=?36.56%), average within-population diversity (H _S?=?0.0990) and genetic diversity within populations (H _pop?=?0.1541), although high levels were reported at species level (PPB?=?81.16%; total genetic diversity for the species, H _T?=?0.0990; and genetic diversity within whole species, H _sp?=?0.2515). The coefficient of genetic differentiation among populations (G _ST) was 0.3869. The analysis of Shannon diversity index in population and for the total data set partitioned (38.72%) and AMOVA (53%) detected a high level of interpopulation diversity, in broad agreement with the result of genetic differentiation analysis. NeighborNet network and principal coordinate analyses clustered the populations in three major groups congruent with geographical regions. Bayesian clustering also confirmed these three distinct genetic clusters. The level of gene flow (Nm) was estimated as 0.3961 individuals per generation among populations, with the genetic identity (I) and genetic distance (D) among populations ranging from 0.8679 to 0.9651 and from 0.0355 to 0.1417, respectively. Therefore, the low levels of genetic variation and high divergence of regional gene pools indicate that there is a need to protect each disjunct region of Western Europe.     

PCR-RFLP analysis of mitochondrial DNA cytochrome b gene among Haruan (Channa striatus) in Malaysia

Haruan (Channa striatus) is in great demand in the Malaysian domestic fish market. In the present study, mtDNA cyt b was used to investigate genetic variation of C. striatus among populations in Peninsular Malaysia. The overall population of C. striatus demonstrated a high level of haplotype diversity (h) and a low-to-moderate level of nucleotide diversity (π). Analysis of molecular variance (AMOVA) results showed a significantly different genetic differentiation among 6 populations (F_ST = 0.37566, P = 0.01). Gene flow (Nm) was high and ranged from 0.32469 to infinity (∞). No significant relationship between genetic distance and geographic distance was detected. A UPGMA tree based on the distance matrix of net interpopulation nucleotide divergence (d_A) and haplotype network of mtDNA cyt b revealed that C. striatus is divided into 2 major clades. The neutrality and mismatch distribution tests for all populations suggested that C. striatus in the study areas had undergone population expansion. The estimated time of population expansion in the mtDNA cyt b of C. striatus populations occurred 0.72-6.19 million years ago. Genetic diversity of mtDNA cyt b and population structure among Haruan populations in Peninsular Malaysia will be useful in fisheries management for standardization for Good Agriculture Practices (GAP) in fish-farming technology, as well as providing the basis for Good Manufacturing Practices (GMP).     

Genetic differentiation in endangered Gynostemma pentaphyllum (Thunb.) Makino based on ISSR polymorphism and its implications for conservation

Studies were performed to investigate the genetic variation of 14 natural populations of Gynostemma pentaphyllum (Thunb.) Makino, an outcrossing clonal plant species in China, using inter-simple sequence repeat (ISSR) markers. Fourteen selected primers were used to amplify DNA samples from 140 individuals, and totally 194 loci were detected. The percentage of polymorphic bands (PPBs) showed that the genetic diversity was pretty high at the species level (PPB = 96.39%) but quite low at the population level (PPB = 1.03–25.26%). Shannon's information index (I) and Nei's gene diversity (h) displayed a similar trend to PPB. According to the hierarchical analysis of molecular variance (AMOVA) and Nei's analysis of gene diversity, the percentages of genetic variation among populations were 88.66 and 88.94%, respectively, indicating a high level of inter-population genetic differentiation. The low levels of genetic diversity within populations and high genetic differentiation among populations were assumed to result from the limited gene flow, the clonal nature and genetic drift. Based on the genetic data, effective conservation strategies were proposed for conserving this traditional Chinese medicinal herb. Concerning the management of G. pentaphyllum, we suggested that in situ conservation be an important and practical measure for maintaining the genetic diversity and that a possibly maximum number of populations be conserved. Populations EMS and HLT, in which particularly low levels of genetic variation were characterized, should be given the priority for ex situ conservation.     

Genetic diversity of a relict plant species, Ligularia sibirica (L.) Cass. (Asteraceae)

Rare plant species can be divided into naturally, ‘old rare’ species and anthropogenically, ‘new rare’ species. Many recent studies explored genetic diversity of ‘new rare’ species. Less is, however, known about genetic diversity of ‘old rare’ species. We examined isozyme genetic variability of 20 populations of an ‘old rare’ plant species, Ligularia sibirica (Asteraceae) in the Czech and Slovak Republic. It is a long-lived perennial herb with mixed-mating breeding system, widely distributed from East Asia to European Russia, with few isolated relict populations in the remaining part of Europe.The results showed high genetic diversity within populations (80.8%) and a low level of genetic differentiation (F_ST = 0.179). Genetic distance between populations correlated significantly with geographic distance. There was also a significant positive correlation between genetic diversity and population size. This is probably caused by destruction of habitats in last centuries and subsequent decrease of population size. Patterns of genetic diversity suggest that the recent distribution is a result of stepwise postglacial migration of the species and subsequent natural fragmentation.We conclude that L. sibirica populations preserve high levels of genetic diversity and are not yet threatened by genetic factors. However, this may change if changes in habitat conditions continue.     

陕西省林麝mtDNA D-loop区序列结构和种群遗传多样性

林麝(Moschus berezovskii)曾广泛分布于中国,由于盗猎和栖息地缩小,秦岭地区野生种群数量迅速下降,圈养繁殖种群已成立了几十年,但大多数圈养种群的遗传背景不清,种群规模增长非常缓慢。为了给这一物种的保护和管理提供有用的信息,调查了陕西省林麝1个圈养种群3个野生种群线粒体DNA(mt DNA)D-Loop 632 bp片段的遗传多样性和种群结构。在69个个体中其碱基组成为A+T的平均含量63.2%高于G+C含量36.8%,共检测到变异位点171个(约占总位点数的27.05%)。核苷酸多样性(Pi)为0.04424,平均核苷酸差异数(K)为19.908。69个个体分属32个单倍型,单倍型间的平均遗传距离(P)为0.070。32个单倍型构建的NJ系统树聚为3个分支,4个林麝群体中的单倍型是随机分布的。4个群体的平均遗传距离为0.043(标准误SE为0.005),凤县养殖场群体与留坝和陇县群体的亲缘关系较远。单倍型间的平均遗传距离为0.043,可见其遗传分化尚未达到种群分化的水平。结果表明,陕西省林麝群体mt DNA D-loop区序列存在着较丰富的变异和遗传多样性,凤县野生群体和凤县养殖场群体的核苷酸多样性和单倍型多样较高,养殖场种群没有出现近亲繁殖及遗传多样性下降的情况。凤县野生群体和凤县养殖场群体两者遗传分化较小,存在着较高的基因流水平。     

ISSR analysis of genetic diversity of the Qinghai-Tibet Plateau endemic Rhodiola chrysanthemifolia (Crassulaceae)

Inter-simple sequence repeat markers (ISSR) were used to estimate genetic diversity within and among 10 populations of Rhodiola chrysanthemifolia along Nianqingtangula Mountains and Brahmaputra, a species endemic to the Qinghai-Tibet Plateau and an endangered medicinal plant. Of the 100 primers screened, 13 produced highly polymorphic DNA fragments. Using these primers, 116 discernible DNA fragments were generated of which 104 (89.7%) were polymorphic, indicating substantial genetic diversity at the species level. Genetic diversity measured by the percentage of polymorphic bands (PPB) at the population level ranged from 21.97% to 48.8%. Analysis of molecular variance (AMOVA) showed that the genetic variation was found mainly among populations (77.3%), but no regional differentiation was discernible. Variance within populations was only 22.7%. The main factor responsible for this high level of differentiation among populations is probably the historical geographical and genetic isolation of populations in a harsh mountainous environment. Concerning the management of R. chrysanthemifolia, the high genetic differentiation of populations indicates the necessity of conserving the maximum possible number of populations.     

Genetic variability and differentiation of <Emphasis Type="Italic">Caragana microphylla</Emphasis> populations as revealed by RAPD markers

Genetic variability in random amplified polymorphic DNA (RAPD) was studied in 90 individuals of Caragana microphylla, an outcrossing perennial shrub species, from five natural populations sampled in Inner Mongolia steppe of China on a small scale. Nineteen selected primers were used to amplify DNA samples, and totally 225 bands were detected. The percentage of polymorphic bands within populations ranged form 58.22% to 63.56%, with an average of 60% at the population level and 71.11% at the species level, indicating relatively high genetic variations in C. microphylla species. Shannon’s information index (l) and Nei’s gene diversity (h) showed the similar trend with each other. According to the analysis of Nei’s gene diversity, the percentage of genetic variation among populations was 7.13%, indicating a low level of genetic differentiation among populations. There existed a strong gene flow (N _m = 3.26) among populations. Although AMOVA analysis also revealed most variation was within populations (Φ_ST = 4.1%), a significant proportion was observed among populations (P < 0.001) in the present study, suggesting genetic differentiation occurred among populations at a certain extent. Based on Mantel’s tests and the results of previous studies, the genetic structure pattern of C. microphylla accorded with the isolation-by-distance model on a very large scale, however, on a small scale, the significant genetic differentiation among populations might be enhanced by the micro-environmental divergence among the sampling sites, rather than by geographic factors. Analysis of the genetic variations of C. microphylla populations provided useful information for the adaptive strategy of Caragana species.     

Genetic diversity of populations of a rare species of Dictamnus gymnostylis Stev. in Bashkir Cis-Urals

We examined genetic diversity in populations of Dictamnus gymnostylis Stev., a rare species growing in the Bashkir Cis-Urals, based on the analysis of 8 gene-enzyme systems and detected a fairly high level of intraspecific genetic diversity and population differentiation. We determined the average number of alleles per locus (A) as 1.57; the portion of polymorphic loci (P ₉₅) as 0.508; the observed heterozygosity (H _o) as 0.139; and the expected heterozygosity (H _e) as 0.169. Of the total genetic diversity, 88.3% stems from variability within populations, and 11.7% is due to variation among populations. The average value for the Nei’s genetic distance (D) constituted 0.028.     

Chloroplast microsatellite diversity of Clintonia udensis (Liliaceae) populations in East Asia

Nineteen populations of Clintonia udensis Trautv. & Mey. were examined to quantify genetic diversity and genetic structure by chloroplast DNA microsatellites (cpSSR). Significant cpSSR genetic diversity (PPB = 63.64%) was detected in this species. Tetraploid populations demonstrated approximately the same level of genetic diversity as the diploid ones. A significant differentiation, however, was found between tetraploids and diploids. Most of the sixteen chloroplast haplotypes were limited to a single population. The level of haplotype diversity was high (Hd = 0.915). AMOVA, PCA and Bayesian clustering analysis revealed that there were significant genetic differences among populations. Inter-population genetic distances among population sites correlated significantly with geographic distances. These results indicate that the mixed-mating – breeding system, limited gene flow, environmental stress, and historical factors may be the main factors causing geographical differentiation in the genetic structure of C. udensis.     

Genetic diversity of Stipa grandis P.Smirn populations across the species’ range in the Inner Mongolia Plateau of China

In order to investigate the genetic diversity of dominant species under the background of climate change and grassland utilization in the Inner Mongolia Plateau of China, we sampled seven Stipa grandis populations along an increasing aridity gradient in the present study. The Nei’s gene diversity of populations (He) was estimated to be 0.15 and the percentage of polymorphic loci (PPL) was 49.28%. The genetic differentiation among populations (Φ_ST) was 0.2431. There was a significant relationship between genetic distance and geographic distance among the S. grandis populations by Mantel's test. The genetic diversity was significantly correlated with longitude and annual mean precipitation, which suggested that a combination of climatic factors affected the genetic diversity. The populations in the marginal habitat should be paid more attention because of their low genetic diversity and its significance for conservation of the whole species.     

High genetic differentiation and low genetic diversity in Incarvillea younghusbandii, an endemic plant of Qinghai-Tibetan Plateau,revealed by AFLP markers

Incarvillea younghusbandii Sprague (Bignoniaceae) is a perennial herbaceous plant endemic to Qinghai-Tibetan Plateau. As a species of medical and horticulture importance, I. younghusbandii is threatened by over exploitation and habitat fragmentation. In this study, we analyze the genetic diversity and population structure of I. younghusbandii using amplified fragment length polymorphism (AFLP) markers. Our data reveal very low levels of genetic diversity in seven natural populations across Tibet. Specifically, at population level, the average Nei's genetic diversity index (H_E) and Shannon's diversity index (I) were 0.063 and 0.096, respectively. In contrast, high genetic differentiation among populations (Gst = 0.6238, Φ_ST = 0.614) is detected. The results of Neighbor-joining cluster, PCO, and STRUCTURE assignment reveal consistent pattern, suggesting seven well-defined genetic groups that are concordant with their geographical origins. The possible mechanisms and implications of these findings for conservation are discussed.     

Molecular and quantitative trait variation within and among small fragmented populations of the endangered plant species Psilopeganum sinense

Background and Aims

Natural selection and genetic drift are important evolutionary forces in determining genetic and phenotypic differentiation in plant populations. The extent to which these two distinct evolutionary forces affect locally adaptive quantitative traits has been well studied in common plant and animal species. However, we know less about how quantitative traits respond to selection pressures and drift in endangered species that have small population sizes and fragmented distributions. To address this question, this study assessed the relative strengths of selection and genetic drift in shaping population differentiation of phenotypic traits in Psilopeganum sinense, a naturally rare and recently endangered plant species.

Methods

Population differentiation at five quantitative traits (Q_ST) obtained from a common garden experiment was compared with differentiation at putatively neutral microsatellite markers (F_ST) in seven populations of P. sinense. Q_ST estimates were derived using a Bayesian hierarchical variance component method.

Key Results

Trait-specific Q_ST values were equal to or lower than F_ST. Neutral genetic diversity was not correlated with quantitative genetic variation within the populations of P. sinense.

Conclusions

Despite the prevalent empirical evidence for Q_ST > F_ST, the results instead suggest a definitive role of stabilizing selection and drift leading to phenotypic differentiation among small populations. Three traits exhibited a significantly lower Q_ST relative to F_ST, suggesting that populations of P. sinense might have experienced stabilizing selection for the same optimal phenotypes despite large geographical distances between populations and habitat fragmentation. For the other two traits, Q_ST estimates were of the same magnitude as F_ST, indicating that divergence in these traits could have been achieved by genetic drift alone. The lack of correlation between molecular marker and quantitative genetic variation suggests that sophisticated considerations are required for the inference of conservation measures of P. sinense from neutral genetic markers.     

High genetic diversity and significant population structure in Cedrus brevifolia Henry, a narrow endemic Mediterranean tree from Cyprus

Endemic island plant species with a narrow distribution are often, but not always, linked to low genetic variation within populations and a lack of differentiation among populations. Cedrus brevifolia is a narrow endemic island tree species of Cyprus. Its range is restricted to a single forest, divided into five neighbouring sites. This study, using biparentally inherited nuclear microsatellites and paternally inherited plastid (chloroplast) microsatellites, assessed the genetic variation of C. brevifolia within its sole population and the level of genetic differentiation among formed sites. The results from both markers showed high diversity (nuclear H _T?=?0.70; plastid H _T?=?0.93), strongly suggesting that the species did not experience severe bottleneck events or extensive genetic drift. Besides, the maintenance of a high genetic diversity in C. brevifolia may suggest that it originates from a widespread congener species. Significant genetic differentiation at nuclear (G _ST?=?0.052) and plastid (G _ST?=?0.119) markers was found among the formed sites. Remarkably, the relatively high genetic differentiation found at plastid markers was comparable to values observed in two widespread congener cedar species. The genetic differentiation probably occurred due to fragmentation of a previously uniform population. This would lead to the shaping of different genetic groups (Bayesian analysis) and to significant population substructure. Furthermore, significant values observed for both isolation by distance and large-scale spatial genetic structure could indicate ineffective gene flow among sites and the early geographical isolation of the more isolated sites from the core population.