# Nitrogen limitation alters the response of specific genes to biotic stress

Mahsa Farjad, Martine Rigault, Stephanie Pateyron, Marie-Laure

Martin-Magniette, Annemarie Krapp, Christian Meyer, Mathilde Fagard

<figure><img src="image_1.png" /></figure>

## To cite this version:

Mahsa Farjad, Martine Rigault, Stephanie Pateyron, Marie-Laure Martin-Magniette, Annemarie Krapp, et al.. Nitrogen limitation alters the response of specific genes to biotic stress. International Journal of Molecular Sciences, 2018, 19 (11), pp.3364. 10.3390/ijms19113364. hal-02625219

HAL Id: hal-02625219

https://hal.inrae.fr/hal-02625219

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Article

Nitrogen Limitation Alters the Response of Specific
Genes to Biotic Stress

Mahsa Farjad <fnref n="1" />, Martine Rigault <fnref n="1" />, Stéphanie Pateyron <fnref n="2" />,<fnref n="3" />, Marie-Laure Martin-Magniette <fnref n="2" />,<fnref n="3" />,<fnref n="4" />, Anne Krapp <fnref n="1" />, Christian Meyer <fnref n="1" /> and Mathilde Fagard <fnref n="1" />,*

Received: 3 October 2018; Accepted: 23 October 2018; Published: 27 October 2018

<figure><img src="image_4.png" /></figure>

check for
updates

**Abstract:** In their natural environment, plants are generally confronted with multiple co-occurring stresses. However, the interaction between stresses is not well known and transcriptomic data in response to combined stresses remain scarce. This study aims at characterizing the interaction between transcriptomic responses to biotic stress and nitrogen (N) limitation. Plants were grown in low or full N, infected or not with *Erwinia amylovora* (*Ea*) and plant gene expression was analyzed through microarray and qRT-PCR. Most *Ea*-responsive genes had the same profile (induced/repressed) in response to *Ea* in low and full N. In response to stress combination, one third of modulated transcripts responded in a manner that could not be deduced from their response to each individual stress. Many defense-related genes showed a prioritization of their response to biotic stress over their response to N limitation, which was also observed using *Pseudomonas syringae* as a second pathosystem. Our results indicate an interaction between transcriptomic responses to N and biotic stress. A small fraction of transcripts was prioritized between antagonistic responses, reflecting a preservation of the plant defense program under N limitation. Furthermore, this interaction also led to a complex and specific response in terms of metabolism and cellular homeostasis-associated genes.

**Keywords:** multistress; bacterial phytopathogen; nitrogen limitation; Arabidopsis; transcriptome; defense

**1. Introduction**

Sustainable protection of crops is a worldwide necessity. Indeed, plants are often continuously exposed to a broad range of biotic and abiotic stresses in their natural habitats [1]. Because biotic and abiotic stresses significantly reduce plant growth and productivity, considerable research has aimed to determine the responses of plants to single stresses [2]. However, the response of plants to a combination of stresses is not the simple addition of the response to each single stress [3]. In addition, it has also been reported that plant responses to different stresses are coordinated by complex and interconnected signaling pathways modulating numerous metabolic networks [4]. Apart from some recent reports, the effects of combined biotic and abiotic stress have been little studied. Indeed, various abiotic stresses can negatively or positively affect plant-pathogen interactions. For example, the exposure of Arabidopsis plants to drought enhances plant susceptibility to an avirulent isolate of *Pseudomonas syringae*, while it decreases the susceptibility of tomato to the fungus *Botrytis cinerea* [5].

<footnote n="1">¹ Institut Jean-Pierre Bourgin, INRA, AgroParisTech, CNRS, Université Paris-Saclay, 78000 Versailles, France; farjad_mahsa@yahoo.fr (M.F.); martine.rigault@inra.fr (M.R.); anne.krapp@inra.fr (A.K.); christian.meyer@inra.fr (C.M.)</footnote>

<footnote n="2">² Institute of Plant Sciences Paris-Saclay (IPS2), CNRS, INRA, Université Paris-Sud, Université d'Evry, Université Paris-Saclay, Bâtiment 630, Plateau de Moulon, 91192 Gif sur Yvette, France; stephanie.pateyron@inra.fr (S.P.); marie_laure.martin-magniette@agroparistech.fr (M.-L.M.-M.)</footnote>

<footnote n="3">³ Institute of Plant Sciences Paris-Saclay (IPS2), CNRS, INRA, Université Paris-Diderot, Sorbonne Paris-Cité, Bâtiment 630, Plateau de Moulon, 91192 Gif sur Yvette, France</footnote>

<footnote n="4">⁴ UMR MIA-Paris, AgroParisTech, INRA, Université Paris-Saclay, 75005 Paris, France</footnote>

<footnote>* Correspondence: mathilde.fagard@inra.fr; Tel.: +33-(0)1-30-83-3200</footnote>
One of the most important abiotic factors is the availability of nutrients that seriously affects plant disease severity [6]. Contradictory reports about the impact of nutrients on plant susceptibility to pathogens indicate that many factors affecting this process are not well understood.

Nitrogen (N) is an essential macronutrient and a major limiting factor of plant growth and development [7]. In addition to growth and developmental effects, many agronomical reports highlight the fact that N fertilizers can impact the plant's ability to cope with biotic stress [8]. However, contradictory data have been reported about the effect of N availability on disease development. On the one hand, the effect of N availability on this process seems partly dependent to the lifestyle of the pathogen. Generally, high N availability increases the susceptibility of plants to biotrophs, while it generally decreases the susceptibility of plants to necrotrophs, although some exceptions have been reported [9]. However, this process seems to be more complex and the impact of N availability can be dependent on the host plant for a given pathogen. For example, in the case of *B. cinerea*, one of most important fungal plant pathogens, high N fertilization enhances infection in strawberry [10], while it reduces susceptibility to this fungus in tomato [11]. The form of N available can also determine the effect of N supply on disease severity. For instance, $NH_4^+$ fertilization enhances the level of symptoms caused by *P. syringae*, while, conversely, $NO_3^-$ reduces plant susceptibility to *P. syringae* [12]. These studies indicate the complexity of the relationship between N metabolism and plant resistance to pathogens.

Although the mechanisms are not well known, it is generally thought that N supply can affect plant-pathogen interactions through its impact on plant defense, pathogen virulence, and the availability of nutrients for the pathogen [9]. However, the effect of N limitation on the expression of defense remains unclear, as the literature harbors contradictory reports concerning the effect of N limitation on the expression of defense. Indeed, some reports highlight that N limitation can influence constitutive or induced defense [13]. Many reports show that limiting N supply drives a decrease in defense [14]. Indeed, Arabidopsis plants grown in low N supply show a reduced basal activity of three defense-associated enzymes (chitinase, chitosanase, and peroxidase). Furthermore, in plants grown in low N and treated with BION®, a chemical elicitor of plant defense, the level of these enzymes is also reduced, compared to plants grown in high N [15]. Recently, it was demonstrated that under low N, there is a decrease in the levels of polyamines, compounds known to increase plant resistance via the triggering of programmed cell death [16,17]. However, other reports suggest that there could be a trade-off between plant growth and defense [18]. Thus, although it is obvious that N availability has an impact on plant defense, it is currently difficult to get a clear-cut idea of the effect of N availability on plant defense.

*Erwinia amylovora* (*Ea*) is the bacterial phytopathogen responsible for fire blight disease. *Ea* is a major concern for apple and pear orchards worldwide, as there is no genetic source of complete resistance. We showed previously that *Ea* can multiply in Arabidopsis and that N supply affects the susceptibility of Arabidopsis to *Ea* [9]. In this paper, we studied the impact of N limitation on the response of Arabidopsis to *Ea*. We analyzed the transcriptome of plants grown in low and full N and subjected to biotic stress. Altogether, our data suggest a preponderance of the plant's response to biotic stress in terms of genes modulation over the response to N limitation. In order to determine the genericity of the effect of N supply on biotic stress, we analyzed the impact of N supply on the response of selected defense marker genes to a virulent and an avirulent strain of the bacterial phytopathogen *P. syringae*.

## 2. Results

### 2.1. *Effect of N Limitation on Gene Expression in Response to Erwinia amylovora (Ea)*

