Introduction

Common oak Quercus robur L. and sessile oak Q. petraea (Matt.) Liebl. are major oaks species in central and eastern European forests. Oak stands occupy about 6% of the total forest area of Poland. Since the early 1900’s, records of powdery mildews affecting European oaks have been mentioned since they are important fungal diseases (Mougou et al. 2008). In Europe, powdery mildew infections of common oak are caused by Erysiphe alphitoides (Braun and Cook 2012) (syn. Microsphaera alphitoides Griffon & Maubl.), which is the most prevalent and damaging foliar fungal species in oak stands, and along with primary insect attacks, they cause periodic dips in oak health (Topalidou and Shaw 2016). Common oak is one of the most susceptible species to E. alphitoides (Desprez-Loustau et al. 2011), and this pathogen is often mentioned as an inciting factor in European oak decline (Keča et al. 2016). Powdery mildew has particularly damaging effects on young seedlings, in nurseries, plantations, and can even damage mature trees under conditions favorable to disease (Marçais and Desprez-Loustau 2014).

Fungi from the genus Trichoderma are commonly found in all climatic zones. The most typical habitats of these fungi include soil and rotting wood (Druzhinina et al. 2006). Trichoderma spp. have been reported to exhibit antagonistic effects against at least 29 species of pathogenic fungi representing 18 genera such as Armillaria, Chondrostereum, Phytophthora, Rhizoctonia, Sclerotinia, and Verticillium (Druzhinina et al. 2011; Contreras-Cornejo et al. 2016). Species of Trichoderma have long been recognized as biological control agents (BCA) of plant diseases (Contreras-Cornejo et al. 2011; Hermosa et al. 2012; Zhang et al. 2016). Isolate T-39 of T. harzianum was the first commercialized fungal BCA (Trichodex®) for plant diseases (Tronsmo and Dennis 1977). Trichoderma sp. IZR D-11 is the basis of Fungilex L (Belarus, State Registration № 06–0078 dating from 28 November 2013) which is used for reduction of Botrytis cinerea on cucumber, parsley, dill, lettuce, but the species name is not known.

The mode of action of Trichoderma spp. as biocontrol agents can be direct by parasitizing the target organisms, or indirect through competition, modifications of the environment, or promoting plant defense mechanisms (Benitez et al. 2004). Trichoderma harzianum T39 applied to roots of greenhouse cucumber was able to reduce powdery mildew (caused by Sphaerotheca fusca) by 75–90% (Elad et al. 1998), and this was attributed to induced resistance. Elad (2000) considered that induced resistance was the main mode of action of T. harzianum against powdery mildew fungi. Another species, T. asperellum, has been used as a biocontrol agent to control Pythium myriotylum, which causes root rot in cocoyam (Mbarga et al. 2012), and Sclerotium sclerotiorum, which causes stem rot in carnation (Vinodkumar et al. 2017). Trichoderma asperellum reduced conidial germination of Leveillula taurica, causing powdery mildew on pepper (López et al. 2019).

Volatile organic compounds (VOCs) produced by plant growth-promoting fungi such as species of Phoma, Cladosporium or Ampelomyces significantly reduced fungal and bacterial disease severity in Arabidopsis by eliciting induced resistance against pathogens (Naznin et al. 2014). Volatiles have been acknowledged for their role in promoting plant growth and disease resistance against pathogens (Thissera et al. 2020; Cheffi et al. 2019; Ben et al. 2018; Belbahri et al. 2017; Mefteh et al. 2017). Trichoderma spp. are also known to have plant-growth promoting activities (Benitez et al. 2004), and volatiles produced by Trichoderma spp. are also known to induce resistance against plant pathogens (Lee et al. 2016). Gas Chromatography-Mass Spectrometry (GC-MS) analysis of VOCs and Liquid Chromatography High Resolution Mass Spectrometry (LC-HRMS) analysis of secondary metabolites have been effective in the discovery of numerous molecules with plant growth promoting abilities such as auxin derivatives and compounds able to interfere with pathogen growth such as antibiotics.

The main goal of this study was to identify the Trichoderma species to which Trichoderma IZR D-11 belongs, and to evaluate its protective effect against a powdery mildew disease of common oak (Q. robur) and to investigate the production of secondary metabolites that may affect the interaction between our isolate and the host Q. robur. The assessment of health of oak seedlings was based on visual disease symptoms and monitoring of physiological changes (chlorophyll fluorescence and water index).

Material and methods

Source of isolates

