About this Journal Submit a Manuscript Table of Contents
International Journal of Forestry Research

Volume 2014 (2014), Article ID 305286, 8 pages

Research Article

Using Random Amplified Polymorphic DNA to Assess Genetic Diversity and Structure of Natural Calophyllum brasiliense (Clusiaceae) Populations in Riparian Forests

1Departamento de Engenharia Florestal, Universidade Federal Rural do Rio de Janeiro, 23890-000 Rio de Janeiro, RN, Brazil

2Departamento de Engenharia Florestal, Universidade de Brasília, 70919-900 Brasília, DF, Brazil

3Unidade Acadêmica Especializada em Ciências Agrárias, Universidade Federal do Rio Grande do Norte, 59280-000 Macaíba, RN, Brazil

4Klabin, Avenida Brasil 26, Bairro Harmonia, 84275-000 Telêmaco Borba, PR, Brazil

5Embrapa Florestas, 74001-970 Goiânia, GO, Brazil

6Departamento de Ciências Florestais, Universidade Federal de Lavras, 37200-000 Lavras, MG, Brazil

Received 28 April 2014; Revised 26 July 2014; Accepted 28 July 2014; Published 21 August 2014

Academic Editor: Robin Reich

Copyright © 2014 Evânia Galvão Mendonça 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.


The objective of this study was to assess the genetic variability in two natural populations of Calophyllum brasiliense located along two different rivers in the state of Minas Gerais, Brazil, using RAPD molecular markers. Eighty-two polymorphic fragments were amplified using 27 primers. The values obtained for Shannon index () were 0.513 and 0.530 for the populations located on the margins of the Rio Grande and Rio das Mortes, respectively, demonstrating the high genetic diversity in the studied populations. Nei’s genetic diversity () was 0.341 for the Rio Grande population and 0.357 for the Rio das Mortes population. These results were not significantly different between populations and suggest a large proportion of heterozygote individuals within both populations. AMOVA showed that 70.42% of the genetic variability is found within populations and 29.58% is found among populations (). The analysis of kinship coefficients detected the existence of family structures in both populations. Average kinship coefficients between neighboring individuals were 0.053 () in Rio das Mortes and 0.040 () in Rio Grande. This could be due to restricted pollen and seed dispersal and the history of anthropogenic disturbance in the area. These factors are likely to contribute to the relatedness observed among these genotypes.

1. Introduction

Of all the ecosystems that constitute the Brazilian Semideciduous Forests, riparian forests contribute significantly to the conservation of biodiversity, mainly due to the relationship between ecological corridors and ecosystem functioning [1]. Because of their interconnectedness, riparian forests play a central role in biogeographical and evolutionary shifts as they facilitate seed dispersal [2]. In Brazil, even though riparian forests are protected by law in the Forest Code (Law 4.771 of 1965), these environments suffer the consequences of human activity, mainly resulting from fragmentation and degradation. The reduction of population sizes and the isolation of populations are direct consequences of anthropogenic activities, resulting in inbreeding and genetic drift [3, 4]. Understanding the distribution of genetic variability of tree populations in these areas is fundamental in management programs aimed at conservation and the survival of riparian forest species.

There are a large number of tree species that occur in riparian forests and as such the choice of the target species for genetic and ecological studies is crucial. Among them, Calophyllum brasiliense Camb. (Clusiaceae), a Neotropical tree commonly known as “guanandi,” stands out as a tree species that occurs preferentially in soils with high humidity saturation and that is sensitive to environmental changes and susceptible to local extinction. The species is widely distributed throughout Central and South America. In Brazil, the wide distribution of C. brasiliense ensures its presence in watersheds of several phytoecological regions.

C. brasiliense presents both male and hermaphroditic flowers. The small “pollen flowers” (<1 cm diameter) are mainly pollinated by halictid bees (Halictidae). Its medium sized, spherical seeds (18–23 mm diameter) are dispersed by gravity and by phyllostomid bats [5]. Bats eat the pulp and discard seeds beneath their feeding roosts [6]. Seeds can also be dispersed secondarily by water, being carried and deposited downstream [7]. Besides the fact that C. brasiliense is a visually attractive species, it is also economically important as it is used for animal feed. It has been widely used in traditional medicine and recent studies show that xanthones from C. brasiliense can potentially decrease the oxidative damage produced by some prooxidants [8]. The species is commonly used in reforestation programs for environmental restoration because its fruit are attractive to fauna and it is also used in tree plantations in areas where the soil is permanently waterlogged [9]. Due to its various uses and importance, C. brasiliense is currently being exploited unsustainably which may eventually alter important genetic parameters in the species’ gene pool.

