Investigation of genotype × environment interaction using GGE biplot method on grain yield in sesame

مقدمه

کنجد (L. Seasamum indicum) گیاهی یکساله،  دولپه‌ای‌ و پیوسته گلبرگ^[1]، متعلق به راسته کاسنی‌داران[2] و خانواده کنجدها^[3] می‌باشد که ارتفاع آن در حدود یک متر است. این خانواده دارای 16 جنس و حدود 60 گونه است. گل‌ها در این گیاه به‌طور منفرد و در کناره برگ‌های قسمت انتهایی ساقه قرار گرفته است و شامل قطعات پنج تایی به هم پیوسته می‌باشد و تعداد پرچم‌ها در آن چهار عدد است. میوه‌ کنجد به صورت کپسول و محتوی دانه‌های کوچک مسطح و بیضوی است که به دلیل دارا بودن روغن قابل استخراج، تنها بخش مورد استفاده گیاه است. سیستم ریشه‌ای کنجد،  کاملا رشد یافته است که سبب مقاومت به خشکی و سازگاری بالایی نسبت به تغییرات محیطی می‌شود  Wikipedia, 2017)). از نظر ژنتیکی، تعداد کرومزومهای موجود در هر سلول کنجد برابر با 2n=2x=26 است (Nyongesa et al.,  2014). اصل و منشاء کنجد به هند و قاره آفریقا برمی‌شود، اما نخستین دست نوشته‌های مربوط به این دانه به حدود 3000 سال قبل از میلاد مسیح در تمدن  آشوری تعلق دارد. بیشترین کاربرد کنجد، به‌صورت مصرف خوراکی و صنایع روغن‌کشی بوده است و در صنعت داروئی نیز کاربرد دارد (Shokohfar & Yaghobi Nejad, 2012; Wikipedia, 2017). بر اساس آمار سازمان خواروبار و کشاورزی ملل متحد، میانگین عملکرد کنجد در دنیا و ایران، به‌ترتیب 803 و 691 کیلوگرم در هکتار بوده است (FAOStat, 2018).

پایداری به تولید ثابت محصول (ثبات عملکرد) در مکان‌ها و سال‌های مختلف گفته می‌شود (Fernandez, 1991). هدف مراکز به‌نژادی، به‌دست آوردن ارقام با سازگاری وسیع، عملکرد بالا، پایدار و مقاوم به انواع تنش‌های محیطی، جهت استفاده در برنامه‌های اصلاحی و توزیع در میان کشاورزان است (Rajaram, 1983). ژنوتیپ‌هایی که دارای اثر متقابل بالایی در محیط‌های مختلف می‌باشند، از ثبات عملکرد دانه کمتری برخوردارند و پایداری پائینی دارند (Agaei, 1993). وجود اثر متقابل ژنوتیپ با محیط، سبب کاهش پایداری عملکرد ارقام در برخی از محیط‌ها می‌شود (Pham & Kang, 1988). این اثر سبب بروز تفاوت‌های قابل ملاحظه‌ای بین تظاهر ژنوتیپ‌ها در محیط‌های مختلف شده و باعث کاهش ارتباط بین بروز صفات فنوتیپی و مقادیر ژنوتیپی صفات می‌شود (DeLacy et al., 1990). اگر اثر متقابل ژنوتیپ و محیط باعث تغییر در رتبه ژنوتیپ‌ها نشود، قابل چشم‌پوشی است، اما تغییر رتبه ژنوتیپ‌ها در محیط‌های مختلف باید مورد بررسی و ارزیابی قرار گیرد (Raiger & Prabhakaran, 2001). اثر متقابل ژنوتیپ با محیط برای پژوهشگران علوم اصلاح نباتات دارای اهمیت ویژه‌ای است و یکی از مسائل پیچیده برنامه‌های به‌نژادی برای تهیه ژنوتیپ پر محصول و پایدار به شمار می‌رود (Cornelius & Crossa, 1999; Yan et al., 2007). وجود اثر متقابل ژنوتیپ و محیط، باعث کاهش بازده روش‌های اصلاحی و سودمندی وسعت کشت ژنوتیپ‌های اصلاح می‌شود و محققین را وادار می‌سازد تا برای مکان‌های مختلف، ژنوتیپ‌های متفاوتی را اصلاح و جهت کشت معرفی نمایند (Becker & Leon, 1988). همچنین وجود تفاوت‌های قابل ملاحظه بین ژنوتیپ‌ها در محیط‌های مختلف، سبب کاهش میزان اثر ژنوتیپ می‌شود ( DeLacy et al., 1990; Abdollahi Nezhad et al., 2005). اثر متقابل ژنوتیپ و محیط، یکی از عوامل کند کننده برنامه‌های اصلاح و آزادسازی ارقام در مناطق مختلف است (Kang, 1997). نتیجه‌گیری از آزمایشات به‌نژادی و به‌زراعی، بدون شناخت و بررسی دقیق اثر متقابل ژنوتیپ و محیط، دارای اعتبار چندانی نمی‌باشد. یکی از روش‌های کاهش اثر متقابل ژنوتیپ و محیط، روش انتخاب ژنوتیپ‌های پایدار است، به‌طوری‌که ژنوتیپ‌های منتخب، حداقل تاثیرپذیری را از محیط دارند (Sadegh Zadeh Ahari et al., 2005).

