Next Article in Journal
Enhancement of Biocontrol Efficacy of Pichia carribbica to Postharvest Diseases of Strawberries by Addition of Trehalose to the Growth Medium
Next Article in Special Issue
Expression of Selected Ginkgo biloba Heat Shock Protein Genes After Cold Treatment Could Be Induced by Other Abiotic Stress
Previous Article in Journal
Antinociceptive and Anti-Inflammatory Effects of Ethanol Extract from Vernonia polyanthes Leaves in Rodents
Previous Article in Special Issue
Recent Advances of Flowering Locus T Gene in Higher Plants
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 13
Issue 3
10.3390/ijms13033900
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Diversity of the Endemic and Medicinally Important Plant Rheum officinale as Revealed by Inter-Simpe Sequence Repeat (ISSR) Markers

by
Xu-Mei Wang
1,*,
Xiao-Qi Hou
2,
Yu-Qu Zhang
2,
Rui Yang
1,
Shi-Fang Feng
2,
Yan Li
1 and
Yi Ren
2
1
College of Medicine, Xi’an Jiaotong University, Xi’an 710061, China
2
College of Life Sciences, Shaanxi Normal University, Xi’an 710062, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2012, 13(3), 3900-3915; https://doi.org/10.3390/ijms13033900
Submission received: 3 February 2012 / Revised: 8 March 2012 / Accepted: 13 March 2012 / Published: 22 March 2012
(This article belongs to the Special Issue Advances in Molecular Plant Biology)
Browse Figures
Versions Notes

Abstract

:
Rheum officinale Baill., an important but endangered medicinal herb, is endemic to China. Inter-simple sequence repeat (ISSR) markers were employed to investigate the genetic diversity and differentiation of 12 populations of R. officinale. Thirteen selected primers yielded 189 bright and discernible bands, with an average of 14.54 per primer. The genetic diversity was low at the population level, but pretty high at the species level (H = 0.1008, I = 0.1505, PPB = 28.95% vs. H = 0.3341, I = 0.5000, PPB = 95.24%, respectively) by POPGENE analysis. Analysis of molecular variance (AMOVA) showed that the genetic variation was found mainly among populations (74.38%), in line with the limited gene flow (Nm = 0.2766) among populations. Mantel test revealed a significant correlation between genetic and geographic distances (r = 0.5381, P = 0.002), indicating the role of geographic isolation in shaping the present population genetic structure. Both Bayesian analysis and UPGMA cluster analysis demonstrated the similar results. Our results imply that the conservation efforts should aim to preserve all the extant populations of this endangered species, and cultivation is proposed in this study.

1. Introduction

To estimate accurately of the genetic diversity of a species is an important element in establishing conservation programs [1–3], because the ability of a species to respond adaptively to environmental changes depends on the level of genetic variability it contains [4]. Therefore, assessing the level and distribution of genetic diversity are crucial for management and the development of effective conservation strategies. Faced with the problem of preserving rare and endangered species, great concerns also should be concentrated on the endemic species with restricted geographic distribution.
Rheum officinale Baill. (Polygonaceae) is a perennial herb, the dried roots and rhizomes of which are called rhubarb (Da Huang in Chinese). Rhubarb is a widely used traditional Chinese medicine having many pharmacological actions, such as purgation, anti-inflammatory, antibacterial, antipyretic, and anticancer effects [5–8]. As described in the Chinese Pharmacopoeia, rhubarb also consists of the roots and rhizomes of R. palmatum L. and R. tanguticum Maxim. ex Balf. [9]. The three genuine species of rhubarb are closely related and monophyletic [10] and endemic to China [11]. As the usage of rhubarb has increased in recent decades, the wild resources have been severely destroyed. Among three genuine species of rhubarb, R. tanguticum has the best pharmacological effect and the most limited distribution in Qinghai, Gansu, Sichuan and Ningxia provinces. Compared with R. tanguticum, R. officinale has a weaker pharmacological effect and wider distribution in western Hubei, western Henan, southern Shaanxi, northern Sichuan, southern Chongqing, northern Guizhou, southeastern of Gansu and northwestern of Yunnan provinces. The distribution of R. palmatum is overlapped with the other two species [11,12]. Because of the overexploitation and the limited distribution, R. tanguticum has become endangered and was listed in the China higher plants endangered list [13,14]. The wild resources of R. officinale face great pressure exacerbated by the reduction in the individual of R. tanguticum. In addition to the increasing medicinal usage of R. officinale, the usage of R. officinale has recently been extended to a functional food and plant protecting agent [15]. The habitat of R. officinale is limited to altitudes ranging from 1100 m to 4600 m, whereby its altitudinal range does not exceed 600 m on any mountain, and it grows generally at the forest edge of the hills and rarely in the forests or in the valleys near the rivers. Because of the deterioration of the habitat of R. officinale, the individuals of this species are decreasing annually, thus the wild resources of it are more easily destroyed. Therefore, it can be predicted that the next original medicinal rhubarb plant which will be endangered will be R. officinale. To date, previous studies of R. officinale have mainly focused on its distribution [12,16], components analysis [17,18] and pharmacological properties [8,19,20]. Although the diversity of R. tanguticum had been studied by a few researchers [13,14], the genetic diversity and population structure of R. officinale remains unknown. In order to conserve and select excellent germplasm for the cultivation in R. officinale, the study on the genetic diversity for this species is becoming necessary and timely.
Of the various molecular markers, ISSR (inter-simple sequence repeat) has a few advantages because ISSR primers anneal directly to simple sequence repeat and thus, unlike SSR markers, no prior knowledge of target sequences is required for ISSR [21,22]. Also, ISSR markers, which have longer primers, allow more stringent annealing temperatures and reveal more polymorphic fragments, can be highly variable within a species and have the advantage over RAPD (random amplified polymorphic DNA) markers [23]. In addition, the cost of the analyses is relatively lower than that of some other markers such as RFLP (restricted fragment length polymorphism) and AFLP (amplified fragment length polymorphism) [24,25]. Therefore, ISSR has been widely used for population genetic studies of various plant species, including many medicinal plants [26–29].
In the present study, ISSR markers were employed to assess the genetic diversity within and among R. officinale populations sampled from all over its distribution. The major objectives of this study were to (1) predict the genetic diversity of the populations of R. officinale; (2) assess the genetic variation within and among populations; and (3) provide basic information for future conservation and management programs in the cultivation of this important medicinal species.

