Rice seeds (O. sativa L. cv. Wuyunjing No. 7, obtained from company of Nanjing Shenzhou Seed) were surface-sterilized with 5% (v/v) sodium hypochlorite (NaClO) for 15 min and thoroughly washed in distilled water. Each treatment was performed in triplicate. For one replicate, 100 seeds were randomly placed on moist filter paper in 200 mm Petri dish. The seeds were germinated in the dark at 25°C with renewal of distilled water every day. After 4 days, these germinating rice seeds were transferred to the mesh over 2.5 L vessel containing different concentrations of Cu sulfate pentahydrate (CuSO₄·5H₂O) solution (0–200μM) for 0–48 h, 1 mM ascorbic acid (Asc) for 12 h or 10 mM H₂O₂ solution for 6 h. A certain number of radicles from each replicate were obtained for the below experiments.

H₂O₂ formation in situ in rice radicles was visually detected based on the infiltration of 3,3′-diaminobenzidine (DAB) as described by Romero-Puertas et al. (2004) with minor modifications. Briefly, six radicles from each replicate (each Petri dish) were immersed in a 1 mg/mL solution of DAB (pH 3.8) and incubated at room temperature for 20 min in the absence of light. After staining, images were captured with a Coolpix 4500 digital camera (Nikon, Tokyo, Japan).

The concentration of H₂O₂ in rice radicles from Cu-treated plants was measured by monitoring the A415 of the titanium-peroxide complex following the method described by Jiang and Zhang (2001). Absorbance values were calibrated to a standard curve established with 0.1–1.0 μM H₂O₂.

Rice seeds germinated for 4 days were treated with 10 mM H₂O₂ for 6 h or with 20 and 100 μM Cu for 12 h, referred as low and high concentration Cu treatment, respectively. Seeds germinated in deionized water without Cu and H₂O₂ were used as controls. Radicles were harvested and ground with a mortar and pestle in liquid nitrogen to obtain a fine powder, and then suspended in four volumes of protein binding buffer (20 mM sodium phosphate, pH 5.8, 500 mM NaCl, 0.1% (w/v) Triton X-100) containing 1 mM phenylmethyl sulfonyl fluoride (PMSF), incubated for 30 min at 4°C, and centrifuged for 30 min at 15,000 g at 4°C. The proteins in the supernatant were used for protein analysis, and the protein concentration was determined according to the Bradford method using bovine serum albumin (0, 0.2, 0.4, 0.6, 0.8, 1.0 mg/mL) as the standard (Bradford, 1976).

Experimental design for proteomic analysis of Cu-binding proteins in rice radicles was shown in Supplementary Figure S1. The used method for separating and isolating Cu-binding proteins was based on the Cu-IMAC method of Song et al. (2014). Protein samples were pre-chromatographed on a column with IDA-Sepharose for removing metal ions from proteins samples before flowing over a Cu-IDA-Sepharose column (referred to as Cu-IMAC) for separating Cu-binding proteins. IDA-Sepharose and Cu-IMAC were connected in tandem with a tube (inner diameter of 0.5 mm). For the pre-chromatography column, IDA-Sepharose with a 2 mL bed volume for each column was poured into a 10-mL glass column with an inner diameter of 10 mm and washed with 10 mL of water at a rate of 0.5 mL/min. For the Cu-IMAC column, IDA-Sepharose with a 2 mL bed volume for each column was poured and washed with 10 mL of water at a rate of 0.5 mL/min, after which the bed volume of 0.2 M CuSO₄ was applied to the column, followed by washing with 15 mL distilled water to remove excess Cu ions at a rate of 0.5 mL/min. Columns were equilibrated with 10 bed volumes of binding buffer at a linear flow rate of 0.5 mL/min after they were connected in tandem.

A total of 20 mg of protein sample solution was loaded onto the pre-chromatography column at a linear flow rate of 0.25 mL/min, after which the column was washed with binding buffer until the ultraviolet absorbance of the effluent from the Cu-IMAC column reached baseline. Binding buffer containing 10 mM imidazole was added to the per-chromatography column at a linear flow rate of 0.5 mL/min until the UV-absorbance of the effluent from the Cu-IMAC column reached baseline, suggesting that non-specific proteins were removed. The Cu-binding proteins were eluted with elution buffer (10 mM sodium acetate, 500 mM NaCl, pH 5.5) containing 40 mM imidazole at a linear flow rate of 0.5 mL/min until the UV-absorbance of the effluent from the Cu-IMAC column reached baseline.

