# Genetic characterization of worldwide Prunus domestica (plum) germplasm using sequence-based genotyping.

Tetyana Zhebentyayeva, Vijay Shankar, Ralph Scorza, Ann Callahan, Michel Ravelonandro, Sarah Castro, Theodore Dejong, Christopher A Saski, Chris Dardick

► To cite this version:

Tetyana Zhebentyayeva, Vijay Shankar, Ralph Scorza, Ann Callahan, Michel Ravelonandro, et al. Genetic characterization of worldwide Prunus domestica (plum) germplasm using sequence-based genotyping. Horticulture research, 2019, 6, pp.12. 10.1038/s41438-018-0090-6. hal-02619720

HAL Id: hal-02619720

https://hal.inrae.fr/hal-02619720

Submitted on 25 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Genetic characterization of worldwide *Prunus domestica* (plum) germplasm using sequence-based genotyping

Tetyana Zhebentyayeva<sup>1,2</sup>, Vijay Shankar<sup>1,2</sup>, Ralph Scorza<sup>3</sup>, Ann Callahan<sup>3</sup>, Michel Ravelonandro<sup>4</sup>, Sarah Castro<sup>5</sup>, Theodore DeJong<sup>5</sup>, Christopher A. Saski<sup>2</sup> and Chris Dardick<sup>3</sup>

**Abstract**

*Prunus domestica* commonly known as European plum is a hexaploid fruit tree species cultivated around the world. Locally it is used for fresh consumption, in jams or jellies, and the production of spirits while commercially the fruit is primarily sold dried (prunes). Despite its agricultural importance and long history of cultivation, many questions remain about the origin of this species, the relationships among its many pomological types, and its underlying genetics. Here, we used a sequence-based genotyping approach to characterize worldwide plum germplasm including the potential progenitor Eurasian plum species. Analysis of 405 DNA samples established a set of four clades consistent with the pomological groups Greengages, Mirabelles, European plums, and d'Agen (French) prune plums. A number of cultivars from each clade were identified as likely clonal selections, particularly among the "French" type prune germplasm that is widely cultivated today. Overall, there was relatively low genetic diversity across all cultivated plums suggesting they have been largely inbred and/or derived from a limited number of founders. The results agree with *P. domestica* having originated as an interspecific hybrid of a diploid *P. cerasifera* and a tetraploid *P. spinosa* that itself may have been an interspecific hybrid of *P. cerasifera* and an unknown Eurasian plum species. The low genetic diversity and lack of true wild-types coupled with the known cultivation history of Eurasian plums imply that *P. domestica* may have been a product of inter-specific cross breeding and artificial selection by early agrarian Eurasian societies.

**Introduction**

European plum *Prunus domestica* L. is a polymorphic allopolyploid (hexaploid) species ($2n = 6x = 48$) commercially grown worldwide for a variety of uses including fresh fruit, prunes, distilling, and as processed additive ingredients. This plum species, commonly referred to as "European plums" or "prune plums", is distinct from the large round diploid "Japanese plums" (*Prunus salicina*) which are widely grown for fresh market consumption.

*P. domestica* plums are critical components of human diets as prunes have been shown to have a broad range of health promoting activities including protection against cardiovascular disease, diabetes, digestive disorders, and osteoporosis¹,². Historical evidence suggests *P. domestica* and other Eurasian plum species including *Prunus cerasifera* (Cherry plum) and *Prunus spinosa* ("Sloe" or "Blackthorn" plum) were important to the development of early European societies. Stones from all three species have been uncovered from Neolithic archeological sites dating back to 4000–6000 BC in Germany and the Ukraine³.

The plum or "prumnon" was first recorded in Archipoledus's "Pollux" written in the 7th century BC³. Around the birth of Christianity, the Roman author Pliny gave the

<figure><img src="image_2.png" /></figure>

<footnote>Correspondence: Chris Dardick (chris.dardick@ars.usda.gov)

¹The Schatz Center for Tree Molecular Genetics, Department of Ecosystem Sciences and Management, The Pennsylvania State University, University Park, PA, 16802, USA

²Genomics and Computational Biology Laboratory, Clemson University, Clemson, SC 29634, USA

Full list of author information is available at the end of the article.</footnote>
first clear account of “Damascus plum” where he descri-
bed a wide variety of cultivars with different fruit types
and colors, their uses, and grafting onto “sorb” rootstocks
(probably *Prunus cerasifera* or *Prunus spinosa*), which is
still a common practice today. Roman depictions of plum
fruit provide further evidence that at least some of the
*P. domestica* varieties grown at that time were very similar
or even possibly the same as those still grown today.
References to different types of plums being imported
from Syria or Persia and the fact that plums were not
native to some of the Neolithic sites indicates that they
were under human selection and cultivation in the Cau-
casus region long before their introduction into Europe.
The propensity of *P. domestica* to form root suckers
would have provided a simple and efficient way to dis-
tribute clonal material across long distances even before
grafting techniques had been developed.

P. domestica is thought to have originated from the
Middle East, specifically, in the area south of the Caucasus
between the Black Sea and the Caspian Sea that encom-
passes Georgia, Armenia, Azerbaijan, and the northern
plateau of Iran⁴. Vavilov⁵ placed the center of its origin in
the region south of Caucasian Mountains through the
Caspian Sea, in the area overlapping the distribution of
P. cerasifera L. and P. spinosa L; which were likewise
concurrently distributed across Europe for cultivation,
consumption, and use as rootstocks. Widespread feral
growth, the apparent lack of natural stands in forests, and
human movement of both cultivated and wild type plums
across the continent has made determining the origin
of P. domestica extremely difficult. The existence of wild
forms of P. domestica in the Caucasus region was recur-
rently reported in the 1930s based on morphological
traits but failed validation from cytogenetic analyses⁵,⁶.
Recent discovery of wild P. domestica stands in Xinjiang
province in northwestern China prompted the hypothesis
that European plum may be of Chinese origin⁷,⁸. How-
ever, narrow molecular diversity of botanical material
from Xinjiang found that the Chinese P. domestica was
likely feral resulting from escape and naturalization⁹.

Consequently, the origin of *P. domestica* has remained
a matter of debate for nearly a century. A major compli-
cation is the wide range of intraspecific variations and
transitional forms¹⁰. Traditionally, plum cultivars have
been divided into a number of different pomological
groups; small fruited mirabelle plums, damsons, small
wild plums or bullaces, greengages, prune plums, and
large-fruited European plums (Fig. 1)¹¹⁻¹³. However, clear
delineations between these groups are extremely difficult
to define due to the inherent phenotypic variability across
the germplasm³,⁴,¹³.

Molecular phylogenetic analyses have more recently
confirmed the close relationships among *P. cerasifera*,
*P. insititia*, *P. domestica* L., *P. spinosa* L., and the

“Marmot” plum *Prunus brigantina* and dated their
divergence from a common Eurasian ancestor to the
Oligocene (31Myr)¹⁴. It is now generally accepted that *P.
domestica* is most likely an interspecific hybrid of *P. cer-
asifera* and *P. spinosa*. However, defining the precise
genetic history of these events has been complicated by
the fact that *P. cerasifera* and *P. spinosa* are found as
diploid and polyploid individuals, sometimes having
morphologies very similar to *P. domestica*. In fact,
hybridization experiments and restriction site analysis of
Ribosomal RNA genes indicated that *P. spinosa* itself may
be an inter-specific hybrid having *P. cerasifera* as one of
the parents¹⁵⁻¹⁹. Crane and Lawrence²⁰ and Rybin⁶ per-
formed direct interspecific hybridization experiments
between *P. cerasifera* (2x), and *P. spinosa* (4x) and found
that, while the vast majority of the resulting offspring
were infertile, a very small number of fertile hybrids that
were morphologically similar to *P. domestica* could be
obtained⁶. Still, many questions remain unresolved as to
the origins of the various pomological plum groups and
whether they represent independent inter-specific hybr-
dization events or arose from the same or limited number
of events³,⁶,¹²,¹³,²⁰.

