Transcriptional profiling of pea ABR17 mediated changes in gene expression in Arabidopsis thaliana

Background Pathogenesis-related proteins belonging to group 10 (PR10) are elevated in response to biotic and abiotic stresses in plants. Previously, we have shown a drastic salinity-induced increase in the levels of ABR17, a member of the PR10 family, in pea. Furthermore, we have also demonstrated that the constitutive expression of pea ABR17 cDNA in Arabidopsis thaliana and Brassica napus enhances their germination and early seedling growth under stress. Although it has been reported that several members of the PR10 family including ABR17 possess RNase activity, the exact mechanism by which the aforementioned characteristics are conferred by ABR17 is unknown at this time. We hypothesized that a study of differences in transcriptome between wild type (WT) and ABR17 transgenic A. thaliana may shed light on this process. Results The molecular changes brought about by the expression of pea ABR17 cDNA in A. thaliana in the presence or absence of salt stress were investigated using microarrays consisting of 70-mer oligonucleotide probes representing 23,686 Arabidopsis genes. Statistical analysis identified number of genes which were over represented among up- or down-regulated transcripts in the transgenic line. Our results highlight the important roles of many abscisic acid (ABA) and cytokinin (CK) responsive genes in ABR17 transgenic lines. Although the transcriptional changes followed a general salt response theme in both WT and transgenic seedlings under salt stress, many genes exhibited differential expression patterns when the transgenic and WT lines were compared. These genes include plant defensins, heat shock proteins, other defense related genes, and several transcriptional factors. Our microarray results for selected genes were validated using quantitative real-time PCR. Conclusion Transcriptional analysis in ABR17 transgenic Arabidopsis plants, both under normal and saline conditions, revealed significant changes in abundance of transcripts for many stress responsive genes, as well as those related to plant growth and development. Our results also suggest that ABR17 may mediate stress tolerance through the modulation of many ABA- and CK-responsive genes and may further our understanding of the role of ABR17 in mediating plant stress responses.


Background
Pathogenesis-related (PR) proteins are part of the plant defense responses that are induced by pathogens as well as by abiotic stresses [1,2]. To date, 17 different families of PR proteins have been identified, based on their specific structural and functional properties [ [3] and references therein]. Among the PR proteins, the PR10 family is composed of intracellular acidic proteins with molecular masses ranging from 15-18 kD and are encoded by multiple genes [1,3]. PR10 genes were first described in Pisum sativum inoculated with Fusarium solani [4] but have been subsequently described in many species [reviewed in [3]]. In addition to their inducible expression in response to stresses, PR10 genes also exhibit constitutive high expressed levels in roots, flowers and pollen during normal growth and development, suggesting additional roles beyond pathogenesis responses [5].
Based on sequence similarities, PR10 proteins have been suggested to be ribonucleases (RNases) [6]. Indeed, PR10 proteins from a variety of species including two pea PR10 proteins have been demonstrated to possess RNase activity [7,8]. Although RNase activities have been detected for many PR10 proteins, they have also been shown to interact with molecules such as cytokinins (CKs), brassinosteroids, fatty acids, and flavonoids [9][10][11]. These observations have led to the suggestion that all PR10 proteins may not be RNases and may be involved during normal plant growth and development as hormone/ligand carriers [10][11][12]. This suggestion is further supported by the fact that CK-specific binding proteins (CSBPs) exhibit amino acid sequence and predicted secondary-structure similarities with PR10 proteins and, for this reason, have been included in the PR10 family [9].
The pea abscisic acid-responsive protein ABR17, induced by the exogenous application of abscisic acid (ABA) is classified as a member of the PR10 family in pea [13]. ABR17 is produced late in seed development, and is homologous to dehydrins and late embryogenesis abundant (LEA)-related proteins [14,15]. ABR17 is also significantly homologous to intracellular pathogenesis related (IPR) proteins and major birch pollen allergen Betv1 proteins [16,17]. Our previous research has demonstrated the expression of ABR17 protein in pea under salt stress [2] and the RNase activity of two members of pea PR10 proteins (PR10.1 and ABR17) [7,8]. Furthermore, we have also demonstrated that the constitutive expression of pea PR10.1 and ABR17 cDNAs enhance germination and early seedling growth under abiotic stress conditions in B. napus and A. thaliana, respectively [18,19]. In addition, the transgenic plants also exhibited phenotypic differences when compared to their WT counterparts, which included precocious flowering, a higher number of lateral branches, and increased numbers of seed pods [8]. Many of these characteristics of ABR17-transgenic A. thaliana are suggestive of a role for CKs in ABR17 action, particularly increased lateral branching and early flowering [20,21]. Our suggestion was further supported by the elevated concentrations of endogenous CKs in PR10.1 transgenic B. napus as well as ABR17-transgenic A. thaliana [7,8].
These observations led us to hypothesize that PR10 proteins, including ABR17, may mediate the observed phenotypic effects through modulation of endogenous CKs. Additional evidence supporting this hypothesis has been provided by the demonstration that exogenous application of CK enhances germination under abiotic stress conditions [8]. In order to further investigate the ABR17mediated changes in A. thaliana we investigated global changes in gene expression using microarrays. Microarray analysis was carried out in an ABR17-transgenic line compared to its WT, salt treated ABR17-transgenic line compared to untreated ABR17-transgenic line, and salt treated WT compared to untreated WT seedlings. Our current findings reveal that, even in the absence of stress, the expression of genes involved in plant growth and development are significantly (and approximately 2-fold) increased in the transgenic line. Salt treated ABR17-transgenic A. thaliana seedlings showed general salt response theme comparable to that of the WT counterpart used in this study. However, both the trend as well as the degree of changes in gene expression of many defense related genes including plant defensins and heat shock proteins was different providing additional insights into the possible ways in which ABR17 may mediate plant responses to stress.