The Cu-binding proteins eluted from Cu-IMAC were precipitated with four volumes of cooled 10% (w/v) acetone (containing 0.07% (w/v) dithiothreitol, DTT) overnight at −20°C for 1 h, followed by centrifugation for 15 min at 10,000 g at 4°C. Afterward, the pellets were dissolved in lysis buffer (7 M urea, 2 M thiourea, 4% (w/v) CHAPS, 0.2% (w/v) Bio-lytes pH 3–10, 65 mM DTT). Protein concentrations were assayed using a Bio-Rad RC DC Protein Assay Kit 1.

For each replicate, 100 μg of total protein extract was loaded onto IPG dry strips (17 cm, pH 4–7 linear gradient; Bio-Rad, Hercules, CA, USA) during the rehydration step (13 h), followed by focusing for a total of 60,000 V·h using a Protean IEF Cell (Bio-Rad). Following isoelectric focusing (IEF), the gel strips were equilibrated for 15 min in 5 mL equilibration buffer containing 0.375 M Tris–HCl (pH 8.8), 6 M urea, 20% (v/v) glycerol, 2% (w/v) sodium dodecyl sulfate (SDS), and 2% (w/v) DTT. The strips were then equilibrated in the same buffer as described above, but including 2.5% w/v iodacetamide instead of DTT. SDS-polyacrylamide gel electrophoresis (PAGE) in the second dimension was performed using 12% SDS-polyacrylamide gels sealed with 0.5% agarose. Electrophoresis was performed at 50 V for the first 30 min, followed by 150 V for 8 h using a Protean Plus Dodeca Cell apparatus (Bio-Rad). Protein spots were visualized using MS-compatible silver staining (Yan et al., 2000). To prevent the gels being overexposed, the developing course was divided into two steps: firstly, the developing solution was drained off after becoming yellow; secondly, the developing course was terminated when the small protein spots begin to become clear.

The gels were scanned using the image scanner UMAX Powerlook III (UMAX Technologies, Dallas, TX, USA) at 300 dpi resolution; image and data analyses of the gels were performed using PDQuest software (version 8.0; Bio-Rad) and a multivariate statistical package (DeCyder EDA, Unscrambler, Samespots), which can automatically deal with missing values during analysis (Valledor and Jorrín, 2011). The abundance of spot mean a summation of the pixel intensities localized within the defined spot area, which obtained by PDQuest (Bio-Rad) image analysis software. Spot quantity was normalized in the “total quantity of valid spot” mode for possible staining differences between gels. Duplicate 2-DE gels were run for each treatment from three independent tissue extractions, only spots with significant and reproducible changes were considered to represent differentially accumulated proteins. The results for the control and Cu- or H₂O₂-treated samples were analyzed for differences using Student's t-test with a significance level of 5%. Protein spots were selected for MS analysis when a difference of 1.5-fold or greater was observed in the level of accumulation between the treatment and control.

Protein spots were excised and destained (Gharahdaghi et al., 1999). The samples were incubated in 50 mM ammonium bicarbonate for 5 min, dehydrated with acetonitrile (ACN), and dried. The peptides were extracted with 60% ACN and 0.1% trifluoroacetic acid after the proteins were digested with trypsin, and were then extracted and desalted with ZipTip C18 columns (Millipore, Bedford, MA, USA). The peptide solution was saturated with α-cyano-4-hydroxycinnamic acid and then air-dried on an MS sample plate.

Peptide mass spectra were obtained using a 4700 Proteomics Analyzer MALDI-TOF/TOF™ mass spectrometer (Applied Biosystems, Framingham, MA, USA) in positive ion reflector mode. The subsequent MS/MS analysis was performed in a data-dependent manner, and the 10 most abundant ions fulfilling certain preset criteria were subjected to high-energy collisional dissociation (CID) analysis. The collision energy was set to 1 keV, and nitrogen was used as the collision gas.

