The Regulatory Network of CMPG1-V in Wheat–Blumeria graminis f. sp. tritici Interaction Revealed by Temporal Profiling Using RNA-Seq

Wheat powdery mildew (Pm), caused by Blumeria graminis f. sp. tritici (Bgt), is a prevalent fungal disease. The diploid wheat relative Haynaldia villosa (H. villosa) showed broad-spectrum resistance (BSR) to Pm. A previous study reported an E3 ligase gene, CMPG1-V from H. villosa, showing BSR to Pm. To elucidate the regulatory network mediated by CMPG1-V, in this study, gene expression profiling of CMPG1-V transgenic plant (CMPG1-VOE) and its receptor Yangmai 158 was analyzed and compared after Bgt inoculation at four infection stages. GO and KEGG analysis revealed obvious reprogramming of SA and ABA signaling, starch/sucrose metabolism, and photosynthesis in CMPG1-VOE, compared with those in Yangmai 158. Transcripts of SA synthesis genes SARD1 and UGT, signaling factors TGA and PRs, and SnRKs in ABA signaling were specifically upregulated in CMPG1-VOE rather than Yangmai 158. Transcripts of LHCII in photosynthesis, GLUC and TPP in starch/sucrose metabolism were also induced distinctly in CMPG1-VOE. WGCNA analysis showed crucial regulatory candidates of CMPG1-V, involving serine/threonine-protein kinase in phosphorylation, glucosyltransferase in flavonoid biosynthesis, defense factor WRKYs, and peroxidase in oxidative stress. Our results facilitate the deciphering of the resistant regulatory network of CMPG1-V and the identification of key candidates which might be employed in breeding programs.


Introduction
When pathogens infect plants, the early defense signals mainly involve the hormone signaling pathways and oxidative homeostasis [1]. Ubiquitination regulated by E3 ubiquitin ligase plays an important role during the process [2]. SCF COI1 type E3 ubiquitin ligase subunit F-box protein COI1 (Coronatine insensitive 1) is the first confirmed ubiquitination system against pathogen infection in JA signaling [3]. Arabidopsis (Arabidopsis thaliana) RING E3 ubiquitin ligase Keep On Going (KEG) interacting with SCF COI1 E3 ubiquitin ligase positively regulates plant immunity [4]. A tobacco (Nicotiana benthamiana) RING-type E3 ubiquitin ligase, NtRFP1, has been shown recently to attenuate disease symptoms caused by bC1 protein through the ubiquitination of bC1 [5].
Salicylic acid (SA), considered to be an important endogenous signaling molecule in plant immunity, can promote the degradation of defense transcription repressors through E3 ubiquitin ligase. SA can promote E3 ubiquitin ligase BTB NPR1 degrading transcriptional repressors TGA2 and NIMIN1 to enhance disease resistance [6]. CUL4-DDB1-based ubiquitin ligase interacted with SA to enhance tomato's resistance to nontumorigenic Agrobacterium tumefaciens [7]. AtSR1 interaction protein 1 (SR1IP1), a CUL3-based E3 ubiquitin ligase, positively regulates SA-related immunity by degrading the defense suppressor AtSR1 in Arabidopsis [8]. Several studies also confirmed the role of E3 ligases in ABA defense signaling. Overexpressing OsDRF1 in tobacco, an F-box type E3 ligase in rice, showed ABA sensitivity and increased resistance to viruses and bacteria [9]. Pepper (Capsicum annuum) RING E3 ubiquitin ligase CaRING1 induced by Xanthomonas oryzae acts as a positive regulator of defense responses [10,11].
Ubiquitin ligase also regulates plant primary metabolism after pathogen infection. Ubiquitin ligase ARABIDOPSIS TOXICOS EN LEVADURA31 (ATL31) and its closest homologue, ATL6, a membrane-associated ubiquitin ligase, are involved in the carbon/nitrogen (C/N) response by regulating the stability of 14-3-3 proteins through ubiquitination in plant immune response [12,13]. Overexpression of LeATL31 and LeATL6 enhanced the Pseudomonas syringae resistance in tomato [14].
Wheat (Triticum aestivum, AABBDD) makes a substantial contribution to human calorie intake, global agricultural sustainability, and food security [15]. Pathogens and pests cause wheat yield losses that account for approximately 21.5% of total losses, reaching up to 28.1% in food-deficit areas [16]. Wheat powdery mildew (Pm), a fungal disease caused by Blumeria graminis f. sp. tritici (Bgt), is one of the most destructive diseases of wheat, causing extensive yield losses worldwide [17]. Due to the host or race specificity of the wheat-Bgt interaction, new virulent races could rapidly evolve and emerge. Thus, the exploration and utilization of durable and broad-spectrum resistance (BSR) genes are of great significance for powdery mildew control.
Although there have many studies on the resistance mechanism of ubiquitin ligase genes, only a few studies performed the transcriptomic changes in response to Pm in a continuous time period. The related metabolism reprogramming in the disease resistance process is also indistinct. In a previous study, an E3 ligase gene, CMPG1-V, was cloned from Haynaldia villosa. Overexpression of CMPG1-V enhanced its BSR to Pm [18]. This provides an ideal system for investigating the temporal process during Bgt infection regulated by CMPG1-V. In this study, CMPG1-V transgenic line (CMPG1-V OE ) and susceptible receptor Yangmai 158 were inoculated with Bgt isolate E26 for RNA-seq. CMPG1-V associated defense pathways were analyzed by the Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment, and Weighted Gene Co-Expression Network (WGCNA). These findings provide insights into the molecular mechanisms of CMPG1-V-associated broad-spectrum powdery mildew resistance in wheat.

