Stability analysis for forage yield in sorghum lines

Authors

1 (Corresponding Author; Tel: 09153510393), Scientific Staff of khorasan Razavi Agricultural and natural Resource Research Center

2 , Scientific Staff of khorasan Razavi Agricultural and natural Resource Research Center

Abstract

The objectives of this study were to analyze genotype by environment(GE) interactions on the forage yield of 15 sorghum lines by
the additive main effects and multiplicative interaction(AMMI) model and to evaluate genotype(G), environment(E) and GE interactions
using stability parameter i.e., AMMI stability value(ASV) and Wrick, s ecovalance (W2
i). The trials were conducted at two research
station(Neyshabur and Torogh Mashhad) for three successive years2007-2009( ). The result’s of AMMI analysis for forage yield indicated
that the Genotype main effects, environment, and GE interactions as well as two first interaction principal components (IPCA1-2)
were significant. AMMI biplot was able to distinguish stable genotypes and environments with high and low genotype discrimination
ability. The genotypes 5 and 15 with higher mean yield than total mean were to be most stable genotypes, while the genotype 1 and 3
with the highest contribution to GE interaction, were to be the most instable. The result of recommended genotypes based on AMMI
analysis showed that the genotypes 1 and 3 were highly adapted to Neyshabur and genotypes 5 and 15 were adapted to Mashhad environments
during this study.

Keywords


مقدمه

سورگوم (Sorghum bicolor (L.) Moench) یکی از گیاهان خانواده غلات (Poaceae)  بوده و عقیده بر این است که قدمت اهلی شدن این گیاه مشابه سایر غلات است. این گیاه از نظر اهمیت مقام پنجم را بعد از گندم،  برنج، ذرت و جو دارد. در سال های اخیر کشت و کار سورگوم به خاطر نیاز کمتر این گیاه به آب در مقایسه با ذرت (3/1 تا 2/1 ذرت) و سازگاری بهتر به شرایط نامتعارف محیطی توسعه زیادی یافته است(کوچکی، 1367). تیپ های علوفه ای معمولا شامل سورگوم علوفه ای، سودان گراس و تلاقی های سودان گراس – سورگوم می باشند. هر کدام از این تیپ ها  خصوصیات رشد ویژه ای دارند که بر نحوه استفاده از آنها تاثیر می گذارد. حتی در داخل هر تیپ ارقام بسیار متفاوتی می تواند وجود داشته باشد. کمبود علوفه در ایران و ظرفیت تولید بالای سورگوم و سازگاری آن با شرایط اقلیمی ایران ایجاب می کند تا ارقام سازگار با شرایط آب و هوایی ایران و در داخل کشور تولید گردد و از هیبرید های سورگوم علوفه ای با ظرفیت تولید بالا که بیشتر از پدیده هتروزیس ناشی می شود بهره بیشتری گرفته شود. این گیاه با شرایط آب و هوایی ایران بخصوص مناطق گرم و خشک مثل سیستان و بلوچستان، کرمان، اصفهان و مناطق جنوبی سازگار است (فومن اجیرلو، 1371). یکی از معیارهای سازگاری ارقام سورگوم دریک ناحیه داشتن عملکرد بالا است و مشاهده شده که عملکرد کمی و کیفی علوفه در شرایط مختلف با صفات مرفولوژیکی مثل قطر ساقه، زمان گلدهی، ارتفاع بوته و نسبت برگ به ساقه بستگی دارد(Dong et al. 1994).

