Research Article

Genetic and phenotypic association of the carnauba palm tree evaluated by inter-simple sequence repeat and biometric traits

Received: May 22, 2018
Accepted: August 03, 2018
Published: August 03, 2018
Genet.Mol.Res. 17(3):


The carnauba palm Copernicia prunifera is the third most important non-timber forest species in Brazil; it is mainly known for the production of carnauba wax. We examined intrapopulation genetic diversity and correlated genotypic and phenotypic traits of 28 trees from a coastal location in Rio Grande do Norte state. Phenotypic variables tree height, diameter at breast height, curvature of the stem, fresh mass of fruits and seeds, and length and diameter of fruits and seeds were evaluated. Genotypic parameters were examined based on inter-simple sequence repeat (ISSR) polymorphisms. The eight ISSR primers detected 79 loci, of which 76% were polymorphic. The polymorphic information content ranged between 0.408 and 0.500, the Nei's diversity index h was 0.327, and Shannon index was 0.470. Although there was no significant association between the matrices of genotypic and phenotypic traits, stepwise multiple regression analysis identified ISSR markers that can be used in a marker-assisted breeding program. The ISSR loci 881(2), 880(1), and M1(4) were significantly and positively associated with plant height. The sMRA also identified a positive association between DBH and 842(5) and a negative association between the curvature of the stem and 857(1). The fresh mass of fruits had a significant negative correlation with ISSR locus 840(1). The other traits showed positive and negative associations when the marker(s) of the previous step(s) were included in the succeeding step in sMRA analyses. These results are relevant for breeding programs, with the perspective of obtaining intra specific hybrids that are superior for desirable characters.


Copernicia prunifera; ISSR; genetic diversity; stepwise multiple regression analysis


Studies on genetic diversity are useful for the evaluation, description and classification of genetic resources for prese

rvation (Govindaraj et al., 2015). The results obtained can subsidize strategies for the selection of genotypes that can be used in genetic improvement. In addition, information about the genetic diversity of populations is necessary for the creation of genetic conservation methods (Kagimbo et al., 2018).

Differentiation of genotypes can be quantified from morphological, biochemical or DNA markers (Adhikari et al., 2017). Characterization of genotypes through DNA molecular markers has become popular because DNA markers are stable and will not be confused with environmental, pleiotropic and epistatic effects (Nybom, 2004; Adhikari et al., 2017). Also, the diversity level in DNA markers is greater. In general, the choice among different types of genetic markers depends on the objectives and their applicability in scientific research (Fajardo et al., 2014). DNA markers are unique because of the simplicity of use, cost, type of technology applied, reproducibility, dominance or codominance, and the ability to distinguish between individuals (Adhikari et al., 2017).

Among molecular markers, the ISSRs (inter-simple sequence repeats) are frequently used and have been shown to be efficient. ISSR markers have a high level of polymorphism, with high reproducibility, and they do not require prior knowledge of the genome, making it a relatively low-cost DNA marker (Desai et al., 2015). The use of ISSRs has been shown to be useful in the analysis of the genetic diversity of palm trees, as in Phoenix dactylifera L. (Mirbahar et al. 2016), Elaies guineensis (Chagas et al., 2015), Attalea vitrivir (Santos et al., 2015), Orbignya phalerata Mart. (Viana et al., 2015), and Mauritia flexuosa (Rossi et al., 2014).

Among species with important economic value in Brazil, we highlight Copernicia prunifera (Arecaceae), known in Brazil as carnauba. This species is native to the Brazilian semi-arid region, with predominance in the states of Ceará, Piauí and Rio Grande do Norte (Leitman et al., 2015). Exploitation of C. prunifera is mainly based on the extraction of awaxy powder from the leaves, which is then processed into carnauba wax, with a production value of US$140 million in the years 2015-2016 (IBGE, 2016). Fibers of C. prunifera are also used in artisanal activities, with enormous economic impact (Sousa et al., 2015).