More recently, genetic diversity of Eurasian plum vari-
eties and species has been studied using bi-parentally
transmitted nuclear and maternal cytoplasmic (chlor-
eloplast, mitochondria) markers. Microsatellite (SSR) and
Inter-simple Sequence Repeat (ISSR) markers were used
for characterization of *P. domestica* germplasm across
Europe in several geographically distant countries<sup>21-24</sup>, in
Turkey<sup>25</sup> and China<sup>9,26</sup>. This approach enabled a recon-
struction of parental lineages among related European
plum cultivars<sup>27</sup> and molecular characterization of “Reine
Claude” types in Spain<sup>28</sup>. Using the allelic variants of the
S-RNase gene, groups of cross-compatible and incompa-
tible cultivars were generated for Latvian plums<sup>29</sup>, and
selected cultivars propagated in Hungary and Slova-
kia<sup>30,31</sup>. In several studies cytoplasmic (chloroplast) and
nuclear SSR markers were used to reveal genetic diversi-
fication within *P. domestica* and to attempt to resolve
its origin. Based on chloroplast markers it was shown
that *P. cerasifera* was a likely progenitor, at least in the
maternal lineage<sup>21,32</sup>. Two main cpDNA haplotypes were
identified, however, major haplotypes were not clearly
associated with pomological groups<sup>21</sup>. Nuclear DNA
markers, on the other hand, indicated a potential con-
tribution of *P. spinosa* to the *P. domestica* genome<sup>14,21</sup>.
Unfortunately, the limited number of molecular markers
used for characterization of plum germplasm to date,
have not yet provided a clear picture of the genetic rela-
tionships among Eurasian plum germplasms.

Restriction site-associated sequencing has been an
effective method for identifying and screening high-
resolution polymorphisms within and between
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Fruit-type diversity of *P. domestica* and potential ancestral species. Illustrations were reprinted from Hedrick¹ with the exception of *P. spinosa* and *P. brigantina*. These early 19th century illustrations were created by the artist Pancrace Bessa and published in "Treats the trees and shrubs grown in France in the open ground" by [H.L.] Duhamel [du Monceau; ed. by Etienne Michel (vol 1–4) and J.L.A. Loiseleur-Deslongchamps (volumes 5–7)]</figcaption></figure>

populations, ecotypes, and species³³,³⁴. In the present study, we generated a set of sequence-based SNP markers densely distributed across the *Prunus* genome in order to: (1) investigate genetic relationships among cultivated plums from different pomological groups that originated under the influence of different ecological factors and were propagated in different geographical regions; (2) estimate the extent of variation among *P. domestica* germplasm and identify diagnostic molecular markers; and (3) estimate the potential contributions of other plum species *P. cerasifera* and *P. spinosa* to the nuclear genome of hexaploid *P. domestica*. The resulting phylogenetic relationships along with our understanding of the origins, history, and future of plum cultivation are discussed.

## Results

### SNP discovery

In total, we sequenced 405 DNA samples representing different pomological groups of *P. domestica* and other Eurasian *Prunus* species including the potential *P. domestica* progenitors *P. cerasifera* and *P. spinosa* (Table 1, Table S1).

Five *Pst1*-digested genomic library pools were sequenced across five Illumina HiSeq2500 lanes producing 2051 million (mln) raw reads. Of these, 97.8% (2005 mln) reads were retained after filtering for base call accuracy, barcode match, and presence of the *Pst1*-restriction site (CTGCA/G). After demultiplexing, the average numbers of reads per DNA sample was 4.93
million reads. Only 5.6% (23 out 405 accessions) had fewer than 1 mln reads. All sequences were aligned

<table><caption>**Table 1** Plum pomological groups and species in this study</caption><thead><tr><td>Plum group/species</td><td>Code</td><td>Accessions</td></tr></thead><tbody><tr><td>European plums</td><td>EUR</td><td>195</td></tr><tr><td>Greengages</td><td>GRG</td><td>46</td></tr><tr><td>Mirabelles</td><td>MIR</td><td>15</td></tr><tr><td>Prunes</td><td>DAP</td><td>107</td></tr><tr><td>P. insititia</td><td>Pi</td><td>9</td></tr><tr><td>P. spinosa</td><td>Psp</td><td>20</td></tr><tr><td>P. cerasifera</td><td>Pc</td><td>10</td></tr><tr><td>P. brigantina</td><td>Pbr</td><td>2<sup>a</sup></td></tr><tr><td>P. simonii</td><td>Psi</td><td>1</td></tr><tr><td>Total</td><td></td><td>405</td></tr></tbody></table>

against the *P. persica* genome assembly v2.1 for SNP identification. The average depth of stacks for variant calling was 221 and varied from 6 to 778. A total of 178,403 unfiltered sites were identified including monomorphic SNPs across all plum accessions. Filtering for SNPs present in more than 80% of accessions resulted in a dataset of 129,110 SNPs that were evenly distributed throughout the *P. persica* pseudochromosomes with little or no large-scale bias (Fig. 2). To keep uncommon variants that are specific to diploid *Prunus* species and potentially contributing to the hexaploid genome, we omitted two common filters, minor allelic frequencies and missing individual rate.

A total of 400 sequence variants derived from the chloroplast (cp) genome were also identified. Of these, 14 SNPs were present in more than 80% of accessions (Table S2). After filtering for Minor Allelic Frequencies (MAF) > 0.05, the only SNP with allelic frequency of 0.324 was a synonymous C/T substitution at 57,543 nt (+) within the large subunit of ribulose 1,5-bisphosphate carboxylase/oxygenase (rbcL) gene.

<figure><img src="image_4.png" /></figure>

<footnote><sup>a</sup>Including *P. brigantina* × *P. cerasifera* hybrid</footnote>
<figure><img src="image_5.png" /></figure>

## Genetic diversity and phylogenetic analysis

We investigated genetic relationships among all 405 plum accessions by performing a hierarchical clustering analysis (between-group average linkage method, or UPGMA) on a dissimilarity matrix calculated in R. In total, 103,382 polymorphic SNPs out of 129,110 sites exported from Stacks were used for phylogenetic classification (Fig. 3; Figure S1). The biological and technical replicates included in the inter-plate control and analysis were used to assess dendrogram quality. Results showed robust and consistent clustering of identical samples and an overall absence of a strong Illumina lane or geographic location effect on genetic signals.

## Prunus domestica germplasm structure

Distinct clades among domesticated hexaploid plums reflected some, but not all, of the known pomological groups including European plums (EUP), Mirabelles (MIR), Greengages (GRG), and d'Agen prune plums (DAP) (Fig. 3; Figure S1; Table S3). Differentiation of other pomological groups including *P. insititia*,

