Scientifica
Volume 2017 (2017), Article ID 3757238, 9 pages
https://doi.org/10.1155/2017/3757238

Research Article

ISSR Marker Based Population Genetic Study of Melocanna baccifera (Roxb.) Kurz: A Commercially Important Bamboo of Manipur, North-East India

Heikrujam Nilkanta, Thoungamba Amom, Leimapokpam Tikendra, Hamidur Rahaman, and Potshangbam Nongdam

Department of Biotechnology, Manipur University, Canchipur, Imphal, Manipur 795003, India

Correspondence should be addressed to Potshangbam Nongdam; purenba@rediffmail.com

Received 6 August 2016; Accepted 4 December 2016; Published 10 January 2017

Academic Editor: Massimo Ciccozzi

Copyright © 2017 Heikrujam Nilkanta et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Melocanna baccifera (Roxb.) Kurz is an economically important bamboo of North-East India experiencing population depletion in its natural habitats. Genetic variation studies were conducted in 7 populations sampled from 5 districts of Manipur using ISSR molecular markers. The investigation was carried out as a primary step towards developing effective conservation strategies for the protection of bamboo germplasm. ISSR marker analysis showed significant level of genetic variation within the populations as revealed by moderately high average values of Nei’s genetic diversity (H 0.1639), Shannon’s diversity index (I 0.2563), percentage of polymorphic bands (PPB 59.18), total genetic variation ( 0.1961), and genetic diversity within population ( 0.1639). The study also divulged a high genetic variation at species level with Shannon’s diversity index (I), Nei’s genetic diversity (H), and percentage of polymorphic band (PPB%) recorded at 0.3218, 0.1939, and 88.37, respectively. Genetic differentiation among the populations () was merely 19.42% leaving 80.58% of genetic variation exhibited within the populations. The low genetic diversity between populations was consistent with AMOVA. The low genetic differentiation among populations coupled with existence of significantly high genetic diversity at species level indicated the urgent necessity of preserving and protecting all the existing natural bamboo populations in the region.

1. Introduction

Bamboos are multipurpose forest tree grasses having more than 1000 different uses [1]. They are extensively used as raw materials in paper and handicraft industries, in house construction, furniture making, water pipes, storage vessels, and other important household items [2, 3]. Melocanna baccifera (Roxb.) Kurz, an evergreen arborescent, non-clump forming bamboo, is a dominant species of Manipur with its multiple uses in making of houses, furniture, handicraft, and other household items [4, 5]. This bamboo also provides important support for local ecosystem and is a priority species involved in protection of local soils and biodiversity of forested regions [6]. Overexploitation of bamboo by locals in its natural reserves may lead to dramatic fall of population. This will result in great environmental degradation due to water loss, soil erosion, and decline in natural biodiversity. Therefore, systematic management of depleting bamboo resource through the adoption of proper conservation strategies with either in situ or ex situ measures is the need of hour in the region. Assessing the level of genetic variation within and among natural bamboo populations is highly crucial for the development of effective conservation methods [7–10]. This is because the ability of a particular plant species to adapt effectively to changing environmental conditions depends on the level of genetic variability it possesses [11, 12]. DNA molecular markers like random amplified polymorphic DNA (RAPD), inter-simple sequence repeats (ISSR), amplified fragment length polymorphism (AFLP), and simple sequence repeat (SSR) have been used as influential molecular tools for determining genetic variation at species or population level in different plants [13–18]. The traditional approaches of employing morphological and vegetative characters for bamboo species identification and genetic variation studies had many shortcomings [19–21]. However, the application of RAPD, ISSR, AFLP, SSR, expressed sequence tag derived-simple sequence repeat (EST-SSR), sequence related amplified polymorphism (SRAP), and restriction fragment length polymorphism (RFLP) has enabled successful investigation of genetic variability in different bamboos [21–29]. ISSR molecular markers are widely used for population genetic analysis of different plants generating more reliable and reproducible bands than RAPD [30, 31]. They are technically simpler as compared to RFLP, SSR, and AFLP markers as no previous sequence information is required for generating DNA amplification products [32–34]. There have been limited studies on the genetic variation of M. baccifera using ISSR markers and no reports are available on the population genetic studies of the bamboo in Manipur. The present study aimed to investigate the genetic diversity and population genetic structure of M. baccifera in 5 districts of Manipur using ISSR markers. This study could contribute to the development of strategies for effective conservation and sustainable utilization of this bamboo for ecological and economic gains by better understanding the genetic diversity profile at the species and population level.

2. Materials and Methods

2.1. Plant Materials and Population Sampling

