نوع مقاله : مقاله پژوهشی
نویسندگان
اﻳﺮان، ﻛﺮﻣﺎن، داﻧﺸﮕﺎه ﺗﺤﺼﻴﻼت ﺗﻜﻤﻴﻠﻲ ﺻﻨﻌﺘﻲ و ﻓﻨﺎوری ﭘﻴﺸﺮﻓﺘﻪ، ﭘﮋوﻫﺸﮕﺎه ﻋﻠﻮم و ﺗﻜﻨﻮﻟﻮژی ﭘﻴﺸﺮﻓﺘﻪ و ﻋﻠﻮم ﻣﺤﻴﻄﻲ، ﮔﺮوه ﺑﻴﻮﺗﻜﻨﻮﻟﻮژی
چکیده
ارزیابی دقیق مقدار و پراکنش تنوع ژنتیکی در گیاهان به منظور تدوین راهبرد حفاظت و استفاده از منابع ژنتیکی از اهمیت ویژهای برخوردار است. در این پژوهش روابط ژنتیکی بین 40 رقم گندم نان (Triticum aestivum L.) با استفاده از 10 آغازگر ISSR و 4 آغازگر RAPD مورد بررسی قرار گرفت. 10 آغازگر ISSR تعداد 92 باند چندشکل با میانگین 2/9 باند چندشکل برای هر آغازگر تولید نمودند و میانگین درصد چندشکلی 79/86% بود. همچنین چهار آغازگر نشانگر RAPD تعداد 19 باند چندشکل با میانگین 75/4 باند چندشکل برای هر آغازگر تولید کردند و میانگین درصد چندشکلی 85/67% بود. بالا بودن معیارهای تعداد آلل مؤثر، تنوع ژنی، شاخص شانون، محتوای اطلاعات چندشکل، در مورد نشانگرهای ISSR برای پرایمرهای UBC813، UBC816 و UBC824 و در مورد نشانگرهای RAPD برای پرایمرهای OH-04 و OPA02 نشاندهنده کارایی بالای این آغازگرها در تمایز ژنوتیپهای گندم مورد مطالعه در این تحقیق بود. تجزیه خوشهای به روش UPGMA و بر اساس ضریب تشابه جاکارد بر اساس کل باندهای چندشکل نشانگر ISSR و RAPD ارقام گندم را در دو گروه جداگانه قرار داد. همچنین تجزیه به مختصات اصلی و تعیین ساختار جمعیت بر اساس روش بیزی با نرمافزار STRUCTURE ارقام را در دو گروه قرار دادند و نتایج تجزیه خوشهای را تایید نمودند. تجزیه واریانس مولکولی نشان داد که تنوع درون گروهها بیشتر از بین گروهها است. اطلاعات به دست آمده از این تحقیق و تنوع موجود در ارقام گندم میتواند در پژوهشهای بهنژادی برای افزایش عملکرد مورد استفاده قرار گیرد.
کلیدواژهها
موضوعات
عنوان مقاله [English]
Evaluation of genetic diversity of different bread wheat (Triticum aestivum L.) varieties using molecular markers ISSR and RAPD
نویسندگان [English]
Dept. of Biotechnology, Institute of Science and High Technology and Environmental Sciences, Graduate University of Advanced Technology, Kerman, I.R. of Iran.
چکیده [English]
Exact evaluation of the amount and distribution of genetic diversity in plants is essential for developing a strategy for conservation and use of genetic resources. Study, the genetic diversity of 40 genotypes of bread wheat (Triticum aestivum L.) was investigated using 10 Inter Simple Sequence Repeat (ISSR) markers and four Random Amplified Polymorphic DNA (RAPD) markers. The 10 ISSR primers produced 92 polymorphic bands with an average of 9.2 polymorphic bands for each primer and the polymorphic percentage was 86.79%. Also, four RAPD primers produced 19 polymorphic bands with an average of 4.75 bands for each primer, and the average polymorphic percentage was 67.85%. High levels of effective number of alleles, Nei index, Shannon index, and the Polymorphism information content of ISSR markers for UBC813, UBC816 and UBC824 primers, and RAPD markers for OH-04 and OPA02 primers was indicated the high efficiency of these primers in differentiating the studied wheat genotypes in this research. Cluster analysis using UPGMA method based on Jaccard's similarity coefficient based on the whole polynomial bands of ISSR and RAPD markers divided wheat genotypes in two separate groups. Also, principal coordinates analysis and determining population structure based on Bayesian method with STRUCTURE software placed the cultivars in two groups and confirmed cluster analysis results. Molecular analysis of variance showed that within groups’ diversity is more than between groups’ diversity. Information obtained from this research and variation in wheat cultivars can be used in the breeding programs for yield improvement.
کلیدواژهها [English]
بررسی تنوع ژنتیکی ارقام مختلف گندم نان (Triticum aestivum L.)با استفاده از نشانگرهای مولکولی ISSR و RAPD
مژده رحمانی، مهدی رحیمی*، مریم عبدلی نسب و محمود ملکی
ایران، کرمان، دانشگاه تحصیلات تکمیلی صنعتی و فناوری پیشرفته، پژوهشگاه علوم و تکنولوژی پیشرفته و علوم محیطی، گروه بیوتکنولوژی
تاریخ دریافت: 4/12/1397 تاریخ پذیرش: 19/1/1398
چکیده
ارزیابی دقیق مقدار و پراکنش تنوع ژنتیکی در گیاهان به منظور تدوین راهبرد حفاظت و استفاده از منابع ژنتیکی از اهمیت ویژهای برخوردار است. در این پژوهش روابط ژنتیکی بین 40 رقم گندم نان (Triticum aestivum L.) با استفاده از 10 آغازگر ISSR و 4 آغازگر RAPD مورد بررسی قرار گرفت. 10 آغازگر ISSR تعداد 92 باند چندشکل با میانگین 2/9 باند چندشکل برای هر آغازگر تولید نمودند و میانگین درصد چندشکلی 79/86 درصد بود. همچنین چهار آغازگر نشانگر RAPD تعداد 19 باند چندشکل با میانگین 75/4 باند چندشکل برای هر آغازگر تولید کردند و میانگین درصد چندشکلی 85/67 درصد بود. بالا بودن معیارهای تعداد آلل مؤثر، تنوع ژنی، شاخص شانون، محتوای اطلاعات چندشکل، در مورد نشانگرهای ISSR برای آغازگرهای UBC813، UBC816 و UBC824 و در مورد نشانگرهای RAPD برای آغازگرهای OH-04 و OPA02 نشاندهنده کارایی بالای این آغازگرها در تمایز ژنوتیپهای گندم مورد مطالعه در این تحقیق بود. تجزیه خوشهای به روش UPGMA و بر اساس ضریب تشابه جاکارد بر اساس کل باندهای چندشکل نشانگر ISSR و RAPD ارقام گندم را در دو گروه جداگانه قرار داد. همچنین تجزیه به مختصات اصلی و تعیین ساختار جمعیت بر اساس روش بیزی با نرمافزار STRUCTURE ارقام را نیز در دو گروه قرار دادند و نتایج تجزیه خوشهای را تأیید نمودند. تجزیه واریانس مولکولی نشان داد که تنوع درون گروهها بیشتر از بین گروههاست. اطلاعات به دست آمده از این تحقیق و تنوع موجود در ارقام گندم میتواند در پژوهشهای بهنژادی برای افزایش عملکرد مورد استفاده قرار گیرد.
واژههای کلیدی: تجزیه خوشهای، تجزیه به مختصات اصلی، ساختار جمعیت، روش UPGMA.
* نوسینده مسئول، تلفن تماس: 09173079874، پست الکترونیکی: me.rahimi@kgut.ac.ir
مقدمه
گندم منبع اصلی کالری و پروتئین برای بخش بزرگی از جمعیت دنیا و همچنین ایران محسوب میشود و سطح زیر کشت گندم در ایران در سال 95-96 برابر با 5437804 هکتار با تولیدی به میزان 12400000 تن بوده است (1) ولی با این وجود کشور سالانه میزان زیادی گندم وارد می کند که میزان واردات در سال 93 برابر با 7122870 تن بوده است (6). گندم دارای گونههای زیادی است ولی بیشترین سطح زیر کشت (90 درصد) و بیشترین میزان تولید (94 درصد) مربوط به گونۀ T. aestivum یا گندم نان است با وجود اینکه کشور ایران از غنای ژنتیکی بالایی برخوردار است و تنوع ژنوتیپها، نژادها و جمعیتهای گیاهان زراعی در گذشته در آن زیاد بوده است، ولی در حال حاضر این تنوع به شدت کاهش یافته است و ارقام زراعی اندکی بخش عمده تولید هر محصول را به خود اختصاص میدهند (8). بنابراین برای کاهش میزان واردات گندم با توجه به کمبود زمینهای قابل کشت باید به دنبال افزایش تولید در هکتار از طریق افزایش عملکرد ارقام جدید بود. این اتفاق زمانی به وجود میآید که امکان انتخاب مواد گیاهی مناسب در برنامههای اصلاحی وجود داشته باشد و این مواد گیاهی نیز دارای تنوع کافی باشند. تنوع مبنای همه گزینشهاست و انتخاب ژنوتیپی نیز نیازمند تنوع است (19). میزان تنوع ژنتیکی موجود در گیاهان زراعی و خویشاوندان وحشی آنها جهت انتخاب والدین مناسب در برنامه اصلاحی کلید موفقیت اصلاحگران است و بنابراین حفظ تنوع ژنتیکی در برنامههای اصلاحی بسیار مهم میباشد (19) و یک برنامه اصلاحی زمانی موفق است که دو عامل تنوع و انتخاب در گیاه مورد آزمایش وجود داشته باشد. روشهای مختلفی مانند انواع نشانگرهای مورفولوژیکی، بیوشیمیایی و مولکولی یا DNA جهت تعیین و برآورد تنوع ژنتیکی گیاهان وجود دارد. در بررسی تنوع ژنتیکی دادههای حاصل از نشانگرهای مولکولی نسبت به دادههای مورفولوژیکی، شجرهنامهها، بیوشیمیایی برتری دارند (17).
به کارگیری نشانگرهای مولکولی در گندم کاربرد گستردهای در دهه گذشته داشته است. از جمله نشانگرهایی که برای بررسی تنوع ژنتیکی در گندم استفاده شده نشانگرهای ISSR و RAPD میباشد (22 و 24). این نشانگرها مبتنی بر PCR بوده و به طراحی آغازگر نیازی ندارند. نشانگر ISSR قطعه DNAی موجود و قابل تکثیر بین دو ریز ماهواره را که در جهت مخالف هم قرار دارند تکثیر میکند (24) و نشانگر RAPD یک نشانگر تصادفی بوده و نواحی تصادفی در ژنوم را تکثیر میکند. این نشانگرها در مطالعات تنوع ژنتیکی، فیلوژنی، نشاندار کردن ژنی، نقشهیابی ژنومی و زیستشناسی تکاملی به کار میروند (21). خلستکنیا و همکاران (14) تنوع ژنتیکی 54 رقم قدیمی و جدید گندم بهاره را با 23 جایگاه ریز ماهواره ارزیابی و در مجموع 151 آلل با تعداد آلل بین 11-3 آلل و میانگین 9/6 آلل برای هر جایگاه گزارش کردند. در تحقیقی تنوع ژنتیکی ژرم پلاسم گندم دوروم (Triticum turgidum ssp. durum) با استفاده از توالیهای ساده تکراری (ریزماهوارهها) ارزیابی شد. یکصد و چهل ژنوتیپ متعلق به مناطق مختلف (10 کشور) با 48 جفت نشانگر ریزماهواره (با پوشش یکنواخت ژنوم) مورد بررسی قرار گرفتند و از این تعداد، 37 جفت نشانگر چند شکلی نشان دادند. در مجموع از 37 جفت نشانگر، 245 آلل چند شکل متفاوت، در دامنهای بین 2 تا 17 آلل و با متوسط 9/6 آلل در هر جایگاه ژنی تکثیر گردید. میزان اطلاعات چندشکل (PIC) از 007/0 (جایگاه ژنی Xgwm274) تا 799/0 (جایگاه ژنی Xgwm427) متفاوت بود (4). در تحقیقی دیگر ارزیابی تنوع ژنتیکی ژنوتیپهای گندم با استفاده از مارکرهای ISSR انجام شد. نتایج نشان داد که بیشترین مقدار شاخص محتوای اطلاعات چندشکلی (PIC) مربوط به پرایمر IS-10 (36/0) و کمترین آن مربوط به پرایمرهای IS-16 و IS-11 (22/0) بود. تجزیه خوشهای به روش UPGMA و بر اساس ضریب تشابه دایس ژنوتیپها را در سه گروه دستهبندی کرد که ژنوتیپ شماره 25 (شاهد گندم نان) در یک گروه کاملاً جداگانه با فاصله زیادی از سایر گروهها قرارگرفت (3). در بررسی الگوی تغییرات ژنتیکی درون و بین دو مجموعه ارقام بومی گندم نگهداری شده در مرکز بین المللی اصلاح گندم و ذرت (CIMMYT) با استفاده از 76 نشانگر ریز ماهواره، تنوع ژنتیکی قابل ملاحظهای درون تودههای مکزیکی و ترکیه مشاهده شد (11). سفالیان و همکاران (23) به منظور بررسی تنوع ژنتیکی 9 رقم گندم بومی بهاره از 11 نشانگر ISSR استفاده کردند که در مجموع 108 باند تشکیل شد که 78 قطعه (22/72 درصد) از آنها چندشکل بودند.
با توجه به اینکه این ارقام با هم و با این آغازگرها مورد مطالعه قرار نگرفتهاند و شباهت و تفاوت این ارقام با هم بررسی نشده است. بنابراین اهداف این مطالعه بررسی تنوع ژنتیکی ارقام گندم نان از طریق نشانگرهای مولکولی ISSR و RAPD و تعیین روابط ژنتیکی آنها برای شناسایی والدین مناسب برای برنامههای اصلاحی و ژنتیکی میباشد. تا بتوان از ارقامی که تفاوت زیادی باهم دارند در برنامههای دورگگیری برای افزایش عملکرد گندم و مناسب کشت در منطقه ماهان استفاده نمود.
مواد و روشها
برای بررسی تنوع ژنتیکی بذور تعداد 40 رقم گندم زراعی که از بانک ژن تهیه گردیده بودند (جدول 1) به مدت 24 ساعت در آب و در تاریکی در دمای اتاق خیسانده شدند و سپس در پتریدیش کشت شدند و پس از 12 روز برگها برای استخراج DNA برداشت شدند.
جدول 1- ارقام گندم نان مورد مطالعه
ردیف |
ژنوتیپ |
ردیف |
ژنوتیپ |
ردیف |
ژنوتیپ |
ردیف |
ژنوتیپ |
1 |
اکبری |
11 |
زرین |
21 |
الوند |
31 |
رسول |
2 |
نیک نژاد |
12 |
تجن |
22 |
شیراز |
32 |
شعله |
3 |
روشن |
13 |
فلات |
23 |
اکسکلیبر |
33 |
سبلان |
4 |
چمران |
14 |
آذر 2 |
24 |
بزاستایا |
34 |
مرودشت |
5 |
هیرمند |
15 |
مغان 3 |
25 |
داراب |
35 |
آرتا |
6 |
مغان 1 |
16 |
امید |
26 |
هامون |
36 |
کرج 2 |
7 |
دز |
17 |
ارگ |
27 |
گاسپارد |
37 |
کویر |
8 |
آستارا |
18 |
سرداری |
28 |
وی ناک |
38 |
بولانی |
9 |
طوس |
19 |
الموت |
29 |
شاه پسند |
39 |
mv-17 |
10 |
سرخ تخم |
20 |
آزادی |
30 |
پیشتاز |
40 |
مهدوی |
جدول 2- خصوصیات آغازگرهای ISSR و RAPD مورد مطالعه
درصد GC |
دمای اتصال |
توالی |
نام آغازگرها |
نوع آغازگر |
47 |
42 |
3´ |
UBC812 |
ISSR |
47 |
42 |
T3´ |
UBC813 |
|
53 |
43 |
G3´ |
UBC815 |
|
47 |
47 |
T3´ |
UBC816 |
|
47 |
47 |
A3´ |
UBC817 |
|
53 |
44 |
C3´ |
UBC823 |
|
53 |
45 |
G3´ |
UBC824 |
|
47 |
48 |
T3´ |
UBC825 |
|
53 |
49 |
C3´ |
UBC826 |
|
47 |
43 |
A3´ |
UBC876 |
|
70 |
34 |
5´GGAAGTCGCC3´ |
OH-04 |
RAPD |
60 |
32 |
5´TCGGCGATAG3´ |
OA12 |
|
60 |
32 |
5´CCGCATCTAC3´ |
OC4 |
|
70 |
30 |
5՜TGCCGAGCTG3՜ |
OPA02 |
برگهای برداشت شده بلافاصله درون نیتروژن مایع قرار گرفته و پس از انجماد سریع به فریزر در 80- درجه سانتی گراد منتقل شدند. استخراج DNA از برگهای جوان به روش دلاپورتا با اندکی تغییرات انجام شد (10). کمیت و کیفیت DNA استخراجی از طریق اسپکتروفتومتری و مطالعه طیف جذبی DNA در طول موجهای 260 و 280 نانومتر و نیز الکتروفورز ژل آگارز یک درصد (13، 16 و 25) تعیین شد. واکنش زنجیرهای پلیمراز با استفاده از 10 آغازگر ISSR و چهار آغازگر RAPD انجام شد (جدول 2). واکنش PCR با استفاده از مواد تهیه شده از شرکت سینا ژن و در حجم 10 میکرولیتر شامل یک میکرولیتر بافر 10X PCR، 25/0 میکرولیتر کلریدمنیزیم 50 میلیمولار، U 25/0 آنزیم Tag DNA Polymerase، 75/0 میکرولیتر از هر آغازگر با غلظت 10 میکرومولار، 25/0 میکرولیتر محلول نوکلئوتیدی (dNTPs) 10 میلیمولار و دو میکرولیتر DNA ژنومی (20 نانوگرم) در دستگاه ترموسایکلر (مدل Biometra) انجام شد.
الگوی دمایی واکنش زنجیرهای پلی مراز برای آغازگر ISSR به این شرح بود: اولین چرخه حرارتی شامل چهار دقیقه واسرشتهسازی اولیه در 94 درجه سانتی گراد انجام گرفت، سپس 35 سیکل بعدی بهصورت 40 ثانیه واسرشتهسازی در 94 درجه، 40 ثانیه مرحله اتصال آغازگر در دمای مناسب (بسته به آغازگر متفاوت بود)، دو دقیقه مرحله بسط در دمای 72 درجه صورت گرفت و آخرین چرخه هم بهصورت پنج دقیقه بسط انتهایی در دمای 72 درجه انجام گرفت، سپس نگهداری در دمای 4 درجه بود. همچنین الگوی دمایی برای آغازگر RAPD نیز به این شرح بود: اولین چرخه حرارتی شامل پنج دقیقه واسرشتهسازی اولیه در 94 درجه، سپس 35 سیکل بهصورت یک دقیقه واسرشتهسازی در 94 درجه، 45 ثانیه مرحله اتصال آغازگر در دمای TM (بسته به آغازگر متفاوت بود)، یک دقیقه مرحله بسط در دمای 72 درجه و پنج دقیقه بسط انتهایی در دمای 94 درجه، سپس نگهداری در دمای در 4 درجه بود. محصولات PCR با استفاده از ژل آگارز یک درصد (13، 16 و 25) با ولتاژ 70 ولت طی یک ساعت الکتروفورز تفکیک شدند. رنگآمیزی با اتیدیوم بروماید به مدت 15 دقیقه انجام شد. در انتها ژل با آب مقطر شسته شده و توسط دستگاه ژل داک (مدل Bio-Rad) عکسبرداری انجام شد. آغازگرهای چند شکل به دست آمده روی 40 ژنوتیپ به صورت دادههای صفر و یک وارد نرم افراز Excel و سپس به وسیله نرم افزارهای Power Marker، STRUCTURE، PAST، PopGene 32 و SAS مورد تجزیه و تحلیل قرار گرفتند. جهت به دست آوردن تعداد آلل مؤثر، شاخص شانون، تنوع ژنی نی و شاخص اطلاعات نشانگری از نرم افزارهای Power Marker و PopGene32 استفاده شد. همچنین تجزیه خوشهای، ضریب همبستگی کوفنتیک و تجزیه به مختصات اصلی با استفاده از نرم افزار PAST انجام شد. تجزیه تابع تشخیص با استفاده از نرمافزار SAS و روش cross-validation انجام شد. همچنین از نرمافزار STRUCTURE برای تأیید نتایج تجزیه خوشهای و همچنین تعیین تعداد گروه و ساختار آن استفاده شد. در نهایت از نرمافزار GenAlEx 6.5 برای بررسی تنوع درون و بین گروههای به دست آمده از نرمافزار STRUCTURE استفاده شد.
نتایج
در شکل 1 تکثیر پرایمر UBC826 (قسمت الف) برای نشانگر ISSR و پرایمر OH-04 (قسمت ب) برای نشانگر RAPD نشان داده شده است. از مجموع 10 آغازگر ISSR مورد استفاده تعداد 92 باند چندشکل بهدست آمد که 4 آغازگر 100 درصد چندشکلی نشان دادند و میانگین باندهای چندشکل 2/9 و درصد چندشکلی 79/86 درصد بهدست آمد و از 4 آغازگر RAPD مورد استفاده تعداد 19 باند چندشکل حاصل شد که یک آغازگر 100 درصد چندشکلی نشان داد و میانگین باندهای چندشکل 75/4 و درصد چند شکلی 85/67 درصد بهدست آمد (جدول 3).
(الف) (ب)
شکل 1- نمونه الگوی باندی پرایمر UBC826 (الف) نشانگر ISSR و پرایمر OH-04 (ب) نشانگر RAPD حاصل از تکثیر ارقام گندم
جدول 3- نتایج تجزیه مولکولی برای آغازگرهای ISSR و RAPD در ارقام گندم
نسبت چندگانه موثر |
شاخص نشانگری |
شاخص شانون |
تنوع ژنی نی |
تعداد آلل مؤثر |
محتوای اطلاعات چندشکل |
درصد چندشکلی |
تعداد باندهای چندشکل |
نام آغازگرها |
8 |
08/2 |
5/0 |
33/0 |
56/1 |
26/0 |
100 |
8 |
UBC812 |
7 |
2/1 |
57/0 |
38/0 |
65/1 |
3/0 |
100 |
7 |
UBC813 |
66/6 |
46/1 |
41/0 |
28/0 |
52/1 |
22/0 |
6/66 |
10 |
UBC815 |
8 |
24/2 |
52/0 |
35/0 |
61/1 |
28/0 |
100 |
8 |
UBC816 |
33/8 |
5/1 |
35/0 |
21/0 |
31/1 |
18/0 |
3/83 |
10 |
UBC817 |
26/11 |
91/1 |
34/0 |
21/0 |
32/1 |
17/0 |
6/86 |
13 |
UBC823 |
14 |
78/3 |
51/0 |
34/0 |
59/1 |
27/0 |
100 |
14 |
UBC824 |
5/4 |
21/1 |
38/0 |
25/0 |
41/1 |
27/0 |
75 |
6 |
UBC825 |
12/6 |
4/1 |
43/0 |
27/0 |
42/1 |
23/0 |
5/87 |
7 |
UBC826 |
36/7 |
76/1 |
45/0 |
3/0 |
52/1 |
24/0 |
8/81 |
9 |
UBC876 |
12/8 |
94/1 |
44/0 |
29/0 |
49/1 |
24/0 |
79/86 |
2/9 |
میانگین (Average) |
12/6 |
77/1 |
54/0 |
38/0 |
71/1 |
29/0 |
5/87 |
7 |
OH-04 |
5 |
65/0 |
26/0 |
14/0 |
17/1 |
13/0 |
100 |
5 |
OA12 |
6/1 |
19/0 |
23/0 |
15/0 |
27/1 |
12/0 |
40 |
4 |
OC4 |
8/1 |
36/0 |
38/0 |
26/0 |
49/1 |
2/0 |
60 |
3 |
OPA02 |
63/3 |
74/0 |
35/0 |
23/0 |
41/1 |
18/0 |
85/67 |
75/4 |
میانگین (Average) |
تجزیه خوشهای بر اساس مجموع باندهای چندشکل نشانگرهای RAPD و ISSR (111 باند چندشکل) بر اساس روشهای مختلف (نزدیکترین همسایه، دورترین همسایه، وارد و UPGMA) و معیارهای تشابه مختلف (جاکارد، دایس، تطابق ساده و ..) انجام شد که در نهایت روش UPGMA با معیار تشابه جاکارد به دلیل داشتن ضریب کوفنتیک بالاتر (78/0) انتخاب گردید. برای تعیین تعداد گروهها روشهای مختلفی انجام شد و براساس پرش ناگهانی در فاصله ادغام برش انجام شد که در نهایت با توجه به ناحیه برش، ارقام در دو گروه قرار گرفتند ( شکل 2). جهت اطمینان از صحیح بودن مکان خط برش و یا صحیح قرار گرفتن هر ژنوتیپ در دسته مربوط به خود در تجزیه خوشهای، از تجزیه تابع تشخیص استفاده شد. بر این اساس گروهبندی بر اساس تعداد باندهای چندشکل هر دو نشانگر، صدرصد تأیید شد و درصد خطای گروهبندی برای تمام کلاسترها صفر درصد بوده و در نتیجه تمام کلاسترها 100 درصد درست گروهبندی شده بودند (جدول 4).
شکل 2- دندروگرام حاصل از تجزیه خوشهای بر اساس دادههای حاصل از نشانگر ISSR و RAPD در ارقام گندم نان
جدول 4- جدول نسبت موفقیت افراد درون گروهها با تابع تشخیص
تعداد افراد در هر گروه |
گروهها |
تعداد افراد پیشبینی شده در هر گروه با تابع تشخیص |
|
1 |
2 |
||
23 |
1 |
23 100% |
0 |
17 |
2 |
0 |
17 100% |
صحت درستی گروهبندی براساس دو گروه |
100% |
برای داشتن دیدگاه بهتر راجع به فواصل ژنتیکی بین ارقام و همچنین برای مشاهده فواصل ژنتیکی بین ارقام بهصورت چندبعدی، تجزیه به مختصات اصلی نیز برای تکمیل نتایج روش تجزیه خوشهای انجام شد. از تجزیه به مختصات اصلی بر روی دادههای حاصل از هر دو نشانگر یک پلات دو بعدی (شکل 3) بهدست آمد.
شکل 3- پلات دوبعدی حاصل از تجزیه به مختصات اصلی با نشانگر ISSR و RAPD
شکل 4- نمودار دو طرفه برای تعیین تعداد بهینه گروه (K) بر اساس آماره ΔK
دو مؤلفه اصلی که بیشترین سهم را در ایجاد تنوع داشتهاند به ترتیب 47/16 و 38/10 درصد از کل تغییرات را توجیه کردند. با توجه به پراکنش ارقام براساس این دو مؤلفه، ارقام در دو گروه قرار گرفتند (شکل 3). همچنین ساختار جمعیت و تعداد گروهها با نرمافزار STRUCTURE براساس روش آماری بیزی مورد بررسی قرار گرفت. در این تحقیق ساختار ژنتیکی جمعیت و چگونگی آن و تعیین تعداد مناسب گروهها (K) بر اساس روش ایوانو و همکاران (12) مشخص شد. این روش بر اساس آمار ΔK استوار است که شیب تابع احتمالی LnP(D) را در نقطهای میشکند که تعداد K فرضی در آن نقطه دارای حداکثر احتمال باشد. بر اساس نتایج تعداد گروه مناسب برای این ژرمپلاسم دو بود (شکل 4) و دو گروه احتمالی و نحوه قرار گرفتن ارقام در آنها در شکل 5 نشان داده شده است. براساس دو گروه به دست آمده تنوع درون و بین دو گروه با تجزیه واریانس مولکولی براساس تعداد باندهای چندشکل هر دو نشانگر مورد بررسی قرار گرفت (جدول 5). نتایج نشان داد که واریانس درون گروهها نسبت به واریانس بین گروهها بیشتر میباشد. بهطوریکه تنوع بین گروهی 16 درصد و تنوع درون گروهی 84 درصد از تنوع کل را توجیه نمودند.
شکل 5- تجزیه کلاستر مبتنی بر مدل Bayesian برای 40 رقم گندم مورد مطالعه بر اساس 111 باند چندشکل.
هر رنگ یک زیرجمعیت یا کلاستر را نشان میدهد. محور عمودی ضریب تعلق هر فرد به هر کلاستر را نشان میدهد.
جدول 5- تجزیه واریانس مولکولی بر اساس نشانگرهای ISSR و RAPD
واریانس نسبی تجمعی (%) |
واریانس |
میانگین مربعات |
مجموع مربعات |
درجه آزادی |
منابع تغییرات |
16% |
325/3 |
004/81 |
004/81 |
1 |
بین گروهها |
84% |
168/17 |
168/17 |
396/652 |
38 |
درون گروهها |
100% |
493/20 |
|
40/733 |
39 |
کل |
162/0 |
PhiPT ارزش آماره |
||||
001/0 |
PhiPT احتمال آماره |
بحث
از آنجایی که یکی از معیارهای مهم در انتخاب آغازگرهای مناسب و سودمند، تعداد آللهای مؤثر است. نتایج حاصل از بررسی تعداد آلل مؤثر برای تمام جایگاههای ژنی نشان داد میانگین تعداد آلل مؤثر در کل جمعیت برای نشانگرهای ISSR، 49/1 به دست آمد. بیشترین تعداد آلل مؤثر به ترتیب مربوط به آغازگرهای UBC813، UBC816 و کمترین برای آغازگر UBC817 در بین کل ارقام بود. همچنین میانگین تعداد آلل مؤثر برای نشانگرهای RAPD در جمعیت موردمطالعه 41/1 و بیشترین مقدار مربوط به آغازگر OH-04 و کمترین برای آغازگر OA12 بهدست آمد. میزان اطلاعات چندشکل (PIC) قدرت تفکیک یک نشانگر را به واسطه تعداد آللهای چندشکل و فراوانی نسبی این آللها در جمعیت تحت مطالعه نشان میدهد. این شاخص در نشانگرهای ISSR برای جایگاه ژنی UBC813 در بیشترین مقدار (30/0) و برای جایگاه ژنی UBC823 در کمترین مقدار (17/0) و برای نشانگر RAPD برای جایگاه ژنی OH-04 بیشترین مقدار (29/0) و برای جایگاه ژنی OC4 کمترین مقدار (12/0) به دست آمد. شاخص اطلاعات شانون (I) و تنوع ژنتیکی نی (Nei) برای هر دو نشانگر نیز محاسبه شد. شاخص شانون معیاری برای اندازهگیری میزان آللهای چندشکل در یک جایگاه است. این شاخص برای نشانگر ISSR در جایگاه ژنی UBC813 بیشترین مقدار (57/0) و برای جایگاه ژنی UBC823 کمترین مقدار (34/0) و برای نشانگرهای RAPD در جایگاه ژنی OH-04 بیشتر از سایر جایگاهها (54/0) و برای جایگاه ژنی OC4 کمتر از سایر جایگاهها (23/0) به دست آمد. یکی دیگر از مهمترین شاخصها برای ارزیابی تنوع ژنی در بین ارقام و جمعیتها، شاخص تنوع ژنی نی میباشد. تنوع ژنتیکی نی برای نشانگر ISSR در جایگاه ژنی UBC813 به میزان 38/0 ( بیشترین مقدار) و برای جایگاه ژنی UBC823 به میزان 21/0 (کمترین مقدار) و برای نشانگر RAPD برای جایگاه ژنی OH-04 به میزان 38/0 (بیشترین مقدار) و برای جایگاه ژنی OA12 به میزان 14/0 (کمترین مقدار) ارزیابی شد. با توجه به نتایج حاصل میتوان بیان داشت که جایگاههای ژنی UBC813 در نشانگر ISSR و OH-04 در نشانگر RAPD به دلیل داشتن بیشترین میزان شاخص شانون، تنوع ژنتیکی نی، تعدا آلل مؤثر و محتوای اطلاعات چندشکل بیشترین تنوع را بر روی جمعیت مورد مطالعه داشته است. در آزمایشی از نشانگرهای SSR برای تنوع ژنتیکی گندم استفاده شد و از مجموع 10 آغازگر انتخابی، 9 آغازگر چند شکلی قابل توجهی نشان داده و 31 نوار چند شکل با میانگین 4/3 نوار به ازای هر آغازگر مشاهده شد (9). همچنین در مطالعه دیگری تنوع ژنتیکی گندم با استفاده از نشانگر ISSR بررسی شد و بیشترین مقدار شاخص محتوای اطلاعات چندشکلی (PIC) مربوط به آغازگر 6 (36/0) و کمترین آن مربوط به آغازگرهای 1 و 5 (22/0) بود و تجزیه خوشهای به روش UPGMA و بر اساس ضریب تشابه دایس ژنوتیپها را در سه گروه دستهبندی کرد که ژنوتیپ شماره 25 (شاهد گندم نان) در یک گروه کاملاً جداگانه با فاصله زیادی از سایر گروهها قرارگرفت (3). همچنین در تحقیقی که بر روی لاینها و ارقام گندم نان و دو رقم گندم دوروم با نشانگر نیمه تصادفی اینترونی-اگزونی انجام شد، کارایی این نشانگرها در بررسی تنوع گندم نشان داده شد و توانست ارقام و لاینهای بهاره و زمستانه، و گندم نان و دوروم را به خوبی از همدیگر تفکیک نماید (2). همچنین کارایی نشانگرهای ISSR در گروهبندی ارقام و لاینهای ذرت در مطالعه دیگری نشان داه شد و توانست این نشانگر به خوبی لاینها را از هم تفکیک نماید (7) و در این تحقیق نیز این نشانگر در تفکیک ارقام توانایی خود را نشان داد.
تجزیه خوشهای بر اساس نشانگرهای ISSR و RAPD ارقام را در دو گروه تقسیم کرد. گروه اول شامل ارقام اکبری، نیک نژاد، روشن، چمران، هیرمند، مغان 1، دز، آستارا، طوس، سرخ تخم، زرین، تجن، فلات، آذر2، مغان3، امید، ارگ، سرداری، الموت، آزادی، الوند، شیراز، رسول و گروه دوم شامل ارقام اکسکلیبر، بزاستایا، داراب، هامون، گاسپارد، وی ناک، شاه پسند، پیشتاز، شعله، سبلان، مرودشت، آرتا، کرج 2، کویر، بولانی، mv-17 و مهدوی بود. نتایج حاصل از میزان شباهت ژنتیکی در مجموع نشان داد که بیشترین میزان شباهت بهترتیب بین ارقام اکسکلیبر و هامون (44/97 درصد) و مرودشت و مهدوی (05/96 درصد) محاسبه گردید. در مقابل بیشترین فاصله ژنتیکی یا کمترین شباهت مربوط به ارقام دز و داراب (6/46 درصد) و دز و شعله (52/48 درصد) محاسبه گردید.
با استفاده از اطلاعات بهدستآمده از پلات دو بعدی تجزیه به مختصات اصلی براساس نشانگر ISSR و RAPD نتیجه گرفته شد که این گروهبندی با گروهبندی حاصل با دندروگرام خوشهای مطابقت بالایی دارد ولی با توجه به اینکه درصد تجمعی این دو مؤلفه (85/26) پایین میباشد ولی توانستهاند نتایج گروهبندی تجزیه خوشهای را تأیید نماید به این احتمال است که این نشانگرها در طول ژنوم به خوبی پراکنده شده و ژنوم را پوشش دادهاند (5). همچنین نتایج گروهبندی به روش بیزی با نرمافزار STRUCTURE نیز نشان داد که دو گروه برای این ارقام وجود دارد و نتایج تجزیه خوشهای را تأیید نمود. تعداد ارقام در دو گروه با تجزیه خوشهای به ترتیب 23 و 17 ژنوتیپ بود درحالی که در روش بیزی در دو گروه به ترتیب 24 و 16 ژنوتیپ قرار گرفته بودند و تنها یک ژنوتیپ در گروهی دیگر قرار گرفته بود و در حدود 5/97 درصد دو گروه باهم شباهت داشتند. تجزیه ساختار با نشانگرهای غالب (مانند ISSR و ...) در مطالعات دیگر هم انجام شده است (15، 18، 20 و 26). تجزیه واریانس مولکولی بر اساس نشانگرهای ISSR و RAPD نشان داد که واریانس درون گروهها نسبت به واریانس بین گروهها بیشتر میباشد بهطوری که تنوع بین جمعیتی 16 درصد و تنوع درون جمعیتی 84 درصد از تنوع کل را توجیه نمودند و این نتیجه نشان میدهد که بر اساس نشانگرهای ISSR و RAPD تفاوت داخل گروهها بیشتر از تفاوت بین گروهها میباشد. همچنین وجود تنوع بالا در درون گروهها نشان میدهد که در درون گروهها هم به اندازه کافی تنوع برای انتخاب وجود دارد. شاخص دیگری که جهت ارزیابی تنوع بین جمعیتها محاسبه گردید شاخص PhiPT بود. این شاخص نشان داد که تنوع بین گروهها از لحاظ نشانگرهای ISSR و RAPD در سطح یک درصد معنیدار میباشد (**162/0) و تأیید در جهت وجود تنوع ژنتیکی معنیدار در بین ارقام گندم مورد مطالعه است.
نتیجهگیری کلی
در این تحقیق تعداد 92 باند چندشکل برای 10 آغازگر نشانگر ISSR و همچنین 19 باند چندشکل برای چهار آغازگر نشانگر RAPD حاصل شد. همچنین ارقام براساس تجزیه خوشهای، تجزیه به مختصات اصلی و تجزیه کلاستر به روش بیزی با نرمافزار STRUCTURE در دو گروه قرار گرفتند و نتایج هم را تأیید نمودند. تجزیه واریانس مولکولی به خوبی تنوع بین گروهها و درون گروهها را نشان داد. بنابراین با توجه به نتایج، بین ارقام گندم نان از لحاظ نشانگرهای ISSR و RAPD تنوع ژنتیکی بسیار خوبی براساس ماتریس شباهت، همچنین نتایج گروهبندی و تجزیه واریانس مولکولی وجود دارد و این امید را میدهد که بتوان با تلاقی ارقام که شباهت کمی باهم دارند و انتخاب در نسلهای در حال تفرق و انتخاب روشهای صحیح اصلاحی باعث انتخاب لاینهای مناسب در جهت بهبود کمیت و کیفیت گندم شد. به منظور ارزیابی دقیق، مطالعات مقایسهای در سطح صفات مورفولوژیکی و مولکولی با نشانگرهای دیگری مانند SSR و همچنین تعداد بیشتری آغازگر با پراکنش کروموزومی بالا پیشنهاد میشود.
تشکر و قدردانی
بدین وسیله از دانشگاه تحصیلات تکمیلی صنعتی و فنآوری پیشرفته به خاطر تأمین هزینه مالی و همچنین امکانات آزمایشگاهی تشکر و قدردانی میشود.