Research Article

Genetic variability analysis of Byrsonima crassifolia germplasm collected in Pará State using ISSR markers

Published: October 17, 2016
Genet. Mol. Res. 15(4): gmr15048887 DOI: 10.4238/gmr15048887

Abstract

Native of the Amazon, the nanche (Byrsonima crassifolia) is a fruit cultivated by family farmers and used in cooking; as such, it represents an opportunity for regional agribusiness. The Embrapa Eastern Amazon set up an active germplasm bank (BAG) consisting of 22 accessions sampled in 11 municipalities of Pará State. Due to its economic potential, there is an interest to advance the genetic breeding program of this species. The aim of this study was to characterize the BAG nanche collection using inter-simple sequence repeat (ISSR) markers. Accessions were genotyped using 23 pre-selected ISSR primers resulting in 109 amplified polymorphic and 51 monomorphic bands. With eight polymorphic bands each, the most polymorphic primers were UBC 809 and UBC 848. An unweighted pair-group method with arithmetic average cluster analysis based on Jaccard’s coefficient indicated that the individuals clustered into two distinct groups. Accessions Igarapé Açu-2 and Augusto Corrêa-Pl 1 were most similar. The genetic dissimilarity values ranged from 0.10 to 0.59. We conclude that the ISSR markers were efficient in detecting polymorphisms in the nanche accessions, and that it is possible to infer the genetic variability among accessions of the collection. This demonstrate the importance of using molecular markers in poorly studied species and the advantages that this information can bring to the genetic improvement of such species.

INTRODUCTION

Nanche (Byrsonima crassifolia (L.) H.B.K.) is considered a semi-domesticated shrub species native of the Amazon and widely distributed in tropical America (Cavalcante, 2010). It is the most important fruit crop of the genus Byrsonima and the most cultivated, since it produces fleshy fruits used in cooking, as well as for making sweets and liqueurs (Cavalcante, 2010; Carvalho and Nascimento, 2013). This species is rich in polyphenols and flavonoid compounds that provide antioxidants. Nanche is therefore considered a functional food (Silva et al., 2007; Souza et al., 2008). It is considered an allogamous plant with auto-compatibility (Sihag, 1995), although low formation of fruits has been observed in self-pollination experiments (Pereira and Freitas, 2002). Cultivation of this species presents an opportunity for diversification in local production structures, expanding the availability of agribusiness products in the Amazon region. It is typically grown by family farmers, has a guaranteed consumer market, and can be produced on a commercial scale due to the development of propagation methods that accelerate and standardize the seed germination (Carvalho and Nascimento, 2008).

The Embrapa Eastern Amazon has assembled an active germplasm bank (BAG) of nanche germplasms and identified promising materials for breeding purposes (Carvalho and Nascimento, 2013). Further experimental evaluations of this germplasm may allow for recommendations of appropriate and stable cultivars for producing regions. Thus, due to this species’ potential, it is necessary to realize a large-scale study of the BAG, which conserves accessions from places of origin of this species. Knowledge about the available germplasm bank can enable traceability and the rescue of accessions by identifying potential replacements in case of loss of plants. Furthermore, the collection may be enriched and kept in order, to preserve its genetic variability (Borém and Miranda, 2013). In addition, the knowledge of the genetic variability of the conserved accessions is indispensable for the development of genetic breeding programs.

Molecular markers have been used in genetic diversity studies in plant species, to determine the variation in collections and germplasm banks, to direct crosses, and to identify markers that describe plant materials and recommended cultivars (Ferreira and Grattapaglia, 1998). Despite the economic importance of B. crassifolia, there are few studies using molecular markers in this species. Inter-simple sequence repeats (ISSRs) are dominant and represent genetic variation within microsatellite regions in the genome. They have been used to estimate genetic diversity and genetic parameters of populations of multiple plant species (Rivas et al., 2013; Lorenzoni et al., 2014; Rossi et al., 2014; Silva et al., 2014; Ghariani et al., 2015; Arrigoni-Blank et al., 2016; Shen et al., 2016). The advantage of ISSRs is the possibility to obtain a high number of polymorphic loci without the need of genome sequence knowledge (Faleiro, 2007).

The aim of this study was to verify the genetic variability of B. crassifolia accessions from the Embrapa Eastern Amazon germplasm bank, using ISSR markers focusing on the identification of superior materials for conventional genetic improvement.

MATERIAL AND METHODS

Plant material

A total of 22 accessions of B. crassifolia were used in this study (Table 1). The accessions were collected in 11 cities of the Pará State (Figure 1), and are conserved as clones at the Embrapa Eastern Amazon germplasm bank (Table 1).