Characterization of ABR17-transgenic plants/seedlings
The appearance of 2-week-old WT and ABR17-transgenic A. thaliana seedlings grown in soil as well as on MS medium (1.5% sucrose, 0.8% agar with pH 5.7) [22] plates are shown in Figure 1. At all growth stages investigated, the ABR17-transgenic line was considerably more developmentally advanced compared to its WT counterpart. For example, in the 5-day-old transgenic seedlings ( Figure 1A) cotyledons were more developed than in their WT counterparts and at 14-days the transgenic seedlings possessed more rosette leaves ( Figure 1A). Similar developmental differences were also observed at 21 days where many transgenic seedlings had started to bolt whereas very few (if any) WT seedlings had advanced to this developmental stage ( Figure 1B). At 28 days, the transgenic seedlings also possessed more lateral branches ( Figure 1C, Table 1). The transgenic seedlings also flowered earlier than WT with an average difference of at least 2.5 days (Table 1). Seedlings for microarray experiments were grown on semi-solid MS media in order to maintain sterility and it was evident that under these growth condi-tions also the transgenic seedlings were more developmentally advanced ( Figure 1D). These results are consistent with our previous observations of seedlings from this and other independently derived ABR17-transgenic lines grown on semi solid MS media [19]. In addition, the ABR17-transgenic seedlings grown on MS media with 100 mM NaCl were greener and their roots appeared to be longer compared to the WT seedlings grown under similar conditions (Figure 2A).
Characteristics like root length, fresh weight, chlorophyll a/b content and carotenoid contents were measured in salt treated-ABR17 and WT seedlings. Roots of ABR17transgenic seedlings were relatively longer in the absence of salt whereas upon salt treatment, the differences in lengths were not that obvious ( Table 1). The fresh weight of ABR17-transgenic seedlings was not significantly different from its WT counterpart in the absence of stress. However, in the presence of 100 mM NaCl, the fresh weights of the transgenic seedlings were significantly (p < 0.05) Appearance of WT and ABR17 transgenic A. thaliana at various growth stages: Seedlings at 5, 14 days (A), 21 days (B), 28 days (C) and MS-grown 14-day-old seedlings (D) are shown higher than their WT counterparts (Table 1). Although the chlorophyll and carotenoid contents were almost similar in ABR17 and WT seedlings without any stress, upon NaCl treatment the transgenic seedlings had significantly (p < 0.05) higher levels of chlorophyll (Table 1). Our results indicate that the NaCl treatment had less deleterious effects on the ABR17-transgenic seedlings compared to the WT.
In order to further characterize the differences between the WT and ABR17-transgenic lines, the ability of both WT and ABR17-transgenic seedlings to germinate in the presence or absence of light at RT was compared. In the dark, 85% of ABR17-transgenic A. thaliana had germinated after one week, whereas only 10% of the WT seeds had germinated under the same conditions (Table 1, Figure 2B). In contrast, in the presence of light, 100% of both ABR17transgenic and WT seeds had germinated in the same period (data not shown). Most Arabidopsis ecotypes require light for germination, which is primarily controlled by a reversible red light dependent equilibrium of the photoreceptors [23]. It is also known that exogenous CKs can substitute for red light and enhance the germination of certain light-requiring species in the dark [24][25][26][27]. Furthermore, A. thaliana detiolated (det) mutants exhibits many characteristics of seeds germinated in the presence of light even when germinated in dark [28], a phenotype that has been attributed to CKs because of the fact that even WT seedlings exhibit the same phenotype when germinated in the dark following exogenous CK application [29]. A role for CKs can also be inferred from the observation that coumarin or far-red light, both of which prevent the formation of CK-nucleosides from storage forms of CKs, inhibit germination of lettuce seeds in the dark [30]. Interestingly, amp1 A. thaliana mutants, which possess higher endogenous CKs, also exhibited a photomorphogenic response that in part similar to our ABR17-transgenic seedlings [31]. Taken together, all these results seem to suggest that endogenous CKs play an important role in the germination of light-sensitive seeds and the elevated endogenous CKs in ABR17-transgenic seedlings previously reported [8] may be responsible for the enhanced germination of this genotype in the dark (Table 1, Figure  2B).