Global Profiling of CMPG1-V OE in Response to Bgt Infection
When Bgt infects wheat, conidia forms rapidly and reaches the epidermal cells from 1 h after inoculation (hai) to 72 hai ( Figure S1) [19]. Based on the above observation, the samples for RNA-seq of CMPG1-V transgenic line (CMPG1-V OE ) and susceptible receptor Yangmai 158 were collected at 1, 8, 18, and 24 hai by moderately Bgt virulent race E26. The mock samples were collected simultaneously at the above corresponding time points. The Q20 and Q30 of RNA-seq sequences for all samples were more than 93% and 84% and the ratios of clean reads accounted for 89.41% to 95.17%, respectively (Table S1). To maximize the number of possible differentially expressed genes (DEGs), the selective parameter was set as Fold Change ≥ 2.00 and FDR ≤ 0.001. In CMPG1-V OE and Yangmai 158, 6363 and 8047 transcripts were upregulated, and 11,475 and 5114 were downregulated, respectively (Figure 1a). The number of upregulated DEGs in CMPG1-V OE was less than that in Yangmai 158 for all the analyzed infection time points. Venn diagram revealed 24% of the upregulated DEGs were specifically expressed in CMPG1-V OE , while for the downregulated DEGs, up to 60.8% were specifically expressed in CMPG1-V OE , and only 12% were specifically expressed in Yangmai 158 (Figure 1b). We have also tested the two materials with a more virulent race, E31, and harvested the samples for RNA-seq at the same time points after inoculation (Table S1 and Figure S2).

Functional Categories Enriched in CMPG1-V OE and Yangmai 158 in Response to Bgt Infection
To understand the biological function of the DEGs, GO enrichment was performed to classify the crucial biological processes and cardinal signaling pathways. Before Bgt infection, no significant difference was observed between CMPG1-V OE and Yangmai 158 ( Figure S3a). The upregulated DEGs in CMPG1-V OE were enriched in normal physiological metabolism, such as "sulfur compound transport", "aminoglycan metabolic process", "glutathione metabolic process", "oxoacid metabolic process", and so on.
After Bgt inoculation, we observed distinct enriched DEGs in the two materials (Figure 2a,b). At 1 hai, most of the upregulated DEGs in the two materials focused on "photosynthesis" and "generation of precursor metabolites and energy", and transcription levels increased in the succeeding infection timepoints. At later infection stages from 18 hai to 24 hai, the upregulated transcripts in "abscisic acid (ABA)-activated signaling pathway", "regulation of salicylic acid (SA) metabolic process", "glucose import and sucrose biosynthetic process", "nitrogen compound", and "nitric oxide metabolic process" were specifically enriched in CMPG1-V OE (Figure 2b). However, at the corresponding infectious timepoints in Yangmai 158, the upregulated transcripts were only enriched in "organonitrogen compound catabolic process", "response to cytokinin", and "carbohydrate metabolic process" (Figure 2a). Similarly, the hormone and energy biological processes were also observed in CMPG1-V OE challenged with E31 ( Figure S4).
Common up-or downregulated DEGs were also present in the two materials after Bgt inoculation. The upregulated transcripts of common defense response in CMPG1-V OE and Yangmai 158 were enriched in "cell wall organization or biogenesis", "glutathione metabolic process", "toxin catabolic process", and "sulfur compound transport" ( Figure S3b). The downregulated transcripts were enriched in "carboxylic acid biosynthetic process", "fatty acid biosynthetic process", "response to acid chemical", and "chitin metabolic process" ( Figure S3c).

Temporal Specificity of CMPG1-V OE in Response to Bgt Infection
To characterize the pathways regulated by CMPG1-V at different time points, temporal specificity of CMPG1-V OE was analyzed (Figure 3a). At 1 hai, "ABA and SA pathways" firstly appeared in CMPG1-V OE . SA signaling pathway was only present in CMPG1-V OE at 1 hai, and ABA-associated processes were sustained during all the tested time points in CMPG1-V OE . "Jasmonic acid (JA) pathway" was induced in CMPG1-V OE at 8 hai, while it was induced at 1 hai and sustained from 18 to 24 hai in Yangmai 158 after Bgt infection. Moreover, transcripts in "photosynthesis" and "generation of precursor metabolites and energy" were significantly upregulated in CMPG1-V OE at 1 hai and thereafter. From 18 hai to 24 hai, the "nitrogen compound transport" was specifically enriched in CMPG1-V OE . The "sucrose biosynthetic process" and "glucose import" were enriched specifically in CMPG1-V OE at 18 hai and 24 hai, respectively. In Yangmai 158 at 18 and 24 hai, there was no significant enrichment of the above three pathways, but "photosynthesis" and "generation of precursor metabolites and energy" were specifically enriched. The "carbohydrate metabolic process" and "organonitrogen compound catabolic process" were activated in Yangmai 158, while not in CMPG1-V OE from 1 hai to 24 hai.
KEGG enrichment analysis was used to classify the crucial biological processes and cardinal signaling pathways (Figure 3b). At 1 hai, the "plant-pathogen interaction" and "plant hormone signal transduction" were obviously upregulated in CMPG1-V OE , and the transcription level in "regulation of autophagy" increased twofold more in CMPG1-V OE than in Yangmai 158. At 8 hai and the later infection stages, the "phenylalanine metabolism" was activated and increased up to threefold more in CMPG1-V OE than in Yangmai 158. Metabolic processes such as "fatty acid elongation" increased more than twofold in CMPG1-V OE . With the infection going on, "photosynthesis" was activated both in CMPG1-V OE and Yangmai 158 at 24 hai, but more significantly in CMPG1-V OE . The downregulated DEGs were mainly enriched in "glycine, serine, and threonine metabolism", "ABC transporters", "isoflavonoid biosynthesis", and "glutathione metabolism" (Figures S5 and S6). Y represents Yangmai 158, C represents CMPG1-V OE . The colors white, pink, and red represent low, medium, and high expression levels, respectively. Upregulated pathways were selected as follows: P value < 0.05; genes with the regulation ratio log ≥ 2 or ≤−2 were selected.

Phytohormone Signaling Was Rapidly Reprogrammed in CMPG1-V OE after Bgt Infection
Transcription analysis indicated obvious biosynthesis and signaling changes in ABA and SA pathways in CMPG1-V OE in response to E26 infection. For the ABA pathway, its degradation gene CYP707A2 decreased more than threefold, while its synthesis gene ABA3 increased more than 4.5-fold at 18 hai in CMPG1-V OE (Figure 4a). More, PP2C, a negative regulatory gene in ABA signaling, was downregulated in CMPG1-V OE for all the analyzed infection processes, while transcripts of SnRK2 and SnRK3 were activated evidently in CMPG1-V OE in all infectious time points (Figure 4a). SnRK2.1 increased more than twofold and SnRK2.5 increased more than fivefold in CMPG1-V OE than in Yangmai 158. Besides, CIPK23 (SnRK3) increased more than sevenfold. The qRT-PCR also validated the RNA-seq data. As shown in Figure 5, the expression level of synthesis gene ABA3 peaked at 18 hai in CMPG1-V OE , not in Yangmai 158; PP2C decreased in CMPG1-V OE from 18 hai to 24 hai, while SnRK2.1 increased to the maximum at 18 hai and SnRK2.5 peaked at 8 hai in CMPG1-V OE ; and the expression level of CIPK23 maximized at 8 hai in CMPG1-V OE , not in Yangmai 158.
For the SA pathway, SA synthesis and transduction pathway were distinctly more upregulated in CMPG1-V OE . Transcript levels of critical genes for SA synthesis, SAR DEFICIENT1 (SARD1) and UDP-glycosyltransferase (UGT), were obviously increased in CMPG1-V OE when infected by Bgt. The UGT was threefold higher in CMPG1-V OE from 1 hai to 18 hai and SARD1 was threefold more in CMPG1-V OE than in Yangmai 158 from 1 hai to 18 hai (Figure 4b). The activation of the SA signal pathway was shown by the TGA1 expression increasing more than sixfold in CMPG1-V OE from 1 hai to 18 hai; the PR was more than fourfold in CMPG1-V OE at 8 hai and increased to ninefold at 24 hai (Figure 4b), while no obvious induction of these genes was observed in Yangmai 158 ( Figure 4b). qRT-PCR of TaTGA1 and TaPR1 validated the results from the RNA-seq data ( Figure 5). The expression level of TaTGA1 peaked at 24 hai in CMPG1-V OE and TaPR1 expression rose in CMPG1-V OE from 18 hai to 24 hai, which were more than those in Yangmai 158.  In addition, transcription levels of COI1 in JA signaling and ETR1 in ET signaling were upregulated in CMPG1-V OE (Figure 4c,d). The expression levels of COI1 and ETR1 maximized at 18 hai in CMPG1-V OE , which was more than those in Yangmai 158 ( Figure 5). The upregulated hormone signaling and associated DEGs were also observed in CMPG1-V OE or Yangmai 158 tested with E31 ( Figure S7).

