Genetic variability of the narrow endemic Rhamnus persicifolia Moris (Rhamnaceae) and its implications for conservation

Rhamnus persicifolia Moris is an endemic small tree belonging to the Rhamnus cathartica group, growing along mountainous streams of Central-Eastern Sardinia (Italy). ISSR markers were used to detect the genetic diversity within and among six populations representative of the species distribution range. In spite of the limited distribution of this endemic taxon, fairly high levels of genetic diversity were detected. Percentage of polymorphic bands (PPB), gene diversity (H_S and H_T) and Shannon information measure (Sh) were calculated both at population (PPB = 30.70%, H_S = 0.1105, Sh = 0.1646) and at species level (PPB = 68.42%, H_T = 0.2066, Sh = 0.3139).The existence of a spatial distribution of genetic diversity in R. persicifolia was revealed by a low gene flow, a fairly high level of genetic differentiation (G_ST = 0.4583) among populations and a positive correlation between genetic and geographic distances (Mantel test, r = 0.71, p = 0.016). The spatial genetic structure was also confirmed with BAPS analysis. Our results show that a certain level of isolation by distance and sex-ratio bias may explain the distribution of genetic diversity among populations.Conservation measures are suggested on the basis of the genetic diversity detected, by implementing an integrated in situ and ex situ conservation program for each population, in order to ensure effective protection for this endemic species.     

Genetic diversity and population structure of Swertia tetraptera (Gentianaceae), an endemic species of Qinghai-Tibetan Plateau

Swertia tetraptera Maxim is an annual alpine herb endemic to the Qinghai-Tibetan Plateau (QTP). Its populations are locally scattered as isolated patches throughout this region. Genetic variation within and among thirty-four populations of this species was assessed using ISSR fingerprinting with 10 primers. High levels of genetic diversity exist within species (P = 98.9%, I = 0.3475; He = 0.2227), while the within-population diversity is low (P = 32.7%, I = 0.177; He = 0.12). High levels of genetic differentiation were detected among populations based on various statistics, including Nei’s genetic diversity analysis (G_ST = 0.4608), Bayesian analysis (θ^B = 0.476) and AMOVA (F_ST = 0.57). That is, populations shared low levels of genetic identity (I = 0.2622–0.0966). This genetic structure was probably due to severe genetic drift, breeding system and limited gene flow. The observed genetic structure of the populations implies that different populations across the distribution range of the species should be sampled to maintain high genetic diversity when a conservation strategy is implemented.     

Highly differentiated populations of the narrow endemic and endangered species Primula cicutariifolia in China,revealed by ISSR and SSR

The perennial herb Primula cicutariifolia Pax is an endangered and endemic species with narrow distribution in eastern Anhui and Zhejiang provinces of China. In this study, the levels of genetic variation and the pattern of genetic structure in five natural populations of P. cicutariifolia were assessed by inter-simple sequence repeats (ISSR) and simple sequence repeat (SSR) markers. Both markers revealed that there was remarkably low genetic variation within populations (e.g., H_e = 0.19 and 0.18, for ISSR and SSR respectively) and high differentiation among populations (G_ST = 0.714 and 0.611; Φ_ST = 0.698 and 0.599, for ISSR and SSR respectively). The level of population genetic diversity was correlated to population size only detected by ISSR markers. The genetic structure of P. cicutariifolia may be explained by limited gene flow that was caused by habitat fragmentation and limited seeds and pollen dispersal ability, self breeding system and biennial life form. To protect and avoid disappearance of P. cicutariifolia, much more attentions should be paid to protect all the remnant populations and their habitats, and three management units, i.e. Tianmushan, Damodao, and Panan units, were proposed.     

Conservation genetics of Heritiera littoralis (Sterculiaceae), a threatened mangrove in China,based on AFLP and ISSR markers

