SciELO - Scientific Electronic Library Online

 
vol.69 issue2Mutualistic interaction of aphids and ants in pepper, Capsicum annuum and Capsicum frutenscens (Solanaceae)Phenology of the endangered palm Ceroxylon quindiuense (Arecaceae) along an altitudinal gradient in Colombia author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

Share


Revista de Biología Tropical

On-line version ISSN 0034-7744Print version ISSN 0034-7744

Rev. biol. trop vol.69 n.2 San José Apr./Jun. 2021

http://dx.doi.org/10.15517/rbt.v69i2.44194 

Artículo

Diversity and genetic structure of Spondias tuberosa (Anacardiaceae) accessions based on microsatellite loci

Diversidad y estructura genética de “accesiones” de Spondias tuberosa (Anacardiaceae) basadas en loci de microsatélites

Viviane Nunes-dos-Santos1 
http://orcid.org/0000-0003-1093-204X

Carlos-Antônio Fernandes-Santos2 
http://orcid.org/0000-0002-6932-6805

Viseldo Ribeiro-de-Oliveira2 
http://orcid.org/0000-0002-5301-7637

Antonio-Elton Da-Silva-Costa3 
http://orcid.org/0000-0003-4653-5660

Fabricio-Francisco Santos-da-Silva4 
http://orcid.org/0000-0001-7938-8308

1. Departamento de Biologia, Universidade Estadual de Feira de Santana, Bahia, Brazil; viviane7414@gmail.com

2. Embrapa Semiárido, Caixa Postal 23. 56300-970, Petrolina, Pernambuco, Brazil; carlos-fernandes.santos@embrapa.br (*Correspondence), viseldo.oliveira@embrapa.br

3. Departamento de Agronomia, Universidade Federal Rural do Pernambuco, Recife, Pernambuco, Brazil; antonioelton.agro@gmail.com

4. Universidade Federal do Vale do São Francisco, Núcleo de Ecologia e Monitoramento Ambiental, Petrolina, Pernambuco, Brazil; fabriciofrancisco2006@gmail.com

Abstract

Introduction:

Spondias tuberosa is a tree endemic to the semiarid region of Brazil with fruticulture potential.

Objective:

To estimate the diversity and genetic structure of S. tuberosa accessions from four areas of the semiarid region of Brazil, in order to facilitate conservation genetic resources studies in this species.

Methods:

DNA was extracted, using the CTAB 2x method, from leaf samples of 24 accessions of S. tuberosa available in the germplasm bank at Embrapa Semiárido, Brazil. Ten microsatellite loci were used in this study.

Results:

The UPGMA dendrogram, generated with a Jaccard coefficient similarity matrix, contains four groups at a 0.44 cutoff point. The similarity coefficient ranged from 0.30 to 0.84, indicating great divergence among the accessions. A Bayesian analysis conducted with the software Structure suggests there are two subpopulations, one formed by accessions from the Januária region and another by accessions from the Juazeiro, Uauá and Petrolina regions. The ΦST value of 0.12 for the analysis of molecular variance indicates moderate genetic differentiation among the four populations, suggesting that the genetic variability is moderately structured in function of region.

Conclusions:

Together, the analyses indicate that the genetic diversity of S. tuberosa is not uniformly distributed in the studied regions. Thus, germplasm from a greater number of populations should be collected to increase the germplasm bank genetic diversity of the species.

Key words: umbu tree; AMOVA; SSR

Resumen

Introducción:

Spondias tuberosa es un árbol endémico de la región semiárida de Brasil con potencial frutícola.

Objetivo:

Estimar la diversidad y caracterizar la estructura genética de accesiones de S. tuberosa en cuatro áreas del semiárido brasileño, para así facilitar estudios de conservación de recursos genéticos de esta especie.

Metodología:

El ADN fue extraído utilizando el método CTAB 2x a partir de muestras de hojas de 24 accesiones de S. tuberosa disponibles en el banco de germoplasma de Embrapa Semiárido, Brasil. Diez loci de microsatélites fueron usados en este estudio.

Resultados:

