Next Article in Journal
Peronosporales Species Associated with Strawberry Crown Rot in the Czech Republic
Previous Article in Journal
Characteristics and Outcomes of Patients with Invasive Pulmonary Aspergillosis and Respiratory Tract Aspergillus Colonization from a Tertiary University Hospital in Thailand
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Fungal Endophyte Penicillium olsonii ML37 Reduces Fusarium Head Blight by Local Induced Resistance in Wheat Spikes

by
Edward C. Rojas
1,2,3,*,
Birgit Jensen
1,2,
Hans J. L. Jørgensen
2,4,
Meike A. C. Latz
2,4,5,
Pilar Esteban
1,2,6 and
David B. Collinge
1,2,*
1
Section for Microbial Ecology and Biotechnology, Department of Plant and Environmental Sciences, University of Copenhagen, Thorvaldsensvej 40, 1871 Copenhagen, Denmark
2
Copenhagen Plant Science Centre, University of Copenhagen, Thorvaldsensvej 40, 1871 Copenhagen, Denmark
3
Chr Hansen A/S, Højbakkegård Alle 30, 2630 Tåstrup, Denmark
4
Section for Plant and Soil Science, Department of Plant and Environmental Sciences, University of Copenhagen, Thorvaldsensvej 40, 1871 Copenhagen, Denmark
5
SciLifeLab, KTH Royal Institute of Technology, 171 65 Solna, Sweden
6
Department of Agricultural, Food and Agro-Environmental Sciences, University of Pisa, Via del Borghetto 80, 56124 Pisa, Italy
*
Authors to whom correspondence should be addressed.
J. Fungi 2022, 8(4), 345; https://doi.org/10.3390/jof8040345
Submission received: 7 March 2022 / Revised: 17 March 2022 / Accepted: 23 March 2022 / Published: 25 March 2022

Abstract

:
The fungal endophyte Penicillium olsonii ML37 is a biocontrol agent of Fusarium head blight in wheat (caused by Fusarium graminearum), which has shown a limited direct inhibition of fungal growth in vitro. We used RNA-seq and LC-MS/MS analyses to elucidate metabolic interactions of the three-way system Penicillium–wheat–Fusarium in greenhouse experiments. We demonstrated that P. olsonii ML37 colonises wheat spikes and transiently activates plant defence mechanisms, as pretreated spikes show a faster and stronger expression of the defence metabolism during the first 24 h after pathogen inoculation. This effect was transient and the expression of the same genes was lower in the pathogen-infected spikes than in those infected by P. olsonii alone. This response to the endophyte includes the transcriptional activation of several WRKY transcription factors. This early activation is associated with a reduction in FHB symptoms and significantly lower levels of the F. graminearum metabolites 15-acetyl-DON and culmorin. An increase in the Penicillium-associated metabolite asperphanamate confirms colonisation by the endophyte. Our results suggest that the mode of action used by P. olsonii ML37 is via a local defence activation in wheat spikes, and that this fungus has potential as a novel biological alternative in wheat disease control.

1. Introduction

Fusarium head blight (FHB) limits wheat productivity around the world [1]. Infection occurs during anthesis, and the disease affects the development of the kernels directly. In addition to reducing yield, the disease affects grain quality due to contamination with mycotoxins such as fumonisins, trichothecenes and zearalenone [2]. These mycotoxins pose a health risk when consumed [3]. The management of FHB remains a challenge due to limited genetic resistance and narrow windows for fungicide applications [4].
Fusarium graminearum and F. culmorum in particular pose problems. Fungal spores germinate and the pathogens colonise the plant’s surface using runner hyphae within 12 h after initial contact with the spikelets. Infection occurs between 24 and 48 h after pathogen spore arrival. Hyphae branch and aggregate, forming appressoria and infection cushions [5]. The fungi initially exhibit a biotrophic phase, where they spread through the intercellular spaces asymptomatically. Nevertheless, molecular interactions take place during this phase. Both mycotoxin production and plant defence mechanisms are highly active during this stage [6]. A later necrotrophic phase is most likely induced by nutrient depletion within plant tissues [7]. During this later phase, fungal hyphae degrade colonised tissues, causing bleached spike symptoms [6]. Subsequently, the pathogens produce aerial sporodochia on the necrotized tissue.
Wheat plants detect Fusarium spp. infection during the early stages of the interaction. The high expression of transcription factors associated with PAMP recognition and signal transduction proteins, such as WRKY and bZIP, takes place during the first phase of infection [8]. In general, this is associated with an activation of the salicylic acid (SA) pathway during the biotrophic phase and the jasmonate-ethylene pathway during the necrotrophic phase, although the roles of these hormones are still under debate [9,10].
Pathogenesis-related genes such as PR-1, PR-2 (β-1,3-glucanase), PR-3 (chitinase), PR-4, PR-5 (thaumatin-like protein) and oxidative stress-related enzymes are highly expressed during the first 12 h after inoculation and peak between 36 and 48 h [11,12,13]. Most of these responses have been associated with QTLs in partially resistant cultivars [14]. However, complete genetic resistance in wheat cultivars remains elusive and partial resistance is polygenic and scattered across the genome [15]. Interestingly, transgenic approaches for conferring resistance against Fusarium are promising, if not accepted globally [16,17,18,19].
Biological control is receiving increased attention as an alternative approach for disease control because it can provide new modes of action whilst enjoying high acceptance among consumers [20,21].
Several studies have shown different levels of FHB biocontrol by living organisms [22]. Bacteria from the genera Bacillus and Pseudomonas have been shown to reduce levels of FHB in greenhouse studies when vegetative cells or spores were sprayed at anthesis [23,24,25,26]. Among fungal candidates, Clonostachys rosea, Cryptococcus sp. and Sarocladium zeae have also been shown to suppress FHB in greenhouse studies when spores were applied at heading [27,28,29]. Most of these studies using microbes have shown an ability to control the disease, but little is known about their mode of action.
Recently, we demonstrated how the fungal endophyte Penicillium olsonii ML37, isolated from wheat leaves, reduced the severity of Septoria tritici blotch [30] and Fusarium head blight [31]. Thus, Rojas et al. [31] found that ML37 reduced FHB symptoms and pathogen biomass in wheat spikes when inoculated before anthesis in greenhouse assays. Interestingly, ML37 did not inhibit F. graminearum growth in vitro and required a lag period before pathogen inoculation to provide a significant biocontrol effect. This suggests that the disease-reducing effect of ML37 is mediated by the plant after inoculation with the pathogen.
The aim of this study is to investigate the mode of action of the biocontrol mediated by P. olsonii ML37 against FHB and, specifically, whether the activation of host defence responses is involved. We perform a transcriptional analysis using the RNA-seq analysis on endophyte-treated and nontreated spikes. Likewise, we evaluate the transcriptomic profile of treated spikes after F. graminearum infection and analyse the content of mycotoxins and other specialised metabolites in the treated spikes. Our findings suggest that biocontrol mediated by ML37 is at least in part triggered by the activation of inherent plant resistance.

2. Materials and Methods

2.1. Plant and Fungal Material

