Original Research ARTICLE
Transcriptome analysis reveals regulatory networks underlying differential susceptibility to Botrytis cinerea in response to nitrogen availability in Solanum lycopersicum
- 1Facultad de Agronomía e Ingeniería Forestal, Pontificia Universidad Católica de Chile, Santiago, Chile
- 2Millennium Nucleus Center for Plant Systems and Synthetic Biology, Santiago, Chile
- 3Departamento de Genética Molecular y Microbiología, Facultad de Ciencias Biológicas, Pontificia Universidad Católica de Chile, Santiago, Chile
- 4Millennium Nucleus for Fungal Integrative and Synthetic Biology, Santiago, Chile
- 5Facultad de Ciencias, Instituto de Bioquímica y Microbiología, Universidad Austral de Chile, Valdivia, Chile
Nitrogen (N) is one of the main limiting nutrients for plant growth and crop yield. It is well documented that changes in nitrate availability, the main N source found in agricultural soils, influences a myriad of developmental programs and processes including the plant defense response. Indeed, many agronomical reports indicate that the plant N nutritional status influences their ability to respond effectively when challenged by different pathogens. However, the molecular mechanisms involved in N-modulation of plant susceptibility to pathogens are poorly characterized. In this work, we show that Solanum lycopersicum defense response to the necrotrophic fungus Botrytis cinerea is affected by plant N availability, with higher susceptibility in nitrate-limiting conditions. Global gene expression responses of tomato against B. cinerea under contrasting nitrate conditions reveals that plant primary metabolism is affected by the fungal infection regardless of N regimes. This result suggests that differential susceptibility to pathogen attack under contrasting N conditions is not only explained by a metabolic alteration. We used a systems biology approach to identify the transcriptional regulatory network implicated in plant response to the fungus infection under contrasting nitrate conditions. Interestingly, hub genes in this network are known key transcription factors involved in ethylene and jasmonic acid signaling. This result positions these hormones as key integrators of nitrate and defense against B. cinerea in tomato plants. Our results provide insights into potential crosstalk mechanisms between necrotrophic defense response and N status in plants.
Nitrogen (N) is an essential macronutrient whose availability significantly impacts plant growth and development. Since natural environments and agricultural fields have limited amounts of N, the production of high-yielding crops relies on the application of large quantities of nitrogenous fertilizers, which come at considerable economic (Good et al., 2004) and environmental costs (Hirel et al., 2011).
The relevance of N for plants is clearly exemplified by its effects on leaf growth (von Wirén et al., 2000), senescence (Vanacker et al., 2006), root system architecture (Zhang et al., 2007; Vidal et al., 2010) and flowering time (Castro Marín et al., 2010). Besides growth and developmental effects, it is also clear that N nutrition can impact the plant's ability to cope with environmental challenges such as plant pathogen attacks (Snoeijers et al., 2000; Walters and Bingham, 2007; Dordas, 2008; Fagard et al., 2014). Different studies have shown that N availability impacts the outcome of plant-pathogen interactions, although the mechanisms underlying this connection are poorly understood, and the effect of N on this process is highly dependent on the crop being studied and on the particular life style of the pathogen involved (Snoeijers et al., 2000; Walters and Bingham, 2007; Dordas, 2008; Fagard et al., 2014). Consequently, it is difficult to derive general rules for the role of N in this process.
The plant defense response is a complex biological process involving numerous changes at the biochemical, physiological, and molecular (transcriptional) level, all governed by an intricate grid of hierarchical and regulatory interactions (Windram et al., 2014). These defense mechanisms are triggered partly by the defense hormones salicylic acid (SA) and jasmonic acid (JA) (Thomma et al., 2001; Glazebrook, 2005; Pieterse et al., 2009; Caarls et al., 2015). In general terms, SA is involved in resistance against biotrophic pathogens (Vlot et al., 2009), while JA participates in the regulation of defense response against necrotrophic pathogens and insects (Farmer et al., 2003).
Since plant defense is an active and energetically costly response mechanism, it is expected that the metabolic state of the plant plays a fundamental role in the outcome of the plant-pathogen interaction. A few agronomic reports indicate that high N availability increases the incidence of crop diseases (Hoffland et al., 2000; Olesen et al., 2003; Ballini et al., 2013). On the other hand, some studies report that a reduction in N fertilization increases disease severity (Hoffland et al., 1999; Long et al., 2007; Linquist et al., 2008). In the model plant Arabidopsis thaliana, proteins involved in plant resistance to infections are up-regulated in response to changes in N levels (Lau and Hamer, 1996; Dietrich et al., 2004). In the case of Botrytis cinerea, one of the most important fungal plant pathogens with regard to both its scientific and agronomic importance (Dean et al., 2012), plant N status can either promote or impede infection, depending on the plant species. For instance, while high N fertilization rates increase disease severity in legumes (Davidson et al., 2007) and lettuce (Lecompte et al., 2013), elevated N concentrations result in reduced susceptibility to this fungus in tomato (Hoffland et al., 1999; Lecompte et al., 2010). As there appears to be a trade-off between plant growth and defense responses (Walters and Heil, 2007), an intricate interconnection between metabolic and stress signaling pathways is required for proper and efficient resource allocation (Hey et al., 2010). The relationship between N metabolism and plant defense responses however, has not been analyzed in detail, although it has recently been recognized that this interconnection may shed new light on the complexity of plant defense strategies (Fagard et al., 2014).
In this study, we analyzed changes in global gene expression patterns during B. cinerea infections in tomato plants grown under contrasting nitrate regimes, with the goal of characterizing the interaction between N supply and defense responses at molecular level. High nitrate availability reduced plant susceptibility to this fungus when measured in leaves and tomato fruits. Global gene expression patterns confirmed that the tomato plant defense response is affected by nitrate availability. To get a broader overview of the interconnection between N metabolism and the plant defense response, a gene network analysis was performed. This strategy identified a transcriptional regulatory network controlling plant susceptibility to this fungus depending on nitrate condition. After validating a network-derived prediction by RT-qPCR in leaves and fruits infected by the fungus, we conclude that the expression of key transcription factors (TFs) involved in ethylene (ET) and JA signaling is modulated by the plant N status when tomato is infected by B. cinerea, suggesting that these hormones play a role in the nitrate-defense response interaction in tomato. Our study identified crosstalk points between N-nutrition, defense response and the ET/JA pathways in plants.
Plant Nitrate Regimes Alter the Progression of B. cinerea Infection in Tomato
As a first step to evaluate a connection between plant nitrate availability and plant response to pathogen infection, we evaluated the growth of S. lycopercicum cv. MicroTom under N (nitrate) conditions that ranged from limiting to sufficient. Tomato plants were grown in pots with vermiculite, an inert growing substrate without N sources, and irrigated with a complete mineral nutrient solution without N supplemented with 2, 4, 6, or 12 mM nitrate (final concentration). As expected, nitrate concentration had a significant impact on tomato plant growth. Maximal growth was attained with 6 mM (N-sufficient condition) under our experimental conditions, and 2 mM and 4 mM nitrate produced severe and mild growth phenotypes, respectively (N-limiting conditions). To evaluate higher N input, tomato plants were grown using 12 mM nitrate. The amount of shoot biomass was reduced upon lowering the amount of nitrate to that normally present in the N-sufficient condition (6 mM), while higher concentrations did not lead to significant changes (Figure S1).
To evaluate whether contrasting nitrate concentrations impact the susceptibility of tomato plants to fungal infection, plants were grown under 2, 4, 6, or 12 mM nitrate as the only N source and challenged with B. cinerea. Leaves in planta were inoculated with an aqueous suspension of 5 × 103 conidia. Typical B. cinerea symptoms, such as necrotic lesions, were observed in leaves under all N regimes. The first visual symptoms of the infections were detected 2 days post inoculation (dpi), at which time darkening of the leaf surface under the inoculum was observed (data not shown). As shown in Figure 1, disease symptoms developed faster in plants grown under N-limiting than in N-sufficient conditions. At 3 dpi, larger lesions were observed under N-limiting conditions, with evident tissue maceration surrounding primary infection sites in leaves (Figure 1A). Conversely, only discrete necrotic lesions were observed under N-sufficient conditions (Figure 1A). Even though, disease symptoms (expanding necrosis, chlorosis, and tissue maceration) were observed in leaves from plants grown under all N regimes at 5 dpi (data not shown), the size of the lesions and the percentage of the leaf exhibiting symptoms were always larger in leaves from plants grown under N-limiting conditions (Figures 1B,C).
Figure 1. Tomato leaves susceptibility to B. cinerea in plants grown under different nitrate concentrations. (A) Representative inoculated leaves (3 dpi) for each nitrate condition used throughout this study. (B) Lesion size, measured as diameter of expanding lesions, for the indicated nitrate concentrations. (C) Disease incidence (percentage of leaf area with disease symptoms). Different letters indicate significant differences among treatments at a given time point (p ≤ 0.05; error bars indicate SEM; n = 6).
To quantitatively assess fungal growth on plant tissue, we employed a quantitative PCR (qPCR) assay based on the relative quantification of fungal and plant DNA in infected plants, as described (Gachon and Saindrenan, 2004). As shown in Figure S2A, and consistent with the results shown in Figure 1, an increase in fungal growth was observed in all N-conditions as disease progressed, with larger quantities in infected plants grown under N-limiting conditions.
Previous studies have shown that nitrate responses are organ and developmental stage-dependent (Wang et al., 2003, 2004; Vidal et al., 2014). To evaluate whether there was a difference in the aforementioned results when fruits instead of vegetative plant tissue was used, and considering that ripening promotes fruit susceptibility to pathogens (Alba et al., 2005; Giovannoni, 2007; Cantu et al., 2009), we evaluated the response of unripe green and ripe red tomato fruits to B. cinerea infection. Figure 2 shows disease symptoms and lesion progression in two developmental stages known as “mature green” (MG) and “red ripe” (RR) fruits obtained from plants grown under the same N conditions described above. Under all nitrate conditions evaluated, MG tomatoes were significantly less susceptible to B. cinerea compared to RR fruits. When MG tomato fruits were inoculated, a small necrotic lesion was observed at the site of infection 3 dpi (Figure 2A). On RR tomato fruits, on the other hand, tissue rotting, and fungal growth were already evident 3 dpi and extended into the pericarp tissue (Figure 2A). The severity of the disease symptoms (Figures 2B,C) and the accumulation of fungal biomass (Figures S2B,C) in infected MG and RR fruits from plants grown under N-limiting conditions were larger than those obtained from a N-sufficient regime. These results indicate that plants grown under N-limiting conditions exhibit enhanced susceptibility to B. cinerea infection.
Figure 2. Tomato fruits susceptibility to B. cinerea obtained from plants grown under different nitrate concentrations. (A) Representative inoculated fruits (3 dpi) for each nitrate concentration and ripening stage (dpa, days post-anthesis). (B,C) Lesion size (diameter) for inoculated MG and RR fruits, respectively. For MG fruits, white arrows indicate the site of infection with minimal necrotic lesion. Different letters indicate significant differences among treatments at a given time point (p ≤ 0.05; error bars indicate SEM; n = 4).
Global Gene Expression Analysis Supports a Connection between Plant N Status and Susceptibly to Fungal Infection in Tomato
To better understand the molecular changes underlying the impact of nitrate availability on plant susceptibility to B. cinerea, we performed transcriptome-profiling assays on mock-treated and infected plants grown under severe limiting (2 mM) and sufficient (6 mM) N conditions, using GeneChip Tomato Genome Arrays (Affymetrix). We also included 12 mM nitrate in our analysis, as a higher N input. These experiments were performed using shoots since most plant defense responses have been described in this particular vegetative tissue (Glazebrook, 2005). As shown in Figure S3, our analysis suggests a significant interaction between the plant N-status and B. cinerea infection (Two-way ANOVA, p < 0.01), where 2110 genes showed significant Nitrogen/B. cinerea (N/B) interaction, representing approximately 4.7% of the analyzed tomato transcriptome (Table S1). Since this result suggests a nitrate and defense signaling crosstalk in response to B. cinerea infection, this set of genes was selected for further analysis.
Principal component analysis (PCA) (Figure 3A) showed that the first principal component (PC1) accounted for 45% of the total variation and segregated mock (M) from infected (B) samples, regardless of N-conditions. In addition, PC1 also allows the differentiation between N-status in B. cinerea-infected samples. The second principal component (PC2) accounted for 13% of the variation, and differentiated limiting vs. sufficient plant N-status in mock-treated samples. Indeed, mock-treated samples from plants grown under N-limiting conditions (2 mM nitrate) cluster more closely with B. cinerea-infected samples, rather than with non-infected samples, under both N-sufficient concentrations (6 mM and 12 mM nitrate). This result suggests that the transcriptome state of infected plants may be similar to N-limiting conditions. Moreover, the total number of differentially expressed genes (DEG) in response to pathogen infection was comparable among all N conditions evaluated (Figure 3B). For DEG analysis, genes with significant N/B model were used for a pairwise comparison between mock and infected samples in each N-condition, using Rank-products (p < 0.05). However, only 18% of up-regulated (22/123) genes were common to all N concentrations analyzed (Figure 3C). Although, the susceptibility phenotypes of plant grown under 2 mM and 6 mM nitrate were different in our experimental conditions, these N-conditions present more common up-regulated and down-regulated genes when compared with plants grown in 12 mM nitrate. These results suggest that higher N-input also affect plant defense response, indicating a complex link between N-metabolism and B. cinerea infections.
Figure 3. The plant N-availability affects global gene expression in response to B. cinerea infection in tomato. (A) Principal component analysis (PCA) of RMA-normalized global gene expression data demonstrate a connection between plant N-status and B. cinerea susceptibility. Colored dots denote each biological replicate. “M” and “B” indicate mock-treated and B. cinerea-infected plants, respectively, while the number next to each dot indicate the N concentration. (B) Number of genes that are either up- or down-regulated by the B. cinerea infection in each N condition analyzed in comparison to mock-inoculated plants (ANOVA; Rank-product, adjusted p < 0.05). (C) Venn diagrams of the data presented in (B) showing the overlap of genes that are either up- or down-regulated in tomato inoculated with B. cinerea in each N condition comparing with corresponding mock-treated samples (B/M).
Differentially Expressed Genes under N-sufficient Condition supports better Defense Response in Tomato
We performed a functional classification of the DEG in each N-condition (Tables S2–S4) using the Generic Gene Ontology (GO) Term Finder tool (Boyle et al., 2004). This analysis showed a wide range of biological processes that were affected due to fungal infection under different N-conditions (Figure 4). Processes associated with metabolism, both primary and secondary, were significantly affected by the pathogen under all N conditions analyzed (p < 0.01). This observation is consistent with previous reports that show that pathogens can reprogram host metabolism, strongly affecting primary, and secondary metabolism in plants (Kliebenstein et al., 2005; Rojas et al., 2014). Figure 4B shows a general overview of metabolic pathways affected by B. cinerea infections in each N-condition, using MapMan software (Thimm et al., 2004). Cell wall metabolisms were strongly affected by the fungus infections in plants grown in N-limiting conditions. In contrast, secondary metabolisms were most affected by the fungus infections in N-sufficient conditions. This finding is consistent with the fact that the accumulation of secondary metabolites is a plant defense mechanism triggered by several fungus pathogens (Kliebenstein et al., 2005; Ward et al., 2010; Smith et al., 2014; Pusztahelyi et al., 2015). Interestingly, we found that genes associated with N transporters (NRT2), N and amino acid metabolism were induced in response to B. cinerea infection depending on the plant's N-status, with it occurring most in N-sufficient conditions (Figure 4B). In agreement with this observation, it has been shown in a recent report that the expression of several genes involved in N metabolism and its transport are strongly affected in response to infections with B. cinerea, Phytophthora infestans, Phytophthora parasitica, and Pseudomonas syringae (Fagard et al., 2014). Other biological functions modified by the B. cinerea infection, and specifically over-represented under N-sufficient conditions, were transport and oxidation-reduction processes (Figure 4, p < 0.01). Redox status significantly impacts both plant defenses and B. cinerea infection (Lamb and Dixon, 1997; Torres et al., 2006).
Figure 4. Biological processes affected in response to B. cinerea infection under different N conditions. (A) The distribution of genes into different GO categories for biological processes is shown. The number of genes induced or repressed (in white and black, respectively) by the fungal infection depending on the N condition used relative to mock-treated plants is presented. Only genes with significant differential expression (p < 0.05) were plotted, corresponding to 132, 151, and 103 genes (2, 6, and 12 mM nitrate, respectively). Overrepresented categories compared to the entire set on the Tomato genome array are indicated with asterisks (**p < 0.01, *p < 0.05). (B) MapMan software was used to provide a snapshot of modulated genes over the main metabolic pathways. DEGs were binned to MapMan functional categories and Log2 fold changes values for each gene upon B. cinerea infection referred to mock-treated samples are represented (B/M in 2 mM, 6 mM, and 12 mM nitrate). Up-regulated and down-regulated transcripts are shown in red and green, respectively.
Our enrichment analysis also identified genes that encode proteins previously implicated in defense responses (AbuQamar et al., 2006; Cantu et al., 2009; Windram et al., 2012; Blanco-Ulate et al., 2013). This group of genes was closely examined to evaluate possible correlations between reduced expression of defense-related genes and susceptibility to B. cinerea infection. Table 1 shows a subset of these genes induced in plants grown under N-sufficient conditions or repressed under N-limiting conditions. These genes include anthranilate synthase 1 (ASA1), 2-oxophytodienoate reductase 3 (OPR3), and pathogenesis-related 4 (PR4). This latter gene encodes a protein with antifungal chitin-binding activity, which is repressed under N-limiting conditions and induced under N-sufficient or higher conditions. In Arabidopsis, this gene has been associated with resistance to necrotrophic pathogens (Catinot et al., 2015) and it is known to be induced by B. cinerea infection (AbuQamar et al., 2006). In addition, genes involved in JA and ET response were also identified. The ACC oxidase gene involved in ET biosynthesis was significantly induced during fungal infection in plants grown under N-sufficient conditions, but not at low nitrate concentrations (Table 1). Thus, a differential response of key defense-related genes under contrasting nitrate concentrations might partially explain the differences in susceptibility to B. cinerea in tomato leaves.
Table 1. Expression profile of defense related genes in response to B. cinerea infections under different N conditions.
Network Analysis Identifies B. cinerea-responsive Regulatory Pathways Associated with the Plant N status
To identify functional relationships and key regulatory TFs among the N/B responsive genes, we performed a gene co-expression network analysis. Due to the limited knowledge of gene interactions in tomato, we decided to incorporate data available from A. thaliana. For ortholog group assignment, tomato sequences (gene models ITAG release 2.4, https://solgenomics.net), matching probe sets on the GeneChip Tomato Genome array, were aligned to the Arabidopsis genome using OrthoMCL (Li et al., 2003). For network inference, we used all Arabidopsis loci identifiers to predict putative protein-DNA interactions, making use of sequence preference information recently obtained for a large number of Arabidopsis TFs (Weirauch et al., 2014) and the gene co-expression database (Obayashi et al., 2014). As shown in Figure 5A, a network comprising 477 genes was obtained. Genes are represented as nodes connected by edges that depict predicted regulatory interactions. In this analysis, 9 co-expression modules were identified using the community cluster (Glay) algorithm in ClusterMaker tools (Morris et al., 2011). This algorithm recognizes functionally related groups and finds densely connected regions in a network (Su et al., 2010). In addition, we found statistically overrepresented biological functions in 7 of these modules using a hypergeometric test in BINGO tools (adjusted p < 0.05 by FDR). Interestingly, network includes functions important for plant defense against necrotrophic fungi, such as stress response, oxide-reduction processes and cell wall and wax biosynthesis processes. Moreover, genes related to ET signaling and the responses to this hormone were also identified as overrepresented in these modules, consistent with the finding describing how ET responses are important for B. cinerea resistance in tomato leaves (Díaz et al., 2002). In addition, we found several TFs belonging to the ERF (Ethylene Response Factor) family, which have been implicated in the plant defense response to B. cinerea (AbuQamar et al., 2006; Windram et al., 2012). Nevertheless, susceptibility of tomato fruit to B. cinerea had been also associated in part with ET pathways (Cantu et al., 2009; Blanco-Ulate et al., 2013), suggesting complex and tissue-specific ET regulatory networks in response to fungus infections.
Figure 5. Network analysis of transcriptomics data predicts a B. cinerea-responsive regulatory gene network affected by the plant N-status. (A) Network analysis of genes with significant N/B ANOVA model. Genes are presented as triangles (TFs) and rectangles (target genes). Colors are used to distinguish each gene network module grouped by topology. The most overrepresented GO term (biological process) and TF in each module are indicated. (B) Network analysis of genes with significant N/B ANOVA model in which genes also present in the B ANOVA model are highlighted. Yellow nodes indicate genes that also respond to the infection alone, irrespective of the N concentration. X and Y axes do not represent any particular scale. (C) Subnetwork of TFs and their putative targets derived from most co-expressed genes. Arrowhead green lines indicate predicted transcriptional activation, while red lines indicate transcriptional repression. Nodes are color-coded based on biological processes: metabolism (gray), unknown processes (white), defense, and stress response (yellow), response to hormone (purple), JA/ET processes (light blue), oxidative-reduction processes (green), and transport (orange).
In order to distinguish the effect of the plant's N-status in the defense response against B. cinerea (Figure 5A; N/B ANOVA model) from the general plant response to fungal infection irrespective of its N condition (B ANOVA model), we highlight genes with a significant B (B. cinerea) ANOVA model within the mentioned network (Figure 5B). When Figures 5A,B were compared, we noticed that only 40% of the mentioned N/B network responded to the infection as a single (B) factor, supporting the link between fungus infection and N response. Interestingly, the majority of the TFs present in the network also link to this interaction (in Figure 5B, note that most TFs at the center of each module are not highlighted in yellow).
To identify TFs that may control relevant functions in the N status-dependent response to B. cinerea, we focused on TFs and those putative targets possessing significant co-expression (Figure 5C). We applied a 90th percentile cut-off in co-expression networks considering the absolute value, such that only the highest scoring 10% of edges were selected. A total of 28 TFs are present in the N/B network, from which 22 were highly connected within its center and might represent important regulators of the N/B interaction. Interestingly, TFs within the network and their respective putative targets are mostly involved in transport, oxide-reduction processes and stress and defense responses. In addition, within this group of regulatory proteins, we also identified overrepresented biological functions related to the response to ET and JA.
The Plant N Status Modulates the Expression of TFs Belonging to Families Involved in Plant Defense
The results depicted in Figure 5 support the importance of plant N metabolism in plant defense pathways. To validate this prediction, we analyzed the expression levels of 18 TF-encoding genes in leaves or fruits infected by B. cinerea, from tomato plants grown under all nitrate concentrations used throughout this work (Table S5). An RT-qPCR analysis of selected TFs supported microarray results and network analysis (Figure 6). Hierarchical clustering analysis of the RT-qPCR data revealed 6 clusters. Cluster #1 (Figure 6, in gray), composed of four TF-encoding genes (NAC-like, WRKY15, BIM, and WRKY70), groups genes that are induced by fungal infection under N-sufficient conditions, but repressed in leaves or marginally induced in fruits of plants grown under N-limiting conditions. In Arabidopsis, WRKY 15 and WRKY 70 (Figure 5, Modules 3 and 8, respectively) have been reported as induced when infected with B. cinerea. Interestingly, wrky70 mutant plants showed enhanced susceptibility to this fungus, in a SA dependent manner (AbuQamar et al., 2006).
Figure 6. Hierarchical clustering of mRNA levels for selected TFs showing a differential gene expression pattern in response to the fungus infection under different N conditions. The dendrogram and colored image were produced using gene expression data after RT-qPCR (n = 3). The color scale ranges from saturated green for log2 ratios −3.0 and below, to saturated red for log2 ratios +3.0 and above. Fold change in gene expression for each gene upon B. cinerea infection are referred to mock-treated samples. Each TF is represented by a single row of colored boxes. Six clusters were identified, and are numbered and denoted with a vertical line (from top to bottom) as follows: cluster #1 in gray, #2 in brown, #3 in light blue, #4 in orange, #5 in green, and #6 in blue.
B. cinerea also triggered the induction of zinc-finger TFs (ZAT) genes in leaves, under all N-conditions (Figure 5, Module 7; Figure 6, Cluster #4). ZAT genes are induced in abiotic stress conditions and are required for cytosolic ascorbate peroxidase expression during oxidative stress in Arabidopsis (Rizhsky et al., 2004). These regulatory genes may be involved in the oxidative burst triggered by the pathogen or in oxidation-reduction processes to protect the plant.
Infection by B. cinerea enhanced the expression of several TFs involved in the ET response (EIN3, ERF5, and ERF13; Figure 5, Module 1; Figure 6, Cluster #3). ERF5 was slightly induced by the fungal infection, most likely in RR fruits. Among the TFs associated with the ET response, several TFs belong to the ERF family, have been associated with stress responses (McGrath et al., 2005) and as modulators of the SA/JA crosstalk (Van der Does et al., 2013). Interestingly, in infected leaves, ERF1 (Ethylene Response Factor 1; Figure 6, Cluster #2) was significantly induced in N-sufficient conditions and repressed in low N conditions. ERF1 has been implicated in JA and ET signaling, and pathogen defense (Berrocal-Lobo et al., 2002; Lorenzo et al., 2003). In addition, the Jasmonate-ZIM-Domain 1 gene (JAZ1; Figure 5, Module 5; Figure 6, Cluster #3) was also induced by the fungus infection in leaves under all analyzed N-conditions. This gene is involved in JA-signaling and, recently, its role in plant defense has been described (Li et al., 2014). These different patterns of expression of TFs related with hormone responses, suggest the existence of several regulatory pathways interacting in a complex and highly dynamic manner, linking plant N status, and the plant defense response.
Although, different agronomic reports indicate that the amount of N fertilization can modulate disease severity caused by pathogens in numerous plant species (Snoeijers et al., 2000; Walters and Bingham, 2007; Dordas, 2008; Fagard et al., 2014), the underlying mechanisms are poorly characterized. In general terms, high N supply decreases the severity of the infections from necrotrophic pathogens, and N-limiting conditions decrease plant susceptibility to biotrophy or hemibiotrophic microorganisms (Walters and Bingham, 2007; Dordas, 2008). Despite these general trends, several reports show contradictory evidence, suggesting a complex relationship linking N metabolism and disease infection processes. In this report, we investigated the role of plant N nutrition in the outcome of the plant-pathogen interaction, employing two important agronomic models: B. cinerea and S. lycopersicum.
We showed that tomato plants grown under contrasting nitrate conditions exhibit differential susceptibility to fungal infection. We used nitrate as a nitrogen source, since it is the most important N source for plants in agricultural soils (von Wirén et al., 2000) and it plays an important role as the major nitrogen resource in tomato plants (Wang et al., 2001; Fu et al., 2015). In agreement with previous studies, tomato plants grown under N-sufficient or higher conditions showed a significantly decreased susceptibility to B. cinerea (Hoffland et al., 1999; Lecompte et al., 2010). We extended previous results to include fruits and showed that there is reduced susceptibility to B. cinerea in fruits at two developmental stages (MG and RR) from plants grown under N-sufficient conditions, suggesting that the nitrogen/disease connection might extend beyond vegetative tissues.
The induction of the plant defense response is an energetically costly process (Berger et al., 2007). As such, different studies have suggested that the major role of the plant's primary metabolism during the plant-pathogen interaction is to support the increased energy requirements brought about by the defense response (Bolton, 2009; Kangasjärvi et al., 2012), and some pathogens in early steps of infections trigger a nutrient-limiting environment (Mathioni et al., 2011). In this scenario, it is reasonable to hypothesize that “well nourished” plants may defend better, as observed in tomato grown under N-sufficient conditions. This observation sharply differs from crop plants that present less susceptibility under low nitrogen availability (Davidson et al., 2007; Lecompte et al., 2013). Notably, as observed from our global gene expression profiling, plant metabolic processes were significantly affected by B. cinerea, irrespective of the nitrogen condition used to grow the plants. This finding is consistent with previous reports showing that the plant-pathogen interaction regulates plant metabolism (Ward et al., 2010; Rojas et al., 2014). In order to establish a favorable energy balance for defense, it has been suggested that the up-regulation of defense-related genes is compensated by the down-regulation of genes involved in other metabolic pathways (Berger et al., 2007; Massad et al., 2012). In agreement with what has been reported for B. cinerea infections in A. thaliana (AbuQamar et al., 2006; Windram et al., 2012), several genes involved in photosynthesis were down regulated in our tomato experiments in response to the fungus infection, under all N conditions.
In the context of plant-pathogens interactions, the nutrition of the plant can have an effect on secondary metabolite production. Plant N-status modifies amino acid accumulation involved in the biosynthesis of defense-associated secondary metabolites, such as flavonoids and phenylpropanoids (Ward et al., 2010; Fagard et al., 2014). Phenylpropanoid pathways regulate the resistance against B. cinerea in the sitiens tomato mutant (Seifi et al., 2013). Accordingly, we found several genes involved in this pathway up-regulated by the fungus infections in N-sufficient conditions. In addition, the amino acid content and its relative concentration are modified during the plant disease process, suggesting that amino acid metabolism can impact plant-pathogen interactions. Recently, reports show that genes involved in amino acid metabolism play a role in plant defense responses to pathogens (Hwang et al., 2011; Seifi et al., 2014). N-metabolism and amino acid homeostasis can modulate the plant redox status affecting the plant disease response (Liu et al., 2010). These results suggest that N metabolism could modulate -at least in part- plant defense responses affecting cellular redox status as well as secondary metabolite pathways.
The relevance of N and amino acid metabolism has been highlighted before in the context of plant-pathogen interactions (Hwang et al., 2011). For instance, several studies have shown that N-related gene expression is altered by pathogens (Pageau et al., 2006; Masclaux-Daubresse and Daniel-Vedele, 2010). Interestingly, the Arabidopsis lht1 (lysine histidine transport 1) mutant is more resistant to a large spectrum of pathogens, an observation that is associated with an altered redox status (Liu et al., 2010). In this regard, we found several genes associated with redox status, N-metabolism and transporters affected by B. cinerea infections, depending on plant N-conditions. NRT2, a putative high-affinity nitrate transporter, was induced only in N-sufficient and higher conditions. Interestingly, in Arabidopsis, two NRT2s have been reported as involved in plant defense (Camañes et al., 2012; Dechorgnat et al., 2012) and the Arabidopsis nrt2.6 mutant exhibits higher sensitivity to the necrotrophic bacterium Erwina amylovora. The susceptibility of this mutant appears to be related to reduced ROS accumulation (Dechorgnat et al., 2012), again linking redox status, N metabolism and defense. Similarly, we found that oxidation-reduction processes were over-represented among the biological functions responding to the fungal infection, under N-sufficient and higher conditions. Recently, nitrate reallocation in plants has been proposed as a regulator of the trade-off between plant growth and environmental adaptation (Zhang et al., 2014). All these results suggest that nitrate metabolism and transport have an impact on the plant's defense response, and that the nitrogen-defense connection goes beyond a direct metabolic relationship. Whether, this can be explained by N-derived alterations in the redox status of the plant requires further investigation.
Different plant hormones may provide an additional layer of regulation, underlying—at least in part—the complex relationship between nitrogen metabolism and plant defense response. It has been described that plant hormones play a role in the defense against B. cinerea (AbuQamar et al., 2006; Windram et al., 2012). Specifically, genetic studies on Arabidopsis and tomato indicate that JA and ET are key regulators of this defense response. Treatments with ET in tomato plants increased resistance to this fungus in leaves, while mutants in JA biosynthesis showed increased susceptibility to B. cinerea in leaves (Díaz et al., 2002; AbuQamar et al., 2008). Consistent with these prior findings, genes involved in ET response were identified by the GO term analysis in our network of B. cinerea-responsive genes associated with the plant N status. Likewise, several TFs associated with the JA/ET response were identified, including ERF1 and ERF5. The latter acts as a positive regulator of JA/ET-responsive defense genes against B. cinerea in Arabidopsis (Berrocal-Lobo et al., 2002). Notably, we found that ERF1 was induced in response to B. cinerea under N-sufficient conditions, while repressed in plants grown under N-limited ones. This TF has been described as a key regulator of JA/ET signaling, activating plant defense responses against B. cinerea in leaves (Berrocal-Lobo et al., 2002; Lorenzo et al., 2003) and necrotrophic Rhizopus nigricans in RR tomato fruit (Pan et al., 2013). Moreover, ERF1 gene expression increases in response to B. cinerea infections in MG tomato fruit but is reduced in infected RR fruit (Blanco-Ulate et al., 2013). MG and RR tomato fruits show different susceptibility responses against this fungus infection, each associated with specific expression profiles of genes involved in ET and others hormone biosynthesis (Blanco-Ulate et al., 2013). More studies are necessary to assess the regulation of this TF in response to infections in different tissues.
Although, the plant hormone SA has been associated with the resistance to biotrophic pathogens, its role in plant defense against B. cinerea is not completely clear. On the one hand, the exogenous application of SA has been reported to result in decreased susceptibility to this fungus in tomato (Audenaert et al., 2002) and Arabidopsis (Ferrari et al., 2003), suggesting that SA affects the resistance against this fungus. On the other hand, the stearoyl-ACP desaturase (ssi2) Arabidopsis mutant, which exhibits increased SA levels and higher expression of the PR1 gene, displayed, increased susceptibility to B. cinerea (Kachroo et al., 2001; Nandi et al., 2005). High resolution temporal transcriptomic analysis of the Arabidopsis defense response against this fungus, revealed complex regulatory networks, with different timings for JA, ET, and SA signaling (Windram et al., 2012) as well as the accumulation of plant phytoalexins and other defense-related proteins and molecules (Van Baarlen et al., 2004; AbuQamar et al., 2006; Scalschi et al., 2015). Our data support a role for these signaling pathways connecting N status and infection. For instance, the silencing of the OPR3 gene enhances susceptibility to B. cinerea by affecting the JA biosynthesis pathway (Scalschi et al., 2015), and we found this gene to be repressed under N-limiting conditions. The WRKY70 TF, which is up-regulated in response to B. cinerea infection, plays a role in SA-JA crosstalk, activating SA-induced genes and repressing JA-responsive genes (Li et al., 2004). Interestingly, the wrky70 mutant is more susceptible to this fungus (AbuQamar et al., 2006), since notably, a gene encoding for a putative WRKY70 TF in tomato is the most repressed gene in response to B. cinerea, as shown here, but only under low N availability, a condition in which higher susceptibility to the fungus was observed.
Taken together, the findings reported here suggest that N metabolism affects different mechanisms, including redox status, secondary metabolites and plant hormones (JA and ET), the latter of which globally alter the expression of defense-related genes. In sum, these mechanisms together modulate disease susceptibility and hence the outcome of the plant-pathogen interaction under different N conditions.
Materials and Methods
Growth Conditions of Tomato Plants
Solanum lycopersicum cv. MicroTom plants were cultivated under different N growth conditions employing 16:8 h light: dark cycles (150 micromoles/m2/s) at 20°C. Plants seeds were spread in plastic pots containing sterilized vermiculite, previously equilibrated with Murashige & Skoog (MS) modified basal salt mixture without nitrogen (PhytoTechnology Laboratories). This nutrient solution was supplemented with different concentrations of nitrate. Six mM nitrate was found to give maximal growth and thus considered as N-sufficient condition, while 2 mM and 4 mM nitrate were used as severe and mild N-limiting conditions, respectively. To evaluate a higher N input, tomato plants were grown using 12 mM nitrate. Four nutrition solutions—each containing the mentioned nitrate concentrations—were prepared and added to each pot every 2 days, employing the same volume (50 ml). Shoot biomass (dry weight) and N concentration (Dumas method) were determined from an average of 8 plants after 4 weeks of growth. Fruits were tagged at 2-day post anthesis (dpa) and harvested at 32 and 40 dpa as mature green and ripening red (MG and RR, respectively). A color chart was employed to confirm ripening stages.
Botrytis cinerea Growth Conditions and Virulence Assays
The necrotrophic fungus used in this study was B. cinerea strain B05.10. This strain was routinely cultivated in Petri dishes containing potato dextrose agar (PDA, AppliChem) with 10% homogenized bean leaves. Conidia employed for virulence assays were obtained from an agar plug followed by glass-wool filtration using Gamborg B5-2% glucose (Duchefa Biochemie) medium.
For virulence assays, 4-week-old tomato plants grown under different nitrate conditions were inoculated in the third and fourth fully expanded leaves in planta as described (Cantu et al., 2008; Canessa et al., 2013). Briefly, conidia were suspended in Gamborg B5-2% glucose medium and adjusted to a final concentration of 5 × 105 conidia/ml in the same medium supplemented with 10 mM KH2PO4/K2HPO4, pH 6.4. Conidial suspensions (10 μl) were used to inoculate leaves of six tomato plants obtained from each N condition analyzed. All plants were incubated inside plastic boxes at 20°C employing the growth conditions mentioned above, under humid environment for the indicated periods of time. Lesions were measured using the ImageJ software employing an external calibration scale (Schneider et al., 2012). On the other hand, fruits were inoculated at the day of harvest, as described (Cantu et al., 2008; Canessa et al., 2013). Briefly, fruits were first disinfected using 10% (v/v) bleach, followed by three consecutive water rinses. Thereafter, fruits were carefully punctured (2 mm depth, 1 mm diameter) employing a sterile needle at five sites per tomato. Subsequently, fruits were inoculated with 10 μl of conidia suspension (5 × 105 conidia/ml). Tomatoes used as mock material were inoculated with the same volume of Gamborg B5-2% glucose-10 mM KH2PO4/K2HPO4 solution. All fruit samples were incubated at 20°C in high humidity, during the indicated periods of time (3, 4, and 5 dpi). At the end of the infection period, tomato lesions were measured as mentioned above. The evaluation of tomato susceptibility to B. cinerea was independently performed three times employing at least 4 fruits for each experimental condition.
For quantification of fungal development on plant tissue, DNA from leaves or fruits for each experimental condition was isolated following a method for genomic DNA described by Edwards et al. (1991). For real time PCR quantification, 25 μl samples were prepared containing 10 μl of DNA solutions and 300 nM of appropriate primers. Real-time PCR reactions were run in triplicates using the Brilliant SYBR Green QPCR Reagents on a StepOnePlus™ Real-Time PCR System (Life Technologies), as described bellow in RT-qPCR section. The sequences of the primers used in this study are detailed in Table S6. Serial dilutions of pure genomic DNA from each species were used to develop a calibration curve, as described (Gachon and Saindrenan, 2004).
Total RNA was extracted from both mock-treated and fungus-infected tomato leaves and fruits, following a CTAB-spermidine extraction protocol (Reid et al., 2006). In the case of leaves, a total of 12 leaves per N treatment were pulled together. This procedure was repeated three times (n = 3) and further used in microarray analysis and RT-qPCR experiments (see below). For fruits, total RNA was prepared from the combined fruit outer pericarp and epidermis, within approximately 0.25 cm radius around the B. cinerea inoculation sites. Each biological replicate (n = 3) consisted of an independent pool of samples from four different fruits. RNA samples were further quantified and analyzed by microfluidic analysis employing the Agilent Technologies' 2100 Bioanalyzer, following the manufacturer's instructions.
For microarray hybridizations, three mock-inoculated and three fungus-infected biological replicates from tomato plants grown under 2, 6, and 12 mM nitrate were selected for global gene expression analysis. Five hundred nanograms of total RNA was processed for microarray hybridization using the GeneChip one-cycle target-labeling kit (Affymetrix), according to the manufacturer's instructions. The fragmented aRNA (7.5 μg) was hybridized on a GeneChip Tomato genome array (containing 10038 tomato probe sets for more than 9200 tomato genes) using standard procedures (45°C for 16 h). The arrays were washed and stained in a Fluidics Station 450 (Affymetrix). Array scanning was carried out with the GeneChip scanner 300 and image analysis was performed using the GeneChip Operating Software. Thereafter, GeneChip array data were quality assessed using a set of standard quality control steps described in the Affymetrix manual “GeneChip Expression Analysis: Data Analysis Fundamentals.” Array data were processed and normalized with RMA (Robust Multi-Array Average) (Irizarry et al., 2003) using the affy R package (Gautier et al., 2004). Approximately 60–75% of probe sets were significantly detected in all microarray hybridizations. To evaluate array reproducibility, spearman rank coefficients were computed and ranged between 0.97 and 0.99. The raw data for all hybridizations can be found in NCBIs Gene Expression Omnibus (Edgar et al., 2002) and are accessible through GEO Series accession number GSE73006.
To determine DEG in response to B. cinerea infection affected by the N-conditions, normalized data were subjected to a Two-way analysis of variance (ANOVA) with adjusted p < 0.01 by Benjamini and Hochberg false discovery rate correction (Benjamini and Hochberg, 1995). For ANOVA analysis, we used a model with expression Y of a given gene i calculated as Yi= β0 + β1B + β2N + β3N/B + ε, where β0 is the global mean, and where β1, β2, and β3 are the effects of the B. cinerea infection, the N-condition and the interaction between these two factors (N/B), respectively (Table S1). The variable ε corresponds to the unexplained variance. In addition, to identify DEG between two conditions (e.g., infected and control samples at each N-condition) the Rank-Product method was used (p < 0.05) (Breitling et al., 2004; see Tables S2–S4).
Functional Annotation of Differentially Expressed Genes (DEG)
For the assignment of functional annotations and categorizations for DEG, we employed the easy to use web tool Generic Gene Ontology (GO) Term Finder (Boyle et al., 2004). For this purpose, we constructed a custom-made GO database employing InterProScan 5 (Jones et al., 2014) and all predicted protein coding genes in the S. lycopersicum genome database (gene models ITAG release 2.4, https://solgenomics.net).
Gene Network Analysis
Genes possessing a significant Nitrogen/B. cinerea (N/B) interaction factor were selected to carry out a gene network construction (Table S1). First, we perform an analysis of ortholog group assignment between tomato and A. thaliana. For this purpose, tomato sequences (gene models ITAG release 2.4), matching to probe sets on the array were aligned to the Arabidopsis genome, using OrthoMCL (Li et al., 2003) orthologs genes were assigned. Thereafter, with all Arabidopsis loci identifiers, a gene co-expression network was constructed using 11,171 microarray experiments obtained from the ATTED-II database (Obayashi et al., 2014), accessed in March 2015. Co-expression values were calculated using weighted Pearson's correlation coefficient, as described in Obayashi et al. (2014). This database collects gene expression data in Arabidopsis from a wide range of microarray experiments. We included protein–DNA interactions (Weirauch et al., 2014), considering at least one TF binding site in the upstream gene region (1000 bp) and over-representation of the TF binding site (two standard deviations) above the mean occurrence in all the upstream sequences in the genome. To improve the regulatory interaction predictions, the protein–DNA interactions were filtered to include only transcription factor/target pairs for which co-expression values are highly (cut-off above the 90th percentile considering absolute value of gene co-expression) (Obayashi et al., 2014) and significantly correlated (p < 0.05) in our microarray experiment. The resulting network was visualized using CYTOSCAPE (Data Sheet S1; Shannon et al., 2003).
RT-qPCR and Clustering Analysis of Gene Expression Data
cDNA synthesis was carried out using the Improm-II reverse transcriptase according to the manufacturer's instructions (Promega). RT-qPCR was carried out using the Brilliant SYBR Green QPCR Reagents on a StepOnePlus™ Real-Time PCR System (Life Technologies). For mRNA levels normalization, the PROTEIN PHOSPHATASE 2A catalytic subunit (PP2Acs, SGN-U567355) was used (Løvdal and Lillo, 2009; Dekkers et al., 2012). The expression of this reference gene was stable in our microarray data (CV = 2%) and in microarrays in response to B. cinerea in tomato fruit at two developmental stages (MG and RR, CV = 3% Cantu et al., 2009).
Total RT-qPCR reaction volume was 25 μl, containing 2.0 μl of cDNA template and 140 nM of each primer. The reactions were performed under the following conditions: 95°C for 10 min, followed by 40 cycles of 95°C, 30 s; 55°C, 30 s; and 72°C, 30 s, followed by a melting curve analysis from 55 to 95°C. Fluorescence values were acquired during the annealing period of the RT-qPCR procedure. All experiments were performed with three biological replicates, with their corresponding three technical replicates. To determine the mRNA levels in response to the B. cinerea infection, infected leaves and fruits from plants grown under 2, 4, 6, and 12 mM nitrate were compared to their corresponding mock-treated control grown under same culture conditions. Genes showing a similar expression pattern were analyzed and visualized using the average-linkage hierarchical clustering performed in Cluster 2.11 software, as described (Eisen et al., 1998).
PC, GH, IR, AV performed and analyzed virulence assays; GH, AV performed and analyzed microarray experiments; PC and AV conducted the GO analysis of gene expression data; TM, JC, AV conducted bioinformatics analysis for gene network construction; IR, GH, JR, AV, PC completed RT-qPCR analysis; AV, PC, and RG analyzed all data and wrote the paper.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Research was funded by the Fondo Nacional de Desarrollo Científico y Tecnológico FONDECYT-Chile 11110095 and 11140678, Millennium Nucleus Center for Plant Systems and Synthetic Biology (NC130030) and Millennium Nucleus for Fungal Integrative and Synthetic Biology (NC120043).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2015.00911
Figure S1. Tomato plants grown under contrasting N conditions. (A) Shoots biomass (average) of MicroTom tomato plants (4-weeks-old) grown under contrasting N conditions. (B) N percentage of plants, determined as dry weight (see Materials and Methods; n = 4). Different letters indicate significant differences among treatments (p ≤ 0.05; error bars indicate SEM).
Figure S2. Fungal growth quantification over a B. cinerea infection time-course. The abundance of the fungus was quantified by qPCR amplification of cutA (see Materials and Methods) in infected leaves samples (A), and MG (B) and RR (C) fruits. Different letters indicate significant differences among treatments (p ≤ 0.05; error bars indicate SEM; n = 3).
Figure S3. ANOVA analysis of microarray gene expression data. The number of genes with significant factors obtained by Two-way ANOVA analysis of global expression data is represented (p < 0.01), employing the SUNGEAR tool (Poultney et al., 2007). The triangle and its vertices represent the analyzed factors: B. cinerea (B), N conditions (N) and the interaction between these two factors (N/B). Circles (and their respective size) represent the number of genes controlled by the different factors, as indicated by the arrows around the circles.
Table S1. Microarray expression data for genes with N/B ANOVA model (adjusted p < 0.01). The table includes the Affymetrix probes, the tomato accessions (ITAG2.4, Sol Genomics Network, https://solgenomics.net) and their annotations, Arabidopsis (TAIR, http://arabidopsis.org) or GeneBank accessions, the p-value adjusted by FDR and the log2-fold changes of the comparison between B. cinerea infected and mock-treated sample (B/M) in each N-conditions (2, 6, and 12 mM nitrate).
Table S2. Microarray expression data for genes that are either up- or down-regulated by the B. cinerea infection in comparison to mock-inoculated plants grown in N-limiting conditions (2 mM nitrate) (ANOVA; Rank-products, p < 0.05). The table includes the Affymetrix probes, the tomato accessions (ITAG2.4, Sol Genomics Network, https://solgenomics.net) and their annotations, Arabidopsis (TAIR, http://arabidopsis.org) or GeneBank accessions, the p-value, Gene Ontology and the log2-fold changes of the comparison between B. cinerea infected and mock-treated sample (B/M) grown in N-limiting conditions (2 mM nitrate).
Table S3. Microarray expression data for genes that are either up- or down-regulated by the B. cinerea infection in comparison to mock-inoculated plants grown in N- sufficient conditions (6 mM nitrate) (ANOVA; Rank-products, p < 0.05). The table includes the Affymetrix probes, the tomato accessions (ITAG2.4, Sol Genomics Network, https://solgenomics.net) and their annotations, Arabidopsis (TAIR, http://arabidopsis.org) or GeneBank accessions, the p-value, Gene Ontology and the log2-fold changes of the comparison between B. cinerea infected and mock-treated sample (B/M) grown in N-sufficient conditions (6 mM nitrate).
Table S4. Microarray expression data for genes that are either up- or down-regulated by the B. cinerea infection in comparison to mock-inoculated plants grown in N- sufficient conditions (12 mM nitrate) (ANOVA; Rank-products, p < 0.05). The table includes the Affymetrix probes, the tomato accessions (ITAG2.4, Sol Genomics Network, https://solgenomics.net) and their annotations, Arabidopsis (TAIR, http://arabidopsis.org) or GeneBank accessions, the p-value, Gene Ontology, and the log2-fold changes of the comparison between B. cinerea infected and mock-treated sample (B/M) grown in N-limiting sufficient (12 mM nitrate).
Table S5. Candidate TF used for RT-qPCR analysis. The table includes the Arabidopsis (TAIR, http://arabidopsis.org) and the tomato accessions (Sol Genomics Network, https://solgenomics.net), the gene names, the module in network analysis and changes in relative expression (log2) between infected and mock-treated samples in leaves and fruits (MG and RR) from plants grown under different N-conditions (2 mM, 4 mM, 6 mM, and 12 mM nitrate).
Table S6. Primer sequences used for qPCR.
Data Sheet S1. Gene network in cytoscape.
AbuQamar, S., Chai, M.-F., Luo, H., Song, F., and Mengiste, T. (2008). Tomato protein kinase 1b mediates signaling of plant responses to necrotrophic fungi and insect herbivory. Plant Cell 20, 1964–1983. doi: 10.1105/tpc.108.059477
AbuQamar, S., Chen, X., Dhawan, R., Bluhm, B., Salmeron, J., Lam, S., et al. (2006). Expression profiling and mutant analysis reveals complex regulatory networks involved in Arabidopsis response to Botrytis infection. Plant J. 48, 28–44. doi: 10.1111/j.1365-313X.2006.02849.x
Alba, R., Payton, P., Fei, Z., McQuinn, R., Debbie, P., Martin, G. B., et al. (2005). Transcriptome and selected metabolite analyses reveal multiple points of ethylene control during tomato fruit development. Plant Cell 17, 2954–2965. doi: 10.1105/tpc.105.036053
Audenaert, K., Pattery, T., Cornelis, P., and Höfte, M. (2002). Induction of systemic resistance to Botrytis cinerea in tomato by Pseudomonas aeruginosa 7NSK2: role of salicylic acid, pyochelin, and pyocyanin. Mol. Plant Microbe Interact. 15, 1147–1156. doi: 10.1094/MPMI.2002.15.11.1147
Berrocal-Lobo, M., Molina, A., and Solano, R. (2002). Constitutive expression of ETHYLENE-RESPONSE-FACTOR1 in Arabidopsis confers resistance to several necrotrophic fungi. Plant J. 29, 23–32. doi: 10.1046/j.1365-313x.2002.01191.x
Blanco-Ulate, B., Vincenti, E., Powell, A. L. T., and Cantu, D. (2013). Tomato transcriptome and mutant analyses suggest a role for plant stress hormones in the interaction between fruit and Botrytis cinerea. Front. Plant Sci. 4:142. doi: 10.3389/fpls.2013.00142
Boyle, E. I., Weng, S., Gollub, J., Jin, H., and Botstein, D. (2004). GO:: TermFinder—open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics 20, 3710–3715. doi: 10.1093/bioinformatics/bth456
Breitling, R., Armengaud, P., Amtmann, A., and Herzyk, P. (2004). Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Lett. 573, 83–92. doi: 10.1016/j.febslet.2004.07.055
Camañes, G., Pastor, V., Cerezo, M., García-Andrade, J., Vicedo, B., García-Agustín, P., et al. (2012). A deletion in NRT2.1 attenuates Pseudomonas syringae-induced hormonal perturbation, resulting in primed plant defenses. Plant Physiol. 158, 1054–1066. doi: 10.1104/pp.111.184424
Canessa, P., Schumacher, J., Hevia, M. A., Tudzynski, P., and Larrondo, L. F. (2013). Assessing the effects of light on differentiation and virulence of the plant pathogen Botrytis cinerea: characterization of the White Collar Complex. PLoS ONE 8:e84223. doi: 10.1371/journal.pone.0084223
Cantu, D., Blanco-Ulate, B., Yang, L., Labavitch, J. M., Bennett, A. B., and Powell, A. L. T. (2009). Ripening-regulated susceptibility of tomato fruit to botrytis cinerea requires NOR but not RIN or Ethylene. Plant Physiol. 150, 1434–1449. doi: 10.1104/pp.109.138701
Cantu, D., Vicente, A. R., Greve, L. C., Dewey, F. M., Bennett, A. B., Labavitch, J. M., et al. (2008). The intersection between cell wall disassembly, ripening, and fruit susceptibility to Botrytis cinerea. Proc. Natl. Acad. Sci. U.S.A. 105, 859–864. doi: 10.1073/pnas.0709813105
Castro Marín, I., Loef, I., Bartetzko, L., Searle, I., Coupland, G., Stitt, M., et al. (2010). Nitrate regulates floral induction in Arabidopsis, acting independently of light, gibberellin and autonomous pathways. Planta 233, 539–552. doi: 10.1007/s00425-010-1316-5
Catinot, J., Huang, J. B., and Huang, P. Y. (2015). ETHYLENE RESPONSE FACTOR 96 positively regulates Arabidopsis resistance to necrotrophic pathogens by direct binding to GCC elements of jasmonate–and …. Plant. doi: 10.1111/pce.12583
Davidson, J. A., Pande, S., Bretag, T. W., Lindbeck, K. D., and Krishna-Kishore, G. (2007). “Biology and management of botrytis spp. in legume crops,” in Botrytis: Biology, Pathology and Control, eds Y. Elad, B. Williamson, P. Tudzynski, and N. Delen (Dordrecht: Springer), 295–318. doi: 10.1007/978-1-4020-2626-3_16
Dean, R., van Kan, J. A. L., Pretorius, Z. A., Hammond-Kosack, K. E., Di Pietro, A., Spanu, P. D., et al. (2012). The Top 10 fungal pathogens in molecular plant pathology. Mol. Plant Pathol. 13, 414–430. doi: 10.1111/J.1364-3703.2011.00783.X
Dechorgnat, J., Patrit, O., Krapp, A., Fagard, M., and Daniel-Vedele, F. (2012). Characterization of the Nrt2.6 gene in Arabidopsis thaliana: a link with plant response to biotic and abiotic stress. PLoS ONE 7:e42491. doi: 10.1371/journal.pone.0042491
Dekkers, B. J. W., Willems, L., Bassel, G. W., van Bolderen-Veldkamp, R. P. M., Ligterink, W., Hilhorst, H. W. M., et al. (2012). Identification of reference genes for RT-qPCR expression analysis in Arabidopsis and tomato seeds. Plant Cell Physiol. 53, 28–37. doi: 10.1093/pcp/pcr113
Dietrich, R., Plo, ß. K., and Heil, M. (2004). Constitutive and induced resistance to pathogens in Arabidopsis thaliana depends on nitrogen supply. Plant Cell Environ. 27, 896–906. doi: 10.1111/j.1365-3040.2004.01195.x
Farmer, E. E., Alméras, E., and Krishnamurthy, V. (2003). Jasmonates and related oxylipins in plant responses to pathogenesis and herbivory. Curr. Opin. Plant Biol. 6, 372–378. doi: 10.1016/S1369-5266(03)00045-1
Ferrari, S., Plotnikova, J. M., and De Lorenzo, G. (2003). Arabidopsis local resistance to Botrytis cinerea involves salicylic acid and camalexin and requires EDS4 and PAD2, but not SID2, EDS5 or PAD4. Plant J. 35, 193–205. doi: 10.1046/j.1365-313X.2003.01794.x
Gachon, C., and Saindrenan, P. (2004). Real-time PCR monitoring of fungal development in Arabidopsis thaliana infected by Alternaria brassicicola and Botrytis cinerea. Plant Physiol. Biochem. 42, 367–371. doi: 10.1016/j.plaphy.2004.04.001
Good, A. G., Shrawat, A. K., and Muench, D. G. (2004). Can less yield more? Is reducing nutrient input into the environment compatible with maintaining crop production? Trends Plant Sci. 9, 597–605. doi: 10.1016/j.tplants.2004.10.008
Hoffland, E., Jeger, M. J., and van Beusichem, M. L. (2000). Effect of nitrogen supply rate on disease resistance in tomato depends on the pathogen. Plant Soil 218, 239–247. doi: 10.1023/A:1014960507981
Hwang, I. S., An, S. H., and Hwang, B. K. (2011). Pepper asparagine synthetase 1 (CaAS1) is required for plant nitrogen assimilation and defense responses to microbial pathogens. Plant J. 67, 749–762. doi: 10.1111/j.1365-313X.2011.04622.x
Irizarry, R. A., Hobbs, B., Collin, F., Beazer-Barclay, Y. D., Antonellis, K. J., Scherf, U., et al. (2003). Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostat. 4, 249–264. doi: 10.1093/biostatistics/4.2.249
Jones, P., Binns, D., Chang, H.-Y., Fraser, M., Li, W., McAnulla, C., et al. (2014). InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236–1240. doi: 10.1093/bioinformatics/btu031
Kachroo, P., Shanklin, J., Shah, J., Whittle, E. J., and Klessig, D. F. (2001). A fatty acid desaturase modulates the activation of defense signaling pathways in plants. Proc. Natl. Acad. Sci. U.S.A. 98, 9448–9453. doi: 10.1073/pnas.151258398
Kangasjärvi, S., Neukermans, J., Li, S., Aro, E.-M., and Noctor, G. (2012). Photosynthesis, photorespiration, and light signalling in defence responses. J. Exp. Bot. 63, 1619–1636. doi: 10.1093/jxb/err402
Kliebenstein, D. J., Rowe, H. C., and Denby, K. J. (2005). Secondary metabolites influence Arabidopsis/Botrytis interactions: variation in host production and pathogen sensitivity. Plant J. 44, 25–36. doi: 10.1111/j.1365-313X.2005.02508.x
Lecompte, F., Abro, M. A., and Nicot, P. C. (2010). Contrasted responses of Botrytis cinerea isolates developing on tomato plants grown under different nitrogen nutrition regimes. Plant Pathol. 59, 891–899. doi: 10.1111/j.1365-3059.2010.02320.x
Lecompte, F., Abro, M. A., and Nicot, P. C. (2013). Can plant sugars mediate the effect of nitrogen fertilization on lettuce susceptibility to two necrotrophic pathogens: Botrytis cinerea and Sclerotinia sclerotiorum? Plant Soil 369, 387–401. doi: 10.1007/s11104-012-1577-9
Li, C., He, X., Luo, X., Xu, L., Liu, L., Min, L., et al. (2014). Cotton WRKY1 mediates the plant defense-to-development transition during infection of cotton by Verticillium dahliae by activating JASMONATE ZIM-DOMAIN1 expression. Plant Physiol. 166, 2179–2194. doi: 10.1104/pp.114.246694
Li, J., Brader, G., and Palva, E. T. (2004). The WRKY70 transcription factor: a node of convergence for jasmonate-mediated and salicylate-mediated signals in plant defense. Plant Cell 16, 319–331. doi: 10.1105/tpc.016980
Linquist, B. A., Byous, E., Jones, G., Williams, J. F., and Six, J. (2008). Nitrogen and potassium fertility impacts on aggregate sheath spot disease and yields of rice. Plant Prod. Sci. 11, 260–267. doi: 10.1626/pps.11.260
Liu, G., Ji, Y., Bhuiyan, N. H., Pilot, G., Selvaraj, G., Zou, J., et al. (2010). Amino acid homeostasis modulates salicylic acid-associated redox status and defense responses in Arabidopsis. Plant Cell 22, 3845–3863. doi: 10.1105/tpc.110.079392
Long, D. H., Lee, F. N., and TeBeest, D. O. (2007). Effect of nitrogen fertilization on disease progress of rice blast on susceptible and resistant cultivars. Plant Dis. 84, 403–409. doi: 10.1094/PDIS.2000.84.4.403
Lorenzo, O., Piqueras, R., Sánchez-Serrano, J. J., and Solano, R. (2003). ETHYLENE RESPONSE FACTOR1 integrates signals from ethylene and jasmonate pathways in plant defense. Plant Cell 15, 165–178. doi: 10.1105/tpc.007468
Løvdal, T., and Lillo, C. (2009). Reference gene selection for quantitative real-time PCR normalization in tomato subjected to nitrogen, cold, and light stress. Anal. Biochem. 387, 238–242. doi: 10.1016/j.ab.2009.01.024
Masclaux-Daubresse, C., and Daniel-Vedele, F. (2010). Nitrogen uptake, assimilation and remobilization in plants: challenges for sustainable and productive agriculture. Ann Bot. 105, 1141–1157. doi: 10.1093/aob/mcq028
Massad, T. J., Dyer, L. A., and Vega, C. G. (2012). Costs of defense and a test of the carbon-nutrient balance and growth-differentiation balance hypotheses for two co-occurring classes of plant defense. PLoS ONE 7:e47554. doi: 10.1371/journal.pone.0047554
Mathioni, S. M., Bel, ó. A., Rizzo, C. J., Dean, R. A., and Donofrio, N. M. (2011). Transcriptome profiling of the rice blast fungus during invasive plant infection and in vitro stresses. BMC Genomics 12:049. doi: 10.1186/1471-2164-12-49
McGrath, K. C., Dombrecht, B., Manners, J. M., Schenk, P. M., Edgar, C. I., Maclean, D. J., et al. (2005). Repressor- and activator-type ethylene response factors functioning in jasmonate signaling and disease resistance identified via a genome-wide screen of Arabidopsis transcription factor gene expression. Plant Physiol. 139, 949–959. doi: 10.1104/pp.105.068544
Morris, J. H., Apeltsin, L., Newman, A. M., Baumbach, J., Wittkop, T., Su, G., et al. (2011). clusterMaker: a multi-algorithm clustering plugin for Cytoscape. BMC Bioinformatics 12:436. doi: 10.1186/1471-2105-12-436
von Wirén, N., Lauter, F. R., Ninnemann, O., Gillissen, B., Walch-Liu, P., Engels, C., et al. (2000). Differential regulation of three functional ammonium transporter genes by nitrogen in root hairs and by light in leaves of tomato. Plant J. 21, 167–175. doi: 10.1046/j.1365-313x.2000.00665.x
Obayashi, T., Okamura, Y., Ito, S., Tadaka, S., Aoki, Y., Shirota, M., et al. (2014). ATTED-II in 2014: evaluation of gene coexpression in agriculturally important plants. Plant Cell Physiol. 55, e6(1–7). doi: 10.1093/pcp/pct178
Olesen, J. E., Jørgensen, L. N., Petersen, J., and Mortensen, J. V. (2003). Effects of rates and timing of nitrogen fertilizer on disease control by fungicides in winter wheat. 2. Crop growth and disease development. J. Agr. Sci. 140, 15–29. doi: 10.1017/S0021859602002897
Pageau, K., Reisdorf-Cren, M., Morot-Gaudry, J.-F., and Masclaux-Daubresse, C. (2006). The two senescence-related markers, GS1 (cytosolic glutamine synthetase) and GDH (glutamate dehydrogenase), involved in nitrogen mobilization, are differentially regulated during pathogen attack and by stress hormones and reactive oxygen species in Nicotiana tabacum L. leaves. J. Exp. Bot. 57, 547–557. doi: 10.1093/jxb/erj035
Pan, X., Zhu, B., Luo, Y., and Fu, D. (2013). Unraveling the protein network of tomato fruit in response to necrotrophic phytopathogenic Rhizopus nigricans. PLoS ONE 8:e73034. doi: 10.1371/journal.pone.0073034
Poultney, C. S., Gutiérrez, R. A., Katari, M. S., Gifford, M. L., Paley, W. B., Coruzzi, G. M., et al. (2007). Sungear: interactive visualization and functional analysis of genomic datasets. Bioinformatics 23, 259–261. doi: 10.1093/bioinformatics/btl496
Reid, K. E., Olsson, N., Schlosser, J., and Peng, F. (2006). An optimized grapevine RNA isolation procedure and statistical determination of reference genes for real-time RT-PCR during berry development. BMC Plant Biol. 6:27. doi: 10.1186/1471-2229-6-27
Rizhsky, L., Davletova, S., Liang, H., and Mittler, R. (2004). The zinc finger protein Zat12 is required for cytosolic ascorbate peroxidase 1 expression during oxidative stress in Arabidopsis. J. Biol. Chem. 279, 11736–11743. doi: 10.1074/jbc.M313350200
Rojas, C. M., Senthil-Kumar, M., Tzin, V., and Mysore, K. S. (2014). Regulation of primary plant metabolism during plant-pathogen interactions and its contribution to plant defense. Front. Plant Sci. 5:17. doi: 10.3389/fpls.2014.00017
Scalschi, L., Sanmartín, M., Camañes, G., Troncho, P., Sánchez-Serrano, J. J., García-Agustín, P., et al. (2015). Silencing of OPR3 in tomato reveals the role of OPDA in callose deposition during the activation of defense responses against Botrytis cinerea. Plant J. 81, 304–315. doi: 10.1111/tpj.12728
Seifi, H., De Vleesschauwer, D., Aziz, A., and Höfte, M. (2014). Modulating plant primary amino acid metabolism as a necrotrophic virulence strategy: the immune-regulatory role of asparagine synthetase in Botrytis cinerea-tomato interaction. Plant Signal. Behav. 9:e27995. doi: 10.4161/psb.27995
Seifi, H. S., Curvers, K., De Vleesschauwer, D., Delaere, I., Aziz, A., and Höfte, M. (2013). Concurrent overactivation of the cytosolic glutamine synthetase and the GABA shunt in the ABA-deficient sitiens mutant of tomato leads to resistance against Botrytis cinerea. New Phytol. 199, 490–504. doi: 10.1111/nph.12283
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Smith, J. E., Mengesha, B., Tang, H., Mengiste, T., and Bluhm, B. H. (2014). Resistance to Botrytis cinerea in Solanum lycopersicoides involves widespread transcriptional reprogramming. BMC Genomics 15:334. doi: 10.1186/1471-2164-15-334
Snoeijers, S. S., Pérez-García, A., Joosten, M. H. A. J., and De Wit, P. J. G. M. (2000). The effect of nitrogen on disease development and gene expression in bacterial and fungal plant pathogens. Eur. J. Plant Pathol. 106, 493–506. doi: 10.1023/A:1008720704105
Thimm, O., Bläsing, O., Gibon, Y., Nagel, A., Meyer, S., Krüger, P., et al. (2004). mapman: a user−driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 37, 914–939. doi: 10.1111/j.1365-313X.2004.02016.x/pdf
Vanacker, H., Sandalio, L., Jiménez, A., Palma, J. M., Corpas, F. J., Meseguer, V., et al. (2006). Roles for redox regulation in leaf senescence of pea plants grown on different sources of nitrogen nutrition. J. Exp. Bot. 57, 1735–1745. doi: 10.1093/jxb/erl012
Van Baarlen, P., Staats, M., and Van Kan, J. A. (2004). Induction of programmed cell death in lily by the fungal pathogen Botrytis elliptica. Mol. Plant Pathol. 5, 559–574. doi: 10.1111/J.1364-3703.2004.00253.X
Van der Does, D., Leon-Reyes, A., Koornneef, A., Van Verk, M. C., Rodenburg, N., Pauwels, L., et al. (2013). Salicylic acid suppresses jasmonic acid signaling downstream of SCFCOI1-JAZ by targeting GCC promoter motifs via transcription factor ORA59. Plant Cell 25, 744–761. doi: 10.1105/tpc.112.108548
Vidal, E. A., Araus, V., Lu, C., Parry, G., Green, P. J., Coruzzi, G. M., et al. (2010). Nitrate-responsive miR393/AFB3 regulatory module controls root system architecture in Arabidopsis thaliana. Proc. Natl. Acad. Sci. U.S.A. 107, 4477–4482. doi: 10.1073/pnas.0909571107
Walters, D. R., and Bingham, I. J. (2007). Influence of nutrition on disease development caused by fungal pathogens: implications for plant disease control. Ann. Appl. Biol. 151, 307–324. doi: 10.1111/j.1744-7348.2007.00176.x
Wang, R., Okamoto, M., Xing, X., and Crawford, N. M. (2003). Microarray analysis of the nitrate response in Arabidopsis roots and shoots reveals over 1,000 rapidly responding genes and new linkages to glucose, trehalose-6-phosphate, iron, and sulfate metabolism. Plant Physiol. 132, 556–567. doi: 10.1104/pp.103.021253
Wang, R., Tischner, R., Gutiérrez, R. A., Hoffman, M., Xing, X., Chen, M., et al. (2004). Genomic analysis of the nitrate response using a nitrate reductase-null mutant of Arabidopsis. Plant Physiol. 136, 2512–2522. doi: 10.1104/pp.104.044610
Wang, Y. H., Garvin, D. F., and Kochian, L. V. (2001). Nitrate-induced genes in tomato roots. Array analysis reveals novel genes that may play a role in nitrogen nutrition. Plant Physiol. 127, 345–359. doi: 10.1104/pp.127.1.345
Ward, J. L., Forcat, S., Beckmann, M., Bennett, M., Miller, S. J., Baker, J. M., et al. (2010). The metabolic transition during disease following infection of Arabidopsis thaliana by Pseudomonas syringae pv. tomato. Plant J. 63, 443–457. doi: 10.1111/j.1365-313X.2010.04254.x
Weirauch, M. T., Yang, A., Albu, M., Cote, A. G., Montenegro-Montero, A., Drewe, P., et al. (2014). Determination and inference of eukaryotic transcription factor sequence specificity. Cell 158, 1431–1443. doi: 10.1016/j.cell.2014.08.009
Windram, O., Madhou, P., McHattie, S., Hill, C., Hickman, R., Cooke, E., et al. (2012). Arabidopsis defense against Botrytis cinerea: chronology and regulation deciphered by high-resolution temporal transcriptomic analysis. Plant Cell 24, 3530–3557. doi: 10.1105/tpc.112.102046
Zhang, G.-B., Yi, H.-Y., and Gong, J.-M. (2014). The Arabidopsis ethylene/jasmonic acid-NRT signaling module coordinates nitrate reallocation and the trade-off between growth and environmental adaptation. Plant Cell 26, 3984–3998. doi: 10.1105/tpc.114.129296
Keywords: nitrate nutrition, tomato, defense mechanisms, jasmonic acid, gene network analysis, ethylene signaling, microarray analysis
Citation: Vega A, Canessa P, Hoppe G, Retamal I, Moyano TC, Canales J, Gutiérrez RA and Rubilar J (2015) Transcriptome analysis reveals regulatory networks underlying differential susceptibility to Botrytis cinerea in response to nitrogen availability in Solanum lycopersicum. Front. Plant Sci. 6:911. doi: 10.3389/fpls.2015.00911
Received: 29 July 2015; Accepted: 12 October 2015;
Published: 04 November 2015.
Edited by:Stéphane Hacquard, Max Planck Institute for Plant Breeding Research, Germany
Reviewed by:Katherine Denby, University of Warwick, UK
Barbara Blanco-Ulate, University of California, Davis, USA
Copyright © 2015 Vega, Canessa, Hoppe, Retamal, Moyano, Canales, Gutiérrez and Rubilar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Andrea Vega, email@example.com
†These authors have contributed equally to this work.