Chromosome-level reference genome of Vitis piasezkii var. pagnucii provides insights into a new locus of resistance to grapevine powdery mildew

Abstract

Grapevine powdery mildew (GPM), caused by Erysiphe necator, poses a significant threat to all green grapevine tissues, leading to substantial economic losses in viticulture. Traditional grapevine cultivars derived from Vitis vinifera are highly susceptible to GPM, whereas the wild Chinese accession Baishui-40 (BS-40) of V. piasezkii var. pagnucii exhibits robust resistance. To illuminate the genetic basis of resistance, we sequenced and assembled the chromosome-level genome of ‘BS-40’, achieving a total mapped length of 578.6 Mb distributed across nineteen chromosomes. A comprehensive annotation identified 897 nucleotide-binding leucine-rich repeat (NLR) genes in the ‘BS-40’ genome, which exhibited high sequence similarity across Vitis genomes. 284 of these NLR genes were differentially expressed upon GPM infection. A hybrid population of ‘BS-40’ and V. vinifera was constructed and 195 progenies were whole-genome re-sequenced. A new GPM-resistant locus, designated Ren17, located within the 0.74-1.23 Mb region on chromosome 1 was identified using genome-wide association study, population selection, and QTL analysis. Recombinant events indicated that an NLR gene cluster between 1 045 489 and 1 089 719 bp on chromosome 1 is possibly the key contributor to GPM resistance in ‘BS-40’. Based on an SNP within this region, a dCAPS marker was developed that can predict the GPM resistance in ‘BS-40’-derived materials with 99.4% accuracy in the progenies of ‘BS-40’ and V. vinifera. This chromosome-level genome assembly of V. piasezkii var. pagnucii provides a valuable resource not only for grapevine evolution, genetic analysis, and pan-genome studies but also a new locus Ren17 as a promising target for GPM-resistant breeding in grapevine.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant No. 32272670) and the National Key Research and Development Program of China (2018YFD1000300). We thank High-Performance Computing (HPC) of Northwest A&F University for providing computing resources. We thank the Horticulture Science Research Center at College of Horticulture, Northwest A&F University for experimental assistance.

Author Contributions

Y.Q.W. and Y.J. conceived the study. L.Z. and Y.H. performed the genome sequencing and analysis. L.Z., Y.H., L.X.G., M.J., X.N.Y., Z.S.J., and M.Z. performed the construction of crossing population, phenotype investigation, and sample collection for resequencing. L.Z., Q.Y.J., and Z.S.J. conducted the dCAPS marker experiment. L.Z. and X.L.D. analyzed transcriptome and resequencing data. L.Z. wrote the manuscript. Y.Q.W., Y.J., and S.Y.X. revised the manuscript. All authors have read and approved the final manuscript.

Data Availability

The genome sequence of ‘BS-40’ has been submitted to the National Center for Biotechnology Information (BioProject ID: PRJNA1104008, accession number: JBDIVC000000000). The raw sequences of the transcriptome of ‘BS-40’ and ‘CS’ were submitted to the SRA database (BioProject ID: PRJNA1105926). The resequencing sequences of ‘PC’ and its’ progenies crossed with ‘BS-40′ were submitted to the SRA database (BioProject ID: PRJNA1105918).

Conflicts of Interest

The authors declare no conflicts of interest.

Supplementary data

Supplementary data is available at Horticulture Research online.

References

[1] Dong Y, Duan S, Xia Q. et al. Dual domestications and origin of traits in grapevine evolution. Science. 2023; 379:892-901

[2] Qiu W, Feechan A, Dry I. Current understanding of grapevine defense mechanisms against the biotrophic fungus (Erysiphe necator), the causal agent of powdery mildew disease. Hortic Res. 2015; 2:15020

[3] Hu Y, Gao YR, Yang LS. et al. The cytological basis of powdery mildew resistance in wild Chinese Vitis species. Plant Physiol Biochem. 2019; 144:244-53

[4] Pap D, Riaz S, Dry IB. et al. Identification of two novel powdery mildew resistance loci, Ren6 and Ren7, from the wild Chinese grape species Vitis piasezkii. BMC Plant Biol. 2016; 16:170

[5] Calonnec A, Cartolaro P, Poupot C. et al. Effects of Uncinula necator on the yield and quality of grapes (Vitis vinifera) and wine. Plant Pathol. 2004; 53:434-45

[6] Gadoury DM, Seem RC, Wilcox WF. et al. Effects of diffuse colonization of grape berries by Uncinula necator on bunch rots, berry microflora, and juice and wine quality. Phytopathology. 2007; 97:1356-65