"Damsons", and "Bullaces" was not supported as such identified individuals did not form discrete clades. Among the four identified clades, subclades consisting of technical replicates, biological replicates and/or known clonal varieties could be distinguished in all pomological groups. With very few exceptions, varieties with known parentage clustered as expected within clades supporting the accuracy of the phylogenetic classification (Table S1). Using the "collapse subtree option" in Geneious³⁵ we estimated a cutoff value for discriminating groups of clonal material using the genetic distances among DAP germplasm including samples from known clones, technical replicates, and seedlings. A subclade comprising 34 individual trees derived from self-pollination of the commercial variety "Improved French" and a subclade of 46 samples that included known DAP clones and technical replicates collapsed at the thresholds of branch lengths 0.392 and 0.373, respectively, while advanced generations of hybrids having d'Agen in their pedigree collapsed at a threshold of 0.454. Consequently, a conservative branch length cutoff value of 0.363 was chosen to distinguish clones across the rest of the dataset. This resulted in the identification of 16 clonal groups that were distributed in all four clades. These results are consistent with the known history of plum domestication as numerous cultivars have been selected by propagation of sports or suckers spontaneously showing new characteristics. It is important to note that this analysis may not have captured all clones as a handful of samples that are likely to be clonal did not meet the threshold (Table S1).

Due to the large number of samples used in this study, a detailed analysis of individual cultivars and their genetic relationships is not reported here. General results from the four distinguishable *P. domestica* clades are given below.

### European plums (EUP)

European plums include many old English and Eastern European cultivars. Members of this group also show a high degree of morphological diversity ranging from small round bullaces to large oblong plums. The EUP clade showed the highest level of genetic diversity. Nine clonal clades were identified among EUP types. The largest consisted predominantly of "Pozegaca" or "Quetsche" types which are round, small fruited plums commonly referred to as "German prunes". These are believed to be ancient cultivars that are still preferred in Eastern Europe but have largely been pushed out of production due to the plum pox virus epidemic³⁶. Other smaller clonal clades consisted of varieties of French origin with mixed but often related names and the so-called "Italian" plums which produce large, oblong shaped fruit (Table S1)³⁷. A set of EUP cultivars with diverse origins including "Ruth Gerstetter", "Topfirst", "Perdrigonne" and
"Cacanska Rana" formed a more distant and poorly defined EUP subclade. "Ruth Gerstetter" is recorded to have been a cross between "Bonne de bry" and "Czar". Our results show that "Bonne de bry" is likely a Greengage type. The other parent, "Czar", was not included in this study, however, other plums with "Czar" parentage including "Edda" and "Herman" fall within the same clade. Likewise, "Topfirst" has "Ruth Gerstetter" as a parent and also falls within this group. One possibility is that this EUP subclade is comprised of plums with mixed heritage.

### Mirabelles (MIR)

Mirabelle plums are small, typically red or yellow in color, and slightly oval in shape. They are often found growing feral and most commonly prepared as jams or jellies. The MIR clade represented a very narrow range of germplasm that was not far removed from EUP plums. Three small clonal groups consisting of 3–6 varieties each were identified (Table S1).

### Greengages

The name "Greengage" comes from Sir William Gage who supposedly imported these small, often green-fruited types from France in the early 1700s³⁸. Many of the cultivars carry their French name "Reine Claude" in honor of Queen Claude who ruled France in the early 1500s. Three groups of clonal material were identified with a single large group (18 members) consisting of cultivars named as "Reine Claude" types. These findings suggest that a significant portion of the commercial "Greengage" germplasm is clonal.

### D'Agen prune plums (DAP)

DAP is the most economically important group of plums and includes the vast majority of commercially growth prunes. The DAP clade (144 accessions) was composed of commercial d'Agen prunes from different geographical regions, progeny from self-pollination of the d'Agen prune "Improved French", cultivars with known d'Agen types in their pedigree and a few unclassified cultivars. Nearly all individuals named as d'Agen types were identified as clones, with the exception of three d'Ente varieties collected from INRA that fell into the EUP clade and d'Ente-jaune which is likely a self-pollinated d'Agen seedling (Table S1). These results indicate that worldwide commercial prune production is predominantly a monoculture.

### Eurasian (EUR) plums and *P. domestica* ancestry

Several diploid plum species were examined including accessions of *P. simonii* (native to China), the interspecific hybrid rootstock *P. cerasifera × P. munsoniana* ('Marianna 2634') (note: *P. munsoniana* is native to North America), several diploid *P. cerasifera* selections, *P. brigantina*, and a hybrid of *P. cerasifera × P. brigantina*.

Tetraploid *P. spinosa* accessions were collected from the U.S., Sweden, and Portugal. Accessions recorded as *P. insititia* included the rootstocks "Saint Julien" as well as several varieties collected from INRA.

The non-Eurasian plum *P. simonii* formed a distinct out-group and, to a lesser extent, so did the *P. cerasifera × P. munsoniana* hybrid "Marianna 2624". In contrast, all *P. cerasifera*, *P. spinosa*, and *P. brigantina* selections displayed relatively short genetic distances to *P. domestica*. These species did not form discrete clades and three potentially clonal groups were identified (Table S1). These include a group of *P. cerasifera* accessions from the UC Davis germplasm repository, a set of *P. spinosa* accessions from Sweden, and a group with apparently mixed identities that included "Krikon" plum and "Mirabelle sans nom" from INRA, a *P. spinosa* accession from Sweden, and two *P. cerasifera* replicates from UC Davis. *P. insititia* samples did not form a distinct clade and were scattered among autochthonous *P. domestica* and EUR germplasm. These results indicate a relatively close genetic similarity among all Eurasian plums regardless of ploidy level. The mixed groupings of the Eurasian plum accessions also indicates a high likelihood that at least some of the material is mis-identified or was subject to collection/ processing errors. Interestingly, three *P. spinosa* genotypes from Portugal grouped with autochthonous *P. domestica* accessions "Wegierka wiedenska" and *P. insititia* "Tersen" from southern Sweden. These findings indicate the potential that at least some *P. domestica* and/or *P. insititia* accessions may represent ancestral germplasm, however, we could not definitively rule out the possibility that these accessions resulted from inter-specific hybridization with *P. domestica*.

### Germplasm stratification—$F_{ST}$ statistics and PCA analysis

We used Wright's F-statistics ($F_{ST}$) to explore the degree of differentiation between pomological groups and Eurasian plum species. $F_{ST}$ estimates the level of reduction in heterozygosity when compared to the expected Hardy-Weinberg equilibrium³⁹. Consistent with random mating, the inbreeding coefficient $F_{IS}$ was low among *P. cerasifera* and negative within *P. domestica* pomological groups. The overall fixation index ($F_{ST}$) for *P. domestica* groups was estimated to be 0.005–0.033 (Table 2) and, in agreement with small pairwise genetic distances, reflecting weak differentiation among the identified *P. domestica* clades.

The $F_{ST}$ values were not uniform across the dataset as differentiation between EUP and DAP was significantly stronger than between EUP and GRG or MIR groups; 0.018, 0.006 and 0.005, respectively. The highest differentiation was found between DAP and GRG and DAP and MIR (0.033 and 0.030, respectively). Overall the
$F_{ST}$ index was largely in the range of 0.05–0.15 for pairwise comparisons between DAP, GRG, MIR, or EUP vs. potential progenitor species *P. cerasifera* and *P. spinosa*. This range represents moderate to low differentiation according to Hartl and Clark⁴⁰ indicating that *P. spinosa* is less differentiated from *P. domestica* groups than *P. cerasifera* (Table 2). The $F_{ST}$ index for the EUP group vs. *P. spinosa* was significantly less than other pomological groups (0.026 vs. 0.05); suggesting that EUP is less diverged from the *P. spinosa* ancestor. This difference

