Resequencing worldwide spinach germplasm for identification of field resistance QTLs to downy mildew and assessment of genomic selection methods

Abstract

Downy mildew, commercially the most important disease of spinach, is caused by the obligate oomycete Peronospora effusa. In the past two decades, new pathogen races have repeatedly overcome the resistance used in newly released cultivars, urging the need for more durable resistance. Commercial spinach cultivars are bred with major R genes to impart resistance to downy mildew pathogens and are effective against some pathogen races/isolates. This work aimed to evaluate the worldwide USDA spinach germplasm collections and commercial cultivars for resistance to downy mildew pathogen in the field condition under natural inoculum pressure and conduct genome wide association analysis (GWAS) to identify resistance-associated genomic regions (alleles). Another objective was to evaluate the prediction accuracy (PA) using several genomic prediction (GP) methods to assess the potential implementation of genomic selection (GS) to improve spinach breeding for resistance to downy mildew pathogen. More than four hundred diverse spinach genotypes comprising USDA germplasm accessions and commercial cultivars were evaluated for resistance to downy mildew pathogen between 2017–2019 in Salinas Valley, California and Yuma, Arizona. GWAS was performed using single nucleotide polymorphism (SNP) markers identified via whole genome resequencing (WGR) in GAPIT and TASSEL programs; detected 14, 12, 5, and 10 significantly associated SNP markers with the resistance from four tested environments, respectively; and the QTL alleles were detected at the previously reported region of chromosome 3 in three of the four experiments. In parallel, PA was assessed using six GP models and seven unique marker datasets for field resistance to downy mildew pathogen across four tested environments. The results suggest the suitability of GS to improve field resistance to downy mildew pathogen. The QTL, SNP markers, and PA estimates provide new information in spinach breeding to select resistant plants and breeding lines through marker-assisted selection (MAS) and GS, eventually helping to accumulate beneficial alleles for durable disease resistance.

References

[1] USDA National Agricultural Statistics Service . NASS - quick stats. USDA National Agricultural Statistics Service. 2020. https://data.nal.usda.gov/dataset/nass-quick-stats (accessed 20 May 2022).

[2] Irish BM, Correll JC, Koike ST et al. Three new races of the spinach downy mildew pathogen identified by a modified set of spinach differentials. Plant Dis. 2007; 91: 1392-6.

[3] Feng C, Saito K, Liu B et al. New races and novel strains of the spinach downy mildew pathogen Peronospora effusa . Plant Dis. 2018; 102: 613-8.

[4] Feng C, Correll JC, Kammeijer KE et al. Identification of new races and deviating strains of the spinach downy mildew pathogen Peronospora farinosa f. sp. spinaciae . Plant Dis. 2014; 98: 145-52.

[5] Plantum. Denomination of Pe : 18 and 19, two new races of downy mildew in spinach. 2021. https://plantum.nl/denomination-of-pe-18-and-19-two-new-races-of-downy-mildew-in-spinach/ (accessed 25 May 2021).

[6] 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.

[7] Irish BM, Correll JC, Feng C et al. Characterization of a resistance locus (Pfs-1) to the spinach downy mildew pathogen (Peronospora farinosa f. sp. spinaciae) and development of a molecular marker linked to Pfs-1 . Phytopathology. 2008; 98: 894-900.

[8] Feng C, Bluhm B, Shi A et al. Development of molecular markers linked to three spinach downy mildew resistance loci. Euphytica. 2018; 214: 174.

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

[10] She H, Qian W, Zhang H et al. Fine mapping and candidate gene screening of the downy mildew resistance gene RPF1 in spinach. Theor Appl Genet. 2018; 131: 2529-41.

[11] 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: 563187.

[12] 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.

[13] Kolmer JA . Genetics of resistance to wheat leaf rust. Annu Rev Phytopathol. 1996; 34: 435-55.

[14] Peressotti E, Wiedemann-Merdinoglu S, Delmotte F et al. Breakdown of resistance to grapevine downy mildew upon limited deployment of a resistant variety. BMC Plant Biol. 2010; 10: 147.