وجود اثرات متقابل ژنوتیپ و محیط مهم ترین چالش فرا روی اصلاحگران گیاه است. تفسیر اثرات متقابل، شناسایی محیط های هدف و معرفی ژنوتیپ های مناسب با سازگاری خصوصی و عمومی برای محیط های مورد مطالعه و تعیین ژنوتیپ های پایدار در سال های مختلف از اهداف مهم در بررسی ژنوتیپ ها در سال ها و مکان های مختلف می باشد. پدیده اثر متقابل ژنوتیپ-محیط به اصلاحگران گیاه کمک می کند تا در ارزیابی ژنوتیپ ها در محیط های مختلف، مکان های غیر ضروری را حذف نموده و در نتیجه موجب کاهش بیشتر هزینه ها شوند(Basford and cooper, 1998; Kang and Magari, 1996;  Shafi et al., 1992).

پارامتر های زیادی برای تجزیه اثر متقابل ژنوتیپ-محیط ارائه شده است. استفاده از روش مبتنی بر رگرسیون

(Finlay and Winkinson, 1963; Eberhart and Russell, 1966 )، در زمره اولین روش های مورد استفاده است. سایر روش های آماری از جمله روش امی(AMMI: Gauch and Zobel, 1996)، روش شمی(SHAMMI) یا مدل ضرب پذیر متغیر(Cornelius et al, 1996 ; Cornelius,1993) و مدل SREG  یا GGE biplot (Yan, 2001) نیز زیاد مورد استفاده قرار گرفته اند، ضمن اینکه برای تعیین مدل پایداری، معیار های ناپایداری هان(Huhn, 1996) که مبتنی بر تغییر مرتبه ارقام می باشند و روش مجموع رتبه کنگ (Kang, 1993 and kang, 1998) که از جمله روش های ناپارامتری هستند می توانند معیارهای معتبری برای انتخاب ژنوتیپ های برتر باشند.

روش امی، یک روش چند متغیره آماری می باشد که اثرات جمع پذیر ژنوتیپ، محیط و اثرات ضرب پذیر ژنوتیپ × محیط را توجیه نموده و تفسیر خوبی از اثر متقابل ژنوتیپ × محیط را ارائه می نماید(Ebdon and Gauch, 2002) . اگر از جزء جمع پذیر مدل، بدون در نظر گرفتن اثر متقابل برای توجیه واریانس آزمایشات استفاده شود، مدل مذکور را امی صفر و اگر از جزء ضرب پذیر امی که اثر متقابل را نیز در نظر می گیرد استفاده شود بسته به این که از کدام مؤلفه  استفاده شود، مدل  امی 1 تا امی F نامیده می شود(Gauch and Zobel, 1996 ).

مؤلفه های امی معیار معتبری برای بررسی پایداری ژنوتیپ ها و ارتباط بین ژنوتیپ ها و محیط می باشند. اگر چندین مؤلفه از مؤلفه های اثر متقابل در مدل امی از لحاظ آماری معنی دار باشند، نشان دهنده وجود اثر متقابل پیچیده ژنوتیپ × محیط می باشد(Gauch and Zobel, 1996). با استفاده از مدل امی، یک پارامتر پایداری معتبر توسط پرکاس(Purchase, 1997) برای پایداری ژنوتیپ ها ارائه شد که آنرا ارزش پایداری امی( ASV1) می نامند و در آن از دو مؤلفه اول امی برای این روش استفاده می گردد. همچنین از اکووالانس ریک به منظور تعیین سهم محیط ها در اثر متقابل می توان استفاده نمود(Isik and Kleinschmit, 2005). آلبرت(Albert, 2004) در بررسی اثر متقابل ژنوتیپ × محیط و برای تعیین پایداری هیبرید های ذرت روش های مختلف پایداری را با هم مقایسه کرد ولی در نهایت مدل  امی را مناسب ترین روش برای تجزیه پایداری معرفی کرد.

سچومن(Schoeman, 2003)، برای بررسی اثر متقابل ژنوتیپ × محیط در آفتابگردان از روش های مختلف پایداری استفاده نمود و بیان داشت که مدل امی نه تنها پایدارترین ژنوتیپ ها را معرفی می کند، بلکه می تواند بیانگر سازگاری خصوصی ارقام نیز باشد.