Trichoderma sp. IZR F 183 D-11 (called here Trichoderma IZR D-11) was obtained from the Institute of Plant Protection, Minsk, Belarus and is available from the Belarusian Collection of Non-pathogenic Microorganisms as Trichoderma sp. БИМ F-457Д. To observe morphological features, the isolate was grown on several media e.g. wort-agar (Merck KGaA, Germany), PDB (HiMedia Laboratories Pvt. Limited, India), Sabouraud agar (Sigma-Aldrich, Milwaukee, WI, USA), and Czapek medium (DifcoTM, MD USA). Colony color and shape and fruiting structures were assessed after 5 days of fungus incubation at 25 °C in the dark. Cultivation of Trichoderma IZR D-11 for inoculation was done in liquid media on PDB at 25 °С and 200 rpm (PSU-20i, e-biogenet, Poland). The microscopic characters of mycelia were investigated on a Nikon Eclipse H550L microscope, using the Imaging Software package NIS-Elements F 2.30, SP3 (Build 374) (Nikon).

Every year, powdery mildew develops on mature Q. robur trees in the Sękocin forest [52°6′21”N, 20°51′47″E] near Warsaw, Poland. In 2016, leaves from oak trees with powdery mildew were used as the source of natural inoculum of the powdery mildew fungus in the first year. In subsequent years, infected seedlings which were also in the greenhouse, provided the source of inoculum. Oak seedlings were infected by E. alphitoides via direct contact and dispersal of ascospores via air in the greenhouse. Infected leaves were attached to healthy ones with rubber band and in addition infected potted seedlings were mixed with healthy ones.

Molecular identification of fungal isolates

Trichoderma IZR D-11 was grown in liquid culture for 72 h, and spores and hyphae were collected. Oak leaves bearing white patches of Erysiphe sp. were collected for species identification. Spores of the powdery mildew were collected by scraping off white powdery material from few heavily infected common oak leaves collected in July 2016.

Genomic DNA for both organisms was extracted using a NucleoSpin Plant II Kit (Macherey-Nagel, Düren, Germany), following the manufacturer’s instructions. For each sample, DNA was eluted in 50 μl H2O, and extracted DNA was stored at −20 °C. For identification of Trichoderma IZR D-11, three markers were used, i.e. internal transcribed spacer region (ITS), translation elongation factor 1 alpha (Tef-1a) and RNA polymerase II subunit 2 (Rpb2). The internal transcribed spacer ITS amplification was done as described by Hermosa et al. (2000) using forward ITS1 (5′- TCCGTAGGTGAACCTGCGG-3′) and reverse ITS4 (5′- TCCTCCGCTTATTGATATGC-3′) primers (White et al. 1990). Translation elongation factor 1 alpha (Tef-1a) amplification was done as described by Carbone and Kohn (1999) using forward EF1-728F (5’-CATCGAGAAGTTCGAGAAGG-3′) and reverse TEF1LLErev (5’-AACTTGCAGGCAATGTGG-3′) primers. RNA polymerase II subunit 2 (Rpb2) amplification was done as described by Liu et al. (1999) using forward fRPB2-5f (5’-GAYGAYMGWGATCAYTTYGG-3′) and reverse fRPB2-7cr (5’-CCCATRGCTTGYTTRCCCAT-3′) primers. PCR primers were synthesized by Sigma-Aldrich (Milwaukee, WI, USA). For identification of the powdery mildew, the ITS1 and ITS4 primers were used with PCR conditions described by White et al. (1990).

PCR reactions were performed in 20 μl volumes each containing 1 μM each primer pair (forward and reverse), 10 μl 2x RedTaq Ready Mix (Sigma-Aldrich, Milwaukee, WI, USA) and 2 μl of genomic DNA as template. Cycling was performed using a Veriti 96-Well Thermal Cycler (ThermoFisher Scientific, Waltham, MA, USA) as follows: an initial denaturation step of 95 °C for 3 min, followed by 35 cycles of 94 °C for 30 s, 55 °C for 30 s, and 72 °C for 30 s, and a final extension step.

All PCR products were visualized by UV illuminator after running through 1% agarose gel and staining with ethidium bromide (Sigma-Aldrich, Milwaukee, WI, USA). Excess dNTPs and unincorporated primers were removed from the PCR product using the Clean-Up Purification Kit (A&A Biotechnology, Gdynia, Poland). Finally, DNA was eluted in 40 μl H2O and the obtained amplicons were sequenced with an ABI 3500xL genetic analyser (ThermoFisher Scientific, Waltham, MA, USA) following the manufacturer’s procedure. For species identification, a criterion of 99% identity over a minimum of 500 base pairs was applied. The assignment of sequences obtained during this study to species was performed by the BOLD Identification System (Ratnasingham and Hebert 2007) and BLAST (Altschul et al. 1990).