2. Results

2.1. Genetic Diversity

Thirteen primers were chosen to amplify a total of 199 individuals of 12 populations of R. officinale (Table 1). The size of all bands ranged from 200 to 2200 bp and the number of bands by each primer was from 11 to 19 with an average of 14.54 per primer (Table 2). Thirteen primers produced a total of 189 clearly identifiable bands of which 180 were polymorphic, i.e., the percentage of polymorphic bands (PPB) for this species was 95.24% (Table 3). While at the population level, the percentage of polymorphic bands (PPB) ranged from 4.00% to 44.57%, with an average of 29.14%. The average effective number of alleles per locus was 1.1740. Assuming Hardy-Weinberg equilibrium, Nei’s gene diversity (H) varied between 0.0139 and 0.1614, with an average of 0.1014, and Shannon’s information index (I) ranged from 0.0213 to 0.2400, with an average of 0.1514. The values of H and I showed a similar trend to PPB. However, the H and I values equalled 0.3341 and 0.5000 respectively at the species level, demonstrating a relatively high level of genetic diversity (Table 3).

2.2. Genetic Differentiation and Relationships

The result of AMOVA showed that the percentages of genetic variation among populations were 64.38% (Gst) and 74.38% (Φst), both of which indicated that the genetic differentiation was found mainly among populations. A significant (P < 0.001) genetic difference was found among and within populations (Table 4). Furthermore, the level of gene flow (Nm) was measured to be 0.2766 individual per generation between populations, suggesting that gene exchange between populations was very low.
POPGENE analysis revealed genetic distances between populations of R. officinale ranged from 0.1610 (between Pop6 and Pop8) to 0.4732 (between Pop7 and Pop9) (Table 5). The UPGMA tree based on Nei’s unbiased genetic distance [30] was depicted in Figure 1, and indicated that the 12 populations were separated into three geographic groups. The populations (from Pop9 to Pop11) sampled from Yunnan province and Pop 12 from Meigu of Sichuan province formed Group I. Group II consisted of the populations from Pop1 to Pop5, and the remaining populations (from Pop6 to Pop8) formed Group III. Mantel test revealed that a significant correlation between matrices of genetic distance and of geographic distance (r = 0.5381, P = 0.002, 999 permutations). The AMOVA analysis also showed that 32.89% of total genetic variability occurred among the three geographic regions (Table 4). In the ISSR admixture analysis using STRUCTURE (Figure 2), the highest likelihood of the data was obtained when samples were clustered into three groups (K = 3). These groupings were entirely consistent with those of the UPGMA clustering results.
Tajima’s D test [31] was positive and not significant (D = 0.513, P = 0.7187). Fu’s Fs statistic [32] was negative and also not significant (Fs = −5.838, P = 0.0566). Thus, no clear signal of demographic expansion by Fu’s Fs test and Tajima’s D test was revealed for all the populations together. Fu’s Fs tests were significantly negative for the three groups (Fs = −6.855, P = 0.0112 for Group I; Fs = −4.575, P = 0.0376 for Group II and Fs = −1.134, P = 0.0037 for Group III) though Tajima’s D tests were positive for Group I and Group II, and negative for Group III (P > 0.05).

3. Discussion

3.1. Genetic Diversity

The results of the present study showed that the genetic diversity of the R. officinale species was PPB = 95.24% and H = 0.3341, higher than the average values of perennial herbaceous (PPB = 39.30%, H = 0.1240) and of species with endemic distribution (PPB = 40.00%, H = 0.096) [33]. The total genetic diversity of R. officinale was similar to that of R. tanguticum (PPB = 92.94%, H = 0.2689 from Hu et al.) [14]. The present study indicated that R. officinale maintained a higher genetic diversity at the species level. However, in contrast, the genetic diversity at the population level was pretty low (PPB = 4.00–44.57%; H = 0.0139–0.1654).
Although many studies have demonstrated that endangered and endemic species tend to possess low levels of genetic diversity based on ISSR data [34,35], some others have showed opposite findings [36–38]. Genetic diversity of a plant species could be affected by many factors such as distribution range, life form, breeding system and the way that its seeds disperse. A species which has had a long life, a high frequency of gene flow and many seeds tends to have high genetic diversity [33]. Species of Rheum is anthophilous [39]. On the other hand, the trigonous achenes of R. officinale have wings, so they are more easily dispersed by wind than other species of Polygonaceae. In addition, R. officinale is a long-living herbaceous and it can grow more than 20 years and produce large numbers of seeds every year. The germination rate of R. officinale is high [40]. Therefore, the long-living herbaceous habit provides more opportunity to accumulate mutant and abundant seeds and may contribute to the high amount of diversity in this species [14,41,42].
Of all the studied populations, the Pop2 population had the highest genetic diversity (PPB = 44.57%, H = 0.1614, Table 3). The Pop2 is located in Shennongjia National Nature Reserve, Hubei province. This area possesses a unique and well-preserved subtropical forest ecosystem and large numbers of rare and endangered species are on the list of UNESCO’s MAB World’s Reserve Network [43,44]. Hence, it is an ideal habitat for R. officinale. In contrast, genetic diversity of the Pop6 population from Mt. Jinfo in Chongqing was obviously lower than other populations. The reason might be that there are only five individuals in this population. In contrast with Shennongjia National Nature Reserve, Mt. Jinfo is a famous tourism area, and human activities frequently happen in this area. Moreover, this area is now a forest fragment surrounded by agricultural lands. Fragmentation may increase levels of inbreeding and/or genetic drift [13]. The Pop6 population is located in a smaller reserved forest compared with the other populations. According to our field investigation, the wild resources have been overexploited greatly in this site and we only found a few individuals in one cluster. From the results of genetic diversity we speculated that five individuals might be propagated from one mother plant. The results also suggested that the sampling strategy is an important factor in genetic diversity.

3.2. Genetic Differentiation