El dendrograma UPGMA generado con una matriz de similitud de coeficientes de Jaccard, formó cuatro grupos con punto de corte en 0.44. El coeficiente de similitud osciló entre 0.30 y 0.84, indicando una gran divergencia entre las accesiones. El análisis Bayesiano realizado en el software Structure sugiere la existencia de dos subpoblaciones, una formada por las accesiones de la región de Januária y otra derivada de las regiones de Juazeiro, Uauá y Petrolina. El valor de ΦST de 0.12 derivado del análisis molecular de la varianza indica moderada variación genética entre las cuatro poblaciones, sugiriendo que la variabilidad genética se estructura moderadamente en función de la región.

Conclusiones:

Los análisis en conjunto indican que la diversidad genética de S. tuberosa no se encuentra distribuida uniformemente en las regiones estudiadas. Por lo tanto, se debe recolectar germoplasma de un mayor número de poblaciones para aumentar la diversidad genética del banco actual de la especie.

Palabras clave: árbol de Umbú; AMOVA; SSR

Introduction

The semiarid region of Brazil has a highly rich flora and the genus Spondias is widely found in this area. Many species of Spondias occur in spontaneous and subspontaneous forms, such as S. mombin L., S. purpurea L., S. cytherea Sonn., Spondias sp., S. tuberosa x S. mombin, and S. tuberosa Arruda) (Pires, 1990). In particular, the umbu tree (S. tuberosa) species has great relevance to the semiarid region due to its production potential and drought resistance. Umbu tree fruits have many uses, for example, they are eaten fresh and used in sweets, juices, jellies and other local products. In addition to its commercial potential, umbu is also notable for its bioactive compounds, such as phenolic, carotenoid, flavonoid and vitamin C compounds that have beneficial health effects (Silva & Alves, 2008).

Agronomic and genetic characterization are important tools for genetic improvement and possible agronomic exploitation. However, there are few studies that focus on the genetic variability of umbu tree. Santos, Rodrigues, and Zucchi (2008) studied the genetic variability of umbu tree in the Brazilian semi-arid region using AFLP markers and found high variability in populations, suggesting that the genetic diversity of umbu tree could be used to improve the species. The Embrapa Semiárido in Brazil maintains an umbu tree germplasm collection, with 80 accessions (Ramos, Queiroz, Romão, & Silva Júnior, 2008).

Molecular markers have been increasingly used due their widely applied in genetic studies (Turchetto-Zolet, Turchetto, Zanella, & Passaia, 2017). Among these biotechnological tools microsatellite markers or single sequence repeats (SSRs) are widely employed due to their ease of use, codominance, multiallelism and high reproducibility. Due to the restriction of SSR developed for species of Spondias sp., studies of marker transferability are common, as in the study of Aguilar-Barajas et al. (2014). Balbino, Martins, Morais, and Almeida (2019) developed 18 polymorphic SSR markers useful for studies of genetic population and conservation and breeding activities. Estimations of genetic diversity parameters, applying SSR, are still rare with umbu tree, mainly among accessions of the Embrapa Semiárido germplasm collection, the most important one.

The objective of the present study was to estimate the genetic diversity and structure of accessions at the Embrapa Semiárido germplasm bank of S. tuberosa from four areas of the semiarid region in Brazil, which will help guide future genetic resource studies of this species.

Materials and methods

Plant material and DNA extraction and quantification: Samples of young and healthy leaves of 24 accessions from four areas of the Brazil semiarid region (Fig. 1) were collected from the umbu tree germplasm bank at Embrapa Semiárido in Petrolina, Pernambuco, Brazil. The DNA was extracted using the CTAB 2X protocol (Doyle & Doyle, 1990) with some modifications. DNA quantification and the integrity were verified in 0.8 % agarose gel, followed by diluting the genomic DNA to 10 ng mL-1.

Fig. 1 Sampling locations of the umbu tree accessions maintained in the germplasm bank at Embrapa Semiárido (BGU), Brazil. Accessions were collected in the following Brazilian municipalities: Petrolina, Pernambuco (PE); Juazeiro and Uauá, Bahia (BA); and Januária, Minas Gerais (MG). 

PCR protocols: The amplification reactions were made using 10 primers, including four developed by Aguilar-Barajas et al. (2014) and six developed by Balbino et al. (2019). The PCRs were adjusted to a final volume of 10 μL containing the following: 1µL of buffer, 2 mM of MgCl2, 0.22 μM of each dNTP, 0.4 μM of each primer (forward and reverse), a unit of Taq DNA polymerase and 10 ng of genomic DNA.

