Transcriptional analysis and histochemistry reveal that hypersensitive cell death and H ₂O ₂ have crucial roles in the resistance of tea plant (Camellia sinensis (L.) O. Kuntze) to anthracnose

Abstract

Anthracnose causes severe losses of tea production in China. Although genes and biological processes involved in anthracnose resistance have been reported in other plants, the molecular response to anthracnose in tea plant is unknown. We used the susceptible tea cultivar Longjing 43 and the resistant cultivar Zhongcha 108 as materials and compared transcriptome changes in the leaves of both cultivars following Colletotrichum fructicola inoculation. In all, 9015 and 8624 genes were differentially expressed between the resistant and susceptible cultivars and their controls (0 h), respectively. In both cultivars, the differentially expressed genes (DEGs) were enriched in 215 pathways, including responses to sugar metabolism, phytohormones, reactive oxygen species (ROS), biotic stimuli and signalling, transmembrane transporter activity, protease activity and signalling receptor activity, but DEG expression levels were higher in Zhongcha 108 than in Longjing 43. Moreover, functional enrichment analysis of the DEGs showed that hydrogen peroxide (H₂O₂) metabolism, cell death, secondary metabolism, and carbohydrate metabolism are involved in the defence of Zhongcha 108, and 88 key genes were identified. Protein–protein interaction (PPI) network demonstrated that putative mitogen-activated protein kinase (MAPK) cascades are activated by resistance (R) genes and mediate downstream defence responses. Histochemical analysis subsequently validated the strong hypersensitive response (HR) and H₂O₂ accumulation that occurred around the hyphal infection sites in Zhongcha 108. Overall, our results indicate that the HR and H₂O₂ are critical mechanisms in tea plant defence against anthracnose and may be activated by R genes via MAPK cascades.

References

[1] Cabrera, C., Artacho, R. & Giménez, R. Beneficial effects of green tea-a review. J. Am. Coll. Nutr. 25, 79-99 (2006).

[2] Yang, C. S., Zhang, J., Zhang, L., Huang, J. & Wang, Y. Mechanisms of body weight reduction and metabolic syndrome alleviation by tea. Mol. Nutr. Food Res. 60, 160-174 (2016).

[3] Wang, Y. C. et al. Diverse Colletotrichum species cause anthracnose of tea plants (Camellia sinensis (L.) O. Kuntze) in China. Sci. Rep. 6, 35287 (2016).

[4] Wang, Y. C. et al. Metabolic changes of caffeine in tea plant (Camellia sinensis (L.) O. Kuntze) as defense response to Colletotrichum fructicola. J. Agric. Food Chem. 64, 6685-6693 (2016).

[5] Kunkel, B. N. & Brooks, D. M. Cross talk between signaling pathways in pathogen defense. Curr. Opin. Plant Biol. 5, 325-331 (2002).

[6] Schwessinger, B. & Zipfel, C. News from the frontline: recent insights into PAMP-triggered immunity in plants. Curr. Opin. Plant Biol. 11, 389-395 (2008).

[7] Thomma, B. P., Nurnberger, T. & Joosten, M. H. Of PAMPs and effectors: the blurred PTI-ETI dichotomy. Plant Cell 23, 4-15 (2011).

[8] Cui, H. T. K. & Parker, J. E. Effector-triggered immunity: from pathogen perception to robust defense. Annu. Rev. Plant Biol. 66, 487-511 (2015).

[9] Greenberg, J. T. & Yao, N. The role and regulation of programmed cell death in plant-pathogen interactions. Cell. Microbiol. 6, 201-211 (2004).

[10] Cannon, P. F., Damm, U., Johnston, P. R. & Weir, B. S. Colletotrichum-current status and future directions. Stud. Mycol. 73, 181-213 (2012).

[11] O’Connell, R. et al. A novel Arabidopsis-Colletotrichum pathosystem for the molecular dissection of plant-fungal interactions. Mol. Plant Microbe Interact. 17, 272-282 (2004).

[12] Bhadauria, V. et al. Identification of Lens culinaris defense genes responsive to the anthracnose pathogen Colletotrichum truncatum. BMC Genet. 14, 31 (2013).

[13] Padder, B. A., Kamfwa, K., Awale, H. E. & Kelly, J. D. Transcriptome profiling of the Phaseolus vulgaris-Colletotrichum lindemuthianum pathosystem. PLoS ONE 11, e0165823 (2016).

[14] Senthilkumar, P., Thirugnanasambantham, K. & Mandal, A. K. Suppressive subtractive hybridization approach revealed differential expression of hypersensitive response and reactive oxygen species production genes in tea (Camellia sinensis (L.) O. Kuntze) leaves during Pestalotiopsis thea infection. Appl. Biochem. Biotechnol. 168, 1917-1927 (2012).

