Skip to main content


Springer Nature is making SARS-CoV-2 and COVID-19 research free. View research | View latest news | Sign up for updates

Tobacco drought stress responses reveal new targets for Solanaceae crop improvement



The Solanaceae are an economically important family of plants that include tobacco (Nicotiana tabacum L.), tomato, and potato. Drought is a major cause of crop losses.


We have identified major changes in physiology, metabolites, mRNA levels, and promoter activities during the tobacco response to drought. We have classified these as potential components of core responses that may be common to many plant species or responses that may be family/species-specific features of the drought stress response in tobacco or the Solanaceae. In tobacco the largest increase in any metabolite was a striking 70-fold increase in 4-hydroxy-2-oxoglutaric acid (KHG) in roots that appears to be tobacco/Solanaceae specific. KHG is poorly characterized in plants but is broken down to pyruvate and glyoxylate after the E. coli SOS response to facilitate the resumption of respiration. A similar process in tobacco would represent a mechanism to restart respiration upon water availability after drought. At the mRNA level, transcription factor gene induction by drought also showed both core and species/family specific responses. Many Group IX Subgroup 3 AP2/ERF transcription factors in tobacco appear to play roles in nicotine biosynthesis as a response to herbivory, whereas their counterparts in legume species appear to play roles in drought responses. We observed apparent Solanaceae-specific drought induction of several Group IId WRKY genes. One of these, NtWRKY69, showed ABA-independent drought stress-inducible promoter activity that moved into the leaf through the vascular tissue and then eventually into the surrounding leaf cells.


We propose components of a core metabolic response to drought stress in plants and also show that some major responses to drought stress at the metabolome and transcriptome levels are family specific. We therefore propose that the observed family-specific changes in metabolism are regulated, at least in part, by family-specific changes in transcription factor activity. We also present a list of potential targets for the improvement of Solanaceae drought responses.


Among the flowering plants, the Solanaceae (nightshade family) ranks third (after grasses and legumes) as the most important crop for human beings [50]. This highly diverse group is comprised of 90 genera and 3000–4000 species, in which half of the species belongs to the genus Solanum [34]. Several members of this genus, such as tomato (S. lycopersicum), potato (S. tuberosum) and eggplant (S. melongena) are important to human diet. Other members of the family are utilized for drug production [34]. This angiosperm family is also interesting because a number of its members such as Nicotiana spp, Solanum spp., Petunia spp., and Datura spp are used as biological model systems [56].

Tobacco is an agriculturally important Solanaceae crop and is one of the most studied plants as biological model system [51]. It is a convenient plant for research because it can be easily transformed and has a relatively short generation time. A cell line (BY-2) is a popular system for functional genomics research because of its fast growth, responsiveness to a variety of plant hormones and ease of transformation. There is also significant interest in understanding genome evolution in tobacco and the Solanaceae [60].

Drought is one of the major constraints in crop production and affects 64% of the global land area [11,48,65]. It is the most common cause of severe food shortage in developing countries [10]. Using candidate gene approaches, there have been a number of reports of genes that improve drought tolerance in tobacco [35]. These include DREB and WRKY transcription factors, genes that alter the levels of trehalose and mannitol, and LEA genes [31,46,63,64,69,70]. The majority of these genes came from bacteria or other plant species and tobacco was used as a model system [31,46,63,64,69,70]. In order to improve drought responses in tobacco and other Solanaceae species, it is desirable to understand how these plants respond to drought stress at multiple levels in the plant. This allows a better understanding of primary and secondary metabolism and the interplay between transcriptional, posttranscriptional, translational and posttranslational regulation. However, there are few reports on metabolic changes in tobacco during drought stress and systems biology data that combine results from multiple different levels in the same samples appear to be lacking.

Here we present extensive data sets on the response of tobacco to drought stress at the physiological, mRNA, metabolite, and promoter levels. We identify novel aspects to the response, such as the accumulation of 4-hydroxy-2-oxoglutaric acid (KHG) in roots, and identify genes, metabolites, transcription factors, and promoters that are potential components of novel strategies to improve drought stress responses in Solanaceae crops. These data provide a framework for crop improvement and are timely because of the recent publication of the tobacco genome sequence [60]. They also provide extensive novel resources for comparative analyses.


The physiological level

Plants typically have to respond to multiple abiotic stresses simultaneously. This makes it difficult to characterize the signalling web that is associated with any one particular stress. To investigate drought responses we therefore performed experiments using hydroponic conditions where temperature, relative humidity, and the light regime were controlled. Tobacco cv ‘Burley 21’ plants were subjected to dehydration through air-drying for 20, 40, 60, 120 and 240 minutes by removing the plants from the hydroponics solution without touching the plants. This strategy ensured a strong and uniform response of the plants to drought stress, rapid harvesting, and a lack of wounding of the tissues. Plants showed symptoms of wilting after 20 minutes of dehydration. Further drought stress up to 240 minutes led to severe wilting. The tobacco plants were still alive after 240 minutes of dehydration because re-watered plants were able to fully recover showing expanded leaves. To determine the water status of the tobacco plants during the time course, osmotic potential and relative water content (RWC) were measured (Figure 1). The kinetics of RWC and osmotic potential changes differ with the greatest change in relative water content occurring in the first 20 minutes whereas the greatest change in osmotic potential occurs in the last 120 min. Stomatal conductance (SC) decreased rapidly (Figure 1) with a dramatic fall from 333 mmol/m/s (unstressed) to 86 mmol/m/s after 20 minutes indicating stomatal closure. Taken together, the physiological data show that stomatal closure is one of the most rapid responses of tobacco plants to drought stress and that the initial fall in osmotic potential occurs more rapidly in leaves than in roots.

Figure 1

Relative water content (RWC, %), osmotic potential (MPa) and stomatal conductance (mmol/m/s) of leaves and osmotic potential of tobacco roots under moderate (20–60 min) and extreme (120–240 min) drought stress. Error bars (standard error) were calculated from three replicates.

The metabolite level

A list of the notable changes in metabolite concentrations during drought is presented in Table 1 and they identify potential metabolic pathways that might be targeted for improving water stress responses. Levels of metabolites were determined by liquid chromatography/mass spectrometry (LC/MS, LC/MS2) and gas chromatography/mass spectrometry (GC/MS) platforms. The resultant dataset comprised 116 named biochemicals in leaves and 140 in root. Following log transformation and imputation with minimum observed values for each compound, Welch’s two-sample t-test was used to identify metabolites that differed significantly between experimental groups. Metabolites were identified that achieved statistical significance (p ≤ 0.05), as well as those approaching significance (0.05 < p < 0.1). An estimate of the false discovery rate (q-value) was also calculated (Additional file 1: Table S1). Our data show that drought stress caused considerable alterations in metabolite profiles in both leaf and root, even at the first time point (20 minutes) (Figure 2 and Additional file 1: Table S1) with changes in the biochemical profile of root tissues far more extensive than leaf (Figure 2). Nearly forty metabolites in the roots increased in level at least 10-fold during drought stress. By comparison, no metabolites in the leaves increased more than 5-fold. Principal Component Analysis (PCA) plots also reflect these trends (Additional file 2: Figure S1). In leaf the 60 min samples grouped more closely to the 0 min control than other groups, consistent with the biphasic nature of the response. In the root plot the time points in the first hour grouped relatively tightly, while the 120 and 240 min samples are well separated.

Table 1 Features of the response of tobacco to drought stress mRNA, metabolite, and promoter levels
Figure 2

Heat map of relative changes in metabolites in leaf and root tissue during drought stress. Samples were analyzed by liquid chromatography/mass spectrometry (LC/MS, LC/MS2) and gas chromatography/mass spectrometry (GC/MS) platforms. The dataset comprised a total of 116 named biochemicals in leaves and 140 in root. Numbers indicate the time after the start of the experiment and the color scales the extent of changes in metabolites. Colored bars on the left of the heat map show the categories of metabolites.

It is striking that the levels of many metabolites increased relative to the time zero control (Additional file 1: Table S1). The few compounds which trend down over time did not appear random. They were mostly in the amino acid category and several were directly related to nitrogen metabolism (e.g. asparagine, aspartate) or were associated with the photorespiration pathway (glycine, serine), which can produce ammonia. Also, the one compound in nucleic acid metabolism which trended down was also a nitrogen storage molecule, allantoin. These patterns in the nitrogen pathways suggest impending nitrogen toxicity.

Compatible solute accumulation

