Research Article |
Corresponding author: Benno A. Augustinus ( benno.augustinus@gmail.com ) Academic editor: Graeme Bourdot
© 2020 Benno A. Augustinus, Suzanne T. E. Lommen, Silvia Fogliatto, Francesco Vidotto, Tessa Smith, David Horvath, Maira Bonini, Rodolfo F. Gentili, Sandra Citterio, Heinz Müller-Schärer, Urs Schaffner.
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:
Augustinus BA, Lommen STE, Fogliatto S, Vidotto F, Smith T, Horvath D, Bonini M, Gentili RF, Citterio S, Müller-Schärer H, Schaffner U (2020) In-season leaf damage by a biocontrol agent explains reproductive output of an invasive plant species. NeoBiota 55: 117-146. https://doi.org/10.3897/neobiota.55.46874
|
One of the biggest challenges in classical biological control of invasive weeds is predicting the likelihood of success. Ambrosia artemisiifolia, a North American plant species that has become invasive in Europe, causes economic losses due to health problems resulting from its huge amount of highly allergenic pollen and as a weed to agricultural crops resulting from high seed densities. Here we assessed whether the pollen and seed output of the annual A. artemisiifolia (at the end of the season) is related to in-season abundance of, or damage by, the accidentally introduced biological control agent Ophraella communa. We monitored the growth and leaf damage of individually labelled A. artemisiifolia plants at four locations in Northern Italy and recorded abundance of different O. communa life stages at regular intervals. We found that the in-season level of leaf damage by O. communa consistently helped to explain seed production in combination with plant volume and site throughout the season. Feeding damage, plant volume and site also explained pollen production by A. artemisiifolia six weeks before male flower formation. At three out of four sites, plants with more than 10% leaf damage in mid-June or early July had a very low likelihood of seed formation. Leaf damage proved to be a better explanatory variable than O. communa abundance. Our results suggest that the monitoring of the in-season leaf damage can help to project the local impact of O. communa on A. artemisiifolia at the end of the season and thus inform management regarding the needs for additional measures to control this prominent invader.
Ambrosia artemisiifolia, biological invasions, classical biological control, common ragweed, herbivory, Ophraella communa
It is now well established that plant species that are introduced into areas outside their native range and become invasive can wreak serious impact on nature and human well-being (
One of the biggest challenges in classical biological control of weeds is predicting the likelihood of success, and thus the necessity for considering additional management practices (
Predicting the impact of herbivore abundance on plants has a long history in crop pest forecasting (
Ambrosia artemisiifolia L. (Asteraceae) is one of the most notorious plant invaders in Europe (
In Europe, O. communa was found for the first time in Northern Italy in 2013, probably also due to an accidental introduction (
In Northern Italy, the yearly peak of O. communa population size is only reached at the time when the first flower buds are produced. Identifying earlier, in-season indicators that are related to the level of biological control at the end of the season could help to project whether in a particular season or location O. communa damages A. artemisiifolia to such an extent that it prevents plants from reproduction, i.e. from producing pollen (which impacts human health) or seeds (which impacts long-term population dynamics and crop yield).
Here we report on a field experiment to assess whether abundance of or damage by O. communa during the season is related to A. artemisiifolia reproduction at the end of the season. We followed individually labelled A. artemisiifolia plants in four locations in Northern Italy during the summer of 2016 to answer the following questions: (1) what is the in-season variation in a) in-season survival of A. artemisiifolia, b) the number of O. communa individuals of, and leaf damage caused by O. communa on individual A. artemisiifolia plants , and (2) what is the effect of in-season O. communa abundance or leaf damage on A. artemisiifolia reproduction at the end of the season?
Ambrosia artemisiifolia is an annual plant that has invaded areas in all continents except Antarctica (
Ophraella communa is a multivoltine leaf beetle which overwinters at the adult stage and lays eggs in egg batches in spring. The beetle then goes through three larval stages, which feed on the green parts of the host plant. It then pupates and starts mating shortly after emergence from the lightly woven cocoon. Adults feed on green parts of the plant as well. In Northern Italy, the beetle can complete up to 4 generations per year (
We selected three former crop fields and one meadow with natural populations of both A. artemisiifolia and O. communa in the Po Plain of the Italian Piedmont and Lombardy regions (see Suppl. material
The study plants were selected between 13 and 18 June 2016, when A. artemisiifolia was between the 4- and the 12-leaf stage. We maximised the variation in initial size of A. artemisiifolia by randomly measuring plants at each site for 10 minutes and separating them into three equally numbered size classes (small, medium, large). We then laid transects of 20 m length through the study plots and selected 20 plants per size class along this transect, with an as homogeneous distribution over the site as possible. Minimum distance between selected plants was 50 cm and the maximum distance away from the transect was 2 m. Plants were individually marked with an aluminium label around the stem and a bamboo stick.
We decided to start our experiment in mid-June to exclude background seedling mortality from the dataset, since seedling establishment can vary considerably within and among sites (
In order to increase intra-site variation in abundance of and damage by O. communa, two subplots of approximately 5 m long along the transect were selected at random for insecticide application. The two subplots contained in total 12 labelled plants (4 plants per size class) per site. These subplots were sprayed twice a month with insecticides, alternating between contact and systemic insecticides. We used Lambda-Cyhalothrin in a dosage of 20g/ha (Syngenta KarateZeon) as contact insecticide, and a combination of Acetamiprid in a dosage of 100g/ha (Sipkam EPIK), and Deltamethrin in a dosage of 20g/ha (Bayer DecisEVO) as systemic insecticides. Insecticides were applied at a spray volume of 1000L/ha using a backpack sprayer. Previous studies revealed that there is no direct effect of this insecticide treatment on the measured plant parameters (
Plant survival and size, O. communa abundance and leaf damage caused by O. communa were assessed on individual plants six times (“censuses”) at three-week intervals from mid-June until mid-September 2016 (see exact dates in Suppl. material
Volume = height * π * (width/4)2
To assess the abundance of O. communa on individual plants, we counted the number of O. communa egg batches, larvae >5 mm long (larger L2 and L3 larvae), and the number of adults on each labelled plant at each census. We disregarded egg batches with less than 5 eggs, because laboratory experiments indicated that eggs from small egg batches are mostly unfertilised (Augustinus, unpublished data). As small larvae are difficult to find since they can hide in buds and flowers, we did not count these to minimize observer errors. In addition, we measured leaf damage per plant by estimating the percent leaf area removed by O. communa from the total leaf area if the plant was intact (plants without leaves were given a value of 100% area removed). We did not score damage that was clearly not caused by O. communa (e.g. with traces of snail mucus). However, we never observed other leaf-chewing insect herbivores on A. artemisiifolia than O. communa, and rarely found traces of molluscs.
To estimate levels of plant competition early in the season, we assessed percent bare soil in a 50×50cm square around each marked plant in early July. A square frame of 50×50cm was laid around a plant and the fraction of that surface covered by bare soil, when projecting the vegetation onto the ground, was estimated by at least two persons and the average taken. Stones or dead leaf material were scored as bare soil as well. In late August, we measured the summed length of all racemes per plant as a proxy for pollen production (
To compare the change in leaf damage over time between sites, we conducted a repeated measures ANOVA with damage as response variable, site as fixed variable, and census as random effect. The fit of the residuals was evaluated graphically, and we took the square root of damage to obtain a better fit.
Because of the highly zero-inflated nature of our data, we applied a hurdle approach to analyse the effect of O. communa numbers on male (i.e. pollen) and female (i.e. seeds) A. artemisiifolia reproduction by first using presence/absence of racemes (pollen-bearing structures) in late August, and of seeds in mid-September to assess the probability of male and female reproduction, respectively. In a second step, we analysed the quantity of male and female reproduction conditional on the probability of reproduction (i.e. only using plants that did produce), using raceme length (as a proxy for the number of pollen produced), and numbers of seeds as response variables, respectively.
In the first part of the hurdle approach, we assessed the effect of O. communa on likelihood of raceme or seed formation in separate analyses by formulating generalised linear models for each of the first four (for raceme formation) or five (for seed formation) censuses. As fixed effects we included site, the natural logarithm of plant volume, as well as none or one of the four O. communa-related variables (number of eggs per plant, number of pupae per plant, number of adults per plant, and percent leaf damage) in each model, as these were inter-correlated. We produced models with and without an interaction term for the O. communa-related variable and plant volume, and with and without percent bare soil. We compared all resulting 18 models for each response variable at each census and selected the model with the lowest conditional Akaike information criterion (AICc) value, which penalizes models with more parameters (
Delta AICc values (upper number) compared to the model with the lowest AICc value, pseudo r-squared (second number), odds ratio for the O. communa related factor (third number), and confidence interval of the odds ratio for the O. communa related factor (lowest number) for models showing correlation between chance of successful raceme formation (left part of table) and total raceme length of raceme-producing plants (right part of table) and explanatory factors at different censuses. Models, where the confidence interval of the odds ratio for the effect size of the O. communa related factor does not cross 0, are shaded. Models including interactions with volume and the explanatory factor are marked with ‘*’. Corresponding p-values can be found in Suppl. material
Factor | Probability of raceme formation dependent on factor | Raceme length dependent on factor … | ||||||
---|---|---|---|---|---|---|---|---|
Mid-June | Early July | Late July | Early August | Mid-June | Early July | Late July | Early August | |
No O. communa parameter | 59 | 58 | 46 | 6.2 | 13 | 16 | 8.1 | 1.5 |
0.40 | 0.40 | 0.37 | 0.40 | 0.724 | 0.737 | 0.738 | 0.713 | |
# egg batches | 59 | 60 | 48 | 8.2 | 15 | 18 | 3.5 | 1.5 |
0.41 | 0.40 | 0.37 | 0.40 | 0.72 | 0.74 | 0.74 | 0.72 | |
0.74 | 0.88 | 0.92 | 1.03 | 1.11 | 1.15 | 1.55 | 1.06 | |
(0.48, 1.15) | (0.58, 1.34) | (0.52,1.65) | (0.88, 1.22) | (0.70, 1.75) | (0.86, 1.55) | (1.12, 2.13) | (0.98, 1.16) | |
# larvae | 61 | 60 | 47 | 7.8 | 14 | 18 | 10 | 2.7 |
0.40 | 0.40 | 0.37 | 0.40 | 0.73 | 0.74 | 0.74 | 0.72 | |
0.82 | 1.01 | 0.74 | 1.13 | 1.82 | 1.08 | 1.20 | 1.06 | |
(0.25, 2.72) | (0.78, 1.32) | (0.42, 1.30) | (0.81, 1.57) | (0.70, 4.73) | (0.91, 1.28) | (0.76, 1.91) | (0.95, 1.19) | |
# adults | 60 | 60 | 48 | 4.5* | 15 | 16 | 10 | 0 |
0.41 | 0.40 | 0.37 | 0.43 | 0.73 | 0.75 | 0.75 | 0.73 | |
0.85 | 1.05 | 0.95 | 24.61 | 1.13 | 1.25 | 1.09 | 1.13 | |
(0.58, 1.25) | (0.65, 1.70) | (0.73, 1.23) | (0.97, 624.27) | (0.72, 1.77) | (0.93, 1.67) | (0.91, 1.30) | (1.00, 1.28) | |
% leaf damage | 59 | 58 | 39* | 0* | 15 | 18 | 11 | 2.6 |
0.41 | 0.41 | 0.43 | 0.46 | 0.73 | 0.74 | 0.74 | 0.72 | |
0.97 | 0.98 | 1.02 | 1.06 | 1.01 | 0.99 | 1.01 | 0.98 | |
(0.94, 1.01) | (0.96, 1.01) | (1.00, 1.04) | (0.96, 1.17) | (0.97, 1.04) | (0.97, 1.02) | (0.98, 1.03) | (0.95, 1.01) |
Delta AICc values (upper number) compared to the model with the lowest AICc value, pseudo r-squared (second number), odds ratio for the O. communa related factor (third number), and confidence interval of the odds ratio for the O. communa related factor (lowest number) for models showing correlation between chance of successful seed formation (left part of table) and total seeds produced (right part of table) and explanatory factors at different censuses. Models, where the confidence interval of the odds ratio for the effect size of the O. communa related factor does not cross 0, are shaded. Models including interactions with volume and the explanatory factor are marked with ‘*’. Corresponding p-values can be found in Suppl. material
Factor | Probability of seed formation dependent on… | Number of seeds produced dependent on… | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
Mid-June | Early July | Late July | Early August | Late August | Mid-June | Early July | Late July | Early August | Late August | |
No O. communa parameter | 69 | 70 | 63 | 63 | 51 | 9.0 | 9.0 | 10 | 1.2 | 9.5 |
0.16 | 0.17 | 0.16 | 0.21 | 0.16 | 0.85 | 0.86 | 0.85 | 0.86 | 0.85 | |
# egg batches | 66* | 66* | 61 | 28* | 51* | 9.5 | 12 | 13 | 1.1 | 13 |
0.21 | 0.23 | 0.16 | 0.32 | 0.19 | 0.85 | 0.86 | 0.85 | 0.85 | 0.85 | |
1.79 | 0.00 | 1.30 | 0.03 | 0 | 1.19 | 1.06 | 1.16 | 0.91 | 1.04 | |
(0.83, 3.87) | (0, 4.51) | (0.76, 2.23) | (0.00, 0.44) | (0, inf) | (0.95, 1.49) | (0.85, 1.33) | (0.77, 1.73) | (0.82, 1.01) | (0.15, 6.99) | |
# larvae | 60* | 71 | 57 | 28* | 26* | 12 | 11 | 13 | 4.4 | 13 |
0.26 | 0.18 | 0.17 | 0.32 | 0.40 | 0.84 | 0.86 | 0.85 | 0.83 | 0.85 | |
0.00 | 1.17 | 1.31 | 0.00 | 0.00 | 1.08 | 1.08 | 1.05 | 1.07 | 1.00 | |
(0, inf) | (0.90, 1.54) | (0.79, 2.18) | (0, 0.79) | (0, 0.27) | (0.54, 2.16) | (0.92, 1.25) | (0.71, 1.55) | (0.68, 1.69) | (0.77, 1.30) | |
# adults | 68* | 64 | 58 | 28* | 31* | 12 | 11 | 10 | 3.9 | 12 |
0.19 | 0.23 | 0.17 | 0.31 | 0.36 | 0.84 | 0.86 | 0.87 | 0.84 | 0.85 | |
0.20 | 2.74 | 1.22 | 0 | 0.14 | 0.95 | 1.15 | 1.12 | 1.04 | 1.02 | |
(0, 109.11) | (1.15, 6.53) | (0.90, 1.65) | (0, 0.58) | (0.03, 0.74) | (0.75, 1.20) | (0.91, 1.45) | (0.97, 1.30) | (0.93, 1.16) | (0.95, 1.09) | |
% leaf damage | 56* | 59* | 48* | 1.9* | 0 | 12 | 6.8 | 12 | 3.2 | 0 |
0.29 | 0.27 | 0.30 | 0.51 | 0.57 | 0.83 | 0.88 | 0.86 | 0.84 | 0.90 | |
0.40 | 0.68 | 0.50 | 0.06 | 0.93 | 1.02 | 0.98 | 0.99 | 0.96 | 0.98 | |
(0.19, 0.83) | (0.51, 0.92) | (0.31, 0.80) | (0.01, 0.35) | (0.89, 0.98) | (0.87, 1.19) | (0.96, 0.99) | (0.97, 1.01) | (0.88, 1.04) | (0.97, 0.99) |
In the second step of the hurdle approach, we assessed the effect of O. communa on total raceme length or number of seeds of those plants that did produce racemes or seeds, respectively. We formulated a set of linear models for the natural logarithm of raceme length and number of seeds, assuming a Gaussian distribution of the response variable. We chose to use a Gaussian distribution over a Poisson distribution since it reduced AICc values of the fitted models by more than 5000 for every case. As fixed effects we included the natural logarithm of volume and site and added none or one of the four O. communa-related variables. To prevent overparameterization, we did not include bare soil and interactions with plant volume in these models, since the sample size of plants that successfully formed racemes and seeds was too low to include more than three fixed effects.
All analyses were conducted in R version 3.5.1 (2018–07–02) --”Feather Spray” (2018). Data were prepared using the readxl (
Until late July, we found less than one egg batch, larva or adult of O. communa per plant (Fig.
Average number of O. communa individuals per plant during the experiment in the four different sites. Different life stages are marked with different lines and symbols. Vertical lines indicate the standard error.
Plant volume steadily increased until late August, and decreased or stayed stable thereafter (Suppl. material
Violin plot of Ambrosia artemisiifolia leaf damage by O. communa feeding. The lines indicate the mean of the leaf damage scored on living plants in the different sites. The distribution of the damage measurements is shown with the grey shapes. Only damage of plants which were not treated with insecticides are displayed.
Models with the lowest delta AICc values (compared to the best performing model) for successful raceme formation included O. communa abundance parameters measured in early August (number of adults), and models with the lowest delta AICc values for raceme length of the plants that successfully formed racemes included O. communa abundance parameters measured in late July (number of egg batches) and early August (number of adults; see Table
In the model with the lowest delta AICc value for successful raceme formation in late July, we found a positive relationship between leaf damage in percent and successful raceme formation (Fig.
The selected models for successful seed formation included O. communa abundance parameters measured in early July (number of adults per plant), early August (number of adults, larvae and egg batches per plant) and late August (number of adults and larvae per plant) (Table
Odds-ratios of effect size of explanatory variables of the models with the lowest AIC per census, explaining successful raceme formation. Red dots/values <1 indicate that the effect is negative, blue dots/values >1 indicate that the effect is positive. The factor “site” with the corresponding site name in square brackets show the effect size of site compared to Busto. Plant volume “vol” (in cm3) is log-transformed for the analysis, leaf damage in percent is abbreviated with “dam”. In models with interaction between leaf damage in percent and volume, the effect size of this factor is described as “vol[log]*dam”.
Effect size of explanatory variables of the models with the lowest AIC per census, explaining successful seed formation. Red dots/values <1 indicate that the effect is negative, blue dots/values >1 indicate that the effect is positive. The factor “site” with the corresponding site name in square brackets shows the effect size of site compared to Busto. Plant volume “vol” (in cm3) is log-transformed for the analysis, leaf damage in percent is abbreviated with “dam”. In models with interaction between leaf damage in percent and volume, the effect size of this factor is described as “vol[log]*dam”.
Including percent leaf damage by O. communa in models for successful raceme formation generated the models with the lowest AICc values for late July and early August, and including percent leaf damage by O. communa measured in late July generated the respective model with the lowest AICc value for raceme length of plants that successfully formed racemes (see Table
Likelihood of successful raceme (A, C) and seed formation (B, D) dependent on plant volume. In A and B the data are presented for four damage classes and in C and D for three O. communa infestation classes (0,1, or ≥2 adults per plant). The large symbols give median volume and mean probability of raceme or seed formation, respectively, together with their associated standard errors. “Only volume” shows the average values without consideration of damage or abundance classes.
With regard to the models for successful seed formation, all models with the lowest AICc values calculated from mid-June to late August included percent leaf damage by O. communa. Of the A. artemisiifolia plants with more than 10% leaf damage in mid-June (19.7% of all plants), none produced seeds at the end of the season (Table
In general, models for the successful formation of racemes and seeds that included percent leaf damage had lower AICc values than those that included O. communa abundance parameters (Tables
In 18 out of 20 cases, the models including an interaction of plant volume and O. communa abundance or damage improved the model fit for successful seed formation, and in the two cases where O. communa abundance or damage improved the models for successful raceme formation, the model included an interaction of abundance or damage with plant volume. To explore the nature of these interaction terms, we displayed the interactions graphically, splitting the data into groups (by level of damage or abundance) and plotted the probability of successful raceme formation against the log of plant volume (Fig.
The models with the lowest AICc values within one census also had the highest pseudo R-squared values. Ophraella communa presence and/or damage explained the likelihood of seed formation better than the likelihood of raceme formation. Including leaf damage increased the pseudo-R2 value of models for successful seed formation much more (max. 40%) than for raceme formation (max. 6%). In contrast, including O. communa abundance parameters hardly improved the pseudo-R2 value of models of seed numbers (max. 4%) or raceme length (max. 2%). For probability of both raceme and seed formation, models had much lower AICc values and higher pseudo R-squared values from early August on. In general, O. communa induced leaf damage and abundance explained more variation the closer it was assessed to the flowering time.
Number of plants with a certain % leaf damage producing seeds at the end of the season. Given are the number of plats within a certain damage category producing seeds / total number of plants within this damage category.
% Damage | Mid-June | Early July | Late July | Early August | Late August | Late September |
---|---|---|---|---|---|---|
0 | 12/83 | 8/59 | 14/66 | 14/35 | 11/18 | 9/14 |
1–10% | 22/92 | 22/103 | 12/76 | 19/90 | 7/11 | 7/9 |
11–20% | 0/23 | 4/29 | 6/20 | 2/22 | 0/2 | 4/8 |
21–30% | 0/8 | 0/7 | 0/7 | 0/4 | 2/2 | 2/4 |
31–40% | 0/4 | 0/3 | 2/3 | 0/4 | 3/3 | 0/3 |
41–50% | 0/1 | 0/2 | 0/2 | 0/3 | 1/3 | 0/1 |
51–60% | 0/3 | 0/4 | 0/3 | 0/4 | 2/3 | 0/1 |
61–70% | 0/3 | 0/0 | 0/0 | 0/3 | 2/4 | 1/4 |
71–80% | 0/1 | 0/1 | 0/1 | 0/1 | 4/5 | 0/6 |
81–90% | 0/0 | 0/1 | 0/2 | 0/2 | 2/17 | 1/5 |
91–100% | 0/0 | 0/4 | 0/3 | 0/3 | 1/89 | 5/61 |
Our study provides evidence that the level of in-season leaf damage by O. communa, in combination with plant volume and site, helps to explain final seed production. Six weeks before flowering, leaf damage by O. communa together with plant volume is correlated to pollen production by A. artemisiifolia at the end of the season. Explanatory power of models improved over the season. Models including leaf damage had generally higher explanatory power than models including O. communa abundance parameters. For successful raceme formation, experimental sites had a much higher explanatory power than leaf damage, but for seed formation, explanatory power of leaf damage was similar to explanatory power of site, with lower variation. This offers possibilities to use in-season leaf damage for developing impact forecast models, which help informing management whether biological control is likely to successfully reduce seed production of this invasive alien plant species in a given region or year, or whether complementary management interventions should be considered to achieve long-term population decrease.
The peak in O. communa abundance in early August coincides with the expected timing of the fourth and last generation in this region (
While damage increased significantly in August at all sites, there was considerable variation in average leaf damage among sites (Table
It should be noted that our study did not cover the very first months of the growing season of A. artemisiifolia. In Northern Italy, gravid O. communa females that have overwintered start laying eggs on A. artemisiifolia seedlings as soon as the temperature is high enough for the beetle to fly (
Significant impact on target weed populations is only expected with high densities of biological control agents (
With regard to the probability of both raceme and seed formation, O. communa leaf damage appears to be a better explanatory variable than O. communa abundance, since AICc values were lower for models including damage than those including abundance for all census dates. This could be due to the behaviour of the beetle; Ophraella communa adults are highly mobile (
All but one model in which O. communa abundance explained the probability of reproductive organ formation contained an interaction with plant volume. Plant volume influences the response of the plant to abundance of herbivores or herbivore-induced leaf damage. In line with
In general, we found more adults on bigger plants, probably explained by a positive effect of plant volume on adult beetle abundance, rather than a positive effect of O. communa abundance on plant size. Caged experiments with varying plant sizes and number of adults could shed some additional light on the potentially interacting effects of plant volume and O. communa impact.
The models for the likelihood of seed formation generally had higher pseudo R-squared values than the models for the likelihood of raceme formation (Tables
Our study provides evidence that the window of impact by O. communa on reproductive output of A. artemisiifolia is relatively narrow (see Fig.
Arthropod demography is strongly influenced by climate, especially temperature, where an increase often results in quicker population growth. Since overall damage is strongly dependent on the number of generations, and as these are expected to increase with temperature in species with a multivoltine life cycles, damage is also expected to increase in a warming climate in the future (
Our findings that average leaf damage from mid-June onwards explained a significant amount of variation in the likelihood of seed formation indicates that O. communa feeding has a direct detrimental effect on female reproduction in A. artemisiifolia. Moreover, while the negative effect of leaf damage on the likelihood of pollen production only was significant in the census made in late July, O. communa exclusion experiments conducted in the same area revealed that O. communa reduces pollen production per unit area by 82% (Lommen et al. unpublished results). These findings are in line with an observed 80% decrease in airborne ragweed pollen counts in the Milano region since the establishment of O. communa (
Hence, our findings suggest that percent leaf damage in mid-June or early July could be used as an indicator for the likelihood that O. communa significantly reduces reproductive output of A. artemisiifolia at the end of the season (see Suppl. material
This study provides evidence that the level of in-season leaf damage by O. communa helps to explain the impact of this biological control agent on seed and – to a lesser extent – pollen production by A. artemisiifolia at the end of the season. Leaf damage measured as early as mid-June partially explains, in combination with plant volume, the likelihood of reproductive output of A. artemisiifolia at the end of the season. For example, none of the plants with more than 10% leaf damage in mid-June formed seeds at the end of the season. It should be noted, though, that at extreme sites where A. artemisiifolia plants grow 2 m and taller (such as at Busto Arsizio), impact of O. communa may be largely explained by plant volume, rather than by average leaf damage in early summer. Our results suggest that in-season assessment of leaf damage and plant volume could be used to develop predictive models for O. communa impact on A. artemisiifolia seed production, similar to the approach used in crop pest forecasting.
We are grateful to the land owner Andrea Airoldi, the communes of Magenta and Magnago, and the University of Torino for allowing us to conduct experiments on their properties. Rudolf Rohr and Peter Stoll provided valuable contributions to the statistical analyses. We thank Yvonne Buckley for reviewing the manuscript, which led to considerable improvements. The study was funded by the EU COST Action FA1203 ‘Sustainable management of Ambrosia artemisiifolia in Europe (SMARTER)’, the Swiss State Secretariat for Education, Research and Innovation (#C14.0063 to U.S. and #C13.0146 to H.M.S.), the Swiss National Science Foundation (#31003A_166448 to H.M.S), and the Pool de Recherche of the University of Fribourg (to H.M.S. and B.A.). Suzanne Lommen has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 786624. Urs Schaffner was supported by CABI with core financial support from its member countries (see http://www.cabi.org/about-cabi/who-we-work-with/key-donors/).
Study sites
Data type: species data
Census dates
Data type: measurement
p-values for O. communa
Data type: statistical data
Figure S1. Mean plant volume ± se of A. artemisiifolia plants measured during the experiment in the four experimental sites
Data type: statistical data
Figure S2
Explanation note: Likelihood of A. artemisiifolia seed formation dependent on O. communa leaf damage in early July. The different line types show the different responses between the sites.
Successful raceme formation
Explanation note: Summaries of selected glms, with successful raceme formation depending on different Ophraella communa abundance parameters, or leaf damage (in percent) inflicted by O. communa, per census.
Raceme length
Explanation note: Summaries of selected linear models, with raceme length depending on O. communa abundance parameters, or leaf damage (in percent) inflicted by O. communa.
Damage ~ abundance
Explanation note: Summaries of selected glms, with leaf damage depending on Ophraella communa abundance parameters.
Successful seed formation
Explanation note: Summaries of selected glm(m)s, with successful seed formation depending on Ophraella communa abundance parameters, or leaf damage (in percent) inflicted by O. communa.
Number of seeds
Explanation note: Summaries of selected lms, with number of seeds produced dependent on Ophraella communa abundance parameters, or leaf damage (in percent) inflicted by O. communa.