Comparative Proteomic Analysis of the Defense Response to Gibberella Stalk Rot in Maize and Reveals That ZmWRKY83 Is Involved in Plant Disease Resistance

Fusarium graminearum is the causal agent of Gibberella stalk rot in maize stem, resulting in maize lodging, yield, quality, and mechanical harvesting capacity. To date, little is known about the maize stem defense mechanism in response to the invasion of F. graminearum. This study represents a global proteomic approach to document the infection by F. graminearum. A total of 1,894 differentially expressed proteins (DEPs) were identified in maize stem with F. graminearum inoculation. Functional categorization analysis indicated that proteins involved in plant-pathogen interaction were inducible at the early stages of infection. We also found that the expression of proteins involved in phenylpropanoid, flavonoid, and terpenoid biosynthesis were upregulated in response to F. graminearum infection, which may reflect that these secondary metabolism pathways were important in the protection against the fungal attack in maize stem. In continuously upregulated proteins after F. graminearum infection, we identified a WRKY transcription factor, ZmWRKY83, which could improve the resistance to plant pathogens. Together, the results show that the defense response of corn stalks against F. graminearum infection was multifaceted, involving the induction of proteins from various immune-related pathways, which had a directive significance for molecular genetic breeding of maize disease-resistant varieties.


INTRODUCTION
Plants are stressed by various pathogens in nature. However, few pathogens can be successfully colonized in host plants, which indicates the existence of recognition and defense mechanisms. Plants have evolved multiple defense strategies for combating invading pathogens. Pathogenassociated molecular patterns (PAMPs)-triggered immunity (PTI) and effector-triggered immunity (ETI) are well-defined modes of plant immunity against pathogens (Thomma et al., 2011). Infected organisms recognize that they are under attack by detecting pathogens directly through PAMPs that bind to pattern recognition receptors. The infected host can identify pathogens by ETI pattern recognition, thereby alerting the host to pathogens through associated damage caused by pathogenic toxins or effectors (Akira et al., 2006;Kazan and Lyons, 2014;Cui et al., 2015). Reactive oxygen species production, MAP kinase activation, which is related to immune responses, can occur in both PTI and ETI (Tsuda and Katagiri, 2010).
Maize is one of the most important crops worldwide. Unfortunately, maize is susceptible to a variety of pathogens during growth and development. Among the fungal pathogens, Fusarium graminearum is the causative agent of ear and stalk rot disease in maize (Santiago et al., 2007). In maize, when the stalks die prematurely, the plant produces light ears with poor grain filling, or when the stalks fall, it causes significant losses and leads to failure in the harvesting of ears (Mesterhazy et al., 2012;Gai et al., 2018). However, the disease resistance mechanisms during maize stalk infection with F. graminearum are not clear, which limits progress in effective disease control. Several quantitative trait loci for resistance have been reported (Pe et al., 1993;Yang et al., 2010;Zhang et al., 2012;Chen et al., 2017;Ma et al., 2017). ZmAuxRP1, which encodes a plastid stroma-localized auxinregulated protein, is demonstrated as a key gene at qRfg2 through map-based cloning (Ye et al., 2018). ZmPRms has been speculated to be involved in resistance to infection by Aspergillus flavus and other pathogens (Majumdar et al., 2017b).
Proteins, as direct executors, participate in plant growth and development, secondary metabolism, and many other processes, reflecting the response of various physiological functions to biotic or abiotic stress stimuli (Hu et al., 2015;Yue et al., 2018;Zhang et al., 2019). Quantitative proteomics focuses on screening and identifying proteome variations among different species or states, revealing, and verifying changes in proteomics (Dunkley et al., 2006;Pechanova et al., 2013). Several studies have used proteomics to explore the response to biotic stress in maize (Huang et al., 2009;Pechanova et al., 2013). Comparative proteomics analysis of maize rachis protein explores resistance and susceptibility to aflatoxin accumulation (Pechanova et al., 2011). Proteome changes induced by F. graminearum inoculation for 48 h were evaluated in developing grains of two inbred lines (Mohammadi et al., 2011). To verify the hypothesis that the primary metabolism of Azospirillum brasilense Sp7 is different from that of C3 and C4 plants, differential proteomics was used to study the protein expression in maize and tomato with A. brasilense treatment (Lade et al., 2018). Proteomics data proved that effector HC toxin is a histone deacetylase inhibitor produced by fungal pathogen Cochliobolus carbonum race 1, which promotes maize virulence by altering protein acetylation (Walley et al., 2018). Comparative proteomics analysis was employed to select candidate proteins related to southern corn rust resistance in susceptible and resistant maize lines (Wang et al., 2019).
Transcription factors play a key role in crosstalk regulation between multiple hormones signaling pathways and signaltransduction-mediated defense gene expression. The WRKY transcription factors form a highly interconnected regulatory subnetwork in Arabidopsis immune response (Rushton et al., 2010;Llorca et al., 2014). WRKY33 is a key regulator in defense against Botrytis cinerea from the WRKY transcription factor family . At the early stages of B. cinerea infection, the expression level of WRKY33 is highly induced (Birkenbihl et al., 2017b;Sham et al., 2017). In the WRKY33 mutant, knock out of WRKY33 increases sensitivity of plants to B. cinerea, whereas overexpression of WRKY33 reduces susceptibility to B. cinerea (Zheng et al., 2006;Birkenbihl et al., 2012). After pathogen infection, WRKY33 can activate defense response mechanisms via MAP-dependent phosphorylation (Wang et al., 2018b;Han et al., 2019). In addition, WRKY33 homologous proteins can regulate plant disease resistance is confirmed in other plants, such as Solanum lycopersicum, Vitis vinifera, and Brassica napus (Wang et al., 2014;Merz et al., 2015;Zhou et al., 2015). Therefore, it is necessary to fully understand the plant defense mechanisms mediated by WRKY33 and its homologous proteins.
Maize stalk rot, caused by members of the F. graminearum species complex, is among the most destructive and economically important diseases in the world (Yang et al., 2010). However, very little is known about the defense responses against F. gramincarum infection in maize stem. In this study, we used high-throughput tandem mass tag (TMT)-based technology for proteomics comparison to explore pathogen-responsive proteins and biological processes in maize stems, which showed moderate resistance to F. gramincarum infection. We identified several defense-related proteins and revealed several secondary metabolism pathways associated with the defense responses. In addition, we also verified that ZmWRKY83 was a key transcription factor in plant disease resistance. This study not only provided several important information to understand the molecular mechanism of the interaction between maize and F. graminearum but also provided important clues for the genetic breeding of disease-resistant maize varieties.

Global Proteomic Changes With F. graminearum Infection in Maize
Fusarium graminearum can cause very serious stalk rot disease in maize. Therefore, it was important to explore the resistance mechanism of stalk rot caused by F. graminearum in maize. To investigate the global proteomic changes in response to F. graminearum infection in maize stem, comparative proteomics-based TMT labeling technology was conducted with three independent biological replicates for each sample ( Figure 1A). A total of 676,400 spectra were generated from TMT analysis, after searching the maize sequence database, yielding 85,386 matched peptides, 67,274 unique peptides, 8,471 matched proteins, and 8,044 quantified proteins ( Figure 1B and Supplementary Table 1), and the biological replications in different infection stage showed good similarity ( Figure 1C and Supplementary Figure 1). According to a criterion of 95% significance and a 1.3-fold cutoff (p-value < 0.05), 1,894 proteins were identified as differentially expressed proteins (DEPs) between 0, 1, and 2 days postinoculation (dpi; Supplementary Table 2). Then, we identified 1,064 DEPs in 1 dpi and 1,376 DEPs in 2 dpi that were significantly up or down expressed compared with 0 dpi (Figure 2). Among the 1,894 DEPs, we found 100 proteins were continuously upregulated ( Supplementary  Figure 2A), while 14 proteins were continuously downregulated after F. graminearum infection (Supplementary Figure 2B).

The Kyoto Encyclopedia of Genes and Genomes (KEGG) Enrichment Analysis of DEPs
To further analyze the key enriched metabolic pathways in response to F. graminearum infection, The KEGG pathway enrichment analysis was performed with pathogen-responsive proteins from the respective proteomes. The KEGG analysis demonstrated that phagosome, plant-pathogen interaction, and ribosome were significantly affected and upregulated by the early stage of F. graminearum infection (1 dpi).
Compared with 2 dpi proteomic data, flavonoid biosynthesis, phenylpropanoid biosynthesis, amino sugar, and nucleotide sugar metabolism, diterpenoid biosynthesis, glyoxylate, and dicarboxylate mentalism, biosynthesis of secondary metabolites, terpenoid backbone biosynthesis, and fatty acid degradation were significantly enriched in upregulated DEPs ( Figure 5). These results suggested that the disease resistance genes could mediate plant-pathogen interaction and the synthesis and metabolism of secondary metabolites in response to pathogen infection.
The KEGG pathway includes glutathione metabolism, starch and sucrose metabolism, oxidative phosphorylation, biotin metabolism, and photosynthesis that were enriched in downregulated DEPs (Figure 5), which suggested that energy metabolism and plant development were downregulated to resist F. graminearum infection. In addition, the results of KEGG analysis were relatively consistent with GO enrichment, and all of these results might reflect how plants resist fungal infection.

Transcriptional Expression Pattern of Some DEPs in Response to F. graminearum Infection
In KEGG pathway enrichment analysis, plant-pathogen interaction was significantly enriched in the upregulated DEPs with F. graminearum early infection, while terpenoid backbone and diterpenoid biosynthesis were enriched in later infection ( Figure 5). To further examine relationships between mRNA expression and protein abundance under F. graminearum infection, the relative transcript level of proteins enriched in plant-pathogen interaction, terpenoid backbone, and diterpenoid biosynthesis was examined by quantitative real-time PCR (qRT-PCR). Transcript levels of these proteins (excepted Zm00001d024903) enriched in plant-pathogen interaction were successfully determined and mostly increased in both 1 and 2 dpi, compared with 0 dpi (Figure 6). The proteins enriched in terpenoid backbone and diterpenoid biosynthesis were also successfully determined. Except for Zm00001d038193 and Zm00001d021709, most of them were observably upregulated during F. graminearum infection, especially in 2 dpi (Figure 7). In addition, we also examined the transcript levels of proteins enriched in the benzoxazinoid biosynthesis pathway which were downregulated, and all transcript levels of these proteins were decreased by F. graminearum infection.
Then, we compared the expression level of the genes enriched in plant-pathogen interaction, terpenoid backbone, and diterpenoid biosynthesis pathway in susceptible (B73) and resistant (Mo17) lines. We found that the expression levels of these genes were also significantly upregulated in Mo17 (Supplementary Table 4). However, most of the genes in Mo17 plants had higher multiples of expression than in B73, especially in 1 dpi. Moreover, we further extend the expression analysis in 0, 1, 2, and 3 dpi (Supplementary Table 5). Compared with 0 dpi, qRT-PCR results showed that most genes enriched in plant-pathogen interaction, terpenoid backbone, and diterpenoid biosynthesis pathways were significantly upregulated in 3 dpi, whereas increase multiple was not as high as 2 dpi. These The colors from yellow to red represent the significance level of the GO terms from high to low. Fisher's exact test, P-value was adjusted using the Benjamini-Yekutieli method, FDR < 0.05. results indicated that these genes played important roles in disease resistance.

ZmWRKY83 Was a Key Transcription Factor in Disease Resistance
In the plant-pathogen interaction pathway, we found that ZmWRKY83 (Zm00001d038023), the homolog of AtWRKY33 in Arabidopsis, was always enriched with high protein abundance. Meanwhile, ZmWRKY83 was continuously upregulated from 0 dpi to 2 dpi both in transcript and protein levels (Figure 6), which indicated that ZmWRKY83 was induced by F. graminearum infection and might play an important role in disease resistance. As a transcription factor, transcriptional activity is necessary. Therefore, we used the yeast two-hybrid system to verify the transcriptional activity of ZmWRKY83. The transformed yeast strain grew well on SD/-Leu/-Trp and SD/-Leu/-Trp/-His medium and appeared blue on the selective SD/-Leu/-Trp/-His/xα-gal medium, which indicated that the ZmWRKY83 transcription factor had the transcriptional activation activity (Supplementary Figure 3).
To investigate the function of ZmWRKY83 in disease resistance, ZmWRKY83 overexpressing transgenic plants were generated in Arabidopsis by Agrobacterium-mediated transformation ( Figure 8A). Then, we conducted a B. cinerea infection assay in Col-0, atwrky33 mutant, and ZmWRKY83 overexpressing transgenic plants. Compared with Col-0 and atwrky33 mutant, the resistibility of ZmWRKY83-OE plants to B. cinerea was relatively strong (Figure 8B and Supplementary Figure 4). In other words, the leaves of Col-0 and atwrky33 mutants were more sensitive to B. cinerea, whereas overexpression of ZmWRKY83 could enhance the resistance of Arabidopsis to this pathogen. Statistical analysis results showed that the size of the lesion area on Col-0 and atwrky33 mutant leaves after B. cinerea-infected was more than the lesion area on B. cinerea-infected ZmWRKY83-OE leaves ( Figure 8C). In addition, we got a Mutator insertion line of ZmWRKY83, in which ZmWRKY83 showed low FIGURE 4 | Gene ontology enrichment of down-regulated differentially expressed proteins between 0, 1, and 2 dpi. Enriched GO terms showed in biological progress (A), molecular function (B), and cellular component (C). The colors from yellow to red represent the significance level of the GO terms from low to high. Fisher's exact test, P-value was adjusted using the Benjamini-Yekutieli method, FDR < 0.05. expression ( Figure 8D). Then, we conducted F. graminearum infection assay in wild type B73 and zmwrky83 mutant. The results showed that the zmwrky83 mutant was more sensitive to F. graminearum infection (Figures 8E,F). These results suggested that ZmWRKY83 plays an important role in the plant defense system.