All protein spectra were submitted for database searching using the online MASCOT program (http://www.matrixscience.com) against NCBInr databases (http://www.ncbi.nlm.nih.gov/protein). The taxonomic category selected was Oryza sativa. The searching parameters were as follows: 0.15 Da mass tolerance for peptides and 0.25 Da mass tolerance for TOF–TOF fragments, one allowed trypsin miscleavage, Cys carbamidomethylation as a fixed modification, and Met oxidation as a variable modification. Only significant hits, as defined by the MASCOT probability analysis (P < 0.05), were accepted.

Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg/ or http://www.kegg.jp/) was used to predict molecular function, biological processes, and significant pathways involved in response to stress.

Total RNA was extracted using the RNA simple Total RNA Kit (LifeFeng, Shanghai, China) according to the manufacturer's instructions and then converted to cDNA after DNase I treatment using a PrimeScript™ RT Master Mix (TaKaRa). Real-time quantitative RT-PCR was performed on a MyiQ Real-Time PCR Detection System (Bio-Rad Hercules, CA, USA) using SYBR Premix Ex Taq (TaKaRa). The primers for protein mRNA were listed in Supplementary Table S1. The PCR protocol included an initial 7-min incubation at 95°C for complete denaturation followed by 40 cycles at 94°C for 30 s, 60°C for 30 s, and 72°C for 30 s. The specificity of the PCR amplification was examined based on a heat dissociation curve (65–95°C) following the final cycle. Normalized relative expression was calculated using the 2^−ΔΔCt (cycle threshold) method.

Data were analyzed using SPSS ver. 16.0 (Statistical Package for Social Science for Windows, SPSS, Inc., Chicago, IL, USA). All values reported in this paper are means ± SE (n = 3) of three separate experiments. Means denoted by the same letter did not significantly differ at P < 0.05 according to Duncan's multiple range test.

The accumulation of H₂O₂ in the radicles of germinating rice seeds was examined using histochemical DAB staining. Our results showed that exposure to excess Cu for 12 h caused an evident accumulation of H₂O₂ in the radicles (Figure 1A). The production of Cu-induced H₂O₂ could be decreased by infiltration with the H₂O₂ scavenger, Asc. In the presence of 100 μM Cu, H₂O₂ accumulation in radicles gradually increased during the first 12 h of exposure and then decreased slightly but remained higher than the control (Figure 1B). The concentrations of H₂O₂ assayed by spectrophotometry were consistent with the results of histochemical detection by DAB staining (Figure 1C).

UV detections to Cu-binding proteins of rice radicles via Cu-IMAC were shown in Supplementary Figure S2. The Cu-binding protein yields from the control, H₂O₂-treated, 20 μM Cu, and 100 μM Cu-treated rice radicles, estimated with the percents of peak area (of total peak area), were not significantly different. Protein maps produced from 2-DE gels showed a high reproducibility among the three independent extractions (Figure 2, Supplementary Figure S3). When analyzed using PDQuest, 780 ± 15, 772 ± 25, 793 ± 13, and 695 ± 24 proteic spots were identified in the range of pH 4–7 and relative molecular masses of 10–120 kDa with the control, H₂O₂-treated, 20 μM Cu and 100 μM Cu-treated rice radicles, respectively. The significantly differential spot patterns between 100 μM Cu treatment and the other treatments might be explained by differing degrees of protein loss resulting from the Cu-binding proteins eluted from the Cu-IMAC column. Among all of the spots, 656 spots were present with all four treatment, other protein spots were treatment-specific. 26 protein spots, exhibited more than 1.5-fold differences in the abundances under at least one treatment (20, 100 μM Cu, or 10 mM H₂O₂) compared to the control. In order to observe more clearly, these protein spots were artificially divided into four regions (A, B, C and D) in gels, and the four regions were enlarged and shown in Figure 3.

Compared with the control, the H₂O₂ treatment increased the abundances of 11 spots and decreased that of nine spots among the 26 protein spots. Among the 11 H₂O₂-increased protein spots, six protein spots were simultaneously increased under both Cu treatments (20 and 100 μM Cu), and five protein spots were increased under 20 μM Cu. Among nine protein spots decreased under H₂O₂, six protein spots were similarly decreased under both Cu treatment, one spot decreased under 20 μM Cu treatment. In addition, the abundances of six protein spots only increased under Cu treatment and unaffected under H₂O₂ treatment (Supplementary Figure S4).

The Cu-binding proteins in 26 spots were analyzed using MALDI-TOF/TOF MS, and all spectra of proteins were submitted to a NCBInr protein database search using the online MASCOT program. Supplementary Table S2 shows the identity of Cu-binding proteins in 26 spot after a database search. The molecular masses (Mr) and isoelectric point (pI) values of each identified protein are listed in Table 1. These identified Cu-binding proteins were found to be involved in different cellular responses and metabolic processes, including antioxidative defense (6 proteins), stress response and detoxification (4 proteins), protein synthesis (5 proteins), protein modification (1 protein in 2 spots), protein metabolism (2 protein in 3 spots), carbohydrate metabolism (3 proteins), nucleotide metabolism (1 protein), and secondary metabolite metabolism (2 proteins).

In this study, among 24 proteins in 26 spots identified, 18 proteins contained one or more of nine metal-binding motifs reported by Smith et al. (2004), and 20 protein species contained one or more of the top six motifs (H–(X)₅–H, H–(X)₇–H, H–(X)₁₂–H, H–(X)₆–M, M–(X)₇–H, and H–(X)₃–C) reported by Kung et al. (2006) in Arabidopsis roots (Table 2). Fifteen proteins contained motifs reported by Smith et al. (2004) and the top six motifs by Kung et al. (2006). However, one protein (spots 18 and 19; putative legumin) contained neither the motifs reported by Smith et al. (2004) nor the top six motifs reported by Kung et al. (2006).

In order to assess the correlation between mRNA expression and protein accumulation, Real-time quantitative RT-PCR was applied to four mRNAs of identified Cu-binding proteins, copper/zinc superoxide dismutase (CuZn-SOD, spot 16), L-ascorbate peroxidase (APX, spot 8), peroxiredoxin (Prx, spot 15), and glutathione S-transferase 2 (GST2, spot 6), involved in antioxidative defense (Figure 4). By compared with the corresponding spots of Figure 3, the expression analyses of four genes were consistent with the proteins accumulation except for CuZn-SOD mRNA change under high Cu treatment, indicate that the accumulation of these proteins have been largely regulated at the transcriptional level.

Numerous studies have shown that excess Cu can induce the formation of ROS (including H₂O₂) and cause oxidative stress. In this study, the formation of H₂O₂ was observed with increasing Cu concentrations and with Cu treatment time in the Cu-treated radicles of rice (Figure 1). Accumulation of H₂O₂ has also been observed in other Cu-exposed plant species using histochemical staining (Tewari et al., 2006; Sgherri et al., 2007; Zhang et al., 2008, 2010). Because H₂O₂ is relatively stable and diffusible through membranes, it is known to modulate gene expression and participate in various physiological processes (Neill et al., 2002; Ahmad et al., 2008). So far, no proteins simultaneously response to H₂O₂ and Cu stress were reported by searching web of science.

In the present study, 24 Cu-IMAC-binding proteins in 26 spots were identified that were differentially accumulated at least by one treatment (20, 100 μM Cu, or 10 mM H₂O₂). The same protein (e.g., protein disulfide isomerase and putative legumin) in varied spots is possible since the spot change its position in the gel due to changes in pI or Mr as a consequence of post-translational modifications. Among these identified proteins, elongation factor EF-2 (EF-2), GST, Prx, APX, quinone-oxidoreductase QR2 (QR2), protein disulfide isomerase (PDI), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), triosephosphate isomerase (TPI), NADPH-dependent mannose 6-phosphate reductase (M6PR), and cytidine/deoxycytidine deaminase-like (CDC) have been identified as Cu-IMAC-binding proteins in Arabidopsis (Kung et al., 2006; Tan et al., 2010), soybean (Wang et al., 2014), microalgae (Smith et al., 2014), and rice (Song et al., 2014). Moreover, the analogs of eukaryotic translation initiation factor 5A (eIF5A), transcription factor (TF), and SOD, such as translation initiation factor Tu, G, 3A, and 4A, iron-dependent transcriptional regulator, and Fe-SOD were identified as Cu-IMAC-binding proteins in hepatoma cells (Smith et al., 2004), in Arabidopsis (Kung et al., 2006) and S. pneumoniae (Sun et al., 2011). SODs are metalloenzymes found in three different molecular forms containing Cu and Zn (CuZn-SOD), Mn (Mn-SOD), or Fe (Fe-SOD) as prosthetic metals. However, to the best of our knowledge, the other eight proteins, including pathogenesis-related protein (PR) known as PR-10a, and Bet v I family protein (PR-b), cytochrome P450-like protein (CYP-L), salt stress-induced protein (SSI), legumin, cathepsin B-like cysteine protease (CBCP), arginine decarboxylase (ADC), and chalcone-flavonone isomerase (CHI), have not been reported as Cu-IMAC-binding proteins in plants. The present data showed that the abundance of 11 spots increased under exogenous H₂O₂ treatment out of 17 spots increased under 20 μM Cu, and that all 20 μM Cu-decreased spots decreased under H₂O₂ (Table 1, Supplementary Figure S3), which may be ascribed to the same amount H₂O₂ produced by low concentration Cu as H₂O₂ treatment. However, H₂O₂ or Cu at high levels can cause oxidative stress and cell damage, which could be the reason that the abundance of five protein spots are decreased under 100 μM Cu stress but increase under H₂O₂ and 20 μM Cu stress.

