Assessment of the Genetic Diversity of Almond (Prunus dulcis) Using Microsatellite Markers and Morphological Traits

Document Type : Research Paper


1 Department of Genomics, Agricultural Biotechnology Research Institute of Iran (ABRII), Mahdasht Road, P.O. Box 31535-1897, Karaj, I.R. Iran

2 Department of Agronomy and Plant Breeding, Faculty of Agriculture, University of Zanjan, P.O. Box 45195-313, Zanjan, I.R. Iran

3 Department of Agronomy and Plant Breeding, Faculty of Agriculture, University of Tehran, P.O. Box 31587-11167, Tehran, I.R. Iran


The genetic diversity among 56 almond (Prunus dulcis) genotypes was analysed using 35 microsatellite markers and 14 morphological traits. Analysis of morphological traits revealed a wide range of variation among the studied genotypes. Out of 35 simple sequence repeats (SSRs) markers, 25 were polymorphic, producing 215 alleles that varied from 2 to 16 with an average of 8.76 alleles per locus. Regression analyses revealed a positive correlation between the CPPCT03 locus and kernel yield, kernel percentage, grain weight, leaf length and tree altitude. The results of analysis of molecular variance (AMOVA) indicated that approximately 4.5% of genetic variance was observed between the collection sites. Based on SSR data, cluster analyses showed that the studied almond genotypes were classified into five main groups. The results of the present study showed that microsatellite markers could be successfully used to assay genetic diversity among Iranian almond landraces/cultivars and to identify informative markers for improving traits in breeding programs.



Almond [Prunus dulcis (Miller) D.A. Webb, syn. Prunus amygdalus Batsch] occupies a very peculiar place among fruit trees (Miller et al., 1989). Because of almond’s tolerance to cold, drought and salinity, it is considered an important tree crop and is cultivated in different climatic regions of Iran. Breeding practices in Prunus face unique challenges resulting from the narrow genetic background of commercial cultivars (Scorza et al., 1985). Morphological traits such as seed and kernel size, kernel yield, and blooming time are usually used for cultivar identification in almond (De iorgio and Polignano, 1999). However, morphological traits are limited because of their environmental fluctuations.
In recent years, molecular markers have been used to study genetic diversity and cultivar identification of peach and almond (Shiran et al., 2007; Sorkheh et al., 2007; Amirbakhtiar et al., 2006; Kadkhodaei et al., 2006; Sanchez-Pérez et al., 2006; Xie et al., 2006; Testolin et al., 2000, 2004; Aranzana et al., 2003). Methods based on knowledge provided by advances in molecular genetics, notably molecular markers, promise faster and more efficient approaches to cultivar improvement. In fact important tools such as molecular markers, maps, DNA sequences, and quantitative trait loci (QTLs) have been developed and made available to researchers, and applications at the breeding program level have already started (Dirlewanger et al., 2004). Recently, DNA microarray-based genome composition analysis has also been used in comparative genomic studies of trees (Martinez-Gomez et al., 2007). The objectives of the present study are to investigate the genetic diversity of major Iranian almond landraces/cultivars, to identify their relationship to important foreign cultivars, and to introduce informative markers for important nut traits using microsatellite markers.

Materials and Methods

Plant materials: Fifty-one Prunus dulcis landraces/cultivars from different provinces of Iran along with three and two registered cultivars from Spain and USA, respectively, were used in this study (Table 1). The trees with similar ages were sown in a randomized complete block design, with four replications, at the experimental field of the Agricultural Biotechnology Research Institute of Iran (ABRII), Isfahan.
Phenotypic analysis: Fourteen independent morphological traits including leaf shape, leaf length (cm), leaf width (cm), petiole length (cm), flowering duration (day), tree altitude (cm), frostbite kernel yield (g), kernel length (cm), kernel width (cm), kernel thickness (cm), nut weight (g), kernel nut weight (g) and kernel percentage were recorded, based on food and agriculture organization (FAO).