Our data showed evidence of drought and oxidative stress responses, as indicated by increased levels of compatible solutes such as sugar alcohols, amino acids, and oligosaccharides, as well as oxidative products of reactive oxygen species (ROS) remediation. The sucrosyl-oligosaccharide pathway, also known as the raffinose pathway, serves to transport and store carbon, and the resulting compounds can serve as osmolytes in situations of low water potential. In the leaf galactinol and raffinose were undetectable until the final 240 minutes time point, suggesting an activation of the pathway due to water stress (Figure 3). A decline in hexose-phosphates (glucose-6-P, fructose-6-P) and an increase in mannitol were also consistent with a net flow of carbon into the compatible solute forms (Figure 3). In roots there was a transient increase of galactinol and myo-inositol, but then a decline at later time points. However, the later time points were marked by a sharp increase in mannitol and trehalose, other compatible solutes often associated with drought stress. It may be that roots follow a different regulatory strategy, or the kinetics of carbon redistribution may be different in roots and leaves.

Figure 3

Changes in the sucrosyl-oligosaccharide/raffinose pathway during drought stress. The pathway is shown together with the relative levels of selected metabolites in the pathway. Error bars represent +/− one standard deviation. Green lines denote changes in leaves and brown lines changes in roots. X-axes numbers represent thedehydration time points in minutes. Y-axes values show the relative levels of the metabolite normalized to the median.

Stress related compounds

Increases in several other stress related compounds were apparent, including oxidized glutathione (GSSG), threonate and dehydroascorbate (both products of ascorbate oxidation) (Additional file 1: Table S1). The glutathione-ascorbate cycle detoxifies hydrogen peroxide which is a reactive oxygen species. The cycle involves the antioxidant metabolites GSSG and dehydroascorbate and our data show that the glutathione-ascorbate cycle is being activated in tobacco as a response to drought stress.

The GABA shunt

The GABA shunt is a well-studied stress response pathway, the functions of which include controlling cytoplasmic pH, maintaining C/N balance by converting glutamate in the cytosol to succinate in the TCA cycle, and aiding in oxidative stress protection by generating NADH and succinate. Here we observed a rapid but transient increase in GABA in leaves, followed by a dramatic reduction after 120 min. In roots GABA slightly increased over the first hour, but then rose sharply after 120 minutes (Figure 4). This contrast in kinetics between tissues was similar to the patterns observed for other oxidative stress markers, such as glutathione, threonate, and dehydroascorbate.

Figure 4

Changes in the GABA Shunt during drought stress in leaves. The pathway is shown together with the relative levels of selected metabolites in the pathway. Error bars represent +/− one standard deviation. X-axes numbers represent dehydration time points in minutes. Y-axes values show the relative levels of the metabolite normalized to the median.

Nitrogen metabolism and maintenance of C/N balance

Nitrogen metabolism and maintenance of C/N balance is a key aspect not only of normal growth and development but also stress management in plants. A strong demand for carbon diversion into oligosaccharides (e.g. raffinose pathway), disaccharides (e.g. trehalose), or sugar alcohols (e.g. mannitol, sorbitol) as osmolytes can lead to excess nitrogen in the system (as toxic ammonia). Photorespiration also generates ammonia, resulting from the conversion of glycine to serine, and of serine to beta-hydroxypyruvate. The normal response to excess ammonia is to capture it through the nitrogen assimilation pathways (glutamine synthase, glutamate dehydrogenase, glutamine oxoglutarate aminotransferase) and sequester and/or transport excess nitrogen as asparagine and other nitrogen rich compounds such as allantoin. Failure of this mechanism can lead to ammonia toxicity and cell damage [44]. In our study, an increase in glutamate and glutamine at the late time points in both leaves and roots, suggested increased assimilation. However, at the same time there was a dramatic drop in aspartate levels, the precursor of asparagine, and a decrease in asparagine and allantoin (Additional file 1: Table S1). This suggests that the conversion of oxaloacetate to aspartate by aspartate aminotransferase (with the amino group donated by glutamate) may be limiting the flow of excess nitrogen into asparagine, and may contribute to ammonia toxicity.

The photorespiration pathway

The photorespiration pathway can generate ammonia, and the two amino acids involved in the pathway (glycine and serine) showed a common pattern contrary to the general amino acid accumulation (Figure 5). They both decreased during drought stress in leaves and roots, most dramatically in leaves. This suggests that tobacco tissues down-regulate photorespiration during drought stress as a mechanism to reduce the accumulation of toxic ammonia. This phenomenon is well documented [4]. The behaviour of another metabolite, 4-hydroxy-2-oxoglutaric acid (KHG, also known as 2-hydroxy-4 oxopentanedioic acid, 4-Hydroxy-2-ketoglutarate, 2-Keto-4-hydroxyglutarate and other names), may be related to this shut down of photorespiration. The levels of KHG show the greatest increase, relative to controls, of any compound in the study (Figure 5 and Additional file 1: Table S1). KHG was undetectable in leaves but there was a rapid 20-fold increase in KHG levels in root tissues between 60 and 120 minutes of drought stress. After four hours, KHG levels reached a striking 70-fold higher than control root tissue. It is not clear what the role of KHG is and further plant-based work is required. KHG can be broken down to pyruvate and glyoxylate by the corresponding aldolase and in other systems the aldolase plays a role in respiratory metabolic pathways, for example during respiration resumption during the termination of the E.coli SOS response [5]. It is possible that KHG accumulates in tobacco during drought stress and is then broken down into pyruvate and glyoxylate when water is available again. This would represent a novel mechanism that is used by tobacco plants to restart respiration upon water availability after drought.

Figure 5

Changes in the Photorespiration pathway during drought stress. The pathway is shown together with the relative levels of selected metabolites in the pathway. Error bars represent +/− one standard deviation. Reactions are shown in the mitochondria, the peroxisomes, and the cytoplasm. Green indicates results for leaf and brown for roots. X-axes numbers represent the dehydration time points in minutes.Y-axes values show the relative levels of the metabolite normalized to the median.

The mRNA level

Gene expression at the mRNA level was monitored using a Nimblegen 12 x plex custom oligo array. The oligo array contained multiple 60mer probes to individual gene sequences from three different sources. Firstly, 40,000 individual genomic survey sequence reads (GSSs) from the Tobacco Genome Initiative project with the highest E-value hits to proteins in the database, secondly, all TOBFAC transcription factors [52], and thirdly the Version 3.0 DFCI Tobacco Gene Index EST sequences. This ensured a wide coverage of tobacco genes (Additional file 3: Table S2 and Additional file 4: Table S3). We previously produced a tobacco MapMan mapping that can be used directly with our NimbleGen oligoarray to visualize the oligo array data [39]. Using these data, we were able to identify genes that showed significant increases or decreases at the mRNA level during drought stress and correlate these data with changes at the metabolite level.

A MapMan visualization of genes involved in secondary metabolism during drought shows that there are major differences between leaf and root tissues at the mRNA level (Figure 6). In roots, genes involved in anthocyanin production were among the most strongly up-regulated genes. Also up-regulated are genes involved in the biosynthesis/catabolism of glucosinolates, chalcones, flavonols, and terpenoids (Figure 6). In leaves, genes involved in anthocyanin, chalcone, phenylpropanoid, lignin and wax metabolism are strongly up-regulated.

Figure 6

MapMan visualization of genes involved in secondary metabolism. Pathway visualization was performed using our tobacco MapMan mapping [39]. Blue denotes genes that are up-regulated and red genes that are down-regulated. The color scheme is shown in a rectangle using a log 2 scale. The most intense color therefore represents changes of 8-fold or more. A. Changes in root after 40 minutes. B. Changes in root after 4 hours. C. Changes in leaf after 40 minutes. D. Changes in leaf after 4 hours.

In leaf, after forty minutes of dehydration, approximately two hundred genes showed at least 8-fold induction in mRNA levels (Additional file 5: Table S4). This number increased to over eight hundred after four hours. The up-regulated genes at the early time points included a large percentage of signalling molecules such as transcription factors, protein kinases, F Box proteins, and phosphatases (Additional file 5: Table S4). Over the time course of drought stress, the most strongly up-regulated gene in tobacco leaves also encodes a signalling molecule, a U-Box E3 ubiquitin ligase that we have called Ubiquitin Protein Ligase Like 1 (NtUPLL1) (Additional file 5: Table S4). The level of NtUPLL1 mRNA increases 52-fold after 40 minutes and over 290-fold by the four hour time point. NtUPLL1 and the closely related NtUPLL2 are similar to the Arabidopsis U-Box E3 ubiquitin ligases AtPUB18 and AtPUB19 that regulate ABA-mediated stomatal closure and drought stress responses [58].

