nuclear and chloroplast dna variability and phylogeny of iranian apples (malus domestica)
TRANSCRIPT
ORIGINAL ARTICLE
Nuclear and chloroplast DNA variability and phylogeny of Iranianapples (Malus domestica)
Abdollah Khadivi-Khub • Shahin Jahangirzadeh •
Elahe Ahadi • Samad Aliyoun
Received: 25 February 2013 / Accepted: 2 February 2014
� Springer-Verlag Wien 2014
Abstract This study was conducted to reveal genetic
diversity among 23 local apple genotypes using nuclear
(RAPD) and chloroplast DNA (PCR-cpRFLP) markers.
Eleven RAPD primers and four cpDNA primer combina-
tions were used in this study. RAPD primers produced a
total of 77 polymorphic fragments with an average of seven
bands per primer. The percentage of polymorphic bands
(68.14 %) showed the efficiency of used RAPD primers in
distinguishing all the genotypes considered. Genetic simi-
larity between studied genotypes varied from 0.38 to 0.72
and cluster analysis showed the abundant diversity, indi-
cating high intraspecific genetic variation between Iranian
apple genotypes. From the four universal chloroplast pri-
mer pairs, three primer pairs amplified the fragments and
their combinations showed polymorphic patterns and
revealed intraspecific chloroplast variation. The informa-
tion will facilitate germplasm identification, conservation
and new cultivar development.
Keywords Nuclear and chloroplast DNAs � Cluster
analysis � Genetic diversity � Phylogenetic relationships �Malus domestica
Introduction
Apple (Malus 9 domestica Borkh.) is a member of Rosa-
ceae and one of the most widespread and popular fruit trees
in the world (Janick et al. 1996). About 59 species and
7500 cultivars were identified in all over the world. Apple
is an ancient fruit crop in Iran (Janick et al. 1996) and there
is an extremely abundant germplasm resource for it. There
is a high level of genetic diversity in Iran’s cultivated apple
due to closely distance to apple origin in Central Asia.
Apple is an economically important fruit crop and has a
broad spectrum uses as food; most apples are eaten raw but
can be made into jellies, pies, cakes and puddings. They
can also be canned, juiced and optionally fermented to
produce apple juice, cider, vinegar and pectin (Jackson
2003).
Knowledge of genetic variation and phylogenetic rela-
tionship among genotypes is an important consideration for
classification, utilization of germplasm resources and
breeding. Traditionally, cultivar identification has relied on
morphological and agronomic characteristics of plant
materials. Although there is substantial intraspecific vari-
ation in vegetative traits, especially leaf and fruit charac-
ters, it is difficult to distinguish genotypes on their external
morphology alone. Further, phenotypic characters are
generally influenced by environmental factors and the
growth stage of the plant. In fruit trees, this requires a
lengthy and expensive evaluation during the whole vege-
tative growth. There is an additional problem to distinguish
among cultivars which come from the same cross. In such
cases, the relatively narrow range of variation of morpho-
logical traits limits cultivar identification, and different
methods must be used (Khadivi-Khub et al. 2008).
The development of randomly amplified polymorphic
DNA (RAPD) markers, generated by the polymerase chain
A. Khadivi-Khub (&)
Department of Horticultural Sciences, Faculty of Agriculture and
Natural Resources, Arak University, 38156-8-8349 Arak, Iran
e-mail: [email protected]; [email protected]
S. Jahangirzadeh
Department of Horticultural Sciences, Faculty of Agriculture,
University of Guilan, Rasht, Iran
E. Ahadi � S. Aliyoun
Department of Horticultural Sciences, Faculty of Agriculture,
University of Tehran, 31587 Karaj, Iran
123
Plant Syst Evol
DOI 10.1007/s00606-014-1007-y
reaction (PCR) using arbitrary primers, has provided a new
tool for the detection of DNA polymorphisms (Welsh and
McClelland 1990; Williams et al. 1990). The RAPD ana-
lysis has already proved to be valuable in apple genome
analysis. These markers have been used for identification
of apple cultivars (Koller et al. 1993; Mulcahy et al. 1993)
and apple rootstocks (Autio et al. 1998) and for paternity
analysis (Harada et al. 1993; Gardiner et al. 1996). Luis
et al. (2001) employed RAPD in discrimination and esti-
mation of genetic similarities among apples cultivars,
while Zhou and Li (2000) used RAPD for the phylogenetic
relationships of the closely related species of cultivated
apple. In most cases, data on genetic similarity obtained by
RAPD analysis matched classifications based on morpho-
logical traits in apple (Landry et al. 1994).
Chloroplast DNA (cpDNA) shows maternal inheritance
in most flowering plants, does not undergo recombination
and has a slow rate of evolution (Palmer 1987); therefore,
cpDNA is a good tool for tracing seed parents and maternal
ancestors. The detection of restriction fragment length
polymorphisms (RFLPs) in specific regions of cpDNA,
amplified by the polymerase chain reaction (PCR), has
been developed as a method for detecting variations in
cpDNA. There are many other horticultural plants whose
backgrounds are unclear, and PCR–RFLP analysis has
been used in horticultural cultivars of chrysanthemum
(Kishimoto et al. 2003), rose hybrids (Takeuchi et al. 1995)
and sweet cherry (Mohanty et al. 2001; Panda et al. 2003).
This method also has been used for analysis of apple
genotypes and phylogenetic relationships in wild and
domesticated apples (Coart et al. 2006).
In this study, we used RAPD and PCR-cpRFLP markers to
characterize nuclear and chloroplast genome in apple geno-
types from different regions of Iran, for both commercial and
selected clones from the breeding program. Molecular
markers were used to study the genetic variation and phylo-
genetic relationships among local apple genotypes.
Materials and methods
Plant material and DNA isolation
Twenty-three local apple genotypes from different geo-
graphical and/or agroclimatic locations in Iran were used in
this study (Table 1). These genotypes are representing
genetic diversity of three main regions of apple in Iran. We
selected these materials with high severity and accuracy. It
is necessary to note that these genotypes have high mor-
phological variation, and their resistance to chilling and
Table 1 The studied apple
genotype names, codes and their
regions
Genotype Region Code of genotype Origin Flavor
T3 Central Alborz T3 Iran Sweet
Ferdos Central Alborz Ferdos Iran Sweet
Shahrod10 Central Alborz Shah10 Iran Sweet
Shahrod19 Central Alborz Shah19 Iran Sweet
Golab-Damavand Central Alborz GolDa Iran Sweet
Golab-Paize Central Alborz GolPA Iran Sweet
Golab-Nemati Central Alborz GolNe Iran Sweet
Azerbaijan6 Azerbaijan Azer6 Iran Sweet
Azerbaijan7 Azerbaijan Azer7 Iran Sweet
Azerbaijan14 Azerbaijan Azer14 Iran Sweet
Beyghi Azerbaijan Beyghi Iran Sweet
Mahali-Urmia Azerbaijan Mahali Iran Sweet–sour
ME 3 Azerbaijan ME 3 Iran Sweet
ME 4 Azerbaijan ME 4 Iran Sweet
Ghara-Yarpagh Azerbaijan Ghara Iran Sweet
Khan-Almasi Azerbaijan Khan Iran Sweet
Meshkin-Shahr Azerbaijan Mesh Iran Sweet
SBA Central Zagros SBA Iran Sweet
Golab-Arak Central Zagros GolAr Iran Sweet
ShahreKord8 Central Zagros ShahKord8 Iran Sweet
Torsh-Kermanshah Central Zagros TorKer Iran Sour
Ghermez1 Central Zagros Gherm1 Iran Sweet
Zodras-Oghlid Central Zagros ZodOgh Iran Sweet
Granny-Smith Australia Granny Australia Sweet–Sour
A. Khadivi-Khub et al.
123
drought is different. Also, flowering and ripening times in
these genotypes are varied. Furthermore, one commercial
cultivar, Granny-Smith, was used as an outgroup. Young
leaves were collected, frozen in liquid nitrogen and subse-
quently stored at -20 �C until processed. Genomic DNA
(nuclear and chloroplast) was isolated from leaf samples
using a CTAB extraction method (Doyle and Doyle 1987).
The extraction buffer contained 2 % CTAB, 1.5 M NaCl,
20 mM EDTA, 100 mM Tris pH 8.0, 2 % soluble poly
vinyl pyrrolidone (PVP), and 2 % 2-mercaptoethanol. DNA
quality was examined by electrophoresis in 0.8 % agarose,
and DNA concentration was quantified on a GeneQuant
(Amersham Pharmacia Biotech).
RAPD analysis
For the RAPD reactions, 20 ng of DNA was used as tem-
plate in a final volume of 25 ll containing 1x reaction
buffer (20 mM (NH4)2SO4, 75 mM Tris–HCl pH 8.8,
0.01 % Tween 20), 1.9 mM MgCl2, 200 lM of each dNTP,
0.4 lM primer and 0.6 U Taq DNA polymerase, overlaid
with a drop of mineral oil. An extra reaction containing all
the components, plus water instead of template DNA, was
included in all the experiments. The DNA amplifications
were performed in a thermocycler (iCycler, Bio Rad Co.,
USA) as follows: 1 cycle of 5 min at 93 �C, 45 cycles of
1 min at 93 �C, 1 min at 36 �C, and 2 min at 72 �C (for
denaturing, annealing and primer extension, respectively).
The last cycle was followed by a final incubation for 10 min
at 72 �C, and the PCR products were stored at 4 �C before
analysis. The DNA amplification products were analyzed
by electrophoresis in 1.5 % agarose gels (Roche Co., Ger-
many) in 1x TBE buffer for 5 h at 3.5 V/cm. The gels were
stained for 20 min with ethidium bromide (1 lg/ml) and
destained in water for 15 min and photographed under UV
light, by a Gel Doc system (UVP, Bio Doc Co., USA).
Fragment length was estimated by comparison with stan-
dard size markers (100 base pair ladder).
Chloroplast analysis (PCR-cpRFLP)
Extracted DNAs were PCR-amplified using four universal
primer pairs for chloroplast DNA (K1K2, CS, HK and TF)
(Dumolin-Lapegue et al. 1997; Grivet et al. 2001). The
amplification of cpDNA regions was performed in 30 ll of
reaction mixture containing 30 ng of total DNA, 10 ng of
each primer pair, 200 lM of each four dNTPs, one unit
Taq DNA polymerase (recombinant, Fermentas, Canada),
2 mM of MgCl2 and 1X PCR buffer with KCl. The ther-
mocycle system used was Gene Amp PCR system 9700.
PCR was carried out as described by Mohanty et al.
(2001), using an initial cycle of 4 min at 94 �C, followed
by 30 cycles of 45 s at 94 �C, 45 s at 54.5 �C to 58.5 �C
and 2 or 3 min at 72 �C (annealing temperature and
extension time depending on the primer and the length of
the fragment to be amplified) and finally a 10 min exten-
sion at 72 �C. Two restriction enzymes (EcoRI and MseI;
Fermentas, Canada) were used for digesting 5 ll PCR
products in a volume of 30 ll containing 10 unit of
restriction enzyme following the manufacturer’s instruc-
tions (each reaction mixture was incubated for 10 h).
cpDNA restriction fragments were separated in 3 % aga-
rose gels cast in a Tris–borate-EDTA buffer (1 9), by
running at 70 V for 4 h, and after that, fragments were
visualized in UV light. Each mixture loaded in agarose gel
contained 10 ll of restriction mixture, 3 ll of loading
buffer and 2 ll Gel Red (Biotium, USA). The size markers
used for analyzing the size of polymorphic bands were
1 kb and 50 bp ladders (Fermentas, Canada).
Data analysis
Data were recorded as presence (1) or absence (0) of
amplified products from the examination of photographic
negatives. The total number of fragments scored, the number
of polymorphic fragments and the percentage of polymor-
phic fragments were determined for each primer used.
Genetic similarity based on RAPD and chloroplast data was
calculated by making a pairwise comparison between all
genotypes according to the Jaccard coefficient using the
Simqual module of NTSYS-pc software version 2.01e
(Rohlf 1998). Similarity matrices were compared using the
Mantel matrix-correspondence test (Mantel 1967) using
MxComp module in NTSYS-pc ver 2.1 (Rohlf 2000).
Bootstrap analysis with 1000 replicates was performed using
FreeTree software (version 0.9.1.50, Pavlicek et al. 1999) to
obtain the confidence of branches of the UPGMA and
Neighbor-joining (NJ) trees and validation of the dendro-
grams. Trees were viewed using Tree View program, version
1.6.6 (Page 1996). The bootstrap analysis estimates the
probability values (P values) for each cluster (the probability
of a true cluster at each edge of the dendrogram). Also, Dollo
Parsimony was used for phylogenetic relationships by PUAP
software. Genetic parameters such as total diversity (Ht),
average intraspecific diversity (Hs) and level of population
subdivision (GSTC) were estimated using POPGENE ver.
1.31 (Yeh et al. 1999) to determine Nei’s genetic diversity
(He) and Shannon’s information index (I) amongregions.
Results and discussion
RAPD analysis
Primer selection is essential for discrimination analysis.
Obviously, the more bands scored and plants studied, the
Nuclear and chloroplast DNA variability and phylogeny
123
higher the statistical significance of the calculation. About
100 bands should be enough to obtain statistically signifi-
cant results (Kocsis et al. 2005), however, Koller et al.
(1993) described 11 apple cultivars using only two RAPD
primers. Tancred et al. (1994) also used only two RAPD
primers to distinguish a new apple genotype from three
commercial apple cultivars. Autio et al. (1998) identified 15
apple rootstocks using only two RAPD primers, and Or-
aguzie et al. (2001) described 155 new and old apple
genotypes using only nine RAPD primers. In the present
study, one-hundred 10-mer primers, corresponding to
TIBMOLBIOL Co., Germany, were initially screened using
three genotypes from different regions to determine the
suitability of each primer for the study. Primers were
selected for further analysis based on their ability to detect
distinct, clearly resolved and polymorphic amplified pro-
ducts between genotypes. To ensure reproducibility, the
primers generating no, weak, or complex patterns were
discarded. In total, 11 primers out of 100 primers produced
clear and reproducible bands. Table 2 shows codes and
sequences of the used primers, total number of the ampli-
fication fragments (bands) resultant from apple samples and
the number of polymorphic fragments for each primer. A
total of 113 RAPD fragments were detected for the studied
genotypes with an average 10.27 bands per primer. Eleven
primers generated 77 polymorphic fragments with an
average 7 bands per primer, thus reflecting the high genetic
diversity within the germplasm, which distinguished
Granny-Smith cultivar from all the Iranian genotypes. Or-
aguzie et al. (2001) reported 42 polymorphic bands for 155
apple genotypes by nine RAPD primers and suggested high
variation and Bayan et al. (2007) reported 82 % polymor-
phism for Syrian local apples using RAPD, while Goulao
et al. (2001) stated that RAPD techniques revealed low level
of polymorphism percentage among apple cultivars.
The amplified fragments per primer ranged from seven
(for TIBMBD-16 primer) to 15 (for TIBMBC-20 primer)
(Table 2). Also, the polymorphic bands varied from four
(for TIBMBD-16 primer) to 10 (for TIBMBC-20 primer).
The relative number of polymorphic fragments to the total
number of amplified fragments ranged between 50 % for
TIBMBB-05 and TIBMBB-17, while TIBMBE-08 primer
showed the highest polymorphism (100 %) that agreed
with the finding of Khadivi-Khub et al. (2008) about this
primer in the study of sweet cherry. In the present study,
the average of polymorphism was 68.14 % for all primers.
Erturk and Akcay (2010) reported 180 polymorphic bands
in 11 apple genotypes by 38 RAPD primers with an
average of 40 % polymorphism. The PIC values for the
used primers were high in this study. The size of amplified
fragments ranged from 300 to 3,000 bp (Fig. 1) which
agreed with the findings of other apple studies by RAPD
(Koller et al. 1993; Mulcahy et al. 1993; Duneman et al.
1994; Autio et al. 1998; Luis et al. 2001). This polymor-
phism showed high genetic diversity in studied genotypes.
Also, previously, these set RAPD primers showed high
polymorphism and diversity in Prunus avium cultivars
(Khadivi-Khub et al. 2008), confirming the suitability of
these RAPD primers for analysis of diversification in fruit
trees that belong to Rosaceae family. In the present study,
studied genotypes were collected from different parts of
Iran, which might be the reasons for the higher percentage
polymorphism (68.14 %) in comparison with the results of
Goulao et al. (2001) about Malus domestica. It was pre-
viously shown that RAPD markers have great potential to
identify relationships among Malus domestica genotypes
(Koller et al. 1993; Mulcahy et al. 1993; Harada et al.
1993; Gardiner et al. 1996; Autio et al. 1998; Luis et al.
2001; Zhou and Li 2000). Some polymorphic bands pro-
duced by RAPD primers seemed to be unique for some
Table 2 List of RAPD primers,
their numbers of total and
polymorphic fragments and
percentage of polymorphism
used in this study
Primer Sequence 50 to 30 Total of
bands (a)
Polymorphic
bands (b)
Polymorphism
%(b/a 9 100)
TIBMBA-06 GGACGACCGT 8 5 62.50
TIBMBA-08 CCACAGCCGA 10 7 70.00
TIBMBA-20 GAGCGCTACC 12 9 75.00
TIBMBB-05 GGGCCGAACA 10 5 50.00
TIBMBB-17 ACACCGTGCC 14 7 50.00
TIBMBC-04 CCACGTGCCA 11 9 81.82
TIBMBC-12 CCTCCACCAG 9 6 66.67
TIBMBC-20 AGCACTGGGG 15 10 66.67
TIBMBD-16 GAACTCCCAG 7 4 57.14
TIBMBD-17 GTTCGCTCCC 8 6 75.00
TIBMBE-08 GGGAAGCGTC 9 9 100.00
Total 113 77
Mean 10.27 7.00 68.14
A. Khadivi-Khub et al.
123
genotypes, indicated that nuclear genome of these geno-
types had some special differences due to mutations (Matus
and Hayes 2002). For example, one band from primer
TIBMBA-08 was unique for Zodras-Oghlid genotype. It
should be noted that this genotype is the most early-rip-
ening fruit between studied genotypes, and may be this
unique band have been related to this traits.
Genetic similarity values between studied genotypes
varied from 0.38 to 0.72 (Table 3), indicating high intra-
specific variation between Iranian apple genotypes. This is
in agreement with a formerly study on Iranian apple
genotypes by Gharghani et al. (2009) and Farrokhi et al.
(2011). They studied genotypes of other regions of Iran and
reported high genetic variation for them, so that Farrokhi
et al. (2011) reported 0.19-0.79 genetic similarity for Ira-
nian apples. Bayan et al. (2007) reported 0.43 to 0.76
diversity between studied Syrian local apples using RAPD,
while, Adebayo et al. (2009) reported 65-87 % genetic
similarity between Nigerian apples by RAPD and sug-
gested a narrow genetic base in the local Nigerian selection
of apple.
The highest genetic similarity was scored between two
genotypes from Central Alborz that included Golab-Da-
mavand and Golab-Paize (72 %) and followed by simi-
larity between two genotypes from Azerbaijan, included
Mahali-Urmia and Azerbaijan6 (71 %), while the lowest
genetic similarity was scored between Zodras-Oghlid from
Central Zagros and Granny-Smith as an outgroup (38 %).
Furthermore, the lowest genetic similarity between Iranian
genotypes was observed between ShahreKord8, from
Central Zagros, and Zodras-Oghlid, from Central Alborz
(46 %) and followed by similarity between ME 4, from
Azerbaijan, and Golab-Arak, from Central Zagros (49 %).
The UPGMA cluster analysis revealed genetic rela-
tionships among genotypes. The cophenetic correlation
coefficient (CCC) indicated high correlation (r = 0.93)
between the similarity matrix and the cophenetic matrix
(obtained from the UPGMA dendrogram), indicating a
good representation of the molecular relationships among
genotypes. The cophenetic correlation coefficient is con-
sidered to be a very good representative of the data matrix
in the dendrogram if it is 0.90 or greater (Romesburg
1990). When the dendrogram determining the rate of rel-
ativity between genotypes on the basis of coefficient of
similarity was examined (Fig. 2), it was observed that the
diversity ratio among the apple genotypes was high. The
dendrogram generated by RAPD analysis showed two
distinct clusters with strong bootstrapping value (100 %)
for these two clusters. The first main cluster included all
Iranian genotypes except Torsh-Kermanshah, while,
Granny-Smith, foreign cultivar, was placed in the second
main cluster. The first cluster was divided into five subcl-
usters but was supported by low to moderate bootstrap
values. This clustering was according to geographical
distribution so that genotypes collected from Central Za-
gros, Central Alborz and Azerbaijan showed high variation
based on their geographical region. These results agreed
with the finding of other studies about apple genotypes
(Royo and Itoiz 2004; Adebayo et al. 2009) that reported
correlation between geographical distribution and intra-
specific variation in apple. Some genetic similarities were
observed between genotypes of different regions. For
example, Azerbaijan14 genotype, from Azerbaijan region,
showed high similarity with Golab-Paize genotype from
Central Alborz. Beside, Ferdos genotype from Central
Alborz was placed with Azerbaijanian genotypes, indicat-
ing possibility of genotype migration or gene flow between
regions. This finding confirmed that the RAPD method can
solve one of the major problems of varietal identification in
pome fruits: the existence of homonyms and synonyms,
particularly with regard to cultivars that were cultivated for
centuries and are widely distributed (Borrego et al. 2002).
Previously, RAPD analysis was used successfully to
reveal genetic variations in many studies (Kocsis et al.
2005; Lisek et al. 2006; Koc et al. 2009). Notwithstanding
the limitations, RAPD markers proved to be a highly
Fig. 1 RAPD profile of 23 Iranian apple genotypes and Granny-Smith cultivar as an outgroup, produced by TIBMBE-08 primer
Nuclear and chloroplast DNA variability and phylogeny
123
Ta
ble
3G
enet
icsi
mil
arit
ym
atri
xb
ased
on
RA
PD
dat
ab
etw
een
stu
die
dap
ple
gen
oty
pes
esti
mat
edac
cord
ing
toJa
ccar
d’s
met
ho
d
T3
Aze
r7S
BA
ME
4M
ahal
iF
erdos
Bey
ghi
ME
3S
hah
l0S
hah
l9G
har
aK
han
GolD
aG
olP
aG
olA
rA
zerl
4M
esh
ShK
ord
8G
oD
Sfe
TorK
erG
her
m1
Aze
r6Z
odO
gh
Gra
nny
T3
1.0
0
Aze
r70.6
81.0
0
SB
A0.5
60.5
11.0
0
ME
40.6
90.6
00.5
81.0
0
Mah
ali
0.6
50.6
80.5
40.6
21.0
0
Fer
dos
0.5
90.6
40.4
80.5
50.6
11.0
0
Bey
ghi
0.4
90.6
20.5
00.6
60.6
40.7
01.0
0
ME
30.5
60.6
50.4
70.5
30.6
20.6
00.6
11.0
0
Shah
lO0.6
00.6
10.4
90.6
40.6
00.5
40.5
60.6
21.0
0
Shah
l90.6
60.6
20.5
50.5
10.6
10.5
70.5
30.5
50.7
11.0
0
Ghar
a0.6
40.5
70.5
50.5
10.6
40.6
00.5
50.5
50.6
20.6
31.0
0
Khan
0.6
70.6
10.5
70.5
70.6
00.6
40.5
40.5
90.6
50.5
90.7
21.0
0
GolD
a0.6
80.5
80.5
40.5
70.6
30.5
90.4
40.5
40.6
80.6
80.6
20.6
31.0
0
GolP
a0.6
00.6
10.5
90.5
70.5
80.5
90.5
20.5
20.5
80.6
60.5
40.6
00.7
21.0
0
GolA
r0.5
50.5
30.5
10.4
90.4
80.6
10.6
20.5
40.5
50.5
60.5
90.6
30.6
00.6
31.0
0
Aze
rl4
0.5
90.5
50.6
50.6
10.5
90.5
60.6
00.5
80.6
60.6
90.5
80.6
60.5
90.6
60.5
91.0
0
Mes
h0.5
80.6
60.5
10.5
40.5
80.6
80.5
60.7
20.6
00.5
40.6
70.6
50.5
50.5
10.6
30.6
81.0
0
ShK
ord
80.6
10.6
10.4
90.5
70.5
50.5
40.5
90.6
80.5
50.4
90.5
70.6
60.5
50.4
60.5
50.6
20.6
31.0
0
GoD
Sfe
0.6
20.5
40.5
50.5
60.5
90.5
30.4
50.5
80.6
20.6
50.5
80.6
20.7
20.6
90.5
10.6
30.5
60.5
71.0
0
TorK
er0.4
90.5
20.6
60.5
60.5
20.5
10.5
30.5
30.5
70.5
10.5
30.5
70.5
40.5
20.5
20.6
50.6
20.4
70.5
11.0
0
Gher
ml
0.5
20.5
20.6
10.5
40.5
90.5
50.6
10.5
30.5
70.6
50.5
80.6
20.5
40.5
20.6
20.6
10.5
40.6
00.6
60.5
41.0
0
Aze
r60.6
50.7
10.6
40.6
20.7
20.5
90.6
40.6
40.6
00.5
90.6
40.6
50.6
70.6
30.6
30.6
00.5
80.6
10.6
00.5
50.6
21.0
0
ZodO
gh
0.5
80.5
90.7
00.5
30.5
80.5
70.5
20.5
20.5
10.4
90.6
50.6
60.6
60.6
60.5
80.6
00.6
10.5
60.6
00.5
30.5
00.6
61.0
0
Gra
nny
0.5
30.5
60.4
90.6
30.5
90.4
90.5
50.5
50.4
50.5
20.4
10.4
90.4
80.5
10.5
10.5
30.4
30.5
40.4
90.5
30.5
30.6
20.3
81.0
0
A. Khadivi-Khub et al.
123
effective and efficient method for the genetic analysis (Luo
et al. 2002; Ulanovsky et al. 2002). Large numbers of data
sets can be generated because different RAPD primers are
commercially available (Fanizza et al. 2000). Our results
confirmed that RAPD analysis can be used to differentiate
apple genotypes and confirmed the potential of RAPD
technology as a reliable, rapid and inexpensive screening
method to discriminate quince genotypes.
Chloroplast DNA analysis
In the present study, from the four universal chloroplast
primer pairs tested on apple studied genotypes (K1K2, CS,
HK and TF), only primer pair TF that amplified fragment
trnT- trnF in tobacco, did not amplify any band in studied
apple genotypes. Similar result was reported by Mohanty
et al. (2001) that used TF primer pairs on Prunus avium
having the same family with genus Malus and did not
obtain any amplification. Other primer pairs amplified the
fragments (Table 4). Chloroplast DNA size in apple has
been reported by Ishikawa et al. (1992) to digest this
molecule by restriction enzyme (BamHI) in Golden
Delicious, McIntosh and Red Delicious cultivars which
are, respectively, 143.4, 143.35 and 143.7 kb. We con-
sidered this molecule (143.48 kb) by averaging the
reported sizes. Therefore, the fragment that primers
amplified (5920 bp) by approximately 4.13 % of the
chloroplast genome was analyzed.
Cluster I
Shahrod19
Azerbaijan14
Golab-Arak
Golab-Paize24
46
Golab-Damavand
SBA48
40
Azerbaijan6
T3
Khan-Almasi
Zodras-Oghlid26
14
8
Shahrod10
Ghermez1
Ghara-Yarpagh
Golab-Nemati11
6
ME3
Meshkinshahr30
4
3
1
ME4
ShahreKord818
1
Azerbaijan7
Ferdos22
Beyghi
Mahali-Urmia34
15
51
Granny-Smith
Torsh-Kermanshah41
100
Cluster II
Fig. 2 Dendrogram of genetic
relationships between apple
genotypes, constructed by
UPGMA based on 11 RAPD
markers. Numbers in the
branches represent bootstrap
values
Nuclear and chloroplast DNA variability and phylogeny
123
From six combinations of primer-restriction enzyme
used in this study for analysis of cpDNA, two combina-
tions, K1K2-MseI and CS-MseI, resolved monomorphic
fragments, and the region amplified by CS primer did not
have any site for EcoRI. The remaining combinations
showed polymorphic pattern. Combinations of K1K2-
EcoRI and HK-MseI had two sites for enzymes. HK-MseI
showed two polymorphic patterns but in combination with
K1K2-EcoRI having two sites; one of the fragments was
monomorphic and the other one was polymorphic
(Table 4). All of these polymorphic patterns resulted
because of insertion–deletion (indel) mutation between 20
and 50 bp ranges. Indel mutations in the combinations of
K1K2-EcoR I and HK- EcoRI are shown in Fig. 3.
This estimation illustrated that the genetic similarity
was varied. In some cases, studied genotypes showed
100 % similarity that confirmed the conservation of the
chloroplast genome within these genotypes (Palmer 1987).
According to polymorphism pattern, all studied genotypes
were placed into two main clusters with strong bootstrap-
ping value (100 %) (Fig. 4). Based on this clustering, in
some case, intraspecific chloroplast variation was observed
within genotypes. This could be due to the existence of
intraspecific cpDNA mutation, which has been previously
reported in sweet cherry (Mohanty et al. 2001) and apple
(Coart et al. 2006). Also, in some cases, genotypes of
different regions showed high chloroplast similarity. The
first main cluster was divided to two subclusters with low
to moderate bootstrapping values that most of Azerbaja-
nian genotypes were placed in this cluster. Also, the sec-
ond cluster was divided to two subclusters that Granny-
Smith cultivar showed high chloroplast similarity to Torsh-
Kermanshah genotype. These two apples have sour fruit
flavor.
Chloroplast DNA is known to be a very conservative
molecule (Palmer 1987), hence chances of detecting
intraspecific cpDNA polymorphism are low. But, in the
present study, intraspecific polymorphism was observed in
most cases. This result agreed with the finding of Coart
et al. (2006) that reported intraspecific chloroplast variation
in apple. But, Verbylait _e et al. (2006) evaluated chloroplast
of Maloideae and reported no intraspecific chloroplast
variation within studied species because of low studied
genotypes for each species.
The results of the present study demonstrated that,
although cpDNA is a very conservative cytoplasmic
molecule, variation exists therein to distinguish even at
the level of intraspecific cultivars. Intraspecific diversity
in plant cpDNA has been believed to be the result of
mutations, genetic drift, selection by man for particular
cpDNA-encoded characters (herbicide resistance), or the
Fig. 3 Restriction patterns in two chloroplast DNA primer-enzyme combinations HK-EcoRI (a) and K1K2-EcoRI (b) resolved on agarose gel
for chloroplast analysis
Table 4 Nei’s analysis of gene diversity in studied apple genotypes
based on chloroplast DNA
cpDNA
primers
Combination Digestion
pattern
Ht Hs Gst I
KIK2 K1K2-EcoRI KEA 0.50 0.49 0.03 0.68
CS K1K2-EcoRI KEB 0.50 0.49 0.03 0.68
HK HK- EcoRI HEA 0.39 0.30 0.23 0.67
TF HK- EcoRI HEB 0.39 0.30 0.23 0.67
HK-MseI HM1A 0.46 0.33 0.29 0.69
HK-MseI HM1B 0.46 0.33 0.29 0.69
HK-MseI HM2A 0.49 0.49 0.01 0.67
HK-MseI HM2B 0.34 0.29 0.15 0.61
HK-MseI HM2C 0.45 0.43 0.05 0.61
Mean – 0.44 0.38 0.14 0.67
St. Dev – 0.01 0.01 – 0.03
Ht: Total diversity; Hs: Average intraspecific diversity; GSTC: level
of population subdivision; I: Shannon’s index
A. Khadivi-Khub et al.
123
introgression of a foreign cpDNA (chloroplast capture)
from related species during the early domestication of the
taxon and even in wild taxa (Harris and Ingram 1991).
Such introgression has been reported to be common for
many plant species (e.g. in Iris species; Arnold et al.
1992).
In some cases, genotypes of Central of Alborz region,
which is the nearest region to Azerbaijan, were beside
Azerbaijanian genotypes. This means that these genotypes
probably have been transferred from Azerbaijan that has a
vast area of culture and history of apple. Meanwhile, the
same maternal parent cannot be completely rejected
because introduction from other areas and selection for
better genotypes have been factors for transmission of
genotypes between regions from the past. One of the main
reasons for the absence of genetic structure among
Golab-Arak
Khan-Almasi
Ghara-Yarpagh
T330
77
Ghermez1
Golab-Paize83
40
ShahreKord8
Azerbaijan14
Ferdos24
88
17
47
Meshkinshahr
Azerbaijan6
Azerbaijan792
Shahrod10
Shahrod1990
65
32
21
Beyghi
ME390
Golab-Damavand
Mahali-Urmia85
48
Golab-Nemati
ME441
30
Granny-Smith
Torsh-Kermanshah53
SBA
Zodras-Oghlid98
13
20
100
Cluster I
Cluster II
Fig. 4 UPGMA dendrogram of studied apple genotypes based on chloroplast data. Numbers in the branches represent bootstrap values
Nuclear and chloroplast DNA variability and phylogeny
123
different populations appears to be cytoplasmic gene flow
among them (that is, seeds migration). The migration rate
of seeds may be higher than expected (Powell et al. 1996).
Phylogenetic analysis revealed the evolution of studied
germplasm (not shown), so that Beyghi, ME3 and Mahali-
Urmia genotypes, from Azerbaijan, had the most close to
their ancestors, while, Ferdos and ShahreKord8 genotypes
showed the lowest close to their ancestors. Granny-Smith
cultivar and Torsh-Kermanshah genotype showed similar
phylogeny. It should be noted that Torsh-Kermanshah has
sour flavor like Granny-Smith.
Mean of total diversity (Ht), average intraspecific diver-
sity (Hs) and level of population subdivision (GSTC) calcu-
lated by POPGENE program were 0.44, 0.38 and 0.14,
respectively (Table 4). According to the genetic differenti-
ation among genotypes of different regions in apple
(GSTC = 0.14) and comparison of it by similar reports such
as in Prunus avium (GSTC = 0.08; Mohanty et al. 2001) and
Prunus spinosa (GSTC = 0.19; Mohanty et al. 2000) which
are from the same family with the genus Malus (Rosaceae), it
could be determined that chloroplast genetic differentiation
among regions in the present study (GSTC) is low and is much
Ferdos
T3
Ghara-Yarpagh
Khan-Almasi54
34
Azerbaijan14
ShahreKord817
4
3
Ghermez1
Golab-Arak27
1
ME3
Meshkinshahr25
Mahali-Urmia
Azerbaijan6
Azerbaijan764
22
Shahrod10
Shahrod1981
2
Golab-Paize
Golab-Damavand
Golab-Nemati47
25
1
1
Torsh-Kermanshah
SBA
Zodras-Oghlid82
16
3
Granny-Smith
Beyghi
ME424
6
100
Cluster I
Cluster II
Fig. 5 UPGMA dendrogram of
studied apple genotypes based
on integrating data of nuclear
and chloroplast DNA
(RAPD ? cpRFLP). Numbers
in the branches represent
bootstrap values
A. Khadivi-Khub et al.
123
lower compared to forest species such as Quercus petraea
(GSTC = 0.90; Petit et al. 1993a), Fagus sylvatica
(GSTC = 0.83; Demesure et al. 1996), Argania spinosa
(GSTC = 0.60; El Mousadik and Petit 1996) and Alnus glu-
tinosa (GSTC = 0.87; King and Ferris 1998). Petit et al.
(1993b) reported that in some species the GSTN (genetic
differentiation among populations using nuclear markers) is
much lower than the GSTC. It could be deduced that cyto-
plasmic genomes which have inheritance from maternal
parents are relatively much more structured than the nuclear
genome because cytoplasmic gene flow is limited just by
seed dispersal (Powell et al. 1996).
Combined analyses of RAPD and cpDNA
The collection and characterization of local germplasm are
of primary importance for the prevention of their loss and
to allow the exploitation of materials that are potentially
useful for future purposes (Palmer 1987). Furthermore, the
characterization of collections of local germplasm can also
help in the clarification of some aspects concerning the
diffusion and evolution of crop plants. Therefore, a char-
acterization of local genotypes and cultivars of apple from
different Iran regions can help us to understand the path-
ways of dissemination and evolution of this species.
0.1
Azerbaijan6
Beyghi
ME422
Granny-Smith
Torsh-Kermanshah27
4
Golab-Damavand
Mahali-Urmia10
Golab-Nemati
Shahrod10
Shahrod1985
3
Ferdos
Azerbaijan7
Meshkinshahr7
2
ME3
ShahreKord814
Golab-Arak
Ghermez1
Khan-Almasi2
1
Ghara-Yarpagh
T314
1
Azerbaijan14
Golab-Paize9
SBA
Zodras-Oghlid81
4
100
Cluster I
Cluster II
Fig. 6 NJ phylogenetic
dendrogram of studied apple
genotypes based on integrating
data of nuclear and chloroplast
DNA (RAPD ? cpRFLP).
Numbers in the branches
represent bootstrap values
Nuclear and chloroplast DNA variability and phylogeny
123
Accordingly, in the present study, 23 apple genotypes were
analyzed using two different types of markers, including
nuclear (RAPD) and chloroplast markers (PCR–RFLP).
The combined use of different types of markers from two
genomes that are likely to be unevenly affected by major
evolutionary forces, such as migration, selection, mutation,
genetic drift and recombination, offers the opportunity to
make inferences from the patterns of genetic diversity, and
thus to define a detailed evolutionary scenario and to fully
describe the genetic makeup of different groups of germ-
plasm (Provan et al. 2001). Our results showed that this
approach is very useful and informative. For instance, the
use of markers showing biparental vs uniparental trans-
mission allowed us to detect the structure of genetic
diversity in the apple. Similarly, we were able to observe
phylogenetic diversity in studied apples based on the
cytoplasmic (chloroplast). Overall, our results argue for a
co-ordinated multiple marker approach for germplasm
characterization, conservation and evaluation. Based on
data, studied apple genotypes showed a higher nuclear
variation (RAPD) than chloroplast variation.
The level of polymorphism revealed by nuclear (RAPD)
and chloroplast markers could determine the relationships
between the DNA profiles and the geographical distribution
of some genotypes. Chloroplast markers are more suitable
for this purpose due to the low rate of mutations and
recombination in the chloroplast genome (Provan et al.
2001) as well as because of their maternal inheritance in
Malus. Also, some cases of similarity were observed
between data based on RAPD and chloroplast data. For
instance, genotypes of Central Alborz region showed
similarity with Azerbaijanian genotypes according to these
two data. The estimation of correlation between matrices of
nuclear data (RAPD) and cpDNA data (r = 0.40) indicated
non-significant correlation between these data and agreed
with the finding of Horvath et al. (2008) that reported an
insignificant correlation coefficient between nuclear DNA
and chloroplast DNA data in Prunus cerasifera (cherry
plum).
Three dendrograms were constructed based on the
combined data of both markers (RAPD ? cpRFLP) using
UPGMA, NJ and Dollo Parsimony methods. The
Fig. 7 Dollo Parsimony phylogenetic dendrogram of studied apple genotypes based on integrating data of nuclear and chloroplast DNA
(RAPD ? cpRFLP). Numbers in the branches represent bootstrap values
A. Khadivi-Khub et al.
123
advantage of combining different datasets is a compre-
hensive taxonomic picture, since each represents a distinct
level of taxonomic differentiation and confirmed view-
point of Sneller et al. (1997) in this case. The UPGMA
dendrogram produced based on combined data presented a
high resolution in separating the different apple genotypes
into two major clusters but with a high bootstrapping
value (100 %) (Fig. 5). The first cluster consists of three
subclusters. Subcluster I contained most genotypes from
Central Alborz. Subcluster II contained two subgroups, so
that subgroup I contained Azerbaijanian genotypes, sug-
gesting a possible genetic relatedness between them.
Another interesting feature was that Golab cultivars which
had much similarity in tree, flower, fruit, flowering and
ripening were placed in the same subgroup by the com-
bined analysis. Subcluster III was consisted of genotypes
of Central Zagros and showed high molecular genetic
similarity between them. Finally, the second major cluster
was consisted of two genotypes of Azerbaijan and
Granny-Smith. These two genotypes have green fruit skin
like Granny-Smith.
Phylogenetic dendrogram according to the integrating
data based on both data (RAPD ? cpRFLP) using NJ
method revealed the evolution of studied germplasm but
with low and moderate bootstrapping values (Fig. 6),
while, phylograms derived from Dollo parsimony analysis
performed using PAUP were similar to those produced by
NJ method in most cases with moderate to high boot-
strapping values (50–100 %) (Fig. 7), so that most of the
bootstrapping values were 50–80 % and some of them
were strong (80–100 %). In both dendrograms of NJ and
Dollo parsimony methods, most cases of similarities were
observed. For instance, genotypes Shahrod10, Shahrod19,
Azerbaijan7 and Azerbaijan6 were placed in the same
group. Also, T3, Ghara-Yarpagh and Khan-Almasi showed
close similarity based on both clusters. Furthermore, both
dendrograms detected close relationships between Mahali-
Urmia, ME3, Beyghi, ME4 and Granny-Smith and also
among ShahreKord8, Ferdos, Meshkinshahr and Azerbai-
jan14. Beside, these similarities were observed for other
genotypes. For the 24 apple genotypes analyzed in the
present study, there existed a theoretical expectation for the
clustering of genotypes within the dendrogram, based on
their botanical pedigree. Genetic similarity dendrograms
were constructed either using the results taken by RAPD or
cpRFLP analysis and by their combination. The prevailing
opinion over the past few years is that nuclear and chlo-
roplast genome-based approaches are complementary
methods to provide the most complete coverage of these
genomes. Hence, in the pedigree analysis, the combined
RAPD ? cpRFLP analysis was found to be the most
accurate system for the identification of genotypes from the
different regions.
Conclusion
We have demonstrated that old varieties are potential
sources of diversity in apples, in addition to the wild
species and current cultivars. Different apple cultivars and
local genotypes have been cultivated in Iran for a long
time, possibly for many centuries and are least frequently
grown in the southeastern, northeastern and central
regions and most widespread in the north and west of Iran.
It can be generalized that genotypes selected from these
regions most likely represent the majority of genotypes
cultivated or local varieties in Iran. The Azerbaijanian
genotypes did not completely separate from other samples,
and this could be due to migration of genotype by mostly
human activity or gene flow. The genotypes from this
collection will be used in the forthcoming cross-breeding
studies. The characterization results will be of great
importance on the selection of the parents and the com-
binations. Information on the origin of cultivated and its
phylogenetic relationship with the closely related species
is now becoming more and more important for future
apple breeding.
References
Adebayo OL, Bola O, Opeyemi W, Gloria M, Temitope O (2009)
Phylogenetic and genomic relationships in the genus Malus
based on RAPDs. Afri J Biotech 815:3387–3391
Arnold ML, Robinson JJ, Buckner CM, Bennet BD (1992) Pollen
dispersal and interspecific gene flow in Louisiana irises.
Heredity 68:399–404
Autio WR, Schupp JR, Ferree DC, Glavin R, Mulcahy DL (1998)
Application of RAPDs to DNA extracted from apple rootstocks.
Hort Science 33:333–335
Bayan MM, Rania AAY, El-Halabi O, Omayma MI (2007) Genetic
identification of some Syrian local apple (Malus sp.) cultivars
using molecular markers. Research Journal of Agriculture and
Biological Sciences 3:704–713
Borrego J, Andres MT, Gomez JL, Ibanez J (2002) Genetic study of
Malvasia and Torrontes groups through molecular markers.
American Journal of Enology and Viticulture. 53:125–130
Coart E, van Glabeke S, de Loose M, Larsen AS, Roldan-Ruiz I
(2006) Chloroplast diversity in the genus Malus: newinsights
into the relationship between the European wild apple (Malus
sylvestris (L.) Mill.) and the domesticated apple (Malus domes-
tica Borkh.). Mol Ecol 15:2171–2182
Demesure B, Comps B, Petit RJ (1996) Chloroplast DNA phyloge-
ography of the common beech (Fagus sylvatica L.) in Europe.
Evol 50:2515–2520
Doyle JJ, Doyle JL (1987) A rapid DNA isolation procedure for small
quantities of fresh leaf tissue. Phytochem. Bull. 19:11–15
Dumolin-Lapegue S, Pemonge MH, Petit RJ (1997) An enlarged set
of consensus primers for the study of organelle DNA in plants.
Mol Ecol 6:393–397
Duneman F, Kahnau R, Schmidt H (1994a) Genetic relationships in
Malus evaluated by RAPD fingerprinting of cultivars and wild
species. Plant Breeding 113:150–159
Nuclear and chloroplast DNA variability and phylogeny
123
Duneman F, Kahnau R, Schmidt H (1994b) Genetic relationships in
Malus evaluated by RAPD fingerprinting of cultivars and wild
species. Plant Breed 113:150–159
El Mousadik A, Petit RJ (1996) Chloroplast DNA phylogeography of
the argan tree of Morocco. Mol Ecol 5:547–555
Erturk U, Akcay MA (2010) Genetic variability in accessions of
‘Amasya’ apple cultivar using RAPD markers. Not Bot Hort
Agro Cluj-Napoca 38:239–245
Fanizza G, Corona MG, Resta P (2000) Analysis of genetic
relationships among Muscat grapevines in Apulia (South Italy)
by RAPD markers. Vitis 39:159–161
Farrokhi J, Darvishzadeh R, Naseri L, MohseniAzar M, Hat-
amiMaleki H (2011) Evaluation of genetic diversity among
Iranian apple (Malus 9 domestica Borkh.) cultivars and land-
races using simple sequence repeat markers. Aust J Crop Sci
5:815–821
Gardiner SE, Bassett HCM, Madie C, Noition DAM (1996) Isozyme,
randomly amplified polymorphic DNA (RAPD), restriction
fragmentlength polymorphism (RFLP) markers to deduce a
putative parent for the ‘Braeburn’ apple. J Amer Soc Hort Sci
121:996–1001
Gharghani A, Zamani Z, Talaie A, Oraguzie NC, Fatahi R, Hajnajari
H, Wiedow C, Gardiner SE (2009) Genetic identity and
relationships of Iranian apples (Malus 9 domestica Borkh.)
cultivars and landraces, wild apple species and representative old
apple cultivars based on SSR markers. Genet Res Crop Evol
56:829–842
Goulao L, Cabrita L, Oliveira CM, Leitao JM (2001) Comparing
RAPD and AFLP analysis in discrimination and estimation of
genetic similarities among apple (Malus domestica Borkh.)
cultivars. Euphytica 119:259–270
Grivet D, Heinze B, Vendramina GG, Petit RJ (2001) Genome
walking with consensus primers: Application to the large single
copy region of chloroplast DNA. Mol Ecol 1:345–349
Harada T, Maksukawa K, Sato T, Ishikawa Niizeki R, Saito KM
(1993) DNA-RAPD detect genetic variation and paternity in
Malus. Euphytica 65:87–91
Harris SA, Ingram R (1991) Chloroplast DNA and biosystematics -
the effects of intraspecific diversity and plastid transmission.
Taxon 40:393–412
Ishikawa S, Kato S, Imakawa S, Mikami T, Shimamoto Y (1992)
Organel DNA polymorphism in apple cultivars and rootstocks.
Theor Appl Genet 83:963–967
Jackson JE (2003) Biology of Apples and Pears. Cambridge
University Press, Cambridge UK 959
Janick J, Cummins JN, Brown SK, Hemmat M (1996) Apples. In:
Janick J, Moore JN (eds) Fruit Breeding, vol I., Tree and
Tropical Fruits. John Wiley and Sons, New York pp, pp 1–77
Khadivi-Khub A, Zamani Z, Bouzari N (2008) Evaluation of genetic
diversity in some Iranian and foreign sweet cherry cultivars by
using RAPD molecular markers and morphological traits. Hortic
Environ Biotechnol 49:188–196
King RA, Ferris C (1998) Chloroplast DNA phylogeography of Alnus
glutinosa (L.). Mol Ecol 7:1151–1161
Kishimoto S, Aida R, Shibata M (2003) Identification of chloroplast
DNA variations by PCR-RFLP analysis in Dendranthema. J Jap
Soc Hort Sci 72:197–204
Koc A, Akbulut M, Orhan E, Celik Z, Bilgener S, Ercisli S (2009)
Identification of Turkish and standard apple rootstocks by
morphological and molecular markers. Genet and Mol Res
8:420–425
Kocsis M, Jaromi L, Putnoky P, Kozma P, Borhidi A (2005) Genetic
diversity among twelve grape cultivars indigenous to the
Carpathian Basin revealed by RAPD markers. Vitis 44:87–91
Koller BA, Lehmann J, Modermott M, Gessier C (1993) Identification
of apple cultivars using RAPD markers. Theor Appl Genet
85:6–7
Landry BS, Li RQ, Cheung WY, Granger RL (1994) Phylogeny
analysis of 25 apple rootstocks using RAPD markers and tactical
gene tagging. Theor Appl Genet 89:847–852
Lisek A, Korbin M, Rozpara E (2006) Simple identification of sweet
cherry rootstocks PHL using RAPD markers. Sod Dar
25:316–321
Luis G, Cabrita L, Oliveira CM, Leitao JM (2001) Comparing RAPD
and AFLP analysis in discrimination and estimation of genetic
similarities among apple (Malus domestica Borkh.) cultivars.
Euphytica 119:250–270
Luo S, He P, Zheng X, Zhou P (2002) Inheritance of RAPD markers
in an interspecific F1 hybrid of grape between Vitis quinquang-
ularis and V. vinifera. Sci Hort 93:19–28
Matus IA, Hayes PM (2002) Genetic diversity in three groups of
barley germplasm assessed by simple sequence repeats. Genome
45:1095–1106
Mohanty A, Martın JP, Aguinagalde I (2000) Chloroplast DNA
diversity within and among populations of the allotetraploid
Prunus spinosa L. Theor Appl Genet 100:1304–1310
Mohanty A, Martin JP, Aguinagaldo I (2001) Chloroplast DNA study
in wild population and some cultivars of Prunus avium. Theor
Appl Genet 103:112–117
Mulcahy DL, Cresti M, Sansavini S, Douglas GC, Linskens HF,
Mulcahy GB, Vignani R, Pancaldi M (1993) The use of random
amplified polymorphic DNAs to fingerprint apple genomes. Sci
Hort 54:89–96
Oraguzie NC, Gardiner SE, Basset HCM, Stefanati M, Ball RD, Bus
VGM, White AG (2001) Genetic diversity and relationship in
Malus sp. Germplasm collections as determined by Random
Amplified Polymorphic DNA. J Amer Soc Hort Sci 126:318–328
Page RDM (1996) TreeView: an application to display phylogenetic
trees on personal computers. Comp Appl Biosci 12:357–358
Palmer JD (1987) Chloroplast DNA evolution and biosystematics
uses of chloroplast DNA variation. Am Nat 130:S6–S29
Panda S, Martin JP, Aguinagalde I, Mohanty A (2008) Chloroplast
DNA variation in cultivated and wild Prunus avium L.: a
comparative study. Plant Breed 122:92–94
Pavlicek A, Hrda S, Flegr J (1999) FreeTree - Freeware program for
construction of phylogenetic trees on the basis of distance data
and bootstrap/jackknife analysis of the tree robustness. Folia
Biol (Praha) 45:97–99
Petit RJ, Kremer A, Wagner DB (1993a) Geographic structure of
chloroplast DNA polymorphisms in European oaks. Theor Appl
Genet 87:122–128
Petit RJ, Kremer A, Wagner DB (1993b) Finite island model for
organelle and nuclear genes in plants. Hered 71:630–641
Powell W, Morgante M, McDevitt R, Vendramin GC, Rafalski JA
(1996) Polymorphic simple sequence repeat regions in chloro-
plast genomes: applications to the population genetics of pines.
Proc Natl Acad Sci USA 92:7759–7763
Provan J, Powell W, Hollingsworth PM (2001) Chloroplast micro-
satellites: new tools for studies in plant ecology and evolution.
Trends Ecol Evol 16:142–147
Rohlf, F.J., 2000. NTSYS-pc Numerical Taxonomy and Multivariate
Analysis System. Version 2. 1. Exeter Software, Setauket, NY
Romesburg HC (1990) Cluster Analysis for Researchers. Krieger
Publishing, Malabar, FL USA
Royo JB, Itoiz R (2004) Evaluation of the discriminance capacity of
RAPD, isoenzymes and morphologic markers in apple (Malus x
domestica Borkh.) and the congruence among classifications.
Genet Res Crop Evol 51:153–160
A. Khadivi-Khub et al.
123
Sneller CH, Miles JW, Hoyt JM (1997) Agronomic performance of
soybean plant introductions and their genetic similarity to elite
lines. Crop Sci 37:1595–1600
Takeuchi S, Suzuki K, Uchiyama H, Yoneda K, Suzuki S (1995)
Identification and phylogeny of wild and cultivated roses using
DNA markers. 3. Maternal inheritance of chloroplast DNA and
inference of mother line in some hybrids. J Jap Soc Hort Sci
64:400–401
Tancred SJ, Zeppa AG, Graham GC (1994) The use of the PCR-
RAPD technique in improving the plant variety rights descrip-
tion of a new Queensland apple (Malus domestica) cultivar. Aust
J ExperAgri 34:665–667
Ulanovsky S, Gogorcena Y, Toda FM, Ortiz JM (2002) Use of
molecular markers in detection of synonymies and homonymies
in grapevines (Vitis vinifera L.). Sci Hort 92:241–254
Verbylaite R, Ford-Lloyd B, Newbury J (2006) The phylogeny of
woody Maloideae (Rosaceae) using chloroplast trnL-trnF
sequence data. Biologija 1:60–63
Welsh J, McClelland M (1990) Fingerprinting genomes using PCR
with arbitrary primers. Nucleic Acids Res 18:7213–7218
Williams JGK, Kubelik AR, Livak KJ, Rafalski JA, Tingey SV
(1990) DNA polymorphisms amplified by arbitrary primers are
useful as genetic markers. Nucleic Acids Res 18:6531–6535
Yeh FC, Yang RC, Boyle T (1999) POPGENE. Microsoft Windows-
based Freeware for Population Genetic Analysis. Release 1.31.
University of Alberta, Edmonton
Zhou ZQ, Li YN (2000) The RAPD evidence for the phylogenetic
relationship of the closely related species of cultivated apple.
Genet Res Crop Evol 47:353–357
Nuclear and chloroplast DNA variability and phylogeny
123