High levels of genetic differentiation of Oryza officinalis Wall. ex Watt. from China

In order to determine the population genetic structure of wild rice (Oryza officinalis Wall. ex Watt.), an endangered tropical and subtropical species, allozyme diversity encoded by 24 loci was analyzed electrophoretically in 145 individuals of eight natural populations from Hainan, Guangxi, and Yunnan provinces, China. A fairly high genetic differentiation (F(ST) = 0.882 and mean I = 0.786) was found among the studied populations. Our results suggest that restricted gene flow may play a significant role in shaping such a population genetic structure. In addition, high genetic differentiation among populations within a geographically limited region may stem from a reduced population size and consequent genetic drift.     

Decoupled mitochondrial and chloroplast DNA population structure reveals Holocene collapse and population isolation in a threatened Mexican-endemic conifer

Chihuahua spruce (Picea chihuahuana Martínez) is a montane subtropical conifer endemic to the Sierra Madre Occidental in northwestern México. Range-wide variation was investigated using maternally inherited mitochondrial (mtDNA) and paternally inherited chloroplast (cpDNA) DNA markers. Among the 16 mtDNA regions analysed, only two mitotypes were detected, while the study of six cpDNA microsatellite markers revealed eight different chlorotypes. The average cpDNA diversity (H = 0.415) was low but much higher than that for mtDNA (H = 0). The distribution of mitotypes revealed two clear nonoverlapping areas (G(ST) = N(ST) = 1), one including northern populations and the second one including the southern and central stands, suggesting that these two regions may represent different ancestral populations. The cpDNA markers showed lower population differentiation (G(ST) = 0.362; R(ST) = 0.230), implying that the two ancestral populations continued to exchange pollen after their initial geographic separation. A lack of a phylogeographic structure was revealed by different spatial analyses of cpDNA (G(ST) > R(ST); and samova), and reduced cpDNA gene flow was noted among populations (Nm = 0.873). Some stands deviated significantly from the mutation-drift equilibrium, suggesting recent bottlenecks. Altogether, these various trends are consistent with the hypothesis of a population collapse during the Holocene warming and suggest that most of the modern P. chihuahuana populations are now effectively isolated with their genetic diversity essentially modelled by genetic drift. The conservation efforts should focus on most southern populations and on the northern and central stands exhibiting high levels of genetic diversity. Additional mtDNA sequence analysis confirmed that P. martinezii (Patterson) is not conspecific with P. chihuahuana, and thus deserves separate conservation efforts.     

The postglacial history of Scots pine (Pinus sylvestris L.) in western Europe: evidence from mitochondrial DNA variation

The geographical structure of mitochondrial (mt)DNA variants (mitotypes) was investigated in 38 western European populations of Scots pine Pinus sylvestris using restriction fragment length polymorphism (RFLP) analysis of total DNA and a homologous cox1 probe. Three major mitotypes (designated a, b and d ) were detected. Within Spain all three major mitotypes were found, gene diversity was high, H_T = 0.586, and this diversity was distributed predominantly among rather than within populations (F_ST(M) = 0.813 for the seven Spanish populations). Mitotype d was present only in the most southerly population from the Sierra Nevada . Elsewhere in Europe, populations showed little or no mtDNA diversity within regions, but there were marked differences between regions. Italian populations were fixed for mitotype b ; populations from northern France, Germany, Poland, Russia and southern Sweden were fixed for mitotype a ; while populations in northern Fennoscandia were fixed for mitotype b . The isolated Scottish populations were predominantly of mitotype a , but mitotype b was present in three of the 20 populations scored. In Scotland, UK gene diversity (H_T = 0.120) and genetic differentiation among populations (F_ST(M) = 0.37) was much lower than in Spain. When interpreted in the light of complementary data from pollen analysis and nuclear genetic markers, the results suggest that present-day populations of P. sylvestris in western Europe have been derived from at least three different sources after glaciation.     

High genetic diversity and population differentiation in Boechera fecunda, a rare relative of Arabidopsis

Conservation of endangered species becomes a critical issue with the increasing rates of extinction. In this study, we use 13 microsatellite loci and 27 single-copy nuclear loci to investigate the population genetics of Boechera fecunda, a rare relative of Arabidopsis thaliana, known from only 21 populations in Montana. We investigated levels of genetic diversity and population structure in comparison to its widespread congener, Boechera stricta, which shares similar life history and mating system. Despite its rarity, B. fecunda had levels of genetic diversity similar to B. stricta for both microsatellites and nucleotide polymorphism. Populations of B. fecunda are highly differentiated, with a majority of genetic diversity existing among populations (F(ST) = 0.57). Differences in molecular diversity and allele frequencies between western and eastern population groups suggest they experienced very different evolutionary histories.     