Transcriptional profiling using microarrays
In order to characterize the molecular changes brought about by the expression of pea ABR17 cDNA in A. thaliana that resulted in the observed phenotypes, we analyzed gene expression by profiling the transcripts of ABR17transgenic plants in the absence and presence of 100 mM NaCl. As described earlier, the first set of microarray analysis was the investigation of the differences in gene expression between ABR17-transgenic and WT A. thaliana in the absence of NaCl (ABR17/WT). The second set of microarray analysis was between 100 mM NaCl-treated WT and untreated WT A. thaliana (100 mM NaCl treated WT/WT). The third set of microarray analysis was between 100 mM NaCl treated ABR17-transgenic versus untreated ABR17transgenic A. thaliana (100 mM NaCl treated ABR17/ ABR17).
Microarrays (70-mer oligonucleotide microarrays) consisting of probes presenting 23,686 unique genes identified by Arabidopsis genome initiative (AGI) locus identifiers were used. We identified transcripts as those with mean signal intensities that differed significantly from 0 at α = 0.05 in a Student's t-test in each set of microarrays. The transcripts were categorized based on shared structural elements and/or inferred function. We selected 12 genes representing different functional categories, which according to our microarray analysis showed enhanced or reduced levels of transcript abundance to validate our microarrays. The results from microarrays and qRT-PCR analysis are discussed below.

First set of transcriptional profiling: genes responsive to ABR17
Of the significantly responsive transcripts due the expression of pea ABR17 in A. thaliana, 124 were observed to be modulated in the transgenic line at least 1.5-fold compared to WT with 83 increasing and 41 decreasing in transcript abundance (Additional file 1). Many of these genes had annotations that were associated with either defense or plant growth and development, or both. A total of 16 genes showed significant differences in transcript abun-Appearance of WT and ABR17 transgenic A. thaliana in response to treatments dance about 2-fold, where 13 genes exhibited increased transcript abundance and 3 genes showed a decrease in transcript abundance (Table 2).
Among the highly induced transcripts in transgenic seedlings that were putatively related to defense responses (Table 2), we detected 5 members of the plant defensin (PDF) family which exhibited an increased in abundancẽ 2-3-fold in the transgenic line. PDFs are small (45-54 amino acids), highly basic cysteine-rich peptides belonging to the large defensin family, and are present throughout the plant kingdom. These proteins are known for their involvement in ancestral non-specific innate immune defense system [32]. In addition to being involved in mediating plant responses to pathogens, defensins may also play an important role in plant growth and development. For example, the constitutive expression of AtPep1 induced the expression of PDF1.2 which resulted in better root development in A. thaliana suggesting that plant defensins may regulate root development [32].
Another interesting transcript that exhibited increased abundance (2-fold; Table 2) in ABR17-transgenic plants was a putative mitogen-activated protein kinase (MAPK). MAPK cascades are known to play crucial roles in physiological processes such as cell growth, cell cycle regulation and developmental control as well as plant defense signaling [33]. They are also known to activate WRKY type transcription factors that are involved in transcriptional activation of disease resistance genes [34]. Indeed, we have observed a modest, but elevated expression of four genes belonging to the WRKY family and disease resistance protein (Additional file 1).
We also observed increased transcript abundance for several genes involved in plant growth and development ( Table 2). For example, expansins were detected as highly induced transcripts in ABR17-transgenic A. thaliana ( Table  2). Expansins are cell wall proteins that are known to induce pH-dependent plant cell wall extension and stress relaxation [35]. The expansins have been related to cell differentiation in tissues such as xylem, leaf primordia and root hairs [36][37][38]. Previous studies on transgenic plants expressing expansin genes have demonstrated precocious leaf development, longer petioles and larger leaf blades [39,40].
Glycine-rich proteins (GRPs) were also detected among growth related genes that whose transcripts increased in abundance in ABR17-transgenic plants ( Table 2). GRPs consist of quasi-repetitive glycine-rich domains, most commonly GGGX, GGXXXGG or GXGX repeats [41]. Some GRPs have been reported as structural components of the plant cell walls based on their co-localization with cell wall [42]. GRPs have also been reported to be activated by osmotic stress [43], cold shock [44] and wounding [45].
The genes that exhibited significant enhanced expression in ABR17-transgenic plants also included genes for proline-rich protein (PRP) family, xyloglucon endotransglycosylase (XTH), glycosyl hydrolase (GH), phytosulfokine precursor 2 (PSK2), No Apical Meristem (NAM) protein family and glutaredoxins (Additional file 1). PRPs represent a family of structural cell wall proteins that have been implicated in various plant developmental processes [46,47]. Similarly, XTH and GH family genes are involved in structuring xyloglucan cross-links in plant cell wall and  [37,48] and [49]]. The PSK2 gene is also involved in cell growth and differentiation [50][51][52]. Similarly, the NAM gene product is required for shoot apical meristem (SAM) formation during embryogenesis as well as for normal flower development [53][54][55]. Glutaredoxins have also been demonstrated to be involved in flower development, probably by mediating post-translational modifications of target proteins required for normal petal organ initiation and morphogenesis [56]. Our current observations that the significantly (albeit modest) higher expression of the above mentioned genes related to growth and development including flowering correlates well with the observed phenotypes which include early flowering, increased lateral branching and seed pods as observed in ABR17-transgenic A. thaliana ( Figure 1C).

