Original Research ARTICLE
Effects of High Toxic Boron Concentration on Protein Profiles in Roots of Two Citrus Species Differing in Boron-Tolerance Revealed by a 2-DE Based MS Approach
- 1Institute of Plant Nutritional Physiology and Molecular Biology, College of Resources and Environment, Fujian Agriculture and Forestry University, Fuzhou, China
- 2College of Horticulture, Fujian Agriculture and Forestry University, Fuzhou, China
- 3Agriculture, Forestry and Water Conservancy Bureau of Xinzhou District, Shangrao, China
- 4Fujian Provincial Key Laboratory of Soil Environmental Health and Regulation, College of Resources and Environment, Fujian Agriculture and Forestry University, Fuzhou, China
- 5The Higher Educational Key Laboratory of Fujian Province for Soil Ecosystem Health and Regulation, Fujian Agriculture and Forestry University, Fuzhou, China
Citrus are sensitive to boron (B)-toxicity. In China, B-toxicity occurs in some citrus orchards. So far, limited data are available on B-toxicity-responsive proteins in higher plants. Thirteen-week-old seedlings of “Sour pummelo” (Citrus grandis) and “Xuegan” (Citrus sinensis) was fertilized every other day until dripping with nutrient solution containing 10 μM (control) or 400 μM (B-toxicity) H3BO3 for 15 weeks. The typical B-toxic symptom only occurred in 400 μM B-treated C. grandis leaves, and that B-toxicity decreased root dry weight more in C. grandis seedlings than in C. sinensis ones, demonstrating that C. sinensis was more tolerant to B-toxicity than C. grandis. Using a 2-dimensional electrophoresis (2-DE) based MS approach, we identified 27 up- and four down-accumulated, and 28 up- and 13 down-accumulated proteins in B-toxic C. sinensis and C. grandis roots, respectively. Most of these proteins were isolated only from B-toxic C. sinensis or C. grandis roots, only nine B-toxicity-responsive proteins were shared by the two citrus species. Great differences existed in B-toxicity-induced alterations of protein profiles between C. sinensis and C. grandis roots. More proteins related to detoxification were up-accumulated in B-toxic C. grandis roots than in B-toxic C. sinensis roots to meet the increased requirement for the detoxification of the more reactive oxygen species and other toxic compounds such as aldehydes in the former. For the first time, we demonstrated that the active methyl cycle was induced and repressed in B-toxic C. sinensis and C. grandis roots, respectively, and that C. sinensis roots had a better capacity to keep cell wall and cytoskeleton integrity than C. grandis roots in response to B-toxicity, which might be responsible for the higher B-tolerance of C. sinensis. In addition, proteins involved in nucleic acid metabolism, biological regulation and signal transduction might play a role in the higher B-tolerance of C. sinensis.
Boron (B) is an essential micronutrient for higher plants (Warington, 1923), where its most important role is associated with cell wall formation, functioning, and strength (Blevins and Lukaszewski, 1998). However, B will become toxic to crops when present in excess (Ben-Gal and Shani, 2003; Chen et al., 2012). B-toxicity is common in areas with high B concentration in underground water mainly resulting from the over-application of B fertilizer (Smith et al., 2013). In China, B-toxicity occurs in some citrus orchards. Up to 74.8 and 22.9% of pummelo (Citrus grandis) orchards in Pinghe, Zhangzhou, China, are excess in leaf B and soil water-soluble B, respectively (Li et al., 2015).
Plants have developed various mechanisms to cope with B-toxicity. Usually, antioxidant system will be activated to defense oxidative damage caused by B-toxicity (Cervilla et al., 2007; Ardic et al., 2009). Antioxidant compounds such as ascorbate and reduced glutathione (GSH) and antioxidant enzymes such as ascorbate peroxidase (APX), superoxide dismutase (SOD), catalase (CAT), and glutathione reductase (GR) are involved in the scavenging of reactive oxygen species (ROS) (Han et al., 2009; Erdal et al., 2014). B-tolerant plant leaves are characterized by a lower B concentration relative to B-sensitive ones, possibly due to a decreased uptake of B into both roots and shoots (leaves) (Camacho-Cristóbal et al., 2008). Sheng et al. (2010) showed that B-tolerant Newhall navel orange trees grafted on Carrizo citrange accumulated more B in roots and leaves than B-sensitive Skagg's Bonanza naval orange trees grafted on Carrizo citrange when exposed to B-toxicity, implying that the former must possess inner mechanisms to tolerate high level of B. Huang et al. (2014) reported that under B-toxicity, total B level was similar between B-tolerant Citrus sinensis and B-sensitive C. grandis roots (leaves), while C. sinensis leaves had lower free B and higher bound B than C. grandis leaves, which might contribute to the higher B-tolerance of C. sinensis. Our recent work with C. sinensis and C. grandis demonstrated that miR397a played a key role in citrus B-tolerance by targeting two laccase genes involved in secondary cell-wall biosynthesis (Huang et al., 2016). Similar result has been obtained on Poncirus trifoliata (Jin et al., 2016). To conclude, the mechanisms for plant B-tolerance are not fully understood yet.
A comprehensive investigation of B-toxicity-responsive proteins will be useful for us to unveil the inner mechanisms of B-tolerance in specific plant species. So far, knowledge on B-toxicity-induced alterations of protein profiles in higher plants is limited. Demiray et al. (2011) used a 2-dimensional electrophoresis (2-DE) based MS approach to identify six B-toxicity-responsive proteins from carrot roots. Atik et al. (2011) used a 2-DE technique to investigate the effects of B-toxicity on protein profiles in barley leaves, suggesting that a B-toxicity-responsive vacuolar H+-ATPase (V-ATPase) subunit E was involved in barley B-tolerance.
In higher plants, citrus are sensitive to B-toxicity (Eaton, 1935; Papadakis et al., 2004). Since B is phloem immobile in citrus (Konsaeng et al., 2005), the typical B-toxic symptom (chlorotic and/or necrotic patches) first develops in the older leaves and extends progressively from the old leaves to the young leaves (Han et al., 2009; Guo et al., 2014; Sang et al., 2015). It was indicated that great differences existed in B-tolerance among citrus species and/or genotypes (Chen et al., 2012). For example, when C. sinensis and C. grandis seedlings were submitted to 400 μM B for 15 weeks, the typical B-toxic symptom only occurred in the latter (Guo et al., 2014; Sang et al., 2015). We previously investigated the differences in B-toxicity-induced alterations of gene expression profiles in roots and leaves and of protein profiles in leaves between B-tolerant C. sinensis and B-sensitive C. grandis and revealed some adaptive responses of citrus to B-toxicity (Guo et al., 2014, 2016; Sang et al., 2015). Thus, B-toxicity-responsive proteins in roots should be different between C. sinensis and C. grandis. In this study, we used a 2-DE based MS approach to investigate comparatively B-toxicity-induced alterations of protein profiles in B-tolerant C. sinensis and B-sensitive C. grandis seedlings roots and corroborated the above hypothesis. For the first time, we demonstrated that the active methyl cycle was upregulated and downregulated in B-toxic C. sinensis and C. grandis roots, respectively, and that C. sinensis roots had a better capacity to keep cell wall and cytoskeleton integrity than C. grandis roots when exposed to B-toxicity, which might be involved in the higher B-tolerance of C. sinensis.
Materials and Methods
Plant Materials and Culture Conditions
This study was conducted at Fujian Agriculture and Forestry University, Fuzhou, China. Seeds of “Sour pummelo” (C. grandis) and “Xuegan” (C. sinensis) were germinated in clean river sand in plastic trays. Five weeks after germination, uniform seedlings with a single stem were transplanted to 6 L pots (two seedlings per pot) filled with clean river sand. Seedlings were grown in a greenhouse under natural photoperiod. Eight weeks after transplanting, each pot was fertilized every other day until dripping with nutrient solution (ca. 500 mL) containing 10 μM (control) or 400 μM (B-toxicity) H3BO3 for 15 weeks as described previously by Guo et al. (2014) and Sang et al. (2015). Thereafter, fully expanded (ca. 7-week-old) leaves were used for all the measurements. Leaf discs (0.2826 cm2 in size) were punched from each leaf using a hole puncher of 0.6 cm in diameter at noon at full sun and immediately frozen in liquid nitrogen. Approximately 5-mm-long white root apices were immediately frozen in liquid nitrogen after they were collected from the same seedlings used for sampling leaves. Both root and leaf samples were stored at −80°C until RNA and protein extraction, and the assay of malondialdehyde (MDA) concentration, H2O2 production and enzyme activities. The remaining seedlings that were not sampled were used to measure root dry weight (DW) and B concentration in fibrous roots, root apices and leaves.
Measurements of Root DW, and B and MDA Concentrations and H2O2 Production in Roots and Leaves
Roots of ten seedlings per treatment from 10 pots were harvested from the remaining seedlings and their DW was measured after being dried at 70°C for 48 h.
Fibrous roots, root apices and ca. 7-week-old fully expanded leaves (midribs and petioles removed) collected from the remaining seedlings were dried at 70°C, then ground to pass a 40-mesh sieve. Root and leaf B concentration was assayed by ICP emission spectrometry after microwave digestion with HNO3 (Wang et al., 2006). There were four replicates per treatment.
Root and leaf MDA was extracted and assayed according to Hodges et al. (1999). There were four replicates per treatment.
Root and leaf H2O2 production was determined according to Chen et al. (2005). About 40 mg frozen roots or 15 frozen leaf discs were incubated in 2 mL reaction mixture containing 50 mM of phosphate buffer (pH 7.0), 0.05% (w/v) of guaiacol and 5 U of horseradish peroxidase (Product No. 77332, lyophilized, powder, beige, ~150 U mg−1, Sigma-Aldrich, Shanghai, China) for 2 h at room temperature in the dark. Then, absorbance was assayed at 470 nm. There were four replicates per treatment.
Root Protein Extraction, 2-DE and Image Analysis
Approximately 1 g frozen roots collected equally from five seedlings (one seedling per pot) were mixed as one biological replicate. There were three biological replicates for each treatment (total of 15 seedlings from 15 pots). Proteins were independently extracted thrice from B-toxic and control samples according to You et al. (2014) using a phenol extraction procedure in order to ensure result reproducibility. Sample protein concentration was assayed according to Bradford (1976). 2-DE and image analysis were made according to Sang et al. (2015) and You et al. (2014). Gel images were obtained using Epson Scanner (Seiko Epson Corporation, Japan) at 300 dpi resolution. Image analysis was performed with PDQuest version 8.0.1 (Bio-Rad, Hercules, CA, USA). The software was used to perform background subtraction, Gaussian fitting, gel alignment, spot detection, matching and normalization. The parameters used to spot detection were as follow: sensitivity 6.05, size scale 3, min peak 600, and local regression model was selected to conduct spot normalization. The spot intensity was expressed as relatively abundant intensity that normalized by total intensities of all spots in one gel. After manual processing, the candidate spots in all triplicate gels were submitted to ANOVA analysis. A protein spot was considered differentially abundant when it had both a P < 0.05 and a fold change of > 1.5.
Protein Identification by MALDI-TOF/TOF-MS and Bioinformatic Analysis
MALDI-TOF/TOF-MS-based protein identification was performed on an AB SCIEX 5800 TOF/TOF (AB SCIEX, Shanghai, China) according to You et al. (2014) and Peng et al. (2015). Briefly, spots were excised from the colloidal Coomassie Brilliant Blue stained gels and plated into a 96-well microtiter plate. Excised spots were first destained twice with 60 μL of 50 mM NH4HCO3 and 50% (v/v) acetonitrile, and then dried twice with 60 μL of acetonitrile. Afterwards, the dried pieces of gels were incubated in ice-cold digestion solution [trypsin (sequencing-grade modified trypsin V5113, Promega, Madison, WI, USA) 12.5 ng/μL and 20 mM NH4HCO3] for 20 min, and then transferred into a 37°C incubator for digestion overnight. Peptides in the supernatant were collected after extraction twice with 60 μL extract solution [5% (v/v) formic acid in 50% (v/v) acetonitrile]. The resulting peptide solution was dried under the protection of N2. Before MS/MS analysis, the pellet was redissolved in 0.8 μL matrix solution [5 mg/mL α-cyano-4-hydroxy-cinnamic acid diluted in 0.1% trifluoroacetic acid (TFA), 50% (v/v) acetonitrile]. Then the mixture was spotted onto a MALDI target plate (AB SCIEX, Shanghai, China). MS analysis of peptide was performed on an AB SCIEX 5800 TOF/TOF. The UV laser was operated at a 400 Hz repetition rate with a wavelength of 355 nm. The accelerated voltage was operated at 20 kV, and mass resolution was maximized at 1,600 Da. Myoglobin digested with trypsin was used to calibrate the mass instrument with internal calibration mode. All acquired spectra of samples were processed using TOF/TOF Explorer™ Software (AB SCIEX, Shanghai, China) in a default mode. The data were searched by GPS Explorer (Version 3.6) with the search engine MASCOT (Version 2.3, Matrix Science Inc., Boston, MA). The search parameters were as follows: viridiplantae database (1,850,050 sequences; 6,42,453,415 residues), trypsin digest with one missing cleavage, MS tolerance was set at 100 ppm, MS/MS tolerance was set at 0.6 Da. At least two peptides were required to match for each protein. Protein identifications were accepted if MASCOT score was not less than 75, and the number of matched peptides was not less than five or the sequence coverage was not less than 20% (Lee et al., 2010; You et al., 2014). Searches were also performed against the C. sinensis databases (https://phytozome.jgi.doe.gov/pz/portal.html#!info?alias$=$Org_Csinensis).
Bioinformatics analysis of proteins was performed according to Yang et al. (2013).
Principal Components Analysis (PCA) of Differentially Abundant Proteins (DAPs)
The ratios of all the DAPs from B-toxic C. sinensis and C. grandis roots were normalized and transformed for the PCA using Princomp function in R circumstance. The first two components were selected and used to visualize two loadings against each other to investigate the relationships between the variables (Mardia et al., 1979). The PCA loading plots were carried out in triplicate.
Approximately 300 mg frozen roots collected equally from five seedlings (one seedling per pot) were pooled as one biological replicate. qRT-PCR analysis was run in three biological and two technical replicates for each treatment (total of 15 seedlings from 15 pots) according to Zhou et al. (2013). In this study, we randomly selected ten DAPs from each citrus species for qRT-PCR analysis. A total of 20 DAPs were selected from B-toxic C. sinensis and C. grandis roots. Specific primers were designed from the corresponding sequences of these selected DAPs in citrus genome (https://phytozome.jgi.doe.gov/pz/portal.html#!info?alias$=$Org_Csinensis) using Primer Primier Version 5.0 (PREMIER Biosoft International, CA, USA). The sequences of the F and R primers used were listed in Table S1. For the normalization of gene expression and reliability of quantitative analysis, two citrus genes [C. sinensis NADP-dependent glyceraldehyde-3-phosphate dehydrogenase (GAPDH; gi|985455672) and C. sinensis DNA-directed RNA polymerase II subunit 4 (RPII; gi|985473508] were selected as internal standards and the roots from control seedlings were used as reference sample, which was set to 1.
Assay of S-Adenosylmethionine (SAM) Synthetase (SAMS) and Adenosine Kinase (ADK)
Both ADK and SAMS were extracted according to Shen et al. (2002) by homogenizing ca. 100 mg of frozen roots in 1 mL extraction buffer including 100 mM of Tris (pH 7.5), 2 mM of ethylenediaminetetraacetic acid (EDTA), 20% (w/v) of glycerol, 20 mM of β-mercaptoethanol, 1 mM of dithiothreitol (DTT) at 4°C. After centrifugation at 10,000 g for 10 min, the supernatant was used immediately for enzyme assay. There were four replicates per treatment.
Total SAMS activity was assayed as described by Kim et al. (1992) and Shen et al. (2002). Briefly, 135 μL of an enzyme extract was incubated in 0.45 mL of a reaction mixture containing 100 mM of Tris (pH 8.0), 30 mM of MgSO4, 10 mM of KCl, 20 mM of ATP, and 5 mM of methionine. Blank contained all reagents except for methionine. Reaction was incubated for 1 h at 25°C and was terminated by adding 0.5 mL of 6% (w/v) sodium dodecyl sulfate (SDS), and the phosphate (Pi) released from the substrate was determined as described by Smith et al. (1984) by adding 0.6 mL of an assay mixture containing 6 parts of 3.6 mM ammonium molybdate in 0.5 M H2SO4, and 1 part of 10% (w/v) ascorbic acid. The sample was incubated at 37°C for 60 min and the absorbance measured at 820 nm.
Root ADK activity was assayed according to Chen and Eckert (1977) and Lindberg et al. (1967) in 1 mL of a reaction mixture containing 20 mM of Tris-maleate (pH 5.8), 0.7 mM of ATP, 0.25 mM of phosphoenolpyruvate (PEP), 0.2 mM of NADH, 0.5 mM of MgCl2, 50 mM of KCl, 0.05 mM of adenosine, 5 U of pyruvate kinase, 5 U of lactate dehydrogenase, and 0.1 mL of enzyme extract. The reaction mixture was always preincubated for 10 min at room temperature (25°C) with all of the regents before starting the reaction by the addition of adenosine.
Experimental Design and Statistical Analysis
There were 20 pots (40 seedlings) per treatment in a completely randomized design. Experiments were performed with 3–10 replicates. Significant differences among four treatments were analyzed by two (species) × two (B levels) ANOVA and four means were separated by the Duncan's new multiple range test at P < 0.05. Significant tests between two means (B-toxicity and control) were performed by unpaired t-test at P < 0.05 level.
Results and Discussion
C. sinensis Was More Tolerant to B-Toxicity than C. grandis
In previous studies, we showed that a concentration of 400 μM B is suitable for the comparative investigation of B-tolerance between B-tolerant C. sinensis and B-sensitive C. grandis. The typical B-toxic symptoms only occurred in C. grandis leaves (Guo et al., 2014; Huang et al., 2014; Sang et al., 2015). We therefore decided to use this B treatment in the present work to reveal specific root proteome signatures in tolerant and sensitive citrus species. As shown in Figures 1A–D, 400 μM B-treatment greatly decreased root DW, increased the concentration of B in leaves, fibrous roots and root apices, and the concentration of B in 400 μM B-treated leaves was far more than the sufficiency range of 30–100 mg kg−1 DW for citrus (Chapman, 1968). Thus, seedlings that received 10 and 400 μM B are considered as B-toxic and B-sufficient (control), respectively.
Figure 1. Effects of B-toxicity on root DW (A), B concentration in leaves (B), fibrous roots (C) and root apices (D), H2O2 production in roots (E) and leaves (F), and MDA concentration in roots (G) and leaves (H). Bars represent means ± SD (n = 10 for root DW and 4 for other parameters). Differences among four treatments were analyzed by two (species) × two (B levels) ANOVA. Different letters above the bars indicate a significant difference at P < 0.05.
Our results showed that the B-toxicity-induced decrease in root DW (Figure 1A) and increase in H2O2 production in roots and leaves (Figures 1E,F) were greater in C. sinensis seedlings than in C. grandis ones, and that B-toxicity increased the concentration of MDA only in C. grandis leaves (Figure 1H). In addition, the typical visible B-toxic symptom only occurred in B-toxic C. grandis leaves, but was not found in B-toxic C. sinensis leaves except for very few seedlings (Figure S1). Previous studies showed that B-toxicity only decreased the concentrations of phosphorus (P) and total soluble proteins in C. grandis roots (Guo et al., 2016). Based on these results, we concluded that C. sinensis had higher B-tolerance than C. grandis.
Protein Yield and DAPs in B-Toxic Roots
Protein yield did not differ among four treatment combinations (Table 1). After Coomassie Brilliant Blue G-250 staining, more than 800 clear and reproducible protein spots were discovered on each gel. The number of protein spots per gel were similar among the four treatment combinations (Table 1 and Figure 2; Figure S2), as obtained on C. sinensis and C. grandis leaves (Sang et al., 2015).
Table 1. Protein yield, number of spots, number of variable spots and number of identified differentially abundant protein spots in Citrus sinensis and Citrus grandis roots.
Figure 2. Representative 2-DE images of proteins extracted from control (A,C) and B-toxic (B,D) roots. (A) Control roots of Citrus sinensis, (B) B-toxic root of C. sinensis, (C) Control roots of Citrus grandis, (D) B-toxic roots of C. grandis.
We detected 43 up- and five down-accumulated, and 35 up- and 20 down-accumulated protein spots from B-toxic C. sinensis and C. grandis roots, respectively. Twenty-seven up- and four down-accumulated, and 28 up- and 13 down-accumulated protein spots were identified from B-toxic C. sinensis and C. grandis roots, respectively after these differentially accumulated protein spots being submitted to the MALDI-TOF/TOF-MS-based identification (Table 1, Figures 2, 3 and Tables S2–S5). These DAPs were mainly involved in protein and amino acid metabolism, stress response, cell wall and cytoskeleton metabolism, carbohydrate and energy metabolism, nucleic acid metabolism, cellular transport, and biological regulation and signal transduction (Tables 2, 3 and Figures 4A,B). Most of B-toxicity-responsive proteins were isolated from B-toxic C. sinensis or C. grandis roots, only nine protein species with the same accession No. were shared by the both. Among the nine overlapping proteins, only five proteins displayed similar change trends in B-toxic C. sinensis and C. grandis roots (Tables 2, 3 and Figure 4C). These results demonstrated that B-toxicity-responsive proteins greatly differed between C. sinensis and C. grandis roots, as obtained on B-toxic C. sinensis and C. grandis leaves (Sang et al., 2015).
Figure 3. Close-up views of the differentially abundant protein spots in control (CK) and B-toxic (BT) roots.
Table 2. Differentially abundant proteins and their identification by MALDI-TOF/TOF-MS in B-toxic Citrus sinensis roots.
Table 3. Differentially abundant proteins and their identification by MALDI-TOF/TOF-MS in B-toxic Citrus grandis roots.
Figure 4. Classification of B-toxicity-responsive protein spots in C. grandis (A), C. sinensis (B) roots, and venn diagram analysis of B-toxicity-responsive protein spots in citrus roots (C).
Principal Component Analysis Loading Plots of DAPs
As shown in Figure 5, 31 and 41 B-toxicity-responsive proteins identified in C. sinensis and C. grandis roots were submitted to PCA procedure. The first two components accounted for 94.6% (70.8% for PC1 and 23.8% for PC2) and 91.8% (69.5% for PC1 and 22.3% for PC2) of total variation in C. sinensis and C. grandis roots, respectively. The DAPs involved in protein and amino acid metabolism and cell wall and cytoskeleton were highly clustered in C. sinensis roots. In contrast, no obvious clustered proteins were observed in C. grandis roots.
Figure 5. Principal component analysis (PCA) loading plots of differentially abundant proteins in B-toxic C. sinensis (A) and C. grandis (B) roots.
qRT-PCR Analysis of Genes for DAPs
The mRNA levels of genes encoding 20 B-toxicity-responsive proteins from C. sinensis (S27, 47, 43, 12, 41, 39, 13, 31, 20, and 36) and C. grandis (G32, 26, 22, 16, 43, 46, 21, 28, 20, and 11) roots (Figures 6A–D) were assayed in order to examine the relationship between the abundances of proteins and the expression levels of genes. The expression levels of all genes except for S39 and G26 matched well with our 2-DE data (Tables 2, 3) regardless of which internal standard was used to calculate the relative expression levels, suggesting that most of B-toxicity-responsive proteins were regulated at the transcriptional level. This is also supported by our analysis that the qRT-PCR data and the 2-DE results were significantly and positively correlated (Figures 6E,F).
Figure 6. Relative expression levels of genes encoding 20 B-toxicity-responsive proteins from C. sinensis (A,B) and C. grandis (C,D) roots using GAPDH (A,C) and RPII (B,D) as internal standards, and the correlation analysis of qRT-PCR results and 2-DE data (E,F). For (A–D), bars represent means ± SD (n = 3); Significant tests between two means were performed by unpaired t-test; Different letters above the bars indicate a significant difference at P < 0.05. For (E,F), 2-DA data from Tables 2, 3.
Proteins Related to Carbohydrate and Energy Metabolism
ADK plays key roles in the maintenance of purine nucleotide pools and in the active methyl cycle (Moffatt et al., 2002; Kettles et al., 2014; Chen et al., 2015). Schoor et al. (2011) observed that silencing of ADK in Arabidopsis caused impaired root growth, small, crinkled rosette leaves, and decreased apical dominance accompanied by an increased concentration of active cytokinin (CK) ribosides, concluding that ADK was responsible for CK homeostasis in vivo. We found that the abundances of ADK2 (S27) and ADK isoform 1T-like protein (S28) were increased in B-toxic C. sinensis roots, while only one down-accumulated ADK isoform 1T-like protein (G24) was identified in B-toxic C. grandis roots (Tables 2, 3). Similarly, B-toxicity increased the abundances of S-adenosylmethionine synthetase 2 [also known as S-adenosyl-L-methionine synthetase 2 (SAMS2); S41] involved in the formation of SAM from methionine and ATP, and 5-methyltetrahydropteroyltriglutamate-homocysteine methyltransferase (also known as methionine synthase; S39) involved in the biosynthesis of methionine in C. sinensis roots and SAMS (G44) in C. grandis roots, and decreased the abundances of two SAMS1 family protein (G33 and 43) in C. grandis roots (Tables 2, 3). In addition, the activities of both SAMS and ADK were increased in B-toxic C. sinensis roots, but decreased in B-toxic C. grandis roots (Figure 7). Therefore, the active methyl cycle might be upregulated in B-toxic C. sinensis roots, but downregulated in B-toxic C. grandis roots. This agrees with the report that the active methyl cycle was induced by drought in drought-resistant rice leaves, but inhibited in drought sensitive rice leaves, and that the cycle played a role in rice drought resistance (Zhang et al., 2012). It is known that SAM not only plays a role in the active methyl cycle but also serves as an intermediate in the biosynthesis of polyamines (PAs) and ethylene (Ravanel et al., 1998). Hassan et al. (2010) reported that the expression of genes encoding SAM decarboxylase [SAMDC, a key enzyme involved in the biosynthesis of PAs (spermidine and spermine)], methinnine synthase 1 and SAMS2 was upregulated in B-tolerant Sahara barley roots, and that an antioxidant mechanisms involving PAs and water-water cycle in Sahara barley might play a role in tolerating high level of soil B. Hassan et al. (2010) also suggested that increased activity of SAMDC on SAM might inhibit ethylene production, hence reducing leaf senescence in Sahara barley. Evidence shows that transgenic plants with elevated levels of PAs have enhanced tolerance to different abiotic stresses (Alcázar et al., 2006). Recently, Tanou et al. (2014) observed that exogenous PAs partially alleviated the NaCl-induced phenotypic and physiological impairments in citrus plants, and systematically upregulated the expression of genes involved in PA biosynthesis (arginine decarboxylase, SAMDC, spermidine synthase, and spermine synthase) and catabolism (diamine oxidase and polyamine oxidase). Also, PAs reprogrammed the oxidative status in salt-stressed citrus plants. Based on these results, we concluded that the B-toxicity-induced upregulation of the active methyl cycle might play a role in the B-tolerance of C. sinensis via enhancing the biosynthesis of PAs.
Figure 7. Effects of B-toxicity on the activities of SAMS (A) and ADK (B) in C. sinensis and C. grandis roots.Bars represent means ± SD (n = 4). Differences among four treatments were analyzed by two (species) × two (B levels) ANOVA. Different letters above the bars indicate a significant difference at P < 0.05.
We found that the abundances of all the four B-toxicity-responsive proteins involved in glycolysis (G22, 42, and 51) and tricarboxylic acid (TCA) cycle (G50) were increased in C. grandis roots and that ATP synthase subunit α (G46) was down-accumulated in B-toxic C. grandis roots (Table 3). Thus, ATP synthase-mediated ATP biosynthesis might be decreased in these roots. This might contribute to the maintenance of ATP balance, when the production of ATP was increased due to upregulated glycolysis and TCA cycle and the consumption of ATP was decreased due to decreased activities of ADK and SAMS. However, the abundance of malate dehydrogenase (MDH, S47) involved in TCA cycle was decreased in B-toxic C. sinensis roots (Table 2).
Stress Response-Related Proteins
Because the production of ROS (H2O2) was increased in B-toxic C. grandis and C. sinensis roots, especially in the former (Figure 1E), antioxidant enzymes might be induced in these roots. As expected, the abundance of Cu/Zn-SOD (G32) was increased in B-toxic C. grandis roots (Table 3). Besides antioxidant enzymes, the abundance of lactoylglutathione lyase (LGL; G26) was augmented in B-toxic C. grandis roots. In addition to the detoxification of methylglyoxal, a cytotoxic compound formed spontaneously from the glycolysis and photosynthesis intermediates glyceraldehyde-3-phosphate and dihydroxyacetone phosphate, LGL also play a role in oxidative stress tolerance (Yadav et al., 2005). The increased abundance of LGL agrees with our data that the abundances of four protein species involved in glycolysis (G22, 42, and 51) and TCA cycle (G50) were elevated in B-toxic C. grandis roots. By contrast, we only obtained one up-accumulated late-embryogenesis abundant protein 2 (LEA-2; S17) from B-toxic C. sinensis roots (Tables 2, 3). To conclude, more proteins related to detoxification were up-accumulated in B-toxic C. grandis roots than in B-toxic C. sinensis roots, which agrees with the increased requirement for detoxification of the more ROS and other toxic compounds such as aldehydes in the former because the production of ROS was higher in B-toxic C. grandis than in B-toxic C. sinensis roots (Figure 1E). We found that the level of MDA did not differ between B-toxic roots and controls (Figure 1C), demonstrating that the upregulation of antioxidant system provided sufficient protection to B-toxic roots against oxidative damage.
Proteins Related to Cell Wall and Cytoskeleton
All of the identified DAPs in cytokeleton were up-accumulated in B-toxic C. sinensis roots, while we isolated two down-accumulated proteins in cytokeleton (actin 1, G34) and polysaccharide biosynthesis (α-1,4-glucan-protein synthase 1; G52), and one up-accumulated α-tubulin in cytokeleton (G25) from B-toxic C. grandis roots (Tables 2, 3). Thus, C. sinensis roots might have a better capacity to keep cytoskeleton and cell wall integrity than C. grandis roots under B-toxicity, which might be responsible for the higher B-tolerance of the former. Similar results have been obtained on B-toxic C. sinensis and C. grandis leaves (Sang et al., 2015).
Proteins Related to Protein, Amino Acid, and Nucleic Acid Metabolisms
Proteasomes are responsible for the degradation of the inactive and futile proteins. Most of proteins degraded by proteasomes are first tagged by ubiquitin (Kurepa and Smalle, 2008). We obtained two up-accumulated proteasomes (G6 and 13) and two up-accumulated ubiquitin-conjugating enzymes (G48 and 54) from B-toxic C. grandis roots, but only one up-accumulated proteasome (S9) from B-toxic C. sinensis (Tables 2, 3), demonstrating that B-toxicity accelerated proteolysis, especially in the former. This agrees with our data that B-toxicity only decreased total soluble protein concentration in B-toxic C. grandis roots (Guo et al., 2016). B-toxicity-induced increase in protein degradation implies that misfolded and damaged proteins were increased in B-toxic C. sinensis and C. grandis roots, especially in the latter. In addition, we identified one up-accumulated α chain of nascent polypeptide associated complex (α-NAC, S8) from B-toxic C. sinensis roots. NAC, including α and β subunits, plays a role in protecting newly synthesized polypeptides on ribosome from proteolysis and in facilitating its folding (Karan and Subudhi, 2012; Kogan and Gvozdev, 2014). Thus, the up-accumulation of α-NAC in B-toxic C. sinensis might alleviate B-toxicity induced protein degradation and misfolding, hence preventing the reduction of proteins. All B-toxicity-responsive proteins (S23, 40, and 36, and G7, 40, 11, and 41) in protein biosynthesis were up-accumulated in C. sinensis and C. grandis roots except for eukaryotic initiation factor 5A (S34 and G30) (Tables 2, 3). Therefore, the lower level of total soluble proteins in B-toxic C. grandis roots might be mainly caused by increased proteolysis rather than by decreased biosynthesis.
As shown in Table 2 and Figure 7A, the abundances of the five DAPs (S43, 41, 39, 42, and 31) involved in amino acid metabolism and the activity of SAMS were increased in B-toxic C. sinensis roots, suggesting that the biosynthesis of some amino acids might be enhanced in these roots. SAM is an allosteric activator of threonine synthase (TS), which is involved in the biosynthesis of branched chain amino acids (BCAAs, valine, leucine and isoleucine; Curien et al., 1998; Ravanel et al., 1998). Thus, TS might be activated in B-toxic C. sinensis roots due to increased SAM biosynthesis resulting from enhanced SAMS activity. Zeh et al. (2001) found that antisense inhibition of TS led to increased level of methionine in transgenic potato plants because of the redirection of carbon flow from the threonine to the methionine branch. Ketol-acid reductoisomerase (KARI) is involved in BCAA biosynthesis. Kochevenko and Fernie (2011) reported that the concentrations of BCAAs in the leaves of transgenic tomato line 7, in which the KARI transcript level was remained at about 70% of the wildtype level, were not lower than in the wildtype leaves, whereas BCAA levels in the leaves of transgenic lines 3 and 6, in which the expression level of KARI was decreased to 27% of the wild-type level, were only 49–79% of the wildtype leaves. Thus, the levels of BCAAs might be enhanced in B-toxic C. sinensis roots due to the activation of TS resulting from increased SAMS activity (Figure 7A) and the increased abundances of KARI (Table 2). However, we obtained two down-accumulated SAMS1 (G33 and 43), one up-accumulated SAMS (G44), and two up-accumulated KARI spots (G27 and 28) from B-toxic C. grandis roots. In addition, SAMS activity was reduced in B-toxic C. grandis roots (Figure 7A). Based on these results, we concluded that BCAA biosynthesis might be disturbed in these roots.
We observed that the abundances of glycine-rich RNA-binding protein (GR-RBP) in C. sinensis (S21) and C. grandis (G21) roots and GR-RBP4 in C. sinensis roots (S13) increased when exposed to B-toxicity (Tables 2, 3), which agrees with the reports that the transcript levels of the genes encoding GR-RBPs were increased in higher plants following exposure to various abiotic stresses (Sachetto-Martins et al., 2000; Kwak et al., 2005). DEAD box RNA helicases, which may actively disrupt misfolded RNA structures by utilizing energy produced from ATP hydrolysis so that correct folding can occur, play key roles in plant response to various stresses (Li et al., 2008; Zhu et al., 2015). Thus, the down-accumulation of DEAD-box RNA helicase-like protein (G35) in B-toxic C. grandis roots might decrease C. grandis stress-tolerance, hence impairing its B-tolerance. However, the abundance of spliceosome RNA helicase BAT1 (G36) was increased in B-toxic C. grandis roots (Table 3).
Proteins Related to Cellular Transport
The abundances of four protein spots involved in cellular transport were increased in B-toxic C. sinensis (S20) and C. grandis (G18, 20, and 54) roots (Tables 2, 3), as reported on B-toxic leaves of barley (Atik et al., 2011), C. sinensis and C. grandis (Sang et al., 2015). However, the mRNA levels of all 13 B-toxicity-responsive genes involved in cellular transport were downregulated in C. grandis and C. sinensis roots except for one upregulated H+-ATPase 4 (Guo et al., 2016). The difference between protein abundances and gene expression levels implies that post-translational modifications (PTMs) might affect protein levels. Atik et al. (2011) observed that heterologous expression of a gene encoding V-ATPase subunit E, a protein induced by B-toxicity in barley leaves, conferred yeast B-tolerance. The up-accumulation of V-ATPase B subunit (S20, and G18 and 20) in the two citrus species might be an adaptive response to B-toxicity by providing energy for compartmentation of excess B in vacuoles (Alemzadeh et al., 2006; Wang et al., 2011). Klychnikov et al. (2007) showed that plant V-ATPase could interact with 14-3-3 proteins. The up-accumulation of V-ATPase B subunit in B-toxic C. sinensis roots agrees with our data that the abundance of 14-3-3 family protein (S14) was increased in these roots. However, the abundance of 14-3-3-like protein GF14 phi (G9 and 10) was decreased in B-toxic C. grandis roots, implying that other factors was involved in the regulation of V-ATPase B subunit.
Proteins Related to Biological Regulation and Signal Transduction
B-toxicity increased the abundance of nucleoside diphosphate kinase 1 (NDPK1; S44) and 14-3-3 family protein (S14) in C. sinensis roots (Table 2), as observed on B-toxic C. sinensis leaves (Sang et al., 2015). However, the abundances of 14-3-3-like protein GF14 phi (G9 and 10) were reduced in B-toxic C. grandis roots (Table 3). Fukamatsu et al. (2003) demonstrated that Arabidopsis NDPK1 played a role in ROS response by interacting with three CATs. Overexpression of NDPKs conferred enhanced tolerance to multiple abiotic stresses in potato, alfalfa and poplar (Fukamatsu et al., 2003; Tang et al., 2008; Wang et al., 2014). 14-3-3 proteins, the master regulators of many signal transduction cascades, have a key role in stress-tolerance in higher plants (Chen et al., 2006). Transgenic potato plants overexpressing 14-3-3 protein genes displayed delayed leaf senescence and enhanced antioxidant activity, while transgenic potato plants with antisense 14-3-3 protein genes displayed early leaf senescence (Wilczynski et al., 1998; Lukaszewicz et al., 2002). Thus, the B-toxicity-induced up-accumulation of NDPK1 and 14-3-3 might contribute to the higher B-tolerance of C. sinensis.
Comparison of B-Toxicity-Responsive Proteins between Roots and Leaves
More B-toxicity-responsive proteins were identified in C. grandis (41) roots than in C. sinensis (31) roots (Tables 2, 3), while Sang et al. (2015) identified 45 and 55 DAPs from B-toxic C. grandis and C. sinensis leaves, respectively. As shown in Tables 2, 3 and Figures 4A,B, we identified more up-accumulated proteins than down-accumulated proteins in B-toxic C. sinensis and C. grandis roots, especially in B-toxic C. sinensis roots, but the reverse was the case in B-toxic C. grandis leaves although the number of up-accumulated proteins (27) in B-toxic C. sinensis leaves was slightly higher than that of down-accumulated proteins (23) (Sang et al., 2015). Furthermore, the vast majority of B-toxicity-responsive proteins were identified only in C. sinensis and C. grandis roots or leaves, only three proteins with the same accession No. were shared by C. grandis roots and leaves (Table 4). In addition, many other differences existed in B-toxicity-responsive proteins between roots and leaves of the two citrus species. For examples, the carbohydrate and energy metabolism-related proteins was the largest category of the B-toxicity-responsive proteins in C. sinensis and C. grandis leaves (Sang et al., 2015). Similar result has been obtained on NaCl-stressed Citrus aurantium leaves (Tanou et al., 2009). However, the protein and amino acid metabolism-related proteins was the most abundant category of B-toxicity-responsive proteins in C. sinensis and C. grandis roots (Tables 2, 3). In the previous study, we isolated similar up- (11) and down-accumulated (9) carbohydrate and energy metabolism-related proteins in B-toxic C. sinensis leaves, and more down- (16) than up-accumulated (9) proteins in B-toxic C. grandis leaves (Sang et al., 2015). However, more up- than down-accumulated proteins involved in carbohydrate and energy metabolism were identified in B-toxic C. sinensis (two up- and one down-accumulated) and C. grandis (four up- and two down-accumulated) roots (Tables 2, 3). As shown in Table 3, all four DAPs related to glycolysis (G22, 42, and 51) and tricarboxylic acid (TCA) cycle (G50) were up-accumulated, and ATP synthase subunit α (G46) involved in ATP was decreased in B-toxic C. grandis roots. By contrast, we isolated four down- and three up-accumulated proteins in glycolysis and TCA cycle and one up-accumulated mitochondrial ATP synthase in B-toxic C. grandis leaves (Sang et al., 2015). In B-toxic C. sinensis leaves, we obtained three up-accumulated malate dehydrogenases (MDHs), while only one down-accumulated MDH (S47) was detected in B-toxic C. sinensis roots (Table 2). Thus, the adaptive responses of carbohydrate and energy metabolism-related proteins to B-toxicity differed between roots and leaves of the two citrus species.
B-toxicity increased the abundances of proteins involved in protein degradation, and decreased the abundances of proteins related to protein biosynthesis in C. grandis and C. sinensis leaves (Sang et al., 2015). By contrast, the abundances of the two kinds of proteins were enhanced in B-toxic C. grandis and C. sinensis roots (Tables 2, 3). Interestingly, total soluble protein level was reduced only in B-toxic C. grandis roots and leaves (Sang et al., 2015; Guo et al., 2016). Thus, the causes for the decrease of total soluble proteins in B-toxic C. grandis roots and leaves were different.
We isolated seven up-accumulated, and one down- and three up-accumulated proteins involved in antioxidation and detoxification from B-toxic C. grandis and C. sinensis leaves, respectively (Sang et al., 2015), and two up-accumulated (G32 and 26) and one up-accumulated (S17) from B-toxic C. grandis and C. sinensis roots, respectively (Tables 2, 3). However, MDA concentration was increased only in B-toxic C. grandis leaves (Figures 1G,H). This might be related to the findings that B mainly accumulated in B-toxic C. grandis and C. sinensis leaves (Figures 1G,H; Jiang et al., 2009; Guo et al., 2014), and that the increased requirement for the detoxification of ROS and other toxic compounds such as reactive aldehtdes was greater in B-toxic C. grandis leaves than in B-toxic C. sinensis leaves (Sang et al., 2015).
To conclude, B-toxicity-induced alterations of protein profiles greatly differed between roots and leaves of the two citrus species.
Using a 2-DE based MS approach, we comparatively investigated the effects of B-toxicity on DAPs in roots of two citrus species with different B-tolerance and obtained 27 up- and four down-accumulated, and 28 up- and 13 down-accumulated proteins from B-toxic C. sinensis and C. grandis) roots, respectively. Most of B-toxicity-responsive proteins only were isolated from C. sinensis or C. grandis roots, only nine proteins were shared by the both. Great differences existed in B-toxicity-induced alterations of protein profiles between C. sinensis or C. grandis roots. Based on our findings, a diagram for the responses of C. grandis and C. sinensis roots to B-toxicity was presented in Figure 8. The higher B-tolerance of C. sinensis might be associated with (a) the B-toxicity-induced upregulation of the active methyl cycle and (b) the better performance in maintaining cell wall and cytoskeleton integrity. In addition, proteins related to nucleic acid metabolism, biological regulation and signal transduction might play a role in the higher B-tolerance of C. sinensis. To conclude, our findings provided some novel cues on the molecular mechanisms of citrus B-toxicity and B-tolerance.
Figure 8. A diagram for the responses of C. grandis and C. sinensis roots to B-toxicity. CG, C. grandis; CGR, C. grandis roots; CS, C. sinensis; CSR, C. sinensis roots; DBRHL, DEAD-box RNA helicase-like protein; DRS, dihydrolipoyllysine-residue succinyltransferase component of 2-oxoglutarate dehydrogenase complex; EF, elongation factor; EIF, eukaryotic initiation factor; ETIF, eukaryotic translation initiation factor; GPS1, α-1,4-glucan-protein synthase 1; MDM, 5-methyltetrahy dropteroyltriglutamate-homocysteine methyltransferase; PGK, Phosphoglycerate kinase; 26SP, 26S proteasome subunit RPN12; TIF, translation initiation factor; TPI, triosephosphate isomerase; TPIL, triosphosphate isomerase-like protein; UBE, ubiquitin-conjugating enzyme; VHA-B, Vacuolar H+-ATPase B subunit.
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD004050.
WS and ZH contributed equally to this works. WS carried most of the experiment and analyzed the data; ZH drafted the manuscript; LY participated in the direction of this study; PG performed the qRT-PCR analysis; XY participated in the analysis of B; LC designed and directed the study and revised the manuscript. All authors have read and approved the final manuscript.
This work was financially supported by the National Natural Science Grant of China (No. 31301740), the earmarked fund for China Agriculture Research System (No. CARS27) and the Provincial Natural Science Grant of Fujian, China (No. 2014J05033).
Conflict of Interest Statement
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.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2017.00180/full#supplementary-material
Alcázar, R., Marco, F., Cuevas, J. C., Patron, M., Ferrando, A., Carrasco, P., et al. (2006). Involvement of polyamines in plant response to abiotic stress. Biotechnol. Lett. 28, 1867–1876. doi: 10.1007/s10529-006-9179-3
Alemzadeh, A., Fujie, M., Usami, S., Yoshizaki, T., Oyama, K., Kawabata, T., et al. (2006). ZMVHA-B1, the gene for subunit B of vacuolar H+-ATPase from the eelgrass Zostera marina L. is able to replace vma2 in a yeast null mutant. J. Biosci. Bioeng. 102, 390–395. doi: 10.1263/jbb.102.390
Ardic, M., Sekmen, A. H., Turkan, I., Tokur, S., and Ozdemir, F. (2009). The effects of boron toxicity on root antioxidant systems of two chickpea (Cicer arietinum L.) cultivars. Plant Soil 314, 99–108. doi: 10.1007/s11104-008-9709-y
Atik, A. E., Bozdağ, G. O., Akinci, E., Kaya, A., Koc, A., Yalcin, T., et al. (2011). Proteomic changes during boron tolerance in barley (Hordeum vulgare) and the role of vacuolar proton-translocating ATPase subunit E. Turk. J. Bot. 35, 379–388. doi: 10.3906/bot-1007-29
Bradford, M. M. (1976). A rapid and sensitive method for quantitation of microgram and quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72, 248–254. doi: 10.1016/0003-2697(76)90527-3
Cervilla, L. M., Blasco, B., Ríos, J. J., Romero, L., and Ruiz, J. M. (2007). Oxidative stress and antioxidants in tomato (Sola um lycopericum) plants subjected to boron toxicity. Ann. Bot. 100, 747–756. doi: 10.1093/aob/mcm156
Chapman, H. D. (1968). “The mineral nutrition of citrus,” in The Citrus Industry, Vol. 2, eds W. Reuther, H. J. Webber, and L. D. Batchelor (California, CA: Division of Agricultural Sciences, University of California), 127–189.
Chen, C., Song, Y., Zhuang, K., Li, L., Xia, Y., and Shen, Z. (2015). Proteomic analysis of copper-binding proteins in excess copper-stressed roots of two rice (Oryza sativa L.) varieties with different Cu tolerances. PLoS ONE 10:e0125367. doi: 10.1371/journal.pone.0125367
Demiray, H., Dereboylu, A. E., Altan, F., and Zeytünlüoğlu, A. (2011). Identification of proteins involved in excess boron stress in roots of carrot (Daucus carota L.) and role of niacin in the protein profiles. Afr. J. Biotechnol. 10, 15545–15551. doi: 10.5897/AJB11.1117
Erdal, S., Genc, E., Karaman, A., Khosroushahi, F. K., Kizilkaya, M., Demir, Y., et al. (2014). Differential responses of two wheat varieties to increasing boron toxicity. Changes on antioxidant activity, oxidative damage and DNA profile. J. Environ. Protect. Ecol. 15, 1217–1229.
Guo, P., Qi, Y. P., Yang, L. T., Ye, X., Huang, J. H., and Chen, L. S. (2016). Long-term boron-excess-induced alterations of gene profiles in roots of two citrus species differing in boron-tolerance revealed by cDNA-AFLP. Front. Plant Sci. 7:898. doi: 10.3389/fpls.2016.00898
Guo, P., Qi, Y. P., Yang, L. T., Ye, X., Jiang, H. X., Huang, J. H., et al. (2014). cDNA-AFLP analysis reveals the adaptive responses of citrus to long-term boron-toxicity. BMC Plant Biol. 14:284. doi: 10.1186/s12870-014-0284-5
Han, S., Tang, N., Jiang, H. X., Yang, L. T., Li, Y., and Chen, L. S. (2009). CO2 assimilation, photosystem II photochemistry, carbohydrate metabolism and antioxidant system of Citrus leaves in response to boron stress. Plant Sci. 176, 143–153. doi: 10.1016/j.plantsci.2008.10.004
Hassan, M., Oldach, K., Baumann, U., Langridge, P., and Sutton, T. (2010). Genes mapping to boron tolerance QTL in barley identified by suppression subtractive hybridization. Plant Cell Environ. 33, 188–198. doi: 10.1111/j.1365-3040.2009.02069.x
Hodges, D. M., DeLong, J. M., Forney, C. F., and Prange, R. K. (1999). Improving the thiobarbituric acid-reactive-substances assay for estimating lipid peroxidation in plant tissues containing anthocyanin and other interfering compounds. Planta 207, 604–611. doi: 10.1007/s004250050524
Huang, J. H., Cai, Z. J., Wen, S. X., Guo, P., Ye, X., Lin, G. Z., et al. (2014). Effects of boron toxicity on root and leaf anatomy in two citrus species differing in boron tolerance. Trees Struct. Funct. 28, 1653–1666. doi: 10.1007/s00468-014-1075-1
Huang, J. H., Qi, Y. P., Wen, S. X., Guo, P., Chen, X. M., and Chen, L. S. (2016). Illumina microRNA profiles reveal the involvement of miR397a in citrus adaptation to long-term boron toxicity via modulating secondary cell-wall biosynthesis. Sci. Rep. 6:22900. doi: 10.1038/srep22900
Jiang, H. X., Tang, N., Zheng, J. G., and Chen, L. S. (2009). Antagonistic actions of boron against inhibitory effects of aluminum toxicity on growth, CO2 assimilation, ribulose-1,5-bisphosphate carboxylase/oxygenase, and photosynthetic electron transport probed by the JIP-test, of Citrus grandis seedlings. BMC Plant Biol. 9:102. doi: 10.1186/1471-2229-9-102
Jin, L. F., Liu, Y. Z., Yin, X. X., and Peng, S. A. (2016). Transcript analysis of citrus miRNA397 and its target LAC7 reveals a possible role in response to boron toxicity. Acta Physiol. Plant. 38, 18. doi: 10.1007/s11738-015-2035-0
Karan, R., and Subudhi, P. K. (2012). Overexpression of a nascent polypeptide associated complex gene (SaβNAC) of Spartina alterniflora improves tolerance to salinity and drought in transgenic Arabidopsis. Biochem. Biophys. Res. Commun. 424, 747–752. doi: 10.1016/j.bbrc.2012.07.023
Kettles, N. L., Kopriva, S., and Malin, G. (2014). Insights into the regulation of DMSP synthesis in the diatom Thalassiosira pseudonana through APR activity, proteomics and gene expression analyses on cells acclimating to changes in salinity, light and nitrogen. PLoS ONE 9:e94795. doi: 10.1371/journal.pone.0094795
Kim, H. J., Balcezak, T. J., Nathin, S. J., McMullen, H. F., and Hansen, D. E. (1992). The use of a spectrophotometric assay to study the interaction of S-adenosylmethionine synthetase with methionine analogues. Anal. Biochem. 207, 68–72. doi: 10.1016/0003-2697(92)90501-W
Klychnikov, O. I., Li, K. W., Lill, H., and de Boer, A. H. (2007). The V-ATPase from etiolated barley (Hordeum vulgare L.) shoots is activated by blue light and interacts with 14-3-3 proteins. J. Exp. Bot. 58, 1013–1023. doi: 10.1093/jxb/erl261
Kwak, K. J., Kim, Y. O., and Kang, H. (2005). Characterization of transgenic Arabidopsis plants overexpressing GR-RBP4 under high salinity, dehydration, or cold stress. J. Exp. Bot. 56, 3007–3016. doi: 10.1093/jxb/eri298
Lee, K., Bae, D. W., Kim, S. H., Han, H. J., Liu, X., Park, H. C., et al. (2010). Comparative proteomic analysis of the short-term responses of rice roots and leaves to cadmium. J. Plant Physiol. 167, 161–168. doi: 10.1016/j.jplph.2009.09.006
Li, D., Liu, H., Zhang, H., Wang, X., and Song, F. (2008). OsBIRH1, a DEAD-box RNA helicase with functions in modulating defence responses against pathogen infection and oxidative stress. J. Exp. Bot. 59, 2133–2146. doi: 10.1093/jxb/ern072
Li, Y., Han, M. Q., Lin, F., Ten, Y., Lin, J., Zhu, D. H., et al. (2015). Soil chemical properties, ‘Guanximiyou’ pummelo leaf mineral nutrient status and fruit quality in the southern region of Fujian province, China. J. Soil Sci. Plant Nutr. 15, 615–628. doi: 10.4067/s0718-95162015005000029
Lukaszewicz, M., Matysiak-Kata, I., Aksamit, A., and Szopa, J. (2002). 14-3-3 protein regulation of the antioxidant capacity of transgenic potato tubers. Plant Sci. 163, 125–130. doi: 10.1016/S0168-9452(02)00081-X
Moffatt, B. A., Stevens, Y. Y., Allen, M. S., Snider, J. D., Pereira, L. A., Todorova, M. I., et al. (2002). Adenosine kinase deficiency is associated with developmental abnormalities and reduced transmethylation. Plant Physiol. 128, 812–821. doi: 10.1104/pp.010880
Papadakis, I. E., Dimassi, K. N., Bosabalidis, A. M., Therios, I. N., Patakas, A., and Giannakoula, A. (2004). Boron toxicity in ‘Clementine’ mandarin plants grafted on two rootstocks. Plant Sci. 166, 539–547. doi: 10.1016/j.plantsci.2003.10.027
Peng, H. Y., Qi, Y. P., Lee, J., Yang, L. T., Guo, P., Jiang, H. X., et al. (2015). Proteomic analysis of Citrus sinensis roots and leaves in response to long-term magnesium-deficiency. BMC Genomics 16:253. doi: 10.1186/s12864-015-1462-z
Ravanel, S., Gakière, B., Job, D., and Douce, R. (1998). The specific features of methionine biosynthesis and metabolism in plants. Proc. Natl. Acad. Sci. U.S.A. 95, 7805–7812. doi: 10.1073/pnas.95.13.7805
Sachetto-Martins, G., Franco, L. O., and de Oliveira, D. E. (2000). Plant glycine-rich proteins: a family or just proteins with a common motif? Biochim. Biophys. Acta 1492, 1–14. doi: 10.1016/s0167-4781(00)00064-6
Sang, W., Huang, Z. R., Qi, Y. P., Yang, L. T., Guo, P., and Chen, L. S. (2015). An investigation of boron-toxicity in leaves of two citrus species differing in boron-tolerance using comparative proteomics. J. Proteomics 123, 128–146. doi: 10.1016/j.jprot.2015.04.007
Schoor, S., Farrow, S., Blaschke, H., Lee, S., Perry, G., von Schwartzenberg, K., et al. (2011). Adenosine kinase contributes to cytokinin interconversion in Arabidopsis. Plant Physiol. 157, 659–672. doi: 10.1104/pp.111.181560
Shen, B., Li, C., and Tarczynski, M. C. (2002). High free-methionine and decreased lignin content result from a mutation in the Arabidopsis S-adenosyl-L-methionine synthetase 3 gene. Plant J. 29, 371–380. doi: 10.1046/j.1365-313X.2002.01221.x
Sheng, O., Zhou, G., Wei, Q., Peng, S., and Deng, X. (2010). Effects of excess boron on growth, gas exchange, and boron status of four orange scion-rootstock combinations. J. Plant Nutr. Soil Sci. 173, 469–476. doi: 10.1002/jpln.200800273
Smith, J. A., Uribe, E. G., Ball, E., Heuer, S., and Lüttge, U. (1984). Characterization of the vacuolar ATPase activity of the crassulacean-acid-metabolism plant Kalanchoë daigremontiana receptor modulating. Eur. J. Biochem. 141, 415–420. doi: 10.1111/j.1432-1033.1984.tb08207.x
Smith, T. E., Grattan, S. R., Grieve, C. M., Poss, J. A., Läuchli, A. E., and Suarez, D. L. (2013). pH dependent salinity-boron interactions impact yield, biomass, evapotranspiration and boron uptake in broccoli (Brassica oleracea L.). Plant Soil 370, 541–554. doi: 10.1007/s11104-013-1653-9
Tang, L., Kim, M. D., Yang, K. S., Kwon, S. Y., Kim, S. H., Kim, J. S., et al. (2008). Enhanced tolerance of transgenic potato plants overexpressing nucleoside diphosphate kinase 2 against multiple environmental stresses. Transgenic Res. 17, 705–715. doi: 10.1007/s11248-007-9155-2
Tanou, G., Job, C., Rajjou, L., Arc, E., Belghazi, M., Diamantidis, G., et al. (2009). Proteomics reveals the overlapping roles of hydrogen peroxide and nitric oxide in the acclimation of citrus plants to salinity. Plant J. 60, 795–804. doi: 10.1111/j.1365-313X.2009.04000.x
Tanou, G., Ziogas, V., Belghazi, M., Christou, A., Filippou, P., Job, D., et al. (2014). Polyamines reprogram oxidative and nitrosative status and the proteome of citrus plants exposed to salinity stress. Plant Cell Environ. 37, 864–885. doi: 10.1111/pce.12204
Wang, J., Nakazato, T., Sakanishi, K., Yamada, O., Tao, H., and Saito, I. (2006). Single-step microwave digestion with HNO3 alone for determination of trace elements in coal by ICP spectrometry. Talanta 68, 1584–1590. doi: 10.1016/j.talanta.2005.08.034
Wang, Z., Li, H., Ke, Q., Jeong, J. C., Lee, H. S., Xu, B., et al. (2014). Transgenic alfalfa plants expressing AtNDPK2 exhibit increased growth and tolerance to abiotic stresses. Plant Physiol. Biochem. 84, 67–77. doi: 10.1016/j.plaphy.2014.08.025
Yadav, S. K., Singla-Pareek, S. L., Reddy, M. K., and Sopory, S. K. (2005). Transgenic tobacco plants overexpressing glyoxalase enzymes resist an increase in methylglyoxal and maintain higher reduced glutathione levels under salinity stress. FEBS Lett. 579, 6265–6271. doi: 10.1016/j.febslet.2005.10.006
Yang, L. T., Qi, Y. P., Lu, Y. B., Guo, P., Sang, W., Feng, H., et al. (2013). iTRAQ protein profile analysis of Citrus sinensis roots in response to long-term boron-deficiency. J. Proteomics 93, 179–206. doi: 10.1016/j.jprot.2013.04.025
You, X., Yang, L. T., Lu, Y. B., Li, H., Zhang, S. Q., and Chen, L. S. (2014). Proteomic changes of citrus roots in response to long-term manganese toxicity. Trees Struct. Funct. 28, 1383–1399. doi: 10.1007/s00468-014-1042-x
Zeh, M., Casazza, A. P., Kreft, O., Roessner, U., Bieberich, K., Willmitzer, L., et al. (2001). Antisense inhibition of threonine synthase leads to high methionine content in transgenic potato plants. Plant Physiol. 127, 792–802. doi: 10.1104/pp.010438
Zhang, X. L., Zhou, J., Han, Z., Shang, Q., Wang, Z. G., Gu, X. H., et al. (2012). Active methyl cycle and transfer related gene expression in response to drought stress in rice leaves. Rice Sci. 19, 86–93. doi: 10.1016/S1672-6308(12)60026-2
Zhou, C. P., Qi, Y. P., You, X., Yang, L. T., Guo, P., Ye, X., et al. (2013). Leaf cDNA-AFLP analysis of two citrus species differing in manganese tolerance in response to long-term manganese-toxicity. BMC Genomics 14:621. doi: 10.1186/1471-2164-14-621
Zhu, M., Chen, G., Dong, T., Wang, L., Zhang, J., Zhao, Z., et al. (2015). SlDEAD31, a putative DEAD-Box RNA Helicase gene, regulates salt and drought tolerance and stress-related genes in tomato. PLoS ONE 10:e0133849. doi: 10.1371/journal.pone.0133849
Keywords: boron-toxicity, Citrus grandis, Citrus sinensis, 2-DE, proteome, roots
Citation: Sang W, Huang Z-R, Yang L-T, Guo P, Ye X and Chen L-S (2017) Effects of High Toxic Boron Concentration on Protein Profiles in Roots of Two Citrus Species Differing in Boron-Tolerance Revealed by a 2-DE Based MS Approach. Front. Plant Sci. 8:180. doi: 10.3389/fpls.2017.00180
Received: 23 September 2016; Accepted: 30 January 2017;
Published: 17 February 2017.
Edited by:Joshua L. Heazlewood, University of Melbourne, Australia
Reviewed by:Athanassios Molassiotis, Aristotle University of Thessaloniki, Greece
Dominique Job, Centre National de la Recherche Scientifique, France
Copyright © 2017 Sang, Huang, Yang, Guo, Ye and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Li-Song Chen, firstname.lastname@example.org