# Primary Selection and Secondary Diversification: Two Key Processes in the History of Olive Domestication

Bouchaib Khadari, Ahmed El Bakkali

► To cite this version:

Bouchaib Khadari, Ahmed El Bakkali. Primary Selection and Secondary Diversification: Two Key Processes in the History of Olive Domestication. International Journal of Agronomy and Agricultural Research, 2018, 2018, 10.1155/2018/5607903. hal-02620833

HAL Id: hal-02620833

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

Submitted on 26 May 2020

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
Research Article

# Primary Selection and Secondary Diversification: Two Key Processes in the History of Olive Domestication

Bouchaib Khadari<sup>1,2</sup> and Ahmed El Bakkali<sup>3</sup>

Correspondence should be addressed to Bouchaib Khadari; khadari@supagro.fr

Received 9 January 2018; Revised 14 May 2018; Accepted 9 August 2018; Published 29 November 2018

Academic Editor: Faruk Toklu

Copyright © 2018 Bouchaib Khadari and Ahmed El Bakkali. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Knowledge on the crop domestication process is important from a cultural and agricultural standpoint since it can shed light on the origin and history of human civilizations as well as the management of genetic resources, while offering guidance for modern breeding. The olive tree (Olea europaea ssp. europaea) is the most iconic of the old crop species of the Mediterranean Basin (MB). Primary domestication from wild olive probably occurred around 6000 BP in the Middle East. However, the question remains as to whether cultivated olive derived from a single domestication event in the Levant, followed by secondary diversification, or whether it was the result of independent domestication events. Here, we analyzed a comprehensive sample collected from 35 wild populations (722 individuals) and 410 cultivars from across the MB using nuclear and plastid DNA markers. Our genetic investigations argue in favor of a single primary domestication event in the eastern MB, followed by diffusion of the first domesticated olive and diversification in the central and western MB as key processes in the olive tree history.

## 1. Introduction

Understanding crop domestication and diversification processes is important to infer the origin of the crop and highlight the history of human civilizations. These investigations can be useful for genetic resource management while offering guidance for modern breeding. Olive (Olea europaea ssp. europaea) is considered to be the most iconic tree in Mediterranean areas. Oliviculture is one of the oldest cropping practices developed in these areas, and olive trees have therefore accompanied the emergence of early Mediterranean civilizations [1]. According to paleobotanical, archaeological, and genetic investigations, the olive tree may have persisted around the Mediterranean Basin (MB) as part of the natural plant community since the Late Tertiary [2]. However, despite the economic, cultural, ecological, and historical importance of the species, its origin and history have yet to be clearly documented. Clarifying the olive domestication and diversification process has therefore long been a focus of active scientific research [3].

According to Carrión et al. [4], during the Middle and Late Pleniglacial (59,000–11,500 yrs BP), Olea europaea had persisted in three thermophilous refugia located in the southern areas of the north MB, the southern Levant, and North Africa. Due to outcrossing of olive species, wild forms (oleasters, O. europaea ssp. europaea var. sylvestris) probably contributed to genetic diversity at the local scale, thus facilitating secondary diversification. It is usually considered that the center of primary olive domestication, from wild progenitors, began roughly around 6000 BP in the Middle East, near the border between Turkey and Syria [1, 5]. This assumption was supported by investigations based on chloroplast DNA polymorphism [5, 6], showing that more than 90% of olive cultivars across the Mediterranean Basin share the same eastern-like haplotype, therefore indicating an east-west human-mediated diffusion of cultivars in the MB. However, the scenario of a single primary domestication center has yet to be demonstrated. Based on paleobotanical and archaeological investigations, early

<footnote n="1"><sup>1</sup>AGAP, University Montpellier, CIRAD, INRA, Montpellier SupAgro, Montpellier, France</footnote>

<footnote n="2"><sup>2</sup>Conservatoire Botanique National Méditerranéen de Porquerolles (CBNMed), UMR AGAP, Montpellier, France</footnote>

<footnote n="3"><sup>3</sup>INRA, UR Amélioration des Plantes et Conservation des Ressources Phytogénétiques, Meknès, Morocco</footnote>
exploitation and use of wild olive trees from the Near East to Spain have been documented since the Neolithic period [7, 8]. The findings of several studies support multiple origins of cultivars across the Mediterranean region [9–13], but it remains unclear whether this reflects secondary diversification or multiple independent primary domestication events [13–15]. Indeed, genetic patterns observed by Diez et al. [13] suggest the occurrence of a second separate olive domestication event in the central MB. The hypothesis of a second independent domestication in the central MB remains to be explored because cultivated olive diversification may also have occurred in this area and in the western MB as the result of local independent domestication [14]. Following the response letter of Diez and Gaut [15], the hypothesis of an independent domestication event in the central MB may ultimately be confirmed or disproven. Therefore, the question as to whether there was a single olive domestication center or multiple ones has yet to be answered.

In the present study, we investigated the history of olive trees through a comprehensive sampling of genuinely wild populations and domesticated forms from across the Mediterranean using nuclear and plastid marker analysis. Our results are discussed based on the question as to whether the current cultivated olive derived from a single domestication event in the Levant followed by secondary diversification or whether it is the result of independent domestication events.

## 2. Materials and Methods

2.1. *Plant Material.* A total of 1,132 distinct genotypes were analyzed in this study, including 410 cultivars from 15 countries and 722 wild olive samples from 35 populations throughout the Mediterranean area (Figure 1; Tables S1 and S2). All cultivars are maintained in the *ex situ* Worldwide Olive Germplasm Bank at the experimental station of Tassaout, INRA, Marrakech, Morocco [6, 16]. Wild populations were sampled in natural areas far from olive agroecosystems in order to minimize admixture with cultivated olives, while taking morphological traits that differ from those of cultivated olives into consideration, such as smaller fruits with less fleshy mesocarp [17]. Thirty of the 35 wild populations were previously described and analyzed [5, 18].