بنابراین هدف از این پژوهش: 1- تجزیه اثرات متقابل ژنوتیپ × محیط و 2- بررسی واکنش ژنوتیپ ها در محیط های مختلف و انتخاب ژنوتیپ سورگوم علوفه ای پایدار به محیط های مورد مطالعه با استفاده از روش های چند متغیره و تک متغیره می باشد.

 

 

مواد و روش‌ها

تعداد 15 لاین امید بخش سورگوم علوفه ای(جدول 1)، در دو منطقه نیشابور(طول جغرافیایی ۵۸ درجه و۱۵ دقیقه شرقی و عرض جغرافیایی ۳۵ درجه و 40 دقیقه شمالی با ارتفاع 1120 متر) و طرق مشهد(با طول جغرافیایی 59 درجه و 38 دقیقه شرقی و عرض جغرافیایی 16 درجه و 36 دقیقه شمالی با ارتفاع 985 متر از سطح دریا) به مدت سه سال زراعی (88-1386)، مورد بررسی و مقایسه قرار گرفتند. منطقه نیشابور در طبقه بندی اقلیمی جزء مناطق معتدل و مشهد معتدل سرد می باشند. هر آزمایش در هرمحیط در قالب طرح بلوک‌های کامل تصادفی در سه تکرار اجراء گردید. کلیه ارقام، حاصل از برنامه های به نژادی سورگوم درکشور بوده و همگی به صورت لاین خالص (Pure Line) از آزمایشات مقدماتی انتخاب شده بودند. هر رقم درهر تکرار در 4 خط به طول 5/6 متر به صورت خطی کشت شد که بعد از سبز شدن و در مرحله 6-4 برگی عملیات تنک مزرعه انجام و فاصله بین بوته ها در روی ردیف به 6 سانتی متر و تعداد بوته در روی هر خط  5/6 متری به 110 بوته تقلیل یافت. عناصر معدنی پر مصرف بر اساس آزمون های خاک و توصیه کودی برای هر دو منطقه برابر با 400 کیلو گرم اوره ،150 کیلو سولفات پتاسیم و 100 کیلوگرم فسفات آمونیوم در هکتار تعیین و اعمال شدند. کلیه مقادیر فسفات و پتاس و یک چهارم اوره قبل از کاشت مصرف شد. کود سرک مرحله اول به مقدار 100 کیلو گرم در هکتار اوره به صورت نواری در ابتدای مرحله رشد سریع بوته کنار پشته ها مصرف و پس از هر چین برداری(دو چین) نیز تکرار گردید.  

در این بررسی به خاطر تنوع ارقام از نظر طول دوره رشد و نمو، افزایش تعداد چین برداری ها و همچنین کاهش خطر مسمومیت توسط ماده اسید پروسیک موجود در ارقام سورگوم و بر اساس دستورالعمل های فنی این گیاه، برداشت محصول زمانی انجام گرفت که متوسط ارتفاع بوته کل ارقام در محدوده 200-160 سانتی متر بود. در شرایط اقلیمی نیشابور و مشهد حداکثر دو چین قابل برداشت می باشد. همچنین برای حفظ سرعت رشد مجدد پنجه ها بعد از برداشت، بوته ها از ارتفاع 15-10 سانتی متری از سطح خاک برداشت شده و بلافاصله توزین شدند. دور آبیاری ثابت و 8 تا 10 روز یکبار بود. برداشت از دو خط وسط با حذف 25 سانتی متر از طرفین هر خط و در مساحت 2/7 متر مربع انجام شد و دو خط کناری به عنوان اثر حاشیه ای حذف گردید.

