Comparative transcriptomic analysis highlights contrasting levels of resistance of Vitis vinifera and Vitis amurensis to Botrytis cinerea

Abstract

Botrytis cinerea is a major grapevine (Vitis spp.) pathogen, but some genotypes differ in their degree of resistance. For example, the Vitis vinifera cultivar Red Globe (RG) is highly susceptible, but V. amurensis Rupr Shuangyou (SY) is highly resistant. Here, we used RNA sequencing analysis to characterize the transcriptome responses of these two genotypes to B. cinerea inoculation at an early infection stage. Approximately a quarter of the genes in RG presented significant changes in transcript levels during infection, the number of which was greater than that in the SY leaves. The genes differentially expressed between infected leaves of SY and RG included those associated with cell surface structure, oxidation, cell death and C/N metabolism. We found evidence that an imbalance in the levels of reactive oxygen species (ROS) and redox homeostasis probably contributed to the susceptibility of RG to B. cinerea. SY leaves had strong antioxidant capacities and improved ROS homeostasis following infection. Regulatory network prediction suggested that WRKY and MYB transcription factors are associated with the abscisic acid pathway. Weighted gene correlation network analysis highlighted preinfection features of SY that might contribute to its increased resistance. Moreover, overexpression of VaWRKY10 in Arabidopsis thaliana and V. vinifera Thompson Seedless enhanced resistance to B. cinerea. Collectively, our study provides a high-resolution view of the transcriptional changes of grapevine in response to B. cinerea infection and novel insights into the underlying resistance mechanisms.

References

[1] Walker, A. S. et al. French vineyards provide information that opens ways for effective resistance management of Botrytis cinerea (grey mould). Pest Manag. Sci. 69, 667-678 (2013).

[2] Bethke, G. et al. Pectin biosynthesis is critical for cell wall integrity and immunity in Arabidopsis thaliana. Plant Cell 28, 537-556 (2016).

[3] Choquer, M. et al. Botrytis cinerea virulence factors: new insights into a necrotrophic and polyphageous pathogen. FEMS Microbiol. Lett. 277, 1-10 (2007).

[4] Castillo, L., Plaza, V., Larrondo, L. F. & Canessa, P. Recent advances in the study of the plant pathogenic fungus Botrytis cinerea and its interaction with the environment. Curr. Protein Pept. Sci. 17, 976-989 (2017).

[5] Govrin, E. M. & Levine, A. The hypersensitive response facilitates plant infection by the necrotrophic pathogen Botrytis cinerea. Curr. Biol. 10, 751-757 (2000).

[6] Veloso, J. & van Kan, J. A. L. Many shades of grey in Botrytis-host plant interactions. Trends Plant Sci. 23, 613-622 (2018).

[7] Wan, R. et al. Resistance evaluation of Chinese wild Vitis genotypes against Botrytis cinerea and different responses of resistant and susceptible hosts to the infection. Front. Plant Sci. 6, 854 (2015).

[8] Rahman, M. U. et al. Screening Vitis genotypes for responses to Botrytis cinerea and evaluation of antioxidant enzymes, reactive oxygen species and jasmonic acid in resistant and susceptible hosts. Molecules 24, 5 (2019).

[9] Song, S. et al. Interaction between MYC2 and ETHYLENE INSENSITIVE3 modulates antagonism between jasmonate and ethylene signaling in Arabidopsis. Plant Cell 26, 263-279 (2014).

[10] Audenaert, K., De Meyer, G. B. & Höfte, M. M. Abscisic acid determines basal susceptibility of tomato to Botrytis cinerea and suppresses salicylic acid-dependent signaling mechanisms. Plant Physiol. 128, 491-501 (2002).

[11] Birkenbihl, R. P., Diezel, C. & Somssich, I. E. Arabidopsis WRKY33 is a key transcriptional regulator of hormonal and metabolic responses toward Botrytis cinerea infection. Plant Physiol. 159, 266-285 (2012).

[12] Han, X. et al. Regulation of GDSL lipase gene expression by the MPK3/MPK6 cascade and its downstream WRKY transcription factors in Arabidopsis immunity. Mol. Plant Microbe Interact. 32, 673-684 (2019).

[13] Liu, S., Kracher, B., Ziegler, J., Birkenbihl, R. P. & Somssich, I. E. Negative regulation of ABA signaling by WRKY33 is critical for Arabidopsis immunity towards Botrytis cinerea 2100. Elife Sci. 4, e07295 (2015).

[14] Jiang, Y. & Yu, D. WRKY57 regulates JAZ genes transcriptionally to compromise Botrytis cinerea resistance in Arabidopsis thaliana. Plant Physiol. 171, 2771-2782 (2016).

[15] Windram, O. et al. Arabidopsis defense against Botrytis cinerea: chronology and regulation deciphered by high-resolution temporal transcriptomic analysis. Plant Cell 24, 3530-3557 (2012).

