Transcriptome reprogramming, epigenetic modifications and alternative splicing orchestrate the tomato root response to the beneficial fungus Trichoderma harzianum

Abstract

Beneficial interactions of rhizosphere microorganisms are widely exploited for plant biofertilization and mitigation of biotic and abiotic constraints. To provide new insights into the onset of the roots–beneficial microorganisms interplay, we characterised the transcriptomes expressed in tomato roots at 24, 48 and 72 h post inoculation with the beneficial fungus Trichoderma harzianum T22 and analysed the epigenetic and post-trascriptional regulation mechanisms. We detected 1243 tomato transcripts that were differentially expressed between Trichoderma-interacting and control roots and 83 T. harzianum transcripts that were differentially expressed between the three experimental time points. Interaction with Trichoderma triggered a transcriptional response mainly ascribable to signal recognition and transduction, stress response, transcriptional regulation and transport. In tomato roots, salicylic acid, and not jasmonate, appears to have a prominent role in orchestrating the interplay with this beneficial strain. Differential regulation of many nutrient transporter genes indicated a strong effect on plant nutrition processes, which, together with the possible modifications in root architecture triggered by ethylene/indole-3-acetic acid signalling at 72 h post inoculation may concur to the well-described growth-promotion ability of this strain. Alongside, T. harzianum-induced defence priming and stress tolerance may be mediated by the induction of reactive oxygen species, detoxification and defence genes. A deeper insight into gene expression and regulation control provided first evidences for the involvement of cytosine methylation and alternative splicing mechanisms in the plant– Trichoderma interaction. A model is proposed that integrates the plant transcriptomic responses in the roots, where interaction between the plant and beneficial rhizosphere microorganisms occurs.

References

[1] Lopez-Bucio, J., Pelagio-Flores, R. & Herrera-Estrella, A. Trichoderma as biostimulant: exploiting the multilevel properties of a plant beneficial fungus. Sci. Hortic. 196 109-123 (2015).

[2] Rouphael, Y., Cardarelli, M., Bonini, P. & Colla, G. Synergistic action of a microbial-based biostimulant and a plant derived-protein hydrolysate enhances lettuce tolerance to alkalinity and salinity. Front Plant Sci. 8, 131 (2017).

[3] Shoresh, M., Harman, G.E. & Mastouri, F. Induced systemic resistance and plant responses to fungal biocontrol agents. Annu. Rev. Phytopathol. 48, 21-43 (2010).

[4] Tucci, M., Ruocco, M., De Masi, L., De Palma, M. & Lorito, M. The beneficial effect of Trichoderma spp. on tomato is modulated by the plant genotype. Mol. Plant Pathol. 12, 341-354 (2011).

[5] Hermosa, R., Viterbo, A., Chet, I. & Monte, E. Plant-beneficial effects of Trichoderma and of its genes. Microbiology 158, 17-25 (2012).

[6] Pieterse, C. M. et al. Induced systemic resistance by beneficial microbes. Annu. Rev. Phytopathol. 52, 347-375 (2014).

[7] Martinez-Medina, A. et al. Shifting from priming of salicylic acid- to jasmonic acid-regulated defences by Trichoderma protects tomato against the root knot nematode Meloidogyne incognita. New Phytol. 213, 1363-1377 (2017).

[8] Brotman, Y. et al. Trichoderma-plant root colonization: escaping early plant defense responses and activation of the antioxidant machinery for saline stress tolerance. PLoS Pathog. 9, e1003221 (2013).

[9] Zhang, S. W., Gan, Y. T. & Xu, B. L. Application of plant-growth-promoting fungi Trichoderma longibrachiatum T6 enhances tolerance of wheat to salt stress through improvement of antioxidative defense system and gene expression. Front. Plant Sci. 7, 1405 (2016).

[10] Fu, J., Liu, Z., Li, Z., Wang, Y. & Yang, K. Alleviation of the effects of saline-alkaline stress on maize seedlings by regulation of active oxygen metabolism by Trichoderma asperellum. PLoS ONE 12, e0179617 (2017).

[11] Alfano, G. et al. Systemic modulation of gene expression in tomato by Trichoderma hamatum 382. Phytopathology 97, 429-437 (2007).