Origin of nanche (Byrsonima crassifolia) accessions used in this study collected from the Pará State and belonging to the Eastern Amazon germplasm bank.

Accession Collection site
1 - Açu (AÇU) Tomé Açu
2 - Santarém-1 (SANT-1) Santarém
3 - Santarém-2 (SANT-2) Santarém
4 - Igarapé Açu-1 (IG.AÇU-1) Igarapé Açu
5 - Tocantins-2 (TO-2) Nova Ipixuna
6 - Cristo Belém
7 - Maracanã-2 (MARAC-2) Maracanã
8 - Guataçara 1.1 (GUAT 1-1) Terra Alta
9 - Tocantins-1 (TO-1) Nova Ipixuna
10 - São José (S. JOSÉ) Igarapé Açu
11 - Aurora 1 (AUR 1) Aurora do Pará
12 - FCAP-7 Belém
13 - Baião-2 Baião
14 - Augusto Corrêa-Pl 1 (AC 1) Augusto Corrêa
15 - Augusto Corrêa-Pl 2 (AC 2) Augusto Corrêa
16 - Kishi-3 Monte alegre
17 - Igarapé Açu-2 (IG.AÇU-2) Igarapé Açu
18 - Igarapé Açu-3 (IG.AÇU-3) Igarapé Açu
19 - Igarapé Açu-4 (IG.AÇU-4) Igarapé Açu
20 - Unnamed (UN) Unknown origin
21 - Aurora-2 (AUR 2) Aurora do Pará
22 - Xininga-5 (XIN-5) Baião

Map showing the collection sites of 22 accessions of nanche (Byrsonima crassifolia) in 11 municipalities in Pará State that were characterized with ISSR markers.

DNA analyses

Leaves were sampled from the accessions and total genomic DNA was extracted according to a procedure similar to that of Doyle and Doyle (1990). The tissue was macerated with liquid nitrogen and polyvinylpyrrolidone (PVP). Subsequently, 3 mL cetyltrimethylammonium bromide (CTAB) extraction buffer (2% CTAB, 5 M NaCl, 0.5 M EDTA, PVP, 1 M Tris-HCl, and sterile water) were added to the macerate. The macerate was homogenized and incubated in hot water bath at 65°C for 1 h. Following this, chloroform:isoamyl alcohol (24:1) was added followed by homogenization, and the samples were centrifuged for 10 min at 11,952 g. To precipitate the DNA, 3 mL 95% ethyl alcohol was added to the supernatant and the samples were centrifuged for 10 min at 11,952 g. Next, the precipitate was washed with 70% ethyl alcohol for 10 min and centrifuged at 2991 g. DNA samples were resuspended in 300 μL TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8.0) and RNAse. DNA was quantified on a 1% agarose gel, using lambda phage DNA as a standard, at different concentrations (50, 100, and 200 ng/μL).

The samples were genotyped using 23 ISSR primers (University of British Columbia, Vancouver, Canada) that had been identified after the primer selection phase, using DNA from five B. crassifolia accessions from the Embrapa Eastern Amazon BAG. A polymerase chain reaction was performed using a final volume of 20 μL, containing 10 ng genomic DNA, 75 µM each dNTP, 2.0 µM primer, 1.0 mg/mL bovine serum albumin, and reaction buffer containing 1.2 mM MgCl2 and 0.2 U Taq DNA polymerase (Applied Biosystems, Carlsbad, CA, USA). For this step, we tested 100 ISSR primers (UBC 801-900), and the reactions were carried out in 0.2-mL microtubes and amplified on an Amplitherm TX96 thermocycler programmed for 35 cycles. First, there was a denaturation phase at 95°C for 5 min. Then, each cycle consisted of DNA denaturation at 95°C for 1 min, primer annealing at temperatures ranging from 50 to 62°C (depending on the primer, Table 2) for 45 s, and elongation at 72°C for 2 min. Following the 35 cycles, there was the final extension at 72°C for 5 min. The reaction products were run on a 1.5% agarose gel (Thermo Fisher Scientific, Boston, MA, USA) prepared with 1.0X TBE buffer (0.45 M Tris-borate and 0.01 M EDTA). Gels were run in a horizontal electrophoresis unit containing 1.0X TBE at a constant voltage of 80 V, were visualized on an ultraviolet light transilluminator (UVB transilluminator, Loccus, LTB-20x20 HE), and images were digitally captured. The fragments were compared with a 1-Kb DNA ladder marker (Applied Biosystems). The annealing temperatures for each primer were chosen according to the result of the reaction that offered the best band sharpness and intensity on the agarose gel for each primer evaluated, and the ability of discrimination of polymorphic bands. The primers were chosen based on their polymorphism.

