Transcriptome analysis reveals absence of unintended effects in drought-tolerant transgenic plants overexpressing the transcription factor ABF3

Abstract

Background

Plants engineered for abiotic stress tolerance may soon be commercialized. The engineering of these plants typically involves the manipulation of complex multigene networks and may therefore have a greater potential to introduce pleiotropic effects than the simple monogenic traits that currently dominate the plant biotechnology market. While research on unintended effects in transgenic plant systems has been instrumental in demonstrating the substantial equivalence of many transgenic plant systems, it is essential that such analyses be extended to transgenic plants engineered for stress tolerance. Drought-tolerant Arabidopsis thaliana were engineered through overexpression of the transcription factor ABF3 in order to investigate unintended pleiotropic effects. In order to eliminate position effects, the Cre/lox recombination system was used to create control plant lines that contain identical T-DNA insertion sites but with the ABF3 transgene excised. This additionally allowed us to determine if Cre recombinase can cause unintended effects that impact the transcriptome.

Results

Microarray analysis of control plant lines that underwent Cre-mediated excision of the ABF3 transgene revealed only two genes that were differentially expressed in more than one plant line, suggesting that the impact of Cre recombinase on the transcriptome was minimal. In the absence of drought stress, overexpression of ABF3 had no effect on the transcriptome, but following drought stress, differences were observed in the gene expression patterns of plants overexpressing ABF3 relative to control plants. Examination of the functional distribution of the differentially expressed genes revealed strong similarity indicating that unintended pathways were not activated.

Conclusions

The action of ABF3 is tightly controlled in Arabidopsis. In the absence of drought stress, ectopic activation of drought response pathways does not occur. In response to drought stress, overexpression of ABF3 results in a reprogramming of the drought response, which is characterized by changes in the timing or strength of expression of some drought response genes, without activating any unexpected gene networks. These results illustrate that important gene networks are highly regulated in Arabidopsis and that engineering stress tolerance may not necessarily cause extensive changes to the transcriptome.

Background

Drought is a major abiotic stress that limits crop productivity [1]. Climate change models predict an increase in summer drying in the midlatitudes, which could contribute to an increase in the number of episodes of drought [2, 3]. Engineering plants with enhanced tolerance of abiotic stresses such as drought is a major objective of plant biotechnology that is expected to be commercialized in the near future [4, 5]. Tolerance to abiotic stress may be achieved through the modification of endogenous plant pathways, often by manipulating important regulatory proteins such as transcription factors. Altering the level of expression of key transcription factors involved in abiotic stress pathways has been shown to enhance tolerance to various abiotic stresses in Arabidopsis [6–9] as well as in important crop species such as rice [10–12], maize [13], and alfalfa [14].

Traits involving tolerance to abiotic stresses are considered to be more complex than those that are currently commercialized due to the large number of genes and pathways that may be affected. Furthermore, the interaction between plants and the environment is an intricate, continuous process that has been difficult to characterize, further adding to the complexity of manipulating abiotic stress tolerance traits. The increased complexity of these traits may correspond with a greater potential for unintended effects to occur in transgenic plants.

In transgenic systems, two different types of unintended effects are generally known to occur [15]. Position effects are attributed to the insertion of a transgene at a particular locus in the genome and the resulting interference this might cause. These effects will vary with the site of integration and will therefore be unique to each independent plant line. Position effects can be easily eliminated by screening for plant lines that have no or little position effects. In contrast, pleiotropic effects are independent of the site of transgene insertion and are the sum of all the phenotypic effects caused by expression of the transgene. While some of these may be the intended trait, others may occur through unexpected interactions of the gene with plant processes and constitute the unintended pleiotropic effects. These effects are of greater interest since they are more difficult to eliminate and more likely to create safety issues.

Engineering more complex traits such as abiotic stress tolerance in plants through the manipulation of transcription factors may uncover cryptic properties of the transcription factor that could produce some of the unintended pleiotropic effects. Many transcription factors are part of large families that have complex evolutionary histories [16, 17]. These families typically arise through gene duplications followed by functional divergence in separate expression domains or through the acquisition of new functions. These processes often result in functional redundancies within the families that can be difficult to detect. Furthermore, some transcription factors may retain ancestral functions that are sometimes only revealed by altering the normal pattern of expression. Therefore the manipulation of transcription factors in engineering complex traits such as abiotic stress tolerance may be likely to produce unintended pleiotropic effects.

The use of non-targeted global profiling technologies, such as microarray analysis, to identify unintended effects in plant systems has proven an effective means of determining the "substantial equivalence" of a transgenic plant to its non-transgenic counterpart. Such approaches have been used to investigate unintended effects in a number of transgenic plant systems [18–27]. To date, these studies have primarily focused on simple, monogenic traits such as those that are currently commercially grown. As transgenic crops with more complex traits involving the modification of endogenous plant pathways will soon be entering the market, it is important to extend these analyses to investigate the potential for unintended pleiotropic effects in such systems.

In order to understand the extent and kinds of unintended effects that could be induced in transgenic plants engineered for complex traits, we conferred drought tolerance on Arabidopsis thaliana by overexpressing the transcription factor ABF3. This system targets drought resistance, a trait that will likely enter the market in the near future. Since transcription factors ultimately function by altering the levels of expression of target genes, we investigated unintended effects using microarray analysis to survey global gene expression profiles. In order to eliminate position effects in our analysis and focus on the pleiotropic unintended effects, we employed the Cre/lox system to excise the ABF3 transgene from the site of insertion, leaving behind the selectable marker, to create control plant lines. Without the ABF3 transgene, the pleiotropic effects will be absent but the site of integration is still interrupted by the selectable marker such that position effects are maintained in these lines.

