Genome-wide association study and genomic prediction of white rust resistance in USDA GRIN spinach germplasm

Abstract

White rust, caused by Albugo occidentalis, is one of the major yield-limiting diseases of spinach (Spinacia oleracea) in some major commercial production areas, particularly in southern Texas in the United States. The use of host resistance is the most economical and environment-friendly approach to managing white rust in spinach production. The objectives of this study were to conduct a genome-wide associating study (GWAS), to identify single nucleotide polymorphism (SNP) markers associated with white rust resistance in spinach, and to perform genomic prediction (GP) to estimate the prediction accuracy (PA). A GWAS panel of 346 USDA (US Dept. of Agriculture) germplasm accessions was phenotyped for white rust resistance under field conditions and GWAS was performed using 13 235 whole-genome resequencing (WGR) generated SNPs. Nine SNPs, chr2_53 049 132, chr3_58 479 501, chr3_95 114 909, chr4_9 176 069, chr4_17 807 168, chr4_83 938 338, chr4_87 601 768, chr6_1 877 096, and chr6_31 287 118, located on chromosomes 2, 3, 4, and 6 were associated with white rust resistance in this GWAS panel. Four scenarios were tested for PA using Pearson’s correlation coefficient (r) between the genomic estimation breeding value (GEBV) and the observed values: (1) different ratios between the training set and testing set (fold), (2) different GP models, (3) different SNP numbers in three different SNP sets, and (4) the use of GWAS-derived significant SNP markers. The results indicated that a 2- to 10-fold difference in the various GP models had similar, although not identical, averaged r values in each SNP set; using GWAS-derived significant SNP markers would increase PA with a high r-value up to 0.84. The SNP markers and the high PA can provide valuable information for breeders to improve spinach by marker-assisted selection (MAS) and genomic selection (GS).

References

[1] Davis W, Lucier G . Vegetable and pulses outlook No. (VGS-366) 68 pp: April 2021. United States Dep Agric Econ Res Serv. 2021.

[2] Dicoteau DR . Vegetable Crops . New Jersey: Printice Hall; 2000.

[3] Morelock TE, Correll JC . Spinach. In: J. Prohens and F. Nuez (eds.), Handbook of Plant Breeding, Vegetables I, Asteraceae, Brassicaceae, Chenopodicaceae, and Cucurbitaceae.New York: Springer, 2008, 189-218.

[4] Brandenberger LP, Correll JC, Morelock TE et al. Characterization of resistance of spinach to white rust (Albugo occidentalis) and downy mildew (Peronospora farinosa f.sp. spinaciae). Phytopathology. 1994;84:431-7.

[5] Correll JC, Bluhm BH, Feng C et al. Spinach: better management of downy mildew and white rust through genomics. Eur J Plant Pathol. 2011; 129: 193-205.

[6] Choi D, Priest MJ . A key to the genus Albugo. Mycotaxon. 1995; 53: 261-72.

[7] Correll JC, Black MC, Koike ST et al. Economically important diseases of spinach. Plant Dis. 1994; 78: 653-60.

[8] Goreta S, Leskovar DI . Screening spinach cultivars for white rust resistance and bolting. HortTechnology. 2006; 16: 162-6.

[9] Vakalounakis DJ, Doulis AG . First record of white rust, caused by Albugo occidentalis, on spinach in Greece. Plant Dis. 2013; 97: 1253.

[10] Correll JC, Feng CD, Liu B . First report of white rust (Albugo occidentalis) of spinach in Mexico . Plant Dis. 2017; 101: 511.

[11] Soylu S, Kara M, Kurt S et al. First report of white blister rust disease caused by Albugo occidentalis on spinach in Turkey. Plant Dis. 2018; 102: 826.

[12] Black MC, Dainello FJ, Kunkel TE . Fungicide evaluations for spinach white rust control on resistant and susceptible cultivars. Fungic Nematic Tests. 1992; 47: 1128.

