Proteome analysis of soybean leaves, hypocotyls and roots under salt stress

Background Salinity is one of the most widespread agricultural problems in arid and semi-arid regions that makes fields unproductive, and soil salinization is a serious problem in the entire world. To determine the effects of salt stress on soybean seedlings, a proteomic technique was used. Results Soybean plants were exposed to 0, 20, 40, or 80 mM NaCl for one week. The effect of treatment at 20 mM NaCl on plant growth was not severe, at 80 mM NaCl was lethal, and at 40 mM NaCl was significant but not lethal. Based on these results, proteins were extracted from the leaves, hypocotyls and roots of soybean treated with 40 mM NaCl. Nineteen, 22 and 14 proteins out of 340, 330 and 235 proteins in the leaves, hypocotyls and roots, respectively, were up- and down-regulated by NaCl treatment. In leaves, hypocotyls and roots, metabolism related proteins were mainly down-regulated with NaCl treatment. Glyceraldehyde-3-phosphate dehydrogenase was down-regulated in the leaf/hypocotyls, and fructokinase 2 was down-regulated in the hypocotyls/root with NaCl treatment. Stem 31 kDa glycoprotein precursor was up-regulated in all three organs with NaCl treatment. Glyceraldehyde-3-phosphate dehydrogenase was specifically down-regulated at the RNA and protein levels by salt stress. Conclusion These results suggest that metabolism related proteins play a role in each organ in the adaptation to saline conditions.


Background
Soybean is an important dicot crop due to the high content of oil and protein in its seeds [1]. However, soybean is subject to abiotic stresses that reduce its yield like many crops. Salinity is one of the most widespread agricultural problems in arid and semi-arid regions that makes fields unproductive, and soil salinization is a serious problem in the entire world [2]. Salt stress severely inhibits plant growth for two reasons: first by an osmotic or water-deficit effect of salinity and second by a salt-specific or ion-excess effect of NaCl. However, plants suffer from composite stresses caused by salinity, including water deficit and ion imbalance [3]. Adaptation to salt stress requires alterations in gene expression and subsequently the protein profile of the plant and is very complicated at the whole plant and cellular levels [4,5].
Some salt-inducible genes have been investigated in soybean. A homologue of oxysterol-binding protein was involved in the salt-stress response and cotyledon senescence of soybean [6]. An acidic isoform of pathogenesisrelated protein group 5 (PR-5) that is responsive to high salt stress and dehydration was located in the extracellular space of soybean roots [7]. A leucine-zipper-like protein was induced under salt stress and acted in mature organs of soybean shoots to counteract water-potential changes [8]. An acid phosphatase was related to the adaptation of soybean to salt stress, and was involved in reactive oxygen species formation or scavenging or in stress-responsive signal transduction pathways [9]. The overexpression of a dehydration responsive element binding protein homologous gene (GmDREB2) in soybean caused the accumulation of a higher level of free proline compared to wild-type plants under salt stress; this gene also was an important transcriptional activator and was useful in improving plant tolerance to salt stress [10]. These salt stress-induced genes may lead to up-regulation or down-regulation of salt stress related proteins [3].
Abbasi and Komatsu [11] studied salt responsive proteins in rice using a proteomic technique, which indicated that an oxygen evolving enhancer protein expressed in the leaf sheath and leaf blade of rice showed a coordinated response to salt stress. Proteomics is a powerful molecular tool for describing the complete proteome at the organelle, cell, organ, or tissue level and for comparing how the proteome is affected by different physiological conditions. Two-dimensional polyacrylamide gel electrophoresis (2-DE) is one of the most sensitive and powerful techniques for resolving hundreds of proteins [12]. 2-DE has been applied to different abiotic treatments for soybean including treatment with salt [13], flooding [14] and exposure to ultraviolet radiation [15]. Proteomic responses of citrus [16], Salicornia europaea [17], Bruguiera gymnorhiza [18], Suaeda Aegyptiaca [19] and rice [20] to salt stress have been reported.
Only one proteomics study has been reported on salt stress in soybean. Aghaei et al. [11] identified salt responsive proteins using a proteomic technique concluding that especially late embryogenesis-abundant protein was involved in the process of adaptation to salt stress at the early seedling stage. Aghaei et al. [11] used 3-day-old seedlings and 100 mM NaCl and proteins were extracted from combined roots and hypocotyls. Soybean seedlings suffer from NaCl and it is very important to improve salt tolerance on seedlings for transition from vegetative to reproductive stage on the farm [1,21]. In this study, proteome analysis was performed on leaf, hypocotyl and root of 7-day-old seedlings to determine the importance of salt-responsive proteins in vegetative stage. Proteins from these organs were separated by 2-DE, and the major differentially expressed proteins were identified by protein sequencer and mass spectrometry (MS).

