A temporal gene expression map of Chrysanthemum leaves infected with Alternaria alternata reveals different stages of defense mechanisms

Abstract

Chrysanthemum (Chrysanthemum morifolium) black spot disease (CBS) poses a major threat to Chrysanthemum cultivation owing to suitable climate conditions and current lack of resistant cultivars for greenhouse cultivation. In this study, we identified a number of genes that respond to Alternaria alternata infection in resistant and susceptible Chrysanthemum cultivars. Based on RNA sequencing technology and a weighted gene coexpression network analysis (WGCNA), we constructed a model to elucidate the response of Chrysanthemum leaves to A. alternata infection at different stages and compared the mapped response of the resistant cultivar ‘Jinba’ to that of the susceptible cultivar ‘Zaoyihong’. In the early stage of infection, when lesions had not yet formed, abscisic acid (ABA), salicylic acid (SA) and EDS1-mediated resistance played important roles in the Chrysanthemum defense system. With the formation of necrotic lesions, ethylene (ET) metabolism and the Ca²⁺ signal transduction pathway strongly responded to A. alternata infection. During the late stage, when necrotic lesions continued to expand, members of the multidrug and toxic compound extrusion (MATE) gene family were highly expressed, and their products may be involved in defense against A. alternata invasion by exporting toxins produced by the pathogen, which plays important roles in the pathogenicity of A. alternata. Furthermore, the function of hub genes was verified by qPCR and transgenic assays. The identification of hub genes at different stages, the comparison of hub genes between the two cultivars and the highly expressed genes in the resistant cultivar ‘Jinba’ provide a theoretical basis for breeding cultivars resistant to CBS.

References

[1] Thomma, B. P. Alternaria spp: from general saprophyte to specific parasite. Mol. Plant Pathol. 4, 225-236 (2003).

[2] Takaoka, S. et al. Complex regulation of secondary metabolism controlling pathogenicity in the phytopathogenic fungus Alternaria alternata. N. Phytol. 202, 1297-1309 (2014).

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

[4] Zipfel, C. Early molecular events in PAMP-triggered immunity. Curr. Opin. Plant Biol. 12, 414-420 (2009).

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

[6] Li, B., Meng, X., Shan, L. & He, P. Transcriptional regulation of pattern-triggered immunity in plants. Cell Host Microbe 19, 641-650 (2016).

[7] Nimchuk, Z., Eulgem, T., Holt, B. F. 3rd & Dangl, J. L. Recognition and response in the plant immune system. Annu. Rev. Genet. 37, 579-609 (2003).

[8] Collier, S. M. & Moffett, P. NB-LRRs work a "bait and switch" on pathogens. Trends Plant Sci. 14, 521-529 (2009).

[9] Le Roux, C. et al. A receptor pair with an integrated decoy converts pathogen disabling of transcription factors to immunity. Cell 161, 1074-1088 (2015).

[10] Vlot, A. C., Dempsey, D. A. & Klessig, D. F. Salicylic acid, a multifaceted hormone to combat disease. Annu. Rev. Phytopathol. 47, 177-206 (2009).

[11] Baggs, E., Dagdas, G. & Krasileva, K. V. NLR diversity, helpers and integrated domains: making sense of the NLR identity. Curr. Opin. Plant Biol. 38, 59-67 (2017).

[12] Zipfel, C. et al. Bacterial disease resistance in Arabidopsis through flagellin perception. Nature 428, 764-767 (2004).

[13] Wang, L. et al. Arabidopsis CaM binding protein CBP60g contributes to MAMP-induced SA accumulation and is involved in disease resistance against Pseudomonas syringae. PLoS Pathog. 5, e1000301 (2009).

[14] van Loon, L. C., Geraats, B. P. & Linthorst, H. J. Ethylene as a modulator of disease resistance in plants. Trends Plant Sci. 11, 184-191 (2006).

[15] Zhang, L. & Xing, D. Methyl jasmonate induces production of reactive oxygen species and alterations in mitochondrial dynamics that precede photosynthetic dysfunction and subsequent cell death. Plant Cell Physiol. 49, 1092-1111 (2008).

[16] Broekgaarden, C., Caarls, L., Vos, I. A., Pieterse, C. M. & Van Wees, S. C. Ethylene: traffic controller on hormonal crossroads to defense. Plant Physiol. 169, 2371-2379 (2015).

[17] Luo, J. et al. Integrated transcriptome analysis reveals plant hormones jasmonic acid and salicylic acid coordinate growth and defense responses upon fungal infection in poplar. Biomolecules 9, https://doi.org/10.3390/biom9010012 (2019).