became more obvious after removing diploid *Prunus* species as well as the interspecific hybrids “Marianna 2634” and the *P. brigantina* × *P. cerasifera* hybrid from the dataset (Table S3).

Next, we conducted Principle Component Analysis (PCA) to estimate the genetic stratification of the germplasm. We compared the patterns of genetic variation using the entire dataset of 129,110 SNPs vs. filtered datasets generated with different minimum minor allelic frequencies (MAF) and missing rates (MAF = 1%, 3%, or 5% and missing rate = 10–30%) (data not shown). Optimal resolution of plum species and pomological groups was achieved at MAF = 5% and a call rate >80% (missing rate = 20%) resulting in an optimized set of 24,978 SNPs which were used for further genetic analysis.

Principal components analysis (PCA) of the 24,978 SNP dataset identified 32 significant principal components (PCs) (explaining 43.72% of total variance). The first six PCs (9.6%, 6.0%, 4.6%, 2.5%, 2.3%, 1.8%–26.8% of variation in total) were selected for classifying/clustering accessions into groups (Fig. 4, PC5 not shown; Table S4).

A large proportion of the genetic variation remained unstructured. In agreement with pairwise $F_{ST}$ statistics, DAP and potential progenitors of hexaploid plum, *P. cerasifera* and *P. spinosa*, were well differentiated along the first two principle components PC1 and PC2. Diploid *Prunus* species (*P. cerasifera* and *P. brigantina*) and

<table><caption>**Table 2** Number of private alleles and pairwise $F_{ST}$ between plum pomological groups and species (interspecific hybrids and underrepresented species deleted)</caption><thead><tr><th></th><th>Private sites<sup>a</sup></th><th>DAP</th><th>GRG</th><th>Pc</th><th>MIR</th><th>Psp</th><th>Pi</th></tr></thead><tbody><tr><td>EUP</td><td>16232</td><td>0.018</td><td>0.006</td><td>0.044</td><td>0.005</td><td>0.026</td><td>0.009</td></tr><tr><td>DAP</td><td>3452</td><td></td><td>0.033</td><td>0.079</td><td>0.030</td><td>0.056</td><td>0.026</td></tr><tr><td>GRG</td><td>1990</td><td></td><td></td><td>0.086</td><td>0.023</td><td>0.054</td><td>0.021</td></tr><tr><td>Pc</td><td>3999</td><td></td><td></td><td></td><td>0.113</td><td>0.057</td><td>0.069</td></tr><tr><td>MIR</td><td>1186</td><td></td><td></td><td></td><td></td><td>0.058</td><td>0.018</td></tr><tr><td>Psp</td><td>12930</td><td></td><td></td><td></td><td></td><td></td><td>0.023</td></tr><tr><td>Pi</td><td>3269</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

<figure><img src="image_6.png" /></figure>

<footnote>a17,339 SNPs in pairwise haplotype file</footnote>
*P. spinosa* were differentiated along PC4 (2.5% of variance). Plotting PC1 against PC3 and PC6 (4.6% and 1.8% of variance, respectively) gave only weak differentiation of GRG and MIR from EUP. This was in contrast to the relatively strong differentiation observed for DAP under all combinations.

## Discriminant analysis of pomological groups and progenitor's species

Random Forest (RF) is a machine-learning approach that allows classification of individuals into groups and the ranking of the relative importance of each SNP based on classification error. Each SNP is essentially removed one-by-one and the significance of its contribution evaluated. This approach was used to confirm whether the GRG and DAP groups can be effectively discriminated from the rest of the Eurasian plum germplasm. First, we applied the RF algorithm and constructed a multidimensional plot for the dataset composed only of three plum groups—DAP, GRG, and combined plums (cPLM) representing the rest of accessions, i.e., EUP, MIR, autochthonous plums, and related species (Fig. 5). We generated seven subsets of the most significant SNPs for separation of GRG and DAP from the rest of germplasm. The out-of-bag error (OOB) was in the range of 4.4–6.7%, supporting the differentiation of GRG and DAP. Orthogonal scattering of GRG and DAP likely reflected selection of different variables (SNPs) associated with breeding for different fruit types. This is in agreement with PCA results and F-statistics.

To extract a subset of the most significant genetic markers for differentiating these groups, we ran a RF analysis for DAP or GRG vs. the rest of the germplasm. After each round of calculations, we extracted the top 20 most significant SNPs and then combined them, retaining only unique marker names. This resulted in a list of 76 significant markers for DAP plums. In a similar manner, we generated a list of the 85 most significant markers for discriminating GRG from other plum groups. Next, we combined these SNPs and verified that they would reproduce the scattering pattern for GRG and DAP on a PCA plot (Figures S2, S3; Table S5). Lastly, Procrustes analysis was used to show the superimposed distributions of accessions for groups between the unreduced SNP set and the subset of 161 most informative SNPs (*P*-value < 0.001 by Monte Carlo permutation procedure, Figure S3). The results show that the subset of 161 markers may have prediction value for classification of plum germplasm by fruit type (GRG vs. DAP vs. EUP vs. MIR) as well as for possible genomic selection for fruit type in plum breeding programs. Next, we mapped the most significant SNPs against the peach genome (Fig. 1; Table S5). The SNPs were non-uniformly distributed on all eight peach pseudochromosomes, indicating

that specific chromosomal regions may have been under selection in the breeding for GRG and DAP types (Fig. 1; Figure S4; Table S5). We analyzed the genes harboring these 161 SNPs to look for potential functional relationships (Table S4). The most abundant gene categories included sugar metabolism and transporters, cell wall modification enzymes, cell cycle regulation, RNA processing, and development.

<figure><img src="image_7.png" /><figcaption>Fig. 5 Multidimensional plot of the proximity matrix by randomForest analysis calculated using dataset of 405 plum accessions and 24,978 SNPs (ntree = 1000). Color indicates main plum groups used in analysis: green—greengages (GRG), red—d'Agen prunes (DAP), blue—remaining plum accessions germplasm (cPLM)</figcaption></figure>

## A rare cp marker is associated with d'Agen prunes

Due to insufficient depth of reads, we failed to establish a cp haplotype for 34 accessions (8.40%), so these were removed from the dataset. In total, 371 samples were retained for analysis. A C/T substitution in the large subunit of the RUBISCO gene was found in 29.9% of the accessions, while 61.7% of cultivars had the cp haplotype C. Interestingly, the haplotype T was almost exclusively found in DAP, among which only 6.50% of samples had haplotype C. BLAST searches against available cp genomes in NCBI revealed that this residue is highly conserved among Angiosperms, suggesting the T haplotype is rare (data not shown). To uncover potential associations between nuclear markers and the cp haplotype we conducted Canonical Correspondence Analysis (CCA) using the vegan package in R⁴¹. The Permutation test (*P* < 0.001) confirmed a significant association with a subset of nuclear markers. These markers were scattered along the CCA1-axis, potentially reflecting their
<figure><img src="image_8.png" /><figcaption>**Fig. 6** Canonical correspondence analysis (CCA) plot. Correlation matrix calculated using dataset of 371 plum accessions genotyped for 24,978 nuclear SNP markers (unconstrained axis). Constrained ordination was done along the CCA1 axis reflecting cp haplotypes C/T (no missing values allowed). *P* ≤ 0.001 (n = 999). Color indicated cp haplotype C (blue) and T (red)</figcaption></figure>