At the later time points in leaf, the list of up-regulated genes includes many biosynthetic enzyme genes that produce osmoprotectants such as raffinose synthase and galactinol synthase and this reflects the observed changes in metabolites (Figure 3). Other well documented drought stress-inducible genes such as dehydrins, aquaporins, and LEA proteins also show induction. This includes glutathione peroxidase genes (Additional file 5: Table S4) which are considered the main enzymatic defence against the oxidative destruction of membranes [3,67]. Glutathione peroxidases may play dual roles, the first in the control of H2O2 homeostasis, and the second in linking ABA and H2O2 signalling in stomatal closure, thereby regulating transpiration [43]. This observation is consistent with an increased flow through the glutathione-ascorbate cycle at the metabolite level. We observed a 33-fold increase in 5-Epiaristolochene 1,3-Dihydroxylase mRNA levels. 5-Epiaristolochene 1,3-Dihydroxylase catalyses the formation of capsidiol, a sesquiterpene produced by Solanaceous plants in response to stresses [68]. Several glutamate decarboxylase (GAD) genes were up-regulated at the mRNA level and this represents a correlation between increases in the metabolite and mRNA levels. GAD catalyses the decarboxylation of glutamate to GABA and CO2 and these data show that the GABA shunt stress response pathway is up-regulated by drought stress in tobacco at the mRNA level and this is subsequently reflected in an increase in GABA at the metabolite level.

In roots at the earliest time point, the up-regulated genes were again predominately signalling molecules including AP2/ERF and MYB transcription factors, MAP kinases and calcium dependent protein kinases (Additional file 5: Table S4). Calmodulin-related genes are also prominent in the list of up-regulated genes. Other strongly induced genes include two CYP94C1 and CYP94B3 genes that catalyse two successive oxidation steps in JA-Ile catabolic turnover. Both JA and ethylene biosynthetic enzyme genes were strongly induced, as were several protein phosphatase 2C genes. Similar to leaf tissues, at later time points, the list of up-regulated genes includes many biosynthetic enzymes genes that produce osmoprotectants. Other notable genes include an anthocyanidin synthase gene (61-fold induced after one hour) that catalyses the penultimate step in the biosynthesis of anthocyanins [6] and a UDP-glycosyltransferase 74B1-like gene (54-fold induced after four hours) involved in glucosinolate biosynthesis. These genes are potential targets for engineering secondary metabolism (anthocyanins and glucosinolates, respectively) and drought tolerance (Table 1).

The phosphoinositide (PI) pathway is an important regulator of cellular functions [12] and the PI pathway and inositol-1,4,5-trisphosphate are implicated in plant responses to stress [45]. An inositol polyphosphate 5-phosphatase-like gene (AM835516) is strongly up-regulated in roots. Arabidopsis plants overexpressing a mammalian inositol polyphosphate 5-phosphatase lost less water and exhibited increased drought tolerance. The onset of drought stress was delayed in these transgenic plants and ABA levels increased less than in wild-type plants [45]. The tobacco inositol-1,4,5-trisphosphate 5-phosphatase enzyme is therefore predicted to modulate the phosphoinositide pathway, ABA levels and drought responses and could be an excellent tool to improve drought stress responses in the Solanaceae (Table 1).

Plant hormones

Both the metabolomics and transcriptomics data show roles for the plant hormones ABA, JA, and ethylene in the response to drought stress in tobacco (Additional file 6: Figure S2). At the metabolite level, ABA (abscisate) increased 8-fold after four hours in root tissue. This observed kinetics of ABA increase suggests that many immediate early genes are activated before any increase in ABA levels (Additional file 5: Table S4 and Additional file 7: Table S5) and this is consistent with observations that early responses may by triggered by a rapid hydraulic signal in an ABA-independent manner [8,9]. There was also a biphasic increase in N-delta-acetylornithine, which has been shown to rise in response to methyl-jasmonate. At the transcriptome level, MapMan visualization revealed that ABA, JA, and ethylene all play major roles in regulating drought stress (Additional file 6: Figure S2).

The biosynthetic enzyme genes in the pathway that produces JA were all up-regulated in tobacco root tissue and to a lesser extent in leaves (Additional file 6: Figure S2). The biosynthetic enzyme genes in the ethylene biosynthetic pathway are also up-regulated (data not shown). We found numerous other differentially regulated genes that provide insights into ABA and JA signalling in drought stressed tobacco plants. A ninja-family protein AFP3/ABI five-binding protein 3 gene is strongly induced in leaves at later time points. NINJA proteins acts as negative regulators of ABA responses and stress responses [19] and are part of a repressor complex that negatively regulates JA signalling [1]. Interestingly, two closely related cytochrome P450 CYP94 genes (CYP94C1 and CYP94B3) are strongly and transiently induced in root tissue (Table 1). In Arabidopsis, these two enzymes catalyse two successive oxidation steps in JA-Ile catabolic turnover [23] and these data suggest that drought stress in tobacco activates the genes for a major catabolic route for turning over JA-Ile. Our data show that not only genes for JA catabolism appear to be up-regulated but also genes for ABA catabolism. ABA 8'-hydroxylase (CYP707A1) genes are strongly and transiently up-regulated in the leaf (Table 1) and encode key enzymes in controlling the level of ABA in plants [36] with the hormonal action of ABA being controlled by the precise balance between its biosynthesis and catabolism.

Transcription factors

Transcriptional reprogramming is a central component of the response to drought stress. For this reason we investigated the transcription factors (TFs) that are up-regulated by drought stress. A MapMan visualization of the TF genes present on the oligo array is presented in Additional file 8: Figure S3. It highlights the tissue-specific differences in the transcriptional regulation of TFs during drought stress in tobacco.

We produced a non-redundant list of the thirty most highly induced TF genes in leaf and root tissues at an early (40 minutes) and late (four hours) time point (Additional file 7: Table S5). All TF genes are named according to the TOBFAC database nomenclature [52]. After 40 minutes of drought stress in root tissues, AP2/ERF TFs were the main group of up-regulated genes (18 genes), followed by MYB, bHLH, WRKY, C2H2 zinc finger, JAZ/TIFY, and bZIP (Additional file 7: Table S5). The most highly induced TF gene was NtERF187 (295-fold) which encodes a protein similar to Arabidopsis AtERF53 (Table 1) which plays a role in drought-regulated gene expression [7,24]. The second most highly induced gene was NtMYB149 which is similar to Arabidopsis AtMYB15. AtMYB15 is involved in ABA-, ethylene-, and JA-mediated signalling pathways, the response to salt stress, and the response to water deprivation. Overexpression of AtMYB15 confers enhanced sensitivity to ABA and improved drought tolerance in Arabidopsis [14] and this suggests that NtMYB149 is a good candidate for improving drought tolerance in the Solanaceae. The next highest induced TF gene was NtERF13 which is in the same clade as CBF1, CBF2, and CBF3 that play roles in cold and dehydration stress. Clearly some of the major TF nodes in drought stress signalling are conserved between tobacco and Arabidopsis and probably across dicot plants as a whole. After four hours, the same three TF genes (NtERF187, NtERF13, and NtMYB149) were still the most highly induced in roots. However, in addition three NAC TFs and three WRKY TFs are also among the most highly induced genes at this time point (Additional file 7: Table S5).

In leaf, there were major differences in the TF mRNA profile compared to root (Additional file 7: Table S5). Although more AP2/ERF genes were up-regulated than other TF families (nineteen genes), WRKY genes (seven) constituted a considerable proportion of the list of thirty most strongly up-regulated genes. The two most strongly up-regulated genes after 40 minutes were the closely related NtERF202 and NtERF114 (also called ACRE111A and ACRE111B) (Table 1). Both are similar to Arabidopsis CBF3/DREB1A. NtWRKY1 is among the most highly induced genes during the early response (Table 1). Prominent among the up-regulated WRKY TF genes, NtWRKY1 is the apparent orthologue of AtWRKY33, which has been shown to play major roles in the response to stresses including abiotic stress [28]. The late time point in leaves revealed a much more diverse profile of differentially regulated TF genes with representatives from the AP2/ERF, WRKY, JAZ/TIFY, NAC, bZIP, MYB, GRAS, homeodomain, and AUX-IAA families.

One of our aims was to identify genes that might be used to improve drought stress responses, not just in tobacco but also in other members of the Solanaceae. We therefore compared drought stress inducible WRKY genes from soybean (Additional file 9: Table S6), Arabidopsis, tobacco, and tomato to determine commonalities that may represent core changes in transcription factor gene expression and differences that could be responsible for family/species-specific responses to drought stress. The drought stress analyses were obtained as follows: tobacco (this study, GEO accession GSE67434), soybean (GEO accession GSE49537), Arabidopsis [33,72] and tomato [20,25,38]. There are three major hotspots of drought stress inducibility in the WRKY phylogenetic tree. The first is in Group III with two sets of tobacco/tomato orthologues (NtWRKY15/NtWRKY95/21/SlWRKY41 and NtWRKY103/SlWRKY54) and inducible Arabidopsis and soybean genes. The second is in Group I in a clade that includes NtWRKY1 and SlWRKY31. In the data sets studied, NtWRKY1 is the most strongly induced WRKY gene after 20 minutes of drought stress in tobacco leaves and SlWRKY31 is the most strongly induced WRKY gene in drought stressed tomato leaves. This suggests that induction of NtWRKY1/SlWRKY31-like genes is a feature of drought responses in the Solanaceae and that these genes are good candidates for manipulating drought responses. The final hotspot is in Group IIa, where NtWRKY86 and SlWRKY39 are induced by drought stress as well as several Arabidopsis and soybean genes (Figure 7). Group IIa is a well-studied group and includes AtWRKY18, AtWRKY40, and AtWRKY60 that play central roles in the responses to multiple stresses [54,65].