A total of 93 individuals of M. baccifera representing 7 populations were collected from different locations spreading across 5 districts of Manipur in North-East India, namely, Bishnupur, Thoubal, Imphal West, Imphal East, and Chandel (Figure 1). The geographical location of each population and its code name and size were presented in Table 1. Fresh leaves were obtained from bamboo plants constituting a particular population. The sample collection was performed from individual plants separated by at least 50 m so as to prevent any possibility of sampling within the same clones. The leaf samples were then stored at −20°C until the DNA extraction was performed.

Table 1: Location of sample collection and sample size for seven populations of M. baccifera in the present study.

Population code
Locality
Sample
size
Latitude (N)
Longitude (E)
Elevation (m)

BISH
Bishnupur district
14
24°37′8.78′′
93°45′42.05′′
828

THBL-A
Areas covering right to Indo-Myanmar road, Thoubal district
14
24°34′7.93′′
94°2′0.86′′
790

IMW-A
Canchipur, Sekmai, Imphal West district
14
24°36′13.03′′
94°4′0.71′′
800

IMPE
Imphal East district
14
24°45′2.48′′
93°56′0.72′′
975

CHDL
Chandel district
14
24°50′9.81′′
93°54′51.22′′
1100

IMW-B
Kangchup areas, Imphal West district
9
24°47′32.38′′
94°2′18.63′′
790

THBL-B
Areas spanning forested regions to the left of Indo-Myanmar road, Thoubal district
14
24°27′10.64′′
94°1′33.04′′
1500

Figure 1: Distribution and location of 7 populations of M. baccifera sampled in 5 districts (blackened area) of Manipur.

2.2. DNA Extraction and ISSR Amplification

The genomic DNA was extracted from the collected leaf samples using the CTAB method with slight modifications [35]. The leaves after being finely ground to fine powder in liquid nitrogen were mixed with freshly prepared CTAB extraction buffer and incubated at 50°C for 15–20 minutes in hot water bath before being subjected to centrifugation at 12000 rpm for 5 minutes. The resultant supernatant was treated with chloroform : isoamyl alcohol (24 : 1) followed by another centrifugation at 13000 rpm for 1-2 minutes. The pellet obtained after 7.5 M ammonium acetate treatment was washed several times with 70% ice-cold ethanol and dried before being resuspended in sterile DNase-free double distilled water. The DNA sample obtained was further incubated at 65°C for 20 minutes to destroy any DNase if present and stored at 4°C for subsequent analysis. DNA quality and quantity were determined through spectrophotometry at 260 and 280 nm, respectively. The purity and integrity were later checked by performing 1.0% agarose gel electrophoresis and comparing the intensity of the resultant bands with 1 kb DNA ladder (Hi-Media). The DNA samples were finally diluted to 50 ng/μL and stored at −20°C for further use. The genetic profile study of M. baccifera was conducted by using 5 ISSR markers, namely, UBC-813, UBC-822, UBC-828, UBC-868, and UBC-878 obtained commercially from the University of British Columbia (Vancouver, Canada). The selected primers showed good, reliable, repetitive, and distinct bands which enabled effective scoring for genetic diversity study within and among the populations.

The DNA amplification mixture of 25 μL contained 25 ng template DNA, 1x PCR buffer, 1.5 mM MgCl₂, 200 mM dNTPs, 1 μM primers (UBC-813, UBC-822, UBC-828, UBC-868, and UBC-878), and 0.6 U Tag DNA polymerase and double distilled sterile water. The PCR components were prepared as master mix for each primer to minimize the pipetting error. The amplification reaction was performed in a thermal cycler (Eppendorf Mastercycler nexus X2) with amplification cycle condition of initial 4 minutes’ strands separation at 94°C followed by 40 cycles of 94°C for 45 secs, 53°C for 1 min, and 72°C for 2 minutes and final extension at 72°C for 10 minutes. The products obtained after PCR amplification were electrophoresed in 2% agarose gel in 0.5x TBE buffer at 100 V for around 3 hours and gel was stained with ethidium bromide (0.5 μg/mL). The fragments after staining were visualized in gel documentation system (Alpha Innotech, USA). A 1000 bp DNA ladder (Hi-media) was used as a size marker for every gel run.

2.3. Data Analysis