2.2. *Molecular Analysis.* Total DNA was extracted from 100 mg of fresh leaf tissue, as described by Khadari et al. [19]. DNA quality was checked on 1% agarose gel, and the concentration was estimated using spectrofluorometry (GENios Plus, TECAN, Grödig, Austria).

A set of 16 nuclear microsatellite loci was chosen and used for sample genotyping [20–23] (Table 1). PCR amplification and product separation using a 3130XL capillary sequencer (Applied Biosystems, Foster City, CA, USA) were conducted as described by El Bakkali et al. [16]. Plastid DNA was characterized using 37 polymorphic simple sequence repeat (SSR) loci and two cleaved amplified polymorphism sites (CAPS-XapI and CAPS-EcoRI), as described by Besnard et al. [24].

2.3. *Data Analysis.* We computed the following genetic diversity parameters for wild and cultivated olives separately and for each genetic group: the number of alleles (Na), expected heterozygosity (He), and observed heterozygosity (Ho) using the Excel Microsatellite Toolkit v3.1 [25]. The inbreeding coefficient (Fis) was calculated using the $F_{STAT}$ program v2.9.3.2b [26], whereas the allelic richness (Ar; [27]) was estimated using the ADZE program [28]. The Mann-Whitney comparison test was used to evaluate the significance of the allelic richness differences.

To investigate the genetic structure pattern within olive samples (wild and cultivated), discriminant analysis of principal components (DAPC; [29]) with the *ADEGENET* 1.3.1 package [30] in the R environment was applied with a *priori* grouping assumptions based on previous studies [6, 13, 18]. Unlike the *STRUCTURE* program, the absence of any assumption about the underlying population genetics model, in particular concerning Hardy-Weinberg equilibrium or linkage equilibrium, is one of the main assets of DAPC [29]. Based on the model-based Bayesian clustering approach implemented in the *STRUCTURE* program [31] as described in previous studies on olive species [6, 13, 18], wild olive was found to be structured in two groups (named western-central and eastern Mediterranean wild), whereas cultivated olive was in three groups (called western, central, and eastern Mediterranean cultivated olive), with one group shared between wild and cultivated olives (eastern Mediterranean). Hence, we set an *a priori* group number of four in the DAPC method for the whole dataset.

Moreover, once the wild and cultivated olive genotypes were assigned to their a posteriori genetic groups, relationships among genotypes and genetic groups were analyzed by principal coordinate analysis (PCoA) based on the simple matching coefficient [32], as implemented in the *DARWIN* v. 6.0.11 program [33]. Pairwise genetic differentiation and significance ($F_{ST}$, [34]) between genetic groups, as revealed by membership assignation using the DAPC method, was estimated using 100,000 permutations with the *GENEPOP* program [35], and the unrooted $F_{ST}$ was plotted using the *POPTREE2* program [36], with 999 bootstrap replicates with the Neighbor-joining method.

## 3. Results

3.1. *Genetic Diversity in Wild and Cultivated Olive.* Based on the analysis of 1,132 genotypes using 16 SSR markers, we identified a total of 427 alleles with an average of 26.69 alleles per locus. The number of alleles observed in wild olive (420) was higher than that in cultivated olive (276). Similarly, the expected heterozygosity (He, diversity index) was greater in wild than in cultivated olive (Table 1).

A total of 33 plastid haplotypes were identified in both wild and cultivated olives. More plastid haplotypes belonging to E1 were identified compared to E2 and E3, i.e., 18, 10, and 4 haplotypes, respectively (Table 2). Otherwise, we revealed more maternal lineages in wild (32) than in cultivated olive (12). In fact, the highest proportion of plastid haplotypes in cultivated olive was observed for E1.1 (79.9%) followed by E1.2 (8.1%), whereas those belonging to E2 and
<figure><img src="image_2.png" /><figcaption>FIGURE 1: Map indicating the location, number of olive samples, and proportion of maternal lineages according to the geographic origin for wild and cultivated olives. Cultivars are shown in green and wild olives in red.</figcaption></figure>