Despite the importance of C. prunifera in the economy of the northeastern region of Brazil, research involving characterization of genotypes is non-existent. Stepwise multiple regression analysis (sMRA) is used to look for associations between phenotypic traits and molecular markers (Khadivi-Khub, 2014; Choudhary et al., 2017). ISSR markers have frequently been utilized in the sMRA approach to evaluate associations different agronomic and commercial traits (Ganopoulos et al., 2011; Ding et al., 2016). Such studies contribute to the adoption of genetic improvement programs and domestication. To this end, we estimated the intra-population genetic diversity of C. prunifera by ISSR markers and examined correlations between phenotypic and genotypic traits among individuals of a local population.

Material and Methods

Study site and sampling

The study was carried out in a planted population at coordinates 5°57’59.14”S and 35°08’34”W in the municipality of Parnamirim, Rio Grande do Norte state, Brazil. The planting was established approximately 25 years ago from seeds of an open pollinated population. The climate of the region is type As, tropical with dry summer and rainy seasons, according to the Köppen classification (Alvares et al., 2013).

Twenty-eight individuals were sampled. Measurements included the height of the plants, DBH (diameter at breast height – approximately 1.3 meters) and curvature of the stem by a qualitative scale (1 - rectilinear, 2 - rectilinear with curvature and 3 - tortuous). Plant height was estimated considering the distance between the soil surface and the last living branch of the plant. The DBH value was obtained through the formula DBH = CBH/π, where CBH is the circumference at breast height. Leaf samples were collected and conditioned in 2 mL plastic tubes containing 2X CTAB (cetyltrimethylammonium bromide).

Ripe fruits were collected from all 11 adult plants that presented fruits during the study period. We collected 15 fruits from each plant, totaling 165 fruits. For each plant, the length (mm), diameter (mm) and fresh mass (g) of the fruits and seeds were recorded with a digital caliper and an electronic scale.

DNA extraction, PCR, and electrophoresis

The DNA was extracted following a standardized protocol for C. prunifera (Vieira et al., 2015). The DNA samples were quantified using an EpochTM spectrophotometer. The PCRs were performed in a Veriti 96-well thermocycler. Twenty-seven ISSR primers were tested, with a length between 13 and 18 nitrogenous bases and with %CG between 44 and 80%. The PCR solution was run in a total volume of 12 μL with Buffer (10 X, IC PHTTM), BSA (1.0 mg.mL-1), MgCl2 (50 mM), dNTP (2.5 mM), primer (2 μM), Taq polymerase (U.μL-1), DNA (50 ng), and ultrapure water. The PCR cycle consisted of denaturation at 94°C for 5 min, followed by 37 cycles at 94°C for 15 s, 47°C for 30 s, and 72°C for 1 min. The process was finalized at 72°C for 7 min and cooled to 4°C.

PCR products were subjected to electrophoresis on a 1.5% agarose gel (v/v) in TAE buffer (1.0 X, EDTA Tris-acetate). We ran the electrophoresis at 110 V and utilized a molecular weight marker (ladder) of 1 kb. The gels were stained with GelRedTM and photographed with ultraviolet light in an E-BoxTM VX2.

Data analysis

Descriptive statistics and genetic diversity

The phenotypic characteristics were analyzed through descriptive statistics, including position and dispersion measures using SPSS 23 software (IBM, 2015). Deviations from the normal distribution were evaluated by skewness (S) and kurtosis (K) values.

The genotypes were evaluated by a binary matrix, according to the presence (1) or absence (0) of the ISSR marker. The polymorphism information content (PIC) was calculated using the formula PIC = equationwhere Pij is the frequency of the allele "j" in the marker "i" (Anderson et al., 1993). The percentage of polymorphic loci (%P), number of alleles observed (na), number of effective alleles (ne), genetic diversity of Nei (h, also called HE) and Shannon Index (I) were calculated using POPGENE 1.32 software (Yeh et al., 1997).