[18] Adie, B. A. et al. ABA is an essential signal for plant resistance to pathogens affecting JA biosynthesis and the activation of defenses in Arabidopsis. Plant Cell 19, 1665-1681 (2007).

[19] Ton, J., Flors, V. & Mauch-Mani, B. The multifaceted role of ABA in disease resistance. Trends Plant Sci. 14, 310-317 (2009).

[20] Cai, S. et al. Evolutionary conservation of ABA signaling for stomatal closure. Plant Physiol. 174, 732-747 (2017).

[21] Ohme-Takagi, M., Suzuki, K. & Shinshi, H. Regulation of ethylene-induced transcription of defense genes. Plant Cell Physiol. 41, 1187-1192 (2000).

[22] Singh, K., Foley, R. C. & Onate-Sanchez, L. Transcription factors in plant defense and stress responses. Curr. Opin. Plant Biol. 5, 430-436 (2002).

[23] Dolmetsch, R. E., Xu, K. & Lewis, R. S. Calcium oscillations increase the efficiency and specificity of gene expression. Nature 392, 933-936 (1998).

[24] Poovaiah, B. W., Du, L., Wang, H. & Yang, T. Recent advances in calcium/calmodulin-mediated signaling with an emphasis on plant-microbe interactions. Plant Physiol. 163, 531-542 (2013).

[25] Gust, A. A. & Felix, G. Receptor like proteins associate with SOBIR1-type of adaptors to form bimolecular receptor kinases. Curr. Opin. Plant Biol. 21, 104-111 (2014).

[26] Yuan, P., Jauregui, E., Du, L., Tanaka, K. & Poovaiah, B. W. Calcium signatures and signaling events orchestrate plant-microbe interactions. Curr. Opin. Plant Biol. 38, 173-183 (2017).

[27] Zipfel, C. & Oldroyd, G. E. Plant signalling in symbiosis and immunity. Nature 543, 328-336 (2017).

[28] Edel, K. H., Marchadier, E., Brownlee, C., Kudla, J. & Hetherington, A. M. The evolution of calcium-based signalling in plants. Curr. Biol. 27, R667-R679 (2017).

[29] Brandt, B. et al. Calcium specificity signaling mechanisms in abscisic acid signal transduction in Arabidopsis guard cells. Elife 4, https://doi.org/10.7554/eLife.03599 (2015).

[30] Raz, V. & Fluhr, R. Calcium requirement for ethylene-dependent responses. Plant cell 4, 1123-1130 (1992).

[31] Reddy, A. S., Ali, G. S., Celesnik, H. & Day, I. S. Coping with stresses: roles of calcium- and calcium/calmodulin-regulated gene expression. Plant Cell 23, 2010-2032 (2011).

[32] Zhao, N. et al. Comparative Transcriptome analysis of waterlogging-sensitive and waterlogging-tolerant Chrysanthemum morifolium cultivars under waterlogging stress and reoxygenation conditions. Int. J. Mol. Sci. 19, https://doi.org/10.3390/ijms19051455 (2018).

[33] Wang, H., Lin, J., Chang, Y. & Jiang, C. Z. Comparative transcriptomic analysis reveals that ethylene/H2O2-mediated hypersensitive response and programmed cell death determine the compatible interaction of sand pear and Alternaria alternata. Front. Plant Sci. 8, 195 (2017).

[34] Wirthmueller, L., Zhang, Y., Jones, J. D. & Parker, J. E. Nuclear accumulation of the Arabidopsis immune receptor RPS4 is necessary for triggering EDS1-dependent defense. Curr. Biol. 17, 2023-2029 (2007).

[35] Heidrich, K. et al. Arabidopsis EDS1 connects pathogen effector recognition to cell compartment-specific immune responses. Science 334, 1401-1404 (2011).

[36] Huang, P. Y., Catinot, J. & Zimmerli, L. Ethylene response factors in Arabidopsis immunity. J. Exp. Bot. 67, 1231-1241 (2016).

[37] Tsuda, K., Sato, M., Stoddard, T., Glazebrook, J. & Katagiri, F. Network properties of robust immunity in plants. PLoS Genet. 5, e1000772 (2009).

[38] Cui, H. et al. Antagonism of transcription factor MYC2 by EDS1/PAD4 complexes bolsters salicylic acid defense in Arabidopsis effector-triggered immunity. Mol. Plant 11, 1053-1066 (2018).

