IIIVmrMLM Provides New Insights into the Genetic Basis of the Agronomic Trait Variation in Chickpea

Chickpea is a staple crop for many nations worldwide. Modeling genotype-by-environment interactions and assessing the genotype’s ability to contribute adaptive alleles are crucial for chickpea breeding. In this study, we evaluated 12 agronomically important traits of 159 accessions from the N.I. Vavilov All Russian Institute for Plant Genetic Resources collection. These included 145 landraces and 13 cultivars grown in different climatic conditions in Kuban (45°18′ N and 40°52′ E) in both 2016 and 2022, as well as in Astrakhan (46°06′ N and 48°04′ E) in 2022. Using the IIIVmrMLM model in multi-environmental mode, we identified 161 quantitative trait nucleotides (QTNs) with stable genetic effects across different environments. Furthermore, we have observed 254 QTN-by-environment interactions with distinct environment-specific effects. Notably, five of these interactions manifested large effects, with R2 values exceeding 10%, while the highest R2 value for stable QTNs was 4.7%. Within the protein-coding genes and their 1 Kb flanking regions, we have discerned 22 QTNs and 45 QTN-by-environment interactions, most likely tagging the candidate causal genes. The landraces obtained from the N.I Vavilov All Russian Institute for Plant Genetic Resources collection exhibit numerous favorable alleles at quantitative trait nucleotide loci, showing stable effects in the Kuban and Astrakhan regions. Additionally, they possessed a significantly higher number of Kuban-specific favorable alleles of the QTN-by-environment interaction loci compared to the Astrakhan-specific ones. The environment-specific alleles found at the QTN-by-environment interaction loci have the potential to enhance chickpea adaptation to specific climatic conditions.