Distinct, reproducible, well resolved fragments were scored as present (1) or absent (0) for each ISSR reaction and were displayed as part of a binary matrix. The data matrices obtained were analyzed using POPGENE version 1.31 [36]. Genetic parameters such as percentage of polymorphic band(s) (PPB), observed number of alleles (), the effective number of alleles (), Nei’s genetic diversity (), Shannon’s information index (), total genetic diversity (), genetic diversity within population (), Nei’s genetic differentiation index among populations (), and gene flow estimates between populations () were determined with POPGENE version 1.31 [36]. Genetic divergence between populations of the bamboo species was also investigated using Nie’s unbiased genetic distances and genetic identities [37]. A study on the correlation (Mantel test) between genetic and geographic distance between the bamboo populations was done using TEPGA version 1.3 with 999 permutations [38]. An analysis of molecular variance (AMOVA) was performed to estimate the variance components and their significance levels of genetic variation within and among populations using GenALEx version 6.5 [39]. The unbiased genetic distance was utilized for the construction of a dendrogram using UPGMA (Unweighted Pair Group Arithmetic Mean Method) in POPGENE program version 1.31. In addition, Nei’s genetic distance matrix was used to construct Neighbor-Joining dendrogram for the 93 individuals belonging to 7 different populations of M. baccifera using Mega 5.10 [40]. Principal coordinate analysis (PCoA) was performed using GenALEx version 6.5 [39] to determine relative genetic distance between individuals and to check the consistency of population genetic differentiation as defined by cluster analysis. The Bayesian model base clustering method of STRUCTURE software version 2.2 was employed to evaluate the genetic population structure and detect the gene pools contributing to bamboo germplasm collection [41]. Ancestry model with admixture and correlated allele frequency model were used to determine the posterior probability of the data and a burn in period of 10,000 was set followed by 10,000 Markov Chain Monte Carlo (MCMC) replications. The number of was set from 1 to 8 with 10 independent runs for every value. The results of the STRUCTURE outfile were utilized to determine the optimum value following the simulation method of Evanno et al. [42] employing software STRUCTURE HARVESTER [43].

3. Results

3.1. Genetic Diversity

The five different ISSR markers amplified a total of 93 individuals representing 7 different populations of M. baccifera and generating high level of genetic polymorphism. High genetic variation at species level was observed in the present investigation with the recording of percentage of polymorphic bands (PPB) at 88.37%. The values for , , , and at species level were also recorded at 1.8837, 1.2751, 0.1939, and 0.3218, respectively, showing a relatively high level of genetic diversity (Table 2). However, the genetic differentiation at population level was relatively low as compared to genetic variation evidenced at species level. This was proved by moderate PPB (%) recorded in the range of 51.16% to 69.77% averaging at 59.18%. Similarly, the effective number of alleles () varied from 1.1626 to 1.3923 with an average value of 1.2625, while the observed number of alleles () ranged from 1.5116 to 1.6977 averaging at 1.5940. Nei’s genetic diversity () also extended from 0.1149 to 0.2282 with an average of 0.1639, and Shannon’s information index () spanned from 0.1927 to 0.3400 with an average value of 0.2563. THBL-A and CHDL exhibited highest genetic diversity at population level while the lowest variation was found in IMW-B.

Table 2: Genetic diversity within populations and genetic differentiation parameters of seven populations of M. baccifera.

Population

PPB (%)

BISH
1.5581
1.1995
0.1320
0.2128
55.81

THBL A
1.6977
1.2411
0.1547
0.2504
69.77

IMW-A
1.5349
1.1626
0.1149
0.1927
53.49

IMPE
1.6279
1.2618
0.1678
0.2651
62.79

CHDL
1.6512
1.3923
0.2282
0.3400
65.12

IMW-B
1.5116
1.2914
0.1720
0.2594
51.16

THBL-B
1.5814
1.2888
0.1780
0.2737
58.14

Average
1.5940
1.2625
0.1639
0.2563
59.18

Species level
1.8837
1.2751
0.1939
0.3218
88.37
Total
0.1961
0.1639
0.1942
2.5455

, observed number of alleles; , effective number of alleles; , Nei’s gene diversity; , Shannon’s information indices; PPB, percentage of polymorphic bands; , total genetic diversity; , genetic diversity within populations; , the relative magnitude of genetic differentiation among populations; , estimate of gene flow among populations.

3.2. Genetic Differentiation and Relationship among Populations

The genetic diversity within population () and the total genetic diversity () of the species were recorded at 0.1639 and 0.1961, respectively. The observed genetic differentiation among populations () was 0.1942 demonstrating the presence of 19.42% of genetic variation among the populations. The finding was consistent with the result of analysis of molecular variance (AMOVA) which revealed 22% of the genetic variation existing among populations and 78% within the populations (Figure 2). The differences exhibited between the populations were found to be highly significant () (Table 3). The result was further substantiated by the existence of significantly high level of gene flow among the populations as shown by high estimate of value (2.5455).

Table 3: Analysis of molecular variance (AMOVA) for five ISSR markers among M. baccifera populations.

Source of variations
Degree of freedom
Sum of squares
Mean square
Variance components
% of total variance
value