Information about the structure and genetic diversity of populations as well as understanding other related factors is crucial for the development and adoption of more adequate management strategies to ensure its genetic conservation. Genetic markers, like random amplified polymorphic DNA (RAPD), are valuable in studies assessing genetic variation in plant species when using carefully scored RAPD markers [10, 11]. This technique is based on the amplification of a DNA fragment, which occurs when a pair of nucleotides hybridizes in opposite directions from the target sequence; this is achieved through denaturation cycles, primer annealing, and enzyme extension of Taq polymerase [12]. Thus, the objective of this study is to evaluate inter- and intrapopulation genetic variability to provide information that can be used in the maintenance and conservation of C. brasiliense.

2. Materials and Methods

2.1. Sampled Populations

The populations sampled in this study are located along two important rivers in Minas Gerais State: the Rio Grande, in the municipality of Itumirim (21°13′08′′S 44°48′14′′W), and Rio das Mortes, in Bom Sucesso (21°07′20′′S 44°48′22′′W), both at an altitude of 850 m. The riparian fragments included in the study are made up of semideciduous forest types. A total of 120 reproductive trees were sampled throughout the study areas. Young leaves were collected and stored at the Tree Species Genetic Conservation Laboratory at the Federal University of Lavras (UFLA), at −80°C, until genomic DNA extraction. All trees had their geographical coordinates recorded with GPS (Garmin 720) to conduct a spatial analysis of the genetic structure of the studied populations.

2.2. DNA Extraction and PCR-RAPD Reactions

Genomic DNA was extracted from 250 mg of young leaves. In order to equalize the concentrations of DNA to 10 ng·mL−1, genotypes were diluted in a buffer according to the formula , where is the volume of the TE solution (µL), is the sample volume, and is the sample concentration. Subsequently, samples were submitted to PCR-RAPD reactions. The DNA was amplified by the RAPD method and reactions were conducted in a GeneAmp PCR system 9700, using 20 “Operon technologies” primers (California, USA). Each RAPD reaction occurred in a volume of 13 µL containing 32 ng of genomic DNA, 0.2 mM of each DNTp, 4 µM of individual primer, 1 unit of Taq polymerase (Phoneutria), PCR reaction buffer (50 mM of Tris pH 8.3, 20 mM of KCl, 2 mM of MgCl2, 5 μμL−1 of BSA, 0.25% of Ficol 400, 10 mM of tartrazine), and ultrapure water. A PCR 9700 thermocycler was used with the following program: 45 cycles with a denaturation step of 94°C for 15 seconds, annealing of the primer at 42°C for 30 seconds, extension at 72°C for 30 seconds, and a final extension step at 72°C for 5 min.

Amplification products were separated in 1.0% (w/v) agarose gel in TBE buffer, running at 200 V for 2 h, and stained with ethidium bromide; bands were estimated by comparison with a 250-bp ladder (Gibco BRL). Negative controls were consistently used to detect amplification of bands from contaminated DNA. Bands with the same molecular weight and mobility were considered identical. The gels were photographed under ultraviolet light using a BioRad transilluminator.

In order to determine the band frequency in each population, images of the gels were carefully analyzed to construct a phenotype matrix composed of 0 and 1, where 0 corresponds to a band’s absence and 1 its presence. The bands that presented weak coloration and minimal definition were discarded. A locus was considered polymorphic when at least two different bands were detected.

2.3. Data Analysis
2.3.1. Identification of the Optimal Number of Markers

In order to verify that the numbers of polymorphic bands generated by the 20 RAPD primers were sufficient to determine genetic similarity, a bootstrap analysis of 10,000 permutations was conducted using GENES software [13]. The squared sum of deviations in relation to the new samples and the stress value (), indicating the fit between the original matrix and the sampled matrix, was calculated. The number of polymorphic bands is considered precise when the stress value is less than 0.05 [14].

