Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

An Integrated Proteomic and Metabolomic Study on the Chronic Effects of Mercury in Suaeda salsa under an Environmentally Relevant Salinity

  • Xiaoli Liu,

    Affiliations Key Laboratory of Coastal Zone Environmental Processes, Yantai Institute of Coastal Zone Research (YIC), Chinese Academy of Sciences (CAS); Shandong Provincial Key Laboratory of Coastal Zone Environmental Processes, YICCAS, Yantai Shandong, P. R. China, The Graduate School of Chinese Academy of Sciences, Beijing, P. R. China

  • Huifeng Wu ,

    hfwu@yic.ac.cn

    Affiliation Key Laboratory of Coastal Zone Environmental Processes, Yantai Institute of Coastal Zone Research (YIC), Chinese Academy of Sciences (CAS); Shandong Provincial Key Laboratory of Coastal Zone Environmental Processes, YICCAS, Yantai Shandong, P. R. China

  • Chenglong Ji,

    Affiliations Key Laboratory of Coastal Zone Environmental Processes, Yantai Institute of Coastal Zone Research (YIC), Chinese Academy of Sciences (CAS); Shandong Provincial Key Laboratory of Coastal Zone Environmental Processes, YICCAS, Yantai Shandong, P. R. China, The Graduate School of Chinese Academy of Sciences, Beijing, P. R. China

  • Lei Wei,

    Affiliations Key Laboratory of Coastal Zone Environmental Processes, Yantai Institute of Coastal Zone Research (YIC), Chinese Academy of Sciences (CAS); Shandong Provincial Key Laboratory of Coastal Zone Environmental Processes, YICCAS, Yantai Shandong, P. R. China, The Graduate School of Chinese Academy of Sciences, Beijing, P. R. China

  • Jianmin Zhao,

    Affiliation Key Laboratory of Coastal Zone Environmental Processes, Yantai Institute of Coastal Zone Research (YIC), Chinese Academy of Sciences (CAS); Shandong Provincial Key Laboratory of Coastal Zone Environmental Processes, YICCAS, Yantai Shandong, P. R. China

  • Junbao Yu

    Affiliation Key Laboratory of Coastal Zone Environmental Processes, Yantai Institute of Coastal Zone Research (YIC), Chinese Academy of Sciences (CAS); Shandong Provincial Key Laboratory of Coastal Zone Environmental Processes, YICCAS, Yantai Shandong, P. R. China

Abstract

As an environmental contaminant, mercury is of great concern due to its high risk to environmental and human health. The halophyte Suaeda salsa is the dominant plant in the intertidal zones of the Yellow River Delta (YRD) where has been contaminated by mercury in some places. This study aimed at evaluating the chronic effects of mercury (Hg2+, 20 µg L−1) and the influence of an environmentally relevant salinity (NaCl, 500 mM) on mercury-induced effects in S. salsa. A total of 43 protein spots with significant changes were identified in response to Hg2+, salinity and combined Hg2+ and salinity. These proteins can be categorized into diverse functional classes, related to metabolic processes, photosynthesis, stress response, protein fate, energy metabolism, signaling pathways and immunosuppression. Metabolic responses demonstrated that Hg2+ could disturb protein and energy metabolisms in S. salsa co-exposed with or without salinity. In addition, both antagonistic and synergistic effects between Hg2+ and salinity were confirmed by differential levels of proteins (magnesium-chelatase and ribulose-l,5-bisphosphate carboxylase/oxygenase) and metabolites (valine, malonate, asparagine, glycine, fructose and glucose) in S. salsa. These findings suggest that a combination of proteomics and metabolomics can provide insightful information of environmental contaminant-induced effects in plants at molecular levels.

Introduction

Mercury is a non-essential and highly toxic heavy metal to organisms [1]. It can be found in cells bound with thiol-containing proteins, glutathione (GSH) and cysteine. In addition, mercury can cause oxidative stresses by producing reactive oxygen species (ROS) such as hydrogen peroxide, superoxide radical and hydroxyl radical in plant cells and lipid peroxides and consequently damages protein, lipid and DNA [2]. In some intertidal zones of the Yellow River Delta (YRD), mercury is a severe heavy metal contaminant and has posed great risk to the intertidal organisms due to its high toxicity [3]. The mercury concentration has been as high as 66 mg kg−1 in the intertidal sediment in some mercury-polluted sites in the YRD [3]. In soils and waters, Hg (II) is the predominant form that can be easily accumulated in plants [1]. Realistically, however, the intertidal zones are of a high salinity which is one of the most common environmental stressors resulting in osmotic pressure, water deficit, ion toxicity and nutrient deficiency in plants [4]. Therefore salinity should be considered in environmental toxicology since it can potentially influence contaminant-induced effects in organisms.

The halophyte Suaeda salsa is the pioneer plant in the intertidal zones of the YRD, where salinity is often up to 3% [5]. This plant species is consumed as an edible vegetable by local residents because of its beneficial nutrient components such as vitamins, proteins and antioxidant ingredients [5]. S. salsa has also demonstrated the virtues in environmental sciences [6], [7]. For example, it has been applied in phyto-remediation for heavy metals and oil in saline soil [6], [7] and deemed as an environmental bioindicator for heavy metals in the intertidal zones because of its substantial tolerance to environmental salinity compared with animals [8]. In this study, we selected S. salsa as the experimental plant to evaluate the influence of environmental salinity on the toxicological effects of environmental relevant mercury in the intertidal zones in the YRD.

Proteomics and metabolomics are two well-established “–omic” techniques in the post-genomic era [9][11]. Proteomics is the study of large-scale proteins in an organism encoded by its genome [12]. Not only is proteomics a powerful tool for describing complete proteomes at organelle, cell, organ or tissue levels, but it can also be used to compare proteomes under various conditions of stress, such as those resulted from the exposure to heavy metals [13]. Metabolomics focuses on a global profile of the low molecular weight (<1 000 Da) metabolites which are the end products of metabolisms in biofluids, tissues and even whole organism [14], [15]. The characterization of endogenous metabolites can then provide the information of metabolic status in an organism to assess the biological responses induced by exogenous factors [16], [17]. Both proteomics and metabolomics are frequently used to characterize the perturbations in metabolic pathways and corresponding enzymes and stress-responsive proteins induced by exogenous factors [18][20]. Obviously, a combination of proteomics and metabolomics could potentially validate and complement one another, when testing the toxicological effects of environmental pollutants in organisms, especially in non-model organisms [21]. To our knowledge, however, no attempt has been made to evaluate the influence of environmental salinity on the toxicological effects of heavy metal contaminants in the intertidal zones using an integrated proteomic and metabolomic approach at molecular levels.

In this study, two “-omic” approaches, traditional 2-DE-based proteomics and proton nuclear magnetic resonance (NMR) spectroscopy-based metabolomics were used to investigate the chronic effects and corresponding proteomic and metabolic biomarkers induced by environmentally relevant mercury (Hg2+, 20 µg L−1) in S. salsa co-exposed without or with an environmentally relevant salinity (NaCl, 500 mM). In addition, the antioxidant status was assessed in S. salsa leaves. Taking all of these outcomes into account, we made an effort to answer three questions. Firstly, what effects can be induced by environmental relevant mercury in S. salsa on the basis of metabolic profiles, protein changes and antioxidant enzyme activities? Secondly, what influence can be caused by environmental salinity on the effects of mercury in S. salsa? Thirdly, are there any stable proteomic and/or metabolomic biomarkers of mercury-induced effects in S. salsa under the experimental salinity?

