A proteomic analysis of seeds from Bt-transgenic Brassica napus and hybrids with wild B. juncea

Transgene insertions might have unintended side effects on the transgenic host, both crop and hybrids with wild relatives that harbor transgenes. We employed proteomic approaches to assess protein abundance changes in seeds from Bt-transgenic oilseed rape (Brassica napus) and its hybrids with wild mustard (B. juncea). A total of 24, 15 and 34 protein spots matching to 23, 13 and 31 unique genes were identified that changed at least 1.5 fold (p < 0.05, Student’s t-test) in abundance between transgenic (tBN) and non-transgenic (BN) oilseed rape, between hybrids of B. juncea (BJ) × tBN (BJtBN) and BJ × BN (BJBN) and between BJBN and BJ, respectively. Eight proteins had higher abundance in tBN than in BN. None of these proteins was toxic or nutritionally harmful to human health, which is not surprising since the seeds are not known to produce toxic proteins. Protein spots varying in abundance between BJtBN and BJBN seeds were the same or homologous to those in the respective parents. None of the differentially-accumulated proteins between BJtBN and BJBN were identical to those between tBN and BN. Results indicated that unintended effects resulted from transgene flow fell within the range of natural variability of hybridization and those found in the native host proteomes.

allergenic, which could have negative impact on human health. Proteomics analyses have been applied to test for unintended effects in GM crops, such as tomato 6,10 , rice 11,14 , maize 8,15,16 , wheat 17 , pea 18,19 and tobacco 20 . Most of these studies showed that transgenic lines did have some changes in the production of proteins-those that were not targets for genetic engineering [15][16][17]20 . An important question to ask is how altered protein production compares with the range of natural variability. Assessing the proteomics of seeds is especially appropriate for edible seed crops, and relatively facile given the relative compact nature of the seed proteome.
In certain hosts, transgenes could be transferred into hybrids through outcrossing of GM crops with crop varieties or wild relatives, and in some cases, introgressed transgenic advanced generations could occur [21][22][23][24][25] . Gene flow from GM crops to their wild relatives is one main environmental regulatory issue 22,24,26 , with the concern of increasing risks, such as increased weediness 27 . In examples of transgenes conferring increased insect, herbicide or virus-resistance, there might be, in turn, a competitive advantage of GM plants when cultivated with non-GM plants 26,28 . Previous studies of transgene flow have focused on effects on plant phenotypic and agronomic characteristics, plant fitness and ecological risks [29][30][31][32] .
Oilseed rape (Brassica napus) has been a widely-used crop to study the ecological consequence of gene flow 25 . Wild brown mustard (B. juncea), an allotetraploid wild relative species of B. napus, is a widespread weed in agricultural fields in China and elsewhere. Hybrids between B. napus and wild B. juncea are successfully obtained by open pollination 21,30 . The aim of this study was, for the first time, to perform a proteomic study in a transgenic B. napus and B. napus/B. juncea hybrid system to better understand potential unintended effects of GM event.

Results
We compared proteomes between conventional (BN) and transgenic (tBN) B. napus seed or among wild B. juncea (BJ) and their hybrid BJ × BN (BJBN) and BJ × tBN (BJtBN) seeds by the 2-D electrophoresis (2-DE, Fig. 1). Proteomic comparison between tBN and BN seeds detected the potential unintended effects of GM event. By comparing BJtBN versus BJBN seed proteomes, it is possible to evaluate the unintended influences by transgene flow. At the same time, the comparison of BJBN versus BJ seed proteomes could investigate the natural variability of hybridization (Fig. 1).

