Grain Yield of Dual-Purpose Barley Genotypes Evaluated by AMMI, BLUP and Non Parametric Measures

Grain Yield of Dual-Purpose Barley Genotypes Evaluated by AMMI, BLUP and Non Parametric Measures

Ajay Verma^* , RPS Verma , J Singh , Lokendra Kumar , Gyanendra Pratap Singh

^*ICAR-Indian Institute of Wheat & Barley Research, Post Bag # 158 AgrasainMarg, Karnal 132001 (Haryana), India

Corresponding Author Email: verma.dwr@gmail.com

DOI : http://dx.doi.org/10.53709/ CHE.2022.v03i01.005

Abstract

Keywords

Download this article as:

Introduction

The nutritious rich green fodder for the livestock and acceptable quality grain for human consumption have been provided by the dual purpose barley.In the recent past, dual purpose barley has been widely utilized as an alternative source of green forage in the arid and semi-arid regions especially owing to the scarcity of green forage availability for the animals. In terms of the world’s most essential crops by production quantity, barley is ranked fourth amongst the cereals after maize, rice and wheat, although eleventh overall, and is widely grown across the world. Barley’s ability to adapt to multiple biotic and abiotic stresses would be desirable for future popularity particularly in the scenario of environmental change. Reliable information about cultivar performance in different environments has been provided by genotype × environment interaction effects with the execution of multi-environment trials [2],[1]. The good number of precise analytic approaches t has been advocated in literature [8]. Visibilities of AMMI based measures (AMMI stability value (ASV), ASV1, Modified AMMI stability value (MASV) & MASV1) have been observed [11].  Best linear unbiased prediction (BLUP) based measures harmonic mean of genotypic values (HMGV), the relative performance of genotypic values (RPGV), and harmonic mean of relative performance of genotypic values (HMRPGV), were also exploited for the stability and adaptability of genotypes [6],[7]. Nonparametric measures have been also utilized S_i¹,S_i²,S_i³,S_i⁴ ,S_i⁵ ,S_i⁶ ,S_i⁷ , NP_i ⁽¹⁾, NP_i ⁽²⁾, NP ⁽³⁾, NP_i ⁽⁴⁾ to interpret the response of genotypes to environmental conditions [8]. Recent analytic measures have been compared to decipher the GxE interactions effects for dual purpose barley genotypes evaluated in northern hills zone of the country.

Materials and Methods

Twenty three promising dual-purpose barley genotypes were evaluated in research field trials at 10centers of the All India Coordinated Research Project across northern hills zone of the country during 2020-21 cropping season. More emphasis had been placed to increase the dual purpose barley production of this zone to augment the fodder and grain production of the country. Field trials were laid out in Randomized block designs with four replications. Recommended practices of packages had followed in total to harvest the good yield. Parentage details and environmental conditions were reflected in table 1 for ready reference.  Various non parametric and parametric measures had been recommended for assessing GxE interaction and stability analysis [8]. For a two-way dataset with k genotypes and n environments X_ij de­notes the phenotypic value of ith genotype in jth environ­ment  where i=1,2, …k, ,j =, 1,2 ,…,n and r_ij  as the rank of the ith genotype in the jth environment, and  as the mean rank across all environments for the ith geno­type. The correction for yield of ith genotype in jth environment as (X*_ij=  X_ij– .+  ) as X*_ij, was the corrected phenotypic value; .was the mean of ith  genotype in all environments and was the grand mean.

	=

Non parametric composite measures NP_i⁽¹⁾, NP_i⁽²⁾, NP_i⁽³⁾ and NP_i⁽⁴⁾ based on the ranks of genotypes as per yield and corrected yield of genotypes. In the formulas, r^*_ij was the rank of X^*_ij, and  and M_di was the mean and median ranks for original (unadjusted) grain yield, where ^* and M^*_di were the same parameters computed from the corrected (adjusted) data.

½