Four identified Cu-binding proteins such as CuZn-SOD (spot 16), QR2 (spot 13), Prx (spot 15), and APX (spot 8), displayed similar behavior under 20 μM Cu and H₂O₂ treatments, which may play important roles in plant antioxidant defense responses. SODs are key players in the antioxidant defense system through the dismutation of O2•− to H₂O₂. SODs as metal chelators may also regulate the intracellular Cu level. In plants, quinones are redox-active compounds that oxidize the thiol groups of proteins and GSH. QR2 catalyzes two electron reductions of quinones to hydroquinones (Malakshah et al., 2007; Vannini et al., 2012). Prx, which consists of many different thiol-disulfide exchange proteins, such as thioredoxins and glutaredoxins, is an H₂O₂-scavenging enzyme that reduces H₂O₂ to H₂O, and Prx possesses a highly reactive Cys that is oxidized to form a disulfide bond coupled with the reduction in H₂O₂ (Dietz et al., 2006). Thus, increase of these Cu-binding proteins may alleviate Cu-induced damage by decreasing free Cu²⁺ activity in the cytoplasm and/or scavenging ROS.

Notably, the abundance changes of specific Cu-binding proteins responded differentially to excess Cu: the abundances of Prxs and QR2 increased, while that of APX decreased. The abundance of CuZn-SOD increased under 20 μM Cu but decreased under higher levels of Cu (100 μM). These different changes of antioxidative enzymes following exposure to excess Cu may be due to their varied functions. Excess Cu increased SOD expression (Sunkar et al., 2006; Cohu et al., 2009; Zhang et al., 2010) and activities (Tewari et al., 2006; Zhang et al., 2008, 2010). However, Cu ions can be dangerous to cellular compartments as free ions. Thus, high Cu treatment (100 μM Cu) decreased the abundance of CuZn-SOD in this study. APX and CAT are two major scavengers of H₂O₂. APX is present throughout the cell and has a higher substrate affinity in the presence of Asc as a reductant. Cu and H₂O₂ have been reported to increase APX expression and activity (Lee et al., 2007). In this study, the abundance of APX decreased under excess Cu and exogenous H₂O₂ treatments. Decrease of APX could lead to the accumulation of H₂O₂ and enhance oxidative stress. Similar decrease of APX was observed in rice leaves (Wan and Liu, 2008) and tobacco cells (Vannini et al., 2012) in response to H₂O₂ stress at high doses (50 mM) or over extended times (6 h). H₂O₂ was suggested to directly inhibit APX activity by causing protein oxidation at concentrations over a threshold value (de Pinto et al., 2006). The abundance changes of GST (spots 6 and 14) responded differentially to excess Cu and exogenous H₂O₂ treatment. A major function of GSTs is to detoxify a variety of hydrophobic and electrophilic compounds by catalyzing their conjugation with GSH (Jwa et al., 2006). Consistent with our results, an increased GSTs was detected in Cu-treated (Song et al., 2013, 2014) and H₂O₂-treated rice (Wan and Liu, 2008). In contrast, a decrease in GST levels in rice exposed to Cu (Ahsan et al., 2007), H₂O₂ (Vannini et al., 2012), and selenium (Se) (Wang et al., 2012) has been observed.

