Impact Factor 4.298

The #1 most cited and #1 largest open-access plant science journal

SUPPLEMENTAL DATA

Original Research ARTICLE

Front. Plant Sci., 24 February 2017 | https://doi.org/10.3389/fpls.2017.00243

Genetic Structure and Selection of a Core Collection for Long Term Conservation of Avocado in Mexico

Luis F. Guzmán1, Ryoko Machida-Hirano2*, Ernesto Borrayo2, Moisés Cortés-Cruz1, María del Carmen Espíndola-Barquera3 and Elena Heredia García4
  • 1Centro Nacional de Recursos Genéticos, Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias, Tepatitlán de Morelos, Mexico
  • 2Gene Research Center, University of Tsukuba, Tsukuba, Japan
  • 3Fundación Salvador Sánchez Colín, CICTAMEX S. C., Ignacio Zaragoza, Col. Centro, Coatepec Harinas, Mexico
  • 4Campo Experimental Bajío, Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias, Celaya, Mexico

Mexico, as the center of origin of avocado (Persea americama Mill.), harbors a wide genetic diversity of this species, whose identification may provide the grounds to not only understand its unique population structure and domestication history, but also inform the efforts aimed at its conservation. Although molecular characterization of cultivated avocado germplasm has been studied by several research groups, this had not been the case in Mexico. In order to elucidate the genetic structure of avocado in Mexico and the sustainable use of its genetic resources, 318 avocado accessions conserved in the germplasm collection in the National Avocado Genebank were analyzed using 28 markers [9 expressed sequence tag-Simple Sequence Repeats (SSRs) and 19 genomic SSRs]. Deviation from Hardy Weinberg Equilibrium and high inter-locus linkage disequilibrium were observed especially in drymifolia, and guatemalensis. Total averages of the observed and expected heterozygosity were 0.59 and 0.75, respectively. Although clear genetic differentiation was not observed among 3 botanical races: americana, drymifolia, and guatemalensis, the analyzed Mexican population can be classified into two groups that correspond to two different ecological regions. We developed a core-collection by K-means clustering method. The selected 36 individuals as core-collection successfully represented more than 80% of total alleles and showed heterozygosity values equal to or higher than those of the original collection, despite its constituting slightly more than 10% of the latter. Accessions selected as members of the core collection have now become candidates to be introduced in cryopreservation implying a minimum loss of genetic diversity and a back-up for existing field collections of such important genetic resources.

Introduction

Mexican avocado (Persea americana Mill.) has been insufficiently studied in its own center of origin despite the wide genetic diversity the region hosts. Understanding its unique population structure and domestication history, as well as exploring strategies for its conservation, has been a long-neglected priority.

Mexico is the world's largest producer of avocados, representing over one-third of global production and harvested area. Avocado is one of the most economically important fruit trees to Mexico and Central America. It is an evergreen subtropical species that has adapted to different climate ranges, extending its production to North and South America, the Caribbean, Asia, Africa, Middle East, and Europe (Food Agriculture Organization of the United Nations, 2014 WEB reference). Such geographical outspread implies adaptation to each specific local environment, which could derive into a reduction of genetic variation due to isolation of populations.

Avocado is a member of a large family called Lauraceae, constituted by 50 genera, including the genus Persea. There are three botanical races: (1) P. americana var. drymifolia (Schlecht. & Cham.) Blake (Mexican race), (2) P. americana var. guatemalensis L. Wms. (Guatemalan race), and (3) P. americana var. americana Mill. (West Indian race) (Table 1). They are recognized as commonly cultivated avocadoes (Ashworth et al., 2011) and their hybrids are currently cultivated worldwide (Borrone et al., 2009). These botanical races are determined by morphological, geographical, physiological, biochemical, molecular, and commercial aspects although some experts have proposed that differences among them are unclear (Ashworth et al., 2011). In addition, within P. americana, other taxa—var. nubigena (L. O. Williams) Kopp (1966) and var. costaricencis Ben-Ya'acob et al. (2003)—are also recognized.

TABLE 1
www.frontiersin.org

Table 1. Summary of main characteristics of botanical races of cultivated avocado.