Results
The effect of 40 mM NaCl on soybean seedlings is significant but not lethal To evaluate the effects of salt stress, soybean seeds were sown and treated with 0, 20, 40 or 80 mM NaCl. Morphological changes were examined ( Figure 1A) and the length and fresh weight of the hypocotyl and root of 7-day-old seedlings were measured ( Figure 1B). After 7 days, first and primary leaves had opened in control plants, whereas they remained closed in treatment with 20 and 40 mM NaCl. There was no first leaf in plants treated with 80 mM NaCl ( Figure 1A). The length of the hypocotyl and root was shorter in NaCl treatment. The length of hypocotyl and root of plant treated with 80 mM NaCl was 60% and 28% shorter, respectively, than control plants. The shorter length of hypocotyl and root was significant for 40 mM NaCl compared to the control plants. The fresh weight of hypocotyls and roots from plants treated with 80 mM NaCl was 42% and 56% of the control plants, respectively. The fresh weight of hypocotyls and roots was significantly lower at 40 mM NaCl compared to the control plants. The effect of treatment at 20 mM NaCl on plant growth was not severe even after 7 days, the effect of treatment at 40 mM NaCl was significant but not lethal, and the effect of treatment at 80 mM NaCl was lethal ( Figure 1B).
Na and K contents in the leaf, hypocotyl and root of soybeans treated with NaCl To investigate the effect of salt stress on the K and Na contents in soybean, the concentrations of these ions were measured in root, hypocotyls and leaf of soybean at 0, 20, 40 and 80 mM NaCl after 2 weeks (see Additional file 1). At 40 mM NaCl, the Na content of the leaf, hypocotyle and root increased, but K content did not changed under NaCl treatment (see Additional file 1).
In the leaf, hypocotyl and root, 19, 22 and 14 proteins are changed, respectively, in response to 40 mM NaCl treatment Soybean seeds were sown and treated with 40 mM NaCl, then proteins were extracted from the leaf, hypocotyl and root, separated by 2-DE, and stained with CBB to evaluate their expression level. Using PDQuest software analysis, 340, 330 and 235 protein spots were reproducibly detected on the 2-DE gels of leaf, hypocotyl and root, respectively (Figures 2, 3 and 4, upper panels), with 11, 12 and 7 proteins visibly up-regulated and 8, 10 and 7 proteins down-regulated by 40 mM NaCl treatment (Figures 2, 3 and 4, lower panels, see Additional files 2, 3 and 4). Aghaei et al. [13] reported that 321 protein spots were detectable from 2-DE gels for 3-day-old soybean seedlings. Among these spots, 7 proteins were up-or down-regulated under salt stress.

