{'fraction_non_alphanumeric': 0.047077162899454406, 'fraction_numerical': 0.04829306313328137, 'mean_word_length': 4.80246020260492, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Proteomic Analysis of UV-B-Induced Virulence-Mutant Strains of Puccinia striiformis f. sp. tritici Based on iTRAQ Technology

The emergence of new physiological races of Puccinia striiformis f. sp. tritici (Pst) causing wheat stripe rust can lead to the loss of resistance of wheat cultivars to stripe rust, thus resulting in severe losses in wheat yield. In this study, after the germination of urediospores of three Pst strains including the original strain (CYR32, a dominant physiological race of Pst in China) and two virulence-mutant strains (CYR32-5 and CYR32-61) acquired from CYR32 via UV-B radiation, proteomic analysis based on isobaric tags for relative and absolute quantification (iTRAQ) technology was performed on the strains. A total of 2,271 proteins were identified, and 59, 74, and 64 differentially expressed proteins (DEPs) were acquired in CYR32-5 vs. CYR32, CYR32-61 vs. CYR32, and CYR32-61 vs. CYR32-5, respectively. The acquired DEPs were mainly involved in energy metabolism, carbon metabolism, and cellular substance synthesis. Furthermore, quantitative reverse transcription PCR assays were used to determine the relative expression of the 6, 7, and 1 DEPs of CYR32-5 vs. CYR32, CYR32-61 vs. CYR32, and CYR32-61 vs. CYR32-5, respectively, at the transcriptional level. The relative expression levels of one, five, and one gene, respectively, encoding the DEPs, were consistent with the corresponding protein abundance determined by iTRAQ technology. Compared with CYR32, the DEPs associated with energy metabolism and stress—including E3JWK6, F4S0Z3, and A8N2Q4—were up-regulated in the mutant strains. The results indicated that the virulence-mutant strains CYR32-5 and CYR32-61 had more tolerance to stress than the original strain CYR32. The results obtained in this study are of great significance for exploring the virulence variation mechanisms of Pst, monitoring the changes in Pst populations, breeding new disease-resistant wheat cultivars, and managing wheat stripe rust sustainably.


INTRODUCTION
Stripe rust (yellow rust), caused by Puccinia striiformis f. sp. tritici (Pst), is a devastating wheat disease in wheat producing regions worldwide (Li and Zeng, 2002;Line, 2002;Chen, 2005;Wan et al., 2007;Chen et al., 2014;Wang et al., 2014). It can have a serious impact on wheat production, and can cause yield losses of 10-30% or even no yield once a disease outbreak or epidemic occurs (Li and Zeng, 2002;Wan et al., 2004Wan et al., , 2007. Pst can mutate in many ways to produce new strains or physiological races (Ma et al., 1993;Li and Zeng, 2002;Jin et al., 2010;Hu et al., 2014;Kang et al., 2015). The emergence of new physiological races can lead to the loss of stripe rust resistance in wheat cultivars and cause periodic epidemics of the disease (Li and Zeng, 2002;Chen, 2007;Chen et al., 2014;Hu et al., 2014;Kang et al., 2015). Hence, exploring the mechanisms of virulence variation of Pst can provide an important basis for wheat breeding and control of stripe rust.
Ultraviolet-B (UV-B) radiation (280-320 nm) is a part of sunlight, and can affect living things on Earth's surface. In most reported studies, UV-B radiation has been used as a stress condition, and many aspects of its effects on living organisms have been investigated, including physiological and biochemical effects (Fargues et al., 1996;Jansen et al., 1998;Fernandes et al., 2007;Wargent and Jordan, 2013;Braga et al., 2015), genetic effects (Kumar et al., 2004;Braga et al., 2015), and effects on proteins Trentin et al., 2015;Gao et al., 2019). As the main propagules of wheat stripe rust, Pst urediospores, in high-altitude areas of northwestern China with high-intensity UV-B radiation and in processes of long-distance dispersal with upper air flows, are affected by UV-B radiation, and virulence-mutant strains may be induced (Li and Zeng, 2002;Cheng et al., 2014;Hu et al., 2014;Kang et al., 2015).
Reported studies on the effects of UV-B radiation on Pst have focused on changes in the epidemiological components of Pst after UV-B radiation (Jing et al., 1993;Cheng et al., 2014), screening of virulence-mutant strains of Pst by UV-B radiation (Shang et al., 1994;Huang et al., 2005;Wang et al., 2009;Zhao et al., 2019), and random amplified polymorphic DNA (RAPD) analysis of the UV-B-induced virulence-mutant strains of Pst (Huang et al., 2005;Wang et al., 2009). An important study on Pst (Zheng et al., 2013) revealed the whole genome information of CYR32, a dominant physiological race of Pst in China, using high-throughput sequencing technology, which provided an important basis for the identification and functional verification of Pst proteins. More recently, a proteomics method based on isobaric tags for relative and absolute quantification (iTRAQ) technology, combining the isotope labeling method and tandem mass spectrometry, has been used to identify and quantify differentially expressed proteins (DEPs) between different proteomes, particularly in studies on plant-pathogen interactions (Fu et al., 2016;Xu et al., 2017;OuYang et al., 2018). Using this methodology, Zhao et al. (2016) investigated the difference between the proteomes of the urediospores of CYR32 before and after germination. The results showed that most of the DEPs were involved in biological processes such as carbon metabolism, energy metabolism, and transport. In our previous study (Zhao et al., 2018), the urediospores of three physiological races of Pst in China-CYR31, CYR32, and CYR33-were irradiated with a dose of UV-B radiation at which the relative lethal rate of the urediospores of each physiological race was 90%. Proteomic analysis of the irradiated urediospores, using methods based on iTRAQ technology, was then performed to explore the effects of UV-B radiation on Pst at the protein level. The results showed that most of the identified DEPs were mainly involved in energy metabolism, substance metabolism, and DNA biosynthesis. However, there are no reports on the differences in the level of protein expression between the original strains of Pst and the UV-B-induced virulence-mutant strains.
In our previous study (Zhao et al., 2019), two UV-Binduced virulence-mutant strains, CYR32-5 and CYR32-61, were screened from the UV-B-irradiated urediospores of CYR32 on the seedlings of the wheat cultivar Guinong 22. In this study, after germination of the urediospores of the original strain (CYR32) and the two UV-B-induced virulence-mutant strains (CYR32-5 and CYR32-61), proteins were extracted from the germinated urediospores and germ tubes, the DEPs among the proteomes of the three Pst strains were screened by using the proteomics method based on iTRAQ technology, and then the acquired DEPs were subjected to COG (Cluster of Orthologous Groups) annotations, KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analyses, and quantitative reverse transcription PCR (qRT-PCR) assays. This study is of great significance for exploring mechanisms of Pst virulence variation, and it can provide a reference for wheat resistance breeding and control of wheat stripe rust.

Materials
In the previous study (Zhao et al., 2019), the urediospores of the Pst physiological race CYR32 were irradiated with UV-B radiation under a radiation dose (the radiation intensity was 250 µw/cm 2 , and the radiation time was 95 min) for which the relative lethal rate of urediospores was approximately 90% (i.e., the relative germination rate of urediospores was approximately 10%); then the irradiated urediospores were inoculated on the seedlings of Guinong 22 to screen virulence-mutant strains. Finally, two virulence-mutant strains named CYR32-5 and CYR32-61, with stable infection types on Guinong 22 throughout four successive generations, were obtained. In this study, the original strain of the Chinese physiological race CYR32 of Pst and the two virulence-mutant strains (CYR32-5 and CYR32-61) were used. The wheat cultivar Mingxian 169, which is susceptible to all known Chinese physiological races of Pst, was used for Pst multiplication. Multiplication of the urediospores of each Pst strain was conducted in an artificial climate chamber (environmental parameters: light time, 12 h/d; light intensity, 10,000 lux; temperature, 11-13 • C; and relative humidity, 60-70%) using the method described by Cheng et al. (2014).

Germination of Pst Urediospores
Fresh urediospores of each of the three strains CYR32, CYR32-5, and CYR32-61 were collected from the diseased leaves of Mingxian 169 seedlings. For each strain, a spore suspension with a concentration of 2.5 mg/mL was prepared with 0.2% Tween-80 solution. Then, 4 mL of the spore suspension was transferred into a Petri dish (16 cm in diameter) containing 1% water agar medium. The urediospores (10 mg) contained in the suspension were evenly scattered in the Petri dish. When the suspension was almost dry, the Petri dish was sealed with plastic wrap and incubated for 10 h at 9 • C in a dark environment. After incubation, a small agar block was picked up and placed on a glass slide for microscopic observation. At least 300 urediospores were checked, and the germination rate of the urediospores was recorded. A urediospore with a germ tube longer than half of the diameter of the urediospore was regarded as germinated. If the germination rate of the urediospores was more than 90%, the urediospores and germ tubes on the surface of the water agar medium were gently collected with a cover glass. For each Pst strain, the urediospores and germ tubes collected from five Petri dishes (approximately 50 mg in total) were treated as a sample.

Protein Extraction and iTRAQ Labeling
In total, two samples of CYR32, three samples of CYR32-5, and three samples of CYR32-61 were acquired for protein extraction. Grinding tools were pre-cooled with liquid nitrogen, and then each of the prepared samples was pulverized with liquid nitrogen in a mortar and pestle. The pulverized sample was suspended in a trichloroacetic acid (TCA) and acetone solution (TCA:acetone = 1:10, w/v; pre-cooled at −20 • C) and then precipitated for 2 h at −20 • C. After centrifugation at 4 • C under 20,000 × g for 30 min, the supernatant was discarded, and the precipitate was suspended in pre-cooled pure acetone, and then precipitated for 30 min at −20 • C prior to centrifugation at 4 • C under 20,000 × g for another 30 min. This process was repeated several times until the precipitate was substantially white. The precipitate was resuspended in a lysis buffer (8 M urea, 30 mM 4-(2-hydroxyethyl)-1piperazineethanesulfonic acid, 1 mM phenylmethanesulfonyl fluoride, 2 mM ethylenediaminetetraacetic acid, and 10 mM dithiothreitol) and sonicated for 5 min (pulse on 2 s, pulse off 3 s, power 180 W), followed by centrifugation at 20,000 × g for 30 min. The supernatant was collected and added to dithiothreitol to a final concentration of 10 mM. After incubation in a water bath at 56 • C for 1 h, iodoacetamide was quickly added, to a final concentration of 55 mM. After incubation for 1 h in a dark environment, the supernatant solution was supplemented with a fourfold volume of pre-chilled acetone, and then precipitated at −20 • C for more than 3 h. After centrifugation at 4 • C under 20,000 × g for 30 min, the precipitate was dissolved in 400 µL of digestion buffer to a final concentration of 50% triethylammonium bicarbonate (TEAB) and 0.1% sodium dodecyl sulfate (SDS). Subsequently, sonication was performed for 3 min with a pulse on time of 2 s and a pulse off time of 3 s at an ultrasonic power of 180 W, and centrifugation was performed for 30 min at 4 • C and 20,000 × g. Finally, the supernatant was obtained, and the protein concentration was quantified using the Bradford assay (Supplementary Table 1).
From each sample solution, 100 µg of protein was transferred to a new clean centrifuge tube. The TEAB solution (0.1% SDS) was added to make the protein solution of each sample in the new tube up to the same volume. After adding 3.3 µL of trypsin (1 µg/µL) to the tube for protein digestion, the solution was incubated in a water bath at 37 • C for 24 h. Subsequently, the tube was supplemented with 1 µL of trypsin (1 µg/µL), and the mixture in the tube was incubated in a water bath at 37 • C for 12 h. After lyophilization of the mixture, 30 µL of TEAB (ddH 2 O:TEAB = 1:1, v/v) was added to the tube to dissolve the peptides. Using an iTRAQ R Reagent-8Plex Multiplex Kit (Applied Biosystems, Foster City, CA, United States), peptides from the eight samples were labeled with the iTRAQ tags as follows: CYR32 (tags 113 and 114), 116,and 117), and CYR32-61 (tags 118, 119, and 121), respectively (Supplementary Table 1).

Proteomic Data Analysis
The acquired raw mass spectrometry data were processed using the PD (Proteome Discoverer 1.3; Thermo Fisher Scientific, San Jose, CA, United States) software with the following parameters: mass range of parent ion, 350-6,000 Da; minimum number of peaks in MS/MS spectrum, 10; signal-to-noise ratio (S/N) threshold, 1.5. The spectra extracted using the PD software were searched against the Basidiomycota_UniProt database using Mascot 2.3.0 (Matrix Science, London, United Kingdom) with the following identification parameters: fixed modification, carbamidomethyl (C); variable modification, oxidation (M), Gln→Pyro-Glu (N-term Q), iTRAQ 8 plex (K), iTRAQ 8 plex (Y), iTRAQ 8 plex (N-term); peptide tolerance, 15 ppm; MS/MS tolerance, 20 mmu; max missed cleavages, 1; enzyme, trypsin. Based on the Mascot search results and the extracted spectra, protein quantitative analysis was performed using the PD software with the following parameters: protein ratio type, median; minimum peptides, 1; normalization method, median; P-value, < 0.05; ratio, ≥ 1.2. In this study, proteins with at least two unique peptides, scores more than 70, P-values less than 0.05, and fold changes more than 1.2 or less than 0.83 were identified as DEPs between the different Pst strains. The fold changes of the up-regulated DEPs were more than 1.2, and those of the down-regulated DEPs were less than 0.83. For the convenience of expression, CYR32-5 vs. CYR32, CYR32-61 vs. CYR32, or CYR32-61 vs. CYR32-5 was used to represent the combination in which the proteome of the former was compared with that of the latter to screen DEPs.
Functional annotations of the proteins were carried out using the Blast2Go program 1 (Conesa et al., 2005). The COG annotations of proteins were conducted using WebMGA 2 . The annotations of the identified DEPs at the biological pathway level were performed using the KEGG database 3 . The KOBAS 2.0 software was used to annotate pathways by comparing similar protein sequences (Xie et al., 2011). Signal peptides were predicted and analyzed using the SignalP 5.0 program 4 (Armenteros et al., 2019).
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium 5 via the iProX partner repository  with the dataset identifier PXD018136.

Relative Quantification of mRNA by qRT-PCR
Using the method described above, 20 mg of the germinated urediospores and germ tubes of each of the three strains CYR32, CYR32-5, and CYR32-61 were collected and treated as a sample. Total RNA from each sample was extracted using the RNeasy kit (Omega) according to the manufacturer's instructions, and then genomic DNA contaminants were removed by DNase I (RNase-free) treatment. The synthesis of the first-strand cDNA was performed by reverse transcription using the PrimeScript TM RT Master Mix kit (Takara, Dalian, China). Specific primers were designed for the selected DEPs according to their corresponding coding sequences using NCBI Primer-BLAST 6 ( Table 1). Pst β-tubulin gene TUBB (GenBank accession No. EG374306) was used as the reference gene, and the corresponding primer reported by Huang et al. (2012) (as shown in Table 1) was used in this study. SYBR Green real-time fluorescence quantitative TABLE 1 | Sequences of the primers used for real-time fluorescence quantitative PCR assays.

Accession
Forward primer (5 -3 ) Reverse primer (5 -3 ) CAGAACCAGGGTGAGATGCT CACCCTCATCGGCGTAAGTA PCR assays were performed using an Applied Biosystems ABI Model 7500 Real Time PCR system. The qRT-PCR reactions were carried out in a total 20 µL volume of reaction mixture containing 10 µL 2 × SYBR Premix DimerEraser, 0.6 µL forward primer (10 µM), 0.6 µL reverse primer (10 µM), 2.0 µL template DNA, 6.4 µL sterile purified water, and 0.4 µL 50 × ROX Reference Dye or Dye II. The amplification procedure was as follows: predenaturation at 95 • C for 30 s, followed by 40 cycles of 95 • C for 5 s, 55 • C for 30 s, and 72 • C for 30 s. The relative gene expression levels were calculated using the 2 − C T method described by Livak and Schmittgen (2001). The experiment was carried out with three biological replicates.

Overview of the Proteomic Data of the Three Pst Strains
A total of 337,731 spectra were acquired from the eight samples of the three Pst strains, and 52,290 spectra were matched by searching against the Basidiomycota_UniProt database using Mascot 2.3.0 (false discovery rate < 1%). In total, 8,882 peptides and 2,271 proteins were obtained (Supplementary Table 2).

Acquired DEPs and COG Annotations
The results of screening the DEPs in the urediospores and germ tubes after germination in three combinations, CYR32-5 vs. CYR32, CYR32-61 vs. CYR32, and CYR32-61 vs. CYR32-5, are shown in Figure 1. As shown in Figure 1, 59 DEPs were obtained in the combination CYR32-5 vs. CYR32, of which 35 proteins were up-regulated and 24 proteins were down-regulated, and the DEPs E3JUB4 and Q9C1C1 (3.39%) were predicted to contain signal peptides. Seventy-four DEPs were obtained in the combination CYR32-61 vs. CYR32, of which 43 proteins were up-regulated and 31 proteins were down-regulated, and the DEPs Q9C1C1, R7SW21, E3JQF7, E3JUB4, E3JSA3, E3KXP2, and E3K3X8 (9.46%) were predicted to contain signal peptides. Sixty-four DEPs were obtained in the combination CYR32-61 vs. CYR32-5, of which 40 were up-regulated and 24 were downregulated, and the DEPs E3K302, E3KDD8, E3JSA3, and E3JQF7 (6.25%) were predicted to contain signal peptides. In total, 144 DEPs were obtained in the three comparison combinations, of which four were common in CYR32-5 vs. CYR32, CYR32-61 vs. CYR32, and CYR32-61 vs. CYR32-5. The DEPs that were predicted to contain signal peptides may be secretary proteins or membrane proteins. All obtained DEPs were subjected to COG functional annotation, and then assigned to 25 categories that were represented by A-Z (Figure 2). Among all the DEPs, seven were annotated with multiple functions, and 22 were not annotated and had no annotation information. The four main functional categories were as follows: [J] translation, ribosomal structure and biogenesis (31, 21.68%), [O] posttranslational modification, protein turnover, chaperones (24, 16.78%), [C] energy production and conversion (17, 11.89%), and [U] intracellular trafficking, secretion, and vesicular transport (11, 7.69%). The detailed functional annotations of the upregulated DEPs and the down-regulated DEPs are shown in Tables 2, 3, respectively. The acquired DEPs were classified into four major functional categories, including "information storage and processing, " "cellular processes and signaling, " "metabolism, " and "poorly characterized." Of the up-regulated DEPs, 14, 43, 27, and 4 were involved in "information storage and processing, " "cellular processes and signaling, " "metabolism, " and "poorly characterized, " respectively, and 18 had no annotation information. Of the down-regulated DEPs, 23, 15, 18, and 5 were involved in "information storage and processing, " "cellular processes and signaling, " "metabolism, " and "poorly characterized, " respectively, and 7 had no annotation information.

KEGG Pathway Enrichment Analysis of DEPs
The biological pathways that the DEPs were involved in were investigated using the KEGG database, and 5, 7, and 7 significant enriched KEGG pathways (P < 0.05) were acquired for the DEPs of CYR32-5 vs. CYR32, CYR32-61 vs. CYR32, and CYR32-61 vs. CYR32-5, respectively (Figure 3 and Supplementary Table 3). The DEPs of CYR32-5 vs. CYR32 were involved in pathways including metabolic pathways, biosynthesis of secondary metabolites, carbon metabolism, ribosome, and pyruvate metabolism. The DEPs of CYR32-61 vs. CYR32 were mainly involved in the carbon metabolism pathway, calcium signaling pathway, and melanogenesis pathway. The DEPs of CYR32-61 vs. CYR32-5 were mainly involved in the ribosome pathway, melanogenesis pathway, and calcium signaling pathway.

Transcriptional Analysis of DEPs
To confirm the differences in protein abundance in each of the three combinations including CYR32-5 vs. CYR32, CYR32-61 vs. CYR32, and CYR32-61 vs. CYR32-5, after germination of the urediospores, qRT-PCR was used to investigate the relative expression levels of the genes encoding the DEPs identified by iTRAQ analysis. The relative expression levels of the genes encoding six proteins (G9B235, F4S8X5, V2XWY1, Q9C1C1, F4S0Z3, and F4RWN9) in CYR32-5 vs. CYR32, seven proteins (G9B235, F4S8X5, V2XWY1, Q9C1C1, E3L0W8, D4QFJ2, and J6EXB0) in CYR32-61 vs. CYR32, and one protein (F4S8X5) in CYR32-61 vs. CYR32-5 were determined, and the results are shown in Figure 4. In CYR32-5 vs. CYR32, the relative expression levels of the genes encoding Q9C1C1, G9B235, F4S8X5, F4RWN9, and F4S0Z3 were not consistent with the corresponding protein abundance based on iTRAQ data, but the relative expression of the gene encoding V2XWY1 was consistent with the corresponding protein abundance. In CYR32-61 vs. CYR32, the relative expression levels of the genes encoding Q9C1C1 and J6EXB0 were inconsistent with the corresponding protein levels, and the genes encoding G9B235, F4S8X5, V2XWY1, E3L0W8, and D4QFJ2 were consistent with the corresponding protein levels. In CYR32-61 vs. CYR32-5, the relative expression level of the gene encoding F4S8X5 at the transcriptional level was consistent with the corresponding protein abundance obtained via iTRAQ analysis.

DISCUSSION
In this study, using the proteomics method based on iTRAQ technology, a total of 144 DEPs were obtained from the germinated urediospores with germ tubes in the CYR32-5 vs. CYR32, CYR32-61 vs. CYR32, and CYR32-61 vs. CYR32-5 conditions. Nine DEPs, including E3JQF7, E3JSA3, E3JUB4, E3K302, E3K3X8, E3KDD8, E3KXP2, Q9C1C1, and R7SW21, were predicted to contain signal peptides. Peptidyl-prolyl cistrans isomerase (PPI, E3JQF7) is a high-efficiency foldase. In biological cells, after ribosome transcription, the peptide chain functions properly only through correct folding and positioning. PPI is associated with the correct folding of a protein (Lu et al., 2007). Heat shock protein 70 (R7SW21) can play an important role in the folding and stretching of newly synthesized polypeptides, can repair denatured proteins, and can prevent protein denaturation by up-regulating protein expression under adverse conditions (Tsan and Gao, 2004). In the present study, both E3JQF7 and R7SW21 were up-regulated in CYR32-61 vs. CYR32, and E3JQF7 was up-regulated in CYR32-61 vs. CYR32-5, indicating that strain CYR32-61 may be more resistant to stress than strains CYR32 and CYR32-5.
The histone chaperone ASF1 (anti-silencing function 1) is a molecular chaperone of histone H3-H4. It has been reported that ASF1 plays important roles in the processes of DNA replication and repair in yeast (Le et al., 1997) and Arabidopsis . In this study, ASF1 (E3L109) was up-regulated in CYR32-61 vs. CYR32-5, indicating that strain CYR32-61 may be more resistant to stress than strain CYR32-5.
G proteins are involved in various cellular activities, including regulation of intracellular Ca 2+ concentration, protein synthesis, and the binding of ribosomes to the endoplasmic reticulum (Lagerstedt et al., 2005). The DEPs G9B235, H6QRQ7, and J9PGK7 acquired in this study were involved in G protein regulation. G9B235 was down-regulated in CYR32-5 vs. CYR32, J9PGK7 was up-regulated in CYR32-61 vs. CYR32-5, and H6QRQ7 and J9PGK7 were up-regulated, and G9B235 was down-regulated in CYR32-61 vs. CYR32. Furthermore, the relative expression level of G9B235 determined by qRT-PCR was consistent with the corresponding protein abundance obtained by iTRAQ technology in CYR32-61 vs. CYR32, but the former was inconsistent with the latter in CYR32-5 vs. CYR32. Therefore, the results indicated that there was little correlation between the expression levels at the transcriptional level and the protein level for some DEPs, which was consistent with a previous report by Zhao et al. (2016).
As an important component of the second messenger system, calmodulin plays a key role in the regulation of the calcium signal system, and is associated with physiological metabolism regulation, gene expression, and normal growth and development of cells (Chin and Means, 2000). In both CYR32-61 vs. CYR32 and CYR32-61 vs. CYR32-5, putative uncharacterized protein (E3L0W8) with the role of calcium ion binding and calmodulin (M5FTW6) were up-regulated, and the relative expression level of E3L0W8 determined by qRT-PCR was consistent with the protein level obtained by iTRAQ technology. Calmodulin (E3KLJ3) was up-regulated in CYR32-61 vs. CYR32-5.
In CYR32-5 vs. CYR32, three DEPs involved in energy production, including ubiquinol-cytochrome c reductase cytochrome c1 subunit (E3JWK6), ubiquinol-cytochrome c reductase iron-sulfur subunit (E3L519), and ATP synthase subunit beta (F4S0Z3), were up-regulated. E3JWK6 and E3L519 are components of complex III of the electron transport chain, and are associated with electron transport. F4S0Z3 is the β subunit of the F1 part of ATP synthase, and contains a site to catalyze ATP synthesis. In living organisms, the electron   transport chain and ATP synthesis are coupled together (Mitchell, 1961). A model for the coupling of electron transport and ATP synthesis in the combination CYR32-5 vs. CYR32 is shown in Figure 5. In this study, the acquired DEPs involved in the electron transport chain and in ATP synthesis were up-regulated, indicating that the processes of energy production during urediospore germination may be promoted in strain CYR32-5 compared to strain CYR32. In the glycolysis pathway, glyceraldehyde-3-phosphate dehydrogenase can catalyze the formation of 1,3diphosphoglycerate from 3-phosphoglyceraldehyde in the presence of NAD + and phosphoric acid; it is, hence, an important enzyme in glycolysis. In this study, D4QFJ2, annotated as glyceraldehyde-3-phosphate dehydrogenase, was up-regulated in CYR32-61 vs. CYR32 according to iTRAQ analysis, which was consistent with the relative expression level determined by qRT-PCR. The results indicated that, during urediospore germination, the glycolysis pathway (as shown in Figure 6) of virulence-mutant strain CYR32-61 may be promoted in comparison with that of strain CYR32. Glycolysis is the main process for the formation of ATP in carbohydrate metabolism. The results obtained in this study showed that the DEPs involved in the glycolysis pathway were up-regulated, further indicating that, during urediospore germination of the virulence-mutant strain CYR32-61 (as compared to the original strain CYR32), polysaccharide may be utilized to form ATP to promote the biological processes of urediospores. In addition, during glycolysis, the final product, pyruvate, can be irreversibly catalyzed by pyruvate dehydrogenase multi-enzyme complex (PDHc) to form acetyl-coenzyme A, which can be oxidatively decomposed in the tricarboxylic acid cycle. Pyruvate dehydrogenase E1   (Danson et al., 1978). In CYR32-5 vs. CYR32, the expression of M5FPE0-annotated as pyruvate dehydrogenase E1 component alpha subunit-was up-regulated, indicating that the formation of ATP during oxidative phosphorylation in the germinated urediospores and germ tubes may be promoted in strain CYR32-5 in comparison with strain CYR32. Inorganic diphosphatase (A8N2Q4) mainly involved in DNA synthesis, can catalyze the conversion of one molecule of pyrophosphate into two molecules of phosphate, accompanied by the generation of high levels of energy (Kukko and Saarento, 1984). A8N2Q4 was up-regulated in CYR32-61 vs. CYR32 and FIGURE 3 | KEGG pathway enrichment of the differentially expressed proteins acquired in the combinations CYR32-5 vs. CYR32, CYR32-61 vs. CYR32, and CYR32-61 vs. CYR32-5 (P < 0.05). FIGURE 5 | Model for the coupling of electron transport and ATP synthesis. The differentially expressed proteins (marked in red) in the electron transport chain complex III, including ubiquinol-cytochrome c reductase cytochrome c1 subunit (E3JWK6) and ubiquinol-cytochrome c reductase iron-sulfur subunit (E3L519), were up-regulated in CYR32-5 vs. CYR32. The β subunit (F4S0Z3) of ATP synthase marked in red was up-regulated in CYR32-5 vs. CYR32.
FIGURE 6 | Glycolysis pathway. The glyceraldehyde-3-phosphate dehydrogenase marked in red was found to be up-regulated in CYR32-61 vs. CYR32.
CYR32-61 vs. CYR32-5 in this study. The results indicated that in comparison with strains CYR32 and CYR32-5, the energy metabolism pathway during urediospore germination may be promoted in strain CYR32-61, which may contribute to the progress in biological processes of urediospores.
It has been reported that UV-B radiation can induce virulence variation in Pst (Shang et al., 1994;Huang et al., 2005;Wang et al., 2009;Zhao et al., 2019). Similarly, in studies on two other important pathogens-Puccinia triticina (Neugebauer et al., 2018) and P. graminis f. sp. tritici (Zhang et al., 2017), causing wheat leaf rust and wheat stem rust, respectively, it has also been reported that mutation may be the main mechanism of virulence variation of the two pathogens, thus resulting in the loss of wheat resistance and severe losses of wheat yield. In this study, after the germination of urediospores, the expression levels of many proteins in UV-B-induced virulence-mutant strains, CYR32-5 and CYR32-61, were different from those in CYR32, the original strain. The changes in the related biological processes in which the DEPs are involved may result in changes in virulence phenotypes of Pst and may be responsible for virulence variation. To further explore the virulence variation mechanisms of Pst, it is of great importance to conduct studies by pathogen sequencing at the genome level and to investigate the metabolomics of Pst after mutations.

CONCLUSION
In this study, using the proteomics method based on iTRAQ technology, a quantitative proteomic analysis of the germinated urediospores (with germ tubes) of the original strain, CYR32, and two UV-B-induced virulence-mutant strains, CYR32-5 and CYR32-61, was undertaken. A total of 2,271 proteins were identified in eight samples of the three Pst strains, and 59, 74, and 64 DEPs were obtained in each of the three combinations, i.e., CYR32-5 vs. CYR32, CYR32-61 vs. CYR32, and CYR32-61 vs. CYR32-5, respectively. The identified DEPs were mainly involved in energy metabolism, carbon metabolism and cellular substance synthesis. The relative expression levels of the genes encoding some DEPs, quantified using qRT-PCR, were consistent with the corresponding protein abundance obtained via iTRAQ analysis. Compared with the original strain, CYR32, the DEPs involved in stress-related and energy metabolism were up-regulated in the virulence-mutant strains, indicating that the virulence-mutant strains, CYR32-5 and CYR32-61, were more tolerant to stress than the original strain. These results are of great significance for further studies on the mechanisms of Pst virulence variation and for implementing control measures for wheat stripe rust.

DATA AVAILABILITY STATEMENT
The datasets generated for this study can be found in the ProteomeXchange Consortium (http://proteomecentral. proteomexchange.org) via the iProX partner repository with the dataset identifier PXD018136.