Abstract

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 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 [913], but it remains unclear whether this reflects secondary diversification or multiple independent primary domestication events [1315]. 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 agro-ecosystems 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 [2023] (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 FSTAT 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 (FST, [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 FST 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 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 FST 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 , five groups could be identified: (i) eastern Mediterranean wild (referred as Wildeast), (ii) eastern cultivated olive (Cultivatedeast), (iii) western and central Mediterranean wild olive (Wildwest-center), (iv) western Mediterranean cultivated olive (Cultivatedwest), and (v) central Mediterranean cultivated olive (Cultivatedcenter). Although they are belonging to the same pool (Figures 2 and 3), Wildeast and Cultivatedeast were considered as two distinct groups to describe the relationships between wild and cultivated olives in the eastern MB. Otherwise, at the assignation level, wild and cultivated olive showed admixture: Wildadmixed and Cultivatedadmixed, 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 (Wildwest-center and Wildeast) were genetically distinct. The second axis separated Cultivatedwest and Cultivatedcenter olives from Wildeast and Cultivatedeast olives. The latter were clustered as one pool. Admixed genotypes for both Wildadmixed and Cultivatedadmixed 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; ; Table 3), indicating a deficit of heterozygotes, as noted for wild genetic groups (Wildwest-center and Wildeast). 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, ; Table 3). Otherwise, a highly significant difference was observed between Wildeast and Cultivatedeast (12.95 vs 8.01), but not between Wildwest-center and Wildeast (10.67 and 12.95, respectively) or between Cultivatedcenter and Cultivatedeast. However, within cultivated olive, Ar for Cultivatedwest was significantly lower than for Cultivatedcenter and Cultivatedeast (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 Cultivatedcenter than in Cultivatedeast and Cultivatedwest 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 (FST) values were significant between all pairs of the five groups (Wildeast and Cultivatedeast treated separately, ; Table 4). The mean FST was 0.105 () for the five groups, indicating that 10% of the total genetic variation resulted from genetic differentiation between groups. The pairwise FST ranged from 0.036 to 0.183 (Table 4). The highest values were observed between Wildwest-center and all the other groups, whereas the lowest values were revealed between Wildeast and Cultivatedeast groups. Relationships between the five groups showed a clear distinction between three main groups: (i) Wildwest-center, (ii) Wildeast and Cultivatedeast, and (iii) Cultivatedwest and Cultivatedcenter, as supported by the high bootstrap values (Figure 3). The group Wildwest-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, archaeological, historical, and molecular data have been accumulated on olive species and the history of its domestication [1, 4, 5, 7, 8, 1113]. A global overview was founded on the basis of several centers of primary selection across the MB [1012]. 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 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 (Wildeast) and eastern cultivars (Cultivatedeast), 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 (FST) between different groups, including wild and cultivated olive trees (Figure 3 and Table 4). Third, the allelic richness revealed highly significant differences between Wildeast and Cultivatedeast 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 Cultivatedwest and Cultivatedcenter 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: [email protected].

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 (FST) 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 . 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)