Plant material and treatment conditions

In spring 2015, 62 acorns of Q. robur collected from the Sękocin forest were planted singly in pots in a 2:1 mixture of horticulture soil (pH 5.5) and vermiculite, and placed outdoors close to adult oaks showing powdery mildew infections. In spring 2016, 42 healthy one-year-old potted oak seedlings were chosen for the experiment with Trichoderma and 20 were left as control seedlings and placed separately in the greenhouse. Diseased leaves from mature oaks infected by powdery mildew were collected in 2016 and rubbed onto leaves of all 42 seedlings to propagate the fungus, prior to their placement in the greenhouse for further growth and infection of leaves.

On 27 July 2016, when the first symptoms of powdery mildew developed on seedling leaves, one group of 21 was sprayed weekly for three weeks with a spore suspension in MiliQ water (1 × 107 spores/ml) of Trichoderma IZR D-11. The other set of 21 seedlings not treated with Trichoderma were randomly interspersed among the treated plants. All 42 seedlings were kept in the greenhouse for the next three years to monitor long term effects of Trichoderma IZR D-11 application on seedlings exposed to powdery mildew.

Rating scale for powdery mildew infection

Powdery mildew was assessed once a year in mid-September on the 42 seedlings using a five-point damage rating scale (Khvasko 2004) as follows:

  • 0 healthy leaves (symptomless);

  • 1 little disease (mycelium coverage of 25% of leaves surface);

  • 2 medium disease (mycelium coverage from 26% to 50% of leaves surface);

  • 3 strong disease (mycelium coverage from 51% to 75% of the leaves surface);

  • 4 severe disease (mycelium coverage more than 76% of the leaves surface).

The formula for calculation of biological efficacy (BE) of Trichoderma IZR D-11 reflecting disease suppression was as follows:

BE (%) = 100 − (100 × disease incidence of treatment/disease incidence of control) (Shin et al. 2017).

In order to evaluate if disease levels differed significantly before and after Trichoderma treatment, an analysis of variance (one-way ANOVA at P = 0.05) was conducted as well as the nonparametric Kruskal-Wallis test of significance between treatments during the three consecutive weeks using STATISTICA 12 (SoftStat, Warszawa, Polska).

Detection of Trichoderma in oak leaves and rhizosphere soil by qPCR

In mid-July 2017, a year after inoculation with Trichoderma IZR D-11, ten oak leaves from 10 seedlings with no symptoms of powdery mildew disease were randomly collected from 21 treated plants to check for the presence of Trichoderma by quantitative PCR. In consecutive years 2018 and 2019, leaves from 10 seedlings and soil from 10 pots were randomly sampled, to be tested for presence of T. asperellum, which had been identified from the sequencing results. From each leaf, about 100 mg of tissue was used for DNA extraction. In parallel, rhizosphere soil from 10 pots (0.5 g per pot) was also subjected to DNA extraction. Genomic DNA from leaves was extracted using a NucleoSpin Plant II Kit, and gDNA from soil with the NucleoSpin Soil Kit (Macherey-Nagel, Düren, Germany), following manufacturer’s instructions. Primers specific for ITS sequences of Q. robur (AH005165.2) and the T. asperellum (MT133310.1) were designed using Primer-BLAST (Ye et al. 2012).

Briefly, qPCR reactions were performed in 20 μl volumes each containing 1 μM of each primer (forward and reverse), 10 μl 2x FastStart SYBR Green Master Mix (Sigma-Aldrich, Milwaukee, WI, USA) and 2 μl of genomic DNA as template. Thermocycling conditions consisted of initial denaturation at 95 °C for 3 min, followed by 40 cycles at 95 °C for 30 s, 55 °C for 30 s and 72 °C for 30 s. Amplifications were performed using a RotorGene 6000 (Qiagen, Hilden, Germany) following manufacturer’s instructions. Ct values below 35 were considered a positive detection response. The negative controls (no DNA) produced no amplification products over 40 cycles.

Analysis of VOCs emitted by Trichoderma asperellum IZR D-11