CMPG1-V Activates Conspicuous Energy Metabolic Signaling during Bgt Infection
Nitrogen assimilation genes nitrate reductase (NR) and glutamine synthetase (GS) were induced apparently in CMPG1-V OE (Figure 6a). The transcript levels of NR and GS increased more than fivefold and threefold in CMPG1-V OE , respectively. Expression levels of NR and GS increased to the maximum at 24 hai in CMPG1-V OE , which were more than those in Yangmai 158 (Figure 7). Glutamate dehydrogenase 1 (GDH1) transcripts rose up to more than sixfold in CMPG1-V OE compared with Yangmai 158. The starch and sucrose metabolism was enriched from 18 hai to 24 hai, and the transcripts of galacturonosyltransferase (GAUT) increased fivefold in CMPG1-V OE . This consisted of the pectinesterase, UDP-glucuronate 4-epimerase (GAE), UDP-glucose 6-dehydrogenase (UGDH), beta-glucosidase (GLUC), and trehalose 6-phosphatephosphatase (TPP), whose transcripts increased fourfold and reached their peaks at 24 hai ( Figure 6b). As shown in Figure 7, expression levels of GLUC and TPP rose to the top at 24 hai in CMPG1-V OE , which were more than those in Yangmai 158 ( Figure 7). CMPG1-V activated the glycolysis/gluconeogenesis pathway. With the infection going on, the transcription levels of ATP-dependent 6-phosphofructokinase (pfkA), pyruvate kinase (PK), pyruvate dehydrogenase E1 alpha subunit (PDHA), and aldehyde dehydrogenase (ALDH) increased prominently in CMPG1-V OE (Figure 6c). The transcripts of PfkA increased fourfold in CMPG1-V OE , and those of PK and PDHA showed a similar pattern. The aldehyde dehydrogenase transcripts were fourfold more from 18 hai until 24 hai. The expression level of pfkA peaked at 24 hai in CMPG1-V OE , while ALDH maximized in CMPG1-V OE at 18 hai, not in Yangmai 158 (Figure 7). In addition, transcription of NAD-dependent malic enzyme (NAD-ME) and tetraspanin-19 (TSPAN) associated with photosynthesis were activated. Their expression levels were more than sixfold in CMPG1-V OE compared with those in Yangmai 158 (Figure 6d). The chlorophyll a-b binding protein (LHCII) increased sevenfold in CMPG1-V OE . The expression level of TSPAN maximized at 18 hai in CMPG1-V OE , LHCII increased to a maximum in CMPG1-V OE at 24 hai, not in Yangmai 158. The NAD-ME had a significant increase in CMPG1-V OE from 1 hai to 8 hai, which was more than that in Yangmai 158 (Figure 7). Similar results were observed in CMPG1-V OE or Yangmai 158 tested with E31 for the induction or inactivation of nitrogen and glycolysis/gluconeogenesis metabolism ( Figure S8).

Weighted Gene Co-Expression Network (WGCNA) of Pm Resistance Regulated by CMPG1-V
To identify crucial regulatory factors or pathways in the CMPG1-V defense network, the co-expression data were analyzed from all DEGs in CMPG1-V OE and Yangmai 158 through WGCNA. The network contained 16,662 DEGs after filtering by the coefficient of variation. As shown in Figure S9, 28 module clusterings were set up and the samples collected at the four infectious time points were corresponded to each module ( Figure S9a,b).
In total, 2259 DEGs enriched in all modules and CMPG1-V were classified into the pink module ( Figure S9c). KEGG enrichment analysis revealed DEGs in the pink module were closely related to plant-pathogen interaction, nitrogen metabolism, phenylalanine metabolism, and ubiquinone and terpenoid-quinone biosynthesis (Figure 8a). Twenty-six DEGs shown as rose red nodes could be categorized into plant-pathogen interaction, including NCED5 and PP2C35 in the ABA pathway, EDS1B in the SA pathway, and RPM1 and PTI1-like tyrosine-protein kinase 1 in disease resistance (Figure 8b). Twenty-four DEGs shown as light blue nodes were associated with phytohormone metabolism, including beta-glucosidase (GLUC), 6-phosphofructokinase 3 (PFK3), and trehalose 6-phosphatephosphatase (TPP) of starch and sucrose metabolism, 4-coumarate-CoA ligase (4CL) of lignin metabolism, and anthocyanidin 5,3-O-glucosyltransferase (GT) of anthocyanin biosynthetic pathway (Figure 8b). In addition, 20 yellow nodes were associated with protein modification and 14 blue nodes involved in oxidation-reduction. Definitively, CMPG1-V interacting proteins previously identified by yeast two-hybrid assay were also identified, including 4-coumarate-CoA ligase-like, beta-glucosidase, protein DETOXIFICATION 19-like, and heavy metal-associated isoprenylated plant protein (Table S2). The DEGs including beta-glucosidase, trehalose 6-phosphatephosphatase, and NAD-ME were confirmed by qRT-PCR (Figure 7). Taken together, Table 1 lists the probable DEGs participating in the CMPG1-V defense network including protein phosphorylation modification, flavonoid metabolism, ABA/SA signaling, and oxidation homeostasis. Rose red nodes are related to signaling, light blue nodes are related to metabolism, blue nodes are related to oxidation-reduction, yellow nodes are related to protein modification. Red label font represents genes interacting with CMPG1-V obtained by yeast two-hybrid. Green label font represents genes having qRT-PCR results in previous study. Genes with the edge weight higher than 0.11 are visualized by Cytoscape.

