- Research article
- Open Access
Genetic variation of naturally growing olive trees in Israel: from abandoned groves to feral and wild?
BMC Plant Biologyvolume 16, Article number: 261 (2016)
Naturally growing populations of olive trees are found in the Mediterranean garrigue and maquis in Israel. Here, we used the Simple Sequence Repeat (SSR) genetic marker technique to investigate whether these represent wild var. sylvestris. Leaf samples were collected from a total of 205 trees at six sites of naturally growing olive populations in Israel. The genetic analysis included a multi-locus lineage (MLL) analysis, Rousset’s genetic distances, Fst values, private alleles, other diversity values and a Structure analysis. The analyses also included scions and suckers of old cultivated olive trees, for which the dominance of one clone in scions (MLL1) and a second in suckers (MLL7) had been shown earlier.
The majority of trees from a Judean Mts. population and from one population from the Galilee showed close genetic similarity to scions of old cultivated trees. Different from that, site-specific and a high number of single occurrence MLLs were found in four olive populations from the Galilee and Carmel which also were genetically more distant from old cultivated trees, had relatively high genetic diversity values and higher numbers of private alleles. Whereas in two of these populations MLL7 (and partly MLL1) were found in low frequency, the two other populations did not contain these MLLs and were very similar in their genetic structure to suckers of old cultivated olive trees that originated from sexual reproduction.
The genetic distinctness from old cultivated olive trees, particularly of one population from Galilee and one from Carmel, suggests that trees at these sites might represent wild var. sylvestris. The similarity in genetic structure of these two populations with the suckers of old cultivated trees implies that wild trees were used as rootstocks. Alternatively, trees at these two sites may be remnants of old cultivated trees in which the scion-derived trunk died and was replaced by suckers. However, considering landscape and topographic environment at the two sites this second interpretation is less likely.
The domestication of crop species started 13,000 to 10,000 years before present by gradual selection of desirable traits and of adaptations to agricultural environments . Such artificial selection of individual plants with desirable traits, e.g., high yield, large fruits, loss of shattering seeds, etc., had an artificial selection effect which resulted in genetic differences between crops and their wild ancestors, both in coding and neutral regions of the genome. However, the long co-existence of crops alongside their wild relatives provided opportunities for hybridization, leading to gene flow between the diverging gene pools. Gene flow between cultivated plants and their wild ancestors has been demonstrated in woody species cultivated for their edible fruits such as almonds (Prunus dulcis and P. orientalis) , grapes (Vitis vinifera subsp. vinifera and V. vinifera subsp. sylvestris) [3, 4] and apples (Malus domestica and M. sylvestris) . In addition to gene flow, dispersal of seeds from cultivated trees into natural surroundings can result in feral populations of natural aspect , as shown for several plants introduced to Australia, including Olea europaea [7, 8]. Both these processes can result in substantial difficulties when trying to identify populations as truly wild.
It is generally accepted that the cultivated olive Olea europaea subsp. europaea var. europaea originated from wild var. sylvestris (Mill) Lehr by artificial selection from wild populations . Recently, analysis of plastid DNA diversity among 1,263 supposedly wild olive trees from 108 localities across the Mediterranean area and 534 cultivars suggested that the north Levant (i.e., the area close to the Syrian/Turkish border) was the primary domestication centre of olives . However, one of the earliest indications of the use of olives and possibly also of its cultivation was found in the southeastern Mediterranean area (i.e., in the area of modern Israel) and dated to 6,500 B.C. .
Wild var. sylvestris, often called ‘oleaster’, resembles cultivated olives except for its shrubby growth and smaller leaves and fruits . These characters, however, are highly variable and do not allow reliable distinction between the wild and cultivated varieties. Thus, the identification of olives growing in natural surroundings in the southeast Mediterranean area as var. sylvestris is often questionable . However, using an ecological niche model based on current climatic parameters, Besnard et al.  could identify the natural distribution range of var. sylvestris and could show that current conditions are suitable for its presence in the southwest Levant, i.e., modern Israel.
Studies employing different molecular marker techniques to investigate the relationship between cultivated and wild olives and to map the distribution of wild olives in the Mediterranean area have been conducted before, e.g. [14–22]. In several cases, genetic similarity between trees growing in natural surroundings and cultivated olives was interpreted as evidence for the feral nature (i.e., descended from cultivated trees) of the former [14, 15]. However, the studies by Baldoni et al.  and Belaj et al.  also revealed the existence of genetically distinct populations in Italy and Spain, respectively, which were interpreted as evidence for the continued existence of isolated populations of wild var. sylvestris in the Mediterranean area. Supporting this hypothesis, other studies using DNA [22–24] and allozyme  variation differentiated between cultivated and wild forms of olives. More recently, a comprehensive Bayesian analysis of microsatellite variation that included cultivated and supposedly wild trees from around the Mediterranean Basin showed that wild trees from the southeastern Mediterranean region were genetically closely similar to Spanish cultivars . The study by Diez et al.  as well as others [15, 26] thus suggest that the identification of naturally growing populations of olives as var. sylvestris requires caution in view of the possibility of gene flow between cultivated and wild populations.
In Israel, naturally growing populations of olive trees can be found in the Mediterranean maquis and garrigues of the Carmel and western Galilee mountain ranges. Considering that it is likely that olives have been cultivated continuously in the area for at least 6,000 years [11, 20, 27–30], and that olive groves occupy large parts of the rural landscape, the continued existence of populations of var. sylvestris in the region perhaps is not likely and needs to be studied. Several studies included samples of naturally growing olive trees from the southeastern Mediterranean to infer the distribution and genetic diversity among population of ‘oleaster’ around the Mediterranean [17, 19, 25]. Higher genetic diversity was found in populations of naturally growing olive trees in the west Mediterranean than in the East Mediterranean area, suggesting the existence of genuine var. sylvestris in the west Mediterranean  but questioning the status of naturally growing olive trees in the southeastern Mediterranean. The genetic variation of populations of var. sylvestris potentially could have enormous importance in breeding programs aiming at the introduction of wild alleles conferring valuable traits that were lost during the domestication process . On this background, knowledge of the status of naturally growing populations of olives is of high importance for developing conservation programs for this valuable germplasm. Conservation efforts should also address the risks of hybridization and introgression from domesticated crops into populations of their wild relatives , as recently shown for fruit trees [2, 5].
To determine the identity of naturally growing olive populations in Israel as wild var. sylvestris, feral (var. europaea) or perhaps as abandoned groves, we used SSR markers for the analysis of six naturally growing olive populations sampled at close to far distances from extant cultivated groves. In a previous study we already used a multi-locus lineage (MLL) analysis with the same SSR markers to infer cultivar identity in the same region . We could show the dominance of one clone in scions, and that another clone is frequent in rootstocks of grafted trees. We used these earlier results to assess genetic similarity between supposedly wild populations and local old cultivated olive trees. More specifically, we hypothesize that genuinely wild populations (i.e., var. sylvestris), if such exist, will be genetically different from cultivated olives and feral populations. On this background, a population of naturally growing olives from outside the hypothetical natural distribution range of wild olives in the region  was included as potential reference as a feral population.
Multi-locus lineage analysis
The genetic analysis included 205 naturally growing olive trees sampled in six populations in Israel (Table 1; Fig. 1). Using 15 SSR markers, the number of alleles per locus in the total of 205 trees ranged from four to 32 (Additional file 1: Table S1). Raw microsatellite data for the 15 markers is available in the Additional file 2: Table S2. Analysis of multi-locus genotypes (MLGs) and grouping of MLGs into multi-locus lineages (MLLs) reduced the probability of mistakes resulting from SSR genotyping errors, thus permitting the comparison of naturally growing populations with grafted old olive trees . The multi-locus lineage analysis (Table 2) showed that 10 trees each of the naturally growing and cultivated trees at BGR belonged to the most commonly cultivated regional clone MLL1 . One of the remaining nine trees belonged to MLL7 and the remaining eight trees were assigned to site-specific single occurrence MLLs of which seven were found in naturally growing trees. In ZUR, only three distinct MLLs were found, with 20 trees assigned to the cultivated MLL1, one to MLL7, and only two to a site-specific MLL. In contrast, trees sampled in IDM (Galilee) and the Carmel populations NOR, BOR and OFR generally belonged to site specific MLLs. Whereas both BOR and OFR contained a small number of trees with MLL7, and the two sampled cultivated trees in OFR belonged to MLL1, neither MLL1 nor MLL7 were found in IDM and NOR. Other than these two common MLLs, no MLLs were shared between populations (Table 2 and Additional file 3: Table S3).
In comparison, samples of suckers (collected at the base of tree trunks) and scions, from presumably grafted trees, belonged to 141 and 18 MLLs, respectively (Table 2; ). Of the total of 269 MLLs (Additional file 3: Table S3), 16 were shared by scions and suckers, two were scion specific (MLL10 and 11) and 125 were specific to suckers, the majority of them as single occurrence MLLs (Table 2 and Additional file 3: Table S3).
Genetic diversity estimates
The genetic diversity values (Table 3) showed that values of allelic richness and mean number of private alleles per locus found in IDM and NOR were higher than those found in the other populations. The average number of private alleles per locus in OFR (0.49) was the lowest among the six populations and in comparison to scions (0.53) and suckers (0.57). No noticeable differences among populations and cultivated trees were found in observed and unbiased expected heterozygosity values (Table 3).
Further analysis of private allelic richness indicated that IDM and NOR had the highest number of private alleles per locus (Fig. 2) when corrected for sample size using ADZE . The number of private alleles in 112 different combinations of populations is presented in the Additional file 4: Figure S1. The highest number of private alleles shared by two populations was found in the combination of IDM and NOR (Additional file 4: Figure S1A). In the combination of three populations, the highest number of private alleles was found in IDM together with BOR and NOR, and in the combination of four populations in IDM, BOR, NOR and OFR. The combination of these four populations with suckers yielded the highest number of private alleles per locus (combination of five populations; Additional file 4: Figure S1D). Any combination of suckers from IDM, BOR, NOR and OFR with scion MLLs resulted in smaller numbers of private alleles per locus (Additional file 4: Figure S1A-D). In combinations of two populations with suckers, the number of private alleles per locus was substantially higher in IDM and NOR with suckers than in BOR and OFR with suckers (Additional file 4: Figure S1B).
Genetic differentiation among wild growing populations and cultivated olives
The pairwise Fst analysis (Table 4) revealed that BGR and OFR are most similar to each other (Fst = 0.016), while the highest genetic differentiation was found between populations ZUR and IDM (Fst = 0.061) and between ZUR and NOR and BOR (Fst = 0.052 and 0.051, respectively). BOR and NOR were very similar to IDM (Table 4).
Rousset’s â values (Fig. 3) indicated close genetic similarity between trees from populations ZUR and BGR and MLL1 (-0.26 and -0.22, respectively) and between ZUR and MLL7 (-0.29). Considering that trees of population ZUR were assigned to only three MLLs (MLL1, 7 and 269; Table 2), these results are not surprising. Trees from OFR also showed relatively high similarity with MLL1 (-0.18), whereas NOR, BOR and IDM were found to be most divergent (-0.05, -0.01 and -0.02, respectively; Fig. 3). Comparisons with MLL7 showed a similar pattern, except for individuals from BOR which are relatively more similar to MLL7 (-0.14) than to MLL1.
Population genetic structure
Confirming the MLL analysis, MLGs of scions were found to be fairly homogenous at all given Ks. The dominance of the scion cluster in the cultivated olive individuals sampled in BGR was also evident at all Ks, and a similar genetic structure was also found in trees from ZUR (Fig. 4). Naturally growing trees at BGR showed evidence of admixture and resembled the genetic structure found in suckers (K = 3 in Fig. 4,). At K = 3, the other naturally growing populations from the Galilee (IDM) and Carmel (NOR, BOR and OFR) also showed an admixed genetic structure resembling that found among suckers.
The existence of wild olive trees (Olea europaea subsp. europaea var. sylvestris) in Israel was here investigated by using SSR variation in naturally growing populations. Based on a Bayesian analysis of SSR markers, it has recently been suggested that 38 trees sampled outside cultivated groves in Israel are presumably feral . However, although that analysis included supposedly wild ‘oleaster’ trees that were collected within (Upper and Lower Galilee, Carmel) and outside (Ashkelon, Coastal Plain and Jerusalem) the putative distribution range of var. sylvestris, it did not include the most common cultivars in the southeast Mediterranean region, which renders the conclusions rather hypothetical. Here, a genetic comparison between grafted old olive trees (i.e. 288 scions and 281 suckers) and naturally growing olive trees from different populations in Israel provided more comprehensive information on their identity as feral, cultivated or genuinely wild, and furthermore allowed us to obtain evidence for the possible source of rootstocks of grafted old trees.
Trees at BGR showed close genetic similarity to scions of old olive trees, most strongly to MLL1 (Table 2 and Figs. 3 and 4). MLL7, the most common lineage among rootstocks of grafted old olive trees  was represented once in this population, and diversity values, estimated as allelic richness and mean number of private alleles per locus, were among the lowest of all populations analyzed (Table 3). Considering the dominance of MLL1 in the supposedly cultivated trees of BGR (Table 2) and the similarity of their genetic structure to scions (Fig. 4), our results confirm our a priori assumption of the existence of an abandoned grove at this site.
Similarly but unexpectedly, most of the trees at ZUR (Galilee), where wild var. sylvestris potentially can grow, were very similar to BGR in their genetic composition and genetic structure. A high number of the sampled trees in ZUR belonged to MLL1 (86.9%), two trees were assigned to site specific MLLs, and MLL7 was present in one individual (Table 2). In addition, genetic differentiation between trees from ZUR and MLL1 and MLL7, as measured by Rousset’s value, was the lowest among all pairwise comparisons (Fig. 3), and the genetic diversity values were lowest of all populations investigated (Table 3). Thus, although the ZUR site is best characterized as Mediterranean garrigue (Fig. 1c), and has no resemblance with a grove in terms of tree spacing, traces of former terraces, etc., our results indicate that most trees at this site represent an old abandoned olive grove.
The remaining four populations (IDM, NOR, BOR, OFR) are very different in their genetic composition from BGR and ZUR. Trees from IDM (Galilee) and the three Carmel populations (BOR, NOR and OFR) were found to be similar in terms of relatively high genetic diversity values (Table 3), and genetic differentiation from both BGR and ZUR (abandoned groves) and from MLL1 and MLL7 is high (Fig. 3). However, while the Carmel populations contain large numbers of single occurrence MLLs (≥69% of the total sample size), higher than those found in suckers of cultivated trees (43%; Table 2), 68% of IDM samples belonged to five site-specific MLLs (Table 2), indicating a high frequency of clonal reproduction and/or inbreeding. Indeed, we found indications of inbreeding in IDM (Ho < He, Table 3). This probably can be interpreted as evidence for small effective population size and a high degree of isolation of this population from others. Beyond these peculiarities of IDM, this population and BOR, NOR and OFR fall into two groups. Whereas NOR and IDM do not contain the common scion (MLL1) or rootstock (MLL7) MLLs, both OFR and BOR contain the rootstock MLL7 (Table 2), explaining their similarity to the common rootstock genotype (mean â = -0.10 and -0.14, respectively; Fig. 3), and OFR also contains the scion MLL1. Furthermore, IDM and NOR showed the highest number of private alleles and of alleles that are private to the combination of two populations (Fig. 2 and Additional file 4: Figure S1), and, among these four populations, are most similar to each other (Table 4). The number of private alleles per locus was higher in a combination of suckers, IDM and NOR than in the combination of suckers, BOR and OFR (Additional file 4: Figure S1B). Finally, considering the Structure analysis at K = 3 (Fig. 4), IDM and NOR appear to be more similar to suckers than BOR and OFR in terms of variation of admixed genotypes. In BOR, genotypes with a high proportion of red and green and a low proportion of blue are common, whereas in OFR a relatively high proportion of blue is common. Genotypes are more variably admixed in IDM and NOR. Taking all evidence together, IDM and NOR are most distinct from BGR and ZUR, and BOR and OFR have a somewhat intermediate genetic structure. This in our opinion allows two interpretations:
First, IDM and NOR should be considered wild populations. This interpretation would confirm previous reports that indicated that supposedly wild populations from Carmel and Galilee genetically resemble wild olive populations from Turkey and Syria [17, 19, 21, 26]. As we had demonstrated before  that the majority of old olive trees in the southeastern Mediterranean were maintained by grafting (>80%), the similarity in genetic structure between IDM and NOR on the one hand and suckers of old cultivated olive trees on the other hand would imply that scions were grafted on wild growing olive trees (var. sylvestris). Similarly, a recent genetic survey of scions and rootstocks of old olive trees in the Iberian Peninsula suggested that old olive trees were grafted on wild growing trees . Based on the distances between grafted trees and their spatial arrangement within the groves, Diez et al.  additionally suggested that natural forests were transformed into olive orchards by grafting. In contrast to the situation in Spain, traditional olive groves in the southeastern Mediterranean are found in terraces with equal distances between trees, suggesting that in this region grafting was practiced in the grove itself. This would imply that scions were grafted on saplings that could have been transplanted from outside the grove or germinated in nurseries within it.
Second, the interpretation of the similarity in genetic structure between IDM and NOR on the one hand and suckers of old cultivated olive trees on the other hand can be reversed: such interpretation would imply that the naturally growing trees at IDM and NOR are remnants of old cultivated trees in which the scion-derived trunk died and was replaced by suckers. If this interpretation of IDM and NOR as essentially feral should be correct, populations BOR and OFR may represent an intermediate stage in the transition of orchards into naturally growing populations of feral origin. However, as naturally growing trees at NOR and IDM grow in conditions that are unsuitable for agriculture (c.f. Fig. 1b), it seems more likely to us that they represent wild than abandoned cultivated trees taken over by their suckers.
The comparison of naturally growing olive tree populations with MLL genotypes of scions and suckers of old cultivated olive trees in Israel allowed us to assess the status of naturally growing populations as abandoned groves, feral or wild var. sylvestris. The interpretation of two of six populations analyzed as wild var. sylvestris implies that grafting in the past used wild plants as rootstocks. In an area where olive cultivation has a history of several thousand years, it is astonishing to have identified naturally growing olive tree populations which are partly well-differentiated from cultivated plants. Considering the high abundance of cultivated and feral olive trees in the region, conservation of this valuable genetic material is of greatest importance.
The studied populations
Naturally growing olive trees can be found sparsely in Israel in natural habitats surrounding cultivated groves and residential areas. Surveys were conducted in the Carmel and Galilee to locate populations of at least 20 olive trees (O. europaea subsp. europaea) growing in natural surroundings and resembling the shrubby phenotype of southeastern Mediterranean var. sylvestris; identification of plants was based on the Analytical Flora of Israel  and Flora Palaestina  and done by Dr. Ori Fragman-Sapir (Head Scientist, Jerusalem Botanical Gardens). Naturally growing olive (var. sylvestris) is not included in the Red List of the Israeli Flora . Nevertheless, sampling in natural reserves was coordinated and approved by the Israel Nature and National Parks Protection Authority (license no. 2014/40360).
Five sites were selected (Table 1; Fig. 1a) within the species’ hypothetical natural distribution range in the region (Galilee and Carmel Mts.), and one population was sampled in the Judean Mts. outside the hypothetical natural distribution range of wild olive . At all sites, olive trees were growing at uneven distances which is untypical for cultivated groves. Two populations were sampled in the Galilee, three in the Carmel Mountain range and one in the Judean Mts. (Table 1; Fig. 1a). Number of samples collected is related to population size. Where possible, trees were sampled randomly along widely spaced transects through the collecting areas in order to represent their genetic diversity; this resulted in a sample of altogether 205 trees. As the Carmel Mountain range is considered the southern limit of the distribution range of wild olive trees in Israel , the populations in the Galilee and the Carmel regions potentially may represent wild var. sylvestris.
At Idmit (Galilee; IDM), the sampled population grows remotely from contemporary olive groves (Additional file 6: Figure S3) as a dense stand on the edge of a cliff with a southwest slope facing the Mediterranean Sea (335 m above sea level; a.s.l). Trees at IDM face strong herbivore pressure, mainly by rock-hyrax (Procavia capensis), and grow in patches as small shrubs (Fig. 1b), different from the trees at the other sites. At IDM, 25 trees were sampled. At the second site, sampled in the western Galilee (ZUR; Zurit, 295 m a.s.l.), closer to an olive cultivation area (~3 km), trees were observed in scattered patches across a large area of about 50 hectares (Fig. 1c). In order to represent the genetic diversity of this population, 23 samples were collected from across the entire area. In the Carmel region, 128 trees were sampled in three locations: (1) on south and north facing slopes of the western part of Nachal Oren (NOR, 140 m a.s.l., 35 trees), a natural habitat that is not suitable for agriculture and has long been used as a model site for biodiversity and speciation studies ; (2) at a higher point of the Carmel Mountain in the vicinity of Beit Oren (BOR, 385 m a.s.l., 54 trees); (3) in the southern part of the Carmel range (Ofer), trees were sampled on north and south facing slopes of the hill (OFR; 170 m a.s.l., 39 trees). At the last site, trees are distributed irregularly, untypical for olive groves, and tree appearance is different from nearby cultivated trees found at the edge of a rural residential area. Hypothesizing that the OFR population might be feral, we included two cultivated trees from the residential area in our analysis; OFR is also relatively close to an old olive grove (32° 37′ 48.00″N, 35° 0′ 0.00″E) which had been sampled in our previous study  (Additional file 6: Figure S3). At Bar Giora (BGR; Judean Mts.; 450 m a.s.l.), trees were sampled in a natural reserve on northeast facing slope terraces in the remains of an abandoned orchard and thus were assumed to represent a population of non-wild olives. Of the total of 29 trees sampled here, 11 resembled cultivated trees; however, since unambiguous distinction of cultivated from naturally growing trees was not possible, the 29 trees were treated as one group. The Galilee (IDM, ZUR) and Carmel (NOR, OFR) areas can be characterized as garrigue (Fig. 1); olive trees in BOR and BGR grow in Mediterranean dense forest (maquis).
DNA was extracted from leaf samples using the Invisorb Plant Mini Kit (Invitek) following the manufacturer’s protocol. Simple Sequence Repeat (SSR) markers used in olive trees [40–48] had previously been screened  resulting in the use of 15 markers with PCR conditions previously described (Additional file 1: Table S1) . SSR products were separated at the Center of Genomic Technologies (The Hebrew University of Jerusalem) on an ABI automated sequencer (Applied Biosystems) as a multiplex of several loci labeled with three different fluorescent dyes (6-FAM, NED and HEX; Applied Biosystems). Electropherograms were scored manually using Genmarker 1.75 (SoftGenetics, State College, Pennsylvania, USA).
Analysis of multi-locus genotypes (MLGs) and grouping of MLGs into multi-locus lineages (MLLs) was done using Genotype 1.2 . To estimate diversities, one sample of the most common MLG of each MLL from each population was used. Genetic diversity values calculated included observed (Ho) and unbiased expected (uHe) heterozygosity using GenAlEx v6.5 [50, 51]. To better account for differing sample sizes, allelic and private allelic richness was calculated with a rarefaction approach using the Allelic Diversity Analyzer ADZE software . ADZE was also used to calculate the number of alleles private to combinations of populations and scions and suckers of old cultivated trees. The diversity values obtained were compared with data from 281 suckers collected at the base of tree trunks and 288 scions of cultivated old olive trees  using the same 15 SSR loci as used here. The old cultivated olive trees were sampled in 32 groves in the southeastern Mediterranean; the location of the groves is listed in Barazani et al. . The ZUR population included only three different MLLs and was excluded from analyses based on MLLs. In addition, Fst values were used to estimate genetic distances among populations, using GenAlEx v6.5. Rousset’s genetic distances (â)  were estimated using Spagedi 1.4c  to compare genetic distances among individuals within and between populations as well as genetic distances to MLL1 and MLL7, the two most commonly cultivated clones in the region .
Structure V.2.3.4  was used for Bayesian clustering with the admixture model to assign each MLG from each of the studied naturally growing populations and from old cultivated trees to K clusters. According to the recommendation by Pritchard et al. , 10 independent runs for given Ks (2 to 8) were performed with a burn-in length of 10,000, followed by 20,000 repetitions. The log likelihoods for a given K were used to choose the best given K based on an ad hoc quantity of ∆K .
Naturally growing olive population at Bar Giora
Naturally growing olive population at Beit Oren
Naturally growing olive population at Idmit
Naturally growing olive population at Nachal Oren
Naturally growing olive population at Ofer
Simple sequence repeat
Naturally growing olive population at Zurit
Purugganan MD, Fuller DQ. The nature of selection during plant domestication. Nature. 2009;457:843–8.
Delplancke M, Alvarez N, Espindola A, Joly H, Benoit L, Brouck E, Arrigo N. Gene flow among wild and domesticated almond species: insights from chloroplast and nuclear markers. Evol Appl. 2012;5:317–29.
De Andres MT, Benito A, Perez-Rivera G, Ocete R, Lopez MA, Gaforio L, Munoz G, Cabello F, Zapater JMM, Arroyo-Garcia R. Genetic diversity of wild grapevine populations in Spain and their genetic relationships with cultivated grapevines. Mol Ecol. 2012;21:800–16.
Di Vecchi-Staraz M, Laucou V, Bruno G, Lacombe T, Gerber S, Bourse T, Boselli M, This P. Low level of pollen-mediated gene flow from cultivated to wild grapevine: consequences for the evolution of the endangered subspecies Vitis vinifera L. subsp silvestris. J Hered. 2009;100:66–75.
Cornille A, Feurtey A, Gelin U, Ropars J, Misvanderbrugge K, Gladieux P, Giraud T. Anthropogenic and natural drivers of gene flow in a temperate wild fruit tree: a basis for conservation and breeding programs in apples. Evol Appl. 2015;8:373–84.
Warwick SI, Beckie HJ, Hall LM. Gene flow, invasiveness, and ecological impact of genetically modified crops. Ann Ny Acad Sci. 2009;1168:72–99.
Spennemann DHR, Allen LR. Feral olives (Olea europaea) as future woody weeds in Australia: a review. Aust J Exp Agr. 2000;40:889–901.
Crossman ND, Bass DA, Virtue JG, Jupp PW. Feral olives (Olea europaea L.) in southern Australia: an issue of conservation concern. Adv Hortic Sci. 2002;16:175–83.
Zohary D, Spiegel-Roy P. Beginnings of fruit growing in the old world. Science. 1975;187:319–27.
Besnard G, Khadari B, Navascues M, Fernandez-Mazuecos M, El Bakkali A, Arrigo N, Baali-Cherif D, Brunini-Bronzini de Caraffa V, Santoni S, Vargas P, et al. The complex history of the olive tree: from Late Quaternary diversification of Mediterranean lineages to primary domestication in the northern Levant. P Roy Soc B. 2013;280:20122833.
Galili E, Stanley D, Sharvit J, Weinstein-Evron M. Evidence for earliest olive-oil production in submerged settlements off the Carmel coast, Israel. J Archeol Sci. 1997;24:1141–50.
Lavee S, Zohary D. The potential of genetic diversity and the effect of geographically isolated resources in olive breeding. Israel J Plant Sci. 2011;59:3–13.
Green PS. A revision of Olea L. (Oleaceae). Kew Bull. 2002;57:91–140.
Baldoni L, Tosti N, Ricciolini C, Belaj A, Arcioni S, Pannelli G, Germana MA, Mulas M, Porceddu A. Genetic structure of wild and cultivated olives in the central Mediterranean basin. Ann Bot-London. 2006;98:935–42.
Belaj A, Munoz-Diez C, Baldoni L, Porceddu A, Barranco D, Satovic Z. Genetic diversity and population structure of wild olives from the North-Western Mediterranean assessed by SSR markers. Ann Bot-London. 2007;100:449–58.
Belaj A, Satovic Z, Rallo L, Trujillo I. Genetic diversity and relationships in olive (Olea europaea L.) germplasm collections as determined by randomly amplified polymorphic DNA. Theor Appl Genet. 2002;105:638–44.
Breton C, Tersac M, Berville A. Genetic diversity and gene flow between the wild olive (oleaster, Olea europaea L.) and the olive: several Plio-Pleistocene refuge zones in the Mediterranean basin suggested by simple sequence repeats analysis. J Biogeogr. 2006;33:1916–28.
De Caraffa VB, Maury J, Gambotti C, Breton C, Berville A, Giannettini J. Mitochondrial DNA variation and RAPD mark oleasters, olive and feral olive from Western and Eastern Mediterranean. Theor Appl Genet. 2002;104:1209–16.
Lumaret R, Ouazzani N, Michaud H, Vivier G, Deguilloux MF, Di Giusto F. Allozyme variation of oleaster populations (wild olive tree) (Olea europaea L.) in the Mediterranean Basin. Heredity. 2004;92:343–51.
Zohary D, Hopf M, Weiss E. Domestication of plants in the old world. 4th ed. Oxford: Oxford Science Publications, Claredon Press; 2012.
Lumaret R, Ouazzani N. Plant genetics: Ancient wild olives in Mediterranean forests. Nature. 2001;413:700.
Vargas P, Kadereit J. Molecular fingerprinting evidence (ISSR, inter-simple sequence repeats) for a wild status of Olea europaea L. (Oleaceae) in the Eurosiberian North of the Iberian Peninsula. Flora. 2001;196:142–52.
Angiolillo A, Mencuccini M, Baldoni L. Olive genetic diversity assessed using amplified fragment length polymorphisms. Theor Appl Genet. 1999;98:411–21.
Besnard G, Baradat P, Breton C, Khadari B, Berville A. Olive domestication from structure of oleasters and cultivars using nuclear RAPDs and mitochondrial RFLPs. Genet Sel Evol. 2001;33:S251–68.
Diez CM, Trujillo I, Martinez-Urdiroz N, Barranco D, Rallo L, Marfil P, Gaut BS. Olive domestication and diversification in the Mediterranean Basin. New Phytol. 2015;206:436–47.
Breton C, Pinatel C, Medail F, Bonhomme F, Berville A. Comparison between classical and Bayesian methods to investigate the history of olive cultivars using SSR-polymorphisms. Plant Sci. 2008;175:524–32.
Neumann FH, Kagan EJ, Schwab M, Stein M. Palynology, sedimentology and palaeoecology of the late Holocene Dead Sea. Quaternary Sci Rev. 2007;26:1476–98.
Neumann FH, Kagan EJ, Leroy SAG, Baruch U. Vegetation history and climate fluctuations on a transect along the Dead Sea west shore and their impact on past societies over the last 3500 years. J Arid Environ. 2010;74:756–64.
Liphschitz N, Gophna R, Hartman M, Biger G. The beginning of olive (Olea europaea) cultivation in the Old World: a reassessment. J Archaeol Sci. 1991;18:441–53.
Kaniewski D, Paulissen E, van Campo E, Bakker J, van Lerberghe K, Waelkens M. Wild or cultivated Olea europaea L. in the eastern Mediterranean during the middle–late Holocene? A pollen-numerical approach. The Holocene. 2009;19:1039–47.
Rugini E, De Pace C, Gutie’rrez-Pesce P, Muleo R. Olea. In: Kole C, editor. Wild Crop Relatives: Genomic and Breeding Resources: Temperate Fruits. Berlin: Springer; 2011. p. 79–117.
Ellstrand NC, Prentice HC, Hancock JF. Gene flow and introgression from domesticated plants into their wild relatives. Annu Rev Ecol Syst. 1999;30:539–63.
Barazani O, Westberg E, Hanin N, Dag A, Kerem Z, Tugendhaft Y, Hmidat M, Hijawi T, Kadereit JW. A comparative analysis of genetic variation in rootstocks and scions of old olive trees - a window into the history of olive cultivation practices and past genetic variation. BMC Plant Biol. 2014;14:146.
Szpiech ZA, Jakobsson M, Rosenberg NA. ADZE: a rarefaction approach for counting alleles private to combinations of populations. Bioinformatics. 2008;24:2498–504.
Diez CM, Trujillo I, Barrio E, Belaj A, Barranco D, Rallo L. Centennial olive trees as a reservoir of genetic diversity. Ann Bot-London. 2011;108:797–807.
Feinbrun-Dothan N, Danin A. Analytical Flora of Israel. Jerusalem, Israel: CANA Publishing House Ltd.; 1998.
Zohary M. Flora Palaestina vol. Part 1. The Israel Academy of Sciences and Humanities: Jerusalem; 1966.
Shmida A, Fragman O, Nathan R, Shamir Z, Sapir Y. The red plants of Israel: A proposal of updated and revised list of plant species protected by the law. Ecol Mediterr. 2002;28:55–64.
Nevo E. “Evolution canyon”: A microcosm of life’s evolution focusing on adaptation and speciation. Isr J Ecol Evol. 2006;52:485–506.
Belaj A, Cipriani G, Testolin R, Rallo L, Trujillo I. Characterization and identification of the main Spanish and Italian olive cultivars by simple-sequence-repeat markers. Hortscience. 2004;39:1557–61.
Carriero F, Fontanazza G, Cellini F, Giorio G. Identification of simple sequence repeats (SSRs) in olive (Olea europaea L.). Theor Appl Genet. 2002;104:301–7.
Cipriani G, Marrazzo M, Marconi R, Cimato A, Testolin R. Microsatellite markers isolated in olive (Olea europaea L.) are suitable for individual fingerprinting and reveal polymorphism within ancient cultivars. Theor Appl Genet. 2002;104:223–8.
De La Rosa R, James CM, Tobutt KR. Isolation and characterization of polymorphic microsatellites in olive (Olea europaea L.) and their transferability to other genera in the Oleaceae. Mol Ecol Notes. 2002;2:265–7.
Diaz A, De la Rosa R, Martin A, Rallo P. Development, characterization and inheritance of new microsatellites in olive (Olea europaea L.) and evaluation of their usefulness in cultivar identification and genetic relationship studies. Tree Genet Genomes. 2006;2:165–75.
La Mantia M, Lain O, Caruso T, Testolin R. SSR-based DNA fingerprints reveal the genetic diversity of Sicilian olive (Olea europaea L.) germplasm. J Hortic Sci Biotech. 2005;80:628–32.
Rallo P, Dorado G, Martin A. Development of simple sequence repeats (SSRs) in olive tree (Olea europaea L.). Theor Appl Genet. 2000;101:984–9.
Saumitou-Laprade P, Vassiliadis C, Epplen JT, Hardt C. Isolation of microsatellite loci for paternity testing in Phillyrea angustifolia L. (Oleaceae). Mol Ecol. 2000;9:112–4.
Sefc KM, Lopes S, Mendonca D, Dos Santos MR, Machado MLD, Machado AD. Identification of microsatellite loci in olive (Olea europaea) and their characterization in Italian and Iberian olive trees. Mol Ecol. 2000;9:1171–3.
Meirmans PG, van Tienderen PH. GENOTYPE and GENODIVE: Two programs for the analysis of genetic diversity of asexual organisms. Mol Ecol Notes. 2004;4:792–4.
Peakall R, Smouse PE. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6:288–95.
Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics. 2012;28:2537–9.
Rousset F. Genetic differentiation between individuals. J Evolution Biol. 2000;13:58–62.
Hardy O, Vekemans X. SPAGeDi: A versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol Ecol Notes. 2002;2:618–20.
Hubisz MJ, Falush D, Stephens M, Pritchard JK. Inferring weak population structure with the assistance of sample group information. Mol Ecology Resour. 2009;9:1322–32.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.
Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14:2611–20.
We thank Mrs. Michal Barzilai, Mr. Isaac Zipori and other colleagues in Israel and abroad who helped at different stages of this study. Helpful comments by an anonymous reviewer are gratefully acknowledged.
This study was supported by the German Research Foundation’s (DFG) trilateral program (Grant no. KA 635/14).
Availability of data and materials
All relevant data supporting our findings is provided in the article and supporting information. Results of our previous SSR analysis of grafted old olive trees can be found in BMC Plant Biol 2014, 14:146.
OB, AD, ZK, and JWK conceived this study. OB, AD, EW, NH, OFS, ZK and YT mapped olive populations and collected the samples. NH did the laboratory work; EW, AKK, GBA and NH analyzed the data. OB and JWK wrote the manuscript with contributions from all co-authors. All co-authors approved submission to BMC Plant Biology. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Olea europaea subsp. europaea var. sylvestris is not included in the Red List of the Israeli flora. Sampling was approved and conducted in accordance with the regulations of the Israel Nature and National Parks Protection Authority (license no. 2014/40360).
SSR markers used, their expected size range, repeated motives and number of alleles found in naturally growing olive populations. Raw microsatellite data is available and enclosed as Additional file 2: Table S2. (PDF 188 kb)
Raw microsatellite data. The fragment sizes (in base pairs) of the two alleles per individual for each locus are given as a and b (0 represents missing data). (XLSX 38 kb)
Number of olive trees assigned to different multi-locus lineages (MLLs). (XLSX 18 kb)
Number of private alleles per locus in combinations of populations. A to D present values for the combination of two to five populations (treating scions and suckers of old olive trees as populations). (PDF 217 kb)
∆K values for the different Ks were calculated according to Evanno et al. , showing that K = 3 is the optimal K for the Structure analysis. (PDF 69 kb)