[15] Palanisamy, S. & Mandal, A. K. Susceptibility against grey blight disease-causing fungus Pestalotiopsis sp. in tea (Camellia sinensis (L.) O. Kuntze) cultivars is influenced by anti-oxidative enzymes. Appl. Biochem. Biotechnol. 172, 216-223 (2014).

[16] Jayaswall, K. et al. Transcriptome analysis reveals candidate genes involved in blister blight defense in tea (Camellia sinensis (L) Kuntze). Sci. Rep. 6, 30412 (2016).

[17] Wang, L. et al. Transcriptome analysis of an anthracnose-resistant tea plant cultivar reveals genes associated with resistance to Colletotrichum camelliae. PLoS ONE 11, e0148535 (2016).

[18] Hao, X. et al. Comprehensive transcriptome analyses reveal differential gene expression profiles of Camellia sinensis axillary buds at para-, endo-, eco-dormancy, and bud flush stages. Front. Plant Sci. 8, 553 (2017).

[19] Grabherr, M. G. et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644-652 (2011).

[20] Li, B. & Dewey, C. N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 12, 323 (2011).

[21] Leng, N. et al. EBSeq-HMM: a Bayesian approach for identifying gene-expression changes in ordered RNA-seq experiments. Bioinformatics 31, 2614-2622 (2015).

[22] Zheng, Q. & Wang, X. J. GOEAST: a web-based software toolkit for Gene Ontology enrichment analysis. Nucleic Acids Res. 36, 358-363 (2008).

[23] Huang, D. W., Sherman, B. T. & Lempicki, R. A. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 37, 1-13 (2008).

[24] Ernst, J. & Barjoseph, Z. STEM: a tool for the analysis of short time series gene expression data. BMC Bioinformatics 7, 191 (2006).

[25] Shannon, P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498-2504 (2003).

[26] Faoro, F., Maffi, D., Cantu, D. & Iriti, M. Chemical-induced resistance against powdery mildew in barley: the effects of chitosan and benzothiadiazole. BioControl 53, 387-401 (2008).

[27] Hao, X. et al. Histochemical studies on the accumulation of H₂O₂ and hypersensitive cell death in the nonhost resistance of pepper against Blumeria graminis f. sp. tritici . Physiol. Mol. Plant Pathol. 76, 104-111 (2011).

[28] Hao, X. et al. Identification and evaluation of reliable reference genes for quantitative real-time PCR analysis in tea plant (Camellia sinensis (L.) O. Kuntze). Int. J. Mol. Sci. 15, 22155-22172 (2014).

[29] Schmittgen, T. D. & Livak, K. J. Analyzing real-time PCR data by the comparative C(T) method. Nat. Protoc. 3, 1101-1108 (2008).

[30] Ferrier-Cana, E. et al. Characterization of expressed NBS-LRR resistance gene candidates from common bean. Theor. Appl. Genet. 106, 251-261 (2003).

[31] Li, J. et al. Genome-wide identification and comparative expression analysis of NBS-LRR-encoding genes upon Colletotrichum gloeosporioides infection in two ecotypes of Fragaria vesca. Gene 527, 215-227 (2013).

[32] Takahara, H. et al. Colletotrichum higginsianum extracellular LysM proteins play dual roles in appressorial function and suppression of chitin-triggered plant immunity. New Phytol. 211, 1323-1337 (2016).

[33] Jones, D. A. & Takemoto, D. Plant innate immunity-direct and indirect recognition of general and specific pathogen-associated molecules. Curr. Opin. Immunol. 16, 48-62 (2004).

[34] Bonardi, V. & Dangl, J. L. How complex are intracellular immune receptor signaling complexes?. Front. Plant Sci. 3, 237 (2012).

[35] Monaghan, J. & Zipfel, C. Plant pattern recognition receptor complexes at the plasma membrane. Curr. Opin. Plant Biol. 15, 349-357 (2012).

[36] Zipfel, C. Plant pattern-recognition receptors. Trends Immunol. 35, 345-351 (2014).

[37] Jones, J. D. & Dangl, J. L. The plant immune system. Nature 444, 323-329 (2006).

[38] Tsuda, K. & Katagiri, F. Comparing signaling mechanisms engaged in pattern-triggered and effector-triggered immunity. Curr. Opin. Plant Biol. 13, 459-465 (2010).

[39] Meng, X. & Zhang, S. MAPK cascades in plant disease resistance signaling. Annu. Rev. Phytopathol. 51, 245-266 (2013).

