Research Article

Effect of Genotype × Environment Interaction on Grain Yield Determinants in Bread Wheat

Malak Atiq Ullah Khan* and Fida Mohammad

Department of Plant Breeding and Genetics, The University of Agriculture, Peshawar, 25130, Khyber Pakhtunkhwa, Pakistan.

Received | October 01, 2017; Accepted | November 09, 2017; Published | January 16, 2018

*Correspondence | Malak Atiq Ullah Khan, Department of Plant Breeding and Genetics, The University of Agriculture, Peshawar, 25130, Khyber Pakhtunkhwa, Pakistan; Email: [email protected]

Citation | Khan, M.A.U. and F. Mohammad. 2018. Effect of genotype × environment interaction on grain yield determinants in bread wheat. Sarhad Journal of Agriculture, 34(1): 54-62.

DOI | http://dx.doi.org/10.17582/journal.sja/2018/34.1.54.62

Keywords | Genotype, Environment, GEI, Yield determinants, Correlation analysis

Introduction

Wheat is the primary food source of people in Pakistan. At national level, wheat was grown on 9.2 million hectares which produced 25.1 million tons with an average yield of 2.7 ton ha-1, while in Khyber Pakhtunkhwa, it was grown on 0.73 million hectares which produced 1.2 million tons with an average yield 1.7 ton ha-1 (Pak. Beurau Stat. 2014-15). The average national wheat yield of Pakistan is far below than agro-technologically advanced countries of the world. In Khyber Pakhtunkhwa province the situation is even worse where yield is lower than the national average yield. Factors for lower grain yield include lack of irrigation water, inadequate rainfall, heat stress, unavailability of high yielding cultivars and absence of quality seed. Demand for wheat being a staple food crop, would rise with the increase in population. Wheat production could be increased either by developing high yielding cultivars or growing more area under cultivation. The choice for wheat cultivation on larger area is partial; however, development of high yielding wheat cultivars with wider adaptability would play significant role.

Cultivars performance largely depends on their genetic makeup (G), environment (E) and their interaction (GEI). Fluctuating response of genotypes across test environments is an usual phenomenon, known as GEI (Akcura et al., 2009; Karimizadeh et al., 2012b; Mohammadi et al., 2012). Yield potential of a cultivar is the result of its performance over locations and years. Therefore, stability analysis of genotypes is required in the presence of GEI to ascertain high yielding and relatively stable genotypes. For several test environments, the GEI governs the credentials of the most stable genotypes that are suitable for specific environment (Annicchiarico, 2002). Thus, the genotypes possessing genetic homeostasis are essential to increase average yield.

Wheat production can boost up through cultivars having broader genetic base and better performance under various agro-climatic conditions. In wheat, genetic improvement is slow process in nature however, the selective process of man can speed it up through appropriate management of environmental factors. Improvement gets complicated when a trait is environment-driven and selection gets more complex (Mohammad, 2011).Keeping in view the importance of GEI in reference to its application for identifying stable genotypes, the present experiment was conducted using 81 recombinant inbred lines (RIL’s) including two check cultivars viz. Janbaz and Atta-Habib. The experiment was planted across nine environments to; i) assess GEI effects on grain yield and its components, ii) identify high yielding and stable genotypes for grain yield, and iii) compute genetic correlation coefficients among yield and yield components in wheat.

Materials and Methods

To interpret genotype by environment interaction (GEI) for grain yield and associated traits, eighty one (81) wheat genotypes were tested over 9 environments of Khyber Pakhtunkhwa, Pakistan during 2013/14, 2014/15 and 2015/2016. Detailed description of the experimental sites, history of breeding material, experimental designs and statistical analysis are given under separate heading as under;

Description of experimental sites

Eighty-one genotypes including 79 F5:8 RILs and two check cultivars “Janbaz” and “Atta-Habib” were evaluated in nine environments during 2013/14, 2014/15 and 2015/16. During 2013/14, the experimental materials were planted at single location i.e. The University of Agriculture Peshawar (E-01) for evaluation and seed multiplication, whereas, the experiments were planted during 2014/15 and 2015/16 at The University of Agriculture Peshawar (E-02 and E-03, respectively), Cereal Crops Research Institute, Pirsabak Nowshehra (E-04 and E-05, respectively), Agricultural Research Station, Charsadda (E-06 and E-07, respectively) and Agricultural Research Station, Swabi (E-08 and E-09, respectively). Hereafter, these will be referred as E-01, E-02, E-03, E-04, E-05, E-06, E-07, E-08, and E-09. Agro-metrological features of test sites/environments including temperature, rainfall etc. are given in Table 1.

