# Genetic diversity and structure of Slovenian native germplasm of plum species (P. domestica L., P. cerasifera Ehrh. and P. spinosa L.)

T. Ternjak, Teresa Barreneche, M. Šiško, A. Ivančič, A. Šušek, José

Quero-Garcia

► To cite this version:

T. Ternjak, Teresa Barreneche, M. Šiško, A. Ivančič, A. Šušek, et al., Genetic diversity and structure of Slovenian native germplasm of plum species (P. domestica L., P. cerasifera Ehrh. and P. spinosa L.). Frontiers in Plant Science, 2023, 14, pp.1150459. 10.3389/fpls.2023.1150459. [https://hal.inrae.fr/hal-04082299](https://hal.inrae.fr/hal-04082299)

HAL Id: hal-04082299

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

Submitted on 26 Apr 2023

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
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN ACCESS

EDITED BY
Xiaohua Qi,
Yangzhou University, China

REVIEWED BY
Monica Marilena Miazzi,
University of Bari Aldo Moro, Italy
Annalisa Marchese,
University of Palermo, Italy

*CORRESPONDENCE

J. Quero-García
[jose.quero-garcia@inrae.fr](mailto:jose.quero-garcia@inrae.fr)
A. Šušek
[andrej.susek@um.si](mailto:andrej.susek@um.si)

SPECIALTY SECTION

This article was submitted to
Plant Breeding,
a section of the journal
Frontiers in Plant Science

RECEIVED 24 January 2023

ACCEPTED 09 March 2023

PUBLISHED 21 March 2023

CITATION

Ternjak T, Barreneche T, Šiško M, Ivančič A, Šušek A and Quero-García J (2023) Genetic diversity and structure of Slovenian native germplasm of plum species (*P. domestica* L., *P. cerasifera* Ehrh. and *P. spinosa* L.). *Front. Plant Sci.* 14:1150459. doi: 10.3389/fpls.2023.1150459

COPYRIGHT

© 2023 Ternjak, Barreneche, Šiško, Ivančič, Šušek and Quero-García. This is an open-access article distributed under the terms of the [Creative Commons Attribution License (CC BY)](https://creativecommons.org/licenses/by/4.0/). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

# Genetic diversity and structure of Slovenian native germplasm of plum species (*P. domestica* L., *P. cerasifera* Ehrh. and *P. spinosa* L.)

T. Ternjak¹, T. Barreneche², M. Šiško¹, A. Ivančič¹, A. Šušek¹* and J. Quero-García²*

¹Faculty of Agriculture and Life Sciences, University of Maribor, Hoče, Slovenia, ²INRAE, Univ. Bordeaux, UMR BFP, Villenave d'Ornon, France

Slovenia has particular climatic, soil, geographic and historical conditions that lead to long tradition of plum cultivation and use. In this work, a set of 11 SSR and three universal cpDNA markers, as well as flow cytometry, were used to (1) evaluate the genetic diversity of 124 accessions of the three *Prunus* species (*P. domestica* L., *P. cerasifera* Ehrh., and *P. spinosa* L.), (2) investigate the possible involvement of *P. cerasifera* and *P. spinosa* species in *P. domestica* origin, (3) study the genetic relationships and variability among the most typical *P. domestica* accessions present in Slovenia. Ten haplotypes of cpDNA were identified and clustered into three groups according to the Neighbor-Joining analysis (NJ). All 11 SSR primer pairs were polymorphic, revealing 116 unique genotypes. A total of 328 alleles were detected with an average value of 29.82 alleles per locus, showing relatively high diversity. Bayesian analysis of genetic structure was used to identify two ancestral populations in the analyses of all three species as well as in a separate set consisting of *P. domestica* material only. Principal Coordinate Analysis (PCoA) showed that accessions clustered largely in agreement with Bayesian analysis. Neighbor-Joining analysis grouped 71 *P. domestica* accessions into three clusters with many subgroups that exhibited complex arrangement. Most accessions clustered in agreement with traditional pomological groups, such as common prunes, mirabelle plums and greengages. In this study, the analyses revealed within *P. domestica* pool valuable local landraces, such as traditional prunes or bluish plums, which seem to be highly interesting from a genetic point of view. Moreover, complementary approaches allowed us to distinguish between the three species and to gain insights into the origin of plum. The results will be instrumental in understanding the diversity of Slovenian plum germplasm, improving the conservation process, recovering local genotypes and enriching existing collections of plant genetic resources.

KEYWORDS

*Prunus* spp., plum, genetic resources, genetic diversity, genetic structure, cpDNA, SSR
# 1 Introduction

Plums are considered as one of the most colourful and delicious temperate fruits, with a wide range of fruit sizes and shapes (Figure S1A), and are well adapted to different growing conditions (Milošević and Milošević, 2018). They are deciduous fruit trees belonging to the stone fruits classified in the Rosaceae family and genus *Prunus*. It is thought that the main existing European plum species, *P. domestica* L. (European plum) (Figure S1B), *P. cerasifera* Ehrh. (Cherry plum, myrobolan) (Figure S1C) and *P. spinosa* L. (blackthorn, sloe) (Figure S1D) descended from an ancestor that migrated from East Asia (Chin et al., 2014).

According to Insero et al. (1991), Eurasian plum is one of the most important and globally distributed *Prunus* group, which includes the hexaploid European plum, *P. domestica* (2n = 6x = 48), and *P. insititia*, which was considered a distinct species by some authors in the past and is now classified as a subspecies of *Prunus domestica* (*P. domestica* subsp. *insititia* C. K. Schneider). European plum cultivars have been generally classified into different pomological groups based on morphological characteristics as well as fruit use, including large-fruited European plums, prunes, egg plums, greengages, mirabelles, damsons, bullaces and St. Julien plums (Zhebentyayeva et al., 2019). Prunes are mainly processed for drying, but can also be used fresh, for distilling, baking and jams. They differ from the other plum groups, as the whole fruit can be dried, without fermenting, due to their high sugar content (Faust and Suranyi, 1999). Traditionally, egg plums produce yellow fruits and are mainly used for canning, while greengages usually produce round, small, green fruits that are considered the sweetest among plums (Topp et al., 2012). Damsons, bullaces, mirabelles and St. Julien plums are mostly used for food processing. St. Julien plums are also commonly used as rootstocks. Very sweet and highly aromatic mirabelle plums have yellow to orange fruits, while damsons and bullaces form yellow or greenish to bluish, small, round to oval fruits (Gaši et al., 2020).

*P. cerasifera* (2n = 2x = 16), the diploid cherry plum, contains wild and cultivated genotypes used as ornamentals, for their fruits, or as rootstocks. Another important species closely related to the cultivated plum is the tetraploid *P. spinosa* (2n = 4x = 32). Some authors consider this wild species as one of the possible ancestors of *P. domestica*. There are still numerous questions about the genetic origin of *P. domestica*. Most authors support the theory that hexaploid *P. domestica* is an interspecific hybrid between diploid *P. cerasifera* and tetraploid *P. spinosa* (Rybin, 1936; Crane and Lawrence, 1938; Werneck, 1961; Decroocq et al., 2004; Bortiri et al., 2009; Ilgin et al., 2014). In contrast, other authors discuss that European plum could have been a result of alloploidization of polyploid forms of *P. cerasifera* (Bajashvili, 1990; Zohary, 1992; Decroocq et al., 2004; Okie and Hancock, 2008; Depypere et al., 2009; Weiss et al., 2012) or *P. spinosa* (Bajashvili, 1990; Woldring, 2003; Liu et al., 2007). Hedrick (1911) hypothesised that *P. domestica* was derived from *P. spinosa*, with *P. insititia* being an intermediate form between the two. More recent results support the possibility that *P. domestica* arose from interspecific hybrids with a hexaploid chromosome complex of *P. cerasifera*, *P. spinosa*, and

potential progenitor species from the Eurasian plum cluster (Zhebentyayeva et al., 2019).

Slovenia, part of the former Republic of Yugoslavia, located in the northwest of the Balkan Peninsula, has a long history of cultivating plum species. It is believed that the Balkan Peninsula is one of the richest regions in Europe in terms of plum fruits, with highly valuable genetic material (Vukojevic et al., 2012). Moreover, according to Vavilov (1949), it belongs to one of the eight independent centres of origin of the most important cultivated plants in the world: the Mediterranean Centre. Mihai et al. (2012) considered *P. domestica* as an indigenous species in the Balkans, represented by thousands of local biotypes that have been cultivated in this area for centuries and are part of the local culture.

On the territory of present-day Slovenia, the first systematic plum cultivation took place during the Roman period (Kazija et al., 2014). In the Middle Ages, plums were grown mainly in village gardens, monasteries and small orchards (Rencelj and Pozrl, 2015). In the nineteenth century, plum was the most prevalent fruit tree species, with the leading *P. domestica* landrace, the common prune ('Domača češplja', syn. 'Bistrica', 'Požegača', 'Dolanka'), also known in the wider area of the Balkan Peninsula. Later, other varieties of European plum spread, and in the last two decades the selection was also enriched by some varieties of Chinese-Japanese origin (Štampar et al., 2005).

Different plum varieties were used both, for fresh consumption and for processing. In the eastern (Styria and Lower Carniola) and western (Goriška brda, Soča and Vipava valleys) parts of the country, plum drying was seasonal, but a very important source of income for rural population, especially women from small farm households (Pockar, 1982). The tradition of peeling and drying plums, called 'prunes' ('prunele'), was brought to Slovenia by French soldiers that deserted Napoleon's army at the beginning of the 19th century (Pockar, 1982). The prevailing genotype used for plums drying was landrace 'Domača češplja'. Distillation of plum brandy (slivovec, slivovica) was also an important source of income for rural folks and usually made with common prune 'Domača češplja' (Rencelj and Pozrl, 2015).

In Styria (NE Slovenia), the sweet and tasty fruits of another important landrace were used. The so-called 'plavkica' (directly translated - small bluish fruits) is an early flowering, round, small fruited, dark blue or violet-blue plum, with greenish to yellowish colour of the flesh that fully or partly adheres to the stone. Written information about this landrace is very scarce. Mišić (1979) stated that 'plavkica' (from this point referred as Bluish plum) was one of the widespread varieties grown in Slovenia, while Honzak (1968) wrote about seedlings of cultivated local varieties used as a rootstocks.

Highly specific climatic, soil, geographic and historical conditions of Slovenia resulted in valuable plum diversity and many unique genotypes could harbour important characteristics for subsequent breeding and food security. Unfortunately, the use of local plum varieties and the number of varieties in cultivation have declined significantly in the recent decades (Usenik et al., 2007). Environmental changes, land use and intensified cultivation practises, introduction of new cultivars, as well as economic and
social changes, have led to significant erosion of indigenous genetic diversity (Benediková and Giovannini, 2011; Luthar, 2018). In addition, *Plum pox virus* disease (Sharka) has greatly affected plum cultivation as well as the distribution of the species in its natural habitat, resulting in plum germplasm depletion (Usenik et al., 2007; Kazija et al., 2014). Climate change is one of the biggest challenges facing food production. Therefore, more sustainable agricultural practises need to be adopted, which go hand in hand with breeding strategies. In order to create new cultivars, well adapted to these challenges, exploring a broader genetic base, including primitive indigenous germplasm, is essential to the success of current breeding programmes (Campoy et al., 2016).

To develop more sustainable solutions for the future it has never been more important, to systematically collect, record, evaluate, and conserve the most valuable genetic material (Vukojevic et al., 2012). In the past, studies on the genetic diversity of *P. domestica* were quite scarce due to the polyploid nature of the species, but recently plum germplasms from different countries have been evaluated (Kazija et al., 2014; Sehic et al., 2015; Urbanovich et al., 2017; Urrestarazu et al., 2018; Manco et al., 2019; Zhebentyayeva et al., 2019). Genotyping of autochthonous genetic resources allows identification of synonyms and homonyms, determination of parentage/origin, and identification/evaluation of important genotypes (Gaši et al., 2020). Additional studies, based on agro-morphological and molecular data are needed to assess the current status of diversity and genetic potential of resources and to prepare effective conservation management as well as increase their use for research, exchange, breeding and cultivation (Kaufmane et al., 2002; Blazek, 2007).

The present research was set up to investigate the status of traditional plum germplasm as well as the genetic relationships and variability among accessions belonging to three *Prunus* species (*P. domestica*, *P. cerasifera* and *P. spinosa*) representing the Slovenian plum gene pool. Based on historical and economic importance, the most frequently present traditional plum genotypes in Slovenia probably belong to the species *P. domestica*. Therefore, our work focused on this species with particular interest on Bluish plum and common prune. In the present study, we used a set of SSR and cpDNA markers to: 1) evaluate the relationship among the three *Prunus* species (*P. domestica*, *P. cerasifera* and *P. spinosa*) 2) examine the possible involvement of *P. cerasifera* and *P. spinosa* species in the origin of *P. domestica* 3) investigate the genetic relationships and structure among *P. domestica* accessions collected *in situ* in the Slovenian landscape.