[39] Lai, Z., Vinod, K., Zheng, Z., Fan, B. & Chen, Z. Roles of Arabidopsis WRKY3 and WRKY4 transcription factors in plant responses to pathogens. BMC Plant Biol. 8, 68 (2008).

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

[41] 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 4, e07295 (2015).

[42] Liu, S., Ziegler, J., Zeier, J., Birkenbihl, R. P. & Somssich, I. E. Botrytis cinerea B05.10 promotes disease development in Arabidopsis by suppressing WRKY33-mediated host immunity. Plant Cell Environ. 40, 2189-2206 (2017).

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

[44] Brown, R. L., Kazan, K., McGrath, K. C., Maclean, D. J. & Manners, J. M. A role for the GCC-box in jasmonate-mediated activation of the PDF1.2 gene of Arabidopsis. Plant Physiol. 132, 1020-1032 (2003).

[45] Meng, X. et al. Phosphorylation of an ERF transcription factor by Arabidopsis MPK3/MPK6 regulates plant defense gene induction and fungal resistance. Plant Cell 25, 1126-1142 (2013).

[46] Li, G. et al. Dual-level regulation of ACC synthase activity by MPK3/MPK6 cascade and its downstream WRKY transcription factor during ethylene induction in Arabidopsis. PLoS Genet. 8, e1002767 (2012).

[47] Dubois, M., Van den Broeck, L. & Inze, D. The pivotal role of ethylene in plant growth. Trends Plant Sci. 23, 311-323 (2018).

[48] Ma, W. et al. Ca2+, cAMP, and transduction of non-self perception during plant immune responses. Proc. Natl Acad. Sci. USA 106, 20995-21000 (2009).

[49] He, X. et al. Structure of a cation-bound multidrug and toxic compound extrusion transporter. Nature 467, 991-994 (2010).

[50] Lu, P. et al. Genome-wide analysis of multidrug and toxic compound extrusion (MATE) family in Gossypium raimondii and Gossypium arboreum and its expression analysis under salt, cadmium, and drought stress. G3 8, 2483-2500 (2018).

[51] Liu, J., Li, Y., Wang, W., Gai, J. & Li, Y. Genome-wide analysis of MATE transporters and expression patterns of a subgroup of MATE genes in response to aluminum toxicity in soybean. BMC Genomics 17, 223 (2016).

[52] Li, N. et al. Genome-wide analysis of MATE transporters and molecular characterization of aluminum resistance in Populus. J. Exp. Bot. 68, 5669-5683 (2017).

[53] Zhang, H. et al. A DTX/MATE-type transporter facilitates abscisic acid efflux and modulates ABA sensitivity and drought tolerance in Arabidopsis. Mol. Plant 7, 1522-1532 (2014).

[54] Estiarte, N., Lawrence, C. B., Sanchis, V., Ramos, A. J. & Crespo-Sempere, A. LaeA and VeA are involved in growth morphology, asexual development, and mycotoxin production in Alternaria alternata. Int. J. Food Microbiol. 238, 153-164 (2016).

[55] Sanchez-Lopez, A. M. et al. Arabidopsis responds to Alternaria alternata volatiles by triggering plastid phosphoglucose isomerase-independent mechanisms. Plant Physiol. 172, 1989-2001 (2016).

[56] Gu, Y. Q. et al. Tomato transcription factors pti4, pti5, and pti6 activate defense responses when expressed in Arabidopsis. Plant Cell 14, 817-831 (2002).

[57] Li, H. et al. RNA-Seq derived identification of differential transcription in the Chrysanthemum leaf following inoculation with Alternaria tenuissima. BMC Genomics 15, 9 (2014).

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

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

[60] Trapnell, C. et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511-515 (2010).

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

[62] Maere, S., Heymans, K. & Kuiper, M. BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 21, 3448-3449 (2005).

[63] Thimm, O. et al. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 37, 914-939 (2004).

[64] Zhang, B. & Horvath, S. A general framework for weighted gene co-expression network analysis. Stat. Appl. Genet. Mol. Biol. 4, Article17 (2005).

[65] Stuart, J. M., Segal, E., Koller, D. & Kim, S. K. A gene-coexpression network for global discovery of conserved genetic modules. Science 302, 249-255 (2003).

[66] Song, A. et al. Comprehensive analysis of mitogen-activated protein kinase cascades in Chrysanthemum. PeerJ 6, e5037 (2018).

[67] Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25, 402-408 (2001).

[68] Li, P. et al. The over-expression of a Chrysanthemum WRKY transcription factor enhances aphid resistance. Plant Physiol. Biochem. 95, 26-34 (2015).