Spring wheat cv. Diskett (SW 45456-Lantmännen) was grown to heading stage in a greenhouse. This cultivar is moderately susceptible to Fusarium head blight. The endophyte Penicillium olsonii ML37 was isolated in our previous study [30]. Agar plugs were cultured on potato dextrose agar (PDA) (Difco™, BD, Franklin Lakes, NJ, USA) plates for 7 days. Spores were harvested by flooding the plate using autoclaved MilliQ water + Tween 20 (0.05% v/v). The spore suspension was filtered using Miracloth® (Merck, Darmstadt, Germany), centrifuged and resuspended. The P. olsonii spore concentration was adjusted to 107 spores per ml using MilliQ water (Merck KGaA, Darmstadt, Germany) + Tween 20 (Thermo Scientific™, Waltham, Massachusetts, USA) (0.05% v/v).
F. graminearum strain WC-091-7 was kindly provided by Dr. Lise N. Jørgensen, Aarhus University. Agar plugs were added to 20 mL RA liquid medium (Imholte and Schramm, 1968) and cultured for 10 days. Macroconidia were collected by filtering the culture through autoclaved Miracloth® (Merck). Spores were centrifuged and resuspended and the F. graminearum spore concentration was adjusted to 105 macroconidia per ml using autoclaved MilliQ water + Tween 20 (0.05% v/v).
In order to identify the transcriptomic and enzymatic changes on wheat spikes induced by P. olsonii ML37 alone (before anthesis) and in the presence of Fusarium infection (after anthesis), four separate experiments were performed, namely, experiments 1, 1.1, 2 and 2.2.
Experiment 1. P. olsoniiML37 transcriptomic effect on heading-stage wheat spikes before anthesis. Heading-stage wheat spikes (BBCH 58) were sprayed with either a P. olsonii ML37 spore suspension or autoclaved MilliQ water + Tween as a mock treatment using a 100 mL spray bottle. Approximately 2 mL was sprayed onto each spike until run-off. Treatments were applied at late afternoon to prevent exposure of the treated spikes to direct sunlight. The spikes were left uncovered and harvested at 48 h after inoculation. Samples were immediately snap-frozen in liquid nitrogen and stored at −80 °C until RNA extraction. A total of three replications (spikes) per treatment was harvested.
Experiment 1.1. β-1,3-glucanase activity of P. olsonii ML37-treated heading-stage wheat spikes before anthesis. Experiment 1 was repeated to evaluate the β-1,3-glucanase activity at three different time points (24, 48 and 72 h), following the protocol of Shetty et al. [32]. Four replications of each condition were used. Total protein was extracted from 200 mg ground tissue using 250 µL 0.05 M sodium acetate buffer (pH 5.2). Specific β-1,3-glucanase activity was calculated based on the absorbance at 540 nm using the absorbance in an Eon™ Microplate Spectrophotometer reader (Winooski, VT, USA). Specific enzyme activity was calculated as μmol released product/min for every mg of initial protein. Data were analysed using a generalized linear model. Comparisons between treatments at each time point were considered significant at p < 0.05.
Experiment 2. P. olsonii ML37 transcriptomic effect on wheat spikes during Fusarium graminearum infection after anthesis. Heading-stage wheat spikes were inoculated with P. olsonii ML37 spores or MilliQ water + Tween as described above. Between three and five days later, as spikes reached mid-anthesis stage (BBCH 65), all the spikes were inoculated with a F. graminearum WC-091-7 macroconidial suspension (105 spores/mL) by spraying approximately 2 mL until run-off. Treated spikes were covered with a previously mist-treated transparent plastic bag and sealed using adhesive tape around the straw to increase humidity around the spike. Plastic bags were maintained until harvest of each time point. Spikes were harvested at 24 and 72 h after inoculation with the pathogen. In total, four replications were harvested for each condition and time point. Spikes were snap frozen in liquid nitrogen immediately after harvest and stored at −80 °C until RNA extraction.
The remaining spikes were inspected for FHB severity at 7 days after inoculation. Severity was calculated as the number of spikelets exhibiting visible FHB symptoms out of the total number of spikelets in each spike. Statistical analysis of FHB severity was performed using a generalized linear model specifying a binomial distribution. All statistical comparisons were performed against FHB controls and were considered significant when p ≤ 0.05.
Experiment 2.1. β-1,3-glucanase activity of P. olsonii ML37-treated wheat spikes during Fusarium graminearum infection. Experiment 2 was repeated to evaluate the β-1,3-glucanase activity at three different time points (24, 48 and 72 h). Four replications of each condition and time point were used. Samples were analysed as described above.
The remaining spikes were harvested 7 days after pathogen inoculation and stored at −20 °C for multi-mycotoxin and specialised (secondary) metabolite analysis.

2.2. RNA Extraction and Illumina Sequencing

Harvested spikes from experiments 1 and 2 were used for RNA extraction. Six spikelets with rachis were dissected from the middle of each treated spike and ground up in liquid nitrogen using a pestle and mortar. Approximately 100 mg of ground tissue was used for RNA extraction. RNA was extracted using the RNeasy Plant Mini Kit (Qiagen, Hilden, Germany). RNA samples were treated with a DNAse treatment using the DNA-free™ DNA removal Kit (Thermo Fisher Scientific, Waltham, MA, USA). RNA quantity and quality were assessed using a NanoDrop 2000 Spectrophotometer (Thermo Scientific™, Waltham, MA, USA) and visualised on a 2% agarose gel. RNA integrity was verified using an Agilent 2100 Bioanalyzer (Thermo Scientific™). In total, 22 samples were used for library preparation and sequencing (6 for experiment 1 and 16 for experiment 2).
Library preparation and Illumina high-throughput sequencing was performed by the NovoGene Sequencing company (Hong Kong, China) using the HiSeq paired-end platform (Illumina, San Diego, CA, USA). Sequencing depth was at least 50 million paired-end reads per sample. Sequence quality and error rate was verified using Qphred scores.

2.3. Bioinformatic Analysis

Raw data were filtered as follows: after demultiplexing, reads containing adapters, reads with Q-score ≤ 5 and reads containing more than 10% N calls were removed. Pseudoalignment of the reads was performed using the kallisto program [33]. Default parameters were used and 50 bootstraps were set up to estimate the technical error in the downstream analysis. Reads were pseudoaligned to the wheat transcriptome (The International Wheat Genome Sequencing Consortium–IWGSC RefSeq v2.0 assembly) [34].
Nonaligned reads from the second experiment were then pseudoaligned against the Fusarium graminearum PH1 transcriptome. Unmapped reads from both experiments were then pseudoaligned to the transcriptome of three closely related Penicillium spp. for reference (P. chrysogenum str. P2niaD18, P. expansum str. CMP-1 and P. griseofulvum str. PG3.). All the genomes are available at http://ensemblgenomes.org (accessed on 6 March 2022).

2.4. Gene Ontology Analysis

In order to identify clusters of transcripts with similar expression patterns, we used the self-organizing map (SOM) competitive learning algorithm from the kohonen package in R [35,36]. Transcripts with fewer than 0.5 tpm were discarded. The mean of each time point and condition was calculated. These values were normalised to 0 with standard deviation from −1 to 1. SOMs were laid down using 100 nodes. Clusters were subdivided using hierarchical clustering. The final number of clusters was decided based on the point were within-group variation reached a low plateau [37]. Individual clusters of genes were then analysed for the enrichment of gene ontology terms (GO terms) with a 0.05 cut-off for false discovery rate using the goseq Bioconductor package [38].

2.5. Differential Gene Expression Analysis of Specific Genes Involved in Fusarium–Wheat Interaction

Clusters of genes, such as WRKY transcription factors, pathogenesis-related proteins (PR), antifungal proteins and specialised metabolites, reported to be active during Fusarium–wheat interactions were examined. Similarly, genes naturally expressed during anthesis and grain filling stage in wheat were inspected. Furthermore, Fusarium genes involved in trichothecene production such as the TRI cluster and CLM1 were also analysed.
Differentially expressed genes were calculated at transcript level using a Wald test on ‘beta’ coefficients (analogous to Fold-change value) using the sleuth tool [39]. This value uses the bootstraps values from kallisto to estimate the ‘inferential variance’ or technical error. p-values were aggregated to determine gene-level results and then adjusted for multiple comparisons using a 0.1 cut-off for false discovery rate and transformed into q-values.

2.6. Multi-Mycotoxin and Specialised Metabolite Analysis

Wheat spikes from experiment 2.2 were used to quantify mycotoxins and specialised metabolites. The plant material was freeze dried for 48 h and homogenized using ceramic beads (10 mm) in a mechanical shaker (SO-40a, Fluid Management Europe B.V.) using a 30 s programme five times. For each sample, 3–5 full spikes were pooled to achieve 1 g of dried material and four replications per treatment were prepared. Five grams of a previously ground sample was extracted for 90 min using 20 mL of acetonitrile/water/acetic acid 79/20/1 (v/v/v). Particles in the raw extracts were removed by sedimentation by keeping the sample suspensions in an upright position for 10 min. An aliquot of the clear extract was diluted 1 + 1 with acetonitrile/water/acetic acid 20/79/1 and, subsequently, analysed (injection volume: 5 μL).
The samples were analysed with a multi-mycotoxin method based on LC-MS/MS analyses using an Agilent 1290 HPLC coupled to an Applied Biosystems 5500 QTrap mass spectrometer. Over 500 different compounds, including major mycotoxins, were measured simultaneously: trichothecenes, zearalenone, aflatoxin, fuminosins, T-2 toxin, depsipeptides, bacterial metabolites, Fusarium metabolites and other unspecific metabolites. This was performed at the Department for Agrobiotechnology in Tulln, part of the University of Natural Resources and Life Sciences (BOKU) in Austria, a research partner of BIOMIN® according to the methods described previously [40]. Mycotoxin concentration data were analysed using a linear model. Comparisons between treatments were considered significant at p < 0.05.

3. Results

3.1. Penicillium olsonii ML37 Activates Defence Mechanisms at Heading-Stage Wheat Spikes