The amplifications were made in a Biometra thermocycler using the program proposed by Aguilar-Barajas et al. (2014): 15 min at 94 ºC, followed by 35 cycles of 30 s at 94 ºC, 1 min 30 s at 59 or 60 ºC and 1 min at 72 ºC, and a final extension of 10 min at 72 ºC. The amplification products were observed using polyacrylamide gel at a concentration of 6 %, according to the methodology described by Costa and Santos (2013), stained with silver nitrate (Creste, Neto, & Figueira, 2001).

Cluster, population structure and AMOVA analysis: the number of base pairs (bp) of each allele was estimated using the inverse mobility method based on the regression of products of known size of a molecular marker with 50 bp (Ludwig Biotec ®). The microsatellites were analyzed for allelic presence (1) and absence (2) to construct a Jaccard index of similarity. A dendrogram was generated using the UPGMA clustering method (unweighted, based on the arithmetic mean). The dendrogram was tested using the cophenetic correlation coefficient. The program NTSYSpc (Rohlf, 2000) was used for these analyses.

The accessions were grouped using the program STRUCTURE 2.3.4 (Pritchard, Stephens, & Donnelly, 2000) and the Markov Monte Carlo chain (MCMC), with 100 000 permutations and 100 000 simulations for cluster inferences. Ten runs were performed for each K value (number of possible clusters). Using STRUCTURE HARVEST (Earl & vonHoldt, 2012), the ΔK value was calculated to detect the probable number of clusters (Evanno, Regnaut, & Goudet, 2005).

The analysis of molecular variance (AMOVA) was conducted by decomposing the total variation of the components between and within populations using the square Euclidean distance (Excoffier, Smouse, & Quattro, 1992). The significance of the genetic parameters was determined by the randomization method (999 permutations). Gene flow (Nm) was estimated by the number of migrants, based on the FST parameter that is analogous to the ΦST, defined as the function of the between-population variance component and the within-population variance component (ΦST = σ2 a/(σ2 a+σ2 b)) (Wright, 1949; Excoffier et al., 1992; Meirmans & Hedrick, 2011). The program GenAlEx 6.5 (Peakall & Smouse, 2006) was used for the AMOVA.

Results

SSR polymorphism: All of the SSR loci had a good amplification pattern in the polyacrylamide gel (6 %). The allelic diversity ranged from two to seven alleles per locus, with an average of 3.5 alleles per locus (Table 1). The polymorphic information content ranged from 0.195 to 0.778 (Table 1). The expected heterozygosity varied from 0.195 to 0.822 and the observed heterozygosity ranged from 0.167 to 0.958.

Table 1 Genetic parameters of ten microsatellite loci evaluated in 24 individuals of Spondias tuberosa 

Locus Primer sequences (5’-3’) Parameters in S. tuberosa
Ta(ºC) NA BS PIC He Ho
SPO4 Forward: CGCTAGTTGTCATTCGCGG Reverse: GCTTAACCTCTGGAAAGTCGC 62 7 548-696 0.778 0.822 0.958
SPO8 Forward: GCAGCAGCCATTTGTGAAC Reverse: CACGTGTTCCCAGTTATGATTTG 60 4 505-529 0.486 0.552 0.458
SPO14 Forward: ACACCAACGTTTGCGGAG Reverse: TCTAGGTAGACAGCGACAAATC 62 3 640-707 0.575 0.662 0.500
SPO18 Forward: TCTATTTGCGTCCAGGTATTTC Reverse: GAATGGGCACGTTCCTTGG 62 3 820-874 0.477 0.580 0.609
TUB78 Forward: TGCTCTGCCTTCCAACATGT Reverse: GTACGTGAGGGACAATGGGG 59 3 558-642 0.388 0.414 0.458
TUB84 Forward: CACCTCCTACGTTACTGCCA Reverse: TCAAACTGGATTCAGGCATGC 59 2 580-613 0.359 0.479 0.333
TUB93 Forward: AGCCTTTTTGAGTCACATGCA Reverse: ACACTGATGGTACGTGAACAA 60 2 570-624 0.195 0.223 0.167
TUB94 Forward: TGTCTGAGGATCGAGACGAGA Reverse: GCACGCGCTTACTTATGTTGG 60 3 689-767 0.522 0.600 0.542
TUB98 Forward: AGCGGAAAAGAATGATGAAGGC Reverse: GTTGGCTTCTTTCTTGCGGC 60 4 647-683 0.474 0.554 0.667
TUB103 Forward: GGAGCAGTGAAACCCCTGAA Reverse: GTCCAGGTCGCCGTAAAGAA 60 4 519-533 0.357 0.392