simultaneous introgression into plum germplasm along with the T haplotype (Fig. 6). In contrast, nuclear markers associated with the C haplotype were distributed along the CA1-axis in an unconstrained fashion, indicating its association with two or more different sources of nuclear markers. The distribution of this cp haplotype mimics the strong PC differentiation observed for DAP relative to GRG, MIR, and EUP (Figure S5).

# Discussion

## Differentiation of *P. domestica* pomological groups

Overlapping distribution of morphological traits among pomological groups of *P. domestica* as well as among seedlings derived from self-pollinated cultivars prompted some early authors to propose that subdivision of *P. domestica* cannot be justified from a genetic point of view<sup>42,43</sup>. Based on the analysis of worldwide germplasm presented here, we report clear genetic evidence that at least some pomological groups of *P. domestica* have distinctive genotypic signatures that can be used for assignment of unknown accessions and paves the way for future genome-enabled prediction of agronomic traits in breeding material.

Two pomological groups, DAP and GRG, were clearly separated from each other and from the rest of the *P. domestica* germplasm. GRG and DAP formed individual clades on the dendrogram derived from the complete dataset of 129,110 polymorphic SNPs and were readily separated by PCA analysis using a representative dataset of 24,498 SNPs. Accordingly, the highest fixation index between all pomological groups ($F_{ST} > 0.03$) indicated that DAP and GRG are the most differentiated among

Eurasian plum germplasm. The relatively strong differentiation of DAP and GRG was also supported by RF analysis, an advanced machine-learning approach<sup>44</sup>. Recently, successful implementation of RF for site-by-site classification of Atlantic and Chinook salmon species allowed the establishment of a relevant, non-redundant and maximally reduced panel of genetic markers<sup>45</sup>. Identification of informative DNA markers useful for differentiating populations was also achieved in cattle breeds<sup>46</sup> and perennial ryegrass<sup>47</sup>. This report provides support for the efficiency of RF in analyzing and unmasking NGS datasets.

We were able to computationally identify two subsets of SNPs that were most informative for assigning unknown accessions to the GRG and DAP groups (subsets 76 and 85, respectively). Based on the uneven distribution of these SNPs across the genome, we hypothesize that these genomic regions may harbor important genes for fruit type or other agronomic traits that are distinct between the DAP and GRG groups. Historically, these pomological groups were selected for different fruit usage such as drying vs. fresh consumption. In many cases these SNPs were within or proximal to orthologues of proteins known to have key roles in sugar metabolism and transport which could underlie important fruit quality differences. Thus, these SNPs may serve as good markers for screening hybrid material in breeding programs aimed at dried or fresh types of fruits.

Despite the ability to clearly resolve some *P. domestica* pomological groups, the results revealed relatively small genetic distances between them. Thus, *P. domestica* was likely derived from a limited set of founders where breeding consisted primarily of self-pollination and/or hybridization among selected siblings. This conclusion was made possible by the inclusion of siblings derived from a self-pollinated "Improved French" individual. This also allowed us to estimate branch length cutoffs to distinguish clonal accessions from siblings. Clonal subclades were identified within all pomological groups and Eurasian plum accessions, consistent with the historical practice of making clonal selections from vegetative buds or sports. Importantly, nearly all commercially grown DAP varieties from France, USA, Argentina, and Australia were found to be clonal, establishing that the plum industry worldwide is predominantly a monoculture. A surprising finding was the unique cp haplotype T that prevailed in DAP plums but was not detected in related *Prunus* species. This suggests that DAP types were historically not cross hybridized to EUP, GRG, or MIR plums and were derived from a distinct hybridization event. Interestingly, DAP plums are very old varieties as they date to before the 12th century where they were reportedly brought to France directly from Syria<sup>37</sup>. Thus, their differentiation from other *P. domestica*
cultivars despite the relatively short genetic distances between them raises the possibility that this cp haplotype is a remnant of one of the inter-specific events that gave rise to *P. domestica* in the Caucasus. Given the lack of evidence for this haplotype in all Eurasian accessions tested, it is possibly derived from a theoretical 3rd Eurasian genomic contributor. Such a possibility is supported by CCA analysis that found sets of nuclear markers were associated with the cp T haplotype. Future genome sequencing of *P. domestica* and its progenitors will be necessary to confirm these hypotheses.

The sequence-based genotyping method used here has limitations caused by the lack of a plum reference assembly and the reliance on the *P. persica* genome for reference assembly of the *P. domestica* sequences. As shown for Sugarcane and Saccharum Complex, the quality of the reference assembly strongly influences GBS performance in polyploids⁴⁸. Consequently, we likely missed an opportunity for discovery of polymorphic markers derived from intergenic regions, as well as SNPs that can discriminate between subgenomes (homo-eoSNPs). Despite this limitation, the results of our analyses were highly consistent with known genetic relationships and support the robustness of the SNP calling pipeline.

## Relationship of *Prunus domestica* to other Eurasian plums

Questions regarding the origin of *P. domestica* have been studied and debated for over a century. The wide range of phenotypic diversity observed in *P. domestica* and its Eurasian relatives along with their overlapping distributions has made resolving these questions particularly difficult. A handful of competing theories have been proposed that are not all mutually exclusive: [1] *P. domestica* is entirely of *P. cerasifera* origin and represents an autopolyploid event of this species, [2] *P. domestica* is a hexaploid of *P. spinosa* which is a discrete species, [3] *P. domestica* is an interspecific hybrid between a diploid *P. cerasifera* and a tetraploid *P. spinosa*, and [4] *P. domestica* is a hybrid of interspecific hybrids having a hexaploid chromosome complement composed of *P. cerasifera*, *P. spinosa*, and possibly contributions from other Eurasian plum species. The results presented here all but rule out models [1] and [2] but do not definitively resolve between [3] and [4]. Consistent with the prevailing theory proposed by early researchers, our data are most consistent with model [4] whereby one subgenome was contributed by diploid *P. cerasifera* and the other from a tetraploid *P. spinosa* that itself is an interspecific hybrid between diploid *P. cerasifera* and a second, as of yet, unknown *Prunus* species such as *Prunus ram burii* Boiss as suggested by Reales et al. in 2010⁶,²⁰,³². The endemic species *P. ramburii* was not included in our study but mixed grouping of *P. brigantina* and its hybrid

*P. brigantina* × *P. cerasifera* with *P. spinosa* and *P. cerasifera* accessions agrees with results reported by these authors and Shi et al. in 2013⁶,²⁰,³²,⁴⁹. Likewise, the lack of the DAP cp haplotype in the *P. cerasifera* or *P. spinosa* accessions tested imply that this haplotype may be derived from a distinct ancestral species. Irrespective of whether a 3rd species contributed to the *P. domestica* genome, the relatively close genetic distances among the germplasm supports the hypothesis by Eryomine, 1990 that some or all of the interspecific hybridization events that gave rise to *P. domestica* did not occur naturally but were artificially selected for by early Eurasian societies¹⁵. Such a scenario is consistent with the lack of known wild *P. domestica*, the low genetic diversity across the germplasm reflecting a more recent origin, and the low polli nation and fertility rates among experimental interspecific hybrids which would presumably have had tremendous difficulty naturally establishing themselves¹⁵,⁵⁰. In addi tion, historical evidence suggesting widespread simulta neous cultivation of *P. domestica*, *P. cerasifera* and *P. spinosa* is consistent with the current state of wild Eurasian germplasm that appears to be comprised by feral populations with relatively low genetic diversity. Whole genome sequencing, accompanied with additional geno typing using a broader range of Eurasian plum germplasm will be necessary to fully resolve these possibilities.