[7] Lowder SR, Neill TM, Peetz AB. et al. A rapid glove-based inocu-lum sampling technique to monitor Erysiphe necator in com-mercial vineyards. Plant Dis. 2023; 107:3096-105

[8] Sapkota S, Zou C, Ledbetter C. et al. Discovery and genome-guided mapping of REN12 from Vitis amurensis, conferring strong, rapid resistance to grapevine powdery mildew. Hortic Res. 2023;10:uhad052

[9] Dry IB, Feechan A, Anderson C. et al. Molecular strategies to enhance the genetic resistance of grapevines to powdery mildew. Aust J Grape Wine Res. 2010; 16:94-105

[10] Barker CL, Donald T, Pauquet J. et al. Genetic and physical map-ping of the grapevine powdery mildew resistance gene, Run1, using a bacterial artificial chromosome library. Theor Appl Genet. 2005; 111:370-7

[11] Riaz S, Tenscher AC, Ramming DW. et al. Using a limited mapping strategy to identify major QTLs for resistance to grapevine pow-dery mildew (Erysiphe necator) and their use in marker-assisted breeding. Theor Appl Genet. 2011; 122:1059-73

[12] Blanc S, Wiedemann-Merdinoglu S, Dumas V. et al. A reference genetic map of Muscadinia rotundifolia and identification of Ren5, a new major locus for resistance to grapevine powdery mildew. Theor Appl Genet. 2012; 125:1663-75

[13] Karn A, Zou C, Brooks S. et al. Discovery of the REN11 locus from Vitis aestivalis for stable resistance to grapevine powdery mildew in a family segregating for several unstable and tissue-specific quantitative resistance loci. Front Plant Sci. 2021; 12: 733899

[14] Dalbó MA, Ye GN, Weeden NF. et al. Marker-assisted selection for powdery mildew resistance in grapes. J Am Soc Hortic Sci. 2001; 126:83-9

[15] Teh SL, Fresnedo-Ramírez J, Clark MD. et al. Genetic dissection of powdery mildew resistance in interspecific half-sib grapevine families using SNP-based maps. Mol Breed. 2016; 37:1

[16] Zyprian E, Ochßner I, Schwander F. et al. Quantitative trait loci affecting pathogen resistance and ripening of grapevines. Mol Gen Genomics. 2016; 291:1573-94

[17] Zendler D, Schneider P, Töpfer R. et al. Fine mapping of Ren3 reveals two loci mediating hypersensitive response against Erysiphe necator in grapevine. Euphytica. 2017; 213:68

[18] Ramming DW, Gabler F, Smilanick J. et al. A single domi-nant locus, ren4, confers rapid non-race-specific resistance to grapevine powdery mildew. Phytopathology. 2011; 101:502-8

[19] Hoffmann S, di Gaspero G, Kovács L. et al. Resistance to Erysiphe necator in the grapevine ’Kishmish vatkana’ is controlled by a single locus through restriction of hyphal growth. Theor Appl Genet. 2008; 116:427-38

[20] Feechan A, Anderson C, Torregrosa L. et al. Genetic dissection of a TIR-NB-LRR locus from the wild north American grapevine species Muscadinia rotundifolia identifies paralogous genes con-ferring resistance to major fungal and oomycete pathogens in cultivated grapevine. Plant Journal. 2013; 76:661-74

[21] Li P, Tan X, Liu R. et al. QTL detection and candidate gene analysis of grape white rot resistance by interspecific grape (Vitis vinifera L. × Vitis davidii Foex.) crossing. Hortic Res. 2023;10:uhad063

[22] Fu P, Wu W, Lai G. et al. Identifying Plasmopara viticola resistance loci in grapevine (Vitis amurensis) via genotyping-by-sequencing-based QTL mapping. Plant Physiol Biochem. 2020; 154:75-84

[23] Possamai T, Wiedemann-Merdinoglu S, Merdinoglu D. et al. Con-struction of a high-density genetic map and detection of a major QTL of resistance to powdery mildew (Erysiphe necator Sch.) in Caucasian grapes (Vitis vinifera L.). BMC Plant Biol. 2021; 21:528

[24] Guo DL, Zhao HL, Li Q. et al. Genome-wide association study of berry-related traits in grape [Vitis vinifera L.] based on genotyping-by-sequencing markers. Hortic Res. 2019; 6:11