Ta(°C): annealing temperature; NA: number of alleles; BS: band size; PIC: polymorphic information content; He: expected heterozygosity; Ho: observed heterozygosity.

Cluster analysis: The clustering of the 24 individuals has a cophenetic correlation coefficient of 1.0 (Fig. 2), indicating that the data are reliable and there is good fit between the genetic distances, original matrix and graphic representation. The similarity matrix ranged from 0.115 to 0.842, indicating high variability among the individuals analyzed (Fig. 2).

Fig. 2 Dendrogram representing the SSR genetic divergence among the 24 accessions of S. tuberosa obtained from the UPGMA method, using the complement of the Jaccard index as the average of dissimilarity. Cophenetic value = 1.0. 

Based on the microsatellite markers and the similarity matrix of the umbu tree accessions, it was possible to separate the accessions into four groups, using the average distance of 0.44 as the cutoff point. The groups are formed by the following accessions: group I by BGU56 (Januária); group II by BGU58, BGU59 and BGU62 (Januária); group III by BGU07, BGU32, BGU60 and BGU61 (Juazeiro, Uauá and Januária, respectively); and group IV by 16 accessions, including BGU01, BGU05, BGU16, BGU41, BGU42 (Juazeiro); BGU28, BGU31, BGU35, BGU37 and BGU38 (Uauá), and BGU12, BGU13, BGU14, BGU22, BGU39 and BGU40 (Petrolina) (Fig. 2).

Genetic structure and gene flow: Two groups were identified based on ΔK (K = 2) (Fig. 3, Fig. 4). Of the clusters obtained from the similarity matrices, there is a group of three accessions (BGU58, BGU59 and BGU62) exclusive to the region of Januária, Minas Gerais that is located in Brazil Southeastern part of the semiarid region.

Fig. 3 Delta K, calculated with the second-order average rate of variation of the K probability divided by the standard deviation of the K probability. 

Fig. 4 Genetic structure of 24 umbu tree accessions from four semiarid regions, based on a Bayesian analysis, considering K = 2, obtained by the ΔK method, from 20 independent simulations for each number of possible clusters (k). 

The Bayesian analysis indicated the existence of only two groups, the three accessions from Januária, Minas Gerais in one group and the remaining accessions from Uauá, Bahia; Juazeiro, Bahia; and Petrolina, Pernambuco in another group (Fig. 4).

The analysis of molecular variance of 24 individuals the umbu tree from four distinct populations revealed that only 12 % of the genetic variability is between populations and 88 % of the variability is within populations (Table 2).

Table 2 Analysis of molecular variance (AMOVA) among four umbu tree populations assessed with ten microsatellite loci 

Source of variation d.f. SS MS p-value Total genetic variation (%) ΦST Nm
Among populations 3 27.833 9.278 < 0.001 12 ΦST = 0.12 1.865
Within population 20 102.833 5.142 < 0.001 88 1- ΦST = 0.88
Total 23 130.667

d.f.: Degree of freedom; SS: Sum of squares; MS: Mean squares; p-value based on 999 permutations; ΦST: Subpopulation genetic variance (S)/total genetic variance (T); Nm: Number of migrants calculated using the method by Wright (1949) [(1-ΦST)/(4 ΦST)].

Discussion

Microsatellite markers used to estimate the genetic diversity and structure of S. tuberosa accessions had good amplification patterns and an average of 3.5 alleles per locus. This is similar to the study of Balbino, Caetano, and Almeida (2018) who found an average of 2.7 alleles per locus for the same species Cristóbal-Pérez, Fuchs, Harvey, and Quesada (2019) evaluated the genetic variability of another species of the genus (S. purpurea) using 24 microsatellites and found an average of 5.88 alleles. Silva et al. (2017) also found high allelic diversity (6.97 alleles per locus) using ISSR markers to characterize the genetic diversity of S. mombin.