Table 1: Description of nine environments used for evaluation of 81 wheat during 2014-2016 cropping season.

Seventy-nine F5:8 RILs were originally developed in the department of Plant Breeding and Genetics during 2002/03 (Table 2). The segregating populations were advanced in bulk till F4 generation. In F5 generation, single heads were selected from the bulk populations based on agronomic fitness and disease resistance. The F5:6 heads were raised as head-to-row for seed multiplication and disease screening against stripe rust in particular. However, lines were found

Table 2: List of wheat RILs with pedigree.

Eighty-one genotypes including 79 F5:8 RILs and two check cultivars Janbaz and Atta-Habib were field-tested in nine environments during 2013/14, 2014/15 and 2015/16. List of genotypes with parentage is given in Table 2. Experimental materials were planted in 9×9 alpha lattice design with two replicates at each environment. Each plot had 6 rows of 5-meter length and a row-to-row space of 30 cm. Standard dose of nitrogen (120 kg ha-1) and phosphorous (80 kg ha-1) was applied. Uniform cultural practices required for wheat crop were followed throughout the growing season.

Traits measurement

Data were recorded on yield and yield components i.e. tillers-2,grains spike-1and 1000-grain weight following the procedure described by Sayre et al. (1997).

Statistical analysis

Analysis of variance: Data recorded on various yield and yield associated traits were statistically analyzed across nine environments appropriate for alpha lattice design using SAS computer program (SAS, 2009). Upon significant differences, means were separated using LSD test at 5% level of probability.

Correlation analysis: Genetic correlations for various morphological and yield traits were computed following the procedure of Singh and Chaudhry (1997) using the following formula:

Results and Discussion

Combined analysis of variance

Combined analysis of variance depicted highly significant differences among genotypes and genotype by environment interactions (GEI) for all traits. (Table 3). Significant GEI exhibited that the performance of genotypes was not stable across different environments and hence justified to proceed for further analysis of genotypes in individual environments. Mean performance of genotypes for various yield and associated traits are discussed below:

Mean Performance

Tillers-2: Productive tillers are considered as one of the important traits in wheat production. It has been estimated that approximately 30 to 50% of the grain yield in wheat comes from the main stem and 50 to 70% comes from the tillers under normal conditions (Elhani et al., 2007). Combined analysis of variance depicted significant differences (P≤0.01) among genotypes, environments and GEI for tillers m-2. Percent variation due to GEI (45.1%) was four times more than the combined effect of environments (6.1%) and genotypes (6.6%), exhibiting that mean performance and ranking of genotypes were mainly due to their interaction with environments (Table 3). Averaged over nine environments, tillers m-2 ranged from 128 to 180 tillers with the mean value of 153 tillers (Table 4). Maximum number of tillers was produced by G-79 (181 tillers), followed by genotypes G-08 (174 tillers) and G-41 (171 tillers), whereas minimum number of tillers were noted for G-67 (129 tillers). Among environments, tillers m-2 ranged from 92 to 229 in E-01; 91 to 229 in E-02; 77 to 224 in E-03; 100 to 198 in E-04; 76 to 220 in E-05; 99 to 275 in E-06; 94 to 229 in E-07; 97 to 178 in E-08 and 92 to 210 in E-09. Genotype G-26 produced maximum tillers in E-01 (229 tillers), G-79 in E-02 (229 tillers), G-39 in E-03 (224 tillers), G-6 (198 tillers) in E-04, G-79 in E-05 (220 tillers), G-31 in E-06 (275 tillers),

G-28 in E-07 (229 tillers) G-6 in E-08 and G-19 in E-09 (210 tillers). None of the genotypes completely dominated rest of the genotypes at all environments, showing site specific performance for tillers m-2. Environment E-06 (163 tillers) and E-09 (145 tillers) were declared as highly productive and less productive environments, respectively for of tillers m-2 (Figure 1).

Table 4: Descriptive statistics of 81 wheat genotypes for various traits across nine environments.

Grains spike-1: Pooled analysis of variance for grains spike-1 unveiled significant differences (P≤0.01) among genotypes, environment and GEI. The GEI explained 48.1% of the total variation, while genotypes and environments captured 18.7% and 0.4%, respectively. Significant differences among genotypes for grain spike-1 provide room for selecting desirable genotype; however, selection should be environment specific due to significant GEI effect (Table 3). Mean over nine environments, grains spike-1 ranged from 51 grains to 71 grains with a mean value of 60 grains (Table 4). Maximum number of grains spike-1 was observed for G-79 (71 grains), followed by genotype G-58 and G-56, while minimum number of grains spike-1 was noticed for G-40 (51 grains). G-79 produced maximum grains spike-1 in all environments except E-06, E-07 and E-09, indicating its superiority and stability across environments. Within environments, grains spike-1 ranged from 42 grains to 83 grains in E-01; 44 grains to 80 grains in E-02; 46 grains to 85 grains in E-03; 42 grains to 83 grains in E-04; 43 grains to 85 grains in E-05; 47 grains to 76 grains in E-06; 45 grains to 78 grains in E-07; 43 grains to 76 grains in E-08 and 47 grains to 73 grains in E-09 (Figure 2).

1000-grain weight: Combined analysis of variance revealed significant differences (P≤0.01) among genotypes, environments and their interactions. The GEI explained maximum variation of 41.1% of total sum of squares, which is four times larger than both genotypes (8.8%) and environment (2.2%) (Table 3). Mean over nine environments, 1000-grain weight ranged from 33.4 to 44.3 g with a mean value of 38.3 g (Table 4).  

Maximum 1000-grain weight was noted for G-72 (44.3 g), followed by G-36 and G-37, whereas minimum value was recorded for genotype G-55 (33.4 g). Among environments, 1000-grain weight ranged from 29.1 to 51.9 g in E-01; 27.1 to 50.2 g in E-02; 25.1 to 65.9 g in E-03; 27.6 to 55.7 g in E-04; 28.0 to 56.0 g in E-05; 27.3 to 51.4 g in E-06; 41.1 to 49.8 g in E-07; 27.6 to 51.1 g in E-08 and 32.9 to 50.5 g in E-09. Maximum 1000-grain weight was recorded for G-36 in E-01, G-08 in E-02 G-40 in E-03, G-14 in E-04, G-79 in E-05, G-46 in E-06, G-36 in E-07, G-37 in E-08 and G-37 in E-09. For 1000-grain weight, E-03 (39.9 g) and E-05 (37.2 g) were declared as highly productive and least productive environments, respectively (Figure 3).

Grain yield : Analysis of combined data exhibited significant differences (P≤0.01) among genotypes, environments and GEI for grain yield. Although, environments and genotypes were significant but they capturedleast sum of squares of 2.2% and 8.8%, respectively. Correspondingly, GEI captured 51.7% of the total variation, exhibiting its role in mean performance and ranking of genotypes across environments (Table 3). Averaged over nine environments, grain yield ranged from 2804 to 4862 kg ha-1 with mean value of 3438 kg ha-1 (Table 4). Overall, 33 out of 81 genotypes had higher yielding than mean yield, while 21 out of 79 had higher yield than both check cultivars. Genotype G-79 produced maximum grain yield (4862 kg ha-1), followed by G-08 (4195 kg ha-1), G-56 (4030 kg ha-1) and G-37 (4003 kg ha-1), whereas minimum value for grain yield was noticed for G-40 (2804 kg ha-1). Grain yield ranged from 2024 to 5272 kg ha-1 in E-01; 1980 to 4975 kg ha-1 in E-02; 1938 to 5102 kg ha-1 in E-03;1995 to 5106 kg ha-1 in E-04; 1931 to 5035 kg ha-1in E-05; 1897 to 5100 kg ha-1 in E-06; 1974 to 5024 kg ha-1 in E-07; 1941 to 4906 kg ha-1 in E-08 and 2110 to 4886 kg ha-1 in E-09. Genotype G-8 (5272 kg ha-1) produced maximum grain yield in E-01, G-9 (4975 kg ha-1) in E-02, G-49 (5071 kg ha-1) in E-03, G-79 (5106 kg ha-1) in E-04, G-79 (5035 kg ha-1) in E-05, G-79 (4920 kg ha-1) in E-06, G-79 (5024 kg ha-1) in E-07, G-79 (4906 kg ha-1) in E-08 and G-79 (4886 kg ha-1) in E-09. E-01 (5272 kg ha-1) and E-09 (4886 kg ha-1) appeared as high and least productive environments, respectively (Figure 4).

Correlation analysis: Correlation coefficients among yield and yield associated traits were computed based on mean data of 81 genotypes over nine environments. Correlation analysis depicted that grain yield had strong positive association with tillers m-2 (rg = 0.72**), grains spike-1 (rg = 0.41**) and 1000-grain weight (rg = 0.30**). Positive relationship of these aforementioned traits with grain yield indicated that these traits had major contribution towards grain yield in wheat (Table 5).Therefore, these traits needs due weightage in selecting wheat genotypes for higher grain yield.

Table 5: Genetic correlation coefficient among yield and yield components of 81 wheat genotypes across nine environments.

Pooled analysis of variance revealed significant differences among genotypes, environment and their interaction for tillers m-2, grains spike-1, 1000-grain weight and grain yield. Previously, Elhani et al. (2007), Ali et al. (2008), Shankarrao et al. (2010), Mohammadi et al. (2012), Mehari et al. (2015), and Ebrahimnejad and Rameeh (2016) also reported similar results in wheat. Contrarily, Khan et al. (2010) and Motamedi et al. (2013) reported non-significant differences among genotypes, environments and their interaction for tillers m-2 and for 1000-grain weight in wheat.Contradictory findings in this regard may be due to difference in genetic makeup of genotypes tested, environmental conditions or both.

In every plant breeding program the integral component is to develop high yielding line/cultivar which determines future of the crop and its growers (Muflin, 2000). For plant breeder there are many challenges to get besides, maintain high productivity and effective practice is to subjecting potential lines to a series of

diverse environmental conditions to identify the best stable yielding genotypes (Kaya et al., 2002; Roozeboom et al., 2004; Loffler et al., 2005). Favorable conditions during flowering and pollination result in higher number of fertile tillers-2 and grains spike-1(Agoston and Pepo, 2005). Tillers m-2, grains spike-1 and 1000-grain weight are very important and direct contributing traits towards grain yield in wheat, therefore stability of these traits leads to grain yield stability (Dreccer et al., 2008). G-79 produced maximum number oftillers m-2and grains spike-1 in almost all environments. Through tillers m-2, grains spike-1and 1000-grain weight considerable genetic progress for higher yield was achieved (Sayre et al., 1997).

Correlation analysis revealed that grain yield had positive associations with tillers plant-1, number of grain spike-1and 1000-grain weight. Mohsen et al. (2012) also reported that grain yield had positive correlations with above mentioned traits and suggested that plant breeders should consider these traits in breeding wheat for grain yield improvement.

Conclusions

Significant genetic differences among genotypes for yield and yield components exhibited the existence of sufficient variability to have an effective selection. Similarly, significant genotype by environment interactions (GEI) indicated inconsistent performance of genotypes across environments for almost all traits. The GEI captured major part of sum of squares exhibiting its greater effects in the phenotypic expression of all traits. Correlation analysis revealed significant positive associations of grain yield with tillers m-2, grains spike-1 and 1000-grain weight indicating that these traits had major contribution towards grain yield in wheat. Genotype G-79 exhibited maximum tillers m-2, grains spike-1 and grain yield across all environments, and was thus identified as high yielding genotype for commercialization in Khyber Pakhtunkhwa.

Author’s Contributions

This research article is an integral part of doctoral study of the first author MAUK. MAUK and FM) formulated the research. MAUK performed the experiments and wrote the first draft of the manuscript. FM made critical corrections in the first draft. Both authors read and approved the final manuscript.

References

Agoston, T. and P. Pepo. 2005. Effect of genetic and ecological factors on yield formation in winter wheat production. Cer. Res. Commun. 33 (1): 37-40. https://doi.org/10.1556/CRC.33.2005.1.8

Akcura. M., Y. Kaya and S. Taner. 2009. Evaluation of durum wheat genotypes using parametric and nonparametric stability statistics. Turk. J. Field Crop. 14 (2): 111-122.

Ali, W., K. Zada, P. Shah, A. Haider, S. Ayaz and S. Malik. 2008. Correlation and regression analysis of agronomic parameters with grain yield of wheat. Sarhad J. Agric. 8(4): 329-332.

Annicchiarico. P. 2002. Genotype × environment interaction: challenges and opportunities for plant breeding and cultivar recommendations, Food and Agriculture Organization of the United Nations.

Dreccer, M.F., S.C. Chapman, F.C. Ogbonnaya, M.G. Borgognone and R.M. Trethowan. 2008. Crop and environmental attributes underpinning genotype by environment interaction in synthetic-derived bread wheat evaluated in Mexico and Australia. Aust. J. Agric. Res. 59: 447-460. https://doi.org/10.1071/AR07220

Ebrahimnejad, S. and V. Rameeh. 2016. Correlation and factor analysis of grain yield and some important component characters in spring bread wheat genotypes. Cer. Agron. 1 (165): 5-15. https://doi.org/10.1515/cerce-2016-0001

Elhani, S., V. Martos, Y. Rharabatti, C. Royo and L.F. Garcia del Moral. 2007. Contribution of main stem and tillers to durum wheat (Triticum turgidum L var durum) grain yield and its components in Mediterranean environments. Field. Crop. Res. 103: 25-35. https://doi.org/10.1016/j.fcr.2007.05.008

Karimizadeh. R., M. Mohammadi, M.K. Shefazadeh, A.A. Mahmoodi, B. Rostami and F. Karimpour. 2012b.

Kaya, Y., C. Palta and S. Taner. 2002. Additive main effects and multiplicative interaction analysis of yield performances in bread wheat genotypes across environments. Tur. J. Agric. For. 26: 275-279.

Khan N., S. Syeed, A. Masood, R. Nazar, N. Iqbal. 2010. Application of salicylic acid increases contents of nutrients and antioxidative metabolism in mungbean and alleviates adverse effects of salinity stress. Int. J. Plant. Biol. 1:1- 10. https://doi.org/10.4081/pb.2010.e1

Loffler, C.M., J. Wei, T. Fast, J. Gogerty, S. Langton, M. Bergman, B. Merrill, and M. Cooper. 2005. Classification of maize environments using crop simulation and geographic information systems. Crop Sci. 45: 1708-1716. https://doi.org/10.2135/cropsci2004.0370

Mehari, M., M. Tesfay, H. Yirga, A. Mesele, T. Abebe, A. Workineh and B. Amare. 2015. GGE biplot analysis of genotype-by-environment interaction and grain yield stability of bread wheat genotypes in South Tigray, Ethiopia. Commun. Biomet. Cr. Sci. 10 (1): 17-26.

Mohammad. F., O.S. Abdalla, S. Rajaram, A. Yaljarouka, N.U. Khan, A.Z. Khan, S.K. Khalil, I.H. Khalil, I. Ahmad and S.A. Jadoon. 2011. Additive main effect and multiplicative analysis of synthetic-derived wheat under varying moisture regimes. Pak. J. Bot. 43 (2): 1205-1210.

Mohammadi. M., R. Karimizadeh, N. Sabaghnia and M.K. Shefazadeh. 2012. Genotype × environment interaction and yield stability analysis of new improved bread wheat genotypes. Turk. J. Field Crop. 17 (1): 67-73.

Mohsen, A., A.S.R. Hegazy and M.H. Taha. 2012. Genotypic and phenotypic interrelationships among yield and yield components in Egyptian bread wheat genotypes. J. Pl br. Crop. Sci. 4(1): 9-16.

Motamedi, M., P. Safari and C. Motamamedi. 2013. Additive main effect and multiplicative interaction analysis of grain yield in bread wheat genotypes across environments. Int. J. Biol. Sci. 3(8): 218-225.

Muflin, B. 2000. Crop improvement 21st century. J. Exp. Bot. 51(342 ): 1-8.

Pakistan Bureau of Statistics. 2014-15. Pakistan economic survey government of Pakistan, finance division, economic advisory wing, Islamabad.

Roozeboom, K.A., J .Fritz. J .Stack, M. Whiteworth, W.B. Classen, W. Gordon, K. Heer, L. Jansen, V. Maddax, P. Martin, J. Evans, T.J. Long, A. Martin, A. Schlegel and M. Witt. 2004. Kanas performance tests with winter wheat varieties. Kansas Agricultural Experiment Station Report of Progress 9030.Kansas Satate Univ., Manhattan.

SAS Institute Inc. 2009. SAS/STAT ® 9.2 User’s Guide, Second Edition. Copyright © 2009, SAS. Inst. Inc. Cary, NC, USA. (https://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/vie wer.htm)

Sayre, K.D., S. Rajaram and R.A. Ficher. 1997. Wheat yield potential in short bread wheat in North Mexico. Crop. Sci. 37:36-42. https://doi.org/10.2135/cropsci1997.0011183X003700010006x

Shankarrao, S.B., S. Mukherjee, A.K. Pal and D.K. De. 2010. Estimation of variability for yield parameters in bread wheat (Triticum aestivum. L) grown in Gangetic West Bengal. Elec. J. Pl. Br. 1 (4): 764-768.

Singh, P. K. and R. D. Choudhary. 1997. Biometrical methods in quantitative genetic analysis. Kalayani publishers, New Delhi. 178-185.