Document Type : Original Article

Authors

1 Department of Agronomy and Plant Breeding, Faculty of Agriculture, Lorestan University, Khorramabad, Iran

2 Assistant Professor. Department of Agronomy and Plant Breeding, Faculty of Agriculture, Lorestan University, Khorramabad, Iran

3 Research assistant, Seed and Plant Improvement Departement, Lorestan Agricultural and Natural Resources, Research and Education Center, AREEO, Khoramabad, Iran

Abstract

Introduction
The main objective in breeding programs is mainly based on evaluation experiments to select superior genotypes and accurate estimates of the mean of varieties of different varieties. Soil fertility, soil water-holding capacity, soil physical and mineral characteristics, wind direction and other unknown agents are often varied across an experimental site. This means that that the residuals of conventional models such as RCB are not independent of their near plots or spatial locations. So, this is an inherent issue which spatial adjustment of errors necessarily improves the precision of estimated treatment. The main objectives of this paper were to (1) explore the characterization of spatial dependence of traits using linear mixed models in large field chickpea trials in normal and drought conditions, (2) compare the efficiency of statistical analysis of spatial model to RCB model, (3) estimation of genetic parameters based on spatial models.
Materials and methods
This research was carried out in order to evaluate the yield and yield components of 64 chickpea genotypes in normal and rainfed conditions in Khorramabad Province, Iran. Experiments were conducted in an 8 × 8 Lattice square design framework in two replications for each condition during the 2018 season. Due to seasonal precipitation, there was no significant statistical difference between rainfed and normal conditions as well genotype by environment interaction in most of the traits. Therefore, rainfed data were considered as replications for the normal environment in all traits and all replications were subjected as one condition. To analyze the data for removing heterogeneity of farm condition in Lattices square design, Federer's spatial models were used to optimize the Lattice square linear model. Finally, the nearest neighbor model was selected as the best model based on the relative advantage and also the statistic -2RLL.
Results and discussion
The results of the combined analysis of variance of chickpea genotypes in two normal and dry conditions were shown in Table 2. These results showed that there was no significant difference between two normal and dry conditions in all traits except for pod weight and root dry weight. These results indicated that the two environments provided the same conditions for the genotypes on average. Also, the interaction of genotype in the environment was not statistically significant in all traits. After integrating the normal and dry data, the analysis of variance was carried out using randomized block design and lattice approach on data. The results showed that the spatial model (4), the nearest neighbor model, is the best linear model for analyzing and comparing the mean of genotypes. This model had the lowest -2RLL (-2 Res Log Likelihood). The lower the value of this statistic, the better the adequacy of the model to justify the response variable. Analysis of variance by the nearest neighbor method for different traits was shown in Table 4. The significant differences existed among genotypes in all traits except for pod weight and root dry weight at the 0.01probability level. Estimation of heritability of different traits showed that most of the studied traits had a heritability ranged from 0.39 to 0.75. For the grain yield traits with 10% intensity selection, the predicted genetic gain was 14% relative to the trait average. Therefore, this genetic gain can be used to advance the development of this trait. Also, genotypes G59, G09, G12, G30, and G41, respectively, in the nearest neighbor model, had the highest ranked estimated grain yields.
 Conclusions
The results of the analysis of the nearest neighbor model showed that there was a significant difference between all genotypes for all traits, as well as high genetic variance and acceptable heritability. Therefore, the variation in the genotypes studied could be used and introduced superior genotypes for breeding programs and even introductions to farmers.

Keywords

