Water Stress and Aphid Feeding Differentially Influence Metabolite Composition in Arabidopsis thaliana (L.)

Abstract

Little is known about how drought stress influences plant secondary metabolite accumulation and how this affects plant defense against different aphids. We therefore cultivated Arabidopsis thaliana (L.) plants under well-watered, drought, and water-logged conditions. Two aphid species were selected for this study: the generalist Myzus persicae (Sulzer) and the crucifer specialist Brevicoryne brassicae (L.). Metabolite concentrations in the phloem sap, which influence aphid growth, changed particularly under drought stress. Levels of sucrose and several amino acids, such as glutamic acid, proline, isoleucine, and lysine increased, while concentrations of 4-methoxyindol-3-ylmethyl glucosinolate decreased. M. persicae population growth was highest on plants under drought stress conditions. However, B. brassicae did not profit from improved phloem sap quality under drought stress and performed equally in all water treatments. Water stress and aphids generally had an opposite effect on the accumulation of secondary metabolites in the plant rosettes. Drought stress and water-logging led to increased aliphatic glucosinolate and flavonoid levels. Conversely, aphid feeding, especially of M. persicae, reduced levels of flavonoids and glucosinolates in the plants. Correspondingly, transcript levels of aliphatic biosynthetic genes decreased after feeding of both aphid species. Contrary to M. persicae, drought stress did not promote population growth of B. brassicae on these plants. The specialist aphid induced expression of CYP79B2, CYP79B3, and PAD3 with corresponding accumulation of indolyl glucosinolates and camalexin. This was distinct from M. persicae, which did not elicit similarly strong camalexin accumulation, which led to the hypothesis of a specific defense adaptations against the specialist aphid.

About the Authors:

Inga Mewis

* E-mail: [email protected]

Affiliation: Department of Quality, Leibniz-Institute of Vegetable and Ornamental Crops Großbeeren/Erfurt e.V., Großbeeren, Germany

Mohammed A. M. Khan

Current address: Department of Entomology, Faculty of Agriculture, Bangladesh Agricultural University, Mymensingh, Bangladesh

Affiliation: Division Urban Plant Ecophysiology, Institute of Horticultural Sciences, Humboldt-Universität zu Berlin, Berlin, Germany

Erich Glawischnig

Affiliation: Lehrstuhl für Genetik, Technische Universität München, Freising, Germany

Monika Schreiner

Affiliation: Department of Quality, Leibniz-Institute of Vegetable and Ornamental Crops Großbeeren/Erfurt e.V., Großbeeren, Germany

Christian Ulrichs

Affiliation: Division Urban Plant Ecophysiology, Institute of Horticultural Sciences, Humboldt-Universität zu Berlin, Berlin, Germany

Introduction

In their natural growing environment, plants encounter multiple environmental stresses that can alter both their chemical composition and the associated herbivore community. Drought is a major environmental stress that not only directly influences plant physiological processes leading to decreased growth, but also changes the allocation of resources and profile of secondary metabolites [1], [2], [3]. Enhanced levels of soluble proteins [1], [4], soluble carbohydrates [5], [6], and free amino acid [7], [8] in plants caused by drought may favor the population growth of phloem-feeding insects. Increased population size on drought-stressed plants is observed for some aphid species [2], [3], [9], [10]. Population growth of the generalist Myzus persicae (Sulzer) was more rapid on drought-stressed plants; however, the crucifer specialist Brevicoryne brassicae (L.) performed equally well on drought and well-watered Brassica oleracea var. italica (broccoli) plants [2], [3].

Besides the general nutritional quality of phloem sap (primary metabolites), defensive compounds of host plants such as glucosinolate may be responsible for such a variation in the aphids’ performance on drought-stressed plants. Decreased levels of glucosinolates in drought-stressed brassicaceous crops in connection with the improved performance of generalist aphids has been reported in previous studies [2], [3], [11]. Unlike generalist aphids such as M. persicae [12], crucifer specialists such as B. brassicae may tolerate to some extent glucosinolates in their host plants. However, the influence of drought stress on secondary metabolite accumulation in relation to the impact this has on plant resistance against aphids with different feeding specialization is still not well understood.

The glucosinolate-myrosinase system is the characteristic chemical defense mechanism against herbivores in plants of the order Brassicales. Tissue damage, such as during herbivory, brings glucosinolates in contact with their hydrolyzing enzymes (myrosinases), which are localized in different cells or compartments [13], [14]. Different biologically active compounds such as isothiocyanates and nitriles are released, which are deterrents to many herbivore insects [15]. Aphids are able to avoid the glucosinolate-myrosinase defense in plants since they only inflict minor wounding on the plant while feeding. Intact glucosinolate, but no corresponding hydrolyses product, has been detected in the honey-dew of aphids [16]. Therefore, glucosinolates in phloem sap, and not their hydrolysis products, may provide a defense against aphids. Although aphids inflict slight wounds during feeding, elicitors released by aphids can activate plants’ responses and modify the glucosinolate profile of host plants [3], [17], [18]. A recent study with broccoli plants has shown that the availability of soil water changes the plant’s glucosinolate response to aphid feeding. Glucosinolate levels in water-stressed plants decreased after M. persicae feeding, but remained unchanged after feeding by B. brassicae [3].

In A. thaliana, the glucosinolate biosynthetic pathway has been well studied and most of the associated genes encoding the required enzymes have been isolated and characterized (Figure 1) [15], [19]. Glucosinolates in plants are synthesized from amino acids and their type depends on the precursor amino acid. The aliphatic glucosinolates are mainly synthesized from methionine, and indolyl glucosinolates from tryptophan. In aliphatic glucosinolate biosynthesis, methionine undergoes chain elongation before entering the core structure pathway. Methionine is initially deaminated by BCAT4 to 2-oxo acid, which enters a cycle of three successive transformations to add a single methylene group, and condensation with acetyl-CoA is catalyzed by MAM1 and MAM3. Chain-elongated methionines are then converted to the corresponding aldoximes by CYP79F1 and CYP79F2, while for the indolyl route CYP79B2 and CYP79B3 accept tryptophan as the substrate. The aldoximes are further metabolized by cytochrome P450 belonging to the CYP83 family (CYP83A1 with a preference for aliphatic and CYP83B1 for indolyl aldoximes) to form S-alkylthiohydroximates with glutathione as the sulfur donor [20]. Thiohydroximates are formed by C-S lyase (SUR1) cleaving S-alkylthiohydroximates. This is followed by glucosylation by glucosyltransferases of the UGT74 family to form desulfo glucosinolates, whereby UGT74C1 glucosylates methionine-derived thiohydroximates; UGT74B1 metabolizes tryptophan-derived thiohydroximates. The final sulfation step is catalyzed by sulfotransferases. Among these, SOT17 and SOT18 metabolize methionine-derived and SOT16 tryptophan-derived precursors. After forming the glucosinolate core structure, secondary modification of the side chain takes place which creates structural biodiversity and changes biological activity [19]. In recent years, transcription factors belonging to the MYB gene family have been identified to regulate glucosinolate biosynthesis. MYB28, MYB29, and MYB76 control aliphatic glucosinolate biosynthesis [19], [21], [22], [23], [24]; MYB34 [25], MYB51, and MYB122 [21] regulate indolyl glucosinolate biosynthesis.

[Figure omitted. See PDF.]

Figure 1. Outline of glucosinolate and camalexin biosynthetic pathways in A. thaliana plants.

(Grey indicates the last steps of the camalexin pathway independent from glucosinolate biosynthesis; large arrows indicate chain elongation of keto-amino acids, each cycle adds one methyl group).

https://doi.org/10.1371/journal.pone.0048661.g001

Cross-talk between the glucosinolate pathway and other metabolic pathways may have an impact on the actual glucosinolate profile produced [26], e.g. since indolyl acetaldoxime serves as a branching point for indolyl glucosinolate and camalexin biosynthesis (Fig. 1) [27], [28]. The phytoalexin camalexin plays a crucial role in the defense against fungal and bacterial pathogens [29] and accumulates at infection sites [30], [31], [32], [33]. A recent study demonstrated that A. thaliana plants accumulate camalexin after infestation with B. brassicae [34]. However, no comparative study has been undertaken to analyze the camalexin response elicited by the feeding of generalist and specialist aphids. Since drought stress alters indolyl glucosinolate metabolism in plants [2], [3], it is possible that biosynthesis of camalexin in A. thaliana plants is also affected by water stress.