High genetic differentiation among populations of R. officinale was detected based on Nei’s gene diversity (Gst = 0.6438) and AMOVA (Φst = 0.7438), which is similar to the higher genetic differentiation in endangered and endemic species in Megacodon stylophorus (Gentianaceae) [45], Rhodiola alsia (Crassulaceae) [46] and Torreya jackii [38]. Correspondingly, the among-population differentiation (Φst) is 0.807, 0.703 and 0.670, respectively. However, the present study showed an opposite finding to the previous studies of R. tanguticum [13,14]. Hu et al. [14] reported that the among-population differentiation coefficients were 0.3585 (Gst) and 0.290 (Φst), and that from Chen et al. were lower (Gst = 0.249, Φst = 0.2118) [13]. The reason for the contradiction between these different studies might be that the populations of R. tanguticum were collected only from Qinghai province in their studies. In fact, R. tanguticum was also distributed in Gansu, Sichuan and Ningxia provinces [12]. The sampled populations of R. officinale in the present study covered the entire distribution. As Pfeifer and Jetschke [47] reported geographic isolation is one major factor influencing genetic differentiation by limiting the amount of gene flow via both pollen and seeds. In the present study, the largest and mean inter-population geographic distances were 1174.91 km (Pop1 vs. Pop11) and 591.854 km, respectively. R. officinale is outcrossing and its seeds are dispersed by wind. Therefore, large geographic barriers, such as Mt. Qinling and Mt. Daba and the human activities between populations greatly hindered gene flow via seed and pollen dispersal among populations. Typically, when Nm < 1, differential selection might be very strong, and population differentiation might be maintained [48]. In R. officinale, Nm estimated from Gst was only 0.2766 and that is far below 1, indicating the gene flow among R. officinale populations with discontinuous distributions was limited and might enhance the genetic differentiation among populations. Moreover, the Mantel test revealed a significant correlation between genetic and geographic distances (r = 0.5381, P = 0.002), displaying a clear isolation of populations by geographic distances in shaping the present genetic structure of R. officinale. This genetic structure was further confirmed by the topology of UPGMA (Figure 1) and the Bayesian analyses (Figure 2). On the other hand, the roots and rhizomes of R. officinale have been overexploited in the past decades, and its habitat gradually limited to small isolated areas. For example, the size of some populations like Pop6 and Pop8 in this study is small, and only a few individuals can be found in the field. The small populations are likely to have been subjected to strong genetic drift. Genetic drift changes the distribution of genetic variation in two ways, by reducing variation within populations and by increasing differentiation between populations [49]. Meanwhile, the theory of population genetics indicates that genetic differentiation should increase in progressively smaller and more isolated populations [49,50].
In the UPGMA tree, populations were mainly divided into three geographic groups (Figure 1). Group I included the populations Pop9, Pop10, Pop11 and Pop12, while the remaining populations were clustered into two other groups. This division reflects geographic distribution pattern of these populations. The populations Pop9, Pop10, Pop11 and Pop12 are geographically distant from the other populations and are located in Hengduan mountain system (southwestern China). Group II included the populations from Pop1 to Pop5 which were sampled from central and eastern part of Mt. Qinling and eastern Mt. Daba, while Pop6, Pop7 and Pop8 sampled from Mt. Dalou (Pop 7 and Pop 8) and Pop6 from western Mt. Daba formed Group III. Populations from Mt. Hengduan which is located in the southeastern part of the Qinghai-Tibet Plateau is in the southwestern part of the distribution of R. officinale, whereas populations from central and eastern Mt. Qinling and eastern Mt. Daba are in the eastern part of the distribution of R. officinale. This finding was consistent with the results of the Bayesian analyses, which was clearly identified the same three main group clusters. Hence, the present genetic structure could be partly explained by the isolation-by-distance model [51], as indicated by the result of the Mantel test. However, the genetic differentiation among the three groups was lower than within groups. The reason of this may be that the weak expansion existed within each of the three groups respectively, as indicated by the neutrality deviation in Fu’s Fs tests.
Combined with the geographic distribution pattern of R. officinale, the genetic diversity of populations which were collected from southwestern part of its distribution (i.e., Pop9-Pop12 in Mt. Hengduan) was lower than that of the populations from the central (e.g., in Mt. Dalou and western Mt. Daba) and northeastern of its distribution (central and eastern Mt. Qinling and eastern Mt. Daba). It is well known that China’s topography is divided into three terraces from the west to the east: Mt. Hengduan lies in the first terrace which mainly includes the Qinghai-Tibet Plateau, with the altitude higher than 3000 m; the second terrace refers to the eastern and northern areas out of the Qinghai-Tibet Plateau, with an average altitude of lower than 3000 m; the third terrace is mainly hills and plains in the east, with an average altitude of lower than 500 m. R. officinale is distributed in the first and the second terrace. Populations Pop9 to Pop12 are located in Mt. Hengduan which belongs to the first terrace, and the altitudes of those populations are higher than 3400 m. The other populations are located in the central and the northeastern area of the species’ distribution in the second terrace, with an altitude of lower than 2900 m. Compared with the second terrace, the first terrace is the youngest orogeny of the mountains in Neozoic era, and the genetic diversity of the populations located in this area is low. The reason may be that accompanying with the uplifts of Mt. Hengduan, the distribution of the species extends to the higher altitude, which might cause the disjunctive distribution of R. officinale. This disjunctive distribution pattern of R. officinale may be another reason for the present genetic structure. These areas with lower genetic diversity are not the suitable habitat for the growth of the species.

3.3. Implications for Conservation and Cultivation

