Research Article |
Corresponding author: Miloň Dvořák ( milon.dvorak@seznam.cz ) Academic editor: Andrea Battisti
© 2023 Miloň Dvořák, Petr Štoidl, Michal Rost.
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:
Dvořák M, Štoidl P, Rost M (2023) Vertical spread of Hymenoscyphus fraxineus propagules. 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: 231-246. https://doi.org/10.3897/neobiota.84.90981
|
Currently, the ash dieback causal agent Hymenoscyphus fraxineus is an established invasive pathogen in most European countries. Its potential to spread quickly among invaded forests is based on its propagules: airborne inoculum composed mainly of ascospores originated in apothecia growing on leaf litter infected during the previous vegetation season. The spread of the inoculum by air masses to distant areas is probable and depends on the availability of the ascospores in higher levels of air. Our study aimed to detect the inoculum in an infected area at heights of more than 20 meters. Our study was conducted in a municipal locality (Boršov nad Vltavou) with tens of infected ash trees (Fraxinus excelsior) in South Bohemia (SW Czechia). The infected trees surround an agricultural silo where five rotating arm spore traps (rotorods) were mounted for ten consequent 48h samplings during the peak of the sporulating season (17th July to 6th August 2020). The spore traps were mounted 48, 37, 25, 14 and 0,3 meters above ground. Samples were quantified by qPCR. Results clearly proved the ability of the spores to reach a height of 48 meters. Furthermore, H. fraxineus DNA was detected from all five spore traps during all ten samplings. Mostly, the amount of detected spores showed a decreasing trend with height, and varied a lot. During some of the samplings, higher spore concetrations were achieved at the top than at the lower traps, which can be explained by horizontal air transfer of the inoculum from other infected areas. Based on GLM analyses, higher spore concentrations were achieved during days without rain, lower air temperatures, after cloudy, humid and rainy weather without strong winds. A combination of rotorod ROTTRAP 52 with qPCR quantification proved to be an efficient technology for a study focused on the vertical spread of H. fraxineus propagules.
airborne inoculum, ash dieback, riseability, rotorod, ROTTRAP, spore trap
Hymenoscyphus fraxineus (T. Kowalski) Baral, Queloz, and Hosoya is the causal agent of the ash dieback (
The symptoms of the ash dieback and mechanism of infection of the host trees have been well described in numerous publications (
The predominant way to spread this pathogen is via the airborne ascospores (
The probability of the long distance spread of the ascospores depends on the height they can reach to be blown with the air masses (
The aim of our study was i) to test a methodological approach for describing the vertical spore dispersal pattern of H. fraxineus, and ii) to prove the presence of the inoculum as high as possible, at such a height where the spread of the aerosols is more likely affected by horizontal movement of air masses rather than by convection from the ground surface.
The sampling point was located in an agricultural silo in Boršov nad Vltavou (South Bohemia; SW Czech Republic); GPS: 48.9244°N, 14.4414°E; 412 m a. s. l. It is an industrial area of a village adjacent to Moldava River. It is located on the SE edge of a large plain called České Budějovice Basin (Českobudějovická pánev), where foothills of Šumava mountains called Blanský les start to rise. Due to the surrounding landscape being characterised by plains (especially in NW, N and E direction), the 52 m high silo of Boršov is probably exposed to wind currents blowing from areas that are at least tens of kilometres away. Trees in this area are mostly represented by mixtures of deciduous middle-European species formed at the riverbanks, gardens and parks, or alleys along railways and roads. Most of the groups of trees comprise a significant portion of ashes (Fraxinus excelsior) with typical symptoms of ash dieback (Fig.
Sampling locality and sampling tools A the sampling point in the focus of ash dieback. Agricultural silo is surrounded by ashes infected with H. fraxineus B sampling spots in the windows of the silo. ROTTRAPs 52 are installed in different heights R2 - R5 in the windows (red circles) of the silo on its NNW side C rotating arm spore trap ROTTRAP 52 installed at R5.
To sample the air inoculum rotating arm spore traps (rotorods) ROTTRAP 52 (Miloň Dvořák, Boršov nad Vltavou, Czech Republic) were employed (Fig.
The rotorods were installed at five different heights (Fig.
Meteorological data were partly measured by automatic meteorological station Signalizátor (AMET, Velké Bílovice, Czech Republic) and partly received from the archive of the Czech Hydrometeorological Institute (CHMI), station České Budějovice – Rožnov. From the sampling point the CHMI meteorological station is positioned 3.6 km to NE. Data taken from Signalizátor were daily means of relative air humidity (further in the text only air humidity). Data received from CHMI were: i) daily mean of air temperature (air temperature); ii) daily duration of sunshine (sunshine); iii) daily mean of air pressure (air pressure); iv) daily amount of precipitation (precipitation) and v) daily mean of wind speed (wind speed).
The genomic DNA from samples was extracted with a DNeasy plant minikit (Düsseldorf, Germany). Each microtube with exposed tape was supplemented with one 3-mm sterile tungsten bead and 20 pcs 2-mm glass beads, 400 µl of AP1 buffer and 4 µl of RNase. This mix was ground twice for 60 seconds using a high speed homogenizer Millmix 20 (Domel, d.o.o., Železniki, Slovenia) set at 30 Hz and incubated for 10 minutes at 65 °C. The microtubes were inverted three times during the incubation. Further steps were following the manufacturer’s instructions; however, the last step (elution) was not repeated to obtain higher concentration of DNA. DNA samples were eluted in 100 µl and stored at -20 °C before further processing.
Direct specific qPCR was performed using a QuantStudio 6 Flex Real-Time PCR System (Life Technologies Holdings Pte. Ltd., Singapore), Light Cycler 480 Probes master (Roche Diagnostics Nederland BV, Almere, the Netherlands) and primers and probes specific to H. fraxineus (
The concentrations of H. fraxineus DNA in the samples were expressed as numbers of copies of the target sequence in 1 µl of template DNA (further only CN). These values were obtained using a standard curve generated from reactions with different CNs (2.5×102 to 2.5×10-2) of plasmid pCR 2.1 TOPO TA vector (Invitrogen, Carlsbad, California, USA) by QuantStudioTM Real-Time PCR System Version 1.3 (Thermo Fisher Scientific). Plasmids contained species - specific insert (PCR products amplified with Cf-F and Cf-S primers). DNA was extracted from pure cultures of H. fraxineus (collection of Mendel University in Brno).
To express the absolute amount of ascospores in every analysed sample, an absolute quantification of ascospore suspension was performed. For that purpose ten apothecia were collected from ash leave rachises and immersed in 1 ml of distilled water in a 2-ml microtube. The following day, the clear liquid upper part without apothecia and other debris was transferred into a clean 2 ml microtube. Vortexed ascospore suspension was quantified in Bürker chamber by microscope. Ten-fold serial dilutions from 18750 to 1.875 ascospores in 100 µl of distilled water were transferred into clean 2 ml microtubes and DNA was extracted with the same protocol as for the spore trap samples. Extracted DNA samples were used as standards for a qPCR absolute quantification of the plasmids previously used for the quantification of the environmental samples. The lowest detectable concentration which turned positive in all three technical repetitions of the sample was 18.75 ascospores per sample (Ct = 36.328, SD = 0.862).
To describe the influence of meteorological factors on the ability of the inoculum to spread vertically, meteorological variables were averaged for three particular days of every sampling. Furthermore, the factor of riseability (FR) has been defined. It is expressed as the ratio of the ascospore concentration recorded at the highest sampler (R5) to the concentration at the lowest sampler (R1). It takes values lower than 1.0 in case the R5 ascospore concentration is lower than R1 concentration.
In order to describe the relationship between CN and the character of the weather described by explanatory variables (air temperature, precipitation, sunshine, air humidity, wind speed and air pressure), generalised linear regression models were constructed with explanatory variables measured on the same day as the dependent variable CN, or with explanatory variables recorded during previous four samplings (i.e. one sampling lag = period of preceding two days before the sampling started, two samplings lag = two to four days before, three samplings lag = four to six days before and four samplings lag = six to eight days before sampling) to simulate the lag of the pathogen’s reaction on the changing weather. Due to the nature of the dependent variable (a strictly positive variable showing a positive skew for R1, R5 and FR), we used the gamma distribution with a logarithmic link function when fitting the regression models. To select variables for individual models, we used the procedure described in
At the same time, we tried to model the relationship between CN and the height of the sampler. Due to convergence problems with an exponential model based on a differential equation, we used a somewhat simpler empirical exponential model with the following analytical form:
CN = β0⋅exp(-β1⋅h)
where β0 and β1 are the estimated regression coefficients, and h is the height in meters. All numerical calculations were performed using the R 4.2.0 programming environment (
Fresh apothecia were observed on infected rachises at the sampling locality from the beginning of July until the end of the sampling (6th August, 2020). Consequently, all ten 48-h samplings showed presence of H. fraxineus in the air (Fig.
Every sampler positively detected inoculum during every sampling. The lowest positively detected spore concentration was detected in samples from the highest sampler R5 from the sampling started on 4th August (CN = 0.018, Ct = 37.428). This least concentrated sample contained 1.89 ascospores; however, this amount was calculated extrapolating the ascospore suspension standard curve, hence the exact amount cannot be taken in consideration. The highest concentration was recorded at R1 during the following sampling started on 6th August (CN = 659.063, Ct = 27.181). Taking in account the sampling rate 52 l/min and sampling period of 48 h, we detected average spore concentration in a range from 0.013 to 462.12 spores per m3.
Most of the detected spore concentrations showed a clear decreasing trend following the height of the sampling point. The resulting nonlinear regression model (Fig.
CN = 150.3243⋅exp(-0.2155⋅h)
GLM analyses resulted in 10 significant models to estimate CNs from meteorological variables measured during the sampling and lagged of 2–4, 4–6 and 6–8 days. Their parameters can be found in supplementary file “Parameters of GLM models”. The calculated p and positive or negative meaning of each parameter are displayed in Table
Meteorological variables as parameters of GLM models and their p during sampling and 0–8 days before sampling (d.b.s.).
R1 | 6–8 d.b.s. | 4–6 d.b.s. | 2–4 d.b.s. | 0–2 d.b.s. | sampling |
Air temperature | .014 | .026 | |||
Sunshine | .009 | .034 | |||
Precipitation | .032 | .002 | |||
Air humidity | .032 | .045 | |||
Wind speed | .01 | .031 | |||
Air pressure | .012 | ||||
R5 | 6–8 d.b.s. | 4–6 d.b.s. | 2–4 d.b.s. | 0–2 d.b.s. | sampling |
Air temperature | .043 | .040 | .018 | ||
Sunshine | .001 | ||||
Precipitation | .001 | .001 | .002 | ||
Air humidity | .006 | .004 | |||
Wind speed | |||||
Air pressure | |||||
FR | 6–8 d.b.s. | 4–6 d.b.s. | 2–4 d.b.s. | 0–2 d.b.s. | sampling |
Air temperature | .016 | ||||
Sunshine | .009 | .029 | |||
Precipitation | .034 | .018 | .009 | .048 | |
Air humidity | .011 | ||||
Wind speed | .016 | .015 | |||
Air pressure | |||||
Legend: | Influence | ||||
positive | negative | ||||
Significant parameter (p < .05) | p | P | |||
Highly significant parameter (p < .01) | p | P |
Our results prove for the first time, that propagules of H. fraxineus are reliably detectable at almost 50 meters above ground, where they have their main source in infected rachises (
The results of this study also confirm statements of other authors (
However, not in every sampling did the CN values descend with increasing height. The sampling spot R4 showed lower values, than R5 for samplings 5, 6 and 7. This was probably due to overloading the sampling rods with dust, which is a critical handicap of rotorods (
Furthermore, the sampling site R3 gave higher CNs than it could be expected to follow the descendent trend during samplings 8 and 9. Similarly, CNs at R5 from samplings 1, 5 and 6 were also not lower than CN of lower sampling points. An explanation for this anomaly could be the effect of horizontal air currents which could bring a higher amount of inoculum from more distant localities, influencing the results of sampling only in these cases of low local concentrations. Generally, abnormalities in the trend of decreasing CN with the height of sampling can be interpreted as a consequence of the fact that inoculum detected in heights above 10 m might be representative for areas within a perimeter of at least tens of kilometres (
Long distance transfers of H. fraxineus by air masses have always been an important issue. It was considered that H. fraxineus had been introduced into Great Britain between 2008 and 2011 via long distance transfer of the air inoculum from mainland Europe. This statement was strongly supported with a model (
Results of the GLM modelling of the determination of the detected spore concentrations by meteorological factors discovered numerous relations. Although the low number of repetitions decreases the reliability of the results’ interpretation, we would like to highlight some of them:
Air temperature proved to have significantly positive effect two days preceding the sampling at both R1 and R5. This confirms previous observations of
The daily amount of precipitation showed different effects at different sampling heights. R1 was significantly positively affected by precipitations two to four days before sampling. Through the enhanced humidity of the ground surface, rain proved to be essential for successful ascospore maturation and release (
Daily duration of sunshine showed highly significantly negative effect on spore concentration at R1 and R5 two days before the sampling and on riseability during the preceding four to six days. However, the sunshine was significantly positive during the sampling at R1 and R5 and during the preceding two to four days it was significantly positively affecting the riseability. This partly confirms and partly neglects results of
Air humidity was found to be an important factor through our sampling. At R5 the humid air was significantly determining the inoculum load after two to six days; similarly, at R1 there was a significantly positive effect on inoculum concentrations lagged by two to four days. The promoting effect of air humidity on the disease spread and establishment has already been emphasised many times (
Wind speed proved to have a significantly negative influence on ascospore concentrations at R1 after two to six days. This effect is probably due to the desiccation of leaf rachises and growing apothecia which are not able to mature and sporulate under such conditions (
Air pressure did not show much importance in our experiment. It was a significant parameter with negative influence on spore concentrations at R1 after two to four days. At the same period and sampling height two other factors had positive influence on the spore concentrations: precipitations and air humidity. Naturally, rainy and humid weather is characterized by low air pressure, which we assume to be a reason for this result.
From a methodological point of view, rotating arm spore trap ROTTRAP 52 proved to be a reliable tool for the detection of H. fraxineus inoculum. All samplers successfully completed all ten 48h samplings at five sampling spots without any blackout even in hot or rainy weather. This reliability has been improved compared to previous experiments (
Our study revealed the permanent presence of the H. fraxineus inoculum up to 48 meters above the ground during the whole sampling period. Its concentration is continuously changing depending on previous and current weather, and decreases with height. It poses a persistent threat to ash trees, either at local or landscape scale. This finding supports a sceptical outlook for the future of ashes in European forests, but also confirms the important role of high height air sampling of the propagules of this invasive alien pathogen to ensure its reliable monitoring.
This work was funded by a HORIZON 2020 project (agreement no. 771271): Holistic Management of Emerging Forest Pests and Diseases (HOMED). We would like to thank to EUROPASTA SE company for cooperation and full access to their silo in Boršov nad Vltavou, namely to the director Ing. Josef Maška, Jaroslav Mrkvička, Ladislav Hakl and Jaroslav Hnilička. Furthermore, we thank Dr. Leticia Botella Sánchez for preparing the plasmid for absolute quantification.