[25] Sargolzaei M, Maddalena G, Bitsadze N. et al. Rpv29, Rpv30 and Rpv31: three novel genomic loci associated with resistance to Plasmopara viticola in Vitis vinifera. Front Plant Sci. 2020; 11:562432

[26] Jang HA, Oo MM, Kim DG. et al. CC-NBS-LRR, a set of VvCRP markers, can distinguish cultivars with ripe rot resistance to Colletotrichum pathogens in grapevine. Horticulture, Environment, and Biotechnology. 2020; 61:915-27

[27] Lian X, Zhang H, Jiang C. et al. De novo chromosome-level genome of a semi-dwarf cultivar of Prunus persica identifies the aquaporin PpTIP2 as responsible for temperature-sensitive semi-dwarf trait and PpB3-1 for flower type and size. Plant Biotechnol J. 2022; 20:886-902

[28] Sun Y, Yang B, Li M. et al. Rpv34: the noval Plasmopara viticola resistance locus in ’Moldova’. Fruit Res. 2024; 4:0

[29] Jaillon O, Aury JM, Noel B. et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007; 449:463-7

[30] Canaguier A, Grimplet J, di Gaspero G. et al. A new version of the grapevine reference genome assembly (12X.v2) and of its annotation (VCost.v3). Genom Data. 2017; 14:56-62

[31] Shi X, Cao S, Wang X. et al. The complete reference genome for grapevine (Vitis vinifera L.) genetics and breeding. Hortic Res. 2023;10:uhad061

[32] Velt A, Frommer B, Blanc S. et al. An improved reference of the grapevine genome reasserts the origin of the PN40024 highly homozygous genotype. G3 (Bethesda). 2023;13:jkad067

[33] Chin C-S, Peluso P, Sedlazeck FJ. et al. Phased diploid genome assembly with single-molecule real-time sequencing. Nat Meth-ods. 2016; 13:1050-4

[34] Zhou Y, Minio A, Massonnet M. et al. The population genetics of structural variants in grapevine domestication. Nat Plants. 2019; 5:965-79

[35] Massonnet M, Cochetel N, Minio A. et al. The genetic basis of sex determination in grapes. Nat Commun. 2020; 11:2902

[36] Howard S, Bouck J, Qiu WP. Assembly and analysis of the chloro-plast genome of the north American hybrid grape Ńorton´. Vitis. 2016; 55:95-8

[37] Patel S, Robben M, Fennell A. et al. Draft genome of the native American cold hardy grapevine Vitis riparia Michx. ´Manitoba 37´. Hortic Res. 2020; 7:92

[38] Minio A, Cochetel N, Massonnet M. et al. HiFi chromosome-scale diploid assemblies of the grape rootstocks 110R, Kober 5BB, and 101-14 Mgt. Sci Data. 2022; 9:660

[39] Park M, Vera D, Kambrianda D. et al. Chromosome-level genome sequence assembly and genome-wide association study of Muscadinia rotundifolia reveal the genetics of 12 berry-related traits. Hortic Res. 2022;9:uhab011

[40] Wang Y, Xin H, Fan P. et al. The genome of Shanputao (Vitis amurensis) provides a new insight into cold tolerance of grapevine. Plant Journal. 2021; 105:1495-506

[41] Wang P, Meng F, Yang Y. et al. De novo assembling a high-quality genome sequence of Amur grape (Vitis amurensis Rupr.) gives insight into Vitis divergence and sex determination. Hortic Res. 2024;11:uhae117

[42] Wang Y, Ding K, Li H. et al. Biography of Vitis genomics: recent advances and prospective. Hortic Res. 2024;11:uhae128

[43] Liu Z, Wang N, Su Y. et al. Grapevine pangenome facilitates trait genetics and genomic breeding. Nat Genet. 2024; 56:2804-14

[44] Cochetel N, Minio A, Guarracino A. et al. A super-pangenome of the north American wild grape species. Genome Biol. 2023; 24:290

[45] Gao YR, Han YT, Zhao FL. et al. Identification and utilization ofanew Erysiphe necator isolate NAFU1 to quickly evaluate powdery mildew resistance in wild Chinese grapevine species using detached leaves. Plant Physiol Biochem. 2016; 98:12-24

[46] Jiao F, Luo R, Dai X. et al. Chromosome-level reference genome and population genomic analysis provide insights into the evo-lution and improvement of domesticated mulberry (Morus alba). Mol Plant. 2020; 13:1001-12

[47] Peng Y, Wang Y, Liu Y. et al. The genomic and epigenomic landscapes of hemizygous genes across crops with contrast-ing reproductive systems. Proc Natl Acad Sci USA. 2025; 122: e2422487122

