Next Article in Journal
Myosin and Other Energy-Transducing ATPases: Structural Dynamics Studied by Electron Paramagnetic Resonance
Next Article in Special Issue
Analysis of Aldo–Keto Reductase Gene Family and Their Responses to Salt, Drought, and Abscisic Acid Stresses in Medicago truncatula
Previous Article in Journal
DNA Methylation Is Correlated with Gene Expression during Diapause Termination of Early Embryonic Development in the Silkworm (Bombyx mori)
Previous Article in Special Issue
Characterization of the Soybean GmIREG Family Genes and the Function of GmIREG3 in Conferring Tolerance to Aluminum Stress
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 2
10.3390/ijms21020673
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Glycerol-Induced Powdery Mildew Resistance in Wheat by Regulating Plant Fatty Acid Metabolism, Plant Hormones Cross-Talk, and Pathogenesis-Related Genes

by
Yinghui Li
1,2,
Lina Qiu
1,
Xinye Liu
1,3,
Qiang Zhang
1,
Xiangxi Zhuansun
1,
Tzion Fahima
2,
Tamar Krugman
2,
Qixin Sun
1 and
Chaojie Xie
1,*
1
Key Laboratory of Crop Heterosis and Utilization (MOE) and State Key Laboratory for Agrobiotechnology, Beijing Key Laboratory of Crop Genetic Improvement, China Agricultural University, Beijing 100193, China
2
Institute of Evolution, University of Haifa, 199 Abba-Hushi Avenue, Mt. Carmel, Haifa 3498838, Israel
3
Ministry of Education Key Laboratory of Molecular and Cellular Biology, College of Life Sciences, Hebei Normal University, Shijiazhuang 050024, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(2), 673; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21020673
Submission received: 8 December 2019 / Revised: 14 January 2020 / Accepted: 17 January 2020 / Published: 20 January 2020
(This article belongs to the Special Issue Abiotic Stress and Gene Networks in Plants 2020)
Browse Figures
Review Reports Versions Notes

Abstract

:
Our previous study indicated that glycerol application induced resistance to powdery mildew (Bgt) in wheat by regulating two important signal molecules, glycerol-3-phosphate (G3P) and oleic acid (OA18:1). Transcriptome analysis of wheat leaves treated by glycerol and inoculated with Bgt was performed to identify the activated immune response pathways. We identified a set of differentially expressed transcripts (e.g., TaGLI1, TaACT1, and TaSSI2) involved in glycerol and fatty acid metabolism that were upregulated in response to Bgt infection and might contribute to G3P and OA18:1 accumulation. Gene Ontology (GO) enrichment analysis revealed GO terms induced by glycerol, such as response to jasmonic acid (JA), defense response to bacterium, lipid oxidation, and growth. In addition, glycerol application induced genes (e.g., LOX, AOS, and OPRs) involved in the metabolism pathway of linolenic and alpha-linolenic acid, which are precursor molecules of JA biosynthesis. Glycerol induced JA and salicylic acid (SA) levels, while glycerol reduced the auxin (IAA) level in wheat. Glycerol treatment also induced pathogenesis related (PR) genes, including PR-1, PR-3, PR-10, callose synthase, PRMS, RPM1, peroxidase, HSP70, HSP90, etc. These results indicate that glycerol treatment regulates fatty acid metabolism and hormones cross-talk and induces the expression of PR genes that together contribute to Bgt resistance in wheat.

1. Introduction

Bread wheat (Triticum aestivum L.) is one of the most important food crops worldwide. Blumeria graminis f.sp. tritici (Bgt), is an obligate biotrophic ascomycete fungus that invades the aerial parts of wheat, causing powdery mildew disease. The damage caused by Bgt can result in yield losses from 30% to 40% in years with severe epidemics [1,2]. Deployment of resistance (R) genes and spraying with foliar fungicides are the common methods for diseases control in wheat. However, many R genes lost their effectiveness due to the appearance of new pathogen races, and fungicides are not always economically feasible and environmentally friendly. Therefore, new sustainable effective ways for disease control in wheat are needed.
Plants have been developing multifaceted innate immunity systems along the co-evolution of plants and pathogens. When plants detect an attempted pathogen invasion, the transmembrane pattern recognition receptors (PRRs) recognize pathogen-associated molecular patterns (PAMPs), resulting in PAMP-triggered immunity (PTI) [3]. The other immune response acts largely within the cell, using NOD-like receptors (NLRs) that specifically recognize fungal effectors, resulting in effector-triggered immunity (ETI) [3,4]. After recognition, several downstream signaling events are elicited, including reactive oxygen species (ROS) accumulation, and transient activation of mitogen-activated protein kinases (MAPK) signaling cascades, as well as interaction with some plant hormones including salicylic acid (SA), jasmonic acid (JA), auxin (IAA), abscisic acid (ABA) and other phytohormones [5,6,7].
Glycerol-3-phosphate (G3P) and oleic acid (OA18:1) are two important signal molecules associated with resistance to fungal pathogens in plants [8,9,10]. Arabidopsis plants overexpressing GLI1/NHO1 (encoding a glycerol kinase) increase G3P levels and enhance resistance against bacterial disease caused by Pseudomonas syringae [11]. Arabidopsis mutation gli1 reduces the G3P level and enhances susceptibility to Colletotrichum higginsianum [12,13]. It was suggested that G3P may contribute to plant resistance. However, ssi2 (encoding a stearoylacyl carrier protein fatty acid desaturase) mutants with lower OA18:1 levels elevate SA and JA levels, induce the expression of pathogenesis-related (PR) proteins, and improve plant resistance [8,14,15,16]. Several reports show that lipids are a major source of organic carbon delivered to the fungus [17,18]. Mutualistic mycorrhizal fungi recruit the plant fatty acid biosynthesis program to facilitate host invasion; with the case of kas1, kar1, and fatb-1 mutant plants of fatty acid biosynthesis enhance disease resistance compared with wild type plants [19]. Therefore, these lipids, as plant signals or fungal organic carbon sources, play an important role in the host–pathogen interactions of plants and pathogenic fungi.
Exogenous glycerol increases G3P level, with a reduction in OA18:1 by acylation of G3P with OA18:1, resulting in the induction of PR gene expression and conferring resistance to different pathogens in Arabidopsis, rice, soybean, and cacao tree [15,20,21,22]. Our previous studies indicated that glycerol can induce resistance of wheat to Bgt with increased G3P and reduced OA18:1 levels, and accumulation of SA, JA, and ROS levels [23]. Moreover, glycerol appears to have some potential to be applied in wheat fields as an environmentally-friendly agricultural chemical to help manage Bgt diseases [22,23]. However, the glycerol-mediated disease resistance pathways have not yet been clarified.
Transcriptome and proteomics technologies provide powerful tools for identification of genes and pathways associated with plant disease resistance. Earlier studies using Affymetrix wheat array reveal 3014 and 2800 wheat genes that are involved in the interaction of Bgt with susceptible cultivar Jingdong 8 and the resistant near-isogenic line carrying Pm30, respectively. The genes’ involvement in lipid degradation and metabolic processes, fatty acid catabolic process, and plant hormone signal transduction pathways have been identified [24]. Using iTRAQ-based quantitative proteomics, 394 proteins were identified in response to Bgt in a resistant wheat line, including PR polypeptides, oxidative stress responsive proteins, and primary metabolism pathways [25]. However, these proteins represent less than 10% of the differentially expressed genes that are identified at different time points post Bgt infection, detected by RNA sequencing [26].
In the current study, wheat plants treated with glycerol were resistant to Bgt, while water-treated plants were susceptible. Differentially expressed genes (DEGs) were identified by using RNAseq analysis. We found that glycerolipid metabolism and fatty acid biosynthesis pathways were in response to Bgt infection, which might contribute to G3P and OA18:1 accumulation. In addition, glycerol regulated hormone cross-talk (increased JA and SA levels, reduced IAA level) and induced PR proteins in wheat, are probably contributing to resistance to Bgt. This study contributes to a better understanding of the glycerol-mediated disease resistance pathways in wheat.

2. Results

2.1. Bgt Resistance Test in Water-Treated and Glycerol-Treated Wheat Leaves

