Application of High-Throughput Sequencing to Evaluate the Genetic Diversity Among Wild Apple Species Indigenous to Shandong,China, and Introduced Cultivars

The Taiyi mountainous region of Shandong province in eastern China has an abundance of wild Malus species. We evaluated the genetic diversity of 88 Malus accessions (45 Asian apple cultivars, 10 American apple cultivars, 12 European apple cultivars, 19 Chinese wild apples, and two apple cultivars with unknown origins) based on single-nucleotide polymorphism (SNP) markers. A total of 38,364 SNPs were obtained with an average of 2256 SNPs per chromosome. The average of the polymorphism information content (PIC), gene diversity, and allele frequency for SNPs was 0.268, 0.306, and 0.364, respectively. A circular phylogenetic tree constructed based on SNP data revealed that the 88 Malus accessions could be divided into three groups. However, a population structure analysis suggested the 88 Malus accessions could be divided into four groups. A principal component analysis (PCA) revealed some population stratification. The first three PCs accounted for 41.62% of the population-wide SNP variation, with PC1 accounting for 33.9%. Moreover, the kinship values of the 88 Malus accessions ranged from 0 to 2.36, with 96.42% of the kinship values between 0 and 0.2. A phylogenetic tree and a PCA indicated the Chinese wild apples widely distributed among the cultivated apples had a diverse genetic background. Characterizing the genetic relationships between cultivated apples and Chinese wild apples is essential for increasing the genetic diversity of the germplasms used by apple breeders.

     

To what extent do wild apples in Kazakhstan retain their genetic integrity?

Kazakhstan belongs to the center of origin of apple. Malus sieversii (Ledeb.) M. Roem., the ancestral progenitor of the cultivated apple is native to this region. Pressure on the natural habitats of this wild apple has been intensified due to agriculture, grazing, and urbanization in the last century. For decades, M. sieversii in Kazakhstan has been subjected to the “Red Book of the Kazakh SSR” and today, this species is threatened with extinction. Wild apple undergoes exceptional losses in habitats, and the risk for losing the genetic integrity becomes worse due to increasing cultivation of cultivated apples and frequently occurring crosspollination events. The present study was focused on the current state of M. sieversii in Kazakhstan, the level of its diversity, its genetic integrity, and the identification of regions where future activities for conservation will have a good chance of success. A total of 311 M. sieversii samples of 12 populations collected in the wild, 16 previously selected wild apple genotypes, and 50 grown cultivars were studied using 16 simple sequence repeat (SSR) markers for genetic analysis. The results suggest that the level of genetic diversity is high. The differentiation between the populations was low, although the within-population heterozygosity was relatively high. A significant number of hybrids (8–95%) between M. sieversii and cultivated apples were found suggesting frequent crop-to-wild gene flow. The percentage of pure wild apple genotypes was highest in Krutoe truct and Tauturgen. These sites should be taken into account for future in situ long-term preservation activities.     

Genetic Diversity and Identity of Chinese Loquat Cultivars/Accessions (<Emphasis Type="Italic">Eriobotrya japonica</Emphasis>) Using Apple SSR Markers

Loquat (Eriobotrya japonica) is an underutilized fruit crop that originated in China and for which only a small number of molecular markers are available. This number can be increased by identifying apple SSRs that are transferable to loquat cultivars/accessions to provide new insight into the level of genetic diversity within loquat and synteny with apple. We evaluated 71 apple SSR markers distributed across 17 linkage groups, and identified 39 SSRs transferable to loquat. Testing 54 loquat accessions, from Japan, Spain, four provinces in China, and two wild species gave a total of 155 different alleles with a mean value of 3.38 per locus. The mean effective number of alleles was 2.21, and the mean observed heterozygosity was 0.47. These values indicate a high degree of genetic diversity in the set of Chinese loquat accessions analyzed. Unweighted pair-group method analysis based on simple matching coefficent clustered the accessions into two groups, cultivated and wild loquat. The cultivated loquat can be subdivided into three subgroups which generally reflect their geographic origin in China. The Spanish cultivars clustered with those of the Jiangsu and Zhejiang provinces. A core set of five SSR markers could distinguish most accessions.     

Polymorphism of SSR alleles in pear cultivars grown in Belarus