[48] Ngou BPM, Ahn HK, Ding P. et al. Mutual potentiation of plant immunity by cell-surface and intracellular receptors. Nature. 2021; 592:110-5

[49] Yuan M, Jiang Z, Bi G. et al. Pattern-recognition receptors are required for NLR-mediated plant immunity. Nature. 2021; 592: 105-9

[50] Wasternack C, Feussner I. The Oxylipin pathways: biochemistry and function. Annu Rev Plant Biol. 2018; 69:363-86

[51] Dorion S, Ouellet JC, Rivoal J. Glutathione metabolism in plants under stress: beyond reactive oxygen species detoxification. Metabolites. 2021; 11:641

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

[53] Long Q, Cao S, Huang G. et al. Population comparative genomics discovers gene gain and loss during grapevine domestication. Plant Physiol. 2024; 195:1401-13

[54] Wicker T, Yahiaoui N, Keller B. Illegitimate recombination is a major evolutionary mechanism for initiating size variation in plant resistance genes. Plant Journal. 2007; 51:631-41

[55] Barragan AC, Weigel D. Plant NLR diversity: the known unknowns of pan-NLRomes. Plant cell. 2021; 33:814-31

[56] Ngou BPM, Ding P, Jones JDG. Thirty years of resistance: zig-Zag through the plant immune system. Plant cell. 2022; 34:1447-78

[57] Zhang W, Yan J, Li X. et al. Transcriptional response of grapevine to infection with the fungal pathogen Lasiodiplodia theobromae. Sci Rep. 2019; 9:5387

[58] Wang J, Tian W, Tao F. et al. TaRPM1 positively regulates wheat high-temperature seedling-plant resistance to Puccinia striiformis f.sp. tritici Front Plant Sci. 2019; 10:1679

[59] Burton JN, Adey A, Patwardhan RP. et al. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat Biotechnol. 2013; 31:1119-25

[60] Tarailo-Graovac M, Chen N. Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinfor-matics. 2009; 25:4.10.1-14

[61] Tang S, Lomsadze A, Borodovsky M. Identification of protein coding regions in RNA transcripts. Nucleic Acids Res. 2015; 43:e78-8

[62] Haas BJ, Salzberg SL, Zhu W. et al. Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments. Genome Biol. 2008;9:R7

[63] Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detec-tion of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997; 25:955-64

[64] Altschul SF, Gish W, Miller W. et al. Basic local alignment search tool. J Mol Biol. 1990; 215:403-10

[65] Conesa A, Götz S, García-Gómez JM. et al. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005; 21:3674-6

[66] Zdobnov EM, Apweiler R. InterProScan - an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001; 17:847-8

[67] Wang Y, Li J, Paterson AH. MCScanX-transposed: detecting transposed gene duplications based on multiple colinearity scans. Bioinformatics. 2013; 29:1458-60

[68] Emms DM, Kelly S. OrthoFinder: phylogenetic orthology infer-ence for comparative genomics. Genome Biol. 2019; 20:238

[69] Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007; 24:1586-91

[70] Han MV, Thomas GWC, Lugo-Martinez J. et al. Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. MolBiolEvol. 2013; 30:1987-97

[71] Nguyen L-T, Schmidt HA, von Haeseler A. et al. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015; 32:268-74

[72] Wu TD, Watanabe CK. GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics. 2005; 21: 1859-75

[73] van Wersch S, Li X. Stronger when together: clustering of plant NLR disease resistance genes. Trends Plant Sci. 2019; 24:688-99

[74] Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015; 12:357-60

[75] Wang, Liu Y, He P. et al. Evaluation of foliar resistance to Uncin-ula necator in chines wild Vitis species. Vitis. 2015; 34:159

[76] Kendig KI, Baheti S, Bockol MA. et al. Sentieon DNASeq variant calling workflow demonstrates strong computational perfor-mance and accuracy. Front Genet. 2019; 10:736

[77] Purcell S, Neale B, Todd-Brown K. et al. PLINK: a tool set for whole-genome association and population-based linkage anal-yses. The American Journal of Human Genetics. 2007; 81:559-75

[78] Broman KW, Wu H, Sen Ś. et al. R/QTL: QTL mapping in experi-mental crosses. Bioinformatics. 2003; 19:889-90

[79] Neff MM, Turk E, Kalishman M. Web-based primer design for single nucleotide polymorphism analysis. Trends Genet. 2002; 18: 613-5