<table><caption>TABLE 1: Genetic parameters of the 16 SSR markers used in this study for wild and cultivated olives.</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Locus</th><th rowspan="2">Size (bp)</th><th colspan="5">Wild (total of 722)</th><th colspan="5">Cultivated olive (total of 410)</th></tr><tr><th>Na</th><th>Ho</th><th>He</th><th>Fis</th><th>Ar</th><th>Na</th><th>Ho</th><th>He</th><th>Fis</th><th>Ar</th></tr></thead><tbody><tr><td>1</td><td>DCA01<sup>(1)</sup></td><td>205-275</td><td>35</td><td>0.47</td><td>0.72</td><td>0.34***</td><td>31.9</td><td>205-275</td><td>21</td><td>0.71</td><td>0.62</td><td>-0.15ns</td><td>20.9</td></tr><tr><td>2</td><td>DCA03<sup>(1)</sup></td><td>222-277</td><td>26</td><td>0.72</td><td>0.81</td><td>0.11***</td><td>25.1</td><td>228-265</td><td>14</td><td>0.91</td><td>0.85</td><td>-0.06ns</td><td>14</td></tr><tr><td>3</td><td>DCA04<sup>(1)</sup></td><td>117-196</td><td>35</td><td>0.72</td><td>0.92</td><td>0.22***</td><td>33.6</td><td>117-200</td><td>33</td><td>0.65</td><td>0.85</td><td>0.23***</td><td>33</td></tr><tr><td>4</td><td>DCA05<sup>(1)</sup></td><td>187-245</td><td>23</td><td>0.82</td><td>0.88</td><td>0.08***</td><td>20.4</td><td>192-214</td><td>12</td><td>0.5</td><td>0.5</td><td>0.01ns</td><td>12</td></tr><tr><td>5</td><td>DCA08<sup>(1)</sup></td><td>118-185</td><td>32</td><td>0.91</td><td>0.93</td><td>0.03***</td><td>29.3</td><td>124-167</td><td>21</td><td>0.92</td><td>0.83</td><td>-0.11*</td><td>21</td></tr><tr><td>6</td><td>DCA09<sup>(1)</sup></td><td>163-225</td><td>29</td><td>0.87</td><td>0.93</td><td>0.06***</td><td>27.5</td><td>163-219</td><td>25</td><td>0.95</td><td>0.89</td><td>-0.07ns</td><td>25</td></tr><tr><td>7</td><td>DCA11<sup>(1)</sup></td><td>124-209</td><td>40</td><td>0.74</td><td>0.82</td><td>0.11***</td><td>37.2</td><td>126-201</td><td>23</td><td>0.88</td><td>0.83</td><td>-0.06*</td><td>23</td></tr><tr><td>8</td><td>DCA14<sup>(1)</sup></td><td>170-194</td><td>25</td><td>0.67</td><td>0.9</td><td>0.25***</td><td>24.2</td><td>170-194</td><td>15</td><td>0.74</td><td>0.71</td><td>-0.04*</td><td>15</td></tr><tr><td>9</td><td>DCA15<sup>(1)</sup></td><td>243-266</td><td>8</td><td>0.64</td><td>0.74</td><td>0.13***</td><td>7.81</td><td>243-266</td><td>6</td><td>0.7</td><td>0.64</td><td>-0.09ns</td><td>5.99</td></tr><tr><td>10</td><td>DCA16<sup>(1)</sup></td><td>122-235</td><td>48</td><td>0.57</td><td>0.95</td><td>0.40***</td><td>45.3</td><td>124-229</td><td>34</td><td>0.95</td><td>0.87</td><td>-0.08ns</td><td>33.9</td></tr><tr><td>11</td><td>DCA18<sup>(1)</sup></td><td>155-205</td><td>26</td><td>0.88</td><td>0.91</td><td>0.03***</td><td>24.7</td><td>157-191</td><td>17</td><td>0.91</td><td>0.83</td><td>-0.09ns</td><td>17</td></tr><tr><td>12</td><td>EMO90<sup>(2)</sup></td><td>183-216</td><td>16</td><td>0.4</td><td>0.63</td><td>0.37***</td><td>15.4</td><td>185-212</td><td>9</td><td>0.71</td><td>0.67</td><td>-0.05ns</td><td>9</td></tr><tr><td>13</td><td>GAPU59<sup>(3)</sup></td><td>206-251</td><td>22</td><td>0.77</td><td>0.84</td><td>0.08*</td><td>19.4</td><td>208-240</td><td>11</td><td>0.59</td><td>0.62</td><td>0.04***</td><td>11</td></tr><tr><td>14</td><td>GAPU71A<sup>(3)</sup></td><td>207-270</td><td>25</td><td>0.66</td><td>0.72</td><td>0.09***</td><td>23.2</td><td>207-252</td><td>14</td><td>0.54</td><td>0.46</td><td>-0.17**</td><td>14</td></tr><tr><td>15</td><td>GAPU 71B<sup>(3)</sup></td><td>132-168</td><td>13</td><td>0.65</td><td>0.69</td><td>0.07***</td><td>12.7</td><td>136-166</td><td>9</td><td>0.9</td><td>0.81</td><td>-0.12**</td><td>8.99</td></tr><tr><td>16</td><td>UDO99-36<sup>(4)</sup></td><td>134-170</td><td>17</td><td>0.76</td><td>0.85</td><td>0.11***</td><td>15.7</td><td>140-168</td><td>12</td><td>0.69</td><td>0.71</td><td>0.023ns</td><td>12</td></tr><tr><td></td><td>Mean</td><td></td><td></td><td>26.3</td><td>0.7</td><td>0.83</td><td>0.15***</td><td>24.6</td><td></td><td>17.3</td><td>0.77</td><td>-0.05***</td><td>17.2</td></tr><tr><td></td><td>Total</td><td></td><td></td><td>420</td><td></td><td></td><td></td><td></td><td></td><td>276</td><td></td><td></td><td></td></tr></tbody></table>

<footnote>(1)Sefc et al. [20], (2)De La Rosa et al. [21], (3)Carriero et al. [22], (4)Cipriani et al. [23]. Significance level at $p \le 0.05^*$, 0.01^**, and 0.001^***. ns: nonsignificant.</footnote>
<table><caption>TABLE 2: Number and proportion of wild and cultivated genotypes assigned to each a posteriori genetic group resulting from DAPC and proportion of major maternal lineages for each group.</caption><thead><tr><th rowspan="2">Genetic groups</th><th rowspan="2">Assigned genotypes* (%)</th><th colspan="4">Maternal lineage (%)</th></tr><tr><th>E1<sup>1</sup></th><th>E2<sup>2</sup></th><th>E3<sup>3</sup></th><th>L1.1</th></tr></thead><tbody><tr><td>Cultivated<sub>west</sub></td><td>50 (12.2)</td><td>50 (100.0)</td><td></td><td></td><td></td></tr><tr><td>Cultivated<sub>center</sub></td><td>138 (33.7)</td><td>129 (93.5)</td><td>3 (2.2)</td><td>6 (4.3)</td><td></td></tr><tr><td>Cultivated<sub>east</sub></td><td>59 (14.4)</td><td>59 (100.0)</td><td></td><td></td><td></td></tr><tr><td>Cultivated<sub>admixed</sub></td><td>150 (36.6)</td><td>129 (86.0)</td><td>14 (9.3)</td><td>6 (4.0)</td><td>1 (0.7)</td></tr><tr><td>Total (cultivated olive)</td><td>410</td><td>377 (91.95)</td><td>18 (4.40)</td><td>14 (3.40)</td><td>1 (0.24)</td></tr><tr><td>Wild<sub>west/center</sub></td><td>398 (55.1)</td><td>55 (13.8)</td><td>307 (77.1)</td><td>36 (9.0)</td><td></td></tr><tr><td>Wild<sub>east</sub></td><td>265 (36.7)</td><td>265 (100.0)</td><td></td><td></td><td></td></tr><tr><td>Wild<sub>admixed</sub></td><td>49 (6.8)</td><td>26 (53.1)</td><td>22 (44.9)</td><td>1 (2.0)</td><td></td></tr><tr><td>Total (wild olive)</td><td>722</td><td>346 (47.92)</td><td>338 (46.82)</td><td>38 (5.26)</td><td></td></tr><tr><td>Total</td><td>1132</td><td>723 (63.87)</td><td>356 (31.45)</td><td>52 (4.59)</td><td>1 (0.09)</td></tr></tbody></table>

