Leaf Eh and pH: A Novel Indicator of Plant Stress. Spatial, Temporal and Genotypic Variability in Rice (Oryza sativa L.)

A wealth of knowledge has been published in the last decade on redox regulations in plants. However, these works remained largely at cellular and organelle levels. Simple indicators of oxidative stress at the plant level are still missing. We developed a method for direct measurement of leaf Eh and pH, which revealed spatial, temporal, and genotypic variations in rice. Eh (redox potential) and Eh@pH7 (redox potential corrected to pH 7) of the last fully expanded leaf decreased after sunrise. Leaf Eh was high in the youngest leaf and in the oldest leaves, and minimum for the last fully expanded leaf. Leaf pH decreased from youngest to oldest leaves. The same gradients in Eh-pH were measured for various varieties, hydric conditions, and cropping seasons. Rice varieties differed in Eh, pH, and/or Eh@pH7. Leaf Eh increases and leaf pH decreases with plant age. These patterns and dynamics in leaf Eh-pH are in accordance with the pattern and dynamics of disease infections. Leaf Eh-pH can bring new insight on redox processes at plant level and is proposed as a novel indicator of plant stress/health. It could be used by agronomists, breeders, and pathologists to accelerate the development of crop cultivation methods leading to agroecological crop protection.


Introduction
pH has been recognized for long as a master variable in plant physiology, from subcellular to cellular and at the whole plant levels [1]. The recognition of the role of redox regulation is more recent, but a wealth of publications in the last decade demonstrate the tremendous progresses achieved by plant physiologists in understanding redox signaling in plants. Reactive oxygen species (ROS), which largely determine redox status in the plant, are no longer seen as toxic by-products of aerobic pathways but as beneficial companions of plants' developmental processes [2]. The consensus view is that redox signaling is intrinsic to many developmental processes and responses to the environment [3][4][5] and that reactive oxygen species (ROS) have a profound influence on almost every aspect of plant biology [6].
ROS are further recognized as important molecules that affect vegetative and pathogenic processes in pathogenic fungi. They are not only secreted during the interaction of host and pathogen but also are involved in tightly controlled intracellular processes [7]. For instance, powdery mildew resistance has been associated to redox balance in pea [8] and the importance of regulation of nitric oxide (NO) and ROS in potatoes' response to pathogens has been highlighted [9]. It is thus expected that a multitude of redox enzymes play an important role in plant immunity regulation [10,11].
Recent research has also pointed to the role of ROS in the functioning of beneficial fungi. For instance, ROS regulate mutualistic interactions between Epichloe festucae and its grass host [12] and play a role in the proliferation of arbuscular mycorrhizae (AM) fungi in plants and in plant/AM interactions [13]. Lucas et al. (2014) also found that a strain of rhizobacteria associated with rice plants primed enzymes related with the detoxification of ROS [14]. Finally, redox signaling has been identified as a possible unifying mechanism for the mode of action of the wide range of chemical-and hormone-based plant bioregulators [15].
Thus, the current increased understanding of ROS acting as signaling molecules has opened new avenues to exploit redox biology for crop improvement required for sustainable food security [16].
However, these works on redox regulations remained largely at cellular and organelle levels. Despites some advances and a growing interest in the understanding of redox mechanisms of plant systemic responses to pathogens and abiotic stresses [4], an understanding of redox regulation at both plant level and field level is still needed. Thus, Noctor et al. [17] stressed the need: (i) to gauge the oxidative stress intensity (or the intensity of the plant response) to establish whether and to what extent a given genetic modification or condition is generating oxidative stress, and (ii) to provide new insight into the mechanisms of oxidative stress and redox regulation in plants. Noctor et al. [17] also emphasized the need for identification of reliable and informative markers, which are crucial to understanding how plants exploit and respond to oxidative stress and how such markers should greatly aid stress diagnostics in agriculture and environmental science.
Developing indicators of plant stress received a lot of attention and various methods have been proposed, as chlorophyll fluorescence [18], photo-oxidative stress markers (including photosynthetic pigments, PSII efficiency, ROS, reactive carbonyl species, antioxidant systems) [19], or digital imaging [20]. Eh (redox potential, measuring the oxidation/reduction state) is often measured at cellular and subcellular levels. However, despites the possible interactions between Eh and pH, studies of redox regulations are generally disconnected from studies of pH, and Eh-pH interactions are ignored Agronomy 2018, 8,209 3 of 24 most of the time [21]. Correcting Eh to pH7 (Eh@pH7) or calculating pe (electronic potential) + pH, which are equivalent notions, could help better define the redox conditions [22].
Thus, it can be expected that measuring both Eh and pH, at plant level, could provide a good indicator of global plant stress.
For agronomists, having such a means to measure directly the reaction of a plant to soil conditions and cropping practices (soil amendment, fertilization, varieties, crop rotation, pesticide application, etc.) could greatly help and speed up the processes of developing new agronomic methods, variety selection, and plant protection approaches.
The aim of the studies presented in this paper was therefore to develop a simple but reliable method for measuring leaf Eh and pH, and produce new insight into redox processes at plant level for agronomists. Rice (Oryza sativa L.) was used as a model plant for redox studies because of its ability to grow in a very wide range of soil redox conditions, under aerobic or anaerobic field conditions. We present here the results of a series of experiments in pots. We first report the temporal variability in Eh-pH of the last fully expanded leaf of the main tiller, according to the time of day, plant age, and across the seasons. We then present the varietal variations for 3 rice varieties, and we finally present the spatial variations of leaf Eh-pH (intra and inter-leaf variations). This knowledge allows us to propose the use of leaf Eh-pH as an indicator of plant health, which could be a useful and powerful tool for agronomist, breeders and plant pathologists.

