Next Article in Journal
Light Quality Affected the Growth and Root Organic Carbon and Autotoxin Secretions of Hydroponic Lettuce
Previous Article in Journal
In Silico Identification of MYB and bHLH Families Reveals Candidate Transcription Factors for Secondary Metabolic Pathways in Cannabis sativa L.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 9
Issue 11
10.3390/plants9111541
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Variation in Damaged Populations of Pistacia atlantica Desf.

by
Amina Labdelli
1,2,
Roberto De La Herrán
3,*,
Rami Arafeh
4,
Francesca Resentini
5,6,*,
Livio Trainotti
6,
Youcef Halis
1,
Ahmed Adda
2 and
Othmane Merah
7,8,*
1
Scientific and Technical Research Centre for Arid Areas (CRSTRA), BP 1682 RP, Biskra 07000, Algeria
2
Laboratory of Agro-Biotechnology and Nutrition in Semi-Arid Areas, Ibn Khaldon University, Tiaret 14000, Algeria
3
Departamento de Genética, Facultad de Ciencias, Universidad de Granada, Avda, Fuentenueva s/n, 18071 Granada, Spain
4
Palestine-Korea Biotechnology Center, Palestine Polytechnic University, P.O. Box (198), Hebron, Palestine
5
Department of Biosciences, Università degli Studi di Milano, Via Giovanni Celoria 26, 20133 Milano, Italy
6
Department of Biology, University of Padova, Viale Giuseppe Colombo 3, 35121 Padova, Italy
7
Laboratoire de Chimie Agro-Industrielle (LCA), Université de Toulouse, INRA, INPT, 31030 Toulouse, France
8
Département Génie Biologique, Université Paul Sabatier, IUT A, 32000 Auch, France
*
Authors to whom correspondence should be addressed.
Plants 2020, 9(11), 1541; https://doi.org/10.3390/plants9111541
Submission received: 28 September 2020 / Revised: 9 November 2020 / Accepted: 9 November 2020 / Published: 11 November 2020
(This article belongs to the Section Plant Molecular Biology)
Browse Figures
Versions Notes

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.

Graphical Abstract

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,3,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, in carpentry, and for firewood [2,3,4,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,7,8,9,10,11,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,7,8,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,14,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,17,18,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,22,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,26,27,28,29,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,26,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).
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).
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.

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.

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; Figure 3 and Figure 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].

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 NaC1, 0.5 M EDTA, 2% CTAB, 2% polyvinylpyrrolidone PVP, 0.2% b-mercaptoethanol, and 0.1% NaHSO3).
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-HCI, 0.l 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 1× 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 (5× Green GoTaq® Reaction Buffer), 0.8 µL of MgCl2 (50 mM) (or DMSO dimethyl sulfoxide + buffer 10× PFU Reaction w/20 Mm MgSO4 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 1× TAE buffer and visualized by UV light.

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 GenAlEx 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 log10 transformed values of Nem 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. [Google Scholar]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar]
  6. El Zerey-Belaskri, A. Taxonomic and botanical retrospective review of Pistacia atlantica Desf. (Anacardiaceae). Arab. J. Med. Aromat. Plants 2019, 5, 47–77. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [CrossRef]
  15. Karimi, H.R.; Kafkas, S. Genetic relationships among Pistacia species studied by SAMPL markers. Plant Syst. Evol. 2011, 297, 207–212. [Google Scholar] [CrossRef] [Green Version]
  16. Kafkas, S.; Perl-Treves, R. Morphological and molecular phylogeny of Pistacia species in Turkey. Theor. Appl. Genet. 2001, 102, 908–915. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [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. [Google Scholar]
  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. [Google Scholar] [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. [Google Scholar] [CrossRef] [PubMed]
  32. Weir, B.S.; Cockerham, C.C. Estimating F-Statistics for the Analysis of Population Structure. Evolution 1984, 38, 1358. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar]
  36. Kafkas, S. Phylogenetic analysis of the genus Pistacia by AFLP markers. Plant Syst. Evol. 2006, 262, 113–124. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [CrossRef]
  39. Peakall, R.; Smouse, P.E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research—An update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar]
  41. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar] [PubMed]
  42. Slatkin, M. Isolation by distance in equilibrium and non-equilibrium populations. Evolution 1993, 47, 264–279. [Google Scholar] [CrossRef] [PubMed]
  43. Mantel, N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27, 209–220. [Google Scholar] [PubMed]