Figure 7

Drought-inducible transcription factor genes in the WRKY, HDZIP, and AP2/ERF families. A. Core responses in the WRKY transcription factor family B. Family/species-specific responses in the AP2/ERF and WRKY families C. Core responses in the AP2/ERF family D. Core responses in the HDZIP family. Drought stress inducibility or association with the GO term “Response to water deprivation” (GO:0009414) is denoted by red circles (soybean), blue circles (Arabidopsis), green circles (tobacco), and pink circles (tomato). Marker domains for Group III and X ERF subfamilies are shown as purple and light blue squares respectively. The evolutionary histories were inferred using the Neighbor-Joining method. The optimal trees are shown. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) is shown next to the branches. The evolutionary distances were computed using the Poisson correction method and are in the units of the number of amino acid substitutions per site. For the HDZIP transcription factors, induction by drought stress and the processes that the transcription factors are involved in are shown. AP2/ERF and homeodomain genes are shown by gene model name except for tobacco where the TOBFAC nomenclature is used. WRKY genes from Arabidopsis are shown by their commonly used names, soybean from our own analysis (Additional file 9: Table S6), tobacco genes were taken from TOBFAC and tomato genes from published sources.

Interestingly, as well as similarities in TF gene induction our data suggest that there may be Solanaceae-specific drought stress inducibility of several Group IId WRKY genes from tobacco and tomato (Figure 7). Tobacco NtWRKY3, NtWRKY10, NtWRKY12, and NtWRKY69, together with SlWRKY10 all show drought stress inducibility whereas under similar experimental conditions, no soybean IId genes showed this induced expression. In addition, the Arabidopsis IId genes AtWRKY7, AtWRKY11, and AtWRKY17 have been shown to modulate transcription in response to pathogen challenge [29] but have also not been shown to play any role in responses to drought.

There are notable similarities in the ERF/AP2 family when we compared drought stress-inducible genes from tobacco and soybean with Arabidopsis genes that are annotated with the GO term “Response to water deprivation” (GO:0009414). As expected, a large number of genes from the DREB subfamily (Groups I, II, IV, and particularly III) are induced by drought stress in tobacco (Additional file 8: Figure S3). Group III (including the Arabidopsis genes AT4G25490/CBF1/DREB1B and AT4G25480/CBF3/DREB1A), is a hotspot with at least 20 out of 36 genes being associated with the response to drought stress (Figure 7). In the ERF subfamily, a small number of genes from Group VIII, a single Group VI gene and no Group VII genes showed inducibility by drought stress. By contrast, many Group IX and X genes showed induction. The functions of many ERF subfamily transcription factors in abiotic stress responses are largely unknown [42] and it was therefore of interest that Group X is a hotspot of drought stress responsive ERF/AP2 genes (Figure 7). Our data suggest that Group X ERF transcription factors play roles in the response to drought stress but unlike the well-studied DREB transcription factors, the exact roles of many remain to be determined.

We again observed apparent family-specific drought stress-inducibility of ERF/AP2 genes. In Subgroup 3 of the Group IX ERF genes there are no fewer than ten soybean ERF genes that are inducible by drought stress. Strikingly, not a single tobacco gene showed induction (Figure 7). This appears to be a family/species-specific role for soybean ERF transcription factors in the response to drought stress. This is particularly noticeable in the area of the phylogenetic tree that includes GLYMA10g00980, GLYMA02g00870, GLYMA10g33060, and GLYMA20g34570. GLYMA20g34570 is induced 35-fold by drought stress in leaves and GLYMA10g00980 is induced 128-fold making it one of the most highly up-regulated genes by drought stress in soybean leaves. By contrast, the most similar tobacco gene NtERF22 (Additional file 10: Figure S4) shows no induction by drought stress whatsoever (Additional file 11: Table S7). This poses the question what the tobacco Subgroup 3 transcription factors may be doing and our previous work sheds light on this [57]. NtERF5, NtERF165, NtERF22 and NtS25XP1 have all been implicated in regulating jasmonate-inducible nicotine biosynthesis in tobacco [57]. Interestingly, nicotine biosynthesis occurs only in the Solanaceae, predominantly in tobacco, and in lower quantities in tomato, potato, eggplant, and green pepper. This suggests that family-specific differences in secondary metabolism are reflected in family-specific roles of some transcription factors.

There are also similarities in drought stress responsiveness in the homeodomain family. In tobacco, two HDZIP genes (HD1 and HD20) are strongly up-regulated by drought stress in leaf and Figure 7 shows that the two orthologues in soybean (Glyma16g02390 and Glyma07g05800) are also very strongly induced by drought stress. Strikingly, the two Arabidopsis orthologues of these genes (AT2G46680/ HOMEOBOX 7 and AT3G61890/HOMEOBOX 12) have been shown to play roles in the responses to drought stress and our data suggest that orthologues of these HDZIP transcription factors regulate drought stress responses across dicot plants and possibly beyond.


To validate and characterize drought stress inducibility, several genes were selected for promoter:reporter analyses in transgenic tobacco plants. We chose five drought stress inducible genes, galactinol synthase (NtGolS), NtUPLL2, NtWRKY69, NtWRKY70, and NtWRKY3. An additional two promoters (NtWRKY95/21 and a raffinose synthase) have previously been validated in transgenic tobacco BY-2 cells [47] where they are inducible both by polyethylene glycol and jasmonate. Our choice of promoters was partly dictated by available promoter sequence. For example, promoter sequences of NtWRKY86, and NtUPLL1 were not available.

Transgenic tobacco plants containing promoter-GUS constructs were tested for reporter gene activity under dehydration, cold and wounding. Figure 8 shows GUS expression from the different promoters during a time course of drought stress. The NtWRKY69 and NtWRKY3 promoters directed early induction by drought stress in leaves with GUS activity observed as early as two hours after the start of dehydration (Figure 8). The kinetics of GUS activity directed by the promoters in planta was similar to the kinetics of mRNA accumulation observed using the oligo array. Both NtWRKY3 and NtWRKY69 show early increases in mRNA levels and these increases are reflected in early promoter activities in planta. In contrast, promoters from NtWRKY70, NtGolS and NtUPLL2 directed later induction by drought stress. Again this correlates with mRNA levels that show major increases during the later stages of dehydration.

Figure 8

Activity of promoter-GUS reporter gene constructs in 4-week old transgenic tobacco plants at different drought stress time points. Inset line graphs show the expression profiles of the genes at varying dehydration time points at the mRNA level from oligo array analysis.

To determine whether these promoters also respond to other types of abiotic stress, the transgenic plants were subjected to cold and wounding (Figure 8). All promoter:reporter transgenic lines showed inducibility by cold treatment (Figure 8). In addition the promoters also showed varying degrees of inducibility by wounding (Figure 8). NtWRKY3, NtWRKY69, NtWRKY70, NtGolS and NtUPLL2 are all therefore responsive to multiple abiotic stresses and are components of interconnected abiotic stress-responsive signalling webs.

We chose to further study the temporal and spatial activation of the NtWRKY69 promoter in more detail because it is one of the Group IId genes that our data suggest may be part of a Solanaceae-specific response to drought stress (Figure 7). Using GFP as the reporter gene enabled monitoring of promoter activity over the course of drought stress and a remarkable activation of the NtWRKY69 promoter was observed (Figure 9). In the absence of drought stress, little or no activity was seen in the leaves. After the start of drought stress, promoter activity started at the base of the stem and then moved up the stem and into the leaf. Promoter activity was subsequently found in the vascular tissue of the leaf and then eventually spread out into the non-vascular tissue. This suggests that a signal from the roots is moving through the stem and then into the leaves through the vascular tissue. Whether this signal is hormonal (such as ABA and/or JA) or hydraulic is unclear. To determine whether the signal might be ABA, 20 μM ABA was sprayed onto the leaves of the transgenic tobacco. The treated leaf areas showed no GFP activity. As a control, the ABA-treated plants were then subjected to dehydration and inducible promoter activity was again seen (data not shown). This suggests that the promoter of NtWRKY69 is not responsive to exogenous ABA and therefore that NtWRKY69 functions in an ABA-independent pathway. Taken together, our data suggest roles for JA-dependent, ethylene-dependent, ABA-dependent and ABA-independent pathways in the water stress induced signalling web in tobacco.

Figure 9