ABF3 belongs to the ABF/AREB subfamily of bZIP transcription factors which consists of thirteen members in Arabidopsis. Several members have been shown to function in ABA signalling either during seed maturation or in response to stress [28]. These factors can bind to ABA-response elements (ABREs), cis-regulatory elements found in the promoters of many ABA- and stress-responsive genes [29–31]. In addition to drought tolerance, overexpression of ABF3 confers tolerance to salt, cold, heat, and oxidative stresses, suggesting that it regulates multiple abiotic stress pathways in Arabidopsis[7, 32]. Three other ABF/AREB transcription factors are predicted to function in ABA-dependent stress signalling based on expression profiling and overexpression studies. Expression of ABF1 is induced by cold-treatment [29]. ABF2/AREB1 is induced by salt-treatment as well as dehydration but not cold and overexpression confers tolerance to a wide range of abiotic stresses, including salt, drought, heat, and oxidative stress [29, 31, 33]. Interestingly, ABF2/AREB1 also appears to function in glucose signalling as well as in the regulation of seedling growth [33]. ABF4/AREB2 is expressed in response to cold, drought, and salt and overexpression renders plants tolerant to drought and salt [7, 29, 31]. Therefore, ABF3 likely shares some redundant functions with other members of the ABF/AREB subfamily.

Plant response to drought involves changes in the expression patterns of a large number of genes [34–36] and, in addition to members of the ABF/AREB family, a number of other transcription factors have been identified that play a role in the drought response in Arabidopsis. These include the AP2/ERF transcription factors DREB2A, DREB2B, and CBF4 [8, 37, 38], AtMYB2, which functions in concert with AtMYC2 [39–41], and the NAC family transcription factors ANAC19, ANAC055/ATNAC3, and ANAC072/RD26 [42], which at least partially function in concert with a zinc finger homeodomain protein ZFHD1 [43]. Therefore, while overexpression of ABF3 affects one of the key drought response pathways, it is not the only pathway mediating the drought response at the gene expression level.

The use of the Cre/lox system to create control lines also created an opportunity to examine the effects of the Cre/lox system on the transcriptome. Site-specific recombination technologies can be used to excise selectable markers or other undesirable genetic elements and can also be used to direct site-specific integration of transgenes [44]. While many studies have employed Cre-mediated recombination in plant systems with no apparent unintended effects [45–48], other studies have observed a range of abnormal phenotypes including growth defects, leaf chlorosis, delayed flowering, and male sterility [49, 50]. In tobacco plants transformed with a chloroplast targeted Cre recombinase, recombination was observed involving cryptic lox sites in the plastid genome, but invariably the second lox site was located within the transgene [51–53]. These recombinations could result in deletions of up to 147 kb, but they did not cause any deleterious effects in the plants [51–53]. Studies in animal systems have similarly revealed that Cre recombinase can have unintended effects, often leading to chromosomal aberrations [54–57].

These studies suggest that the Cre/lox system has the potential to cause unintended effects in plant systems, by mediating recombination with cryptic lox sites that may be present in the genome resulting in large deletions. Such cryptic lox sites are difficult to identify since they may deviate substantially from conventional loxP sites [52] and not all of the unintended effects may produce readily apparent phenotypic abnormalities, so studying the unintended effects of Cre recombinase using a non-targeted approach such as microarray analysis is essential for establishing the utility and safety of this technology.

In this study, we performed microarray analysis on Arabidopsis plants engineered to be drought-tolerant through overexpression of the transcription factor ABF3 with the goal of identifying unintended pleiotropic effects. The results suggest that overexpression of ABF3 has a minimal impact on the transcriptome, with differences in the gene expression pattern only detectable in response to drought and then being suggestive of transcriptional reprogramming as opposed to the activation of novel pathways. In addition, we examined the impact of Cre recombinase on the transcriptome to detect any unintended effects of this technology and found that it had minimal effects on gene expression patterns in plants following transgene excision.

Results

Phenotype of plants overexpressing ABF3

A construct was created containing the CaMV 35S promoter followed by the coding sequence of ABF3 and then the nopaline synthase transcriptional terminator, all of which was flanked by two loxP sites (Figure 1A). This construct was transformed into Arabidopsis thaliana to generate 35S:ABF3 plants. Fifty-nine independent transformants were recovered and three high-expressing lines containing single insertions were selected for further analysis (35S:ABF3-48, -57, and -59).

Figure 1

35S:ABF3 and control T-DNA regions. Schematic map of the 35S:ABF3 and control T-DNA regions. ABF3, coding region of ABF3 gene; 35S, CaMV 35S promoter; nos, nopaline synthase transcriptional terminator; bar, bialaphos resistance; lox, loxP sites.

Control plant lines were generated by excising the 35S:ABF3 transgene, leaving only the selectable marker transgene at the site of insertion (Figure 1B). This was achieved by crossing 35S:ABF3 plants with plants expressing the Cre recombinase gene and then backcrossing the progeny to wild-type Arabidopsis to remove the Cre recombinase gene, generating three control lines (Control-48, -57, and -59). The loss of both the ABF3 transgene and the Cre recombinase transgene was confirmed by PCR (Additional file 1).

Previously, Arabidopsis plants overexpressing ABF3 were found to be more tolerant to drought conditions, which could at least partially be attributed to a lower rate of transpiration [7]. To confirm that the 35S:ABF3 lines show a similar phenotype, the transpiration rate was determined. Excised leaves from 4-week-old 35S:ABF3 plants left at room temperature for 1 day lost less fresh weight than both wild-type and control plants, indicating a lower rate of transpiration (Figure 2A). In addition, 35S:ABF3 plants showed a mild growth retardation that became more evident as plants matured (Figure 2B, C), which is also consistent with previous results [7].

Figure 2

Phenotype of 35S:ABF3 plants. (A) Transpiration rate of leaves from 35S:ABF3, control, and wild-type plants. Leaves from 4-week-old plants were excised and left at room temperature for one day and the percentage of leaf weight remaining was measured (n = 5). For each measurement, five leaves were pooled. Error bars indicate standard deviation. (B, C) Growth retardation of 35S:ABF3 plants. Three-day-old seedlings were transplanted to fresh MS plates and fresh weight was measured after 1 (B) (n = 6) or 4 weeks (C) (n = 5) of growth. For each measurement, eight (B) or 5 (C) seedlings were pooled. Error bars indicate standard deviation.

Identification of position effects and impact of Cre recombinase on the transcriptome