Metabolic related proteins are down-regulated in soybean roots by NaCl treatment
In the root, 13 proteins were identified by MALDI-TOF MS (see Additional files 4 and 7). These proteins were putative fructokinase 2 (R01), dienelactone hydrolase family protein (R02), caffeoyl-CoA-O-methyltransferase (R06), NADPH:isoflavone reductase (R07), putative cinnamyl alcohol dehydrogenase (R09), putative quinone oxidoreductase (R10), stem 31 kDa glycoprotein precursor (R11) and ripening related protein (R13). In the root, proteins involved in defense (2 proteins) were up-regulated, whereas proteins involved in plant metabolism (5 proteins) were down-regulated by NaCl treatment ( Figure 5). Figure 5 Functional distribution of identified proteins expressed in leaf, hypocotyl and root of soybean seedlings after NaCl treatment. After identification of proteins up-(black columns) and down-(white columns) regulated with NaCl treatment, the proteins were classified into 9 functional groups, which are: metabolism, photosynthesis, defense, cell transfer, protein destination and storage, protein synthesis, transcription, not hit and not determined. Proteins in common among leaves, hypocotyls and roots show organ and developmental stage specificity at the mRNA level A total of 19, 22 and 14 proteins were changed by 40 mM NaCl treatment in the leaf, hypocotyl and root, respectively. Three of these proteins were involved in at least two of these organs with salt stress. These proteins were glyceraldehyde-3-phosphate dehydrogenase (H15, L04), stem 31 kDa glycoprotein precursor (H13, L04, R11) and putative fructokinase 2 (H09, R01) ( Figure 6). The remaining proteins were organ specific, suggesting that salt stress affects organ specific proteins by up-or down-regulation. These results indicate that three proteins are mostly affected by NaCl treatment. To identify whether these genes coding proteins are expressed with organ specificity, RNA expression levels were examined. Total RNA was extracted from 3-and 5-day-old hypocotyls and from 7-day-old leaves, hypocotyls and roots of soybean seedlings and quantitative real-time PCR was then performed.
The relative amount of glyceraldehyde-3-phosphate dehydrogenase mRNA was decreased in the leaves, hypocotyls and roots relative to controls by NaCl treatment. The decrease was significant for hypocotyls and roots ( Figure 7A). The expression of this protein was down-regulated in leaves and hypocotyls (see Additional files 5 and 6). These results indicate that the down-regulation of glyceraldehyde-3-phosphate dehydrogenase is caused at both the mRNA and protein level by NaCl treatment. On the other hand, fructokinase 2 was differentially expressed in the hypocotyl and root under salt stress ( Figure 6). The relative amount of putative fructokinase 2 mRNA was increased in the leaf, hypocotyl and root, but not significantly ( Figure 7C). At the protein level, this protein was down-regulated in the hypocotyl and root (see Additional files 6 and 7). Fructokinase catalyzes the phosphorylation of fructose to fructose-6-phosphate. Fructose-6-phosphate is a major substrate for many sugar metabolic pathways including glycolysis, starch biosynthesis, and the oxidative pentose pathway [22]. Under salt stress, by inhibiting photosynthesis, the amount of fructose-6-phosphate decreases, because of the decrease in photosynthetic sucrose. Fructokinase was down-regulated in this study, presumably as a result of negative feedback caused by a decrease in the amount of fructose-6-phosphate.
The relative amount of stem 31 kDa glycoprotein precursor mRNA was increased in the hypocotyl and root but decreased in the leaf by NaCl treatment. The changes for stem 31 kDa glycoprotein precursor mRNA were not significant among the three organs ( Figure 7B). The expression of this protein was up-regulated in all three organs (see Additional files 5, 6 and 7). Aghaei et al. [13] reported that the expression of stem 31 kDa Figure 6 Venn diagram analysis showing up-or down-regulated proteins that overlapped between the leaf, hypocotyl and root of soybean seedlings. Effects of salt stress on these three parts of soybean seedlings were analyzed. Numbers correspond to the protein spots in the 2-DE patterns from the leaf, hypocotyl and root. glycoprotein precursor was down-regulated in hypocotyls and roots of 3-day-old soybean seedlings by NaCl treatment. In this study, the amount of stem 31 kDa glycoprotein precursor mRNA was increased in the hypocotyl of 3-day-old to 5-day-old soybean seedlings ( Figure  7B). Based on this study, after 3 days there is a developmental increase in the amount of stem 31 kDa glycoprotein precursor mRNA. Stem 31 kDa glycoprotein precursor may function as a somatic storage protein during early seedling development and it accumulates mainly in the stems of developing soybean seedlings [23,24]. Its differential change in different organs of soybean seedlings under salt stress indicates its rate of consumption from storage and de novo biosynthesis during development. At the protein level, this protein was upregulated in all three organs examined, which indicates the decrease in consumption of this storage protein caused by the decrease in plant growth under salt stress. Figure 7 Relative amount of mRNA for salt responsive genes in leaf, hypocotyl and root of soybean seedlings. Total RNA was extracted from 3-and 5-day-old hypocotyls and from 7-day-old leaves, hypocotyls and roots of soybean seedlings after 40 mM NaCl treatment. cDNA was synthesized by reverse transcription, and quantitative real-time PCR was performed using the cDNA products. An 18S rRNA cDNA clone was used as a template to produce a standard curve and the relative amount of mRNA was calculated by the ΔCt method. Values are the mean ± SE from 3 independent experiments and asterisks indicate significant (P < 0.05) differences between control (white columns) and NaCl treatments (black columns). GAPDH, glyceraldehyde-3-phosphate dehydrogenase; S31, stem 31 kDa glycoprotein precursor; FK, putative fructokinase 2; L, leaf; H, hypocotyl; R, root; H3, 3-day-old hypocotyl and H5, 5-day-old hypocotyl.

Discussion
Growth in length and fresh weight of the hypocotyls and roots was completely abolished by 80 mM NaCl treatment. It was reported by Aghaei et al. [13] that growth was inhibited by 100 mM NaCl treatment of 3-day-old soybean seedlings. Based on the present study and the report of Aghaei et al. [13], it is improved that soybean is a relatively NaCl sensitive crop and that these deleterious effects of NaCl are mainly due to up-or down-regulation of genes and their corresponding proteins following NaCl treatment. In this study, Na content was increased under salt stress, but K content did not decrease (see Additional file 1). This result suggests that Na toxicity leads to damaging effects of NaCl in soybean. In order to investigate the effects of NaCl on 7-day-old soybean seedlings, the leaf, hypocotyl and root of 40 mM NaCl-treated seedlings were used for proteome analysis.
Proteins changed in soybean leaves by salt stress show that photosynthesis related proteins are mainly downregulated and suggest that NaCl affects photosynthesis and leads to energy reduction inside the plant and consequent reduction in plant growth. Calreticulin is an important calcium-binding protein with chaperone functions and plays a pivotal role in regulating calcium homeostasis and protein folding in the endoplasmic reticulum of plants [25,26]. Calreticulin was down-regulated in rice under osmotic stress [27]. This indicates the main role of calcium as a main secondary messenger for soybean seedlings under salt stress. The principal role of RuBisCO activase is to release inhibitory sugar phosphates such as ribulose-1,5-biphosphate from the active sites of RuBisCO so that CO 2 can activate the enzyme by carbamylation [28]. RuBisCO activase also has another role as a chaperone during stress [29]. RuBisCO activase was down-regulated by NaCl and this led to the main inhibitory effect of NaCl on soybean photosynthesis [5]. The 20-kDa chaperonin functions only as a cochaperone, along with cpn60, and in certain cases is essential for the discharge of biologically active proteins from cpn60 [30]. Chaperones act to repair the potential damage caused by misfolding of proteins. Most newly synthesized proteins can fold in the absence of chaperones, but a minority strictly requires them [31]. The 20-kDa chaperonin was up-regulated in this study, suggesting that protection of proteins by the chaperone in soybean is very important to prevent misfolding of proteins under salt stress. The large 50S ribosomal subunit catalyses the peptidyl-transfer reaction of messenger RNA-directed protein biosynthesis [32]. Down-regulation of the 50S ribosomal protein indicates the inhibitory effect of NaCl on soybean protein biosynthesis and presumably leads to the consequent reduction in plant growth (Figure 1).
Proteins changed in soybean hypocotyls by salt stress show that metabolic related proteins in the hypocotyl of soybean seedlings are mainly affected by down-regulation under salt stress. Kinesin proteins are a large family of plus-or minus-end-directed microtubule motors important in intracellular transport, mitosis, meiosis, and development [33]. Control and maintenance of the cell cycle and cell integrity are critically important for cell-to-cell communication and signaling, especially under stress [34]. Kinesin motor protein was up-regulated, suggesting that it maintains the cell cycle and cell integrity under salt stress. NADH dehydrogenase 1 beta subcomplex subunit 8 is an accessory subunit of the mitochondrial membrane respiratory chain NADH dehydrogenase, which is not involved in catalysis. This protein functions in the transfer of electrons from NADH to the respiratory chain. The immediate electron acceptor for the enzyme is ubiquinone [35]. In this study, NADH dehydrogenase 1 beta subcomplex subunit 8 was down-regulated, suggesting that this led to a decrease in the ATP pool of the soybean seedling cells and a subsequent decrease in plant growth under salt stress. Trypsin inhibitor is capable of reducing dehydroascorbate when in the reduced form but acquires trypsin-inhibiting activity in the oxidized state [36]. Proteolysis would be accelerated under stress and in response to stress conditions, accumulation of trypsin inhibitors will occur [37]. Trypsin inhibitor was up-regulated and this indicates its main role in prevention of protein proteolysis and in H 2 O 2 detoxification in soybean seedlings under salt stress. The role of alcohol dehydrogenase in tolerance of crops to flooding stress has already been identified [38]. Under stress, plant respiration shifts from aerobic to anaerobic conditions to produce ATP via the fermentative pathway. Alcohol dehydrogenase was up-regulated in the hypocotyl in this study, indicating the main role of this enzyme produce ATP and consume glycolytic products under salt stress. The role of annexin is important for integrity of the cell membrane and cell signaling, and especially under stress conditions, this protein family can help the plant to tolerate the stress and maintain the integrity of the cells [39,40]. Annexin was up-regulated in this study, which indicates the main role of this protein in tolerance of salt stress in the soybean seedling hypocotyl.
Proteins changed in soybean roots by salt stress show that metabolic related proteins in the roots of soybean seedling are mainly affected by down-regulation under salt stress. Dienelactone hydrolase hydrolyzes the conversion of dienelactone to maleylacetate, which both intermediates for aerobic degradation of haloaromatic compounds [41,42]. A dienelactone hydrolase family protein was down-regulated in this study and this could cause a decrease in degradation of these secondary metabolites under salt stress. Caffeoyl-CoA-O-methyltransferase can catalyze the conversion of caffeoyl-CoA to sinapoyl-CoA [43]. These products are the intermediates for lignification of the cell wall and down-regulation of this enzyme may indicate a reduction in cell wall lignification and a consequent decrease in soybean seedling growth under salt stress (Figure 1). Isoflavone reductase catalyzes the NADPH-dependent reduction of 2'-hydroxyformononetin to vestitone, which is the penultimate step in the synthesis of medicarpin in general flavonoid biosynthesis [44]. Down-regulation of this enzyme indicates that flavonoid compounds do not have a main role in tolerance to salt stress in soybean seedlings. Quinone is the product of polyphenol oxidase and is converted to melanin after polymerization [45]. Quinone oxidoreductase was down-regulated and this could be due to a reduction in quinone production or in the amount of phenol compounds in soybean seedling roots under salt stress. Ripening related proteins are involved in plant maturation and ripening [46]. A relationship between the stress hormone ethylene and plant ripening has already been reported [47]. Ethylene induces ripening related proteins and leads the plant to mature and ripen fruit sooner than usual under stress. Up-regulation of ripening related protein could be due to an increase in ethylene concentration of soybean seedlings under salt stress.
Proteome responses of soybean as a salt-sensitive crop to salt stress are different compared to other plants that are mentioned in present study. In soybean, metabolism related proteins were mainly down-regulated with NaCl treatment, but for Salicornia europaea, energy production and ion homeostasis associated proteins played important roles for this plant salt tolerance ability [17] or overlapping of hydrogen peroxide and nitric oxide caused the acclimation of citrus plants to salt stress [16]. For Suaeda aegyptiaca as a halophytes plant, proteins associated with salt tolerance such as proteins involved in glycinebetain synthesis and oxidative stress tolerance had main roles to tolerate salt stress [19]. Caffeoyl-CoA-O-methyltransferase was up-regulated by salt stress at rice [20], but for present study it was downregulated at roots (Table 1C). Lignification changes are important in protecting plants from abiotic stress. This indicates different strategies and subsequently different salt tolerant potential between soybean and other plants.
Based on this study, stem 31 kDa glycoprotein precursor shows developmental specificity and putative fructokinase 2 shows organ specificity in response to NaCl treatment. Glyceraldehyde-3-phosphate dehydrogenase was down-regulated at both the mRNA and protein levels in response to NaCl treatment, suggesting that it plays a role in salt stress and can be used as a target gene in soybean seedlings. Jeong et al. [48] improved salt tolerance of potato by transfer of the glyceraldehyde-3 phosphate dehydrogenase gene. This gene has a main role in tolerance to salt stress and there are several reports related to its relationship to improved salt tolerance in plants [49]. Under salt stress, which inhibits photosynthesis, the substrates for glycolysis decrease and there is a resulting decrease in the rate of the glycolytic reactions. At the step catalyzed by glyceraldehyde-3-phosphate dehydrogenase in glycolysis, one NADH is produced from NAD + . It is indicated that the ATP production will be reduced by down-regulation of glyceraldehydes-3-phosphate dehydrogenase and consequently there will be a decrease in plant growth under salt stress.