The Chromosomal Distribution of Candidate Genes
The chromosome distributions of the DEGs related to plant hormone pathways, energy metabolic signaling pathways, and in the network of the CMPG1-V module were surveyed through in silico mapping using Chinese Spring genomic sequences (http://www.wheatgenome.org/) as reference ( Figure S10 and Table S3). In addition, all of them have homologous genes in Arabidopsis, with a similarity reached up to 80% (Table S3). As is shown in Figure S10, the 22 DEGs in plant hormone pathways were assigned to 11 different chromosomes; 29 DEGs in energy metabolism were assigned to 14 different chromosomes; and the remaining 23 DEGs in the network of CMPG1-V were assigned to 16 different chromosomes. Chromosome distribution analysis of the DEGs might facilitate the association study between these DEGs and previously reported genes in resistance to powdery mildew.

Discussion
Wheat powdery mildew is one of the most destructive diseases of wheat, causing extensive yield losses worldwide [17]. Due to the host or race specificity of the wheat-Bgt interaction, a new virulent race could rapidly evolve and emerge. Here, CMPG1-V, an E3 ligase gene from H. villosa, showed BSR to powdery mildew when overexpressed in wheat [18], thus the transgenic line CMPG1-V OE and its receptor variety Yangmai 158 provide an ideal system for studying BSR during wheat-Bgt interaction. We identified 3542 DEGs, which were only specifically upregulated in CMPG1-V OE after Bgt infection.

Phytohormone SA and ABA Play Important Roles in CMPG1-V-Associated Powdery Mildew Resistance
Plants are continuously exposed to diverse phytopathogenic microorganisms and have elaborated a variety of defense mechanisms to successfully avoid infection by limiting pathogen invasion and multiplication [35,36]. Plant hormone SA plays an important role in the interaction between plants and pathogens [37]. The exogenous application of SA conferred resistance against Magnaporthe oryzae [38]. SA-deficient NahG tobacco was hypersusceptible to powdery mildew (Euoidium longipes), as judged by significantly more severe powdery mildew symptoms and enhanced pathogen accumulation [39]. In a previous study, exogenous SA could increase the expression level of CMPG1-V, and TaPR1 and TaPR2 were upregulated significantly in CMPG1-V OE when infected by Bgt [18]. RNA-seq showed obvious transcription level increase of the SA synthesis genes SARD1 and UGT in CMPG1-V OE , hinting at the important role of the SA pathway in Pm resistance mediated by CMPG1-V. Arabidopsis AtSARD1 could enhance the resistance against Verticillium dahliae [40]. Overexpressing TaUGT3 enhanced fusarium head blight resistance of wheat [41]. What's more, TGA TFs family positively regulates defense responses against biotrophic and necrotrophic pathogens [42], consistent with the transcription level of TGA accelerated in CMPG1-V OE . Thus, it will be important to elucidate the synthesis of SA during CMPG1-V regulating Bgt resistance. And what is more, the strategies targeting TAG TFs are more crucial for enhancing powdery mildew resistance during this process.
ABA, originally described for their function in response to abiotic stresses, has emerged as a crucial player in plant-pathogen interactions [43,44]. Exogenous ABA quickly induced CMPG1-V, suggesting that CMPG1-V may be involved in ABA-dependent defense response [18]. Here, we found transcripts of SnRK2 were upregulated in all infectious periods. In Arabidopsis, SnRK2.8 phosphorylated NPR1, enhancing systemic acquired resistance (SAR) [45]. OsSAPKs (SnRK2) play positive roles in response to pathogen infection in rice [46]. ABA-inducible SnRK2-type kinase OsSAPK10 phosphorylated WRKY72, released its suppression on AOS1, and enhanced resistance to Xanthomonas oryzae pv. oryzae in rice [47]. Besides, PP2C, a negative regulator in ABA signaling, was significantly downregulated in CMPG1-V OE after Bgt inoculation compared with Yangmai 158. Some rice PP2Cs subgroup K negatively regulate X. oryzae pv. oryzae infection [48]. Thus, PP2Cs are thought to play a specific role in CMPG1-V-associated defense pathways, and ABA-dependent defense signaling possibly plays a positively regulatory role in CMPG1-V resistance network.
Chromosome distribution analysis of the DEGs provided information to identify whether these DEGs were associated with previously identified genes or QTL. By gene expression comparison between resistant and susceptible durum wheat lines, a WAK2 gene in the mapped region of QFhb.mgb-2A was identified and responsible for Fusarium Head Blight resistance [49,50]. The DEGs which might be involved in CMPG1-V regulatory network have been mapped in silico onto chromosome regions of wheat (Table S3). We have found that a CIPK29 gene (Unigene49646_All) which was located in the terminal region 2BS might be associated with powdery-mildew-resistant gene Ml5323 in the interval of 2BS FL0.84-1.00 [49]. CIPKs were reported to play roles in the ABA signaling pathway [51]. This finding might provide clues to the cloning and elucidation of the resistant mechanism mediated by Ml5323.

CMPG1-V Reprogrammed Starch and Sucrose Metabolism/Photosynthesis in Response to Bgt Infection
When a plant defends against pathogen invasion, levels of resistance in the whole plant are influenced by systemic signals mediated by plant hormones [52]. Plant hormones can coordinate various signal pathways to mediate metabolism under pathogen infection [1]. NR-mediated NO generation plays a key role in protecting plants from abiotic stresses through activating antioxidant enzymes [53]. GS1 can be induced during development of the infectious process [54]. We observed that, in response to Bgt infection, the nitrogen assimilation genes NR and GS were upregulated apparently in CMPG1-V OE . What is more, typical stress metabolite trehalose plays an important role in regulation of defense response against pathogens. Some selected SlTPPs could be inducted by Botrytis cinerea and Pseudomonas syringae pv. tomato (Pst) DC3000 as well as defense signaling hormones [55]. We found TPP involving starch and sucrose metabolism and pfkA, PK, PDHA involved in glycolysis/gluconeogenesis were induced notably in CMPG1-V OE . The GAPDHs, multifunctional enzymes in glycolysis, have been proven to be involved in regulation of ROS, autophagy, and plant immune responses [56]. Photosynthesis also modulates plant defense responses induced by pathogen infection and abiotic signals such as light, circadian rhythm, and temperature [57]. During infection with Colletotrichum higginsianum, loss-of-function mutants of NADP-ME2 (nadp-me2) showed enhanced susceptibility [58]. The transcription levels of NAD-dependent malic enzyme, chlorophyll a-b binding protein, and tetraspanin-19 were upregulated remarkably in CMPG1-V OE , compared to those in Yangmai 158. Thus, the reprogramming of starch/sucrose metabolism and photosynthesis pathways played an important role in the CMPG1-V defense response to Bgt infection.

CMPG1-V Exerts Its Distinct Defense Response by Reprogramming a Specific Network
Plants have evolved innate immune systems that recognize the presence of potential pathogens and initiate effective defense responses. E3 ubiquitin ligase plays a crucial role in the plant immunity process [2]. There are, however, a few notable highlights or differences in resistance mechanisms against Bgt in CMPG1-V. WGCNA of CMPG1-V module revealed its close relation with phosphorylation modification, flavonoid metabolism, oxidation homeostasis, and WRKY. Crosstalk between ubiquitination and phosphorylation modification plays a key role in controlling the defense signaling events to ensure plant pathogen resistance [2]. Crosstalk between posttranslational modifications, such as ubiquitination and phosphorylation, play key roles in controlling the duration and intensity of signaling events to ensure cellular homeostasis [59]. Phosphorylation events contribute additively to the stabilization of PUB22 in response to the perception of pathogen-associated molecular patterns (PAMPs) [59]. Phosphorylation of NPR1 could switch autoubiquitination activity of BTB NPR1 to substrate ubiquitination by enhancing its affinity to targets such as TGA2 [60]. For oxidation homeostasis, E3 ubiquitin ligase can induce ROS to enhance the resistance level. Knockdown of stress inducible OsSRFP1, encoding an E3 ubiquitin ligase with transcriptional activation activity, conferred abiotic stress tolerance through enhancing antioxidant protection in rice [61]. Arabidopsis E3 ubiquitin ligase PUB13 regulated chitin receptor LYSIN MOTIF RECEPTOR KINASE5 and induced rapid responses, such as the production of ROS [62]. Moreover, other E3 ligases can interact with WRKY to enhance the resistance. The E3 ligase likely to be involved in plant immunity is UPL5, which was first identified by a yeast two-hybrid screen. UPL5 interacted with WRKY53, a transcription factor acting positively in leaf senescence [63]. The Chinese wild grapevine (Vitis pseudoreticulata) E3 ligase EIRP1 activated plant defense responses by inducing proteolysis of the VpWRKY11 transcription factor [64]. However, there are few studies about ubiquitination in plant immunity related with flavonoid metabolism. In this study, we found flavonoid metabolism may have an influence on ubiquitination in plant immunity. This will give us a new perspective to clarify the resistance function of CMPG1-V in wheat powdery mildew resistance. On all accounts, these candidates will provide extensive insights into molecular mechanisms of CMPG1-V-associated broad-spectrum powdery mildew resistance in wheat.

Plant Material and Fungal Isolates
Yangmai 158 is a moderate powdery-mildew-susceptible wheat variety and it is susceptible to Bgt isolates E26 or E31. Overexpression of CMPG1-V in Yangmai 158 enhanced its broad-spectrum resistance. Bgt isolates E26 and E31 were collected from Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing, China. They are maintained on seedlings of susceptible wheat variety Sumai 3 in a spore-proof greenhouse. All seedlings were grown in a growth chamber with 20 • C/16 • C (day/night), 16 h/8 h (light/dark). CMPG1-V transgenic line (CMPG1-V OE ) and the receptor variety Yangmai 158 were inoculated with Bgt isolates E26 and E31 at the two-leaf stage and RNAs were isolated at 1 h, 8 h, 18 h, and 24 h before and after Bgt inoculation, followed by freezing in liquid nitrogen for subsequent RNA extraction.

RNA-Seq Library Construction and Sequencing
Total RNA extraction was performed using the Trizol reagent (Invitrogen, Waltham, MA, USA), according to the manufacturer's instructions. RNA concentration was measured using a NanoDrop spectrophotometer and 1.2% agarose gel electrophoresis.
Sample detection and sequencing were performed by BGI-Shenzhen, Shenzhen Beijing. After extracting total RNA and treating with DNase I, Oligo(dT) was used to isolate mRNA. Mixed with the fragmentation buffer, the mRNA were fragmented. Then, cDNA was synthesized using the mRNA fragments as templates. Short fragments were purified and resolved with EB buffer for end reparation and single nucleotide A (adenine) addition. After that, the short fragments were connected with adapters. The suitable fragments were selected for the PCR amplification. During the QC steps, Agilent 2100 Bioanaylzer and ABI Step One Plus Real-Time PCR System were used in quantification and qualification of the sample library. Then, the library was sequenced using Illumina HiSeq 4000 (Illumina, San Diego, CA, USA). The raw data and processed data could be obtained from the NCBI website address. The initial short reads data sets were available at the NCBI Short Read Archive (SRA) with the accession number SRR12416499-SRR12416522.

Sequencing Reads Filtering and De Novo Assembly
To obtain clean reads, raw reads, including reads with adaptors, reads in which unknown bases represented more than 5% of the total bases, and low-quality reads (percentage of low-quality bases with a quality value ≤ 10 in more than 20% of a read), were removed. We used Trinity to perform de novo assembly with clean reads that PCR duplication removed (in order to improve the efficiency), then used Tgicl to cluster transcripts to Unigenes. Trinity

Detection and Analysis of Differentially Expressed Gene
We detected DEGs with PossionDis as requested and parameters were set as Fold Change ≥ 2.00 and FDR ≤ 0.001. With the GO annotation result, we classified DEGs according to official classification, and we also performed GO functional enrichment using phyper, a function of R (version: v3.4.2). With the KEGG annotation result, we classified DEGs according to official classification, and we also performed pathway functional enrichment using phyper, a function of R. We calculated false discovery rate (FDR) for each p-value; in general, the terms with FDR not larger than 0.001 were defined as significantly enriched. Hierarchical clustering and heat map generation were performed in R. The gene expression data were log2-transformed and then quantile-normalized prior to generating the heat map for direct comparison of the data.

Weighted Gene Co-Expression Network Analysis and Visualization
Co-expression networks were constructed using the WGCNA package in R. We used the log2-transformed FPKM values, and replaced values smaller than one by zero. In total, 16,662 genes were used for the WGCNA analysis. The modules were obtained using the automatic network construction function block wise Modules with default settings. Parameters were set as follows: the power was 20, TOM-Type was unsigned, min Module Size was 30, and merge Cut Height was 0.25. We inferred an undirected, weighted network. Within the network, clusters of genes with similar expression patterns, called modules, were then inferred using a clustering of the Topological Overlay Matrix. The eigengene value was calculated for each module and 16,662 genes were clustered into 27 specific modules. The networks were visualized using Cytoscape _v.3.5.1.

Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) Assay
All selected genes were used for SYBR green real-time RT-PCR; primers are listed in Table S4. The Tubulin gene was used as the internal control gene. The qRT-PCR reaction for target gene transcript amplification was carried out in a final volume of 20 µL containing 20 ng cDNA, 0.2 µM of each primer, 2× AceQ qPCR SYBR Green Master Mix (Vazyme, Nanjing, China). The PCR reaction conditions were denaturation at 95 • C for 5 min followed by 40 cycles of 95 • C for 10 s, annealing at the appropriate temperature (from 57 to 61 • C) for 30 s, extension at 72 • C for 30 s, followed by 95 • C for 15 s, 60 • C for 1 min, then 95 • C for 15 s to obtain melt curves to ensure primer specificity. All reactions were done in triplicate. Reactions were performed on a Bio-Rad IQ single-color Real-Time PCR detection System (Bio-Rad, Hercules, CA, USA). The data were compiled from the mean Ct values of all the replicates after normalizing with the Ct values of the endogenous control. The relative expression level in terms of fold change was calculated using the 2 −∆∆CT method [69]. All the qRT-PCRs had three biological replicates.

Statistical Analysis
All statistical analyses were performed using SPSS version 13.0, and one-way ANOVA was performed with a homogeneity of variance test, followed by an LSD test to check for quantitative differences between treatments. p < 0.05 was set as the significance cut-off.

Conclusions
In this study, we analyzed the temporal specificity of Pm resistance regulated by CMPG1-V at the whole genome level, represented by four time points after Bgt infection. We suspected that when the pathogen attacked, the plant activated CMPG1-V, then enhanced hormone defense signaling through SA and ABA signaling. Activation of C/N and photosynthesis metabolism provided energy during wheat-Bgt interaction, which affected the expression of downstream PR genes ( Figure 9). Moreover, phosphorylation modification probably regulated different transcription of genes in defense signaling and related metabolisms. Our findings point to strengthening the utility of CMPG1-V in enhancing BSR to powdery mildew.

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