E3 (total of 6 haplotypes) were detected only in 32 cultivars (7.8%; Table S2).

### 3.2. Genetic Clustering.
The genetic structure of Mediterranean olive was investigated using discriminant analysis of principal components (DAPC). The “find.clusters” function was used to determine the number of clusters maximizing the variation between clusters [30]. To avoid the loss of information, the function was performed with 200 principal components, accounting for more than 98% of the variance (Figure S1(a)). The Bayesian information criterion (BIC) was used to identify the optimal number of clusters, i.e., 11 clusters (Figure S1(b)). Based on these 11 clusters as a first analysis, DAPC clustering was represented according to the origin of olives classified in four *a priori* groups: western-central Mediterranean wild olive, eastern Mediterranean wild and cultivated olives, western cultivated olive, and central cultivated olive (Figure S1(c)). Western cultivated olive showed narrow genetic diversity (cluster 7), whereas those from the central Mediterranean Basin displayed high diversity included in 3 clusters (clusters 1, 3, and 6). Similarly, eastern and western-central wild olive displayed high diversity, i.e., 3 and 4 clusters, respectively (Figures S1(c) and S1(d)). Pairwise *F*<sub>ST</sub> values among the 11 predefined clusters resulting from DAPC ranged from 0.017 (cluster 5-cluster 8) to 0.129 (cluster 4-cluster 7) (Table S3).

Based on the assignation membership probability resulting from DAPC at *p* ≥ 0.8, five groups could be identified: (i) eastern Mediterranean wild (referred as Wild<sub>east</sub>), (ii) eastern cultivated olive (Cultivated<sub>east</sub>), (iii) western and central Mediterranean wild olive (Wild<sub>west-center</sub>), (iv) western Mediterranean cultivated olive (Cultivated<sub>west</sub>), and (v) central Mediterranean cultivated olive (Cultivated<sub>center</sub>). Although they are belonging to the same pool (Figures 2 and 3), Wild<sub>east</sub> and Cultivated<sub>east</sub> were considered as two distinct groups to describe the relationships between wild and cultivated olives in the eastern MB. Otherwise, at the *p* < 80% assignation level, wild and cultivated olive showed admixture: Wild<sub>admixed</sub> and Cultivated<sub>admixed</sub>, respectively (Figure 2; Tables 2 and S4). For admixed wild and cultivated forms, we noted the occurrence of the three plastid lineages with a high

proportion of E1 for cultivated olive (86.0%), whereas for admixed wild olive, close proportions of E1 and E2 were observed (53.1% and 44.9%, respectively, Table 2). When considering the *a priori* groups, more wild admixed genotypes were noted in the western-central part of MB. Similarly, more admixed cultivars were observed in the western-central MB area compared to the east (Table S5).

To investigate genetic relationships between the five groups defined above, principal coordinate analysis (PCoA) was performed (Figure S2). Most of the variation (13.06%) was explained by the first two axes. For both wild and cultivated olives, the first axis corresponded to the east-west spatial distribution at the MB scale, where both wild olive groups (Wild<sub>west-center</sub> and Wild<sub>east</sub>) were genetically distinct. The second axis separated Cultivated<sub>west</sub> and Cultivated<sub>center</sub> olives from Wild<sub>east</sub> and Cultivated<sub>east</sub> olives. The latter were clustered as one pool. Admixed genotypes for both Wild<sub>admixed</sub> and Cultivated<sub>admixed</sub> were plotted midway between the five genetic groups (Figure S2).

### 3.3. Genetic Variation and Relationships among Genetic Groups.
The mean heterozygosity (Ho = 0.701 observed) noted for wild olive was less than expected (He = 0.826) based on the Hardy-Weinberg equilibrium findings (Fis = 0.151; *p* < 0.001; Table 3), indicating a deficit of heterozygotes, as noted for wild genetic groups (Wild<sub>west-center</sub> and Wild<sub>east</sub>). These results may be explained by the subdivision of local populations into isolated and differentiated units (Wahlund effect), as revealed by the 11 predefined DAPC clusters (Figures S1(c) and S1(d)).

Allelic richness (Ar) was estimated and revealed a highly significant difference between wild olive and cultivars (24.58 vs 17.23; Mann-Whitney test, *p* < 0.001; Table 3). Otherwise, a highly significant difference was observed between Wild<sub>east</sub> and Cultivated<sub>east</sub> (12.95 vs 8.01), but not between Wild<sub>west-center</sub> and Wild<sub>east</sub> (10.67 and 12.95, respectively) or between Cultivated<sub>center</sub> and Cultivated<sub>east</sub>. However, within cultivated olive, Ar for Cultivated<sub>west</sub> was significantly lower than for Cultivated<sub>center</sub> and Cultivated<sub>east</sub> (4.19, 7.52, and 8.01, respectively). When focusing on the maternal lineage, more maternal lineages belonging to E2 and E3 were revealed in

<footnote>*Under the assignation probability *p* > 80%, ¹E1: E1.1, E1.2, E1.3, E1.4, E1.5, E1.6, E1.7, E1.8, E1.9, E1.10, E1.11, E1.12, E1.13, E1.14, E1.15, E1.20, E1.21, E1.27, ²E2: E2.1, E2.2, E2.3, E2.4, E2.5, E2.6, E2.7, E2.8, E2.12, E2.14, ³E3: E3.1, E3.2, E3.3, E3.5.</footnote>
<figure><img src="image_3.png" /><figcaption>FIGURE 2: Scatterplot from a discriminant analysis of principal components (DAPC) (a) using the four *a priori* groups. The sample origins are labeled within their 95% inertia ellipses, and individuals are connected to the corresponding group centroids. The number of PCs retained and discriminant functions are shown. Bar plot from DAPC highlights the membership of individuals in different groups according to geographical zones for both wild and cultivated olives (b).</figcaption></figure>