Genetic variation of Trigonobalanus verticillata,a primitive species of Fagaceae,in Malaysia revealed by chloroplast sequences and AFLP markers

The genetic variation of Trigonobalanus verticillata, the most recently described genus of Fagaceae, was studied using chloroplast DNA sequences and AFLP fingerprinting. This species has a restricted distribution that is known to include seven localities in tropical lower montane forests in Malaysia and Indonesia. A total of 75 individuals were collected from Bario, Kinabalu, and Fraser's Hill in Malaysia. The sequences of rbcL, matK, and three non-coding regions (atpB-rbcL spacer, trnL intron, and trnL-trnF spacer) were determined for 19 individuals from these populations. We found a total of 30 nucleotide substitutions and four length variations, which allowed identification of three haplotypes characterizing each population. No substitutions were detected within populations, while the tandem repeats in the trnL -trnF spacer had a variable repeat number of a 20-bp motif only in Kinabalu. The differentiation of the populations inferred from the cpDNA molecular clock calibrated with paleontological data was estimated to be 8.3 MYA between Bario and Kinabalu, and 16.7 MYA between Fraser's Hill and the other populations. In AFLP analysis, four selective primer pairs yielded a total of 431 loci, of which 340 (78.9%) were polymorphic. The results showed relatively high gene diversity (H(S) = 0.153 and H(T) = 0.198) and nucleotide diversity (pi(S) = 0.0132 and pi(T) = 0.0168) both within and among the populations. Although the cpDNA data suggest that little or no gene flow occurred between the populations via seeds, the fixation index estimated from AFLP data (F(ST) = 0.153 and N(ST) = 0.214) implies that some gene flow occurs between populations, possibly through pollen transfer.     

Population structure of the African savannah elephant inferred from mitochondrial control region sequences and nuclear microsatellite loci

Two hundred and thirty-six mitochondrial DNA nucleotide sequences were used in combination with polymorphism at four nuclear microsatellite loci to assess the amount and distribution of genetic variation within and between African savannah elephants. They were sampled from 11 localities in eastern, western and southern Africa. In the total sample, 43 haplotypes were identified and an overall nucleotide diversity of 2.0% was observed. High levels of polymorphism were also observed at the microsatellite loci both at the level of number of alleles and gene diversity. Nine to 14 alleles per locus across populations and 44 alleles in the total sample were found. The gene diversity ranged from 0.51 to 0.72 in the localities studied. An analysis of molecular variance showed significant genetic differentiation between populations within regions and also between regions. The extent of subdivision between populations at the mtDNA control region was approximately twice as high as shown by the microsatellite loci (mtDNA F(ST) = 0.59; microsatellite R(ST) = 0.31). We discuss our results in the light of Pleistocene refugia and attribute the observed pattern to population divergence in allopatry accompanied by a recent population admixture following a recent population expansion.     

Microsatellite analysis of genetic diversity and population genetic structure of a wild rice (<Emphasis Type="Italic">Oryza rufipogon</Emphasis> Griff.) in China

Genetic diversity and population genetic structure of natural Oryza rufipogon populations in China were studied based on ten microsatellite loci. For a total of 237 individuals of 12 populations collected from four regions, a moderate to high level of genetic diversity was observed at population levels with the number of alleles per locus ( A) ranging from 2 to 18 (average 10.6), and polymorphic loci ( P) from 40.0% to 100% (average 83.3%). The observed heterozygosity ( H(O)) varied from 0.163 to 0.550 with the mean of 0.332, and the expected heterozygosity ( H(E)) from 0.164 to 0.648 with the mean of 0.413. The level of genetic diversity for Guangxi was the highest. These results are in good agreement with previous allozyme and RAPD studies. However, it was unexpected that high genetic differentiation among populations was found ( R(ST) = 0.5199, theta = 0.491), suggesting that about one-half of the genetic variation existed between the populations. Differentiation (pairwise theta) was positively correlated with geographical distance ( r = 0.464), as expected under the isolation by distance model. The habitat destruction and degradation throughout the geographic range of O. rufipogon may be the main factor attributed to high genetic differentiation among populations of O. rufipogon in China.     