[12] Shoresh, M . & Harman, G. E. The molecular basis of shoot responses of maize seedlings to Trichoderma harzianum T22 inoculation of the root: a proteomic approach. Plant Physiol. 147, 2147-2163 (2008).

[13] Mathys, J. et al. Genome-wide characterization of ISR induced in Arabidopsis thaliana by Trichoderma hamatum T382 against Botrytis cinerea infection. Front. Plant Sci. 3, 108 (2012).

[14] De Palma, M. et al. Suppression subtractive hybridization analysis provides new insights into the tomato (Solanum lycopersicum L.) response to the plant probiotic microorganism Trichoderma longibrachiatum MK1. J. Plant Physiol. 190, 79-94 (2016).

[15] Hermosa, R. et al. The contribution of Trichoderma to balancing the costs of plant growth and defense. Int Microbiol 16, 69-80 (2013).

[16] Ruocco, M. et al. Multiple roles and effects of a novel Trichoderma hydrophobin. Mol. Plant-Microbe Interact. 28, 167-179 (2015).

[17] Vos, C. M., De Cremer, K., Cammue, B. P. & De Coninck, B. The toolbox of Trichoderma spp. in the biocontrol of Botrytis cinerea disease. Mol. Plant Pathol. 16, 400-412 (2015).

[18] Lace, B. et al. Gate crashing arbuscular mycorrhizas: in vivo imaging shows the extensive colonization of both symbionts by Trichoderma atroviride. Environ. Microbiol. Rep. 7, 64-77 (2015).

[19] Contreras-Cornejo, H. A., Macias-Rodriguez, L., Cortes-Penagos, C. & Lopez-Bucio, J. Trichoderma virens, a plant beneficial fungus, enhances biomass production and promotes lateral root growth through an auxin-dependent mechanism in Arabidopsis. Plant Physiol. 149, 1579-1592 (2009).

[20] Howard, B. E. et al. High-throughput RNA sequencing of pseudomonas-infected Arabidopsis reveals hidden transcriptome complexity and novel splice variants. PLoS ONE 8, e74183 (2013).

[21] Zogli, P. & Libault, M. Plant response to biotic stress: Is there a common epigenetic response during plant-pathogenic and symbiotic interactions?. Plant Sci. 263, 89-93 (2017).

[22] Da, K., Nowak, J. & Flinn, B. Potato cytosine methylation and gene expression changes induced by a beneficial bacterial endophyte, Burkholderia phytofirmans strain PsJN. Plant Physiol. Biochem. 50, 24-34 (2012).

[23] Chamala, S., Feng, G., Chavarro, C. & Barbazuk, W. B. Genome-wide identification of evolutionarily conserved alternative splicing events in flowering plants. Front Bioeng. Biotechnol. 3, 33 (2015).

[24] Yang, S., Tang, F. & Zhu, H. Alternative splicing in plant immunity. Int J. Mol. Sci. 15, 10424-10445 (2014).

[25] Tan, G. X. et al. Transcriptome analysis of the compatible interaction of tomato with Verticillium dahliae using RNA-sequencing. Front. Plant Sci. 6, 428 (2015).

[26] Ruzicka, D. et al. Inside arbuscular mycorrhizal roots - molecular probes to understand the symbiosis. Plant Genome 6, https://doi.org/10.3835/plantgenome2012.06.0007 (2013).

[27] Sun, Y. & Xiao, H. Identification of alternative splicing events by RNA sequencing in early growth tomato fruits. BMC Genom. 16, 948 (2015).

[28] Yedidia, I. I., Benhamou, N. & Chet, I. I. Induction of defense responses in cucumber plants (Cucumis sativus L.) by the biocontrol agent Trichoderma harzianum. Appl. Environ. Microbiol. 65, 1061-1070 (1999).

[29] Chacon, M. R. et al. Microscopic and transcriptome analyses of early colonization of tomato roots by Trichoderma harzianum. Int. Microbiol. 10, 19-27 (2007).