قبل از تجزیه امی، یکنواختی واریانس خطاهای آزمایش با استفاده از آزمون لون(Leven,s test) بررسی گردید و سپس تجزیه امی بر روی مجموع عملکرد علوفه تازه ژنوتیپ‌ها در محیط‌های مختلف با استفاده از نرم افزار GenStat  انجام شد. به منظور تجزیه پایداری عملکرد علوفه لاین‌ها و ارقام مورد بررسی از مدل امی و از مؤلفه‌های اثر متقابل اول و دوم امی(IPCA1, IPCA2) به عنوان پارامتر‌های پایداری برای ژنوتیپ‌ها و محیط‌ها استفاده گردید(Annicchiarico, 1997). همچنین از مدل بای پلات امی جهت بررسی واکنش ژنوتیپ‌ها در محیط‌ها استفاده شد. بای پلات ها به دلیل نمایش گرافیکی واکنش ژنوتیپ‌ها و محیط ها در پدیده اثر متقابل ابزارهای مفیدی جهت شناسایی ژنوتیپ‌های سازگار به محیط می باشند(Suadric et al., 2006). به منظور تحلیل بهتر اثرات متقابل از آماره های پایداری امی (ASV, Purchase et al., 2000) و اکووالانس ریک (: Wrick, 1962 W2i) استفاده گردید. آماره ASV به دلیل این که اثرات دو مؤلفه اول اثر متقابل را به طور هم زمان مورد بررسی قرار می‌دهد نقش مهمی در تفسیر نتایج تجزیه امی دارد. ژنوتیپ‌ها و محیط‌های با مقادیر کم پارامتر ASV پایدار و با مقادیر بیشتر ناپایدار خواهند بود(Purchase et al., 2000). پارامتر پایداری اکووالانس ریک(W2i) که کمتر بودن آن برای ژنوتیپ‌ها و محیط‌ها نشان دهنده پایداری بیشتر آنها و سهم کمتر در اثر متقابل می‌باشد جهت کمک در تفسیر نتایج مورد استفاده قرار گرفت.

 

 نتایج و بحث

جدول تجزیه واریانس مرکب علوفه تر ژنوتیپ ها در جدول 2 ارائه شده است. نتایج جدول 2 نشان داد که بین سال ها، مکان های آزمایش و همچنین اثرات متقابل ژنوتیپ و محیط تفاوت معنی داری وجود دارد(p≤ 0.01). لذا تجزیه پایدارای ژنوتیپ ها امکان پذیر می باشد.

نتایج تجزیه واریانس برای اثرات اصلی افزایشی و ضرب پذیر در جدول 3 آمده است. تجزیه واریانس اثرات افزایشی نشان داد که تفاوت معنی‌داری بین ژنوتیپ ها، محیط و اثر متقابل ژنوتیپ و محیط در سطح احتمال یک درصد وجود دارد. در این مدل (امی)، بزرگی اثرات اصلی افزایشی برای ژنوتیپ، محیط و اثر متقابل ژنوتیپ ومحیط به ترتیب 67/9 درصد، 56/83 درصد و 78/6 درصد مجموع مربعات کل بود. بزرگی اثرات محیط بیانگر متفاوت بودن محیط‌ها بوده که باعث ایجاد تنوع در عملکرد علوفه ی ژنوتیپ‌ها شده است. نتایج آزمون معنی‌داری مؤلفه‌های اثر متقابل نیز نشان داد که دو مؤلفه اول اثر متقابل مدل امی در سطح احتمال یک درصد معنی‌دار می باشند(جدول 3). نتایج تجزیه واریانس غیر افزایشی نشان داد که اولین مؤلفه اثر متقابل 01/70 درصد و دومین مؤلفه اثر متقابل 17/16 درصد از تغییرات ژنوتیپ و محیط را توجیه می‌نمایند. این دو مؤلفه از لحاظ آماری در توجیه اثرات متقابل معنی‌دار بودند. بنابراین مدل دوم امی(AMMI2 ) که شامل دو مؤلفه اول اثر متقابل و اثرات افزایشی ژنوتیپ و محیط می‌باشد می‌تواند مورد استفاده قرار گیرد. به عبارت دیگر استفاده از مدل دوم امی به خوبی در تفسیر نتایج مفید می‌باشد.