# 2 Materials and methods

## 2.1 Plant material

The study included 124 accessions (Table S1). Among these, 54 accessions were collected *in situ*: in the wild, along roads, in home gardens, on abandoned estates and in orchards. The material was collected mostly throughout Slovenia with the emphasis on the north-eastern part, which is considered one of the typical plum-growing regions in Slovenia (Pockar, 1982), while two samples were found in Styria region of Austria and one accession was brought

from Croatia. In addition, 70 samples were obtained from various *ex situ* collections: 38 from the Slovene Plant Gene Bank (SPGB) of the Faculty of Agriculture and Life Sciences (Maribor), 21 from the Plum French National Collection (PFNC) of the 'Centre de Ressources Biologiques - CRB *Prunus Juglans*', French National Research Institute for Agriculture, Food and Environment, (INRAE) (Nouvelle Aquitaine Bordeaux), and eleven from different botanical gardens or private collections from Slovenia (Maribor and Ljubljana) and Austria (Vienna and Graz). Of the 124 studied accessions, a set of 14 internationally known cultivars was used as reference. The accessions from the PFNC Collection were added in order to compare the genetic resources of Slovenian plum with a sample representing a large diversity within the species. This collection has been well described morphologically in INRAE multispecies integrative information system dedicated to plant and fungi pests (GnpIS, https://urgi.versailles.inra.fr/gnpis), as well as studied with molecular markers by different authors (Horvath et al., 2011; Urrestarazu et al., 2018; Zhebentyayeva et al., 2019). Prior to analysis, all accessions were classified into three different species groups (*P. domestica*, *P. cerasifera* or *P. spinosa*) based on data from the existing collections and our morphological descriptions.

The studied material was also classified according to its status, e.g., wild genotype, landrace, or improved variety. Based on passport data from the existing *ex situ* collections and according to Campoy et al. (2016) and the bibliography, improved varieties were considered "early selections" (such as 'President') or modern breeding cultivars (such as 'Topper'), depending on the period when they were selected.

## 2.2 DNA extraction

Total genomic DNA was isolated following the CTAB protocol described by Doyle and Doyle (1987) with some modifications: DNA was precipitated by adding 60 µl sodium acetate, and samples were incubated at -20°C for 20 min without addition of resuspension buffer or RNAse. DNA concentration was quantified using a fluorimeter (Hoefer DQ 300, Holliston, USA), and quality was checked by electrophoresis on a 3% agarose gel. Products were visualised under UV light (302 nm), using a Benchtop UVP UV Transilluminator (Analytik Jena, Jena, Germany).

## 2.3 Assessment of ploidy level

Relative DNA content was estimated by flow cytometry using DAPI (4',6'-diamidino-2-phenylindole) staining in Otto buffers. Samples were prepared according to Bohanec (2003) using *Trifolium repens* L. as an internal standard with a 2C nuclear DNA content of 2.07 pg (Arumuganathan and Earle, 1991). Ploidy analysis was performed with a CyFlow® Space Flow Cytometer (Sysmex Partec, Goerlitz, Germany) using a linear scale. Histograms were analysed using Flomax software 2.0 (Sysmex Partec, Goerlitz, Germany). The absolute amount of DNA in a sample was calculated according to the protocols proposed by Dolezel (1991), Marie and Brown (1993) and Dolezel and Bartos (2005).
The ploidy level was already determined for the accessions from the PFNC Collection (INRAE, Horvath, unpublished data), accordingly some of the material was used as reference [e.g.: 170 'D'Ente double' (2n = 6x), 181 blackthorn (2n = 4x), 187 myrobolan 'Agdzadzor' (2n = 2x)]. For each reference accession, the absolute amount of DNA was calculated and the values were compared between all accessions.

## 2.4 PCR-RFLP analyses

Three universal primers were used for chloroplast DNA (cpDNA) analyses [i.e. HK, K1K2 (Demesure et al., 1995) and VL (Dumolin-Lapegue et al., 1997)]. These primers were previously used in *Prunus* sp. (Bouhadida et al., 2007; Horvath et al., 2008; Horvath et al., 2011).

Amplified fragments were digested with three primer pair-restriction enzyme combinations for K1K2 (*Hinf I*, *Taq I*, *Alu I*) and with two for VL and HK (*Hinf I*, *Taq I*). Analyses were performed by BioGEVES platform (GEVES, Beaucouzé, France) following the protocol described by Horvath et al. (2011). The size of each fragment was evaluated and scored as binary data (presence/absence). Bands were coded as unique combinations of letters (A, B, C, D) and assigned to a specific haplotype. Restriction fragment data were transformed into a binary matrix and the Unweighted Neighbor-Joining dendrogram was created using DARwin 6.0.21 software (Perrier and Jacquemoud-Collet, 2006).

## 2.5 SSR analyses

All individuals were analysed with a set of 11 SSR primer pairs (Table 1) developed on different *Prunus* species (peach, sweet cherry and Japanese plum) (Cipriani et al., 1999; Wünsch, 2009; Mnejja et al., 2010). In addition to their position on the *Prunus* reference map, markers were also selected based on their amplification quality and polymorphism, as reported in previously published works. Primers were labelled with Cy5 or Cy5.5 dye according to Sigma-Aldrich (Saint Louis, USA). Polymerase chain reactions (PCR) were performed using HotStarTaq® Plus Master Mix Kit (Qiagen, Hilden, Germany), and the annealing temperature was adjusted accordingly for each marker (approximately 5°C below the Tm). Amplifications were performed in a Biometra TProfessional thermal cycler (Analytik Jena, Jena, Germany), and the quality of the fluorescently-labelled PCR products was checked. Fragment size analysis was performed using a Beckman CEQ 8000 DNA sequencer (Beckman Coulter Inc., Brea, USA) according to the manufacturer's instructions. Alleles appear as peaks in CEQ 8000 DNA sequencer and are compared to the standards during scoring to determine their size. Allele fragments were scored for each accession. First, we started with *P. cerasifera* accessions since they are diploids and therefore easier to evaluate. Then, each SSR was characterized by analyzing the peak profiles for all accessions. Following this approach, the accessions of the other two species were evaluated considering the corresponding profile observed in *P. cerasifera*. When *P. cerasifera* accessions had more than two peaks within the theoretical range of

<table><caption>TABLE 1 SSR loci, linkage group (LG), repeat motifs and references of the 11 primer pairs used in this study.</caption><thead><tr><td>Locus</td><td>LG</td><td>Source species</td><td>Primer sequence (5'-3')</td><td>Repeat motif</td><td>Reference</td></tr></thead><tbody><tr><td>UDP96-005</td><td>1</td><td><i>P. persica</i> (L.) Batsch</td><td>For: GTAACGCTCGCTACCACAAA<br>Rev: CCTGCATATCACCACCCAG</td><td>(AC)<sub>16</sub>TG(CT)<sub>2</sub>CA(CT)<sub>11</sub></td><td>Cipriani et al. (1999)</td></tr><tr><td>BPPCT034</td><td>2</td><td><i>P. persica</i> (L.) Batsch</td><td>For: CTACCTGAAATAAGCAGAGCCAT<br>Rev: CAATGGAGAATGGGGTGC</td><td>(GA)<sub>19</sub></td><td>Dirlewanger et al. (2002)</td></tr><tr><td>EMPAS12</td><td>3</td><td><i>P. avium</i> L.</td><td>For: TGTGCTAATGCCAAAAATACC<br>Rev: ACATGCATTTCAACCCACTC</td><td>(TG)<sub>10</sub>A/GA)<sub>10</sub>AA(GA)<sub>13</sub></td><td>Vughan and Russell (2004)</td></tr><tr><td>UCD-CH17</td><td>4</td><td><i>P. avium</i> L.</td><td>For: TGGACTTCACTCATTTTCAGAGA<br>Rev: ACTGCAGAGAATTTCCACAACCA</td><td>(CT)<sub>11</sub></td><td>Struss et al. (2003)</td></tr><tr><td>EMPAS06</td><td>4</td><td><i>P. avium</i> L.</td><td>For: AAGCGGAAAGCACAGGTAG<br>Rev: TTGCTAGCATAGAAAAGAATTGTAG</td><td>(CT)<sub>12</sub></td><td>Vughan and Russell (2004)</td></tr><tr><td>EMPAS11</td><td>5</td><td><i>P. avium</i> L.</td><td>For: ACCACTTTGAGGAACTTGGG<br>Rev: CTGCCTGGAAGAGCAATAAC</td><td>(TC)<sub>25</sub></td><td>Vughan and Russell (2004)</td></tr><tr><td>EMPAS14</td><td>5</td><td><i>P. avium</i> L.</td><td>For: TCCGCCATATCACAATCAAC<br>Rev: TTCCACACAAAAAACCAATCC</td><td>(TC)<sub>10</sub>CCAT(TC)<sub>5</sub>CCAT(TC)<sub>8</sub></td><td>Vughan and Russell (2004)</td></tr><tr><td>BPPCT014</td><td>5</td><td><i>P. persica</i> (L.) Batsch</td><td>For: TTGTCTGCCTCTCATCTTAACC<br>Rev: CATCGCAGAGAACTGAGAGC</td><td>(AG)<sub>23</sub></td><td>Dirlewanger et al. (2002)</td></tr><tr><td>BPPCT025</td><td>6</td><td><i>P. persica</i> (L.) Batsch</td><td>For: TCCTGCGTAGAAGAAGGTAGC<br>Rev: CGACATAAAAGTCCAAATGGC</td><td>(GA)<sub>29</sub></td><td>Dirlewanger et al. (2002)</td></tr><tr><td>CPSCT026</td><td>7</td><td><i>P. salicina</i> Lindl.</td><td>For: TCTCACACGCTTTCGTCAAC<br>Rev: AAAAAGCCAAAAGGGGTTGT</td><td>(CT)<sub>16</sub></td><td>Mnejja et al. (2004)</td></tr><tr><td>CPPCT006</td><td>8</td><td><i>P. persica</i> (L.) Batsch</td><td>For: AATTAACTCCAACAGCTCCA<br>Rev: ATGGTTGCTTAATTCAATGG</td><td>(CT)<sub>16</sub></td><td>Aranzana et al. (2002)</td></tr></tbody></table>
allele size expected for a given SSR, we assumed that the corresponding primer pair amplified two loci. We considered this information when scoring the alleles of the polyploid species. Sometimes for the polyploid material the number of peaks did not match the ploidy level, in these cases, allele dosage was based on the height of the peak. Examples of electropherograms with recorded allele fragments are provided in supplementary material (Figures S2, S3).

Due to the polyploid nature of the *Prunus* species discussed in this work, calculation of classical genetic diversity parameters is not possible because of the ambiguous allele dosage. Hence, distribution of SSR alleles in the studied *Prunus* species was estimated as follows: number of detected alleles, alleles found in all species, and alleles common only to *P. domestica*, *P. cerasifera* or *P. spinosa*.

Bayesian model-based analysis was applied using the software package STRUCTURE V2.3.4 (Pritchard et al., 2000) to investigate the genetic structure of accessions and estimate the number of subpopulations. This software is suitable for data sets with different ploidy levels. Since the studied individuals were polyploids with ambiguous genotypes, the recessive allele approach was used (Falush et al., 2007). When all species were analysed together (diploids, tetraploids, and hexaploids), genotypes were coded according to Horvath et al. (2011). Among the options offered by the software, we chose the admixture model with correlated allele frequencies (Falush et al., 2007). The parameter K was set between 1 and 20 inferred clusters, with 20 independent iterations for each simulation. For each run, 100.000 burn-in periods followed by 750.000 MCMC (Markov Chain Monte Carlo) replicates, were performed.

The most relevant number of clusters (K) for the analysed data was estimated by calculating ΔK based on the method of Evanno (Evanno et al., 2005) implemented in the Structure Harvester V0.6.94 application (Earl and von Holdt, 2012). The threshold of 90% was chosen to assign a given individual to a population. The raw data (Q-matrix) obtained with the STRUCTURE software were used as a basis for the R package ‘tidyverse’ to visualise the results in the form of bar plots. In a second step, the Structure software was applied to a separate data set of *P. domestica* accessions with the same parameters.

Principal Coordinate Analysis (PCoA) was used to confirm the genetic structure of the three studied species and independently of the *P. domestica* accessions using DARwin 6.0.21 software (Perrier and Jacquemoud-Collet, 2006). This statistical method examines dissimilarities by converting data on distances between individuals into a location in low-dimensional space (e.g., 2D or 3D graphs). The software option single data format was used and allele data for each accession were converted into a binary matrix of presence (1) and absence (0) (Table S2). Despite the co-dominant nature of the SSR markers, the results were scored in a dominant manner. As Rouger and Jump (2014) noted, a part of genetic information is lost as a result, nevertheless this method has been successfully used in diversity studies of polyploid *Prunus* species (Horvath et al., 2011; Sehic et al., 2015; Halász et al., 2017; Urrestarazu et al., 2018).

First, dissimilarities (30.000 bootstraps) were calculated with Sokal and Michener index (simple matching):

$$d_{ij} = \frac{u}{m + u}$$

where $d_{ij}$ represents dissimilarity between units i and j, u is the number of unmatched variables and m is the number of matched variables. Information about the presence of the allele is rather confused in polyploids because allele frequencies are not clear. Therefore, the presence and absence modalities must be considered of equal weight, which minimises the loss of information (Noyer et al., 2003). However, some authors (Cordeiro et al., 2003) have used the Jaccard coefficient, which disregards the joint absence of bands in pairwise comparison and therefore most likely underestimates dissimilarity. On that ground, dissimilarities were also calculated using the Jaccard coefficient and the results were compared. The calculated dissimilarity was transformed into Euclidean distance using the 0.5 power transformation, and the graphical representation of the PCoA results was complemented using the by R package ‘scatterplot3d’. In addition, the calculated dissimilarity was used to estimate the genetic relationships among *P. domestica* accessions. Unweighted Neighbor-Joining (NJ) method was used to generate a tree with DARwin 6.0.21 software (Perrier and Jacquemoud-Collet, 2006). This method has been successfully used in similar studies (Campoy et al., 2016).

# 3 Results

## 3.1 Ploidy level determination of the plum species

According to relative DNA content estimated by flow cytometry, analysis of 124 studied accessions revealed 56 hexaploid, seven tetraploid and 41 diploid accessions (Table S3), with 2C nuclear DNA (pg) ranging from 1.35-1.698, 0.85-1.14 and 0.44-0.52, respectively (Table S4). The measuring was not successful for six samples, so their ploidy level was considered unknown.

Thirty-three accessions that were classified as *P. domestica* and considered hexaploid prior to analysis, based on recorded data from existing collections, along with morphological description, were found to be diploids (Table S3).

## 3.2 Diversity analyses of the plum species

### 3.2.1 Chloroplast DNA markers analysis

The combination of three universal cpDNA primers resulted in distinctive haplotypes representing the studied accessions. Among the analysed material, eight accessions did not show clear PCR amplifications. All seven-primer pair-restriction enzyme combinations (K1K2-Hinf I, K1K2-Taq I, K1K2-Alu I, VL-Hinf I, VL-Taq I, HK-Hinf I, and HK-Taq I) amplified polymorphic
patterns. After enzymatic digestion, ten polymorphic sites with 22
scorable bands were detected and ten haplotypes were identified
(H1-H10) (Table S5). The most prevalent haplotype was H1, which
was shared by 62 accessions (53%), followed by H4, which was
present in 27 accessions (23.1%). H2, H6 and H7 included 15
(12.8%), two (1.7%), and six accessions (5.1%), respectively. The
remaining haplotypes, H3, H5, H8, H9, and H10, were represented
by only one genotype each. Haplotype H1 represented hexaploid
material (P. domestica) with the exception of Blackthorn accession
(42), of unknown ploidy level. Most of the reference accessions
displayed affiliation to haplotype H1. Haplotypes H5 and H6
included solely accessions of hexaploid plum (P. domestica),
whereas H2 and H3 were specific to the diploid cluster (P.
kerasifera). Haplotype H4 was shared by these two Prunus
species, with P. cerasifera (22 accessions) being the most
abundant compared to the five P. domestica accessions.
Haplotypes H7 to H10 were found only in P. spinosa accessions.

The Neighbor-Joining tree (Figure 1) divided the haplotypes
into three clusters supported by a high bootstrap value (>50).
Cluster I was further divided into two sub-clusters, representing
predominantly *P. domestica* haplotypes (H1 and H6) and haplotype
H2 found only in *P. cerasifera*. Haplotype H5 consisted of only one
*P. domestica* accession: ‘Pitestean’ (69), which was not positioned
with the other hexaploid plums in the first sub-cluster. The cluster II
grouped together all *P. spinosa* accessions, which were further
subdivided into two sub-clusters, comprising haplotypes H7-H8
and H9-H10, respectively. In the cluster III, diploid plums were
mainly represented (H3 and H4, in two sub-clusters). H3 was
represented by only one *P. cerasifera* accession whereas sub-
cluster H4 was associated with diploid plums and five hexaploid
accessions: local Plum (52), ‘Empress’ (80), ‘Čačanska lepotica’ (87),
‘Quetsche du Carmel’ (169) and ‘Krikon’ (179).

**3.2.2 Polymorphism of SSR**

Based on 11 SSR markers, the analysis revealed 116 unique
genotypes. Identification of genotypes based on SSR profiles

revealed the presence of three duplicate groups accounting for
6.5% (8) of the studied accessions: Plum (150) and ‘Spindel plum’
(153); common prune accessions 108 and 110; and a group of four
Bluish plum accessions (112, 113, 114, and 145).

The 11 primer pairs used were polymorphic and amplified a single
locus with the exception of UDP96-005, EMPAS12 and BPPCT014,
which amplified two loci. The allele range for all primers extended
from 131-173 (UCD-CH17) to 133-219 bp (BPPCT025), whereas the
microsatellite markers (UDP96-005, EMPAS12 and BPPCT014),
which amplified two loci each, had a range of alleles of 94-174 bp,
94-156 bp and 187-283 bp, respectively (Table 2).

The allelic richness, generated by the 11 polymorphic primers,
varied from 22 (UCD-CH17, EMPAS14, CPPCT006) to 40
(BPPCT025) alleles per locus. The total number of alleles for all
accessions was 328 with an average of 29.82 alleles per locus
(Table 2). Within identified alleles, 81 were present in the three
studied species. The number and range of alleles was highest in P.
domestica, followed by P. spinosa and P. cerasifera. The distribution
of unique/specific alleles was the richest for P. domestica accessions
(28.6%), while P. cerasifera (12.7%) and P. spinosa (8.6%) had lower
values. We also calculated alleles found in only two species: P.
domestica and P. spinosa shared 83 alleles; 49 alleles were found in
P. domestica and P. cerasifera, while P. cerasifera and P. spinosa had
only 20 alleles in common (Table 2).

3.3 Genetic structure analyses of the
plum species

Bayesian analysis was used to examine the structural patterns of
the three species. Regarding the most likely number of clusters (K)
identified using the Evanno method, the maximum value for ΔK
was two (1213.8, Figure S4, Table S6), corresponding to two
ancestral populations (Figure 2). With a threshold of qI ≥ 0.90,
42 accessions representing diploid plums (P. cerasifera) were
assigned to Cluster 1 (green), with the exception of the hexaploid

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

*Accession Blackthorn (42) was initially classified as *P. spinosa*, but it is probably Bluish plum or an interspecific hybrid.

FIGURE 1

Neighbor-Joining tree based on cpDNA marker data for the three plum species (*P. cerasifera*, *P. domestica*, and *P. spinosa*). The complete list of studied accessions with the corresponding haplotypes and their codes can be found in Table S5.
<table><caption>TABLE 2 Distribution of SSR alleles in plum species (*P. domestica*, *P. cerasifera* and *P. spinosa*).</caption><thead><tr><th>SSR</th><th>BPPCT014_1*</th><th>BPPCT014_2*</th><th>BPPCT025</th><th>BPPCT034</th><th>CPPCT006</th><th>CPSCT026</th><th>EMPAS06</th><th>EMPAS11</th></tr></thead><tbody><tr><td>Linkage group</td><td>5</td><td>5</td><td>6</td><td>2</td><td>8</td><td>7</td><td>4</td><td>5</td></tr><tr><td>Total number of alleles</td><td>11</td><td>27</td><td>40</td><td>30</td><td>22</td><td>29</td><td>31</td><td>23</td></tr><tr><td>Number of alleles found in all three species</td><td>4</td><td>0</td><td>9</td><td>6</td><td>6</td><td>8</td><td>7</td><td>9</td></tr><tr><td>Range of alleles P. domestica</td><td>187-209</td><td>213-283</td><td>149-219</td><td>215-273</td><td>156-224</td><td>161-221</td><td>205-261</td><td>57-103</td></tr><tr><td>Total number of alleles in P. domestica</td><td>9</td><td>26</td><td>32</td><td>26</td><td>22</td><td>23</td><td>26</td><td>23</td></tr><tr><td>Number of alleles specific to P. domestica</td><td>0</td><td>16</td><td>15</td><td>2</td><td>6</td><td>7</td><td>11</td><td>5</td></tr><tr><td>Range of alleles P. cerasifera</td><td>189-199</td><td>0</td><td>133-197</td><td>217-257</td><td>172-202</td><td>171-211</td><td>179-245</td><td>57-79</td></tr><tr><td>Total number of alleles in P. cerasifera</td><td>5</td><td>0</td><td>18</td><td>11</td><td>7</td><td>17</td><td>13</td><td>9</td></tr><tr><td>Number of alleles specific to P. cerasifera</td><td>1</td><td>0</td><td>6</td><td>1</td><td>0</td><td>5</td><td>1</td><td>0</td></tr><tr><td>Range of alleles P. spinosa</td><td>187-209</td><td>211-251</td><td>143-193</td><td>211-263</td><td>174-210</td><td>171-215</td><td>197-254</td><td>57-99</td></tr><tr><td>Total number of alleles in P. spinosa</td><td>9</td><td>11</td><td>16</td><td>15</td><td>15</td><td>13</td><td>15</td><td>17</td></tr><tr><td>Number of alleles specific to P. spinosa</td><td>1</td><td>1</td><td>2</td><td>3</td><td>0</td><td>0</td><td>3</td><td>0</td></tr><tr><td>Number of alleles common to P. domestica and P. cerasifera</td><td>1</td><td>0</td><td>3</td><td>3</td><td>1</td><td>2</td><td>4</td><td>9</td></tr><tr><td>Number of alleles common to P. domestica and P. spinosa</td><td>5</td><td>10</td><td>1</td><td>5</td><td>9</td><td>3</td><td>4</td><td>17</td></tr><tr><td>Number of alleles common to P. cerasifera and P. spinosa</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td><td>9</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 2 Continued</caption><thead><tr><th>SSR</th><th>EMPAS12_1*</th><th>EMPAS12_2*</th><th>EMPAS14</th><th>UCD-CH17</th><th>UDP96-005_1*</th><th>UDP96-005_2*</th><th>TOTALS</th><th>MEAN</th></tr></thead><tbody><tr><td>Linkage group</td><td>3</td><td>3</td><td>5</td><td>4</td><td>1</td><td>1</td><td>/</td><td>/</td></tr><tr><td>Total number of alleles</td><td>14</td><td>18</td><td>22</td><td>22</td><td>18</td><td>21</td><td>328</td><td>29.82</td></tr><tr><td>Number of alleles found in all three species</td><td>8</td><td>5</td><td>4</td><td>5</td><td>5</td><td>10</td><td>81</td><td>24.7</td></tr><tr><td>Range of alleles P. domestica</td><td>94-120</td><td>122-156</td><td>165-199</td><td>131-173</td><td>94-130</td><td>132-174</td><td>/</td><td>/</td></tr><tr><td>Total number of alleles in P. domestica</td><td>13</td><td>18</td><td>13</td><td>21</td><td>18</td><td>20</td><td>290</td><td>26.36</td></tr><tr><td>Number of alleles specific to P. domestica</td><td>3</td><td>5</td><td>1</td><td>4</td><td>3</td><td>4</td><td>83</td><td>28.62</td></tr><tr><td>Range of alleles P. cerasifera</td><td>94-118</td><td>122-148</td><td>171-219</td><td>135-167</td><td>96-118</td><td>132-174</td><td>/</td><td>/</td></tr><tr><td>Total number of alleles in P. cerasifera</td><td>9</td><td>10</td><td>9</td><td>12</td><td>9</td><td>13</td><td>142</td><td>12.9</td></tr><tr><td>Number of alleles specific to P. cerasifera</td><td>0</td><td>0</td><td>3</td><td>1</td><td>0</td><td>0</td><td>18</td><td>12.68</td></tr><tr><td>Range of alleles P. spinosa</td><td>94-118</td><td>122-154</td><td>167-207</td><td>139-171</td><td>100-130</td><td>140-172</td><td>/</td><td>/</td></tr><tr><td>Total number of alleles in P. spinosa</td><td>10</td><td>7</td><td>16</td><td>10</td><td>11</td><td>10</td><td>175</td><td>15.9</td></tr><tr><td>Number of alleles specific to P. spinosa</td><td>0</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td><td>15</td><td>8.57</td></tr><tr><td>Number of alleles common to P. domestica and P. cerasifera</td><td>0</td><td>5</td><td>4</td><td>6</td><td>4</td><td>7</td><td>49</td><td>15.65</td></tr><tr><td>Number of alleles common to P. domestica and P. spinosa</td><td>1</td><td>2</td><td>6</td><td>10</td><td>6</td><td>4</td><td>83</td><td>26.77</td></tr><tr><td>Number of alleles common to P. cerasifera and P. spinosa</td><td>1</td><td>0</td><td>2</td><td>5</td><td>0</td><td>1</td><td>20</td><td>8.16</td></tr></tbody></table>

<footnote>* Primers pairs of UDP96-005, EMPAS12 and BPPCT014 amplified two loci.</footnote>
<figure><img src="image_4.png" /><figcaption>FIGURE 2</figcaption></figure>

Bar plot of Bayesian analysis results (K = 2) for plum species genotypes. *P. cerasifera* accessions are shown in green (Cluster 1), and *P. domestica* accessions are arranged in blue (Cluster 2). *P. spinosa* accessions were admixed.

accession Plum (52). Sixty-nine *P. domestica* accessions were assigned to Cluster 2 (blue), which were mainly identified as hexaploid, with the exception of the accessions: ‘Stanley’ (73), and Plum (144), which were found to be diploid in flow cytometry analysis. In addition, twelve genotypes showed mixed ancestry (membership values lower than 90% in both clusters). Among the admixed accessions, the highest number of accessions (10) belonged to *P. spinosa* with the remaining two genotypes representing a diploid accession, Plum (162) and one hexaploid, ‘Krikon’ (179).

Before performing the Principal Coordinate Analysis, dissimilarities were estimated using two distinctive indices: Sokal and Michener index and Jaccard coefficient. Comparison of the results showed no significant differences between the two indices, and only the results of the first index are presented. PCoA analysis revealed three main clusters (Figure 3). Cluster 1, corresponding to *P. cerasifera*, formed a denser arrangement in space. In contrast,

Cluster 2, corresponding to *P. domestica* accessions, was more dispersed. The admixed accessions (mainly *P. spinosa*) were distributed between the two clusters but closer to Cluster 1 (*P. cerasifera*).

Based on the results of various analyses (flow cytometry, cpDNA, Structure, PCoA and Neighbor-Joining), thirty-one accessions were identified, among the material collected in Slovenia, that were assumed to be hexaploid *P. domestica* but actually belonged to the diploid species *P. cerasifera*.

## 3.4 Genetic structure analyses and relationships of *P. domestica*

One of the main objectives of this study was to assess the genetic structure of *P. domestica* collected in Slovenia.

<figure><img src="image_5.png" /><figcaption>FIGURE 3</figcaption></figure>

3D graph of Principal Coordinate Analysis (PCoA) based on polymorphism at 11 SSR loci of the 124 genotypes. The colors of the accessions indicate the most relevant parameter K = 2: green for *P. cerasifera* accessions, blue for mainly *P. domestica* material, and black for admixed accessions.
<figure><img src="image_6.png" /><figcaption>FIGURE 4
Bar plot of Bayesian analysis results (K = 2) for *P. domestica* accessions. Twenty-two accessions representing mainly local traditional material collected in Slovenia, were assigned to the Cluster 1 (blue). Cluster 2 (violet) consisted of 32 accessions of mixed origin (*in situ* and *ex situ*) representing different pomological groups. Seventeen accessions were classified admixed.</figcaption></figure>

When investigating patterns of structure among *P. domestica* accessions with Bayesian analysis, two ancestral populations were identified using the ΔK method (Figure S5, Table S7). The threshold value of qI ≥ 0.90 assigned twenty-two accessions to the Cluster 1 (blue). They represented mainly the local traditional *P. domestica* accessions of Bluish plum and common prunes collected in Slovenia, and additionally two traditional German prunes: ‘HZW Meschenmoser’ (74) and ‘HZW typ Mare’ (89) (Figure 4). Cluster 2 (violet) consisted of 32 accessions of mixed origin, including: five accessions collected *in situ* in Slovenia, Austria and Croatia; four improved cultivars [‘Stanley’ (73), ‘Hanita’ (75), ‘Ruth Gerstetter’ (81) and ‘Topper’ (86)]; eleven accessions from the SPGB collection, the Botanical Garden of the University of Maribor, the Botanical Garden of the Karl-Franzens University of Graz, and twelve

accessions from the PFNC collection. The latter are reference accessions representing different pomological groups (such as mirabelles and greengages), some French landraces [such as ‘Quetsche du Carmel’ (169), ‘D’Ente double’ (170), ‘Damas blanc’ (171), ‘Perdrigon’ (177), ‘St. Julien’ (180)], and genotypes such as damson ‘Krikon’ (179). Seventeen accessions with estimated membership coefficients below 0.9 were considered admixed, including five accessions from Bluish plum group (44, 45, 46, 50 and 189). The average distance (expected heterozygosity) between individuals was higher in Cluster 2 (0.8627) than in Cluster 1 (0.6617).

Examining the structure of the plum collection for K ranging between 3 and 6 (Figure S6A–D) allowed us to detect interesting dynamics among the different plum groups. At K=3 (Figure S6A),

<figure><img src="image_7.png" /><figcaption>FIGURE 5
3D graph of Principal Coordinate Analysis (PCoA) based on polymorphism at 11 SSR loci of the 71 *P. domestica* genotypes. The colors of the accessions indicated the clusters obtained by Bayesian analysis for K=2: Bluish plum accessions and common prunes collected in Slovenia (blue), plums genetically very diverse (violet). Accessions located between the two clusters (in black) correspond to the admixed accessions.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>FIGURE 6</figcaption></figure>

Neighbor-Joining tree based on the dissimilarity matrix calculated from the data set of 11 SSRs for 71 P. domestica accessions using the Sokal and Michener index. The colors of the accessions reflect the type of studied material: *in situ* material (blue) and *ex situ* material (black).

Cluster 1 split into two sub-clusters, Bluish plum group (olive green) and the common prunes (blue), while the original Cluster 2 (violet) and admixed accessions maintained the same profile or remained unchanged, except for the accession 'Kriecherl Innerhofer' (194), which shared a large part of its genome with the Bluish plum group. When K=4 (Figure S6B) was considered, 14 accessions (red) separated from Cluster 2 (K=2). They corresponded to known pomological groups such as mirabelle plums (140, 175, 176) and greengage accessions (172, 173, 174), as well as 'Damas blanc' (171), 'St. Julien' (180) and 'Italian prune' (57). At K=5 (Figure S6C), Bluish plum accessions 44, 45, 46 (before admixed) formed a new cluster (yellow) that also included accessions Plum (150) and 'Spindel plum' (153). At K=6 (Figure S6D), Plum (150) and 'Spindel plum' (153) were associated with five other accessions: 'Plum red' (65), 'Stanley' (73), 'Čačanska lepotica' (87), 'President' (88) and 'Jojo' (127), forming an additional ancestral population (cyan).

PCoA confirmed the results obtained with Structure for K=2 (Figure 5). Cluster 1 was represented by Bluish plum accessions and common prunes collected in Slovenia. Cluster 2 consisted of 32 genetically very diverse materials: some plum accessions collected in situ, international improved cultivars, and reference accessions (described in detail above). It was more dispersed and divided into several sub-clusters, demonstrating a higher genetic diversity. The accessions located between the two clusters corresponded to the admixed accessions identified by the Bayesian analysis.

Neighbor-Joining analysis allowed us to evaluate the genetic relationships between 71 P. domestica accessions (Figure 6, Figure S7). Three clusters were observed in the dendrogram. Cluster I was the greatest (44 accessions), quite complex and divided into two

sub-clusters. The first sub-cluster included plum accessions that were very different morphologically. This sub-cluster could be related to Cluster 2 identified by Bayesian analysis. First, all French greengages were established: 'Reine Claude' (172), 'Reine Claude Violette' (173), and 'Reine Claude Dorée' (174), together with an improved cultivar, 'Plum green' (56), from Slovenia. The French landraces 'Perdrigonne Violette' (167), 'Damas Blanc' (171) and 'Saint Julien' (180) were associated to this group, but with low bootstrap values. In the second group, two Slovenian accessions, Mirabelle (140) and supposed Myrobolan (141), clustered very strongly with French mirabelle reference cultivars, e.g., 'Mirabelle de Nancy' (175) and 'Mirabelle de Metz' (176). Two accessions were placed near the mirabelle group, Plum (144) and an 'Italian Prune' (57). Below, two pairings stood out: 'Zibate blue' (191) and 'Zibate yellow' (192), native small-fruited subsp. *insititia* genotypes from Austria; 'Krikon' (179) clustered together with a local accession Plum (52). At the end of the first part of this sub-cluster, a pair of duplicate genotypes, Plum (150) with 'Spindel plum' (153), was found. The second half of the first sub-cluster consisted mainly of well-known improved cultivars: 'Ruth Gerstetter' (81), which appeared to be genetically close to the landrace 'Californian Prune' (66), 'Pitestean' (69), and 'Topper' (86). Associated with the latter were an improved variety 'Plum Blue' (53) maintained in the SPGB and a very old variety from Croatia, Plum (85), which is a prune type. The Slovenian small fruited genotype 'Plum Green' (37) was connected to this sub-cluster, although with a low bootstrap value. The second sub-cluster of Cluster I consisted of common prunes, collected *in situ* in Slovenia, well-known cultivars such as 'Bistrica' (79, 84) or 'Quetsche commune' (168), or derived from clonal selection ['HZW Meschenmoser' (74) and 'HZW Typ Mare']
(89)] or from crosses involving common prunes [‘Ersinger’ (72), ‘Hanita’ (75) and ‘Presenta’ (78)].