Among populations
6
136.87
22.81
1.35
22.00
<0.001

Within populations
86
419.92
4.88
4.88
78.00
<0.001

Total
92
556.79

6.23
100

Figure 2: Analysis of molecular variance (AMOVA) revealed 22% of the genetic variation existing among populations and 78% within the populations.

Mantel test was performed to find out if there was any correlation between genetic relatedness between populations and geographic distance separating them. The test result indicated that there was no significant correlation between the genetic and geographic distance (; ) (Figure 3). The genetic distance obtained between the 7 populations of M. baccifera varied from 0.011 between BISH and IMW-A to 0.076 between CHDL and THBL-B, respectively (Table 4). The dendrogram display based on Nei’s [44] unbiased measures of genetic distances showed only one major cluster comprising of BISH, IMW-A, THBL-A, IMPE, CHDL, and THBL-B. The CHDL population tended to separate from the cluster indicating its genetic isolation from other populations in the group (Figure 4). But CHDL exhibited more genetic proximity to THBL-B as compared to other remaining populations in the same cluster with the recording of significantly high value of Nei’s unbiased measures of genetic identity (0.967). The existing greater genetic identity between CHDL and THBL-B in spite of individuals of CHDL being sampled from hilly dominated region may be attributed to geographic closeness of the two populations. The lone IMW-B population diverged from the major cluster in the dendrogram depicting its less genetic relatedness to remaining populations existing in the major cluster. Neighbor-Joining radial tree obtained for the individual samples based on Nei’s genetic distance showed most of the plants belonging to different population origins separated distinctly though some individuals were partly mixed clustered (Figure 5). Individual genotypes belonging to same population were more closely related genetically than those of other populations.

Table 4: Nei’s unbiased measures of genetic identity (above diagonal) and genetic distance (below diagonal) of the seven populations of M. baccifera.

BISH
THBL-A
IMW-A
IMPE
CHDL
IMW-B
THBL-B
Population

0.986
0.989
0.985
0.958
0.948
0.984
BISH

0.014

0.983
0.982
0.967
0.949
0.962
THBL-A

0.011
0.017

0.985
0.942
0.945
0.977
IMW-A

0.015
0.018
0.015

0.952
0.940
0.980
IMPE

0.043
0.033
0.060
0.049

0.936
0.927
CHDL

0.054
0.052
0.056
0.062
0.066

0.939
IMW-B

0.016
0.039
0.023
0.021
0.076
0.063

THBL-B

Figure 3: Correlation between genetic and geographic distance among 7 populations of M. baccifera.

Figure 4: UPGMA dendrogram based on Nei’s (1972) unbiased measures of genetic distance among 7 populations of M. baccifera.

Figure 5: Neighbor-Joining dendrogram clustering pattern for 93 individuals belonging to 7 populations of M. baccifera.

The principal coordinate analysis (PCoA) was conducted to determine the spatial representation of genetic distances observed among individuals of different populations and also to check the consistency of population genetic differentiation as defined by cluster analysis. The two-dimensional PCoA plot showed the first principal coordinate accounting for 16.32% while the second coordinate produced 10.90% of the total genetic variation. The individuals of seven bamboo populations were distributed in the plot in accordance with the cluster analysis. The individuals belonging to CHDL and IMW-B were scattered separately from the rest of the populations (Figure 6). The remaining populations were clustered as one group with individuals more or less intermixed as was depicted in the cluster analysis.

Figure 6: Two-dimensional plot of principal coordinate analysis (PCoA) showing clustering of individual samples belonging to 7 populations of M. baccifera.

The highest delta which determined the best value of was obtained from STRUCTURE HARVESTER and was found to be shown at (Figure 7). The result indicated that all the 93 individuals of seven populations of M. baccifera shared three genetic pools and were assigned into three different clusters (Figure 8). The seven populations displayed significant degree of mixed ancestry as represented in the structure diagram.

Figure 7: Result of Bayesian assignment analysis suggesting as most likely number of clusters as delta value maximum at .

Figure 8: Population structure pattern of 7 populations of M. baccifera generated using STRUCTURE program when .

4. Discussion