To determine to what extent the nutritional status of a plant affects its response to a bacterial pathogen, we grew Arabidopsis plants in low (0.5 mM) or full (5 mM) $NO_3^-$ in plugs of peat moss substrate for 5 weeks (Figure 1A), and infected them with the bacterial phytopathogen *Ea*. In these conditions, the rosette diameter was reduced by 38% (Figure 1B) and nitrate content was reduced by
16% (Figure S1), indicating that NO₃⁻ was limiting for plant growth in the 0.5 mM NO₃⁻ condition compared to 5 mM NO₃⁻ (full N). In order to determine if plant susceptibility to bacteria was affected by these growth conditions, plants were inoculated with *Ea* wild-type strain and in planta bacterial cell numbers was analyzed 6 and 24 h post inoculation (hpi). Bacterial cell numbers were lower in leaves of Arabidopsis plants grown in low N compared to those grown in full N (Figure 1C).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Impact of nitrogen (N) on physiological and transcriptional characteristics of Arabidopsis plants. (A) Five-week-old Arabidopsis rosettes grown under low and full N. (B) Maximum rosette diameter quantified by ImageJ. (C) Bacterial cell numbers of *Erwinia amylovora* (*Ea*) CFBP1430 in Arabidopsis rosette leaves at 6 and 24 h post inoculation (hpi). (B,C) Different letters indicate significant differences according to the Mann-Whitney test (*p*-value < 0.05). All experiments were repeated twice or more with similar results: For (B), the results of two independent experiments were pooled (*n* = 20); for (C), a representative experiment is shown (*n* = 3). (D,E) Transcriptomic analysis of Arabidopsis plants grown in contrasted N regimes and inoculated with *Ea*. Plants were sampled 6 hpi; two biological replicates were performed for each condition. (D) Venn diagrams illustrating the overlap of upregulated (red) and downregulated (green) Arabidopsis genes in response to *Ea* between low and full N. (E) Distribution of functional categories according to the FunCatDB. The pie charts represent: (a) 2982 genes induced and (b) 3248 genes repressed by *Ea* in plants grown in low N; (c) 3017 genes induced and (d) 3401 genes repressed by *Ea* in plants grown under full N. Asterisks (*) indicate significant differences according to hypergeometric distribution (*p*-value < 0.05).</figcaption></figure>

In order to test the hypothesis that lower in planta bacterial cell numbers are due to a difference in
the expression of defense in plants grown in low N, we analyzed the transcriptome of plants grown in
low and full N and infected or not with *Ea*. Leaves of five-week-old plants were *Ea*- or mock-inoculated and sampled 6 hpi. Gene expression was analyzed using the CATMA microarray. We first compared gene expression in noninfected plants to a previous microarray dataset from Arabidopsis plants grown in low nitrate [19]. As expected, several genes previously described as repressed in low nitrate conditions, such as *ASN2* [20], were less expressed in plants grown in low N than in plants grown in full N. Then, log2 of the ratio between expression in infected and mock plants was calculated. In order to validate our dataset, we compared our data to previous transcriptomic analysis performed in response to *Ea* [21] and found that previously identified *Ea*-responsive genes showed, for the most part, a similar profile in the present dataset (Figure S2). The Pearson correlation coefficient between the two sets of data was high (R² = 0.95; Figure S3), indicating that response to infection of plants grown in contrasted N supply conditions was very close. Statistical analysis of the data indicated that the expression of 2982 and 3017 genes was significantly induced (Bonferroni, *p*-value < 0.05) following infection with bacteria, in plants grown in low and full N, respectively. Altogether, 2602 genes showed an increase in expression following infection in both N regimes (Figure 1D). We also found that 3248 and 3401 genes were significantly repressed (Bonferroni, *p*-value < 0.05) following *Ea* infection in plants grown in low and full N, respectively. 2799 genes showed a decrease in expression following infection in both N conditions (Figure 1D). Several hundreds of genes were modulated in response to *Ea* only in one N condition (Figure 1D). However, only a small number of these genes, described later, showed a strong difference in expression profile between the two N regimes (Table S1 and Figure S3).

We then analyzed the functional categories of the infection-modulated genes using the Munich Information Center for Protein Sequences database [22]. For *Ea*-induced and *Ea*-repressed genes, the main functional categories represented were found both for plants grown in low and full N (Figure 1E), which is not surprising, given the similarity of the datasets obtained for plants grown at in low and full N (Figure S3). The largest categories represented unclassified and metabolism-related genes (Figure 1E). Several functional categories were significantly overrepresented in our datasets compared to their representation in the whole genome (categories with an asterisk in Figure 1E). In most cases, their overrepresentation was different between *Ea*-induced and *Ea*-repressed genes, but was not affected by the N regime of the plants. Only the “cell fate” category, which mainly contains genes associated with cell growth, showed a difference according to the $NO_3^-$ level, as it was significantly overrepresented only in plants grown in full N for bacteria-induced genes. Among bacteria-induced genes, the “metabolism”, “protein fate”, cellular transport”, “cellular communication”, and “cell rescue and defense” functional categories were overrepresented. Among bacteria-repressed genes, the “metabolism” and “biogenesis of cellular components” categories were overrepresented. These categories are consistent with the response of the plant to biotic stress, which is known to lead to strong metabolic readjustment, defense setup and a shutting down of the photosynthesis apparatus.

Our data indicate that the very large transcriptional reprogramming previously observed in Arabidopsis leaves infected with *Ea* [21] occurs whether plants are grown in low or full N. Despite the strong impact of the N limitation on the growth of the plants, these plants showed a response to biotic stress that was mostly similar to plants grown in full N, with only strong differences in expression for a small subset of genes.

## 2.2. *N Availability Modulates Specific Defense-Related Genes*

Little is known on the impact of plant nutrition on the expression of genes following biotic stress. To analyze the impact of N supply on known defense-associated gene expression, we compared the gene expression ratio (log2) between infected and mock plants grown in low or full N (Table 1). Most Ethylene- (ET) and Salicylic acid (SA)-related genes analyzed were highly induced by *Ea* infection in plants grown in low and full N. In general, the level of modulation by infection of the SA- and ET-associated genes analyzed here was not strongly affected by N availability. However, two SA-responsive genes, *PR2* and *PR5*, were induced only in plants grown upon N limitation (Table 1). Conversely, most Jasmonic acid (JA)-associated genes were either not modulated following
bacterial infection or were repressed by bacteria. However, we found that two genes involved in JA biosynthesis, AOS (allene oxide synthase) and AOC2 (allene oxide cyclase), and the JA-responsive gene JR1 (Jacalin lectin family protein) were repressed only in plants grown in full N. Only the JA-responsive PR4 gene was induced by Ea, and this was restricted to plants grown under low N conditions. Altogether, these data indicate that N availability modulates, at least in part, the known plant defense response to pathogens, with a strong impact on the JA-signaling pathway.