Soil, Nutrients, Plants, and Water Management
In all experiments but one, the soil used was produced by Africa Rice farm unit, as a thoroughly homogenized clayey soil coming from the "bas-fond" area of the experimental station. The 11 dm 3 pots were filled-in with 8 dm 3 of soil. Only the third season of the "Intra leaf variability" experiment was conducted with a thoroughly homogenized mixture of sandy soil (Arenosol) and organic matter (forest top soil).
In the first experiment of each series, no fertilizer was added. In the next experiments, 1 g of NPK fertilizer (15-15-15) diluted in 1 dm 3 of water was added in each pot the day before sowing. Finally, in the last experiment of each series, Hoagland solution [23] was used instead of NPK fertilizer to prevent any risk of nutrient deficiency (1 dm 3 /pot before sowing, 1 dm 3 /pot 30 days after sowing (DAS)). Hoagland's solution was prepared with Sigma Aldrich (Saint-Louis, MO, USA) products suitable for plant cell culture.
Genetically pure rice seeds were obtained from the Genetic Resources unit of Africa Rice in Cotonou.
Before planting rice seeds pre-germinated for 72 h, soils were first saturated with water by addition of 3 dm 3 /pot for all treatments. In the "lowland management" treatments, pots were kept submerged with 2 cm of water above the soil surface during the complete cultivation period. In the "upland management" treatments, two hours after saturation with water, pots were opened in the bottom by 8 mm diameter holes. They were then watered every day with 1 dm 3 /pot, to prevent any water stress. Soil water content was kept permanently between 18% and 35% on a volume basis, and was between 21% and 27%, close to field capacity, during all the periods of measurement.