The primary objective in nature conservation is to preserve as much as possible of the evolutionary potential of species through maintaining as much genetic diversity as possible [34]. The maintenance of genetic variation is a major objective within conservation plans for endangered species [52,53]. Information obtained in the present study provides significant implications for conservation strategies of R. officinale. The results reported here revealed low genetic diversity at population level and high genetic diversity at species level in R. officinale. High genetic differentiation that occurred among populations and might be due to limited gene flow and genetic drift. In general, this species is not genetically depauperate. The main factor responsible for this threatened species may be recent over-collection of the medicinal organs from the field, rather than a lack of overall genetic diversity. Therefore, the exploitation of wild resources should be forbidden. Considering the high genetic differentiation of R. officinale, preservation of any population would be insufficient to conserve all the variation in the species. Thus, the priority must be to protect all the existing populations in situ, especially those extant populations with high levels of genetic variation of different regions such as Pop1, Pop2 and Pop5, we suggest that their habitats be protected and the exploitation of wild resources be prevented. If ex situ conservation is required, samples should be collected from as many populations as possible from the whole distribution, and populations with small sizes (e.g., Pop6 and Pop8) should receive more attention. On the other hand, to meet the bulk commercial demand for this traditional medicinal plant, cultivation facilities can be established as an alternative source of raw materials. Therefore, the high genetic diversity of the materials (e.g., Pop1, Pop2 and Pop5) should be used in cultivation in order to reduce the inbreeding depression.

4. Experimental Section

4.1. Plant Sampling

Twelve populations were sampled throughout the distribution of R. officinale, including Henan, Hubei, Shaanxi, Sichuan, Chongqing, Guizhou and Yunnan provinces, China. 16–20 (only five and ten individuals in Chongqing and Guizhou populations, respectively) fresh leaves were collected randomly in each population, depending on accessibility and population size. Each population was positioned by a GPS and the detailed locations of the studied populations are listed in Table 1. The young leaves were stored and dried in ziplock bags with silica gel and transported back to our laboratory for DNA extraction. The voucher specimens were deposited in the Shaanxi Normal University Herbarium (SANU).

4.2. DNA Extraction and PCR Amplification

Total DNA was extracted from the silica gel-dried leaves using the modified 2× CTAB procedure [54]. The quality and quantity of DNA were performed by UV-spectrophotometer (ND-2000, NanoDrop, USA). DNA concentration and purity were also determined by electrophoresis on 1.0% agarose gels based on the intensities of band when compared with l kb plus DNA ladder as marker. The DNA samples were diluted to the concentration of 50 ng/μL and stored at −20 °C for use.
One hundred ISSR primers synthesized by Shanghai Sangon Biological Engineering Technology & Service (China), according to the primer set published by University of British Columbia, Canada (UBC set No.9) were used for amplification to standardize the PCR conditions. Thirteen of 100 ISSR primers produced clear, reproducible and relatively high polymorphism bands were selected for all samples of R. officinale (Table 2). The effects of Mg2+, dNTPs, DNA templates, primers concentrations and annealing temperature on the amplification were tested, and the final amplification was carried out according to our previous study for the optimized ISSR-PCR reactions of R. officinale [55]. PCR products were electrophoresed on 1.6% (w/v) agarose gels, in 1× TBE Buffer at 110 V for 1.5 h and stained with ethidium bromide (0.5 μg/mL). Gels with amplification fragments were visualized and photographed in UV light by using Bio-Rad Gel Documentation System (Bio-Rad Laboratories, UK). DL2000 ladder (TaKaRa Biotechnology, China) was used as DNA molecular weight.

4.3. Data Analysis

The amplified fragments, with the same mobility according to their molecular weight (bp), were scored in terms of a binary code as present (1) or absent (0). Only those consistently reproducible bands were scored, and smeared and weak bands were excluded. Genetic parameters including the percentage of polymorphic bands (PPB), observed number of alleles (Na), the effective number of alleles (Ne), Nei’s gene diversity (H) [56], Shannon’s index (I) [57], Nei’s genetic identity and genetic distance [30], Nei’s genetic differentiation index among populations (GST) [56] and gene flow (Nm) were calculated using the computer program POPGENE version 1.32 [58]. An estimate of Nm among populations was computed using the formula of Nm =0.5(1 − GST)/GST [59]. The obtained genetic distance matrix was then used to perform the cluster analysis and construct the unweighted pair-group method with arithmetic average (UPGMA) dendrogram using MEGA 4.0 [60]. In addition, an analysis of molecular variance (AMOVA) procedure was used to estimate the partitioning of genetic variance among and within populations. Input data files for the AMOVA 1.55 [61] were generated using DCFA version 1.1 [62]. The significance of variance components were tested statistically by nonparametric randomization tests using 1000 permutations. Geographic distances were interpreted by the latitudes and longitudes with Mapinfo 8.0 Program. The Mantel test of genetic and geographic distances was carried out to evaluate the correlation between the two data matrices using TFPGA software [63] (computing 999 permutations). Tajima’s D test [31] and Fu’s Fs test [32] were conducted to test the neutral mutation hypothesis using ARLEQUIN Version 3.0 [64]. The demographic history of a population could be inferred by comparing such neutrality tests, given that a range expansion is suggested when Tajima’s D and Fu’s Fs are significantly negative [32].
A Bayesian analysis of ISSR population structure was performed on the entire data set using the program STRUCTURE (version 2.3) [65] to detect population structure and estimate the number of populations (K) in a sample and to assign individuals to one or more of these populations (K). The number of genetically distinct clusters (K) was set to vary from 1 to 12 (the total number of populations). The model was run for eight independent simulations for each K, used a burn-in length of 50,000 and a run length of 100,000 iterations. Following the program’s dominant marker settings, the “no admixture” model was used, and uncorrelated allele frequencies among populations were assumed. The most likely number of clusters was estimated according to the model value (ΔK) based on the second order rate of change, with respect to K, of the likelihood function, following the procedure described by Evanno et al. [66].

5. Conclusions

In summary, our results indicated that the genetic diversity of R. officinale was high at the species level, while the genetic diversity was low at the population level. High genetic differentiation was found mainly among populations, which may be attributed to geographic isolation by the limited gene flow and genetic drift. Cluster analysis using the UPGMA method grouped all populations into three geographic groups. Based on these findings, strategies are proposed for the conservation and cultivation of the species.

Acknowledgments

This work was supported by National Nature Science Foundation of China (No. 81001602) and Fundamental Research Funds for the Central Universities (No. xjj20100038).

