Volume 7, Issue 5, October Issue - 2019, Pages:468-476 |
Authors: Ajay Verma, J. Crossa, M. Vargas, A.K. Joshi, G.P. Singh |
Abstract: Best linear unbiased predictors of wheat genotypes were analysed as per analytic measures of adaptability for the North Eastern Plains zone of the country under mixed model framework. Among the tested genotypes, genotype K1006 had better adaptation while HD2733 & HD2967 were of specific adaptability during the cropping season of 2016-17. Biplot analysis based on first two significant principal components, accounted for 72.2 % of total GxE interaction sum of squares, observed stable performance of HD2967 and DBW187 genotypes while K1007 and HD2733 would be suitable for specific adaptations. Further, genotypes HD2967, DBW187 and HD3219 had expressed specific adaptation to Varanasi, Faizabad and Dumka, whereas K0307 and DBW39 were better adapted to conditions of Coochbehar & Chianki. Bikramganj with Manichak, Coochbehar with Chianki, Varanasi with Faizabad, Banka with Sabour and Bardwan maintained acute angles therefore would show similar performance of genotypes. Better adaptable genotypes with higher yield were HD2733, PBW762, HD3249 & DBW223 as per year cropping season of 2017-18. Biplot analysis identified HD2967, DBW233, HD3254 along with PBW762 and DBW39 genotypes of stable performance. Genotypes WH1218, DBW223, HD2733, HD3249 and DBW221 expressed unstable performance. Harmonic means of the relative performance of predicted genetic values by BLUP is suitable to identify the better adaptive genotypes with high yield. |
[Download PDF] |
Full Text: 1 Introduction Major objective of a wheat improvement program of any country is to identify promising genotypes which represent yield advantage relative to established check varieties (Crespo et al., 2017). Efficiency of a breeding program is judged by positive changes in yield performance over time (de Oliveira et al., 2017). Main concern in the selection process is to exploit G×E interaction effects so as to reduce the degree of uncertainty in the promotion process for broad or specific adaptability (Burgueño et al., 2011; Silva et al., 2014). Recently linear mixed model analysis of plant breeding experiments under multi environments trials has proved efficient and widely used (Friesen et al., 2016). Focus of MET analysis is to select the better performing genotypes, based on their yield rankings, owing to random performance of genotypes and consequent use of BLUPs is justified (Gogel et al., 2018). BLUPs maximize the correlation between the realized values and true values (Piepho et al., 2008). Full multivariate BLUP model utilizes all information along with variance heterogeneity (Kleinknecht et al., 2011). Factor analytic (FA) model with sufficient multiplicative terms is computationally robust (Smith & Cullis, 2018) and superiority of the Factor Analytic model in a breeding program had been demonstrated by Nuvunga et al. (2018). North eastern plains zone of India comprises eastern Uttar Pradesh, Bihar, Jharkhand, Assam and plains of West Bengal. Wheat is cultivated under highly diverse situations in around 8 million ha area. Among different wheat growing zones, this zone occupies 27% of total wheat area and accounts for 22% of the total wheat production in the country. The prime objective was to study adaptability of wheat genotypes by analytic measures under mixed models methodology. 2 Materials and Methods Seven advanced wheat genotypes were evaluated under multi environments trials at 14 major locations of the zone and fifteen advanced wheat genotypes were at thirteen major locations of the zone during 2016-17 and 2017-18 cropping seasons respectively (Figure 1). Field trials were laid out in Randomized Block Design with three replications and recommended agronomical interventions followed to harvest the good crop. More over the yield were analysed further to generate useful information. Details of genotypes and locations were reflected in tables 1 & 2 for ready reference. Mixed model to estimate fixed effects of blocks (b), random effects of genotypes (g) and interaction GxE effects (c) along with random errors (e) described as Y = Xb + Zg +Wc + e Where incidence matrices were X, Z and W respectively (Resende & Duarte, 2007). Simple and effective measure for adaptability is calculated as the relative performance of genetic values (PRVG) across environments (Resende, 2007). Moreover measure MHVG method (harmonic mean of genetic values), based on the harmonic mean of the genotypic values across in different environments (Resende, 2007). Lower the standard deviation of genotypic performance across environments, the greater is the harmonic mean of its genotypic values. Another measure based on the harmonic mean of the relative performance of the genotypic values (MHPRVG) under the mixed models analysis (Resende, 2007) for the simultaneous analysis of stability, adaptability and yield. Consequently, the selection for higher values of the harmonic mean results in selection for both yield and stability. PRVGij = VGij / VGi MHVGi = Number of environments / i=1k1Xi MHPRVGi. = n / j=1k1PRVGij VGij is the genotypic value of the i genotype, in the j environment, expressed as a proportion of the average in this environment. To facilitate interpretation, PRVG and MHPRVG values were multiplied by the general mean (GM), to provide results in the same magnitude as of the average wheat yield (Verardi et al., 2009). Estimation of the variance components were carried out by ASReml-R package using residual maximum likelihood (REML) along with estimation / prediction of the fixed as well as random effects (Cullis et al., 2014). 3 Results and discussion 3.1 First year of study (2016-17) Average yield of genotypes as per BLUP’s identified in K1006 & DBW39 as desirable while HD2967 along with HD2733 expressed low yields (Table 3). Harmonic mean of genotypes values selected K1006 & DBW39 as high yielders at the same lower yield of HD2733 & HD2967. Mean and Harmonic yield of studied genotypes based on BLUE’s pointed towards DBW187 and HD3219. PRVG as well as PRVG*GM pointed out K1006 & DBW39 as of better adaptable genotypes whereas HD2733 & HD2967 with low adaptability across locations of North Eastern Plains Zone of the country. Recent analytic measures HMPRVG and HMPRVG*GM marked K1006 & K0307 genotypes of high yield and better adaptability across zone. In general an agreement was observed in the classification of wheat genotypes, based on BLUP, MHVG, PRVG, MHPRVG and average yield (Table 3). Genotypic values can be predicted by measures based on Harmonic means that cares of adaptability along with yield (Verardi et al., 2009). Under mixed model methodology MHPRVG measures observed as an efficient for wheat improvement programs as mentioned by Mendes et al 2012 for other crops (Verardi et al., 2009; Mendes et al., 2012). Large variation in average yield of wheat genotypes had observed on the BLUE and BLUP across zone (Figure 2). More over consistent lower values had seen for genotypes based on BLUE except for HD3219 & DBW187. However the heights standard error of genotypes were at par with fixed and random effects assumptions of genotypes behavior. First two significant principal components accounted for 72.2 % of total GxE interaction sum of squares (Figure 4). Biplot analysis pointed out stable performance of genotypes or environments occupy positions near the cross section of coordinate system and vice versa for unstable behavior (Durate & Vencovsky, 1999). According to Yang & Kang (2003) specific adaptability of genotypes justified by close proximity near to environment. Genotypes HD2967 and DBW187 would be of stable performance while K1007 and HD2733 would be suitable for specific adaptations (Figure 4). This interpretation is based on distance of genotypes from origin of biplot analysis. Environments of Bardwan, Ranchi, Kanpur and Kaul would be the larger contributors to the G × E interaction of genotypes. On the other hand, Sabour and Shillongani environments would be suitable for stable performance of wheat genotypes. Genotypes and environments placed close to each other in the biplot have positive associations; this behavior would enable to identify specific zones for relatively similar performance of genotypes. HD2967, DBW187 and HD3219 had a specific adaptation to Varanasi, Faizabad and Dumka environments, whereas K0307 and DBW39 were better adapted genotypes to environments Coochbehar & Chianki. Genotype HD2733 expressed suitability to Ranchi & Ghazipur locations. |
Burgueño J, Crossa J, Cotes JM, Vicente FS, Das B (2011) Prediction assessment of linear mixed models for multi environment trials. Crop Science 51: 944-954. Crespo-Herrera LA, Crossa J, Huerta-Espino J, Autrique E, Mondal S, Velu G (2017) Genetic yield gains in CIMMYT’s International Elite Spring Wheat Yield Trials by modeling the genotype × environment interaction. Crop Science 57:789–801. Cullis BR, Jefferson P, Thompson R, Smith AB (2014) Factor analytic and reduced animal models for the investigation of additive genotype by environment interaction in out crossing plant species with application to a pinus radiata breeding program. Theoretical and Applied Genetics 127:2193–2210. Duarte JB, Vencovsky R (1999) Interação genótipos x ambientes: uma introdução à análise AMMI. Ribeirão Preto: Soc Bras Genét. Pp. 60. Friesen LF, Brule-Babel AL, Crow GH, Rothenburger PA (2016) Mixed model and stability analysis of spring wheat genotype yield evaluation data from Manitoba, Canada. Canadian Journal of Plant Science 96 : 305–320. Gogel BJ, Smith AB, Cullis BR (2018) Comparison of a one and two-stage mixed model analysis of Australia’s National Variety Trial Southern Region wheat data. Euphytica 214: 44-64. Kleinknecht K, Laidig F, Piepho HP, Möhring J (2011) Best linear unbiased prediction (BLUP): Is it beneficial in official variety performance trials?. Biuletyn Oceny Odmian 33:21–33. Mendes FF, Guimarães LJM, Souza JC, Guimarães PEO, Pacheco CAP, de Machado J RA, Meirelles W F, da Silva A R, Parentoni S (2012) Adaptability and stability of maize varieties using mixed model methodology. Crop Breeding and Applied Biotechnology 12: 111-117. Murakami DM, Cruz CD (2004) Proposal of methodologies for environment stratification and analysis of genotype adaptability. Crop Breeding and Applied Biotechnology 4:7-11. Nuvunga JJ, Oliveira LA, Silva CP, Pamplona AKA, Silva AQ, Moura EG, Maleia MP, Balestre M (2018) Adaptability and stability of cotton cultivars (Gossypium hirsutum L. race latifolium H.) using factor analytic model. Genetics and Molecular Research 17 :1-10. de Oliveira IJ, Atroch AL, Dias MC, Guimarães LJ, de Paulo EOG (2017) Selection of corn cultivars for yield, stability and adaptability in the state of Amazonas, Brazil. Pesquisa Agropecuária Brasileira 52:455-463. Piepho HP, Möhring J, Melchinger AE, Büchse A (2008) BLUP for phenotypic selection in plant breeding and variety testing. Euphytica 161: 209-228. Resende MDV (2007) Selegen-reml/blup: statistical system and computerized genetic selection by linear mixed models. Embrapa Florestas, Colombo, PR, Brazil. Resende MDV, Duarte JB (2007) Precision and quality control in variety trials. Pesquisa Agropecuária Tropical 37: 182-194. Silva P da, Bisognin DA, Locatelli AB, Storck L (2014) Adaptability and stability of corn hybrids grown for high grain yield. Acta Scientiarum. Agronomy 36:175-181. Smith AB, Cullis BR (2018) Plant breeding selection tools built on factor analytic mixed models for multi-environment trial data. Euphytica 214:143-161. Verardi CK, Resende MDZV, Costa RB, Gonçalves PS (2009) Adaptabilidade e estabilidade da produção de borracha e seleção em progênies de seringueira. Pesquisa Agropecuária Brasileira 44: 1277-1282. Yang W, Kang MS (2003) GGE Biplot Analysis: A Graphical Tool for Breeders, Geneticists, and Agronomists. CRC Press, Boca Raton, FL. |