- Research article
- Open Access
cDNA-AFLP analysis reveals the adaptive responses of citrus to long-term boron-toxicity
BMC Plant Biologyvolume 14, Article number: 284 (2014)
Boron (B)-toxicity is an important disorder in agricultural regions across the world. Seedlings of `Sour pummelo' (Citrus grandis) and `Xuegan' (Citrus sinensis) were fertigated every other day until drip with 10 μM (control) or 400 μM (B-toxic) H3BO3 in a complete nutrient solution for 15 weeks. The aims of this study were to elucidate the adaptive mechanisms of citrus plants to B-toxicity and to identify B-tolerant genes.
B-toxicity-induced changes in seedlings growth, leaf CO2 assimilation, pigments, total soluble protein, malondialdehyde (MDA) and phosphorus were less pronounced in C. sinensis than in C. grandis. B concentration was higher in B-toxic C. sinensis leaves than in B-toxic C. grandis ones. Here we successfully used cDNA-AFLP to isolate 67 up-regulated and 65 down-regulated transcript-derived fragments (TDFs) from B-toxic C. grandis leaves, whilst only 31 up-regulated and 37 down-regulated TDFs from B-toxic C. sinensis ones, demonstrating that gene expression is less affected in B-toxic C. sinensis leaves than in B-toxic C. grandis ones. These differentially expressed TDFs were related to signal transduction, carbohydrate and energy metabolism, nucleic acid metabolism, protein and amino acid metabolism, lipid metabolism, cell wall and cytoskeleton modification, stress responses and cell transport. The higher B-tolerance of C. sinensis might be related to the findings that B-toxic C. sinensis leaves had higher expression levels of genes involved in photosynthesis, which might contribute to the higher photosyntheis and light utilization and less excess light energy, and in reactive oxygen species (ROS) scavenging compared to B-toxic C. grandis leaves, thus preventing them from photo-oxidative damage. In addition, B-toxicity-induced alteration in the expression levels of genes encoding inorganic pyrophosphatase 1, AT4G01850 and methionine synthase differed between the two species, which might play a role in the B-tolerance of C. sinensis.
C. sinensis leaves could tolerate higher level of B than C. grandis ones, thus improving the B-tolerance of C. sinensis plants. Our findings reveal some novel mechanisms on the tolerance of plants to B-toxicity at the gene expression level.
Althought boron (B) is a micronutrient element required for normal growth and development of higher plants, it is harmful to plants when present in excess. Whilst of lesser importance than B-deficiency (a widespread problem in many agricultural crops), B-toxicity is also an important problem in agricultural regions across the world, which citrus trees are cultivated -. Despite the importance of B-toxicity for crop productivity, the mechanisms by which plants respond to B-toxicity are poorly understood yet. Recently, increasing attention has been paid to plant B-toxicity as a result of the increased demand for desalinated water, in which the B level may be too high for healthy irrigation of crops .
Alteration of gene expression levels is an inevitable process of plants responding to environmental stresses. Kasajima and Fujiwara first investigated high B-induced changes in gene expression in Arabidopsis thaliana roots and rosette leaves using microarray, and identified a number of high B-induced genes, including a heat shock protein and a number of the multi-drug and toxic compound extrusion (MATE) family transporters . Hassan et al. preformed suppression subtractive hybridization on root cDNA from bulked B-tolerant and -intolerant doubled haploid barley lines grown under moderate B-stress and identified 111 upregulated clones in the tolerant bulk under B-stress, nine of which were genetically mapped to B-tolerant quantitative trait loci. An antioxidative response mechanism was suggested to provide an advantage in tolerating high level of soil B . Recently, Aquea et al. found that B-toxicity upregulated the expression of genes related to ABA signaling, ABA response and cell wall modification, and downregulated the expression of genes involved in water transporters in Arabidopsis roots, concluding that root growth inhibition was caused by B-toxicity-induced water-stress . Most research, however, has focused on roots and herbaceous plants (i.e., barley, A. thaliana), very little is known about the differential expression of genes in response to B-toxicity in leaves and woody plants.
Citrus belongs to evergreen subtropical fruit trees. In China, B-toxicity often occurs in citrus orchards from high level of B in soils and/or irrigation water and from inappropriate application of B fertilizer especially under low-rainfall conditions ,. During 1998-1999, Huang et al. investigated the nutrient status of soils and leaves from 200 `Guanximiyou' pummelo (Citrus grandis) orchards located in Pinghe, Zhangzhou, China. Up to 61.5% and 17.0% of orchards were excess in leaf B and soil water-soluble B, respectively . Previous studies showed that B-toxicity disturbed citrus plant growth and metabolism in multiple way, including interference of nutrient uptake , ultrastructural damage of roots and leaves -, inhibition of CO2 assimilation, photosynthetic enzymes and photosynthetic electron transport, decrease of chlorophyll (Chl), carotenoid (Car) and total soluble protein levels, affecting leaf carbohydrate metabolism and antioxidant system ,. However, our understanding of the molecular mechanisms underlying these processes in citrus is very limited. To our best knowledge, no high B-toxicity-induced changes in gene expression profiles have been reported in citrus plants to date. Here we investigated the effects of B-toxicity on growth, leaf CO2 assimilation, leaf concentrations of malondialdehyde (MDA), pigments and total soluble protein, root and leaf concentration of B, leaf concentration of phosphorus (P), and leaf gene expression profiles using cDNA-amplified fragment length polymorphism (cDNA-AFLP) in Citrus grandis and Citrus sinensis seedlings differing in B-tolerance . The aims of this study were to elucidate the adaptive mechanisms of citrus plants to B-toxicity and to identify B-tolerant genes.
Effects of B-toxicity on seedlings growth, B concentration in roots and leaves, and P concentration in leaves
Because B is phloem immobile in citrus plants, B-toxic symptoms first developed in old leaves. The typical visible symptom produced in B-toxic leaves was leaf burn (chlorotic and/or necrotic), which only occurred in C. grandis plants. In the later stages, B-toxic leaves shed premature. By contrast, almost no visible symptoms occurred in C. sinensis plants except for very few plants (Additional file 1).
B-toxicity-induced decreases in root, shoot and whole plant dry weights (DWs) were more pronounced in C. grandis than in C. sinensis seedlings (Figure 1A-C). Root DW decreased to a larger extent than shoot DW in response to B-toxicity, and resulted in a decrease in root DW/shoot DW ratio of both C. grandis and C. sinensis seedlings (Figure 1A-B and D).
B-toxicity increased B concentration in roots and leaves, especially in leaves and decreased P concentration in C. grandis leaves. No significant differences were found in root and leaf B concentration and leaf P concentration between the two species at each given B treatment except that B concentration was higher in B-toxic C. sinensis leaves than in B-toxic C. grandis ones (Figure 2).
Effects of B-toxicity on leaf gas exchange, pigments, total soluble protein and MDA
B-toxicity-induced decreases in both CO2 assimilation and stomatal conductance were higher in C. grandis than in C. sinensis leaves. Intercellular CO2 concentration increased in C. grandis leaves, but did not significantly change in C. sinensis leaves in response to B-toxicity. CO2 assimilation and stomatal conductance in control leaves did not differ between the two species, but were higher in B-toxic C. sinensis leaves than in B-toxic C. grandis ones. Intercellular CO2 concentration in control leaves was higher in C. sinensis than in C. grandis, but the reverse was the case in B-toxic leaves (Figure 3A-C).
B-toxicity decreased concentrations of Chl a + b and Car and ratio of Chl a/b in C. grandis and C. sinensis leaves. In control leaves, all the three parameters did not differ between the two species, but Chl a + b and Car concentrations were higher in B-toxic C. sinensis leaves than in B-toxic C. grandis ones (Figure 3E-G).
Leaf concentrations of total soluble protein and MDA were decreased and increased by B-toxicity in C. grandis leaves, respectively, but were not significantly affected in C. sinensis ones (Figure 3D and H).
B-toxicity-induced differentially expressed genes revealed by cDNA-AFLP
Here we used a total of 256 selective primer combinations to isolate the differentially expressed transcript-derived fragments (TDFs) from B-toxic leaves of two citrus species differing in B-tolerance. A representative picture of a silver-stained cDNA-AFLP gel showing B-toxicity-induced genes in C. grandis and C. sinensis leaves was presented in Additional file 2. As shown in Table 1, a total of 6050 clear and unambiguous TDFs were detected from the B-toxic leaves, with an average of 25.7 (15-40) TDFs for each primer combination. Among these TDFs, 932 TDFs only presented in C. grandis, 631 TDFs only presented in C. sinensis, and 4587 TDFs presented in the two species.
A total of 218 and 104 differentially expressed and reproducible TDFs were successfully obtained from B-toxic C. grandis and C. sinensis leaves, respectively. All these TDFs were re-amplified, cloned and sequenced. For C. grandis, 183 of fragments yielded usable sequence data. Aligment analysis showed 132 TDFs were homologous to genes encoding known, putative predicted, uncharacterized, hypothetical or unnamed proteins, and the remaining 51 TDFs showed no significant matches (Tables 1 and 2). Among these matched TDFs, 67 (50.8%) TDFs were up-regulated and 65 (49.2%) were down-regulated by B-toxicity. These TDFs were related to different biological processes such as cell transport (12.9%), lipid metabolism (2.3%), nucleic acid metabolism (12.9%), carbohydrate and energy metabolism (12.1%), protein and amino acid metabolism (25.0%), stress responses (6.1%), cell wall and cytoskeleton modification (6.8%), signal transduction (2.3%), other and unknown processes (19.7%) (Figure 4A). For C. sinensis leaves, 90 differentially expressed TDFs produced readable sequences (Tables 1 and 2), 68 of which displayed homology to genes encoding known, putative, hypothetical, uncharacterized or unnamed proteins. The remaining 22 TDFs had no database matches. Of these matched TDFs, 31 (45.6%) TDFs increased and 37 (54.4%) decreased in response to B-toxicity. These TDFs were involved in cell transport (8.8%), lipid metabolism (4.4%), nucleic acid metabolism (13.2%), carbohydrate and energy metabolism (20.6%), protein and amino acid metabolism (25.0%), stress responses (7.4%), cell wall and cytoskeleton modification (2.9%), signal transduction (1.5%), other and unknown processes (16.2%) (Figure 4B).
Validation of cDNA-AFLP data using qRT-PCR
In this study, nine TDFs from C. sinensis leaves and nine TDFs from C. grandis ones were selected for qRT-PCR analysis in order to validate their expression patterns obtained by cDNA-AFLP analysis. Except for two TDFs (i.e., TDFs #187_1 and 195_1), the expression profiles of all the TDFs obtained by qRT-PCR were in agreement with the expression patterns produced by cDNA-AFLP (Figure 5). This technique was thus validated in 88.9% of cases. In addition to gene family complexity, the changes in the intensity of individual bands in the cDNA-AFLP gels might be responsible for the discrepancies between qRT-PCR and cDNA-AFLP analysis.
C. sinensis displayed higher B-tolerance than C. grandis
Our results showed that the effects of B-toxicity on plant growth (Figure 1A-C), and leaf gas exchange, pigments, total soluble protein, MDA (Figure 3) and P (Figure 2C) were more pronounced in C. grandis than in C. sinensis seedlings, meaning that C. sinensis has higher B-tolerance than C. grandis. The present work, like that of the previous workers ,,, indicates that the major of B in B-toxic citrus plants was accumulated in the leaves (Figure 2A and B). As shown in Figure 2B, B concentration was not lower in C. sinensis than in C. grandis leaves regardless of B concentration in the nutrient solution, indicating that C. sinensis leaves may tolerate higher level of B. Similar result has been obtained by Huang et al. . Here we isolated 67 up-regulated and 65 down-regulated TDFs from B-toxic C. grandis leaves, whilst only 31 up-regulated and 37 down-regulated TDFs from B-toxic C. sinensis ones (Figure 4), suggesting that B-toxicity affects C. sinensis leaves gene expression less than C. grandis ones. These data also support above inference that C. sinensis leaves may tolerate higher level of B.
We found that CO2 assimilation was lower in toxic leaves than in control leaves, while stomatal conductance was not lower in the former (Figure 3A-C), implies that B-toxicity-induced inhibition of CO2 assimilation in two citrus species is primarily due to non-stomatal factors. Similar results have been obtained on B-toxic C. grandis and C. sinensis,, `Navelina' orange and `Clementine' mandarin plants grafted on sour orange and Swingle citrumelo rootstocks ,, Newhall and Skagg's Bonanza navel orange plants grafted on Carrizo citrange and trifoliate orange .
Leaf carbohydrate and energy metabolism
Since B-toxicity decreased CO2 assimilation (Figure 3A), genes involved in photosynthesis and related biological processes might be affected by B-toxicity. As expected, 16 TDFs in C. grandis leaves and 14 TDFs in C. sinensis ones related to carbohydrate and energy metabolism were altered under B-toxicity (Table 2 and Figure 4). We found that B-toxicity decreased the transcript level of ribulose-1,5-bisphosphate (RuBP) carboxylase/oxygenase (Rubisco) small subunit precursor (TDF #143_2) gene in C. grandis leaves (Table 2), which agrees with the previous report that B-toxicity decreased the activity of Rubisco in C. grandis leaves . Hudson et al. showed that the reduction of Rubisco concentration by anti-small subunit led to decreased photosynthesis in transgenic tobacco plants, but unchanged stomatal conductance . Also, the mRNA abundances of photosystem II (PSII) 32 kDa protein (PsbA, TDF #251_1), chloroplast PSII oxygen-evolving complex 23 kDa polypeptide (TDF #112_2) and NifU-like protein (TDF #239_4) genes were down-regulated in B-toxic C. grandis leaves (Table 2). Khan et al. reported that PsbA knockout tobacco plants lacked PSII activity, accompanied by promoted senescence . By using differential RNA interference (RNAi), Ishihara et al. demonstrated that PSII activity was linearly correlated with the total amount of PsbP (PSII 23 kDa protein) . Ifuku et al. reported that PsbP is essential for the regulation and stabilization of PSII in higher plants . Yabe et al. proposed that Arabidopsis chloroplastic NifU-like protein, which can act as a Fe-S cluster scaffold protein, was required for biogenesis of ferredoxin and photosystem I (PSI) . B-toxicity-induced decreases in the transcript levels of PsbA, chloroplast PSII oxygen-evolving complex 23 kDa polypeptide and NifU-like protein genes agree with our report that B-toxicity impaired the whole photosynthetic electron transport from PSII donor side up to the reduction of end acceptors of PSI in C. grandis leaves . By contrast, B-toxicity increased the transcript levels of chloroplast PSII oxygen-evolving complex 23 kDa polypeptide (TDF #112_2) and glyceraldehyde-3-phosphate dehydrogenase B (TDE #23_2) in C. sinensis leaves (Table 2). NADP-glyceraldehyde-3-phosphate dehydrogenase is one of the two chloroplast enzymes which catalyze the reduction of 3-phosphoglycerate to triose phosphate . However, the expression of Rubisco activase (TDF #6_4) gene in C. sinensis leaves decreased in response to B-toxicity (Table 2). Generally speaking, B-toxic C. sinensis leaves had higher expression levels of photosynthetic genes than B-toxic C. grandis ones. This might be responsible for the greater decrease in CO2 assimilation in B-toxic C. grandis leaves compared with B-toxic C. sinensis ones. It is noteworthy that the mRNA level of gene encoding sedoheptulose-1,7-bisphosphatase (SBPase, TDF #249_3), a key factor for the RuBP regeneration, was up-regulated in B-toxic leaves of the two citrus species (Table 2). Harrison et al. showed that a small decrease in SBPase activity caused a decline in CO2 assimilation by reducing the capacity for RuBP regeneration . Lefebvre et al. observed that transgenic tobacco plants over-expressing SBPase had enhanced photosynthesis and growth from an early stage in development . Wang reported that transgenic tomato plants over-expressing SBPase were more tolerance to low temperature and had higher photosynthetic capacity under low temperature . Therefore, the up-regulation of SBPase might be an adaptive response to B-toxicity.
As shown in Table 2, B-toxicity decreased leaf expression levels of three genes [i.e., ADP-glucose pyrophosphorylase (TDF #235_2) in C. sinensis, starch branching enzyme I (TDF #42_1) in C. grandis and glucose-1-phosphate adenylyltransferase large subunit 1 (TDF #59_2) in the two citrus species] related to starch biosynthesis, which agrees with the previous report that B-toxicity decreased starch concentration in C. grandis leaves .
B-toxicity increased the mRNA levels of three genes encoding citrate synthase (TDF #75-2), pyruvate dehydrogenase E1 component subunit beta (TDF #87_1) and aconitate hydratase 3 (TDF #33-2) in C. grandis leaves (Table 2), indicating that tricarboxylic acid cycle might be up-regulated in B-toxic C. grandis leaves. Similarly, the transcript level of a glycolysis gene encoding 2,3-bisphosphoglycerate-independent phosphoglycerate mutase (TDF #161_3) was enhanced in B-toxic C. sinensis leaves (Table 2). However, the mRNA levels of plastidial pyruvate kinase 3 (TDF #35_1) and aconitate hydratase 1 (TDF #33_2) genes were reduced in B-toxic C. sinensis leaves (Table 2). There is evidence showing that plastidic pyruvate kinase plays a key role in fatty acid synthesis by controlling the supply of ATP and pyruvate for de novo fatty acid synthesis in plastids . Thus, the fatty acid metabolism in B-toxic C. sinensis leaves might be impaired due to decreased plastidic pyruvate kinase.
In Arabidopsis, three NADPH: protochlorophyllide oxidoreductases (PORs), denoted as PORA, PORB, and PORC participate in mediating the light-dependent protochlorophyllide reduction . Pattanayak and Tripathy showed that over-expression of PORC in Arabidopsis led to coordinated up-regulation of gene/protein expression of several Chl biosynthetic pathway enzymes, thus enhancing Chl synthesis, and that the 1O2-mediated photo-oxidative damage in transgenic plants overexpressing PORC was minimal under high light stress . The observed lower transcript level of PORC (TDF #171_2) in B-toxic C. grandis and C. sinensis leaves (Table 2) agrees with the results that B-toxicity decreased the concentration of Chl a + b in citrus leaves (Figure 3E).
Cytochrome P450s play a key role in biotic and abiotic stresses. Transgenic tobacco and potato plants expressing cytochrome P450 with increased monooxygenase activity tolerated better oxidative stress after herbicide treatment . We found that B-toxicity increased the expression levels of genes encoding cytochrome P450 (TDF #5_1) and cytochrome P450 like protein (TDF #76-1) in C. grandis leaves (Table 2), which agrees with the previous report that some of the 49 cytochrome P450 genes in Arabidopsis were upregulated by biotic (i.e., Alternaria brassicicola and Alternaria alternata) and abiotic [i.e., drought, high salinity, low temperature, hormones, paraquat, rose bengal, UV stress (UV-C), mechanical wounding and heavy metal stress (CuSO4)] stresses . Thus, the up-regulation of cytochrome P450s in B-toxic C. grandis leaves might be an adaptive response. However, B-toxicity decreased the expression of cytochrome P450 in Arabidopsis roots .
Taken all together, we isolated eight up-regulated and eight down-regulated TDFs from B-toxic C. grandis leaves, and five up-regulated and nine down-regulated from B-toxic C. sinesnsis ones. Among these differentially expressed TDFs, only SBPase (TDF #249_3) and PORC (TDF #171_2) were similarly affected by B-toxicity in the two species (Table 2). These results demonstrated that the transcript profiles in the two species were differentially altered under B-toxicity.
Leaf lipid metabolism
Allene oxide synthase (AOS) and hydroperoxide lyase (HPL) branches of the oxylipin pathway, which are responsible for the production of jasmonates and aldehydes, respectively, participate in a range of stresses. Recently, Liu et al. showed that depletion of rice OsHPL3 greatly stimulated the jasmonic acid-governed defense response . Therefore, the AOS pathway and jasmonate level might be up-regulated in the B-toxic C. sinensis leaves due to decreased expression of fatty acid HPL (TDF #10_1; Table 2), thus contributing to B-tolerance. In addition, B-toxicity also affected the transcript levels of three genes [i.e., plastidial pyruvate kinase 3 (TDF #35_1), sugar-dependent1 (TDF #195_1) and 3-oxoacyl-reductase (TDF #233_3)] related to lipid metabolism in C. sinensis leaves (Table 2). Thus, lipid metabolism might be altered in B-toxic C. sinensis leaves.
Tang et al. reported that transgenic tobacco plants over-expressing acyl carrier protein (ACP)-1 (or expressing antisense ACP1) exhibited an increase (or decrease) in leaf concentrations of total lipids and the main fatty acids, and were more tolerant (or sensitive) to cold stress . Branen et al. showed that reduction of ACP4 by antisense RNA led to a decrease in total leaf lipids and decreased photosynthetic efficiency, and concluded that ACP4 might play a major role in the biosynthesis of fatty acids for chloroplast membrane development . The lower transcript level of gene encoding ACP1, chloroplastic-like (TDF #8_1) in B-toxic C. grandis leaves (Table 2) means that fatty acid biosynthesis in these leaves might be impaired. However, the expression of α/β-hydrolase domain-containing protein (TDF #194_1) and phospholipase-like protein (PEARLI 4) domain-containing protein (TDF #186_4) genes were up-regulated in B-toxic C. grandis leaves (Table 2).
Leaf nucleic acid metabolism
As shown in Table 2, eight up-regulated genes (TDFs #52_1, 49_1, 72_4, 120_1, 44_1, 159_2, 164_1 and 73_2) and nine down-regulated genes (TDFs #250_3, 157_2, 11_1, 71_3, 67_4, 10_4, 22_3, 104_1 and 68_2) were isolated from B-toxic C. grandis leaves, while only five up-regulated genes (TDFs #73_2, 250_3, 157_2, 60_1 and 131_1) and four down-regulated genes (TDFs #164_1, 71_3, 186_1 and 108_1) were identified in B-toxic C. sinensis leaves. Obviously, B-toxicity affected nucleic acid metabolism more in C. grandis leaves than in C. sinensis ones. This agrees with our inference that C. sinensis may tolerate higher level of B.
Leaf protein and amino acid metabolism
All these differentially expressed TDFs encoding chloroplatic translation initiation factor IF-2 (TDF #236_1) involved in promoting the binding of formylmethionyl-tRNA to 30'S ribosomal subunits, eukaryotic release factor 1-3 (TDF #117_4) involved in the termination step of protein synthesis, EMB1241 (At5g17710; TDF #93_3 ) related to protein folding and stabilization, Ankyrin repeat domain-containing protein (TDF #73_3) involved mainly in mediating protein-protein interactions, and ribosomal proteins [i.e., 50S ribosomal protein L15 (TDF #179_4), 30S ribosomal protein S17 (TDF #105_1), putative 60S ribosomal protein L6 (TDF #99_6) and 60S ribosomal protein L4_1 (TDF #186_2) ] related to mature ribosome assembly and translation processes except for SHEPHERD (TDF #93_2) involved in the correct folding and/or complex formation of CLAVATA (CLV) proteins , 60S ribosomal protein L23 (TDF #129_2) and 60S ribosomal protein L10B (TDF #161_1), were down-regulated in B-toxic C. grandis leaves (Table 2), indicating that B-toxicity impairs protein biosynthesis in C. grandis leaves ,. By contrast, only three down-regulated genes [30S ribosomal protein S17 (TDF #105_1), chaperonin 20 (TDF #98_1) involved in protein folding and stabilization and AT5G47880 (TDF #69_3) involved in the termination step of protein synthesis] were detected in B-toxic C. sinensis leaves (Table 2). These results demonstrated that B-toxicity affected protein biosynthesis more in the former than in the latter. This agrees with our data that B-toxicity only decreased total soluble protein concentration in C. grandis leaves (Figure 3H).
Here we observed four down-regulated genes [i.e., mitogen-activated protein (MAP) kinase (TDF #23_4), putative leucine-rich repeat receptor-like protein kinase (TDF #139-4), CBL-interacting protein kinase 19 (TDF #72_1) and At1g25390/F2J7_14 (TDF #39_3)] and three up-regulated genes [i.e., CDK activating kinase (TDF #12_2), serine/threonine protein kinase ATR (TDF #22_2) and receptor-like protein kinase (TDF #235_3) ] involved in phosphorylation and one up-regulated gene [i.e., protein phosphatase 2C (TDF #99_1)] involved in dephosphorylation in B-toxic C. grandis leaves, while only one down-regulated gene [i.e., receptor-like protein kinase (TDF #110_1)] and one up-regulated gene [i.e., protein phosphatase 2C (TDF #99_1)] in B-toxic C. sinensis leaves (Table 2). This means that C. sinensis leaves might achieve a better balance between phosphorylation and dephosphorylation than C. grandis ones under B-toxicity, which might contribute to the B-tolerance of C. sinensis.
Inactive (i.e., incorrect folding) and futile proteins for cell are tagged by ubiquitin for proteolysis . In this study, we found four up-regulated genes [i.e., C3H4 type zinc finger protein (TDF #99_2), AT5g57360/MSF19_2 (TDF #54_1), E3 ligase SAP5 (TDF #57_1) and root phototropism protein 2 (TDF #234_1)] and three down-regulated genes [i.e., E3 ubiquitin-protein ligase BRE1-like protein (TDF #96_1), Skp1-like protein 1 (TDF #187_1) and polyubiquitin (TDF #120_2)] involved in ubiquitination in B-toxic C. grandis leaves, and one up-regulated gene [i.e., putative E3 ubiquitin-protein ligase XBAT31 isoform 2 (TDF #158_2)] and two down-regulated genes [i.e., F-box family protein (TDF #73_1) and F-box with WD-40 2 (TDF #112_1)] involved in ubiquitination in B-toxici C. sinensis leaves. This indicates that ubiquitination might be involved in the adaptive response of citrus leaves to B-toxicity. Plant proteases has been shown to play key roles in controlling strict protein quality and degrading specific sets of proteins in response to environmental stresses . As expected, several genes (TDFs #38_3, 81_1, 38_4, 73_4, 240_1, 39_1, 145_1, 67_1 and 75_1) involved in proteolysis were altered in B-toxic C. grandis and C. sinensis leaves (Table 2).
S-adenosylmethionine (AdoMet) participates in a number of essential metabolic pathways in plants and is the principal biological methyl donor. AdoMet-dependent methylation is essential for keeping cellular functions in plants . Methionine synthase, which catalyzes the last reaction in de novo methionine synthesis, also serves to regenerate the methyl group of AdoMet. As shown in Table 2, B-toxicity increased the expression of AT4G01850 (TDF #138_4) involved in AdoMet biosynthesis in C. sinensis leaves, but decreased Methionine synthase expression (TDF #245_1) in C. grandis leaves, which might contribute to the higher tolerance of C. sinensis leaves to B-toxicity than that of C. grandis ones.
N-carbamoylputrescine amidase (TDF #213_4) involved in polyamine (putrescine) biosynthesis were down-regulated in B-toxic C. grandis leaves (Table 2). This means that the biosynthesis of polyamine might be inhibited in B-toxic C. grandis leaves, which disagrees with the previous report that 1000 μM B increased leaf concentration of putrescine in B-sensitive barley cultivar, but decreased its concentration in B-tolerant one .
The up-regulation of 2-oxoglutarate-dependent dioxygenase gene (TDF #61_2) in B-toxic C. grandis leaves (Table 2) agrees with the reports that B-toxicity stimulated the general amino acid control system in Saccharomyces cerevisiae and that the concentration of total amino acids in tomato leaves increased under B-toxicity . Evidence shows that 2-oxoglutarate-dependent dioxygenase participates in glucosinolate biosynthesis . Thus, the concentration of glucosinolates might be enhanced in B-toxic C. grandis leaves.
There is evidence showing that a few cystathionine-β-synthase (CBS) domain-containing proteins (CDCPs) play a role in plant stress response/tolerance and development . Overexpression of OsCBSX4 improved tobacco plant tolerance to salinity, oxidative, and heavy metal stresses . We observed that B-toxicity decreased the transcript level of CDCP (TDF #251_3) in C. sinensis leaves (Table 2), as obtained on manganese (Mn)-toxic C. grandis leaves . However, B-deficient C. sinensis roots had higher level of CBS family protein . Singh et al. observed that the expression of OsCBSX4 was up-regulated under high salinity, heavy metal, and oxidative stresses at seedling stage of a salt tolerant (Pokkali) rice cultivar, whilst its expression was upregulated only under NaCl stress, downregulated under heavy metal stress and kept unchanged under oxidative stress in a salt sensitive (IR64) rice one . Taken all together, the influence of stresses on expression of CDCP genes deponds on the kinds of stresses and plant species/cultivars.
Leaf stress responses
Inorganic pyrophosphatase (PPase), which cleaves pyrophosphate molecules to liberate two molecules of inorganic phosphate, are essential for the viability of organisms, because the removal of pyrophosphate, a by-product of a host of biosynthetic reactions, is required for preventing the inhibition of thermodynamically unfavorable reactions ,. George et al. observed that Nicotiana benthamiana plants lacking plastidial soluble PPase exhibited reduced drought tolerance as a result of the impaired leaf anabolic pathways . The up-regulation of PPase 1 (TDF #118_1) in B-toxic C. sinensis leaves (Table 2) might be an adaptive response to B-toxicity. By contrast, its expression (TDF #118_1) was down-regulated in B-toxic C. grandis leaves (Table 2).
Because leaf CO2 assimilation was decreased in B-toxic leaves (Figure 3A), less of the absorbed light energy was utilized in photosynthetic electron transport in these leaves, particularly under high light. Thus, reactive oxygen species (ROS) production might be enhanced in B-toxic leaves because of more excess absorbed photon flux . In addition to various ROS scavenger enzymes, "house-keeping" enzymes such as Nudix hydrolases (NUDXs) also play a role in ROS scavenging. Ogawa et al.  and Ishikawa et al.  showed that transgenic Arabidopsis plants overexpressing AtNUDX2 and AtNUDX7 exhibited higher tolerance to oxidative stress than wild type plants. Therefore, the higher expression level of NUDX19 (TDF #148_2) in B-toxic C. sinensis leaves might be an adaptive response to B-toxicity (Table 2). However, its expression level (TDF #148_2) in C. grandis leaves decreased in response to B-toxicity (Table 2).
Up to 10% of the ascorbate content of the whole leaf is localized in the apoplast, where it forms the first line of defense against external oxidants . In the apoplast, ascorbate oxidase (AO) oxidizes ascorbate to the unstable radical monodehydroascorbate which rapidly disproportionates to yield dehydroascorbate and ascorbate, thus participating in the regulation of the redox state of ascorbic acid pool. AO has been suggested to play a role in cell expansion via the modulation of redox control of the apoplast . Pignocchi et al.  showed that enhanced AO activity decreased the concentration and the redox state of ascorbic acid pool in the apoplast, whereas reduced AO activity increased its amount and redox state in the apoplast. Overexpression of AO in the apoplast of tobacco resulted in lowered capacity for scavenging ROS in the leaf apoplast accompanied by increased sensitivity to ozone . Fotopoulos et al.  observed that AO-overexpressing transgenic tobacco plants had increased sensitivity to various oxidative stress-promoting agents accompanied by a general suppression of the plant antioxidative metabolism. By contrast, a diminution in AO activity improved tomato yield under water deficit . The down-regulation of gene encoding Fe (II)/ascorbate oxidase family protein SRG1 (TDF #59_1) in B-toxic C. sinensis leaves (Table 2) might increase the amount and the redox state of AA pool in the apoplast, thus enhancing the B-tolerance.
Thioredoxins, which participates in supplying reducing power to reductases required for detoxifying lipid hydroperoxides or repairing oxidized proteins, play key roles in plant tolerance of oxidative stress . We found that the expression level of thioredoxin superfamily protein (TDF #137_2) was up-regulated in B-toxic C. grandis leaves (Table 2), indicating that thioredoxins might be involved in the ROS detoxification. However, the transcript level of thioredoxin superfamily protein (TDF #68_3) gene was down-regulated in B-toxic C. grandis leaves.
Our finding that B-toxicity increased the expression level of group 5 late embryogenesis abundant protein (LEA5, TDF #2_1) in C. grandis leaves (Table 2) agrees with the previous report that drought, heat and salt stresses stimulated the expression of LEA5 in citrus leaves . Accumulation of AtRAB28 (LEA5) protein in Arabidopsis through transgenic approach improved the germination rate under standard conditions or salt and osmotic stresses and the cation toxicity tolerance . Also, B-toxicity increased the transcript level of thaumatin-like protein 1 (TLP1, TDF #125_1) in C. grandis leaves (Table 2). The family of thaumatin-like proteins (also designated PR-5), which comprises proteins with various functions, is induced by biotic and abiotic factors in plants . Therefore, the up-regulation of LEA5 and TLP1 in B-toxic C. grandis leaves might be an adaptive response.
Protein sodium-and lithium-tolerant 1 (SLT1) gene isolated from tobacco (NtSLT1) and A. thaliana (AtSLT1) has been implicated in mediating salt tolerance by regulating Na+ homeostasis via the calcineurin (CaN) and SPK1/HAL4 (SPK1/HAL4 which encodes a serine-threonine kinase) signal transduction . Later, Antoine et al.  showed that rice OsSLT1 had molecular chaperone activity in vitro, and that OsSLT1 could be an important component of the cell immediate defenses against possible protein denaturation and aggregation. The down-regulation of SLT1 (TDF #99_5) in B-toxic C. grandis leaves (Table 2) means that Na+ homeostasis or related processes mediated by SLT1 are impaired in B-toxic C. grandis leaves.
Plant autophagy plays a role in various stress responses, pathogen defense, and senescence . Xiong et al. , showed that AtATG18a was necessary for the formation of autophagosomes during nutrient stress and senescence in A. thaliana and that autophagy participated in the degradation of oxidized proteins under oxidative stress conditions in Arabidopsis. AtATG18a RNAi plants usually senesce earlier and have lower tolerance to various stresses including drought, salt and oxidative stresses compared with wild-type plants ,. Our result showed that the transcript level of transducin/WD40 domain-containing protein (ATG18a, TDF #104_3) in C. sinensis leaves decreased in response to B-toxicity (Table 2), indicating that autophagy is impaired in C. sinensis leaves.
As shown in Table 2, B-toxicity down-regulated the expression of "cold-regulated" gene (cold regulated 314 thylakoid membrane 2, TDF # 109_1) in C. sinensis leaves and universal stress protein A-like protein (TDF #150_2) in C. grandis leaves (Table 2), indicating that B-toxicity might affect the tolerance of plants to other stresses.
Leaf signal transduction
Here four genes involved in signal transduction were altered by B-toxicity (Table 2 and Figure 4). Evidence shows that that signal recognition particle 54 kDa protein (SRP54) plays important roles in chloroplast development ,. The down-regulation of signal recognition particle 54 kDa protein 2 (TDF #182_2) in B-toxic C. grandis leaves (Table 2) means that the biosynthesis of Chl is impaired in these leaves. This agrees with our results that B-toxicity affected Chl more in C. grandis leaves than in C. sinensis ones (Figure 3E).
Increasing evidence shows that 14-3-3 proteins play an important role in plant stress responses ,. The most direct evidence for the role of 14-3-3 proteins in stress responses comes from transgenic rice plants over-expressing ZmGF14-6 encoding a maize 14-3-3 protein  and cotton plants over-expressing Arabidopsis 14-3-3 . These transgenic plants displayed enhanced tolerance to drought stress. Heterotrimeric GTP-binding proteins (G proteins, consisting of subunits Gα, Gβ, and Gγ) are signaling molecules required for various eukaryotic organisms. Joo et al.  observed that A. thaliana mutant plants losing the Gβ protein were less tolerant to O3 damage than wild-type plants. Thus, the B-tolerance of C. grandis leaves might be down-regulated due to decreased transcript level of genes encoding 14-3-3 protein (TDF #108_2) and heterotrimeric GTP-binding protein subunit beta 1 (TDF #200_1) (Table 2).
In higher plants, the endogenous circadian clock is involved in the manipulation of different various cellular processes ranging from photosynthesis to stress responses ,. It also confers plants with competitive advantages, including improved photosynthesis, growth and survival . Nakamichi et al.  observed that A PRR9, 7 and 5 triple mutant of Arabidopsis had higher tolerance against drought, salt and cold stresses compared to wild type, demonstrating the involvement of the three genes in abiotic stress responses as negative regulators. The up-regulation of pseudo-response regulator 5 (PRR5; TDF #70_2) in B-toxic C. sinensis leaves (Table 2) agrees with the previous reports that PRR5 was induced by cold treatment in apical shoots of cassava  and in Arabidopsis leaves . Fukushima et al.  showed that PRR9, 7 and 5 negatively regulated the biosynthetic pathways of Chl, Car, ABA and α-tocopherol. This agrees with our results that B-toxici C. sinensis leaves had decreased concentrations of Chl a + b and Car (Figure 3E and H).
Leaf cell transport
As shown in Table 2 and Figure 4, the number of differentially expressed TDFs involved in cell transport was far less in B-toxic C. sinensis leaves than in B-toxic C. grandis ones, meaning that cell transport is less affected in the former than in the latter, which agrees with our inference that C. sinensis leaves may tolerate higher level of B.
Most of the differentially expressed TDFs (TDFs #26_1, 124_3, 66_1, 97_1, 53_1, 210_1, 178_1, 49_3, 137_1, 63_1, 51_1, 250_2, 79_2 and 67_3) associated with cell transport were up-regulated in B-toxic C. grandis leaves except for AT5g24810/F6A4_20 (TDF #89_2), protein transport protein SEC61 γ subunit (TDF #6_1) and putative β-subunit of adaptor protein complex 3, PAT2 (TDF #249_2) (Table 2), indicating that cell transport might be enhanced in B-toxic C. grandis leaves. Plasma-membrane H+-ATPase plays a crucial role in the plant response to environmental stresses, such as salt stress, aluminum (Al) stress, P and potassium (K) deficiencies . Wu et al.  reported that pumping of Ca2+ and Mn2+ by an endoplasmic reticulum-type Ca2+-ATPase (ECA1) into the endoplasmic reticulum was necessary for maintaining plant growth under calcium (Ca)-deficiency or Mn-toxicity. The PIB-ATPases (also known as heavy metal ATPases), which are involved in heavy metal transport across cellular membranes, play a crucial role in metal homeostasis and detoxification in plants . Proton pump interactor 1 (PPI1), an interactor of plasma-membrane H+-ATPase, stimulates its activity in vitro. The up-regulation of PPI1 (TDF #97_1) in B-toxic leaves agrees with our data that the transcript level of H+-ATPase 6 (TDF #26_1) in C. grandis leaves increased in response to B-toxicity (Table 2) and with the report that the expression of PPI1 in potato tuber was up-regulated by salt stress and cold .
ATP-binding cassette (ABC) transporters are involved in metal ion efflux from the plasma-membrane. AtPDR8, an ABC transporter localized in the plasma-membrane of A. thaliana root hairs and epidermal cells, confers metal tolerance . Our finding that the expression of ABC transporter G family member 40 (TDF #53_1) gene was up-regulated in B-toxic C. grandis leaves agrees with the reports that AtPDR8 in Arabidopsis roots and shoots was induced when exposed to copper (Cu), cadmium (Cd) and lead (Pb) , and that ABC transporter G family member 40 gene and ABC transporter A family member 7 gene were induced in drought-sensitive and -tolerant genotypes of Gossypium herbaceum, respectively under drought stress . However, the expression of AT5g24810/F6A4_20 (TDF #89_2) was down-regulated in B-toxic C. grandis leaves (Table 2).
Cu transporters (COPTs/Ctrs) are involved in the maintenance of Cu homeostasis in plants. Generally speaking, COPTs/Ctrs are up-regulated by Cu deprivation and down-regulated by Cu excess . COPT1 antisense Arabidopsis plants have decreased Cu level due to decreased Cu uptake and display sensitivity to Cu chelators . The up-regulation of COPT (TDF #210_1) in B-toxic C. grandis leaves might play a role in the maintenance of leaf Cu homeostasis.
Plant cyclic nucleotide gated channels (CNGCs) paly a role in heavy metal homeostasis. Previous study showed that transgenic tobacco plants overexpressing a truncated NtCBP4 (tobacco CNGC) had higher tolerance to Pb compared with wild type . Chan et al.  reported that cngc2 Arabidopsis mutants were hypersensitive to increased soil Ca. However, transgenic tobacco plants overexpressing NtCBP4 were hypersensitivity to Pb . B-toxicity-induced increase in transcript level of CNGC1 (TDF #178_1) in C. grandis leaves (Table 2) agrees with the report that the expression of AtCNGC2 was induced during Arabidopsis leaf senescence and AtCNGC2 might be involved in programmed cell death .
Membrane traffic is not only required for plant normal cellular function and maintenance of cellular viability, but also plays an important roles in plant responses to the environment ,. The transcript levels of genes [i.e., vacuolar-sorting receptor 3 (TDF #49_3), vacuolar protein-sorting-associated protein 37-1 (TDF #137_1), vesicle-associated membrane protein-associated protein (TDF #63_1), secY protein transport family protein (TDF #51-1), fat-free-like protein (TDF #250_2) and non-specific lipid-transfer protein (TDF #79_2)] involved in membrane traffic increased in B-toxic C. grandis leaves except for genes encoding protein transport protein SEC61 γ subunit (TDF #6_1) and putative β-subunit of adaptor protein complex 3, PAT2 (TDF # 249_2) (Table 2). This indicates that membrane traffic might be enhanced in B-toxic C. grandis leaves.
Plant sieve element occlusion (SEO) genes have been shown to encode the common phloem proteins (P-proteins) that plug sieve plates after wounding. Tobacco SEO-RNA interference lines were essentially devoid of P-protein structures and lost photoassimilates more rapidly after injury than control plants . Therefore, the up-regulation of sieve element occlusion protein 1 gene (TDF #67_3) in B-toxic C. grandis leaves (Table 2) might be of advantage to prevent the loss of photoassimilates. Recently, Huang et al. observed that many electron-dense particles deposited near sieve plates of B-toxic C. grandis and C. sinensis leaves . In conclusion, the up-regulation of cell transport in B-toxic C. grandis leaves might be an adaptive response of plants to B-toxicity.
By contrast, we isolated three down-regulated [i.e., putative β-subunit of adaptor protein complex 3, PAT2 (TDF #249_2), sugar transporter ERD6-like 5 (TDF #61_1) and metal tolerance protein (MTP, TDF #179_2)] and three up-regulated [i.e., sieve element occlusion protein 1 (TDF #67_3), kinesin-related protein (TDF #51_4) and bidirectional sugar transporter SWEET7 (TDF #36_2) TDFs from B-toxic C. sinensis leaves (Table 2). Generally speaking, cell transport might be not enhanced in B-toxicity leaves.
In plants, kinesins are involved in a variety of cellular processes including intracellular transport, spindle assembly, phragmoplast assembly, chromosome motility, MAP kinase regulation and microtubule stability . Li et al.  reported that mutation of rice BC12/GDD1 encoding a kinesin-like protein led to dwarfism with impaired cell elongation. Nishihama et al.  demonstrated that the expansion of the cell plate in tobacco plant cytokinesis required kinesin-like proteins (i.e., NACK1 and NACK2) to regulate the activity and localization of MAP kinase kinase kinase. Therefore, the up-regulation of kinesin-like protein (TDF #51_4) in C. sinensis leaves (Table 2) might be an adaptive responsive to B-toxicity. However, the transcript level of putative β-subunit of adaptor protein complex 3, PAT2 (TDF #249_2) in C. sinensis leaves decreased in response to B-toxicity (Table 2).
Plant SWEETs function as facilitators involved in the influx and the efflux of sugar into and out of cells . We found that the expression level of SWEET7 (TDF #36_2) in C. sinensis leaves increased in response to B-toxicity (Table 2), which agrees with the previous report that SWEET15/SAG29 was enhanced in senescing Arabidopsis leaves . However, the expression of gene encoding sugar transporter ERD6-like 5 (TDF #61_1), a passive facilitator for the diffusion of glucose across the tonoplast membrane, was down-regulated in B-toxic C. sinensis leaves (Table 2). This disagrees with the previous report that the expression of AtESL1 (ERD six-like 1) was induced by various stresses including drought, high salinity and ABA in Arabidopsis plants .
MTPs are a subfamily of the cation diffusion facilitator (CDF) family found in plants. So far, most studied CDF family members confer heavy metal tolerance by affecting heavy metal efflux from the cytoplasm . The down-regulation of MTP (TDF #179_2) in C. sinensis leaves (Table 2) means that the tolerance of plants to heavy metal might be reduced in B-toxic plants. This agrees with our previous report that the tolerance of C. grandis plants to Al-toxicity was higher under adequate B supply than under excess B .
Leaf cell wall and cytoskeleton modification
Eleven TDFs associated with cell wall and cytoskeleton modification were altered by B-toxicity (Table 2 and Figure 4). O-methyltransferase (OMT) genes are involved in lignin biosynthesis. Fu et al.  showed that down-regulation of the caffeic acid 3-O-methyltransferase (COMT) gene in switchgrass lowered lignin level in whole tillers and stems of transgenic plants and enhanced forage quality. Transgenic Leucaena leucocephala plants expressing antisense OMT displayed decreased activity of OMT activity and concentration of lignin . Therefore, the biosynthesis of lignin in B-toxic C. grandis and C. sinensis leaves might be reduced due to decreased expression of COMT (TDF #49_4) and COMT3 (TDF #125_2) (Table 2). In addition, the biosynthesis of chitin in C. grandis and C. sinensis leaves and cellulose in C. grandis leaves might be down-regulated under B-toxicity due to the down-regulation of chitinase (TDF #10_3) and cellulose synthase (TDF #249_4) (Table 2). These results demonstrated that B-toxicity might impair citrus cell wall metabolism, which agrees with the previous suggestion that leaf cupping, a specific visible B-toxic symptom in some species might be due to the inhibition of cell wall expansion, through disturbance of cell wall crosslinks . However, the transcript levels of genes encoding OMT1 (TDF #33_3), LIM domain-containing protein (TDF #241_1), UDP-glucose flavonoid 7-O-glucosyltransferase (TDF #124_2), UDP-glucosyltransferase family 1 protein (TDF #3_3), limonoid UDP-glucosyltransferase (TDF #70_4) and putative glucosyltransferase (TDF #63_2) in C. grandis increased in response to B-toxicity (Table 2).
Evidence shows that lily LIM1  and all Arabidopsis LIM domain proteins  participate in regulating actin cytoskeleton organization and dynamics. Tobacco LIM1 protein acts in the cytoplasm as an actin binding and bundling protein  and in the nucleus as a transcription factor regulating the expression of genes related to lignin biosynthesis . Recently, Moes et al.  demonstrated the involvement of tobacco LIM2 in actin-bundling and histone gene transcription. The up-regulation of LIM domain-containing protein (TDF #241_1) in B-toxic C. grandis leaves (Table 2) agrees with the report that the expression of LIM domain-containing protein in Physcomitrella patens increased under cold acclimation .
Glycosyltransferases (GTs), which catalyze the formation of glycosidic bonds between donor sugars and acceptor molecules, participate in many aspects of a plant life, including cell wall biosynthesis ,. In Arabidopsis, up to 10 or 12 GT2 family members form the cellulose synthase catalytic subunit and callose synthase gene families . In plants, UDP-glucosyltransferases (UGTs) have been suggested to play important roles in keeping cell homeostasis, regulating plant growth and improving their tolerance to environmental stresses . Overexpression of UGT74E2 conferred tolerance to salinity and drought stresses in A. thaliana. Transgenic tobacco plants overexpressiong UGT85A5 exhibited enhanced salt tolerance . Therefore, the up-regulation of UDP-glucose flavonoid 7-O-glucosyltransferase (TDF #124_2), UGT family 1 protein (TDF #3_3), limonoid UGT (TDF #70_4) and putative GT (TDF #63_2) genes in B-toxic C. grandis leaves (Table 2) might play a role in B-tolerance of plants. However, loss of function of a UGT73B2 alone or in conjunction with UGT73B1 and UGT73B3 resulted in enhanced oxidative stress tolerance in Arabidopsis, whilst transgenic Arabidopsis plants overexpressing UGT73B2 displayed decreased oxidative stress tolerance .
Overexpression of bacterial or plant gene encoding phytoene synthase (PSY), a key regulatory enzyme in Car biosynthesis, led to enhanced level of total Car in various higher plants ,. Transgenic Arabidopsis plants overexpressing PSY from euhalophyte Salicornia europaea had higher tolerance to salt stress than wild type plants by enhanced photosynthetic efficiency and antioxidative capacity . Cidade et al.  showed that ectopic expression of PSY from Citrus paradisi fruit conferred abiotic stress tolerance in transgenic tobacco, which was correlated with the increased endogenous ABA level and expression of stress-responsive genes. Our finding that B-toxic C. grandis leaves had lower transcript of PSY (TDF #229_4; Table 2) means that the biosynthesis of Car and the antioxidative capacity may be decreased in B-toxic leaves. This agrees with our data that B-toxicity affected Car more in C. grandis leaves than in C. sinensis one (Figure 3G) and the inference that C. grandis may tolerate lower level of B.
Strictosidine synthase (Str), a key enzyme in alkaloid biosynthesis, catalyzes the condensation of tryptamine and secologanin leading to the synthesis of numerous monoterpenoid indole alkaloids in higher plants . The up-regulation of Str family protein gene (TDF #231_1) in B-toxic C. sinensis leaves (Table 2) agrees with the previous report that Str in Catharanthus roseus leaves was enhanced under dehydration, salt and UV stresses  and that B-toxicity decreased IAA level in Triticum durum seedlings , because the expression of Str was inhibited by auxin . B-toxicity-induced up-regulation of Str family protein gene (TDF #231_1) also agrees with our reference that the AOS pathway and jasmonate level might be up-regulated in the B-toxic C. sinensis leaves due to decreased expression of fatty acid HPL gene (TDF #10_1) (Table 2), because Str has been shown to be induced by jasmonate . By contrast, the expression of Str family protein gene (TDF #231_1) was down-regulated in B-toxic C. grandis leaves (Table 2), which agrees with the previous report that cold stress led to Str down-regulation in C. roseus leaves .
B-toxicity affected C. grandis seedling growth, leaf CO2 assimilation, pigments, total soluble protein, MDA and P more than C. sinensis, indicating that C. sinensis have higher B-tolerance than C. grandis ones. Under B-toxicity, C. sinensis leaves accumulated more B than C. grandis ones, meaning that the former may tolerate higher level of B. Using cDNA-AFLP, we successfully isolated 67 up-regulated and 65 down-regulated TDFs from B-toxic C. grandis leaves, whilst only 31 up-regulated and 37 down-regulated TDFs from B-toxic C. sinensis ones. This indicates that gene expression is less affected in B-toxic C. sinensis leaves than in C. grandis ones, which might be associated with the fact that C. sinensis leaves can tolerate higher level of B. The higher B-tolerance of C. sinensis might be related to the findings that B-toxic C. sinensis leaves had higher expression levels of genes involved in photosynthesis, which might contribute to the higher photosynthesis and light utilization and less excess light energy compared to the B-toxic C. grandis ones, and in ROS scavenging, thus preventing them from photo-oxidative damage. In addition, B-toxicity-induced alteration in the expression levels of genes encoding inorganic PPase 1, AT4G01850 and methionine synthase differed between the two species, which might also contribute to the B-tolerance of C. sinensis. In this study, a total of 174 differentially expressed TDFs were isolated from two citrus species, only 26 TDFs presented in the two citrus, the remaining TDFs presented only in C. grandis or C. sinensis, demonstrating that the B-toxicity-responsive genes differ between the two citrus species. For example, cell transport were up-regulated in B-toxicity C. grandis leaves, whilst this did not occur in B-toxic C. sinensis ones.
This study was conducted from February to December, 2011 at Fujian Agriculture and Forestry University. Plant culture and B treatments were performed according to Han et al. . Briefly, 5-week-old uniform seedlings of `Xuegan' (Citrus sinensis) and `Sour pummelo' (Citrus grandis) were transplanted to 6 L pots containing fine river sand. Plants, two per pot, were grown in a greenhouse under natural photoperiod at Fujian Agriculture and Forestry University. Eight weeks after transplanting, each pot was supplied every other day until dripping with nutrient solution containing 10 μM (control) or 400 μM (B-toxic) H3BO3 and 6 mM KNO3, 4 mM Ca (NO3)2, 2 mM NH4H2PO4, 1 mM MgSO4, 10 μM H3BO3, 2 μM MnCl2, 2 μM ZnSO4, 0.5 μM CuSO4, 0.065 μM (NH4)6Mo7O24 and 20 μM Fe-EDTA for 15 weeks. At the end of the experiment, fully expanded leaves from different replicates and treatments were used for all the measurements. Leaves were collected at noon under full sun and immediately frozen in liquid nitrogen and were stored at −80°C until extraction.
Measurements of plant DW, root and leaf B, leaf P, total soluble protein, MDA and pigments
Ten plants per treatment from different pots were harvested and divided into their parts (roots and shoots). The plant parts were then dried at 75°C for 48 h and their DWs measured. B concentration in roots and leaves was assayed by ICP emission spectrometry after microwave digestion with HNO3. Leaf P concentration was measured according to Ames . Leaf total soluble protein was measured according to Bradford  using bovine serum albumin as standard after being extracted with 50 mM Na2HPO4-KH2PO4 (pH 7.0) and 5% (w/v) insoluble polyvinylpyrrolidone. Extraction and determination of leaf MDA were performed according to Hodges et al. . Chl, Chl a, Chl b and Car were assayed according to Lichtenthaler  after being extracted with 80 (v/v) actetone.
Measurements of leaf gas exchange
Leaf gas exchange was measured using a CIARS-2 portable photosynthesis system (PP systems, Herts, UK) at ambient CO2 concentration under a controlled light intensity of 990-1010 μMol m-2 s-1 between 9:00 and 11:00 on a clear day. During measuring, leaf temperature and air relative humidity were 32.2 ± 0.2°C and 66.6 ± 0.8%, respectively.
Leaf RNA extraction, cDNA synthesis and cDNA-AFLP analysis
Total RNA was extracted from ca. 300 mg of frozen mixed leaves from B-toxic and control plants of C. grandis and C. sinensis using Recalcirtant Plant Total RNA Extraction Kit (Centrifugal column type, Bioteke Corporation, China). There were three biological replicates for each treatment. Leave of 4-5 plants from different pots were mixed as a biological replicate. Equal amounts of leaves were collected from each plant. cDNA synthesis and cDNA-AFLP analysis were performed according to Zhou et al. .
Validation of cDNA-AFLP data using qRT-PCR
Total RNA was extracted from the frozen leaves as described above. qRT-PCR analysis was performed according to Zhou et al. . Specific primers were designed from the sequences of 16 differentially expressed TDFs using Primer Primier Version 5.0 (PREMIER Biosoft International, CA, USA). The sequences of the F and R primers used were listed in Additional file 3. Samples for qRT-PCR were run in 3 biological replicates with 3 technical replicates. Leave of 4-5 plants from different pots were mixed as a biological replicate. Relative gene expression was calculated using ddCt algorithm. For the normalization of gene expression, citrus actin (GU911361.1) was used as an internal standard and the leaves from control plants were used as reference sample, which was set to 1.
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. Results represented the mean ± SE. Statistical analyses of data were carried out by ANOVA tests. Means were separated by the least significant difference test at P <0.05 level.
Acyl carrier protein
Allene oxide synthase
CBS domain-containing protein
Cation diffusion facilitator
cDNA-amplified fragment length polymorphism
Cyclic nucleotide gated channel
Caffeic acid 3-O-methyltransferase
Group 5 late embryogenesis abundant protein
Multi-drug and toxic compound extrusion
Metal tolerance protein
Proton pump interactor 1
Pseudo-response regulator 5
PSII 32 kDa protein
PSII 23 kDa protein
Reactive oxygen species
Sieve element occlusion
Protein sodium-and lithium-tolerant 1
Nable RO, Banuelos GS, Paull JG: Boron toxicity. Plant Soil. 1997, 193: 181-198.
Papadakis IE, Dimassi KN, Therios IN: Response of two citrus genotypes to six boron concentrations: concentration and distribution of nutrients, total absorption, and nutrient use efficiency. Aust J Agr Res. 2003, 54: 571-580.
Cervilla LM, Blasco B, Ríos J, Romero L, Ruiz J: Oxidative stress and antioxidants in tomato (Solanum lycopericum) plants subjected to boron toxicity. Ann Bot. 2007, 100: 747-756.
Parks JL, Edwards M: Boron in the environment. Crit Rev Environ Sci Technol. 2005, 35: 81-114.
Kasajima I, Fujiwara T: Identification of novel Arabidopsis thaliana genes which are induced by high levels of boron. Plant Biotech. 2007, 24: 355-360.
Hassan M, Oldach K, Baumann U, Langridge P, Sutton T: Genes mapping to boron tolerance QTL in barley identified by suppression subtractive hybridization. Plant Cell Environ. 2010, 33: 188-198.
Aquea F, Federici F, Moscoso C, Vega A, Jullian P, Haseloff J, Arce-Johnson P: A molecular framework for the inhibition of Arabidopsis root growth in response to boron toxicity. Plant Cell Environ. 2012, 35: 719-734.
Chen LS, Han S, Qi YP, Yang LT: Boron stresses and tolerance in citrus. Afr J Biotech. 2012, 11: 5961-5969.
Sheng O, Zhou GF, Wei QJ, Peng SA, Deng XX: Effects of excess boron on growth, gas exchange, and boron status of four orange scion-rootstock combinations. J Plant Nutr Soil Sci. 2010, 173: 469-476.
Huang YZ, Li J, Wu SH, Pang DM: Nutrition condition of the orchards in the main production areas of Guanxihoney pomelo trees (Pinhe county). J Fujian Agri Univ. 2001, 30: 40-43.
Papadakis IE, Dimassi KN, Bosabalidis AM, Therios IN, Patakas A, Giannakoula A: Effects of B excess on some physiological and anatomical parameters of `Navelina' orange plants grafted on two rootstocks. Environ Exp Bot. 2004, 51: 247-257.
Papadakis IE, Dimassi KN, Bosabalidis AM, Therios IN, Patakas A, Giannakoula A: Boron toxicity in `Clementine' mandarin plants grafted on two rootstocks. Plant Sci. 2004, 166: 539-547.
Huang JH, Cai ZJ, Wen SX, Guo P, Ye X, Lin GZ, Chen LS: Effects of boron toxicity on root and leaf anatomy in two citrus species differing in boron tolerance. Trees Struct Funct 2014. doi:10.1007/s00468-014-1075-1.,
Han S, Tang N, Jiang HX, Yang LT, Li Y, Chen LS: CO2 assimilation, photosystem II photochemistry, carbohydrate metabolism and antioxidant system of citrus leaves in response to boron stress. Plant Sci. 2009, 176: 143-153.
Sheng O, Song SW, Chen YJ, Peng SA, Deng XX: Effects of exogenous B supply on growth, B accumulation and distribution of two navel orange cultivars. Trees Struct Funct. 2009, 23: 59-68.
Hudson GS, Evans JR, von Caemmerer S, Arvidsson YBC, Andrewset TJ: Reduction of ribulose-1,5-bisphosphate carboxylase/oxygenase content by antisense RNA reduces photosynthesis in transgenic tobacco plants. Plant Physiol. 1992, 98: 294-302.
Khan MS, Hameed W, Nozoe M, Shiina T: Disruption of the psbA gene by the copy correction mechanism reveals that the expression of plastid-encoded genes is regulated by photosynthesis activity. J Plant Res. 2007, 120: 421-430.
Ishihara S, Yamamoto Y, Ifuku K, Sato F: Functional analysis of four members of the PsbP family in photosystem II in Nicotiana tabacum using differential RNA interference. Plant Cell Physiol. 2005, 46: 1885-1893.
Ifuku K, Yamamoto Y, Ono TA, Ishihara S, Sato F: PsbP protein, but not PsbQ protein, is essential for the regulation and stabilization of photosystem II in higher plants. Plant Physiol. 2005, 139: 1175-1184.
Yabe T, Morimoto K, Kikuchi S, Nishio K, Terashima I, Nakai M: The Arabidopsis chloroplastic NifU-Like protein CnfU, which can act as an iron-sulfur cluster scaffold protein, is required for biogenesis of ferredoxin and photosystem I. Plant Cell. 2004, 16: 993-1007.
Price GD, Evans JR, Voncaemmerer S, Yu JW, Badger MR: Specific reduction of chloroplast glyceraldehyde-3-phosphate dehydrogenase-activity by antisense RNA reduces CO2 assimilation via a reduction in ribulose-bisphosphate regeneration in transgenic tobacco plants. Planta. 1995, 195: 369-378.
Harrison EP, Olcer H, Lloyd JC, Long SP, Raines CA: Small decreases in SBPase cause a linear decline in the apparent RuBP regeneration rate, but do not affect Rubsico carboxylation capacity. J Exp Bot. 2001, 52: 1779-1784.
Lefebvre S, Lawson T, Fryer M, Zakhleniuk OV, Lloyd JC, Raines CA: Increased sedoheptulose-1,7-bisphosphatase activity in transgenic tobacco plants stimulates photosynthesis and growth from an early stage in development. Plant Physiol. 2005, 138: 451-460.
Wang ML: Molecular Cloning and Transformation of Sedoheptulose-1,7-Bisphosphatase in Lycopersicon esculentum, PhD thesis. China: Shandong Agricultural University; 2011.
Baud S, Wuill me S, Dubreucq B, de Almeida A, Vuagnat C, Lepiniec L, Miquel M, Rochat C: Function of plastidial pyruvate kinases in seeds of Arabidopsis thaliana. Plant J. 2007, 52: 405-419.
Frick G, Su Q, Apel K, Armstrong GA: An Arabidopsis porB porC double mutant lacking light-dependent ADPH:protochlorophyllide oxidoreductases B and C is highly chlorophyll-deficient and developmentally arrested. Plant J. 2003, 35: 141-153.
Pattanayak GK, Tripathy BC: Overexpression of protochlorophyllide oxidoreductase C regulates oxidative stress in Arabidopsis. PLoS One. 2011, 6: e26532-
Gorinova N, Nedkovska M, Atanassov A: Cytochrome P450 monooxygenases as a tool for metabolization herbicides in plants. Biotechnol Biotechnol Eq. (Special issue): 105-115.
Narusaka Y, Narusaka M, Seki M, Umezawa T, Ishida J, Nakajima M, Enju A, Shinozaki K: Crosstalk in the responses to abiotic and biotic stresses in Arabidopsis: analysis of gene expression in cytochrome P450 gene superfamily by cDNA microarray. Plant Mol Biol. 2004, 55: 327-342.
Liu X, Li F, Tang J, Wang W, Zhang F, Wang G, Chu J, Yan C, Wang T, Chu C, Li C: Activation of the jasmonic acid pathway by depletion of the hydroperoxide lyase OsHPL3 reveals crosstalk between the HPL and AOS branches of the oxylipin pathway in rice. PLoS One. 2012, 7: e50089-
Tang GY, Wei LQ, Liu ZJ, Bi YP, Shan L: Ectopic expression of peanut acyl carrier protein in tobacco alters fatty acid composition in the leaf and resistance to cold stress. Biol Plant. 2012, 56: 493-501.
Branen JK, Shintani DK, Engeseth NJ: Expression of antisense acyl carrier protein-4 reduces lipid content in Arabidopsis leaf tissue. Plant Physiol. 2003, 132: 748-756.
Ishiguro S, Watanabe Y, Ito N, Nonaka H, Takeda N, Sakai T, Kanaya H, Okada K: SHEPHERD is the Arabidopsis GRP94 responsible for the formation of functional CLAVATA proteins. EMBO J. 2002, 21: 898-908.
Chen M, Mishra S, Heckathorn SA, Frantz JM, Krause C: Proteomic analysis of Arabidopsis thaliana leaves in response to acute boron deficiency and toxicity reveals effects on photosynthesis, carbohydrate metabolism, and protein synthesis. J Plant Physiol. 2014, 171: 235-242.
Uluisik I, Alaattin Kaya A, Fomenko DE, Karakaya HC, Carlson BA, Gladyshev VN, Koc A: Boron stress activates the general amino acid control mechanism and inhibits protein synthesis. PLoS One. 2012, 6: e27772-
Mazzucotelli E, Mastrangelo AM, Crosatti C, Guerra D, Stanca AM, Cattivelli L: Abiotic stress response in plants: when post-transcriptional and post-translational regulations control transcription. Plant Sci. 2008, 174: 420-431.
García-Lorenzo M, Sjödin A, Jansson S, Funk C: Protease gene families in Populus and Arabidopsis. BMC Plant Biol. 2006, 6: 30-
Sauter M, Moffatt B, Saechao MC, Hell R, Wirtz M: Methionine salvage and S-adenosylmethionine: essential links between sulfur, ethylene and polyamine biosynthesis. Biochem J. 2013, 451: 145-154.
Roessner U, Patterson JH, Forbes MG, Fincher GB, Langridge P, Bacic A: An investigation of boron toxicity in barley using metabolomics. Plant Physiol. 2006, 142: 1087-1101.
Cervilla LM, Blasco B, Ríos JJ, Rosales MA, Rubio-Wilhelmi MM, Sánchez-Rodríguez E, Romero L, Ruiz JM: Response of nitrogen metabolism to boron toxicity in tomato plants. Plant Biol. 2009, 11: 671-677.
Baskar V, Gururani MA, Yu JW, Park SW: Engineering glucosinolates in plants: current knowledge and potential uses. Appl Biochem Biotechnol. 2012, 168: 1694-1717.
Kushwaha HR, Singh AK, Sopory SK, Singla-Pareek SL, Pareek A: Genome wide expression analysis of CBS domain containing proteins in Arabidopsis thaliana (L.) Heynh and Oryza sativa L. reveals their developmental and stress regulation. BMC Genomics. 2009, 10: 200-
Singh AK, Kumar R, Pareek A, Sopory SK, Singla-Pareek SL: Overexpression of rice CBS domain containing protein improves salinity, oxidative, and heavy metal tolerance in transgenic tobacco. Mol Biotechnol. 2012, 52: 205-216.
Zhou CP, Qi YP, You X, Yang LT, Guo P, Ye X, Zhou XX, Ke FJ, Chen LS: Leaf cDNA-AFLP analysis of two citrus species differing in manganese tolerance in response to long-term manganese-toxicity. BMC Genomics. 2013, 14: 621-
Yang LT, Qi YP, Lu YB, Guo P, Sang W, Feng H, Zhang HX, Chen LS: iTRAQ protein profile analysis of Citrus sinensis roots in response to long-term boron-deficiency. J Proteomics. 2013, 93: 179-206.
George GM, van der Merwe MJ, Nunes-Nesi A, Bauer R, Fernie AR, Kossmann J, Lloyd JR: Virus-induced gene silencing of plastidial soluble inorganic pyrophosphatase impairs essential leaf anabolic pathways and reduces drought stress tolerance in Nicotiana benthamiana. Plant Physiol. 2010, 154: 55-66.
Hernández-Domíguez EE, Valencia-Turcotte LG, Rodríguez-Sotres R: Changes in expression of soluble inorganic pyrophosphatases of Phaseolus vulgaris under phosphate starvation. Plant Sci. 2012, 187: 39-48.
Ogawa T, Ishikawa K, Harada K, Fukusaki E, Yoshimura K, Shigeoka S: Overexpression of an ADP-ribose pyrophosphatase, AtNUDX2, confers enhanced tolerance to oxidative stress in Arabidopsis plants. Plant J. 2009, 57: 289-301.
Ishikawa K, Ogawa T, Hirosue E, Nakayama Y, Harada K, Fukusaki E, Yoshimura K, Shigeoka S: Modulation of the poly (ADP-ribosyl) ation reaction via the Arabidopsis ADP-ribose/NADH pyrophosphohydrolase, AtNUDX7, is involved in the response to oxidative stress. Plant Physiol. 2009, 151: 741-754.
Fotopoulos V, Sanmartin M, Kanellis AK: Effect of ascorbate oxidase over-expression on ascorbate recycling gene expression in response to agents imposing oxidative stress. J Exp Bot. 2006, 57: 3933-3943.
Kato N, Esaka M: Expansion of transgenic tobacco protoplasts expressing pumpkin ascorbate oxidase is more rapid that that of wild type protoplasts. Planta. 2000, 210: 1018-1022.
Pignocchi C, Kiddle G, Hernández I, Foster SJ, Asensi A, Taybi T, Barnes J, Foyer CH: Ascorbate oxidase-dependent changes in the redox state of the apoplast modulate gene transcript accumulation leading to modified hormone signaling and orchestration of defense processes in tobacco. Plant Physiol. 2006, 141: 423-435.
Sanmartin M, Drogoudi PA, Lyons T, Pateraki I, Barnes J, Kanellis AK: Over-expression of ascorbate oxidase in the apoplast of transgenic tobacco results in altered ascorbate and glutathione redox states and increased sensitivity to ozone. Planta. 2003, 216: 918-928.
Garchery C, Gest N, Do PT, Alhagdow M, Baldet P, Menard G, Rothan C, Massot C, Gautier H, Aarrouf J, Fernie AR, Stevens R: A diminution in ascorbate oxidase activity affects carbon allocation and improves yield in tomato under water deficit. Plant Cell Environ. 2013, 36: 159-175.
Vieira Dos Santos C, Rey P: Plant thioredoxins are key actors in the oxidative stress response. Trends Plant Sci. 2006, 11: 329-334.
Naot D, Ben-Hayyim G, Eshdat Y, Holland D: Drought, heat and salt stress induce the expression of a citrus homologue of an atypical late-embryogenesis Lea5 gene. Plant Mol Biol. 1995, 27: 619-622.
Borrell A, Cutanda MC, Lumbreras V, Pujal J, Goday A, Culiáñez-Macià FA, Pagès M: Arabidopsis thaliana atrab28: a nuclear targeted protein related to germination and toxic cation tolerance. Plant Mol Biol. 2002, 50: 249-259.
Ahmed NU, Park JI, Jung HJ, Chung MY, Cho YG, Nou IS: Characterization of thaumatin-like gene family and identification of Pectobacterium carotovorum subsp. carotovorum inducible genes in Brassica oleracea. Plant Breed Biotech. 2013, 1: 111-121.
Matsumoto TK, Pardo JM, Takeda S, Bressan RA, Hasegawa PM: Tobacco and Arabidiopsis SLT1 mediate salt tolerance of yeast. Plant Mol Biol. 2001, 45: 489-500.
Antoine W, Stewart JM, Reyes BG D l: The rice homolog of the sodium/lithium tolerance (SLT1) gene functions as molecular chaperone in vitro. Physiol Plant. 2005, 125: 299-310.
Liu Y, Bassham DC: Autophagy: pathways for self-eating in plant cells. Annu Rev Plant Biol. 2012, 63: 215-237.
Xiong Y, Contento AL, Bassham DC: AtATG18a is required for the formation of autophagosomes during nutrient stress and senescence in Arabidopsis thaliana. Plant J. 2005, 42: 535-546.
Xiong Y, Contento AL, Nguyen PQ, Bassham DC: Degradation of oxidized proteins by autophagy during oxidative stress in Arabidopsis. Plant Physiol. 2007, 143: 291-299.
Richter CV, Bals T, Schünemann D: Component interactions, regulation and mechanisms of chloroplast signal recognition particle-dependent protein transport. Eur J Cell Biol. 2010, 89: 965-973.
Zhang F, Luo X, Hu B, Yong Wan Y, Xie J: YGL138(t), encoding a putative signal recognition particle 54 kDa protein, is involved in chloroplast development of rice. Rice. 2013, 6: 7-
Chevalier D, Morris ER, Walker JC: 14-3-3 and FHA domains mediate phosphoprotein interactions. Annu Rev Plant Biol. 2009, 60: 67-91.
Chen F, Li Q, Sun L, He Z: The rice 14-3-3 gene family and its involvement in responses to biotic and abiotic stress. DNA Res. 2006, 13: 53-63.
Campo S, Peris-Peris C, Montesinos L, Peñas G, Messeguer J, San Segundo B: Expression of the maize ZmGF14-6 gene in rice confers tolerance to drought stress while enhancing susceptibility to pathogen infection. J Exp Bot. 2012, 63: 983-999.
Yan J, He C, Wang J, Mao Z, Holaday SA, Allen RD, Zhang H: Overexpression of the Arabidopsis 14-3-3 protein GF14λ in cotton leads to a "stay-green" phenotype and improves stress tolerance under moderate drought conditions. Plant Cell Physiol. 2004, 45: 1007-1014.
Joo JH, Wang SY, Chen JG, Jones AM, Fedoroff NV: Different signaling and cell death roles of heterotrimeric G protein a and b subunits in the Arabidopsis oxidative stress response to ozone. Plant Cell. 2005, 17: 957-970.
Dodd AN, Salathia N, Hall A, Kévei E, Tóth R, Nagy F, Hibberd JM, Millar AJ, Webb AA: Plant circadian clocks increase photosynthesis, growth, survival, and competitive advantage. Science. 2005, 309: 630-633.
Nakamichi N, Kusano M, Fukushima A, Kita M, Ito S, Yamashino T, Saito K, Sakakibara H, Mizuno T: Transcript profiling of an Arabidopsis PSEUDO RESPONSE REGULATOR arrhythmic triple mutant reveals a role for the circadian clock in cold stress response. Plant Cell Physiol. 2009, 50: 447-462.
An D, Yang J, Zhang P: Transcriptome profiling of low temperature-treated cassava apical shoots showed dynamic responses of tropical plant to cold stress. BMC Genomics. 2012, 13: 64-
Byun YJ, Kim HJ, Lee DH: LongSAGE analysis of the early response to cold stress in Arabidopsis leaf. Planta. 2009, 229: 1181-1200.
Fukushima A, Kusano M, Nakamichi N, Kobayashi M, Hayashi N, Sakakibara H, Mizuno T, Saito K: Impact of clock-associated Arabidopsis pseudo-response regulators in metabolic coordination. Proc Natl Acad Sci USA. 2009, 106: 7251-7256.
Shen H, Chen J, Wang Z, Yang C, Sasaki T, Yamamoto Y, Matsumoto H, Yan X: Root plasma membrane H+-ATPase is involved in the adaptation of soybean to phosphorus starvation. J Exp Bot. 2006, 57: 1353-1362.
Wu Z, Liang F, Hong B, Young JC, Sussman MR, Harper JF, Sze H: An endoplasmic reticulum-bound Ca2+/Mn2+ pump, ECA1, supports plant growth and confers tolerance to Mn2+ stress. Plant Physiol. 2002, 130: 128-137.
Takahashi R, Bashir K, Ishimaru Y, Nishizawa NK, Nakanishi H: The role of heavy-metal ATPases, HMAs, in zinc and cadmium transport in rice. Plant Signal Behav. 2012, 7: 1605-1607.
García MNM, País SM, Téllez-Iñón MT, Capiati DA: Characterization of StPPI1, a proton pump interactor from Solanum tuberosum L. that is up-regulated during tuber development and by abiotic stress. Planta. 2011, 233: 661-674.
Kim DY, Bovet L, Maeshima M, Martinoia E, Lee Y: The ABC transporter AtPDR8 is a cadmium extrusion pump conferring heavy metal resistance. Plant J. 2007, 50: 207-218.
Ranjan A, Pandey N, Lakhwani D, Dubey NK, Pathre UV, Sawant SV: Comparative transcriptomic analysis of roots of contrasting Gossypium herbaceum genotypes revealing adaptation to drought. BMC Genomics. 2012, 13: 680-
Yuan M, Li X, Xiao J, Wang S: Molecular and functional analyses of COPT/Ctr-type copper transporter-like gene family in rice. BMC Plant Biol. 2011, 11: 69-
Sancenón V, Puig S, Mateu-Andrés I, Dorcey E, Thiele DJ, Peñarrubia L: The Arabidopsis copper transporter COPT1 functions in root elongation and pollen development. J Biol Chem. 2004, 279: 15348-15355.
Sunkar R, Kaplan B, Bouche N, Arazi T, Dolev D, Talke I, Maathuis FJM, Sanders D, Bouchez D, Fromm H: Expression of a truncated tobacco NtCBP4 channel in transgenic plants and disruption of the homologous Arabidopsis CNGC1 gene confer Pb2+ tolerance. Plant J. 2000, 24: 533-542.
Chan CW, Schorrak LM, Smith RK, Bent AF, Sussman MR: A cyclic nucleotide-gated ion channel, CNGC2, is crucial for plant development and adaptation to calcium stress. Plant Physiol. 2003, 132: 728-731.
Arazi T, Sunkar R, Kaplan B, Fromm H: A tobacco plasma membrane calmodulin-binding transporter confers Ni2+ tolerance and Pb2+ hypersensitivity in transgenic plants. Plant J. 1999, 20: 171-182.
K hler C, Merkle T, Roby D, Neuhaus G: Developmentally regulated expression of a cyclic nucleotide-gated ion channel from Arabidopsis indicates its involvement in programmed cell death. Planta. 2001, 213: 327-332.
Ohno H: Membrane traffic in multicellular systems: more than just a housekeeper. J Biochem. 2006, 139: 941-942.
Wang HJ, Wan AR, Jauh GY: An actin-binding protein, LILIM1, mediates calcium and hydrogen regulation of actin dynamics in pollen tubes. Plant Physiol. 2008, 147: 1619-1636.
Ernst AM, Jekat SB, Zielonka S, Müller B, Neumann U, Rüping B, Twyman RM, Krzyzanek V, Prüfer D, Noll GA: Sieve element occlusion (SEO) genes encode structural phloem proteins involved in wound sealing of the phloem. Proc Natl Acad Sci USA. 2012, 109: E1980-E1989.
Shen Z, Collatos AR, Bibeau JP, Furt F, Vidali L: Phylogenetic analysis of the kinesin superfamily from Physcomitrella. Front Plant Sci. 2012, 3: 230-
Li J, Jiang J, Qian Q, Xu Y, Zhang C, Xiao J, Du C, Luo W, Zou G, Chen M, Huang Y, Feng Y, Cheng Z, Yuan M, Chong K: Mutation of rice BC12/GDD1, which encodes a kinesin-like protein that binds to a GA biosynthesis gene promoter, leads to dwarfism with impaired cell elongation. Plant Cell. 2011, 23: 628-640.
Nishihama R, Soyano T, Ishikawa M, Araki S, Tanaka H, Asada T, Irie K, Ito M, Terada M, Banno H, Yamazaki Y, Machida Y: Expansion of the cell plate in plant cytokinesis requires a kinesin-like protein/MAPKKK complex. Cell. 2002, 109: 87-99.
Slewinski TL: Diverse functional roles of monosaccharide transporters and their homologs in vascular plants: a physiological perspective. Mol Plant. 2011, 4: 641-662.
Quirino BF, Reiter WD, Amasino RD: One of two tandem Arabidopsis genes homologous to monosaccharide transporters is senescence-associated. Plant Mol Biol. 1999, 46: 447-457.
Yamada K, Osakabe Y, Mizoi J, Nakashima K, Fujita Y, Shinozaki K, Yamaguchi-Shinozaki K: Functional analysis of an Arabidopsis thaliana abiotic stress-inducible facilitated diffusion transporter for monosaccharides. J Biol Chem. 2010, 285: 1138-1146.
Ricachenevsky FK, Menguer PK, Sperotto RA, Williams LE, Fett JP: Roles of plant metal tolerance proteins (MTP) in metal storage and potential use in biofortification strategies. Front Plant Sci. 2013, 4: 144-
Jiang HX, Tang N, Zheng JG, Chen LS: 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. 2009, 9: 102-
Fu C, Mielenz JR, Xiao X, Ge Y, Hamilton CY, Rodriguez M, Chen F, Foston M, Ragauskas A, Bouton J, Dixon RA, Wang ZY: Genetic manipulation of lignin reduces recalcitrance and improves ethanol production from switchgrass. Proc Natl Acad Sci USA. 2011, 108: 3803-3808.
Rastogi S, Dwivedi UN: Down-regulation of lignin biosynthesis in transgenic Leucaena leucocephala harboring O-methyltransferase gene. Biotechnol Prog. 2006, 22: 609-916.
Loomis WD, Durst RW: Chemistry and biology of boron. Biofactors. 1992, 3: 229-239.
Papuga J, Hoffmann C, Dieterle M, Moes D, Moreau F, Tholl S, Steinmetz A, Thomas C: Arabidopsis LIM proteins: a family of actin bundlers with distinct expression patterns and modes of regulation. Plant Cell. 2010, 22: 3034-3052.
Thomas C, Hoffmann C, Dieterle M, Van Troys M, Ampe C, Steinmetz A: Tobacco WLIM1 is a novel F-Actin binding protein involved in actin cytoskeleton remodeling. Plant Cell. 2006, 18: 2194-2206.
Kawaoka A, Kaothien P, Yoshida K, Endo S, Yamada K, Ebinuma H: Functional analysis of tobacco LIM protein Ntlim1 involved in lignin biosynthesis. Plant J. 2000, 22: 289-301.
Moes D, Gatti S, Hoffmann C, Dieterle M, Moreau F, Neumann K, Schumacher M, Diederich M, Grill E, Shen WH, Steinmetz A, Thomas C: A LIM domain protein from tobacco involved in actin-bundling and histone gene transcription. Mol Plant. 2013, 6: 483-502.
Sun MM, Li LH, Xie H, Ma RC, He YK: Differentially expressed genes under cold acclimation in Physcomitrella patens. J Biochem Mol Biol. 2007, 40: 986-1001.
Harholt J, Sørensen I, Fangel J, Roberts A, Willats WGT, Scheller HV, Petersen BL, Banks JA, Ulvskov P: The glycosyltransferase repertoire of the spikemoss Selaginella moellendorffii and a comparative study of its cell wall. PLoS One. 2012, 7: e35846-
Scheible WR, Pauly M: Glycosyltransferases and cell wall biosynthesis: novel players and insights. Curr Opin Plant Biol. 2004, 7: 285-295.
Bowles D, Lim EK, Poppenberger B, Vaistij FE: Glycosyltransferases of lipophilic small molecules. Annu Rev Plant Biol. 2006, 57: 567-597.
Tognetti VB, Van Aken O, Morreel K, Vandenbroucke K, van de Cotte B, De Clercq I, Chiwocha S, Fenske R, Prinsen E, Boerjan W, Genty B, Stubbs KA, Inzé D, Van Breusegem F: Perturbation of indole-3-butyric acid homeostasis by the UDP-glucosyltransferase UGT74E2 modulates Arabidopsis architecture and water stress tolerance. Plant Cell. 2010, 22: 2660-2679.
Sun YG, Wang B, Jin SH, Qu XX, Li YJ, Hou BK: Ectopic expression of Arabidopsis glycosyltransferase UGT85A5 enhances salt stress tolerance in tobacco. PLoS One. 2013, 8: e59924-
Kim IA, Heo JO, Chang KS, Lee SA, Lee MH, Lim CE, Lim J: Overexpression and inactivation of UGT73B2 modulate tolerance to oxidative stress in Arabidopsis. J Plant Biol. 2010, 53: 233-239.
Rodríguez-Villalón A, Gas E, Rodríguez-Concepción M: Phytoene synthase activity controls the biosynthesis of carotenoids and the supply of their metabolic precursors in dark-grown Arabidopsis seedlings. Plant J. 2009, 60: 424-435.
Zhang J, Tao N, Xu Q, Zhou W, Cao H, Xu J, Deng X: Functional characterization of citrus PSY gene in Hongkong kumquat (Fortunella hindsii Swingle). Plant Cell Rep. 2009, 28: 1737-1746.
Han H, Li Y, Zhou S: Overexpression of phytoene synthase gene from Salicornia europaea alters response to reactive oxygen species under salt stress in transgenic Arabidopsis. Biotechnol Lett. 2008, 30: 1501-1507.
Cidade LC, de Oliveira TM, Mendes AFS, Macedo AF, Floh EIS, Gesteira AS, Soares-Filho WS, Costa MGC: Ectopic expression of a fruit phytoene synthase from Citrus paradisi Macf. promotes abiotic stress tolerance in transgenic tobacco. Mol Biol Rep. 2012, 39: 10201-10209.
Dutta A, Sen J, Deswal R: New evidences about strictosidine synthase (Str) regulation by salinity, cold stress and nitric oxide in Catharanthus roseus. J Plant Biochem Biotech. 2013, 22: 124-131.
Gemici M, Aktaş LY, Türkyilmaz B, Güven A: The effects of the excessive boron applications on indole-3-acetic acid levels in Triticum durum Desf cv. Gediz seedlings. Cumhuriyet Üniversitesi Fen Bilimleri Dergisi. 2002, 23 (2): 17-24.
Pasquali G, Goddijn OJ, de Waal A, Verpoorte R, Schilperoort RA, Hoge JH, Memelink J: Coordinated regulation of two indole alkaloid biosynthetic genes from Catharanthus roseus by auxin and elicitors. Plant Mol Biol. 1992, 18: 1121-1131.
Menke FL, Champion A, Kijne JW, Memelink J: A novel jasmonate- and elicitor-responsive element in the periwinkle secondary metabolite biosynthetic gene Str interacts with a jasmonate- and elicitor-inducible AP2-domain transcription factor, ORCA2. EMBO J. 1999, 18: 4455-4463.
Wang J, Nakazato T, Kinya Sakanishi K, Yamada O, Tao H, Saito I: Single-step microwave digestion with HNO3 alone for determination of trace elements in coal by ICP spectrometry. Talanta. 2006, 68: 1584-1590.
Ames BN: Assay of inorganic phosphate, total phosphate and phosphatase. Methods Enzymol. 1966, 8: 115-118.
Bradford MM: A rapid and sensitive method for quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976, 72: 248-254.
Hodges DM, DeLong JM, Forney CF, Prange RK: Improving the thiobarbituric acid-reactive-substances assay for estimating lipid peroxidation in plant tissues containing anthocyanin and other interfering compounds. Planta. 1999, 207: 604-611.
Lichtenthaler HK: Chlorophylls and carotenoids: pigments of photosynthetic biomembranes. Methods Enzymol. 1987, 148: 350-382.
This study was financially supported by the earmarked fund for China Agriculture Research System.
The authors declare that they have no competing interests.
PG carried out most of the experiments and drafted the manuscript. YPQ participated in the design of the study. LTY participated in the design of the study and coordination. XY carried out the measurement of B and P. HXJ performed the statistical analysis. JHH carried out the cultivation of seedlings. LSC designed and directed the study and revised the manuscript. All authors have read and approved the final manuscript.