AFLP and ISSR markers were used to determine the genetic variations in eight mangrove and non-mangrove populations of Heritiera littoralis (Sterculiaceae), a threatened species in China. Our results showed a moderate to high level of genetic variation in this species (P = 63.69%, H_T = 0.20 for AFLP; P = 76.07%, H_T = 0.22 for ISSR), and a relatively high level of genetic differentiation among populations (G_ST = 0.24 for AFLP and 0.27 for ISSR). Life history traits, long-distance dispersal of floating seeds, and local environments may play important roles in shaping the genetic diversity and genetic structure of this species. Investigation of the plant’s reproductive capacity showed that the natural seed production of H. littoralis was low, as was the germination rate and the transformation rate from juvenile to adult. H. littoralis is seriously threatened and is in urgent need of conservation in China.     

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.     

Effect of Yangtze River on population genetic structure of the relict plant Parrotia subaequalis in eastern China

Parrotia subaequalis (Hamamelidaceae) is a Tertiary relic species endemic in eastern China. We used inter‐simple sequence repeat (ISSR) markers to access genetic diversity and population genetic structure in natural five populations of P. subaequalis. The levels of genetic diversity were higher at species level (H = 0.2031) but lower at population level (H = 0.1096). The higher genetic diversity at species levels might be attributed to the accumulation of distinctive genotypes which adapted to the different habitats after Quaternary glaciations. Meanwhile, founder effects on the early stage, and subsequent bottleneck of population regeneration due to its biological characteristics, environmental features, and human activities, seemed to explain the low population levels of genetic diversity. The hierarchical AMOVA revealed high levels (42.60%) of among‐population genetic differentiation, which was in congruence with the high levels of Nei's genetic differentiation index (G_ST = 0.4629) and limited gene flow (N_m = 0.5801) among the studied populations. Mantel test showed a significant isolation‐by‐distance, indicating that geographic isolation has a significant effect on genetic structure in this species. Unweighted pair‐group method with arithmetic average clustering, PCoA, and Bayesian analyses uniformly recovered groups that matched the geographical distribution of this species. In particular, our results suggest that Yangtze River has served as a natural barrier to gene flow between populations occurred on both riversides. Concerning the management of P. subaequalis, the high genetic differentiation among populations indicates that preserving all five natural populations in situ and collecting enough individuals from these populations for ex situ conservation are necessary.     

Assessment of variability among morphological and molecular characters in wild populations of mint [Mentha longifolia (L.) L.] germplasm

Mentha longifolia is an important medicinal and aromatic perennial herb that exhibits wide distribution range from sub-tropical to temperate regions. In the present study, agro-morphological traits and genetic differences in 19 different populations of M. longifolia were studied to evaluate the level and extent of its diversity. Analysis of variance (ANOVA) showed that the different phenotypic characters show considerable differences among various populations and was significant at p < 0.05. Molecular diversity analysis performed by using arbitrary amplified eleven ISSR primers generated a total of 121 amplicons that range within the size of 200–2500 base pairs (bp). Each primer on average generated 11 amplicons with percentage polymorphism being 100. The analysis of molecular variance (AMOVA) showed more (64%) among population genetic diversity and less (36%) within the populations. Greater genetic differentiation (Gst = 0.6852) among these populations occurs due to low gene flow (Nm = 0.2297) and greater habitat variability. Geographic and genetic distances were positively correlated according to Mantel’s test. In order to remove any kind of biases, we used R software to perform cluster and redundancy analysis to analyse the extent of relatedness among studied populations. In terms of morphological and molecular aspects, the populations were grouped into four and five clusters respectively based on hierarchical clustering method. The results demonstrated that M. longifolia displays a great degree of morphological and genetic variation and can be utilized in breeding, genetic improvement, and gene bank conservation programmes in future.     

Genetic diversity and population structure of the endemic Azorean juniper,Juniperus brevifolia (Seub.) Antoine,inferred from SSRs and ISSR markers