[13] Dainello F, Black M, Disease TKP et al. Control of white rust of spinach with partial resistance and multiple soil applications of metalaxyl granules. Plant Dis. 1990; 74: 913-6.

[14] Dainello F, Heineman R . Relative white rust resistance and adaptability of spinach varieties in Southwest Texas. Texas Agricultural Experiment Station. 1981; 3878: 7 pages.

[15] Dohm JC, Minoche AE, Holtgrawe D et al. The genome of the recently domesticated crop plant sugar beet (Beta vulgaris) . Nature. 2014; 505: 546-9.

[16] Hirakawa H, Toyoda A, Itoh T et al. A spinach genome assembly with remarkable completeness, and its use for rapid identification of candidate genes for agronomic traits. DNA Res. 2021; 28.

[17] Hulse-Kemp AM, Bostan H, Chen S et al. An anchored chromosome-scale genome assembly of spinach improves annotation and reveals extensive gene rearrangements in euasterids. Plant Genome. 2021; 14: e20101.

[18] Xu C, Jiao C, Sun H et al. Draft genome of spinach and transcriptome diversity of 120 Spinacia accessions. Nat Commun. 2017; 8: 15275.

[19] Cai X, Sun X, Xu C et al. Reference genome and resequencing of 305 accessions provide insights into spinach evolution, domestication and genetic basis of agronomic traits. Nat Commun. 2021; 12: 7246.

[20] Davey JW, Hohenlohe PA, Etter PD et al. Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nat Rev Genet. 2011; 12: 499-510.

[21] Li H, Peng Z, Yang X et al. Genome-wide association study dissects the genetic architecture of oil biosynthesis in maize kernels. Nat Genet. 2013; 45: 43-50.

[22] Yano K, Yamamoto E, Aya K et al. Genome-wide association study using whole-genome sequencing rapidly identifies new genes influencing agronomic traits in rice. Nat Genet. 2016; 48: 927-34.

[23] Ma J, Shi A, Mou B et al. Association mapping of leaf traits in spinach (Spinacia oleracea L.) . Plant Breed. 2016; 135: 399-404.

[24] Chitwood J, Shi A, Mou B et al. Population structure and association analysis of bolting, plant height, and leaf erectness in spinach. HortScience. 2016; 51: 481-6.

[25] Shi A, Mou B, Cheng ZM . Genetic diversity and association analysis of leafminer (Liriomyza langei) resistance in spinach (Spinacia oleracea) . Genome. 2016; 59: 581-8.

[26] Shi A, Mou B, Correll JC . Association analysis for oxalate concentration in spinach. Euphytica. 2016; 212: 17-28.

[27] Shi A, Mou B, Correll J et al. SNP association analysis of resistance to Verticillium wilt (Verticillium dahliae Kleb.) in spinach . Aust J Crop Sci. 2016; 10: 1188-96.

[28] Shi A, Mou B, Correll J et al. Association analysis and identification of SNP markers for Stemphylium leaf spot (Stemphylium botryosum f. sp. spinacia) resistance in spinach (Spinacia oleracea) . Am J Plant Sci. 2016; 07: 1600-11.

[29] Qin J, Shi A, Mou B et al. Genetic diversity and association mapping of mineral element concentrations in spinach leaves. BMC Genomics. 2017; 18: 941.

[30] Awika HO, Marconi TG, Bedre R et al. Minor alleles are associated with white rust (Albugo occidentalis) susceptibility in spinach (Spinacia oleracea). Hortic Res. 2019;6:129.

[31] Awika HO, Bedre R, Yeom J et al. Developing growth-associated molecular markers via high-throughput phenotyping in spinach. Plant Genome. 2019; 139: 402-18.

[32] Awika HO, Cochran K, Joshi V et al. Single-marker and haplotype-based association analysis of anthracnose (Colletotrichum dematium) resistance in spinach (Spinacia oleracea) . Plant Breed. 2020; 139: 402-18.

