Genome-wide transcriptome and functional analysis of two contrasting genotypes reveals key genes for cadmium tolerance in barley
© Cao et al.; licensee BioMed Central Ltd. 2014
Received: 25 April 2014
Accepted: 14 July 2014
Published: 19 July 2014
Cadmium (Cd) is a severe detrimental environmental pollutant. To adapt to Cd-induced deleterious effects, plants have evolved sophisticated defence mechanisms. In this study, a genome-wide transcriptome analysis was performed to identify the mechanisms of Cd tolerance using two barley genotypes with distinct Cd tolerance.
Microarray expression profiling revealed that 91 genes were up-regulated by Cd in Cd-tolerant genotype Weisuobuzhi and simultaneously down-regulated or non-changed in Cd-sensitive Dong17, and 692 genes showed no change in Weisuobuzhi but down-regulated in Dong17. Novel genes that may play significant roles in Cd tolerance were mainly via generating protectants such as catalase against reactive oxygen species, Cd compartmentalization (e.g. phytochelatin-synthase and vacuolar ATPase), and defence response and DNA replication (e.g. chitinase and histones). Other 156 up-regulated genes in both genotypes also included those encoding proteins related to stress and defence responses, and metabolism-related genes involved in detoxification pathways. Meanwhile, biochemical and physiological analysis of enzyme (ATPase and chitinase), phytohormone (ethylene), ion distribution and transport (Cd, Na+, K+, Ca2+, ABC transporter) demonstrated that significantly larger Cd-induced increases of those components in Weisuobuzhi than those in Dong17. In addition, Cd-induced DNA damage was more pronounced in Dong17 than that in Weisuobuzhi.
Our findings suggest that combining microarray, physiological and biochemical analysis has provided valuable insights towards a novel integrated molecular mechanism of Cd tolerance in barley. The higher expression genes in Cd tolerant genotype could be used for transgenic overexpression in sensitive genotypes of barley or other cereal crops for elevating tolerance to Cd stress.
KeywordsATPase Cadmium fluorescent indicator Cadmium tolerance Microarray Hordeum vulgare L Integrated molecular mechanism
Cadmium (Cd) in soil represents a direct contact risk to both humans and ecological recipients due to its high toxicity and ready uptake by plants [1, 2]. At present, Cd has become one of the most harmful and widespread pollutants in agricultural soils. This situation has resulted primarily from industrial emissions, application of Cd containing phosphate fertilisers and municipal waste disposal [1, 3]. Once released into the soils, moderate Cd pollution could result in considerable Cd accumulation in edible parts of crops. Such levels of Cd not only affect the quality and yield of crops, but pose a great threat to human health . Accordingly, these developments raise serious concerns for both the environment and human health. Therefore, there is an urgent need to elucidate the mechanisms of Cd tolerance in plants and to develop crop varieties with high Cd tolerance and yield.
Cd affects many important physiological processes and inhibits plant growth and development [4, 5]. To minimise the detrimental effects of Cd toxicity, plants have evolved a range of detoxification mechanisms, including Cd exclusion, chelation and compartmentalisation in vacuoles . For example, phytochelatins (PCs) and Cd have been found to form PC-Cd complexes in cytosol, which are subsequently transported into vacuoles. Thereby, it can protect plants from the deleterious effects of Cd . The elevated expression of heavy metal transporting ATPases4 (HMA4) P1B-type ATPase furnishes an efficient mechanism for increasing Cd tolerance in plants under Cd toxicity via the maintenance of low cellular Cd in the cytoplasm . However, regulatory mechanisms in Cd-tolerance, which is still a focal point in plant research, is a complex process that contains many genes regulated by a variety of physiological pathways .
Identification of stress-induced genes and proteins are fundamental in understanding the molecular mechanisms of stress responses and in developing transgenic plants with enhanced tolerance . Studies have been conducted to identify plant defences against Cd toxicity and investigate Cd-specific genes. For example, the up-regulation of well-known Cd-detoxifying proteins, such as phytochelatin synthase (PCS), antioxidative enzymes and glutathione S-transferases (GST), were observed in the response of plants to Cd stress using proteomic and metabolomic approaches . Uraguchi and Fujiwara  summarised several Cd transport-related genes, such as OsLCT1, OsHMA, OsNramp1 and OsIRT1, involved in Cd transport and tolerance in rice. Despite the identification of those genes, the underlying knowledge of molecular mechanisms for plant Cd tolerance is still fragmental.
Knowledge about networks of global gene expression is imperative for further understanding the molecular mechanisms in plant Cd tolerance. Microarray analysis is a powerful technique for analysing the profiles of gene expression related to abiotic stress in plants [11, 12]. Concerning Cd toxicity and detoxification in plants, the genome-wide transcriptome profiling has been explored mostly in herbaceous plants such as Arabidopsis thaliana, Arabidopsis halleri, and Thlaspi caerulescen[13–16]. However, the regulatory system for genes conferring Cd tolerance in many cereal crops is largely unknown and remains an essential issue to be addressed.
Barley, the fourth most important cereal in the world, is as an ideal model for genetic and physiological study . Our previous studies have demonstrated distinctive genotypic difference between Cd-tolerant (Weisuobuzhi) and Cd-sensitive (Dong17) genotypes in response to Cd stress including plant growth, photosynthesis, antioxidant enzyme activities [1, 18, 19]. We, therefore, hypothesised that there were large differences between the two genotypes in their genome-wide response to Cd stress. Candidate genes associated with Cd tolerance in the tolerant genotype were identified using the Affymertix barley GeneChip whole genome arrays, and the expression and enzyme activity of many key genes were further validated by quantitative real-time PCR (qRT-PCR) and physiological and biochemical analysis. We also proposed an integrated schematic diagram of the mechanisms involved in Cd tolerance and adaptation. These results are very important to provide novel clues for understanding the mechanisms in Cd tolerance of barley and open prospective for excavating novel genes and for the genetic improvement of plant tolerance to Cd stress.
Plant material, growth conditions and Cd treatments
A greenhouse hydroponic experiment was carried out at Huajiachi Campus of Zhejiang University, Hangzhou, China. Two barley genotypes, Cd-tolerant Weisuobuzhi and Cd-sensitive Dong17 , were used throughout the experiments.
Barley seeds were surface sterilised by soaking in 2% H2O2 for 30 min and fully rinsed with deionised water. After soaking in deionised water at room temperature for 4 h, seeds were germinated in sterilised moist sand in an incubator at 22 ± 1°C. Ten-day-old healthy and uniform plants were selected and transplanted to 5-L containers filled with 4.5 L basal nutrient solution. The composition of nutrient solution was described by Chen et al. [18, 19]. A week after transplanting, Cd (as CdCl2) was added to the corresponding containers to form 4 treatments: basal nutrient solution (control, without Cd) and 5, 50 and 500 μM Cd. The experiment was laid in a split-plot design with Cd-treatment as main-plot with six replicates. The nutrient solution was continuously aerated with pumps and renewed every 5 d.
Plant samples were harvested 15 d after Cd treatment, and root tips were used for fluorescence imaging of Cd. The second fully expanded leaves of control and 5 μM Cd treatment were collected for microarray, physiological and biochemical analysis.
Fresh root tips were immersed in 20 mM disodium ethylenediamine tetra-acetic acid (Na2-EDTA) for 15 min and then gently rinsed for three times with deionised water. The specimen sections were then immediately immersed in the Cd probe Leadmium™ Green AM solution (Molecular Probes, Life Technologies, California, USA) for 45 min in the dark and then washed three times (5 min each time) with deionised water. A stock solution of Leadmium™ Green AM was made by adding 50 μL of dimethyl sulfoxide (DMSO) to one vial of the dye. This stock solution was then diluted 1:10 with 0.85% NaCl. Roots were immersed in this solution for 90 min in the dark. The sections were examined with a laser confocal scanning microscope (Leica TCS SP5; Berlin, Germany) with excitation and emission wavelengths at 488 and 515 nm, respectively. The fluorescence density of Cd was calculated by selecting the root sections and measuring the total Integrated Density using “Analyse and Measure” function of the Image J software (NIH, Bethesda, MD, USA).
RNA isolation, probe preparation and microarray hybridisation
Total RNA was initially extracted from leaf tissue of Weisuobuzhi and Dong17 under control and 5 μM Cd treatment using Trizol reagent protocol (P/N 15596–018, Life Technologies, Carlsbad, CA, USA). RNA was purified on RNeasy spin columns (RNeasy MinElute Cleanup Kit, Qiagen, Dusseldorf, Germany) and with on-column DNase1 treatment. The eluted total RNAs were quantified on an Agilent 2100 bioanalyser (Palo Alto, CA, USA) and adjusted to a final concentration of 1 μg μL-1. Sample processing, cDNA synthesis, biotin-labelled cRNA synthesis, hybridisation, washing, staining and scanning of Affymetrix Barley 1 GeneChip were performed following the standard protocol (Affymetrix Inc., Santa Clara, CA, USA). The Barley 1.0 Affymetrix microarray GeneChip array consists of 22,795 probe sets designed from an exemplar set of barley sequences derived from 350,000 high-quality expressed sequence tag (EST) contigs. Total RNA was used to prepare double strand cDNA. Labelled cRNA preparation and hybridisation on GeneChip and scanning was done following the standard Affymetrix procedures (http://media.affymetrix.com/support/downloads/manuals/3_ivt_express_kit_manual.pdf). Three replications of each sample were conducted to test the quality and reproducibility of the chip hybridisation. Each treatment had three replications.
Microarray data analysis
Data analysis was conducted using Refiner and Analyst, two analytical tools in Affymetrix GeneChip Operating Software Version 1.4. The Refiner tool condenses and normalises the raw signal. The Analyst tool provides statistical analysis and data visualisation capability. Detection, signal condensation and normalisation were conducted using Robust Multichip Analysis (RMA) of the Refiner tool . Further analysis of transcript abundance was conducted using the Analyst tool. The correlation of expression signals between replications for each genotype was 0–0.99 across all probe sets tested on the GeneChip. Each Barley1.0 GeneChip probe set was first tested for barley transcript specificity on the basis of quantitative and qualitative transcript abundance differences between barley genotype Betzes and wheat cultivar Chinese Spring (CS).
To detect barley transcripts in the two lines (Weisuobuzhi and Dong17), differences in the abundance of barley transcript signals in the two lines were tested against the background signal in CS with a t-test (P-value < 0.001). Only those showing significantly higher transcript levels both in Betzes and in the two lines compared to those in CS were subjected to qualitative analysis to detect barley-specific transcripts in the two lines (presence/absence test, P < 0.001 using MAS 5.0). All data from the Affymetrix scanner have been deposited at BarleyBase (http://www.barleybase.org) in the form of DAT, CEL, EXP, and CHP files. We considered a change of at least two folds as an indication of a significant change in gene expression for up- and down-regulation. To perform hierarchical clustering, the differentially regulated genes were clustered according to similarities in expression profiles. The following nine groups were used: Weisuobuzhi (W) up-regulated and Dong17 (D) down-regulated (Group I) or no change (Group II), W no change and D down-regulated (Group III), both up-regulated (Group IV), W down-regulated and D up-regulated (Group V), both down-regulated (Group VI), W down-regulated and D no change (Group VII), W no change and D up-regulated (Group VIII), and both no change (Group IX).
The hydroponic experiment was carried once again using Weisuobuzhi and Dong17 under control and 5 μM Cd treatment with four replicates. Total RNA was isolated from leaves after 15 days of Cd treatment using the TRIzol reagent (Invitrogen, Karlsruhe, Germany). Residual DNA was removed using purifying columns. One microgram of each RNA sample was subsequently employed for cDNA synthesis with 0.5 μg of oligo (dT) 12–18 and 200 units of Superscript II (Invitrogen, Karlsruhe, Germany). cDNA samples used for GeneChip analysis and this experiment were assayed by quantitative real time PCR (qRT-PCR) in an iCycler iQ™ Real-time PCR Detection System (Bio-Rad, Hercules, CA, USA) using the SYBR Green PCR Master Mix (Applied Biosystems, Life Technologies, CA, USA). The PCR conditions consisted of denaturation at 95°C for 3 min, followed by 40 cycles of denaturation at 94°C for 1 min, annealing at 58°C for 30 s and extension at 72°C for 30 s, and continued extension at 72°C for 5 min. The barley ACTIN (accession no. AY145451) was used as a reference gene, which were selected from a number of candidates by Chen et al. [18, 19]. The primers are listed in Additional file 1: Table S1. Two independent experiments and six biological replications in total were conducted.
Enzyme activity measurements
For the determination of GST (EC 22.214.171.124) and ATPase (EC 126.96.36.199) activities, plant tissue was homogenised in 8 ml 50 mM sodium phosphate buffer (PBS, pH 7.8,) using a pre-chilled mortar and pestle, then centrifuged at 10000 × g for 20 min at 4°C. The supernatant was used for enzyme assay. GST and ATPase activities were determined with an enzyme assay kit according to the manufacturers’ protocol (Jiancheng Bio Co., Nanjing, China). Chitinase (EC 188.8.131.52) were extracted and determined with an ELISA kit according to the manufacturers’ protocol (Becton, Dickinson and Company, FranklinLakes, NJ, USA).
Determination of DNA damage
DNA damage was assayed according to Menke et al.  with minor modification. The experiment was performed in a darkroom with dim red light. About 150 mg leaf samples were sliced in 1 mL PBS buffer (160 mM NaCl, 8 mM Na2HPO4, 4 mM NaH2PO4, pH 7) containing 50 mM EDTA on ice with a fresh razor blade. The 30 μL suspension was then taken on regularly used microscopic slides (pre-coated with 1% of agarose in double distilled H2O and dried over night at room temperature) followed by addition of 30 μL 1% agarose solution (42°C). DNA damage was analysed by the comet assay according to the alkali-alkali (A/A) method as described by Menke et al. . For comet assay, unwinding was done in high alkali for 5 min, and then electrophoresis for 10 min with 21 V, 300 mA in a chamber cooled on ice, followed by a short neutralisation of 3 min in 100 mM Tris–HCl (pH 7.5). To remove the starch grains, the slides were kept for 10 min in 1% Triton prior to dehydration in 70% (2 min) and 96% (5 min) ethanol and air-drying. The gels were then stained with 15 μL ethidium bromide (5 μg ml-1) and immediately used for evaluation. Images were taken using a fluorescence microscope (BX50WI, Olympus, Japan) equipped with a digital CCD camera (Olympus, Japan).
Ethylene production in the leaves of barley seedling was measured as described by Chen et al.  with minor modification. In brief, leaf tissues (0.5 g fresh weight) were collected, immediately weighed, and then placed in sealed glass vials containing water-saturated filter paper. Samples were incubated for 1 h under illumination. One millilitre of gas was collected using a gas-tight syringe and injected into a gas chromatography (Focus GC, Thermo Fisher Scientific, MA, USA) equipped with a capillary column (CP-carboPLOT P7, Varian, CA, USA) and flame-ionisation detector for ethylene determination. Ethylene production was calculated on the basis of known ethylene standards and leaf fresh weight.
Leaf element analysis
Leaf samples were dried at 80°C, digested with HNO3/HClO4 (4:1, v/v) for 3 h, and then diluted to 25 ml by adding de-ionized water. Calcium (Ca2+), sodium (Na+), and potassium (K+) concentrations were determined using an Inductively Coupled Plasma Optical Emission Spectrometer (ICP-OES, Optima 8000DV, PerkinElmer, USA).
All data are the averages of three to six biological replicates in each experiment. Statistical analyses were performed with Data Processing System (DPS) statistical software. An ANOVA followed by a Duncan’s multiple range test (DMRT) were used to evaluate treatment effects at significance level of P < 0.05.
Cd is highly accumulated in root tips of a Cd-tolerant genotype
Cd toxicity induces large scale changes in gene expression
The distinct cellular Cd accumulation pattern (Figure 1) and tolerance [18, 19] of the two genotypes has led us to explore the expression pattern of their Cd-induced genes. The microarray data showed that, compared with control, the gene expression profiles of the two genotypes changed significantly after exposing to 5 μM Cd for 15 d. Based on at least ±2.0-fold changes (P < 0.05), a total of 1,750 genes were differentially expressed between Cd-stressed and control plants. Of these genes, only 247 and 103 genes were up- and down-regulated in Weisuobuzhi, while there were 3.6- (898) and 6.9-fold (710) more genes up-regulated and down-regulated in Dong17, respectively (Additional file 2: Figure S1 and Additional file 3: Table S2).
Key genes for Cd tolerance are highly expressed only in the Cd-tolerant genotype
These genes were up-regulated in Cd-tolerant genotype Weisuobuzhi (W-up) and down- (D-down) or un-regulated (D-none) in Cd sensitive genotype Dong17 (Additional file 4: Table S3 to Additional file 5: Table S4, Figure 3A). Group I included 7 genes encoding key enzymes for nitrogen catabolism such as asparaginase and C13 endopeptidase NP1 precursor. The 84 genes in Group II were up-regulated after Cd treatment in Weisuobuzhi but showed no change in Dong17 (Additional file 5: Table S4). Of these, there were 6 signal transduction related genes (e.g. the gene encoding 23 KD jasmonate-induced protein 1), 11 stress and defence response related genes (such as the genes encoding catalase isozyme 2, peroxidase, and chitinase 2), 10 carbohydrate and fat metabolism related genes (e.g. the gene encoding acetyl-CoA synthetase). Additionally, the expression of 10 transporter genes (e.g. the genes encoding vacuolar ATP synthase, 16 kDa proteolipid subunit and putative lipid transfer protein) under Cd stress was up-regulated in Weisuobuzhi but showed no change in Dong17.
Key genes for Cd-tolerance are down-regulated only in the Cd-sensitive genotype
The majority of the 692 genes in Group III were responsible for stress and defence response, transport, transcription and signal transduction (Additional file 6: Table S5, Figure 3B). For instance, the salicylic acid (SA), jasmonate (JA), ethylene (ET) and Ca2+ induced genes were all significantly down-regulated in Cd-sensitive genotype Dong17, but were not affected in Weisuobuzhi after 15 d of Cd exposure, with the exception of a 23 kDa JA-induced protein that was up-regulated in Weisuobuzhi (Additional file 5: Table S4 and Additional file 6: Table S5). There were 2 isoforms of GST showing differential expression in both genotypes (Additional file 6: Table S5). Genes encoding nitrate (NRaT) and nitrite (NRiT) transporters, as well as a gene encoding dehydrins (late embryogenesis abundant, LEA) were down-regulated in Dong17 but remained unchanged in Weisuobuzhi. Heat shock proteins (HSPs) were highly down-regulated under Cd stress in Dong17.
qRT-PCR confirms expression pattern of the microarray
Biochemical and physiological validation of the roles of key genes in Cd tolerance
Effects of 15 d of Cd treatment on ethylene emission and the activities of glutathione S-transferases (GST) and chitinase in leaves of two barley genotypes
Ethylene emission (nL g-1FW h-1)
GST (unit g-1FW)
Chitinase (unit g-1FW)
5 μM Cd
5 μM Cd
Effects of 15 d of Cd treatment on Ca 2+ , Na + , and K + content in leaves of two barley genotypes
5 μM Cd
5 μM Cd
In comparison to the control, Cd stress resulted in significant increases in the activity of H+-, Na+K+-, Ca2+Mg2+- and total-ATPase of 217%, 140%, 104%, and 143% in Weisuobuzhi, but those numbers were only 80%, 74%, 69% and 72% for Dong17, respectively (Figure 5). Moreover, Cd stress caused a 110% increase in GST activity in Weisuobuzhi but no statistically significant change was observed in Dong17 (Table 1). The significant increase of chitinase activity under Cd stress was 1.7-fold higher in Weisuobuzhi than that in Dong17 (Table 1).
DNA appeared as densely condensed structured resembling a bead and no DNA damage was observed in both genotypes in the control (Figure 6). However, Cd stress induced a marked increase in DNA damage in the Cd-sensitive genotype Dong17. Moreover, Cd stress led to a significant increase in ethylene emission in both genotypes. However, the rate of ethylene emission was significant higher in Weisuobuzhi than that in Dong17 (Table 1).
Cd treatment resulted in a significant decrease of leaf Ca2+ concentration, which was 2.0-fold higher for Dong17 in contrast to Weisuobuzhi. Cd stress caused an 11.0% decrease of K+ content in Dong17, but no significant difference was observed in Weisuobuzhi. Surprisingly, Na+ content was also markedly reduced under Cd stress in both genotypes (Table 2).
Comparative transcriptome analysis reveals key genes associated with Cd tolerance
Detoxification genes are largely responsible for Cd tolerance in barley
One aspect of Cd tolerance and adaptation is the ability to scavenge Cd-induced reactive oxygen species (ROS) in order to protect membranes and other macromolecules . Antioxidant enzymes, such as catalase (CAT) and CuZn-superoxide dismutase (SOD), play key roles in scavenging ROS under heavy metal stress . GST includes a family of potent detoxification enzymes, and is thought to contribute to the primary cellular defence against oxidative stress . In this study, the expression levels of the genes encoding CAT2 and CuZn-SOD confirmed the higher enzyme activity in Weisuobuzhi in our previous publication . Also, higher transcripts of GST found in both microarray and qRT-PCR were matched by a significantly elevated GST activity in Weisuobuzhi as compared to a small decrease in Dong17 (Table 1; Figure 4; Additional file 5: Table S4 and Additional file 6: Table S5). As a result, Dong17 accumulated more O2 . , H2O2 and MDA than Weisuobuzhi .
Membrane bound enzyme ATPases provide energy and an H+ gradient for the co-transport of copper, cobalt, lead, and Cd with protons for the detoxification of these metal ions in plants [25, 26]. Cation antiporter activity driven by vacuolar ATPase-dependent proton motive force contributes significantly to the detoxification of Cd via vacuolar compartmentalisation in plants . In this study, transcripts of V-ATPase were significantly elevated in Weisuobuzhi but showed no change in Dong17 (Additional file 5: Table S4). In addition, we validated the microarray data using biochemical analysis with significantly higher Cd-induced H+-, and Ca2+Mg2+-ATPase and much larger Cd-induced increase in Na+K+- and total-ATPase activities in Weisuobuzhi (Figure 5).
DNA methylation is an important modification of DNA that plays a key role in gene regulation, DNA replication and repair, and chromatin determination . In the present study, Cd significantly inhibited DNA methyltransferase in Dong17 but showed no change in Weisuobuzhi, which represented a 2.6-fold higher relative expression of this gene in Weisuobuzhi (Additional file 6: Table S5). In addition, Cd stress induced a marked increase in DNA damage to the Cd-sensitive genotype Dong17 (Figure 6). Phytochelatins are widely accepted as a major agent of plant detoxification and tolerance to Cd stress  and a lack of PC synthase activity resulted in an increased sensitivity to Cd . Our previous research has shown that genotypic differences in Cd tolerance were positively linked to the elevation of PCs in rice . In this study, transcripts of PC synthase were significantly inhibited in Dong17, while no change was found in Weisuobuzhi (Additional file 6: Table S5). In summary, the results demonstrated that Cd-tolerant Weisuobuzhi is more capable to scavenge Cd-induced ROS by increasing the activity of antioxidant enzymes, to maintain DNA structural stability and protect the normal DNA methylation, and to sequester more Cd2+ to vacuole to reduce Cd toxicity.
Membrane transport genes modulate ion homeostasis for Cd tolerance
Cd affects the distribution of nitrogen, but nitrogen could be recycled and be translocated as a Cd protection and storage strategy . Li et al.  reported that nitrate transporter NRT1.8-regulated nitrate distribution plays an important role in Cd tolerance. In addition, inwardly rectifying K+ channels (IRK) contributes to cellular K+ homeostasis in higher plants . In the present study, the genes encoding transporters for nitrogen (NRaT and NRiT) and potassium (IRK) were significantly decreased in Dong17 but remained unchanged in Weisuobuzhi (Additional file 6: Table S5). Cd stress also induced a significant decrease of K+ content in Dong17 but no change in Weisuobuzhi (Table 2), indicating a role for IRK in the reduced K+ uptake in Dong17. Moreover, ABC transporters are involved in the homeostasis of organic anions, heavy metals, xenobiotics and lipids , including vacuolar compartmentalisation of Cd . The ABC transporters of Arabidopsis, AtMRP3 and AtATM3, have been shown to confer Cd resistance [35, 36]. AtPDR8 is a Cd extrusion pump conferring Cd and Pb resistance . However, the genes encoding ABC transporters, in this study, were significantly inhibited in Dong17 while not affected in Weisuobuzhi (Additional file 6: Table S5). Therefore, a higher Cd accumulation in Weisuobuzhi (Figure 1) did not affect its overall Cd tolerance, but rather supported the hypothesis that more Cd is transported into the vacuoles of Weisuobuzhi, alleviating Cd toxicity to the cytoplasm.
As a result, these data enable us to generate a specific model on key components of vacuolar compartmentalisation for Cd tolerance (Additional file 13: Figure S3). Interestingly, ten key genes found to be down-regulated in the Cd-sensitive Dong17, were unchanged in Weisuobuzhi. All these highlighted that the high Cd accumulation in Weisuobuzhi (Figure 1; Chen et al. ) is linked with expression patterns of those specific genes involved in vacuolar compartmentalization, which is one of the most crucial strategies for Cd tolerance in plants.
Hormonal and Ca2+signal transduction related genes are important for Cd tolerance
Phytohormones such as SA, JA and ET play fundamental roles via signalling crosstalk in abiotic stresses in plants . For instance, Cd stress triggered an accumulation of ET in bean . It has also been suggested that ET signalling pathways affect the early phase of Cd stress response, and 1-aminocyclopropane-1-carboxylic acid oxidase (ACO) catalyses the last step of ET biosynthesis in Arabidopsis . Transcripts of ACO was markedly increased in Weisuobuzhi but was no change in Dong17 (Additional file 5: Table S4). ET activates stress-responsive genes in the hormonal signalling cascade against Cd toxicity. However, transcripts of ET-induced proteins were all significantly decreased in Dong17 (Additional file 6: Table S5). Cd-induced ET emission was also lower in Dong17 compared to that in Weisuobuzhi, which was consistent with the transcripts of ACO (Table 1; Additional file 5: Table S4). These results may suggest that Weisuobuzhi is more likely to defend itself against Cd stress through ET signalling. Meanwhile, Ca2+-binding protein calmodulin (CAM) transduces second messenger signals into a wide range of cellular responses . Accumulation of Cd2+ may compete with cellular Ca2+ for CAM binding sites . Under Cd stress, there was significantly less leaf Ca2+ in Dong17 than that in Weisuobuzhi (Table 2). Our results showed that Cd induces more signal molecules and activate Cd responsive genes more rapidly in Cd-tolerant Weisuobuzhi than those in Dong17. This difference again contributes to the high Cd-tolerance of Weisuobuzhi.
Carbohydrate metabolism related genes regulate cell wall structure for Cd tolerance
Cell wall, consisting of cellulose, hemicelluloses and pectin, which contain carboxyl, hydroxyl and aldehyde, can sequester a substantial amount of Cd under Cd stress . For instance, Cosio et al.  found Cd binding in leaf cell wall could play a major role in Cd tolerance and hyperaccumulation in Thlaspi caerulescens and Arabidopsis halleri. Xiong et al.  found that 200 μM Cd markedly decreased the pectin and hemicellulose content in rice root cell wall. Glycosyl hydrolases (GH) are also essential for the modification of cell wall polysaccharides . In the present study, GH and GT were found to be down-regulated in Dong17, but showed no change in Weisuobuzhi. Hence, Weisuobuzhi could maintain normal cell wall synthesis, remodelling and modification, and consequently have the capacity to accumulate much higher Cd at the cell wall (Figure 1), therefore reducing Cd toxicity to the cytoplasm.
Defence and DNA replication related genes are critical for Cd tolerance
Chitinases are components of plant defences against high concentrations of heavy metals such as Cd and As . Chitinase 2 genes showed significant up-regulation under Cd stress in Weisuobuzhi (Additional file 5: Table S4), which was confirmed by the elevated chitinase activity (Table 1). Moreover, Cd affects cell cycle progression, differentiation, DNA replication and repair . Histones such as H2B may be able to repair heavy metal-induced DNA damage in plant cells [46, 47]. In the current study, transcripts of H2A, H2B, H3 and H4 were all significantly down-regulated in Dong17 but were not modified in Weisuobuzhi (Additional file 6: Table S5). Also, genes encoding DNA replication related proteins (e.g. DNA binding protein, polymerase and origin recognition complex subunit 4) were all significantly inhibited in Dong17 under Cd stress, while no change was found in Weisuobuzhi (Additional file 6: Table S5). Considering the significant Cd-induced DNA damage found in Dong 17 (Figure 6), we thus suggest that defence and DNA replication related proteins like Chitinases and histones could be key determinants for Cd tolerance in barley.
Cd-responsive miRNAs showed a negative link to the expression of key genes in this study
miRNAs are a large family of small non-coding RNAs that negatively regulate mRNA at the post-transcriptional level . Establishing a link between published miRNAs in the literature and our microarray data could provide a better understanding and validation to further investigate the key genes in Cd tolerance in plants. We identified 13 genes, showed significant negative correlation (r2 = 0.441; P < 0.001) between the microarray data (Additional file 14: Table S11) and different miRNA families [48–51]. The results demonstrated combination of microarray and miRNA analysis could narrow down the number of candidate genes conferring Cd tolerance from thousands to tens, providing a promising outlook for future functional analysis of these genes.
The use of genome-wide transcriptome analysis highlights novel integrated molecular mechanisms associated with Cd-tolerance. Our results are potentially important for the characterisation of molecular mechanisms underlying Cd tolerance in barley. We demonstrated that Cd-tolerant Weisuobuzhi (1) is more capable to scavenge Cd-induced ROS; (2) is able to maintain ion homeostasis and sequester more Cd into the vacuoles via the ABC transporters and ATPase; (3) has higher efficiency in ET and Ca2+ signal transduction; (4) maintains normal cell wall function; and (5) expresses defence and DNA replication related proteins for Cd tolerance in barley. These distinct differences in gene expression profiles, biochemical and physiological functions between Cd-tolerant and -sensitive genotypes will provide critical information for extending our knowledge and guiding our future investigations into the candidate genes and proteins underlying Cd tolerance in barley.
This project was funded by National Natural Science Foundation of China (30571097, 31071365), the National 863 program (2012AA101105), the Ph.D. Programs Foundation of Ministry of Education of China (20110101110088), and the Key Research Foundation of Science and Technology Department of Zhejiang Province of China (2012C12902-2). Z.H.C. is supported by a Research Lectureship of University of Western Sydney and a Discovery Early Career Researcher Award (DE140101143) of the Australian Research Council. We deeply appreciate Ms Mei Li from 985-Institute of Agrobiology and Environmental Science (985-IAES) of Zhejiang University, for her kind help with the experiment.
- Chen F, Wang F, Zhang GP, Wu FB: Identification of barley varieties tolerant to cadmium toxicity. Biol Trace Elem Res. 2008, 121: 171-179.PubMedView ArticleGoogle Scholar
- Uraguchi S, Fujiwara T: Cadmium transport and tolerance in rice: perspectives for reducing grain cadmium accumulation. Rice. 2012, 5: 5-PubMed CentralPubMedView ArticleGoogle Scholar
- Wu FB, Zhang GP, Dominy P: Four barley genotypes respond differently to cadmium: lipid peroxidation and activities of antioxidant capacity. Environ Exp Bot. 2003, 50: 67-78.View ArticleGoogle Scholar
- Clijsters H, Vanassche FV: Inhibition of photosynthesis by metals. Photosynth Res. 1985, 7: 31-40.PubMedView ArticleGoogle Scholar
- Wu FB, Dong J, Qian QQ, Zhang G: Subcellular distribution and chemical form of Cd and Cd–Zn interaction in different barley genotypes. Chemosphere. 2005, 60: 1437-1446.PubMedView ArticleGoogle Scholar
- Hall JL: Cellular mechanisms for heavy metal detoxification and tolerance. J Exp Bot. 2002, 53: 1-11.PubMedView ArticleGoogle Scholar
- Cobbett CS: Phytochelatins and their roles in heavy metal detoxification. Plant Physiol. 2000, 123: 825-832.PubMed CentralPubMedView ArticleGoogle Scholar
- Courbot M, Willems G, Motte P, Arvidsson S, Roosens N, Saumitou-Laprade P, Verbruggen N: A major quantitative trait locus for cadmium tolerance in Arabidopsis halleri colocalizes with HMA4, a gene encoding a heavy metal ATPase. Plant Physiol. 2007, 144: 1052-1065.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhao CR, Ikka T, Sawaki Y, Kobayashi Y, Suzuki Y, Hibino T, Sato S, Sakurai N, Shibata D, Koyama H: Comparative transcriptomic characterization of aluminum, sodium chloride, cadmium and copper rhizotoxicities in Arabidopsis thaliana. BMC Plant Biol. 2009, 9: 32-PubMed CentralPubMedView ArticleGoogle Scholar
- Villiers F, Ducruix C, Hugouvieux V, Jarno N, Ezan E, Garin J, Junot C, Bourguignon J: Investigating the plant response to cadmium exposure by proteomic and metabolomic approaches. Proteomics. 2011, 11: 1650-1663.PubMedView ArticleGoogle Scholar
- Rabbani MA, Maruyama K, Abe H, Khan MA, Katsura K, Ito Y, Yoshiwara K, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: Monitoring expression profiles of rice genes under cold, drought, and high-salinity stresses and abscisic acid application using cDNA microarray and RNA gel-blot analyses. Plant Physiol. 2003, 133: 1755-1767.PubMed CentralPubMedView ArticleGoogle Scholar
- Sasidharan R, Mustroph A, Boonman A, Akman M, Ammerlaan AMH, Breit T, Schranz ME, Voesenek LACJ, van Tienderen PH: Root transcript profiling of two Rorippa species reveals gene clusters associated with extreme submergence tolerance. Plant Physiol. 2013, 163: 1277-1292.PubMed CentralPubMedView ArticleGoogle Scholar
- Bovet L, Feller U, Martinoia E: Possible involvement of plant ABC transporters in cadmium detoxification: a cDNA sub-microarray approach. Environ Int. 2005, 31: 263-267.PubMedView ArticleGoogle Scholar
- Herbette S, Taconnat L, Hugouvieux V, Piette L, Magniette ML, Cuine S, Auroy P, Richaud P, Forestier C, Bourguignon J, Renou JP, Vavasseur A, Leonhardt N: Genome-wide transcriptome profiling of the early cadmium response of Arabidopsis roots and shoots. Biochimie. 2006, 88: 1751-1765.PubMedView ArticleGoogle Scholar
- Kramer U: Metal hyperaccumulation in plants. Annu Rev Plant Biol. 2010, 61: 517-534.PubMedView ArticleGoogle Scholar
- Mendoza-Cózatl DG, Jobe TO, Hauser F, Schroeder JI: Long distance transport, vacuolar sequestration, tolerance, and transcriptional responses induced by cadmium and arsenic. Curr Opin Plant Biol. 2011, 14: 554-562.PubMed CentralPubMedView ArticleGoogle Scholar
- Barley: Production, Improvement, and Uses. Edited by: Ullrich SE. 2010, New Jersey: Wiley-BlackwellGoogle Scholar
- Chen F, Wang F, Wu FB, Mao WH, Zhang GP, Zhou MX: Modulation of exogenous glutathione in antioxidant defense system against Cd stress in the two barley genotypes differing in Cd tolerance. Plant Physiol Bioch. 2010, 45: 663-672.View ArticleGoogle Scholar
- Chen F, Wang F, Sun HY, Cai Y, Mao WH, Zhang GP, Vincze E, Wu FB: Genotype-dependent effect of exogenous nitric oxide on Cd-induced changes in antioxidative metabolism, ultrastructure, and photosynthetic performance in barley seedlings (Hordeum vulgare). J Plant Growth Regul. 2010, 29: 394-408.View ArticleGoogle Scholar
- Irizarry RA, Hobbs B, Collin F, Beazer-Barclay , Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264.PubMedView ArticleGoogle Scholar
- Menke M, Chen IP, Angelis KJ, Schubert I: DNA damage and repair in Arabidopsis thaliana as measured by the comet assay after treatment with different classes of genotoxins. Mutat Res. 2001, 493: 87-93.PubMedView ArticleGoogle Scholar
- Chen L, Dodd IC, Davies WJ, Wilkinson S: Ethylene limits abscisic acid- or soil drying-induced stomatal closure in aged wheat leaves. Plant Cell Environ. 2013, 36: 1850-1859.PubMedView ArticleGoogle Scholar
- Rodriguez-Serrano M, Romero-Puertas MC, Zabalza A, Corpas FJ, Gomez M, Del Rio LA, Sandalio LM: Cadmium effect on oxidative metabolism of pea (Pisum sativum L.) roots. Imaging of reactive oxygen species and nitric oxide accumulation in vivo. Plan Cell Environ. 2006, 29: 1532-1544.View ArticleGoogle Scholar
- Roxas VP, Lodhi SA, Garrett DK, Mahan JR, Allen RD: Stress tolerance in transgenic tobacco seedlings that overexpress glutathione S-transferase/glutathione peroxidase. Plant Cell Physiol. 2000, 41: 1229-1234.PubMedView ArticleGoogle Scholar
- Serrano R: Structure and function of plasma membrane ATPase. Annu Rev Plant Biol. 1989, 40: 61-94.View ArticleGoogle Scholar
- Williams LE, Mills RF: P1B-ATPases - an ancient family of transition metal pumps with diverse functions in plants. Trends Plant Sci. 2005, 10: 491-502.PubMedView ArticleGoogle Scholar
- Dietz KJ, Tavakoli N, Kluge C, Mimura T, Sharma SS, Harris GC, Chardonnens AN, Golldack D: Significance of the V-type ATPase for the adaptation to stressful growth conditions and its regulation on the molecular and biochemical level. J Exp Bot. 2001, 52: 1969-1980.PubMedView ArticleGoogle Scholar
- Finnegan EJ, Peacock WJ, Dennis ES: Reduced DNA methylation in Arabidopsis thaliana results in abnormal plant development. Proc Natl Acad Sci U S A. 1996, 93: 8449-8454.PubMed CentralPubMedView ArticleGoogle Scholar
- Fusco N, Micheletto L, Corso GD, Borgato L, Furini A: Identification of cadmium-regulated genes by cDNA-AFLP in the heavy metal accumulator Brassica juncea L. J Exp Bot. 2005, 56: 3017-3027.PubMedView ArticleGoogle Scholar
- Cai Y, Cao FB, Cheng WD, Zhang GP, Wu FB: Modulation of exogenous glutathione in phytochelatins and photosynthetic performance against Cd stress in the two rice genotypes differing in Cd tolerance. Biol Trace Elem Res. 2011, 143: 1159-1173.PubMedView ArticleGoogle Scholar
- Chaffei C, Pageau K, Suzuki A, Gouia H, Ghorbel MH, Masclaux-Daubresse C: Cadmium toxicity induced changes in nitrogen management in Lycopersicon esculentum leading to a metabolic safeguard through an amino acid storage strategy. Plant Cell Physiol. 2004, 45: 1681-1693.PubMedView ArticleGoogle Scholar
- Li JY, Fu YL, Pike SM, Bao J, Tian W, Chen CZ, Zhang Y, Li HM, Huang J, Li LG, Schroeder JI, Gassmann W, Gong JM: The Arabidopsis nitrate transporter NRT1.8 functions in nitrate removal from the xylem sap and mediates cadmium tolerance. Plant Cell. 2010, 22: 1633-1646.PubMed CentralPubMedView ArticleGoogle Scholar
- Schachtman DP, Schroeder JI, Lucas WJ, Anderson JA, Gaber RF: Expression of an inward-rectifying potassium channel by the Arabidopsis KAT1 cDNA. Science. 1992, 258: 1654-1658.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Bovet L, Eggmann T, Meylan-Bettex M, Polier J, Kammer P, Marin E, Feller U, Martinoia E: Transcript levels of AtMRPs after cadmium treatment: induction of AtMRP3. Plant Cell Environ. 2003, 26: 371-381.View ArticleGoogle Scholar
- Kim DY, Bovet L, Kushnir S, Noh EW, Martinoia E, Lee Y: AtATM3 is involved in heavy metal resistance in Arabidopsis. Plant Physiol. 2006, 140: 922-932.PubMed CentralPubMedView ArticleGoogle Scholar
- Fujita M, Fujita Y, Noutoshi Y, Takahashi F, Narusaka Y, Yamaguchi-Shinozaki K, Shinozaki K: Crosstalk between abiotic and biotic stress responses: a current view from the points of convergence in the stress signaling networks. Curr Opin Plant Biol. 2006, 9: 436-442.PubMedView ArticleGoogle Scholar
- Fuhrer J: Ethylene biosynthesis and cadmium toxicity in leaf tissue of beans (Phaseolus vulgaris L.). Plant Physiol. 1982, 70: 162-167.PubMed CentralPubMedView ArticleGoogle Scholar
- Zielinski RE: Calmodulin and calmodulin-binding proteins in plants. Annu Rev Plant Biol. 1998, 49: 697-725.View ArticleGoogle Scholar
- Garnier L, Simon-Plas F, Thuleau P, Agnel JP, Blein JP, Ranjeva R, Montillet JL: Cadmium affects tobacco cells by a series of three waves of reactive oxygen species that contribute to cytotoxicity. Plant Cell Environ. 2006, 29: 1956-1969.PubMedView ArticleGoogle Scholar
- Cosio C, Martinnoia E, Keller C: Hyperaccumulation of cadmium and zinc in Thlaspi caerulescens and Arabidopsis halleri at the leaf cellular level. Plant Physiol. 2004, 134: 716-725.PubMed CentralPubMedView ArticleGoogle Scholar
- Xiong J, An LY, Lu H, Zhu C: Exogenous nitric oxide enhances cadmium tolerance of rice by increasing pectin and hemicellulose contents in root cell wall. Planta. 2009, 230: 755-765.PubMedView ArticleGoogle Scholar
- Lee EJ, Matsumura Y, Soga K, Hoson T, Koizumi N: Glycosyl hydrolases of cell wall are induced by sugar starvation in Arabidopsis. Plant Cell Physiol. 2007, 48: 405-413.PubMedView ArticleGoogle Scholar
- Békésiová B, Hraška Š, Libantová J, Moravčíková J, Matušíková I: Heavy-metal stress induced accumulation of chitinase isoforms in plants. Mol Biol Rep. 2008, 35: 579-588.PubMedView ArticleGoogle Scholar
- Bertin G, Averbeck D: Cadmium: cellular effects, modifications of biomolecules, modulation of DNA repair and genotoxic consequences (a review). Biochimie. 2006, 88: 1549-1559.PubMedView ArticleGoogle Scholar
- Jenuwein T, Allis CD: Translating the histone code. Science. 2001, 293: 1074-1080.PubMedView ArticleGoogle Scholar
- Sobkowiak R, Deckert J: Proteins induced by cadmium in soybean cells. J Plant Physiol. 2006, 163: 1203-1206.PubMedView ArticleGoogle Scholar
- Ding YF, Chen Z, Zhu C: Microarray-based analysis of cadmium-responsive microRNAs in rice (Oryza sativa). J Exp Bot. 2011, 60: 3563-3573.View ArticleGoogle Scholar
- Zhou ZS, Song JB, Yang ZM: Genome-wide identification of Brassica napus microRNAs and their targets in response to cadmium. J Exp Bot. 2012, 63: 4597-4613.PubMed CentralPubMedView ArticleGoogle Scholar
- Huang SQ, Peng J, Qiu CX, Yang ZM: Heavy metal-regulated new microRNAs from rice. J Inorg Biochem. 2009, 103: 282-287.PubMedView ArticleGoogle Scholar
- Xu L, Wang Y, Zhai LL, Xu YY, Wang LJ, Zhu XW, Gong YQ, Yu RG, Limera C, Liu LW: Genome-wide identification and characterization of cadmium-responsive microRNAs and their target genes in radish (Raphanus sativus L.) roots. J Exp Bot. 2013, 64: 4271-4287.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.