P. olsonii ML37 inoculation changed the transcriptomic profile of wheat spikes at 48 h after inoculation. The SOM clustering and gene ontology analysis detected a group of transcripts highly expressed in ML37-treated plants. This cluster was significantly enriched with genes involved in defence responses, specifically defence to fungi, as well as stress-related and respiratory burst (Figure 1A).
The differential expression analysis showed that P. olsonii ML37 inoculation upregulated 105 genes and downregulated 48 genes (Supplementary Materials Table S1). Among the upregulated genes, several transcription factors involved in biotic and abiotic stress responses such as WIR1A, WRKY51, WRKY10 [41,42,43] were detected. Some key pathogenesis-related genes, such as PR-1 and PR-2, showed increase expression, although they were not statistically significantly upregulated (Figure 2). While some others such as PR1-16 and chitinase IV (Cht4) were significantly upregulated (Figure S1).
In addition, a key regulator of abscisic acid (ABA) catabolism, the ABA 8′-hydroxylase ABA8 gene [44], was significantly overexpressed.
Specific β-1,3-glucanase activity showed a significant increase at 48 h on P. olsonii ML37-treated spikes compared to the mock-treated ones, but the enzymatic activity declined and was not significant at 72 h (Figure S2A). This confirmed that ML37 inoculation activates microbe-related defence responses in the plant after inoculation.
Among the downregulated genes, there was a significant enrichment of GO terms related to cell proliferation, cell cycle, DNA regulation and the heterochromatin structure (Figure 1B). Several H2 and H3 histone genes were downregulated in ML37-treated plants. Notably, the largest part of transcripts in the samples (approximately 50%) remained unaffected by the treatments (Figure 1C). Among these, GO terms associated with photosynthesis and biological processes, such as translation, were over-represented. Several genes reported to be highly active during anthesis, such as α-amylase and ATP-dependent 6-phosphofruktinase [45,46] involved in carbohydrate transport, were not differentially expressed (Figure S1).

3.2. P. olsonii ML37-Treated Plants Show Higher Levels of Defence Responses during the First 24 h of FHB Attack

Wheat response to Fusarium infection varied drastically over time between ML37-treated and mock-treated plants. At 1 dpi, 293 genes were differentially expressed between treatments, while at 3 dpi, this number increased to 18,442 genes (Tables S2 and S3). A large cluster of genes significantly enriched with GO terms involved in response to chitin, cell death, pathogenesis, defence response pathways and others showed a higher expression levels in ML37-inoculated plants at 1 dpi and a steep decrease at 3 dpi (Figure 3A).
Within this cluster, key marker genes for the pathogenesis and activation of defence responses, such as PR-1-1, PR-1-9, PR-2, PR-4A [47,48] and others with antifungal activity such as chitinases, endochitinases and specialised metabolite production, were significantly overexpressed in ML37-treated plants compared to mock-treated plants at 1 dpi (Figure S3). These same genes showed the opposite trend at 3 dpi, as mock-treated plants showed a higher expression than the ML37-treated ones (Figure 4).
ML37-treated wheat spikes under FHB infection showed an increased specific β-1,3-glucanase activity as early as 24 dpi after pathogen inoculation compared to the mock-treated ones. This effect was greater at 48 h and declined at 72 dpi (Figure S2B).
An equally large second cluster of transcripts showed lower levels of expression at 1 dpi in ML37-treated plants and a higher expression at 3 dpi, while remaining constantly expressed in mock-treated plants across both time points (Figure 3B). This cluster was significantly enriched with GO terms associated with thylakoid, chloroplast relocation and photosynthesis, as well as floral development, cell division, starch metabolic process and postembryonic development.
A third, smaller cluster with constantly high levels of expression at both 1 and 3 dpi in the ML37-treated plants while also being expressed at low levels in the mock-treated plants was detected. This was significantly enriched with GO terms involved in cell wall formation and maintenance such as pectin esterase activity (Figure 3C).

3.3. Fusarium graminearum Shows a Reduced Metabolism during P. olsonii ML37-Mediated Biocontrol

Only 2–15% of the unmapped reads were successfully aligned to the Fusarium graminearum transcriptome. Such low read quantities impeded a clearer visualization of the transcriptomic state of the pathogen during infection using Gene Ontology tools. Still, SOM clustering identified several groups of genes that showed reduced expression patterns on ML37-treated plants compared to the mock-treated ones at both time points (Figure S4). None of these clusters were significantly enriched with any GO terms after p-value adjustment, but the most abundant were related to ribosomal function and mRNA translation.
On the other hand, 112 genes were detected as differentially expressed in F. graminearum on ML37-treated spikes compared to mock-treated spikes at 1 dpi. This number increased to 228 at 3 dpi (Tables S4 and S5). Interestingly, some of the genes significantly overexpressed in the ML37-treated spikes were involved in stress-related responses. Among these were ABC transporter multidrug resistance protein 1, cytochrome P45O monooxygenase and antifungal defensin (Figure S4C).
We further inspected the expression of several fungal genes involved in trichothecene biosynthesis, namely, TRI-6, TRI-5 and TRI-12 [49]. These genes showed low expression levels at 3 dpi, but they were not significantly affected by ML37 inoculation at this early stage. Interestingly, CLM1, a gene involved in the synthesis of culmorin-like metabolites, was significantly repressed in ML37-treated plants at 3 dpi (Figure 5).
Spikes were assessed for FHB severity at 7 dpi and a significant reduction in infection was observed on the ML37-treated plants (Figure 6). The LC-MS/MS analysis of over 500 microbial specialised metabolites in full spikes at 7 dpi revealed lower concentrations of trichothecenes in ML37-treated spikes compared to mock-treated spikes, including statistically significant reductions in 15-acetyl-deoxynivalenol and culmorin (Table 1).

3.4. Expression of Penicillium Genes Not Detected, but Fungal Metabolites Were Detected on Treated Spikes

Unmapped reads from experiment two were pseudoaligned to the genomes of three closely related Penicillium spp., i.e., P. chrysogenum, P. expansum and P. griseofulvum [50]. On average, only 0.5 of the reads used aligned against these genomes. Moreover, there were no differences in the number of reads mapped between the ML37-treated and nontreated plants. Thus, we concluded that the detection of ML37 transcripts was unsuccessful.
The multi-mycotoxin method detected significantly higher concentrations of the fungal metabolites N-benzoyl-phenylalanine and asperphenamate on the ML37-treated spikes compared to mock-treated spikes after Fusarium infection (Table 1). These molecules were not detected on the untreated controls.

4. Discussion

