# Genetic Variation in Damaged Populations of Pistacia atlantica Desf.

Labdelli Amina, Roberto de La Herran, Rami Arafeh, Francesca Resentin, Livio Trainotti, Youcef Halis, Ahmed Haddad, Othmane Merah

► To cite this version:

Labdelli Amina, Roberto de La Herran, Rami Arafeh, Francesca Resentin, Livio Trainotti, et al. Genetic Variation in Damaged Populations of Pistacia atlantica Desf.. Plants, 2020, 9 (11), pp.1541. 10.3390/plants9111541. [hal-03009509](https://hal.inrae.fr/hal-03009509)

HAL Id: hal-03009509

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

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

Article

# Genetic Variation in Damaged Populations of *Pistacia atlantica* Desf.

Amina Labdelli<sup>1,2 ID</sup>, Roberto De La Herrán<sup>3,* ID</sup>, Rami Arafeh<sup>4</sup>, Francesca Resentini<sup>5,6,* ID</sup>, Livio Trainotti<sup>6 ID</sup>, Youcef Halis<sup>1</sup>, Ahmed Adda<sup>2</sup> and Othmane Merah<sup>7,8,* ID</sup>

Received: 28 September 2020; Accepted: 9 November 2020; Published: 11 November 2020

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

check for
updates

**Abstract:** The Atlas Pistachio tree, *Pistacia atlantica* Desf., has great importance in the ecological landscape of North Africa, due to its adaptive plasticity, as well as its use as a rootstock in the cultivation of the economically important species, *Pistacia vera* L. The conservation and valuation of this species require sampling and an assessment of its genetic variability. For the first time in North Africa, the inter-simple sequence repeats (ISSR) molecular marker has been used in genetic-diversity assessment and in the population relationships of *P. atlantica* subsp. atlantica. The ISSR markers tested showed 74.1% polymorphism, while molecular variance (AMOVA) analysis revealed a high percentage of the total genetic diversity of 55.7% among the four populations studied. Cluster analysis with neighbor-joining (NJ) and principal coordinate analysis (PCO) divided the study sites into four distinct groups according to their geographical locations (Tiaret, Batna, Djelfa, and Bechar). Isolation by distance or Mantel test gave a positive correlation of *r* = 0.86 between geographical and genetic distances. The results in this study indicate an absence of gene flow, implying that conservation efforts should be taken separately for each population.

**Keywords:** Atlas pistachio; population analyses; endangered species; genetic conservation; ISSR variation

## 1. Introduction

The genus *Pistacia*, Anacardiaceae family consists of at least 11 dioecious species [1], four of these being found in North Africa, including Algeria (*Pistacia vera* L., *P. terebinthus* L., *P. lentiscus* L., and *P. atlantica* Desf.). Pistachio (*P. vera* L.) is known for its economic value in agriculture as a nut crop [2–4]. The other three species possess ecological importance in addition to their use as medicinal plants and sources of feed for animals and birds [4]. The Atlas pistachio (*P. atlantica* Desf.) is used for its leaves and fruits that have high nutritional value, and the wood is also used in handicrafts,

<footnote n="1">¹ Scientific and Technical Research Centre for Arid Areas (CRSTRA), BP 1682 RP, Biskra 07000, Algeria; aminalabdelli@yahoo.fr (A.L.); youcefhal@gmail.com (Y.H.)</footnote>

<footnote n="2">² Laboratory of Agro-Biotechnology and Nutrition in Semi-Arid Areas, Ibn Khaldon University, Tiaret 14000, Algeria; adda2ahmed@yahoo.fr</footnote>

<footnote n="3">³ Departamento de Genética, Facultad de Ciencias, Universidad de Granada, Avda, Fuentenueva s/n, 18071 Granada, Spain</footnote>

<footnote n="4">⁴ Palestine-Korea Biotechnology Center, Palestine Polytechnic University, P.O. Box (198), Hebron, Palestine; arafeh@ppu.edu</footnote>

<footnote n="5">⁵ Department of Biosciences, Università degli Studi di Milano, Via Giovanni Celoria 26, 20133 Milano, Italy</footnote>

<footnote n="6">⁶ Department of Biology, University of Padova, Viale Giuseppe Colombo 3, 35121 Padova, Italy; livio trainotti@unipd.it</footnote>

<footnote n="7">⁷ Laboratoire de Chimie Agro-Industrielle (LCA), Université de Toulouse, INRA, INPT, 31030 Toulouse, France</footnote>

<footnote n="8">⁸ Département Génie Biologique, Université Paul Sabatier, IUT A, 32000 Auch, France</footnote>

<footnote>* Correspondence: rherran@ugr.es (R.D.L.H.); francesca.resentini@unimi.it (F.R.); othmane.merah@ensiacet.fr (O.M.); Tel.: +33-534323523 (O.M.)</footnote>
in carpentry, and for firewood [2–5]. *Pistacia atlantica* Desf. is an Irano-Turanian tree distributed from south-west Asia to north-west Africa [1]. There are four subspecies in *P. atlantica* have been recognized, subsp. *atlantica*, *kurdica*, *mutica*, and *cabulica*. The subsp. *atlantica* in North Africa is found in a wide distribution range from sub-humid environments to extremely arid Saharan sites [6]. Being xerophyte makes *P. atlantica* Desf. a prime choice for establishing pistachio cultivation in Algeria because it is used as a rootstock for the common pistachio crop [1]. This species plays important ecological roles because it combats soil erosion and prevents landslides when used for the reforestation of arid and steep slopes. Unfortunately, in many regions, *P. atlantica* is facing genetic erosion after being subjected to overgrazing, bush fires, deforestation, desertification, pollution, and other adverse factors. Consequently, it is listed among the endangered species in many regions [7,8]. The genetic relationships based on morphological or molecular markers have been performed on the genus *Pistacia* to determine their taxonomy [1,6–12]. Among the species of this genus, the Atlas pistachio tree is characterized by a remarkable morphological variability that allows it to adapt to different ecological conditions [6–9]. Based on this fact, it is urgent to proceed with an assessment of the genetic variability of this species to adopt the appropriate strategies for its preservation. It is becoming of greater importance to explore genetic diversity and population relationships of endangered plant species, especially those with slow growth and high morphological and biochemical variation [3,7,9], as in the case of Atlas pistachio.

Despite the importance of this species, very few genetic studies have been carried out on its populations or in relation to other *Pistacia* species [13–15]. Several molecular markers have been used to determine genetic diversity and relations among *Pistacia* species and cultivars. The randomly amplified polymorphic DNA (RAPD) was among the most commonly used methods in pistachio to study genetic diversity [13,16–19] for being quick and easy to develop. However, it lacks reproducibility—so today, the method has fallen into disuse. Ibrahim Basha et al. [20] detected high diversity in *P. vera* L. by amplified fragment-length polymorphisms (AFLP) markers in Syria—a country considered one of the centers of origin for this species. Other markers used in different studies in *Pistacia* ssp. [19,21–23] include the simple sequence repeats (SSR), which are specific and highly polymorphic markers [24], but require prior knowledge of the genomic sequence to design specific primers, and therefore, its development is limited to more economically important species, and in many cases, cross-amplification is ineffective. Furthermore, other methods have been used, such as the sequence-related amplified polymorphism (SRAP), the selectively amplified microsatellite polymorphic loci (SAMPL), and the inter-simple sequence repeats (ISSR) as well [15,25–30].

The molecular marker, called the inter-simple sequence repeat (ISSR), which provides high repeatability, has been used in systematic comparisons below the species level and varietal identification [20,27]. The technique does not require prior information about the genome sequence and leads to multi-locus and highly polymorphic patterns [25–27]. Also, it involves longer primers (16–18 nucleotides) encoding microsatellite elements that amplify DNA segment intra-microsatellite repeats [31]. Here, such a study is conducted for the first time on *P. atlantica* Desf. Algerian populations, based on the ISSR marker.

The aim of the present study was to examine the intraspecific genetic diversity of four populations of *P. atlantica* Desf. growing in different zones of Algeria. This is intended to provide insight into the relationships between the study populations and to provide data for the proper conservation of this tree, as well as for future selection programs of rootstock variability for pistachio cultivation.

## 2. Results

Two primer pairs (ISSR1 and ISSR 2) did not produce any polymorphic fragments with different annealing temperatures in the polymerase chain reaction (PCR) for any of the different *P. atlantica* Desf. samples. For the other two primer pairs, a total of 112 DNA fragments were amplified, 72 of which (64.3%) were polymorphic. The number of polymorphic fragments varied between the populations and the ISSR primers, showing a range from 4 (ISSR4 in Bechar) to 14 (ISSR3 in Tiaret). The percentage of polymorphic loci was 91% on average, and the highest degree of polymorphism
was observed among wild populations of Tiaret (76.9%), whereas the lowest was found in the Djelfa population (50%; Table 1).

<table><caption>**Table 1.** The number of amplified and polymorphic bands observed at the different locations. ISSR, inter-simple sequence repeats.</caption><thead><tr><th>Primer</th><th></th><th>Batna</th><th>Djelfa</th><th>Tiaret</th><th>Bechar</th></tr></thead><tbody><tr><td rowspan="3">ISSR3</td><td>Total number of bands</td><td>19</td><td>14</td><td>19</td><td>13</td></tr><tr><td>Number of polymorphic bands</td><td>13</td><td>7</td><td>14</td><td>9</td></tr><tr><td>% polymorphic bands</td><td><strong>68.42</strong></td><td><strong>50</strong></td><td><strong>73.68</strong></td><td><strong>69.23</strong></td></tr><tr><td rowspan="3">ISSR4</td><td>Total number of bands</td><td>14</td><td>14</td><td>13</td><td>6</td></tr><tr><td>Number of polymorphic bands</td><td>8</td><td>7</td><td>10</td><td>4</td></tr><tr><td>% polymorphic bands</td><td><strong>57.14</strong></td><td><strong>50</strong></td><td><strong>76.92</strong></td><td><strong>66.67</strong></td></tr></tbody></table>

In the four natural populations of *P. atlantica* Desf. the number of effective alleles proved similar, ranging from 1.440 (Djelfa) to 1.64 (Tiaret). Expected heterozygosis and the genetic diversity parameter by Shannon's information index (SI) values ranged from 0.233 (Djelfa) to 0.344 (Tiaret) and 3.6 (Bechar) to 5.34 (Djelfa), respectively (Table 2).

<table><caption>**Table 2.** Genetic diversity of the four populations of *P. atlantica* Desf. inferred by the two ISSR markers.</caption><thead><tr><th>Populations</th><th>Na</th><th>%P</th><th>Nd</th><th>Ne</th><th>He</th><th>SI</th></tr></thead><tbody><tr><td>Batna</td><td>83</td><td>63.64%</td><td>1.636 ± 0.085</td><td>1.504 ± 0.076</td><td>0.272 ± 0.039</td><td>5.26 ± 0.29</td></tr><tr><td>Djelfa</td><td>83</td><td>50%</td><td>1.500 ± 0.096</td><td>1.440 ± 0.086</td><td>0.233 ± 0.045</td><td>5.34 ± 0.36</td></tr><tr><td>Tiaret</td><td>83</td><td>75%</td><td>1.750 ± 0.078</td><td>1.640 ± 0.069</td><td>0.344 ± 0.036</td><td>5.32 ± 0.17</td></tr><tr><td>Bechar</td><td>33</td><td>68.42%</td><td>1.684 ± 0.110</td><td>1.565 ± 0.071</td><td>0.305 ± 0.051</td><td>3.6 ± 0.19</td></tr></tbody></table>

Na: Number of alleles; %P: Percentage of polymorphic loci; Nd: Number of different alleles; Ne: Number of effective alleles; He: Expected heterozygosity; and SI: Shannon's information index.

The analysis of the molecular variance (AMOVA) based on ISSR data showed high percentages of the total genetic diversity of *P. atlantica* Desf. populations distributed at this small spatial scale—with 55.7% among populations and 44.3% within populations. The pairwise PhiST value (coefficient: standard Jaccard distance transformation: d = 1 – s) was 0.56, indicating high differentiation among populations. Significant correlation between different populations was found in this study (p ≤ 0.001; pairwise PhiST values, standard Jaccard coefficient, distance transformation d = 1 – s) (B-group, D-group, T-group: 0.45 ***, 0.57 ***, 0.59 ***). An analysis using Structure software, performing runs from K = 2 up to K = 6, showed the genetic structure of the populations analyzed.

The AMOVA approach also generates FST values (one of the F-statistics, such as Weir and Cockerham's 98 measures [32], this parameter is calculated as the diversity among populations/total diversity. Values close to zero indicate little differentiation among populations (signifying mostly the genetic diversity within populations), while values close to one indicate high differentiation among populations. The principal coordinate analysis (PCO) clearly showed the four main groups from distinct geographical locations, Batna (B), Djelfa (D), Bechar (A), and Tiaret (T) (Figure 1A). Positions of the individual samples showed an overall homogeneity within the sampled site. On the other hand, isolation by distance or the Mantel test gave a positive correlation (r = 0.86, p = 0.131; Figure 1B). This positive r-value would indicate a certain degree of genetic isolation within a very narrow scale of the entire distribution range and also an indicator for genetic isolation between the four populations studied with an absence of gene flow between them, although they represent a narrow area at the eastern portion of the *P. atlantica* Desf. distribution range (Figure 1B). The ISSR data, used to find genetic relatedness among the genotypes, clustered the samples into four groups. Our results also indicated that *P. atlantica* Desf. subsp. atlantica clustering was related to the geographical locations and climatic conditions of its growing regions (Figure 1C). As a complementary study, an analysis
using Structure software, performing runs from *K* = 2 up to *K* = 6, showed the genetic structure of the populations analyzed. The estimated membership coefficients for each individual (Figure 2) is represented according to differentiation among populations, where different genetic structures for the four populations in different degrees can be observed.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** PCO Analysis, Correlation between genetic and geographic distances, and Midpoint rooted neighbor-joining tree. (A) Three-dimensional PCO analysis of ISSR data; (B) Correlation between genetic and geographic distances plotted with *Z* = 2119.05, *r* = 0.867 as calculated by the Mantel test; (C) Midpoint rooted neighbor-joining tree of samples from four populations of *P. atlantica* Desf. generated with ISSR data. Individuals from A1–A10 belong to the Batna population, A11–A21 to the Djelfa population, A22–A33 to the Tiaret population, and A34–A39 to the Bechar population.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Graphical representation of the estimated membership coefficients for each individual obtained from the Bayesian clustering analysis of genetic structure for *K* = 2–4 computed from ISSR. Each individual is represented by a bar broken into *K* colored segments. The populations of Batna, Djelfa, Tiaret, and Bechar are indicated with the numbers 1 to 4, respectively.</figcaption></figure>

**3. Discussion**

Due to the lack of genomic information of this species, these markers are highly useful, since they require no prior knowledge of the sequences, and therefore, the primers need not be site-specific because they are directed to a genomic region containing a complementary microsatellite motif [33,34]. By means of the ISSR analysis, which has proved to be simple and efficient with high reproducibility, we found pronounced genetic variation in Atlas pistachio populations among the samples studied, i.e., enough to characterize the different Algerian populations analyzed. Thus, the ISSR technique proved adequate for discriminating different *P. atlantica* Desf. populations. Both of the ISSR primers used here generated polymorphisms and proved sufficient to discriminate all four populations. The principal coordinate analysis (PCO) and neighbor-joining (NJ) analysis yielded congruent patterns in the determination of genetic differences between the populations that are grown in different regions in Algeria. This finding suggests the influence of climatic conditions and geographic distances on the diversity of different sites of *P. atlantica* Desf.

Moreover, El Zerey-Belaskri et al. [19], studying populations of this species in north-western Algeria, demonstrated that a certain genetic similarity was established between them through evaluation by the SSR marker. According to these authors, this genetic linkage would occur through the dissemination of pollen grains and/or seeds between geographically close populations. This data supports the results obtained in our study, since the populations studied by El Zerey-Belaskri et al. [19], were distributed along a northwest-northeast transect of about 230 km, while the populations studied here are at much greater geographic distances, preventing the most cases a genetic flow between them (Figure 3).

The genetic relatedness among the genotypes was investigated by neighbor-joining (NJ) cluster analysis. This grouped the genotypes into four major groups with increasing aridity: low and medium aridity (T: Tiaret, B: Batna, D: Djelfa) and high aridity (A: Bechar) (Figure 1C). It bears noting that the most severe environmental conditions found in Bechar correlated with the highest genetic diversity within individuals of this population. For this population, the genetic diversity indexes also were similar to the other populations when they were compared (Table 2). These results corroborate that the ISSRs are good markers to detect genetic variation among closely related variants and in crop classification [29]. Mediterranean forests, in general, tend to be poor in endemic species [35]. In Algeria, the climate is transitional between maritime in the north and semi-arid to arid in the middle and the south [6]. Bechar is located in the south-west of Algeria, characterized by an arid climate with Saharan tendency and irregular precipitation. The interaction of Mediterranean and
Saharan flora includes very high diversity, and therefore, deserves special conservation [8]. The high morphological and biochemical plasticity of *P. atlantica* Desf. in response to aridity may explain its wide ecological distribution in North Africa [6,7]. In this region, the ecological upheaval caused by deforestation and overgrazing has given rise to a serious dysfunction, which affects high-altitude forests in particular. Conservation of forests and forest species in the Mediterranean basin is a complex issue, given the wide range of ecological situations, diverse forest uses, and the pressure exerted by different cultural groups [35]. Based on the results of this study, conservation efforts should be focused on each population separately, due to the genetic isolation and absence of gene flow between the four populations. Additionally, the information on the genetic variability of these populations may be of interest at the time of selecting the specimens used as rootstocks in the cultivation of *P. vera* L.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Geographical locations of the four sampled populations of *P. atlantica* Desf. collected in Algeria. The distance between them are: Batna (B)-Tiaret (T) ~387 Km; Batna (B)-Djelfa (D) ~300 Km; Batna (B)-Bechar (A) ~899 Km; Tiaret (T)-Djelfa (D) ~198 Km; Tiaret (T)-Bechar (A) ~563 Km; and Djelfa (D)-Bechar (A) ~604 Km.</figcaption></figure>

# 4. Materials and Methods

## 4.1. Plant Materials

Leaf samples of 39 *P. atlantica* Desf. trees were collected from four different locations in Algeria (Table 3; Figures 3 and 4). A voucher specimen was deposited in the herbarium of the Scientific and Technical Research Centre for Arid Areas (CRSTRA), Biskra, Algeria (006CRSTRA0002). Details of the ecological factors for each collection site are reported in Labdelli et al. [3].

<table><caption>**Table 3.** Geographical locations and numbers of samples collected from each *P. atlantica* Desf. population analyzed.</caption><thead><tr><th>Assigned Code</th><th>Location</th><th>Number of Samples</th><th>Altitude (m)</th><th>Latitude (N)</th><th>Longitude</th></tr></thead><tbody><tr><td>B</td><td>Batna Semi-arid</td><td>10</td><td>1027</td><td>35°37'10"</td><td>6°22'13" E</td></tr><tr><td>T</td><td>Tiaret Semi-arid</td><td>12</td><td>808</td><td>35°22'33"</td><td>02°09'5" W</td></tr><tr><td>D</td><td>Djelfa Arid</td><td>11</td><td>630</td><td>34°02'11"</td><td>03°40'22" E</td></tr><tr><td>A</td><td>Bechar Hyperarid</td><td>06</td><td>979</td><td>32°04'6"</td><td>02°18'5" W</td></tr></tbody></table>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** *Pistacia atlantica* Desf. *atlantica* tree and leaves in regions of Batna (A) and Bechar (B) (September 2015). Photos: Labdelli A.</figcaption></figure>

## 4.2. DNA Extraction

The leaf samples were dried by silica gel and stored at -80 °C until used. Genomic DNA was extracted from leaf tissue by the CTAB (hexadecyltrimethyl ammonium bromide) method with the modifications adopted by Kafkas [36] and Al-Sousli et al. [37]. Young leaf tissues (100 mg) was manually ground into a fine powder in liquid nitrogen with a plastic pestle in a 2.0-mL Eppendorf tube and mixed with 900 µL of CTAB extraction buffer (1 M TRIS-HCl pH 8, 5 M NaCl, 0.5 M EDTA, 2% CTAB, 2% polyvinylpyrrolidone PVP, 0.2% b-mercaptoethanol, and 0.1% NaHSO₃).

Samples were incubated at 65 °C for 1 h, mixed with an equal volume of chloroform-isoamyl alcohol (24:1), and centrifuged for 5 min at 13,000 rpm. The aqueous phase was transferred to a new 2-mL Eppendorf tube, mixed with equal volume of ice-cold isopropanol to precipitate the DNA, and left at -80 °C for 1 h. The tubes were gently inverted several times. The nucleic acid pellet was recovered by centrifugation at 13,000 rpm for 5 min, washed with 1.0 mL of 10.0 mM ammonium acetate in 70% ethanol for a few min, dried for 15 min at room temperature, and resuspended in 100 µL modified TE buffer pH 8.0 (10 mM TRIS-HCl, 0.1 mM EDTA) added with RNAase (5 µg/mL). The concentration of the DNA was estimated by comparing band intensity with a DNA of known concentrations on 0.8% agarose gel in 1x TAE buffer. The DNA was diluted to 20 ng/µL for ISSR reactions.

## 4.3. ISSR-PCR Amplification

The samples were genotyped for ISSR markers using GC- and CA-rich primers. Amplifications were performed in 25 µL reaction mixture containing 1.7 µL of genomic DNA (20 ng), 0.8 µL primer (Table 4), 4.0 µL of reaction buffer (5x Green GoTaq® Reaction Buffer), 0.8 µL of MgCl₂ (50 mM) (or DMSO dimethyl sulfoxide + buffer 10x PFU Reaction w/20 Mm MgSO₄ Promega), 0.8 µL dNTPs, 0.2 µL of Taq DNA Polymerase (PFU DNA Polymerase 5 U/µL), 11.7 µL of high-performance liquid chromatography (HPLC) grade water. The thermocycler program included 1 cycle of initial DNA denaturation at 94 °C for 10 min, followed by 35 cycles of 45 s at 94 °C, 45 s annealing at 54 °C and 2 min extension at 72 °C, followed by a final extension for 10 min at 72 °C [37]. PCR reactions were run in agarose gel (3.0% w/v) in 1x TAE buffer and visualized by UV light.
<table><caption>**Table 4.** List of the ISSR primers used in this study, their sequences, and the annealing temperature applied to each. (W= A or T, Y= C or T, B= T, or C or G).</caption><thead><tr><th>Primer Name</th><th>Sequence 5'-3'</th><th>Annealing Temperature (°C)</th></tr></thead><tbody><tr><td>ISSR 1</td><td>(CG)<sub>9</sub>W</td><td>From 50 to 60</td></tr><tr><td>ISSR 2</td><td>(GC)<sub>9</sub>W</td><td>From 50 to 60</td></tr><tr><td>ISSR 3</td><td>(AG)<sub>9</sub>B</td><td>54</td></tr><tr><td>ISSR 4</td><td>(GA)<sub>9</sub>Y</td><td>54</td></tr></tbody></table>

## 4.4. Data Analysis

The ISSR amplicons (bands) were scored as present (1) or absent (0) according to their size, and then a binary data matrix was constructed. Only clear and reproducible bands were scored and used in the data analysis. Bands were assumed to be independent, and those of identical size were assumed to have identical sequences. The genetic diversity indices, including the number of alleles (Na), percentage of polymorphic loci (%P), number of different alleles (Nd), number of effective alleles (Ne), expected heterozygosity (He), and Shannon's information index (SI), were assessed at each population levels using FAMD Software version 1.31 [38] and GenAIEx Software [39].

The similarity matrix was also used to perform a hierarchical analysis of molecular variance (AMOVA), according to Excoffier et al. [40] using fingerprinting analysis with missing data (FAMD) Software version 1.31 [38]. This analysis enables the partitioning of the total ISSR variation into within and among geographical region variation components (called Phi statistics).

The genetic structure of populations was simulated by a Bayesian clustering analysis with the software STRUCTURE, version 2.3 [41]. The genotypic data of populations was analyzed under a clustering model of admixture of individuals with correlated allele frequencies among populations. Analysis of ISSR data included a model of clusters from $K = 2$ to $K = 4$.

The genetic distance matrix was estimated based on Jaccard's similarity coefficient using the multi-locus fingerprinting data sets containing missing data [FAMD software, version 1.31]. A cluster analysis was performed using the neighbor joining (NJ) method with arithmetic averages and viewed with Tree View software (Win32, version 1.6.6). PCO analysis was made with FAMD version 1.3 (Fingerprinting analysis with missing data) using the Jaccard's similarity index [42].

To test for isolation by distance [40], pairwise log₁₀ transformed values of Nₑ^m among populations were linearly regressed against their geographical (=aerial map) distances (in km). The significance of the regression slopes was evaluated by Mantel's test [43] with 10,000 random permutations using the Isolation by Distance Web Service version (3.23) under the link http://ibdws.sdsu.edu/~ibdws.

## 5. Conclusions

This work shows that Algeria still contains a rich genetic diversity for the Atlas Pistachio, which is likely to be valued. Despite the strict allogamy of this species and its highly dynamic gene flow, the geographical remoteness of the study areas indicates that this species still maintains a diversity of genetic entities. This inter-population diversity would certainly be remodeled by natural selection under the diversity of the environmental conditions of the different prospected areas. The results also indicate that the rate of intra-population genetic variability is negatively correlated with the increase in density of the populations studied.

**Author Contributions:** A.L., F.R., and L.T. designed the methodology; R.A., A.L. performed the data analysis; A.L., F.R., R.D.L.H., Y.H., A.A., and O.M. software; R.D.L.H., Y.H., A.A., and O.M. participated in writing and original draft preparation; A.A. and O.M. project administration. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Conflicts of Interest:** The authors declare no conflict of interest.
## References

1. Zohary, M. A monographical study of the genus *Pistacia*. *Palestine J. Bot.* **1952**, *5*, 187-228.

2. Labdelli, A.; Zemour, K.; Simon, V.; Cerny, M.; Adda, A.; Merah, O. *Pistacia atlantica* Desf., a Source of Healthy Vegetable Oil. *Appl. Sci.* **2019**, *9*, 2552. [CrossRef]

3. Labdelli, A.; Rebiai, A.; Tahirine, M.; Adda, A.; Merah, O. Nutritional Content and Antioxidant Capacity of the Seed and the Epicarp in Different Ecotypes of *Pistacia atlantica* Desf. subsp. atlantica. *Plants* **2020**, *9*, 1065. [CrossRef] [PubMed]

4. Ben Ahmed, Z.; Yousfi, M.; Viaene, J.; Dejaegher, B.; Demeyer, K.; Heyden, Y.V. Four *Pistacia atlantica* subspecies (atlantica, cabulica, kurdica and mutica): A review of their botany, ethnobotany, phytochemistry and pharmacology. *J. Ethnopharmacol.* **2020**, *265*, 113329. [CrossRef] [PubMed]

5. Belhadj, S.; Derridj, A.; Gauquelin, T. Pistachio use in Algeria. In: Avanzato D, Vassallo I. (eds) Following pistachio footprints (*Pistacia vera* L.). Cultivation and culture, folklore and history, traditions and uses. *Script. Hortic.* **2008**, *7*, 15-21.

6. El Zerey-Belaskri, A. Taxonomic and botanical retrospective review of *Pistacia atlantica* Desf. (Anacardiaceae). *Arab. J. Med. Aromat. Plants* **2019**, *5*, 47-77. [CrossRef]

7. Said, S.A.; Fernandez, C.; Greff, S.; Derridj, A.; Gauquelin, T.; Mevy, J.-P. Inter-population variability of leaf morpho-anatomical and terpenoid patterns of *Pistacia atlantica* Desf. ssp. atlantica growing along an aridity gradient in Algeria. *Flora Morphol. Distrib. Funct. Ecol. Plants* **2011**, *206*, 397-405. [CrossRef]

8. Amara, M.; Bouazza, M.; Al Saghir, B.A.M. Anatomical and Adaptation Features of *Pistacia atlantica* Desf. to Adverse Climate Conditions in Algeria. *Am. J. Plant Sci.* **2017**, *8*, 137-153. [CrossRef]

9. Zerey-Belaskri, E.; Benhassaini, H. Morphological leaf variability in natural populations of *Pistacia atlantica* Desf. subsp. atlantica along climatic gradient: New features to update *Pistacia atlantica* subsp. atlantica key. *Int. J. Biometeorol.* **2015**, *60*, 577-589. [CrossRef]

10. Al-Saghir, M.G.; Porter, D.M. Taxonomic Revision of the Genus *Pistacia* L. (Anacardiaceae). *Am. J. Plant Sci.* **2012**, *3*, 12-32. [CrossRef]

11. Shahghobadi, H.; Shabanian, N.; Rahmani, M.-S.; Khadivi-Khub, A. Genetic characterization of *Pistacia atlantica* subsp. kurdica from northern Zagros forests in Iran. *Trees* **2018**, *33*, 481-490. [CrossRef]

12. Xie, L.; Yang, Z.-Y.; Wen, J.; Li, D.-Z.; Yi, T.-S. Biogeographic history of Pistacia (Anacardiaceae), emphasizing the evolution of the Madrean-Tethyan and the eastern Asian-Tethyan disjunctions. *Mol. Phylogenet. Evol.* **2014**, *77*, 136-146. [CrossRef] [PubMed]

13. Golan-Goldhirsh, A.; Barazani, O.; Wang, Z.S.; Khadka, D.; Saunders, J.A.; Kostiukovsky, V.; Rowland, L.J. Genetic relationships among Mediterranean *Pistacia* species evaluated by RAPD and AFLP markers. *Plant Syst. Evol.* **2004**, *246*, 9-18. [CrossRef]

14. Arefi, H.; Abdi, A.; Saydian, S.; Nasirzadeh, A.; Nadushan, H.; Rad, M.; Golbabaii, H.; Azdoo, Z.; Ziedabadi, D. Genetics and breeding of pistacia atlantica in Iran. *Acta Hortic.* **2006**, *726*, 77-84. [CrossRef]

15. Karimi, H.R.; Kafkas, S. Genetic relationships among Pistacia species studied by SAMPL markers. *Plant Syst. Evol.* **2011**, *297*, 207-212. [CrossRef]

16. Kafkas, S.; Perl-Treves, R. Morphological and molecular phylogeny of Pistacia species in Turkey. *Theor. Appl. Genet.* **2001**, *102*, 908-915. [CrossRef]

17. Katsiotis, A.; Hagidimitriou, M.; Drossou, A.; Pontikis, C.; Loukas, M. Genetic relationships among species and cultivars of *Pistacia* using RAPDs and AFLPs. *Euphytica* **2003**, *132*, 279-286. [CrossRef]

18. Mirzaei, S.; Bahar, M.; Sharifnabi, B. A phylogenetic study of iranian wild pistachio species and some cultivars using rapd markers. *Acta Hortic.* **2006**, *726*, 39-44. [CrossRef]

19. El Zerey-Belaskria, A.; Ribeirob, T.; Alcarazc, M.L.; El Zereyd, W.; Castroe, S.; Loureiroe, J.; Benhassainia, H.; Hormaza, J.I. Molecular characterization of *Pistacia atlantica* Desf. subsp. atlantica (Anacardiaceae) in Algeria: Genome size determination, chromosome count and genetic diversity analysis using SSR markers. *Sci. Hortic.* **2018**, *227*, 278-287. [CrossRef]

20. Basha, A.I.; Padulosi, S.; Chabane, K.; Hadj-Hassan, A.; Dulloo, E.; Pagnotta, M.A.; Porceddu, E. Genetic diversity of Syrian pistachio (*Pistacia vera* L.) varieties evaluated by AFLP markers. *Genet. Resour. Crop Evol.* **2007**, *54*, 1807-1816. [CrossRef]
21. Pazouki, L.; Mardi, M.; Shanjani, P.S.; Hagidimitriou, M.; Pirseyedi, S.M.; Naghavi, M.R.; Avanzato, D.; Vendramin, E.; Kafkas, S.; Ghareyazie, B.; et al. Genetic diversity and relationships among Pistacia species and cultivars. *Conserv. Genet.* **2010**, *11*, 311–318. [CrossRef]

22. Arabnezhad, H.; Bahar, M.; Pour, A.T. Evaluation of genetic relationships among Iranian pistachios using microsatellite markers developed from Pistacia khinjuk Stocks. *Sci. Hortic.* **2011**, *128*, 249–254. [CrossRef]

23. Motalebipour, E.Z.; Kafkas, S.; Khodaeiaminjan, M.; Çoban, N.; Gözel, H. Genome survey of pistachio (*Pistacia vera* L.) by next generation sequencing: Development of novel SSR markers and genetic diversity in Pistacia species. *BMC Genom.* **2016**, *17*, 998. [CrossRef] [PubMed]

24. Khadivi-Khub, A.; Esmaeili, A.; Mardani, N. Genetic diversity of cultivated pistachio as revealed by microsatellite molecular markers. *Biotechnol. Biotechnol. Equip.* **2018**, *32*, 602–609. [CrossRef]

25. Kafkas, S.; Ozkan, H.; Ak, B.E.; Acar, I.; Atli, H.S.; Koyuncu, S. Detecting DNA Polymorphism and Genetic Diversity in a Wide Pistachio Germplasm: Comparison of AFLP, ISSR, and RAPD Markers. *J. Am. Soc. Hortic. Sci.* **2006**, *131*, 522–529. [CrossRef]

26. Pourian, M.A.; Bakhshi, D.; Aalami, A.; Hokmabadi, H. Assessment of Genetic Relationships among Cultivated and Wild Pistachios (*Pistacia vera* L.) using Molecular Markers. *J. Hortic. Res.* **2019**, *27*, 37–46. [CrossRef]

27. Wu, W.; Chen, F.; Yeh, K.-W.; Chen, J. ISSR Analysis of Genetic Diversity and Structure of Plum Varieties Cultivated in Southern China. *Biology* **2018**, *8*, 2. [CrossRef]

28. Fares, K.; Guasmi, F.; Touil, L.; Triki, T.; Ferchichi, A. Genetic Diversity of Pistachio Tree using Inter-Simple Sequence Repeat Markers ISSR Supported by Morphological and Chemical Markers. *Biotechnology* **2009**, *8*, 24–34. [CrossRef]

29. Tagizad, A.; Ahmadi, J.; Haddad, R.; Zarrabi, M. A comparative analysis of ISSR and RAPD markers for studying genetic diversity in Iranian pistachio cultivars. *IJGPB* **2010**, *1*, 6–16.

30. Turkeli, Y.; Kafkas, S. First genetic linkage map in pistachio constructed using an interspecific cross between *Pistacia vera* L. and monoecious *Pistacia atlantica* Desf. *Sci. Hortic.* **2013**, *151*, 30–37. [CrossRef]

31. Gupta, M.; Chyi, Y.S.I.; Romero-Severson, J.; Owen, J.L. Amplification of DNA markers from evolutionary diverse genomes using single primers of simple sequence repeats. *Theor. Appl. Genet.* **1994**, *89*, 998–1006. [CrossRef] [PubMed]

32. Weir, B.S.; Cockerham, C.C. Estimating F-Statistics for the Analysis of Population Structure. *Evolution* **1984**, *38*, 1358. [CrossRef] [PubMed]

33. Arcade, A.; Anselin, F.; Rampant, P.F.; Lesage, M.C.; Pâques, L.E.; Prat, D. Application of AFLP, RAPD and ISSR markers to genetic mapping of European and Japanese larch. *Theor. Appl. Genet.* **2000**, *100*, 299–307. [CrossRef]

34. Goulão, L.; Valdiviesso, T.; Santana, C.; Oliveira, C.M. Comparison between phenetic characterization using RAPD and ISSR markers and phenotypic data of cultivated chestnut (*Castanea sativa* Mill.). *Genet. Resour. Crop Evol.* **2001**, *48*, 329–338. [CrossRef]

35. Quézel, P.; Médail, F.; Loisel, R.; Barbero, M. Biodiversity and conservation of forest species in the Mediterranean basin. *Unasylva* **1999**, *50*, 21–28.

36. Kafkas, S. Phylogenetic analysis of the genus *Pistacia* by AFLP markers. *Plant Syst. Evol.* **2006**, *262*, 113–124. [CrossRef]

37. Al-Sousli, M.; Faory, H.; Nakar, M.; Zaid, S.; Al-Safadi, B.; Al-Saghir, M. Genetic relationships among some Pistacia species (Anacardiaceae) in Syria. *Middle East J. Sci. Res.* **2014**, *21*, 1487–1496. [CrossRef]

38. Schlüter, P.M.; Harris, S.A. Analysis of multilocus fingerprinting data sets containing missing data. *Mol. Ecol. Notes* **2006**, *6*, 569–572. [CrossRef]

39. Peakall, R.; Smouse, P.E. GenAIEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research—An update. *Bioinformatics* **2012**, *28*, 2537–2539. [CrossRef]

40. Excoffier, L.; Smouse, P.E.; Quattro, J.M. Analysis of Molecular Variance Inferred from Metric Distances among DNA Haplotypes: Application to Human Mitochondrial DNA Restriction Data. *Genetics* **1992**, *131*, 479–491.

41. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. *Genetics* **2000**, *155*, 945–959. [PubMed]

42. Slatkin, M. Isolation by distance in equilibrium and non-equilibrium populations. *Evolution* **1993**, *47*, 264–279. [CrossRef] [PubMed]
43. Mantel, N. The detection of disease clustering and a generalized regression approach. *Cancer Res.* **1967**, *27*, 209–220. [PubMed]

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).