2.3.2. Genetic Structure

To assess genetic variability, the software POPGENE (version 1.32) was used to calculate the statistics of genetic variation considering the use of dominant markers and diploid data [15]. Keeping in mind that the software takes into account the dominance of the alleles, this program assumes that for each estimate of allelic frequencies, the loci are in Hardy-Weinberg equilibrium. The allelic frequencies were estimated from the square root of the number of occurrences of the null genotype (recessive).

Analysis of molecular variance, AMOVA [15], was carried out using the software ARLEQUIN version [16]. The total variance is partitioned into covariance components as they relate to intra- and interindividual and/or population differences [17].

2.3.3. UPGMA

The estimate of genetic similarity between each pair of genotypes was calculated using PC-ORD 4.14 software [17]. The values of genetic distance were calculated based on Jaccard’s coefficient [18]. We used dendrograms of the unweighted pair-group method with arithmetic average (UPGMA) to represent genetic similarities [19]. Group consistency was analyzed by calculating the cophenetic correlation coefficient among matrices based on the arithmetic complement and cophenetic values, known as Mantel test () [20], using NTSYS-pc software, version 2.1 [21]. This test compares two matrices, element by element, resulting in a correlation value () that represents the level of relation among the matrices, with the significance of tested by permutations [22].

2.3.4. Fine-Scale Spatial Genetic Structure (SGS)

Kinship coefficients () were calculated from RAPD data according to the dominant marker estimator described by Hardy [23]. This coefficient is the probability that two genes drawn randomly from two individuals are identical by descent; that is, the coefficient estimates the genetic coancestry between individuals. The estimate of with dominant markers requires the calculation of the inbreeding coefficient, [23]. We estimated based on allozymic data [24]. The inbreeding coefficient () was −0.046 for the Rio das Mortes and 0.311 for the Rio Grande. Spatial genetic patterns were further quantified by the statistic [25] for each population based on kinship analyses, using the program SPAGeDi v. 1.2g [26]. The value was estimated as , where is the regression slope and is the mean kinship coefficient between individuals belonging to the first distance class. A jackknife procedure (across loci) was used to estimate standard errors for each distance class and 1,000 randomizations of spatial location were conducted to test for overall spatial structure under the null hypothesis ().

3. Results

3.1. Profile of RAPD Markers

Of the total 27 primers selected, interpretation was possible for only 20 primers, resulting in a total 91 and 92 markers for the Rio Grande and Rio das Mortes populations, respectively (Table 1). Each primer revealed from 1.0 to 7.0 polymorphic loci, producing an average of 5.2 polymorphic loci per primer.

Table 1: Primers, nucleotide sequence, total number of loci (NL), and number of polymorphic loci (NPL) for two analyzed Calophyllum brasiliense populations.
3.2. Evaluating the Optimal Number of Markers

In order to evaluate if the number of amplified polymorphic fragments was able to represent the species’ existent genetic diversity, we conducted a bootstrap resampling analysis. In this analysis, we used the 92 generated fragments from the Rio Grande population (Figure 1) and 91 fragments generated for the Rio das Mortes population (Figure 2). The optimal number of bands was 81 in both populations. With this number of fragments, the correlation’s magnitude value was very close to the maximum for the trees sampled in Rio Grande () and in Rio das Mortes ().

Figure 1: Resampling analysis (bootstrap) and the correlation () with number of fragments from the Rio Grande population. SQd is the squared sum of the deviations in relation to the new samples; is the stress value.
Figure 2: Resampling analysis (bootstrap) and the correlation () with number of fragments from the Rio das Mortes population. SQd is the squared sum of the deviations in relation to the new samples; is the stress value.
3.3. Interpopulation Genetic Diversity

Nei’s genetic diversity () was 0.341 for the Rio Grande population and 0.357 for the Rio das Mortes population, and they were not significantly different (Table 2). The estimated values for the Shannon index () suggest high genetic diversity within and among C. brasiliense populations, with a value of 0.513 for Rio Grande, 0.530 for Rio das Mortes, and 0.588 for the entire population.

Table 2: Genetic diversity and Shannon index for two Calophyllum brasiliense populations.
3.4. Population Genetic Structure