A role for cytokinins in ABR17-induced changes in gene expression?
Interestingly, members of most of the gene families described above (defensin, expansin, MAPK, NAM, WRKY, GRP, PSK2 and Glutaredoxin) that are involved in plant defense as well as growth and development, have been previously reported to be regulated by CKs. For example, genome-wide expression profiling of immediate-early and delayed CK-response genes of A. thaliana, has led to the identification of many genes that are upregulated by CKs including members of expansin (At1g69530), GRPs (At2g21060), NAM (At4g27410), Fbox protein (At3g61060), ERBFs, putative ring zinc finger protein (At1g76410), a member of the bHLH family (At2g18300), blue copper binding protein (At5g20230) and PSK2 [57]. The blue copper binding protein (At5g20230) and putative ring zinc finger protein (At1g76410) identified by Brenner et al. (2005) [57] as CK-induced were observed to be up-regulated in our microarrays analysis. Similarly, gene expression analysis of transgenic A. thaliana seedlings transformed with bacterial isopentenyl transferase (IPT) [58] gene revealed increased transcript abundance for many members of MAPK and WRKY gene families, which included the specific WRKY gene -At1g80840 that has been detected in our microarray studies as being induced by ABR17 expression (Additional file 1). Another investigation into CK action in Arabidopsis has demonstrated increased expression of genes for cytochrome P450, PDF, expansin, patatin, WRKY members and putative disease resistance protein in response to CKs [59]. Therefore, it is apparent that several genes whose transcript levels were modulated by ABR17 expression in A. thaliana have been previously reported in the literature as being CK-responsive, there by suggesting an important role for CK-mediated gene expression in ABR17 action in planta.

Second set of transcriptional profiling: genes responsive to salt stress in WT A. thaliana
Microarray-based analyses of the salt responses in Arabidopsis have been published in several reports. However, most of these studies have investigated responses to very short-term exposure to salt. In this study, we report the transcriptional changes in A. thaliana as a result of longterm, continuous exposure to 100 mM NaCl. Here, we allowed A. thaliana seeds to germinate and grow on semisolid medium in the presence of 100 mM NaCl for 2 weeks, and the RNA extracted from whole seedlings were used for cDNA synthesis and subsequent microarray analysis. The results from microarray analysis of salt treated WT Arabidopsis seedlings (Additional file 2) agreed with previous studies using similar approaches [60,61]. We identified 163 genes that showed more than four fold changes in transcript abundance which have been previously reported as being responsive to salt. Our results, therefore, indicate that both short-term "shock" treatments with NaCl as well as long-term treatment used in this study elicit similar responses in A. thaliana at the transcript level (Additional file 2).
Members of protease inhibitor/lipid transfer protein (LTP) family were seen among highly up regulated and/or down regulated genes. At least five members showed increase in transcript abundance and 1 member showed decrease in transcript abundance of more than 4 fold (Additional file 2). LTP genes contain ABA-responsive (ABRE) element (GTACGTGG) and are induced by abscisic acid (ABA) [62,63]. It has been reported in the literature that NaCl, mannitol or ABA treatments induce the expression of a gene encoding an LTP-like protein in tomato [62,64]. In addition, the changes in the expression of LTP genes during salt stress have been previously reported [60,61]. Although, most of the LTP genes were up regulated after short term treatment with salt, they were found to be down regulated after 24 h of salt treatment [60,61]. From our studies, it appears that many of the LTP genes will get up regulated in response to long term stress as a result of expected increase in ABA levels.
Other major groups of genes with increase in transcript abundance following NaCl treatment included two members of glycosyltransferases (GTs) and five members of glycoside hydrolases (GHs). GTs and GHs are major families of carbohydrate-active that play a primary role in structuring xyloglucan cross-links in plant cell wall [48,49]. They have been previously reported to be induced by salinity stress in plants and this has been implicated in drought and salt tolerance in A. thaliana [49,60]. Other genes exhibiting increased transcript abundance included ribonuclease RNS1, osmotin-like protein, hydroxycinnamoyl benzoyltransferase-related, oxidoreductase, 2OG-Fe(II) oxygenase family, glutathione trans-ferase and zinc finger (C3HC4-type RING finger) protein family (Additional file 2). Similarly, the genes which showed decrease in transcript abundance more than 4 fold included many photosynthetic genes, plant defensins, heat shock proteins, auxin-induced proteins, disease resistance protein, Bet v I allergen family and bHLH protein. These results are once again consistent with the previously reported results from microarraybased investigation into salinity stress responses [60,61].

