ارزیابی تنوع ژنتیکی25 ژنوتیپ پسته ایرانی با استفاده از نشانگر مولکولی ISSR

نوع مقاله : مقاله پژوهشی

نویسندگان

چکیده

در مطالعه حاضر، به منظور بررسی تنوع ژنتیکی 25 ژنوتیپ‌ ماده پسته از 17 نشانگر ISSR  استفاده گردید.  با استفاده از نشانگرهای ISSR، در مجموع148 قطعه DNA تکثیر شد که از بین آن‌ها تعداد 123 قطعه (83درصد) چند شکل بودند که متوسط چندشکلی ایجاد شده برای هر آغازگر، 23/7 باند بود. گروه­بندی ژنوتیپ­های مورد مطالعه با استفاده از تجزیه خوشه­ای  به روش UPGMA و بر اساس ماتریس تشابه جاکارد آنها را در 7 گروه مجزا قرار داد. مطابق با نتایج حاصل از ماتریس تشابه بیشترین شباهت بین ژنوتیپ‌های خنجری دامغان و فندقی48 (63%) و کمترین شباهت بین ژنوتیپ‌های ابراهیمی و شستی (2%) مشاهده شد  که نشان دهنده سطح بالای چندشکلی در بین ژنوتیپ‌های مورد مطالعه است. نتایج این تحقیق نشان داد که نشانگر ISSR برای شناسایی ژنوتیپ‌های ماده پسته و تهیه­ی شناسنامه مولکولی کامل برای توده­های پسته ایرانی مناسب است.

کلیدواژه‌ها


عنوان مقاله [English]

Evaluation of genetic diversity in 25 Iranian pistachio genotypes using ISSR markers

نویسندگان [English]

  • Ramin Arjmand Qahestani
  • Iraj Tavassolian
  • Qasem Mohammadi Nezhad
چکیده [English]

In this study 17 ISSR molecular markers were employed to investigate the genetic diversity among 25 female pistachio genotypes. The markers produced  148 bands in total, of which 123 (83%) was polymorphic and the average of polymorphism for each primer was 7.23 band per primer. Cluster analyses of genotypes were estimated based on Jaccard's similarity algorithm and UPGMA clustering method and placed the genotypes in seven groups. The results obtained by similarity matrix revealed the highest similarity between cultivars "Khanjari Damghan" and "Fandoghi 48" (63%) and the lowest similarity between cultivars "Ebrahimi" and "Shasti" (2%) which is due to high amount of polymorphism among the genotypes. This experiment shows that ISSR markers are reliable for assessment of pistachio female genotypes and genetic diversity study. Moreover, these markers can be used to provide molecular genetics identity for Iranian pistachio populations.

کلیدواژه‌ها [English]

  • Pistachio
  • genetic diversity
  • Jaccard similarity coefficient
  • Cluster analysis
  • ISSR

 

ارزیابی تنوع ژنتیکی25 ژنوتیپ پسته ایرانی با استفاده از نشانگر مولکولی  ISSR

 

رامین ارجمند قهستانی1، ایرج توسلیان*2، قاسم محمدی نژاد2

 

1دانش آموخته کارشناسی ارشد علوم باغبانی، دانشگاه شهید باهنر کرمان، کرمان، ایران.

2پژوهشکده باغبانی، دانشگاه شهید باهنر کرمان، کرمان، ایران.

 

تاریخ دریافت: 19/011/1391، تاریخ پذیرش: 25/11/1393

 

چکیده

در مطالعه حاضر، به منظور بررسی تنوع ژنتیکی 25 ژنوتیپ‌ ماده پسته از 17 نشانگر ISSR  استفاده گردید.  با استفاده از نشانگرهای ISSR، در مجموع148 قطعه DNA تکثیر شد که از بین آن‌ها تعداد 123 قطعه (83درصد) چند شکل بودند که متوسط چندشکلی ایجاد شده برای هر آغازگر، 23/7 باند بود. گروه­بندی ژنوتیپ­های مورد مطالعه با استفاده از تجزیه خوشه­ای  به روش UPGMA و بر اساس ماتریس تشابه جاکارد آنها را در 7 گروه مجزا قرار داد. مطابق با نتایج حاصل از ماتریس تشابه بیشترین شباهت بین ژنوتیپ‌های خنجری دامغان و فندقی48 (63%) و کمترین شباهت بین ژنوتیپ‌های ابراهیمی و شستی (2%) مشاهده شد  که نشان دهنده سطح بالای چندشکلی در بین ژنوتیپ‌های مورد مطالعه است. نتایج این تحقیق نشان داد که نشانگر ISSR برای شناسایی ژنوتیپ‌های ماده پسته و تهیه­ی شناسنامه مولکولی کامل برای توده­های پسته ایرانی مناسب است.

کلید واژه: پسته، تنوع ژنتیکی، ضریب تشابه جاکارد، تجزیه خوشه­ای، ISSR.

 



مقدمه

پسته به عنوان یکی از مهمترین محصولات باغی و سومین کالای صادراتی کشور به لحاظ ارزآوری، از اهمیت ویژه ای برخوردار است. گسترش روزافزون کشت پسته در دیگر مناطق جهان و میانگین تولید بالاتر در واحد سطح در مقایسه با ایران، تهدید جدی برای صادرات این محصول در کشور می باشد. از طرف دیگر برای احداث باغ های یکدست و با عملکرد و کیفیت بالای هر محصول از طریق بهنژادی، در مرحله اول نیازمند بررسی صفات مهم مورفولوژیکی وشناسایی ژنوتیپ‌های برتر هستیم (Al-saghir (and Porter, 2006. جنس پسته متعلق به خانواده پسته‌سانان[1] و راسته Sapindales است (Zohary, 1952). پسته به عنوان گیاهی که نقش مهمی در تغذیه و اقتصاد کشورهایی مانند ایران، ترکیه و سوریه دارد دارای تنوع ژنتیکی بالایی می‌باشد (Ozden-Tokatli et al., 2010). این جنس دارای بیش از 13 گونه است که به صورت درخت یا درختچه‌های خزان پذیر یا همیشه‌سبز و دو‎پایه می‌باشند. گل آذین در گونه‌های این جنس به‌صورت خوشه یا پانیکول و میوه از نوع شفت می‎باشد. در بین گونه‌های پسته فقط گونه ورا[2] دارای خشک‌میوه‌ی خندان است و ارزش تجاری دارد (Zohary, 1952). سایر گونه‌ها و زیرگونه‌ها خشک میوه‌های کوچکتری تولید می‌کنند و به‌طور عمده به عنوان پایه و یا استخراج روغن، جنگل‌کاری و تولید الوار استفاده می‌شوند (Barghchi & Alderson, 1989). همه گونه‌های جنس پسته دوپایه هستند و با باد گرده افشانی می‌شوند. مجموعه‌ای از درختان آتلانتیکا (P. atlantica) در کوه‌های یانت[3] ترکیه یافت شده است که به صورت یک‌پایه بوده است ( Kafkas et al.,2000).اعتقاد بر این است که قدیمی‌ترین گونه موجود در جنس پسته، گونه ورا می‌باشد و احتمالاً دیگر گونه‌ها از آن مشتق شده‌اند (Zohary, 1952). تا کنون جامع‌ترین رده‌بندی جنس پسته توسط زهری گزارش شده است که جنس پسته را به 4 بخش و 11 گونه بر اساس خصوصیات برگ و مورفولوژی میوه تقسیم کرده است (Zohary, 1952). اولین طبقه بندی گونه‌های پسته در سطح مولکولی با استفاده ازDNA کلروپلاستی انجام شد و گزارش شد که جنس پسته در دو بخش تربینتوس Terebinthus)) و لنتیسکوس (Lentiscus) به ترتیب خزان‌پذیر و همیشه‌سبز تقسیم می‌شود (Parfitt & Badenes, 1997). ایران منشا سه گونه ورا، خنجوک و موتیکا (P. vera, P. mutica, P. khinjuk) از جنس پسته می‌باشد.