[15] McDonald BA, Linde C . The population genetics of plant pathogens and breeding strategies for durable resistance. Euphytica. 2002; 124: 163-80.

[16] Parlevliet JE . Durability of resistance against fungal, bacterial and viral pathogens; present situation. Euphytica. 2002; 124: 147-56.

[17] Palloix A, Ayme V, Moury B . Durability of plant major resistance genes to pathogens depends on the genetic background, experimental evidence and consequences for breeding strategies. New Phytol. 2009; 183: 190-9.

[18] Quenouille J, Montarry J, Palloix A et al. Slower, stronger: how the plant genetic background protects a major resistance gene from breakdown. Mol Plant Pathol. 2013; 14: 109-18.

[19] Poland J, Rutkoski J . Advances and challenges in genomic selection for disease resistance. Annu Rev Phytopathol. 2016; 54: 79-98.

[20] Cai X, Sun X, Xu C et al. Genomic analyses provide insights into spinach domestication and the genetic basis of agronomic traits. Nat Commun. 2021; 12: 1-12.

[21] Simko I, Peng H, Sthapit Kandel J et al. Genome-wide association mapping reveals genomic regions frequently associated with lettuce field resistance to downy mildew. Theor Appl Genet. 2022; 135: 2009-24.

[22] Wang X, Bao K, Reddy UK et al. The USDA cucumber (Cucumis sativus L.) collection: genetic diversity, population structure, genome-wide association studies, and core collection development . Hortic Res. 2018; 5: 1-13.

[23] Colque-Little C, Abondano MC, Lund OS et al. Genetic variation for tolerance to the downy mildew pathogen Peronospora variabilis in genetic resources of quinoa (Chenopodium quinoa) . BMC Plant Biol. 2021; 21: 1-19.

[24] Talukder ZI, Ma G, Hulke BS et al. Linkage mapping and genome-wide association studies of the Rf gene cluster in sunflower (Helianthus annuus L.) and their distribution in world sunflower collections. Front Genet. 2019; 10: 216.

[25] Nemri A, Atwell S, Tarone AM et al. Genome-wide survey of arabidopsis natural variation in downy mildew resistance using combined association and linkage mapping. Proc Natl Acad Sci U S A. 2010; 107: 10302-7.

[26] Meuwissen THE, Hayes BJ, Goddard ME . Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001; 157: 1819-29.

[27] Jannink JL, Lorenz AJ, Iwata H . Genomic selection in plant breeding: from theory to practice. Brief Funct Genomics. 2010; 9: 166-77.

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

[29] Bernardo R. Genomewide selection with minimal crossing in self-pollinated crops. Crop Sci. 2010; 50: 624-7.

[30] Sehgal D, Rosyara U, Mondal S et al. Incorporating genome-wide association mapping results into genomic prediction models for grain yield and yield stability in CIMMYT spring bread wheat. Front Plant Sci. 2020; 11: 197.

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

[32] Lorenzana RE, Bernardo R . Accuracy of genotypic value predictions for marker-based selection in biparental plant populations. Theor Appl Genet. 2009; 120: 151-61.

[33] 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.

[34] Gezan SA, Osorio LF, Verma S et al. An experimental validation of genomic selection in octoploid strawberry. Hortic Res. 2017; 4: 16070.

[35] Islam MS, Fang DD, Jenkins JN et al. Evaluation of genomic selection methods for predicting fiber quality traits in upland cotton. Mol Gen Genomics. 2020; 295: 67-79.

[36] Shi A, Bhattarai G, Xiong H et al. Genome-wide association study and genomic prediction of white rust resistance in USDA GRIN spinach germplasm. Hortic Res. 2022; 9: uhac069.

[37] Habier D, Fernando RL, Garrick DJ . Genomic BLUP decoded: a look into the black box of genomic prediction. Genetics. 2013; 194: 597-607.

[38] Daetwyler HD, Calus MPL, Pong-Wong R et al. Genomic prediction in animals and plants: simulation of data, validation, reporting, and benchmarking. Genetics. 2013; 193: 347-65.

[39] Asoro FG, Newell MA, Beavis WD et al. Accuracy and training population design for genomic selection on quantitative traits in elite north American oats. The Plant Genome. 2011; 4: 132-44.

