- Research article
- Open Access
Expression profiling of Crambe abyssinicaunder arsenate stress identifies genes and gene networks involved in arsenic metabolism and detoxification
BMC Plant Biology volume 10, Article number: 108 (2010)
Arsenic contamination is widespread throughout the world and this toxic metalloid is known to cause cancers of organs such as liver, kidney, skin, and lung in human. In spite of a recent surge in arsenic related studies, we are still far from a comprehensive understanding of arsenic uptake, detoxification, and sequestration in plants. Crambe abyssinica, commonly known as 'abyssinian mustard', is a non-food, high biomass oil seed crop that is naturally tolerant to heavy metals. Moreover, it accumulates significantly higher levels of arsenic as compared to other species of the Brassicaceae family. Thus, C. abyssinica has great potential to be utilized as an ideal inedible crop for phytoremediation of heavy metals and metalloids. However, the mechanism of arsenic metabolism in higher plants, including C. abyssinica, remains elusive.
To identify the differentially expressed transcripts and the pathways involved in arsenic metabolism and detoxification, C. abyssinica plants were subjected to arsenate stress and a PCR-Select Suppression Subtraction Hybridization (SSH) approach was employed. A total of 105 differentially expressed subtracted cDNAs were sequenced which were found to represent 38 genes. Those genes encode proteins functioning as antioxidants, metal transporters, reductases, enzymes involved in the protein degradation pathway, and several novel uncharacterized proteins. The transcripts corresponding to the subtracted cDNAs showed strong upregulation by arsenate stress as confirmed by the semi-quantitative RT-PCR.
Our study revealed novel insights into the plant defense mechanisms and the regulation of genes and gene networks in response to arsenate toxicity. The differential expression of transcripts encoding glutathione-S-transferases, antioxidants, sulfur metabolism, heat-shock proteins, metal transporters, and enzymes in the ubiquitination pathway of protein degradation as well as several unknown novel proteins serve as molecular evidence for the physiological responses to arsenate stress in plants. Additionally, many of these cDNA clones showing strong upregulation due to arsenate stress could be used as valuable markers. Further characterization of these differentially expressed genes would be useful to develop novel strategies for efficient phytoremediation as well as for engineering arsenic tolerant crops with reduced arsenic translocation to the edible parts of plants.
Arsenic (As) is ubiquitous in the earth's crust and has widely been known as an acute poison and carcinogen [1, 2]. Chronic As ingestion from contaminated food supplies and drinking water is strongly associated with an increased risk of liver, kidney, skin, and lung cancer . Inorganic forms of arsenic, pentavalent arsenate (AsV) and trivalent arsenite (AsIII), are more toxic than the organic forms and are highly abundant in the environment. High levels of As in soil and drinking water have been reported around the world including Bangladesh, India, USA, and South America, where As in drinking water exceeds 10 ppb, the permissible limit established by the World Health Organization . Industrial and agricultural practices, including the use of As in pesticides and herbicides , wood preservatives , and irrigation with contaminated groundwater [6, 7], have significantly increased As levels in agricultural soils. Several studies indicated that food crops such as rice and vegetables grown on As contaminated soil can accumulate high levels of As in roots, shoots and seeds [7–9]. Thus As uptake by crop plants plays an important role in transfer of this toxic element into the food chain. Additionally, inorganic As species are phytotoxic and the elevated concentration of As in the soil causes a significant reduction in crop yield [7, 10, 11]. Therefore, developing strategies for remediation of As contaminated agricultural soil and water and decreasing the As contents in the food crops for preventing As contamination of the food chain are highly desirable. Furthermore, the phytotoxic effects suffered by crops grown in soil with As residues could be overcome by developing crops resistant to As. Progress towards developing such genetics-based strategies, however, has been hindered by the lack of a basic understanding of the molecular and biochemical mechanisms of As uptake and detoxification in higher plants.
Recently, plants have been engineered to detoxify and accumulate As in the aboveground tissues. Dhankher et al  engineered Arabidopsis with bacterial arsenate reductase, ArsC, and γ-glutamylecysteine synthatase, γ-ECS, thereby increasing the As accumulation and tolerance. This innovative study showed that the pathways for As metabolism could be engineered in plants either for hyperacumulation and tolerance or preventing the uptake of As in the food chain once we better understand the details of As metabolism in plants. While the molecular mechanism of As detoxification and tolerance has yet to be fully determined in plants, it has been shown that plants detoxify As by reducing AsV to AsIII, which is subsequently detoxified via forming complexes with thiol-reactive peptides such as γ-glutamylcysteine (γ-EC), glutathione (GSH), and phytochelatins (PCs) [13–16]. The AsIII-thiol complexes are suggested to be sequestered into vacuoles by glutathione-conjugating pumps, GCPs, although this still remains to be proven.
Recently, an arsenate reductase in Arabidopsis (AtACR2), which is a CDC25 phosphatase homologue, has been isolated and characterized [17, 18]. The RNAi-mediated silencing of AtACR2 increased the sensitivity of Arabidopsis to AsV but not to AsIII, and the ACR2-deficient Arabidopsis plants translocated significantly higher amount of As in the aboveground tissues as compared to wild type plants . Arsenate reductases have also been isolated from P. vitatta and rice that complemented the AsV sensitivity of yeast mutants deficient in arsenate reductase, ACR2 [19, 20]. In addition, a glutaredoxin, PvGrx5, has been reported to have a role in regulating the intracellular levels of AsIII in P. vittata .
A comparative transcriptome analysis of AsV and AsIII stress in rice seedlings showed differential expression of a wide array of transcripts encoding transporters, regulatory genes, and genes involved in various stress responses . Recently, two silicon transporters, LsiI and LsiII, which belong to the nodulin 26-like intrinsic proteins (NIPs) subfamily of the plant aquaporins, have been shown to co-transport AsIII and thus playing a major role in the entry of AsIII into rice roots . Other members of NIP subfamily from Arabidopsis, rice, and lotus have been shown to facilitate the bidirectional transport of AsIII in yeast lacking the endogenous aquaglyceroporin transporter, FPS1 .
Several plant species, especially members of Brassicaseae and pteridophytes such as P. vittata, have been identified as promising candidates for phytoremediation [25, 26]. In spite of all the recent As-related studies, we are still far from a comprehensive understanding of As uptake, detoxification, and sequestration in plants including P. vittata. Therefore, a better understanding of the molecular and biochemical mechanism of As uptake and detoxification in plants is needed to develop cost-effective strategies for phytoremediation as well as for engineering of As resistant food crops with reduced As accumulation. C. abyssinica, a member of Brassicaceae commonly known as 'abyssinian mustard', is a non-food, high biomass crop with short life cycle that holds potential for efficient phytoremediation of arsenic contaminated soils and sediments. C. abysinnica, referred to here as crambe, can be grown wherever oilseed rape is grown and it does not out-cross with canola or any other agriculture crop [27, 28]. Moreover, crambe is highly resistant to As and other heavy metals, and accumulates significantly higher levels of As as compared to other Brassica species . For these reasons, we selected crambe to study the molecular mechanism of As accumulation and detoxification and identify genes that are differentially regulated by AsV exposure. In this study, we isolated differentially expressed cDNAs from crambe exposed to AsV employing the PCR-Select Suppression Subtraction Hybridization (SSH) approach. We report here many uncharacterized novel transcripts as well as those participating in known biochemical pathways that are differentially expressed due to arsenate exposure.
Results and Discussion
Arsenic accumulation and tolerance in crambe
In order to select a plant suitable for studying As tolerance, we compared As accumulation in crambe plants to that of five other Brassica family members, namely Brassica napus, B. compestris cv Turkey, B. compestris cv Yellow Sarson, B. oleracea and B. rapa. Crambe plants were more tolerant to As than the various Brassica species examined and accumulated significantly higher amounts of As in aboveground shoot tissues followed by B. compestris and B. rapa (Figure 1). Similar results have been shown for uptake of As and Cd in crambe . This suggests that crambe has an inherent capacity to accumulate As and other toxic metals at higher levels than other Brassica species. Thus, crambe plants are highly suitable for examining AsV stress response. To determine the optimum concentration of AsV required for inducing plant response, we treated 10 days old crambe seedlings with varying concentrations of sodium arsenate (0 to 350 μM) for seven days. The fresh weight of plants decreased moderately at 200 and 250 μM sodium arsenate concentrations with no symptoms of severe cellular damage (Figure 2). At higher concentrations (300 and 350 μM), however, AsV was toxic to plants as they were severely stunted with highly reduced biomass and had visible injury symptoms on leaves. Low levels of AsV (100 and 150 μM) have little effect on fresh weight or plant morphology. Therefore, an effective concentration of 250 μM sodium arsenate was selected to treat crambe seedlings to induce AsV stress.
Differentially expressed crambe transcripts in response to arsenate stress
We used a PCR-Select SSH approach to isolate the differentially expressed transcripts from crambe in response to AsV exposure. As the ratio of driver to tester affects the stringency of the subtraction, two parallel subtractions were performed varying the driver/tester ratio and thus a total of 435 subtracted cDNA clones were obtained. A differential screening procedure finally identified 105 clones containing the differentially expressed mRNA sequences eliminating a large number of false positives as visible on the blots (Figure 3). We subjected all 105 subtracted cDNA sequences to homology search against the NCBI nucleotide database and found that they represented 38 different coding transcripts (Table 1). These crambe cDNA sequences were highly similar to that of Arabidopsis, as both species belong to the Brassicaceae family. A number of transcripts in our subtracted cDNA library were homologous to well characterized proteins, however several transcripts encode uncharacterized proteins that could be involved in novel pathways of As metabolism in plants.
The identity of the subtracted cDNAs suggested that AsV stress significantly affects pathways related to oxidative stress, defense, ion transport, sulfur assimilation, signal transduction, and photosynthesis and metabolism. Among the subtracted cDNA sequences (Table 1), glutathione-S-transferase (GSTs, both tau and phi subfamilies) was represented by the maximum number of clones (16%) followed by the heat shock proteins (15%), various reductases (15%), and proteins in the ubiquitin pathway (8%) indicating the significance of AsV as an acute cytotoxin. An iron-binding oxidoreductase (FeBOx) that participates in phenylpropanoid metabolism and has a functional domain similar to that of anthocyanidin synthase  represented 12% of cDNA clones in the library. Most of the remaining sequences encode proteins involved in various metabolic pathways including sulfate assimilation, carbohydrate and nitrogen metabolism, and amino acid synthesis. AsV stress also induced several general stress and defense-related proteins. Transcripts that code for a cation transporter associated protein, ChaC, and a secondary multidrug and toxin extrusion pump (MATE) along with many proteins of unknown function were also present in the library.
An arsenic-induced transcription profile has recently been reported in Arabidopsis  and rice  and identified GSTs as highly represented proteins in response to AsV exposure, followed by antioxidant synthesis and proteins involved in sulfur transport and metabolism. Our results for the transcriptional response of crambe to AsV support this consensus. In addition, our library contained the transcripts encoding glucanases, RNA binding proteins, proteins involved in amino acid biosynthesis as reported in Arabidopsis , as well as HSPs, MATE, glucosidases, jasmonic acid synthesis, and glyceraldehyde 3-phosphate dehydrogenase (GAPDH) as reported in rice [22, 32]. However, almost 50% of the sequences identified in our subtraction library have not been reported in As related expression studies in other plant species. Important among them, were the serine palmitoyl transferase (SPT), adenylyl phosphosulfokinase (APSK), ChaC, MATE, pathogenesis and defense related proteins (PDR), transport protein particle component Bet3p like (TRAPP beta3) involved in vesicular trafficking and the proteins with unknown functions. Based on the sequence identity and their functions in various metabolic pathways, we clustered the transcripts into different functional groups as discussed below.
Oxidative stress and antioxidants
Proteins involved in oxidative stress were highly represented in our subtracted cDNA library. There were four transcripts with 75% identity to the known aldo/keto reductases (AKR) and five transcripts similar to peptide methionine sulfoxide reductases (PMSR). AKR and PMSR have established roles in mitigating oxidative protein damage. Arabidopsis plants over-expressing PMSR3 were tolerant to methyl viologen, ozone, and metals such as Zn and Mn, illustrating its role as a vital antioxidant . Recently, members of the NADPH dependent AKR family have been portrayed as antioxidants as well as cellular detoxification enzymes . Arsenic has been shown to interact with phospholipid bilayers and thus modifying the properties of biomembranes by altering membrane fluidity . Serine palmitoyl transferase (SPT), which catalyses the initial step of sphingolipid biosynthesis was present in our library, and was upregulated by AsV (Figure 4) and ozone treatments (Table 1 and 2). The lethal phenotype of SPT knockouts and reduced growth of Arabidopsis lines with RNAi downregulated SPT, illustrates the vital role of these molecules in membrane integrity .
The jasmonic acid signaling pathway is thought to mediate biotic as well as abiotic stress signaling. The 12-oxophytodienoate (OPDA) with structural resemblance to jasmonate is the first metabolic intermediate in the jasmonate synthesis pathway. It has been shown that the OPDA was capable of inducing the OPR, several GSTs and genes containing xenobiotic responsive elements . Our subtracted library contains multiple copies of OPDA reductase (OPR), which was upregulated both in arsenate (Table 1, Figure 4) and ozone treatment (Table 2). Recently, members of the OPR family (OPR1, OPR2, and OPR3) have been shown to reduce the xenobiotic 2,4,6-trinitrotoluene (TNT) in Arabidopsis . These studies and the upregulation of OPR in response to AsV exposure in crambe suggest the significance of OPRs in signaling and detoxification of toxic pollutants.
Glutathione-S-Transferases and sulfate assimilation
Transcripts belonging to the glutathione-S-transferases (GSTs) form the largest group in the subtracted cDNA library. Most of these transcripts (12%) fall in the Tau subfamily (GST-Tau), while the remaining sequences (4%) are similar to the Phi subfamily. GSTs are considered to be involved in detoxification of both endogenous and xenobiotic compounds with electrophilic centers by the nucleophilic addition of glutathione . Members of the ABC transporter family, human multidrug resistant proteins (MRPs) and yeast cadmium factor1 (YCF1), were found to export glutathione conjugated AsIII out of the cytoplasm [40, 41]. However, GST mediated conjugation of glutathione with inorganic ions has not been proven. GSTs might be mitigating oxidative stress by other means as ectopic expression of tomato GST significantly improved the tolerance of yeast to H2O2 induced oxidative stress .
Another sequence present in our library, monodehydroascorbate reductase (MDAR) participates in the glutathione-ascorbate redox cycle that is responsible for the active oxygen species (AOS) scavenging process . Ascorbate reacts with AOS and produces an unstable compound, monodehydroascorbate (MDHA) that dissociates into ascorbate and dehydroascorbate (DHA). The later is eventually reduced to ascorbate by DHA reductase using reduced glutathione as the source of electrons . However, the increased regeneration of ascorbate directly from MDHA due to higher levels of MDAR could decrease the demand for glutathione . It has been shown that tobacco seedlings overexpressing a tobacco GST exhibited a concomitant increase in the expression of MDAR, possibly due to a lower amount of reduced glutathione . Surprisingly, both the GST-Tau and MDAR in our subtracted library displayed identical time-dependent expression patterns in response to arsenate (Figure 4), which is in agreement with a possible role of diminishing cellular concentration of GSH on induction of MDAR.
Our library also contains genes related to biosynthesis of cysteine, a precursor of glutathione biosynthesis, thus emphasizing the role of GSH in arsenic detoxification. Cysteine biosynthesis represents the reductive sulfate assimilation pathway. This is initiated by adenylylation of sulfate by sulfate adenylyl transferase or adenosine phosphosulfurylase (APS) followed by its reduction to sulfite by APS reductase (APR) and then to sulfide by sulfite reductase (SiR). The sulfide is incorporated into O-acetyl serine by cysteine synthase . Considering the occurrence of elevated levels of GSH in plants overexpressing cysteine synthase , the presence of highly upregulated SiR and APS in our library prompted us to speculate that there is an overall increase in glutathione synthesis under AsV stress. Further studies are needed to establish a clear role for these enzymes in As metabolism.
Genes involved in the proteasomal degradation system constitutes 8% of the subtracted cDNAs in our library. This includes the proteolytic beta-subunit of the 20S particle, polyubiquitin 14 (UBQ14) and an ubiquitin-associated protein. The elevated expression of 20S proteasome subunits and UBQ14 indicated that, most likely, the oxidized or the arsenite bound proteins in AsV stressed plant cells were ubiquitinated for proteasome-mediated degradation. Besides, the ubiquitin-associated protein present in our subtraction library contained a zinc ion-binding domain, which could be a crambe homologue of an Ubrp1 type of ubiquitin ligase (E3) that transfers the activated ubiquitin to target proteins . Therefore, we suggest that the proteasome-mediated protein degradation machinery is vital to counteract the toxicity induced by As. However, further studies are needed to elucidate the exact roles of the genes involved in the proteasomal degradation system in As metabolism and detoxification.
Two membrane transporters, a multidrug and toxin extrusion pump (MATE) and a putative cation transporter associated protein (ChaC), were also found in the pool of our subtracted cDNA clones. MATE belongs to a large family of secondary transporters with 56 ORFs in the A. thaliana genome. An Arabidopsis MATE transporter, AtDTX1, was reported to transport alkaloids, antibiotics and other toxic compounds . In addition, ectopic expression of AtDTX1 conferred Cd resistance to bacterial cells and two members of this family were upregulated in As treated rice seedlings . In prokaryotes, the putative transporter regulatory protein ChaC is thought to be associated with a Ca2+/H+ cation transporter, ChaA. A ChaC homologue from rice, OsARP, conferred salt tolerance when overexpressed in tobacco , suggesting a role for ChaC in abiotic stress tolerance. These transporters have not yet been characterized in plants and further studies are needed to understand their role in arsenic metabolism.
Unknown novel proteins
There were five subtracted cDNA sequences that did not match to any of the characterized proteins in the database. However, they have homologs in Arabidopsis with very high sequence similarity, but with unknown function, and these transcripts were named as UK (unknown) with their colony ID. For example, UK325, one of the novel sequences, has a methyl transferase and a Suppressor of variegation-Enhancer of zeste-Trithorax (SET) domains for transferring a methyl group to its substrates from S-adenosyl methionine. Previously characterized arsenic methyl transferases in rabbit liver demonstrated a GSH-dependent transfer of methyl group from S-adenosyl methionine to As . Further characterization of these unknown novel sequences will shed light on their role in As metabolism.
Validation of PCR-Select subtraction data by RT-PCR
In order to confirm the upregulation of the subtracted cDNA clones in response to AsV exposure, a semi-quantitative RT-PCR was performed for each gene present in the subtracted cDNA library. Nearly 85% of the transcripts present in our library responded within 6 hrs after the AsV treatment (Figure 4). Both AKR and PMSR proteins were upregulated within 6 hrs after the As treatment and showed a similar expression pattern up to 24 hrs. MDAR and OPR responded within 6 hrs of AsV treatment and transcripts of MDAR continued to increase several folds higher at 24 hrs (Figure 4). A similar pattern of transcriptional response was observed in the case of SPT, known to participate in the ROS response. Transcripts of both subfamilies of GSTs were upregulated within 6 hrs of AsV treatment and the mRNA levels remained high throughout the period of As treatment.
Among the enzymes in sulfate assimilation, SiR showed strong induction throughout the duration of AsV treatments. The transcripts level of APSK increased consistently, however, APS was upregulated only at 12 hrs and then decreased after 24 hrs. Since the energetics of the activation of sulfate favors the reverse reaction , a delay in the induction of APS was expected as compared to the downstream enzymes in sulfate assimilation. The UBQ14 and 20S beta-subunit were significantly upregulated at 6 hrs and this higher expression level remained unaltered at all the time points tested. The expression pattern of the ubiquitin-associated protein was noticeably different, however, with a sharp increase in mRNA levels at 6 hrs and then decreased at the following time intervals. The secondary transporter, MATE, was induced in As treated seedlings only after 24 hrs of exposure to AsV, whereas, ChaC was induced at 6 hrs and continued the higher level of expression at 12 and 24 hrs. The unknown sequences, UK17 and UK231 were consistently upregulated at 6, 12 and 24 hrs after AsV treatment, whereas, UK21 and UK325 showed maximum mRNA levels at 6 hrs and then decreased steadily over time. The unknown protein with a domain similar to the luminal binding protein (UK-lum) was upregulated only after 12 hrs with further increase in expression at 24 hrs.
The expression levels of the sequences corresponding to the glycoprotein EDGP, FeBOx, HSP21, β-tubulin, and TRAPP beta3 were upregulated and maintained many fold higher expression level at 6, 12, and 24 hrs after AsV treatment. A similar pattern but at lower expression levels was observed in case of thioglucosidase, tryptophan synthase, cysteine proteinase, and aldose epimerase. Contrary to this trend, both trypsin inhibitor and defense related proteins were strongly induced only after 24 hrs of AsV exposure. The differential regulation of most of these subtracted cDNA suggests their role in As metabolism and detoxification processes.
Comparison with phosphate and ozone stress
A comparison of the expression of the homologous transcripts in Arabidopsis in response to phosphate deficiency and ozone-induced oxidative stress revealed an overlapping effect of these related stress conditions (Table 2). Expression levels of one third of the sequences in the subtracted cDNA library were affected by varying phosphate level in the media and strikingly, all of them were found to have higher mRNA levels in phosphorus deficiency. This trend was verified in crambe by semi-quantitative RT-PCR in seedlings grown with or without phosphate for 6, 12, and 24 hrs as short-term and 3 and 7 days as long-term exposure. In crambe, only a few transcripts such as cysteine proteinase, thioglucosidase, aldose epimerase, and PMSR showed increased transcription within 24 hrs of phosphate deficiency (Figure 5). However, the transcripts encoding EDGP, AKR, trypsin inhibitor, FBPA, GSTs, UK17, tryptophan synthase, defense related protein, UK-lum, and UK21 were upregulated after 7 days of phosphorus deficiency.
Ozone is known to induce oxidative stress and transcripts encoding enzyme with antioxidant property were upregulated in response to ozone stress. One-third of the transcripts related to the subtracted cDNAs in our library showed elevated expression levels when Arabidopsis seedlings were subjected to oxidative stress by ozone (Table 2). This supported the previously reported observations of As-induced ROS resulting in lipid peroxidation and higher SOD activity in Holcus lanatus .
Present study proves that crambe plants have the potential for As phytoremediation because they are highly tolerant to As and in addition can accumulate high levels of this toxic metal in the above ground parts. Moreover, our study identifies several genes that are highly upregulated in response to AsV stress and that encode novel proteins and proteins involved in various vital metabolic pathways. These arsenate-induced proteins may play a crucial role in As metabolism in plants. By further characterizing these differentially expressed genes, we may be able to acquire knowledge essential for developing novel strategies for efficient phytoremediation as well as engineering As tolerant crop plants with reduced translocation of this deadly metalloid to the edible parts of plants.
Plant growth and arsenate treatments
Crambe abyssinica (cv. BelAnn) seeds were surface disinfected using 70% ethanol for five minutes followed by a treatment with 30% clorox for 30 minutes. Seeds were washed four times with sterile deionized water and then germinated on 0.5 × MS (Murashige and Skoog) agar media. Germinated seeds were aseptically transferred to 0.5 × MS liquid media in 250 ml flasks and grown under controlled conditions (16/8 hour light/dark periods at 22°C) with constant shaking. After 10 days, sodium arsenate (Na2HAsO4.7H2O) at various concentrations (0, 100, 150, 200, 250, 300 and 350 μM) was added to the medium and plants were grown for an additional 7 days to determine the optimum concentration for As treatment. A final concentration of 250 μM AsV, where plants showed significant reduction in biomass without exhibiting severe toxicity symptoms, was then used for treating plants for cDNA subtraction experiments. Ten days old seedlings were exposed to AsV for 24 hrs and the seedlings were harvested, frozen in liquid nitrogen, ground in to a powder and then stored at -80°C for further processing. For inducing the phosphate (P) deficiency, plants were grown in 0.5 × MS phosphate dropout liquid media (pH 7.0) in 250 ml flasks under 16/8 hour light/dark periods at 22°C with constant shaking. After the various time periods (6, 12, and 24 hrs as short-term and 3 and 7 days as long-term exposure), seedling samples were harvested, washed in deionized water and stored at -80°C till further use.
In order to compare the As accumulation in different Brassica species (Brassica napus, B. compestris cv Turkey, B. compestris cv Yellow Sarson, B. oleracea, B. rapa) and C. abyssinica, the seeds were germinated and grown for 14 days on 0.5 × MS agar media supplemented with various concentrations (0, 200 and 400 μM) of sodium arsenate. Total As in the shoot tissues was analyzed using Inductively Coupled Plasma-Mass Spectroscopy (ICP-MS).
RNA isolation, cDNA synthesis and PCR-Select subtraction
Total RNA was extracted from control and AsV treated samples using RNeasy Plant Mini kit (Qiagen) followed by the mRNA isolation using NucleoTrap mRNA isolation kit. Double stranded cDNA synthesis and cDNA subtractions were performed using the PCR select cDNA subtraction kit (Clontech, CA) according to the manufacturer's protocol. In short, the cDNAs from driver (control) and tester (As-treated) samples were digested with endonuclease RsaI and the tester cDNAs were ligated to two different adapters at either end of the cDNAs separately . The twin samples of tester differing in the adapters were denatured and allowed to hybridize in separate tubes with excess denatured driver cDNAs in each tube. Apart from enriching the differentially expressed cDNAs, this step normalized the high and low abundant sequences. The tester samples were then mixed without denaturation in presence of the fresh denatured driver cDNAs and allowed to anneal overnight. Subsequent PCR amplifications, first with primers corresponding to both adapters then with nested primers, selectively amplified and enriched the differentially expressed cDNAs. A reverse subtraction, by switching the tester and driver cDNAs was also performed for the purpose of differential screening. In order to construct the subtracted cDNA library, the subtracted and PCR amplified products were ligated to the pGEM-T-Easy T/A cloning vector (Promega) and transformed into DH5α cells.
A differential screening procedure was used to eliminate the false positives in the subtracted cDNA library. This step was carried out using a differential screening kit (Clontech) following the manufacturer's protocol. In this procedure, the subtracted cDNAs blotted onto the nylon membranes were hybridized with forward and reverse subtracted probes and the signals were compared for each clone. Colonies containing plasmids with cloned cDNAs were grown in liquid LB medium overnight and an equal amount of each culture was spotted in duplicate on two replicate nylon membranes, which were placed on LB plates and incubated at 37°C overnight. Membranes containing the colonies were transferred to Whatmann filters pre-saturated with denaturing solution (0.5 M NaOH, 1.5 M NaCl), and then to Whatmann filters pre-saturated with neutralizing solution (0.5 M Tris-HCl, 1.5 M NaCl). Each membrane was air-dried and the DNA was fixed by baking the membranes at 80°C for 2 hours. These membranes were then hybridized with either 32P labeled forward or reverse subtracted probes overnight as described previously  and washed with low stringency (2 × SSC, 0.5% SDS) and high stringency (0.2 × SSC, 0.5% SDS) wash buffer before exposing to X-ray film at -80°C overnight. The hybridization signal strength for each cDNA clone on membrane probed with forward subtracted cDNA was compared with cDNA clones on the duplicate membrane probed with reverse subtracted cDNA.
DNA sequencing and gene identification
A total of 105 differentially expressed cDNA clones were sequenced both in forward and reverse orientations. The DNA sequences of these clones were used to search the NCBI Genebank database to identify homologues sequences and their identity using the BLAST algorithm . All the subtracted cDNA sequences were used for BLAST search against both DNA and protein sequences on the database and their homologous gene sequences were identified both from Prokaryotes and Eukaryotes.
For RT-PCR, 10 days old crambe seedlings were treated with 250 μM sodium arsenate and the samples were harvested after 6, 12 or 24 hours treatment. Total RNA was extracted using the RNAeasy Plant mini kit (Qiagen), and cDNA synthesized from total RNA using Thermoscript RT-PCR kit (Invitrogen) following the manufacturers protocol. Primers for each subtracted cDNA were designed using Primer3 . Threshold of the number of cycles was determined by an optimization step in which the PCR was performed from a known amount of cDNA for a range of cycles. The number of PCR cycles selected for the experiments were chosen from the plot where the trend line exhibited highest correlation to exponential amplification. Once the number of cycles for a few genes was optimized, it was extrapolated to all other genes in the library based on their microarray expression (ArrayExpress) values, as it represented the relative expression level of each mRNA. In addition, the constitutively expressed actin2 gene, ACT2, was used as an internal control to verify that an equal amount of cDNA template was used in each PCR reaction. A touchdown PCR was performed for all the transcripts with annealing temperature varying from 64°C to 52°C. Similarly, to examine the effect of phosphate deficiency on gene expression, seedlings were grown in half strength MS medium with or without phosphate and after 6, 12 or 24 hrs treatment samples were collected. In order to evaluate the long-term effects of phosphate deficiency, RNA was isolated from seedlings grown in phosphate deficient media for 3 and 7 days.
Statistical analysis and comparison with microarray data
All statistical analyses were performed using SAS software (SAS version 9.1). For comparison, Arabidopsis microarray data on phosphate [, ArrayExpress ID: E-MEXP-791] and ozone treatments [, ArrayExpress ID: E-MEXP-342], were downloaded from the ArrayExpress , processed by Bioconductor in R programming environment and annotated using a Perl script. ANOVA was performed on the log-transformed expression values for genes which were present in the subtracted library (25 transcripts in Affymetrix array) identified by AGI locus and followed by Dunnett's procedure  for ozone treatments. Duncan's New Multiple range test  was performed for the pair-wise comparison of expression levels in high phosphate and low phosphate treatments.
OPD is an Assistant Professor and BP is a PhD candidate in the Department of Plant, Soil, and Insect Sciences, University of Massachusetts Amherst. SK is currently an undergraduate student at Cornell University.
Kaiser J: Second look at arsenic finds higher risk. Science. 2001, 293: 2189-10.1126/science.293.5538.2189b.
Bhattacharjee Y: A sluggish response to humanity's biggest mass poisoning. Science. 2007, 315: 1659-1661. 10.1126/science.315.5819.1659.
World Health Organization (WHO): Arsenic in drinking water, Fact sheet No. 210. [http://www.who.int/int-fs/en/fact210.html]
Baker RS, Arle HF, Miller JH, Holstun JT: Effects of organic arsenical herbicides on cotton response and chemical residues. Weed Sci. 1969, 17: 37-40.
US Department of Agriculture: Wood Preservatives. The Pesticide Review. Washington DC. 1974
Zhao FJ, McGrath SP, Meharg AA: Arsenic as a Food Chain Contaminant: Mechanisms of Plant Uptake and Metabolism and Mitigation Strategies. Annu Rev Plant Biol. 2010,
Zhu YG, Williams PN, Meharg A: Exposure to inorganic arsenic from rice: a global health issue?. Environ Pollut. 2008, 154: 169-171. 10.1016/j.envpol.2008.03.015.
Bennett JP, Chiriboga E, Coleman J, Waller DM: Heavy metals in wild rice from Northern Wisconsin. Sci Total Environ. 2000, 246: 261-269. 10.1016/S0048-9697(99)00464-7.
Das HK, Mitra AK, Sengupta PK, Hossain A, Islam F, Rabbani GH: Arsenic concentrations in rice, vegetables, and fish in Bangladesh: a preliminary study. Environ Int. 2004, 30: 383-387. 10.1016/j.envint.2003.09.005.
Meharg AA: Arsenic in rice: understanding a new disaster for South-East Asia. Trends Plant Sci. 2004, 9: 415-417. 10.1016/j.tplants.2004.07.002.
Marin AR, Masscheleyn PH, Patrik WH:jr Soil redox-pH stability of arsenic and its influence on arsenic uptake in rice. Plant Soil. 1993, 152: 245-253. 10.1007/BF00029094.
Dhankher OP, Li Y, Rosen BP, Shi J, Salt D, Senecoff JF, Sashti NA, Meagher RB: Engineered tolerance and hyperaccumulation of arsenic in plants by combining arsenate reductase and γ-glutamylcysteine synthetase expression. Nature Biotechnol. 2002, 20: 1140-1145. 10.1038/nbt747.
Dhankher OP: Arsenic metabolism in plants: An inside story. New Phytol. 2005, 168: 503-505. 10.1111/j.1469-8137.2005.01598.x.
Pickering IJ, Prince RC, George MJ, Smith RD, George GN, Salt DE: Reduction and coordination of arsenic in Indian mustard. Plant Physiol. 2000, 122: 1171-1177. 10.1104/pp.122.4.1171.
Li Y, Dhankher OP, Carreira L, Balish R, Meagher R: Engineered over-expression of-glutamylcysteine synthetase in plants confers high levels arsenic and mercury tolerance. Environ Toxicol Chem. 2005, 24: 1376-1386. 10.1897/04-340R.1.
Raab A, Schat H, Meharg A, Feldmann J: Uptake, translocation, and tansformation of arsenate and arsenite in sunflower Helianthus annuus: formation of arsenic pytochelatin complexes during exposure to high arsenic concentrations. New Phytol. 2005, 168: 551-558. 10.1111/j.1469-8137.2005.01519.x.
Dhankher OP, Rosen BP, McKinney EC, Meagher RB: Enhanced arsenic uptake in Arabidopsis plants by suppressing endogenous arsenate reductase AtACR2 gene. Proc Nat Acad Sci USA. 2006, 103: 5413-5418. 10.1073/pnas.0509770102.
Bleeker PM, Hakvoort HW, Bliek M, Souer E, Schat H: Enhanced arsenate reduction by a CDC25-like tyrosine phosphatase explains increased phytochelatin accumulation in arsenate-tolerant Holcus lanatus. Plant J. 2006, 45: 917-929. 10.1111/j.1365-313X.2005.02651.x.
Ellis DR, Gumaelius L, Indriolo E, Pickering IJ, Banks JA, Salt DE: A novel arsenate reductase from the arsenic hyperaccumulating fern Pteris vittata. Plant Physiol. 2006, 141: 1544-1554. 10.1104/pp.106.084079.
Duan Y, Zhou YP, Mukhopadhyay R, Rosen BP, Zhu YG: A CDC25 homologous from rice functions as an arsenate reductase. New Phytol. 2007, 174: 31-321. 10.1111/j.1469-8137.2007.02009.x.
Sundaram S, Rathinasabapathi B, Ma LQ, Rosen BP: An arsenate-activated glutaredoxin from the arsenic hyperaccumulator fern Pteris vittata L. regulates intracellular arsenite. J Biol Chem. 2008, 283: 6095-6101. 10.1074/jbc.M704149200.
Chakrabarty D, Trivedi P K, Misra P, Tiwari M, Shri M, Shukla D, Kumar S, Rai A, Pandey A, Nigam D, Tripathi R, Tuli R: Comparative transcriptome analysis of arsenate and arsenite stresses in rice seedlings. Chemosphere. 2009, 74: 688-702. 10.1016/j.chemosphere.2008.09.082.
Ma JF, Yamaji N, Mitani N, Xu XY, Su YH, McGrath SP, Zhao FJ: Transporters of arsenite in rice and their role in arsenic accumulation in rice grain. Proc Nat Acad Sci USA. 2008, 105: 9931-9935. 10.1073/pnas.0802361105.
Bienert GP, Thorsen M, Schüssler MD, Nilsson HR, Wagner A, Tamás MJ, Jahn TP: A subgroup of plant aquaporins facilitate the bi-directional diffusion of As (OH)3 and Sb(OH)3 across membranes. BMC Biol. 2008, 6: 26-10.1186/1741-7007-6-26.
Ma LQ, Komar KM, Tu C, Zhang WH, Cai Y, Kennelley ED: A fern that hyperaccumulates arsenic: a hardy versatile fast-growing plant helps to remove arsenic from contaminated soils. Nature. 2001, 409: 579-579. 10.1038/35054664.
Salido AL, Hasty KL, Lim JM, Butcher DJ: Phytoremediation of arsenic and lead in contaminated soil using Chinese brake ferns Pteris vittata and Indian mustard Brassica juncea. Int J Phytorem. 2003, 5: 89-103. 10.1080/713610173.
Anon: Crambe production and utilization. BioOptions, Newsletter of the Centre for Alternative Plant and Animal Products, University of Minnesota. 1993, 1-3.
European Union Strategic Research Agenda, part II: Plant for the Future. Stakeholder proposal for a strategic research agenda 2005 including Draft Action Plan 2010. 2005
Artus NN: Arsenic and cadmium phytoextraction potential of crambe compared with Indian mustard. J Plant Nutr. 2006, 29: 667-679. 10.1080/01904160600564444.
Saito K, Kobayashi M, Gong Z, Tanaka Y, Yamazaki M: Direct evidence for anthocyanidin synthase as a 2-oxoglutarate-dependent oxygenase: molecular cloning and functional expression of cDNA from a red forma of Perilla frutescens. Plant J. 1999, 17: 181-189. 10.1046/j.1365-313X.1999.00365.x.
Abercrombie JM, Halfhill MD, Ranjan P, Rao MR, Saxton AM, Yuan JS, Stewart CN: jr Transcriptional responses of Arabidopsis thaliana plants to AsV stress. BMC Plant Biol. 2008, 8: 87-10.1186/1471-2229-8-87.
Ahsan N, Lee D-G, Alam I: Comparative proteomic study of arsenic-induced differentially expressed proteins in rice roots reveals glutathione plays a central role during As stress. Proteomics. 2008, 8: 3561-3576. 10.1002/pmic.200701189.
Romero HM, Berlett BS, Jensen PJ, Pell EJ, Tien M: Investigations into the role of the plastidial peptide methionine sulfoxide reductase in response to oxidative stress in Arabidopsis. Plant Physiol. 2004, 136: 3784-3794. 10.1104/pp.104.046656.
Chang Q, Petrash JM: Disruption of aldo-keto reductase genes leads to elevated markers of oxidative stress and inositol auxotrophy in Saccharomyces cerevisiae. Biochem Biophys Acta. 2008, 1783: 237-245. 10.1016/j.bbamcr.2007.08.008.
Jemiola-Rzeminska M, Rivera C, Suwalsky M, Strzalka K: Interaction of arsenic compounds with model phospholipid membranes. Thermochim Acta. 2007, 458: 132-137. 10.1016/j.tca.2007.01.004.
Dietrich CR, Han G, Chen M, Berg RH, Dunn TM, Cahoon EB: Loss-of-function mutations and inducible RNAi suppression of Arabidopsis LCB2 genes reveal the critical role of sphingolipids in gametophytic and sporophytic cell viability. Plant J. 2008, 54: 284-98. 10.1111/j.1365-313X.2008.03420.x.
Mueller S, Hilbert B, Dueckershoff K, Roitsch T, Krischke M, Mueller MJ, Berger S: General detoxification and stress responses are mediated by oxidized lipids through TGA transcription factors in Arabidopsis. Plant Cell. 2008, 20: 768-785. 10.1105/tpc.107.054809.
Beynon ER, Symons ZC, Jackson RG, Lorenz A, Rylott EL, Bruce NC: The role of oxophytodienoate reductases in the detoxification of the explosive 2, 4, 6- trinitrotoluene by Arabidopsis. Plant Physiol. 2009, 151: 253-261. 10.1104/pp.109.141598.
Wagner U, Edwards R, Dixon DP, Mauch F: Probing the diversity of the Arabidopsis glutathione S-transferase gene family. Plant Mol Biol. 2002, 49: 515-532. 10.1023/A:1015557300450.
Zaman GJR, Lankelma J, van Tellingen O, Beijnen J, Dekker H, Paulusma C, Oude Elferink RPJ, Bass F, Borst P: Role of glutathione in the export of compounds from cells by the multidrug and resistance-associated protein. Proc Nat Acad Sci USA. 1995, 92: 7690-7694. 10.1073/pnas.92.17.7690.
Ghosh M, Shen J, Rosen BP: Pathways of As III detoxification in Saccharomyces cerevisiae. Proc Nat Acad Sci USA. 1999, 96: 5001-5006. 10.1073/pnas.96.9.5001.
Kampranis SC, Damianova R, Atallah M, Toby G, Kondi G, Tsichlis PN, Makris AM: A novel plant glutathione S-transferase/peroxidase suppresses Bax lethality in yeast. J Biol Chem. 2000, 275: 29207-29216. 10.1074/jbc.M002359200.
Noctor G, Foyer CH: Ascorbate and Glutathione: keeping active oxygen under control. Annu Rev Plant Physiol. 1998, 49: 249-279. 10.1146/annurev.arplant.49.1.249.
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. 10.1093/pcp/pcd051.
Leustek T: Sulfate Metabolism. The Arabidopsis Book. American Society of Plant Biologists; Cold Spring Harbor, Cold Spring Harbor Laboratory PressCold Spring Harbor, Cold Spring Harbor Laboratory Press; 2002. 1-16. 10.1199/tab.0017.
Noji M, Saito M, Nakamura M, Aono M, Saji H, Saito K: Cysteine synthase overexpression in tobacco confers tolerance to sulfur-containing environment pollutants. Plant Physiology. 2001, 126: 973-980. 10.1104/pp.126.3.973.
Callis J, Vierstra RD: Protein degradation in signaling. Curr Opin Plant Biol. 2000, 3: 381-386. 10.1016/S1369-5266(00)00100-X.
Li L, He Z, Pandey GK, Tsuchiya T, Luan S: Functional cloning and characterization of a plant efflux carrier for multidrug and heavy metal detoxification. J Biol Chem. 2002, 277: 5360-5368. 10.1074/jbc.M108777200.
Uddin Md I, Qi Y, Yamada S, Shibuya I, Deng X, Kwak S, Kaminaka H, Tanaka K: Overexpression of a new rice vacuolar antiporter regulating protein OsARP improves salt tolerance in tobacco. Plant Cell Physiol. 2008, 49: 880-890. 10.1093/pcp/pcn062.
Zakharyan R, Wu Y, Bogdan GM, Aposhian HV: Enzymatic methylation of arsenic compounds: assay partial purification and properties of arsenite methyltransferase and monomethylarsonic acid methyltransferase from rabbit liver. Chem Res Toxicol. 1995, 8: 1029-1038. 10.1021/tx00050a006.
Hartley-Whitaker J, Ainsworth G, Meharg AA: Copper and arsenate induced oxidative stress in Holcus lanatus L. clones with differential sensitivity. Plant Cell Environ. 2001, 24: 713-722. 10.1046/j.0016-8025.2001.00721.x.
Diatchenko L, Lau YC, Campbell AP, Chenchik A, Moqadam F, Hung B, Lukyanov S, Lukyanov K, Gurskay N, Sverdlov ED, Siebert PD: Suppression subtractive hybridization: a method for generating differentially regulated or tissue-specific cDNA probes and libraries. Proc Nat Acad Sci USA. 1996, 93: 6025-6030. 10.1073/pnas.93.12.6025.
Maniatis T, Fritsch EF, Sambrook J: Molecular cloning: a laboratory manual. Cold Spring Harbor, Cold Spring Harbor Laboratory Press; 1989.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search program. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Misson J, Thibaud MC, Nussaume L, Raghothama K, Jain A, Jouhet J, Block M, Bligny R, Ortet P, Herrerra-Estrella L, Somerville S: A genome wide transcriptional analysis using Arabidopsis thaliana Affymetrix gene chips determined plant responses to phosphate deprivation. Proc Nat Acad Sci USA. 2005, 102: 11934-11939. 10.1073/pnas.0505266102.
Short E, McAinsh M, Shirras A: Array Express ID: E-MEXP-342. Microarray data on Arabidopsis treated with ozone. 2005, [http://www.ebi.ac.uk/microarray-as/aer/#ae-main]
ArrayExpress microarray database: [http://www.ebi.ac.uk/microarray-as/ae/]
Dunnett CW: A multiple comparison procedure for comparing several treatments with a control. J Am Stat Assoc. 1955, 50: 1096-1121. 10.2307/2281208.
Duncan DB: Multiple range tests for correlated and Heteroscedastic means. Biometrics. 1957, 13: 164-176. 10.2307/2527799.
Authors are thankful to Evan Vaine, Drs Rob Wick and John Burand for their critical reading of the manuscript and suggestions. This research was supported from the start-up funds and a Healey Endowment faculty research grant (#21644) to OPD from the University of Massachusetts, Amherst.
BP performed the subtraction procedure, cDNA cloning, sequence analysis, RT-PCR, analyzed data and participated in writing the manuscript. SK grew crambe and Brassica species in arsenic media, performed the arsenic uptake assays and analyzed the data. OPD designed and supervised this whole project and participated in writing the manuscript. All authors read and approved this manuscript.
About this article
Cite this article
Paulose, B., Kandasamy, S. & Dhankher, O.P. Expression profiling of Crambe abyssinicaunder arsenate stress identifies genes and gene networks involved in arsenic metabolism and detoxification. BMC Plant Biol 10, 108 (2010) doi:10.1186/1471-2229-10-108
- Suppression Subtraction Hybridization
- Subtract cDNA Library
- Sulfate Assimilation