References

  1. Cote, C.T. Genetic variation in rare and common plants. Annu. Rev. Ecol. Syst 2003, 34, 213–237. [Google Scholar]
  2. Qian, W.; Ge, S.; Hong, D.Y. Genetic variation within and among populations of a wild rice Oryza granulata from China detected by RAPD and ISSR markers. Theor. Appl. Genet 2001, 102, 440–449. [Google Scholar]
  3. Kim, S.C.; Lee, C.; Santos-Guerra, A. Genetic analysis and conservation of the endangered Canary Island woody sow-thistle, Sonchus gandogeri (Asteraceae). J. Plant Res 2005, 118, 147–153. [Google Scholar]
  4. Ayala, F.J.; Kiger, J.A. Modern Genetics, 2nd ed; Benjamin/Cumings: Menlo Park, CA, USA, 1984. [Google Scholar]
  5. Xiao, P.; He, L.; Wang, L. Ethnopharmacologic study of Chinese rhubarb. J. Ethnopharmacol 1984, 10, 275–293. [Google Scholar]
  6. Tseng, S.H.; Lee, H.H.; Chen, L.G.; Wu, C.H.; Wang, C.C. Effects of three purgative decoctions on inflammatory mediators. J. Ethnopharmacol 2006, 105, 118–124. [Google Scholar]
  7. Liu, L.; Guo, Z.J.; Lv, Z.G.; Sun, Y.; Cao, W.; Zhang, R.; Liu, Z.G.; Li, C.; Cao, S.S.; Mei, Q.B. The beneficial effect of Rheum tanguticum polysaccharide on protecting against diarrhea, colonic inflammation and ulceration in rats with TNBS-induced colitis: The role of macrophage mannose receptor in inflammation and immune response. Int. J. Immunopharmacol 2008, 8, 1481–1492. [Google Scholar]
  8. Li, W.Y.; Chan, S.W.; Guo, D.J.; Chung, M.K.; Leung, T.Y.; Yu, P.H. Water extract of Rheum officinale Baill. induces apoptosis in human lung adenocarcinoma A549 and human breast cancer MCF-7 cell lines. J. Ethnopharmacol 2009, 124, 251–256. [Google Scholar]
  9. Chinese Pharmacopoeia Committee, Pharmacopoeia of the People’s Republic of China; China Medical Science and Technology Press: Beijing, China, 2010; pp. 22–23.
  10. Wang, A.L.; Yang, M.H.; Liu, J.Q. Molecular phylogeny, recent radiation and evolution of gross morphology of the rhubarb genus Rheum (Polygonaceae) inferred from chloroplast DNA trnL-F sequences. Ann. Bot 2005, 96, 489–498. [Google Scholar]
  11. Bao, B.J.; Grabovskaya-Borodina, A.E. Rheum. In Flora of China; Li, A.R., Bao, B.J., Eds.; Science Press: Beijing, China; Missouri Botanical Garden: St. Louis, USA, 2003; Volume 5, pp. 341–350. [Google Scholar]
  12. Wang, X.M.; Hou, X.Q.; Zhang, Y.Q.; Li, Y. Distribution pattern of genuine species of rhubarb as traditional Chinese medicine. J. Med. Plants Res 2010, 4, 1865–1876. [Google Scholar]
  13. Chen, F.J.; Wang, A.L.; Chen, K.M.; Wan, D.S.; Liu, J.Q. Genetic diversity and population structure of the endangered and medically important Rheum tanguticum (Polygonaceae) revealed by SSR markers. Biochem. Syst. Ecol 2009, 37, 613–621. [Google Scholar]
  14. Hu, Y.P.; Wang, L.; Xie, X.L.; Yang, J.; Li, Y.; Zhang, H.G. Genetic diversity of wild populations of Rheum tanguticum endemic to China as revealed by ISSR analysis. Bioch. Syst. Ecol 2010, 38, 264–274. [Google Scholar]
  15. Yang, X.J.; Ma, X.X.; Yang, L.J.; Yu, D.Z.; Qian, Y.X.; Ni, H.W. Efficacy of Rheum officinale liquid formulation on cucumber powdery mildew. Crop Prot 2009, 28, 1031–1035. [Google Scholar]
  16. Shen, S.L.; Li, Z.G.; Zhang, J.Y. Resources distribution and usage of medicinal plants of Rheum in Gansu. J. Gansu Sci 1996, 8, 20–26. [Google Scholar]
  17. Yang, F.Q.; Zhang, T.Y.; Tian, G.L.; Cao, H.F.; Liu, Q.H.; Ito, Y. Preparative isolation and purification of hydroxyanthraquinones from Rheum officinale Baill. by high-speed counter-current chromatography using pH-modulated stepwise elution. J. Chromatogr. A 1999, 858, 103–107. [Google Scholar]
  18. Xu, Z.X.; Sun, Y.; Zhang, H.Q. Study on the anti-aging effect of rhapontinum. J. Anhui Med. Pharmaceu 2002, 6, 9–10. [Google Scholar]
  19. Tsai, J.C.; Tsai, S.; Chang, W.C. Effect of ethanol extracts of three Chinese medicinal plants with laxative properties on ion transport of the rat intestinal epithelia. Biol. Pharm. Bull 2004, 27, 162–165. [Google Scholar]
  20. Moon, M.K.; Kang, D.G.; Lee, J.K.; Kim, J.S.; Lee, H.S. Vasodilatory and anti-inflammatory effects of the aqueous extract of rhubarb via a NO-cGMP pathway. Life Sci 2006, 78, 1550–1557. [Google Scholar]
  21. Zietkiewicz, E.; Rafalski, A.; Labuda, D. Genome fingerprinting by simple sequence repeat (SSR)-anchored polymerase chain reaction amplification. Genomics 1994, 20, 176–183. [Google Scholar]
  22. Godwin, I.D.; Aitken, E.A.; Smith, L.W. Application of inter simple sequence repeat (ISSR) markers to plant genetics. Electrophoresis 1997, 18, 1524–1528. [Google Scholar]
  23. Esselman, E.J.; Li, J.Q.; Crawford, D.J.; Winduss, J.L.; Wolfe, A.D. Clonal diversity in the rare Calamagrostis porteri ssp. insperata (Poaceae): Comparative results for allozymes and randomamplified polymorphic DNA (RAPD) and inter simple sequence repeat (ISSR) markers. Mol. Ecol 1999, 8, 443–451. [Google Scholar]
  24. Reddy, M.P.; Sarla, N.; Siddiq, E.A. Inter simple sequence repeat (ISSR) polymorphism and its application in plant breeding. Euphytica 2002, 128, 9–17. [Google Scholar]
  25. Song, Z.Q.; Li, X.F.; Wang, H.G.; Wang, J.H. Genetic diversity and population structure of Salvia miltiorrhiza Bge. in China revealed by ISSR and SRAP. Genetica 2010, 138, 241–249. [Google Scholar]
  26. Camacho, F.J.; Liston, A. Population structure and genetic diversity of Botrychium pumicola (Ophioglossaceae) based on inter-simple sequence repeats (ISSR). Am. J. Bot 2001, 88, 1065–1070. [Google Scholar]
  27. Lei, Y.D.; Gao, H.; Tsering, T.; Shi, S.H.; Zhong, Y. Determination of genetic variation in Rhodiola crenulata from the Hengduan Mountains Region, China using inter-simple sequence repeats. Genet. Mol. Biol 2006, 29, 339–344. [Google Scholar]
  28. Xiao, M.; Li, Q.; Wang, L.; Guo, L.; Li, J.; Tang, L.; Chen, F. ISSR Analysis of the genetic diversity of the endangered species Sinopodophyllum hexandrum (Royle) Ying from western Sichuan Province, China. J. Integr. Plant Biol 2006, 48, 1140–1146. [Google Scholar]
  29. Qiu, Y.Q.; Zong, M.; Yao, H.; Chen, B.L.; Zhou, X.L.; Chen, Z.L.; Fu, C.X. Genetic variation in wild and cultivated Rhizoma Corydalis revealed by ISSRs markers. Planta Med 2008, 74, 1–5. [Google Scholar]
  30. Nei, M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 1978, 89, 583–590. [Google Scholar]
  31. Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 1989, 123, 585–595. [Google Scholar]
  32. Fu, Y.X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 1997, 147, 915–925. [Google Scholar]
  33. Hamrick, J.L.; Godt, M.J.W. Allozyme Diversity in Plant Species. In Plant Population Genetics, Breeding, and Genetic Resources; Brown, A.H.D., Clegg, M.T., Kahler, A.L., Weir, B.S., Eds.; Sinauer Associates: Sunderland, MA, USA, 1990; pp. 43–63. [Google Scholar]
  34. Xiao, L.Q.; Ge, X.J.; Gong, X.; Hao, G.; Zheng, S.X. ISSR variation in the endemic and endangered plant Cycas guizhouensis (Cycadaceae). Ann. Bot 2004, 94, 133–138. [Google Scholar]
  35. Li, F.G.; Xia, N.H. Population structure and genetic diversity of an endangered species, Glyptostrobus pensilis (Cupressaceae). Bot. Bull. Acad. Sin 2005, 46, 155–162. [Google Scholar]
  36. Ge, Y.Q.; Qiu, Y.X.; Ding, B.Y.; Fu, C.X. An ISSR analysis on population genetic diversity of the relict plant Ginkgo biloba. Biodiv. Sci 2003, 11, 276–287. [Google Scholar]
  37. Luan, S.; Chiang, T.Y.; Gong, X. High genetic diversity vs. low genetic differentiation in Nouelia insignis (Asteraceae), a narrowly distributed and endemic species in China, revealed by ISSR fingerprinting. Ann. Bot 2006, 98, 583–589. [Google Scholar]
  38. Li, J.M.; Jin, Z.X. Genetic variation and differentiation in Torreya jackii Chun, an endangered plant endemic to China. Plant Sci 2007, 172, 1048–1053. [Google Scholar]
  39. Knuth, P. Handbook of Flower Pollination Based upon Hermann Müller’s Work ‘The Fertilisation of Flowers by Insects’, 2nd ed; Clarendon Press: Oxford, UK, 1906; Volume 3, p. 340. [Google Scholar]
  40. Xiao, S.P.; Chen, M.; Huang, L.Q.; Gao, F. Primary study on shapes of fruits and germination characters of seeds of Radix et Rhizoma Rhei. China J. Chin. Mat. Med 2007, 32, 195–199. [Google Scholar]
  41. Nybom, H. Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants. Mol. Ecol 2004, 13, 1143–1155. [Google Scholar]
  42. Zhao, Y.; Chen, X.Y.; Wang, X.R.; Pian, R.Q. ISSR analysis of genetic diversity among Lespedeza bicolor populations. J. Plant Genet. Resour 2007, 8, 195–199. [Google Scholar]
  43. Chen, L.Z. Present Situation of Biodiversity and Its Protection in China; Science Press: Beijing, China, 1993; pp. 12–43. [Google Scholar]
  44. Fu, C.J.; Liu, Y.H.; Xiong, W.J. Regeneration characteristics of gaps and their effect on Abies fargesii forests in the Shennongjia National Nature Reserve, China. For. Stud. China 2010, 12, 206–212. [Google Scholar]
  45. Ge, X.J.; Zhang, L.B.; Yuan, Y.M.; Hao, G.; Chang, T.Y. Strong genetic differentiation of the East-Himalayan Megacodon stylophorus (Gentianaceae) detected by inter-simple sequence repeats (ISSR). Biodiv. Conserv 2005, 14, 849–861. [Google Scholar]
  46. Xia, T.; Chen, S.L.; Chen, S.Y.; Ge, X.J. Genetic variation within and among populations of Rhodiola alsia (Crassulaceae) native to the Tibetan Plateau as detected by ISSR markers. Biochem. Genet 2005, 43, 87–101. [Google Scholar]
  47. Pfeifer, M.; Jetschke, G. Influence of geographical isolation on genetic diversity of Himantoglossum hircinum (Orchidaceae). Folia Geobot 2006, 41, 3–20. [Google Scholar]
  48. Slatkin, M. Gene flow in natural populations. Annu. Rev. Ecol. Syst 1985, 16, 393–430. [Google Scholar]
  49. Ellstrand, N.C.; Elam, D.R. Population genetic consequences of small population size: Implications for plant conservation. Annu. Rev. Ecol. Syst 1993, 24, 217–242. [Google Scholar]
  50. Templeton, A.R.; Shaw, K.; Routman, E.; Davis, S.K. The genetic consequences of habitat fragmentation. Ann. Mo. Bot. Gard 1990, 77, 13–27. [Google Scholar]
  51. Wright, S. Isolation by distance. Genetics 1943, 28, 114–138. [Google Scholar]
  52. Milligan, B.G.; Leebens-Mack, J.; Strand, A.E. Conservation genetics: Beyond the maintenance of marker diversity. Mol. Ecol 1994, 12, 844–855. [Google Scholar]
  53. Avise, J.C.; Hamrick, J.L. Conservation Genetics: Case Histories from Nature; Chapman and Hall: New York, NY, USA, 1996. [Google Scholar]
  54. Doyle, J.J.; Doyle, J.L. A rapid DNA isolation procedure for small quantities of fresh leaf material. Phytochem. Bull 1987, 19, 11–15. [Google Scholar]
  55. Wang, X.M. Optimization of DNA isolation, ISSR-PCR system and primers screening of genuine species of rhubarb, an important herbal medicine in China. J. Med. Plants Res 2010, 4, 904–908. [Google Scholar]
  56. Nei, M. Analysis of gene diversity in subdivided populations. Proc. Natl. Acad. Sci. USA 1973, 70, 3321–3323. [Google Scholar]
  57. Lewontin, R.C. Apportionment of human diversity. Evol. Biol 1972, 6, 381–398. [Google Scholar]
  58. Yeh, F.; Yang, R.; Boyle, T. POPGENE Version 1.32; Microsoft Window-based Freeware for Population Genetic Analysis; Molecular Biology and Biotechnology Center, University of Alberta: Edmonton, AB, Canada, 1999. [Google Scholar]
  59. McDermott, J.; McDonald, B. Gene flow in plant pathosystems. Annu. Rev. Phytopathol 1993, 31, 353–373. [Google Scholar]
  60. Tamura, K.; Dudley, J.; Nei, M.; Kumar, S. MEGA 4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol. Biol. Evol 2007, 24, 1596–1599. [Google Scholar]
  61. Excoffier, L.; Smouse, P.E.; Quattro, J.M. Analysis of molecular variance inferred from metric distances among DNA haplotypes: Applications to human mitochondrial DNA restriction data. Genetics 1992, 131, 479–491. [Google Scholar]
  62. Zhang, F.M.; Ge, S. Data analysis in population genetics. I. Analysis of RAPD data with AMOVA. Biodivers. Sci 2002, 10, 438–444. [Google Scholar]
  63. Miller, M.P. Tools for Population Genetic Analyses (TFPGA) Version 1.3: A Windows program for the Analysis of Allozyme and Molecular Population Genetic Data; Utah State University: Logan, UT, USA, 1997. [Google Scholar]
  64. Excoffier, L.; Laval, G.; Schneider, S. Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol. Bioinform 2005, 1, 47–50. [Google Scholar]
  65. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar]
  66. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol 2005, 14, 2611–2620. [Google Scholar]