In the present study, the model plant Arabidopsis thaliana (L.) ecotype Columbia (Col-0) was used to study the genetic regulation of glucosinolate biosynthesis as a response to water stress and aphid feeding. Three water regimes - well-watered, drought, and water-logged - and two important aphid pests - the specialist B. brassicae and the generalist M. persicae [35] - were used to elucidate the molecular and biochemical response of plants cultivated under different water status conditions to insects. Besides major glucosinolate biosynthetic genes and transcription factors, we analyzed genes related to major plant signaling pathways, and PAD3 involved in camalexin biosynthesis [36], [37]. The main research objective was to investigate factors that can vary the performance of aphid species with different specializations. Unlike Khan et al. (2), who focused on water stress, and its effects on both aphid performance and accumulation of glucosinolates and flavonoids in broccoli, we undertook a wider approach by using A. thaliana and included the molecular level for defensive metabolite biosynthesis. Water stress mediated changes in the nutritional components of phloem sap, such as levels of sugars and free amino acids were analyzed along with flavonoids as secondary metabolites and defensive compounds such as glucosinolate and camalexin.

Materials and Methods

Plant Material and Cultivation

Arabidopsis thaliana (L.) Columbia wild type (Col-0) was used to study the effect of water stress-induced changes in plant metabolites, corresponding gene expression, and aphid performance. Seeds were sown in pots (10 cm×10 cm) filled with potting substrate (Einheitserde Werkverband e.V., Germany) and stratified for three days at 4°C in the dark. Seven days after germination, the seedlings were transferred to single pots (7 cm×6 cm). Plants were cultivated at 22±1°C, 60±5% relative humidity and a photoperiod of 10 hours with 200 µmol m−2 s−1 light intensity. The photoperiod was changed from 10 hours to 12 hours with the same light intensity when the experiment with the aphids began. Plants were watered three times per week, with quantities as required, before applying the different water conditions.

Rearing the Aphids

B. brassicae and M. persicae were reared on Brassica rapa ssp. chinensis (pak-choi) plants in the insectarium at the Urban Plant Ecophysiology Division of Humboldt University Berlin, Germany, at 22±1°C temperature and a 14-hour photoperiod. To maintain the aphid population, new plants were added roughly every two weeks, and old plants were removed after the aphids had settled on the new plants.

Water Stress Treatment and Experimental Procedures

25-day-old A. thaliana Col-0 plants were assigned to three different irrigation regimes (treatment): 160 ml water/week (water-logged), 80 ml/week (well-watered) and 40 ml/week (drought stress), each delivered three times per week. Plants were placed on 4 cm deep plastic trays (30 cm wide×60 cm long) to retain any excess water and twenty plants were arranged per tray. After one week of different water treatments, when the plants were 32 days old, ten plants per treatment were harvested in pairs by cutting at the base, then flash frozen in liquid nitrogen and stored at −80°C for chemical analysis. The bioassay with aphids was started in parallel and 15 adult B. brassicae or M. persicae were released on each A. thaliana plant. Each plant was caged individually in a transparent plastic cylinder (6 cm diameter×10 cm height) covered with mesh gauze. The experiment with aphid species was performed with 13 replicates. Control plants of each water treatment were also caged. After one week, ten control plants and 10 aphid-treated plants were harvested in pairs for biochemical analysis, resulting in five replicates per treatment. The remaining three plants per treatment were harvested separately for molecular biological studies. Control plants were cut directly at the base and immediately flash-frozen in liquid nitrogen, whereas aphids were brushed off and counted in the insect treatments. Plant material was stored at −80°C before freeze drying samples for biochemical analysis (ALPHA 2- 4 LD plus, Martin Christ Gefriertrocknungsanlagen GmbH, Germany). Samples were ground using a mixer mill MM 301 (Retsch GmbH & Co. KG, Germany).

An additional experiment, in parallel, with similar water treatments and controlled conditions, as described previously, was conducted to collect phloem sap after one and two weeks of water stress. Collected phloem sap samples were analyzed for sugars, amino acids, and glucosinolates. Furthermore, ten plants per water treatment were cultivated to determine the water content of the soil and the plants using the gravimetric method. The fresh weights of soil and plants were recorded, along with the dry weights after 48 h drying at 105°C in an oven. The water content of soil, expressed by weight as the ratio of the water mass present to the dry weight of the soil sample, determined for well-watered, drought, and water-logged treatments was 77.21%, 59.12%, and 81.93%, respectively. Also the dry weight of plants was recorded after two weeks to estimate plant growth under different water conditions.

Collection of Phloem Sap

Phloem exudates were collected from A. thaliana leaves using the EDTA method according to King and Zeevaart [38]. Upper, completely expanded leaves (4th to 5th) were taken from plants cultivated under different water conditions. Phloem exudates were collected twice, after one or two weeks of water stress. Briefly, each leaf was cut at the base, weighed and incubated in 2 ml vials containing 200 µl 15 mM EDTA-disodium salt solution (pH 7.5) for four hours in a humid climate chamber at 20°C. The phloem exudates were then stored at −80°C.

Analysis of Phloem Amino Acid

Phloem exudates from ten leaves were combined for each replicate and mixed gently. Four replicates per treatment were analyzed for amino acid concentrations. The phloem samples were run on an ARACUS amino acid analyzer (membraPure GmbH, Germany) and the concentration of amino acid was calculated using amino acid standards according to the manufacturers’ protocol (Sigma Aldrich, USA). Tryptophan was detected in the exudates, but levels were too low for reproducible quantitative analysis.

Analysis of Phloem Sugar

For determination of phloem sugar, ten phloem sap samples were pooled for each replicate, lyophilized and filled up to 200 µl. Four replicates were processed per treatment. Sucrose, D-glucose, and D-fructose were determined in the phloem sap samples using a UV spectrometric method with an Enzymatic BioAnalysis kit (Boehringer Mannheim/R-Biopharm) and following the manufacturer’s instructions.

Glucosinolate Analysis

Glucosinolates were extracted as described in detail by Mewis et al. [39]. Briefly, a freeze-dried sample (20 mg) was extracted with 70% (v/v) methanol at 80°C in a heat block for 5 min and centrifuged for 5 min at 4000 g. The pellet was re-extracted twice and supernatants were combined. To quantify the glucosinolate content, 60 µl of 1 mM p-hydroxybenzyl GS (isolated from Sinapis alba seeds) was initially added as an internal standard. To convert the methanol-extracted glucosinolate to the desulfo-form, the extract was loaded onto 9 cm Poly-Prep® Chromatography columns (Biorad, CA, USA) filled with 500 µl of a 10% (w/v) suspension of DEAE Sephadex A-25 in 2 M acetic acid and preconditioned with 6 M imidazole-formate solution in 30% formic acid. The column was washed with 0.02 M sodium acetate buffer (pH 4.0) before adding 75 µl of aryl sulfatase solution (Sigma-Aldrich; H-1 from Helix pomatia). The column was incubated overnight. Desulfated glucosinolates were eluted with 1 ml of ultrapure water (Milli-Q quality). Extracts were run on a Dionex Summit P680A HPLC system equipped with an ASI-100 auto sampler and a PDA-100 photodiode array detector. Glucosinolates were separated on a narrow bore 5 µm column (RP18 Acclaim ™ 120, 250×2.1 mm, Dionex) at a flow rate of 0.4 ml min−1 and a column temperature of 25°C. A 43 min gradient program was used with the eluents ultrapure water and 40% (v/v) acetonitrile (HPLC grade). Glucosinolate peaks were monitored at 229 nm. Glucosinolates were identified using standards, retention time, and UV spectra (see Table S1). The glucosinolate amount was calculated from the HPLC peak areas using response factors from desulfo glucosinolate at 229 nm.

Peak identity of the 13 desulfo GS detected in extracts was confirmed with LC-MS analysis on a HPLC-1100 series chromatograph (Agilent Technologies, Böblingen, Germany) coupled to an Esquire 6000 ESI-ion trap mass spectrometer (Bruker Daltonics, Bremen, Germany) operated in negative ion mode in the range m/z 50–700 as described in Mewis et al. (39). Following mass spectra of protonated molecular ions, [M+H]+ were obtained at certain retention times: 1) 3-methylsulfinylpropyl GS, m/z 344 and 183, 7 min, 2) 4-methylsulfinylbutyl GS, m/z 358 and m/z 197, 10 min, 3) 5-methylsulfinylpentyl GS, m/z 372 and 211, 11 min, 4) 6-methysulfinylhexyl GS, m/z 386 and 225, 13 min, 5) 4-hydroxy-indol-3-ylmethyl GS, m/z 385 and 224, 14 min, 6) 7-methylsulfinylheptyl GS m/z 400 and 239, 18 min, 7) 4-methylthiobutyl GS, m/z 342 and 181, 19 min, 8) 8-methylsulfinyloctyl GS, m/z 414 and 253, 20 min 9) indol-3-ylmethyl GS, m/z 369 and 207, 21 min, 10) 4-methoxy-indol-3-ylmethyl GS, m/z 399 and 237, 23 min, 11) 1-methoxy-indol-3-ylmethyl GS, m/z 399 and 237, 27 min, 12) 7-methylthioheptyl GS, m/z 344 and 223, 30 min, and 13) 8-methylthiooctyl GS m/z 398 and 237, 32 min.