Materials and Methods

Plant Materials and Experimental Conditions

The seeds of S. salsa were collected from the wetland located in the Yellow River Delta near the Nature Reserve in November, 2010 and stored in a refrigerator at 4°C. No specific permissions were required for these locations/activities. The location is not privately-owned or protected in any way. The field studies also did not involve endangered or protected species. After surface sterilization in 0.5% HgCl2, forty seeds were washed in sterilized double distilled water for three times and sown in the sterilized sands in 4 replicate polyvinylchloride cylinders (20 cm diameter, 20 cm depth). The sown seedlings (n = 10) irrigated with Hoagland’s solution (containing 0.1% sodium chloride, pH = 6.0) were divided into four groups: (1) control, (2) Hg2+ (20 µg L−1), (3) NaCl (500 mM) and (4) combined Hg2+ (20 µg L−1) and NaCl (500 mM). The experimental concentrations of mercury and salinity were environmentally relevant to the real situation of polluted intertidal zones in the YRD [3], [22]. Mercury concentrations were ranged from 19.3 to 20.3 µg L−1 determined by an atomic fluorescence spectrometer (SK-2002A, Jinsuokun LTD., Beijing, China). The seedlings were cultivated at 28±4°C with the photoperiod of 12 h light/12 h darkness, relative humidity 70% and photo-synthetically active radiation 600 µmol m−2 s−1. After exposure for 30 days, the seedlings (n = 6) of S. salsa from both control and exposed groups were randomly harvested. After quick measure of fresh weight and total length of seedlings, all plant leaf tissues were collected and flash-frozen in liquid nitrogen and stored at −80°C before extraction for proteins, total RNA and metabolites and enzymatic assay.

Protein Extraction

Denaturing protein extraction was performed based on the method of Hajheidari et al and Saravanan et al [23], [24]. The frozen leaf tissue was homogenized in liquid nitrogen, resuspended in cold buffer (10% w/v trichloroacetic acid in acetone contained 0.07% w/v DTT). The samples were precipitated overnight at −20°C, centrifuged at 15 000 g for 30 min at 4°C, and the supernatant was discarded. The pellets were resuspended in cold acetone containing 0.07% w/v DTT, then kept at −20°C for 1 h and centrifuged at 15 000 g for 30 min at 4°C. This procedure was repeated twice. The pellets were solubilized in the lysis buffer (7 M urea; 2 M thio urea; 4% m/V CHAPS; 65 mM DTT and 0.2% W/V Bio-lyte buffer) and then incubated for 3 h at room temperature [25]. The homogenate was centrifuged at 15 000 g for 10 min and the supernatant was used for electrophoresis. The total protein was determined by Protein Assay Kit of Tiangen (Beijing, China) using bovine serum albumin as a standard.

Two-dimensional Gel Electrophoresis

The first dimension (IEF) was performed using the Immobiline Drystrip (24 cm, pH 4–7, linear). One hundred and seventy microgram (170 µg) of proteins was loaded. Isoelectric focusing gel solution containing 7 M urea, 2 M thio urea, 4% m/v CHAPS, 65 mM DTT, 0.001% m/v Bromophenol blue and 0.2% w/v Bio-lyte buffer. IEF was conducted at 20°C with an Etan IPGphor3 system for a total of 90 000 Vh. After the first dimension, strips were placed in the equilibration buffer (0.05 M Tris-HCl, pH 8.8; 6 M urea; 30% glycerol; 2% w/v SDS; containing 1% w/v DTT) and were slowly shaken for 13 min. Subsequently, the strips were incubated for 13 min in the same equilibration buffer with 2.5% (w/v) iodoacetamide without DTT [26]. The second dimension was performed on 10% SDS-PAGE gels using the Ettan DALTsix system (GE Healthcare). After electrophoresis, the gels were silver stained by following the method of Mortz et al and Gharahdaghi et al [27], [28]. Images were captured by ImageScanner III and spots were quantitatively analyzed by using ImageMaster 2D Platinum 7.0. The abundance of each protein spot was estimated by the percentage volume (% vol).

In Gel Digestion and MS Analysis

In gel digestion was performed according to Katayama et al [29]. Protein spots were washed three times with ultrapure water and de-stained with 25 mmol/L NH4HCO3 in 50% v/v acetonitrile at room temperature for 30 min. The gels were dried in 50% acetonitrile for 30 min and 100% acetonitrile for another 30 min. The samples were rehydrated in 10 µL cover solution (0.02 g L−1 w/v trypsin, 25 mM NH4HCO3 and 10% acetonitrile) for 30 min, and then covered with the same solution but without trypsin for digestion overnight at 37°C. The supernatants were extracted with 5% TFA in 67% acetonitrile at 37°C for 30 min, then centrifuged at 5 000 g for 5 min, so the peptide extracts and the supernatant of the gel spot were combined.

After being completely dried, the samples were resuspended with 5 µL 0.1% TFA followed by mixing in 1∶1 ratio with a saturated solution of α-cyano-4-hydroxy-trans-cinnamic acid in 50% acetonitrile [30]. One microliter of mixture was analyzed by an ABI 4800 MALDI-TOF/TOF Plus mass spectrometer (Applied Biosystems, Foster City, USA), and data were acquired in a positive MS reflector using a CalMix5 standard to calibrate the instrument (ABI4800 Calibration Mixture). Both the MS and MS/MS data were integrated and processed using the GPS Explorer V3.6 software (Applied Biosystems, USA) with default parameters. Proteins were successfully identified based on 95% or higher confidence interval of their scores in the MASCOT V2.4 search engine (Matrix Science Ltd., London, U.K.). The following parameters were used in the search: NCBInr Viridiplantae (Green Plants) database; trypsin as the digestion enzyme; one missed cleavage site; partial modifications of cysteine carbamidomethylation and methionine oxidization; no fixed modifications; 0.15 Da for precursor ion tolerance and 0.25 Da for fragment ion tolerance. Individual ions scores >40 indicated identity or extensive homology (P<0.05).

Metabolite Extraction

Polar metabolites were extracted from leaf tissues using an extraction system of methanol/water (1/1) as described previously [31]. Briefly, the tissue sample (approximately 300 mg) was ground in a liquid N2-cooled mortar and pestle. The tissue powder was thoroughly homogenized in 3.33 mL g−1 methanol/water (1/1) using a high throughput homogenizer, Precellys 24 (Bertin, France). After homogenization, the sample was transferred to an Eppendorf tube and centrifuged for 10 minutes (3 000 g, 4°C). The supernatant was removed and then lyophilized. It was subsequently re-dissolved in 600 µL phosphate buffer (0.1 M Na2HPO4 and NaH2PO4, containing 0.5 mM TSP, pH 7.0) in D2O. The mixture was vortexed and then centrifuged at 3 000 g for 5 min at 4°C. The supernatant substance (550 µL) was pipetted into a 5 mm NMR tube for 1H NMR spectroscopic analysis.

NMR Spectroscopy