AMOVA showed that 70.42% of the genetic variability of C. brasiliense is found within populations and 29.58% among populations, with an value of 0.2958 (Table 3).

Table 3: Analysis of molecular variance (AMOVA) for 120 individuals sampled from two populations of Calophyllum brasiliense using RAPD molecular markers.
3.5. Genetic Similarity

We observed that only individuals 22, 23, 37, and 53 sampled from the Rio Grande population were grouped with individuals from the Rio das Mortes population (Figure 3). This result underscores the level of divergence among trees sampled in both studied populations. The results confirmed that the polymorphic bands obtained in these two populations, using 20 RAPD primers, were able to assess the genetic similarity of C. brasiliense trees. We can therefore assume that the grouping of the individuals based on the cophenetic correlation is highly reliable (Figure 3).

Figure 3: Hierarchical cluster analysis of Calophyllum brasiliense genotypes in populations from Rio Grande (RG) and Rio das Mortes (RM) using UPGMA, based on data obtained with RAPD. RC: cophenetic correlation.
3.6. Fine-Scale Spatial Genetic Structure

Kinship coefficient analysis detected the existence of family structures in both studied populations (Figure 4). The average kinship coefficient between neighboring individuals (represented by ) was 0.053 () in the Rio das Mortes population and 0.040 () in the Rio Grande population, indicating that neighboring individuals have a higher level of genetic relatedness than random pairs of individuals. Negative values of occurred with increased distances between individuals for both populations. Overall slopes () of two correlograms were significantly different from the null hypothesis of no spatial genetic structure ():  (s.e.) for Rio das Mortes and  (s.e.) for Rio Grande. values revealed strong patterns of SGS for the total sample and they were similar for both populations ( for Rio das Mortes and for Rio Grande).

Figure 4: Correlograms of average coancestry coefficients () per distance class in two Calophyllum brasiliense populations. Confidence intervals around each -value were obtained using a jackknife procedure across all loci.

4. Discussion

4.1. Genetic Diversity

The number of RAPD markers used to evaluate genetic variability in plants is highly variable and not dependent on the species’ level of domestication. Pigato and Lopes [27] when evaluating the genetic variability of four generations of Eucalyptus urophylla identified 86 amplified fragments. Zimback et al. [28] detected 72 polymorphic fragments when analyzing the genetic diversity of Trichilia pallida, a Meliaceae climax species. Wang et al. [29], investigating the genetic structure of Neolitsea sericea, found 50 polymorphic fragments. Torezan et al. [30], studying Aspidosperma polyneuron using RAPD markers, encountered 98 polymorphic fragments, while Phong et al. [31], in studying the genetic variability of Dalbergia oliveri, a leguminous tree of the Fabaceae family, found 24 polymorphic fragments when using 29 RAPD primers. Considering the number of polymorphic fragments described in the many published studies on natural tree species populations, the number of fragments obtained in this study of 92 and 91 (Table 2) can be considered satisfactory, enabling an accurate evaluation of genetic diversity for this species.

The sufficiency of the number of polymorphic bands for genetic similarity, calculated by resampling, showed low stress values (), 0.045 and 0.040 at Rio Grande and Rio das Mortes, respectively (Figures 1 and 2). According to Kruskal [14], stress values lower than 0.05 are an excellent indicator of precision. Therefore, the values obtained in our study of 0.045 for Rio Grande and 0.040 for Rio das Mortes, respectively, validate the results obtained for number of fragments.

The interpopulation genetic divergence indexes obtained in our analysis demonstrate that the values of He (0.341 and 0.357) were similar to the values obtained in other studies of natural tree species populations. Torezan et al. [30], studying six natural populations of Aspidosperma polyneuron, obtained an He of 0.299 for adults and seedlings. Estopa et al. [32] found values ranging from 0.299 to 0.333 in Eremanthus elytropappus, while Zimback et al. [28] found a variation of 0.27 to 0.33 for when studying Trichilia pallida. The high genetic diversity within natural populations of C. brasiliense is consistent with other studies of tree species when comparing the values of Shannon’s index (), where similar methodologies were used. Lacerda et al. [33] observed values ranging from 0.301 to 0.367 for Plathymenia reticulada. Torezan et al. [30] obtained an value equal to 0.410 for adult individuals of Aspidosperma polyneuron, and Zang et al. [34] found diversity values of 0.4382 using RADP markers in natural populations of Larix gmelinii. Comparatively, these studies suggest that there is strong evidence that C. brasiliense presents high levels of genetic variability with values of 0.341 for the Rio Grande population and 0.357 for the Rio das Mortes population. This can be related to factors such as insect-pollinated flowers and a predominantly outcrossed mating system [24]. In the study area, however, populations of C. brasiliense have been greatly reduced and disturbed; thus, the distribution of genetic variability is likely concentrated in small groups of individuals.

