Abstract
Vegetable soybean is an important economic and nutritious crop. In this study, 48 differentially expressed proteins were identified from filling seeds of soybean (Glycine max) cv. Mindou 6 by using two-dimensional electrophoresis (2-DE) combined with liquid chromatography electrospray ionization tandem mass spectrometry (LC-ESI-MS/MS). Among them, 25% were related to protein destination and storage, 42% to energy and metabolism, 15% to disease/defense, 6% to transporters, 4% to secondary metabolism, 4% to transcription, 2% to protein synthesis, and 2% to cell growth/division. Along with the maturity of seeds, the number of unchanged abundance proteins decreased, while that of both upregulated and downregulated proteins increased. Both downregulated expression of caffeic acid O-methyltransferase (COMT) and upregulated expression of sucrose-binding protein (SBP) 2 precursor may contribute to increase in digestibility, nutritional value, and eating quality of vegetative seeds at suitable picking period. The pattern of unchanged proteins during the whole seed-filling stage may be also beneficial to the quality of vegetable soybean.
Vegetable soybean, called “Edamame” in Japan or “Mao Dou” in China, is one of the oldest cultivated vegetables and an excellent source of proteins, carbohydrates, dietary fibers, vitamins, minerals, and phytoestrogens (Kim et al., 2014). Therefore, the immature seeds of such a soybean are directly consumed as a kind of vegetable or snack. Now, it is an important economic and nutritious crop in Southeast Asian countries or areas including China, Japan, and Taiwan, and it increasingly grows in the Western Hemisphere (Arikit et al., 2011; Juwattanasomran et al., 2011). Because immature seeds of vegetable soybean, compared with seeds of normal soybean varieties, have added-value traits, such as high sweet taste, distinct fragrance, and large size, they have high commercial price and demand. So soybean breeders always give high priority to these characters in breeding.
In recent years, a great progress in improving quality traits of vegetable soybean has been made. Normally, protein and oil contents of fresh vegetable soybean seeds range from 33% to 39% and 13% to 16%, respectively (Rao et al., 2002). The different vegetable soybeans show considerable variations in the contents of sugar (15.131–33.979 mg·g−1 dry matter), total free amino acids (4.581–10.180 mg·g−1 dry matter), and organic acids (3.754–6.752 mg·g−1 dry matter) (Song et al., 2013). Zhang et al. (2015) reported that the cumulative contribution to eating quality of sucrose, fructose, glucose, raffinose, stachyose, proteins, free amino acids, and oil was 67.8%. Li et al. (2012) reported that sucrose was the highest content (≈70%) of soluble sugar at the edible stage of vegetable soybean. It is found that there are a significantly positive correlation between eating quality score and sucrose content and a significantly negative correlation between eating quality score and protein content (Zhang et al., 2015). Such an organoleptic quality of vegetable soybean is proved to be dependent on various factors, including different varieties, the harvest stage, the duration between the harvest and processing date, and storage conditions (Czaikoski et al., 2013; Li et al., 2012; Mozzoni et al., 2009; Song et al., 2013).
Despite the completion of the soybean genome sequencing (Schmutz et al., 2010), functional genomics study of soybean is still challenging, and therefore proteomics approaches could alternatively be a powerful tool (Arruda et al., 2013; Komatsu and Ahsan, 2009). Hajduch et al. (2005) reported that both 2-DE and MALDI-TOF/TOF overall decrease in metabolism-related proteins and increase in proteins associated with destination and storage during seed filling. Comparative proteomics of soybean and rapeseed revealed increased expression of proteins related to glycolytic and fatty acid metabolism in rapeseed, suggesting that a possible mechanism for higher oil rapeseed is involved in the concerted commitment of hexoses to glycolysis and eventually to de novo fatty acid synthesis pathways (Agrawal et al., 2008).
Until now, little is known about the protein expression changes in vegetable soybean, especially in the aspects of proteomics. To understand protein profiles related to organoleptic properties of the vegetable soybean, in this study, we carried out a detailed proteomic analysis of soybean seeds at the seed-filling stage using 2-DE followed by LC-MS/MS.
Materials and Methods
Plant materials and growth conditions.
The vegetable soybean (Glycine max) Mindou 6 was newly bred by the Institute of Crop Sciences, Fujian Academy of Agricultural Sciences, and approved by Crop variety approval committee of Fujian Province in 2013. In fall 2014, soybean Mindou 6 seeds were grown in the experimental field of Institute of Crop Sciences, Fujian Academy of Agricultural Sciences. According to morphological features of soybean, its reproductive growth can be separated into eight stages: R1 and R2 depicting flowering; R3 and R4 depicting pod development; R5 and R6 depicting seed development; and R7 and R8 depicting plant mature (Mcwilliams et al., 1999). Soybean seeds stayed green at R5 and R6 and began to turn yellow at R7 stage (Tekrony et al., 1979), therefore R5 and R6 are termed as the stage of seed filling. To investigate the dynamic protein profiles of seeds at the filling stage, the flowering days of 100 plants were tagged. Only the similar sized pods from the same part of the tagged plants were harvested and a total of 20 pods from 20 different plants were mixed as a sample. Pods including seeds were sampled at 6-d intervals from initiation of seed formation to suitable to eat as vegetable soybean. Pods appeared light yellow at 67 days after flowering (DAF67). At the whole seed-filling period (DAF49–67), the color of seeds had no significant difference. At DAF49, 55, 61, and 67, developing pods including seeds were harvested and stored in liquid nitrogen until used.
Protein sample preparation.
A total of 10 similar sized pods from 20 pods at any stages were screened and stripped the pods, then only similar sized seeds excluding pods were used for protein sample preparation. Developing soybean seed proteins were extracted using a trichloroacetic acid/acetone procedure following Zhang et al. (2011). The protein concentration of the final supernatant was measured by Bradford approach, using bovine serum albumin as the standard (Bradford, 1976).
2-DE, gel staining, and image analysis.
A 1-DE was performed in a Bio-Rad PROTEAN IEF Cell. The protein extract was diluted to a final concentration of 800 μg·mL−1 with an IEF rehydration solution [7 m urea, 2 m thiourea, 4% CHAPS (w/v), 50 mm DTT, and 0.5% (v/v) IPG buffer (pH 5–8)]. After 15 min of centrifugation at 10,000 gn, 300 μL supernatant was added to each commercially available precast 17-cm-long IPG strip with a 5–8 pH gradient and then actively rehydrated at 50 V for 13 h at 20 °C. Focusing was performed in a Bio-Rad PROTEAN IEF Cell under the following conditions: 100 V for 1 h, 500 V for 1 h, 1000 V for 2 h, 8000 V for 4 h, and then 8000 V for 7.5 h and ≈60,000 Vh was achieved. Before SDS-PAGE, the strips were equilibrated for 15 min in 10 mL reducing equilibration buffer [6 m urea, 0.375 m Tris-HCl (pH 8.8), 2% (w/v) SDS, 20% glycerol (v/v), and 2% (w/v) DTT]. The strips were subsequently placed for another 15 min in alkylating equilibration buffer containing 2.5% (w/v) iodoacetamide instead of 2% DTT.
The second dimension ran on a 12% polyacrylamide SDS gel using an Ettan™ DALT SIX System (GE Healthcare). Gel electrophoresis was carried out at 20 °C with a 1.0 w/gel for 40 min and then with 10 w/gel until the dye front reached about 1 cm from the bottom of the gel.
The gels were stained with silver nitrate as described by Yan et al. (2000). The silver-stained 2-DE gels were digitized with a UMAX PowerLook 2100XL scanner (UMAX Systems GmbH, Willichm, Germany), and analyzed using PDQuest software (version 7.1; Bio-Rad, Hercules, CA). Only those changes that were reproducible and showed at least a 2-fold increase or decrease, as well as being statistically significant (using Student’s t test, at P < 0.05), were considered to be differentially expressed protein spots.
Protein in-gel digestion and LC-ESI-MS/MS.
The in-gel digestion of protein spots was carried out following Gharahdaghi et al. (1999). The AB Sciex Triple TOF 5600 System is a hybrid quadruple TOF-MS equipped with a DuoSpray ionization source coupled to a Shimadzu 30 series high-performance liquid chromatography system. Full scan mass spectrometric analysis and product ion multiple tandem mass spectra (MS/MS) analysis using information-dependent acquisition experiments were performed using the 5600 TF on the reduced and reduced/alkylated injected samples. The LC separation was achieved using a Thermo C18 4.6 × 150 mm column at a linear 1.3% b [90% acetonitrile/0.1% formic acid (aq)] per minute gradient with a flow rate of 0.3 mL·min−1 over 60 min (Dutertre et al., 2013). A cycle of one full scan of the mass range (MS) (300–2000 m/z) followed by MS/MS was applied using a rolling collision energy relative to the m/z and charge state of the precursor ion up to a maximum of 80 eV. The ProteinPilot™ 5.0 software (AB SCIEX, Framingham, MA) was used for sequence identification by searching the LC-ESI-MS/MS mass lists obtained at a mass tolerance of 0.05 Da for precursor ions using the reduced and reduced/alkylated samples. Cysteine carbamidomethylation and methionine oxidation were set as fixed and variable modifications, respectively, accepting two missed cleavages per peptide. The search was conducted against the subset of soybean protein sequences downloaded from the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/) and concatenated with the reverse one. The threshold “Conf” value for accepting identified spectra was set to 95. Identified peptide sequences were inspected manually to confirm assignment.
Functional assignment and cellular localization of proteins.
The identified proteins were annotated as their biological function according to Bevan et al. (1998). The subcellular localization prediction of 48 differently abundant proteins was based on Plant-PLoc (Chou and Shen, 2008) (http://www.csbio.sjtu.edu.cn/bioinf/plant/).
Results
Bioinformatic analysis of 48 identified protein spots.
To understand proteomic changes of soybean-developing seeds, 2-DE coupled with LC-ESI-MS/MS was used to analyze soybean seeds at DAF49, 55, 61, and 67. About 800 spots were detected on each gel for different stages (Fig. 1; Supplemental Fig. 1) and analyzed by the PDQuest 7.1 software. A total of 48 spots showed significant changes in abundant (Fig. 2) and successfully identified by LC-ESI-MS/MS. With Plant-PLoc software (http://www.csbio.sjtu.edu.cn/bioinf/plant/) (Chou & Shen, 2008), these proteins were categorized into eight subcellular compartments: chloroplast (10), cytoplasm (8), extracellular (4), mitochondrion (6), nucleus (2), plasma membrane (10), plastid (3), and vacuole (5) (Supplemental Table 1). No proteins found to localize in cell wall, endoplasmic reticulum, Golgi apparatus, and peroxisome. The prediction of protein functions was categorized according to the method of Beavan et al. (1998). The 48 identified proteins were related to nine metabolic pathways: protein destination and storage (25%), energy (23%), metabolism (19%), disease/defense (15%), transporter (6%), secondary metabolism (4%), transcription (4%), protein synthesis (2%), and cell growth/division (2%).

Representative spot maps of filling seeds of soybean (Glycine max) cv. Mindou 6 at different stages. Protein samples were loaded onto 17-cm immobilized strips (pH 5–8) for IEF, followed by 12% SDS-PAGE. Four different stages of seed development are shown by 49 days after flowering (DAF49), 55, 61, and 67.
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16

Representative spot maps of filling seeds of soybean (Glycine max) cv. Mindou 6 at different stages. Protein samples were loaded onto 17-cm immobilized strips (pH 5–8) for IEF, followed by 12% SDS-PAGE. Four different stages of seed development are shown by 49 days after flowering (DAF49), 55, 61, and 67.
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16
Representative spot maps of filling seeds of soybean (Glycine max) cv. Mindou 6 at different stages. Protein samples were loaded onto 17-cm immobilized strips (pH 5–8) for IEF, followed by 12% SDS-PAGE. Four different stages of seed development are shown by 49 days after flowering (DAF49), 55, 61, and 67.
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16

Identification of 48 differently expressed protein spots from filling seeds by two-dimensional electrophoresis and liquid chromatography electrospray ionization tandem mass spectrometry. Arrow and numbers are indicating the differential spots at days after flowering DAF55, 61, and 67 in comparison with DAF49.
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16

Identification of 48 differently expressed protein spots from filling seeds by two-dimensional electrophoresis and liquid chromatography electrospray ionization tandem mass spectrometry. Arrow and numbers are indicating the differential spots at days after flowering DAF55, 61, and 67 in comparison with DAF49.
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16
Identification of 48 differently expressed protein spots from filling seeds by two-dimensional electrophoresis and liquid chromatography electrospray ionization tandem mass spectrometry. Arrow and numbers are indicating the differential spots at days after flowering DAF55, 61, and 67 in comparison with DAF49.
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16
In comparison of protein profiles among four different stages by Venn diagram, a reference map between differentially expressed proteins is clearly shown in (Fig. 3). Compared with DAF49, DAF55, 61, and 67 had abundant protein spots identified and 18 protein spots commonly expressed both DAF61 and 67, suggesting that seeds expressed more proteins along with the developmental progress. Only three proteins are shared by three stages (Fig. 3), inferring that some proteins accumulated at different stages.

Differences of protein abundance in filling seeds at days after flowering (DAF55), 61, and 67 in comparison with DAF49.
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16

Differences of protein abundance in filling seeds at days after flowering (DAF55), 61, and 67 in comparison with DAF49.
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16
Differences of protein abundance in filling seeds at days after flowering (DAF55), 61, and 67 in comparison with DAF49.
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16
As for the expression level of 48 proteins, 15% of proteins were upregulated and 13% of proteins were downregulated, while 73% of proteins were unchanged in DAF55 compared with DAF49. In the DAF61 and 67 samples, 21% and 31% of proteins were overexpressed, whereas 27% and 44% of the proteins were repressed, respectively (Table 1).
Expressed analysis of 48 identified protein spots at days after flowering (DAF55), 61, and 67 in comparison with DAF49.


Principal component analysis for 48 identified protein spots.
Principal component analysis for 48 protein spots was carried out between DAF49 and DAF55, 61, or 67 using Student’s t test (P < 0.05, at least 2-folds). The results showed that a 2-D (plotting PC1 and PC2) representation showed the effective separation of samples, using identified spots, the abundance, and the number of developing soybean seeds proteins (Fig. 4). The results indicated that the 48 identified protein spots were centralized into two principal components (PCs), of which PC1 and PC2 explained 60.04% and 14.30% of the variance, respectively.

Presentation of the samples based on main principal components found after principal component analysis. Two-dimensional plot of main principal components (PC1 and PC2).
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16

Presentation of the samples based on main principal components found after principal component analysis. Two-dimensional plot of main principal components (PC1 and PC2).
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16
Presentation of the samples based on main principal components found after principal component analysis. Two-dimensional plot of main principal components (PC1 and PC2).
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16
Discussion
Proteins involved in storage.
Soybean storage proteins account for ≈70% to 80% of total seed proteins. These storage proteins, such as glycinin and β-conglycinin, are largely responsible for the nutritional and physicochemical properties of soybean seeds. The three major allergen proteins of soybean are GlymBd 60K, GlymBd 30K, and GlymBd 28K (Mona et al., 2013). Various seed maturation proteins are also synthesized during the later stages of seed development (Mona et al., 2014). Saldivar et al. (2011) reported the level of proteins decreased 2–6% during the first 3–5 weeks after flowering and gradually increased thereafter until maturity, which attributes to the rapid synthesis of oil and starch in early seed development. The best suitable pattern of protein accumulation in breeding is that the protein content remains unchanged until fully maturity (at R6–R7 stages). In our study, among the five protein spots, glycinin (spot 10, 20) were downregulated, whereas β-conglycinin (spot 13, 18, and 41) were upregulated. The results show that soybeans make a large adjustment in the proteome during seed filling and selectively increase accumulation of some proteins.
Proteins involved in transporters.
Soluble sugar is an important eating-quality trait in vegetable soybeans. It is important to increase its content for improving the taste and flavor. Sucrose is the main carbohydrate. It not only functions as a transport metabolite, but also contributes to the osmotic driving force for phloem translocation and serves as a signal to activate or repress expression of specific genes in different tissues (Yang et al., 2004). The SBP has been implicated in the sucrose uptake system in Nicotiana tabacum (Pedra et al., 2000), pea (Castillo et al., 2000), and Vicia faba (Heim et al., 2001). Two identified protein spots with increasing amount (36 and 44) were mapped to the same NCBInr protein (SBP2, gi|351722438), which was in agreement with the previous studies in Arabidopsis microarray and soybean proteome. Transcriptome investigation of Arabidopsis seeds indicated increasing expression of SBP through seed filling (Ruuska et al., 2002). The systematic proteome study of soybean seed filling identified SBP with increased tendency (Martin et al., 2005). Similarly, as members of nutrient transport, outer envelope protein (OEP16) solute channels (OEP of 16 kDa) are found to be specifically expressed in embryo and seed development (Pudelski et al., 2010). Due to the alternating expression pattern in Pisum sativum seeds, OEP16.1 is involved in accumulation processes during early seed maturation. When a transcriptional and metabolic switch induces the ABA-controlled seed desiccation and dormancy, OEP16.2 comes into a player, bridging late seed development with early germination (Pudelski et al., 2012). In consequence, the loss of OEP16 causes metabolic imbalance, in particular that of amino acids during seed development and early germination. In our study, spot 16 was indentified as OEP16–2. The expression ratio at DAF61 and 67 vs. DAF49 was 3.15- and 2.17-folds, respectively.
Proteins involved in secondary metabolism.
Only two identified protein spots (COMT, spot 17 and isoflavone reductase homolog 1, spot 25) were involved in secondary metabolism. Soybean isoflavone synthase (IFS) genes encode the key enzymes involved in the phenylpropanoid pathway. Phenylpropanoid products, such as isoflavonoids and flavonoids, play diverse roles in response to different biotic and abiotic stresses (Yu and Mcgonigle, 2005). COMT gene encodes an enzyme, which catalyzes O-methylation at the C5 position of 5-hydroxyconiferaldehyde and 5-hydroxyconiferyl alcohol to yield sinapaldehyde and sinapyl alcohol, respectively (Jung et al., 2013; Louie et al., 2010). Downregulation of the COMT gene in rapeseed alters the lignin content and composition, improve forage quality, and increases digestion efficiency (Bhinu et al., 2009; Oraby and Ramadan, 2015). The protein level of COMT (spot 17) was downregulated in seed developmental process of Mindou 6. The result suggests that downregulated expression is contributed to increase digestibility, nutritional, and energetic values of vegetable soybean.
Proteins involved in cell growth/division.
Protein spot 3 was identified as cell division cycle protein 48 (CDC48) homolog. CDC48 is an essential and conserved ATP-driven chaperone in eukaryotic cells that functions in diverse cellular processes (Annette et al., 2013). Exogenous GA3 upregulates CDC48 and enhances rice leaf expansion (Wang et al., 2013). CDC48 homolog (spot 3) was upregulated in soybean seed filling process in our study, indicating that CDC48 may be responsible for seed cell expansion in soybean.
Proteins involved in disease/defense.
Proteolysis is fundamental for the normal function of multicellular organisms and plays key roles in a variety of processes. The serine proteases are one of the best characterized groups of proteolytic enzymes in higher organisms. Subtilisin-like proteases are serine proteases (Dodson and Wlodawer, 1998), which fulfill highly specific functions in plant development cascades (Figueiredo et al., 2014). The role of downregulation of subtilisin precursor (spot 14, 21) in the process of soybean seed filling remains to be further studied. Catalases, such as superoxide dismutase, CAT, and APX are key antioxidant enzymes in the defense pathway (Molassiotis et al., 2005; Muneer et al., 2012). Upregulation of superoxide dismutase [Mn] (mitochondrial-like, spot 27), methylecgonone reductase-like (spot 43), 1-Cys peroxiredoxin-like (spot 45), and lambda class glutathione S-transferase (GST, spot 7) during soybean seed filling suggests that the defense system is also enhanced in soybean seed development as a previous report (Pellinen et al., 2002).
Conclusions
This study applied a proteomic approach to understand the protein change of developing seeds of vegetable soybean, and provided informative data for vegetable-soybean breeding. Comparative proteomics of soybean filling seeds identified 48 proteins by 2-DE combined with LC-ESI-MS/MS, which is involved in various biological processes. Along with the maturity of seeds, the number of proteins with unchanged in abundance decrease, but the number of both upregulated and downregulated proteins increase; such tendencies are all beneficial to the quality of vegetable soybean. Especially, both downregulation of COMT and upregulation of SBP 2 precursor may contribute to increase digestibility, nutritional value, and eating quality.
Literature Cited
Agrawal, G.K., Hajduch, M., Graham, K. & Thelen, J.J. 2008 In-depth investigation of the soybean seed-filling proteome and comparison with a parallel study of rapeseed Plant Physiol. 148 504 518
Annette, N., Khalid, A. & Manfred, H. 2013 CDC48 function during TMV infection: Regulation of virus movement and replication by degradation? Plant Signal. Behav. 8 e22865
Arikit, S., Yoshihashi, T., Wanchana, S., Tanya, P., Juwattanasomran, R., Srinives, P. & Vanavichit, A. 2011 A PCR-based marker for a locus conferring aroma in vegetable soybean (Glycine max L.) Theor. Appl. Genet. 122 311 316
Arruda, S.C., Barbosa, H.S., Azevedo, R.A. & Arruda, M.A. 2013 Comparative studies focusing on transgenic through cp4EPSPS gene and non-transgenic soybean plants: An analysis of protein species and enzymes J. Proteomics 93 107 116
Bevan, M., Bancroft, I., Bent, E., Love, K., Goodman, H., Dean, C., Bergkamp, R., Dirkse, W., Van Staveren, M., Stiekema, W., Drost, L., Ridley, P., Hudson, S.-A., Patel, K., Murphy, G., Piffanelli, P., Wedler, H., Wedler, E., Wambutt, R., Weitzenegger, T., Pohl, T.M., Terryn, N., Gielen, J., Villarroel, R., De Clerck, R., Van Montagu, M., Lecharny, A., Auborg, S., Gy, I., Kreis, M., Lao, N., Kavanagh, T., Hempel, S., Kotter, P., Entian, K.-D., Rieger, M., Schaeffer, M., Funk, B., Mueller-Auer, S., Silvey, M., James, R., Montfort, A., Pons, A., Puigdomenech, P., Douka, A., Voukelatou, E., Milioni, D., Hatzopoulos, P., Piravandi, E., Obermaier, B., Hilbert, H., Düsterhöft, A., Moores, T., Jones, J.D.G., Eneva, T., Palme, K., Benes, V., Rechman, S., Ansorge, W., Cooke, R., Berger, C., Delseny, M., Voet, M., Volckaert, G., Mewes, H.-W., Klosterman, S., Schueller, C. & Chalwatzis, N. 1998 Analysis of 1.9 Mb of contiguous sequence from chromosome 4 of Arabidopsis thaliana Nature 391 485 488
Bhinu, V.S., Li, R., Huang, J., Kaminskyj, S., Sharpe, A. & Hannoufa, A. 2009 Perturbation of lignin biosynthesis pathway in Brassica napus (canola) plants using RNAi Can. J. Plant Sci. 89 441 453
Bradford, M.M. 1976 A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding Anal. Biochem. 72 248 254
Castillo, J., Rodrigo, M.I., Márquez, J.A., Zúñiga, Á. & Franco, L. 2000 A pea nuclear protein that is induced by dehydration belongs to the vicilin superfamily Eur. J. Biochem. 267 2156 2165
Chou, K.C. & Shen, H.B. 2008 Cell-PLoc: A package of Web servers for predicting subcellular localization of proteins in various organisms Nat. Protoc. 3 153 162
Czaikoski, K., Leite, R.S., Mandarino, J.M.G., Carrão-Panizzi, M.C., Silva, J.B.D. & Ida, E.I. 2013 Canning of vegetable-type soybean in acidified brine: Effect of the addition of sucrose and pasteurisation time on color and other characteristics Ind. Crops Prod. 45 472 476
Dodson, G. & Wlodawer, A. 1998 Catalytic triads and their relatives Trends Biochem. Sci. 23 347 352
Dutertre, S., Jin, A.H., Kaas, Q., Jones, A., Alewood, P.F. & Lewis, R.J. 2013 Deep venomics reveals the mechanism for expanded peptide diversity in cone snail venom Mol. Cell. Proteomics 12 312 329
Figueiredo, A., Monteiro, F. & Sebastiana, M. 2014 Subtilisin-like proteases in plant-pathogen recognition and immune priming: A perspective Front. Plant Sci. 5 739
Gharahdaghi, F., Weinberg, C.R., Meagher, D.A., Imai, B.S. & Mische, S.M. 1999 Mass spectrometric identification of proteins from silver-stained polyacrylamide gel: A method for the removal of silver ions to enhance sensitivity Electrophoresis 20 601 605
Hajduch, M., Ashwin, G., Stein, J.W. & Thelen, J.J. 2005 A systematic proteomic study of seed filling in soybean. Establishment of high-resolution two-dimensional reference maps, expression profiles, and an interactive proteome database Plant Physiol. 137 1397 1419
Heim, U., Wang, Q., Kurz, T., Borisjuk, L., Golombek, S., Neubohn, B., Adler, K., Gahrtz, M., Sauer, N., Weber, H. & Wob, U. 2001 Expression patterns and subcellular localization of a 52 kDa sucrose-binding protein homologue of Vicia faba (VfSBPL) suggest different functions during development Plant Mol. Biol. 47 461 474
Jung, J.H., Vermerris, W., Gallo, M., Fedenko, J.R., Erickson, J.E. & Altpeter, F. 2013 RNA interference suppression of lignin biosynthesis increases fermentable sugar yields for biofuel production from field-grown sugarcane Plant Biotechnol. J. 11 709 716
Juwattanasomran, R., Somta, P., Chankaew, S., Shimizu, T., Wongpornchai, S., Kaga, A. & Srinives, P. 2011 A SNP in GmBADH2 gene associates with fragrance in vegetable soybean variety “Kaori” and SNAP marker development for the fragrance Theor. Appl. Genet. 122 533 541
Kim, S.Y., Kim, B.M., Kim, J.B., Shanmugavelan, P., Kim, H.W., Kim, S.Y., Kim, S.N., Cho, Y.S., Choi, H.S. & Park, K.M. 2014 Effect of steaming, blanching, and high temperature/ high pressure processing on the amino acid contents of commonly consumed Korean vegetables and pulses Prev. Nutr. Food Sci. 19 220 226
Komatsu, S. & Ahsan, N. 2009 Soybean proteomics and its application to functional analysis J. Proteomics 72 325 336
Li, Y., Du, M., Zhang, Q.Y., Wang, G.H., Hashemi, M. & Liu, X.B. 2012 Greater differences exist in seedprotein, oil, total soluble sugar and sucrosecontentof vegetable soybean genotypes [Glycine max (L.) Merrill] in Northeast China Austral. J. Crop Sci. 6 1681 1686
Louie, G.V., Bowman, M.E., Tu, Y., Mouradov, A., Spangenberg, G. & Noel, J.P. 2010 Structure-function analyses of a caffeic acid O-methyltransferase from perennial ryegrass reveal the molecular basis for substrate preference Plant Cell 22 4114 4127
Martin, H., Ganapathy, A., Stein, J.W. & Thelen, J.J. 2005 A systematic proteomic study of seed filling in soybean. Establishment of high-resolution two-dimensional reference maps, expression profiles, and an interactive proteome database Plant Physiol. 137 1397 1419
Mcwilliams, D., Berglund, D. & Endres, G. 1999 Soybean growth and management: Quick guide. North Dakota State University, Fargo, ND
Molassiotis, A.N., Diamantidis, G.C., Therios, I.N., Tsirakoglou, V. & Dimassi, K.N. 2005 Oxidative stress, antioxidant activity and Fe(III)-chelate reductase activity of five Prunus rootstocks explants in response to Fe deficiency Plant Growth Regulat. 46 69 78
Mona, T., Alkharouf, N.W., Khan, F.H. & Savithiry, N. 2013 SoyProDB: A database for the identification of soybean seed proteins Bioinformation 9 165 167
Mona, T., Alkharouf, N.W. & Natarajan, S.S. 2014 SoyProLow: A protein database enriched in low abundant soybean proteins Bioinformation 10 599 601
Mozzoni, L.A., Chen, P., Morawicki, R.O., Hettiarachchy, N.S., Brye, K.R. & Mauromoustakos, A. 2009 Quality attributes of vegetable soybean as a function of boiling time and condition Intl. J. Food Sci. Technol. 44 2089 2099
Muneer, S., Kim, T.H. & Qureshi, M.I. 2012 Fe modulates Cd-induced oxidative stress and the expression of stress responsive proteins in the nodules of Vigna radiata Plant Growth Regulat. 68 421 433
Oraby, H.F. & Ramadan, M.F. 2015 Impact of suppressing the caffeic acid O-methyltransferase (COMT) gene on lignin, fiber, and seed oil composition in Brassica napus transgenic plants Eur. Food Res. Technol. 240 931 938
Pedra, J.H.F., Delú-Filho, N., Pirovani, C.P., Contim, L.A.S., Dewey, R.E., Otoni, W.C. & Fontes, E.P.B. 2000 Antisense and sense expression of a sucrose binding protein homologue gene from soybean in transgenic tobacco affects plant growth and carbohydrate partitioning in leaves Plant Sci. 152 87 98
Pellinen, R.I., Minna-Sisko, K., Tauriainen, A.A., Tapio, P.E. & Jaakko, K.R. 2002 Hydrogen peroxide activates cell death and defense gene expression in birch Plant Physiol. 130 549 560
Pudelski, B., Kraus, S., Soll, J. & Philippar, K. 2010 The plant PRAT proteins - preprotein and amino acid transport in mitochondria and chloroplasts Plant Biol. 12 Suppl. 1 42 55
Pudelski, B., Schock, A., Hoth, S., Radchuk, R., Weber, H., Hofmann, J., Sonnewald, U., Soll, J. & Philippar, K. 2012 The plastid outer envelope protein OEP16 affects metabolic fluxes during ABA-controlled seed development and germination J. Expt. Bot. 63 1919 1936
Rao, M.S.S., Bhagsari, A.S. & Mohamed, A.I. 2002 Fresh green seed yield and seed nutritional traits of vegetable soybean genotypes Crop Sci. 42 1950 1958
Ruuska, S.A., Thomas, G., Christoph, B. & Ohlrogge, J.B. 2002 Contrapuntal networks of gene expression during Arabidopsis seed filling Plant Cell 14 1191 1206
Saldivar, X., Wang, Y.J., Chen, P. & Hou, A. 2011 Changes in chemical composition during soybean seed development Food Chem. 124 1369 1375
Schmutz, J., Cannon, S.B., Schlueter, J., Ma, J., Mitros, T., Nelson, W., Hyten, D.L., Song, Q., Thelen, J.J. & Cheng, J. 2010 Genome sequence of the palaeopolyploid soybean Nature 463 178 183
Song, J., Liu, C., Li, D. & Gu, Z. 2013 Evaluation of sugar, free amino acid, and organic acid compositions of different varieties of vegetable soybean (Glycine max [L.] Merr) Ind. Crops Prod. 50 743 749
Tekrony, D.M., Egli, D.B., Balles, J., Pfeiffer, T. & Fellows, R.J. 1979 Physiological maturity in soybean 1 Agron. J. 71 771 775
Wang, X., Han, F., Yang, M., Yang, P. & Shen, S. 2013 Exploring the response of rice (Oryza sativa) leaf to gibberellins: A proteomic strategy Rice (N. Y.) 6 1 12
Yan, J.X., Wait, R., Berkelman, T., Harry, R.A., Westbrook, J.A., Wheeler, C.H. & Dunn, M.J. 2000 A modified silver staining protocol for visualization of proteins compatible with matrix-assisted laser desorption /ionization and electrospray ionization-mass spectrometry Electrophoresis 21 3666 3672
Yang, Z., Lei, Z., Diao, F., Huang, M. & Wu, N. 2004 Sucrose regulates elongation of carrot somatic embryo radicles as a signal molecule Plant Mol. Biol. 54 441 459
Yu, O. & Mcgonigle, B. 2005 Metabolic engineering of isoflavone biosynthesis Adv. Agron. 86 147 190
Zhang, Q., Li, Y., Liu, C., Tian, B., Xu, B. & Mao, J. 2015 Key components of eating quality and their dynamic accumulation in vegetable soybean varieties [Glycine max (L.) Merr.] Acta Agron. Sin. 41 1692 1700 (in Chinese with English abstract)
Zhang, Y.M., Zhao, J.M., Xiang, Y., Bian, X.C., Zuo, Q.M., Qi, S., Gai, J.Y. & Xing, H. 2011 Proteomics study of changes in soybean lines resistant and sensitive to Phytophthora sojae Proteome Sci. 9 52

Three repeats of spot maps of filling seeds of Mindou 6 at 49 days after flowering (DAF49), 55, 61, and 67.
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16

Three repeats of spot maps of filling seeds of Mindou 6 at 49 days after flowering (DAF49), 55, 61, and 67.
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16
Three repeats of spot maps of filling seeds of Mindou 6 at 49 days after flowering (DAF49), 55, 61, and 67.
Citation: HortScience 51, 11; 10.21273/HORTSCI10911-16
Identification of 48 differently expressed protein spots from filling seeds of cv. Mindou 6 (Glycine max) at DAF55, 61 and 67 in comparison with DAF49.