Microsatellite analysis: Total genomic DNA was extracted according to the method described by Doyle and Doyle (1987), with minor modifications. Thirty-five simple sequence repeat (SSR) markers, isolated from peach and almond were used in this study (Testolin et al., 2004; Dirlewanger et al., 2002). Amplification reaction products were separated on a 6% (w/v) denaturing polyacrylamide gel using a Sequi-Gen GT Sequencing Cell 30 cm gel apparatus (BioRad Laboratories Inc., Hercules, CA, USA). The amplified fragments were detected by the silver staining method as described by Bassam et al. (1991). The gels were visually scored by two independent observations.

Data analysis: Each polymorphic  fragment was scored as either present (1) or absent (0) across all genotypes. The data were used to calculate the similarity matrix among cultivars employing simple matching coefficients. The similarity matrix was then used to construct dendrograms using the unweighted pair group method with arithmetic averages (UPGMA). This was achieved by employing the sequential, agglomerative, hierarchical, and nested clustering (SAHN) using the numerical taxonomy and multivariate analysis system (NTSYS-PC), version 2.00 (Rohlf, 1998). Observed heterozygosity (Ho) and expected heterozygosity (He) were calculated using the POPGENE version 1.32 (Yeh et al., 1997). The degree of polymorphism was quantified using the polymorphic information content (PIC). Probability of identity (PI) was estimated according to Paetkau et al. (1995). Analysis of molecular variance (AMOVA) was performed using the Arlequin version 2.00 (Schneider et al., 2000) to determine genetic variation (Nei, 1972). Average value of the Shannon index was also measured (Shannon and Weaver, 1949). Informative markers were determined by stepwise regression using the SPSS software version 10.0 for windows (SPSS Inc., Chicago, IL).

Morphological trait analysis: Mean, maximum, minimum and the percentage of coefficient of variation (CV%) of 14 morphological characters are shown in Table 2. A large diversity in the characters was observed, indicating a high level of variation in the studied plant materials.

SSR marker analysis: The results of this study showed cross amplification ability of microsatellite markers among the studied almond genotypes. Out of 35 SSR markers, Out of 35 SSR markers, 25 were polymorphic and produced 215 alleles.  The number of alleles per locus ranged from 2 to 16, with an average of 8.76 (Table 3). Average value of the Shannon index was 1.79, which varied from 0.35 in UDP96-008 to 2.6 in CPPCT3. Mean He across microsatellite loci ranged from 0.92 in CPPCT3 to 0.17 in UDP96-008. The highest level of observed heterozygosity was found in XAM18 and CPPCT22 and the lowest in UDP96-008. According to PI, the most informative loci were UDP98-412 and CPPCT3 with values of 0.041 and 0.042, respectively. PIC for these two loci was greater (0.7) than others. The least informative locus was XAM04 with PI of 0.98 and PIC of 0.159, followed by XAM18 with PI and PIC values of 0.494 and 0.0018, respectively. The average of PI and PIC values for all loci were 0.258 and 0.475, respectively (Table 3). Rare polymorphic alleles (i.e. those with a frequency of £ 0.005) and their weights were determined for the purpose of rapid cultivar identification (Table 4). Regression analyses revealed that there was a positive correlation between the CPPCT03 locus and kernel yield (b = 0.424), kernel percentage (b = 0.49), grain weight (b = 0.35), leaf length (b = 0.32) and tree altitude (b = 0.327) (Table 5).
     Based on sampling sites, average He was 0.697 and the largest heterozygosity was observed for cultivars from Hamadan (0.731). The results of AMOVA indicated that approximately 4.5% of genetic variance belonged to between collected sites (Table 6). Based on SSR data, the studied almond genotypes were classified into five main groups (Fig. 1). The first cluster included some landraces and cultivars from the Shiraz, Isfahan, Hamadan and Arak provinces. The second cluster included two sub-clusters: the first sub-cluster contained 4 landraces from the Shiraz province and the second sub-cluster contained registered cultivars from Spain, USA and Azerbaijan. Two landraces from Shiraz and Arak provinces were gathered into cluster III. One registered cultivar from USA (HO) and one registered cultivar from Azerbaijan (Harir) were located in two distinct clusters (IV and V).

