Skim resequencing finely maps the downy mildew resistance loci RPF2 and RPF3 in spinach cultivars whale and Lazio

Abstract

Commercial production of spinach (Spinacia oleracea L.) is centered in California and Arizona in the US, where downy mildew caused by Peronospora effusa is the most destructive disease. Nineteen typical races of P. effusa have been reported to infect spinach, with 16 identified after 1990. The regular appearance of new pathogen races breaks the resistance gene introgressed in spinach. We attempted to map and delineate the RPF2 locus at a finer resolution, identify linked single nucleotide polymorphism (SNP) markers, and report candidate downy mildew resistance (R) genes. Progeny populations segregating for RPF2 locus derived from resistant differential cultivar Lazio were infected using race 5 of P. effusa and were used to study for genetic transmission and mapping analysis in this study. Association analysis performed with low coverage whole genome resequencing-generated SNP markers mapped the RPF2 locus between 0.47 to 1.46 Mb of chromosome 3 with peak SNP (Chr3_1,221,009) showing a LOD value of 61.6 in the GLM model in TASSEL, which was within 1.08 Kb from Spo12821, a gene that encodes CC-NBS-LRR plant disease resistance protein. In addition, a combined analysis of progeny panels of Lazio and Whale segregating for RPF2 and RPF3 loci delineated the resistance section in chromosome 3 between 1.18–1.23 and 1.75–1.76 Mb. This study provides valuable information on the RPF2 resistance region in the spinach cultivar Lazio compared to RPF3 loci in the cultivar Whale. The RPF2 and RPF3 specific SNP markers, plus the resistant genes reported here, could add value to breeding efforts to develop downy mildew resistant cultivars in the future.

Author contributions

AS, GB conceived the study. AS managed funds and resources. GB planned and performed the experiment, maintained the downy mildew pathogen, performed inoculation, phenotyping, genotyping, and data analysis, and wrote and revised the manuscript. BM provided the progeny population. JC provided P. effusa inoculum. All authors made edits and approved the final version of the manuscript.

Data availability

The genetic datasets generated in this study are publicly available in FigShare: https://doi.org/10.6084/m9.figshare.20441916.v2.

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary data

Supplementary data is available at Horticulture Research online.

References

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

[2] Morelock TE, Correll JC . Spinach. In: Prohens, J., Nuez, F. (eds) Vegetables I. Springer: New York, NY, 2008, 189-218.

[3] Bhattarai G, Shi A, Department of Horticulture, University of Arkansas, Fayetteville, AR 72701, USA. Research advances and prospects of spinach breeding, genetics, and genomics. Vegetable Research. 2021; 1: 1-18.

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

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

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

[8] 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: 740437.

[9] Ribera A, Bai Y, Wolters AMA et al. A review on the genetic resources, domestication and breeding history of spinach (Spinacia oleracea L.). Euphytica. 2020; 216: 48.

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

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

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

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

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

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

[16] Bhattarai G, Olaoye D, Mou B et al. Mapping and selection of downy mildew resistance in spinach cv. Whale by low coverage whole genome sequencing. Front Plant Sci. 2022; 13: 1012923.

[17] Gao S, Lu T, She H et al. Fine mapping and identification of a candidate gene of downy mildew resistance, RPF2, in spinach (Spinacia oleracea L.). Int J Mol Sci. 2022; 23: 14872.

[18] Bhattarai G, Shi A, Mou B et al. Resequencing worldwide spinach germplasm for identification of field resistance QTLs to downy mildew and assessment of genomic selection methods. Horticulture Research. 2022; 9: uhac205.

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

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

[21] Bayer PE, Ruperao P, Mason AS et al. High-resolution skim genotyping by sequencing reveals the distribution of crossovers and gene conversions in Cicer arietinum and Brassica napus. Theor Appl Genet. 2015; 128: 1039-47.

[22] Gao ZY, Zhao SC, He WM et al. Dissecting yield-associated loci in super hybrid rice by resequencing recombinant inbred lines and improving parental genome sequences. Proc Natl Acad Sci U S A. 2013; 110: 14492-7.

[23] Hu Z, Deng G, Mou H et al. A re-sequencing-based ultra-dense genetic map reveals a gummy stem blight resistance-associated gene in Cucumis melo. DNA Res. 2018; 25: 1-10.

[24] Malmberg MM, Barbulescu DM, Drayton MC et al. Evaluation and recommendations for routine genotyping using skim whole genome re-sequencing in canola. Front Plant Sci. 2018; 9: 1809.

[25] Bhattarai G, Feng C, Dhillon B et al. Detached leaf inoculation assay for evaluating resistance to the spinach downy mildew pathogen. Eur J Plant Pathol. 2020; 158: 511-20.

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

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

[28] Jones JDG, Dangl JL . The plant immune system. Nature. 2006; 444: 323-9.

[29] Gururani MA, Venkatesh J, Upadhyaya CP et al. Plant disease resistance genes: current status and future directions. Physiol Mol Plant Pathol. 2012; 78: 51-65.

[30] Innes RW . Guarding the goods. New insights into the central alarm system of plants. Plant Physiol. 2004; 135: 695-701.

[31] Sinapidou E, Williams K, Nott L et al. Two TIR:NB:LRR genes are required to specify resistance to Peronospora parasitica isolate Cala2 in Arabidopsis . Plant J. 2004; 38: 898-909.

[32] Ashikawa I, Hayashi N, Yamane H et al. Two adjacent nucleotide-binding site-leucine-rich repeat class genes are required to confer Pikm-specific rice blast resistance. Genetics. 2008; 180: 2267-76.

[33] Buckley RM, Harris AC, Wang GD et al. Best practices for analyzing imputed genotypes from low-pass sequencing in dogs. Mamm Genome. 2022; 33: 213-29.

[34] Browning BL, Browning SR . Genotype imputation with millions of reference samples. Am J Hum Genet. 2016; 98: 116-26.

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

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

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

[38] Gogarten SM, Sofer T, Chen H et al. Genetic association testing using the GENESIS R/bioconductor package. Bioinformatics. 2019; 35: 5346-8.