For the analysis of glucosinolate in phloem sap, twenty phloem sap samples were pooled for each replicate. Analysis was carried out in five replications. Phloem sap samples were loaded directly on to the DEAE Sephadex A-25 mini columns for desulfation of glucosinolates as described above. Desulfated glucosinolates were analyzed by HPLC.

Camalexin Quantification

Homogenized plant material, 20 mg, was extracted twice with methanol/ultrapure water (4/1, v/v) for 30 min at 65°C. Samples were analyzed by reverse phase HPLC using a RP-18 LiChroCART column (4×250 mm, 5 µm, Merck) and at a flow of 1 ml·min−1. The gradient program started with methanol/ultrapure water at 1∶1 for 2 min, followed by a 10 min linear gradient to 100% methanol and holding for 2 min at 100%, before returning to the initial run conditions. The peak at 10.2 min was identified as camalexin by comparison with verified standards with respect to retention time and UV spectrum (photodiode array detector PDA-100, Dionex) and quantified using a Shimadzu F-10AXL fluorescence detector (318 nm excitation, 370 nm emission).

Chemical Analysis of Flavonoids

Flavonoids in A. thaliana samples were analyzed as aglycones. A 20 mg freeze-dried, ground sample was dissolved in 750 µl 80% (v/v) methanol (pH = 4) and sonicated at a temperature of 4°C for 20 min. After 5 min centrifugation at 4000 g, the extract was collected and the pellet re-extracted twice. Combined supernatants were concentrated in a centrifugation evaporator (Speed Vac, SC 110) to 250 µl and then an equal volume of 2M HCl was added. The extract was heated at 70°C for two hours in a thermo mixer at 400 rpm (Eppendorf, Germany) to hydrolyze the flavonoid glucosides to aglycones. The extract was cooled down and 100% methanol (HPLC grade) was added to the sample to a final volume of 1 ml. The extract was filtered through Spin-X centrifuge tube filters (Sigma-Aldrich, Inc.) and the flow-through was transferred to HPLC vials.

Samples were run on a Dionex Summit P680A HPLC system, equipped with an ASI-100 auto sampler and a PDA-100 photodiode array detector. The aglycones of flavonoids were separated on a narrow bore 3 µm column (AcclaimPA C16, 2.1×150 mm, Dionex) at a flow rate of 0.4 ml min−1 and a column temperature of 25°C. A 32 min gradient program was used with the eluents ultrapure water (pH 3, with acetic acid) and 100% acetonitrile (HPLC grade). Flavonoid aglycone peaks were monitored at 370 nm. Flavonoids were identified using standards, retention time, and UV spectra. The quantity of flavonoids was calculated from HPLC peak areas using the calibration curve from aglycone standards at 370 nm.

RNA Isolation, cDNA Synthesis and Semi-quantitative RT-PCR Analysis

Total RNA was isolated from 100 mg of frozen A. thaliana rosettes with the RNeasy Plant Mini Kit (Qiagen GmbH, Germany) following the standard protocol. RNA was quantified by UV spectroscopy (Nanodrop ND1000, Technology Inc., USA) and its integrity was visually assessed on ethidium bromide stained agarose gels. RNA was first converted to cDNA using the Moloney murine leukemia virus reverse transcriptase according to the Promega protocol (Promega Corporation, USA).

Gene-specific primers were designed for PCR (Table S2) and the PCR reactions were optimized to ensure that the reaction was in the linear range for quantification by diluting the template. The PCR reaction was performed in a total volume of 20 µl, containing 5XGoTaq PCR buffer 4 µl (Promega), 0.2 mM dNTPs, 2.1 mM MgCl2, 0.5 µmol of the forward and reverse primer, 0.1 µl Taq DNA polymerase (Promega), and 1.5 µl of diluted cDNA. Then the PCR program was run on a thermo cycler (Applied Biosystems): 2 min at 96°C for initial denaturation, 30 cycles of 0.5 min at 94°C for denaturation, 0.5 min at 54°C for annealing, and 0.5 min at 72°C for extension, followed by a final extension for 7 min at 72°C. Actin8 (AC8) was used as a reference gene. The forward (f) primer for AC8 was 5′-ATGAAGATTAAGGTCGTGGCAC and the reverse (r) primer 5′-GTTTTTATCCGAGTTTGAAGAGGC. To quantify the expression profiles, PCR products were run in gel electrophoresis (1.5% agarose gel containing SYBR® Safe DNA gel stain) and band intensities were visualized on a Biostep trans-illuminator. Phoretix Total Lab Quant software was used to measure band intensities and calculate the amount of PCR product using the low mass ladder (Invitrogen) for calibration. Expression of genes was normalized by dividing the band intensities with the corresponding Actin8 intensities. The average gene expression levels in drought and water-logged treatments were compared to well-watered treatments, and gene expression levels in insect treatments were related to their controls in each water treatment to estimate the fold changes in transcript levels. The PCR analysis was carry out in 2 technical and 3 biological replications for each gene. We have also included drought stress marker genes such as P5CS1 (DELTA1-PYRROLINE-5-CARBOXYLATE SYNTHASE 1), PP2CA (PROTEIN PHOSPHATASE 2CA), and CDSP32 (CHLOROPLASTIC DROUGHT-INDUCED STRESS PROTEIN OF 32 KD), which have been shown to increase in microarray analysis (40) and also increased in response to drought stress in our analysis.

All primers successfully amplified a band of the correct size and PCR products were cloned into the T-overhang vector pCR2.1 with the TOPO TA cloning kit (Invitrogen). PCR products were fully sequenced to confirm their accuracy.

Results

Effects of the Water Status of Plants on Aphid Performance

Populations of the specialist B. brassicae increased between six and seven-fold within one week on A. thaliana Col-0 plants cultivated under well-watered, drought, and water-logged conditions, with no significant difference in population size between treatments (Fig. 2). In contrast, aphid numbers of the generalist M. persicae on A. thaliana plants differed between water treatments, with significantly higher numbers found on drought-stressed plants (Fig. 2).

[Figure omitted. See PDF.]

Figure 2. Drought stress favors the population growth of Myzus persicae on A. thaliana.

Average population size of B. brassicae and M. persicae after feeding for 7 days on A. thaliana plants cultivated under different water conditions (n = 10, error bars: SE, different letters indicate significant differences in aphid species numbers between treatments, Tukey’s HSD test P<0.05).

https://doi.org/10.1371/journal.pone.0048661.g002

Effects of Water Availability on Plant Health and Metabolites: Plant Growth

Plant’s access to water can influence plant growth and the metabolite composition of phloem sap, which are important for the plant’s tolerance against aphids. We therefore investigated the relevant metabolites and plant parameters that could influence aphid infestation levels.

Plant Growth

The water content of A. thaliana plants depended on the water status of the soil. Water content decreased under both drought and water-logged conditions, although significantly lower water content levels were determined for water-logged plants compared to well-watered plants (Fig. 3A). Also plant growth of A. thaliana was affected by water stress. Biomass production measured as dry weight was significantly lower in plants grown under drought and water-logged conditions than in well-watered plants (Fig. 3B).

[Figure omitted. See PDF.]

Figure 3. Water stress reduced plant water content and dry weight.

Mean water content (A) and final dry weight (B) of A. thaliana plants after cultivation under different water conditions (n = 10, error bars: SD, different letters indicate significant differences in values obtained between treatments, Tukey’s HSD test P<0.05).

https://doi.org/10.1371/journal.pone.0048661.g003

Free Amino Acids Levels in Phloem Sap

Eighteen different amino acids were detected in the phloem exudates of A. thaliana plants, showing that non-essential amino acids such as aspartic acid, serine, asparagines, glutamic acid, alanine, and proline dominate (Table 1). The concentration of amino acids in the phloem sap of plants that had experienced different water treatments varied significantly after one and two weeks of stress application. Total amino acid levels were higher in drought-stressed plants and lower in water-logged plants compared to well-watered plants; however, only the differences between amino acid concentrations in phloem sap of drought-stressed and water-logged plants were significant. Of the non-essential amino acids, particularly proline, glutamic acid, and aspartic acid were significantly higher in the phloem sap of drought-treated A. thaliana plants compared to water-logged plants (Table 1). Among the essential amino acids, concentrations of arginine, isoleucine, leucine, and lysine, were influenced the most by water availability, with levels increasing under drought stress and partially decreasing under water-logged conditions.

[Figure omitted. See PDF.]

Table 1. Concentrations of amino acid in the phloem sap of A. thaliana plants cultivated under different water stress treatments.

https://doi.org/10.1371/journal.pone.0048661.t001

Sugar Levels in Phloem Sap

After one week, the sugars levels (glucose, fructose, and sucrose) in the phloem sap of A. thaliana plants were not influenced by the different water stress treatments. However, after continued cultivation under water stress, significantly higher amounts of sucrose were detected in the phloem sap of drought-stressed plants compared to well-watered and water-logged plants, which both had little or no detectable amounts of this carbohydrate (Fig. 4).

[Figure omitted. See PDF.]