Ijms 13 03900f1 1024
Figure 1. UPGMA dendrogram based on Nei’s (1978) genetic distances among populations.
Figure 1. UPGMA dendrogram based on Nei’s (1978) genetic distances among populations.
Ijms 13 03900f1
Ijms 13 03900f2 1024
Figure 2. Genetic relationships among the 12 populations of R. officinale estimated using STRUCTURE program based on ISSR data. The model with K = 3 showed the highest ΔK value.
Figure 2. Genetic relationships among the 12 populations of R. officinale estimated using STRUCTURE program based on ISSR data. The model with K = 3 showed the highest ΔK value.
Ijms 13 03900f2
Table
Table 1. Sampling details of R. officinale populations in the present study.
Table 1. Sampling details of R. officinale populations in the present study.
PopulationLocalityLongitude (E)Latitude (N)Altitude (m)Sample sizeVoucher
Pop1Baotianman Reserve, Mt. Funiu, Neixiang county, Henan33°29.179′111°54.96′110020Xu-mei Wang 08082417
Pop2Muyu, Mt. Shennongjia, Xingshan county, Hubei31°27.088′110°16.172′290818Xu-mei Wang and Xiao-qi Hou 09071101
Pop3Doumugong, Mt. Taibai, Mei county, Shannxi34°02.286′107°42.869′284119Xiao-qi Hou 10072001
Pop4Taibaimiao, Ningshan county, Shaanxi33°25.406′108°31.833′187820Xiao-qi Hou 10071301
Pop5Mt. Hualong, Pingli county, Shaanxi32°01.393′109°21.499′291920Xiao-qi Hou 2010072801
Pop6Daping, Mt. Jinfo, Nanchuan county, Chongqing28°58.414′107°11.023′14125Xu-mei Wang and Xiao-qi Hou 09072711
Pop7Nanjiang county, Sichuan32°35.668′107°06.78′180919Yu-qu Zhang 10080801
Pop8Hailongtun, Gaoping, Zunyi county, Guizhou27°48.766′106°49.097′125210Xu-mei Wang and Xiao-qi Hou 20090731
Pop9Gelachang, Haba, Sanba, Xianggelila county, Yunan27°23.219′100°02.754′399516Xiao-qi Hou 09080718
Pop10Parch house, Haba, Sanba, Xianggelila county, Yunan27°23.741′100°02.257′372716Xiao-qi Hou 09080720
Pop11Xiaozhongdian, Xianggelila county, Yunan27°35.118′99°50.835′344120Xiao-qi Hou 1010202
Pop12Hongxi, Meigu county, Sichuan28°40.243′102°58.341′362316Yu-qu Zhang and Xiao-qi Hou 09091824
Table
Table 2. Inter-simple sequence repeat (ISSR) primers used for ISSR analysis in the present study, Y = (C, T); B = (C, G, T); D = (A, G, T); H = (A, G, T); V = (A, C, G).
Table 2. Inter-simple sequence repeat (ISSR) primers used for ISSR analysis in the present study, Y = (C, T); B = (C, G, T); D = (A, G, T); H = (A, G, T); V = (A, C, G).
Primer codeSequence (5′→3′)Annealing temperature (°C)No. of amplified bandsNo. of polymorphic bands
UBC807(AG)8T511212
UBC811(GA)8C531311
UBC816(CA)8T521615
UBC825(AC)8T521211
UBC834(AG)8YT521414
UBC835(AG)8YC52109
UBC836(AG)8YA521616
UBC841(GA)8YC521514
UBC842(GA)8YG561111
UBC888BDB(CA)7521918
UBC889DBD(AC)7521717
UBC890VHV(GT)7561817
UBC891HVH(TG)7521615
Total--189180
Table
Table 3. Genetic diversity within populations of R. officinale, Na: observed number of alleles; Ne: effective number of alleles; H: Nei’s (1973) gene diversity; I: Shannon’s information index; PPB: percentage of polymorphic bands.
Table 3. Genetic diversity within populations of R. officinale, Na: observed number of alleles; Ne: effective number of alleles; H: Nei’s (1973) gene diversity; I: Shannon’s information index; PPB: percentage of polymorphic bands.
PopulationsNaNeHIPPB (%)
Pop11.44571.26260.15550.233244.57
Pop21.44571.27710.16140.240044.57
Pop31.39431.22100.13050.197039.43
Pop41.36571.21560.12510.187736.57
Pop51.43431.27300.15630.231843.43
Pop61.04001.02130.01390.02134.00
Pop71.32001.11290.06990.111532.00
Pop81.11431.07660.04350.064011.43
Pop91.24571.17370.09900.144724.57
Pop101.22291.16140.09140.132922.29
Pop111.17711.11900.06880.101417.71
Pop121.26861.16130.09400.140226.86
Average1.28951.17300.10080.150528.95
Species level1.95431.56960.33410.500095.24
Table
Table 4. Analysis of molecular variance (AMOVA) from 12 populations of R. officinale using 13 inter-simple sequence repeat markers, d.f., degree of freedom; SSD, sum of squares; MSD, mean squared deviations; VC, variance component; TVP, total variance percentage;
Table 4. Analysis of molecular variance (AMOVA) from 12 populations of R. officinale using 13 inter-simple sequence repeat markers, d.f., degree of freedom; SSD, sum of squares; MSD, mean squared deviations; VC, variance component; TVP, total variance percentage;
Source of variationd.f.SSDMSDVCTVP (%)P-value *
Among populations112076.62207.6624.3774.38<0.001
Within populations187663.008.398.3925.62<0.001
Among geographic regions2712.42356.2111.4232.89%<0.001
Within geographic regions1962027.2023.3023.3067.11%<0.001
*Significance tests after 1000 permutations.
Table
Table 5. Nei’s (1978) unbiased measures of genetic distance (below diagonal) and genetic identity (above diagonal) between R. officinale populations.
Table 5. Nei’s (1978) unbiased measures of genetic distance (below diagonal) and genetic identity (above diagonal) between R. officinale populations.
PopulationPop1Pop2Pop3Pop4Pop5Pop6Pop7Pop8Pop9Pop10Pop11Pop12
Pop1****0.76740.80840.74910.79280.69600.71650.72380.69440.72790.70720.7026
Pop20.2648****0.74550.77830.75130.66900.69040.69200.68900.70540.67780.6903
Pop30.21270.2937****0.80220.81960.69540.72580.75200.70270.69220.67840.7258
Pop40.28890.25070.2205****0.81080.74200.74080.77700.67700.71400.72890.7044
Pop50.23220.28590.19900.2098****0.73470.73130.76470.70840.71730.73500.6925
Pop60.36250.40190.36330.29840.3083****0.74990.85130.63490.66270.63830.6522
Pop70.33340.37040.32040.30000.31290.2878****0.79970.62300.69610.65430.6810
Pop80.32320.36820.28500.25230.26830.16100.2235****0.68510.71190.66420.6888
Pop90.36470.37250.35280.39000.34470.45420.47320.3783****0.80640.71080.7865
Pop100.31750.34900.36780.33690.33230.41140.36230.33980.2152****0.73640.7803
Pop110.34640.38900.38800.31620.30790.44900.42410.40920.34130.3059****0.7040
Pop120.35290.37070.32040.35040.36740.42750.38420.37280.24010.24810.3510****

