Research Article |
Corresponding author: Alex Stemmelen ( alex.stemmelen@inrae.fr ) Academic editor: Manuela Branco
© 2023 Alex Stemmelen, Bastien Castagneyrol, Quentin Ponette, Simone Prospero, Gilles San Martin, Salome Schneider, Hervé Jactel.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Stemmelen A, Castagneyrol B, Ponette Q, Prospero S, San Martin G, Schneider S, Jactel H (2023) Tree diversity reduces co-infestation of Douglas fir by two exotic pests and pathogens. In: Jactel H, Orazio C, Robinet C, Douma JC, Santini A, Battisti A, Branco M, Seehausen L, Kenis M (Eds) Conceptual and technical innovations to better manage invasions of alien pests and pathogens in forests. NeoBiota 84: 397-413. https://doi.org/10.3897/neobiota.84.94109
|
The number of non-native invasive pests and pathogens has increased dramatically in recent years, with disastrous consequences for the health of forests worldwide. Multiple studies have shown that mixed forests may suffer less damage from insect pests than single species forest. This “associational resistance” can be notably explained by the fact that heterospecific neighbours make it more difficult for herbivores to locate and then exploit their host tree. However, the validity of these findings in the case of non-native, invasive pests and pathogens remains to be demonstrated. In this study, we monitored over two hundred Douglas firs in pure and mixed plots of a tree diversity experiment to assess the damage from the non-native gall midge Contarinia pseudotsugae and the non-native needle cast Nothophaeocryptopus gaeumannii. The probability of Swiss needle-cast infection was lower in Douglas fir trees surrounded by heterospecific neighbours. Gall midge damage was lower on Douglas firs surrounded by taller neighbours, consistent with the hypothesis of reduced host Apparency. Douglas fir trees that were more damaged by C. pseudotsugae were also more often infected by N. gaeumannii. Our study thus provides partial support of the associational resistance hypothesis of mixed forests against exotic pests and pathogens. Promoting forest species diversity at the stand level could, therefore, offer interesting prospects for reducing the impact of biological invasions, especially those involving both pests and pathogens.
Douglas fir, exotic pathogen, exotic pest, resource concentration, tree apparency, tree diversity
In the last decades, climate change and increasing global trade and travel have led to a dramatic increase in invasion by non-native forest insects and pathogens (
Studies from agricultural systems have shown that plant diversity can lead to a decrease in insect damage (
The mechanisms driving associational resistance effects are further complicated by overlooked horizontal interactions between herbivores and between herbivores and plant pathogens. For example, a meta-analysis by
Contarinia pseudotsugae Condrashoff (Cecidomyiidae) is an exotic invasive forest pest in western Europe, originating from North America where it causes damage on Douglas fir (Pseudotsuga menziesii Franco), a coniferous tree species also originating from North America (
One of the main pathogens frequently co-occurring with C. pseudotsugae on Douglas fir is the fungus Nothophaeocryptopus gaeumannii (T. Rohde) Videira, C. Nakash., which causes a foliar disease called Swiss needle cast. Although discovered in Switzerland in 1925, this pathogen is native to North America – hence its misleading common name – and is nowadays present in Europe and Australasia (
In this study, we assessed the effect of forest mixtures on damage caused to Douglas fir by the exotic pest C. pseudotsugae and the exotic pathogen N. gaeumannii, as well as the potential effect of C. pseudotsugae damage on the co-occurring infection by N. gaeumannii. We used a long-term tree diversity experiment to sample 207 Douglas firs of various height, located in plots of increasing tree diversity from monoculture to four species mixtures. We measured needle damage by the Douglas fir needle midge and infection by the Swiss needle cast to test the following hypotheses: (1) Needle damage by C. pseudotsugae and N. gaeumannii increases with increasing proportion of Douglas fir in forest plot (host concentration hypothesis); (2) Douglas firs taller than neighbouring trees suffer more needle damage than smaller ones (host Apparency hypothesis); (3) Damage by C. pseudotsugae increases the probability of infection by N. gaeumannii. In doing so, our study aimed at a better understanding of the ecological factors driving primary and subsequent invasion by exotic pests and pathogens in pure vs. mixed forests.
The study was conducted in Belgium, six years after the first detection of C. pseudotsugae in the country, in the tree diversity experiment FORBIO (
Map of Belgium with the location and experimental set-up of the FORBIO Gedinne sub-sites. The tree species diversity per plot ranged from one species (white) to four species (dark grey).
Tree species were planted at an equivalent proportion in mixtures (e.g. 50:50 in two-species mixtures). Even in mixture, trees were planted in monospecific patches of 3 × 3 trees, with patches arranged in a checkerboard pattern in the two-species mixtures and randomly distributed in the three- and four-species mixtures. Therefore, individual trees were surrounded by a varying number of conspecific neighbours, even in mixed plots in which the central tree of each monospecific patch was always surrounded by conspecific neighbours. Finally, four subplots were delimited in each plot and consisted in a 4 × 4 tree patches, where tree height and circumference are measured yearly.
In each of the 20 plots containing Douglas fir (10 plots in each sub-site), we randomly selected three of the four sub-plots in which tree height and circumference had been measured in the current year (2021). We sampled every Douglas fir present in those subplots, for a total of 207 trees. Finally, we calculated the mean height difference between each selected Douglas fir and their direct neighbours, regardless of whether they were other Douglas firs or trees of the other three species (henceforth called “Apparency”), setting tree height to zero when neighbours of sampled Douglas fir were missing (dead).
On each sampled Douglas fir, we selected a branch, at mid-height of the tree crown and collected five current-year shoots, starting with the terminal shoot, then picking up every second lateral shoot down to the base of the branch. Shoots were then placed in paper bags, sealed and sent to Bordeaux (France) for damage assessment by a single observer (AS), blind to treatment (plot) identity.
We estimated “needle damage” as the percentage of needles on a shoot that turned dark-brown because of the presence of Contarinia pseudotsugae galls, following the methodology developed by the Walloon Forest Health Observatory (OWSF). We used an ordinal scale of six percentage classes of damage: 0%; 1 – 10%; 11 – 20%; 21 – 40%; 41 – 60%; 61 – 90%; 91 – 100%, assigning each of the five shoots to one damage class. Finally, we estimated needle damage at the level of individual tree by averaging the median values of defoliation class across the five shoots. Five needles per trees were collected and used to detect the presence or absence of N. gaeumannii.
For DNA extraction, the collected needles were placed in 2ml Eppendorf Tubes containing a sterile metal bead (4mm diameter). Tubes were closed with an AirPore tape, stored at -20 °C and lyophilised in an Alpha 2-4 LD plus freeze dryer (Christ, Osterode am Harz, Germany) for 1.5 days. After lyophilisation, the needles were crushed in a MM400 Retsch Mill (Retsch, Haan, Germany) at 30 Hz for 2 min. Further steps in DNA extraction were done using the Sbeadex Plant kit (LGC Genomics GmbH, Berlin, Germany). Lysis was performed by adding 300 μl Lysis Buffer PVP (LGC Genomics GmbH) containing 40 μl/ml Debris capture beads (LGC Genomics GmbH), 0.3 μg/μl Proteinase K (LGC Genomics GmbH) and 1.3% 1‐thioglycol (Sigma‐Aldrich, Buchs, Switzerland) to each sample and subsequent mixing until the samples were well suspended. After incubation at 60 °C for 1 hr, samples were centrifuged at 5,700 g for 5 min and 200 μl of lysate added to a 96-well deep-well plate (Thermo Fisher Scientific, Waltham, USA) containing 420 μl Binding Buffer PN (LGC Genomics GmbH) and 10 μl Sbeadex Particles (LGC Genomics GmbH). The following steps were conducted using the automated KingFisherTM Flex Purification 96 System (Thermo Fisher Scientific). After mixing for 30 min, the Sbeadex Particles with the DNA attached were collected and transferred to 400 μl Wash Buffer PN1 (LGC Genomics GmbH) and mixed for 6 min. This step was repeated before the particles were transferred into 400 μl Wash Buffer PN2 (LGC Genomics GmbH). After mixing for 6 min, the particles were transferred to 50 μl Elution Buffer AMP (LGC Genomics GmbH) and mixed for 7 min to elute the DNA attached to the particles. All DNA purification steps were performed at room temperature. The DNA extracts were stored at -20 °C before dilution and further analyses.
To detect N. gaeumannii, we performed a multiplex quantitative real-time PCR (qPCR) with a specific primer pair targeting the beta tubulin gene and the corresponding dual-labelled probe (
To test the hypothesis that taller trees are more damaged by C. pseudotsugae or more exposed to infection by N. gaeumannii, we used Douglas fir absolute (Height) and relative (Apparency) height. We also included the density of Douglas fir in the plot and in the direct neighbourhood – number of Douglas fir trees present amongst the eight nearest neighbors – of a focal Douglas fir tree to investigate the potential validity of the resource concentration hypothesis. Density of Douglas fir corresponded to the number of Douglas fir either in the plot or in the direct neighbourhood of a focal Douglas fir. Finally, as damage by C. pseudotsugae could alter the biology of the needle or weaken Douglas fir prior an infection by the Swiss needle cast, we included C. pseudotsugae damage as an additional explanatory variable only in the model built to explain the probability of infection by N. gaeumannii.
We built two different models to test for the effect of explanatory variables on needle damage by C. pseudotsugae and infection by P. gaeumannii.
First, we used a linear mixed model to test the effect of Douglas fir absolute (Height) and relative (Apparency) height and Douglas fir density in the plot (Density plot) and in the direct neighbourhood (Density neigh.) as well as the interaction between height and Apparency on mean C. pseudotsugae damage per tree. Subplots (Subplot_ID), nested within Plot (Plot_ID), were included as random factors to account for spatial autocorrelation. Response variable was square-rooted to satisfy model assumptions of normality and homogeneity of residuals and predictors were scaled. Variance inflation factors were checked for every explanatory variables used in the model and were never > 5, the usual cut-off values used to check for multicollinearity issues (
Second, we used a generalised linear model with a binomial error distribution family to analyse sources of variation in the probability of Douglas fir infection by the Swiss needle cast. We used the same model structure as in Eq. 1, adding C. pseudotsugae. damage as an additional covariate. Again, the interaction between height and apparency had no significant effect and we removed it from the final model during the model simplification process. Statistical analyses were performed using R software version 4.1.0 (
The height of the focal Douglas fir studied ranged from 175 to 954 cm, with a mean (± SE) of 559.6 ± 13.0 cm. Mean height (± SE) of European beech, sycamore maple, sessile oak and hybrid larch were 357.8 ± 9.2, 364.0 ± 16.3, 407.6 ± 14.2 and 948.8 ± 13.6 cm, respectively (Supplementary material). Height difference between a Douglas fir and its closest neighbours was on average (± SE) -14.2 ± 13.8 cm, i.e. Douglas firs were on average shorter than their neighbouring trees. Maximum differences were recorded in plots 20 and 22 (Species composition: Beech, maple, larch and Douglas fir, in the Gribelle and Gouverneurs blocks, respectively), where Douglas fir was 400 cm shorter and 520 cm taller than their neighbours, respectively.
Mean needle damage by C. pseudotsugae was on average 16.25 ± 1.25% needles attacked, ranging from no damage at all on 12 Douglas firs to 91.5% of damage on the most impacted Douglas fir. Infection by the Swiss needle cast was confirmed on 118 (57%) of the 207 sampled Douglas firs.
Taller Douglas firs (absolute height) suffered significantly more damage by C. pseudotsugae than smaller ones (average standardised model coefficient parameter estimate ± standard error: 1.10 ± 0.24, n = 192, Fig.
Summary of model coefficients for the two models presented in this study. Bold characters indicate that parameters are significant. R2m and R2c represent R2 of fixed and fixed plus random factors, respectively.
Models / Parameter | Estimate | Standard error | 95% Cl | P value | R2m (R2c) |
---|---|---|---|---|---|
Model 1 – Mean needle damage | 0.121 (0.197) | ||||
(Intercept) | 3.47 | 0.19 | (3.10; 3.83) | < 0.001 | |
Apparency | 0.59 | 0.22 | (0.16; 1.04) | 0.008 | |
Height | 1.10 | 0.24 | (0.63; 1.60) | < 0.001 | |
Density neigh. | -0.01 | 0.15 | (-0.31; 0.29) | 0.919 | |
Density plot | 0.08 | 0.19 | (-0.29 0.45) | 0.668 | |
Model 2 – Infection rate | 0.196 (0.518) | ||||
(Intercept) | 0.29 | 0.36 | (-0.48; 1.09) | 0.420 | |
Apparency | -0.12 | 0.34 | (-0.83; 0.53) | 0.711 | |
Height | -0.54 | 0.42 | (-1.46; 0.26) | 0.201 | |
Density neigh. | 0.70 | 0.22 | (0.27; 1.17) | 0.002 | |
Density plot | 0.61 | 0.37 | (-0.15; 1.45) | 0.105 | |
Contarinia damage | 0.65 | 0.23 | (0.21; 1.16) | 0.006 |
Effect of tree height of a focal Douglas fir on mean damage by C. pseudotsugae. Green dots represent individual trees. Solid and dashed lines represent prediction and adjusted standard error of model 1, with values of other model variables set at average (Table
Effect of Douglas fir apparency on the mean damage by C. pseudotsugae. Apparency is the mean height difference between a focal Douglas fir and its direct neighbours. Dots at the right and at the left of the vertical dashed line are Douglas fir trees that are taller and smaller than their neighbours, respectively. Solid and dashed lines represent prediction and adjusted standard error of model 1, with values of other model parameters set at average (Table
The probability of presence of N. gaeumannii increased significantly with the density of Douglas firs at the neighbouring scale (0.70 ± 0.22, n = 192, Fig.
Effects of (A) the density of Douglas fir in the direct neighbourhood of a focal Douglas fir and (B) of mean damage by C. pseudotsugae on the probability of a Douglas fir being infected by Nothophaeocryptopus gaeumannii, the causal agent of Swiss needle cast. Solid and dashed lines represent prediction and adjusted standard error of model 2 (Table
Our study shows that tree height and Apparency were two main drivers of needle midge C. pseudotsugae on Douglas fir. Douglas firs taller than their neighbours exhibited more needle damage. We found no effect of Douglas fir concentration, nor at the plot level, nor at the direct neighbours’ scale, on damage by C. pseudotsugae. However, there was a significant, positive effect of Douglas concentration at the local scale on the probability of infection by the Swiss needle cast fungus N. gaeumannii. Douglas fir suffering more needle damage by C. pseudotsugae also had higher probability of being infected by the fungal pathogen.
Tree Apparency is a component of tree diversity. The concept was initially coined to describe the probability of a plant being identified by its herbivores (
In addition to the effect of tree Apparency, we also found a positive effect of absolute tree height, with taller Douglas firs suffering from more needle damage than smaller ones. This result confirms the outcome of a study conducted by
Unexpectedly, we did not find any effect of Douglas fir density on mean damage by C. pseudotsugae. Plots with more abundant resource for the invasive pest were not more infested than plots with a lower concentration of host trees, contradicting the resource concentration hypothesis. This might be due to the small size of the plots in the FORBIO experiment (42 m × 42 m) (
By contrast, the probability of a Douglas fir to be infected by N. gaeumannii was positively correlated with the density of Douglas fir in its direct neighbourhood. This result is consistent with recent studies that linked host tree proportion and pathogen transmission amongst neighbouring individuals (
Finally, we found that damage by C. pseudotsugae increased the probability of an infection by N. gaeumannii. Interactions between two or more invasive species in an area can result in invasion meltdown (
Our study provides indirect support of the associational resistance hypothesis of mixed forests against exotic pests and pathogens. In the case of Douglas-fir Swiss needle cast, the effect was mediated by host dilution amongst heterospecific, non-host tree neighbours. In the case of Douglas-fir midge, the resistance effect was due to a reduction in host tree Apparency by taller neighbours. It should be noted, however, that the reduction of infestations by these exotic organisms was not general in mixed stands and might mainly depend on the species composition of the mixtures. In particular, the need to incorporate tree species growing faster than Douglas-fir, such as hybrid larch, could be investigated as a mean of reducing pest damage. We also demonstrated for the first time a positive interaction between Swiss needle cast and Douglas fir needle midge. This underlines the necessity to develop a holistic approach to the problem of biological invasions in forests and, in particular, to look for generic prevention methods, amongst which the increase of forest diversity appears to be promising.