The Cluster II included all accessions of the Bluish plum group, separated into two groups according to the results of Bayesian analysis. One corresponded to Structure Cluster 1 and the other combined the accessions previously identified as admixed in the genetic structure analysis. The first group included material from different locations in the Styria region, associated to them, but clearly separated, was accession ‘Kriecherl Innerhofer’ (194) from Austria. The second group, much smaller, consisted of three accessions collected in the Prekmurje region. Related to both groups were found a *Prunus insititia* accession ‘Tersen’ (178), originating from Sweden and maintained in the PFNC collection, and the last Bluish plum (189), collected in eastern Slovenia. Cluster III contained several very well-known modern cultivars such as ‘Elena’ (68), ‘Čačanska lepotica’ (87), ‘President’ (88), ‘Jojo’ (127) and ‘Stanley’ (73).

# 4 Discussion

## 4.1 Relationship among *P. spinosa*, *P. cerasifera* and *P. domestica* species

In this study we tried to collect mainly *P. domestica* accessions/ germplasm, but we discovered that approximately a third of the collected material actually belonged to *P. cerasifera*. This surprising result might be related to the fact that the Slovenian term ‘ringlo’, traditionally used for greengages, can also refer to any plum that is round, which is a typical characteristic of *P. cerasifera* fruits. When we carried out the initial prospection, we did not aim at species identification, but at categorizing accessions based on fruit description or other interesting agronomic traits.