[30] Salas-Marina, M. A. et al. Colonization of Arabidopsis roots by Trichoderma atroviride promotes growth and enhances systemic disease resistance through jasmonic acid/ethylene and salicylic acid pathways. Eur. J. Plant Pathol. 131, 15-26 (2011).

[31] Ho, C. L. et al. De novo transcriptome analyses of host-fungal interactions in oil palm (Elaeis guineensis Jacq.). BMC Genom. 17, 66 (2016).

[32] Alonso-Ramirez, A. et al. Salicylic acid prevents Trichoderma harzianum from entering the vascular system of roots. Mol. Plant Pathol. 15, 823-831 (2014).

[33] Martinez-Medina, A., Appels, F. V. W. & van Wees, S. C. M. Impact of salicylic acid- and jasmonic acid-regulated defences on root colonization by Trichoderma harzianum T-78. Plant Signal Behav. 12, e1345404 (2017).

[34] Ahmad, P. et al. Role of Trichoderma harzianum in mitigating NaCl stress in Indian mustard (Brassica juncea L.) through antioxidative defense system. Front. Plant Sci. 6, 868 (2015).

[35] Mastouri, F., Bjorkman, T. & Harman, G. E. Trichoderma harzianum enhances antioxidant defense of tomato seedlings and resistance to water deficit. Mol. Plant Microbe Interact. 25, 1264-1271 (2012).

[36] Vinale, F. et al. Trichoderma-plant-pathogen interactions. Soil Biol. Biochem. 40, 1-10 (2008).

[37] Moran-Diez, E. et al. Transcriptomic response of Arabidopsis thaliana after 24 h incubation with the biocontrol fungus Trichoderma harzianum. J. Plant Physiol. 169, 614-620 (2012).

[38] Suzuki, N. et al. Respiratory burst oxidases: the engines of ROS signaling. Curr. Opin. Plant Biol. 14, 691-699 (2011).

[39] Kimura, S., Waszczak, C., Hunter, K. & Wrzaczek, M. Bound by fate: the role of reactive oxygen species in receptor-like kinase signaling. Plant Cell 29, 638-654 (2017).

[40] Navazio, L. et al. Calcium-mediated perception and defense responses activated in plant cells by metabolite mixtures secreted by the biocontrol fungus Trichoderma atroviride. BMC Plant Biol. 7, 41 (2007).

[41] Contreras-Cornejo, H. A., Macias-Rodriguez, L., Beltran-Pena, E., Herrera-Estrella, A. & Lopez-Bucio, J. Trichoderma-induced plant immunity likely involves both hormonal- and camalexin-dependent mechanisms in Arabidopsis thaliana and confers resistance against necrotrophic fungi Botrytis cinerea. Plant Signal Behav. 6, 1554-1563 (2011).

[42] Asai, T. et al. MAP kinase signalling cascade in Arabidopsis innate immunity. Nature 415, 977-983 (2002).

[43] Higashi, K. et al. Modulation of defense signal transduction by flagellin-induced WRKY41 transcription factor in Arabidopsis thaliana. Mol. Genet. Genom. 279, 303-312 (2008).

[44] Shoresh, M., Yedidia, I. & Chet, I. Involvement of jasmonic acid/ethylene signaling pathway in the systemic resistance induced in cucumber by Trichoderma asperellum T203. Phytopathology 95, 76-84 (2005).

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

[46] Van der Ent, S. et al. MYB72 is required in early signaling steps of rhizobacteria-induced systemic resistance in Arabidopsis. Plant Physiol. 146, 1293-1304 (2008).

[47] Segarra, G., Van der Ent, S., Trillas, I. & Pieterse, C. M. MYB72, a node of convergence in induced systemic resistance triggered by a fungal and a bacterial beneficial microbe. Plant Biol. (Stuttg.) 11, 90-96 (2009).

[48] Ruzicka, K. et al. Ethylene regulates root growth through effects on auxin biosynthesis and transport-dependent auxin distribution. Plant Cell 19, 2197-2212 (2007).

[49] Harman, G. E., Howell, C. R., Viterbo, A., Chet, I. & Lorito, M. Trichoderma species -opportunistic, avirulent plant symbionts. Nat. Rev. Microbiol 2, 43-56 (2004).

[50] Carol, R. J. & Dolan, L. The role of reactive oxygen species in cell growth: lessons from root hairs. J. Exp. Bot. 57, 1829-1834 (2006).

[51] Gibbs, D. J. et al. AtMYB93 is a novel negative regulator of lateral root development in Arabidopsis. New Phytol. 203, 1194-1207 (2014).

[52] Segarra, G. et al. Proteome, salicylic acid, and jasmonic acid changes in cucumber plants inoculated with Trichoderma asperellum strain T34. Proteomics 7, 3943-3952 (2007).

[53] Saijo, Y. & Reimer-Michalski, E.-M. in Epigenetic Memory and Control in Plants (eds Gideon Grafi & Nir Ohad) 57-76 (Springer, Berlin, Heidelberg, 2013).

[54] Ding, B. & Wang, G. L. Chromatin versus pathogens: the function of epigenetics in plant immunity. Front. Plant Sci. 6, 675 (2015).

[55] Espinas, N. A., Saze, H. & Saijo, Y. Epigenetic control of defense signaling and priming in plants. Front. Plant Sci. 7, 1201 (2016).

[56] Satgé, C. et al. Reprogramming of DNA methylation is critical for nodule development in Medicago truncatula. Nat. Plants 2, 16166 (2016).

[57] Gutzat, R. & Mittelsten Scheid, O. Epigenetic responses to stress: triple defense?. Curr. Opin. Plant Biol. 15, 568-573 (2012).

[58] Conrath, U., Beckers, G. J., Langenbach, C. J. & Jaskiewicz, M. R. Priming for enhanced defense. Annu. Rev. Phytopathol. 53, 97-119 (2015).

[59] Gassmann, W. Alternative splicing in plant defense. Curr. Top. Microbiol. Immunol. 326, 219-233 (2008).

[60] Shang, X., Cao, Y. & Ma, L. Alternative splicing in plant genes: a means of regulating the environmental fitness of plants. Int. J. Mol. Sci. 18, 432 (2017). https://doi.org/10.3390/ijms18020432

[61] Pearson, W. R., Wood, T., Zhang, Z. & Miller, W. Comparison of DNA sequences with protein sequences. Genomics 46, 24-36 (1997).

[62] Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114-2120 (2014).

[63] Kim, D. et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14, R36 (2013).

[64] Anders, S., Pyl, P. T. & Huber, W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics 31, 166-169 (2015).

[65] Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139-140 (2010).

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

[67] Alexa, A., Rahnenfuhrer, J. & Lengauer, T. Improved scoring of functional groups from gene expression data by decorrelating GO graph structure. Bioinformatics 22, 1600-1607 (2006).

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

[69] Ramsak, Z. et al. GoMapMan: integration, consolidation and visualization of plant gene annotations within the MapMan ontology. Nucleic Acids Res. 42, D1167-D1175 (2014).

[70] Herrero, J., Valencia, A. & Dopazo, J. A hierarchical unsupervised growing neural network for clustering gene expression patterns. Bioinformatics 17, 126-136 (2001).

[71] Pfaffl, M.W. A new mathematical model for relative quantification in real-time PCR. Nucleic Acids Res. 29, e45 (2001)

[72] Scherm, B., Schmoll, M., Balmas, V., Kubicek, C. P. & Migheli, Q. Identification of potential marker genes for Trichoderma harzianum strains with high antagonistic potential against Rhizoctonia solani by a rapid subtraction hybridization approach. Curr. Genet. 55, 81-91 (2009).

[73] Mittal, V. K. & McDonald, J. F. R-SAP: a multi-threading computational pipeline for the characterization of high-throughput RNA-sequencing data. Nucleic Acids Res. 40, e67 (2012).

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

[75] Tian, T. et al. agriGOv2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Res. 45, W122-W129 (2017).

[76] Supek, F., Bosnjak, M., Skunca, N. & Smuc, T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS ONE 6, e21800 (2011).

[77] Schenk, S. T. et al. N-Acyl-homoserine lactone primes plants for cell wall reinforcement and induces resistance to bacterial pathogens via the salicylic acid/oxylipin pathway. Plant Cell 26, 2708-2723 (2014).