Plant Physiol. Illumina
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


First published online June 5, 2009; 10.1104/pp.109.135327

Plant Physiology 150:1972-1980 (2009)
© 2009 American Society of Plant Biologists

OPEN ACCESS ARTICLE
This Article
Free via Open Access: OA
Right arrow OA Abstract
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrowOA All Versions of this Article:
150/4/1972    most recent
pp.109.135327v2
pp.109.135327v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Google Scholar
Right arrow Articles by Maruyama, K.
Right arrow Articles by Yamaguchi-Shinozaki, K.
PubMed
Right arrow PubMed Citation
Right arrow Articles by Maruyama, K.
Right arrow Articles by Yamaguchi-Shinozaki, K.
Agricola
Right arrow Articles by Maruyama, K.
Right arrow Articles by Yamaguchi-Shinozaki, K.
ENVIRONMENTAL STRESS AND ADAPTATION TO STRESS

Metabolic Pathways Involved in Cold Acclimation Identified by Integrated Analysis of Metabolites and Transcripts Regulated by DREB1A and DREB2A1,[W],[OA]

Kyonoshin Maruyama, Migiwa Takeda, Satoshi Kidokoro, Kohji Yamada, Yoh Sakuma, Kaoru Urano, Miki Fujita, Kyouko Yoshiwara, Satoko Matsukura, Yoshihiko Morishita, Ryosuke Sasaki, Hideyuki Suzuki, Kazuki Saito, Daisuke Shibata, Kazuo Shinozaki and Kazuko Yamaguchi-Shinozaki*

Biological Resources Division, Japan International Research Center for Agricultural Sciences, Tsukuba, Ibaraki 305–8686, Japan (K.M., S.K., K.Y., Y.S., K.Y., S.M., K.Y.-S.); Kazusa DNA Research Institute, Kisarazu, Chiba 292–0818, Japan (M.T., Y.M., R.S., H.S., K. Saito, D.S.); Laboratory of Plant Molecular Physiology, Graduate School of Agricultural and Life Sciences, University of Tokyo, Bunkyo-ku, Tokyo 113–8657, Japan (S.K., K.Y., Y.S., S.M., K.Y.-S.); Gene Discovery Research Group, RIKEN Plant Science Center, Tsukuba, Ibaraki 305–0074, Japan (K.U., M.F., K. Shinozaki); and Department of Molecular Biology and Biotechnology, Graduate School of Pharmaceutical Sciences, Chiba University, Chiba 263–8522, Japan (K. Saito)


    ABSTRACT
 TOP
 ABSTRACT
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 LITERATURE CITED
 
DREB1A/CBF3 and DREB2A are transcription factors that specifically interact with a cis-acting dehydration-responsive element (DRE), which is involved in cold- and dehydration-responsive gene expression in Arabidopsis (Arabidopsis thaliana). Overexpression of DREB1A improves stress tolerance to both freezing and dehydration in transgenic plants. In contrast, overexpression of an active form of DREB2A results in significant stress tolerance to dehydration but only slight tolerance to freezing in transgenic plants. The downstream gene products for DREB1A and DREB2A are reported to have similar putative functions, but downstream genes encoding enzymes for carbohydrate metabolism are very different between DREB1A and DREB2A. We demonstrate that under cold and dehydration conditions, the expression of many genes encoding starch-degrading enzymes, sucrose metabolism enzymes, and sugar alcohol synthases changes dynamically; consequently, many kinds of monosaccharides, disaccharides, trisaccharides, and sugar alcohols accumulate in Arabidopsis. We also show that DREB1A overexpression can cause almost the same changes in these metabolic processes and that these changes seem to improve freezing and dehydration stress tolerance in transgenic plants. In contrast, DREB2A overexpression did not increase the level of any of these metabolites in transgenic plants. Strong freezing stress tolerance of the transgenic plants overexpressing DREB1A may depend on accumulation of these metabolites.


Low temperature and dehydration are adverse environmental conditions that affect plant growth and productivity. Many genes have been described that respond to both stresses at the transcriptional level, and their gene products are thought to function in stress tolerance and response even though these stresses are quite different (Thomashow, 1999Go; Zhu, 2002Go; Shinozaki and Yamaguchi-Shinozaki, 2003Go; Bartels and Sunkar, 2005Go). These genes include key metabolic enzymes, late embryogenesis-abundant (LEA) proteins, detoxification enzymes, chaperones, protein kinases, and transcription factors (Thomashow, 1999Go; Zhu, 2002Go; Shinozaki and Yamaguchi-Shinozaki, 2003Go; Bartels and Sunkar, 2005Go). The cis-acting elements that function in stress-responsive gene expression have been analyzed to elucidate the molecular mechanisms of gene expression in response to these stresses (Yamaguchi-Shinozaki and Shinozaki, 2006Go). The dehydration-responsive element (DRE), containing the core sequence A/GCCGAC, was identified as a cis-acting element that regulates gene expression in response to both dehydration and cold in Arabidopsis (Arabidopsis thaliana; Yamaguchi-Shinozaki and Shinozaki, 1994Go). A similar motif was identified as the C-repeat and low temperature-responsive element in promoter regions of low temperature-inducible genes (Baker et al., 1994Go; Jiang et al., 1996Go).

Arabidopsis cDNAs encoding the ethylene-responsive element-binding factor/APETALA2 (ERF/AP2)-type DRE-binding (DREB) proteins CBF1, DREB1A, and DREB2A have been isolated by yeast one-hybrid screening (Stockinger et al., 1997Go; Liu et al., 1998Go). They specifically bind to the DRE/C-repeat sequence and activate transcription of genes driven by the DRE sequence in Arabidopsis. Three DREB1/CBF proteins are encoded by genes that lie in tandem on chromosome 4 in the order DREB1B/CBF1, DREB1A/CBF3, and DREB1C/CBF2 (Gilmour et al., 1998Go; Liu et al., 1998Go). Arabidopsis also has two DREB2 proteins, DREB2A and DREB2B (Liu et al., 1998Go). Although expression of all three DREB1/CBF genes is induced by cold but not by dehydration, both DREB2 genes are induced by dehydration and high salinity (Liu et al., 1998Go). Both DREB1/CBF and DREB2 proteins bind to DRE, but DREB1/CBFs are thought to function in cold-responsive gene expression, whereas DREB2s are involved in dehydration-responsive gene expression.

Overexpression of DREB1/CBFs driven by the cauliflower mosaic virus 35S promoter increases stress tolerance to freezing, dehydration, and high salinity in transgenic Arabidopsis (Jaglo-Ottosen et al., 1998Go; Liu et al., 1998Go; Kasuga et al., 1999Go). More than 40 downstream targets of DREB1A/CBF3 have been identified by microarrays (Seki et al., 2001Go; Fowler and Thomashow, 2002Go; Maruyama et al., 2004Go; Vogel et al., 2005Go). Overexpression of the constitutively active form of DREB2A (35S:DREB2A-CA) significantly increases dehydration tolerance but only slightly increases freezing tolerance. Microarray analyses of the 35S:DREB2A-CA plants revealed that DREB2A regulates the expression of many dehydration-responsive genes. However, some genes regulated by DREB1A are not regulated by DREB2A (Sakuma et al., 2006aGo).

Promoter analysis of DREB1A- and DREB2A-regulated genes and gel mobility shift assays revealed that DREB1A and DREB2A have different DNA-binding specificities. DREB1A has highest affinity for A/GCCGACNT, whereas DREB2A preferentially binds ACCGAC (Maruyama et al., 2004Go; Sakuma et al., 2006aGo). This affinity controls the induction of different downstream genes of DREB1A and DREB2A. Recently, Sakuma et al. (2006b)Go reported that DREB2A expression is also regulated by heat shock (HS) stress and that thermotolerance is significantly increased in 35S:DREB2A-CA plants. HS-inducible genes such as HS proteins and HS transcription factors are also up-regulated in 35S:DREB2A-CA plants but not in 35S:DREB1A plants. Although DREB2A-regulated genes are involved in dehydration, high salinity, and HS stress tolerance, they are not sufficient to withstand freezing stress. DREB1A and DREB2A downstream gene products mostly have similar functions, but downstream genes encoding enzymes for carbohydrate metabolism are different between DREB1A and DREB2A. Several large-scale analyses have investigated metabolites in cold-exposed and 35S:DREB1A plants (Cook et al., 2004Go; Kaplan et al., 2004Go, 2007Go; Hannah et al., 2006Go). The levels of many metabolites (e.g. carbohydrates, amines, and organic acids) increase significantly in both cold-exposed and 35S:DREB1A plants. However, there has been no reported effort to perform large-scale metabolite analyses in dehydration-exposed and 35S:DREB2A-CA plants.

Here, we performed integrated analysis of both metabolites and transcripts in dehydration-exposed and 35S:DREB2A-CA plants and compared them with those in cold-exposed and 35S:DREB1A plants. We also analyzed the expression of genes involved in carbohydrate metabolism under both cold and dehydration conditions.


    RESULTS
 TOP
 ABSTRACT
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 LITERATURE CITED
 

Identification of Metabolites in Various Plants

We measured metabolites using gas chromatography-time of flight-mass spectrometry (GC-TOF-MS), capillary electrophoresis-mass spectrometry (CE-MS), and liquid chromatography-ion trap-mass spectrometry (LC-IT-MS) in eight kinds of plants: untreated, cold exposed (1 and 4 d), dehydration exposed (2 and 3 d), 35S:DREB1A, 35S:DREB2A-CA, and control plants containing pGreen. In GC-TOF-MS and CE-MS analyses, 254 independent metabolites were identified by retention time indices and specific mass fragments (Supplemental Fig. S1). The levels of 155, 162, 37, and 28 metabolites increased significantly in cold-exposed, dehydration-exposed, 35S:DREB1A, and 35S:DREB2A-CA plants, respectively (Benjamini and Hochberg false discovery rate [FDR]; P < 0.01). The levels of 50, 85, nine, and 38 metabolites decreased significantly in cold-exposed, dehydration-exposed, 35S:DREB1A, and 35S:DREB2A-CA plants, respectively (FDR; P < 0.01; Supplemental Table S1).

Venn diagrams illustrate the identified metabolites in the plants (Fig. 1 ; Supplemental Fig. S2). The levels of 102 metabolites increased in both cold-exposed and dehydration-exposed plants. The levels of 17 metabolites increased in both 35S:DREB1A and 35S:DREB2A-CA plants. Most metabolites (89%) that increased in 35S:DREB1A plants also increased in cold-exposed plants. These metabolites contained sugars such as Suc, galactinol, myoinositol, raffinose, and unknown metabolites. Likewise, most metabolites (89%) that increased in 35S:DREB2A-CA plants also increased in dehydration-exposed plants. In this case, these metabolites contained organic acids such as argininosuccinate, fumarate, malic acid, and unknown metabolites. In addition, most metabolites that were almost unchanged in 35S:DREB1A plants were also unchanged in 35S:DREB2A-CA plants; unfortunately, these metabolites are unknown (Fig. 1).