Share and Cite

MDPI and ACS Style

Wang, X.-M.; Hou, X.-Q.; Zhang, Y.-Q.; Yang, R.; Feng, S.-F.; Li, Y.; Ren, Y. Genetic Diversity of the Endemic and Medicinally Important Plant Rheum officinale as Revealed by Inter-Simpe Sequence Repeat (ISSR) Markers. Int. J. Mol. Sci. 2012, 13, 3900-3915. https://doi.org/10.3390/ijms13033900

AMA Style

Wang X-M, Hou X-Q, Zhang Y-Q, Yang R, Feng S-F, Li Y, Ren Y. Genetic Diversity of the Endemic and Medicinally Important Plant Rheum officinale as Revealed by Inter-Simpe Sequence Repeat (ISSR) Markers. International Journal of Molecular Sciences. 2012; 13(3):3900-3915. https://doi.org/10.3390/ijms13033900

Chicago/Turabian Style

Wang, Xu-Mei, Xiao-Qi Hou, Yu-Qu Zhang, Rui Yang, Shi-Fang Feng, Yan Li, and Yi Ren. 2012. "Genetic Diversity of the Endemic and Medicinally Important Plant Rheum officinale as Revealed by Inter-Simpe Sequence Repeat (ISSR) Markers" International Journal of Molecular Sciences 13, no. 3: 3900-3915. https://doi.org/10.3390/ijms13033900

APA Style

Wang, X. -M., Hou, X. -Q., Zhang, Y. -Q., Yang, R., Feng, S. -F., Li, Y., & Ren, Y. (2012). Genetic Diversity of the Endemic and Medicinally Important Plant Rheum officinale as Revealed by Inter-Simpe Sequence Repeat (ISSR) Markers. International Journal of Molecular Sciences, 13(3), 3900-3915. https://doi.org/10.3390/ijms13033900

Article Metrics

Back to TopTop