# Powdery mildew decreases the radial growth of oak trees with cumulative and delayed effects over years

Didier Bert, Jean-Baptiste Lasnier, Xavier Capdevielle, Aline Dugravot, Marie

Laure Desprez-Loustau

► To cite this version:

Didier Bert, Jean-Baptiste Lasnier, Xavier Capdevielle, Aline Dugravot, Marie Laure Desprez-Loustau. Powdery mildew decreases the radial growth of oak trees with cumulative and delayed effects over years. PLoS ONE, 2016, 11 (5), 19 p. 10.1371/journal.pone.0155344. [hal-02634265](https://hal.inrae.fr/hal-02634265v1)

HAL Id: hal-02634265

https://hal.inrae.fr/hal-02634265v1

Submitted on 27 May 2020

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 CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Powdery Mildew Decreases the Radial Growth of Oak Trees with Cumulative and Delayed Effects over Years

Didier Bert¹*, Jean-Baptiste Lasnier², Xavier Capdevielle¹, Aline Dugravot³, Marie-Laure Desprez-Loustau¹

1 BIOGECO, INRA, Univ. Bordeaux, 33610 Cestas, France, 2 DRIAAF ILE DE France, F-94234 Cachan, France, 3 Inserm U1018—F-94800 Villejuif, France

* didier.bert@pierroton.inra.fr

## Abstract

Quercus robur and Q. petraea are major European forest tree species. They have been affected by powdery mildew caused by Erysiphe alphitoides for more than a century. This fungus is a biotrophic foliar pathogen that diverts photosynthetate from the plant for its own nutrition. We used a dendrochronological approach to investigate the effects of different levels of infection severity on the radial growth of young oak trees. Oak infection was monitored at individual tree level, at two sites in southwestern France, over a five-year period (2001–2005). Mean infection severity was almost 75% (infected leaf area) at the end of the 2001 growing season, at both sites, but only about 40% in 2002, and 8%, 5% and 2% in 2003, 2004 and 2005, respectively. Infection levels varied considerably between trees and were positively related between 2001 and 2002. Increment cores were taken from each tree to assess annual ring widths and increases in basal area. Annual radial growth was standardised to take the effect of tree size into account. Annual standardised radial growth was significantly and negatively correlated with infection severity in the same year, for both 2001 and 2002, and at both sites. The decrease in growth reached 70–90% for highly infected trees. The earlywood width was poorly correlated with infection severity, but the proportion of latewood in tree rings was lower in highly infected trees (60%) than in less heavily infected trees (85%). Infection in 2001 and 2002 was found to have a cumulative effect on radial growth in these years, together with a delayed effect detectable in 2003. Thus, even non-lethal pathogens like powdery mildew can have a significant impact on tree functioning. This impact should be taken into account in growth and yield models, to improve predictions of forest net primary production.

## Introduction

The impact of abiotic factors on forest productivity constitutes a highly active broad field of investigation, particularly in the context of climate change [1, 2]. Dendrochronology analyses of tree-ring width provide an integrative approach to the assessment of this impact [3–7].

Far fewer studies have investigated the effects of biotic factors. Strong impacts of biotic factors have been reported for seedling or tree survival (e.g. [8–10]), and for tree growth in severe

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

OPEN ACCESS

**Citation:** Bert D, Lasnier J-B, Capdevielle X, Dugravot A, Desprez-Loustau M-L (2016) Powdery Mildew Decreases the Radial Growth of Oak Trees with Cumulative and Delayed Effects over Years. PLoS ONE 11(5): e0155344. doi:10.1371/journal.pone.0155344

**Editor:** Mark Gijzen, Agriculture and Agri-Food Canada, CANADA

**Received:** November 27, 2015

**Accepted:** April 8, 2016

**Published:** May 13, 2016

**Copyright:** © 2016 Bert et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** Tree-ring width data are available from the NOAA National Centers for Environmental Information www.ncdc.noaa.gov/data-access/paleoclimatology-data/datasets/tree-ring. Basal Area Increments are included in S4 Table together with powdery mildew infection per tree.

**Funding:** The study was funded by French institute INRA. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
**Competing Interests:** The authors have declared that no competing interests exist.

outbreaks of pest and pathogen infestations [11, 12]. More limited or chronic effects of biotic factors on tree growth have been investigated in only a few pathosystems, such as *Dothistroma* needle blight on *Pinus radiata* in New Zealand [13–14], Swiss needle-cast caused by *Phaeo-cryptopus gaeumannii* on *Pseudotsuga menziesii* in New Zealand [15] or North West America [16], *Armillaria* root rot on various conifers [14, 17, 18] or *Mycosphaerella* leaf blights on *Euca-lyptus globulus* [19]. As regards *Quercus robur*, only the effects of defoliation by insects have been investigated on mature trees [20].

Relating pathogen infection to tree growth is especially important when dealing with forest plantations where disease translates into yield loss [13–15, 19]. However, decreased growth caused by pathogens can more generally affect all ecological functions and services provided by forests, in particular their role in biogeochemical cycles. Up to now, long-term growth reductions due to pathogens are not included in ecosystem models which could lead to inaccurate estimations of carbon uptake and storage [21]. Several authors have emphasized the need of “adding pests into the equation” for forest tree growth modelling [21–23]. Among various types of pathogen damage, foliar infections are likely the easiest to model a quantitative impact on tree growth since they directly reduce the photosynthetically active area of trees [21]. However, predicting the outcome of foliar attack on host growth is not so straightforward, because non-linear relationships between foliar infection and productivity losses are expected to occur at different scales [22–24]. Studies on whole leaves have shown that the decrease in photosynthesis resulting from leaf damage by parasites is generally greater than could be accounted for by the direct effects of the loss of leaf area alone [25–27]. On the other hand, compensatory effects may occur at the whole-plant level [23, 26]. More investigations are thus needed to improve our understanding of the relationships linking foliage damage and biomass production, and to enable pathogen damage to be included in forest tree growth modelling [21–23].

This study focuses on oak powdery mildew. Oaks, including the pedunculate oak *Quercus robur* L. and the sessile oak *Q. petraea* (Matt.) Liebl., are major tree species in European forests. They are widely distributed and are considered to be of high economic, cultural and environmental value [28]. Since the beginning of the 20th century, European oaks have been affected by powdery mildew caused by *Erysiphe alphitoides* (previously known as *Microsphaera alphitoides*) and related species [29]. Powdery mildew fungi are obligate biotrophs that absorb nutrients from the living cells of their host plant, without killing them, through differentiated structures called haustoria [30]. They thus divert photosynthetate from the plant for their own nutrition [31]. Studies on the physiological effects of powdery mildew infection on oak leaves have reported negative effects of the disease on net CO₂ assimilation rates [32], although the effects were moderate, with photosynthesis maintained even when the leaf surface was almost completely covered with mycelium [33]. Powdery mildew has particularly damaging effects on young seedlings, in nurseries, plantations and natural regeneration conditions, and on mature trees, in conjunction with other factors leading to decline [34]. Decreases and changes in growth have been reported in young pedunculate oak seedlings [35]. However, no estimates of growth loss in older moderately infected trees are currently available. The aim of this study was to assess the impact of powdery mildew on the radial growth of oak trees, by studying the relationship between levels of foliar infection and growth increments, using a dendrochronological approach.

# Materials and Methods

## Study Site and Oak Stands

Data was obtained from two *Quercus robur* L. field experiments in the Landes de Gascogne region in south-western France. The experiments were in community-owned forests and an agreement allowing samples to be taken from was signed with INRA. One experiment was
located near Biscarosse ("B1 site" at 44°26N-1°04W, 40m a.s.l, S1 Fig) and the other near Pon-
tenx-les-Forges ("P8 site" at 44°17N-1°04W). Oaks were planted in November 1994, after two
years in a nursery. Each experiment site consisted of six rows of 17 oaks with 2 m between adja-
cent trees in each row and 3 m separating the rows. The rows corresponded to three geographic
provenances of Q. robur: two rows were planted with seedlings raised from acorns collected in
the Vouzeron forest (47°15N-2°13E) located 400 km northeast of the study site, two rows were
of local provenance, from the Hermitage forest (44°44N-0°46W), and two rows were of local
provenance originating from a forest 75 km further south, at Bias (44°08N-1°13W). Provenances
were planted in unreplicated blocks at each site meaning that provenance effects and main-plot
effects were confounded. However, the small size of each plot and the homogeneity of the field
could be expected to reduce this disadvantage. This region has a temperate-maritime climate,
with an annual mean temperature of 13.6°C and about 940 mm precipitation per year (mean:
1980-2010). P8 is drier than B1, because the hydromorphic sandy spodosol at this site has a rela-
tively cemented spodic horizon (dense podzol) over 100% of the area at a depth of 43 cm,
whereas the iron pan covers only 44% of the surface at B1 and is located at a depth of 58 cm.

Monitoring of Infection and Tree Growth

Oak growth and powdery mildew infection were monitored at tree level at both sites in 2001 and 2002, using additional observations in 2003–2005 for the P8 site. Shoot growth phenology was monitored by weekly assessments in spring, with a seven-point scale: 0 = dormant bud, 1 = swollen bud, 2 = leaf tips extended beyond bud scales, 3 = first leaves unfolded, 4 = leaves unfolded but not yet full-size, 5 = leaves fully expanded but still tender and light green, 6 = mature leaves (dark green and tough). In 2001 at the B1 site and in 2002 at both sites, late frost caused leaf damage (mid-May assessment). The percentage of leaves affected was recorded on each tree.

Powdery mildew infection was assessed at tree level during the growing season, by trained observers. Initial (primary) infections occur in early spring as small whitish circular lesions on a few leaves, mostly visible on the upper leaf surface. Progressively, an increasing number of leaves are infected, and, as secondary infections occur, the whole leaf surface may be covered by mycelium and spores, resulting in a whitish appearance. The percentage of the leaf area affected by oak powdery mildew (infection severity) was estimated as follows. The leaves of each tree were examined by eye and assigned to one of four damage groups: "0" for leaves with no powdery mildew, "A" if less than 50% of the surface area of the leaves was infected, "B" if the infected area covered more than 50%, "C" for severely distorted and/or necrotic or dead leaves. The whole tree infection index (Inf) was then calculated as follows:

$$
\text{Inf} = 0.25 (\% \text{ leaves in group A}) + 0.75 (\% \text{ leaves in group B}) + 1 (\% \text{ leaves in group C})
$$

Inf therefore varies from 0 to 1, corresponding to a percentage of infected leaf area from 0 to
100%.

In 2001, both sites were monitored nine times, from May 18th to August 23rd. A mean infec-
tion index was calculated for the first flush, from the first four assessments until mid-June
("spring infection"). We also calculated another mean infection index for the end of June to the
end of September ("summer infection"). In the other years, infection levels in the spring were
much lower (see below) so assessments were carried out only during the summer.

Dendrochronology

In April 2013, an increment core was obtained 0.30 m above ground level, from all the surviv-
ing oaks, on the main stem and up to two additional stems; individuals with a bushy
appearance, with no main stem and more than three stems were discarded. The cores were planed with a cutter and digitalised on a scanner at a resolution of 1600 dpi. Ring widths were measured to the nearest 0.01 mm on the image, with dendrochronology software (Windendro™, Regent Instruments Inc., Canada). The series were cross-dated with the Gleichläufigkeit index in Windendro software (i.e., the percentage of cases of agreement in the signs of the first-difference of two time series). We rejected the few samples with inner tree rings dating from 1998 or later due to apex death followed by re-sprouting, or for which the core missed the pith. This selection procedure resulted in the use of trees with at least three measured tree rings dating from before 2001. The final dataset consisted of cores from 67 oaks at the B1 site and 61 oaks at the P8 site.

The ring widths (in mm) were then converted into basal area increment (BAI in cm²) values, to provide a better estimate of total radial growth. For each year, the basal area (BA) was calculated as

$$BA_n = \pi r_n^2$$

with $r_n$ the sum of ring widths from the pith to year n inclusive. Basal area increments (BAI) were then calculated as:

$$BAI_n = BA_n - BA_{n-1}$$

For multi-stem trees, the sum of $BA_n$ for the different stems was calculated per year, and $BAI_n$ was then calculated as if the tree had one stem. BA values ranged from 9 to 190 cm² in 2012. The largest oak sampled therefore had a BA 21 times larger than that of the smallest oak.

## Data Analysis

The summer powdery mildew infection index was analysed using a linear mixed model including site, provenance, year (2001 and 2002) and their interactions as fixed effects and individual trees as a random factor. A preliminary analysis showed that there was no significant interaction between year and site or year and provenance. The infection model was therefore defined as follows:

$$Inf_{ijkt} = \mu + \alpha Site_i + \beta Provenance_j + \gamma Site_i * Provenance_{ij} + \delta Year_t + Tree_{ijk} + \varepsilon_{ijkt} \quad (1)$$

where $Inf_{ijkt}$ is the summer mean infection index for site "i", provenance "j", tree "k", in year "t", $\mu$ is the overall mean. $\alpha$, $\beta$, $\gamma$, $\delta$ are the parameters corresponding to the fixed effects of site, provenance, their interaction, and year, $Tree_{ijk}$ is the random tree effect and $\varepsilon_{ijkt}$ is the residual error such as $Tree_{ijk} \sim N(0, \sigma^2_{Tree})$ and $\varepsilon_{ijkt} \sim N(0, \sigma^2_{\varepsilon})$, and $Tree_{ijk}$ are independent of $\varepsilon_{ijkt} \forall i,j, k,t$ and $\varepsilon_{ijkt}$ are independent $\forall i,j,k,t$.

Multiple comparisons of means were carried out between provenances at a single site, for a given year, with the Bonferroni test. Linear models were used to investigate the relationship between growth and powdery mildew infection. Initial tree size was included in the model because large trees generally grow faster than small ones.

We first analysed effects of infection on radial increments on an annual basis, for 2001 and 2002. Infection levels were too low in subsequent years to perform such analyses. The $BAI_n/BA_{n-1}$ ratio (expressed as a percentage) describing the standardised radial growth of the tree was used as the independent variable. It was assumed that there was a general linear relationship with no intercept between the basal area increment of one year and the basal area of the stem in the previous year (S2 Fig, [36]). This ratio takes into account the wide range of sizes of the trees studied and makes it possible to compare their growth during a given year (i.e., 2001 or 2002). We assumed that infection had a linear effect on this ratio, as described by the
following model:

$$ \frac{BAI_n}{BA_{n-1}} = \beta_0 + \beta_1 Inf_n + \varepsilon \quad (2) $$

where $\beta_0$ is the intercept and $\beta_1$ is the slope. As the value of the slope of the regression line depends on the range of BAI$_n$/BA$_{n-1}$, comparisons between years and between sites required the use of a standardised slope coefficient (SSC). When Inf increases by 1 unit, BAI$_n$/BA$_{n-1}$ decreases by $\beta_1$ (actually its absolute value, assuming $\beta_1$ is negative) and this decrease can be expressed as a percentage of $\beta_0$:

$$ SSC = 100 * \frac{|\beta_1|}{\beta_0} \quad (3) $$

For a given year, when SSC is 0.75, for example, oaks with an infection severity of 100% display 75% less growth than uninfected oaks. We used the delta method to computed 95% CI of SSC. The delta method approximates the standard errors using a first-order Taylor approximation [37].

We also analysed the cumulative effects of infection in 2001 and 2002 on growth over these two years, as follows:

$$ \frac{BAI_{2001} + BAI_{2002}}{BA_{2000}} = \beta_0 + \beta_1 Inf_{2001+2002} + \varepsilon \quad (4) $$

with $Inf_{2001+2002}$ the mean of the summer infection indices for 2001 and 2002

We further analysed the cumulative effects of infection by partitioning the whole sample of trees for the 2002 analyses into two subsamples with the median value of infection severity in 2001 used as the cut-off point. Model (2) for 2002 was then tested for the two subsamples.

Finally, we explored delayed effects of infection, by looking at effects on growth in the years following a season with infection. Our study period consisted of two successive years (2001 and 2002) with high powdery mildew infection followed by several years with low levels of infection (see below). We therefore assessed the cumulative effects of infection in 2001 and 2002 on radial growth in 2003 and subsequent years, as follows:

$$ \frac{BAI_n}{BA_{2000}} = \beta_0 + \beta_1 Inf_{2001+2002} + \varepsilon \quad (5) $$

with $n = 2003, 2004$ and $2005$.

For model (2), we determined whether the mid-May frost had an effect on the relationship between powdery mildew infection and radial growth, by running the model in 2001 for the B1 site for two equally distributed classes with and without frost damage.

Statistical analyses were carried out and graphs were produced using the lm(), lme() and deltamethod() functions in the nlme and msn libraries of R [38].

# Results

## Infection Monitoring

Budburst occurred during late April in all years (see Fig 1 for the P8 site). Very similar dynamics were observed at the two sites for all variables and years. Late frost damage was observed in mid-May in 2001 and 2002. In 2001, only the B1 site was affected, with 53% of trees and 30% of leaves affected on average. More generalised frost damage to leaves was observed in 2002, with 99% of trees affected at both sites, and a mean of 41% of the leaves frozen at B1, and 36% at P8.
In 2001, powdery mildew infections increased very rapidly after budburst and almost all leaves were infected at the end of May. Mean infection severity was estimated at about 50% of the leaf surface area during June-July and reached almost 75% at the end of the season, at both sites (Fig 1). Infection dynamics were slower, resulting in the persistence of infection at much lower levels, in subsequent years. In 2002, infection was delayed until refoliation after frost damage and progressed slowly through the summer, until a mean of about 40% of the leaf surface area was infected. No powdery mildew epidemics were observed in 2003, 2004 and 2005, resulting in very low final infection levels: 8%, 5%, and 2%, respectively (Fig 2A).

The linear mixed model of infection (model 1) showed the following significant effects: site ($p$ = 0.0157), provenance ($p$ = 0.0002), year ($p$<0.0001), interaction between site and provenance ($p$ = 0.0159), and tree effect (variance of the random tree effect, $p$< 0.0001). No significant difference in infection was observed between sites in 2001 with mean summer infection rates of 62% at B1 and 59% at P8 ($t$ = 0.939, DF = 118, $p$ = 0.349). However, a significant difference occurred in 2002, as infection rates were slightly lower at P8 than at B1 with a mean of 32% at B1 and 25% at P8 ($t$ = 2.935, DF = 124, $p$ = 0.004). The provenance effect resulted mostly from higher rates of infection among trees originating from the Vouzeron forest 400 km north of the study site than in trees with a local provenance, particularly those from Bias, with significant differences recorded at B1 (Fig 3; Table 1). The ranking of the Bias and Hermitage provenances differed between sites, resulting in a significant interaction between site and provenance in the analysis of variance. The strongest effect in the model was that of year, due to a large difference in mean infection levels between years (Figs 2A and 3).

Tree infection levels were positively related between years (Fig 4). The Pearson coefficient was 0.764 ($p$<0.0001; df = 61) for the B1 site and 0.333 ($p$ = 0.0087; df = 59) for the P8 site. At the P8 site, trees with lower infection levels in 2001 also showed low levels of infection in 2002, whereas highly infected trees in 2001 displayed a wide range of infection level in 2002.

## Effects of Infection on Growth

Soil characteristics differed between the sites, with a lower potential for growth at the P8 site. The oaks at P8 had a mean basal area only 50–70% of that at B1, from 2000 to 2005. The basal area increment (BAI) varied from 2 to 3 cm² over years and sites, with large differences between trees. It was significantly lower at the P8 site than at the B1 site in 2000, 2003 and 2004 (Fig 2B).

Standardised radial growth was significantly negatively correlated with infection severity in 2001 and 2002, at both sites (Fig 5). In 2001, the correlation with spring infection severity (Fig 5A and 5B) was weaker than that with summer infection severity (Fig 5C and 5D), at both sites (Table 2). For this year, the earlywood and latewood widths were distinguished within the annual ring, and the latewood/total ring width ratio was then calculated as a percentage. The correlation between this ratio and powdery mildew infection in 2001 was negative and highly significant: B1 site r = -0.509, $p$<0.0001 (S3 Fig) and P8 site r = -0.397, $p$ = 0.0009. The observed growth reduction concerned principally the latewood, as earlywood width remained fairly constant. The correlation between earlywood width and summer infection was significant for the B1 site (r = -0.260, $p$ = 0.0348) but not for the P8 site (r = -0.141, $p$ = 0.259). Tree rings consisted of 85% latewood in trees with lower levels of infection but only 60% latewood in highly infected trees.

The data were highly dispersed, but the relationship between infection and radial growth could be approximated by a linear regression over the range of infection severity, with no obvious threshold (Fig 5). The SSC was 0.7–0.9, with a mean of 0.785 in the case of severe summer infection in 2001 and 2002 (Table 2). The loss of radial growth may therefore be close to 78.5% at very high levels of powdery mildew infection, regardless of the site and the year studied.
<figure><img src="image_4.png" /><figcaption>Fig 1. At the P8 site: the green line indicates phenology (median score, from 0 = dormant bud to 6 = fully developed leaves), and the red line indicates powdery mildew severity (mean % of the leaf surface infected). The blue arrow indicates the date of the late frost.</figcaption></figure>

doi:10.1371/journal.pone.0155344.g001
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Final infection severity and basal area increment (BAI), by year and site. The boxplot shows the median in black, the mean in red, the first and third quartile, and the outliers. *n* = 67 oaks at B1 (grey) and *n* = 61 oaks at P8 (white).</figcaption></figure>

doi:10.1371/journal.pone.0155344.g002

The cumulated effects of infections over 2001 and 2002 on cumulated radial growth for these two years were even stronger than the annual effects estimated separately. Very high standardised slope coefficients of 0.90±0.25 (95% C.I.) at Site B1 and 1.27±0.25 at Site P8 were obtained, corresponding to more than 90% growth reduction for 100% infection over the two years (**Fig 6A and 6B; Table 2**).

The cumulative effect on radial growth of infection in 2001 and 2002 was also demonstrated by the observation that the effect of infection in 2002 on radial growth in the same year depended on the level of infection in the previous year at the P8 site (Fig 7; Table 2). Only trees highly infected in 2001 displayed a significant relationship between infection and radial growth in 2002. A similar analysis was not possible at the B1 site because the strong correlation between infection levels in 2001 and 2002 meant that trees highly infected in 2002 were invariably also highly infected in 2001.

<figure><img src="image_6.png" /><figcaption>**Fig 3.** Mean infection index, by site, year and provenance. The black lines and dots are for the B1 site, the blue lines and squares are for the P8 site. The continuous and dotted lines are for 2001 and 2002, respectively.</figcaption></figure>

doi:10.1371/journal.pone.0155344.g003
<table><caption>**Table 1.** Number of trees, mean infection index and standard deviation, by site, provenance and year. The means in bold differed significantly (*p*<0.05) for the year concerned.</caption><thead><tr><th rowspan="2">Year</th><th rowspan="2">Provenance</th><th rowspan="2">n</th><th colspan="3">B1 site</th><th colspan="2">P8 site</th></tr><tr><th>Mean</th><th>r</th><th>n</th><th>Mean</th><th>σ</th></tr></thead><tbody><tr><td>2001</td><td>Bias</td><td>22</td><td><strong>50.49</strong></td><td>16.24</td><td>24</td><td>59.02</td><td>15.39</td></tr><tr><td>2001</td><td>Hermitage</td><td>21</td><td>63.71</td><td>20.26</td><td>13</td><td>52.14</td><td>15.82</td></tr><tr><td>2001</td><td>Vouzeron</td><td>23</td><td><strong>71.87</strong></td><td>19.27</td><td>24</td><td>63.15</td><td>13.02</td></tr><tr><td>2002</td><td>Bias</td><td>22</td><td><strong>24.66</strong></td><td>9.22</td><td>24</td><td>23.49</td><td>9.72</td></tr><tr><td>2002</td><td>Hermitage</td><td>21</td><td>32.86</td><td>12.54</td><td>14</td><td>24.89</td><td>12.17</td></tr><tr><td>2002</td><td>Vouzeron</td><td>21</td><td><strong>37.35</strong></td><td>10.98</td><td>24</td><td>27.28</td><td>13.66</td></tr></tbody></table>

doi:10.1371/journal.pone.0155344.t001

We also observed a cumulative and delayed effect of infection in 2001 and 2002 on the radial growth achieved in 2003 (Fig 8). No further effects on radial growth were observed in 2004 and 2005 (Table 2).

Discussion

The impact of pathogens (and pests) on forests has been well documented for severe epidemics or outbreaks [9, 39, 40]. In addition to pathogen-induced mortality, impacts on tree growth have been assessed by comparing epidemic and non-epidemic years [41], or affected and healthy reference trees in the same diameter class [16] or infected trees and fungicide-treated controls [42, 43]. However, the more chronic impact on tree growth of moderate levels of infection has been investigated in less detail. In particular, functions describing the growth response to various intensities of infection at stand or tree level have been produced for only a few pathosystems, such as Swiss needle cast on Douglas fir in particular [14, 19, 44, 45].

We carried out detailed monitoring of powdery mildew infection in two young *Quercus robur* plantations in which natural epidemics occurred in two successive years, producing considerable variation in the level of infection between trees. Infection was assessed at individual tree level, by trained observers using a combination of disease prevalence (frequency) and severity (classes) criteria. *Erysiphe* species are epiphyllous. Infection is therefore conspicuous

<figure><img src="image_7.png" /><figcaption>**Fig 4.** Infection severity in 2002, as a function of severity in 2001, at the two sites.</figcaption></figure>

doi:10.1371/journal.pone.0155344.g004
<figure><img src="image_8.png" /><figcaption>Fig 5. Growth index BAI<sub>n</sub>/BA<sub>n-1</sub> as a function of infection index for two years and two sites.</figcaption></figure>

doi:10.1371/journal.pone.0155344.g005
<table><caption>**Table 2. Statistics for two sites and two years of powdery mildew infection.** Regression coefficients (intercept and slope of regression lines in figures), slope probability (p), adjusted r², 95% confidence intervals, standardised slope coefficient (SSC, see Eq 3), 95% confidence interval of SSC assessed with Delta method [37], and Pearson correlation coefficient (r). Note: r and the slope have equal probabilities.</caption><thead><tr><th>Fig</th><th>Dependent variable</th><th>Infection period</th><th>Site</th><th>Intercept</th><th>Slope</th><th>p slope</th><th>Adj r<sup>2</sup></th><th>SSC</th><th>±95% CI</th><th>r</th></tr></thead><tbody><tr><td>5A</td><td>BAI2001/BA2000</td><td>Spring 2001</td><td>B1</td><td>40.948</td><td>-0.305</td><td>&lt;0.0001</td><td>0.202</td><td>0.745</td><td>0.242</td><td>-0.463</td></tr><tr><td>5B</td><td>BAI2001/BA2000</td><td>Spring 2001</td><td>P8</td><td>80.536</td><td>-0.620</td><td>0.0129</td><td>0.085</td><td>0.770</td><td>0.382</td><td>-0.317</td></tr><tr><td>5C</td><td>BAI2001/BA2000</td><td>Summer 2001</td><td>B1</td><td>47.645</td><td>-0.320</td><td>&lt;0.0001</td><td>0.260</td><td>0.672</td><td>0.161</td><td>-0.521</td></tr><tr><td>5D</td><td>BAI2001/BA2000</td><td>Summer 2001</td><td>P8</td><td>111.133</td><td>-0.997</td><td>0.0001</td><td>0.204</td><td>0.897</td><td>0.212</td><td>-0.466</td></tr><tr><td>5E</td><td>BAI2002/BA2001</td><td>Summer 2002</td><td>B1</td><td>22.889</td><td>-0.185</td><td>0.0385</td><td>0.052</td><td>0.806</td><td>0.562</td><td>-0.259</td></tr><tr><td>5F</td><td>BAI2002/BA2001</td><td>Summer 2002</td><td>P8</td><td>34.511</td><td>-0.264</td><td>0.0343</td><td>0.058</td><td>0.764</td><td>0.560</td><td>-0.271</td></tr><tr><td>6A</td><td>(BAI2001+BAI2002)/BA2000</td><td>Mean 2001&amp;2002</td><td>B1</td><td>87.268</td><td>-0.788</td><td>&lt;0.0001</td><td>0.218</td><td>0.903</td><td>0.247</td><td>-0.480</td></tr><tr><td>6B</td><td>(BAI2001+BAI2002)/BA2000</td><td>Mean 2001&amp;2002</td><td>P8</td><td>207.088</td><td>-2.622</td><td>&lt;0.0001</td><td>0.277</td><td>1.266</td><td>0.246</td><td>-0.538</td></tr><tr><td>7</td><td>BAI2002/BA2001 Oaks with Inf 2001 ≤ 60.5</td><td>2002</td><td>P8</td><td>33.680</td><td>0.020</td><td>0.93 NS</td><td>0.035</td><td>0.060</td><td>1.333</td><td>0.017</td></tr><tr><td>7</td><td>BAI2002/BA2001 Oaks with Inf 2001 > 60.5</td><td>2002</td><td>P8</td><td>28.175</td><td>-0.230</td><td>0.0386</td><td>0.110</td><td>0.817</td><td>0.576</td><td>-0.373</td></tr><tr><td>8A</td><td>BAI2003/BA2000</td><td>Mean 2001&amp;2002</td><td>B1</td><td>51.945</td><td>-0.436</td><td>0.0046</td><td>0.110</td><td>0.839</td><td>0.349</td><td>-0.352</td></tr><tr><td>8B</td><td>BAI2003/BA2000</td><td>Mean 2001&amp;2002</td><td>P8</td><td>69.318</td><td>-0.730</td><td>0.0054</td><td>0.109</td><td>1.052</td><td>0.405</td><td>-0.352</td></tr><tr><td>NA</td><td>BAI2004/BA2000</td><td>Mean 2001&amp;2002</td><td>B1</td><td>52.959</td><td>-0.296</td><td>0.145 NS</td><td>0.019</td><td>0.559</td><td>0.553</td><td>-0.186</td></tr><tr><td>NA</td><td>BAI2004/BA2000</td><td>Mean 2001&amp;2002</td><td>P8</td><td>21.615</td><td>1.090</td><td>0.072 NS</td><td>0.038</td><td>5.043</td><td>17.181</td><td>0.232</td></tr><tr><td>NA</td><td>BAI2005/BA2000</td><td>Mean 2001&amp;2002</td><td>B1</td><td>25.719</td><td>-0.093</td><td>0.376 NS</td><td>0.003</td><td>0.363</td><td>0.665</td><td>-0.113</td></tr><tr><td>NA</td><td>BAI2005/BA2000</td><td>Mean 2001&amp;2002</td><td>P8</td><td>54.987</td><td>-0.127</td><td>0.629 NS</td><td>0.013</td><td>0.230</td><td>0.840</td><td>-0.063</td></tr></tbody></table>

doi:10.1371/journal.pone.0155344.t002

and visible symptoms are assumed to be strongly correlated with the amounts of tissue infected. We were therefore able to obtain a quantitative estimate of the percentage of the leaf area affected by the fungus for each tree. The validity of this assessment method was evaluated in preliminary trials (unpublished data). The almost Gaussian distribution of the tree disease severity values obtained for years with high levels of infection in this study and the significant

<figure><img src="image_9.png" /><figcaption>**Fig 6. Growth index as a function of mean infection index in 2001–2002 and regression lines for two sites.**</figcaption></figure>

doi:10.1371/journal.pone.0155344.g006
<figure><img src="image_10.png" /><figcaption>Fig 7. Growth index according to infection index in 2002, split into two groups on the basis of infection index for 2001 (low or high) at the P8 site.</figcaption></figure>

doi:10.1371/journal.pone.0155344.g007

positive correlations between trees obtained for independent assessments within and between years also suggest that this assessment method was appropriate. We used a dendrochronological approach to study the impact of infection on growth. This approach is useful for the reconstruction of growth dynamics over several years. Previous studies based on such an approach have compared the temporal growth dynamics of a particular species periodically affected by a pest or pathogen with those of a reference species growing in the same location but not

<figure><img src="image_11.png" /><figcaption>Fig 8. Growth index in 2003 as a function of mean infection index in 2001–2002 for two sites.</figcaption></figure>

doi:10.1371/journal.pone.0155344.g008
susceptible to the pest or pathogen concerned [46, 47]. However, these studies focused mainly on temporal patterns of infection and growth. Other studies have addressed the effect of a pathogen on the radial growth of trees by taking into account the stem size distribution or growth history [48]. The healthy or infected trees were often paired together on the basis of stem size or competition level, resulting in discrete classes, particularly for root pathogens [49, 50]. We aimed to assess the effect of infection over a continuous range. We therefore adjusted for initial tree size in each year, making it possible to use the whole range of stem size for the sample studied without the need to pair individuals or to form classes.

This study of more than one hundred trees at two sites clearly demonstrated a significant negative effect of powdery mildew infection on the radial growth of oak trees. This negative effect was demonstrated in two years (2001 and 2002), differing in infection level and climate. In 2001, severe oak powdery mildew epidemics occurred throughout southwestern France [51], but the spring and summer precipitation was high which was beneficial to growth (Fig 2). Such conditions proved favourable for showing a significant effect of infection. In 2002, the growing season was rather dry, resulting in lower growth, and infection severity was lower, but between-tree variability remained high at our experimental sites and there was also a significant infection effect. A potential cumulative effect of leaf frost damage and powdery mildew infection might partly account for the observed growth patterns in our experiment, because late frost and powdery mildew occurred together in three of the four combinations of site and year situations studied, and the oaks with the highest levels of infection were also those displaying frost damage (S4 Fig). However, in 2001, the powdery mildew-growth relationship was similar between P8, where no frost damage was observed, and B1, where significant frost damage occurred. High levels of powdery mildew infection associated with oak refoliation after frost damage have frequently been reported [34]. In our experiment, frost damage may therefore exert its deleterious effects principally by increasing rates of powdery mildew infection, without modifying the effect of infection on tree growth.

The quantitative impact of *E. alphitoides* on oak growth has previously been described only for seedlings (reviewed in [34]). For example, height increment losses of 30–40% were observed in an experiment in south-west France in which naturally infected seedlings (about 30% of leaf area infected) were compared with control seedlings treated with fungicide [35]. Soutrenon [52] also reported large benefits of experimental fungicide treatments in naturally regenerating oak, in terms of seedling survival and biomass (greater increases in height and diameter than for untreated seedlings). The impact of foliar pathogens on the growth of older trees has been documented for other pathosystems. A 50–70% decrease in growth was reported following *Phaeocryptopus gaeumannii* infections in Douglas fir stands in southern Germany (Merkle 1951 in [14]), and a 27% decrease in basal area increment was reported in New Zealand [14]. Eucalypt growth following attacks by several foliar pests and pathogens has also been shown to be highly variable, with no effect until approximately 50% defoliation was reached in some cases, but a significant effect with as little as 10% defoliation in other situations [21].

In this study, a linear model fitted the data well for both years and sites, suggesting a proportional effect of infection (expressed as the percentage of the leaf area affected) on radial growth and the absence of a threshold effect. This contrasts with the findings of other studies, in which only high levels of infection were shown to affect growth [21, 53]. For example, the effects of the foliar pathogen *Dothistroma pini* on the growth of radiata pines were observed only at infection severities exceeding 20–30% in three- to four-year-old *Pinus radiata* trees [53, 54]. The growth reduction caused by two other needle cast diseases, *Lophodermium* sp. on pines and *Phaeocryptopus gaumannii* on Douglas fir, has also been shown to be less than 1:1 with respect to foliar disease indices [53]. However a more detailed analysis of our data also suggested that moderate levels of infection in a single year may have a limited effect. Linear models
fitted the data particularly well (high R²) at both sites in 2001, when infection levels ranged from 30% to 100%. In 2002, with infection severity in the range 10–50%, an impact on growth was observed only on trees already severely affected in 2001 at the P8 site.

This cumulative effect of infection in successive years was clearly shown by even stronger infection-growth relationships in the analyses based on a two-year period than in separate estimations of annual effects. Very high standardised slope coefficients were obtained, exceeding 0.9 and corresponding to a 90% growth reduction for 100% infection over the two years. A similar effect of successive defoliations by oak herbivore insects on the relationship between increment loss and defoliation index has been reported [20]. The cumulative effect of infection in two successive years may be accounted for by delayed effects of the infection observed in 2001 on growth in the following year. First-order autocorrelations have been shown to be strong in tree-ring series in *Quercus*, because the climate in year$_n$ also determines some of the growth in year$_{n+1}$ [55]. Similarly, we observed carry-over effects of infection on growth in subsequent years. The strong correlations between infection rates in 2001 and 2002 made it difficult to separate the effects of these two years, although we were able to go some way towards this goal at P8, as discussed above. Nevertheless, use of the whole sequence of years, from 2001 to 2005, showed that infection in 2001–2002 had an effect until the end of the 2003 growing season. The absence of a relationship to tree growth in subsequent years suggests that trees recover after one or two years, as often reported for insect defoliations [12]. These findings also suggest that growth loss was a direct effect of powdery mildew infection and not an artefact due to inherent growth differences between trees displaying different levels of disease severity.

The cumulative effects of infection are probably favoured by the consistent behaviour of trees between years, as shown by the strong correlation between the infection levels in 2001 and 2002 for individual trees, and the significant tree effect in the infection model. Several factors may explain such individual effects. We observed a provenance effect, with higher levels of infection on trees originating from Vouzeron, the provenance furthest away from the experimental site. The greater resistance of local provenances may be accounted for by better local adaptation of oaks, although the pathogen has been shown to benefit more from local adaptation than the host in this pathosystem [56, 57]. Genetic effects on tree susceptibility to disease are likely, in the form of direct effects [35] or through host phenology, a highly heritable trait [58] that also influences vulnerability to late frost [59].

The effects of disease on growth depend on the mechanisms connecting infection to growth processes at various temporal and spatial scales, from the physiological processes in infected tissues to processes at tree level. For example, the impact of conifer foliar pathogens is highly dependent on the age classes of the needles infected, conditioning their contribution to photosynthesis [53]. The pattern of infection/defoliation within a tree is also very important for determining growth responses and may contribute to the differences in such responses between trees of different ages [21, 41]. Infection also impairs growth by other more complex and diverse mechanisms, such as foliage retention, the occlusion of stomata, changes in light penetration, foliar nutrient translocation, gaseous exchanges, and carbohydrate dynamics [21, 26, 45, 53, 60–62]. The impact of broadleaf foliar biotrophic pathogens, such as powdery mildews, should be directly related to the percentage of the leaf area infected, particularly in young trees, because these pathogens develop on photosynthetically active leaves and hijack photosynthate for their own development via their haustoria, which tap into living leaf cells [30, 63, 64]. E. alphitoides has been shown to decrease the net CO₂ assimilation rate in infected leaves [31, 33], through several mechanisms, including decreases in stomatal conductance [33] and the quantum efficiency of photosystem II [65]. Changes in carbohydrate translocation patterns, lower rates of photorespiration and higher rates of dark respiration have also been demonstrated [31–33]. The effect of foliar pathogens on the net photosynthetic rate of leaves usually exceeds
that predicted from the lesion area, because the physiology of the remaining green leaf tissue is also affected [24]. Compensatory mechanisms in the healthy parts of the leaf may account for the impact of lesion size being smaller than expected in some cases ([25] for *Chrysomyxa rho-dodendri* on *Picea abies*). By contrast, changes in translocation due to heavy infection can lead to cell death and leaf shedding, resulting in an even greater carbon requirement for crown restoration [64].

Interannual variation in the radial growth of oaks in relation to climate has been shown to be due mostly to a decrease in latewood production during dry periods, with the production of earlywood remaining fairly constant [55]. Rubtsov [20] also reported a closer relationship of latewood increment than of earlywood increment to insect defoliation. Similarly, we show here that the effect of powdery mildew infection was principally due to a decrease in latewood width in the 2001 tree ring. Moreover, the impact on annual growth was more strongly related to summer infection (overall severity) than to spring infection (precocity of infection). This may be because the large vessels of the early wood are formed two weeks before budburst and are probably constructed from the carbon pool stored during the previous growing season [66]. The intensity of infection over the whole season therefore seems to be more important than the timing of infection, in terms of the impact on annual tree growth. The change in earlywood/ latewood ratio in infected trees also indicates that powdery mildew has both qualitative and quantitative effects on wood production, as wood density depends on the proportion of latewood [67].

## Conclusion–Implications

Unlike root pathogens, which spread throughout a forest stand and persist indefinitely, foliar pathogens, such as *E. alphitoides* cause peaks of infection only in some years. Cyclic infestations have been reported for several pest insects [68]. In pathogens, between-year variations have often been shown to be related to climate variability, favouring fungus and disease development to different extents [69]. This is particularly true for oak powdery mildew, epidemics of which have been shown to be correlated with mild winters, such as that in 2001 [34]. In the context of global warming, favourable years for infection are therefore likely to become more frequent, increasing the impact of powdery mildew. Furthermore, this impact is likely to be increased by interactions with abiotic stresses, such as frost and drought, which are also likely to increase in frequency. The impact of pathogens should therefore be included in growth and yield models, to improve predictions of forest net primary production and carbon balance [22, 27, 39].

## Supporting Information

S1 Fig. Map of France showing the two study sites (yellow and blue dots) and three oak provenances (red dots). Map plotted with R software.
(TIF)

S2 Fig. For both sites, Basal Area Increment in 2001 according to Basal Area of the tree at the end of year 2000. The oaks were split in two groups of low or high infection when they were less or more infected than 63.8 or 62.9%, for B1 or P8 respectively. The regression lines give the relationship between the annual radial growth and the size of the tree. The model was: $BAI = \mu + \alpha BA + \beta(BA * \text{Infection}) + \varepsilon$. BA effect and the interaction between BA and Infection were significant (p<0.0001) for both sites, i.e., the slopes were different because the infection decreased significantly the radial growth all along the tree size range.
(TIF)
**S3 Fig. Percentage of latewood in the 2001 tree ring as a function of infection levels in 2001 for the B1 site.**

(TIF)

**S4 Fig. Growth index according to infection in 2001, for oaks with and without frost damage in spring 2001.** The crosses indicate the means for the two groups.

(TIF)

**S1 Table. Data: one line per tree.** Columns: Site, Provenance, Line number in the plot, Oak number in the line, Frost index in 2001 and 2002, Mean Infection index from May 18 to June 22 in 2001, Mean Infection from June 29 to August 23 in 2001 (taken as Infection index in 2001 in the data analyses), Infection index per year over 2002 to 2005, Basal Area per year over 1994 to 2012, Basal Area Increment per year over 1994 to 2012.

(XLSX)

**Acknowledgments**

We would like to thank the *Unité expérimentale Forêt Pierroton*, UE 0570, INRA, 69 Route d'Arcachon, 33612 CESTAS and H. Jactel and I. Van Halder, for the installation of the experimental sites, L. Riols and F. Bernier for soil characterisation and dendrometric measurements, G. Capron for his contribution to powdery mildew monitoring, J. Sappa of Alex Edelman & Associations for language editing, and two reviewers for their suggestions.

**Author Contributions**

Conceived and designed the experiments: MLDL DB JBL. Performed the experiments: DB JBL XC AD MLDL. Analyzed the data: DB JBL MLDL AD. Contributed reagents/materials/analysis tools: XC DB MLDL. Wrote the paper: DB MLDL AD.

**References**

1. Lindner M, Maroschek M, Netherer S, Kremer A, Barbati A, Garcia-Gonzalo J, et al. (2010) Climate change impacts, adaptive capacity, and vulnerability of European forest ecosystems. Forest Ecology and Management 259: 698–709.

2. Bontemps JD, Bouriaud O (2014) Predictive approaches to forest site productivity: recent trends, challenges and future perspectives. Forestry 87: 109–128.

3. Fritts HC (1976) Tree rings and climate. London: Academic Press.

4. Cook ER (1987) The decomposition of tree-ring series for environmental studies. Tree-Ring Bulletin 47: 37–59.

5. Bert D (1993) Impact of ecological factors, climatic stresses and pollution on growth and health of silver fir (Abies alba Mill.) in the Jura mountains: an ecological and dendrochronological study. Acta Oecologica 14(2): 229–246.

6. Mérian P, Bontemps JD, Bergès L, Lebourgeois F (2011) Spatial variation and temporal instability in climate-growth relationships of sessile oak (*Quercus petreae* [Matt.] Liebl.) under temperate conditions. Plant Ecology 212: 1855–1871.

7. Guiot J, Boucher E, Gea-Izquierdo G (2014) Process models and model-data fusion in dendroecology. Frontiers in Ecology and Evolution 2(52): 1–12.

8. Augspurger CK (1984) Seedling survival of tropical tree species: interactions of dispersal distance, light-gaps, and pathogens. Ecology 65: 1705–1712.

9. Meentemeyer RK, Rank NE, Shoemaker DA, Oneal CB, Wickland AC, Frangioso KM, et al. (2008) Impact of sudden oak death on tree mortality in the Big Sur ecoregion of California. Biological Invasions 10: 1243–1255.

10. Chanthorn W, Caughlin T, Dechka S, Brockelman WY (2013) The relative importance of fungal infection, conspecific density and environmental heterogeneity for seedling survival in a dominant tropical tree. Biotropica doi: 10.1111/btp.12044: 1–7.
11. Iqbal M, MacLean DA, Kershaw JA (2011) Impacts of hemlock looper defoliation on growth and survival of balsam fir, black spruce and white birch in Newfoundland, Canada. Forest Ecology and Management 261: 1106-1114.

12. Sangüesa-Barreda G, Camarero JJ, Garcia-Martin A, Hernandez R, de la Riva J (2014) Remote-sensing and tree-ring based characterization of forest defoliation and growth loss due to the Mediterranean pine processionary moth. Forest Ecology and Management 320: 171-181.

13. Bulman LS (1993) *Cyclaneusma* needle-cast and *Dothistroma* needle blight in NZ pine plantations. New-Zealand Forestry Aug: 21-24.

14. Shaw CG III, Toes EHA (1977) Impact of *Dothistroma* needle blight and *Armillaria* root rot on diameter growth of *Pinus radiata*. Phytopathology 67: 1319-1323.

15. Kimberley MO, Hood IA, Knowles RL (2011) Impact of Swiss needle-cast on growth of Douglas fir. Phytopathology 101: 583-593. doi: 10.1094/PHYTO-05-10-0129 PMID: 20923368

16. Black BA, Shaw DC, Stone JK (2010) Impact of Swiss needle cast on overstory Douglas-fir forests of the Western Oregon Coast Range. Forest Ecology and Management 259: 1673-1680.

17. Battaglia M, Pinkard EA, Sands PJ, Bruce JL, Quentin A (2011) Modelling the impact of defoliation and leaf damage on forest plantation function and production. Ecological Modelling 222: 3193-3202.

18. Westwood AR, Conciatori F, Tardif JC, Knowles K (2012) Effects of *Armillaria* root disease on the growth of *Picea mariana* trees in the boreal plains of central Canada. Forest Ecology and Management 266: 1-10.

19. Cruickshank MG, Morrison DJ, Lalumiere A (2011) Site, plot, and individual tree yield reduction of interior Douglas-fir associated with non-lethal infection by *Armillaria* root disease in southern British Columbia. Forest Ecology and Management 261: 297-307.

20. Rubtsov VV (1996) Influence of repeated defoliations by insects on wood increment in common oak (*Quercus robur* L.). Annales des Sciences Forestières 53: 407-412.

21. Dietze MC, Matthes JH (2014) A general ecophysiological framework for modelling the impact of pests and pathogens on forest ecosystems. Ecology Letters doi: 10.1111/ele.12345: 1-9.

22. Pinkard EA, Battaglia M, Roxburgh S, O'Grady AP (2011) Estimating forest net primary production under changing climate: adding pests into the equation. Tree Physiology 31: 686-699. doi: 10.1093/treephys/tpr054 PMID: 21746746

23. Whitehead D (2011) Modelling the impacts of pests on forest productivity: a pathway through complexities and conundrums. Tree Physiology 31: 683-685. doi: 10.1093/treephys/tpr071 PMID: 21849589

24. Pinkard EA, Quentin A, Battaglia M, Wardlaw TJ, Mohammed CL (2007) The physiology of host responses to pest attack in *Eucalyptus globulus* plantations. Proceedings IUFRO Working Group 2.08.03 'Eucalypts and diversity: Balancing productivity and sustainability. Durban, South-Africa.

25. Bastiaans L (1991) Ratio between virtual and visual lesion size as a measure to describe reduction in leaf photosynthesis of rice due to leaf blast. Phytopathology 81: 611-615.

26. Mayr S, Siller C, Kriss M, Oberhuber W, Bauer H (2001) Photosynthesis in rust-infected adult Norway spruce in the field. New Phytologist 151: 683-689.

27. Pinkard EA, Mohammed CL (2006) Photosynthesis of *Eucalyptus globulus* Labill. with *Mycosphaerella* leaf disease. New Phytologist 170: 119-127. PMID: 16539609

28. Ducousso A, Bordacs S (2004) Technical guidelines for genetic conservation and use for pedunculate and sessile oaks (*Quercus robur* and *Q. petraea*). International Plant Genetic Resources Institute, Rome, Italy. Available from: http://www.bioversityinternational.org/e-library/publications

29. Mougou A, Dutech C, Desprez-Loustau ML (2008) New insights into the identity and origin of the causal agent of oak powdery mildew in Europe. Forest Pathology 38: 275-287.

30. Glawe DA (2008) The powdery mildews: a review of the world's most familiar (yet poorly known) plant pathogens. Annual Review of Phytopathology 46: 27-51. doi: 10.1146/annurev.phyto.46.081407.104740 PMID: 18680422

31. Hewitt HG, Ayres PG (1976) Effect of infection by *Microsphaera alphitoides* (powdery mildew) on carbohydrate levels and translocation in seedlings of *Quercus robur*. New Phytologist 77: 379-390.

32. Hewitt HG, Ayres PG (1975) Changes in CO₂ and water vapour exchange rates in leaves of *Quercus robur* infected by *Microsphaera alphitoides* (powdery mildew). Physiological Plant Pathology 7: 127-137.

33. Hajji M, Dreyer E, Marçais B (2009) Impact of *Erysiphe alphitoides* on transpiration and photosynthesis in *Quercus robur* leaves. Journal of Plant Pathology 125: 63-72.

34. Marçais B, Desprez-Loustau ML (2014) European oak powdery mildew: impact on trees, effects of environmental factors, and potential effects of climate change. Annals of Forest Science 71: 633-642.
35. Desprez-Loustau ML, Saint-Jean G, Barrès B, Dantec CF, Dutech C. (2014) Oak powdery mildew changes growth patterns in its host tree: host tolerance response and potential manipulation of host physiology by the parasite. Annals of Forest Science 71(5): 563-573.

36. Bowen AT Jr (1977) The relation of tree and stand characteristics to basal area growth of red spruce trees in partially cut stands in eastern Maine. Maine Agricultural Experiment Station Bulletin 627: 1-33.

37. Oehlert GW (1992) A note on the Delta method. The American Statistician 46: 27-29.

38. R Core Team (2013) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org/.

39. Hicke JA, Allen CD, Desai AR, Dietze MC, Hall RJ, Hogg EH, et al. (2012) Effects of biotic disturbances on forest carbon cycling in the United States and Canada. Global Change Biology 18: 7-34.

40. Boyd IL, Freer-Smith PH, Gilligan CA, Godfray HCJ (2013) The consequence of tree pests and diseases for ecosystem services. Science 342: 1235773. doi: 10.1126/science.1235773 PMID: 24233727

41. Hanso M, Drenkhan R (2012) *Lophodermium* needle cast, insect defoliation and growth responses of young Scots pines in Estonia. Forest Pathology 42: 124-135.

42. Carnegie AJ, Ades PK (2002) *Mycosphaerella* leaf disease reduces growth of plantation-grown *Euca-lyptus globulus*. Australian Forestry 66: 113-119

43. Gastine F, Berthelot A, Bouvet A, Roy B, Servant H (2003) La protection phytosanitaire du cultivar « Beaupré » est-elle efficace? Informations-Forêt 2(667).

44. Maguire DA, Kanaskie A, Voelker W, Johnson R, Johnson G (2002) Growth of young Douglas-fir plantations across a gradient in Swiss needle cast severity. Western Journal of Applied Forestry 17: 86-95.

45. Saffel BJ, Meinzer FC, Voelker SL, Shaw DC, Brooks JR, Lachenbruch B, et al. (2014) Tree-ring stable isotopes record the impact of a foliar fungal pathogen on CO₂ assimilation and growth in Douglas-fir. Plant, Cell and Environment 37: 1536-1547. doi: 10.1111/pce.12256 PMID: 24330052

46. Mott DG, Nairn LD, Cook JA (1957) Radial growth in forest trees and effects of insect defoliation. Forest Science 3: 286-304.

47. Lee EH, Beedlow PA, Waschmann RS, Burdick CA, Shaw DC (2013) Tree-ring analysis of the fungal disease Swiss needle cast in western Oregon coastal forests. Canadian Journal of Forest Research 43: 677-690.

48. Rentch J, Fajvan MA, Evans RA, Onken B (2009) Using dendrochronology to model hemlock woolly adelgid effects on eastern hemlock growth and vulnerability. Biological Invasions 11: 551-563.

49. Thies WG (1983) Determination of growth reduction in Douglas-fir infected by *Phelinus weirii*. Forest Science 29: 305-315.

50. Bendz-Hellgren M, Stenlid J (1995) Long term reduction in the diameter growth of butt rot affected Norway spruce, *Picea abies*. Forest Ecology and Management 74: 239-243.

51. Département Santé des Forêts (2002) Les Cahiers du D.S.F., La santé des forêts en France en 2001, Min. Agri. Pêche (DERF), Paris.

52. Soutrenon A (1998) Une experimentation pluri-annuelle confirme l'impact de l'oïdium sur de jeunes sujets. [A pluriannual experiment confirms the impact of powdery mildew on young seedlings]. Les cahiers du DSF 1-2000 (la santé des forets [France] en 1997, Paris):93-94.

53. Maguire DA, Mainwaring DB, Kanaskie A (2011) Ten-year growth and mortality in young Douglas-fir stands experiencing a range in Swiss needle cast severity. Canadian Journal of Forest Research 41: 2064-2076.

54. Gibson IAS (1974) Impact and control of *Dothistroma* blight of pines. European Journal of Forest Pathology 4: 89-100.

55. Lebourgeois F, Cousseau G, Ducos Y (2004) Climate-tree-growth relationships of *Quercus petraea* Mill. stand in the Forest of Bercé ("Futaie des Clos", Sarthe, France). Annals of Forest Science 61: 361-372.

56. Roslin T, Laine AL, Gripenberg S (2007) Spatial population structure in an obligate plant pathogen colonizing oak *Quercus robur*. Functional Ecology 21: 1168-1177.

57. Desprez-Loustau ML, Feau N, Mougou-Hamdane H, Dutech C (2011) Interspecific and intraspecific diversity in oak powdery mildews in Europe: coevolution history and adaptation to their hosts. Mycoscience 52: 165-173.

58. Alberto F, Bouffier L, Louvet JM, Lamy JB, Delzon S, Kremer A (2011) Adaptive responses for seed and leaf phenology in natural populations of sessile oak along an altitudinal gradient. Journal of Evolutionary Biology 24: 1442-1454. doi: 10.1111/j.1420-9101.2011.02277.x PMID: 21507119
59. Dantec CF, Ducasse H, Capdevielle X, Fabreguettes O, Delzon S, Desprez-Loustau M (2015) Escape of spring frost and disease through phenological variations in oak populations along elevation gradients. Journal of Ecology 103: 1044–1056.

60. Farrar JF (1992) Beyond photosynthesis: the translocation and respiration of diseased leaves. In: Ayres EPG, editor. Pests and pathogens. Plants responses to foliar attack. Oxford: BIOS Scientific Publishers. p. 107–128.

61. Hansen EM, Stone JK, Capitano BR, Rosso P, Sutton W, Kanaskie A, et al. (2000) Incidence and impact of Swiss needle cast in forest plantations of Douglas-fir in coastal Oregon. Plant Disease 84: 773–779.

62. Manter DK, Bond BJ, Kavanagh KL, Rosso PH, Filip GM (2000) Pseudothecia of the Swiss needle cast fungus, *Phaeocryptopus gaeumannii*, physically block stomata of Douglas-fir, reducing CO₂ assimilation. New Phytologist 148: 481–491.

63. Panstruga R, Dodds PN (2009) Terrific protein traffic: the mystery of effector protein delivery by filamentous plant pathogens. Science 324: 748–750. doi: 10.1126/science.1171652 PMID: 19423815

64. Oliva J, Stenlid J, Martinez-Vilalta J (2014) The effect of fungal pathogens on the water and carbon economy of trees: implications for drought-induced mortality. New Phytologist 203: 1028–1035. doi: 10.1111/nph.12857 PMID: 24824859

65. Newsham KK, Oxborough K, Greenslade PD, McLeod AR (2000) UV-B radiation constrains the photo-synthesis of *Quercus robur* through impacts on the abundance of *Microsphaera alphitoides*. Forest Pathology 30: 265–275.

66. Michelot A, Simard S, Rathgeber C, Dufrêne E, Damesin C (2012) Comparing the intra-annual wood formation of three European species (*Fagus sylvatica*, *Quercus petraea* and *Pinus sylvestris*) as related to leaf phenology and non-structural carbohydrate dynamics. Tree Physiology 32: 1033–1043. doi: 10.1093/treephys/tps052 PMID: 22718524

67. Zhang SY, Owoundi RE, Nepveu G, Mothe F, Dhote JF (1993) Modelling wood density in European oak (*Quercus petraea* and *Quercus robur*) and simulating the silvicultural influence. Canadian Journal of Forest Research 23: 2587–2593.

68. Li S, Daudin JJ, Piou D, Robinet C, Jactel H (2015) Periodicity and synchrony of pine processionary moth outbreaks in France. Forest Ecology and Management 354: 309–317.

69. Coakley SM (1988) Variation in climate and prediction of disease in plants. Annual Review of Phytopathology 26: 163–81.