Volatile secondary metabolites of the identified strain IZR D-11 as Trichoderma asperellum were analyzed with headspace solid phase microextraction and gas chromatography-mass spectrometry (HS-SPME/GC-MS) methods. The research was performed using previously developed methodology (Stocki et al. 2018, 2019; Isidorov et al. 2019; Oszako et al. 2019). Petri dishes with Trichoderma asperellum IZR D-11 were sealed with parafilm and warmed for 30 min at 40 °C. Afterwards, the parafilm was pierced by a needle protecting the SPME fiber (Supelco, Bellefonte, PA, USA) and the DVB/CAR/PDMS (divinylbenzene/carboxen/polydimethylsiloxane) fiber was exposed to the headspace gas phase. After 30 min of exposure, the SMPE fiber was placed for 10 min in an injection port of the GC–MS apparatus for thermal desorption of volatile substances from fiber. The GC-MS analyses were done on an Agilent 7890A gas chromatograph with an Agilent 5975C mass spectrometer (Agilent Technologies Inc., CA, USA). The chromatographic separation was performed using a capillary column HP-5MS (30 m × 0.25 mm × 0.25 μm) at a helium flow rate of 1 ml/min. The injector worked in a splitless mode at 250 °C, with initial column temperature at 35 °C, rising to 250 °C at 5 °C/min. The ion source and quadrupole temperatures were 230 °C and 150 °C, respectively. Electron ionization mass spectra were obtained at an ionization energy 70 eV. The detection was performed in a full scan mode from 29 to 600 a.m.u. After integration, the content of each component in the total ion current (% TIC) was calculated. Mass spectral data and calculated retention indices were used to identify components. Mass spectrometric identification was carried out with an automatic system of GC–MS data processing supplied by NIST and Wiley mass spectra libraries. Retention indices of compounds were determined taking into account C5-C40 n-alkanes retention times. To confirm reliable identification of compounds, experimental retention indices (RIexp.) were compared with literature values of retention indices (RIlit.).

Chlorophyll fluorescence and water content in leaves

In August 2016, chlorophyll fluorescence was measured directly on oak leaves with a Handy PEA (Handy Plant Efficiency Analyser) fluorometer from Hansatech Instruments Ltd. (King’s Lynn, Nortfolk, Great Britain). These measurements of physiological stress of investigated oaks were done on the central part of each leaf blade in an area of 4 mm2. In total, 10 measurements were made on randomly chosen leaves for each treated and untreated plant. Prior to measurements, seedlings were kept for 45–60 min in the dark at room temperature. The measurements were made following manufacturer instructions, with a 3500 μmol m−2 s1 light pulse for 1 s. Several photosynthesis parameters were assessed, namely initial fluorescence (Fo), maximal fluorescence (FM), Fv = FM - Fo, Tfm - maximum time value of fluorescence (Fm), and total Performance Index PI. The mean data were tested by LSD (0.05).

Water content was calculated for a middle leaf from each plant including Trichoderma-treated plants as non-treated controls following Mazets et al. (2012). The water content in the leaves was calculated as follows:

L = 100 × [(Mf – Md)/Md],

where : L – mean of leaf moisture (%),

$$ {M}_{\mathrm{f}}\ \mathrm{is}\ \mathrm{the}\ \mathrm{fresh}\ \mathrm{mass},\mathrm{and}\ {M}_{\mathrm{d}}\ \mathrm{is}\ \mathrm{the}\ \mathrm{dry}\ \mathrm{mass}. $$

Results

Cultural and morphological features of Trichoderma asperellum IZR D-11

When the diameter of the colony grown on wort-agar reached 90 mm, its macromorphological and microscopic characters were examined. The colonies were circular, edges were smooth, and the colony profiles were convex. The aerial mycelium was white in the initial period of growth, but became green with sporulation. The base substrate mycelium was white. The color of the reverse side of the colony was light yellow.

Colony zone lines were noticeable in the form of concentric rings. On PDB, the colonies were similar with regular round shape and clearly defined colony margins. The profile was also convex, and the texture was tomentose. The color of the colony was yellow-green, with the reverse side light yellow (Fig. S1). The center was crater-like, with concentric rings toward the margins. On the Sabouraud medium, colonies had the regular round shape with clearly defined margins. The profile was convex and the texture was woolly. The color of the colony was yellow-green, with a light yellow reverse side. The center was not noticeably different. On the Czapek medium colonies had regular round shape, the colony margins were irregular with diffuse hyphae. The profile was convex with a velvety texture. The color of the colony was intense green, with base light yellow on reverse side. The center was crater-like with concentric zone lines toward the margins. On all media, conidiophores were produced. The conidiophores were branched with bottle-shaped phialides, located with whorls 2–3 or singly, and 6.4 × 1.4 μm. Sporulation was abundant and collected in spherical heads. Conidia were ellipsoid, 3.1 × 2.6 μm, and greenish in mass. Chlamydospores were unicellular, terminal, and roughly spherical, measuring 20.5 × 21.2 μm (Fig. S2). These measurements were based on 100 observations.