Penicillium olsonii ML37 has been identified as a new potential biocontrol agent against wheat pathogens [30,31]. Here, we present the transcriptomic response of wheat spikes after inoculation with P. olsonii ML37 alone and during Fusarium infection. Our results highlighted the importance of timing behind the biological control of FHB as ML37-treated plants showed a reduced disease severity and mycotoxin accumulation.
The transcriptomic data showed that P. olsonii ML37 spores that landed on the wheat spikes were recognized by the plant. Hence, the transcription factors WRKY10, WIR1A and WRKY51, which are associated with the response to stress [51], were upregulated by ML37. These genes are also highly responsive to the exogenous application of purified chitin in Arabidopsis [52]. Moreover, the expression of WRKY10 and WIR1A has been associated with the deployment of basal resistance in rice [53] as well as the downstream activation of pathogen responsive genes in the cell wall [42]. As far as we are aware, this was the first time WRKY gene induction was seen following colonisation with an endophytic fungus.
The recognition of P. olsonii ML37 triggered plant defence mechanisms that led to the local expression of PR-proteins and a significant increase in β-1,3-glucanase activity in treated spikes. This is a common response after interactions with beneficial organisms. It is part of a basal defence mechanism called induced resistance [54]. Similar responses have been observed in other crops after inoculation with different microorganisms, such as Trichoderma in maize [55,56] and Bacillus spp. in soybean [57], and they have also been reported during interactions between plants and endophytic fungi [58,59].
The basal defence response to P. olsonii ML37 caused an active reprogramming of metabolic processes in the wheat plants. Histone-related genes and chromatin configuration genes were significantly downregulated after ML37 inoculation. Epigenetic regulation has been suggested to play a key role in the deployment of defence responses in general and, specifically, in the onset of defence responses elicited by beneficial microbes [60,61]. Simultaneously, this activation reduced photosynthetic activity on wheat spikes. This is a common consequence of plant defence activation, and its effects on plant fitness have been discussed previously [62]. However, it has been shown that changes in photosynthesis and growth-related genes do not change the expression of normal anthesis developmental processes [45,46].
Fusarium-inoculated spikes that were pretreated with P. olsonii ML37 spores showed a higher expression of genes involved in pathogen defence, the oxidative burst, PR proteins, chitinases, glycosyltransferases as well as a higher specific β-1,3-glucanase activity at 1 dpi than in the mock-treated spikes. Interestingly, several PR-genes (PR1-1, PR1-9, PR1-16, PR4A and PR4B) showed a greater expression in the mock-treated spikes compared to the ML-37-treated ones at 3 dpi, suggesting a delayed response for these key defence genes. By having a higher expression of these genes at 24 h, the pretreated plants were able to minimize the initial infection, while the mock-treated plants showed a relatively lower expression at this point. This pattern was reversed at 3 dpi, where the mock-treated plants showed a higher expression of these genes. This suggested a stronger defence response to a high infection level of Fusarium. This showed that activating defence responses during the first 24 h after pathogen arrival is the key driver of the P. olsonii ML37-mediated effect. Thus, our data strongly suggest that the timing of the defence response deployment is a key factor in the mode of action of ML37-mediated biocontrol, as the early activation of defence mechanisms was associated with FHB suppression.
The rapid activation and deployment of defence responses has been associated with FHB resistance in several cultivars. A greater activation of the defence genes PR-1 and PR-2 at 6-12 hpi was observed in the tolerant cultivar Sumai 3 [11]. In wheat, these genes have been reported to be highly responsive to fungal pathogens [11,63,64].
We showed that specific β-1,3-glucanase activity peaked at 48 hpi, but decreased at 72 hpi in the ML37-treated compared to the mock-treated plants. These proteins are produced after both the SA and JA pathway activation, representing the activity of several different enzymes. They inhibit pathogen growth by hydrolysing the main chain of β-1,3-glucans in fungal cell walls [65,66]. They have been observed in transcriptomic and proteomic studies in FHB-tolerant cultivars [12,67].
We observed a high expression of glycosyltransferase and sulfotransferase proteins in ML37-treated compared to mock-treated spikes. These enzymes use glucosinolates, brassinosteroids, jasmonates and flavonoids as substrates to generate specialised metabolites. These types of enzymes, especially uridinediphosphate glycosyltransferases and glutathione S-transferases, have been associated with major FHB-disease resistance QTLs [68], presumably by the detoxification of DON into less toxic forms such as D3G [15]. In our study, the simultaneous activation of antifungal enzymes and specialised metabolite production was associated with a reduced translation and ribosomal activity in the pathogen, and likely to be the cause of the reduced FHB symptoms in ML37-treated plants.
The expression of three key trichothecene pathway genes during the first 72 h was not affected by the P. olsonii ML37 treatment. However, there was a significant reduction in the accumulation of trichothecenes and other Fusarium metabolites in whole wheat spikes after 7 days. The sesquiterpenoid mycotoxins, 15-Acetyl-DON and culmorin, were significantly reduced in ML37-treated plants. 15-Acetyl-DON is the acetylated derivate of DON and possesses equally regulatory importance, whereas culmorin is a co-occurring Fusarium metabolite, which has been suggested to be an “emerging mycotoxin” [69]. It enhances DON accumulation by inhibiting DON detoxification via glycosylation in plants [70]. Interestingly, Wipfler et al. [71] reported that culmorin biosynthesis genes were more highly to be expressed than the TRI5 trichothecene biosynthetic gene during the first week after infection. This suggests that ML37-mediated defence activation in wheat spikes prevents culmorin production during the first 72 h after infection, consequently, reducing trichothecene production and the onset of the necrotrophic phase of FHB. This was further validated by the RNA-seq data showing the repression of the gene CLM1 involved in the synthesis of culmorin.
Previous studies have suggested an indirect mode of action of some biocontrol agents. For example, Jochum et al. [72] showed that Lysobacter enzymogenes C3 heat-inactivated cultures reduced FHB to the same extent as fresh cultures when sprayed 8 days before F. graminearum inoculation. Other studies have used plant hormone accumulation as an indication of an increased defence activation. Palazzini et al. [23] suggested the phytohormone induction of SA and JA by Bacillus velezensis RC-218 when co-inoculated with Fusarium graminearum, but this effect was not consistently higher than the F. graminearum control during the first 72 h, while the bacteria alone showed no differences compared to the mock (water) control. Similarly, Kemp et al. [29] showed that spikes pre-inoculated with the endophytic yeast Sarocladium zeae showed higher levels of both SA and JA during FHB compared to the F. graminearum-only control. However, this effect was only significant at 1 h after pathogen inoculation for SA and up to 8 h for JA. In neither of these studies was there direct evidence of defence mechanisms known to inhibit Fusarium spp, and/or their dynamics before and after pathogen inoculation. Our study was the first to fully describe the effective biological control of FHB through plant defence activation.
P. olsonii ML37 has shown significant biocontrol effects against two major pathogens in wheat, namely, Z. tritici and F. graminearum [30,31]. In this study, we were unable to detect P. olsonii gene expression. However, we were able to detect high amounts of the linear amino acid ester asperphenamate and one of its components, N-benzoylphenylalanine, only in the ML37-treated spikes, further confirming an active role of ML37. These molecules are produced by fungi such as Aspergillus and Penicillium and they are particularly common among endophytic fungi [73]. They are biologically active and have been linked to several chemical and pharmaceutical applications [74], as well as displaying antimicrobial properties [75]. However, P. olsonii ML37 does not inhibit the growth of Z. tritici nor F. graminearum in vitro. It is possible that asperphenamate and N-benzoylphenylalanine are only produced in planta. In fact, asperphenamate has been detected in diverse groups of unrelated plants, which could suggest a common generalist fungal endophyte origin [74,76].
Multiple and simultaneous modes of action are a common trait of biocontrol agents [77], especially in fungal endophytes [78]. For example, the biocontrol agent Clonostachys rosea is a mycoparasite [79] with the ability to detoxify the Fusarium mycotoxin zearalenone [80]. Similarly, mycoparasitism and competition were shown in biological control by Trichoderma gamsii T6085 against Fusarium spp. [81]. Although it is likely that ML37 possesses several modes of action, this remains to be verified using targeted experiments.
Recently, it was shown how fhb7, a gene that confers the resistance gene against Fusarium head blight in wheat, has a fungal endophyte origin and was transmitted to wild relatives of wheat via a horizontal gene transfer from Epichloë fungi [82]. This highlights the degree of coevolution between fungal endophytes and their host plants. Indeed, naturally occurring fungal endophytes have been observed to be associated with healthy wheat spikes and negatively correlated with Fusarium abundance [83]. Here, we showed how P. olsonii ML37 increased wheat tolerance to biotic stress caused by F. graminearum. This capacity has the potential to be used as a novel biocontrol strategy against FHB in wheat production.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jof8040345/s1, Figure S1: transcriptome changes in wheat spikes at 48 h after P. olsonii ML37 inoculation, Figure S2: specific β-1,3-glucanase activity in wheat spikes inoculated with combinations of P. olsonii ML37, Fusarium or water at 24, 48 and 72 hpi., Figure S3: expression of five different chitinase genes and five genes involved in specialised metabolite production in wheat in P. olsonii ML37-treated plants at 24 h after Fusarium inoculation, Figure S4: transcriptome changes in Fusarium graminearum during P. olsonii ML37-mediated biocontrol of FHB at 24 and 72 hpi, Table S1: differentially expressed genes in wheat at 48 h after P. olsonii ML37 inoculation, Table S2: differentially expressed genes in wheat—P. olsonii ML37-mediated biocontrol of FHB (24 hpi), Table S3: differentially expressed genes in wheat—P. olsonii ML37-mediated biocontrol of FHB (72 hpi), Table S4: differentially expressed genes in F. graminearum during P. olsonii ML37-mediated biocontrol of FHB (24 hpi), Table S5: differentially expressed genes in F. graminearum during P. olsonii ML37-mediated biocontrol of FHB (72 hpi).

Author Contributions

Conceptualization, E.C.R., B.J., H.J.L.J. and D.B.C.; methodology, E.C.R., B.J., H.J.L.J., M.A.C.L. and D.B.C.; formal analysis, E.C.R. and H.J.L.J.; investigation, E.C.R. and P.E.; resources, D.B.C., E.C.R. and M.A.C.L.; data curation, E.C.R.; writing—original draft preparation, E.C.R.; writing—review and editing, E.C.R., B.J., H.J.L.J., M.A.C.L. and D.B.C.; visualization, E.C.R.; supervision, E.C.R., B.J., H.J.L.J. and D.B.C.; project administration, D.B.C.; funding acquisition, D.B.C., B.J. and H.J.L.J. All authors have read and agreed to the published version of the manuscript.

Funding

This project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 674964. Pilar Esteban was supported by the Erasmus+ programme of the EU.
Jof 08 00345 i001

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The sequencing data that support the findings of this study are available in the NCBI repository under SRA submission SUB11165978. Other experimental data are available upon request.

Acknowledgments