$\text{Cultivated}_{\text{center}}$ than in $\text{Cultivated}_{\text{east}}$ and $\text{Cultivated}_{\text{west}}$ olives (Table 2; Figure 3). Moreover, for cultivars identified as admixed on the basis of DAPC (membership assignation <0.8), E2 and E3 haplotypes were found in higher proportion than for other cultivated groups (Table S5).

The genetic differentiation ($F_{ST}$) values were significant between all pairs of the five groups (Wild$_{east}$ and Cultivated$_{east}$ treated separately, $p < 0.001$; Table 4). The mean $F_{ST}$ was 0.105 ($p < 0.0001$) for the five groups, indicating that 10% of the total genetic variation resulted from genetic differentiation between groups. The pairwise $F_{ST}$ ranged from 0.036 to 0.183 (Table 4). The highest values were observed between Wild$_{west-center}$ and all the other groups, whereas the lowest values were revealed between Wild$_{east}$ and Cultivated$_{east}$ groups. Relationships between the five

groups showed a clear distinction between three main groups: (i) Wild$_{west-center}$, (ii) Wild$_{east}$ and Cultivated$_{east}$, and (iii) Cultivated$_{west}$ and Cultivated$_{center}$, as supported by the high bootstrap values (Figure 3). The group Wild$_{west-center}$ was highly separated from the others based on both nuclear and plastid polymorphism, with the highest proportion of maternal lineages belonging to E2 and E3 haplotypes (86.1%; Table 2; Figure 3).

## 4. Discussion

Over the last two decades, substantial paleobotanical, archeological, historical, and molecular data have been accumulated on olive species and the history of its domestication [1, 4, 5, 7, 8, 11–13]. A global overview was
<figure><img src="image_4.png" /><figcaption>FIGURE 3: Unrooted neighbor-joining dendrogram based on the genetic differentiation index ($F_{ST}$), showing relationships among the five genetic groups as identified by DAPC and the proportion of the three maternal lineages (E1, E2, and E3) within each group.</figcaption></figure>

<table><caption>TABLE 3: Summary of genetic parameters for different groups.</caption><thead><tr><th>Groups</th><th>n</th><th>A</th><th>A%</th><th>Ho</th><th>He</th><th>Fis</th><th>Ar*</th></tr></thead><tbody><tr><td>Wild olive</td><td>722</td><td>420</td><td>98.3</td><td>0.701</td><td>0.826</td><td>0.151***</td><td>24.58a<sup>1</sup></td></tr><tr><td>Cultivated olive</td><td>410</td><td>276</td><td>64.6</td><td>0.765</td><td>0.731</td><td>-0.047***</td><td>17.23b<sup>1</sup></td></tr><tr><td>Wild<sub>west/center</sub></td><td>398</td><td>331</td><td>77.5</td><td>0.664</td><td>0.737</td><td>0.094***</td><td>10.67ac<sup>2</sup></td></tr><tr><td>Wild<sub>east</sub></td><td>265</td><td>362</td><td>84.7</td><td>0.737</td><td>0.806</td><td>0.084***</td><td>12.95a<sup>2</sup></td></tr><tr><td>Cultivated<sub>west</sub></td><td>50</td><td>79</td><td>18.5</td><td>0.770</td><td>0.589</td><td>-0.310***</td><td>4.19d<sup>2</sup></td></tr><tr><td>Cultivated<sub>center</sub></td><td>138</td><td>195</td><td>45.6</td><td>0.751</td><td>0.704</td><td>-0.067***</td><td>7.52b<sup>2</sup></td></tr><tr><td>Cultivated<sub>east</sub></td><td>59</td><td>160</td><td>37.4</td><td>0.752</td><td>0.712</td><td>-0.057***</td><td>8.01bc<sup>2</sup></td></tr></tbody></table>

n, size of the group; A, number of alleles; A%, percentage of alleles with respect to the total number of alleles observed within wild and cultivated samples (427 alleles); Ho, observed heterozygosity; He, unbiased expected heterozygosity; Fis, inbreeding coefficient; Ar, allelic richness; $^1$estimated based on standardized sample size of G = 410; $^2$estimated based on standardized sample size of G = 50. *Letters (a, b, c, and d) represent significance groupings after a Mann-Whitney pairwise comparison test ($p < 0.05$).

founded on the basis of several centers of primary selection across the MB [10–12]. However, the question remains unclear as to whether cultivated olives derived from a single primary domestication center followed by secondary diversification events or whether they are the result of independent primary selection events. A scenario of at least two independent primary selection centers in the eastern and central Mediterranean was proposed by Diez et al. [13]. Investigations on wild olive in the eastern Mediterranean were, however, limited to few sampled populations that were likely feral, as assumed by Diez et al. [13]. The lack of genuinely eastern wild populations has drastically limited the possibility of testing alternative complex domestication scenarios, as pointed out by Besnard and Robio de Casas [14] and Diez and Gaut [15]. Hence, investigating a comprehensive

sample of wild olives throughout the MB could bring insight to help solve questions related to primary domestication and secondary diversification centers. Moreover, contrary to the findings of Diez et al. [13], our eastern wild populations were clearly distinct from the western-central wild olive populations, thus indicating their genuine status.

