Genetic variability and population structure of Bergenia ciliata (Saxifragaceae) in the Western Himalaya inferred from DAMD and ISSR markers

Genetic variability and population structure of Bergenia ciliata (Haw.) Sternb., commonly known as “Pashanbheda” (Stone-breaker), collected from the Western Himalayan region of India were estimated using two DNA fingerprinting methods viz., directed amplification of minisatellite DNA (DAMD) and inter simple sequence repeats (ISSR). The cumulative data analysis of DAMD and ISSR markers for 74 accessions from eight populations showed 86.1% polymorphism. Analysis of molecular variance (AMOVA) showed highest percentage of variation within individuals of populations (73.6%) and 21.7% among populations. STRUCTURE and PCoA analyses on the hierarchical partitioning of genetic diversity showed strong admixture of individuals among the eight assumed geographical populations of B. ciliata. The data suggests that high genetic flow is one of the major factors responsible for low genetic differentiation. Preservation of genetic diversity of B. ciliata is important, both to promote adaptability of the populations to changing environment as well as to preserve a large gene pool for future prospection. The present study using DAMD and ISSR markers, therefore, provide the means of rapid characterization of accessions within the populations, and thus enable the selection of appropriate accessions for further utilization in conservation and prospection programmes.     

Estimation of genetic variability and population structure in Sapindus trifoliatus L., using DNA fingerprinting methods

Genetic variability and population structure of Sapindus trifoliatus L. (Sapindaceae), collected from Gujarat, Karnataka and Uttar Pradesh states were estimated using three DNA fingerprinting methods viz., random amplified polymorphic DNA (RAPD), directed amplification of minisatellite DNA (DAMD) and inter-simple sequence repeats (ISSR). The cumulative data analysis carried out for all three markers showed 69.42 % polymorphism. The intra-population genetic diversity analysis revealed the highest values of Nei’s genetic diversity (0.16), Shannon information index (0.24) and polymorphic loci (43.99 %) among Bhavnagar (BH) population, whereas lowest values were found in Junagarh (JU) population. The maximum inter-population average genetic distance (0.20) was between Allahabad (AL) and JU populations. Analysis of molecular variance (AMOVA) showed highest percentage of variation among individuals of populations (56 %) followed by 25 % among populations and 19 % among regions. Principal coordinate analysis and UPGMA dendrogram revealed that genetic diversity was in congruence with the geographical diversity. The data strongly suggest that low genetic flow, geographic isolation and to some extent genetic drift are the major factors responsible for high genetic differentiation. Preservation of genetic diversity of S. trifoliatus is important, both to promote adaptability of the populations to changing environment as well as to preserve a large gene pool for future genetic improvement. The present study using RAPD, DAMD and ISSR profiles of S. trifoliatus provide the means of rapid characterization of accessions within the populations, and thus enable the selection of appropriate accessions for further utilization in conservation and prospection programs of this important plant genetic resource.     

Genetic variability and population structure in Sapindus emarginatus Vahl from India

Sapindus emarginatus is an economically important tropical tree species sparsely distributed in different geographical provinces like Gangetic Plains, Western Ghats, and Deccan Plateau in India. In the present paper estimation of genetic variability within and among 41 accessions representing five populations was carried out using 3 single primer amplification reaction (SPAR) methods viz. RAPD, DAMD and ISSR. The cumulative data analysis was carried out for all three SPAR methods, and showed 82.32% polymorphism across all the accessions of S. emarginatus. Jaccard's similarity values among 41 accessions ranged from 0.15 to 0.49 with an average value of 0.37. The intra-population genetic diversity revealed highest values of Nei's genetic diversity (0.19,) Shannon information index (0.29) and polymorphic loci (55.18%), among the accessions of Gujarat (GJ) population, while the corresponding lowest values were (0.10), (0.15) and (26.40%) respectively among the accessions of Rajasthan (RJ) population. The maximum inter-population average genetic distance (0.20) was between Karnataka (KA) and RJ, while the corresponding least genetic distance (0.06) was between Allahabad (AL) and Varanasi (VS) populations. The analysis of molecular variance (AMOVA) revealed maximum percentage of variation among individuals of populations (72%) followed by 16% among regions and 12% among populations. Principal coordinate analysis (PCA) of cumulative data also supported the clustering pattern in the UPGMA dendrogram. These results suggest that genetic diversity is corroborating with the geographical diversity. Mantel's test was performed which revealed a highly significant correlation between cumulative vs RAPD, and showed the maximum (0.93) correlation coefficient, followed by cumulative vs ISSR (0.78) and cumulative vs DAMD (0.91) respectively, and this clearly indicates that the SPAR methods (RAPD, DAMD and ISSR) are sufficiently informative and are suitable to analyze the genetic variability within and among the populations of S. emarginatus.     