به این منظور طی دهه‌های گذشته عمده ژنوتیپ‌های پسته کشورهای ایران، آمریکا، ایتالیا، سوریه و ... شناسایی و از لحاظ مورفولوژیکی و مولکولی مورد بررسی قرار گرفته‌اند. در مطالعات اولیه مربوط به تنوع ژنتیکی پسته، نشانگرهای مورفولوژیکی مانند شکل میوه و برگ مورد استفاده محققین قرار گرفتند (Zohary, 1952). در مراحل بعدی مطالعات آیزوزایمی متعددی با هدف تمایز بین گونه‌ها و ارقام پسته انجام شد (Dyszel & Petitt, 1990) که در تمامی موارد، ناکافی بودن چندشکلی آیزوزایم‌ها در بین ارقام نزدیک به هم، قابلیت آن‌ها را در تشخیص  و بررسی دقیق روابط خویشاوندی مورد تردید قرار داد. بنابراین به تدریج مطالعات مولکولی پسته به سمت استفاده از نشانگرهای مبتنی بر DNA متمایل گشت ( Golan-Goldhirsh et al., 2004). تحقیقات مختلفی با روش­های متفاوتی چون RAPD و AFLP (Katsiotis et al., 2003; Golan-Goldhrish et al., 2004)، روش SSRs (Vandramin et al., 2009) و ایزوآنزیم­ها (Baron et al., 1996) به منظور بررسی تفاوت­های ژنتیکی بین گونه­های مختلف پسته انجام شده است. تنوع ژنتیکی 30 ژنوتیپ پسته شامل 24 رقم پسته ورا، واریته سرخس، گونه‌های خنجوک، بنه، بنه باغی، اینتگریما و پالستینا را با سه سیستم ایزوزایمی (استراز، پراکسیداز و مالات دهیدروژناز) مورد بررسی قرار گرفت ( .(Alami et al., 2003 براساس پژوهش آنها تجزیه کلاستر ژنوتیپ‌ها را به هشت گروه اصلی و 20 گروه فرعی تقسیم نمود به طوری‌که پنج گروه اول، به ارقام پسته اختصاص یافت و سه گروه پایانی گونه‌ها را شامل شد. نتایج حاصل از مطالعه ایشان نشان داد که خنجوک بیشترین قرابت را با ارقام رایج پسته دارد. بررسی روابط درون گونه‎ای تعدادی از گونه‎های جنس پسته با استفاده از نشانگر RAPD نشان داد که گونه‎های مورد بررسی تشکیل دو گروه مجزا می‌دهد که گروه اول شامل گونه‌های ورا، آتلانتیکا، خنجوک، یوریکارپا و انتیگریما[4] و گروه دوم شامل گونه‌های تربینتوس، مکزیکانا[5]، پالستینا[6] و لنتیسکوس[7] می‌باشد (Kafkas & Perl-Treves, 2002). در پژوهش ایشان، 20 آغازگر مورد استفاده توانست 228 باند چندشکل در سطح درون گونه‌ای ایجاد نماید.

از نشانگرهای RAPD وAFLP برای مطالعه روابط ژنتیکی گونه‌های پسته موجود در یونان استفاده شد و بر اساس داده‌های مولکولی بدست آمده، گونه‌های مورد بررسی به دو دسته همیشه سبز و خزان پذیر تفکیک شدند (Katsiotis et al., 2003).

الصغیر و پرتر مطالعه فیلوژنی 4 گونه پسته (خنجوک، لنتیسکوس، تربینتوس و ورا)  بوسیله نشانگر RAPD راگزارش کردند. در این گزارش23 آغازگر مورد استفاده، 248 باند تولید کرد که 139 عدد آن در سطح درون گونه‌ای چند شکلی نشان داد. تجزیه کلاستر داده‌ها نشان داد که خنجوک نزدیک‌ترینگونه به ورا می‌باشد (Al-Saghir, & Porter, 2006).

کفکاس روابط فیلوژنی بین گونه‌های پسته موجود در ترکیه را بوسیله نشانگر AFLP بررسی نمود و گزارش کرد که ژنوتیپ‎های خنجوک در ترکیه رابطه ژنتیکی نزدیک‌تری به گونه‌های یوریکارپا و آتلانتیکا دارد که دلیل آن به خاطر تعیین هویت اشتباه نمونه‌های خنجوک می‌باشد (Kafkas, 2006).

میرزایی و همکاران با استفاده از نشانگر مولکولی RAPD روابط فیلوژنی گونه‌های وحشی و برخی ارقام پسته موجود در ایران را بررسی نمودند. از مجموع 77 آغازگر مورد استفاده در پژوهش ایشان، 15 آغازگر تولید الگوهای چند شکل نمودند که مجموعاً 57 قطعه چند شکل تولید شد (Mirzaei et al., 2002). تجزیه کلاستر بر اساس ضریب تشابه جاکارد و به روش [8]UPGMA ژنوتیپ‌های مورد مطالعه را در سه گروه گونه ورا، فرم وحشی ورا و گونه وحشی بنه قرار داد. نتایج حاصل نشان داد، ارقام پسته مورد بررسی سطح تشابه بالایی دارند و واریته سرخس بیشترین رابطه ژنتیکی را با ارقام اهلی پسته دارد، در حالی‎که خنجوک و بنه کم‎ترین رابطه را نشان می‌دهند. آنها گزارش دادند که تفاوت‌های موجود در ارقام رایج باغی به علت موتاسیون در اندام‌های رویشی و زایشی و همچنین به علت تکثیر رویشی می‌باشد.

بررسی تنوع ژنتیکی ارقام و گونه‎های وحشی پسته موجود در ایران به همراه 22 رقم خارجی با استفاده از نشانگر SSR انجام شد. نتایج حاصل از مطالعه نشان داد که تنوع ژنتیکی در زیر گونه کردیکا[9] پایین‎تر از گونه ورا و خنجوک است. تجزیه به مولفه های اصلی[10] به‎طور واضحی ارقام و گونه‌های ایرانی را از غیر ایرانی تفکیک نمود و نتایج حاصل از تجزیه واریانس مولکولی نشان داد که، جریان ژنی بین ارقام ایرانی و جمعیت‌های وحشی پسته ورا وجود دارد (    et al., 2010 Pazouki).

بررسی تعیین جنسیت  چهار رقم پسته ماده (اکبری، احمدآقایی، فندقی، کله قوچی) و یک ژنوتیپ پسته نر با استفاده از نه آغازگر ISSR (تکثیر بین ردیف‌های تکراری ساده) نشان داد که دو آغازگر (AC)8CG و (AC)8TA قادر به تعیین هویت گیاهان ماده و نر بوسیله تولید باندهای DNA وابسته به جنسیت در گیاهان ماده شدند و گزارش شد که نشانگر ISSR جهت تعیین جنسیت گیاهان پسته در مرحله نونهالی مناسب است (Ehsanpour et al., 2008). با توجه به موارد فوق و این که استفاده از نشانگرهای ISSR نیازی به اطلاعات توالی ژنوم ندارد و منجر به ایجاد الگوهای چندجایگاهی و بسیار چندشکل می­شود (Askari et al., 2011; Ghasemi et al., 2010; Zamani et al., 2011; Zamani et al., 2015) هدف از اجرای این پژوهش ارزیابی تنوع ژنتیکی و شناسایی چند شکلی ژنتیکی بین ارقام پسته با استفاده از نشانگر مولکولی ISSR می­باشد.


مواد و روش ها

در این مطالعه نمونه برگ 25 ژنوتیپ ماده پسته گونه ورا، از ایستگاه شماره دو مؤسسه تحقیقات پسته کشور واقع در رفسنجان جمع‌آوری گردید و بر روی یخ به آزمایشگاه اصلاح نباتات دانشگاه شهید باهنر کرمان انتقال داده شد. از هر ژنوتیپ، برگ‌های جوان و سالم به طور مجزا در پاکت‌های نایلونی زیپ‌دار قرار داده شد، سپس نام آن بر روی پاکت درج گردید و بلافاصله نمونه‌ها به فریزر 20- درجه سانتی‌گراد منتقل شد (جدول 1).

DNA ژنومی از بافت برگ تازه، به روش CTAB[11]  مطابق مقاله دویل و دویل با کمی تغییرات استخراج گردید (Doyle & Doyle, 1987). ترکیب بافر استخراج شامل Tris-HCl، NaCl، EDTA، CTAB، Na2S2O5، PVP، β-mercaptoethanol طبق جدول 2 تهیه شده و در ادامه کیفیت و کمیت DNA استخراج شده با استفاده از روش های اسپکتوفتومتری و الکتروفورز با  ژل آگارز 8/0% تعیین گردید.

واکنش تکثیر قطعات DNA با استفاده از 17 آغازگر ISSR به طول 18 تا 24 نوکلئوتیدی (جدول3) انجام گرفت. هر واکنش در حجم 15 میکرولیتر و با استفاده ازمواد جدول 4 صورت گرفت.

پس از اتمام واکنش PCR به هر میکروتیوب، 5 میکرولیتر بافر بارگذاری اضافه شد و مقدار 10 میکرولیتر از مخلوط حاصله به عنوان پیش بار استفاده شد. جهت پیش لود کردن، از غلظت 8/0 درصد آگارز و جهت بارگذاری نهایی، از غلظت 5/1 درصد آگارز و ولتاژ ثابت 100 ولت به مدت 3 ساعت استفاده شد. ژل 5/1 درصد آگارز، در داخل قالب ریخته و پس از سرد شدن کامل ژل، شانه را با دقت از ژل جدا نموده و DNA مربوط به‌ همراه بافر لودگذاری در چاهک‌ها لود شد و محصولات تکثیر توسط دستگاه الکتروفورز تفکیک شد. سپس ژل­ها در محلول اتیدیوم بروماید رنگ‌آمیزی شد. پس از شستشو ژل با آب مقطر به دستگاه ژل داک جهت مشاهده باندهای تکثیر شده زیر نور فرابنفش منتقل شد و از ژل­ها عکس برداری شد.

پس از انجام آزمایشات ISSR، الگوهای باندی، بر اساس حضور و عدم حضور باند، به صورت یک و صفر امتیازدهی شدند. بعد از تشکیل ماتریس صفر و یک ، ماتریس تشابه ژنوتیپ‌ها با استفاده از ضریب تشابه جاکارد (Jaccard, 1908) و دندروگرام براساس ماتریس تشابه و به روش UPGMA بدست آمد. تجزیه داده‌ها با استفاده از نرم‌افزار.NTSYS Ver 2.02 انجام گرفت. به منظور انجام تجزیه به مختصات اصلی بر اساس نشانگرهای ISSR  که هدف از این تجزیه (PCoA)[12]، یافتن ترکیباتی از P متغیر جهت ایجاد شاخص­های مستقل به نام مؤلفه­های اصلی (PC) می­باشد. نیز از نرم افزارNTSYS Ver 2.02  استفاده گردید.


 

 

 

جدول 1- لیست ژنوتیپ های ماده پسته مورد مطالعه با نشانگرهای ISSR.

Table 1- List of the female genotypes were tested in this study by ISSR markers.

نام رقم

Cultivar

شماره

Number

نام رقم

Cultivar

شماره

Number

نام رقم

Cultivar

شماره

Number

 (Italian)ایتالیایی

19

سیریزی (Cirizi)

10

هراتی (Harati)

1

بادامی راور (Badami Ravar)

20

سبزپسته نوق  (Sabzpesteh Nogh)

 

11

خنجری دامغان (Khanjari Damghan)

2

 (Shasti)شستی

21

ابراهیم آبادی (Ebrahim Abadi)

12

ممتاز (Momtaz)

3

فندقی زودرس (Fandoghi Zodras)

22

ممتاز تاج آبادی (Momtaz Taj Abadi)

13

ابراهیمی (Ebrahimi)

4

فندقی ریز (Fandoghi Riz)

23

 (Jandaghi) جندقی

 

14

فندقی48 (Fandoghi48)

5

سفیدپسته نوق (Sefidpesteh Nogh)

24

خنجری راور (Khanjari Ravar)

 

15

شاهپسند (Shah Pasand)

6

 (Amiri)امیری

25

قزوینی (Ghazvini)

 

 

16

بادامی زرند (Badami Zarand)

7

   

اوحدی (Ohadi)

 

17

موسی آبادی (Mosa Abadi)

8

 

 

بادامی نیش کلاغی (Badami Nishklaghi)

18

حسن زاده (Hassan Zadeh)

9

 


جدول 2- آغازگرهای ISSR مورد استفاده در این تحقیق.

Table 2- ISSR Primers used in the present study.

Length (Mer)

Tm ºC

Oligo Sequence

Oligo Name

No

18 Mer

551

GAGAGAGAGAGAGAGARC

(GA)8-RC

1

18  Mer

54

AGAGAGAGAGAGAGAGYT

(AG)8YT

2

18 Mer

54

ACACACACACACACACYG

(AC)8YG

3

18 Mer

51

TGTGTGTGTGTGTGTGRT

(TG)8RT

4

18 Mer

51

CACACACACACACACART

(CA)8RT

5

18 Mer

54

CTCTCTCTCTCTCTCTRG

(CT)8RG

6

24  Mer

65

BDBTCCTCCTCCTCCTCCTCCTCC

BDB(TCC)7

7

19  Mer

56

CCACTCTCTCTCTCTCTCT

CCA(CT)8

8

18  Mer

56

AGAGAGAGAGAGAGAGSC

(AG)8SC

9

20 Mer

55

GACAGACAGACAGACAGACA

(GACA)5

10

18 Mer

47

GAAGAAGAAGAAGAAGAA

(GAA)6

11

18 Mer

54

CACACACACACACACAGT

(CA)8-GT

12

18Mer

55

GTGTGTGTGTGTGTGTYC

(GT)8-YC

13

17Mer

52

TGTGTGTGTGTGTGTGG

(TG)8G

14

18  Mer

55

ACACACACACACACACYG

(AC)8YT

15

17  Mer

50

AGAGAGAGAGAGAGAGT

(AG)8T

16

18  Mer

52

CACACACACACACACART

(CA)8RT

17

 

جدول 3- مواد لازم و مقادیر مربوطه برای تهیه محلول واکنش ISSR-PCR.

Table 3- Reagents and their concentrations  for ISSR-PCR reactions.

مواد

Materials

غلظت پایه

Basic level

غلظت نهایی

Final level

یک واکنش

One reaction

24واکنش×1/1*

24 reactions

آب دوبار تقطیر استریل

(Double Distillated water)

_

 

_

 

 8.4µl

 

221.76 µl

 

 

بافر PCR

(PCR buffer)

x 10

 

1x

 

1.5 µl

 

39.6 µl

 

MgCl2

mM 50

1.75mM

 0.6 µl

15.84 µl

dNTPs

mM 10

mM 10

 0.3 µl

7.92 µl

آغازگر Primers

10µM

µM 0.2

 2 µl

52.8 µl

DNA الگو

(Template DNA)

ng/µl 30

0 ng/µl 30

 2 µl

-

آنزیم Taq  (Taq DNA Polymerase)

5 Unit/µl

Unit/µl 1

 0.2 µl

5.28 µl

 

جمع کل (Total volume)

 

 

15 µl

 

* ضریب خطای نمونه برداری  ( Sampling error coefficient)

 

جدول 4- برنامه زمانی و دمای مورد نیاز برای انجام مراحل واکنش ISSR- PCR.

Table 4- PCR protocol for the amplification of ISSR-PCR.

مرحله

Step

مرحله انجام شده

Phase

تعدادسیکل

Number of cycle

دما

Temperature

زمان

Time (minute)

1

شروع باز شدن DNA کروموزومی

(Initial Denaturation)

1

ºC 94

4

 

تک رشته‌ای شدن DNA

(Denaturation)

 

ºC 94

2

2

اتصال آغازگر (Annealing)

40

سیکل (40 cycles)

متفاوت است (Variable)

1

 

بسط آغازگر (Extension)

 

ºC 72

2

3

تکمیل بسط (final extension)

1

ºC 72

7

 

 

 

ºC 4

1

 


نتایج و بحث

برای بررسی چندشکلی DNA بین ژنوتیپ‌های پسته مورد آزمایش از 17 آغازگر ISSR استفاده گردید. اندازه قطعات در تمام آغازگرها بین 24- 18 نوکلئوتید بود و فقط باندهای با وضوح بالا مورد تجزیه و تحلیل قرار گرفت. 17 آغازگر مذکور، مجموعاً 148 قطعه DNA را تکثیر کردند که از بین آن‌ها تعداد 123 قطعه چند شکل و 25 قطعه یک شکل بودند. متوسط باند تولید شده برای هر آغازگر، 70/8 باند و متوسط چندشکلی ایجاد شده برای هر آغازگر، 23/7 باند بود (جدول 6). در پژوهشی مطالعه تنوع ژنتیکی برخی از ژنوتیپ های پسته با استفاده از آغازگرهای ISSR گزارش گردید که از 114 باند تولید شده، 73 باند (64درصد) چندشکلی نشان دادند (Taghizadeh et al, 2010). دلیل چندشکلی کمتر ایجاد شده در پژوهش آن‌ها احتمالا، تعداد ژنوتیپ‌های کمتر (19) و همچنین متفاوت بودن آغازگرهای به‌کار رفته می‌باشد. تعداد باند تکثیر شده با آغازگرهای مختلف، متفاوت بود، به‌طوری‌که بیشترین تعداد باندهای تولید شده متعلق به آغازگر (AC)8YG 15 باند و کمترین باند متعلق به آغازگر (TG)8G چهار باند بود. بیشترین درصد چندشکلی مربوط به آغازگرهای (CT)8RC، (AG)8YT، (AC)8YG، BDB(TCC)7، (GACA)5، (CA)8-GT، (AC)8YT، (100%) و کمترین درصد چندشکلی مربوط به آغازگر (AG)8T (20%) بود (جدول5).

الگوهای باندی نشانگرهای چندشکل ISSR بر اساس وجود یا عدم وجود باند به‌صورت اعداد یک و صفر، برای بدست آوردن ماتریس تشابه ژنوتیپ‌ها بر اساس ضریب تشابه جاکارد استفاده شد. ضریب تشابه جاکارد، بیشترین استفاده را برای بررسی میزان تنوع ژنتیکی بدست آمده از نشانگرهای غالب است، زیرا شباهت‌ها را براساس وجود و یا عدم وجود باند گروه‌بندی می‌کند (Naghavi et al, 2009). ماتریس تشابه میزان تشابه دو به دوی بین ژنوتیپ‌ها را نشان می‌دهد. بنابراین هر چه دو ژنوتیپ از نظر نشانگرهای مختلف الگوی شبیه‌تری  با هم‌ داشته باشند، تشابه ژنتیکی آنها بیشتر خواهد بود و بالعکس. 

 

 

جدول 5- لیست آغازگرهای مورد استفاده و درصد چندشکلی حاصل از آن‌ها.

Table 5- List of the ISSR Primers used and their polymorphism percentage.

ردیف

Number

آغازگر

Primers name

تعداد کل قطعات تکثیر شده (a)

Number of amplified fragments (a)

تعداد قطعات چندشکل(b)

Number of polymorphic fragments (b)

درصد چندشکلی

(b/a)

Percentage of polymorphism

1

(GA)8-RC

4

4

100

2

(CT)8RC

8

8

100

3

(AG)8YT

8

8

100

4

(AC)8YG

5

5

100

5

(TG)8RT

7

1

100

6

(CA)8RT

8

7

87.5

7

BDB(TCC)7

8

8

100

8

CCA(CT)8

8

7

87.5

9

(AG)8SC

7

6

85.71

10

(GACA)5

6

6

100

11

(GAA)6

4

3

75

12

(CA)8-GT

7

7

100

13

(GT)8-YC

3

3

100

14

(TG)8G

4

3

75

15

(AC)8YT

7

7

100

16

(AG)8T

5

4

80

17

(CA)8RT

6

5

83.3

کل (Total)

-

115

107

-

میانگین(Means)

-

5

4.65

93.20

 

   

   

 

براساس نتایج حاصل از ماتریس تشابه، بیشترین شباهت (63%) بین ژنوتیپ‌های خنجری دامغان و فندقی48 و کمترین شباهت (2%) بین ژنوتیپ‌های ابراهیمی با شستی و ممتازتاج آبادی و خنجری دامغان و فندقی48 وجود داشت که نشان دهنده سطح بالای چندشکلی در بین ژنوتیپ‌های مورد مطالعه است. همچنین شباهت بین ژنوتیپ‌های خنجری راور و اوحدی (62%)، بادامی راور و خنجری راور (60%)، بدست آمد. در بررسی 19 ژنوتیپ پسته اهلی با استفاده از 20 آغازگر ISSR، 64 درصد چندشکلی گزارش گردید (Taghizadeh et al, 2010). در مطالعات مشابهی میرزایی و همکاران (2002) ضمن مطالعه تنوع ژنتیکی گونه‌های وحشی و ارقام پسته موجود در ایران با استفاده از نشانگر RAPD، درصد چندشکلی ایجاد شده را 73/74 درصد گزارش کردند (Mirzaei et al., 2002).  در پژوهش دیگری کمترین فاصله بین ارقام کله­قوچی و احمدآقایی با میزان تشابه 83/0 و بیشترین فاصله ژنتیکی بین ارقام کله­قوچی، حسین­خانی و ابراهیم­آبادی با میزان تشابه 54/0 مشاهده شد (Taghizadeh et al, 2012).

 

تجزیه خوشه ای حاصل از داده‌های مولکولی

دندروگرام بر اساس ماتریس تشابه و به روش UPGMA بدست آمد. ضریب همبستگی کوفنتیک به‌دست آمده از دندروگرام و ماتریس تشابه 94/0 بود که نشان‌دهنده برازش خیلی خوب بین ماتریس تشابه و دندروگرام است.  در مطالعه باقی­زاده و همکاران (2010) دندروگرام حاصل از داده­های بدست آمده بر اساس ترکیب هر سه مارکر RAPD، ISSR و SSR با دندروگرام حاصل از تجزیه خوشه­ای جداگانه داده­های هر مارکر به صورت جداگانه مطابقت داشت.

در فاصله تشابه 36/0 ژنوتیپ‌ها در هفت گروه مجزا قرار گرفتند (شکل 1)، که به شرح زیر می‌باشند.

گروه اول: به دو زیر گروه تقسیم شد. زیرگروه اول شامل ژنوتیپ‌های هراتی، خنجری دامغان، فندقی48، موسی آبادی و زیرگروه دوم شامل ژنوتیپ های فندقی زودرس و سفیدپسته نوق می‌باشد. در این گروه بیشترین تشابه بین ژنوتیپ‌های خنجری دامغان، فندقی48 (63%) و کمترین تشابه بین ژنوتیپ‌های هراتی و سفیدپسته نوق (24%) مشاهده شد.

گروه دوم: شامل ژنوتیپ های ممتاز، ابراهیم آبادی، ممتازتاج آبادی، که بیشترین تشابه (51/0) بین ژنوتیپ ابراهیم آبادی و ممتازتاج آبادی و کمترین تشابه (42/0) بین ارقام ممتاز و ابراهیم­آبادی مشاهده شد. نتایج فوق با نتایج میرزایی و همکاران مبنی بر قرار گرفتن ارقام ممتاز تاج آبادی و ممتاز در یک گروه مطابقت دارد (Mirzaei et al., 2002).

گروه سوم: شامل ژنوتیپ های شاهپسند، بادامی زرند، سیریزی، سبزپسته نوق، قزوینی، که بیشترین تشابه (48/0) بین ژنوتیپ سبزپسته نوق و قزوینی و کمترین تشابه (34/0) بین ارقام شاهپسند و قزوینی مشاهده شد. تشابه ژنوتیپ های سبزپسته نوق و بادامی زرند و تشابه بادامی زرند و سیریزی طی گزارشی منتشر گردیده است (Mirzaei et al., 2002 ; Baghizadeh et al., 2010). گروه چهارم: در برگیرنده ژنوتیپ‌های خنجری راور، اوحدی، شستی، ایتالیایی، بادامی راور، بادامی نیش کلاغی، امیری می باشد که بادامی نیش کلاغی و امیری جدا از بقیه در زیر گروه دیگری قرار گرفته اند. در این گروه بیشترین تشابه بین ژنوتیپ‌های خنجری راور و اوحدی (62%) و کمترین تشابه بین ژنوتیپ‌های اوحدی و بادامی نیش کلاغی (32%) مشاهده شد. شباهت بالای ژنوتیپ های خنجری راور و اوحدی توسط نوروزی و همکاران، نیز گزارش شده است (Noroozi et al., 2009). در پژوهشی تشابه ژنوتیپ اوحدی و بادامی نیش کلاغی گزارش شد (Baghizadeh et al., 2010). همچنین ارقام اوحدی و ایتالیایی در مطالعه­ای در یک گروه واقع شدند (Alami et al., 2003). ژنوتیپ های ایتالیایی و بادامی نیش کلاغی طی بررسی با استفاده از نشانگر مولکولی SSR، در یک گروه قرار گرفتند (Arabnezhad et al., 2008).

گروه پنجم: این گروه شامل ژنوتیپ حسن زاده و جندقی می باشد.

گروه ششم: ژنوتیپ فندقی ریز در این گروه جداگانه قرار گرفت.

 گروه هفتم: در این گروه ژنوتیپ ابراهیمی به تنهایی قرار گرفت که در گزارشی آمده است رقم ابراهیمی از سایر ارقام پسته مجزا شده و در یک گروه جدا، واقع می‌شود (Baghizadeh et al., 2010). بای‌پلات و تری­پلات (شکل 2و3) براساس دو و سه مولفه اصلی اولیه ترسیم شده­اند و پراکندگی ارقام را نشان می­دهند.

 

 

 شکل 1- دندروگرام حاصل از داده­های نشانگرهای ISSR مربوط به 25 ژنوتیپ ماده پسته به روش UPGMA.

Figure 1- UPGMA dendrogram revealed by 25 female genotypes of pistachio based on ISSR markers.

 

 

شکل 2- نتایج تجزیه ‌پلات دو بعدی ژنوتیپ‌های پسته مورد مطالعه با استفاده از نشانگر مولکولی ISSR.

Figure 2- Results of two dimentional plot showing relationship among 25 female genotypes of pistachio using ISSR markers.

 

شکل 3- نتایج تجزیه ‌پلات سه بعدی ژنوتیپ‌های پسته مورد مطالعه با استفاده از نشانگر مولکولی ISSR.

Figure 3- Results of three dimentional plot showing relationship among 25 female genotypes of pistachio using ISSR markers.

 



نتیجه­گیری

در مطالعه حاضر درصد چندشکلی کل به دست آمده 10/83 درصد بود. بیشترین تعداد باندهای تولید شده متعلق به آغازگر (AC)8YG 15 باند و کمترین باند متعلق به آغازگر (TG)8G چهار باند بود. بیشترین شباهت (63%) بین ژنوتیپ‌های خنجری دامغان و فندقی48 و کمترین شباهت (2%) بین ژنوتیپ‌های ابراهیمی با شستی و ممتازتاج آبادی و خنجری دامغان و فندقی48 وجود داشت. با توجه به مقایسات انجام شده می‌توان گفت که نشانگر ISSR به خصوص نشانگرهای مورد استفادهی این پژوهش که بالاترین درصد چند شکلی را نشان دادند((CT)8RC، (AG)8YT، (AC)8YG، BDB(TCC)7، (GACA)5، (CA)8-GT، (AC)8YT)، برای شناسایی ژنوتیپ‌های پسته مناسب می‌باشد.


 

منابع

Ahmadiafzadi M, Tabatabaei S B, Mohammadi SA, Tajabadipour A (2007).  Comparison of genetic diversity in species and cultivars of Pistachio (Pistacia sp.) based on Amplified Fragment Length Polymorphism (AFLP) markers. Iranian Journal of Agricultural Biotechnology 3: 147-152.

Alami E, Taeb M, Lotfi E, Sadeghian Motahar I (2003). Study of polymorphic esterase isosymes, peroxidase and mallat dehydrogenase in Iranian pistachio cultivars and species. Agricultural and environmental science and technology 7: 107 – 113.

Al-saghir MG, Porter DM (2006). Random Amplified Polymorphic DNA (RAPD) markers of Pistacia species (Anacardiaceae). Asian Journal of Plant Sciences 5: 1002-1006.

Arabnezhad HM, Bahar A, Tajabadi Pour (2008). Assessment of genetic diversity among Iranian pistachios using microsatellites isolated from Pistacia khinjuk. Journal of Science and Technology of Agriculture and Natural Resources 45: 218-228.

Askari N, Mohammad Abadi MR, Baghizadeh A (2011). ISSR markers for assessing DNA polymorphism and genetic characterization of cattle, goat and sheep populations. Iranian Journal of Biotechnology 9: 222–229.

Baghizadeh A, Noroozi Sh, Jalili Javaran M (2010). Study on genetic diversity of some Iranian pistachio (Pistacia vera L.) cultivars using Random Amplified Polymorphic DNA (RAPD), Inter Sequence Repeat (ISSR) and Simple Sequence Repeat (SSR) markers: A comparative study. African Journal of Biotechnology 9: 7632-7640.

Barghchi M, Alderson PG (1989). Pistachio (Pistacia vera L.). In: Biotechnology in Agriculture and Forestry. Bajaj, Y.P.S. (ed.). Springer Verlag, Berlin 5: 68-98.

Baron, E, Marco L.D, Marra F.P,  Sidari M (1996). Isozymes and canonical discriminant analysis to identify pistachio (Pistacia vera L.) germplasm. Horticulture Science 31, 134-138.

Doyle JJ, Doyle JL (1987). A rapid isolation procedure for small quantities of fresh leaf tissue. Phytochem Bulletin 19:11–15.

Dyszel SM, Petitt BC (1990). Determination of the country of origin of pistachio nut by DSC and HPLC. Journal of the American Oil Chemists' Society 67:947-951.

Ehsanpour AA, Tavassoli M, Arab L (2008). Sex determination of Pistacia vera L. using ISSR markers. Malaysian Applied Biology Journal 37: 25-28.

Ghasemi M, Baghizadeh A, Mohammadabadi MR (2010). Determination ofgenetic polymorphism in Kerman Holstein and Jersey cattle population using ISSR markers. Australian Journal of Basic Applied Science 4: 5758–5760.

Golan-Goldhirsh A, Barazani O, Wang Z S , Khadka DK, Saunders JA, Kostiukovsky V, Rowland LJ (2004). Genetic relationships among Mediterranean Pistacia species evaluated by RAPD and AFLP markers. Plant Systematics and Evolution 246: 9–18.

Jaccard P (1908). Nouvellesrecherchessur la distribution.florale. Bulletin de la Société vaudoise des sciences naturelles. 44: 223-270.

Kafkas S, Perl-Treves RN, Kaska (2000). Unusual Pistacia atlantica Desf. monoecious sex types in the Yunt Mountains of the Manisa province of Turkey. Israel  Journal  of Plant Sciences 48: 277-280.

Kafkas S, Ebru K, Perl-Treves R (2002). Morphological diversity and germplasm survey of three wild Pistacia species in Turkey. Genetic Resources and Crop Evolution 49: 261-270.

Kafkas S (2006). Phylogenetic analysis of the genus Pistacia by AFLP markers. Plant Systematics and Evolution 262: 113-124.

Katsiotis A, Hagidimitriou M, Drossou A, Pontikis C, Loukas M (2003). Genetic relationships among species and cultivars of Pistaciausing RAPDs and AFLPs. Euphytica 132, 279–86.

Mirzaei S, Bahar M, Sharif Nabi B,Vatan Pour Azghadi E, Bani Nasab B  (2002). Genetic diversity among Iranian pistachio cultivars using RAPD molecular marker. MSc. Thesis. Isfahan University of Technology.

Naghavi, M. Gharyazi, B. Hossini Salekdeh, G. (2009). Molecular Markers). Tehran   University Publication. 340 Pages. Noroozi Sh, Baghizadehand AM, Javaran J (2009). The genetic diversity of Iranian pistachio (Pistacia vera L.) cultivars revealed by ISSR markers. Bio.Di Con 2: 50-56.

Ozbek S, Ayfer M (1957). Pistacia tuleriuzerin desitolojikarastirmalar. Ankara Univeristies iZiraat Facultesi Yilliji. 3: 203-222.

Ozden-Tokatli Y, Akdemir H, Tilkat E, Onay A (2010). Current status and conservation of Pistacia germplasm. Biotechnology Advances  28, 130–141.

Parfitt DE, Badenes ML (1997). Phylogeny of the genus Pistacia as determined from analysis of the chloroplast genome. Proceedings of the National Academy of Sciences 94: 7987- 7992.

Pazouki L, Mardi M, SalehiShanjani P, Hagidimitrio M, Pirseyedi S.M, Naghavi MR, Avanzato D, Vendramin E, Kafkas S, Ghareyazie B, Ghaffari MR,  KhayamNekoui SM (2010). Genetic diversity and relationship among Pistacia species and cultivars. Conservation Genetics  11: 311-318.

Rohlf FJ (1992). NTSYS-pc (Numerical Taxonomy and Multivariate Analysis System).Version 1.70.Exeter, Setauket, NY.

Taghizadeh A, Ahmadi J, Haddad R, Zarrabi M. (2010). A comparative analysis of ISSR and RAPD markers for studying genetic diversity in Iranian pistachio cultivars. Iranian Journal of Genetics and plant Breeding 1: 6-16.

Vendramin E, Dettori M T, Verde I, Micali S, Giovinazzi J, Mardi M (2009). Molecular characterization  of  Pistacia genus by microsatellite markers. Acta Horticulturae 825: 55–61.

Zamani P, Akhondi M, Mohammadabadi MR (2015). Associations of Inter-Simple Sequence Repeat loci with predicted breeding values of body weight in Sheep. Small Ruminant Research 132: 123–127.

 Zamani P, Akhondi M, Mohammadabadi MR, Saki AA, Ershadi A, Banabazi MH, Abdolmohammadi AR (2011). Genetic variation of Mehraban sheepusing two inter-simple sequence repeat (ISSR) markers. African Journal of Biotechnology 10: 1812–1817.

Zohary M (1952). A monographic study of the genus Pistacia Palestine.  American Journal of Botany  5: 187-228.

 

 

 

 


Evaluation of genetic diversity in 25 Iranian pistachio genotypes using ISSR markers

 

Arjmand Ghahestani R.1, Tavassolian I.*2, Mohammadi Nejad Gh.2

 

1Department of Horticulture, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, Iran.

2Horticultural Research Institute, Shahid Bahonar University of Kerman, Kerman, Iran.

 

Abstract

In this study 17 ISSR molecular markers were employed to investigate the genetic diversity among 25 female pistachio genotypes. The markers produced  148 bands in total, of which 123 (83%) was polymorphic and the average of polymorphism for each primer was 7.23 band per primer. Cluster analyses of genotypes were estimated based on Jaccard's similarity algorithm and UPGMA clustering method and placed the genotypes in seven groups. The results obtained by similarity matrix revealed the highest similarity between cultivars "Khanjari Damghan" and "Fandoghi 48" (63%) and the lowest similarity between cultivars "Ebrahimi" and "Shasti" (2%) which is due to high amount of polymorphism among the genotypes. This experiment shows that ISSR markers are reliable for assessment of pistachio female genotypes and genetic diversity study. Moreover, these markers can be used to provide molecular genetics identity for Iranian pistachio populations.

 

Keywords: Pistachio, Genetic diversity, Jaccard similarity coefficient, Cluster analysis, ISSR.

 

 

 



* نویسنده مسئول: ایرج توسلیان                                تلفن: 03431322633                    Email: itavasoli@uk.ac.ir      

[1]- Anacardiaceae

[2]-P.vera L.

[3]- Yunt

[4]-P. integerrima

[5]-P. mexicana

[6]-P. palaestina

[7]-P. lentiscus

[8]- Unnerighted Paired Group Method Using Aritmatic Average

[9]- P. atlantica subsp. Kurdica

[10]-Principle coordinate analysis

[11]- cetyltrimethylammonium bromide

[12]principle coordinate analysis

* Corresponding Author: Tavassolian I.           Tel: 03431322633                  Email: itavasoli@uk.ac.ir

Ahmadiafzadi M, Tabatabaei S B, Mohammadi SA, Tajabadipour A (2007).  Comparison of genetic diversity in species and cultivars of Pistachio (Pistacia sp.) based on Amplified Fragment Length Polymorphism (AFLP) markers. Iranian Journal of Agricultural Biotechnology 3: 147-152.
Alami E, Taeb M, Lotfi E, Sadeghian Motahar I (2003). Study of polymorphic esterase isosymes, peroxidase and mallat dehydrogenase in Iranian pistachio cultivars and species. Agricultural and environmental science and technology 7: 107 – 113.
Al-saghir MG, Porter DM (2006). Random Amplified Polymorphic DNA (RAPD) markers of Pistacia species (Anacardiaceae). Asian Journal of Plant Sciences 5: 1002-1006.
Arabnezhad HM, Bahar A, Tajabadi Pour (2008). Assessment of genetic diversity among Iranian pistachios using microsatellites isolated from Pistacia khinjuk. Journal of Science and Technology of Agriculture and Natural Resources 45: 218-228.
Askari N, Mohammad Abadi MR, Baghizadeh A (2011). ISSR markers for assessing DNA polymorphism and genetic characterization of cattle, goat and sheep populations. Iranian Journal of Biotechnology 9: 222–229.
Baghizadeh A, Noroozi Sh, Jalili Javaran M (2010). Study on genetic diversity of some Iranian pistachio (Pistacia vera L.) cultivars using Random Amplified Polymorphic DNA (RAPD), Inter Sequence Repeat (ISSR) and Simple Sequence Repeat (SSR) markers: A comparative study. African Journal of Biotechnology 9: 7632-7640.
Barghchi M, Alderson PG (1989). Pistachio (Pistacia vera L.). In: Biotechnology in Agriculture and Forestry. Bajaj, Y.P.S. (ed.). Springer Verlag, Berlin 5: 68-98.
Baron, E, Marco L.D, Marra F.P,  Sidari M (1996). Isozymes and canonical discriminant analysis to identify pistachio (Pistacia vera L.) germplasm. Horticulture Science 31, 134-138.
Doyle JJ, Doyle JL (1987). A rapid isolation procedure for small quantities of fresh leaf tissue. Phytochem Bulletin 19:11–15.
Dyszel SM, Petitt BC (1990). Determination of the country of origin of pistachio nut by DSC and HPLC. Journal of the American Oil Chemists' Society 67:947-951.
Ehsanpour AA, Tavassoli M, Arab L (2008). Sex determination of Pistacia vera L. using ISSR markers. Malaysian Applied Biology Journal 37: 25-28.
Ghasemi M, Baghizadeh A, Mohammadabadi MR (2010). Determination ofgenetic polymorphism in Kerman Holstein and Jersey cattle population using ISSR markers. Australian Journal of Basic Applied Science 4: 5758–5760.
Golan-Goldhirsh A, Barazani O, Wang Z S , Khadka DK, Saunders JA, Kostiukovsky V, Rowland LJ (2004). Genetic relationships among Mediterranean Pistacia species evaluated by RAPD and AFLP markers. Plant Systematics and Evolution 246: 9–18.
Jaccard P (1908). Nouvellesrecherchessur la distribution.florale. Bulletin de la Société vaudoise des sciences naturelles. 44: 223-270.
Kafkas S, Perl-Treves RN, Kaska (2000). Unusual Pistacia atlantica Desf. monoecious sex types in the Yunt Mountains of the Manisa province of Turkey. Israel  Journal  of Plant Sciences 48: 277-280.
Kafkas S, Ebru K, Perl-Treves R (2002). Morphological diversity and germplasm survey of three wild Pistacia species in Turkey. Genetic Resources and Crop Evolution 49: 261-270.
Kafkas S (2006). Phylogenetic analysis of the genus Pistacia by AFLP markers. Plant Systematics and Evolution 262: 113-124.
Katsiotis A, Hagidimitriou M, Drossou A, Pontikis C, Loukas M (2003). Genetic relationships among species and cultivars of Pistaciausing RAPDs and AFLPs. Euphytica 132, 279–86.
Mirzaei S, Bahar M, Sharif Nabi B,Vatan Pour Azghadi E, Bani Nasab B  (2002). Genetic diversity among Iranian pistachio cultivars using RAPD molecular marker. MSc. Thesis. Isfahan University of Technology.
Naghavi, M. Gharyazi, B. Hossini Salekdeh, G. (2009). Molecular Markers). Tehran   University Publication. 340 Pages. Noroozi Sh, Baghizadehand AM, Javaran J (2009). The genetic diversity of Iranian pistachio (Pistacia vera L.) cultivars revealed by ISSR markers. Bio.Di Con 2: 50-56.
Ozbek S, Ayfer M (1957). Pistacia tuleriuzerin desitolojikarastirmalar. Ankara Univeristies iZiraat Facultesi Yilliji. 3: 203-222.
Ozden-Tokatli Y, Akdemir H, Tilkat E, Onay A (2010). Current status and conservation of Pistacia germplasm. Biotechnology Advances  28, 130–141.
Parfitt DE, Badenes ML (1997). Phylogeny of the genus Pistacia as determined from analysis of the chloroplast genome. Proceedings of the National Academy of Sciences 94: 7987- 7992.
Pazouki L, Mardi M, SalehiShanjani P, Hagidimitrio M, Pirseyedi S.M, Naghavi MR, Avanzato D, Vendramin E, Kafkas S, Ghareyazie B, Ghaffari MR,  KhayamNekoui SM (2010). Genetic diversity and relationship among Pistacia species and cultivars. Conservation Genetics  11: 311-318.
Rohlf FJ (1992). NTSYS-pc (Numerical Taxonomy and Multivariate Analysis System).Version 1.70.Exeter, Setauket, NY.
Taghizadeh A, Ahmadi J, Haddad R, Zarrabi M. (2010). A comparative analysis of ISSR and RAPD markers for studying genetic diversity in Iranian pistachio cultivars. Iranian Journal of Genetics and plant Breeding 1: 6-16.
Vendramin E, Dettori M T, Verde I, Micali S, Giovinazzi J, Mardi M (2009). Molecular characterization  of  Pistacia genus by microsatellite markers. Acta Horticulturae 825: 55–61.
Zamani P, Akhondi M, Mohammadabadi MR (2015). Associations of Inter-Simple Sequence Repeat loci with predicted breeding values of body weight in Sheep. Small Ruminant Research 132: 123–127.
 Zamani P, Akhondi M, Mohammadabadi MR, Saki AA, Ershadi A, Banabazi MH, Abdolmohammadi AR (2011). Genetic variation of Mehraban sheepusing two inter-simple sequence repeat (ISSR) markers. African Journal of Biotechnology 10: 1812–1817.
Zohary M (1952). A monographic study of the genus Pistacia Palestine.  American Journal of Botany  5: 187-228.