Next Article in Journal
Polar Lipids in Starch-Rich Commodities to be Analyzed with LC-MS-Based Metabolomics—Optimization of Ionization Parameters and High-Throughput Extraction Protocols
Previous Article in Journal
Microbial Transformations of Organically Fermented Foods
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 9
Issue 8
10.3390/metabo9080166
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Differential Proteomics Based on TMT and PRM Reveal the Resistance Response of Bambusa pervariabilis × Dendrocalamopisis grandis Induced by AP-Toxin

by
Qianqian He
,
Xinmei Fang
,
Tianhui Zhu
,
Shan Han
,
Hanmingyue Zhu
and
Shujiang Li
*
College of Forestry, Sichuan Agricultural University, No. 211, Huimin Road in Wenjiang District, Chengdu 611130, China
*
Author to whom correspondence should be addressed.
Metabolites 2019, 9(8), 166; https://doi.org/10.3390/metabo9080166
Submission received: 30 June 2019 / Revised: 7 August 2019 / Accepted: 8 August 2019 / Published: 10 August 2019
Browse Figures
Review Reports Versions Notes

Abstract

:
Bambusa pervariabilis McClure × Dendrocalamopsis grandis (Q.H.Dai & X.l.Tao ex Keng f.) Ohrnb. blight is a widespread and dangerous forest fungus disease, and has been listed as a supplementary object of forest phytosanitary measures. In order to study the control of B. pervariabilis × D. grandis blight, this experiment was carried out. In this work, a toxin purified from the pathogen Arthrinium phaeospermum (Corda) Elli, which causes blight in B. pervariabilis × D. grandis, with homologous heterogeneity, was used as an inducer to increase resistance to B. pervariabilis × D. grandis. A functional analysis of the differentially expressed proteins after induction using a tandem mass tag labeling technique was combined with mass spectrometry and liquid chromatography mass spectrometry in order to effectively screen for the proteins related to the resistance of B. pervariabilis × D. grandis to blight. After peptide labeling, a total of 3320 unique peptides and 1791 quantitative proteins were obtained by liquid chromatography mass spectrometry analysis. Annotation and enrichment analysis of these peptides and proteins using the Gene ontology and Kyoto Encyclopedia of Genes and Genomes databases with bioinformatics software show that the differentially expressed protein functional annotation items are mainly concentrated on biological processes and cell components. Several pathways that are prominent in the Kyoto Encyclopedia of Genes and Genomes annotation and enrichment include metabolic pathways, the citrate cycle, and phenylpropanoid biosynthesis. In the Protein-protein interaction networks four differentially expressed proteins-sucrose synthase, adenosine triphosphate-citrate synthase beta chain protein 1, peroxidase, and phenylalanine ammonia-lyase significantly interact with multiple proteins and significantly enrich metabolic pathways. To verify the results of tandem mass tag, the candidate proteins were further verified by parallel reaction monitoring, and the results were consistent with the tandem mass tag data analysis results. It is confirmed that the data obtained by tandem mass tag technology are reliable. Therefore, the differentially expressed proteins and signaling pathways discovered here is the primary concern for subsequent disease resistance studies.

1. Introduction

When a plant is induced by a biologically inducing factor (pathogenic toxin, etc.), the plant will change the expression of the proteins in its cells, thereby changing the activity of related enzymes to respond to the induction. In order to understand the interactions between biologically inducible factors and host plants, we need to conduct in-depth research based on proteomics. Currently, many proteomic studies focus on the response of plants to abiotic factors, such as the effects of drought and salt content, etc. [1,2].
Proteomics is the study of all proteins expressed by a cell under specific conditions. This process includes the study of protein expression levels, post-translational modifications, and interactions between proteins and other proteins or with other biomolecules. Proteomics can thus reveal the molecular mechanisms involved in a biological process or disease at the protein level [3]. At present, isobaric tags for relative and absolute quantification (iTRAQ)/tandem mass tag (TMT) technology includes a wide range of methods for proteomics research [3,4]. The general technique utilizes four or eight isotopically encoded tags that specifically label the amino group of a peptide and then uses tandem mass spectrometry to simultaneously compare the relative or absolute content of the protein in four or eight different samples [3]. González-Fernández et al. [5] reviewed the progress of proteomics in plant fungal pathogens research, which is an excellent tool that can give us a great deal of information about fungal pathogenicity by high-throughput studies. This approach has allowed the identification of new fungal virulence factors by characterizing signal transduction or biochemical pathways and studying the fungal life cycle and their life-style. Similarly, in the plant-fungus intricacies, proteomics provides rapid insight and is expected to be one of the imminent and integrative tools in biological research [6]. Rustagi et al. [6] summarized more than 17 important plant-fungal pathogen studies that used proteomics by methods including Two-Dimensional Polyacrylamide Gel Electrophoresis (2D-PAGE), Fluorescent Two-Dimensional Difference Gel Electrophoresis (2D-DIGE), Isotope-Coded Affinity Tags (ICAT), iTRAQ, Multidimensional Protein Identification Technology (MudPIT), and Mass Spectrometry (MS).
With the developing breadth of research in plant-induced disease resistance, the idea of using a pathogen toxin as an induction factor for disease resistance has garnered more attention and recognition. Many studies around the world have incorporated this idea. Yang et al. [7] explored the toxin extracted from Phytophthora infestans to induce potato resistance to Phytophthora infestans. Hui et al. [8] used the crude toxin of leaf mold to induce the changes of defense enzymes and reactive oxygen species in tomato seedlings. Ma et al. [9] reported a toxin-induced maize resistance to Bipolaris maydis. They found an increase in the amount of various defense enzymes, such as peroxidase (POD) and phenylalanine ammonia-lyase (PAL) in the maize seedlings after treatment with the inducer, indicating that the plant resistance is improved. Conrath et al. [10] successfully induced partial potato resistance using the incompatible races of Phytophthora infestans and showed that the mechanism of action of disease resistance is that the potato produces phytoalextin. Therefore, this phytoalexin has provided the basis for the study of plant induced disease resistance. However, in the area of forest devastating disease research, there are few studies into the use of toxins as inducers for disease resistance, and the mechanism of induction, studied by using proteomics, has not been established.
B. pervariabilis × D. grandis blight is a widespread and dangerous forest fungus disease. In many provinces of China, large-scale B. pervariabilis × D. grandis death has occurred due to blight, in an area spanning more than 3000 hm2. This disease greatly threatens the process of returning farmland to forests and the construction of ecological barriers [11]. B. pervariabilis × D. grandis blight has been listed as a supplementary object of forest phytosanitary measures [12]. The pathogen A. phaeospermum belongs to Fungi, Dikarya, Ascomycota, Pezizomycotina, Sordariomycetes, Xylariomycetidae, Xylariales, Apiosporaceae, Arthrinium [13].
In this study, the purified toxin of A. phaeospermum (AP-Toxin) was obtained by Li et al. [14,15,16]. Through the results of a previous study, the best inducer was the purified pathogen toxin inactivated at 60 °C, and the concentration was 40 μg/mL [17]. The TMT technique combined with liquid chromatography mass spectrometry (LC-MS/MS) were used to explore the differential expression of proteins and metabolic pathways induced by the toxin and then inoculated by A. phaeospermum in B. pervariabilis × D. grandis. This work provides a theoretical basis for the in-depth study of the plant’s resistance under the action of biological inducers and lays a theoretical foundation for the management of forest diseases.

2. Results

2.1. Symptoms of B. Pervariabilis × D. Grandis after AP-Toxin Induction

B. pervariabilis × D. grandis that was induced with AP-toxin and inoculated with A. phaeospermum shows significant visual differences when compared with the control (sterile water + A. phaeospermum) (Figure 1). After treatment with AP-Toxin + A. phaeospermum, the leaves and branches of B. pervariabilis × D. grandis show a healthy state without symptoms of yellowing and wilting (Figure 1A). Most of the B. pervariabilis × D. grandis leaves treated with sterile water + A. phaeospermum have yellowing symptoms, and the branches also show signs of withering (Figure 1B). It is thus shown that the inducer can enhance the ability of the B. pervariabilis × D. grandis to resist A. phaeospermum.

2.2. Quantitative Proteomic Response to Inducer

TMT technology for protein enzymatic hydrolysis, peptide quantification and labeling, strong cation exchange (SCX) fractionation, and mass spectrometry were used to identify peptides, quantify proteins, and analyze differentially expressed proteins. According to statistical analyses of the results (Table 1), 3320 unique peptides were identified, and 1791 proteins were quantified. The percentage of total protein covered by two or more different peptides reached 45.5%, indicating that the selected TMT markers combined with mass spectrometry can effectively separate and identify proteins from the B. pervariabilis × D. grandis experimental and control groups. A differentially expressed protein was defined to be one that had a change in expression greater than 1.2-fold (up-and-down), and the significance P value levels less than 0.05 were used as a standard for screening differentially expressed proteins. The screened differentially expressed proteins were blasted by National Center for Biotechnology Information (NCBI). There were 163 differentially expressed proteins, of which 75 were up-regulated and 88 were down-regulated.

2.3. GO and KEGG Analysis of Differentially Expressed Proteins

The Gene Ontology (GO) functional analysis of proteins was quantified for B. pervariabilis × D. grandis considering three aspects: biological processes, molecular functions, and cell components [18] (Figure 2). The results show that a significant enrichment of biological processes mainly affects the single-organism metabolic process (GO:0044710; 80 proteins; p value_0.00000000141) and carbohydrate metabolic process (GO:0005975; 38 proteins; p value_0.000000339). They all appear to activate the up-regulation. Significant enrichment of the molecular function mainly concerns catalytic activity (GO:0003824; 101 proteins; p value_0.000319) and cation binding (GO:0043169; 55 proteins; p value_0.00000316). Finally, the affected cellular components are mainly composed of an intracellular organelle part (GO:0044446; 61 proteins; p value_0.0000000309) and chloroplast (GO:0009507; 59 proteins; p value_0.000000000000000162) (Table S1). The bioaccumulated entries show that the multiple biological pathways of B. pervariabilis × D. grandis that were significantly changed after induction agent treatment are related to plant stress tolerance: including response to osmotic stresses, external biotic stimulus, oxidative stresses, and so on. In addition, the carbohydrate metabolic process is closely related to energy metabolism. They all appear to activate the up-regulation. Carbohydrate and energy adjustment respond to plants’ responses to pathogen stress, providing the necessary dynamic support through the up-regulation of energy metabolism-related proteins to the resistance stress pathway.
Differential proteins are classified into A–F groups, including A (Metabolism), B (Genetic information processing), C (Environmental information processing), D (Cellular processes), E (Organismal systems), and F (Human diseases). The Kyoto Encyclopedia of Genes and Genomes (KEGG) classification results indicate that the differentially expressed proteins are mainly concentrated in metabolic processes (A metabolism) (Figure 3). The metabolic pathways with the highest number of differentially expressed proteins are metabolic pathways (map01100; 53 proteins; p value_0.000123), biosynthesis of secondary metabolites (map01110; 39 proteins; p value_0.0000117), and phenylpropanoid biosynthesis (map00940; 13 proteins; p value_0.0000346) (Figure 3, Table S2). The phenylpropanoid biosynthesis is related to plant defense reactions. Through a series of reactions, polyphenols such as lignin, flavonoids, and flavonols, which can be synthesized can effectively prevent the invasion of various pathogens and participate in the interaction process between plants and pathogens. Moreover, this signaling is associated with plant immune responses.

2.4. PPI Analysis

In the protein-protein interaction (PPI) network map (used for establishing protein networks) (Figure 4), the most significantly enriched proteins, with high connectivity, are shown in red, and the significant enrichment of the metabolic pathway is shown in dark blue. The three metabolic pathways that were significantly enriched were Biosynthesis of secondary metabolites, Phenylpropanoid biosynthesis, and Metabolic pathways. There are 16 differential proteins directly related to the Biosynthesis of secondary metabolites pathway, nine of which are up-regulated; three differential proteins are directly linked to the Phenylpropanoid biosynthesis pathway, two of which are up-regulated; 22 differential proteins are linked to it in the Metabolic pathways pathway, 11 of which are up-regulated. A total of 11 differential proteins were up-regulated in this network map, and it was found that PH01000123G1100, PH01000159G0130, PH01000713G0340, PH01000761G0570, and PH01003309G0170 were associated with plant resistance expression, both of which showed significant up-regulation and at the same time affected multiple significantly enriched metabolic pathways (Table 2).
Round nodes represent proteins/genes (for fold change analysis, red shows up-regulation and green shows down-regulation). Rectangular nodes represent KEGG pathways/biological processes.
A significant p value is represented by a yellow–blue gradient, where yellow indicates a small p value and blue indicates a large p value.

2.5. Mass Spectrometry Analysis of PRM Validation Candidate Peptides

Parallel reaction monitoring (PRM, a sort of Multiple Reaction Monitoring for peptides) mass spectrometry was used to verify the target peptides found in the TMT analysis of the AP-toxin treated B. pervariabilis × D. grandis, and the 14 candidate peptides of the target protein were subjected to LC-PRM/MS analysis. The Skyline analysis results of each candidate peptide are shown in Figure S1, Table S3 and include the chromatographic peak contrast map (Skyline analysis map) of each peptide segment in different samples. After the Skyline analysis of 3–5 sub-ions with high abundance and as continuous as possible in the secondary mass spectrometry of candidate peptide segments, the peak area of each target peptide segment (Table S4) was obtained. This included the name of the target protein, sequence of peptide segments, charge number of parent ions, selected sub-ions, charge number of the sub-ions, and the original peak area of each sub-ion for quantification. The sum of all ion peak areas of each peptide segment is the ion peak area of the peptide segment. The ionic peak area of the normalized (N) corresponding peptide segment is presented in Table S5. The normalized peak area of the peptide segment was used to quantitatively analyze the target peptide segments in different samples, and the results are shown in Table S6.

2.6. PRM Verification of the Top Candidates for Differentially Expressed Proteins

A quantitative analysis using LC-PRM/MS of the expression levels of 14 candidate proteins associated with response-induced resistance was used to validate the results of TMT. The PRM quantitative results (Table 3) show similar trends to TMT for 13 candidate proteins, indicating that the data obtained by this experimental TMT technique combined with LC/MS is reliable. PH01000087G1730, PH01000860G0570, PH01000898G0600, PH01001064G0040, PH01001724G0160, PH01001918G0090, PH01002295G0170, and PH01100083G0010 exhibit PRM results close to the TMT results; however, PH01000159G0130, PH01000713G0340, PH01000761G0570, and PH01003309G0170 show multiple differences from TMT. It is possible that LC-PRM/MS quantitative analysis technology has higher sensitivity and resolution for these proteins [19].

3. Discussion

TMT-labeled nanoscale liquid chromatography tandem mass spectrometry (LC-MS/MS), as an emerging quantitative method technology, overcomes the shortcomings of traditional methods that cannot quantify macromolecules and proteins [20]. We used TMT technology combined with LC/MS analysis for the first time to study forest-induced resistance. The main results identified 163 differentially expressed proteins by three independent experiments, of which 75 were up-regulated and 88 were down-regulated. GO annotation and enrichment analysis of the differentially expressed proteins determined the biological process, cellular component, or molecular function of the various proteins/genes. The results indicate that the differentially expressed proteins are mainly concentrated in biological processes and cellular components. KEGG pathway annotation and enrichment analysis indicates that the differential proteins induced by AP-toxin in B. pervariabilis × D. grandis mainly enriched the metabolic pathways such as starch and sucrose metabolism, phenylpropanoid biosynthesis, the citrate cycle (TCA cycle), and amino acid transport and metabolism. The key differential proteins for screening response inducers are sucrose synthase, ATP-citrate synthase beta chain protein 1, peroxidase, and phenylalanine ammonia-lyase. By using PRM combined with LC/MS quantitative analysis to verify TMT results, 14 related proteins were screened from KEGG and GO annotation and functional classification based on the influence of plant resistance. The results showed that 13 candidate proteins showed similar trends to TMT, indicating TMT experiment. The data is reliable.
Inducing disease resistance mainly defends against disease by activating the resistance mechanism in plants and enhances preventive measures and controllability. It is an important method for current plant disease prevention and control [21,22]. Ping [23] used iTRAQ to comprehensively analyze the protein profile induced in Arabidopsis thaliana root, where an abundance of proteins related to the phenylpropanoid pathway, oxidative stress, and respiratory processes were highly increased. Jing et al. [24] studied induced sweet orange and found that the proteins that were up-regulated were related to the defense response. Lippert et al. [25] studied the protein profile of Norwegian spruce under chitosan induction and found proteins associated with oxidative stress response. Stresses in plants after induction include significant changes in various protein levels, such as pathogenesis-related proteins, protease inhibitors, polyphenol oxidase, peroxidase, lipoxygenase, and phenylalanine ammonia-lyase [26,27,28,29]. Furthermore, differentially expressed proteins such as sucrose synthase in the energy-related glucose metabolism and citrate synthase in the citric acid cycle are also significantly upregulated [30,31]. Similar to the above results, in this research, differentially expressed proteins involved in carbohydrate and energy formation processes were up-regulated, and differential proteins involved in signal sensing and amino acid biosynthesis were also significantly expressed. The adjustment of carbohydrates and energy is associated with plant response to pathogen stress [32,33]. Signal-aware metabolic pathways respond to plant resistance to pathogens [34]. The synthesis of amino acids is also considered to be an important mechanism of plant stress response and regulation [35,36]. In summary, B. pervariabilis × D. grandis was induced to activate the defense and energy metabolism-related pathways, and the necessary protein support was provided by the up-regulation of energy metabolism-related proteins. Therefore, the major differential proteins involved in these metabolic pathways are considered to be key proteins in response to inducer action.
Among the starch and sucrose metabolism and citrate cycle (TCA cycle) pathways, the differentially expressed proteins are sucrose synthase and ATP-citrate synthase beta chain protein 1. Sucrose synthase can catalyze the conversion of sucrose to uridine diphosphate glucose (UDPG) and fructose [37], a widely occurring UDP-glycosyltransferase and a key enzyme that catalyzes sucrose metabolism [38,39]. Glycogen is an important energy storage material for synthetic bio-organisms [40]. It is involved in the biosynthesis of tannins, cellulose, starch, glycolipids and glycoproteins [37], providing the energy and carbon skeleton needed for plant growth and regulating the gene expression of the signaling molecule [38]. Thus, glycogen participates in the response of plant stress against the infection of pathogenic bacteria [41]. The TCA cycle is linked to the glycolysis pathway and participates in the carbon metabolism pathway, which plays a vital role in the energy metabolism of organisms. Citrate synthase, which is differentially expressed in the TCA cycle pathway, is a key enzyme in the sugar metabolism pathway and has a high specificity. It specifically catalyzes the condensation reaction of acetyl-CoA with oxaloacetate to form citric acid and coenzyme. Citrate synthase from Arabidopsis mitochondria [42], releasing more energy, could accelerate physiological metabolic reactions against adverse environments. This conclusion is the same as in our study. B. pervariabilis × D. grandis significantly increases the expressions of sucrose synthase and ATP-citrate synthase beta chain protein 1 after the inducer is introduced. The reason for this may be that in the starch and sucrose metabolism and TCA cycle pathways, B. pervariabilis × D. grandis uses sugar and other substances to oxidize and obtain energy. The various metabolites produced by the cyclic process are the raw materials for the biosynthesis of many important substances in the organism, which may be related to the plant immune response [43].
In the process of resisting stress, plants can regulate the activity of various antioxidant enzymes and reduce the damage of reactive oxygen species by regulating the changes of peroxidase activity in the metabolic pathway of phenylpropanoid biosynthesis [44]. Sugarcane enhances its ability to respond to low temperatures by increasing the activity of antioxidant enzymes in the body [45]. Studies on the cold resistance of rapeseed found that peroxidase isoenzyme activity was strong and stable in cold-resistant varieties [46]. The above conclusions are similar to the conclusions of this experiment. After the bamboo was induced by AP-toxin, the activity of peroxidase significantly changed. Perxidase (POX) is related to respiration, photosynthesis, and the oxidation of auxin. It removes excess reactive oxygen species from B. pervariabilis × D. grandis, maintains the active oxygen balance of the B. pervariabilis × D. grandis, and avoids damage to the B. pervariabilis × D. grandis cells by O2− [47].
The differentially enriched protein in the amino acid transport and metabolism pathway is phenylalanine ammonia-lyase. Phenylalanine ammonia-lyase catalyzes the deamination of phenylalanine to trans-cinnamic acid in the phenylpropane metabolic pathway, and then synthesizes coumaroyl-CoA through a series of catalytic reactions to finally synthesize polyphenolic compounds such as lignin, flavonoids, and flavonols [48]. PAL is a key enzyme in the specific synthesis pathway of lignin monomers and plays a crucial role in lignin biosynthesis [49]. Lignin [50] is not only the basic component of cell walls, but also participates in the process of cell wall lignification. The natural barrier formed by lignin and cellulose can effectively prevent the invasion of various pathogens [51] and participate in the signal of plant-pathogen interactions. The conduction process [52] is associated with plant immune responses [53]. Flavonoids and flavonols have also been considered as screening indicators for host disease resistance in previous studies [54]. Cinnamic acid can also form benzoic acid by oxidation and then hydroxylate to form salicylic acid. Studies have shown that PAL, as a key enzyme and rate-limiting enzyme in the phenylpropane synthesis pathway, is a branching point that controls primary metabolism to secondary metabolism and can respond to a variety of biotic and abiotic stresses, such as pathogen infection and biological factor induction [55]. In this study, phenylalanine ammonia-lyase was significantly up-regulated after the action of the B. pervariabilis × D. grandis inducer, which was consistent with previous studies.
Although the study found many important related proteins for B. pervariabilis × D. grandis response to AP-toxin, there is still a lack of understanding regarding its intrinsic mechanism of action. The function of the corresponding proteins and metabolic pathways can only be clarified by further combining physiological and biochemical research and verification of gene function. This would provide a clearer understanding of the mechanism of action for the corresponding proteins. Therefore, future work should use targeted proteomics technology to mine specific anti-reverse related proteins, verify their function and structure, and clone the corresponding genes. Such studies could use transgenic means to elucidate the expression regulation mechanism of the related genes/proteins and reveal the molecular regulatory network of B. pervariabilis × D. grandis to AP-toxin induction. Finally, such work would provide technical guidance for the breeding and genetic improvement of blight-resistant B. pervariabilis × D. grandis varieties, as well as better applications in production practice.

4. Materials and Methods

4.1. Materials

Plants: one-year-old B. pervariabilis × D. grandis planted in the cultivation area of Sichuan Agricultural University (103°01′N, 29°54′E, elevation 515.9 m, annual average temperature 16 °C, average annual precipitation 1400–1700 mm) in Sichuan, China.
Strain: Arthrinium phaeospermum (Corda) Ellis, isolated from B. pervariabilis × D. grandis blight, was provided by the Forest Protection Laboratory of Sichuan Agricultural University of Sichuan Province, China [11]. The pathogen suspension stored at −20 °C was activated on a PDA plate and cultured at 25 °C for 5 days. Dilute pathogen spores were added at a concentration of 10−7 to sterile water.
Purified AP-toxin was prepared by the Forest Conservation Laboratory of Sichuan Agricultural University [14,15,16] and configured as the inducer (diluted to 40 μg/mL with sterile distilled water and inactivated at 60 °C) [17].

4.2. Plant Sample Processing

Thirty plants of one-year-old B. pervariabilis × D. grandis with uniform growth were selected, and the pH of the induction was adjusted to 6–8. Fifteen plants were randomly selected and sprayed on the upper eight shoots by the inducer (inactivated AP-toxin) in 4.1. The remaining 15 plants were treated with sterile water, subjected to bagging moisturizing, and sprayed once every 12 h, 3 times in total. After three days, all of the B. pervariabilis × D. grandis plants underwent acupuncture inoculation with the A. phaeospermum suspension (the acupuncture site was at the tip of the twig, not pierced). Dropping the A. phaeospermum suspension at the wound was suitable to fill the shoots without flow, at about 0.05 mL. The plant underwent further bagging and moisturizing; and samples were collected after about 30–40 days. The samples were quickly stored in liquid nitrogen for future use.

4.3. Protein Extraction, Enzymatic Hydrolysis and Peptide Quantification

We took a 300 μg aliquot of each sample, froze them with liquid nitrogen, and ground them into powder. We took an appropriate amount of powder from each sample, added 200 μL of lysate, sonicated them, and subjected them to TCA-acetone precipitation at −20 °C overnight. A 16,000 g centrifugation was conducted for fifteen minutes at −4 °C, the precipitate was washed twice with cold acetone and dried in air. Then, 150 μL of lysate was added to each tube sample, sonicated, centrifuged at 16,000 g for 15 min, and the supernatant was taken. Protein quantification was performed using the bicinchoninic acid (BCA) method (Beyotime Institute of Biotechnology, Shanghai, China) [56].
We took 300 μg aliquot of each sample for enzymatic hydrolysis [57], added DL-Dithiothreitol (DTT) to 100 mM of the sample, boiled in water for 5 min, and cooled to room temperature. We added 200 μL of Urea (UA) buffer (8 M urea, 150 mM Tris-HCl, pH 8.0), mixed, transferred it to a 10 KD ultrafiltration centrifuge tube, and centrifuged the sample at 12,000 g for 15 min. We added 200 μL of UA buffer and centrifuged it at 12,000 g for 15 min, then discarded the filtrate. We added 100 μL of IAA (50 mM IAA in UA), shook it at 600 rpm for 1 min, avoiding light, left it at room temperature for 30 min, and centrifuged at 12,000 g for 10 min. We added 100 μL of UA buffer, centrifuged at 12,000 g for 10 min, repeating these two times. We added 100 μL of NH4HCO3 buffer (Sigma-Aldrich, St. Louis, USA), centrifuged at 14,000 g for 10 min, repeating these two times. We added 40 μL trypsin buffer (6 μg Trypsin in 40 μL NH4HCO3 buffer), shook it at 600 rpm for 1 min, and then kept it at 37 °C for 16–18 h. We replaced the collection tube, centrifuged at 12,000 g for 10 min, collected the filtrate, added an appropriate amount of 0.1% trifluoroacetic acid solution, desalinated the C18 Cartridge, quantified the peptide, and determined the peptide concentration by OD280.

4.4. TMT Peptide Labeling and Peptide Grading

One hundred milligrams of peptide were taken from each sample and labeled according to the Thermo Fisher TMT Labeling Kit instructions. The experimental induction group and the control group were labeled as TMT-126 and TMT-127, and three replicates were set. The labeled peptides were mixed in equal amounts, and the dried peptides were fractionated using High-pH (Pierce™ High pH Reversed-Phase Peptide Fractionation Kit, Thermo Fisher Scientific, Waltham, USA). The sample collection was eventually combined into ten components. The peptides of each component were dried and reconstituted with 0.1% FA for LC-MS analysis. TMT markers and LC-MS data acquisition were determined by Shanghai Baipu Biotechnology Co., Ltd. (Shanghai, China).

4.5. Liquid Chromatography Mass Spectrometry (LC-MS/MS) Analysis

The reconstituted peptide solution was subjected to LC-MS/MS analysis, and each fraction was injected once for a total of 10 mass spectrometry analyses. Separation was performed using a nanoscale flow rate HPLC liquid phase system Easy-nLC1200 (Thermo Fisher Scientific, Waltham, USA). Buffer: solution A was 0.1% formic acid in water, and solution B was 0.1% formic acid in acetonitrile (98% acetonitrile). The column was equilibrated with 95% A solution. The sample was loaded onto a Trap column (2 cm × 100 μm, 5 μm: C18) and separated by an analytical Thermo scientific EASY column (75 μm × 120 mm, 3 μm: C18) at a flow rate of 300 nL/min. The relevant liquid phase gradient was as follows: 0–2 min, with a linear gradient of B liquid from 4–7%; 2–57 min, with a linear gradient of B liquid from 7–30%; 57–62 min, with a linear gradient of B liquid from 30–45%; 62–67 min, with a linear gradient of B liquid from 45–90%; 67–75 min, with B liquid maintained at 90%. The peptides were chromatographed and analyzed by tandem mass spectrometry using a Q-Exactive Plus Mass Spectrometer (Thermo Scientific). Analysis duration: 75 min. Detection method: positive ion. Parent ion scanning range: 300–1800 m/z. The mass-to-charge ratio of peptide and polypeptide fragments was collected as follows: 20 fragment maps (MS2 scan, higher-energy collisional dissociation (HCD)) were acquired after each full scan. Primary mass spectrometry resolution: 70,000 @ m/z 200, AGC target: 1e6, Level 1 Maximum IT: 50 ms. Secondary MS resolution: 17,500 @ m/z 200. AGC target: 1e5. Level 2 Maximum IT: 50 ms, MS2 Activation Type: HCD. Isolation window: 1.6 Th. Normalized collision energy: 35.

4.6. Mass Spectrometry Data Analysis

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the iProX partner repository [58] with the dataset identifier PXD014363, the subject ID is IPX0001646002.
The original RAW file obtained by LC-MS/MS was imported into MaxQuant software (version 1.0.6.16) for database searching. In this experiment, the protein database translated by the B. pervariabilis × D. grandis transcriptome (All-gene.fa) file was selected for retrieval and analysis. The main library parameters were set as in Table 4.

4.7. Bioinformatics Analysis

Bioinformatics analysis of the B. pervariabilis × D. grandis induction group and control group was performed. The functional annotation of proteins was performed using the Blast2 GO database program (http://geneontology.org/), while the KEGG database (http://www.genome.jp/kegg/) was used to classify and group identified proteins. The Fisher’s exact test was used to analyze the significance level of each pathway and certain GO term protein enrichments, and the target proteins were analyzed by the GO [18] and KEGG [59] pathways. We searched for the protein genes in the IntAct molecular interaction database to obtain the mutual interaction data between proteins (PPI). The results were downloaded and imported into the Cytoscape package (version 3.2.1) in extensible graph markup and modelling language (XGMML) format to visualize and further analyze the function of the differentially expressed proteins. In addition, the importance of the proteins in the PPI network was assessed by calculating the degree of interaction for each protein.

4.8. PRM Verification

In order to verify the accuracy of TMT and LC-MS/MS data analysis, important target proteins were selected for targeted shotgun mass spectrometry, in order to directly compare the abundance of the differentially expressed proteins. We imported the PRM analysis protein information into the software Xcalibur for PRM method settings. The PRM mass spectrometry RAW file of the hybrid bamboo induction group and the experimental group was opened using the software Xcalibur 4.0.
We took 2 μg of peptide from each sample for LC-PRM/MS analysis, and performed chromatographic separation using a nanoliter flow rate Easy nLC 1200 chromatography system (Thermo Scientific) after loading. Buffer: solution A was 0.1% aqueous formic acid, and solution B was 0.1% formic acid, acetonitrile, and water (95% acetonitrile). The column was equilibrated with 95% solution A. The sample was injected into a trap column (100 μm × 20 mm, 5 μm, C18, Dr. Maisch GmbH, Ammerbuch-Entringen, Germany) and subjected to gradient separation through a chromatography column (75 μm × 150 mm, 3 μm, C18, Dr. Maisch GmbH) at a flow rate of 300 nL/min. The liquid phase separation gradient was as follows: 0–5 min, with a linear gradient of liquid B from 2% to 5%; 5–45 min, with a linear gradient of liquid B from 5% to 23%; 45–50 min, with a linear gradient of liquid B from 23% to 40%; 50–52 min, with a linear gradient of liquid B from 40% to 100%; 52–60 min, with liquid B maintained at 100%. Peptide fragmentation and targeted PRM mass spectrometry were performed using a Q-Exactive Plus Mass Spectrometer (Thermo Scientific). The analysis time was 75 min. Detection mode: positive ion. Parent ion scanning range: 350–1500 m/z. Primary mass spectrometry resolution: 70,000 @ m/z 200. AGC target: 3e6. First-order mass spectrometer maximum IT: 200 ms. Peptide secondary mass spectrometry was performed as follows: after each full scan (Full MS scan), the precursor m/z of 18 target peptides was sequentially selected according to the inclusion list for secondary mass spectrometry analysis (MS2), in which MS2 was resolved. The rate was 17,500 @ m/z 200. AGC target: 3e6, secondary mass spectrometry maximum IT: 100 ms. MS2 activation type: HCD. Isolation window: 2.0 Th. Normalized collision energy: 27. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the iProX partner repository [58] with the dataset identifier PXD014364, the subject ID is IPX0001651001. The resulting mass spectrum RAW file was analyzed using the software Skyline 4.1 [60] for PRM data.

4.9. Data Analysis Software

The data of the experimental induction group and the control group were collected three times, and the data were statistically analyzed using SPSS v20.0. The significance of differences was determined by Student’s t test, and the significance level was p < 0.05. All data shown are the means ± standard deviations (SDs) (n = 3).

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-1989/9/8/166/s1. Table S1: Gene Ontology (GO) functional analysis of proteins of B. pervariabilis × D. grandis; Table S2: Kyoto Encyclopedia of Genes and Genomes (KEGG) functional analysis of proteins of B. pervariabilis × D. grandis; Table S3: Skyline analysis of candidate peptide fragments of target proteins in different samples; Table S4: Analysis of target peptide PRM quantitative Skyline date-Transition results; Table S5: Analysis of target peptide PRM quantitative Skyline date-Normalized AUC; Table S6: Analysis of target peptide PRM quantitative Skyline data-Peptide quantitative; Figure S1: Secondary spectrum of candidate peptides.

Author Contributions

Conceptualization, S.L. and T.Z.; Methodology, T.Z.; Software, S.L.; Validation, Q.H.; Formal analysis, Q.H.; Investigation, T.Z.; Resources, X.F.; Data curation, Q.H. and H.Z.; Writing—original draft, Q.H. and X.F.; Writing—review and editing, S.L. and H.Z.; Visualization, S.H.; Supervision, T.Z.; Project administration, S.H.; Funding acquisition, S.L.

Funding

This work was supported by the National Natural Science Foundation of China [grant number 31700568]; the Certificate of China Postdoctoral Science Foundation [grant number 2018T110989]; the China Postdoctoral Science Foundation [grant number 2016M602705]; and the College Students’ Innovation and Entrepreneurship Training Program of Sichuan Province [grant number 201810626117].

Acknowledgments

We are grateful to MDPI English editing who provided language help.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bronwyn, J.B.; Thelma, C.C.; Jose, L.D.; Andrea, M.; Hans, P.M.; Francisco, P.A.; Ghasem, H.S.; Katja, W.; Christian, Z. Elucidation of salt stress defense and tolerance mechanisms of crop plants using proteomics-current achievements and perspectives. Proteomics 2013, 13, 1885–1900. [Google Scholar]
  2. Wang, Q.G.; Li, Z.; Liu, W.; Pan, J. Studies on drought-responsive proteins of foxtail millet by applying quantitative proteomic. Shandong Agric. Sci. 2018, 50, 1–8. [Google Scholar]
  3. Zhong, Y. Research of Transcription Genomics and Proteomics of Citrus Induced by Candidatus Liberibacter Asiaticus. Ph.D. Thesis, Hunan Agricultural University, Changsha, China, 2012. [Google Scholar]
  4. Chen, Q.Z.; Guo, W.S.; Feng, L.Z.; Miao, S.H.; Lin, J.S. Pathogenic Identification of Cylindro Cladium Leaf Blight Collected from Fujian Province and Studies on Transcriptomics and Proteomics of Eucalyptus Induced by Calonectria Pseudoreteaudii. Ph.D. Thesis, Fujian Agriculture and Forestry University, Fuzhou, China, 2013. [Google Scholar]
  5. González-Fernández, R.; Prats, E.; Jorrín-Novo, V.J. Proteomics of plant pathogenic fungi. J. Biomed. Biotechnol. 2010, 2010, 932527. [Google Scholar] [CrossRef] [PubMed]
  6. Rustagi, A.; Singh, G.; Agrawal, S.; Gupta, P.K. Proteomic studies revealing enigma of plant–pathogen interaction. In Molecular Aspects of Plant-Pathogen Interaction; Springer: Berlin, Germany, 2018; pp. 239–264. [Google Scholar]
  7. Yang, Y.L.; Xiao, L.T.; Hu, X.Q. Study on the relationship between the toxin of Phytophthora infestans (Mont.) de bary and resistance of potato. Agric. Sci. China 2011, 10, 238–245. [Google Scholar] [CrossRef]
  8. Hui, Q.; Liu, H.P.; Han, J.C. Induce of crude toxin from Cladosporium fulvum on defense enzyme and activity oxygen of tomato seedlings. North China Agric. J. 2009, 24, 215–218. [Google Scholar]
  9. Ma, C.H.; Zhai, C.X.; Wang, L.A.; Chen, X.; Li, Y.C.; Guo, X.L.; Cui, S.P.; Li, G.M. Induced resistance by the toxin filtrate of Bipolaris maydis race T cultivation. Agric. Sci. China 2006, 5, 678–684. [Google Scholar] [CrossRef]
  10. Conrath, U.; Thulke, O.; Katz, V.; Schwindling, S.; Kohler, A. Priming as a mechanism in induced systemic resistance of plant. Eur. J. Plant Pathol. 2001, 107, 113–119. [Google Scholar] [CrossRef]
  11. Zhu, T.H.; Huang, Z.C.; Gao, Q.Z.; Mou, Y.M.; Luo, Y.Q.; Mao, X.T.; Li, X.D. Research on survey and space distribution of Bambusa pervariabilis × Dendrocalamopsis daii blight. J. Sichuan For. Sci. Technol. 2009, 30, 29–31, 37. [Google Scholar]
  12. Zhu, J.H. Preliminarily approaching the relationship between the universal incidence and gravity of Ceratosphaeria phyllostachydis. J. Fujian For. Sci. Tech. 1997, 24, 43–46. [Google Scholar]
  13. Ma, G.L.; Hu, G.B.; Yu, C.Z.; Wu, J.L.; Xu, B.C. Phyllostachys prominens plum shoot wilt pathogenic fungoid and its biological characteristics. J. Zhejiang For. Coll. 2003, 20, 44–48. [Google Scholar]
  14. Li, S.J. Protein toxin produced from Arthrinium Phaeospermum Causing Blight in Bambusa Pervariabilis × Dendrocalamopsis Grandis and Its Accurate Pathogenic Mechanism. Ph.D. Thesis, Sichuan Agricultural University, Chengdu, China, July 2013. [Google Scholar]
  15. Li, S.J.; Zhu, T.H. Purification of the toxin protein Pc from Arthrinium phaeospermum and its effect on the defence enzymes of Bambusa pervariabilis × Dendrocalamopsis grandis varieties. For. Pathol. 2014, 44, 96–106. [Google Scholar] [CrossRef]
  16. Li, S.J.; Zhu, T.H.; Zhu, H.M.Y.; Liang, M.; Qiao, T.M.; Han, S.; Che, G.N. Purification of protein AP-toxin from Arthrinium phaeospermum causing blight in Bambusa pervariabilis × Dendrocalamopsis grandis and its metabolic effects on four bamboo varieties. Phytopathology 2013, 103, 135–145. [Google Scholar] [CrossRef]
  17. He, Q.Q.; Liu, Y.X.; Fang, X.M.; Zhu, T.H.; Qiao, T.M.; Han, S.; Li, S.J. Screening of induced factor of Bambusa pervariabilis × Dendrocalamopisis grandis against Arthrinium phaeospermum and its duration. Plant Prot. 2019, 45, 81–89. [Google Scholar]
  18. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene ontology: Tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef]
  19. Peterson, A.C.; Russell, J.D.; Bailey, D.J.; Westphall, M.S.; Coon J, J. Parallel reaction monitoring for high resolution and high mass accuracy quantitative, targeted proteomics. Mol. Cell. Proteom. 2012, 11, 1475–1488. [Google Scholar] [CrossRef]
  20. Overall, C.M. Can proteomics fill the gap between genomics and phenotypes? J. Proteom. 2014, 100, 1–2. [Google Scholar] [CrossRef]
  21. Shoresh, M.; Harman, G.E.; Mastouri, F. Induced systemic resistance and plant responses to fungal biocontrol agents. Annu. Rev. Phytopathol. 2010, 48, 21–43. [Google Scholar] [CrossRef]
  22. Yu, Z.M.; Li, J.Y.; Lin, Z.H.; Zhang, Q.; He, H.B. Research advances on induced-defense of plant resistance. J. Agric. Sci. 2013, 34, 69–76. [Google Scholar]
  23. Lan, P.; Li, W.; Wen, T.N.; Shiay, J.Y.; Wu, Y.C.; Lin, W.; Schmidt, W. iTRAQ protein profile analysis of Arabidopsis roots reveals new aspects critical for iron homeostasis. Plant Physiol. 2011, 155, 821–834. [Google Scholar] [CrossRef]
  24. Fan, J.; Chen, C.; Yu, Q.; Brlansky, R.H.; Li, Z.J.; Gmitter, G.G.J. Comparative iTRAQ proteome and transcriptome analyses of sweet orange infected by “Candidatus liberibacter asiaticus”. Physiol. Plantarum 2011, 143, 235–245. [Google Scholar] [CrossRef]
  25. Lippert, D.N.; Ralph, S.G.; Phillips, M.; Phillips, M.; White, R.; Smith, D.; Hardie, D.; Gershenzon, J.; Ritland, K.; Borchers, C.H.; et al. Quantitative iTRAQ proteome and comparative transcriptome analysis of elicitor-induced Norway spruce (Picea abies) cells reveals elements of calcium signaling in the early conifer defense response. Proteomics 2009, 9, 350–367. [Google Scholar] [CrossRef]
  26. Duffey, S.; Stout, M. Antinutritive and toxic components of plant defense against insects. Arch. Insect Biochem. Physiol. 1996, 32, 3–37. [Google Scholar] [CrossRef]
  27. Felton, G.W. Nutritive quality of plant protein; sources of variation and insect herbivore responses. Arch. Insect Biochem. Physiol. 1996, 32, 107–130. [Google Scholar] [CrossRef]
  28. Wang, J.; Constabel, C. Polyphenol oxidase overexpression in transgenic Populus enhances resistance to herbivory by forest tent caterpillar (Malacosoma disstria). Planta 2004, 220, 87–96. [Google Scholar] [CrossRef]
  29. Henrik, U.S.; Barry, R.P.; Jürgen, K.; Kerstin, W.; Jacqueline, F.; Antje, B.; Thomas, M.O. Induced plant defense responses against chewing insects. Ethylene signaling reduces resistance of Arabidopsis against Egyptian cotton worm but not diamondback moth. Plant Physiol. 2000, 124, 1007–1017. [Google Scholar]
  30. Timpa, J.D.; Burke, J.J.; Quiseberry, J.E.; Wendt, C.W. Effects of water stress on the organic acid and carbohydrate compositions of cotton plants. Plant Physiol. 1986, 3, 724–728. [Google Scholar] [CrossRef]
  31. Francoise, F.; Daniel, L.R.; John, G.S. Effects of salt stress on amino acid, organic acid, and carbohydrate composition of roots, bacteroids and cytosol of alfalfa (Medicago sativa L.). Plant Physiol. 1991, 96, 1228–1236. [Google Scholar]
  32. Ghosh, D.; Xu, J. Abiotic stress responses in plant roots: A proteomics perspective. Front. Plant Sci. 2014, 5, 6. [Google Scholar] [CrossRef]
  33. Wang, X.; Fan, P.; Song, H.; Chen, X.; Li, X.; Li, Y. Comparative proteomic analysis of differentially expressed proteins in shoots of Salicornia europaea under different salinity. J. Proteome Res. 2009, 8, 3331–3345. [Google Scholar] [CrossRef]
  34. Zhao, Q.; Zhang, H.; Wang, T.; Chen, S.X.; Dai, S.J. Proteomics-based investigation of salt-responsive mechanisms in plant roots. J. Proteomics 2013, 82, 230–253. [Google Scholar] [CrossRef]
  35. Waduwara-Jayabahu, I.; Oppermann, Y.; Wirtz, M.; Hull, Z.T.; Schoor, S.; Plotnikov, A.N.; Hell, R.; Sauter, M.; Moffatt, B.A. Recycling of methylthioadenosine is essential for normal vascular development and reproduction in Arabidopsis. Plant Physiol. 2012, 158, 1728–1744. [Google Scholar] [CrossRef]
  36. Zheng, B.B.; Fang, Y.N.; Pan, Z.Y.; Sun, L.; Deng, X.X.; Grosser, J.W.; Guo, W.W. iTRAQ-based quantitative proteomics analysis revealed alterations of carbohydrate metabolism pathways and mitochondrial proteins in a male sterile cybrid pummelo. J. Proteome Res. 2014, 13, 2998–3015. [Google Scholar] [CrossRef]
  37. Coleman, H.D.; Yan, J.; Mansfield, S.D. Sucrose synthase affects carbon partitioning to increase cellulose production and altered cell wall ultrastructure. Proc. Natl. Acad. Sci. USA 2009, 106, 13118–13123. [Google Scholar] [CrossRef] [Green Version]
  38. Ruan, Y.L. Sucrose metabolism: Gateway to diverse car-bon use and sugar signaling. Annu. Rev. Plant. Biol. 2014, 65, 33–67. [Google Scholar] [CrossRef]
  39. Xu, S.M.; Brill, E.; Llewellyn, D.J.; Llewellyn, D.J.; Furbank, R.T.; Ruan, Y.L. Overexpression of a potato sucrose synthase gene in cotton accelerates leaf expansion, reduces seed abortion, and enhances fiber production. Mol. Plant 2012, 5, 430–441. [Google Scholar] [CrossRef]
  40. Shewmaker, C.K.; Boyer, C.D.; Wiesenborn, D.P.; Thompson, D.B.; Boersig, M.R.; Oakes, J.V.; Stalker, D.M. Expression of Escherichia coli glycogen synthase in the tubers of transgenic potatoes (Solanum tuberosum) results in a highly branched starch. Plant Physiol. 1994, 104, 1159–1166. [Google Scholar] [CrossRef]
  41. Li, Y.Q.; Zhao, S.B.; Ni, K.C.; Chen, W.R.; Xin, D.D.; Guo, W.D. Studies on differentially expressed genes in flower bud during dormancy release in Prunus pseudocerasus Lindl. ‘duanbing’. Acta Hortic. Sin. 2012, 39, 1253–1262. [Google Scholar]
  42. Schmidtmann, E.; Konig, A.C.; Orwat, A.; Leister, D.; Hartl, M.; Finkemeier, I. Redox regulation of Arabidopsis mitochondrial citrate synthase. Mol. Plant 2014, 7, 156–169. [Google Scholar] [CrossRef]
  43. Yang, J.L. Effect of l-Glutamate on Inhibiting Postharvest Diseases by Inducing Resistance in Fruit and the Possible Defense Mechanisms Involved. Ph.D. Thesis, Zhejiang Univercity, Hangzhou, China, 2017. [Google Scholar]
  44. Mittler, R.; Vanderauwera, S.; Gollery, M.; Van Breusegem, F. Reactive oxygen gene network of plants. Trends. Plant Sci. 2004, 9, 490–498. [Google Scholar] [CrossRef]
  45. Sun, F. Effects of Chilling Stress on Physiological Metabolism Protein Expression in Chloroplasts from Settlings of Sugarcane. Ph.D. Thesis, Guangxi university, Nanning, China, 2011. [Google Scholar]
  46. Song, L.T. Study on the relations of esterase, peroxidase and cold hardiness of winter rape (B. napu L.). J. Anhui Agric. Sci. 1995, 23, 333–336. [Google Scholar]
  47. Ma, X.J.; Zhu, D.H. Functional roles of the plant superoxide dismutase. Hereditas 2003, 25, 225–231. [Google Scholar]
  48. Avdiushko, S.A.; Ye, X.; Kuc, J. Detection of several enzymatic activities in leaf prints of cucumber plants. Physiol. Mol. Plant P 1993, 42, 441–454. [Google Scholar] [CrossRef]
  49. Zhu, J.X.; Sun, J.J.; Yuan, X.L.; Wang, J.; Yang, Y.M.; Wang, Y. Cloning and Sequence Analysis of Cinnamoyl Alcohol Dehydrogenase Gene from Guadua paraguayana. J. West China For. Sci. 2017, 1, 6–11. [Google Scholar]
  50. Whetten, R.W.; Mackay, J.J.; Sederoff, R.R. Recent advances in understanding lignin biosynthesis. Annu. Rev. Plant Physiol. Plant Mol Biol. 1998, 1, 585–609. [Google Scholar] [CrossRef]
  51. Frey, M.; Chomet, P.; Glawischnig, E.; Stettner, C.; Grun, S.; Winklmair, A.; Eisenreich, W.; Bacher, A.; Meeley, R.B.; Briggs, S.P.; et al. Analysis of a chemical plant defense mechanism in grasses. Science 1997, 277, 696–699. [Google Scholar] [CrossRef]
  52. Yoda, H.; Fujimura, K.; Takahashi, H.; Munemura, I.; Uchimiya, H.; Sano, H. Polyamines as a common source of hydrogen peroxide in host-and nonhost hypersensitive response during pathogen infection. Plant Mol. Biol. 2009, 70, 103–112. [Google Scholar] [CrossRef]
  53. Bell, A.A. Biochemical mechanism of disease resistance. Ann. Rev. Plant Physiol. 1981, 32, 21–81. [Google Scholar] [CrossRef]
  54. Qiu, Y.X.; Ke, Y.Q.; Dai, H.J.; Pan, T.G. Phenolics metabolism in sweet potato infected by Fusarium oxysporum Schlecht. f. sp. Chin. J. Eco Agric. 2007, 15, 167–170. [Google Scholar]
  55. Ritter, H.; Schulz, G.E. Structural basis for the entrance into the phenylpropanoid metabolism catalyzed by phenylalanine ammonia-lyase. Plant Cell 2004, 16, 3426–3436. [Google Scholar] [CrossRef]
  56. Smith, P.K.; Krohn, R.I.; Hermanson, G.T.; Mallia, A.K.; Gartner, F.H.; Provenzano, M.D.; Fujimoto, E.K.; Goeke, N.M.; Olson, B.J.; Klenk, D.C. Measurement of protein using bicinchoninic acid. Anal. Biochem. 1985, 150, 76–85. [Google Scholar] [CrossRef]
  57. Wiśniewski, J.R.; Zougman, A.; Nagaraj, N.; Mann, M. Universal sample preparation method for proteome analysis. Nat. Methods 2009, 6, 359–362. [Google Scholar] [CrossRef]
  58. Ma, J.; Chen, T.; Wu, S.F.; Yang, C.Y.; Bai, M.Z.; Shu, K.X.; Li, K.L.; Zhang, G.Q.; Jin, Z.; He, F.C.; et al. iProX: An integrated proteome resource. Nucleic Acids Res. 2018, 47, D1211–D1217. [Google Scholar] [CrossRef]
  59. Kanehisa, M.; Goto, S.; Sato, Y.; Furumichi, M.; Tanabe, M. KEGG for integration and interpretation of largescale molecular data sets. Nucleic Acids Res. 2012, 40, D109–D114. [Google Scholar] [CrossRef]
  60. MacLean, B.; Tomazela, D.M.; Shulman, N.; Chambers, M.; Finney, G.L.; Frewen, B.; Kern, R.; Tabb, D.L.; Liebler, D.C.; MacCoss, M.J. Skyline: An open source document editor for creating and analyzing targeted proteomics experiments. Bioinformatics 2010, 26, 966–968. [Google Scholar] [CrossRef]
Metabolites 09 00166 g001 550
Figure 1. Photos of characteristic samples of B. pervariabilis × D. grandis after A. phaeospermum (AP) infection. (A) Treated with AP-toxin, and (B) treated with water before infection with A. phaeospermum.
Figure 1. Photos of characteristic samples of B. pervariabilis × D. grandis after A. phaeospermum (AP) infection. (A) Treated with AP-toxin, and (B) treated with water before infection with A. phaeospermum.
Metabolites 09 00166 g001
Metabolites 09 00166 g002a 550Metabolites 09 00166 g002b 550
Figure 2. Gene ontology (GO) analysis results for significant protein enrichment between the B. pervariabilis × D. grandis induced by AP-toxin and the control group. (A) The top 10 most significant enriched entries in the three categories biological processes (BP), cell component (CC), and molecular function (MF). The p value is set to 0.05, and the items in the same category are sorted by the p value. The left y-axis represents the percentage, and the right y-axis represents the number of genes/proteins enriched under a certain category. (B) The results for the biological processes that are significantly enriched, (C) the results for the significantly enriched cellular components, and (D) the results for the molecular functions that are significantly enriched. In (B,D), the horizontal axis indicates the percentage of genes/proteins enriched in each case, and each bar is followed by the number of proteins and the corresponding p value.
Figure 2. Gene ontology (GO) analysis results for significant protein enrichment between the B. pervariabilis × D. grandis induced by AP-toxin and the control group. (A) The top 10 most significant enriched entries in the three categories biological processes (BP), cell component (CC), and molecular function (MF). The p value is set to 0.05, and the items in the same category are sorted by the p value. The left y-axis represents the percentage, and the right y-axis represents the number of genes/proteins enriched under a certain category. (B) The results for the biological processes that are significantly enriched, (C) the results for the significantly enriched cellular components, and (D) the results for the molecular functions that are significantly enriched. In (B,D), the horizontal axis indicates the percentage of genes/proteins enriched in each case, and each bar is followed by the number of proteins and the corresponding p value.
Metabolites 09 00166 g002aMetabolites 09 00166 g002b
Metabolites 09 00166 g003a 550Metabolites 09 00166 g003b 550
Figure 3. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways with significant protein enrichment between the B. pervariabilis × D. grandis induced by AP-toxin and the control group. (A) The top 10 KEGG signaling pathways with the most significant p values. Two limits of p values, 0.05 and 0.01, are indicated. (B) KEGG pathway classification indicates the signal pathways that are significantly enriched. Each bar is followed by the number of proteins and the corresponding p value. The KEGG pathway is divided into the following categories: A: metabolism, B: genetic information processing, C: environmental information processing, D: cellular processes, E: organismal systems, F: human diseases.
Figure 3. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways with significant protein enrichment between the B. pervariabilis × D. grandis induced by AP-toxin and the control group. (A) The top 10 KEGG signaling pathways with the most significant p values. Two limits of p values, 0.05 and 0.01, are indicated. (B) KEGG pathway classification indicates the signal pathways that are significantly enriched. Each bar is followed by the number of proteins and the corresponding p value. The KEGG pathway is divided into the following categories: A: metabolism, B: genetic information processing, C: environmental information processing, D: cellular processes, E: organismal systems, F: human diseases.
Metabolites 09 00166 g003aMetabolites 09 00166 g003b
Metabolites 09 00166 g004 550
Figure 4. Protein-protein interaction (PPI) network map of the differentially expressed proteins in the B. pervariabilis × D. grandis between the induced group and the control group.
Figure 4. Protein-protein interaction (PPI) network map of the differentially expressed proteins in the B. pervariabilis × D. grandis between the induced group and the control group.
Metabolites 09 00166 g004
Table
Table 1. Statistics on the protein identification results.
Table 1. Statistics on the protein identification results.
Identification ResultUnique PeptidesQuantified ProteinsUp-RegulatedDown-RegulatedSignificantly Different Proteins
Total332017917588163
Table
Table 2. Differentially expressed protein information for high connectivity in the PPI analysis.
Table 2. Differentially expressed protein information for high connectivity in the PPI analysis.
Protein IDDescriptionKEGG Pathways
PH01000123G1100porphobilinogen deaminaseBiosynthesis of secondary metabolites
Metabolic pathways
PH01000159G0130sucrose synthaseStarch and sucrose metabolism
PH01000713G0340Adenosine triphosphate (ATP)-citrate synthase beta chain protein 1Citrate cycle (Tricarboxylic acid (TCA) cycle)
Energy production and conversion
PH01000761G0570PeroxidasePhenylpropanoid biosynthesis
Carbohydrate transport and metabolism
PH01003309G0170Phenylalanine ammonia-lyasePhenylalanine metabolism
Phenylpropanoid biosynthesis
Secondary metabolites biosynthesis, transport and catabolism
Table
Table 3. Quantitative results for 14 candidate proteins determined using the PRM and tandem mass tag (TMT) methods.
Table 3. Quantitative results for 14 candidate proteins determined using the PRM and tandem mass tag (TMT) methods.
Protein NamePeptide SequenceAP-Toxin + A. phaeospermum Sterile Water + A. phaeospermumPRM Fold ChangeTMT Fold ChangeConsistency between PRM and TMT
PH01000087G1730AVAHQPVSVAIEAGGR1,756,917.4221,409,646.381.251.52Yes
PH01000123G1100TTGDMILDKPLADIGGK3,128,973.5763,315,437.0470.941.62No
PH01000133G0340LVENDEVVR6,132,972.8984710,029.371.302.01Yes
PH01000159G0130QQGLNITPR30,652,687.985,578,763.1545.493.27Yes
PH01000713G0340FGGAIDDAAR21,475,599.318,257,753.8412.601.63Yes
PH01000761G0570NNPSDIDPSLNPSYAK4,349,530.1581,633,319.5232.661.59Yes
PH01000860G0570DVDLSTYK2,738,630.0871,653,255.4041.661.54Yes
PH01000898G0600MGNINPLTGTAGQIR66,407,505.35444,831,682.90.150.35Yes
PH01001064G0040EHLIAGWAPK7,994,610.9084,179,165.6951.911.44Yes
PH01001724G0160VNVYYNEASCGR4,996,071.9073,097,515.3461.611.57Yes
PH01001918G0090APDFEAEAVFDQEFIK1,768,290.2643,727,484.20.470.64Yes
PH01002295G0170VVVSSCGHDGPFGATGVK11,037,762.066,282,695.561.761.45Yes
PH01003309G0170VGQVAAVAQAK63,790,370.9819,648,147.323.252.24Yes
PH01100083G0010YFSAAASQALDTAER2,771,584.72316,486,693.420.170.46Yes
Table
Table 4. MaxQuant search library parameter settings.
Table 4. MaxQuant search library parameter settings.
ItemValue
TypeReporter ion MS2
Isobaric labelsTMT 6plex
EnzymeTrypsin
Reporter mass tolerance0.005 Da
Max missed cleavages2
Main search peptide tolerance4.5 ppm
First search peptide tolerance20 ppm
Mass spectrometry (MS) /MS tolerance20 ppm
Fixed modificationsCarbamidomethyl (C)
Variable modificationsOxidation (M), Acetyl (Protein N-term)
DatabaseBamboo.fasta
Database patternTarget-Reverse
Peptide spectrum matches false discovery rate (PSM FDR)≤0.01
Protein FDR≤0.01
Protein quantificationRazor and unique peptides were used for protein quantification.

Share and Cite

MDPI and ACS Style

He, Q.; Fang, X.; Zhu, T.; Han, S.; Zhu, H.; Li, S. Differential Proteomics Based on TMT and PRM Reveal the Resistance Response of Bambusa pervariabilis × Dendrocalamopisis grandis Induced by AP-Toxin. Metabolites 2019, 9, 166. https://doi.org/10.3390/metabo9080166

AMA Style

He Q, Fang X, Zhu T, Han S, Zhu H, Li S. Differential Proteomics Based on TMT and PRM Reveal the Resistance Response of Bambusa pervariabilis × Dendrocalamopisis grandis Induced by AP-Toxin. Metabolites. 2019; 9(8):166. https://doi.org/10.3390/metabo9080166

Chicago/Turabian Style

He, Qianqian, Xinmei Fang, Tianhui Zhu, Shan Han, Hanmingyue Zhu, and Shujiang Li. 2019. "Differential Proteomics Based on TMT and PRM Reveal the Resistance Response of Bambusa pervariabilis × Dendrocalamopisis grandis Induced by AP-Toxin" Metabolites 9, no. 8: 166. https://doi.org/10.3390/metabo9080166

APA Style

He, Q., Fang, X., Zhu, T., Han, S., Zhu, H., & Li, S. (2019). Differential Proteomics Based on TMT and PRM Reveal the Resistance Response of Bambusa pervariabilis × Dendrocalamopisis grandis Induced by AP-Toxin. Metabolites, 9(8), 166. https://doi.org/10.3390/metabo9080166

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

Article Metrics

Back to TopTop