The results of this study support those of Sosinski et al. (2000), regarding the cross amplification ability of microsatellite markers across the Prunus species. High level of heterozygosity for all loci (0.697) can be attributed to cross pollination and the self-incompatibility nature of almond. The high values of polymorphic loci (71%), average number of alleles per locus (8.76), He (0.775), average polymorphism information content (0.475) and PI (0.258) observed in this study indicate that SSR markers are able to identify genetic variation among the studied almond genotypes. According to PI and PIC values, CPPCT3, UDP98-412, UDP96-409, XAM05, XAM08, XAM09, XAM15 and XAM19 are the best loci for further studies of almond genetic diversity. The percentage of polymorphic SSR loci (71%) in this study was much higher than that estimated for RFLPs (21.9%), suggesting that SSRs can act as better systems for almond cultivar identification (Eldredge et al., 1992).
During this research, alleles were identified that correlated with yield-related traits. The allele belonging to the XAM09 locus had a positive correlation with blooming duration (0.418) (Table 5). In addition, CPPCT17 was found to be an informative marker for nut weight, average kernel thickness and leaf width (Table 5).
In this investigation, cluster analyses showed that most Iranian landraces are well separated from the Spanish and American (USA) cultivars, indicating that they may be native to Iran. However, Shiraz almond landraces are assigned to the same group as the Spanish and American cultivars. A possible explanation is that they might carry a common genetic background. According to the results of this study, SSR data failed to separate genotypes based on their sampling sites. Germplasm migration or insufficient SSR markers can explain this incomplete separation. The results show that Iranian registered cultivars including Yalda, Shokofe and Sahand are similar to the foreign cultivars.
Informative markers are most applicable for breeding purposes. These markers have previously been used in the identification of peach and nectarine varieties (Manubens et al., 1999). A combination of molecular and morphological data is the best choice to find informative markers. In summary, results of the present study reveal that microsatellite markers can be successfully used to assay genetic diversity among Iranian almond landraces/cultivars and to identify informative markers for breeding of important traits.


We would like to thank the Agricultural Biotechnology Research Institute of Iran (ABRII) for funding this study.