Third set of transcriptional profiling: genes responsive to salt stress in presence of ABR17
The results from microarray analysis of salt treated ABR17 transgenic A. thaliana seedlings are presented in Tables 3  and 4. We identified 129 genes showing either increase or decrease in transcript abundance more than 4-fold, which included transcription factors (15), stress responsive genes (16), carbohydrate and cell wall metabolism (8), electron transport and oxidoreductases (6), lipid metabolism (3), protein and amino acid metabolism (9), proteins involved in transport across membranes (12) and 60 unknown or unclassified genes.
Transcriptional factors are necessary for the proper transcriptional regulation in response environmental cues [65] and those exhibiting significant increases in transcript abundance included bHLH, 4 members of AP2 related, 2 members of NAM, zinc finger (C3HC4-type RING finger) protein family, ATMYB74 (MYB domain protein 74), ATHB-7 (Arabidopsis thaliana HOMEOBOX 7), and WRKY families. bHLH092 has been indicated among the highly induced transcripts in response to NaCl treatment in the previous transcriptomic studies and are suggested to be important regulators of the NaCl-stress response in Arabidopsis [60]. The APETALA2 (AP2) domain defines a large family of transcription factors which play important roles in plant growth, development as well as stress tolerance [66][67][68][69][70][71].
Similarly, as previously stated, NAM genes have been found to be induced by abiotic stresses implying roles in stress responses in addition to those in plant growth and development [72,73]. NAM/NAC proteins belong contain highly conserved NAC (for NAM, ATAF1, 2, and CUC2) domains in their N-terminal regions, that specifically binds target DNA [54]. It has also been demonstrated that NAC transcription factors are ABA-responsive [74,75] and are also induced by other plant hormones like NAA and ethylene [75][76][77][78]. Overexpression of NAC genes has been shown to result in an increase in lateral roots, and tolerance to abiotic stresses like drought and salt stress. NAC genes are believed to exert their stress ameliorating activity through the regulation of stress-inducible genes [77][78][79]. Similarly, the WRKY family TF genes and myb family genes are known to be biotic and/or abiotic stress respon-sive [60,80]. Thus, it is possible that the increased tolerance of ABR17-transgenic seedlings to NaCl is the combined effect of the modulation of the levels of abundance of transcripts for these transcription factors with demonstrated roles in stress tolerance.
The highest transcript abundance of any gene observed in salt treated ABR17 plants was XTR-6, which showed 4.7fold increase, compared to the untreated ABR17-transgenic line. Xyloglucan endo-transglycosylase (XET) has been suggested to be a key enzyme involved in the modification of the xyloglucan cross-links that controls the strength and extensibility of the plant cell wall [81]. Three members of GH family were also seen among genes which up regulated more than 4-fold change ( Table 3). The importance of GHs genes in plant stress [49,60] has already been discussed in the previous section.
Others salt responsive genes in the ABR17-transgenic line included osmotin, mannitol dehydrogenase, steroid sulfotransferases and RD20 (Table 3) which are known to be regulated by ABA, are expressed in salt-stressed plants and have been used to engineer salinity tolerance [82][83][84][85][86]. In addition, we also observed increase in transcript abundance for ribonuclease-RNS1, peroxidases, copper/zinc superoxidase dismutase (CSD1), cytochrome p450 family, MATE efflux protein and protein kinases which have been previously demonstrated to accumulate in salt treated tissues by others [60,61]. From our microarray results, it appears that many genes involved in mediating responses to salinity stress are increased in transcript abundance as would be expected.