The pairwise genetic similarities between genotypes were evaluated using a Dice coefficient. The clustering procedure was the method UPGMA (Unweighted Pair-Group Method with Arithmetic Mean). The cophenetic correlation (rC) and bootstrap support values were calculated considering 1,000 replicates. This analysis was performed using the software PAST version 3.18 (Hammer et al. 2001). The resulting dendrogram was drawn using TreeView 1.6.6 software (Page, 1996).

Association among traits and stepwise multiple regression analysis

The phenotypic traits were standardized to Z scores prior to obtaining the Euclidean distance matrix in the SPSS software. The genetic distance matrix was obtained based on Nei's standardized genetic distance (Ds) using FreeTree software (Hampl et al., 2001). The correlation between the phenotypic matrices (Euclidean distance) and the genotypic matrix (Ds) was evaluated by the Mantel (rM) test using GenAlex 6.503 software (Peakall and Smouse, 2012). The significance value of the rM was obtained through the Monte Carlo test with 999 permutations.

Association between ISSR markers (as independent variables) and the phenotypic traits (as dependent variables) was estimated through stepwise multiple regression analysis (sMRA) according to Virk et al. (1996), using the software SPSS. The stepping criteria employed to select independent variables for entry and removal was F values with 0.045 and 0.099 probabilities, respectively.


Phenotypic traits

The height and DBH ranged, respectively, from 3.20 to 6.25 m, and 31.21 to 45.22 cm. The mean fruit average length, diameter, and fresh mass were, respectively, 26.12 mm, 16.25 mm and 2.64 g (Table 1). The standard error values (s.e. < 1.0) indicated good accuracy in the sampling of phenotypic traits.

Traits n Minimum Mean (s.e.) Maximum CV (%) S K
Plants’ height (m) 28 3.20 4.69 (0.16) 6.25 17.94 -0.09 -0.66
DBH (cm) 28 31.21 38.85 (0.77) 45.22 10.61 -0.31 -1.06
Fruits (FM) (g) 165 0.61 2.64 (0.13) 7.15 67.44 1.07 -0.18
Fruits (L) (mm) 165 20.20 26.12 (0.23) 33.80 11.52 0.48 -0.25
Fruits (D) (mm) 165 6.60 16.25 (0.21) 22.00 16.65 -0.0003 -0.04
Seeds (FM) (g) 165 0.17 1.25 (0.07) 3.41 66.76 0.80 -0.55
Seeds (L) (mm) 165 11.10 16.22 (0.20) 22.90 15.62 0.33 -0.62
Seeds (D) (mm) 165 4.20 10.55 (0.23) 17.00 27.89 0.15 -1.00

Table 1: Descriptive statistics of the phenotypic traits of Copernicia prunifera.

The CV values indicated that the fresh mass (FM) of fruits and seeds varied more than other traits, with a high CV of 67.44 and 66.76%, respectively. The kurtosis values indicated that all the traits had a platykurtic distribution (K < 0), indicating wide dispersion of the data around the mean.

ISSR markers and genetic diversity

The eight selected ISSR primers detected 79 loci (Table 2), of which 76% were polymorphic. The number of loci ranged from 6 to 13 per primer. The PIC ranged between 0.408 and 0.500.

Primer ISSR Sequence (5’ – 3’) NL PIC
Mean   9.88 0.461

Table 2: Nucleotide sequence of the ISSR primers, number of loci (NL), and the polymorphic information content (PIC) for each primer.

The sampled population had an observed number of alleles na = 1.76 (s.e. 0.05), number of effective alleles ne = 1.59 (s.e. 0.04), Nei’s genetic diversity h = 0.327 (s.e. 0.02), and Shannon index I = 0.470 (s.e. 0.03). Figure 1 presents the clustering analysis generated for all genotypes.

Figure 1: Cluster analysis of genotypes of Copernicia prunifera, based on Dice coefficient similarity. Bootstrap values greater than 15% are shown.The cophenetic correlation was rC = 0.8.

Association among traits and stepwise multiple regression analysis