According to archeological evidence, avocado consumption began more than 7,000 years ago (Smith, 1966) and it has been considered as one of the trees to be domesticated previous to annual crops in Mesoamerica (Galindo-Tovar et al., 2008). Furthermore, ancient ethnic groups may have contributed to some extent to the divergence of avocado into different botanical races through selection of the variants that best adapted to different environmental conditions (Cañas-Gutiérrez et al., 2015). Founded on genetic markers, there is a wide variety of studies regarding the genetic diversity and domestication history of avocado (AFLP's, Cañas-Gutiérrez et al., 2015; RAPD's, Sharon et al., 1997; RFLP's and VNTR's, Chanderbali et al., 2008; SNP's, Chen et al., 2008, 2009). Among these, microsatellites are commonly used for genetic diversity studies due to their reproducibility, high variability, codominance inheritance, abundance, and distribution across the avocado genome (Ashworth et al., 2004; Gross-German and Viruel, 2013). Although genetic diversity in wild avocado varieties has been reported (Chen et al., 2008), this is not the case for an extensive genetic diversity evaluation and phylogenetic relationship that may focus specifically in the avocado germplasm present in Mexico.

Genetic diversity plays an important role not only in preservation of biodiversity, but also of ecological and cultural aspects, such as diet and economy (Sarr et al., 2008; Rincón-Hernández et al., 2011). Projects have been launched with the aim to avoid the loss of such important resources. The conservation of plant species depends greatly on whether their seeds are of orthodox or recalcitrant nature. The latter species are often conserved in field collections, as is the case for the local collection in Celaya Experimental Station of the National Forestry, Crops, and Livestock Research Institute (CEBAJ-INIFAP), which dates back to 1972 and whose main objective is to collect and conserve varieties of avocado local to Mexico. At present time, this collection has been used as a source for cultivation and breeding materials retrieval. Recently, an initiative has arisen for the establishment of the “National Avocado Germplasm Depository (BNGA),” under a joint effort between the Salvador Sanchez Colin Foundation, CICTAMEX and CEBAJ-INIFAP. The collection—consisting of more than 350 accessions, mainly of Mexican genotypes (P. americana var. drymifolia) collected from different localities across Mexico—aims to conserve, as backup, these accessions in different institutions.

For the long-term conservation and/or safety backup of field collection of recalcitrant species, as avocado is, two main approaches are usually taken: in vitro and cryopreservation (Efendi and Litz, 2003). Specific protocols for target species must be available for these approaches to be successful. In addition, preparation and maintenance of in vitro cultures require skilled technicians and a maintenance budget. Cryopreservation has been recognized as the safest alternative for long-term conservation of plant genetic resources, as it does not require continuous manipulations (González-Arnao et al., 2014). To ensure the conservation and sustainable utilization of recalcitrant species under limited resources, we propose: a combination of cryopreservation with Core-Collection (CC) selection, i.e., a reduction to a small-number subset of the larger germplasm collection representing the maximum possible genetic diversity with minimum repetitiveness (Frankel, 1984). The CC represents an ideal set of accessions to use for the establishment of cryopreservation protocols and for the initial introduction of accessions in cryostorage. The combination of CC and cryopreservation approaches will ensure a feasible long-term conservation of avocado by focusing on a limited number of representative elements with the safest conservation method.

In this study, we carried out a genetic diversity evaluation of avocado germplasm conserved in the CEBAJ-INIFAP genebank by means of microsatellite markers, in order to infer variation and phylogenetic relationships which will in turn provide the relevant information to deduce the domestication history of avocado and to identify a CC that contained a representative subset -in terms of genetic diversity- of the original collection to ensure efficient long-term germplasm preservation.

Materials and Methods

Plant Material and Genomic DNA Extraction

A total of 319 accessions were sampled at Campo Experimental Bajío, INIFAP, in Mexico. The collection includes 313 of Persea americana. Among these, information on botanical races [var. drymifolia (104), var. guatemalensis (15), var. americana (12), var. costareicensis (3), and hybrids (31)] was available for 165 accessions, while information about geographical origin was available for 299 accessions. Three other Persea species [P. longipes (2), P. nubigena (2), and P. schiedeana (1) and an accession from other genus (Beilmiedia anay)] are also included in the collection. Fresh leaves collected from individual trees in the field were transferred to the laboratory and freeze-dried on the same day. Total DNA was extracted by a modified tropical-plant-specific protocol (Huang et al., 2013) to overcome polysaccharide and phenolic compound complications. The extraction buffer contained 2 M of NaCl, 25 mM of ethylenediaminetetraacetic acid, 200 mM of Tris, 2% of cetyltrimethylammonium bromide, 2% of polyvinylpolypyrrolidone, 1% of lauroyl sarcosine, 20 mM of borax, and 140 mM of β-mercaptorthanol. During incubation (65°C for 45 min), 700 μl of dichloromethane was used instead of the original-protocol chloroform: phenol: isoamyl alcohol (25:24:1). Incubation and precipitation were performed twice, followed by the previously reported unmodified protocol.

Microsatellite Analysis

A total of 47 EST-SSR and the genomic SSR (Gross-German and Viruel, 2013) were tested using an aleatory generated number (4) of genotypes randomly selected from the BNGA collection. Within all the 47 sets of primers, we discarded primer pairs which produced an ambiguous allelic pattern, and either monomorphic or more than two alleles. Finally, 28 primer pairs (9 EST-SSR and 19 genomic SSR) were used for genotyping. A combination of Universal tailed primers with multiplex PCR approach was applied for genotyping. A specific DNA sequence was attached to the 5′ end of each forward primer (tailed forward primer, Blacket et al., 2012). For microsatellite amplification, combinations of a tailed forward primer; a tail primer labeled with one of 6-FAM, VIC, PET, or NED (Life Technologies) dyes consisting of the complementary tail sequence; and a reverse primer were used. The multiplex PCR reactions were performed with four different marker sets in a total volume of 10 μL containing approximately 10 ng of template DNA using 1 × QIAGEN Type-it Microsatellite PCR Kit (QIAGEN) according to the manufacturer's standard protocol. PCR cycles consisted of a denaturing step of 15 min at 95°C; followed by 30 cycles of 94°C for 30 s, 57°C for 90 s, and 72°C for 60 s; and a final elongation step of 30 min at 60°C. PCR products were determined by capillary sequencer (3500xl Genetic Analyzer, Life Technologies). The allelic composition of each marker was determined for each accession, and putative alleles were indicated by the estimated size in base pair by GeneMapper (Life Technologies).

Genetic Diversity Analysis

A total of 319 accessions of the BNGA collection were genotyped. One accession was not included due to low yield of amplicon. The botanical race information on the 165 accessions had been previously determined by a genebank curator regarding the fruit-skin characteristics and anise-like odor of the leaf. Number of total alleles scored (A), observed heterozygosity (Ho), unbiased expected heterozygosity (uHe), and fixation index (F) were calculated to determine the genetic diversity and the past inbreeding. Allele frequencies of each botanical race and deviation from Hardy-Weinberg equilibrium (HWE) of loci were calculated for the subset of genotyping data on 165 accessions with botanical race information. All these values were calculated using GenAlEx 6.5 (Peakall and Smouse, 2012). Allelic richness (AR) and private allelic richness (PAR) for each botanical race were measured via rarefaction by HP-Rare v.1.0 (Kalinowski, 2005), calculated based on a minimum sample size of 12 (sample number of americana). Linkage disequilibrium (LD) between pairs of markers was tested by Genpop on the Web (Raymond and Rousset, 1995; Rousset, 2008; http://genepop.curtin.edu.au/).

Race Assignment and Genetic Group Classification of P. americana in BNGA

The membership compositions of each cluster, inferred by STRUCTURE analysis (Pritchard et al., 2000), were compared with the subset of 165 accessions with botanical race classification defined by morphological characteristics. The analysis was carried out to infer the most-likely number of genetic clusters (K) for the whole set of accessions showing both admixture and non-mixture by correlated allele frequency models.

Firstly, we screened numbers of clusters (K) from 1 to 8; each K was simulated 5 times with a burn-in of 10,000 iterations before collecting data, and 10,000 iterations of the Markov-Chain Monte Carlo (MCMC) method. Secondly, to obtain accurate results, parameter estimates were increased to a burn-in of 100,000 iterations and 100,000 iterations of MCMC with K from 1 to 5 with 5 iterations. The optimum K value was determined from log probability of data at each step of the MCMC, Pr(X|K), and ad hoc statistic ΔK of Evanno et al. (2005) using STRUCTURE HARVESTER (Earl and von Holdt, 2012). Genetic groups determined by STRUCTURE were then projected on a map using data from the 299 accessions with information on geographical origin. Geographical coordinate data was obtained from passport data, either as a georeferenced position or based on the localities of collection. For accessions that only had country of origin, the information was represented as the capital city's geographical coordinates.

Selection and Validation of Core-Collection

PCA-Kmeans method (Borrayo et al., 2016) was used for CC selection. The CC is intended to select a small subset of individuals of the National Avocado Germplasm Depository (BNGA) for long-term conservation by in vitro culture and cryopreservation. Only genotyping data were used as input in order to retain as much genetic variation as possible. We arbitrarily set the CC-number to 36 accessions (12.1% of the original collection), in what we consider to be a manageable sub-collection and in agreement with previous reports that establish that CC should contain ≥10% elements from the original collection (Guo et al., 2014), and then evaluated along with different CC sizes: 12, 24, 48, 72, and 96. Input data were depurated by the elimination of monomorphic alleles and accessions that lacked complete information. Twenty accessions were eliminated due to missing data from a total of 318 accessions genotyped, resulting in 298 accessions with complete genotyping data that served as the original collection for the CC selection procedure. Subsequently, the genotyping matrix with fragment sizes was transformed into a 0-to-1 scale and the PCA-Kmeans CC selection method was applied to the new data matrix. All CC selections and mathematical evaluations were performed with CorColv2.1-beta (program is available on request), an MS-Windows executable file built from a Python implementation of original FREEMAT codes published by Borrayo et al. (2016). A phenetic tree was constructed based on the same data using the neighbor-joining method (Saitou and Nei, 1987) in PowerMarker 3.25 (Liu and Muse, 2005) with 1,000 bootstrap replications. Distributions of CC members on the dendrogram were visualized by MEGA 5.2 (Tamura et al., 2011). The selected CC was validated by both mathematical evaluation parameters and their distribution patterns on the dendrogram; as well as by race information, geographical origin, and ratio of genetic group determined by STRUCTURE analysis.

Results

Genetic Diversity of BNGA Collection

The 28 markers that were analyzed (9 EST-SSR and 19 SSRs) detected a total of 547 alleles ranging from 4 to 34 alleles with an average of 19.1 per locus, out of which, 393 (71.7%) were rare alleles present in frequency less than 0.05. Except for LMAV20, all the markers tested were not under HWE (Table 2). Deviations from HWE of each marker were also tested by botanical race groups. Twenty-seven loci were not in HWD in drymifolia, 16 in guatemalensis, one in americana, and none in constericensis (Table 3). Out of 378 observed marker combinations, we detected significant (p < 0.01) LD between pairs of markers in 344 (91.0%) combinations in all 318 accessions. The botanical race-wise values for significant LD were 205 (54.2%), 125 (33.1%), and 5 (1.3%) combinations for drymifolia, guatemalensis, and americana, respectively. For costaricensis, only six marker combinations were valid to calculate LD due to the small number of samples, none of which was significant.

TABLE 2
www.frontiersin.org

Table 2. Diversity parameters associated with the 318 accessions of avocado analyzed.

TABLE 3
www.frontiersin.org

Table 3. Result of HWE statistics of each botanical race.

Analysis of the data subset of 165 accessions with botanical race information provided genetic diversity for each botanical race (Supplementary Figure 1). Values of observed heterozygosity, expected heterozygosity and number of private alleles are shown in Table 4. Allelic richness (AR) and private allele richness (PAR) in the population were: 5.95 and 0.63 for drymifolia, 6.13 and 0.89 for guatemalensis, 6.22 and 0.69 for americana, and 3.82 and 0.90 for costaricensis, respectively (Table 4).

TABLE 4
www.frontiersin.org

Table 4. Summary statistics for each botanical race.

Race Assignment

Genetic assignment methods implemented in STRUCTURE (Pritchard et al., 2000) revealed that K = 2 showed the highest likelihood in Mexican avocado by both the log likelihood [LnP(D)] and the Evanno's ΔK. Estimate of population divergence of allele frequencies in these two inferred populations from ancestral frequencies was 0.14 (Cluster A) and 0.05 (Cluster B), respectively.

The values represent the calculated average of 5 repeated runs. Average of a net nucleotide difference between clusters was 0.12. Average distance between individuals was 0.79 (Cluster A) and 0.71 (Cluster B). Accessions were assigned to one of the clusters when more than 80% of the genetic background belonged to such cluster. The number of accessions assigned to each cluster was 71 (22.3%) for Cluster A and 188 (59.1%) for Cluster B. Fifty-nine (18.6%) individuals were considered to be of admixed origin (80> A and B >20) (Supplementary Table 1).

Comparison of these genetic groups assigned by STRUCTURE and their botanical race information did not coincide with each other (Supplementary Figure 2). Therefore, we could not assign the botanical races based on genetic groups determined by STRUCTURE. Precise comparisons were impossible due to the disproportional presence of drymifolia compared with other botanical races. Nevertheless, individuals belonging to Cluster B (72) occur about five times more often than those in Cluster A (15) for drymifolia. Geographical distribution of Clusters A and B was represented on a map (Figure 1). The distribution of both clusters overlapped with some differences: distribution of Cluster A stretches wider in longitudinal range from the Yucatan Peninsula to the Pacific Coast, whereas distribution of Cluster B is more concentrated in the Central Highland.

FIGURE 1
www.frontiersin.org

Figure 1. Geographical distribution of assignment results of STRUCTURE (K = 2) analysis. Only the accessions assigned more than 80% to one of the clusters (A or B) were plotted on the map. The red circle represents Cluster A and the blue triangle, Cluster B. Distribution of the accessions originated in Mexico and Guatemala are shown.

Core Collection

As an adequate evaluation threshold was archived at CC36 and above (CC48, 72, and 96), and being 36 the smaller CC to meet this criterion, the core reference set of 36 individuals was selected from the 298 individual accessions. Scores obtained from the evaluation of CC36 were compared with those of other CCs as well as with the original collection (Table 5). In addition to individuals from three botanical races of P. americana, two other individuals were selected to CC: one belonging to the same sub-genus, P. americana var. nubigena, and the other from P. longipes, which belongs to sub-gender P. eriodaphne.

TABLE 5
www.frontiersin.org

Table 5. Core collection evaluation values for different K Core collection sizes.

A smaller ANE (average distance between each original collection and the nearest CC sample) means the diversity in the original collection is homogeneously represented by the CC. On the other hand, ENE (average distance between each CC sample and nearest CC sample) and E (average distance between CC samples) indicate the dispersion of data for the CC, where higher values stand for a better representation of extreme values. Values of ANE decrease as size of the CC increases; this is because members of a larger CC fill the gap of the smaller CC decreasing the average distance between the selected elements of the CC and those of the original collection.

ENE's largest value in CC36 indicates that this CC set is the one with the largest dispersion. The E-value increases as the number of CC increases; this indicates that extreme individuals are being included as the extreme values begin to form a single cluster.

Allele coverage (CA) reaches over 80% when CC size exceeds 36 and will continue to increase as CC becomes larger. This is explained by the large proportion of rare alleles in the total number of alleles detected (71.7% of total alleles); the larger the CC, the higher the number of rare alleles included in it. The frequency of rare alleles in CC36 was 197/352 (56.0%). Members of CC36 were evenly distributed on the dendrogram drawn by neighbor-joining method using the same dataset for CC selection. Compared with the distribution of a smaller number of the CC, CC36 members filled the non–represented branches by those CC12 and CC24 selected members (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. Position on the NJ dendrogram of accessions selected for core collection (K = 12, K = 24, and K = 36).

The distribution patterns of the 298-original collection and CC36 were compared with geographical origin, botanical race distribution, and ratio of genetic group determined by STRUCTURE analysis. The original collection consists of accessions collected from 17 States in Mexico and 7 other countries. Elements selected in CC36 were collected from 9 Mexican States and one other country (Costa Rica), encompassing the most representative area of distribution of the geographic origin of accessions in the original collection. Components of botanical races and species in the original collection and CC36 were 104 (32.6%) and 11 (30.6%) for drymifolia; 15 (4.7%) and 3 (8.3%) for guatemalensis; 12 (3.8%) and 1 (2.8%) for americana; 3 (0.9%) and 2 (5.6%) for costaricensis; 31 (9.7%) and 3 (8.3%) for hybrids; as well as 148 (46.4%) and 14 (39%) for unknown. The ratio of genetic groups determined by STRUCTURE analysis of the original collection and CC36 were almost the same. A list of the 36 accessions is found in Supplementary Table 1.

Discussion

Genetic Characteristics and Diversity of the BNGA Collection in CEBAJ-INIFAP

Mexico is the center of origin of avocado and there exist several germplasm collections and conservation activities regarding the species (Ashworth et al., 2011). However, there is a limited number of studies about genetic diversity and background of Mexican varieties. Among these studies, a global evaluation of avocado germplasm conserved in USA (Schnell et al., 2003) included 51 accessions of drymifolia. Cuiris-Pérez et al. (2009) used ISSR markers from 77 accessions of a nationwide avocado collection in Uruapan, Michoacan and reported high variation among varieties while focusing specifically on the Mexican race (drymifolia). Galindo-Tovar et al. (2011) also studied genetic relationships of avocado in Mexico using microsatellites and claimed to have found two genetic groups; however, the samples analyzed belonged to a considerably limited geographic area.

In this work, a considerable number of accessions (318 accessions) from the BNGA collection were used to illustrate the genetic diversity and background of avocado present in Mexico. As the BNGA collection in CEBAJ-INIFAP contains accessions from a wide geographical range in the country, it served as an optimal data source to infer genetic diversity and background of avocado present in Mexico. The result of this study demonstrated that this collection has achieved as much genetic diversity as other international avocado collections. Expected genetic diversity (He) values of the BNGA collection (He = 0.75) were smaller than those of 42 accessions of avocado varieties preserved in several institutes in Spain (He = 0.831), which consist of a mixture of drymifolia, guatemalensis, and americana collected from a broad geographical range (Gross-German and Viruel, 2013). The total average of the detected number of alleles at the BNGA collection (Na = 18.9, Table 2) exceeds that of the Spanish collection (A = 11.4). This high genetic diversity of the BNGA collection has derived mainly from the presence of rare alleles, which account for 71.7% (393/548) of total alleles. These rare alleles give significant values to Mexican avocado and emphasize the importance of conservation.

High Linkage Disequilibrium in Avocado in Mexico

Human imposed selection, intentional, or otherwise, modifies genetic diversity as it limits the number of lineages that are maintained for propagation. Although avocado is a subtropical tree species which predominantly outcrosses, we observed a deviation from HWE as well as a significant level of LD between pairs of markers in the materials used in this study. This could be explained by the extensive number of alleles being analyzed, yet it may suggest a domestication bottleneck process experienced differently by each avocado botanical race. Use of unlinked or weakly-linked genetic regions is recommended for STRUCTURE analysis to make meaningful inferences (Falush et al., 2003). Our results, especially drymifolia, demonstrated LD and substantially high deviation in all loci from HWE, which should be considered while interpreting these data. Chen et al. (2008) studied LD in wild avocado by sequencing 4 nuclear loci. They reported that significant excess of interlocus LD was observed when the three botanical races (drymifolia, guatemalensis and americana) were analyzed together, but not when the analysis was performed within each botanical race. They concluded that the LD arose from the genetic structure of the different botanical races used in their study.

Although in small proportion, compared to drymifolia (104), our initial analysis included a small number of individuals of guatemalensis (15), americana (12), and costarricensis (3), as well as some non-P. americana species. We contrasted this analysis with one that only included drymifolia and no significant difference was found in this comparison (data not shown). The same approach was carried out with CC36 and, although the LD decreases at several loci, the overall high LD tendency prevailed. Microsatellite markers used in this study consisted in est-SSR and genomic SSR markers. Both markers suggested interlocus LD, and showed no evidence of selection at the est-SSR loci corresponding to the markers analyzed in this study. Currently, genome sequencing of Mexican avocado is in process (Ibarra-Laclette et al., 2015). The complete genomic information will allow us to carry out further studies regarding interlocus LD and association mapping.

Implication of Avocado Domestication by Population Structure Analysis

Previous studies succeeded in grouping based on the botanical races using genetic classification by STRUCTURE analysis (Gross-German and Viruel, 2013). In this study, STRUCTURE analysis and clustering pattern based on utilized markers suggest that there are two genetic groups within the BNGA collection that did not coincide with determination of botanical races. A previous study based on ISSR analysis also reported two genetic groups in botanical race drymifolia (Cuiris-Pérez et al., 2009). Chen et al. (2009) evaluated haplotype analysis by sequencing 4 nucleic genes and reported that two genetic groups, which were determined by longitude and altitude, exist in wild avocados of Mexico. Torres-Gurrola et al. (2009) evaluated foliar chemical diversity using 35 accessions conserved in BNGA in CEBAJ-INIFAP, which partly coincide with materials used in this study, but the chemical profile did not demonstrate clear correlation between geographic distributions. The grouping pattern of our study partly coincides with results of the Martínez-Villagomez et al. (2016) analysis of the distribution of genus Persea based on climate parameters. The distribution of our Cluster A coincides with their group I (humid semi-warm to humid semi-cold), and our Cluster B with their group II (humid semi-warm to hot semi-dry). This congruence between genetic clustering and climate distribution clustering might imply that avocado in Mexico may have evolved separately in two climatic regions, which could have also led to differences in both est- and genomic microsatellite markers.

Other than the geographical/environmental factors, this grouping may have some relation with possible multiple domestication of avocado in Mexico. Clegg et al. (1993) suggest that several domestication events for avocado occurred in the past. Gama-Campillo and Gomez-Pompa (1992) concluded that avocado is a semi-domesticated tree and is still under the process of domestication by frequent exchange of trees/seedlings between wild habitat and home garden/orchard. Very low values of divergence of allele frequencies, especially found in Cluster B (0.05), suggests that the materials used in this study have also undergone such process.

Since genetic flow persists both among cultivated botanical races and non-cultivated wild avocadoes, correlation between race-identification based on morphological characteristics (Table 1) and genotypic information is difficult to achieve in Mexican avocado. Based on their morphological characteristics, 31 accessions were preliminarily identified as hybrids between botanical races. However, results rendered by STRUCTURE analysis identified that 59 accessions possess admixed genetic background. It is known that the complexity of the hybrid status (e.g., multiple backcrossing), and segregation of a race-specific trait (e.g., green or black skin) in hybrid origin progenies have made determination of botanical assignments extremely difficult (Ashworth and Clegg, 2003).

Molecular studies in clustering analysis have provided general agreement between botanical-race classification and the employed molecular markers with some exceptions (Ashworth and Clegg, 2003; Schnell et al., 2003; Alcaraz and Hormaza, 2007; Gross-German and Viruel, 2013). However, these studies are primarily based on improved varieties that have trackable pedigree records. Our study, on the other hand, is mainly based on local cultivars in the center of origin of the species, where gene flow is to be expected among their wild counterparts. Therefore, genealogical relationships are harder to determine and may also explain why botanical races do not coincide with the selected genotype-based clustering.

Core Collection and Long-Term Conservation of the BNGA Collection

Long-term conservation of commercially important tree species is one of the most urgent and challenging strategies aimed to assure the availability of such resources for future generations. This becomes particularly difficult with recalcitrant species. In the present work, we developed a core collection from 298 genotypes conserved in the BNGA in CEBAJ-INIFAP by PCA-K means clustering method with 28 microsatellite markers.

The objective of the study was to identify a core-collection that contained genetic diversity representative of the original collection. Agro-morphological traits were excluded as a selection criterion since their inclusion would have impacted the genotypic selection process, possibly affecting the desired optimal allelic representation. Despite this exclusion, the selected CC ended up successfully representing the original collection in terms of geographical distribution, representation of botanical races, and genetic group assignment determined by STRUCTURE analysis. However, specifically designed CCs may be selected for other purposes based on different selection criteria.

Currently, the BNGA collection is maintained as field trees. This conservation method presents possible threats by natural disaster, as pest and diseases are inevitable. Although a backup conservation approach has been already implemented among field collections in different localities, such collections are still at potential risk. As an alternative, recent advances in in vitro and cryopreservation techniques have proven effective in safeguarding recalcitrant species with a biotechnological approach (González-Arnao et al., 2014). Cryopreservation is an efficient and less expensive alternative, becoming even more so when applied to core collections, hence the importance of an adequate CC selection process.

Conclusion

The analysis of the BNGA avocado collection by molecular markers did not establish a clear difference among avocado botanical races. Our results suggest that the two genetic groups inferred were admixed and have contributed to the development of the current genetic structure of such populations, which might present evidence for putative environmental and/or domestication history of the species. This information, combined with agro-morphological characteristics, could prove useful when used for association studies, important in the acceleration of breeding procedures.

The selected core collection efficiently represents the diversity of the original collection and is a suitable candidate set for long-term cryopreservation conservation. This strategy may serve as a model for the conservation of other important recalcitrant species.

Author Contributions

RM conceived and designed the experiments; EH, ME, MC, and EB provided reagents/materials/analysis tools; LG and RM performed the laboratory procedures; RM, LG, EB, and ME analyzed the data; LG and RM wrote the paper; EH, ME, MC, and EB edited and provided critical review of the manuscript. All authors read and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We acknowledge Prof. Dr. Kazuo Watanabe, Dr. José Fernando de la Torre Sánchez, and Dr. Makoto Kawase for their project framework planning. We also thank Mr. Maharshi Ledezma Rodriguez and Mr. Juan Ramon Reynoso for their laboratory assistance. This research was supported by JST/JICA, Science and Technology Research Partnership for Sustainable Development (SATREPS) “Diversity Assessment and Development of Sustainable Use of Mexican Genetic Resources,” and Project INIFAP Number 13401032577.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.00243/full#supplementary-material

References

Alcaraz, M. L., and Hormaza, J. I. (2007). Molecular characterization and genetic diversity in an avocado collection of cultivars and local Spanish genotypes using SSRs. Hereditas 144, 244–253. doi: 10.1111/j.2007.0018-0661.02019x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashworth, V. E., Chen, H., and Clegg, M. T. (2011). “Persea,” in Wild Crop Relatives: Genomic and Breeding Resources Tropical and Subtropical Fruits, ed C. Kole (Berlin: Springer-Verlag), 173–189.

Google Scholar

Ashworth, V. E., and Clegg, M. T. (2003). Microsatellite markers in avocado (Persea americana Mill.): genealogical relationships among cultivated avocado genotypes. J. Hered. 94, 407–415. doi: 10.1093/jhered/esg076

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashworth, V. E., Kobayashi, M., De La Cruz, M., and Clegg, M. (2004). Microsatellite markers in avocado (Persea americana Mill.): development of dinucleotide and trinucleotide markers. Sci. Hortic. 101, 255–267. doi: 10.1016/j.scienta.2003.11.008

CrossRef Full Text | Google Scholar

Ayala, T., and Ledesma, N. (2014). “Chapter 8. Avocado history, biodiversity and production,” in Sustainable Horticultural Systems, Sustainable Development and Biodiversity, ed D. Nandwani (Springer International Publishing), 157–205.

Google Scholar

Ben-Ya'acob, A., Solis-Molina, A., and Bulfler, G. (2003). “The mountain avocado of Costa-Rica. Persea Americana var. Costericensis, A new sub-species,” in Proceedings V World Avocado Congress (Granada-Málaga), 27–33.

Blacket, M. J., Robin, C., Good, R. T., Lee, S. F., and Miller, A. D. (2012). Universal primers for fluorescent labelling of PCR fragments—an efficient and cost-effective approach to genotyping by fluorescence. Mol. Ecol. Resour. 12, 456–463. doi: 10.1111/j.1755-0998.2011.03104.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Borrayo, E., Machida-Hirano, R., Takeya, M., Kawase, M., and Watanabe, K. (2016). Principal components analysis - K-means transposon element based foxtail millet core collection selection method. BMC Genet. 17:42. doi: 10.1186/s12863-016-0343-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Borrone, J. W., Brown, J. S., Tondo, C. L., Mauro-Herrera, M., Kuhn, D. N., Violi, H. A., et al. (2009). An EST-SSR-based linkage map for Persea americana Mill. (avocado). Tree Genet. Genomes 5, 553–560. doi: 10.1007/s11295-009-0208-y

CrossRef Full Text | Google Scholar

Cañas-Gutiérrez, G. P., Galindo-López, L. F., Arango-Isaza, R., and Saldamando-Benjumea, C. I. (2015). Diversidad genetica de cultivares de aguacate (Persea americana) en Antioquia, Colombia. Agron. Mesoam. 26, 129–143. doi: 10.15517/am.v26i1.16936

CrossRef Full Text | Google Scholar

Chanderbali, A. S., Albert, V. A., Ashworth, V. E., Clegg, M. T., Litz, R. E., Soltis, D. E., et al. (2008). Persea americana (avocado): bringing ancient flowers to fruit in the genomics era. BioEssays 30, 386–396. doi: 10.1002/bies.20721

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H., Morrell, P. L., Ashworth, V. E., de la Cruz, M., and Clegg, M. T. (2009). Tracing the geographic origins of major avocado cultivars. J. Hered. 100, 56–65. doi: 10.1093/jhered/esn068

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H., Morrell, P. L., de la Cruz, M., and Clegg, M. T. (2008). Nucleotide diversity and linkage disequilibrium in wild avocado (Persea americana Mill.). J. Hered. 99, 382–389. doi: 10.1093/jhered/esn016

PubMed Abstract | CrossRef Full Text | Google Scholar

Clegg, M. T., Brandon, G. S., Duval, M. R., and Davis, J. (1993). Inferring plant evolutionary history from molecular data. N.Z. J. Bot. 3, 307–316. doi: 10.1080/0028825X.1993.10419508

CrossRef Full Text | Google Scholar

Cuiris-Pérez, H., Guillén-Andrade, H., Pedraza-Santos, M. E., López-Medina, J., and Vidales-Fernández, I. (2009). Genetic variability within Mexican race avocado (Persea americana Mill.) Germplasm collections determined by ISSRs. Rev. Chapingo Ser. Hortic. 15, 169–175. doi: 10.5154/r.rchsh.2009.15.023

CrossRef Full Text | Google Scholar

Earl, A., and von Holdt, B. M. (2012). STRUCTURE HARVESTER: a website and program for visualizing, STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 4, 359–361. doi: 10.1007/s12686-011-9548-7

CrossRef Full Text | Google Scholar

Efendi, D., and Litz, R. E. (2003). “Cryopreservation of avocado,” in Proceedings V World Avocado Congress (Granada-Málaga), 111–114.

Google Scholar

Evanno, G., Regnaut, S., and Goudet, J. (2005). Detecting the number of clusters of individuals using the software structure: a simulation study. Mol. Ecol. 14, 2611–2620. doi: 10.1111/j.1365-294X.2005.02553.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Falush, D., Stephens, M., and Pritchard, J. K. (2003). Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164, 1567–1587.

PubMed Abstract | Google Scholar

Frankel, O. H. (1984). “Genetic perspectives of germplasm conservation,” in Genetic Manipulation: Impact on Man and Society, eds W. K. Arber, K. Llimensee, W. J. Peacock, and P. Starlinger (Cambridge: Cambridge University Press), 161–170.

Google Scholar

Food Agriculture Organization of the United Nations (2014). Data from FAOSTAT Database. Available online at: http://www.fao.org/faostat/en/#data (Accessed on December 30, 2016).

Galindo-Tovar, M. E., Milagro-Pérez, P. A., Alejandro-Rosas, J. A., Leyva-Ovalle, O. R., Landero-Torres, I., Lee-Espinosa, H., et al. (2011). Genetic relationships within avocado (Persea americana Mill.) in seven municipalities of central Veracruz, characterized using microsatellite markers. Trop. Subtrop. Agroecosyst. 13, 339–346.

Galindo-Tovar, M. E., Ogata-Aguilar, N., and Arzate-Fernández, A. M. (2008). Some aspects of avocado (Persea americana Mill.) diversity and domestication in Mesoamerica. Genet. Resour. Crop Evol. 55, 441–450. doi: 10.1007/s10722-007-9250-5

CrossRef Full Text | Google Scholar

Gama-Campillo, L., and Gomez-Pompa, A. (1992). “An ethnoecological approach for the study of persea: a case study in the maya area,” in Proceedings Second World Avocado Congress (Orange City, FL), 11–17.

Google Scholar

González-Arnao, M. T., Martínez-Montero, M. E., Cruz-Cruz, C. A., and Engelmann, F. (2014). “Advances in cryogenic techniques for the long-term preservation of plant biodiversity,” in Biotechnology and Biodiversity, Sustainable Development and Biodiversity, eds M. R. Ahuja and K. G. Ramawat (Springer International Publishing), 120–170.

Google Scholar

Gross-German, E., and Viruel, M. A. (2013). Molecular characterization of avocado germplasm with a new set of SSR and EST-SSR markers: genetic diversity, population structure, and identification of race-specific markers in a group of cultivated genotypes. Tree Genet. Genomes 9, 539–555. doi: 10.1007/s11295-012-0577-5

CrossRef Full Text | Google Scholar

Guo, Y., Li, Y., Hong, H., and Qiu, L. J. (2014). Establishment of the integrated applied core collection and its comparison with mini core collection in soybean (Glycine max). Crop J. 2, 38–45. doi: 10.1016/j.cj.2013.11.001

CrossRef Full Text | Google Scholar

Huang, Q. X., Wang, X. C., Kong, H., Guo, Y. L., and Guo, A. P. (2013). An efficient DNA isolation method for tropical plants. Afr. J. Biotechnol. 12, 2727–2732. doi: 10.5897/AJB12.524

CrossRef Full Text | Google Scholar

Ibarra-Laclette, E., Méndez-Bravo, A., Pérez-Torres, C. A., Albert, V. A., Mockaitis, K., Kilaru, A., et al. (2015). Deep sequencing of the Mexican avocado transcriptome, an ancient angiosperm with a high content of fatty acids. BMC Genomics 16:599. doi: 10.1186/s12864-015-1775-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalinowski, S. T. (2005). HP-RARE 1.0: a computer program for performing rarefaction on measures of allelic richness. Mol. Ecol. Notes 5, 187–189. doi: 10.1111/j.1471-8286.2004.00845.x

CrossRef Full Text | Google Scholar

Kopp, L. E. (1966). A taxonomic revision of the genus Persea in the western hemisphere (Perseae–Lauraceae). Mem. N.Y. Bot. Gard. 14, 1–120.

Google Scholar

Liu, K., and Muse, S. V. (2005). PowerMarker: an Integrated analysis environment for genetic marker analysis. Bioinformatics 21, 2128–2129. doi: 10.1093/bioinformatics/bti282

PubMed Abstract | CrossRef Full Text | Google Scholar

Martínez-Villagomez, M., Campos-Rojas, E., Ayala-Arreola, J., Barrientos-Priego, A. F., and Espíndola-Barquera, M. C. (2016). Diversidad y distribución del género Persea Mill. en México. Agroproductividad 9, 72–77.

Google Scholar

Peakall, R., and Smouse, P. E. (2012). GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research–an update. Bioinformatics 28, 2537–2539. doi: 10.1093/bioinformatics/bts460

PubMed Abstract | CrossRef Full Text | Google Scholar

Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics 155, 945–959.

PubMed Abstract | Google Scholar

Raymond, M., and Rousset, F. (1995). GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J. Hered. 86, 248–249.

Google Scholar

Rincón-Hernández, C. A., Sánchez Pérez, J. D., and Espinosa-García, F. J. (2011). Caracterización química foliar de los árboles de aguacate criollo (Persea americana var. drymifolia) en los bancos de germoplasma de Michoacán, México. Rev. Mex. Biodivers. 82, 395–412.

Google Scholar

Rousset, F. (2008). Genepop'007: a complete reimplementation of the Genepop software for Windows and Linux. Mol. Ecol. Resour. 8, 103–106. doi: 10.1111/j.1471-8286.2007.01931.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Saitou, N., and Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4, 406–425.

Google Scholar

Sarr, M., Goeschl, T., and Swanson, T. (2008). The value of conserving genetic resources for R&D: a survey. Ecol. Econ. 67, 184–193. doi: 10.1016/j.ecolecon.2008.03.004

CrossRef Full Text | Google Scholar

Schnell, R. J., Brown, J. S., Olano, C. T., Power, E. J., Krol, C. A., Kuhn, D. N., et al. (2003). Evaluation of avocado germplasm using microsatellite markers. J. Am. Soc. Hort. Sci. 128, 881–889.

Google Scholar

Scora, R. W., Wolstenholme, B. N., and Lavi, U. (2002). “Chapter 2. Taxonomy and botany” in The Avocado: Botany, Production and Uses eds A. W. Whiley, B. Schaffer, and B. N. Wolstenholme (Wallingford: CAB International), 15–37.

Sharon, D., Cregan, P. B., Mhameed, S., Kusharska, M., Hillel, J., Lahav, E., et al. (1997). An integrated genetic linkage map of avocado. Theor. Appl. Genet. 95, 911–921. doi: 10.1007/s001220050642

CrossRef Full Text | Google Scholar

Smith, C. E. J. (1966). Archeological evidence for selection in avocado. Econ. Bot. 20, 169–175. doi: 10.1007/BF02904012

CrossRef Full Text | Google Scholar

Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., and Kumar, S. (2011). MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony method. Mol. Biol. Evol. 28, 2731–2739. doi: 10.1093/molbev/msr121

CrossRef Full Text | Google Scholar

Torres-Gurrola, G., Montes-Hernández, S., and Espinosa-García, F. J. (2009). Patrones de variación y distribución geográfica en fenotipos químicos foliares de Persea americana var. drymifolia. Rev. Fitotec. Mex. 32, 19–30.

Google Scholar

Keywords: Persea americana, microsatellites, botanical races, core collection, PCA-Kmeans

Citation: Guzmán LF, Machida-Hirano R, Borrayo E, Cortés-Cruz M, Espíndola-Barquera MdC and Heredia García E (2017) Genetic Structure and Selection of a Core Collection for Long Term Conservation of Avocado in Mexico. Front. Plant Sci. 8:243. doi: 10.3389/fpls.2017.00243

Received: 11 November 2016; Accepted: 08 February 2017;
Published: 24 February 2017.

Edited by:

Rodomiro Ortiz, Swedish University of Agricultural Sciences, Sweden

Reviewed by:

Robert Lawrence Jarret, Agricultural Research Service (USDA), USA
Kaye Enid Basford, University of Queensland, Australia
Arun Jagannath, University of Delhi, India
Marinus J. M. Smulders, Wageningen University and Research Centre, Netherlands

Copyright © 2017 Guzmán, Machida-Hirano, Borrayo, Cortés-Cruz, Espíndola-Barquera and Heredia García. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ryoko Machida-Hirano, hiranoryoko@gmail.com

These authors have Contributed equally to this work.