Using SSR markers designed for Malus × domestica Borkh. genetic polymorphism of 43 pear accessions cultivated in Belarus was examined. A total of 217 alleles were identified with the mean number of 12.8 alleles per marker. The mean PIC value was 0.81; the mean number of informative alleles, 6.49. The heterozygosity level ranged from 0.30 to 0.84. Genetic diversity of SSR alleles in pear and apple genomes was compared. A method of identification of commercial pear cultivars using a set of six SSR markers was suggested.     

Unravelling genetic diversity and cultivar parentage in the Danish apple gene bank collection

Characterization of apple germplasm is important for conservation management and breeding strategies. A set of 448 Malus domestica accessions, primarily of local Danish origin, were genotyped using 15 microsatellite markers. Ploidy levels were determined by flow cytometry. Special emphasis was given to pedigree reconstruction, cultivar fingerprinting and genetic clustering. A reference set of cultivars, mostly from other European countries, together with a private nursery collection and a small set of Malus sieversii, Malus sylvestris and small-fruited, ornamental Malus cultivars, was also included. The microsatellite markers amplified 17–30 alleles per loci with an average degree of heterozygosity at 0.78. We identified 104 (23%) duplicate genotypes including colour sports. We could infer first-degree relationships for many cultivars with previously unknown parentages. STRUCTURE analysis provided no evidence for a genetic structure but allowed us to present a putative genetic assembly that was consistent with both PCA analysis and parental affiliation. The Danish cultivar collection contains 10% duplicate genotypes including colour sports and 22% triploids. Many unique accessions and considerable genetic diversity make the collection a valuable resource within the European apple germplasm. The findings presented shed new light on the origin of Danish apple cultivars. The fingerprints can be used for cultivar identification and future management of apple genetic resources. In addition, future genome-wide association studies and breeding programmes may benefit from the findings concerning genetic clustering and diversity of cultivars.     

Genetic Diversity and Characterization of a Core Collection of Malus Germplasm Using Simple Sequence Repeats (SSRs)

Simple sequence repeats (SSRs) were used to assess genetic diversity and study genetic relatedness in a large collection of Malus germplasm. A total of 164 accessions from the Malus core collection, maintained at the University of Illinois, were genotyped using apple SSR markers. Each of the accessions was genotyped using a single robust SSR marker from each of the 17 different linkage groups in Malus. Data were subjected to principal component analysis, and a dendrogram was constructed to establish genetic relatedness. As expected, this diverse core collection showed high allelic diversity; moreover, this allelic diversity was higher than that previously reported. Cluster analysis revealed the presence of four distinct clusters of accessions in this collection.     

Genetic diversity and population structure of Moringa oleifera

Moringa is a genus of the tropical flowering plant family Moringaceae containing 13 diverse species. Among the different species, only Moringa oleifera L. is cultivated. This species has great potential in serving as a high-value crop for food, medicinal products, as well as fodder for animals, particularly in developing tropical regions of the world. In this study, the genetic diversity and population structure of world-wide collections of M. oleifera were investigated using DNA markers. A total of 19 microsatellite or simple sequence repeat (SSR) markers along with a partial sequence of the chloroplast gene atpB were used to study genetic diversity within 161 accessions of M. oleifera collected from Asia, Africa, North and South America, and the Caribbean. On average, 8.3 alleles/per SSR were amplified in each accession. A total number of 158 alleles were detected in 131 accessions collected from the wild in Pakistan and from 30 accessions obtained from ECHO (Florida). Observed heterozygosity varied from 0.16 to 0.86, with an average of 0.58, while the average PIC value was 0.59. Partial sequencing of chloroplast genes of 43 of 161 plants generated mixed patterns. These findings have demonstrated that there is a large genetic diversity present in wild collections of M. oleifera collected in Pakistan; whereas low genetic diversity is detected in cultivated accessions obtained from ECHO. Taken together, these results agree with previous reports that M. oleifera is native to the Indo-Pakistan ecological region, and provides sufficient diversity for genetic exploration as well as for genetic improvement efforts.     

Assessment of molecular diversity and evolutionary relationship of kenaf (Hibiscus cannabinus L.), roselle (H. sabdariffa L.) and their wild relatives