[16] AbuQamar, S. F., Moustafa, K. & Tran, L. S. P. ‘Omics’ and plant responses to Botrytis cinerea. Front. Plant Sci. 7, 1658- 1658 (2016).

[17] Agudelo-Romero, P. et al. Transcriptome and metabolome reprogramming in Vitis vinifera cv. Trincadeira berries upon infection with Botrytis cinerea. J. Exp. Bot. 66, 3714-3717 (2015).

[18] Blanco-Ulate, B. et al. Developmental and metabolic plasticity of white-skinned grape berries in response to Botrytis cinerea during noble rot. Plant Physiol. 169, 2422-2443 (2015).

[19] Haile, Z. M. et al. Molecular analysis of the early interaction between the grapevine flower and Botrytis cinerea reveals that prompt activation of specific host pathways leads to fungus quiescence. Plant Cell Environ. 40, 1409-1428 (2017).

[20] Kelloniemi, J. et al. Analysis of the molecular dialogue between gray mold (Botrytis cinerea) and grapevine (Vitis vinifera) reveals a clear shift in defense mechanisms during berry ripening. Mol. Plant Microbe Interact. 28, 1167-1180 (2015).

[21] Ramirez, V., Garcia-Andrade, J. & Vera, P. Enhanced disease resistance to Botrytis cinerea in myb46 Arabidopsis plants is associated to an early down-regulation of CesA genes. Plant Signal. Behav. 6, 911-913 (2011).

[22] Sharma, P., Jha, A. B., Dubey, R. S. & Pessarakli, M. Reactive oxygen species, oxidative damage, and antioxidative defense mechanism in plants under stressful conditions. J. Bot. 2012, 217037 (2012).

[23] Wang, P. et al. Genome-wide characterization of bHLH genes in grape and analysis of their potential relevance to abiotic stress tolerance and secondary metabolite biosynthesis. Front. Plant Sci. 9, 64 (2018).

[24] Pattison, R. J., Csukasi, F., Zheng, Y. & Fei, Z. Comprehensive tissue-specific transcriptome analysis reveals distinct regulatory programs during early tomato fruit development. Plant Physiol. 168, 1684-1701 (2015).

[25] Buchanan, B., Wilhelm, G. & Russell, L. J. Plant Biochemistry and Molecular Biology (Chinese version) (Science press, 2004).

[26] Seifi, H. S., Van, B. J., Angenon, G. & Höfte, M. Glutamate metabolism in plant disease and defense: friend or foe? Mol. Plant Microbe Interact. 26, 475-485 (2013).

[27] Langfelder, P. & Horvath, S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 9, 559 (2008).

[28] Guo, C. et al. Evolution and expression analysis of the grape (Vitis vinifera L.) WRKY gene family. J. Exp. Bot. 65, 1513-1528 (2014).

[29] Xu, X., Chen, C., Fan, B. & Chen, Z. Physical and functional interactions between pathogen-induced Arabidopsis WRKY18, WRKY40, and WRKY60 transcription factors. Plant Cell 18, 1310-1326 (2006).

[30] Moffett, P., Farnham, G., Peart, J. & Baulcombe, D. C. Interaction between domains of a plant NBS-LRR protein in disease resistance-related cell death. EMBO J. 21, 4511-4519 (2002).

[31] Tian, S., Yin, X., Fu, P., Wu, W. & Lu, J. Ectopic Expression of grapevine gene VaRGA1 in Arabidopsis improves resistance to downy mildew and Pseudomonas syringae pv. tomato DC3000 but increases susceptibility to Botrytis cinerea. Inter. J. Mol. Sci. 21, 193 (2020).

[32] Ye, J. et al. Knockdown of SlNL33 accumulates ascorbate, enhances disease and oxidative stress tolerance in tomato (Solanum lycopersicum). Plant Growth Regul. 89, 49-58 (2019).

[33] Massimo, A. et al. Current and novel aspects on the non-lysosomal β-glucosylceramidase GBA2. Neurochem. Res. 41, 210-220 (2016).

[34] Michaelson, L. V., Napier, J. A., Molino, D. & Faure, J.-D. Plant sphingolipids: their importance in cellular organization and adaption. Bba. Mol. 1861, 1329-1335 (2016).

[35] Reiss, U. et al. Sphingosine-phosphate lyase enhances stress-induced ceramide generation and apoptosis. J. Biol. Chem. 279, 1281-1290 (2004).

[36] Deighton, N., Muckenschnabel, I., Goodman, B. A. & Williamson, B. Lipid peroxidation and the oxidative burst associated with infection of Capsicum annuum by Botrytis cinerea. Plant J. 20, 485-492 (1999).

