# Dynamics of Foliar Responses to O3 Stress as a Function of Phytotoxic O3 Dose in Hybrid Poplar

Benjamin Turc, Pierre Vollenweider, Didier Le Thiec, Anthony Gandin,

Marcus Schaub, Mireille Cabané, Yves Jolivet

► To cite this version:

Benjamin Turc, Pierre Vollenweider, Didier Le Thiec, Anthony Gandin, Marcus Schaub, et al. Dynamics of Foliar Responses to O3 Stress as a Function of Phytotoxic O3 Dose in Hybrid Poplar. Frontiers in Plant Science, 2021, 12, pp.679852. 10.3389/fpls.2021.679852. hal-03981723

HAL Id: hal-03981723

https://hal.univ-lorraine.fr/hal-03981723

Submitted on 10 Feb 2023

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

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

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

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

Check for
updates

# Dynamics of Foliar Responses to O₃ Stress as a Function of Phytotoxic O₃ Dose in Hybrid Poplar

Benjamin Turc¹,²*, Pierre Vollenweider², Didier Le Thiec¹, Anthony Gandin¹, Marcus Schaub², Mireille Cabané¹ and Yves Jolivet¹

¹University of Lorraine, AgroParisTech, INRAE, SILVA, Nancy, France, ²Section Forest Dynamics, Swiss Federal Institute for Forest, Snow and Landscape Research WSL, Birmensdorf, Switzerland

With background concentrations having reached phytotoxic levels during the last century, tropospheric ozone (O₃) has become a key climate change agent, counteracting carbon sequestration by forest ecosystems. One of the main knowledge gaps for implementing the recent O₃ flux-based critical levels (CLs) concerns the assessment of effective O₃ dose leading to adverse effects in plants. In this study, we investigate the dynamics of physiological, structural, and morphological responses induced by two levels of O₃ exposure (80 and 100 ppb) in the foliage of hybrid poplar, as a function of phytotoxic O₃ dose (POD₀) and foliar developmental stage. After a latency period driven by foliar ontological development, the gas exchanges and chlorophyll content decreased with higher POD₀ monotonically. Hypersensitive response-like lesions appeared early during exposure and showed sigmoidal-like dynamics, varying according to leaf age. At current POD₁ SPEC CL, notwithstanding the aforementioned reactions and initial visible injury to foliage, the treated poplars had still not shown any growth or biomass reduction. Hence, this study demonstrates the development of a complex syndrome of early reactions below the flux-based CL, with response dynamics closely determined by the foliar ontological stage and environmental conditions. General agreement with patterns observed in the field appears indicative of early O₃ impacts on processes relevant, e.g., biodiversity ecosystem services before those of economic significance – i.e., wood production, as targeted by flux-based CL.

**Keywords:** ozone, poplar, hypersensitive response-like, accelerated cell senescence, foliar response

## INTRODUCTION

The ground-level concentrations of ozone (O₃) have increased during the past century (Maas and Grennfelt, 2016), and are predicted to remain stable or increase during the 21st century (Revell et al., 2015; Fu and Tian, 2019). They have already reached levels negatively affecting crop plants and the natural vegetation (Wittig et al., 2009; Jolivet et al., 2016; Proietti et al., 2016; Li et al., 2017), and steady or increasing impacts are expected over the course of next decades (Karlsson et al., 2017).

**OPEN ACCESS**

**Edited by:**
Nam-Chon Paek,
Seoul National University,
South Korea

**Reviewed by:**
Neil J. Willey,
University of the West of England,
United Kingdom
Ling Li,
Mississippi State University,
United States

**\*Correspondence:**
Benjamin Turc
benjamin.turc@univ-lorraine.fr

**Specialty section:**
This article was submitted to
Plant Physiology,
a section of the journal
Frontiers in Plant Science

**Received:** 25 March 2021

**Accepted:** 03 June 2021

**Published:** 28 June 2021

**Citation:**
Turc B, Vollenweider P, Le Thiec D, Gandin A, Schaub M, Cabané M and Jolivet Y (2021) Dynamics of Foliar Responses to O₃ Stress as a Function of Phytotoxic O₃ Dose in Hybrid Poplar. *Front. Plant Sci.* 12:679852.
doi: 10.3389/fpls.2021.679852
Once entering the leaf through stomata, O₃ degradation causes the formation of reactive oxygen species (ROS), the accumulation of which triggers rapid oxidative bursts (Schraudner et al., 1998; Pasqualini et al., 2003; Moura et al., 2018). ROS can also act as elicitors of programed cell death (PCD) reminiscent of plant responses during defensive plant/pathogen interactions which are subsequently designated as hypersensitive response-like (HR-like; Vollenweider et al., 2002; Bhattacharjee, 2005; Günthardt-Goerg and Vollenweider, 2007; Moura et al., 2018). In parallel, an acceleration of cell senescence (ACS), with distinct apparent mechanisms, can be observed (Pell et al., 1999; Günthardt-Goerg and Vollenweider, 2007; Vollenweider et al., 2019). The characteristic symptoms thus include marked degenerative injuries in chloroplasts, in apparent relation to an increase in the constitutive ROS load resulting from the daily photosynthetic activity. As a consequence, these latter organelles are particularly sensitive to O₃ stress (Joo et al., 2005; Kangasjarvi et al., 2005). However, the sequence of plant reactions in response to O₃ stress remains unclear, especially given the driving – but still partially understood – effects of interacting environmental conditions and ontological development. In field vs. climate chamber conditions, for example, the high vs. low-intensity illumination can lead to contrasted symptom expression, with clear synergies between photooxidative and O₃ stress in the former case only (Günthardt-Goerg and Vollenweider, 2007; Paoletti et al., 2009; Moura et al., 2018; Vollenweider et al., 2019). Hence, the dynamics of responses to O₃ stress as a function of environmental conditions needs further research.

Although the effects of O₃ stress have been observed in both mature and developing foliage, their intensity is strongly related to the leaf ontology, the mature leaves being more sensitive than those still in expansion. However, the younger vs. older leaves can show higher rates of stomatal conductance and O₃ uptake (Reich, 1983; Strohm et al., 1999; Bagard et al., 2008; Zhang et al., 2010; Guerrero et al., 2013), suggesting an enhanced detoxification capacity (Bellini and De Tullio, 2019). Still, the mechanisms underlying the higher O₃ tolerance in developing foliage remain largely obscure (Strohm et al., 2002) and the differences in response dynamics as a function of leaf ontogenetic development require further investigations.