Microarray analysis of control plant lines was performed to identify position effects as well as to determine if use of the Cre/lox system can cause unintended effects in plant systems. Position effects should be specific to each of the three control lines, as each line is expected to have a unique insertion site. In contrast, unintended effects resulting from excision at cryptic lox sites are more likely to affect independent plant lines similarly, since cryptic sites that may be present in the Arabidopsis genome will equally be found in each plant line.

Using a P-value cut-off of 0.05, only a small number of genes were differentially expressed in control lines compared to wild-type plants (Table 1). In the Control-48 line, only a single gene was differentially expressed. In the Control-57 and -59 lines, 10 and 4 genes were differentially expressed and two of these, rps7 (AtCg00900/AtCg01240) and rps12.1 (transplice part 1 of 2) (AtCg00065), were common to both lines. These genes could represent either position effects or natural background variation in gene expression. The rps7 and rps12.1 genes that were differentially expressed in both Control-57 and -59 lines may reflect unintended effects of the Cre/lox system.

Table 1 Genes differentially expressed (P < 0.05) in control lines following Cre-mediated excision of the ABF3 transgene

Effect of ABF3 overexpression on the transcriptome in the absence of drought

Microarray analysis of 35S:ABF3 plants was performed to identify unintended effects resulting from overexpression of a transcription factor in the absence of stress. Expression of ABF3 in the absence of stress is generally very low but it is rapidly induced in response to ABA [7, 29]. Overexpression of ABF3 would therefore be expected to initiate those pathways that are typically activated in response to stress via ABA-mediated signalling. This has been observed for other transcription factors and is generally used as a means of identifying targets of that particular transcription factor [39, 42, 58].

Three independent plant lines containing the 35S:ABF3 transgene were compared to the corresponding control plant lines from which the 35S:ABF3 transgene was excised. Based on microarray analysis, using a cut-off of P < 0.05, only a small number of genes (7, 1, and 8) were differentially expressed in the three 35S:ABF3 transgenic plant lines (Table 2). The only gene that was differentially expressed in all three plant lines was ABF3.

Table 2 Genes differentially expressed (P < 0.05) in 35S:ABF3 transgenic plants

Drought response of 35S:ABF3 plants

Overexpression of ABF3 in Arabidopsis confers enhanced drought tolerance [7]. Since ABF3 is a transcription factor, it is likely that this is achieved at the level of gene expression. Therefore, microarray analysis should provide insight into the mechanism of drought tolerance. Since the expression pattern of ABF3 during the drought response is achieved by constitutive overexpression from the CaMV 35S promoter, it is possible that unintended effects may also be generated downstream. By comparing the transcriptional profile of 35S:ABF3 plants to control plants, it may be possible to identify these unintended effects. Two time points were examined in order to consider both early and late responses to drought. For this analysis, the Control-48 and 35S:ABF3-48 lines were selected as the Control-48 lines showed the smallest number of differentially expressed genes (Table 1).

Following 2 h of drought stress, using a cut-off of P < 0.05 and a fold-change of 2, 508 genes were differentially expressed in the control plants exposed to drought stress compared to unstressed control plants and 468 genes were differentially expressed in 35S:ABF3 plants exposed to drought stress compared to unstressed 35S:ABF3 plants. Following 24 h of drought stress, 592 genes were differentially expressed in control plants exposed to drought stress compared to unstressed control plants and 797 genes were differentially expressed in 35S:ABF3 plants exposed to drought stress compared to unstressed 35S:ABF3 plants.