New evidence from mitochondrial DNA of a progenitor-derivative species relationship between black spruce and red spruce (Pinaceae)

Mitochondrial DNA (mtDNA) markers were used to assess the genetic diversity in allopatric populations of black spruce (Picea mariana [Mill.] BSP) and red spruce (P. rubens Sarg.). Patterns of mitochondrial haplotypes (mitotypes) were strikingly different between the two species. All mtDNA markers surveyed were polymorphic in black spruce, revealing four different mitotypes and high levels of mtDNA diversity (P(p) = 100%, A = 2.0, H = 0.496). In contrast, populations of red spruce had only two mitotypes and harbored low levels of ggenetic diversity (P(p) = 13.2%, A = 1.1, H = 0.120). When the southernmost allopatric populations of red spruce were considered, only one mitotype was detected. As previously reported for nuclear gene loci, the diversity observed for mtDNA in red spruce was a subset of that found in black spruce. Comparison of present and previously published data supports the hypothesis of a recent progenitor-derivative relationship between these species, red spruce presumably being derived by allopatric speciation of an isolated population of black spruce during the Pleistocene.     

Levels and patterns of genetic variation in the endangered species Abronia macrocarpa

Genetic variation was evaluated in the federally endangered species Abronia macrocarpa (large-fruited sand-verbena), an herbaceous perennial restricted to deep sandy soils and endemic to three counties of east-central Texas. Seven of the ten known populations were sampled and analyzed using starch gel electrophoresis of eight enzymes coded by 18 interpretable loci. Duplicate gene expression was observed for four loci, suggesting polyploid ancestry for the lineage that includes A. macrocarpa. Values for estimators of genetic polymorphism within populations (ranges: P = 38.9%-61.1%, A = 1.7-2.1, H = 0.122-0.279) exceeded average values for seed plants (P = 34.2%, A = 1.53, H = 0.113). Genotype proportions at most loci in most populations were in Hardy-Weinberg equilibrium, consistent with obligate outcrossing previously documented for this species; exceptions could be attributed to population substructure. Values of F(ST) tended to be high, ranging from 0.021 to 0.481 for individual loci (mean F(ST) = 0.272), indicating substantial divergence and limited gene flow among populations, despite their close geographic proximity. Pairwise values of Nei's genetic identity between populations ranged from 0.799 to 0.975 and tended to be influenced by geographic proximity of population pairs. Collectively, these data suggest a long history of isolation among populations that have not been subjected to bottlenecks. Isolation of A. macrocarpa populations apparently results from the disjunct occurrence of suitable habitat and perhaps has been accentuated by human disturbance.     

Selection of high heterozygosity popcorn varieties in Brazil based on SSR markers

We analyzed genetic structure and diversity among eight populations of popcorn, using SSR loci as genetic markers. Our objectives were to select SSR loci that could be used to estimate genetic diversity within popcorn populations, and to analyze the genetic structure of promising populations with high levels of heterozygosity that could be used in breeding programs. Fifty-seven alleles (3.7 alleles per locus) were detected; the highest effective number of alleles (4.21) and the highest gene diversity (0.763) were found for the Umc2226 locus. A very high level of population differentiation was found (F(ST) = 0.3664), with F(ST) for each locus ranging from 0.1029 (Umc1664) to 0.6010 (Umc2350). This analysis allowed us to identify SSR loci with high levels of heterozygosity and heterozygous varieties, which could be selected for production of inbred lines and for developing new cultivars.     

Phylogeography of Camellia taliensis (Theaceae) inferred from chloroplast and nuclear DNA: insights into evolutionary history and conservation