Assessment of genetic diversity and population structure of Bergenia stracheyi (Saxifragaceae) in the Western Himalaya (India)

Genetic diversity and population structure in Bergenia stracheyi, a threatened medicinal herb in the Western Himalaya of India was analysed using directed amplification of minisatellite DNA (DAMD) and inter simple sequence repeats (ISSR) markers. A total of 41 accessions of B. stracheyi representing three populations (Khillenmarg –KLM, Jalori Pass-JLP and Rohtang-RTG) were considered in the present study. The cumulative data analysis for 26 (10 DAMD + 16 ISSR) markers revealed 87.1% polymorphism. The maximum inter-population genetic distance was found between KLM and JLP, whereas the minimum genetic distance was found between RTG and JLP populations. The analysis of molecular variance (AMOVA) revealed maximum percentage of variation among individuals within populations (75%) than among the populations (25%). Clustering pattern of the three sample populations in STRUCTURE and PCoA analyses showed high genetic variation at population level. The present study revealed that distribution patterns, high altitudinal ranges, high habitat specificity, relatively high gene flow, small and isolated population size have shaped the current population structure of B. stracheyi in the Western Himalayan region. DAMD and ISSR markers have provided significant insights into characterization of B. stracheyi populations, and facilitate selection of appropriate accessions for further utilization in conservation and bioprospecting programmes.     

Population Genetic Study of Fagopyrum tataricum from Western Himalaya Using ISSR Markers

Inter-simple sequence repeat (ISSR) markers were used to analyze genetic diversity and relatedness of 15 germplasms of Fagopyrum tataricum. Samples representing 75 individuals were collected from a range of altitudes in the Western Himalaya. The 13 ISSR primers revealed 98.1% polymorphism among populations, whereas average polymorphism was extremely low (2.18%) within populations. The coefficient of population differentiation was 0.9750, with limited gene flow (N _m) of 0.0128. The average PIC value of the ISSR markers was high (0.812), with a marker ratio of 0.65 and marker index of 6.66. The genetic diversity of F. tataricum significantly correlated with altitude and gene diversity, Shannon’s index, and the percentage of polymorphic bands. The genetic diversity among populations showed broad genetic base and provided a developmental strategy for crop improvement.     

都支杜鹃遗传多样性的ISSR分析

通过ISSR分子标记对都支杜鹃(Rhododendron shanii)5个居群114个个体的遗传多样性进行了分析。8个引物共检测到69个位点,其中多态位点36个,多态位点百分率 （P） 为52.17%;物种水平的Shannon多样性指数（I）为0.2536,Nei指数（H）为0.1659;居群水平的P、I和H分别平均为32.46%、0.1955和0.1344;Nei’s基因分化系数GST=0.1845,基因流Nm=2.2104;遗传结构的AMOVA分析显示,总的遗传变异中88.3%存在于居群内,11.7%存在于居群间;PCA分析和Splite Tree构建的NJ树均显示居群间分化小;遗传距离与地理距离显著相关。较强的基因流、多年生木本及异交与分布区狭窄等可能是导致都支杜鹃遗传多样性较低的主要原因。目前该种的保护应以就地保护为主,尤其要注意对大居群的保护。     

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.     

Analysis of genetic variability and population structure of the endemic medicinal Limonium sinense using molecular markers

Limonium sinense is an endemic medicinal herb used to treat fever, hemorrhage and other disorders. In the present study, population genetic diversity was elucidated using random amplified polymorphic DNA (RAPD), inter simple sequence repeat (ISSR) and amplified fragment length polymorphism (AFLP) primers. Percentage of polymorphic bands, Nei's gene diversity and Shannon's information index revealed a high level of genetic diversity at species level. The analysis of molecular variance revealed that 69.88% (RAPD), 71.19% (ISSR) and 70.97% (AFLP) of variability were partitioned among individuals within populations, which indicated the coherent trend by Gst (0.3849/0.3577/0.3670). Gene flow number (Nm) was 0.581/0.618/0.612, which indicated that there was a limited gene exchange between populations. The UPGMA clustering results showed that the genetic distance had no significant correlation with geographic distance. These results indicate that these markers were reliable tools for the differentiation and determination of the genetic diversity among the populations of L. sinense and the conservation of existing natural population is necessary.     

Natural population structure and genetic differentiation for heterodicogamous plant: <Emphasis Type="Italic">Cyclocarya paliurus</Emphasis> (Batal.) Iljinskaja (Juglandaceae)