## Materials and methods

### Plant material

Plant material for this study included samples from 251 cultivars and forms of *P. domestica* held at the germplasm repository at INRA (Bordeaux, France), 66 cultivars from germplasm collections maintained at Uni versity of California at Davis and two USDA Agricultural Research Service (ARS) sites—National Clonal Germ plasm Repository (NCGR) and the Appalachian Fruit Research Station (AFRS, Kearneysville, WV). Also, 31 accessions representing related Eurasian plum species and cultivated varieties were obtained from Argentina, Aus tralia, Portugal, Sweden, and Ukraine. A complete list of accessions, their origin and pomological characteristics is provided in Supplemental Materials (Table S1). Classifi cation for pomological groups was conducted according to Neumüller¹² who separated plums, prunes, Reine Claudes (Greengages), Mirabelles, primitive forms and autochthonous forms. Finally, 34 genotypes derived from self-pollination of the cultivar “Improved French” main tained at AFRS (Kearneysville, WV) were used to aid in distinguishing self-pollinated genotypes from clonal material.

To verify the reproducibility of genotyping and dendo gram construction, we randomly selected and sequenced 23 of the DNA samples twice. Technical replicates were included in analyses under cultivar names followed by “.1”
and “.2”. Sixteen sets of biological replicates were also included, consisting of independently extracted DNAs from the same cultivar present in different locations, and listed under cultivar names followed by numbers 1 or 2 without a period.

## DNA extraction, library construction, and sequencing

Genomic DNAs were isolated using either the Qiagen DNeasy DNA extraction kit (Qiagen, Inc.) following manufacturer's instruction or a modified CTAB protocol by Kubisiak et al.⁵¹. The quality and integrity of the DNA was examined using a NanoDrop ND-8000 (Thermo Fisher Scientific Inc., USA) followed by electrophoresis on 1% agarose gels. Quantification of DNA samples was done using QuantiFluor® dsDNA kit (Promega, Inc.) and a Synergy H1 microplate reader (Biotek, USA).

Libraries for sequencing were constructed following protocols by Elshire et al.⁵² with a few modifications. Briefly, 100 ng of genomic DNAs were digested in the presence of *Pst*I (New England Biolabs, Inc.). Custom barcoded adapters were ligated to the *Pst*I-compatible overhang, and samples were pooled and then purified with a commercial PCR purification kit (Qiagen, Inc.), and finally quantified with a NanoDrop ND-1000. For test-amplification, 46 ng of DNA was amplified with common primers and NEBNext® High-Fidelity 2X PCR Master Mix (New England Biolabs, Inc.) as outlined in Elshire et al.⁵². The total volume of the PCR mix (50 µl) was divided into five tubes (10 µl each). Typically, 12, 14, 16, and 18 cycles were set up for testing against a non-amplified control. DNA fragments were resolved on 1.5% agarose gels against the Hi-Lo™ DNA marker (Bionexus, USA) to select optimal number of cycles for fragment enrichment in the range of 200–700 bp. Final libraries for genotyping-by-Sequencing (GBS) were amplified in two replicates (combined volume 100 µl) and purified using Mag-Bind® Total Pure NGS magnetic beads (Omega Bio-Tek, GA, USA). Two-step size selection with 0.4× and 0.8× volumes of magnetic beads depleted fragments larger than 1.5 kb and smaller than 121 bp, respectively. Quality checks for GBS libraries were done using a 2100 BioAnalyzer (Agilent Technologies, CA, USA). Fragment size distribution agreed with that predicted from “in silico” digestion of the *P. persica* reference genome with the *Pst*I restriction enzyme (not shown). Amplified GBS libraries were quantified using a Qubit® 2.0 fluorometer (ThermoFisher Scientific, Inc.) and were submitted for sequencing. Four 96-plexed and one 48-plexed sets of libraries were pair-end sequenced (2 × 125 bp) on a single lane each using the Illumina HiSeq 2500 instrument at the core facilities of the Medical University of South Carolina (Charleston, SC). Raw reads for all 405 accessions were deposited into NCBI's SRA database (BioProject PRJNA436025, SRA accession: SRP134093).

## Raw sequence data processing and genotyping

Data processing and SNP genotyping were performed using Stacks v.1.44-v.1.45³⁹,⁵³. Briefly, raw paired-end (PE) sequences were de-multiplexed according to barcodes, trimmed to remove adapters and low quality (<phred33) sequences, and filtered for the presence of the *Pst*I restriction sites. Typical success rate was in a range of 95–97%. Using the GSNAP software package⁵⁴ reads were aligned to the *P. persica* v 2.1 genome⁵⁵,⁵⁶. An inclusive catalog of tags and SNP genotypes (bi-allelic sites only) were generated using the “ref_map.pl” command with default parameters. The SNP genotypes were generated in the population mode of Stacks, recorded in the Variant Call Format (VCF) and filtered using VCFtools⁵⁷. To simplify analyses and limit errantly called variants, all indels and multi-nucleotide variants were deleted, leaving only bi-allelic SNPs and invariant sites. Individual genotypes were filtered to remove calls supported by less than 5 reads. Only SNPs present in more than 80% of all accessions (330 out of 405) were kept in the dataset for genetic analyses.

Structural and functional annotation of SNPs in coding and non-coding regions on the *P. persica* pseudochromosomes was done using the SnpEff 4.3e (build 2016-11-19) variant annotation and effect prediction tool⁵⁸. The physical positions (bp) of reference-based SNPs were correlated with the GFF file of the peach genome annotation v2.1 deposited at the Genome Database for Rosaceae (GDR)⁵⁶,⁵⁹. The genomic distributions of SNPs across all eight peach chromosomes (pseudomolecules) were plotted individually based on their physical positions (100K screening window) and visualized using the Circos visualization tool⁶⁰.

The same analytical pipeline was used for isolation of genomic fragments derived from the peach chloroplast (cp) genome. Demultiplexed reads were aligned against the *P. persica* cp sequence deposited at GenBank (HQ336405.1). The cp-derived polymorphic sites were called using “ref_map” command adjusted for ploidy level 1³⁹,⁵³.

## Genetic diversity and population structure analysis

The SNPRelate package in R⁶¹ was used to investigate genetic diversity of the plum germplasm. A distance matrix derived from the Euclidian metric was calculated and used in hierarchical cluster analysis. A dendrogram was constructed using the between-group average linkage method, or UPGMA. The confidence interval of the genetic relationships among the accessions was determined by performing 5000 bootstrap replicates. A graphical representation of the resulting tree (in Newick format) was visualized in Geneious R10³⁵. Wright's fixation index ($F_{ST}$) was used to explore the degree of differentiation between pomological groups and species.
$F_{ST}$ values and standard errors for all pairwise comparisons were calculated using a “population” mode in Stacks³⁹.

Principal component analysis (PCA) was performed using the R package SNPRelate⁶¹. The VCF file with genotypic data was converted into the genomic data structure (GDS) file format, filtered for MAF ≥ 0.05 and missing rate ≤0.20 and used for PCA analysis with the snpgdsPCA function. The percentage of explained variation was calculated for the first six PCs used for plotting the genotypes on a two-dimensional scale. All plots were generated using ggplot2 included with the library.