We would like to thank technicians Ulla Rasmussen and Monika Dragsted for their dedicated work in the lab and in the greenhouse, respectively. Special thanks to Cristobal Uauy from the John Innes Centre, UK, and Thomas Baldwin from USDA-ARS, USA, for their advice during data analysis, as well as to Lise Nistrup Jørgensen from Aarhus University for the Fusarium graminearum isolates. Additionally, thanks to Michael Sulyok, University of Natural Resources and Life Sciences, Department for Agrobiotechnology, Tulln, Austria, via BIOMIN Holding GmbH for metabolite analysis.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Savary, S.; Willocquet, L.; Pethybridge, S.J.; Esker, P.; McRoberts, N.; Nelson, A. The global burden of pathogens and pests on major food crops. Nat. Ecol. Evol. 2019, 3, 430–439. [Google Scholar] [CrossRef] [PubMed]
  2. Bryła, M.; Waśkiewicz, A.; Ksieniewicz-Woźniak, E.; Szymczyk, K.; Jędrzejczak, R. Modified Fusarium Mycotoxins in Cereals and Their Products—Metabolism, Occurrence, and Toxicity: An Updated Review. Molecules 2018, 23, 963. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Antonissen, G.; Martel, A.; Pasmans, F.; Ducatelle, R.; Verbrugghe, E.; Vandenbroucke, V.; Li, S.; Haesebrouck, F.; Van Immerseel, F.; Croubels, S. The impact of Fusarium mycotoxins on human and animal host susceptibility to infectious diseases. Toxins 2014, 6, 430–452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Rojas, E.C.; Jørgensen, H.J.L.; Jensen, B.; Collinge, D.B. Fusarium diseases: Biology and management perspectives. In Integrated Disease Management of Wheat and Barley; Oliver, R.P., Ed.; Burleigh Dodds Science Publishing: Cambridge, UK, 2019; pp. 23–45. [Google Scholar]
  5. Boenisch, M.J.; Schäfer, W. Fusarium graminearum forms mycotoxin producing infection structures on wheat. BMC Plant Biol. 2011, 11, 110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Brown, N.A.; Evans, J.; Mead, A.; Hammond-Kosack, K.E. A spatial temporal analysis of the Fusarium graminearum transcriptome during symptomless and symptomatic wheat infection. Mol. Plant Pathol. 2017, 18, 1295–1312. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Kazan, K.; Gardiner, D.M. Transcriptomics of cereal–Fusarium graminearum interactions: What we have learned so far. Mol. Plant Pathol. 2018, 19, 764–778. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Erayman, M.; Turktas, M.; Akdogan, G.; Gurkok, T.; Inal, B.; Ishakoglu, E.; Ilhan, E.; Unver, T. Transcriptome analysis of wheat inoculated with Fusarium graminearum. Front. Plant Sci. 2015, 6, 867. [Google Scholar] [CrossRef] [Green Version]
  9. Foroud, N.A.; Pordel, R.; Goyal, R.K.; Ryabova, D.; Eranthodi, A.; Chatterton, S.; Kovalchuk, I. Chemical activation of the ethylene signaling pathway promotes Fusarium graminearum resistance in detached wheat heads. Phytopathology 2019, 109, 796–803. [Google Scholar] [CrossRef]
  10. Sarowar, S.; Alam, S.T.; Makandar, R.; Lee, H.; Trick, H.N.; Dong, Y.; Shah, J. Targeting the pattern-triggered immunity pathway to enhance resistance to Fusarium graminearum. Mol. Plant Pathol. 2019, 20, 626–640. [Google Scholar] [CrossRef] [Green Version]
  11. Pritsch, C.; Muehlbauer, G.J.; Bushnell, W.R.; Somers, D.A.; Vance, C.P. Fungal development and induction of defense response genes during early infection of wheat spikes by Fusarium graminearum. Mol. Plant-Microbe Interact. 2000, 13, 159–169. [Google Scholar] [CrossRef] [Green Version]
  12. Ding, L.; Xu, H.; Yi, H.; Yang, L.; Kong, Z.; Zhang, L.; Xue, S.; Jia, H.; Ma, Z. Resistance to hemi-biotrophic F. graminearum Infection is associated with coordinated and ordered expression of diverse defense signaling pathways. PLoS ONE 2011, 6, e19008. [Google Scholar] [CrossRef] [Green Version]
  13. Khaledi, N.; Taheri, P.; Falahati-Rastegar, M. Reactive oxygen species and antioxidant system responses in wheat cultivars during interaction with Fusarium species. Australas. Plant Pathol. 2016, 45, 653–670. [Google Scholar] [CrossRef]
  14. Buerstmayr, M.; Steiner, B.; Buerstmayr, H. Breeding for Fusarium head blight resistance in wheat—Progress and challenges. Plant Breed. 2020, 139, 429–454. [Google Scholar] [CrossRef]
  15. He, Y.; Wu, L.; Liu, X.; Jiang, P.; Yu, L.; Qiu, J.; Wang, G.; Zhang, X.; Ma, H. TaUGT6, a novel UDP-Glycosyltransferase gene enhances the resistance to FHB and DON accumulation in wheat. Front. Plant Sci. 2020, 11, 1549. [Google Scholar] [CrossRef]
  16. Makandar, R.; Essig, J.S.; Schapaugh, M.A.; Trick, H.N.; Shah, J. Genetically engineered resistance to Fusarium head blight in wheat by expression of Arabidopsis NPR1. Mol. Plant-Microbe Interact. 2006, 19, 123–129. [Google Scholar] [CrossRef] [Green Version]
  17. Li, H.P.; Zhang, J.B.; Shi, R.P.; Huang, T.; Fischer, R.; Liao, Y.C. Engineering Fusarium Head Blightresistance in wheat by expression of a fusion protein containing a fusarium-specific antibody and an antifungal peptide. Mol. Plant-Microbe Interact. 2008, 21, 1242–1248. [Google Scholar] [CrossRef] [Green Version]
  18. Collinge, D.B.; Jørgensen, H.J.L.; Lund, O.S.; Lyngkjær, M.F. Engineering pathogen resistance in crop plants—Current trends and future prospects. Annu. Rev. Phytopathol. 2010, 48, 269–291. [Google Scholar] [CrossRef] [Green Version]
  19. Collinge, D.B.; Sarrocco, S. Transgenic approaches for plant disease control: Status and prospects 2021. Plant Pathol. 2022, 71, 207–225. [Google Scholar] [CrossRef]
  20. Collinge, D.B.; Jensen, D.F.; Rabiey, M.; Sarrocco, S.; Shaw, M.W.; Shaw, R. Biological control of plant diseases—What has been achieved and what is the direction? Plant Pathol. 2022, 10, 596. [Google Scholar] [CrossRef]
  21. Collinge, D.B.; Jørgensen, H.J.L.; Latz, M.A.C.; Manzotti, A.; Ntana, F.; Rojas, E.C.; Jensen, B. Searching for novel fungal biological control agents for plant disease control among endophytes. In Endophytes: For a Growing World; Hodkinson, T.R., Doohan, F.M., Saunders, M., Murphy, B.R., Eds.; Cambridge University Press: Cambridge, UK, 2019; pp. 25–51. [Google Scholar]
  22. Legrand, F.; Picot, A.; Cobo-Díaz, J.F.; Chen, W.; Le Floch, G. Challenges facing the biological control strategies for the management of Fusarium Head Blight of cereals caused by F. graminearum. Biol. Control 2017, 113, 26–38. [Google Scholar] [CrossRef]
  23. Palazzini, J.; Roncallo, P.; Cantoro, R.; Chiotta, M.; Yerkovich, N.; Palacios, S.; Echenique, V.; Torres, A.; Ramirez, M.; Karlovsky, P.; et al. Biocontrol of Fusarium graminearum sensu stricto, reduction of deoxynivalenol accumulation and phytohormone induction by two selected antagonists. Toxins 2018, 10, 88. [Google Scholar] [CrossRef] [Green Version]
  24. Baffoni, L.; Gaggia, F.; Dalanaj, N.; Prodi, A.; Nipoti, P.; Pisi, A.; Biavati, B.; Di Gioia, D. Microbial inoculants for the biocontrol of Fusarium spp. in durum wheat. BMC Microbiol. 2015, 15, 242. [Google Scholar] [CrossRef]
  25. Schisler, D.A.; Khan, N.I.; Boehm, M.J. Biological control of fusarium head blight of wheat and deoxynivalenol levels in grain via use of microbial antagonists. In Mycotoxins and Food Safety; Advances in Experimental Medicine and Biology; DeVries, J.W., Trucksess, M.W., Jackson, L.S., Eds.; Springer: Boston, MA, USA, 2002; Volume 504. [Google Scholar]
  26. Da Luz, W.C.; Stockwell, C.A.; Bergstrom, G.C. Biological control of Fusarium graminearum. In Fusarium Head Blight of Wheat and Barley; Canty, S., Lewis, J., Ward, R., Eds.; APS Press: St. Paul, MN, USA, 2003. [Google Scholar]
  27. Khan, N.I.; Schisler, D.A.; Boehm, M.J.; Slininger, P.J.; Bothast, R.J. Selection and evaluation of microorganisms for biocontrol of fusarium head blight of wheat incited by Gibberella zeae. Plant Dis. 2001, 85, 1253–1258. [Google Scholar] [CrossRef] [Green Version]
  28. Xue, A.G.; Chen, Y.; Voldeng, H.D.; Fedak, G.; Savard, M.E.; Längle, T.; Zhang, J.; Harman, G.E. Concentration and cultivar effects on efficacy of CLO-1 biofungicide in controlling Fusarium head blight of wheat. Biol. Control 2014, 73, 2–7. [Google Scholar] [CrossRef]
  29. Kemp, N.D.; Vaughan, M.M.; McCormick, S.P.; Brown, J.A.; Bakker, M.G. Sarocladium zeae is a systemic endophyte of wheat and an effective biocontrol agent against Fusarium head blight. Biol. Control 2020, 149, 104329. [Google Scholar] [CrossRef]
  30. Latz, M.A.C.; Jensen, B.; Collinge, D.B.; Lyngs Jørgensen, H.J. Identification of two endophytic fungi that control Septoria tritici blotch in the field, using a structured screening approach. Biol. Control 2020, 141, 104128. [Google Scholar] [CrossRef]
  31. Rojas, E.C.; Jensen, B.; Jørgensen, H.J.L.; Latz, M.A.C.; Esteban, P.; Ding, Y.; Collinge, D.B. Selection of fungal endophytes with biocontrol potential against Fusarium head blight in wheat. Biol. Control 2020, 144, 104222. [Google Scholar] [CrossRef]
  32. Shetty, N.P.; Jensen, J.D.; Knudsen, A.; Finnie, C.; Geshi, N.; Blennow, A.; Collinge, D.B.; Jørgensen, H.J.L. Effects of b-1,3-glucan from Septoria tritici on structural defence responses in wheat. J. Exp. Bot. 2009, 60, 4287–4300. [Google Scholar] [CrossRef] [Green Version]
  33. Bray, N.L.; Pimentel, H.; Melsted, P.; Pachter, L. Near-optimal probabilistic RNA-seq quantification. Nat. Biotechnol. 2016, 34, 525–527. [Google Scholar] [CrossRef]
  34. Zimin, A.V.; Puiu, D.; Hall, R.; Kingan, S.; Clavijo, B.J.; Salzberg, S.L. The first near-complete assembly of the hexaploid bread wheat genome, Triticum aestivum. GigaScience 2017, 6, gix097. [Google Scholar] [CrossRef] [Green Version]
  35. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  36. Wehrens, R.; Kruisselbrink, J. Flexible Self-Organizing Maps in kohonen 3.0. J. Stat. Softw. 2018, 87, 1–18. [Google Scholar] [CrossRef] [Green Version]
  37. Vesanto, J.; Alhoniemi, E. Clustering of the self-organizing map. IEEE Trans. Neural Netw. 2000, 11, 586–600. [Google Scholar] [CrossRef]
  38. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. Gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biol. 2010, 11, R14. [Google Scholar] [CrossRef] [Green Version]
  39. Pimentel, H.; Bray, N.L.; Puente, S.; Melsted, P.; Pachter, L. Differential analysis of RNA-seq incorporating quantification uncertainty. Nat. Methods 2017, 14, 687–690. [Google Scholar] [CrossRef]
  40. Sulyok, M.; Stadler, D.; Steiner, D.; Krska, R. Validation of an LC-MS/MS-based dilute-and-shoot approach for the quantification of > 500 mycotoxins and other secondary metabolites in food crops: Challenges and solutions. Anal. Bioanal. Chem. 2020, 412, 2607–2620. [Google Scholar] [CrossRef] [Green Version]
  41. Coram, T.E.; Huang, X.; Zhan, G.; Settles, M.L.; Chen, X. Meta-analysis of transcripts associated with race-specific resistance to stripe rust in wheat demonstrates common induction of blue copper-binding protein, heat-stress transcription factor, pathogen-induced WIR1A protein, and ent-kaurene synthase transcripts. Funct. Integr. Genom. 2010, 10, 383–392. [Google Scholar] [CrossRef]
  42. Wang, C.; Deng, P.; Chen, L.; Wang, X.; Ma, H.; Hu, W.; Yao, N.; Feng, Y.; Chai, R.; Yang, G.; et al. A wheat WRKY transcription factor TaWRKY10 confers tolerance to multiple abiotic stresses in transgenic tobacco. PLoS ONE 2013, 8, e65120. [Google Scholar] [CrossRef]
  43. Hu, Z.; Wang, R.; Zheng, M.; Liu, X.; Meng, F.; Wu, H.; Yao, Y.; Xin, M.; Peng, H.; Ni, Z.; et al. TaWRKY51 promotes lateral root formation through negative regulation of ethylene biosynthesis in wheat (Triticum aestivum L.). Plant J. 2018, 96, 372–388. [Google Scholar] [CrossRef] [Green Version]
  44. Chono, M.; Matsunaka, H.; Seki, M.; Fujita, M.; Kiribuchi-Otobe, C.; Oda, S.; Kojima, H.; Kobayashi, D.; Kawakami, N. Isolation of a wheat Triticum aestivum L. mutant in ABA 8′-hydroxylase gene: Effect of reduced ABA catabolism on germination inhibition under field condition. Breed. Sci. 2013, 63, 104–115. [Google Scholar] [CrossRef] [Green Version]
  45. Ma, J.; Li, R.; Wang, H.; Li, D.; Wang, X.; Zhang, Y.; Zhen, W.; Duan, H.; Yan, G.; Li, Y. Transcriptomics analyses reveal wheat responses to drought stress during reproductive stages under field conditions. Front. Plant Sci. 2017, 8, 592. [Google Scholar] [CrossRef] [Green Version]
  46. Rangan, P.; Furtado, A.; Henry, R.J. The transcriptome of the developing grain: A resource for understanding seed development and the molecular control of the functional and nutritional properties of wheat. BMC Genom. 2017, 18, 766. [Google Scholar] [CrossRef] [Green Version]
  47. Bertini, L.; Leonardi, L.; Caporale, C.; Tucci, M.; Cascone, N.; Di Berardino, I.; Buonocore, V.; Caruso, C. Pathogen-responsive wheat PR4 genes are induced by activators of systemic acquired resistance and wounding. Plant Sci. 2003, 164, 1067–1078. [Google Scholar] [CrossRef]
  48. Lu, S.; Edwards, M.C. Molecular characterization and functional analysis of PR-1-like proteins identified from the wheat head blight fungus Fusarium graminearum. Phytopathology 2018, 108, 510–520. [Google Scholar] [CrossRef] [Green Version]
  49. Kimura, M.; Tokai, T.; Takahashi-Ando, N.; Ohsato, S.; Fujimura, M. Molecularand genetic studies of Fusarium trichothecene biosynthesis: Pathways, genes, and evolution. Biosci. Biotechnol. Biochem. 2007, 71, 2105–2123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Houbraken, J.; Samson, R.A. Phylogeny of Penicillium and the segregation of Trichocomaceae into three families. Stud. Mycol. 2011, 70, 1–51. [Google Scholar] [CrossRef]
  51. Phukan, U.J.; Jeena, G.S.; Shukla, R.K. WRKY Transcription factors: Molecular regulation and stress responses in plants. Front. Plant Sci. 2016, 7, 760. [Google Scholar] [CrossRef] [Green Version]
  52. Wan, J.R.; Zhang, S.Q.; Stacey, G. Activation of a mitogen-activated protein kinase pathway in Arabidopsis by chitin. Mol. Plant Pathol. 2004, 5, 125–135. [Google Scholar] [CrossRef]
  53. Choi, N.; Im, J.H.; Lee, E.; Lee, J.; Choi, C.; Park, S.R.; Hwang, D.-J. WRKY10 transcriptional regulatory cascades in rice are involved in basal defense and Xa1-mediated resistance. J. Exp. Bot. 2020, 71, 3735–3748. [Google Scholar] [CrossRef]
  54. De Kesel, J.; Conrath, U.; Flors, V.; Luna, E.; Mageroy, M.H.; Mauch-Mani, B.; Pastor, V.; Pozo, M.J.; Pieterse, C.M.J.; Ton, J.; et al. The induced resistance lexicon: Do’s and don’ts. Trends Plant Sci. 2021, 26, 685–691. [Google Scholar] [CrossRef]
  55. Saravanakumar, K.; Li, Y.; Yu, C.; Wang, Q.-q.; Wang, M.; Sun, J.; Gao, J.-x.; Chen, J. Effect of Trichoderma harzianum on maize rhizosphere microbiome and biocontrol of Fusarium Stalk rot. Sci. Rep. 2017, 7, 1771. [Google Scholar] [CrossRef]
  56. Carella, P. Xylem-mobile oxylipins are critical regulators of induced systemic resistance in maize. Plant Cell 2019, 32, 13–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Xing, Z.; Wu, X.; Zhao, J.; Zhao, X.; Zhu, X.; Wang, Y.; Fan, H.; Chen, L.; Liu, X.; Duan, Y. Isolation and identification of induced systemic resistance determinants from Bacillus simplex Sneb545 against Heterodera glycines. Sci. Rep. 2020, 10, 11586. [Google Scholar] [CrossRef]
  58. Chen, X.; Marszałkowska, M.; Reinhold-Hurek, B. Jasmonic acid, not salicyclic acid restricts endophytic root colonization of rice. Front. Plant Sci. 2020, 10, 1758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Ntana, F.; Johnson, S.R.; Hamberger, B.; Jensen, B.; Jørgensen, H.J.L.; Collinge, D.B. Regulation of tomato specialised metabolism after establishment of symbiosis with the endophytic fungus Serendipita indica. Microorganisms 2022, 10, 194. [Google Scholar] [CrossRef] [PubMed]
  60. Ramirez-Prado, J.S.; Piquerez, S.J.M.; Bendahmane, A.; Hirt, H.; Raynaud, C.; Benhamed, M. Modify the histone to win the battle: Chromatin dynamics in plant–pathogen interactions. Front. Plant Sci. 2018, 9, 355. [Google Scholar] [CrossRef]
  61. Alonso, C.; Ramos-Cruz, D.; Becker, C. The role of plant epigenetics in biotic interactions. New Phytol. 2019, 221, 731–737. [Google Scholar] [CrossRef] [PubMed]
  62. Huot, B.; Yao, J.; Montgomery, B.L.; He, S.Y. Growth–defense tradeoffs in plants: A balancing act to optimize fitness. Mol. Plant 2014, 7, 1267–1287. [Google Scholar] [CrossRef] [Green Version]
  63. Yang, F.; Li, W.; Jørgensen, H.J.L. Transcriptional reprogramming of wheat and the hemibiotrophic pathogen Septoria tritici during two phases of the compatible interaction. PLoS ONE 2013, 8, e81606. [Google Scholar] [CrossRef] [Green Version]
  64. Zhang, J.; Wang, F.; Liang, F.; Zhang, Y.; Ma, L.; Wang, H.; Liu, D. Functional analysis of a pathogenesis-related thaumatin-like protein gene TaLr35PR5 from wheat induced by leaf rust fungus. BMC Plant Biol. 2018, 18, 76. [Google Scholar] [CrossRef] [Green Version]
  65. Akiyama, T.; Jin, S.; Yoshida, M.; Hoshino, T.; Opassiri, R.; Ketudat Cairns, J.R. Expression of an endo-(1,3;1,4)-β-glucanase in response to wounding, methyl jasmonate, abscisic acid and ethephon in rice seedlings. J. Plant Physiol. 2009, 166, 1814–1825. [Google Scholar] [CrossRef]
  66. Yoshimi, A.; Miyazawa, K.; Abe, K. Function and biosynthesis of cell wall α-1,3-glucan in fungi. J. Fungi 2017, 3, 63. [Google Scholar] [CrossRef]
  67. Zhou, W.C.; Kolb, F.L.; Riechers, D.E. Identification of proteins induced or upregulated by Fusarium head blight infection in the spikes of hexaploid wheat (Triticum aestivum). Genome 2005, 48, 770–780. [Google Scholar] [CrossRef]
  68. Venske, E.; dos Santos, R.S.; da Rosa Farias, D.; Rother, V.; da Maia, L.C.; Pegoraro, C.; Costa de Oliveira, A. Meta-Analysis of the QTLome of Fusarium Head Blight resistance in bread wheat: Refining the current puzzle. Front. Plant Sci. 2019, 10, 727. [Google Scholar] [CrossRef] [Green Version]
  69. Spanic, V.; Katanic, Z.; Sulyok, M.; Krska, R.; Puskas, K.; Vida, G.; Drezner, G.; Šarkanj, B. Multiple fungal metabolites including mycotoxins in naturally infected and Fusarium-inoculated wheat samples. Microorganisms 2020, 8, 578. [Google Scholar] [CrossRef] [Green Version]
  70. Woelflingseder, L.; Warth, B.; Vierheilig, I.; Schwartz-Zimmermann, H.; Hametner, C.; Nagl, V.; Novak, B.; Šarkanj, B.; Berthiller, F.; Adam, G.; et al. The Fusarium metabolite culmorin suppresses the in vitro glucuronidation of deoxynivalenol. Arch. Toxicol. 2019, 93, 1729–1743. [Google Scholar] [CrossRef] [Green Version]
  71. Wipfler, R.; McCormick, S.P.; Proctor, R.; Teresi, J.; Hao, G.; Ward, T.; Alexander, N.; Vaughan, M.M. Synergistic phytotoxic effects of culmorin and trichothecene mycotoxins. Toxins 2019, 11, 555. [Google Scholar] [CrossRef] [Green Version]
  72. Jochum, C.C.; Osborne, L.E.; Yuen, G.Y. Fusarium head blight biological control with Lysobacter enzymogenes strain C3. Biol. Control 2006, 39, 336–344. [Google Scholar] [CrossRef]
  73. Subko, K.; Wang, X.; Nielsen, F.H.; Isbrandt, T.; Gotfredsen, C.H.; Ramos, M.C.; Mackenzie, T.; Vicente, F.; Genilloud, O.; Frisvad, J.C.; et al. Mass spectrometry guided discovery and design of novel asperphenamate analogs from Penicillium astrolabium reveals an extraordinary NRPS flexibility. Front. Microbiol. 2021, 11, 3510. [Google Scholar] [CrossRef]
  74. Frisvad, J.C.; Houbraken, J.; Popma, S.; Samson, R.A. Two new Penicillium species Penicillium buchwaldii and Penicillium spathulatum, producing the anticancer compound asperphenamate. FEMS Microbiol. Lett. 2013, 339, 77–92. [Google Scholar] [CrossRef] [Green Version]
  75. Pomini, A.M.; Ferreira, D.T.; Braz-Filho, R.; Saridakis, H.O.; Schmitz, W.; Ishikawa, N.K.; Faccione, M. A new method for asperphenamate synthesis and its antimicrobial activity evaluation. Nat. Prod. Res. 2006, 20, 537–541. [Google Scholar] [CrossRef]
  76. Macabeo, A.P.G.; Tudla, F.A.; Alejandro, G.J.D.; Kouam, S.F.; Hussain, H.; Krohn, K. Benzoylated derivatives from Uvaria rufa. Biochem. Syst. Ecol. 2010, 38, 857–860. [Google Scholar] [CrossRef]
  77. Köhl, J.; Kolnaar, R.; Ravensberg, W.J. Mode of action of microbial biological control agents against plant diseases: Relevance beyond efficacy. Front. Plant Sci. 2019, 10, 845. [Google Scholar] [CrossRef] [Green Version]
  78. Latz, M.A.C.; Jensen, B.; Collinge, D.B.; Jørgensen, H.J.L. Endophytic fungi as biocontrol agents: Elucidating mechanisms in disease suppression. Plant Ecol. Divers. 2018, 11, 555. [Google Scholar] [CrossRef] [Green Version]
  79. Gimeno, A.; Kägi, A.; Drakopoulos, D.; Bänziger, I.; Lehmann, E.; Forrer, H.-R.; Keller, B.; Vogelgsang, S. From laboratory to the field: Biological control of Fusarium graminearum on infected maize crop residues. J. Appl. Microbiol. 2020, 129, 680–694. [Google Scholar] [CrossRef] [Green Version]
  80. Kosawang, C.; Karlsson, M.; Vélëz, H.; Rasmussen, P.H.; Collinge, D.B.; Jensen, B.; Jensen, D.F. Zearalenone detoxification by zearalenone hydrolase is important for the antagonistic ability of Clonostachys rosea against mycotoxigenic Fusarium graminearum. Fungal Biol. 2014, 118, 364–373. [Google Scholar] [CrossRef]
  81. Sarrocco, S.; Valenti, F.; Manfredini, S.; Esteban, P.; Bernardi, R.; Puntoni, G.; Baroncelli, R.; Haidukowski, M.; Moretti, A.; Vannacci, G. Is exploitation competition involved in a multitrophic strategy for the biocontrol of fusarium head blight? Phytopathology 2019, 109, 560–570. [Google Scholar] [CrossRef] [Green Version]
  82. Wang, H.; Sun, S.; Ge, W.; Zhao, L.; Hou, B.; Wang, K.; Lyu, Z.; Chen, L.; Xu, S.; Guo, J.; et al. Horizontal gene transfer of Fhb7 from fungus underlies Fusarium head blight resistance in wheat. Science 2020, 368, eaba5435. [Google Scholar] [CrossRef]
  83. Rojas, E.C.; Sapkota, R.; Jensen, B.; Jørgensen, H.J.L.; Henriksson, T.; Jørgensen, L.N.; Nicolaisen, M.; Collinge, D.B. Fusarium head blight modifies fungal endophytic communities during infection of wheat spikes. Microb. Ecol. 2020, 79, 397–408. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Transcriptome changes induced by P. olsonii ML37 in wheat spikes at 48 h after inoculation. (A) (left) Normalised transcripts per million (tpm) of a cluster of genes with positive regulation after ML37 inoculation. Grey lines represent mean expression of transcripts, blue line represents total cluster mean. (A) (right) Significantly enriched gene ontology terms within each cluster. Bar plots represent −log(p-value) of each GO term. (B) (left) Normalised tpm of a cluster of genes with negative regulation after ML37 inoculation. (B) (right) Significantly enriched GO terms. (C) (left) Normalised tpm of a cluster of genes unaffected by ML37 inoculation. (C) (right) Significantly enriched GO terms. Gene ontology conventions—BP: biological process; CC: cell component; MF: molecular function.