Cyclocarya paliurus, a Chinese monotypic genus belonging to the Juglandaceae, exhibits a widespread distribution in southern China. To provide resource support for its future breeding and medicinal utilization, much information involving natural population succession and genetic diversity of C. paliurus should be obtained. In this study, the structures of 30 natural populations throughout its distribution were analyzed based on the investigation, showing the ratio of individuals with DBH < 40 cm accounting for 71.7%. Inter-simple sequence repeat (ISSR) and simple sequence repeat (SSR) were used for evaluation of 244 individuals from 29 natural populations. High-level polymorphism (100%) indicates high efficiency of two markers. Analysis demonstrated that genetic diversity (value of 0.18 for ISSR and 0.09 for SSR) was low among populations. While gene flow among populations reached 1.182 and 1.335 for ISSR and SSR, respectively, indicating genetic variations mainly existing within population. Twenty-six populations (size ≥5 individuals) were grouped into five clusters based on the combined ISSR and SSR loci. The value of Nei’s coefficient varying from 0.87 to 0.96 indicates small differentiation among groups. Furthermore, we speculate that C. paliurus is an expanding species in subtropical China (most of populations are increasing type (IP)) and discuss the strategy of germplasm conservation.     

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 Relationships Among Wild and Cultivated Accessions of Curry Leaf Plant (Murraya koenigii (L.) Spreng.), as Revealed by DNA Fingerprinting Methods

Murraya koenigii (L.) Spreng. (Rutaceae), is an aromatic plant and much valued for its flavor, nutritive and medicinal properties. In this study, three DNA fingerprinting methods viz., random amplification of polymorphic DNA (RAPD), directed amplification of minisatellite DNA (DAMD), and inter-simple sequence repeat (ISSR), were used to unravel the genetic variability and relationships across 92 wild and cultivated M. koenigii accessions. A total of 310, 102, and 184, DNA fragments were amplified using 20 RAPD, 5 DAMD, and 13 ISSR primers, revealing 95.80, 96.07, and 96.73% polymorphism, respectively, across all accessions. The average polymorphic information content value obtained with RAPD, DAMD, and ISSR markers was 0.244, 0.250, and 0.281, respectively. The UPGMA tree, based on Jaccard’s similarity coefficient generated from the cumulative (RAPD, DAMD, and ISSR) band data showed two distinct clusters, clearly separating wild and cultivated accessions in the dendrogram. Percentage polymorphism, gene diversity (H), and Shannon information index (I) estimates were higher in cultivated accessions compared to wild accessions. The overall high level of polymorphism and varied range of genetic distances revealed a wide genetic base in M. koenigii accessions. The study suggests that RAPD, DAMD, and ISSR markers are highly useful to unravel the genetic variability in wild and cultivated accessions of M. koenigii.     

Assessment of genetic diversity through RAPD,ISSR and AFLP markers in Podophyllum hexandrum: a medicinal herb from the Northwestern Himalayan region

Total synthesis of podophyllotoxin is an expensive process and availability of the compound from the natural resources is an important issue for pharmaceutical companies that manufacture anticancer drugs. In order to facilitate reasoned scientific decisions on its management and conservation for selective breeding programme, genetic analysis of 28 populations was done with 19 random primers, 11 ISSR primers and 13 AFLP primer pairs. A total of 92.37 %, 83.82 % and 84.40 % genetic polymorphism among the populations of Podophyllum were detected using RAPD, ISSR and AFLP makers, respectively. Similarly the mean coefficient of gene differentiation (Gst) were 0.69, 0.63 and 0.51, indicating that 33.77 %, 29.44 % and 26 % of the genetic diversity resided within the population. Analysis of molecular variance (AMOVA) indicated that 53 %, 62 % and 64 % of the genetic diversity among the studied populations was attributed to geographical location while 47 %, 38 % and 36 % was attributed to differences in their habitats using RAPD, ISSR and AFLP markers. An overall value of mean estimated number of gene flow (Nm) were 0.110, 0.147 and 0.24 from RAPD, ISSR and AFLP markers indicating that there was limited gene flow among the sampled populations.     

珍稀濒危植物长叶红砂种群遗传多样性的ISSR分析