The gene products of four identified antioxidant proteins (CuZn-SOD, APX, Prx and GST2) showed similar changes obtained from proteomics experiments except for CuZn-SOD change under high Cu treatments (Figure 4). Previous studies showed that Cu availability is the major factor that determines whether Fe-SOD or CuZn-SOD are expressed (Cohu et al., 2009), the CuZn-SOD accumulation is mediated by a microRNA, miR398, which targets CuZn-SOD mRNA for degradation under some condition (Sunkar et al., 2006), and CuZn-SOD proteins accumulated only when Cu ions were available for final assembly and stability. Thus, it is possible that the abundance decrease of CuZn-SOD protein spot under high Cu treatment was not consistent with the results of gene expression analyses.

Three proteins including PR-b, CYP-L, and SSI out of four Cu-binding proteins involved in the stress response and detoxification displayed similar behavior under low concentration Cu and H₂O₂ treatment. In plants, CYP proteins are involved in the synthesis of fatty acids, lignin, hormones, and flavonoids, as well as xenobiotic metabolism in higher plants (Schuler and Werck-Reichhart, 2003). In this study, the abundance of one CYP-L (spot 17) increased slightly under 20 μM Cu and H₂O₂ but significantly decreased under 100 μM Cu treatment. Li et al. (2009) observed an increase of the CYP-like protein in soybean at 2 h post inoculation. In contrast, decrease of CYP proteins was observed in Cu-treated rice germinating embryos (Zhang et al., 2009) and a Cd-treated Phytolacca americana leaf (Zhao et al., 2011). Wan and Liu (2008) observed that one putative salt-induced protein increased under H₂O₂ in rice leaves. In this study, Cu and H₂O₂ treatments significantly increased the abundance of SSI (spot 23), but its function in Cu-stressed plants remains unknown. PR proteins play a role in a wide range of cell functions, including cell wall rigidification, signal transduction, and antimicrobial activity (Markovic-Housley et al., 2003). Elevated levels of ROS have been reported to induce PR proteins in rice (Jwa et al., 2006). The increase of PR proteins has also been observed in Cu-treated Phaseolus vulgaris (Cuypers et al., 2005), Elsholtzia splendens (Li et al., 2009), and rice (Zhang et al., 2009). In this study, treatment with H₂O₂ and Cu decreased or did not affect the abundance of PR-10a (spot 1) and PR-b (spot 2), although 100 μM Cu increased that of PR-10a. The opposite change patterns of PR proteins by excess Cu suggest that they have different roles.