Kenaf (Hibiscus cannabinus L.) and roselle (H. sabdariffa L.) are valuable fibre crop species with diverse end use. Phylogenetic relationship of 73 accessions of kenaf, roselle and their wild relatives from 15 countries was assessed using 44 inter-simple sequence repeat (ISSR) and jute (Corchorus olitorius L.) specific simple sequence repeats (SSR) markers. A total of 113 alleles were identified of which 61.95 % were polymorphic. Jute specific SSR markers exhibited high polymorphism and resolving power in kenaf, although ISSR markers exhibited higher resolving power than SSR markers. Number of polymorphic alleles varied from 1 to 5 for ISSR and 1 to 6 for SSR markers. Cultivated species exhibited higher allele polymorphism (57 %) than the wild species (35 %), but the improved cultivars exhibited lower genetic diversity compared to germplasm accessions. Accessions with common genetic lineage and geographical distribution clustered together. Indian kenaf varieties were distinct from cultivars bred in other countries and shared more genetic homology with African accessions. High genetic diversity was observed in the Indian (J = 0.35–0.74) and exotic kenaf germplasm collections (J = 0.38–0.79), suggesting kenaf might have been introduced in India from Africa through Central Asia during early domestication. Genetic similarity-based cluster analysis was in close accordance with taxonomic classification of Hibiscus.     

Genetic diversity and structure of local apple cultivars from Northeastern Spain assessed by microsatellite markers

A set of 493 old and local Spanish accessions of apple (Malus x domestica Borkh) maintained at three collections in Northeastern Spain was studied using 16 simple sequence repeats in order to estimate their genetic diversity and to identify the genetic structure and relationships among their accessions. An additional diverse set of 45 apple cultivars, including old Spanish and international cultivars, was added as reference. Genetic analyses performed by Bayesian model-based clustering revealed a very strong differentiation of two major groups. The first one clustered 159 individuals (52?% of unique genotypes) including local accessions and six old Spanish cultivars. The second major group was formed by 145 individuals, including 38 international reference cultivars and one old Spanish cultivar. Nested Bayesian clustering was applied to those two groups and two and four sub-groups were found at each one, respectively. The identification of private and unique alleles, and the remarkable differences in allelic richness among groups and sub-groups constitute further evidence of a clear genetic structure. The results obtained through the factorial correspondence and analyses of molecular variance confirmed those obtained by Bayesian analyses, revealing moderate but significant differentiation among the two major groups (F _ST?=?0.076) and the six sub-groups (F _ST?=?0.111). Our results highlight that the genetic diversity encompassed by currently cultivated apple accounts only for a small fraction of that existing within the species, and that an important part (??60?%) of the local material analyzed constitutes a good example of genetic distinctness with respect to the main cultivars used in European orchards.     

Genetic Diversity of Some Pear Cultivars and Genotypes Using Simple Sequence Repeat (SSR) Markers

The genetic diversity and relationships among 47 pear cultivars and genotypes (Pyrus spp.), including 4 Japanese pears (Pyrus pyrifolia), 40 European pears (Pyrus communis), 1 Chinese pear (Pyrus bretschneideri) as well as 2 wild relatives (Pyrus salicifolia and Pyrus mazandaranica) were studied using 28 microsatellite primer pairs. A total of 174 alleles were produced at the 28 SSR loci with their sizes ranging from 81 to 290?bp. The number of observed alleles for each locus ranged from 3 (TsuENH014 and TsuENH046) to 12 (NB103a), with an average of 6.21 alleles per locus. In some SSR loci, more than two alleles were amplified in some cultivars and genotypes, suggesting that duplication has occurred in those accessions. This information suggests that at least two genomic regions exist for these loci in the pear genome. The observed heterozygosity (H _o) values of amplified loci ranged from 0.17 (TsuENH006) to 0.97 (NB103a). Shannon's information index (I) value was observed to be highest (2.14) in the NB103a locus, while the TsuENH006 locus had the lowest value with an average of 1.37 among SSR loci. The Dice genetic similarity coefficient ranged from 0.29 (??Nijisseiki?? and P. mazandaranica) to 0.91 (??Chojuro?? and ??Nijisseiki??) among samples. UPGMA cluster analysis showed two major groups corresponding to the Japanese and European pears.     