Tissue extracts of S. salsa were analyzed on a Bruker AV 500 NMR spectrometer at 500.18 MHz (at 298 K) [31]. One dimensional (1D) 1H NMR spectra were obtained using a 11.9 µs pulse, 6009.6 Hz spectral width, mixing time 0.1 s, and 3.0 s relaxation delay with standard 1D NOESY pulse sequence, with 128 transients collected into 16 384 data points. Datasets were zero-filled to 32 768 points, and exponential line-broadenings of 0.3 Hz were applied before Fourier transformation. All 1H NMR spectra were phased, baseline-corrected, and calibrated (TSP at 0.0 ppm) manually using TopSpin (version 2.1, Bruker). The metabolites were identified and quantified following tabulated chemical shifts and using the software, Chenomx (Evaluation Version, Chenomx Inc., Canada) [32]. The metabolite concentrations were normalized to the mass of leaf tissue by calculating the concentration of metabolites in each NMR tube.

Data Pre-processing and Multivariate Analysis

One dimensional proton NMR spectra were converted to a format for multivariate analysis using custom-written ProMetab software in Matlab (version 7.0; The MathsWorks, Natick, MA) [33]. Each spectrum was segmented into 0.005 ppm bins between 0.2 and 10.0 ppm with bins from 4.70 to 5.20 ppm (water) excluded from all the NMR spectra. The total spectral area of the remaining bins was normalized to unity that could facilitate the comparison between the spectra. All the NMR spectra were generalized log transformed (glog) with transformation parameter λ = 1.0×10−9 [33] to stabilize the variance across the spectral bins and to increase the weightings of the less intense peaks. Data were mean-centered before principal components analysis (PCA).

PCA is an exploratory unsupervised pattern recognition method of analysis which is blind to the status of each sample, and serves to reduce the dimensionality of the data and summarize the similarities and differences between control and treatment groups [8]. The algorithm of this pattern recognition method calculates the highest amount of correlated variation along PC1, with subsequent PCs containing correspondingly smaller amounts of variance. For each model built, the loading vector for the PC was examined to identify the metabolites which contributed to the clusters. For proteomic data, PCA was conducted on the percentage volume (% vol) of common protein spots (matched protein spots in all 2-DE gels) to discover the proteomic differences between different treatments using software PLS Toolbox (version 4.0, Eigenvector Research, Manson, WA). One way analysis of variance (ANOVA) was conducted on the PC scores from each group to test the statistical significance (P<0.05) of separations. SAM software was then used to find significant metabolic differences from various groups with appropriate false discovery rate (FDR <0.01) cutoffs [34]. Especially, only those proteins with significant changes (>2.0 folds, P<0.05) were considered to be differentially expressed proteins.

Measurement of Antioxidant Enzyme Activities

The antioxidant enzymes for the activity measurement included superoxide dismutase (SOD, EC 1.15.1.1), peroxidase (POD, EC 1.11.1.7), and catalase (CAT, EC 1.11.1.6). The leaf tissue of S. salsa was homogenized in cold buffer (10% w/v phosphate buffer, pH = 7.8), and the antioxidant enzyme activities were assayed by a Multiskan spectrum microplate spectrophotometer (Infinite M200, TECAN) according to the manufacturer’s protocols using the enzyme kits (Jiancheng, Nanjing, China). All the enzyme activities were expressed as units per milligram protein (U/mg protein).

Total mRNA Extraction and Gene Quantification

The expression of the housekeeping genes (Table 1) in S. salsa was determined by qRT-PCR. The data were analyzed with geNorm to calculate the expression stability (M values) of potential reference genes required for accurate normalization (V values) [35]. GeNorm identified β-actin as the most stable gene with a V2/3 value 0.139 less than the proposed geNorm cutoff value of 0.15, which was then followed by 18S ribosomal RNA, helicase, polyubiquitin, ubiquitin, histone. Therefore β-actin was used as the internal control for gene expression normalization.

thumbnail
Table 1. The list of primers used for the determination of internal control and quantification of gene expressions by qRT-PCR.

https://doi.org/10.1371/journal.pone.0064041.t001

Total RNA from leaf tissues of S. salsa was isolated using the Trizol reagent (Invitrogen). Single-strand cDNA was synthesized from total RNA with M-MLV reverse transcriptase (Promega, USA). Gene-specific primers for heat shock protein 70, aspartate aminotransferase, S-adenosylmethionine synthase, transketolase, peroxiredoxin, formate dehydrogenase, glutamine synthetase GS2, ribulose-1,5-bisphosphate carboxylase/oxygenase (rubisco), chlorophyll A/B binding protein, carbonic anhydrase and ATP synthase delta chain were used to amplify amplicons specific for S. salsa samples. The sequences of primers are given in Table 1.

The quantitative real time RT-PCR amplifications were carried out in triplicate with a total volume of 50 µL containing 25 µL of 2×SYBR Green PCR Master Mix, 20 µL of the diluted cDNA, 1 µL of each of primers (10 µmol/L), and 3 µL of DEPC-treated water. The quantitative real time RT-PCR program was as follows: 50°C for 2 min, 95°C for 10 min followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. After PCR program, the data were analyzed with the ABI 7500 SDS software (Applied Biosystems). The comparative CT method (2−ΔΔCT) was used to analyze the expression levels of genes [36].

Statistical Analysis

Data of growth parameters, enzyme activities, gene expressions and metabolite concentrations were expressed as the mean ± standard deviation (n = 6). The Minitab software (Version 15, Minitab Inc., USA) was used for the statistical analysis. One way analysis of variance (ANOVA) with Tukey’s test was conducted on the data, and a P value at 0.05 was considered significant.

Results

The Growth Data of S. salsa from Control and Treatment Groups

The average lengths and fresh weights of S. salsa seedlings were shown in Figure 1. Obviously Hg2+ exposure (20 µg L−1) did not significantly (P>0.05) affect the growth of S. salsa. However, both salinity and Hg2++salinity treatments significantly inhibited the growth of S. salsa after exposure for 1 month.

thumbnail
Figure 1. The average length (A) and fresh weight (B) of S. salsa seedlings from control, Hg2+−, salinity- and Hg2++salinity-treated groups.

https://doi.org/10.1371/journal.pone.0064041.g001

Proteomic Responses in Leaf Tissues of S. salsa

Comparative proteomic analysis was used to determine the proteomic profiles in S. salsa exposed to Hg2+, salinity and combined Hg2+ and salinity (designated as Hg2++salinity). A total of 60 spots resolved in 2-DE gels were differentially expressed (>2 folds, P<0.05). Figure 2 represented the differential protein spots in S. salsa with different exposures. In the present study, the proteins spots observed in all three biological replicates of silver-stained gels were analyzed by MALDI-TOF/TOF mass spectrometry and 43 (∼70%) proteins were successfully identified. The details of matched proteins including their protein names, accession numbers and other MS data are summarized in Table 2.

thumbnail
Figure 2. Representative 2-DE images of proteins extracted from S. salsa leaves.

Proteins were submitted to isoelectric focusing on 4–7 IPG strips (24 cm) followed by electrophoresis on 10% SDS-PAGE. Gels were stained by silver stain. Gels are from the exposed groups of control (A), Hg2+ (20 µg L−1) (B), salinity (NaCl, 500 mM) (C) and Hg2++salinity (D).

https://doi.org/10.1371/journal.pone.0064041.g002

thumbnail
Table 2. List of protein spots which were differentially expressed in S. salsa leaves under mercury, salinity treatment after 30 days.

https://doi.org/10.1371/journal.pone.0064041.t002

In Hg2+-treated group, a total of 17 differentially expressed spots were discovered, including 4 up-regulated and 13 down-regulated. These proteins were related to metabolism, protein biosynthesis, stress and defense, signaling pathway and energy metabolism (Table 2). Nine differentially expressed proteins were observed in salinity-treated S. salsa samples, including 3 up-regulated and 6 down-regulated proteins that were involved in metabolism, photosynthesis, energy metabolism and stress and defense. For the Hg2++salinity-treated group, a total of 22 proteins (7 up-regulated and 15 down-regulated) were significantly altered. All these differentially expressed proteins from different treatments were basically involved in metabolism, photosynthesis, protein folding and assembly, protein biosynthesis, stress and defense, signaling pathway, immunosuppression and energy metabolism. Figure 3 revealed the total numbers of protein spots differentially expressed from S. salsa samples with various treatments.

thumbnail
Figure 3. Number of protein spots from S. salsa samples (A) and the differentially expressed proteins in S. salsa from Hg2+ (20 µg L−1), salinity (NaCl, 500 mM) and Hg2++salinity (B).

https://doi.org/10.1371/journal.pone.0064041.g003

The Relationship between Gene Expressions and Protein Abundances

To further verify the results of 2-DE and compare the correlation between gene expression and protein abundances, eleven representative genes related to regulated proteins were quantified using quantitative real-time RT-PCR technique. The results indicated that the expression levels of the genes had similar alteration tendency with corresponding proteins (Figure 4).

thumbnail
Figure 4. Expression levels of eleven genes in leaf tissues of Suaeda salsa from control, Hg2+ (20 µg L−1), salinity (NaCl, 500 mM) and Hg2++salinity-treated groups after exposure for 30 days.

Data (n = 6) were expressed as mean ± standard deviation. Significant difference between control and exposed groups was tested by one-way analysis of variance with Tukey’s test and indicated by *(P<0.05) and **(P<0.01).

https://doi.org/10.1371/journal.pone.0064041.g004

Metabolic Responses in Leaf Tissues of S. salsa

Figure 5 shows one representative 1H NMR spectrum of leaf tissue extracts from the control group. Different classes of metabolites were identified in leaf tissue of S. salsa, including amino acids (branched chain amino acids: valine, leucine and isoleucine, alanine, threonine, glutamate, aspartate, asparagine, glycine, etc.), sugars (fructose and glucose), an organic osmolyte (betaine), organic acids (acetate and malonate) and intermediates in the tricarboxylic acid (TCA) cycle (succinate, malate and fumarate). Especially, all NMR spectra were dominated by the osmolyte, betaine that is synthesized from choline by choline monooxygenase [37]. In the higher plant, halophyte S. salsa, betaine is the most important secondary metabolite playing a pivotal role in the protection from osmotic stress [38].

thumbnail
Figure 5. One representative 1-dimensional 500 MHz 1H NMR spectrum of metabolite extracts of the leaf tissue from S. salsa in control group using extraction solvent system of methanol/water (1/1).

Metabolite assignments: (1) branched chain amino acids: leucine, isoleucine and valine, (2) ethanol, (3) threonine, (4) alanine, (5) acetate, (6) glutamate, (7) glutamine, (8) succinate, (9) malate, (10) aspartate, (11) asparagine, (12) malonate, (13) choline, (14) phosphocholine, (15) betaine, (16) glycine, (17) fructose, (18) glucose, (19) unknown (5.39 ppm) and (20) fumarate.

https://doi.org/10.1371/journal.pone.0064041.g005

PCA resulted in significant separation (P<0.05) control and each treatment (data not shown), which indicated that Hg2+, salinity and Hg2++salinity induced significant effects in S. salsa. Especially, the inhibited plant growth data confirmed the effects by salinity and Hg2++salinity exposures (Figure 1). For further analysis, one way ANOVA was directly conducted on the metabolite concentrations to explore the significant metabolic responses induced by these treatments. For Hg2+-treated S. salsa group, the metabolic profiles of leaf tissue extracts indicated significant (P<0.05) increases of branched chain amino acids (valine, leucine and isoleucine) and phosphocholine, and significant (P<0.05) decreases of ethanol, acetate, succinate, malonate and choline (Table 3). The significant metabolic differences induced by salinity in the leaf tissues of S. salsa included the depleted threonine, alanine and acetate, and elevated malate, choline, betaine, fructose and glucose (Table 3). For Hg2++salinity-treated S. salsa samples, the significant metabolic changes resulted in decreases of ethanol, acetate, succinate, malonate and glycine, and increases of valine, asparagine, phosphocholine and betaine.

thumbnail
Table 3. Metabolite concentrations (µmol g−1 wet tissue) in leaf tissues of S. salsa exposed to Hg2+ (20 µg L−1), salinity (NaCl, 500 mM) and combined Hg2+ and salinity.

https://doi.org/10.1371/journal.pone.0064041.t003

Antioxidant Enzyme Activities

After exposure for 30 days, the average activities of antioxidant enzymes including SOD, POD and CAT were significantly (P<0.05) increased in both salinity and Hg2++salinity-exposed samples of S. salsa (Figure 6). However, Hg2+ treatment induced significant (P<0.05) decreases in antioxidant enzyme activities in the S. salsa samples.

thumbnail
Figure 6. The activities (U mg−1 protein) of SOD (A), POD (B) and CAT (C) in leaf tissues from control, Hg2+ (20 µg L−1), salinity (NaCl, 500) and Hg2++salinity-exposed S. salsa samples (n = 6).

Data were expressed as mean ± standard deviation. Significant difference between control and exposed groups was tested by one-way analysis of variance with Tukey’s test and indicated by *(P<0.05) and **(P<0.01).

https://doi.org/10.1371/journal.pone.0064041.g006

Discussion

Effects of Hg2+, Salinity and Hg2++Salinity on the Proteome of S. salsa

Figure 7A summarized the pathways involved in the response of S. salsa to Hg2+. Some proteins related to glycolysis, pentose phosphate pathway, fatty acid synthesis, TCA cycle were changed (decreased or increased) in Hg2+-treated S. salsa group (Table 2), which clearly demonstrated the disturbance in metabolism caused by Hg2+ exposure in S. salsa [39][45]. One 14-3-3 protein was decreased (>4 folds) in mercury-treated group. 14-3-3 proteins in plant cells participate in signal transduction, regulate the activity of enzymes, and are associated with the ATP synthases, representing a mechanism for plant adaptation to environmental changes [46][49]. Especially, ATP synthase CF1 epsilon subunit was dramatically decreased (>34 folds) in Hg2+-treated S. salsa group. Bunney et al. (2001) reported that 14-3-3 protein is a regulator of the mitochondrial and chloroplast ATP synthase (Figure 7A) [46], therefore, the consistent down-regulation of 14-3-3 protein and ATP synthase CF1 epsilon subunit confirmed the disturbance in energy metabolism caused by Hg2+ exposure in S. salsa. Salt tolerance proteins can increase salt tolerance in some plants and can be over expressed when plants are exposed to salinities [49]. Interestingly, one salt tolerance protein was significantly down-regulated in Hg2+-treated S. salsa samples, which might suggest a contrary effect of Hg2+ to salinity in S. salsa. Heat shock protein 70 (Hsp 70) has multiple functions in preventing aggregation, assisting refolding, protein import and translocation, signal transduction, anti-oxidative stress and transcriptional activation under different environmental stress conditions such as heat, cold and drought [50]. Mercury can induce oxidative stress in organisms. Therefore, S. salsa may use various protective antioxidants to resist oxidative stress and probably up-regulates the contents of antioxidants, such as heat shock proteins and anti-oxidative enzyme activities, such as SOD and CAT [51]. Elbaz et al. (2010) reported that mercury (>200 µg L−1) could induce significant oxidative stress in Chlamydomonas reinhardtii [51]. However, Hsp 70 was down-regulated under Hg2+ exposure, together with decreased SOD, CAT and POD activities which were significantly increased in S. salsa. Clearly, the environmentally relevant concentration (20 µg L−1) of Hg2+ induced a contrary result in S. salsa compared with that in Chlamydomonas reinhardtii Hg2+ (>200 µg L−1), which probably implied that this relatively low concentration (20 µg L−1) of Hg2+ induced a hormesis phenomenon of antioxidants (Hsp 70, SOD, CAT and POD) in S. salsa. In both salinity- and Hg2++salinity-treated groups, these antioxidant enzymes were significantly increased which apparently indicated the oxidative stress induced by salinity in S. salsa but Hg2+.