All Cu-binding proteins involved in protein synthesis displayed similar behavior under low concentration Cu and H₂O₂ treatments (Table 1). The abundances of susceptibility homeodomain transcription factor (SHTF, spot 3), EF-2, and (spot 21) increased under excess Cu and H₂O₂, excluding the high Cu treatment that did not affect that of SHTF and decreased that of legumin. In contrast, the abundances of eIF5A (spot 11) and eIF5A-2 (spots 9 and 10) decreased under both Cu and H₂O₂. eIF5A was also thought to play a role in translation elongation (Saini et al., 2009) and other aspects of RNA metabolism such as RNA export (Liu et al., 2008). The expression of eIF5A in plants usually increases in response to abiotic stress (Li et al., 2009; Xu et al., 2011; Meng et al., 2014; Parkash et al., 2014). In agreement with our result, a significant decrease of eIF5A was observed in rice after a long-term salt stress (Parker et al., 2006), which may be associated with premature senescence. EFs (EF1A, EF1B, and EF-2) are fundamental regulatory proteins of the translational elongation step in higher plants, as well as other eukaryotic organisms. EF-2 catalyzes GTP-dependent translocation of peptidyl-tRNA from the A site to the P site of the ribosome during peptide chain elongation (Browne and Proud, 2002). In this study, the abundances of EF-2 increased by 30.0-, 11.4-, and 6.1-fold in the presence of H₂O₂, low and high concentration Cu, respectively. Similar increase of EF-2 was observed in Schizosaccharomyces pombe in response to H₂O₂ stress (Weeks et al., 2006). In contrast, decrease of the EF-2 protein was observed in Cu-treated E. splendens roots (Li et al., 2009), Cd-treated P. americana (Zhao et al., 2011), and B-deficient Brassica napus (Wang et al., 2010).