Microsatellite (SSR) markers reveal genetic identities, genetic diversity and relationships in a Malus×domestica borkh. core subset collection

 A collection of 66 Malus×domestica Borkh. accessions from the USDA-ARS Plant Genetic Resources Unit’s core collection was screened with a set of eight SSR (simple sequence repeat) primers developed at the PGRU in order to determine genetic identities, estimate genetic diversity, and to identify genetic relationships among these accessions. All eight primer pairs generated multiple fragments when used in amplification reactions with DNA from these accessions. High levels of variation were detected with a mean of 12.1 alleles per locus and a mean heterozygosity across all eight loci of 0.693. The eight primer pairs utilized in this study unambiguously differentiated all but seven pairs of accessions in this collection of 66 M.×domestica Borkh. genotypes. The probability of matching any two genotypes at all eight loci in this study was approximately 1 in 1 billion. The markers detected two misnamed accessions in the collection. Genetic-identity data produced a genetic-relatedness phenogram which was concordant with geographic origins and/or known pedigree information. These SSR markers show great promise as tools for managing Malus ex situ germplasm collections as well as for collection and preservation strategies concerning wild Malus populations in situ. Received: 28 March 1998 / Accepted: 29 April 1998     

Characterisation of the virescent locus controlling a recessive phenotype in apple rootstocks (Malus pumila Mill.)

Certain progenies of Malling apple rootstocks (Malus pumila) have been reported to segregate for a virescent trait: leaves are chlorotic at germination or bud break but turn green as the season progresses. The M432 rootstock mapping progeny, from which a linkage map has recently been elaborated with 323 simple sequence repeat (SSR) markers and 3,069 single nucleotide polymorphism (SNP) markers, also segregates for this phenotype. In this investigation, 188 seedlings were scored and, on the basis of a 3:1 segregation, virescence was attributed to the recessive gene (vir) for which the two parents, M.27 and M.116, are heterozygous. At least seven of 28 Malling rootstocks are heterozygous for this apparently deleterious trait. With the published marker data the gene was mapped to linkage group 12, tightly flanked by the SSR CH01g12 and the SNP marker 475880474, and was located in a physical interval of 2.36 Mb on the Golden Delicious genome sequence. A PCR-based marker was developed from the SNP and along with the SSR was scored in a set of Malus rootstock accessions. The screening of this collection demonstrated that those accessions known to be heterozygous at the vir locus all carried the 152 allele of the SSR and the G allele of the SNP, whilst a virescent accession was homozygous for the alleles. The results we present here could help predict the genotype of apple rootstocks at the vir locus, assist in the fine mapping of the vir locus to identify potential candidate genes for the trait and also aid rootstock breeding.     

基于TP-M13-SSR指纹图谱的中国原产苹果属植物分子身份证的建立

利用TP-M13-SSR分子标记方法,构建27份中国原产苹果属植物在12个SSR位点的指纹图谱,运用条码技术生成其分子身份证。12对引物共获得251个等位基因,平均21个。引物多态性好,仅用引物CH05b06即可区分全部供试材料。27份苹果材料在12个SSR位点遗传多样性、多态性信息含量和位点杂合度的变化范围为0.6620～0.9455、0.6327～0.9211和0.6538～0.9319。基于CH05b06位点处获得的指纹谱图即可得到每份供试材料独有的分子身份证。TP-M13-SSR分子标记技术适用于苹果属植物种质资源的指纹图谱构建,利于分子基础数据库的积累。基于苹果种质资源TP-M13-SSR指纹图谱可获得每份苹果种质资源独有的分子身份证。     

The southwestern origin and eastward dispersal of pear (<Emphasis Type="Italic">Pyrus pyrifolia</Emphasis>) in East Asia revealed by comprehensive genetic structure analysis with SSR markers

Pyrus pyrifolia is considered one of the most important cultivated Pyrus species. Hundreds of landraces and bred cultivars have been developed through the natural and artificial hybridizations necessary due to self-incompatibility. In this study, the genetic diversity of 478 Pyrus accessions, including Chinese landraces, bred cultivars, and wild samples, as well as introduced pear cultivars from Japan and Korea, was investigated with a set of 17 simple sequence repeat (SSR) markers distributed across all 17 linkage groups of the pear genome. A total of 121 alleles were detected, including 4 rare alleles with a frequency lower than 5%. Diversity statistics indicated a high level of genetic variation as quantified by the average values of the observed heterozygosity, the expected heterozygosity, and Wright’s fixation index, at 0.76, 0.78, and 0.02, respectively. Population structure and discriminant analysis of principal component analysis implied extensive genetic communication between sand pears in China and revealed four contiguous geographical clusters with overlapping geographical regions. The diversity of the four clusters and approximate Bayesian computation (ABC) indicated that sand pear spread from west to east along the Pearl River and Yangtze River valleys. High diversity and polyphyletic genetic components of cultivars in southwestern China further support southwestern China as the probable center of divergence for Pyrus species. A core collection of 80 out of 470 cultivars was selected, accounting for about 17% of accessions, and capturing 91% of all alleles, including all rare alleles. Our research provides a comprehensive understanding of sand pear germplasm in East Asia and constructs a preliminary core collection, which will be useful for association genetics studies, germplasm conservation, and breeding programs.     