Amirbakhtiar N, Shiran BH, Moradi H, Sayed-Tabatabaei BE (2006). Molecular characterization of almond cultivars using microsatellte markers. ISHS Acta Horticulture 726: IV International Symposium on Pistachios and Almonds. Tehran, Iran.
Aranzana MJ, Pineda A, Cosson P, Dirlewanger E, Ascasibar J, Cipriani G, Ryder CD, Testolin R, Abbott A, King GJ, Iezzoni AF, Arus P (2003). A set of simple-sequence repeat (SSR) markers covering the Prunus genome. Theor Appl Genet. 106: 819-825.
Bassam BJ, Caetano-Anollés G, Gresshoff PM (1991). Fast and sensitive silver staining of DNA in polyacrylamide gels. Anal Biochem. 196: 81-84.
De iorgio D, Polignano G (1999). Evaluating almond biodiversity of cultivars grown a germoplasm collection field in south of Italy. Proc. ISCO Conference, May 23-28, Purdue University, West Lafayette, Indiana, USA. 
Dirlewanger E, Cosson P, Tavaud M, Aranzana MJ, Poizat C, Zanetto A, Arus P, Laigret F (2002). Development of microsatellite markers in peach [Prunus persica (L.) Batsch] and their use in genetic diversity analysis in peach and sweet cherry (Prunus avium L.). Theor Appl Genet. 105: 127-138.
Dirlewanger E, Graziano E, Joobeur T, Garriga-Caldere F, Cosson P, Howad W, Arus P (2004). Comparative mapping and marker-assisted selection in Rosaceae fruit crops. PNAS. 101: 9891-9896.
Doyle JJ, Doyle JL (1987). A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull. 19: 11-15.
Eldredge L, Ballard R, Baird W, Abbott A, Morgens P, Callahan A, Scorza R, Monet R (1992). Application of RFLP analysis to genetic mapping in peaches. HortScience 27: 160-163.
Kadkhodaei S, Aghdaei SRT, Grigorian V, Moghadam M, Hashemi SMM (2006). A study on genetic variation among some wild almond species using RAPD markers. ISHS Acta Horticulture 726: IV International Symposium on Pistachios and Almonds. Tehran, Iran.
Manubens A, Lobos S, Jadue Y, Toro M, Messina R, Liadser M, Seelenfreund D (1999). DNA isolation and AFLP fingerprinting of nectarine and peach varieties (Prunus persica). Plant Mol Biol Rpt. 17: 255-267.
Martinez-Gomez P, Sanchez-Perez R, Dicenta F, Howad W, Arus P, Gradziel TM (2007). Genome mapping and molecular breeding in plants. Springer Berlin Heidelberg. 4: 229-242.
Miller PJ, Parfitt DE, Weinbaum SA (1989). Outcrossing in peach. HortScience 24: 359-360.
Nei M (1972). Genetic distance between collections. Amer Naturalist. 106: 283-292.
Paetkau D, Calvert W, Stirling I, Strobeck C (1995). Microsatellite analysis of population structure in Canadian polar bears. Mol Ecol. 4: 347-354.
Rohlf FJ (1998). NTSYSpc: Numerical Taxonomy and Multivariate Analysis System, Version 2.02, Exeter Software, Setauket, New York.
Sanchez-Pérez R, Dicenta F, Martinez-Gomez P, Howad W, Arus P (2006). Construction of linkage map and QTL analysis of agronomic traits in almon using SSR markers. ISHS Acta Horticulture 726: IV International Symposium on Pistachios and Almonds. Tehran, Iran.
Schneider S, Roessli D, Excoffier L (2000). Arlequin: A Software for Population Genetics Data Analysis, Version 2.000 Genetics and Biometry Laboratory, Dept. of Anthropology, University of Geneva, Switzerland.
Shannon CE, Weaver W (1949). The mathematical theory of communication. Urbana, Ill: Univer. of Illinois Press.
Shiran B, Amirbakhtiar N, Kiani S, Mohammadi SH, Sayed-Tabatabaei BT, Moradi H (2007). Molecular characterization and genetic relationship among almond cultivars assessed by RAPD and SSR markers. Sci Hortic. 111: 280-292.
Scorza R, Mehlenbacher SA, Lighner GW (1985). Inbreeding and coancestry of freestone peach cultivars of the eastern United States and implications for peach germplasm improvement. J Am Soc Hortic Sci. 110 : 547-552.
Sorkheh K, Shiran B, Gradziel, TM,  Epperson BK,  Mrtinez-Gomez P, Asadi E (2007). Amplified fragment length polymorphism as a tool for molecular characterization of almond germplasm: genetic diversity among cultivated genotypes and related wild species of almond, and its relationships with agronomic traits. Euphytica 156: 327-344.
Sosinski B, Gannavarapu M, Hager LD, Beck E, King GJ, Ryder CD, Rajapakse S, Baird WV, Ballard RE, Abbott AG (2000). Characterization of microsatellite markers in peach [Prunus persica (L.) Batsch]. Theor Appl Genet. 101: 421-428.
Testolin R, Marrazzo T, Cipriani G, Quarta R, Verde I, Dettori M, Pancaldi M, Sansavini S (2000). Microsatellite DNA in peach (Prunus persica L. Batsch) and its use in fingerprinting and testing the genetic origin of cultivars. Genome 43: 512-520.
Testolin R, Messina R, Lain O, Marrazzo MT, Huang WG, Cipriani G (2004). microsatellites isolation in almond from an AC-repeat enriched library. Mol Ecol Notes. 4: 459-461.
Xie H, Sui Y, Chang FQ, Xu Y, Ma RC (2006). SSR allelic variation in almond (Prunus dulcis Mill.). Theor Appl Genet. 112: 366-72.
Wright S (1951). The genetical structure of collections. Ann Eugen. 15: 323.
Yeh FC, Yang RC, Boyle T (1997). POPGENE version 1.21. (URL