جهت برآورد اثر متقابل ژنوتیپ و محیط، از تجزیه واریانس مرکب استفاده می‌گردد (Miller et al., 1959). بدین منظور آزمایشات یکنواختی طی سال‌ها و مکان‌های مختلف اجرا و پس از تجزیه آماری، واریانس فنوتیپی به اجزای اصلی و اثرات متقابل آنها تجزیه می‌شود و به‌نژادگر از این تفکیک واریانس، به میزان اثر متقابل ژنوتیپ با محیط پی می‌برد. در صورت عدم وجود اثر متقابل ژنوتیپ و محیط، متوسط اختلاف بین ژنوتیپ‌ها در محیط‌های مختلف ثابت خواهد بود (Xie & Mosjidis, 1996). مدل این روش به صورت زیر می‌باشد:

رابطه (1)

در این رابطه، : عملکرد ژنوتیپ gام در محیط eام و با تکرار rام، α و β : به‌ترتیب اثرات اصلی ژنوتیپ و محیط،  : اثر متقابل ژنوتیپ و محیط و  : اثر خطا می‌باشند (Zobel et al., 1988).

در روش بای‌پلات، با بهره‌گیری از نمایش گرافیکی اثر متقابل ژنوتیپ و محیط، پایداری ژنوتیپ‌‌های برتر و همچنین روابط بین محیط‌ها و شناسایی محیط‌های مناسب در برنامه‌های اصلاحی آینده را بررسی می‌شود (Yan, 2001b). منشاء بسیاری از تغییرات در عملکرد، محیط می‌باشد، اما ارزیابی ژنوتیپ‌ها، ارتباطی با اثر اصلی محیط ندارد و نیاز است تا با حذف آن، به متغیرهای ژنوتیپ و اثر متقابل ژنوتیپ با محیط توجه شود (Yan & Kang, 2002).

روش بای‌پلات، نوعی تجزیه به مؤلفه‌های اصلی، برای اثر اصلی ژنوتیپ و اثر متقابل ژنوتیپ در محیط می‌باشد که که نمودارهای گرافیکی ترسیمی، تصویری از روابط بین ژنوتیپ و محیط ارائه و اطلاعات ذیل را به محقق ارائه می‌نمایند:

1- بررسی روابط متقابل بین محیط‌ها: کوسینوس زاویه بین دو بردار محیطی را ضریب همبستگی بین محیط‌ها تعریف می‌کنند (Yan, 2001b). نتایج آزمایشات در محیط‌های دارای همبستگی مثبت، قابل تعمیم به یکدیگر است و می‌توان آزمایشات را در یکی از آن محیط‌ها انجام داد تا هزینه آزمایشات کاهش و کارآیی اصلاح افزایش یابد. همچنین طول بردار محیط‌ها، تقریبی از انحراف معیار درون هر محیط است و ژنوتیپ‌ها در محیط‌های دارای بردارهای بلندتر، از قابلیت تمایز بیشتری برخوردار می‌باشند ( Yan & Kang, 2002; Baraki & Gebremariam, 2018). علاوه بر موارد فوق، با استفاده از این بای‌پلات می‌توان محیط‌های مختلف آزمایش را گروه‌بندی و ابرمحیط‌‌^[4]ها را مشخص نمود.

