Document Type : Research Paper
Authors
1 Department of Crop Biotechnology and Breeding, Faculty of Agriculture, Ferdowsi University of Mashhad, Mashhad, Iran
2 Department of Agrotechnology, Faculty of Agricultureو Ferdowsi University of Mashhad, Iran
3 Ph.D Student of Agrotechnology of Ferdowsi University of Mashhad, Iran
Abstract
Keywords
Main Subjects
مقدمه
با افزایش جمعیت جهان نیاز به مواد غذایی ازجمله دانههای روغنی رو به افزایش است (Desire et al., 2021). دانههای روغنی پس از غلات دومین ذخایر غذایی جهان را تشکیل میدهند و از نظر اسیدهای چرب و پروتئین از جایگاه ویژهای برخوردار هستند (Zanetti et al., 2013). در ایران مصرف روغنهای گیاهی باتوجهبه افزایش جمعیت کشور و تغییر الگوی غذایی مردم در حال افزایش است. از طرفی نیاز به واردات بخش اعظم روغن خوراکی و تداوم این روند، ضرورت سنجش امکانات و پتانسیلها و برنامهریزی جامع برای توسعه کشت دانههای روغنی را بیشازپیش مورد تأکید قرار داده است (Helali, 2018). علاوه بر این، کنجاله کنجد (Sesamum indicum L.) بهعنوان یک منبع پروتئینی در تغذیه دام حایز اهمیت میباشد. کنجاله کنجد تقریباً حاوی 42 درصد پروتئین خام است و میتواند بهعنوان جایگزین مناسب برای کنجاله سویا (Glycine max L.) در خوراک دام در نظر گرفته شود (Shahizad et al., 2019).
کنجد از خانواده Pedaliaceae بهدلیل سهولت استخراج، پایداری روغن و تحمل بالا به خشکی و گرما بهعنوان کشت تابستانه اهمیت زیادی در توسعه کشاورزی مناطق خشک و نیمهخشک دارد (Masoudi & Ahmadi, 2019 Yadav et al., 2022; Rezvani Moghadam et al., 2010; Weiss, 2000). وجود اسیدهای چرب اشباعنشده بالا، حضور برخی از آنتیاکسیدانها و میزان پایین کلسترول، کیفیتی مطلوب به روغن کنجد بخشیده است ( .(Gholinezhad & Darvishzadeh, 2018این موارد در کنار مقادیر بالای پروتئین، ویتامینها و لیگنانهای مختلف منجر شده است تا گیاه کنجد نقش بسیار مهمی در حفظ امنیت غذایی دارا باشد (Li et al., 2018) .طبـق آمار فائو در سال 2020 میزان سطح زیر کشت کنجد در ایران 42 هزار هکتار، میانگین عملکرد دانه 690 کیلوگرم در هکتـار و میزان تولید دانه کنجد 29 هزار تن بوده برآورد شده است (FAO, 2021) که این وضعیت میتواند با انتخاب ارقام باکیفیت و پتانسیل سازگاری بالا به شرایط آبوهوایی مختلف بهبـود یابـد (Yol & Uzun, 2012). عدم بهکارگیری روشهای نوین کشت و نیز فقدان مواد ژنتیکی اصلاحشده (ارقام جدیـد) سـازگار بـه منـاطق مختلف کشور، منجر به کاهش عملکرد شده است. از طرف دیگر، عمده زراعت کنجد در ایران با استفاده از تودههای بومی صورت میگیرد (Rmazani & Mansouri, 2017; Dossa et al., 2017) که از معایب آن میتوان به دیررسـی، غیر یکنـواختی و حساسیت به بیماریها اشاره کرد (Yadav et al., 2022; Zeinalzadeh et al., 2021).
در این راستا، استفاده از ارقام اصلاحشده میتواند موجب افزایش کمی و کیفی روغن در گیاه کنجد شود و موفقیت در یک برنامه بهنژادی بستگی به تنوع ژنتیکی و انتخاب مؤثر ژنوتیپهای برتر دارد (Masoudi & Ahmadi, 2019). گزینش ژنوتیپها توسط روشهای آماری چندمتغیره همچون همبستگی و تجزیه علیت، رگرسیون چندمتغیره، تجزیه به عاملها و تجزیه خوشهای انجام میشود. بهطورکلی مطالعـات همبستگی و اسـتفاده از روشهـای آمـاری چندمتغیره ایـن امکان را فراهم میسازد تا صفات مهم و تعیینکننده عملکـرد و میزان سهم نسبی هریک بر عملکرد مشـخص شود
(Askari et al., 2016). نتایج پژوهشها نشان داده است که عملکرد دانه، زیستتوده، شاخص برداشت و تعداد کپسول در بوته از صفات مهم مرتبط با عملکرد کنجد محسوب میشوند (Ranjithkumar et al, 2022;
Roy et al., 2022). علاوهبراین، همبستگی و تجزیه مسیر نشان داده است که صفات تعداد دانه در کپسول و ارتفاع بوته اثر مستقیم مثبت و همبستگی مثبت معنیداری با عملکرد دانه دارند؛ بنابراین میتوان این صفات را بهعنوان معیاری برای بهبود عملکرد دانه در برنامههای اصلاحی کنجد در نظر گرفت (Kalaiyarasi et al., 2019). نتایج مطالعه 705 توده کنجد نشان داد تعداد کپسول در بوته، تعداد دانه در بوته و وزن هزار دانه بهمنظور افزایش عملکرد میتوانند بهطور مستقیم مورد گزینش قرار گیرند (Zhou et al., 2018). نتایج پژوهشهای پیشین در بررسـی سـازگاری و روابط صـفات کمـی در 14 لایـن برتـر کنجـد بههمراه رقم اولتـان و تـوده محلـی بیرجند بـا اسـتفاده از تجزیـه بـه مؤلفـههـای اصـلی و تجزیـه خوشهای نشان داد که صفات تعداد شاخه فرعی، وزن هزار دانـه و عملکـرد کـل بیشـترین وزن را در گروهبندی ارقام داشتند و بر اساس صفات ذکرشده برای کشت کنجد در منطقه بیرجند، توده محلی بیرجند توصیه شده است (Ramazani & Mansouri, 2017). در پژوهشی دیگر 65 ژنوتیپ کنجد ازنظر صفات مورفولوژیک و زراعی ازجمله ارتفاع بوته، تعداد انشعاب در بوته، تعداد کپسول در بوته، وزن دانه و عملکرد دانه در پنج گروه توسط تجزیه خوشهای گروهبندی شدند (Srikanth & Ghodke, 2022). شناسایی ژنوتیپهای برتر حاصل از آزمایشهای سازگاری میتواند منجر بـه معرفـی ارقام جدید و جایگزینی آنها با تودههای بومی فعلی شود. هدف از این پژوهش بررسی عملکرد و اجزای عملکرد ژنوتیپهای کنجد بر اساس گـزینش صـفات مطلوب در شـرایط آبوهوایی مشهد بود.
مواد و روشها
این مطالعه در قالب طرح بلوکهای کامل تصادفی با سه تکرار با 20 ژنوتیپ کنجد در مزرعه تحقیقاتی دانشکده کشاورزی دانشگاه فردوسی مشهد واقع در 10 کیلومتری جنوب شرق مشهد (عرض جغرافیایی 36 درجه و 15 دقیقه شمالی و طول جغرافیایی 56 درجه و 28 دقیقه شرقی و ارتفاع 985 متر از سطح دریا) در سال زراعی 1400 اجرا شد. قبل از کاشت، عملیات آمادهسازی زمین و خاکورزی انجام شد و هر کرت شامل چهار ردیف با فاصله 50 سانتیمتر بود و فاصله کرتها و بلوکها از یکدیگر یک و نیم متر در نظر گرفته شد. فاصله بوتهها روی ردیف 18 سانتیمتر با تراکم 111 هزار بوته در هکتار کشت شد. در این آزمایش آبیاری به روش جوی و پشتهای انجام شد. اولین آبیاری بلافاصله پس از کاشت و آبیاریهای بعدی با فواصل هر هفت روز یکبار انجام شد. کنترل علفهای هرز در طول فصل رشد دو بار بهصورت وجین دستی صورت گرفت و از هیچگونه کود و سم شیمیایی استفاده نشد. مشخصات آبوهوایی و ویژگیهای فیزیکوشیمیایی خاک محل آزمایش بهترتیب در شکل 1 و جدول 1 ارائه شده است. تاریخ کاشت دوم خردادماه و تاریخ برداشت 19 مهرماه 1400 (132 روز) بود.
شکل 1- شرایط آبوهوایی رشد کنجد از کاشت تا برداشت (محدوده بین پیکانها دوره رشد کنجد را نشان میدهد).
Figure 1. Weather conditions during the growing season, from the sowing date to harvest time (The range between the arrows shows the period of sesame growth).
جدول 1- خصوصیات فیزیکی و شیمیایی خاک مزرعه آزمایشی (0-30 سانتیمتر).
Table 1. Physical and chemical properties of experimental field soil (0-30 cm).
Soil texture |
Electrical conductivity (dS.m-1) |
Soil pH |
Mineral carbon (%) |
Nitrogen (%) |
Phosphorus (mg.Kg-1) |
Potassium (mg.Kg-1) |
Loam |
2.75 |
7.39 |
1.05 |
0.063 |
14 |
117 |
در پایان فصل رشد، پس از رسیدگی فیزیولوژیک، ارتفاع بوته، ارتفاع اولین کپسول، عملکرد دانه، زیستتوده و اجزای عملکرد شامل تعداد شاخه فرعی، تعداد کپسول در بوته و وزن هزار دانه در پنج بوته بهصورت تصادفی اندازهگیری شد.
تجزیه به عاملها با استفاده از روش تجزیه به مؤلفه اصلی برآورد شدند. سپس بهمنظور گروهبندی ژنوتیپهای مورد بررسی، تجزیه خوشهای مبتنی بر روش وارد (Ward) مورد استفاده قرار گرفت. بهمنظور تأیید صحت گروهبندی انجامشده، از تجزیه واریانس چندمتغیره و تجزیه تابع تشخیص استفاده شد. توابع تشخیص بهدستآمده از نوع توابع خطی فیشر بودند. گروهبنـدی ژنوتیپها بـا اسـتفاده از دو عامـل اصلی اول که بیشترین درصد تغییرات را توجیـه کردنـد، صـورت گرفت. برای بررسی تفاوت گروهها ازلحاظ صفات مختلف، مقایسه میانگین گروهها برای صفات مورد بررسی انجام شد. مقایسه میانگین دادهها با استفاده از آزمون دانکن در سطح احتمال پنج درصد انجام گرفت. همچنین، برای آزمون نرمالبودن دادهها و یکنواختی واریانسها از نرمافزار، Minitab16 و برای تجزیه واریانس چندمتغیره، تجزیه تابع تشخیص، تجزیه خوشهای و تجزیه به مؤلفههای اصلی و ترسیم نمودارهای دوبعدی از نرمافزارهای SPSS19 و STATISTICA8 و برای آنالیز همبستگی صفات از نرمافزار JMP4 استفاده شد.
نتایج و بحث
عملکرد و اجزای عملکرد
در آزمایش حاضر صفات مورد ارزیابی در ژنوتیپهای کنجد در سطح احتمال یک درصد معنیدار شدند (جدول 2).
جدول 2- تجزیه واریانس (میانگین مربعات) صفات مورد مطالعه در کنجد در شرایط آب و هوایی مشهد.
Table 2. Analysis of variance (mean of squares) of different sesame traits under Mashhad climatic conditions.
S.O.V |
df |
Plant height |
Lowest capsule height |
Branch no. |
Capsuleno. plant-1 |
Biomass plant-1 |
Grain weight plant-1 |
1000- grain weight |
Biological yield |
Grain yield |
Harvest index |
Block |
2 |
1321** |
140** |
2.02** |
469* |
817** |
4.08* |
0.22** |
9404199** |
244414** |
195** |
Genotypes |
19 |
133** |
195** |
2.34** |
1580** |
279** |
60.68** |
0.72** |
3316189** |
715835** |
166** |
Error |
38 |
45.05 |
14.50 |
0.24 |
193.3 |
60.52 |
1.91 |
0.072 |
818232 |
23192 |
21.75 |
CV (%) |
-- |
7.97 |
10.04 |
14.03 |
14.03 |
18.86 |
29.84 |
8.85 |
20.07 |
29.91 |
44.41 |
** و * بهترتیب بیانگر اختلاف معنیدار در سطوح احتمال یک درصد و پنج درصد میباشد.
* and ** are significant at 5% and 1% probability levels, respectively.
نتایج حاصل از مقایسه میانگین ژنوتیپها (جدول 3) نشان داد که ژنوتیپ داورزن، افغان1 و بشرویه دارای بیشترین ارتفاع بوته بودند که با هم تفاوت معنیداری نشان ندادند و کمترین ارتفاع بوته مربوط به ژنوتیپ طبس بود. ژنوتیپ کهنوج از ارتفاع اولین کپسول بیشتری نسبت به سایر ژنوتیپها برخوردار بود. ارتفاع بوته مناسب، بهمنظور برداشت مکانیزه کنجد مؤثر میباشد؛ لذا ارتفاع بوته در کنجد حائز اهمیت بوده و وجود تنوع ژنتیکی برای این صفت، امکان انتخاب ژنوتیپهای مناسب ازلحاظ ارتفاع را فراهم میآورد (Takele et al., 2022).
بیشترین تعداد شاخه فرعی مربوط به ژنوتیپ افغان1 بود و کمترین مقدار این صفت مربوط به ژنوتیپ داورزن با میانگین 93/1 بود. تعداد شاخه در بوته در برنامههای اصلاحی کنجد، صفتی مهم محسوب میشود؛ زیرا انواع تک ساقه زودرستر هستند و معمولاً یکنواختی رسیدگی بیشتری نسبت به انواع منشعب دارند
(Khajehpour, 2013). ژنوتیپ داورزن دارای بیشترین تعداد کپسول در بوته بود و پس از آن ژنوتیپ افغان1 بیشترین تعداد کپسول در بوته را به خود اختصاص داد.
ژنوتیپ بشرویه بیشترین وزن هزار دانه و پسازآن ژنوتیپ افغان1 و افغان5 بیشترین مقدار این صفت را به خود اختصاص دادند. ازنظر صفت وزن خشک اندام هوایی تکبوته و زیستتوده افغان1 بیشترین مقدار را به خود اختصاص داد و بهترتیب ژنوتیپهای کهنوج و وحدتی2 کمترین مقدار زیستتوده را دارا بودند.
عملکرد دانه تنوع گستردهای را نشان داد؛ بهطوریکه ژنوتیپ افغان1 با تولید 2109 کیلوگرم در هکتار و ژنوتیپ کهنوج با عملکرد 46 کیلوگرم در هکتار بهترتیب بیشترین و کمترین عملکرد دانه را بهخود اختصاص دادند. همچنین ژنوتیپ افغان1 از بیشترین شاخص برداشت برخوردار بود و کمترین مقدار این صفت مربوط به ژنوتیپ کهنوج بود (جدول 3). ژنوتیپ افغان1 با شاخص برداشت بـالا، میـزان بیـشتری از مـاده خشک را به دانهها اختصاص داده است. بـه عبـارت دیگـر بـا توزیع بیشتر مواد فتوسنتزی به مخازن، قـسمت زیـادی از ماده خشک به عملکرد دانـه تخصیص یافته است.
عملکرد دانه در گیاهان زراعی یک ویژگی کمّی و پیچیده اسـت که تحت تأثیر اجزای عملکـرد و خصوصیات زراعی است؛ لذا آگاهی از نقـش و تـأثیر ایـن اجـزا و صـفات زراعی موثر بر عملکرد و شناخت روابط بین آنها برای عملیات اصلاحی و بهنژادی ضروری است
(Koochaki & Banayan, 1994; Bayat & Vaezi, 2016). تولید شاخ و برگ بیشتر در گیاه موجب افزایش جـذب نـور و فتوسنتز در گیاه شده که این امر موجب میشود تا گیاه در تولید عملکرد نهایی و افزایش وزن دانهها نیز موفق عمل کنـد (Torabi et al., 2007). بهطوریکه ژنوتیپ افغان1 با بیشترین تعداد شاخه فرعی بیشترین وزن دانه و عملکرد دانه را به خود اختصاص داد. دراینرابطه، نتایج پژوهشهای پیشین نشان داد ژنوتیپهایی با تعداد شاخه بارور، تعداد کپسول در بوته، تعداد دانه در کپسول و عملکرد دانه بیشتر میتواند بهعنوان ژنوتیپهای مطلوب معرفی شوند
(Vivek et al., 2022). همچنین در مطالعه دیگری صفات تعداد کپسول در بوته، تعداد دانه در کپسول، ارتفاع بوته و عملکرد دانه را بهعنوان معیاری برای انتخاب ژنوتیپ برتر در نظر گرفتند (Kalaiyarasi et al., 2019).
جدول 3- مقایسه میانگین صفات مورد مطالعه در ژنوتیپهای کنجد در شرایط آب و هوایی مشهد.
Table 3. Mean comparison of different genotypes of sesame under Mashhad climatic conditions.
Genotypes |
Plant height (cm) |
Lowest capsule height (cm) |
Branch no. |
Capsule no. plant-1 |
Biomass (g.plant-1) |
Grain weight (g.plant-1) |
1000- grain weight (g) |
Biological yield (kg.ha-1) |
|
Grain yield (kg.ha-1) |
Harvest index (%) |
||||||
Afghan 1 |
95.0ab |
36.4bc |
4.73a |
91.8ab |
66.5a |
19.17a |
3.56a-c |
7313a |
|
2109a |
30.3a |
||||||
Afghan 2 |
80.9a-c |
28.5c |
2.20de |
69.3a-g |
38.0bc |
2.69de |
2.83c-g |
4180b |
|
296de |
8.31cd |
||||||
Afghan 3 |
87.7a-c |
35.9bc |
4.60a |
83.3a-d |
51.3a-c |
8.18bc |
3.10b-g |
5646ab |
|
899bc |
15.9a-c |
||||||
Afghan 4 |
85.9a-c |
43.3b |
4.06ab |
78.1a-f |
48.7a-c |
3.17de |
3.04c-g |
5359ab |
|
350de |
6.47cd |
||||||
Afghan 5 |
85.9a-c |
35.5bc |
2.40de |
83.6a-d |
54.7ab |
9.02b |
3.90ab |
6017ab |
|
991b |
16.6a-c |
||||||
Afghan 6 |
83.6a-c |
32.9bc |
2.46de |
56.6a-g |
30.0c |
1.24de |
2.26g |
3300b |
|
137de |
4.13cd |
||||||
Afghan 7 |
76.7bc |
39.8bc |
3.67a-d |
54.9a-h |
34.0bc |
1.57de |
2.80c-g |
3740b |
|
173de |
4.61cd |
||||||
Afghan 8 |
79.3b-c |
36.3bc |
3.60a-d |
89.9a-c |
52.8a-c |
8.99b |
3.00c-g |
5809ab |
|
989b |
17.2a-c |
||||||
Boshrooyeh |
91.5a-c |
38.1bc |
4.56a |
46.4d-h |
34.9bc |
2.99de |
4.20a |
3839b |
|
328de |
9.48cd |
||||||
Davarzan |
101a |
44.6b |
1.93e |
97.8a |
44.0a-c |
3.29de |
3.10b-g |
4817ab |
|
362de |
7.83cd |
||||||
Eshghabad |
81.3a-c |
38.5bc |
4.20ab |
49.9c-h |
35.3bc |
2.82de |
3.40a-d |
3886b |
|
310de |
9.07cd |
||||||
Feizabad |
80.1b-c |
32.9bc |
4.20ab |
40.8e-h |
34.1bc |
2.05de |
3.00c-g |
3748b |
|
225de |
6.24cd |
||||||
Gorgan |
87.1a-c |
29.5c |
2.53de |
61.7a-g |
47.1a-c |
9.53b |
3.43a-d |
4600ab |
|
1048b |
24.9ab |
||||||
Kahnooj |
90.5a-c |
66.7a |
2.73b-e |
12.3h |
31.3bc |
0.42e |
2.45fg |
3446b |
|
46.0e |
1.29d |
||||||
Kalat |
88.1a-c |
33.6bc |
3.93a-d |
51.6a-h |
41.3bc |
1.47de |
3.16b-f |
4546ab |
|
161de |
3.20cd |
||||||
Sarayan |
79.8bc |
35.1bc |
4.46a |
37.4f-h |
40.7bc |
2.88de |
2.48fg |
4473b |
|
316de |
7.45cd |
||||||
Tabas |
74.1c |
33.0bc |
3.53a-d |
41.9e-h |
32.0bc |
1.88de |
2.55fg |
3520b |
|
206de |
5.82cd |
||||||
Torbat-e-Jam |
80.9a-c |
33.5bc |
3.72a-d |
42.8e-h |
40.0bc |
5.40bd |
3.13b-f |
4400b |
|
594bd |
14.5bc |
||||||
Vahdati1 |
76.5bc |
39.7bc |
4.00a-c |
54.0a-h |
36.7bc |
4.02a-e |
2.66d-g |
4033b |
|
442c-d |
11.7bd |
||||||
Vahdati 2 |
83.8a-c |
44.1b |
2.80a-d |
33.9g-h |
31.3bc |
1.79de |
2.59e-g |
3446b |
|
197de |
4.86cd |
||||||
مقادیر دارای حروف مشترک در هر ستون فاقد تفاوت معنیدار آماری در سطح احتمال پنج درصد بر اساس آزمون آماری دانکن میباشند.
Values in each column followed by similar letters are not significantly different at 5% probability level, using Duncan’s multiple range test, No: Number.
تعیین ضرایب همبستگی صفات
شناخت رابطه بین عملکرد دانه و صفات مورفولوژیک در انتخاب ژنوتیپ برتر، از اهمیت زیادی برخوردار است. بر اساس میانگین ژنوتیپها ضرایب همبستگی بین صفات مورد بررسی محاسبه شدند. نتایج نشان داد که عملکرد دانه بهترتیب با تعداد شاخه فرعی (**53/0r=)، تعداد کپسول در دانه (**58/0r=) و وزن هزار دانه (**42/0r=) همبستگی مثبت و معنیداری دارد. بیشترین ضرایب همبستگی بهترتیب بین زیستتوده با وزن خشک اندام هوایی تکبوته (**97/0r=)، عملکرد دانه با وزن دانه تکبوته (**90/0r=)، شاخص برداشت با وزن دانه تکبوته (**89/0r=) و شاخص برداشت با عملکرد دانه (**89/0r=) مشاهده شد (جدول 4). در بررسـی انجـامشـده، روی 12 ژنوتیپ کنجد هـمبسـتگی مثبت و معنیداری بـین تعداد شاخه فرعی، وزن کپسـول در بوته و وزن هزار دانه مشاهده شد (Askari et al., 2016; Zhou et al., 2018). در دیگر پژوهشها هـمبسـتگی بین تعداد شاخه فرعی، تعداد کپسول در بوته و وزن هزار دانه با عملکرد مثبـت و معنیدار گزارش شده است (Patil & Lokesha, 2018; Kumar et al., 2022). همبستگی این صفات با عملکرد دانه میتواند بهدلیل وجود پیوستگی قوی بین ژنهای کنترلکننده این صفات با عملکرد دانه باشد
(Tahmasebi et al., 2021). بر مبنای مطالعات یادشده عوامـل مهـم مؤثر در عملکرد کنجـد شـامل وزن هـزار دانـه، تعـداد شـاخه فرعی و تعـداد کپسـول در بوتـه هستند و لـذا مـیتواننـد معیــارهــای مناســبی بــرای انتخاب ژنوتیپ برتر محسوب شوند. در آزمایش حاضر، همبستگی مثبت و معنیدار عملکرد دانه با زیستتوده نشان میدهد که فرایند انتقال و تسهیم فراوردههای فتوسنتزی به مخزن بهطور مستقیم با منبع در ارتباط است. بنابراین توانایی بالا در تولید زیستتوده میتواند به افزایش تولید کپسول در بوته و درنهایت افزایش عملکرد دانه منجر شود.
جدول 4- ضرایب همبستگی بین عملکرد و اجزای عملکرد ژنوتیپهای کنجد در شرایط آب و هوایی مشهد.
Table 4. Correlation coefficient between yield and yield components of sesame genotypes under Mashhad climatic conditions.
Traits |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
1. Plant height |
1 |
|
|
|
|
|
|
|
|
|
2. Lowest capsule height |
0.30** |
1 |
|
|
|
|
|
|
|
|
3. Branch no. |
-0.04ns |
-0.17ns |
1 |
|
|
|
|
|
|
|
4. Capsule no. plant-1 |
-0.24ns |
0.27* |
0.03ns |
1 |
|
|
|
|
|
|
5. Biomass plant-1 |
-0.01ns |
0.43** |
0.24ns |
0.69** |
1 |
|
|
|
|
|
6. Grain plant-1 |
-0.19ns |
0.12ns |
0.25ns |
0.58** |
0.67** |
1 |
|
|
|
|
7. 1000- grain weight |
-0.25ns |
0.29* |
0.14ns |
0.28* |
0.35** |
0.42** |
1 |
|
|
|
8. Biological yield |
-0.01ns |
0.43** |
0.25ns |
0.69** |
0.97** |
0.63** |
0.33** |
1 |
|
|
9. Grain yield |
-0.19ns |
0.12ns |
0.53** |
0.58** |
0.67** |
0.90** |
0.42** |
0.63** |
1 |
|
10. Harvest index |
-0.26* |
-0.03ns |
0.18ns |
0.40** |
0.37** |
0.89** |
0.41** |
0.30** |
0.89** |
1 |
ns، * و ** بهترتیب بیانگر عدم اختلاف معنیدار و اختلاف معنیدار در سطوح احتمال پنج درصد و یک درصد میباشد.
ns, * and ** are non-significant and significant at 5% and 1% probability levels
تجزیه به عاملها
روش تجزیه به عاملها یکی از روشهای آماری چندمتغیره است که در آن میتوان تعداد زیادی از متغیرهای همبسته را به تعداد کمتری عامل اصلی کاهش داد. پارامترهای حاصل از تجزیه به مؤلفههای اصلی شامل مقادیر ویژه، درصد واریانس توجیهشده، درصد واریانس تجمعی در جدول 5 نمایش داده شده است. در این مطالعه تجزیه به عاملها با استفاده از دادههای استانداردشده نشان داد که هر سه عامل دارای مقادیر ویژه بالاتر از یک بود و درمجموع 1/83 درصد تغییرات کل را شامل شدند (جدول 5). عامل اول 7/57 درصد از تغییرات را توجیه کرد که شامل تعداد کپسول در بوته، وزن خشک اندام هوایی، وزن دانه تکبوته، وزن هزار دانه، زیستتوده، عملکرد دانه و شاخص برداشت با بار منفی بود. عامل اول نشاندهنده ژنوتیپهایی با صفات تعداد کپسول در بوته، وزن خشک اندام هوایی در تکبوته، وزن دانه تکبوته، وزن هزار دانه، زیستتوده، عملکرد دانه و شاخص برداشت پایین است. این عامل نشاندهنده ژنوتیپهایی با تعداد کپسول در بوته، وزن خشک اندام هوایی در تک بوته، وزن دانه تکبوته، وزن هزار دانه، زیستتوده، عملکرد دانه و شاخص برداشت کمتری هستند. عامل دوم حدود 2/15 درصد تغییرات را توجیه کرد که شامل ارتفاع بوته و ارتفاع اولین کپسول با بار منفی بود. عامل سوم 2/10 درصد از تغییرات را توجیه کرد که تعداد شاخه فرعی با بار مثبت بیشترین تأثیر را در این عامل داشت؛ این عامل نشاندهنده ژنوتیپهایی با تعداد شاخه فرعی بیشتری است (جدول 5). ازآنجاییکه عملکرد دانه مهمترین ویژگی در بهگزینی ژنوتیپها است، گزینش بر اساس عامل اول بهعلت ضریب بیشتر این ویژگی، به غربال ژنوتیپهایی با عملکرد دانه بالا منجر خواهد شد.
جدول 5- ضرایب حاصل از تجزیه به عاملها در ژنوتیپهای کنجد.
Table 5. Factor analysis for sesame genotypes.
Factor 1 |
Factor 2 |
Factor 3 |
|
Plant height |
-0.391 |
-0.824 |
0.114 |
Lowest capsule height |
0.306 |
-0.678 |
0.512 |
Branch no. |
-0.229 |
0.559 |
0.738 |
Capsule no. plant-1 |
-0.779 |
-0.128 |
-0.382 |
Biomass plant-1 |
-0.960 |
-0.061 |
0.002 |
Grain plant-1 |
-0.957 |
0.066 |
0.051 |
1000- grain weight |
-0.626 |
-0.102 |
0.207 |
Biological yield |
-0.942 |
-0.060 |
0.034 |
Grain yield |
-0.957 |
0.066 |
0.051 |
Harvest index |
-0.913 |
0.153 |
-0.022 |
Eigen value |
5.77 |
1.52 |
1.02 |
Cumulative variance % |
57.7 |
72.9 |
83.1 |
باتوجهبهاینکه دو عامل اصلی اول و دوم بهویژه عامل اول بیشترین تغییرات واریانس دادهها را توجیه کردند از صفات تعداد کپسول در بوته، وزن خشک اندام هوایی در تکبوته، وزن دانه تکبوته، وزن هزار دانه، زیستتوده، عملکرد دانه و شاخص برداشت برای بهدستآوردن پراکنش و شناسایی ژنوتیپهای برتر در دستگاه مختصات استفاده شد (شکل 2). همانطور که ملاحظه میشود ژنوتیپهای طبس، افغان2، سرایان، فیضآباد، افغان6، وحدتی1، افغان8، عشقآباد ازنظر عامل اول از ضرایب پایینی برخوردارند و بهعنوان ژنوتیپهای با عملکرد پایین معرفی میشوند (شکل 2). علاوهبراین، نمودار دوبعدی نشان داد که ژنوتیپ افغان1 از ضرایب بالایی برای عامل اول برخوردار است و بهصورت ژنوتیپی با عملکرد و اجزای عملکرد بالا شناسایی شد (شکل 2). همچنین در بین ژنوتیپهای مورد مطالعه، ژنوتیپ کهنوج با کمترین صفات از نظر تعداد کپسول در بوته، وزن هزار دانه و عملکرد زیستتوده به عنوان ژنوتیپی با کمترین عملکرد دانه معرفی شد (شکل 2).
تجزیه خوشهای
بهمنظور تعیین قرابت ژنوتیپها و گروهبندی آنها برمبنای صفات مورد بررسی، از تجزیه خوشهای بهروش Ward و استفاده از فاصله ماهالانوبیس بهمنظور تعیین خویشاوندی ژنوتیپها و گروهبندی آنها بر اساس صفات مهم زراعی استفاده شد. برایناساس، ژنوتیپهای مورد بررسی در سه گروه تقسیمبندی شدند که در گروه اول هشت ژنوتیپ، در گروه دوم هفت ژنوتیپ و در گروه سوم پنج ژنوتیپ قرار گرفتند که خوشه سوم از میانگینهای بالاتری نسبت به سایر گروهها ازنظر عملکرد و اجزای عملکرد برخوردار بود (شکل 3).
شکل 2- پراکنش ژنوتیپهای کنجد بر اساس عوامل اول و دوم.
Figure 2. The distribution of different genotypes of sesame based on the operating results of the factor analysis.
بهمنظور تأیید تفاوت بین گروهها، تجزیه واریانس چندمتغیره بر پایه طرح کاملاً تصادفی نامتعادل برای صفات مورد نظر انجام شد که در آن، آماره ویلکس لامبدا (Wilks' Lambda) در سطح احتمال یک درصد و پنج درصد معنیدار بود (جدول 6)؛ بنابراین، میتوان نتیجه گرفت بین بردار میانگینها اختلاف معنیداری وجود دارد. بهاینترتیب، ژنوتیپهای قرارگرفته درون گروهها نسبت به ژنوتیپهای قرارگرفته در گروههای متفاوت ازنظر صفات مورد بررسی، شباهت بیشتری با هم داشته و گروهبندی، بهطور صحیح انجام شده است.
تجزیه تابع تشخیص
بررسی صحت گروهبندیهای بهدستآمده از روش تجزیه خوشهای، از طریق تابع تشخیص انجام شد (جدول 7). نتایج تجزیه تابع تشخیص نشان داد که تمامی ژنوتیپها بهطور صحیح گروهبندی شدهاند و میزان موفقیت تابع تشخیص برای تمام گروهها 100 درصد است که این مقدار را میزان موفقیت تابع تشخیص میگویند. میزان موفقیت نشان میدهد که تابع تشخیص تا چه حد در گروهبندی یا تشخیص بین گروهها موفق بوده است. در این رابطه، پژوهشگران 94 ژنوتیپ کنجد را ازنظر پنج صفت توسط تجزیه خوشهای به پنج گروه دستهبندی کردند و با تابع تشخیص نشان دادند که تمامی ژنوتیپهای مورد بررسی بهطور صحیحی گروهبندی شدند و میزان موفقیت تابع تشخیص برای تمام گروهها 100 درصد بود
(Masoudi & Ahmadi, 2019). همچنین مطالعه 20 ژنوتیپ یونجه (Medicago sativa L.) ازلحاظ صفات زراعی و مورفولوژیک، در سه گروه دستهبندی شد و نتایج تابع تشخیص نشان داد 100 درصد ژنوتیپها به گروه خود تعلق داشتند (Khodarahmpour & Motamedi, 2016).
شکل 3- دندوگرام مربوط به گروهبندی ژنوتیپهای کنجد با استفاده از صفات مورد مطالعه.
Figure 3. Dendrogram of clustering of studied sesame genotypes based on traits.
جدول 6– تجزیه واریانس چندمتغیره بر پایه طرح کاملاً تصادفی نامتعادل، آمارۀ ویلکس لامبدا در ژنوتیپهای کنجد.
Table 6. Analysis variation of multi variables based on unbalanced completely randomized design (CRD) Wilks' Lambda in sesame genotypes.
Function |
Df |
Wilks' Lambda |
Chi-square |
Probability level |
1 |
18 |
0.010 |
59.42 |
0.000 |
2 |
8 |
0.296 |
15.81 |
0.045 |
در تجزیه تابع تشخیص کانونیکی، متغیر اول کانونیک که مقادیر ویژه بالاتر از یک داشت، درمجموع 92 درصد واریانس موجود را تبیین کردند که میتواند بهعنوان معیاری مطمئن بهمنظور انتساب ژنوتیپ جدید به گروه صحیح مورد استفاده قرار گیرد (جدول 8). همبستگی کانونیکی معنیدار بین ژنوتیپها با اولین متغیر کانونیک (**98/0=R) و دومین متغیر کانونیک (**83/0=R) نشاندهنده این است که متغیرهای کانونیک تفاوت بین
ژنوتیپها را بهخوبی توجیه میکنند (جدول 8).
جدول 7- نتایج تابع تشخیص برای صحت گروهبندی ژنوتیپهای کنجد
Table 7. The results of discriminant function for clustering validity of sesame genotypes.
Group |
Group Membership |
||||
1 |
2 |
3 |
Total |
||
Total |
1 |
8 |
0 |
0 |
8 |
2 |
0 |
7 |
0 |
7 |
|
3 |
0 |
0 |
5 |
5 |
|
Percent |
1 |
100 |
0 |
0 |
100 |
2 |
0 |
100 |
0 |
100 |
|
3 |
0 |
0 |
100 |
100 |
100 درصد گروهها بهطور صحیح گروهبندی شدهاند. 100.0% of original grouped cases correctly classified.
ضرایب تشخیص استانداردشده کانونیکی همبستگی خطی ساده بین متغیرهای اصلی و متغیرهای کانونیکی را محاسبه میکند؛ لذا ضرایب تشخیص استانداردشده کانونیکی، منعکسکننده واریانس مشترکی است که متغیرهای کانونیک دارند و میتواند در ارزیابی توجیه نسبی هر متغیر در هر مرحله کانونیک مورد تفسیر قرار گیرد (Rencher, 2002). درواقع در تجزیه تشخیص کانونیکی، تفاوت گروهها بر اساس همبستگی میان متغیرهای مستقل (صفات اندازهگیریشده) و ارتباط آنها با متغیرهای وابسته (ژنوتیپها) میباشد
(Vaylay & Van Santen, 2002).
بر اساس ضرایب استانداردشده کانونیکی، صفات زیستتوده، شاخص برداشت، وزن خشک اندام هوایی در تکبوته و تعداد کپسول در بوته در اولین معادله تشخیص کانونیکی بالا بود. علاوه بر این صفات وزن دانه تکبوته و وزن هزار دانه در دومین معادله تشخیص کانونیکی قابل توجه بود (جدول 8). بر مبنای نتایج بهدستآمده این صفات بیشترین تأثیر را در تنوع بین ژنوتیپها مورد مطالعه دارند و به عبارتی سهم بیشتری از تنوع بین ژنوتیپها بر اساس این صفات قابل توجیه است.
در این بررسی، از متغیرهای کانونیکی معنیدار اول و دوم برای گروهبندی ژنوتیپها استفاده شد (شکل 4). بر این اساس سه گروه کاملاً مجزا بهدست آمد و همانطور که مشاهده میشود بیشترین فاصله بین گروه اول و سوم و کمترین فاصله بین گروه دوم و سوم مشاهده شد (شکل 4). در برنامههای اصلاحی، انتخاب والدین مناسب برای برنامههای هیبریداسیون بسیار پراهمیت است که بهرهبرداری از حداکثر تنوع ژنتیکی را تسهیل میکند. تلاقی ژنوتیپهایی که در یک گروه قرار دارند، نمیتواند پاسخگوی انتظارات اصلاحگران برای بهگزینی آنها باشد (Bertan et al., 2007)؛ بنابراین ژنوتیپهای مربوط به گروه دوم باوجود عملکرد بالا بهدلیل فاصله ژنتیکی کم، بهتر است با گروه سوم تلاقی داده نشود. در برنامههای بهنژادی بهمنظور ایجاد رقم هیبرید برای حصول حداکثر هتروزیس بهتر است از افراد مربوط به جمعیتهایی که در گروههای متفاوت و دارای فاصله بیشتر (گروه 1 و 3 در این پژوهش) هستند استفاده شود.
جدول 8- ضرایب استاندارد کانونیکی صفات مورد مطالعه در ژنوتیپهای کنجد در شرایط آب و هوایی مشهد.
Table 8. Standardized canonical discriminant function coefficients measured groups in sesame genotypes under Mashhad climatic conditions.
Traits |
Canonical Function |
|
1 |
2 |
|
Plant height |
0.002* |
-1.79 |
Lowest capsule height |
0.067 |
1.19* |
Branch no. |
0.507 |
0.33* |
Capsule no. plant-1 |
0.449* |
1.09 |
Biomass plant-1 |
1.91* |
2.88 |
Grain plant-1 |
8.43 |
6.73* |
1000- grain weight |
1.43 |
1.55* |
Biological yield |
-5.63* |
-3.79 |
Harvest index |
-6.79* |
-6.45 |
Eigenvalue |
27.6 |
2.38 |
Cumulative % |
92.1 |
100 |
Canonical correlation |
0.982** |
0.839** |
*: بالاترین همبستگی مشاهدهشده بین هر صفت و متغیر کانونیکی
*: Largest absolute correlation between each variable and any discriminant function
شکل 4- گروهبندی ژنوتیپهای کنجد بر اساس متغیرهای کانونیک معنیدار
Figure.4. Cluster grouping of sesame genotypes based on significant canonical variable
بهمنظور مطالعه دقیقتر گروهها برای هر یک از صفات مورد مطالعه بهصورت جداگانه تجزیه واریانس یکطرفه انجام شد (جدول 9). نتایج نشان داد صفات تعداد کپسول در بوته، وزن خشک اندام هوایی تکبوته، وزن دانه تکبوته، عملکرد دانه، زیستتوده و شاخص برداشت بین گروهها تفاوت معنیداری داشتند و در سایر صفات تفاوت
معنیداری مشاهده نشد (جدول 9).
جدول 9- تجزیه واریانس (میانگین مربعات) گروهها بر اساس صفات مورد مطالعه کنجد در شرایط آب و هوایی مشهد
Table 9. Analysis of variance (mean square) based on measured groups in sesame genotypes under Mashhad climatic conditions.
Traits |
Between Groups |
Within Groups |
Df |
2 |
17 |
Plant height |
17.95ns |
47.75ns |
Lowest capsule height |
62.97ns |
65.46ns |
Branch no. |
0.678ns |
0.796ns |
Capsule no. plant-1 |
3013** |
234** |
Biomass plant-1 |
739** |
17.19** |
Grain plant-1 |
96.14** |
10.73** |
1000- grain weight |
0.282ns |
0.237ns |
Biological yield |
8985754** |
178295** |
Grain yield |
1163490** |
129796** |
Harvest index |
210* |
37.37* |
ns، * و ** بهترتیب بیانگر عدم اختلاف معنیدار و اختلاف معنیدار در سطوح احتمال پنج درصد و یک درصد میباشد.
ns, * and ** are non-significant and significant at 5% and 1% probability levels
مقایسه میانگین گروهها با استفاده از آزمون چنددامنهای دانکن برای ارتفاع بوته نشان داد که گروه سوم بیشترین میانگین را دارا بود و با گروه اول و دوم از نظر آماری تفاوت معنیداری نداشت (جدول 10). ازنظر ارتفاع اولین کپسول گروه اول در بالاترین رتبه قرار گرفت و با گروه اول و دوم از نظر آماری تفاوت معنیداری نشان نداد. تعداد شاخه فرعی در دو گروه اول و سوم بیشترین مقدار را به خود اختصاص دادند (جدول 10). ازنظر تعداد کپسول در بوته گروه سوم بیشترین مقدار را به خود اختصاص داد که تفاوت میانگین این گروه با میانگین کل 45 درصد بود و با گروه اول و سوم از نظر آماری تفاوت معنیداری را نشان داد (جدول 10).
صفات وزن خشک اندام هوایی در تک بوته، وزن دانه تک بوته و وزن هزار دانه در گروه سوم بیشترین مقدار را به خود اختصاص داد که نسبت به میانگین کل بهترتیب 32، 109 و 9 درصد برتری داشت. ازنظر صفات زیستتوده و عملکرد دانه گروه سوم در بالاترین رتبه قرار گرفت که نسبت به میانگین کل به ترتیب 109 و 32 درصد برتری داشت. صفت شاخص برداشت در گروه سوم بیشترین مقدار را به خود اختصاص داد که نسبت به میانگین کل 64 درصد برتری داشت و با گروه دوم در یک گروه آماری یکسان قرار گرفت (جدول 10). گروه سوم ازنظر صفات ارتفاع، تعداد شاخه فرعی، تعداد کپسول در بوته، وزن خشک اندام هوایی در تکبوته، وزن دانه تکبوته، وزن هزار دانه، زیستتوده، عملکرد دانه و شاخص برداشت بیشترین مقدار را به خود اختصاص داد و از میانگین بالاتری نسبت به میانگین کل برخوردار بود. بر اساس نتایج مقایسه میانگین صفات برای گروهها (جدول 10) ژنوتیپهای گروه سوم ازنظر اکثر صفات زراعی مورد بررسی در مقایسه با گروه اول و دوم، از مقادیر بالاتری برخوردار بودند. بدیهی است اگر میانگین یک صفت در یک گروه، از میانگین آن صفت در سایر گروهها و همچنین میانگین کل بالاتر باشد، بدین مفهوم است که ژنوتیپهای آن صفت برای بهگزینی مناسب هستند.
بر اساس این نتایج میتوان عنوان کرد که ژنوتیپهای گروه سوم ازنظر صفات مؤثر مشخصشده در تعیین ضرایب استانداردشده کانونیکی شامل تعداد کپسول در بوته، وزن خشک اندام هوایی در بوته، وزن دانه تکبوته، زیستتوده و شاخص برداشت برتر هستند. باتوجهبه نقش تعداد کپسول در بوته بهعنوان عامـل مهـم مؤثر در عملکرد کنجـد میتواند معیــار مناســبی بــرای انتخاب ژنوتیپهای برتر محسوب شوند.
جدول 10- مقایسه میانگین در صفات مورد بررسی در ژنوتیپهای کنجد در شرایط آب و هوایی مشهد.
Table 10– Cluster analysis in seasem genotypes in mashhad condition.
Group 1 |
Eshgabad, Boshrooyeh, Feizabad, Afghan 8, Vahdati 2, Tabas, Kahnooj, Afghan 6 |
|
||
Group 2 |
Gorgan, Vahdati 1, Afghan 2, Torbat-e-Jam, Sarayan, Kalat, Davarzan |
|
||
Group 3 |
Afghan 1, Afghan 4, Afghan 5, Afghan 7, Afghan 3 |
|
||
Group number |
Group 1 |
Group 2 |
Group 3 |
Total |
Number of genotypes |
8 |
7 |
5 |
20 |
Traits |
Group mean |
Group mean |
Group mean |
Total mean |
Plant height (cm) |
82.7a±6.09 |
84.8a±8.36 |
86.0a±6.72 |
84.3±6.68 |
Lowest capsule height (cm) |
40.7a±11.20 |
34.9a±5.64 |
37.5a±10.1 |
37.9±8.08 |
Branch no. |
3.52a±0.78 |
3.20a±0.970 |
3.88a±0.870 |
3.50±0.885 |
Capsule no. plant-1 |
42.1b±14.20 |
59.2b±20.1 |
86.4a±20.5 |
59.2±23.0 |
Biomass (g.plant-1) |
32.9c±1.94 |
41.1b±3.54 |
54.8a±10.2 |
41.2±9.65 |
Grain (g.plant-1) |
1.84b±0.82 |
4.18b±2.65 |
9.70a±5.22 |
4.63±4.44 |
1000- grain weight (g) |
2.90a±0.62 |
2.97a±0.320 |
3.32a±0.570 |
3.04±0.492 |
Biological yield (kg.ha-1) |
3615c±214 |
4435b±263 |
6029a±1117 |
4506±1051 |
Grain yield (kg.ha-1) |
203b±90.7 |
460b±291 |
1068a±574 |
510±489 |
HI (%) |
5.69b±2.97 |
11.1ab±7.02 |
17.3a±7.73 |
10.5±7.45 |
حروف متفاوت نشاندهنده تفاوت معنیدار میان گروههای مختلف است.
Means with the same letter are not significantly different from each other.
نتیجهگیری کلی
نتایج این پژوهش نشان داد بین ژنوتیپهای مورد بررسی ازنظر اکثر صفات زراعی تنوع فراوانی مشاهده شد. صفات تعداد کپسول در بوته، تعداد شاخه فرعی و وزن هزار دانه ازجمله صفات بسیار مؤثر در عملکرد کنجد بوده و ازآنجاییکه عملکرد بالا یکی از اهداف مهم اصلاحی در کنجد است؛ بنابراین ژنوتیپهایی با این ویژگیها را میتوان بهعنوان ژنوتیپهای برتر معرفی کرد. همچنین استفاده از روش تجزیه خوشهای در تمایز ژنوتیپها به زیرگروههای مشابه بر اساس صفات زراعی نتیجه مطلوبی داشت. بهطوریکه ژنوتیپهای کنجد گروه سوم از عملکرد و اجزای عملکرد بالاتری برخوردار بودند؛ بنابراین با درنظرگرفتن نتایج حاصل از روشهای آماری مختلف، ژنوتیپهای افغان1، افغان3 و افغان5 بهعنوان ژنوتیپهای برتر بهویژه ازنظر عملکرد دانه شناخته شدند که دارای سازگاری بالایی به شرایط آبوهوایی مشهد بودند.
REFERENCES
REFERENCES