Identification of fungal isolates

PCR amplifications of Trichoderma IZR D-11 strain DNA yielded the following amplicon sizes for each genomic region: 603 bp for ITS, 500 bp for Tef-1a, and 1069 bp for Rpb2. The BOLD Identification System placed the queried ITS sequence as Trichoderma asperellum with 99.17% identity. With BLAST, 93 records were found with 100% identity: 79 belonged to T. asperellum, 6 to Trichoderma sp., 3 to T. asperelloides, 2 to T. viridae, 2 to T. yunnanense and 1 to T. koningiosis. With BLAST, among 10 records of Rpb2 with 100% identity to the IZR D-11 sequence, three belonged to T. asperellum, three to T. asperelloides and four to Trichoderma sp. From eight records of Tef-1a (with 99% identity to the IZR D-11 sequence) three belonged to T. asperellum and five to Trichoderma sp. Based on these results, we concluded that IZR D-11 belongs to Trichoderma asperellum. The obtained sequences were deposited in GenBank and were assigned by the following accession numbers: MT133310.1, MT197117.1 and MT197116.1 for ITS, Tef-1a and Rpb2, respectively.

Sequencing of the powdery mildew ITS-amplicon resulted in a 622 bp fragment which matched Erysiphe alphitoides (MK357383.1) at 99% identity over the entire length. The obtained sequence has been deposited in GenBank (accession number KY346982).

Disease suppression with Trichoderma asperellum IZR D-11

Seedlings treated with T. asperellum IZR D-11 on 27 July 2017, showed smaller powdery mildew severity (0 - healthy, 4 - highly diseased) compared to the untreated controls (Fig. 1).

Fig. 1
figure 1

Condition of common oak leaves in August 2016 which had been treated with Trichoderma asperellum IZR-D11 in July 2016 (a), in comparison to the untreated control (b).

On a scale of 0 = healthy to 4 = diseased, the first rating before biocontrol application in in 2016 gave disease rating values of 0.57 for the treatment set of plants compared to 0.62 for the control plants which were not significantly different (Table 1). A week later after the first T. asperellum IZR D-11 treatment with continued powdery mildew development, the difference became more apparent and was statistically significant with disease rating at 1.00 for treated plants and 1.48 for untreated plants. After the third treatment, the difference was statistically significant with 2.10 for the treated plants and 2.81 for untreated control plants (Table 1). Based on ANOVA, the group of treated plants significantly differ from the control (p = 0.032), independently of the treatment conditions (days).

Table 1 Effect of T. asperellum IZR D-11 treatment on oak powdery mildew caused by E. alphitoides with disease ratings of 0 = healthy and 4 = highly diseased.

The T. asperellum IZR D-11 treatment significantly reduced levels of powdery mildew on the leaves of oak seedlings for the subsequent three years in the greenhouse, where they were maintained and monitored annually (Table 2).

Table 2 Residual and biological efficacy of Trichoderma asperellum IZR D-11 inoculated three times in the first year on the development of oak powdery mildew E. alphitoides in the subsequent three years where disease was rated as 0 = healthy and 4 = highly diseased.

Detection of Trichoderma asperellum in oak leaves and rhizosphere soil by qPCR

Primers specific to Q. robur (AH005165.2) ITS were designed by Primer-BLAST: the forward primer, 5’-ACCTGCACAGCGGAACG-3′, has a Tm of 60.3 °C, and the reverse primer, 5′- AATTCATTAGACGCCGACCG −3′, has a Tm of 58.4 °C. The predicted amplicon size was 191 bp. The same way were designed primers specific to T. asperellum (MT133310.1) ITS: 5’-CCTGCGGAGGGATCATTACC-3′, has a Tm of 60.0 °C, and the reverse primer, 5’-AAGAAATACGTCCGCGAGGG-3′, has a Tm of 55.0 °C. The predicted amplicon size was 169 bp. Accuracy of designed primers were tested on DNA extracted from a T. asperellum IZR D-11 culture and Q. robur seedlings.