۲- انتخاب ژنوتیپ برتر جهت هر محیط: در این نمودار دورترین ژنوتیپ‌ها در بای‌پلات، به‌وسیله خطوطی به هم متصل و یک چند ضلعی تشکیل می‌دهند. بر هر ضلع، یک خط عمود از مبدأ رسم و چند ضلعی به چندین بخش مختلف تقسیم می‌شود. ژنوتیپی‌هایی که در هر بخش قرار دارند، دارای بهترین عملکرد در درون مکان‌های همان بخش هستند (Burgueño et al., 2001; Yan, 2001a; Ansarifard et al., 2020).

3- میانگین عملکرد ژنوتیپ‌ها و پایداری آن‌ها: دو محور میانگین عملکرد و محور اثر متقابل ژنوتیپ و محیط در این نمودار قابل مشاهده است. ژنوتیپ‌های پایدار، دارای میانگین عملکرد بالا و حداقل فاصله با محور اثر متقابل ژنوتیپ و محیط هستند ( Yan & Hunt, 2002; Yan & Tinker, 2006; Saremi-Rad et al., 2020).

4- رتبه‌بندی محیط‌ها بر اساس محیط ایده‌آل فرضی: هدف، شناسایی محیطی است که بهینه‌ترین الگوی پاسخ ژنوتیپی را دارا باشد (Saremi-Rad et al., 2020).

5- رتبه‌بندی ژنوتیپ‌ها بر اساس ژنوتیپ ایده‌آل: از نظر تئوری، ژنوتیپی ایده‌آل شناخته می‌شود که هم میانگین عملکرد بالا و هم پایداری بالا داشته باشد. چنین ژنوتیپ ایده‌آلی را می‌توان به عنوان مرجعی برای مقایسه با سایر ژنوتیپ‌ها استفاده نمود (Karimzadeh et al., 2013; Kindeya et al., 2020).

اهداف این تحقیق شامل بررسی وجود اثر متقابل ژنوتیپ و محیط، شناسایی ژنوتیپ‌های پایدار، توصیه ژنوتیپ‌های پایدار برای هر یک از مکان‌های آزمایش و همچنین کاهش هزینه‌های اصلاح از طریق گروه‌بندی ابر محیط‌ها بود.

مواد و روش‌ها

به‌منظور بررسی سازگاری و پایداری ارقام مورد بررسی کنجد در شرایط آب و هوایی مختلف، 15 ژنوتیپ مختلف که فهرست آن‌ها در جدول 1 آورده شده است، در آزمایشی در چهار منطقه آب و هوایی شامل اراک، بیرجند، شیراز و کرج، با مشخصات گفته شده در جدول 2، طی سال‌های 96 و 97 مورد ارزیابی قرار گرفتند.

جدول 1- نام ژنوتیپ‌های کنجد کشت شده در این تحقیق

Table 1. Sesame genotypes used in this study

جدول 2 – مناطق آزمایشی و مشخصات آن‌ها

Table 2. Characteristics of the study sites