The polymorphic information content (PIC) values were moderate (0.5 > PIC > 0.25) to highly (PIC > 0.5) informative, according to the classification by Botstein, White, Skolnick, and Davis (1980), except for TUB93 that had a less-informative PIC value of 0.195, and the SPO4, SPO14 and TUB94 loci were the most informative (Table 1; Silva et al., 2017) estimated the genetic diversity of S. mombin using ISSR markers and found PIC values above 0.250 for most of the markers used, similar to our study.

For the SPO8, SPO14, TUB84, TUB93, TUB94 and TUB103 markers, the expected heterozygosity was higher than the observed heterozygosity, meaning high genetic variability and mixing of populations. Cristóbal-Pérez et al. (2019) also found higher expected heterozygosity values when evaluating 139 individuals of S. purpurea, from three Mexican localities, based on 10 polymorphic SSR loci.

The expected and observed heterozygosity values were similar to those observed by Balbino et al. (2018), who found values between 0.158 and 0.607, and 0.170 and 0.781, respectively, in a study about the phylogeographic pattern of S. tuberosa using accD-psaI plastid sequences and SSR markers of individuals from 20 localities of Brazil Northeastern.

In the cluster analysis, the cophenetic correlation coefficient of 1.0 indicates the confidence of the data and shows there is a good fit between the genetic distances, original matrix, and graphic representation. Dendrogram of the 24 individuals analyzed, based on the UPGMA method, showed four groups. Reduced cluster number was reported by Santana et al. (2011), when analyzed the genetic diversity among 17 Spondias sp. accessions from Brazil Northeast, with 25 ISSR markers, a king of universal SSR. Based on ISSR, Silva et al. (2017) reported nine groups with a sampling of 126 individuals of S. mombin from three populations, collected in the Mato Grosso State, Brazil.

The BGU56 accession comprised the first group diverging from the other studied accessions. The second group only contains the accessions from Januária-MG, municipality located in the Brazilian Southeast Region, geographically distant, around 2 000 Km, from the other analyzed regions. The fourth group has the most accessions and comprises individuals from the Juazeiro, state of Bahia, Uauá, state of Bahia, and Petrolina, state of Pernambuc regions, which are relatively close to each other, maximum of 150 km, in the same ecogeographic region. Thus, we can infer that their proximity justifies this cluster. Based on phenotypic characters, Santos (1997) concluded that variability in umbu tree is uniformly distributed in the semiarid region of Brazil. Differently, Santos et al. (2008), based on AFLP markers, concluded that the genetic variability of umbu tree is not uniformly distributed in this region and that geographic barriers or edaphoclimatic conditions have limited the crossing and frequency of the alleles among populations. The present study also indicates that the variability of umbu tree is not uniformly distributed in the semiarid region, since individuals from Januária form nearly exclusive group, while individuals from Uauá, Juazeiro and Petrolina are almost all in the same group.

Based on Bayesian statistics, two genetic groups (K) were found in this study. Of these two groups obtained for K, one includes the three accessions (BGU58, BGU59 and BGU62) from Januária-MG, in the southeast part of the Brazilian semiarid region, and the other contains the remaining accessions from Uauá, Juazeiro, and Petrolina. Balbino et al. (2018) also found K = 2 in a study of the phylogeographic pattern of S. tuberosa using sequences of the accD-psaI plastid region and six SSR markers for individuals from 20 localities of the Brazil North-eastern. By Bayesian analysis, the two groups found in the present study can be seen at a cutoff point of 0.40 in the UPGMA dendrogram (Fig. 1), in three groups. Costa and Santos (2017) also reported concordances between UPGMA and Bayesian analyses when studying accessions of Psidium (guava) with SNPs.

The analysis of molecular variance indicated moderate genetic differentiation diversity among population (12 % of the variability). A similar result was found in the study of Balbino et al. (2018), where the authors detected 13 % genetic variability among populations of S. tuberosa from Brazil North-eastern regions. These data go to what was reported by Paiva (1998) who noted that in natural plant populations in tropical regions most genetic variability is preserved within populations. Still about, according to Wright (1965), FST (=ΦST) values above 0.25 indicate high levels of genetic differentiation and an FST value of 0.12 indicates moderate differentiation.