ABSTRACT: BACKGROUND: As one of the most important but seriously endangered wild relatives of the cultivated tea, Camellia taliensis harbors valuable gene resources for tea tree improvement in the future. The knowledge of genetic variation and population structure may provide insights into evolutionary history and germplasm conservation of the species. RESULTS: Here, we sampled 21 natural populations from the species' range in China and performed the phylogeography of C. taliensis by using the nuclear PAL gene fragment and chloroplast rpl32-trnL intergenic spacer. Levels of haplotype diversity and nucleotide diversity detected at rpl32-trnL (h = 0.841; pi = 0.00314) were almost as high as at PAL (h = 0.836; pi = 0.00417). Significant chloroplast DNA population subdivision was detected (GST = 0.988; NST = 0.989), suggesting fairly high genetic differentiation and low levels of recurrent gene flow through seeds among populations. Nested clade phylogeographic analysis of chlorotypes suggests that population genetic structure in C. taliensis has been affected by habitat fragmentation in the past. However, the detection of a moderate nrDNA population subdivision (GST = 0.222; NST = 0.301) provided the evidence of efficient pollen-mediated gene flow among populations and significant phylogeographical structure (NST > GST; P < 0.01). The analysis of PAL haplotypes indicates that phylogeographical pattern of nrDNA haplotypes might be caused by restricted gene flow with isolation by distance, which was also supported by Mantel's test of nrDNA haplotypes (r = 0.234, P < 0.001). We found that chlorotype C1 was fixed in seven populations of Lancang River Region, implying that the Lancang River might have provided a corridor for the long-distance dispersal of the species. CONCLUSIONS: We found that C. taliensis showed fairly high genetic differentiation resulting from restricted gene flow and habitat fragmentation. This phylogeographical study gives us deep insights into population structure of the species and conservation strategies for germplasm sampling and developing in situ conservation of natural populations.     

Population genetic structure of Arabidopsis lyrata in Europe

Population genetic theory predicts that the self-incompatible and perennial herb, Arabidopsis lyrata, will have a genetic structure that differs from the self-fertilizing, annual Arabidopsis thaliana. We quantified the genetic structure for eight populations of A. lyrata ssp. petraea in historically nonglaciated regions of central Europe. Analysis of 20 microsatellite loci for 344 individuals demonstrated that, in accordance with predictions, diploid populations had high genome-wide heterozygosity (H(O) = 0.48; H(E) = 0.52), high within-population diversity (83% of total) compatible with mutation-drift equilibrium, and moderate differentiation among populations (F(ST) = 0.17). Within a single population, the vast majority of genetic variability (92%) was found at the smallest spatial scale (< 3 m). Although there was no evidence of biparental inbreeding or clonal propagation at this scale (F(IS) = 0.003), significant fine-scale spatial autocorrelation indicated localized gene flow presumably due to gravity dispersed seeds (Sp = 0.018). Limited gene flow between isolated population clusters (regions) separated by hundreds of kilometres has given rise to an isolation by distance pattern of diversification, with low, but significant, differentiation among regions (F(ST) = 0.05). The maintenance of geographically widespread polymorphisms and uniformly high diversity throughout central Europe is consistent with periglacial survival of A. lyrata ssp. petraea north of the Alps in steppe-tundra habitats during the last glacial maximum. As expected of northern and previously glaciated localities, A. lyrata in Iceland was genetically less diverse and highly differentiated from central Europe (H(E) = 0.37; F(ST) = 0.27).     

Microsatellite and mitochondrial haplotype diversity reveals population differentiation in the tiger shrimp (Penaeus monodon) in the Indo-Pacific region

The black tiger shrimp (Penaeus monodon) is an ecologically and economically important penaeid species and is widely distributed in the Indo-Pacific region. Here we investigated the genetic diversity of P. monodon (n = 355) from eight geographical regions by genotyping at 10 microsatellite loci. The average observed heterozygosity at various loci ranged from 0.638 to 0.743, indicating a high level of genetic variability in this region. Significant departures from Hardy-Weinberg equilibrium caused by heterozygote deficiency were recorded for most loci and populations. Pairwise F(ST) and R(ST) values revealed genetic differentiation among the populations. Evidence from the assignment test showed that the populations in the West Indian Ocean were unique, whereas other populations examined were partially admixed. In addition, the non-metric multidimensional scaling analysis indicated the presence of three geographic groups in the Indo-Pacific region, i.e. the African populations, a population from western Thailand and the remaining populations as a whole. We also sequenced and analysed the mitochondrial control region (mtCR) in these shrimp stocks to determine whether the nuclear and mitochondrial genomes show a similar pattern of genetic differentiation. A total of 262 haplotypes were identified, and nucleotide divergence among haplotypes ranged from 0.2% to 16.3%. Haplotype diversity was high in all populations, with a range from 0.969 to 1. Phylogenetic analysis using the mtCR data revealed that the West Indian Ocean populations were genetically differentiated from the West Pacific populations, consistent with the microsatellite data. These results should have implications for aquaculture management and conservation of aquatic diversity.     

银杉花粉生命力及其变异