thumbnail
Figure 7. Schematic presentation of molecular responsive-mechanisms in S. salsa to Hg2+ co-exposed without (A) or with (B) the environmentally relevant salinity (NaCl, 500 mM) according to Kyoto Encyclopedia of Genes and Genomes (http://www.genome.jp/kegg/) and Uniprot (http://www.uniprot.org/).

The identified proteins and metabolites were shown by marking the protein names in red (up-regulated) or blue (down-regulated).

https://doi.org/10.1371/journal.pone.0064041.g007

In response to salinity, the identified proteins relating to metabolism were involved in pentose phosphate pathway (transketolase), polyamine biosynthesis (spermidine synthase), amino-acid biosynthesis (cysteine synthase) and isoprenoid biosynthetic process (1-deoxy-D-xylulose 5-phosphate reductoisomerase). The down-regulated ATP synthase CF1 alpha subunit and ATPase subunit implied the disturbance in energy metabolism caused by salinity exposure in S. salsa [52]. The similar result of down-regulated ATP synthase was also observed in another halophyte Thellungiella exposed to salinities, which could imply the similar responsive mechanism of energy metabolism halophytes to salinity exposures [52].

Figure 7B summarized the proteomic and metabolomic responses involved in the pathways, including Calvin cycle, ROS scavenging, TCA cycle, chlorophyll synthesis, nucleotide metabolism, immunosuppression, signaling pathway and energy metabolism, in S. salsa to Hg2++salinity. The proteins associated with metabolism were involved in Calvin cycle, pyruvate metabolism, sulfur metabolism, cysteine and methionine metabolism, methane metabolism, glutamate metabolism, pantothenate and CoA biosynthesi [53][57], as shown in Table 2 and Figure 7B. Auxins play a vital role in division and expansion of plant cells, and its action is mediated by the putative receptor called auxin-binding protein [58] which was dramatically down-regulated (∼9 folds) in Hg2++salinity-treated samples. Sauer et al. (2006) found that auxin-binding protein was significantly up-regulated in maize during shoot-borne root initiation, which confirmed that auxin-binding protein was important in plant growth [59]. Therefore, the down-regulation of auxin-binding protein indicated that the co-exposure of Hg2+ and salinity inhibited plant growth in S. salsa, which was indicated by the reduced (P<0.001) plant weight and length of Hg2++salinity-treated samples (Figure 1). Of particular interest is that cyclophilin was found to be decreased under Hg2++salinity stress. Cyclophilin belongs to a family of immunosuppressant receptor proteins presented in all subcellular compartments which are involved in a wide variety of processes such as protein trafficking and maturation and receptor complex stabilization [60]. The down-regulation of cyclophilin probably indicated that multiple cellular pathways were affected by Hg2++salinity exposure. Peroxiredoxin, an antioxidant enzyme reducing hydrogen peroxide (H2O2) and alkyl hydroperoxides treatments, regulates peroxide-mediated signaling cascades [61], [62]. In another native halophyte Thellungiella in the YRD, peroxiredoxin was down-regulated the seedlings of Thellungiella was exposed to NaCl (150 mM) [52]. In our case, the single exposure of Hg2+ or salinity did not induce obvious alteration of peroxiredoxin in S. salsa, however, the combined exposure of Hg2+ and salinity significantly down-regulated peroxiredoxin. It might be the synergistic effects between salinity and Hg2+ in S. salsa. As to protein folding and assembly, chloroplast RNA binding protein and protein disulfide isomerase (PDI) were down-regulated. PDI is a multifunctional enzyme catalyzes a set of disulfide-exchange reaction that includes formation, reduction or isomerization of protein disulfide bonds [63], [64]. In addition, as a member of the superfamily of thioredoxin, PDI has been found in response to oxidative stress [65]. Proteins are the main nonwater components of a cell and they are the major targets attacked by intracellular reactive oxygen species (ROS). As a member of thioredoxin family, PDI can reverse these oxidative modifications. Therefore, the down-regulation of PDI could imply that the Hg2++salinity affected the antioxidant system in S. salsa, which is evidenced by the increased activities of CAT, SOD and POD in this treatment.

50S ribosomal protein L12 was down-regulated in both Hg2+- and Hg2++salinity-treated groups. The synthesis of this protein occurs in large macromolecular particles called ribosomes [66]. And it is presumed to be involved in the binding of translation factors, stimulating factor-dependent GTP hydrolysis [67]. Transketolase was down-regulated and dehydroascorbate reductase-like protein was up-regulated in both Hg2+- and salinity-treated group. Dehydroascorbate reductase (DHAR) plays an important role in maintaining the normal level of ascorbic acid by recycling oxidized ascorbic acid [68], [69]. The increasing dehydroascorbate reductase-like protein is critical for protection of cellular components against mercury and salinity stresses. Proteins linked to the function of photosynthesis (magnesium-chelatase and ribulose-l,5-bisphosphate carboxylase/oxygenase) were decreased in both salinity- and Hg2++salinity-treated groups. Magnesium-chelatase catalyzes the insertion of Mg2+ into protoporphyrin IX, the first dedicated step in chlorophyll biosynthesis [70], [71]. Ribulose-l,5-bisphosphate carboxylase/oxygenase acts as a carboxylase with the substrates RuBP and CO2 and as an oxygenase with the substrates RuBP and O2. The degradation of Rubisco can be a key factor that determined the photosynthetic activity in plants [72], [73]. The down-regulation of magnesium-chelatase and ribulose-l,5-bisphosphate carboxylase/oxygenase demonstrated that Hg2+ and salinity stresses could weaken the photosynthesis system of S. salsa. Interestingly, the level of magnesium-chelatase in salinity-treated group was significantly lower than that in Hg2++salinity-treated group, which suggested the antagonistic effect between salinity and Hg2+ in S. salsa. Compared to Hg2++salinity-treated group, however, the relatively high level of ribulose-l,5-bisphosphate carboxylase/oxygenase in salinity-treated group implied the additive effects of Hg2+ to salinity.

Essentially, gene expression means the tendency of the corresponding encoded protein that is the down-stream product of the gene. In some cases, however, gene expressions are not always correlated with the protein abundances. The lack of correlation between gene and protein expressions is not strange since there could be different regulatory mechanisms of gene expressions and protein expressions to contaminant exposures. In addition, the protein translation of mRNA transcription does not always happen due to the posttranscriptional modifications. In our case, the expression levels of selected genes were consistent with alteration tendency of corresponding proteins, which could partly confirm the regulations of proteins in S. salsa exposed to Hg2+, salinity and Hg2++salinity.

Effects of Hg2+, Salinity and Hg2++Salinity on the Metabolome of S. salsa

For the leaf tissues of S. salsa seedlings exposed to Hg2+, the levels of branched chain amino acids (valine, isoleucine and leucine) were approximately 3 times higher than those in the samples of control group (Table 3), which was probably associated with enhanced protein bio-degradation or inhibited protein biosynthesis caused by Hg treatment in S. salsa [19]. Ethanol is the end product of anaerobic metabolism through fermentation pathway. Subsequently, ethanol can be catalyzed into acetate by alcohol dehydrogenase and aldehyde dehydrogenase [74]. In this study, both ethanol and acetate were significantly (P<0.05) decreased in Hg2+-exposed S. salsa samples, which indicated that Hg2+ exposure inhibited anaerobic metabolism in the leaf of S. salsa. Succinate, one of the intermediates involved in tricarboxylic acid (TCA) cycle, was significantly (P<0.01) depleted in S. salsa samples, indicating the disturbance in TCA cycle that is related to energy metabolism [19]. Interestingly, succinyl-CoA-synthetase that is involved in TCA cycle, was down-regulated as well in Hg2+-exposed S. salsa samples. Succinyl-CoA-synthetase is an enzyme that catalyzes the reversible reaction of succinyl-CoA to succinate (as shown in Figure 7A). Therefore the consistent alteration of both succinyl-CoA and succinate confirmed that proteomics and metabolomics might validate one another in special metabolic pathways. Choline and ATP can be converted into phosphocholine and ADP catalyzed by choline kinase. Therefore the elevated phosphocholine and reduced choline might indicate the enhanced conversion of choline to phosphocholine and hence enhanced energy cost (ATP→ADP) to prevent S. salsa from Hg2+-induced damage.

Salinity is a common environmental stressor to plants. Although S. salsa is a halophyte that can tolerate high salinities, the extremely high concentration (∼500 mM) of environmental salinity can lead to significant inhibition of plant growth, which was observed in this study (data not shown). At metabolite level, the experimental concentration of salinity (500 mM, also realistically relevant) induced differential metabolic responses in leaf tissues of S. salsa compared with those induced by Hg2+ exposure. The decreased threonine and alanine indicated possible inhibition of protein bio-degradation or enhanced protein biosynthesis caused by salinity stress in S. salsa. The significant decrease of acetate could be related to the inhibited anaerobic metabolism combined with the reduced level of ethanol (although not significant, P = 0.11). The intermediate in TCA cycle, malate, was significantly (P<0.05) increased, which suggested the disturbance in energy metabolism. Both choline and betaine were significantly increased in salinity-stressed S. salsa leaf tissues. The pathway of betaine synthesis is catalyzed by choline monooxygenase converting choline to betaine aldehyde that is subsequently converted into betaine by betaine aldehyde dehydrogenase [37], [75]. Sugars, such as glucose and fructose, are produced from gluconeogenesis and photosynthesis. In plants, sugars are also compatible osmolytes to maintain osmotic balance between plant and environment. Therefore the elevated choline, betaine, fructose and glucose were the indicators of osmotic stress induced by salinity in S. salsa.

In realistic situation, S. salsa seedlings are exposed to Hg and salinity in a combined manner. After exposure of combined Hg and salinity for 30 days, several metabolites including valine, ethanol, succinate, malonate and phosphocholine were consistently altered with the metabolic responses from Hg-treated S. salsa samples, which demonstrated that these metabolic biomarkers induced by Hg were relatively stable under the experimental salinity (NaCl, 500 mM). However, other metabolic responses including isoleucine, leucine, threonine, alanine, choline, malate, fructose and glucose were at control level. It apparently confirmed the antagonistic effects between Hg and salinity. However, asparagine and glycine were uniquely altered in Hg2++salinity-exposed S. salsa samples, which suggested that there were synergistic effects between Hg and salinity as well. All these findings indicated that salinity could influence Hg-induced effects in S. salsa.

Conclusions

In this study, the protein changes, metabolic profiles and antioxidant enzyme activities were characterized in the halophyte Suaeda salsa with environmentally relevant exposures of Hg2+ (20 µg L−1), salinity (NaCl, 500 mM) and combined Hg2+ (20 µg L−1) and salinity (NaCl, 500 mM) for 30 days. Both altered proteins and metabolites of Hg2+ exposure were different from those of Hg2++salinity exposure in S. salsa. It clearly demonstrated that environmental salinity could influence the chronic effects of mercury in S. salsa. Figure 7 draw a broader picture of the effects characterized by both proteomics and metabolomics in S. salsa with the exposures of Hg2+ and combined Hg2+ and salinity. In addition, both proteomic and metabolomic responses indicated the antagonistic and synergistic effects between Hg2+ and salinity. Interestingly, six metabolites including valine, ethanol, acetate, succinate, malonate and phosphocholine and one protein, 50S ribosomal protein L12, were consistently regulated in both Hg2+ and Hg2++salinity-exposed S. salsa samples. These responses could potentially be used as biomarkers of mercury in S. salsa, even under an extreme environmental salinity (∼500 mM NaCl). Our results suggest that a combination of proteomics and metabolomics can provide more insightful information of environmental contaminant-induced effects in plants at molecular levels.

Author Contributions

Conceived and designed the experiments: HW. Performed the experiments: XL CJ LW. Analyzed the data: XL HW JZ. Contributed reagents/materials/analysis tools: HW JZ JY. Wrote the paper: XL.

References

  1. 1. Clarkson TW (1997) The toxicology of mercury. Critical Reviews in Clinical Laboratory Sciences 34: 369–403.
  2. 2. Patra M, Bhowmik N, Bandopadhyay B, Sharma A (2004) Comparison of mercury, lead and arsenic with respect to genotoxic effects on plant systems and the development of genetic tolerance. Environmental and Experimental Botany 52: 199–223.
  3. 3. Zhang X (2001) Investigation of pollution of Pb, Cd, Hg, As in sea water and deposit of the Bohai Sea area (In Chinese). Heilongjiang Environmental Journal 25: 87–90.
  4. 4. Nazarbeygi E, Yazdi HL, Naseri R, Soleimani R (2011) The Effects of different levels of salinity on proline and A-, B-chlorophylls in Canola. American-Eurasian Journal of Agricultural and Environmental Science 10: 70–74.
  5. 5. Zhao HL, Ma YZ, Li JX, Piao L (2010) A study on edible value of Suaeda salsa (In Chinese). Journal Anhui Agricultural Science 38: 14350–14351.
  6. 6. Li X, Xiang XM, Zhou JT, Qu BC, Zhou XB (2007) Advance of study on soil remediation and sewage treatment by halophyte Suaede salsa. Jiangsu Environmental Science and Technology 20: 53–54.
  7. 7. Xu CY, Liu XB, Liu ZG, Wang J, Jiang ZP, et al. (2007) Remedial effect of Suaeda salsa Pall planting on the oil polluted coastal zones. Journal of Safety and Environment 1: 37–39.
  8. 8. Liu X, Yang C, Zhang L, Li L, Liu S, et al. (2011) Metabolic profiling of cadmium-induced effects in one pioneer intertidal halophyte Suaeda salsa by NMR-based metabolomics. Ecotoxicology 20: 1422–1431.
  9. 9. Fernandez-Garcia N, Hernandez M, Casado-Vela J, Bru R, Elortza F, et al. (2011) Changes to the proteome and targeted metabolites of xylem sap in Brassica oleracea in response to salt stress. Plant, Cell & Environment 34: 821–836.
  10. 10. Galindo Gonzalez LM, Kayal WE, Ju CJ-T, Allen CCG, King-Jones S, et al. (2012) Integrated transcriptomic and proteomic profiling of whites pruce stems during the transition from active growth to dormancy. Plant, Cell & Environment 35: 682–701.
  11. 11. Wu H, Wang W-X (2010) NMR-based metabolomic studies on the toxicological effects of cadmium and copper on green mussels Perna viridis. Aquatic Toxicology 100: 339–345.
  12. 12. Anderson NL, Anderson NG (1998) Proteome and proteomics: new technologies, new concepts, and new words. Electrophoresis 19: 1853–1861.
  13. 13. Ahsan N, Renaut J, Komatsu S (2009) Recent developments in the application of proteomics to the analysis of plant responses to heavy metals. Proteomics 9: 2602–2621.
  14. 14. Lindon JC, Nicholson JK, Everett JR (1999) NMR spectroscopy of biofluid. Annual Reports on NMR Spectroscopy 38: 1–88.
  15. 15. Davis B (2005) ‘Growing pains for metabolomics’. The Scientist 19: 25–28.
  16. 16. Pedras MSC, Zheng QA (2010) Metabolic responses of Thellungiella halophila/salsuginea to biotic and abiotic stresses: Metabolite profiles and quantitative analyses. Phytochemistry 71: 581–589.
  17. 17. Wu H, Zhang X, Li X, Li Z, Wu Y, et al. (2005) Comparison of metabolic profiles from serum from hepatotoxin-treated rats by nuclear-magnetic-resonance-spectroscopy-based metabonomic analysis. Analytical Biochemistry 340: 99–105.
  18. 18. Poynton HC, Taylor NS, Hicks J, Colson K, Chan S, et al. (2011) Metabolomics of microliter hemolymph samples enables an improved understanding of the combined metabolic and transcriptional responses of Daphnia magna to cadmium. Environmental Science and Technology 45: 3710–3717.
  19. 19. Zhang JT, Zhang Y, Du YY, Chen SY, Tang HR (2011) Dynamic metabonomic responses of tobacco (Nicotiana tabacum) plants to salt stress. Journal of Proteome Research 10: 1904–1914.
  20. 20. Zhang W, Liu Y, Zhang H, Dai J (2012) Proteomic analysis of male zebrafish livers chronically exposed to perfluorononanoic acid. Environment International 42: 20–30.
  21. 21. Zhang H, Ding L, Fang X, Shi Z, Zhang Y, et al. (2011) Biological responses to perfluorododecanoic acid exposure in rat kidneys as determined by integrated proteomic and metabonomic studies. PLoS One 6: e20862.
  22. 22. Li R, Li H, Li Y, Zhang S, Qian Z, et al. (2001) Study of the heavy metals, nitrogen and phosphorus contaminants in the sediments of the Yellow River Delta (In Chinese). Acta Sedimentologica Sinica 19: 622–629.
  23. 23. Hajheidari M, Abdollahian-Noghabi M, Askari H, Heidari M, Sadeghian SY, et al. (2005) Proteome analysis of sugar beet leaves under drought stress. Proteomics 5: 950–960.
  24. 24. Saravanan RS, Rose JK (2004) A critical evaluation of sample extraction techniques for enhanced proteomic analysis of recalcitrant plant tissues. Proteomics 4: 2522–2532.
  25. 25. Tullius MV, Phillips NJ, Scheffler NK, Samuels NM, Munson JRS, et al. (2002) The lbgAB gene cluster of haemophilus ducreyi encodes a β-1,4-galactosyltransferase and an α-1,6-DD-heptosyltransferase involved in lipooligosaccharide biosynthesis. Infection and Immunity 70: 2853–2861.
  26. 26. Lee SC, Chan J, Clement M-V, Pervaiz Dr S (2006) Functional proteomics of resveratrol-induced colon cancer cell apoptosis: Caspase-6-mediated cleavage of lamin A is a major signaling loop. Proteomics 6: 2386–2394.
  27. 27. Mortz E, Krogh TN, Vorum H, Görg A (2001) Improved silver staining protocols for high sensitivity protein identification using matrix-assisted laser desorption/ionization-time of flight analysis. Proteomics 1: 1359–1363.
  28. 28. Gharahdaghi F, Weinberg CR, Meagher DA, Imai BS, Mische SM (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.
  29. 29. Katayama H, Nagasu T, Oda Y (2001) Inprovement of in-gel digestion protocal for peptide mass fingerprinting by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Rapid Communications in Mass Spectrometry 15: 1416–1421.
  30. 30. Shevchenko A, Wilm M, Vorm O, Mann M (1996) Mass spectrometric sequencing of proteins from silver-stained polyacrylamide gels. Analytical Chemistry 68: 850–858.
  31. 31. Liu X, Zhang L, You L, Wu H, Zhao J, et al. (2011) Metabolomic study on the halophyte Suaeda salsa in the Yellow River Delta. Clean – Soil, Air, Water 39: 720–727.
  32. 32. Fan WMT (1996) Metabolite profiling by one- and two-dimensional NMR analysis of complex mixtures. Progress in Nuclear Magnetic Resonance Spectroscopy 28: 161–219.
  33. 33. Purohit PV, Rocke DM, Viant MR, Woodruff DL (2004) Discrimination models using variance-stabilizing transformation of metabolomic NMR data. OMICS: A Journal of Integrative Biology 8: 118–130.
  34. 34. Santos EM, Ball JS, Williams TD, Wu H, Ortega F, et al. (2010) Identifying health impacts of exposure to copper using transcriptomics and metabolomics in a fish model. Environmental Science & Technology 44: 820–826.
  35. 35. Dang W, Sun L (2011) Determination of internal controls for quantitative real time RT-PCR analysis of the effect of Edwardsiella tarda infection on gene expression in turbot (Scophthalmus maximus). Fish Shellfish Immunology 30, 720–728.
  36. 36. Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25: 402–408.
  37. 37. Peel GJ, Mickelbart MV, Rhodes D (2010) Choline metabolism in glycine betaine accumulating and non-accumulating near-isogenic lines of Zea mays and Sorghum bicolor. Phytochemistry 71: 404–414.
  38. 38. Moghaieb REA, Saneoka H, Fujita K (2004) Effect of salinity on osmotic adjustment, glycine betaine accumulation and the betaine aldehyde dehydrogenase gene expression in two halophytic plants, Salicornia europaea and Suaeda maritime. Plant Science 166: 1345–1349.
  39. 39. Graiia X, de Lecea L, El-Maghrabi MR, Ureiia JM, Caellas C, et al. (1992) Cloning and sequencing of a cDNA encoding 2,3-bisphosphoglycerate-independent phosphoglycerate mutase from maize. The Journal of Biological Chemistry 267: 12797–12803.
  40. 40. Sirover MA (1999) New insights into an old protein: the functional diversity of mammalian glyceraldehyde-3-phosphate dehydrogenase. Biochimica et Biophysica Acta 1432: 159–184.
  41. 41. Schultz CJ, Hsu M, Miesak B, Coruzzi GM (1998) Arabidopsis mutants define an in vivo role for isoenzymes of aspartate aminotransferase in plant nitrogen assimilation. Genetics 149: 491–499.
  42. 42. Wadsworth GJ, Marmaras SM, Matthews BF (1993) Isolation and characterization of a soybean cDNA clone encoding the plastid form of aspartate aminotransferase. Plant Molecular Biology 21: 993–1009.
  43. 43. Przybyla-Zawislak B, Dennis RA, Zakharkin SO, Mccammon MT (1998) Genes of succinyl-CoA ligase from Saccharomyces cerevisiae. European Journal of Biochemistry 258: 736–743.
  44. 44. Matsushika A, Goshima T, Fujii T, Inoue H, Sawayama S, et al. (2012) Characterization of non-oxidative transaldolase and transketolase enzymes in the pentose phosphate pathway with regard to xylose utilization by recombinant Saccharomyces cerevisiae. Enzyme and microbial technology 51: 16–25.
  45. 45. Li M-J, Li A-Q, Xia H, Zhao C-Z, Li C-S, et al. (2009) Cloning and sequence analysis of putative type II fatty acid synthase genes from Arachis hypogaea L. Journal of Biosciences. 34: 227–238.
  46. 46. Bunney TD, van Walraven HS, de Boer AH (2001) 14–3-3 protein is a regulator of the mitochondrial and chloroplast ATP synthase. Proceedings of the National Academy of Sciences 98: 4249–4254.
  47. 47. van Hemert MJ, Steensma HY, van Heusden GPH (2001) 14–3-3 proteins: key regulators of cell division, signalling and apoptosis. Bioessays 23: 936–946.
  48. 48. Lalle M, Camerini S, Cecchetti S, Sayadi A, Crescenzi M, et al. (2012) Interaction network of the 14–3-3 protein in the ancient protozoan parasite Giardia duodenalis. Journal of Proteome Research 11: 2666–2683.
  49. 49. Nagaoka S, Takano T (2003) Salt tolerance-related protein STO binds to a Myb transcription factor homologue and confers salt tolerance in Arabidopsis. Journal of Experimental Botany 54: 2231–2237.
  50. 50. Wang W, Vinocur B, Shoseyov O, Altman A (2004) Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends in plant science 9: 244–252.
  51. 51. Elbaz A, Wei YY, Meng Q, Zheng Q, Yang ZM (2010) Mercury-induced oxidative stress and impact on antioxidant enzymes in Chlamydomonas reinhardtii. Ecotoxicology 19: 1285–1293.
  52. 52. Pang Q, Chen S, Dai S, Chen Y, Wang Y, et al. (2010) Comparative proteomics of salt tolerance in Arabidopsis thaliana and Thellungiella halophile. Journal of Proteome Research 9: 2584–2599.
  53. 53. Ghosh S, Bagchi S, Majumder AL (2001) Chloroplast fructose-1,6-bisphosphatase from Oryza differs in salt tolerance property from the Porteresia enzyme and is protected by osmolytes. Plant Science 160: 1171–1181.
  54. 54. Yadav SK, Singla-Pareek SL, Kumar M, Pareek A, Saxena M, et al. (2007) Characterization and functional validation of glyoxalase II from rice. Protein Expression and Purification 51: 126–132.
  55. 55. Hoshida H, Tanaka Y, Hibino T, Hayashi Y, Tanaka A, et al. (2000) Enhanced tolerance to salt stress in transgenic rice that overexpresses chloroplast glutamine synthetase. Plant Molecular Biology 43: 103–111.
  56. 56. Tilton GB, Wedemeyer WJ, Browse J, Ohlrogge J (2006) Plant coenzyme A biosynthesis: characterization of two pantothenate kinases from Arabidopsis. Plant Molecular Biology 61: 629–642.
  57. 57. Raman SB, Rathinasabapathi B (2004) Pantothenate synthesis in plants. Plant Science 167: 961–968.
  58. 58. Chen J-G, Shimomura S, Sitbon F, Sandberg G, Jones AM (2001) The role of auxin-binding protein 1 in the expansion of tobacco leaf cells. The Plant Journal 28: 607–617.
  59. 59. Sauer M, Jakob A, Nordheim A, Hochholdinger F (2006) Proteomic analysis of shoot-borne root initiation inmaize (Zea mays L.). Proteomics 6: 2530–2541.
  60. 60. Romano PG, Horton P, Gray JE (2004) The Arabidopsis cyclophilin gene family. Plant Physiology 134: 1268–1282.
  61. 61. Wood ZA, Poole LB, Karplus PA (2003) Peroxiredoxin evolution and the regulation of hydrogen peroxide signaling. Science 300: 650–653.
  62. 62. Horling F, Lamkemeyer P, Konig J, Finkemeier I, Kandlbinder A, et al. (2003) Divergent light-, ascorbate-, and oxidative stress-dependent regulation of expression of the peroxiredoxin gene family in Arabidopsis. Plant Physiology 131: 317–325.
  63. 63. Ado K, Takeda N, Kikuchi M, Taniguchi Y (2006) The pressure effect on the structure and functions of protein disulfide isomerase. Biochimica et Biophysica Acta 1764: 586–592.
  64. 64. Huang D-J, Chen H-J, Lin Y-H (2005) Isolation and expression of protein disulfide isomerase cDNA from sweet potato (Ipomoea batatas [L.] Lam ‘Tainong 57′) storage roots. Plant Science 169: 776–784.
  65. 65. McDonagh B, David S (2007) Effect of oxidative stress on protein thiols in the blue mussel Mytilus edulis: Proteomic identification of target proteins. Proteomics 7: 3395–3403.
  66. 66. Marquez V, Frohlich T, Armache JP, Sohmen D, Donhofer A, et al. (2011) Proteomic characterization of archaeal ribosomes reveals the presence of novel archaeal-specific ribosomal proteins. Journal of Molecular Biology 405: 1215–1232.
  67. 67. Wahl MC, Bourenkov GP, Bartunik HD, Huber R (2000) Flexibility, conformational diversity and two dimerization modes in complexes of ribosomal protein L12. The EMBO Journal 19: 174–186.
  68. 68. Qin A, Shi Q, Yu X (2011) Ascorbic acid contents in transgenic potato plants overexpressing two dehydroascorbate reductase genes. Molecular Biology Reports 38: 1557–1566.
  69. 69. Urano J, Nakagawa T, Maki Y, Masumura T, Tanaka K, et al. (2000) Molecular cloning and characterization of a rice dehydroascorbate reductase. FEBS Letters 466: 107–111.
  70. 70. Ikegami A, Yoshimura N, Motohashi K, Takahashi S, Romano PG, et al. (2007) The CHLI1 subunit of Arabidopsis thaliana magnesium chelatase is a target protein of the chloroplast thioredoxin. The Journal of Biological Chemistry 282: 19282–19291.
  71. 71. Baginsky S, Gruissem W (2009) The chloroplast kinase network: new insights from large-scale phosphoproteome profiling. Molecular Plant 2: 1141–1153.
  72. 72. Fukayama H, Abe R, Uchida N (2010) SDS-dependent proteases induced by ABA and its relation to Rubisco and Rubisco activase contents in rice leaves. Plant Physiology and Biochemistry 48: 808–812.
  73. 73. Brooks A, Farquhar GD (1985) Effect of temperature on the CO2/O2 specificity of ribulose-1,5-bisphosphate carboxylase/oxygenase and the rate of respiration in the light. Planta 165: 397–406.
  74. 74. MacDonald RC, Kimmerer TW (1993) Metabolism of transpired ethanol by eastern cottonwood (Populus deltoides Bartr.) Plant Physiology. 102: 173–179.
  75. 75. Greenway H, Osmond CB (1972) Salt responses of enzymes from species differing in salt tolerance. Plant Physiology 49: 256–259.