The overall qPCRs analysis performed with the ITS primers specific for T. asperellum and Q. robur yielded positive results with DNA for T. asperellum (Ct 23.18–25.26) and Q. robur (Ct 20.04–23.12) while primers specific T. asperellum not yielded any positive signal for DNA of Q. robur. Similarly, primers designed for Q. robur not yielded any signal with DNA of T. asperellum (Table S1). None of leaves tested gave positive result of Trichoderma (Ct >40), however the fungus was present in all sampled rhizosphere soil of treated oaks (Ct 26.83–33.08). To exclude possibility that DNA extracted from leaves contains inhibitors, the qPCR analysis was performed with primers specific for Q. robur. qPCR and this yielded positive results for all analyzed samples (Ct 20.04–24.12). This indicated that extracted DNA was free of inhibitors and that the lack of signals with primers specific for T. asperellum was associated with lack of T. asperellum DNA in analyzed samples.

Volatile secondary metabolites of Trichoderma asperellum

To better understand the mechanism of biocontrol activity, the secondary metabolites produced by the fungus were assessed. A chromatogram with HS-SPME/GC-MS analysis of VOCs emitted by T. asperellum is presented in the Fig. 2. The volatiles emitted by T. asperellum contained 49 compounds (Table S2). Chemical analysis revealed VOC’s including ones unique to T. asperellum (Table S2). Sesquiterpenes were the largest group of volatiles emitted by T. asperellum, and their content (% of TIC) reached nearly 80%. Daucane type sesquiterpenes e.g. daucene, dauca-4(11),8-diene and isodaucene dominated among them.

Fig. 2
figure 2

Chromatogram with HS-SPME/GC-MS analysis of VOCs emitted by Trichoderma asperellum.

The composition of T. asperellum VOCs was also rich in aliphatic alcohols, especially 3-hepten-1-ol. Aliphatic carbonyls such as acetone, 3-octanone and acetoin were found, as well as esters, including ethyl tiglate, ethyl 3-hydroxybutanoate and methyl 2-ethylhexanoate. T. asperellum also emitted phenylethyl alcohol, linalool and 1,2-dimethylsulfide (Table S2).

Fluorometric analysis of oak leaves and water content

Fluorometric analysis in oak seedling leaves treated with T. asperellum IZR-D1 showed that the total performance index of chlorophyll fluorescence (PI) was higher (nearly double) than in the control (Table S3). Total performance index showed statistically significant difference (p < 0.05) between control (2.37) and treatment with Trichoderma (4.13); (Table S3). Similar significant differences were observed for other parameters, Fm (Fp), Tfm and Area.

After treatment with T. asperellum IZR D-11, all investigated parameters of chlorophyll fluorescence increased, but the PI increased the most (by 75%) in seedlings treated with Trichoderma compared to the control (Fig. 3). The area parameter above the fluorescence transient for the T. asperellum IZR D-11 treatment was also greater.

Fig. 3
figure 3

A spider plot of the chlorophyll fluorescence parameters of oak leaves. Abbreviations: Fo – minimum chlorophyll fluorescence; Fm (Fp) – maximum chlorophyll fluorescence; Fv/fm – maximum quantum efficiency of photosystem II; Tfm - maximum time value of fluorescence (Fm); Area - the area above the chlorophyll fluorescence curve between Fo and Fm; PI - total performance index of chlorophyll fluorescence (also called productivity index).

The water content measured in 42 leaves via wet and dry weight assessment, was significantly higher in oaks treated with T. asperellum IZR D-11 than in the control plants (Table S4).

Discussion

Biological control efficacy

The test results with efficacy over three years after application of the T. asperellum IZR D-11 as a BCA has not been previously reported in Poland or elsewhere as far as we know. In 2017 (after BCA treatment in 2016), disease rating in the control plants reached 3.22 out of 4 which was the highest disease scoring reflecting heavy infection of leaves caused by E. alphitoides, where no symptoms of powdery mildew occurred on the leaves of plants treated with T. asperellum. After three years, molecular tests (qPCR) revealed that in asymptomatic leaves, Trichoderma was not detected, but it was present in the rhizosphere soil of treated oak seedlings.