[33] Bhattarai G. Genetic resistance to the downy mildew pathogen and mapping the RPF resistance loci in spinach. In: Doctoral Dissertation . University of Arkansas: Fayetteville, USA, 2019. Available at https://scholarworks.uark.edu/etd/3442.

[34] Bhattarai G, Shi A, Feng C et al. Genome wide association studies in multiple spinach breeding populations refine downy mildew race 13 resistance genes. Front Plant Sci. 2020; 11.

[35] Bhattarai G, Yang W, Shi A et al. High resolution mapping and candidate gene identification of downy mildew race 16 resistance in spinach. BMC Genomics. 2021; 22: 478.

[36] Olaoye D. Resistance Screening and Association Mapping for Resistance to the Downy Mildew Pathogen of Spinach. In: MS Thesis . University of Arkansas: Fayetteville, USA, 2021. Available at https://scholarworks.uark.edu/etd/4107.

[37] Brandenberger LP, Correll JC, Morelock TE . Identification of and cultivar reactions to a new race (race 4) of Peronospora farinosa f. sp. spinaciae on spinach in the United States . Plant Dis. 1991; 75: 630-4.

[38] Albrecht T et al. Genome-based prediction of testcross values in maize. Theor Appl Genet. 2011; 123: 339-50.

[39] Bernardo R, Yu J . Prospects for genome-wide selection for quantitative traits in maize. Crop Sci. 2007; 47: 1082-90.

[40] Zhang J, Song Q, Cregan P, Jiang GL . Genome-wide association study, genomic prediction and marker-assisted selection for seed weight in soybean (Glycine max). Theor Appl Genet. 2016; 129(1): 117-30.

[41] Heffner EL, Jannink JL, Iwata H et al. Genomic selection accuracy for grain quality traits in biparental wheat populations. Crop Sci. 2011; 51: 2597-606.

[42] Heffner EL, Sorrells ME, Jannink JL . Genomic selection for crop improvement. Crop Sci. 2009; 49: 1-12.

[43] Onogi A et al. Toward integration of genomic selection with crop modelling: the development of an integrated approach to predicting rice heading dates. Theor Appl Genet. 2016; 129: 805-17.

[44] Poland J, Endelman J, Dawson J et al. Genomic selection in wheat breeding using genotyping-by-sequencing. Plant Genome. 2012; 5.

[45] Poudel HP, Sanciangco MD, Kaeppler SM et al. Genomic prediction for winter survival of lowland switchgrass in the northern USA. G3 Genes Genomes Genet. 2019; 9: 1921-31.

[46] Würschum T, Reif JC, Kraft T et al. Genomic selection in sugar beet breeding populations. BMC Genet. 2013; 14: 85-8.

[47] Hernandez CO, Wyatt LE, Mazourek MR . Genomic prediction and selection for fruit traits in winter squash. G3 Genes|Genomes|Genetics. 2020; 10: 3601- 10.

[48] Xavier A, Muir WM, Rainey KM et al. Assessing predictive properties of genome-wide selection in soybeans. G3 (Bethesda). 2016; 6: 2611-6.

[49] Ravelombola W, Shi A, Huynh BL . Loci discovery, network-guided approach, and genomic prediction for drought tolerance index in a multi-parent advanced generation intercross (MAGIC) cowpea population. Hortic Res. 2021; 8: 24.

[50] Ravelombola WS, Qin J, Shi A et al. Genome-wide association study and genomic selection for tolerance of soybean biomass to soybean cyst nematode infestation. PLoS One. 2020; 15: e0235089.

[51] Ravelombola WS, Qin J, Shi A et al. Genome-wide association study and genomic selection for soybean chlorophyll content associated with soybean cyst nematode tolerance. BMC Genomics. 2019; 20.