The genetic variation analysis conducted in 7 populations of M. baccifera revealed the existence of high genetic diversity within population which was similar to earlier finding in Dendrocalamus membranaceus using ISSR markers [28]. Plants which are long lived, outcrossing, and self-incompatible have higher genetic variation at species level and lower differentiation among populations [11, 45]. The high genetic variation within the population of M. baccifera was expected as being a long lived woody bamboo with long vegetative phase [46]. Self-incompatibility is another important factor in maintaining high genetic variability in population [47] and M. baccifera is known to be self-incompatible and outcrossing bamboo thereby displaying significant genetic diversity within population. Nybom [48] also reported retention of most of genetic variability within populations of long lived and outcrossing plant species when many RAPD and sequence tagged microsatellite site (STMs) based analysis were performed. The plants with high geographical ranges tend to maintain higher genetic diversity than geographically localized species [49]. M. baccifera is known to show diverse distribution pattern along the various geographical locations of Manipur [6]. Genetic diversity within populations is also influenced by many factors such as mating system, population size, extended time period with low number of individuals, genetic drift, and gene flow [50]. High genetic diversity within small populations can also be exhibited if reduction of population size had taken place very recently, especially when it occurred within a generation or two for the concerned species [51]. In such cases the surviving individuals are effectively sampled from the populations that have existed before. The existence of significant genetic variation within the populations of M. baccifera under study may also be due to sudden reduction in population size in short span of time. This presumption might hold true as more accessible individuals growing in the region have been exploited by locals for construction of houses and other household products along with the utilization of their tender shoots as popular vegetable items.

The coefficient of genetic differentiation among populations () was 0.1942 which was higher than the mean value ( 0.073) observed for 121 woody species examined using allozyme markers [52]. Yang et al. [28] also recorded low (0.252) in ISSR analysis of Dendrocalamus membranaceus populations in Yunnan province of China. The low value recorded for M. baccifera indicated low level of genetic differentiation among the populations. Low genetic differentiation between populations had also been confirmed by AMOVA (analysis of molecular variance). The finding was also consistent with high (2.5455) recorded which was much higher than (0.0101) recorded for Dendrocalamus giganteus by Tian et al. [21]. The genetic differentiation between populations is negatively correlated with gene flow [53]. Gene flow between the populations might have been induced by pollens and seeds highly influencing plant evolution [54]. High gene flow ( 2.545) was observed in spite of poor seed dispersal mechanism of M. baccifera. The flower cycle of this bamboo is relatively long and large seeds which fall near the foot of plants are mostly eaten by rats and rodents. However, M. baccifera being one of the economically important bamboos of the region might have undergone human mediated movement of genotypes through local people, farmers, and entrepreneurs. This may lead to enhancement of gene flow among populations and produce overlapping and intermixing of bamboo genotypes from different populations. The easy genetic movement might also be possible as the regions of population sampling in the present study were confined mostly to plain valley area. In fact majority of the populations under study were located in valley districts except for the one population sampled from hill district of Chandel. The STRUCTURE result also revealed high degree of admixture among the 7 populations in three genetic clusters. The existence of great extent of admixture may be attributed to cultivation and utilization pattern of the bamboo species. Geographic isolation is another major factor influencing the genetic differentiation among populations by limiting the amount of gene flow via both pollen and seeds and also through human activities [55]. Theoretically gene flow of more than four migrants per generation is sufficient to prevent genetic differentiation between populations due to drift alone. In our study, the gene flow estimate was high ( 2.5455) ruling out the possibility of inducing the genetic differentiation among populations of M. baccifera due to geographical isolations. The results from the Mantel test also supported the above observation as it showed no correlation between the genetic and geographic distances between different populations. The ultimate goal of conservation is to ensure the continuous survival of population and to maintain their evolutionary potential by preserving natural levels of genetic diversity within and between populations [9, 56, 57]. The natural M. baccifera resources in the selected sites have been exploited by locals for many years leading to loss of plant resources, habitat destruction, and fragmentation. The dramatic fall of bamboo natural population due to increasing demand for its multiutility characters ensures that there is urgent requirement of effective conservation. The population size of M. baccifera sampled from different locations was small requiring population enhancement to maintain standard effective size for high genetic diversity preservation. Genetic drift in very small population might also cause rapid genetic erosion and increased the risk of extinction of the bamboo species. Considering the M. baccifera maintaining significant diversity at species level but low genetic variation among populations and that many native habitats and ecosystems are being destroyed by human interferences, great effort should be made to preserve individuals of the existing populations. Exploitation of bamboos and habitat destruction should be prohibited in the region till effective in situ and ex situ conservation strategies are established. This will prevent further reduction of population size and conserve the overall genetic base and structure of the bamboo species. The in situ efforts to conserve the remaining habitats should be combined with ex situ approach through seed and vegetative propagation of culm cutting, division, and tissue culture technology. These activities should be performed with a view to establishing a new generation of plants both cultivable and wild. A more effective protection strategy may also be framed especially for BISH, IMW-A, IMW-B, and THBL-B which displayed low genetic diversity at population level. Sensitization of locals about the prevailing genetic scenario of alarming bamboo populations and involvement of local communities in framing protection policies will be a highly effective approach to bamboo conservation. The presence of high genetic variation within populations of M. baccifera emphasizes the necessity of preserving and conserving all the existing 7 natural populations and their habitats in Manipur.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of the paper.