<table><caption>**Table 1.** Effect of N limitation on defense-related genes. Values are log2 signal ratios between infected and water-treated control plants. BS: Genes involved in hormone biosynthesis, R: Hormone-responsive genes, S: Genes involved in hormone signaling. Two levels of significance threshold were considered according to the Bonferroni statistical test (a: *p*-value < 0.05; b: *p*-value < 10⁻⁸).</caption><thead><tr><th>Gene Number</th><th>Name</th><th>Function</th><th>Low N<br>Ea vs. Mock</th><th>Full N<br>Ea vs. Mock</th><th></th><th></th></tr></thead><tbody><tr><td colspan="7" style="text-align: center; font-weight: bold;">SA biosynthesis and responsive genes</td></tr><tr><td>AT1G74710</td><td>ICS1</td><td>BS</td><td>6.14</td><td>b</td><td>5.73</td><td>b</td></tr><tr><td>AT3G52430</td><td>PAD4</td><td>BS</td><td>5.15</td><td>b</td><td>4.51</td><td>b</td></tr><tr><td>AT3G48090</td><td>EDS1</td><td>S</td><td>3.85</td><td>b</td><td>3.82</td><td>b</td></tr><tr><td>AT4G39030</td><td>EDS5</td><td>S</td><td>6.33</td><td>b</td><td>5.67</td><td>b</td></tr><tr><td>AT1G64280</td><td>NPR1</td><td>S</td><td>2.14</td><td>a</td><td>1.66</td><td>a</td></tr><tr><td>AT2G14610</td><td>PR1</td><td>R</td><td>5.10</td><td>b</td><td>3.23</td><td>b</td></tr><tr><td>AT3G57260</td><td>PR2</td><td>R</td><td>1.50</td><td>a</td><td>0.03</td><td></td></tr><tr><td>AT1G75040</td><td>PR5</td><td>R</td><td>2.00</td><td>b</td><td>0.88</td><td></td></tr><tr><td colspan="7" style="text-align: center; font-weight: bold;">ET biosynthesis and responsive genes</td></tr><tr><td>AT1G05010</td><td>ACO</td><td>BS</td><td>2.64</td><td>b</td><td>1.62</td><td>a</td></tr><tr><td>AT3G04580</td><td>EIN4</td><td>R</td><td>2.76</td><td>b</td><td>2.66</td><td>b</td></tr><tr><td>AT4G17500</td><td>ERF-1</td><td>R</td><td>1.93</td><td>a</td><td>1.92</td><td>a</td></tr><tr><td>AT2G40940</td><td>ERS1</td><td>R</td><td>1.42</td><td>a</td><td>1.05</td><td>a</td></tr><tr><td>AT1G66340</td><td>ETR1</td><td>R</td><td>1.19</td><td></td><td>1.14</td><td></td></tr><tr><td colspan="7" style="text-align: center; font-weight: bold;">JA biosynthesis and responsive genes</td></tr><tr><td>AT3G25770</td><td>AOC2</td><td>BS</td><td>-1.17</td><td></td><td>-2.12</td><td>b</td></tr><tr><td>AT5G42650</td><td>AOS</td><td>BS</td><td>-1.13</td><td></td><td>-2.13</td><td>b</td></tr><tr><td>AT2G46370</td><td>JAR1</td><td>BS</td><td>-0.23</td><td></td><td>-0.51</td><td></td></tr><tr><td>AT3G45140</td><td>LOX2</td><td>BS</td><td>-0.35</td><td></td><td>-0.97</td><td></td></tr><tr><td>AT2G06050</td><td>OPR3</td><td>BS</td><td>0.16</td><td></td><td>0.26</td><td></td></tr><tr><td>AT1G19640</td><td>JMT</td><td>R</td><td>-0.60</td><td></td><td>-0.23</td><td></td></tr><tr><td>AT3G16470</td><td>JR1</td><td>R</td><td>0.08</td><td></td><td>-2.01</td><td>a</td></tr><tr><td>AT2G26020</td><td>PDF1.2b</td><td>R</td><td>-1.77</td><td>a</td><td>-2.06</td><td>a</td></tr><tr><td>AT3G12500</td><td>PR3</td><td>R</td><td>0.06</td><td></td><td>-0.06</td><td></td></tr><tr><td>AT3G04720</td><td>PR4</td><td>R</td><td>2.43</td><td>b</td><td>0.35</td><td></td></tr><tr><td>AT5G24770</td><td>VSP2</td><td>R</td><td>0.14</td><td></td><td>-0.02</td><td></td></tr><tr><td>AT1G32640</td><td>ATMYC2</td><td>S</td><td>-1.70</td><td>a</td><td>-1.82</td><td>a</td></tr><tr><td>AT2G39940</td><td>COI1</td><td>S</td><td>0.25</td><td></td><td>0.10</td><td></td></tr></tbody></table>

We also analyzed the impact of N supply on the response to bacterial infection of known defense-related regulator genes. For example, WRKY transcription factors (TF) are well known regulators of the biotic stress response in plants. The majority of WRKY TFs were indeed modulated by Ea infection, but in most cases, there was no impact of N supply on this modulation (Table S2). Only five WRKY TFs showed a differential modulation in response to bacteria under different N regimes (Table S2): Four were induced only in full N (WRKY42, WRKY47, WRKY64, WRKY67) and two were repressed only in full N (WRKY3 and WRKY69). Concerning MYB TFs, which have been in some cases associated with the defense response [23], most genes were not modulated by Ea or repressed independent of the N regime. Only three genes showed an N-regime specific profile (AT1G25550, AT5G17300 and AT1G74840).
We hypothesized that the genes showing the strongest difference in the amplitude of the modulation by bacterial infection could be good candidates to explain the difference in susceptibility to *Ea* of Arabidopsis plants grown under different N regimes. We thus looked for genes with the highest differential in log2 ratio in response to infection between plants grown in low and full N. The thirty genes with the highest differential in log2 ratio showed mostly two types of profiles (Table S1): Genes repressed by bacterial infection in low N but not in full N, and genes induced by bacterial infection specifically or more highly in low N. Interestingly, genes which were repressed by bacteria only in low N were mainly associated to metabolism (Table S1), while genes which were more or only induced in low N were mainly linked to plant defense responses. For instance, two kelch repeat-containing F-box family proteins (AT1G80440, AT2G44130) were involved in the regulation of the phenylpropanoid biosynthesis pathway [24], an ankyrin repeat protein (AT5G54610) was involved in regulation of immunity [25], and a leucine-rich repeat family protein (AT3G11010) was involved in defense signaling [26].

## 2.3. N Limitation and Biotic Stresses Interact

In order to test the hypothesis of an interaction between the response of plants to N limitation and the response to biotic stress, we analyzed the response of plants to the combination of stresses in more detail. For this, we compared each single stress (N limitation or bacteria) to the combination of the two stresses (Figure 2A). We used the categories defined by the authors of [1] to classify the genes according to their pattern of expression, using a log2 (ratio) of 1 as a cutoff. Furthermore, the responses were considered in a simplified fashion as “induced/not modulated/repressed”, without considering the level of expression. Comparison of the responses of single versus combined stresses showed that around two thirds of the genes modulated in our experiments show an independent response (Figure 2B). Indeed, these genes show a response to the stress combination (N + B) that corresponds to the response of the gene to one of the single stresses. This was expected, since the overlap between infection-modulated genes between plants grown in low and full N is very high (Figure 1D). However, one third of the genes showed a response to the combination of the stresses that could not be deduced from their response to each individual stress, suggesting an interaction between the two single stress responses. These nondeducible gene patterns fall into three categories, defined previously [1], which are combinatorial, cancelled, and prioritized (Figure 2C). Only a few transcripts responded in the prioritized manner (Figure 2C). The prioritized category corresponds exclusively to genes that are repressed by N limitation, induced by biotic stress in full N conditions and that remain induced by biotic stress in low N conditions. The majority of the genes showing a nondeducible pattern showed either a combinatorial pattern (16.4% of total modulated genes) or a cancelled pattern (13.4%, of total modulated genes). Furthermore, the expression profiles corresponded to a small number of specific expression patterns (Figure 2C). In the cancelled category, we found five subcategories (C1–C5; Figure 2C) corresponding to genes induced or repressed by one of the single stresses and not modulated in response to the combination. Among these five categories, the most abundant correspond to genes induced by biotic stress and not modulated in response to the combination (C3) and genes repressed by N limitation and not modulated in response to the combination (C4). The most surprising category corresponds to the combinatorial pattern (Co1–Co2; Figure 2C). These genes correspond to genes that are not modulated in response to single stresses, but that show repression (the most abundant subcategory; Co1) or induction (Co2) in response to the combination.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Gene expression patterns in response to single and combined stresses. (A) Schematic representation of single stresses (N:N limitation; B: *Ea* infection) and combined stresses (N + B:N limitation and *Ea* infection); arrows represent the dataset comparisons performed that are represented as "expression pattern" in (C). (B) Number of genes in the different categories of responses to stress combination; for each category, the percentage among the total modulated genes is indicated. The dotted line indicates the three nonpredictable categories, which represent one third of *Ea*-responsive transcripts. (C) Detail of the expression patterns among the different nonpredictable categories. Expression patterns: The dotted line represents transcript level in control plants (mock/5 mM NO₃⁻); the full lines correspond to N, B, and N + B (from left to right). For each subcategory, only the type of response is indicated with a full line (induced, repressed, no response). Enrichment in Gene ontology (GO) terms in each subcategory was determined according to the FunCatDB. Asterisks (*) indicate a significant difference relative to the Arabidopsis genome (*p*-value < 0.05). P: Prioritized, C: Cancelled, and CO: Combinatorial. Met: Metabolism, Tr: Transcription, Pf: Protein fate, Ct: Cellular transport, Cc: Cellular communication, Cr: Cell rescue, defense, Cf: Cell fate, Bc: Biogenesis of cellular components, Uc: Unclassified.</figcaption></figure>

We then determined whether transcripts of each specific response mode could be linked to particular biological functions (Figure 2C). The prioritized mode was primarily associated to "cell rescue, defense" and "cellular communication". In the canceled and combinatorial categories, the largest subcategories (C3 and Co1) corresponded to a more diverse set of functional categories. In the Co1 pattern, the "metabolism" was very important, while in the C3 category, the "protein fate and cellular transport" functional categories were very important (Figure 2C). This may reflect a complex and specific adaptation of the plant to the combination of stresses.

