Figures
Abstract
Wheat (Triticum aestivum L.) is an important crop worldwide. The physiological deterioration of seeds during storage and seed priming is closely associated with germination, and thus contributes to plant growth and subsequent grain yields. In this study, wheat seeds during different stages of artificial ageing (45°C; 50% relative humidity; 98%, 50%, 20%, and 1% Germination rates) and priming (hydro-priming treatment) were subjected to proteomics analysis through a proteomic approach based on the isobaric tandem mass tag labeling. A total of 162 differentially expressed proteins (DEPs) mainly involved in metabolism, energy supply, and defense/stress responses, were identified during artificial ageing and thus validated previous physiological and biochemical studies. These DEPs indicated that the inability to protect against ageing leads to the incremental decomposition of the stored substance, impairment of metabolism and energy supply, and ultimately resulted in seed deterioration. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the up-regulated proteins involved in seed ageing were mainly enriched in ribosome, whereas the down-regulated proteins were mainly accumulated in energy supply (starch and sucrose metabolism) and stress defense (ascorbate and aldarate metabolism). Proteins, including hemoglobin 1, oleosin, agglutinin, and non-specific lipid-transfer proteins, were first identified in aged seeds and might be regarded as new markers of seed deterioration. Of the identified proteins, 531 DEPs were recognized during seed priming compared with unprimed seeds. In contrast to the up-regulated DEPs in seed ageing, several up-regulated DEPs in priming were involved in energy supply (tricarboxylic acid cycle, glycolysis, and fatty acid oxidation), anabolism (amino acids, and fatty acid synthesis), and cell growth/division. KEGG and protein-protein interaction analysis indicated that the up-regulated proteins in seed priming were mainly enriched in amino acid synthesis, stress defense (plant-pathogen interactions, and ascorbate and aldarate metabolism), and energy supply (oxidative phosphorylation and carbon metabolism). Therefore, DEPs associated with seed ageing and priming can be used to characterize seed vigor and optimize germination enhancement treatments. This work reveals new proteomic insights into protein changes that occur during seed deterioration and priming.
Citation: Lv Y, Zhang S, Wang J, Hu Y (2016) Quantitative Proteomic Analysis of Wheat Seeds during Artificial Ageing and Priming Using the Isobaric Tandem Mass Tag Labeling. PLoS ONE 11(9): e0162851. https://doi.org/10.1371/journal.pone.0162851
Editor: Pingfang Yang, Wuhan Botanical Garden, CHINA
Received: April 29, 2016; Accepted: August 28, 2016; Published: September 15, 2016
Copyright: © 2016 Lv et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All the relevant mass spectrometry proteomics data files are available from ProteomeXchange database (accession number(s) PXD004564, 10.6019/PXD004564).
Funding: This research was financially supported by grants from the Natural Science Foundation of China grants (31371850, 31501575, and 31401542), Natural Science Foundation of Henan Provincial Education Department (13A210164) and Natural Science Foundation of Youth Support Plan of Henan University of Technology (2015QNJH06).
Competing interests: The authors have declared that no competing interests exist.
Introduction
Wheat (Triticum aestivum L.), one of the most important, oldest and widely cultivated crops, is a staple food source for humans and livestock feed worldwide because of its high nutritional value [1, 2]. As orthodox type seeds, wheat seeds undergo desiccation after maturation, which enables to survive for a long time in a metabolic standstill situation [3]. As storage time is prolonged, seed vigor gradually decreases, and the germination rate eventually diminishes; as a consequence, commercial and genetic losses occur [4, 5]. Hence, seed ageing and germination mechanisms should be understood to develop new measures for seed conservation and production.
Seed ageing causes the physiological deterioration of seeds, which includes a reduced germination rate and an increased post-germination growth time [6, 7]. At present, the altered physiological and biochemical characteristics of seeds have been extensively investigated to elucidate seed aging mechanisms [8, 9]. Seed deterioration is mainly influenced by the accumulation of reactive oxygen species, lipid peroxidation mediated by free radicals, disruption of cellular membranes, and damage to proteins and nucleic acids [7, 8, 10–15]. Proteomic studies on artificially-aged Arabidopsis and Zea mays (maize) seeds indicated that differentially expressed proteins (DEPs) are mainly involved in oxidative stress, metabolism, and energy supply, which indicated that the proteomic changes can occur during deterioration at the dry state of aged seeds [4, 7]. Das and Sen-Mandi [16] further revealed that the physiological deterioration of wheat begins in its embryo, and this phenomenon is correlated with germination. Nevertheless, the mechanism underlying artificial ageing of wheat seeds remains unknown.
Crop seed germinability is a vital factor that contributes to seedling performance, plant establishment, and subsequent crop development and growth. Seed germination is controlled by both internal and external factors, including genetics, seed structure, seed chemistry, humidity, and temperature [17]. To improve and synchronize seed germination and emergence, researchers apply seed invigoration treatments called seed priming. Seed priming involves pretreatments with water and various chemical reagents, including polyethylene glycol, ascorbic acid, hormones, and vitamins [18, 19, 20, 21]. Proteomic investigations have been conducted during the seed germination of several plant seeds, such as wheat [22, 23], alfalfa [21], Arabidopsis [3], and maize [24]. These proteomic studies, conducted using two-dimensional (2-D) electrophoresis [21, 22] and 2-D differential gel electrophoresis [23], have provided critical information on the metabolic process of seed germination. However, 2-D-gel-based approaches suffer from low reproducibility and under-representation of low abundance and hydrophobic proteins [25]. These limitations can be overcome by a non-gel-based quantitative proteomic approach using isobaric tagging reagents. Isobaric tagging reagents, such as tandem mass tags (TMT) and isobaric tags for relative and absolute quantification (iTRAQ), have been developed for mass spectrometry (MS)-based protein detection and quantification in complicated samples [26, 27]. For instance, iTRAQ has been applied to conduct a quantitative proteomics study on wheat grain development and drought response [28, 29]. However, quantitative proteomics studies on wheat seed priming have yet to be reported.
The Chinese wheat cultivar ‘Aikang58’, a medium-hard wheat widely cultivated in the main wheat production areas of China, exhibits excellent characteristics, including drought tolerance, freeze resistance, wide adaptability, and high yield [30]. In this work, the first TMT-based quantitative proteome analysis of elite Chinese wheat cultivar ‘Aikang58’ seeds was conducted during artificial ageing and priming. We uncovered new information on the proteomic changes during ageing and priming in wheat seeds that might provide new insights into metabolic pathways, as well as into adverse defense mechanisms during seed deterioration and priming.
Materials and Methods
Wheat seeds, artificial ageing and priming treatment
The elite Chinese bread wheat cultivar ‘Aikang 58’ seeds used in this study were purchased from the Henan Academy of Agricultural Science. Seeds with similar sizes and weights were selected, and the germination rate (Gr) was calculated in accordance with the methods proposed by Dong et al. [23]. The original Gr was 98.0% (designated as WH98), and the seed moisture content was 9.79%.
Seeds were artificially aged by sealing them in air-tight plastic bottles and then stored at 45°C (± 1°C) and 50% relative humidity in a constant-temperature- and humidity-controlled cabinet (Binder KMF720, Tuttingen, Germany) in accordance with previous described methods [4] with minor modification. The vigor of the wheat seeds was determined regularly. The seed samples with three biological replicates were collected when Gr were 50%, 20%, and 1% (designated as WH50, WH20, and WH01 respectively). The collected samples were stored at -80°C prior to analysis. Unaged seeds (WH98) were used as control specimens.
Hydro-priming treatment was based on a previous method [3]. In brief, hydro-primed seeds (designated as WH100) were prepared by immersing dry mature seeds in sterilized water for 8 h at 25°C, packed in a wet gauze, and incubated for 12 h at 25°C.
Scanning electron microscopy (SEM)
SEM analysis was performed on the basis of a previous method [31]. WH98, WH50, WH20, and WH01 were halved vertically to the ventral side by using a cryostat (Leica CM1950, Solms, Germany). The cut side of embryo of the aged seeds was examined under a SEM (Quanta 250 FEG, FEI, Hillsboro, OR, USA).
Protein extraction and trypsin digestion
Embryo samples were collected by the dissection of WH98, WH50, WH20, WH01, and WH100 seeds as previously described [32], with each having three biological replicates. Afterward, the embryo samples of the wheat seeds were sonicated thrice on ice using a high intensity ultrasonic processor (Scientz, Ningbo, China) in lysis buffer (8 M urea, 1% Triton-100, 10 mM dithiothreitol and 1% protease inhibitor cocktail VI). After the samples were centrifuged at 20,000 × g for 10 min at 4°C, the supernatant was precipitated with 15% cold trichloroacetic acid for 2 h at -20°C and then centrifuged for 10 min at 4°C. The obtained precipitate was washed with cold acetone thrice, and re-dissolved in buffer (8 M urea and 100 mM triethylammonium bicarbonate, pH 8.0). Protein concentration was determined by using a 2-D Quant kit (GE Healthcare) according to the manufacturer’s instructions, and then stored at -80°C for further use.
Approximately 100 μg of proteins for each sample was digested with trypsin for the subsequent experiments. In briefly, proteins from wheat seed embryos were reduced with 10 mM dithiothreitol for 1 h at 37°C and alkylated with 20 mM indole-3-acetic acid for 45 min at room temperature in the dark. Finally, trypsin was added at 1:50 trypsin:protein mass ratio for the first digestion overnight and at a 1:100 trypsin:protein mass ratio for a second 4-h digestion.
TMT labeling and high-performance liquid chromatography (LC) fractionation
After trypsin digestion was completed, peptides were desalted using a Strata X C18 SPE column (Phenomenex, Torrance, CA, USA) and vacuum-dried. Peptides were reconstituted in 0.5 M triethylammonium bicarbonate and labeled. In brief, one unit of TMT reagent (labeled 100 μg of protein) was thawed and reconstituted in 24 μl acetonitrile. The peptide mixtures were incubated for 2 h at room temperature, pooled, desalted, and dried through vacuum centrifugation. Five samples with three biological replicates were labeled with TMT tags. WH100 and WH98 were labeled with 130 and 129 respectively. WH50, WH20, and WH01 were labeled with 128, 127 and 126respectively.
The TMT labeled samples were then fractionated through high pH reverse-phase high-performance liquid chromatography (LC) by using Agilent 300 Extend C18 columns (5 μm particles, 4.6 mm ID, 250 mm length). The LC gradient was run with 2% to 60% acetonitrile in 10 mM ammonium bicarbonate (pH 10) for 80 min to generate 80 fractions. Afterward, which all of the fractions were combined into 18 fractions. The fractionated samples were dried through vacuum centrifugation and stored at -20°C.
LC-MS/MS analysis
LC-MS/MS was performed on the basis of a previous report [33]. In brief, the peptides were dissolved in 0.1% formic acid, directly loaded onto a reversed-phase pre-column (Acclaim PepMap 100, Thermo Fisher Scientific, Waltham, MA, USA), and separated using a reversed-phase analytical column (Acclaim PepMap RSLC, Thermo Fisher Scientific). The gradient was increased from 7% to 20% solvent B (0.1% formic acid in 98% acetonitrile) for 24 min, from 20% to 35% in 8 min, and from 35% to 80% in 5 min. The gradient was then maintained at 80% for the last 3 min. A constant flow rate of 280 nl/min was set in an EASY-nLC 1000 UPLC system. The resulting peptides were analyzed by using Q Exactive™ hybrid quadrupole-Orbitrap mass spectrometer (Thermo Fisher Scientific).
The peptides were subjected to a nanospray ion source followed by tandem mass spectrometry (MS/MS) in Q Exactive™ (Thermo Fisher Scientific) coupled online to ultra performance LC. The peptides were then selected for MS/MS by using the NCE settings of 27, 30, and 33. A data-dependent procedure that alternated between one MS scan followed by 20 MS/MS scans was applied to the top 20 precursor ions above a threshold ion count of 1.0E4 in the MS survey scan with a 30.0 s dynamic exclusion. The fixed first mass was set at 100 m/z. The electrospray voltage and m/z scan range were described previously [33].
Database search and bioinformatics analysis
The resulting MS/MS data were processed using Mascot search engine (v.2.3.0 http://www.matrixscience.com/). Tandem mass spectra were used as query in the Uniprot_Triticum_aestivum database (100,981 sequences). Trypsin/P was specified as the cleavage enzyme that allows up to two missing cleavage sites. The mass error was set to 10 ppm for precursor ions and 0.02 Da for fragment ions. Carbamidomethyl on cysteine, TMT-6plex (N-term), and TMT-6 plex (K) were specified as fixed modifications, and oxidation on methionine was specified as a variable modification. The false discovery rate was adjusted to < 1%, and the peptide ion score was set at > 20. The proteins displaying a 1.2 fold change between artificially aged and normal seeds (WH01 vs WH98, WH20 vs WH98 and WH50 vs WH98) and between priming and normal seeds (WH100 vs WH98) were considered as DEPs if p < 0.05.
Protein annotation and functional analysis
Gene ontology (GO) (http://www.geneontology.org) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg/pathway.html) analyses were conducted in accordance with previously reported methods [34]. Proteins were subjected to a eukaryotic orthologous group (KOG) analysis by performing a homology search of the KOG database with the following parameters: E value < 1e-5, identities > 80%, and percent of match length > 60%. A domain annotation was performed using InterProScan on the InterPro (http://www.ebi.ac.uk/interpro/) domain database via Web-based interfaces and services [35]. Subcellular localization was determined by Wolfpsort (version of PSORT/PSORT II, http://psort.hgc.jp/).
GO, KEGG pathway, and protein domain enrichment analyses were performed, and a two-tailed Fisher’s exact test was employed to examine the enrichment of the DEPs against all of the identified proteins. Multiple hypothesis testing was corrected by using standard false discovery rate control methods, and domains with a corrected p-value < 0.05 were considered significant.
For hierarchical clustering based on different protein functional classifications (GO, domain, pathway), proteins from the categories were obtained after enrichment and then the categories were filtered to identify those that were at least enriched in one of the clusters with p-value < 0.05. This filtered p-value matrix was transformed by using the function x = −log10 (p-value), and x values were z-transformed for each functional category. Cluster membership was visualized by using a heat map via the “heatmap.2” function from “gplots” R-package.
Protein-protein interaction (PPI) network
All of the DEPs were used as query against the STRING database (version 10.0, http://string-db.org) [36] to identify protein-protein interactions. STRING defines a metric called “confidence score” to describe the interaction confidence. We obtained all of the interactions with a confidence score ≥ 0.7 (high confidence). The interaction network from STRING was visualized with Cytoscape (http://www.cytoscape.org/) [37]. Molecular complex detection (MCODE) was utilized to analyze densely connected regions. MCODE is part of the plug-in toolkit of the network analysis and visualization software Cytoscape.
RNA isolation and quantitative real-time PCR (qRT-PCR) analysis
Total RNA from wheat embryos of WH98, WH50, WH20, and WH01 were extracted by using RNAiso Plus reagent (Takara, Tokyo, Japan), and genomic DNA was removed by treating with DNase I (Takara) following the manufacturer’s protocol. Reverse transcription was performed by using a PrimeScriptRT reagent kit (Takara, RR047A) with a random primer mix in accordance with the manufacturer’s instructions. Primer pairs for qRT-PCR analysis (S3 Table) were designed by Primer5 and checked by querying the primer sequences, using the BLAST algorithm against the NCBI database. All of the primers were consistent with their respective target gene sequences. qRT-PCR of the translationally-controlled tumor protein, (TCTP; Q8LRM8), asparagine synthetase (AS; W5CK94) and catalase (CAT; W5HND1) genes was performed using Mastercycler ep realplex (Eppendorf, Hamburg, Germany) and SYBR Premix Ex Taq (Takara, RR420A) with ADP-ribosylation factor as a reference gene [38].
Results
Germinability and microstructural changes during artificial ageing of wheat seeds embryo
To investigate seed vigor and structural changes in embryos, Gr and SEM were employed to evaluate the alterations. There was no obvious difference in the Gr of seeds stored at room temperature for nearly 210 days, while the artificial ageing treatment resulted in a decreased germinability from 98% to 50% (65 d), 20% (131 d), and 1% (210 d). Gr was decreased rapidly from 30 (85%) to 90 d (41%), and almost no seeds could germinate after being stored 210 d (1%) at 45°C (Fig 1A). Gr also indicated that the deterioration rate accelerated after 30 d at 45°C. After priming was performed, the radicle protruded through the seed coat and the embryo was collected for further analysis (S1A Fig).
(A) Seed germination during artificial ageing. (B) Scanning electronic microscope observations of embryos from artificial ageing grains.
The SEM images (magnification was 4000×) of the cut side of embryos during artificial ageing were shown in Fig 1B. The results indicated that artificially aged seeds (WH50, WH20, and WH01) showed a degradation of granules and fold change on the surface compared with those of WH98. The altered microstructure and decreased Gr suggested that various physiological reactions might occur during this phase.
Overview of quantitative proteomics analysis
In this study, the TMT-based quantitative proteomic characteristics of the ‘Aikang58’ cultivar were investigated to uncover the altered proteins involved in artificial deterioration and priming. A global profile of the quantitative proteome was obtained by using the embryos from WH98, WH50, WH20, WH01 and WH100, with three biological replicates each. The SDS-PAGE analysis of proteins from the embryo was shown in S1B Fig. Each lane was loaded with 10 μg of proteins. When samples were subjected to an LC-MS/MS analysis, the data validation was evaluated. The distribution of mass error is near zero and most errors are less than 0.02 Da. These findings indicated that the MS data’s accuracy met the requirement (S1C Fig). The lengths of most of the peptides were distributed between 8 and 16 aa, which was consistent with the properties of tryptic peptides (S1D Fig), indicating that the sample preparation met the standard.
A total of 6281 proteins were identified, of which 3574 proteins with quantitative information were elicited from T. aestivum in the three biological replicates (S1 Table), and the coverage of identified proteins among the three biological repeats are also presented (S2A Fig). The protein expression levels were comparatively analyzed and divided into two groups: protein changes occurring during artificial ageing (WH01 vs WH98, WH20 vs WH98, and WH50 vs WH98), and protein changes occurring during priming (WH100 vs WH98). Pearson’s correlation coefficient was obtained to evaluate the repeatability of protein relative quantitation. The pair-wise Pearson’s correlation coefficients of the samples (protein changes during artificial ageing) are presented in a red-white-green heat map format (S2B Fig) and the reproducibility of the protein quantitation (protein changes during priming) of the samples is presented (S2C Fig). A 45°-diagonal line was obtained throughout the detection range, and this finding indicated the expected distribution without obvious changes among the three biological replicates (artificial ageing). The correlation coefficient was > 0.82 among the three biological replicates (priming), which indicated good reproducibility.
All the identified proteins were annotated (S2 Table), including GO, KEGG, KOG, and domain annotations. Of the 3574 quantified proteins, 162 were identified as DEPs, containing 36 up-regulated (≥1.2-fold, p-value ≤ 0.05) and 126 down-regulated (≤ 0.83-fold, p-value ≤ 0.05) in at least one of the artificial ageing stages compared with unaged seeds (Fig 2A, S3 Table). The 36 up-regulated proteins included 16 DEPs at WH01, 18 DEPs at WH20, and 10 DEPs at WH50, of which 8 DEPs were at two stages, meanwhile the 126 down-regulated proteins included 80 DEPs at WH01, 40 DEPs at WH20, and 57 at WH50, of which 10 DEPs were shared by all of the artificial ageing stages and 32 DEPs were shared by two stages (12 at WH01 and WH20, 12 at WH01 and WH50, and 8 at WH20 and WH50).
(A) Differentially expressed proteins during artificial ageing compared with unaged seeds; (B) sub-cellular localization analysis; (C) eukaryotic orthologous group (KOG) analysis.
Subcellular localization revealed the following: 31% of DEPs in the chloroplast, 22% of DEPs in the cytosol, 21% of DEPs in the nucleus, 11% of DEPs in the extracellular, and 9% of DEPs in the mitochondria (Fig 2B). The KOG annotation showed that DEPs were mainly involved in translation, ribosomal structure and biogenesis (13%), post-translational modification, protein turnover, chaperones (10%), energy production and conversion (9%), and signal transduction mechanisms (9%) (Fig 2C). GO analysis demonstrated that DEPs were mainly involved in binding, metabolic process, catalytic activity, and cellular process (S3A Fig).
To further understand the DEPs during artificial ageing, a functional enrichment analysis was performed (S4 Table). A Fisher’s exact test p-value was obtained to evaluate the enrichment analysis with those of the quantified proteins as references. Our results indicated that the proteins mainly involved in the extracellular region, serine-type endopeptidase inhibitor activity, negative regulation of peptidase, and endopeptidase and peptidase activities were down-regulated (S4A Fig), while the DNA-binding complex, nucleosome, DNA-packaging complex, protein heterodimerization and asparagine metabolism were accumulated in the GO enrichment analysis (S4B Fig). The KEGG pathway enrichment analysis showed that proteins involved in the ribosome were up-regulated, whereas proteins in starch, sucrose, ascorbate and aldarate metabolism were down-regulated (S4C, S5 and S6 Figs). Domain enrichment analysis indicated that the proteins involved in the histone core, histone fold, and asparagines synthase were over-represented in the up-regulated proteins, while proteins in the bifunctional inhibitor/plant lipid transfer protein/seed storage helical domain and bifunctional trypsin/alpha-amylase inhibitor helical domain were significantly over-represented in the down-regulated proteins (S4D and S4E Fig).
To better understand the dynamics of DEPs during artificial ageing, a hierarchical clustering analysis was performed to obtain the dynamic expression patterns. The up-regulated proteins, including carbohydrate derivatives, organonitrogen compounds, amino sugars, and aminoglycans, in the biological process category (Fig 3A) were mainly enriched in energy catabolic processes at the beginning of artificial ageing (from WH98 to WH50). With an extended artificial ageing time (WH01), these proteins were mainly enriched in amide biosynthesis. The down-regulated proteins were mainly implicated in defense- and stress-related pathways, and over-represented when Gr was nearly zero (WH01). The proteins related to lipid metabolic processes were accumulated when Gr was 20% (WH20). The results also showed that proteins involved in the negative regulation of metabolic processes were enriched in the later stages of artificial ageing (WH20 and WH01). The up-regulated DEPs in the cellular component category (Fig 3B) were mostly implicated in chromatin and nucleus, and were enriched at WH20 and WH01. By contrast, the proteins related to mitochondria and lipids were significantly down-regulated when Gr were 50% (WH50) and 1% (WH01), respectively. The up-regulated DEPs in the molecular function category (Fig 3C) were mainly related to chitinase and hydrolase activities, and clustered at WH50. By comparison, asparagine synthase was accumulated at WH01 and protein dimerization/heterodimerization activity was concentrated at WH20. The down-regulated proteins related to alpha-amylase inhibitor, lipid binding, and nutrient reservoir and enzyme inhibitor activities were enriched at WH20, while endopeptidase and its regulatory activity accumulated at WH20 and WH01. For DEPs in the KEGG pathway (Fig 3D), the up-regulated proteins related to the ribosome accumulated significantly at WH20, whereas the down-regulated proteins related to starch and sucrose metabolism were observed during the whole artificial ageing process. The proteins associated with ascorbate and aldarate metabolism were significantly enriched at WH50, and the amino sugar and nucleotide sugar metabolism were enriched at WH01. For DEPs in the protein domain cluster analysis (Fig 3E), the up-regulated proteins related to translation and histone were clustered at WH20, glutamine amidotransferase and asparagine synthase at WH01, and glycoside hydrolase at WH50. The down-regulated proteins concerned with plant lipid transfer and alpha-amylase inhibitor were enriched at WH20, whereas UDP-glucose/GDP-mannose dehydrogenase accumulated at WH50 and WH01 (Fig 3E).
(A) Biological process; (B) Cellular component; (C) Molecular function; (D) KEGG pathway; (D) Protein domain.
To illuminate the mechanisms of DEPs that mediated artificial ageing, the PPI network for DEPs were constructed with STRING and Cytoscape. The PPI networks of up and down-regulated proteins were also drawn separately (S7 Fig and S5 Table).
During seed priming, 531 DEPs, including 312 and 219 up-regulated and down-regulated proteins respectively, were presented (S3 Table). Subcellular localization revealed the following: 33% of DEPs were located in the cytosol, 32% of DEPs in the chloroplast, 14% of DEPs in the nucleus, 7% of DEPs in the mitochondria, and 5% of DEPs in the extracellular (Fig 4A). DEPs were mainly classified into five main categories, namely, posttranslational modification, protein turnover, and chaperones (18%); energy production and conversion (10%); amino acid transport and metabolism (9%); carbohydrate transport and metabolism (7%); and intracellular trafficking, secretion, and vesicular transport (7%) (Fig 4B). GO analysis revealed that DEP categories were mainly involved in metabolic process, binding, catalytic activity and cellular process (S3B Fig).
(A) Sub-cellular localization; (B) KOG categorization.
To better investigate the DEPs during seed priming, functional enrichment (GO, KEGG pathway, and protein domain) analyses were performed. A Fisher’s exact test p-value (-log10 [p-value]) was performed to evaluate the enrichment analysis, and the larger the value, the more DEPs were enriched in this category (Fig 5).
The red bar represents up-regulated proteins, and the green bar represents down-regulated proteins. (A, B) GO enrichment analysis; (C) KEGG enrichment analysis; (D) Protein domain enrichment analysis.
The up-regulated proteins significantly accumulated in the anion binding and biosynthetic process categories, including organic acid, carboxylic acid, small molecule, and metabolic processes involving oxoacid and carboxylic acid (Fig 5A). The proteins enriched in each category are listed in detail in S6 Table. In Fig 5B, the down-regulated proteins accumulated in the categories of nutrient reservoir activity, lipid particle and monolayer-surrounded lipid storage body. In Fig 5C, the down-regulated proteins mainly clustered in the protein processing sites in the endoplasmic reticulum (KEGG pathway: S8 Fig) in accordance with the KEGG pathway enrichment analysis. Meanwhile, the enriched KEGG pathways for up-regulated proteins included the phagosome (KEGG pathway: S9 Fig), biosynthesis of amino acids, amino sugar and nucleotide sugar metabolism, plant-pathogen interaction (KEGG pathway: S10 Fig), phenylalanine, tyrosine and tryptophan biosynthesis (KEGG pathway: S11 Fig), ascorbate and aldarate metabolism (KEGG pathway: S12 Fig), and pentose and glucuronate interconversions (KEGG pathway: S13 Fig). The proteins clustered in each pathway are listed in S6 Table. The domain enrichment analysis demonstrated that the down-regulated proteins contained seed maturation protein domains, alpha crystallin/Hsp20 domains, and HSP20-like chaperones. By contrast, the up-regulated proteins contained mainly the pyridoxal phosphate-dependent transferase, 90 kD heat shock protein, and 70 kD heat shock protein peptide-binding domains. The proteins clustered in each domain are listed in S6 Table.
To better elucidate the protein mechanisms involved in seed priming, we constructed the PPI network of DEPs using STRING and Cytoscape. With the MCODE plug-in toolkit, three enriched interaction clusters were associated with oxidative phosphorylation, carbon metabolism, and ATP binding (Fig 6 and S7 Table).
Differential protein expression during seed artificial ageing
On the basis of protein functions, these DEPs (S3 Table) could be grouped into four categories according to Xin et al. [4], except for the uncharacterized proteins. The four categories included protein destination and storage, metabolism/energy supply, defense/stress, and unclear classification (Table 1).
In Table 1, five storage proteins related with wheat flour quality (P10385, Q2A784, P08453, Q43659, and D2KFG9) and one regulatory protein (W5A8E0), were involved in protein destination and storage. The results indicated that storage proteins, except glutenin (P10385), were significantly changed at WH20 in all of the artificial ageing stages (WH50, WH20, and WH01) compared with WH98. Three kinds of proteins related to metabolism/energy, including tricarboxylic acid (TCA) cycle (malic enzyme, W5A810), sucrose synthase (W5B5R3 and W5I774) and starch branching enzyme 1 (O04074) were present and distributed in all of the stages. A total of 24 DEPs that significantly accumulated at the last artificial ageing stages WH20 and WH01 were categorized into defense/stress. These DEPs could be divided into three groups, namely, storage material degradation inhibitors (four alpha-amylase inhibitor and three oleosin proteins); response to oxidative, osmotic, and temperature pressures (catalase, UDP-glucose 6-dehydrogenase, TCTP, and multiprotein bridging factor), and response to biotic stresses (hemoglobin, agglutinin isolectin, five non-specific lipid-transfer proteins, ribosomal protein L3, and two asparagine synthetases). The results also indicated that a considerable number of DEPs involved in defense/stress were down-regulated.
Differential protein expression during seed priming
According to Dong et al. [23] and Ma et al. [28], DEPs (S3 Table) during seed priming, except the uncharacterized proteins, could be grouped into five categories, namely, energy supply and metabolism, signal transduction/cell growth, protein destination and storage, defense/stress, and unclear classification (Table 2).
In Table 2, proteins involved in energy supply could be mainly divided into four functional groups, the pentose phosphate pathway (PPP; W5AU57, W5B201, and W5HFZ6), glycolysis (W5FIN0 and W5CVS7), fatty acid oxidase (W5DVJ7), and TCA cycle (W5A810). The DEPs in metabolism could be mainly divided into three groups, fatty acid synthesis (B2ZGL2, W5AXD5, W5DZ64, W5F3T3, and W5FD02), protein synthesis (W5H9X1, D8L9B5, F4Y593, Q0Q0I7, W5F285, and W5GCA3), and amino acid synthesis (W5AVC9, W5CM54, W5DYX9, W5I8U3, Q45NB6, W5ARG6, and W5ECP3). The DEPs in signal transduction/cell growth could be mainly classified as those involved in phosphorylation (L0GED8), temperature response (P04464 and A4K4Y1), ubiquitin hydrolase (W5B4E9 and W5CWD6) and photosynthesis (P08823, W5B777, and W5ECJ8). The DEPs in defense/stress could be mainly divided into three groups, storage material degradation inhibitors (P16347, P01083, P16159, P17314, P04568, P22701, P42755, Q9ZR70, W5AW62, W5BE38, W5E9F6, W5FHA6, and W5FYW1); response to osmotic, temperature, and oxidative pressure (W5GDM5, C1K737, W5GEI6, Q6W8Q2, W5FJN1, W5FQ55, and W5FDW8); and response to biotic stresses (A1YE31, Q58A30, W5FEK8, Q1XHC6, O64393, Q8GZB0, P82901, and W5FSX7).
Transcriptional expression analysis as revealed by qRT-PCR
To provide further information on the correspondence between proteins and their mRNA expression patterns, qRT-PCR was performed to investigate the dynamic transcriptional expression patterns of three representative DEPs during artificial ageing. In Fig 7, the gene encoding asparagine synthetase displayed the same expression pattern as its protein. However, the two genes encoding TCTP and catalase showed opposite expression patterns than their proteins. Previous work demonstrated that the expression of TCTP is regulated at both the transcriptional and post-transcriptional levels [39], and such inconsistencies here might be due to post-transcriptional modification. The expression pattern of catalase is similar to that involved in wheat grain development [28], which is regulated by post-transcriptional regulation [40].
Solid lines represent mRNA expression patterns, and dotted lines represent protein expression patterns.
Discussion
The wheat embryo is a significant tissue where seed physiological deterioration and germination begin. In our study, a TMT labeling-based quantitative proteome was employed to detect protein changes during seed artificial ageing and priming. A total of 162 and 531 DEPs were identified during artificial ageing and seed priming, respectively. The DEPs were mainly involved in energy supply, metabolism, and stress response.
Effects of artificial ageing on wheat seed vigor and proteome
Wheat seed viability was gradually reduced as the Gr decreased from 98% to 1% after 210 d under the artificial ageing conditions used in this study. This effect was similar to the response of maize seeds [4]. In addition, the reduction in Gr was accompanied by structural changes in the embryo. Seed ageing affected the proteome of dry wheat seeds, and such effect is similar to that of maize seeds; this finding confirmed that proteome variation is possible in lowly hydrated seeds [4, 7].
DEPs involved in seed deterioration
Seed deterioration is a key developmental process during storage. When subjected to artificial ageing, proteins related to energy, storage nutrients, stress, and defense were altered. Several proteins previously reported in other processes were identified for the first time in the aged wheat seeds in this study. This study helps enhance our understanding of biochemical and molecular processes underlying seed deterioration.
Metabolism and energy supply.
Proteins related to glycolysis, TCA cycle, and PPP significantly contribute to seed vigor [41,42]. In our study, the abundance of 6-phosphofructokinase 1 (W5CP16 and W4ZRX8; the third enzyme and major point of regulation in glycolysis) and glyceraldehyde-3-phosphate dehydrogenase (W5BLN6; the sixth enzyme in glycolysis) decreased. Decreased glyceraldehyde-3-phosphate dehydrogenase levels are also observed in aged Arabidopsis seeds [7]; however, no protein participating in the PPP was detected in our experiment, which was similar to the results of studies on aged Arabidopsis and maize seeds [4,7]. The malic enzyme (malate dehydrogenase, W5A810), responsible for the TCA cycle replenishment pathway, is involved in catalyzing the conversion of pyruvate into malate; this enzyme was down-regulated by seed ageing (Table 1). The changes in malic enzyme was consistent with that in maize and Arabidopsis seeds during artificial ageing, indicating that the TCA cycle was disturbed by artificial ageing [4]. In our work, sucrose synthase (W5B5R3 and W5I774) participating in cell respiration and storage deposit [43], and starch branching enzyme 1 (O04074) involved in starch synthesis in plants [44], were also down-regulated. These results signified that the proteins involved in oligosaccharide metabolism were activated in a manner similar to those of maize seeds [4]. No protein involved in the PPP was detected in our study and therefore supported previous works on maize and aged Arabidopsis seeds [4,7]. The down-regulated KEGG pathways of energy supply were also enriched (Fig 3D). Kondoh et al. [45] found that enhanced glycolytic activation could inhibit cell senility by protecting against oxidative damage caused by reactive oxygen species, whereas the suppressed glycolytic pathway would result in aged cells. In maize (Z. mays) seeds, a controlled deterioration treatment inhibited the glycolytic pathway and promoted reactive oxygen species production and accumulation. Then, cell aging or death in seed embryos was accelerated, ultimately leading to seed deterioration and vigor loss [46]. Thus our results indicated that artificially aged wheat seeds triggered responses through the modification of the glycolytic pathway.
Storage protein decomposition.
Wheat grain proteins accumulate during seed development and decrease with seed ageing. Dell’Aquila previously reported protein change patterns in aged soybean and barley seeds. However, the specific proteins were not identified [47]. In this work, five kinds of storage proteins, particularly glutenin (P10385), avenin-like a1 (Q2A784), gamma-gliadin (P08453), 15 kD grain softness protein (Q43659), and gliadin/avenin-like seed protein (D2KFG9), were down-regulated. These storage proteins are closely related to the quality of wheat flour and thus provide energy for germination [48]. The break-down of these storage proteins might result in an inefficient material supply for the synthesis of new proteins essential for seed germination [16]. A deleterious effect on germination is associated with the protein-synthesizing system [47]. Our results also showed that 60S ribosomal protein L18 (W5A8E0) was up-regulated in aged seeds. Double-stranded RNA-activated protein kinase, a fundamental control step in the regulation of protein synthesis initiation, is negatively regulated by the 60S ribosomal subunit protein L18 in human cells [49]. Determining whether the regulation mechanism in wheat seeds varies will require further study. The KEGG enrichment analysis indicated that up-regulated proteins were associated with the ribosome (Fig 3D) where protein processing occurs. The results might be in accordance with the decreased proteins associated with artificial ageing.
Response to biotic and abiotic stimuli.
To protect themselves against various stresses, including a biotic or abiotic stimulus, such as oxidative stress, humidity, high temperatures, and pathogen penetration, a number of stress- and defense-related proteins changed during artificial ageing. The decrease in Gr was correlated with the increase of reactive oxygen species, the degradation of storage nutrients, and decreased activity levels of antioxidant enzymes and hydrolases inhibitors, as well as the content of antioxidants, such as L-ascorbic acid and glutathione, after the ageing treatment [5, 9, 50]. In our work, two proteins, catalase (W5HND1) and L-ascorbate peroxidase (W5DYW7), which is related to L-ascorbic acid, were down-regulated. This was consistent with the reduced mitochondrial antioxidant enzyme activities reported in aged soybean seeds [5]. In addition, the KEGG analysis indicated that two kinds of UDP-glucose 6-dehydrogenase (W5FQ55 and W5FDW8) involved in starch, sucrose, and ascorbate metabolism were also down-regulated (S5 and S6 Figs). Previous research indicated that UDP-glucose 6-dehydrogenase could catalyze the biosynthetic oxidation of UDP-glucose into UDP-glucuronic acid [51]. However, its function in oxidative responses has not been reported previously. The alpha-amylase inhibitors may protect stored starch against endogenous amylase activity and protect plants against wounding [23]. In this work, four proteins of alpha-amylase/trypsin inhibitors (P16159, P16851, P17314, and P01083) were down-regulated, indicating that the storage material began to degrade.
Hemoglobin is widely distributed in plants, and plays significant roles in improving the energy metabolism of hypoxic cells and resisting extreme external environmental conditions, such as low temperatures and pathogen attacks [52]. Mirelman et al. [53] reported that wheat germ agglutinin could inhibit fungal growth. In our research, hemoglobin 1 (M5AJV9) and agglutinin isolectin 3 (P10969) were both down-regulated, which signified a decrease in defensive responses against pathogens.
To date, only non-specific lipid-transfer proteins (nsLTPs) were found in plants and their seeds, whereas both LTPs and nsLTPs were discovered in animal tissues [54]. nsLTPs belong to the defense-associated proteins in plants because of their strong activity against plant pathogens [55]. nsLTPs in onion seeds inhibit the growth of Fusarium solani [56]; nsLTPs from barley and maize leaves inhibit the growth of Alternaria brassicola, Ascochyta pisi, Botrytis cinerea, Fusarium culmorum, and Verticillium dahliae [57]; and nsLTPs from Vigna unguiculata seeds inhibit the growth of Fusarium oxysporum and F. solani [58]. In our present work, five nsLTPs including P26913, Q2PCC3, W5FSX7, P82901 and W5D2W6, were down regulated, which indicated that the defense capability against pathogen penetration might decrease gradually.
TCTP, a highly conserved protein during evolution, is widely located in eukaryotic cells, including those of plants [59]. TCTP is an important component of the target of rapamycin signaling pathway, the major regulator of cell growth in animals and fungi. Despite its relevance, the molecular functions of plant TCTP homologs remain to be studied. Previous studies showed that the expression of TCTP was enhanced by stress stimuli, including high temperatures and salt stresses in cabbage (Brassica oleracea L.) [60], whereas in Arabidopsis, TCTP functions as an important growth regulator [59]. Additionally, in animal cells, TCTP functions as a novel heat shock protein with a chaperone-like activity [61]. The multiprotein bridging factor, encoded by TaMBF1c in wheat, plays a pivotal role in plant thermo-tolerance and may be useful in improving heat tolerance in other crops [62]. In our study, a TCTP homolog (Q8LRM8) was up-regulated, while the multiprotein bridging factor was down-regulated during the artificial ageing treatment at 45°C. The PPI network analysis indicated that TCTP might interact with other up-regulated proteins during artificial ageing. (S7 Fig and S5 Table). These results suggest that thermo-response factors were triggered when subjected to a high temperature and that the tolerance to high temperature in wheat seeds was balanced by various regulators.
In higher plants, asparagine synthesis is catalyzed by asparagine synthetase, which contributes to the protein synthesis and nitrogen metabolism [63]. Asparagine synthetase is also involved in responding to various factors, such as pathogen infection and oxidative stress [64]. In wheat, deoxynivalenol (DON) is a virulence factor of Fusarium, and resistance against DON is considered to be part of Fusarium resistance [65]. The evidence for resistance to DON by ribosomal protein L3 (RPL3) was claimed by Miller and Ewen [66] and Miller and Arnison [67]. In wheat, six RPL3 genes, three homologous versions of both paralogs, RPL3-A and RPL3-B, were expected to be found in T. aestivum [65]. In our present work, two AS (W5CK94 and W5FDU3) and two RPL3 (RPL3, A1YE31 and A1YE34) were up-regulated during artificial ageing, indicating that the protection system was triggered in wheat seeds to avoid pathogen penetration.
In plant seeds, storage triacylglycerols (TAG) are present in small discrete intracellular organelles called oil bodies. An oil body has a matrix of TAG, which is surrounded by phospholipids (PL) and alkaline proteins, termed oleosins [68]. Previous research showed that oleosins play important roles in maintaining the structure of the oil body during seed dehydration to facilitate seed germination, promoting the decomposition rate of the oil body during germination, and providing a binding site for lipase during germination [69,70]. In this work, three oleosins (W5AW62, W5FHA6, and W5FYW1) were down-regulated, which might be in accordance with the reduced Gr.
DEPs participate in seed priming
The priming treatment could accelerate seed germination and improve seedling uniformity. Upon hydro-priming, embryonic cells switch from a quiescent state to a highly active metabolic state. In the present study, many identified proteins might aid in characterizing germination vigor and optimizing germination enhancement treatments.
Activated energy supply.
The main physiological characteristics during seed priming and subsequent germination are storage degradation, active metabolism, and cell growth [71]. Upon priming with water, storage substances begin to decompose and the central carbon metabolism becomes highly active. Enzymes in glycolysis, the PPP and the TCA cycle are activated during this process. Triosephosphate isomerase (TPI, W5FIN0) is essential for efficient energy production; the enzyme catalyzes the reversible interconversion of dihydroxyacetone phosphate and D-glyceraldehyde 3-phosphate. Meanwhile, pyruvate kinase (W5CVS7), the final enzyme in glycolysis, catalyzes phosphoenolpyruvate to pyruvate, yielding one molecule of ATP [72]. Glucose-6-phosphate-dehydrogenase (W5AU57 and W5B201), the first and also rate-limiting enzyme in PPP, catalyzes the formation of glucose 6-phosphate to form 6-phosphoglucono-δ-lactone, an intramolecular ester that undergoes oxidation. Decarboxylation catalyzed by 6-phosphogluconate dehydrogenase (W5HFZ6), the second enzyme up-regulated in PPP, produces the ketopentose ribulose 5-phosphate (Table 2). The up-regulated malic enzyme (W5A810) in priming catalyzes pyruvate into malate, which is then converted into oxaloacetate, which replenishes the TCA cycle. Our data indicated that the proteins involved in glycolysis, PPP, and TCA were all up-regulated, which was in agreement with those of wheat seed germination [22, 23]. In addition, fructan 1-exohydrolase w2 (1-FEH w2, Q84LA1) was up-regulated, being considerably higher during the fructan breakdown phase [73].
Besides the proteins indicated above, fatty acid oxidation is another energy supply pathway. The first step of β-oxidation in the peroxisome is catalyzed by the enzyme acyl-CoA oxidase (ACOX, W5DVJ7). This process is also the rate-limiting step in peroxisomal β-oxidation [74]. The up-regulated level of acyl-CoA oxidase indicated that energy was also provided by fatty acid decomposition, which had not been previously reported. These observations demonstrated that energy metabolism was triggered after imbibition to provide energy for subsequent seed germination.
The storage proteins, including globulin-3A (I6QQ39), gamma-gliadin (P04730), and glutenin, and the low molecular weight subunit (P10385) were down-regulated. This finding indicated that they were decomposed to produce energy for seed germination.
The PPI network analysis indicated that the enriched clusters were involved in oxidative phosphorylation and carbon metabolism. Oxidative phosphorylation is the culmination of energy yielding metabolism. All of the oxidative steps in the degradation of carbohydrates, fats, and amino acids converge at this final stage of cellular respiration, in which the energy of oxidation drives the synthesis of ATP [72]. The PPI results indicated that carbohydrates and their oxidative phosphorylation were highly active, and provided energy during seed priming.
Physiological metabolism and organizational development.
Various physiological metabolic pathways, as well as organizational development, became more active, with fatty acid synthesis, amino acid synthesis, protein translation, and protein folding occurring during this phase. In both eukaryotes and prokaryotes, acetyl-CoA carboxylase is a biotinylated enzyme that catalyzes the first committed step of de novo fatty acid biosynthesis by the carboxylation of acetyl-CoA to malonyl-CoA. In plants, acetyl-CoA carboxylase activity is found in the plastids where primary fatty acid biosynthesis occurs [75]. In plants, monounsaturated fatty acids maintain membrane fluidity and are also sources of storage triacylglycerols in tissues [76]. Acyl-acyl carrier protein (ACP) desaturases are responsible for the conversion of saturated fatty acyl-ACPs into monounsaturated-ACPs. This process is the penultimate step in fatty acid biosynthesis and produces the main source of monounsaturated fatty acids in plant tissues [72,77]. ACP transports the growing fatty acid chain between enzymatic domains of fatty acid synthase during fatty acid biosynthesis [78, 79]. In this work, the accumulation of acetyl-CoA carboxylase (W5DVJ7) involved in fatty acid synthesis was consistent with results of previous research [22]. Acyl-ACP (W5AXD5) and ACP (W5DZ64, W5F3T3, and W5FD02) involved in fatty acid synthesis have not been reported previously during wheat seed priming.
During seed priming and germination, amino acids and proteins accumulate. The aromatic amino acids, as well as other aromatic compounds, which function in plant defense, electron transport, signaling, communication, and wound responses, were synthesized via the shikimate pathway [80, 81]. In our present work, 3-phosphoshikimate 1-carboxyvinyltransferase and chorismate synthase were accumulated. The enzyme 5-enolpyruvylshikimate- 3-phosphate synthase, encoded by aroA (3-phosphoshikimate 1-carboxyvinyltransferase), is the sixth enzyme in the shikimate pathway and catalyzes the conversion of phosphoenol-pyruvate and shikimate-3-phosphate into 5-enolpyruvylshikimate 3-phosphate. Meanwhile, chorismate synthase, the seventh enzyme in the shikimate pathway catalyzes the transformation of 5-enolpyruvylshikimate 3-phosphate to the last common precursor, chorismate, in the biosynthesis of numerous aromatic compounds in bacteria, fungi, and plants [72,80,82]. The up-regulated proteins that participated in phenylalanine, tyrosine, and tryptophan biosynthesis are shown in S11 Fig. In addition, enzymes involved in histidine biosynthesis (W5AVC9; histidinol dehydrogenase, last step in histidine biosynthesis) [83], cysteine biosynthesis (W5CM54; cysteine synthase) [84], proline biosynthesis (Q45NB6 and W5ARG6; glutamine synthetase and glutamate dehydrogenase, respectively) [85], and glutamate and asparagine biosynthesis (W5ECP3; N-acetyl-gamma-glutamyl-phosphate reductase) [86] were also accumulated. The results suggested that the level of activated amino acid synthesis was similar to that during wheat germination [23].
The active amino acid synthesis pathway provided raw materials for protein synthesis, which might result in the induction of the unfolded protein response (UPR) pathway. In our study, the translation controlling protein eukaryotic translation initiation factor 3 (W5F285 and W5GCA3) [87] was up-regulated. UPR-relevant marker proteins, such as endoplasmic reticulum chaperone BiP (W5GCC9) and protein disulfide isomerase (D8L9B5) were also accumulated during priming. In yeast, the dissociation of Kar2p/BiP from an endoplasmic reticulum sensory molecule, Ire1p, triggers the UPR. The accumulation of BiP suggested that the UPR response might be triggered after priming. Protein disulfide isomerase protein, a multifunctional protein that resides in the endoplasmic reticulum’ lumen, was up-regulated. This enzyme may hold two resulting functions, the formation of molecules with disulfide bonds in wheat, and response to protein stress in the endoplasmic reticulum [88, 89]. Moreover, the highly conserved molecular chaperone 90 kD heat shock protein (F4Y593 and Q0Q0I7), which fulfills a housekeeping function by contributing to protein folding, and a regulatory function by preventing unfolded proteins from aggregating [90], was also up-regulated. Reticulon-like protein, required for endoplasmic reticulum organization [91], was down-regulated in our study, and its function remains to be studied further. The DEPs that participated in protein processing in the endoplasmic reticulum are presented in S8 Fig. In our present work (Table 2), adenylosuccinate synthetase (W5DWA6; involved in purine biosynthesis) and lipoxygenase (W5F9D7; involved in the determination of wheat grain quality) [92] were also up-regulated.
Signal transduction/ cell growth.
The abundance of 14-3-3 proteins (L0GED8 in Table 2) increased after priming, which was in accordance with a previous study on wheat germination [22]. These proteins are involved in signal transduction chains and facilitate the phosphorylation of target proteins [93]. A calmodulin, which is a type of protein involved in heat shock signal transduction in wheat [94], was found in this study. The α- and β-tubulins are considered important structural elements in cell growth, morphogenesis, and response to cold acclimation [95, 96]. β-Tubulin has been observed in wheat during early germination [93]. In our study, an up-regulated α-tubulin-2A was present, which was in agreement with the report of Gallardo et al. [3]. Two up-regulated ubiquitin carboxyl-terminal hydrolases (W5B4E9 and W5CWD6) were discovered in our study. The ubiquitin carboxyl-terminal hydrolases are a subset of de-ubiquitinating proteases that control biological activity, and regulatory and structural proteins in both plants and animals. In Arabidopsis, it contributes to the modification of plant shoot architecture [42]. Another accumulated protein associated with radicle protrusion, S-adenosylmethionine synthase (B0LXM0 and W5G6F1), is synthesized from the pool of stored mRNAs and is involved in the formation of S-adenosyl-methionine, which is essential for plant growth and embryo development [97]. A ribulose-1, 5-bisphosphatecarboxylase/oxygenases large subunit (P08823), two serine hydroxymethyl transferases (W5B777 and W5ECJ8), plasma membrane ATPase (P83970) and fructose-bisphosphate aldolase (W5D043), which are involved in dark reaction/photorespiration/phototropic growth, were found [23, 98–100]. In addition, an up-regulated ferritin (W5ETC7; involved in mineral transport) was also observed [101,94]. These proteins contributed to the onset of the subsequent development of the wheat embryo.
Stress response and defense.
During seed priming and subsequent germination, various defense proteins stored in mature dry seeds are activated by environmental conditions, including moisture, oxygen, temperature, and oxidative stress [23]. Unlike the down-regulated proteins of L-ascorbic acid, and UDP-glucose 6-dehydrogenase (W5DP16, W5FDW8, and W5FQ55) participating in ascorbate and aldarate metabolisms in aged seeds, they were up-regulated to defend against oxidative stress during priming (S12 Fig and S6 Table). Mak [22] observed marked decreases in the abundances of hydrolytic enzyme inhibitors, including alpha-amylase inhibitor and trypsin inhibitor, in embryos during germination. By contrast, Dong [23] reported that the alpha-amylase inhibitor was abundantly expressed to protect plants against wounding and is beneficial for organ development. In our present work, endogenous alpha-amylase/subtilisin inhibitor (P16347) and alpha-amylase inhibitor 0.28 (P01083) were down-regulated, while alpha-amylase/trypsin inhibitor CM16 (P16159) and CM3 (P17314) were up-regulated. The Em proteins, synthesized at the very early stage in the wheat embryo, may provide a matrix of bound water, which opposes the protein denaturation [102]. In our study, the four Em proteins (P04568, P22701, P42755, and Q9ZR70) were decreased, indicating that protein degradation might be activated. Kibinza [103] indicated that priming induces the synthesis of catalase, which is involved in seed recovery during priming. Up-regulated catalase (W5GEI6) levels were also observed in our study. Previous studies reported that 1-cys peroxiredoxin (Q6W8Q2), which functions as a molecular chaperone under oxidative stress conditions, is down-regulated during seed germination [22]. In our work, 1-cys peroxiredoxin was also down-regulated. Seeds may be exposed to various microorganisms during germination, such as bacterial and fungi. In our present work, several DEPs associated with biotic stresses were identified. These DEPs included ribosomal protein L3-A3 (A1YE31; resistance to DON) [66], P450 (Q58A30; resistance to Fusarium head blight) [104], methylenetetrahydrofolate reductase (W5FEK8; infection-related morphogenesis) [105], multidomain cystatin (Q1XHC6; inhibit growth of the snow mold fungus) [106] and Wheatwin-2 (O64393; antifungal activity toward B. cinerea, F. culmorum and Fusarium graminearum) [107], and nsLTPs (Q8GZB0, P82901, and W5FSX7) [57]. Unlike in artificial ageing, a decreased expression level of multiprotein bridging factor 1 (C1K737), a plant thermo-tolerance protein [62], was observed at normal temperature (25°C) compared with that at higher artificial temperature (45°C). The TCTP was up-regulated during priming, which is consistent with the observation in Arabidopsis, in which it functions as an important growth regulator [59]. Chalcone-flavanone isomerase (W5FJN1) in flavonoid biosynthesis has a well-recognized role in the stress response of bread wheat [108–110] and was down-regulated in our research. In addition, proteins, including the up-regulated 3-ketoacyl-CoA synthase (W5GDM5) [111], which is involved in osmotic stress, were also observed in our study. Oleosins, which can promote the decomposition rate of the oil body during germination [69], were down-regulated (W5AW62, W5BE38, W5E9F6, W5FHA6, and W5FYW1), suggesting that oil body decomposition might occur.
Different physiological states during wheat seed artificial ageing and priming
Seed ageing and seed priming are two distinguishable physiological states of wheat seed. During seed ageing, wheat seeds undergo a metabolic pause and wheat vigor progressively decreases. Upon hydro-priming, embryonic cells switch from a quiescence state to a highly active metabolism state and seed vigor increases rapidly. Therefore, it was important to compare the biochemical behaviors of these two distinct states.
Our present proteome analysis revealed that the enzymes that participated in energy supply displayed opposite changes in the artificial aged and hydro-primed seeds. Enzymes involved in glycolysis (6-phosphofructokinase and glyceraldehyde-3-phosphate dehydrogenase) and the TCA cycle (malic enzyme) were strongly reduced during artificial ageing (Table 1). However, the opposite was observed for proteins involved in glycolysis (triosephosphate isomerase), the PPP pathway (glucose-6-phosphate-dehydrogenase), and the TCA cycle (malic enzyme). Such proteins exhibited strongly increased levels during hydro-priming (Table 2).
Other significant changes observed between ageing and priming involved the stress and response patterns of the seed proteome. Because of the accumulation of reactive oxygen species during seed ageing, antioxidant enzymes, as well as the enzymes participating in antioxidant biosynthesis, including catalase, L-ascorbate peroxidase and UDP-glucose 6-dehydrogenase, were down-regulated. During priming, defense proteins were activated, including UDP-glucose 6-dehydrogenase and catalase. Although our data enumerate the DEPs between ageing and priming, nsLTPs, TCTP, and oleosin remains to be studied further.
Finally, it is remarkable that catabolism was present in the aged seeds, while anabolism was observed in the primed seeds. The abundances of the five kinds of storage proteins including glutenin, avenin-like a1, gamma-gliadin, 15 kD grain softness protein, and gliadin/avenin-like seed protein were decreased. In the primed seeds, several enzymes involved in fatty acid (acetyl-CoA carboxylase, ACP and acyl-ACP), amino acid (shikimate pathway, histidine and proline) and protein synthesis were up-regulated. The up-regulated levels of the endoplasmic reticulum chaperone BiP and protein disulfide isomerase in the UPR pathway confirmed the active status of protein synthesis.
Conclusions
In summary, wheat embryo during artificial ageing and seed priming was subjected to quantitative proteomic analyses. A total of 162 and 531 DEPs were respectively identified during artificial ageing and seed priming of the ‘Aikang58’ cultivar, and the quantitative expression characterization, functional analysis and protein-protein interactions of these DEPs were investigated.
The DEPs present during artificial ageing participated in various cellular processes, such as protein destination and storage, metabolism and energy supply, and defense/stress, suggesting that the artificial ageing affected these pathways. The reduced ability to protect against ageing may lead to an increase in amylase and protease levels, decomposition of storage substances, impairment of metabolism and energy supply, and ultimately seed deterioration. A considerable number of proteins involved in metabolism and energy supply (malic enzyme, glyceraldehydes-3phosphate dehydrogenase, and 6-phosphofructokinase 1) and defense/stress (catalase, L-ascorbate peroxidase, hemoglobin 1, oleosin, and nsLTPs) were down-regulated. These proteins may play important roles in wheat seed deterioration and might be considered new protein markers for labeling seed ageing. Unlike in seed ageing, most of the DEPs involved in metabolism, energy supply, and signal transduction/ cell growth, were activated. A considerable number of proteins involved in energy supply (glycolysis, the TCA cycle, and the PPP), metabolism (fatty acid oxidation and synthesis, and amino acid synthesis), signal transduction/cell growth (14-3-3 proteins, heat and cold responses), and photorespiration were up-regulated. This pattern is consistent with that of highly active metabolic state. Our results provide comprehensive proteome insights into protein changes that occur during seed deterioration and priming.
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [112] partner repository with the dataset identifier PXD004564 and 10.6019/PXD004564.
Supporting Information
S1 Fig. Morphology of seed priming (A) SDS-PAGE analysis of seed embryos proteins during different stages; (B) Mass error distribution of all identified peptides (C); Peptide length distribution (D).
https://doi.org/10.1371/journal.pone.0162851.s001
(TIF)
S2 Fig. Venn diagrams of the identified proteins between three biological repeats (A); Pair-wise pearson’s correlation coefficient of three biology replicates during artificial ageing (B); Protein quantitation reproducibility of three repeat experiment during priming (C).
https://doi.org/10.1371/journal.pone.0162851.s002
(TIF)
S3 Fig. Gene Ontology analysis of DEPs during artificial ageing (A) and priming (B).
https://doi.org/10.1371/journal.pone.0162851.s003
(TIF)
S4 Fig. Functional enrichment analysis of DEPs in artificial ageing.
Red bars indicate the up-regulated proteins; green bars indicate down-regulated proteins. GO Ontology enrichment (A, B); KEGG enrichment (C); Domain enrichment (D, E).
https://doi.org/10.1371/journal.pone.0162851.s004
(TIF)
S5 Fig. Down-regulated proteins (highlighted green boxes) participated in starch and sucrose metabolism during artificial ageing.
https://doi.org/10.1371/journal.pone.0162851.s005
(TIF)
S6 Fig. Down-regulated proteins (highlighted green boxes) participated in ascorbate and aldarate metabolism during artificial ageing.
https://doi.org/10.1371/journal.pone.0162851.s006
(TIF)
S7 Fig. The protein-protein interaction network analysis of up-regulated proteins (A) and down-regulated proteins (B) identified by TMT-labeling during seed artificial ageing.
https://doi.org/10.1371/journal.pone.0162851.s007
(TIF)
S8 Fig. DEPs participated in protein processing in endoplasmic reticulum during seed priming.
Red boxes indicate the up-regulated proteins; green boxes indicate down-regulated proteins
https://doi.org/10.1371/journal.pone.0162851.s008
(TIF)
S9 Fig. Up-regulated proteins (red highlighted boxes) involved in phagosome during priming.
https://doi.org/10.1371/journal.pone.0162851.s009
(TIF)
S10 Fig. Up-regulated proteins (red highlighted boxes) involved in plant-pathogen interaction during priming.
https://doi.org/10.1371/journal.pone.0162851.s010
(TIF)
S11 Fig. Up-regulated proteins (red highlighted boxes) involved in phenylalanine, tyrosine and tryptophan biosynthesis during priming.
https://doi.org/10.1371/journal.pone.0162851.s011
(TIF)
S12 Fig. Up-regulated proteins (red highlighted boxes) involved in ascorbate and aldarate metabolism during priming.
https://doi.org/10.1371/journal.pone.0162851.s012
(TIF)
S13 Fig. Up-regulated proteins (red highlighted boxes) involved in pentose and glucuronate interconversions during priming.
https://doi.org/10.1371/journal.pone.0162851.s013
(TIF)
S1 Table. The identified proteins of three biological replicates and combined data.
https://doi.org/10.1371/journal.pone.0162851.s014
(XLS)
S2 Table. Annotation of the identified proteins.
https://doi.org/10.1371/journal.pone.0162851.s015
(XLS)
S3 Table. Differentially expressed proteins (DEPs) during artificial ageing and priming and primers of DEPs encoding genes for qRT-PCR.
https://doi.org/10.1371/journal.pone.0162851.s016
(XLS)
S4 Table. Functional enrichment analysis (GO, KEGG and Domain enrichment) of DEPs during artificial ageing.
https://doi.org/10.1371/journal.pone.0162851.s017
(XLS)
S5 Table. DEPs involved in protein-protein interaction networks during artificial ageing.
https://doi.org/10.1371/journal.pone.0162851.s018
(XLS)
S6 Table. Functional enrichment analysis (GO, KEGG and Domain enrichment) of DEPs during priming.
https://doi.org/10.1371/journal.pone.0162851.s019
(XLS)
S7 Table. DEPs involved in protein-protein interaction networks during priming.
https://doi.org/10.1371/journal.pone.0162851.s020
(XLS)
Acknowledgments
This research was financially supported by grants from the Natural Science Foundation of China (31371850, 31501575, and 31401542), Natural Science Foundation of Henan Provincial education Department (13A210164) and Natural Science Foundation of Youth Support Plan of Henan University of Technology (2015QNJH06).
Author Contributions
- Validation: YL.
- Writing – original draft: YL SZ.
- Writing – review & editing: JW YH.
References
- 1. Šramková Z, Gregová E, Šturdík E (2009) Chemical composition and nutritional quality of wheat grain. Acta Chimica Slovaca, 2: 115–138.
- 2. Brenchley R, Spannagl M, Pfeifer M, Barker GL, D’Amore R, Allen AM, et al. (2012) Analysis of the bread wheat genome using whole-genome shotgun sequencing. Nature, 491: 705–710. pmid:23192148
- 3. Gallardo K, Job C, Groot SP, Puype M, Demol H, Vandekerckhove J, et al. (2001) Proteomic analysis of Arabidopsis seed germination and priming. Plant Physiol, 126:835–848. pmid:11402211
- 4. Xin X, Lin XH, Zhou YC, Chen XL, Liu X, Lu XX (2011) Proteome analysis of maize seeds: the effect of artificial ageing. Physiol plantarum, 143: 126–138.
- 5. Xin X, Tian Q, Yin G, Chen X, Zhang J, Ng S, et al. (2014) Reduced mitochondrial and ascorbate-glutathione activity after artificial ageing in soybean seed. J Plant Physiol, 171: 140–147. pmid:24331429
- 6. Ellis RH (1992) Seed and seedling vigour in relation to crop growth and yield. Plant Growth Regul, 11: 249–255.
- 7. Rajjou L, Lovigny Y, Groot SP, Belghazi M, Job C, Job D (2008) Proteome-wide characterization of seed aging in Arabidopsis: a comparison between artificial and natural aging protocols. Plant Physiol, 148: 620–641. pmid:18599647
- 8. Bailly C, Kranner I (2011) Analyses of reactive oxygen species and antioxidants in relation to seed longevity and germination. Methods Mol Biol, 773:343–67. pmid:21898265
- 9. Morscher F, Kranner I, Arc E, Bailly C, Roach T (2015) Glutathione redox state, tocochromanols, fatty acids, antioxidant enzymes and protein carbonylation in sunflower seed embryos associated with after-ripening and ageing. Ann Bot, 116: 669–678. pmid:26346716
- 10. Walters C (1998) Understanding the mechanisms and kinetics of seed aging. Seed Sci Res, 8: 223–244.
- 11. McDonald MB (1999) Seed deterioration: physiology, repair and assessment. Seed Sci Technol, 27:177–237.
- 12. Bellani LM, Salvini L, Dell’Aquila A, Scialabba A (2012) Reactive oxygen species release, vita-min E, fatty acid and phytosterol contents of artificially aged radish (Raphanussativus L.) seeds during germination. Acta Physiol Plant, 34:1789–99.
- 13. Hu D, Ma G, Wang Q, Yao J, Wang Y, Pritchard HW, et al. (2012) Spatial and temporal nature of reactive oxygen species production and programmed cell death inelm (Ulmus pumila L.) seeds during controlled deterioration. Plant Cell Environ, 35:2045–59. pmid:22582978
- 14. Yao Z, Liu L, Gao F, Rampitsch C, Reinecke DM, Ozga JA, et al. (2012) Developmental and seed aging mediated regulation of antioxidative genes and differential expres-sion of proteins during pre-and post-germinative phases in pea. Plant Physiol, 169:1477–88.
- 15. Lee J, Welti R, Roth M, Schapaugh WT, Li J, Trick HN (2012) Enhanced seed viability and lipid compositional changes during natural ageing by suppressing phospholipase Dα in soybean seed. Plant Biotechnol J, 10: 164–173. pmid:21895945
- 16. Das G, Sen-Mandi S (1988) Root formation in deteriorated (aged) wheat embryos. Plant physiol, 88: 983–986. pmid:16666491
- 17. Lei L, Zhu X, Wang S, Zhu M, Carver BF, Yan L (2013) TaMFT-A1 is associated with seed germination sensitive to temperature in winter wheat. PloS one, 8(9), e73330. pmid:24069187
- 18. Heydecker W, Higgins J, Gulliver RL (1973) Accelerated germination by osmotic seed treatment. Nature, 246: 42–44.
- 19. Chhetri DR, Rai AS, Bhattacharjee A (1993) Chemical manipulation of seed longevity of four crop species in an unfavourable storage environment. Seed Sci Technol, 21: 31–44.
- 20. Powell AA, Yule LJ, Jing H-C, Groot SPC, Bino RJ, Pritchard HW (2000) The influence of aerated hydration seed treatment on seed longevity as assessed by the viability equations. J Exp Bot 51: 2031–2043. pmid:11141177
- 21. Yacoubi R, Job C, Belghazi M, Chaibi W, Job D (2011) Toward characterizing seed vigor in alfalfa through proteomic analysis of germination and priming. J Proteome Res, 10: 3891–3903. pmid:21755932
- 22. Mak Y, Willows RD, Roberts TH, Wrigley CW, Sharp PJ, Copeland L (2009) Germination of wheat: a functional proteomics analysis of the embryo. Cereal Chem, 86: 281–289.
- 23. Dong K, Zhen S, Cheng Z, Cao H, Ge P, Yan Y (2015) Proteomic analysis reveals key proteins and phosphoproteins upon seed germination of wheat (Triticum aestivum L.). Frontiers in Plant Science, 6:1017. pmid:26635843
- 24. Fu Z, Jin X, Ding D, Li Y, Fu Z, Tang J (2011) Proteomic analysis of heterosis during maize seed germination. Proteomics, 11: 1462–1472. pmid:21365753
- 25. Ge P, Hao P, Cao M, Guo G, Lv D, Subburaj S (2013). iTRAQ-based quantitative proteomic analysis reveals new metabolic pathways of wheat seedling growth under hydrogen peroxide stress. Proteomics, 13: 3046–3058. pmid:23929510
- 26. van Ulsen P, Kuhn K, Prinz T, Legner H, Schmid P, Baumann C (2009) Identification of proteins of Neisseria meningitidis induced under iron-limiting conditions using the isobaric tandem mass tag (TMT) labeling approach. Proteomics, 9: 1771–1781. pmid:19288525
- 27. Karp NA, Huber W, Sadowski PG, Charles PD, Hester SV, Lilley KS (2010) Addressing accuracy and precision issues in iTRAQ quantitation. Mol Cell Proteomics, 9: 1885–1897. pmid:20382981
- 28. Ma C, Zhou J, Chen G, Bian Y, Lv D, Li X, et al. (2014) iTRAQ-based quantitative proteome and phosphoprotein characterization reveals the central metabolism changes involved in wheat grain development. BMC genomics, 15:1029. pmid:25427527
- 29. Ford KL, Cassin A, Bacic A (2011). Quantitative proteomic analysis of wheat cultivars with differing drought stress tolerance. Front Plant Sci, 2: 44. pmid:22639595
- 30. Zeng J, Li G, Gao H, Ru Z (2011) Comparison of A and B starch granules from three wheat varieties. Molecules, 16: 10570–10591. pmid:22183883
- 31. Dong K, Ge P, Ma C, Wang K, Yan X, Gao L, et al. (2012) Albumin and globulin dynamics during grain development of elite Chinese wheat cultivar Xiaoyan 6. J Cereal Sci, 56: 615–622.
- 32. Yu Y, Zhen S, Wang S, Wang Y, Cao H, Zhang Y, et al. (2016) Comparative transcriptome analysis of wheat embryo and endosperm responses to ABA and H2O2 stresses during seed germination. BMC genomics, 17: 1.
- 33. Jin W, Wu F (2016) Proteome-Wide Identification of lysine succinylation in the proteins of tomato (Solanum lycopersicum). PloS one, 11(2), e0147586. pmid:26828863
- 34. Wu J, Mao X, Cai T, Luo J, Wei L (2006) KOBAS server: a web-based platform for automated annotation and pathway identification. Nucleic Acids Res, 34: W720–W724. pmid:16845106
- 35. Zdobnov EM, Apweiler R (2001) InterProScan-an integration platform for the signature-recognition methods in InterPro. Bioinformatics, 17: 847–848. pmid:11590104
- 36. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. (2014) STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res, gku1003.
- 37. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. (2003) Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res, 13: 2498–2504. pmid:14597658
- 38. Paolacci AR, Tanzarella OA, Porceddu E, Ciaffi M (2009) Identification and validation of reference genes for quantitative RT-PCR normalization in wheat. BMC Mol Biol, 10: 1–27.
- 39. Xu A, Bellamy AR, Taylor JA (1999) Expression of translationally controlled tumour protein is regulated by calcium at both the transcriptional and post-transcriptional level. Biochem J, 342: 683–689. pmid:10477280
- 40. Ni W, Trelease RN (1991) Post-transcriptional regulation of catalase isozyme expression in cotton seeds. The Plant Cell, 3: 737–744. pmid:12324611
- 41. Kim ST, Wang Y, Kang SY, Kim SG, Rakwal R, Kim YC, et al. (2009) Developing rice embryo proteomics reveals essential role for embryonic proteins in regulation of seed germination. J Proteome Res, 8: 3598–3605. pmid:19472976
- 42. Yang P, Smalle J, Lee S, Yan N, Emborg TJ, Vierstra RD (2007) Ubiquitin C-terminal hydrolases 1 and 2 affect shoot architecture in Arabidopsis. The Plant J, 51: 441–457. pmid:17559514
- 43. Abid G, Silue S, Muhovski Y, Jacquemin JM, Toussaint A, Baudoin JP (2009) Role of myo-inositol phosphate synthase and sucrose synthase genes in plant seed development. Gene, 439: 1–10. pmid:19306919
- 44. Nair RB, Baga M, Scoles GJ, Kartha KK, Chibbar RN (1997) Isolation, characterization expression analysis of a starch branching enzyme cDNA from wheat. Plant Sci, 122: 153–163.
- 45. Kondoh H, Lleonart ME, Gil J, Wang J, Degan P, Peters G, et al. (2005) Glycolytic enzymes can modulate cellular life span. Cancer Res, 65: 177–185. pmid:15665293
- 46. Yang WF, Zhang JL, Lü WZ, Cao GC, Chen JY (2014) Study on the differential genes expression in maize embryo treated by a controlled deterioration treatment. Scientia Agricultura Sinica, 47: 1878–1893.
- 47. Dell’Aquila A (1994) Wheat seed ageing and embryo protein degradation. Seed Sci Res 4: 293–298
- 48. Song Y, Zheng Q (2007) Dynamic rheological properties of wheat flour dough and proteins. Trends Food Sci Tech, 18:132–138.
- 49. Kumar KU, Srivastava SP, Kaufman RJ (1999) Double-stranded RNA-activated protein kinase (PKR) is negatively regulated by 60S ribosomal subunit protein L18. Mol Cell Biol, 19: 1116–1125. pmid:9891046
- 50. Goel A, Goel AK, Sheoran IS (2003) Changes in oxidative stress enzymes during artificial ageing in cotton (Gossypium hirsutum L.) seeds. J Plant Physiol, 160: 1093–1100. pmid:14593811
- 51. Weyler C, Bureik M, Heinzle E (2016) Selective oxidation of UDP-glucose to UDP-glucuronic acid using permeabilized Schizosaccharomyces pombe expressing human UDP-glucose 6-dehydrogenase. Biotechnol Lett, 38: 477–481. pmid:26582015
- 52. Qu ZL, Wang HY, Xia GX (2005) GhHb1: a nonsymbiotic hemoglobin gene of cotton responsive to infection by Verticillium dahliae. BBA-Gene Struct Exp, 1730: 103–113.
- 53. Mirelman D, Galum E, Sharon N, Lotan R (1975) Inhibition of fungal growth by wheat germ agglutinin. Nature, 256: 414–416. pmid:1170502
- 54. Kader JC (1996) Lipid-transfer proteins in plants. Annu Rev Plant Biol, 47: 627–654.
- 55. García-Olmedo F, Molina A, Segura A, Moreno M (1995) The defensive role of nonspecific lipid-transfer proteins in plants. Trends Microbiol, 3: 72–74. pmid:7728389
- 56. Cammue BPA, Thevissen K, Hendriks M, Eggermont K, Goderis IJ, Proost P, et al. (1995) A potent antimicrobial protein from onion seeds showing sequence homology to plant lipid transfer proteins. Plant Physiol, 109(2), 445–455. pmid:7480341
- 57. Molina A, Segura A, García-Olmedo F (1993) Lipid transfer proteins (nsLTPs) from barley and maize leaves are potent inhibitors of bacterial and fungal plant pathogens. Febs Lett, 316: 119–122. pmid:8420795
- 58. Carvalho AO, Machado OLT, Da Cunha M, Santos IS, Gomes VM (2001) Antimicrobial peptides and immunolocalization of a LTP in Vigna unguiculata seeds. Plant Physiol Bioch, 39: 137–146.
- 59. Berkowitz O, Jost R, Pollmann S, Masle J (2008) Characterization of TCTP, the translationally controlled tumor protein, from Arabidopsis thaliana. The Plant Cell, 20: 3430–3447. pmid:19060111
- 60. Cao B, Lu Y, Chen G, Lei J (2010) Functional characterization of the translationally controlled tumor protein (TCTP) gene associated with growth and defense response in cabbage. Plant Cell Tiss Org, 103: 217–226.
- 61. Gnanasekar M, Dakshinamoorthy G, Ramaswamy K (2009) Translationally controlled tumor protein is a novel heat shock protein with chaperone-like activity. Biochem Bioph Res Co, 386: 333–337.
- 62. Qin D, Wang F, Geng X, Zhang L, Yao Y, Ni Z, et al. (2015) Overexpression of heat stress-responsive TaMBF1c, a wheat (Triticum aestivum L.) Multiprotein bridging factor, confers heat tolerance in both yeast and rice. Plant Mol Biol, 87: 31–45. pmid:25326264
- 63. Seebauer JR, Moose SP, Fabbri BJ, Crossland LD, Below FE (2004) Amino acid metabolism in maize earshoots. Implications for assimilate preconditioning and nitrogen signaling. Plant Physiol, 136: 4326–4334. pmid:15531710
- 64. Olea F, Pérez-García A, Cantón FR, Rivera ME, Cañas R, Ávila C, et al. (2004) Up-regulation and localization of asparagine synthetase in tomato leaves infected by the bacterial pathogen Pseudomonas syringae. Plant Cell Physiol, 45: 770–780. pmid:15215512
- 65. Lucyshyn D, Busch BL, Abolmaali S, Steiner B, Chandler E, Sanjarian F, et al. (2007) Cloning and characterization of the ribosomal protein L3 (RPL3) gene family from Triticum aestivum. Mol Genet Genomics, 277: 507–517. pmid:17216491
- 66. Miller JD, Arnison PG (1986) Degradation of deoxynivalenol by suspension cultures of the Fusarium Head Blight resistant wheat cultivar Frontana. Can J Plant Pathol 8:147–150.
- 67. Miller JD, Ewen MA (1997) Toxic effects of deoxynivalenol on ribosomes and tissues of the spring wheat cultivars Frontana and Casavant. Nat Toxins 5:234–237 pmid:9615311
- 68. Tzen JT, Huang AH (1992) Surface structure and properties of plant seed oil bodies. The Journal of cell biology, 117: 327–335. pmid:1560029
- 69. Anthony H, Huang C (1992) Oil bodies and oleosins in seeds. Annual Review of Plant Physiology and Plant Molecular Biology, 43: 177–200.
- 70. Tzen JT, Cao YZ, Laurent P, Ratnayake C, Huang AH (1993) Lipids, proteins, and structure of seed oil bodies from diverse species. Plant physiology, 101: 267–276. pmid:12231682
- 71. Bewley J D (1997) Seed germination and dormancy. The Plant Cell, 9: 1055. pmid:12237375
- 72.
Lehninger AL, Nelson DL, Cox MM (2005) Principles of biochemistry. Fourth Edition, Worth Publishers, New York.
- 73. Van den Ende W, Clerens S, Vergauwen R, Van Riet L, Van Laere A, Yoshida M, et al. (2003) Fructan 1-exohydrolases. β-(2, 1)-trimmers during graminan biosynthesis in stems of wheat? Purification, characterization, mass mapping, and cloning of two fructan 1-exohydrolase isoforms. Plant Physiol, 131: 621–631. pmid:12586886
- 74. He AY, Liu CZ, Chen LQ, Ning LJ, Zhang ML, Li EC, et al. (2014) Identification, characterization and nutritional regulation of two isoforms of acyl-coenzyme A oxidase 1 gene in Nile tilapia (Oreochromis niloticus). Gene, 545: 30–35. pmid:24802117
- 75. Liu W, Harrison DK, Chalupska D, Gornicki P, O'Donnell CC, Adkins SW, et al. (2007) Single-site mutations in the carboxyltransferase domain of plastid acetyl-CoA carboxylase confer resistance to grass-specific herbicides. P Natl Acad Sci, 104: 3627–3632.
- 76. Whittle E, Cahoon EB, Subrahmanyam S, Shanklin J (2005) A multifunctional acyl-acyl carrier protein desaturase from Hedera helix L.(English ivy) can synthesize 16-and 18-carbon monoene and diene products. J Biol Chem, 280: 28169–28176. pmid:15939740
- 77. Ohlrogge J, Browse J (1995). Lipid biosynthesis. The Plant Cell, 7: 957. pmid:7640528
- 78. Chan DI, Vogel HJ (2010) Current understanding of fatty acid biosynthesis and the acyl carrier protein. Biochem J, 430: 1–19. pmid:20662770
- 79. Nguyen C, Haushalter RW, Lee DJ, Markwick PR, Bruegger J, Caldara-Festin G, et al. (2014) Trapping the dynamic acyl carrier protein in fatty acid biosynthesis. Nature, 505: 427–431. pmid:24362570
- 80. Macheroux P, Schmid J, Amrhein N, Schaller A (1999) A unique reaction in a common pathway: mechanism and function of chorismate synthase in the shikimate pathway. Planta, 207: 325–334. pmid:9951731
- 81. Bentley R, Haslam E (1990) The shikimate pathway-a metabolic tree with many branches. Crit. Rev. Biochem Mol Biol, 25: 307–384. pmid:2279393
- 82. Han J, Tian YS, Xu J, Wang LJ, Wang B, Peng RH, et al. (2014) Functional characterization of aroA from Rhizobium leguminosarum with significant glyphosate tolerance in transgenic Arabidopsis. J Microbiol Biotechn, 24: 1162–1169.
- 83. Wong YS, Mazelis M (1981) Detection and properties of L-histidinol dehydrogenase in wheat germ. Phytochemistry, 20: 1831–1834.
- 84. Feldman-Salit A, Wirtz M, Lenherr ED, Throm C, Hothorn M, Scheffzek K, et al. (2012) Allosterically gated enzyme dynamics in the cysteine synthase complex regulate cysteine biosynthesis in Arabidopsis thaliana. Structure, 20: 292–302. pmid:22325778
- 85. Wang ZQ, Yuan YZ, Ou JQ, Lin QH, Zhang CF (2007) Glutamine synthetase and glutamate dehydrogenase contribute differentially to proline accumulation in leaves of wheat (Triticum aestivum) seedlings exposed to different salinity. J Plant Physiol, 164: 695–701. pmid:16777263
- 86. Dupont FM (2008) Metabolic pathways of the wheat (Triticum aestivum) endosperm amyloplast revealed by proteomics. BMC Plant Biol, 8: 1.
- 87. Martineau Y, Derry MC, Wang X, Yanagiya A, Berlanga JJ, Shyu AB, et al. (2008) Poly (A)-binding protein-interacting protein 1 binds to eukaryotic translation initiation factor 3 to stimulate translation. Mol Cell Biol, 28: 6658–6667. pmid:18725400
- 88. Lin XQ, Liang SL, Han SY, Zheng SP, Ye YR, Lin Y (2013) Quantitative iTRAQ LC-MS/MS proteomics reveals the cellular response to heterologous protein overexpression and the regulation of HAC1 in Pichia pastoris. J Proteomics, 91: 58–72. pmid:23851310
- 89. Kimura S, Higashino Y, Kitao Y, Masuda T, Urade R (2015) Expression and characterization of protein disulfide isomerase family proteins in bread wheat. BMC Plant Biol, 15: 1.
- 90. Picard D (2002) Heat-shock protein 90, a chaperone for folding and regulation. Cell Mol Life Sci, 59: 1640–1648. pmid:12475174
- 91. O'Sullivan NC, Jahn TR, Reid E, O'Kane CJ (2012) Reticulon-like-1, the Drosophila orthologue of the hereditary spastic paraplegia gene reticulon 2, is required for organization of endoplasmic reticulum and of distal motor axons. Hum Mol Genet, 21: 3356–3365. pmid:22543973
- 92. Permyakova MD, Trufanov VA, Pshenichnikova TA, Ermakova MF (2010) Role of lipoxygenase in the determination of wheat grain quality. Applied Biochem Micro, 46: 87–92.
- 93.
Trewavas A (2000) Signal perception and transduction. Pages 930–987 in: Biochemistry and Molecular Biology of Plants. Buchanan B. B., Gruissem W., and Jones R. L., eds. Am. Soc. Plant Physiologists: Rockville, MD.
- 94. Liu HT, Li B, Shang ZL, Li XZ, Mu RL, Sun DY, et al. (2003) Calmodulin is involved in heat shock signal transduction in wheat. Plant Physiol, 132:1186–1195. pmid:12857801
- 95. Christov NK, Imai R, Blume Y (2008) Differential expression of two winter wheat alpha-tubulin genes during cold acclimation. Cell Biol Int, 32: 574–578. pmid:18162419
- 96. Ridha Farajalla M, Gulick PJ (2007) The α-tubulin gene family in wheat (Triticum aestivum L.) and differential gene expression during cold acclimation. Genome, 50: 502–510. pmid:17612619
- 97. Mathur M, Saluja D, Sachar RC (1991) Post-transcriptional regulation of S-adenosylmethionine synthetase from its stored mRNA in germinated wheat embryos. BBA-Protein Struct Mol M, 1078: 161–170.
- 98. Rebeille F, Neuburger M, Douce R (1994) Interaction between glycine decarboxylase, serine hydroxymethyltransferase and tetrahydrofolate polyglutamates in pea leaf mitochondria. Biochem J, 302: 223–228. pmid:7520695
- 99. Hohm T, Demarsy E, Quan C, Petrolati LA, Preuten T, Vernoux T, et al. (2014) Plasma membrane H+-ATPase regulation is required for auxin gradient formation preceding phototropic growth. Mol Syst Biol, 10:751. pmid:25261457
- 100. Uematsu K, Suzuki N, Iwamae T, Inui M, Yukawa H (2012) Increased fructose 1, 6-bisphosphate aldolase in plastids enhances growth and photosynthesis of tobacco plants. J Exp Bot, 63: 3001–3009. pmid:22323273
- 101. Neal AL, Geraki K, Borg S, Quinn P, Mosselmans JF, Brinch-Pedersen H, et al. (2013) Iron and zinc complexation in wild-type and ferritin-expressing wheat grain: implications for mineral transport into developing grain. J Biol Inorg Chem, 18: 557–570. pmid:23624806
- 102. McCubbin WD, Kay CM, Lane BG (1985) Hydrodynamic and optical properties of the wheat germ Em protein. Can J Biochem Cell Biol, 63: 803–811.
- 103. Kibinza S, Bazin J, Bailly C, Farrant JM, Corbineau F, El-Maarouf-Bouteau H (2011) Catalase is a key enzyme in seed recovery from ageing during priming. Plant Sci, 181: 309–315. pmid:21763542
- 104. Li X, Zhang JB, Song B, Li HP, Xu HQ, Qu B, et al. (2010) Resistance to Fusarium head blight and seedling blight in wheat is associated with activation of a cytochrome P450 gene. Phytopathology, 100:183–191. pmid:20055652
- 105. Yan X, Que Y, Wang H, Wang C, Li Y, Yue X, et al. (2013) The MET13 methylenetetrahydrofolate reductase gene is essential for infection-related morphogenesis in the rice blast fungus Magnaporthe oryzae. PloS one, 8(10), e76914. pmid:24116181
- 106. Christova PK, Christov NK, Imai R (2006) A cold inducible multidomain cystatin from winter wheat inhibits growth of the snow mold fungus, Microdochium nivale. Planta, 223: 1207–1218. pmid:16320069
- 107. Caruso C, Caporale C, Chilosi G, Vacca F, Bertini L, Magro P, et al. (1996) Structural and antifungal properties of a pathogenesis-related protein from wheat kernel. J Protein Chem, 15: 35–44. pmid:8838588
- 108. Khlestkina E (2013) The adaptive role of flavonoids: emphasis on cereals. Cereal Res Commun, 41: 185–198.
- 109. Liu Z, Liu Y, Pu Z, Wang J, Zheng Y, Li Y, et al. (2013) Regulation, evolution, and functionality of flavonoids in cereal crops. Biotechnol Lett, 35: 1765–1780. pmid:23881316
- 110. Shoeva OY, Khlestkina EK, Berges H, Salina EA (2014) The homoeologous genes encoding chalcone-flavanone isomerase in Triticum aestivum L.: Structural characterization and expression in different parts of wheat plant. Gene, 538:334–341. pmid:24480448
- 111. Lee SB, Jung SJ, Go YS, Kim HU, Kim JK, Cho HJ, et al. (2009) Two Arabidopsis 3-ketoacyl CoA synthase genes, KCS20 and KCS2/DAISY, are functionally redundant in cuticular wax and root suberin biosynthesis, but differentially controlled by osmotic stress. The Plant J, 60:462–475. pmid:19619160
- 112. Vizcaíno JA, Csordas A, del-Toro N, Dianes JA, Griss J, Lavidas I, et al. (2016) 2016 update of the PRIDE database and related tools. Nucleic Acids Res 44(D1): D447–D456. pmid:26527722