The nonparametric tree-based ensemble random forest (RF) analysis⁴⁴ was executed using the R package Random Forest with 1000 permutation trees⁶². The variable (i.e., SNP) importance was estimated based on the Mean Decrease Accuracy (MDA) metric that measures contribution of individual SNPs to an accuracy of prediction when all other SNPs are constant. Higher MDA values indicate large changes to the prediction power, and hence, the importance of the SNP. In total, 7 output datasets of the most significant SNPs were generated with the “ntree = ” parameter settings of 250, 500, 1000, 2000, 5000, 8000, and 15,000. The distribution of accessions on a PCA plot (with most significant SNPs) and RF plot was compared using the Procrustes technique implemented in the R package vegan 2.4-4¹<fnref n="1" />. Association between nuclear markers and cp haplotypes was estimated by the CCA using vegan in R. Significance of constrain was estimated with the *anova.cca* function by permutation test under reduced model (n = 999).

## Acknowledgements

This work was supported in part by a National Institute of General Medical Sciences Grant: P20GM109094 and funding from the California Dried Plum Board. Also, we are thankful to Joe Turkovich for coordinating commercial sample collections and Drs. Rita Costa (INIAV, Portugal), Larisa Gustavsson-Garkava (SLU, Sweden), Craig Hardner (Australia), and Humberto Prieto (Chile) for providing DNA or leaf material that made this study possible. The authors would like to thank Pascal Briard (INRA-Bordeaux) for sampling leaves and contributing to DNA extraction, Bernard Lafargue (INRA-Bordeaux) for providing valuable information about the sources of *Prunus* collections and Marine Blouin in charge of the Centre de Ressources Genetiques, INRA-CRG-Prunus, Domaine de la Tour de Rance, 47320 Bourran.

## Author details

## Conflict of interest

The authors declare that they have no conflict of interest.

## Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Supplementary Information** accompanies this paper at (https://doi.org/10.1038/s41438-018-0090-6).

Received: 30 May 2018 Revised: 24 August 2018 Accepted: 12 September 2018

Published online: 01 January 2019

## References

1. Stacewicz-Sapuntzakis, M. Dried plums and their products: composition and health effects-an updated review. *Crit. Rev. Food Sci. Nutr.* **53**, 1277-1302 (2013).

2. Wallace, T. C. Dried plums, prunes and bone health: a comprehensive review. *Nutrients* **9**, E401 (2017).

3. Faust, M. & Surányi, D. Origin and dissemination of plums. *Hortic. Rev.* **23**, 179-231 (1999).

4. Hedrick, U. P. The Plums of New York. pp 844 (J. B. Lyon Co., State Printers, Albany, 1911).

5. Vávilov, N. Wild progenitors of fruit trees and shrubs in Turkestan and the Caucasus and the problem of the origin of fruit trees. In *Ninth International Horticultural Congress Report*, London (1930).

6. Rybin, W. A. Spontaneous and experimentally produced hybrids between blackthorn and cherry plum and the descent problem of the cultivated plum. *Planta* **25**, 22-58 (1936).

7. Lin, P. & Shi, L. The discovery and distribution of Ili wild *Prunus domestica* (*P. communis* Fritsch) in Xinjiang, pp. 282-286. In *Proceedings of International Symposium on Horticultural Germplasm, Cultivated and Wild*. Chinese Society for Horticultural Science, Beijing, 5-9 Sept. 1988 (1990).

8. Lin, S. H., Pu, F. S., Zhang, J. H., Gao, X. Y. & Li, X. J. Observation on chromosomes numbers of *Prunus*. *China Fruits* **2**, 8-10 (1991).

9. Liu, W. et al. Genetic diversity and phylogenetic relationships among plum germplasm resources in China assessed with inter-simple sequence repeat markers. *J. Am. Soc. Hort. Sci.* **132**, 619-628 (2007).

10. Woldring, H. On the origin of plums: a study of sloe, damson, cherry plum, domestic plums and their intermediates. *Palaeohistoria* **3**, 535-562 (2000).

11. Zohary, D. & Hopf, M. *Domestication of plants in the Old World*. 3rd edn, 316 pp. (Oxford University Press, New York, 2000).

12. Neumüller, M. Fundamental and applied aspects of plum (*Prunus domestica*) breeding. *Fruit*, *Veg. Cereal Sci. Biotechnol.* **5**, 139-156 (2011).

13. Topp, B. L., Russell, D. M., Neumüller, M., Dalbó, M. A., & Liu, W. in *Fruit Breeding, Vol. 8 of Handbook of Plant Breeding* (eds Badenes, M. L. & Byrne, D. H.) pp. 571-621 (Springer, New York, NY, 2012).

14. Chin, S. W., Shaw, J., Haberle, R., Wen, J. & Potter, D. Diversity of almonds, peaches, plums and cherries-molecular systematics and biogeographic history of *Prunus* (*Rosaceae*). *Mol. Phylogenet. Evol.* **76**, 34-48 (2014).

15. Eryomine, G. V. New data on origin of *Prunus domestica* L. *Acta Hortic.* **283**, 27-30 (1990).

16. Reynders, S. & Salesses, G. Study on the relationships within the subgenus Prunophora. Restriction maps of the ribosomal genes in *P. cerasifera* and *P. spinosa*. *Acta Hortic.* **283**, 27-29 (1991).

17. Reynders-Aloisi, S. & Grellet, E. Characterization of the ribosomal DNA units in two related *Prunus* species (*P. Cerasifera* and *P. Spinosa*). *Plant Cell Rep.* **13**, 641-646 (1994).

18. Bortiri, E. S. H. et al. Phylogeny and systematics of *Prunus* (*Rosaceae*) as determined by sequence analysis of ITS and the chloroplast trnL-trnF spacer DNA. *Syst. Bot.* **26**, 797-807 (2001).

19. Potter, D. in *Wild Crop Relatives: Genomic and Breeding Resources* (ed. Kole, C.) pp. 129-145 (Springer, Berlin, Heidelberg, 2011).

20. Crane, M. & Lawrence, W. Genetics of Garden Plants p. 236 (The MacMillan Co, London, 1934).

21. Horvath, A. et al. Phenotypic variability and genetic structure in plum (*Prunus domestica* L.), cherry plum (*P. cerasifera* Ehrh.) and sloe (*P. spinosa* L.). *Sci. Hortic.* **129**, 283-293 (2011).

22. Athanasiadis, I., Nikoloudakis, N. & Hagidimitrion, M. Genetic relatedness among cultivars of the Greek plum germplasm. *Not. Bot. Horti. Agrobo*. **41**, 491-498 (2013).

<footnote n="1">¹The Schatz Center for Tree Molecular Genetics, Department of Ecosystem Sciences and Management, The Pennsylvania State University, University Park, PA, 16802, USA. ²Genomics and Computational Biology Laboratory, Clemson University, Clemson, SC 29634, USA. ³USDA Appalachian Fruit Research Laboratory, Kearneysville, WV 25430, USA. ⁴UMR BFP1332 - INRA-Bordeaux, Bordeaux II University, 33882 Villenave d'Ornon, France. ⁵Department of Plant Sciences, University of California, Davis, CA 95616, USA</footnote>
23. Kazija, D. H. et al. Plum germplasm in Croatia and neighbouring countries assessed by microsatellites and DUS descriptors. *Tree. Genet. Genom.* **10**, 761–778 (2014).