Four groups of wheat plants were prepared for the identification of the transcriptional responses to glycerol treatment of wheat leaves, with and without Bgt infection: H0 and G0 groups represented water or glycerol-treated plants without Bgt infection; H24 and G24 groups represented water or glycerol-treated plants with Bgt infection (Figure 1A). The morphology of conidiophores and hyphae was monitored in water and glycerol-treated wheat leaves infected with Bgt (H24 and G24 groups). The obtained results showed that ~90% of the spores successfully germinated on water-treated leaves, while only about 10.7% of the spores germinated on the glycerol-treated leaves (Figure 1B,C). These results suggest that glycerol treatment induced effective resistance against Bgt infection. Therefore, changes in gene expression induced by both glycerol treatment and Bgt infection were explored using whole transcriptome analysis.

2.2. Transcriptome Sequencing, Assembly, and Functional Annotation

Transcriptome sequencing revealed a total of 681,705,618 clean reads (85.21 G clean data) from twelve cDNA libraries, constituting 41,214,380–77,987,478 clean reads (5.15–9.74 G clean data) for each pool, with a quality score of Q30 ≥89.12%. Following assembly, 81.99%~86.69% of the reads were mapped on the wheat genome reference sequence IWGSC RefSeq v1.0 (Table S1). In total, we annotated 106,755 distinct assembled unigenes after blast searching in six databases (Table S2, Figure 1D).

2.3. Identification of Differentially Expressed Genes (DEGs)

In the water-treated group (H24 vs. H0), 5616 DEGs were identified of which 3489 were upregulated and 2127 were downregulated in response to Bgt. In the glycerol-treated group (G24 vs. G0), only 2494 DEGs were identified, of which 1723 were upregulated and 771 were downregulated in response to Bgt. A total of 414 DEGs (376 up, 38 down) were identified in the G0 vs. H0 groups, and 463 DEGs (308 up; 155 down) were identified in the H24 vs. G24 groups (Table 1). A total of 1317 upregulated and 467 downregulated DEGs were overlapped in H24 vs. H0 and G24 vs. G0 groups, respectively. More than half of the DEGs in G0 vs. H0 groups overlapped with G24 vs. G0 or H24 vs. H0 groups (Figure 1D). DEGs identified in different biological replicates were clustered together in a heat map of expression levels, indicating good reproducibility between replicates (Figure S1). The annotated DEGs of different groups (H24 vs. H0, G24 vs. G0, G0 vs. H0, G24 vs. H24) are attached in Table S3.

2.4. Expression Patterns in Response to Bgt Infection (H24 vs. H0 and G24 vs. G0)

To better understand the response to Bgt infection in both water- and glycerol-treated plants, we annotated DEGs into Gene Ontology (GO) terms and KEGG pathways (Tables S4 and S5). DEGs were mainly classified into cell killing, immune system process, response to stimulus, signaling, and metabolic process GOs (e.g., biological process terms, Table S4). Those upregulated DEGs both in H24 vs. H0 and G24 vs. G0 groups showed significant enrichment in plant–pathogen interaction, sphingolipid metabolism, and phenylalanine metabolism pathways (Figure 2A,B); while the linoleic acid metabolism pathway was only markedly enriched in upregulated DEGs of H24 vs. H0 groups (Figure 2A). However, those downregulated DEGs both in H24 vs. H0 and G24 vs. G0 groups showed significant enrichment in photosynthesis-antenna proteins, steroid biosynthesis, fatty acid elongation, starch, and sucrose metabolism, and sesquiterpenoid and triterpenoid biosynthesis showed significant enrichment (Figure 2C,D). KEGG-enriched analysis showed that 36 DEGs in H24 vs. H0 groups and 22 DEGs in G24 vs. G0 groups belong to plant hormone signal transduction pathways. Those DEGs are involved in the signaling pathways of auxin (AUX/IAA), cytokinin (AHP and A-ARR), abscisic acid (PP2C and SnPK2), ethylene (ETR, EIN2, EIN3, and ERF1/2), and salicylic acid (NPR1, PR-1) (Figures S2 and S3). Overall, 66 DEGs in H24 vs. H0 groups and 40 DEGs in G24 vs. G0 groups were identified in the plant–fungus interaction pathway. Among them, PPRs (CEBiP), FLS2, MKK1/2, WRKY33, pti1, NHO1, and PR1 were induced by Bgt in the H24 vs. H0 groups, and CaMCML, MKK1/2, WRKY33, NHO1, and PR1 were induced by Bgt in the G24 vs. H0 groups, which might be associated with the PAMP-triggered immunity. Moreover, RPS2 and HSP90 were induced by Bgt in H24 vs. H0 groups, and RPM1, RPS,2 and HSP90 genes were induced in response to Bgt in G24 vs. G0 groups, which may be involved in effector-triggered immunity (Figures S4 and S5).

2.5. Changes in Expression Patterns of Glycerolipid Metabolism and Fatty Acid Biosynthesis Pathways in Response to Bgt Infection

We analyzed the DEGs of glycerolipid metabolism and fatty acid biosynthesis pathways, which responded to Bgt infection in H24 vs. H0. We found that 13 genes from 265 genes which were involved in glycerolipid metabolism pathways were regulated by Bgt infection (Figure 3, Tables S5 and S6). Among them, three TaGLI1 homologous genes (2.7.1.30), identified as glycerol kinase from A, B, and D genomes of wheat, were upregulated by Bgt infection. In addition, 1.1.1.21 identified as an aldehyde reductase, which transfers D-glycer-aldehyde to glycerol, was also upregulated by Bgt infection. While, 2.3.1.15 (named as TaACT1) identified as glycerol-3-phosphate acyltransferase was downregulated by Bgt infection. Moreover, 11 DEGs from 265 genes which were involved in fatty acid biosynthesis pathways were regulated by Bgt infection (Figure 3, Tables S5 and S6). We found that three TaSSI2 genes (1.14.192) from A, B, and D wheat genome, and some upstream genes identified as FabI, FabZ, FabF, and FabD were induced by Bgt infection (Figure 4, Table S6). Moreover, those genes were also regulated by Bgt in G24 vs. G0 groups and showed the same expression pattern with H24 vs. H0 groups (Table S6). The regulation of those genes in response to Bgt infection might contribute to the accumulation of G3P and several fatty acids. Moreover, it was demonstrated in our previous study that the levels of G3P, OA18:1, linoleic acid (LA18:2), and α-linolenic acid (ALA18:3) are induced by Bgt infection in wheat leaves [23].

2.6. Expression Patterns in Response to Glycerol Application (G0 vs. H0)

A total of 414 DEGs (376 up, 38 down) were identified in G0 vs. H0 groups in response to glycerol. Among them, 77 upregulated DEGs and seven downregulated DEGs were overlapped with those identified in H24 vs. H0 and G24 vs. G0 groups (Figure 1D, Tables S7 and S8). Those 414 DEGs were annotated in some GO terms such as cell killing, immune system process, and response to stimulus (Table S4). GO enrichment analysis revealed some important GO terms, including response to JA, response to SA, defense response to bacterium, wounding, lipid oxidation, growth, lateral root formation, oxylipin metabolic process, linoleate 9S-lipoxygenase activity, allene oxide synthase activity, linoleate 13S-lipoxygenase activity, and hydroperoxide dehydratase activity (Figure 5A). KEGG pathway analysis revealed that 21.1% (19/90) of genes in the linoleic acid metabolism, 2.3% (3/131) genes in cutin, suberine, and wax biosynthesis, and 1.9% (5/265) genes in glycerolipid metabolism could be regulated by glycerol (Table S5). KEGGenrichment analysis showed that the most highly enriched pathways were linoleic acid metabolism and alpha-linolenic acid metabolism (Figure 5B).
The following DEGs induced by glycerol were found in the linoleic acid metabolism pathway: 16 DEGs identified as linoleate 9S-lipoxygenase (LOX1_5, 1.13.11.58), and three DEGs identified as lipoxygenase (LOX2S, 1.13.11.12) (Figure S6, Table 2). In the alpha-linolenic acid metabolism pathway: two DEGs identified as lipoxygenase (LOX2S, 1.13.11.12), seven DEGs identified as allene oxide synthase (hydroperoxide dehydratase, AOS, 4.2.1.92), and two DEGs identified as 12-oxophytodienoic acid reductase (OPRs, 1.3.1.42) were induced by glycerol (Figure S7, Table 2). Among them, the AOS family genes were the most significantly-induced genes from 8.5–20.7 folds by glycerol application, and most were also induced by Bgt infection. Traes_2DL_FE7B99D58, identified as probable linoleate 9S-lipoxygenase 5, showed higher expression than other homologous gene members, and it was induced by both glycerol and Bgt (Table 2).
Glycerol-induced plant hormone signal transduction pathway was indicated by the upregulation of AUX/IAA, PP2C, and JAZ, annotated in the IAA, ABA, and JA pathways (Figure S8). Glycerol also induced plant–pathogen interaction pathways, such as HSP90 (two DEGs) (Figure S9). In addition, many pathogenesis-related proteins were induced by glycerol, such as genes encoding PR10, PR1, PRMS, disease resistance protein RPM1, chitinase 1, chitinase 8, HSP70, and almost all of those PR genes were also induced by Bgt (Table 2 and Table S3). KEGG enrichment analysis showed that down-regulated DEGs by glycerol were enriched in those pathways of photosynthesis and monoterpenoid biosynthesis (Figure S10).

2.7. Expression Patterns and DEGs in G24 vs. H24 Groups

A total of 308 DEGs were upregulated and 155 DEGs were downregulated in the G24 vs. H24 groups. Among them, 120 upregulated DEGs and two downregulated DEGs overlapped with G0 vs. H0 groups (Figure 6A, Table S9). The most enriched GO terms of those 308 DEGs in G24 vs. H24 groups were similar to G0 vs. H0 groups, such as response to jasmonic acid, oxylipin metabolic process, lipid oxidation, linoleate 9S-lipoxygenase activity, allene oxide synthase activity, linoleate 13S-lipoxygenase activity, and hydroperoxide dehydratase (Figure 6B). The two pathways of linoleic acid metabolism and alpha-linolenic acid metabolism also was identified from the enriched KEGG pathways in the G24 vs. H24 groups (Figure 6C).
In the alpha-linolenic acid metabolism pathway, 1.13.11.12 identified as a lipoxygenase (LOX2S), eight DEGs identified as allene oxide synthase (AOS, 4.2.1.92), and 5.3.99.6 identified as an allene oxide cyclase (AOC) were upregulated in the G24 pool compared with H24 pool. In the linoleic acid metabolism pathway (ko00591), six DEGs identified as linoleate 9S-lipoxygenase (LOX1_5, 1.13.11.58), and 1.13.11.12 identified as a lipoxygenase (LOX2S), were upregulated in the G24 pool compared with H24 pool (Table S10).
In the plant hormone signal transduction pathway, SnPK2, and JAZ (two DEGs), were upregulated in the G24 pool compared with H24 pool (Figure S11). PR-1 was downregulated in the G24 pool. In the plant–pathogen interaction pathway, HSP90 (two DEGs) and FLS2 were upregulated in the G24 pool (Figure S12). Moreover, those genes encoding peroxygenase, chitinase 8, peroxidase, beta-glucosidase, UDP-glycosyltransferase, glucan endo-1,3-beta-glucosidase, and xylanase inhibitor were upregulated in the G24 pool (Tables S10 and S3). We found that most of AOC, LOX, and PR genes were also induced in the G0 pool compared with H0 pool (Table S10). It suggested that those genes were significantly changed by glycerol application (G24 vs. H24; G0 vs. H0) and might take part in the resistance pathways in the G24 group.

2.8. Plant Hormones Cross-Talk in Response to Glycerol Application

AUX/IAA, PP2C, and JAZ were induced by glycerol, which was annotated in the IAA, ABA, and JA pathways, indicating activation of plant hormone signal transduction pathways (Figure S8). Moreover, the enriched GO terms in response to jasmonic acid and salicylic acid were identified in the G0 vs. H0 groups (Figure 5A). These results suggest that glycerol was involved in the regulation of SA, JA, IAA, and ABA pathways. Quantification of plant hormones validated glycerol induction of SA and JA levels, and reduced IAA level, while ABA level showed no change (Figure 7).

3. Discussion

Glycerol-induced resistance was reported in different plants to different pathogens as a typical non-host resistance, with upregulation of several PR genes and ROS accumulation [15,20,21,22]. ETI is the main mechanism for host resistance, with a race-specific resistance involving cell death. However, some overlapping defense signaling pathways or genes have been identified between non-host and host resistance, such as PR genes, oxidative burst-associated genes, and cell defense genes [27]. In this study, to better understand the mechanism, the whole transcriptome analysis of wheat plants infected by Bgt after glycerol or water treatment was used to reveal extensive alteration in gene expression. Our results showed that glycerol induced some defense-related pathways or genes which associated with ROS accumulation, cell death, and PR genes, which might contribute to resistance to Bgt in wheat (Figure 8).
We identified genes involved in the accumulation of G3P and OA18:1, including the upregulation of three TaGLI1 and three TaSSI2 homologs and downregulation of TaACT1, in response to Bgt infection (Figure 3 and Figure 4; Table S5). Previous reports show that the expression of SSI2 and GLI1 genes are induced by powdery mildew both in wheat and barley [23]. G3P and OA18:1 are important signal molecules involved in plant resistance. Mutants defective in G3P synthesis are compromised in systemic acquired resistance (SAR), while exogenous G3P induces SAR [9,13]. Overexpressing GLI1 increases G3P levels and enhances resistance against bacterial disease in Arabidopsis [11]. Hence, the regulation of genes in the glycerol metabolism pathway, TaGLI1 and TaACT1, which was followed by Bgt infection, might contribute to G3P accumulation and induction of disease resistance (Figure 3). However, numerous studies showed that downregulation of OA18:1 improves plant resistance [14,15,28]. For example, the ssi2 mutants or RNAi-mediated knockdown of SSI2 show enhanced resistance to different pathogens, with a lower OA18:1 level. Our previous results showed that OA18:1, LA18:2, and ALA18:3 levels were induced by Bgt infection [23]. It was demonstrated that unsaturated fatty acids (OA18:1, LA18:2, and ALA18:3) can inhibit innate immunity responses related to callose deposition in wheat [29]. Plant fatty acids can also be transferred to the pathogenic fungus and are required for colonization by pathogens as a fungal organic carbon source [19]. In this situation, the role of upregulated TaSSI2, FabI, FabZ, FabF, and FabD, following Bgt infection, might not contribute to wheat resistance (Figure 4). In addition, lipids encompassing fatty acids, fatty acid-based polymers, and fatty acid derivatives are part of the membrane structures of cells and tissues. Upon microbe infestation, the membrane is modified and membrane lipid defense-signaling molecules, including free fatty acids, oxylipins, JA, and the potent second messenger phosphatidic acid are released [30]. However, the role of membrane lipid is still unclear for disease resistance in plants, and the impact on membrane lipid after glycerol application also needs further research.
JA is derived from the fatty acid ALA18:3, and the first step is catalyzed by lipoxygenases (LOXs), as the formation of FA hydroperoxides, then formation of the catalyzed reaction by the AOS, AOC, and OPR3 [31,32]. In our analysis of glycerol application, linolenic acid and alpha-linolenic acid (ALA18:3) metabolism were the two most enriched KEGG pathways (Figure 5B). The upregulations of LOX2S, LOX1_5, AOS, and OPRs by glycerol application from those two pathways, might contribute to JA accumulation. Further support for these results is provided in the current study by quantification of plant hormones that showed JA accumulation induced by glycerol treatment. However, JA is involved in the defense against necrotrophic pathogens, preventing plant cell death, and inducing defense responses to restrict further pathogen infection [33]. Treatment with JA is shown to protect plants against herbivore attack and reduce the severity of infection by necrotrophic fungi [34,35]; while the application of MeJA does not induce resistance to powdery mildew (biotrophic fungi) in wheat [36]. Furthermore, we found the JASMONATE ZIM DOMAIN (JAZ) in JA signal transduction pathway was induced by glycerol. JAZ is a key target of microbes in their effort to suppress JA-dependent signaling [37]. Previous results suggest that tomato JAZ proteins regulate the progression of cell death during host and non-host interactions with different bacterial strains [38]. So, the upregulation of JAZ by glycerol might contribute to resistance to Bgt in wheat.
Furthermore, we found that glycerol also reduced IAA and induced SA levels. IAA induces the expression of expansins, proteins that loosen the cell wall, which is key for plant growth but may also make the plant vulnerable to biotic intruders. Exogenous application of IAA aggravates some pathogenic diseases progression [39,40]. Auxin-resistant Arabidopsis by mutation of genes functioning in auxin signaling shows enhanced disease resistance [41]. Previous reports showed that Aux/IAAs repressor proteins confer an improved disease resistance against pathogens, regulate the transcript levels of PRs, ROS accumulation, and callose development in the plant defense response, and act as transcriptional repressors of ARF genes to regulate downstream auxin-regulated genes [41,42,43]. Thus, the low level of IAA and upregulated Aux/IAAs in response to glycerol may improve plant wheat resistance to Bgt. SA is recognized as a very important SAR signal molecule, hence, the upregulation of SA by glycerol may also contribute to plant resistance to Bgt. SA causes global repression of auxin-related genes, including the TIR1 receptor gene, resulting in stabilization of the Aux/IAA repressor proteins and inhibition of auxin responses; the inhibitory effect on auxin signaling may be a part of the SA-mediated disease resistance mechanism [41]. Moreover, the amount of ABA was not changed after glycerol application, and other plant hormones-related pathways were not significantly changed with transcriptome analysis. Altogether, our analysis indicated that glycerol application altered plant hormone biosynthesis, signaling, and levels involved in Bgt resistance mechanisms.
Some pathogenesis related genes were induced by glycerol application such as encoding, PR10, PR1, chitinase-1, chitinase-8, peroxidase, glucan endo-1,3-beta-glucosidase, beta-glucosidase, UDP-glycosyltransferase, pathogenesis-related protein PRMS, disease resistance protein RPM1, heat shock 70 kDa protein, and heat shock 90 kDa protein. Previous reports showed that the downstream genes of SAR induce the expression of some PR protein genes, such as PR1, PR2, chitinase (PR3, PR8, and PR11), peroxidase (PR9), and oxalate oxidase (PR15 and PR16), and those PR genes can be used as potential candidate genes for improvement of the pathogen resistance of wheat and barley [44]. For example, PR genes encoding hydrolytic enzymes chitinases and β-1,3-glucanases, are very important in plants for invading pathogens, and overexpression of some PRs genes improves the resistance to different pathogens in different plants [45,46]. These results suggest that up-expression of those PRs by glycerol application might contribute to wheat resistance, and those PRs could be used as candidate genes for improving wheat resistance.