به منظور بررسی روابط ژنوتیپ‌ها و محیط‌ها از نمایش گرافیکی بای پلات استفاده شد. در بای پلات شکل 1 محور افقی نمایانگر اثرات اصلی جمع پذیر یا میانگین عملکرد علوفه بر حسب تن درهکتار و محور عمودی اثرات متقابل ضربی یا مقادیر اولین مؤلفه اصلی(IPCA1 ) یعنی ضرایب عاملی، برای ژنوتیپ‌ها و محیط‌ها به طور جداگانه می باشد. در بای پلات مذکور دو جفت از داده‌ها روی محور‌ها نمایش داده شده‌اند. اولین جفت مربوط به میانگین عملکرد هر رقم (محور افقی) و مقادیر اولین مؤلفه اصلی هر ژنوتیپ(محور عمودی) و دومین جفت مربوط به میانگین عملکرد هر محیط (محور افقی) و مقادیر اولین مؤلفه اصلی هر محیط(محور عمودی) می باشد.

به منظور تجزیه واکنش ژنوتیپی (استفاده هم زمان از تفکیک‌های دسته بندی و بردار یابیPurchase et al.,2000;

 Croosa et al., 1996; Van Eeuwijk, 1992)، و بررسی دقیق تر پایداری ژنوتیپ ها و محیط ها، تجزیه خوشه ای بر اساس مقادیر اولین مؤلفه اصلی ژنوتیپ ها و همچنین محیط ها(جدول 4) انجام شد. تجزیه خوشه ای مقادیر اولین مؤلفه اصلی برای ژنوتیپ ها سه گروه ژنوتیپی را مشخص نمود که گروه اول شامل ژنوتیپ های شماره 12، 10، 8، 7 و 4 با مقادیر بالا و مثبت IPCA1، گروه سوم شامل ژنوتیپ های 1 و 3 با مقادیر بالا و منفی IPCA1، و گروه دوم شامل سایر ژنوتیپ ها بودند که کمترین مقادیر IPCA1 را به خود اختصاص داده بودند.(شکل 2) همچنین تجزیه خوشه ای روی مقادیر اولین مؤلفه اصلی محیط ها سه گروه را مشخص نمود. گروه اول شامل محیط های 1 و 5  بود که مقادیر مثبت و نزدیک به صفر را دارا بودند. گروه دوم شامل محیط های 4 و 6 بودند که مقادیر بالا و مثبت  IPCA1  را داشتند. گروه سوم شامل محیط های 2 و 3 بودند که مقادیر بالا و منفی IPCA1   را به خود اختصاص داده بودند.(شکل 4)

بررسی بای پلات شکل 1 نشان می دهد که ژنوتیپ های شماره 1 و 3 و همچنین محیط های E3 و E4 دارای اثرات متقابل بزرگ بوده و بیشترین تاثیر را در ایجاد اثر متقابل دارند. ژنوتیپ هایی که در مرکز بای پلات قرار گرفته اند، اثر متقابل نزدیک به صفر را دارند و دارای پایداری عمومی بیشتری هستند، ژنوتیپ های شماره 15، 5، 11، 13 و 9 دارای اثر متقابل کم می باشند، ولی ژنوتیپ های شماره 5 و 15 به علت داشتن میانگین عملکرد بالاتر از میانگین کل می توانند به عنوان ژنوتیپ های با پایداری مطلوب مورد توجه قرار گیرند.

 

تجزیه پایداری ژنوتیپ ها