[52] Bao Y, Vuong T, Meinhardt C et al. Potential of association mapping and genomic selection to explore PI 88788 derived soybean cyst nematode resistance. Crop Science. 2014; 7.

[53] Shi A, Gepts P, Song Q et al. Genome-wide association study and genomic prediction for soybean cyst nematode resistance in USDA common bean (Phaseolus vulgaris) Core collection. Front Plant Sci. 2021; 12: 624156.

[54] Wang J, Zhang Z . GAPIT version 3: boosting power and accuracy for genomic association and prediction. Genomics Proteomics Bioinformatics. 2021; 19: 629-40.

[55] Zeng J, Pszczola M, Wolc A et al. Genomic breeding value prediction and QTL mapping of QTLMAS2011 data using Bayesian and GBLUP methods. BMC Proc. 2012; 6: S7.

[56] Van Dijk EL, Jaszczyszyn Y, Thermes C . Library preparation methods for next-generation sequencing: tone down the bias. Exp Cell Res. 2014; 322: 12-20.

[57] Collins K, Zhao K, Jiao C et al. SpinachBase: a central portal for spinach genomics. Database (Oxford). 2019; 2019: baz072.

[58] Li H, Durbin R . Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009; 25: 1754-60.

[59] McKenna A, Hanna M, Banks E et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20: 1297-303.

[60] Lipka AE, Tian F, Wang Q et al. GAPIT: genome association and prediction integrated tool. Bioinformatics. 2012; 28: 2397-9.

[61] Bradbury PJ, Zhang Z, Kroon DE et al. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007; 23: 2633-5.

[62] Zhang Z, Ersoz E, Lai CQ et al. Mixed linear model approach adapted for genome-wide association studies. Nat Genet. 2010; 42: 355-60.

[63] Wang Q, Tian F, Pan Y et al. A SUPER powerful method for genome wide association study. PLoS One. 2014; 9: e107684.

[64] Liu X, Huang M, Fan B et al. Iterative usage of fixed and random effect models for powerful and efficient genome-wide association studies. PLoS Genet. 2016; 12.

[65] Huang M, Liu X, Zhou Y et al. BLINK: a package for the next level of genome-wide association studies with both individuals and markers in the millions. Gigascience. 2019; 8.

[66] Endelman JB . Ridge regression and other kernels for genomic selection with R package rrBLUP. Plant Genome. 2011; 4: 250-5.

[67] Legarra A, Robert-Granié C, Croiseau P et al. Improved lasso for genomic selection. Genet Res (Camb). 2011; 93: 77-87.

[68] Ogutu JO, Piepho HP, Schulz-Streeck T . A comparison of random forests, boosting and support vector machines for genomic selection. BMC Proc. 2011; 5: 1-5.

[69] Shikha M, Kanika A, Roa AR et al. Genomic selection for drought tolerance using genome-wide SNPs in maize. Front Plant Sci. 2017; 8.

[70] Brandenberger L . Studies to quantify disease resistance in spinach to the white rust (Albugo occidentalis) and downy mildew (Peronospora farinosa f. sp. spinaciae) pathogens. In: Doctoral Dissertation.University of Arkansas: Fayetteville, USA, 1992.

[71] Keller B, Ariza-Suarez D, de la Hoz J et al. Genomic prediction of agronomic traits in common bean (Phaseolus vulgaris L.) under environmental stress . Front Plant Sci. 2020; 11.

[72] Qin J, Shi A, Song Q et al. Genome wide association study and genomic selection of amino acid concentrations in soybean seeds. Front Plant Sci. 2019; 10.

[73] Spindel JE, Begum H, Akdemir D et al. Genome-wide prediction models that incorporate de novo GWAS are a powerful new tool for tropical rice improvement. Heredity (Edinb). 2016; 116: 395-408.

[74] Bhattarai G, Shi A . Research advances and prospects of spinach breeding, genetics, and genomics. Vegetable Research. 2021; 1: 1-18.