- Research article
- Open Access
Characterization of the cytokinin-responsive transcriptome in rice
BMC Plant Biology volume 16, Article number: 260 (2016)
Cytokinin activates transcriptional cascades important for development and the responses to biotic and abiotic stresses. Most of what is known regarding cytokinin-regulated gene expression comes from studies of the dicotyledonous plant Arabidopsis thaliana. To expand the understanding of the cytokinin-regulated transcriptome, we employed RNA-Seq to analyze gene expression in response to cytokinin in roots and shoots of the monocotyledonous plant rice.
We identified over 4,600 and approximately 2,400 genes differentially expressed in response to cytokinin in roots and shoots respectively. There were some similarities in the sets of cytokinin-regulated genes identified in rice and Arabidopsis, including an up-regulation of genes that act to reduce cytokinin function. Consistent with this, we found that the preferred DNA-binding motif of a rice type-B response regulator is similar to those from Arabidopsis. Analysis of the genes regulated by cytokinin in rice revealed a large number of transcription factors, receptor-like kinases, and genes involved in protein degradation, as well as genes involved in development and the response to biotic stress. Consistent with the over-representation of genes involved in biotic stress, there is a substantial overlap in the genes regulated by cytokinin and those differentially expressed in response to pathogen infection, suggesting that cytokinin plays an integral role in the transcriptional response to pathogens in rice, including the induction of a large number of WRKY transcription factors.
These results begin to unravel the complex gene regulation after cytokinin perception in a crop of agricultural importance and provide insight into the processes and responses modulated by cytokinin in monocots.
Phytohormones regulate many aspects of plant growth and development as well as responses to the environment in part through an integrated modulation of the transcriptome . Cytokinins are N 6-substituted adenine derivatives that were discovered based on their ability to promote cell division in cultured cells [2, 3]. Since their discovery, cytokinins have been implicated in almost all facets of plant growth and development, including leaf senescence, meristem maintenance, sink/source relationships, and biotic/abiotic interactions [4–6]. Much progress has been made in understanding cytokinin biosynthesis, signaling, and gene regulation, mostly from genetic and molecular studies in the dicot Arabidopsis thaliana, including the identification of the key elements of the cytokinin signaling pathway , which enabled studies of the function of cytokinin in other plants, including its role in nodulation in legumes  and phylotaxy in maize .
The cytokinin signaling pathway is comprised of a ‘two-component’ phosphorelay that culminates in the activation of the type-B response regulators (RRs), transcription factors that modulate downstream gene expression . The first step in the signaling pathway is mediated by the hybrid histidine-kinase receptors (HKs), which are embedded in the endoplasmic reticulum membrane such that the amino-terminal cytokinin-binding domain is localized within the lumen of the endoplasmic reticulum [9, 10]. Upon cytokinin binding, the HKs autophosphorylate on a conserved histidine, and then transfer the phosphate to an Asp residue within their C-terminal receiver domain [11–14]. This phosphate is subsequently transferred to the histidine phosphotransfer proteins (HPs), which in turn transfer the phosphate to the type-A and type-B RRs [15, 16]. The type-B RRs are activated by phosphorylation on a conserved Asp and act as positive elements of signaling. They contain a C-terminal DNA-binding domain and initiate a complex transcriptional cascade to drive the appropriate gene expression changes downstream of cytokinin perception [17–20]. The type-A RRs are cytokinin primary response genes that negatively regulate the signaling pathway [21–24].
Rice (Oryza sativa) is a staple crop for almost half of the world’s population. The global consumption of rice has risen steadily over the past 50 years and is projected to continue to increase as the world’s population rises through 2050 , pointing toward an urgent need to increase production. Cytokinin plays a key role in determining yield from rice; a major quantitative trait loci underlying yield was identified as a cytokinin oxidase gene, which encodes a protein that degrades cytokinin . Phylogenetic analysis has identified rice homologs of Arabidopsis cytokinin signaling pathway components [27–29] and functional studies of a few of these genes have confirmed their role in cytokinin signaling. For example, the OsHK6 gene encodes an HK that binds cytokinin and complements an Arabidopsis cytokinin receptor mutant . Similarly, a rice type-B RR (OsRR22) complements a type-B RR mutant of Arabidopsis . Disruption of expression of OsHP1 and OsHP2 by RNAi results in reduced sensitivity to exogenous cytokinin and phenotypes consistent with reduced cytokinin function . Overexpression of several type-A OsRRs in rice blocks shoot regeneration in response to cytokinin and results in reduced root sensitivity to cytokinin, suggesting that they negatively regulate cytokinin signaling, similar to what is found in Arabidopsis [32, 33]. While these results indicate that, at a basic level, the backbone of the cytokinin signaling pathway likely operates in a similar manner in monocots and dicots, the downstream processes regulated by cytokinin in rice have yet to be established.
Global gene expression responses to cytokinin have been extensively studied in Arabidopsis using microarray and RNA-Seq analyses [34–38]. Many genes that are differentially expressed as early as 15 min after the treatment encode transcription factors, suggesting that cytokinin not only triggers immediate gene expression changes, but also activates complex transcriptional cascades. Here, we use RNA-Seq to identify genes regulated by cytokinin in the roots and shoots of rice seedlings. Identifying these changes in response to exogenous cytokinin defines the distinct patterns of expression in response to cytokinin in the two different tissues. Comparing the differentially expressed genes in rice to a similar experiment in Arabidopsis reveals similarities and differences in the role of cytokinin between these monocot and dicot species. This study begins to unravel the complex gene regulation after cytokinin perception in a crop of agricultural importance and provides insight into the processes and responses modulated by cytokinin in monocots.
Results and discussion
Identification of cytokinin-responsive genes in rice
To investigate cytokinin regulation of gene expression in a monocot, we performed high throughput cDNA sequencing (RNA-Seq) of libraries prepared from rice seedlings treated for two hours with the cytokinin benzyladenine (BA). Rice seedlings were grown hydroponically and cytokinin delivered via addition to the hydroponic media. Twelve libraries were prepared in total, comprised of three replicates each of BA and mock-treated roots and shoots. Libraries were sequenced on the Illumina HiSeq platform, yielding approximately 30 to 50 million single-end, 100-bp reads per library. In each library, at least 90% of reads had a mean Phred score of ≥ 28 and more than 95% could be mapped to a single location in the rice genome. Altogether, more than 447 million reads mapped to a unique genomic location. Thus coverage of the rice transcriptome was deep enough to provide a detailed view of how cytokinin affected gene expression in both roots and shoots of rice seedlings. To facilitate re-use of the data in other studies, we configured the Integrated Genome Browser (IGB)  to offer access to RNA-Seq alignment files, pre-computed coverage graphs, and splice junction files. To view the data, readers should download the browser, select the latest rice genome, and then browse and select data in the Data Access tab.
Read alignments were compared to rice gene models from the Michigan State University rice annotation project’s MSU7 release . Comparing read alignments to annotated genes in MSU7 identified approximately 30,000 genes with 20 mapped reads or more across all samples (Additional file 1: Table S1). Using this as a minimal threshold for calling a gene expressed, we detected expression for 53% of the 55,987 annotated rice genes. For comparisons between genes, expression values were calculated as the number of reads per kilobase of expressed sequence per million mapped reads (RPKM; Additional file 2: Table S2). Other gene model collections are available, such as annotations from the Rice Annotation Project Database , but we chose to use the MSU7 release in part because of the availability of informatics tools needed for functional interpretation of the data, such as GO annotations and Arabidopsis ortholog assignments. In general, we have found that the MSU7 and RGAP-DB annotation collections are congruent in that genes annotated in MSU7 are typically present in the RGAP-DB annotations, and vice versa; tools that map gene names between models are available, and to further facilitate comparisons, we configured IGB to provide both sets of annotations together with the RNA-Seq data.
Many of the same genes were expressed in shoots and roots of rice, but the overall profile of gene expression was different between the two tissues. The similarities and differences between gene expression in roots and shoots is clear when visualized on a chromosome or region-wide scale using IGB. Figure 1a shows an example of RNA-Seq coverage graphs from the six control libraries depicting read density by chromosome position for a portion of chromosome 1. Several of these genes show differential expression in the roots relative to shoot. Such regions are widespread throughout the rice genome. In total, more than 17,000 genes were differentially expressed between the roots and shoots in control conditions (Additional file 3: Table S3), with about 9,000 of them having a higher expression level in the root.
Consistent with RNA-Seq studies in other plant species [37, 42, 43], a relatively small number of highly expressed genes in rice gave rise to most of the sequence reads (Additional file 1: Tables S1 and Additional file 2: Table S2). In each of the 12 libraries, more than 75% of reads mapped to the top 10% of the most highly expressed genes. Overall diversity of gene expression in roots and shoots was similar, but the types of genes that were most highly expressed in the roots versus shoots reflected tissue-specific functions. In shoots, genes involved in photosynthesis were typically the most highly expressed; RuBISCO and chlorophyll a/b binding protein genes were six of the ten most highly expressed genes, exceeding 7,000 RPKM in the case of RuBISCO homolog (LOC_Os12g19381). In roots, three of the top ten most highly expressed genes encoded metallothioneins, cysteine-rich metal chelators associated with heavy metal tolerance , and two of these (LOC_Os03g17870, LOC_Os12g38300) were also highly enriched in roots compared to shoots. The most highly expressed gene in rice roots was metallothionein-encoding gene LOC_Os12g38300, which averaged 5,425 RPKM in roots and only 10 RPKM in shoots, a more than 500-fold difference in relative expression. Other genes that were highly expressed in roots but not shoots included several expressed proteins of no known function, a putative aquaporin, several putative protease inhibitors, and a rapid alkalinization protein family member (RALF). Interestingly, RALF protein genes were also among the most highly expressed genes in Arabidopsis pollen , and yet the function of RALF proteins remains largely unknown.
To identify cytokinin-responsive genes in rice, we performed differential expression analysis of the RNA-Seq data, analyzing root and shoot samples separately. In roots, the BA treatment caused differential expression (DE) of approximately 4,700 genes (FDR ≤ 0.0001), of which approximately 60% were up-regulated (Additional file 4: Table S4). In shoots, approximately 2,400 genes were differentially expressed, of which 64% were up-regulated (Fig. 1b; Additional file 5: Table S5). There were 609 genes that were DE in both roots and shoots (Fig. 1c), and most of these changed in same direction. However, the magnitude of the change was usually larger in roots than shoots Fig. 1d). The larger number of genes induced in the roots and their greater induction could reflect the mode of application of the hormone: BA was added directly to the hydroponic media bathing the roots and thus the added cytokinin reached the shoots primarily via the transpiration stream.
To confirm results from the RNA-Seq analysis, we used the NanoString nCounter assay [45, 46] to re-test the response of a subset of genes identified by RNA-Seq as DE in response to cytokinin using three biological replicates distinct from those used for the RNA-Seq analysis. This subset included genes identified as induced or repressed specifically in roots, induced in both tissues, and oppositely regulated in roots and shoots; the subset included genes involved in hormone function (cytokinin, auxin, GA, ethylene, and brassinosteroid), biotic interactions (including five WRKY transcription factors), transcription factors, and various other genes (Table 1). Of the 32 genes retested as being DE in roots, 30 were confirmed as DE by cytokinin using the NanoStrings assay (Fig. 2a). Of the eight genes retested as being DE in shoots, six were confirmed using the NanoStrings assay (Fig. 2b). Furthermore, the log2 fold changes obtained from NanoStrings and RNA-Seq were similar (Fig. 2). For the comparison of genes identified in roots as DE, the Pearson’s correlation coefficient was 0.93, and for the shoots comparison, the correlation was 0.86. Overall, the NanoStrings assay indicates high confidence in the genes identified as DE in the RNA-Seq dataset. There were 477 genes that were DE in both roots and shoots in the same direction (Additional file 6: Table S6). The list included several genes involved in modulating the response to cytokinin, including eight type-A RRs, four cytokinin oxidases, two cytokinin cytokinin-O-glucosyltransferases, and two genes (OsHK3 and OsHK4) encoding cytokinin receptor histidine kinases (Additional file 6: Table S6 and Additional file 7: Figure S1). The list also included 46 genes annotated as transcription factors  from diverse families, with an enrichment for members of the MYB family and two annotated with functions related to ethylene and gibberellin signaling pathways (the gibberellin-associated factor LOC_Os07g39470 was up-regulated; the putative ethylne-responsive transcription factor2 LOC_Os01g54890 was down-regulated). The DE genes also included genes encoding a putative gibberellin receptor, a gibberellin 2-beta-dioxygenase, an ethylene-responsive protein, two ACC oxidase proteins (involved in ethylene biosynthesis) and eight genes annotated with auxin-related functions, including 3 Aux-IAA genes. Two genes annotated as brassinosteroid-insensitive 1-associated receptor kinases were differentially expressed; one was up-regulated and the other down-regulated. Overall, the modulation of expression of these genes suggests that application of cytokinin triggers changes in other hormone-related pathways, highlighting how plant hormone pathways intersect and influence each other. Other notable genes regulated in both roots and shoots included six peroxidases, nine cytochrome P450s, four aquaporins, and two HKT Na+ and Na+/K+ transporters (OsHKT1;5 and OsHKT2;1).
Most genes that were DE in both roots and shoots changed in the same direction. However, there were 131 genes in which the direction of the change was different in the two tissues (Additional file 8: Table S7). Cytokinin affects root and shoots in different ways , and so it was interesting to investigate these genes in detail. Eleven of the 132 genes encoded transcriptional regulators, including the type-B RR OsRR22 (LOC_Os06g08440), which was down-regulated in roots and up-regulated in shoots. Only one other type-B RR (OsRR26, LOC_Os01g67770) was regulated by cytokinin, being down-regulated in roots and unaffected in shoots. Other genes that were regulated in opposite directions included ACT domain repeat protein 5 (ACR5), likely involved in glutamine-sensing , OsSAUR21, the phosphatidylethanolamine-binding protein RCN1, which is involved in the transition of rice plants to the reproductive phase , and several protein kinases and other genes with annotated functions related to signal transduction, including genes involved in calcium signaling.
We examined if the DE genes from rice roots were similar to those differentially regulated by cytokinin in Arabidopsis. To this end, we examined a list of genes robustly regulated by cytokinin in Arabidopsis across several studies, the so-called “golden list”  and identified rice counterparts of these genes using orthology annotations from the MSU web site that were generated using OrthoMCL with the default parameters [40, 51]. Approximately half of the Arabidopsis genes on the golden list were annotated with one or more orthologs from rice (116/227). Of these, 55% were differentially expressed in the rice data set in the same direction as in Arabidopsis, and 17% were differentially expressed in opposite directions (Additional file 7: Figure S2). The similarly regulated genes included a number known to negatively regulate cytokinin function, including type-A ARRs and cytokinin oxidases. Thus, a substantial portion of the transcriptional response is conserved between monocots and dicots.
OsRR22 binds upstream motifs similar to those bound by Arabidopsis type-B RRs
The type-B response regulators are positive elements of cytokinin signaling and are required to modulate transcription of genes in response to the hormone [17, 18, 20, 52]. In Arabidopsis, a common DNA motif sufficient for binding of ARR1, ARR2, ARR10, and ARR11 has been identified as AGAT(T/C), with an extended ARR1 binding motif (AAGAT[T/C]TT) being elucidated due to its enrichment in promoters of known primary response genes [52–55]. Repeats of the extended motif have been incorporated into cytokinin reporters and successfully used to visualize cytokinin activity in plant tissues . More recently, protein binding microarrays have been used to identify transcription factor binding motifs, including those for a subset of the type-B RRs of Arabidopsis (ARR1, ARR2, ARR11, ARR14, ARR18) [57, 58], an approach we chose to use for characterizing type-B RRs of rice.
To determine binding sites of the rice type-B OsRRs, we expressed the DNA-binding domains of OsRR22, OsRR27 and OsRR29 as fusions to maltose-binding protein in vitro and used the protein binding microarray technology (PBM11) to assay binding on all possible 8-mers followed by computational analysis to identify binding motifs . We characterized OsRR22 because it belongs to the same type-B RR subfamily as the Arabidopsis type-B RRs implicated in regulating the cytokinin response, OsRR22 being most similar to AtARR10 and AtARR12. Functional analysis has shown that OsRR22 is able to complement the mutant phenotype of an Arabidopsis arr1;arr12 loss-of-function line, as well as transactivate AtARR6-LUC in transient expression assays, providing evidence of conserved function across monocot and dicot species . Three related binding motifs were identified for OsRR22 (Fig. 3a), two of these with an AGAT core sequence and the third with a palindromic sequence built around the AT core. Significantly, the same motifs have been identified for Arabidopsis type-B RRs of the same subfamily , indicating a conservation of binding between the rice and Arabidopsis type-B RRs implicated in regulating the cytokinin response. In contrast to OsRR22, OsRR27 and OsRR29 belong to monocot-specific subfamilies of type-B RRs and contain substantially diverged Myb-like DNA-binding domains . No significantly enriched motifs were identified for OsRR27 and OsRR29, although both produced fusion proteins at similar levels to those found with OsRR22 (Additional file 7: Figure S3). OsRR27 and OsRR29 may require additional sequences outside the Myb-like domain, additional cofactors and/or plant-specific post-translational modifications to allow for binding to DNA.
We searched for instances of the OsRR22 binding motifs within the putative regulatory regions (1 kb upstream of the transcriptional start site) of cytokinin-regulated genes, the prediction being that direct targets of the type-B RRs would be enriched for the binding motifs. Of the 2,890 genes that were up-regulated in roots, 27.8% contained one or more of the OsRR22-like extended binding motifs. To determine if this percentage was unusually large or small, we generated 1,000 random samples of 2,890 rice genes and for each sample, calculated the percentage of genes that contained an OsRR22 binding site. Out of 1,000 random samples of 2,890 genes, only 2.2% contained a higher percentage of genes with OsRR22 binding sites than the set of 2,890 genes that were up-regulated in roots (Fig. 3b). By contrast, of the 1,780 genes that were down-regulated in roots, only 23.7% contained an OsRR22 binding site (Additional file 7: Fig. 3c). To assess the significance of this, we repeated the random sampling experiment, selecting 1,000 random samples of 1,780 genes and for each sample, calculating the percentage of genes with one or more OsRR22 sites. Of the 1,000 random samples, only 1% contained a lower percentage of genes with an OsRR22 binding site than the set of 1,780 genes that were down-regulated by cytokinin. We observed similar, but less significant trends for differentially expressed genes in shoots. Thus, OsRR22 binding sites were unusually prevalent among up-regulated but not down-regulated genes, consistent with OsRR22 being primarily a positive regulator of cytokinin signaling. However, many DE genes (more than 70%) lacked an extended OsRR22 binding motif. These genes may contain only the core type-B RR binding motif (which, although sufficient for binding, is too short to be diagnostic for enrichment), a type-B RR binding site diverged from this consensus, a type-B RR binding site outside of the 1-kb window examined, or may be dependent on other factors besides the type-B RRs for regulation by cytokinin.
Gene Ontology enrichment analysis of cytokinin-regulated genes in rice
To explore the biological processes in which cytokinin participates in rice, we conducted a Gene Ontology enrichment analysis on genes DE in response to cytokinin. We first examined genes induced or repressed by cytokinin in both root and shoot tissues (Additional file 7: Figures S4 and S5). Among the induced genes, as expected, we saw an enrichment of genes involved in two-component signaling and signal transduction (Additional file 7: Figures S1 and S4). Additionally, we observed a significant enrichment in genes involved in adenine and purine salvage, processes that have been implicated in the interconversion of cytokinin (CK) bases, ribosides, and nucleotides . Similar to GO term enrichment studies in Arabidopsis treated tissues [34, 37, 38, 61], genes involved in secondary metabolic processes were enriched among the genes induced in both roots and shoots. Among the genes repressed in both roots and shoots, the most significant enrichment was for GO terms related to the response to oxidative stress (Additional file 7: Figure S5).
We next examined GO term enrichment among the genes DE in the root tissue. Consistent with cytokinin regulating a transcriptional cascade, there is an enrichment of genes with a role in the regulation of transcription among the genes induced by cytokinin in the root (Additional file 7: Figure S6), with 550 genes encoding transcription factors showing DE in response to cytokinin out of a total of 1,611 transcription factors encoded by the rice genome . Thus, remarkably, more than one-third of rice transcription factors are differentially expressed in response to cytokinin in rice roots. This suggests large scale changes in the rice root transcriptome following cytokinin treatment, consistent with the large number of genes we observe as differentially expressed in response to cytokinin.
There is a strong enrichment for genes involved in post-translational modification of proteins among the genes induced by cytokinin in the root, specifically in ubiquitination and phosphorylation (Additional file 7: Figure S6), the former perhaps indicating a role for cytokinin in regulating protein turnover in rice. Several studies have linked cytokinin to the regulation of protein turnover in Arabidopsis, although such a strong GO enrichment as we find in rice has not been previously noted in other cytokinin transcriptome studies. In Arabidopsis, the RPN12 subunit of the 26S proteasome plays a role in regulating cytokinin sensitivity . In addition, the cytokinin signaling pathway itself is subject to regulation by protein turnover. The turnover of a subset of type-A ARR proteins occurs through a Rub protein modification-dependent pathway  and is regulated by cytokinin via differential phosphorylation of the conserved Asp residue . Cytokinin transiently stabilizes degradation of ARR1 in Arabidopsis , a type-B RR, whose turnover is also regulated by the KISS ME DEADLY family of F-box-proteins . In contrast, cytokinin is proposed to promote the 26S proteasome-dependent turnover of ARR2 via regulation of phosphorylation of the conserved Asp target . Cytokinin also regulates the turnover of other proteins. For example, cytokinin regulates the turnover of ACC synthase, which catalyzes the key step involved in ethylene biosynthesis [68, 69]. Further, during lateral root development in Arabidopsis, cytokinin targets dephosphorylated PINs present preferentially on anticlinal membranes for lytic degradation . Together, these studies from Arabidopsis indicate that cytokinin regulates the turnover of its own signaling pathway, as well as additional downstream targets. The large number of genes involved in protein degradation regulated by cytokinin in rice suggests that the modulation of protein turnover may be a prominent mechanism by which cytokinin regulates cell function in rice.
One intriguing category of genes enriched among those induced by cytokinin in the root encode proteins involved in the synthesis of siderophores (Additional file 7: Figure S6). Siderophores are small, high-affinity iron chelating compounds secreted by graminaceous plants to increase the uptake of these minerals [71–73]. In rice these compounds are secreted at lower levels than in other grasses and may be more involved in the uptake of zinc rather than iron . While cytokinin has been linked to the uptake of other minerals, including nitrate, phosphate, sulfur, and iron , it has not previously been linked to zinc uptake.
Finally, we examined GO enrichment among the genes DE by cytokinin in the shoot. In addition to the previously mentioned two-component signaling and adenine and purine salvage, we observed an enrichment of genes induced by cytokinin that are involved in the regulation of cell growth, lipid metabolism, and secondary metabolism (Additional file 7: Figure S7). The GO term secondary metabolic processes is also enriched in cytokinin-treated Arabidopsis tissues [37, 38, 61]. Genes involved in the synthesis of siderophores are enriched among the down-regulated genes (Additional file 7: Figure S8), which is distinct from the induction of these genes by cytokinin in roots (see above). Other enriched terms include flower development and protein phosphorylation.
To further investigate the functions of cytokinin-regulated genes, we analyzed differentially expressed genes in rice roots (Fig. 4a-c) and shoots (Additional file 7: Figure S9) using functional annotations from MapMan . In general, the enrichment analysis using MapMan categories correlated well with results observed in GO analysis. As expected, there were a large number of genes linked to development differentially regulated by cytokinin in both roots and shoots (Figs. 4a and Additional file 7: Figure S8). Many genes involved in hormone signaling and metabolism were regulated in both the root and shoot, most markedly those related to ethylene (Figs. 4a and Additional file 7: Figure S8), the majority of which were induced by cytokinin. There was also a consistent up-regulation of genes involved in jasmonate function, and general down-regulation of genes involved in gibberellin function. Cytokinin also affected many genes involved in auxin, ABA, and salicylic acid function in the root.
As observed with the GO enrichment analysis described above, there were a large number of genes involved in the regulation of transcription and genes linked to protein degradation regulated by cytokinin in rice roots (Fig. 4a). Interestingly, many kinases, especially receptor-like kinases, were up-regulated in roots, which is consistent with the over-representation of the GO term protein phosphorylation noted above and which suggests that cytokinin may regulate cellular function through modification of Ser/Thr phosphorylation of many downstream targets. Consistent with a role as a negative regulator of cell proliferation in roots, multiple genes involved in cell division and the cell cycle were down-regulated by cytokinin (Fig. 4b). MapMan analysis also identified an unusually large percentage of cytokinin-regulated genes related to stress-response functions (Figs. 4b and 5). These genes included those involved in the response to biotic stimuli and abiotic stimuli (heat, cold, and drought/salt) with the majority found to be up-regulated by cytokinin. There was also an enrichment of genes involved in calcium and G-protein signaling (Fig. 4a), suggesting that cytokinin may crosstalk with these key signaling pathways. Other enriched categories of genes included, peroxidases and GSTs, and genes involved in cell wall function (Fig. 4c).
Using statistical analysis of MapMan annotations, we identified MapMan categories that contained an unusually large number of differentially expressed genes (Fig. 5 and Additional file 9: Tables S8 and Additional file 10: Table S9). Members of the WRKY transcription factor family contained an unusually large number of DE genes in roots (Fig. 5), consistent with a modulation of a genes involved in biotic stress response. Several categories of kinases were overrepresented, including LRKs, S-locus receptor kinases. Other notable categories of genes enriched in cytokinin-induced genes in roots were cytochrome P450s and genes involved in calcium signaling. Among the down-regulated genes in the roots, genes involved in cell wall function, chromatin modification, and peroxidases were enriched. In the shoots, genes involved in protein function (folding, initiation of protein synthesis and elongation) and the ribosome were overrepresented.
Cytokinin-regulated genes overlap with genes driving defense responses
Our analysis suggests a substantial overlap in the genes regulated by cytokinin in rice and those involved in biotic stress responses. To further explore this, we examined publicly available expression data related to pathogen interactions in rice for overlap with the cytokinin-related gene changes found here by RNA-Seq. We first examined gene expression changes in response to benzothiadiazole (BTH) , which activates the salicylic acid pathway, mimicking a plant’s natural defense response. More than half of BTH-regulated genes were also regulated by cytokinin in roots (Fig. 6a), suggesting that cytokinin activates a substantial portion of the BTH targets in roots. We also examined the overlap of the cytokinin-regulated transcriptome in roots with changes in the transcriptome in response to the plant pathogens Xanthomonas oryzae pv. oryzicola (Xoo)  and Magnaporthe oryzae , the causative agents of rice blight and rice blast respectively. Similar to the BTH results, nearly half of genes affected by the bacterial pathogen Xoo or by the fungal pathogen M. oryzae were also differentially expressed in the cytokinin RNA-Seq experiment (Fig. 6b-c). This indicates that cytokinin signaling alters a major fraction of the pathogen-regulated transcriptome, suggesting that cytokinin plays a crucial role in the regulation of genes driving the response to biotic stress.
Infection by M. oryzae has recently been shown to elevate iP-type cytokinins in rice leaf blades ; altered cytokinin levels have been observed in many other plant pathogen interactions . Further, M. oryzae itself has been shown to synthesize various cytokinin species . As was found in Arabidopsis , cytokinin acts synergistically with salicylic acid to activate defense responses in rice, in a manner dependent on both WRKY45 and OsNPR1 . The authors suggest that, similar to the interaction of Arabidopsis with the pathogen Hyaloperonospora arabidopsidis , M. oryzae likely elevates cytokinin in order to increase its proliferation by modification of host processes such as increasing sink activity; rice perceives this inappropriate modulation of cytokinin levels and activates defense responses .
Many resistance genes in plants encode nucleotide-binding site leucine-rich repeat (NBS-LRR) proteins, which are involved in the detection of diverse pathogens, including bacteria, viruses, fungi, nematodes, insects and oomycetes . Consistent with a link between cytokinin and pathogen response signaling, many NBS-LRR genes were regulated by cytokinin in rice, including 26 genes related to the Arabidopsis RPM1/(RPS3) gene and several other RPS homologs. Further, the induction of a large fraction of the receptor-like kinases (RLKs) in rice by cytokinin (Fig. 4a) may reflect the intersection with pathogen signaling as it has been postulated that the expansion of the rice RLK gene family after the rice/Arabidopsis split involves defense-related genes .
The WRKY family has been closely linked to biotic stress response. Many rice WRKYs are rapidly regulated transcriptionally upon pathogen infection [78, 83] and many have been shown to directly affect susceptibility to pathogens . In our RNA-Seq experiment, a large number of OsWRKYs were up-regulated by cytokinin, particularly in the root tissue (Fig. 7). The rice genome contains 98 to 102 OsWRKY transcription factors  and almost half (50/102) were differentially expressed in BA-treated rice roots, and most of these (45/50) were up-regulated. In rice, multiple WRKY loci have been linked to pathogen responses , and several of these are induced by cytokinin, including OsWRKY28, OsWRKY45, OsWRKY45, OsWRKY53, OsWRKY62, OsWRKY71, OsWRKY76 (Fig. 7). For example, overexpression of OsWRKY53 and OsWRKY71 increased resistance to the fungal and bacterial pathogens, respectively [87, 88]. Overexpression of OsWRKY28, OsWRKY62, OsWRKY71, and OsWRKY76 resulted in increased resistance to Xoo . The WRKY45 gene, which is induced 9.5-fold in response to cytokinin treatment in rice roots (Fig. 7), plays key role in the response to both Xoo and M. oryzae: overexpression of WRKY45 conferred strong resistance to the Xoo and M. oryzae ; OsWRKY45 is essential for BTH-primed plant immunity to M. oryzae and Xoo . WRKY45-2, WRKY13, and WRKY42 form a sequential transcriptional regulatory cascade required for resistance to M. oryzae . Further, the role of cytokinin in defense responses has been linked to WRKY45 as cytokinin and salicylic acid act synergistically to elevate genes encoding the enzymes involved in the biosynthesis of diterpenoid phytoalexin defense compounds in a WRKY45-dependent manner . Together, these data suggest that cytokinin plays an important role in facilitating the response to pathogen in rice, particularly by the regulation of the NBS-LRR receptors and the WRKY family of transcription factors. Analysis of rice lines disrupted for cytokinin metabolism and signaling should help elucidate the complex interactions among cytokinin, defense signaling and the NBS-LRR/WRKYs.
Here we demonstrate that there are differences in effects of cytokinin on the transcriptome of rice and Arabidopsis, including both a large increase in the number of genes induced by cytokinin in rice and differences in the suites of genes regulated, which suggests distinct cytokinin outputs in monocots. One of the most enriched sets of genes regulated by cytokinin in rice are linked to biotic responses, and there is substantial overlap between genes regulated by cytokinin and by various pathogens, suggesting a close link between cytokinin and pathogen responses in rice. These results begin to unravel the complex gene regulation after cytokinin perception in a crop of agricultural importance and provide insight into the processes and responses modulated by cytokinin in monocots.
Rice RNA-Seq and NanoString analysis
Rice seeds (Nipponbare) obtained from plants grown in our research greenhouse were surface sterilized with 2.5% sodium hypochlorite for 15 min and germinated on filter paper containing distilled water at 37 °C for one day. Germinated seeds were transferred to hydroponic culture containing Kimura B solution  in a growth chamber set to 14 h light (28 °C) and 8 h of dark (23 °C) with light intensity 700 μmol m−2 s−1. Plants were grown in two pots. At day 10 of hydroponic cultivation, culture media was replaced with new media containing 5 μM BA or 0.05 mN NaOH as a control. Following 120 min of the treatment, roots and shoots were harvested separately. Libraries were made from each of the four rice sample types (treated root, treated shoot, untreated root, and untreated shoot) using the ThruPLEX DNA-seq kit (Rubicon Genomics). The multiplexed samples were combined and sequenced on three lanes of an Illumina HiSeq flowcell for 100 cycles, yielding 100 base, single end reads. Sequencing was done at the David H. Murdock Research Institute in Kannapolis. Each lane contained one library from each of the four sample types.
The NanoString nCounter analysis was performed by the University of North Carolina Genomics and Bioinformatics Core Facility. 160 ng of total RNA was directly hybridized with gene specific color-coded probes and the data collection is carried out in the nCounter Digital Analyzer as described by the manufacturer (NanoStrings Technologies; Seattle, WA). The NanoString Codeset (Additional file 11: Table S10) was designed and synthesized by NanoString Technologies (Seattle, WA), which including three reference genes, LOC_Os01g16414 (Actin), LOC_Os02g32030 (Elongation factor), and LOC_Os05g01600 (Actin7). In addition, six positive-control and eight negative-control probes were added to each reaction to produce a standard curve for normalization. All the reaction counts were within the linear dynamic range of the standard curve. For each gene analyzed, the average plus two times standard deviation of the negative controls was subtracted from the raw data and then normalized to the standard curve within each reaction and three reference genes.
Sequences from rice were aligned onto the O. sativa japonica genome assembly Os-Nipponbar-Reference-IRGSP-1.0 (released October 2011)  using tophat (v.2.0.5)  and bowtie . For the tophat alignment, the default parameters were used except that the maximum intron size was set to 5000. The default settings used allowed up to two mismatches per read. The analysis was limited to reads with exactly one reported alignment, using the SAM format NH tag to distinguish such reads. Alignment, splice junction, and coverage graph files were deployed on an IGB QuickLoad site to facilitate visualization in Integrated Genome Browser . The number of single-mapping reads that overlapped each annotated protein-coding gene in the MSU Rice Genome Annotation Project (Release 7) genome annotations (RGAP7) was calculated using the samtools view –c function and supplied as inputs for differential expression analysis.
Differential expression analysis
Differentially expressed genes were identified using edgeR generalized linear modeling and exactTest methods. Differential expression was examined in roots and shoots (separately) between treated and untreated samples. In all cases, only genes with at least one read in one of the compared samples were considered (i.e. if a gene had no reads at all in any sample being compared, we did not test it for differential expression). An FDR cutoff of 0.0001 was used to determine differential expression of individual genes. Code used to process and analyze data is available from https://bitbucket.org/lorainelab/ricecyto.
Determination of rice Type-B DNA binding motifs
The DNA binding domains for the three rice type-B RRs (RR22, RR27, RR29) were fused to MBP by cloning into pMAL-c2x using restriction enzymes BamH I and Sal I. Primers used for cloning were as follows: CCCGGATCCTCAGCTGCAAAGAAGCCC and GTCGACTCAACCTAGTCTCTTGAGGTAAAG for RR22; CCCGGATCCAGGTTCACATGGACGACG and GTCGACTCATCTGTATTTCTGTAGATGGCT for RR27; CCCGGATCCACTAAGAAGAAATATTATCTCATG and GTCGACTCACAAATCCTTTGTTAGCCGTAG for RR29. Constructs were verified by sequencing and plasmids were introduced into the BL-21 strain of E. coli for protein expression. Expression of recombinant proteins was performed using the Zymo Dual Media Kit (Zymo Research). Induction was carried out using 0.1 mM IPTG for 8 h at 37 °C. Total soluble protein extracts from E. coli cultures expressing the recombinant proteins were analyzed by gel electrophoresis and immunoblot to confirm induction efficiency. OsRR22, OsRR27 and OsRR29 DNA binding specificities were then determined using protein binding microarrays (PBM11) as described by . The complete list of 25 8-mers and their E-scores can be found in the project repository PbaAnalysis/data/OsRR22_results.txt under PbaAnalysis/data/OsRR22_results.txt.
Gene Ontology and MapMan gene category enrichment analysis
Gene families with unusually many or few differentially expressed genes were detected using GOSeq, which accounts the tendency of genes with larger transcripts to be preferentially detected as differentially expressed in RNA-Seq data. GO term trees were visualized using Virtual Plant. MapMan gene to category assignments for MSU7 rice annotations were downloaded from the MapMan Web site at http://mapman.gabipd.org/ as plain text files. Gene Ontology annotations were obtained from the Gene Ontology Web site. Categories containing unusually many differentially expressed genes were identified using GOSeq; the code is available in the project repository PbaAnalysis/data/OsRR22_results.txt. The MapMan desktop application was used to visualize.
Histidine phosphotransfer proteins
Integrated Genome Browser
Reads per kilobase of expressed sequence per million mapped reads
Nemhauser JL, Hong F, Chory J. Different plant hormones regulate similar processes through largely nonoverlapping transcriptional responses. Cell. 2006;126:467–75.
Miller CO, Skoog F, Okomura FS, von Saltza MH, Strong FM. Isolation, structure and synthesis of kinetin, a substance promoting cell division. J Am Chem Soc. 1956;78:1345–50.
Miller CO, Skoog F, Von Saltza MH, Strong F. Kinetin, a cell division factor from deoxyribonucleic acid. J Am Chem Soc. 1955;77:1392.
Argueso CT, Ferreira FJ, Kieber JJ. Environmental perception avenues: the interaction of cytokinin and environmental response pathways. Plant Cell Environ. 2009;32:1147–60.
Mok DW, Mok MC. Cytokinin metabolism and action. Annu Rev Plant Physiol Plant Mol Biol. 2001;89:89–118.
Kieber JJ, Schaller GE. Cytokinins. Arabidopsis Book. 2014;12:e0168.
Frugier F, Kosuta S, Murray JD, Crespi M, Szczyglowski K. Cytokinin: secret agent of symbiosis. Trends Plant Sci. 2008;13:115–20.
Giulini A, Wang J, Jackson D. Control of phyllotaxy by the cytokinin-inducible response regulator homologue ABPHYL1. Nature. 2004;430:1031–4.
Caesar K, Thamm AMK, Witthöft J, Elgass K, Huppenberger P, Grefen C, Horak J, Harter K. Evidence for the localization of the Arabidopsis cytokinin receptors AHK3 and AHK4 in the endoplasmic reticulum. J Exp Bot. 2011;62:5571–80.
Lomin SN, Yonekura-Sakakibara K, Romanov GA, Sakakibara H. Ligand-binding properties and subcellular localization of maize cytokinin receptors. J Exp Bot. 2011;62:5149–59.
Inoue T, Higuchi M, Hashimoto Y, Seki M, Kobayashi M, Kato T, Tabata S, Shinozaki K, Kakimoto T. Identification of CRE1 as a cytokinin receptor from Arabidopsis. Nature. 2001;409:1060–3.
Ueguchi C, Sato S, Kato T, Tabata S. The AHK4 gene involved in the cytokinin-signaling pathway as a direct receptor molecule in Arabidopsis thaliana. Plant Cell Physiol. 2001;42:751–5.
Suzuki T, Miwa K, Ishikawa K, Yamada H, Aiba H, Mizuno T. The Arabidopsis sensor His-kinase, AHK4, can respond to cytokinins. Plant Cell Physiol. 2001;42:107–13.
Yamada H, Suzuki T, Terada K, Takei K, Ishikawa K, Miwa K, Yamashino T, Mizuno T. The Arabidopsis AHK4 histidine kinase is a cytokinin-binding receptor that transduces cytokinin signals across the membrane. Plant Cell Physiol. 2001;41:1017–23.
Hutchison CE, Li J, Argueso C, Gonzalez M, Lee E, Lewis MW, Maxwell BB, Perdue TD, Schaller GE, Alonso JM, et al. The Arabidopsis histidine phosphotransfer proteins are redundant positive regulators of cytokinin signaling. Plant Cell. 2006;18:3073–87.
Punwani JA, Hutchison CE, Schaller GE, Kieber JJ. The subcellular distribution of the Arabidopsis Histidine Phosphotransfer proteins is independent of cytokinin signaling. Plant J. 2010;62:473–82.
Mason MG, Mathews DE, Argyros DA, Maxwell BB, Kieber JJ, Alonso JM, Ecker JRS, Schaller GE. Multiple type-B response regulators mediate cytokinin signal transduction in Arabidopsis. Plant Cell. 2005;17:3007–18.
Ishida K, Yamashino T, Mizuno T. Expression of the cytokinin-induced type-A response regulator gene ARR9 is regulated by the circadian clock in Arabidopsis thaliana. Biosci Biotechnol Biochem. 2008;72:3025–9.
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.
Argyros RD, Mathews DE, Chiang YH, 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–16.
To JP, Deruère J, Maxwell BB, Morris VF, Hutchison CE, Schaller GE, Kieber JJ. Cytokinin regulates type-A Arabidopsis response regulator activity and protein stability via two-component phosphorelay. Plant Cell. 2007;19:3901–14.
To JPC, Haberer G, Ferreira FJ, Deruère J, Mason MG, Schaller GE, Alonso JM, Ecker JR, Kieber JJ. Type-A ARRs are partially redundant negative regulators of cytokinin signaling in Arabidopsis. Plant Cell. 2004;16:658–71.
Brandstatter I, Kieber JJ. Two genes with similarity to bacterial response regulators are rapidly and specifically induced by cytokinin in Arabidopsis. Plant Cell. 1998;10:1009–20.
D'Agostino I, Deruère J, Kieber JJ. Characterization of the response of the Arabidopsis ARR gene family to cytokinin. Plant Physiol. 2000;124:1706–17.
Mohanty S. Trends in global rice consumption. Rice Today. 2013;12:44–5.
Ashikari M, Sakakibara H, Lin S, Yamamoto T, Takashi T, Nishimura A, Angeles ER, Qian Q, Kitano H, Matsuoka M. Cytokinin oxidase regulates rice grain production. Science. 2005;309:741–5.
Tsai YC, Weir N, Hill K, Zhang W, Kim HJ, Shiu SH, Schaller E, Kieber JJ. Characterization of genes involved in cytokinin signaling and metabolism from rice. Plant Physiol. 2012;158:1666–84.
Du L, Jiao F, Chu J, Jin G, Chen M, Wu P. The two-component signal system in rice (Oryza sativa L.): A genome-wide study of cytokinin signal perception and transduction. Genomics. 2007;89:697–707.
Schaller GE, Doi K, Hwang I, Kieber JJ, Khurana JP, Kurata N, Mizuno T, Pareek A, Shiu S-H, Wu P, et al. Nomenclature for two-component signaling elements of rice. Plant Physiol. 2007;143:555–7.
Choi J, Lee J, Kim K, Cho M, Ryu H, An G, Hwang I. Functional identification of OsHk6 as a homotypic cytokinin receptor in rice with preferential affinity for iP. Plant Cell Physiol. 2012;53:1334–43.
Sun L, Zhang Q, Wu J, Zhang L, Jiao X, Zhang S, Zhang Z, Sun D, Lu T, Sun Y. Two rice authentic histidine phosphotransfer proteins, OsAHP1 and OsAHP2, mediate cytokinin signaling and stress responses in rice. Plant Physiol. 2014;165:335–45.
Hirose N, Makita N, Kojima M, Kamada-Nobusada T, Sakakibara H. Overexpression of a type-A response regulator alters rice morphology and cytokinin metabolism. Plant Cell Physiol. 2007;48:523–39.
Cheng X, Jiang H, Zhang J, Qian Y, Zhu S, Cheng B. Overexpression of type-A rice response regulators, OsRR3 and OsRR5, results in lower sensitivity to cytokinins. Genet Mol Res. 2010;9:348–59.
Rashotte AM, Carson SD, To JP, Kieber JJ. Expression profiling of cytokinin action in Arabidopsis. Plant Physiol. 2003;132:1998–2011.
Brenner WG, Romanov GA, Kollmer I, Burkle 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–33.
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–55.
Bhargava A, Clabaugh I, To JP, Maxwell BB, Chiang YH, Schaller GE, Loraine A, Kieber JJ. Identification of cytokinin-responsive genes using microarray meta-analysis and RNA-Seq in Arabidopsi. Plant Physiol. 2013;162:272–94.
Brenner WG, Schmülling T. Transcript profiling of cytokinin action in Arabidopsis roots and shoots discovers largely similar but also organ-specific responses. BMC Plant Biol. 2012;12:112.
Nicol JW, Helt GA, Blanchard SGJ, Raja A, Loraine AE. The Integrated Genome Browser: free software for distribution and exploration of genome-scale datasets. Bioinformatics. 2009;25:2730–1.
Kawahara Y, de la Bastide M, Hamilton J, Kanamori H, McCombie W, Ouyang S, Schwartz D, Tanaka T, Wu J, Zhou S, et al. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice. 2013;6:4.
Sakai H, Lee SS, Tanaka T, Numa H, Kim J, Kawahara Y, Wakimoto H, Yang CC, Iwamoto M, Abe T, et al. Rice annotation project database (RAP-DB): An integrative and interactive database for rice genomics. Plant Cell Physiol. 2013;54:e6.
Ramsköld D, Wang ET, Burge CB, Sandberg R. An abundance of ubiquitously expressed genes revealed by tissue transcriptome sequence data. PLoS Comput Biol. 2009;5:e1000598.
Loraine AE, McCormick S, Estrada A, Patel K, Qin P. RNA-seq of Arabidopsis pollen uncovers novel transcription and alternative splicing. Plant Physiol. 2013;162:1092–109.
Cobbett C, Goldsbrough P. Phytochelatins and metallothioneins: roles in heavy metal detoxification and homeostasis. Annu Rev Plant Biol. 2002;53:159–82.
Geiss GK, Bumgarner RE, Birditt B, Dahl T, Dowidar N, Dunaway DL, Fell HP, Ferree S, George RD, Grogan T, et al. Direct multiplexed measurement of gene expression with color-coded probe pairs. Nat Biotech. 2008;26:317–25.
Malkov V, Serikawa K, Balantac N, Watters J, Geiss G, Mashadi-Hossein A, Fare T. Multiplexed measurements of gene signatures in different analytes using the Nanostring nCounterTM Assay System. BMC Res Notes. 2009;2:80.
Jin J, Zhang H, Kong L, Gao G, Luo J. PlantTFDB 3.0: a portal for the functional and evolutionary study of plant transcription factors. Nucleic Acids Res. 2014;42:D1182–7.
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–92.
Hayakawa T, Kudo T, Ito T, Takahashi N, Yamaya T. ACT domain repeat protein 7, ACR7, interacts with a chaperone HSP18.0-CII in rice nuclei. Plant Cell Physiol. 2006;47:891–904.
Nakagawa M, Shimamoto K, Kyozuka J. Overexpression of RCN1 and RCN2, rice TERMINAL FLOWER 1/CENTRORADIALIS homologs, confers delay of phase transition and altered panicle morphology in rice. Plant J. 2002;29:743–50.
Li L, Stoeckert CJ, Roos DS. OrthoMCL: Identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89.
Hosoda K, Imamura A, Katoh E, Hatta T, Tachiki M, Yamada H, Mizuno T, Yamazaki T. Molecular structure of the GARP family of plant Myb-related DNA binding motifs of the Arabidopsis response regulators. Plant Cell. 2002;14:2015–29.
Sakai H, Aoyama T, Oka A. Arabidopsis ARR1 and ARR2 response regulators operate as transcriptional activators. Plant J. 2000;24:703–11.
Imamura A, Kiba T, Tajima Y, Yamashino T, Mizuno T. In vivo and in vitro characterization of the ARR11 response regulator implicated in the His-to-Asp phosphorelay signal transduction in Arabidopsis thaliana. Plant Cell Physiol. 2003;44:122–31.
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–77.
Zürcher E, Tavor-Deslex D, Lituiev D, Enkerli K, Tarr PT, Müller B. A robust and sensitive synthetic sensor to monitor the transcriptional output of the cytokinin signaling network in planta. Plant Physiol. 2013;161:1066–75.
Franco-Zorrilla JM, López-Vidriero I, Carrasco JL, Godoy M, Vera P, Solano R. DNA-binding specificities of plant transcription factors and their potential to define target genes. Proc Natl Acad Sci U S A. 2014;111:2367–72.
Weirauch Matthew T, Yang A, Albu M, Cote AG, Montenegro-Montero A, Drewe P, Najafabadi Hamed S, Lambert Samuel A, Mann I, Cook K, et al. Determination and inference of eukaryotic transcription factor sequence specificity. Cell. 2014;158:1431–43.
Godoy M, Franco-Zorrilla JM, Pérez-Pérez J, Oliveros JC, Lorenzo Ó, Solano R. Improved protein-binding microarrays for the identification of DNA-binding specificities of transcription factors. Plant J. 2011;66:700–11.
Schoor S, Farrow S, Blaschke H, Lee S, Perry G, von Schwartzenberg K, Emery N, Moffatt B. Adenosine kinase contributes to cytokinin interconversion in Arabidopsis. Plant Physiol. 2011;157:659–72.
Brenner WG, Ramireddy E, Heyl A, Schmülling T. Gene regulation by cytokinin in Arabidopsis. Front Plant Sci. 2012;3:8.
Xiong Y, Liu T, Tian C, Sun S, Li J, Chen M. Transcription factors in rice: A genome-wide comparative analysis between monocots and eudicots. Plant Mol Biol. 2005;59:191–203.
Smalle J, Kurepa J, Yang P, Babiychuk E, Kushnir S, Durski A, Vierstra RD. Cytokinin growth responses in Arabidopsis involve the 26S proteasome subunit RPN12. Plant Cell. 2002;14:17–32.
Li Y, Kurepa J, Smalle J. AXR1 promotes the Arabidopsis cytokinin response by facilitating ARR5 proteolysis. Plant J. 2013;74:13–24.
Kurepa J, Li Y, Smalle JA. Cytokinin signaling stabilizes the response activator ARR1. Plant J. 2014;78:157–68.
Kim HJ, Chiang Y-H, Kieber JJ, Schaller GE. SCF(KMD) controls cytokinin signaling by regulating the degradation of type-B response regulators. Proc Natl Acad Sci U S A. 2013;110:10028–33.
Kim K, Ryu H, Cho YH, Scacchi E, Sabatini S, Hwang I. Cytokinin-facilitated proteolysis of ARABIDOPSIS RESPONSE REGULATOR 2 attenuates signaling output in two-component circuitry. Plant J. 2012;69:934–45.
Vogel JP, Woeste KW, Theologis A, Kieber JJ. Recessive and dominant mutations in the ethylene biosynthetic gene ACS5 of Arabidopsis confer cytokinin insensitivity and ethylene overproduction, respectively. Proc Natl Acad Sci U S A. 1998;95:4766–71.
Chae HS, Faure F, Kieber JJ. The eto1, eto2 and eto3 mutations and cytokinin treatment elevate ethylene biosynthesis in Arabidopsis by increasing the stability of the ACS5 protein. Plant Cell. 2003;15:545–59.
Marhavý P, Duclercq J, Weller B, Feraru E, Bielach A, Offringa R, Friml J, Schwechheimer C, Murphy A, Benková E. Cytokinin controls polarity of PIN1-dependent auxin transport during lateral root organogenesis. Curr Biol. 2014;24:1031–7.
Takagi S, Nomoto K, Takemoto T. Physiological aspect of mugineic acid, a possible phytosiderophore of graminaceous plants. J Plant Nutr. 1984;7:469–77.
Römheld V, Marschner H. Evidence for a specific uptake system for iron phytosiderophores in roots of grasses. Plant Physiol. 1986;80:175–80.
Nozoye T, Nagasaka S, Kobayashi T, Takahashi M, Sato Y, Sato Y, Uozumi N, Nakanishi H, Nishizawa NK. Phytosiderophore efflux transporters are crucial for iron acquisition in graminaceous plants. J Biol Chem. 2011;286:5446–54.
Ptashnyk M, Roose T, Jones DL, Kirk GJD. Enhanced zinc uptake by rice through phytosiderophore secretion: a modelling study. Plant Cell Environ. 2011;34:2038–46.
Thimm O, Bläsing O, Gibon Y, Nagel A, Meyer S, Krüger P, Selbig J, Müller 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–39.
Shimono M, Sugano S, Nakayama A, Jiang CJ, Ono K, Toki S, Takatsuji H. Rice WRKY45 plays a crucial role in benzothiadiazole-inducible blast resistance. Plant Cell. 2007;19:2064–76.
Zhou Y-L, Xu M-R, Zhao M-F, Xie X-W, Zhu L-H, Fu B-Y, Li Z-K. Genome-wide gene responses in a transgenic rice line carrying the maize resistance gene Rxo1 to the rice bacterial streak pathogen, Xanthomonas oryzae pv. oryzicola. BMC Genomics. 2010;11:78.
Wei T, Ou B, Li J, Zhao Y, Guo D, Zhu Y, Chen Z, Gu H, Li C, Qin G, et al. Transcriptional profiling of rice early response to Magnaporthe oryzae identified OsWRKYs as important regulators in rice blast resistance. PLoS One. 2013;8:e59720.
Jiang C-J, Shimono M, Sugano S, Kojima M, Liu X, Inoue H, Sakakibara H, Takatsuji H. Cytokinins act synergistically with salicylic acid to activate defense gene expression in rice. Mol Plant Microbe Interact. 2012;26:287–96.
Argueso CT, Ferreira FJ, Hutchison CE, To JPC, Epple P, Mathews DE, Schaller GE, Dangl JL, Kieber JJ. Two-component elements mediate interactions between cytokinin and salicylic acid in plant immunity. PLoS Genet. 2012;8:e1002448.
McHale L, Tan X, Koehl P, Michelmore RW. Plant NBS-LRR proteins: adaptable guards. Genome Biol. 2006;7:212.
Shiu SH, Karlowski WM, Pan R, Tzeng YH, Mayer KF, Li WH. Comparative analysis of the receptor-like kinase family in Arabidopsis and rice. Plant Cell. 2004;16:1220–34.
Ryu H-S, Han M, Lee S-K, Cho J-I, Ryoo N, Heu S, Lee Y-H, Bhoo S, Wang G-L, Hahn T-R, et al. A comprehensive expression analysis of the WRKY gene superfamily in rice plants during defense response. Plant Cell Rep. 2006;25:836–47.
Eulgem T, Somssich IE. Networks of WRKY transcription factors in defense signaling. Curr Opin Plant Biol. 2007;10:366–71.
Ross CA, Liu Y, Shen QJ. The WRKY gene family in rice (Oryza sativa). J Integr Plant Biol. 2007;49:827–42.
Cheng H, Liu H, Deng Y, Xiao J, Li X, Wang S. The WRKY45-2 WRKY13 WRKY42 transcriptional regulatory cascade is required for rice resistance to fungal pathogen. Plant Physiol. 2015;167:1087–99.
Chujo T, Takai R, Akimoto-Tomiyama C, Ando S, Minami E, Nagamura Y, Kaku H, Shibuya N, Yasuda M, Nakashita H, et al. Involvement of the elicitor-induced gene OsWRKY53 in the expression of defense-related genes in rice. Biochim Biophys Acta. 2007;1769:497–505.
Liu X, Bai X, Wang X, Chu C. OsWRKY71, a rice transcription factor, is involved in rice defense response. J Plant Physiol. 2007;164:969–79.
Peng Y, Bartley LE, Canlas P, Ronald PC. OsWRKY IIa transcription factors modulate rice innate immunity. Rice. 2010;3:36–42.
Shimono M, Koga H, Akagi AYA, Hayashi N, Goto S, Sawada M, Kurihara T, Matsushita A, Sugano S, Jiang C-J, et al. Rice WRKY45 plays important roles in fungal and bacterial disease resistance. Mol Plant Pathol. 2012;13:83–94.
Nakayama A, Fukushima S, Goto S, Matsushita A, Shimono M, Sugano S, Jiang C-J, Akagi A, Yamazaki M, Inoue H, et al. Genome-wide identification of WRKY45-regulated genes that mediate benzothiadiazole-induced defense responses in rice. BMC Plant Biol. 2013;13:150.
Akagi A, Fukushima S, Okada K, Jiang C-J, Yoshida R, Nakayama A, Shimono M, Sugano S, Yamane H, Takatsuji H. WRKY45-dependent priming of diterpenoid phytoalexin biosynthesis in rice and the role of cytokinin in triggering the reaction. Plant Mol Biol. 2014;86:171–83.
Ma JF, Goto S, Tamai K, Ichii M. Role of root hairs and lateral roots in silicon uptake by rice. Plant Physiol. 2001;127:1773–80.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
National Science Foundation Plant Genome Research Program award number IOS-1238051 to JJK, AEL, and GES and the Spanish Ministry for Science and Innovation grant BIO2013-44407-R and Fundacion UAM Grant 2015007 to RS funded this work. NIH award 2R01GM103463-06A1 to AL supports development of the IGB software.
Availability of data and materials
Illumina sequence data are available from the Short Read Archive under accession SRP049054 (rice data) and SRP059384 (the Arabidopsis data). Processed data files are available for visualization in Integrated Genome Browser via IGB QuickLoad data source http://www.igbquickload.org/cyto. Code used for statistical analysis is version-controlled in git repositories located at http://www.bitbucket.org/lorainelab/ricecyto, https://bitbucket.org/lorainelab/atcytokinin2 and https://bitbucket.org/lorainelab/pbamotif.
ICB, TR, AEL, and JJK did data analysis. YCT did RNA-Seq and NanoStrings experiments. JMW, JMFZ, and RS did the OsRR22 experiments. TR, ICB, JMW, JJK, GES and AEL wrote the manuscript. All authors have read and approved this manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Number of rice RNA-Seq reads mapping to genes by library. (XLSX 12 kb)
Rice RNA-Seq gene expression values in RPKM. (XLSX 13 kb)
Genes that were DE between roots and shoots in rice. (XLSX 3612 kb)
Genes differentially regulated in response to cytokinin in roots. (XLSX 10749 kb)
Genes differentially regulated in response to cytokinin in shoots. (XLSX 1485 kb)
Genes that were DE in roots and shoots in the same direction. (XLSX 419 kb)
Regulation of genes involved in cytokinin biosynthesis, signaling, response, and degradation. Figure S2. Cytokinin regulation of orthologous genes in rice and Arabidopsis. Figure S3. Characterization of rice type-B RRs OsRR27 and OsRR29 through use of protein binding microarrays. Figure S4. GO categorization of genes induced by cytokinin in both the root and the shoot. Figure S5. GO categorization of genes repressed by cytokinin in both the root and the shoot. Figure S6. GO categorization of genes induced by cytokinin in the root. Figure S7. GO categorization of genes induced by cytokinin in the shoot. Figure S8. GO categorization of genes repressed by cytokinin in the shoot. Figure S9. MapMan overviews of metabolism, cellular response and regulation of cytokinin-regulated genes in shoots in rice. (XLSX 220 kb)
Genes that were DE in roots and shoots in opposite directions. (XLSX 58 kb)
MapMan categories and DE genes from roots. (PDF 1088 kb)
MapMan categories and DE genes from shoots. (XLSX 26 kb)
Primers used in the NanoString analysis. (XLSX 44 kb)
About this article
Cite this article
Raines, T., Blakley, I.C., Tsai, Y. et al. Characterization of the cytokinin-responsive transcriptome in rice. BMC Plant Biol 16, 260 (2016) doi:10.1186/s12870-016-0932-z
- Plant hormone