Arora, P.P., Jeena, A.S., 2001. Genetic variability studies in chickpea. Legume Research. 25, 137-138.
Arshad, M., Bakhsh, A., Ghafoor, A., 2004. Path coefficient analysis in chickpea (Cicer arietinum L.) under rainfed conditions. Pakistan Journal of Botany. 36, 75-82.
Bowman, D., 1990. Trend analysis to improve efficiency of agronomic trials in flue-cured tobacco. Agronomy Journal. 82, 499-501.
Burgueño, J., Cadena, A., Crossa, J., Banziger, M., Gilmour, A., Cullis, B. 2000. User's guide for spatial analysis of field variety trials using ASREML. Cimmyt, Mexico.
Chauhan, R.S., 2011. Studies on genetic variability for yield and quality traits of chickpea (Cicer arietinum) grown under late sown condition. Phd Thesis. Jnkvv, Jabalpur. India.
Coombes, N., 2009. DiGGer design search tool in R. New South Wales Department of Primary Industry) Available at http://nswdpibiom. org/austatgen/software/[Verified 29 August 2017].
Cullis, B.R., Smith, A.B., Coombes, N.E., 2006. On the design of early generation variety trials with correlated data. Journal of Agricultural, Biological, and Environmental Statistics. 11, 381-393.
Dutkowski, G.W., Costa e Silva, J., Gilmour, A. R., Wellendorf, H., Aguiar, A., 2006. Spatial analysis enhances modelling of a wide variety of traits in forest genetic trials. Canadian Journal of Forest Research. 36, 1851-1870.
FAO, 2017. FAOSTAT. Food and Agriculture Organization of the United Nations. http://faostat.fao.org/default.aspx
Federer, W.T., 1996. Recovery of interblock, intergradient, and intervariety information for incomplete block and lattice rectangle designed experiments. BU-1315-M in the Technical Report Series of the Biometrics Unit, Cornell University, Ithaca, NY 14853.
Federer, W.T., 1995. SAS for ANOVAs and recovery of interblock, intercolumn, and intergradient information. BU-1295 in the Technical Report Series of the Biometrics Unit, Cornell University, June.
Federer, W.T., Schlottfeldt, C.S., 1954. The use of covariance to control gradients in experiments. Biometrics. 10, 282-290.
Federer, W.T., Newton, E.A., Altman, N.S., 1997. Combining standard block analyses with spatial analyses under a random effects model. Pages 373-386 Modelling Longitudinal and Spatially Correlated Data. Springer, Netherland.
Holland, J., 2006. Estimating genotypic correlations and their standard errors using multivariate restricted maximum likelihood estimation with SAS Proc MIXED. Crop Science. 46, 642-654.
Iran Meteorological Organization, 2018. Iran Meteorological Organization, Khorramabad Station. Temperature and precipitation data for year 2018. [In Persian]
Ismaili, A., Karami, F., Akbarpour, O., Rezaeinejad, A., 2016. Estimation of genotypic correlation and heritability of apricot traits, using restricted maximum likelihood in repeated measures data. Canadian Journal of Plant Science. 96, 439-447.
Kirk, H.J., Haynes, F.L., Monroe, R.J., 1980. Application of trend analysis to horticultural field trials. Journal of the American Society of Horticultural Science. 105, 189 – 193.
Kumar, J., Abbo. S., 2001. Genetics of flowering time in chickpea and its bearing on productivity in semiarid environments. Advances in Agronomy. 72, 107-138.
Kumar, S., Arora, P., Jeena, A., 2001. Genetic variability studied for quantitative traits in chickpea. Agricultural Science Digest. 21, 263-264.
Littel, R., Milliken, G., Stroup, W., Wolfinger, R., 2006. SAS system for mixed models. SAS Institute Inc, Cary, NC, USA.
Malhotra, R.S., Singh, M., Erskine, W., 2004. Application of spatial variability models in enhancing precision and efficiency of selection in Chickpea trials. Journal of the Indian Society of Agricultural Statistics. 57, 71–83.
Ministry of Agriculture, 2017. Agriculture - Statistical Year Book of Iran. First volume: Crops. Information and Communication Technology Center, Department of Planning and Economy, Ministry of Agriculture, Tehran. Iran. [In Persian]
Nayyar, H., Chander, K., Kumar, S., Bains, T., 2005. Glycine betaine mitigates cold stress damage in chickpea. Agronomy for Sustainable Development. 25, 381-388.
Nyquist, W.E., Baker, R., 1991. Estimation of heritability and prediction of selection response in plant populations. Critical Reviews in Plant Sciences. 10, 235-322.
Papadakis, J., 1937. Methode statistique pour des experiences sur champ, Institut d'Amelioration des Plantes a Thessaloniki (Thessalonika, Greece). Bulletin Scientifique. 23, 1–30.
Parsa M, Bagheri A., 2008. Pulses. Mashhad University Press, Mashhad, Iran. [In Persian].
Piepho, H.P., Richter, C., Williams, E., 2008. Nearest neighbour adjustment and linear variance models in plant breeding trials. Biometrical Journal. 50,164-189.
Purushothaman, R., Upadhyaya, H., Gaur, P., Gowda, C., Krishnamurthy, L., 2014. Kabuli and desi chickpeas differ in their requirement for reproductive duration. Field Crops Research. 163, 24-31.
R Development Core Team, 2011. R: A Language and Environment for Statistical Computing. R Development Core Team, Vienna, Austria, http://www.r- project.org/.
Rezvani Moghaddam, P., Sadeghi Samarjan, R., 2008. Effect of sowing dates and different irrigation regimes on morphological characteristics and grain yield of chickpea (Cicer arietinum L.). Journal of Iranian Field Crop Research. 6, 315-325. [In Persian with English summary].
Roff, D. A., 2012. Evolutionary quantitative genetics. Springer Science & Business Media, Netherland.
Sandhu, T., Gumber, R., Bhullar, B., 1991. Correlated response of grain yield and protein content in chickpea. Legume Research, 14, 45-49.
Sandhu, T.S., Gumber, R.K., Bhullar, B.S., 1991. Correlated response of grain yield and protnin content in chickpea (Cicer arietinum L.). Legume Research. 14, 45-46.
Sarker, A., Singh, M., 2015. Improving breeding efficiency through application of appropriate experimental designs and analysis models: a case of lentil (Lens culinaris Medikus subsp. culinaris) yield trials. Field Crops Research. 179, 26-34.
Serraj, R., Krishnamurthy, L. Kashiwagi, J., Kumar, J., Chandra, S., Crouch, J., 2004. Variation in root traits of chickpea (Cicer arietinum L.) grown under terminal drought. Field Crops Research. 88,115-127.
Singh, M., Chaubey, Y., Sarker, A., Sen, D., 2010. Modeling unstructured heterogeneity along with spatially correlated errors in field trials. Journal of the Indian Society of Agricultural Statistics. 64, 313-321.
Singh, M., Malhotra, R., Ceccarelli, S., Sarker, A., Grando, S., Erskine, W., 2003. Spatial variability models to improve dryland field trials. Experimental Agriculture. 39, 151-160.
Stringer, J., Smith, A.B., Cullis, B.R., 2011. Spatial analysis of agricultural field experiments. In: Hinkelmann, K., (eds.), Design and analysis of experiments: special designs and applications, Volume 3 (pp. 109-136). Hoboken, N.J: Wiley-Interscience. USA.
Tamura, R.N., Nelson, L.A., Naderman, G.C., 1988. An investigation of the validity and usefulness of trend analysis for field plot data. Agronomy Journal. 80,712-718.
Warren, J., Mendez, I., 1982. Methods for Estimating Background Variation in Field Experiments. Agronomy Journal. 74, 1004-1009.
Wilkinson, G., Eckert, S., Hancock, T., Mayo, O., 1983. Nearest neighbour (NN) analysis of field experiments. Journal of the Royal Statistical Society. Series B (Methodological). 45(2), 151-211.
Zhang, H., Pala, M., Oweis, T., Harris, H., 2000. Water use and water-use efficiency of chickpea and lentil in a Mediterranean environment. Australian Journal of Agricultural Research. 51, 295-304.