Assessment of genetic diversity using PCR-RFLP analyses showed that the total number of haplotypes (10) was lower compared to similar studies: *P. spinosa* (4), *P. domestica* (3), *P. cerasifera* (2), and, as mentioned above, one haplotype was shared by the latter two species. Horvath et al. (2011) discovered the highest total number of haplotypes (32): *P. cerasifera* (15), *P. spinosa* (12) and *P. domestica* (5). On the other hand, Urrestarazu et al. (2018) reported the highest number of *P. domestica* haplotypes (8), which is consistent with the fact that they studied the highest number of *P. domestica* accessions (166). The most predominant haplotype in our study contained mainly *P. domestica* accessions. Urrestarazu et al. (2018) also reported the presence of a predominant haplotype, whereas two major haplotypes were identified in the case of Horvath et al. (2011). Both authors argue that a reduction in cpDNA haplotypes of *P. domestica* could be related to historical events (geological events, climatic oscillations, and independent dispersal of the species across continents), human interference and a narrow genetic base of material when plum was introduced into Europe. However, the relatively small number of accessions in our study could also explain this observation. Concerning *P. spinosa*, Horvath et al. (2011) observed more

complex clustering because although most haplotypes of *P. spinosa* clustered together (8), four were linked to haplotypes of *P. domestica* and *P. cerasifera*. Because several identical accessions were examined in both studies (21) and the same primer pair-restriction enzyme combinations were used for cpDNA marker analysis, this discrepancy is rather unexpected. For 17 accessions, we confirmed the same trend for the corresponding haplotypes in both studies. However, three *P. spinosa* (181, 183, 184) accessions did not form individual haplotypes as was the case in Horvath et al. (2011). We suspect that this discrepancy is due to the difference in resolution of the amplified products.