Conclusions
The growth of 7-day-old soybean seedlings was affected and their protein expression was changed in the leaf, hypocotyl and root by salt stress. Glyceraldehyde-3phosphate dehydrogenase was down-regulated at both the protein and mRNA levels in leaves and hypocotyls of soybean seedling by salt stress (Figure 7, see Additional files 5, 6 and 7). ATP production will be reduced by down-regulation of glyceraldehydes-3-phosphate dehydrogenase and consequently there will be a decrease in plant growth under salt stress. Based on previous and present study, this gene can be one of the target genes to improve salt tolerance of soybean in the future.

Plant materials and treatment
Seeds of soybean (Glycin max [L.] Merr. cv. Enrei) were sown on the sand and then treated with 0, 20, 40 or 80 mM NaCl. They were grown in a growth chamber under white fluorescent light (600 μmol m -2 s -1 ; 16 h light/8 h dark) at 25°C and 70% relative humidity. Length and fresh weight of hypocotyls and roots of healthy seedlings were measured 7 days after sowing. Physiologically independent experiments were carried out 3 times, and 5 soybean plants were used for each experiment.

Analysis of Na and K contents
Freeze-dried hypocotyls, roots and leaves were used for analysis of Na and K contents. The samples were powdered by grinding with pestl and mortar. A portion of the powered samples were digested with concentrated HNO 3 at 110°C for 2 h. Na and K contents in the digested samples were determined by atomic absorption.