Plants 09 01541 g001 550
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.
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.
Plants 09 01541 g001
Plants 09 01541 g002 550
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.
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.
Plants 09 01541 g002
Plants 09 01541 g003 550
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.
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.
Plants 09 01541 g003
Plants 09 01541 g004 550
Figure 4. Pistacia atlantica Desf. atlantica tree and leaves in regions of Batna (A) and Bechar (B) (September 2015). Photos: Labdelli A.
Figure 4. Pistacia atlantica Desf. atlantica tree and leaves in regions of Batna (A) and Bechar (B) (September 2015). Photos: Labdelli A.
Plants 09 01541 g004
Table
Table 1. The number of amplified and polymorphic bands observed at the different locations. ISSR, inter-simple sequence repeats.
Table 1. The number of amplified and polymorphic bands observed at the different locations. ISSR, inter-simple sequence repeats.
Primer BatnaDjelfaTiaretBechar
ISSR3Total number of bands19141913
Number of polymorphic bands137149
% polymorphic bands68.425073.6869.23
ISSR4Total number of bands1414136
Number of polymorphic bands87104
% polymorphic bands57.145076.9266.67
Table
Table 2. Genetic diversity of the four populations of P. atlantica Desf. inferred by the two ISSR markers.
Table 2. Genetic diversity of the four populations of P. atlantica Desf. inferred by the two ISSR markers.
PopulationsNa%PNdNeHeSI
Batna8363.64%1.636 ± 0.0851.504 ± 0.0760.272 ± 0.0395.26 ± 0.29
Djelfa8350%1.500 ± 0.0961.440 ± 0.0860.233 ± 0.0455.34 ± 0.36
Tiaret8375%1.750 ± 0.0781.640 ± 0.0690.344 ± 0.0365.32 ± 0.17
Bechar3368.42%1.684 ± 0.1101.565 ± 0.0710.305 ± 0.0513.6 ± 0.19
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.
Table
Table 3. Geographical locations and numbers of samples collected from each P. atlantica Desf. population analyzed.
Table 3. Geographical locations and numbers of samples collected from each P. atlantica Desf. population analyzed.
Assigned CodeLocationNumber of SamplesAltitude (m)Latitude (N)Longitude
BBatna
Semi-arid		10102735°37′10″6°22′13″ E
TTiaret
Semi-arid		1280835°22′33″02°09′5″ W
DDjelfa
Arid		1163034°02′11″03°40′22″ E
ABechar
Hyperarid		0697932°04′6″02°18′5″ W
Table
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).
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).
Primer NameSequence
5′–3′		Annealing Temperature
(°C)
ISSR 1(CG)9WFrom 50 to 60
ISSR 2(GC)9WFrom 50 to 60
ISSR 3(AG)9B54
ISSR 4(GA)9Y54
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Labdelli, A.; De La Herrán, R.; Arafeh, R.; Resentini, F.; Trainotti, L.; Halis, Y.; Adda, A.; Merah, O. Genetic Variation in Damaged Populations of Pistacia atlantica Desf. Plants 2020, 9, 1541. https://doi.org/10.3390/plants9111541

AMA Style

Labdelli A, De La Herrán R, Arafeh R, Resentini F, Trainotti L, Halis Y, Adda A, Merah O. Genetic Variation in Damaged Populations of Pistacia atlantica Desf. Plants. 2020; 9(11):1541. https://doi.org/10.3390/plants9111541

Chicago/Turabian Style

Labdelli, Amina, Roberto De La Herrán, Rami Arafeh, Francesca Resentini, Livio Trainotti, Youcef Halis, Ahmed Adda, and Othmane Merah. 2020. "Genetic Variation in Damaged Populations of Pistacia atlantica Desf." Plants 9, no. 11: 1541. https://doi.org/10.3390/plants9111541

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop