Introduction
Wheat is one of the most important food crops in the world. Low-temperature injury during winter and spring can be particularly destructive on wheat production. The winter hardiness of modern cultivars and good management practices have reduced winterkilling of wheat, so injury during winter is not a big problem nowadays. While, with the frequent occurrence of climatic anomalies, SRF has become a significant limiting factor on wheat production. SRF injury frequently occurs when low temperatures coincide with sensitive plant growth stages [
1-
3], from pistil and stamen differentiation stage to anther differentiation stage [
4-
7]. The research into cold tolerance at the vegetative stages of wheat development has revealed that genetic variation does exist which was controlled by few genetic loci with major effects, such as
VRN, Fr and
CBF transcription factors [
8-
12]. To date, there have been few studies on the underlying mechanisms involved in the SRF response at the reproductive stages of common wheat [
13-
15].
Transcriptome profiling has contributed markedly to our understanding of cold stress in species including Arabidopsis and rice. However, it has some limitations because mRNA levels are not always correlated to those of corresponding proteins, due in part to post-transcriptional regulation. Recently, proteomic analysis has become a preferred strategy to reveal the dynamics of expression under cold stress.
Protein synthesis and accumulation in developing wheat grains has been extensively investigated [
16-
18], with most studies focusing on plant proteomic analysis under abiotic stress conditions [
19,
20]. Also comparative proteomic analysis of wheat seedling roots, leaves, and grains, as well as embryos, between cultivars differing in tolerance to drought and salt has been reported [
21-
23]. Some work has been conducted to uncover proteomic changes in wheat under cold stress. Danyluk et al. [
24] found that 18 proteins were transiently induced in the three wheat cultivars after cold treatment. There are 53 kinds of proteins sustaining high levels of expression in the induction of freezing tolerance required for four weeks, in which 34 kinds of proteins of the best freezing tolerance cultivars also had a higher level of expression. Sara et al. [
25] analyzed the wheat plants grown at 20 or 4°C until entering the reproductive stage and a cross-comparison on the leaf proteomes experiments, which found a series of proteins related to cold tolerance. However, most of the above mentioned studies were focused on the cold responding at the vegetative stage.
Our previous studies revealed that Jimai22 has almost the same development process and frost resistance with another wheat cultivar Luyuan301 at the vegetative stage, but in the event of SRF, Jimai22, at pistil and stamen differentiation stage, was less affected, while Luyuan301 was injured to bald seriously. In the present study, we performed comparative proteome analyses using two-dimensional electrophoresis (2-DE) and MALDI-TOF MS. The objectives were to determine the changes of protein expression in metabolic pathways of wheat (T. aestivum) after cold treatment at the reproductive stage and improve our understanding the effects of SRF on plant growth and metabolism, which will provide valuable information for the breeding commercial cultivars with SRF tolerant.
Materials and methods
Plant materials and stress treatment
Experiments were carried out on two wheat cultivars, Jimai22 and Luyuan301. Seeds were germinated and grown in at 2°C for vernalization. After 40 days, the seedlings were transplanted to greenhouse and grown under 22/16°C day/night cycle. Forty days after transplanting, most of the young spikes had developed to pistil and stamen differentiation stage, and then the seedlings were divided into two groups. One was used as the control sample at 22 /16°C in green house, while the other was moved to a chilling chamber for frost treatment at 0°C for 48 h. The young spikes and leaves were harvested and mixed together, stored in liquid nitrogen for protein extraction.
Physiologic analysis
The fresh weights (FW) of second leaves were measured immediately after harvesting. The leaves were floated on deionized water for 24 h at 4°C under low irradiance. The turgid leaves were then quickly weighed (TW) and dry weights were determined after oven-drying at 70°C for 48 h. The relative water content (RWC) was calculated as: RWC/% = (FW-DW)/(TW-DW) × 100 [
26].
Total soluble sugars of wheat leaves (about 100 mg) were extracted in boiling water for 30 min and determined by anthrone reagent using glucose as standard. The absorbance was read at 630 nm and the sugar concentration was determined from a glucose standard curve [
27].
Protein extraction and quantification
Soluble proteins from leaf tissues were homogenized separately to a fine powder in liquid nitrogen, and were transferred into 50 mL tubes. The pellets were added to 25 mL tricholoroacetic acid/acetone (1:9), 65 mmol·L-1 DTT and incubated at -20°C for 1 h and centrifuged for 45 min at 12000 r·min-1. The supernatant was removed, 25 mL pure acetone was added and incubated at -20°C overnight. This step was repeated three times and the pellets were lyophilized. Powdered pellet (200 mg) was dissolved in 0.5 mL 2D protein extraction buffer and ultrasonic for 3 min in an ice bath. After centrifugation for 45 min at 12000 rpm, the supernatant was filtered (0.22 μm). Finally, the protein concentration was determined with a 2-D Quant Kit (Sangon, Shanghai, China) using BSA (2 mg·mL-1) as the standard following the manufacturer’s instructions. The protein samples were kept at -80°C until used.
2-DE and image analysis
Protein samples of 80 μg were loaded onto an Immobiline DryStrip (GE, pH 4-7) and were then placed on an IPG strip to perform IEF. Isoelectric focusing was conducted in five steps: 30 V for 12 h; 500 V for 1 h; 1000 V for 1 h; 8000 V for 8 h; and 500 V for 4 h. The focused strips were then treated twice in equilibration solution for 15 min. For the second dimension, the proteins were separated on 12.5% SDS-polyacrylamide gels. After SDS-PAGE, gels were stained with silver staining. The images of the stained gels were acquired using a GS-800 calibrated densitometer (Bio-Rad) and analyzed using PDQuest software, version 7.2.0 (Bio-Rad). After background subtraction and spot detection, spots were matched and normalized using the method of total density in the gel image. Student’s t-tests (n = 3, P<0.05) were performed to determine if the proteins differentially accumulated.
In-gel digestion and MS analysis
Protein spots showing at least a 2-fold difference (P<0.05) in abundance between control and treatment were excised manually into 1.5 mL microtubes. Gel pieces were first washed with 200 μL of bleaching solution (100 mmol·L-1 NH4HCO3/30% acetonitrile), repeated 2-3 times until colorless. The dry particles were digested with 5 μL of 10 ng·μL-1 trypsin enzyme (Promega, USA), and incubated at 37°C for 20 h.
Tryptic peptides were analyzed with a 4800PlusMALDI TOF/TOFTM Analyzer (Applied Biosystems, USA). The spots were analyzed with parameters of laser intensity 5200, acceleration voltage 8 kV, grid voltage 91.25% and delay time 250 ns, UV nitrogen laser 355 nm, and 40-50 times of single scan of MS signal accumulation. The MS/MS ion search was performed using MASCOT (http://www.matrixscience.com) searching the NCBI (National Center for Biotechnology Information) non-redunant database. The parameters for searching were: monoisotoic peptide masses, ± 0.2 Da peptide mass tolerance; one missed cleavage site, enzyme of trypsin, modifications allowed for oxidation of methionine and carboxy-amidomethylation of cysteine. Only significant hits as defined by the MASCOT probability analysis (P<0.05) were accepted. BLAST homology and similarity searches were conducted with the protein Viridiplantae index of the non-redundant NCBI database.
Results
Physiologic responses induced by low temperature stress
The relative water content assay (RWC) was measured and results showed that the RWC of the two cultivars were both decreased after the cold treatment (Fig. 1a). This indicated that the leaves of wheat were subjected to chilling injury. RWC of Jimai22 was higher than Luyuan301, which did not the significant differences.
Soluble sugars content was detected. Both of the two wheat cultivars were relatively low and at similar levels as shown in Fig. 1b. After the stress treatment, both the sugar increased, and the contents showed significant differences between Jimai22 and Luyuan301 plants.
2-DE analysis of cold-responsive total proteins in Jimai 22 and Luyuan 301
2-DE analysis was carried out and a representative gel is shown in Fig. 2. Approximately 1100-1400 protein spots were detected on silver-stained gels in both Jimai 22 and Luyuan 301, and proteins were better separated in both directions of the isoelectric and molecular weight. We detected approximately 1200 proteins in total and 24 (2%) of them were considered to have changed significantly between cold-treated and untreated plants. The results showed that the occurrence of cold-treated resulted in both up- and downregulation of a series of genes, which may relate to the morphological and physiologic changes of cold-treated plants.
Identification of the differentially expressed proteins (DEPs) by MALDI-TOF MS
A total of 49 protein spots showed at least 2-fold differences in abundance between the two cultivars in response to cold stress. Forty-three DEPs from the 49 protein spots between cold-treated and untreated plants were identified by MALDI-TOF MS. The results revealed that 21 proteins in Jimai22 differed in the protein map, including 13 upregulated and 8 downregulated (Table 1); 22 proteins in Luyuan301 differed including 9 upregulated and 13 downregulated (Table 2). These DEPs were mainly involved in six biologic processes, carbohydrate metabolism, amino acid metabolism, resistance proteins and antioxidant enzymes, photosynthesis and cellular respiration proteins, cell-wall related proteins, protein translation/processing/degradation and signal transduction. The protein levels indicated that wheat cold tolerance was involved in the changes in cell metabolism and signal or biochemical pathways. The peptide sequence data for the 43 different expression proteins identified by MS was listed in Table S1.
Discussion
Carbon metabolism proteins
A number of proteins important in sugar metabolic pathways had the same expression pattern in cold-treated plants compared to the untreated plants. Fructose-bisphosphate aldolase (spot 682 and 689), an important enzyme associated both in glycolysis and pentose phosphate pathway, which could provide ATP and substrate that is essential for the organism. The upregulated of this enzyme implied the sugar metabolism alteration in Jimai22, while in Luyuan301 was downregulated (spot 428) compared to their untreated plants [
28].
RuBisCO (ribulose-1,5-bisphosphate carboxylase/oxygenase) is a bifunctional enzyme that catalyzes both the carboxylation and oxygenation of ribulose-1,5-bisphosphate (RuBP), thus fixing carbon dioxide as the first step of the Calvin cycle. Compared to untreated plants, the relative contents of Rubisco proteins (spots 76, 142, and 366) and Rubisco activase1 (spot 444) in the total proteins of Jimai22 in cold-treated plants were lower than untreated plants, which suggested that the capability of carbon fixation and assimilation might be relatively decreased in cold-treated plants. And in Luyuan301 plants, Rubisco proteins (spot 667) and Rubisco activase1 (spot 640) were upregulated, while Rubisco activase (spot 394) was downregulated compared to untreated plants.
Malate dehydrogenase (MDH) (spot 375) was commonly found in animals, plants, bacteria and other organisms, is a key enzyme in the glucose metabolism, which catalyzes the reversible conversion between malic acid and oxaloacetate.
Differential accumulation of these proteins implied a shift of sugar metabolic direction after cold stress and then resulted in abundance changes of components in sugar pool, which is likely to be an important reason for drastic morphological and physiologic modification caused by cold.
Amino acid metabolism proteins
The abundance of some enzymes involved in amino acid metabolism showed significant differences between cold-treated and untreated plants. Cysteine synthase (spot 316) was significantly accumulated in cold-treated plants. The result indicated that the effect of cold stress was also reflected in a dynamic modification of both amino acid catabolism and anabolism, and also suggested the specificity for synthesis and/or degradation of amino acids after cold treatment.
Resistance proteins and antioxidant enzymes
Plant resistance to cold stress involves a series of inducible defensive responses, which include the synthesis of various cold-related proteins and multiple enzymes with antioxidant activity.
Quinone oxidoreductase-like protein (spot 692) belongs to zinc-containing alcohol dehydrogenase family and quinone oxidoreductase subfamily. It has recently been demonstrated that purified NAD(P)H:quinone oxidoreductase 1 (NQO1) is able to scavenge superoxide (O2(.-)) though the rate of reaction of O2(.-) with NQO1 is much lower than the rate of enzymatic dismutation catalyzed by superoxide dismutase (SOD). This gene was confirmed to play an important role as the force generators for fish in a low water environment. In addition, NADH-quinone oxidoreductase 3 subunit gene was associated with cold adaptation of common carp to some extent in oxidation of respiratory chain electron transport system [
29].
Vitamin B6 is an essential coenzyme for numerous metabolic enzymes and is a potent antioxidant. The de novo pathway of vitamin B6 biosynthesis has only been described recently and involves the protein Pyridoxal phosphate synthase protein (PDX1, spot 345).
Cold-shock proteins (CSPs) (spot 149) are nucleic acid binding proteins containing a cold shock domain (CSD), which are highly conserved in bacteria, plants and animals. They play a critical role for transcription, translation, development and stress responses by the RNA chaperones function.
Photosynthesis and cellular respiration proteins
A light-harvesting complex (spot 127) is a complex of subunit proteins that may be part of a larger supercomplex of a photosystem, the functional unit in photosynthesis. The increase of LCHI aggregates in the study may suggest that LCHI aggregates may represent a protective mechanism to cold-stress plants.
ATPases (spot 455 and 447) are membrane-bound enzyme complexes/ion transporters that combine ATP synthesis and/or hydrolysis with the transport of protons across a membrane. ATPases can harness the energy from a proton gradient, using the flux of ions across the membrane via the ATPase proton channel to drive the synthesis of ATP. The plasma membrane proton pump ATPase [
30] is important in transport across the plasma membrane. As a primary transporter, it mediates ATP-dependent H(+) extrusion to the extracelllar space, thus creating pH and potential differences across the plasma membrane that activate a large set of secondary transporters. ATPases (spot 455 and 447) of Jimai22 were downregulated while spot 455 of Luyuan301 was upregulated when they compared to non-stressed plants.
Cell-wall related proteins
Germin and Germin-like protein have various proposed roles in plant developmental stages and stress-related processes. It may function through the controlling of the cell wall extensibility in the early plants development, but their expression is also regulated by pathogens and abiotic stress factors such as photoperiod [
31].
Glycine-rich RNA binding proteins (GRPs; spot 109) have been found in the cell walls of many higher plants and form a group of structural protein components of the wall in addition to extensins and proline-rich proteins. GRPs are involved in post-transcriptional regulation of gene expression in plants under various conditions and, in most cases, they are accumulated in the vascular tissues and their synthesis is part of defense mechanisms in plants [
32].
Protein translation, processing and degradation
Regulation of gene expression is achieved at several levels; transcriptional, post-transcriptional, translational, and post-translational. Pentatricopeptide repeat proteins (PPR; spot 138) are encoded by a large gene family, which are involved in plant development, organelle biogenesis, restoring of cytoplasmic male sterilities, RNA processing and editing in mitochondrial and chloroplasts, plastid-to-nucleus retrograde signaling, and responses to environmental stresses. PPR proteins have been implicated in RNA editing, RNA processing, RNA splicing and translational activation. Many miRNAs were processed and edited from PPR proteins, which may be directly or indirectly involved in plant resistance to environmental stresses.
Changes of translation activity were considered to be primary stress responses [
33]. Spots 57 and 66 were 50S ribosomal protein which involved in the protein translation and synthesis.
Eukaryotic translation initiation factor (eIF-5A) (spot 84) is important in the regulation of biologic growth, aging and environmental responses. It’s upregulation can enhance the plants stress tolerance.
Heat shock proteins (spot 49), also known as stress proteins, are a group of highly structurally conserved peptide, which are widely found in prokaryotic, eukaryotic cells and are involved in cell damage and repair.
Signal transduction
Ras GTP binding protein (spot 653) is an evolutionary conserved eukaryotic GTPase, which is likely to be involved in nuclear translocation of proteins and cell cycle progression [
34]. However, little is known about the function of Ras in plant stress response.
Conclusions
Cold tolerance is a complex phenomenon in most plants, and involves numerous mechanisms at the cellular, tissue, organ and whole plant levels. In the present study, a comparative proteomic analysis of wheat subjected to cold treatment for 48 h at the pistil and stamen differentiation stage was performed. The physiologic results showed the relative higher RWC and soluble sugars may have increased the cold stress tolerance of Jimai22 compared to Luyuan301. Forty-three different cold-responsive proteins were identified (Table 1 and Table 2). The proteins identified were implicated in six physiologic process; including carbohydrate metabolism, amino acid metabolism, resistance proteins and antioxidant enzymes, photosynthesis and cellular respiration proteins, cell-wall related proteins, protein translation/processing/degradation and signal transduction.
The same proteins involved in the two wheat cultivars may had different expression mode, such as Fructose-bisphosphate aldolase, Rubisco proteins and ATPase. Other proteins were cultivar specific and involved in many different biochemistry pathways, for example, quinone oxidoreductase-like protein, cold-shock protein, germ-like protein and 50S ribosome protein. The different expression of these proteins in various pathways may cause the different cold tolerance for the two wheat cultivars. The unknown proteins may represent novel mechanisms for cold tolerance in wheat. Further study is needed to reveal their structures and functions.
The proteins identified in this study represent only a small part of the wheat proteome responsive to cold treatment, and may other cold-responsive proteins still need to be identified. The identification of cold-responsive proteins provides not only new insights to cold stress responses but also a good starting point for further dissection of their functions using genetic and approaches.
Higher Education Press and Springer-Verlag Berlin Heidelberg