- Research article
- Open Access
Transcript profiling of cytokinin action in Arabidopsis roots and shoots discovers largely similar but also organ-specific responses
BMC Plant Biology volume 12, Article number: 112 (2012)
The plant hormone cytokinin regulates growth and development of roots and shoots in opposite ways. In shoots it is a positive growth regulator whereas it inhibits growth in roots. It may be assumed that organ-specific regulation of gene expression is involved in these differential activities, but little is known about it. To get more insight into the transcriptional events triggered by cytokinin in roots and shoots, we studied genome-wide gene expression in cytokinin-treated and cytokinin-deficient roots and shoots.
It was found by principal component analysis of the transcriptomic data that the immediate-early response to a cytokinin stimulus differs from the later response, and that the transcriptome of cytokinin-deficient plants is different from both the early and the late cytokinin induction response. A higher cytokinin status in the roots activated the expression of numerous genes normally expressed predominantly in the shoot, while a lower cytokinin status in the shoot reduced the expression of genes normally more active in the shoot to a more root-like level. This shift predominantly affected nuclear genes encoding plastid proteins. An organ-specific regulation was assigned to a number of genes previously known to react to a cytokinin signal, including root-specificity for the cytokinin hydroxylase gene CYP735A2 and shoot specificity for the cell cycle regulator gene CDKA;1. Numerous cytokinin-regulated genes were newly discovered or confirmed, including the meristem regulator genes SHEPHERD and CLAVATA1, auxin-related genes (IAA7, IAA13, AXR1, PIN2, PID), several genes involved in brassinosteroid (CYP710A1, CYP710A2, DIM/DWF) and flavonol (MYB12, CHS, FLS1) synthesis, various transporter genes (e.g. HKT1), numerous members of the AP2/ERF transcription factor gene family, genes involved in light signalling (PhyA, COP1, SPA1), and more than 80 ribosomal genes. However, contrasting with the fundamental difference of the growth response of roots and shoots to the hormone, the vast majority of the cytokinin-regulated transcriptome showed similar response patterns in roots and shoots.
The shift of the root and shoot transcriptomes towards the respective other organ depending on the cytokinin status indicated that the hormone determines part of the organ-specific transcriptome pattern independent of morphological organ identity. Numerous novel cytokinin-regulated genes were discovered which had escaped earlier discovery, most probably due to unspecific sampling. These offer novel insights into the diverse activities of cytokinin, including crosstalk with other hormones and different environmental cues, identify the AP2/ERF class of transcriptions factors as particularly cytokinin sensitive, and also suggest translational control of cytokinin-induced changes.
Cytokinin is a plant hormone regulating numerous developmental and physiological activities [1, 2]. Many molecular details of cytokinin metabolism and signal transduction have been discovered during recent years. It is known that cytokinin acts locally as well as over a distance, that it is synthesized and degraded in different root and shoot tissues and that its signal is transduced through a complex two-component system. In Arabidopsis thaliana, three membrane-located sensor histidine kinases – AHK2, AHK3 and CRE1/AHK4 – perceive the signal and transmit it through a phosphorelay via phosphotransmitter proteins (AHP) to the nucleus where B-type response regulators (ARR), a class of transcription factors, become activated and induce the transcription of numerous target genes [3–5]. It has been shown that the B-type ARRs are involved in mediating most if not all of the transcriptional responses to cytokinin [6–8]. Genome-wide analyses of the Arabidopsis transcriptome have identified many genes that show a rapid up- or downregulation of their steady-state mRNA level in response to cytokinin [6, 7, 9–18], reviewed in . Among these are numerous genes encoding transcription factors suggesting that transcriptional cascades operate downstream of cytokinin and are involved in realizing the biological output reactions. Indeed, functional studies have confirmed that beside the B-type ARRs, other factors are involved in mediating the transcriptional response to cytokinin. These include a subgroup of the ERF/AP2 transcription factor family named CYTOKININ RESPONSE FACTOR (CRF) , members of the GeBP transcription factor family  and GATA22 [21, 22].
An interesting feature of cytokinin is its opposite activity in regulating growth of roots and shoots. Cytokinin is a negative regulator of root growth and branching. Plants with a reduced cytokinin status – generated either by reduction of the endogenous cytokinin content [23–26] or by disruption of cytokinin signalling [6–8, 27] show enhanced root growth. In the root apical meristem, cytokinin promotes the transition from cells to differentiation [24, 28] acting through the auxin-response factor IAA3/SHY2 . In contrast, cytokinin retards cellular differentiation in the shoot apical meristem by interacting with the WUS/CLV pathway, as well as by a separate pathway [30–33]. Similarly, cytokinin inhibits root branching [34, 35], but stimulates the growth of lateral shoots . Moreover, a high cytokinin:auxin ratio favours shoot induction in undifferentiated callus tissue, while a low cytokinin:auxin ratio leads to the formation of roots . In addition to these developmental changes, it has been shown that cytokinin-deficiency has different consequences on the primary metabolism of root and shoot tissue .
Different activities of cytokinin in the two organs could be partly realized through a different transcriptional response. Shoots and roots express specific sets of genes . Here we describe the results of a genome-wide expression profiling of roots and shoots of seedlings treated with cytokinin for different periods of time. In addition, we compared the root and shoot transcriptomes of wild-type plants with those of cytokinin-deficient plants overexpressing the cytokinin oxidase/dehydrogenase gene CKX1 under the control of the 35S promoter . We looked for similarities and differences of cytokinin-dependent gene expression in roots and shoots to identify genes with organ-specific regulation, as well as to assign a root- or shoot-specific response to already known cytokinin-regulated genes. The results of this transcript profiling approach are a starting point for further research into the differential responses of roots and shoots to cytokinin.
Results and discussion
In order to search for short-term and long-term effects on the transcriptome, we treated 5-d-old wild-type seedlings with cytokinin for 15 min (named BA15), 2 h (BA120), and 18 h (BA1080) and harvested roots and shoots separately. To contrast the cytokinin induction, we analyzed root and shoot samples of cytokinin-deficient 35S:CKX1 expressing seedlings of the same age (named CKX1). These seedlings contain about 38% of the cytokinin levels of wild-type seedlings in the roots and about 12% in the shoots . The reference samples for the whole experiment were roots and shoots of mock-treated wild-type seedlings (BA0). mRNA preparations were obtained from two biological replicates and subjected to microarray analysis as described in Methods. The hybridization strategy used for this study yielded absolute expression values instead of the fold-changes usually obtained from two-colour microarrays. To identify cytokinin-regulated genes, we considered only those transcripts which were detected on at least 25% of the microarrays (i. e. 8 microarrays for cytokinin induction and 4 microarrays for cytokinin deficiency), because inclusion of weakly expressed genes caused a large number of false positives as determined by re-examination by qRT-PCR (data not shown). In total, 8,549 and 7,557 genes were above the cutoff on the arrays used for the analysis of gene regulation by cytokinin treatment or cytokinin-deficiency, respectively. These genes were then classified as described in Methods and Additional file 1: Figure S1. Genes with more than 2.5-fold change of expression compared to the control samples are listed in Additional file 2: Table S1 (cytokinin induction, 1,450 genes) and Additional file 3: Table S2 (cytokinin deficiency, 1,473 genes), respectively. Selected examples of these genes are listed in Table 1 (cytokinin induction) and Table 2 (cytokinin deficiency) and will be commented on further below.
Global expression analysis by principal component analysis
In order to get an insight into the global structure of the complex dataset of this microarray investigation, we performed a principal component analysis (PCA) on a reduced dataset (see Methods) using the steady-state mRNA levels determined by the application of a common reference hybridization strategy (see Methods). This type of analysis applied to the samples of a microarray experiment typically results in clusters of samples with a similar expression pattern, thereby revealing the main factors that lead to changes in gene expression. Additionally, we performed a support tree (ST) clustering, which served as a confirmation of the PCA.
Cytokinin deficiency has a large and specific influence on the transcriptome
The results of the PCA show a separation of the dataset into distinct clusters (Figure 1b-e). As expected, the largest distance was found between root and shoot samples confirming the functionality of the PCA (Figure 1b, c, e). The samples of the cytokinin-deficient plants were separated from the wild-type samples on a different axis than the organ effect (Figure 1b – d). Additionally, they form their own cluster in ST clustering (Figure 1a). The large distance of the cytokinin-deficient samples from the wild-type clusters means that cytokinin deficiency causes a global transcript change to an extent comparable to the difference between roots and shoots.
Taken together, PCA analysis groups the samples into four clusters — wild-type roots, wild-type shoots, cytokinin-deficient roots, and cytokinin-deficient shoots — each of which occupies one quadrant of the clustering space (Figure 1b-e).
The cytokinin status influences the organ-specific transcriptome pattern
Closer examination of the distribution of the samples within the clustering space revealed that a long-term cytokinin treatment of roots and constitutively lowered cytokinin content in shoots causes distinct changes of the gene expression pattern shifting the location of the root and shoot samples in the PCA towards the respective other organ.
Firstly, the long-term cytokinin-treated root sample is located more closely to the shoot coordinates than any other root sample (Figure 1b, c, e). This change indicates that after 18 h of cytokinin treatment, a shift in gene expression towards a shoot-like pattern has started without gross visible developmental changes having occurred in the tissue. Next, we looked for affected genes and identified these based on their differential expression between untreated roots and shoots (≥ 2.5-fold), their altered expression level after 18 h of cytokinin treatment (≥ 2.5-fold), and by their p-values of the organ, cytokinin and/or interaction effects (≤ 0.03). The transcripts of 106 genes are regulated in this fashion (Figure 2, Additional file 4: Table S3). Among the 59 genes that are upregulated towards a more shoot-like transcript level in long-term cytokinin-treated roots, we identified genes encoding chloroplast-localized proteins as the major overrepresented gene ontology class (70% of the regulated genes) (Additional file 5: Figure S2a, Additional file 4: Table S3a). These include genes encoding proteins necessary for plastid proliferation, such as the plastid developmental protein DAG (AT2G33430) and the chloroplast division protein FtsZ (AT2G36250). Conversely, 47 predominantly root-expressed genes, such as AT5G47450 encoding the tonoplast intrinsic protein TIP2;3 transporting ammonium across the tonoplast (Table 1, Additional file 4: Table S3b, Additional file 6: Figure S3a, b, Additional file 7: Figure S4a, b) , are repressed by long-term cytokinin treatment in roots. However, no overrepresented functional GO category could be found among this group of genes. In comparison, only 30 genes acquired a more root-like expression upon cytokinin treatment of shoots (Figure 2) and no shift of the respective sample towards root-like expression was observed in the PCA (Figure 1b-e). No GO category could be found to be overrepresented among these genes.
Secondly, a similar, albeit opposite, developmental shift effect was observed in the shoot sample of cytokinin-deficient plants, which is located more closely to the root coordinates than any other shoot sample (Figure 1a - c, e). This global shift of the transcript profile to a more root-like pattern suggests that, beside the visible reduction of shoot growth, a low cytokinin status reduces the “shootyness” of the global transcription pattern. We identified genes displaying a corresponding expression pattern as described above. 246 genes show a more root-like transcript level in shoots of cytokinin-deficient seedlings (Figure 2, Additional file 8: Table S4). Also in this case, the set of 169 genes that was downregulated to a more root-like transcript level contains a high fraction (52%) of chloroplast-related genes (Additional file 5: Figure S2b, Additional file 8: Table S4b). The gene encoding the PsbR subunit of photosystem II (AT1G79040), for example, is downregulated to about one tenth of its expression level in shoots of cytokinin-deficient plants. Conversely, 77 otherwise predominantly shoot-expressed genes are stronger expressed in cytokinin-deficient roots (Figure 2, Additional file 8: Table S4a), such as PIN2, encoding a root-specific auxin efflux carrier discussed further below, but no functional GO category was found to be overrepresented among them.
Interestingly, 63 of all regulated genes showing an enhanced transcript level in cytokinin-treated roots were also repressed in cytokinin-deficient shoots and vice versa indicating that these are particularly cytokinin-sensitive genes (Table 3). 34 of the 44 genes enhanced by cytokinin in roots and repressed in cytokinin-deficient shoots encode plastid proteins with functions e.g. in the light reaction or the calvin cycle as well as several plastidal ribosomal proteins (Table 3). Examples of the 19 genes repressed by cytokinin in roots and enhanced in cytokinin-deficient shoots encode the predominantly root-expressed sodium transporter HKT1 discussed further below, and the signalling protein CBL-interacting protein kinase 23 (CIPK23, AT1G30270).
The fact that cytokinin deficiency results in downregulation of chloroplast genes in green tissue, whereas excess cytokinin causes chloroplast genes expressed in the root is consistent with a primary role of cytokinin in regulating plastid development . It is noteworthy that only a fraction (3.5%) of the ca. 3,000 nuclear genes encoding plastidic proteins  were affected by an altered cytokinin status, suggesting that cytokinin acts only on part of the genes and/or that there is a sequential order of cytokinin-regulated gene expression of which we detected only a part at the analyzed time point. Consistent with our data others found also a strong increase of genes encoding plastid proteins in cytokinin-treated root callus tissue, albeit at much later time points [10, 43]. The expression of a comparatively small number of genes was shifted towards a more shoot-like profile in cytokinin-deficient roots (Figure 2), and correspondingly no shift of the cytokinin-deficient root sample towards the shoot coordinates was observed in the PCA (Figure 1b-e).
Taken together, these findings indicate that the cytokinin status apparently determines a more root- or shoot-like transcriptomic pattern that is at least partially independent of morphologically discernable altered organ identity. We used the category of nuclear genes encoding plastid proteins as a diagnostic tool to assess the “shootyness” of the transcriptome but it is impossible to determine whether the transcriptomic shift in “shootyness” is cause or consequence of a shift in organ identity. However, several studies dissecting the temporal pattern of gene expression during cytokinin-induced shoot development from root calli found that shoot marker genes start to be upregulated only several days after the onset of the hormonal treatment [10, 44–46] supporting the notion that during the time period analyzed here no shoot identity was obtained by the cytokinin-treated root tissue.
Similarity of the overall transcriptional response of roots and shoots to cytokinin
The effect of cytokinin treatment at different time points was not resolved in the PCA analysis, because the large organ and genotype effects prevailed. Therefore, we removed the largest effect — the organ effect — by normalizing the root and shoot samples for each experimental condition as described in Methods. As a result, support tree clustering showed distinct root-shoot-pairs at each experimental condition indicating that the overall transcriptional response to cytokinin is similar in roots and shoots (Figure 1f). The two early time points of cytokinin induction showed the closest relation, while the long-term treated samples separated from the other wild-type samples.
PCA analysis using this dataset normalized for the organ effect resolved two temporal phases of the transcriptomic response to cytokinin. For the early time points of cytokinin induction, the separation is approximately parallel to the z-axis and for the later time points parallel to the y-axis, suggesting that the late response affects a different set of genes than the early response (Figure 1j). Both root and shoot samples follow a similar pattern (shown as dashed lines in Figure 1j), indicating that the transcriptomic response to cytokinin treatment is similar in both organs. This conclusion was confirmed by an additional PCA without the cytokinin-deficient samples, where the arrangement of the wild-type samples in one plane remained stable (data not shown).
Given the fundamental differences of the growth response of roots and shoots to cytokinin, we initially hypothesized that this could be reflected by large differences at the transcriptome level. In fact, cytokinin response genes are regulated by the hormone in a context-dependent manner , which is likely due to the presence of multiple cis-regulatory elements in their promoters. Indeed, the only known cytokinin response element is not present in the promoters of all cytokinin responsive genes , suggesting that a variety of possibilities exist to regulate a gene by cytokinin. It might well have been possible that large specific subsets of genes are regulated by cytokinin in an organ-specific manner. However, that appears not to be the case. Therefore it is remarkable and to some extent unexpected that the global transcriptional responses in both organs appear to be rather similar. This indicates that it is not the differential bulk expression determining the rather different cytokinin-dependent growth responses in both organs. This would also imply that cytokinin signaling to the nucleus is to a large extent similar and that organ-specific features are limited to few changes, which regulate the growth response. It would be very interesting to find out how these comparably small differences can drive completely different growth responses in roots and shoots. Furthermore, because of these relatively small differences one can expect a relatively low number of genes that are regulated by cytokinin in an organ-specific manner.
It should be noted that the absolute expression level of the genes, although showing qualitatively similar cytokinin responses, can be quite different in roots and shoots (Additional file 2: Table S1 and Additional file 3: Table S2). For instance, the transcripts encoding the auxin transporter PIN2 or the tonoplast intrinsic protein TIP2;3 are more than hundred-fold more abundant in the root than in the shoot, while the phosphate transporter gene PHT2;1 shows a much higher expression level in the shoot (Additional file 6: Figure S3, Additional file 7: Figure S4). Thus the differential response of the tissues to cytokinin may be affected by additional factors that regulate the absolute abundances of the different messages.
Chronic cytokinin deficiency and cytokinin induction affect different sets of genes
The PCA of the normalized dataset (Figure 1g-j) shows that the CKX samples separate from the wild-type samples along the x-axis while the effects of cytokinin treatment on the transcriptome are distributed over the y- and z-axes and remain within the same plane (Figure 1g-j, purple markings). This means that the transcriptomic changes resulting from chronic cytokinin deficiency are not simply the opposite of those resulting from early or late cytokinin treatment, but that another (third) set of genes is affected. Consequently, we analyzed cytokinin-deficiency separately from the dataset derived from the cytokinin-induced plants.
Regulation of cytokinin metabolism and signalling genes
To search for potential homeostatic mechanisms within the cytokinin system, we studied the influence of cytokinin treatment and cytokinin deficiency on the transcript level of cytokinin metabolism and signalling genes (Figure 3). Unfortunately, the CATMA microarray lacks a part of these cytokinin-related genes and expression of another part was below the detection limit. Therefore, only a limited number of the relevant genes could be analyzed. A subset of these showed significant differences in organ specificity and cytokinin response of transcript abundance.
Among the cytokinin-synthesizing IPT genes, IPT5 showed a higher expression in the shoot under all conditions, while IPT7 (and less strongly also IPT9) was consistently more highly expressed in the root (Figure 3a). IPT5 was strongly downregulated in long-term cytokinin-treated roots and upregulated in cytokinin-deficient shoots (see also Figure 3a). In contrast, a strong upregulation in long-term cytokinin-treated roots was detected for IPT9, and IPT3 was apparently neither differentially expressed nor regulated by cytokinin. The gene encoding the cytokinin hydroxylase CYP735A2 was strongly upregulated by cytokinin in a root-specific fashion (see below). Among the cytokinin-degrading CKX genes, strong differences were found for CKX4 and CKX5, confirming previous data . Both genes were expressed more strongly in the root than in the shoot and both were upregulated by cytokinin and downregulated under conditions of cytokinin deficiency, the response in the root being stronger than in the shoot (Figure 3a). CKX4 has been consistently described as a primary target gene of B-type response regulator ARR1 . In contrast to these two genes, CKX6 was expressed more strongly in the shoot and was hardly regulated by the cytokinin status in our study. Please note that the strong upregulation of CKX1 found in cytokinin-deficient roots results from detection of the transgene. Of the putative cytokinin transporter genes, PUP4 and PUP11 were more strongly expressed in the root and both showed a transient upregulation in the shoot but not in the root in response to cytokinin (Figure 3a).
The known higher expression of the cytokinin receptor gene CRE1/AHK4 in the root [51, 52] and its induction by cytokinin  was confirmed. The transcript level of the two B-type Arabidopsis response regulator genes (ARR1, ARR2) were more strongly expressed in the root and – as expected – not regulated by cytokinin (Figure 3a). In contrast, all known primary cytokinin response genes encoding A-type response regulators (ARR4, ARR6, ARR7, ARR9, and ARR16) detected were equally responsive to cytokinin in both roots and shoots (Figure 3a, Additional file 6: Figure S3a, Additional file 7: Figure S4a, b). Figure 3b shows the reduced expression levels of four exemplary A-type ARR genes (ARR3, ARR4, ARR6, and ARR7) in 35S:CKX1 transgenic roots and shoots as compared to the corresponding wild-type tissues and their induction pattern following cytokinin treatment in roots and shoots of the wild type (see also Additional file 6: Figure S3c and Additional file 7: Figure S4c, d). Despite a generally similar response to an altered cytokinin status, the kinetics and amplitudes of changes in steady state transcript levels differ for the four genes.
Together, these data confirm that the transcript abundance of several cytokinin metabolism and signalling genes depends on the cytokinin status [6, 9, 13, 15, 17] indicating regulatory loops between cytokinin metabolism and signalling. The analysis also showed that the reaction patterns are not uniform but complex, with differences in kinetics and organ specificity.
Known and novel genes responding to cytokinin
In the following sections we will analyse separately the sets of genes that show altered expression patterns following either a treatment of seedlings with cytokinin or under conditions of constitutive cytokinin deficiency (Additional file 2: Table S1 and Additional file 3: Table S2). The organ-specific regulation of some known cytokinin response genes will be addressed and an arbitrary selection of novel cytokinin-regulated genes, which may have escaped previous discovery due to unspecific sampling, is discussed. We will explore similar and differential changes of steady-state transcript levels in both organs, as well as changes that are restricted to either roots or shoots. This organizing principle is followed, although in part genes belonging to a similar context (e.g. related to auxin or brassinosteroid, meristem regulation, nutrient transport, etc.) were categorized into different groups of organ-specific cytokinin-regulation.
Cytokinin-induced changes of transcript abundance
1,450 of the 8,549 genes (i.e. ~17%) that were detected by the arrays were responsive to a short-term or long-term cytokinin treatment (Figure 4). Out of these, 670 genes could be assigned unambiguously to one specific category as follows: the majority, i.e. 396 genes, were regulated in both organs with either a similar or differential expression pattern. The former genes (355) were regulated in both organs in the same direction, while the latter genes (41) were regulated in both organs in the opposite direction. 177 genes were regulated by cytokinin in a root-specific fashion, and 97 genes were regulated in a shoot-specific fashion. A large number of genes (780) were regulated significantly by cytokinin under at least one of the experimental conditions, but lacked a significant organ effect. Most of these genes showed a different steady-state transcript level at only one time point (Additional file 2: Table S1). They were not further analysed, as the focus of this work was on organ specificity. 7099 genes were classified as non-regulated under our conditions (Figure 4 and Additional file 2: Table S1f).
Assignment of organ-specificity to known cytokinin response genes
Firstly, we examined whether we could assign a root- or shoot-specific response to known cytokinin-regulated genes which had been discovered previously by genome-wide microarray analyses [6, 9, 13, 15, 17]. In fact, a large number of these genes were detected by the CATMA array and most, but not all, of them were found to be regulated in both organs (Table 4 and Additional file 9: Table S5).
A noteworthy example of a root-specific cytokinin response gene encodes the cytokinin hydroxylase CYP735A2 which responded strongly to cytokinin within minutes. Its transcript level in roots after 120 min of cytokinin treatment was strongly increased, while it was non-responsive in shoot tissue (Table 4, Additional file 6: Figure S3a). Significantly, the basal level of CYP735A2 expression was similar in both roots and shoots (Additional file 6: Figure S3b, Additional file 7: Figure S4a, b). This regulation is interesting because the hydroxylase catalyzes the conversion of isopentenyl riboside phosphates to trans-zeatin (tZ) riboside phosphates , which are inactive precursors of tZ-type cytokinins. tZ-type cytokinins are the main cytokinins transported in the xylem from roots to shoots [5, 55]. The exclusive induction of CYP735A2 by cytokinin in roots, which probably results in a locally increased production of tZ-type cytokinins, lends support to tZ-type cytokinins having a distinct role in root-to-shoot communication [55–57].
An example of a shoot-specific regulated known cytokinin response gene is CDKA;1, encoding a cyclin-dependent kinase gene. CDKA;1 is an important cell cycle regulator and is involved in a variety of developmental processes. Because cells expressing more CDKA;1 are more competent for cell division, it is thought to be a prerequisite for post-embryonic cell division [58, 59]. It was shown recently that the CDK activity level in the shoot apical meristem is important for cellular differentiation , a process that is under cytokinin control [30–32]. CDKA;1 regulation in the shoot by cytokinin could be part of the activities orchestrated by cytokinin in the shoot apical meristem.
The organ-specificity of some other known cytokinin response genes is also described in the following sections.
Genes with a similar response in roots and shoots upon cytokinin treatment
Of the 355 genes displaying a similar cytokinin response in roots and shoots, 133 genes showed an increased transcript level following cytokinin treatment, while the transcript level of 222 genes was reduced (Additional file 2: Table S1). 66 genes showed a rapid change of transcript abundance within 15 min and thus belong to the immediate-early response genes. The vast majority of the changes became apparent 120 min or 1080 min after the treatment. Interestingly, the transcript level of many genes that was increased after 120 min remained increased also after 1080 min. As already mentioned above, it should be noted that there are substantial differences in the basal expression levels of numerous of these genes in roots and shoots which could be of functional relevance as well. An example is the CLV1 gene discussed below, which is approx. 4-fold higher expressed in the untreated wild-type shoots compared to the roots (Additional file 6: Figure S3b, Additional file 7: Figure S4a, b).
Among the upregulated developmental genes, a particularly interesting one encodes the CLAVATA (CLV) formation protein SHEPHERD (SHD, AT4G24190). SHD is a member of the HSP90 protein family and has been suggested to act as a chaperone for the correct folding of nascent CLV1 and/or CLV3-like polypeptides in the endoplasmatic reticulum . Consistently, shd mutant plants are phenotypically reminiscent of clv mutants. The CLV1 and CLV3 proteins are part of the regulatory loop controlling meristem size and activity in concert with WUSCHEL (WUS) . Cytokinin regulates the expression of CLV1 and WUS, thereby controlling meristematic activity [30, 31, 33]. Transcriptional regulation of the SHD gene possibly represents an additional activity of cytokinin in regulating meristematic activity.
Members of the ERF/AP2 family were noticeably frequent among the cytokinin-responsive transcription factor genes. The transcript levels of eight members of the family (AT5G07580, AT1G77640, AT5G61600, AT1G28370, AT3G15210, AT5G51190, AT5G47220, and AT2G44840) were strongly downregulated by the hormone, but none was upregulated (Table 1, Additional file 6: Figure S3a, Additional file 7: Figure S4a, b). This shows that, in addition to the ERF/AP2 subclass of cytokinin-regulated factors (CRF) with its current six members , other members of other subclades of the ERF/AP2 family are also cytokinin-regulated (see also below).
Interestingly, numerous ribosomal genes were identified as cytokinin responsive in this study as well as in previous studies [9, 13, 63]. Several developmental processes react sensitively to mutations in ribosomal genes, including embryogenesis, root growth and leaf development, underpinning the functional relevance of the specific assembly of ribosomal proteins in addition to their fundamental role in translation, possibly because they are required to translate specific proteins [64–66]. According to the Gene Ontology database of May 2011, the Arabidopsis genome encodes 474 ribosomal proteins [67, 68] and at least 87 transcripts of them responded to cytokinin (Additional file 10: Table S6). The cytokinin-responsive ribosomal genes encode proteins of both plastidic and cytoplasmic localisation. Noteworthy, most of the transcripts of nuclear genes encoding plastidic ribosomal proteins reflect the organ shift discussed above, i.e. induction in long-term cytokinin-treated roots and downregulation in cytokinin-deficient shoots. Ribosomal genes encoding cytosolic proteins show diverse cytokinin response patterns. The altered composition of ribosomes and the translational changes that this may cause are a largely unexplored part of the cytokinin response.
The flavonol synthase gene FLS1 was strongly upregulated, particularly in the shoot. FLS1 is part of the flavonol biosynthesis pathway  and it is likely that its strong upregulation is part of the induction of anthocyanin accumulation by cytokinin . Cytokinin and flavonoids are involved in regulating several developmental processes: both, for example, are positive regulators of nodule organogenesis  but have opposite roles in the control of shoot branching [72–74]. Flavonoids are inhibitors of auxin transport and thus may mediate crosstalk between cytokinin and auxin action. In any case, there are multiple links between cytokinin and flavonoids and effects on the FLS1 transcript level could be part of that. FLS1 was also found among the strongly upregulated genes in cytokinin-deficient plants (see below), illustrating nicely that cytokinin deficiency does not simply induce the opposite changes caused by cytokinin treatment, as already shown by the PCA discussed above.
Five members of the nodulin MtN21 gene family are downregulated by cytokinin: AT2G28120, AT4G08290, AT4G01450, AT4G08300, and AT3G25190 (Additional file 6: Figure S3a, Additional file 7: Figure S4a, b). The role of nodulin MtN21-like genes is not yet understood, but the presence of seven transmembrane domains and structural homologies with bacterial multidrug exporters might suggest a role in a transport function. It was proposed that the maize nodulin MtN21-like gene is probably involved in the transport of a component related to vascular tissue assembly , which is interesting in view of the role of cytokinin in vascular development .
Several other regulated genes can be associated with known functions of cytokinin, such as crosstalk with other hormones including auxin (e.g. downregulation of IAA7 and IAA13) and gibberellin (downregulation of the gibberellin 2-oxidase gene AT1G02400) or its functions in mineral uptake (e.g. IRT2, ZIP4) (Additional file 2: Table S1).
Root-specific changes in gene expression upon cytokinin treatment
Among the 177 genes that were regulated only in roots, the strongest differences in expression level were detected for a gene encoding a member of the legume lectin family protein, its transcript abundance increased >150-fold (Additional file 2: Table S1). As has been discussed above, the largest specific group of regulated genes encoded plastid proteins, and a particularly interesting gene showing root-specific regulation was cytokinin hydroxylase gene CYP735A2. A root-specific regulation by cytokinin was assigned to two other cytochrome P450 genes for which cytokinin-regulation has been previously described [14, 15]. One of them, AT2G34490, encodes the C22 sterol desaturase CYP710A2 involved in brassinosteroid metabolism. It catalyzes the final step in the biosynthesis of brassicasterol and stigmasterol [77, 78]. This is interesting as cytokinin and brassinosteroids have been recently shown to act in common pathways in regulating root architecture . The molecular basis for this crosstalk is not known but it is conceivable that regulation of metabolism genes may be part of it. Noticeably, the related gene encoding CYP710A1 and the brassinosteroid biosynthesis gene DIM/DWF are deregulated in cytokinin-deficient plants (see below) indicating that brassinosteroid metabolism is regulated by cytokinin on different levels.
The largest fold-change in transcript abundance between roots of the cytokinin-deficient and long-term cytokinin-treated seedlings was noted for the LBD4 gene, encoding a member of the LOB domain protein family  and already noted before in relation to cytokinin . The function of LBD4 is unknown, but members of this plant-specific family of DNA-binding transcription factors are involved in various developmental processes, such as leaf polarity, and floral or inflorescence architecture in different plant species [80–84], as well as lateral root formation and the regulation of anthocyanin and nitrogen metabolism .
Again, a number of genes were identified that are relevant for auxin action, e.g. encoding the IAA-amido synthase DFL1 (AT5G54510) and PIN4 (AT2G01420), suggesting that cytokinin-auxin crosstalk is realized in many different ways in a tissue-specific fashion. Similarly, genes encoding proteins involved in mineral uptake are also regulated in a root-specific fashion, including iron transporter gene IRT1 as well as phosphate (PHT2;1, AT3G26570) and sulfate transporter genes (Additional file 2: Table S1, Additional file 6: Figure S3a, Additional file 7: Figure S4a, b).
Shoot-specific changes in gene expression upon cytokinin treatment
97 genes were exclusively regulated in the shoot, of which 23 were upregulated (up to 6-fold), and 74 were downregulated, in both cases most consistently after 120 and 1080 min. One of the most prominent genes that was rapidly and strongly downregulated in the shoot is CLV1 (Table 1, Additional file 6: Figure S3a, Additional file 7: Figure S4a, b). Its negative regulation by cytokinin has been previously reported [30, 31, 33], but the shoot-specific nature of this was not known. Similarly, cytokinin regulation of CRF3, encoding a member of the subfamily B-5 of the ERF/AP2 transcription factor family (AT5G53290), had been shown before , but not its shoot-specificity (Table 1). Another member of the ERF/AP2 transcription factor family (subfamily B-3), ERF5 (AT5G47230), is one of the most strongly downregulated genes in the shoot (Table 1). ERF5 is an activator of GCC box-dependent transcription and is also regulated by ethylene and abiotic stress .
The gene encoding the digalactosyldiacylglycerol synthase 1 (DGD1, AT3G11670) was upregulated almost 6-fold by cytokinin treatment. DGD1 catalyzes the synthesis of the galactolipid digalactosyldiacylglycerol (DGDG), requiring monogalactosyldiacylglycerol (MGDG) as a substrate, which is synthesized by mongalactosyldiacylglycerol synthase (MGD) [87, 88]. Both galactolipids constitute the predominant and functionally relevant component of thylakoid membranes . Previously, it was shown that the accumulation of MGDG in cucumber seedlings required cytokinin in addition to light, and that both factors play a co-operative role . The response of DGD1 shown here indicates that the cytokinin regulation also comprises this step of galactolipid synthesis, adding another example of the influence of the hormone on chloroplast function.
The serin/threonine protein kinase gene PINOID (PID), encoding a regulator of the auxin response, is downregulated approximately 3-fold. PID phosphorylates the auxin efflux carrier proteins PIN, thereby changing their apical-basal polar localization and thus increasing basipetal polar auxin transport [92–96]. Although it is not known whether transcriptional regulation of PID plays a role in regulating its kinase activity, it is possible that the downregulation of PID by cytokinin causes an increase in basipetal auxin transport. In any case, in whatever functional context, regulation of PID by cytokinin could be a mechanism in addition to the known regulation of PIN genes [29, 97, 98] that contributes to the manifold regulatory influences of cytokinin on auxin function.
Genes with differential expression patterns in root and shoot upon cytokinin treatment
41 genes were upregulated in shoots and downregulated in roots or vice versa. The existence of such a differential organ-specific response pattern could be, for example, achieved by transcription factors that act depending on the cellular or developmental context either as activator or repressor . However, given the low number of genes in this class and the frequently transient nature of their differential regulation, this does not seem to be a general feature of differential cytokinin activity in roots and shoots. Two genes encoding cytochrome P450 enzymes belong in this class and show altered transcript levels at several time points (Additional file 2: Table S1). One of these genes encodes CYP83A1, also called REF2 (AT4G13770), catalyzing the initial conversion of aldoximes to thiohydroximates in the synthesis of glucosinolates not derived from tryptophan. It has also a role in auxin homeostasis [100–102]. F5H2 (AT5G04330) encoding a ferulate 5-hydroxylase involved in lignification  is regulated in the opposite way, namely up in the root and down in the shoot. The functional significance of this type of regulation is not known.
Genes differentially expressed under conditions of constitutive cytokinin deficiency
In cytokinin-deficient tissues of 35S:CKX1 plants, 7,557 genes were detected, of which 1,473 (~20%) showed an altered transcript abundance compared to the wild type. 718 of these genes could be assigned to one specific category. The major effect of cytokinin deficiency was manifested in the shoot, where 639 genes were differently expressed. 363 of these genes were exclusively affected in this organ and 277 were affected in the root as well. A smaller number of 79 genes were regulated only in the root and 2 genes were identified showing differential regulation in root and shoot (Additional file 3: Table S2, Figure 4). One outcome of the analysis was the relatively small overlap of genes affected by cytokinin-deficiency or cytokinin-treatment as described above. There are examples of opposite gene regulation in both conditions, e.g. the A-type ARR genes, which are upregulated by cytokinin and downregulated under cytokinin deficiency. Conversely, there are several cases of a change of transcript level in the same direction under both conditions (e.g. FLS1, several AP2/ERF transcription factor genes) illustrating that opposite hormonal conditions may have a similar effect on transcript abundance. The relative shift of the shoot expression pattern towards a more root-like pattern, especially concerning genes encoding plastid proteins, was particularly noteworthy in the light of the effects of long-term cytokinin treatment on the roots described above. In addition, the data show that roots and shoots clearly react with different trancriptomic changes to cytokinin deficiency. Exceptions are general indicators of the cytokinin status, such as the A-type ARR genes, which are lowered in both organs (Figure 3). The alteration of the transcript level of several regulated genes may contribute to the establishment of the cytokinin deficiency syndrome.
Genes with similar reaction in root and shoot under constitutive cytokinin deficiency
275 genes showed a similar response in roots and shoots: 113 thereof were upregulated and 162 were downregulated (Additional file 3: Table S2). There are several regulated genes that provide insight into the biological functions of cytokinin, such as the strong downregulation of nitrate reductase genes NR1 and NR2 reflecting the functional link between cytokinin and nitrogen usage .
A very strongly altered expression level was shown by the sodium transporter gene HKT1, which was upregulated 24-fold in roots and 16-fold in shoots (Additional file 3: Table S2, Additional file 6: Figure S3c, Additional file 7: Figure S4c, d). The main function of HKT1 is the control of Na+ entry into plant roots and the root/shoot distribution of Na+[105, 106]. Thus, HKT1 has an important role in mediating salt tolerance, as has been shown in ecotypes as well as by forward genetics, and it has been utilized to generate salt-tolerance using a transgenic approach [107–110]. The regulation found here is consistent with the previous finding that HKT1;1 is repressed by cytokinin treatment but showed significantly elevated expression in the B-type response regulator double mutant arr1 arr12. These data suggested that cytokinin, acting through the transcription factors ARR1 and ARR12, regulates sodium accumulation in the shoots by controlling the expression of HKT1 in the roots. Cytokinin receptors AHK2 and AHK3 were consistently identified as negative regulators of the osmotic stress response and receptor mutants as well as mutants with a reduced cytokinin content were found to be salt resistant [112, 113]. Another cytokinin-regulated transporter gene was ZIP1 (AT3G12750), which encodes a Zn2+ transporter and was upregulated 3.2- and 4.4-fold in shoots and roots, respectively . Its upregulation by cytokinin deficiency is in accord with the finding that cytokinin-deficient plants accumulate higher amounts of Zn2+ in their shoot tissue .
Among the transcription factor genes, four ERF/AP2 transcription factor genes (AT1G28370, AT3G16770, AT5G05410, and AT3G15210) are downregulated to about one third of their expression level (Additional file 6: Figure S3c, Additional file 7: Figure S4c, d). Interestingly, three of them are identical to the cytokinin-regulated AP2/ERFs mentioned above (see Table 1), indicating again that increase and decrease of the cytokinin status may have similar consequences.
MYB12 is another strongly regulated transcription factor gene. MYB12 belongs to subgroup 7 of the R2R3-MYB family and is a flavonol-specific activator of flavonoid biosynthesis . Together with its two other subgroup members, MYB11 and MYB111, it strongly activates the promoters of genes encoding chalcone synthase (CHS), flavanone 3-hydroxylase (F3H), flavonol synthase (FLS1), and – to a lesser extent – chalcone flavanone isomerase (CHI), which are all involved in the formation of flavonols . The upregulation of MYB12 consistently resulted in an enhanced transcript level of CHS (AT5G13930) and FLS1 (Additional file 3: Table S2). We hypothesize that cytokinin-deficient plants contain higher levels of flavonols and/or flavonoids, which have a variety of functions in plant stress defence. Interestingly, auxin and ethylene induce flavonol accumulation through partly identical transcriptional networks .
Two well known auxin-related genes, PIN2 and AXR1, were upregulated under cytokinin-deficiency. PIN2 belongs to those genes which show a much higher shoot expression in response to cytokinin deficiency (55-fold increase compared to wild-type shoots) while it is normally predominantly expressed in roots (Table 2 and Additional file 3: Table S2, Additional file 6: Figure S3c, d, Additional file 7: Figure S4, Additional file 6: Figure S3, d). PIN2 encodes a member of the auxin efflux carrier proteins [118–121]. Its regulation by cytokinin has been reported in the context of root meristem size regulation [97, 122], but it might also be involved in co-ordinating the activities of cytokinin and auxin in regulating root gravitropism . AXR1, which was upregulated about 5-fold in both organs, encodes part of a dimeric E1 ubiquitin-activating enzyme of the ubiquitin-proteasome protein degradation pathway [124–127] mainly involved in regulating the auxin response, but also other signalling pathways dependent on ubiquitin-dependent protein degradation [128–133]. The upregulation of AXR1 may enhance the potential for neddylation of proteins to be targeted for degradation by the proteasome. The Aux/IAA proteins, transcriptional repressors of auxin-inducible genes, are primary targets of this pathway. Higher expression of AXR1 under cytokinin-deficiency suggests an enhanced auxin signalling in cytokinin-deficient plants. The identification of AXR1 as a regulatory target of cytokinin adds another potential regulatory level to cytokinin-auxin interaction, which, because of the involvement of AXR1 in auxin signalling in many tissues, might be of general relevance for the whole plant.
Root-specifically altered transcript levels under constitutive cytokinin deficiency
The transcript level of 79 genes was changed in a root-specific fashion: 35 genes were upregulated and 44 genes were downregulated (Additional file 3: Table S2). The strongest upregulation in roots showed a gene for an amino acid transporter family protein (AT5G65990) which was increased 18-fold (Table 2, Additional file 6: Figure S3a, b, Additional file 7: Figure S4a, b).
The transcript of the CLV3-related gene CLE1 (AT1G73165) was upregulated 5.7-fold under cytokinin deficiency, however, the p-value (0.049) was slightly above our threshold (0.03). There is not much known about the actual function of CLE1. However, CLE1 complemented the clv3 mutant well and also requires CLV1 for signalling, and thus it may be considered a functional homologue . The role of CLV3 and related CLE peptides in the shoot apical meristem has been investigated in great detail . However, they are equally important in root meristem maintenance. Overexpression, as well as the exogenous application of several CLE peptides – including CLV3 – to roots, results in the consumption of root meristematic cells [136–140] and both CLV3 and CLE40, its functional orthologue in roots, promote differentiation in the distal root meristem. CLE40 forms a ligand-receptor pair with ACR4 to regulate WOX5 expression, resulting in a similar regulatory loop for the balance between proliferation and differentiation as the CLV/WUS system in the shoot meristem . It could be that the root-specific upregulation of CLE1 under cytokinin deficiency is functionally relevant for the altered differentiation behaviour of cytokinin-deficient roots [24, 26].
Shoot-specifically altered transcript levels under constitutive cytokinin deficiency
363 genes were deregulated in a shoot-specific fashion in 35S:CKX1 transgenic seedlings, which is a more complex response compared to roots. 189 of the genes were upregulated and many showed very strong changes (~80-fold) of their expression level and a strong shift towards root-like expression (Additional file 3: Table S2) which has caused the shift of the cytokinin-deficient root sample towards a more root-like position in the PCA analysis discussed above. Interestingly, the three most strongly downregulated genes and two more genes with strongly reduced transcript levels encode AP2/ERF transcription factors (Additional file 3: Table S2, Additional file 6: Figure S3c, Additional file 7: Figure S4c, d), four of which are identical to those downregulated by cytokinin treatment (AT5G51190, AT5G61600, AT2G44840, and AT5G47220) mentioned above.
Several components of the light signalling network, in particular linked to phytochrome A, were found to be deregulated in cytokinin-deficient shoots. SPA1 and COP1 belong to the upregulated genes, while the phyA transcript itself is downregulated to one third of its level in wild-type plants (Additional file 3: Tables S2b and e). PhyA is involved in two light responses: the very-low-fluence response (VLFR) and the far-red-light-dependent high-irradiance response (HIR), the former promoting seed germination with very low quantities of red or far-red light and the latter promoting photomorphogenesis . SPA1 (an atypical bHLH transcription factor) and COP1 (an E3 ubiquitin ligase) are negative regulators of partial aspects of both VLFR and HIR. COP1 and SPA1 are able to form a complex which targets PhyA for proteolytic degradation, thereby promoting photomorphogenesis . The upregulation of the COP1 and SPA1 transcripts by cytokinin may contribute to the rapid degradation of PhyA in its P fr form and the transcriptional downregulation of PhyA by cytokinin is likely to attenuate PhyA accumulation. Together, the gene regulation is consistent with altered light responses of cytokinin-deficient plants  and the known links between cytokinin and light responses [1, 2].
Finally, the transcript of the DIM/DWF1 gene encoding an enzyme involved in brassinosteroid biosynthesis  was downregulated to one third of its wild-type level in cytokinin-deficient shoots. dim/dwf1 mutants are brassinosteroid-deficient and show the typical dwarfed phenotype. It would be interesting to study in what way an altered brassinosteroid level is causally involved in establishing the cytokinin deficiency syndrome of which reduced shoot growth is a hallmark . In any case, together with the regulation of CYP710A and CYP710B genes mentioned above our results suggest an influence of cytokinin on brassinosteroid biosynthesis. Reciprocal influences of the two hormones on the expression of their metabolism genes have been observed before; however, the regulation of brassinosteroid metabolism genes by cytokinin was described as being ambiguous and not allowing a clear prediction whether this results in an enhanced or decreased brassinosteroid level .
This study reports hitherto unknown differences of transcriptomic changes in root and shoot tissues in response to cytokinin treatment or cytokinin deficiency. One interesting result was the shift of long-term treated wild-type root samples to a more “shoot-like” profile and the opposite shift of cytokinin-deficient shoot samples towards “root-like” expression, both primarily affecting nuclear genes encoding plastid proteins. It would be interesting to analyze whether the whole transcriptome reacts in a rheostat-like fashion to a change in cytokinin status. There are some noteworthy examples of tissue-specific cytokinin responses, the most prominent probably being the cytokinin hydroxylase gene CYP735A2, which is strongly upregulated in root tissue but lacks a response in shoot tissue. This implies that strong negative regulators suppress the cytokinin response of CYP735A2 in shoot tissue. Numerous novel cytokinin response genes were discovered which were missed previously because of unspecific sampling and it can be anticipated that studies with an even higher spatial resolution will continue to discover novel cytokinin-regulated genes in their functional context . Several novel regulatory influences of cytokinin on metabolism and signalling genes of other hormonal regulators of plant development such as auxin and brassinosteroids were found. Similarly, novel connections between cytokinin and environmental cues operating through transcriptional changes were unravelled supporting the notion that the hormone is involved in regulating a vast array of plant functions . Together these findings generate new hypotheses that need to be tested experimentally in order to fit them into the emerging picture of regulatory transcriptional networks [145, 147].
Growth, treatment and harvest of Arabidopsis thaliana
Arabidopsis thaliana Col-0 seeds and seeds of ectopic overexpressors of the cytokinin oxidase/dehydrogenase gene CKX1 were surface-sterilized with chlorine gas . Seedlings were grown before treatment for five days in liquid medium (half-strength MS medium, 1 g L–1 sucrose, 0.5 g L–1 MES, pH 5.7) under a 16 h light/8 h dark regime at 22°C. The whole seedlings were then treated with either 5 μM 6-benzyladenine (BA) for 15, 120 or 1080 min, or for 120 min in a solvent control, which were both added to the medium. To minimize circadian effects on gene expression, the treatment was timed in such a way that the harvest could be carried out for all samples between 6 h and 8 h after the onset of light. To achieve this, both the induction and the harvesting process for each sample covered the whole period of 2 h. At the end of the incubation period, roots and shoots were harvested separately by taking each seedling individually with a soft forceps at the hypocotyl, breaking off the roots of shock-frozen seedlings at the inner wall of a tube swimming in liquid nitrogen and collecting the remaining shoot in another tube. In this way, the breaking point was always very close to the root-hypocotyl junction and very pure root and shoot samples were obtained. The roots and shoots of at least 2,400 seedlings were pooled for one sample. Two samples were collected per condition.
RNA extraction, preparation of fluorescently labelled probes, hybridisation of microarrays and generation of raw data was performed as described previously . Each sample was hybridized onto two microarrays, resulting in four microarrays per condition. The CATMA V2.3 microarray  of the Microarray Facility of the University Hospital Leuven (MAF, Leuven, Belgium, ArrayExpress accession no. A-MEXP-120) was used for this study. We followed the hybridization strategy used in the European Compendium of Arabidopsis Gene Expression project (http://www.cagecompendium.org), employing a Cy3-labeled mixture of oligonucleotides complementary to the secondary primers for probe amplification as a common reference sample and Cy5 labeling for the biological samples. This allowed for determination of the steady-state mRNA levels as opposed to the ratios usually obtained by two-color microarrays. The raw data are accessible at ArrayExpress (http://www.ebi.ac.uk/arrayexpress/) under the accession E-CAGE-111.
The raw data were normalized and pre-processed using the pre-processing pipeline developed with Bioconductor  in R  during the course of the CAGE project [152, 153]. After pre-processing, linear models were fitted for each gene [154, 155]. Linear modelling was performed using two factors: (1) organ with the levels “root” and “shoot”, and (2) cytokinin with the levels “CKX1”, “BA0”, “BA15”, “BA120”, and “BA1080”. From the individual p-values calculated in linear modelling, FDR-corrected p-values (“q-values”) were calculated . Genes which were above background on less than four arrays were called “not expressed” and excluded from the analysis. Genes with less than ten spots above background exhibited a very high percentage of false positives, so only genes with at least ten spots above background were chosen for further analysis. The resulting dataset of (log) expressions, (log) ratios, p-values, q-values, and spots above background (AbB) was used to generate all further data derived.
Interaction plots were generated in R and reformatted with standard image editing software (CorelDraw).
Colour-coded expression maps (Figure 2a) were generated using Mapman [157, 158]. Tables with log-ratios for the conditions compared were generated in R for this purpose and exported as tab-delimited text files in a Mapman-compatible format. Genes which were detected on less than two microarrays in the respective more highly expressed condition were flagged as not detected. The Mapman-generated images were exported and the colour-coded boxes were aligned using standard image editing software. Genes were categorized according to their expression patterns by an algorithm shown in Additional file 1: Figure S1a. Two sorting algorithms adapted to the datasets of the cytokinin-induced and the cytokinin-deficient plants, respectively, were implemented in R. This was done by splitting the dataset up into the subsets “root-specific”, “shoot-specific”, “differential”, “similar”, “uncategorized”, and “not regulated” using an algorithm outlined in Additional file 1: Figure S1.
Principal Component Analysis (PCA) based on the steady-state mRNA levels was carried out using the implementation in TIGR MEV 4  with the standard settings using a restricted dataset of 911 genes with ≥ 10 spots above background and an FDR-corrected p-value of ≤ 0.0001. The graphic output was re-formatted using standard image editing software.
GO term enrichment was performed using the Gene Ontology website http://amigo.geneontology.org/cgi-bin/amigo/term_enrichment. The graphic result was downloaded from the website and adapted to printing format. Normalization of the dataset for the organ effect was carried out using the following formula on the log-transformed expression values for the five shoot samples of each gene:
norm, normalized value; raw, native value; , average of the native values.
Real-time quantitative PCR analysis
Gene expression profiles were verified by real-time quantitative PCR analysis using the following protocol. Each reaction contained 1× reaction buffer (10× stock, 160 mM (NH4)2SO4, 1 M Tris–HCl, pH 8.3, 0.1% Tween-20), 2 mM MgCl2, 100 μM each dNTP, 0.1× SYBR Green I (Fluka, St. Louis, MO, USA, 10,000× stock), 50 nM 5-carboxy-X-rhodamine triethylamine salt (ROX, Sigma, St. Louis, MO, USA), 50 – 300 nM of each primer, 10 mU/μL immolase (BioLine, Luckenwalde, Germany) in a total volume of 20 μL. Each sample-gene combination was run as triple replicates. Real-time PCR runs were performed in a 7500 Fast system (Applied Biosystems), controlled and evaluated by the 7500 Software V2.0.1 RC. In order to exclude primers binding to unspecific sequences, they were designed using GENOPLANTE™S.P.A.D.S (http://www.psb.ugent.be/SPADS/)  or NCBI Primer-BLAST (http://www.ncbi.nlm.nih.gov/tools/primer-blast/) . The primer sequences are listed in Additional file 11: Table S7.
Argueso CT, Ferreira FJ, Kieber JJ: Environmental perception avenues: the interaction of cytokinin and environmental response pathways. Plant Cell Env. 2009, 32: 1147-1160.
Werner T, Schmülling T: Cytokinin action in plant development. Curr Opin Plant Biol. 2009, 12: 527-538.
Heyl A, Schmülling T: Cytokinin signal perception and transduction. Curr Opin Plant Biol. 2003, 6: 480-488.
Müller B, Sheen J: Advances in cytokinin signaling. Science. 2007, 318: 68-69.
Sakakibara H: CYTOKININS: Activity, biosynthesis, and translocation. Annu Rev Plant Biol. 2006, 57: 431-449.
Argyros RD, Mathews DE, Chiang Y-H, Palmer CM, Thibault DM, Etheridge N, Argyros DA, Mason MG, Kieber JJ, Schaller GE: Type B response regulators of Arabidopsis play key roles in cytokinin signaling and plant development. Plant Cell. 2008, 20: 2102-2116.
Heyl A, Ramireddy E, Brenner WG, Riefler M, Allemeersch J, Schmülling T: The transcriptional repressor ARR1-SRDX suppresses pleiotropic cytokinin activities in Arabidopsis. Plant Physiol. 2008, 147: 1380-1395.
Ishida K, Yamashino T, Yokoyama A, Mizuno T: Three type-B response regulators, ARR1, ARR10 and ARR12, play essential but redundant roles in cytokinin signal transduction throughout the life cycle of Arabidopsis thaliana. Plant Cell Physiol. 2008, 49: 47-57.
Brenner WG, Romanov GA, Köllmer I, Bürkle L, Schmülling T: Immediate-early and delayed cytokinin response genes of Arabidopsis thaliana identified by genome-wide expression profiling reveal novel cytokinin-sensitive processes and suggest cytokinin action through transcriptional cascades. Plant J. 2005, 44: 314-333.
Che P, Gingerich DJ, Lall S, Howell SH: Global and hormone-induced gene expression changes during shoot development in Arabidopsis. Plant Cell. 2002, 14: 2771-2785.
Hoth S, Ikeda Y, Morgante M, Wang X, Zuo J, Hanafey MK, Gaasterland T, Tingey SV, Chua NH: Monitoring genome-wide changes in gene expression in response to endogenous cytokinin reveals targets in Arabidopsis thaliana. FEBS Lett. 2003, 554: 373-380.
Kiba T, Aoki K, Sakakibara H, Mizuno T: Arabidopsis response regulator, ARR22, ectopic expression of which results in phenotypes similar to the wol cytokinin-receptor mutant. Plant Cell Physiol. 2004, 45: 1063-1077.
Kiba T, Naitou T, Koizumi N, Yamashino T, Sakakibara H, Mizuno T: Combinatorial microarray analysis revealing Arabidopsis genes implicated in cytokinin responses through the His → Asp phosphorelay circuitry. Plant Cell Physiol. 2005, 46: 339-355.
Lee DJ, Park JY, Ku SJ, Ha YM, Kim S, Kim MD, Oh MH, Kim J: Genome-wide expression profiling of ARABIDOPSIS RESPONSE REGULATOR 7(ARR7) overexpression in cytokinin response. Mol Genet Genomics. 2007, 277: 115-137.
Rashotte AM, Carson SD, To JP, Kieber JJ: Expression profiling of cytokinin action in Arabidopsis. Plant Physiol. 2003, 132: 1998-2011.
Rashotte AM, Mason MG, Hutchison CE, Ferreira FJ, Schaller GE, Kieber JJ: A subset of Arabidopsis AP2 transcription factors mediates cytokinin responses in concert with a two-component pathway. Proc Natl Acad Sci USA. 2006, 103: 11081-11085.
Taniguchi M, Sasaki N, Tsuge T, Aoyama T, Oka A: ARR1 directly activates cytokinin response genes that encode proteins with diverse regulatory functions. Plant Cell Physiol. 2007, 48: 263-277.
Yokoyama A, Yamashino T, Amano Y, Tajima Y, Imamura A, Sakakibara H, Mizuno T: Type-B ARR transcription factors, ARR10 and ARR12, are implicated in cytokinin-mediated regulation of protoxylem differentiation in roots of Arabidopsis thaliana. Plant Cell Physiol. 2007, 48: 84-96.
Brenner WG, Ramireddy E, Heyl A, Schmülling T: Gene regulation by cytokinin. Frontiers in Plant Science. 2012, 3: 8-
Chevalier F, Perazza D, Laporte F, Le Henanff G, Hornitschek P, Bonneville JM, Herzog M, Vachon G: GeBP and GeBP-like proteins are noncanonical leucine-zipper transcription factors that regulate cytokinin response in Arabidopsis. Plant Physiol. 2008, 146: 1142-1154.
Köllmer I, Werner T, Schmülling T: Ectopic expression of different cytokinin-regulated transcription factor genes of Arabidopsis thaliana alters plant growth and development. J Plant Physiol. 2011, 168: 1320-1327.
Richter R, Behringer C, Müller IK, Schwechheimer C: The GATA-type transcription factors GNC and GNL/CGA1 repress gibberellin signaling downstream from DELLA proteins and PHYTOCHROME-INTERACTING FACTORS. Genes Dev. 2010, 24: 2093-2104.
Miyawaki K, Matsumoto-Kitano M, Kakimoto T: Expression of cytokinin biosynthetic isopentenyltransferase genes in Arabidopsis: tissue specificity and regulation by auxin, cytokinin, and nitrate. Plant J. 2004, 37: 128-138.
Werner T, Motyka V, Laucou V, Smets R, van Onckelen H, Schmülling T: Cytokinin-deficient transgenic Arabidopsis plants show multiple developmental alterations indicating opposite functions of cytokinins in the regulation of shoot and root meristem activity. Plant Cell. 2003, 15: 2532-2550.
Werner T, Motyka V, Strnad M, Schmülling T: Regulation of plant growth by cytokinin. Proc Natl Acad Sci U S A. 2001, 98: 10487-10492.
Werner T, Nehnevajova E, Köllmer I, Novak O, Strnad M, Krämer U, Schmülling T: Root-specific reduction of cytokinin causes enhanced root growth, drought tolerance, and leaf mineral enrichment in Arabidopsis and tobacco. Plant Cell. 2010, 22: 3905-3920.
Riefler M, Novak O, Strnad M, Schmülling T: Arabidopsis cytokinin receptor mutants reveal functions in shoot growth, leaf senescence, seed size, germination, root development, and cytokinin metabolism. Plant Cell. 2006, 18: 40-54.
Dello Ioio R, Linhares FS, Scacchi E, Casamitjana-Martinez E, Heidstra R, Costantino P, Sabatini S: Cytokinins determine Arabidopsis root-meristem size by controlling cell differentiation. Curr Biol. 2007, 17: 678-682.
Dello Ioio R, Linhares FS, Sabatini S: Emerging role of cytokinin as a regulator of cellular differentiation. Curr Opin Plant Biol. 2008, 11: 23-27.
Bartrina I, Strnad M, Werner T, Schmülling T: Cytokinin regulates the activity of reproductive meristems, flower organ size, ovule formation and thus seed yield in Arabidopsis thaliana. Plant Cell. 2011, 23: 69-80.
Gordon SP, Chickarmane VS, Ohno C, Meyerowitz EM: Multiple feedback loops through cytokinin signaling control stem cell number within the Arabidopsis shoot meristem. Proc Natl Acad Sci U S A. 2009, 106: 16529-16534.
Leibfried A, To JP, Busch W, Stehling S, Kehle A, Demar M, Kieber JJ, Lohmann JU: WUSCHEL controls meristem function by direct regulation of cytokinin-inducible response regulators. Nature. 2005, 438: 1172-1175.
Lindsay DL, Sawhney VK, Bonham-Smith PC: Cytokinin-induced changes in CLAVATA1 and WUSCHEL expression temporally coincide with altered floral development in Arabidopsis. Plant Sci. 2006, 170: 1111-1117.
Cary AJ, Liu W, Howell SH: Cytokinin action is coupled to ethylene in its effects on the inhibition of root and hypocotyl elongation in Arabidopsis thaliana seedlings. Plant Physiol. 1995, 107: 1075-1082.
Laplaze L, Benkova E, Casimiro I, Maes L, Vanneste S, Swarup R, Weijers D, Calvo V, Parizot B, Herrera-Rodriguez MB, et al: Cytokinins act directly on lateral root founder cells to inhibit root Initiation. Plant Cell. 2007, 19: 3889-3900.
Faiss M, Zalubìlová J, Strnad M, Schmülling T: Conditional transgenic expression of the ipt gene indicates a function for cytokinins in paracrine signaling in whole tobacco plants. Plant J. 1997, 12: 401-415.
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.
Werner T, Holst K, Pors Y, Guivarc'h A, Mustroph A, Chriqui D, Grimm B, Schmulling T: Cytokinin deficiency causes distinct changes of sink and source parameters in tobacco shoots and roots. J Exp Bot. 2008, 59: 2659-2672.
Schmid M, Davison TS, Henz SR, Pape UJ, Demar M, Vingron M, Schölkopf B, Weigel D, Lohmann JU: A gene expression map of Arabidopsis thaliana development. Nat Genet. 2005, 37: 501-506.
Loqué D, Ludewig U, Yuan L, von Wirén N: Tonoplast intrinsic proteins AtTIP2;1 and AtTIP2;3 facilitate NH3 transport into the vacuole. Plant Physiol. 2005, 137: 671-680.
Parthier B: The role of phytohormones (cytokinin) in chloroplasts development. Biochem Physiol Pfl. 1979, 174: 173-214.
Race HL, Herrmann RG, Martin W: Why have organelles retained genomes?. Trends Genet. 1999, 15: 364-370.
Che P, Lall S, Nettleton D, Howell SH: Gene expression programs during shoot, root, and callus development in Arabidopsis tissue culture. Plant Physiol. 2006, 141: 620-637.
Cary A, Uttamchandani SJ, Smets R, Van Onckelen HA, Howell SH: Arabidopsis mutants with increased organ regeneration in tissue culture are more competent to respond to hormonal signals. Planta. 2001, 213: 700-707.
Cary AJ, Che P, Howell SH: Developmental events and shoot apical meristem gene expression patterns during shoot development in Arabidopsis thaliana. Plant J. 2002, 32: 867-877.
Atta R, Laurens L, Boucheron-Dubuisson E, Guivarc’h A, Carnero E, Giraudat-Pautot V, Rech P, Chriqui D: Pluripotency of Arabidopsis xylem pericycle underlies shoot regeneration from root and hypocotyl explants grown in vitro. Plant J. 2009, 57: 626-644.
Müller B, Sheen J: Cytokinin and auxin interaction in root stem-cell specification during early embryogenesis. Nature. 2008, 453: 1094-1097.
Werner T, Köllmer I, Bartrina I, Holst K, Schmülling T: New insights into the biology of cytokinin degradation. Plant Biol. 2006, 8: 371-381.
Taniguchi M, Kiba T, Sakakibara H, Ueguchi C, Mizuno T, Sugiyama T: Expression of Arabidopsis response regulator homologs is induced by cytokinins and nitrate. FEBS Lett. 1998, 429: 259-262.
Bürkle L, Cedzich A, Döpke C, Stransky H, Okumoto S, Gillissen B, Kühn C, Frommer WB: Transport of cytokinins mediated by purine transporters of the PUP family expressed in phloem, hydathodes, and pollen of Arabidopsis. Plant J. 2003, 34: 13-26.
Higuchi M, Pischke MS, Mähönen AP, Miyawaki K, Hashimoto Y, Seki M, Kobayashi M, Shinozaki K, Kato T, Tabata S, et al: In planta functions of the Arabidopsis cytokinin receptor family. Proc Natl Acad Sci U S A. 2004, 101: 8821-8826.
Ueguchi C, Koizumi H, Suzuki T, Mituno T: Novel family of sensor histidine kinase genes in Arabidopsis thaliana. Plant Cell Physiol. 2001, 42: 231-235.
de León BG-P, Zorrilla JMF, Rubio V, Dahiya P, Paz-Ares J, Leyva A: Interallelic complementation at the Arabidopsis CRE1 locus uncovers independent pathways for the proliferation of vascular initials and canonical cytokinin signalling. Plant J. 2004, 38: 70-79.
Takei K, Yamaya T, Sakakibara H: Arabidopsis CYP735A1 and CYP735A2 encode cytokinin hydroxylases that catalyze the biosynthesis of trans-zeatin. J Biol Chem. 2004, 279: 41866-41872.
Hirose N, Takei K, Kuroha T, Kamada-Nobusada T, Hayashi H, Sakakibara H: Regulation of cytokinin biosynthesis, compartmentalization and translocation. J Exp Bot. 2008, 59: 75-83.
Heyl A, Riefler M, Romanov GA, Schmülling T: Properties, functions and evolution of cytokinin receptors. Eur J Cell Biol. 2011, 10.1016/j.ejcb.2011.1002.1009. in press
Romanov GA, Lomin SN, Schmülling T: Biochemical characteristics and ligand-binding properties of Arabidopsis cytokinin receptor AHK3 compared to CRE1/AHK4 as revealed by a direct binding assay. J Exp Bot. 2006, 57: 4051-4058.
Adachi S, Nobusawa T, Umeda M: Quantitative and cell type-specific transcriptional regulation of A-type cyclin-dependent kinase in Arabidopsis thaliana. Dev Biol. 2009, 329: 306-314.
Hemerly AS, Ferreira P, de Almeida Engler J, Van Montagu M, Engler G, Inze D: cdc2a expression in Arabidopsis is linked with competence for cell division. Plant Cell. 1993, 5: 1711-1723.
Gaamouche T, Manes C-LO, Kwiatkowska D, Berckmans B, Koumproglou R, Maes S, Beeckman T, Vernoux T, Doonan JH, Traas J, et al: Cyclin-dependent kinase activity maintains the shoot apical meristem cells in an undifferentiated state. Plant J. 2010, 64: 26-37.
Ishiguro S, Watanabe Y, Ito N, Nonaka H, Takeda N, Sakai T, Kanaya H, Okada K: SHEPHERD is the Arabidopsis GRP94 responsible for the formation of functional CLAVATA proteins. EMBO J. 2002, 21: 898-908.
Barton MK: Twenty years on: The inner workings of the shoot apical meristem, a developmental dynamo. Dev Biol. 2010, 341: 95-113.
Rashotte A, Goertzen L: The CRF domain defines cytokinin response factor proteins in plants. BMC Plant Biol. 2010, 10: 74-
Byrne ME: A role for the ribosome in development. Trends Plant Sci. 2009, 14: 512-519.
Horiguchi G, Mollá-Morales A, Pérez-Pérez JM, Kojima K, Robles P, Ponce MR, Micol JL, Tsukaya H: Differential contributions of ribosomal protein genes to Arabidopsis thaliana leaf development. Plant J. 2011, 65: 724-736.
Weijers D, Franke-van Dijk M, Vencken R-J, Quint A, Hooykaas P, Offringa R: An Arabidopsis Minute-like phenotype caused by a semi-dominant mutation in a RIBOSOMAL PROTEIN S5 gene. Development. 2001, 128: 4289-4299.
Swarbreck D, Wilks C, Lamesch P, Berardini TZ, Garcia-Hernandez M, Foerster H, Li D, Meyer T, Muller R, Ploetz L, et al: The Arabidopsis Information Resource (TAIR): gene structure and function annotation. Nucleic Acids Res. 2008, 36: D1009-D1014.
Gene Ontology Consortium.: Creating the gene ontology resource: design and implementation. Genome Res. 2001, 11: 1425-
Prescott AG, Stamford NPJ, Wheeler G, Firmin JL: In vitro properties of a recombinant flavonol synthase from Arabidopsis thaliana. Phytochemistry. 2002, 60: 589-593.
Deikman J, Hammer PE: Induction of anthocyanin accumulation by cytokinins in Arabidopsis thaliana. Plant Physiol. 1995, 108: 47-57.
Frugier F, Kosuta S, Murray JD, Crespi M, Szczyglowski K: Cytokinin: secret agent of symbiosis. Trends Plant Sci. 2008, 13: 115-120.
Böhner S, Gatz C: Characterisation of novel target promoters for the dexamethasone-inducible/tetracycline-repressible regulator TGV using luciferase and isopentenyl transferase as sensitive reporter genes. Mol Gen Genet. 2001, 264: 860-870.
Borevitz JO, Xia Y, Blount J, Dixon RA, Lamb C: Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell. 2000, 12: 2383-2394.
Tamas IA: Plant hormones and their role in plant growth and development. Hormonal regulation of apical dominance. Edited by: Davies PJ. Dordrecht: Martinus Nijhoff Publishers; 1987: 393-410.
Guillaumie S, Goffner D, Barbier O, Martinant J-P, Pichon M, Barrière Y: Expression of cell wall related genes in basal and ear internodes of silking brown-midrib-3, caffeic acid O-methyltransferase (COMT) down-regulated, and normal maize plants. BMC Plant Biol. 2008, 8: 71-
Bishopp A, Help H, Helariutta Y, Kwang WJ: International Review of Cell and Molecular Biology. Volume 276. Cytokinin signaling during root development. Edited by: Jeon K. San Diego: Academic; 2009: 1-48.
Morikawa T, Mizutani M, Aoki N, Watanabe B, Saga H, Saito S, Oikawa A, Suzuki H, Sakurai N, Shibata D, et al: Cytochrome P450 CYP710A encodes the sterol C-22 desaturase in Arabidopsis and tomato. Plant Cell. 2006, 18: 1008-1022.
Morikawa T, Mizutani M, Ohta D: Cytochrome P450 subfamily CYP710A genes encode sterol C-22 desaturase in plants. Biochem Soc Trans. 2006, 34: 1202-1205.
Vercruyssen L, Gonzalez N, Werner T, Schmülling T, Inze D: Combining enhanced root and shoot growth reveals cross talk between pathways that control plant organ size in Arabidopsis. Plant Physiol. 2011, 155: 1339-1352.
Shuai B, Reynaga PCG, Springer PS: The LATERAL ORGAN BOUNDARIES gene defines a novel, plant-specific gene family. Plant Physiol. 2002, 129: 747-761.
Husbands A, Bell EM, Shuai B, Smith HMS, Springer PS: LATERAL ORGAN BOUNDARIES defines a new family of DNA-binding transcription factors and can interact with specific bHLH proteins. Nucleic Acids Res. 2007, 35: 6663-6671.
Ueno Y, Ishikawa T, Watanabe K, Terakura S, Iwakawa H, Okada K, Machida C, Machida Y: Histone deacetylases and ASYMMETRIC LEAVES2 are involved in the establishment of polarity in leaves of Arabidopsis. Plant Cell. 2007, 19: 445-457.
Bortiri E, Chuck G, Vollbrecht E, Rocheford T, Martienssen R, Hake S: ramosa2 encodes a LATERAL ORGAN BOUNDARY domain protein that determines the fate of stem cells in branch meristems of maize. Plant Cell. 2006, 18: 574-585.
Xu B, Li Z, Zhu Y, Wang H, Ma H, Dong A, Huang H: Arabidopsis genes AS1, AS2, and JAG negatively regulate boundary-specifying genes to promote sepal and petal development. Plant Physiol. 2008, 146: 566-575.
Majer C, Hochholdinger F: Defining the boundaries: structure and function of LOB domain proteins. Trends Plant Sci. 2011, 16: 47-52.
Fujimoto SY, Ohta M, Usui A, Shinshi H, Ohme-Takagi M: Arabidopsis ethylene-responsive element binding factors act as transcriptional activators or repressors of GCC box–mediated gene expression. Plant Cell. 2000, 12: 393-404.
Härtel H, Lokstein H, Dörmann P, Grimm B, Benning C: Changes in the composition of the photosynthetic apparatus in the galactolipid-deficient dgd1 mutant of Arabidopsis thaliana. Plant Physiol. 1997, 115: 1175-1184.
Reifarth F, Christen G, Seeliger AG, Dörmann P, Benning C, Renger G: Modification of the water oxidizing complex in leaves of the dgd1 mutant of Arabidopsis thaliana deficient in the galactolipid digalactosyldiacylglycerol. Biochemistry. 1997, 36: 11769-11776.
Block MA, Dorne AJ, Joyard J, Douce R: Preparation and characterization of membrane fractions enriched in outer and inner envelope membranes from spinach chloroplasts. II. Biochemical characterization. J Biol Chem. 1983, 258: 13281-13286.
Yamaryo Y, Kanai D, Awai K, Shimojima M, Masuda T, Shimada H, Takamiya K-, Ohta H: Light and cytokinin play a co-operative role in MGDG synthesis in greening cucumber cotyledons. Plant Cell Physiol. 2003, 44: 844-855.
Christensen SK, Dagenais N, Chory J, Weigel D: Regulation of auxin response by the protein kinase PINOID. Cell. 2000, 100: 469-478.
Friml J, Yang X, Michniewicz M, Weijers D, Quint A, Tietz O, Benjamins R, Ouwerkerk PBF, Ljung K, Sandberg G, et al: A PINOID-dependent binary switch in apical-basal PIN polar targeting directs auxin efflux. Science. 2004, 306: 862-865.
Huang F, Kemel Zago M, Abas L, van Marion A, Galván-Ampudia CS, Offringa R: Phosphorylation of conserved PIN motifs directs Arabidopsis PIN1 polarity and auxin transport. Plant Cell. 2010, 22: 1129-1142.
Kleine-Vehn J, Huang F, Naramoto S, Zhang J, Michniewicz M, Offringa R, Friml J: PIN auxin efflux carrier polarity is regulated by PINOID kinase-mediated recruitment into GNOM-independent trafficking in Arabidopsis. Plant Cell. 2009, 21: 3839-3849.
Michniewicz M, Zago MK, Abas L, Weijers D, Schweighofer A, Meskiene I, Heisler MG, Ohno C, Zhang J, Huang F, et al: Antagonistic regulation of PIN phosphorylation by PP2A and PINOID directs auxin flux. Cell. 2007, 130: 1044-1056.
Zhang J, Nodzyński T, Pĕnčík A, Rolčík J, Friml J: PIN phosphorylation is sufficient to mediate PIN polarity and direct auxin transport. Proc Natl Acad Sci U S A. 2010, 107: 918-922.
Růžička K, Šimášková M, Duclercq J, Petrášek J, Zažímalová E, Simon S, Friml J, Van Montagu MCE, Benková E: Cytokinin regulates root meristem activity via modulation of the polar auxin transport. Proc Natl Acad Sci U S A. 2009, 106: 4284-4289.
Scacchi E, Salinas P, Gujas B, Santuari L, Krogan N, Ragni L, Berleth T, Hardtke CS: Spatio-temporal sequence of cross-regulatory events in root meristem growth. Proc Natl Acad Sci U S A. 2010, 107: 22734-22739.
Ikeda M, Mitsuda N, Ohme-Takagi M: Arabidopsis WUSCHEL is a bifunctional transcription factor that acts as a repressor in stem cell regulation and as an activator in floral patterning. The Plant Cell. 2009, 21: 3493-3505.
Bak S, Feyereisen R: The involvement of two P450 enzymes, CYP83B1 and CYP83A1, in auxin homeostasis and glucosinolate biosynthesis. Plant Physiol. 2001, 127: 108-118.
Hemm MR, Ruegger MO, Chapple C: The Arabidopsis ref2 mutant is defective in the gene encoding CYP83A1 and shows both phenylpropanoid and glucosinolate phenotypes. Plant Cell. 2003, 15: 179-194.
Naur P, Petersen BL, Mikkelsen MD, Bak S, Rasmussen H, Olsen CE, Halkier BA: CYP83A1 and CYP83B1, two nonredundant cytochrome P450 enzymes metabolizing oximes in the biosynthesis of glucosinolates in Arabidopsis. Plant Physiol. 2003, 133: 63-72.
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.
Kiba T, Kudo T, Kojima M, Sakakibara H: Hormonal control of nitrogen acquisition: roles of auxin, abscisic acid, and cytokinin. J Exp Bot. 2011, 62: 1399-1409.
Rus A, Yokoi S, Sharkhuu A, Reddy M, B-h L, Matsumoto TK, Koiwa H, Zhu J-K, Bressan RA, Hasegawa PM: AtHKT1 is a salt tolerance determinant that controls Na+ entry into plant roots. Proc Natl Acad Sci U S A. 2001, 98: 14150-14155.
Mäser P, Eckelman B, Vaidyanathan R, Horie T, Fairbairn DJ, Kubo M, Yamagami M, Yamaguchi K, Nishimura M, Uozumi N, et al: Altered shoot/root Na+ distribution and bifurcating salt sensitivity in Arabidopsis by genetic disruption of the Na+ transporter AtHKT1. FEBS Lett. 2002, 531: 157-161.
Jha D, Shirley N, Tester M, Roy SJ: Variation in salinity tolerance and shoot sodium accumulation in Arabidopsis ecotypes linked to differences in the natural expression levels of transporters involved in sodium transport. Plant Cell Env. 2010, 33: 793-804.
Koiwa H, Bressan RA, Hasegawa PM: Identification of plant stress-responsive determinants in Arabidopsis by large-scale forward genetic screens. J Exp Bot. 2006, 57: 1119-1128.
Møller IS, Gilliham M, Jha D, Mayo GM, Roy SJ, Coates JC, Haseloff J, Tester M: Shoot Na+ exclusion and increased salinity tolerance engineered by cell type-specific alteration of Na+ transport in Arabidopsis. Plant Cell. 2009, 21: 2163-2178.
Rus A, Baxter I, Muthukumar B, Gustin J, Lahner B, Yakubova E, Salt DE: Natural variants of AtHKT1 enhance Na+ accumulation in two wild populations of Arabidopsis. PLoS Genet. 2006, 2: e210-
Mason MG, Jha D, Salt DE, Tester M, Hill K, Kieber JJ, Eric Schaller G: Type-B response regulators ARR1 and ARR12 regulate expression of AtHKT1;1 and accumulation of sodium in Arabidopsis shoots. Plant J. 2010, 64: 753-763.
Nishiyama R, Watanabe Y, Fujita Y, Le DT, Kojima M, Werner T, Vankova R, Yamaguchi-Shinozaki K, Shinozaki K, Kakimoto T, et al: Analysis of cytokinin mutants and regulation of cytokinin metabolic genes reveals important regulatory roles of cytokinins in drought, salt and abscisic acid responses, and abscisic acid biosynthesis. The Plant Cell. 2011, 23: 2169-2183.
Tran L-SP, Urao T, Qin F, Maruyama K, Kakimoto T, Shinozaki K, Yamaguchi-Shinozaki K: Functional analysis of AHK1/ATHK1 and cytokinin receptor histidine kinases in response to abscisic acid, drought, and salt stress in Arabidopsis. Proc Natl Acad Sci U S A. 2007, 104: 20623-20628.
Grotz N, Fox T, Connolly E, Park W, Guerinot ML, Eide D: Identification of a family of zinc transporter genes from Arabidopsis that respond to zinc deficiency. Proc Natl Acad Sci U S A. 1998, 95: 7220-7224.
Mehrtens F, Kranz H, Bednarek P, Weisshaar B: The Arabidopsis transcription factor MYB12 is a flavonol-specific regulator of phenylpropanoid biosynthesis. Plant Physiol. 2005, 138: 1083-1096.
Stracke R, Ishihara H, Huep G, Barsch A, Mehrtens F, Niehaus K, Weisshaar B: Differential regulation of closely related R2R3-MYB transcription factors controls flavonol accumulation in different parts of the Arabidopsis thaliana seedling. Plant J. 2007, 50: 660-677.
Lewis DR, Ramirez MV, Miller ND, Vallabhaneni P, Ray WK, Helm RF, Winkel BSJ, Muday GK: Auxin and ethylene induce flavonol accumulation through distinct transcriptional networks. Plant Physiol. 2011, 156: 144-164.
Chen R, Hilson P, Sedbrook J, Rosen E, Caspar T, Masson PH: The Arabidopsis thaliana AGRAVITROPIC 1 gene encodes a component of the polar-auxin-transport efflux carrier. Proc Natl Acad Sci U S A. 1998, 95: 15112-15117.
Luschnig C, Gaxiola RA, Grisafi P, Fink GR: EIR1, a root-specific protein involved in auxin transport, is required for gravitropism in Arabidopsis thaliana. Genes Dev. 1998, 12: 2175-2187.
Müller A, Guan C, Gälweiler L, Tänzler P, Huijser P, Marchant A, Parry G, Bennett M, Wisman E, Palme K: AtPIN2 defines a locus of Arabidopsis for root gravitropism control. EMBO J. 1998, 17: 6903-6911.
Utsuno K, Shikanai T, Yamada Y, Hashimoto T: AGR, an Agravitropic locus of Arabidopsis thaliana, encodes a novel membrane-protein family member. Plant Cell Physiol. 1998, 39: 1111-1118.
Dello Ioio R, Nakamura K, Moubayidin L, Perilli S, Taniguchi M, Morita MT, Aoyama T, Costantino P, Sabatini S: A genetic framework for the control of cell division and differentiation in the root meristem. Science. 2008, 322: 1380-1384.
Aloni R, Aloni E, Langhans M, Ullrich CI: Role of cytokinin and auxin in shaping root architecture: regulating vascular differentiation, lateral root initiation, root apical dominance and root gravitropism. Ann Bot. 2006, 97: 883-893.
del Pozo JC, Estelle M: The Arabidopsis cullin At CUL1 is modified by the ubiquitin-related protein RUB1. Proc Natl Acad Sci U S A. 1999, 96: 15342-15347.
del Pozo JC, Timpte C, Tan S, Callis J, Estelle M: The ubiquitin-related protein RUB1 and auxin response in Arabidopsis. Science. 1998, 280: 1760-1763.
Leyser HMO, Lincoln CA, Timpte C, Lammer D, Turner J, Estelle M: Arabidopsis auxin-resistance gene AXR1 encodes a protein related to ubiquitin-activating enzyme E1. Nature. 1993, 364: 161-164.
Timpte C, Lincoln C, Pickett FB, Turner J, Estelle M: The AXR1 and AUX1 genes of Arabidopsisfunction in separate auxin-response pathways. Plant J. 1995, 8: 561-569.
Estelle M, Klee HJ: Arabidopsis. Auxin and Cytokinin. New York: Cold Spring Harbor Laboratory Press; 1994,
Ruegger M, Dewey E, Gray WM, Hobbie L, Turner J, Estelle M: The TIR1 protein of Arabidopsis functions in auxin response and is related to human SKP2 and yeast Grr1p. Genes Dev. 1998, 12: 198-207.
Schwechheimer C, Serino G, Deng X-W: Multiple ubiquitin ligase-mediated processes require COP9 signalosome and AXR1 function. Plant Cell. 2002, 14: 2553-2563.
Tiryaki I, Staswick PE: An Arabidopsis mutant defective in jasmonate response is allelic to the auxin-signaling mutant axr1. Plant Physiol. 2002, 130: 887-894.
Xu L, Liu F, Lechner E, Genschik P, Crosby WL, Ma H, Peng W, Huang D, Xie D: The SCFCOI1 ubiquitin-ligase complexes are required for jasmonate response in Arabidopsis. Plant Cell. 2002, 14: 1919-1935.
Nakamura A, Shimada Y, Goda H, Fujiwara MT, Asami T, Yoshida S: AXR1 is involved in BR-mediated elongation and SAUR-AC1 gene expression in Arabidopsis. FEBS Lett. 2003, 553: 28-32.
Ni J, Clark SE: Evidence for functional conservation, sufficiency, and proteolytic processing of the CLAVATA3 CLE domain. Plant Physiol. 2006, 140: 726-733.
Wang G, Fiers M: CLE peptide signaling during plant development. Protoplasma. 2010, 240: 33-43.
Casamitjana-Martínez E, Hofhuis HF, Xu J, Liu C-M, Heidstra R, Scheres B: Root-specific CLE19 overexpression and the sol1/2 suppressors implicate a CLV-like pathway in the control of Arabidopsis root meristem maintenance. Curr Biol. 2003, 13: 1435-1441.
Fiers M, Golemiec E, Xu J, van der Geest L, Heidstra R, Stiekema W, Liu C-M: The 14-amino acid CLV3, CLE19, and CLE40 peptides trigger consumption of the root meristem in Arabidopsis through a CLAVATA2-dependent pathway. Plant Cell. 2005, 17: 2542-2553.
Fiers M, Hause G, Boutilier K, Casamitjana-Martínez E, Weijers D, Offringa R, van der Geest L, van Lookeren Campagne M, Liu C-M: Mis-expression of the CLV3/ESR-like gene CLE19 in Arabidopsis leads to a consumption of root meristem. Gene. 2004, 327: 37-49.
Hobe M, Müller R, Grünewald M, Brand U, Simon R: Loss of CLE40, a protein functionally equivalent to the stem cell restricting signal CLV3, enhances root waving in Arabidopsis. Dev Genes Evol. 2003, 213: 371-
Strabala TJ, O'Donnell PJ, Smit A-M, Ampomah-Dwamena C, Martin EJ, Netzler N, Nieuwenhuizen NJ, Quinn BD, Foote HCC, Hudson KR: Gain-of-function phenotypes of many CLAVATA3/ESR genes, including four new family members, correlate with tandem variations in the conserved CLAVATA3/ESR domain. Plant Physiol. 2006, 140: 1331-1344.
Stahl Y, Wink RH, Ingram GC, Simon R: A signaling module controlling the stem cell niche in Arabidopsis root meristems. Curr Biol. 2009, 19: 909-914.
Wang H, Deng XW: Dissecting the phytochrome A-dependent signaling network in higher plants. Trends Plant Sci. 2003, 8: 172-178.
Saijo Y, Zhu D, Li J, Rubio V, Zhou Z, Shen Y, Hoecker U, Wang H, Deng XW: Arabidopsis COP1/SPA1 complex and FHY1/FHY3 associate with distinct phosphorylated forms of phytochrome A in balancing light signaling. Mol Cell. 2008, 31: 607-613.
Klahre U, Noguchi T, Fujioka S, Takatsuto S, Yokota T, Nomura T, Yoshida S, Chua N-H: The Arabidopsis DIMINUTO/DWARF1 gene encodes a protein involved in steroid synthesis. Plant Cell. 1998, 10: 1677-1690.
Nemhauser JL, Hong F, Chory J: Different plant hormones regulate similar processes through largely nonoverlapping transcriptional responses. Cell. 2006, 126: 467-475.
Zhao Z, Andersen SU, Ljung K, Doležal K, Miotk A, Schultheiss SJ, Lohmann JU: Hormonal control of the shoot stem-cell niche. Nature. 2010, 465: 1089-1092.
Petricka JJ, Benfey PN: Reconstructing regulatory network transitions. Trends Cell Biol. 2011, 21: 442-451.
Clough SJ, Bent AF: Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 1998, 16: 735-743.
Hilson P, Allemeersch J, Altmann T, Aubourg S, Avon A, Beynon J, Bhalerao RP, Bitton F, Caboche M, Cannoot B, et al: Versatile gene-specific sequence tags for Arabidopsis functional genomics: transcript profiling and reverse genetics applications. Genome Res. 2004, 14: 2176-2189.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-
R Development Core Team: R: A language and environment for statistical computing. Vienna, Austria: R Foundation for statistical computing; 2006. ISBN 3-900051-07-0, http://www.R-project.org
Hilson P, Small I, Kuiper MT: European consortia building integrated resources for Arabidopsis functional genomics. Curr Opin Plant Biol. 2003, 6: 426-429.
Allemeersch J: Statistical analysis of microarray data: applications in platform comparison, compendium data, and array CGH. Katholieke Universiteit Leuven, Electrical engineering: Doctoraat; 2006.
Kerr MK, Martin M, Churchill GA: Analysis of variance for gene expression microarray data. J Comput Biol. 2000, 7: 819-837.
Wolfinger RD, Gibson G, Wolfinger ED, Bennett L, Hamadeh H, Bushel P, Afshari C, Paules RS: Assessing gene significance from cDNA microarray expression data via mixed models. J Comput Biol. 2001, 8: 625-637.
Benjamini Y, Hochberg Y: Controlling the false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc B Methods. 1995, 9: 289-300.
Thimm O, Blaesing O, Gibon Y, Nagel A, Meyer S, Krueger P, Selbig J, Mueller LA, Rhee SY, Stitt M: MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004, 37: 914-939.
Usadel B, Nagel A, Thimm O, Redestig H, Blaesing OE, Palacios-Rojas N, Selbig J, Hannemann J, Piques MC, Steinhauser D, et al: Extension of the visualization tool MapMan to allow statistical analysis of arrays, display of coresponding genes, and comparison with known responses. Plant Physiol. 2005, 138: 1195-1204.
Pavlidis P, Noble WS: Analysis of strain and regional variation in gene expression in mouse brain. Genome Biol. 2001, 2: research0042–-research0042.15.
Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S, Hub A, Group WPW: AmiGO: online access to ontology and annotation data. Bioinformatics. 2009, 25: 288-289.
Thareau V, Dehais P, Serizet C, Hilson P, Rouze P, Aubourg S: Automatic design of gene-specific sequence tags for genome-wide functional studies. Bioinformatics. 2003, 19: 2191-2198.
Rozen S, Skaletsky HJ: Bioinformatics Methods and Protocols: Methods in Molecular Biology. Primer3 on the WWW for general users and for biologist programmers. Edited by: Krawetz S, Misener S. Totowa, NJ: Humana Press; 2000: 365-386.
We thank Joke Allemeersch for help and advice in bioinformatics and Camilla Liput and Petra-Michaela Hartmann for performing transcript analyses by real-time PCR. We thank Stefanie Zintl for discussing the links between cytokinin and light. This work was funded by the DFG in the frame of the Arabidopsis Functional Genomics Network (AFGN) and the CAGE project (http://www.cagecompendium.org, EU project No. QLK3-CT-2002-02035).
WGB carried out all the experimental and bioinformatical work, designed the research and bioinformatic approaches and contributed to the analysis of the data and writing the paper. Both read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1. Flow chart representing the algorithm to categorize genes detected in this microarray study. Schematic representation of the algorithm used for gene categorization. Only genes, which were detected on at least 25% of the microarrays were fed into the algorithm. (PDF 86 KB)
Additional file 2: Table S1. Gene lists from cytokinin-treated samples. Lists of all genes detected on at least 25% of the microarrays (i.e. 8 arrays) containing expression data of wild-type samples sorted into the different categories as described in Methods and Additional file 1: Figure S1. The spreadsheet is divided into sections a-f containing genes which are root-specifically (a), shoot-specifically (b), differentially (c), or similarly (d) regulated by cytokinin. Section e contains cytokinin-regulated genes which are not unambiguously classifiable into one of the other categories. Section f contains genes not regulated by cytokinin. CATMA ID, unique identifier of the probe on the CATMA microarray; FDR p-value, p-value corrected for multiple testing; AGI, unique gene identifier assigned to Arabidopsis thaliana genes by TAIR. BA0, mock treatment; BA15, 15 min cytokinin treatment; BA120, 2 h cytokinin treatment; BA1080, 18 h cytokinin treatment. The sorting criteria were (a) the sum of all ratios in the root; (b) the sum of all ratios in the shoot; (c) the sum of the absolute value of the log2 of all ratios in root and shoot; (d) the sum of all ratios in root and shoot; (e) the same as for (c); (f) CATMA ID. (XLS 7 MB)
Additional file 3: Table S2. Gene lists from cytokinin-deficient samples. Lists of all genes detected on at least 25% of the microarrays (i.e. 4 arrays) containing expression data of mock-treated wild-type samples (BA0) and samples of cytokinin-deficient plants (CKX1), sorted into the different categories as described in Methods and Additional file 1: Figure S1. The spreadsheet is divided into sections a-f containing genes which are root-specifically (a), shoot-specifically (b), differentially (c), or similarly (d) regulated. Section e contains cytokinin-regulated genes which are not unambiguously classifiable into one of the other categories. Section f contains genes not found regulated by cytokinin. CATMA ID, unique identifier of the probe on the CATMA microarray; FDR p-value, p-value corrected for multiple testing; AGI, unique gene identifier assigned to Arabidopsis thaliana genes by TAIR. The sorting criteria were (a) the ratio CKX1 vs. Col-0 root; (b) the ratio CKX1 vs. Col-0 shoot; (c) the sum of the absolute value of the log2 of both ratios; (d) the sum of both CKX1 vs. Col-0 ratios; (e) the same as for (c); (f) CATMA ID. (XLS 4 MB)
Additional file 4: Table S3. Genes expressed in a more shoot-like fashion in cytokinin-treated roots. (a) List of genes, which are higher expressed in the shoot than in the root and upregulated by 18 h of cytokinin treatment in the root. (b) List of genes, which are lower expressed in the shoot than in the root and downregulated by 18 h of cytokinin treatment in the root. FDR-corrected p-values ≤ 0.03 marked in blue indicate whether a gene is significantly regulated by cytokinin (cytokinin effect), significantly differently expressed roots and shoots (organ effect), or whether its regulation by cytokinin differs significantly between roots and shoots (organ effect). The threshold set for differential expression is 2.5-fold. (XLS 128 KB)
Additional file 5: Figure S2. GO term enrichment for genes indicative of a developmental shift. (a) GO term enrichment of genes which are more shoot-like expressed in cytokinin-induced roots. The figure shows that the only overrepresented class of genes encodes proteins located in the chloroplast. The dataset of Additional file 4: Table 3a was used to generate this figure. Darker boxes indicate a higher significance of overrepresentation of the GO term indicated. (b) GO term enrichment of genes which are more root-like expressed in cytokinin-deficient shoots. The figure shows that the only overrepresented class of genes encodes proteins located in the chloroplast. The dataset of Additional file 4: Table 3b was used to generate this figure using the AmiGO tool (see Methods). (PDF 114 KB)
Additional file 6: Figure S3. Expression of selected organ-specific cytokinin-regulated genes mentioned in the text as revealed by microarray analysis and real-time quantitative RT-PCR. (a) Differentially expressed genes after cytokinin induction. (b) Expression ratios of untreated shoots and roots of the genes shown in (a). (c) Differentially expressed genes under cytokinin deficiency. (d) Expression ratios in wild-type shoots and roots of the genes shown in (c). Cytokinin regulation of selected genes with different organ specificity of the response was tested by microarrays (indicated by diamonds) and qRT-PCR (indicated by bars) using independent biological material. BD, below detection limit of microarray; bd, below detection limit of qRT-PCR. (PDF 1 MB)
Additional file 7: Figure S4. Expression levels of selected organ-specific cytokinin-regulated genes mentioned in the text as revealed by microarray analysis and real-time quantitative RT-PCR. (a, b) Expression levels of differentially expressed genes after cytokinin treatment. (c, d) Expression levels of differentially expressed genes under cytokinin-deficiency. The expression levels of selected genes with different organ specificity of the response were tested under the conditions indicated in the legend by microarray analysis (a, c) and qRT-PCR (b, d) using independent biological material. Open diamonds indicate that less than two of the four spots on the microarray were significantly above background. nd, not detected. (PDF 1 MB)
Additional file 8: Table S4. Genes expressed in a more root-like fashion in cytokinin-deficient shoots. (a) List of genes, which are higher expressed in the root than in the shoot and upregulated by chronic cytokinin deficiency in the shoot. (b) List of genes, which are lower expressed in the root than in the shoot and downregulated by chronic cytokinin deficiency in the shoot. FDR-corrected p-values ≤ 0.03 marked in blue indicate whether a gene is significantly regulated by cytokinin (cytokinin effect), significantly differently expressed roots and shoots (organ effect), or whether its regulation by cytokinin differs significantly between roots and shoots (organ effect). The threshold set for differential expression is 2.5-fold. (XLS 247 KB)
Additional file 9: Table S5. Root- and shoot-specific changes of transcript abundance of known cytokinin-regulated genes of Arabidopsis thaliana. A selection of previously published cytokinin-responsive genes from various publications were analyzed for their transcriptional regulation in response to cytokinin induction and cytokinin deficiency in root and shoot samples. Expression ratios are colour-coded as described for Table 1. Significance codes indicate p-values calculated as described in Methods: 0.1 > · > 0.05 > * 0.01 > ** 0.001 > ***. Genes that were published in at least two previous publications are listed here. In addition, selected examples of known regulated genes that were only published once were added. A short list of genes can be found in Table 2. The publications are coded by letters: A, Rashotte et al., 2003 ; B, Brenner et al., 2005 ; C, Kiba et al., 2005 ; D, Rashotte et al., 2006 ; E, Taniguchi et al., 2007 ; F, Argyros et al., 2008 . AGI, unique gene identifier assigned to Arabidopsis thaliana genes by TAIR. (PDF 336 kb) (PDF 336 KB)
Additional file 10: Table S6. Cytokinin-responsive genes encoding ribosomal proteins. The table shows all cytokinin-regulated genes encoding ribosomal protein according to the GO localization database. Columns and colour-codes are the same as for Table 2. The additional column localization shows whether the protein is part of a plastid or cytosolic ribosome, according to the GO database (http://www.arabidopsis.org/tools/bulk/go/index.jsp). (PDF 231 KB)
About this article
Cite this article
Brenner, W.G., Schmülling, T. Transcript profiling of cytokinin action in Arabidopsis roots and shoots discovers largely similar but also organ-specific responses. BMC Plant Biol 12, 112 (2012). https://doi.org/10.1186/1471-2229-12-112
- Shoot Apical Meristem
- Shoot Sample
- Cytokinin Response
- Cytokinin Treatment
- Cytokinin Response Factor