NtWRKY69 promoter activity in 4-week old transgenic tobacco plants during drought stress. ProNtWRKY69:GFP-containing plants illustrate the progression of promoter activation in stems and then leaves during drought stress.


Core and family/species-specific features of the drought stress response in tobacco

The responses of many plant species to drought stress have been extensively studied [13,30,33,59,61,66,71]. Several hormones play major roles in drought responses, with ABA being a major player. ABA plays a major role in drought signalling and water stress causes an immediate hydraulic signal that triggers ABA biosynthesis [49]. ABA perception occurs by a mechanism where binding of ABA to the ABA receptors RCARs/PYR1/PYLs results in inactivation of type 2C protein phosphatases including ABI1 and ABI2. These type 2C protein phosphatases appear to function as co-receptors and they inactivate SNF1-type kinases. This results in activation of ABA-dependent gene expression [49]. To understand the mechanisms involved in the response to drought stress, it is necessary to understand the changes in the plant cell at many interconnected levels. We therefore sought to correlate changes at different levels by using the same plant samples for physiological, mRNA, and metabolite analyses. We have also been able to suggest common features of plant drought stress responses and additonally features that may be unique to tobacco or members of the Solanaceae. Comparative analyses of large omics data sets are hampered by the varying experimental setup [32]. However, the data set from soybean (GEO accession GSE49537) was produced under very similar conditions of hydroponics/dehydration to the tobacco data set and facilitates direct comparisons between a member of the Solanaceae and a legume. Our comparisons, together with recent reports from other groups [13,30,59,61,66,71] lead us to suggest that the following represent some of the core metabolic changes of plants to drought stress.

  1. 1)

    The production of glutathione and tocopherol as anti-oxidants.

  2. 2)

    Ammonia detoxification.

  3. 3)

    Increases in amino acids that act as osmolytes.

  4. 4)

    Activation of the raffinose pathway to produce osmolytes.

  5. 5)

    Regulation of nitrogen metabolism (asparagine, allantoin and glutamine).

  6. 6)

    Activation of the GABA Shunt to control cytoplasmic pH, maintain C/N balance and protect from oxidative stress.

The suggestion concerning the detoxification of ammonia is based mainly on the present study and the study by Oliver et al. [44] on desiccation tolerance in Sporobolus stapfianus. These studies point to ammonia detoxification as a core metabolic process during water stress. The ammonia itself would appear to come from a strong demand for carbon diversion into oligosaccharides (e.g. raffinose pathway), disaccharides (e.g. trehalose), or sugar alcohols (e.g. mannitol, sorbitol) and these osmolytes can lead to excess nitrogen in the system (as toxic ammonia). Photorespiration also generates ammonia, resulting from the conversion of glycine to serine, and of serine to beta-hydroxypyruvate. The exact role of ammonia detoxification in water stress responses and the mechanisms involved require further study.

In addition to these core responses, there are other responses that appear species/family specific. For example KHG, the most strongly induced metabolite in tobacco is not induced by drought stress in all plants. Our data show that although KHG is detected in drought stressed soybean roots, it is not induced by drought stress (Table 2). This shows that major responses to drought stress are family specific and not part of the core responses to drought stress.

Table 2 Relative increases in 4-hydroxy-2-oxoglutaric acid levels in tobacco and soybean during dehydration

Similar to the metabolite level, at the mRNA level transcription factor genes show both core changes and family specific responses in transcription factor gene expression. Core changes include up-regulation of members of the Group I, IIa, and III WRKY, Group III and X AP2/ERF, and HDZIP (Figure 7) transcription factors. Family-specific changes include Group IX AP2/ERF genes from soybean that are up-regulated by drought stress and Group IId WRKY genes from tobacco and tomato that likewise show drought stress inducibility. Taken together, it is likely that the family-specific changes in metabolism are regulated, at least in part, by family-specific changes in transcription factor activity. A challenge for the future is to further establish the roles of individual transcription factors in regulating these family-specific changes in metabolism. However, our data has already provided insights into this. Many Group IX Subgroup 3 ERF transcription factors in soybean appear to play roles in the response to drought stress (Figure 7), whereas their counterparts in tobacco appear not to regulate drought stress responses but rather the Solanaceae-specific increase in nicotine biosynthesis as a response to herbivory [57]. These observations provide a framework for the parallel study of family-specific metabolic changes and family-specific differences in transcription factor function. The large nature of our data set has precluded a study of all transcription factor families, but we are confident that similar core and family-specific changes are to be found outside of the AP2/ERF, WRKY, and homeodomain families. Indeed, an elegant recent study has suggested that there are core environmental stress response genes that are coordinately regulated not only as a result of drought stress but also more widely as a response to abiotic stresses [22].

The very high accumulation of KHG in tobacco roots strongly suggests that it is playing a role in combating drought stress but it is not currently clear what this mechanism is because KHG is poorly studied in plants. KHG can be broken down to pyruvate and glyoxylate by the corresponding aldolase (Figure 5) and in other systems the aldolase plays a role in respiratory metabolic pathways. One hypothesis comes from studies of the E. coli SOS response where KHG accumulates to promote respiration resumption during the termination of the E. coli SOS response [5]. We suggest a similar mechanism where KHG accumulates in tobacco during drought stress and is then broken down into pyruvate and glyoxylate when water is available again. This would represent a novel mechanism that is used by tobacco plants to restart respiration upon water availability after drought. Alternatively, the rise in activity of glutamine oxoglutarate aminotransferase (GOGAT) and it is possible that KHG serves as a scavenger for ferredoxin, NADP or the detoxification of NO2, which could have accumulated due to the dehydration. This and/or other functions of KHG in plants require further study.

Targets and new strategies for Solanaceae crop improvement

Table 1 lists some of the noteworthy features of the tobacco response to drought stress at the gene, metabolite, plant hormone, transcription factor, and promoter levels. Many of these represent potential targets for the improvement of Solanaceae crop plants such as tomato and potato. These systems biology data will prove valuable when designing strategies for crop improvement. For example, a drought stress-inducible promoter can be used to drive a transgene that affects the accumulation of identified key metabolites. An inducible promoter will restrict transgene expression to when and where it is needed and may reduce unwanted side effects caused by constitutive activation of stress responses. One notable example of this type of strategy comes from a member of the Solanaceae with transgenic potato tubers that expressed the Arabidopsis AP2/ERF gene AtDREB1A under the control of the drought stress-inducible RD29A promoter or the strong constitutive CaMV35S promoter [26]. Metabolite profiling of the transgenic lines revealed elevated levels of glutathione and GABA in addition to β-cyanoalanine, which is a biosynthesis by product of ethylene. Both glutathione and GABA are components of what we have described as the plant core metabolite response to drought stress and it appears that use of transcription factors similar to AtDREB1A (such as NtERF218 and NtERF228) as a transgene holds promise, especially in conjunction with inducible promoters such as the NtWRKY69 promoter (Figure 9).

Systems of reduced complexity have served plant scientists well in the discovery of many important components of the responses of plants to different abiotic stresses. Our data have likewise revealed new potential targets for crop improvement (Table 1). However, we are well aware that field conditions are much more complex that these hydroponics/growth room based studies and that our discoveries will need to be extended to field conditions. For example, studies have revealed that the response of plants to a combination of two different abiotic stresses is unique and cannot be directly extrapolated from the response of plants to each of the different stresses applied individually [41]. Drought tolerance is a complex quantitative and multigenic trait with a significant environmental component [37,40]. As a result, the genetic control of traits associated with tolerance to drought often shows low heritability. Therefore, one of the major hurdles in using the toolbox of genes, promoters and metabolites described in Table 1 to improve drought responses in Solanaceae crops lies not in producing these plants but in accurately defining them. The ongoing development of new phenomics technologies promises to greatly facilitate the characterization of transgenic lines, especially under field conditions. For this reason, field phenotyping systems may give us our best opportunity to determine whether strategies to improve drought stress responses can actually produce plants with increased drought tolerance.


We propose components of a core metabolic response to drought stress in plants but also show that major responses to drought stress at the metabolome and transcriptome levels are family specific.


Plant materials, growth conditions and drought stress treatments

Tobacco cv. ‘Burley 21’ seedlings were germinated on agar plates. When the plants had grown to almost reach the top of the petri dish, they were then transferred to a mini hydroponics set up in sterile MK-5 polycarbonate vessels with half-strength MS liquid medium. Plants were supported by plastic inserts with holes in them for the plants and grown for 2 weeks at 25°C. They were then subjected to dehydration in a growth room for 0 (control), 20, 40, 60, 120 and 240 minutes by removing them from the liquid using the supports. Each time point consisted of three replicates with 20 plants per replicate. Roots and leaves were harvested by flash freezing in liquid nitrogen. All samples from the time course of water stress were taken during the morning.