## 2.4. Defense-Associated Genes Are Prioritized in Response to the Combination of N Limitation with Different Pathogens

Since our analysis of the transcriptome data indicated that defense-associated genes were overrepresented in the prioritized mode, we analyzed the expression profile of specific defense-associated genes known to be expressed in response to biotic stress. We found a large number with a prioritized pattern. This is particularly true for the WRKY family of transcription factors [27] that have largely been described as being involved in the response of plants to stress (Table 2). Another
typical defense gene, PR1, was also found to be regulated in the prioritized mode. We compared these results to N metabolism-associated genes (Table 2). Most N-related genes showed profiles that were independent in the stress combinations (Table S3), including three members of the NRT2 family of putative nitrate transporters, NRT2.1 and NRT.6, known to be involved in the response to biotic stress [28,29]. Only five genes related to N metabolism showed a specific profile in response to the combination of the two stresses (Table 2). Three of these genes showed a prioritized pattern, while two showed a cancelled pattern. Interestingly, the three genes showing a prioritized pattern have been linked to defense responses. Indeed, *Lysine Histidine Transporter 1* (LHT1) and *Ammonium Transporter 1* (AMT1) have been shown to be involved in defense against pathogens [30,31] while *Wound-responsive gene 3* (WR3/NAR1), encoding a component of the high-affinity nitrate transporter system, is involved in JA-independent wound signal transduction [32]. Thus, our data show that genes known to play a role in Arabidopsis defense against pathogens were regulated in a prioritized manner, which is consistent with the fact that N limitation affected the response of the plant to biotic stress only for a limited subset of genes (Figure 2).

<table><caption>**Table 2.** Expression profiles of selected genes in response to single and combined stresses. Nondeducible gene profiles of selected defense-associated genes and N metabolism are presented. Values represent log2 signal ratios of the fold-changes upon single stresses (N:N limitation; B: *Ea* infection) or combined stresses (N + B:N limitation and *Ea* infection). The column at right indicates the transcriptional response mode according to the categories described by Rasmussen et al. 2013. N:N limitation, B: Bacteria, N + B: Combined stresses (N limitation and bacteria).</caption><thead><tr><th>Gene Name</th><th>Gene Number</th><th>N</th><th>B</th><th>N + B</th><th>Response Mode</th></tr></thead><tbody><tr><td colspan="6" style="text-align: center;">Response to salicylic acid</td></tr><tr><td>PR5</td><td>AT1G75040</td><td>-2.07</td><td>0.88</td><td>-0.06</td><td>cancelled</td></tr><tr><td>EDS5</td><td>AT4G39030</td><td>-1.09</td><td>5.67</td><td>5.24</td><td>prioritized</td></tr><tr><td>PR1</td><td>AT2G14610</td><td>-2.71</td><td>3.23</td><td>2.40</td><td>prioritized</td></tr><tr><td colspan="6" style="text-align: center;">ET/JA response and biosynthesis</td></tr><tr><td>JR1</td><td>AT3G16470</td><td>-0.60</td><td>-2.01</td><td>-0.52</td><td>cancelled</td></tr><tr><td>ETR1</td><td>AT1G66340</td><td>-0.41</td><td>1.14</td><td>0.78</td><td>cancelled</td></tr><tr><td>LOX3</td><td>AT1G17420</td><td>-1.15</td><td>2.44</td><td>1.96</td><td>prioritized</td></tr><tr><td colspan="6" style="text-align: center;">EDS1 pathway</td></tr><tr><td>FMO1</td><td>AT1G19250</td><td>-1.28</td><td>5.85</td><td>5.02</td><td>prioritized</td></tr><tr><td>AtNUDT7</td><td>AT4G12720</td><td>-1.10</td><td>3.28</td><td>2.97</td><td>prioritized</td></tr><tr><td colspan="6" style="text-align: center;">WRKY transcription factors</td></tr><tr><td>WRKY75</td><td>AT5G13080</td><td>-1.12</td><td>4.20</td><td>3.72</td><td>prioritized</td></tr><tr><td>WRKY51</td><td>AT5G64810</td><td>-1.50</td><td>5.93</td><td>5.47</td><td>prioritized</td></tr><tr><td>WRKY38</td><td>AT5G22570</td><td>-1.93</td><td>3.68</td><td>3.12</td><td>prioritized</td></tr><tr><td>WRKY25</td><td>AT2G30250</td><td>-1.23</td><td>2.53</td><td>1.96</td><td>prioritized</td></tr><tr><td>WRKY54</td><td>AT2G40750</td><td>-1.04</td><td>2.22</td><td>1.88</td><td>prioritized</td></tr><tr><td>WRKY70</td><td>AT3G56400</td><td>-1.30</td><td>1.69</td><td>1.46</td><td>prioritized</td></tr><tr><td>WRKY53</td><td>AT4G23810</td><td>-1.42</td><td>2.24</td><td>1.22</td><td>prioritized</td></tr><tr><td>WRKY60</td><td>AT2G25000</td><td>-1.15</td><td>1.77</td><td>1.34</td><td>prioritized</td></tr><tr><td>WRKY50</td><td>AT5G26170</td><td>-1.02</td><td>3.45</td><td>2.64</td><td>prioritized</td></tr><tr><td>WRKY30</td><td>AT5G24110</td><td>-1.15</td><td>6.11</td><td>4.90</td><td>prioritized</td></tr><tr><td>WRKY33</td><td>AT2G38470</td><td>-1.43</td><td>4.63</td><td>3.79</td><td>prioritized</td></tr><tr><td colspan="6" style="text-align: center;">Resistance genes</td></tr><tr><td>"LRR family protein"</td><td>AT5G45510</td><td>-1.24</td><td>1.87</td><td>1.43</td><td>prioritized</td></tr><tr><td colspan="6" style="text-align: center;">N metabolism</td></tr><tr><td>GLN1.3</td><td>AT3G17820</td><td>-0.38</td><td>-1.02</td><td>-0.90</td><td>cancelled</td></tr><tr><td>GDH3</td><td>AT3G03910</td><td>-0.06</td><td>1.21</td><td>0.24</td><td>cancelled</td></tr><tr><td>AMT1.1</td><td>AT4G13510</td><td>-1.13</td><td>2.40</td><td>1.98</td><td>prioritized</td></tr><tr><td>WR3</td><td>AT5G50200</td><td>-1.28</td><td>3.38</td><td>2.63</td><td>prioritized</td></tr><tr><td>LHT1</td><td>AT5G40780</td><td>-1.30</td><td>3.04</td><td>2.72</td><td>prioritized</td></tr></tbody></table>
In order to determine how generic the interaction between stresses is, we analyzed the effect of N limitation on the response of Arabidopsis to another bacterial phytopathogen, *P. syringae* pv. *tomato*. We used the virulent strain DC3000 and the avirulent strain DC3000 *avrrpm1*. Arabidopsis plants were grown for five weeks in peat moss in low N or full N, as described above. After 5 weeks, rosette leaves were inoculated with the virulent or the avirulent strain of *P. syringae*. As a control, we also inoculated plants with *Ea*. In order to study the kinetic of gene expression, plants were harvested at 6 and 24 hpi and gene expression was analyzed by qRT-PCR as above. We selected three defense-related genes, *PR1*, *WRKY33* and *WRKY60*, which exhibited a prioritized response and one, *PR5*, which exhibited a cancelled response. Interestingly, we not only found that *PR1* was more highly expressed in response to *Ea* in plants grown under low N, but also that *PR1* was induced by *Ea* as early as 6 hpi only in plants grown in full N (Figure 3A). Furthermore, qRT-PCR analysis of these genes confirmed the prioritized response of these genes in combination of N limitation and *P. syringae*, both for the virulent and the avirulent strain tested (Figure 3B).

<figure><img src="image_7.png" /></figure>