Acknowledgments

The authors are thankful to SERB (Science Engineering and Research Board), New Delhi, India, for providing financial support.

References

A. Lewington, Plant for People, Natural History Museum Publications, London, UK, 1990.

M. Das, S. Bhattacharya, P. Singh, T. S. Filgueiras, and A. Pal, “Bamboo taxonomy and diversity in the era of molecular markers,” Advances in Botanical Research, vol. 47, pp. 225–268, 2008.

P. Nongdam and L. Tikendra, “The nutritional facts of bamboo shoots and their usage as important traditional foods of northeast India,” International Scholarly Research Notices, vol. 2014, Article ID 679073, 17 pages, 2014.

H. B. Naithani, “Diversity of bamboo species in North-east India,” in Ecology, Diversity and Conservation of Plants and Ecosystems in India, H. N. Pandey and S. K. Barik, Eds., pp. 312–322, Regency Publication, New Delhi, India, 2006.

J. M. S. Tomar, D. K. Hore, and A. Annadurai, “Bamboo and their conservation in North east India,” Indian Forester, vol. 6, pp. 817–824, 2009.

H. B. Naithani, N. S. Bisht, and S. Singsit, Distribution of Bamboo Species in Manipur, Principle Chief Conservation of Forest, Manipur, India, 1st edition, 2010.

J. Liu, S. Shi, E. Chang, W. Yang, and Z. Jiang, “Genetic diversity of the critically endangered Thuja sutchuenensis revealed by ISSR markers and the implications for conservation,” International Journal of Molecular Sciences, vol. 14, no. 7, pp. 14860–14871, 2013.

Z. Q. Mei, S. Y. Fu, H. Q. Yu et al., “Genetic characterization and authentication of Dimocarpus longan Lour. using an improved RAPD technique,” Genetics and Molecular Research, vol. 13, no. 1, pp. 1447–1455, 2014.

P. Desai, B. Gajera, M. Mankad et al., “Comparative assessment of genetic diversity among Indian bamboo genotypes using RAPD and ISSR markers,” Molecular Biology Reports, vol. 42, no. 8, pp. 1265–1273, 2015.

L. Yeasmin, M. N. Ali, S. Gantait, and S. Chakraborty, “Bamboo: an overview on its genetic diversity and characterization,” 3 Biotech, vol. 5, no. 1, pp. 1–11, 2015.

J. L. Hamrigk and M. J. W. Godt, “Effects of life history traits on genetic diversity in plant species,” Philosophical Transactions of the Royal Society B: Biological Sciences, vol. 351, no. 1345, pp. 1291–1298, 1996.

R. Frankham, “Genetic adaptation to captivity in species conservation programs,” Molecular Ecology, vol. 17, no. 1, pp. 325–333, 2008.

S.-M. Chung, J. E. Staub, and J.-F. Chen, “Molecular phylogeny of Cucumis species as revealed by consensus chloroplast SSR marker length and sequence variation,” Genome, vol. 49, no. 3, pp. 219–229, 2006.

D. Ng'uni, M. Geleta, and T. Bryngelsson, “Genetic diversity in sorghum (Sorghum bicolor (L.) Moench) accessions of Zambia as revealed by simple sequence repeats (SSR),” Hereditas, vol. 148, no. 2, pp. 52–62, 2011.

G. Della Rocca, T. Osmundson, R. Danti et al., “AFLP analyses of california and mediterranean populations of Seiridium cardinale provide insights on its origin, biology and spread pathways,” Forest Pathology, vol. 43, no. 3, pp. 211–221, 2013.

J. Sabir, M. Mutwakel, A. El-Hanafy et al., “Applying molecular tools for improving livestock performance: from DNA markers to next generation sequencing technologies,” Journal of Food, Agriculture and Environment, vol. 12, no. 2, pp. 541–553, 2014.

M. Basahi, “ISSR-based analysis of genetic diversity among sorghum landraces growing in some parts of Saudi Arabia and Yemen,” Comptes Rendus Biologies, vol. 11, no. 338, pp. 723–727, 2015.

J. Gao, W. Zhang, J. Li, H. Long, W. He, and X. Li, “Amplified fragment length polymorphism analysis of the population structure and genetic diversity of Phoebe zhennan (Lauraceae), a native species to China,” Biochemical Systematics and Ecology, vol. 64, pp. 149–155, 2016.

D. Ohrnberger and J. Goerrings, The Bamboos of the World, Elsevier, Odenthal, Germany, 1986.

