Document Type : Research Paper
Authors
Department of Plant Production and Genetics, Faculty of Agriculture, Urmia University, Urmia, Iran.
Abstract
Keywords
Main Subjects
. مقدمه
آفتابگردان گیاه دانه روغنی مهم با مصارف متعدد است (Radanovíc et al., 2018; Rauf, 2019; Fernandez et al., 2019). روغن آن با توجه به غنای اسیدهای چرب غیر اشباع و توکوفرول بهعنوان سالمترین روغن در رژیم غذایی جایگاه ویژهای دارد
(Rauf et al., 2017, 2020). با توجه به تغییرات اقلیمی و گسترش مناطقی که با تنشهای کمآبی و شوری مواجه هستند، پیشبینی میشود سطح زیر کشت آفتابگردان که واکنش نسبتاً متحمل (نیمهمتحمل) به تنشهای فوق نشان میدهد، بهطور قابل توجهی کاهش یابد (Dufresne et al., 2006; Shi et al., 2023). بنابراین، تولید و توسعه ژنوتیپهای متحمل پُرمحصول و پایدار هدف اصلی برنامههای بهنژادی به شمار میرود (Alizadeh et al., 2021). موفقیت هر برنامه بهنژادی به توانایی بهنژادگر در انتخاب ژنوتیپهای برتر برای صفات متعدد بهطور همزمان بستگی دارد. بایپلات ژنوتیپ-صفت، یک برنامه کاربردی از تکنیک بایپلات GGE است. از این بایپلات برای بررسی ژنوتیپها بر اساس دادههای صفات استفاده میشود؛ بهطوریکه انتخاب بر اساس چندین صفت با بالاترین بازده و موفقیت همراه باشد (Koç et al., 2018; Shojaei et al., 2022). این نوع از بایپلات اولین بار توسط یان معرفی شد (Yan, 2001). بایپلات ژنوتیپ-صفت (Genotype-Trait; GT) یک راهنمای گرافیکی مفیدی برای مطالعه و گزینش برترین ژنوتیپها بر اساس ترکیب صفات مختلف است (Motamedi & Safari, 2019). با GT هم رابطه بین صفات و هم رابطه بین ژنوتیپها قابل بررسی است. چنین بایپلاتی میتواند با تجسم همبستگیهای ژنتیکی بین صفات، شناسایی ژنوتیپها و همچنین صفاتی که میتوان از آنها در انتخاب غیر مستقیم برای بهبود (بهنژادی) یک صفت خاص استفاده کرد را آسانتر کند.
استفاده از بایپلات ژنوتیپ × صفت جهت ملموسکردن مزایا و معایب ژنوتیپهای خاص نیز پیشنهاد میشود؛ این امر هم برای ارزیابی تنوع و هم برای انتخاب والدین مهم است (Yan & Tinker, 2006; Kendal, 2022). بایپلات GT جهت ارزیابی برخی از ژنوتیپهای آفتابگردان در شرایط مختلف محیطی بهکار گرفته شده است (Yasar et al., 2023). همچنین از GT برای انتخاب ژنوتیپهای باقلا با عملکرد دانه بالا از طریق وزن صد دانه و ارتفاع بوته استفاده شده است (Koç et al., 2018). در مطالعه پیشِ رو از بایپلات GT، جهت ارزیابی 100 ژنوتیپ آفتابگردان دانه روغنی تحت شرایط آبیاری معمول و محدود استفاده شده است تا با در نظر گرفتن روابط بین صفات، صفات مطلوب برای گزینش و همچنین ژنوتیپهای با پاسخ مطلوب از نظر صفات در هر دو شرایط معرفی شوند.
تعداد 100 ژنوتیپ آفتابگردان که از مراکز مختلف تحقیقاتی تهیه شدند (جدول 1) در قالب طرح لاتیس ساده 10×10 در دو شرایط آبیاری معمول و محدود در روستای قزلجه سلماس با موقعیت جغرافیایی 38 درجه و 11 دقیقه و 50 ثانیه شمالی و 44 درجه و 45 دقیقه و 55 ثانیه شرقی و 2100 متر ارتفاع از سطح دریا، در سالهای 1392 و 1393 ارزیابی شدند. روستای قزلجه با توجه به گزارشهای آب و هوایی، دارای میزان متوسط بارندگی سالیانه 388 میلیمتر و میانگین دمای سالیانه 5/10 درجه سانتیگراد میباشد (http://en.climate-data.org/location/1784/). هر واحد آزمایش از یک خط کاشت به طول پنج متر تشکیل شد و بذور به فواصل 50 سانتیمتر روی خطوط کاشت و با فواصل 60 سانتیمتر بین خطوط کشت شدند. فاصله دو تکرار از هم 5/1 متر و فاصله دو آزمایش از هم پنج متر انتخاب شد. در هر دو سال مورد مطالعه، تا مرحله هشتبرگی آبیاری مزرعه در هر دو آزمایش (آبیاری معمول و محدود) بعد از 90 میلیمتر تبخیر از تشتک تبخیر کلاس A انجام گرفت. از مرحله هشتبرگی آبیاری در آزمایش آبیاری معمول تا انتهای فصل رشد به منوال قبل ادامه یافت. در آزمایش آبیاری محدود، از مرحله هشتبرگی به بعد، آبیاری بعد از 180 میلیمتر تبخیر از تشتک تبخیر کلاس A انجام شد. در هر کرت زراعی در هر تکرار تعداد 13 صفت شامل ارتفاع بوته (PH, cm)، تعداد برگ (LN)، طول برگ (LL, cm)، عرض برگ (LW, cm)، طول دمبرگ (PL, cm)، قطر ساقه (SD, cm)، کلروفیل (Ch)، روز تا گلدهی (DF, day)، روز تا رسیدگی (DM, day)، قطر طبق (HD, cm)، محتوای آب نسبی (RWC, %)، محتوای روغن دانه (Oil, %)، و عملکرد (Yield, g) اندازهگیری شدند.
جدول 1. مشخصات ژنوتیپهای آفتابگردان مورد مطالعه.
Country |
Line Name |
Code |
Country |
Line Name |
Code |
Country |
Line Name |
Code |
France |
H156A/RHA274 |
69 |
France |
H603R |
35 |
France |
H100A/83HR4 |
1 |
U.S.A |
SDB1 |
70 |
Iran |
4 |
36 |
France |
H209A/LC1064 |
2 |
U.S.A |
HAR-4 |
71 |
France |
703-CHLORINA |
37 |
France |
H205A/H543R |
3 |
France |
AS5305 |
72 |
France |
NSF1-A4*R5 |
38 |
France |
AS5306 |
4 |
U.S.A |
RHA274 |
73 |
Iran |
28 |
39 |
U.S.A |
RHA858 |
5 |
France |
H158A/H543R |
74 |
Iran |
30 |
40 |
France |
H209A/83HR4 |
6 |
France |
H100A/RHA274 |
75 |
Hungary |
F1250/03 |
41 |
France |
AS3211 |
7 |
France |
H209A/H566R |
76 |
U.S.A |
SDR18 |
42 |
France |
254-ENSAT |
8 |
France |
ASO-1-POP-A |
77 |
France |
LP-SCYB |
43 |
France |
AS5304 |
9 |
France |
AS6305 |
78 |
Serbia |
803-1 |
44 |
France |
1009329.2(100K) |
10 |
France |
B-FIPOPB |
79 |
France |
1009370-1(100K) |
45 |
France |
270-ENSAT |
11 |
U.S.A |
D34 |
80 |
France |
CSWW2S |
46 |
France |
AS613 |
12 |
France |
CAY |
81 |
France |
1009370-3(100K) |
47 |
France |
A-FLPOPA |
13 |
Iran |
346 |
82 |
France |
H158A/H543R |
48 |
France |
OES |
14 |
France |
NS-F1-A5*R5 |
83 |
France |
H100A |
49 |
France |
H100A/LC1064 |
15 |
Iran |
36 |
84 |
France |
15031 |
50 |
U.S.A |
RHA266 |
16 |
Iran |
38 |
85 |
France |
H250A/83HR4 |
51 |
France |
PAC2 |
17 |
France |
SDB2 |
86 |
U.S.A |
RHA265 |
52 |
France |
H157/LC1064 |
18 |
|
H158A/LC1064 |
87 |
|
PM1-3 |
53 |
France |
5DES20QR |
19 |
France |
H156A/H543R |
88 |
France |
RT948 |
54 |
France |
15038 |
20 |
France |
H543R/H543R |
89 |
|
283-ENSAT |
55 |
France |
1009337(100K) |
21 |
France |
H543R |
90 |
|
QHP-1 |
56 |
France |
AS3232 |
22 |
France |
SF076 |
91 |
U.S.A |
SDR19 |
57 |
|
12AASB3 |
23 |
|
B-FIPOPB |
92 |
U.S.A |
HA337B |
58 |
|
8ASB2 |
24 |
France |
SF085 |
93 |
France |
H100B |
59 |
|
9CSA3 |
25 |
|
SF092 |
94 |
Hungary |
B454/03 |
60 |
|
H049+FSB |
26 |
|
A-CONTROLPLASTIPIC |
95 |
U.S.A |
HA304 |
61 |
|
SSD-580 |
27 |
Iran |
10-59 |
96 |
France |
RT931 |
62 |
France |
5AS-F1/A2*R5AS-29- F1/A2*R2 |
28 |
France |
H-100A-90RL8 |
97 |
U.S.A |
HA335B |
63 |
France |
7CR1=PRH6 |
29 |
|
SF109 |
98 |
France |
NS_B5 |
64 |
France |
ENSAT699 |
30 |
France |
SF105 |
99 |
U.S.A |
SDB3 |
65 |
France |
SSD-581 |
31 |
|
SF-023 |
100 |
France |
LC1064C |
66 |
France |
TMB-51 |
32 |
|
|
|
France |
NS-R5 |
67 |
Iran |
11*12 |
33 |
|
|
|
U.S.A |
DM-2 |
68 |
Iran |
110 |
34 |
Num: The numerator degrees of freedom; Den: The denominator degrees of freedom.
2-1. تجزیههای آماری
برای تجزیه واریانس مرکب دادهها در قالب مدل آماری طرح پایه از نرمافزار SAS 9.4 و اسکریپت زیر استفاده شد.
Data A;
Input Year Rep Block Entry X;
Cards;
...
;
Ods Rtf;
Proc Print;
Run;
%Include 'C:\Mult.Sas';
Ods Output Diffs=Diffs Lsmeans=Lsmeans;
Proc Mixed Data=A;
Class Year Rep Block Entry;
Model X=Entry/Ddfm=Kr;
Random Year Rep(Year) Block(Year*Rep) Entry*Year;
Lsmeans Entry/Pdiff Adjust=Tukey;
%Mult (Trt=Entry);
Run;
Ods Rtf Close;
برای تجزیه بایپلات GT از مدل تجزیه به مؤلفههای اصلی تصحیحشده بر پایه میانگین صفت و استاندارد شده توسط انحراف استاندارد استفاده شد:
در این مدل، مقدار مورد انتظار برای ژنوتیپ i در ترکیب با صفت j، میانگین کل تمام ترکیبهای ژنوتیپ × صفت؛ اثر اصلی صفت j؛ و بهترتیب اثرهای اولیه برای ژنوتیپ i و صفت j؛ و بهترتیب اثرهای ثانویه برای ژنوتیپ i و صفت j؛ dj انحراف استاندارد فنوتیپی در صفت j و عبارت از مقدار باقیمانده توجیهنشده بهوسیله اثرهای اولیه و ثانویه است. بایپلات مورد نظر به وسیله رسم در مقابل و در مقابل در یک نمودار پراکنش ایجاد شد. برای بهکارگیری معادله فوق از رایجترین روش، یعنی تجزیه به مقادیر منفرد (Singular value decomposition; SVD) استفاده شد:
که در آن و به ترتیب مقادیر ویژه اولین و دومین مؤلفه اصلی بزرگ (PC1 و PC2) هستند. و به ترتیب بردارهای ویژه ژنوتیپ i برای PC1 و PC2 و و به ترتیب بردارهای ویژه صفت j برای PC1 و PC2 هستند. این روش آماری با جزئیات بیشتر توسط Yan & Kang (2003) و Yan (2014) توصیف شده است. تجزیه و تحلیل GT-Biplot با برنامه metan در نرمافزار R انجام گرفت. برای این منظور ماتریس دو طرفه ژنوتیپ-صفت بر اساس میانگین تکرارها در دو سال برای هر یک از شرایط تنظیم و به نرمافزارمعرفی شد.
بر اساس نتایج تجزیه واریانس بین ژنوتیپهای آفتابگردان در هر دو شرایط آبیاری معمول و محدود در رابطه با اغلب صفات مورد مطالعه اختلاف آماری معنیدار وجود دارد که حاکی از تنوع ژنتیکی برای صفات فوق است (جدول 2، 3 و 4).
3-1. شناسایی بهترین ژنوتیپها بر اساس صفات مورد مطالعه
در تحلیل بایپلات ژنوتیپ-صفت در شرایط آبیاری معمول، دو مؤلفه اول در مجموع 97/49% از تنوع کل را توجیه میکنند؛ از این مقدار 68/36 % توسط PC1 و 29/13% توسط PC2 توجیه میشود. این درصد نشان از پیچیدگی روابط بین صفات دارد
(Yan & Rajcan, 2002; Sheikh et al., 2022). در شرایط آبیاری محدود، 58/68 % از تنوع کل؛ 4/47 % توسط PC1 و 18/21 % توسط PC2، توجیه میشوند. در نمودارهای چندضلعی، ژنوتیپهای واقع در یک بخش که در کنار یک یا چند صفت خاص قرار دارند، از عملکرد خوبی نسبت به آن صفات برخوردار هستند. ژنوتیپهایی که نزدیک به مبدأ نمودار قرار میگیرند تغییرات آنچنانی در صفات نشان نمیدهند (Thomason & Phillips, 2006; Koç et al., 2018)؛ ولی بهدلیل دارا بودن سطح متوسط برای تمام صفات از پتانسیل بالقوه برخوردارند (Yan & Fregeau-Reid, 2008). ژنوتیپهای قرار گرفته در رأس در شرایط نرمال، شامل ژنوتیپهای شماره 8 (254-ENSAT)، 11 (270-ENSAT)، 21 ((K100)1009337)، 72 (AS5305)،71 (HAR-4)، 39 (ایرانی 28)، 24 (8ASB2)، 22 (AS3232) و 60 (B454/03) بودند. این ژنوتیپها به دلیل فاصله طولانیتری که از مبدأ بایپلات دارند، در برخی یا همه صفات جزو بهترین یا ضعیفترین ژنوتیپها محسوب میشوند (شکل 1).
ژنوتیپهای قرار گرفته در رأس شامل ژنوتیپهای شماره 8 (254-ENSAT)، 11 (270-ENSAT)، 22 (AS3232) و 60 (B454/03) در اکثر صفات بجز محتوای نسبی آب، محتوای روغن و کلروفیل برترین ژنوتیپها بودند و ژنوتیپهای موفقی از نظر اکثر صفات بهویژه عملکرد بهشمار میروند. در این میان، ژنوتیپ شماره 8 (254-ENSAT)، بهدلیل فاصله طولانیتری که از مبدأ بایپلات دارد نسبت به بقیه ژنوتیپها از واکنش مطلوبتری در صفات فوق برخوردار میباشد. تحت شرایط نرمال در ناحیه مربوط به ژنوتیپهایِ رأسیِ 72 (AS5305)، 71 (HAR-4) و 39 (ایرانی 28) صفت شاخصی قرار نداشت و ژنوتیپهای واقع در این ناحیه بهویژه ژنوتیپهای رأسی فوق از نظر هیچ صفتی دارای برتری خاص نسبت به دیگر ژنوتیپها نبودند؛ بنابراین در جمیع صفات نامطلوب ارزیابی شدند. ژنوتیپهای قرار گرفته در رأس در شرایط محدودیت آبیاری شامل ژنوتیپهای شماره 8 (254-ENSAT)، 10 [1009329(100K)]، 16 (RHA266)، 70 (SDB1)، 35 (H603R)، 84 (ایرانی 36)، 6 (H209A/83HR4)، 17 (PAC2)، 21 ( (K100)1009337)، 43 (LP-SCYB)، 20 (15038) و 4 (AS5306) بودند. در این میان، ژنوتیپ شماره 8 (254-ENSAT)، از نظر صفاتِ کلروفیل، روز تا رسیدگی، روز تا گلدهی، طول دمبرگ، تعداد برگ و ارتفاع بوته مطلوب بود. ژنوتیپهای رأسی 35 (H603R)، 84 (ایرانی 36) و سایر ژنوتیپهای موجود در ناحیه مربوطه واکنش خوبی از نظر صفت عملکرد و محتوای نسبی آب داشتند (شکل 1). بنابراین ژنوتیپهای رأسی فوق، از نظر صفات عملکرد و محتوای نسبی آب مطلوبترین بودند؛ در حالیکه، در ناحیه مربوط به ژنوتیپهای رأسی 17 (PAC2)، 6 (H209A/83HR4)، 20 (15038)، 43 (LP-SCYB)، 4 (AS5306)، 21 ( (K100)1009337) و 10 ((K100)1009329) هیچ صفت شاخصی حضور نداشت و ژنوتیپهای واقع در این ناحیه از نظر مجموع صفات ضعیف ارزیابی شدند. در مجموع تحت شرایط نرمال و تنش خشکی (محدودیت آبیاری) ژنوتیپ 8 (254-ENSAT)، از مقادیر بالایی در اکثر صفات برخوردار بود. استفاده از بایپلات چندضلعیِ ژنوتیپ در صفت برای مقایسه، انتخاب ژنوتیپ و در تعیین رابطه بین صفات و ژنوتیپ، در گیاهان متعددی همچون سویا (Yan & Rajcan, 2002)، کلزا (Hosseini, 2016)، ذرت (Shojaei, 2020)، باقلا
(Sheikh et al., 2022)، توتون (Porkabiri et al., 2019a; 2019b) و سورگوم (Kaplan et al., 2017) گزارش شده است.
در نمودار بایپلات، روابط متقابل میان صفات از طریق بررسی زاویه بین بردارهای صفات امکانپذیر میشود. در واقع ضریب همبستگی بین جفت صفات با کسینوس زاویه بین بردارها برآورد میشود (Kendal, 2022). بر این اساس زاویه بین 0-90 نشاندهنده وجود رابطه مثبت بین دو صفت است. افزایش مقدار زاویه از 90 به180 حاکی از وجود رابطه منفی بین صفات است
(Kendal, 2022). همچنین طول بردار نیز تعیینکننده میزان تأثیر صفت است (Yan & Tinker, 2005). هر چقدر طول بردار اتصالدهنده صفت به مبدأ بیشتر باشد، گویای انحراف معیار بیشتر صفت و تنوع بالا بین ژنوتیپها است (Sheikh et al., 2022; Yan et al., 2007). تحلیل روابط از طریق بایپلات در شرایط نرمال (شکل 2)، نشان از وجود ارتباط مثبت قوی بین اکثر صفات مورد بررسی است. این صفات توسط زوایای منفرجه کوچک بین بردارهای آنها قابل تمیز هستند (r = cos 0 ° = +1). مماسبودن بردار طول دمبرگ با قطر طبق و بردار قطر ساقه با عملکرد نشانگر وجود همبستگی مثبت و قوی بین آنهاست که حاکی از ارتباط قوی ژنتیکی و فنوتیپی بین این صفات میباشد (Morovati et al., 2019). بین عملکرد با همه صفات بجز صفت کلروفیل همبستگی مثبت مشاهده شد. در این میان همبستگی مثبت و قوی بین طول برگ، عرض برگ، قطر ساقه، طول دمبرگ، قطر طبق، درصد روغن دانه با عملکرد مشاهده شد؛ بنابراین انتخاب در جهت افزایش این صفات با افزایش عملکرد همراه خواهد بود. در نتیجه از این صفات مورفولوژیک میتوان جهت انتخاب غیر مستقیم برای افزایش عملکرد تحت شرایط نرمال استفاده کرد.
Sources |
Conditions |
PH |
|
|
LN |
|
|
LL |
|
|||
of variation |
|
F |
Num |
Den |
|
F |
Num |
Den |
|
F |
Num |
Den |
Genotype |
Normal |
3.20*** |
99 |
93 |
|
1.92 |
99 |
95.1 |
|
2.23*** |
99 |
93.5 |
Stress |
2.30*** |
99 |
93.5 |
|
1.54* |
99 |
94.3 |
|
1.86** |
99 |
95.2 |
|
|
|
WL |
|
|
|
PL |
|
|
|
SD |
|
|
|
|
F |
Num |
Den |
|
F |
Num |
Den |
|
F |
Num |
Den |
Genotype |
Normal |
1.87** |
99 |
91.3 |
|
1.86** |
99 |
92.7 |
|
2.09*** |
99 |
93.1 |
Stress |
1.49* |
99 |
94.4 |
|
1.88** |
99 |
91.8 |
|
1.51* |
99 |
93.1 |
|
|
|
Ch |
|
|
DF |
|
|
DM |
|
|||
|
|
F |
Num |
Den |
|
F |
Num |
Den |
|
F |
Num |
Den |
Genotype |
Normal |
2.01*** |
99 |
297 |
|
2.68*** |
99 |
87.3 |
|
3.80*** |
99 |
85 |
Stress |
1.77*** |
99 |
269 |
|
2.74*** |
99 |
84.5 |
|
1.97*** |
99 |
91.4 |
|
|
|
HD |
|
|
RWC |
|
|
OIL |
|
|||
|
|
F |
Num |
Den |
|
F |
Num |
Den |
|
F |
Num |
Den |
Genotype |
Normal |
0.03 |
99 |
99 |
|
0.88 |
99 |
83.1 |
|
6.76*** |
99 |
98.8 |
Stress |
1.12 |
99 |
95.8 |
|
1.36 |
99 |
90.7 |
|
7.89*** |
99 |
98.6 |
|
|
|
Yield |
|
|
|
|
|
|
|
|
|
|
|
|
F |
Num |
Den |
|
|
|
|
|
|
|
|
Genotype |
Normal |
2.1*** |
99 |
88.9 |
|
|
|
|
|
|
||
Stress |
0.95 |
99 |
94.9 |
|
|
|
|
|
|
|
|
*، ** بهترتیب معنیدار در سطح احتمال پنج و یک درصد، ارتفاع بوته (PH, cm)، تعداد برگ (LN)، طول برگ (LL, cm)، عرض برگ (LW, cm)، طول دمبرگ (PL, cm)، قطر ساقه (SD, cm)، کلروفیل (Ch)، روز تا گلدهی (DF, day)، روز تا رسیدگی (DM, day)، قطر طبق (HD, cm)، محتوای نسبی آب (RWC, %)، محتوای روغن دانه (Oil, %)، و عملکرد (Yield, g).
جدول 3. برآورد واریانس اثرات تصادفی در طرح لاتیس ساده تحت شرایط آبیاری معمول طی دو سال زراعی
Experimental error |
Genotype × Year |
Block (Year × Replication) |
Replication / Year |
Year |
Character
|
151.89 |
41.854 |
63.825 |
7.101 |
77.90 |
Plant height |
8.423 |
2.764 |
0.943 |
0 |
6.381 |
Number of leaf |
4.776 |
1.699 |
0.281 |
0.0525 |
3.517 |
Leaf length |
6.836 |
2.395 |
0.214 |
0.039 |
4.759 |
Leaf width |
2.146 |
0.790 |
0.108 |
0.158 |
4.243 |
Petiole length |
0.541 |
0.197 |
0.051 |
0.025 |
0.823 |
Stem diameter |
21.378 |
0 |
0 |
0.314 |
3.999 |
Chlorophyll |
14.68 |
1.926 |
0.509 |
0.422 |
8.895 |
Day to flowering |
20.207 |
1.235 |
0.137 |
0.277 |
0.204 |
Day to maturity |
1.443 |
8.557 |
0.360 |
0.241 |
1.920 |
Head diameter |
95.274 |
4.416 |
1.577 |
0 |
49.54 |
RWC |
0.480 |
3.449 |
0.012 |
0.005 |
1.862 |
Oil performance |
RWC; Relative water content
قویترین همبستگی مثبت بین صفت عملکرد با قطر ساقه مشاهده شد و با توجه به طول بردار قطر ساقه این صفت بیشترین تأثیر را بر عملکرد بعد از روز تا گلدهی، عرض برگ، و طول برگ دارا میباشد و دارای تنوع بالا بین ژنوتیپها است. زاویه بین بردارهای دو صفت محتوای نسبی آب با صفت عرض برگ عمود برهم است (r = cos 90 ° = 0) که حاکی از عدم وجود همبستگی بین این دو صفت میباشد (شکل 2(. محتوای نسبی آب با صفات روز تا گلدهی، روز تا رسیدگی، ارتفاع بوته، تعداد برگ، طول برگ، طول دمبرگ، قطر طبق و کلروفیل همبستگی منفی نشان داد. همچنین با توجه به طول بلند بردارهای صفات عملکرد، قطر ساقه، طول برگ، عرض برگ، طول دمبرگ، قطر طبق، ارتفاع بوته، روز تا رسیدگی و روز تا گلدهی، به نظر ژنوتیپها طیف گستردهای از پراکندگی را برای این صفات داشتند؛ بهعبارتی برای این صفات تنوع بالا بین ژنوتیپها تحت شرایط نرمال مشاهده شد.
جدول 4. برآورد واریانس اثرات تصادفی در طرح لاتیس ساده تحت شرایط آبیاری محدود طی دو سال زراعی.
Experimental error |
Genotype × Year |
Block (Year × Replication) |
Replication / Year |
Year |
Character |
67.746 |
66.62 |
66.775 |
0 |
248 |
Plant height |
6.742 |
5.896 |
3.189 |
0.193 |
17.364 |
Number of leaf |
2.705 |
1.379 |
1.295 |
0.064 |
2.589 |
Leaf length |
2.961 |
1.628 |
1.210 |
0.042 |
1.825 |
Leaf width |
0.906 |
0.308 |
0.389 |
0.042 |
0 |
Petiole length |
0.33 |
0.131 |
0.056 |
0 |
0.485 |
Stem diameter |
23.108 |
0 |
1.033 |
0 |
0 |
Chlorophyll |
12.542 |
0.670 |
0.177 |
0.125 |
9.321 |
Day to flowering |
15.638 |
7.126 |
0.030 |
0 |
0 |
Day to maturity |
3.439 |
1.823 |
1.122 |
0.32 |
2.312 |
Head diameter |
76.561 |
7.646 |
14.143 |
26.725 |
3.145 |
RWC |
0.737 |
3.086 |
0.067 |
0.073 |
3.058 |
Oil performance |
107.18 |
59.731 |
16.191 |
0.8689 |
0 |
Yield |
RWC; Relative water content
بر پایه نتایج جدول همبستگی بیشترین همبستگی تحت شرایط نرمال با عملکرد بهترتیب مربوط به صفات قطر طبق (601/0)، عرض برگ (571/0) و قطر ساقه (533/0) است؛ درحالیکه در بایپلات ژنوتیپ× صفت بین عملکرد با صفت قطر ساقه همبستگی مثبت و قوی مشاهده شد که نشان میدهد بین نتایج بایپلات با جدول همبستگی ساده (جدول 5) مطابقت کامل وجود ندارد و این بهدلیل نحوه عملکرد متفاوت آنها است (Yan & Rajcan, 2002). به نظر انتخاب بر اساس نتایج بایپلات با موفقیت بیشتری همراه باشد.
تحت شرایط تنش خشکی (محدودیت آبیاری)، بین اکثر صفات همبستگی وجود داشت. بین عملکرد با اکثر صفات بجز محتوای کلروفیل، تعداد برگ، روز تا رسیدگی و روز تا گلدهی همبستگی مثبت وجود داشت. بین روز تا رسیدگی، روز تا گلدهی، طول برگ، طول دمبرگ، کلروفیل، ارتفاع بوته و تعداد برگ با محتوای نسبی آب همبستگی منفی مشاهده شد. در تنش خشکی (محدودیت آبیاری) بهدلیل همبستگی مثبت موجود بین صفت قطر طبق با اکثر صفاتی که همبستگی مثبت با عملکرد دارند، انتخاب غیر مستقیم بر اساس صفت قطر طبق مطلوبتر خواهد بود. تحت شرایط تنش خشکی (محدودیت آبیاری) صفات عملکرد، قطر ساقه، طول برگ، عرض برگ، قطر طبق، روز تا رسیدگی و روز تا گلدهی بیشترین تأثیرپذیری یا تنوع را در ژنوتیپها نشان دادند. در این بین بیشترین تأثیر مثبت مشاهدهشده بین صفات با عملکرد مربوط به صفات قطر طبق، قطر ساقه، عرض برگ، و طول برگ و کمترین تأثیر مثبت مربوط به صفت درصد روغن دانه بود. بیشترین تأثیر منفی مربوط به صفات روز تا گلدهی و روزتا رسیدگی و کمترین تأثیر منفی بر عملکرد مربوط به صفت محتوای کلروفیل بود. با توجه به روابط موجود در شرایط نرمال با صفت قطر ساقه و در شرایط تنش خشکی با صفت قطر طبق، فقط با یک صفت (حذف اندازهگیریهای اضافی) میتوان اقدام به انتخاب کرد (Mohammadi, 2019). همبستگی بین صفت قطر ساقه و قطر طبق با عملکرد پیشتر نیز گزارش شده بود
(Gholizadeh et al., 2021; Zabet, 2017; Zabet et al., 2016).
میانگین صفات اندازهگیریشده ژنوتیپها با ترسیم محور مختصات میانگین تستر (Average Tester Coordination; ATC) و تعیین ژنوتیپآیدهال که شامل ترکیبی از صفات خوب است، مورد ارزیابی قرار گرفتند. بر این اساس ژنوتیپی برتر است که از نظر اکثر صفات مطلوب بوده و کمترین فاصله را با دایره مرکزی روی مختصات افقی میانگین تستر (ATC) داشته باشد. با توجه به محور قائم ATC، ژنوتیپهایی که در سمت راست قرار میگیرند از مقادیر متوسط به بالا در صفات واقعشده در این بخش برخوردارند (شکل 3 و 4). همچنین صفات واقع در این بخش در برتر بودن ژنوتیپها نقش بیشتری دارند (Najafi et al., 2022). بر این اساس تحت هر دو شرایط آبیاری معمول و محدود تمام صفات بجز محتوای نسبی آب بر شاخصبودن ژنوتیپها نقش دارند. تحت شرایط نرمال بهترتیب ژنوتیپ شماره 8، 22، 57 ،60 و 11 ژنوتیپهای برتر از نظر تمام صفات هستند و نزدیکترین فاصله را با ژنوتیپ ایدهآل دارند و ژنوتیپهای شماره 59، 42 و 72 بهترتیب بدترین ژنوتیپها از نظر تمام صفات بر اساس ژنوتیپ ایدهآل میباشند یا بهعبارتی در تمام صفات در سطح پایینتر از متوسط قرار دارند (شکل 3 و 4). ژنوتیپ شماره 57 از پایداری بیشتر یا تنوع کمتر در همه صفات مورد بررسی از جمله عملکرد تحت شرایط نرمال برخوردار است (شکل3). تحت شرایط تنش خشکی (محدودیت آبیاری) بهترتیب ژنوتیپهای شماره 8، 70، 11، 13و 69 ژنوتیپهای برتر از نظر تمام صفات معرفی میشوند. در این میان ژنوتیپ شماره 70 از نظر تمام صفات از بالاترین میانگین و پایداری برخوردار است؛ بنابراین بهعنوان ژنوتیپ ایدهآل در شرایط تنش خشکی (محدودیت آبیاری) معرفی میشود که چندین صفت خوب را در ساختار ژنتیکی خود دارا میباشد (Mohaddesi et al., 2017; Badu-Apraku & Akinwale, 2011). ژنوتیپ ایدهآل (70) یا ژنوتیپهای برتر که در نزدیکترین فاصله با ژنوتیپ ایدهآل قرار دارند را میتوان بهعنوان ژنوتیپهای بالقوه مطلوب در مطالعات آتی مورد استفاده قرار داد (Mohaddesi et al., 2017). ژنوتیپهای 19، 59 و 17 بهترتیب بدترین ژنوتیپها از نظر تمام صفات بر اساس ژنوتیپ ایدهآل میباشند (شکل 3 و 4).
شکل 1. چندضلعی GT-Biplot برای تعیین ژنوتیپ برتر در شرایط نرمال (سمت راست) و تنش خشکی (محدودیت آبیاری)(سمت چپ).
جدول 5. ضرایب همبستگی فنوتیپی بین صفات مختلف مورفولوژیکی در ژنوتیپهای آفتابگردان تحت شرایط نرمال و تنش خشکی (محدودیت آبیاری).
Character |
PH |
LN |
LL |
LW |
PL |
SD |
Ch |
DF |
DM |
HD |
RWC |
Oil |
Yield |
PH |
1 |
0.408 |
0.550 |
0.532 |
0.447 |
0.408 |
0.063 |
0.311 |
0.371 |
0.345 |
-0.105 |
0.137 |
0.369 |
LN |
0.382 |
1 |
0.233 |
0.229 |
0.104 |
0.238 |
0.060 |
0.171 |
0.208 |
0.091 |
-0.010 |
0.060 |
0.121 |
LL |
0.409 |
0.285 |
1 |
0.873 |
0.588 |
0.703 |
-0.047 |
0.329 |
0.381 |
0.504 |
-0.135 |
0.163 |
0.496 |
LW |
0.357 |
0.175 |
0.838 |
1 |
0.577 |
0.740 |
0.021 |
0.254 |
0.290 |
0.485 |
-0.05 |
0.012 |
0.571 |
PL |
0.322 |
0.134 |
0.530 |
0.491 |
1 |
0.493 |
-0.038 |
0.277 |
0.281 |
0.445 |
-0.105 |
0.130 |
0.460 |
SD |
0.369 |
0.201 |
0.733 |
0.771 |
0.482 |
1 |
0.011 |
0.204 |
0.205 |
0.506 |
-0.034 |
-0.013 |
0.533 |
Ch |
0.196 |
0.135 |
0.047 |
0.038 |
0.141 |
0.148 |
1 |
0.242 |
0.165 |
0.096 |
-0.024 |
0.01 |
0.131 |
DF |
0.278 |
0.228 |
0.242 |
0.166 |
0.234 |
0.148 |
0.172 |
1 |
0.808 |
0.273 |
-0.307 |
0.060 |
0.172 |
DM |
0.291 |
0.220 |
0.256 |
0.221 |
0.297 |
0.178 |
0.157 |
0.836 |
1 |
0.264 |
-0.253 |
-0.025 |
0.243 |
HD |
0.305 |
-0.028 |
0.550 |
0.588 |
0.295 |
0.656 |
0.017 |
-0.007 |
0.052 |
1 |
-0.062 |
0.149 |
0.601 |
RWC |
-0.099 |
-0.117 |
-0.172 |
0.004 |
-0.098 |
0.017 |
0.003 |
-0.263 |
-0.233 |
-0.030 |
1 |
0.026 |
0.077 |
Oil |
0.073 |
0.131 |
0.061 |
0.039 |
-0.054 |
0.068 |
-0.078 |
-0.043 |
-0.003 |
0.054 |
0.065 |
1 |
0.095 |
Yield |
0.220 |
-0.200 |
0.358 |
0.499 |
0.057 |
0.549 |
0.003 |
-0.154 |
-0.107 |
0.790 |
0.112 |
0.071 |
1 |
همبستگیهای بزرگتر یا مساوی 193/0 در سطح پنج درصد و همبستگیهای بزرگتر یا مساوی 251/0 در سطح یک درصد معنیدار میباشند. در جدول همبستگی بالای قطر اصلی همبستگی تحت شرایط نرمال بین صفات و پایین قطر اصلی همبستگی بین صفات تحت شرایط تنش خشکی (محدودیت آبیاری) را نشان میدهد. ارتفاع بوته (PH, cm)، تعداد برگ (LN)، طول برگ (LL, cm)، عرض برگ (LW, cm)، طول دمبرگ (PL, cm)، قطر ساقه (SD, cm)، کلروفیل (Ch)، روز تا گلدهی (DF, day)، روز تا رسیدگی (DM, day)، قطر طبق (HD, cm)، محتوای نسبی آب (RWC, %)، محتوای روغن دانه(Oil, %)، و عملکرد (Yield, g).
شکل 2. بایپلات نقشه همبستگی بین صفات مورد مطالعه در شرایط نرمال (راست) و تنش خشکی (محدودیت آبیاری) (چپ).
شکل 3. بایپلات رتبهبندی ژنوتیپهای آفتابگردان دانه روغنی در مقایسه با ژنوتیپ ایدهآل در شرایط نرمال (راست) و تنش خشکی (محدودیت آبیاری) (چپ).
شکل 4. بایپلات مقایسه ژنوتیپهای آفتابگردان دانه روغنی با ژنوتیپ ایدهآل از نظر صفات و پایداری در شرایط نرمال (راست) و تنش خشکی (محدودیت آبیاری) (چپ).
استفاده از بایپلات ژنوتیپ در صفت (GT-biplot)، ضمن صرفهجویی در زمان، به واسطه تفسیر مطلوب حجم وسیعی از دادهها کمک شایان توجهی به برنامههای بهنژادی میکند. در این مطالعه تحت هر دو شرایط آبیاری معمول و محدودیت آبیاری مطلوب، بهترتیب ژنوتیپهای شماره 8 و 59 بهعنوان بهترین و بدترین ژنوتیپها معرفی شدند. ژنوتیپهای 59، 42 و 72 تحت شرایط آبیاری معمول و ژنوتیپهای 19، 59 و 17 تحت شرایط محدودیت آبیاری بدترین ژنوتیپها از نظر جمیع صفات مرتبط با عملکرد شناسایی شدند. ژنوتیپ 70 بهعنوان ژنوتیپ ایدهآل تحت تنش خشکی (محدودیت آبیاری) ارزیابی شد. شناسایی ژنوتیپهای کاملاً متمایز کمک شایانی به برنامههای بهنژادی تولید ارقام هیبرید در برنامههای آتی میکند. با توجه به روابط صفات تحت شرایط نرمال؛ صفت قطر ساقه و تحت محدودیت آبیاری؛ صفت قطر طبق نقش مهم در گزینش غیر مستقیم برای عملکرد خواهند داشت. بایپلات GT، امکان بررسی چندجانبه و تلخیص اطلاعات جهت رسیدن به یک نتیجه جامع و کامل را مقدور میسازد.
Alizadeh, B., Rezaizad, A., Hamedani, M.Y., Shiresmaeili, G., Nasserghadimi, F., Khademhamzeh, H.R., & Gholizadeh, A. (2021). Genotype-environment interactions and simultaneous selection for high seed yield and stability in winter rapseed (Brassica napus) multi-environment trials. Agricultural Research, 10, 1–10.
Badu-Apraku, B., & Akinwale, R.O. (2011). Cultivar evaluation and trait analysis of tropical early maturing maize under strigainfested and striga-free environments. Field Crops Research, 121, 186-194.
Dufresne, J.L., Salas, Y., Mélia, D., Denvil, S., Tyteca, S., Arzel, O., Bony, S., Braconnot, P., Brockmann, P., Cadule, P., Caubel, A., & Chauvin, F. (2006). Simulation du climat récent et future par les modèles du CNRM et de l’IPSL. La Météorologie, SHF-29 Hydraulics Days Congress: Climate variations and hydrology, Lyon, 27-28 March 2007 – JL Dufresne- JF Royer. 55, 45-59.
Fernandez, O., Urrutia, M., Berton, T., Bernillon, S., Deborde, C., Jacob, D., Maucourt, M., Maury, P., Durufé, H., Gibon, Y., Langlade, N.B., & Moing, A. (2019). Metabolomics characterization of sunflower leaf allows discriminating genotype groups or stress levels with a minimal set of metabolic markers. Metabolomics, 15, 56.
Gholizadeh, A., Ghaffari, M., Payghamzadeh, K., & Kia, S. (2021). Study on relationships between yield and other agronomic traits and identification of the superior hybrids in sunflower (Helianthus annuus L.). Journal of Crop Breeding, 13(37), 171-184. (In Persian).
Hosseini, S.Z. (2016). Evaluation of drought tolerance in canola (Brassica napus L.) genotypes, using biplot analysis. Journal of Crop Breeding, 8, 192-202. (In Persian).
Kendal, E. (2022). Using biplot analysis technique to selection in tritikale breeding studies. Yuzuncu Yıl University. Journal of Agricultural Sciences, 32(1), 186-198.
Koç, S., Orak, A., Tenikecier, H.S., & Nezihi Sağlam, N. (2018). Relationship between seed yield and yield characteristics in faba bean (Vicia faba L.) by GGE-biplot analysis. Journal of Life Sciences, 12, 105-110.
Kaplan, M., Arslan, M., Kale, H., Kara, K., & Kökten, K. (2017). GT biplot analysis for silage potential, nutritive value, gas and methane production of stay-green grain sorghum shoots. International Journal of Agriculture and Natural Resources, 44(3), 2019.
Motamedi, M, & Safari, P. (2019). Evaluation of water deficient stress tolerance in some wheat cultivars and their hybrids using canonical discriminant analysis and genotype by trait biplot. Journal of Crop Breeding, 11(29), 104-116. (In Persian).
Mohaddesi, A., Erfani, R., Sharifi, P., Aminpanah, H., & Abouzar Abbasian, A. (2017). Studying the relationships between yield and yield components and stability of some of rice genotypes using biplot method. Cereal Research, 6(4), 411-421. (In Persian).
Mohammadi, R. (2019). Genotype by yield × trait biplot for genotype evaluation and trait profiles in durum wheat. Cereal Research Communications, 47(3), 541-551.
Morovati, Z., Farshadsfar, E., & Romena, M.H. (2019). Genetic evaluation of physiological traits related to drought tolerance in some bread wheat genotypes under rain-fed conditions. Agronomy and Plant Breeding, 15(2), 35-50. (In Persian).
Najafi, S., Mohammadi, R., Shooshtari, L., Etminan, A., & Mehrabi, A.M. (2022). Evaluation of morpho-physiological diversity of durum wheat genotypes using genotype × trait biplot method. Journal of Crop Breeding, 14(44), 211-226. (In Persian).
Porkabiri, Z., Sabaghnia, N., Ranjbar, R., & Maleki, H. (2019a). Genetic variation of some tobacco (Nicotiana tabacum L.) genotypes by morphological traits. Scientia Agriculturae Bohemica, 50(1), 1-7.
Porkabiri, Z., Sabaghnia, N., Ranjbar, R., & Maleki, H.H. (2019b). Morphological traits and resistance to Egyptian broomrape weed (Orobanche aegyptiaca Pers.) in tobacco under greenhouse condition. Australian Journal of Crop Science, 13(2), 287.
Radanovíc, A., Miladinovíc, D., Cvejíc, S., Jockovíc, M., & Jocíc, S. (2018). Sunflower genetics from ancestors to modern hybrids—A Review. Genes, 9, 528.
Rauf, S., Jamil, N., Tariq, S.A., Khan, M., Kausar, M., & Kaya, Y. (2017). Progress in modification of sunflower oil to expand its industrial value. Journal of the Science of Food and Agriculture, 97, 1997–2006.
Rauf, S., Ortiz, R., Shehzad, M., Haider, W., & Ahmed, I. (2020). The exploitation of sunflower (Helianthus annuus L.) seed and other parts for human nutrition, medicine and the industry. Helia, 43(73), 167–184.
Rauf, S. (2019). Breeding strategies for sunflower (Helianthus annuus L.) genetic improvement. In: Al-Khayri, J.M., Jain, S.M., & Johnson, D.V. (Eds.), Advances in plant breeding strategies. Industrial and Food Crops, 6, 637–673.
Sheikh, F., Nazari, H., & Fanaii, H. (2022). Graphic analysis of trait relations and stability of faba bean genotypes using the biplot method. Crop Production Journal, 15(2), 117-136. (In Persian).
Shi, H., Wu, Y., Yi, L., Hu, H., Su, F., Wang, Y., Li, D., & Hou, J. (2023). Analysis of QTL mapping for germination and seedling response to drought stress in sunflower (Helianthus annuus L.). Peer Journal, 11, e15275 http://doi.org/10.7717/peerj.15275.
Shojaei, S.H., Ansarifard, I., Mostafavi, K., Bihamta, M.R., & Zabet, M. (2022). GT biplot analysis for yield and related traits in some sunflower (Helianthus annuus L.) genotypes. Journal of Agriculture and Food Research, 10, 100370.
Shojaei, S.H., Mostafavi, K.H., Khosroshahli, M., & Bihamta, M.R. (2020). Assessment of genotype-trait interaction in maize (Zea mays L.) hybrids usi ng GGT biplot analysis. Food Science and Nutrition, 8, 5340-5351.
Thomason, W.E., & Phillips, S.B. (2006). Methods to evaluate wheat cultivar testing environments and improve cultivar selection protocols. Field Crops Research, 99(2–3), 87–95.
Yan, W. (2001). GGE biplot - a windows application for graphical analysis of multi-environment trial data and other types of two-way data. Agronomy Journal, 93, 1111-1118.
Yan, W., & Rajcan, I. (2002). Biplot analysis of test sites and trait relations of soybean in Ontario. Crop Science, 42, 11-20.
Yan, W., & Kang, M.S. (2003). GGE biplot analysis: A graphical tool for breeders, geneticists, and agronomist, CRC press, Boca Raton, FL, p. 288. https://doi.org/10.1201/ 9781420040371.
Yan, W. (2014). Crop variety trials: Data management and analysis. Wiley-Blackwell.
Yan, W., & Tinker, N.A. (2006). Biplot analysis of multi-environment trial data: Principles and applications. Canadian Journal of Plant Science, 86(3), 623-645.
Yan, W., & Tinker, N.A. (2005). An integrated biplot analysis system for displaying, interpreting, and exploring genotype × environment interaction. Crop Science, 45, 1004-1016.
Yan, W., Kang, M.S., Ma, B., Woods, S., & Cornelius, P.L. (2007). GGE biplot vs. AMMI analysis of genotype-by-environment data. Crop Science, 47, 643-653.
Yan, W., & Frégeau-Reid, J. (2008). Breeding line selection based on multiple traits. Crop Science, 48, 417-423.
Yașar, M., Çil, A.N., & Çil, A. (2023). Investigation of genotype × environment interaction in some sunflower (Helianthus annuus L.) genotypes in different environmental conditions. MAS Journal of Applied Sciences, 8(1), 41-55.
Zabet, M. (2017). Identification of superior genotypes of rapeseed by GT biplot and GGE biplot methodology in normal and stressed conditions. Iranian Journal of Field Crop Science, 48, 207-220. (In Persian).
Zabet, M., Samadzadeh, A.R., & Shorvarzi, A. (2016). Selection of the most effective traits on yield of sunflower under normal irrigation and drought stress conditions in Birjand region. Environmental Stresses in Crop Sciences, 8(2), 217-231. (In Persian).