Pre-selection of ISSR primers used in the molecular characterization of nanche (Byrsonima crassifolia) accessions collected in Pará State.

Primer No. of polymorphic bands Primer No. of polymorphic bands
808 5 844 3
809 4 845 5
811 5 846 6
812 5 848 6
813 4 850 4
814 5 855 4
815 2 856 4
825 5 857 6
826 4 858 3
835 6 888 5
836 4 889 4
843 4

Total polymorphic bands = 103. The number of polymorphic bands shown had an annealing temperature of 47°C.

Data analyses

The data obtained from the analysis of amplified bands using the 23 primers according to the quality of the polymorphic bands visualized on the gel. Bands with the same run pattern were considered to be the same locus. Presence of a band was scored as (1) and band absence as (0), generating a binary matrix. Only polymorphic bands were analyzed.

The matrix of genetic similarity was generated using PAST (Hammer et al., 2001), based on Jaccard’s coefficient:

sgij=aa+b+c

(Equation 1)

where a is the number of events where the band occurs in both genotypes; b is the number of events where the band occurs only in genotype i; and c is the number of events where the band occurs only in genotype j.

Based on the genetic similarity matrix, a dendrogram was generated using the unweighted pair group method with arithmetic mean (UPGMA). The relationship between the similarity matrix and the dendrogram was estimated using the coefficient of cophenetic correlation (CCC), according to Sokal and Rohlf (1962).

The assignment of genotypes to groups and relatedness among groups were assessed using STRUCTURE v. 2.3.3 (Pritchard et al., 2000). Admixture and correlated allele models without prior population information were used. Following a burn-in period of 100,000 iterations, ten independent runs were carried out for each number of clusters (K, from 1 to 10), each with 100,000 iterations. The choice of the most likely K was carried out by calculating the ∆K statistic. This statistic is based on the rate of change in the log probability of the data among successive K values. The selection was done as described by Evanno et al. (2005), using STRUCTURE HARVESTER (Earl and VonHoldt, 2012). Among the 10 runs per value of K, the run with the highest maximum likelihood was used to assign individual genotypes to clusters. STRUCTURE uses a Bayesian approach and indicates the number of genetic groups (K value) that best fits the data

RESULTS

The pre-selection result of primers indicated polymorphisms in 23 ISSR using an annealing temperature of 47°C. These primers showed an average polymorphism of 4.4 bands per primer. We noted that some primers resulted in fewer bands using the annealing temperature of 47°C, compared to DNA amplification using higher annealing temperatures. This was observed using primers UBC 809, UBC 835, UBC 845, UBC 846, UBC 848, and UBC 850. No polymorphism was observed in 77 primers for the five clones evaluated at a temperature of 47°C. The 23 polymorphic primers were used for the characterization of B. crassifolia accessions from the Embrapa Eastern Amazon.

The 23 selected primers had an average annealing temperature of 52.43°C. The highest annealing temperature used was 60°C, which was used for UBC 888 and UBC 889, whereas the lowest annealing temperature (47°C) was used for UBC 848 and UBC 850 (Table 3). The annealing temperature of 53°C showed the highest amplification quality for the largest number of primers (UBC 811, UBC 813, UBC 843, UBC 855, UBC 857, and UBC 858), followed by 52°C annealing temperature, which was optimal for five primers (UBC 825, UBC 826, UBC 835, UBC 846, and UBC 856). Only the ISSR-845 primer obtained satisfactory amplification products at the 48°C annealing temperature.

Ideal annealing temperature, number of polymorphic bands, and percent polymorphism for the selected ISSR primers used in the molecular characterization of nanche (Byrsonima crassifolia) accessions collected in Pará State.

Primer 5' to 3' sequence Annealing temperature (°C) No. of bands No. of polymorphic bands Polymorphism (%)
808 (AG)8 C 50 7 5 71.43
809 (AG)8 G 57 11 8 72.73
811 (GA)8 C 53 9 5 55.56
812 (GA)8 A 49 7 5 71.43
813 (CT)8 T 53 6 4 66.67
814 (CT)8 A 49 7 5 71.43
815 (CT)8 G 49 3 2 66.67
825 (AC)8 T 52 7 5 71.43
826 (AC)8 C 52 6 4 66.67
835 (AG)8 YC 52 4 2 50.00
836 (AG)8 YA 56 7 5 71.43
843 (CT)8 RA 53 5 4 80.00
844 (CT)8 RC 56 5 3 60.00
845 (CT)8 RG 48 10 7 70.00
846 (CA)8 RT 52 7 5 71.43
848 (CA)8 RG 47 12 8 66.67
850 (GT)8 YG 47 9 6 66.67
855 (GT)8 YT 53 6 4 66.67
856 (AC)8 YA 52 5 4 80.00
857 (AC)8 YG 53 8 6 75.00
858 (TG)8 RT 53 4 3 75.00
888 BDB (CA)7 60 8 5 62.50
889 DBD (AC)7 60 7 4 57.14
Total Mean = 52.43 160 109 Mean = 68.125