Using the AFLP molecular marker, Santos et al. (2008) studied the distribution of the genetic variability of umbu tree in the semiarid region of Brazil and found high genetic differentiation (FST = 0.3138), suggesting that this species has restricted flow, with less than one migrant per generation (Nm = 0.567), and high variability between populations. Using an isoenzymatic polymorphism analysis Silva, Martins, and Oliveira (2009) estimated the genetic diversity and structure of S. lutea populations in the forest zone in Pernambuco State, in Northeastern Brazil, and found an Nm value of 5.27, which differs from that found in the present study.

The cluster analysis, AMOVA and Bayesian analysis of the present study indicate that the genetic diversity of S. tuberosa is not uniformly distributed in the Januária-MG, Juazeiro-BA, Uauá-BA and Petrolina-PE regions. Thus, germplasm from a greater number of populations should be collected in other Brazilian regions to increase the genetic diversity of the germplasm collection maintained at Embrapa Semiárido, Brazil.

The FST value of 0.12 indicates moderate genetic differentiation among the S. tuberosa populations from Januária-MG, Juazeiro-BA, Uauá-BA and Petrolina-PE, suggesting that the genetic variability of the accessions of the Embrapa germplasm collection is moderately structured in function of origin. The genetic diversity of S. tuberosa is not uniformly distributed in the four studied Brazilian semiarid regions and germplasm expedition should consider sampling in other regions to increase the collection variability.

Ethical statement: authors declare that they all agree with this publication and made significant contributions; that there is no conflict of interest of any kind; and that we followed all pertinent ethical and legal procedures and requirements. All financial sources are fully and clearly stated in the acknowledgements section. A signed document has been filed in the journal archives.

Acknowledgments

The authors would like to thank Embrapa (Brazilian Agricultural Research Corporation) for the infrastructure and financial support for carrying out the experiments; CAPES (Coordination for the Improvement of Higher Education Personnel) and CNPq (National Council for Scientific and Technological Development) for the grant of the postgraduate scholarship. We also thank Mr. Geraldo Freire dos Santos for his support in the field and Ms. Tatiana Ayako Taura for the support on the map.

References

Aguilar-Barajas, E., Sork, V.L., González-Zamora, A., Rocha-Ramírez, V., Arroyo-Rodríguez, V., & Oyama, K. (2014). Isolation and characterization of polymorphic microsatellite loci in Spondias radlkoferi (Anacardiaceae). Applications in Plant Sciences, 2(11), 1400079. [ Links ]

Balbino, E., Caetano, B., & Almeida, C. (2018). Phylogeographic structure of Spondias tuberosa Arruda Câmara (Anacardiaceae): seasonally dry tropical forest as a large and continuous refuge. Tree Genet Genomes, 14(5), 1-8. [ Links ]

Balbino, E., Martins, G., Morais, S., & Almeida, C. (2019). Genome survey and development of 18 microsatellite markers to assess genetic diversity in Spondias tuberosa Arruda Câmara (Anacardiaceae) and cross-amplification in congeneric species. Molecular Biology Reports, 46(3), 3511-3517. [ Links ]

Botstein, D., White, R.L., Skolnick, M., & Davis, R.W. (1980). Construction of a genetic linkage map in man using restriction fragment length polymorphisms. American Journal of Human Genetics, 32(3), 314-331. [ Links ]

Costa, S.R., & Santos, C.A.F. (2013). Allelic database and divergence among Psidium accessions by using microsatellite markers. Genetics and Molecular Research, 12(4), 6802-6812. [ Links ]

Costa, S.R., & Santos, C.A.F. (2017). Genetic divergence among Psidium accessions based on single nucleotide polymorphisms developed for Eucalyptus. Genetics and Molecular Research, 16(2), gmr16029566. [ Links ]

Creste, S., Neto, A.T., & Figueira, A. (2001). Detection of single sequence repeat polymorphisms in denaturing polyacrylamide sequencing gels by silver staining. Plant Molecular Biology Reporter, 19(4), 299-306. [ Links ]

Cristóbal-Pérez, E.J., Fuchs, E.J., Harvey, N., & Quesada, M. (2019). Isolation and characterization of microsatellites loci in Spondias purpurea (Anacardiaceae) and cross amplification in congeneric species. Molecular Biology Reports, 46(5), 5581-5585. [ Links ]

Doyle, J.J., & Doyle, J.L. (1990). Isolation of plant DNA from fresh tissue. Focus, 12(13), 39-40. [ Links ]