Comparison of salt responses in WT and ABR17 transgenic seedlings
Although transcriptional changes were almost similar both in salt treated-ABR17 and WT seedlings, the transcript abundance of some genes exhibited significant differences in both the trend as well as the degree of modulation of transcript abundance (Table 5). For instance, as mentioned previously, transcript abundance of xyloglucan endotransglycosylase (XTR-6) (At4g25810) increased 4.7-fold in salt treated ABR17 seedlings, whereas it showed only a 2.4-fold increase in salt-treated WT seedlings (Table 5). Similarly, AP2 domain related transcription factor RAP2.6 (At1g43160) increased 4.5fold in salt treated ABR17 compared to 1.67-fold in treated WT plants. The expressed proteins-ABA-responsive protein-related (At3g02480) and unknown protein (At5g24640) also showed increase in transcript abundance of at least 4-fold in salt-treated ABR17 transgenic line compared to the 2-fold increase observed for the WT in response to salt. Other genes which exhibited increase in transcript abundance of more than 2 fold in salt-treated ABR17 transgenic but showed less abundance in treated WT included pathogenesis-related protein, (At4g33720) and glutamine-dependent asparagine synthetase (At3g47340). On the other hand, the retroelement pol polyprotein (At2g04460) with unknown function showed a down regulation of more than 2-fold in salt-treated transgenic ABR17 line, but was up regulated in salt treated WT Arabidopsis plants.
Interestingly, many members of heat shock protein (Hsp) family and PDF family showed the opposite response in salt-treated ABR17-transgenic seedlings compared to the salt-treated WT counterparts, with an increase of transcript abundance in salt-treated ABR17-transgenic (Table 5). This difference in the direction of the response in gene expression i.e. induction in the transgenic seedlings versus the repression in the WT may have important consequences with respect to the ability to tolerate salinity (and perhaps other) stress. For example, the Hsp family contains chaperones, which have important roles in protein folding, assembly and in the disposal of unwanted nonfunctional proteins. Hsps are usually induced by environmental stress, and the accumulation of Hsps coincides with enhanced stress tolerance [60,[87][88][89][90][91]. In addition, transgenic Arabidopsis plants overexpressing AtHSP17.7 accumulate high levels of AtHSP17.7 protein and show enhanced tolerance to drought and salinity [89,90]. The abundance of Hsps in plants and their functional characteristics suggest that Hsps play an important role in plant stress tolerance. Thus, the increased abundance of HSP transcripts in the ABR17-transgenic seedlings may be important for the increased tolerance of this line to the imposed stress.
The up regulation of PDFs in ABR17-transgenic A. thaliana grown under normal conditions ( Table 2) and their importance in growth and development has already been discussed earlier. The literature also supports a role for PDFs in stress tolerance [32,91]. Most of the previously characterized PDFs exhibit anti-fungal, antibacterial, antiinsect and protease inhibitor activity [92]. However, the halophyte salt cress (Thellungiella halophila), a relative of Arabidopsis over expresses PDFs under normal conditions and hence defensins are believed to play a role in salt tolerance [93]. It is therefore possible that the observed relatively tolerant phenotype of ABR17-transgenic seedlings could be due, at least in part, to the elevated expression of XTR6, RAP2.6 transcription factors, unknown proteins (At3g02480, At5g24640), Hsp and PDF gene(s). Additional studies are underway in our laboratory to precisely Spermine synthase (ACL5) -2.02 0.14 2.72E-05 All expression ratios are significant (α = 0.05) and are in a log2 scale where fold change is salt treated ABR17/control ABR17. AGI a -Arabidopsis Genome Initiative SE b -Standard error  characterize the role(s) of the above-mentioned genes in stress tolerance of ABR17-transgenic A. thaliana seedlings.

qRT-PCR validation of microarray observations
In order to confirm the fact that CK-responsive genes were indeed up-regulated in the ABR17-transgenic lines, we performed qRT-PCR experiments with the following genes: plant defensin protein (PDF1.2a, At5g44420), expansin (EXPL1, At3g45970), GRP (At1g07135) and putative MAPK 11 (Atg01560) using qRT-PCR. Among the CK-inducible genes identified from our first set of microarray experiments we chose the above 4 genes for qRT-PCR as their transcripts were observed to be at least ~2fold (1 in log2 ratio) or higher in the transgenic line compared to WT. Our microarray analysis revealed increase in transcripts for defensin, expansin, GRP and MAPK in pea ABR17 seedlings by 2.6, 2.5, 1.9 and 1.9-fold, respectively ( Figure 3). Our qRT-PCR results were consistent with the microarray data and showed up-regulation of defensin, expansin, GRP and MAPK by 3.6, 2.5, 2.5 and 2.7-folds, respectively ( Figure 3). From these results it is apparent that all the four genes that were up-regulated in our microarray analysis also demonstrated up-regulation in the qRT-PCR relative expression experiments (Figure 3).
Because of the fact that the specific members of gene families whose transcripts were detected to be modulated by ABR17 cDNA expression in A. thaliana were not exactly identical to those specific members of these families identified by other studies, we wanted to investigate whether those specific members detected in our studies were indeed CK-inducible/repressible. In these experiments, we used WT A. thaliana tissue germinated and grown for 14-days medium supplemented with 5 μM zeatin for additional qRT-PCR experiments. This concentration of zeatin was chosen based on our earlier observations that it induced the largest phenotypic responses in A. thaliana when exogenously applied [8]. It must also be noted that Footnote: All expression ratios are significant (α = 0.05) and in a log2 scale where fold change is salt treated ABR17/control ABR17 and salt treated WT/control WT. ABR17-WT = Difference in log2 ratio of salt treated ABR17/control ABR17 and salt treated WT/control WT. AGI a -Arabidopsis Genome Initiative SE b -Standard error even though 5 μM zeatin was used in our experiments, it is difficult to estimate how much of this exogenously supplied CK actually gets into the seed in order to exert a physiological effect. From the results shown in Figure 3, it is apparent that the expression of EXPL1 (At3g45970), putative MAPK (Atg01560) and GRP (At1g07135) was upregulated in response to exogenous zeatin by 1.6, 3.8, and 2-fold, respectively ( Figure 3). In contrast, the expression of defensin gene was observed to be down-regulated (0.3fold) in response to the exogenous application of CK (Figure 3). The results for expansin, MAPK and GRP are consistent with our microarray and qRT-PCR results with respect to increased transcript abundance in ABR17-transgenic A. thaliana previously shown to possess higher endogenous CK concentrations [8]. However, in the case of defensin, even though our microarray and qRT-PCR experiments revealed that this gene was up-regulated in the ABR17-transgenic line (Figure 3), its expression was not induced by exogenous CK (Figure 3). The reason for this discrepancy is not immediately clear; however, this may be due to the concentration as well as type of CK used for our exogenous experiments. Furthermore, as indicated previously, the amount of the exogenously supplied CK entering the seed to exert physiological affect may also be different from the concentrations required to elicit induction of this gene.
In order to confirm results from our second and third set of microarray analysis, we performed qRT-PCR experiments with the following 12 genes: unknown proteins (At3g02480; At5g24640; At1g14880), XTR6 (At4g25810), bHLH (At5g43650), AP2 domain transcription factor RAP 2.6 (At1g43160), ATNAC3 (At3g15500), ACD6 (At4g14400), PDF1.2a (At5g44420), EXPL1 (At3g45970), GRP (At1g07135) and MAPK 11 (Atg01560). The unknown proteins were chosen because expression of two of them (At3g02480 and At5g24640) were among highly induced transcripts in salt treated ABR17-trangenic line and also showed comparatively less but high level of transcript abundance in salt treated WT A. thaliana lines ( Table  5, Figure 4). Two of them (At1g14880 and At4g14400) were among highly down regulated genes in both salt treated ABR17-transgenic line and WT A. thaliana lines (Table 3 and 4, Figure 4). Our qRT-PCR data showed similar trend as observed by microarrays for all the abovementioned genes in both salt-treated ABR17-transgenic and salt-treated WT microarrays (Figure 4). The gene XTR6 (At4g25810) was selected because it was among one of the most highly induced transcripts of any gene on our salt treated ABR17-transgenic A. thalina microarray ( Table 3, Figure 4). The genes bHLH (At5g43650), RAP 2.6 (At1g43160) and ATNAC3 (At3g15500) were chosen because their expression was the highest among any other transcription factors identified in response to salt in ABR17-transgenic line ( Table 3

Relative expression of CK-biosynthetic genes (IPT and CKX) in ABR17-transgenic A. thaliana
As discussed earlier, our observations indicated that many of the genes identified in the transgenic plants as being up-regulated are from families that contain CK-responsive members. We have also previously reported higher endogenous concentrations of CK in this line [8], which suggested the possibility that this may be due to either enhanced de novo CK biosynthesis or decreased degradation. Specifically, the endogenous concentration of total CK in the transgenic line used in this study was ~1-3-fold higher, with the concentration of zeatin (cis and trans combined) being ~1.4-fold and IP being ~2-fold higher in these transgenic lines. However, we did not detect any IPT (involved in CK biosynthesis) or CKX (involved in CK degradation) genes as being significantly up-or down-regulated genes in our microarray experiments suggesting that the elevated endogenous CK concentrations previously reported may not be the result of increased or decreased activities of IPT and CKX genes, respectively. In order to confirm our microarray results and to lend additional support to our above-mentioned hypothesis with respect to the roles (or lack thereof) of IPT and CKX expression in ABR17-transgenic A. thaliana, we also performed qRT-PCR analysis of the expression of IPT and CKX genes using qRT-PCR. There are 9 known IPT genes and 7 known CKX genes but sequence of CKX5 and 7 are very similar therefore we performed qRT-PCR analysis on the 9 IPT and 6 CKX genes. The results from these experiments are summarized in Table 6 and it is apparent that most of the IPT genes exhibit similar expression patterns in both transgenic and WT seedlings. The only exception appears to be IPT 8 where only 0.5-fold expression of this gene was observed in the transgenic line (Table 6). Similarly, CKX expression in the transgenic line was also quite similar to its expression in the WT (Table 6). Our results suggest that the differences in endogenous CK concentrations previously observed in the ABR17-transgenic line may not be the result of increased IPT or decreased CKX levels. However, frequently, there is no correlation between transcript abundance and protein levels and therefore it is possible that IPT and/or CKX protein concentrations may have been affected in the transgenic line resulting in increased endogenous CKs as a result of posttranslational processes. However, our previously reported proteome studies on this transgenic line did not reveal any differences between transgenic and WT seedlings with respect to the levels of these proteins [19].
It is possible that the activity of neither IPT nor CKX is responsible for the increased endogenous concentrations of CKs in the ABR17-transgenic lines and the increased endogenous CKs previously reported in the ABR17-trans-genic lines may be the result of tRNA degradation by the previously demonstrated RNase activity of pea ABR17 protein [8]. Thus, an increase in free cellular CK would not necessarily involve enhanced IPT or reduced CKX activity; rather it may reflect an increased access to existing, yet tRNA-bound, CK.

Conclusion
We have demonstrated that pea ABR17 cDNA expression modulates the level of a number of transcripts related to plant defense, growth and development, which may explain the observed phenotypic differences between WT and ABR17-transgenic A. thaliana. The gene expression of many transcription factors and defense responsive genes like Hsps and PDFs showed different degree and kind of response between salt treated-ABR17 transgenic and WT A. thaliana, which explains the observed enhanced germination and early seedling vigor in ABR17 transgenic lines, compared to its WT counterpart. Many of the genes exhibiting a 2-fold or higher increase in transcript abundance are known CK-responsive genes providing additional evidence of a role for CKs in ABR17 function. Furthermore, a detailed expression analysis of IPTs and CKXs revealed that the levels of these transcripts were similar in both WT and transgenic seedlings, suggesting the possibility that ABR17 modulates endogenous CKs through an, as of yet, uncharacterized mechanism including the possible degradation of tRNAs which contain CK moieties [94]. These possibilities are currently being investigated in our laboratory.

Plant material and growth conditions
Transformation of A. thaliana with the pea ABR17 cDNA and the generation of homozygous ABR17-transgenic A. thaliana (line 6.9) have been previously described [19]. This line (6.9) was one of the three independently derived transgenic lines that were characterized in that earlier study. The WT (ecotype WS) and transgenic A. thaliana plants were grown in the green house for observations as previously described [8]. Lateral branches were counted on plants from three independent biological replicates with at least 72 plants per replicate. Average number of days required for the opening of the first flower was also recorded on plants from three biological replicates with 36 plants in each replicate.
In order to measure root lengths of seedlings seeds of A. thaliana (line 6.9) and the WT were surface sterilized [8] and placed on half strength Murashige & Skoog (MS) medium [22] with or without salt (75 mM NaCl or 100 mM NaCl) in square dishes with grids. These dishes were placed vertically in a growth chamber (at 21°C and with light intensity of 250 μmol m -2 s -1 ) and root lengths were measured after 10 days. The seeds of ABR17 and WT seeds were also grown on half strength MS medium with 0 or 100 mM NaCl to determine their fresh weight and chlorophyll and carotenoid contents in order to assess their ability to grow in the presence of salt. The length of the primary roots of 10-day-old seedlings from three independent biological replicates with at least sixty seedlings per replicate, were calculated using the Image J software (Image J, NIH, MD, USA).
Chlorophyll and carotenoids were extracted from the pooled 2-week-old tissue grown on MS media, using the procedure as described by Srivastava et al., 2006 [7]. Total chlorophyll was estimated using a nomogram [95] and total carotenoid was measured using the formula: where A is the absorbance and CAR is the carotenoid content [96]. The fresh weight, chlorophyll and carotenoid were calculated using pooled tissue from three independent biological replicates. Percent germination after one week for ABR17-transgenic and WT seeds in the dark and in the presence of light (fluorescent light, 30 μmol m -2 s -1 ) were compared in Petri dishes at RT. This experiment included three independent biological replicates with at least 45 seeds per replicate. All statistical analyses were performed using the Student's t-test procedure in SAS version 8e (Statistical Analysis System 1985).
Tissue for microarray analysis was obtained by placing surface sterilized seeds of A. thaliana (line 6.9) and the WT on half strength MS medium in Petri dishes with or without 100 mM NaCl at RT (21 ± 2°C) under continuous fluorescent light 30 μmol m -2 s -1 for 14 days. Seedlings (14day-old) from three independently grown biological replicates in all three set of experiments (comparison of ABR17-transgenic with WT without any stress; comparison of salt treated WT with untreated WT; comparison of salt treated ABR17 transgenic with untreated transgenic) were removed from the MS plates, flash frozen in liquid nitrogen and stored at -80°C until used for RNA extraction.

RNA extraction, cDNA synthesis and microarray analysis
In order to investigate the ABR17-induced gene expression changes under normal and salinity stress conditions, we conducted microarray analysis in three separate hybridization experiments. The first set (set I), consisted of comparison of cDNA samples prepared from ABR17transgenic and WT seedlings, which were grown in the absence of any stress. Set II consisted of cDNA obtained from salt-treated samples of WT and untreated WT; and set III, cDNA samples of salt treated ABR17-and untreated ABR17-transgenic seedlings for hybridization to the oligonucleotide arrays. Each microarray experiment consisted of six hybridizations according to the principles of dyeswap design [97] on tissues across three biological replicates of the experiments.
RNA was isolated using the QIAGEN RNeasy Plant Mini Kit (Qiagen Inc., Mississauga, ON, Canada) from 2-weekold WT and ABR17 seedling tissue from all three set of experiments and the integrity of all RNA samples assessed by agarose gel (1.2%) electrophoresis. For microarray hybridization, 6 μg of total RNA was used to synthesize cDNAs using SuperScript ® II RT (Invitrogen Inc., Burlington, ON, Canada) with RT polyA-capture primers in 3D Array 900TM (Genisphere Inc., Hatfield, PA, USA). In these microarray experiments, 70-mer oligonucleotide arrays were used which contained 26,090 probes (Array-Ready Oligo Set for Arabidopsis genome Version 1.0, Qiagen Operon, Alameda, CA, USA), plus additional probes for quality control. Oligonucleotide arrays were spotted on superamine aminosilane-coated slides (TeleChem International Inc., Sunnyvale, CA, USA). Each pair of samples within each of the three biological replicates was labeled in a reciprocal dye-swap design, for a total of 18 hybridizations (overnight, at 52°C) in all three sets of experiments. Slides were scanned using ArrayWoRx e (Applied Precision, Issaquah, WA, USA) and spot intensities were measured, quantified, normalized and analyzed using TM4 [98]. Spots with intensity ratios that differed significantly from 0 (log2 scale) were identified by Student's t-test. This procedure highlights the spots that demonstrated statistically significant differential expression between the different samples. The raw microarray data of 18 hybridizations as well as the protocols used to produce the data were deposited in the ArrayExpress database [ArrayExpress: E-MEXP1024 and E-MEXP1566].

Quantitative real-time PCR (qRT-PCR) validation of microarray data
Primers for qRT-PCR were designed using the Primer Express software (Applied Biosystems Inc., Foster City, CA, USA) to ensure that PCR products of approximately 70-80 bp were generated (Additional file 3). cDNA synthesis and qRT-PCR analysis of gene expression of 19 genes were performed using the Taqman system as described previously [8] on an ABI Prism 7700 Sequence detector (Applied Biosystems Inc., Foster City, CA, USA) and the SNP RT template program or using the SYBR green system as described by Yang and others [99] was used to validate the expression of 8 genes. In both cases, the deltadelta method [100] was used to calculate relative gene expression using actin as the endogenous control. The relative transcript abundance in the controls was normalized to 1 and was used as a basis for comparison to the treatments. Plant tissue from three biological replicates was used in qRT-PCR experiments and reactions for each biological replicate were performed in duplicate (n = 6).
Publish with Bio Med Central and every scientist can read your work free of charge