Genome-wide association study on stem rust resistance in Kazakh spring barley lines

Yerlan Turuspekov, Danara Ormanbekova, Aralbek Rsaliev, Saule Abugalieva

Stem rust (SR) is one of the most economically devastating barley diseases in Kazakhstan, and in some years it causes up to 50 % of yield losses. Routine conventional breeding for resistance to stem rust is almost always in progress in all Kazakhstan breeding stations. However, molecular marker based approach towards new SR resistance genes identification and relevant marker-assisted selection had never been employed in Kazakhstan yet. In this study, as a preliminary step the GWAS (genome-wide association study) mapping was applied in attempt to identify reliable SNP markers and quantitative trait loci (QTL) associated with resistance to SR. Barley collection of 92 commercial cultivars and promising lines was genotyped using a high-throughput single nucleotide polymorphism (9,000 SNP) Illumina iSelect array. 6,970 SNPs out of 9,000 total were polymorphic and scorable. 5,050 SNPs out of 6,970 passed filtering threshold and were used for association mapping (AM). All 92 accessions were phenotyped for resistance to SR by observing adult plants in artificially infected plots at the Research Institute for Biological Safety Problems in Dzhambul region of Kazakhstan. GLM analysis allowed the identification of 15 SNPs associated with the resistance at the heading time (HA) phase, and 2 SNPs at the seed's milky-waxy maturity (SM) phase. However, after application of 5 % Bonferroni multiple test correction, only 2 SNPs at the HA stage on the same position of chromosome 6H can be claimed as reliable markers for SR resistance. The MLM analysis after the Bonferroni correction did not reveal any associations in this study, although distribution lines in the quantile-quantile (QQ) plot indicates that overcorrection in the test due to both Q and K matrices usage. Obtained data suggest that genome wide genotyping of 92 spring barley accessions from Kazakhstan with 9 K Illumina SNP array was highly efficient. Linkage disequilibrium based mapping approach allowed the identification of highly significant QTL for the SR resistance at the HA phase of growth on chromosome 6H. On the other hand, no significant QTL was detected at the SM phase, assuming that for a successful GWASmapping experiment a larger size population with more diverse genetic background should be tested. Obtained results provide additional information towards better understanding of SR resistance in barley.