Juniperus brevifolia is an important woody species endemic of the Azores archipelago (Portugal), found from coastal to mountain environments. Due to colonization and grazing pressure this species has suffered fragmentation, leading to extinction in one island and being threatened in others. The genetic diversity and population structure of J. brevifolia populations was studied to provide guidelines for restoration and conservation programmes. Nuclear Single Sequence Repeats (nSSR) from Juniperus communis and Juniperus przewalskii and Inter Simple Sequence Repeats (ISSR) were tested and results compared to test the transferability of the microsatellites to J. brevifolia. Samples from ten populations over three islands, divided as coastal, mid-altitude and mountain were analysed. Both marker systems revealed results statistically and strongly correlated with each other, and not dependent on population sample size. We observed positive fixation indexes, moderate to high levels of genetic diversity (h = 0.415 for nSSR and h = 0.245 for ISSR), low to moderate φ_pt genetic differentiation (0.070 for nSSR and 0.129 for ISSR) and high gene flow (N_m > 2.432). Regarding Nei's genetic distance the coastal communities clustered together reflecting the phenotypic plasticity, but no specific clustering was observed regarding φ_pt values. Therefore no populations with substantial genetic differentiation were identified, once the diversity is mostly observed within populations. However it is advised the continuous monitoring of J. brevifolia.     

Genetic diversity of the threatened Chinese endemic plant,Sinowilsonia henryi Hemsi. (Hamamelidaceae), revealed by inter-simple sequence repeat (ISSR) markers

Sinowilsonia henryi Hemsi., the only representative of the monotypic genus Sinowilsonia Hemsi. (Hamamelidaceae), is a threatened plant endemic to China with high phylogenetical, ecological and economical values. In the present study, inter-simple sequence repeat (ISSR) markers were employed to investigate the genetic diversity and differentiation of 214 individuals sampled from 14 populations. Fifteen selected primers yielded a total of 178 bright and discernible bands. The genetic diversity was low at the population level (h = 0.1025; I = 0.1506; PPL = 26.7%), but quite high at the species level (h = 0.2449; I = 0.3690; PPL = 72.5%). In line with the limited gene flow (N_m = 0.3537), the hierarchical analysis of molecular variance (AMOVA) revealed pronounced genetic differentiation among populations (Φ_ST = 0.6639). Furthermore, the Mantel test revealed a significant correlation between genetic and geographic distances among populations (r = 0.688, P = 0.001), indicating the role of geographic isolation in shaping its present population genetic structure. The present patterns of genetic diversity of S. henryi were assumed to result largely from its evolutionary history and geographic factors. Based on these findings, conservation strategies were proposed to preserve this threatened plant.     

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 variation within the sand-fixation species Caragana microphylla (Leguminosae) in Horqin sandy land detected by inter-simple sequence repeats analysis

Caragana microphylla is the most dominant and constructive shrub species in the Horqin sandy land in the northeast of China. We evaluated it's the level of genetic variation within and among populations sampled from two different populations types in Horqin sandy land by using inter-simple sequence repeat polymorphism (ISSR) molecular markers. The results showed that eight ISSR primers generated 106 bands, of which 87 were polymorphic. At the species level, genetic diversity was relatively high (P = 82.08%, h = 0.2831, I = 0.4233). Genetic variation in natural populations (h = 0.2152, I = 0.3169) was more than that in plantation populations (h = 0.2021, I = 0.3040). Based on Nei's G_ST value, more genetic differentiation among plantation populations was detected (G_ST = 0.7787). Six populations of C. microphylla clustered into two clades. These results have important implications for restoring and managing the degraded ecosystem in arid and semi-arid areas.     

Genetic,Epigenetic, and HPLC Fingerprint Differentiation between Natural and Ex Situ Populations of Rhodiola sachalinensis from Changbai Mountain,China