<table><thead><tr><th>Gene</th><th>Treatment</th><th>N</th><th>B</th><th>N+B</th><th>Response type</th></tr></thead><tbody><tr><td rowspan="3">PR1</td><td>Ea</td><td>-2.52a</td><td>3.09b</td><td>5.96c</td><td>prioritized</td></tr><tr><td>Ps vir</td><td>-2.43a</td><td>0.81b</td><td>1.89c</td><td>prioritized</td></tr><tr><td>Ps avir</td><td>-2.43a</td><td>0.68b</td><td>3.64c</td><td>prioritized</td></tr><tr><td rowspan="3">PR5</td><td>Ea</td><td>-1.51a</td><td>0.17b</td><td>0.95c</td><td>cancelled</td></tr><tr><td>Ps vir</td><td>-1.02a</td><td>-0.51a</td><td>0.46b</td><td>cancelled</td></tr><tr><td>Ps avir</td><td>-1.02a</td><td>-1.02a</td><td>0.63b</td><td>cancelled</td></tr><tr><td rowspan="3">WRKY33</td><td>Ea</td><td>-2.53a</td><td>2.27b</td><td>4.09c</td><td>prioritized</td></tr><tr><td>Ps vir</td><td>-0.57a</td><td>0.61b</td><td>1.67b</td><td>prioritized</td></tr><tr><td>Ps avir</td><td>-0.57a</td><td>3.06b</td><td>3.52c</td><td>prioritized</td></tr><tr><td rowspan="3">WRKY60</td><td>Ea</td><td>-1.15a</td><td>1.77b</td><td>1.34b</td><td>prioritized</td></tr><tr><td>Ps vir</td><td>-0.51a</td><td>1.48b</td><td>2.33c</td><td>prioritized</td></tr><tr><td>Ps avir</td><td>-0.51a</td><td>1.42b</td><td>1.92b</td><td>prioritized</td></tr></tbody></table>

**Figure 3.** Impact of N supply on defense-related genes expression in response to different bacterial pathogens. (A) Time-course of *PR1* gene expression at 6 and 24 hpi following mock or *Ea* treatment. (B) Response of selected defense-related genes in response to N, B, and N + B at 6 hpi. (A,B) Col-0 plants were grown in full (5 mM) or low (0.5 mM) NO₃⁻, and infiltrated with bacteria (*Ea* wild type, virulent and a virulent strain of *P. syringae*) or mock-inoculated. Expression is normalized to the *PP2a3* (At1g69960) constitutive gene. Values are log2 signal ratios between infected and mock plants. Similar results were obtained in a minimum of three independent experiments, including three biological replicates each; a representative experiment is shown. The bars represent standard deviation. (A,B): Different letters correspond to values that are significantly different according to the Mann-Whitney statistical test (*p*-value < 0.05).
Our data show that defense-related genes respond in the same manner when N limitation is combined to different bacterial pathogens.

### 3. Discussion

Several studies have demonstrated that the availability of nutrients, in particular of nitrogen (N), influences the outcome of plant-pathogen interactions. Nevertheless, the mechanisms underlying this connection are poorly understood, in part because the effect of N availability on this biotic stress is dependent on the plant-pathogen interaction considered [6,9,14,33]. Therefore, it is currently difficult to define general rules for the impact of N availability on the response of plants to biotic stress. Previous reports have suggested that defense activation in plants grown under N limitation is reduced [15]. However, these data concerned a limited number of defenses and this did not allow to determine whether N limitation affected the response of a plant to biotic stress on a large scale. More recently, a transcriptomic analysis showed that the response of tomato to the fungus *B. cinerea* is affected by N supply [11]. However, data concerning combinations of biotic and abiotic stresses at the transcriptomic level remain scarce [27] and no study has analyzed the impact of N supply on the response to biotic stress in Arabidopsis.

Bacterial infection leads to large modifications in the transcriptomic profile of Arabidopsis plants at early time-points post inoculation [30,31]. This was shown for virulent and avirulent *P. syringae* strains, the response to which mainly differs in timing and intensity [32], as well as for necrotrophic bacterial pathogens, such as *Ea* [6,9,14,33]. To test the hypothesis that the response of plants to *Ea* infection is affected by N limitation, we analyzed the transcriptome of rosette leaves of plants grown in low or full N and exposed to the phytopathogenic bacterium *Ea*. Our data indicated, as expected, a large transcriptomic reprogramming following plant infection by *Ea*, with an important induction of defense-associated genes in both low and full N growth conditions. Close examination of the transcriptomic data indicated a large overlap in the plants' response to bacteria in plants grown in low or full N. Indeed, in our dataset, there was an 86% overlap between genes modulated in low and full N, suggesting that, at least at the qualitative level, the response of plants to biotic stress was close in low and full N. This was rather surprising given the important impact of N limitation on plant size in our experiments. Indeed, the commonly held idea of a tradeoff between plant growth and plant defense has led to the widespread idea that plants have a lower capacity to defend themselves when growth is optimal, thus when nitrogen is abundant [34].

We further analyzed our transcriptome data in response to single (N limitation or bacteria) and combined (N limitation and bacteria) stresses, according to a previous analysis of transcriptomic data performed on several combinations of stresses [15]. As in this previous study, we found that three categories (independent, combinatorial, and cancelled), among the five types of expression profiles identified, represent the most abundant transcript response modes (>95% of the total transcripts). Furthermore, as in this previous study, the deducible profiles (independent and similar) were much more frequent than the nondeducible profiles (combinatorial, canceled, and prioritized), which constituted a third of the total transcripts. In addition, we found that the response of plants to the combination of stresses was closer to its response to bacterial treatment than to N limitation alone. This reflects a dominance of the response to biotic stress over the response to the abiotic stress in our experiments. On the other hand, a recent transcriptomic study on Arabidopsis response to sequential double stresses indicated that plants first subjected to drought or herbivory stress and then infected by *B. cinerea* responded similarly to *B. cinerea* treatment alone. The authors proposed that when two stresses are applied in sequence, plants display a transcriptome profile, which is very similar to the second stress, regardless of the nature of the first stress [11]. Since our experimental setup involves applying N limitation before bacterial infection, one could imagine that the dominance of biotic stress is, at least partially, a result of the experimental setup.

The prioritized category corresponded exclusively to genes repressed by N limitation and induced by bacteria and combined stress. The "cell rescue, defense" was overrepresented among these genes.
This indicates that signaling pathways regulating defense against bacteria are negatively regulated by abiotic stress. Previous studies showed that defense genes activated by *B. cinerea* were repressed under drought stress [35], suggesting that repression of defense expression occurs in response to different abiotic stresses and prioritized in response to different pathogens. We confirmed the prioritized mode for some selected defense genes in response to N limitation and another bacterial phytopathogen, *P. syringae*. The expression of N-related genes has been described to be affected by bacterial infection and recent reports show that some of these genes are indeed involved in plant defense responses to pathogens [1, 8, 10, 36]. Our transcriptomic data indicated that several genes related to N metabolism, such as *LHT1* and *NRT2.6*, were strongly induced by bacteria. Interestingly, we found that these genes showed a prioritized profile, like many defense genes, and not an independent mode, like most N-associated genes. This suggests that although these genes are involved in N metabolism and/or transport, they are regulated by the plant as defense genes independently of N supply. On the other hand, *NRT2.1* was upregulated by bacteria only in plants grown in full N. Interestingly, a recent transcriptomic study on tomato indicated that *NRT2.1* is also upregulated by *B. cinerea* only in plants grown in full N [11]. This suggests that induction of *NRT2.1* is N-dependent in response to different pathogens, suggesting a conserved mechanism. *NRT2.1* has been shown to be a down-regulator of salicylic acid-dependent defenses in response to *P. syringae*. These results suggest that *NRT2.1* could play a role in negatively controlling defense activation in response to pathogens in full N in different pathosystems.

Infection with *Ea* activates defense-related hormonal signaling pathways. The SA-signaling pathway is indeed strongly induced following *Ea* infection in a T3SS-dependent manner, both in host and nonhost plants [37,38]. Regulation of the JA-signaling pathway in response to *Ea* infection seems more complex. Indeed, several genes involved in JA biosynthesis are repressed following *Ea* infection in Arabidopsis, but other JA-dependent genes are induced following infection by *Ea* ([21] and this study). In parallel, it was shown that T3SS-dependent downregulation of the JA pathway is a critical element in the infection process of Malus spp. by *Ea*, since the addition of methyl-jasmonate to susceptible plants increases their resistance to *Ea*. In contrast, the SA pathway was similarly induced in both resistant and susceptible Malus spp by *Ea* [39]. Study of hormonal pathway related genes in our transcriptomic data indicated that SA- and ET-associated genes were highly induced by *Ea*, while most of the JA-related genes were repressed by *Ea*. Generally, in response to *Ea*, the ET and SA pathways were not significantly affected by N supply; however, some genes involved in JA biosynthesis and JA-responsive genes were repressed only in plants grown in full N, suggesting that N supply modulates plant-pathogen interaction through the JA pathway signaling. Thus, the higher bacterial cell numbers in plants grown in full N could be linked to higher repression of JA-associated defense in these plants. Interestingly, these results are reminiscent of those observed with the *B. cinerea*-tomato interaction, for which expression of the JA pathway was clearly identified as being associated with the lower symptoms observed in plants grown in high N [40].