Soybean W-82 seeds were grown in hydroponics using 0.5 × Hoagland’s solution, pH 5.8 in a growth chamber with a 16 hour/8 hour day/night cycle at 25°C and 50% relative humidity. After 30 days, plants were subjected to drought stress for 30 minutes, one hour, two hours, three hours and five hours by removing them by means of the plastic lid supports without touching the plants. Leaves and roots were harvested by flash freezing in liquid nitrogen. Nine plants were utilized for each time-point (three replicates per time-point and three plants per replicate).

Osmotic potential, Stomatal conductance, and relative water content

A Vapro 5520 osmometer (Wescor, Logan, UT, USA) was used to measure the osmotic potential using 10 μL cell sap samples. Stomatal conductance was determined using an SC-1 leaf porometer (Decagon, Pullman, WA, USA). Relative water content was measured following the protocol of Gonzalez and Gonzalez-Vilar [21].

Total RNA extraction

RNA isolation was performed using an RNeasy kit (Qiagen, Valencia, CA, USA) following the manufacturer’s protocol. Genomic DNA contamination was removed using an Ambion DNA-free and DNase removal kit (Life Technologies, Carlsbad, CA, USA). 10 μg total RNA from each sample was used for micro-array analysis.


A Nimblegen custom oligo array was used for transcriptome analyses. The oligo array contained 385,000 probes with approximately seven 60-mer probes for each gene sequence. The sequences came from three different sources. Firstly, 40,000 individual genomic survey sequence reads (GSSs) from the Tobacco Genome Initiative that had the highest E-value hits to proteins in the NCBI nr database (A data set of 1,159,022 genomic survey sequences was downloaded from the TGI in 2008), secondly all TOBFAC transcription factors [52], and thirdly the Version 4.0 DFCI Tobacco Gene Index EST sequences (NTGI4;; file name: Source files/NTGI.071508.fasta). The genome survey sequences have subsequently been deposited at The National Center for Biotechnology Information; file name: Source files/TGSS_expressed.fasta although these deposited sequences were the result of re-reading to extend the length of the sequences and have more errors. Oligoarray experiments were performed at MOgene LC (St Louis, MO) using their standard protocols. Data analysis was performed using ArrayStar v4. Differential expression was calculated using 90% confidence (FDR Benjamini Hochberg) and 8-fold change as the cut off. We chose a high fold inducibility because under the experimental conditions very high inducibilities were obtained and a cut off value of 4-fold gave several thousand genes at the later time points. 8-fold was chosen to focus on only the most highly induced or repressed genes. 90% confidence was chosen because at early time points some genes failed the 95% confidence limit but as the time course progressed, these and other genes typically passed the 95% or even the 99% confidence limit. Pathway visualization was performed using our tobacco MapMan mapping [39].

For soybean, a custom made 12 × plex array was designed by Roche NimbleGen, Inc. containing multiple 60mer oligomers to all high and low confidence genes from the GLYMAv1.0 release of the soybean genome. This was used for soybean expression analyses. Hybridizations and data analyses were performed in an identical way to tobacco.


The same set of tobacco samples that was used for transcriptomic analyses was also used for metabolomics analysis at Metabolon, Inc. (North Carolina) using their analysis pipeline. The global unbiased metabolic profiling platform was based on a combination of three independent platforms: UHLC/MS/MS2 optimized for basic species, UHLC/MS/MS2 optimized for acidic species, and GC/MS. This platform has been described in detail [16]. Compounds were identified by comparison to library entries of purified standards or recurrent unknown entities, a total of over one thousand compounds. Following log transformation and imputation with minimum observed values for each compound, Welch’s two-sample t-test was used to identify biochemicals that differed significantly between experimental groups. Biochemicals that achieved statistical significance (p ≤ 0.05), as well as those approaching significance (0.05 < p < 0.1), were highlighted in the data set. An estimate of the false discovery rate (q-value) was also calculated to take into account the multiple comparisons that normally occur in metabolomic-based studies. The q-value describes the false discovery rate; a low q-value (q < 0.10) is an indication of high confidence in a result.

Phylogenetic analyses