Proteomic analysis of GM event effects: tBN vs. BN.
Approximately 800 protein spots were detected in 2-D gels of tBN and BN seeds after Coomassie brilliant blue R-250 (CBB) staining (representative image in Fig. 2A,B). We considered the proteomic differences between tBN and BN seeds significant if they had greater than a 1.5-fold change in abundance at a P-level of p < 0.05 using Student's t-test, which resulted in a total of 31 protein spots meeting these criteria ( Fig. 2A,B, Supplemental Table S1). Proteins were successfully identified in 28 spots by searching the PMFs of each protein spot against the NCBInr database. Among these spots, 24 were comprised of just one protein, which could be assigned to 23 respective unique genes ( Table 1). The other four genes were matched to two proteins (Supplemental Table S2). If one spot consisted of two proteins, it is difficult to deconvolute protein identity. Therefore, the only spots that were analyzed herein are those with a one-to-one match.
The proteins were classified into six functional categories according to Bevan et al. (1998) 33 (Table 1) and mainly involved in four functional categories, storage proteins (33%), cell defense and rescue (25%), energy (17%), and metabolism (17%) (Fig. 2C). Eight spots accumulated in higher abundance in tBN than in BN seeds (Fig. 2D, Table 1). Among these spots, six were storage proteins, and the other two involved in energy production and cell defense and rescue (Fig. 2D, Table 1). The other 16 protein spots were accumulated in lower abundance in tBN than in BN seeds (Fig. 2D, Table 1). Five of them involved in cell defense and rescue, four in metabolism and three in energy production (Fig. 2D, Table 1). Four protein spots, BnaC06g06810D (spot 12), Cruciferin storage protein (spot 8), BnaC03g41580D (spot 15) Proteomic analysis of hybridization effects: BJBN vs. BJ. Thirty-seven protein spots changed more than 1.5-fold in abundance in hybrid BJBN seeds compared to the seeds of its the maternal BJ parent ( Fig. 4A,B, Supplemental Table S4). Proteins were identified in 35 spots, of which 34 spots matched to only one protein and derived from 31 unique genes (Table 3). One spot (spot 56) matched to two proteins (Supplemental Table S2) and was excluded in the following analysis.
The identified protein spots that varied in abundance between BJBN and BJ were annotated to cell defense and rescue (41%), storage protein (26%) and metabolic functions (21%) (Fig. 4C). Twenty-four spots, involved mainly in cell defense and rescue (seven spots), metabolism (seven spots) and storage (seven spots), had higher abundance in BJBN than in BJ seeds (Fig. 4D, Table 3). The abundance of the other ten protein spots was lower in BJBN than in BJ seeds, among which, cell defense and rescue (seven spots) were the most abundant (Fig. 4D, Table 3).
Comparison of proteome changes associated with GM event, transgene flow and hybridization. One protein, BnaC01g09900D exhibited a higher accumulation in tBN than in BN seeds (spot 26) and in BJtBN than in BJBN seeds (spot 32) ( Fig. 5; Tables 1 and 2). However, this protein accumulated in spots 26 and 32 had completely different experimental molecular mass and pI (Tables 1 and 2), indicating that they have been differentially modified after translation. Therefore, none of protein spots accumulated differentially in abundance between BJtBN and BJBN seeds were substantially identical to those between tBN and BN seeds.
Seven protein spots (spot 41,43,73,70,34,61,59) varied similarly in abundance between BJtBN/ BJBN and BJBN/BJ comparisons (  representative CBB R-250 stained gels of BN and tBN seeds, respectively. A total of 400 μ g of proteins was loaded onto the gel strip and separated by 2-D gel electrophoresis and visualized with CBB R-250. The red and blue arrows with numbers indicate the spots whose abundance were more than 1.5 fold higher and lower in tBN than in BN seeds, respectively. The spots were excised, and proteins were extracted and finally identified by MS/MS. Only protein spots successfully identified by MS/MS and matched to one protein were used in the functional analysis (C,D). (C) functional distribution of identified differentially accumulated protein spots between tBN and BN seeds; (D) number and distribution of identified protein spots accumulated in higher or lower abundance in tBN than in BN seeds.  Discussion Proteomic comparison between tBN and BN seeds did not detect new known toxic or nutritionally-harmful compounds. Twenty four proteins were identified to change in abundance in oilseed rape as a result of transgenic events (Fig. 2, Table 1). This number fell within the range of previous proteomic studies on seeds of transgenic crops. A proteomic analysis revealed that 43 protein spots increased or decreased in abundance in GM maize seeds compared with its non-GM isoline 16 . Gong et al. (2012) found that 17 protein spots varied in abundance between insect-resistant transgenic rice harboring Cry1Ac and its non-transgenic control, and 12 spots between herbicide-resistant transgenic rice carrying bar and non-transgenic control 11 . A proteomic study on transgenic tomato reported that transgene insertion had no influence on seed proteomes 10 .
To date no proteomic studies on transgenic plants have raised any new safety concerns 3,5 ; our study did not change this paradigm. In the seeds of oilseed rape, no new protein was found in tBN seeds when compared to BN seeds ( Table 1). Abundance of eight protein spots was higher in tBN than in BN seeds ( Fig. 2D). They are cruciferin gene family proteins (spots 25, 26, 27, 29, 30 and 31), glyceraldehyde-3-phosphate dehydrogenase (spot 24) and catalase (spot 23) ( Table 1). Cruciferin is one type of storage proteins in oilseed rape seeds. Storage proteins are importantly nutritional compounds in seeds and play essential roles in seed germination and seedling growth [34][35][36] . Glyceraldehyde-3-phosphate dehydrogenase is a well-known glycolytic enzyme, and catalase is one of the antioxidant enzymes for removal of reactive oxygen species (ROS) 37,38 . None of these proteins has been reported to be toxic or nutritionally harmful to human health. The result is what was expected, since these two Brassica species are not known to produce toxic proteins in seeds. Nonetheless, our study was about unintended effects, and it is conceivable that there could have been novel protein biosynthesis in the seeds; we found none. Regardless, the spontaneous production of novel allergens and toxins in transgenic is a long-held point of concern by members of the public.
Accumulation of several protein spots in BN seeds were completely inhibited after transgenic modification ( Fig. 2A,B, Table S1). This phenomenon can also be found in other proteomic evaluation of GM events 16,18,19 . It may be due to the random insertion of transgene into plant genome, which inhibits the gene expression and accumulation of these proteins.
In seeds, it appears that storage protein is most susceptible to be altered by transgene insertion. Storage proteins are the common proteins found to vary in abundance between GM and non-GM seeds in many species, such as rice 11 , maize 16 and pea 18,19 . In oilseed rape, we found that about 33% of the differentially accumulated proteins in tBN seeds were storage proteins when compared to BN seeds (Fig. 2C). Transgenic events seemed to influence not only storage protein biosynthesis, but also their post-translational modification. In oilseed rape seeds, storage proteins varying in abundance between transgenic and non-transgenic seeds had different experimental molecular weights (MW), pI, or both, respective to their theoretical MW and pI (Table 1). Similar results have also been found in proteomic studies in seeds of transgenic maize 16 and pea 18,19 . Storage proteins are the most abundant in seeds and account for as much as 60% of total proteins 39 . This finding suggests that storage proteins might be among the most susceptible proteins for alteration in transgenic plants.
Besides storage proteins, some other proteins accumulated differentially in abundance between tBN and BN also exhibited different experimental MW and pI as compared with theoretical MW and pI (Table 1). Similar results were also found in the comparison of BJtBN and BJBN seeds ( Table 2) and in previous proteomic studies on GM plants 9,11,16,18 . These results suggest that a GM event might affect both protein synthesis and post-translational modification.
No protein products derived from the inserted transgenes (Cry1Ac and gfp) were identified in the present proteomic analysis. Similar results have also been found in other proteomic studies on transgenic crops 11,15,16,40 . This may be due to that expression of transgene protein products in tBN seeds is under the sensitivity of CBB staining for detection of protein spots in 2-D gel. In addition, these and other low-produced protein would be occluded by highly synthesized storage proteins, which dominated the 2-D gel (for example, Fig. 2A Proteomic changes in transgenic hybrids were within the bounds of natural variability. As in tBN seeds, proteomics analysis found no transgene protein products in BJtBN seeds. However, a Bt ELISA detected accumulation of Bt-toxin proteins in the BJtBN seeds (Supplemental Figure S1), showing gene flow to the wild relative 41,42 . Transgene flow has potential effects on plant genome, morphology traits and plant growth and reproduction, which may confer altered host ecology 22,26,27 . Table 1. Identified protein spots that are differentially accumulated between transgenic (tBN) and nontransgenic (BN) Brassica napus seeds. Spot ID is the spot number shown in Fig. 2; accumulation pattern indicates the spot whose abundance is higher (H) or lower (L) in tBN seeds than BN seeds; Ratio tBN/ BN, normalized spot volume in tBN seeds divided by the normalized volume in BN seeds; No. of sequence peptides, the peptides matched by the MS/MS spectra (Ion score > 20 Transgene flow is likely to result in unintended effects in transgenic hybrids. Compared with non-transgenic hybrid seeds of BJBN, fifteen protein spots varied in abundance in the transgenic hybrid seeds of BJtBN (Fig. 3A,C, Table 2). BJBN and BJtBN plants were cultivated together in the same environment. Therefore, the observed variation was likely attributed primarily to transgene movement into wild B. juncea. Among these proteins, storage protein was the most abundant (Fig. 3C), which support that storage protein is most sensitive to transgene insertion as discussed above. Proteins involved in cell defense and rescue also accounted for a large difference (Fig. 3C), in which three spots (70, 34 and 38) were identified as late embryogenesis abundant (LEA) proteins (Table 2). LEA proteins are thought to play an essential role in desiccation tolerance and vigor [43][44][45] , by replacing water, sequestering ions, removing ROS and/or stabilizing protein and membrane structure 45,46 . In this category, one protein, BnaC03g55840D, which is homologous to the 17.4 kDa class I heat shock protein (HSP), was newly accumulated in BJtBN seeds (Table 2). HSPs are well known to function as chaperones, which stabilize newly synthesized proteins to ensure correct folding or helping refold damaged proteins from stress.
It appears that there were few unintended effects of the GM event relative to natural hybridization effects, and data fell within the range of natural variability of hybridization. The number of protein spots (15) varied in abundance between BJtBN and BJBN seeds was far less than those (34) found between BJBN and BJ seeds. Among these 15 protein spots in BJtBN/BJBN, seven spots exhibited a same accumulation pattern between BJBN and BJ seeds (Fig. 5), and five spots (40, 32, 58, 72 and 38) were homologous to some of the differentially accumulated protein spots in BJBN/BJ (Tables 2 and 3). For example, spot 40 accumulated in higher abundance in BJtBN than in BJBN seeds is homologous to spots 41 and 42 whose abundance were higher in BJBN than in BJ seeds (Tables 2 and 3). Spots 32, 58 and 72, identified to be cruciferin storage protein accumulated in higher abundance in BJtBN seeds than BJBN seeds ( Table 2). Similar accumulation pattern of cruciferin storage protein (spots 43, 46, 65, 66 and 73) could be observed between BJBN and BJ seeds (Table 3). Spot 38 (LEA protein family protein) in BJtBN/ BJBN is homologous to spot 54 in BJBN/BJ, and they showed a same accumulation pattern in these two pair-wise comparisons (Tables 2 and 3). Herman and Price (2013) reviewed 20 years of research on unintended compositional changes in GM crops and found that the changes were small compared with those produced from traditional breeding and environmental factors 47 . Our results agree with this conclusion and extend it to the next generation of GM crops, namely hybrids with wild relatives.
Some of the differentially accumulated proteins identified by tBN/BN have similar functions as those identified by BJtBN/BJBN. For instances, some cruciferin storage proteins were differentially accumulated between BJtBN and BJBN seeds and between tBN and BN seeds (Tables 1 and 2). BnaC04g48420D (spot 5, Table 1) and LEA protein (spot 18, Table 1) in comparison of tBN vs. BN and At2g42560 (spots 34 and 38, Table 2) in comparison of BJtBN vs. BJBN had lower abundance in transgenic seeds than in non-transgenic seeds; these were members of the LEA protein family. Except for the functional similarity, however, none of the differentially accumulated protein spots were substantially identical between BJtBN/BJBN and tBN/BN seeds (Fig. 5). This indicates that the unintended effects resulting from insertion of exogenous gene in crops were not inherited in the next generation.
According to our experimental design and proteomic methods in this study, there were a few differentially accumulated proteins, but no clear patterns were detected between crops and between crop-wild  napus. (A,B) representative CBB R-250 stained gels of BJBN and BJtBN seeds, respectively. A total of 400 μ g of proteins was loaded onto the gel strip and separated by 2-D gel electrophoresis and visualized with CBB R-250. The red and blue arrows with numbers indicate the spots whose abundance were more than 1.5 fold higher and lower in BJtBN than in BJBN seeds, respectively. The spots were excised, and proteins were extracted and finally identified by MS/ MS. Only protein spots successfully identified by MS/MS and matched to one protein were used in the functional analysis (C,D). (C) functional distribution of identified differentially accumulated protein spots between BJtBN and BJBN seeds; (D) number of identified protein spots accumulated in higher or lower abundance in BJtBN than in BJBN seeds.
Scientific RepoRts | 5:15480 | DOi: 10.1038/srep15480 relative generations. We concluded that potential unintended effects of this transgenic-wild system are negligible and within bounds of natural variation.

Materials and Methods
Plants. Seeds of wild brown mustard (Brassica juncea, 2n = 36, AABB) originating from a local field collection (Nanjing, China) were provided by Prof. S. Qiang, Nanjing Agricultural University. Transgenic oilseed rape (B. napus cv. Westar, 2n = 38, AACC) was produced by transforming with a pSAM12 plasmid containing genetically linked gfp (encoding a green fluorescent protein) and Bt Cry1Ac cassettes that are regulated by independent CaMV 35S promoters 48 .
Wild mustard, non-transgenic and transgenic B. napus were planted in three segregated greenhouses (natural light, average daily temperature varied between 20 and 30 °C) at the Chinese Research Academy of Environmental Sciences (Beijing, China) respectively. After emasculating the stamens of the maternal wild plants, one wild mustard plant (BJ) was hybridized with the pollen from one conventional B. napus (BN) to obtain F 1 hybrid seeds (BJBN), and another wild mustard was crossed with the pollen of one  transgenic B. napus (tBN) to form F 1 hybrid seeds (BJtBN) (Fig. 1). The remainder of BN, tBN and BJ plants were self-pollinated (Fig. 1).
After maturation, seeds of BN, tBN, BJ, BJBN and BJtBN were harvested separately and dried in an air-conditioned room (28 °C, 45% humidity) until the water content reached a constant level (about 0.09 g H 2 O g −1 dry weight). The accumulation of Cry1Ac proteins in BN, tBN, BJ, BJBN and BJtBN seeds were detected using enzyme-linked immunosorbent assays (ELISA) from Agidia (EnviroLogix, USA), and the ELISA was performed following the manufacturer's instructions. The Cry1Ac protein accumulated both in tBN and BJtBN (Supplemental Fig. S1). The negative controls (BN, BJ and BJBN) had no measurable Bt protein.

Preparation of protein samples. Fifty seeds each of BN, tBN, BJ and their hybrids BJBN and BJtBN
were homogenized in a total of 1.5 ml precooled extraction buffer containing 50 mM Tris-HCl (pH 7.5), 30% (w/v) sucrose, 10 mM EGTA, 1 mM PMSF, 1 mM DTT, and 1% (v/v) Triton X-100. The homogenate was centrifuged at 16 000 g for 10 min followed by at 32 000 g for 20 min at 4 °C. The resulting supernatant was mixed with two-fold volume of ice-cold Tris-HCl (0.1 M, pH 7.5) saturated phenol and shaken on ice for 30 min. After centrifugation at 16 000 g for 20 min, the phenol phase was collected and mixed with five times volume of precooled methanol saturated with (NH 4 ) 2 SO 4. Samples were incubated at -20 °C for over 6 hours and then centrifuged at 16 000 g for 5 min. The resulting pellets were rinsed four times with ice-cold acetone (100%) containing 13 mM DTT, and then lyophilized. Lyophilized protein samples were solubilized in lysis buffer composed of 7 M urea, 2 M thiourea, 2% (w/v) CHAPS, 20 mM DTT, and 0.5% (v/v) immobilized pH gradient (IPG) buffer (pH 3-10) and then used for determination of protein concentration by Bradford method 49 using bovine serum albumin as the standard.

2-DE.
Isoelectrofocusing (IEF) was performed using a PROTEAN i12 IEF system (Bio-Rad; Hercules, Piscataway, NJ, USA) and 17 cm Immobiline Dry Strips with a nonlinear pH gradient of 3-10 (Bio-Rad; Hercules, USA). The 400 μ g protein sample dissolved in 300 μ l lysis buffer was loaded onto the gel strip by passive rehydration: the gel strip was incubated at 20 °C for 16 h and then used for IEF. IEF was performed by applying a voltage of 250 V for 1 h, ramping to 500 V over 1 h, 2 000 V for 2 h, 10 000 V for 4 h and held at 10 000 V until a total of 60 kVh was reached. After IEF, the gel strip was reduced for 15 min with 65 mM DTT dissolved in 3 ml of equilibration buffer (6 M urea, 30% (v/v) glycerol, 2% (w/v) SDS, 50 mM Tris-HCl (pH 8.8) and 0.01% (w/v) bromophenol blue) and then alkylated with 2.5% (w/v) iodoacetamide in the same buffer for 15 min. The reduced and alkylated strip was placed onto a vertical SDS-polyacrylamide gel (12% resolving and 5% stacking), and the low-molecular-range marker (Bio-Rad) was loaded at one end of the strip. The strip was sealed with 0.5% (w/v) low-melting agarose in SDS buffer containing bromophenol blue before electrophoresis. Electrophoresis was performed at 15 °C in SDS electrophoresis buffer (pH 8.   . (A,B) representative CBB R-250 stained gels of BJ and BJBN seeds, respectively. A total of 400 μ g of proteins was loaded onto the gel strip and separated by 2-D gel electrophoresis and visualized with CBB R-250. The red and blue arrows with numbers indicate the spots whose abundance were more than 1.5 fold higher and lower in BJBN than in BJ seeds, respectively. The spots were excised, and proteins were extracted and finally identified by MS/MS. Only protein spots successfully identified by MS/MS and matched to one protein were used in the functional analysis (C,D). (C) functional distribution of identified differentially accumulated protein spots between BJBN and BJ seeds; (D) number of identified protein spots accumulated in higher or lower abundance in BJBN than in BJ seeds. Image analysis, in-gel digestion with trypsin and protein identification by MALDI-TOF-TOF mass spectrometry. The 2-D gels were scanned at a 300 dpi resolution with a UMAX Power Look 2100XL scanner (Maxium Tech., Taipei, China). Spot detection and gel comparisons were made with ImageMaster 2D Platimum, version 5.01 (GE Healthcare Bio-Science, Little Chalfont, UK). After automatic detection and matching, manual editing was carried out to correct the mismatched and unmatched spots.
Well-separated gels of the three independent biological replicates were used for proteomic comparisons. Spots were considered reproducible if they were detected in all the biological replicates. The normalized volume (based on the total spot volumes) of each spot was assumed to represent its accumulation abundance. Protein spot was considered to be differentially accumulated when the change was more than 1.5-fold and significant (p < 0.05, Students' t-test). Protein spots, which varied in abundance were excised from the stained gels and digested according to Shevchenko et al. 50 with minor modifications. Gel spots were washed and destained in 100 μ l water followed by 100 μ l 50% acetonitrile (ACN) (Fisher Scientific; Fair Lawn, NJ, USA) twice, respectively and then dehydrated with 50 μ l 100% ACN at room temperature (RT). Dehydrated gel spots were first incubated in 10 mM DTT/100 mM NH 4 HCO 3 for 45 min at 56 °C, followed by incubation with 55 mM iodoacetamide/100 mM NH 4 HCO 3 for 30 min at RT in darkness. After a series of washes (as described above), the gel pieces were dehydrated in 100% ACN. Samples were subsequently rehydrated in digestion buffer containing 10 ng trypsin in 50 mM NH 4 HCO 3 for 45 min at 4 °C. After removal of the excess of digestion buffer, the gel pieces were incubated in 50 mM NH 4 HCO 3 for overnight at 37 °C.
Peptide mixtures were desalted on in-house made microcolumns packed with reverse phase material POROS 20R2 (Applied Biosystems, US) before MS analysis. The desalted samples were eluted from the column with 3 μ l matrix solution (5 μ g/μ l α -cyano-4-hydroxycinnamic acid in 70% ACN, 0.1% TFA) directly onto the MALDI target plate (Opti-TOF Fig. 4; accumulation pattern indicates the spot whose abundance is higher (H) or lower (L) in BJBN seeds than BJ seeds; Ratio BJBN/BJ, normalized spot volume in BJBN seeds divided by the normalized volume in BJ seeds; No. of sequence peptides/matched peptides, exp. protein mass and theo. protein mass: theoretical protein mass are described in Table 1. A: appeared. Proteins in the bracket are the homologous proteins of the identified proteins.
spectra in the range 700-3500 m/z were acquired automatically with external calibration to a standard β -lactoglobulin tryptic digest. Collision-induced dissociation was used for fragmentation of the 10 most intense precursor ions and was performed automatically with default calibration. Protein identification was performed using Biotools 3.2 (Bruker Daltonics, Billerica, MA, US) by searching NCBInr database with Mascot 2.2.04 (Matrix Science). The following search parameters were used: Database -NCBInr (NCBInr 20140927), Taxonomy-Green plants (49886901 sequences; 17905752166 residues), maximum one missed cleavage, cysteine carbamidomethylation as a fixed modification, methionine oxidation and N-terminal acetylation as variable modifications, mass tolerance of 100 ppm in MS mode and 0.6 Da for MS/MS. Protein scores greater than 76 (in NCBInr database) were significant (p < 0.05). The identified protein name, accession number, Mascot score, sequence coverage, number of sequenced peptides (the peptides matched by the MS/MS spectra (ion score > 20) and of matched peptides (the peptides matched by the peptide mass fingerprinting (PMF) and theoretical protein mass and pI found in the databases are shown in Tables 1-3.