B: A, G, T; Y: C, T; R: A, G; D: A, G, T.

After the first evaluation, the 23 primers were selected and revealed variation in the number of amplified fragments, totaling 160 DNA fragments. We were able to identify 109 polymorphic bands (Table 3), with an average of 4.74 markers per ISSR primer and 68.13% polymorphism. The number of polymorphic bands per primer varied from two (UBC 815) to eight (UBC 809 and UBC 848), and the polymorphism ranged from 50 (UBC 835) to 80% (UBC 843 and UBC 856).

The values of genetic dissimilarities based on Jaccard’s coefficient, varied from 0.10 to 0.59, with an average of 0.33 (Figure 2A). The mean among individual genetic similarity was 30.0% with an average standard deviation of 0.08. Most combinations were found in the range of 0.199-0.296, whereas only three combinations were found in the range of 0.494-0.59 (Figure 2B).

Analysis of genetic identity between accessions of Byrsonima crassifolia, obtained by ISSR markers. A. Matrix of genetic identity. Maximum and minimum values are shown in bold face. AÇ: AÇU; S1: SANT-1; S2: SANT-2; I1: IG.AÇU-1; T2: TO-2; CR: CRISTO; MR: MARAC-2; GT: GUAT 1.1; T1: TO-1; SJ: S. JOSÉ; A1: AUR 1; F7: FCAP-7; B2: BAIÃO-2; C1: AC 1; C2: AC 2; KS: KISHI-3; I2: IG.AÇU-2; I3: IG.AÇU-3; I4: IG.AÇU-4; UN: UN; A2: AUR 2; XN: XIN-5. B. Distribution of genetic identity among accessions.

Based on the analysis of genetic similarity, the individuals could be allocated to groups (Figure 3). The highest genetic similarity was detected among an accession from Igarapé Açu (Igarapé Açu-2) and one from Augusto Corrêa (Augusto Corrêa-2), whereas the most divergent accessions were Nova Ipixuna (Tocantins-2) and Igarapé Açu (São José) (Figure 3). The dendrogram based on the UPGMA method and Jaccard’s genetic similarities among the 22 accessions indicated the formation of one group of individuals collected in the northeast of Pará, and another group of individuals containing accessions mainly from southern and western Pará. However, Baião-2 and Xininga-5, collected in the northeast region of the state, were grouped in the latter group, and the accession Tocantins-2 did not cluster with any other accession (Figure 3).

Dendrogram representing the genetic dissimilarity between 22 Byrsonima crassifolia accessions collected in Pará State obtained using the UPGMA method and the arithmetic complement of Jaccard’s index as a measure of dissimilarity. The analysis was based on ISSR molecular characterization. Bootstrap values are shown.

The analysis based on the dendrogram confirmed a clustering concordant with the accessions’ geographic origin. However, the sample FCAP-7, from Belém, was closer to Baião-2 (Baião locality). The sample of unknown origin (UN) clustered with a sample from Igarapé Açu, while the accession collected in Tomé Açu (AÇU) was genetically distant from all other individuals of the collection. The result of the cophenetic correlation analysis showed a 72% association between the distances obtained using the Jaccard coefficient (the dissimilarity matrix) and the dendrogram (the cophenetic matrix). Cruz and Carneiro (2003) stated that the higher the CCC value, the smaller the distortion caused by grouping individuals using the UPGMA. The CCC value was considered appropriate, since values of r ≥ 0.56 are considered ideal, which reflects a good agreement with the values’ genetic similarity (Patto et al., 2004).

The STRUCTURE analysis resulted in the highest ΔK value for K = 2, in addition to a few other suboptimal K values (Figure 4A). This suggests that there are two genetically homogeneous accession groups. The structure bar plot shown in Figure 4B also shows two distinct groups.