4. Material and Methods

4.1. Plant and Fungal Materials

The Bgt susceptible bread wheat cultivar Xuezao was used in this experiment. Powdery mildew isolate E09 was provided by Xiayu Duan, Institute of Plant Protection, Chinese Academy of Agricultural Sciences.

4.2. Plants Growth and Glycerol Treatments

Plants potted and maintained in a growth chamber with humidity of 75%, 26/20 °C day/night temperature regime, 12:12 h light/dark cycle, and light intensity of 40 µmol m−2 s−1. Wheat leaves at seedling stage (about one-week-old) were sprayed with 3% glycerol (containing 0.02% Silwett L-77) when the first leaves were fully expanded. Control plants were sprayed with water, containing 0.02% Silwett L-77. Glycerol and water treatments were preformed two times (once a day) until liquid was dropping from the leaves, to make sure sufficient treatments. Each treatment included three pots of plants, and each pot (12 × 12 cm) included 16 plants.

4.3. Pathogen Maintenance and Inoculation

Isolate E09 was maintained on the susceptible wheat line Xuezao through weekly transfer to new plants in the growth chamber with the same controls as described in Section 4.2. The fresh conidia were collected on the dark paper and papered before inoculation. After one day of glycerol or water treatments, wheat plants were put into a vaccination tower; conidia were slowly blown into the vaccination tower by using a blower with an estimated density of 100–150 conidia/mm2.

4.4. Coomassie Blue Staining

For microscopic observations of fungal development, Bgt-infected leaf segments were collected at 24 h post powdery mildew infection (hpi) for Coomassie blue staining as described previously [23]. Wheat leaf segments (5 cm in length) were stored in 50% glycerol and examined under an Olympus BX-43 microscope (Olympus Corporation, Tokyo, Japan). Three plants in each treatment and three leaves per plant were used for Comassie staining and the germination/penetration rates were calculated. The germination/penetration rates of conidiospores (number of germinated spores or penetrating spores relative to the total number of spores) were calculated and presented as the means of three independent replicated experiments. Statistical significance was determined by paired Student’s t test.

4.5. RNA Extraction, cDNA Library Construction, RNA-Seq, and Data Analysis

Four groups of wheat plants were prepared for RNA-seq: H0 and G0 groups represent water or glycerol-treated plants without Bgt infection; H24 and G24 groups represent water or glycerol-treated plants infected with Bgt one day after glycerol or water treatments. Each group included three biological repetitions. Ten leaves were pooled together at 24 hpi as one biological repetition for RNA extraction. The leaf samples of H0 and G0 groups were also collected at the same time point as controls of H24 and G24 groups at 24 hpi. The total RNA was extracted using RNA pure Plant Kit (Tiangen Biotech, Beijing, China). Sequencing libraries were generated using NEB Next UltraTM RNA Library Prep Kit for Illumina (NEB, Beijing, China) following manufacturer’s recommendations. The parried-end reads were generated on Illumina Hiseq 2500 platform.
Sequencing data were analyzed by using BMKCloud. Trimmomatic was used to remove the adaptors. TopHat2 tools software were used to map and assemble the reads to reference genome IWGSC RefSeq v1.0, [47]. Only reads with a perfect match or one mismatch were further analyzed and annotated based on the reference genome. Gene function was annotated based on the following databases: Nr (NCBI non-redundant protein sequences); Nt (NCBI non-redundant nucleotide sequences); Pfam (protein family); KOG/COG (clusters of orthologous groups of proteins); Swiss-Prot (a manually annotated and reviewed protein sequence database); KO (KEGG ortholog database); GO (Gene Ontology). Differential expression analysis of four groups was performed using the DESeq R package (1.10.1). Genes with p-value ≤ 0.05 were assigned as differentially expressed, with |log2 fold change| ≥1.8. GO enrichment analysis of the differentially expressed genes (DEGs) was implemented by the GOseq R packages based on the Wallenius non-central hyper-geometric distribution, which can adjust for gene length bias in DEGs [48]. KOBAS software was used to test the statistical enrichment of differential expression genes in KEGG pathways [49].

4.6. Measurements of Endogenous Plant Hormones

The leaf tissues of both the water treatment (control) and glycerol treatment were collected at the seeding stages one day after glycerol and water treatments. Plant hormones (ABA, IAA, JA, and SA) were extracted and measured from leaf tissue (approximately 50 mg of fresh weight) by simplified icELISA method as described previously [50], slightly modified by using a different antibody which was developed in the Wang Baomin’s lab, China Agricultural University.

5. Concluding Remarks

Bgt infection might contribute to G3P and OA18:1 accumulation in wheat. Glycerol application could improve G3P levels and reduced OA18:1 levels. Furthermore, glycerol altered plant hormones signaling pathways, which were confirmed by increase in JA and SA and reduction of IAA levels. Moreover, glycerol induced the expression of many PR proteins. These glycerol-induced changes probably contributed to wheat resistance to Bgt (Figure 8).

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/21/2/673/s1, Supporting Information items should be referred to in the text as follows: Supporting figures: Figures S1–S12. Supporting tables: Tables S1–S10.

Author Contributions

C.X., Y.L. and Q.S. conceived the project. Y.L. and L.Q. performed most of the experiments. X.L. provided help for using software and analysis of RNA-seq data. Q.Z. and X.Z. provided data curation. Y.L. wrote the manuscript. T.F., T.K. and C.X. improved and revised this manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Natural Science Foundation of China (31671676), the National Key Research and Development Project of China (2016YFD0101802), and the Israel Science Foundation (1366/18).

Acknowledgments