Legumin is a major storage protein in plant seeds, including α and basic polypeptides of 40 and 20 kDa, respectively, bound by a disulfide bridge (Sabir et al., 1973). This protein contained neither the motifs reported by Smith et al. (2004) nor the top six motifs reported by Kung et al. (2006), but contained 8 of 117 potential metal-binding motifs (C-(X)n-C, C-(X)n-H, C-(X)n-M, H-(X)n-C, H-(X)n-H, H-(X)n-M, M-(X)n-C, M-(X)n-H, and M-(X)n-M, where n = 0–12) reported by Kung et al. (2006). Cu caused a reduction in the germination rate of bean, which increased the level of storage proteins compared to the control (Karmous et al., 2011). In this study, the abundances of legumin (spots 18 and 19) increased under excess Cu and H₂O₂. It is unknown whether increase of legumin protein abundance may alleviate Cu-induced damage by decreasing free Cu²⁺ activity in the cytoplasm or be a consequence of Cu toxicity and oxidative stress.

Two enzymes (GAPDH and M6PR) involved in carbohydrate metabolism and two enzymes (ADC and CHI) involved in secondary metabolism displayed similar behavior under Cu and H₂O₂ treatments. Treatments with both Cu and H₂O₂ decreased the abundance of GAPDH (spot 12) and increased M6PR (spot 20). Arabidopsis Cytosolic GAPDH may be a potential target of H₂O₂-dependent oxidation in plant protein extractions (Hancock et al., 2005). TPI and GAPDH are important enzymes in the glycolytic pathway. Here, the decrease of GAPDH and increase of TPI may favor the accumulation of glyceraldehyde 3-phosphate under stress conditions. M6PR, a key enzyme in mannitol biosynthesis, catalyzes the conversion of mannose 6-phosphate into mannitol 1-phosphate. Overexpression of M6PR genes from celery and Arabidopsis result in increased tolerance to salt stress (Sickler et al., 2007; Chan et al., 2011). Treatments with Cu and H₂O₂ increased abundance of ADC (spot 4) and decreased that of CHI (spot 5). Cu-induced increases in ADC activity were also observed in previous reports (Groppa et al., 2008; Xu et al., 2011). The roles of ADC and CHI as Cu-binding proteins in the tolerance to Cu and oxidative stresses are still unknown.

Whereas five identified Cu-binding proteins including GST (spot 6), PDI (spots 25 and 26), CBCP (spot 24), TPI (spot 7), and CDC (spot 22) were increased only under low and high Cu but unaffected under H₂O₂ stress, which hint Cu ions can regulate the Cu-binding proteins accumulation by no H₂O₂ pathway. PDI has been identified as a Cu-binding protein in previous reports (Smith et al., 2004; Song et al., 2014). PDI is a thioredoxin superfamily oxidoreductase from the endoplasmic reticulum, and catalyzes a wide range of thiol-disulfide exchange reactions, including oxidation, reduction, and isomerization, and also displays chaperone, calcium-binding, and Cu-binding activity (Hatahet and Ruddock, 2009; Laurindo et al., 2012). Overexpression of PDI gene from Methanothermobacter thermoautotrophicum enhances mercury tolerance in transgenic rice (Chen et al., 2012). Proteomic analyses showed that PDI accumulation increased in rice roots in the presence of Cu (Song et al., 2014), in rice leaves in the presence of H₂O₂ (Wan and Liu, 2008), and in soybean leaves during salt stress (Ma et al., 2012). Here, the abundance of PDIs was markedly higher in Cu-treated rice, which can enhance Cu tolerance in germinating rice seed by binding Cu and thiol-disulfide exchange reactions. In contrast, PDIs accumulation were down-regulated by Cu stress in roots of the tolerant plant E. splendens (Li et al., 2009), by H₂O₂ in rice root apoplasts (Zhou et al., 2011), and by flooding stress in soybean roots (Khatoon et al., 2012).

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.