ΔK analysis using STRUCTURE (Pritchard et al., 2000) based on ISSR data of 22 Byrsonima crassifolia accessions. A. Indication of the optimal group number present in the sample (K = 2) as determined by the method of Evanno et al. (2005). B. Structure bar plot illustrating the two genetic groups identified by the analysis. Accessions numbers are represented by vertical colored lines. Dark gray columns correspond to group 1 and light gray columns represents group 2. 1: AÇU; 2: SANT-1; 3: SANT-2; 4: IG.AÇU-1; 5: TO-2; 6: CRISTO; 7: MARAC-2; 8: GUAT 1.1; 9: TO-1; 10: S. JOSÉ; 11: AUR 1; 12: FCAP-7; 13: BAIÃO-2; 14: AC 1; 15: AC 2; 16: KISHI-3; 17: IG.AÇU-2; 18: IG.AÇU-3; 19: IG.AÇU-4; 20: UN; 21: AUR 2; 22: XIN-5.

DISCUSSION

B. crassifolia is considered a semi-domesticated plant and need of additional attention. Advancing the genetic breeding programs of the species may be achieved using molecular markers. Our evaluation of nanche accessions of the Embrapa Eastern Amazon BAG using ISSR markers revealed the need to identify specific temperatures for each of the primers used, as has previously been observed in other species (Pessoni, 2007; Souza et al., 2013; Pepineli et al., 2014). In a similar study conducted for Platonia insignis, 72 genotypes from the Embrapa Mid-North germplasm collection were evaluated using 18 primers. An average of 13 loci per primer was identified (Souza et al., 2013). Among the primers used, seven primers (UBC 811, UBC 825, UBC 826, UBC 855, UBC 856, UBC 888, and UBC 890) were also used in the current study. However, different annealing temperatures were used for P. insignis (Souza et al., 2013).

In the present study, 160 bands were amplified, resulting in an average of 7.0 loci per primer. The cashew germplasm collection of Embrapa Tropical Agroindustry, consisting of 91 samples, was characterized using 11 polymorphic primers (Pessoni, 2007). UBC 812 was analyzed on a polyacrylamide gel and 92.30% polymorphic bands were identified (Pessoni, 2007). In the present study, UBC 812 resulted in 71.43% polymorphic bands. There are several studies in the literature using dominant markers (Carvalho et al., 2012; Lorenzoni et al., 2014), including ISSRs, that consider the percentage polymorphic loci a genetic diversity measure.

Martínez et al. (2013) analyzed B. crassifolia trees found in two different locations in the state of Tabasco, Mexico, using six random amplification of polymorphic DNA (RAPD) primers, three of which were polymorphic. The within-population genetic diversity was higher (85.20%) than the between-population genetic diversity (14.80%), and it was concluded that dominant markers may be considered useful for identifying genetic variation of this species (Martínez et al., 2013). The low polymorphism found was taken to indicate that the genetic variability was small for the measured material, and that the RAPD primers used in the study showed low efficiency to recognize homologous DNA regions (Martínez et al., 2013).

Lorenzoni et al. (2014) studied the genetic diversity among accessions of Rollinia mucosa using ISSR markers. They identified dissimilarity values ranging from 0.0909 to 0.5147, showing high genetic diversity among the accessions. Thus, it can be inferred that the results of genetic dissimilarity in our study reflect a considerable genetic variability among the nanche accessions. This is justified by the fact that the evaluated nanche collection consisted of a small number of individuals.

Nanche is considered an allogamous and auto-compatible species (Sihag, 1995), although Pereira and Freitas (2002) observed low formation of fruits in flowers subjected to self-pollination. Our results reflect the genetic variability of the accessions studied, which highlights the importance of collecting more materials in the area represented in this study. Such an effort may identify genotypes with good prospect of commercial exploitation that can be used in genetic breeding programs.

Our results demonstrate the importance of using molecular markers for discrimination and identification of genetic materials. This makes it possible to quantify the genetic variability among individuals of the same species, even in species that are not considered domesticated and that have been little studied, as in the case of nanche. It also makes it possible to group genetically similar individuals. ISSR markers were efficient in detecting molecular polymorphisms in nanche. However, in future studies it will be necessary to expand the collection area, increase the number of prospected genotypes, and to include natural populations. It will then be possible to obtain information regarding the evolutionary and population mechanisms that may be used to decide on maintenance strategies and species conservation, particularly in the choice of individuals and areas to be sampled. Understanding the population dynamics of this understudied species could facilitate the choice of contrasting genotypes to be used in genetic breeding and in maintenance strategies. Such genetic information can be used together with morpho-agronomical information, in order to direct and accelerate the improvement of the species, in addition to making decisions about collecting new plant materials, conservation of the germplasm, and management of B. crassifolia genetic resources.