- Research article
- Open Access
Downy mildew resistance induced by Trichoderma harzianum T39 in susceptible grapevines partially mimics transcriptional changes of resistant genotypes
BMC Genomics volume 13, Article number: 660 (2012)
Downy mildew, caused by Plasmopara viticola, is one of the most severe diseases of grapevine and is commonly controlled by fungicide treatments. The beneficial microorganism Trichoderma harzianum T39 (T39) can induce resistance to downy mildew, although the molecular events associated with this process have not yet been elucidated in grapevine. A next generation RNA sequencing (RNA-Seq) approach was used to study global transcriptional changes associated with resistance induced by T39 in Vitis vinifera Pinot Noir leaves. The long-term aim was to develop strategies to optimize the use of this agent for downy mildew control.
More than 14.8 million paired-end reads were obtained for each biological replicate of T39-treated and control leaf samples collected before and 24 h after P. viticola inoculation. RNA-Seq analysis resulted in the identification of 7,024 differentially expressed genes, highlighting the complex transcriptional reprogramming of grapevine leaves during resistance induction and in response to pathogen inoculation. Our data show that T39 has a dual effect: it directly modulates genes related to the microbial recognition machinery, and it enhances the expression of defence-related processes after pathogen inoculation. Whereas several genes were commonly affected by P. viticola in control and T39-treated plants, opposing modulation of genes related to responses to stress and protein metabolism was found. T39-induced resistance partially inhibited some disease-related processes and specifically activated defence responses after P. viticola inoculation, causing a significant reduction of downy mildew symptoms.
The global transcriptional analysis revealed that defence processes known to be implicated in the reaction of resistant genotypes to downy mildew were partially activated by T39-induced resistance in susceptible grapevines. Genes identified in this work are an important source of markers for selecting novel resistance inducers and for the analysis of environmental conditions that might affect induced resistance mechanisms.
Plasmopara viticola (Berk. and Curt.) Berl. and de Toni is a biotrophic oomycete that causes downy mildew in grapevine. This devastating disease occurs worldwide, particularly in regions with warm and wet conditions during the growing season. P. viticola mainly infects leaves and clusters of young berries and produces oil spot lesions on the adaxial leaf surface accompanied by massive sporulation on the abaxial surface. Although downy mildew can be controlled by frequent applications of chemical fungicides, concerns about the environmental impact of pesticide overuse and the development of resistant P. viticola populations have sparked an interest in alternatives to chemical treatments.
The grapevine industry relies predominantly on Vitis vinifera, which is susceptible to downy mildew. Resistance traits have been identified in wild grapevine species, and the mechanisms of resistance to downy mildew have been characterized in resistant genotypes. Transcriptomic analysis supports the view that downy mildew resistance is mainly a post-infection phenomenon and highlights the importance of transcriptional reprogramming in both resistant and susceptible genotypes in response to P. viticola inoculation[4–8]. Transcriptional changes associated with P. viticola infection of susceptible grapevines have been related to a weak defence response and to the establishment of a compatible interaction[5, 7, 9, 10]. The response of resistant genotypes has been characterized by strong and rapid transcriptional reprogramming of processes related to defence, signal transduction, and secondary metabolism, which are either not induced or induced to a lesser extent in susceptible grapevines[4, 8, 11–14]. In particular, downy mildew resistance has been correlated with enhanced expression of genes encoding pathogenesis-related (PR) proteins and enzymes of phenylpropanoid biosynthesis, and with specific modulation of signal transduction components and markers of hypersensitive response (HR) in resistant grapevines[4, 8, 11–13].
Downy mildew symptoms can be significantly reduced in susceptible grapevines by preventive application of resistance inducers, such as chitosan, laminarins[16–18], β-aminobutyric acid (BABA)[19, 20], acibenzolar-S-methyl (BTH)[21, 22], and thiamine. Treatments with plant extracts or microbial extracts[22, 25] have also been found to increase grapevine resistance to downy mildew. The ascomycete Trichoderma harzianum strain T39 (T39) significantly reduces downy mildew symptoms by activating grapevine resistance both locally and systemically, although the molecular events responsible for resistance induction have not yet been clarified. Trichoderma spp. strains have been characterized in model systems based on their ability to induce plant resistance against pathogens[26, 27] by reprogramming the plant transcriptome[28–32]. Specific strains of beneficial microorganisms can improve plant performance by activating a plant-mediated defence mechanism known as induced systemic resistance (ISR). Through root or leaf interactions[21, 34], beneficial microorganisms are recognized by the plant, which results in a mild but effective activation of the plant immune responses in all tissues. ISR confers broad-spectrum resistance to various types of pathogens and abiotic stresses[27, 36] and is usually regulated by jasmonic acid (JA)- and ethylene (ET)-dependent signalling pathways. Rather than directly activating plant defences, beneficial microorganisms that induce resistance usually prime the plant so that it responds more strongly upon exposure to the stress condition[35, 36]. Primed plants display faster and/or stronger activation of the defence responses after pathogen inoculation. Because plant defences are activated only when they are really needed, priming provides advantages in terms of energy costs for the plant, and it is probably evolved to save energy under pathogen-free conditions. The benefits of priming outweigh its costs when disease does occur, and priming is seen as a promising strategy in modern disease management. The absence of apparent energy costs associated with T39-induced resistance in grapevine suggested a priming state activation. However, the molecular mechanisms underlying resistance induction have been only partially revealed[26, 43], and information regarding induced resistance in grapevine is particularly scarce[20, 44].
In this study, we used Illumina RNA-Seq analysis to characterize the global transcriptional dynamics associated with T39-induced resistance in grapevine. To the best of our knowledge, this study is the first to use high-throughput sequencing technology to investigate molecular events underlying induced resistance in plants, and it is also the first transcriptome-wide characterization of resistance induced by a beneficial microbe in a non-model plant. Our analysis revealed that the reduction of downy mildew symptoms is related to a complex transcriptional reprogramming in T39-treated plants, both before and 24 h after pathogen inoculation. In particular, the reaction of T39-treated plants to pathogen inoculation is associated with enhanced expression of P. viticola-responsive genes and specific modulation of some genes related to defence in resistant grapevines. Our study has identified genes that could be valuable as markers of ISR activation for subsequent selection of new resistance inducers with improved ability to stimulate plant defences.
RNA-Seq sequencing and mapping of reads to the grapevine genome
The ability of T39 to activate local and systemic resistance processes has been previously reported in grapevine[21, 42], and leaf treatment was used to study responses to T39 treatment and pathogen inoculation in the same plant organ. Leaf samples were collected before and 24 h after P. viticola inoculation from T39-treated and control plants, and resistance induction was confirmed by phenotypic analysis (Figure1). Four different treatments were analysed by RNA-Seq: control (C), T39-treated (T39), P. viticola-inoculated control (C+P.v.) and P. viticola-inoculated T39-treated (T39+P.v.) plants. Three biological replicates of each treatment were analysed, with each biological replicate sequenced twice in separate lanes (sequencing replicate). Between 6.1 and 18.8 million paired-end reads of 100 nucleotides were obtained for each sequencing replicate, and an average of 94% of these passed the quality control test (filtered reads; Additional files1 and2). Summing the reads of two sequencing replicates, more than 14.8 million filtered reads were obtained for each biological replicate (Table1), corresponding to a coverage of at least 32× the grapevine transcriptome (Table2).
An average of 91% of filtered reads mapped to the grapevine genome (mapped reads) for each biological replicate (Table1), and similar percentages mapped for each sequencing replicate (Additional file1). An average of 96% of mapped reads matched to unique locations (unique reads), 4% displayed multiple matches (multi-reads) to the grapevine genome (Table1 and Additional files1 and3), and about 77% of mapped reads matched to grapevine genes (Additional files3 and4).
Gene expression estimation by RNA-Seq
Gene expression level was assessed on the basis of both unique and multi-reads (Additional file3) to improve evaluation of members of multigene families[45, 46], and it was expressed as fragments per kilobase of transcript per million fragments mapped (FPKM) values using Cufflinks. Between one and nine transcript isoforms were recognized for each grapevine gene, and 3,679 to 7,548 novel genes were identified in each sequencing replicate (Additional file5), representing important new information for genome annotation. About 66% of grapevine genes were expressed in each sequencing replicate (Additional files5 and6). High correlations between sequencing replicates were obtained (Pearson’s correlation coefficients and R2 values greater than 0.98 and 0.95, respectively; Additional files7,8 and9), and read counts were summed to obtain better coverage and improve variance estimation of each biological replicate. Considering the expression of all grapevine genes, the Pearson’s correlation coefficient between control and T39-treated plants (T39 vs. C) was greater than that between control and P. viticola-inoculated plants (C+P.v. vs. C and T39+P.v. vs. C; Table3), suggesting that few transcriptional changes were caused by T39 treatment and that major transcriptional reprogramming occurred after pathogen inoculation.
Grapevine genes differentially expressed during Trichoderma harzianum T39-induced resistance
DeSeq statistical analysis revealed that 7,024 genes were differentially expressed with respect to the control with a false discovery rate (FDR) of 5% and a minimum fold-change of two in at least one pairwise comparison (Additional file10). About 90% of these genes were expressed in all treatments (Figure2a), indicating that T39 and P. viticola caused a significant modulation (up- or down-regulation greater than 2-fold) of genes normally expressed in control leaves. Relatively few genes were specifically expressed in response to the T39 treatment (12 genes) and after P. viticola inoculation in control (36 genes) and T39-treated plants (22 genes).
The expression profiles of the differentially expressed genes were determined by cluster analysis, and genes were grouped into 10 clusters based on their expression modulation (Figure2b and Additional file10). Among genes modulated by T39, 343 had comparable expression levels before and after P. viticola inoculation in T39-treated plants (cluster 1), while 233 genes were initially modulated by T39 and then further modulated by P. viticola inoculation in T39-treated plants (cluster 2). In addition, the expression of 479 genes was affected exclusively by T39 treatment and not by P. viticola (cluster 3). P. viticola inoculation resulted in modulation of 3,454 genes with comparable expression levels in control and T39-treated plants at 24 h after inoculation; of these, 948 genes were directly modulated by T39 treatment (cluster 4) and 2,506 genes were not (cluster 5). Interestingly, 868 genes showed reinforced modulation in T39-treated plants compared to control plants after P. viticola inoculation (ISR-primed genes), indicating enhancement of the grapevine defence reaction to the pathogen. Among these genes, 518 and 350 genes were modulated (cluster 6) or not (cluster 7) by T39 treatment, respectively. A total of 888 genes were modulated by P. viticola exclusively in control plants; of these, 267 genes were directly modulated by T39 treatment (cluster 8) and 621 genes were not (cluster 9). Conversely, 759 genes were modulated by P. viticola exclusively in T39-treated plants (cluster 10, ISR-responsive specific genes), and they represent the specific reaction to P. viticola of T39-treated plants. ISR-responsive specific genes were mainly induced (63%), while genes modulated by P. viticola exclusively in control plants were mainly repressed (55%).
Functional annotation of differentially expressed genes
Differentially expressed genes were automatically annotated and then grouped into 15 functional categories of Gene Ontology (GO) biological processes (Figure3); genes that could not be associated with any biological process category were assigned to the GO root (biological process). Although the group of genes with unknown function was the largest (35%), 11 categories were represented to different degrees in differentially expressed genes compared to the grapevine transcriptome. DNA (5.0%) and protein (8.4%) metabolic processes were significantly underrepresented among the differentially expressed genes, but response to stress (5.9%), response to stimulus (9.8%), and signal transduction (2.7%) categories were significantly overrepresented. Likewise, transport (6.2%) and metabolism of carbohydrates (3.4%), lipids (2.9%), secondary compounds (1.1%), and energy (1.2%) were significantly overrepresented. The large fraction of differentially expressed genes with unknown function highlighted the identification of novel processes potentially relevant for the induction of plant resistance and for the response of grapevine to P. viticola inoculation. However, additional studies are required to better characterize the relevance of these genes in the mechanisms of T39-indiuced resistance.
After genes with unknown function were discounted, a specific distribution of up- and down-regulated functional categories was observed in the various clusters (Figure4). Genes of cluster 1 involved in carbohydrate, protein, and secondary metabolism were mainly induced by T39 (69, 76 and 83%, respectively), whereas those related to lipid metabolism and transport were mainly repressed (86 and 65%, respectively; Figure4a). Genes of cluster 2 related to protein metabolism and transcription were mainly repressed by T39 and subsequently induced by P. viticola inoculation of T39-treated plants. More specifically, T39 treatment directly induced the expression of genes encoding protein kinases and receptor-like protein kinases, MYB transcription factors, hormone-responsive genes, defence-related proteins, and enzymes of secondary metabolism. Annotation of genes modulated exclusively by T39 (cluster 3) revealed down-regulation of processes related to response to stimulus, cellular component organization, and protein metabolism (Figure4b).
Genes that were similarly modulated by P. viticola in control and T39-treated plants (clusters 4 and 5) were mainly related to primary metabolic processes and signal transduction (Figure4c). In particular, genes involved in cellular component organization and carbohydrate metabolism were mainly induced (74 and 65%, respectively). More specifically, P. viticola induced six cellulose synthase-like genes, six expansin genes, as well as genes encoding enzymes of glucan biosynthesis and hexose metabolism. Conversely, lipid metabolic process, signal transduction, response to stress, and response to stimulus were mainly repressed by P. viticola (60, 58, 56 and 55%, respectively); this included genes encoding protein kinases, receptor-like protein kinases, disease resistance proteins, and PRs.
The categories of response to stress, secondary metabolism, and DNA metabolism were significantly overrepresented in clusters 6 and 7, and the genes in these categories were mainly up-regulated following P. viticola inoculation (Figure4d). In particular, expression of protein kinase, disease resistance protein, chitinase (Chit), stilbene synthase (STS), resveratrol o-methyltransferase, phenylalanine ammonia-lyase (PAL), and transcription factor genes was enhanced after P. viticola inoculation in T39-treated plants compared with control plants.
Clusters of genes modulated by P. viticola exclusively in control plants revealed global down-regulation of genes related to protein metabolic process (62 and 71% in clusters 8 and 9, respectively), response to stress (70 and 60%, respectively), and response to stimulus (66 and 57%, respectively; Figure4e). Interestingly, the same categories repressed by P. viticola in control plants were globally induced in T39-treated plants (cluster 10; Figure4f). In particular, genes related to protein metabolic process, response to stress, response to stimulus, and transcription were mainly up-regulated by P. viticola in T39-treated plants (77, 67, 61, and 67%, respectively). This active reaction to the pathogen in T39-treated plants included expression of protein kinases, transcription factors, auxin and JA/ET signals, phenylpropanoid biosynthesis, and defence-related genes.
Validation of RNA-Seq analysis
To validate the RNA-Seq results, we used real-time RT-PCR to analyse the expression levels of 24 differentially expressed genes (Figure5). We selected genes with different expression profiles and expression levels, including genes associated with different functional categories (Additional file11). Although the extent of modulation revealed by real-time RT-PCR and RNA-Seq may differ[6, 48, 49], the real-time RT-PCR expression profiles in our analysis were in complete agreement with the RNA-Seq data for 19 genes. The expression profiles generated by real-time RT-PCR and RNA-Seq differed for three genes in one treatment and for two genes in all treatments. These differences could be due to differences in sensitivity, particularly in distinguishing members of multigene families. Specificity for RNA-Seq analysis was observed in two grapevine chitinases (Table4) that are known to be strongly (Chit1a) or weakly (Chit1b) up-regulated after P. viticola inoculation.
We also validated the RNA-Seq method by confirming the expression of previously identified markers of plant response to Trichoderma spp. and P. viticola (Table4). Our RNA-Seq data confirmed P. viticola-dependent up-regulation of PR 1, PR 2, PR 4, and PR 10 genes in susceptible grapevines[5, 7, 11, 19, 51, 52]. Expression profiles of PR 2, PR 4, PR 5, PR 10, and lipoxygenase 9 (LOX 9) were consistent with previous real-time RT-PCR data obtained during T39-induced resistance. The RNA-Seq analysis revealed T39-dependent induction of genes that are known to be involved in Trichoderma-induced resistance in other systems[26, 29, 53, 54]. Moreover, RNA-Seq analysis confirmed the modulation of genes known to be affected by P. viticola in grapevines[4–7, 9, 11, 17, 52].
RNA-Seq analyses of leaf samples collected before and 24 h after P. viticola inoculation from T39-treated and control plants resulted in the identification of 7,024 differentially expressed genes. These genes formed 10 clusters of different expression profiles, highlighting the complex transcriptional reprogramming of T39-induced resistance. T39 treatment directly affected the expression of grapevine genes and to a greater extent enhanced grapevine response to P. viticola inoculation, indicating a dual effect of T39. At the sampling time points selected in this analysis, we showed that a limited number of changes in gene expression were caused by T39 treatment and that more intense transcriptional reprogramming took place after pathogen inoculation. In particular, opposing modulation of genes related to response to stress and protein metabolism was observed in T39-treated and control plants at 24 h after P. viticola inoculation, indicating inhibition of disease-related processes and induction of active defence in T39-treated grapevines. Based on the expression profiles, genes directly modulated by T39, as well as genes with reinforced or specific modulation in T39-treated plants after pathogen inoculation are strong candidates for activation of plant self-protection and consequent inhibition of disease-related processes and symptoms development.
Grapevine processes directly affected by Trichoderma harzianum T39
Analysis of the expression profiles revealed a set of genes directly modulated by T39 (182 up- and 161 down-regulated) and whose expression was not affected by the subsequent pathogen inoculation (cluster 1). Several genes of cluster 1 were involved in signal transduction processes, indicating that they may be related to the initial events of recognition of the beneficial microorganism and induction of resistance, as recently demonstrated by the proteomic analysis of T39-induced resistance. Enzymes known to mediate microbial recognition and to trigger defence responses in plant species were up-regulated by T39, and these included 66 receptor-like protein kinase genes, three protein kinase genes, and one protein phosphatase gene. In particular, the serine/threonine kinase receptor (STKR) was similar to the Arabidopsis FRK1 induced by bacterial flagellins, and a protein kinase was homologous to the Trichoderma-induced kinase (TIPK) of cucumber. Furthermore, the STKR gene (JG391826) is induced by P. viticola in resistant grapevines, suggesting that it plays a crucial role in the activation of specific defence processes against downy mildew.
Our data suggest that resistance induction was also mediated by hormone signalling and transcriptional reprogramming. Genes related to ET metabolism (two 1-aminocyclopropane-1-carboxylate oxidases) were induced by T39 treatment, together with two MYB (including MYB72), one FAMA, and one NAC transcription factors. Activation of ET metabolism was in agreement with previous data showing that JA/ET signals are involved in T39-induced resistance[42, 58]. In Arabidopsis, MYB72 plays a crucial role in Trichoderma-induced resistance and is required in early signalling of rhizobacteria-induced resistance in Arabidopsis. The NAC gene is involved during oxidative stress, and response to oxidative stress after T39 treatment was suggested by modulation of two peroxidases, three glutathione S-transferases (GSTs), and one thioredoxin. Moreover, defence-related genes (two STSs, a Chit, and a cellulose synthase) were also directly induced, indicating that the maximum expression level of these genes was probably reached after T39 treatment and that this was sufficient to contribute to defence against subsequent P. viticola inoculation. Other defence-related genes were pre-induced by T39 and further induced after P. viticola inoculation (cluster 2), indicating reinforcement of grapevine defence after pathogen inoculation.
Enhanced reaction of T39-treated plants to Plasmopara viticola
Plant resistance induced by beneficial microorganisms has been associated with faster and/or stronger activation of defence responses after pathogen inoculation. P. vitiocla-responsive genes with enhanced expression in T39-treated plants were clustered to distinguish between those directly modulated (cluster 6) and those not modulated (cluster 7) by T39 treatment. These expression profiles provide strong support for the view that Trichoderma spp. may have a dual effect: it directly stimulates induction of some genes and further reinforces modulation of these and other genes after pathogen inoculation. The dual effect was also reported for defence-gene modulation during resistance induced by sulfated laminarin and for phytoalexin accumulation during resistance induced by Rheum palmatum extracts.
Enhanced reaction of T39-treated plants included induction of genes related to response to stimulus and response to stress categories, suggesting improvement of signalling pathways and activation of defence reactions in response to pathogen inoculation. In particular, 50 receptor-like kinase genes (35 in cluster 6 and 15 in cluster 7) and six protein kinase genes (one in cluster 6 and five in cluster 7) of signal transduction were primed in T39-treated plants. In agreement with our results, priming was associated with increased expression of mitogen-activated protein kinases (MPK3 and MPK6) in Arabidopsis plants. Up-regulation of MPK genes was also associated with the reaction of resistant grapevines to P. viticola[4, 6], and post-translational modification may be the additional mechanism for activating MPK s in response to pathogens. Moreover, expression of stress-related genes was enhanced in inoculated T39-treated plants, and these genes included PR 2, PR 4, PR 5, PR 10, Chit3, and osmotin (OSM1). Similar profiles were observed for Arabidopsis PR genes during T. asperelloides T203 and Pseudomonas syringae interactions, indicating common mechanisms of Trichoderma-induced resistance against biotrophs. Primed profiles were found in resistance-related genes (12 genes in cluster 6 and eight in cluster 7). Peronospora parasitica resistance genes (RPP) must be expressed at optimal levels to function against downy mildew in Arabidopsis, and specific profiles of grapevine RPP genes suggest fine tuning during T39-induced resistance.
A key role of transcriptional regulation and secondary metabolic processes in T39-induced resistance was also indicated by the priming profiles of four MYB, two WRKY, and 12 ET-responsive transcription factors (e.g., AP2), and of PAL and STSs genes. Moreover, the category DNA metabolic process was significantly overrepresented in ISR-primed genes. In particular, expression of the histone-lysine N-methyltransferases SUVR4 and ATX2 was enhanced after P. viticola inoculation, reflecting possible involvement of epigenetic modifications[62–64] in T39-induced resistance.
Genes with enhanced expression in T39-treated plants have been previously related to defence against downy mildew in resistant grapevines. This was the case for PR genes (PR 2, PR 4, PR 10, and OSM1), invertase genes (INV1 and INV2), and genes related to secondary metabolism (PAL, STS s, STKR, copper-containing amine oxidase, polyphenol oxidase, and three resveratrol O-methyltransferases), for which up-regulation was greater in resistant than in susceptible grapevines upon P. viticola inoculation[4, 8, 11–14]. The expression profiles of these marker genes suggest that the defence processes usually activated against downy mildew in resistant grapevines are partially stimulated in susceptible plants by T39-induced resistance.
Disease-related processes inhibited in T39-treated plants
Another important aspect of T39-induced resistance was evidenced by genes modulated by P. viticola exclusively in control plants (clusters 8 and 9). These genes were mainly down-regulated (57 and 55%, respectively), and they reflect exploitation of cellular resources and/or suppression of defence responses during the compatible interaction. Down-regulation caused by P. viticola in control plants involved categories of response to stress and primary metabolic processes. Many signalling components (kinase, phosphatase, calmodulin, and calcium signalling), transcription factors (WRKY and MYB), and disease resistance proteins were repressed, supporting the view that suppression of endogenous signals is required to establish the compatible interaction. Three ABC transporters (e.g., ABC C) were also repressed, and suppression of some ABC transporters increases the susceptibility to oomycete pathogens. Specific alteration of carbohydrate metabolism by P. viticola in control plants was highlighted by modulation of glucosidase, galactosidase, mannosidase, and sucrose synthase (SS) genes. Moreover, repressed genes of cluster 9 were classified into the categories of energy metabolism (a phytochrome C, two malic enzymes, and a ribokinase) and defence response (five chitinases, three glucanases, three superoxide dismutases, and a callose synthase), reflecting disease-related process employed by P. viticola only in control plants. P. viticola might need to actively suppress plant defences during leaf colonisation through microbial effectors, as demonstrated in other oomycetes[66, 67]. Particularly, hemibiotrophic and biotrophic species establish intimate associations with plants[66, 67]. To establish infection, these pathogens must suppress the plant defence and manipulate the host metabolism by microbial effectors (virulence factors) that are translocated inside the plant cell or secreted into the extracellular space within plant tissue[66, 68]. P. viticola effector genes have been recently described, and modulation of host cell defences through virulence factors in susceptible grapevines has been indicated by histochemical and transcriptomic analyses. In agreement with the phenotypic observations, the specific modulation of grapevine genes in inoculated control but not in inoculated T39-treated plants indicates that T39-induced resistance acts by inhibiting some disease-related processes and probably by interfering with some pathogen-induced processes.
Specific transcriptional response of T39-treated plants to Plasmopara viticola
In contrast to broad down-regulation of genes in control plants (clusters 8 and 9), genes specifically modulated by P. viticola in T39-treated plants (cluster 10) were mainly induced (63%). These opposing reactions to P. viticola are particularly evident in genes related to protein metabolism, response to stimulus, and response to stress, which were mainly induced in T39-treated plants and mainly repressed in control plants. Up-regulation of genes associated with these categories has been observed in resistant grapevines[4, 13, 71], indicating that the defence processes of resistant genotypes could be partially activated in susceptible varieties by T39-induced resistance. In particular, all NBS-encoding resistance (NBS R) genes modulated by P. viticola in control plants (cluster 9) were repressed, whereas those with ISR-responsive specific profiles (cluster 10) or ISR-primed profiles (cluster 7) were mainly up-regulated. Opposing modulation of NBS R genes probably reflects suppression of plant defence in control plants and activation of defence responses in T39-treated plants. Interestingly, the NBS R genes of clusters 7 and 10 belong mainly to the Va component genome of grapevine, indicating subgenome-dependent regulation of gene expression in grapevine.
Defence signals specifically activated in T39-treated plants included those mediated by auxin (two auxin transporters, two auxin-induced proteins, and two indole-3-acetic acid amido synthetases), ET (ACC oxidase and five ERF transcription factors), and JA (three lipoxygenases and two fatty acid desaturases). The role of JA/ET signalling pathways in T39-induced resistance has also been demonstrated by phytohormone-affected Arabidopsis mutants and by expression analysis of grapevine marker genes. ISR is commonly regulated by JA/ET-dependent signalling pathways, and it is especially active against pathogens deterred by defences that are controlled by JA and ET[33, 36]. The auxin response pathway is connected to the SA and JA/ET signalling networks, and crosstalk between hormonal pathways enables the fine tuning of defence mechanisms so that the plant can tailor its response to the specific invader. ET exerts its resistance-stimulating activity in concert with JA, and JA pathways are involved in the reaction to P. viticola in resistant grapevines[4, 6]. Thus, enhancement of JA/ET signals in T39-treated plants supports the view that increased resistance to downy mildew is mediated by partial activation of extant defence mechanisms normally activated in resistant genotypes. Reaction to the pathogen was also mediated by specific up-regulation of 59 receptor kinases, 10 protein kinases, two bHLHs genes, one MYB gene, and the NPR1.1 gene.
Our results also suggest that the cell redox balance is altered in T39-treated plants after pathogen inoculation. The reaction of T39-treated plants to pathogen inoculation included the induction four peroxidases and a GST. Antioxidant enzymes are often induced in response to pathogens, and alteration of oxidative-stress metabolism has a prominent role in the T39-induced resistance of grapevine to downy mildew. Peroxidases play several important roles in pathogen resistance by contributing to the production of reactive oxygen species, the reinforcement of cell walls, and the production of phytoalexins. Accumulation of stilbene phytoalexin is one of the most important defence processes activated by resistant grapevines in response to P. viticola; genes of phenylpropanoid biosynthesis (flavanone 3-dioxygenase, laccase, and dihydroflavonol-4-reductase) were specifically induced in T39-treated plants, confirming activation of pathways known in resistant genotypes. However, additional defence mechanisms against P. viticola are activated in resistant genotypes. In particular, HR-related genes (Avr/Cf9, Hin1, and Hsr203j) were not induced and localised HR necrosis was not observed in T39-treated plants.
Common transcriptional response of control and T39-treated plants to Plasmopara viticola
Although specific transcriptional reprogramming of T39-treated plants was observed, 3,454 genes had comparable expression levels in control and T39-treated plants after P. viticola inoculation (clusters 4 and 5). The pathogen-responsive processes not affected by resistance induction were mainly related to primary metabolism and signal transduction. In particular, expression profiles of genes related to starch metabolism (up-regulation of α-amylase and sugar transporters, down-regulation of β-amylase, glucose-1-phosphate adenylyltransferases, sucrose phosphate synthase, and phosphoglucan water dikinase) indicated the source-to-sink transition of P. viticola-infected leaves. Likewise, genes related to photosynthesis (two quinone oxidoreductases, a chlorophyllide oxygenase, a protochlorophyllide transporter, a chlorophyllase-2 and chlorophyll a-b binding proteins) and to the Calvin cycle (two rubisco genes and a phosphoglycerate kinase gene) were similarly modulated in inoculated control and T39-treated plants, possibly reflecting the establishment of a compatible interaction. P. viticola inoculation also resulted in the down-regulation of genes involved in the signal transduction processes (20 receptor-like protein kinase genes, 11 protein kinase genes, and five protein phosphatase genes) and defence response (23 probable disease resistance genes, six unspecified PR genes, a GST gene and a thaumatin-like gene), indicating a pathogen-dependent suppression of the host reaction mechanisms.
Other markers of P. viticola infection in susceptible grapevines showed comparable modulation in T39-treated and control plants, such as enzymatic resistance protein, zeaxanthin epoxidase, catalase, isoforms of chalcone synthase, and dihydroflavonol 4-reductase. These results suggested that transcriptional changes associated with the compatible interaction[5, 7, 10, 13] were not completely inhibited in T39-treated plants, which is consistent with the observation that downy mildew symptoms were reduced but not completely blocked by T39 treatment.
The transcriptome analysis reported here represents a major contribution to the characterization of induced resistance mechanisms in a non-model plant. We used the RNA-Seq approach to characterize the transcriptional changes associated with the reduction of downy mildew symptoms in T39-treated grapevines and thereby established a foundation for a more detailed time-course analysis of induced resistance in grapevines. Our data show that T39 directly activates microbial recognition mechanisms in absence of pathogen infection and enhances the expression of defence processes after P. viticola inoculation. Reduction of downy mildew symptoms is related to inhibition of disease-related processes and activation of defence mechanisms after P. viticola inoculation of T39-treated plants. T39-induced resistance is associated with enhanced expression of specific genes related to resistance against downy mildew in wild grapevines, indicating that induced resistance can partially mimic defence processes of resistant genotypes. In addition, the genes identified in this work represent an important source of markers of resistance induction and can be used to select beneficial microorganisms with an improved ability to induce plant resistance. These markers can also be used to clarify how environmental conditions affect induced resistance.
Two-year-old plants of the susceptible grapevine (V. vinifera) cultivar Pinot Noir grafted onto Kober 5BB were individually planted in 2.5 L pots containing a mixture of peat and pumice (3:1). Plants were grown for two months in a greenhouse at 25 ± 1°C with a photoperiod of 16 h light and relative humidity (RH) of 70 ± 10%. T. harzianum T39 (Trichodex, Makhteshim Ltd., Israel) was applied to grapevine leaves at 8 g L-1 in water, corresponding to a conidial suspension of 105 colony-forming units (cfu) mL-1. A P. viticola isolate was collected from an untreated vineyard in the Trentino region (northern Italy) and maintained by subsequent inoculations on V. vinifera Pinot Noir plants under controlled greenhouse conditions. Plants with oil spot symptoms were incubated overnight in the dark at 99–100% RH. P. viticola sporangia were collected by washing the abaxial surfaces bearing freshly sporulating lesions with cold distilled water, and the inoculum concentration was adjusted to 105 sporangia mL-1 with a haemocytometer.
Bioassay of induced resistance in grapevine
All leaves of each shoot were treated three times with T39 (at one, two and three days before pathogen inoculation) to induce the greatest phenotypic resistance response. Leaves directly treated with the resistance inducer were analysed to study local effects of T39-induced resistance, and untreated plants were used as control. The suspensions were applied to the abaxial and adaxial leaf surfaces using a compressed air hand sprayer (20–30 mL for each plant).
One day after the last T39 treatment, the abaxial surfaces of all leaves of each plant were inoculated with a P. viticola suspension (105 sporangia mL-1) using a compressed-air hand sprayer (20–30 mL for each plant). All plants were then incubated overnight in the dark at 25 ± 1°C with 99–100% RH and then kept under controlled greenhouse conditions.
Six days after inoculation with the pathogen, all plants were incubated overnight in the dark at 25°C with 99–100% RH. Disease severity was visually assessed as the percentage of abaxial leaf area covered by sporulation. Each treatment was carried out on six replicates (plants). An ANOVA analysis was performed using the Statistica 9 software (StatSoft, Tulsa, OK) followed by Tukey’s test to detect significant differences (P < 0.05) in disease severity between treatments.
Sample collection and RNA isolation
Leaf samples were collected from control and T39-treated plants immediately before and 24 h after P. viticola inoculation. This time point was chosen because it is associated with leaf colonization by primary hyphae[24, 55, 75] and with modulation of defence-related genes[4, 17, 19, 42] for the establishment of defence responses in resistant genotypes. Four treatments were analysed by RNA-Seq: control (C), T39-treated (T39), P. viticola-inoculated control (C+P.v.), and P. viticola-inoculated T39-treated (T39+P.v.) plants. For each treatment, leaf samples from three replicates (plants) were collected at each time point, i.e., three plants per treatment were sampled before inoculation and three plants per treatment were sampled 24 h after inoculation; different plants were sampled at each time point to avoid the effects of wounding stress. Each sample comprised three leaves taken from the same plant, and only leaves of the 4th-5th node from the top of the shoot were collected to avoid ontogenic resistance effects. Samples were immediately frozen in liquid N2 and stored at −80°C. Total RNA was extracted using the Spectrum Plant total RNA kit (Sigma-Aldrich, St. Louis, MO) and quantified using the Nanodrop 8000 (Thermo Fisher Scientific, Wilmington, DE). The quality of the RNA extracts was checked using the Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA).
Library construction and Illumina sequencing
Three biological replicates for each treatment were subjected to RNA-Seq library construction, and each library was sequenced twice in separate lanes (sequencing replicate). In a preliminary experiment, a technical replicate of one sample was run, and the results showed that the library preparation had a high level of technical reproducibility (R2 = 0.93, Pearson’s correlation = 0.96). Libraries were prepared using the TruSeq SBS v5 protocol (Illumina, San Diego, CA), and paired-end reads of 100 nucleotides were obtained using an Illumina HiSeq 2000 at Fasteris (Plan-les-Ouates, Switzerland). Briefly, poly (A)-containing mRNA was isolated from 10 μg of total RNA in two rounds of purification using poly-T oligo-attached magnetic beads. Purified mRNA was fragmented using Zn-catalysed hydrolysis and converted into double-stranded cDNA by random priming. Following end repair, indexed adapters were ligated and cDNA fragments of 200 ± 25 bp were purified. Purified cDNA was amplified by PCR and validated by Sanger sequencing, after which mRNA-Seq libraries were multiplexed (six libraries per lane) and sequenced according to the manufacturer’s instructions. The sequences have been deposited at the Sequence Read Archive of the National Center for Biotechnology under BioProject number PRJNA168987.
Mapping of sequenced reads and assessment of gene expression
For sequence quality control, in-house python scripts were used for quality trimming and sequence filtering: sequencing adapters, k-mers, and bases with a Phred quality score lower than 30 were removed from the read ends, and reads shorter than 50 nucleotides were then discarded. Filtered reads were mapped to the Pinot Noir grapevine genome Release 3 using TopHat 1.2.0 release with default settings; up to 100 hits for each read were allowed to improve gene expression evaluation of members of multigene families. Gene expression values were determined using Cufflinks 1.0.3 release, and the FPKM values were calculated for each transcript. Cufflinks default settings were adopted and gene FPKM values were computed by summing the FPKM values of different transcripts of the same gene. The option for fragment bias correction was applied, and the RABT assembly option was used to improve the identification of novel genes and transcripts. Sequencing replicates were first analysed separately, and then the read counts of two sequencing replicates were summed and FPKM values recalculated for each biological replicate.
Differential gene expression analysis
The number of reads falling into each transcript was estimated from the FPKM values according to the formula reported by Zenoni et al., and gene counts were computed by summing counts of different transcripts of the same gene. Differentially expressed genes were identified by the DESeq package, with an FDR of the Benjamini-Hochberg multiple tests of 5% (P < 0.05) based on read counts and a minimum fold-change of two in at least one pairwise comparison: T39 versus C, C+P.v. versus C, or T39+P.v. versus C. Fold-changes were calculated from FPKM expression values with the minimum expression value (10-4 FPKM) imposed on non-expressed genes. Genes modulated by pathogen inoculation exclusively in T39-treated plants were classified as ISR-responsive specific genes. Among genes modulated by P. viticola in control and T39-treated plants, genes that showed an up- or down-regulation greater than 1.5-fold in inoculated T39-treated versus inoculated control plants were selected as ISR-primed genes. This priming cut-off value is based on previous expression data of Arabidopsis genes that showed a primed expression profile after pathogen attack. ISR-primed genes were further distinguished as genes directly modulated (cluster 6) or not (cluster 7) by T39 treatment. Principal component analysis (PCA) was performed with Statistica 9 software (StatSoft, Tulsa, OK) on FPKM expression values of all grapevine genes, and the expression profiles of differentially expressed genes were determined by cluster analysis using the T-MeV 4.8.1 software.
The grapevine genes were automatically annotated using the ARGOT2 function prediction tool[83, 84], which was able to annotate grapevine protein sequences by blastp search against the UniProtKB database (downloaded from January 2012) and by association of Gene Ontology biological process terms weighted according to both semantic similarity relations and associated scores[83, 84]. Annotation of differentially expressed genes and their distribution over 15 selected functional categories was compared with the entire grapevine transcriptome using GOstat statistical analysis. For manual annotation, grapevine proteins were aligned against the Swiss Prot database (downloaded from February 2012) using blastp, and the best hits (E-value lower than 1E-5) were used to select protein descriptions. The corresponding genes of the grapevine PN40024 genotype were identified by blastp search against the 12× release with an imposed E-value lower than 1E-20.
Gene expression analysis by quantitative real-time RT-PCR
Total RNA was treated with DNase I (Invitrogen, Carlsbad, CA), and the first-strand cDNA was synthesized from 1.0 μg of total RNA using Superscript III (Invitrogen) and oligo-dT. Reactions were carried out with Platinum SYBR Green qPCR SuperMix-UDG (Invitrogen) and specific primers (Additional file11) using the Light Cycler 480 (Roche Diagnostics, Germany). The PCR conditions were: 50°C for 2 min and 95°C for 2 min as initial steps, followed by 50 cycles of 95°C for 15 s and 60°C for 1 min. Each sample was examined in three technical replicates, and dissociation curves were analysed to verify the specificity of each amplification reaction. Light Cycler 480 SV1.5.0 software (Roche) was used to extract Ct values using the second derivative calculation, and LinReg software was used to calculate reaction efficiencies. The relative expression of each gene was then calculated according to the Pfaffl equation using control plants as the calibrator. The actin gene (glimmer.VV78X114914.6_2) was used as the constitutive gene for normalization[4, 52] because its expression was not significantly affected by the treatments (Additional file12), and comparable results were obtained with VATP16. Mean expression and standard error of three biological replicates were calculated for each treatment.
Gessler C, Pertot I, Perazzolli M: Plasmopara viticola, the causal agent of downy mildew of grapes. Phytopathol Mediterr. 2011, 50: 3-44.
Pimentel D, Acquay H, Biltonen M, Rice P, Silva M, Nelson J, Lipner V, Giordano S, Horowitz A, D'Amore M: Environmental and economic costs of pesticide use. Bioscience. 1992, 42: 750-760. 10.2307/1311994.
Chen WJ, Delmotte F, Richard-Cervera S, Douence L, Greif C, Corio-Costet MF: At least two origins of fungicide resistance in grapevine downy mildew populations. Appl Environ Microbiol. 2007, 73: 5162-5172. 10.1128/AEM.00507-07.
Polesani M, Bortesi L, Ferrarini A, Zamboni A, Fasoli M, Zadra C, Lovato A, Pezzotti M, Delledonne M, Polverari A: General and species-specific transcriptional responses to downy mildew infection in a susceptible (Vitis vinifera) and a resistant (V. riparia) grapevine species. BMC Genomics. 2010, 11: 117-10.1186/1471-2164-11-117.
Polesani M, Desario F, Ferrarini A, Zamboni A, Pezzotti M, Kortekamp A, Polverari A: cDNA-AFLP analysis of plant and pathogen genes expressed in grapevine infected with Plasmopara viticola. BMC Genomics. 2008, 9: 142-10.1186/1471-2164-9-142.
Wu J, Zhang Y, Zhang H, Huang H, Folta KM, Lu J: Whole genome wide expression profiles of Vitis amurensis grape responding to downy mildew by using Solexa sequencing technology. BMC Plant Biol. 2010, 10: 234-10.1186/1471-2229-10-234.
Legay G, Marouf E, Berger D, Neuhaus JM, Mauch-Mani B, Slaughter A: Identification of genes expressed during the compatible interaction of grapevine with Plasmopara viticola; through suppression subtractive hybridization (SSH). Eur J Plant Pathol. 2011, 129: 281-301. 10.1007/s10658-010-9676-z.
Malacarne G, Vrhovsek U, Zulini L, Cestaro A, Stefanini M, Mattivi F, Delledonne M, Velasco R, Moser C: Grapevine resistance to Plasmopara viticola is associated to stilbenoids accumulation and to specific transcriptional responses as revealed by metabolic and gene expression profiling of resistant and susceptible individuals in a segregating population. BMC Plant Biol. 2011, 11: 114-10.1186/1471-2229-11-114.
Hayes MA, Feechan A, Dry IB: Involvement of abscisic acid in the coordinated regulation of a stress-inducible hexose transporter (VvHT5) and a cell wall invertase in grapevine in response to biotrophic fungal infection. Plant Physiol. 2010, 153: 211-221. 10.1104/pp.110.154765.
Gamm M, Heloir MC, Bligny R, Vaillant-Gaveau N, Trouvelot S, Alcaraz G, Frettinger P, Clement C, Pugin A, Wendehenne D, Adrian M: Changes in carbohydrate metabolism in Plasmopara viticola infected grapevine leaves. Mol Plant Microbe Interact. 2011, 24: 1061-1073. 10.1094/MPMI-02-11-0040.
Kortekamp A: Expression analysis of defence-related genes in grapevine leaves after inoculation with a host and a non-host pathogen. Plant Physiol Biochem. 2006, 44: 58-67. 10.1016/j.plaphy.2006.01.008.
Casagrande K, Falginella L, Castellarin SD, Testolin R, Di Gaspero G: Defence responses in Rpv3-dependent resistance to grapevine downy mildew. Planta. 2011, 234: 1097-1109. 10.1007/s00425-011-1461-5.
Figueiredo A, Monteiro F, Fortes AM, Bonow-Rex M, Zyprian E, Sousa L, Pais MS: Cultivar-specific kinetics of gene induction during downy mildew early infection in grapevine. Funct Integr Genomics. 2012, 12: 379-386. 10.1007/s10142-012-0261-8.
He M, Xu Y, Cao J, Zhu Z, Jiao Y, Wang Y, Guan X, Yang Y, Xu W, Fu Z: Subcellular localization and functional analyses of a PR10 protein gene from Vitis pseudoreticulata in response to Plasmopara viticola infection. Protoplasma. 2012, 10.1007/s00709-012-0384-8
Aziz A, Trotel-Aziz P, Dhuicq L, Jeandet P, Couderchet M, Vernet G: Chitosan oligomers and copper sulfate induce grapevine defense reactions and resistance to gray mold and downy mildew. Phytopathology. 2006, 96: 1188-1194. 10.1094/PHYTO-96-1188.
Aziz A, Poinssot B, Daire X, Adrian M, Bezier A, Lambert B, Joubert JM, Pugin A: Laminarin elicits defense responses in grapevine and induces protection against Botrytis cinerea and Plasmopara viticola. Mol Plant Microbe Interact. 2003, 16: 1118-1128. 10.1094/MPMI.2003.16.12.1118.
Trouvelot S, Varnier AL, Allègre M, Mercier L, Baillieul F, Arnould C, Gianinazzi-Pearson V, Klarzynski O, Joubert JM, Pugin A, Daire X: A beta-1,3 glucan sulfate induces resistance in grapevine against Plasmopara viticola through priming of defense responses, including HR-like cell death. Mol Plant Microbe Interact. 2008, 21: 232-243. 10.1094/MPMI-21-2-0232.
Allègre M, Héloir MC, Trouvelot S, Daire X, Pugin A, Wendehenne D, Adrian M: Are grapevine stomata involved in the elicitor-induced protection against downy mildew?. Mol Plant Microbe Interact. 2009, 22: 977-986. 10.1094/MPMI-22-8-0977.
Hamiduzzaman MM, Jakab G, Barnavon L, Neuhaus JM, Mauch-Mani B: Beta-aminobutyric acid-induced resistance against downy mildew in grapevine acts through the potentiation of callose formation and jasmonic acid signaling. Mol Plant Microbe Interact. 2005, 18: 819-829. 10.1094/MPMI-18-0819.
Dubreuil-Maurizi C, Trouvelot S, Frettinger P, Pugin A, Wendehenne D, Poinssot B: Beta-aminobutyric acid primes an NADPH oxidase-dependent reactive oxygen species production during grapevine-triggered immunity. Mol Plant Microbe Interact. 2010, 23: 1012-1021. 10.1094/MPMI-23-8-1012.
Perazzolli M, Dagostin S, Ferrari A, Elad Y, Pertot I: Induction of systemic resistance against Plasmopara viticola in grapevine by Trichoderma harzianum T39 and benzothiadiazole. Biol Control. 2008, 47: 228-234. 10.1016/j.biocontrol.2008.08.008.
Harm A, Kassemeyer H-H, Seibicke T, Regner F: Evaluation of chemical and natural resistance inducers against downy mildew (Plasmopara viticola) in grapevine. Am J Enol Vitic. 2011, 62: 184-192. 10.5344/ajev.2011.09054.
Boubakri H, Wahab MA, Chong J, Bertsch C, Mliki A, Soustre-Gacougnolle I: Thiamine induced resistance to Plasmopara viticola in grapevine and elicited host-defense responses, including HR like-cell death. Plant Physiol Biochem. 2012, 57: 120-133.
Godard S, Slacanin I, Viret O, Gindro K: Induction of defence mechanisms in grapevine leaves by emodin- and anthraquinone-rich plant extracts and their conferred resistance to downy mildew. Plant Physiol Biochem. 2009, 47: 827-837. 10.1016/j.plaphy.2009.04.003.
Thuerig B, Binder A, Boller T, Guyer U, Jiménez S, Rentsch C, Tamm L: An aqueous extract of the dry mycelium of Penicillium chrysogenum induces resistance in several crops under controlled and field conditions. Eur J Plant Pathol. 2006, 114: 185-197. 10.1007/s10658-005-4512-6.
Segarra G, van der Ent S, Trillas I, Pieterse CMJ: MYB72, a node of convergence in induced systemic resistance triggered by a fungal and a bacterial beneficial microbe. Plant Biol. 2009, 11: 90-96. 10.1111/j.1438-8677.2008.00162.x.
Shoresh M, Harman GE, Mastouri F: Induced systemic resistance and plant responses to fungal biocontrol agents. Annu Rev Phytopathol. 2010, 48: 21-43. 10.1146/annurev-phyto-073009-114450.
Bailey BA, Bae H, Strem MD, Roberts DP, Thomas SE, Crozier J, Samuels GJ, Choi IY, Holmes KA: Fungal and plant gene expression during the colonization of cacao seedlings by endophytic isolates of four Trichoderma species. Planta. 2006, 224: 1449-1464. 10.1007/s00425-006-0314-0.
Alfano G, Ivey ML, Cakir C, Bos JI, Miller SA, Madden LV, Kamoun S, Hoitink HA: Systemic modulation of gene expression in tomato by Trichoderma hamatum 382. Phytopathology. 2007, 97: 429-437. 10.1094/PHYTO-97-4-0429.
Brotman Y, Lisec J, Michaël M, Chet I, Willmitzer L, Viterbo A: Transcript and metabolite analysis of the Trichoderma-induced systemic resistance response to Pseudomonas syringae in Arabidopsis thaliana. Microbiology. 2012, 158: 139-146. 10.1099/mic.0.052621-0.
Mathys J, De Cremer K, Timmermans P, Van Kerckhove S, Lievens B, Vanhaecke M, Cammue BP, De Coninck B: Genome-wide characterization of ISR induced in Arabidopsis thaliana by Trichoderma hamatum T382 against Botrytis cinerea infection. Front Plant Sci. 2012, 3: 108-
Morán-Diez E, Rubio B, Domínguez S, Hermosa R, Monte E, Nicolás C: Transcriptomic response of Arabidopsis thaliana after 24h incubation with the biocontrol fungus Trichoderma harzianum. J Plant Physiol. 2012, 169: 614-620. 10.1016/j.jplph.2011.12.016.
Pieterse CMJ, Leon-Reyes A, Van der Ent S, Van Wees SC: Networking by small-molecule hormones in plant immunity. Nat Chem Biol. 2009, 5: 308-316. 10.1038/nchembio.164.
Pieterse CMJ, Van Pelt JA, Ton J, Parchmann S, Mueller MJ, Buchala AJ, Métraux JP, Van Loon LC: Rhizobacteria-mediated induced systemic resistance (ISR) in Arabidopsis requires sensitivity to jasmonate and ethylene but is not accompanied by an increase in their production. Physiol Mol Plant Pathol. 2000, 57: 123-134. 10.1006/pmpp.2000.0291.
Van Wees SCM, Van der Ent S, Pieterse CMJ: Plant immune responses triggered by beneficial microbes. Curr Opin Plant Biol. 2008, 11: 443-448. 10.1016/j.pbi.2008.05.005.
van der Ent S, Koornneef A, Ton J, Pieterse CMJ: Annual plant reviews - molecular aspects of plant disease resistance. Edited by: Parker JE. 2009, Oxford, UK: Wiley-Blackwell, 334-370. Induced resistance-orchestrating defence mechanisms through crosstalk and priming, 34.
Conrath U, Beckers GJM, Flors V, Garcia-Agustin P, Jakab G, Mauch F, Newman MA, Pieterse CMJ, Poinssot B, Pozo MJ, Pugin A, Schaffrath U, Ton J, Wendehenne D, Zimmerli L, Mauch-Mani B: Priming: getting ready for battle. Mol Plant Microbe Interact. 2006, 19: 1062-1071. 10.1094/MPMI-19-1062.
Verhagen BW, Glazebrook J, Zhu T, Chang HS, van Loon LC, Pieterse CMJ: The transcriptome of rhizobacteria-induced systemic resistance in Arabidopsis. Mol Plant Microbe Interact. 2004, 17: 895-908. 10.1094/MPMI.2004.17.8.895.
van Hulten M, Pelser M, van Loon LC, Pieterse CMJ, Ton J: Costs and benefits of priming for defense in Arabidopsis. Proc Natl Acad Sci USA. 2006, 103: 5602-5607. 10.1073/pnas.0510213103.
Walters DR, Heil M: Costs and trade-offs associated with induced resistance. Physiol Mol Plant Pathol. 2007, 71: 3-17. 10.1016/j.pmpp.2007.09.008.
Beckers GJM, Conrath U: Priming for stress resistance: from the lab to the field. Curr Opin Plant Biol. 2007, 10: 425-431. 10.1016/j.pbi.2007.06.002.
Perazzolli M, Raoatti B, Bozza E, Pertot I: Trichoderma harzianum T39 induces resistance against downy mildew by priming for defence without costs for grapevine. Biol Control. 2011, 58: 74-82. 10.1016/j.biocontrol.2011.04.006.
Beckers GJ, Jaskiewicz M, Liu Y, Underwood WR, He SY, Zhang S, Conrath U: Mitogen-activated protein kinases 3 and 6 are required for full priming of stress responses in Arabidopsis thaliana. Plant Cell. 2009, 21: 944-953. 10.1105/tpc.108.062158.
Verhagen BW, Trotel-Aziz P, Couderchet M, Hofte M, Aziz A: Pseudomonas spp.-induced systemic resistance to Botrytis cinerea is associated with induction and priming of defence responses in grapevine. J Exp Bot. 2010, 61: 249-260. 10.1093/jxb/erp295.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 621-628. 10.1038/nmeth.1226.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28: 511-515. 10.1038/nbt.1621.
Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.
Mizuno H, Kawahara Y, Sakai H, Kanamori H, Wakimoto H, Yamagata H, Oono Y, Wu J, Ikawa H, Itoh T, Matsumoto T: Massive parallel sequencing of mRNA in identification of unannotated salinity stress-inducible transcripts in rice (Oryza sativa L.). BMC Genomics. 2010, 11: 683-10.1186/1471-2164-11-683.
Zenoni S, Ferrarini A, Giacomelli E, Xumerle L, Fasoli M, Malerba G, Bellin D, Pezzotti M, Delledonne M: Characterization of transcriptional complexity during berry development in Vitis vinifera using RNA-Seq. Plant Physiol. 2010, 152: 1787-1795. 10.1104/pp.109.149716.
Chong J, Le Henanff G, Bertsch C, Walter B: Identification, expression analysis and characterization of defense and signaling genes in Vitis vinifera. Plant Physiol Biochem. 2008, 46: 469-481. 10.1016/j.plaphy.2007.09.010.
Mzid R, Marchive C, Blancard D, Deluc L, Barrieu F, Corio-Costet MF, Drira N, Hamdi S, Lauvergeat V: Overexpression of VvWRKY2 in tobacco enhances broad resistance to necrotrophic fungal pathogens. Physiol Plant. 2007, 131: 434-447. 10.1111/j.1399-3054.2007.00975.x.
Milli A, Cecconi D, Bortesi L, Persi A, Rinalducci S, Zamboni A, Zoccatelli G, Lovato A, Zolla L, Polverari A: Proteomic analysis of the compatible interaction between Vitis vinifera and Plasmopara viticola. J Proteomics. 2012, 75: 1284-1302. 10.1016/j.jprot.2011.11.006.
Shoresh M, Gal-On A, Leibman D, Chet I: Characterization of a mitogen-activated protein kinase gene from cucumber required for trichoderma-conferred plant resistance. Plant Physiol. 2006, 142: 1169-1179. 10.1104/pp.106.082107.
Moreno CA, Castillo F, González A, Bernal D, Jaimes Y, Chaparro M, González C, Rodriguez F, Restrepo S, Cotes AM: Biological and molecular characterization of the response of tomato plants treated with Trichoderma koningiopsis. Physiol Mol Plant Pathol. 2009, 74: 111-120. 10.1016/j.pmpp.2009.10.001.
Palmieri MC, Perazzolli M, Matafora V, Moretto M, Bachi A, Pertot I: Proteomic analysis of grapevine resistance induced by Trichoderma harzianum T39 reveals specific defence pathways activated against downy mildew. J Exp Bot. 2012, 63: 6237-6251. 10.1093/jxb/ers279.
Lukasik E, Takken FL: STANDing strong, resistance proteins instigators of plant defence. Curr Opin Plant Biol. 2009, 12: 427-436. 10.1016/j.pbi.2009.03.001.
Thilmony R, Underwood W, He SY: Genome-wide transcriptional analysis of the Arabidopsis thaliana interaction with the plant pathogen Pseudomonas syringae pv. tomato DC3000 and the human pathogen Escherichia coli O157:H7. Plant J. 2006, 46: 34-53. 10.1111/j.1365-313X.2006.02725.x.
Korolev N, Rav David D, Elad Y: The role of phytohormones in basal resistance and Trichoderma-induced systemic resistance to Botrytis cinerea in Arabidopsis thaliana. BioControl. 2008, 53: 667-683. 10.1007/s10526-007-9103-3.
Van der Ent S, Verhagen BW, Van Doorn R, Bakker D, Verlaan MG, Pel MJ, Joosten RG, Proveniers MC, Van Loon LC, Ton J, Pieterse CM: MYB72 is required in early signaling steps of rhizobacteria-induced systemic resistance in Arabidopsis. Plant Physiol. 2008, 146: 1293-1304. 10.1104/pp.107.113829.
Wu A, Allu AD, Garapati P, Siddiqui H, Dortay H, Zanor MI, Asensi-Fabado MA, Munne-Bosch S, Antonio C, Tohge T, Fernie AR, Kaufmann K, Xue GP, Mueller-Roeber B, Balazadeh S: JUNGBRUNNEN1, a reactive oxygen species-responsive NAC transcription factor, regulates longevity in Arabidopsis. Plant Cell. 2012, 24: 482-506. 10.1105/tpc.111.090894.
Mohr TJ, Mammarella ND, Hoff T, Woffenden BJ, Jelesko JG, McDowell JM: The Arabidopsis downy mildew resistance gene RPP8 is induced by pathogens and salicylic acid and is regulated by W box cis elements. Mol Plant Microbe Interact. 2010, 23: 1303-1315. 10.1094/MPMI-01-10-0022.
Conrath U: Molecular aspects of defence priming. Trends Plant Sci. 2011, 16: 524-531. 10.1016/j.tplants.2011.06.004.
Jaskiewicz M, Conrath U, Peterhansel C: Chromatin modification acts as a memory for systemic acquired resistance in the plant stress response. EMBO Rep. 2011, 12: 50-55. 10.1038/embor.2010.186.
Slaughter A, Daniel X, Flors V, Luna E, Hohn B, Mauch-Mani B: Descendants of primed Arabidopsis plants exhibit resistance to biotic stress. Plant Physiol. 2012, 158: 835-843. 10.1104/pp.111.191593.
Bultreys A, Trombik T, Drozak A, Boutry M: Nicotiana plumbaginifolia plants silenced for the ATP-binding cassette transporter gene NpPDR1 show increased susceptibility to a group of fungal and oomycete pathogens. Mol Plant Pathol. 2009, 10: 651-663. 10.1111/j.1364-3703.2009.00562.x.
Kamoun S: A catalogue of the effector secretome of plant pathogenic oomycetes. Annu Rev Phytopathol. 2006, 44: 41-60. 10.1146/annurev.phyto.44.070505.143436.
Stassen JH, Van den Ackerveken G: How do oomycete effectors interfere with plant life?. Curr Opin Plant Biol. 2011, 14: 407-414. 10.1016/j.pbi.2011.05.002.
Birch PR, Rehmany AP, Pritchard L, Kamoun S, Beynon JL: Trafficking arms: oomycete effectors enter host plant cells. Trends Microbiol. 2006, 14: 8-11. 10.1016/j.tim.2005.11.007.
Mestre P, Piron MC, Merdinoglu D: Identification of effector genes from the phytopathogenic oomycete Plasmopara viticola through the analysis of gene expression in germinated zoospores. Fungal Biol. 2012, 116: 825-835. 10.1016/j.funbio.2012.04.016.
Díez-Navajas AM, Wiedemann-Merdinoglu S, Greif C, Merdinoglu D: Nonhost versus host resistance to the grapevine downy mildew, Plasmopara viticola, studied at the tissue level. Phytopathology. 2008, 98: 776-780. 10.1094/PHYTO-98-7-0776.
Malacarne G, Perazzolli M, Cestaro A, Sterck L, Fontana P, Van de Peer Y, Viola R, Velasco R, Salamini F: Deconstruction of the (paleo)polyploid grapevine genome based on the analysis of transposition events involving NBS resistance genes. PLoS One. 2012, 7: e29762-10.1371/journal.pone.0029762.
Schnable JC, Springer NM, Freeling M: Differentiation of the maize subgenomes by genome dominance and both ancient and ongoing gene loss. Proc Natl Acad Sci USA. 2011, 108: 4069-4074. 10.1073/pnas.1101368108.
Perazzolli M, Bampi F, Faccin S, Moser M, De Luca F, Ciccotti AM, Velasco R, Gessler C, Pertot I, Moser C: Armillaria mellea induces a set of defense genes in grapevine roots and one of them codifies a protein with antifungal activity. Mol Plant Microbe Interact. 2010, 23: 485-496. 10.1094/MPMI-23-4-0485.
EPPO: European and Mediterranean Plant Protection Organization (EPPO): Guidelines for the efficacy evaluation of fungicides: Plasmopara viticola. EPPO Bull. 2001, 31: 313-317.
Unger S, Buche C, Boso S, Kassemeyer HH: The course of colonization of two different Vitis genotypes by Plasmopara viticola indicates compatible and incompatible host-pathogen interactions. Phytopathology. 2007, 97: 780-786. 10.1094/PHYTO-97-7-0780.
Sequence read archive (SRA).http://www.ncbi.nlm.nih.gov/sra/,
Velasco R, Zharkikh A, Troggio M, Cartwright DA, Cestaro A, Pruss D, Pindo M, Fitzgerald LM, Vezzulli S, Reid J, Malacarne G, Iliev D, Coppola G, Wardell B, Micheletti D, Macalma T, Facci M, Mitchell JT, Perazzolli M, Eldredge G, Gatto P, Oyzerski R, Moretto M, Gutin N, Stefanini M, Chen Y, Segala C, Davenport C, Dematte L, Mraz A, et al: A high quality draft consensus sequence of the genome of a heterozygous grapevine variety. PLoS One. 2007, 2: e1326-10.1371/journal.pone.0001326.
Pinot noir grapevine genome browser.http://genomics.research.iasma.it,
Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25: 1105-1111. 10.1093/bioinformatics/btp120.
Roberts A, Trapnell C, Donaghey J, Rinn JL, Pachter L: Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. 2011, 12: R22-10.1186/gb-2011-12-3-r22.
Roberts A, Pimentel H, Trapnell C, Pachter L: Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics. 2011, 27: 2325-2329. 10.1093/bioinformatics/btr355.
Saeed AI, Bhagabati NK, Braisted JC, Liang W, Sharov V, Howe EA, Li J, Thiagarajan M, White JA, Quackenbush J: TM4 microarray software suite. Methods Enzymol. 2006, 411: 134-193.
Fontana P, Cestaro A, Velasco R, Formentin E, Toppo S: Rapid annotation of anonymous sequences from genome projects using semantic similarities and a weighting scheme in gene ontology. PLoS One. 2009, 4: e4619-10.1371/journal.pone.0004619.
Falda M, Toppo S, Pescarolo A, Lavezzo E, Di Camillo B, Facchinetti A, Cilia E, Velasco R, Fontana P: Argot2: a large scale function prediction tool relying on semantic similarity of weighted Gene Ontology terms. BMC Bioinforma. 2012, 13 (Suppl 4): S14-10.1186/1471-2105-13-S4-S14.
The gene ontology.http://www.geneontology.org/,
Beissbarth T, Speed TP: GOstat: find statistically overrepresented Gene Ontologies within a group of genes. Bioinformatics. 2004, 20: 1464-1465. 10.1093/bioinformatics/bth088.
Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, Vezzi A, Legeai F, Hugueney P, Dasilva C, Horner D, Mica E, Jublot D, Poulain J, Bruyere C, Billault A, Segurens B, Gouyvenoux M, Ugarte E, Cattonaro F, Anthouard V, Vico V, Del Fabbro C, Alaux M, Di Gaspero G, Dumas V, et al: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449: 463-467. 10.1038/nature06148.
Grape genome browser.http://www.genoscope.cns.fr/externe/GenomeBrowser/Vitis/,
Luu-The V, Paquet N, Calvo E, Cumps J: Improved real-time RT-PCR method for high-throughput measurements using second derivative calculation and double correction. Biotechniques. 2005, 38: 287-293. 10.2144/05382RR05.
Ruijter JM, Ramakers C, Hoogaars WM, Karlen Y, Bakker O, van den Hoff MJ, Moorman AF: Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 2009, 37: e45-10.1093/nar/gkp045.
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.
Gamm M, Héloir MC, Kelloniemi J, Poinssot B, Wendehenne D, Adrian M: Identification of reference genes suitable for qRT-PCR in grapevine and application for the study of the expression of genes involved in pterostilbene synthesis. Mol Genet Genomics. 2011, 285: 273-285. 10.1007/s00438-011-0607-2.
This research was supported by the Post-Doc Project 2006 Resistevite, Envirochange and Accordo di programma funded by the Autonomous Province of Trento. The authors thank Dr. Yigal Elad (The Volcani Center, Bet Dagan, Israel) for supplying Trichodex and Dr. Laurent Farinelli (Fasteris SA, Plan-les-Ouates, Switzerland) for technical support with the RNA-Seq Illumina sequencing.
The authors declare that they have no competing interests.
MP performed the grapevine treatments, P. viticola inoculations, RNA extractions, and real-time RT-PCR experiments. MP, MM, PF, and AF analysed the RNA-Seq sequences by mapping, gene expression estimation, and statistical analysis. MP, CM, RV, MD, and IP contributed to data interpretation and manuscript writing. MP, CM, MD, and IP conceived the study, designed the experiment, and coordinated all research activities. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Numbers of RNA-Seq reads passing the quality check and mapping to the grapevine genome are reported for each sequencing replicate (named A and B) of each biological replicate (numbered from 1 to 3) for control (C), Trichoderma harzianum T39-treated (T39), Plasmopara viticola-inoculated control (C+P.v.), and P. viticola-inoculated T39-treated (T39+P.v.) plants. (PDF 11 KB)
RNA-Seq sequencing and coverage of the grapevine transcriptome for each biological and sequencing replicate.
Additional file 2: Total bases (Mbp) sequenced by RNA-Seq analysis and coverage of the grapevine transcriptome are reported for each sequencing replicate (named A and B) of each biological replicate (numbered from 1 to 3) for control (C), Trichoderma harzianum T39-treated (T39), Plasmopara viticola-inoculated control (C+P.v.), and P. viticola-inoculated T39-treated (T39+P.v.) plants. (PDF 7 KB)
Additional file 3: (A) Proportions of reads mapping to unique locations (unique reads, yellow), reads mapping to multiple locations with 2–100 matches (multi-reads, red), and reads not mapping or ambiguously mapping to more than 100 locations (unused reads, white) to the grapevine genome. The mean number of sequenced reads and standard errors for three biological replicates are presented for each treatment. Percentages (%) are calculated with respect to the total filtered reads. Unique reads and multi-reads were used for evaluating gene expression. (B) Proportions of reads mapping to grapevine genes (green) and to intergenic regions (red) of the grapevine genome. The mean number of sequenced reads and standard errors for three biological replicates are presented for each treatment. Percentages (%) are calculated with respect to the total mapping reads. Grapevine treatments: control (C), Trichoderma harzianum T39-treated (T39), Plasmopara viticola-inoculated control (C+P.v.), and P. viticola-inoculated T39-treated (T39+P.v) plants. (PDF 5 KB)
Additional file 4: Numbers of RNA-Seq reads mapping to grapevine genes are reported for each sequencing replicate (named A and B) of each biological replicate (numbered from 1 to 3) for control (C), Trichoderma harzianum T39-treated (T39), Plasmopara viticola-inoculated control (C+P.v.), and P. viticola-inoculated T39-treated (T39+P.v.) plants. (PDF 8 KB)
Additional file 5: Numbers of expressed grapevine genes and novel genes with respect to the grapevine annotation are reported for each sequencing replicate (named A and B) of each biological replicate (numbered from 1 to 3) for control (C), Trichoderma harzianum T39-treated (T39), Plasmopara viticola-inoculated control (C+P.v.), and P. viticola-inoculated T39-treated (T39+P.v.) plants (PDF 7 KB)
Additional file 6: Gene expression values (FPKM) and standard errors are reported for all (33,514) Pinot Noir grapevine predicated genes  Release 3  in control (C), Trichoderma harzianum T39-treated (T39), Plasmopara viticola-inoculated control (C+P.v.), and P. viticola-inoculated T39-treated (T39+P.v.) plants. Expression values in each biological replicate (numbered from 1 to 3) and each sequencing replicate (named A and B) are also listed. (XLS 13 MB)
Additional file 7: Comparison of the expression levels of all grapevine genes, expressed as fragments per kilobase of transcript per million fragments mapped (FPKM), in the two sequencing replicates (named A and B) of (A) control (C) biological replicate no. 1; (B) C biological replicate no. 2 (B); (C) C biological replicate no. 3; (D) Trichoderma harzianum T39-treated (T39) biological replicate no. 1; (E) T39-treated biological replicate no. 2; (F) T39-treated biological replicate no. 3; (G) Plasmopara viticola-inoculated control (C+P.v.) biological replicate no. 1; (H) C+P.v. biological replicate no. 2; (I) C+P.v. biological replicate no. 3; (J) P. viticola-inoculated T39-treated plants (T39+P.v) biological replicate no. 1; (K) T39+P.v biological replicate no. 2; and (L) T39+P.v biological replicate no. 3. (PDF 2 MB)
Additional file 8: Principal component analysis (PCA) is based on the expression values (FPKM) of all grapevine genes for each sequencing replicate (named A and B) of each biological replicate (numbered from 1 to 3) for control (C), Trichoderma harzianum T39-treated (T39), Plasmopara viticola-inoculated control (C+P.v.), and P. viticola-inoculated T39-treated (T39+P.v.) plants. (PDF 6 KB)
Pearson’s correlation coefficients between sequencing and biological replicates of RNA-Seq analysis.
Additional file 9: Pearson’s correlation coefficients are based on the gene expression values (FPKM) of all grapevine genes in two sequencing replicates (A and B) of each biological replicate (numbered from 1 to 3) for control (C), Trichoderma harzianum T39-treated (T39), Plasmopara viticola-inoculated control (C+P.v.), and P. viticola-inoculated T39-treated (T39+P.v.) plants. (PDF 7 KB)
Additional file 10: Gene expression values (FPKM) and standard errors are reported for differentially expressed (7,024) genes in control (C), Trichoderma harzianum T39-treated (T39), Plasmopara viticola- inoculated control (C+P.v.), and P. viticola-inoculated T39-treated (T39+P.v.) plants. Differentially expressed genes were identified by the DESeq package with a false discovery rate (FDR) of 5% and a fold-change greater than two in at least one pairwise comparison. Fold-changes, clustering results, and functional annotation are reported for each gene. (XLS 4 MB)
Additional file 11: Forward and reverse primer sequences of real-time RT-PCR analysis are reported for 24 selected grapevine genes. Grapevine Actin (TC81781) [4, 52] and VATP16 (XM_002269086.1)  were used as constitutive genes for normalising the real-time RT-PCR data. (XLS 32 KB)
Additional file 12: actin gene family. Gene expression values (FPKM) and standard errors are reported for grapevine genes belonging to the Actin gene family in control (C), Trichoderma harzianum T39-treated (T39), Plasmopara viticola-inoculated control (C+P.v.), and P. viticola-inoculated T39-treated (T39+P.v) plants. (PDF 8 KB)
About this article
Cite this article
Perazzolli, M., Moretto, M., Fontana, P. et al. Downy mildew resistance induced by Trichoderma harzianum T39 in susceptible grapevines partially mimics transcriptional changes of resistant genotypes. BMC Genomics 13, 660 (2012). https://doi.org/10.1186/1471-2164-13-660
- Induced resistance
- Next generation sequencing
- Gene expression
- Vitis vinifera
- Plant-pathogen interactions