Figure 1. Transcriptome changes induced by P. olsonii ML37 in wheat spikes at 48 h after inoculation. (A) (left) Normalised transcripts per million (tpm) of a cluster of genes with positive regulation after ML37 inoculation. Grey lines represent mean expression of transcripts, blue line represents total cluster mean. (A) (right) Significantly enriched gene ontology terms within each cluster. Bar plots represent −log(p-value) of each GO term. (B) (left) Normalised tpm of a cluster of genes with negative regulation after ML37 inoculation. (B) (right) Significantly enriched GO terms. (C) (left) Normalised tpm of a cluster of genes unaffected by ML37 inoculation. (C) (right) Significantly enriched GO terms. Gene ontology conventions—BP: biological process; CC: cell component; MF: molecular function.
Jof 08 00345 g001
Figure 2. Expression of three WRKY transcription factors and two pathogenesis-related genes in wheat spikes after inoculation with P. olsonii ML37 alone before anthesis. Boxplot shows range and median of mean tpm for each gene. Statistically significant differences (q < 0.01 using Walt test) are denoted as ‘**’. Deep green colour denotes differentially expressed genes. Light green colour denotes non-differentially expressed genes.
Figure 2. Expression of three WRKY transcription factors and two pathogenesis-related genes in wheat spikes after inoculation with P. olsonii ML37 alone before anthesis. Boxplot shows range and median of mean tpm for each gene. Statistically significant differences (q < 0.01 using Walt test) are denoted as ‘**’. Deep green colour denotes differentially expressed genes. Light green colour denotes non-differentially expressed genes.
Jof 08 00345 g002
Figure 3. Transcriptome changes during P. olsonii ML37-mediated biocontrol of FHB at 1 and 3 dpi. (A) (left) Normalised transcripts per million (tpm) of a cluster of genes with positive regulation in ML37-treated plants after F. graminearum inoculation. Grey lines represent mean expression of transcripts, blue line represents total cluster mean. (A) (right) Significantly enriched gene ontology terms within each cluster. Bar plots represent –log(P-value) of each GO term. (B) (left) Normalised tpm of a cluster of genes with negative regulation in ML37-treated plants after pathogen inoculation. (B) (right) Significantly enriched GO terms. (C) (left) Normalised tpm of a cluster of genes highly expressed at both time points in ML37-treated plants. (C) (right) Significantly enriched GO terms. Gene ontology conventions—BP: biological process; CC: cell component; MF: molecular function.
Figure 3. Transcriptome changes during P. olsonii ML37-mediated biocontrol of FHB at 1 and 3 dpi. (A) (left) Normalised transcripts per million (tpm) of a cluster of genes with positive regulation in ML37-treated plants after F. graminearum inoculation. Grey lines represent mean expression of transcripts, blue line represents total cluster mean. (A) (right) Significantly enriched gene ontology terms within each cluster. Bar plots represent –log(P-value) of each GO term. (B) (left) Normalised tpm of a cluster of genes with negative regulation in ML37-treated plants after pathogen inoculation. (B) (right) Significantly enriched GO terms. (C) (left) Normalised tpm of a cluster of genes highly expressed at both time points in ML37-treated plants. (C) (right) Significantly enriched GO terms. Gene ontology conventions—BP: biological process; CC: cell component; MF: molecular function.
Jof 08 00345 g003
Figure 4. Expression of five different Pathogenesis-related (PR) genes in P. olsonii ML37 and mock-treated spikes at 1 and 3 dpi with F. graminearum. Boxplot shows range and median of mean tpm for each gene. Statistically significant differences (q < 0.01 using Walt test) are denoted as ‘**’.
Figure 4. Expression of five different Pathogenesis-related (PR) genes in P. olsonii ML37 and mock-treated spikes at 1 and 3 dpi with F. graminearum. Boxplot shows range and median of mean tpm for each gene. Statistically significant differences (q < 0.01 using Walt test) are denoted as ‘**’.
Jof 08 00345 g004
Figure 5. Expression of three genes from the TRI cluster and CLM1 involved in deoxynivalenol and culmorin production, respectively, on endophyte-treated spikes and FHB control at 3 dpi. Boxplot shows range and median of mean tpm for each gene. Light red colour denotes nondifferentially expressed genes. Boxplot shows range and median of mean tpm for each gene. Statistically significant differences (q < 0.01 using Walt test) are denoted as ‘**’.
Figure 5. Expression of three genes from the TRI cluster and CLM1 involved in deoxynivalenol and culmorin production, respectively, on endophyte-treated spikes and FHB control at 3 dpi. Boxplot shows range and median of mean tpm for each gene. Light red colour denotes nondifferentially expressed genes. Boxplot shows range and median of mean tpm for each gene. Statistically significant differences (q < 0.01 using Walt test) are denoted as ‘**’.
Jof 08 00345 g005
Figure 6. Effect of Penicillium olsonii ML37 on FHB severity at 5 dpi with F. graminearum under greenhouse conditions. Boxplot shows range and median severity for each treatment. Points below and up boxes are outlier data points. Statistically significant differences at p < 0.001 compared to FHB control are denoted as ‘***’.
Figure 6. Effect of Penicillium olsonii ML37 on FHB severity at 5 dpi with F. graminearum under greenhouse conditions. Boxplot shows range and median severity for each treatment. Points below and up boxes are outlier data points. Statistically significant differences at p < 0.001 compared to FHB control are denoted as ‘***’.
Jof 08 00345 g006
Table 1. HPLC-ESI-MS/MS-based multi-mycotoxin analysis of wheat spikes treated with P. olsonii ML37 at 7 dpi with Fusarium graminearum.
Table 1. HPLC-ESI-MS/MS-based multi-mycotoxin analysis of wheat spikes treated with P. olsonii ML37 at 7 dpi with Fusarium graminearum.
Concentration ppb (µg/kg)
TypeMoleculeMockP. olsonii ML37
TrichothecenesDeoxynivalenol106,602 a84,390 a
15-Acetyildeoxynivalenol7618 a4784 b
DON-3-glucoside10,511 a7929 a
Nivalenol115 a89 a
Fusarium-specific metabolitesCulmorin183,035 a133,627 b
15-Hydroxyculmorin64,928 a54,355 a
15-Hydroxyculmoron8188 a6133 a
5-hydroxyculmorin11,725 a9484 a
Unspecific metabolitesN-Benzoyl-Phenylalanine0 b34 a
Asperphenamate35 b1055 a
Altersetin8 a0 b
a and b denote statistical differences (p-value < 0.05) using Fisher’s exact test.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rojas, E.C.; Jensen, B.; Jørgensen, H.J.L.; Latz, M.A.C.; Esteban, P.; Collinge, D.B. The Fungal Endophyte Penicillium olsonii ML37 Reduces Fusarium Head Blight by Local Induced Resistance in Wheat Spikes. J. Fungi 2022, 8, 345. https://doi.org/10.3390/jof8040345

AMA Style

Rojas EC, Jensen B, Jørgensen HJL, Latz MAC, Esteban P, Collinge DB. The Fungal Endophyte Penicillium olsonii ML37 Reduces Fusarium Head Blight by Local Induced Resistance in Wheat Spikes. Journal of Fungi. 2022; 8(4):345. https://doi.org/10.3390/jof8040345

Chicago/Turabian Style

Rojas, Edward C., Birgit Jensen, Hans J. L. Jørgensen, Meike A. C. Latz, Pilar Esteban, and David B. Collinge. 2022. "The Fungal Endophyte Penicillium olsonii ML37 Reduces Fusarium Head Blight by Local Induced Resistance in Wheat Spikes" Journal of Fungi 8, no. 4: 345. https://doi.org/10.3390/jof8040345

APA Style

Rojas, E. C., Jensen, B., Jørgensen, H. J. L., Latz, M. A. C., Esteban, P., & Collinge, D. B. (2022). The Fungal Endophyte Penicillium olsonii ML37 Reduces Fusarium Head Blight by Local Induced Resistance in Wheat Spikes. Journal of Fungi, 8(4), 345. https://doi.org/10.3390/jof8040345

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop