# Diversity of Algerian oases date palm (Phoenix dactylifera L., Arecaceae) : Heterozygote excess and cryptic structure suggest farmer management had a major impact on diversity

Souhila Moussouni, Jean-Christophe Pintaud, Yves Vigouroux, Nadia Bouguedoura

► To cite this version:

Souhila Moussouni, Jean-Christophe Pintaud, Yves Vigouroux, Nadia Bouguedoura. Diversity of Algerian oases date palm (Phoenix dactylifera L., Arecaceae): Heterozygote excess and cryptic structure suggest farmer management had a major impact on diversity. PLoS ONE, 2017, 12 (4), pp.e0175232. 10.1371/journal.pone.0175232. [https://doi.org/10.1371/journal.pone.0175232](https://doi.org/10.1371/journal.pone.0175232). [hal-03346207](https:// HAL Id: hal-03346207)

Submitted on 16 Sep 2021

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Diversity of Algerian oases date palm (*Phoenix dactylifera* L., *Arecaceae*): Heterozygote excess and cryptic structure suggest farmer management had a major impact on diversity

Souhila Moussouni¹*, Jean-Christophe Pintaud², Yves Vigouroux², Nadia Bouguedoura¹

1 Université des Science et de la Technologie Houari Boumediene (USTHB), Faculté des Sciences Biologiques, Laboratoire de Recherche sur les Zones Arides (LRZA), Alger, Algeria, 2 Institut de Recherche pour le Développement, UMR DIADE, DYNADIV, Montpellier, France

* moussouni.souhila@gmail.com

## Abstract

Date palm (*Phoenix dactylifera* L.) is the mainstay of oasis agriculture in the Saharan region. It is cultivated in a large part of the Mediterranean coastal area of the Sahara and in most isolated oases in the Algerian desert. We sampled 10 oases in Algeria to understand the structure of date palm diversity from the coastal area to a very isolated desert location. We used 18 microsatellite markers and a chloroplast minisatellite to characterize 414 individual palm trees corresponding to 114 named varieties. We found a significant negative inbreeding coefficient, suggesting active farmer selection for heterozygous individuals. Three distinct genetic clusters were identified, a ubiquitous set of varieties found across the different oases, and two clusters, one of which was specific to the northern area, and the other to the drier southern area of the Algerian Sahara. The ubiquitous cluster presented very striking chloroplast diversity, signing the frequency of haplotypes found in Saudi Arabia, the most eastern part of the date palm range. Exchanges of Middle Eastern and Algerian date palms are known to have occurred and could have led to the introduction of this particular chlorotype. However, Algerian nuclear diversity was not of eastern origin. Our study strongly suggests that the peculiar chloroplastic diversity of date palm is maintained by farmers and could originate from date palms introduced from the Middle East a long time ago, which since then, has been strongly introgressed. This study illustrates the complex structure of date palm diversity in Algerian oases and the role of farmers in shaping such cryptic diversity.

## Introduction

The date palm (*Phoenix dactylifera* L., *Arecaceae*) is a perennial monocotyledon (2n = 36). It is an ecologically, culturally and economically important crop, widely cultivated in arid and semi-arid Mediterranean regions, in the Sahara, and in the Middle East [1, 2, 3]. More than

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

### OPEN ACCESS

**Citation:** Moussouni S, Pintaud J-C, Vigouroux Y, Bouguedoura N (2017) Diversity of Algerian oases date palm (*Phoenix dactylifera* L., *Arecaceae*): Heterozygote excess and cryptic structure suggest farmer management had a major impact on diversity. PLoS ONE 12(4): e0175232. https://doi.org/10.1371/journal.pone.0175232

**Editor:** Tzen-Yuh Chiang, National Cheng Kung University, TAIWAN

**Received:** August 8, 2016

**Accepted:** March 22, 2017

**Published:** April 14, 2017

**Copyright:** © 2017 Moussouni et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was financed by the AUF MeRSi project (6313PS001) and the Algerian Ministère de l'Enseignement Supérieur et de la Recherche Scientifique. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
**Competing interests:** The authors have declared that no competing interests exist.

3,000 cultivars are estimated to be used for date production worldwide, of which around 60 are widely grown and have important national and international markets [4].

Early cultivation of the date palm is recorded in the eastern part of its cultivation area, in southern Mesopotamia in the 5th Millennium BC. In the western part of its cultivation area, evidence for domesticated date palm has been found in Egypt in the 4th century BC [5]. Genetic studies of date palm clearly separate the eastern and western group [6, 7, 8, 9, 10]. The two primary gene pools (eastern and western) are also observed in the maternally-inherited chloroplast genome [7] and the paternally-inherited Y chromosome [11]. Two major alleles were found in the chloroplast, an occidental and oriental haplotype [7].

In Algeria, chloroplast diversity includes roughly 70% of the eastern chloroplast [10]. In neighboring Egypt, Tunisia and Morocco, the proportion of eastern haplotype only ranges from 11% to 42%, but Algerian nuclear diversity is similar to that found in its neighboring countries, Tunisia and Morocco [10]. The contrast between nuclear and chloroplast diversity observed in Algeria thus remains largely unexplained, perhaps because of the only fragmentary analysis of the diversity of the 1,000 varieties described in Algeria.

Date palms are cultivated in Algerian oases in most of the regions south of the Saharan Atlas Mountains. In 2002, Algerian date palm groves contained 13.5 million trees occupying 120,830 ha, while in 2015, 18 million date palms occupied 169,380 ha [12]. Nearly 1,000 cultivars clonally propagated from offshoots have been inventoried and their distribution shows a very marked breakdown into eastern, central and western parts of the country. Some cultivars are found in two or three regions but most are restricted to their area of origin [13]. What is more, cultivars are not evenly distributed across oases, as they are adapted to slightly different types of soil, ranges of temperature and humidity, and often do not produce a satisfactory yield when cultivated outside their place of origin [14, 12]. Palms grown from seeds occur randomly in the oases and are called "khalts" or "dgouls". Khalts represent up to 10% of a population of date palm sand area valuable resource for new selection by farmers.

The aims of the present study were to (1) genetically characterize date palm agrobiodiversity in Algeria using a set of representative female fruiting date palm cultivars and (2) investigate the structure of genetic diversity in the oases south of the Saharan Atlas based on nuclear and chloroplast genotyping.

# Materials and methods

## Ethics statement

No specific permissions are required for the activities conducted in this study. The study did not involve endangered or protected species.

## Plant material, DNA extraction and quantification

A total of 414 samples of date palm (*Phoenix dactylifera* L.) were sampled from 10 oases distributed in the north-western (3 oases), southern (1 oases), north-central (2 oases) and north-eastern (4 oases) date palm regions of Algeria (Fig 1). The plant material consisted in portions of young leaflets collected from 12–15 year-old palms whenever possible. Fresh material was dried in silica gel and stored at room temperature. The 414 samples corresponded to 114 named female varieties, 14 unidentified female varieties and 10 male plants. DNA was extracted from 45 mg of dried leaf tissue. The samples were ground to fine powder in a Tissue-Lyser II homogenizer (Qiagen), in the presence of steel balls. DNA was extracted using the DNeasy Plant MiniKit (Qiagen) following the manufacturer's instructions. DNA quality was checked by electrophoresis on 1% agarose gel [15]. The DNA concentration was determined using a NanoDrop spectro-photometer.
<figure><img src="image_4.png" /><figcaption>**Fig 1. Map showing the distribution of the 10 oases sampled in Algeria.** The 414 samples came from the following sources: Biskra 13 cultivars; Oued Souf 22 cultivars; Touggourt 30 cultivars; Ouargla 5 cultivars; Ghardaia 17 cultivars; El Menia 10 cultivars; Timimoun 16 cultivars; Adrar 18 cultivars; Beni Abbes 14 cultivars, and Tamanrasset 4 cultivars.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0175232.g001

## SSR genotyping

One dodecanucleotide plastid minisatellite located in the *trnG-trnfM* inter gene space was genotyped [7]. The occidental chlorotype is characterized by three repetitions of the 12 bp motif (type 3 chlorotype), while the oriental chlorotype has four repetitions (type 4 chlorotype) [16, 17]. This polymorphism provides the easiest access to the differentiation of the occidental and oriental chloroplast genomes, which is otherwise extensive. Indeed, sequence alignment of ca. 50 kb of an occidental chloroplast from Elche, Spain and an oriental one (cv Khalass) from the Saudi Arabian peninsula, revealed 18 SNPs along with a number of other mutations, including microsatellites, indels and small inversions [18]. The comparison of two complete sequences of oriental chloroplast genomes of cv Khalass and Aseel [19] resulted in only 3 SNPs. The two amplified alleles at this locus are 242 bp and 254 bp in length.

The following reaction mixture was used for genotyping with *trnG-trnfM*-primers: 1 µl of DNA, PCR buffer, 5 µl of FailSaif premix, 0.1 µl of Taq DNA polymerase, 0.1 µl forward primer, 0.1 µl of antisense primer, and water qs 3.6 µl Promega to reach a final volume of 10 µl. The amplification protocol comprised a denaturation step at 95°C for 3 minutes followed by 30 cycles each comprising a denaturation step at 94°C for 30 seconds, an annealing step at 56°C for 1 minute 30 seconds, and an elongation step at 72°C for 1 minute 30 seconds. Genotyping was performed with the Qiaaxel DNA analyzer.

Eighteen nuclear microsatellite loci were also used for this study (Table 1). These loci were analyzed on 192 samples including 178 female clones representing 114 identified varieties, 11 unidentified female clones, one seed-grown female genotype and 2 male plants. The individuals sampled were collected from the following locations: Biskra (13 cultivars / 20 individuals), Oued Souf (22 cultivars / 30 individuals), Touggourt (30 cultivars / 32 individuals), Ouargla (5 cultivars / 5 individuals), Ghardaia (17 cultivars / 29 individuals), El Menia (10 cultivars / 17 individuals), Timimoun (16 cultivars / 17 individuals), Adrar (18 cultivars / 20 individuals), Beni Abbes (14 cultivars / 15 individuals) and Tamanrasset (4 cultivars / 7 individuals) (S1 Table).

For all DNA extracts, the concentration was homogenized at 5 ng/µl. Amplification reactions were performed in a final volume of 20 µl containing 15 ng of template DNA, 10×
<table><caption>Table 1. Summary data for 18 microsatellite loci.</caption><thead><tr><th>Marker name</th><th>Short name</th><th>Forward primer sequence (5'→3')</th><th>Reverse primer sequence (5'→3')</th><th>SSR motifs</th><th>Length (pd)</th><th>At (C')</th><th>Reference</th><th>GeneBank accessions</th></tr></thead><tbody><tr><td>mPdCIR010</td><td>Pd10</td><td>ACCCCGGACGTGAGGTG</td><td>CGTCGATCTCCTCCTTIGT</td><td>(GA)<sub>22</sub></td><td>164</td><td>56</td><td>[20]</td><td>AJ571673</td></tr><tr><td>mPdCIR015</td><td>Pd15</td><td>AGCTGCTCCTCCCTTCTTA</td><td>GCTCGGTTGGACTTTGTCT</td><td>(GA)<sub>15</sub></td><td>135</td><td>52</td><td>[20]</td><td>AJ571674</td></tr><tr><td>mPdCIR016</td><td>Pd16</td><td>AGCGGAATGAAAGGTAT</td><td>ATGAAACGTGCCAAATGTC</td><td>(GA)<sub>14</sub></td><td>134</td><td>52</td><td>[20]</td><td>AJ571675</td></tr><tr><td>mPdCIR025</td><td>Pd25</td><td>GCACGAGAAGGCTTATAGT</td><td>CCCCTCATTAGGATTCCTAC</td><td>(GA)<sub>22</sub></td><td>230</td><td>49</td><td>[20]</td><td>AJ571676</td></tr><tr><td>mPdCIR032</td><td>Pd32</td><td>CAAATCTTTGCCSTGAG</td><td>GGTGTGGAGTAATCATGTAGTAG</td><td>(GA)<sub>19</sub></td><td>300</td><td>52</td><td>[20]</td><td>AJ571677</td></tr><tr><td>mPdCIR035</td><td>Pd35</td><td>ACAAACGGCATGGGATTAC</td><td>CCGCAGCTCACCCTTCTTAT</td><td>(GA)<sub>15</sub></td><td>188</td><td>54</td><td>[20]</td><td>AJ571678</td></tr><tr><td>mPdCIR057</td><td>Pd57</td><td>AAGCAGCAGCCCTTCCGTAG</td><td>GTTCCTCACTCGCCCAAAATAC</td><td>(GA)<sub>20</sub></td><td>282</td><td>55</td><td>[20]</td><td>AJ571682</td></tr><tr><td>mPdCIR063</td><td>Pd63</td><td>CTTTTATGTGGTCTGAGAGA</td><td>TCTCTGATCTTTGGGTTCTGT</td><td>(GA)<sub>17</sub></td><td>167</td><td>50</td><td>[20]</td><td>AJ571683</td></tr><tr><td>mPdCIR078</td><td>Pd78</td><td>TGGATTTCCATTGTGAG</td><td>CCCGAAGAGACGCTATT</td><td>(GA)<sub>13</sub></td><td>121</td><td>50</td><td>[20]</td><td>AJ571685</td></tr><tr><td>mPdCIR085</td><td>Pd85</td><td>GAGAGGGTGGTGTATT</td><td>TTCATCCAGAACCACAGTA</td><td>(GA)<sub>29</sub></td><td>179</td><td>50</td><td>[20]</td><td>AJ571686</td></tr><tr><td>PdAG1-ssr</td><td>AG1</td><td>TCTGATTTCGTTTACTTCTTAGGA</td><td>TTCATATTCAAGTTGTCGGGTGA</td><td>(GA)</td><td>260</td><td>52</td><td>[21]</td><td></td></tr><tr><td>PdCUC3-ssr1</td><td>CUC3-1</td><td>CGTGGACTCATGACTCGCATGTCC</td><td>GGTCCTTGCCGGTGGCCTTC</td><td>(GT)<sub>14</sub></td><td>330</td><td>60</td><td>[22]</td><td>HM622273</td></tr><tr><td>PdCUC3-ssr2</td><td>CUC3-2</td><td>ACATTGCTCTTTTGCCATGGGCT</td><td>CCAGCAGTGGGGTTCCGGT</td><td>(GA)<sub>22</sub></td><td>350</td><td>59</td><td>[22]</td><td>HM622273</td></tr><tr><td>PdAP3-ssr</td><td>AP3</td><td>GAGAAATAGAGCTGTGCAAG</td><td>CTGCAGTACTCGGAGAACTTG</td><td>(GA)<sub>25</sub></td><td>331</td><td>57</td><td>[23]</td><td>KC188337</td></tr><tr><td>mPdIRD013</td><td>P13</td><td>GGCGAGACAGGAGATCGTAA</td><td>CTTGACTGCTTCTGCTGCTG</td><td>(CAC)<sub>6</sub></td><td>204</td><td>60</td><td>[24]</td><td>PDK_20s1496731g002</td></tr><tr><td>mPdIRD031</td><td>P31</td><td>GCAGGTGGACTGCAAAATCT</td><td>CTATTGGGGTGCTGATCCAT</td><td>(CCA)<sub>7</sub></td><td>198</td><td>60</td><td>[24]</td><td>PDK_20s1419261g003</td></tr><tr><td>mPdIRD033</td><td>P33</td><td>GGAGCATACAGTGGGTTTGC</td><td>CAGCCTGGGAATGAGGATAG</td><td>(CAG)<sub>7</sub></td><td>199</td><td>60</td><td>[24]</td><td>PDK_20s1569281g001</td></tr><tr><td>mPdIRD040</td><td>P40</td><td>GAGAGATGCGTCAGGGAATC</td><td>CCAGAATCTTCCAAGCAAGC</td><td>(CCAGT)<sub>4</sub></td><td>193</td><td>60</td><td>60</td><td>PDK_20s1327401g002</td></tr></tbody></table>

Information concerning the microsatellite markers used to genotype date palm cultivars (identified from the genome of *Phoenix dactylifera*) using the primer sequences. The SSR motifs, PCR product size, primer annealing temperatures, referenced and GeneBank accession numbers are given. PCR = polymerase chain reaction; SSR = simple sequence repeat; At = annealing temperature.

https://doi.org/10.1371/journal.pone.0175232.t001
reaction buffer, 5 pmol each of forward and reverse primer, 0.2 mM of each deoxynucleotide, 2 mM MgCl₂, and 1 unit of Taq polymerase (Sigma). The forward primers were 5' labelled with one of three fluorescent compounds (6-FAM, NED or HEX) to enable analysis with automated sequencers. PCR was carried out using an Eppendorf Mastercycler pro equipped with vapor protect technology (AG, Hamburg, Germany). The PCR conditions were: an initial denaturation step at 95°C for 5 min, followed by 35 cycles each consisting of denaturation at 95°C for 30 seconds, hybridization at 51-57°C for 60 seconds and elongation at 72°C for 30 seconds, and a final step at 60°C for 30 min. Amplified products were run on an ABI 3130XL Genetic Analyzer (Applied Biosystems, USA) and alleles were scored using the GeneMapper V3.7 software (Applied Biosystems).

## Diversity and statistical analysis

Genetic diversity parameters, i.e. the number of alleles per locus (Na), allelic frequencies, heterozygosity and polymorphism information content value (PIC value) were estimated with Power Marker version 3.25 software [25].

$F_{IS}$ and $F_{ST}$ indices were calculated according to Weir and Cockram [26] using Fstat [27] and Genetix 4.05 software [28]. Allelic richness was computed with Fstat [27].

Anon-parametric Wilcoxon test was used to determine the significance of differences in diversity between groups using R [29]. Probability values (P-value) of the test are given for each comparison between pairs of populations.

Cluster analysis was conducted by generating a genetic distance matrix by calculating the shared allele frequencies using the neighbor-joining algorithm implemented in Mega version 5.05 software [30]. We also assessed clustering based on the Bayesian approach implemented in STRUCTURE version 2.3 [31]. We assessed different numbers of groups (K) ranging from 1 to 10, with 10 repetitions for each given value of K. For each run, we used a burn in period of 10,000 iterations, and a post burn in simulation length of 1,000,000. The most probable number of clusters was assessed using both the likelihood and an *ad hoc* quantity based on the second order rate of change in the log probability of data between different K values [32]. Diversity (number of alleles, heterozygosity) was assessed in the different clusters and compared using the Wilcoxon test. We assigned an individual to a group if its ancestry in this group was > 75% [33]. To check for a difference in chlorotype frequency between groups, we performed a Chi-square test, using SPSS Statistics 20.

## Results

### Chloroplast minisatellite and nuclear SSR genotyping

The two genotyping methods used (either Qiaxel or the ABI genetic analyzer) were consistent and revealed two chlorotypes, previously reported as occidental and oriental [7] (S2 Table). The occidental (western) chlorotype was found in 31.3% and the oriental (eastern) in 68.7% of the total sample. The proportion varied among oases (Fig 2C).

Linkage disequilibrium was very low between SSR markers (S3 Table). A total of 143 alleles were identified with the 18 nuclear markers. The average number of alleles was 7.9, average expected heterozygosity was 0.60, and the average PIC was 0.57 (Table 2). The average number of alleles per locus at the individual oases ranged from 3.83 (Ouargla) to 5.78 (Oued Souf). Allelic richness varied slightly between oases with values ranging from 3.36 (Touggourt) to 3.85 (Beni Abbes). The same was true of expected heterozygosity, which ranged from 0.53 (Ouargla) to 0.59 (Oued Souf and Beni Abbes), and of observed heterozygosity, which ranged from 0.56 (Tamanrasset) to 0.65 (Ghardaia).
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Bayesian cluster analysis using the STRUCTURE program: Results for K = 3. (A) Plots of (a) maximum log likelihood over the 10 runs and (b) delta k from the structure analysis was calculated according to the method of Evanno et al. (B) Estimated population structure inferred from all the individual date palms for K = 3. Each individual is represented by a thin vertical line divided into K colored segments representing the fraction of the individual's estimated membership of the K clusters. Pie charts show the frequencies of the haplotypes belonging to the three structure groups. (C) Geographic distributions of the 10 oases and two haplotypes. The pie chart shows the proportions of haplotypes in each oasis. (D) Sampling location of the date palms. Pie charts show the proportion of membership of each sampled population inferred by structure analysis (K = 3).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0175232.g002

No significant difference in diversity (allelic richness, heterozygosity observed or expected) between oases was observed using a Bonferroni corrected p-value for multiple testing (S4, S5, S6 Tables). Significant negative fixation indices (FIS) were obtained for two out of 10 oases (Table 3). Only three isolated southern oases: Tamanrasset, Adrar and Beni Abbes, had a deficit of heterozygous individuals (Table 3).

## Genetic differentiation between populations

Overall differentiation was weak between oases, with an overall average differentiation of 0.0119. Overall, there was isolation by geographical distance (Mantel test, p-value = 0.039). Geographically close populations were generally weakly differentiated while greater differentiation was observed for populations sampled further apart (Table 4).

## Variety name and diversity

Based on an individual phylogenetic tree, relationships between varieties with similar names are easy to see. Varietal names need to be compared because of the many synonyms, such as: Al Kayed and Tanteboucht; Tgaza and Takerboucht; Halimi and Masri; Aharthan and Harthan Oumazer; Chikh, Mhammed Chikh and Hamuri; Bent Cherk and Cherka; Tati and
<table><caption>Table 2. Descriptive genetic parameters for 18 microsatellite loci analyzed on 192 individual date palms.</caption><thead><tr><th>Marker</th><th>N° of genotypes</th><th>N° of alleles</th><th>Major allele frequency</th><th>Gene diversity</th><th>Expected heterozygosity He</th><th>PIC</th></tr></thead><tbody><tr><td>mPdCIR010</td><td>31</td><td>10</td><td>0.28</td><td>0.80</td><td>0.86</td><td>0.78</td></tr><tr><td>mPdCIR015</td><td>19</td><td>8</td><td>0.50</td><td>0.70</td><td>0.68</td><td>0.67</td></tr><tr><td>mPdCIR016</td><td>8</td><td>4</td><td>0.51</td><td>0.61</td><td>0.56</td><td>0.53</td></tr><tr><td>mPdCIR025</td><td>19</td><td>10</td><td>0.29</td><td>0.77</td><td>0.87</td><td>0.74</td></tr><tr><td>mPdCIR032</td><td>25</td><td>9</td><td>0.41</td><td>0.75</td><td>0.77</td><td>0.72</td></tr><tr><td>mPdCIR035</td><td>15</td><td>7</td><td>0.70</td><td>0.48</td><td>0.27</td><td>0.46</td></tr><tr><td>mPdCIR057</td><td>18</td><td>7</td><td>0.40</td><td>0.72</td><td>0.77</td><td>0.68</td></tr><tr><td>mPdCIR063</td><td>14</td><td>6</td><td>0.34</td><td>0.75</td><td>0.74</td><td>0.71</td></tr><tr><td>mPdCIR078</td><td>41</td><td>16</td><td>0.31</td><td>0.84</td><td>0.83</td><td>0.82</td></tr><tr><td>mPdCIR085</td><td>30</td><td>10</td><td>0.25</td><td>0.83</td><td>0.84</td><td>0.80</td></tr><tr><td>PdAG1-ssr</td><td>54</td><td>19</td><td>0.27</td><td>0.87</td><td>0.86</td><td>0.86</td></tr><tr><td>PdCUC3-ssr1</td><td>1</td><td>1</td><td>1.00</td><td>0.00</td><td>0.00</td><td>0.00</td></tr><tr><td>PdCUC3-ssr2</td><td>42</td><td>12</td><td>0.17</td><td>0.88</td><td>0.85</td><td>0.87</td></tr><tr><td>PdAP3-ssr</td><td>20</td><td>9</td><td>0.37</td><td>0.74</td><td>0.83</td><td>0.70</td></tr><tr><td>mPdIRD013</td><td>2</td><td>2</td><td>0.98</td><td>0.05</td><td>0.05</td><td>0.05</td></tr><tr><td>mPdIRD031</td><td>5</td><td>3</td><td>0.83</td><td>0.29</td><td>0.31</td><td>0.26</td></tr><tr><td>mPdIRD033</td><td>5</td><td>4</td><td>0.91</td><td>0.16</td><td>0.16</td><td>0.16</td></tr><tr><td>mPdIRD040</td><td>9</td><td>5</td><td>0.70</td><td>0.47</td><td>0.49</td><td>0.43</td></tr><tr><td>Mean</td><td>19.89</td><td>7.94</td><td>0.51</td><td>0.59</td><td>0.60</td><td>0.57</td></tr></tbody></table>

The number of genotypes, the number of alleles, the frequency of the main alleles, gene diversity, heterozygosity and the polymorphism information content (PIC) of 18 loci in *Phoenix dactylifera* are given. PIC = polymorphism information content.

https://doi.org/10.1371/journal.pone.0175232.t002

Tacherwint; Kesba and Sokrya; Azizaou and Adam Zrak; Alig, Bu'Rus and Our'Rous. In a few cases, we also observed varieties that shared the same name but were relatively distant. For example, three samples of the variety Hamraya collected at three different oases (Biskra, Tougourt and Oued Souf) differed genetically.

<table><caption>Table 3. Summary statistics for 18 microsatellite loci in 10 oases.</caption><thead><tr><th>Oasis</th><th>N°</th><th>Na</th><th>Ar</th><th>Ho</th><th>He</th><th>F<sub>IS</sub></th><th>P value F<sub>IS</sub></th></tr></thead><tbody><tr><td>Biskra</td><td>18</td><td>4.72</td><td>3.45</td><td>0.58</td><td>0.56</td><td>-0.02</td><td>0.3216</td></tr><tr><td>Touggourt</td><td>32</td><td>5.33</td><td>3.36</td><td>0.59</td><td>0.54</td><td>-0.06</td><td>0.0097</td></tr><tr><td>Oued Souf</td><td>30</td><td>5.78</td><td>3.73</td><td>0.60</td><td>0.59</td><td>-0.01</td><td>0.2457</td></tr><tr><td>Ouargla</td><td>5</td><td>3.83</td><td>3.83</td><td>0.61</td><td>0.53</td><td>-0.05</td><td>0.2563</td></tr><tr><td>Ghardaia</td><td>29</td><td>5.22</td><td>3.53</td><td>0.65</td><td>0.56</td><td>-0.13</td><td>0.0000</td></tr><tr><td>Tamanrasset</td><td>7</td><td>4.11</td><td>3.68</td><td>0.56</td><td>0.57</td><td>0.09</td><td>0.9603</td></tr><tr><td>El Menia</td><td>17</td><td>4.78</td><td>3.47</td><td>0.58</td><td>0.55</td><td>-0.03</td><td>0.1719</td></tr><tr><td>Timimoun</td><td>16</td><td>5.11</td><td>3.67</td><td>0.64</td><td>0.57</td><td>-0.04</td><td>0.1694</td></tr><tr><td>Adrar</td><td>20</td><td>5.06</td><td>3.53</td><td>0.58</td><td>0.56</td><td>0.004</td><td>0.5505</td></tr><tr><td>Beni Abbes</td><td>14</td><td>5.22</td><td>3.85</td><td>0.58</td><td>0.59</td><td>0.05</td><td>0.9013</td></tr></tbody></table>

The number of samples (N°), the number of alleles per locus (Na), allelic richness (Ar), observed heterozygosity (Ho), expected heterozygosity (He), theinbreeding coefficient (F<sub>IS</sub>) and the p-value of F<sub>IS</sub> are given for each population (averaged across the 18 loci).

https://doi.org/10.1371/journal.pone.0175232.t003
<table><caption>Table 4. Genetic distance between oases (F<sub>ST</sub>).</caption><thead><tr><th></th><th>Biskra</th><th>Touggourt</th><th>Oued Souf</th><th>Ouargla</th><th>Ghardaia</th><th>Tamanrasset</th><th>El Menia</th><th>Timimoun</th><th>Adrar</th><th>Beni Abbes</th></tr></thead><tbody><tr><td>Biskra</td><td></td><td><strong>0.0288*</strong></td><td>0.0288</td><td>0.0062</td><td><strong>0.0701*</strong></td><td>0.0205</td><td><strong>0.0668*</strong></td><td><strong>0.0581*</strong></td><td><strong>0.0559*</strong></td><td><strong>0.0623*</strong></td></tr><tr><td>Touggourt</td><td></td><td></td><td>0.0098</td><td>-0.0288</td><td><strong>0.0379*</strong></td><td><strong>0.0473*</strong></td><td><strong>0.0420*</strong></td><td><strong>0.0479*</strong></td><td><strong>0.0492*</strong></td><td><strong>0.0555*</strong></td></tr><tr><td>Oued Souf</td><td></td><td></td><td></td><td>-0.0065</td><td><strong>0.0240*</strong></td><td>0.0270</td><td>0.0289</td><td><strong>0.0322*</strong></td><td><strong>0.0334*</strong></td><td><strong>0.0313*</strong></td></tr><tr><td>Ouargla</td><td></td><td></td><td></td><td></td><td>0.0197</td><td>0.0049</td><td>0.0016</td><td>0.0063</td><td>0.0079</td><td>0.0132</td></tr><tr><td>Ghardaia</td><td></td><td></td><td></td><td></td><td></td><td>0.0322</td><td><strong>0.0414*</strong></td><td>0.0143</td><td><strong>0.0365*</strong></td><td><strong>0.0372*</strong></td></tr><tr><td>Tamanrasset</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.0245</td><td>0.0141</td><td>0.0129</td><td>0.0172</td></tr><tr><td>El Menia</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.0179</td><td>0.0072</td><td>0.0250</td></tr><tr><td>Timimoun</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.0119</td><td>0.0056</td></tr><tr><td>Adrar</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.0086</td></tr><tr><td>Beni Abbes</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

The F<sub>ST</sub> component was determined using a method that calculates the genetic distance between two oases, the differences being in pairs. A Bonferroni correction was applied, the significance level passing p < 0.001.

https://doi.org/10.1371/journal.pone.0175232.t004

## Population structure and differentiation

Analysis of population structure led us to keep three sub-populations (Fig 2) based on Evanno statistics (Fig 2b). Among the 192 accessions with distinct genotypes, 186 were assigned to a single group (Fig 2). Only six accessions remained "unclassified" and were considered "admixed". The first group included 51 varieties from the north-eastern and north-central part of the Algerian Sahara. The second group, which was the smallest with 48 varieties, was present in every oasis. The third group was the largest with 87 varieties, including varieties from the north-western, north-central and southern part of the Algerian Sahara, and 20% from the north-eastern part.

Allelic richness was lowest in the northern group (Table 5). No difference in observed and expected heterozygosity was found between the three groups (Table 5, Wilcoxon signed-rank test, p-value > 0.05). The inbreeding coefficient F<sub>IS</sub> was negative and significant in the first and third group (Table 5). Only group 2 had no excess of heterozygotes. Pairwise F<sub>ST</sub> values in the three groups ranged from 0.066 to 0.092 (Table 6). The highest differentiation was found between group 1 and group 2 (F<sub>ST</sub> = 0.0924). The three groups are thus genetically different (significantly different at p < 1%).

A relationship was found between the three genetic structure clusters and the two chloro-type frequencies (Fig 2, χ<sup>2</sup> = 28.1, dof = 2, p<0.001).

<table><caption>Table 5. Description and variability parameters of the 18 polymorphic microsatellite loci in the three STRUCTURE populations.</caption><thead><tr><th>Population</th><th>N'</th><th>Na</th><th>Ar</th><th>Ho</th><th>He</th><th>F<sub>IS</sub></th><th>P value F<sub>IS</sub></th></tr></thead><tbody><tr><td>Population 1</td><td>51</td><td>5.28</td><td>5.15</td><td>0.60</td><td>0.53</td><td>-0.13</td><td>0.000002</td></tr><tr><td>Population 2</td><td>48</td><td>6</td><td>5.88</td><td>0.57</td><td>0.56</td><td>-0.004</td><td>0.426301</td></tr><tr><td>Population 3</td><td>87</td><td>6.94</td><td>6.32</td><td>0.61</td><td>0.58</td><td>-0.05</td><td>0.000101</td></tr></tbody></table>

The number of samples (N'), the number of alleles per locus (Na), allelic richness (Ar), expected heterozygosity (He), observed heterozygosity (Ho), the inbreeding coefficient (F<sub>IS</sub>) and the p-value of F<sub>IS</sub> are given for each population and each locus.

https://doi.org/10.1371/journal.pone.0175232.t005

<footnote>* P < 0.001 (significant).</footnote>
<table><caption>Table 6. Genetic distance between Structure populations (F<sub>ST</sub>).</caption><thead><tr><th></th><th>population 1</th><th>population 2</th><th>population 3</th></tr></thead><tbody><tr><th>population 1</th><td></td><td>0.0924*</td><td>0.0742*</td></tr><tr><th>population 2</th><td></td><td></td><td>0.0659*</td></tr><tr><th>population 3</th><td></td><td></td><td></td></tr></tbody></table>

The F<sub>ST</sub> component was determined using a method that calculates the genetic distance between two assumed populations, the differences being in pairs. A Bonferroni correction was applied, the significance level passing p < 0.01.

https://doi.org/10.1371/journal.pone.0175232.t006

## Genetic distance and dendrogram construction

The matrix of estimated distances between individuals of the total population was used to construct a phylogenetic tree of individuals in the total population. The distances ranged from 0 to 0.75, underlining the wide genetic variability of our study population. Zero distance between two individuals suggests a clonal relationship (Fig 3).

Analysis of the phylogenetic tree showed that the individuals were grouped independent of their geographic origin and ethnic name ecotype. The distribution on the phylogenetic tree

<figure><img src="image_6.png" /><figcaption>**Fig 3. Neighbor-joining tree with microsatellite genotypes using shared allele distance.** The neighbor-joining tree shows the genetic relationships among the date palm genotypes included in the study, where each branch represents a single individual. The individuals in the population STRUCTURE1 are in red, those in population 2 in green, and those in population 3 in blue.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0175232.g003

<footnote>* p < 0.01 (significant)</footnote>
can be explained by the existence of a common genetic basis between different populations
and ecotypes, despite geographic distance and phenotypic divergence.

Discussion

Genetic diversity and population structure of date palm

The diversity of Algerian date palm was shown to be high with an average number of alleles per locus of 7.94, similar to the number observed in Qatar, Iraq or Tunisia [34, 35, 36, 37, 38]. Seven populations had negative fixation indices (FIS), suggesting an excess of heterozygosity at these loci. Two populations, Touggourt (FIS = -0.06) and Ghardaia (FIS = -0.13), had significant positive fixation indices (FIS). A negative FIS value suggests that different heterozygous genotypes have been maintained. Negative and significant FIS value were also found when individuals were grouped based on Bayesian analysis of population structure. This particular signature is not very common in plants, and may indicate direct selection for individuals with the best performance, i.e. those with high heterozygosity. It suggests farmers had a direct impact on maintaining the most heterozygous plants [39]. As date palms are also propagated by stem (offshoots), it should be noted that if all the plants were at Hardy-Weinberg equilibrium, propagation by stem would not change or increase FIS through the production of clones. So this excess of heterozygous plants could only have been be maintained by choosing the fittest individual [39].

The difference in F<sub>ST</sub> between oases was relatively low, ranging from 0.0016 to 0.0701.
However, slightly higher differentiation was observed between oases in the northern and
southern Algerian Sahara. This could reflect the use of different varieties in different ecological
conditions. In the Algerian Sahara, the climate is characterized by aridity and heat, and gradu-
ally becomes hotter and drier from north to south. The northern regions of the Sahara are
characterized by a semi-arid climate, which is the case of Biskra located at the foot of the Saha-
ran Atlas mountain range, the natural boundary between northern and southern Algeria. Simi-
larly, Touggourt, Oued Souf, Ouargla and Ghardaia havea semi-arid to arid climate. Southern
Saharan regions, i.e. Beni Abbes, Adrar, Timimoun and Tamanrasset are characterized by an
extremely hot dry desert climate. The analyses of structure also support the northern and
southern grouping of diversity.

Explanation for the genetic history of diffusion of date palm cultivars in
Algeria

Recent studies based on the global diversity of date palm suggest a geographic differentiation between African and Middle Eastern date palms [10, 8]. Chloroplastic diversity and Pairwise the non-recombinant portion of the segment of the Y chromosome both reveal differentiation between African and Middle Eastern cultivars [11]. A hypothesis of two origins, one in the east and the other in the west has been proposed [10, 8]. However, it is likely that this differentiation was caused by successive bottlenecks during the diffusion of cultivated date palms. However that may be, the diversity of the chloroplast is structured with major eastern and western alleles. Chloroplast genomes help trace the path of the female lines [7], and our chloroplast data provided insight into this maternal structure.

Interestingly, we observed a high frequency (almost 75%) of the eastern (oriental) chloro-
plast in a specific group of varieties. In this particular group, the frequency was similar to that
found further east. Conversely, in the two other structure groups, chloroplastic allele frequency
was similar to that observed further west.
The 18 microsatellite loci used in the study differed between eastern and western date palms [10]. Algeria does not differ particularly from other western countries and falls in the western genetic group. While chloroplasts inherited maternally suggest some proximity to the east, the nuclear markers clearly show a proximity to the west. One hypothesis proposed to explain this peculiar pattern, is that seeds and offshoots were imported from the Middle East to North Africa, and, after generations of crossing with local males (western nuclear genome), varieties acquired a western nuclear genome although an eastern chlorotype was maintained. Such varieties therefore carry the eastern haplotype. Another possible explanation is that both eastern and western haplotypes were widespread in both the eastern and western Sahara, but because farmers selected a specific genetic group of plants, chloroplastic diversity varied considerably from one local genetic group to the other in the west. Whatever the scenario, this pattern is only compatible with human selection that maintained the maternal lineage separately. The fact that we also found that in Algeria, nuclear SSR is structured in three groups shows that the groups were also to a certain extent kept separate by controlled pollen flow.

In conclusion, humans shaped the diversity of date palms by selecting heterozygous individuals and maintaining them. Moreover, cryptic diversity was observed both at genome and chloroplast level. Taken together, our results suggest that human selection played a major role in maintaining the diversity and lineage of date palms in Algeria.

**Supporting information**

S1 Fig. **Clustering of individuals by STRUCTURE at K = 2 and K = 3.** Individuals are represented by vertical colored lines. Individuals of the same color belong to the same cluster. Individuals with several different colors show the percentage of the genome that was inherited from each cluster.

(TIF)

S1 Table. **List of samples used in nuclear SSR genotyping.**

(PDF)

S2 Table. **Genetic assignment of the sample cultivars according to chlorotype and STRUCTURE grouping for K = 3.**

(PDF)

S3 Table. **Linkage disequilibrium results for 18 markers.**

(PDF)

S4 Table. **P-value of allelic richness between oases calculated by the Wilcoxon test.**

(PDF)

S5 Table. **P-value of expected heterozygosity between oases calculated by the Wilcoxon test.**

(PDF)

S6 Table. **P-value of observed heterozygosity calculated by the Wilcoxon test between oases.**

(PDF)

S7 Table. **P-value of Fis calculated by the Wilcoxon test between oases.**

(PDF)

S8 Table. **P-value of allelic richness between Structure populations calculated by the Wilcoxon test.**

(PDF)
**S9 Table.** P-value of expected heterozygosity between Structure populations calculated by the Wilcoxon test.
(PDF)

**S10 Table.** P-value of observed heterozygosity between Structure populations calculated by the Wilcoxon test.
(PDF)

**S11 Table.** P-value of FIS between Structure populations calculated by the Wilcoxon test.
(PDF)

**S12 Table.** Geographic distance (in km.) between oases.
(PDF)

**Acknowledgments**

JC Pintaud passed away during the course of the study. This work is dedicated to his memory.
The publication of this manuscript was partly supported by the Agropolis Resource Center for Crop Conservation, Adaptation and Diversity (ARCAD) funded by the Agropolis Fondation, Montpellier, France.

**Author Contributions**

**Conceptualization:** SM JCP NB.

**Data curation:** SM YV.

**Formal analysis:** SM YV.

**Funding acquisition:** SM.

**Investigation:** SM JCP NB.

**Methodology:** SM JCP.

**Project administration:** SM NB JCP.

**Resources:** SM NB JCP YV.

**Software:** SM YV JCP.

**Supervision:** SM NB YV JCP.

**Validation:** SM JCP YV.

**Visualization:** SM YV.

**Writing – original draft:** SM YV.

**Writing – review & editing:** SM YV NB.

**References**

1. Munier P. Le palmier-dattier. Maisonneuve et Larose, Coll. Techniques Agricoles et Productions Tropicales n° XXIV, Paris. 1973.

2. Barrow S. A monograph of *Phoenix* L. (Palmae, Coryphoideae). Kew Bulletin. 1998; 53: 513–575.

3. Dransfield J, Uhl NW, Asmussen CB, Baker WJ, Harley MM, Lewis CE. Genera Palmarum. The evolution and classification of palms. Royal Botanic Gardens, Kew, Richmond, U.K. 2008; 732 p.
4. Chao CCT, Krueger RR. *Phoenix dactylifera*, date palm. 2008; pp. 138–151 in Janick J, Paull RE (eds.): The encyclopedia of fruit and nuts. CAB International, Cambridge University Press, Cambridge, U.K.

5. Terral JF, Newton C, Ivorra S, Gros-Balthazard M, Tito de Morais C, Picq S, et al. Insights into the historical biogeography of the date palm (*Phoenix dactylifera* L.) using geometric morphometry of modern and ancient seeds. Journal of Biogeography. 2012; 39: 929–941.

6. Arabnezhad H, Bahar M, Mohammadi HR, Latifian M. Development, characterization, and use of micro-satellite markers for germplasm analysis in date palm (*Phoenix dactylifera* L.). Scientia Horticulturae. 2012; 134: 150–156.

7. Pintaud JC, Ludeña B, Aberlenc-Bertossi F, Zehdi S, Gros-Balthazard M, Ivorra S, et al. Biogeography of the date palm (*Phoenix dactylifera* L., Arecaceae): insights on the origin and on the structure of modern diversity. Acta Horticulturae. 2013; 994: 19–38.

8. Mathew LS, Seidel MA, George B, Mathew S, Spannagl M, Haberer G, et al. Genome-Wide Survey of Date Palm Cultivars Supports Two Major Subpopulations in *Phoenix dactylifera*. G3. 2015; 5:1429–1438. https://doi.org/10.1534/g3.115.018341 PMID: 25957276

9. Zehdi S, Cherif E, Rhouma S, Santoni S, Salhi HA, Pintaud JC. Molecular polymorphism and genetic relationships in date palm (*Phoenix dactylifera* L.): The utility of nuclear microsatellite markers. Scientia Horticulturae. 2012; 148: 255–263.

10. Zehdi S, Cherif E, Moussouni S, Gros-Balthazard M, Abbas Naqvi S, Ludeña B, et al. Genetic structure of the date palm (*Phoenix dactylifera* L.) in the Old World reveals a strong differentiation between eastern and western populations. Annals of Botany. 2015; 116: 101–112. https://doi.org/10.1093/aob/mcv068 PMID: 26113618

11. Cherif E, Zehdi S, Castillo K, Chabrillange N, Abdoulkader S, Pintaud JC, et al. Male-specific DNA markers provide genetic evidence of an XY chromosome system and recombination arrest, and allow to trace paternal lineages in date palm. New Phytologist. 2013; 197: 409–415. https://doi.org/10.1111/nph.12069 PMID: 23231423

12. Bouguedoura N, Bennaceur M, Babahani S, Benziouche SE. Date Palm Status and Perspective in Algeria. 2015; 4: 125–168. In: Al-Khayri JM, Jain SM, Johnson DV: Date Palm Genetic Resources and Utilization. Volume 1: Africa and the Americas.

13. Bouguedoura N, Bennaceur M, Benkhalifa A. Le palmier dattier en Algérie Situation, contraintes et apports de la recherche. Actes du 3e Séminaire du réseau AUF-BIOVEG «Biotechnologie du palmier du palmier dattier » Montpellier, France. 2008.

14. Bouguedoura N. Le palmier dattier: développement morphogénétique des appareils végétatif et reproducteur. Editions universitaires européennes. Book Editions Universitaires et Internationales. 2012.

15. Sambrook J, Fritsch EF, Maniatis T. Molecular cloning: a laboratory manual. Cold Spring Harbor Laboratory Press: Cold Spring Harbor, NY. 1989.

16. Henderson SA, Billotte N, Pintaud JC. Genetic isolation of Cape Verde Island *Phoenixatlantica* (Arecaceae) revealed by microsatellite markers. Conservation Genetics. 2006; 7: 213–223.

17. Pintaud JC, Zehdi S, Couvreur T, Barrow S, Henderson S, Aberlenc-Bertossi F, et al. Species delimitation in the genus *Phoenix* (Arecaceae) based on SSR markers, with emphasis on the identity of the Date Palm (*Phoenix dactylifera* L.). 2010; 267–286. In: Seberg O, Petersen G, Barfod A, Davis J (eds.): Diversity, phylogeny, and evolution in the Monocotyledons, Aarhus University Press, Denmark.

18. Scarcelli N, Barnaud A, Eiserhardt W, Treier U, Seveno M, d'Anfray A, et al. A set of 100 DNA primer pairs to study population genetics and phylogeny in Monocotyledons. PLoS ONE. 2011; 6(5): e19954. https://doi.org/10.1371/journal.pone.0019954 PMID: 21637837

19. Khan A, Khan IA, Heinze B, Azim MK. The chloroplast genome sequence of date palm (*Phoenix dactylifera* L. cv. 'Aseel'). Plant Molecular Biology Reporter. 2012; 30: 666–678.

20. Billotte N, Marseillac N, Brottier P, Noyer JL, Jacquemoud-Collet JP, Moreau C, et al. Nuclear microsatellite markers for the date palm (*Phoenix dactylifera* L.): characterization and utility across the genus *Phoenix* and in other palm genera. Molecular Ecology Notes. 2004; 4: 256–258.

21. Ludeña B, Chabrillange N, Aberlenc-Bertossi F, Adam H, Tregear JW, Pintaud JC. Phylogenetic utility of the nuclear genes Agamous 1 and phytochrome B in palms (Arecaceae): an example within Bactridiaceae. Annals of Botany. 2011; 108: 1433–1444. https://doi.org/10.1093/aob/mcr191 PMID: 21828068

22. Adam H, Marguerettaz M, Qadri R, Adroher B, Richaud F, Collin M, et al. Divergent expression patterns of miR164 and CUP-SHAPED COTYLEDON genes in palms and other Monocots: Implication for the evolution of meristem function in Angiosperms. Molecular Biology and Evolution. 2010; 28:1439–1454. https://doi.org/10.1093/molbev/msq328 PMID: 21135149

23. Chabrillange N, Aberlenc F, Pintaud JC, Ludena B. Microsatellite markers for the date palm (Phoenix dactylifera L.) diversity characterization. Unpublished.
24. Aberlenc-Bertossi F, Castillo K, Tranchant-Dubreuil C, Chérif E, Ballardini M, Abdoulkader S, et al. In silico mining of microsatellites in coding sequences of the date palm (Arecaceae) genome, characterization, and transferability. Applications in Plant Sciences. 2014; 2: 1300058.

25. Liu K, Muse SV. Power Marker: an integrated analysis environment for genetic marker analysis. Bioinformatics. 2005; 21: 2128–2129. https://doi.org/10.1093/bioinformatics/bti282 PMID: 15705655

26. Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984; 38: 1358–1370.

27. Goudet J. Fstat Software,v 2.9.3.2. 2002. http://www2.unil.ch/popgen/softwares/fstat.htm.

28. Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Montpellier, France: Laboratoire Génome, Populations, Interactions, CNRS UMR 5171, Université de Montpellier II. 2004.

29. R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2008. ISBN 3-900051-07-0, URL http://www.R-project.org.

30. Kumar S, Tamura K, Nei M. MEGA3: integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004; 5: 150–163. PMID: 15260895

31. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000; 155: 945–959. PMID: 10835412

32. Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software Structure: A simulation study. Molecular Ecology. 2005; 14: 2611–2620. https://doi.org/10.1111/j.1365-294X.2005.02553.x PMID: 15969739

33. Vigouroux Y, Glaubitz JC, Matsuoka Y, Goodman MM, Sánchez JG, Doebley J. Population structure and genetic diversity of New World maize races assessed by DNA microsatellites. American Journal of Botany. 2008; 95(10): 1240–1253. https://doi.org/10.3732/ajb.0800097 PMID: 21632329

34. Zehdi S, Sakka H, Rhouma A, Ould Mohamed Salem A, Marrakchi M, Trifi M. Analysis of Tunisian date palm germplasm using simple sequence repeat primers. African journal of Biotechnology. 2004; 3: 215–219.

35. Khierallah HSM, Bader SM, Baum M, Hamwieh A. Genetic diversity of Iraqi Date palms revealed by microsatellite polymorphism. Journal of the American Society for Horticultural Science. 2011; 136: 282–287.

36. Elmeer K, Sarwath H, Malek J, Baum M, Hamwieh A. New microsatellite markers for assessment of genetic diversity in date palm (*Phoenix dactylifera* L.). 3 Biotec. 2011; 1:91–97.

37. Elmeer K, Mattat I. Marker-assisted sex differentiation in date palm using simplesequence repeats. 3 Biotech. 2012; 2:241–247.

38. Hamza H, Ben Abederrahim MA, Elbekkay M, Ferchichi A. Comparison of the effectiveness of ISSR and SSR markers in determination of date palm (*Phoenix dactylifera* L.) agronomic traits. Australian Journal of Crop Science. ISSN: 1835. 2013; 7(6):763–769.

39. Pujol B, Patrice David P, McKey D. Microevolution in agricultural environments: how a traditional Amer-Indian farming practice favours heterozygosity in cassava (Manihot esculenta Crantz, Euphorbiaceae). Ecology Letters. 2005; 8: 138–147.