Figure 1
View larger version (49K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 1. Venn diagrams of identified metabolites that are increased relative to the controls. The diagrams illustrate the number of identified metabolites in six kinds of plants: cold exposed (1 and 4 d), dehydration exposed (2 and 3 d), 35S:DREB1A, and 35S:DREB2A-CA.

 
LC-IT-MS analyses showed that accumulation of three secondary metabolites increased significantly in the cold-exposed, 35S:DREB1A, and 35S:DREB2A-CA plants but not in the dehydration-exposed plants (Supplemental Fig. S3). Mass fragment/retention time comparisons revealed these three secondary metabolites as kaempferol 3-[6'-(glucosyl) rhamnoside] 7-rhamnoside, kaempferol 3-glucoside 7-rhamnoside, and kaempferol 3-rhamnoside 7-rhamnoside.


Seventeen Metabolites Were Unique to Both Cold-Exposed and 35S:DREB1A Plants

We compared the metabolite profiles of the eight analyzed plants by principal component analysis (PCA; Fig. 2A ; Supplemental Tables S2–S6). The cumulative contribution ratio of the PCA reached 61.4% up to the second principal component (PC2), indicating correct application of PCA to the metabolite profiles of the analyzed plants. The PCA showed the following. First, the first principal component (PC1) reflected increases of metabolites in the eight kinds of plants. The PC1 value of plants exposed to 3 d of dehydration was the highest in the eight kinds of plants. Plants exposed to 4 d of cold had the second highest PC1 value. The PC1 value of untreated plants was similar to that of control plants, and that value was the lowest in the eight kinds of plants. Second, PC2 reflected differences in the variety of metabolites according to different environmental conditions. The PC2 value of cold-exposed plants was positive, but PC2 for dehydration-exposed plants was negative. The PC2 values of both untreated and control plants were nearly zero. Third, the metabolite profile of the 35S:DREB1A plants was similar to that of cold-exposed plants, and the PC2 values of the plants were positive. Fourth, the metabolite profile of 35S:DREB2A-CA plants was similar to that of dehydration-exposed plants, and PC2 values were negative.


Figure 2
View larger version (29K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 2. Statistical analyses of metabolite profiles. We analyzed two independent lines of 35S:DREB1A ({alpha} and β) and 35S:DREB2A-CA ({alpha} and β) plants. The levels of metabolites for both DREB1A and DREB2A in each β line were higher than those in each {alpha} line. A, PCA of metabolites. The y and x axes represent PC1 and PC2, respectively. The circles indicate untreated, cold-exposed, and dehydration-exposed plants. The diamonds represent control, 35S:DREB1A, and 35S:DREB2A-CA plants. B, Representative metabolites for which the eigenvector values were the first and second highest or lowest. In each case, the maximum level of the metabolite was set to 100. Error bars indicate SD for three experiments. A star indicates that the metabolite was not detected. Metabolites in the {alpha} and β lines of each transgenic plant are shown by the left and right bars, respectively. C, Selected metabolites that were increased in both cold-exposed and 35S:DREB1A plants but minimal in 35S:DREB2A-CA plants.

 
We selected representative metabolites for which the eigenvector values were the first and second highest or lowest, and those values are displayed in bar charts (Fig. 2B). The PC1 eigenvector value of metabolite 190 was the highest. The level of metabolite 190 increased significantly in plants exposed to 3 d of dehydration or 4 d of cold. This metabolite was also detected in 35S:DREB1A plants but not in 35S:DREB2A-CA plants. Raffinose had the second highest PC1 eigenvector value. The level of raffinose increased significantly in plants exposed to 3 d of dehydration or 4 d of cold and in 35S:DREB1A plants. The level of raffinose in the 35S:DREB1A plants was the highest of those in the eight kinds of plants. We did not detect raffinose in either untreated or control plants. Raffinose was detectable in 35S:DREB2A-CA plants, but its level was very low compared with 35S:DREB1A plants. Metabolite 70 had the highest PC2 eigenvector value. The level of metabolite 70 increased significantly in both cold-exposed and 35S:DREB1A plants but not in dehydration-exposed or 35S:DREB2A-CA plants. Metabolite 73 had the second highest PC2 eigenvector value. The level of metabolite 73 increased significantly in cold-exposed, 35S:DREB1A, and 35S:DREB2A-CA plants. Metabolite 72 had the lowest PC2 eigenvector value. Metabolite 72 increased significantly in plants exposed to 3 d of dehydration and decreased significantly in both cold-exposed and 35S:DREB1A plants. Metabolite 4 had the second lowest PC2 eigenvector. The level of metabolite 4 increased significantly in plants exposed to 3 d of dehydration.

To identify candidates that play vital roles under low-temperature conditions, we tried to select metabolites that were unique to both cold-exposed and 35S:DREB1A plants but were not affected in 35S:DREB2A-CA plants. Seventeen such metabolites were selected, and their values are displayed in bar charts (Fig. 2, B and C). Among the 17 metabolites, 13 are unknown and four are myoinositol (metabolite 247), galactinol (metabolite 239), raffinose (metabolite 251), and Suc (metabolite 253). The levels of metabolites 55, 70, 92, 155, 156, and 189 increased significantly in both cold-exposed and 35S:DREB1A plants but not in dehydration-exposed plants. The levels of raffinose, galactinol, myoinositol, Suc, and metabolites 19, 80, 85, 100, 193, 200, and 209 increased significantly in not only cold-exposed and 35S:DREB1A plants but also in dehydration-exposed plants. Among the 17 metabolites, raffinose, galactinol, myoinositol, Suc, and metabolite 70 were detectable in 35S:DREB2A-CA plants at very low levels; however, the other 12 metabolites were undetectable in 35S:DREB2A-CA plants.


Identification of Cold- or Dehydration-Responsive Genes and Downstream Genes of DREB1A and DREB2A

To identify cold- or dehydration-responsive genes and downstream genes of DREB1A and DREB2A, we performed array analyses using an Agilent Arabidopsis 2 Oligo Microarray. In the cold-exposed, dehydration-exposed, 35S:DREB1A, and 35S:DREB2A-CA plants, 945, 2,641, 259, and 373 genes were significantly up-regulated, respectively (FDR; P < 0.05 and fold change [FC] > 2; Supplemental Tables S7–S12). In contrast, 732, 2,999, 174, and 234 genes were down-regulated significantly in the cold-exposed, dehydration-exposed, 35S:DREB1A, and 35S:DREB2A-CA plants, respectively (FDR; P < 0.05 and FC < 1/2; Supplemental Tables S7–S12). Of the 259 up-regulated genes in 35S:DREB1A plants, 125 were unknown and the other gene products included dehydrin, LEA protein, starch-degrading enzymes, transcription factors, and protein kinases. Of the 373 up-regulated genes in 35S:DREB2A-CA plants, 158 were unknown and the other gene products included dehydrin, LEA protein, HS protein family, enzymes for toxin catabolic processes, transcription factors, and protein kinases (Fig. 3 ).


Figure 3
View larger version (46K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 3. Functional categorization of DREB1A and DREB2A-CA downstream genes. Shown are 20 functional categories of DREB1A and DREB2A-CA downstream genes.

 

Expression Patterns of Genes for Starch-Degrading Enzymes and Sugar Alcohol Synthases in 35S:DREB1A Plants Were Similar to Those in Plants Exposed to Cold for 4 d

We reconstructed a metabolic enzyme database for Arabidopsis. First, we downloaded all of the amino acid sequence data from The Arabidopsis Information Resource and searched enzyme motifs using the HMMER program (Eddy, 1998Go). Then, we downloaded data for enzymes involved in metabolism from KEGG, AraCyc, and Kappa-view and performed similarity searches of the downloaded data against all Arabidopsis proteins using the PSI-BLAST program. Based on these results, we selected 2,232 metabolic enzymes from all Arabidopsis proteins (Supplemental Table S13) as a database. Among the genes for the selected 2,232 enzymes, 2,170 were detectable in the microarray (Supplemental Table S13). Of the 2,170 genes, 115, 385, 42, and 44 were significantly up-regulated in cold-exposed, dehydration-exposed, 35S:DREB1A, and 35S:DREB2A-CA plants, respectively (FDR; P < 0.05 and FC > 2; Supplemental Fig. S4; Supplemental Tables S14–S19). In contrast, 143, 414, 35, and 50 genes were down-regulated significantly in cold-exposed, dehydration-exposed, 35S:DREB1A, and 35S:DREB2A-CA plants, respectively (FDR; P < 0.05 and FC < 1/2; Supplemental Fig. S4; Supplemental Tables S14–S19).

Expression levels of genes encoding starch-degrading enzymes, Suc metabolism enzymes, and sugar alcohol synthases correlated positively with the accumulation of monosaccharides, disaccharides, trisaccharides, and sugar alcohols in cold- and dehydration-exposed plants (Fig. 4 ). Genes encoding starch-degrading enzymes, such as {alpha}-amylase, β-amylase, and glucan-water dikinase, were up-regulated significantly under both cold and dehydration conditions. The downstream products of these enzymes, such as Glc, Fru, and Suc, accumulated under both cold and dehydration conditions. In addition, genes for galactinol synthase were up-regulated significantly under both cold and dehydration conditions. The levels of both galactinol and raffinose increased under both cold and dehydration conditions. On the other hand, these sugars except Glc and Fru increased in the 35S:DREB1A plants, but none of them increased in the 35S:DREB2A-CA plants as compared with control plants (Fig. 4).


Figure 4
View larger version (38K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 4. Map of starch degradation and Suc metabolism pathways. Each small square indicates the expression level of the gene that shows the highest expression in each gene family. Red squares show more than 4-fold increase relative to the control. Orange squares show between 2- and 4-fold increase. Octagons indicate metabolite accumulations. Each pink section indicates increased metabolites in each plant. We could not measure starch, maltooligosaccharides, Glc-1-P, Fru-6-P, Suc-6-P, UDP-Glc, and UDP-Gal. C, D, 1A, and 2A indicate cold-exposed, dehydration-exposed, 35S:DREB1A, and 35S:DREB2A-CA plants, respectively. We analyzed two independent lines of 35S:DREB1A ({alpha} and β) and 35S:DREB2A-CA ({alpha} and β) plants. GWD, Glucan-water dikinase.

 
We then focused on genes for starch-degrading enzymes and sugar alcohol synthase (Fig. 5A ). We confirmed expression levels of these genes using quantitative reverse transcription (qRT)-PCR (Fig. 5B). The genes for starch-degrading enzymes were divided into two groups: one included cold-inducible genes, and the other included dehydration-inducible genes. Expression of AMY3 and BAM3/BMY8 increased specifically under cold conditions. In contrast, expression of AMY2, BAM1/BMY7, BAM2/BMY9, BAM6, and BAM9/BMY3 increased specifically under dehydration conditions. Furthermore, most of these genes were regulated oppositely under dehydration and cold conditions (Fig. 5). Whereas expression of AMY3, BAM3/BMY8, PHS1/PHO1, PHS2/PHO2, ISA3, DPE1, DPE2, and GWD1/SEX1 increased under cold conditions, these transcripts decreased under dehydration conditions. In addition, expression of BAM1/BMY7, BAM6, and BAM9/BMY3 increased under dehydration conditions but decreased under cold conditions. Within the galactinol synthase gene family, expression of AtGolS3 and AtGolS2 increased significantly in cold- and dehydration-exposed plants, respectively (Fig. 5). Expression patterns of genes for starch-degrading enzymes and sugar alcohol synthase in 35S:DREB1A plants were similar to those in plants exposed to cold for 4 d (Fig. 5). However, expression of these genes did not increase in 35S:DREB2A-CA plants, except for AtGolS3 and AtGolS2 (Fig. 5).


Figure 5
View larger version (33K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 5. Expression of genes for starch degradation- and Suc metabolism-related enzymes. A, Each small square indicates the level of gene expression. Red squares show more than 4-fold increase relative to the control. Orange squares show between 2- and 4-fold increase. Light blue squares show between 1/2- and 1/4-fold increase. Dark blue squares show less than 1/4-fold increase. C1, C4, D2, D3, 1A, and 2A indicate cold-exposed (1 d), cold-exposed (4 d), dehydration-exposed (2 d), dehydration-exposed (3 d), 35S:DREB1A, and 35S:DREB2A-CA plants, respectively. We analyzed two independent lines of 35S:DREB1A ({alpha} and β) and 35S:DREB2A-CA ({alpha} and β) plants. The levels of BAM3, BAM1, AtGolS3, AtGolS2, DIN10, and SIP transcripts were detected by qRT-PCR. B, The levels of transcripts for genes encoding β-amylase, galactinol, and raffinose synthase determined by qRT-PCR. Panels 1 and 2 show the levels of transcripts for BAM3 and BAM1 encoding β-amylase, respectively. Panels 3 and 4 show the levels of transcripts for AtGolS3 and AtGolS2 encoding galactinol synthase, respectively. Panels 5 and 6 show the levels of transcripts for DIN10 and SIP1 encoding raffinose synthase, respectively. Error bars indicate SD for three experiments. A star indicates that the transcript was not detected.

 

    DISCUSSION
 TOP
 ABSTRACT
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 LITERATURE CITED
 
Previously, we showed that overexpression of DREB1A improved both dehydration and freezing stress tolerance, whereas overexpression of DREB2A-CA improved dehydration stress tolerance but not freezing tolerance in transgenic plants (Sakuma et al., 2006aGo). In this study, we used microarrays to analyze genes that are regulated by DREB1A and DREB2A under the same conditions as used for metabolite analyses and found that 259 and 373 genes were significantly increased in 35S:DREB1A and 35S:DREB2A-CA plants, respectively, in comparison with vector control plants (Fig. 1). The putative functions of most of the DREB1A- and DREB2A-targeted genes were similar, but certain genes were specific to DREB2A, such as molecular chaperones and enzymes for toxin catabolic processes. We also found that most putative functions of DREB1A downstream gene products were included in those of the DREB2A downstream gene products, but the number of DREB1A downstream genes encoding enzymes for carbohydrate metabolism was much higher than that of DREB2A downstream genes. Thus, expression of these DREB1A-specific downstream genes encoding enzymes for carbohydrate metabolism may be involved in freezing tolerance in 35S:DREB1A plants.

We aimed to detect many kinds of metabolites that are regulated by abiotic stress using GC-TOF-MS, LC-IT-MS, and CE-MS. GC-TOF-MS and LC-IT-MS are suitable for analysis of polar metabolites and secondary metabolites, respectively. CE-MS was used to analyze charged metabolites. We identified 254 independent metabolites using GC-TOF-MS and CE-MS. Among them, 61.0%, 63.7%, 14.5%, and 11.0% of metabolites increased significantly in cold-exposed, dehydration-exposed, 35S:DREB1A, and 35S:DREB2A-CA plants, respectively (Fig. 1). Detected metabolites were compared by PCA (Fig. 2A). The metabolite profiles were classified into three groups according to plant growth conditions: cold, dehydration, and untreated/control. The metabolite profile of the 35S:DREB1A plants was similar to that of the cold class but not to that of the dehydration class. The metabolite profile of the 35S:DREB2A-CA plants resembled that of the dehydration class more than that of the cold class. Because 35S:DREB1A plants, but not 35S:DREB2A-CA plants, were freezing tolerant, the metabolites for which the level increased in the cold-exposed and 35S:DREB1A plants but not in the 35S:DREB2A-CA plants probably participate in improvement of freezing stress tolerance. Seventeen metabolites were selected as candidates. Among them, 13 are unknown and four are myoinositol, Suc, galactinol, and raffinose. Although both transgenic plants showed strong dehydration tolerance, the metabolite profiles in these plants were very different. These results suggest that the dehydration-specific metabolites appearing on the negative side of the PC2 axis may not be important for dehydration-specific tolerance in 35S:DREB1A plants. Several unknown metabolites that increased under both cold and dehydration conditions and also in both transgenic plants might be involved in dehydration tolerance in the transgenic plants (Supplemental Table S5). We determined using LC-IT-MS that accumulation of only three secondary metabolites increased significantly in the cold-exposed and 35S:DREB1A plants (Supplemental Fig. S3). These secondary metabolites also increased in the 35S:DREB2A-CA plants but not in the dehydration-exposed plants. Korn et al. (2008)Go also reported that accumulation of several secondary metabolites similar to these metabolites increased in cold-treated Arabidopsis plants. However, our comparative analyses suggest that these secondary metabolites may not be important for freezing tolerance in 35S:DREB1A plants, because they also increased significantly in the 35S:DREB2A-CA plants.

Our integrated analysis of metabolites and transcripts indicated that expression of many genes encoding starch-degrading enzymes, Suc metabolism enzymes, and sugar alcohol synthases changed dynamically and resulted in the accumulation of monosaccharides, disaccharides, trisaccharides, and sugar alcohols, including Suc, myoinositol, galactinol, and raffinose, in cold-exposed and dehydration-exposed plants. Because most of these genes are targets of DREB1A, expression of these genes also changed dynamically and the metabolites also accumulated in the 35S:DREB1A plants. In contrast, because these genes are not regulated by DREB2A (except the genes for galactinol synthase and raffinose synthase), none of these metabolites accumulated in 35S:DREB2A-CA plants (Fig. 4). Expression of these genes and accumulation of these metabolites correlated with differences in freezing tolerance between the 35S:DREB1A and 35S:DREB2A-CA plants. These results suggest that dynamic transcriptional regulation of the carbohydrate network is necessary for the accumulation of specific carbohydrates, such as Suc, galactinol, myoinositol, and raffinose, and that the accumulation of these carbohydrates may be important for improving freezing tolerance in transgenic plants overexpressing DREB1A. Moreover, it is likely that DREB1A is one of the major factors regulating carbohydrate metabolism under cold conditions, whereas DREB2A is not.

A cold-inducible galactinol synthase, AtGolS3, plays a key role in the accumulation of galactinol and raffinose under cold conditions (Taji et al., 2002Go; Cook et al., 2004Go). Although the level of the AtGolS3 transcript in the 35S:DREB2A-CA plants was similar to that in cold-exposed and 35S:DREB1A plants, galactinol and raffinose content in the 35S:DREB2A-CA plants did not increase compared with cold-exposed and 35S:DREB1A plants (Figs. 4 and 5). Dynamic changes in gene expression may be necessary for plants to accumulate these metabolites, and strong freezing stress tolerance of 35S:DREB1A plants may depend on the accumulation of these metabolites. It has been reported that an increase of the raffinose content in transgenic petunia (Petunia hybrida) plants improves freezing tolerance (Pennycooke et al., 2003Go). By contrast, using transgenic Arabidopsis, Zuther et al. (2004)Go determined that raffinose is not an essential component of basic freezing tolerance and cold acclimation. Future studies may be necessary to establish the importance of raffinose in the freezing tolerance of Arabidopsis plants. Recently, Nishizawa et al. (2008)Go reported that galactinol and raffinose act as scavengers of the hydroxyl radical in vitro. Under low-temperature conditions, raffinose probably acts not only as an osmoprotectant and a stabilizer of cellular membranes but also as a scavenger of reactive oxygen species to protect the photosynthesis complex in chloroplasts of the cold-exposed and 35S:DREB1A plants.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 LITERATURE CITED
 

Plant Materials and Growth Conditions

Both wild-type and transgenic Arabidopsis (Arabidopsis thaliana) plants were grown in plastic pots filled with peat moss for 3 weeks (principal growth stage 1.07–1.08) under a 16-h-light/8-h-dark regimen (40 ± 10 µmol photons m–2 s–1) at 22°C. The light intensity was the same at 22°C and 4°C. Plants were harvested 2 h later after lights were turned on. For cold treatment, 3-week-old plants were transferred from 22°C to 4°C and were grown for 1 or 4 d. For dehydration treatment, 3-week-old plants were grown for 2 or 3 d without watering. To obtain accurate results, we carefully raised single plants in petri dishes, each containing an equal amount of soil. Soil moisture contents were calculated from soil dry weight. Untreated, the soil moisture content was 84.3%. Under dehydration, on the 2nd d, the soil moisture content was 51.1%. Under dehydration, on the 3rd d, the moisture content was 11.6% (Supplemental Fig. S5).


Metabolite Extraction and Detection

GC-TOF-MS Analysis
Metabolites were extracted from aerial parts of Arabidopsis plants (50 mg) with methanol. Extraction and derivatization were performed as described (Kusano et al., 2007Go). Prepared metabolites were detected by an Agilent 6890 GC apparatus (Agilent Technologies) coupled to a TOF-MS device (Leco) on a 30-m DB-17ms column (J&W Scientific; 0.25 mm i.d., 0.25-mm film). Ribitol was used as an internal standard (Kusano et al., 2007Go). The reproducibility of GC-TOF-MS analysis was determined using three biological replicates in each experiment.

CE-MS Analysis
Metabolites were extracted from aerial parts of Arabidopsis plants (50 mg) with chloroform/methanol solution using a mixer and a 5-kD cutoff filter (Millipore). Prepared metabolites were detected by CE-MS (Agilent 1100 series MSD mass spectrometer, Agilent 1100 series isocratic HPLC pump, Agilent G1603A CE-MS adapter kit, and Agilent G1607A CE-MS sprayer kit). Met sulfone was used as an internal standard. The reproducibility of CE-MS analysis was determined using three biological replicates in each experiment. Unknown metabolites were measured by CE-MS.

LC-IT-MS Analysis
Metabolites were extracted from aerial parts of Arabidopsis plants (50 mg) with methanol:acetic acid:water (9:1:10) as described (Tohge et al., 2005Go). Prepared metabolites were detected by an Agilent HPLC 1100 series coupled to a Finnigan LTQ mass spectrometer (Thermo Electron). Kaempferol 3-O-rutinoside was used as a reference compound for the standard curve. The reproducibility of LC-IT-MS analysis was determined using five biological replicates in each experiment.

Microarray and qRT-PCR Analyses
Total RNA was isolated from aerial parts of Arabidopsis plants with Trizol reagent (Invitrogen) and used for the preparation of Cy5- and Cy3-labeled complementary RNA probes. All microarray experiments, including the data analysis, were performed according to Agilent methods. The reproducibility of microarray analysis was assayed by biological and technical (dye-swap) replicates in each experiment (Invitrogen), and qRT-PCR was performed as described (Sakuma et al., 2006aGo). Microarray analysis and data mining were carried out as described (Qin et al., 2008Go). Array design and data from this article have been deposited at MIAMExpress under accession numbers E-MEXP-2173, E-MEXP-2174, and E-MEXP-2175.


Functional Categorization of Genes

Cold- or dehydration-responsive genes and downstream genes of DREB1A and DREB2A were annotated according to results of motif (HMMER) and similarity (PSI-BLAST) searches and were classified according to Gene Ontology annotation at The Arabidopsis Information Resource, KEGG, and Kappa-view.


Supplemental Data

The following materials are available in the online version of this article.

Supplemental Figure S1. Venn diagrams of identified metabolites.
Supplemental Figure S2. LC-IT-MS analysis.
Supplemental Figure S3. Functional categorization of DREB1A and DREB2A-CA downstream genes.
Supplemental Figure S4. Venn diagrams of metabolite-related gene expression.
Supplemental Figure S5. Dehydration treatment.
Supplemental Table S1. Mean and SD.
Supplemental Table S2. Sum mean SD.
Supplemental Table S3. Variance covariance matrix.
Supplemental Table S4. Eigenvalue.
Supplemental Table S5. Eigenvector.
Supplemental Table S6. PC score.
Supplemental Table S7. Cold-exposed plants (1 d), all.
Supplemental Table S8. Cold-exposed plants (4 d), all.
Supplemental Table S9. Dehydration-exposed plants (2 d), all.
Supplemental Table S10. Dehydration-exposed plants (3 d), all.
Supplemental Table S11. 35S:DREB1A transgenic plants, all.
Supplemental Table S12. 35S:DREB2A-CA transgenic plants, all.
Supplemental Table S13. List of metabolic enzymes.
Supplemental Table S14. Cold-exposed plants (1 d), metabolic enzymes.
Supplemental Table S15. Cold-exposed plants (4 d), metabolic enzymes.
Supplemental Table S16. Dehydration-exposed plants (2 d), metabolic enzymes.
Supplemental Table S17. Dehydration-exposed plants (3 d), metabolic enzymes.
Supplemental Table S18. 35S:DREB1A transgenic plants, metabolic enzymes.
Supplemental Table S19. 35S:DREB2A-CA transgenic plants, metabolic enzymes.


    ACKNOWLEDGMENTS
 
We are grateful to E. Ohgawara, K. Murai, E. Kishi, K. Amano, and H. Sado for their excellent technical support and to M. Toyoshima for skillful editorial assistance. We also thank the Rice Genome Resource Center at the National Institute of Agrobiological Sciences for use of the 22 K microarray analysis system and the technical support provided by Dr. Y. Nagamura and R. Motoyama.

Received January 5, 2009; accepted June 3, 2009; published June 5, 2009.


    FOOTNOTES
 
1 This work was supported by project grants from the Ministry of Agriculture, Forestry, and Fisheries of Japan, grants-in-aid from the Ministry of Education, Culture, Sports, Science, and Technology of Japan, and project grants from the New Energy and Industrial Technology Development Organization of Japan. Back

The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: Kazuko Yamaguchi-Shinozaki (kazukoys{at}jircas.affrc.go.jp).

[W] The online version of this article contains Web-only data. Back

[OA] Open Access articles can be viewed online without a subscription. Back

www.plantphysiol.org/cgi/doi/10.1104/pp.109.135327

* Corresponding author; e-mail kazukoys{at}jircas.affrc.go.jp.


    LITERATURE CITED
 TOP
 ABSTRACT
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 LITERATURE CITED
 
Baker SS, Wilhelm KS, Thomashow MF (1994) The 5'-region of Arabidopsis thaliana cor15a has cis-acting elements that confer cold-, drought- and ABA-regulated gene expression. Plant Mol Biol 24: 701–713[CrossRef][Web of Science][Medline]

Bartels D, Sunkar R (2005) Drought and salt tolerance in plants. Crit Rev Plant Sci 24: 23–58[CrossRef]

Cook D, Fowler S, Fiehn O, Thomashow MF (2004) A prominent role of the CBF cold response pathway in configuring the low temperature metabolome of Arabidopsis. Proc Natl Acad Sci USA 101: 15243–15248[Abstract/Free Full Text]

Eddy SR (1998) Profile hidden Markov models. Bioinformatics 14: 755–763[Abstract/Free Full Text]

Fowler S, Thomashow MF (2002) Arabidopsis transcriptome profiling indicates multiple regulatory pathways are activated during cold acclimation in addition to the CBF cold-response pathway. Plant Cell 14: 1675–1690[Abstract/Free Full Text]

Gilmour SJ, Zarka DG, Stockinger EJ, Salazar MP, Houghton JM, Thomashow MF (1998) Low temperature regulation of the Arabidopsis CBF family of AP2 transcriptional activators as an early step in cold-induced COR gene expression. Plant J 16: 433–442[CrossRef][Web of Science][Medline]

Hannah MA, Wiese D, Freund S, Fiehn O, Heyer AG, Hincha DK (2006) Natural genetic variation of freezing tolerance in Arabidopsis. Plant Physiol 142: 98–122[Abstract/Free Full Text]

Jaglo-Ottosen KR, Gilmour SJ, Zarka DG, Schabenberger O, Thomashow MF (1998) Arabidopsis CBF1 overexpression induces COR genes and enhances freezing tolerance. Science 280: 104–106[Abstract/Free Full Text]

Jiang C, Iu B, Singh J (1996) Requirement of a CCGAC cis-acting element for cold induction of the BN115 gene from winter Brassica napus. Plant Mol Biol 30: 679–684[CrossRef][Web of Science][Medline]

Kaplan F, Kopka J, Haskell DW, Zhao W, Schiller KC, Gatzke N, Sung DY, Guy CL (2004) Exploring the temperature stress metabolome of Arabidopsis. Plant Physiol 136: 4159–4168[Abstract/Free Full Text]

Kaplan F, Kopka J, Sung DY, Zhao W, Popp M, Porat R, Guy CL (2007) Transcript and metabolite profiling during cold acclimation of Arabidopsis reveals an intricate relationship of cold-regulated gene expression with modifications in metabolite content. Plant J 50: 967–981[CrossRef][Web of Science][Medline]

Kasuga M, Liu Q, Miura S, Yamaguchi-Shinozaki K, Shinozaki K (1999) Improving plant drought, salt, and freezing tolerance by gene transfer of a single stress-inducible transcription factor. Nat Biotechnol 17: 287–291[CrossRef][Web of Science][Medline]

Korn M, Peterek S, Mock HP, Heyer AG, Hincha DK (2008) Heterosis in the freezing tolerance, and sugar and flavonoid contents of crosses between Arabidopsis thaliana accessions of widely varying freezing tolerance. Plant Cell Environ 31: 813–827[Medline]

Kusano M, Fukushima A, Arita M, Jonsson P, Moritz T, Kobayashi M, Hayashi N, Tohge T, Saito K (2007) Unbiased characterization of genotype-dependent metabolic regulations by metabolomic approach in Arabidopsis. BMC Syst Biol 21: 1–53

Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-Shinozaki K, Shinozaki K (1998) Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-responsive gene expression, respectively, in Arabidopsis. Plant Cell 10: 1391–1406[Abstract/Free Full Text]

Maruyama K, Sakuma Y, Kasuga M, Ito Y, Seki M, Goda H, Shimada Y, Yoshida S, Shinozaki K, Yamaguchi-Shinozaki K (2004) Identification of cold-inducible downstream genes of the Arabidopsis DREB1A/CBF3 transcriptional factor using two microarray systems. Plant J 38: 982–993[CrossRef][Web of Science][Medline]

Nishizawa A, Yabuta Y, Shigeoka S (2008) Galactinol and raffinose constitute a novel function to protect plants from oxidative damage. Plant Physiol 147: 1251–1263[Abstract/Free Full Text]

Pennycooke JC, Jones ML, Stushnoff C (2003) Down-regulating {alpha}-galactosidase enhances freezing tolerance in transgenic petunia. Plant Physiol 133: 901–909[Abstract/Free Full Text]

Qin F, Sakuma Y, Tran LS, Maruyama K, Kidokoro S, Fujita Y, Fujita M, Umezawa T, Sawano Y, Miyazono K, et al (2008) Arabidopsis DREB2A-interacting proteins function as RING E3 ligases and negatively regulate plant drought stress-responsive gene expression. Plant Cell 20: 1693–1707[Abstract/Free Full Text]

Sakuma Y, Maruyama K, Osakabe Y, Qin F, Seki M, Shinozaki K, Yamaguchi-Shinozaki K (2006a) Functional analysis of an Arabidopsis transcription factor, DREB2A, involved in drought-responsive gene expression. Plant Cell 18: 1292–1309[Abstract/Free Full Text]

Sakuma Y, Maruyama K, Qin F, Osakabe Y, Shinozaki K, Yamaguchi-Shinozaki K (2006b) Dual function of an Arabidopsis transcription factor DREB2A in water-stress-responsive and heat-stress-responsive gene expression. Proc Natl Acad Sci USA 103: 18822–18827[Abstract/Free Full Text]

Seki M, Narusaka M, Abe H, Kasuga M, Yamaguchi-Shinozaki K, Carninci P, Hayashizaki Y, Shinozaki K (2001) Monitoring the expression pattern of 1300 Arabidopsis genes under drought and cold stresses by using a full-length cDNA microarray. Plant Cell 13: 61–72[Abstract/Free Full Text]

Shinozaki K, Yamaguchi-Shinozaki K (2003) Regulatory network of gene expression in the drought and cold stress responses. Curr Opin Plant Biol 6: 410–417[CrossRef][Web of Science][Medline]

Stockinger EJ, Gilmour SJ, Thomashow MF (1997) Arabidopsis CBF1 encodes an AP2 domain-containing transcriptional activator that binds to the C-repeat/DRE, a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and water deficit. Proc Natl Acad Sci USA 94: 1035–1040[Abstract/Free Full Text]

Taji T, Ohsumi C, Iuchi S, Seki M, Kasuga M, Kobayashi M, Yamaguchi-Shinozaki K, Shinozaki K (2002) Important roles of drought- and cold-inducible genes for galactinol synthase in stress tolerance in Arabidopsis. Plant J 29: 417–426[CrossRef][Web of Science][Medline]

Thomashow MF (1999) Plant cold acclimation: freezing tolerance genes and regulatory mechanisms. Annu Rev Plant Physiol Plant Mol Biol 50: 571–599[CrossRef][Web of Science]

Tohge T, Nishiyama Y, Hirai MY, Yano M, Nakajima J, Awazuhara M, Inoue E, Takahashi H, Goodenowe DB, Kitayama M, et al (2005) Functional genomics by integrated analysis of metabolome and transcriptome of Arabidopsis plants over-expressing an MYB transcription factor. Plant J 42: 218–235[CrossRef][Web of Science][Medline]

Vogel JT, Zarka DG, Van Buskirk HA, Fowler SG, Thomashow MF (2005) Roles of the CBF2 and ZAT12 transcription factors in configuring the low temperature transcriptome of Arabidopsis. Plant J 41: 195–211[CrossRef][Web of Science][Medline]

Yamaguchi-Shinozaki K, Shinozaki K (1994) A novel cis-acting element in an Arabidopsis gene is involved in responsiveness to drought, low-temperature, or high-salt stress. Plant Cell 6: 251–264[Abstract]

Yamaguchi-Shinozaki K, Shinozaki K (2006) Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. Annu Rev Plant Biol 57: 781–803[CrossRef][Medline]

Zhu JK (2002) Salt and drought stress signal transduction in plants. Annu Rev Plant Biol 53: 247–273[CrossRef][Medline]

Zuther E, Büchel K, Hundertmark M, Stitt M, Hincha DK, Heyer AG (2004) The role of raffinose in the cold acclimation response of Arabidopsis. FEBS Lett 576: 169–173[CrossRef][Web of Science][Medline]





This Article
Free via Open Access: OA
Right arrow OA Abstract
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrowOA All Versions of this Article:
150/4/1972    most recent
pp.109.135327v2
pp.109.135327v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Google Scholar
Right arrow Articles by Maruyama, K.
Right arrow Articles by Yamaguchi-Shinozaki, K.
PubMed
Right arrow PubMed Citation
Right arrow Articles by Maruyama, K.
Right arrow Articles by Yamaguchi-Shinozaki, K.
Agricola
Right arrow Articles by Maruyama, K.
Right arrow Articles by Yamaguchi-Shinozaki, K.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
ASPB Publications PLANT PHYSIOLOGY® THE PLANT CELL
Copyright © 2009 by the American Society of Plant Biologists