Rhodiola sachalinensis is an endangered species with important medicinal value. We used inter-simple sequence repeat (ISSR) and methylation-sensitive amplified polymorphism (MSAP) markers to analyze genetic and epigenetic differentiation in different populations of R. sachalinensis, including three natural populations and an ex situ population. Chromatographic fingerprint was used to reveal HPLC fingerprint differentiation. According to our results, the ex situ population of R. sachalinensis has higher level genetic diversity and greater HPLC fingerprint variation than natural populations, but shows lower epigenetic diversity. Most genetic variation (54.88%) was found to be distributed within populations, and epigenetic variation was primarily distributed among populations (63.87%). UPGMA cluster analysis of ISSR and MSAP data showed identical results, with individuals from each given population grouping together. The results of UPGMA cluster analysis of HPLC fingerprint patterns was significantly different from results obtained from ISSR and MSAP data. Correlation analysis revealed close relationships among altitude, genetic structure, epigenetic structure, and HPLC fingerprint patterns (R² = 0.98 for genetic and epigenetic distance; R² = 0.90 for DNA methylation level and altitude; R² = –0.95 for HPLC fingerprint and altitude). Taken together, our results indicate that ex situ population of R. sachalinensis show significantly different genetic and epigenetic population structures and HPLC fingerprint patterns. Along with other potential explanations, these findings suggest that the ex situ environmental factors caused by different altitude play an important role in keeping hereditary characteristic of R. sachalinensis.     

Genetic Diversity and Structure of Sinopodophyllum hexandrum (Royle) Ying in the Qinling Mountains,China

Sinopodophyllum hexandrum is an important medicinal plant whose genetic diversity must be conserved because it is endangered. The Qinling Mts. are a S. hexandrum distribution area that has unique environmental features that highly affect the evolution of the species. To provide the reference data for evolutionary and conservation studies, the genetic diversity and population structure of S. hexandrum in its overall natural distribution areas in the Qinling Mts. were investigated through inter-simple sequence repeats analysis of 32 natural populations. The 11 selected primers generated a total of 135 polymorphic bands. S. hexandrum genetic diversity was low within populations (average H_e = 0.0621), but higher at the species level (H_e = 0.1434). Clear structure and high genetic differentiation among populations were detected by using the unweighted pair group method for arithmetic averages, principle coordinate analysis and Bayesian clustering. The clustering approaches supported a division of the 32 populations into three major groups, for which analysis of molecular variance confirmed significant variation (63.27%) among populations. The genetic differentiation may have been attributed to the limited gene flow (N_m = 0.3587) in the species. Isolation by distance among populations was determined by comparing genetic distance versus geographic distance by using the Mantel test. Result was insignificant (r = 0.212, P = 0.287) at 0.05, showing that their spatial pattern and geographic locations are not correlated. Given the low within-population genetic diversity, high differentiation among populations and the increasing anthropogenic pressure on the species, in situ conservation measures were recommended to preserve S. hexandrum in Qinling Mts., and other populations must be sampled to retain as much genetic diversity of the species to achieve ex situ preservation as a supplement to in situ conservation.     

Genetic diversity and structure of Pinus dabeshanensis revealed by expressed sequence tag-simple sequence repeat (EST-SSR) markers

Assessing patterns of genetic variation in rare endangered species is critical for developing both in situ and ex situ conservation strategies. Pinus dabeshanensis Cheng et Law is an endangered species endemic to the Dabieshan Mountains of eastern China. To obtain fundamental information of genetic diversity, population history, effective population size, and gene flow in this species, we explored patterns of genetic variation of natural populations, in addition to an ex situ conserved population, using expressed sequence tag-simple sequence repeats (EST-SSR) markers. Our results revealed moderate levels of genetic diversity (e.g., H_E = 0.458 vs. H_E = 0.423) and a low level of genetic differentiation (F_ST = 0.028) among natural and conserved populations relative to other conifers. Both contemporary and historical migration rates among populations were high. Bayesian coalescent-based analyses suggested that 3 populations underwent reductions in population size ca. 10,000 yr ago, and that two populations may have experienced recent genetic bottlenecks under the TPM. Bayesian clustering revealed that individuals from the ex situ population were largely assigned to the ‘red’ cluster. Additionally, our results identified private alleles in the natural populations but not in the ex situ population, suggesting that the ex situ conserved population insufficiently represents the genetic diversity present in the species. Past decline in population size is likely to be due to Holocene climate change. Based on the genetic information obtained for P. dabeshanensis, we propose some suggestions for the conservation and efficient management of this endangered species.     