When studying SSR markers polymorphism, we obtained an average value of 29.82 alleles per locus, similar to the value of 29.25 alleles for eight SSR loci reported by Gaši et al. (2020). In contrast, Horvath et al. (2011) found an average of 41 alleles for five loci, while Urrestarazu et al. (2018) determined an average of 23.36 alleles after analysing 11 SSRs. Slightly lower values of 22.7 and 18.7 alleles for nine SSR loci were determined by Sehic et al. (2015) and Kazija et al. (2014), respectively. These studies differ in the evaluated material. While Sehic et al. (2015), Urrestarazu et al. (2018), and Gaši et al. (2020) studied exclusively *P. domestica* accessions, Horvath et al. (2011) and our work included accessions of *P. domestica*, *P. cerasifera* and *P. spinosa* species. Our results showed lower diversity than those of Horvath et al. (2011), which could be due to the fact that our study focused more on local landraces such as the Bluish plum group and common prunes from Slovenia, while Horvath et al. (2011) included material from different countries across Europe, especially for the *P. spinosa* and *P. cerasifera* species.

When comparing all three species, the number, range of alleles, and distribution of unique/specific alleles were highest for *P. domestica*, which is not surprising given its polyploid nature as well as the largest number of accessions studied. The number of *P. spinosa* accessions for this study was the lowest of all three species, nevertheless the number and range of alleles were higher compared to *P. cerasifera*. This demonstrates the originality of the collected *P. spinosa* material.