Increased total performance index (PI) reflects plant tolerance to adverse environmental conditions, meaning that plants can cope better with the stress. A lower area above the fluorescence curve between Fo and Fm in the control treatment indicated that the transfer of electrons from the reaction centers to the quinine pool is disrupted. Area reflects any changes in the shape of the induction kinetics between Fo and Fm, that is not evident on the basis of other parameters such as Fo, Fm, Fv/fm, which express only change amplitude extreme values Fo and Fm. The Tfm value increased by 6%, which suggested that the plants needed the prolonged time to reach a full reduction state of the primary electron acceptors and the significant reduction of PSII performance. Generally, T. asperellum IZR D-11 treatment improved the photosynthetic efficiency compared to the control seedlings.

The antagonistic fungus T. asperellum positively influenced not only efficiency of oak seedlings photosynthesis but also increased the total water content in oak leaves. We supposed that BCA can serve as a preventive measure to reduce the energy losses in the process of water transpiration in infected plants. Water stress is one of important factors, which can dramatically change physiological condition of plants. Prosyannikova (2002) noted that infection of oaks Q. petraea by the pathogenic fungus E. alphitoides is accompanied by pathological changes in the water regime of plants. She also showed an increase of transpiration in diseased leaves. Enhanced water consumption in the process of transpiration is not fully compensated by absorption of water by roots, increases the water deficit in leaves of infected plants, and reduce their overall water content. Changes in cell metabolism of water Q. petraea, arising under the influence of the microspheres, lead to the inhibition of growth and, ultimately, can cause the death of the undergrowth part of plants. Our findings revealed that the use of antagonistic fungus positively affected the total content of water in oak leaves, and BCA can serve as a preventive measure to reduce energy losses in the process of water transpiration in infected plants.

VOCs produced by Trichoderma

It has been reported that the VOCs released by Trichoderma spp. play an important role in its antifungal activity (Vinodkumar et al. 2017; Sunpapao et al. 2018); they can induce systemic resistance (Yi et al. 2009) or promote plant growth (Vinale et al. 2006). The Trichoderma genus employs several modes of action contributing to biocontrol activity, and these modes vary much depending on the strain and environmental conditions. Some strains act directly against pathogens via mycoparasitism, competition, or antibiosis (Vinale et al. 2008; Monfil and Casas-Flores 2014; Vos et al. 2015); others promote plant growth (Aly and Manal 2009; Christopher et al. 2010; Yadav et al. 2011) or induce natural plant protection at the site of infection as well as at a distance from the pathogen, resulting in induced plant systemic resistance (López-Bucio et al. 2015; Vos et al. 2015; Zhang et al. 2016). We speculate that the curative effect (the lightly infected leaf blades went from a low level to even less disease) in 2016 on common pedunculate oak seedlings was achieved because of mycoparasitism and plant growth promotion, but the residual efficacy over the next three years might have been a consequence of the Trichoderma asperellum being able to establish in the plant or on the plant surface, and induce systemic resistance in subsequent seasons.

Furthermore some of the volatiles produced by T. asperellum are known to have antifungal activity and this might contribute to the decreased disease. To our knowledge, some of the volatile compounds identified were recorded for the first time in T. asperellum and they are inhibitory to powdery mildew fungi (Contreras-Cornejo et al. 2016; Li et al. 2018; Guo et al. 2019). We found aliphatic carbonyls such as acetone, 3-octanone and acetoin, as well as esters, including ethyl tiglate, ethyl 3-hydroxybutanoate and methyl 2-ethylhexanoate. Acetoin is a known activator of induced systemic resistance (ISR) (Rudrappa et al. 2010). Confirmation of this speculation can be done in future research to reveal the mechanism of Trichoderma secondary metabolites and their role in the disease biocontrol process. The volatile bioactive metabolites of Trichoderma such as 6-pentyl-α-pyrone are known to exhibit multiple actions such as inhibiting growth of fungal mycelium, spore germination, and pigmentation of plant pathogenic fungi (Salwan et al. 2019). Also, Bhardwaj and Kumar (2017) reported that T. asperellum can produce many important VOCs that allow diverse functions ranging from anti-pathogenic to plant growth promotion. Our observations of Trichoderma in leaves where it had not been applied suggests that it can spread and establish. As well, the finding of Trichoderma in rhizosphere soil is consistent with induction of plant resistance perhaps mediated by VOCs. Defense responses elicited by Trichoderma in Arabidopsis are complex and involve the canonical defense hormones salicylic acid (SA) and jasmonic acid (JA), as well as phytoalexin camalexin, which may be important factors in boosting plant immunity. Co-cultivation of plant roots with these fungi also induces hydrogen peroxide accumulation which correlates with induction of pathogenesis-related reporter markers pPr1a:uidA and pLox2:uidA (Contreras-Cornejo et al. 2011). Another strain of T. asperellum T34 conferred protection against Pseudomonas syringe through regulation of proteins involved in adaptation to stress, isoprenoid and ethylene biosynthesis, photorespiration and carbohydrate metabolism (Segarra et al. 2009; Yedidia et al. 1999).

Volatile organic compounds play an important role in the protection of plants against pests and pathogens. Volatiles produced by T. asperellum were tested against Fusarium oxysporum, Fusarium solani, Macrophomina phaseolina and Pythium ultimum (Boat et al. 2020). As a result, it was found that fungus produced volatile and non-volatile inhibitory metabolites as well as hydrolytic enzymes (chitinase, cellulase, protease and lipase). Li et al. (2018) found that T. asperellum, in contrast to T. virens or T. harzianum, did not secrete greater amount of volatiles in response to F. oxysporum infection. Moreover, Wonglom et al. 2019 demonstrated that VOC’s of T. asperellum can mediate antifungal activity, defense response and promote growth of lettuce (Lactuca sativa).

In recent years, VOC’s emitted by different strains of T. asperellum, such as: T1 (Wonglom et al. 2020), LF15 (Sim et al. 2019), LU1370 (Nieto-Jacobo et al. 2017) and IsmT5 (Kottb et al. 2015), have been analyzed. These investigations revealed the emission of 4 (Kottb et al. 2015) to maximum 22 compounds (Kottb et al. 2015; Wonglom et al. 2019), whereas we detected 49 volatiles emitted by T. asperellum IZR D-11. Our results confirmed results from previous studies showing the presence of aliphatic alcohols e.g. 2-methyl-1-butanol and 2-ethyl-1-hexanol in T. asperellum emissions (Sim et al. 2019; Wonglom et al. 2019), as well as sesquiterpenes and monoterpenes (Nieto-Jacobo et al. 2017; Sim et al. 2019; Wonglom et al. 2020), and aliphatic esters and carbonyls (Kottb et al. 2015). Moreover, in volatile emissions of T. asperellum, we identified compounds previously detected in other Trichoderma species but not previous reported for T. asperellum including sesquiterpenes such as β-farnesene found in T. atroviride (Stoppacher et al. 2010) and T. asperelloides (Gal-Hemed et al. 2011), β-bisabolene found in T. atroviride (Stoppacher et al. 2010; Polizzi et al. 2011) and T. asperelloides (Gal-Hemed et al. 2011), and β-elemene and germacrene D in T. virens (Nieto-Jacobo et al. 2017). It is quite likely that the compounds mentioned play crucial roles in the plant defense. For instance, Watanabe et al. (1990) found that sesquiterpenes of the daucane type exhibited antifungal activity against various yeast and dermatophytes, and effectively inhibit growth of Candida albicans.

Aforementioned compounds recorded in our strain of T. asperellum IZR D-11 have not been previously reported for other Trichoderma species (Siddiquee 2014), therefore future work should investigate the role of these compounds in plant disease resistance. Our results suggest that antagonistic microorganism T. asperellum IZR D-11 has a potential as BCA in the control of oak powdery mildew in forest nurseries.

Conclusions

  1. 1)

    DNA sequencing of the Belarusian strain Trichoderma IZR D-11 identified it as Trichoderma asperellum which can act as a potential Biological Control Agent against Erysiphe alphitoides, the cause of oak powdery mildew.

  2. 2)

    The volatile secondary metabolites identified below were reported for the first time in T. asperellum and may play an important role in suppression of disease caused by E. alphitoides: sesquiterpenes - daucene, dauca-4(11),8-diene isodaucene; aliphatic alcohols - 3-hepten-1-ol; aliphatic carbonyls - acetone, 3-octanone and acetoin, as well as esters, including: ethyl tiglate, ethyl 3-hydroxybutanoate and methyl 2-ethylhexanoate also phenylethyl alcohol, linalool and 1,2-dimethylsulfide.

  3. 3)

    The antagonistic fungus T. asperellum not only increased oak seedling photosynthesis but also increased the total water content of oak leaves, which can serve preventively to reduce energy losses during water transpiration in infected plants.

  4. 4)

    Applications of T. asperellum in the first year was able to suppress disease in the subsequent three years and hence is promising for longer-lasting control of oak powdery mildew particularly in nurseries.