We are grateful to Wang Baomin for the analysis of plant hormones.

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. Mehta, Y.R. Wheat Diseases and Their Management; Springer: Cham, Switzerland, 2014; p. 256. [Google Scholar]
  2. Singh, R.P.; Singh, P.K.; Rutkoski, J.; Hodson, D.P.; He, X.; Jørgensen, L.N.; Hovmøller, M.S.; Huerta-Espino, J. Disease impact on wheat yield potential and prospects of genetic control. Annu. Rev. Phytopathol. 2016, 54, 303–322. [Google Scholar] [CrossRef] [Green Version]
  3. Jones, J.D.; Dangl, J.L. The plant immune system. Nature 2006, 444, 323. [Google Scholar] [CrossRef] [Green Version]
  4. Dangl, J.L.; Horvath, D.M.; Staskawicz, B.J. Pivoting the plant immune system from dissection to deployment. Science 2013, 341, 746–751. [Google Scholar] [CrossRef] [Green Version]
  5. Truman, W.; Bennett, M.H.; Kubigsteltig, I.; Turnbull, C.; Grant, M. Arabidopsis systemic immunity uses conserved defense signaling pathways and is mediated by jasmonates. Proc. Natl. Acad. Sci. USA 2007, 104, 1075–1080. [Google Scholar] [CrossRef] [Green Version]
  6. Truman, W.M.; Bennett, M.H.; Turnbull, C.G.; Grant, M.R. Arabidopsis auxin mutants are compromised in systemic acquired resistance and exhibit aberrant accumulation of various indolic compounds. Plant Physiol. 2010, 152, 1562–1573. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Ton, J.; Flors, V.; Mauch-Mani, B. The multifaceted role of ABA in disease resistance. Trends Plant Sci. 2009, 14, 310–317. [Google Scholar] [CrossRef] [PubMed]
  8. Kachroo, A.; Kachroo, P. Fatty Acid-derived signals in plant defense. Annu. Rev. Phytopathol. 2009, 47, 153–176. [Google Scholar] [CrossRef] [PubMed]
  9. Chanda, B.; Xia, Y.; Mandal, M.K.; Yu, K.; Sekine, K.T.; Gao, Q.M.; Selote, D.; Hu, Y.; Stromberg, A.; Navarre, D.; et al. Glycerol-3-phosphate is a critical mobile inducer of systemic immunity in plants. Nat. Genet. 2011, 43, 421–427. [Google Scholar] [CrossRef]
  10. Yang, Y.; Zhao, J.; Liu, P.; Xing, H.; Li, C.; Wei, G.; Kang, Z. Glycerol-3-phosphate metabolism in wheat contributes to systemic acquired resistance against Puccinia striiformis f. sp. tritici. PLoS ONE 2013, 8, e81756. [Google Scholar] [CrossRef]
  11. Kang, L.; Li, J.; Zhao, T.; Xiao, F.; Tang, X.; Thilmony, R.; He, S.; Zhou, J.M. Interplay of the Arabidopsis nonhost resistance gene NHO1 with bacterial virulence. Proc. Natl. Acad. Sci. USA 2003, 100, 3519–3524. [Google Scholar] [CrossRef] [Green Version]
  12. Chanda, B.; Venugopal, S.C.; Kulshrestha, S.; Navarre, D.A.; Downie, B.; Vaillancourt, L.; Kachroo, A.; Kachroo, P. Glycerol-3-phosphate levels are associated with basal resistance to the hemibiotrophic fungus Colletotrichum higginsianum in Arabidopsis. Plant Physiol. 2008, 147, 2017–2029. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Venugopal, S.C.; Chanda, B.; Vaillancourt, L.; Kachroo, A.; Kachroo, P. The common metabolite glycerol-3-phosphate is a novel regulator of plant defense signaling. Plant Signal Behav. 2009, 4, 746–749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Kachroo, P.; Shanklin, J.; Shah, J.; Whittle, E.J.; Klessig, D.F. A fatty acid desaturase modulates the activation of defense signaling pathways in plants. Proc. Natl. Acad. Sci. USA 2001, 98, 9448–9453. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Jiang, C.J.; Shimono, M.; Maeda, S.; Inoue, H.; Mori, M.; Hasegawa, M.; Sugano, S.; Takatsuji, H. Suppression of the rice fatty-acid desaturase gene OsSSI2 enhances resistance to blast and leaf blight diseases in rice. Mol. Plant Microbe Interact. 2009, 22, 820–829. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Savchenko, T.; Walley, J.W.; Chehab, E.W.; Xiao, Y.; Kaspi, R.; Pye, M.F.; Mohamed, M.E.; Lazarus, C.M.; Bostock, R.M.; Dehesh, K. Arachidonic acid: An evolutionarily conserved signaling molecule modulates plant stress signaling networks. Plant Cell 2010, 22, 3193–3205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Keymer, A.; Pimprikar, P.; Wewer, V.; Huber, C.; Brands, M.; Bucerius, S.L.; Delaux, P.M.; Klingl, V.; von Roepenack-Lahaye, E.; Wang, T.L.; et al. Lipid transfer from plants to arbuscular mycorrhiza fungi. Elife 2017, 6, 29107. [Google Scholar] [CrossRef]
  18. Luginbuehl, L.H.; Menard, G.N.; Kurup, S.; Van Erp, H.; Radhakrishnan, G.V.; Breakspear, A.; Oldroyd, G.E.; Eastmond, P.J. Fatty acids in arbuscular mycorrhizal fungi are synthesized by the host plant. Science 2017, 356, 1175–1178. [Google Scholar] [CrossRef] [Green Version]
  19. Jiang, Y.; Wang, W.; Xie, Q.; Liu, N.; Liu, L.; Wang, D.; Zhang, X.; Yang, C.; Chen, X.; Tang, D.; et al. Plants transfer lipids to sustain colonization by mutualistic mycorrhizal and parasitic fungi. Science 2017, 356, 1172–1175. [Google Scholar] [CrossRef] [Green Version]
  20. Kachroo, A.; Venugopal, S.C.; Lapchyk, L.; Falcone, D.; Hildebrand, D.; Kachroo, P. Oleic acid levels regulated by glycerolipid metabolism modulate defense gene expression in Arabidopsis. Proc. Natl. Acad. Sci. USA 2004, 101, 5152–5157. [Google Scholar] [CrossRef] [Green Version]
  21. Kachroo, A.; Fu, D.Q.; Havens, W.; Navarre, D.; Kachroo, P.; Ghabrial, S.A. An oleic acid-mediated pathway induces constitutive defense signaling and enhanced resistance to multiple pathogens in soybean. Mol. Plant Microbe Interact. 2008, 21, 564–575. [Google Scholar] [CrossRef] [Green Version]
  22. Zhang, Y.; Smith, P.; Maximova, S.N.; Guiltinan, M.J. Application of glycerol as a foliar spray activates the defense response and enhances disease resistance of Theobroma cacao. Mol. Plant Pathol. 2015, 16, 27–37. [Google Scholar] [CrossRef] [PubMed]
  23. Li, Y.; Song, N.; Zhao, C.; Li, F.; Geng, M.; Wang, Y.; Liu, W.; Xie, C.; Sun, Q. Application of glycerol for induced powdery mildew resistance in Triticum aestivum L. Front. Physiol. 2016, 7, 413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Xin, M.; Wang, X.; Peng, H.; Yao, Y.; Xie, C.; Han, Y.; Ni, Z.; Sun, Q. Transcriptome comparison of susceptible and resistant wheat in response to powdery mildew infection. Genom. Proteom. Bioinform. 2012, 10, 94–106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Fu, Y.; Zhang, H.; Mandal, S.N.; Wang, C.; Chen, C.; Ji, W. Quantitative proteomics reveals the central changes of wheat in response to powdery mildew. J. Proteom. 2016, 130, 108–119. [Google Scholar] [CrossRef]
  26. Zhang, H.; Yang, Y.; Wang, C.; Liu, M.; Li, H.; Fu, Y.; Wang, Y.; Nie, Y.; Liu, X.; Ji, W. Large-scale transcriptome comparison reveals distinct gene activations in wheat responding to stripe rust and powdery mildew. BMC Genom. 2014, 15, 898. [Google Scholar] [CrossRef] [Green Version]
  27. Oh, S.K.; Lee, S.; Chung, E.; Park, J.M.; Yu, S.H.; Ryu, C.M.; Choi, D. Insight into Types I and II nonhost resistance using expression patterns of defense-related genes in tobacco. Planta 2006, 223, 1101–1107. [Google Scholar] [CrossRef]
  28. Song, N.; Hu, Z.; Li, Y.; Li, C.; Peng, F.; Yao, Y.; Peng, H.; Ni, Z.; Xie, C.; Sun, Q. Overexpression of a wheat stearoyl-ACP desaturase (SACPD) gene TaSSI2 in Arabidopsis ssi2 mutant compromise its resistance to powdery mildew. Gene 2013, 524, 220–227. [Google Scholar] [CrossRef]
  29. Blümke, A.; Falter, C.; Herrfurth, C.; Sode, B.; Bode, R.; Schäfer, W.; Feussner, I.; Voigt, C.A. Secreted fungal effector lipase releases free fatty acids to inhibit innate immunity-related callose formation during wheat head infection. Plant Physiol. 2014, 165, 346–358. [Google Scholar] [CrossRef] [Green Version]
  30. Lee, H.J.; Park, O.K. Lipases associated with plant defense against pathogens. Plant Sci. 2019, 279, 51–58. [Google Scholar] [CrossRef]
  31. Wasternack, C. Jasmonates: An update on biosynthesis, signal transduction and action in plant stress response, growth and development. Ann. Bot. 2007, 100, 681–697. [Google Scholar] [CrossRef] [Green Version]
  32. Lim, G.H.; Singhal, R.; Kachroo, A.; Kachroo, P. Fatty acid-and lipid-mediated signaling in plant defense. Annu. Rev. Phytopathol. 2017, 55, 505–536. [Google Scholar] [CrossRef] [PubMed]
  33. Pandey, D.; Rajendran, S.R.C.K.; Gaur, M.; Sajeesh, P.K.; Kumar, A. Plant defense signaling and responses against necrotrophic fungal pathogens. J. Plant Growth Regul. 2016, 35, 1159–1174. [Google Scholar] [CrossRef]
  34. Baldwin, I.T. Jasmonate-induced responses are costly but benefit plants under attack in native populations. Proc. Natl. Acad. Sci. USA 1998, 95, 8113–8118. [Google Scholar] [CrossRef] [Green Version]
  35. Zalewski, K.; Lahuta, L.B.; Martysiak-Żurowska, D.; Okorski, A.; Nitkiewicz, B.; Zielonka, Ł. The effect of exogenous application of methyl jasmonate on the lipid and carbohydrate content and composition of winter triticale (Triticosecale Wittm.) grain and the severity of fungal infections in triticale plants and grain. J. Agric. Food Chem. 2019, 67, 5932–5939. [Google Scholar] [CrossRef] [PubMed]
  36. Xiang, Y.; Song, M.; Wei, Z.; Tong, J.; Zhang, L.; Xiao, L.; Ma, Z.; Wang, Y. A jacalin-related lectin-like gene in wheat is a component of the plant defence system. J. Exp. Bot. 2011, 62, 5471–5483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Guo, Q.; Yoshida, Y.; Major, I.T.; Wang, K.; Sugimoto, K.; Kapali, G.; Havko, N.E.; Benning, C.; Howe, G.A. JAZ repressors of metabolic defense promote growth and reproductive fitness in Arabidopsis. Proc. Natl. Acad. Sci. USA 2018, 115, 10768–10777. [Google Scholar] [CrossRef] [Green Version]
  38. Ishiga, Y.; Ishiga, T.; Uppalapati, S.R.; Mysore, K.S. Jasmonate ZIM-domain (JAZ) protein regulates host and nonhost pathogen-induced cell death in tomato and Nicotiana benthamiana. PLoS ONE 2013, 8, e75728. [Google Scholar] [CrossRef]
  39. Ding, X.; Cao, Y.; Huang, L.; Zhao, J.; Xu, C.; Li, X.; Wang, S. Activation of the indole-3-acetic acid–amido synthetase GH3-8 suppresses expansin expression and promotes salicylate-and jasmonate-independent basal immunity in rice. Plant Cell 2008, 20, 228–240. [Google Scholar] [CrossRef] [Green Version]
  40. Fu, J.; Liu, H.; Li, Y.; Yu, H.; Li, X.; Xiao, J.; Wang, S. Manipulating broad-spectrum disease resistance by suppressing pathogen-induced auxin accumulation in rice. Plant Physiol. 2011, 155, 589–602. [Google Scholar] [CrossRef] [Green Version]
  41. Wang, D.; Pajerowska-Mukhtar, K.; Culler, A.H.; Dong, X. Salicylic acid inhibits pathogen growth in plants through repression of the auxin signaling pathway. Curr. Biol. 2007, 17, 1784–1790. [Google Scholar] [CrossRef] [Green Version]
  42. Fousia, S.; Tsafouros, A.; Roussos, P.A.; Tjamos, S.E. Increased resistance to Verticillium dahliae in Arabidopsis plants defective in auxin signalling. Plant Pathol. 2018, 67, 1749–1757. [Google Scholar] [CrossRef]
  43. Fan, S.; Chang, Y.; Liu, G.; Shang, S.; Tian, L.; Shi, H. Molecular functional analysis of auxin/indole-3-acetic acid proteins (Aux/IAAs) in plant disease resistance in cassava. Physiol. Plant. 2019. [Google Scholar] [CrossRef] [PubMed]
  44. Wang, X.; Bi, W.; Jing, G.; Yu, X.; Wang, H.; Liu, D. Systemic acquired resistance, NPR1, and pathogenesis-related genes in wheat and barley. J. Integr. Agric. 2018, 17, 2468–2477. [Google Scholar] [CrossRef] [Green Version]
  45. Ebrahim, S.; Usha, K.; Singh, B. Pathogenesis related (PR) proteins in plant defense mechanism. Sci. Against Microb. Pathog. 2011, 2, 1043–1054. [Google Scholar]
  46. Ali, S.; Ganai, B.A.; Kamili, A.N.; Bhat, A.A.; Mir, Z.A.; Bhat, J.A.; Tyagi, A.; Islam, S.T.; Mushtaq, M.; Yadav, P.; et al. Pathogenesis-related proteins and peptides as promising tools for engineering plants with multiple stress tolerance. Microbiol. Res. 2018, 212, 29–37. [Google Scholar] [CrossRef] [PubMed]
  47. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14, 36. [Google Scholar] [CrossRef] [Green Version]
  48. 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]
  49. Mao, X.; Cai, T.; Olyarchuk, J.G.; Wei, L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 2005, 21, 3787–3793. [Google Scholar] [CrossRef]
  50. Zhao, J.; Li, G.; Yi, G.X.; Wang, B.M.; Deng, A.X.; Nan, T.G.; Li, Z.H.; Li, Q.X. Comparison between conventional indirect competitive enzyme-linked immunosorbent assay (icELISA) and simplified icELISA for small molecules. Anal. Chim. Acta 2006, 571, 79–85. [Google Scholar] [CrossRef]