در فصل بهار و پس از گاورو شدن مزرعه و تأمین دمای پایه لازم برای جوانه‌زنی کنجد، با انجام شخم عمیق، دیسک‌زنی، ایجاد ردیف‌های کشت و کرت و نیز ایجاد نهرها، نسبت به عملیات آماده سازی زمین اقدام شد. مطابق بررسی وضعیت شیمیایی خاک، کود­های حاوی نیتروژن، فسفر و پتاسیم به‌ترتیب از منابع اوره (100 کیلو­گرم در هکتار)، سوپر فسفات تریپل (70 کیلو­گرم در هکتار) و سولفات پتاسیم (60 کیلو­گرم در هکتار) استفاده شدند. تلاش شد تا مشخصات زراعی اجرای طرح در همه مناطق یکسان باشد. آزمایش به‌صورت طرح بلوک‌های کامل تصادفی با سه تکرار اجرا شد. طول هر کرت دو متر بود و درون هر کرت، چهار خط با فاصله  ۵۰ سانتیمتر قرار داشت. فواصل کاشت 10 سانتیمتر  و عمق کاشت در حدود 5/2 سانتیمتر در نظر گرفته شد. تاریخ‌های کاشت و برداشت با توجه به شرایط آب و هوایی منطقه متفاوت بودند که تاریخ تقریبی آن طی دو سال انجام آزمایش در جدول شماره 2 آورده شده است. کاشت به صورت خطی و با دست انجام شد و بلافاصله پس از کاشت، آبیاری صورت پذیرفت. دوره‌های آبیاری بر اساس رژیم حرارتی منطقه و پتانسیل تبخیر آب صورت گرفت. تنک کردن در مرحله دو برگی گیاه و مبارزه با علف‌های هرز به صورت دستی صورت انجام شد.

برای ارزیابی صفت عملکرد دانه، بوته‌های یک متر میانی هر کرت برداشت و و وزن آن‌ها اندازه‌گیری شد و داده‌ها بر اساس کیلوگرم در هکتار محاسبه شدند.

از نرم افزار SAS برای تجزیه و تحلیل داده‌ها جهت وجود اثر متقابل ژنوتیپ در محیط و از برنامه آماری GenStat 12^th جهت ترسیم بای‌پلات‌ها استفاده شد.

نتایج و بحث

با توجه به نتایج جدول تجزیه واریانس مرکب (جدول شماره 3)، اثرات مستقیم ژنوتیپ و مکان غیرمعنی‌دار ارزیابی شدند که به معنی عدم تفاوت معنی‌دار ژنوتیپ‌ها در صفت عملکرد دانه و همچنین اختلاف ناچیز ژنوتیپ‌های مورد آزمون در چهار مکان مورد آزمایش می‌باشد؛ همچنین اثر متقابل ژنوتیپ در مکان نیز معنی‌دار نشد. اما اثر ساده زمان (سال) و متقابل سال در مکان، ژنوتیپ در سال و ژنوتیپ در مکان در سال در سطح احتمال یک درصد معنی‌دار شد که این موضوع نشان دهنده تاثیر محیط بر صفت عملکرد دانه می‌باشد. در جدول 4 نیز نتایج مربوط به عملکرد دانه هر ژنوتیپ در چهار مکان مورد آزمایش آورده شده است.

جدول 3- تجزیه واریانس مرکب صفت عملکرد دانه در دو سال و چهار مکان

Table 3. Combined variance analysis of grain yield in 2 years and 4 locations.

ns و ^**: به ترتیب غیرمعنی‌دار و معنی‌دار در سطح احتمال یک درصد.

^ns and ^**: Non significant and significant at 1% of probability level, respectively.

جدول 4- میانگین عملکرد دانه هر ژنوتیپ در هر یک از مکان‌های مورد آزمایش

Table 4 . Genotypes mean grain yield in each location

مقادیر برحسب کیلوگرم بر هکتار می‌باشند                                        Values are base on kg/ha

بر اساس روش GGE biplot، دو مؤلفه اصلی اول و دوم به‌ترتیب 06/84 و 40/8 درصد از تغییرات را توجیه نمودند. این بدان معنی است که بای‌پلات این دو مؤلفه، 46/92 درصد از تغییرات ژنوتیپ و اثر متقابل ژنوتیپ و محیط را توجیه می‌نماید که بیانگر اعتبار نسبتاً بالای نمودار بای‌پلات حاصل از این مطالعه در تبیین تغییرات ژنوتیپ همراه اثر متقابل ژنوتیپ در محیط است. به‌طور کلی، عدم تبیین بالای تغییرات موجود مجموع مؤلفه‌های اصلی اول و دوم، نشانگر ماهیّت پیچیده اثر متقابل ژنوتیپ و محیط خواهد بود (Yan & Tinker, 2005)، اما این امر به معنای عدم اعتبار بای‌پلات نمی‌باشد (Yan et al., 2007). اگر نمودار بای‌پلات قادر به توجیه حداقل 60 درصد از واریانس داده‌ها باشد، می‌توان برای تعیین ابر محیط‌ها از آن استفاده نمود (Yang et al., 2009).