Genetic Diversity Analysis of Tibetan Wild Barley Using SSR Markers

One hundred and six accessions of wild barley collected from Tibet, China, including 50 entries of the two-rowed wild barley Hordeum vulgare ssp. spontaneum (HS), 29 entries of the six-rowed wild barley Hordeum vulgare ssp. agriocrithon (HA), and 27 entries of the six-rowed wild barley Hordeum vulgare ssp. agriocrithon var. lagunculiforme (HL), were analyzed using 30 SSR markers selected from the seven barley linkage groups for studying genetic diversity and evolutionary relationship of the three subspecies of Tibetan wild barley to cultivated barley in China. Over the 30 genetic loci that were studied, 229 alleles were identified among the 106 accessions, of which 70 were common alleles. H. vulgare ssp. spontaneum possesses about thrice more private alleles (2.83 alleles/locus) than HS (0.93 alleles/locus), whereas almost no private alleles were detected in HL. The genetic diversity among-subspecies is much higher than that within-subspecies. Generally, the genetic diversity among the three subspecies is of the order HS > HL > HA. Phylogenetic analysis of the 106 accessions showed that all the accessions of HS and HA was clustered in their own groups, whereas the 27 accessions of HL were separated into two groups (14 entries with group HS and the rest with group HA). This indicated that HL was an intermediate form between HS and HA. Based on this study and previous works, we suggested that Chinese cultivated barley might evolve from HS via HL to HA.     

Development and characterization of microsatellite markers in Indian sesame (Sesamum indicum L.)

The genetic characterization of Indian sesame cultivars and related wild species was analysed using 102 simple sequence repeat (SSR; microsatellite) markers. Of these, 62 were novel sesame-specific microsatellites isolated in the course of the present investigation by constructing genomic libraries. Characterization of the 68 sesame accessions and three related wild species using 72 polymorphic SSR primers resulted in the detection of 170 alleles. The number of alleles ranged from two to four with an average of 2.5 alleles per locus. Polymorphic information content of the markers ranged from 0.43 to 0.88 with an average of 0.66. UPGMA cluster analysis grouped all the accessions into two major clusters with a genetic similarity ranging from 0.40 to 0.91. A moderate to high level of genetic variability was observed. The three wild accessions used in the study formed separate clades and distant genetic relationships were observed between the cultivar lines and wild species. Differentiation of genotypes according to geographical region was not observed. Analysis of molecular variance (AMOVA) analysis revealed that a high percentage of variation was within populations (87.1 %). An overall F _st of 0.11 among the populations indicated low population differentiation. The SSR markers developed will be useful for further genetic analysis, linkage mapping and selection of parents in future breeding programmes.     

Rosaceae conserved orthologous sequences marker polymorphism in sweet cherry germplasm and construction of a SNP-based map

The Rosaceae Conserved Orthologous Set (RosCOS) provides a gene-based genome-wide set of markers that have been used in comparative analyses of peach (Prunus persica), apple (Malus × domestica), and strawberry (Fragaria spp.). In order to extend the use of these RosCOS to sweet cherry (Prunus avium L.), we identified markers that are polymorphic in breeding germplasm. Ninety-five percent (595/627) of previously designed RosCOS primer pairs amplified a product in six sweet cherry cultivars predicted to represent the range of genetic diversity in breeding germplasm. A total of 45% (282/627) RosCOS were polymorphic among the six cultivars, and allele number ranged from 2 to 6, with a genome-wide mean of 2.35. A subset of 92 genome-wide single nucleotide polymorphisms (SNPs) corresponding to 76 RosCOS was analyzed in 36 founder accessions and progeny. The expected and observed heterozygosity suggested that 83% of the RosCOS were in Hardy–Weinberg equilibrium, implying that most RosCOS behave as neutral markers. Principal coordinate analysis (PCO) identified one wild accession and two Spanish landraces that clustered differently from the other accessions. The relatively high number of unique alleles found in the three differentially clustered selections suggested that their use as parents has potential to increase the genetic diversity in future US-bred cultivars. Of the 92 RosCOS SNPs, 81 SNPs that represented 68 genome-wide RosCOS segregated in four mapping populations. These RosCOS were mapped in four F₁ populations, thereby greatly improving the genetic linkage map of sweet cherry.     