4.2. Population Differentiation and Fine-Scale Spatial Genetic Structure

Gillies et al. [35], studying the genetic diversity of Swietenia macrophylla using RAPD markers, observed that 87.43% of the genetic variability is found within populations. Lacerda et al. [33] observed that 12.3% of the genetic variation of Plathymenia reticula is related to differences among populations, based on AMOVA. Trindade and Chaves [36], analyzing the genetic structure of natural populations of Eugenia dysenterica, observed that 8.63% of the genetic variability is found among populations. Zucchi et al. [37], analyzing ten populations of Eugenia dysenterica, found 27.03% of variability among populations and 72.97% within populations. Cardoso et al. [38], studying five remaining populations of Caesalpinia echinata, observed that 28.4% of the genetic variability is due to geographical differences between groups, with 29.6% of the variability occurring among groups and 42% among individuals. The results obtained in our study are consistent with that expected for tree species, indicating that the studied populations of C. brasiliense exhibit considerable levels of genetic variability within populations (70.42%) and high levels of divergence (29.58%). The high genetic divergence among populations may be related to the spatial aggregation of trees and geographical distance between the studied populations (approximately 11 km). Although seeds are generally dispersed by bats [39], the genetic differentiation found in this study may be an indication of restricted seed dispersal.

C. brasiliense can present an aggregated distribution depending on the context in which it occurs, for example, small rock outcroppings. Moreover, the landscape in the study region is highly disturbed and as such the remaining C. brasiliense populations are restricted to certain areas. Therefore, the current gene flow in such habitats can be limited. According to Wright [40], values of genetic differentiation among populations ranging between 0.15 and 0.05 indicate moderate differentiation and between 0.15 and 0.25 elevated differentiation, and values greater than 0.25 can be considered very high differentiation. Thus, the genetic differentiation between the two C. brasiliense populations can be considered to be very high as our analysis produced a value of 0.29.

In this study, the observation of significantly positive kinship coefficients over a short distance indicates a clear family structure, that is, genetic cluster, within the two populations. A similar pattern between populations was found based on the correlograms. The significant SGS found in both populations suggests restricted seed dispersal of C. brasiliense, a primary factor responsible for the observed SGS of this species [41]. The significant genetic clustering at short distances in both populations likely reflects seedling recruitment around the maternal plant [42, 43]. According to Wright [44], if seed dispersal is localized at the scale of investigation, it will result in spatial clustering of genetically related individuals and the development of significant SGS.

The Sp statistic is a robust method to compare fine-scale genetic patterns among plant species, independent of sampling strategies, life histories, and population densities [25]. In comparison to mean values for other species, the estimated Sp values in our study (average ) are within the range of those reported in predominantly mixed mating (), pollen animal-dispersed (), and seed gravity-dispersed () species [25].

4.3. Implications for Genetic Conservation

Appropriate conservation and management strategies are necessary in order to ensure the long-term genetic variability of C. brasiliense. The distribution of genetic variation across the landscape is a prime factor to consider in the conservation and management of natural plant populations [45]. Our results show the existence of family structures and a detailed depiction of fine-scale genetic variation patterns. This spatial information can be used to determine future sampling strategies that ensure appropriate distances between sampled genotypes. For example, sampling methods for seed banks (ex situ conservation) and definitions of the size of an area for in situ conservation can benefit greatly from this information. Based on the correlogram of average in local C. brasiliense populations, individuals located within an area of 535 m (at Rio das Mortes) and 206 m (at Rio Grande) are likely to be genetically similar. Therefore, to optimize sampling design and limit redundant genotypes, only individuals located at greater distances than 535 m or 206 m should be collected. Furthermore, samples should be collected from both populations since our results show that they are significantly different.