طی آزمایشی که بر روی شانزده ژنوتیپ کنجد در سه مکان و طی چهار سال در کشور اوگاندا انجام شد، اثرات مستقیم سال و ژنوتیپ و همچنین اثرات متقابل سال در مکان، تیمار در مکان و تیمار در سال در مکان، بسیار معنی‌دار برآورد شد، درصورتی‌که اثر مستقیم مکان و متقابل تیمار در سال غیر معنی‌دار بودند؛ همچنین مؤلفه‌های اصلی اول و دوم 21/47 و 25/14 درصد ارزیابی شدند (Okello-Anyanga et al., 2016).

بر اساس تحقیقی که بر روی لاین‌های پیشرفته گیاه علوفه‌ای خللر انجام شد، کلیه منابع تغییر اعم از اثرات مستقیم و متقابل، معنی‌دار ارزیابی شدند که نشان‌گر وجود اثرات متقابل ژنوتیپ (لاین‌در × محیط می‌باشند (Ahmadi et al., 2016).

نمودارهای بای‌پلات بر اساس چهار الگوی زیر ترسیم و ارائه شد:

در شکل 1، بین مکان‌های اراک، بیرجند و شیراز، همبستگی مثبت و بالایی وجود داشت، بنابراین ژنوتیپ‌ها در این مکان‌ها، پاسخ مشابهی داشتند و تفاوت چندانی بین رتبه عملکرد ژنوتیپ‌ها مشاهده نشد و می‌توان از یکی از این سه مکان برای کاهش هزینه‌ها برای تحقیقات آینده استفاده نمود. همچنین محیط‌های بیرجند، شیراز و اراک با محیط کرج، تقریباً‌ فاقد همبستگی بودند و به عبارتی، ژنوتیپ‌ها در این دو مجموعه، روند عملکردی مستقلی از یکدیگر داشتند. محیط‌های کرج و شیراز، طول بردارهای بلندتری داشتند و ژنوتیپ‌ها در این محیط‌ها، از قابلیت تمایز بالاتری برخوردار بودند. همچنین مناطق مورد بررسی به دو ابر محیط کرج و ابر محیط اراک،‌ بیرجند و شیراز تقسیم شدند.

در شکل 2، ژنوتیپ‌های TS-3 (G2)، داراب 14 (G8)، محلی سیرجان (G11)، PF3 (G12)، محلی فارس (G14) و محلی خنداب (G15) در رئوس چندضلعی قرار گرفتند. ژنوتیپ‌‌های برتر (از نظر عملکرد و حداقل تاثیرپذیری از محیط) در محیط کرج که در رأس چند ضلعی بخش کرج قرار دارند، شامل محلی فارس (G14) و محلی خنداب (G15) و پس از آن ژنوتیپ قابل توصیه برای محیط مذکور، داراب 1 (G7) بود.

در محیط‌های اراک، بیرجند و شیراز، داراب 14 (G8) برترین ژنوتیپ و پس آن ژنوتیپ‌های صفی آباد 1 (G9)، محلی اهواز (G10) و محلی اصفهان (G13) عملکرد بیشتر و پایداری بالاتری از سایر ژنوتیپ‌ها از خود نشان دادند. ژنوتیپ‌های TS-3 (G2)، Yellow White (G3)، جیرفت 13 (G4)، اولتان (G5)، داراب 2 (G6)، محلی سیرجان (G11) و PF3 (G12) در هیچ محیطی برتری نداشتند و جزو ژنوتیپ‌های ضعیف در اکثر محیط‌ها بودند. همچنین ژنوتیپ‌های دشتستان 5 (G1) و داراب 2 (G6)، تقریباً در مرکز بای‌پلات قرار گرفتند که نشان‌دهنده عکس‌العمل یکسان و ضعیف آن‌ها در بیشتر محیط‌های مورد مطالعه بود.

بر مبنای شکل 3، ژنوتیپ داراب 14 (G8)، بالاترین میزان عملکرد و ژنوتیپ TS-3 (G2)، کمترین میزان میانگین عملکرد دانه را دارابودند. همچنین ژنوتیپ محلی سیرجان (G11)، بیشترین ناپایداری را از خود نشان داد. ژنوتیپ‌های صفی آباد 1 (G9) و محلی خنداب (G15)، برترین (از نظر عملکرد و پایداری) ژنوتیپ‌ها شناخته شدند، زیرا علاوه بر این‌که در سمت راست نمودار و خط رسم شده از محور PC1 قرار داشتند، دارای کمترین فاصله با خط ترسیمی از محور PC2 بودند. همچنین محیط شیراز به علت قرار گرفتن در نزدیکترین دایره نسبت به مرکز دوایر متحدالمرکز، به

عنوان بهترین ایستگاه شناخته شد و پس از آن به‌ترتیب اراک، کرج و بیرجند قرار گرفتند. محیط بیرجند به علت بیشترین فاصله از مرکز دوایر متحدالمرکز، به عنوان ضعیف‌ترین مکان در این تحقیق بود.مطابق بای‌پلات ترسیمی شماره 4، ژنوتیپ داراب 14 (G8) در مرکز بای‌پلات قرار گرفت و به عنوان ایده‌آل‌ترین ژنوتیپ شناخته شد و پس از آن ژنوتیپ‌های محلی اصفهان (G13)، صفی آباد 1 (G9) و محلی خنداب (G15) به عنوان ژنوتیپ‌های ایده‌آل شناخته شدند؛ ژنوتیپ‌های TS-3 (G2) و Yellow White (G3) نامناسبترین ژنوتیپ‌ها بودند.

در تحقیقی که بر روی کنجد در چند مکان مختلف صورت گرفت، اثر مستقیم ژنوتیپ و متقابل ژنوتیپ و محیط بسیار معنی‌دار ارزیابی شدند. همچنین دو مؤلفه اصلی اول و دوم، به‌ترتیب 50/92 و 94/4 درصد (در مجموع 44/97 درصد) از تغییرات اثر مستقیم ژنوتیپ و اثر متقابل ژنوتیپ و محیط را توجیه نمود. به‌علاوه با کمک نمودار بای‌پلات چند ضلعی تعیین بهترین ژنوتیپ‌ها برای هر یک از مکان‌ها و بای‌پلات میانگین صفت و پایداری ژنوتیپ‌ها، ارقام برتر و ایده‌آل شناسایی شدند (Jahanfar et al., 2014).

در بررسی بای‌پلات‌های رسم شده اثر متقابل ژنوتیپ و محیط 16 ژنوتیپ کنجد در سه مکان طی دو سال، مکان‌های آزمون به دو ابر محیط تقسیم و با رسم چند ضلعی ژنوتیپ‌های مناسب برای هر محیط مشخص شد. به‌علاوه بای‌پلات میانگین صفت و پایداری ژنوتیپ‌ها و محیط‌ها نیز ترسیم شد. بای‌پلات‌های دیگر ترسیمی این آزمون شامل بای‌پلات روابط بین محیط‌ها، رتبه‌بندی محیط‌ها بر اساس محیط ایده‌آل فرضی و رتبه‌بندی ژنوتیپ‌ها بر اساس ژنوتیپ ایده‌آل فرضی بودند (Okello-Anyanga et al., 2016).

نتیجه‌گیری کلی

در کل، برترین ژنوتیپ دارای بیشترین عملکرد و پایداری بالا در کلیه محیط‌هایی مورد تحقیق، داراب 14 (G8) بود. ژنوتیپ‌های برتر قابل توصیه برای محیط کرج شامل محلی فارس (G14) و محلی خنداب (G15) و برای محیط‌های اراک، بیرجند و شیراز داراب 14 (G8) بودند. بر اساس بای‌پلات روابط بین محیط‌ها، میزان تاثیرپذیری هر یک ژنوتیپ‌ها در سه محیط اراک، بیرجند و شیراز، تشابه زیادی با آن ژنوتیپ در محیط دیگر داشت و بنابراین جهت تحقیقات آینده می‌توان تنها یکی از این سه محیط را برای مراحل اولیه اصلاحی انتخاب کرد و هزینه‌ها را کاهش داد.

REFERENCES

1. Abdollahi Nezhad, K., Alishah, O. & Sirani, S. (2005). The assessment of genotype environment interaction and yield stability on new cotton hybrids using parametric methods. Agricultural Sciences and Natural Resources, 12(4), 71-79 (In Persian).
2. Agaei, M. (1993). Study of genotype× environment interaction in barely cultivars on Tabriz. Journal of Agricultural Science, 1(2), 28–40.
3. Ahmadi, j., Vaezi, B. & Pour-Aboughadareh, A.R. (2016). Evaluation of Forage Yield Stability of Advanced lines of Grass pea (Lathyrus sativa) by Parametric and Non- parametric Methods. Journal of Crop Breeding, 8(17), 149-159 (In Persian).
4. Ansarifard, I., Mostafavi, K., Khosroshahli, M., Reza Bihamta, M. & Ramshini, H. (2020). A study on genotype–environment interaction based on GGE biplot graphical method in sunflower genotypes (Helianthus annuus L.). Food Science & Nutrition, 8, 3327–3334.
5. Baraki, F. & Gebremariam, G. (2018). Genotype environment interaction and stability of oil content of sesame (Sesamum indicum) genotypes in Northern Ethiopia. Journal of Cereals and Oilseeds, 9(3), 20–28.
6. Becker, H. C. & Leon, J. (1988). Stability analysis in plant breeding. Plant Breeding, 101(1), 1–23.
7. Burgueño, J., Crossa, J. & Vargas, M. (2001). SAS programs for graphing GE and GGE biplots. Mexico: CIMMYT Int.
8. Cornelius, P. L. & Crossa, J. (1999). Prediction assessment of shrinkage estimators of multiplicative models for multi‐environment cultivar trials. Crop Science, 39(4), 998–1009.
9. DeLacy, I. H., Eisemann, R. L. & Cooper, M. (1990). The importance of genotype-by-environment interaction in regional variety trials. Genotype-by-Environment Interaction and Plant Breeding. Louisiana State University, Baton Rouge, LA, 287–301.
10. Food and Agriculture Organization. (2018). FAOStat: Food and agriculture data. Retrieved June 15, 2020, from http://www.fao.org/faostat/en/#data/QC
11. Fernandez, G. C. J. (1991). Analysis of genotype× environment interaction by stability estimates. HortScience, 26(8), 947–950.
12. Jahanfar, V., Mostafavi, K. & Bihamta, M. R. (2014). Evaluation of seed yield stability of sesame genotypes using GGE biplot method. In 1st Iternational and 13th Iranian Genetics Congress, 24-26 May, Tehran, (pp. 1–5), (in Persian).
13. Kang, M. S. (1997). Using genotype-by-environment interaction for crop cultivar development. In Advances in Agronomy , 62, 199–252).
14. Karimzadeh, R., Mohammadi, M., Sabaghni, N., Mahmoodi, A. A., Roustami, B., Seyyedi, F. & Akbari, F. (2013). GGE biplot analysis of yield stability in multi-environment trials of lentil genotypes under rainfed condition. Notulae Scientia Biologicae, 5(2), 256–262.
15. Kindeya, Y. B., Mekbib, F. & Alle, E. A. (2020). AMMI and GGE Bi-plot analysis for seed yield and oil content of sesame (Sesamum indicum) genotypes in Tigray, Northern Ethiopia. Journal of Agriculture and Crops, 6(4), 58–67.
16. Miller, P. A., Williams, J. C. & Robinson, H. F. (1959). Variety✕ environment interactions in cotton variety tests and their implications on testing methods 1. Agronomy Journal, 51(3), 132–134.
17. Nyongesa, B. O., Were, B. A., Gudu, S., Dangasuk, O. G. & Onkware, A. O. (2014). Genetic relationship between sesame (Sesamum indicum) and related wild species based on chromosome counts and isozyme markers. African Journal of Agricultural Research, 9(12), 1052–1060.
18. Okello-Anyanga, W., Rubaihayo, P., Gibson, P. & Okori, P. (2016). Genotype by environment interaction in sesame (Sesamum indicum) cultivars in Uganda. African Journal of Plant Science, 10(10), 189–202.
19. Pham, H. N. & Kang, M. S. (1988). Interrelationships among and repeatability of several stability statistics estimated from international maize trials. Crop Science, 28(6), 925–928.
20. Raiger, H. L. & Prabhakaran, V. T. (2001). A study on the performance of a few non-parametric stability measures using pearl-millet data. Indian J. Genet, 61(1), 7–11.
21. Rajaram, S. (1983). Adaptation, stability and high yield potential of certain 1B/1R CIMMYT wheats. In proceedings of 6th international wheat genetics , November 28-December 3, Held at Kyoto, Japan, pp. 613–621.
22. Sadegh Zadeh Ahari, D., Hosseini, S. K., HosseinPour, T., Jafarbay, J., Khalilzade, G. & Alizade Disaj, K. (2005). Staudy on adaptability and stability of grain yield in durum wheat lines in warm and semi-warm dryland areas. Seed and Plant, 21(4), 561-576 (In Tehran).
23. Saremi-Rad, A., Mostafavi, K. & Mohammadi, A. (2020). Genotype-environment interaction study base GGE biplot method for kernel yield in Sunflower (Helianthus annuus) cultivars. Journal of Crop Breeding, 12(34), 43–53.
24. Shokohfar, A. & Yaghobi Nejad, S. (2012). The effect of drought stress on yield components of Sesame (Sesumum Indicum) cultivars. Agronomy & Plant Breeding of Iran, 8(4), 19-29 (In Tehran).
25. (2017). Sesame. Retrieved January 7, 2021 from https://www.en.wikipedia.org/wiki/Sesame.
26. Xie, C. & Mosjidis, J. A. (1996). Selection of stable cultivars using phenotypic variances. Crop Science, 36(3), 572–576.
27. Yan, W. (2001a). A study on the methodolgy of cultivar evaluation based on yield trial data with special reference to winter wheat in Ontario.D. Thesis. Faculty of Agriculture University of Guelph, Canada.
28. Yan, W. (2001b). GGEbiplot—A Windows application for graphical analysis of multienvironment trial data and other types of two‐way data. Agronomy Journal, 93(5), 1111–1118.
29. Yan, W. & Hunt, L. A. (2002). 19 Biplot analysis of multi-environment trial data. http://www.ggebiplot.com/Yan&Hunt2003.pdf, accessed July 9, 2020.
30. Yan, W. & Kang, M. S. (2002). GGE biplot analysis: A graphical tool for breeders, geneticists, and agronomists. CRC press.
31. 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(2), 643–653.
32. Yan, W. & Tinker, N. A. (2005). An integrated biplot analysis system for displaying, interpreting, and exploring genotype× environment interaction. Crop Science, 45(3), 1004–1016.
33. Yan, W. & Tinker, N. A. (2006). Biplot analysis of multi-environment trial data: Principles and applications. In Canadian Journal of Plant Science, 86, 623–645).
34. Yang, R., Crossa, J., Cornelius, P. L. & Burgueño, J. (2009). Biplot analysis of genotype× environment interaction: Proceed with caution. Crop Science, 49(5), 1564–1576.
35. Zobel, R. W., Wright, M. J. & Gauch Jr, H. G. (1988). Statistical analysis of a yield trial. Agronomy Journal, 80(3), 388–393.

^[1] - Sympetalae

^[2] - Tubiflorae

^[3]-^- Pedaliaceae

 ^[4] -Mega Environment