Genetic diversity and population structure in the US Upland cotton (Gossypium hirsutum L.)

Key message

Genetic diversity and population structure in the US Upland cotton was established and core sets of allelic richness were identified for developing association mapping populations in cotton.

Abstract

Elite plant breeding programs could likely benefit from the unexploited standing genetic variation of obsolete cultivars without the yield drag typically associated with wild accessions. A set of 381 accessions comprising 378 Upland (Gossypium hirsutum L.) and 3 G. barbadense L. accessions of the United States cotton belt were genotyped using 120 genome-wide SSR markers to establish the genetic diversity and population structure in tetraploid cotton. These accessions represent more than 100 years of Upland cotton breeding in the United States. Genetic diversity analysis identified a total of 546 alleles across 141 marker loci. Twenty-two percent of the alleles in Upland accessions were unique, specific to a single accession. Population structure analysis revealed extensive admixture and identified five subgroups corresponding to Southeastern, Midsouth, Southwest, and Western zones of cotton growing areas in the United States, with the three accessions of G. barbadense forming a separate cluster. Phylogenetic analysis supported the subgroups identified by STRUCTURE. Average genetic distance between G. hirsutum accessions was 0.195 indicating low levels of genetic diversity in Upland cotton germplasm pool. The results from both population structure and phylogenetic analysis were in agreement with pedigree information, although there were a few exceptions. Further, core sets of different sizes representing different levels of allelic richness in Upland cotton were identified. Establishment of genetic diversity, population structure, and identification of core sets from this study could be useful for genetic and genomic analysis and systematic utilization of the standing genetic variation in Upland cotton.     

A genetic linkage map of an apple rootstock progeny anchored to the Malus genome sequence

An apple rootstock progeny raised from the cross between the very dwarfing ??M.27?? and the more vigorous ??M.116?? (??M.M.106???×???M.27??) was used for the construction of a linkage map comprising a total of 324 loci: 252 previously mapped SSRs, 71 newly characterised or previously unmapped SSR loci (including 36 amplified by 33 out of the 35 novel markers reported here), and the self-incompatibility locus. The map spanned the 17 linkage groups (LG) expected for apple covering a genetic distance of 1,229.5?cM, an estimated 91% of the Malus genome. Linkage groups were well populated and, although marker density ranged from 2.3 to 6.2?cM/SSR, just 15 gaps of more than 15?cM were observed. Moreover, only 17.5% of markers displayed segregation distortion and, unsurprisingly in a semi-compatible backcross, distortion was particularly pronounced surrounding the self-incompatibility locus (S) at the bottom of LG17. DNA sequences of 273 SSR markers and the S locus, representing a total of 314 loci in this investigation, were used to anchor to the ??Golden Delicious?? genome sequence. More than 260 of these loci were located on the expected pseudo-chromosome on the ??Golden Delicious?? genome or on its homeologous pseudo-chromosome. In total, 282.4?Mbp of sequence from 142 genome sequence scaffolds of the Malus genome were anchored to the ??M.27???×???M.116?? map, providing an interface between the marker data and the underlying genome sequence. This will be exploited for the identification of genes responsible for traits of agronomic importance such as dwarfing and water use efficiency.     

Phylogeny of the Malus (Apple Tree) Species,Inferred from the Morphological Traits and Molecular DNA Analysis

Some debated issues of the genus Malus (apple) taxonomy were examined using a variety of species from the collection of the Maikop Experimental Station, Vavilon Research Institute of Plant Industry (Krasnodar krai). Phylogenetic relationships among these species were studied using traditional analysis of morphological traits, RAPD, and complete sequencing of the 5"- internal transcribed spacer (ITS1), 5.8S rRNA, 3"- internal transcribed spacer (ITS2) (constituting a cluster of the rRNA genes), and the terminal fragment of the matK gene encoding chloroplast maturase. The results showed that the Sorbomalussection was polyphyletic; the American apple M. fusca was closely related to the species contributing to the East Asian center of the genus origin, and the American speciesM. angustifolia, M. coronaria, and M. ioensis were closely related to the M. trilobata relict species, whose assignment to the genus Malus is debated by some authors. Molecular analysis of the species relationships showed that the Middle Asian apple M. sieversii is the species from which apple domestication started.