Ijms 21 00673 g001 550
Figure 1. Application of glycerol enhanced Bgt resistance in wheat and transcriptome analysis. (A) Experimental design. (B) Morphology of conidiospore and hyphae in untreated, water-treated, and glycerol-treated Xuezao leaves at 24 h post powdery mildew infection (hpi). (C) Germination rates of powdery mildew on the water-treated and glycerol-treated wheat leaves at 24 hpi. Each value is the mean ± standard error (SE) of three independent biological repetitions. Asterisk indicates a significant difference from the germination rates of treated leaves at p ≤ 0.05 by Student’s t-test. (D) Differentially expressed genes (DEGs) overlapped in different groups comparison. H0 and G0 groups presented water or glycerol-treated plants without Bgt infection; H24 and G24 group presented water or glycerol-treated plants with Bgt infection.
Figure 1. Application of glycerol enhanced Bgt resistance in wheat and transcriptome analysis. (A) Experimental design. (B) Morphology of conidiospore and hyphae in untreated, water-treated, and glycerol-treated Xuezao leaves at 24 h post powdery mildew infection (hpi). (C) Germination rates of powdery mildew on the water-treated and glycerol-treated wheat leaves at 24 hpi. Each value is the mean ± standard error (SE) of three independent biological repetitions. Asterisk indicates a significant difference from the germination rates of treated leaves at p ≤ 0.05 by Student’s t-test. (D) Differentially expressed genes (DEGs) overlapped in different groups comparison. H0 and G0 groups presented water or glycerol-treated plants without Bgt infection; H24 and G24 group presented water or glycerol-treated plants with Bgt infection.
Ijms 21 00673 g001
Ijms 21 00673 g002 550
Figure 2. The top 20 enriched KEGG pathways of upregulated (in H24) DEGs in the (A) H24 vs. H0 and (B) G24 vs. G0 groups. The top 20 enriched KEGG pathways of downregulated DEGs (in H24) in the (C) H24 vs. H0 and (D) G24 vs. G0 groups. The color represents the Q-value as shown in the legend. Q-values are the p-values corrected for multiple hypothesis testing and range from 0 to 1. The closer the Q-value is to zero, the more significant the enrichment. The horizontal axis indicates the rich factor, which means that the ratio of the DEGs number and the number of genes are annotated in this pathway. The greater the rich factor, the greater degree of enrichment. The size of each circle indicates the number of DEGs in that pathway. The larger the circle, the greater number of DEGs.
Figure 2. The top 20 enriched KEGG pathways of upregulated (in H24) DEGs in the (A) H24 vs. H0 and (B) G24 vs. G0 groups. The top 20 enriched KEGG pathways of downregulated DEGs (in H24) in the (C) H24 vs. H0 and (D) G24 vs. G0 groups. The color represents the Q-value as shown in the legend. Q-values are the p-values corrected for multiple hypothesis testing and range from 0 to 1. The closer the Q-value is to zero, the more significant the enrichment. The horizontal axis indicates the rich factor, which means that the ratio of the DEGs number and the number of genes are annotated in this pathway. The greater the rich factor, the greater degree of enrichment. The size of each circle indicates the number of DEGs in that pathway. The larger the circle, the greater number of DEGs.
Ijms 21 00673 g002
Ijms 21 00673 g003 550
Figure 3. Gene expressions in the glycerolipid metabolism pathway (Ko00561) in response to Bgt infection. The red color indicates gene expression induction in H24 group; the green color indicates gene expression reduction in H24 group; the blue color indicates that the expression patterns of some annotated genes were induced and some homologous genes were reduced in H24 group.
Figure 3. Gene expressions in the glycerolipid metabolism pathway (Ko00561) in response to Bgt infection. The red color indicates gene expression induction in H24 group; the green color indicates gene expression reduction in H24 group; the blue color indicates that the expression patterns of some annotated genes were induced and some homologous genes were reduced in H24 group.
Ijms 21 00673 g003
Ijms 21 00673 g004 550
Figure 4. Gene expressions in the fatty acid biosynthesis pathway (Ko00061) in response to Bgt infection. The red color indicates gene expression induction by Bgt; the green color indicates gene expression reduction by Bgt.
Figure 4. Gene expressions in the fatty acid biosynthesis pathway (Ko00061) in response to Bgt infection. The red color indicates gene expression induction by Bgt; the green color indicates gene expression reduction by Bgt.
Ijms 21 00673 g004
Ijms 21 00673 g005 550
Figure 5. The top 20 enriched GO terms and KEGG pathways. (A) GO enrichment analysis of genes differentially expressed in response to glycerol. Data are presented according to the p-value. The names of the 20 most highly enriched GO terms are arranged on the vertical axis according to the −log 10 (p-value); the horizontal axis represents the −log 10 (p-value). ** indicates highly significant differences from H0 group at p ≤ 0.01 by Student’s t-test. (B) Statistical scatterplot showing the pathways enriched in genes differentially regulated in response to glycerol. The color represents the Q-value as shown in the legend. Q-values are p-values corrected for multiple hypothesis testing and range from 0 to 1. The closer the Q-value is to zero, the more significant the enrichment. The horizontal axis indicates the rich factor, which means that the ratio of the DEGs number and the number of genes are annotated in this pathway. The greater the rich factor, the greater degree of enrichment. The size of each circle indicates the number of DEGs in that pathway. The bigger the circle, the greater number of DEGs.
Figure 5. The top 20 enriched GO terms and KEGG pathways. (A) GO enrichment analysis of genes differentially expressed in response to glycerol. Data are presented according to the p-value. The names of the 20 most highly enriched GO terms are arranged on the vertical axis according to the −log 10 (p-value); the horizontal axis represents the −log 10 (p-value). ** indicates highly significant differences from H0 group at p ≤ 0.01 by Student’s t-test. (B) Statistical scatterplot showing the pathways enriched in genes differentially regulated in response to glycerol. The color represents the Q-value as shown in the legend. Q-values are p-values corrected for multiple hypothesis testing and range from 0 to 1. The closer the Q-value is to zero, the more significant the enrichment. The horizontal axis indicates the rich factor, which means that the ratio of the DEGs number and the number of genes are annotated in this pathway. The greater the rich factor, the greater degree of enrichment. The size of each circle indicates the number of DEGs in that pathway. The bigger the circle, the greater number of DEGs.
Ijms 21 00673 g005
Ijms 21 00673 g006 550
Figure 6. The DEGs in G24 vs. H24 groups. (A) DEGs overlapped in different groups comparison. (B) GO enrichment analysis of DEGs in G24 vs. H24 groups. (C) Statistical scatterplot showing the pathways enriched in genes differentially regulated in response to glycerol in the G24 vs. H24 groups.
Figure 6. The DEGs in G24 vs. H24 groups. (A) DEGs overlapped in different groups comparison. (B) GO enrichment analysis of DEGs in G24 vs. H24 groups. (C) Statistical scatterplot showing the pathways enriched in genes differentially regulated in response to glycerol in the G24 vs. H24 groups.
Ijms 21 00673 g006
Ijms 21 00673 g007 550
Figure 7. Content of abscisic acid (ABA), auxin (IAA), jasmonic acid (JA), and salicylic acid (SA) levels in water-treated and glycerol-treated Xuezao leaves. Each value is the mean ± SE of three independent biological repetitions. Asterisks indicate a significant difference from the water-treated mock control at p ≤ 0.05 determined by Student’s t-test.
Figure 7. Content of abscisic acid (ABA), auxin (IAA), jasmonic acid (JA), and salicylic acid (SA) levels in water-treated and glycerol-treated Xuezao leaves. Each value is the mean ± SE of three independent biological repetitions. Asterisks indicate a significant difference from the water-treated mock control at p ≤ 0.05 determined by Student’s t-test.
Ijms 21 00673 g007
Ijms 21 00673 g008 550
Figure 8. The glycerol-induced resistance pathway. Upon Bgt infection, the TaGLI1 and TaSSI2 were induced, and TaACT1 was reduced, which might contribute to G3P and OA18:1 accumulation. G3P accumulation might be an active defense process; OA18:1 accumulation might be not good for plant resistance. Exogenous application of glycerol increased the G3P levels and decreased the level of OA18:1. Glycerol induced some genes involved in JA synthesis (e.g., LOX, AOS, and OPRs), then induced the JA level. Glycerol reduced IAA and induced SA levels with cross-talk of plant hormones. Glycerol induced ROS accumulation, PR and Hsp90 genes expression, and cell death, then contributed resistance to Bgt in wheat.
Figure 8. The glycerol-induced resistance pathway. Upon Bgt infection, the TaGLI1 and TaSSI2 were induced, and TaACT1 was reduced, which might contribute to G3P and OA18:1 accumulation. G3P accumulation might be an active defense process; OA18:1 accumulation might be not good for plant resistance. Exogenous application of glycerol increased the G3P levels and decreased the level of OA18:1. Glycerol induced some genes involved in JA synthesis (e.g., LOX, AOS, and OPRs), then induced the JA level. Glycerol reduced IAA and induced SA levels with cross-talk of plant hormones. Glycerol induced ROS accumulation, PR and Hsp90 genes expression, and cell death, then contributed resistance to Bgt in wheat.
Ijms 21 00673 g008
Table
Table 1. Comparison of the number of transcripts identified in the DEGs.
Table 1. Comparison of the number of transcripts identified in the DEGs.
DEG SetAll DEGsUpregulatedDownregulated
H24 vs. H0561634892127
G24 vs. G024941723771
G0 vs. H041437638
G24 vs. H24463308155
Table
Table 2. DEGs in the G0 vs. H0 groups.
Table 2. DEGs in the G0 vs. H0 groups.
Gene #IDH0G0H24G24Annotation
1.13.11.12Traes_5BS_06078574032.93123.1246.6897.83Lipoxygenase 2.1
Traes_5DS_E8892706A93.92330.19155.71290.37
Traes_6DS_7CA5A8F1241.6273.3730.5138.32Lipoxygenase 2.3
4.2.1.92Traes_4AS_41FB87D392.3721.675.5624.73Allene oxide synthase 2
Traes_4AS_9F1B2A7DD1.6614.561.007.24
Traes_4BL_523D155E212.6343.796.7949.98
Traes_4BL_C01E043B61.1010.831.467.70
Traes_4BL_DD6DD748717.33119.9138.96170.11
Traes_4BL_EC6D200265.55115.0829.68122.97
Traes_4DL_B3E978E9F1.7614.961.026.22
1.3.1.42Traes_6DL_94DCF0B700.291.182.482.3512-oxophytodienoate reductase
Traes_7BS_62CC4CA590.636.320.192.01
1.13.11.58Traes_2DL_17913EE21109.35358.45275.74312.34Linoleate 9S-lipoxygenase 1
Traes_2DL_CE85DC5C035.10133.1893.19117.98
Traes_2DL_D4BCDAA7624.1270.5541.8565.25
Traes_4BS_63DD9D0362.135.625.7510.62
Traes_4BS_71CB57A0D4.3315.908.6017.88
Traes_4BS_939C791841.459.814.199.49
Traes_4BS_9DDF3D7C62.056.425.439.57
Traes_4DS_7868A8C2E2.1311.587.1618.30
Traes_2DL_B5B62EE1131.1695.4265.1585.79Probable linoleate 9S-Lipoxygenase 5
Traes_2BL_77148B8D818.7276.3638.6561.36
Traes_2DL_FE7B99D583700210951654415919
Traes_6AS_9557563D18.1770.756.1088.52Putative linoleate 9S-Lipoxygenase 3
Traes_6BS_B26FD03C87.0862.704.0141.90
Triticum_aestivumLinn_newGene_117321.4413.990.6817.84
Traes_5BL_304FAFA260.091.780.153.64Linoleate 9S-lipoxygenase 2
Traes_2AL_5BAB2682727.3975.7478.1076.00Seed linoleate 9S-lipoxygenase-3
1.11.2.3Traes_2AL_6A8D574C40.323.420.252.85Peroxygenase
Traes_2BL_5880DAAC30.111.310.000.41
Traes_2DL_00618F3150.091.910.100.83
AUX/IAATriticum_aestivumLinn_newGene_1260.701.470.871.62Auxin-responsive protein IAA18
PP2CTraes_4BL_3403452C00.612.360.300.88Serine/threonine protein phosphatase 2C 30
JAZTriticum_aestivumLinn_newGene_185790.501.320.520.92Protein TIFY 10B
HSPTraes_2AS_67EFE0FAE0.241.530.290.79Heat shock protein 90
Traes_2DS_3B16D81730.091.580.420.83
Traes_1AL_51CED3DBF4.027.935.778.92Heat shock cognate 70 kDa protein 4
Traes_1AS_5DE9A16CD0.861.902.532.77Heat shock 70 kDa protein 4L
Traes_3B_B67388A960.331.890.350.87Heat shock cognate 70 kDa protein
Traes_3B_FB10B725B0.542.530.461.36Heat shock cognate 70 kDa protein
Traes_5BL_C318204D21.583.142.163.18Hsp70 nucleotide exchange factor fes1-like
PR-1Traes_3DL_85AC9E60D0.7122.84152.38109.17Glucan endo-1,3-beta-glucosidase GII
Traes_3DL_F5930F58D2.1336.99344.10210.39
Traes_3DL_48F92563F5.7918.03258.95136.05Glucan endo-1,3-beta-glucosidase GIII
Triticum_aestivumLinn_newGene_52030.922.3235.7822.19beta-1,3-glucanase precursor
Traes_3AL_BCFD5F3030.5416.40152.58103.74beta-1,3-glucanase
Traes_3B_9F3320C780.123.5938.2126.41Glucan endo-1,3-beta-D-glucosidase
Traes_5DL_ED441B7EB1.623.543.345.37Cell wall beta-glucosidase
Triticum_aestivumLinn_newGene_16520.131.5035.0717.75Glucan endo-1,3-beta-glucosidase 13
Traes_1BL_A6F7A9A540.601.390.451.03Glucan endo-1,3-beta-glucosidase 14
PR-3Traes_7AL_FAE816A853.6829.0060.3180.11Chitinase 1
Traes_7DL_24AA718601.5511.6728.8735.39
Traes_1AL_E96C0662D8.9715.92297.95169.09Chitinase 8
Traes_5BL_DE0C53CE21.4218.7627.6025.49
Traes_5DL_EEF38A7E40.233.516.4113.68
Callose synthaseTraes_3DL_9021F9E750.511.121.111.07Callose synthase 12
Traes_6BL_2D48C932A0.130.410.420.16Callose synthase 3
Traes_7DS_B481462CF0.561.191.161.06Callose synthase 10
PeroxidaseTraes_2AS_4576043590.444.4467.4150.64Peroxidase 2
Traes_2BS_EAB2C09D00.141.8835.6236.93
Traes_2DS_708F03DA320.2550.88301.61183.45
Traes_3DL_B38DFDDFF0.3911.672.5422.57
Traes_2DS_D76AB139C0.251.08141.2660.39Peroxidase 1
Traes_6BS_0BDACE20520.8842.56197.18181.62Peroxidase 6
Traes_5BL_36EBD512B0.200.940.771.91Peroxidase 35
Disease-related geneTraes_2BL_19C5224BE0.412.350.411.24Pathogenesis-related protein 1
Traes_2BL_B657F7F3A0.001.430.120.77
Traes_4AS_C5AE1BBDD2.145.2716.468.89
Traes_4DL_5C688784F3.189.0022.6116.89
Traes_6DL_F504536D54.649.1317.0224.47Putative disease resistance RGA4
Traes_7BS_94EB3B3D611.5442.41484.99343.31Pathogenesis-related protein PRMS
Triticum_aestivumLinn_newGene_118332.825.116.739.12Putative disease resistance RPM1
Triticum_aestivumLinn_newGene_229340.371.801.283.06Putative LRR receptor-like serine/threonine-protein kinase

