Submit or Track your Manuscript LOG-IN

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





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.

Abstract | Knowledge of genotype by environment interaction (GEI) is essential for identifying environment-specific and widely adapted genotypes. Multi-environment trials were conducted to evaluate the magnitude of genotype (G), environment (E) and GEI effects on yield and yield components of wheat. Eighty-one wheat genotypes were evaluated under nine environments within Khyber Pakhtunkhwa, Pakistan during 2013/14, 2014/15 and 2015/16 cropping seasons. Combined analyses of variance for G and GEI were significant for almost all traits. Significant GEI suggested that the performance of genotypes was not consistent across environments, revealing strong influence of environments. Generally, GEI explained major portion of total variation and thus had larger effect than G and E on the expression of phenotype for all traits. Average over nine environments, genotype G-79 exhibited maximum values for tillers m-2 (181 tillers); grains spike-1 (71 grains); grain yield (4862 kg ha-1) and thus appeared as leading genotype for yield and production traits. Similarly, G-79 had higher grain yield in all test environments except E-03 i.e. E-01 (4537 kg ha-1); E-02 (4840 kg ha-1); E-04 (5035 kg ha-1); E-05 (4976 kg ha-1); E-06 (4797 kg ha-1); E-07 (5024 kg ha-1); E-08 (4767 kg ha-1); E-09 (4886 kg ha-1) signifying it as high yielding and widely adaptable genotype. However, G-08 with higher grain yield in E-01 (5272 kg ha-1) and E-07 (4840 kg ha-1), G-37 in E-04 (5106 kg ha-1) and E-08 (4906 kg ha-1), G-41 in E-05 (5034 kg ha-1) and E-06 (5272 kg ha-1), G-9 in E-02 (4975 kg ha-1) and G-49 in E-03 (5271 kg ha-1), established their specific adaptability in respective environments. Among environments, E-01, E-02, E-03, E-06 and E-07 were identified as highly productive environments in terms of grain yield. Correlation analysis revealed significant positive association of grain yield with tillers m-2 (rg = 0.72**), grains spike-1 (rg = 0.41**) and 1000-grain weight (rg = 0.30**). Based on mean yield performance, G-79 was found as high yielding genotype and thus could be recommended for commercialization in Khyber Pakhtunkhwa.

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:

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.


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


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.

Environments Growing season Geographical Positon Altitude (m.a.s.l) Average rainfall (mm) Temperature (0c)

Latitude Longitude Min. Max.
E-01 2014 (AUP) 34.0150° N

71.5805° E

359 238 20.1 34.8
E-02 2015 (AUP) --do-- --do-- --do-- 415 19.5 35.4
E-03 2016 (AUP) --do-- --do-- --do-- 189 17.8 38.2
E-04 2015 (CCRI) 34.0159° N

71.9755° E

288 220 10.1 28.6
E-05 2016 (CCRI) --do-- --do-- --do-- 112 16.3 35.9
E-06 2015 (ARSS) 34.1442° N

72.3785° E

321 263 18.0 36.7
E-07 2016 (ARSS) --do-- --do-- --do-- 312 14.5 32.1
E-08 2015 (ARSC) 34.1494° N

71.7428° E

381 460 10.4 28.5
E-09 2016 (ARSC) --do-- --do-- --do-- 392 17.4


Source: Meteorological Department, Khyber Pakhtunkhwa, Pakistan.

Breeding history of plant material

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.

Code Pedigree Code Pedigree
G-01 Takbir × Khatakwal-3-1 G-42 Tatara x Inqilab-26-7
G-02 Takbir × Khatakwal-3-5 G-43 Tatara × Inqilab-26-11
G-03 Takbir × Khatakwal-3-7 G-44 Tatara × Inqilab-26-15
G-04 Takbir × Khatakwal-3-8 G-45 Tatara × Inqilab-26-20
G-05 Takbir × Khatakwal-3-9 G-46 Tatara × Ghaznavi 98-31-1
G-06 Takbir × Khatakwal-3-16 G-47 Tatara × Ghaznavi 98-31-2
G-07 Takbir × Khatakwal-3-18 G-48 Tatara × Ghaznavi 98-31-4
G-08 Tatara × Inqilab-4-3 G-49 Tatara × Ghaznavi 98-31-7
G-09 Tatara × Inqilab-4-6 G-50 Ghaznavi 98 × Khatakwal -33-5
G-10 Tatara × Inqilab-4-9 G-51 Ghaznavi 98 × Khatakwal -33-7
G-11 Tatara × Inqilab-4-10 G-52 Ghaznavi 98 × Khatakwal -33-10
G-12 Tatara × Inqilab-4-11 G-53 Ghaznavi 98 × Khatakwal -33-15
G-13 Tatara × Inqilab-4-13 G-54 Tatara × Ghaznavi 98-37-15
G-14 Tatara × Inqilab-4-16 G-55 Tatara × Margala-43-2
G-15 Wafaq × Ghaznavi 98 G-56 Tatara × Margala-43-4
G-16 Wafaq × Ghaznavi 98 G-57 Tatara × Margala-43-11
G-17 Wafaq × Ghaznavi 98 G-58 Tatara × Margala-43-12
G-18 Tatara × Takbir-9-8 G-59 Tatara × Inqilab -45-10
G-19 Tatara × Takbir-9-10 G-60 Takbir × Inqilab -45-12
G-20 Tatara × Takbir-9-12 G-61 Tatara × Ghaznavi 98-48-2
G-21 Tatara × Takbir-9-813 G-62 Tatara × Ghaznavi 98-48-3
G-22 Tatara × Inqilab-18-15 G-63 Tatara × Ghaznavi 98-48-13
G-23 Tatara × Inqilab-18-19 G-64 Tatara × Ghaznavi 98-48-15
G-24 Tatara × Inqilab-18-20 G-65 Tatara × Ghaznavi 98-48-19
G-25 Tatara × Takbir-19-3 G-66 Wafaq × Ghaznavi 98-49-2
G-26 Tatara × Takbir-19-4 G-67 Wafaq × Ghaznavi 98-49-4
G-27 Tatara × Takbir-19-8 G-68 Wafaq × Ghaznavi 98-49-5
G-28 Tatara × Takbir-19-11 G-69 Wafaq × Ghaznavi 98-49-6
G-29 Tatara × Takbir-19-16 G-70 Wafaq × Ghaznavi 98-49-9
G-30 Tatara × Takbir-19-18 G-71 Wafaq × Ghaznavi 98-49-10
G-31 Tatara × Ghaznavi 98-22-1 G-72 Wafaq × Ghaznavi 98-49-12
G-32 Tatara × Ghaznavi 98-22-2 G-73 Wafaq × Ghaznavi 98-49-13
G-33 Tatara × Ghaznavi 98-22-6 G-74 Wafaq × Ghaznavi 98-49-15
G-34 Tatara × Ghaznavi 98-22-8 G-75 Wafaq × Ghaznavi 98-49-16
G-35 Tatara × Ghaznavi 98-22-9 G-76 Wafaq × Ghaznavi 98-49-19
G-36 Tatara × Ghaznavi 98-22-12 G-77 Wafaq × Ghaznavi 98-49-20
G-37 Tatara × Ghaznavi 98-22-13 G-78 Tatara × Takbir-19-17
G-38 Tatara × Ghaznavi 98-22-19 G-79 Tatara × Takbir-19-18
G-39 Tatara × Ghaznavi 98-22-20 Check Janbaz
G-40 Tatara × Inqilab-26-4 Check Atta-Habib
G-41 Tatara × Inqilab-26-6    

segregating for stripe rust resistance. Therefore, twenty heads were hand harvested and planted as head-to-row during 2012/13. The heavy infestation of stripe and leaf rusts fungi offered opportunity for screening F5:7 RILs. For the current study, 79 F5:8RILs were selected for further evaluation in multi-location trials over years.

Experimental design and procedure

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.

  Parameter Tiller-2

Grains spike-1

1000-grain weight Grain yield
Across nine environments Range 128-180 50.7-70.7 33.4-44.3 2804-4862
Desirable RIL G-79 G-79 G-72 G-79
Mean 153 60 38.3 3438
E-01 Range 91-229 42-83 29.1-51.7 1938-5272
Desirable RIL G-26 G-79 G-36 G-7,G-79
Mean 157 60 39.2 3597
E-02 Range 91-229 42-83 29.1-51.7 1630-4975
Desirable RIL G-79 G-51 G-08 G-9,G-79
Mean 161 60 39.2 3670
E-03 Range 77-224 46-85 25.1-65.9 1939-5071
Desirable RIL G-39,G-63 G-79 G-40 G-49
Mean 158 60 39.9 3620
E-04 Range 100-198 42-83 27.6-55.7 1931-5106
Desirable RIL G-06 G-79 G-14 G-37,G-79
Mean 148 60 37.2 3302
E-05 Range 76-220 43-85 28.0-56.0 1873-5035
Desirable RIL G-79 G-41 G-79


Mean 147 60 37.2 3284
E-06 Range 99-275 47-76 27.3-51.4 2023-4920
Desirable RIL G-31,G-41 G-78,G-79 G-46 G-41,G-79
Mean 163 61 37.4 3608
E-07 Range 94-229 47-76 31.4-49.8 1938-5024
Desirable RIL G-19,G-28 G-12,G-26 G-36,G-79 G-08,G-79
Mean 161 61 38.1 3615
E-08 Range 97-178 43-76 27.6-51.1 1931-4906
Desirable RIL G-10,G-36 G-25,G-79 G-37 G-37,G-79
Mean 137 59 37.9 3089
E-09 Range 92-210 47-73 32.9-50.5 2109-4886
Desirable RIL G-19,G-79 G-54 G-77,G-79 G-79
Mean 145 59 38.8


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.

Traits Tiller-2

Grains spike-1

1000-grain weight Grain yield



-0.11 ns

0.20 ns


Grains spike-1




1000-grain weight     -


Grain yield       -


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.


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.


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.

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.

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.

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.

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.

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.

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. ( 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.

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.

Sarhad Journal of Agriculture


Vol.39, Iss. 4, Pages 781-1002


Click here for more

Subscribe Today

Receive free updates on new articles, opportunities and benefits

Subscribe Unsubscribe