24. Sehic, J., Nybom, H., Hjeltnes, S. H. & Gaši, F. Genetic diversity and structure of Nordic plum germplasm preserved ex situ and on-farm. *Sci. Hortic.* **190**, 195–202 (2015).

25. Öz, M. H. et al. Molecular analysis of East Anatolian traditional plum and cherry accessions using SSR markers. *Genet. Mol. Res.* **12**, 5310–5320 (2013).

26. Sun, Q. et al. Analysis of genetic relationships in *Prunus domestica* L. in Xinjiang Using ISSR markers. *Agric. Sci. Technol.* **16**, 447–453 (2015).

27. Decroocq, V., Hagen, L. S., Favé, M. G., Eyquard, J. P. & Pierronnet, A. Micro-satellite markers in the hexaploid *Prunus domestica* species and parentage lineage of three European plum cultivars using nuclear and chloroplast simple-sequence repeats. *Mol. Breed.* **13**, 135–142 (2004).

28. Gharbi, O., Wunsch, A. & Rodrigo, J. Characterization of accessions of 'Reine Claude Verte' plum using Prunus SRR and phenotypic traits. *Sci. Hortic.* **169**, 57–65 (2014).

29. Kota-Dombrovska, I. & Lacis, G. Evaluation of self-incompatibility locus diversity of domestic plum (*Prunus domestica* L.) using DNA-based S-genotyping. *Proc. Latvian Acad. Sci.* **67**, 109–115 (2013).

30. Halász, J., Kurilla, A. & Hegedűs, A. Preliminary characterization of the self-incompatibility genotypes of European plum (*Prunus domestica* L.) cultivars. *International. J. Hortic. Sci.* **20**, 23–26 (2014).

31. Halász, J., Makovics-Zsohár, N., Szőke, F., Ercisli, S. & Hegedűs, A. Simple sequence repeat and S-locus genotyping to explore genetic variability in polyploid *Prunus spinosa* and *P. insititia*. *Biochem. Genet.* **55**, 1–12 (2017).

32. Reales, A., Sargent, D. J., Tobutt, K. R. & Rivera, D. Phylogenetics of Eurasian plums, *Prunus* L. section *Prunus* (Rosaceae), according to coding and non-coding chloroplast DNA sequences. *Tree Genet. Genomes* **6**, 37–45 (2010).

33. Clevenger, J., Chavarro, C., Pearl, S. A., Ozias-Akins, P. & Jackson, S. A. Single nucleotide polymorphism identification in polyploids: a review, example, and recommendations. *Mol. Plant* **8**, 831–846 (2015).

34. Andrews, K. R, Good, J. M., Miller, M. R., Luikart, G. & Hohenlohe, P. A. Harnessing the power of RADseq for ecological and evolutionary genomics. *Nat. Rev. Genet.* **17**, 81–92 (2016).

35. Kearse, M. et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. *Bioinformatics* **28**, 1647–1649 (2012).

36. Németh, M. History and importance of plum pox in stone-fruit production. *EPBO Bull.* **24**, 525–536 (1994).

37. Doyle, J. F., DeBuse, C. J. & DeJong, T. M. Varieties. in Prune Production Manual, Vol. 3507 (ed. Buchner, R.) pp. 75–92 (University of California Agriculture and Natural Resources Publication, 2012).

38. Britten, J. Greengages. *Ir. Mon.* **17**, 635–639 (1889).

39. Catchen, J., Amores, A., Hohenlohe, P., Cresko, W. & Postlethwait, J. Stacks: building and genotyping loci de novo from short-read sequences. *G3: Genes Genom. Genet.* **1**, 171–182 (2011).

40. Hartl, D. L. & Clark, G. C. Principles of PopulationGenetics. 4th edn, pp. 95–105 (Sinauer Associates, Inc. Publishers Sunderland, Massachusetts, 1997).

41. Oksanen, J. Vegan: an introduction to ordination. Available at https://github.com/vegandevs/vegan/vignettes (2008).

42. Schmidt, M. Contributions to breeding research in plums. *Breeder* **24**, 157–161 (1954).

43. Johansson, J. & Oldén, E. J. in *Handbook of Plant Breeding, Vol. 6: Breeding of Vegetables, Fruits, Vines and Forest Plants* (eds Kappert, H. & Rudorf, W.) pp. 602–624 (Paul Parey Verlag, Berlin, 1962).

44. Breiman, L. Random forests. *Mach. Learn.* **45**, 5–32 (2001).

45. Sylvester, E. V. A. et al. Applications of random forest feature selection for fine-scale genetic population assignment. *Evol. Appl.* **11**, 153–165 (2018).

46. Bertolini, F. et al. Combined use of principal component analysis and Random Forests identify population-informative single nucleotide polymorphisms: application in cattle breeds. *J. Anim. Breed. Genet.* **132**, 346–356 (2015).

47. Byrne, S. L. et al. Using variable importance measures to identify a small set of SNPs to predict heading date in perennial ryegrass. *Sci. Rep.* **7**, 3566 (2017).

48. Yang, X. et al. Mining sequence variations in representative polyploid sugarcane germplasm accessions. *BMC Genom.* **18**, 594 (2017).

49. Shi, S., Li, J., Sun, J., Yu, J. & Zhou, S. Phylogeny and classification of *Prunus sensu lato* (*Rosaceae*). *J. Integr. Plant Biol.* **55**, 1069–1079 (2013).

50. Oldén, E. J. Interspecific plum crosses. Balsgard Fruit Breeding Inst. Fjalkestad. *Sweden* **1**, 57 (1965).

51. Kubisiak, T. L. et al. A transcriptome-based genetic map of Chinese chestnut (*Castanea mollissima*) and identification of segmental homology with peach (*Prunus persica*). *Tree Genet. Genomes* **9**, 557–571 (2013).

52. Elshire, R. J. et al. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. *PLoS ONE* **6**, e19379 (2011).

53. Rochette, N. & Catchen, J. Deriving genotypes from RAD-seq short-read data using Stacks. *Nat. Prot.* **12**, 2640–2659 (2017).

54. Wu, T. D. & Nacu, S. Fast and SNP-tolerant detection of complex variants and splicing in short reads. *Bioinformatics* **26**, 873–881 (2010).

55. Verde, I. et al. The high-quality draft genome of peach (*Prunus persica*) identifies unique patterns of genetic diversity, domestication and genome evolution. *Nat. Genet.* **45**, 487–494 (2013).

56. Verde, I. et al. The peach v2.0 release: high-resolution linkage mapping and deep resequencing improve chromosome-scale assembly and contiguity. *BMC Genom.* **18**, 225 (2017).

57. Danecek, P. et al. 1000 Genomes Project Analysis Group: The variant call format and VCFtools. *Bioinformatics* **27**, 2156–2158 (2011).

58. Cingolani, P. et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of *Drosophila melanogaster* strainw1118; iso-2; iso-3. *Fly* **6**, 80–92 (2012).

59. Jung, S. et al. The Genome Database for Rosaceae (GDR): year 10 update. *Nucleic Acids Res.* **42**, D1237–D1244 (2014).

60. Krzywinski, M. et al. Circos: an information aesthetic for comparative genomics. *Genome Res.* **19**, 1639–1645 (2009).

61. Zheng, X. et al. A high-performance computing toolset for relatedness and principal component analysis of SNP data. *Bioinformatics* **28**, 3326–3328 (2012).

62. Liaw, A. et al. Classification and regression by random forest. *R. News* **2**, 18–22 (2002).