[40] Steinhorst, L. & Kudla, J. Calcium and reactive oxygen species rule the waves of signaling. Plant Physiol. 163, 471-485 (2013).

[41] Li, J. et al. Proteomic analysis of the compatible interaction of wheat and powdery mildew (Blumeria graminis f. sp. tritici). Plant Physiol. Biochem. 111, 234-243 (2017).

[42] Lumbreras, V. et al. MAPK phosphatase MKP2 mediates disease responses in Arabidopsis and functionally interacts with MPK3 and MPK6. Plant J. 63, 1017-1030 (2010).

[43] Vilela, B., Pagès, M. & Lumbreras, V. Regulation of MAPK signaling and cell death by MAPK phosphatase MKP2. Plant Signal. Behav. 5, 1497-1500 (2010).

[44] Tenhaken, R., Doerks, T. & Bork, P. DCD-a novel plant specific domain in proteins involved in development and programmed cell death. BMC Bioinformatics 6, 169 (2005).

[45] Shi, L. et al. Involvement of sphingoid bases in mediating reactive oxygen intermediate production and programmed cell death in Arabidopsis. Cell Res. 17, 1030-1040 (2007).

[46] Gechev, T. S., Van Breusegem, F., Stone, J. M., Denev, I. & Laloi, C. Reactive oxygen species as signals that modulate plant stress responses and programmed cell death. Bioessays 28, 1091-1101 (2006).

[47] Daudi, A. et al. The apoplastic oxidative burst peroxidase in Arabidopsis is a major component of pattern-triggered immunity. Plant Cell 24, 275-287 (2012).

[48] Li, W. et al. A natural allele of a transcription factor in rice confers broad-spectrum blast resistance. Cell 170, 114-126 (2017).

[49] Choi, H. W. & Hwang, B. K. The pepper extracellular peroxidase CaPO₂ is required for salt, drought and oxidative stress tolerance as well as resistance to fungal pathogens . Planta 235, 1369-1382 (2012).

[50] Ascencio-Ibáñez, J. T. et al. Global analysis of Arabidopsis gene expression uncovers a complex array of changes impacting pathogen response and cell cycle during geminivirus infection. Plant Physiol. 148, 436-454 (2008).

[51] Rodríguez-Serrano, M., Romero-Puertas, M. C., Sanz-Fernández, M., Hu, J. & Sandalio, L. M. Peroxisomes extend peroxules in a fast response to stress via a reactive oxygen species-mediated induction of the peroxin PEX11a. Plant Physiol. 171, 1665-1674 (2016).

[52] Vellosillo, T., Vicente, J., Kulasekaran, S., Hamberg, M. & Castresana, C. Emerging complexity in reactive oxygen species production and signaling during the response of plants to pathogens. Plant Physiol. 154, 444-448 (2010).

[53] Kubicek, C. P., Starr, T. L. & Glass, N. L. Plant cell wall-degrading enzymes and their secretion in plant-pathogenic fungi. Annu. Rev. Phytopathol. 52, 427-451 (2014).

[54] Trouvelot, S. et al. Carbohydrates in plant immunity and plant protection: roles and potential application as foliar sprays. Front. Plant Sci. 5, 592 (2014).

[55] Cui, L. et al. Identification of UDP-glycosyltransferases involved in the biosynthesis of astringent taste compounds in tea (Camellia sinensis). J. Exp. Bot. 67, 2285-2297 (2016).

[56] Figueiredo, A. et al. Transcriptional and metabolic profiling of grape (Vitis vinifera L.) leaves unravel possible innate resistance against pathogenic fungi. J. Exp. Bot. 59, 3371-3381 (2008).

[57] Trépanier, M. et al. Dependence of arbuscular-mycorrhizal fungi on their plant host for palmitic acid synthesis. Appl. Environ. Microbiol. 71, 5341-5347 (2005).

[58] Jiang, Y. et al. Plants transfer lipids to sustain colonization by mutualistic mycorrhizal and parasitic fungi. Science 356, 1172-1175 (2017).

[59] Rubio, M. et al. Analysis of gene expression changes in peach leaves in response to Plum pox virus infection using RNA-Seq. Mol. Plant Pathol. 16, 164-176 (2015).

[60] Guo, W. et al. An ethylene response-related factor, GbERF1-like, from Gossypium barbadense improves resistance to Verticillium dahliae via activating lignin synthesis. Plant Mol. Biol. 91, 305-318 (2016).

[61] Becker, M. G. et al. Transcriptome analysis of the Brassica napus-Leptosphaeria maculans pathosystem identifies receptor, signaling and structural genes underlying plant resistance. Plant J. 90, 573-586 (2017).