D. Fan, X. Jiarong, Y. Yuming et al., “Study on flowering phenomenon and its type of bamboo in Yunnan in past fifteen years,” Scientia Silvae Sinica, vol. 3, pp. 57–68, 2000.

B. Tian, H.-Q. Yang, K.-M. Wong, A.-Z. Liu, and Z.-Y. Ruan, “ISSR analysis shows low genetic diversity versus high genetic differentiation for giant bamboo, Dendrocalamus giganteus (Poaceae: Bambusoideae), in China populations,” Genetic Resources and Crop Evolution, vol. 59, no. 5, pp. 901–908, 2012.

E. Friar and G. Kochert, “A study of genetic variation and evolution of Phyllostachys (Bambusoideae: Poaceae) using nuclear restriction fragment length polymorphisms,” Theoretical and Applied Genetics, vol. 89, no. 2-3, pp. 265–270, 1994.

C. Lai and J. Y. Hsiao, “Genetic variation of Phyllostachys pubescens (Bambusoideae, Poaceae) in Taiwan based on DNA polymorphisms,” Botanical Bulletin of Academia Sinica, vol. 38, no. 3, pp. 145–152, 1997.

T. R. Hodkinson, S. A. Renvoize, G. N. Chonghaile, C. M. A. Stapleton, and M. W. Chase, “A comparison of ITS nuclear rDNA sequence data and AFLP markers for phylogenetic studies in Phyllostachys (Bambusoideae, Poaceae),” Journal of Plant Research, vol. 113, no. 3, pp. 259–269, 2000.

N. A. Barkley, M. L. Newman, M. L. Wang, M. W. Hotchkiss, and G. A. Pederson, “Assessment of the genetic diversity and phylogenetic relationships of a temperate bamboo collection by using transferred EST-SSR markers,” Genome, vol. 48, no. 4, pp. 731–737, 2005.

S. Bhattacharya, M. Das, R. Bar, and A. Pal, “Morphological and molecular characterization of bambusa tulda with a note on flowering,” Annals of Botany, vol. 98, no. 3, pp. 529–535, 2006.

R. K. Sharma, P. Gupta, V. Sharma, A. Sood, T. Mohapatra, and P. S. Ahuja, “Evaluation of rice and sugarcane SSR markers for phylogenetic and genetic diversity analyses in bamboo,” Genome, vol. 51, no. 2, pp. 91–103, 2008.

H.-Q. Yang, M.-Y. An, Z.-J. Gu, and B. Tian, “Genetic diversity and differentiation of Dendrocalamus membranaceus (Poaceae: Bambusoideae), a declining bamboo species in Yunnan, China, as based on Inter-Simple Sequence Repeat (ISSR) analysis,” International Journal of Molecular Sciences, vol. 13, no. 4, pp. 4446–4447, 2012.

S. Zhu, T. Liu, Q. Tang, L. Fu, and S. Tang, “Evaluation of bamboo genetic diversity using morphological and SRAP analyses,” Russian Journal of Genetics, vol. 50, no. 3, pp. 267–273, 2014.

T. Nagaoka and Y. Ogihara, “Applicability of inter-simple sequence repeat polymorphisms in wheat for use as DNA markers in comparison to RFLP and RAPD markers,” Theoretical and Applied Genetics, vol. 94, no. 5, pp. 597–602, 1997.

L.-J. Zhang and S.-L. Dai, “Genetic variation within and among populations of Orychophragmus violaceus (Cruciferae) in China as detected by ISSR analysis,” Genetic Resources and Crop Evolution, vol. 57, no. 1, pp. 55–64, 2010.

S. P. Joshi, V. S. Gupta, R. K. Aggarwal, P. K. Ranjekar, and D. S. Brar, “Genetic diversity and phylogenetic relationship as revealed by inter simple sequence repeat (ISSR) polymorphism in the genus Oryza,” Theoretical and Applied Genetics, vol. 100, no. 8, pp. 1311–1320, 2000.

B. Bornet and M. Branchard, “Nonanchored inter simple sequence repeat (ISSR) markers: reproducible and specific tools for genome fingerprinting,” Plant Molecular Biology Reporter, vol. 19, no. 3, pp. 209–215, 2001.

M. X. Wang, H. L. Zhang, D. L. Zhang et al., “Genetic structure of Oryza rufipogon Griff. in China,” Heredity, vol. 101, no. 6, pp. 527–535, 2008.

J. J. Doyle and J. L. Doyle, “A rapid DNA isolation procedure for small quantities of fresh leaf material,” Phytochemical Bulletin, vol. 19, pp. 11–15, 1987.