# 4. Materials and Methods

## 4.1. Growth Conditions of Arabidopsis Plants

Seeds of *Arabidopsis thaliana* Col-0 were obtained from the INRA Versailles collection. Plants were grown for 5 weeks in 4 cm plugs of peat moss substrate (70% blond peat, 20% perlite, and 10% vermiculite) wrapped in a nonwoven film [41] and were subjected to an 8 h-light (150 µmol·m⁻²·s⁻¹ irradiation) and 16 h-dark cycle at 21 °C (day)/18 °C (night), with 65% relative humidity. Nitrogen limitation was performed as described in Reference [9]. Briefly: Plants were grown for five weeks in soil and watered to a final humidity of 60%, with a nutrient solution containing full (5 mM NO₃⁻) or low nitrogen (0.5 mM NO₃⁻). In all cases, 5-week-old stressed or control plants were mock or pathogen-inoculated as described in the text.
## 4.2. Pathogen Infections

Rosette leaves of 5-week-old plants were infiltrated with *Ea* CFBP1430 using a needleless syringe. Bacterial suspensions were prepared in sterile water (10⁷ CFU·mL⁻¹). Six and twenty-four hours after infection (hpi), we performed bacterial counting by grinding infected leaves using glass beads in a TissueLyser (Qiagen/Retsch, Hilden, Germany). The bacterial suspensions were used to prepare serial dilutions, which were plated on an LB medium, and after 1 or 2 days the colonies formed were counted to evaluate the initial number of bacteria.

## 4.3. RNA Isolation and qRT-PCR Analysis

For RNA extraction, twelve leaves of three plants (pathogen- or mock-treated) were collected at the indicated time-point after treatment, pooled, and immediately frozen in liquid nitrogen. The experiment was repeated twice independently, thus *n* = 6. Total RNA was extracted from 100 mg of frozen ground leaves using Trizol® reagent (Invitrogen Life Technologies, Saint-Aubin, France). RNA quality was evaluated by electrophoretic run on 1% agarose gel. For the qRT-PCR analysis, first-strand cDNA was synthesized using Superscript reverse transcriptase SSII (Invitrogen, Saint-Aubin, France) from 1 µg of DNase-treated (Invitrogen) total RNA in a 20 µL reaction volume. qPCR reactions were performed using SYBR® Selected MasterMix 2x (Applied Biosystem, Villebon Sur Yvette, France), following the manufacturer's protocol. The cycling conditions consisted of an initial 5 min at 95 °C, followed by 40 three-step cycles at 94 °C for 15 s, 60 °C for 30 s, and 72 °C for 30 s. Melting curve analysis was performed after cycle completion to validate amplicon identity. Relative expression levels were calculated following the standard curve based method [42]. Expression of the *Protein Phosphatase 2a Subunit A3 (PP2A3)* reference gene [43] was used for normalization of every target gene studied. For each treatment, three biological replicates, corresponding to a pool of 4 leaves from a single plant, were analyzed and each qRT-PCR reaction was carried out in duplicate; the complete experiment was conducted twice independently and one representative experiment is presented in Figure 3. The gene-specific primers used in this analysis are indicated in Table S4.

## 4.4. Transcriptome Studies