采用ISSR分子标记技术,对濒危小灌木长叶红砂(Reaumuria trigyna) 集中分布的5个种群的遗传多样性水平和遗传结构进行了研究。14条引物共检测到114个位点,其中99个为多态位点,多态位点比率为86.84％,长叶红砂种群具有较高的遗传多样性。物种水平上Shannon多样性指数(I)为0.468 8,Nei基因多样性指数(H)为0.308 4;种群水平上,多态位点比率P为77.89%,I为0.410 6,H为0.260 9,基因分化系数Gst为0.106 9,揭示了长叶红砂种群遗传变异多存在于种群内,种群间的遗传分化较小,占10.69%。 基因流(Nm)为4.178 7＞1,说明种群间的基因交流,防止了由于遗传漂变导致的遗传分化。聚类分析表明长叶红砂种群遗传距离与地理距离之间无显著的相关性。研究结果说明遗传多样性水平与物种本身特性和所处不同群落有关,濒危植物并不一定表现为遗传变异水平的降低。     

ISSR and DAMD markers revealed high genetic variability within Flavoparmelia caperata in Western Himalaya (India)

Flavoparmelia caperata (L.) Hale is medicinally very important and possesses antifungal and antibacterial activities. F. caperata is the only species found in India. Inter simple sequence repeat (ISSR) and Directed amplification of minisatellite DNA (DAMD) methods were used to analyze the genetic variability within F. caperata from the Western Himalayan region of India. Eleven ISSR and 10 DAMD primers produced 139 and 117 polymorphic bands, and detected 91.44 and 82.34 % polymorphisms, respectively. Cumulative band data generated for ISSR and DAMD markers resulted in 86.86 % polymorphism across all the accessions of F. caperata. The average Polymorphic information content (PIC) value obtained with ISSR, DAMD, and cumulative band data were 0.28, 0.27, and 0.27, respectively. The clustering of the F. caperata accessions in the UPGMA dendrogram showed that these accessions are intermingled with each other in different subclusters irrespective of their geographical affiliations. The pattern of genetic variations within F. caperata accessions could be due to free exchange of spores that might have taken place among these accessions in the wild. ISSR and DAMD markers efficiently and reliably resulted in discrete banding patterns and polymorphic profiles. These markers despite targeting different regions of genome, revealed almost similar levels of polymorphism across all the accessions. The wide range of genetic distance and high level of polymorphism detected by ISSR and DAMD reflected a high genetic variability among the different accessions of F. caperata.     

Single primer amplification reaction (SPAR) reveals intra-specific natural variation in Prosopis cineraria (L.) Druce

Prosopis cineraria, an important multipurpose tree and vital component of the otherwise fragile ecosystem of arid and semiarid regions of India. It is highly drought tolerant and sprouts profusely during the extreme dry summer months when most other trees are leafless. P. cineraria is known to exhibit comparable genetic variations at intra-specific and inter-population levels reflected through morphological and cytogenetical diversity in regions, where this plant grows naturally. In the present study, single primer amplification reaction (SPAR) methods have been used for determination of diversity at DNA level in 30 accessions of P. cineraria collected from different districts of Rajasthan. The analyses include the use of six minisatellite core sequence primers for direct amplification of minisatellite DNA (DAMD), eight inter simple sequence repeats (ISSR) and 20 arbitrary primed decamer sequences for random amplification (RAPD) reactions. Upon analysis of the data generated, all the three SPAR methods, either independently and/or in combination, revealed wide range of genetic variation among accessions. Comparison of matrix of individual SPAR method using MxComp component of NTSYS-pc 2.02 K software proving that analysis of natural genetic variation using combination of SPAR methods particularly ISSR and DAMD, rather than an isolated approach, is very effective. Such an approach also yields better information and reflection of the relatedness and affinities at intra-species and inter-population levels. Therefore, it is opined that in order to reveal the intrinsic intra-specific variation, SPAR approaches involving more than one DNA marker may reveal more authentic genetic variation in tropical tree species like P. cineraria.     

Genetic diversity analysis of Perinereis aibuhitensis based on ISSR and SRAP markers of Chinese coast populations

The objective of this study was to obtain an overview of the genetic relationships within Perinereis aibuhitensis using Inter-Simple Sequence Repeat (ISSR) and Sequence-Related Amplified Polymorphism (SRAP) markers that were derived from related populations residing in the Chinese coasts. The percentage of polymorphic bands, Nei's gene diversity and Shannon's information index revealed a high level of genetic diversity at the species level. The analysis of molecular variance revealed that 81.22% (ISSR) and 76.29% (SRAP) of variability were partitioned among individuals within populations, which indicated the coherent trend by Nei's genetic differentiation (Gst) (0.2568/0.2876). The gene flow number (Nm) was 1.4470/1.2385, which indicated that there was limited gene exchange between populations. The phylogenetic tree of the ten P. aibuhitensis populations was separated into four major clusters using the neighbor-joining (NJ) method. These results provide a simple and useful basis for P. aibuhitensis germplasm research and aquaculture breeding.     

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.     