مقادیر پارامتر های ژنوتیپی و رتبه بندی ژنوتیپ ها بر اساس این پارامتر ها در جدول 5 آمده است. ضرایب دو مولفه اول اثر متقابل به عنوان ساده ترین پارامتر های پایداری جهت انتخاب ژنوتیپ ها قبلا مورد استفاده قرار گرفته است(Annicchiarico, 1997; Purchase et al., 2000)) گراس گروبر و همکاران (Grausgruber et al., 2000) در گندم نان و محمدی,و همکاران(2008) در گلرنگ از ضرایب مولفه های اصلی اثر متقابل(IPCA) جهت گزینش ژنوتیپ های پایدار استفاده نمودند. درپژوهش حاضر، کمترین مقدار IPCA1 مربوط به ژنوتیپ های 11، 5، 15، 9، 14 و 13 و کمترین مقدار IPCA2 مربوط به ژنوتیپ های 15، 10، 13، 12، 7 و 5 بود. بر اساس مقادیر IPCA1   وIPCA2 ژنوتیپ های 5، 15 با میزان عملکرد بالاتر از میانگین کل پایدارترین ژنوتیپ ها بودند.

بیشترین میزان IPCA1 متعلق به ژنوتیپ های 1و 3 با میانگین عملکرد بیشتر از میانگین کل(32/98 تن در هکتار) بود. همچنین بیشترین ضرایب IPCA2 مربوط به ژنوتیپ های 3، 1و2 با میانگین عملکرد بیشتر از میانگین کل(32/98 تن در هکتار) بود. آماره پایداری امی(ASV) نیز نشان داد که ژنوتیپ های 5، 15 و 11 با کمترین مقدار ASV پایدارترین ژنوتیپ ها بودند. در حالی که ژنوتیپ های 1 و 3 با بالاترین عملکرد های علوفه، ناپایدارترین ژنوتیپ ها در مجموع محیط ها بودند. بر اساس پارامتر پایداری اکووالانس ریک، ژنوتیپ های 5، 13 و 11 پایدارترین ژنوتیپ ها بوده و به ترتیب با 55/0، 16/1 و 33/1 درصد کمترین نقش را در اثر متقابل ژنوتیپ و محیط داشتند. در صورتی که ژنوتیپ های 1 و 3 به ترتیب با 34/26 و 11/17 درصد بیشترین نقش را در این اثر متقابل دارا بودند(جدول 5).

      

 تجزیه پایداری محیط ها

محمدی و همکاران (2008) به منظور تجزیه بهتر اثرات متقابل ژنوتیپ و محیط از پارامتر های پایداری برای محیط ها استفاده نمودند. مقادیر این پارامتر ها و رتبه بندی محیط ها بر اساس این پارامتر ها در جدول 6 آمده است. جهت گزینش محیط های مناسب با قدرت بالا در تفکیک ژنوتیپ ها، محیط ها بایستی دارای مقادیر IPCA1   بالا و IPCA2  پائین باشند. بر اساس این دو پارامتر، یان و همکاران(Yan et al., 2000) درگندم و یان و راجاکان(Yan and Rajcan, 2002 ) در سویا، ژنوتیپ های پایدار به محیط ها و نیز محیط های با قدرت تفکیک ژنوتیپی بالا از سایر محیط ها را شناسایی نمودند. در تحقیق حاضر بر اساس IPCA1، محیط E1 با پتانسیل عملکرد بالا و محیط E5 با پتانسیل عملکرد پائین، دارای بیشترین پایداری و کمترین نقش در اثر متقابل بودند. محیط های E3 و E4 با بیشترین IPCA1  بیشترین نقش را در ایجاد اثر متقابل ژنوتیپ و محیط داشتند. کمترین IPCA2  محیطی به E6 و بیشترین آن به محیط های E2 و E3 تعلق داشت. محیط های ایده ال بر اساس IPCA1  بالاتر وIPCA2 کمتر عبارت بودند از E4 و E6 (با مقادیر مثبت IPCA1 ). بر اساس آماره پایداری محیطی امی(ASVj )، محیط های E1 با پتانسیل عملکرد بالا و E5 با پتانسیل عملکرد پائین دارای بیشترین پایداری و کمترین نقش در پدیده اثر متقابل بودند در صورتی که محیط های E3 با بیشترین پتانسیل عملکرد و E4 با کمترین پتانسیل عملکرد، بیشترین نقش در ناپایداری را داشتند. بر اساس پارامتر اکووالانس ریک، محیطE3 با 88/27 درصد بیشترین نقش و محیط E1 با 19/5 درصد کمترین نقش را در ایجاد اثر متقابل داشتند(جدول 6).

 