Earl, D.A., & vonHoldt, B.M. (2012). STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conservation Genetics Resources, 4(2), 359-361. [ Links ]

Evanno, G., Regnaut, S., & Goudet, J. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Molecular Ecology, 14(8), 2611-2620. [ Links ]

Excoffier, L., Smouse, P.E., & Quattro, J.M. (1992). Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics, 131(2), 479-491. [ Links ]

Meirmans, P.G., & Hedrick, P.W. (2011). Assessing population structure: FST and related measures. Molecular Ecology Resources, 11(1), 5-18. [ Links ]

Paiva, J.R. (1998). Melhoramento genético de espécies agroindustriais na Amazônia: estratégias e novas abordagens. Fortaleza, CE: Embrapa-CNPAT. [ Links ]

Peakall, R.O.D., & Smouse, P.E. (2006). GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes, 6(1), 288-295. [ Links ]

Pires, M.D.G.D.M. (1990). Estudo taxonômico e área de ocorrência de Spondias tuberosa Arr. Cam. (umbuzeiro) no Estado de Pernambuco-Brasil (Dissertação de Mestrado). Universidade Federal Rural de Pernambuco, Recife, PE, Brasil. [ Links ]

Pritchard, J.K., Stephens, M., & Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics, 155(2), 945-959. [ Links ]

Ramos, S., Queiroz, M.D., Romão, R.L., & Silva Júnior, J.D. (2008). Germoplasma vegetal conservado no nordeste brasileiro: situação atual, prioridades e perspectivas. Magistra, Cruz das Almas, 20(3), 205-319. [ Links ]

Rohlf, F.J. (2000). NTSYS-pc numerical taxonomy and multivariate analysis system (Version 2.10). Setauket: Exeter Software. [ Links ]

Santana, I.B.B., Oliveira, E.J.D., Soares Filho, W.D.S., Ritzinger, R., Amorim, E.P., Costa, M.A.P.D.C., & Moreira, R.F.C. (2011). Variabilidade genética entre acessos de umbu-cajazeira mediante análise de marcadores ISSR.Revista Brasileira de Fruticultura, 3(3), 868-876. [ Links ]

Santos, C.A.F. (1997). Dispersão da variabilidade fenotípica do umbuzeiro no semi-árido brasileiro. Pesquisa Agropecuária Brasileira, 32(9), 923-930. [ Links ]

Santos, C.A.F., Rodrigues, M.A., & Zucchi, M.I. (2008). Variabilidade genética do umbuzeiro no Semi-Árido brasileiro, por meio de marcadores AFLP.Pesquisa Agropecuária Brasileira, 43(8), 1037-1043. [ Links ]

Silva, B.M., Rossi, A.A.B., Tiago, A.V., Schmitt, K.F.M., Dardengo, J.F.E., & Souza, S.A.M. (2017). Genetic diversity of Cajazeira (Spondias mombin L.) in three geographic regions. Genetics and Molecular Research, 16(1), gmr16018946. [ Links ]

Silva, E.F.D., Martins, L.S.S., & Oliveira, V.R.D. (2009). Diversity and genetic structure in caja tree (Spondias mombin L.) populations in Northeastern Brazil. Revista Brasileira de Fruticultura, 31(1), 171-181. [ Links ]

Silva, S.M., & Alves, R.E. (2008). Desenvolvimento e fisiologia da maturação de frutos do gênero Spondias. In I.E. Lederman, J.S.D. Lira Júnior, & J.F.D. Silva Júnior (Eds.), Spondias no Brasil: umbú, cajá e espécies afins (pp. 23-33). Recife, Brasil: Empresa Pernambucana de Pesquisa Agropecuária, IPA/UFRPE. [ Links ]

Turchetto-Zolet, A.C., Turchetto, C., Zanella, C.M., & Passaia, G. (2017). Marcadores Moleculares na Era Genômica: Metodologias e Aplicações. Ribeirão Preto, São Paulo: Sociedade Brasileira de Genética. [ Links ]

Wright, S. (1949). Genetical structure of populations. Annals of Eugenics, 15(1), 323-354. [ Links ]

Wright, S. (1965). The interpretation of population structure by F-statistics with special regard to systems of mating. Evolution, 19(3), 395-420. [ Links ]

Received: October 13, 2020; Revised: March 25, 2021; Accepted: May 10, 2021

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License