The amino acid sequences of the DNA binding domains of WRKY, and AP2/ERF transcription factors, and the complete amino acid sequences of HDZIP transcription factors, were used to construct phylogenetic trees. Alignments were constructed using MUSCLE [15] and the following parameters; Gap Penalties: Gap open −2.9, Gap Extended 0, Hydrophobicity multiplier 1.2 Memory/Iterations: Max Memory in MB 4095, Max Iterations 8; Clustering Method Iteration 1, 2 (UPGMB), Clustering Method (Other Iterations (UPGMB), Min. Diag. Length (Lambda) 24. Alignments are presented as Additional file 12: Table S8, Additional file 13: Table S9, Additional file 14: Table S10. In each case, the evolutionary history was inferred using the Neighbor-Joining method [55]. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) were determined [17]. The evolutionary distances were computed using the Poisson correction method [73] and are in the units of the number of amino acid substitutions per site. All ambiguous positions were removed for each sequence pair. Evolutionary analyses were conducted in MEGA6 [62]. All positions containing alignment gaps and missing data were eliminated in pairwise sequence. The amino acid sequences of the DNA binding domains were taken from the TOBFAC database (tobacco), TAIR (Arabidopsis), Phytozome (soybean), and (tomato). Expression data came from tobacco (GEO accession GSE67434), soybean (GEO accession GSE49537), Arabidopsis [33] and tomato [20,25,38]. Data on Arabidopsis gene function were taken from the TAIR website (

Promoter analyses

The choice of promoters was dictated by available promoter sequence. The blast search and contig builder functions in the TOBFAC data base were used to extend genomic sequences to include the promoters. Where promoter sequences were obtained, the promoter regions were PCR amplified, verified by sequencing and inserted as HindIII-SacI fragments the pGPTV-GUS-KAN [2,53] or pGPTV-GFP-KAN binary vectors. Tobacco transformation was performed following the protocol of Gallois and Marinho [18]. Four-week old seedlings were subjected to different treatments. Dehydration was imposed by removing the plants from liquid medium into an empty MK-5 polyethylene vessel. The vessels with the plants were then placed into a large plastic container with water to allow for slow dehydration of the plant. Wounding was achieved by pressure exerted by serrated tweezers. Cold treatments were performed by transferring the plants to 4°C. GFP activity was monitored at 5, 12 and 24 h while GUS activity was observed after 24 h. Histochemical staining of GUS activity was performed as described by Jefferson et al. [27]. GFP visualization of transgenic plants was performed using an Olympus SZX16 Epi-Fluorescent stereo microscope. The promoter sequences of NtWRKY69, NtWRKY70, NtWRKY3, UPLL2, and galactinol synthase are available under the GenBank accession numbers x-x.

Protein motif analysis

Analysis of selected full length AP2/ERF transcription factors for protein domain architecture and putative conserved domains was carried out using MEME ( The settings were; any number of repetitions of a single motif, minimum width of a motif ten amino acids, maximum width of a motif sixty amino acids, maximum number of motifs to find nine.

Availability of supporting data

The tobacco oligo array data is available at the Gene Expression Omnibus online repository (NCBI GEO) as GEO accession GSE67434.

The soybean oligo array data is available at the Gene Expression Omnibus online repository (NCBI GEO) as GEO accession GSE49537 (


  1. 1.

    Acosta IF, Gasperini D, Chetelat A, Stolz S, Santuari L, Farmer EE. Role of NINJA in root jasmonate signaling. Proc Natl Acad Sci U S A. 2013;110:15473–8.

  2. 2.

    Becker D, Kemper E, Schell J, Masterson R. New plant binary vectors with selectable markers located proximal to the left T-DNA border. Plant Mol Biol. 1992;20:1195–7.

  3. 3.

    Brigelius-Flohe R, Aumann KD, Blocker H, Gross G, Kiess M, Kloppel KD, et al. Phospholipid-hydroperoxide glutathione peroxidase. Genomic DNA, cDNA, and deduced amino acid sequence. J Biol Chem. 1994;269:7342–8.

  4. 4.

    Brodribb T. Dynamics of Changing Intercellular CO2 Concentration (ci) during Drought and Determination of Minimum Functional ci. Plant Physiol. 1996;111:179–85.

  5. 5.

    Cayrol C, Petit C, Raynaud B, Capdevielle J, Guillemot JC, Defais M. Recovery of respiration following the SOS response of Escherichia coli requires RecA-mediated induction of 2-keto-4-hydroxyglutarate aldolase. Proc Natl Acad Sci U S A. 1995;92:11806–9.

  6. 6.

    Cheng L, Xu Y, Grotewold E, Jin Z, Wu F, Fu C, et al. Characterization of Anthocyanidin Synthase (ANS) Gene and anthocyanidin in rare medicinal plant-Saussurea medusa. Plant Cell Tissue Org Cult. 2007;89:63–73.

  7. 7.

    Cheng M-C, Hsieh E-J, Chen J-H, Chen H-Y, Lin T-P. Arabidopsis RGLG2, functioning as a RING E3 ligase, interacts with AtERF53 and negatively regulates the plant drought stress response. Plant Physiol. 2012;158:363–75.

  8. 8.

    Christmann A, Grill E, Huang J. Hydraulic signals in long-distance signaling. Curr Opin Plant Biol. 2013;16:293–300.

  9. 9.

    Christmann A, Weiler EW, Steudle E, Grill E. A hydraulic signal in root-to-shoot signalling of water shortage. Plant J. 2007;52:167–74.

  10. 10.

    Cominelli E, Tonelli C. Transgenic crops coping with water scarcity. N Biotechnol. 2010;27:473–7.

  11. 11.

    Cramer GR, Urano K, Delrot S, Pezzotti M, Shinozaki K. Effects of abiotic stress on plants: a systems biology perspective. BMC Plant Biol. 2011;11:163.

  12. 12.

    Dieck CB, Boss WF, Perera IY. A role for phosphoinositides in regulating plant nuclear functions. Front Plant Sci. 2012;3:50.

  13. 13.

    Dinakar C, Bartels D. Desiccation tolerance in resurrection plants: new insights from transcriptome, proteome and metabolome analysis. Front Plant Sci. 2013;4:482.

  14. 14.

    Ding Z, Li S, An X, Liu X, Qin H, Wang D. Transgenic expression of MYB15 confers enhanced sensitivity to abscisic acid and improved drought tolerance in Arabidopsis thaliana. J Genet Genomics. 2009;36:17–29.

  15. 15.

    Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.

  16. 16.

    Evans AM, DeHaven CD, Barrett T, Mitchell M, Milgram E. Integrated, nontargeted ultrahigh performance liquid chromatography/electrospray ionization tandem mass spectrometry platform for the identification and relative quantification of the small-molecule complement of biological systems. Anal Chem. 2009;81:6656–67.

  17. 17.

    Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39:783–91.

  18. 18.

    Gallois P, Marinho P. Leaf disk transformation using Agrobacterium tumefaciens-expression of heterologous genes in tobacco. In: Jones H, editor. Plant Gene Transfer and Expression Protocols. New York: Springer; 1996. p. 39–48.

  19. 19.

    Garcia ME, Lynch T, Peeters J, Snowden C, Finkelstein R. A small plant-specific protein family of ABI five binding proteins (AFPs) regulates stress response in germinating Arabidopsis seeds and seedlings. Plant Mol Biol. 2008;67:643–58.

  20. 20.

    Gong P, Zhang J, Li H, Yang C, Zhang C, Zhang X, et al. Transcriptional profiles of drought-responsive genes in modulating transcription signal transduction, and biochemical pathways in tomato. J Exp Bot. 2010;61:3563–75.

  21. 21.

    González L and González-Vilar M. Determination of Relative Water Content: Handbook of Plant Ecophysiology Techniques, (Reigosa Roger MJ ed), Springer Netherlands. p^pp 207–212. 2003.

  22. 22.

    Hahn A, Kilian J, Mohrholz A, Ladwig F, Peschke F, Dautel R, et al. Plant core environmental stress response genes are systemically coordinated during abiotic stresses. Int J Mol Sci. 2013;14:7617–41.

  23. 23.

    Heitz T, Widemann E, Lugan R, Miesch L, Ullmann P, Desaubry L, et al. Cytochromes P450 CYP94C1 and CYP94B3 catalyze two successive oxidation steps of plant hormone Jasmonoyl-isoleucine for catabolic turnover. J Biol Chem. 2012;287:6296–306.

  24. 24.

    Hsieh EJ, Cheng MC, Lin TP. Functional characterization of an abiotic stress-inducible transcription factor AtERF53 in Arabidopsis thaliana. Plant Mol Biol. 2013;82:223–37.

  25. 25.

    Huang S, Gao Y, Liu J, Peng X, Niu X, Fei Z, et al. Genome-wide analysis of WRKY transcription factors in Solanum lycopersicum. Mol Genet Genomics. 2012;287:495–513.

  26. 26.

    Iwaki T, Guo L, Ryals JA, Yasuda S, Shimazaki T, Kikuchi A, et al. Metabolic profiling of transgenic potato tubers expressing Arabidopsis dehydration response element-binding protein 1A (DREB1A). J Agric Food Chem. 2013;61:893–900.

  27. 27.

    Jefferson RA. The GUS reporter gene system. Nature. 1989;342:837–8.

  28. 28.

    Jiang Y, Deyholos MK. Functional characterization of Arabidopsis NaCl-inducible WRKY25 and WRKY33 transcription factors in abiotic stresses. Plant Mol Biol. 2009;69:91–105.

  29. 29.

    Journot-Catalino N, Somssich IE, Roby D, Kroj T. The transcription factors WRKY11 and WRKY17 act as negative regulators of basal resistance in Arabidopsis thaliana. Plant Cell. 2006;18:3289–302.

  30. 30.

    Kang Y, Han Y, Torres-Jerez I, Wang M, Tang Y, Monteros M, et al. System responses to long-term drought and re-watering of two contrasting alfalfa varieties. Plant J. 2011;68:871–89.

  31. 31.

    Kasuga M, Miura S, Shinozaki K, Yamaguchi-Shinozaki K. A combination of the Arabidopsis DREB1A gene and stress-inducible rd29A promoter improved drought- and low-temperature stress tolerance in tobacco by gene transfer. Plant Cell Physiol. 2004;45:346–50.

  32. 32.

    Kilian J, Peschke F, Berendzen KW, Harter K, Wanke D. Prerequisites, performance and profits of transcriptional profiling the abiotic stress response. Biochim Biophys Acta. 2012;1819:166–75.

  33. 33.

    Kilian J, Whitehead D, Horak J, Wanke D, Weinl S, Batistic O, et al. The AtGenExpress global stress expression data set: protocols, evaluation and model data analysis of UV-B light, drought and cold stress responses. Plant J. 2007;50:347–63.

  34. 34.

    Knapp S, Bohs L, Nee M, Spooner DM. Solanaceae–a model for linking genomics with biodiversity. Comp Funct Genomics. 2004;5:285–91.

  35. 35.

    Kolodyazhnaya YS, Kutsokon NK, Levenko BA, Syutikova OS, Rakhmetov DB, Kochetov AV. Transgenic plants tolerant to abiotic stresses. Cytol Genet. 2009;43:132–49.

  36. 36.

    Kushiro T, Okamoto M, Nakabayashi K, Yamagishi K, Kitamura S, Asami T, et al. The Arabidopsis cytochrome P450 CYP707A encodes ABA 8'-hydroxylases: key enzymes in ABA catabolism. EMBO J. 2004;23:1647–56.

  37. 37.

    Langridge P, Fleury D. Making the most of ‘omics’ for crop breeding. Trends Biotechnol. 2011;29:33–40.

  38. 38.

    Li JB, Luan YS, Jin H. The tomato SlWRKY gene plays an important role in the regulation of defense responses in tobacco. Biochem Biophys Res Commun. 2012;427:671–6.

  39. 39.

    Ling MHT, Rabara RC, Tripathi P, Rushton PJ, Ge X. Extending MapMan ontology to tobacco for visualization of gene expression. Dataset Pap Biol. 2013;2013:7.

  40. 40.

    Mir RR, Zaman-Allah M, Sreenivasulu N, Trethowan R, Varshney RK. Integrated genomics, physiology and breeding approaches for improving drought tolerance in crops. Theor Appl Genet. 2012;125:625–45.

  41. 41.

    Mittler R. Abiotic stress, the field environment and stress combination. Trends Plant Sci. 2006;11:15–9.

  42. 42.

    Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K. AP2/ERF family transcription factors in plant abiotic stress responses. Biochim Biophys Acta. 2012;1819:86–96.

  43. 43.

    Navrot N, Collin V, Gualberto J, Gelhaye E, Hirasawa M, Rey P, et al. Plant glutathione peroxidases are functional peroxiredoxins distributed in several subcellular compartments and regulated during biotic and abiotic stresses. Plant Physiol. 2006;142:1364–79.

  44. 44.

    Oliver MJ, Guo L, Alexander DC, Ryals JA, Wone BWM, Cushman JC. A sister group contrast using untargeted global metabolomic analysis delineates the biochemical regulation underlying desiccation tolerance in Sporobolus stapfianus. Plant Cell. 2011;23:1231–48.

  45. 45.

    Perera IY, Hung CY, Moore CD, Stevenson-Paulik J, Boss WF. Transgenic Arabidopsis plants expressing the type 1 inositol 5-phosphatase exhibit increased drought tolerance and altered abscisic acid signaling. Plant Cell. 2008;20:2876–93.

  46. 46.

    Pilon-Smits EAH, Terry N, Sears T, Kim H, Zayed A, Hwang S, et al. Trehalose-producing transgenic tobacco plants show improved growth performance under drought stress. J Plant Physiol. 1998;152:525–32.

  47. 47.

    Rabara RC, Tripathi P, Lin J, Rushton PJ. Dehydration-induced WRKY genes from tobacco and soybean respond to jasmonic acid treatments in BY-2 cell culture. Biochem Biophys Res Commun. 2013;431:409–14.

  48. 48.

    Rabara RC, Tripathi P, Rushton PJ. The Potential of Transcription Factor-Based Genetic Engineering in Improving Crop Tolerance to Drought. Omics: a journal of integrative biology. 2014.

  49. 49.

    Raghavendra AS, Gonugunta VK, Christmann A, Grill E. ABA perception and signalling. Trends Plant Sci. 2010;15:395–401.

  50. 50.

    Rigano MM, De Guzman G, Walmsley AM, Frusciante L, Barone A. Production of pharmaceutical proteins in solanaceae food crops. Int J Mol Sci. 2013;14:2753–73.

  51. 51.

    Rushton PJ, Bokowiec MT, Han S, Zhang H, Brannock JF, Chen X, et al. Tobacco transcription factors: novel insights into transcriptional regulation in the Solanaceae. Plant Physiol. 2008;147:280–95.

  52. 52.

    Rushton PJ, Bokowiec MT, Laudeman TW, Brannock JF, Chen X, Timko MP. TOBFAC: the database of tobacco transcription factors. BMC Bioinformatics. 2008;9:53.

  53. 53.

    Rushton PJ, Reinstadler A, Lipka V, Lippok B, Somssich IE. Synthetic plant promoters containing defined regulatory elements provide novel insights into pathogen- and wound-induced signaling. Plant Cell. 2002;14:749–62.

  54. 54.

    Rushton PJ, Somssich IE, Ringler P, Shen QJ. WRKY transcription factors. Trends Plant Sci. 2010;15:247–58.

  55. 55.

    Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.

  56. 56.

    Sarkinen T, Bohs L, Olmstead RG, Knapp S. A phylogenetic framework for evolutionary study of the nightshades (Solanaceae): a dated 1000-tip tree. BMC Evol Biol. 2013;13:214.

  57. 57.

    Sears MT, Zhang H, Rushton PJ, Wu M, Han S, Spano AJ, et al. NtERF32: a non-NIC2 locus AP2/ERF transcription factor required in jasmonate-inducible nicotine biosynthesis in tobacco. Plant Mol Biol. 2014;84:49–66.

  58. 58.

    Seo JS, Sohn HB, Noh K, Jung C, An JH, Donovan CM, et al. Expression of the ArabidopsisAtMYB44 gene confers drought/salt-stress tolerance in transgenic soybean. Mol Breed. 2012;29:601–8.

  59. 59.

    Shao HB, Chu LY, Jaleel CA, Manivannan P, Panneerselvam R, Shao MA. Understanding water deficit stress-induced changes in the basic metabolism of higher plants - biotechnologically and sustainably improving agriculture and the ecoenvironment in arid regions of the globe. Crit Rev Biotechnol. 2009;29:131–51.

  60. 60.

    Sierro N, Battey JN, Ouadi S, Bakaher N, Bovet L, Willig A, et al. The tobacco genome sequence and its comparison with those of tomato and potato. Nat Commun. 2014;5:3833.

  61. 61.

    Silvente S, Sobolev AP, Lara M. Metabolite adjustments in drought tolerant and sensitive soybean genotypes in response to water stress. PLoS One. 2012;7:e38554.

  62. 62.

    Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.

  63. 63.

    Tarczynski MC, Jensen RG, Bohnert HJ. Expression of a bacterial mtlD gene in transgenic tobacco leads to production and accumulation of mannitol. Proc Natl Acad Sci U S A. 1992;89:2600–4.

  64. 64.

    Tarczynski MC, Jensen RG, Bohnert HJ. Stress protection of transgenic tobacco by production of the osmolyte mannitol. Science. 1993;259:508–10.

  65. 65.

    Tripathi P, Rabara RC, Rushton PJ. A systems biology perspective on the role of WRKY transcription factors in drought responses in plants. Planta. 2014;239:255–66.

  66. 66.

    Urano K, Maruyama K, Ogata Y, Morishita Y, Takeda M, Sakurai N, et al. Characterization of the ABA-regulated global responses to dehydration in Arabidopsis by metabolomics. Plant J. 2009;57:1065–78.

  67. 67.

    Ursini F, Bindoli A. The role of selenium peroxidases in the protection against oxidative damage of membranes. Chem Phys Lipids. 1987;44:255–76.

  68. 68.

    Vogeli U, Chappell J. Induction of sesquiterpene cyclase and suppression of squalene synthetase activities in plant cell cultures treated with fungal elicitor. Plant Physiol. 1988;88:1291–6.

  69. 69.

    Wang Y, Jiang J, Zhao X, Liu G, Yang C, Zhan L. A novel LEA gene fromTamarix androssowii confers drought tolerance in transgenic tobacco. Plant Sci. 2006;171:655–62.

  70. 70.

    Wei W, Zhang Y, Han L, Guan Z, Chai T. A novel WRKY transcriptional factor from Thlaspi caerulescens negatively regulates the osmotic stress tolerance of transgenic tobacco. Plant Cell Rep. 2008;27:795–803.

  71. 71.

    Zhang JY, Cruz DECMH, Torres-Jerez I, Kang Y, Allen SN, Huhman DV, et al. Global reprogramming of transcription and metabolism in Medicago truncatula during progressive drought and after rewatering. Plant Cell Environ. 2014

  72. 72.

    Zimmermann P, Hirsch-Hoffmann M, Hennig L, Gruissem W. GENEVESTIGATOR. Arabidopsis microarray database and analysis toolbox. Plant Physiol. 2004;136:2621–32.

  73. 73.

    Zuckerkandl E, Pauling L. Molecules as documents of evolutionary history. J Theor Biol. 1965;8:357–66.

Download references


We wish to thank John Ryals at Metabolon Inc., and Naveen Kumar, Mani Kant Choudhary, Malini Rao, Marissa Miller and Nikhil Kesarla in the Rushton lab. This project was supported by National Research Initiative grants 2008-35100-04519 and 2008-35100-05969 from the USDA National Institute of Food and Agriculture.

Author information

Correspondence to Paul J Rushton.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

RCR, PT, and PJR conceived the experiments, generated data, analyzed the data and drafted the manuscript. DA and DLR analyzed the data and drafted the manuscript. RNR conceived the experiments, generated data, and analyzed the data. MPT and QJS conceived the experiments and drafted the manuscript. All authors edited and approved the final manuscript.

Roel C Rabara and Prateek Tripathi contributed equally to this work.

Additional files

Additional file 1: Table S1.

Changes in metabolites in tobacco root and leaf tissue during drought.

Additional file 2: Figure S1.

Principal Component Analysis (PCA) plots of all independent data points.

Additional file 3: Table S2.

Genomic survey sequence reads from the Tobacco Genome Initiative project and TOBFAC transcription factors present on the oligo array together with their identification numbers on the oligo array.

Additional file 4: Table S3.

Version 3.0 DFCI Tobacco Gene Index EST sequences that were present on the oligo array.

Additional file 5: Table S4.

Genes that are up- or down-regulated at least 8-fold during drought.

Additional file 6: Figure S2.

Roles for the plant hormones ABA, JA, and ethylene in the response to drought stress in tobacco.

Additional file 7: Table S5.

Early and late up-regulated transcription factor genes during drought.

Additional file 8: Figure S3.

MapMan visualization of changes in transcription factor mRNA levels during drought stress in tobacco.

Additional file 9: Table S6.

The WRKY transcription factor family in soybean together with apparent pseudogenes.

Additional file 10: Figure S4.

Protein architecture of selected Group IX AP2/ERF transcription factors as determined by MEME.

Additional file 11: Table S7.

Drought inducibility of Subgroup 3 Group IX ERF transcription factor genes from tobacco and soybean.

Additional file 12: Table S8.

MUSCLE alignment used by MEGA6 to create the Neighbor-Joining phylogenetic tree of the WRKY family.

Additional file 13: Table S9.

MUSCLE alignment used by MEGA6 to create the Neighbor-Joining phylogenetic tree of the AP2/ERF family.

Additional file 14: Table S10.

MUSCLE alignment used by MEGA6 to create the Neighbor-Joining phylogenetic tree of part of the homeodomain family.

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Rabara, R.C., Tripathi, P., Reese, R.N. et al. Tobacco drought stress responses reveal new targets for Solanaceae crop improvement. BMC Genomics 16, 484 (2015).

Download citation


  • Drought Stress
  • Relative Water Content
  • WRKY Gene
  • Drought Stress Response
  • Gaba Shunt