In order to compare the profile of genes expressed in the two plant lines in response to drought stress, functional categorizations of the differentially expressed genes were obtained from The Arabidopsis Information Resource (TAIR; http://www.arabidopsis.org) and compared. Although the profile of genes differentially expressed in control and 35S:ABF3 plants show some differences, overall they show a similar pattern of distribution among the different functional classes (Additional file 2). The greatest difference between 35S:ABF3 and control plant lines in the percentage of genes belonging to a particular functional category was 1.3% for the 'other metabolic processes' category at the 2 h time point. This suggests that the overall functional response of 35S:ABF3 and control plant lines at the gene expression level was similar.

In total, 1234 genes were differentially expressed in at least one plant line during at least one time point. The overlap in genes expressed in the two plant lines at the two different time points is depicted with a four-way Venn diagram in Figure 3. These genes can be subdivided into three categories. There are 564 genes that were differentially expressed in both control and 35S:ABF3 lines at the same time points suggesting that they were commonly regulated in both lines. There are 407 genes that are regulated differently in control and 35S:ABF3 plant lines that show an enhanced response in the 35S:ABF3 line. Finally, there are 263 genes that are regulated differently in control and 35S:ABF3 plant lines that show an attenuated response in the 35S:ABF3 line. In the latter two categories, these genes are either uniquely differentially expressed in one line or the other, or they are differentially expressed in one line at one time point but not in the other line at that time point. Since only two time points were examined, it is difficult to determine if the observed differences in the two lines are due to differences in the timing of gene expression or in the magnitude of gene expression or some combination of both of these factors. Additional file 3 contains a list of all of the differentially expressed genes found in each of the three categories.

Figure 3

Overlap of differentially expressed genes in 35S:ABF3 and control plants in response to drought. Four-way Venn diagram showing the overlap of differentially expressed genes in 35S:ABF3 and control plant lines in response to drought at 2 h and 24 h. Regions corresponding to genes that are regulated commonly in both plant lines are shaded in purple. Regions corresponding to genes that show enhanced regulation in 35S:ABF3 plants are shaded in red and those that show attenuated regulation in 35S:ABF3 plants are shaded in blue.

In order to confirm the microarray results, RT-PCR was performed on 32 genes from the enhanced and attenuated categories. Twenty-nine of the examined genes exhibited expression patterns that were consistent with the microarray results, confirming the reliability of the microarray data (Figure 4). Amplification of one gene, At2g22760, produced two bands, one corresponding to the expected size of 526 bp. The primers for this gene were designed around an intron and the size of the second band is similar to the 664 bp that would be expected if the intron were not spliced, suggesting that this band may represent a splice variant.

Figure 4

RT-PCR confirmation of microarray data. RT-PCR was performed using RNA from unstressed (-) or drought stressed (+) control and 35S:ABF3 plants at either the 2 h or 24 h time point. Genes showing either enhanced or attenuated regulation in 35S:ABF3 plants in response to drought were amplified using gene-specific primers. Elongation factor 1-α (EF-1) was amplified as an internal control.

Genes commonly regulated in both 35S:ABF3 and control lines

There are 564 genes that are commonly regulated in both 35S:ABF3 and control lines. Of these, 172 show differential expression in both lines at both time points, 120 are only differentially expressed at 2 h and 272 are only differentially expressed at 24 h.

A number of genes in this category are known to act in pathways that are upstream or independent of ABF3, which is consistent with their common pattern of regulation in both 35S:ABF3 and control plant lines (Table 3). Since ABF3 is an ABA-dependent transcription factor, ABA biosynthesis should occur upstream of ABF3 activity. Consistent with this, AtNCED3 (At3g14440), which encodes a 9-cis-epoxycarotenoid dioxygenase enzyme involved in ABA biosynthesis, was upregulated in both lines at both time points in response to drought treatment as was the gene CYP707A3 (At5g45340), which encodes a cytochrome P450 monooxygenase involved in ABA catabolism.

Table 3 Genes commonly regulated in 35S:ABF3 and control plant lines that act in pathways upstream or independent from ABF3

Several transcription factors have been identified that mediate drought response pathways independent from the ABF3 pathway. Members of the DREB family of transcription factors function in abiotic stress signalling and members of the DREB2 subfamily are known to function in ABA-independent drought stress signalling [59]. DREB2A (At5g05410) and DREB2B (At3g11020), as expected, were upregulated in response to drought in both plant lines at both time points (Table 3). In addition, DREB1B/CBF1 (At4g25490) and DREB1C/CBF2 (At4g25470) were also upregulated in both lines at both time points and DDF1 (At1g12610) was upregulated at 24 h (Additional file 3). While the DREB1 subfamily is primarily associated with cold stress signalling [59], there is some evidence that members also function in dehydration stress [8, 38, 60, 61].

Three NAC transcription factors also function in a drought signalling pathway that is independent from ABF3 [42]. ANAC19 (At1g52890), ANAC055/ATNAC3 (At3g15500), and ANAC072/RD26 (At4g27410) are all expressed at both time points in both plant lines (Table 3). Another six members of the NAC family of transcription factors (At1g01720, At1g69490, At1g77450, At2g02450, At3g49530, and At5g63790) are also commonly regulated in both 35S:ABF3 and control plant lines (Additonal file 3).

Also in this category are many genes that are known to be regulated in response to drought stress. This includes genes involved in the biosynthesis of osmolytes, late embryogenesis abundant proteins (LEA), kinases, phosphatases, and transcription factors, as well as several different types of transporters (Table 4). This further suggests that some drought signalling pathways are unaffected by ABF3 overexpression

Table 4 Examples of typical drought-responsive genes regulated similarly in 35S:ABF3 and control plant lines

Genes with an enhanced response in 35S:ABF3 plants

There are 407 genes that demonstrated an enhanced response in 35S:ABF3 plants. Of these, 332 are uniquely regulated in 35S:ABF3 plants with 9 differentially expressed at both time points, 87 differentially expressed only at 2 h and 236 differentially expressed at 24 h. In addition, 47 genes are differentially expressed in the control line at only 24 h and 27 genes are differentially expressed in the control line at only 2 h while in the 35S:ABF3 line they are differentially expressed at both time points. Finally, there is a transposable element gene that is differentially expressed in 35S:ABF3 plants at 2 h but is only differentially expressed in control plants at 24 h.

A number of genes in this category could contribute to the enhanced drought tolerance of 35S:ABF3 plants (Table 5). Several genes involved in the biosynthesis of the osmolyte trehalose show enhanced upregulation while two genes involved in osmolyte catabolism show enhanced downregulation. Two genes involved in detoxification show enhanced upregulation in 35S:ABF3 as do three LEA proteins. A number of transporters also show enhanced regulation in 35S:ABF3 lines.

Table 5 Genes showing an enhanced response in 35S:ABF3 plants that may contribute to drought-tolerance

Direct targets of ABF3 are likely found in this category of genes. In addition to showing enhanced regulation in 35S:ABF3 plants, it is expected that these genes possess at least one ABRE in their promoter. In addition, these genes are likely to be significantly differentially expressed at the 2 h time point. Amongst those genes showing an enhanced response in 35S:ABF3 plants that are significantly differentially expressed at the 2 h time point, 24 contain at least one ABRE according to the in silico analysis performed by Gómez-Porras et al [62] and are therefore identified as putative ABF3 targets (Table 6). Included in this group are four transcription factors, one choline kinase, one trehalose biosynthetic enzyme, one gene involved in ubiquitin-mediated protein degradation, and two transporters as well as seven other genes with undefined roles in the drought response and eight unknown genes.

Table 6 Potential targets of ABF3

Amongst those genes that show an enhanced response in 35S:ABF3 plants, there is an enrichment of genes that function in RNA processing pathways (Table 7). There are 8 genes involved in RNA processing that are similarly regulated in both control and 35S:ABF3 plants lines, 8 genes that show an attenuated response in 35S:ABF3, while 18 genes show an enhanced response. Many of these appear to be uniquely downregulated in 35S:ABF3 plants at 24 h. Eleven of these are predicted to function in RNA splicing or to be associated with the RNA splicing machinery while the others have roles in nucleocytoplasmic transport, deadenylation, or the function is not specifically known. Several genes that function in RNA processing have been found to function in ABA and abiotic stress signalling pathways [63–66]. This suggests that the differential regulation of RNA processing genes in 35S:ABF3 lines could contribute to its enhanced drought tolerance.

Table 7 Genes showing enhanced regulation in 35S:ABF3 with predicted function in RNA processing

Interestingly, a member of the DREB family of transcription factors is also included in this group. DREB1D/CBF4 (At5g51990) was upregulated in 35S:ABF3 plants at both time points but was only upregulated in control plants at 24 h. DREB1D/CBF4 functions in both cold and drought stress signalling and, unlike the DREB2 subfamily, its expression is induced by ABA [38]. This suggests that CBF4 may act downstream of ABF3 in the ABA-dependent drought signalling pathway.

Genes with an attenuated response in 35S:ABF3 plants

There are 263 genes that show an attenuated response in 35S:ABF3 plants. Of these, 224 are uniquely regulated in control lines with 3 differentially expressed at both time points, 147 differentially expressed only at 2 h and 74 differentially expressed only at 24 h. In addition, 18 genes are differentially expressed in the 35S:ABF3 line at only 24 h and 5 genes are differentially expressed in the 35S:ABF3 line at only 2 h while in the control line they are differentially expressed at both time points. Finally, there are 16 genes that are differentially expressed in the control line at 2 h and in the 35S:ABF3 line at 24 h.

Surprisingly, a number of genes that typically function in conferring drought tolerance show an attenuated response in 35S:ABF3 plants. This includes several peroxidases and glutaredoxins involved in detoxification, two genes involved in biosynthesis of the osmolytes raffinose and spermine, a heat shock protein, and a number of transporters (Additional file 3). Similarly a number of transcription factors and other signalling components also show an attenuated response in 35S:ABF3 plants. The reason for the attenuated response of these genes in 35S:ABF3 plants is not clear but might reflect differences in the physiological state of 35S:ABF3 plants due to their enhanced drought resistance.

A number of genes encoded by the chloroplast and mitochondrial genomes show an attenuated response in 35S:ABF3 plants (Table 8). Some of these genes encode proteins with electron transport activity or NADH dehydrogenase activity while others are predicted to function in transcription and translation processes. Most of these genes are upregulated in control lines only at the 2 h time point. The reason for the exclusive upregulation of these genes in the control line is not clear. The chloroplast NADH dehydrogenase (NDH) complex is predicted to be involved in cyclic electron transport around photosystem I, thereby dissipating energy and maintaining ATP supply under conditions of low CO2 availability following stomatal closure in response to stress [67]. The NDH complex may also be involved in detoxification in the chloroplast [68]. The upregulation of genes encoding NADH dehydrogenases may therefore reflect the greater sensitivity of the control plant line to stress. Rice plants overexpressing ABF3 were able to maintain higher photochemical efficiency during drought stress [11]. This suggests that one aspect of the drought tolerance conferred by ABF3 overexpression may be a minimization of the negative impact of drought stress on photosynthesis, which is reflected by differences in gene expression in the chloroplast. Similar effects may also occur in the mitochondria.

Table 8 Genes encoded by the chloroplast and mitochondrial genomes showing attenuated regulation in 35S:ABF3 plants

Four transposable element genes were also included in the attenuated response category (Table 9). Only one other transposable element gene was significantly differentially expressed and it was included in the enhanced regulation category. Many transposable elements are transcriptionally activated in response to stress conditions [69, 70]. The attenuated regulation of transposable element genes in 35S:ABF3 plants might suggest that they are experiencing a lower level of stress that is insufficient to activate the transposable element genes, consistent with the drought tolerance of these plants.

Table 9 Transposable element genes showing attenuated regulation in 35S:ABF3 plants

Interestingly, DREB1A/CBF3 showed an attenuated response in 35S:ABF3 plants. DREB1A/CBF3 was upregulated in control plant lines at both time points but was only upregulated in 35S:ABF3 plants at 2 h. This may again reflect differences in the tolerance of 35S:ABF3 and control plants to the drought stress, with control plants requiring a stronger or longer activation of DREB1A/CBF3 expression in response to the increased stress.

Discussion

The Impact of Cre recombinase on the transcriptome is minimal

In order to eliminate position effects and focus on unintended pleiotropic effects of transcription factor overexpression, the Cre/lox recombination system was employed to create a series of control plant lines that contain the selectable marker at the site of transgene insertion but from which the ABF3 transgene was excised. The use of the Cre/lox recombination system also allowed us to determine the impact of Cre recombinase on the transcriptome.

In tomato, petunia, tobacco, and to a lesser extent Arabidopsis, expression of Cre recombinase has resulted in abnormal phenotypes, including leaf chlorosis, stunted growth, and sterility [49, 50]. Similarly, expression of Cre recombinase resulted in reduced proliferation and chromosomal abnormalities in cultured embryonic mouse cells [55, 57], toxicity in dividing cells of Drosophila melanogaster[54], and chromosomal rearrangements in mouse spermatids leading to male sterility [56]. These abnormal effects of Cre recombinase are suspected to result from Cre-mediated recombination using cryptic lox sites that may be found in eukaryotic genomes. Cryptic lox sites that can be recognized by Cre recombinase have been identified in the genomes of yeast and humans as well as the chloroplast genome of tobacco [52, 53, 71, 72].

If cryptic lox sites exist in plant genomes, expression of Cre recombinase could induce deletions or inversions of genome segments or even chromosome translocations that could adversely impact the plant. These deletions or inversions may alter the expression of genes found within this segment, which would be detectable by microarray. In three independent plants lines in which Cre recombinase was employed to excise the ABF3 transgene from the T-DNA insertion, only a small number of genes were found to be differentially expressed (Table 1). Of these, two genes were found to be differentially expressed in two out of the three control lines. The two genes are the chloroplast encoded rps7 (AtCg00900 and AtCg01240) and rps12.1 (AtCg00065) genes. In tobacco, one of the cryptic lox sites identified in the chloroplast genome is found just downstream of the start site of the rps12.2/rps7 operon [52, 53]. The altered expression of rps7 might suggest that this cryptic lox site is conserved in Arabidopsis and may have undergone Cre-mediated recombination. Cre recombinase may therefore also be able to act on cryptic lox sites in the Arabidopsis chloroplast genome, resulting in a change in expression of the affected genes. It is, however, unclear how Cre recombinase is targeted to the chloroplast since it is only predicted to be targeted to the nucleus. Cre-mediated chloroplast genome deletions are not likely to be of great concern since chloroplast genomes containing deletions in essential genes are typically rapidly lost due to selection pressures [73–75], especially once the Cre recombinase has been removed. The impact of Cre recombinase on the nuclear transcriptome was negligible, which demonstrates that in Arabidopsis this technology does not produce unintended effects.

The activity of ABF3 is strictly controlled

Microarray analysis of Arabidopsis plants overexpressing the transcription factor ABF3 suggests that alterations to the transcriptome are minimal when position effects are eliminated as a source of variation. In the absence of stress, a small number of genes were differentially expressed in three 35S:ABF3 plant lines (Table 2), but no genes were differentially expressed in more than one independent line.

Members of the ABF/AREB family of transcription factors bind to ABREs found in the promoters of ABA-responsive genes [29, 31]. If the genes identified by microarray analysis are actual downstream targets of ABF3, they would be expected to contain at least one ABRE. An in silico analysis of the Arabidopsis nuclear genome has identified 3829 genes containing one or more ABREs [62]. None of the nuclear genes identified by microarray analysis of 35S:ABF3 transgenic plants are predicted to contain an ABRE. Other members of the ABF/AREB subfamily of transcription factors localize to the nucleus [76, 77], and it is likely that ABF3 similarly functions in the nucleus. Therefore, it is unlikely that the chloroplast genes identified by microarray analysis are functional targets of ABF3. This suggests that overexpression of ABF3 alone is not sufficient to alter the transcriptome.

This result was unexpected as previous work identified a number of genes with altered expression in Arabidopsis[7] and rice [11] overexpressing ABF3. Similarly, overexpression studies of many other transcription factors have revealed alterations in gene expression and this approach is typically used to identify the gene network controlled by that particular transcription factor [39, 42, 58]. The absence of differentially expressed genes in 35S:ABF3 transgenic plants suggests that an additional signal is required to activate ABF3 that is not present in unstressed plants.

There is accumulating evidence that members of the ABF/AREB family of transcription factors are regulated by phosphorylation. ABF2/AREB1 transactivation of a reporter gene in the presence of ABA was inhibited by the addition of the protein kinase inhibitor staurosporine [31] and ABI5 is phosphorylated following ABA treatment [78]. Several studies have suggested a role for members of the SnRK2 family of protein kinases in the phosphorylation of ABF/AREB transcription factors [79–84]. ABF3 and ABF4/AREB2 interact with the calcium-dependent protein kinase AtCPK32 and evidence suggests that it phosphorylates a highly conserved serine residue in ABF4/AREB2 that is necessary for activity [85]. The protein kinases CPK4 and CPK11 are also likely to phosphorylate ABF1 and ABF4/AREB2 and their activity is enhanced by ABA [86]. It is possible that in the absence of stress, ABF3 is not phosphorylated and therefore cannot activate gene expression.

Furthermore, other factors necessary for the activity of ABF3 may not be expressed in the absence of abiotic stress. Members of the ABF/AREB family have been shown to interact with diverse proteins that are predicted to modulate their transcriptional activity. ABF2/AREB1 interacts with an arm-repeat protein that is predicted to positively regulate its activity [87]. ABI5, ABF1, ABF3, and ABF4/AREB2 can interact with the transcription factor ABI3 [88, 89]. Furthermore, the rd29a promoter contains both an ABRE as well as a dehydration-responsive elements that is bound by members of the DREB/CBF family of transcription factors and the two elements function interdependently to activate expression of rd29a[90]. Members of the ABF/AREB can also heterodimerize [91], suggesting that other members of this family may need to be expressed in order for ABF3 to be functional. Therefore, it is possible that other components of the stress response pathway are necessary in order for ABF3 to be active, preventing ABF3 from altering gene expression in the absence of stress.

While the 35S:ABF3 plants did not show any changes in transcription in the absence of drought stress, there were some phenotypic differences compared to control plants. Most notably, the 35S:ABF3 plants were smaller in size than control plants of the same age, with the difference becoming more pronounced with increased age (Figure 2B, C). This is a common observation for plants overexpressing transcription factors [7, 33, 39, 60] and is often overcome by using tissue-specific or inducible promoters [32, 60]. At least some of the growth retardation may be attributable to reduced transpiration rates of 35S:ABF3 plants compared to control plants (Figure 2A), which is consistent with the observation that Arabidopsis plants overexpressing ABF3 typically have stomata with smaller openings than do wild-type plants [7]. This would suggest that ABF3 may govern gene networks involved in stomatal closure. Consistent with this, ABF3 is expressed in guard cells and its expression is further induced in these cells in response to ABA [7, 92]. Since our analysis was performed on whole plants, it is likely that changes in the transcriptional network of guard cells would not be readily detectable.

Overexpression of ABF3 results in transcriptional reprogramming of the drought response

Overexpression of ABF3 confers drought tolerance to Arabidopsis plants and since ABF3 is a transcription factor, it can be predicted that this will occur through changes to the transcriptional network of the plants. Consistent with this, the expression profile of Arabidopsis plants overexpressing ABF3 differed from that of control plants. As might be expected, there were a number of genes with expression patterns that appeared to be enhanced in 35S:ABF3 plants compared to control plants. Whether this occurred through alterations in the timing or strength of expression could not be established with the two time points considered in this study. Many of the genes with enhanced expression are known to function in mitigating drought stress, suggesting that they could contribute incrementally to the enhanced drought tolerance of 35S:ABF3 plants (Table 5). Those genes showing enhanced regulation in 35S:ABF3 plants likely include some direct targets of ABF3. In particular, those genes that are differentially expressed at 2 h and contain at least one ABRE are the most likely targets of ABF3 (Table 6).

Interestingly, there seemed to be a number of genes involved in RNA processing that showed enhanced expression in 35S:ABF3 plant lines. Most of these were downregulated at the 24 h time point. A number of RNA processing mutants impaired in ABA response have demonstrated the importance of RNA processing to ABA signalling [93]. Loss-of-function mutations in two genes encoding subunits of a nuclear cap-binding complex cause ABA hypersensitivity [63, 66] as do mutations in a pre-mRNA splicing factor that is important for both mRNA splicing and turnover [94], a poly(A)-specific ribonuclease that is predicted to function in mRNA degradation [65], and a protein with homology to the Sm-like small nuclear ribonucleoproteins (snRNPs) that function in mRNA splicing, export, and degradation [95]. Two phosphatases that belong to the family of proteins that dephosphorylate the C-terminal domain of RNA polymerase II negatively regulate stress-responsive gene transcription [96] and a mutation in one of these results in decreased ABA sensitivity [97]. The enrichment of genes involved in RNA processing that show enhanced expression in the 35S:ABF3 line might suggest that ABF3 plays an important role in the regulation of ABA-responsive RNA processing events. The downregulation of many of the RNA processing genes is consistent with the negative regulatory role observed for several of the RNA processing proteins previously identified in the ABA-signalling pathways.

In addition to those genes showing enhanced regulation in 35S:ABF3 plants, there were also a number of genes with attenuated expression. Many of these are known to encode proteins with roles in minimizing drought stress and their attenuated expression is not consistent with the drought tolerance of the 35S:ABF3 plants. These genes may be reflective of a greater transcriptional reprogramming in 35S:ABF3 plants than merely enhancing the rate or level of expression of a subset of genes. The strong activation of the ABF3 pathway may result in co-ordinated feedback that modulates other drought responsive pathways, resulting in attenuated gene expression in some cases. This might reflect cooperativity between some of the drought signalling pathways. In many cases, drought-responsive transcription factors have been shown to function in concert to activate gene expression [39, 40, 43, 90]. Furthermore, it has been observed that downregulation of the phosphoinositide pathway in Arabidopsis results in an upregulation of the DREB2A gene as well as several DREB2A-regulated genes, suggesting a negative interaction between two drought signalling pathways [98].

The degree of stress experienced by the control plants should be greater than that experienced by the 35S:ABF3 plants and this may ultimately have secondary consequences on the transcriptional network that is reflected by the genes showing an attenuated transcriptional response. This is consistent with the observation that several transposable element genes show attenuated expression in 35S:ABF3 plants (Table 9). Many transposons are activated in response to stress [69, 70] and their delayed activation in 35S:ABF3 plants could be a result of the drought tolerance of these plants.

Similarly, the enrichment of genes encoded by the chloroplast and mitochondrial genomes in the attenuated category (Table 8) may also reflect the greater drought sensitivity of the control plant lines. Drought has a significant and complex impact on the activities of both the mitochondria and chloroplasts and this is reflected in changes in gene expression [99, 100]. Several of the chloroplast and mitochondria-encoded genes showing attenuated expression encode NADH dehydrogenases. The chloroplast NDH complex is predicted to function in cyclic electron transport around photosystem I to help dissipate excess energy during abiotic stresses such as drought and thereby to alleviate oxidative stress [67]. The elevated expression of chloroplast NDH complex subunits and other chloroplast and mitochondria-encoded genes in control plants may therefore be indicative of the increased stress of these plants.

Drought stress can change the composition of a plant, altering levels of oils, proteins and other constituents, which can affect the commercial and nutritional value of the plant [101–107]. In some cases, drought can also initiate the accumulation of higher levels of dangerous toxins and anti-nutrients [102, 105, 108]. The altered expression of transposable element genes and genes that are encoded by the chloroplast and mitochondrial genomes in control plants suggests that these plants are experiencing a higher level of drought stress than 35S:ABF3 plants. It is therefore possible that control plants exposed to drought will exhibit more compositional changes compared to unstressed plants than 35S:ABF3 plants. This is an important consideration because compositional analysis is one of the parameters used to determine the substantial equivalence of transgenic plants to their non-transgenic comparators during the risk assessment process. These results might suggest that 35S:ABF3 plants are better able to maintain compositional standards under drought stress than their non-transgenic counterparts.

Overexpression of ABF3 does not activate unintended gene networks during the drought response

While differences were observed in the patterns of gene expression of 35S:ABF3 and control plant lines, at the functional level the response was very similar. Functional categorization of the significantly differentially expressed genes demonstrated that the percentage of genes in each category was relatively similar between the two plant lines (Additional file 2). This suggests that overexpression of ABF3 does not activate new gene networks but simply functions to modify existing gene networks that function in drought response. Although there were genes that were uniquely regulated in each of the two plant lines, many of these genes did observe a change in expression in the other plant line that was not of a high enough magnitude to meet our criteria for differential expression (Additional file 3). This suggests that these genes are showing stronger and/or earlier differential regulation as opposed to being uniquely regulated.

Among the genes that were differentially regulated between the 35S:ABF3 and control plant lines, there was no indication that overexpression of ABF3 activated any unintended gene networks. Several closely related members of the ABF/AREB transcription factor family function in seed germination and seed and early seedling developmental pathways, including ABI5 whose role in these pathways has been best characterized [109–111]. Although ABF3 is primarily associated with abiotic stress signalling in vegetative tissues, there is some evidence that it may also function in seed and early seedling developmental processes, although its role may be relatively minor. Microarray data from GENEVESTIGATOR [112] indicates that ABF3 is expressed during seed development, although at low levels compared to ABI5. Also, levels of ABF3 are enhanced by dehydration and salinity stresses or by ABA treatment in germinating embryos [113] and an alternative splice form of the ABF3 gene was identified from a cDNA library prepared from immature seed [91]. Double mutant analysis has also revealed several redundancies between ABF3 and ABI5, including sensitivity to ABA during germination, stress sensitivity of root growth, resistance to glucose, and regulation of the ABA-induced vegetative expression of RAB18 and RD29B[88]. ABI5 and ABF3 also appear to antagonistically cross-regulate each other [88, 114].

If ABF3 does function in seed developmental pathways, it is possible that altering its pattern of expression could ectopically activate gene networks involved in those pathways. Another member of the ABF/AREB family of transcription factors, ABF2/AREB1, which functions in abiotic stress signalling and glucose response, is also expressed in embryonic axes in dry siliques [33, 77]. Overexpression of a phosphorylated active form of ABF2/AREB1 led to the activation of several seed storage protein genes in vegetative tissues, many of which also have binding sites for ABI3, another transcription factor involved in seed development [81]. However, none of the seed storage genes activated by the phosphorylated active form of ABF2/AREB1 were differentially expressed in the 35S:ABF3 plants, nor were AtEm1 and AtEm6, two LEA-class genes that are ABI5 targets expressed during seed maturation (data not shown). This further demonstrates that overexpression of ABF3, while modifying patterns of gene expression, did not activate unintended pathways in response to drought stress in Arabidopsis. This suggests that overexpression of a transcription factor to confer an abiotic stress tolerance trait may not necessarily produce unintended pleiotropic effects.

Conclusions

The Cre/lox recombination system allowed us to create paired plant lines with identical T-DNA insertion sites either with or without the ABF3 transgene in order to eliminate position effects in our analysis of unintended effects. This approach also allowed us to examine unintended effects resulting from the expression of Cre recombinase. We found that Cre recombinase had a minimal impact on the transcriptome, which suggests that it produces few unintended effects in Arabidopsis. Microarray analysis of Arabidopsis plants overexpressing ABF3 demonstrated that the impact on the transcriptome is minimal. In the absence of drought stress, there were no differentially expressed genes. In response to drought stress, a reprogramming of the drought response was observed, suggestive of changes in the timing or strength of expression of some genes in 35S:ABF3 plants. Some of these changes may be directly related to the action of ABF3 while others may reflect an altered physiological state as a result of the enhanced drought tolerance of 35S:ABF3 plants. Amongst the differentially expressed genes, no unintended pathways appeared to be activated as a result of ABF3 overexpression. These results are significant because they demonstrate that plant responses to abiotic stresses such as drought may be strictly coordinated at multiple regulatory steps and this limits the extent of unintended pleiotropic effects. This demonstrates that engineering stress tolerance through manipulation of endogenous plant pathways may not necessarily produce unintended pleiotropic effects despite the complexity of such traits. This is an important finding for establishing the safety of such traits as they begin to enter the market in the near future.

Methods

Generation of transgenic plants

To generate 35S:ABF3 plant lines, Arabidopsis thaliana Col-0 plants were transformed by the floral dip method [115] using Agrobacterium tumefaciens strain GV3101 harbouring the pCAMBIA3300 vector containing the Cauliflower mosaic virus 35S promoter, ABF3 coding region, and nopaline synthase transcriptional terminator, bordered on either side loxP sites. Homozygous single insert transformants were identified by Southern blot and segregation analysis and three lines with high levels of ABF3 expression, as determined by RT-PCR analysis, were selected for further analysis.

To generate plant lines expressing Cre recombinase, Arabidopsis thaliana Col-0 plants were transformed as above with a pCAMBIA1200 vector containing the Cauliflower Mosaic Virus 35S promoter, Cre recombinase coding region, and nopaline synthase transcriptional terminator.

To generate control plant lines, 35S:ABF3 plants were crossed with plants expressing Cre recombinase in order to excise the 35S-ABF3-nos transgene. The F1 plants were PCR genotyped to identify those that underwent successful excision, using the primers ABF-F (5'-ATGGGGTCTAGATTAAAC-3') and Nos-R (5'-CCGATCTAGTAACATAGATG-3'), which are specific for the ABF3 gene and the nopaline synthase transcriptional terminator, respectively. One plant from each cross, Control-48; Cre 1.1.2, Control-57; Cre 1.2.2 and Control-59; Cre 2.2.1, was selected from which loss of the 35S-ABF3-nos transgene was confirmed by PCR (Additional File 1). F2 progeny from the selected F1 plants were then backcrossed to wild-type plants to eliminate Cre recombinase. The F1 plants from this cross were again PCR genotyped to identify those from which the Cre gene was lost, using primers CRE.F (5'-CCAGGCGTTTTCTGAGCATACCTG-3') and CRE.R (5'-CTCTGACCAGAGTCATCCTTAGCG-3'), which are both specific for the Cre recombinase gene. One plant from each cross, Control-48 1.1.2.5, Control-57 1.2.2.5 and Control-59 2.2.1.4, was selected from which loss of the Cre recombinase gene was confirmed by PCR (Additional File 1). The selected F1 plants represent the control plant lines and the F2 progeny of these plants were used in subsequent experiments.

Growth of Arabidopsis plants

Seeds were surface sterilized by soaking for 20 min in 25% (v/v) commercial Clorox (final concentration of 1.3% sodium hypochlorite) and 0.05% (v/v) Triton X-100 (Fisher Scientific, Hampton, NH, USA), and then rinsing four times with distilled water. Seeds were germinated on MS media and grown at 22°C with a 16 h photoperiod. To determine the growth rate of plants, three-day-old seedlings were transplanted onto fresh MS plates and after one week or four weeks of growth the fresh weight of the seedlings was measured. To measure transpiration rate, leaves of a similar developmental stage were excised from four-week-old plants and the loss of weight over a 24 h-period was measured. For microarray analysis, seeds were germinated on MS media and one-week-old seedlings were collected. For the drought stress, seedlings were transferred onto paper towels and harvested after 2 and 24 h.

Microarray Analysis

For each sample, three biological replicates were prepared. RNA was extracted from seedlings using the RNeasy Plant Mini Kit (Qiagen, Germantown, MD, USA), following the manufacturer's protocol. A 2100 Bioanalyser (Agilent Technologies, Santa Clara, CA, USA) was used to determine the quality of the RNA. Microarray analysis was performed using the GeneChip^® Arabidopsis ATH1 Genome Array (Affymetrix Inc., Santa Clara, CA, USA). Standard RNA processing, hybridization, and scanning protocols were followed as recommended by the GeneChip^® Expression Analysis Technical Manual (Affymetrix Inc., Santa Clara, CA, USA). Hybridization, and scanning were performed at Agriculture and Agri-Food Canada (Winnipeg, MB, Canada) by Mark Jordan. RMA procedure [116] was performed to normalize data using the AffylmGUI R software package from Bioconductor (http://www.bioconductor.org/)[117]. Analysis of differential expression was done using a moderated t-test with empirical Bayes smoothing [118]. Microarray data from this study have been deposited at ArrayExpress (accession number E-MEXP-2435).

Reverse transcription-PCR (RT-PCR)

DNaseI-treated RNA was used for first strand cDNA synthesis using Superscript III reverse transcriptase (Invitrogen, Carlsbad, CA, USA) and oligo (dT)₁₈ primers according to the manufacturer's protocol. Between 24 and 30 cycles of PCR amplification was performed using gene-specific primers. As an internal control, elongation factor 1-α (EF-1) was amplified. Sequences of all primers used in RT-PCR analysis can be found in Additional file 4.