There was not significant correlation between the Euclidean distance matrix (for plant height, DBH, and stem curvature) and the genetic distance matrix (Mantel test, rM = 0.096; P = 0.151) and between the Euclidean distance matrix (for fruits and seeds) and the genetic distance matrix (rM = 0.006; P = 0.366).

However, that the ISSR loci 881(2), 880(1), and M1(4) were significantly and positively correlated with plant height (sMRA analysis, Table 3). The sMRA also identified a positive association between DBH and 842(5) and a negative association between stem curvature and 857(1). The fresh mass of fruits significantly and negatively correlated with the ISSR locus 840(1). The other traits showed positive and negative associations when the marker(s) of the previous step(s) were included in the succeeding step in the sMRA analyses.

Traits ISSR Markers(Locus) r Adjusted R2 FANOVA(P< 0.05) ß t value P value
Plant height (m) 881(2) 0.547 0.272 10.7 0.481 3.528 0.002
  +880(1) 0.472 0.451 11.7 0.370 2.659 0.014
  +M1(4) 0.465 0.527 10.6 0.308 2.194 0.039
DBH (cm) 842(5) 0.454 0.174 6.5 0.454 2.545 0.017
Stem 857(1) -0.466 0.186 6.9 -0.466 -2.635 0.014
Fruits (FM) (g) 840(1) -0.790 0.577 13.3 -1.097 -15.806 0.000
  +840(5) 0.187 0.917 50.9 0.612 9.102 0.000
  +881(10) 0.018 0.966 86.5 -0.212 -3.330 0.016
Fruits (L) (mm) 873(10) 0.846 0.680 20.1 0.936 17.452 0.000
  +825(5) -0.023 0.843 25.2 -0.551 -10.067 0.000
  +840(4) -0.478 0.931 41.4 -0.272 -5.311 0.003
  +825(8) 0.269 0.980 110.5 0.202 3.950 0.011
Fruits (D) (mm) 880(7) -0.745 0.499 10.0 -1.060 -15.446 0.000
  +880(8) 0.276 0.825 22.3 0.709 10.327 0.000
  +873(7) -0.023 0.966 86.4 0.363 5.477 0.002
Seeds (FM) (g) 840(1) -0.755 0.517 10.6 -1.092 -14.695 0.000
  +840(5) 0.187 0.827 22.5 0.595 8.267 0.000
  +881(10) -0.112 0.961 75.3 -0.343 -5.024 0.002
Seeds (L) (mm) 840(1) -0.736 0.484 9.5 -0.827 -7.585 0.000
  +881(11) -0.550 0.800 19.0 -0.661 -5.945 0.001
  +857(3) 0.113 0.902 28.6 -0.333 -2.878 0.028
Seeds (D) (mm) 840(1) -0.720 0.458 8.6 -1.072 -14.433 0.000
  +840(5) 0.167 0.708 11.9 0.553 7.683 0.000
  +881(10) -0.236 0.961 75.3 -0.466 -6.833 0.000

Table 3: ISSR markers associated with phenotypic traits in Copernicia prunifera as revealed by multiple regression analysis.


Phenotypic traits

Copernicia prunifera typically reaches 7-10 m in height and can reach up to 15 m (Lorenzi, 2010). In the study population, the plants were relatively small, with an average height of 4.69 m. This difference is probably due to environmental factors. The plant height and DBH variables presented S < 0, indicating that tall individuals and with higher DAP predominate in this population.

There was significant variability in the fresh mass of fruit and seed. This is probably due to environmental effect (e.g. availability of water), which is relevant for the production of fleshy fruits (Chenet al., 2017). The mass and the length of the fruits showed positive skewness (S > 0), indicating that fruits with smaller length and mass predominate in this sample. The diameter of the fruit presented S < 0, evidencing asymmetric distribution to the left, meaning a tendency of fruits with a larger diameter in the sample. The seeds presented S > 0, indicating a predominance of seeds with less mass, length and diameter in the sample.