[37] Cremer, K. D. et al. RNAseq-based transcriptome analysis of Lactuca sativa infected by the fungal necrotroph Botrytis cinerea. Plant Cell Environ. 36, 1992-2007 (2013).

[38] Wang, X. et al. CRISPR/Cas9-mediated efficient targeted mutagenesis in grape in the first generation. Plant Biotechnol. J. 16, 844-855 (2018).

[39] Rubio, S. et al. Triple loss of function of protein phosphatases type 2C leads to partial constitutive response to endogenous abscisic acid. Plant Physiol. 150, 1345-1355 (2009).

[40] Breia, R. et al. VvSWEET7 is a mono- and disaccharide transporter up-regulated in response to Botrytis cinerea infection in grape berries. Front. Plant. Sci. 10, 1753 (2020).

[41] Seifi, H. S. et al. Concurrent overactivation of the cytosolic glutamine synthetase and the GABA shunt in the ABA-deficient sitiens mutant of tomato leads to resistance against Botrytis cinerea. New Phytol. 199, 490-504 (2013).

[42] Hatmi, S. et al. Osmotic stress-induced polyamine oxidation mediates defence responses and reduces stress-enhanced grapevine susceptibility to Botrytis cinerea. J. Exp. Bot. 65, 75-88 (2014).

[43] Lecompte, F., Abro, M. A. & Nicot, P. C. Contrasted responses of Botrytis cinerea isolates developing on tomato plants grown under different nitrogen nutrition regimes. Plant Pathol. 59, 891-899 (2010).

[44] Li, D. et al. Ectopic expression of rice OsBIANK1, encoding an ankyrin repeat-containing protein, in Arabidopsis confers enhanced disease resistance to Botrytis cinerea and Pseudomonas syringae. J. Phytopathol. 161, 27-34 (2013).

[45] Yin, X. et al. Rice copine genes OsBON1 and OsBON3 function as suppressors of broad-spectrum disease resistance. Plant Biotechnol. J. 16, 1476-1487 (2018).

[46] Love, A. J., Milner, J. J. & Sadanandom, A. Timing is everything: Regulatory overlap in plant cell death. Trends Plant Sci. 13, 589-595 (2008).

[47] Birkenbihl, R. P., Kracher, B., Roccaro, M. & Somssich, I. E. Induced genome-wide binding of three Arabidopsis WRKY transcription factors during early MAMP-triggered immunity. Plant Cell 29, 20-38 (2017).

[48] Cheng, C. et al. Gibberellin-induced changes in the transcriptome of grapevine (Vitis labrusca × V. vinifera) cv. Kyoho flowers. BMC Genomics 16, 128 (2015).

[49] Gao, J. Experimental Instruction In Plant Physiology (Higher Education Press, 2006).

[50] Zhong, S. et al. High-throughput illumina strand-specific RNA sequencing library preparation. Cold Spring Harb. Protoc. 8, 940-949 (2011).

[51] Jiao, C., Gao, M., Wang, X. & Fei, Z. Transcriptome characterization of three wild Chinese Vitis uncovers a large number of distinct disease related genes. BMC Genomics 16, 223 (2015).

[52] Jaillon, O. et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature 449, 463-467 (2007).

[53] Staats, M. & van Kan, J. A. L. Genome update of Botrytis cinerea Strains B05.10 and T4. Eukaryot. Cell 11, 12 (2012).

[54] Kent, W. J. BLAT - The BLAST-like alignment tool. Genome Res. 12, 656-664 (2002).

[55] Anders, S. & Huber, W. Differential expression analysis for sequence count data. Genome Biol. 11, R106 (2010).

[56] Joung, J. G. et al. Plant MetGenMAP: an integrative analysis system for plant systems biology. Plant Physiol. 151, 1758-1768 (2009).

[57] Sturn, A., Quackenbush, J. & Trajanoski, Z. Genesis: cluster analysis of microarray data. Bioinformatics 18, 207-208 (2002).

[58] Hardison, R. C. & Taylor, J. Genomic approaches towards finding cis-regulatory modules in animals. Nat. Rev. Genet. 13, 469-483 (2012).

[59] Wasserman, W. W. & Sandelin, A. Applied bioinformatics for the identification of regulatory elements. Nat. Rev. Genet. 5, 276-287 (2004).

[60] Singh, V. K., Mangalam, A. K., Dwivedi, S. & Naik, S. Primer premier: program for design of degenerate primers from a protein sequence. Biotechniques 24, 318-319 (1998).

[61] Guo, R. et al. Ectopic expression of a grape aspartic protease gene, AP13, in Arabidopsis thaliana improves resistance to powdery mildew but increases susceptibility to Botrytis cinerea. Plant Sci. 248, 17-27 (2015).

[62] Clough, S. J. & Bent, A. F. Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 16, 735-743 (1998).