In addition, cpDNA analysis complemented with SSR markers allowed us to reveal some examples of complex relationships between the three species. Firstly, the supposedly *P. spinosa* accession Blackthorn (42), shared haplotype H1 with *P. domestica* material. Furthermore, the results of a Neighbor-Joining analysis for the three species (data not shown) indicated that Blackthorn (42) clustered together with a group of Bluish plum accessions (44, 45, and 46). However, Bayesian analysis classified this accession as admixed, together with other *P. spinosa* accessions. Hence, genetic relationship analyses suggest that Blackthorn (42) accession could be either a Bluish plum or an interspecific hybrid. Secondly, ‘Krikon’ (179) is another accession that exhibited a putative hybrid origin. When all species were examined together (data not shown), ‘Krikon’ clustered with the *P. cerasifera* and *P. spinosa* material. Interestingly, Zhebentyayeva et al. (2019) found that the same ‘Krikon’ accession from the PFNC collection clustered close to a *P. cerasifera* accession from the UC Davis collection and a *P. spinosa* accession from Sweden.
## 4.2 Involvement of *P. cerasifera* and *P. spinosa* in *P. domestica* origin

The combined results presented in this study allowed us to investigate the potential involvement of *P. cerasifera* and *P. spinosa* species in the origin of *P. domestica*. Different analyses provided somehow contradictory information, which is not surprising given the unresolved question of the genetic origin of *P. domestica*. Assessment of genetic diversity using chloroplast DNA markers showed that *P. spinosa* clearly clustered separately, whereas *P. domestica* and *P. cerasifera* clustered together. The haplotypes of *P. domestica* and *P. cerasifera* appeared to be genetically closer compared to the haplotypes of *P. spinosa*, suggesting that *P. cerasifera* may have contributed to the maternal chloroplast DNA of *P. domestica*. This was also discussed by Bortiri et al. (2009), Reales et al. (2009), and Horvath et al. (2011). On the contrary, SSR markers revealed a possible involvement of *P. spinosa* in the genome of *P. domestica*, sharing almost twice as many alleles as *P. domestica* with *P. cerasifera*, which is consistent with results of Horvath et al. (2011). This could support a hybrid origin of the European plum. Moreover, the genetic structure analyses revealed *P. domestica* and *P. cerasifera* as independent groups, while the accessions of *P. spinosa* showed admixed ancestry shared with both species.

## 4.3 Genetic structure and relationships between *P. domestica* accessions

In this study, different approaches allowed us to establish genetic relationships among the local plum material collected in Slovenia and the traditional plum groups, such as common prunes, greengages and mirabelles.

Duplicates, represented 11.3% of the studied *P. domestica* accessions. Kazija et al. (2014) found a similar proportion of duplicates (11%), when they examined the plum germplasm in Croatia. The percentage of identified duplicates was slightly higher in Urrestarazu et al. (2018) at 18%, while Sehic et al. (2015) detected 30% duplicates.

The STRUCTURE results were confirmed by the representation of PCoA analysis based on a genetic distance matrix. The studied *P. domestica* accessions formed two distinct ancestral populations, which was also reported for plum by Gaši et al. (2020), as well as for other germplasm collections of important fruit tree species such as apple (Urrestarazu et al., 2012), cherry (Campoy et al., 2016), walnut (Bernard et al., 2018) and apricot (Bourguiba et al., 2020). The average distance between individuals (expected heterozygosity) was higher in Cluster 2, which could be explained by the diverse origin of the accessions belonging to this cluster (landraces from the different pomological groups, accessions from botanical gardens, and international modern cultivars) and by their higher number. On the contrary, Cluster 1 consisted mainly of two plum groups (Bluish plum group and common prunes) collected in Slovenia.

The landraces belonging to the Bluish plum group are historically and genetically important plums that are slowly disappearing from the Slovenian environment. Since there is very scarce information about this plum group (except in oral tradition),

we have tried to find reference material, corresponding to the morphological description of this varietal group in order to link them with some known genetic material. The French landrace 'Perdrigon' is a small, round-fruited, bluish-purple plum that was cultivated in Slovenia in the 19th century (Zalokar, 1854). However, two accessions of this landrace group, 'Perdrignonne Violette' (167) and 'Perdrigon' (177), as well as the accession 'Saint Julien' (180), which are morphologically similar to the description of Bluish plum landraces, were not genetically close to the Bluish plum group. The same is true for 'Krikon' (179), which is probably of Swedish origin and represents a group of sour-tasting feral plums with small, round, and bluish (bluish-red) fruits (Sehic et al., 2015). On the other hand, *Prunus insititia* 'Tersen' (178), which also originates from Sweden, clustered with several Bluish plum accessions. The fruits of 'Tersen' are small, blue-black in colour, and traditionally used for jams. Finally, 'Kriecherl Innerhofer' (194) is a clonal selection within subsp. *insititia*, a small, round, bluish-purple coloured landrace from Styria region in Austria. It is traditionally used for the production of plum brandy. Because its geographic origin, morphological appearance, and use overlap with our Bluish plum landraces, we suspected that they might be genetically close, which was confirmed by our analyses. The SSR analyses divided the Bluish plum accessions into two distinct populations, which were established in agreement with the region of origin. Similar to the first population, the Bluish plum accession (189) was also collected in the same region (Styria), although approximately 75 km of air distance apart and 100 m difference in elevation. It was in the same cluster with other Bluish plum populations, but separated from them and connected with *Prunus insititia* 'Tersen' (178). Since it seems to be a rather original material, it should be further investigated, by studying a larger sample at the same location. While Mišić (1979) described Bluish plum to be one of the widespread landraces in Slovenia, this did not correspond to our experience while collecting material. Bluish plum has low to medium susceptibility to some important plum diseases like *Plum pox virus*, *Monilia laxa* and *Monilia fructigena*. Perhaps this is the reason why they are slowly disappearing from the natural environment. Nevertheless, it might be interesting to test genetically well-identified Bluish plum accessions for tolerance to these diseases. There is also a possibility that Slovenian fruit producers and local farmers have tried to introduce more modern cultivars in recent years and somehow 'abandoned' the traditional material.

Another important locally collected group of plums are the common prunes, which clearly formed a unique group and seemed to include similar, but not identical accessions. When structure was examined (K=3), common prunes originating from different parts of Europe, such as Slovenia, France or Germany, formed a separate cluster, which was also confirmed by the Neighbor-Joining analysis. This result corroborates the fact that common prunes are widespread in Europe, with a large number of types that are distributed under different names such as 'Hauszwetsche' (Germany), 'Quetsche Commune' (France), 'Casalinga' (Italy), 'Besztercei' (Hungary) and 'Vinete romanesti' (Romania) (Mišić, 1979). Linked by parentage to the common prune cluster were early selections or modern cultivars: 'Ersinger' (72), 'Hanita' (75)
comprehensive study is the common prunes. Moreover, this type of research could be a next step for a comprehensive evaluation of plant genetic resources in the Balkan Peninsula, which harbors potentially valuable and unique plum genotypes.

