- Research article
- Open Access
Comparative transcriptome analysis of Arabidopsis thaliana infested by diamond back moth (Plutella xylostella) larvae reveals signatures of stress response, secondary metabolism, and signalling
© Ehlting et al; licensee BioMed Central Ltd. 2008
- Received: 17 November 2007
- Accepted: 09 April 2008
- Published: 09 April 2008
Plants are exposed to attack from a large variety of herbivores. Feeding insects can induce substantial changes of the host plant transcriptome. Arabidopsis thaliana has been established as a relevant system for the discovery of genes associated with response to herbivory, including genes for specialized (i.e. secondary) metabolism as well as genes involved in plant-insect defence signalling.
Using a 70-mer oligonulceotide microarray covering 26,090 gene-specific elements, we monitored changes of the Arabidopsis leaf transcriptome in response to feeding by diamond back moth (DBM; Plutella xylostella) larvae. Analysis of samples from a time course of one hour to 24 hours following onset of DBM feeding revealed almost three thousand (2,881) array elements (including 2,671 genes with AGI annotations) that were differentially expressed (>2-fold; p[t-test] < 0.05) of which 1,686 also changed more than twofold in expression between at least two time points of the time course with p(ANOVA) < 0.05. While the majority of these transcripts were up-regulated within 8 h upon onset of insect feeding relative to untreated controls, cluster analysis identified several distinct temporal patterns of transcriptome changes. Many of the DBM-induced genes fall into ontology groups annotated as stress response, secondary metabolism and signalling. Among DBM-induced genes associated with plant signal molecules or phytohormones, genes associated with octadecanoid signalling were clearly overrepresented. We identified a substantial number of differentially expressed genes associated with signal transduction in response to DBM feeding, and we compared there expression profiles with those of previously reported transcriptome responses induced by other insect herbivores, specifically Pieris rapae, Frankliniella occidentalis, Bemisia tabaci,Myzus persicae, and Brevicoryne brassicae.
Arabidopsis responds to feeding DBM larvae with a drastic reprogramming of the transcriptome, which has considerable overlap with the response induced by other insect herbivores. Based on a meta-analysis of microarray data we identified groups of transcription factors that are either affected by multiple forms of biotic or abiotic stress including DBM feeding or, alternatively, were responsive to DBM herbivory but not to most other forms of stress.
- Differentially Express
- Plutella Xylostella
- Signal Transduction Component
- Methylerythritol Phosphate
- Auxin Related Gene
Arabidopsis thaliana has emerged as a useful system for genomic studies of plant insect-interactions [1–6]. Because of the large amount of genomic information available for Arabidopsis, it is possible to perform comparisons of gene expression profiles across many different conditions or treatments including various forms of interactions with herbivores and pathogens [1, 2]. With regard to specific pathways involved in plant defence against insects, the Arabidopsis genomic resources have much advanced, for example, the discovery of genes and proteins of secondary metabolism (specifically glucosinolate, phenolic, and terpenoid metabolism) [7–11], as well as genes involved in plant-insect defence signalling [12–14].
Previous large-scale gene expression microarray analyses of Arabidopsis-herbivore interactions involved plants affected by Pieris rapae (larvae of cabbage white butterfly), Spodoptera littoralis (larvae of mediterranean brocade), Frankliniella occidentalis (western flower thrip), Bemisia tabaci (silverleaf whitefly nymphs), Brevicoryne brassicae (cabbage aphid), and Myzus persicae (green peach aphid) [1–5, 15]. These insects represent leave-chewing larvae (P. rapae and S. littoralis) as well as cell-sucking (F. occidentalis) or phloem sap-feeding (M. persicae, B. brassicae, B. tabaci) adults with P. rapae and B. brassicae being specialist herbivores adapted to members of the Brassicaceae as their hosts. The present study complements previous work with an analysis of Arabidopsis rosette leaves fed upon by larvae of a different leave-chewing specialist herbivore, Plutella xylostella (diamond back moth – DBM). DBM larvae feed on several crucifer plants and are a frequent pest of agricultural crops including cabbage, broccoli, cauliflower, and rape .
Overall, our findings from a fully replicated time-course transcriptome analysis of Arabidopsis challenged by DBM larvae identified almost three thousand (2,881 array elements; 2,671 genes with AGI annotations) differentially expressed genes (>2-fold; p[t-test] < 0.05) and several distinct temporal patterns of changes of transcript abundance with prominent changes of transcripts associated with stress response, secondary metabolism, and signalling. In addition, we provide a first comprehensive meta-analysis of array data of herbivore-induced Arabidopsis transcription factors, which identified insect-induced transcription factors that are also affected by other forms of biotic or abiotic stress as well as transcription factors that appear to be more specific to the insect-induced response.
Overall changes of the Arabidopsis leaf transcriptome in response to DBM
Arabidopsis rosette leaves (ecotype Ler) were challenged with feeding DBM larvae (third to fifth instars). For microarray gene expression profiling, rosette leaves were harvested after 1 h, 4 h, 8 h, and 24 h of continuous DBM feeding. For each time point leaves were also harvested from unchallenged control plants. In order to obtain enough plant material for RNA isolation (leave material became limiting in particular at the later time points of DBM feeding) each treatment and time point consisted of four or five plants grown together in one pot and exposed to a group of seven DBM larvae. For each treatment and control and for each time point two independent biological replicate experiments were performed and RNA derived from each biological replicate was used for two separate array hybridizations each using reversed fluorescence labels (dye-flip). This experimental design thus resulted in four replicate microarray hybridizations per time point and treatment with two biological and two technical replicates comparing RNA derived from treated plants with the corresponding control harvested in parallel.
Overall summary of differentially expressed genes
p(t-test) < 0.01
p(t-test) < 0.05
p(t-test) < 0.05, FCb>2
In at least 1 time point
p(ANOVA) < 0.01
p(ANOVA) < 0.05
p(ANOVA) < 0.05, FCc>2
Temporal patterns of the Arabidopsis leaf transcriptome affected by DBM
Similarly, many down-regulated genes displayed transient expression profiles (cluster G and cluster H), although a majority of down-regulated genes maintained lower expression levels over the time course analyzed (cluster I). Interestingly, a portion of genes placed in the down-regulated clusters displayed reversed expression ratios at different time points, e. g. some genes in cluster H were transiently down-regulated at early time-points but were up-regulated later in the experiment. Likewise, many genes in cluster I were transiently up-regulated early in the treatment, but were down-regulated 24 h after the onset of DBM feeding. In reverse, many genes in cluster F were up-regulated 8 h into the time course, but displayed repressed expression at the 24 h time point. This cluster also contains genes that display a biphasic expression pattern, with repressed expression at early (1 h) and late (24 h) time points.
In summary, it is noteworthy that despite continuous insect feeding over the time period analyzed, a majority (60%) of up-regulated genes displayed a transient pattern of change of transcript abundance.
Annotation and expression profiles of DBM induced stress-related genes
Most of these DBM-affected stress-related genes were strongly up-regulated within 1 to 4 h after the onset of DBM feeding (Figure 3), and the majority were associated with wound response. In contrast, only five genes were up-regulated only at the 24 h time point, including two pathogenesis-related (PR) genes (PR1/At2g14610 and PRB1/At2g14580) that are associated with salicylic acid dependent pathogen defence. These genes are down-regulated early during DBM feeding before being up-regulated at 24 h. Many of the wound-response genes that are strongly up-regulated by DBM are involved in octadecanoid biosynthesis. All known enzymatic steps of the octadecanoid pathway were represented in the cluster of DBM-induced genes (Figure 3A). These include the lipoxygenase LOX2 (At3g45140), two other putative lipoxygenases (At1g17420 and At1g72520), the single copy allene oxid synthase AOS (CYP74A, At5g42650), the allene oxide cyclases AOC1 (At3g2576), AOC2 (At3g25780), and AOC4 (At1g13280), the 12-oxophytodienoate reductases OPR2 (At1g76690) and OPR3 (At2g06050), and 3-oxo-2-(2'-[Z]-pentenyl)cyclopentane-1-octanoic acid C oA L igase1 OPCL1 (At1g20510) [20, 21]. Other genes in this group of up-regulated wound-response genes are involved in the shikimate pathway (anthranilate synthase, ASB, At1g25220), and tocopherol biosynthesis (tyrosine aminotransferase, TAT, At2g24850) [22, 23].
We further analyzed expression patterns of genes of secondary metabolite pathways that are known to be affected by herbivory, namely glucosinolate, phenolic and terpenoid metabolism (Figure 3B–D). Many transcripts of these pathways were differentially expressed upon DBM feeding. For the glucosinolate pathway, two genes involved in the chain elongation of methionin (MAM1 and MAML, At5g23020 and At5g23010) were up-regulated, while two other members of the MAM gene family  were not detectibly differentially expressed. Also up-regulated were cytochrome P450 monooxygenases involved in the biosynthesis of various glucosinolates: CYP79A2 (At5g05260) catalyzing the conversion of phenylalanine to the corresponding oxime in benzylglucosinolate biosynthesis ; CYP79B2 (At4g39950) converting tryptophan and tryptophan analogs to the oxime in indole glucosinolate biosynthesis ; CYP79F2 (At1g16400) involved in the synthesis of long chain aliphatic glucosinolates ; and CYP83B1 (At4g31500) catalyzing the oxidation of indole-3-acetyldoxime in indole glucosinolate biosynthesis . CYP79B2 is also involved in camalexin and auxin biosynthesis [29, 30]. DBM feeding further induced the C-S-lyase (SUR1, At2g20610) that converts S-(alkylacetohydroximoyl)-L-cysteines to the corresponding thiohydroximic acids . While none of the three myrosinase encoding genes present in Arabidopsis  were differentially expressed at detectable levels, two myrosinase associated proteins (MAG) were affected. Of these MAG2 (At3g14210), which has been characterized as a quantitative trait locus (termed epithiospecifier 1; ESM1) controlling the ratio of nitrile to isothiocyanate production during glucosinolate hydrolysis , was suppressed by DBM feeding, while the related gene MAGL (At1g54030) was induced.
Relatively few genes that have been functionally characterized to encode enzymes of the phenylpropanoid pathway are induced by DBM feeding (Figure 3C), but many of the genes of this pathway just failed the call to be significantly differentially expressed at stringent conditions. In contrast, eight genes with similarity to known phenylpropanoid genes but otherwise of unknown function were transcriptionally up-regulated (Figure 3C).
Related to terpenoid biosynthesis, two differentially expressed genes encoding enzymes of the methylerythritol phosphate pathway, 1-deoxy-D-xylulose 5-phosphate synthase (DXPS1, At3g21500) and 2-C-methyl-D-erythritol 4-phosphate cytidyltransferase (MCT, At2g02500)  were down-regulated by DBM feeding. In contrast the only gene of the mevalonate pathway that was detected as affected by DBM feeding encodes acetoacetyl-CoA thiolase (AACT1, At5g47720) and was up-regulated. In addition, a farnesyl diphosphate synthase (FPPS1, At4g17190) and three of the more than thirty Arabidopsis terpene synthase genes (TPS04, At1g61120; TPS10, At2g24210; TPS21, At5g23960)  were up-regulated by DBM feeding (Figure 3D).
Expression profiles of phytohormone related genes affected by DBM feeding
Genes associated with gibberellic acid, brassinosteroids, and abscisic acid were also not over-represented among the differentially expressed genes (Figure 4). However, among those genes that were differentially expressed, most genes associated with abscisic acid were rapidly up-regulated upon herbivory, while genes associated with gibberellic acid were predominantly down-regulated (Figure 4 and Figure 5).
In addition to genes associated with jasmonate signalling, genes associated with the hormones auxin and cytokinin were also over-represented in the DBM affected transcriptome (Figure 4). While three cytokinin oxidase family members involved in cytokinin catabolism were up-regulated late in the feeding experiment, several cytokinin response regulators are transiently down-regulated during herbivory (Figure 5). Likewise, a large group of auxin-induced genes were transcriptionally down-regulated starting at 4 h after the onset of herbivory feeding (Figure 4). However, a smaller group of genes in the same group displayed a reverse expression pattern, as did some genes involved in auxin metabolism and the auxin-response transcription factor MONOPTEROUS (Figure 5) . Overall, our results highlight the importance of jasmonate in herbivore induced signalling, and may also suggest roles for cytokinin and auxin as well.
Expression profiles of genes associated with signal transduction
Of the 902 protein kinases present on the microarray, 98 were differentially expressed, 71 were up-regulated, while 27 were down-regulated (Additional File 2). Fourteen kinases were differentially expressed in at least two time points (Figure 6). Most genes in this group code for receptor-like kinases such as leucin rich repeat (LRR) and peptido-glucan (LysM) binding domain containing kinases. In addition, two mitogen activated protein kinase kinase kinases (MAPKKK) were transcriptionally up-regulated: ANP2 (At3g46160) is a MAPKKK protein related to Nicotiana protein kinase 1 (NPK1) which may negatively regulate stress responses ; and Raf27 (At4g18950) contains an ankyrin domain but has not been further characterized. A calcium-dependent protein kinase, At2g3036 (CIPK11, SnRK3.22) which is a member of a plant specific protein kinase family that specifically interacts with the calcium sensor protein CalcineurinB-like , was also up-regulated by DBM feeding. Finally, two AGC kinases (protein kinases A, G, and C), which belong to a family of effectors of the intracellular second messengers cAMP, cGMP, phospholipids, and Ca2+ , were up-regulated upon herbivory.
Targeted protein degradation via the ubiquitin/26S-proteasome pathway is another important regulatory process . Among the 1,403 Arabidopsis genes annotated to be involved in this pathway, 1,230 are present on the array and 82 were differentially expressed upon DBM feeding (Additional File 2). Among these are a 75 putative E3-ubiquitin-protein-ligases that were affected by DBM, in addition to five differentially expressed 26S-proteasome components, a single ubiquitin-like gene and a E2-ubiquitin activating enzyme (UBC10, At5g53300) (the latter two were down-regulated). Figure 6 shows the expression profile of ten different E3-ubiquitin-protein-ligases that were differentially expressed in at least two time points. In summary, we identified a large number of signal transduction components affected by DBM feeding. In particular, members of the AP2-EREBP family of transcription factors stand out as being rapidly induced by herbivory suggesting roles for this family in DBM induced signal transduction networks.
Infestation by feeding DBM larvae induces substantial overall changes in the Arabidopsis leaf transcriptome, with 1,854 array elements representing 1,664 annotated genes that were significantly induced and 1,007 elements representing 913 annotated genes that were repressed significantly. Despite continuous feeding the majority of differentially expressed genes displayed a transient expression profile with a maximum transcript abundance level at 8 h after onset of feeding and were down to their initial transcript level after 24 h (clusters B and C, Figure 2). It is not known what proportion of the induced changes of the transcriptome result in any downstream changes; and it is possible that some defence responses only require a transient burst at the level of the transcriptome to become effective and a large part of the initial response could return to the pre-attack steady-state level of gene expression. Repressor proteins, such as the recently discovered JAZ proteins involved in the mediation of jasmomate signalling [13, 14], may be involved in shaping such rapid and transient responses.
Although only a few plant species have been studied for the impact of insect attack on large-scale transcriptome changes, their diversity ranges from relatively short-lived herbaceous plants to long-lived trees, including angiosperms and gymnosperms. Results obtained with these systems support the general notion that insect feeding induces massive changes of the host plant transcriptome [1–6, 43, 44]. A few general patterns have emerged from these studies suggesting that herbivory can results in down-regulation of primary metabolic processes while at the same time activating defence related processes including secondary defence metabolism. These findings are well supported by our analysis of the Arabidopsis transcriptome affected by DBM feeding. In addition, the massive reprogramming of primary and secondary metabolic processes as part of the insect-induced defence response involves rapid changes in signalling and other regulatory processes. The present study establishes a signature of DBM-induced changes of the signalling transcriptome of Arabidopsis leaves.
Surprisingly, in the meta-analysis of all microarray data represented in this comparative Arabidopsis-herbivory transcriptome study (Figure 8 and Additional File 3), we found only one gene being up-regulated in nine of the ten experiments compared: The cytochrome P450 monooxygenase CYP79B2, which catalyzes the conversion of tryptophan to indole-3-acetaldoxime, the precursor of indole glucosinolates, camalexin, and also auxin [26, 29, 30, 45]. Six additional genes were found up-regulated in eight of the ten experiments, including three other genes related to tryptophan metabolism and glucosinolate biosynthesis as well as the jasmonate inducible tyrosine aminotransferase TAT, which is involved in tocopherol biosynthesis [22, 46]. An additional 40 genes were found induced in at least six of the experiments compared (Additional File 3), and almost half of these (18 genes) were found in the GeneOntology category 'response to stress' including genes encoding enzymes of the shikimate pathway (in particular the tryptophan branch), phenylpropanoid metabolism, glucosinolate bioynthesis, glutathione metabolism, and chlorophyll degradation (Additional File 3). Genes encoding functions in the octadecanoid pathway were also found induced in most of the individual datasets of this comparative Arabidopsis-herbivory transcriptome meta-analysis. The role of octadecanoids in mediating herbivore-induced responses is well established, and it has been estimated that up to 80% of all herbivore-induced Arabidopsis genes are octadecanoid regulated .
In the present study we identified a large number of transcripts that are affected by DBM feeding and are broadly associated with signal transduction components (examples are highlighted in Figure 6 and the complete data set is provided in Additional File 2). Transcriptome patterns associated with signalling in insect attack have not been well established. In other plant species this is largely due to the smaller array platforms commonly available that may not include many transcription factors or other signalling features. Also, lack of relevant reference datasets such as those obtained from array analyses of other stress treatments in the same plant species make meta-analyses of signalling transcriptomes currently a difficult task for most plant species, except for Arabidopsis. We established a first meta-analysis of DBM-induced regulatory proteins (protein kinases, transcription factors, ubiquitin ligases) to integrate results obtained in the present microarray analysis of DBM-induced Arabidopsis leaves with previously published Arabidopsis microarray data from various treatments related to biotic and abiotic stress response (Figure 9, Additional File 4). Specifically, we selected microarray data from a total of 295 Arabidopsis samples treated with a wide range of biotic and abiotic stresses, phytohormones and biological elicitors, or other chemicals (Figure 9, Additional File 4). The data analyzed in this comparison are based on the Affymetrix ATH1 platform, were retrieved from the 'Genevestigator' database , and were processed as described in Additional File 4. Based on this meta-analysis, the DBM-affected transcription factors and other signalling components identified in our experiments can be divided into two large classes: (i) those that are responsive to a wide range of different forms of biotic and abiotic stress [Figure 9, group I (induced) and group III (repressed)]; and (ii) those that are not effected by most other treatments of this comparison (Figure 9 and Additional File 4, group II).
Signalling components in group-I and group-III thus constitute candidates involved generally in stress perception and transduction. Group-I and group-III transcription factors include all DBM-affected AP2-EREBP genes, the bHLH PIF1, WRKY80, WRKY40, the C2H2 At5g04340, MYB15, and all DBM-effected NAC type transcription factors. Most of these genes were also found to be induced by other herbivores, although with notable exceptions (Figure 6). For example, the AP2-EREBP At2g2088 is rapidly and strongly induced by DBM feeding, but appears unaffected by the other herbivore treatments, while treatments with some bacterial or fungal pathogens, drought, and abscisic acid lead to transiently induced expression (Additional File 4).
Signalling components in group-II appear to be specific to DBM treatment, as they were not found induced in most other stress treatment (Figure 9). However, lack of induction in other treatments needs to be interpreted cautiously, since different expression profiling platforms were used. The transcription factors JAG (At1g68480) and bHLH121 (At3g19860), which are members of group II, were also found induced upon herbivory with P. rapae and M. persica e, respectively (Figure 6) but were found induced in only three and five, respectively, of the other 297 stress treatments (Figure 9 and Additional File 4). JAG had previously been shown to be necessary for the development and shaping of lateral organs such as leaves . Taken together, we identified a large set of signal transduction components that likely orchestrate a rapid and general response to a wide range of external stresses, but also a large set of signaling components that may mediate responses more specific to plant-herbivore responses.
Very few of the protein kinases, transcription factors and ubiquitin ligases that are affected by DBM feeding (Figure 9, Additional File 4) have well characterized functions. A notable exception is the transcription factor MYB34/ATR1 (At5g60890), which encodes a positive regulator of indole glucosinolate biosynthesis. ATR1 has previously been implied in insect-induced signalling  and is induced in particular late into the DBM feeding experiment. ATR1 was also found induced upon herbivory by P. rapae , but appears repressed upon treatment with F. occidentalis and M. persicae  (Figure 6). Three additional transcription factors, the two AP2-EREBP RAP2.6 (At1g43160) and At2g20880 as well as MYB15 (At3g23250), were previously identified as induced by wounding, methyl jasmonate, various pathogens, and several forms of abiotic stress .
The Arabidopsis transcriptome changes substantially in response to leaf feeding DBM larvae. Major DBM induced changes are involved in specialized (secondary) defence metabolism and in signalling. The DBM induced response shows considerable overlap with the response induced by other insect herbivores. A first large-scale meta-analysis of Arabidopsis microarray data obtained for a large number of biotic and abiotic interactions revealed groups of transcription factors and other signalling components that are similarly affected by multiple forms of biotic or abiotic stress including DBM feeding or, alternatively, appear more specifically responsive to DBM herbivory.
Plant and insect materials
Arabidopsis thaliana plants (ecotype Landsberg erecta) were grown in plastic pots (10 cm wide × 8 cm tall) containing Terra-lite Redi-earth (W.R. Grace and Co., Ajax, Ontario, Canada). Each pot contained four or five plants, which were grown in a growth chamber at 20°C constant temperature, 8 hr/16 hr Light/Dark photoperiod at 50–60% ambient humidity, for 8 to 9 weeks. Short day conditions prevented the onset of flowering and the plants were thus maintained in growth stage 1 (leaf production) with 13 to 15 rosette leaves larger than 1 mm (stage 1.13 to 1.14). Diamondback moth (DBM, Plutella xylostella) larvae were provided by Dr. Murray Isman (Faculty of Agricultural Sciences, University of British Columbia, Vancouver, Canada) and maintained on cabbage (Brassica oleracea) plants in a climate-controlled room at 25°C, 12 hr photoperiod with 50%-60% relative humidity. Two days before exposing A. thaliana plants to herbivore treatment, plants were transferred to a climate-controlled room (22°C, 50–60% humidity, 12 hr photoperiod). For insect treatment, seven DBM larvae (third to fifth instars) were placed on a group of four or five plants until time of harvest, for each time point separately. As control, Arabidopsis plants were maintained under the same condition except without exposure to DBM larvae. Rosette leaves from DBM-exposed and control plants were harvested at 1 h, 4 h, 8 h, 12 h and 24 h after onset of herbivory. For each treatment group and time point, all rosette leaves were harvested from the four or five plants per treatment group and flash frozen in liquid nitrogen.
Microarrays, RNA isolation, cDNA labelling, and microarray hybridization
The design and production of the A. thaliana 26,090 element 70-mer oligonucleotide microarray was previously described with oligonucleotides designed in collaboration with and purchased from Operon (Huntsville AL, USA) . All procedures for RNA isolation, RNA labelling and microarray hybridization were performed as described . Microarray experiments involved two independent biological replicates for each time point and treatment with each replicate consisting of four or five plants to provide enough plant material for RNA isolation. In addition, microarray hybridizations for each time point and treatment were performed with two technical replicates (dye-flip labelling) for each of the biological replicates for a total of four dual channel microarray hybridizations per time point comparing treatment with control.
Microarray data analysis
Microarrays were scanned with a ScanArray Express (Perkin Elmer, Woodbridge ON, Canada) scanner with laser power set to 95% and photo-multiplier-tube set to 54 to 64. We identified and quantified hybridization signals using ImaGene software (BioDiscovery, Marina Del Rey CA, USA). Grids were manually placed and spot finding was performed using the 'Auto adjust' spot function repeated three times. Spot finding was subsequently verified by visual inspection and manually adjusted when necessary. Poor spots were manually flagged (flag 1) and were not used in further data analyses. For all analyses, the median pixel intensities for each spot were used. All microarray expression data were submitted to the GEO database  under the accessions series GSE10681. Further analyses were performed using customized scripts for R and Bioconducter . For background correction, we defined the mean of the lowest 10% of spot intensities from a particular subgrid as the background for that subgrid. This mean was subtracted from each spot in the subgrid. Signal intensities that did not exceed the background plus 3 standard deviations thereof were defined as not detectable and were excluded from further analyses. We normalized using loess curves . For each array element, we first used the data from the four replicate array hybridizations (two biological replicates each with two technical replicates) for each time point and treatment to perform a paired Student's t-test using the Welch approximation to degrees of freedom. Subsequently, an analysis of variance (ANOVA) using data from all experimental samples (four normalized log2-expression ratios per time point for a total of 20 data points) was performed for each element. In order to assess the type I error rate, we calculated q-values estimating the false discovery rate based on the parametric p-values . Genes were first placed into one of three expression groups: Group A) up-regulated genes displaying a significant (p[t-test] < 0.05) and more than twofold higher signal in insect treated leaves compared to control leaves in at least one time point; Group B) down-regulated genes displaying a significant (p < 0.05) and more than twofold lower signal in insect treated compared to control plants in at least one time point; and Group C) genes with mixed expression using the same definition as in A and B. For clustering, mean log2-expression ratios for genes identified as differentially expressed (DE) in at least one time point were used. To derive at a reliable dataset, we defined genes as DE only if they met all of the following criteria: (1) significant (p[t-test] < 0.05) and more than twofold difference of transcript abundance between treatment and control for at least one time point, and (2) change of expression of more than twofold between at least two time points of the treatment time course with p(ANOVA) < 0.05. K-means clustering of DE genes was performed using Genesis v1.5  defining nine clusters with a maximum of 50 iterations. The normalized expression ratios and the results for all statistical analyses are summarized in Additional File 1.
Analyses of genes of interest
Gene lists containing selected genes of interests were retrieved from 'The Arabidopsis Resource Information database' (TAIR)  or from published gene family compilations. Lists of genes involved in "response to pathogens or wounding" were retrieved from TAIR (status December 2004). We selected genes placed in the GeneOntology (GO) categories "involved in" the "biological process" "response to pest, pathogen, or parasite" and/or "response to wounding". Because of the large number of genes of uncharacterized functions associated with the GO-terms, we only selected genes if they were also curator annotated based on experimental evidence by TAIR. Children terms of these GO categories were also included in the selection of these genes. Lists of complete gene families involved in "Arabidopsis secondary metabolism of glucosinolates, phenylpropanoids, or terpenoids" were compiled based on published surveys of the Arabidopsis genome [9–11, 32]. Complete lists of putative "protein kinases, transcription factors", and genes involved in the "26S proteasome pathway" were retrieved from the PlantsP (protein kinase families ), AtTFDB (transcription factor families ), and PlantsUBQ (26S proteasome gene families ) databases, respectively. For analysis of genes associated with "phytohormones or signalling compounds", curator annotated genes placed in the GO terms "involved in the metabolism of", "involved in the signalling mediated by", or "involved in the response to" the phytohormones or signalling molecules "auxin", "abscisic acid", "brassinosteroid", "cytokinin", "ethylene", "gibberellic acid", "jasmonic acid", and "salicylic acid" were retrieved from TAIR. Each of these gene lists was filtered to avoid multiple entries per list of the same gene.
Expression data for members of these gene lists that were found differentially expressed in our experiments according to the DE definition described above were visualized as heatmaps using Genesis v1.5 . To assess if any of the groups of genes associated with phytohormones or signalling compounds was significantly overrepresented in the insect-effected Arabidopsis transcriptome, the frequency of differentially expressed genes associated with each of the different phytohormones or signalling compounds in the groups A, B, and C (see above) was compared to the frequency of all genes associated with these categories represented on the microarray using a hypergeometric distribution. Phytohormone or signalling compound GO groups were defined as significantly over-represented in the group of differentially expressed genes, when p[hyper] < 0.01.
Quantitative real time RT-PCR (qRT-PCR)
Total RNA (15 μg) was digested with 15U DNAse in 1 × buffer (Invitrogen, Carlsbad CA, USA) for 15 min at room temperature. The reaction was stopped with EDTA (2.3 mM final concentration) and heat-inactivation (65°C, 10 min). RNA was precipitated with a 1/10 volume of 3 M sodium acetate and 2.5 volumes of ethanol at -20°C overnight, and subsequently pelleted at 20,000 × g for 30 min at 4°C. The precipitate was washed with 70% ethanol, centrifuged, and resuspended in RNAse free water to an approximate concentration of 1 μg/μl. Actual RNA concentration was determined spectrophotometrically. DNAse-treated total RNA (10 μg) of was used for reverse transcription with 0.27 μM T17VN primer, 0.15 mM dNTPs, 40 U RNAseOut, and 400 U SuperscriptII (Invitrogen) in 10 mM DTT and 1 × first strand buffer in a total volume of 40 μl. Prior to addition of enzymes the solution was heated to 65°C for 5 min and for primer annealing cooled to 42°C. Following an incubation at 42°C for 2.5 h, the RNA was degraded with 8μl 1 M sodium hydroxide for 15 min at 65°C, neutralized with 8 μl 1 M hydrochloric acid and buffered with 4 μl 1 M Tris-pH 7.5. Synthesized cDNA was purified using the Quiagen (Hilden, Germany) PCR-purification kit prior to quantitative PCR reaction using the ABI TaqMan system. The Custom TaqMan Gene Expression Assays (consisting of gene-specific TaqMan® MGB probe and primer sets, supplied as 20× stocks) used for quantitative real time PCR were from the Applied Biosystems (Foster City CA, USA) Custom Oligonucleotide Synthesis Service. The gene-specific probe and primer sets were designed using the Primer Express software from Applied Biosystems. Oligonucleotide sequences of all primer pairs and the respective probes are given in Additional File 5. Multiplex PCR reaction in triplicate (20 μl) containing cDNA equivalent to 100 ng RNA were performed in 96-well plates with the Opticon 2 (BioRad, Hercules CA, USA) using 1 μl of the 20× Custom TaqMan Gene Expression Assay (consisting 900 μmole of primer/250 μmole FAM-labeled probe, final concentration) for each of the genes analyzed; 2 μl primer pair and probe combination (300 μmole of primer/125 μmole of VIC-labeled probe) of the endogenous control (β-tubulin); 10 μl of 2× TaqMan Universal PCR Master Mix (containing all necessary components for fast reaction set-up for 5' nuclease reactions, including AmpliTaq Gold DNA polymerase, and AmpErase UNG). After an initial hold at 50°C for 2 min for the activation of AmpErase UNG, and denaturing at 95°C for 10 min, 40 cycles at 95°C for 15 sec and 60°C for 1 min, followed by fluorescence reading were performed. Data analysis was done according to a protocol by Applied Biosystems. Briefly, the baseline was set such that the amplification curve growth began at a cycle number that was greater than the highest baseline number. The threshold value was set within the exponential phase of the logarithmic scale amplification plot. Relative quantification of gene expression was calculated from the threshold cycle (CT) values for each replicate well on the reaction plate. Relative gene expression levels were calculated manually from the exported results file. Briefly, the VIC CT values were subtracted from the FAM CT values to calculate ΔCT for the control and samples at each of the time points for each of the transcription factors [ΔCT = CT (FAM) - CT (VIC)]. The ΔCT values for the triplicate wells of the Control samples at each time point for each of the transcription factor were averaged to obtain the mean ΔCT (Control). The mean ΔCT (Control) for a gene at a particular time point was then subtracted from the ΔCT values of this gene at that time point to calculate its ΔΔCT (Sample) [ΔΔCT (Sample) = ΔCT (Sample) - mean ΔCT (Control)]. The ΔΔCT (Sample) values were then averaged for the triplicate wells of each gene at each time point to normalize the target mRNA quantity to the internal control (β-Tubulin). The Average ΔΔCT (Sample) for each sample was then used to calculate the relative quantification values [2-mean ΔΔCT].
We gratefully acknowledge financial support (to JB) from the Natural Science and Engineering Research Council (NSERC) of Canada, Genome Canada, and Genome British Columbia. JB is a UBC Distinguished University Scholar and NSERC E.W.R. Steacie Memorial Fellow.
- De Vos M, Van Oosten VR, Van Poecke RMP, Van Pelt JA, Pozo MJ, Mueller MJ, Buchala AJ, Metraux JP, van Loon LC, Dicke M, Pieterse CMJ: Signal signature and transcriptome changes of Arabidopsis during pathogen and insect attack. Mol Plant Microbe In. 2005, 18: 923-937. 10.1094/MPMI-18-0923.View ArticleGoogle Scholar
- Kempema LA, Cui XP, Holzer FM, Walling LL: Arabidopsis transcriptome changes in response to phloem-feeding silverleaf whitefly nymphs. Similarities and distinctions in responses to aphids. Plant Physiol. 2007, 143: 849-865. 10.1104/pp.106.090662.PubMedPubMed CentralView ArticleGoogle Scholar
- Little D, Gouhier-Darimont C, Bruessow F, Reymond P: Oviposition by pierid butterflies triggers defense responses in Arabidopsis. Plant Physiol. 2007, 143: 784-800. 10.1104/pp.106.090837.PubMedPubMed CentralView ArticleGoogle Scholar
- Reymond P, Weber H, Damond M, Farmer EE: Differential gene expression in response to mechanical wounding and insect feeding in Arabidopsis. Plant Cell. 2000, 12: 707-720. 10.1105/tpc.12.5.707.PubMedPubMed CentralView ArticleGoogle Scholar
- Reymond P, Bodenhausen N, Van Poecke RMP, Krishnamurthy V, Dicke M, Farmer EE: A conserved transcript pattern in response to a specialist and a generalist herbivore. Plant Cell. 2004, 16: 3132-3147. 10.1105/tpc.104.026120.PubMedPubMed CentralView ArticleGoogle Scholar
- Van Poecke RM, Dicke M: Indirect defence of plants against herbivores: using Arabidopsis thaliana as a model plant. Plant Biol. 2004, 6: 387-401. 10.1055/s-2004-820887.PubMedView ArticleGoogle Scholar
- Wittstock U, Halkier BA: Glucosinolate research in the Arabidopsis era. Trends Plant Sci. 2002, 7: 263-270. 10.1016/S1360-1385(02)02273-2.PubMedView ArticleGoogle Scholar
- Raes J, Rohde A, Christensen JH, Van de Peer Y, Boerjan W: Genome-wide characterization of the lignification toolbox in Arabidopsis. Plant Physiol. 2003, 133: 1051-1071. 10.1104/pp.103.026484.PubMedPubMed CentralView ArticleGoogle Scholar
- Ehlting J, Mattheus N, Aeschliman DS, Li E, Hamberger B, Cullis IF, Zhuang J, Kaneda M, Mansfield SD, Samuels L, Ritland K, Ellis BE, Bohlmann J, Douglas CJ: Global transcript profiling of primary stems from Arabidopsis thaliana identifies candidate genes for missing links in lignin biosynthesis and transcriptional regulators of fiber differentiation. Plant J. 2005, 42: 618-640. 10.1111/j.1365-313X.2005.02403.x.PubMedView ArticleGoogle Scholar
- Aubourg S, Lecharny A, Bohlmann J: Genomic analysis of the terpenoid synthase ( AtTPS) gene family of Arabidopsis thaliana. Mol Genet Genomics. 2002, 267: 730-745. 10.1007/s00438-002-0709-y.PubMedView ArticleGoogle Scholar
- Lange BM, Ghassemian M: Genome organization in Arabidopsis thaliana: a survey for genes involved in isoprenoid and chlorophyll metabolism. Plant Mol Biol. 2003, 51: 925-948. 10.1023/A:1023005504702.PubMedView ArticleGoogle Scholar
- Li C, Zhao J, Jiang H, Wu X, Sun J, Zhang C, Wang X, Lou Y, Li C: The wound response mutant suppressor of prosystemin-mediated responses6 (spr6) is a weak allele of the tomato homolog of CORONATINE-INSENSITIVE1 (COI1). Plant Cell Physiol. 2006, 47: 653-663. 10.1093/pcp/pcj034.PubMedView ArticleGoogle Scholar
- Thines B, Katsir L, Melotto M, Niu Y, Mandaokar A, Liu GH, Nomura K, He SY, Howe GA, Browse J: JAZ repressor proteins are targets of the SCFCO11 complex during jasmonate signalling. Nature. 2007, 448: 661-6U2. 10.1038/nature05960.PubMedView ArticleGoogle Scholar
- Chini A, Fonseca S, Fernandez G, Adie B, Chico JM, Lorenzo O, Garcia-Casado G, Lopez-Vidriero I, Lozano FM, Ponce MR, Micol JL, Solano R: The JAZ family of repressors is the missing link in jasmonate signalling. Nature. 2007, 448: 666-6U4. 10.1038/nature06006.PubMedView ArticleGoogle Scholar
- Kusnierczyk A, Winge P, Midelfart H, Armbruster WS, Rossiter JT, Bones AM: Transcriptional responses of Arabidopsis thaliana ecotypes with different glucosinolate profiles after attack by polyphagous Myzus persicae and oligophagous Brevicoryne brassicae. J Exp Bot. 2007, 58: 2537-2552. 10.1093/jxb/erm043.PubMedView ArticleGoogle Scholar
- Global Crop Pest Identification and Information Services in Integrated Pest Management from The Cornell International Institute for Food, Agriculture and Development - Diamondback Moth (DBM). 1995, [http://www.nysaes.cornell.edu/ent/hortcrops/english/dbm.html]
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003, 100: 9440-9445. 10.1073/pnas.1530509100.PubMedPubMed CentralView ArticleGoogle Scholar
- Rhee SY, Beavis W, Berardini TZ, Chen G, Dixon D, Doyle A, Garcia-Hernandez M, Huala E, Lander G, Montoya M, Miller N, Mueller LA, Mundodi S, Reiser L, Tacklind J, Weems DC, Wu Y, Xu I, Yoo D, Yoon J, Zhang P: The Arabidopsis Information Resource (TAIR): a model organism database providing a centralized, curated gateway to Arabidopsis biology, research materials and community. Nucl Acids Res. 2003, 31: 224-228. 10.1093/nar/gkg076.PubMedView ArticleGoogle Scholar
- Consortium GO: The Gene Ontology Project in 2006. Nucl Acids Res. 2006, 34: D322-D326. 10.1093/nar/gkj021.View ArticleGoogle Scholar
- Creelman RA, Mulpuri R: The Oxylipin Pathway in Arabidopsis. The Arabidopsis book. Edited by: Somerville CR and Meyerowitz EM. 2002, Rockville, MD, American Society of Plant Biologists,, 1-24. 10.1199/tab.0012. [http://www.bioone.org/perlserv/?request=get-document&issn=1543-8120&volume=36&issue=1&page=1]Google Scholar
- Koo AJK, Chung HS, Kobayashi Y, Howe GA: Identification of a peroxisomal acyl-activating enzyme involved in the biosynthesis of jasmonic acid in Arabidopsis. J Biol Chem. 2006, 281: 33511-33520. 10.1074/jbc.M607854200.PubMedView ArticleGoogle Scholar
- Niyogi KK, Last RL, Fink GR, Keith B: Suppressors of trp1 fluorescence identify a new Arabidopsis gene, TRP4, encoding the anthranilate synthase beta subunit. Plant Cell. 1993, 5: 1011-1027. 10.1105/tpc.5.9.1011.PubMedPubMed CentralGoogle Scholar
- Lopukhina A, Dettenberg M, Weiler EW, Hollander-Czytko H: Cloning and characterization of a coronatine-regulated tyrosine aminotransferase from Arabidopsis. Plant Physiol. 2001, 126: 1678-1687. 10.1104/pp.126.4.1678.PubMedPubMed CentralView ArticleGoogle Scholar
- Field B, Cardon G, Traka M, Botterman J, Vancanneyt G, Mithen R: Glucosinolate and amino acid biosynthesis in Arabidopsis. Plant Physiol. 2004, 135: 828-839. 10.1104/pp.104.039347.PubMedPubMed CentralView ArticleGoogle Scholar
- Wittstock U, Halkier BA: Cytochrome P450 CYP79A2 from Arabidopsis thaliana L. catalyzes the conversion of L-phenylalanine to phenylacetaldoxime in the biosynthesis of benzylglucosinolate. J Biol Chem. 2000, 275: 14659-14666. 10.1074/jbc.275.19.14659.PubMedView ArticleGoogle Scholar
- Hull AK, Vij R, Celenza JL: Arabidopsis cytochrome P450s that catalyze the first step of tryptophan-dependent indole-3-acetic acid biosynthesis. Proc Natl Acad Sci U S A. 2000, 97: 2379-2384. 10.1073/pnas.040569997.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen SX, Glawischnig E, Jorgensen K, Naur P, Jorgensen B, Olsen CE, Hansen CH, Rasmussen H, Pickett JA, Halkier BA: CYP79F1 and CYP79F2 have distinct functions in the biosynthesis of aliphatic glucosinolates in Arabidopsis. Plant J. 2003, 33: 923-937. 10.1046/j.1365-313X.2003.01679.x.PubMedView ArticleGoogle Scholar
- Hansen CH, Du L, Naur P, Olsen CE, Axelsen KB, Hick AJ, Pickett JA, Halkier BA: CYP83B1 Is the oxime-metabolizing enzyme in the glucosinolate pathway in Arabidopsis. J Biol Chem. 2001, 276: 24790-24796. 10.1074/jbc.M102637200.PubMedView ArticleGoogle Scholar
- Glawischnig E, Hansen BG, Olsen CE, Halkier BA: Camalexin is synthesized from indole-3-acetaldoxime, a key branching point between primary and secondary metabolism in Arabidopsis. Proc Natl Acad Sci U S A. 2004, 101: 8245-8250. 10.1073/pnas.0305876101.PubMedPubMed CentralView ArticleGoogle Scholar
- Ljung K, Hull AK, Celenza J, Yamada M, Estelle M, Nonmanly J, Sandberg G: Sites and regulation of auxin biosynthesis in Arabidopsis roots. Plant Cell. 2005, 17: 1090-1104. 10.1105/tpc.104.029272.PubMedPubMed CentralView ArticleGoogle Scholar
- Mikkelsen MD, Naur P, Halkier BA: Arabidopsis mutants in the C-S lyase of glucosinolate biosynthesis establish a critical role for indole-3-acetaldoxime in auxin homeostasis. Plant J. 2004, 37: 770-777. 10.1111/j.1365-313X.2004.02002.x.PubMedView ArticleGoogle Scholar
- Kliebenstein DJ, Kroymann J, Mitchell-Olds T: The glucosinolate-myrosinase system in an ecological and evolutionary context. Curr Opin Plant Biol. 2005, 8: 264-271. 10.1016/j.pbi.2005.03.002.PubMedView ArticleGoogle Scholar
- Zhang Z, Ober JA, Kliebenstein DJ: The gene controlling the quantitative trait locus EPITHIOSPECIFIER MODIFIER1 alters glucosinolate hydrolysis and insect resistance in Arabidopsis. Plant Cell. 2006, 18: 1524-1536. 10.1105/tpc.105.039602.PubMedPubMed CentralView ArticleGoogle Scholar
- Hardtke CS, Berleth T: The Arabidopsis gene MONOPTEROS encodes a transcription factor mediating embryo axis formation and vascular development. EMBO J. 1998, 17: 1405-1411. 10.1093/emboj/17.5.1405.PubMedPubMed CentralView ArticleGoogle Scholar
- Gribskov M, Fana F, Harper J, Hope DA, Harmon AC, Smith DW, Tax FE, Zhang G: PlantsP: a functional genomics database for plant phosphorylation. Nucl Acids Res. 2001, 29: 111-113. 10.1093/nar/29.1.111.PubMedPubMed CentralView ArticleGoogle Scholar
- Davuluri RV, Sun H, Palaniswamy SK, Matthews N, Molina C, Kurtz M, Grotewold E: AGRIS: Arabidopsis gene regulatory information server, an information resource of Arabidopsis cis-regulatory elements and transcription factors. BMC Bioinformatics. 2003, 4: 25-10.1186/1471-2105-4-25.PubMedPubMed CentralView ArticleGoogle Scholar
- Paterson A, Bowers J, Kresovich S, Hash CT, Messing J, Peterson D, Schmutz J, Rokhsar D: Phytozome Sorghum bicolor. 2007, [http://www.phytozome.net/sorghum]Google Scholar
- Krysan PJ, Jester PJ, Gottwald JR, Sussman MR: An Arabidopsis mitogen-activated protein kinase kinase kinase gene family encodes essential positive regulators of cytokinesis. Plant Cell. 2002, 14: 1109-1120. 10.1105/tpc.001164.PubMedPubMed CentralView ArticleGoogle Scholar
- Batistic O, Kudla J: Integration and channeling of calcium signaling through the CBL calcium sensor/CIPK protein kinase network. Planta. 2004, 219: 915-924. 10.1007/s00425-004-1333-3.PubMedView ArticleGoogle Scholar
- Bögre L, Ökrész L, Henriques R, Anthony RG: Growth signalling pathways in Arabidopsis and the AGC protein kinases. Trends Plant Sci. 2003, 8: 424-431. 10.1016/S1360-1385(03)00188-2.PubMedView ArticleGoogle Scholar
- Celenza JL, Quiel JA, Smolen GA, Merrikh H, Silvestro AR, Normanly J, Bender J: The Arabidopsis ATR1 Myb transcription factor controls indolic glucosinolate homeostasis. Plant Physiol. 2005, 137: 253-262. 10.1104/pp.104.054395.PubMedPubMed CentralView ArticleGoogle Scholar
- Glickman MH, Ciechanover A: The ubiquitin-proteasome proteolytic pathway: destruction for the sake of construction. Physiol Rev. 2002, 82: 373-428.PubMedView ArticleGoogle Scholar
- Ralph SG, Yueh H, Friedmann M, Aeschliman D, Zeznik JA, Nelson CC, Butterfield YS, Kirkpatrick R, Liu J, Jones SJ, Marra MA, Douglas CJ, Ritland K, Bohlmann J: Conifer defence against insects: microarray gene expression profiling of Sitka spruce (Picea sitchensis) induced by mechanical wounding or feeding by spruce budworms (Choristoneura occidentalis) or white pine weevils (Pissodes strobi) reveals large-scale changes of the host transcriptome. Plant Cell Environ. 2006, 29: 1545-1570. 10.1111/j.1365-3040.2006.01532.x.PubMedView ArticleGoogle Scholar
- Ralph S, Oddy C, Cooper D, Yueh H, Jancsik S, Kolosova N, Philippe RN, Aeschliman D, White R, Huber D, Ritland CE, Benoit F, Rigby T, Nantel A, Butterfield YS, Kirkpatrick R, Chun E, Liu J, Palmquist D, Wynhoven B, Stott J, Yang G, Barber S, Holt RA, Siddiqui A, Jones SJ, Marra MA, Ellis BE, Douglas CJ, Ritland K, Bohlmann J: Genomics of hybrid poplar (Populus trichocarpa x deltoides) interacting with forest tent caterpillars (Malacosoma disstria): normalized and full-length cDNA libraries, expressed sequence tags, and a cDNA microarray for the study of insect-induced defences in poplar. Mol Ecol. 2006, 15: 1275-1297. 10.1111/j.1365-294X.2006.02824.x.PubMedView ArticleGoogle Scholar
- Mikkelsen MD, Hansen CH, Wittstock U, Halkier BA: Cytochrome P450 CYP79B2 from Arabidopsis catalyzes the conversion of tryptophan to indole-3-acetaldoxime, a precursor of indole glucosinolates and indole-3-acetic acid. J Biol Chem. 2000, 275: 33712-33717. 10.1074/jbc.M001667200.PubMedView ArticleGoogle Scholar
- Sandorf I, Holländer-Czytko H: Jasmonate is involved in the induction of tyrosine aminotransferase and tocopherol biosynthesis in Arabidopsis thaliana. Planta. 2002, 216: 173-179. 10.1007/s00425-002-0888-0.PubMedView ArticleGoogle Scholar
- Skoog F, Miller CO: Chemical regulation of growth and organ formation in plant tissues cultured in vitro. Symp Soc Exp Biol. 1957, 54: 118-130.Google Scholar
- Coenen C, Lomax TL: Auxin-cytokinin interactions in higher plants: old problems and new tools. Trends Plant Sci. 1997, 2: 351-356. 10.1016/S1360-1385(97)84623-7.PubMedView ArticleGoogle Scholar
- Eklöf S, Astot C, Sitbon F, Moritz T, Olsson O, Sandberg G: Transgenic tobacco plants co-expressing Agrobacterium iaa and ipt genes have wild-type hormone levels but display both auxin- and cytokinin-overproducing phenotypes. Plant J. 2000, 23: 279-284. 10.1046/j.1365-313x.2000.00762.x.PubMedView ArticleGoogle Scholar
- Cowan AK, Freeman M, Bjorkman PO, Nicander B, Sitbon F, Tillberg E: Effects of senescence-induced alteration in cytokinin metabolism on source-sink relationships and ontogenic and stress-induced transitions in tobacco. Planta. 2005, 221: 801-814. 10.1007/s00425-005-1489-5.PubMedView ArticleGoogle Scholar
- Thornburg RW, Li X: Wounding Nicotiana tabacum leaves causes a decline in endogenous indole-3-acetic acid. Plant Physiol. 1991, 96: 802-805.PubMedPubMed CentralView ArticleGoogle Scholar
- Schmelz EA, Engelberth J, Alborn HT, O'Donnell P, Sammons M, Toshima H, Tumlinson JH: Simultaneous analysis of phytohormones, phytotoxins, and volatile organic compounds in plants. Proc Natl Acad Sci U S A. 2003, 100: 10552-10557. 10.1073/pnas.1633615100.PubMedPubMed CentralView ArticleGoogle Scholar
- Baldwin IT, Zhang ZP, Diab N, Ohnmeiss TE, McCloud ES, Lynds GY, Schmelz EA: Quantification, correlations and manipulations of wound-induced changes in jasmonic acid and nicotine in Nicotiana sylvestris. Planta. 1997, 201: 397-404. 10.1007/s004250050082. [http://www.springerlink.com/content/x7bcdtruq56cvbbx/?p=c52d64ad4dfa45c89810af3222b44295&pi=1]View ArticleGoogle Scholar
- Cheong YH, Chang HS, Gupta R, Wang X, Zhu T, Luan S: Transcriptional profiling reveals novel interactions between wounding, pathogen, abiotic stress, and hormonal responses in Arabidopsis. Plant Physiol. 2002, 129: 661-677. 10.1104/pp.002857.PubMedPubMed CentralView ArticleGoogle Scholar
- Zimmermann P, Hirsch-Hoffmann M, Hennig L, Gruissem W: GENEVESTIGATOR. Arabidopsis microarray database and analysis toolbox. Plant Physiol. 2004, 136: 2621-2632. 10.1104/pp.104.046367.PubMedPubMed CentralView ArticleGoogle Scholar
- Ohno CK, Reddy GV, Heisler MG, Meyerowitz EM: The Arabidopsis JAGGED gene encodes a zinc finger protein that promotes leaf tissue development. Development. 2004, 131: 1111-1122. 10.1242/dev.00991.PubMedView ArticleGoogle Scholar
- Chen W, Provart NJ, Glazebrook J, Katagiri F, Chang HS, Eulgem T, Mauch F, Luan S, Zou G, Whitham SA, Budworth PR, Tao Y, Xie Z, Chen X, Lam S, Kreps JA, Harper JF, Si-Ammour A, Mauch-Mani B, Heinlein M, Kobayashi K, Hohn T, Dangl JL, Wang X, Zhu T: Expression profile matrix of Arabidopsis transcription factor genes suggests their putative functions in response to environmental stresses. Plant Cell. 2002, 14: 559-574. 10.1105/tpc.010410.PubMedPubMed CentralView ArticleGoogle Scholar
- NAtional Centre for Biotechnology Information (NCBI): GEO Gene Expression Omnibus. 2007, [http://www.ncbi.nlm.nih.gov/geo/]Google Scholar
- Team RDC: R: A language and environment for statistical computing. 2006, Vienna, Austria, R Foundation for Statistical Computing, [http://www.R-project.org]Google Scholar
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucl Acids Res. 2002, 30: e15-10.1093/nar/30.4.e15.PubMedPubMed CentralView ArticleGoogle Scholar
- Sturn A, Quackenbush J, Trajanoski Z: Genesis: cluster analysis of microarray data. Bioinformatics. 2002, 18: 207-208. 10.1093/bioinformatics/18.1.207.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.