Figure 4. Sucrose levels in phloem sap increased under drought stress.

Sugar levels in the phloem sap of A. thaliana plants cultivated under different water conditions (n = 3, different letters indicate significant differences between treatments, * indicates significant differences among treatments of sucrose only, Tukey’s HSD test P<0.05).

https://doi.org/10.1371/journal.pone.0048661.g004

Glucosinolates in Phloem Sap

Only two glucosinolates (GS) were detected in the phloem sap of mature A. thaliana plants, namely 1-methoxyindol-3-ylmethyl and 4-methoxyindol-3-ylmethyl glucosinolate. These levels did not change significantly within one week of water stress (ANOVA: F2,11 = 1.3, P = 0.337). But as the period of stress continued, 4-methoxyindol-3-ylmethyl glucosinolate concentrations decreased significantly in drought-stressed plants (Fig. 5, Tukey’s HSD test P<0.05).

[Figure omitted. See PDF.]

Figure 5. Drought stress reduced indolyl glucosinolate concentrations in the phloem sap of plants.

Glucosinolate (GS) levels in the phloem sap of A. thaliana plants cultivated for one week under different water conditions (n = 3, different letters indicate significant differences in total glucosinolate levels between water treatments, and different Greek letters indicate significant differences in single compounds between the water treatments, Tukey’s HSD test P<0.05).

https://doi.org/10.1371/journal.pone.0048661.g005

Correlation Results of Aphid Performance to Compounds in Phloem Sap

Simple correlations were used to explore the likelihood that constitutive metabolite levels in the phloem sap of plants under different water conditions might be a basis of plant suitability for aphid species. The number of aphids and increasing indolyl glucosinolate concentrations in the phloem sap were inversely related, whereby the correlation coefficient was stronger for M. persicae with R = −0.802 than for B. brassicae with R = −0.617. In contrast aphid performance was positive related to primary metabolite content in the phloem sap. M. persicae population increased with a somewhat stronger positive correlation to carbohydrate and amino acid concentrations in the phloem sap, with R = 0.718 and R = 0.754 compared to B. brassicae with R = 0.613 and R = 0.542.

Effects of Water Availability on Plant Secondary Metabolites and the Response to Aphid Feeding

Besides analyzing water stress-induced changes in compounds found in phloem sap, which directly influences phloem feeding insects, secondary metabolite accumulation in water-stressed plants was investigated to elucidate changes within the general plant defense response.

Glucosinolate Response

Thirteen different glucosinolates were identified in mature A. thaliana Col-0 rosettes (Table S1). Levels of aliphatic glucosinolates, but not indolyl glucosinolates, changed significantly upon water stress (ANOVA: F2,14 = 5.55, P = 0.021). Higher aliphatic glucosinolate amounts were determined in drought-stressed plants compared to well-watered plants, with intermediate levels in water-logged plants. Significant increases in 4-methylsulfinylbutyl and 5-methylsulfinylpentyl glucosinolate concentrations were responsible for an increase in total aliphatic glucosinolate levels in drought stressed plants (Tukey’s HSD test P<0.05). Furthermore, 4-methoxyindol-3-ylmethyl glucosinolate significantly decreased under drought stress conditions Tukey’s HSD test P<0.05.

Feeding of aphids changed the glucosinolate profile of A. thaliana after one week (Fig. 6). The induced plant changes were aphid species-specific, and especially feeding of M. persicae lead to significantly decreased aliphatic glucosinolates in all water treatments (Fig. 6A). Contrary to the generalist aphid, feeding by the specialist B. brassicae reduced aliphatic glucosinolates significantly only in water-logged plants (Fig. 6A). In response to aphid feeding especially the portion of methylsulfnylalkyl glucosinolates decreased. Indolyl glucosinolates in plants decreased only under well-watered and water-logged conditions in response to feeding of M. persicae (Fig. 6B). Furthermore, under drought conditions feeding by the specialist B. brassicae significantly induced levels of indolyl glucosinolates, especially of 4-methoxyindol-3-ylmethyl glucosinolate (Fig. 6B).

[Figure omitted. See PDF.]

Figure 6. Especially M. persicae reduced aliphatic glucosinolate levels under all water conditions.

Glucosinolate (GS) levels in A. thaliana plants cultivated under different water conditions and after one week feeding of M. persicae (Mp) and B. brassicae (Bb) (n = 5, different letters indicate significant differences in total glucosinolate levels of control and aphid-treated plants within a water treatment, different Greek letters indicate significant differences in single compounds within each water treatments for the main compounds, Tukey’s HSD test P<0.05; minor compounds not significantly different).

https://doi.org/10.1371/journal.pone.0048661.g006

Camalexin Response

Water stress altered levels of camalexin in A. thaliana plants (Fig. 7). Camalexin content was determined as being twice as high in water-logged plants compared to well-watered plants, approximately four-fold less in plants under drought stressed conditions. However, camalexin levels were only significant different between water-logged and drought-stressed plants (Fig. 7). A similar pattern was obtained under infestation with M. persicae, albeit camalexin levels in general were slightly higher. Interestingly, camalexin was induced up to 90-fold in A. thaliana after feeding by B brassicae, regardless of the plant’s water availability.

[Figure omitted. See PDF.]

Figure 7. The aphid species elicited different camelexin accumulation in plants under different water treatments.

Camalexin levels in control and aphid infested A. thaliana plants after two weeks cultivation under different water conditions (n = 5, different letters indicate significant differences between water treatments and bars with an asterisk indicate significant differences in B. brassicae infestation compared to M. persicae and controls, Tukey’s HSD test, P<0.05; WC: well-watered control, DC: drought control, LC: water-logged control, WM: well-watered + M. persicae, DM: drought + M. persicae, LM: water-logged + M. persicae, WB: well-watered + B. brassicae, DB: drought + B. brassicae, LB: water-logged + B. brassicae).

https://doi.org/10.1371/journal.pone.0048661.g007

Flavonoid Response

Mature A. thaliana plants contained isorhamnetin, kaempferol, and quercetin as flavonoids (Fig. 8). Level of flavonoids were not significantly different in the plants after one week of treatment (one way ANOVA, F2,14 = 1.7, P = 0.222). However, after two weeks of water stress, significantly higher levels of flavonoids were detected in water-logged plants, followed by drought-stressed plants, when compared to well-watered plants (Tukey’s HSD test P<0.05). Aphid feeding did not elicit overall changes in the flavonoid content of plants under well-watered conditions. In contrast, plants under drought, and particularly water-logged stress conditions had overall higher levels of flavonoids, which after feeding of both aphid species decreased substantially in (Fig. 8.) Among the above-mentioned flavonoids, only quercetin concentrations decreased significantly following feeding of either aphid under all the water treatments.

[Figure omitted. See PDF.]

Figure 8. Flavonoid levels are influenced by aphids and water stress.

Flavonoid content of A. thaliana after two weeks under different water conditions and one week of infestation with M. persicae (MP) and B. brassicae (BB) compared to non-insect treated plants (n = 5, different lower case letters indicate significant differences in total flavonoid content within the water treatments; different Greek letters indicate significant differences in single compounds within the water treatments and listed on the right; Bars with an asterisk indicate differences within the flavonoid content of water-logged plants compared to respective treatments on well-watered and drought-stress plants for single and total levels; Tukey’s HSD Test P<0.05).

https://doi.org/10.1371/journal.pone.0048661.g008

Water Stress and Aphid-mediated Changes in Glucosinolate Biosynthetic and Regulatory Gene Transcript Levels

Since glucosinolate levels of A. thaliana were influenced by water treatment and aphid feeding, we analyzed the expression of selected genes associated with glucosinolate biosynthesis and signaling pathways. Under drought stress, MAM1 and two transcription factors involved in aliphatic glucosinolate biosynthesis - MYB29 and MYB76 - were up-regulated two-fold, while genes related to indolyl GS remained unchanged, except for CYP79B2, also involved in camalexin biosynthesis, which decreased (Table 2). Also the transcript level of PAD3 responsible for the final step in camalexin biosynthesis, decreased about two-fold. Under water-logged conditions, only MYB76 transcription increased. From genes associated with signaling pathways, only the expression level of PR1, which is associated with the salicylic acid pathway, decreased under drought stress.

[Figure omitted. See PDF.]

Table 2. Transcript levels of glucosinolate biosynthetic, regulatory and camalexin biosynthesis genes in A. thaliana after two weeks’ water stress and one week of aphid feeding (plants39 days old, n = 3 biological replicates).

https://doi.org/10.1371/journal.pone.0048661.t002

Although aliphatic glucosinolates decreased in A. thaliana only after feeding of M. persicae in all water treatments, transcript levels of biosynthetic genes MAM1, CYP79F1, CYP83A1, and UGT74C1 decreased after infestation with either aphid species (Table 2). Feeding of M. persicae did not alter the expression of genes involved in earlier indolyl glucosinolate biosynthesis steps. Furthermore, expression of genes encoding aliphatic glucosinolate biosynthesis transcription factors, i.e. MYB28, MYB29, and MYB76, was mostly down-regulated by M. persicae and B. brassicae in both water stress treatments. Notably, expression of genes encoding indolyl and aliphatic glucosinolate biosynthesis transcription factors responded differentially to aphid feeding in A. thaliana, with increased MYB51 and decreased MYB76. Transcript levels of MYB122 were up-regulated as a result of feeding by B. brassicae in all water treatments.

Interestingly, infestation with B. brassicae increased transcripts of CYP79B2 and CYP79B3 in drought and water-logged plants, which corresponds to the accumulation of indolyl glucosinolates and camalexin in such challenged plants. Feeding of both, M. persicae and B. brassicae, up-regulated the expression of PAD3 in A. thaliana, but induction was stronger in response to B. brassicae. One gene encoding for myrosinase1 in A. thaliana, TGG1, was included in our study. Water stress did not change its expression levels, but plant infestation with M. persicae or B. brassicae down-regulated TGG1 in all treatments.

From the signaling pathway genes, aphid feeding did not change expression levels of CaEF, a gene that responds to ethylene (Table 2). B. brassicae up-regulated PR1 only under drought stress, while the salicylic acid responsive gene remained unchanged upon aphid feeding in both well-watered and water-logged conditions. BGL1, commonly induced by JA, was down-regulated after feeding by either M. persicae or B. brassicae in all water treatments.

[Figure omitted. See PDF.]

Figure 9. Overall effects of aphid species and water stress on secondary metabolite accumulation in A. thaliana.

Scheme depicting the effect of aphid infestation (biotic stress) or water stress (abiotic stress) on inducing (+), reduction (−), or not changing (+/−) plant metabolites.

https://doi.org/10.1371/journal.pone.0048661.g009

Discussion

Water Mediated Effects on Metabolite Composition in Phloem Sap Promoted Only Growth of Myzus Persicae

Aphids feed exclusively on phloem sap and their association with plants is believed to depend on the quality (metabolite composition) of the sap. Besides the negative effects of water stress on plant growth, the present study revealed distinct changes within the chemical metabolite composition in the phloem sap. Drought stress enriched the nutritional quality of phloem sap in A. thaliana Col-0 plants by increasing the amounts of free amino acids and sugars (Table 1 and Fig. 4). Several plant families, including Brassicaceae, translocate sucrose or raffinose together with one polyol in the phloem sap for long-distance allocation of photo-synthetically fixed carbon [40]. As efficient osmolytes, many plant species accumulate polyols under salt or drought stress. This is in agreement with the remarkable accumulation of sucrose in the phloem sap of drought-stressed A. thaliana plants observed in this study. It is also hypothesized that elevated free amino acids act as osmolytes, protecting cells from dehydration. Correspondingly, several essential amino acids (e.g. leucine and isoleucine) and non-essential amino acids (e.g. proline and glutamic acid) increased in plants following drought stress in our study, but declined under water-logged conditions in A. thaliana. Proline, one of the known markers for water stress, increased dramatically in drought stressed plants when compared to their respective controls, whereby proline has been shown to increase in several different plant species and is thought to provide an osmoprotective function [41], [42]. Accumulation of free isoleucine and other branched-chain amino acids in response to drought stress was also reported by Joshi and Jander [43].

Higher sugar content in their diet, especially sucrose, has been shown to increase the performance of aphids [44], and we found a positive correlation of the two aphid species to increasing sugar levels in the phloem sap of plants under different water condition. However, in other studies high sugar content in the sap had no effect on aphids or was considered as barrier to aphid feeding [45]. Also amino acid concentrations, e.g. tyrosine, alanine, leucine, and glutamic acid, in Brassica species were reported to account for 43% of variations in the intrinsic rate of increasing populations of the generalist M. persicae and the specialist B. brassicae [46]. However, the supply of essential amino acids in phloem sap may be comparatively less important for aphid performance, since they are able to obtain these amino acids from their symbiotic bacteria, such as Buchnera aphidicola, in exchange for providing the bacteria with non-essential amino acids, e.g. aspartic/glutamic acid, as well as sucrose [45], [47]. Our observed greater population growth of M. persicae on drought stressed A. thaliana compared to well-watered and water-logged plants is in line with the enriched quality (increased amino acids) of phloem sap in water deficient plants (Fig. 2, Fig. 4, and Table 1). Karley et al. [48] reported that M. persicae and Macrosiphum euphorbiae (Thomas) perform best on young potato plants and with an artificial diet enriched with non-essential amino acids, particularly glutamine and asparagine, which is consistent with our findings. That increasing amino acids in phloem sap is beneficial for M. persicae is also confirmed in the present study, since a positive correlation of aphid population growth to increasing amino acid levels in differently water treated plant was determined. However, the performance of the specialist aphid, B. brassicae was not influenced by the plant’s water availability. Drought stress can make the plant more or less suitable for attacking herbivores, depending on drought intensity and insect species. In fact, there was only a weak relationship of B. brassicae population growth to increasingg amino acid levels indicating other metabolites make drought stress plants less suitable. Water deficiency in Brassica napus also negatively influenced B. brassicae and had no effect on the population growth of Lipaphis erysimi (Kaltenbach) [10].

Secondary metabolites present in the phloem sap can also alter host-plant suitability for aphids and could contribute to the observed different effects of drought stress on aphid performance. The concentration of glucosinolates in the phloem sap of A. thaliana differed between the various water treatments, and we determined lower amounts of 4-methoxyindol-3-ylmethyl glucosinolate under drought stress than in well-watered and water-logged plants (Fig. 5). Using artificial diets it was previously shown that indolyl glucosinolate breakdown products have an anti-feeding effect on M. persicae and reduced aphids’ reproduction [12]. The same authors showed that atr1D mutant plants overproducing indolyl glucosinolates are more resistant to M. persicae, whereas cyp79B2cyp79B3 double mutants lacking indolyl glucosinolates are more suitable to M. persicae growth. A stronger effect of 4-methoxyindol-3-ylmethyl glucosinolate compared to indole-3-ylmethyl glucosinolate was determined for M. persicae in the absence of myrosinase [16]. Our observed lower amount of 4-methoxyindol-3-ylmethyl glucosinolate in the phloem sap of drought-stressed plants, together with higher amino acid and sucrose content, may cause the observed improvement in population growth of M. persicae (Fig. 2). However, glucosinolate levels may not explain the performance of B. brassicae on plants receiving the water stress treatments. In contrast to M. persicae, this specialist can tolerate glucosinolates to some extent and sequesters these compounds to protect itself against enemies [49], [50]. The improvement in phloem sap quality in drought stressed plants had little effect on B. brassicae, indicating a possible alteration in further growth relevant metabolites. Equivalent performance of B. brassicae on broccoli plants cultivated under different water conditions has been reported by Khan et al. [2], [3]. High glucosinolate concentrations in the well-watered and water-logged plants may have influenced the feeding habits of B. brassicae, causing them to ingest more phloem sap, as reported by Gabrys et al. [51]. Consequently, they would ingest amounts of nutrients similar to normal feeding on drought-stressed plants.

Water Stress Alters Plant Defense Metabolites and Response to Aphids

Different water availability may change both the plant’s constitutive levels of secondary metabolites and its induced defense response to herbivory. Plant growth and defenses are restricted by their internal resources and according to the growth-differentiation balance hypothesis, plants must strike a balance between growth processes and production of defensive compounds [52]. Considering these relationships, low stress conditions with unlimited growth will result in no increase in secondary metabolites, whereas moderate water stress may result in elevating secondary metabolites; however, secondary metabolites should decrease under strong and significant water limitation.

With this hypothesis in mind, we found reduced growth in plants under both water stress conditions along with increased levels of aliphatic glucosinolates and flavonoids under drought stress; flavonoids increased most under water-logged conditions (Fig. 3 and Fig. 8). Similarly, elevated levels of flavonoids under drought and water-logged conditions were observed for broccoli plants [3]. In several studies, a drought-mediated glucosinolate elevating effect was determined for Brassica species [53], [54], [55]. However, glucosinolate accumulation was shown to depend on drought intensity, and glucosinolate levels remained unchanged [52] or decreased under drought conditions in some species [3], [56]. Changes in glucosinolate levels in the plant tissue should influence chewing insects, as reported by Mewis et al. [57] and Gutbrodt et al. [56], but may not affect aphids which feed on phloem sap.

It has been proposed that glucosinolates are produced in drought tissue at low turgor and that under these conditions glucosinolate precursors are produced for later use [53]. Indeed, we noticed a trend towards increasing levels of branch chain amino acids, including methionine, in the phloem sap of A. thaliana under drought stress, but levels of 4-methoxy-indol-3ylmethyl glucosinolate decreased. Along with increased aliphatic glucosinolates in drought-stressed plants, transcript levels of MAM1 and two transcription factors involved in aliphatic glucosinolate biosynthesis, MYB29 and MYB76, increased. MYB29 and MYB76 over-expression was reported to noticeably increase levels of 3-methylsulfinylpropyl, 4-methylsulfinylbutyl, and 5-methylsulfinylpentyl glucosinolate [24]. MAM1 catalyzes only the short chain methionine derivatives, resulting in accumulation of short-chained aliphatic glucosinolates [58]. Although MYB76 was slightly up-regulated in water-logged plants, the glucosinolate profile of these plants did not change. MYB76 was shown to integrate different stress responses with aliphatic glucosinolate regulation, while MYB29 is more important for maintaining a basic level of aliphatic glucosinolates, not only in non-stress situations but also in response to environmental cues [24]. Interestingly, both drought and water-logged conditions had no overall-effect on indolyl glucosinolate accumulation and corresponding gene expression.

The two aphids induced partially different plant responses in the present study, and this was influenced by the water available to the plants. Aliphatic glucosinolate decreased after M. persicae feeding on A. thaliana in all water treatments, while feeding of B. brassicae only reduced levels of aliphatic glucosinolate in drought-stressed plants (Fig. 6A). A similar decrease in total aliphatic glucosinolate content in Col-0 plants in response to M. persicae feeding was observed by Kim et al. [12]. Contrary to their results, in our study M. persicae did not influence levels of 4-methoxyindol-3-ylmethyl glucosinolate in Col-0; however this indolyl defense compound accumulated in response to B. brassicae feeding (Fig. 6B). Decreased glucosinolate contents in M. persicae-infested, water-stressed plants were associated with down-regulation of aliphatic glucosinolate biosynthetic (MAM1, CYP79F1, CYP83A1, and UGT74C1) and transcription factor genes (MYB28, MYB29, and MYB76, Table 2). Response to B. brassicae only reduced aliphatic glucosinolates in water-logged plants; however biosynthetic core structure genes as early response were down-regulated in response to feeding under all water conditions. Kuśnierczyk et al. [34] reported that the transcript levels of genes associated with the aliphatic glucosinolate biosynthesis pathway were either unchanged or slightly down-regulated upon infestation by B. brassicae and M. persicae, but in their study, aphids fed on the plant for only three days and not for seven days, as in our study.

Interestingly, the generalist M. persicae and the specialist aphid B. brassicae modulated the transcript profile of indolyl glucosinolate biosynthetic genes differently. While the feeding of M. persicae induced no changes in transcripts involved in earlier steps of glucosinolate biosynthesis in A. thaliana, B. brassicae induced expression of CYP79B2 and CYP79B3 under drought and water-logged conditions, along with increased levels of 4-methoxyindol-3ylmetyl glucosinolate (Table 2). Furthermore, two of the three transcription factors which induce indolyl glucosinolate biosynthesis, MYB51 and MYB122 [22], were only up-regulated in differently water treated plants only in the case of B. brassicae. Unlike our study, Kuśnierczyk et al. [59] found increased expression indolyl glucosinolate biosynthesis genes such as CYP79B2 and CYP79B3 after a short feeding period of both M. persicae and B. brassicae.

TGG1, encoding for the enzyme myrosinase, was down-regulated after infestation by M. persicae and B. brassicae in plants subjected to all water conditions. This observation is consistent with the studies of Kuśnierczyk et al. [59], where both myrosinase encoding genes TGG1 and TGG2 are down-regulated after feeding of B. brassicae and M. persicae. Down-regulation of the myrosinase encoding genes may have been irrelevant for M. persicae and B. brassicae, since they can avoid the myrosinase-glucosinolate defense system of plants. Indeed, knocking out TGG1 and TGG2 did not influence the performance of M. persicae [60]. Rather, this change might increase the plants’ susceptibility to chewing insects and pathogens, since glucosinolate breakdown products have been shown to suppress microbial growth [15].

Feeding of aphids decreased flavonoid accumulation under drought and water-logged stress (Fig. 8). Flavonoids benefit the plant since they are scavengers of reactive oxygen species (ROSs) under stress conditions and can influence insect herbivory [61], [62], but maybe not aphid growth. Lattanzio et al. [63] and Kim et al. [64] suggested that elicitors released by aphids during feeding manipulate flavonoid resistance mechanisms.

Our study also provides evidence that camalexin accumulation depends on the water available to the plant. Transcript levels of the camalexin biosynthetic genes CYP79B2 and PAD3 were down-regulated by drought stress, resulting in less camalexin in plants under drought conditions than in plants under excessive water. The feeding of B. brassicae, but not M. persicae, increased expression of CYP79B2, CYP79B3, and PAD3 resulting in a 50 to 100-fold increase in the concentration of camalexin, known as a phytoalexin active against pathogens, content in plants under all water conditions (Fig. 7). In comparison, the camalexin response elicited by the generalist M. persicae was marginal and statistically not significantly different from the corresponding water treatment controls. Consistent with our finding, a previous study by Kuśnierczyk et al. [34] found increased camalexin in B. brassicae-infested plants. However, their study only included one aphid species and no water stress.

In theory, the improved phloem sap quality in drought-stressed plants should also favor the specialist aphid, and not just the generalist aphid. Since the highest camalexin levels were found in drought-stressed, B. brassicae-infested A. thaliana plants one could speculate whether an alternative chemical defense mechanism exists against this specialist, reasonably glucosinolate-tolerant aphid. The aphid fitness experiment conducted by Kuśnierczyk et al. [34] found that B. brassicae was significantly more fecund on camalexin-lacking pad3-1 mutants than on wild-type plants. This supports our hypothesis that high camalexin levels in drought-treated plants reduce their suitability as hosts. Further studies are required, however, to prove that camalexin is present in phloem sap. Interestingly, Pegadaraju et al. [65] observed no differences in population growth of M. persicae on pad3-1 or wild-type plants. They concluded that camalexin has no effect on this generalist aphid. However, our study rather suggests that the similar performance of M. persicae on the two genotypes could be attributed to the very low induction of camalexin in Col-0 by this aphid species.

Genes associated with signaling pathways were included in this study to elucidate cross talk in plant responses to different stressors. In a drought assay using A. thaliana mutants modified in abscisic acid and jasmonate, Harb et al. [66] found that the associated network of these plant hormones is crucial for the plants’ response to drought. In our study, BGL1 (responsive to jasmonic acid) and CaEF (ethylene signaling) expression did not change [67], but lower levels of PR1 transcript indicate down-regulation of salicylic acid signaling (Table 1). The plants’ response to aphids is partially associated with the salicylic acid pathway, and PR1 is induced by aphids, as reported by Moran and Thomson [17] and Mewis et al. [18]. In this study, PR1 was only up-regulated by B. brassicae in drought-stressed plants, along with accumulation of 4-methoxyindol-3-ylmethyl glucosinolate. Exogenous application of salicylic acid has been shown to specifically promote the accumulation of 4-methoxyindol-3-ylmethyl glucosinolate in A. thaliana [67], [68]. With this in mind, higher salicylic acid levels in response to B. brassicae infestation could lead to the accumulation of 4-methoxyindol-3-ylmethyl glucosinolate in the plant. Drought treatment alone down-regulated the transcripts of PR1 and significantly lowered 4-methoxyindol-3-ylmethyl glucosinolate contents in drought-stressed plants. This finding suggests that the salicylic acid signaling pathway is a key element in regulating 4-methoxyindol-3-ylmethyl glucosinolate production from the indol-3-ylmethyl glucosinolate precursor.

In addition to the importance of salicylic acid signaling, accumulation of glucosinolates is often linked to the jasmonic acid pathway. Exogenous application of methyl jasmonate and jasmonic acid was reported to facilitate the induction of glucosinolate biosynthesis genes, with corresponding increases of indolyl glucosinolates in A. thaliana [57], [68]. In our study, aphid feeding decreased transcripts of BGL1 in plants under all conditions, while levels of CaEF remained unchanged. Corresponding to the down-regulation of jasmonic acid signaling, aliphatic and indolyl glucosinolates decreased after feeding of M. persicae, irrespective of the water treatment. The observed B. brassicae induced accumulation of indolyl glucosinolate under drought stress seems to be jasmonate-independent (Fig. 6B) and might be also one explanation for the similar performance on drought and well-water plants, which might be contrary to the initial lower constitutive levels of indolyl glucosinolates in the phloem sap.

Conclusions

The drought stress-mediated enrichment of phloem sap with free amino acids, along with the decreased concentration of 4-methoxyindol-3-ylmethyl glucosinolate promoted population growth of the generalist aphid M. persicae on A. thaliana. The specialist B. brassicae did not particularly profit from the increased phloem sap quality, since it performed equally well on plants cultivated under different water conditions. Plants responded to infestation with B. brassicae, but not with M. persicae, by producing high levels of camalexin in plant tissues under all water conditions. This may lead to decreased host plant suitability for this aphid species. However, further investigations are required to determine possible transport, and, therefore, presence of camalexin in phloem sap.