ASV	ASV = [
ASV1	ASV1 = [
Modified AMMI stability Value
HMGVi	=  Number of environments /  genetic value of ith genotype in jth environments
Relative performance of genotypic values across environments	RPGVij =  /
Harmonic mean of Relative performance of genotypic values	HMRPGVi. =  Number of environments /
Geometric Adaptability Index	GAI =

AMMISOFT version 1.0 software utilized for AMMI analysis of data sets and SAS software version 9.3 for further analysis.

Results and Discussion

AMMI analysis

Highly significant variations due to environments, GxE interactions, and genotypes were observed by AMMI analysis (Table 2). This analysis also revealed about 71.8% of the total sum square of variation for yield was due to environments followed by 21.5%  of GxE interactions,  whereas genotypes accounted for marginally 3.3% [3]. Interaction effects are further portioned into  seven Interaction principal components totalled more than 98.6% of interactions sum of square variations. AMMI1 explained a total variation of 34.3%, followed by 28.1% for AMMI2, 11.5% for AMMI3, AMMI4 accounted for 8.3% and followed by 6.6%, 4.5% & 3.3% respectively. The first two AMMI components in total showed 62.8% of the total variation indicating the two AMMI components well fit and confirm the use of AMMI model [7]. Estimated sums of squares for G×E signal and noise were 88.5%  and 11.4% of total G×E respectively.  Early IPCs selectively capture signal, and late ones noise[4]. Note that the sum of squares for GxE-signal is 5.83 times that for genotypes main effects. Hence, narrow adaptations are important for this dataset [12]. Even just IPC1 alone is 2.26 times the genotype‘s main effects. Also, note that GxE-noise is 0.76 times the genotypes effects. Discarding noise improves accuracy, increases repeatability, simplifies conclusions, and accelerates progress.

Ranking of genotypes as per on AMMI based measures

Since the genotypes yield, the mean yield was considered an important measure to assess the genotypes. Mean yield of genotypes expressed highly significant variations and yield potential selected G5, G9, G14 and the lowest yield for G23, G17 (Table 3). Stability or adaptability of genotypes judged by values of IPCA’s in the AMMI analysis. The, greater the IPCA scores reflect the specific adaptation of genotype to certain locations. While, the values approximate to zero were recommended for in general adaptations of the genotype [9].  Absolute IPCA-1 scores pointed for G17, G16, G11 as per IPCA-2, G2, G8, G18  genotypes would be of choice. Values of IPCA-3 favored G23,G19, G20 genotypes. As per IPCA-4, G20, G17 G9 genotypes would be of stable performance. As per IPCA5, G9, G15, G21 while IPCA6 favoured G9, G21, G22 whereas IPCA7 settled for G19,G6, G20 First two IPCAs in ASV & ASV1 measures utilized 68.2% of G×E interaction sum of squares. The two IPCAs have different values and meanings and the ASV and ASV1 parameters using the Pythagoras theorem and to get estimated values between IPCA1 and IPCA2 scores to produce a balanced measure between the two IPCA scores. Also, ASV parameter of this investigation used advantages of cross validation due to computation from first two IPCAs [10]. Using first two IPCAs in stability analysis could benefits dynamic concept of stability in the identification of the stable high yielder genotypes. ASV1 measures recommended (G20, G5, G19) and ASV pointed towards (G20 G5 G19) as of stable performance. Adaptability measures MASV and MASV1considered all seven significant IPCAs of the AMMI analysis using 98.6% of GxE interactions sum of squares [5]. Values of MASV1 identified G19, G22, G20 genotypes would express stable yield whereas genotypes G19, G22, G20 be of stable yield performance by MASV measure respectively.

Ranking of genotypes as per BLUP based & Non parametric measures

Major advantages of BLUP based measures are to account for the random nature of the genotype behavior in changes climatic conditions. At the same time allow ranking genotypes in relation to their performance based on the genetic effects [11]. An average yield of genotypes pointed towards G5, G9,G14 as high yielders. The consistent yield of G18, G11,G17  as per least values of standard deviation more over the values of CV identified G18, G11, G20 genotypes for the consistent yield performance for northern hills zone of the country. More over the values of GM favoured G18,G5,and G9. The BLUP-based simultaneous selections, such as HMGV identified G18, G9,G10 whilevalues of RPGV favored G18, G5, G9 and HMRPGV settled G18, G9, G5 genotypes. The evaluation of adaptability and stability of wheat genotypes through these BLUP-based indices was reported by Pour-Aboughadareh [8]. The estimates of HMGV, RPGV, and HMRPGV had the same genotype ranking that was reported Anuradha [2].

Non parametric measures ranked the genotypes as per their corrected yield across environments S_i¹ values pointed for G19,G15,G17 while S_i² selected G19,G15,G20  and values of S_i³favoured G19,G20,G15  as desirable genotypes (Table 4). G19,G20,G15 selected by values of S_i⁴, S_i⁵ for G1,9 G20,G15, S_i⁶ forG19,G20,G9  and lastly S_i⁷ for G19, G20,G15 (Table 4). The mentioned strategy determines the stability of genotype over environment if its rank is similar over other environments (biological concept). Nonparametric measures of phenotypic stability were associated with the biological concept of stability [12]. Non parametric composite measures NP_i⁽¹⁾  to NP_i⁽⁴⁾, consider the ranks of genotypes as per their yield and corrected yield across environments simultaneously [9]. NP_i ⁽¹⁾ measure the observed suitability ofG19,G20,G15 whereas as per NP_i⁽²⁾,   genotypes G15,G19,G20 would be of choice while NP_i⁽³⁾  identified G15,G19,G17. Last composite measure NP_i⁽⁴⁾  found G15,G19,G20 as genotypes of choice for this zone.

Biplot analysis

The first two significant PC’s has explained about 60.1% of the total variation in the AMMI, BLUP and non parametric measures considered for this study in biplot analysis (Table 5) with respective contributions of 37.8% &22.2% by first and second principal components respectively [1]. Measures S_i¹, S_i², S_i³, S_i⁴, S_i⁵,S_i⁶ ,S_i⁷,NP_i ⁽¹⁾, NP_i ^(2), NP_i ⁽³⁾, NP_i ^(4), MASV, MASV1 accounted more of share in first principal component whereas Mean, Average, HMPRVG, PRVG, IPC5, GM, HM were major contributors  in PC2. The association analysis among measures had been explored with the biplot analysis.  In the biplot vectors of measures expressed acute angles would be positively correlated whereas those achieved obtuse or straight line angles would be negatively correlated. Independent types of relationships had expressed by right angles between vectors. Very tight positive relationships observed for NP_i⁽²⁾ , NP_i⁽³⁾ , NP_i⁽⁴⁾with Stdev, IPC7, IPC4.BLUP based measures GM, HM, PRVG, HMPRVG, Averagealso expressed tight relationships among themselves in the same quadrant.AMMI based measures MASV, MASV1 closely associated with Stdev, S_i² , S_i³ S_i⁵ , S_i^{6 ,} ASV  ASV1 values ^. BLCV observed withS_i¹ , S_i⁴ S_i⁷. This group of measures expressed no relationship with group consisted of BLCV ,S_i¹ , S_i⁴ S_i⁷ .BLCV alsoexpressed no relationships IPC3, NP_i⁽³⁾, NP_i⁽⁴⁾. Group  MASV, MASV1 also expressed right angles with BLUP based measures (Fig. 1).  Seven clusters of small and moderate sizes observed in biplot analysis. AMMI analysis based measures IPC1 & IPC5 formed first cluster as placed in first quadrant. Three clusters were observed in second quadrant. BLUP based measures GM, HM, PRVG, HMPRVG grouped with average and adjacent cluster comprised of NP_i⁽²⁾, NP_i⁽³⁾, NP_i⁽⁴⁾ whileStdevgrouped with IPC4, IPC7 in smaller cluster. Next quadrant occupied by two AMMI based measures i.e. one of ASV, ASV1 with MASV1 whereas other consisted of NP_i⁽¹⁾, S_i² ,S_i⁵S_i² ,S_i⁵ ,S_i⁶ with MASV Last cluster of CV with IPC2 & IPC6 placed in last quadrant (Fig. 2).

Acknowledgements

The training by Dr. J Crossa and financial support by Dr. A.K Joshi & Dr. RP Singh CIMMYT, Mexico sincerely acknowledged along with the hard work of the staff to carry out the field evaluation of genotypes at coordinating centers.

Conflict of Interests

No conflict of interests reported by the authors

References

[1]Ahakpaz F., Abdi H., Neyestani E., Hesami A., Mohammadi B., Nader Mahmoudi K., Abedi-Asl G., JazayeriNoshabadi M.R., Ahakpaz F., Alipour H .(2021). Genotype-by-environment interaction analysis for grain yield of barley genotypes under dry land conditions and the role of monthly rainfall.Agric Water Manag 245:10665

[2]Anuradha N., Patro T.S.S.K., Singamsetti A., Sandhya Rani Y., Triveni U., NirmalaKumari A., Govanakoppa N., Lakshmi Pathy T. and Tonapi V.A .(2022). Comparative Study of AMMI- and BLUP-Based Simultaneous Selection for Grain Yield and Stability of Finger Millet [Eleusinecoracana (L.)Gaertn.]Genotypes. Front. Plant Sci. 12:786839. DOI: 10.3389/fpls.2021.786839

[3]Bocianowski J., Tratwal A., NowosadK .(2021). Genotype by environment interaction for main winter triticale varieties characteristics at two levels of technology using additive main effects and multiplicative interaction model.Euphytica 217:26

[4]George N., Lundy M .(2019). Quantifying genotype x Environment effects in long-term common wheat yield trials from an agroecologically diverse production region. Crop Science 59:1960–1972

[5]Gerrano A.S., Rensburg W.S.J.V., Mathew I., Shayanowako A.I.T., Bairu M.W., Venter S.L., Swart W., Mofokeng A., Mellem J., Labuschagne M.(2020). Genotype and genotype x environment interaction effects on the grain yield performance of cowpea genotypes in the dry land farming system in South Africa. Euphytica 216:80

[6]Gonçalves G. de M. C. , Gomes R. L. F., Lopes Â. C. de A. and Vieira P. Fe. de M. J. (2020). Adaptability and yield stability of soybean genotypes by REML/BLUP and GGE Biplot. Crop Breeding and Applied Biotechnology, 20(2): e282920217.

[7]Pour‑Aboughadareh, A.  , Ali B., Ali K. S., Mehdi J., Akbar M., Ahmad G., Kamal S.H., Hassan Z., PoodinehOmid  andMasoome, K. (2022).  Dissection of genotype‑by‑environment interaction and yield stability analysis in barley using AMMI model and stability statistics. Bulletin of the National Research Centre 46:19

[8]Pour-Aboughadareh, A., Yousefian M., Moradkhani H.,  Poczai P., and Siddique K.H. (2019). STABILITYSOFT: A  new online program to calculate parametric and non-  parametric  stability statistics for crop traits.  Applications in Plant Sciences 7(1): e1211

[9]Mehraban R. A., Hossein-Pour T., Koohkan E., Ghasemi S., Moradkhani H., Siddique K.H .(2019). Integrating different stability models to investigate genotype × environment interactions and identify stable and high-yielding barley genotypes. Euphytica 215:63

[10]Silva E. M. da, Nunes E. W. L. P., Costa J. M. da, Ricarte A. de O., Nunes G. H. de S. and Aragão Fernando Antonio Souza de .(2019). Genotype x environment interaction, adaptability and stability of ‘Piel de Sapo’ melon hybrids through mixed models Crop Breeding and Applied Biotechnology, 19(4): 402-411.

[11]Sousa, A.M.C.B., Silva V.B., Lopes A.C.A., Ferreira-Gomes R.L. and CarvalhoL.C.B .(2020). Prediction of grain yield, adaptability, and stability in landrace varieties of lima bean (Phaseoluslunatus L.) Crop Breeding and Applied Biotechnology, 20: e295120115

[12]Vaezi, B., A. Pour-Aboughadareh, A. Mehraban, T. Hossein-Pour, Mohammadi R.  Armion,  M. and Dorri, M. (2018). The use of parametric and  non-  parametric measures for selecting stable and adapted barley lines.  Archives of Agronomy and Soil Science, 64: 597–611

Table 1: Parentage vis-a-vis location details of dual purpose barley genotypes

Code	Genotype	Parentage	Locations	Latitude	Longitude	Altitude
G1	HBL873	P.STO/3/LBIRAN/UNA80//LIGNEE640/4/BLLU/5/PETUNIA1/6/ P.STO/3/ LBIRAN/UNA80//LIGNEE640/4/BLLU/ 5/PETUNIA 1 (6th GSBON-2018-19-Ent 87)	Almora	29° 35 ‘ N	79° 39 ‘E	1610
G2	HBL870	VLB 118 x HBL 712	Berthein	31°50’N	77°9’E	1103
G3	VLB170	VB 1709 INBYT-HI (2016)-12 (CHAMICO/TOCTE //CONGONA/3/PETUNIA 2/4/ PENCO/CHEVRON-BAR)	Kangra
G4	BHS483	BHS352/BHS366	Katrain
G5	UPB1093	RD2784/RD2035	Majhera	29° 16′ N	80° 5′ E	1532
G6	VLB118	14th EMBSN-9313	Khudwani	33° 70′ N	75°10′ E	1590
G7	BHS487	BBM593/ BHS169	Malan	32°08 ‘ N	76°35’E	846
G8	BHS400	34th IBON-9009	Rajauri	31.01	75.92
G9	BHS486	HBL276/BHS365	Ranichauri	28° 43′ N	81°02′ E	2200
G10	VLB173	P.STO/3/LBIRAN/UNA80//LIGNEE640/4/BLLU/5/ PETUNIA1/6/ GLORIA- BAR/COPAL (IBON-HI-18-91)	Shimla	31°10 ‘ N	77°17’E	‎2276
G11	BHS352	HBL240/BHS504//VLB129
G12	HBL869	DWR 81 x BH 936
G13	VLB172	ZIGZIG/3/PENCO/CHEVRON-BAR//PETUNIA 1 (INBYT-HI-15-16-20)
G14	HBL113	SELECTION FROM ZYPHYZE
G15	BHS485	HBL276/BHS369
G16	BHS484	BHS352/BHS 169
G17	HBL872	P.STO/3/LBIRAN/UNA80//LIGNEE640/4/BLLU/5/PETUNIA1/6/P.STO/3/LBIRAN/ UNA80//LIGNEE640/4/BLLU/5/PETUNIA 1 (6th GSBON-2018-19 -Ent 86)
G18	UPB1092	RD2828/K551
G19	VLB171	BISON 110.3//CANELA/ZHEDAR#2 (IBON-HI-18-36)
G20	HBL871	TRADITION/6/VMorales/7/LEGACY//PENCO/CHEVRON-BAR  (IBON 16-17-Ent72 or EIBGN 2017-18, Ent-49)
G21	BHS380	VOILET/MJA/7/ABN-B6/BA/GAL// FZA-B /5/DG/DC-B/ PT-BAR /3/RA-B/ BA /3/4/TRYIGAL…
G22	VLB174	LIMON/BICHY2000//DEFRA/DESCONOCIDA-BAR (IBON-HI-18-83)
G23	UPB1091	RD2828/RD2552

Table 2: AMMI analysis of dualpurpose barley genotypes evaluated under coordinated trials

Source	Degree of  freedom	Mean Sum of Squares	Significance level	Proportional contribution of factors	GxE interaction Sum of Squares (% )	Cumulative Sum of Squares (% ) by IPCA’s
Treatments	229	202.1488	***	96.69
Genotype (G)	22	71.30654	***	3.28
Environment ( E )	9	3821.847	***	71.85
GxE interactions	198	52.15502	***	21.57
IPC1	30	118.0639	***		34.30	34.3
IPC2	28	103.7608	***		28.13	62.43
IPC3	26	45.89658	***		11.56	73.99
IPC4	24	36.0393	***		8.38	82.37
IPC5	22	31.15171	***		6.64	89.00
IPC6	20	23.69014	***		4.59	93.59
IPC7	18	19.23269	***		3.35	96.94
Residual	30	10.5298	*
Error	230	6.8866
Total	459	104.305

Table 3: AMMI based measures of dual purpose barley genotypes 

Genotype	Mean	IPC1	IPC2	IPC3	IPC4	IPC5	IPC6	IPC7	ASV1	ASV	MASV1	MASV	Average	Stdev	CV
G1	24.88	1.0684	-1.1416	0.7808	0.4632	-1.2261	1.4952	0.5252	1.73	1.64	4.940	4.130	24.89	10.10	40.58
G2	23.74	1.4543	-0.0193	-0.4993	-0.5389	1.1851	-1.1202	-0.6250	1.77	1.61	3.600	3.260	23.60	9.64	40.86
G3	23.75	-0.4967	-1.2817	1.6095	0.9338	0.7936	-0.8146	0.2235	1.42	1.39	5.043	4.152	23.74	9.99	42.08
G4	23.29	2.1107	1.5634	0.2267	0.9424	-0.8746	-1.1796	0.8676	3.01	2.81	5.755	4.677	23.33	9.42	40.37
G5	26.57	0.5222	-0.6414	-1.1806	0.3345	-0.9192	-2.3898	0.2889	0.90	0.86	5.167	4.582	26.17	11.82	45.16
G6	23.56	-1.3447	0.3324	-1.6548	0.5188	1.5086	0.3013	0.1003	1.67	1.52	4.406	3.937	23.36	9.43	40.36
G7	21.16	1.9748	-0.9179	-0.7627	1.1159	0.3093	0.8837	-1.1796	2.58	2.37	4.526	3.915	21.38	9.53	44.56
G8	24.47	-2.0286	0.0549	-1.7048	-1.5721	-0.3126	-0.1843	0.5930	2.47	2.24	4.662	4.267	24.17	12.57	52.01
G9	25.56	1.8762	-0.8786	0.7747	-0.1855	-0.0001	0.0279	-0.3911	2.45	2.25	3.545	2.933	25.24	10.11	40.03
G10	23.62	-1.9996	1.3721	1.8523	1.5860	-0.4011	-0.2855	-0.4764	2.80	2.60	6.035	5.100	23.50	9.44	40.17
G11	21.71	0.2550	1.5066	1.4432	-2.1051	1.3184	0.1556	-0.4427	1.54	1.53	6.245	5.231	21.94	8.43	38.44
G12	23.24	0.6587	1.5720	-1.5094	-0.3315	-0.2813	0.7206	0.3597	1.77	1.73	5.151	4.030	23.22	9.09	39.15
G13	24.69	-2.1222	-2.9801	1.1983	-0.3100	0.1564	-0.3450	0.6427	3.95	3.79	8.571	6.355	24.61	13.81	56.11
G14	25.21	0.9586	-1.5315	0.3211	-2.0016	-1.0634	0.4593	0.3072	1.93	1.86	5.701	4.674	25.08	12.27	48.92
G15	20.23	0.8213	1.2798	0.7231	-1.0888	-0.0155	-0.6287	-0.4576	1.62	1.57	4.276	3.394	20.64	9.05	43.84
G16	25.36	-0.1706	-1.6502	-1.4259	0.7134	-0.7127	0.1366	-0.8037	1.66	1.66	5.329	4.161	25.08	11.55	46.06
G17	19.99	0.1713	2.2338	0.6881	0.0856	-0.9129	0.4785	0.4482	2.24	2.24	6.280	4.593	20.44	8.56	41.88
G18	24.78	1.4385	0.0697	-0.2428	1.0582	1.7025	0.8191	1.6556	1.76	1.59	4.451	4.078	24.67	7.67	31.08
G19	23.61	-0.3015	-0.8463	-0.0720	-0.4296	-0.6738	0.3059	0.0852	0.92	0.91	2.696	2.083	23.62	10.59	44.84
G20	22.94	-0.2757	-0.6192	0.1625	-0.0197	1.3053	0.4004	-0.2531	0.70	0.69	2.939	2.467	22.93	8.91	38.85
G21	20.73	-2.1470	1.0305	-0.3826	0.2324	0.0886	0.0730	1.0735	2.81	2.58	3.996	3.305	21.15	10.13	47.90
G22	22.86	-1.0207	0.5192	-0.3268	0.2930	0.1763	-0.1352	-1.7111	1.35	1.24	2.651	2.386	22.81	9.27	40.64
G23	19.89	-1.4029	0.9732	-0.0186	0.3055	-1.1507	0.8258	-0.8304	1.97	1.83	4.059	3.378	20.28	9.87	48.68

Table 4: BLUP based and Non parametric measures of dual purpose barley genotypes

	Genotype	GM	HM	PRVG	HMPRVG	Si1	Si2	Si3	Si4	Si5	Si6	Si7	NPi1	NPi2	NPi3	NPi4
HBL873	G1	22.58	19.82	1.0636	1.0291	8.69	52.68	4.28	7.26	6.04	4.91	38.54	5.70	0.6000	0.7331	0.8777
HBL870	G2	21.84	20.22	1.0235	1.0000	7.24	38.40	3.00	6.20	5.36	4.19	27.00	5.20	0.4727	0.5296	0.6192
VLB170	G3	22.04	20.60	1.0347	1.0088	8.20	48.99	3.52	7.00	6.10	4.39	31.72	6.10	0.4880	0.5555	0.6508
BHS483	G4	21.58	19.90	1.0228	0.9778	10.07	76.94	6.16	8.77	7.80	6.24	55.40	7.70	0.5133	0.7074	0.8118
UPB1093	G5	23.42	20.63	1.1011	1.0655	7.33	39.29	3.07	6.27	5.20	4.06	27.63	5.20	0.8667	0.8247	0.9649
VLB118	G6	21.78	20.41	1.0211	0.9978	7.51	40.00	3.33	6.32	5.00	4.17	30.00	5.00	0.4167	0.5406	0.6420
BHS487	G7	18.93	15.83	0.9003	0.8517	7.69	43.60	3.82	6.60	5.04	4.42	34.42	4.80	0.2824	0.4316	0.5025
BHS400	G8	20.99	17.47	0.9974	0.9449	9.20	61.16	4.93	7.82	6.80	5.48	44.39	6.80	0.6800	0.6982	0.8214
BHS486	G9	23.26	21.15	1.0875	1.0684	6.73	32.99	2.73	5.74	4.32	3.57	24.54	4.30	0.7167	0.7557	0.8860
VLB173	G10	22.02	20.78	1.0539	0.9803	9.64	73.96	7.55	8.60	7.56	7.71	67.92	6.80	1.2364	0.9149	1.0260
BHS352	G11	20.67	19.59	0.9760	0.9393	9.56	66.93	6.31	8.18	7.40	6.98	56.83	7.40	0.6167	0.6818	0.7963
HBL869	G12	21.55	19.88	1.0104	0.9862	7.33	39.16	3.69	6.26	4.92	4.64	33.25	4.80	0.5333	0.5587	0.6548
VLB172	G13	21.47	18.56	1.0178	0.9737	9.11	64.49	4.13	8.03	6.56	4.21	37.21	6.20	0.4133	0.6130	0.6955
HBL113	G14	21.36	16.20	1.0314	0.9313	8.67	52.49	3.92	7.24	5.80	4.33	35.25	5.80	0.5800	0.7393	0.8844
BHS485	G15	19.01	17.59	0.8876	0.8744	5.71	26.68	2.28	5.17	4.36	3.73	20.52	3.90	0.2229	0.3056	0.3379
BHS484	G16	21.56	16.54	1.0371	0.9449	8.29	52.77	4.43	7.26	6.28	5.28	39.91	5.70	0.8769	0.7489	0.8545
HBL872	G17	18.80	17.21	0.8878	0.8524	7.40	38.06	3.62	6.17	5.10	4.86	32.62	5.10	0.3091	0.3954	0.4744
UPB1092	G18	23.57	22.49	1.1133	1.0720	9.36	60.77	5.02	7.80	6.52	5.39	45.20	6.50	0.7222	0.7954	0.9546
VLB171	G19	21.09	18.12	0.9860	0.9661	4.29	13.83	1.11	3.72	3.00	2.40	9.96	2.70	0.2348	0.3263	0.3762
HBL871	G20	21.35	19.83	0.9941	0.9843	5.98	26.94	2.16	5.19	3.70	2.96	19.40	3.70	0.2846	0.3993	0.4598
BHS380	G21	19.23	17.67	0.9105	0.8683	8.13	47.60	4.10	6.90	6.00	5.17	36.93	6.00	0.3636	0.4539	0.5351
VLB174	G22	21.14	19.59	0.9895	0.9694	7.76	41.88	3.77	6.47	5.30	4.77	33.95	5.30	0.4818	0.5531	0.6629
UPB1091	G23	17.99	15.66	0.8526	0.8117	7.84	42.99	4.34	6.56	5.30	5.35	39.08	5.30	0.3533	0.4047	0.4842

Table 5:  Loadings of AMMI, BLUP and Non parametric measures

Measure	Principal Component 1	Principal Component  2	Measure	Principal Component 1	Principal Component  2
Mean	0.1258	-0.3182	GM	0.1295	-0.3446
IPC1	-0.0239	-0.1254	HM	0.0936	-0.2578
IPC2	-0.0077	0.1737	PRVG	0.1540	-0.3287
IPC3	0.0660	0.0701	HMPRVG	0.0980	-0.3506
IPC4	0.0497	-0.0368	Si1	0.2751	0.0950
IPC5	-0.0169	-0.0797	Si2	0.2808	0.0982
IPC6	-0.0505	0.0945	Si3	0.2585	0.1368
IPC7	0.1007	-0.0472	Si4	0.2798	0.1000
ASV1	0.1547	0.1940	Si5	0.2747	0.1073
ASV	0.1542	0.2004	Si6	0.2321	0.1624
MASV1	0.1983	0.1216	Si7	0.2585	0.1368
MASV	0.2295	0.0965	NPi1	0.2673	0.1033
Average	0.1274	-0.3162	NPi2	0.2225	-0.1363
Stdev	0.0616	-0.0308	NPi3	0.2552	-0.1625
CV	-0.0047	0.1493	NPi4	0.2514	-0.1680
60.05	37.86	22.20

PC1=37.8%;PC2=22.2%;TOTAL =60.1%

Figure 1: Biplot analysis of AMMI, BLUP and Non parametric measures

Figure 2: Clustering pattern of AMMI, BLUP and Non parametric measures