Finally, the obtained molecular data and the characterization of the collected material will be crucial in order to understand the diversity of Slovenian germplasm, to improve the conservation process, to recover local genotypes and to enrich the current SPGB collection. Furthermore, a complete agro-morphological study would be essential to evaluate prospective material in an attempt to integrate valuable genotypes into breeding programs. In addition, a thorough evaluation of commercial value for use in plum production would also be required.

## Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

## Author contributions

AI, MŠ, AŠ and TT designed the initial research, TB and JQ-G suggested complementary research and methods in data analysis. AI, MŠ, AŠ and TT conducted the field work and sampling. TT carried out SSR molecular analyses. TT and TB evaluated the SSR analyses results. TT analyzed the data. TB and JQ-G contributed to the data analysis. TT, TB and JQ-G, wrote the manuscript, AI, MŠ, AŠ edited and reviewed the original manuscript. All authors contributed to the article and approved the submitted version.

## Funding

TT was financially supported by the following grants: French Government Scholarship Program (BGF) and WFS National scholarship for Slovenian PhD students, offered by the Slovenian Science Foundation.

## References

* Aranzana, M. J., Garcia-Mas, J., Carbó, J., and Arús, P. (2002). Development and variability analysis of microsatellite markers in peach. *Plant Breed.* 121, 87–92. doi: 10.1046/j.1439-0523.2002.00656.x

* Arumuganathan, K., and Earle, E. D. (1991). Nuclear DNA content of some important plant species. *Plant Mol. Biol. Rep.* 9, 208–218. doi: 10.1007/bf02672069

* Bajashvili, E. I. (1990). Studies of some species of *Prunus* mill. *Genus. Acta Hortic.* 283, 31–34. doi: 10.17660/ActaHortic.1990.283.3

* Benediková, D., and Giovannini, D. (2011). Review on genetic resources in the ECPGR prunus working group. *Acta Hortic.* 981, 43–51. doi: 10.17660/ActaHortic.2013.981.2

* Bernard, A., Barreneche, T., Lheureux, F., and Dirlewanger, E. (2018). Analysis of genetic diversity and structure in a worldwide walnut (*Juglans regia* L.) germplasm using SSR markers. *PloS One* 13, 19. doi: 10.1371/journal.pone.0208021

* Blazek, J. (2007). A survey of the genetic resources used in plum breeding. *Acta Hortic.* 734, 31–45. doi: 10.17660/ActaHortic.2007.734.2

## Acknowledgments

The authors would like to thank the INRAE's 'Prunus Juglans Biological Resources Centre' of the Fruit Tree Experimental Unit of INRAE-Nouvelle Aquitaine-Bordeaux (UE0393), Slovene Plant Gene Bank, the Botanical Garden of the University of Vienna, the Botanical Garden of the Karl-Franzens University of Graz, the Botanical Garden of the University of Ljubljana, the Botanical garden of the University of Maribor, Sadjarstvo Mirosan d.o.o and Versuchsstation Obst- und Weinbau Haidegg for providing valuable genetic material used in this study. We would like to express our gratitude to Borut Bohanec and Ester Stajič from the Biotechnical Faculty (University of Ljubljana) for their assistance in determining the ploidy level.

## Conflict of interest

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

## Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

## Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1150459/full#supplementary-material

* Bohanec, B.M. (2003). "Ploidy determination using flow cytometry," in *Doubled haploid production in crop plants: A manual*. Eds. Maluszynski, K. J. Kasha, B. P. Forster and I. Szarejko (Netherland: Speinger Science Buisiness media LLC), 397–403.

* Bortiri, E., Oh, S.-H., Jiang, J., Baggett, S., Granger, A., Weeks, C., et al. (2009). 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. doi: 10.1043/0363-6445-26.4.797

* Bouhadida, M., Martin, J. P., Eremin, G., Pinochet, J., Moreno, M.Á., and Gogorcena, Y. (2007). Chloroplast DNA diversity in *Prunus* and its implication on genetic relationships. *J. Am. Soc Hortic. Sci.* 132, 670–679. doi: 10.21273/jashs.132.5.670

* Bourguiba, H., Scotti, I., Sauvage, C., Zhebentyayeva, T., Ledbetter, C., Krška, B., et al. (2020). Genetic structure of a worldwide germplasm collection of *Prunus armeniaca* L. reveals three major diffusion routes for varieties coming from the species' center of origin. *Front. Plant Sci.* 11. doi: 10.3389/fpls.2020.00638
Campoy, J. A., Lerigoleur-Balsemin, E., Christmann, H., Beauvieux, R., Girollet, N., Quero-García, J., et al. (2016). Genetic diversity, linkage disequilibrium, population structure and construction of a core collection of *Prunus avium* l. landraces and bred cultivars. *BMC Plant Biol.* 16, 49. doi: 10.1186/s12870-016-0712-9

Chin, S. W., Shaw, J., Haberle, R., Wen, J., and Potter, D. (2014). Diversification of almonds, peaches, plums and cherries - molecular systematics and biogeographic history of *Prunus* (Rosaceae). *Mol. Phylogenet. Evol.* 76, 34–48. doi: 10.1016/j.ympev.2014.02.024

Cipriani, G., Lot, G., Huang, W.-G., Marrazzo, M. T., Peterlunger, E., and Testolin, R. (1999). AC/GT and AG/CT microsatellite repeats in peach [*Prunus persica* (L.) batsch]: isolation, characterisation and cross-species amplification in *Prunus*. *theor. Appl. Genet.* 99, 65–72. doi: 10.1007/s001220051209

Cordeiro, G., Pan, Y.-B., and Henry, R. (2003). Sugarcane microsatellites for assessment of genetic diversity in sugarcane germplasm. *Plant Sci.* 165, 181–189. doi: 10.1016/S0168-9452(03)00157-2

Crane, M. B., and Lawrence, W. J. C. (1938). *The genetics of garden plants* (London: Macmillan and CO., Limited).

Decroocq, V., Hagen, L. S., Favé, M. G., Eyquard, J. P., and Pierronnet, A. (2004). Microsatellite markers in the hexaploid *Prunus domestica* species and parentage lineage of three European plum cultivars using nuclear and chloroplast simple-sequence repeats. *Molec. Breed.* 13, 135–142. doi: 10.1023/B:MOLB.0000018761.04559.b3

Demesure, B., Sodzi, N., and Petit, R. (1995). A set of universal primers for amplification of polymorphic non-coding regions of mitochondrial and chloroplast DNA in plants. *Mol. Ecol.* 4, 129–132. doi: 10.1111/j.1365-294X.1995.tb00201.x

Depypere, L., Chaerle, P., Vander Mijnsbrugge, K., and Goetghebeur, P. (2009). Classification trees and plant identification: A case study of European *Prunus* section *Prunus* taxa. *Belg. J. Bot.* 142, 163–176.

Dirlewanger, E., Cosson, P., Tavaud, M., Aranzana, J., Poizat, C., Zanetto, A., et al. (2002). Development of microsatellite markers in peach [*Prunus persica* (L.) batsch] and their use in genetic diversity analysis in peach and sweet cherry [*Prunus avium* l.]. *Theor. Appl. Genet.* 105, 127–138. doi: 10.1007/s00122-002-0867-7

Dolezel, J. (1991). Flow cytometric analysis of nuclear DNA content in higher plants. *Phytochem. Anal.* 2, 143–154. doi: 10.1002/pca.2800020402

Dolezel, J., and Bartos, J. (2005). Plant DNA flow cytometry and estimation of nuclear genome size. *Ann. Bot.* 95, 99–110. doi: 10.1093/aob/mci005

Doyle, J. J., and Doyle, J. L. (1987). A rapid DNA isolation procedure for small amounts of fresh leaf tissue. *Phytochem. Bull.* 19, 11–15.

Dumolin-Lapegue, S., Pemonge, M. H., and Petit, R. J. (1997). An enlarged set of consensus primers for the study of organelle DNA in plants. *Mol. Ecol.* 6, 393–397. doi: 10.1046/j.1365-294X.1997.00193.x

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

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

Falush, D., Stephens, M., and Pritchard, J. (2007). Inference of population STRUCTURE using multilocus genotype data: dominant markers and null alleles. *Mol. Ecol. Notes.* 7, 574–578. doi: 10.1111/j.1471-8286.2007.01758.x

Faust, M., and Suranyi, D. (1999). Origin and dissemination of plums. *Hortic. Rev.* 23, 179–231. doi: 10.1002/9780470650752.ch4

Gaši, F., Sehic, J., Grahic, J., Hjeltnes, S. H., Ordidge, M., Benedikova, D., et al. (2020). Genetic assessment of the pomological classification of plum *Prunus domestica* l. *accessions sampled across Europe*. *Genet. Resour. Crop Evol.* 67, 1137–1161. doi: 10.1007/s10722-020-00901-y

Halász, J., Makovics-Zsóhár, N., Szőke, F., Ercisli, S., and Hegedűs, A. (2017). Simple sequence repeat and s-locus genotyping to explore genetic variability in polyploid *Prunus spinosa* and p. *insititia*. *Biochem. Genet.* 55, 22–33. doi: 10.1007/s10528-016-9768-3

Hartmann, W., and Eckhart, F. (2008). *Farbatlas alte obstsorten* (Stuttgart: Eugen Ulmer).

Hedrick, U. P. (1911). *The plums of new York* (Albany: J. B. Lyon company, State printers).

Heinkel, R., Hartmann, W., and Stösser, R. (2000). On the origin of the plum cultivars 'Cacaks beauty', 'Cacaks best', 'Cacaks early' and 'Cacaks fruitful' as investigated by the inheritance of random amplified polymorphic DNA (RAPD) fragments. *Sci. Hortic.* 83, 149–155. doi: 10.1016/S0304-4238(99)00074-6

Honzak, D. (1968). *Sadjarstvo* (Ljubljana: Tehniška založba Slovenije).