Protein extraction
A portion (200 mg and 500 mg, respectively) of hypocotyls and roots from soybean was homogenized in phosphate saline buffer (pH 7.6) containing 65 mM K 2 HPO 4 , 2.6 mM KH 2 PO 4 , 400 mM NaCl and 3 mM NaN 3 at 4°C using a glass mortar and pestle on ice. The homogenate was centrifuged at 15,000 × g for 10 min, and then cold trichloroacetic acid was added to the supernatant to a final concentration of 10%. The solution was kept on ice for 30 min and then centrifuged for 10 min at 15,000 × g at 4°C. The resultant precipitate was homogenized with 100 μL of lysis buffer containing 8.5 M urea, 2.5 M thiourea, 5% CHAPS, 100 mM dithiothreitol and 0.5% Bio-Lyte (pH 3.0-10.0 and pH 5.0-8.0) in a glass mortar and pestle on ice. A portion (200 mg) of leaf was directly homogenized in lysis buffer containing 7 M urea, 0.2 M thiourea, 0.2 mM TBP, 0.4% CHAPS, 0.2% Bio-Lyte (pH 3.0-10.0) and 5% polyvinylpyrrolidone 40. All homogenates were centrifuged twice at 15,000 × g at room temperature for 10 min. The supernatant was subjected to electrophoresis.

Two-dimensional polyacrylamide gel electrophoresis
Proteins were separated by 2-DE in the first dimension by isoelectric focusing tube gel (11 cm) and in the second dimension by SDS-PAGE [12]. For the first dimension, the isoelectric focusing gel consisted of 8 M urea, 3.5% polyacrylamide, 2% Nonidet P-40, 2% Bio-Lyte (pH 3.0-10.0 and pH 5.0-8.0), ammonium peroxodisulfate and tetramethylethylenediamine. Electrophoresis was carried out at 200 V for 30 min, followed by 400 V for 16 h and 600 V for 1 h. After separation, SDS-PAGE as the second dimension was performed using a 15% polyacrylamide gel (160 × 140 × 1 mm) with a 5% stacking gel. The gels were stained with Coomassie brilliant blue (CBB) and image analysis was performed. 2-DE analysis was done 5 times for independent samples. For each sample 2-DE gels were casted at least 2 times.

Image acquisition and data analysis
Three gels out of 5 gels were used for image acquisition and data analysis. Spot detection, spot measurement, background subtraction and spot matching were performed after CBB staining of the gels using PDQuest software (version 7.1, Bio-Rad, Hercules, CA, USA). Following automatic spot detection, gel images were carefully edited. Before spot matching, one of the gel images was selected as a reference gel. The amount of a protein spot was expressed as the volume of that spot, which was defined as the sum of the intensities of all the pixels that make up the spot. In order to correct the variability due to CBB staining and to reflect the quantitative variation in intensity of protein spots, the spot volumes were normalized as a percentage of the total volume in all spots in the gel. The resulting data from image analysis were transferred to PDQuest software for querying protein spots that showed quantitative and qualitative variations. The pI and Mr of each protein were determined using 2-DE markers (Bio-Rad).

N-terminal amino acid sequence analysis and homology searches
Following separation by 2-DE, proteins were electroblotted onto a polyvinylidine difluoride (PVDF) membrane (Pall, Port Washington, NY, USA) using a semidry transfer blotter (Nippon Eido, Tokyo, Japan), and detected by CBB staining. The stained protein spots were excised from the PVDF membrane and applied to the upper glass block of the reaction chamber in a gasphase protein sequencer (Procise cLC, Applied Biosystems, Foster City, CA, USA). Edman degradation was performed according to the standard program supplied by the manufacturer. The amino acid sequences obtained were compared with those of known proteins in the Swiss-Prot, PIR, GenPept and PDB databases with the Web-accessible search program FASTA http://www. dna.affrc.go.jp.

Analysis using matrix-assisted laser desorption ionization time-of-flight mass spectrometry
Proteins were excised from 2-DE gels stained by CBB and then alkylation and protein digestion with trypsin were done using a DigestPro96 robotic system (Intavis AG, Koeln, Germany). The generated peptides were purified using NuTip C-18 columns (Glygen Corp., Columbia, MD, USA). The purified peptides were added to an α-cyano-4-hydroxycinamic acid matrix and dried onto a plate for analysis using a matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) MS (Voyager-DE RP, Applied Biosystems). Calibration was external and data were collected in the reflector mode. Data were searched on the Internet using an in-house licensed Mascot search engine (version 2.2.18, Matrix Science, Ltd., London, UK) against all entries in the soybean genome database (version 4, 62,199 sequences), which was especially constructed for this research based on preliminary soybean genome sequences from the Department of Energy (DOE) Joint Genome Institute and the Soybean Genome Sequencing Consortium. Genome sequences were downloaded from the DOE database (http://www.phytozome.net, release data 24 January, 2008) and then converted into FASTA format. Carbamidomethylation of cysteines was set as a fixed modification and oxidation of methionines was set as a variable modification. Trypsin was specified as the proteolytic enzyme and one missed cleavage was allowed. Identified proteins with a peptide mass fingerprint were denoted as having an unambiguous identification by the following criteria: (1) the deviation between the experimental and theoretical peptide masses needed to be less than 50 ppm; (2) at least five different predicted peptide masses were needed to match the observed masses for an identification to be considered valid; (3) the matching peptides needed to cover at least 10% of the known protein sequence; and (4) protein scores needed to have > 60 identity for soybean DOE database and > 81 identity for NCBI database (p < 0.05).
RNA extraction and quantitative RT-PCR from the leaf, hypocotyl and root of soybean seedlings Total RNA was extracted from 100 mg each of leaves, hypocotyls and roots of soybean seedlings using an RNeasy plant mini kit with DNAse (Qiagen, Valencia, CA, USA) treatment. A portion (800 ng) of total RNA from each sample was reverse-transcribed to cDNA in a 20 μL reaction volume using iScript cDNA synthesis kit (Bio-Rad) according to the manufacturer's protocol. Real-time quantitative PCR was performed using the cDNA product corresponding to 20 ng of total RNA in a 20 μL reaction volume using the iQ SYBR Green supermix (Bio-Rad) on a MyiQ single-color real-time PCR detection system (Bio-Rad). The PCR conditions were as follows: 95°C for 3 min, then 45 cycles of 95°C for 30 sec, 58°C for 30 sec and 72°C for 30 sec. To normalize gene expression, 18S rRNA (X02623) was used as an internal control. A dilution series of 18S rRNA cDNA clone and the cDNA products were used as templates to confirm amplification effieciency. The relative amount of mRNA was calculated by the ΔCt calculation [50]. The calculation was performed by following calculation formula: 2 (threshold cycle of 18S rRNA-threshold cycle of target gene) ·10 8 . The primers were designed using the Pri-mer3 web interface http://frodo.wi.mit.edu/ [51]. Primer sets for the genes examined had the following sequences: for the stem 31 kDa glycoprotein precursor, the forward sequence was ACCTTCACCTCTCTCAA-CAATC and the reverse sequence was AAAGCCGCG-TAAGAGACAAC; for glyceraldehyde-3-phosphate dehydrogenase, the forward sequence was GAGGGG-GATCGAGTTTTTG and the reverse sequence was CACGAATCCATGATTCAAAGC; for the putative fructokinase 2, the forward sequence was AGAG-GAGGCTGCACTGAAAC and the reverse sequence was GCGTTATTTGGCAAAACCTG and for 18S rRNA, the forward sequence was TGATTAACAGGGA-CAGTCGG and the reverse sequence was ACGGTA TCTGATCGTCTTCG.
Additional file 1: Na and K contents in the leaf, hypocotyl and root of soybeans treated with NaCl. Soybeans were spwn on the sand and treated with 0, 20, 40 or 80 mM NaCl. They were grown for 2 weeks, and Na and K contents of leaf, hypocotyls and root were measured. Ten plants in each treatment were used. The experiments were repeated 3 times and the results show the average ± SE.
Additional file 2: Relative levels of protein expression of leaf of soybeans treated with NaCl. Relative intensities of leaf were obtained using PDQuest software. Values are the mean ± SE. Error bars show SE of 3 spots from 3 independent experiments. The numbers of spots are the same as in Figure 2. White and black columns are the control and NaCl treatments, respectively.
Additional file 3: Relative levels of protein expression of hypocotyl of soybeans treated with NaCl. Relative intensities of hypocotyl were obtained using PDQuest software. Values are the mean ± SE. Error bars show SE of 3 spots from 3 independent experiments. The numbers of spots are the same as in Figure 3. White and black columns are the control and NaCl treatments, respectively.
Additional file 4: Relative levels of protein expression of root of soybeans treated with NaCl. Relative intensities of root were obtained using PDQuest software. Values are the mean ± SE. Error bars show SE of 3 spots from 3 independent experiments. The numbers of spots are the same as in Figure 4. White and black columns are the control and NaCl treatments, respectively.
Additional file 5: Salt stress responsive proteins in leaves of soybean seedlings. a) Spot No, Spot number; b) The sequence shown is the N-terminal amino acid sequence determined by protein sequencing; c) Accession No, Accession number; d) Exp. Mr/pI shows experimental molecular weight and isoelectric point; e) Theor. Mr/pI shows theoretical molecular weight and pH isoelectric; f) PM, Number of matched peptides; g) SC, Sequence coverage by peptide mass fingerprinting using MALDI-TOF MS; h) U & D, up-regulated and down-regulated spots based on significant (p < 0.05) differences between control and NaCl treatments; i) CV ± SE, Spot volume of control ± standard error; j) TV ± SE, Spot volume of treatment ± standard error; k) T/C ratio, Treatment spot volume/control spot volume ratio; l) Category shows functional classification; m) ND, Not determined; M, metabolism; P, photosynthesis; D, defence; CT, cell transfer; PD, protein destination and storage; PS, protein synthesis; T, transcription.
Additional file 6: Salt stress responsive proteins in hypocotyls of soybean seedlings. a) Spot No, Spot number; b) The sequence shown is the N-terminal amino acid sequence determined by protein sequencing; c) Accession No, Accession number; d) Exp. Mr/pI shows experimental molecular weight and isoelectric point; e) Theor. Mr/pI shows theoretical molecular weight and pH isoelectric; f) PM, Number of matched peptides; g) SC, Sequence coverage by peptide mass fingerprinting using MALDI-TOF MS; h) U & D, up-regulated and down-regulated spots based on significant (p < 0.05) differences between control and NaCl treatments; i) CV ± SE, Spot volume of control ± standard error; j) TV ± SE, Spot volume of treatment ± standard error; k) T/C ratio, Treatment spot volume/control spot volume ratio; l) Category shows functional classification; m) ND, Not determined; M, metabolism; P, photosynthesis; D, defence; CT, cell transfer; PD, protein destination and storage; PS, protein synthesis; T, transcription. Additional file 7: Salt stress responsive proteins in roots of soybean seedlings. a) Spot No, Spot number; b) The sequence shown is the Nterminal amino acid sequence determined by protein sequencing; c) Accession No, Accession number; d) Exp. Mr/pI shows experimental molecular weight and isoelectric point; e) Theor. Mr/pI shows theoretical molecular weight and pH isoelectric; f) PM, Number of matched peptides; g) SC, Sequence coverage by peptide mass fingerprinting using MALDI-TOF MS; h) U & D, up-regulated and down-regulated spots based on significant (p < 0.05) differences between control and NaCl treatments; i) CV ± SE, Spot volume of control ± standard error; j) TV ± SE, Spot volume of treatment ± standard error; k) T/C ratio, Treatment spot volume/control spot volume ratio; l) Category shows functional classification; m) ND, Not determined; M, metabolism; P, photosynthesis; D, defence; CT, cell transfer; PD, protein destination and storage; PS, protein synthesis; T, transcription.