Phenotypic studies including biometric observations of fruits and seeds are useful for understanding the intra-population variation (Vieira and Carvalho, 2009), and can support conservation and genetic improvement programs. Freitas et al. (2015) studied Spondias sp. genotypes and generated valuable information for selection based on biometric parameters of the fruits. Júnior et al. (2011) studied Prunus persica and provided relevant data for selection of superior genotypes regarding fruit quality. Yadav et al. (2017) biometrically characterized fruits and seeds of Annona squamosa and indicated that this approach is useful for tree selection in genetic improvement programs.

Genetic diversity

The PIC demonstrates the efficiency of the ISSR primers used to detect the polymorphism. All the markers were moderately informative, according to the Botstein et al. (1980) classification. Vieira et al. (2015) studied C. prunifera genotypes and found PIC values ranging from 0.057 to 0.444. The PIC observed here (0.46) was higher than in other studies of palm trees using ISSR, including Borassus flabellifer (0.23) (Vinayagam et al., 2009), Phoenix dactylifera (0.32) (Sabir et al., 2014), and Attalea vitrivir (0.45) (Santos et al., 2015).

The h value in this study (0.327) was higher than expected for long-lived perennial and out crossing species (0.25 and 0.27, respectively) (Nybom, 2004). Our results were similar to Pinheiro et al. (2017) that found h = 0.341 and I = 0.505 for adults individuals of C. prunifera, but superior to those of Vieira et al. (2016), who found h = 0.267 and I = 0.427 in a natural population of C. prunifera.

Association among traits and stepwise multiple regression analysis

Silva et al. (2007) also did not identify an association between the biometric characteristics of the fruits and genetic distance of the palm tree Geonoma schottiana. Given that the phenotype is the result of the action of the genotype under the influence of the environment, the phenotypic characteristics of the fruits are probably affected mainly by the availability of water, which is the primary determinant of fruit mass (Chen et al., 2017).

Although we did not find a significant correlation between genotypic and phenotypic traits based on the Mantel test, the sMRA analysis identified ISSR markers that can be used in a marker-assisted breeding program. These results are relevant for breeding programs, with the perspective of obtaining intraspecific hybrids for desirable characters. Integrating phenotypic traits with molecular markers has recently become a useful tool for estimating genetic divergence and consequently select plants for breeding (Choudhary et al., 2017). Additional studies using assemblages of DNA-based markers, biological and agronomic traits, are necessary to create a database for ex situ conservation and domestication of C. prunifera genotypes.

In conclusion, the population of C. prunifera showed considerable genetic diversity among the plants, based on the ISSR markers. The phenotypic traits evaluated were also diverse, especially the fresh mass of fruits and seeds. There was no correlation between the matrices of phenotypic and genotypic traits, indicating significant environmental influence in the set of observations. However, a positive and significant correlation was observed between some ISSR markers and phenotypic traits, indicating that such characteristics can be determined by genetic factors and thus should be considered in a breeding program.

Conflicts of interest

The authors declare no conflict of interest.


This work was supported by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), process no.471099/2012-0. FA Vieira received a Post-doctoral Scholarship Abroad from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), process no.88881.118979/2016-01. We thank LI Santos for laboratory assistance. We acknowledge Dr. Alistair Jump for providing office facilities at the University of Stirling (UK).

About the Authors

Corresponding Author

C.G. Fajardo