Horvath, A., Balsemin, E., Barbot, J.-C., Christmann, H., Manzano, G., Reynet, P., et al. (2011). 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. doi: 10.1016/j.sciencia.2011.03.049

Horvath, A., Christmann, H., and Laigret, F. (2008). Genetic diversity and relationships among *Prunus cerasifera* (cherry plum) clones. *Botany.* 86, 1311–1318. doi: 10.1139/B08-097

Ilgin, M., Kafkas, S., and Ercisli, S. (2014). Molecular characterization of plum cultivars by AFLP markers. *Biotechnol. Biotechnol. Equip.* 23, 1189–1193. doi: 10.1080/13102818.2009.10817636

Insero, O., Monastra, F., Guerriero, R., Fontanazza, G., Godini, A., Tamponi, G., et al. (1991). *Frutticoltura speciale* (Roma: Reda).

Kaufmane, E., Ikase, L., Trajkovski, V., and Lacis, G. (2002). Evaluation and characterization of plum genetic resources in Sweden and Latvia. *Acta Hortic.* 577, 207–213. doi: 10.17660/ActaHortic.2002.577.34

Kazija, D. H., Jelačić, T., Vujević, P., Milinović, B., Čiček, D., Biško, A., et al. (2014). Plum germplasm in Croatia and neighboring countries assessed by microsatellites and DUS descriptors. *Tree Genet. Genomes.* 10, 761–778. doi: 10.1007/s11295-014-0721-5

Liu, W., Liu, D., Zhang, A., Feng, C., Yang, J., Yoon, J., et al. (2007). Genetic diversity and phylogenetic relationships among plum germplasm resources in China assessed with inter-simple sequence repeat markers. *J. Am. Soc Hortic. Sci.* 132, 619–628. doi: 10.21273/JASHS.132.5.619

Luthar, Z. (2018). "Buckwheat genetic resources in central Europe," in *Buckwheat germplasm in the world*. Eds. M. Zhou, I. Kreft, G. Suvorova, Y. Tang and S. H. Woo (Cambridge MA: Academic Press), 127–143.

Manco, R., Basile, B., Capuozzo, C., Scognamiglio, P., Forlani, M., Rao, R., et al. (2019). Molecular and phenotypic diversity of traditional European plum (*Prunus domestica* l.) germplasm of southern Italy. *Sustainability* 11, 14. doi: 10.3390/su11154112

Marchese, A., Bošković, R., Caruso, T., and Tobutt, K. (2010). Intra-allelic variation in introns of the S13-RNase allele distinguishes sweet, wild and sour cherries. *Tree Genet. Genomes.* 6, 963–972. doi: 10.1007/s11295-010-0305-y

Marie, D., and Brown, S. (1993). A cytometric exercise in plant DNA histograms, with 2C values for 70 species. *Bio. Cell.* 78, 41–51. doi: 10.1016/0248-4900(93)01131-S

Mihai, B., Lidija, T., Frida, C., Miljan, C., Viktor, G., Biljana, L., et al. (2012). *The Balkan pomology - plums* (Alnarp: Exactaprinting AB).

Milošević, T., and Milošević, N. (2018). "Plum (Prunus spp.) breeding," in *Advances in plant breeding strategies: Fruits*. Eds. J. M. Al-Khayri, S. M. Jain and D. V. Johnson (Cham, Switzerland: Springer Nature Switzerland AG), 165–215.

Mišić, P. D. (1979). *Šljiva*. Beograd: Nolit.

Mnejja, M., Garcia-Mas, J., Audergon, J.-M., and Arús, P. (2010). Prunus microsatellite marker transferability across rosaceous crops. *Tree Genet. Genomes.* 6, 689–700. doi: 10.1007/s11295-010-0284-z

Mnejja, M., Garcia-Mas, J., Howad, W., Badenes, M. L., and Arús, P. (2004). Simple-sequence repeat (SSR) markers of Japanese plum (*Prunus salicina* lindl.) are highly polymorphic and transferable to peach and almond. *Mol. Ecol. Notes.* 4, 163–166. doi: 10.1111/j.1471-8286.2004.00603.x

Noyer, J. L., Billot, C., Weber, A., Brottier, P., Quero-García, J., and Lebot, V. (2003). "Genetic diversity of taro (Colocasia esculenta (L.) schott) assessed by SSR markers," in *3rd taro symposium*. Eds. L. Guarino, M. Taylor and T. Osborn (Nadi, Fidji Islands: Secretariat of the Pacific Community), 174–180.

Okie, W. R., and Hancock, J. F. (2008). "Plums," in *Temperate fruit crop breeding: Germplasm to genomics*. Ed. J. F. Hancock (Dordrecht: Springer Netherlands), 337–358.

Perrier, X., and Jacquemoud-Collet, J. P. (2006). *DARwin software* (France: Cirad).

Pockar, I. (1982). *Slivarji* (Brezice: Posavski muzej).

Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. *Genetics.* 155, 945. doi: 10.1093/genetics/155.2.945

Reales, A., Sargent, D., Tobutt, K., and Rivera, D. (2009). Phylogenetics of Eurasian plums, *Prunus* l. section *Prunus* (Rosaceae), according to coding and non-coding chloroplast DNA sequences. *Tree Genet. Genomes.* 6, 37–45. doi: 10.1007/s11295-009-0226-9

Rencelj, S., and Pozrl, J. (2015). *Kraski brinjevec in brkinski slivovec* (Sezana: Drustvo proizvajalcev Brkinskega slivovca, Drustvo proizvajalcev Kraskeke brinjevca).

Rouger, R., and Jump, A. (2014). A seascape genetic analysis reveals strong biogeographical structuring driven by contrasting processes in the polyploid saltmarsh species *Puccinellia maritima* and *Triglochin maritima*. *mol. Ecol.* 23, 3158–3170. doi: 10.1111/mec.12802

Rybin, W. A. (1936). Spontane und experimentell erzeugte bastarde zwischen schwarzdorn und kirschpflaume und das abstammungsproblem der kulturpflaume. *Planta.* 25, 22–58. doi: 10.1007/BF01909303

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

Štampar, F., Lešnik, M., Veberič, R., Solar, A., Koron, D., Usenik, V., et al. (2005). *Sadjarstvo* (Ljubljana: Kmečki glas).

Struss, D., Ahmad, R., Southwick, S., and Boritzki, M. (2003). Analysis of sweet cherry (*Prunus avium* l.) cultivars using SSR and AFLP markers. *J. Am. Soc Hortic. Sci.* 128, 904–909. doi: 10.21273/JASHS.128.6.0904

Topp, B. L., Russell, D. M., Michael, N., A., D. M., and Weisheng, L. (2012). "Plum," in *Fruit breeding*. Eds. M. L. Badenes and D. H. Byrne (New York: Springer), 571–621.
Urbanovich, O., Kuzmitskaya, P., and Kilchevsky, A. (2017). Identification and genetic diversity of plum cultivars grown in Belarus. *Russ. J. Genet.* 53, 775–784. doi: 10.1134/S1022795417070134

Urrestarazu, J., Errea, P., Miranda, C., Santesteban, L. G., and Pina, A. (2018). Genetic diversity of Spanish *Prunus domestica* l. germplasm reveals a complex genetic structure underlying. *PloS One* 13, 20. doi: 10.1371/journal.pone.0195591

Urrestarazu, J., Miranda, C., Santesteban, L. G., and Royo, J. B. (2012). Genetic diversity and structure of local apple cultivars from northeastern Spain assessed by microsatellite markers. *Tree Genet. Genomes*. 8, 1163–1180. doi: 10.1007/s11295-012-0502-y

Usenik, V., Stampar, F., and Fajt, N. (2007). Pomological and phenological characteristics of some autochthonous Slovenian plum cultivars. *Acta Hortic*. 734, 53–59. doi: 10.17660/ActaHortic.2007.734.4

Vavilov, N. I. (1949). *The origin, variation, immunity and breeding of cultivated plants* (Waltham, MA: Chronica Botanica).

Vughan, S. P., and Russell, K. (2004). Characterization of novel microsatellites and development of multiplex PCR for large-scale population studies in wild cherry, prunus avium. *Mol. Ecol. Notes*. 4, 429–431. doi: 10.1111/j.1471-8286.2004.00673.x

Vukojevic, D., Simic, J., Dragisic, N., Sevo, D., Misimovic, M., Zavisic, N., et al. (2012). "Evaluation of the quality of autochthonous plum cultivars in the area of

bosanski petrovac,” in *Book of proceedings of the third international scientific symposium, agrosym jahorina*. Ed. D. Kovacevic (Jahorina), 161–166.

Weiss, E., Zohary, D., and Hopf, M. (2012). *Domestication of plants in the old world - the origin and spread of domesticated plants in south-west Asia, Europe, and the Mediterranean basin* (Oxford: Oxford University Press).

Werneck, H. I. (1961). *Die wurzel-und kernechten stammformen der pflaumen in oberösterreich* (Linz: Naturkundliches Jahrbuch der Stadt Linz).

Woldring, H. (2003). On the origin of plums: A study of sloe, damson, cherry plum, domestic plums and their intermediate forms. *Palaeohistoria* 39/40, 535–562.

Wünsch, A. (2009). Cross-transferable polymorphic SSR loci in *Prunus* species. *Sci. Hortic.* 120, 348–352. doi: 10.1016/j.scienta.2008.11.012

Zalokar, J. (1854). *Umno kmetovanje in gospodarstvo* (Ljubljana: C.K. Krajnska kmetijska družba).

Zhebentyayeva, T., Shankar, V., Scorza, R., Callahan, A., Ravelonandro, M., Castro, S., et al. (2019). Genetic characterization of worldwide *Prunus domestica* (plum) germplasm using sequence-based genotyping. *Hortic. Res.* 6, 12. doi: 10.1038/s41438-018-0090-6

Zohary, D. (1992). Is the European plum, *Prunus domestica* l., a p. cerasifera ehrh. x *P. spinosa* l. allo-polyploid? *Euphytica*. 60, 75–77. doi: 10.1007/BF00022260