[40] Li H . A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. 2011; 27: 2987-93.

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

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

[43] Shi A, Qin J, Mou B et al. Genetic diversity and population structure analysis of spinach by single-nucleotide polymorphisms identified through genotyping-by-sequencing. PLoS One. 2017; 12: e0188745.

[44] Zhang C, Dong SS, Xu JY et al. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics. 2019; 35: 1786-8.

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

[46] 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. Giga Science. 2019; 8: 1-12.

[47] 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: e1005767.

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

[49] Liaw A, Wiener M . Classification and regression by randomForest. R News. 2002; 2: 18-22.

[50] Pérez P, De Los CG. Genome-wide regression and prediction with the BGLR statistical package. Genetics. 2014; 198: 483-95.

[51] Correll JC, Feng C, Matheron ME et al. Evaluation of spinach varieties for downy mildew resistance, Monterey County, CA, 2017. Plant Disease Management Reports. 2018; 12: V132.

[52] Dhillon B, Feng C, Bhattarai G et al. Evaluation of spinach varieties for downy mildew resistance, Yuma, AZ 2019. Plant Disease Management Reports. 2020; 14: V025.

[53] Dhillon B, Feng C, Bhattarai G et al. Evaluation of spinach varieties for downy mildew resistance, San Juan Bautista, CA 2018. Plant Disease Management Reports. 2019; 13: V017.

[54] Gyawali S, Bhattarai G, Shi A et al. Genetic diversity, structure, and selective sweeps in Spinacia turkestanica associated with the domestication of cultivated spinach . Front Genet. 2021; 12: 2469.

[55] Feng C, Villarroel-Zeballos M, Ficheus P et al. New races and novel strains of the spinach downy mildew pathogen. Phytopathology. 2020; 110: S2.1-S2.206.

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

[57] Mundt CC . Pyramiding for resistance durability: theory and practice. Phytopathology. 2018; 108: 792-802.

[58] Johnson R . Past present and future opportunities in breeding for disease resistance, with examples from wheat. Euphytica. 1992; 63: 3-22.

[59] Lindhout P. The perspectives of polygenic resistance in breeding for durable disease resistance. Euphytica. 2002; 124: 217-26.

[60] He L, Xiao J, Rashid KY et al. Genome-wide association studies for pasmo resistance in flax (Linum usitatissimum L.) . Front Plant Sci. 2018; 9: 1982.

[61] Liu X, Lu H, Liu P et al. Identification of novel loci and candidate genes for cucumber downy mildew resistance using GWAS. Plan Theory. 2020; 9: 1-15.

[62] Enciso-Rodriguez F, Douches D, Lopez-Cruz M et al. Genomic selection for late blight and common scab resistance in tetraploid potato (Solanum tuberosum) . G3: Genes, Genomes, Genetics. 2018; 8: 2471-81.

[63] Merrick LF, Burke AB, Chen X et al. Breeding with major and minor genes: genomic selection for quantitative disease resistance. Front Plant Sci. 2021; 12: 713667.

[64] Daetwyler HD, Pong-Wong R, Villanueva B et al. The impact of genetic architecture on genome-wide evaluation methods. Genetics. 2010; 185: 1021-31.

[65] Heslot N, Yang HP, Sorrells ME et al. Genomic selection in plant breeding: a comparison of models. Crop Sci. 2012; 52: 146-60.

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

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

[68] Rutkoski J, Singh RP, Huerta-Espino J et al. Efficient use of historical data for genomic selection: a case study of stem rust resistance in wheat. The Plant Genome. 2015; 8.

[69] Muleta KT, Bulli P, Zhang Z et al. Unlocking diversity in germplasm collections via genomic selection: a case study based on quantitative adult plant resistance to stripe rust in spring wheat. The Plant Genome. 2017; 10: 1-15.

[70] McElroy MS, Navarro AJR, Mustiga G et al. Prediction of cacao (Theobroma cacao) resistance to Moniliophthora spp. diseases via genome-wide association analysis and genomic selection. Front Plant Sci. 2018; 9: 343.