Conflict of Interests

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


The authors would like to thank the Fundação de Amparo à Pesquisa do Estado de Minas Gerais (FAPEMIG) for financial support. The authors also thank Dr. Evelyn Nimmo for editing the English of the paper.


  1. G. Viegas, C. Stenert, U. H. Schulz, and L. Maltchik, “Dung beetle communities as biological indicators of riparian forest widths in southern Brazil,” Ecological Indicators, vol. 36, pp. 703–710, 2014.
  2. E. Imbert and F. Lefèvre, “Dispersal and gene flow of Populus nigra (Salicaceae) along a dynamic river system,” Journal of Ecology, vol. 91, no. 3, pp. 447–456, 2003. View at Publisher · View at Google Scholar · View at Scopus
  3. P. R. Aldrich and J. L. Hamrick, “Reproductive dominance of pasture trees in a fragmented tropical forest mosaic,” Science, vol. 281, no. 5373, pp. 103–105, 1998. View at Publisher · View at Google Scholar · View at Scopus
  4. D. Couvet, “Deleterious effects of restricted gene flow in fragmented populations,” Conservation Biology, vol. 16, no. 2, pp. 369–376, 2002. View at Publisher · View at Google Scholar · View at Scopus
  5. E. Fischer and F. A. M. dos Santos, “Demography, phenology and sex of Calophyllum brasiliense (Clusiaceae) trees in the Atlantic forest,” Journal of Tropical Ecology, vol. 17, no. 6, pp. 903–909, 2001. View at Scopus
  6. I. Sazima, W. A. Fischer, M. Sazima, and E. A. Fischer, “The fruit bat Artibeus lituratus as a forest and city dweller,” Ciência e Cultura, vol. 46, pp. 164–168, 1994.
  7. M. C. M. Marques and C. A. Joly, “Germinação e crescimento de Calophyllum brasiliense (Clusiaceae), uma espécie típica de florestas inundadas,” Acta Botânica Brasílica, vol. 14, pp. 113–120, 2000.
  8. T. Blanco-Ayala, R. Lugo-Huitrón, E. M. Serrano-López et al., “Antioxidant properties of xanthones from Calophyllum brasiliense: prevention of oxidative damage induced by FeSO4,” BMC Complementary and Alternative Medicine, vol. 13, p. 262, 2013.
  9. R. B. Torres, L. A. F. Matthes, R. R. Rodrigues, and H. F. Leitão-Filho, “Espécies florestais nativas para plantio em áreas de brejo,” O Agronômico, vol. 44, pp. 13–16, 1992.
  10. M. A. González-Pérez, P. A. Sosa, and F. J. Batista, “Genetic variation and conservation of the endangered endemic Anagyris latifolia Brouss. ex Willd. (Leguminosae) from the Canary Islands,” Plant Systematics and Evolution, vol. 279, no. 1–4, pp. 59–68, 2009. View at Publisher · View at Google Scholar · View at Scopus
  11. Ö. Özbek and A. Kara, “Genetic variation in natural populations of Capparis from Turkey, as revealed by RAPD analysis,” Plant Systematics and Evolution, vol. 299, pp. 1911–1933, 2013. View at Publisher · View at Google Scholar · View at Scopus
  12. J. G. K. Williams, A. R. Kubelik, K. J. Livak, J. A. Rafalski, and S. V. Tingey, “DNA polymorphisms amplified by arbitrary primers are useful as genetic markers,” Nucleic Acids Research, vol. 18, no. 22, pp. 6531–6535, 1990. View at Publisher · View at Google Scholar · View at Scopus
  13. C. D. Cruz, Programa GENES—versão Windows, UFV, Viçosa, Minas Gerais, Brazil, 2001.
  14. J. B. Kruskal, “Multidimensional scaling by optimizing goodness of fit to a nonmetric hypothesis,” Psychometrika, vol. 29, pp. 1–27, 1964. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  15. A. Young, D. Boshier, and T. Boyle, Forest Conservation Genetics, CSIRO Publishing, Melbourne, Australia, 2000.
  16. D. S. Schineider, L. Roesslli, and L. Excoffier, Arlequin ver: A Software for Population Genetics Data Analysis, Genetics and Biometry Laboratory /University of Geneva, Geneva, Switzerland, 2000.
  17. B. S. Weir, Genetic Data Analysis: 2—Methods for Discrete Population Genetic Data, North Carolina State University/Sinauer Associates, Suderland, UK, 1996.
  18. M. R. Anderberg, Cluster Analysis for Applications, Academic Press, New York, NY, USA, 1973. View at MathSciNet
  19. P. H. A. Sneath and R. R. Sokal, Numerical Taxonomy: The Principles and Practice of Numerical Classification, W. H. Freeman, San Francisco, Calif, USA, 1973. View at MathSciNet
  20. B. F. J. Manly, Randomization, Bootstrap and Monte Carlo Methods in Biology, Chapman & Hall, London, UK, 2nd edition, 1997. View at MathSciNet
  21. F. J. Rohlf, “Numerical taxonomy and multivariate analysis system,” New York, NY, USA, version 1.70, 470 pages , 1992.
  22. N. Mantel, “The detection of disease clustering and a generalized regression approach,” Cancer Research, vol. 27, no. 2, pp. 209–220, 1967. View at Scopus
  23. O. J. Hardy, “Estimation of pairwise relatedness between individuals and characterization of isolation-by-distance processes using dominant genetic markers,” Molecular Ecology, vol. 12, no. 6, pp. 1577–1588, 2003. View at Publisher · View at Google Scholar · View at Scopus
  24. A. M. Souza, Estrutura genética de populações naturais de Calophyllum brasiliense Camb. na bacia do alto Rio Grande, UFLA, Tese (Doutorado), Lavras, Brazil, 2006.
  25. X. Vekemans and O. J. Hardy, “New insights from fine-scale spatial genetic structure analyses in plant populations,” Molecular Ecology, vol. 13, no. 4, pp. 921–935, 2004. View at Publisher · View at Google Scholar · View at Scopus
  26. O. J. Hardy and X. Vekemans, “SPAGeDI: a versatile computer program to analyse spatial genetic structure at the individual or population levels,” Molecular Ecology Notes, vol. 2, no. 4, pp. 618–620, 2002. View at Publisher · View at Google Scholar · View at Scopus
  27. S. M. P. C. Pigato and C. R. Lopes, “Avaliação da variabilidade genética em quatro gerações de Eucalyptus urophylla S.T. Blake por meio do marcador molecular RAPD,” Scientia Florestalis, vol. 60, pp. 119–113, 2001.
  28. L. Zimback, E. S. Mori, P. Y. Kageyama, R. F. A. Veiga, and J. R. S. Mello Júnior, “Estrutura genética de populações de Trichillia pallida Swartz (Meliaceae) por marcadores RAPD,” Scientia Florestalis, vol. 65, pp. 114–119, 2004.
  29. Z. Wang, S. An, H. Liu, X. Leng, J. Zheng, and Y. Liu, “Genetic structure of the endangered plant Neolitsea sericea (Lauraceae) from the Zhoushan archipelago using RAPD markers,” Annals of Botany, vol. 95, no. 2, pp. 305–313, 2005. View at Publisher · View at Google Scholar · View at Scopus
  30. J. M. D. Torezan, R. F. de Souza, P. M. Ruas, C. de Fátima Ruas, E. H. Camargo, and A. L. L. Vanzela, “Genetic variability of pre and post-fragmentation cohorts of Aspidosperma polyneuron Muell. Arg. (Apocynaceae),” Brazilian Archives of Biology and Technology, vol. 48, no. 2, pp. 171–180, 2005. View at Publisher · View at Google Scholar · View at Scopus
  31. D. T. Phong, V. T. Hien, T. T. Thanh, and D. V. Tang, “Comparison of RAPD and ISSR markers for assessment of genetic diversity among endangered rare Dalbergia oliveri (Fabaceae) genotypes in Vietnam,” Genetic and Molecular Research, vol. 10, no. 4, pp. 2382–2393, 2011. View at Scopus
  32. R. A. Estopa, A. M. Souza, C. O. M. Moura, M. C. G. Botrel, E. G. Mendonça, and D. Carvalho, “Diversidade genética em populações naturais de candeia (Eremanthus erytropapus (DC.) MacLeish),” Scientia Florestalis, vol. 70, pp. 97–106, 2006.
  33. D. R. Lacerda, M. D. P. Acedo, J. P. L. Filho, and M. B. Lovato, “Genetic diversity and structure of natural populations of Plathymenia reticulata (Mimosoideae), a tropical tree from the Brazilian Cerrado,” Molecular Ecology, vol. 10, no. 5, pp. 1143–1152, 2001. View at Publisher · View at Google Scholar · View at Scopus
  34. L. Zhang, H. G. Zhang, and X. F. Li, “Analysis of genetic diversity in Larix gmelinii (Pinaceae) with RAPD and ISSR markers,” Genetics and Molecular Research, vol. 12, no. 1, pp. 196–207, 2013. View at Publisher · View at Google Scholar · View at Scopus
  35. A. C. M. Gillies, C. Navarro, A. J. Lowe et al., “Genetic diversity in mesoamerican populations of mahogany (Swietenia macrophylla), assessed using RAPDs,” Heredity, vol. 83, no. 6, pp. 722–732, 1999. View at Publisher · View at Google Scholar · View at Scopus
  36. M. D. G. Trindade and L. J. Chaves, “Genetic structure of natural Eugenia dysenterica DC (Myrtaceae) populations in northeastern Goiás, Brazil, accessed by morphological traits and RAPD markers,” Genetics and Molecular Biology, vol. 28, no. 3, pp. 407–413, 2005. View at Publisher · View at Google Scholar · View at Scopus
  37. M. I. Zucchi, J. B. Pinheiro, L. J. Chaves et al., “Genetic structure and gene flow of Eugenia dysenterica natural populations,” Pesquisa Agropecuaria Brasileira, vol. 40, no. 10, pp. 975–980, 2005. View at Publisher · View at Google Scholar · View at Scopus
  38. M. A. Cardoso, J. Provan, W. Powell, P. C. G. Ferreira, and D. E. de Oliveira, “High genetic differentiation among remnant populations of the endangered Caesalpinia echinata Lam. (Leguminosae-Caesalpinioideae),” Molecular Ecology, vol. 7, no. 5, pp. 601–608, 1998. View at Publisher · View at Google Scholar · View at Scopus
  39. M. A. R. Mello, N. O. Leiner, P. R. Guimarães Jr., and P. Jordano, “Size-based fruit selection of Calophyllum brasiliense (Clusiaceae) by bats of the genus Artibeus (Phyllostomidae) in a Restinga area, southeastern Brazil,” Acta Chiropterologica, vol. 7, no. 1, pp. 179–182, 2005. View at Publisher · View at Google Scholar · View at Scopus
  40. S. Wright, Evolution and the Genetics of Populations: Variability within and among Natural Populations, vol. 4, University of Chicago Press, Chicago, Ill, USA, 1978.
  41. J. L. Hamrick and J. D. Nason, “Consequences of dispersal in plants,” in Population Dynamics in Ecological Space and Time, O. E. Rhodes Jr., R. K. Chesser, and M. H. Smith, Eds., pp. 203–236, University of Chicago Press, Chicago, Ill, USA, 1996.
  42. B. A. Loiselle, V. L. Sork, J. Nason, and C. Graham, “Spatial genetic structure of a tropical understory shrub, Psychotria officinalis (Rubiaceae),” The American Journal of Botany, vol. 82, no. 11, pp. 1420–1425, 1995. View at Scopus
  43. S. Kalisz, J. D. Nason, F. M. Hanzawa, and S. J. Tonsor, “Spatial population genetic structure in Trillium grandiflorum: the roles of dispersal, mating, history, and selection,” Evolution, vol. 55, no. 8, pp. 1560–1568, 2001. View at Scopus
  44. S. Wright, “Isolation by distance,” Genetics, vol. 28, pp. 114–138, 1943.
  45. B. K. Epperson, “Spatial patterns of genetic variation within plant populations,” in Plant Population Genetics, Breeding, and Genetic Resources, A. H. D. Brown, M. T. Clegg, A. L. Kahler, and B. S. Weir, Eds., pp. 229–253, Sinauer Associates, Sunderland, Mass, USA, 1990.