Experimental Settings
Three series of experiments were conducted concomitantly, at the Africa Rice research station in Cotonou, Bénin (6 • 25'13" N; 2 • 19'38" E, 12 m above mean sea level), as summarized in Table 1. For all series of experiments, the pots were arranged at 1 × 1 m, directly on the soil, in an open field, and the treatments were randomly arranged in each block (first series) or fully randomized (second and third series).

Diurnal Change Experiments
In the first series of experiments, we studied the diurnal changes of the middle part of the last fully expanded leaf (i.e., the last fully developed leaf) of the main tiller in three rice (Oryza sativa L.) varieties, with a measurement every 30 min for each variety, from 6:00 a.m. (before sunrise) to sunset. Plants were grown outside, in 11 dm 3 plastic pots, in aerobic conditions ("upland management"). For each of the three varieties, five pre-germinated rice seeds were planted per pot. Fifteen pots were planted for each variety, arranged in 5 blocks with 3 pots of each variety. The day of the measurement, plants with the same development stage were selected (with the same number of tillers and leaves per tiller) and measured in a randomized order, with one replication of each variety every 30 min.
This experiment was conducted 4 times (= four "seasons"), with variations in fertilization (No fertilization during the first season, NPK fertilizer in the second and third season, and Hoagland solution in the last season) and varieties. In the first two seasons, three upland rice varieties were used, with an equivalent cycle of 90-95 days from sowing to 50% maturity: In the next seasons, IDSA6 was replaced by another O. sativa sub. japonica upland rice variety: Azucena, frequently used by plant physiologists.
In total, measurements were conducted for 14 days and 3 varieties per day. Measurements on 10 October 2016 and 31 January 2017 were abandoned due to strong storms. Thus, we report here the results for 36 days × varieties.

Intra-Leaf Variability Experiments
In the second series of experiments, we studied the Eh and pH variations within the leaves in four experiments. Plants were grown in 11 dm 3 pots, with 1 plant per pot, and 4 pots for each treatment. For each leaf measured (last fully expanded leaf in the first three seasons, last and third last fully expanded leaves in the last season), the base, the middle, and the tip of each leaf were measured in a randomized order.
Nerica4 variety was grown: • in "upland" conditions, without fertilizer in the first season. • in "upland" (aerobic) and in "lowland" (permanently submerged) conditions with NPK fertilizer in the second season. • in "upland" and in "lowland" conditions with Hoagland solution in the third season.
IR64 variety was grown in lowland conditions and Hoagland solution in the last season.

Inter-Leaf Variability Experiments
In the third series, we studied the spatial variability according to leaf position on the plant, measuring Eh and pH in the middle part of each leaf, for each tiller. Plants were grown in 11 dm 3 pots, with one plant per pot to avoid competition for light. Six pots were planted for each treatment and measurement date, and randomly arranged outside in a field. Measurements were conducted plant after plant. All the leaves of one plant were measured in a randomized order and all the measurements were conducted between 11:00 a.m. and 4:00 p.m. This experiment was conducted four times in pots, using Nerica4 variety, with some variations in the growing conditions: In the first and the second seasons, the experiments were conducted in aerobic conditions ("upland management"), without fertilizer during the first season and with NPK application during the second season.
In the third season, as Husson et al. [24] have shown that soil Eh-pH gradient could be inverted by cropping practices. The experiment was conducted in aerobic conditions in two treatments creating inversed Eh and pH gradients in soil, to assess the impact of Eh-pH soil gradient on plants. Instead of mixing the arenosol and the organic soil in the pots, a "natural soil gradient" was created by filling the bottom half of the pots with the arenosol and the upper part with the organic soil, and a "reversed soil gradient" was created doing the opposite (organic soil in the bottom and arenosol in the upper half). To prevent nutrient deficiency, Hoagland solution was applied in both treatments.
In the fourth season, the experiment was conducted with Nerica4 variety grown in anaerobic conditions ("lowland management") and a lowland variety (IR64, O. sativa sub. Indica, well-known by physiologists) grown both in aerobic ("upland management") and in anaerobic ("lowland management") conditions, with application of Hoagland solution.
In a fifth season, Nerica4 variety was grown in a field, in upland conditions at Africa Rice Station in Cotonou, on an arenosol, with application of 200 kg/ha NPK fertilizer (15-15-15) at sowing.

Soil Eh and pH
Soil Eh was measured with WTW (Weilhem, Germany) ProfiLine pH 3110 pH meters (>5 × 10 12 Ohm input resistance), using radiometer analytical Ag/AgCl (KCl 3M) reference electrodes Ref 321 and custom-made Paleo terra (Amsterdam, The Netherlands) platinum electrodes. These electrodes were 30 cm long, 8 mm in diameter, and made from fiberglass-reinforced epoxy with 3 platinum rings (width: 0.4 mm; surface: 5 mm 2 ) inserted at 5 cm distance into the pot soil. Two electrodes were directly inserted into the soil, to locate the platinum rings at 3, 8, and 13 cm deep. For each depth, the two electrodes were connected in parallel to a WTW ProfiLine pH 3110 m, and a reference-electrode connected to the voltmeter was inserted into the soil surface. Measurements were taken after stabilization of the voltage (less than 1 mV drift in one minute). The electrodes were cleaned and tested each day, before and after each set of measurements, according to the procedure indicated by Husson et al. [22] with a Hach (Loveland, Colorado, CO, USA) 220 mV redox buffer solution (ref LZW9400.99), and in the same solution diluted (1%) in 0.1 M KCl solution.
All redox potential data presented are actual Eh (mV) vs. standard hydrogen electrode (SHE), The measured E (mV vs. Ag/AgCl 3M) values have been corrected to Eh (mV) as a function of the temperature through the equation: Eh (mV) = E (mV) + 225.8 − 0.73 t ( • C) derived from the table provided by Radiometer Analytical (Vaux-en-Velin, France) for the REF 321 Ag/AgCl KCl 3M electrode.
Soil pH and temperature were taken with a Hach 5014T pH combined electrode (LZW5014T.97.002) directly inserted in the humid soil, and connected to WTW 3110 pHmeter. Calibration of the pH meter was performed every day, with Hach pH buffer solution at pH 7.00 (LZW9464.98) and pH 4.01 (LZW9463.99).
Soil humidity was measured with a WET 2 sensor connected to a Delta HH2 moisture Meter (Delta-T Devices, Cambridge, UK), calibrated for "clayey" soil.

Leaf Eh and pH
The method for measuring leaf Eh has been adapted from Benada [25,26]. It aimed at a direct measurement in a living organism, in aerobic conditions, recognizing that respiration and photosynthesis strongly impact Eh and that Eh of living organisms should be measured when they are still alive, in their natural conditions i.e., aerobic conditions in the case of plants.
The choice of the location for conducting the measurements is crucial for leaf Eh measurement, as it is influenced by electromagnetic waves as shown for soil Eh measurement by Husson et al. [22]. These effects are even stronger when measuring in plant leaves, having a high water content and being less poised than soils. All measurements were therefore conducted outdoor, in an environment identified as being free of electromagnetic interference.
As for soil measurement, the electrodes were cleaned and checked before each set of measurements, in the 220 mV buffer solution and in its dilution at 1% in 0.1 M KCl.
A filter paper was placed in a Petri dish and watered with 0.1 M KCl solution (used as an electric bridge). The reference electrode was placed vertically on a support, standing with the porous plug on the filter paper. The middle part of the leaf was taken from the living plant and rapidly (within 1 min) rolled in order to make several layers, sufficiently thick to allow penetration of the whole platinum plate (5 mm, requiring usually twelve layers for rice leaf). The platinum electrode was then pinched into the rolled leaf, parallel to the main vein, at equal distance between the main vein and the edge of the leaf. The leaf and the electrode were then placed with the electrode standing vertically on the leaf, placed on the filter paper at approx. 2-3 cm from the reference electrode, as seen in Figure S1. The electric potential between the two electrodes showed a decrease, reaching a plateau within a few minutes, and then started rising. This lower turn point value was taken as the measured value, reflecting the Eh of the living leaf. The platinum electrode was then removed from the leaf and pinched again, on the other side of the main vein of the leaf. In some cases, the value on the voltmeter immediately raised. In such case, the first measurement was considered as the measured E (vs. Ag/AgCl) value for the leaf. In other cases, it showed again a decrease, reaching a plateau, usually a few mV lower than the first one. In this case, the lowest of the two measurements was considered as the E (vs. Ag/AgCl) of the plant.
All the measured values of E (vs. Ag/AgCl) were corrected to Eh (vs. SHE), as for soils. Eh@pH7, which is equivalent to pe + pH, can be used in addition to Eh and pH to better define the redox conditions of a system, characterizing the electron activity independently of the proton activity [22]. Eh@pH7 was calculated as Eh@pH7 = Eh + (pH−7) × R × T × Ln10/F, where Eh is measured in volts, R is the perfect gas constant (8.314471 JK −1 mol −1 ), F is the Faraday constant (96485.3383 C mol −1 ) and T is the temperature in Kelvin.
The leaf was then rinsed with a few drops of distilled water, dried by vigorous shaking and placed into a ceramic mortar. It was gently hand-crushed for 1 min and placed in a 5 mL plastic syringe. Pressing the plunger of the syringe allowed the extraction of a few drops which were used for direct measurement of pH and temperature with an ISFET type pH meter (Hach H160 pH Meter with a PHW37-SS electrode), calibrated each day before measurement with Hach pH buffer solution at pH 7.00 (LZW9464.98) and pH 4.01 (LZW9463.99).
Solar radiation was measured at each sampling time with a LI-250A Radiometer/Photometer and a LIC190RBN02 Quantum captor from LI-COR (Lincoln, NE, USA).

Statistical Analysis
Scatter plots were performed for the relationship between Eh, pH, Eh@pH7 and time or leaf length showing different types of relationships. Quadratic regressions between Eh, pH, Eh@ph7 and time of the day were performed for the different DAS × varieties of measurements. Linear regressions were performed between and Eh, pH, Eh@pH7 and the leaf length for the first three age classes.
Analysis of variance (ANOVA) was conducted on Eh, pH, Eh@ph7, temperature, and solar radiation for comparison of growing seasons, plant age (DAS) and varieties. ANOVAs were followed by a comparison of means using the Ryan-Einot-Gabriel-Welsh test Q (REGW-Q) as it provides the best compromise between the need to have a powerful test and the need to limit the familywise error rate at α, according to Howell (2009) [27]. All statistical analyses were performed with XLSTAT (XLSTAT Premium, v. 2018.5, Addinsoft, Paris, France) [28].

Diurnal Variations
Eh of the last fully expanded leaf of the main tiller showed the same diurnal variations, for all the varieties and for each of the 12 days of measurement. Eh was maximum at the end of the night and decreased during the morning until 10:00 to 11:00 a.m. It was stable in the middle of the day, as seen in Figure 1, until 4:00 to 5:00 p.m., when it started to increase again. Significant quadratic correlations between Eh and time of the day were measured for the 36 days × varieties of measurement, shown in Table S1.1, with R 2 ranging from 0.30 to 0.85, and a median R 2 of 0.65.

Diurnal Variations
Eh of the last fully expanded leaf of the main tiller showed the same diurnal variations, for all the varieties and for each of the 12 days of measurement. Eh was maximum at the end of the night and decreased during the morning until 10:00 to 11:00 a.m. It was stable in the middle of the day, as seen in Figure 1, until 4:00 to 5:00 p.m., when it started to increase again. Significant quadratic correlations between Eh and time of the day were measured for the 36 days × varieties of measurement, shown in Table S1.1, with R 2 ranging from 0.30 to 0.85, and a median R 2 of 0.65.  pH of the last fully expanded leaf appeared to be more stable than Eh and was weakly correlated with time in most of the days (significant correlation in only 13 days × variety), with a R 2 above 0.5 in only 11% of the cases, a median R 2 of 0.27, and a maximum R 2 of 0.71 over the 36 days × varieties of measurement ( Figure 1 and Table S1.1).
Eh@pH7 (Eh corrected to pH 7) showed the same diurnal variations as Eh, with a significant correlation in 30 over 36 days × variety, a maximum R 2 of 0.86 and a median R 2 of 0.51 ( Figure 1 and Table S1.1).
Eh and Eh@pH7 were also correlated with: • temperature, with a median R 2 of 0.48 and 0.39 respectively, as seen in Table S1.2, and • solar radiation, with median R 2 of 0.58 and 0.47 respectively, as seen in Table S1.3.

Plant Age
Plants that had grown for 60 days after sowing (DAS) appeared to be significantly more oxidized (higher Eh) and acidic (lower pH) than plants after 40 DAS, for all the varieties and the different seasons experienced during the experiment, as shown in Figures 1 and 2A and Table 2. Importantly, this cannot be explained by differences in temperature or solar radiation at the day of measurement, as measurements at 40 DAS were carried out at higher temperature and lower solar radiation than at 60 DAS during season III, and the opposite during season IV. Interestingly, the increase in Eh with plant age was not related to an increase in Eh@pH7: depending on cropping season and fertilization, as seen in Figures 1 and 2A pH of the last fully expanded leaf appeared to be more stable than Eh and was weakly correlated with time in most of the days (significant correlation in only 13 days × variety), with a R 2 above 0.5 in only 11% of the cases, a median R 2 of 0.27, and a maximum R 2 of 0.71 over the 36 days × varieties of measurement ( Figure 1 and Table S1.1).
Eh@pH7 (Eh corrected to pH 7) showed the same diurnal variations as Eh, with a significant correlation in 30 over 36 days × variety, a maximum R 2 of 0.86 and a median R 2 of 0.51 ( Figure 1 and Table S1.1).
Eh and Eh@pH7 were also correlated with: • temperature, with a median R 2 of 0.48 and 0.39 respectively, as seen in Table S1.2, and • solar radiation, with median R 2 of 0.58 and 0.47 respectively, as seen in Table S1.3.

Plant Age
Plants that had grown for 60 days after sowing (DAS) appeared to be significantly more oxidized (higher Eh) and acidic (lower pH) than plants after 40 DAS, for all the varieties and the different seasons experienced during the experiment, as shown in Figures 1 and 2A and Table 2. Importantly, this cannot be explained by differences in temperature or solar radiation at the day of measurement, as measurements at 40 DAS were carried out at higher temperature and lower solar radiation than at 60 DAS during season III, and the opposite during season IV. Interestingly, the increase in Eh with plant age was not related to an increase in Eh@pH7: depending on cropping season and fertilization, as seen in Figures 1 and 2A,B, Table 2 and Tables S2, S3.1, S3.2 and S4, Eh@pH7 could decrease with plant age as in season III with NPK fertilizer, or increase as in season IV with Hoagland solution.

Cropping Season
Highly significant differences in Eh and significant differences in pH and Eh@pH7 were measured for the fully expended leaf (last fully expanded leaf) of the Nerica4 rice variety, as seen in Table 3, for plants grown in the same conditions (same variety, soil, fertilization, water management, and plant age) and at different periods of the year (second and third "seasons" of the series). The same trend was observed for IRBLTA-2Pi variety, but differences were not significant. Temperature was not significantly different between cropping seasons, as seen in Table 3, while solar radiation was significantly different. Higher pH and lower Eh and Eh@pH7 were observed when solar radiation was high (Table 3  and Tables S5.1 and S5.2). The same trends in Eh and pH were measured, with significant differences, for: (i) the two varieties tested during the second, third and fourth seasons (Table S5.1), and (ii) the three varieties tested during the third and fourth seasons (Table S5.2) although the fertilization during the fourth season was different: NPK in seasons II and III was replaced by Hoagland's solution [23] in the fourth season.

Varietal Differences in Eh, pH and Eh@pH7
The measurements made on the last fully expanded leaf of the main tiller between 11:00 a.m. and 4:00 p.m., during a period of low fluctuation of the 3 parameters, were used to compare varieties. They showed significant differences in Eh, pH and/or Eh@pH7, shown in Figure 2A,B, Table 4 and  Table S4 for all cropping seasons and plant ages.
Nerica4 had a significantly lower Eh than the other varieties. IRBLTA-2Pi had a higher pH than Azucena. Nerica4 had a significantly lower Eh@pH7 than Azucena, and IRBLTA-2Pi had significantly the highest Eh@pH7 (Figure 1; Table 4 and Table S4).

Intra-Leaf Variability
In most cases, the base of the leaf had higher Eh and pH than the tip of the leaf, as seen in Figure 3 and Table S2, and the middle part being most of the time at an intermediate Eh and pH. This was observed over the four seasons, on Nerica4 (upland variety, Oryza sativa type japonica x O. glaberrima) and IR 64 (lowland variety, Oryza sativa sub. indica), in aerobic ("upland" water management) and anaerobic ("lowland" management, with submersion) conditions, for the last fully expanded leaves of all primary tillers, and for older leaves in the main tiller, as seen in Table S2.   Eh of the last (youngest) leaf was high and highly variable, as seen in Figure 4A, and correlated with leaf length. Eh of the last and second last fully expanded leaves (second and third leaf from the top of the tiller) was the lowest in the plant. Eh of the photosynthetically active leaves were either equal, as seen in Figures 5 and 6A, or, in most cases, increasing from the last fully expanded leaf to older leaves, from the second leaf from top to bottom of the plant, with a strong increase in old, dying leaves, as shown in Tables S3.1-S3.7. Eh of the last (youngest) leaf was high and highly variable, as seen in Figure 4A, and correlated with leaf length. Eh of the last and second last fully expanded leaves (second and third leaf from the top of the tiller) was the lowest in the plant. Eh of the photosynthetically active leaves were either equal, as seen in Figures 5 and 6A, or, in most cases, increasing from the last fully expanded leaf to older leaves, from the second leaf from top to bottom of the plant, with a strong increase in old, dying leaves, as shown in Tables S3.1-S3.7.       pH showed a very clear and significant pattern, decreasing from top to bottom of the plant. pH was always maximum in the youngest leaves, and minimum in the oldest leaves, as seen in Figures  5 and 6B and Tables S3.1-S3.7.
Eh@pH7, as a consequence, was significantly the highest. It was correlated with leaf length in the youngest leaf, as seen in Figure 4B. It was slightly decreasing from top to bottom (from young to older leaves) in most of the cases, as shown in Figures 5 and 6C, and Tables S3.1-S3.7.
The same trends were observed for all seasons and fertilization levels, in all tillers, and for all plant ages, as shown in Tables S5.1-S5.7. Interestingly, leaves of the same age on a single rice plant showed similar Eh-pH in all the primary tillers, as seen in Figure 5 and Tables S3.1-S3.7, indicating that leaf Eh-pH is related to leaf age.
These gradients in leaf Eh and pH were not related to Eh-pH gradients in soil horizons, as reversing the soil gradient did not affect Eh-pH gradients in leaves, as recorded in Tables S3.4 and S3.5.
The same gradients were also measured for upland (Nerica4) and lowland (IR64) varieties, in aerobic ("upland" water management) and anaerobic ("lowland" management, with submersion) conditions as shown in Table S3.6.
Finally, rice plants grown in natural field conditions presented the same Eh-pH gradients as those measured in pots, as seen in Table S3.7. pH showed a very clear and significant pattern, decreasing from top to bottom of the plant. pH was always maximum in the youngest leaves, and minimum in the oldest leaves, as seen in Figures 5 and 6B and Tables S3.1-S3.7.
Eh@pH7, as a consequence, was significantly the highest. It was correlated with leaf length in the youngest leaf, as seen in Figure 4B. It was slightly decreasing from top to bottom (from young to older leaves) in most of the cases, as shown in Figures 5 and 6C, and Tables S3.1-S3.7.
The same trends were observed for all seasons and fertilization levels, in all tillers, and for all plant ages, as shown in Tables S5.1-S5.7. Interestingly, leaves of the same age on a single rice plant showed similar Eh-pH in all the primary tillers, as seen in Figure 5 and Tables S3.1-S3.7, indicating that leaf Eh-pH is related to leaf age.
These gradients in leaf Eh and pH were not related to Eh-pH gradients in soil horizons, as reversing the soil gradient did not affect Eh-pH gradients in leaves, as recorded in Tables S3.4 and S3.5.
The same gradients were also measured for upland (Nerica4) and lowland (IR64) varieties, in aerobic ("upland" water management) and anaerobic ("lowland" management, with submersion) conditions as shown in Table S3.6.
Finally, rice plants grown in natural field conditions presented the same Eh-pH gradients as those measured in pots, as seen in Table S3.7.

Spatial and Temporal Variability
Photosynthesis is the fundamental reduction reaction and thus the electron accumulation process at the origin of the energetic functioning of all living organisms on earth. It affects both plant Eh and pH.
As it could be expected, the spatial and temporal variations measured in leaf Eh, pH and Eh@pH7 are in accordance with the well-known variations in photosynthetic activity: Eh, pH and Eh@pH7 (and their variability) were maximum in the top, youngest leaf, in which photosynthetic activity had not fully started yet, and Eh was maximum in the lowest, oldest leaves, in which photosynthesis was probably not fully functional any longer. Inversely, the lowest Eh and Eh@pH7 were observed in the youngest fully developed and photosynthetically active leaves, intercepting high solar radiation and in which photosynthetic activity is expected at its highest intensity [29].
The lower Eh, pH, and Eh@pH7 in the tip of the leaf might also be attributed to higher photosynthetic activity, this part of the leaf capturing more light than the base, which is more in a vertical and shaded position [29].
Similarly, the temporal variations measured in leaf Eh, pH and Eh@pH7 are in accordance with the variations in photosynthetic activity: On a daily basis, Eh and Eh@pH7 decreased when two major factors impacting photosynthesis, light and temperature [29], increased.
Variations across seasons might also be explained by differences in photosynthetic activity. Lower Eh, pH, and Eh@pH7 have been measured when solar radiation and temperature were higher, as seen in Table S5.2.
To confirm the relation of photosynthetic activity and Eh-pH, the stomatal conductance, another major determinant of photosynthetic rate in rice [30], should be measured in future experiments.
Interestingly, aging seems to be related to oxidation and strong acidification of plant leaves. However, despite a global trend always showing a decrease in pH with age, increase in Eh is not always significant, and as a consequence, changes in Eh@pH7 are variable. Variations in leaf Eh, pH, and Eh@pH7 with age also seem to depend on the environment (soil and climate) and cropping practices (water management, fertilization, and cropping season). This can be seen when comparing Eh and pH across the different experiments and seasons of our study. Breaking down the complexity of genotype × environment × management (G × E × M) interactions [31] into variations of Eh-pH-Eh@pH7 depending on genotype, environmental factors and cropping practices would be an extremely useful tool to decipher these complex dynamics and patterns.

Validation of The Measurement Method
The measurement method can certainly be improved (especially reducing the variability). However, the repeatability of the results, the moderate variability, and the significance of the differences measured in space, in time, and between varieties indicate that the method developed here can confidently be used and developed in further studies.
This study also provides essential information to define the conditions of measurement for these parameters to be used in various disciplines like agronomy, plant breeding, or pathology. It confirmed the importance to measure Eh-pH in living plants and to assure that the measuring takes place in an environment free of electromagnetic fields interfering with the measurement.
More important, this study showed that: • As photosynthesis and respiration strongly affect redox level in plants, the measurements should be done in a time window during which solar radiation and temperature variations, hence when Eh variations, are minor: between 11:00 a.m. and 16:00 p.m. in our experimental conditions. • Because of the high spatial and temporal variability of leaf Eh-pH, the measurement should be always conducted on the same part of the same leaf: in the case of upright monocotyledon plants such as rice, we propose to use the middle part of the last fully expanded leaf, as it has the lowest Eh of the plant and the second highest pH, and presents a low variability with time. The last, youngest leaf, which is in rapid development, is not recommended for measurement as its Eh is very variable. Old, senescent leaves are not recommended either, as their Eh also evolves rapidly with leaf age and thus leads to very variable measurements.

•
As plant age impacts leaf Eh and pH, this parameter should also be taken into account when comparing plants.

Leaf Eh-pH and Eh@pH7 as an Indicator of Plant Stresses
Eh and pH regulation is central in plant physiology and phenology and the vast majority of biotic and abiotic stresses are translated into redox signals in plants [21]. Therefore, it can be hypothesized that a mean Eh-pH value at leaf level can be a good indicator, integrating and synthetizing the various stresses as well as beneficial factors faced by the whole plant.
The method proposed here simply aims at measuring such a "mean" Eh and pH at the "coarse" scale of plant organs, not taking into account the important compartmentation of pH and Eh [3]. This study shows that, despites the apparent limitations, a mean value of Eh-pH at plant and organ level can bring important information for agronomy and related disciplines and can be proposed as an indicator of plant stress. It can bring new perspectives, coming as complementary information and adding to the wealth of knowledge on molecular redox processes at cell and organelle levels developed by physiologists and pathologists in recent years [17,32].
Interestingly, differences between varieties were often more significant according to Eh@pH7 than Eh or pH separately (Table 4 and Table S4). This confirms that Eh and pH are interacting in plants. Unfortunately, most of the studies on redox and pH in plants have been conducted independently up to now. Integrating oxidation-reduction and acid-base conditions in future studies would certainly help in deciphering these complex, interacting processes.
It would be interesting in further studies to assess the impact of various kinds of biotic and abiotic stresses on leaf Eh and pH and relate these indicators to stress intensity. Overproduction of reactive oxygen species (ROS) in plants is induced by various environmental perturbations like drought, heat, high light intensity, salinity, osmotic stress, chilling, herbicides, heavy metals, pathogens, wounding, ozone (O 3 ) and atmospheric pollutants [33][34][35]. It can be expected that these ROS affect Eh and pH at the plant level. Could Eh-pH be a unifying indicator, integrating a multiplicity of plant stresses? If so, it would be a powerful tool to analyze genotype × environment × management × bioagressors interactions (G × E × M × B).
Next studies could also use leaf Eh-pH measurements to identify rapidly beneficial effects of e.g., soil amendments, fertigation blends, plant defense elicitors, pest treatments or new land management practices on plant health and growth.
Follow-up studies could also link the variations of Eh and pH to the physiological changes within the leaf during ageing, and how Eh-pH signals interact with programmed cell death [36].
Furthermore, extending the cropping area of a given plant species over a range of latitudes and climatic zones has been a major challenge during plant domestication [37]. Follow-up studies on the relationship between leaf Eh-pH levels at different latitudes and the genetic diversity or plasticity within plant species would deepen our understanding of the adaptability of plants to climate. It could bring new insights in the context of adaptation to climate change.
Finally, plant and canopy architecture influence growth and response to stress [38]. Follow-up studies on monocotyledons with contrasted plant stand (such as rosette-like) or dicotyledons would deepen our understanding of the impact of plant architecture on plant growth.

Spatial and Temporal Variations of Leaf Eh-pH and Disease Patterns and Dynamics
It is striking that the pattern and dynamics of Eh and pH in plants and leaves revealed in this study are in accordance with the pattern and dynamic of development of the major diseases in rice.

Fungi
Fungi in soil exhibit wider pH range for optimal growth than bacteria [39]. However, in plant many pathogenic fungi develop preferably at a slightly acidic pH and seem to be hindered by low Eh and favored by high Eh [21,40]. Rice blast for instance has a good growth on medium at pH ranging from 6 to 7, with an optimum growth at pH 6.5. Its growth is reduced at very low pH [41].
It is remarkable that rice susceptibility/resistance to blast is in accordance with leaf Eh-pH: • Roumen [42] reports that at the 6th leaf stage, blast lesion density was lower in one week-old leaves than in younger, 1-2 days-old leaves, i.e., blast development was higher in the high Eh conditions of immature leaves than in photosynthetically active leaves, at lower Eh.

•
The decline in blast lesion density with leaf age was faster under higher radiation and temperature conditions, i.e., lower Eh and under more active photosynthesis [42]. • Susceptibility of rice to blast decreases from the vegetative stage to the reproductive stage [43], which could be related to decreasing leaf pH (going below 6) and Eh@pH7.
The fact that fungi develop their appressoria in the early morning hours as observed by Benada [42] for Erysiphe graminis might also be related to the high redox potential at this period of the day, when photosynthesis did not yet set in again.

Viruses
Several studies have shown that virus development is slowed down by low soil and plant pH, and favored by high pH and high Eh [21]. From our study, at the plant level, high Eh and pH (and thus high Eh@pH7) are found in rather young plants, and in the youngest leaf. It is also remarkable that the rice yellow mottle virus (RYMV) inoculation is the most efficient on young rice plants and disease severity is the highest in the flag leaf as compared to older leaves [44,45].

Bacteria
Bacteria have a narrow pH range for optimal growth [46]. Pathogenic bacteria seem to be favored by neutral pH. For instance, optimum pH for Xanthomonas oryzae pv. oryzae is 7.0 [47]. This Gram-negative bacteria is obligate aerobic, using oxygen as terminal electron acceptor [48] but as many pathogenic bacteria is sensitive to ROS, the major defense of the host plant [49,50].
In rice leaves, neutral pH with a moderate Eh are found in young plants but not in the oldest leaves where both pH and Eh are lower. Again, this is in accordance with the resistance of rice to bacterial blast (Xanthomonas campestris pv. oryzae), which increases considerably with plant age, and is higher in older, mature leaves [51].
These examples suggest that leaf Eh and pH, reflecting the physiological state of plants, could explain the variable resistance of cereals to obligates parasites, and especially (i) the disease gradients on plants, (ii) the change of susceptibility of organs during the ontogeny and growth, (iii) the difference in resistance in individual plant cells and (iv) the rapid change of resistance in a couple of hour [52].
The possible role of Eh-pH homeostasis in plant health is also reflected in the importance of major rice diseases in different rice ecosystems [53]: rice blast is found mainly in upland conditions (with often oxidized and acidic soil conditions) and more rarely observed in irrigated conditions, especially in Sahelian environment (with high light intensity). Inversely, RYMV and bacterial blight are dominant in irrigated and lowland environment, i.e., in submerged conditions, in which soil pH is increased and soil Eh is decreased by submersion [54].

Eh-pH Homeostasis as a Mechanism of Resistance/Tolerance to Major Diseases?
This accordance between spatial and temporal variations of leaf Eh-pH and disease patterns and dynamics suggests that Eh-pH homeostasis could be an overlooked mechanism of resistance/tolerance to major diseases. Evidence of the role of redox regulation in host plant-pathogen interactions is accumulating [55][56][57][58]. Sclerotinia sclerotiorum (via oxalic acid) generates a reducing environment in host cells that suppress host defense responses including the oxidative burst [59]. The genome of the generalist plant pathogen Fusarium avenaceum is enriched with genes involved in redox signaling [60].
In the rice blast fungus Magnaporthae oryzae, robust anti-oxidant defenses confer tolerance to the host oxidative burst [61]. Magnaporthe grisea also undergoes an oxidative burst of its own during plant infection, which is associated with its development of appressoria [62,63]. The same authors report that ascorbic acid, a scavenger of oxygen radicals, significantly delayed spore germination and appressorium formation in a dose-dependent manner. All these authors mention that these processes are particularly important in the early stage of infection, in the apoplasm. pH also impacts M. oryzae's biology. Conidia production was 5-fold higher at pH 8 than at pH 5, while conversely, conidia germination decreased from pH 8 to pH 3 [64]. As for Eh, M. oryzae also modulates the pH of its environment [64].
Furthermore, virulence of M. oryzae has also been related to Eh-pH regulation [64,65]. The stability of RYMV is pH dependent, with a swollen, unstable form of the virus at basic pH [66]. Furthermore, reversible protein modification by redox alterations could be a key regulatory mechanism for the control of RNA silencing and such a RYMV-encoded viral suppressor of RNA silencing protein (P1) has been identified, with redox dependent flexibility [67].
A final example of how Eh-pH regulation may be an important factor of control of pathogens is the protection of rice against Xanthomonas oryzae by salicylic-acid induced accumulation of superoxide anion [68].
All these studies brought together suggest that each plant pathogen is adapted to specific Eh-pH conditions. They may be able to develop only at certain stages of plant development and when plants are unbalanced and not able to sustain Eh-pH conditions or modify them to prevent the pathogen development.
Similarly, oxidation-reduction processes are also involved in interactions between herbivorous insects and plants [21] and it can be hypothesized that plant Eh-pH conditions affect insects attacks/plant resistance as recently shown for aphids [69].

Leaf Eh-pH and Eh@pH7 as an Indicator of Plant Health
Thus, it is proposed to use leaf Eh-pH not only as an indicator of plant stress, but also as an indicator of plant health. A global Eh-pH perspective at the organ and plant levels could therefore bring new insights for agronomists, pathologists, entomologists, and epidemiologists. For instance, this Eh-pH perspective may help in deciphering:

1.
the spatial and temporal changes in host phenology related to plant and leaf age as observed by Farber and Mundt [70] and the induction of resistance to disease during plant development [71]; 2.
the complex interactions between Genotype, Environment and Management practices and the development of pests and diseases, including the interactions between plant nutrition and diseases [72][73][74] and the interactions between light and diseases [75]; 3.
the variable effects of soil amendment, especially biochar addition on plant diseases [76], and how biochar addition influences severity of diseases [77,78]; 4.
the co-infection and pathogen-pathogen-host interactions, which has been recently recognized as deserving more attention [79], and the suppression of diseases by growth promoting bacteria [80]; 5.
the trade-off between disease resistance and plant fertility [81] and between disease resistance and abiotic stress tolerance [82]; 6.
the processes of co-evolution of plants and their pathogens [83].

Conclusions
Agronomists could use plant Eh-pH to rapidly assess the impact of cropping practices on plant health and relate it to soil health. This could be a very useful tool for the design of innovative cropping systems and the steering of cropping practices.
An Eh-pH approach could also bring original perspectives for breeders. For instance, breeding for resistance/tolerance to specific pathogens could be facilitated by creating environments favoring the development of these specific pathogens through management practices and thus creating a high pressure for selection. It can also be a useful variable to measure to further optimize the inoculation efficiency by adapting the location and time of inoculation to specific pathogens in addition to the current empirical knowledge. Such an Eh-pH approach might also make it possible to develop methods for early and rapid screening of varieties for their resistance/sensitivity to major pathogens.
A major challenge for the future is therefore to develop an intelligence of how environment and practices, particularly through soil characteristics, interact with genotype and impact plant Eh and pH and, as a consequence, pathogen development. This would pave the way towards agroecological crop protection driven by Eh and pH, reducing the environmental impact of agriculture, which is essential [84].