نتیجه گیری نهایی

تفاوت در رتبه ژنوتیپ ها در محیط ها نشان دهنده اثر متقابل ژنوتیپ و محیط می باشد، به طوری که این اثر متقابل با معنی دار شدن آن از لحاظ آماری در سطح یک درصد قابل توجه می باشد(جدول 3). در جدول 7، چهار ژنوتیپ توصیه شده بر اساس مدل امی برای هر محیط نشان داده شده است. بر اساس نتایج حاصل، ژنوتیپ های 15 و 5 بیشترین سازگاری را به شرایط مشهد در هر سه سال زراعی داشتند(جدول 7). برای نیشابور ژنوتیپ های 1 و 3 در هر سه سال زراعی بهترین ژنوتیپ ها با سازگاری ویژه بودند. بر اساس نتایج این آزمایش در مجموع سال ها و دو منطقه، ژنوتیپ های 15 و 5 قدرت پایداری و سازگاری عمومی بالایی با محیط های مورد مطالعه را دارا بودند.

پاورقی ها

ASV= AMMI Stability Value

 

 

 

فومن، عزیز. (1371). بررسی و مقایسه عملکرد ارقام سورگوم دانه ای. مجله نهال و بذر، 8 (3-4). ص. 32-28 .
کوچکی، علیرضا. (1367). جنبه های از مقاومت به خشکی در سورگوم. مجله علوم و صنایع کشاورزی. 2(2). ص.77-81.
Albert, M.j. A. (2004). A comparison of statistical methods to describe genotype × environment interaction and yield stability in multi- location maize trials. MSc. Thesis. Department of Plant Sci. The University of the Free State, Bloemfontein.
Annicchiarico, p.(1997).Joint regression vs AMMI analysis of genotype-environment interactions for cereals in Italy. Euphytica. Vol, 94. pp: 53-62.
Basford, K.E., and Cooper, M.(1998). Genotype by environment interaction and some considerations of their implication for wheat breeding in Australia. Aust. J. Agric. Res . Vol, 49. pp: 154-175.
Cornelius, P.I. (1993). Statistical tests and retention of terms in the additive main effects and multiplicative interaction model for cultivar trials. Crop Sci. Vol, 33. pp: 1186-1193.
Cornelius, P.L., Crossa, J. and Seyedsadr, M.S. (1996). Statistical tests and estimates of multiplicative models for GE interaction. In: kang, M.S. and H.G. Jr. Gauch (eds). Genotype-By- environment interaction. (pp.199-234).CRC Press, Boca Raton, Florida.
Croosa, J., Gauch, G, H. and Zobell, R.W.(1990). Additive main effects and multiplicative interaction analysis of two international maize cultivar trials. Crop Sci. Vol, 30. pp: 493- 500.
Dong, O.P., LTN, C. S., Ander, O., Johnson, J.J., and Hamilton, R.I.(1994). Selecting stable and high-yielding sorghum cultivars for semi-arid tropics. International Sorghum and Millet Newsletter. 74: 59-762.
Ebdon, J.S. and Gauch, H. G.(2002). Additive main effect and multiplicative interaction analysis of national turf grass performance trials: I Interpretation of Genotype × environment interaction.  Crop Sci. Vol,42. pp: 489- 496.
Eberhart, S. A. and Russell, W.A.(1966).Stability parameters for comparing varieties. Crop Sci. Vol, 6. pp: 36-40.
Finlay, K. W. and Wilkinson, G.N.(1963).The analysis of adaptation in a plant breeding program. Aust. J. Agric. Res. Vol, 14. pp: 742-754.
Gauch, H.G. and Zobel, R.W.(1996).AMMI analysis of yield trials. In: Kang, M.S. and H.G. Jr. Gauch(eds), Genotype- by- environment interaction. (pp. 85-122). CRC Press, Boca Raton, Florida.
Grausgruber, H., Oberforster, M., Werteker, M., Ruckenbauer, P., and Vollmann,j.(2000).Stability of quality traits in Australian grown winter wheat.  Field rops Res. Vol, 66. pp: 257-267.
Hayward, M., Bosemard, D. and Romagosa, L.(1993).Plant breeding. Chapman and Hall, UK.
Huhn, M.(1996).Nonparametric analysis of genotype × environment interaction by ranks. In: Kang, M. S. and H. G. Jr. Gauch(eds),  Genotype- by- environment interaction. (pp. 235-271). CRC Press, Boca Raton. Florida.
Isik, K., and Kleinschmit. J.(2005).Similarities and effectiveness of test environments in selecting and deploying desirable genotypes. Theor Apple Genet. Vol, 110. pp: 311-322.
Kang, M. S., and Magari, R.(1996).New developments in selecting for phenotypic stability in crop breeding. In: M.S. Kang, and H. G. Zobel (eds), Genotype- by- Environment interaction, 1-14. CRC Press, Boca Raton.
Kang, M.S.(1993).Simultaneous selection for yield and stability in crop performance genotype × environment interaction 239 trials: consequences for growers.  Agronomy Journal, Vol, 85. pp: 754-757.
Kang, M.S.(1998).Using genotype × environment interaction for crop cultivar development. Advances in Agronomy. Vol, 62. pp: 199-252.
Mohammadi, R., Pourdad S.S. and Amri, A.(2008).Grain yield stability of spring safflower (Carthamus tinctorius  L.).  Aust. J. of Agric. Res. Vol, 59. pp: 546-553.
Purchase, J.(1997). Parametric analysis to describe Genotype × environment  interaction and yield stability in winter wheat. Ph.D.University of the Free State, South Africa.
Purchase, J.L., Hatting, H., and Van Deventer, C.S.(2000).Genotype × environment interaction of winter wheat in south Africa: II. Stability analysis of yield performance.  South Africa Journal of plant and soil. Vol, 17(3)pp: 101-107
Schoeman, L. J.(2003).Genotype × environment interaction in sunflower(Helianthus annuus) in south Africa. MSc.Thesis, Department of Agronomy, University of the Free State, Bloemfontein.
Shafi, B., Mahler, K.A., Price, W.J., and Auld, D.L.(1992).Genotype × environment interaction effects on  winter rapseed yield and oil content.  Rop Science,Vol, 32. pp: 922-927.
Suadric, A., Simic, D., and Vratric, M.(2006).Characterization of genotype by environment interactions in soybean breeding program of southeast Europe.  Plant Breed. Vol,125. pp: 125-191.
Van Eeuwijk, F. F. (1992). Multiplicative models for genotype. Environment interaction in plant breeding. Stat. Appl. Vol, 4. pp: 393-406.
Wricke, G. (1962) Uber eine method zur refassung der okologischen streubretite in feldversuchen,  
Yan, W.(2001).GGEbiplot- a Windows application for graphical analysis of multi-environment trial data and other types of two- way data.  Agron. J. Vol, 93(5)pp:  1111-1118.
Yan, W., and Rajcan, I.(2002). Biplot analysis of test sites and trait relations of soybean in Ontario. Crop Sci. Vol, 42. pp: 11-20.
Yan, W., Hunt, L.A. Sheng Q., and Szlavnics, Z.(2000).Cultivar evaluation and mega- environment investigation based on the GGE biplot. Crop Sci. Vol, 40. pp: 597-605.