Share and Cite

MDPI and ACS Style

Li, Y.; Qiu, L.; Liu, X.; Zhang, Q.; Zhuansun, X.; Fahima, T.; Krugman, T.; Sun, Q.; Xie, C. Glycerol-Induced Powdery Mildew Resistance in Wheat by Regulating Plant Fatty Acid Metabolism, Plant Hormones Cross-Talk, and Pathogenesis-Related Genes. Int. J. Mol. Sci. 2020, 21, 673. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21020673

AMA Style

Li Y, Qiu L, Liu X, Zhang Q, Zhuansun X, Fahima T, Krugman T, Sun Q, Xie C. Glycerol-Induced Powdery Mildew Resistance in Wheat by Regulating Plant Fatty Acid Metabolism, Plant Hormones Cross-Talk, and Pathogenesis-Related Genes. International Journal of Molecular Sciences. 2020; 21(2):673. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21020673

Chicago/Turabian Style

Li, Yinghui, Lina Qiu, Xinye Liu, Qiang Zhang, Xiangxi Zhuansun, Tzion Fahima, Tamar Krugman, Qixin Sun, and Chaojie Xie. 2020. "Glycerol-Induced Powdery Mildew Resistance in Wheat by Regulating Plant Fatty Acid Metabolism, Plant Hormones Cross-Talk, and Pathogenesis-Related Genes" International Journal of Molecular Sciences 21, no. 2: 673. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21020673

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