Water stress and aphids generally influenced the accumulation of secondary metabolites an opposite ways and participated differentially in signaling pathways as summarized in Figure 9. Water-logging and drought stress led to increased aliphatic glucosinolate and flavonoid levels in the plants. Conversely, aphid feeding, especially of M. persicae, reduced the accumulation of glucosinolates and flavonoids in most water treatments, along with down-regulation of the jasmonic acid pathway. How the plants responded with secondary metabolite accumulation was obviously very different after feeding of M. persicae and B. brassicae and depended on the water status of the plants. For example, B. brassicae induced the salicylic acid pathway only under drought stress, meaning the plants accumulated 4-methoxyindol-3-ylmethyl glucosinolate. Different chemical elicitor present in the saliva might mediate the discriminative plants’ response to M. persicae and B. brassicae, a question that remains to be investigated in the future. Our findings suggest a complex network of regulating secondary metabolite biosynthesis in response to multiple environmental factors and a fine-tuned plant response to stressors, excluding a general prediction. Further experimentations should include GUS transformed A. thaliana mutants to visualize camalexin induction by the different aphid species. Also a direct comparison of M. persicae and B. brassicae performance on indolyl GS lacking cyp79B2cyp79B3 mutants and camalexin deficient pad3-1 mutants under different water conditions may be worth further investigation.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

UV-spectra of desulfo GS from Arabidopsis (Col-0).

https://doi.org/10.1371/journal.pone.0048661.s001

(PDF)

Table S2.

Primers list of genes associated with glucosinolate, camalexin biosynthesis or signaling pathways.

https://doi.org/10.1371/journal.pone.0048661.s002

(PDF)

Acknowledgments

We thank Dr. Sonja Löffler, LFE, Brandenburg, Germany and Martina Hitzenbichler, TU München, Germany for their excellent technical assistance in analyzing amino acids and camalexin.

Author Contributions

Conceived and designed the experiments: MK IM CU. Performed the experiments: MK IM EG. Analyzed the data: MK IM. Contributed reagents/materials/analysis tools: IM MS EG CU. Wrote the paper: IM MK MS EG CU.

Citation: Mewis I, Khan MAM, Glawischnig E, Schreiner M, Ulrichs C (2012) Water Stress and Aphid Feeding Differentially Influence Metabolite Composition in Arabidopsis thaliana (L.). PLoS ONE 7(11): e48661. https://doi.org/10.1371/journal.pone.0048661

References

1. Garg BK, Kathju S, Burman U (2001) Influence of water stress on water relations, photosynthetic parameters and nitrogen metabolism of moth bean genotypes. Biologia Plantarum 44: 289–292.

2. Khan MAM, Ulrichs C, Mewis I (2010) Influence of water stress on the glucosinolate profile of Brassica oleracea var. italica and the performance of Brevicoryne brassicae and Myzus persicae.. Entomologia Experimentalis et Applicata 137: 229–236.

3. Khan MAM, Ulrichs C, Mewis I (2011) Water stress alters aphid induced glucosinolate response in Brassica oleracea var. italica differently. Chemoecology 21: 235–242.

4. Chandra A, Bhatt RK, Misra LP (1998) Effect of water stress on biochemical physiological characteristics of oat genotypes. Journal of Agronomy and Crop Science 181: 45–48.

5. Pinheiro C, Chaves RM, Ricardo CP (2001) Alteration of carbon and nitrogen metabolism induced by water deficit in the stems and leaves of Lupinus albus L. Journal of Experimental Botany. 52: 1063–1070.

6. Streeter JG, Lohnes DG, Fioritto RJ (2001) Paterns of pinitrol accumulation in soybean plants and relationships to drought tolerance. Plant Cell and Environment 24: 429–438.

7. Ingram J, Bartels D (1996) The molecular basis of dehydration tolerance in plants. Annual Review of Plant Physiology and Plant Molecular Biology 47: 377–403.

8. Bray EA (1997) Plant response to water deficit. Trends of Plant Science 2: 48–54.

9. Wearing CH (1972) Responses of Myzus persicae and Brevicoryne brassicae to leaf age and water stress in Brussels sprouts grown in pots. Entomologia Experimentalis et Applicata 15: 61–80.

10. King C, Jacob HS, Berlandier F (2006) The influence of water deficiency on the relationship between canola (Brassica napus L.), and two aphid species (Hemiptera: Aphididae), Lipaphis erysimi (Kaltenbach) and Brevicoryne brassicae (L.). Australian Journal of Agricultural Research 57: 439–445.

11. Robbins RJ, Keck AS, Banuelos G, Finley JW (2005) Cultivation conditions and selenium fertilization alter the phenolic profile, glucosinolate, and sulforaphane content of broccoli. Journal of Medicinal Food 8: 204–214.

12. Kim JH, Lee BW, Schroeder FC, Jander G (2008) Identification of indole glucosinolate breakdown products with antifeedant effects on Myzus persicae (green peach aphid). The Plant Journal 54: 1015–1026.

13. Kelly PJ, Bones AM, Rossiter JT (1998) Sub-cellular immunolocalization of the glucosinolate singrin in seedlings of Brassica juncea. Planta 206: 370–377.

14. Koroleva OA, Davies A, Deeken R, Thorpe MR, Tomos AD, et al. (2000) Identification of a new glucosinolate-rich cell type in Arabidopsis flower stalk. Plant Physiology 124: 599–608.

15. Halkier BA, Gershenzon J (2006) Biology and biochemistry of glucosinolates. Annual Review of Plant Biology 57: 303–333.

16. Kim JH, Jander G (2007) Myzus persicae (green peach aphid) feeding on Arabidopsis induces the formation of a deterrent indole glucosinolate. The Plant Journal 49: 1008–1019.

17. Moran PJ, Thompson GA (2001) Molecular responses to aphid feeding in Arabidopsis in relation to plant defense pathways. Plant Physiology 125: 1074–1085.

18. Mewis I, Tokuhisa JG, Schultz JC, Appel HM, Ulrichs C, et al. (2006) Gene expression and glucosinolate accumulation in Arabidopsis thaliana in response to generalist and specialist herbivores of different feeding guilds and the role of defense signaling pathways. Phytochemistry 67: 2450–2462.

19. Sønderby IE, Geu-Flores F, Halkier BA (2010) Biosynthesis of glucosinolates – gene discovery and beyond. Trends in Plant Science 15: 283–290.

20. Geu-Flores F, Møldrup ME, Böttcher C, Olsen CE, Scheel D, et al. (2011) Cytosolic {gamma}-glutamyl peptidases process glutathione conjugates in the biosynthesis of glucosinolates and camalexin in Arabidopsis. The Plant Cell 23: 2456–2469.

21. Hirai MY, Sugiyama K, Sawada Y, Tohge T, Obayashi T, et al. (2007) Omics-based identification of Arabidopsis Myb transcription factors regulating aliphatic glucosinolate biosynthesis. Proceedings of the National Academy of Sciences USA 104: 6478–6483.

22. Gigolashvili T, Berger B, Mock HP, Müller C, Weisshaar B, et al. (2007) The transcription factor HIG1/MYB51 regulates indolic glucosinolate biosynthesis in Arabidopsis thaliana. The Plant Journal 50: 886–901.

23. Gigolashvili T, Yatusevich R, Berger B, Müller C, Flügge UI (2007) The R2R3-MYB transcription factor HAG1/MYB28 is a regulator of methionine-derived glucosinolate biosynthesis in Arabidopsis thaliana. The Plant Journal 51: 247–261.

24. Gigolashvili T, Engqvist M, Yatusevich R, Müller C, Flügge UI (2008) HAG2/MYB76 and HAG3/MYB29 exert a specific and coordinated control on the regulation of aliphatic glucosinolate biosynthesis in Arabidopsis thaliana. New Phytologist 177: 627–642.

25. Celenza JL, Quiel JA, Smolen GA, Merrikh H, Silvestro AR, et al. (2005) The Arabidopsis ATR1 myb transcription factor controls indolic glucosinolate homeostasis. Plant Physiology 137: 253–262.

26. Yan X, Chen S (2007) Regulation of plant glucosinolate metabolism. Planta 226: 1343–1352.

27. Zhao Y, Hull AK, Gupta NR, Goss KA, Alonso J, et al. (2002) Trp-dependent auxin biosynthesis in Arabidopsis: involvement of cytochrome P450s CYP79B2 and CYP79B3. Genes and Development 16: 3100–3112.

28. Glawischnig E, Hansen BG, Olsen CE, Halkier BA (2004) Camalexin is synthesized from indole-3-acetaldoxime, a key branching point between primary and secondary metabolism in Arabidopsis. Proceedings of the National Academy of Sciences USA 101: 8245–8250.

29. Glawischnig E (2007) Camalexin. Phytochemistry 68: 401–406.