DISCUSSION
Plant diseases caused by pathogenic bacteria, such as fungi, bacteria, etc., have caused significant economic losses to important crops in the world (Ahuja et al., 2012). Fusarium species, such as F. graminearum, Fusarium verticillioides, and Fusarium proliferatum, are fungal pathogens associated with maize diseases (Dorn et al., 2011;Nguyen et al., 2016;Shu et al., 2017). Fusarium spp. can infect maize through all stages of growth and development and cause significant damage to different tissues of maize. For example, F. graminearum can cause ear rot and stem rot of maize and F. graminearum causes corn ear rot and stalk rot . Previous studies focused more on-ear and kernel rot in maize, as mechanized harvesting has attracted much more attention, stalk rot has become more important.
In a previous study, global proteomics using iTRAQ technology were employed to study the early infection by F. graminearum in maize kernels. A subset of 878 proteins was identified between mock and treatment. Many proteins related to defense response after infection are induced, including chitinases, pathogenesis-related (PR) proteins, and some differentially expressed proteins (DEPs) between maize-resistant and sensitive-inbred lines (Mohammadi et al., 2011). This study aimed to identify the defense response proteins of maize to F. graminearum infection. Proteomics analysis of maize infected by F. graminearum in different time points was performed using the TMT technique (Figure 1). As a result, a total of 1,894 DEPs were identified, which included 1,177 upregulated and 717 downregulated proteins (Figure 2). Function analysis of the DEPs during F. graminearum infection showed that might modulate the immune system by slowing down plant growth and energy metabolism (Figures 3-5).
Phenylpropanoids are secondary metabolites of plants derived from aromatic amino acid phenylalanine in most plants or tyrosine in some monocotyledons (Deng and Lu, 2017). Phenylpropanoids can participate in local and systemic signal transduction induced by defense response genes (Dixon et al., 2002). In the process of protecting plants from pathogen infection, the key enzymes of the phenylalanine pathway, such FIGURE 5 | Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis of upregulated and downregulated DEPs between 0, 1, and 2 dpi. The colors from blue-to-red represent the significance level of the KEGG pathways from low to high. as phenylalanine ammonia-lyase (PAL), peroxidase (POD), and polyphenol oxidase (PPO), are involved in the biosynthesis of lignin metabolites, which provides a structural barrier against fungi (Dixon et al., 2002;Mohammadi and Kazemi, 2002). In this study, we found that 29 upregulated proteins during F. graminearum infection were enriched in phenylpropanoid biosynthesis, including PODs, PALs, and other key enzymes (Supplementary Table 2).
In Arabidopsis, it is proved that calcium-dependent protein kinases (CDPKs) can mediate calcium-directed phosphorylation during plant defense activation and calmodulin-like (CML) proteins are important mediators of Ca 2+ -dependent signals during the immune response to plant pathogens (Harmon et al., 2000;Harper and Harmon, 2005;Leba et al., 2012). In plant-pathogen interaction pathway from upregulated DEPs enrichment analysis, two CMLs (Zm00001d029028, Zm00001d010649) and three CDPKs (Zm00001d021835, Zm00001d009708, and Zm00001d013030) were specifically upregulated during F. graminearum infection, which emphasizes the possible role of Ca 2+ -dependent signaling in the resistance to F. graminearum infection.
Production of PR proteins is a defense strategy often adopted by plants to resist the invasion of pathogens (Sels et al., 2008). PR protein has direct antibacterial properties, and some PR proteins can also inhibit pathogens by regulating key genes in the host defense system (van Loon et al., 2006). ZmPR10 plays a role in maize host resistance to A. flavus infection and aflatoxin contamination . The PR maize seed (PRms, Zm00001d009296) gene is proven to be involved in resistance to A. flavus and F. verticillioides infection (Shu et al., 2015;Majumdar et al., 2017a). In this study, ZmPR10, ZmPRms, and several PR proteins were upregulated both in transcript and protein levels during F. graminearum infection (Figure 6 and Supplementary Table 2), which suggested that ZmPRms played a key role in maize ear and stalk rot diseases and was a broadspectrum disease resistance gene.
Terpenoids are the most diverse in structure, and they act as plant antitoxins in direct defense of plants or signal in indirect defense reactions (Cheng et al., 2007). Terpenoids, monoterpenes, sesquiterpenes, semiterpenoids, and diterpenoids can participate in the protection of plants against abiotic and biotic stresses (Loreto et al., 2014). In this study, 13 proteins were detected and enriched in terpenoids and diterpenoids biosynthesis during F. graminearum's later infection (Figure 6 and Supplementary Table 2). These results indicated that the proteins involved in terpenoids and diterpenoids biosynthesis might be important to the resistance of fungal infection. Xylanase inhibitor proteins (XIP) are reported to be involved in the plant defense mechanisms against fungal pathogens and has been fully characterized in several model plant species (Misas-Villamil and van der Hoorn, 2008;Vasconcelos et al., 2011;Sun et al., 2018). Two XIPs were detected during F. graminearum infection and upregulated (Supplementary  Table 2), which suggested that these XIPs may act as a barrier to prevent fungal pathogens secreting xylanase from degrading cell walls. Several studies confirm that WRKY transcription factors are key regulators in response to abiotic stresses (Rushton et al., 2010). In maize, ZmWRKY17 is involved in regulating the expression of some ABA and stress response genes and negatively reduces ABA sensitivity, and regulates salt tolerance (Cai et al., 2017). ZmWRKY106 participates in a variety of abiotic stress response pathways and plays an active role under drought and high-temperature stress (Wang et al., 2018a). ZmWRKY33 can be induced by salt stress, drought stress, cold treatment, and ABA treatment, and overexpression in Arabidopsis plants can enhance salt tolerance (Li et al., 2013). ZmWRKY4 can be induced by Cd, and its role in regulating the expression and activity of SOD and APX in maize was studied by transient expression analysis based on maize protoplasts (Hong et al., 2017). However, the role of WRKY transcription factors in biotic stress is few studied in maize.
AtWRKY33 is involved in regulating gene transcription activity and it plays an important role in plant defense against B. cinerea (Birkenbihl et al., 2012(Birkenbihl et al., , 2017aSham et al., 2017). In this study, we found a WRKY transcription factor ZmWRKY83 (homolog of AtWRKY33) that may play a key role in maize resistance to F. graminearum infection by proteomics. A recent study shows that ZmWRKY52, ZmWRKY71, and ZmWRKY83 genes display significantly differential expression levels among the QTL-NILs, which indicated that these maize WRKY transcription factors may be involved in maize resistance to A. flavus (Parish et al., 2019). However, there is no direct experimental evidence that ZmWRKY83 was involved in resistance to pathogen infection. In this study, ZmWRKY83 protein was first identified to be associated with F. graminearum infection (Figure 6). Through overexpression of ZmWRKY83 in Arabidopsis and B. cinerea infection assay, it is shown that this gene was closely related to plant disease resistance (Figure 8), which indicated that ZmWRKY83 could be applied to molecular genetic breeding of maize disease-resistant varieties.
In summary, we used quantitative proteomics technology to compare the differential protein abundance after F. graminearum infection. Through a TMT-based proteomic approach, we identified 1,894 DEPs from maize stems challenged with the necrotrophic fungus pathogen F. graminearum. GO and KEGG pathway enrichment analysis showed that proteins related to plant-pathogen interaction and secondary metabolites biosynthesis were upregulated, whereas energy metabolism and plant development were downregulated to resist F. graminearum infection. In addition, we identified a disease-resistant transcription factor and proved that ZmWRKY83 plays an important role in plant immunity against fungal infection. Collectively, our results would also provide compelling evidence that the plant immune-related proteins are involved in the response to F. graminearum infection, which complements the latest knowledge about the maize resistance to F. graminearum infestation.

Plant Materials and Treatments
The maize inbred line B73 was grown in the greenhouse in Baoding (Hebei Agricultural University, Hebei province, China). Maize zmwrky83 mutant was obtained from the ChinaMu project (insertion site: chr 6, 145531076 1 ) (Liang et al., 2019). F. gramincarum PH-1 was grown on PDA plates at 28 • C for 6 days. Conidial suspensions were harvested by adding sterile distilled water containing 0.05% (v/v) Triton X-100 and scraping the plates using a glass spreader. The concentration of conidia was quantified using a hemocytometer and diluted to 1 × 10 6 spores/ml for inoculation. Maize stem was inoculated at the second or third internode using a sterile micropipet tip (10 mm hole depth) above the soil line, followed by injection of 20 µl macroconidia suspension. The plants were kept growing with the wounds being covered by sterile gauze to maintain moisture and avoid contamination from other organisms.

Protein Extraction, Quantification, and Digestion
The sample was ground into powder by liquid nitrogen and then transferred to a 5 ml centrifuge tube. Then, lysis buffer (including 1% TritonX-100, 10 mM dithiothreitol, and 1% protease inhibitor cocktail, 50 µM PR-619,3 µM TSA,50 mM NAM, and 2 mM EDTA) was added to the powder, followed by sonication three times on ice using a high-intensity ultrasonic processor (Scientz). After an equal volume of Tris-saturated phenol (pH 8.0) was added, the mixture was further vortexed for 5 min. After centrifugation (4 • C, 10 min, 5,000 g), the upper phenol phase was transferred to a new centrifuge tube. Proteins were precipitated by adding at least four volumes of ammonium sulfate-saturated methanol and incubated at −20 • C for at least 6 h. After centrifugation at 4 • C for 10 min, the supernatant was discarded. The remaining precipitate was washed with icecold methanol once, followed by ice-cold acetone three times. The protein was resuspended in 8 M urea (Sigma), and then the redissolved protein concentration was detected with a BCA kit (Beyotime) according to the instructions of the manufacturer.
The protein solution was reduced with 5 mM dithiothreitol (Sigma) for 30 min at 56 • C and alkylated with 11 mM iodoacetamide (Sigma) for 15 min in darkness at room temperature. The protein sample was then diluted by adding 100 mM TEAB to urea concentration less than 2 M. Trypsin (Promega) was added at 1:50 trypsin-to-protein mass ratio for the first digestion overnight and 1:100 trypsin-to-protein mass ratio for a second 4-h-digestion.

TMT Labeling and Peptide Fractionation
The peptide was desalted by Strata X C18 SPE column (Phenomenex) and vacuum-dried after trypsin digestion. According to the manufacturer's protocol, peptide was reconstituted in 0.5 M TEAB and processed with TMT kit (Thermo Fisher Scientific 2 ). Briefly, 1 unit of TMT reagent was thawed and reconstituted in acetonitrile. Then, the peptide mixture was incubated at room temperature for 2 h, and the pool was centrifuged in a vacuum, desalted, and dried. The tryptic peptides were fractionated into fractions by high pH reverse-phase HPLC using Agilent 300 Extend C18 column (5 µm particles, 4.6 mm ID, and 250 mm length). Then, the peptides were combined into 18 components with a gradient of 8-32% acetonitrile (pH 9.0) and dried by vacuum centrifugation.

LC-MS-MS Analysis
Trypsin peptides are dissolved in 0.1% formic acid and directly loaded on a self-made reversed-phase analytical column (15 cm length, 75 µm i.d.). The gradient of solvent B (0.1% formic acid in 98% acetonitrile) was gradually increased from 6% to 23% through 26 min, increased from 23% to 35% over 8 min, climbing to 80% in 3 min, and then maintained at 80% for 3 min, all at a constant flow rate of 400 nl/min on an EASY-nLC 1000 UPLC system. The separated peptides were subjected to NSI source followed by tandem mass spectrometry (MS/MS) in Q Exactive TM Plus (Thermo) coupled online to the UPLC.
A 2.0 kV electrospray voltage was applied. At 70,000 resolutions, the intact peptides were detected in the Orbitrap, with 350-1800 m/z full scan range. Up to 20 most abundant precursors were then selected for further MS/MS analyses with 30-s dynamic exclusion. The HCD fragmentation was performed at a normalized collision energy of 28%. The fragments were detected in the Orbitrap at a resolution of 17,500. The fixed first mass was set as 100 m/z. The automatic gain control target was set at 5E4, with an intensity threshold of 1E4 and a maximum injection time of 200 ms.

Database Search and Bioinformatics Analysis
Tandem mass spectrometry data were processed by the Maxquant search engine (version 1.5.2.8). MS/MS were searched against MaizeGDB (version 4, 131,585 sequences) database concatenated with reverse decoy database. Trypsin/P was designated as lyase, allowing up to two missing cleavages. The mass tolerance of precursor ions is set to 20 ppm in the first search, 5 ppm in the main search, and 0.02 Da in the fragment ions. The aminoformyl methyl on Cys was designated as a fixed modification, and the oxidation on Met was designated as a variable modification. The minimum score for peptides was set > 40 and FDR was adjusted to < 1%.
Functional annotations of DEPs were performed using agriGOv2 (Tian et al., 2017). The KEGG was used to predict the metabolic pathways and biochemical signals transduction pathways (Kanehisa and Goto, 2000). A p-value < 0.05 (Fisher's exact test) was used as the threshold to determine the significant enrichment of GO and KEGG pathways.

RNA Extraction and Quantitative Real-Time PCR
The same samples used in proteomics were homogenized in liquid nitrogen before RNA isolation. Total RNA of the samples was isolated using TRIZOL (Invitrogen) and purified using Qiagen RNeasy columns. qRT-PCR was conducted using actin as an internal reference and cDNAs from samples collected at different time points as the template with TransStart Tip Green qPCR SuperMix according to the manufacturer's instructions. Furthermore, comparative Ct analysis (2-Ct) of each gene in B73 and its relative expression levels at different time points were employed, and quantitative data were expressed as mean ± SEM. Primers used in this study were designed with Primer 6.0 software and listed in Supplementary Table 3.

Detection of ZmWRKY83 Transcription Activation Activity
To detect the transcriptional activation of ZmWRKY83, the CDS of ZmWRKY83 were amplified and cloned into the yeast vector pGBKT7 (Clontech) to obtain BD-ZmWRKY83 construct. BD-ZmWRKY83 construct and empty vector pGBKT7 were, respectively, transformed into yeast strain AH109. A transformed single colony (2 mm diameter) grown on selection medium for 3 days was resuspended in 100 µl of autoclaved distilled H 2 O, and 10 µl of resuspended cells was plated on SD/-Leu-Trp, SD/-Leu-Trp-His/x-α-gal for 3-5 days.

Construction of ZmWRKY83 Overexpression Plants in Arabidopsis
A full-length ZmWRKY83 CDS sequence was amplified by PCR with a specific primer (Supplementary Table 3). After sequence verification, the fragment was inserted into the binary vector pCAMBIA-1300. The recombinant plasmid was introduced into Columbia-0 (Col-0) plants by Agrobacterium using floral dipping. The transformed seeds were selected on a half-strength MS medium containing Basta.

B. cinerea Inoculation
Botrytis cinerea strain B05.10 was cultured on PDA (potato dextrose agar) media for 5-10 days at 18-22 • C. Conidia inoculums of B. cinerea were prepared by collecting conidia in sterile water, filtered to remove mycelia, and suspended to a final concentration of 10 6 conidia/ml. For resistance analysis, the single leaf of plants was inoculated with 5 µl conidia suspension of B. cinerea. Inoculated plants were kept in a growth chamber with high humidity for 24 h and then transferred to normal conditions. Lesion formation was monitored at 7 days postinoculation (dpi). Trypan blue staining was performed as described by the previous study (Hao et al., 2013).

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: The datasets generated and analyzed during the current study are available in the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD024342.

AUTHOR CONTRIBUTIONS
KZ, JX, and JD designed the experiments, analyzed the data, and wrote the study. HB, HS, and JZ analyzed the data, performed phenotype analysis, and qRT-PCR assays. XP, HC, and LY performed transcription activation activity assays. HB and HC constructed the ZmWRKY83 overexpression plants in Arabidopsis. KZ, HB, and XP identified the WRKY33 mutants and performed B. cinerea inoculation. All authors read and approved the contents of this study.