Fragmented habitat drives significant genetic divergence in the Chinese endemic plant,Urophysa henryi (Ranuculaceae)

Urophysa henryi (Oliv.) Ulbr., endemic to China with small populations, is known as a medicinal plant. In this study, ISSR markers were used to assess the genetic diversity and population structure throughout its entire distribution areas. Twelve primers revealed high genetic diversity at the species level (PPB = 95.6%; H = 0.3441; I = 0.5111), as well as high level of genetic differentiation (F_ST = 0.659, p < 0.001; G_ST = 0.677) and restricted gene flow (Nm = 0.239) among populations. According to the UPGMA and PCoA analysis, the 9 populations were clustered into three main groups, which were roughly in accordance with their geographical regions. In addition, a significant correlation between the genetic difference and geographic distances among populations was detected from the IBD analysis (r = 0.516, p = 0.003). These results indicated that the habitat heterogeneity and physical barriers play important roles in the modern distribution pattern and population divergence of U. henryi. However, human activities have posed serious threat to its living environment and continued survival. It is necessary to adopt some measures to restrict anthropogenic disturbances and preserve the existing populations.     

Population genetic structure and conservation implications of Ceriops decandra in Malay Peninsula and North Australia

Comprehensive information of mangrove genetic resources is requisite for developing strategies for their effective conservation and sustainable use. Genetic diversity within and among populations of a widespread viviparous mangrove Ceriops decandra was determined using inter-simple sequence repeat (ISSR). Ten natural populations were collected from Malay Peninsula and North Australia. At the species level, high genetic variation was detected (P = 72%, H_E = 0.253, and I = 0.379). The estimate of G_ST was 0.882, indicating a high level of genetic differentiation among populations. When populations were grouped according to geographic regions, i.e., East Malaya, West Malaya, Southmost Malaya, and North Australia, AMOVA suggested that most of the total variation (87%) was accounted for by differentiation between regions, with only 4% accounting for variation among populations within regions, and a further 9% partitioned among individuals within a population. A UPGMA dendrogram based on genetic distance revealed a deep split between populations from the eastern Indian Ocean and all others from the western Pacific Ocean, which may result from the historical lowering of sea level at these regions during the recent Pleistocene glaciations. An understanding of the genetic structure of C. decandra provides insight for the conservation and management of this species.     

Genetic diversity of endangered Manglietia patungensis assessed by inter simple sequence repeat and sequence-related amplified polymorphism markers

Manglietia patungensis Hu is an endangered plant native to China. Knowledge of its genetic diversity and structure would aid its conservation. This study assessed nine natural populations of M. patungensis using two methods: inter simple sequence repeat (ISSR) and sequence-related amplified polymorphism (SRAP) markers. Using 10 ISSR primer pairs, 334 bands were generated, and 10 SRAP primer pairs generated 276 bands. The percent of polymorphic bands (91.32% and 93.48%), Nei's genetic diversity (0.3448 and 0.3323), and Shannon's information index (0.5075 and 0.4935) revealed a high level of genetic diversity at the species level. Total heterozygosity was 0.3439 by ISSR and 0.3281 by SRAP. The mean heterozygosity was 0.2323 by ISSR and 0.2521 by SRAP. The coefficient of genetic differentiation among natural populations was 0.3245 by ISSR and 0.2316 by SRAP. These data indicated higher levels of genetic diversity of M. patungensis within, rather than among, populations. Estimates of gene flow among natural populations were 1.0411 and 1.0589, which implied a certain amount of gene exchange among populations. A Mantel test revealed no significant correlation between genetic and geographic distance. ISSR and SRAP markers are both effective for genetic diversity research in M. Patungensis. Based on these results, conservation of M. patungensis should be performed both in situ and ex situ.     