Olive tree history is complex, as previously highlighted by several studies (see review in Besnard et al. [3]). Instead of multiple primary domestication centers, we argue in favor of a single primary domestication in the Levant, followed by human-mediated diffusion of the first domesticated forms and admixtures with wild olives in the central and western Mediterranean Basin. However, we cannot exclude the occurrence of minor domestication centers in western and central parts of the MB, as some varieties have been found to
<table><caption>TABLE 4: Pairwise genetic differentiation $F_{ST}$ values between the five groups as revealed by DAPC method for wild and cultivated olives.</caption><thead><tr><th></th><th>Cultivated<sub>center</sub></th><th>Cultivated<sub>west</sub></th><th>Cultivated<sub>east</sub></th><th>Wild<sub>east</sub></th></tr></thead><tbody><tr><td>Cultivated<sub>west</sub></td><td>0.075<sup>***</sup></td><td></td><td></td><td></td></tr><tr><td>Cultivated<sub>east</sub></td><td>0.052<sup>***</sup></td><td>0.140<sup>***</sup></td><td></td><td></td></tr><tr><td>Wild<sub>east</sub></td><td>0.060<sup>***</sup></td><td>0.142<sup>***</sup></td><td>0.036<sup>***</sup></td><td></td></tr><tr><td>Wild<sub>west-center</sub></td><td><strong>0.168<sup>***</sup></strong></td><td><strong>0.176<sup>***</sup></strong></td><td><strong>0.183<sup>***</sup></strong></td><td>0.144<sup>***</sup></td></tr></tbody></table>

Pairwise significance after sequential Bonferroni corrections: ***, significance at *p* < 0.001. Numbers in bold are the highest values observed and underlined are the lowest.

harbor maternal E2 or E3 lineages specific to local genetic resources, indicating their ancient local selection heritage (Figure 1, Table S4). Moreover, morphometric olive stone and charcoal analyses have revealed the use of wild olive before the Neolithic period, suggesting local domestication could have occurred in the western MB area [7]. Here, by investigating current varieties using both nuclear and plastid markers, we obtained evidence of primary selection and secondary diversification as two key processes in the history of olive domestication based on the following arguments. First, we used DAPC and identified a single group including both eastern wild olive (Wild<sub>east</sub>) and eastern cultivars (Cultivated<sub>east</sub>), thus indicating direct selection from wild olive populations, as suggested by Gurbuz-Veral et al. [37]. Second, most cultivated olives have an eastern-like maternal haplotype as a signature of the diffusion of the first domesticated olives from the eastern to western Mediterranean Basin [5]. Note that the above two arguments are supported by the genetic differentiation index ($F_{ST}$) between different groups, including wild and cultivated olive trees (Figure 3 and Table 4). Third, the allelic richness revealed highly significant differences between Wild<sub>east</sub> and Cultivated<sub>east</sub> olives. Contrary to other perennial fruit species such as apple [38], a substantial reduction in allelic diversity was observed between domesticated and wild olives across the MB (up to 30%), especially from the eastern MB (up to 38.1%; Table 3). This finding is in line with the selection pattern during the domestication process, as reviewed by Gaut et al. [39] and Besnard et al. [3] and references therein. Fourth, the genetic pattern of the Cultivated<sub>west</sub> and Cultivated<sub>center</sub> groups indicated a diversification process based on selection from crosses between the first domesticated olive forms and local olives, thus supporting the assumption of human-mediated diffusion of cultivars. Indeed, among varieties from the central MB, 35% (81 varieties) were admixed with limited gene flow from western-central wild populations (Figure 2; Table S4) and displaying the three maternal lineages (Table S5). Among varieties from the central MB, we found 11% (24) harboring maternal lineages E2 and E3. Moreover, as reported by Belaj et al. [40] and Klepo et al. [41], some central Mediterranean varieties retain wild-like phenotypic characteristics, such as low endocarp weight and a smooth endocarp surface. These findings suggest a second center of domestication, as reported by Diez et al. [13], but evidence to back this assumption has yet to be documented [3].

We argue here in favor of a diversification process occurring in the western and central MB. Fifth, Diez et al. [13] found that most first-degree relationships were from the same

genetic group (i.e., western cultivated olive; 96.3%) in which two cultivars from Spain (i.e., Gordal Sevillana and Lechin de Granada) had more than 60 first-degree relationships. Varieties harboring E2 and E3 from the western MB were found to be closely related within the western cultivated group such as Lechin de Sevilla with E2.3 maternal lineage (Table S2) and five first-degree relationships [13]. Regardless of the maternal lineages, the presence of highly related varieties indicated diversification based on crosses between cultivated olives as a key olive domestication process in the central and western MB.

## 5. Conclusion

Beyond a single primary olive tree domestication event, our investigation underlines the importance of admixtures within cultivated olive groups from the central and western Mediterranean Basin. Clarifying the evolutionary processes responsible for these groups will help gain important insight into accurately identify the genes under selection. This will also help to design methods for sampling of Mediterranean olive germplasm, including wild olives suitable for genome-wide association studies and genomic selection under the impact of climate change and within the sustainable oliviculture setting. Moreover, identifying olive diversity hotspots in the MB could also help to develop cost-effective diversity-prioritized approaches for *in situ* olive genetic resource conservation and management.

## Data Availability

The complete dataset is available upon request to the corresponding author: khadari@supagro.fr.

## Disclosure

First results from the present study have been presented as an oral communication within the Fourth International American Moroccan, Agriculture Science Conference “AMAS Conference IV” held in May 9–11, 2018, at Agropolis, Meknès, Morocco (see http://www.amas-conference.org/wp-content/uploads/2018/05/AMAS-IV-FinalProgramBooklet-5-3-18-RegularPrint.pdf).

## Conflicts of Interest

The authors declare that they have no conflicts of interest.

## Authors’ Contributions

Bouchaib Khadari and Ahmed El Bakkali contributed equally to this work.
## Acknowledgments

We thank A. Al Ibrahim, G. Besnard, M. Gurbuz, H. Haouane, A. Moukhli, and J. F. Terral for providing olive samples. This work was conducted at UMR AGAP. We thank L. Essalouh, S. Santoni, and Ch. Tollon for laboratory assistance. This research was supported by the project OliveMed/Agropolis Fondation no. 1202-066 through the *Investissements d'avenir/Labex Agro ANR-10-Labex-0001-01* managed by the French National Research Agency (ANR) and by the BeFOre project "Bioresources for Oliviculture" 2015-2019, H2020-MSCA-RISE-Marie Skłodowska-Curie Research and Innovation Staff Exchange, Grant Agreement no. 645595.

## Supplementary Materials

Table S1: list of the 35 wild olive populations. The sampling locations, the number of populations, individuals per location (size), and the GPS coordinates are given. Table S2: list of the 410 Mediterranean olive cultivars analyzed in the present study, along with their origins and maternal lineages. Table S3: pairwise genetic differentiation ($F_{ST}$) among the 11 subclusters, as identified by DAPC using the "find.clusters" function. Table S4: number of wild and cultivated olives per region of origin and the number and proportion of individuals assigned to each group based on the DAPC findings with a membership probability of 0.8. Table S5: proportion of maternal lineages according to the *a priori* grouping clusters for both wild and cultivated olives identified as admixed genotypes by DAPC with a membership assignation of *p* < 0.8. Figure S1: discriminant analysis of principal component (DAPC) results. Cumulative variance explained by the principal component analysis (PCA) relative to the number of principal components (PCs) retained in the analysis (a). Selection of the optimal number of clusters in the DAPC using the lowest Bayesian information criterion (BIC; b)). Comparison of clustering performed by DAPC (K = 11) and the *a priori* wild and cultivated olive groups (c). Squares represent the number of individuals in each pairwise comparison. Scatterplot from a DAPC of olive genotypes showing the relationships between the 11 identified clusters (d). Figure S2: principal coordinate analysis (PCoA) based on the simple matching coefficient showing the relationships among cultivated and wild olive genotypes according to groups resulting from DAPC. (*Supplementary Materials*)

## References

[1] D. Kaniewski, E. Van Campo, T. Boiy, J. F. Terral, B. Khadari, and G. Besnard, "Primary domestication and early uses of the emblematic olive tree: palaeobotanical, historical and molecular evidences from the middle east," *Biological Reviews*, vol. 87, no. 4, pp. 885-899, 2012.

[2] E. Palamarev, "Paleobotanical evidences of the tertiary history and origin of the Mediterranean sclerophylldendroflora," *Plant Systematics and Evolution*, vol. 162, no. 1-4, pp. 93-107, 1989.

[3] G. Besnard, J. F. Teral, and A. Cornille, "on the origins and domestication of the olive: a review and perspectives," *Annals of Botany*, vol. 121, no. 3, pp. 385-403, 2018.

[4] Y. Carrión, M. Ntinou, and E. Badal, "Olea europaea L. in the North Mediterranean Basin during the Pleniglacial and the Early-Middle Holocene," *Quaternary Science Reviews*, vol. 29, no. 7-8, pp. 952-968, 2010.

[5] G. Besnard, B. Khadari, M. Navascués et al., "The complex history of the olive tree: from late quaternary diversification of Mediterranean lineages to primary domestication in the northern Levant," *Proceedings of the Royal Society B: Biological Sciences*, vol. 280, no. 1756, article 20122833, 2013.

[6] H. Haouane, A. El Bakkali, A. Moukhli et al., "Genetic structure and core collection of the World Olive Germplasm Bank of Marrakech: towards the optimised management and use of Mediterranean olive genetic resources," *Genetica*, vol. 139, no. 9, pp. 1083-1094, 2011.

[7] J. F. Terral, N. Alonso, R. B. Capdevila et al., "Historical biogeography of olive domestication (Olea europaea L.) as revealed by geometrical morphometry applied to biological and archaeological material," *Journal of Biogeography*, vol. 31, no. 1, pp. 63-77, 2004.

[8] D. Zohary, M. Hopf, and E. Weiss, *Domestication of Plants in the Old World: The Origin and Spread of Cultivated Plants in Southwest Asia, Europe, and the Mediterranean Basin*, Oxford University Press, Oxford, UK, 2012.

[9] J. F. Terral, *La domestication de l'olivier (Olea europaea L.) en Méditerranée nord-occidentale: approche morphométrique et implications paléoclimatiques*, Ph.D. thesis, Université Montpellier II, Montpellier, France, 1997.

[10] G. Besnard, P. Baradat, C. Breton, B. Khadari, and A. Berville, "Olive domestication from structure of oleasters and cultivars using nuclear RAPDs and mitochondrial RFLPs," *Genetics Selection Evolution*, vol. 33, pp. S251-S268, 2001.

[11] L. Baldoni, N. Tosti, C. Ricciolini et al., "Genetic structure of wild and cultivated olives in the central Mediterranean basin," *Annals of Botany*, vol. 98, no. 5, pp. 935-942, 2006.

[12] C. Breton, M. Tersac, and A. Berville, "Genetic diversity and gene flow between the wild olive (oleaster, Oleaeuropaea L.) and the olive: several Plio-Pleistocene refuge zones in the Mediterranean basin suggested by simple sequence repeats analysis," *Journal of Biogeography*, vol. 33, no. 11, pp. 1916-1928, 2006.

[13] C. M. Diez, I. Trujillo, N. Martinez-Urdiroz et al., "Olive domestication and diversification in the Mediterranean Basin," *New Phytologist*, vol. 206, no. 1, pp. 436-447, 2015.

[14] G. Besnard and R. Rubio de Casas, "Single vs multiple independent olive domestications: the jury is (still) out," *New Phytologist*, vol. 209, no. 2, pp. 466-470, 2016.

[15] C. M. Diez and B. S. Gaut, "The jury may be out, but it is important that it deliberates: a response to Besnard and Rubio de Casas about olive domestication," *New Phytologist*, vol. 209, no. 2, pp. 471-473, 2016.

[16] A. El Bakkali, H. Haouane, A. Moukhli, E. Costes, P. Van Damme, and B. Khadari, "Construction of core collections suitable for association mapping to optimize use of Mediterranean olive (Olea europaea L.) genetic resources," *PLoS One*, vol. 8, no. 5, Article ID e61265, 2013.