Conservation genetics of endangered medicinal plant Commiphora wightii in Indian Thar Desert

To ascertain the conservation priorities and strategies for Commiphora wightii, an endangered medicinal plant of Indian Thar Desert, genetic diversity was estimated within and among different populations. The total of 155 amplification products were scored using ten each of RAPD and ISSR primers, exhibiting an overall 86.72% polymorphism across 45 individuals representing eight populations. The cumulative data of two markers were used to compute pair-wise distances. The Neighbor-Joining tree revealed high genetic differentiation among populations except Kiradu population. Nei's gene diversity (h) ranged between 0.082 and 0.193 with total diversity at species level is 0.294. Shannon's information index (I) ranged between 0.118 and 0.275 with an overall diversity of 0.439. Analysis of molecular variance showed more diversity among population level (56.65%) than at within population level (43.35%). The low gene flow value (Nm = 0.349) and high coefficient of genetic differentiation (G_ST = 0.589) and high fixation index (F_ST = 0.566) demonstrated elevated genetic differentiation among the population and can be predicted that these populations are not in Hardy–Weinberg proportions. Principal Co-ordinate Analysis confirms that Akal population has become phylogenetically more distinct and less diverse than the rest of the samples. Mantel's test revealed no correlation between genetic and geographical distances of populations (R² = 0.122). Overall highest diversity was observed in the population of Machiya Safari Park and Kiradu, while lowest in Akal population, later may constitute an evolutionary significant unit, having merit for special management.     

Genetic diversity of Indian jujube cultivars using SCoT,ISSR, and rDNA markers

Genetic variation and relationships among 37 cultivars of Ziziphus mauritiana (Lamk.) native of India were analyzed using start codon targeted (SCoT), inter-simple sequence repeats (ISSR), and ribosomal DNA (rDNA) markers. High level of polymorphism among SCoT (61.6%) and ISSR (61%) primers with higher PIC values ranging from 63.1 to 90.4% of SCoT and 47.3 to 88.8% of ISSR primers was recorded. SCoT and ISSR dendrograms revealed similarity coefficients ranging from 0.80 to 0.92 and 0.79 to 0.96, respectively, and clearly delineated all the cultivars of Z. mauritiana into well-supported distinct clusters. Greater Gst signifies higher amount of differentiation observed over multiple loci among seven Z. mauritiana populations. On the other hand, higher gene flow demonstrating a very high migration rate between Z. mauritiana populations indicated higher rates of transfer of alleles or genes from one population to another. The genetic diversity of population 1 (Rajasthan) was the richest among all the seven populations. The largest genetic distance was measured between Maharashtra and West Bengal and the least between Rajasthan and Punjab cultivars. Most of the genetic diversity exists within population rather than among populations. Substantial variation in the ITS-1 region signifies its phylogenetic utility specifically in assessing genetic diversity in Z. mauritiana. The clustering patterns using three molecular marker systems vis-à-vis place of origin exhibited no consistency in grouping of Z. mauritiana cultivars as cultivars from the same place of origin were genetically cataloged into different SCoT, ISSR, and ITS phylogram clusters indicating wide genetic diversity and distribution across agro-climatic zones validating the robustness of marker systems tested.     

Population Genetic Structure of a Widespread Bat-Pollinated Columnar Cactus

Bats are the main pollinators and seed dispersers of Stenocereus thurberi, a xenogamous columnar cactus of northwestern Mexico and a good model to illustrate spatial dynamics of gene flow in long-lived species. Previous studies in this cactus showed differences among populations in the type and abundance of pollinators, and in the timing of flowering and fruiting. In this study we analyzed genetic variability and population differentiation among populations. We used three primers of ISSR to analyze within and among populations genetic variation from eight widely separated populations of S. thurberi in Sonora, Mexico. Sixty-six out of 99 of the ISSR bands (P = 66.7%) were polymorphic. Total heterozygosity for all populations sampled revealed high genetic diversity (H_sp = 0.207, H_BT = 0.224). The AMOVA showed that most of the genetic variation was within populations (80.5%). At the species level, estimates of population differentiation, θ = 0.175 and θ^B = 0.194, indicated moderate gene flow among populations. The absence of a significant correlation between genetic and geographic distances indicated little isolation by geographic distance. The large genetic variation and diversity found in S. thurberi is consistent with its open reproductive system and the high mobility of bats, a major pollinator. However, small changes in number or kind of pollinators and seed dispersal agents, in the directionality of migratory routes, and/or in the timing of flowering and fruiting among populations, can critically affect gene flow dynamics.