中国特有植物银杉(Cathay argyrophylla)是最濒危的松杉类植物之一。研究表明,受精前胚珠的高败育率是银杉生殖效率低的重要原因之一,但迄今对银杉花粉活力及其变异仍不得而知。由于花粉活力对授粉、受精、种子产量和质量以及后代适合度都有直接影响,该研究采用TTC染色和体外萌发法测定了来自大瑶山(Dayaoshan, DYS)和花坪(Huaping, HP)两个地区7个种群16个个体52份银杉花粉样品的生活力。结果表明银杉花粉的活力水平(93.3%)与其它裸子植物相当,干燥低温条件下银杉的花粉活力比较稳定,体外萌发的适宜蔗糖浓度在13%左右。ANOVA分析揭示种群内个体间的花粉活力差异不大(p>0.05),而地区内的种群间及地区间的花粉活力有显著差异(p<0.05)。银杉花粉活力与其生殖成功没有相关性。在上述基础上,从控制授粉的角度讨论了银杉的保护策略和复壮的可能性。     

濒危物种巴东木莲的等位酶遗传多样性及其保护策略

为了确定中国特有濒危植物巴东木莲(Manglietia patungensis)的就地保护优先单元和制定迁地保护的取样策略,采用超薄平板聚丙烯酰胺等电聚焦电泳技术对巴东木莲的7个野生居群的等位酶遗传多样性及遗传结构进行了分析。通过对8个酶系统19个酶位点的分析结果来看,巴东木莲具有较高的遗传多样性和一定程度的遗传分化。每位点平均等位基因数（A）为1.57,多态位点百分率（P）为48.1%,平均预期杂合度（H_e）为0.192。巴东木莲居群间的遗传分化系数（G_ST）为0.165,说明其16.5％的遗传变异存在于居群间;基于遗传分化系数计算的基因流（N_m）为1.27,说明居群间存在适度的基因流;固定指数(F)为-0.191,显示巴东木莲居群杂合体轻微过量,纯合体略显不足。本研究表明,目前巴东木莲仍具有较高的遗传多样性并且在其适宜生境中可以完成自然更新,因此应采取以就地保护为主的综合保育策略。     

Allozyme variation in the endangered insular endemic Castilleja grisea

BACKGROUND AND AIMS: Genetic diversity in Castilleja grisea, an endangered, perennial herb endemic to San Clemente Island, California was investigated. Subsequent to the elimination of goats from the island in 1992, many populations of C. grisea have reappeared and have been increasing in size. METHODS: Nineteen populations were surveyed for their genotype at 19 allozyme loci. KEY RESULTS: At the taxon level, 57.9 % of loci are polymorphic with A(P) = 3.09 and H(E) = 0.137. Populations averaged 33.0 % polymorphic loci with A(P) = 2.43 and H(E) = 0.099. Most variation is found within rather than among populations (G(ST) = 0.128), although differentiation among populations is significant. Genetic identities range from I = 0.960 to I = 1.000 with mean I = 0.990. There is no significant relationship between genetic and geographic distance. Gene flow among populations is Nm = 2.50 based on private alleles and Nm = 1.70 based on F(ST). Outcrossing rates based on fixation indices average t = 1.01, indicating a primarily out-crossed mating system. CONCLUSIONS: The observed genetic variation is moderately high, unusually so for an insular endemic species, suggesting that C. grisea may not have lost substantial genetic variation during 150 years of overgrazing, and indicating that it is unlikely to be endangered by genetic factors.     

Conservation genetics and population history of the threatened European mink Mustela lutreola, with an emphasis on the west European population

In species of great conservation concern, special attention must be paid to their phylogeography, in particular the origin of animals for captive breeding and reintroduction. The endangered European mink lives now in at least three well-separated populations in northeast, southeast and west Europe. Our aim is to assess the genetic structure of these populations to identify 'distinct population segments' (DPS) and advise captive breeding programmes. First, the mtDNA control region was completely sequenced in 176 minks and 10 polecats. The analysis revealed that the western population is characterized by a single mtDNA haplotype that is closely related to those in eastern regions but nevertheless, not found there to date. The northeast European animals are much more variable (pi = 0.012, h = 0.939), with the southeast samples intermediate (pi = 0.0012, h = 0.469). Second, 155 European mink were genotyped using six microsatellites. The latter display the same trends of genetic diversity among regions as mtDNA [gene diversity and allelic richness highest in northeast Europe (H(E) = 0.539, R(S) = 3.76), lowest in west Europe (H(E) = 0.379, R(S) = 2.12)], and provide evidences that the southeast and possibly the west populations have undergone a recent bottleneck. Our results indicate that the western population derives from a few animals which recently colonized this region, possibly after a human introduction. Microsatellite data also reveal that isolation by distance occurs in the western population, causing some inbreeding because related individuals mate. As genetic data indicate that the three populations have not undergone independent evolutionary histories for long (no phylogeographical structure), they should not be considered as distinct DPS. In conclusion, the captive breeding programme should use animals from different parts of the species' present distribution area.     