30. Rogers EE, Glazebrook J, Ausubel FN (1996) Mode of action of the Arabidopsis thaliana phytoalexin camalexin and its role in Arabidopsis–pathogen interactions. Molecular Plant-Microbe Interaction 9: 748–757.

31. Ferrari S, Galletti R, Denoux C, De Lorenzo G, Ausubel FM, et al. (2007) Resistance to Botrytis cinerea induced in Arabidopsis by elicitors is independent of salicylic acid, ethylene, or jasmonate signaling but requires PHYTOALEXIN DEFICIENT3. Plant Physiology 144: 367–379.

32. Sellam A, Iacomi-Vasilescu B, Hudhomme P, Simoneau P (2007) In vitro antifungal activity of brassinin, camalexin and two isothiocyanates against the crucifer pathogens Alternaria brassicicola and Alternaria brassicae. Plant Pathology 56: 296–301.

33. Schuhegger R, Rauhut T, Glawischnig E (2007) Regulatory variability of camalexin biosynthesis. Plant Physiology 164: 636–644.

34. Kuśnierczyk A, Winge P, Jørstad TS, Troczyńska J, Rossiter JT, et al. (2008) Towards global understanding of plant defense against aphids – timing and dynamics of early Arabidopsis defense responses to cabbage aphid (Brevicoryne brassicae) attack. Plant Cell and Environment 31: 1097–1115.

35. Ahuja I, Rohloff J, Bones AM (2010) Defense mechanisms of Brassicaceae: implications for plant-insect interactions and potential for integrated pest management. A review. Agronomy for Sustainable Development 30: 311–348.

36. Schuhegger R, Nafisi M, Mansourova M, Petersen BL, Olsen CE, et al. (2006) CYP71B15 (PAD3) catalyzes the final step in camalexin biosynthesis. Plant Physiology 141: 1248–1256.

37. Böttcher C, Westphal L, Schmotz C, Prade E, Scheel D, et al. (2009) The multifunctional enzyme CYP71B15 (PHYTOALEXIN DEFICIENT3) converts cysteine-indole-3-acetonitrile to camalexin in the indole-3-acetonitrile metabolic network of Arabidopsis thaliana. The Plant Cell 21: 1830–1845.

38. King RW, Zeevaart JAD (1974) Enhancement of phloem exudation from cut petioles by chelating agents. Plant Physiology 53: 96–103.

39. Mewis I, Schreiner M, Nguyen CN, Krumbein A, Ulrichs Ch, et al. (2012) UV-B irradiation changes specifically the secondary metabolite profile in broccoli sprouts – Induced signalling overlaps with the plant response to biotic stressors. Plant Cell Physiology 53: 1546–1560.

40. Pommerrenig B, Papini-Terzi FS, Sauer N (2007) Differential regulation of sorbitol and sucrose loading into the phloem of Plantago major in response to salt stress. Plant Physiology 144: 1029–1038.

41. Witt S, Galicia L, Lisec J, Cairns J, Tiessen A, et al. (2012) Metabolic and phenotypic responses of greenhouse-grown maize hybrids to experimentally controlled drought stress: Molecular Plant. 5: 401–417.

42. Bowne JB, Erwin TA, Juttner J, Schnurbusch T, Langridge P, et al. (2012) Drought responses of leaf tissues from wheat cultivars of differing drought tolerance at the metabolite level. Molecular Plant 5: 418–429.

43. Joshi V, Jander G (2009) Arabidopsis methionine γ-lyase is regulated according to isoleucine biosynthesis need but plays a subordinate role to threonine deaminase. Plant Physiology 151: 367–378.

44. Simpson J, Abisgold JD, Douglas AE (1995) Response of the pea aphid (Acyrthosiphon pisum) to variation in dietary levels of sugar and amino acids: the significance of amino acid quality. Journal of Insect Physiology 41: 71–75.

45. Douglas AE (2006) Phloem-sap feeding by animals: problems and solutions. Journal of Experimental Botany 57: 747–754.

46. Cole RA (2003) The relative importance of glucosinolates and amino acids to the development of two aphid pests Brevicoryne brassicae and Myzus persicae on wild and cultivated brassica species. Entomologia Experimentalis et Applicata 85: 121–133.

47. Gündüz EA, Douglas AE (2009) Symbiotic bacteria enable insect to use a nutritionally inadequate diet. Proceeding of the Royal Society. Series B 276: 987–991.

48. Karley AJ, Douglas AE, Parker WE (2002) Amino acid composition and nutritional quality of potato leaf phloem sap for aphids. Journal of Experimental Biology 205: 3009–3018.

49. Kazana E, Pope TW, Tibbles L, Bridges M, Pickett JA, et al. (2007) The cabbage aphid: a walking mustard oil bomb. Proceeding of the Royal Society London Series B 274: 2271–2277.

50. Pratt C, Pope TW, Powell G, Rossiter JT (2008) Accumulation of glucosinolates by the cabbage aphid Brevicoryne brassicae as a defense against two Coccinellid species. Journal of Chemical Ecology 34: 323–329.

51. Gabrys B, Tjallingii WF, van Beek TA (1997) Analysis of EPG recorded probing by cabbage aphid on host plant parts with different glucosinolate contents. Journal of Chemical Ecology 23: 1661–1673.

52. Mattyssek R, Schnyder H, Elster EF, Munch JC, Pretsch H, et al. (2002) Growth and parasite defense in plants: the balance between resource sequestration and retention. Plant Biology 4: 133–136.

53. Jensen CR, Mogensen VO, Mortensen G, Fieldsen JK, Thage JH (1996) Seed glucosinolate, oil and protein contents of field-grown rape (Brassica napus) affected by soil drying and evaporative demand. Field Crops Research 47: 93–105.

54. Schreiner M (2005) Vegetable crop management strategies to increase the quantity of phytochemicals. European Journal of Nutrition 44: 85–94.

55. Schreiner M, Beyene B, Krumbein A, Stützel H (2009) Ontogenetic changes of 2-propenyland 3-indolylmethyl glucosinolates in Brassica carinata leaves as affected by water supply. Journal of Agriculture and Food Chemistry 57: 7259–7263.

56. Gutbrodt B, Mody K, Dorn S (2011) Drought changes plant chemistry and causes contrasting response in lepidopteran herbivores. Oikos 120: 1732–1740.

57. Mewis I, Appel HM, Hom A, Raina R, Schultz JC (2005) Major signaling pathways modulate Arabidopsis glucosinolate accumulation and response to both phloem-feeding and chewing insects. Plant Physiology 138: 1149–1162.

58. Textor S, de Kraker JW, Hause B, Gershenzon J, Tokuhisa JG (2007) MAM3 catalyzes the formation of all aliphatic glucosinolate chain lengths in Arabidopsis thaliana.. Plant Physiology 144: 60–71.

59. Kuśnierczyk A, Winge P, Midelfart H, Armbruster WS, Rossiter JT, et al. (2007) Transcriptional responses of Arabidopsis thaliana ecotypes with different glucosinolate profiles after attack by polyphagous Myzus persicae and oligophagous Brevicoryne brassicae. Journal of Experimental Botany 58: 2537–2552.

60. Barth C, Jander G (2006) Arabidopsis myrosinases TGG1 and TGG2 have redundant function in glucosinolate breakdown and insect defense. The Plant Journal 46: 549–562.

61. Mittler R (2002) Oxidative stress, antioxidants and stress tolerance. Trend Plant Science 7: 405–410.

62. Treutter D (2006) Significance of flavonoids in plant resistance: a review. Environmental Chemistry Letters 4: 147–157.

63. Lattanzio V, Arpaia S, Cardinali A, Di Venere D, Linsalata V (2000) Role of endogenous flavonoids in resistance mechanism of Vigna to aphids. Journal of Agriculture and Food Chemistry 48: 5316–5320.

64. Kim KS, Bilgin DD, Diers BW (2011) Identification of flavonoids associated with soybean aphid resistance in soybean. Plant and Animal Genomes XIX Conference. January 15–19, San Diego, CA. p 398.

65. Pegadaraju V, Knepper C, Reese J, Shah J (2005) Premature leaf senescence modulated by the Arabidopsis PHYTOALEXIN DEFICIENT4 gene is associated with defense against the phloem-feeding green peach aphid. Plant Physiology 139: 1927–1934.

66. Harb A, Krishnan A, Ambavaram MMR, Pereira A (2010) Molecular and physiological analysis of drought stress in Arabidopsis reveals early responses leading to acclimation in plant growth. Plant Physiology 154: 1254–1271.

67. Mikkelsen MD, Petersen BL, Glawischnig E, Jensen AB, Andreasson E (2003) Modulation of CYP79 genes and glucosinolate profiles in Arabidopsis by defense signaling pathways. Plant Physiology 131: 298–308.

68. Kliebenstein DJ, Figuth A, Mitchell-Olds T (2002) Genetic architecture of plastic methyl jasmonate responses in Arabidopsis thaliana. Genetics 161: 1685–1696.