C. F. Yeh, R. Yang, and T. Boyle, POPGENE VERSION 1.31. Microsoft Windows-Based Freeware for Population Genetics Analysis, University of Alberta and Centre for International Forestry Research, Alberta, Canada, 1999.

M. Nei, “Estimation of average heterozygosity and genetic distance from a small number of individuals,” Genetics, vol. 89, no. 3, pp. 583–590, 1978.

M. P. Miller, Tools for Populations Genetic Analyses (TFPGA) Version 1.3: A Windows Program for the Analysis of Allozymes and Molecular Population Genetic Data, Utah State University, Longan, Utah, USA, 1997.

R. Peakall and P. E. Smouse, “GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update,” Bioinformatics, vol. 28, no. 19, pp. 2537–2539, 2012.

K. Tamura, D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. Kumar, “MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods,” Molecular Biology and Evolution, vol. 28, no. 10, pp. 2731–2739, 2011.

J. K. Pritchard, M. Stephens, and P. Donnelly, “Inference of population structure using multilocus genotype data,” Genetics, vol. 155, no. 2, pp. 945–959, 2000.

G. Evanno, S. Regnaut, and J. Goudet, “Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study,” Molecular Ecology, vol. 14, no. 8, pp. 2611–2620, 2005.

D. A. Earl and B. M. VonHoldt, “Structure Harvester: a website and program for visualizing Structure output and implementing the Evanno Method,” Conservation Genetics Resources, vol. 4, no. 2, pp. 359–361, 2012.

M. Nei, “Genetic distance between populations,” The American Naturalist, vol. 106, no. 949, pp. 283–292, 1972.

G. Zawko, S. L. Krauss, K. W. Dixon, and K. Sivasithamparam, “Conservation genetics of the rare and endangered Leucopogon obtectus (Ericaceae),” Molecular Ecology, vol. 10, no. 10, pp. 2389–2396, 2001.

S. Jeeva, S. Kiruba, H. Lahruaitlung et al., “Flowering of Melocanna baccifera (Bambusaceae) in North-eastern India,” Current Science, vol. 96, pp. 1165–1166, 2009.

E. L. Borba, J. Semir, and G. J. Shepherd, “Self-incompatibility, inbreeding depression and crossing potential in five Brazilian Pleurothallis (Orchidaceae) species,” Annals of Botany, vol. 88, no. 1, pp. 89–99, 2001.

H. Nybom, “Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants,” Molecular Ecology, vol. 13, no. 5, pp. 1143–1155, 2004.

L. Chen, F. Chen, S. He, and L. Ma, “High genetic diversity and small genetic variation among populations of Magnolia wufengensis (Magnoliaceae), revealed by ISSR and SRAP markers,” Electronic Journal of Biotechnology, vol. 17, no. 6, pp. 268–274, 2014.

S. George, J. Sharma, and V. L. Yadon, “Genetic diversity of the endangered and narrow endemic Piperia Yadonii (Orchidaceae) assessed with ISSR polymorphisms,” American Journal of Botany, vol. 96, no. 11, pp. 2022–2030, 2009.

X. Zhao, Y. Ma, W. Sun, X. Wen, and R. Milne, “High genetic diversity and low differentiation of Michelia coriacea (Magnoliaceae), a critically endangered endemic in southeast Yunnan, China,” International Journal of Molecular Sciences, vol. 13, no. 4, pp. 4396–4411, 2012.

J. L. Hamrick, M. J. W. Godt, and S. L. Sherman-Broyles, “Factors influencing levels of genetic diversity in woody plant species,” New Forests, vol. 6, no. 1, pp. 95–124, 1992.

V. Grant, The Evolutionary Process: A Critical Study of Evolutionary Theory, Columbia University Press, New York, NY, USA, 1991.

J. L. Hamrick, “Gene flow and distribution of genetic variation in plant populations,” in Differentiation Patterns in Higher Plants, K. Urbanska, Ed., pp. 53–76, Academic Press, New York, NY, USA, 1987.

M. Pfeifer and G. Jetschke, “Influence of geographical isolation on genetic diversity of Himantoglossum hircinum (orchidaceae),” Folia Geobotanica, vol. 41, no. 1, pp. 3–20, 2006.

Z.-Y. Zhang, Y.-Y. Chen, and D.-Z. Li, “Detection of low genetic variation in a critically endangered Chinese pine, Pinus squamata, using RAPD and ISSR markers,” Biochemical Genetics, vol. 43, no. 5-6, pp. 239–249, 2005.

H. Wu, Z. Li, and H. Huang, “Genetic differentiation among natural populations of Gastrodia elata (Orchidaceae) in Hubei and germplasm assessment of the cultivated populations,” Biodiversity Science, vol. 14, no. 4, pp. 315–326, 2006.