To assess and prevent O₃ injury on vegetation and forest trees, a concentration-based index, namely, the accumulated O₃ exposure threshold over 40 ppb (AOT40), was initially proposed (Fuhrer et al., 1997). Given the dependency of O₃ phytotoxicity on stomatal conductance, the biologically (Karlsson et al., 2007; Mills et al., 2011) and environmentally (Musselman et al., 2006; Dizengremel et al., 2013; Büker et al., 2015) more relevant flux-based approach has been increasingly implemented. Nowadays, the O₃ critical level (CL) is defined for given vegetation types or plant species and calculated as the Phytotoxic O₃ Dose over a Y threshold for a specific species or group of species (POD<sub>Y SPEC</sub> (Mills et al., 2017). Based on empirical evidence from risk assessment studies – linking POD<sub>Y SPEC</sub> values to tree biomass loss or foliar injury – the current CL typically targets 4% maximum, i.e., growth reductions by oxidative stress. However, such markers represent some late

O₃ stress effects, at least partly resulting from earlier processes in foliage (i.e., reduced physiological activity/extensive cellular injury) which dynamics primarily depends on detoxification processes (Dghim et al., 2013; Dumont et al., 2014; Dusart et al., 2019a). With a view to the larger implementation and acceptance of flux-based approach, there is then an important knowledge gap regarding the dynamics and effective PODₓ of first effective O₃ stress effects, prior to the appearance of current risk assessment markers.

In this study, our main objective was to characterize the dynamics of early physiological and structural responses to O₃ stress in poplar trees as a function of flux-based O₃ dose and before, e.g., growth reduction and extensive foliar injury, the primary markers of O₃ stress for defining O₃ CL (Sanz and Catalayud, 2011; Mills et al., 2017). The tested hypotheses (H) included: (H1) the development of injury and growth response to O₃ stress, as well as physiological and structural changes, proceeds in sequential order, with each response showing specific dynamics; (H2a) O₃ elicits different injury responses within the foliage of trees (H2b) with ACS occurring before the development of HR-like lesions (Vollenweider et al., 2019); (H3a) at comparable O₃ dose and irrespective of the applied O₃ concentration, leaves show similar responses and (H3b) response dynamics; (H4) the dynamics of responses depends on the leaf developmental stage (Moura et al., 2018). Therefore, rooted cuttings of hybrid poplar (Populus tremula x alba) were exposed to three O₃ concentrations in fully controlled conditions for a month. The leaf physiology, development of ACS and HR-like lesions, and appearance of visible injuries were monitored over the course of 29 days. The interaction between foliar response dynamics and leaf ontological development was evaluated by assessing the responses to treatments at two distinct leaf positions.

# MATERIALS AND METHODS

## Plant Material and Controlled O₃ Exposure

Young trees from a hybrid Populus tremula x alba clone (INRAE 717-1b4) were cultivated similarly to Cabane et al. (2004). Before the experimental exposure, micro propagated cuttings were grown for 2 weeks in 0.5 L pots containing compost (Gramoflor Universel) and perlite [1:1 (v/v)], and placed in containers covered with transparent acrylic hoods inside a growth chamber. The environmental conditions were set at 22°C/18°C day/night temperature, 350 µmol m⁻² s⁻¹ photosynthetic active radiation (PAR, 1 m below lamps) during a 14-h photoperiod (Philips Son-T Agro 400 W lamps), 75%/85% relative humidity (day/night). The young trees were then transplanted into 10 L pots filled with compost (Gramoflor Universel) and fertilized with 3 g l⁻¹ of slow-release Nutricot T 100 granules (13:13:13:2 N:P:K:MgO, Fertil, Boulogne-Billancourt, France). They were further cultivated for 1 month in the same growth chamber and watered to field capacity every day. The trees retained (n = 48), with a view to the forthcoming O₃ exposure experiments, were 29.5 ± 0.2 cm high, with 13.1 ± 0.1 leaves. During experiments, all foliar assessments were repeated at the third and tenth leaf position
from the tree base, thus selecting the youngest fully expanded leaf and that still in expansion at treetop by the start of exposure.

Before exposure, the selected trees were randomly distributed among six ventilated phytotron chambers (1 air change min⁻¹; 120 cm × 117 cm and 204 cm high) within the O₃-exposure facility of the PEPLor platform (Faculty of Sciences and Technologies, University of Lorraine). Within each chamber (N = 8 trees), the plant position was randomized by each assessment. The transferred poplars were then left to acclimate for 1 week, with environmental conditions similar to those in the growth chamber. The O₃ exposure experiment included three treatments [charcoal-filtered (CF) air; CF + 80 ppb O₃; CF + 100 ppb O₃] replicated in two chambers each and performed for 30 days (N = 16 trees per treatment). O₃ was generated from pure oxygen using an O₃ generator (Innovatec II, Rheinbach, Germany), and provided to the chambers during the daytime period in the form of a 13 h square wave, starting 1 h after the light was switched on. The O₃ concentrations within each phytotron chamber were monitored twice an hour using a computer-assisted automatic O₃ analyzer (O341M, Environment SA, Paris).

## Dynamics of Leaf Physiology Responses and Estimation of Phytotoxic O₃ Dose

The effect of treatments on the dynamics of leaf gas exchanges was assessed by measuring the net CO₂ assimilation rates (Aₘₐt) and stomatal conductance to water vapor (gₗ) every 2 days, 3 h after starting the O₃ exposure. Selecting six trees per treatment, the measurements were performed using a Li-6400XT portable photosynthesis system (LiCor, Inc., Lincoln, NE, United States), with cuvette temperature set at 22°C, light intensity (PAR) at 300 and 320 µmol m⁻² s⁻¹ for measurements at the third and tenth leaf position, respectively, airflow at 300 µmol s⁻¹, CO₂ concentration at 400 ppm, and leaf vapor pressure deficit (VPDₜₐf) < 1 kPa. The values were recorded once gₗ and Aₘₐt remained stable for 30 s.

The gₗ estimates (mol H₂O m⁻² s⁻¹) were used to calculate the instantaneous O₃ uptake into the leaf under environmental stable conditions (Fₐ₇), according to Bagard et al. (2015):

$$F_{O_3} = [O_3]_{atm} * g_{O_3}$$

with Fₐ₇ as the O₃ flux (nmol O₃ m⁻² s⁻¹), [O₇]ₐₘₐt as the O₃ concentration (ppb) in the phytotron chamber, and gₐ₇ (O₃ m⁻² s⁻¹) as the stomatal conductance to O₃, according to (Lamaud et al., 2009):

$$g_{O_3} = \frac{D_{O_3}}{D_{H_2O}} * g_w$$

with Dₐ₇ and Dₐ₇₀ as the O₃ and water molecular diffusivity (cm⁻² s⁻¹) respectively (Massman, 1998). The hourly O₃ uptake (mmol O₃ m⁻² h⁻¹) was calculated by integrating Fₐ₇ over an hour and the POD₀ (mmol O₃ m⁻²), by cumulating the hourly O₃ uptake since the beginning of experiment. Missing gₗ measurements were estimated based on values from flanking days (Bagard et al., 2008).

The effect of treatments on the dynamics of leaf chlorophyll content was assessed by measuring estimates of surface-based concentrations of chlorophylls (total chlorophyll index) every day, 1 h after switching the light on and before the start of O₃ treatment. Selecting six trees per treatment, the estimates were obtained averaging 10 measurements per leaf performed with a leaf clamp sensor device (Dualex Force-A, Orsay, France).

## Dynamics of Microscopic and Visible Leaf Injury

The development of HR-like lesions within the mesophyll and that of O₃ symptoms throughout foliage was monitored using completing microscopic assessments and visible injury observations. For microscopic assessments, two discs (diameter = 6 mm) per leaf position in two trees per chamber were sampled every 2 days, until HR-like lesions were detected in the 100 ppb O₃ treatment at both leaf positions; the sampling interval was then extended (3–6 days). The harvested discs were processed immediately after sampling.

Necrotic cells within mesophyll as a consequence of HR-like lesions were evidenced using the Trypan blue assay (Pasqualini et al., 2003; Joo et al., 2005; Faoro and Iriti, 2009). Briefly, the leaf discs were stained for 3 min in a hot lactophenol Trypan blue mixture (60 ml staining solution: 10 g phenol, 10 mg Trypan blue, 30 ml ethanol, 10 ml glycerol, 10 ml lactic acid, and 10 ml distilled water) and the necrotic cells contrasted for 20 min against a clear background using 2.5 g ml⁻¹ hot chloral hydrate, before mounting in 60% glycerol (Pasqualini et al., 2003). The preparations were then transferred to WSL where all microscopy assessments were performed. The disk's central part, free of staining artifacts, was observed using the 5× objective of a Leica microscope (Leitz DM/RB). Given the disk thickness (>200 µm) and to create high contrast pictures, the preparations were imaged after inserting the 10x condenser and removing most filters and diaphragms, using the INFINITY 2-1R camera and Lumenera Infinity Analyze (release 6.4) software (Lumenera Corp., Ottawa, ON, Canada). The center of each disk preparation was photographed, creating composite images made of nine tiles each. The percentage area, particle size, and shape properties of HR-like lesions inside of stitched images were quantified using computer-assisted color image analysis (software WinCELL™ 2004, Regent Instruments Inc., Québec, QC, Canada). Briefly, the software attributed the whole lesion or part of it to one of two color classes (non-oxidized: violet hue; oxidized: dark blue hue) made of 10 shades each, defined based on a representative batch of images and contrasting with the background color class (grayish hue, based on 10 white to gray shades). The quantified parameters characterized the size and shape properties of total and individual lesion particles.

The HR-like lesions and oxidation diagnosis were verified based on hand, and semi-thin sections from samples collected in all treatments at the two leaf positions during the whole study and subsequently processed and observed as described previously (Moura et al., 2018). Briefly, supplementary leaf discs were infiltrated upon sampling with EM-grade 2.5%
glutaraldehyde buffered at pH 7.0 with 0.067 M Soerensen phosphate buffer, renewed after vacuum infiltration. Sections (60 µm) obtained using a custom-made hand microtome and kept unstained were used for visualizing chlorophylls and the oxidation of HR-like lesions. Technovit-embedded 1.5 µm sections, obtained using a Supercut Reichert 2050 microtome and stained with Toluidine blue (Vollenweider et al., 2016), were used to identify HR-like markers within necrotic mesophyll cells, after observation with phase contrast illumination in bright field microscopy using the 5–100× objectives of the Leica microscope and imaged using the Infinity camera, as mentioned above.

The development of O₃ injury in response to treatments was surveyed in the morning, before the start of O₃ exposure, daily. Upon appearance, the development of visible injury was monitored with pictures of symptomatic leaves. The percentage area of necrosis per leaf within the latter material was quantified using color image analysis, using the Color Segmentation plugin in Fiji freeware (ver. 2.0.0; Schindelin et al., 2012).

## Morphological Assessments

After 30 days of exposure, all trees were harvested and biometric assessments were conducted. Tree height was recorded and stem diameter 1.5 cm above root collar was measured, using a hand caliper. The number of leaves per tree, shed or still attached, was recorded before harvest. Leaf and stem material was oven-dried to constant weight, before determining the dry mass of each fraction.

## Statistical Analysis

The dynamics of physiological and structural responses to treatments in foliage and the differences in whole-tree morphology and biomass between groups by the end of the experiment were analyzed using linear mixed-effects models (lmem). The fixed-effect factors included the O₃ treatment, leaf position, time or POD₀ and interactions, whereas the tree nested in the chamber (leaf data; with the leaf position as the statistical unit) or the chamber (morphology/biomass data; with the tree as the statistical unit) were introduced in models as random terms. Homoscedasticity and normality of residuals were verified graphically, and the dependent variables were log- or square-transformed to meet the model assumptions, as needed. The differences between treatments at given assessment dates were tested using post-hoc tests (Tukey's honest significant difference). All statistical analyses were performed using R

statistical software, version 3.5.0 (R Development Core Team, 2017), with the packages lme4 (Bates et al., 2015) for linear mixed-effects models, and emmeans (Lenth, 2016) for post-hoc testing.

# RESULTS

## Morphological Responses

After 30 days of exposure, no change in tree height, stem diameter, or foliar dry mass in response to O₃ exposure was observed (Table 1). However, the stem biomass (*p* = 0.014), amount of leaves (*p* = 0.003), and leaf shedding (*p* = 1.2 × 10⁻⁶) were increased, with significant differences between the 80 and 100 ppb O₃ treatments in the case of leaf shedding.

## Dynamics of Stomatal Responses and Changes in the Phytotoxic O₃ Dose

At both leaf positions, the 100 and 80 ppb O₃ treatments significantly reduced *g*<sub>w</sub> (Figure 1A; O₃ treatment: *p* < 0.001) and accelerated its leaf ontology-driven decrease (O₃ treatment*Time: *p* < 0.001). This reduction was delayed at the tenth vs. third leaf position (O₃ treatment*leaf position, O₃ treatment*time*leaf position: *p* < 0.001), with a 50% decrease in *g*<sub>w</sub> reached in 15 vs. 6 days, respectively, in the 100 ppb treatment. As indicated by increasing *g*<sub>w</sub> in maturing leaves (tenth leaf position) during the first 10 days of exposure irrespective of O₃ exposure (leaf position: *p* < 0.001), the O₃ treatment affected *g*<sub>w</sub> only once the ontological development had been achieved (latency phase, Figure 1A). By the end of the experiment and only at the tenth leaf position, the differences in *g*<sub>w</sub> between the 100 and 80 ppb treatments were significant.

By the end of the experiment, trees in the 100 vs. 80 ppb O₃ treatment showed a larger POD₀, as a consequence of their higher O₃ concentrations and mostly similar *g*<sub>w</sub> (Figure 1B; O₃ treatment, O₃ treatment*time: *p* < 0.001). After 30 days of exposure, the POD₀ at the third and tenth leaf positions was thus 1.4 and 1.2 times higher in the 100 vs., 80 ppb O₃ treatments. The POD₀ was also higher in leaves at the tenth vs. third leaf position (leaf position: *p* < 0.001), as a consequence of the delayed leaf ontogeny and higher *g*<sub>w</sub> (O₃ treatment*leaf position: *p* < 0.001). After 10 days of exposure, the POD₀ levels in younger foliage thus exceeded those in older material by approximately 25% and outpaced them by 40% by the end of exposure (O₃ treatment*time*leaf position: *p* < 0.001).

<table><caption>TABLE 1 | Morphological responses to O₃ treatments in hybrid poplar (Populus tremula x alba) at the end of the experiment.</caption><thead><tr><th>Treatment</th><th>Tree height (cm)</th><th>Stem diameter (mm)</th><th>Foliage biomass (g)</th><th>Stem dry mass (g)</th><th>Leaf shedding (%)</th><th>Leaf number</th></tr></thead><tbody><tr><td>Charcoal-filtered</td><td>100.31 ± 0.95</td><td>9.93 ± 0.18</td><td>30.09 ± 0.86</td><td>15.72 ± 0.53a</td><td>1.59 ± 0.81a</td><td>33.31 ± 0.54a</td></tr><tr><td>80 ppb ozone</td><td>104.38 ± 1.66</td><td>10.16 ± 0.33</td><td>31.54 ± 1.50</td><td>17.51 ± 0.91b</td><td>7.41 ± 1.99b</td><td>35.00 ± 0.61b</td></tr><tr><td>100 ppb ozone</td><td>101.75 ± 0.81</td><td>10.23 ± 0.12</td><td>28.83 ± 0.61</td><td>16.50 ± 0.39ab</td><td>16.94 ± 2.21c</td><td>35.44 ± 0.39b</td></tr><tr><td>Treatment</td><td>ns</td><td>ns</td><td>ns</td><td>*</td><td>***</td><td>**</td></tr></tbody></table>

<footnote>model: lmer(variable) - O₃ treatment + (1 | chamber); *** *p* ≤ 0.001; ** *p* ≤ 0.01; * *p* ≤ 0.05; ns, not significantly different. Values represent means ± SE, n = 16. Different letters indicate significant differences between treatments by the end of experiment (post-hoc Tukey's honest significant difference, *p* ≤ 0.05).</footnote>
<figure><img src="image_3.png" /></figure>

<table><thead><tr><th></th><th>g<sub>w</sub></th><th></th><th>POD<sub>0</sub></th></tr></thead><tbody><tr><td>O<sub>3</sub> treatment</td><td>***</td><td>O<sub>3</sub> treatment</td><td>***</td></tr><tr><td>Leaf Position</td><td>***</td><td>Leaf Position</td><td>***</td></tr><tr><td>Time</td><td>***</td><td>Time</td><td>***</td></tr><tr><td>O<sub>3</sub> treatment:Time</td><td>***</td><td>O<sub>3</sub> treatment:Time</td><td>***</td></tr><tr><td>O<sub>3</sub> treatment:Leaf Position</td><td>***</td><td>O<sub>3</sub> treatment:Leaf Position</td><td>***</td></tr><tr><td>Leaf Position:Time</td><td>***</td><td>Leaf Position:Time</td><td>***</td></tr><tr><td>O<sub>3</sub> treatment:Time:Leaf Position</td><td>***</td><td>O<sub>3</sub> treatment:Time:Leaf Position</td><td>***</td></tr></tbody></table>

**FIGURE 1** | Dynamics of changes in the (A) stomatal conductance to water (g<sub>w</sub>) and (B) phytotoxic O<sub>3</sub> dose (POD<sub>0</sub>) of hybrid poplar leaves (*Populus tremula* x *alba*), as a function of O<sub>3</sub> treatment (charcoal-filtered air ●, 80 ppb O<sub>3</sub> ▲, 100 ppb O<sub>3</sub> ■), leaf position, time of assessment and interactions (model: lmer[sqrt(variable)] ~ O<sub>3</sub> treatment * leaf position * time + (1 | tree / chamber); *** * p* ≤ 0.001). Values represent means ± SE, *n* = 6. Different letters indicate significant differences between treatments for a given assessment date (*post-hoc* Tukey's Honest significant difference, *p* ≤ 0.05).

## Dynamics of Leaf Physiology Responses

Irrespective of the leaf position, A<sub>net</sub> showed responses to treatments and response dynamics similar to g<sub>w</sub> (Figure 2A vs. Figure 1A). After 30 days of exposure, A<sub>net</sub> in the 100 ppb O<sub>3</sub> vs. CF treatment was decreased by 70 and 35% at the third and tenth leaf position, respectively. Once the leaf ontogenetic differentiation achieved, A<sub>net</sub> decreased with POD<sub>0</sub> in a constant and monotonic manner (Figure 2B; POD: *p* < 0.001), regardless of the treatment or leaf position. The significant effects of O<sub>3</sub> treatment and O<sub>3</sub> treatment*POD factors (*p* < 0.001) could then be related to the less affected A<sub>net</sub> values in the 80 vs. 100 ppb treatment at the highest POD<sub>0</sub> reached by the end of the exposure. The observed reduction in A<sub>net</sub> as a function of POD<sub>0</sub> was stronger at the third vs. tenth leaf position (Figure 2B; leaf position *p* < 0.001) but the dynamics at both leaf positions was similar (leaf position*POD: ns). Hence, at POD<sub>0</sub> of 5 mmol O<sub>3</sub> m<sup>-2</sup>, A<sub>net</sub> at the tenth leaf position did not show any reduction relative to CF treatment yet, vs. 50% A<sub>net</sub> loss in older leaves, irrespective of the O<sub>3</sub> treatment. The response differences between the two leaf positions were further observed at higher POD<sub>0</sub>. With POD<sub>0</sub> above 9 mmol O<sub>3</sub> m<sup>-2</sup>, as recorded in younger leaves only, and also as a consequence of the aforementioned latency effect (Figure 1B), A<sub>net</sub> never dropped to levels observed at the third leaf position for lower POD<sub>0</sub>. Consequently, the photosynthetic activity in younger vs. older foliage appeared less sensitive to the absorbed O<sub>3</sub> dose.

An O<sub>3</sub> impact on the chlorophyll content index of leaves was detected after 24 days of the experiment. The impact was

restricted to the third leaf position (Figure 3A; O<sub>3</sub> treatment: ns; O<sub>3</sub> treatment*time, O<sub>3</sub> treatment*leaf position: *p* < 0.001), showing a decrease of 30% for the total chlorophyll index in the 100 ppb O<sub>3</sub> vs. CF treatment. These findings primarily related to latency effects due to leaf ontological maturation, which were observed at the third as well as the tenth leaf position in the case of this parameter, lasting 7 and 17 days, respectively. Accordingly, a larger latency peak was observed in younger than older foliage. The total chlorophyll index decreased with higher POD<sub>0</sub> (Figure 3B; POD: *p* < 0.05), irrespective of the O<sub>3</sub> treatment (O<sub>3</sub> treatment: ns, O<sub>3</sub> treatment*POD: ns). At low POD<sub>0</sub>, the decline was rather monotonic, but accelerated with values exceeding 8 mmol m<sup>-2</sup> at the third leaf position, thus contrasting with the nearly linear drop observed in younger leaves. Confirming a higher O<sub>3</sub> tolerance leaves at the tenth leaf position showed smaller (leaf position: *p* < 0.05) and slower (leaf position*POD<sub>0</sub> *p* < 0.001) drops with higher POD<sub>0</sub>.

## Dynamics of HR-Like Lesion Spread and Development

The microscopic necrosis observed in mesophyll using the Trypan blue assay was diagnosed as being caused by HR-like processes (Figures 4A–D), based on several typical O<sub>3</sub>-stress markers (Paoletti et al., 2009; Vollenweider et al., 2019). They included (1) the characteristic intercostal distribution of lesions (Figure 4E), (2) the development of injury first in older leaves (Figure 5A), or (3) the multiple HR-like events restricted to cells or small groups of cells within mesophyll (Figures 4E–G,I vs. Figure 4H). Collapsed dead cells were mainly observed in
<figure><img src="image_4.png" /></figure>

<table><thead><tr><th></th><th>A<sub>net</sub></th><th></th><th>A<sub>net</sub></th></tr></thead><tbody><tr><td>O<sub>3</sub> treatment</td><td>***</td><td>O<sub>3</sub> treatment</td><td>***</td></tr><tr><td>Leaf Position</td><td>***</td><td>Leaf Position</td><td>***</td></tr><tr><td>Time</td><td>***</td><td>POD<sub>0</sub></td><td>***</td></tr><tr><td>O<sub>3</sub> treatment:Time</td><td>***</td><td>O<sub>3</sub> treatment:POD<sub>0</sub></td><td>***</td></tr><tr><td>O<sub>3</sub> treatment:Leaf Position</td><td>***</td><td>O<sub>3</sub> treatment:Leaf Position</td><td>ns</td></tr><tr><td>Leaf Position:Time</td><td>***</td><td>Leaf Position:POD<sub>0</sub></td><td>ns</td></tr><tr><td>O<sub>3</sub> treatment:Time:Leaf Position</td><td>***</td><td>O<sub>3</sub> treatment:POD<sub>0</sub>:Leaf Position</td><td>ns</td></tr></tbody></table>

**FIGURE 2** | Dynamics of changes in the net CO₂ assimilation (A<sub>net</sub>) of hybrid poplar leaves (*Populus tremula* x *alba*), as a function of assessment time (A), phytotoxic O₃ dose (POD₀; **B**), O₃ treatment (charcoal-filtered air ●, 80 ppb O₃ ▲, 100 ppb O₃ ■), leaf position and interactions [model: Imer(variable) ~ O₃ treatment * leaf position * time + (1 | tree / chamber)]; *** *p* ≤ 0.001; ns, not significantly different]. Values represent means ± SE, *n* = 6. Different letters indicate significant differences between treatments at a given assessment date (*post-hoc* Tukey's Honest significant difference, *p* ≤ 0.05).

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

<table><thead><tr><th></th><th>Chlorophyll Content<br>Index</th><th></th><th>Chlorophyll Content<br>Index</th></tr></thead><tbody><tr><td>O<sub>3</sub> treatment</td><td>ns</td><td>O<sub>3</sub> treatment</td><td>ns</td></tr><tr><td>Leaf Position</td><td>***</td><td>Leaf Position</td><td>*</td></tr><tr><td>Time</td><td>***</td><td>POD<sub>0</sub></td><td>***</td></tr><tr><td>O<sub>3</sub> treatment:Time</td><td>***</td><td>O<sub>3</sub> treatment:POD<sub>0</sub></td><td>ns</td></tr><tr><td>O<sub>3</sub> treatment:Leaf Position</td><td>***</td><td>O<sub>3</sub> treatment:Leaf Position</td><td>*</td></tr><tr><td>Leaf Position:Time</td><td>***</td><td>Leaf Position:POD<sub>0</sub></td><td>***</td></tr><tr><td>O<sub>3</sub> treatment:Time:Leaf Position</td><td>***</td><td>O<sub>3</sub> treatment:POD<sub>0</sub>:Leaf Position</td><td>ns</td></tr></tbody></table>

**FIGURE 3** | Dynamics of changes in the surface-based concentration of chlorophylls (total chlorophyll content index of Dualex) within hybrid poplar leaves (*Populus tremula* x *alba*), as a function of assessment time (A), phytotoxic O₃ dose (POD₀; **B**), O₃ treatment (charcoal-filtered air ●, 80 ppb O₃ ▲, 100 ppb O₃ ■), leaf position and interactions [model: Imer(variable) ~ O₃ treatment * leaf position * time + (1 | tree / chamber)]; *** *p* ≤ 0.001; * *p* ≤ 0.05; ns, not significantly different]. Values represent means ± SE, *n* = 6. Different letters indicate significant differences between treatments at a given assessment date (*post-hoc* Tukey's Honest significant difference, *p* ≤ 0.05).

the lower palisade parenchyma (Figures 4I–L). Non-oxidized lesions (Figure 4G) showed up first (Figures 5A,C, 6C), with, for instance, still green chloroplasts visible within collapsed

dead cells (Figures 4I,J). The dark hues of oxidized lesions (Figure 4L) were enhanced by staining with Trypan blue (Figures 4E,G). Oxidized cells showed sharp wall angles indicative
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | Distribution, morphology, and structural properties of HR-like lesions within hybrid poplar leaves (*Populus tremula* x *alba). (A–D) Foliar discs excised at low leaf position (third) from leaves exposed to 100 ppb O₃ during 2 (A), 8 (B), 13 (C), and 27 (D) days. (A) Central area within each disk was photographed each time (frame in D). (E–G) Image analysis of HR-like lesion after 27 days of treatment, as framed in (D). (E) Synthetic digital image, made of 9 stitched micrographs (5x magnification). The lesions are separated by veinlets and non-symptomatic tissues. (F) The binary image of total HR-like lesion vs. background (white). (G) Classification of HR-like injuries into oxidized and non-oxidized lesion groups, based on color classes. (H–L) Changes in the mesophyll tissue and cell structure underlying the HR-like lesions. (H) Asymptomatic leaf tissues in a leaf sample from the filtered air treatment. I–L: necrotic cells within the upper (UPP) and lower (LPP) palisade parenchyma underlying the HR-like lesions. (I,J) Within mesophyll cells having recently undergone HR-like necrosis (#), the chloroplasts (chl) were still visible and had retained their green color (J). (K,L) at a later stage, the HR-like lesions (§) showed cell-content disruption and oxidation (L). Other structures: UE, upper epidermis; SP, spongy parenchyma; LE, lower epidermis; v, vacuole; *, nucleus. Technical specifications: staining with Trypan blue (A–E) and Toluidine blue (H,I,K); observations in bright field microscopy (A–E,H–L) using phase-contrast (H,I,K); (J,L) fresh, unfixed and unstained leaf sample preparations.</figcaption></figure>

of breaks and disrupted cell content (Figures 4K,L). All these typical HR-like traits showed little variation, regardless of the O₃ treatment or leaf position.

The first HR-like lesions at the third and tenth leaf position were observed after 6 and 13 days of exposure, respectively. This was much earlier than reductions in the chlorophyll content index (Figure 5A). Despite large response variability among trees, the effect of O₃ treatment was significant (Figure 5A; O₃ treatment: *p* < 0.05; O₃ treatment*time*: *p* < 0.001). A larger leaf percentage area showing HR-like lesions was observed for the 100 vs. 80 ppb O₃ treatment, with differences between the two treatments at the third leaf position becoming significant after 20 days of exposure (O₃*time*; *p* < 0.001). After 27 days of treatment, the percentage area of lesions in the 100 vs.

80 ppb O₃ treatment was two and five times higher at the third and tenth leaf positions, respectively. However, each leaf position showed specific response dynamics (O₃ treatment*leaf position: *p* < 0.01, O₃ treatment*time*leaf position: *p* < 0.05), rather sigmoidal-like vs. linear – once lesions appeared – in older vs. younger foliage (Figure 5A). Moreover, the HR-like lesions in response to the two O₃ concentrations showed up simultaneously at the third leaf position whereas a 7-day delay was observed at the tenth leaf position (Figure 5A).

When expressed as a function of POD₀, the differences between the two O₃ treatments in the leaf percentage area showing HR-like lesions were leveled out, especially at the third leaf position (Figure 5B; O₃ treatment: ns). However, the dependency on POD₀ was lessened at the tenth vs. third
leaf position (leaf position*POD: *p* < 0.001), and distinctly higher lesion percentage areas in response to similar POD₀ were observed in the 100 vs. 80 ppb O₃ treatment in younger leaves (O₃ treatment*POD*leaf position: *p* < 0.05). Hence, not only the O₃ dose but also the O₃ absorption rate then determined the lesion severity. The oxidized HR-like lesions, expressed as a function of time or POD₀, showed responses and response dynamics similar to HR-like lesions taken globally (Figures 5C,D vs. Figures 5A,B). The main differences included a smaller percentage of the injured area and a weaker symptom dynamics. At the tenth leaf position, oxidized HR-like lesions in the 80 ppb O₃ treatment were observed only occasionally.

Analyzing single HR-like injuries, the shape (data not shown) and size of lesions remained stable over time or with increasing POD₀. Furthermore, they did not respond to the O₃ treatment, leaf position, or interaction factors (Figures 6A,B; all factors: ns). The only change observed was increasing oxidation with longer exposures and at higher POD₀ (Figures 6C,D;

Time, POD: *p* < 0.05). Hence, the observed increases in the percentage area of HR-like lesions with time, or higher POD₀ and in response to the O₃ treatment (Figures 5A,B) resulted as a consequence of the multiplication of single HR-like reactions and higher lesion density, rather than from increased growth of already developed injuries. However, both the higher lesion density and the growing oxidation of individual HR-like lesions could contribute to the observed increase in the leaf percentage area showing oxidation (Figures 5C,D).

## Emergence of Visible Symptoms

The first visible symptoms in leaves exposed to the O₃ treatments appeared by the end of the experiment, that is, 23 days after the start of exposure. The first HR-like lesions had been detected more than 2 weeks earlier, whereas the observed drops in foliar chlorophyll content index were rather synchronous (Figure 7A vs. Figures 3A, 5A). These visible symptoms consisted of intercostal necrotic dark spots spread in the leaf

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

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

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

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

<table><thead><tr><th></th><th>Total HR-like</th><th>Oxidized HR-like</th><th></th><th>Total HR-like</th><th>Oxidized HR-like</th></tr></thead><tbody><tr><td>O₃ treatment</td><td>*</td><td>ns</td><td>O₃ treatment</td><td>ns</td><td>ns</td></tr><tr><td>Leaf Position</td><td>ns</td><td>ns</td><td>Leaf Position</td><td>ns</td><td>ns</td></tr><tr><td>Time</td><td>***</td><td>***</td><td>POD₀</td><td>***</td><td>***</td></tr><tr><td>O₃ treatment:Time</td><td>***</td><td>***</td><td>O₃ treatment:POD₀</td><td>**</td><td>***</td></tr><tr><td>O₃ treatment:Leaf Position</td><td>**</td><td>ns</td><td>O₃ treatment:Leaf Position</td><td>ns</td><td>ns</td></tr><tr><td>Leaf Position:Time</td><td>*</td><td>ns</td><td>Leaf Position:POD₀</td><td>***</td><td>***</td></tr><tr><td>O₃ treatment:Time:Leaf Position</td><td>*</td><td>ns</td><td>O₃ treatment:POD₀:Leaf Position</td><td>*</td><td>ns</td></tr></tbody></table>

**FIGURE 5** | Development dynamics of HR-like lesions (A,B) and lesion oxidation (C,D) in hybrid poplar leaves (Populus tremula x alba), as a function of the assessment time (A,C), phytotoxic O₃ dose (POD₀; B,D), O₃ treatment (charcoal-filtered air ●, 80 ppb O₃ ▲, 100 ppb O₃ ■), leaf position and interactions (model: lmer[sqrt variable]) ~ O₃ treatment * leaf position * time + (1 | tree / chamber); ***p ≤ 0.001; **p ≤ 0.01; *p ≤ 0.05; ns, not significantly different). Values represent means ± SE of percentage area of leaf discs showing microscopic injury, n = 4. Different letters indicate significant differences between treatments at a given assessment date (post-hoc Tukey's honest significant difference, *p* ≤ 0.05).
<figure><img src="image_11.png" /></figure>

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

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

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

<table><thead><tr><th></th><th>Lesion size</th><th>Lesion oxidation</th><th></th><th>Lesion size</th><th>Lesion oxidation</th></tr></thead><tbody><tr><td>O<sub>3</sub> treatment</td><td>ns</td><td>ns</td><td>O<sub>3</sub> treatment</td><td>ns</td><td>ns</td></tr><tr><td>Leaf Position</td><td>ns</td><td>ns</td><td>Leaf Position</td><td>ns</td><td>ns</td></tr><tr><td>Time</td><td>ns</td><td>*</td><td>POD<sub>0</sub></td><td>ns</td><td>*</td></tr><tr><td>O<sub>3</sub> treatment:Time</td><td>ns</td><td>ns</td><td>O<sub>3</sub> treatment:POD<sub>0</sub></td><td>ns</td><td>ns</td></tr><tr><td>O<sub>3</sub> treatment:Leaf Position</td><td>ns</td><td>ns</td><td>O<sub>3</sub> treatment:Leaf Position</td><td>ns</td><td>ns</td></tr><tr><td>Leaf Position:Time</td><td>ns</td><td>ns</td><td>Leaf Position:POD<sub>0</sub></td><td>ns</td><td>ns</td></tr><tr><td>O<sub>3</sub> treatment:Time:Leaf Position</td><td>ns</td><td>ns</td><td>O<sub>3</sub> treatment:POD<sub>0</sub>:Leaf Position</td><td>ns</td><td>ns</td></tr></tbody></table>

**FIGURE 6** | Dynamics of changes in the size (A,B) and degree of oxidation (C,D) of single HR-like lesions in hybrid poplar leaves (Populus tremula x alba), as a function of the assessment time (A,C), phytotoxic O₃ dose (POD₀; B,D), O₃ treatment (charcoal-filtered air ●, 80 ppb O₃ ▲, 100 ppb O₃ ■), leaf position and interactions [model: lmer(variable) ~ O₃ treatment * leaf position * time + (1 | tree/chamber); * p ≤ 0.05; ns, not significantly different, na not tested]. (A,B) The size of HR-like lesions did not respond to the experimental factors and interactions. Values represent means ± SE of distinct lesion size (A,B) and percentage area showing oxidation (C,D), n = 4.

blade, as previously observed in poplar (Figure 7A; Cabane et al., 2004; Giacomo et al., 2010; Dghim et al., 2013). By the end of exposure, only low levels of injury could develop (O₃ treatment: ns; O₃ treatment*Time *p* < 0.001), with significantly higher percentages in the 100 vs. 80 ppb O₃ treatment at the third leaf position only. In younger leaves, the visible injury appeared 4 days later than at the third leaf position, and differences between treatments throughout the experiment remained non-significant (leaf position, O₃ treatment*time*leaf position: *p* < 0.05; O₃ treatment*leaf position: *p* < 0.01).

When expressed as a function of POD₀, the visible injuries appeared at a lower O₃ threshold at the third than tenth leaf position (9 and 12 mmol O₃ m⁻² in the case of 100 ppb treatment; Figure 7B). With only the beginning of injury development assessed in a 30-day experiment, the O₃ treatment effects could not reach any significance, and only preliminary information on the dynamics of visible symptom development

was thus obtained. The simultaneous detection of early injury in the two O₃ treatments (Figure 7A) thereby resulted in higher injury values for similar POD₀ in the 80 vs. 100 ppb O₃ treatment (Figure 7B; O₃ treatment*POD: *p* < 0.05). Similarly, differences between younger and older leaves were detected as a trend only, with still very low injury values recorded at the tenth leaf position and for the higher POD₀ values only (leaf position: ns; leaf position*POD: *p* < 0.001).

DISCUSSION

Dynamics of Physiological and Structural Responses to O₃ Stress

The physiological and structural responses detected during 30 days of O₃ exposure developed mostly before and in some cases at the same time as the initial visible injury and first
<figure><img src="image_15.png" /></figure>

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

<table><thead><tr><th></th><th>Leaf injury</th><th></th><th>Leaf injury</th></tr></thead><tbody><tr><td>O<sub>3</sub> treatment</td><td>ns</td><td>O<sub>3</sub> treatment</td><td>ns</td></tr><tr><td>Leaf Position</td><td>*</td><td>Leaf Position</td><td>ns</td></tr><tr><td>Time</td><td>***</td><td>POD<sub>0</sub></td><td>***</td></tr><tr><td>O<sub>3</sub> treatment:Time</td><td>***</td><td>O<sub>3</sub> treatment:POD<sub>0</sub></td><td>*</td></tr><tr><td>O<sub>3</sub> treatment:Leaf Position</td><td>**</td><td>O<sub>3</sub> treatment:Leaf Position</td><td>ns</td></tr><tr><td>Leaf Position:Time</td><td>**</td><td>Leaf Position:POD<sub>0</sub></td><td>***</td></tr><tr><td>O<sub>3</sub> treatment:Time:Leaf Position</td><td>*</td><td>O<sub>3</sub> treatment:POD<sub>0</sub>:Leaf Position</td><td>*</td></tr></tbody></table>

**FIGURE 7** | Development dynamics of visible O₃ injury in hybrid poplar leaves (*Populus tremula* x *alba*), as a function of the assessment time (**A**), phytotoxic O₃ dose (POD₀; **B**), O₃ treatment (charcoal-filtered air ●, 80 ppb O₃ ▲, 100 ppb O₃ ■), leaf position and interactions (lmer[log(variable)] ~ O₃ treatment * leaf position * time + (1 | tree / chamber); ***p ≤ 0.001; **p ≤ 0.01; *p ≤ 0.05; ns, not significantly different). The insert in (**A**) shows typical O₃ symptoms, in the form of intercostal necrotic dark spots, within a leaf at third leaf position exposed for 30 days (POD₀ = 9.3 mmol O₃ m⁻², 100 ppb O₃ treatment) vs. an asymptomatic leaf from the CF treatment, sampled at the same leaf position and the end of treatment. Values represent means ± SE of leaf percentage area showing visible injury, *n* = 6. Letters indicate significant differences between treatments at a given assessment date (*post-hoc* Tukey's Honest significant difference, *p* ≤ 0.05).

leaf shedding, whereas no O₃ effect on the gross morphology of trees was observed. The two O₃ treatments accelerated the ontological decline of leaf gas exchange at the two-leaf positions, once leaf physiology had reached maturity. O₃-induced reductions in the stomatal conductance and net CO₂ assimilation are well-documented in the case of various species, including poplars (Pell et al., 1992; Bagard et al., 2015; Dusart et al., 2019a) but more rarely with a leaf ontogenetic perspective. The decrease in net CO₂ assimilation could result from a smaller stomatal aperture, limiting CO₂ availability. Indeed, the O₃ effects on stomata are well-established (Kangasjarvi et al., 2005), including slower movements of cell guard cells upon exposure referred to as stomatal sluggishness, as observed in different species but with higher measurement frequency than in our case (Paoletti and Grulke, 2010; Dumont et al., 2013; Dusart et al., 2019b). Lower mesophyll conductance could also contribute to the observed acceleration of Aₘₐt reduction, as compared with the ontogenetic decrease observed in CF trees (Xu et al., 2019). Although not significant in older foliage, the drop in Aₘₐt tended to be more expressed in the 100 vs. 80 ppb treatment after 20 days, suggesting an additional reduction in carboxylation efficiency and rate of electron transport (Bagard et al., 2008; Shang et al., 2017). A possible cause could be the starting degradation of the photosynthetic machinery, as suggested by the concomitant reduction in the chlorophyll content index.

The reduction in leaf chlorophyll content is another well-documented leaf physiological response to elevated O₃ (Reich, 1983; Bagard et al., 2015; Dusart et al., 2019b), indicative of

ACS, together with other markers of chloroplast degeneration (Mikkelsen and HeideJorgensen, 1996; Günthardt-Goerg and Vollenweider, 2007; Moura et al., 2018). Also typical of ACS and degenerative processes was the progressive and mostly monotonic reductions observed in the case of all physiology parameters (g<sub>w</sub>, Aₘₐt, and chlorophyll index). This contrasted with abrupt responses upon exceedance of a threshold, as observed in the case of HR-like reactions. In the field, the ontogeny-driven ACS latency observed in both younger and older foliage (total chlorophyll index) may also contribute to delaying the onset of degenerative events in response to O₃-stress, as suggested by the development of ACS traits and visible injury rarely occurring before summer, once the foliage has fully matured (Vollenweider et al., 2019).

The new cell necrosis assay, using computer-assisted color image analysis, allowed us to monitor the emergence and development of validated HR-like reactions for the first time. It provided unprecedented capacity for quantitative assessments of cell death reactions in experimental conditions, overcoming limitations and uncertainties regarding the observation of visible injury only. Typical HR-like markers were detected in the lesions (Paoletti et al., 2009; Vollenweider et al., 2013, 2019; Feng et al., 2016). However, there were marked differences in classical traits as well, including the mid- rather than the upper-mesophyll location of HR-like lesions or a missing intra- and intercellular gradient of injury. Such features indicated low levels of photo-oxidative stress (Foyer et al., 1994; Günthardt-Goerg and Vollenweider, 2007;
Guerrero et al., 2013). Given maximum PAR above 2000 vs. 350 µmol m⁻² s⁻¹ with high (Ritchie, 2010; Poorter et al., 2019) vs. low light conditions, HR-like reaction peculiarities – together with the late onset of ACS – can be attributed to specifics in the environmental conditions, especially regarding PAR supply. This finding thus provides further confirmation of the close dependency relating the O₃ symptom expression in foliage and precise experimental and exposure conditions of tested material (Paoletti et al., 2009; Moura et al., 2018; Vollenweider et al., 2019).

The rather sigmoidal-like injury dynamics observed in older foliage was in good agreement with already existing molecular and trait evidence on HR-like processes. The 6 days/3 mmol O₃ m⁻² s⁻¹ delay between the start of exposure and occurrence of first lesions was thus indicative of the O₃ dose-dependent onset of genetically controlled PCD (Rao et al., 2000; Overmyer et al., 2005). The steep increase in injury percentage area reflected the rapid cell death completion once PCD started (Overmyer et al., 2005; Günthardt-Goerg and Vollenweider, 2007). Finally, the plateau reached was indicative of lesion containment, blocking its further spread (Overmyer et al., 2003; Kangasjarvi et al., 2005; Marchica et al., 2019). In younger foliage, the experiment was terminated before a plateau could be reached, with values of lesion percentage area which would probably have been sizably higher than at the third leaf position. Further suggesting the genetic control of PCD, in older foliage, the first HR-like reactions occurred 17 days/at 5.3 mmol O₃ m⁻² s⁻¹ before any evidence of biochemical limitation and chloroplast injury, as indicated by low levels of chlorophyll content indexes. The early HR-like responses, and their antecedence concerning ACS and first visible injury, contrasted with field evidence (Vollenweider et al., 2019), further outlining how important the environmental conditions can be regarding response order and dynamics.

Image analysis in WinCELL based on three color classes allowed us to quantify the structurally contrasted non-oxidized and oxidized lesions, based on constitutive and stained-color characteristics. REDOX changes during oxidative stress and cell death form an important cell physiology process (Foyer and Noctor, 2005) that is well-documented in the case of O₃ stress (Ranieri et al., 2000; Baier et al., 2005; Chen and Gallie, 2005; Bellini and De Tullio, 2019), detectable with different structural and ultrastructural markers (Moura et al., 2018; Vollenweider et al., 2019) and underlying changes in visible symptom expression. Expressed as leaf area percentages, the oxidized and total HR-like lesions showed similar dynamics and responses to increasing POD₀ (same results for non-oxidized lesions, data not shown). The main differences in oxidized vs. non-oxidized lesions included their (1) lower percentage area, (2) delay in development, and (3) higher severity (i.e., cell wall breaks, cell content disruption). Non-oxidized and oxidized lesions may thus correspond to two types or two stages of HR-like reactions. However, the first hypothesis appears unlikely, given lacking molecular evidence for alleged PCD severity variation. In favor of the second, oxidized vs. non-oxidized lesions appeared later, and the oxidation degree of a lesion increased with time. However, it implies the further evolution of

HR-like lesions after cell collapse and death, which needs further structural and ultrastructural confirmation.

The visible injury was detected only once 2–5% of the leaf percentage area showed oxidized lesions, thus with a detection delay and a POD₀ gap compared to the onset of HR-like reactions amounting to 18 days and 4.9 mmol O₃ m⁻² s⁻¹. Similarly, risk assessment studies using visible injury markers rather target the late and final structural evolution of responses to O₃ stress in foliage, with possible interspecific variation, instead of the injury appearance in foliage.

## Leaf Position Dependency of Responses to O₃ Stress

Reductions in leaf gas exchanges or the development of HR-like lesions and visible symptoms at the tenth vs. third leaf position occurred later and for larger POD₀. Given the higher stomatal conductance and POD₀ in younger leaves, their greater physiological activity and lower levels of injury suggest higher O₃ tolerance, while a contribution by enhanced stomatal closure can be excluded. This finding is confirmed by similar reports on enhanced O₃ tolerance in maturing leaves (Reich, 1983; Paakkonen et al., 1996; Strohm et al., 2002; Bagard et al., 2008; Zhang et al., 2010; Guerrero et al., 2013). This may be related to sink functional properties and larger resource availability for defense and repair (Coleman, 1986). Resource availability in young leaves could be increased by the supply of nutrients (such as nitrogen, potassium, and phosphorus) coming from senescent leaves (Maillard et al., 2015; Have et al., 2017). Hence, the concentration of phenolics with antioxidant properties and other antioxidative capacities decline during the sink-to-source transition in maturing foliage, thus increasing leaf susceptibility to oxidative stress (Coleman, 1986; Strohm et al., 2002; Blokhina and Fagerstedt, 2006; Bellini and De Tullio, 2019). However, and in contrast to older foliage, the development of HR-like lesions as a function of POD₀ at the tenth leaf position depended on the O₃ treatment, with a higher O₃ tolerance in the 80 ppb O₃ treatment. Given the high O₃ dose in maturing leaves, this finding highlights the importance of the O₃ absorption rate given the saturation of the antioxidative system. The higher O₃ tolerance in younger vs. older foliage was further confirmed by their still comparable leaf percentage areas showing HR-like lesions in the 100 ppb O₃ treatment despite 1.8 times higher POD₀ at the tenth leaf position.

## Reaction Gradient in Foliage Concerning Critical O₃ Levels

In our experiment, the current CL (POD<sub>Y SPEC</sub> for beech and birch = 5.2 mmol O₃ m⁻²; Mills et al., 2017) was equivalent to a POD₀ of 5.7 mmol O₃ m⁻². By the end of exposure, this CL had thus been exceeded by 1.54 and 2.35 times at the third and tenth leaf positions, respectively. If any impairment of tree morphology and biomass was still lacking, reductions in leaf gas exchange, development of structural injury, and the emergence of visible symptoms at the third leaf position had already been observed for O₃ dose, amounting to 0.82, 0.69, and 1.46 times the current CL, respectively. At the tenth
leaf position, these responses were detected for POD~V SPEC~ 1.83, 1.18, and 2.32 times above CL. These findings highlight the high dependency of sensitivity evaluations on the selected parameters and scale of observation. They also outline the within-tree gradient of sensitivity to O~3~ stress, given the large size of such organisms, which, as a result, complicates O~3~ risk assessment. They finally indicate that below CL, significant effects in the foliage of trees, such as in the impairment of leaf physiology and development of microscopic necrosis in extended parts of mesophyll, can be expected. These responses may already contribute to reduced carbon uptake and storage in foliage and other tree organs before reaching CL thresholds.

## CONCLUSION

In this study, we characterized the dynamics of physiological, structural, and morphological responses to two levels of O~3~ exposure and as a function of time, POD~0~ and leaf position, in fully controlled conditions. We observed contrasting dynamics, monotonic or sigmoidal-like, as a function of plant responses but irrespective of leaf position, before any visible symptoms and effects on the gross morphology of trees. The first microscopic necrosis developed weeks before the appearance of visible symptoms and at half the O~3~ dose. Concerning experimental hypotheses (H), the sequential development and distinct dynamics of physiological, structural, and morphological responses to O~3~ stress was confirmed (confirmation of H1); both HR-like and ACS responses were elicited, the former occurring first (confirmation of H2a, rejection of H2b). When expressed as a function of POD~0~, leaf responses did not depend on the O~3~ treatment (confirmation of H3a), except for the development of structural injury that depended on the O~3~ absorption rates in younger foliage (partial rejection of H3b). Finally, response dynamics were strongly related to leaf age as a function of time or POD~0~, showing delay in younger foliage (confirmation of H4). This study thus sheds light on the syndrome of early reactions to O~3~ stress and disentangles the specific dynamics of distinct but co-occurring plant responses, before CL exceedance. The resulting variety of symptoms, as observed by the end of the experiment, provides an exemplary experimental demonstration for integrative injury display, as found in the field late in summer. Given ACS

and HR-like timing inversion, compared to field conditions, the ontogenetical and environmental drivers also appear to have a prevailing effect over the sensitivity of affected markers, regarding the timing and dynamics of each cellular response. Whatever the exact sequence order of early reactions to O~3~ stress below CL – in the field or controlled conditions, they will modify key foliage properties. Such impacts can be relevant for some, e.g., biodiversity ecosystem services before those of economic significance – i.e., wood production, as targeted by flux-based CL.

## DATA AVAILABILITY STATEMENT

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

## AUTHOR CONTRIBUTIONS

BT, YJ, MC, and PV: conception or design of the work and final approval of the version to be published. BT: data collection. BT, PV, AG, DT, YJ, and MC: data analysis and interpretation. BT and PV: drafting the article. BT, PV, AG, DT, YJ, MC, and MS: critical revision of the article. All authors contributed to the article and approved the submitted version.

## FUNDING

The research was funded by grants from the French National Research Agency (ANR, "Investissements d'Avenir" from the program Lab of Excellence ARBRE: ANR-11-LABX-0002-01), Grand-Est region, France, and WSL (201701N1428).

## ACKNOWLEDGMENTS

Technical support by Terry Menard (microscopy), Christophe Robin (leaf physiology), and Jean-Charles Olry and Stéphane Martin (members of the experimental phytotrons platform of Lorraine PEPLor, University of Lorraine, France) was gratefully acknowledged.

## REFERENCES

Bagard, M., Jolivet, Y., Hasenfratz-Sauder, M. P., Gerard, J., Dizengremel, P., and Le Thiec, D. (2015). Ozone exposure and flux-based response functions for photosynthetic traits in wheat, maize and poplar. *Environ. Pollut.* 206, 411–420. doi: 10.1016/j.envpol.2015.07.046

Bagard, M., Le Thiec, D., Delacote, E., Hasenfratz-Sauder, M. P., Banvoy, J., Gerard, J., et al. (2008). Ozone-induced changes in photosynthesis and photorespiration of hybrid poplar in relation to the developmental stage of the leaves. *Physiol. Plant.* 134, 559–574. doi: 10.1111/j.1399-3054.2008.01160.x

Baier, M., Kandlbinder, A., Golldack, D., and Dietz, K. J. (2005). Oxidative stress and ozone: perception, signaling and response. *Plant Cell Environ.* 28, 1012–1020. doi: 10.1111/j.1365-3040.2005.01326.x

Bates, D., Machler, M., Bolker, B. M., and Walker, S. C. (2015). Fitting linear mixed-effects models using lme4. *J. Stat. Softw.* 67, 1–48. doi: 10.18637/jss.v067.i01

Bellini, E., and De Tullio, M. C. (2019). Ascorbic acid and ozone: novel perspectives to explain an elusive relationship. *Plants* 8:122. doi: 10.3390/plants8050122

Bhattacharjee, S. (2005). Reactive oxygen species and oxidative burst: roles in stress, senescence and signal transduction in plants. *Curr. Sci.* 89, 1113–1121.

Blokhina, O., and Fagerstedt, K. (2006). *Oxidative Stress and Antioxidant Defenses in Plants*. Covent Garden: Imperial College Press.

Büker, P., Feng, Z., Uddling, J., Briolat, A., Alonso, R., Braun, S., et al. (2015). New flux based dose-response relationships for ozone for European forest tree species. *Environ. Pollut.* 206, 163–174. doi: 10.1016/j.envpol.2015.06.033

Cabane, M., Pireaux, J. C., Leger, E., Weber, E., Dizengremel, P., Pollet, B., et al. (2004). Condensed lignins are synthesized in poplar leaves exposed to ozone. *Plant Physiol.* 134, 586–594. doi: 10.1104/pp.103.031765

Chen, Z., and Gallie, D. R. (2005). Increasing tolerance to ozone by elevating foliar ascorbic acid confers greater protection against ozone than increasing avoidance. *Plant Physiol.* 138, 1673–1689. doi: 10.1104/pp.105.062000
Coleman, J. S. (1986). Leaf development and leaf stress: increased susceptibility associated with sink-source transition. *Tree Physiol.* 2, 289-299. doi: 10.1093/treephys/2.1-2-3.289

Dghim, A. A., Dumont, J., Hasenfratz-Sauder, M. P., Dizengremel, P., Le Thiec, D., and Jolivet, Y. (2013). Capacity for NADPH regeneration in the leaves of two poplar genotypes differing in ozone sensitivity. *Physiol. Plant.* 148, 36-50. doi: 10.1111/j.1399-3054.2012.01686.x

Dizengremel, P., Jolivet, Y., Tuzet, A., Ranieri, A., and Le Thiec, D. (2013). "Integrative leaf-level phytotoxic ozone dose assessment for forest risk modelling," in *Climate Change, Air Pollution and Global Challenges: Understanding and Perspectives from Forest Research*. eds. R. Matyssek, N. Clarke, P. Cudlin, T. N. Mikkelsen, J. P. Tuovinen, G. Wieser et al (Netherlands: Elsevier), 267-288.

Dumont, J., Keski-Saari, S., Keinanen, M., Cohen, D., Ningre, N., Kontunen-Soppela, S., et al. (2014). Ozone affects ascorbate and glutathione biosynthesis as well as amino acid contents in three Euramerican poplar genotypes. *Tree Physiol.* 34, 253-266. doi: 10.1093/treephys/tpu004

Dumont, J., Spicher, F., Montpied, P., Dizengremel, P., Jolivet, Y., and Le Thiec, D. (2013). Effects of ozone on stomatal responses to environmental parameters (blue light, red light, CO2 and vapour pressure deficit) in three *Populus deltoides* x *Populus nigra* genotypes. *Environ. Pollut.* 173, 85-96. doi: 10.1016/j.envpol.2012.09.026

Dusart, N., Gerard, J., Le Thiec, D., Collignon, C., Jolivet, Y., and Vaultier, M. N. (2019a). Integrated analysis of the detoxification responses of two Euramerican poplar genotypes exposed to ozone and water deficit: focus on the ascorbate-glutathione cycle. *Sci. Total Environ.* 651, 2365-2379. doi: 10.1016/j.scitotenv.2018.09.367

Dusart, N., Vaultier, M. N., Olry, J. C., Bure, C., Gerard, J., Jolivet, Y., et al. (2019b). Altered stomatal dynamics of two Euramerican poplar genotypes submitted to successive ozone exposure and water deficit. *Environ. Pollut.* 252, 1687-1697. doi: 10.1016/j.envpol.2019.06.110

Faoro, F., and Iriti, M. (2009). Plant cell death and cellular alterations induced by ozone: key studies in Mediterranean conditions. *Environ. Pollut.* 157, 1470-1477. doi: 10.1016/j.envpol.2008.09.026

Feng, G., Calatayud, V., García-Breijo, F., Reig-Armiñana, J., and Feng, Z. (2016). Effects of elevated ozone on physiological, anatomical and ultrastructural characteristics of four common urban tree species in China. *Ecol. Indic.* 67, 367-379. doi: 10.1016/j.ecolind.2016.03.012

Foyer, C. H., Lelandais, M., and Kunert, K. J. (1994). Photooxidative stress in plants. *Physiol. Plant.* 92, 696-717. doi: 10.1111/j.1399-3054.1994.tb03042.x

Foyer, C. H., and Noctor, G. (2005). Redox homeostasis and antioxidant signaling: a metabolic interface between stress perception and physiological responses. *Plant Cell* 17, 1866-1875. doi: 10.1105/tpc.105.033589

Fu, T. M., and Tian, H. (2019). Climate change penalty to ozone air quality: review of current understandings and knowledge gaps. *Curr. Pollut. Rep.* 5, 159-171. doi: 10.1007/s40726-019-00115-6

Fuhrer, J., Skarby, L., and Ashmore, M. R. (1997). Critical levels for ozone effects on vegetation in Europe. *Environ. Pollut.* 97, 91-106. doi: 10.1016/S0269-7491(97)00067-5

Giacomo, B., Forino, L. M. C., Tagliasacchi, A. M., Bernardi, R., and Durante, M. (2010). Ozone damage and tolerance in leaves of two poplar genotypes. *Caryologia* 63, 422-434. doi: 10.1080/00087114.2010.10589755

Guerrero, C. C., Gunthardt-Goerg, M. S., and Vollenweider, P. (2013). Foliar symptoms triggered by ozone stress in irrigated holm oaks from the city of Madrid, Spain. *PLoS One* 8:e69171. doi: 10.1371/journal.pone.0069171

Günthardt-Goerg, M. S., and Vollenweider, P. (2007). Linking stress with macroscopic and microscopic leaf response in trees: new diagnostic perspectives. *Environ. Pollut.* 147, 467-488. doi: 10.1016/j.envpol.2006.08.033

Have, M., Marmagne, A., Chardon, F., and Masclaux-Daubresse, C. (2017). Nitrogen remobilization during leaf senescence: lessons from *Arabidopsis* to crops. *J. Exp. Bot.* 68, 2513-2529. doi: 10.1093/jxb/erw365

Jolivet, Y., Bagard, M., Cabane, M., Vaultier, M. N., Gandin, A., Afif, D., et al. (2016). Deciphering the ozone-induced changes in cellular processes: a prerequisite for ozone risk assessment at the tree and forest levels. *Ann. For. Sci.* 73, 923-943. doi: 10.1007/s13595-016-0580-3

Joo, J. H., Wang, S. Y., Chen, J. G., Jones, A. M., and Fedoroff, N. V. (2005). Different signaling and cell death roles of heterotrimeric G protein alpha and beta subunits in the *Arabidopsis* oxidative stress response to ozone. *Plant Cell* 17, 957-970. doi: 10.1105/tpc.104.029603

Kangasjarvi, J., Jaspers, P., and Kollist, H. (2005). Signaling and cell death in ozone-exposed plants. *Plant Cell Environ.* 28, 1021-1036. doi: 10.1111/j.1365-3040.2005.01325.x

Karlsson, P., Braun, S., Broadmeadow, M., Elvira, S., Emberson, L., Gimeno, B. S., et al. (2007). Risk assessments for forest trees: the performance of the ozone flux versus the AOT concepts. *Environ. Pollut.* 146, 608-616. doi: 10.1016/j.envpol.2006.06.012

Karlsson, P. E., Klingberg, J., Engardt, M., Andersson, C., Langner, J., Karlsson, G. P., et al. (2017). Past, present and future concentrations of ground-level ozone and potential impacts on ecosystems and human health in northern Europe. *Sci. Total Environ.* 576, 22-35. doi: 10.1016/j.scitotenv.2016.10.061

Lamaud, E., Loubet, B., Irvine, M., Stella, P., Personne, E., and Cellier, P. (2009). Partitioning of ozone deposition over a developed maize crop between stomatal and non-stomatal uptakes, using eddy-covariance flux measurements and modelling. *Agric. For. Meteorol.* 149, 1385-1396. doi: 10.1016/j.agrformet.2009.03.017

Lenth, R. V. (2016). Least-squares means: the R package lsmeans. *J. Stat. Softw.* 69, 1-33. doi: 10.18637/jss.v069.i01

Li, P., Feng, Z. Z., Catalayud, V., Yuan, X. Y., Xu, Y. S., and Paoletti, E. (2017). A meta-analysis on growth, physiological and biochemical responses of woody species to ground-level ozone highlights the role of plant functional types. *Plant Cell Environ.* 40, 2369-2380. doi: 10.1111/pce.13043

Maas, R., and Grennfelt, P. (2016). Towards cleaner air. Scientific assessment report. Emap steering body and working group on effects of the convention on long-range transboundary air pollution, Oslo. 12-16.

Maillard, A., Diquelou, S., Billard, V., Laine, P., Garnica, M., Prudent, M., et al. (2015). Leaf mineral nutrient remobilization during leaf senescence and modulation by nutrient deficiency. *Front. Plant Sci.* 6:317. doi: 10.3389/fpls.2015.00317

Marchica, A., Lorenzini, G., Papini, R., Bernardi, R., Nali, C., and Pellegrini, E. (2019). Signaling molecules responsive to ozone-induced oxidative stress in *Salvia officinalis*. *Sci. Total Environ.* 657, 568-576. doi: 10.1016/j.scitotenv.2018.11.472

Massman, W. J. (1998). A review of the molecular diffusivities of H₂O, CO₂, CH₄, CO, O₃, SO₂, NH₃, N₂O, NO, and NO₂ in air, O₂ and N₂ near STP. *Atmos. Environ.* 32, 1111-1127. doi: 10.1016/S1352-2310(97)00391-9

Mikkelsen, T. N., and Heidejorgensen, H. S. (1996). Acceleration of leaf senescence in *Fagus sylvatica* L by low levels of tropospheric ozone demonstrated by leaf colour, chlorophyll fluorescence and chloroplast ultrastructure. *Trees-Struct. Funct.* 10, 145-156. doi: 10.1007/BF02340766

Mills, G., Harmens, H., Hayes, F., Pleijel, H., Buker, P., and González, I. (2017). "Mapping critical levels for vegetation," in *Manual on Methodologies and Criteria for Modelling and Mapping Critical Loads and Levels and Air Pollution Effects. Risks and Trends*. International Cooperative Programme on Effects of Air Pollution on Natural Vegetation and Crops, 1-66.

Mills, G., Pleijel, H., Braun, S., Buker, P., Bermejo, V., Calvo, E., et al. (2011). New stomatal flux-based critical levels for ozone effects on vegetation. *Atmos. Environ.* 45, 5064-5068. doi: 10.1016/j.atmosenv.2011.06.009

Moura, B. B., Alves, E. S., Marabesi, M. A., De Souza, S. R., Schaub, M., and Vollenweider, P. (2018). Ozone affects leaf physiology and causes injury to foliage of native tree species from the tropical Atlantic Forest of southern Brazil. *Sci. Total Environ.* 610, 912-925. doi: 10.1016/j.scitotenv.2017.08.130

Musselman, R. C., Lefohn, A. S., Massman, W. J., and Heath, R. L. (2006). A critical review and analysis of the use of exposure- and flux-based ozone indices for predicting vegetation effects. *Atmos. Environ.* 40, 1869-1888. doi: 10.1016/j.atmosenv.2005.10.064

Overmyer, K., Brosche, M., and Kangasjarvi, J. (2003). Reactive oxygen species and hormonal control of cell death. *Trends Plant Sci.* 8, 335-342. doi: 10.1016/S1360-1385(03)00135-3

Overmyer, K., Brosche, M., Pellinen, R., Kuittinen, T., Tuominen, H., Ahlfors, R., et al. (2005). Ozone-induced programmed cell death in the *Arabidopsis* radical-induced cell death1 mutant. *Plant Physiol.* 137, 1092-1104. doi: 10.1104/pp.104.055681

Paakkonen, E., Metsarinne, S., Holopainen, T., and Karenlampi, L. (1996). The ozone sensitivity of birch (*Betula pendula*) in relation to the developmental stage of leaves. *New Phytol.* 132, 145-154. doi: 10.1111/j.1469-8137.1996.tb04520.x

Paoletti, E., Contran, N., Bernasconi, P., Gunthardt-Goerg, M. S., and Vollenweider, P. (2009). Structural and physiological responses to ozone in
manna ash (Fraxinus ornus L.) leaves of seedlings and mature trees under controlled and ambient conditions. *Sci. Total Environ.* 407, 1631–1643. doi: 10.1016/j.scitotenv.2008.11.061

Paoletti, E., and Grulke, N. E. (2010). Ozone exposure and stomatal sluggishness in different plant physiognomic classes. *Environ. Pollut.* 158, 2664–2671. doi: 10.1016/j.envpol.2010.04.024

Pasqualini, S., Piccioni, C., Reale, L., Ederli, L., Della Torre, G., and Ferranti, F. (2003). Ozone-induced cell death in tobacco cultivar Bel W3 plants. The role of programmed cell death in lesion formation. *Plant Physiol.* 133, 1122–1134. doi: 10.1104/pp.103.026591

Pell, E. J., Eckardt, N., and Enyedi, A. J. (1992). Timing of ozone stress and resulting status of ribulose bisphosphate carboxylase oxygenase and associated net photosynthesis. *New Phytol.* 120, 397–405. doi: 10.1111/j.1469-8137.1992.tb01080.x

Pell, E. J., Sinn, J. P., Brendley, B. W., Samuelson, L., Vinten-Johansen, C., Tien, M., et al. (1999). Differential response of four tree species to ozone-induced acceleration of foliar senescence. *Plant Cell Environ.* 22, 779–790. doi: 10.1046/j.1365-3040.1999.00449.x

Poorter, H., Niinemets, U., Ntagkas, N., Siebenkas, A., Maenpaa, M., Matsubara, S., et al. (2019). A meta-analysis of plant responses to light intensity for 70 traits ranging from molecules to whole plant performance. *New Phytol.* 223, 1073–1105. doi: 10.1111/nph.15754

Proietti, C., Anav, A., De Marco, A., Sicard, P., and Vitale, M. (2016). A multi-sites analysis on the ozone effects on gross primary production of european forests. *Sci. Total Environ.* 556, 1–11. doi: 10.1016/j.scitotenv.2016.02.187

Ranieri, A., Petacco, F., Castagna, A., and Soldatini, G. F. (2000). Redox state and peroxidase system in sunflower plants exposed to ozone. *Plant Sci.* 159, 159–167. doi: 10.1016/S0168-9452(00)00352-6

Rao, M. V., Lee, H., Creelman, R. A., Mullet, J. E., and Davis, K. R. (2000). Jasmonic acid signaling modulates ozone-induced hypersensitive cell death. *Plant Cell* 12, 1633–1646. doi: 10.1105/tpc.12.9.1633

R Development Core Team (2017). "R: A language and environment for statistical computing." Vienna, Austria: R Foundation for Statistical Computing.

Reich, P. B. (1983). Effects of low concentrations of O3 on net photosynthesis, dark respiration, and chlorophyll contents in aging hybrid poplar leaves. *Plant Physiol.* 73, 291–296. doi: 10.1104/pp.73.2.291

Revell, L. E., Tummon, F., Stenke, A., Sukhodolov, T., Coulon, A., Rozanov, E., et al. (2015). Drivers of the tropospheric ozone budget throughout the 21st century under the medium-high climate scenario RCP 6.0. *Atmospheric Chem. Phys.* 15, 5887–5902. doi: 10.5194/acp-15-5887-2015

Ritchie, R. J. (2010). Modelling photosynthetic photon flux density and maximum potential gross photosynthesis. *Photosynthetica* 48, 596–609. doi: 10.1007/s11099-010-0077-5

Sanz, M., and Catalayud, V. (2011). *Ozone Injury in European Forest Species*. Available at: http://www.ozoneinjury.org/ (Accessed December 12, 2020).

Schindelin, J., Arganda-Carreras, I., Frise, E., Kaynig, V., Longair, M., Pietzsch, T., et al. (2012). Fiji: an open-source platform for biological-image analysis. *Nat. Methods* 9, 676–682. doi: 10.1038/nmeth.2019

Schraudner, M., Moeder, W., Wiese, C., Van Camp, W., Inze, D., Langebartels, C., et al. (1998). Ozone-induced oxidative burst in the ozone biomonitor plant, tobacco Bel W3. *Plant J.* 16, 235–245. doi: 10.1046/j.1365-313x.1998.00294.x

Shang, B., Feng, Z. Z., Li, P., Yuan, X. Y., Xu, Y. S., and Calatayud, V. (2017). Ozone exposure- and flux-based response relationships with photosynthesis, leaf morphology and biomass in two poplar clones. *Sci. Total Environ.* 603, 185–195. doi: 10.1016/j.scitotenv.2017.06.083

Strohm, M., Eiblmeier, M., Langebartels, C., Jouanin, L., Polle, A., Sandermann, H., et al. (1999). Responses of transgenic poplar (*Populus tremula* x *P-alba*) overexpressing glutathione synthetase or glutathione reductase to acute ozone stress: visible injury and leaf gas exchange. *J. Exp. Bot.* 50, 365–374. doi: 10.1093/jxb/50.332.365

Strohm, M., Eiblmeier, M., Langebartels, C., Jouanin, L., Polle, A., Sandermann, H., et al. (2002). Responses of antioxidative systems to acute ozone stress in transgenic poplar (*Populus tremula* x *P-alba*) over-expressing glutathione synthetase or glutathione reductase. *Trees-Struct. Funct.* 16, 262–273. doi: 10.1007/s00468-001-0157-z

Vollenweider, P., Fenn, M. E., Menard, T., Günthardt-Goerg, M., and Bytnerowicz, A. (2013). Structural injury underlying mottling in ponderosa pine needles exposed to ambient ozone concentrations in the San Bernardino Mountains near Los Angeles, California. *Trees* 27, 895–911. doi: 10.1007/s00468-013-0843-7

Vollenweider, P., Günthardt-Goerg, M. S., Menard, T., Baumgarten, M., Matyssek, R., and Schaub, M. (2019). Macro- and microscopic leaf injury triggered by ozone stress in beech foliage (*Fagus sylvatica* L.). *Ann. For. Sci.* 76:71. doi: 10.1007/s13595-019-0856-5

Vollenweider, P., Menard, T., Arend, M., Kuster, T. M., and Gunthardt-Goerg, M. S. (2016). Structural changes associated with drought stress symptoms in foliage of central European oaks. *Trees-Struct. Funct.* 30, 883–900. doi: 10.1007/s00468-015-1329-6

Vollenweider, P., Ottiger, M., and Günthardt-Goerg, M. S. (2002). Validation of leaf ozone symptoms in natural vegetation using microscopical methods. *Environ. Pollut.* 124, 101–118. doi: 10.1016/s0269-7491(02)00412-8

Wittig, V. E., Ainsworth, E. A., Naidu, S. L., Karnosky, D. F., and Long, S. P. (2009). Quantifying the impact of current and future tropospheric ozone on tree biomass, growth, physiology and biochemistry: a quantitative meta-analysis. *Glob. Change Biol.* 15, 396–424. doi: 10.1111/j.1365-2486.2008.01774.x

Xu, Y. S., Feng, Z. Z., Shang, B., Dai, L. L., Uddling, J., and Tarvainen, L. (2019). Mesophyll conductance limitation of photosynthesis in poplar under elevated ozone. *Sci. Total Environ.* 657, 136–145. doi: 10.1016/j.scitotenv.2018.11.466

Zhang, J., Schaub, M., Ferdinand, J. A., Skelly, J. M., Steiner, K. C., and Savage, J. E. (2010). Leaf age affects the responses of foliar injury and gas exchange to tropospheric ozone in *Prunus serotina* seedlings. *Environ. Pollut.* 158, 2627–2634. doi: 10.1016/j.envpol.2010.05.003

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Turc, Vollenweider, Le Thiec, Gandin, Schaub, Cabané and Jolivet. *This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.*