Site description and plant material
The research was conducted on a sun exposed area in Limassol, Cyprus (34°42'N, 32°59'E; elevation: 100 m a.s.l.), during 22 rainless days in May 2018. The climate is Mediterranean, with hot and dry summers. Supplementary Fig. S1 presents the climatic data, recorded with an on-site data logger (Kistock KH 250; Kimo). On an average day during the experiment, a maximum temperature of 38.3 °C was achieved in full shade, corresponding with 24% RH. On an average night, the minimum temperature dropped to 19.5 °C and the relative humidity reached 70%.
This study comprised two cultivars of Vitis vinifera, Xynisteri and Chardonnay. Xynisteri is the main white grape cultivar grown in Cyprus, while Chardonnay, one of the most planted white grape cultivars internationally, has been introduced in Cyprus. In 2014, they respectively covered 30.2% and 1.6% of the ca. 6,142 ha viticultural area of Cyprus (Statistical service of Republic of Cyprus 2016). Of each cultivar, 60 self-rooted cuttings were planted in 5-liter polyethylene pots containing soil, originating from the traditional vineyard area in Limassol. The soil properties were previously described by Tzortzakis et al. (2020) and briefly, the soil had a clay-loam texture, organic matter of 2.19%; total CaCO3 of 66.9%; pH of 7.42; electrical conductivity (EC) 0.28 of mS cm− 1. The plants were grown in field conditions and were automatically irrigated at field capacity using a drip irrigation system. Three months after planting, the plants were uniformly distributed over 12 treatment groups. For each treatment, five replicates were used per cultivar. The experimental set-up is shown in Fig. 1. Each group was treated with one of four abiotic stress treatments (7 or 14 days of full or deficit irrigation) to assess the effect of short and prolonged drought stress. Two groups were sampled destructively at 7 and 14 days of treatment (dot). In vitro inoculations were performed on disks of these leaves. In the evening, the remaining intact plants were inoculated with either pathogen or water. For these plants, the irrigation regime was maintained until disease evaluation. Some leaves were sampled at 9 and 16 dot to establish the effect of pathogen attack at 1.5 days post inoculation (dpi) on drought-stressed plants.
Abiotic stress
Plants were either well-watered in the full irrigation control treatment or exposed to drought stress by deficit irrigation. Fully irrigated plants received irrigation at field capacity from an automatic drip system, every 6 hours for 5 min. Deficit irrigation was maintained at 40% of the full irrigation, based on the volumetric water content of the soil (VWC). The deficit-irrigated plants were irrigated manually every two days. To verify and accurately adjust the irrigation, the VWC was measured daily in 8 randomly chosen pots using a portable Time-Domain Reflectometer (TDR) (FieldScout TDR 300 Soil Moisture Probe; Spectrum Technologies) with 4.7 inch rods (Supplementary Fig. S2).
Biotic stress
To examine the combined effect of abiotic and biotic stress on the vine, a pathogen stress was imposed on the intact plants after 7 or 14 days of drought stress. Plasmopara viticola isolate Fcpv1, obtained from Chardonnay in France, was grown for 10 days at 22 °C on detached Chardonnay leaves on water agar (0.65%). Sporangia were collected with distilled water and the suspension was adjusted to 2.5 × 104 sporangia mL− 1. The artificial inoculation was performed in the evening. The abaxial sides of all leaves were sprayed until run-off with 3 mL sporangia suspension of P. viticola. The control plants were sprayed with distilled water. The irrigation regimes were maintained until the disease evaluation. Since P. viticola needs 95–100% relative humidity during the night for an optimal infection and sporulation, each plant was equipped with a container of water and a humid plastic cover in the evening. To prevent extreme temperature development within the cover, the cover was removed in the morning and a light shade was created using a shadow mesh.
Sampling and disease evaluation took place at 1.5 and 7 dpi respectively. Each plant was evaluated according to the following classes: 0, no symptoms; 1, few oil spots with little to no sporulation; 2, moderate symptoms and non-spreading sporulation; 3, clearly diseased with spreading sporulation; 4, severe symptoms with dense sporangiophore carpets.
Field measurements
At 3, 7, 9, 14 and 16 dot, stomatal conductance, chlorophyll fluorescence and chlorophyll content were recorded. The measurements were conducted on the 4th or 5th leaf starting from the apical meristem on randomly chosen plants at mid-morning, 4 h after onset of light. The stomatal conductance to water vapor (gs) was measured on three to five plants, using a transient state diffusion porometer (AP4; Delta-T Devices). The chlorophyll fluorescence (Fv/Fm), an indicator of the maximum quantum efficiency of Photosystem II, was monitored on three or four plants after exposure to darkness for 20 minutes with a dark adaptation pin using a chlorophyll fluorometer (OS30p; Opti-Sciences). The chlorophyll content per unit leaf area was estimated using a non-destructive Soil Plant Analysis Development (SPAD) meter (SPAD 502 Plus; Spectrum Technologies). SPAD measurements were conducted twice on the same leaf of five or six plants.
In vitro assessment of disease susceptibility
The 3rd and 4th leaf, counted from the apex, sampled at 7 and 14 dot before application of biotic stress treatment, were used to investigate the effect of the previous exposure to drought stress on the susceptibility to P. viticola. Leaf disks (11 mm diameter) were treated with 20 µL distilled water or 20 µL P. viticola sporangia suspension containing 2.5 × 104 sporangia mL− 1 and incubated on water agar (0.65%) at 22 °C. At 5 dpi, the number of sporangiophores was counted to assign each disk to one of the following classes: 0, 0 sporangiophores; 1, 1–6 sporangiophores; 2, 7–20 sporangiophores; 3, > 20 sporangiophores; 4, dense sporangiophore carpet. An average of 60 disks were evaluated per treatment.
Quantification of phytohormones
The levels of abscisic acid (ABA), indole-3-acetic acid (IAA), jasmonic acid (JA) and salicylic acid (SA) were determined in leaves sampled at 7, 9, 14 and 16 dot. For each of the five replicates per treatment, two leaves were pooled, immediately frozen in liquid N2 and kept at -80 °C until analysis. The procedure for quantification of phytohormones is described in detail by Haeck et al. (2018). The ground tissue (100 mg) was incubated with 5 mL modified Bieleski extraction solvent (methanol/water/formic acid 75:20:5, v/v/v) for 20–24 h at -80 °C. After this cold extraction, filtration (30 kDa Amicon® Ultra centrifugal filter unit, Merck Millipore, Overijse, Belgium) and evaporation (TurboVap® LV, Biotage, Uppsala, Sweden), the extracts were reconstituted in 0.5 mL methanol/water/formic acid (20:80:0.1, v/v/v). Chromatographic separation was performed on an ultra-high-performance liquid chromatography system (U-HPLC, Thermo Fisher Scientific) equipped with a Nucleodur C18 column (50 × 2 mm; 1.8 µm particle diameter). Mass spectrometric analysis was achieved in targeted single ion monitoring mode on a Q-Exactive™ quadrupole-Orbitrap mass spectrometer (Thermo Fisher Scientific), equipped with a heated electrospray ionization source, at a resolution of 70,000 full width at half maximum. In negative ionization mode, SA, ABA and JA were measured using an elution gradient (300 µL min− 1) of (A) methanol and (B) water, both with 0.01% formic acid. The formic acid concentration of solvent B was adjusted to 0.1% for measurement of IAA in positive ionization mode. The following linear gradient was applied (solvent A): 0–1 min at 20%, 1-2.5 min from 20 to 45%; 2.5-9 min from 45 to 100%; 9–10 min at 100%; 10–14 min at 20%. External and deuterated internal (d4-SA at 200 µg L− 1, d6-ABA and d5-IAA at 1 µg L− 1) standards were used for accurate quantification of the hormone content.
Quantification of photosynthetic pigments
Leaf samples were collected at 7, 9, 14 and 16 dot with five replications per treatment, each consisting of a pool of two leaves. Leaf tissue (100 mg) was incubated in a heat bath at 65 °C for 30 min, with 10 mL dimethyl sulfoxide (DMSO). The absorbance of the extract was measured at 645 nm and 663 nm, using a microplate spectrophotometer (Thermo Scientific, Multiskan GO) and chlorophyll a (Chl a) and chlorophyll b (Chl b) concentrations were calculated as described by Richardson et al. (2002).
Quantification of hydrogen peroxide content, lipid peroxidation and proline
For the quantification of hydrogen peroxide (H2O2), lipid peroxidation, in terms of malondialdehyde (MDA) content, and proline, two leaves were sampled and pooled for each of the five plants per treatment at 7, 9, 14 and 16 dot. Fresh leaves were immediately frozen in liquid N2 and kept at -80 °C until analysis. Before analysis, ground leaf tissue (200 mg) was homogenized with ice-cold 0.1% trichloroacetic acid (TCA). The extract was centrifuged and the supernatant was used for the quantification of H2O2 and MDA (Chrysargyris et al. 2017). For the quantification of H2O2, 0.5 mL of the supernatant was mixed with 0.5 mL of 10 mM potassium phosphate buffer (PPB) (pH 7.0) and 1 mL of 1M potassium iodide (KI). The content of H2O2 was calculated using standards of 5 to 500 µM of H2O2 and a calibration curve was plotted accordingly. The absorbance was measured at 390 nm. For the MDA content, 0.5 mL of the supernatant was incubated with 1.5 mL of 0.5% thiobarbituric acid (TBA) in 20% TCA at 95 °C for 25 min. The reaction was stopped in an ice bath and the absorbance was measured at 532 nm and 600 nm. The MDA content was calculated using the extinction coefficient of 155 mM cm− 1.
Proline content was also determined using this frozen ground tissue. Leaf tissue (200 mg) was homogenized in 2 mL of 3% aqueous sulfosalicylic acid (SSA). Extracts were then centrifuged and 1 mL of the supernatant was incubated with 1 mL of acid ninhydrin and 1 mL of glacial acetic acid, for 1 h at 100 °C. Then, the formed chromogen was extracted with toluene and the absorbance was measured at 520 nm, using toluene as blank. The proline concentration was determined using serial dilutions (0-100 µg mL− 1) of D-proline (Khedr 2003).
Quantification of antioxidant enzymes
The ground leaf samples were also used for the determination of the activity of the antioxidant enzymes. The tissue (200 mg) was homogenized with 3 mL ice cold 50 mM potassium phosphate buffer (pH 7.0), including 1 mM ethylenediaminetetraacetic acid (EDTA), 1% w/v polyvinylpolypyrrolidone (PVPP), 1 mM phenylmethylsulfonyl fluoride (PMSF) and 0.05% polyethylene glycol tert-octylphenyl ether (Triton X-100). The homogenate was centrifuged at 16,000 g for 20 min, at 4 °C. The supernatant was collected and an aliquot was first used to determine the protein content via the Bradford method (1976), with bovine serum albumin (BSA) as the protein standard.
Catalase (CAT, EC 1.11.1.6) activity was determined by following the consumption of H2O2 (extinction coefficient 39.4 mM cm− 1) at 240 nm for 3 min, as assayed by Jiang and Zhang (2002). The reaction mixture contained 100 mM PPB (pH 7.0), plant extract and 200 µL of 75 mM of H2O2. Results were expressed as CAT units per milligram of protein. One unit of enzyme decomposed 1 µmol of H2O2 per min.
Superoxide dismutase (SOD, EC 1.15.1.1) was assayed using the photochemical method. The reaction mixture (1.5 mL) contained 50 mM PPB (pH 7.5), 13 mM methionine, 75 µM nitro blue tetrazolium (NBT), 0.1 mM EDTA, 2 µM riboflavin and an enzyme aliquot. Reaction started after the addition of riboflavin. Tubes containing the reaction were then placed under a light source of two 15-watt fluorescent lamps for 15 min. The reaction was stopped by placing the tubes in the dark. Reaction without the extract developed maximal color (control) and non-irradiated mixture was used as a blank. The absorbance was determined at 560 nm and activity was expressed as SOD units per mg of protein. One unit of SOD activity was defined as the amount of enzyme required to cause 50% inhibition of the NBT photoreduction rate (Chrysargyris et al. 2018).
Peroxidase activity (POD, EC 1.11.1.7) was determined according to the method used by Tarchoune et al. (2012). POD activity was assayed using pyrogallol, following the increase in absorbance at 430 nm, after the oxidation to purpurgallin. The reaction mixture of 2 mL contained 1,665 µL of 100 mM PPB (pH 6.5), 200 µL of 100 mM pyrogallol and 50 µL of extract. The reaction started with the addition of 85 µL of 40 mM H2O2. The increase in absorbance at 430 nm was measured on a kinetic cycle for 3 min. Calculations were performed using the coefficient of extinction of 2.47 mM cm− 1. One POD unit was defined as the amount of enzyme needed to decompose 1 µmol of H2O2 per min.
Statistical analysis
All statistical analysis were conducted using R, version 3.6.1 (R Core Team 2019). Normality and homoscedasticity were verified with the Shapiro-Wilk’s and Levene’s test (p = 0.05). Data meeting the assumptions, were treated with one-way analysis of variance (ANOVA), followed by Tukey’s HSD (honestly significant difference) to compare the means (p = 0.05). Nonparametric data were analyzed using the Kruskal-Wallis test, followed by the Mann-Whitney U-test (p = 0.05). For the analysis of the interactions between the cultivar, drought and pathogen stress, linear regression analysis was performed. A generalized least squares (GLS) model was improved by eliminating interaction terms until the lowest Akaike Information Criterion score (AIC) was reached.