[email protected]


  • Adhikari S, Saha S, Biswas A, Rana TS, et al. (2017). Application of molecular markers in plant genome analysis: a review. The Nucleus 60: 283-297.
  • Alvares CA, Stape JL, Sentelhas PC, de Moraes G, et al. (2013). Köppen’s climate classification map for Brazil. Meteorol. Z. 22: 711-728.
  • Anderson JA, Churchill GA, Autrique JE, Tanksley SD, et al. (1993). Optimizing parental selection for genetic linkage maps. Genome 36: 181-186.
  • Botstein D, White RL, Skolnick M, Davis RW (1980). Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am. J. Hum. Genet. 32: 314-331.
  • Chagas KPT, Sousa RF, Fajardo CG, Vieira FA (2015). Seleção de marcadores ISSR e diversidade genética em uma população de Elaies guineensis. Agrária 10: 147-152.
  • Chen SC, Cornwell WK, Zhang HX, Moles AT (2017). Plants show more flesh in the tropics: variation in fruit type along latitudinal and climatic gradients. Ecography 40: 531-538.
  • Choudhary SB, Sharma HK, Kumar AA, Chowdhury I, et al. (2017). Genetic diversity spectrum and marker trait association for agronomic traits in global accessions of Linumusitatissimum L. Ind. Crops Prod. 108: 604-615
  • Desai P, Gajera B, Mankad M, Shah S, et al. (2015). Comparative assessment of genetic diversity among Indian bamboo genotypes using RAPD and ISSR markers. Mol. Biol. Rep. 42: 1265-1273.
  • Ding J, Ruan CJ, Guan Y, Shan JY, et al. (2016). Characterization and identification of ISSR markers associated with oil content in sea buckthorn berries. Genet. Mol. Res. 15: gmr.15038278.
  • Fajardo CG, Vieira FA, Molina WF (2014). Interspecific genetic analysis of orchids in Brazil using molecular markers. Plant. Syst. Evol. 300: 1825-1832.
  • Freitas WES, Almeida MLB, Morais PLD, Silva JR (2015). Different ‘umbu-cajá’ clones fruit quality. Científica 43: 236-240. 2015 v 43 n3 p 236-240
  • Ganopoulos IV, Kazantzis K, Chatzicharisis I, Karayiannis I, et al. (2011). Genetic diversity, structure and fruit trait associations in Greek sweet cherry cultivars using microsatellite based (SSR/ISSR) and morpho-physiological markers. Euphytica 181: 237-251.
  • Govindaraj M, Vetriventhan M, Srinivasan M (2015). Importance of genetic diversity assessment in crop plants and its recent advances: an overview of its analytical perspectives. Genet. Res. Int. 431487:14
  • Hammer Ø, Harper DAT, Ryan PD (2001). Past: paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4(1): 1-9.
  • Hampl V, Pavlícek A, Flegr J (2001). Construction and bootstrap analysis of DNA fingerprinting-based phylogenetic trees with a freeware program FreeTree: Application to trichomonad parasites. Int. J. Syst. Evol. Microbiol. 51: 731-735.
  • IBGE- Instituto Brasileiro de Geografia Estatística. (2016). SIDRA – Banco de tabelas estatísticas. Available at  []. Accessed February 2, 2018.
  • IBM Corp. Released 2015. IBM SPSS Statistics for Windows, Version 23.0. Armonk, NY: IBM Corp. Júnior AW, Bruckner CH, Cantín CM, Sánchez MAM, et al. (2011). Seleção de progênies e genitores de pessegueiro com base nas características dos frutos. Rev. Bras. Fruticultura 33: 170-179. 29452011005000044
  • Khadivi-Khub A (2014). Regression association analysis of fruit traits with molecular markers in cherries. Plant. Syst. Evol. 300:1163-1173.
  • Kagimbo FM, Shimelis H, Sibiya J (2018). Diversity assessment of sweetpotato germplasm collections for yield and yield-related traits in western Tanzania. Acta Agric. Scand. B. 68: 121-129.
  • Leitman P, Soares K, Henderson A, Noblick L, et al. Arecaceae in Lista de Espécies da Flora do Brasil. Jardim Botânico do Rio de Janeiro, 2015. Available at  [] Accessed February 20, 2018.
  • Lorenzi H (2010). Flora brasileira: Arecaceae (palmeiras). Instituto Plantarum de Estudos da Flora. Nova Odessa, 384 p.
  • Mirbahar AA, Khan S, Markhand GS, Kauser N, et al. (2016). DNA fingerprinting of some pakistani date palm (Phoenixdactylifera L.) cultivars using ISSR markers. Pak. J. Bot. 48: 2005-2010.
  • Nybom H (2004). Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants. Mol. Ecol. 13: 1143-1155.
  • Page RDM (1996). TreeView: An application to display phylogenetic trees on personal computers. Comput. Appl. Biosci. 12: 357-358.
  • Peakall R, Smouse PE (2012). GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics 28: 2537-2539.
  • Pinheiro LG, Chagas KPT, Freire ASM, Ferreira MC, et al. (2017). Anthropization as a determinant factor in the genetic structure of Copernicia prunifera (Arecaceae). Genet. Mol. Res.16: 1-14.
  • Rossi FS, Rossi AAB, Dardengo JFE, Brauwers LR, et al. (2014). Diversidade genética em populações naturais de Mauritia flexuosa L. f. (Arecaceae) com uso de marcadores ISSR. Sci. For. 42(104): 631-639.
  • Sabir J, Abo-Aba S, Bafeel S, Edris S, et al. (2014). Characterization of ten date palm (Phoenix dactylifera L.) cultivars from Saudi Arabia using AFLP and ISSR markers. C. R. Biol. 337: 6-18.
  • Santos RRM, Cavallari MM, Pimenta MAS, Abreu AG, et al. (2015). Population genetic structure of Attalea vitrivir Zona  (Arecaceae)  in  fragmented  areas  of  southeast  Brazil.  Genet.  Mol.  Res.  14(2):  6472-648.
  • Silva MS, Vieira FA, Carvalho D (2007). Biometria dos frutos e divergência genética em uma população de Geonoma schottiana Mart. Rev. Bras. Biocienc. 5(S1): 582-584.
  • Sousa RF, Silva RAR, Rocha TGF, Santana JAS, et al. (2015). Etnoecologia e etnobotânica da palmeira carnauba no semiárido brasileiro. Cerne 21: 587-594.
  • Viana JPG, Borges ANC, Lopes ACA, Gomes RL, et al. (2015). Comparison of eight methods of genomic DNA extraction from babassu. Genet. Mol. Res.14(14): 18003-18008.
  • Vieira FA, Carvalho D (2009). Maturação e morfometria dos frutos de Miconia albicans (Swartz) Triana (Melastomataceae) em um remanescente de floresta estacional semidecídua montana em Lavras, MG. Rev. Árvore 33: 1015-1023.
  • Vieira FA, Sousa RF, Silva RAR, Fajardo CG, et al. (2015). Diversidade genética de Copernicia prunifera com o uso de marcadores moleculares ISSR. Agrária 10(4): 525:531. 10.5039/agraria.v10i4a5040
  • Vieira FA, Sousa RF, Fajardo CG, Brandão MM (2016). Increased relatedness among the neighboring plants from seedling to adult stages in carnauba wax palm. Genet. Mol. Res. 15: 1-10.
  • Vinayagam P, Dhandapani J, Raman P, Alagarsamy RK,et al. (2009). Molecular characterization for genetic diversity of Palmyrah (Borassus flabellifer) accessions using inter simple sequence repeat (ISSR) markers. Asian Australas. J. Plant. Sci. Biotechnol. 3(1): 11-5.
  • Virk PS, Ford-Lloyd BV, Jackson MT, Pooni HS, et al. (1996). Predicting quantitative variation within rice germplasm using molecular markers. Heredity 76: 296–304.
  • Yadav V, Singh AK, Singh S, Rao VVA (2017). Variability in custard apple (Annona squamosa) genotypes for quality characters from Gujarat. Indian J.Agr. Sci. 87(12): 1627-1632.
  • Yeh FC, Yang RC, Boyle TBJ, Ye ZH, et al. (1997). Pop-gene, the user-friendly shareware for population genetic analysis molecular biology and biotechnology center. Molecular Biology and Biotechnology Center, University of Alberta, Edmonton, Alberta, Canada.

Full PDF