Low genetic diversity and significant structuring in the endangered Mentha cervina populations and its implications for conservation

Eighteen populations of the endangered aromatic and medicinal plant Mentha cervina (Lamiaceae) were sampled across its natural range, in the western half of the Iberian Peninsula, and inter-simple sequence repeats (ISSRs) markers were used to assess genetic diversity and population structure. M. cervina populations exhibited a relatively low genetic diversity (percentage of polymorphic loci PPB = 14.2–58.3%, Nei's genetic diversity H_e = 0.135–0.205, Shannon's information index I = 0.08 − 0.33). However, the genetic diversity at species level was relatively high (PPB = 98.3%; H_e = 0.325; I = 0.23). The results of the analysis of molecular variance indicated very structured populations, with 50% of the variance within populations, 44% among populations and 6% between regions defined by hydrographic basins, in line with the gene differentiation coefficient (G_ST = 0.532). A Mantel test did not find significant correlation between genetic and geographic distance matrices (r = 0.064), indicating that isolation by distance is not shaping the present genetic structure. The levels and patterns of genetic diversity in M. cervina populations were assumed to result largely from a combination of evolutionary history and its unique biological traits, such as breeding system, low capacity of dispersion, small effective size and habitat fragmentation. The high genetic differentiation among populations indicates the necessity of conserving the maximum possible number of populations. The results also provide information to select sites for ex situ conservation. Optimal harvesting strategies, cultivation and tissue culture should also be developed as soon as possible to guarantee sustainable use of the species under study.     

Genetic diversity and structure of the noxious alien grass Praxelis clematidea in southern China

Praxelis clematidea (Asteraceae), a plant species native to South America, is a noxious weed in southern China. We examined the genetic variation and population structure of 12 populations (76 individuals) of P. clematidea from Fujian, Guangdong, and Hainan Provinces in China using inter-simple sequence repeat (ISSR) analysis. From an initial set of 69 ISSR primers, 10 were selected which yielded 80 reproducible bands. Polymorphic bands (P) were 100%, Shannon's information index (I) was 0.4226, and Nei's gene diversity (H) was 0.2791. We infer that the high levels of genetic diversity exhibited by P. clematidea may have contributed to its invasiveness. Gene flow among populations was 2.4930, which has led to homogenization. The coefficient of population differentiation (Gst = 0.1671) indicated low levels of genetic variation among populations and high levels of genetic polymorphism within populations. There was a negative correlation between population elevation and genetic diversity, while there was a significant positive correlation between genetic distance and geographic distance based on a Mantel test (r = 0.5820, P < 0.01). Some populations from different provinces clustered together in principal coordinate and UPGMA analyses indicating that human-mediated events may have contributed to the dispersal of the species.     

Genetic diversity within and among populations of Neopicrorhiza scrophulariiflora (Scrophulariaceae) in China,an endangered medicinal plant

Neopicrorhiza scrophulariiflora Pennel (Hong) (Scrophulariaceae) is an endangered medicinal plant, endemic to the Eastern Himalayas and the Hengduan Mountains. Levels of genetic variation and genetic structure of seven populations of N. scrophulariiflora in China were studied using inter-simple sequence repeat (ISSR) markers. Thirteen primers amplified 82 total loci from 7 populations composed of 136 individuals. All 82 loci were polymorphic, showing a percentage of polymorphic bands (PPB) of 100%, indicating considerable genetic variation at the species level. In contrast, a low level of genetic diversity within populations was detected with a mean PPB of 30.56%. In addition, Nei’s genetic diversity analysis (0.4073) and Shannon’s diversity index (0.5917) revealed similar genetic structure. High levels of genetic differentiation (Gst = 0.6955) and restricted gene flow (Nm = 0.2198) among populations were also detected. Anthropologic impacts, together with clonal propagation, genetic drift and geographical isolation might be the reasons which had shaped the genetic structure of this species.