Microarray analysis was carried out at the Institute of Plant Sciences Paris-Saclay (IPS2, Evry, France), using the CATMAv7 array [34] based on AGILENT technology. The CATMAv7 design of *Arabidopsis thaliana* genome was made with gene annotations included in FLAGdb++ (http://tools.ips2.u-psud.fr/FLAGdb), an integrative database around plant genome [1]. The single high density CATMAv7 microarray slide contains four chambers, each containing 149916 primers. Each 60 bp primer is triplicate in each chamber for robust analysis and in both strands. As part of all probes, 35,754 in triplicate correspond to gene TAIRv8 (among which, 476 probes correspond to mitochondrial and chloroplast genes), together with 1289 probes corresponding to EUGENE software predictions and 658 probes for miRNA/MIR, and finally 240 controls. Two independent biological replicates were produced. For each biological replicate, RNA samples were obtained by pooling RNAs from more than three leaves. Leaves were collected on plants at 3.90 developmental growth stages [35] cultivated in short day conditions. Total RNA was extracted using Trizol®, followed by a purification step on RNeasy column (Qiagen, Courtaboeuf, France) according to the supplier's instructions. For each comparison, one technical replicate with fluorochrome reversal was performed for each biological replicate (i.e., four hybridizations per comparison). The labeling of cRNAs with Cy3-dUTP or Cy5-dUTP was performed as described in Two-Color Microarray-Based Gene Expression Analysis Low Input Quick Amp Labeling manual (©Agilent Technologies, Inc., Les Ulis, France). The hybridization and washing were performed according to Agilent Microarray Hybridization Chamber User Guide instructions (©Agilent Technologies, Inc.). Two-micron scanning was performed with InnoScan900 scanner (InnopsysR, Carbonne, France) and raw data were extracted using MapixR software (version 7.1.0, InnopsysR, Carbonne, France).
Microarray data from this article were deposited in the international repository GEO, Gene Expression Omnibus (Edgar R. 2002, http://www.ncbi.nlm.nih.gov/geo/), accession No. GSE97582) and all steps of the experiment, from growth conditions to bioinformatic and statistical analyses, were detailed in CATdb [28,29] (http://tools.ips2.u-psud.fr/CATdb/; Project: RA14-05_Multipass) according to the "Minimum Information About a Microarray Experiment" standards.

## 4.5. Statistical Analysis of Microarray Data

Experiments were designed with the Genomic networks team of IPS2. For each array, the raw data comprised the logarithm of median feature pixel intensity at wavelengths 635 nm (red) and 532 nm (green). For each array, a global intensity-dependent normalization using the loess procedure [37,38] was performed to correct the dye bias. The differential analysis is based on the log-ratios averaging over the duplicate probes and over the technical replicates. Hence, the number of available data for each gene equals the number of biological replicates, and these were used to calculate the moderated *t*-test [39].

Under the null hypothesis, no evidence that the specific variances vary between probes was highlighted by Limma and, consequently, the moderated *t*-statistic was assumed to follow a standard normal distribution. To control the false discovery rate, adjusted *p*-values found using the optimized FDR approach [44] were calculated. We considered as being differentially expressed the probes with an adjusted *p*-value ≤ 0.05. Analysis was done with the R software (R Development Core Team, 2005, https://cran.r-project.org). The function SqueezeVar of the Limma library was used to smooth the specific variances by computing empirical Bayes posterior means. The library kerfdr was used to calculate the adjusted *p*-values.

## 4.6. Transcriptional Response Modes

Transcript sets were created by grouping genes exhibiting similar expression patterns under single (N limitation or bacteria) and combined stress (N limitation and bacteria) treatments. The responses were considered in a simplified fashion as "induced/not modulated/repressed", without taking into account the level of expression: Log-fold changes between -1 and 1 were considered as nonregulated genes, and log-fold changes higher than 1 and lower than -1 were considered as induced and repressed genes, respectively (Table 2). Plant response to N limitation corresponds to ratios (log2) between control plants grown in low (0.5 mM) and full (5 mM) NO₃⁻; plant response to bacteria corresponds to ratios (log2) between infected and control plants grown in full (5 mM) NO₃⁻; plant response to combined stresses corresponds to ratios (log2) between infected plant grown in low (0.5 mM) NO₃⁻ and control plants grown in full (5 mM) NO₃⁻.

Genes were grouped into 20 subcategories, each representing a specific expression pattern. The 20 subcategories were assembled into five larger categories (cancelled, combinatorial, prioritized, independent, and similar) according to [1]. GO terms associated with each specific transcriptional response profile using the FunCatDB [22].

## 4.7. Data Availability

The datasets generated and analyzed during the current study are available in the CATdb repository: http://tools.ips2.u-psud.fr/cgi-bin/projects/CATdb/consult_project.pl?project_id=402.

# 5. Conclusions

Our work shows that while a large proportion of the Arabidopsis genes showed the same modulation profile (induced or repressed) in response to bacterial infection in two contrasted N regimes, a small number of transcripts showed specific responses to the combination of stresses, including known defense-related genes, thus maybe reflecting a preservation of the plant defense program under N limitation treatment. We found that N limitation had a similar effect on defense gene expression in response to two different bacterial pathogens, indicating at least partial conservation of
the interaction between the response to biotic stress and to N limitation. Finally, our work suggests a main role for the JA-signaling pathway in the impact of N supply on the response of plants to biotic stress, as suggested by a previous report [11]. It is probably of great interest to further study this link in future work.

**Supplementary Materials:** Supplementary materials can be found at http://www.mdpi.com/1422-0067/19/11/3364/s1. Table S1. Genes with the most contrasted response to *E. amylovora* in low and high N. Table S2. Modulation of WRKY TFs by *E. amylovora* in low and high N. Table S3. Response to single and combined stresses of selected N-related genes showing an independent profile. Table S4. Sequence of the gene-specific primers used in this analysis. Figure S1. Nitrate content in Arabidopsis rosette leaves grown under low (0.5 mM) or high (5 mM) nitrate [45]. Figure S2. Modulation profile of previously identified Arabidopsis *E. amylovora*-responsive genes [21]. Figure S3. Correlation of transcriptome data.

**Author Contributions:** Conceptualization, A.K., C.M., and M.F. (Mathilde Fagard); Formal analysis, M.-L.M.-M.; Funding acquisition, A.K., C.M., and M.F. (Mathilde Fagard); Investigation, M.F. (Mahsa Farjad), M.R., S.P., and M.F. (Mathilde Fagard); Supervision, M.F. (Mathilde Fagard); Writing—original draft, M.F. (Mahsa Farjad) and M.F. (Mathilde Fagard); Writing—review & editing, M.F. (Mahsa Farjad), M.-L.M.-M., A.K., C.M., and M.F. (Mathilde Fagard).

**Funding:** This research was funded by INRA BAP, grants Multipass & Nitropath BAP2014_63. The IJPB benefits from the support of the LabEx Saclay Plant Sciences-SPS (ANR-10-LABX-0040-SPS).

**Acknowledgments:** We thank O. Loudet for help in experiment setup. We thank H. Ferry and E. Gilbault for technical assistance. We thank R. Le Hir, D. Expert, M.-C. Soulié, and L. Soubigou-Taconnat for useful comments and critical reading of the manuscript.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

**Abbreviations**

B Bacteria

ET Ethylene

GO Gene Ontology

HPI Hours post inoculation

JA Jasmonic acid

MYB MYeloBlastosis

N Nitrogen

SA Salicylic acid

TF Transcription factor

**References**

1. Rasmussen, S.; Barah, P.; Suarez-Rodriguez, M.C.; Bressendorff, S.; Friis, P.; Costantino, P.; Bones, A.M.; Nielsen, H.B.; Mundy, J. Transcriptome Responses to Combinations of Stresses in Arabidopsis. *Plant Physiol.* **2013**, *161*, 1783–1794. [CrossRef] [PubMed]

2. Chew, Y.H.; Halliday, K.J. A stress-free walk from Arabidopsis to crops. *Curr. Opin. Biotechnol.* **2011**, *22*, 281–286. [CrossRef] [PubMed]

3. Suzuki, N.; Rivero, R.M.; Shulaev, V.; Blumwald, E.; Mittler, R. Abiotic and biotic stress combinations. *New Phytol.* **2014**, *203*, 32–43. [CrossRef] [PubMed]

4. Nakashima, K.; Ito, Y.; Yamaguchi-Shinozaki, K. Transcriptional Regulatory Networks in Response to Abiotic Stresses in Arabidopsis and Grasses. *Plant Physiol.* **2009**, *149*, 88–95. [CrossRef] [PubMed]

5. Atkinson, N.J.; Lilley, C.J.; Urwin, P.E. Identification of Genes Involved in the Response of Arabidopsis to Simultaneous Biotic and Abiotic Stresses. *Plant Physiol.* **2013**, *162*, 2028–2041. [CrossRef] [PubMed]

6. Dordas, C. Role of nutrients in controlling plant diseases in sustainable agriculture. A review. *Agron. Sustain. Dev.* **2008**, *28*, 33–46. [CrossRef]

7. Masclaux-Daubresse, C.; Daniel-Vedele, F.; Dechorgnat, J.; Chardon, F.; Gaufichon, L.; Suzuki, A. Nitrogen uptake, assimilation and remobilization in plants: Challenges for sustainable and productive agriculture. *Ann. Bot.* **2010**, *105*, 1141–1157. [CrossRef] [PubMed]
8. Ballini, E.; Nguyen, T.T.; Morel, J.-B. Diversity and genetics of nitrogen-induced susceptibility to the blast fungus in rice and wheat. *Rice* **2013**, *6*, 32. [CrossRef] [PubMed]

9. Fagard, M.; Launay, A.; Clément, G.; Courtial, J.; Dellagi, A.; Farjad, M.; Krapp, A.; Soulié, M.C.; Masclaux-Daubresse, C. Nitrogen metabolism meets phytopathology. *J. Exp. Bot.* **2014**, *65*, 5643–5656. [CrossRef] [PubMed]

10. Daugaard, H.; Sørensen, L.; Løschenkohl, B. Effect of Plant Spacing, Nitrogen Fertilisation, Post-Harvest Defoliation and Finger Harrowing in the Control of *Botrytis cinerea* Pers. in Strawberry. *Eur. J. Hortic. Sci.* **2016**, *68*, 77–82.

11. Vega, A.; Canessa, P.; Hoppe, G.; Retamal, I.; Moyano, T.C.; Canales, J.; Gutiérrez, R.A.; Rubilar, J. Transcriptome analysis reveals regulatory networks underlying differential susceptibility to *Botrytis cinerea* in response to nitrogen availability in *Solanum lycopersicum*. *Front. Plant Sci.* **2015**, *6*, 1964. [CrossRef] [PubMed]

12. Mur, L.A.J.; Simpson, C.; Kumari, A.; Gupta, A.K.; Gupta, K.J. Moving nitrogen to the centre of plant defence against pathogens. *Ann. Bot.* **2017**, *119*, 703–709. [CrossRef] [PubMed]

13. Ward, J.L.; Forcat, S.; Beckmann, M.; Bennett, M.; Miller, S.J.; Baker, J.M.; Hawkins, N.D.; Vermeer, C.P.; Lu, C.; Lin, W.; et al. The metabolic transition during disease following infection of *Arabidopsis thaliana* by *Pseudomonas syringae* pv. *tomato*. *Plant. J.* **2010**, *63*, 443–457. [CrossRef] [PubMed]

14. Snoeijers, S.; Pérez-García, A.; Joosten, M.; De Wit, P. The effect of nitrogen on disease development and gene expression in bacterial and fungal plant pathogens. *Eur. J. Plant. Pathol.* **2000**, *106*, 493. [CrossRef]

15. Dietrich, C.R.; Ploß, K.; Heil, K. Constitutive and induced resistance to pathogens in *Arabidopsis thaliana* depends on nitrogen supply. *Plant. Cell Environ.* **2004**, *27*, 896–906. [CrossRef]

16. Gupta, K.J.; Brotman, Y.; Segu, S.; Zeier, T.; Zeier, J.; Persijn, S.T.; Cristescu, S.M.; Harren, F.J.; Bauwe, H.; Fernie, A.R.; et al. The form of nitrogen nutrition affects resistance against *Pseudomonas syringae* pv. *phaseolicola* in tobacco. *J. Exp. Bot.* **2012**, *64*, 553–568. [CrossRef] [PubMed]

17. Tiburcio, A.F.; Altabella, T.; Bitrián, M.; Alcázar, R. The roles of polyamines during the lifespan of plants: From development to stress. *Planta* **2014**, *240*, 1–18. [CrossRef] [PubMed]

18. Huot, B.; Yao, J.; Montgomery, B.L.; He, S.Y. Growth-Defense Tradeoffs in Plants: A Balancing Act to Optimize Fitness. *Mol. Plant.* **2014**, *7*, 1267–1287. [CrossRef] [PubMed]

19. Peng, M.; Hudson, D.; Schofield, A.; Tsao, R.; Yang, R.; Gu, H.; Bi, Y.M.; Rothstein, S.J. Adaptation of Arabidopsis to nitrogen limitation involves induction of anthocyanin synthesis which is controlled by the NLA gene. *J. Exp. Bot.* **2008**, *59*, 2933–2944. [CrossRef] [PubMed]

20. Lothier, J.; Gaufichon, L.; Sormani, R.; Lemaître, T.; Azzopardi, M.; Morin, H.; Chardon, F.; Reisdorf-Cren, M.; Avice, J.C.; Masclaux-Daubresse, C. The cytosolic glutamine synthetase GLN1;2 plays a role in the control of plant growth and ammonium homeostasis in Arabidopsis rosettes when nitrate supply is not limiting. *J. Exp. Bot.* **2011**, *62*, 1375–1390. [CrossRef] [PubMed]

21. Moreau, M.; Degrave, A.; Vedel, R.; Bitton, F.; Patrit, O.; Renou, J.P.; Barny, M.A.; Fagard, M. EDS1 contributes to nonhost resistance of *Arabidopsis thaliana* against *Erwinia amylovora*. *Mol. Plant-Microbe Interact.* **2012**, *25*, 421–430. [CrossRef] [PubMed]

22. Ruepp, A. The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. *Nucleic Acids Res.* **2004**, *32*, 5539–5545. [CrossRef] [PubMed]

23. Alves, M.; Dadalto, S.; Gonçalves, A.; de Souza, G.; Barros, V.; Fietto, L. Transcription Factor Functional Protein-Protein Interactions in Plant Defense Responses. *Proteomes* **2014**, *2*, 85–106. [CrossRef] [PubMed]

24. Zhang, X.; Gou, M.; Liu, C.J. Arabidopsis Kelch Repeat F-Box Proteins Regulate Phenylpropanoid Biosynthesis via Controlling the Turnover of Phenylalanine Ammonia-Lyase. *Plant Cell.* **2014**, *25*, 4994–5010. [CrossRef] [PubMed]

25. Yang, Y.; He, M.; Zhu, Z.; Li, S.; Xu, Y.; Zhang, C.; Singer, S.D.; Wang, Y. Identification of the dehydrin gene family from grapevine species and analysis of their responsiveness to various forms of abiotic and biotic stress. *BMC Plant Biol.* **2012**, *12*, 140. [CrossRef] [PubMed]

26. Galon, Y.; Nave, R.; Boyce, J.M.; Nachmias, D.; Knight, M.R.; Fromm, H. Calmodulin-binding transcription activator (CAMTA) 3 mediates biotic defense responses in Arabidopsis. *FEBS Lett.* **2008**, *582*, 943–948. [CrossRef] [PubMed]

27. Birkenbihl, R.P.; Liu, S.; Somssich, I.E. Transcriptional events defining plant immune responses. *Curr. Opin. Plant Biol.* **2017**, *38*, 1–9. [CrossRef] [PubMed]
28. Camanes, G.; Pastor, V.; Cerezo, M.; García-Andrade, J.; Vicedo, B.; García-Agustín, P.; Flors, V. A deletion in NRT2.1 attenuates *Pseudomonas syringae*-induced hormonal perturbation, resulting in primed plant defenses. *Plant Physiol.* **2012**, *158*, 1054–1066. [CrossRef] [PubMed]

29. Dechorgnat, J.; Patrit, O.; Krapp, A.; Fagard, M.; Daniel-Vedele, F. Characterization of the Nrt2.6 Gene in Arabidopsis thaliana: A Link with Plant Response to Biotic and Abiotic Stress. *PLoS ONE* **2012**, *7*, e42491. [CrossRef] [PubMed]

30. Liu, G.; Ji, Y.; Bhuiyan, N.H.; Pilot, G.; Selvaraj, G.; Zou, J.; Wei, Y. Amino Acid Homeostasis Modulates Salicylic Acid-Associated Redox Status and Defense Responses in Arabidopsis. *Plant Cell* **2010**, *22*, 3845–3863. [CrossRef] [PubMed]

31. Pastor, V.; Gamir, J.; Camañes, G.; Cerezo, M.; Sanchez-Bel, P.; Flors, V. Disruption of the ammonium transporter AMT1.1 alters basal defenses generating resistance against *Pseudomonas syringae* and *Plectosphaerella cucumerina*. *Front. Plant Sci.* **2014**, *5*, 838. [CrossRef] [PubMed]

32. Titarenko, E.; Rojo, E.; Leon, J.; Sanchez-Serrano, J.J. Jasmonic acid-dependent and -independent signaling pathways control wound-induced gene activation in *Arabidopsis thaliana*. *Plant Physiol.* **1997**, *115*, 817–826. [CrossRef] [PubMed]

33. Walters, D.R.; Bingham, I.J. Influence of nutrition on disease development caused by fungal pathogens: Implications for plant disease control. *Ann. Appl. Biol.* **2007**, *151*, 307–324. [CrossRef]

34. Massad, T.J.; Dyer, L.A.; Vega, C.G. 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* **2012**, *7*, e47554. [CrossRef] [PubMed]

35. Coolen, S.; Proietti, S.; Hickman, R.; Davila Olivas, N.H.; Huang, P.P.; Van Verk, M.C.; Van Pelt, J.A.; Wittenberg, A.H.; De Vos, M.; Prins, M.; et al. Transcriptome dynamics of Arabidopsis during sequential biotic and abiotic stresses. *Plant. J.* **2016**, *86*, 249–267. [CrossRef] [PubMed]

36. Atkinson, N.J.; Urwin, P.E. The interaction of plant biotic and abiotic stresses: From genes to the field. *J. Exp. Bot.* **2012**, *63*, 3523–3543. [CrossRef] [PubMed]

37. Degrave, A.; Fagard, M.; Perino, C.; Brisset, M.N.; Gaubert, S.; Laroche, S.; Patrit, O.; Barny, M.A. *Erwinia amylovora* type three-secreted proteins trigger cell death and defense responses in *Arabidopsis thaliana*. *Mol. Plant Microb. Interact.* **2008**, *21*, 1076–1086. [CrossRef] [PubMed]

38. Venisse, J.-S.; Malnoy, M.; Faize, M.; Paulin, J.-P.; Brisset, M.-N. Modulation of defense responses of *Malus* spp. during compatible and incompatible interactions with *Erwinia amylovora*. *Mol. Plant-Microbe Interact.* **2002**, *15*, 1204–1212. [CrossRef] [PubMed]

39. De Bernonville, T.D.; Gaucher, M.; Flors, V.; Gaillard, S.; Paulin, J.P.; Dat, J.F.; Brisset, M.N. T3SS-dependent differential modulations of the jasmonic acid pathway in susceptible and resistant genotypes of *Malus* spp. challenged with *Erwinia amylovora*. *Plant Sci.* **2012**, *188–189*, 1–9. [CrossRef] [PubMed]

40. Eulgem, T. Regulation of the Arabidopsis defense transcriptome. *Trends Plant Sci.* **2005**, *10*, 71–78. [CrossRef] [PubMed]

41. Tao, Y.; Xie, Z.; Chen, W.; Glazebrook, J.; Chang, H.S.; Han, B.; Zhu, T.; Zou, G.; Katagiri, F. Quantitative Nature of Arabidopsis Responses during Compatible and Incompatible Interactions with the Bacterial Pathogen *Pseudomonas syringae*. *Plant Cell* **2003**, *15*, 317–330. [CrossRef] [PubMed]

42. Larionov, A.; Krause, A.; Miller, W. A standard curve based method for relative real time PCR data processing. *BMC Bioinform.* **2005**, *6*, 62. [CrossRef] [PubMed]

43. Ceccato, L.; Masiero, S.; Roy, D.S.; Bencivenga, S.; Roig-Villanova, I.; Ditengou, F.A.; Palme, K.; Simon, R.; Colombo, L. Maternal control of PIN1 is required for female gametophyte development in Arabidopsis. *PLoS ONE* **2013**, *8*, e66148. [CrossRef] [PubMed]

44. Storey, J.D.; Tibshirani, R. Statistical significance for genomewide studies. *Proc. Natl. Acad. Sci. USA* **2003**, *100*, 9440–9445. [CrossRef] [PubMed]

45. Miranda, K.M.; Espey, M.G.; Wink, D.A. A Rapid, Simple Spectrophotometric Method for Simultaneous Detection of Nitrate and Nitrite. *Nitric Oxide-Biol. Chem.* **2001**, *5*, 62–71. [CrossRef] [PubMed]

<figure><img src="image_8.png" /></figure>