中国6个国家级保护鸭种群遗传多样性分析

为探讨国家级保护鸭种群的遗传多样性和遗传分化关系,利用17个微卫星标记,对我国6个国家级保护鸭品种资源的遗传多样性进行了分析,统计了平均遗传杂合度(H)、多态信息含量(PIC)、基因多样度(FST)、平均遗传分化系数(GST)和遗传距离等。结果表明,各鸭品种的杂合度较高,除两个群体表现为显著的杂合子缺失外,其他群体均处于Hard-Weinberg平衡状态。6个鸭种群间平均FST值为17.0%,平均遗传分化系数GST为14.7%,各品种平均杂合度为0.706～0.604,平均多态信息含量(PIC)为0.561～0.663。本研究说明我国6个国家级保护鸭品种资源各品种内和品种间的遗传变异大,遗传多样性丰富。     

Mitochondrial and chloroplast phylogeography of Picea crassifolia Kom. (Pinaceae) in the Qinghai-Tibetan Plateau and adjacent highlands

The disjunct distribution of forests in the Qinghai-Tibetan Plateau (QTP) and adjacent Helan Shan and Daqing Shan highlands provides an excellent model to examine vegetation shifts, glacial refugia and gene flow of key species in this complex landscape region in response to past climatic oscillations and human disturbance. In this study, we examined maternally inherited mitochondrial DNA (nad1 intron b/c and nad5 intron 1) and paternally inherited chloroplast DNA (trnC-trnD) sequence variation within a dominant forest species, Picea crassifolia Kom. We recovered nine mitotypes and two chlorotypes in a survey of 442 individuals from 32 populations sampled throughout the species' range. Significant mitochondrial DNA population subdivision was detected (G(ST) = 0.512; N(ST) = 0.679), suggesting low levels of recurrent gene flow through seeds among populations and significant phylogeographical structure (N(ST) > GST, P < 0.05). Plateau haplotypes differed in sequence from those in the adjacent highlands, suggesting a long period of allopatric fragmentation between the species in the two regions and the presence of independent refugia in each region during Quaternary glaciations. On the QTP platform, all but one of the disjunct populations surveyed were fixed for the same mitotype, while most populations at the plateau edge contained more than one haplotype with the mitotype that was fixed in plateau platform populations always present at high frequency. This distribution pattern suggests that present-day disjunct populations on the QTP platform experienced a common recolonization history. The same phylogeographical pattern, however, was not detected for paternally inherited chloroplast DNA haplotypes. Two chlorotypes were distributed throughout the range of the species with little geographical population differentiation (G(ST) = N(ST) = 0.093). This provides evidence for highly efficient pollen-mediated gene flow among isolated forest patches, both within and between the QTP and adjacent highland populations. A lack of isolation to pollen-mediated gene flow between forests on the QTP and adjacent highlands is surprising given that the Tengger Desert has been a geographical barrier between these two regions for approximately the last 1.8 million years.     

Genetic Structure Within and Among Populations of Saruma henryi, an Endangered Plant Endemic to China

The endangered perennial plant Saruma henryi Oliv. is endemic to China and has phylogenetic, ecological, and medicinal value. We used 10 microsatellite (SSR) loci to investigate genetic diversity and differentiation in 16 natural populations. Genetic diversity was high at the species level (H(E) = 0.9427, h = 0.9410, I = 3.0213) and low at the population level (H(E) = 0.4441, h = 0.4307, I = 0.6822). Pronounced genetic differentiation was detected among populations (G(ST) = 0.5428, F(ST) = 0.5524), in line with the limited among-population gene flow (Nm = 0.21). The significant isolation-by-distance pattern revealed by a Mantel test (r = 0.311, P = 0.001) suggested a clear geographic tendency in the distribution of genetic variability. Bayesian assignment and principal coordinates analyses supported the clustering of 16 populations into three groups. The present SSR results were also compared with previously published ISSR results. These results have significant implications for conservation of the species.