Skip to main content
  • Research article
  • Open access
  • Published:

Phosphite-induced changes of the transcriptome and secretome in Solanum tuberosum leading to resistance against Phytophthora infestans

Abstract

Background

Potato late blight caused by the oomycete pathogen Phytophthora infestans can lead to immense yield loss. We investigated the transcriptome of Solanum tubersoum (cv. Desiree) and characterized the secretome by quantitative proteomics after foliar application of the protective agent phosphite. We also studied the distribution of phosphite in planta after application and tested transgenic potato lines with impaired in salicylic and jasmonic acid signaling.

Results

Phosphite had a rapid and transient effect on the transcriptome, with a clear response 3 h after treatment. Strikingly this effect lasted less than 24 h, whereas protection was observed throughout all time points tested. In contrast, 67 secretome proteins predominantly associated with cell-wall processes and defense changed in abundance at 48 h after treatment. Transcripts associated with defense, wounding, and oxidative stress constituted the core of the phosphite response. We also observed changes in primary metabolism and cell wall-related processes. These changes were shown not to be due to phosphate depletion or acidification caused by phosphite treatment. Of the phosphite-regulated transcripts 40% also changed with β-aminobutyric acid (BABA) as an elicitor, while the defence gene PR1 was only up-regulated by BABA. Although phosphite was shown to be distributed in planta to parts not directly exposed to phosphite, no protection in leaves without direct foliar application was observed. Furthermore, the analysis of transgenic potato lines indicated that the phosphite-mediated resistance was independent of the plant hormones salicylic and jasmonic acid.

Conclusions

Our study suggests that a rapid phosphite-triggered response is important to confer long-lasting resistance against P. infestans and gives molecular understanding of its successful field applications.

Background

Potato late blight disease caused by the oomycete pathogen Phytophthora infestans is among the most severely damaging diseases of the potato crop. The disease is spread by sporangia and zoospores of the oomycete. Under suitable conditions, the encysted zoospores germinate, penetrate the leaf and form haustoria. The initial biotrophic phase leads into a necrotrophic phase characterized by colonization and sporulation on the leaf tissue which in turn gives rise to zoospores for a new cycle of infection [1]. Although there is naturally occurring resistance in wild potato relatives, sustainable resistance has been difficult to achieve at least partly due to rapid adaptation of the oomycete [2],[3]. The predominant method to control for late blight disease is through frequent use of fungicides. Because of the overall harmful impact of continual fungicide spray and high costs incurred due to spraying there is a need to develop alternative methods to control late blight disease [4].

One of the alternative methods in pest control that has been pursued in a variety of patho-systems is induced resistance [5],[6]. Several reports have shown that plants treated with resistance inducing agents like β-aminobutyric acid (BABA), thiamine (vitamin B1), thiadiazole-7-carbothioc acid S-methyl ester (BTH) and phosphite (H2PO3-) show enhanced resistance after pathogen attack. Recent molecular studies indicate activation of a broad range of defense responses during induced resistance. For example, Ton, et al. [7] showed that BABA-activated defense in Arabidopsis is dependent on a cyclin-dependent kinase-like protein. BABA has also been shown to induce salicylic acid (SA)-dependent induced resistance against P. infestans in potato [8],[9]. Similarly, metabolically inert phosphite-based compounds that have a direct inhibitory effect on the mycelial growth of P. infestans[10]-[13] have also been shown to induce resistance in potato against late blight disease [14]-[16]. In agriculture, phosphite-based compounds are marketed as fertilizers, activators of natural resistance or systemic fungicides and are also widely used in some developing countries where they because of the low risk to human health and environment have been identified as potential alternatives to conventional fungicides [17].

However, little is known about the molecular mechanisms behind phosphite-mediated induced resistance and a better understanding of the underlying molecular mechanisms could assist the development of new plant protection strategies. Studies have shown that phosphite enters the cell via phosphate transporters and interferes with phosphate signaling mechanisms because of its close steric resemblance, which potentially could lead to indirect induction of resistance [18]-[21]. Massoud, et al. [22] demonstrated the importance of the SA pathway in phosphite induced resistance in Arabidopsis against Hyaloperonospora arabidopsidis, while Eshraghi, et al. [23] and Machinandiarena, et al. [16] reported that potassium phosphite generates resistance via excessive accumulation of hydrogen peroxide and PR1 expression in Arabidopsis and potato, respectively. In a study of soluble proteins from potato leaf treated with the phosphite product Confine™, Lim, et al. [24] found increased abundance of proteins involved in SA-dependent defense responses, reactive oxygen species (ROS) and calcium-dependent pathways, whereas proteins involved in amino acid and starch metabolism were down-regulated. They also showed induction of hypersensitive response and callose formation after pathogen attack in phosphite treated leaves. This indicates that complex and multiple processes are involved in phosphite-induced resistance.

No previous studies have combined genome-wide transcriptomics with quantitative proteomics on phosphite treated plants. In this study a transcriptomic analysis was complemented by proteomic investigation of secreted proteins in the apoplast, which is regarded as an important interface during plant-pathogen interactions [25]. In parallel, phosphite and phosphate levels were measured in leaves to determine the effect of phosphite on phosphate uptake and metabolism. The level of phosphite-mediated resistance against P. infestans was also investigated with transgenic lines deficient in salicylic and jasmonic acid signaling.

Results

Phosphite-induced protection against P. infestans is observed in detached leaflet assay throughout the time series

A clear reduction of P. infestans infection was observed at all tested time points after phosphite treatment, i.e. 3, 6, 11, 24 and 120 h (Figure 1). Furthermore, the phosphite treated leaflets showed “hR-like” (“hypersensitive response like”) symptoms at the site of P. infestans inoculation compared to extensive sporulation observed on control leaflets 7 days after infection (Additional file 1: Figure S1A). To test whether phosphite was transported in the potato plant and triggered a systemic induced resistance against P. infestans, some leaves were covered during treatment. These “covered leaves” were equally infected as the water controls (Figure 1). In order to test if there is a possible direct-effect of phosphite on the pathogen, a set of plants were washed with tap water to remove phosphite from the leaf surface and dried for a minimum of 5 h before sampling for detached leaflet assay. These leaflets, hereafter referred to as “washed”, had similar protection as observed in phosphite treated leaflets (Figure 1).

Figure 1
figure 1

Detached leaflet assay of potato plants. Potatoes (cv. Desiree) were foliar sprayed either with 36 mM proalexin (Potassium phosphite; phosphite treated) or tap water (Water treated). “Covered leaves” leaflets were obtained by covering two leaves per plant during phosphite spray. Washed leaflets were obtained by spraying leaves with 36 mM proalexin, washing and drying away the phosphite present on the leaves. Infection was measured as lesion size 7 days after inoculation with P. infestans. Data corresponds to mean ± SD obtained from 12 biological replicates.

Phosphite is distributed rapidly in planta

An enzymatic assay was used to measure phosphite concentrations in potato leaflets. Accumulation of phosphite was observed in both phosphite treated and “covered leaves” sampled already 3 h after treatment (Figure 2A). Also at 6, 24 and 120 h after treatment, phosphite was detected in both types of leaflets with no significant difference. In order to estimate phosphite levels on the leaf surface, “washed” leaflets were compared to phosphite-treated leaflets. However, no significant difference in phosphite levels were seen between “washed” and “phosphite treated” leaflets (p > 0.05), indicating that phosphite is either taken up completely or it is not possible to remove phosphite from the leaf surface by excessive washing.

Figure 2
figure 2

Phosphite and phosphate measurements. Phosphite (A) and phosphate measurement (B) of water sprayed, phosphite treated and “covered leaves” leaflets based on an enzymatic activity assay. Data presented here corresponds to mean ± SD obtained from 3 biological replicates.

No significant difference (p = 0.53) in phosphate levels was detected in phosphite treated leaflets and “covered leaves” leaflets (Figure 2B). In addition, no significant difference (p = 0.19) was observed in phosphate levels of phosphite treated leaflets across time-points. Also phosphate levels did not significantly differ between phosphite treated leaflets and water-sprayed control leaflets (Figure 2B). Similar range of phosphite levels (100-1200 ug/g FW vs. 100-900 ug/g FW by Borza, et al. [26] and phosphate (>1000 ug/g FW) was identified in leaf tissue after Confine™, a phosphite-based fungicide, application on potato when measured by ion chromatography.

Phosphite treatment has a rapid and transient effect on the leaf transcriptome

Analysis of the microarray data revealed that significant differential expression of transcripts was observed only at 3, 6 and 11 h after phosphite treatment (Figure 3A). In the result section and discussion we only highlight transcripts up- or down-regulated more than two times compared to the control. No transcript was observed to be differentially expressed in leaflets sampled at 24 h, while at 48 h after treatment only one transcript annotated as encoding reticuline oxidase (DMP400036592) was differentially expressed. Since altered pH can effect transcript levels, leaves were sprayed with acidified water (pH 5.2), but no transcript changes were seen at the two time points tested by microarrays, 3 and 11 h after treatment (data not shown). Phosphite has a rapid effect on the leaf transcriptome with 738 transcripts differentially expressed in leaves sampled at 3 h after treatment. The number of differentially expressed transcripts increased to 5788 at 6 h after treatment and decreased to 4418 at 11 h after treatment (Figure 3A). A comparison of the number of up-regulated and repressed transcripts at each time point revealed that 87% of significantly differentially expressed transcripts at 3 h are up-regulated, while this number drops to 57% and 49% at 6 and 11 h, respectively (Figure 3A). A larger overlap of differentially expressed transcripts at time points 6 and 11 h was observed and a “core” of 207 transcripts were differentially expressed at 3, 6 and 11 h (Figure 3B). The expression levels of eight genes were validated by qPCR and a high agreement between microarrays and qPCR was found (R = 0.9; Additional file 2: Figure S2).

Figure 3
figure 3

Differentially expressed genes. A comparison of number of transcripts induced and repressed at each time point (A), area proportional Venn diagram depicting overlap of transcripts significantly altered at all the time points (B).

The “core” phosphite-induced transcripts are associated with biotic and abiotic stress responses

GO analysis of the 207 transcripts significantly up-regulated across all time points identified 174 enriched GO terms. Based on semantic grouping of the enriched GO terms using ReviGO [27], two major clusters were identified (Figure 4A). One cluster contained GO terms such as response to: wounding (GO:0009611), defense (GO:0006952), chitin (GO:0010200), chemical stimulus (GO:0042221) and SA stimulus (GO:0009751). The second cluster was associated with metabolic activity and contained GO terms such as purine nucleotide metabolic process (GO:0006163), nucleoside biosynthetic process (GO:0009163). GO terms that were enriched but did not belong to either of the two clusters were associated with, e.g., phosphate-containing compound metabolic process (GO:0006796) and organophosphate biosynthesis (GO:0019637) (Figure 4A). Transcripts encoding pentacyclic triterpene synthase (PEN1; DMP400036965), Jasmonate ZIM-domain protein 1 (DMP400005281), WRKY transcription factor-30 (DMP400010347), non-race specific disease resistance (NDR1/NHL25;DMP400037182), cytochrome BC1 synthesis protein (BCS1; DMP400036856), wall associated kinase (WAK; DMP400053153) associated with biotic stress were up-regulated across all the time points in addition to transcripts associated with abiotic stress such as ethylene-responsive transcriptional co-activator (DMP400051868), ZPT2-13 (DMP400027219) and salt responsive protein (DMP400000646) (Additional file 3: Table S1).

Figure 4
figure 4

Gene ontology (GO) analysis. Clusters of enriched GO terms based on semantic similarity identified among transcripts expressed at all the time points (A), representation of processes and associated example transcripts (in brackets) significantly regulated at each time point (B).

Phosphite has a complex effect on stress pathways in addition to affecting processes associated with primary and secondary metabolism

GO analysis of differentially expressed transcripts at 3 h identified enriched GO terms associated with biotic and abiotic stress (Figure 4B, Additional file 4: Figure S3). An increase in expression of defense related transcripts encoding ascorbate oxidase (DMP400009720), cell wall peroxidase (DMP400026523) and phenyl ammonia lyase (PAL; DMP400037388) were observed at this time point. Abiotic stress responsive transcripts encoding salt responsive protein 2 (DMP400000646), early response to dehydration 7 (ERD 7;DMP400010986) associated with dehydration, and transcription factor ZPT2-13 (DMP400027220) were also induced. Transcripts encoding glucan endo-1,3- β -glucosidase (DMP400036382) involved in cell wall biosynthesis and serine/threonine-protein kinase bri1 (DMP400033050) involved in brassinosteroid hormone perception were instead repressed. In addition to activation of defense processes, analysis of differentially expressed transcripts at 6 h revealed significant changes in lipoprotein, amino acid biosynthesis and polysaccharide metabolism (Figure 4B, Additional file 4: Figure S3). Increased expression of a neutral invertase (DMP400023171) transcript involved in sucrose breakdown and phospholipase PLDb1 (DMP400039419) involved in cell wall phospholipid biosynthesis was observed. Significant up-regulation was also observed in genes encoding alternative oxidase (DMP400013470), citrate synthase (DMP400023850), malic enzyme (DMP400004672) and phosphofructokinase (DMP400030437), all of which have been shown to be involved in stress and primary metabolism. Significant up-regulation was also observed in transcripts belonging to amino acid metabolism pathways such as tryptophan synthase (DMP400019982) and chorismate synthase (DMP400042200). Chloroplast omega-6 fatty acid desaturase (DMP400019659) associated with lipid metabolism, phosphatidylinositol-4-phosphate kinase (PIP; DMP400032031) associated with cell wall modification and cellulose synthase (DMP400007078) associated with cell wall biosynthesis were repressed. Analyses of transcripts differentially expressed at 11 h suggest that in addition to activation of stress responses phosphite also has an effect on cell wall related processes (Figure 4B, Additional file 4: Figure S3). Stress related transcripts encoding UDP-Glucosyltransferase (DMP400021191), Glutathione S transferase (DMP400003866), fatty acid desaturase (DMP400041239) and linalool synthase (DMP400048541) were up-regulated. In addition, ontology terms such as response to chemical stimulus (GO:0042221), response to cyclopentenone (GO:0010583) and response to stress (GO:0006950) were enriched. Interestingly, numerous transcripts associated with cell wall related processes were repressed, e.g., pectinesterase (DMP400009250), polygalactouranase (DMP400023907), glycine rich cell wall protein (DMP400050455), chitin binding lectin (DMP400055565) and β-1-3 glucan synthase (DMP400049943). In coherence with this, GO terms such as cell wall modification (GO: 0042545), phosphatidylinositol 3-kinase activity (GO: 0035004) and external encapsulating structure (GO: 0030312) were significantly enriched for the repressed transcripts at 11 h after phosphite treatment (Figure 4B, Additional file 4: Figure S3).

Phosphite has an effect similar to BABA on the leaf transcriptome

A large overlap of differentially expressed genes was observed between phosphite and 48 h after treatment with 10 mM BABA (Figure 5; [28]). Over 300 GO terms were significantly enriched among transcripts that were common to both BABA and phosphite 3 h after treatment (Additional file 5: Table S2). Numerous stress and defense related transcripts such as jasmonate ZIM-domain protein 1 (DMP400005280), thaumatin (DMP400007586), wall-associated kinase (DMP400010164) were common for both treatments. In addition to the large overlap, a high correlation in expression levels of transcripts affected by both BABA and phosphite was observed (Figure 5).

Figure 5
figure 5

Phosphite-BABA comparison. Scatter plots displaying correlation of expression levels of transcripts 48 h after BABA and 3, 6 and 11 h after phosphite (Proalexin) treatment.

Principal component analysis (PCA) of transcriptomic datasets obtained from microarray experiments with potato clones either susceptible or resistant to P. infestans (Ali et al. submitted) revealed that the phosphite responsive transcriptome clustered closer to BABA than to uninfected late blight resistant potato clones (cv. Sarpo Mira and cv. SW93-1015) or uninfected susceptible clone (cv. Des; Additional file 6: Figure S4). This suggests that phosphite treatment triggers induced resistance responses on the transcriptomic level that is similar to BABA treatment.

Secreted proteins involved in cell wall related processes and defense differ in abundance 48 h after treatment

Since no major changes in gene expression was observed from 24 h and onwards, we decided to explore the secretome at a later time point, 48 h, to see whether the marked changes in gene expression were persistent at the protein level. We detected a total of 67 proteins with significant different abundance, whereas only 4 proteins were repressed (Table 1). Mapman visualization showed that proteins associated with cell wall and stress related processes had increased abundance (Additional file 7: Figure S5). Proteins such as ceramidase (DMP400007784), aspartic proteinase nepenthesin-1 (DMP400009572), β-D-glucan exohydrolase (DMP400010541), alpha-galactosidase (DMP400018078), mannose binding lectin (DMP400012540), fasciclin-like arabinogalactan protein 13 (DMP400022582) were observed to change significantly in abundance (Table 1). Significant changes were also observed in stress related proteins such as Kunitz trypsin inhibitor (DMP400046981), peroxidase (DMP400043335), polygalacturonase inhibiting protein (DMP400014905) and Class III peroxidase (DMP400046178). 20 of the proteins identified in this secreted fraction were also observed to change in abundance in the secretome 48 h after BABA treatment [28]. At least 3 proteins encoding chitinases (DMP400002757, DMP400015232 and DMP400015454) were identified in proteins that change in abundance both after phosphite and BABA treatment. However, PR1 induced by BABA was not seen to be induced by phosphite.

Table 1 Secreted proteins changing in abundance 48 h post phosphite treatment

Similar protection observed in phosphite treated jasmonic and salicylic acid potato transgenes

In order to test if phosphite-mediated induced resistance was dependent on SA as shown previously by Massoud, et al. [22] in Arabidopsis, a whole plant infection assay of phosphite and water treated transgenic lines impaired in SA and JA hormone signaling was performed. Analyses of late blight disease symptoms 3, 5 and 7 days post infection revealed that phosphite conferred protection regardless of transgenic line and very few or no symptoms were observed in phosphite-treated transgenic lines in comparison to water treated controls (Figure 6). In untreated plants the NahG line impaired in SA signaling had a lower resistance to P. infestans as reported previously (Figure 6A; [29]).

Figure 6
figure 6

Whole-plant P. infestans infection. Progress of disease (% infection) 3, 5 and 7 days post infection in water (A) and phosphite treated salicylic acid impaired (NaHG), jasmonic acid impaired knockout (RNAi AOC, RNAi OPR3 and RNAi COIX5) and wild type (cv. Desiree) plants (B). Data presented here corresponds to mean ± SE obtained from 3 biological replicates.

Discussion

Phosphite provides efficient disease protection in many different plant pathosystems. In the present study we observed resistance against P. infestans already at 3 h after phosphite treatment and resistance was sustained even on the 5th day after treatment (Figure 1). The phosphite treatment was accompanied by altered expression of hundreds of transcripts. This study shows for the first time that massive phosphite-induced transcriptomic changes occur rapidly and result in early and sustained protection. Intriguingly, phosphite was detected in parts of treated plants not directly exposed to phosphite already at 3 h after treatment (Figure 2A), demonstrating instantaneous mobility of the molecule inside the plant. Still, phosphite-induced protection was absent in leaves not directly treated with phosphite (Figure 1), suggesting that phosphite treatment leads to local resistance and requires direct foliar application in potato. The in planta mobility of phosphite observed is consistent with a recent study by Borza et al. [26], who reported an active translocation of phosphite from leaves to tubers in field-grown potato.

Previous studies have shown that the presence of phosphite impacts phosphate sensing, uptake and metabolism especially when plants are grown in low phosphate conditions [19],[21],[30]-[33]. However, in our study no clear differences in phosphate levels were seen within 5 days of phosphite treatment and levels of phosphate were constantly higher than phosphite (Figure 2B) suggesting that phosphite does not change the plant’s phosphate status with our fertilization regime and that this is not the likely cue behind the increase in resistance.

Phosphite is also known to be directly toxic to P. infestans[11],[13] depending on the concentration used [12],[22]. In field applications, the dual nature of the phosphite molecule both being an inducer of plant resistance and having a direct toxic effect on oomycetes might explain the high efficacy. We used the recommended concentration of phosphite that has also resulted in a significant reduction in potato late blight disease in field studies (Liljeroth et al., unpublished data). Although observation of “hR-like” symptoms 7 dpi on phosphite treated leaflets (Additional file 1: Figure S1a), there could still be a direct toxic effect of phosphite. With this in mind, we set out to test if the presence of phosphite on the surface of the leaflets had a direct effect on P. infestans infection by washing phosphite-treated leaves before infection. These “washed” leaflets showed the same level of resistance as treated leaflets (Figure 1), indicating that a direct interaction on the leaf surface between P. infestans and phosphite was not necessary for resistance. Furthermore, potassium phosphite causes acidification (in our case pH 5.2) and external acidification to pH 4.5 has previously been shown to change the expression of hundreds of genes in Arabidopsis[34]. Therefore we tested whether foliar spraying with acidified tap water (pH 5.2) caused an acid-mediated stress leading to increased protection. However, no increased resistance was observed by this treatment (Additional file 1: Figure S1b) and microarray analysis did not show a change in gene expression compared to the control at two time points tested, 3 and 11 h after acidified water treatment. This indicates that phosphite-mediated protection cannot be attributed to its presence on the leaf surface or to its acidifying property. We could also measure phosphite in “covered leaves” without any protection against P. infestans, indicating that a strong direct toxic effect is not at play in our experiments.

Machinandiarena, et al. [16] and Massoud, et al. [22] studied a smaller selection of defense related genes first at 72 h after phosphite treatment and did not observe any significant induction of genes then. This is consistent with our observations where we only detected significant changes in the transcriptome at earlier time points with the exception of one transcript at 48 h. In a pre-study, the lack of changes in the transcriptome was confirmed with Aliette, an aluminum-phosphite compound, at 48 h (data not shown). In a proteomic study directed towards soluble proteins, Lim, et al. [24] identified significant changes in abundance of numerous proteins after phosphite treatment in potato. However, after infection of phosphite treated leaves with P. infestans, only 4 defense related proteins significantly changed in abundance in comparison to infected control, this led them to hypothesize that phosphite induces pre-activation of defense prior to pathogen inoculation. Indeed, the transcriptomic data analysis presented here reveals that phosphite has an early and transient effect on the transcriptome, which disappears already at 24 h after treatment. Our secretome data also underpins a strong effect of phosphite treatment.

The rapid effect of phosphite 3 h after treatment on the transcriptome with hundreds of transcripts affected includes activation of genes associated to both biotic and abiotic stress response (Figure 4B, Additional file 4: Figure S3). A transcript encoding PAL, a classic marker in plant-pathogen interaction studies associated with SA synthesis, was observed to be up-regulated more than 10-times. BCS1 and NDR1, also associated with SA mediated immunity, were induced as well [35]-[37]. Of special interest is increased expression of a wall associated kinase (WAK) transcript since WAKs previously have been implicated in mediating defense response during pathogen and wound induced cell-wall damage [38]. Numerous abiotic stress responsive transcripts were also markedly up-regulated. For example, multiple ZPT2-13 transcripts that belong to C2H2 type transcription factor family, which share high sequence homology with Arabidopsis ZAT transcription factors that have been shown to be induced in response to abiotic stress [39].

In addition to activation of stress responses, significant changes in transcripts relating to primary metabolism were observed 6 h after phosphite treatment (Figure 4B, Additional file 4: Figure S3). Reprogramming of primary metabolic processes is one of the hallmarks of plants stress response exemplified by the observed induction of neutral invertases which are known to be rapidly induced during incompatible interactions [40]. Expression of TCA cycle related enzymes such as malic enzyme and citrate synthase also increased, which is an additional indication of heightened defense responses related to metabolism [40]. Fructose bisphosphate aldolase transcript was up-regulated and increased abundance of this protein has been reported in potato leaves after treatment with different inducing agents [41]. On the contrary, Lim, et al. [24] reported decreased abundance of this protein in potato leaves after phosphite treatment by quantitative proteomics. In addition, they observed decreased abundance of sucrose synthase and other proteins related to carbon fixation, carbohydrate metabolism and energy production, whereas we observe an induction of the sucrose synthase transcripts 6 h after phosphite treatment and generally see an up-regulation of transcripts related to these processes. These differences in observations can be due to divergence in transcript and protein regulation, but also due to the differences in experimental set-ups; in the study by Lim, et al. [24] field-grown samples were collected after a series of 5 applications of phosphite 3 days after the final application. The effect of sampling time is also evident in our transcriptomic study as significant repression in processes associated with metabolism were observed 11 h after treatment when compared in samples obtained at 6 h (Additional file 8: Figure S6). In addition we also observed changes on the levels of secreted proteins at 48 h after treatment.

Increased expression of transcripts encoding omega-6 fatty acid desaturases was observed, suggesting heightened defense response after phosphite treatment. These genes have been shown to be repressed during compatible potato-P. infestans interactions [42]. Changes in transcripts related to cell wall related processes were also observed (Figure 4B), for example UDP-forming cellulose synthase involved in cellulose biosynthesis. Knock-outs of Arabidopsis homologs of cellulose synthase (IRX5, IRX3 and IRX1) have previously been shown to confer SA, JA and ethylene-independent resistance to Ralstonia solanaceareum and Plectospharella cucmerina[43]. Previously, Yaeno, et al. [44] demonstrated that avirulence factor Avr3a of Phytophthora capsici requires host PIP kinase to stabilize and cause cell death in planta in tobacco, suggesting that PIP could be a susceptibility factor and its down-regulation might play a role in early phosphite-mediated induced resistance. Indeed, in this study a PIP kinase associated with various cell wall related processes was repressed both at 6 and 11 h after phosphite treatment.

After 11 h of phosphite treatment a much more pronounced effect on several cell wall modification processes was observed (Figure 4B, Additional file 4: Figure S3). By clustering Arabidopsis microarray data obtained under different stress and chemical treatments, Ma, et al. [45] observed over-representation of cell wall modification processes in clusters with down-regulated genes. Among the cell wall related transcripts, repression of β-1,3-glucan synthase associated to callose biosynthesis. Previous reports have revealed that callose formation can negatively regulate SA-mediated response by physically secluding the site of pathogen penetration [46],[47].

Since no major changes in gene expression was observed from 24 h and onwards, we decided to explore the secretome at a later time point to see whether the marked change in gene expression was persistent at the protein level. 63 secreted proteins changed in abundance indicating a prolonged change in protein abundance compared to gene expression. In silico functional analysis of these secreted proteins revealed increased abundance of proteins related to cell wall and defense processes (Additional file 7: Figure S5). For example, an aspartic proteinase nepenthesin-1 (Table 1) was found to increase in abundance after phosphite treatment. Numerous secreted aspartic proteases such as constitutive disease resistance 1 (CDR1) have been previously implicated in enhancing resistance in Arabidopsis [48]. A mannose-binding lectin (Table 1) also significantly increased in abundance. These have been shown to play a crucial role during plant pathogen interaction by aiding in recognition of specific glyco-conjugates present on the surface of bacteria and fungi [49]. In addition, a polygalacturonase inhibiting protein, a cell-wall associated protein with well-documented roles in plant defense, and which if over-expressed has been shown to change cell-wall properties prior to pathogen infection, increased in abundance [50],[51] further indicating that phosphite effects secreted proteins related to defense and cell wall at later time points. A total of 20 proteins that changed in abundance after phosphite treatment were also found to change after BABA treatment. Proteins encoding chitinases that are members of pathogenesis-related proteins and are produced in response to biotic and abiotic stress [52] were identified among these proteins. This probably indicates the reason behind efficient antifungal activity by phosphite based compounds as mentioned by Deliopoulos, et al. [53] and others since fungal cell walls are primarily composed of chitin.

Massoud, et al. [22] previously demonstrated the importance of the SA pathway in what they described as “priming” when infecting phosphite treated Arabidopsis with the oomycete Hyaloperonospora arabidopsidis. Conversely, we observed similar protection on both phosphite-treated transgenic lines impaired in SA and JA signaling, either suggesting that SA and JA act in complementation to mediate induced resistance after phosphite treatment or that the induced resistance response is SA and JA-independent (Figure 6). We observed expression in transcripts associated to both JA- and SA-metabolism and signaling, such as allene oxide synthase and 12-oxophytodienoate reductase 3 both involved in JA biosynthesis and phenylalanine ammonia-lyase and phytoalexin-deficient 4-2 protein involved in SA synthesis and signaling (Additional file 9: Table S3), pointing towards importance of both SA and JA during phosphite-mediated response.

Surprisingly, ~40% overlap among transcripts induced by both phosphite and BABA [28] was observed, a similar trend was also observed from large scale clustering analysis of microarray data, with phosphite effect on the transcriptome clustering closer to BABA than to potato lines resistant to late blight disease (Additional file 6: Figure S4). Interestingly though, we observed significant induction of PR1 after BABA treatment [9], whereas PR1 transcript was not induced after phosphite treatment despite similar protection levels conferred against P. infestans.

Conclusion

This is the first detailed investigation of transcriptomic and proteomic changes after phosphite treated leaves prior to pathogen infection. It seems that multiple defense pathways are rapidly induced by phosphite treatment that causes heightened defense leading to enhanced resistance after pathogen infection in local tissue. Our results also indicate that phosphite treatment influences primary metabolism and cell wall associated metabolic processes, and detailed investigation of these processes, e.g. study of cell wall composition and structure, will deepen the understanding of induced resistance mediated by phosphite. Recently, Arabidopsis plants engineered to use phosphite as their source of phosphate were produced in order to outcompete weeds in a low-phosphate environment [54]. In addition use of phosphite or a combination of phosphite and fungicide has also been shown to be comparatively less expensive fungicides alone in downy mildew disease management in grapevines [55]. A future agricultural system with the relatively inexpensive phosphite as the phosphorus source would have the additional benefit in that phosphite reduces the spread of oomycete and other plant pathogens since it triggers the plant’s innate immune system. It is also evident from research conducted by us and others that there are benefits using phosphite-based salts in various plant protection strategies and improved understanding of induced resistance by phosphite could facilitate expansion of its usage.

Methods

Plant material

Solanum tuberosum cv. Desiree was grown in 3.5 liter pots in climate chambers for four to five weeks with 16 h light and 8 h dark at 22°C. Plants were foliar sprayed with 40 mL tap water (control) or with 40 mL 1.25% (v/v) Proalexin (LMI AB, Helsingborg, Sweden) diluted in tap water corresponding to 36 mM phosphite as per the recommendation of the manufacturer. Sampling was done at 3, 6, 11, 24, 48 and 120 h. To demonstrate phosphite redistribution in the plant and possible systemic effects on resistance, two leaves of phosphite treated plants were covered with transparent plastic bags to during Proalexin treatment, and are referred to as “Covered leaves”. To test whether the acidity of Proalexin could trigger induced resistance, plants were sprayed with acidified tap water adjusted with hydrochloric acid (HCl) to pH 5.4 (equivalent to the pH of the final Proalexin solution). In order to remove Proalexin from the leaf surface, Proalexin-treated plants were washed in tap water and then let to dry for at least 5 h before harvest. Sampling for detached leaflet assay and determination of phosphite levels of so called “washed” leaflets was done for the 6, 11, 24 and 120 h time points. Each treatment was done in three biological replicates at each time point. For the detached leaflet assay, four leaflets from each replicate were sampled while three leaflets were sampled and pooled from each replicate for RNA extraction. Secreted protein fraction or secretome were obtained by pooling five leaflets from each replicate. Six leaflets were pooled from each replicate for phosphate and phosphite measurements.

Transgenic lines [56],[57] impaired in salicylic acid (SA) production (NAHGD2) and jasmonic acid (JA) RNAi silenced transgenes of allene oxidase cyclase gene (AOCZ4), 12-Oxophytodienoate reductase 3 gene (OPR3A5, OPR3Z2) and Coronatine insensitive gene (COIX5) were cultivated the same way as the control (cv. Desiree).

P. infestans infection assay

P. infestans strain 88069 (kindly provided by Francine Grovers, Wageningen University) was used for detached leaflet and whole plant infection assays performed as in Ali, et al. [58]. For the detached leaf assay two 20 μl (15,000 zoospores/ml) drops were spotted on the abaxial side of the leaflet. Lesion size for the detached leaflet assay was measured seven days post inoculation. For the whole plant infection assay comparing phosphite treated transgenic lines and background control (cv. Desiree), the plants were either sprayed with tap water or 36 mM Proalexin, and 24 h after treatment plants were sprayed with a solution containing 15,000 zoospores/ml. Percentage infection was scored visually as percentage of infected leaf area.

Phosphite and phosphate measurements

For the 3, 6, 24 and 120 h time points, phosphate and phosphite were extracted from approximately 40 mg freeze-dried leaf tissue with 1% (v/v) acetic acid at a 1:10 mass:volume ratio using a tissue lyser (Qiagen Australia Pty. Ltd., Doncaster, Australia) for 1 min at 25 hz. After clearing the lysate by centrifugation at 14,000 g and 4°C for 15 min, the supernatant was assayed for phosphate and phosphite as described [59],[60]. Statistically significant differences between phosphite and phosphate levels was determined using ANOVA (general linear model) with leaflet sample type (phosphite treated/“Covered leaves”/ “washed”) and time-point as fixed factors in Minitab 16.

RNA extraction, microarray analysis and verification by qPCR

RNA was isolated with Qiagen RNeasy mini kits (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. RNA quality was determined on a Nanodrop spectrophotometer (Thermo, Saveen Werner, Malmö, Sweden) and integrity was tested on a BioRad Experion (BioRad, Herecules, CA) before further analysis. For mRNA expression analysis a custom-made expression array (Agilent JHI Solanum tuberosum 60 k v1) based on the predicted transcripts in the Solanum phureja genome (version 3.4) was run according to the supplier’s (Agilent) instructions. The complete microarray design is available in ArrayExpress (A-MEXP-2272).

The resulting probe intensities were background corrected and normalised using the quantile method in the Limma R-package (Smyth, Gordon K. 2004). Fold changes and standard errors were obtained by fitting a linear model to each gene and standard errors were smoothed by empirical Bayes. Genes with p < 0.05 (Benjamini-Hochberg adjusted) were regarded as statistically significant. The microarray data was deposited in the ArrayExpress database (accession number E-MTAB-2243). Venn diagrams were drawn using BioVenn [61]. Gene ontology (GO) terms for probes were constructed by clustering 26 plant genomes using a parallelized version of OrthoMCL according [28]. GO term enrichment was performed using GOEast [62] and gene ontology (GO) terms with p < 0.05 (Benjamini-Yekutieli adjusted) were regarded as significantly enriched, clusters from significantly enriched GO term lists were created in ReviGO [27] with default semantic cluster settings (SimRel allowed similarity = 0.7). Cytoscape (v 2.8.3) [63] was used for visualization of GO term clusters. Using Qlucore Omics explorer v 2.2 (Qlucore AB, Lund, Sweden) with variance filter adjusted to 0.25, an unsupervised principal component analysis of gene expression data produced with the Agilent JHI Solanum tuberosum 60 k v1 microarray from 10 mM BABA after 48 h (cv. Desiree; [28]), uninfected late blight resistant clones (cv. Sarpo Mira and SW93-1015), uninfected late blight susceptible cultivar (cv. Desiree) [64] was performed. Expression levels of significantly differentially expressed transcripts 48 h after 10 mM BABA treatment and 3, 6 and 11 h after Proalexin treatment were compared using a linear regression model in ggplot2 in R.

The expression levels of eight genes (StPEN1, StAOS, StFAD3, StMLO1, StWRKY8, StLX-3R, StWRKY1, StNOD084) were validated by qPCR. Primers were designed with the help of Primer 3 [65] according to criteria which include a predicted melting temperature of 57-62°C, a primer length of 18-24 nucleotides, a product size of 100-250 base pairs (bp) and a GC content of 30-70%. Primer sequences are given in Additional file 10: Table S4. For cDNA synthesis 500 ng of total RNA was transcribed to cDNA using SuperScript® III Reverse Transcriptase including degradation with RNase H according to the manufacturer’s protocol (InvitroGen). qPCR was performed with a CFX96 (ABI) using Power SYBR® Master Mix (InvitroGen) and PCR cycles ran according to the manufacturer’s recommendations. The comparative CT method was used for relative quantification of transcripts [66]. A high agreement between microarrays and qPCR was found (R = 0.9; Additional file 8: Figure S6).

Secretome isolation and mass spectrometry of secreted proteins

Secretome isolation was performed using vacuum infiltration with phosphate buffer saline according to a previously described protocol [58]. The secreted protein fraction was dissolved in 6x SDS-PAGE buffer containing DTT, and denatured at 95°C for 3 minutes. Each of the pooled samples (30 μl) was loaded on polyacrylamide gels and separated for 2°Cm with SDS-PAGE. After staining with Coomassie, the gel lane from each sample was cut into about 1 mM2 pieces. Samples were then subjected to in-gel tryptic digestion with incubation (modified sequencing grade; Promega, Madison, WI, USA) overnight at 37°C. Peptides were extracted in 50-80% acetonitrile and excess acetonitrile was vapourised using centrifugation under vacuum. De-salting was performed using UltraMicro spin columns (Nest group).

LC-MS/MS analysis was performed on a LTQ Orbitrap XL ETD with an Eksigent nano-LC system (Eksigent technologies, Dublin, CA, USA). A 5 μl sample was injected and separated at a flow rate of 300 nl/min with a 90 minute gradient. The four most intense ions were selected in data-dependent mode and fragmented in the linear ion trap, with settings as in Ali, et al. [58]. Files were converted to mzML [67] and Mascot Generic Format (MGF) using ProteoWizard [68] and uploaded to the Proteios Software Environment, ProSE [69]. MGF files were used for MS/MS identification, and mzML files for feature detection using msInspect [70]. Peptide and protein identification were performed in Mascot (http://www.matrixscience.com) and X!Tandem (http://thegpm.org/tandem/) in a database consisting of all Solanum proteins in UniProt (http://www.uniprot.org) and all annotated proteins from the potato genome project [71], extended with an equal amount of decoy (reverse sequence) proteins for false discovery rate (FDR) estimation. The MS mass tolerance was set to 5 ppm and MS/MS fragment tolerance to 0.5 Da, with one potential missed cleavage allowed. Cysteine carbamidomethylation was set as fixed and methionine oxidation as variable modification. Peptide cutoff for the combined searched were set to an FDR rate of 1% as described previously [72] within ProSE [69]. Label-free quantification of peptides was performed using a precursor intensity-based strategy [73]. To quantify possible peptides, msInspect [70] feature detection was performed from ProSE using default settings. The features were matched to MS/MS identifications with a retention time tolerance of 0.2 minutes and an m/z tolerance of 0.005 Da as well as a requirement of same charge and LC-MS/MS run. Alignment of peptide features between LC-MS/MS runs was performed within ProSE using the built-in algorithm described by Sandin, et al. [74]. A report of the features corresponding between runs was exported for further analysis. The proteomics data was deposited in PRIDE with the project accession: PXD001031.

Secretome data analysis

Peptides with a FDR of < 0.01 were selected for further analysis. For normalization, we used the Eigen MS method incorporated in DanteR (v0.2) that uses Eigenvalues to find trends in the data for normalization [75],[76]. In DanteR, data was filtered, missing values imputed and an ANOVA was run on the peptide level as before. After Benjamini-Hochberg adjustment, differentially expressed peptides with p < 0.05 were selected for further analysis. The median fold-change of peptides associated to the same proteins was calculated. Functional analysis of identified proteins was performed using MapMan [77]. Signal sequence was predicted with SignalP 4.1 [78] using default parameters.

Availability of supporting data

The microarray data was deposited in ArrayExpress, accession number: E-MTAB-2243. The proteomics data was deposited in PRIDE with the project accession: PXD001031.

Additional files

Abbreviations

BABA:

β-aminobutyric acid

BTH:

S-methyl ester

SA:

Salicylic acid

JA:

Jasmonic acid

ROS:

Reactive oxygen species

“hR-like”:

hypersensitive response like

FW:

Fresh weight

PEN1:

Pentacyclic triterpene synthase

NDR1/NHL25:

Non-race specific disease resistance

BCS1:

Cytochrome BC1 synthesis protein

PAL:

Phenyl ammonia lyase

ERD7:

Early response to dehydration 7

PIP:

Phosphatidylinositol-4-phosphate

PCA:

Principal component analysis

PR1:

Pathogenesis responsive protein 1

AOC:

Allene oxidase cyclase

OPR3:

12-Oxophytodienoate reductase

COI:

Coronatine insensitive

GFP:

Green fluorescent protein

ANOVA:

Analysis of Variance

SDS-PAGE:

Sodium dodecyl sulphate - Polyacrylamide gel electrophoresis

LTQ:

Linear trap quadropole

MS:

Mass spectrometry

LC:

Liquid chromatography

dpi:

Days post infection

GO:

Gene ontology

References

  1. Nowicki M, Foolad MR, Nowakowska M, Kozik EU: Potato and tomato late blight caused by Phytophthora infestans: an overview of pathology and resistance breeding. Plant Dis. 2012, 96 (1): 4-17. 10.1094/PDIS-05-11-0458.

    Article  Google Scholar 

  2. Vleeshouwers VG, Raffaele S, Vossen JH, Champouret N, Oliva R, Segretin ME, Rietman H, Cano LM, Lokossou A, Kessel G, Pel M, Kamoun S: Understanding and exploiting late blight resistance in the age of effectors. Annu Rev Phytopathol. 2011, 49: 507-531. 10.1146/annurev-phyto-072910-095326.

    Article  CAS  PubMed  Google Scholar 

  3. Whisson SC, Avrova AO, Boevink PC, Armstrong MR, Seman ZA, Hein I, Birch PRJ: Exploiting knowledge of pathogen effectors to enhance late blight resistance in potato. Potato Res. 2011, 54 (4): 325-340. 10.1007/s11540-011-9197-y.

    Article  Google Scholar 

  4. Haverkort AJ, Boonekamp PM, Hutten R, Jacobsen E, Lotz LAP, Kessel GJT, Visser RGF, Vossen EAG: Societal costs of late blight in potato and prospects of durable resistance through cisgenic modification. Potato Res. 2008, 51 (1): 47-57. 10.1007/s11540-008-9089-y.

    Article  Google Scholar 

  5. Mishra AK, Sharma K, Misra RS: Elicitor recognition, signal transduction and induced resistance in plants. J Plant Interact. 2012, 7 (2): 95-120. 10.1080/17429145.2011.597517.

    Article  Google Scholar 

  6. Van Loon LC: Induced resistance in plants and the role of pathogenesis-related proteins. Eur J Plant Pathol. 1997, 103 (9): 753-765. 10.1023/A:1008638109140.

    Article  CAS  Google Scholar 

  7. Ton J, Jakab G, Toquin V, Flors V, Iavicoli A, Maeder M, Métraux J, Mauch-Mani B: Dissecting the beta-aminobutyric acid-induced priming phenomenon in Arabidopsis. Plant Cell. 2005, 17 (3): 987-999. 10.1105/tpc.104.029728.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Eschen-Lippold L, Altmann S, Rosahl S: dl-β-aminobutyric acid-induced resistance of potato against Phytophthora infestans requires salicylic acid but not oxylipins. Mol Plant-Microbe Interact. 2010, 23 (5): 585-592. 10.1094/MPMI-23-5-0585.

    Article  CAS  PubMed  Google Scholar 

  9. Bengtsson T, Holefors A, Witzell J, Andreasson E, Liljeroth E: Activation of defence responses to Phytophthora infestans in potato by BABA. Plant Pathol. 2014, 63 (1): 193-202. 10.1111/ppa.12069.

    Article  CAS  Google Scholar 

  10. McDonald AE, Grant BR, Plaxton WC: Phosphite (Phosphorous Acid): its relevance in the environment and agriculture and influence on plant phosphate starvation response. J Plant Nutr. 2001, 24 (10): 1505-1519. 10.1081/PLN-100106017.

    Article  CAS  Google Scholar 

  11. King M, Reeve W, Van der Hoek MB, Williams N, McComb J, O’Brien PA, Hardy GE: Defining the phosphite-regulated transcriptome of the plant pathogen Phytophthora cinnamomi. Mol Genet Genomics. 2010, 284 (6): 425-435. 10.1007/s00438-010-0579-7.

    Article  CAS  PubMed  Google Scholar 

  12. Lobato M, Olivieri F, Daleo G, Andreu A: Antimicrobial activity of phosphites against different potato pathogens. J Plant Dis Prot. 2010, 3 (117): 102-109.

    Google Scholar 

  13. Smillie R, Grant B, Guest D: The mode of action of phosphite: evidence for both direct and indirect modes of action on three Phytophthora spp. in plants. Phytopathology. 1989, 79 (9): 921-926. 10.1094/Phyto-79-921.

    Article  CAS  Google Scholar 

  14. Thao HTB, Yamakawa T: Phosphite (phosphorous acid): fungicide, fertilizer or bio-stimulator?. Soil Sci Plant Nutr. 2009, 55 (2): 228-234. 10.1111/j.1747-0765.2009.00365.x.

    Article  CAS  Google Scholar 

  15. Olivieri FP, Feldman ML, Machinandiarena MF, Lobato MC, Caldiz DO, Daleo GR, Andreu AB: Phosphite applications induce molecular modifications in potato tuber periderm and cortex that enhance resistance to pathogens. Crop Prot. 2012, 32: 1-6. 10.1016/j.cropro.2011.08.025.

    Article  CAS  Google Scholar 

  16. Machinandiarena MF, Lobato MC, Feldman ML, Daleo GR, Andreu AB: Potassium phosphite primes defense responses in potato against Phytophthora infestans. J Plant Physiol. 2012, 169 (14): 1417-1424. 10.1016/j.jplph.2012.05.005.

    Article  CAS  PubMed  Google Scholar 

  17. Kromann P, Pérez WG, Taipe A, Schulte-Geldermann E, Sharma BP, Andrade-Piedra JL, Forbes GA: Use of phosphonate to manage foliar potato late blight in developing countries. Plant Dis. 2012, 96 (7): 1008-1015. 10.1094/PDIS-12-11-1029-RE.

    Article  CAS  Google Scholar 

  18. Conrath U, Klessig DF, Bachmair A: Tobacco plants perturbed in the ubiquitin-dependent protein degradation system accumulate callose, salicyclic acid, and pathogenesis-related protein 1. Plant Cell Rep. 1998, 17 (11): 876-880. 10.1007/s002990050501.

    Article  CAS  Google Scholar 

  19. Danova-Alt R, Dijkema C, DE Waard P, Kock M: Transport and compartmentation of phosphite in higher plant cells--kinetic and P nuclear magnetic resonance studies. Plant Cell Environ. 2008, 31 (10): 1510-1521. 10.1111/j.1365-3040.2008.01861.x.

    Article  CAS  PubMed  Google Scholar 

  20. Berkowitz O, Jost R, Kollehn DO, Fenske R, Finnegan PM, O’Brien PA, Hardy GE, Lambers H: Acclimation responses of Arabidopsis thaliana to sustained phosphite treatments. J Exp Bot. 2013, 64 (6): 1731-1743. 10.1093/jxb/ert037.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Ticconi CA, Delatorre CA, Abel S: Attenuation of phosphate starvation responses by phosphite in Arabidopsis. Plant Physiol. 2001, 127 (3): 963-972. 10.1104/pp.010396.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Massoud K, Barchietto T, Le Rudulier T, Pallandre L, Didierlaurent L, Garmier M, Ambard-Bretteville F, Seng JM, Saindrenan P: Dissecting phosphite-induced priming in Arabidopsis infected with Hyaloperonospora arabidopsidis. Plant Physiol. 2012, 159 (1): 286-298. 10.1104/pp.112.194647.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Eshraghi L, Anderson J, Aryamanesh N, Shearer B, McComb J, Hardy GES, O’Brien PA: Phosphite primed defence responses and enhanced expression of defence genes in Arabidopsis thaliana infected with Phytophthora cinnamomi. Plant Pathol. 2011, 60 (6): 1086-1095. 10.1111/j.1365-3059.2011.02471.x.

    Article  CAS  Google Scholar 

  24. Lim S, Borza T, Peters RD, Coffin RH, Al-Mughrabi KI, Pinto DM, Wang-Pruski G: Proteomics analysis suggests broad functional changes in potato leaves triggered by phosphites and a complex indirect mode of action against Phytophthora infestans. J Proteomic. 2013, 93: 207-223. 10.1016/j.jprot.2013.03.010.

    Article  CAS  Google Scholar 

  25. Alexandersson E, Ashfaq A, Resjö S, Andreasson E: Plant secretome proteomics. Front Plant Sci. 2013, 4: 9-10.3389/fpls.2013.00009.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Borza T, Schofield A, Sakthivel G, Bergese J, Gao X, Rand J, Wang-Pruski G: Ion chromatography analysis of phosphite uptake and translocation by potato plants: dose-dependent uptake and inhibition of Phytophthora infestans development. Crop Prot. 2014, 56: 74-81. 10.1016/j.cropro.2013.10.024.

    Article  CAS  Google Scholar 

  27. Supek F, Bošnjak M, Škunca N, Šmuc T: REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011, 6 (7): e21800-10.1371/journal.pone.0021800.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Bengtsson T, Weighill D, Proux-Wera E, Levander F, Resjo S, Burra DD, Moushib LI, Hedley PE, Liljeroth E, Jacobson D, Alexandersson E, Andreasson E: Proteomics and transcriptomics of the BABA-induced resistance response in potato using a novel functional annotation approach. BMC Genomics. 2014, 15 (1): 315-10.1186/1471-2164-15-315.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Halim VA, Eschen-Lippold L, Altmann S, Birschwilks M, Scheel D, Rosahl S: Salicylic acid is important for basal defense of Solanum tuberosum against Phytophthora infestans. Mol Plant Microbe Interact. 2007, 20 (11): 1346-1352. 10.1094/MPMI-20-11-1346.

    Article  CAS  PubMed  Google Scholar 

  30. Thao HTB, Yamakawa T, Myint AK, Sarr PS: Effects of phosphite, a reduced form of phosphate, on the growth and phosphorus nutrition of spinach (Spinacia oleraceaL.). Soil Sci Plant Nutr. 2008, 54 (5): 761-768. 10.1111/j.1747-0765.2008.00290.x.

    Article  CAS  Google Scholar 

  31. Thao HTB, Yamakawa T: Phosphate absorption of intact komatsuna plants as influenced by phosphite. Soil Sci Plant Nutr. 2010, 56 (1): 133-139. 10.1111/j.1747-0765.2009.00431.x.

    Article  CAS  Google Scholar 

  32. Ávila FW, Faquin V, Araújo JL, Marques DJ, Ribeiro Júnior PM, Lobato AKDS, Ramos SJ, Baliza DP: Phosphite supply affects phosphorus nutrition and biochemical response in maize plants. Aust J Crop Sci. 2011, 5 (6): 646-653.

    Google Scholar 

  33. Varadarajan DK, Karthikeyan AS, Matilda PD, Raghothama KG: Phosphite, an analog of phosphate, suppresses the coordinated expression of genes under phosphate starvation. Plant Physiol. 2002, 129 (3): 1232-1240. 10.1104/pp.010835.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Lager I, Andreasson O, Dunbar TL, Andreasson E, Escobar MA, Rasmusson AG: Changes in external pH rapidly alter plant gene expression and modulate auxin and elicitor responses. Plant Cell Environ. 2010, 33 (9): 1513-1528.

    PubMed Central  CAS  PubMed  Google Scholar 

  35. Lee J, Nam J, Park HC, Na G, Miura K, Jin JB, Yoo CY, Baek D, Kim DH, Jeong JC, Kim D, Lee SY, Salt DE, Mengiste T, Gong Q, Ma S, Bohnert HJ, Kwak SS, Bressan RA, Hasgawa PM, Yun DJ: Salicylic acid-mediated innate immunity in Arabidopsis is regulated by SIZ1 SUMO E3 ligase. Plant J. 2007, 49 (1): 79-90. 10.1111/j.1365-313X.2006.02947.x.

    Article  CAS  PubMed  Google Scholar 

  36. Ho LH, Giraud E, Uggalla V, Lister R, Clifton R, Glen A, Thirkettle-Watts D, Van Aken O, Whelan J: Identification of regulatory pathways controlling gene expression of stress-responsive mitochondrial proteins in Arabidopsis. Plant Physiol. 2008, 147 (4): 1858-1873. 10.1104/pp.108.121384.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Studham ME, MacIntosh GC: Phytohormone signaling pathway analysis method for comparing hormone responses in plant-pest interactions. BMC Res Notes. 2012, 5: 392-10.1186/1756-0500-5-392.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Kohorn BD, Kohorn SL: The cell wall-associated kinases, WAKs, as pectin receptors. Front Plant Sci. 2012, 3: 88-10.3389/fpls.2012.00088.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Kodaira KS, Qin F, Tran LS, Maruyama K, Kidokoro S, Fujita Y, Shinozaki K, Yamaguchi-Shinozaki K: Arabidopsis Cys2/His2 zinc-finger proteins AZF1 and AZF2 negatively regulate abscisic acid-repressive and auxin-inducible genes under abiotic stress conditions. Plant Physiol. 2011, 157 (2): 742-756. 10.1104/pp.111.182683.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Bolton MD: Primary metabolism and plant defense-fuel for the fire. Mol Plant-Microbe Interact. 2009, 22 (5): 487-497. 10.1094/MPMI-22-5-0487.

    Article  CAS  PubMed  Google Scholar 

  41. Arasimowicz-Jelonek M, Kosmala A, Janus L, Abramowski D, Floryszak-Wieczorek J: The proteome response of potato leaves to priming agents and S-nitrosoglutathione. Plant Sci. 2013, 198: 83-90. 10.1016/j.plantsci.2012.10.004.

    Article  CAS  PubMed  Google Scholar 

  42. Restrepo S, Myers KL, del Pozo O, Martin GB, Hart AL, Buell CR, Fry WE, Smart CD: Gene profiling of a compatible interaction between phytophthora infestans and Solanum tuberosum suggests a role for carbonic anhydrase. Mol Plant-Microbe Interact. 2005, 18 (9): 913-922. 10.1094/MPMI-18-0913.

    Article  CAS  PubMed  Google Scholar 

  43. Hernandez-Blanco C, Feng DX, Hu J, Sanchez-Vallet A, Deslandes L, Llorente F, Berrocal-Lobo M, Keller H, Barlet X, Sanchez-Rodriguez C, Anderson LK, Somerville S, Marco Y, Molina A: Impairment of cellulose synthases required for Arabidopsis secondary cell wall formation enhances disease resistance. Plant Cell. 2007, 19 (3): 890-903. 10.1105/tpc.106.048058.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Yaeno T, Li H, Chaparro-Garcia A, Schornack S, Koshiba S, Watanabe S, Kigawa T, Kamoun S, Shirasu K: Phosphatidylinositol monophosphate-binding interface in the oomycete RXLR effector AVR3a is required for its stability in host cells to modulate plant immunity. Proc Natl Acad Sci. 2011, 108 (35): 14682-14687. 10.1073/pnas.1106002108.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Ma S, Bohnert HJ: Integration of Arabidopsis thaliana stress-related transcript profiles, promoter structures, and cell-specific expression. Genome Biol. 2007, 8 (4): R49-10.1186/gb-2007-8-4-r49.

    Article  PubMed Central  PubMed  Google Scholar 

  46. Nishimura MT, Stein M, Hou B-H, Vogel JP, Edwards H, Somerville SC: Loss of a callose synthase results in salicylic acid-dependent disease resistance. Science. 2003, 301 (5635): 969-972. 10.1126/science.1086716.

    Article  CAS  PubMed  Google Scholar 

  47. Chen X-Y, Kim J-Y: Callose synthesis in higher plants. Plant Signal Behav. 2009, 4 (6): 489-492. 10.4161/psb.4.6.8359.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Xia Y, Suzuki H, Borevitz J, Blount J, Guo Z, Patel K, Dixon RA, Lamb C: An extracellular aspartic protease functions in Arabidopsis disease resistance signaling. EMBO J. 2004, 23 (4): 980-988. 10.1038/sj.emboj.7600086.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Hwang IS, Hwang BK: The pepper mannose-binding lectin gene CaMBL1 is required to regulate cell death and defense responses to microbial pathogens. Plant Physiol. 2011, 155 (1): 447-463. 10.1104/pp.110.164848.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. De Lorenzo G, D’Ovidio R, Cervone F: The role of Polygalacturonase-Inhibiting Proteins (PGIPS) in defense against pathogenic fungi. Annu Rev Phytopathol. 2001, 39 (1): 313-335. 10.1146/annurev.phyto.39.1.313.

    Article  CAS  PubMed  Google Scholar 

  51. Alexandersson E, Becker JV, Jacobson D, Nguema-Ona E, Steyn C, Denby KJ, Vivier MA: Constitutive expression of a grapevine polygalacturonase-inhibiting protein affects gene expression and cell wall properties in uninfected tobacco. BMC Res Notes. 2011, 4: 493-10.1186/1756-0500-4-493.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Grover A: Plant chitinases: genetic diversity and physiological roles. Crit Rev Plant Sci. 2012, 31 (1): 57-73. 10.1080/07352689.2011.616043.

    Article  CAS  Google Scholar 

  53. Deliopoulos T, Kettlewell PS, Hare MC: Fungal disease suppression by inorganic salts: a review. Crop Prot. 2010, 29 (10): 1059-1075. 10.1016/j.cropro.2010.05.011.

    Article  CAS  Google Scholar 

  54. Lopez-Arredondo DL, Herrera-Estrella L: Engineering phosphorus metabolism in plants to produce a dual fertilization and weed control system. Nat Biotechnol. 2012, 30 (9): 889-893. 10.1038/nbt.2346.

    Article  CAS  PubMed  Google Scholar 

  55. Pinto KMS, Do Nascimento LC, De Souza Gomes EC, da Silva HF, Dos Reis Miranda J: Efficiency of resistance elicitors in the management of grapevine downy mildew Plasmopara viticola: epidemiological, biochemical and economic aspects. Eur J Plant Pathol. 2012, 134 (4): 745-754. 10.1007/s10658-012-0050-1.

    Article  Google Scholar 

  56. Halim VA, Altmann S, Ellinger D, Eschen-Lippold L, Miersch O, Scheel D, Rosahl S: PAMP-induced defense responses in potato require both salicylic acid and jasmonic acid. Plant J. 2009, 57 (2): 230-242. 10.1111/j.1365-313X.2008.03688.x.

    Article  CAS  PubMed  Google Scholar 

  57. Halim VA, Hunger A, Macioszek V, Landgraf P, Nürnberger T, Scheel D, Rosahl S: The oligopeptide elicitor Pep-13 induces salicylic acid-dependent and -independent defense reactions in potato. Physiol Mol Plant Pathol. 2004, 64 (6): 311-318. 10.1016/j.pmpp.2004.10.003.

    Article  CAS  Google Scholar 

  58. Ali A, Moushib LI, Lenman M, Levander F, Olsson K, Carlson-Nilson U, Zoteyeva N, Liljeroth E, Andreasson E: Paranoid potato: phytophthora-resistant genotype shows constitutively activated defense. Plant Signal Behav. 2012, 7 (3): 400-408. 10.4161/psb.19149.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  59. Ames BN: Assay of inorganic phosphate, total phosphate and phosphatase. Methods Enzymol. 1966, 8: 115-118. 10.1016/0076-6879(66)08014-5.

    Article  CAS  Google Scholar 

  60. Berkowitz O, Jost R, Pearse SJ, Lambers H, Finnegan PM, Hardy GE, O’Brien PA: An enzymatic fluorescent assay for the quantification of phosphite in a microtiter plate format. Anal Biochem. 2011, 412 (1): 74-78. 10.1016/j.ab.2011.01.014.

    Article  CAS  PubMed  Google Scholar 

  61. Hulsen T, de Vlieg J, Alkema W: BioVenn-a web application for the comparison and visualization of biological lists using area-proportional Venn diagrams. BMC Genomics. 2008, 9 (1): 488-10.1186/1471-2164-9-488.

    Article  PubMed Central  PubMed  Google Scholar 

  62. Zheng Q, Wang XJ: GOEAST: a web-based software toolkit for Gene Ontology enrichment analysis. Nucleic Acids Res. 2008, 36 (Web Server issue): W358-W363. 10.1093/nar/gkn276.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  63. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13 (11): 2498-2504. 10.1101/gr.1239303.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  64. Ali A, Alexandersson E, Sandin M, Resjö S, Lenman M, Hedley P, Levander F, Andreasson E: Quantitative proteomics and transcriptomics of potato in response to Phytophthora infestans in compatible and incompatible interactions. BMC Genomics. 2014, 15 (1): 497-10.1186/1471-2164-15-497.

    Article  PubMed Central  PubMed  Google Scholar 

  65. Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Meth Protoc Springer. 1999, 6: 365-386. 10.1385/1-59259-192-2:365.

    Article  Google Scholar 

  66. Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008, 3 (6): 1101-1108. 10.1038/nprot.2008.73.

    Article  CAS  PubMed  Google Scholar 

  67. Martens L, Chambers M, Sturm M, Kessner D, Levander F, Shofstahl J, Tang WH, Römpp A, Neumann S, Pizarro AD, Palazzi LM, Tasman N, Coleman M, Reisinger F, Souda P, Hermjakob H, Binz PA, Deutsch PA: mzML-a community standard for mass spectrometry data. Mol Cell Proteomics. 2011, 10 (1): R110-000133-10.1074/mcp.R110.000133.

    Article  PubMed Central  PubMed  Google Scholar 

  68. Kessner D, Chambers M, Burke R, Agus D, Mallick P: ProteoWizard: open source software for rapid proteomics tools development. Bioinformatics. 2008, 24 (21): 2534-2536. 10.1093/bioinformatics/btn323.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  69. Häkkinen J, Vincic G, Månsson O, Wårell K, Levander F: The proteios software environment: an extensible multiuser platform for management and analysis of proteomics data. J Proteome Res. 2009, 8 (6): 3037-3043. 10.1021/pr900189c.

    Article  PubMed  Google Scholar 

  70. Bellew M, Coram M, Fitzgibbon M, Igra M, Randolph T, Wang P, May D, Eng J, Fang R, Lin C, Chen J, Goodlett D, Whiteaker J, Paulovich A, McIntosh M: A suite of algorithms for the comprehensive analysis of complex protein mixtures using high-resolution LC-MS. Bioinformatics. 2006, 22 (15): 1902-1909. 10.1093/bioinformatics/btl276.

    Article  CAS  PubMed  Google Scholar 

  71. Genome sequence: and analysis of the tuber crop potato. Nature. 2011, 475 (7355): 189-195. 10.1038/nature10158.

  72. Levander F, Krogh M, Wårell K, Gärdén P, James P, Häkkinen J: Automated reporting from gel-based proteomics experiments using the open source Proteios database application. Proteomics. 2007, 7 (5): 668-674. 10.1002/pmic.200600814.

    Article  CAS  PubMed  Google Scholar 

  73. Sandin M, Teleman J, Malmström J, Levander F: Data processing methods and quality control strategies for label-free LC-MS protein quantification. Biochim Biophys Acta. 2014, 1844 (1): 29-41. 10.1016/j.bbapap.2013.03.026.

    Article  CAS  PubMed  Google Scholar 

  74. Sandin M, Ali A, Hansson K, Mansson O, Andreasson E, Resjo S, Levander F: An adaptive alignment algorithm for quality-controlled label-free LC-MS. Mol Cell Proteomics. 2013, 12 (5): 1407-1420. 10.1074/mcp.O112.021907.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  75. Karpievitch YV, Taverner T, Adkins JN, Callister SJ, Anderson GA, Smith RD, Dabney AR: Normalization of peak intensities in bottom-up MS-based proteomics using singular value decomposition. Bioinformatics. 2009, 25 (19): 2573-2580. 10.1093/bioinformatics/btp426.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  76. Taverner T, Karpievitch YV, Polpitiya AD, Brown JN, Dabney AR, Anderson GA, Smith RD: DanteR: an extensible R-based tool for quantitative analysis of -omics data. Bioinformatics. 2012, 28 (18): 2404-2406. 10.1093/bioinformatics/bts449.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  77. Thimm O, Bläsing O, Gibon Y, Nagel A, Meyer S, Krüger P, Selbig J, Müller LA, Rhee SY, Stitt M: mapman: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004, 37 (6): 914-939. 10.1111/j.1365-313X.2004.02016.x.

    Article  CAS  PubMed  Google Scholar 

  78. Petersen TN, Brunak S, von Heijne G, Nielsen H: SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Meth. 2011, 8 (10): 785-786. 10.1038/nmeth.1701.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

This work was supported by Swedish Foundation for Strategic Research, the Crafoord foundation, Mistra Biotech, Swedish Research Council Formas and the Australian Research Council Linkage Project (LP0776252) and its partners. We thank Estelle Proux-Wéra (PlantLink) for bioinformatics support. We thank Sabine Rosahl for providing the transgenic potato lines, Francine Grovers for contributing with the P. infestans strain, Marit Lenman for qPCR primer design and Mia Mogren for excellent assistance in the laboratory.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Erik Alexandersson.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

DDB, EOA, EL, EA planned the experiments and wrote the paper. OB performed the enzymatic assay to measure phosphite. PH and JM performed the microarray analysis. SR and FL performed the proteomics and helped with the analysis. DDB and EOA performed the other lab experiments and did the bioinformatics analysis. All authors read and approved the final manuscript.

Electronic supplementary material

12870_2014_254_MOESM1_ESM.pptx

Additional file 1: Figure S1.: (A) “HR like” symptoms observed at the site of P. infestans infestation in phosphite treated leaflets (top panel) while extensive sporulation of P. Infestans observed on water treated leaflets (bottom panel) 7 dpi in the detached leaflet assay, (B) detached leaflet assay of Water sprayed (control) and acidified water sprayed leaflets. (PPTX 2 MB)

12870_2014_254_MOESM2_ESM.docx

Additional file 2: Figure S2.: Correlation between microarrays and qPCR for eight different genes specified in Additional file 10: Table S4. (DOCX 16 KB)

12870_2014_254_MOESM3_ESM.docx

Additional file 3: Table S1.: The. List of “Core” of phosphite induced transcripts that are differentially expressed at 3, 6 and 11 h post treatment. (DOCX 28 KB)

12870_2014_254_MOESM4_ESM.pptx

Additional file 4: Figure S3.: Selection of significantly enriched gene ontology terms at each time point (Benjamini-Yekutieli adj. p-value <0.05). (PPTX 162 KB)

12870_2014_254_MOESM5_ESM.docx

Additional file 5: Table S2.: List of significantly enriched gene ontology terms among transcripts regulated both by BABA 48 h after treatment and phosphite 3 h after treatment. (DOCX 31 KB)

12870_2014_254_MOESM6_ESM.pptx

Additional file 6: Figure S4.: Principal component analysis of transcriptomic data obtained from 3 biological replicates of BABA treated plants (cv. Desiree, 10 mM BABA 48 h post treatment) referred to as BABA, phosphite treated plants (cv. Desiree, 36 mM phosphite) sampled 3 (Phi3), 6 (Phi6) and 11 (Phi11) h after phosphite treatment, uninfected late blight resistant clones [58] cv Sarpo Mira (SMC) and clone SW93-1015 (SWC) and uninfected late blight susceptible clone cv. Desiree (DesC). (PPTX 42 KB)

12870_2014_254_MOESM7_ESM.pptx

Additional file 7: Figure S5.: Functional analysis using MapMan of secreted proteins with significantly changed abundance detected 48 h after phosphite treatment [77]. (PPTX 88 KB)

12870_2014_254_MOESM8_ESM.pptx

Additional file 8: Figure S6.: Mapman visualization of metabolic processes in transcripts significantly changing 6 and 11 h after treatment. (PPTX 209 KB)

12870_2014_254_MOESM9_ESM.docx

Additional file 9: Table S3.: Significant differences observed in transcripts associated with salicylic and jasmonic acid pathways, the markers were selected as suggested in Studham, et al. [37]. (DOCX 20 KB)

Additional file 10: Table S4.: Primer sequences for qPCR. (DOCX 14 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Burra, D.D., Berkowitz, O., Hedley, P.E. et al. Phosphite-induced changes of the transcriptome and secretome in Solanum tuberosum leading to resistance against Phytophthora infestans. BMC Plant Biol 14, 254 (2014). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-014-0254-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-014-0254-y

Keywords