[17] H. Hannachi, C. Breton, M. Msallem, S. Ben El Hadj, M. El Gazzah, and A. Bervillé, "Differences between native and introduced olive cultivars as revealed by morphology of drupes, oil composition and SSR Polymorphisms: a case study in Tunisia," *Scientia Horticultural*, vol. 116, no. 3, pp. 280-290, 2008.
[18] G. Besnard, A. El Bakkali, H. Haouane, D. Baali-Cherif, A. Moukhli, and B. Khadari, "Population genetics of Mediterranean and Saharan olives: geographic patterns of differentiation and evidence for early generations of admixture," *Annals of Botany*, vol. 112, no. 7, pp. 1293-1302, 2013.

[19] B. Khadari, J. Charafi, A. Moukhli, and M. Ater, "Substantial genetic diversity in cultivated Moroccan olive despite a single major cultivar: a paradoxical situation evidenced by the use of SSR loci," *Tree Genetic and Genomes*, vol. 4, no. 2, pp. 213-221, 2008.

[20] K. M. Sefc, M. S. Lopes, D. Mendonca et al., "Identification of microsatellite loci in olive (Olea europaea) and their characterization in Italian and Iberian olive trees," *Molecular Ecology*, vol. 9, no. 8, pp. 1171-1173, 2000.

[21] R. De la Rosa, C. M. James, and K. R. Tobutt, "Isolation and characterization of polymorphic microsatellites in olive (Olea europaea L.) and their transferability to other genera in the Oleaceae," *Molecular Ecology Notes*, vol. 2, no. 3, pp. 265-267, 2002.

[22] F. Carriero, G. Fontanazza, F. Cellini, and G. Giorio, "Identification of simple sequence repeats (SSRs) in olive (Olea europaea L.)," *Theoretical and Applied Genetics*, vol. 104, no. 2-3, pp. 301-307, 2002.

[23] G. Cipriani, M. T. Marrazzo, R. Marconi, A. Cimato, and R. Testolin, "Microsatellite markers isolated in olive (Olea europaea L.) are suitable for individual fingerprinting and reveal polymorphism within ancient cultivars," *Theoretical and Applied Genetics*, vol. 104, no. 2-3, pp. 223-228, 2002.

[24] G. Besnard, P. Hernandez, B. Khadari, G. Dorado, and V. Savolainen, "Genomic profiling of plastid DNA variation in the Mediterranean olive tree," *BMC Plant Biology*, vol. 11, p. 80, 2011.

[25] S. D. E. Park, *Trypanotolerance in West African cattle and the population genetic effects of selection*, Ph.D. thesis, University of Dublin, Dublin, Ireland, 2001.

[26] J. Goudet, "FSTAT (Version 1.2): a computer program to calculate F-statistics," *Journal of Heredity*, vol. 86, no. 6, pp. 485-486, 1995.

[27] R. J. Petit, A. El Mousadik, and O. Pons, "Identifying populations for conservation on the basis of genetic markers," *Conservation Biology*, vol. 12, no. 4, pp. 844-855, 1998.

[28] Z. A. Szpiech, M. Jakobsson, and N. A. Rosenberg, "ADZE: a rarefaction approach for counting alleles private to combinations of populations," *Bioinformatics*, vol. 24, no. 21, pp. 2498-2504, 2008.

[29] T. Jombart, S. Devillard, and F. Balloux, "Discriminant analysis of principal components: a new method for the analysis of genetically structured populations," *BMC Genetics*, vol. 11, no. 1, p. 94, 2010.

[30] T. Jombart, "adegenet: a R package for the multivariate analysis of genetic markers," *Bioinformatics*, vol. 24, no. 11, pp. 1403-1405, 2008.

[31] J. K. Pritchard, M. Stephens, and P. Donnelly, "Inference of population structure from multilocus genotype data," *Genetics*, vol. 155, no. 2, pp. 945-959, 2000.

[32] R. R. Sokal and C. D. Michener, *A Statistical Method for Evaluating Systematic Relationships*, vol. 38 of University of Kansas Science Bulletin, University of Kansas, Lawrence, KS, USA, 1958.

[33] X. Perrier and J. P. Jacquemoud-Collet, "DARwin software," 2006, http://darwin.cirad.fr/darwin.

[34] B. S. Weir and C. C. Cockerham, "Estimating F-statistics for the analysis of population structure," *Evolution*, vol. 38, no. 6, pp. 1358-1370, 1984.

[35] M. Raymond and F. Rousset, "GENEPOP (version 1.2): population genetics software for exact tests and ecumenism," *Journal of Heredity*, vol. 86, no. 3, pp. 248-249, 1995.

[36] N. Takezaki, M. Nei, and K. Tamura, "POPTREE2: software for constructing population trees from allele frequency data and computing other population statistics with Windows-interface," *Molecular Biology and Evolution*, vol. 27, no. 4, pp. 747-752, 2010.

[37] M. Gurbuz-Veral, A. El Bakkali, L. Essalouh et al., "New insights on olive domestication in Turkey," *Acta Horticulturae*, vol. 1199, pp. 15-20, 2018.

[38] A. Cornille, P. Gladieux, M. J. M. Smulders et al., "New insight into the history of domesticated apple: secondary contribution of the European wild apple to the genome of cultivated varieties," *PLoS Genetics*, vol. 8, no. 5, article e1002703, 2012.

[39] B. S. Gaut, C. M. Díez, and P. L. Morrell, "Genomics and the contrasting dynamics of annual and perennial domestication," *Trends in Genetics*, vol. 31, no. 12, pp. 709-719, 2015.

[40] A. Belaj, L. Leon, Z. Satovic, and R. de la Rosa, "Variability of wild olives (Olea europaea subsp. europaea var. sylvestris) analyzed by agro-morphological traits and SSR markers," *Scientia Horticulturae*, vol. 129, no. 4, pp. 561-569, 2011.

[41] T. Klepo, R. De la Rosa, Z. Satovic, L. Leon, and A. Belaj, "Utility of wild germplasm in olive breeding," *Scientia Horticulturae*, vol. 152, pp. 92-101, 2013.
<figure><img src="image_5.png" /></figure>

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

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

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

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

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

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

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

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

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

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

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

Hindawi

Submit your manuscripts at
www.hindawi.com

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

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

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

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

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

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

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

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