Gene Expression Analysis of Rice Seedling under Potassium Deprivation Reveals Major Changes in Metabolism and Signaling Components

Plant nutrition is one of the important areas for improving the yield and quality in crops as well as non-crop plants. Potassium is an essential plant nutrient and is required in abundance for their proper growth and development. Potassium deficiency directly affects the plant growth and hence crop yield and production. Recently, potassium-dependent transcriptomic analysis has been performed in the model plant Arabidopsis, however in cereals and crop plants; such a transcriptome analysis has not been undertaken till date. In rice, the molecular mechanism for the regulation of potassium starvation responses has not been investigated in detail. Here, we present a combined physiological and whole genome transcriptomic study of rice seedlings exposed to a brief period of potassium deficiency then replenished with potassium. Our results reveal that the expressions of a diverse set of genes annotated with many distinct functions were altered under potassium deprivation. Our findings highlight altered expression patterns of potassium-responsive genes majorly involved in metabolic processes, stress responses, signaling pathways, transcriptional regulation, and transport of multiple molecules including K+. Interestingly, several genes responsive to low-potassium conditions show a reversal in expression upon resupply of potassium. The results of this study indicate that potassium deprivation leads to activation of multiple genes and gene networks, which may be acting in concert to sense the external potassium and mediate uptake, distribution and ultimately adaptation to low potassium conditions. The interplay of both upregulated and downregulated genes globally in response to potassium deprivation determines how plants cope with the stress of nutrient deficiency at different physiological as well as developmental stages of plants.


Introduction
Potassium is one of the essential macronutrients required for plant growth and development. It plays a major role in different physiological processes like cell elongation, stomatal movement, turgor regulation, osmotic adjustment, and signal transduction by acting as a major osmolyte and component of the ionic environment in the cytosol and subcellular organelles [1][2][3][4][5][6][7]. Potassium is also required for balancing the electrical charge of membranes, energy generation by proton pump activity, longdistant transport of ions from root to shoot, protein synthesis, enzyme activation, and metabolism of sugars and nitrogen [2,8,9]. Since potassium is one of the major plant macronutrients (cytosolic K + concentration is approximately 100 mM), potassium deficiency poses a severe agricultural challenge and requires the use of large quantities of chemical fertilizers for sustainable agricultural practices.
Previous studies report that potassium acts as an activator or cofactor in several enzyme systems [10]. Enzymatic activity of pyruvate kinase, starch synthase, nitrate reductase, and rubisco are all directly connected to metabolic changes under potassium deficiency [11][12][13][14]. One of the hallmarks of potassium deficiency is chlorosis (yellowing) in older leaves, a consequence of mobilization of potassium from older leaves to younger growing tissues [2,9,15].
Potassium uptake takes place in the roots of plants; potassium is subsequently redistributed to plant tissues and organs and stored in abundance in vacuoles. Plant roots tolerate short-term potassium deprivation by utilizing potassium stored in the vacuole when available. When plants grow in potassium-deficient soil, the root cells sense the low concentrations of K + and initiate a series of physiological reactions [16,17]. The detailed physiological role of potassium absorption and uptake has been studied in several plant species, and the molecular mechanisms of potassium transport have been largely elucidated in Arabidopsis. A large number of transporters and channels in Arabidopsis have been implicated in the uptake and mobilization of potassium from root to other parts of the plant [4,[18][19][20][21]. To adjust fluctuation of potassium levels in the soil, plants have adopted two modes of potassium uptake in their roots, namely high-affinity and low-affinity uptake [4,22,23]. Recently, studies implicated calcium-mediated CBL-CIPK signaling in regulating the shaker family K + channels AKT1 and AKT2, however detailed mechanistics of potassium sensing remain elusive [6,24,25]. Research related to the molecular mechanisms of K + sensing, uptake, distribution, and homeostasis in cereal and non-cereal crops is still miniscule. Although considerable work has been done in the model plant Arabidopsis, an extensive amount of work is still needed in crop plants to understand the detail mechanisms of K + nutrition and signaling.
In this study, we used whole genome microarrays to determine the transcriptomic profile of rice seedlings exposed to short-term K + deficiency followed by K + resupply. We applied Benjamini-Hochberg correction to filter the differentially expressed genes in different conditions. We also performed PCA (Principal Component Analysis) and Pearson correlation coefficient analysis to ensure reliability of the data. According to our microarray data, potassium deficiency affects the expression of various genes, which were grouped into different categories such as metabolism, transcription factor, transporter, signal transduction. The objective of this study is to determine the expression pattern of genes under low potassium stress conditions and to determine the role of these genes in potassium homeostasis and adaptation in potassiumdeficient soil.

Phenotypic Analysis of K + Deficiency on Rice Seedling Growth
For transcriptional analysis under potassium deficiency, 5 days old hydroponically grown indica rice IR64 seedlings were transferred to a potassium-deficient medium for 5 days followed by resupply of potassium for 6 h. Seedlings grown in potassiumdeficient media had reduced shoot and root growth and leaf blade height ( Figure 1A, B). However, there was no phenotypic difference between K + minus or starved (KM) and K + resupply (KR) seedlings after treatment for 6 h. In addition, we observed a higher number of root hairs in KM seedlings while more lateral roots were seen in K + plus or normal growth (KP) seedlings ( Figure 1A, B). The biomass and ion content (K + , Na + ) of seedlings were also measured for each condition. The fresh and dry weight of KM seedlings was lower after 5 days of K + deficiency compared to normal growth conditions, whereas no significant difference was observed between KM and KR ( Figure 1C). In order to determine the total K + and Na + content, single-channel flame photometry (see Material and Methods) was performed for all the three conditions. KM seedlings contained about 4 times less K + than KP seedlings while KR showed slightly higher K + than KM plants ( Figure 1D). Na + content was lower in KP condition as compared to KM and KR condition, however no significant difference was found in KM and KR seedlings.

Global Expression Analysis in Response to Potassium Deficiency and Resupply Condition
To obtain insights into changes in rice gene expression profiles under potassium deprivation and resupply of potassium, we performed rice whole genome microarrays using Affymetrix gene chip (57 K). Rice seedlings were subjected to three different conditions, viz. K + plus (KP), K + minus (KM), and K + resupply (KR) after five days of normal growth as described in the material and method. Total RNA was isolated from treated seedlings and used for whole genome microarray analysis. The microarray gene expression data was normalized against data obtained for samples grown on normal media for five days (untreated seedlings) in order to eliminate potassium-responsive genes involved in early growth and development of seedlings. To determine the accuracy and reliability of the data obtained, microarray data of three biological replicates for each treatment was analysed for their correlation coefficient, and Principal Component Analysis (PCA) was also performed. Correlation between the biological replicates was calculated using Pearson's correlation coefficient based on the signal intensities. The three replicates of each treatment showed correlation coefficient of greater than 0.98 (Table 1). The values of the correlation coefficients of the three replicates were further reinforced by the result of the PCA, wherein the biological replicates for each treatment showed a high degree of correlation. Interestingly, the PCA cluster for the third treatment i.e. resupply of K + -deficient medium with K + for 6 h (KR) was positioned closer to the PCA cluster of the K + -plus condition (KP) as compared to the PCA cluster for the K + -minus condition (KM) (Figure 2A). This indicated that resupplying potassium externally for 6 h restored expression of many genes to their normal levels [i.e. K + plus condition (KP)].
After normalization of the microarray data, differential expression analysis was performed for the following pairwise combinations of the three test conditions: KP and KM (KP/KM) and KM and KR (KM/KR). Differentially expressed genes were identified on the basis of following two criteria: first, the fold change, wherein genes with a fold change of $2 were considered for analysis and second, the unpaired t-test, wherein genes with pvalue #0.05 were selected for analysis. False discovery rate (FDR) correction (Benjamini-Hochberg Correction) was applied, which added greater significance to the data. 722 genes were differentially expressed under potassium deficient (KM) condition while 1867 genes were differentially expressed upon potassium resupply (KR). The details (probe set Ids, signal intensity value, p-value, fold change and regulation) of 722 and 1867 genes is given in Table S1 and Table S2 respectively. Among these 722 differentially expressed genes, 281 genes showed up regulation and 441 genes were down regulated ( Figure 3A). Out of the 1867 genes differentially expressed upon potassium resupply, 866 showed up regulation, while 1001 genes were down regulated ( Figure 3B).
A total of 307 differentially expressed genes were common for both KP and KR condition ( Figure 2B, Figure 3C), details of which are given in Table S3. For the validation of microarray data, 19 significantly and differentially expressed genes representing different categories were selected for quantitative realtime PCR (qPCR) analysis. The results of qPCR of the 19 tested genes were in accordance with the microarray results, which authenticated our microarray data and established its reliability ( Figure 4, Table S4). The raw data sets (CEL) and the normalized expression data sets have been deposited in the Gene Expression Omnibus at the National Center for Biotechnology Information with the accession number GSE44250 (http://www.ncbi.nlm.nih.gov/ geo/).

Functional Classification of the Differentially Expressed Genes
The differentially expressed genes in potassium-deficient conditions were functionally classified by homology search against the Gene Ontology (GO) and NCBI Non-redundant (NR) databases using BLAST through NCBI. Genes encoding hypothetical proteins were classified as genes of unknown function. Among the 722 genes that were differentially expressed under potassium deprivation, 15% of the genes belonged to the unknown function category ( Figure 5). The remaining 85% genes were categorised into 17 comprehensive subdivisions corresponding to the following functions: primary metabolism, secondary metabolism, nucleic acid metabolism, transporters, transcription factors, auxin signaling components, cell wall metabolism, cell death, growth and development, photosynthesis, stress responses, etc. ( Figure 5). Genes that could not be classified into any of the above-mentioned categories were placed in the category of ''others''. Overall, this analysis indicate that potassium deficiency affects the expression of genes with diverse functions in metabolism, transcriptional regulation, stress responses, molecular transport and signal transduction. In all, nearly 56% of total differentially expressed genes are involved in these processes according to their annotation.
18% of differentially expressing genes (DEG) were involved in various metabolic pathways, including 4% in nucleic acid metabolism, 3% in secondary metabolism and 2% in protein metabolism, carbohydrate metabolism (primary metabolism), lipid metabolism, glycolytic enzymes and other biosynthetic pathways (B) Length of rice shoot, root and leaf blade under three different conditions. (C) Fresh and dry weight comparison of rice seedling. (D) K + and Na + content of rice seedlings. Differences between mean values of treatments and controls were compared using t -tests (* P,0.05, ** P,0.01). doi:10.1371/journal.pone.0070321.g001 (Table 2). Of all the DEG, 7% were transporters, including three potassium transporters HAK1, OsHKT2;3 and OsHKT2;4; ABC transporter; 3 peptide transporter PTR2 and PTR3-A, and many other diverse transporters. Diverse group of kinases (CIPKs and MAPKs), phosphatases (PP2C and other Ser/Thr protein phosphatases), and calcium sensors (CMLs and CBLs) known to be involved in signal transduction mechanism constituted 5% of DEG. Significant percentages of the total DEG were comprised of genes annotated as transcription factors (9%) and stress responsive genes (8%). A small group of genes related to photosynthesis (3%), growth and development (3%), cell wall (2%), defence (2%), and auxin (1%) were also differentially expressed under potassiumdeficient conditions.

K-means Cluster Analysis
To elucidate potassium-specific gene expression profiles, coregulated genes were further analysed by K-Means Cluster analysis. Generally, K-means cluster analysis is based on the assumption that genes involved in either a similar function or a common pathway will have similar expression profiles and can likely to be grouped together. The analysis was performed for genes that were significantly expressed during potassium starvation. A total of 722 DEG under potassium-deficient conditions were clustered into 16 primary clusters and were then grouped into 7 different groups based on significant changes in expression pattern ( Figure 6). Group 1, the largest group, consisted of 151 genes that were further divided into four sub-groups (i.e. 1A, 1B, 1C and 1D). All the genes in this group exhibited high expression levels under potassium-deficient conditions (KM), while their expression was low in both potassium-sufficient (KP) and potassium resupply (KR) conditions. Stress-related genes such as oxidoreductase, GST, heat shock protein and a few genes related to metabolism such as oligosaccharyl transferase, terpene synthase, ribulose 1,5 biphosphate were included in this group. 32 and 51 genes from group 2 a and 2 b, respectively, showed similar expression pattern as group 1 under KM condition, but they exhibited slightly elevated expression during potassium resupply (KR) compared to K + plus (KP) conditions. This group included genes related to plant growth, plant development, and transcription factors. In group 3 (133 genes) and group 4 (94 genes), the overall expression pattern of genes was similar, but the difference in expression levels of genes under KM and KR was greater in group 4 as compared to group 3. Group 3 included 133 genes; all of them showed elevated expression in KP and low expression in KM. After resupply of potassium, they regained expression levels comparable to KP. Most of the photosynthesis, signal transduction, and transporter-related genes along with flavonol synthase, expansin precursor, and RNA polymerase sigma factor, comprised group 3 and 4. 163 genes of group 5 had elevated expression in K + -sufficient (KP) conditions and constant expression in K + -deficient (KM) and K + -resupply condition (KR). CHIT10, WRKY 50, WRKY 65, a MYB family transcription factor were also contained in this group. Group 6 genes were not differentially expressed in any of the three conditions tested and maintained a constant low level of expression. On the other hand, genes placed in group 7 showed high expression in all three conditions.

Genes Involved in Metabolism
A large number of DEG was found to be involved in metabolic processes, particularly primary metabolism (carbohydrate and lipid), secondary metabolism, protein and nucleic acid metabolism. Potassium is a crucial ion for plants as it acts as a cofactor for various enzymes in metabolic pathways [10]. From this gene expression analysis, many genes related to metabolism showed altered expression in potassium-deficient conditions (KM). Upon resupply of potassium, the expression was restored to similar levels as seen in K + -sufficient (KP) conditions. As mentioned earlier, 18% of differentially regulated genes were involved in various metabolic and biosynthetic processes. Among these 18% of genes, 30 genes are related to carbohydrate metabolism, of which 10 genes were downregulated in KM condition and upregulated upon external resupply of potassium (KR). Another 14 genes showed the reverse trend of upregulation in KM condition and downregulation in KR condition (Table S5). Carboxyvinyl-carboxyphosphonate phosphorylmutase (LOC_Os12g08760) showed downregulation in potassium-deficient conditions, as well as glyoxalase (LOC_Os07g06660), pyruvate kinase (LOC_Os11g10980), and lactate/malate dehydrogenase (LOC_Os07g43700). UDP-glucoronosyl and UDP-glucosyl transferases (LOC_Os01g08440, LO-C_Os02g11640), and a glycosyltransferase (LOC_Os01g07530) showed downregulation in both deficient and resupply condition. Phosphoglycerate kinase (LOC_Os06g45710) and indole-3-glycerol phosphate synthase (LOC_Os09g08130) were up-and downregulated in KM and KR condition, respectively. Alpha- amylase (LOC_Os09g28430, LOC_Os02g52700) showed the highest expression in deficient condition ( Table 2). Eleven of the lipid metabolism-associated genes were downregulated in potassium-deficient conditions and showed elevated expression in resupply conditions (Table S5). Beta-amylase (LOC_Os03g04770), omega-3 fatty acid desaturase (LO-C_Os03g18070), and lipase (LOC_Os04g56240) exhibited low expression in deficient condition and high expression upon resupply. Three genes related to lipid metabolism showed the reverse trend and were upregulated in KM and downregulated in KR condition. A number of other genes related to metabolism such as transferases (10 genes), dehydrogenases (5 genes) other kinases/phosphatases (4 genes) were also differentially regulated (Table S5).

Genes Involved in Secondary Metabolism
Eleven K + -responsive genes related to secondary metabolism were also downregulated upon KM treatment and upregulated in KR conditions. This category included many enzymes such as two polyprenyl synthetases (LOC_Os05g50550, LOC_Os04g56230) and a protein containing a FAD-binding domain (LOC_Os02g51080).

Nucleic Acid and Protein Metabolism
In potassium-deficient conditions, the majority of DEG not only belonged to carbohydrate and lipid metabolism, but several genes were identified that are annotated as being involved with nucleic acid and protein metabolism. Among 22 genes associated with nucleic acid metabolism, only 8 showed downregulation in KM condition and elevated expression in KR condition, whereas the remaining exhibited reversed expression pattern (Table S5). Transcripts presenting this profile encoded, for example, two AMP deaminases (LOC_Os05g28180, LOC_Os07g46630); endonuclease/exonuclease/phosphatase domain containing (LO-C_Os01g58690); RNA polymerase sigma factor (LO-C_Os08g06630, LOC_Os05g51150); and 39-59 exonuclease (LOC_Os04g14810). Most of the protein metabolism-related transcripts shown in Table 2 were downregulated upon potassium treatment, and their expression increased upon KR treatment. The putative C-terminal processing peptidase OsCttP3 (LO-C_Os06g21380), the putative prolyl oligopeptidase OsPOP9 (LOC_Os04g47360), aspartic proteinase nepenthesin-2 (LO-C_Os06g03120), and 8 other genes were included in this category.

Genes Related to Photosynthesis and Plant growth
Photosynthesis is a process by which green plants capture light energy to synthesize organic compounds from carbon dioxide and water. Seventeen photosynthesis-related genes were expressed differentially under KM and KR conditions. Surprisingly, all 17 genes were downregulated in KM condition and upregulated in KR condition. Some of these genes (9) encoded chlorophyll a/b binding protein, photosystem I reaction centre subunit (LO-C_Os09g30340), and magnesium-protoporphyrin IX monomethyl ester cyclase (LOC_Os01g17170) ( Table 2). Sixteen transcripts (9 upregulated and 7 downregulated) that showed altered expression in potassium-deficient conditions relative to their expression under normal growth conditions were related to plant growth and development. The expression of these genes was reversed upon resupply of potassium (Table S5).

Transcription Factors and Related Genes
In higher plants, gene expression in response to various stresses is regulated by transcription factors [26][27][28]. In many biological pathways, transcription factors are considered to be master regulators, since they can control the expression of suites of multiple genes at a given time in response to a particular stress condition. Our microarray results showed that 55 genes encoding transcription factors were induced under potassium deficiency  (Table 2).

Transporters
Potassium deficiency affected the expression of a large number of transcripts encoding transporters representative of several distinct protein families. Similar to transcription factors, some transporters also showed downregulation under potassium deficiency and upregulation upon resupply. Of the 40 differentially expressed transporters, 9 were metal transporters; 7 were peptide transporters; and 5 were lipid transporters. In this study, sodium transporters (OsHKT2;3 and OsHKT2;4 2 ), sulphate transporters (2 genes), and ABC transporters (2 genes) were highly downregulated in KM conditions and upregulated under KR conditions ( Table 2). Thirteen other transporters (annotated as aquaporins, MATE efflux transporters, and AAA-type ATPases) also displayed similar expression pattern as described above. Surprisingly, we did not find any potassium transporters that showed differential expression in our microarray data with the criteria of more than 2-fold change. But when the fold change criterion was reduced to 1.5, two potassium transporters (OsHAK1), one potassium channel (OsKAT1) and one putative potassium transporter were found to be differentially expressed (Table S6).

Signal Transduction Related Genes
Based on the available literature for analyses of gene expression changes in response to various nutrient deficiencies, several proteins involved in signal transduction networks such as calcium sensors, kinases, and phosphatases were found to be differentially regulated under those conditions [24,25,29,30]. Kinases and phosphatases are involved in a large number of distinct signaling pathways, such as the signal transduction pathways controlling cell growth, differentiation, and death. In our analysis, 31 differentially expressed genes were classified as signal transduction-related genes. Three CIPKs (CIPK7, CIPK8 and CIPK29) and three MAPKs display distinct patterns of downregulation and upregulation in potassium-deficient and resupply conditions, respectively (Table S5). Among the other signaling related genes, 6 phosphatases, 2 PP2Cs (LO-C_Os01g46760, LOC_Os04g33080), and a phosphocholine phosphatase exhibited low expression under KM condition and elevated expression under KR condition, whereas one of the Ser/Thr protein phosphatases and two of the pyrophosphatases exhibited the reverse expression pattern. Some of the kinases such as a CDPK (LOC_Os03g48270) and the serine/ threonine protein kinase SNT7 (LOC_Os05g47560) also showed downregulation in KM condition, but BRASSINOSTEROID INSENSITIVE 1 (LOC_Os11g31530) was upregulated 12-fold under KM conditions and down regulated 10-fold under KR conditions ( Table 2).

Stress Related Genes
Reactive oxygen species (ROS) are known to be involved in signaling pathways specific to low potassium stress conditions [31,32]. Oxidoreductase and glutathione S-transferase (GST) play important roles during oxidative stress in plants [33]. In our study, a total of 50 genes grouped as stress-related were differentially expressed under low potassium conditions. These included 28 upregulated and 22 downregulated genes. Among these, 5 oxidoreductase genes and 4 GST were identified that showed upregulation in KM conditions and downregulation in KR conditions (Table S5), which we view as unsurprising given the well established role of potassium in stress responses [34][35][36]. A total of 11 heat stress-related genes (HSPs) and 5 dehydrationrelated genes were identified. Most of the HSPs (9 genes) were downregulated in potassium-deficient conditions, and the expression pattern was reversed upon resupply of potassium. Other drought-responsive genes showed a qualitatively similar expression pattern (Table S2).
Cell death-related transcripts constituted only a small category in the microarray results. This category contained 7 genes, including 5 genes encoding ubiquitin, which mostly showed upregulation under low potassium conditions (KM) ( Table S5).

Phytohormones and Defense Related Genes
Changes in potassium supply affect the transcription of genes related to phytohormones like jasmonic acid (JA), auxin, gibberellin and others [5,37]. Seven genes were found to be related to auxin, such as OsSAUR21, OsSAUR53, OsIAA2 (Table 2, S5). Interestingly, two jasmonic acid-related genes, lipoxygenase (LOC_Os02g10120) and jacalin-like lectin (LOC_Os12g14440) were upregulated 16-and 65-fold in potassium-deficient conditions, respectively. Their expression levels declined rapidly after resupply. Two GA-related genes also showed similar expression patterns as the jasmonic acid related genes, however GID1, exhibited the opposite expression pattern.

cis-Regulatory Element Analysis of Potassium Deficiency Responsive Genes
When analysing the expression of a group of genes that respond similarly to a particular condition, it is anticipated that these genes might have some common features or elements driving their expression. The promoters of co-expressed genes might share some common regulatory elements or be regulated by a common set of transcription factors. Thus, the identification of shared cis-regulatory elements in the promoter regions of genes co-expressed in response to altered levels of potassium availability could provide new insights into the transcriptional regulatory networks that are triggered by potassium deficiency in plants. To identify common cis-regulatory elements, 1 kb regions directly upstream of potassium-responsive genes in Arabidopsis and rice were analysed in PLACE database [38,39]. Genes that exhibited high levels of expression under potassium deficiency in this study and earlier studies on Arabidopsis including CIPK9 [40], HAK5 [41], peroxidase, glycosyl hydrolase, AP2 domain-containing transcription factor, rice alpha-amylase, glycosyl hydrolase, peroxidase, potassium transporter HAK1, GST, dehydration stress-induced protein, and ethylene-responsive transcription factor were compared. A total of 31 common cis-regulatory elements were found, including ABRELATERD1 (ACGTG), ARR1AT (NGATT), CAAT-BOX1 (CAAT), DOFCOREZM (AAG), GATABOX (GATA), GT1CONSENSUS (GRWAAW), POLLEN1LELAT52 (AGAAA), TAAAGSTKST1 (TAAAG), WRKY71OS (TGAC) and many others that were prominently present in the 1 kb putative regulatory regions of the aforementioned genes. A list of all of the 31 cis-regulatory elements with their sequences is provided in Table S7.

Comparative Analysis of Rice and Arabidopsis in Potassium and Other Nutrient Deprivation
Transcriptome data from this study were compared to earlier published results on the transcriptome response to potassium deficiency in Arabidopsis and rice, particularly an earlier transcriptomic profile of 2 week old rice roots grown under potassiumdeficient conditions [37]. However, the transcriptome profile of potassium-deficient plants is comparatively far better characterized in Arabidopsis, where both two week old shoot and root tissues were analyzed under potassium deprivation [5,42]. We compared the dataset; we obtained from seedlings grown under potassiumdeficient conditions under these two published transcriptome datasets. In earlier reports, far fewer genes were reported to be differentially regulated with a 2-fold or greater expression change in response to potassium deficiency in both rice (722 in rice whole seedling [this study] versus 356 in rice root [37]) and Arabidopsis (119 in shoot, 299 in root [5]) (Table 3). In an attempt to identify the similarities and differences among potassium-responsive genes in rice and Arabidopsis, we categorised and compared 356 DEG from rice root [37] with 722 DEG from whole rice seedling (this study) and 418 genes from Arabidopsis shoot and root transcriptome [5] ( Figure S1). This comparative analysis revealed that the majority of DEG was related to metabolism and signal transduction in response to potassium deficiency in both rice and Arabidopsis ( Figure S1). The percentage of genes related to defense and cell wall were significantly greater in number in Arabidopsis than rice. In whole rice seedlings, the number of differentially expressed transcription factor genes were greater (this study, Figure S1) compared to the roots of potassium-deficient rice seedlings [37]. The present transcriptome profile of potassium-deficient rice seedlings was further compared to transcriptome profiles of rice under other nutrient deprivation (i.e. nitrogen, phosphorus, and iron) ( Figure S2). From the above studies, it was inferred that genes differentially expressed in response to deprivation of distinct nutrients were associated with various metabolic processes and stress adapatations and were classified as transporters, transcription factors, and signal transduction components. This suggest that there might be a functional overalap of these genes in responses to deprivation of nutrients other than potassium in both rice and Arabidopsis.

Relationship of Potassium Ion Response Related Genes with Abiotic Stresses
We analysed potassium responsive DEGs (mainly related to metabolism, signal transduction and stress) under abiotic stresses like cold, drought, heat and salt stress, using Genevestigator database for their transcript status. Most of the genes showed differential expression pattern under these stresses. Metabolism related carboxyvinyl-carboxyphosphonate phosphorylmutase (LO-C_Os12g08760), heparanase (LOC_Os06g08090), UDP glucoronosyl (LOC_Os01g08440) showed downregulation in all stresses except heat stress. But phosphoglycerate kinase (LO-C_Os06g45710) and beta-amylase (LOC_Os03g04770) were upregulated in all stress conditions ( Figure S3A). Genes related to signal transduction, such as OsCML14 (LOC_Os05g50180) and BRI I, (LOC_Os11g31530) showed an upregulation in both cold and drought stresses whereas no differential expression was found in heat and salt stresses. MAP3K4 (LOC_Os01g50370) expression was elevated in response to above mentioned four stresses. In contrast, EF-hand containing protein (LO-C_Os09g28510) was showing downregulation in stressed state ( Figure S3B). Early responsive to dehydration protein (LO-C_Os12g43720), oxidoreductase (LOC_Os11g07930), HSP (LO-C_Os2g54140), HSF (LOC_Os02g35960) genes were downregulated in cold and heat stresses while upregulated under drought stress. Dehydrin (LOC_Os02g44870) showed elevated expression under cold, drought, heat and salt stresses whereas 2-oxoglutarate (LOC_Os10g40960) and oxidoreductase (LOC_Os02g53510) were downregulated under these conditions ( Figure S3C).

Discussion
Nutritional deficiencies alter the physiology and metabolism of plants during both short-term and long-term periods of growth. In response to nutrient deficiency, plants activate a plethora of genes and gene networks encoding a large number of proteins involved in acclimation and adaptation. Significant research has contributed to our understanding of the uptake, distribution and homeostasis of potassium ions in Arabidopsis by employing advanced tools for genetics, biochemistry, cell biology, and physiology. Our knowledge of adjustment and adaptation to short-term and long-term nutrient deficiency is still insufficient to establish a basic model of how plants manage to thrive physiologically despite nutrient deficiencies in soil. As mentioned in the introduction, research related to potassium nutrition in plants has drawn the attention of several leading research groups who seek to enhance our understanding of the mechanisms of plant tolerance to potassium deficiency. In this study, we have contributed to this effort by performing detailed transcriptomic profiling of whole rice seedlings exposed to short-term (5 days) potassium deprivation.
It is quite possible that different plant species might have different mechanism to absorb, transport, distribute, and utilize potassium [4,6,7,[43][44][45]. Biomass is one important parameter to measure plant growth under low nutrient conditions, and it is a significant criterion to judge the crucial role played by particular nutrient in plant growth and development.
In order to complement our transcriptomic analysis of rice seedlings, we also performed a physiological experiment to determine the growth symptoms of potassium deficiency. The experimental results indicate that rice grew significantly more poorly under low potassium conditions than it did under higher concentrations of potassium that we believe are more reflective of optimal growth conditions for rice. This effect is clearly seen in the qualitative and quantitative assessment of seedling biomass under both conditions ( Figure 1A, B, C). At the molecular level, the content of potassium ions in plants grown in low potassium media was also appreciably lower than the potassium content of plants grown in media with sufficient levels of potassium, which demonstrates that rice seedlings are indeed able to take up higher levels of potassium when external levels are sufficiently higher and further establishes a direct link between potassium supply and plant growth ( Figure 1D). Similar results were previosuly described in Arabidopsis plants starved of potassium [5]. Analysis of the growth of Arabidopsis under long-term potassium starvation (14 days) by Armengaud et al. (2004) [5] revealed chlorosis (bleaching) of older leaves, as well as a reduction in lateral root growth of seedlings. In contrast to these reported results in Arabidopsis, we observed an obvious reduction in the growth of rice seedlings by the fifth day of growth on low potassium containing media without any significant signs of chlorosis or bleaching of seedling leaves. We also observed a very similar, severe impairment of growth in rice seedlings maintained on potassium-deficient media for longer period of time (i.e. 14 days; data not shown). After observing growth differences of rice seedlings exposed to short-term or longterm potassium deprivation, we decided to investigate gene expression profiles under short-term (5 days) potassium starvation.
During the preparation of this manuscript, we discovered a report on the transcriptome profile of rice roots deprived of potassium [37]. The major difference between this study and our study is that we chose to use whole seedlings rather than just roots. Moreover, this study is more comprehensive due to our transcriptional profiling of seedlings resupplied with potassium after short-term potassium deprivation. These growth conditions closely mimicked a previous study by Armengaud et al. 2004 [5]. The usage of higher level of stringency (i.e. 2-fold or greater change) and false discovery rate (FDR) in the present study also strengthens the reliability of our data on differential gene Table 3. Differentially expressed genes from present microarray data and data obtained from 2 week K + -starved Arabidopsis shoot-root (Armengaud et al. 2004), roots of rice plants after 6  expression. Despite several differences, some of our findings are consistent with differentially expressing genes previously reported by Ma et al. 2012 [37]. The most significant contribution of the present work is the identification of genes implicated in the growth and development of rice seedling exposed to the short-term stress of potassium deficiency and the resupply of potassium to deprived seedlings. This study provides a platform for the investigation of the mechanisms of regulatory responses to low potassium conditions. From our microarray analysis, a large number of genes (722) were found to be differentially expressed during potassium deprivation and resupply (1867 DEG). Only 307 of the affected genes, however, exhibited reversal of gene expression (i.e. restoration of comparable expression levels as seen during optimal growth conditions) upon resupply of potassium to the deficient media ( Figure 2B), therefore this study clearly identifies a subset of the differentially expressed genes that are specifically expressed only during potassium deprivation. This difference in expression is also seen in the results of our PCA analysis, where the three biological replicates of the potassium resupply treatment were more similar to those grown with sufficient levels of potassium (Figure 2A). This recognizable difference in gene expression between the plants during potassium deficient conditions and following resupply highlights the flexibility of biological pathways to adapt quickly as per changes in the environment. Through this study, we reveal adaptive responses of plants to fluctuations in potassium ion concentration. As previously reported, genes involved in various processes like primary and secondary metabolism, ion transport, signal transduction, and hormone signaling were linked to potassium deficiency [5,[44][45][46][47]. We have also identified differentially expressed genes related to carbohydrate metabolism, carboxyvinyl-carboxyphosphonate phosphorylmutase, glyoxalase, lactate/malate dehydrogenase, phosphoglycerate kinase, UDPglucosyl transferase and many other enzymes and biological processes ( Table 2). Many genes related to carbohydrate metabolism have also been reported to play important roles in response to phosphorus deficiency in rice roots and potassium deficiency [46,48]. One of the genes identified in this study, UDPglucosyl transferase, was also found to be influenced by iron toxicity in rice [49]. Similarly, alpha amylase transcript that shows the greatest magnitude of upregulation in response to potassium deficiency was also reported to be involved in starch hydrolysis [50].
Our microarray data suggests a number of proteins involved in diverse aspects of carbohydrate metabolism were differentially regulated and thereby indicates that potassium deficiency has widespread effects on carbohydrate metabolism. Carbohydrate metabolism directly impacts plant biomass and photosynthetic yield and, based on our findings, could be a significant component of the stress imposed by potassium deprivation.
Lipases are enzymes involved in the cleavage of lipids and are crucial components of lipid metabolism responsible for the mobilization of lipids during stress conditions. In this study, several lipases were identified that exhibited downregulation during potassium deprivation. In Arabidopsis, Armengaud et al. (2009) [5] reported the involvement of lipases and their similar differential expression in potassium deficiency treatments.
Based on earlier physiological analysis of growth symptoms under potassium deficiency, photosynthesis was identified as one of the major processes directly affected by potassium deficiency [51][52][53]. Photosynthesis serves as the major energy source of plants, and its efficiency is reduced under nutrient deficiencies [54]. The efflux of potassium ion from guard cells causes decreased stomatal conductance and thereby reduces rates of photosynthesis as a tradeoff for reduced water loss [51,[55][56][57][58]. In our differential gene expression analysis, we have identified several genes involved in photosynthesis. In particular, chlorophyll A-B binding protein and magnesium-protoporphyrin IX monomethyl ester cyclase were downregulated in KM conditions and strongly upregulated upon KR conditions. Similarly, several genes involved in jasmonic acid (JA), auxin, ethylene, and gibberellin signaling pathways are known to be directly influenced by potassium deficiency in plants [5,37,59].
With extensive research in the field of stress biology in plants, it is widely appreciated that many of these stresses are inter-related and produce similar secondary messenger molecules such as reactive oxygen species (ROS), calcium (Ca 2+ ), and others [32,54,60]. Research in Arabidopsis pertaining to deprivation of multiple nutrients (N, S, P, K, Fe and others) suggests the activation of signaling pathways implicated in other abiotic stresses such as drought or toxic ion exposure [32,61]. Presumably due to the interconnected nature of diverse plant stresses, plants grown under low potassium conditions have also been shown to be particularly susceptibile to drought and cold stress [9,34,36]. Moreover, many of the abiotic or biotic stresses affect metabolism and hence affect growth [62]. Different cellular processes like growth, photosynthesis, carbohydrate and lipid metabolism, osmotic homeostasis, protein synthesis are affected by both biotic and abiotic stresses [63][64][65][66].
Reactive oxygen species (ROS) are a byproduct of normal cell metabolism in plants. Accumulation of ROS is damaging to cells, however ROS also act as important signaling molecules. The balance between production and detoxification of ROS is tightly regulated by many proteins and enzymes [67][68][69]. ROS are generated as a signal of several biotic and abiotic stress conditions [68,[70][71][72]. ROS level were also elevated in response to deficiency of several nutrients (N, S, P, Fe), and ROS are a known signal of potassium starvation in roots [31,32,45,73]. Under potassiumdeficient conditions, many ROS related genes such as oxidoreductase and glutathione S-transferase (GST) are involved in detoxifying the ROS to maintain homeostasis of these molecules in the cytosol. We observed that transcripts encoding these enzymes were upregulated in response to potassium deprivation and subsequently downregulated following resupply of potassium.
Recently in the field of calcium signaling in plants, calcium has often been placed downstream of ROS [26,31,32,40,74,75]. In previous reports in Arabidopsis, calcium-binding proteins were described to be differentially regulated in response to potassium starvation conditions [45]. Calcium sensors such as calmodulin (CaM), calcineurin B-like proteins (CBLs), CBL interacting protein kinases (CIPKs), and Ca 2+ -dependent protein kinases (CDPKs), were differentially expressed in potassium-deficient conditions [5]. In the rice root transcriptome, OsCIPK29 was found to be downregulated in response to potassium deficiency [37]. Three other CIPKs (OsCIPK7, OsCIPK8, and OsCIPK29) were found to be downregulated under potassium starvation. Their expression increased after resupply of potassium in this study, which suggests that these kinases might be acting as negative regulators of potassium nutrition and signaling in rice. One recently discovered, important aspect of potassium uptake and nutrition that has been elucidated in Arabidopsis is regulation of the high affinity K + channel AKT1 by physical interactions with the calcium-binding proteins CBL1 and CBL9 and the Ser/Thr kinase CIPK23 in the root to regulate uptake of potassium [24,25]. The downregulation of two PP2C protein phosphatases in the transcriptomic data suggests that these phosphatases might be acting as counteracting molecules in protein phosphorylation pathways that signal potassium-deficient conditions. Previously, it was reported that several members of the PP2C family are involved in regulating responses to multiple abiotic stresses as well as developmental processes in rice [76]. Some MAP kinases were also differentially regulated in potassium-deficient conditions, which suggest other signaling components are involved in the regulation of potassium starvation responses. MAPKs have also been implicated in both biotic and abiotic stress responses in addition to potassium starvation, development, cell proliferation, and hormone physiology [47,61,77,78]. The involvement of phytohormones like JA and auxin in potassium-deficient conditions was shown earlier by identification of differentially regulated genes in the metabolism or signaling of these phytohormones [5,37]. Surprisingly, we identified only a few candidate genes involved in JA and auxin metabolism and signaling during potassium-deficient conditions in rice seedling. It is quite possible that the altered expression of these signal transduction, stress and hormone-related genes might also play roles in enabling plants to adapt under low potassium stress conditions in rice.
Genome-wide analysis has led to the identification of a large family of genes encoding K + transporters/channels in both Arabidopsis and rice [79][80][81][82][83]. Moreover, several detailed gene expression profiles by northern blots and RT-PCR analyses identified multiple K + transporters that are differentially regulated under potassium starvation conditions [19,61,84]. However, in this transcriptomic analysis, we could not identify any K + transporters with at least two fold change in gene expression, in contrast to the findings of Ma et al 2012 [37]. Similar to our findings, two previous reports on transcriptomic responses to potassium starvation in Arabidopsis did not identify many K + transporters with significant changes in expression [5,42]. Upon relaxing our analysis stringency to 1.5 fold with FDR, p,0.05, three transporters including HAK1 (LOC_Os04g32920) and KAT1 (LOC_Os01g11250) were identified. OsHAK1 is a member of the KUP/HAK/KT family of putative high-affinity K + transporters that have been known to be involved in potassium uptake under low potassium conditions [37,85]. KAT1 is a shakertype K + channel that is specifically expressed in guard cells and mediates potassium fluxes for turgor-dependent regulation of the stomatal aperture [86,87]. As reported earlier, HKT class I transporters conduct Na + uptake under KM condition in rice. In contrast, class II HKT transporters showed higher permeability for K + ions. [88][89][90][91][92][93]. In transcriptomic analysis, the expression of sulfate transporters such as SULTR1;2 was found to be altered under potassium-deficient condition in Arabidopsis [5]. Approximately 30 DEG annotated as transporters that conduct molecules other than K + were identified in this study, and these are believed to be involved in the transport of metals, lipids, sulphate, and other molecules. We identified representatives of ABC transporters, aquaporins, MATE efflux transpirters, and AAA-ATPase-type transporters. The differential expression of these transporters in potassium-deficient conditions might be related to the homeostasis of other molecules to maintain proper osmoticum and neutralize the charges across the membranes and in the cytosol [19,32,94].
The identification of fewer K + transporters led to an argument for the possible existence of post-translational regulatory mechanisms of K + transporters that direct uptake and distribution of potassium from the root to shoot. In fact, this mode of regulation has been recently discovered in Arabidopsis, where a high affinity K + channel was regulated by a calcium-mediated phosphorylation-dephosphorylation network in roots [24,25,95]. However, no such phosphorylation-dephosphorylation signaling network has been identifed yet in rice. We speculate that similar posttranslational modes of regulation are likely present in rice for regulation of K + transporters during low potassium conditons.
In the present transcriptomic analysis, genes for a large number of transcription factors were differentially regulated in KM condition and their expression returned to baseline levels upon resupply of potassium. The maximum effect was shown by members of the zinc-finger transcription factor family, followed by MYB family transcription factors (Table 2). Zinc-finger family and basic helix-loop-helix (BHLH038, BHLH039, BHLH100, and BHLH101) transcription factors were previously reported to be involved in responses to iron deficiency in Arabidopsis [96]. Multiple reports in plants suggest that MYB family transcription factors regulate a wide range of processes in plants, including biotic and abiotic stress responses, plant growth and development, and secondary metabolite production [97][98][99][100][101]. In addition to potassium and iron, this family is also involved in phosphate starvation signaling in both vascular plants and unicellular algae [102].
Regulatory cis elements present in promoter regions are one of the salient aspects of gene regulation and often provide a mechanistic explanation for coexpression of a group of coexpressed genes. Several common cis-regulatory elements have been identified in multiple genes regulated similarly at transcriptional level [103]. Upon analysis of cis-regulatory elements in some of the highly differentially expressed genes in the present study, we found that they contained previously identifed drought responsive, sugar responsive, and light responsive regulatory cis-elements including ABRELATERD1, ACGTABOX, GT1CONSENSUS, GATA-BOX and others [104][105][106][107]. Interestingly, an earlier study concluded that ABREs functions Ca 2+ -responsive cis-elements in response to calcium transients induced by biotic and abiotic stresses [108][109][110]. Also the expressions of many potassium transporters and channels (KCO) have been shown to be activated by cytosolic Ca 2+ [111]. The TAAAGSTKST1 element was found in promoter of the StKST1 gene encoding a K + channel, and this gene is a known target site for trans-acting StDof1 protein controlling guard cell-specific gene expression [112]. Thus, the presence of these cis-regulatory elements in the promoter sequences of genes involved in nutrient stress might be associated with several abiotic and biotic stresses that lead to changes in the expression of signal transduction components and might also regulate signaling pathway crosstalk to enable plants to acclimatize to a wide range of conditions. Despite the presence of several common cis-regulatory elements in the promoters of highly differentially expressed genes under potassium deprivation, we could not identify a motif or cis-element specific for potassium responsiveness.
The comparative transcriptome analysis of rice and Arabidopsis ( Figure S1, S2) suggests a common regulatory mechanism for genes whose expression is highly dependent on potassium availability. Such genes are implicated in a variety of cellular processes including metabolism, molecular transport, transcriptional regulation and signal transduction. We believe the concerted interplay of this diverse suite of genes coordinates plant adaptation to potassium starvation.
Cold, drought, and high salinity are common stress conditions that adversely affect plant growth and crop productivity. Due to abiotic stresses, the plant metabolism and homeostasis is also affetced. In order to achieve adjustment of metabolic pathways, different signaling networks also worked together [113]. Plant adaptation to environmental stresses is dependent upon the activation of various cascades of molecular networks involved in stress perception, signal transduction, and the expression of specific stress-related genes and metabolites [114]. Genes found in potassium deprivation condition were also involved in plant stress management. Carboxyvinyl-carboxyphosphonate phosphor-ylmutase showed 62 fold downregulation in potassium deficient condition, was also downregulated in different stress treatment. Also potassium responsive genes such as heparanase, UDP glucoronosyl, phosphoglycerate kinase and beta-amylase were also directly involved in plant stress management [115,116,117,118]. Various reports revealed that genes involved in signal transduction like CML14, BRI I, MAPK were also responsive to several environmental stresses [119,120,121].

Conclusions
In the present study, comprehensive gene expression analysis was performed using Affymetrix rice genome arrays during shortterm potassium deficiency and resupply of potassium to deprived seedlings. A large number of genes were found to be differentially expressed. Taking into account the many potential consequences of the gene expression changes we observed, our study emphasizes the major metabolic adjustments associated with adaptation to the stress of nutrient deficiency. Potassium nutrition impacts a number of metabolic, signaling and physiological processes; together these processes coordinate plant adaptation to low potassium conditions.
The overlap of low potassium responses with responses to biotic and abiotic stresses attests to the complex physiological changes triggered by nutrient starvation conditions. Under low potassium conditions, a large array of genes are activated, starting from sensing of nutrient deficiency in the soil to the transducing of this signal by activation of a plethora of signaling networks to, ultimately, the regulation of downstream proteins including transporters, enzymes, and transcription factors. We expect that the transcriptional changes in several diverse physiological and metabolic pathways are critical for adaptive responses to plant potassium deprivation conditon. Besides transcriptional regulation, post-translational modulation could certainly constitute another important mode of regulation in response to low potassium conditions. Despite the complexity of responses to potassium deficiency, we have made a hypothetical model to summarize the results of this study ( Figure 7). As described in (Figure 7), low potassium influences metabolic processes and impacts genes related to different stress responses. Although, both processes are interconnected by various signaling pathways, the differential expression of individual components might affect the net signaling output of these highly interrelated processes. The changes in expression patterns of these components might govern the acclimation of the plant to such conditions. Further verification of these differentially expressed genes in planta by genetic and biochemical approaches will provide insights into the mechanistic regulation of these gene (s) during low potassium stress conditions, which ultimately can be translated to develop plants that can grow on soil with relatively low levels of potassium without compromising yield and productivity.

Plant Material and Growth Conditions
Seeds of indica rice (Oryza sativa subsp. indica cv. IR64) were surface sterilized and grown hydroponically in culture solution prepared as described [122]. The solution contained 1.14 mM NH 4 SO 4 ) and other one half to medium without potassium (KM, 0 mM). After 5 days, potassium (0.409 mM K 2 SO 4 ) was resupplied to the plants grown in KM medium for 6 h. The nutrient medium was changed at the interval of every 2 days. The seedlings were then harvested and immediately frozen in liquid nitrogen and stored at 280uC until total RNA isolation.

Measurement of Total Ion Content
Seedlings treated with the above mentioned conditions of K + concentration were collected and dried in oven at 80 degree for 24-48 h. To 100 mg of the dried tissue, 1 ml of spectroscopic grade concentrated HNO 3 was added for the solubilisation. Samples were then incubated at room temperature (RT) for 2 days to ensure complete solubilisation of the tissue. The solubilized tissue was digested in the digester at 80u celcius temperature to vaporize almost 95% of the acid. The completely digested tissue was then transferred to a fresh micro-centrifuge tube and the volume was adjusted to a total of 1 ml using triple distilled milliQ water. The content was centrifuged at 10,000 g for 5 min and the supernatant was filtered using a 0.2 micron filter (Mdi, India). Three biological replicates were used for the statistical analysis. The concentrations of Na + and K + ions in seedling tissues were measured using a single-channel flame photometer (Digital Flame Analyzer model 2655-00; Cole-Parmer). In each analytical batch, reagent blanks and spiked samples were included in the acid digestion to assess the accuracy of the chemical analysis. The recovery of spike was 94-95% (n = 6). 100 mM and 40 mM solutions of NaCl and KCl were used as reference standards during analysis. Three biological replicates were analysed for each treatment.

Measurement of Fresh and Dry Weight
A total of 15 rice seedlings were grown in normal (KP) and potassium deficient (KM) condition followed by weighing (fresh weight) separately and dried at 80uC for 48 h, and then again weighed (dry weight). DW/FW ratio was calculated and t -test (*P,0.05, **P,0.01) was used to analyze the statistical significance. Three biological replicates were used for fresh and dry weight measurement.

RNA Extraction and Microarray Experiments
Total RNA was isolated from whole seedlings, except seeds, using the TRIzol method (Invitrogen Inc., USA), and RNA was further purified by Nucleospin RNA clean-up columns (Macherey-Nagel, Germany). The total RNA samples were quantified using nanodrop (ND-1000 Spectrophotometer (NanoDrop Technologies, USA ) and then analysed on denaturing MOPS agarose gel (2%) as well as by using Agilent Bioanalyser (Agilent Technologies, CA) for determining integrity and quality of RNA. After verifying the quality, 500 nanograms of RNA with 260:280 ratios of 1.9-2.0 and 260:230 ratios more than 2.0 was used for cDNA synthesis. Affymetrix GeneChipH Rice Genome Arrays representing 57,381 probe sets have been used to study the transcriptome profiles of potassium responsive genes in potassium deficient stress in rice. Labeling and hybridizations were carried out according to Affymetrix manual for one-cycle target labeling (Affymetrix, Santa Clara, CA). Hybridization was performed in a GeneChipH Hybridization oven 640 for 16 hours at 45uC and 60 rpm. GeneChips were washed using the fluidics protocol Eu-kGE_WS2V5_450 in Affymetrix fluidic station model 450 and the chips were scanned using the Gene-ChipH Scanner 3000 (Affymetrix, Santa Clara, CA). Three biological replicates were used for each treatment.

Microarray Data Analysis
The CEL files created by genechip operating software (GCOS) were imported into array assist (version 5.0) software (Stratagene, La Jolla, CA) for further analysis. GCRMA algorithm and log transformation were used for the normalization of microarray data. Student's t-test was performed on log transformed value with the Benjamini-Hochberg correction to identify the differentially expressed genes. Benjamini-Hochberg correction ensured that transcripts showing significant changes were only analysed. Genes with fold change $2 and p,0.05 were considered to be significantly expressed. Oligonucleotide sequences of the probes were then mapped to Osa1 Rice Genome Annotation Project release 6 (RGAP: http://rice.plantbiology.msu.edu/) to annotate the differentially expressed genes. For hierarchical clustering, Euclidean distance metric and Wards linkage rule was used. Additional K-means clustering was carried out to identify the expression patterns of all differentially expressed genes in potassium-deficient conditions.

Real-time PCR Analysis for the Validation of Microarray Data
Nineteen significantly differentially expressed genes under potassium deficiency were selected for quantitative real time PCR. The same RNA samples that were used for microarrays were also used for real time PCR. For all the selected genes primers were designed mainly from 39 end of the gene using PRIMER EXPRESS version 2.0 (PE Applied Biosystems, USA) with default parameters. The primer sequences are listed in Table  S4. Specificity of the primers was checked using BLAST tool of NCBI. First strand cDNA was prepared from 2 mg of DNase treated total RNA in 50 ml of reaction volume using high-capacity cDNA Archive kit (Applied Biosystems, USA). One ml of the first strand cDNA reaction was used for quantitative real time PCR. KAPA SYBR FAST Master Mix (KAPABIOSYSTEMS, USA) was used to determine the expression level for the genes in ABI Prism 7000 Sequence detection System (Applied Biosystems, This condition alters plant metabolic as well as signaling responses. During K + deprivation conditons, the initial reaction of plant is to sense the K + deficiency followed by activation of several genes and gene networks, which might be working in tandem. As depected several genes encoding metabolic enzymes, protein involved in stress signaling, transporters, transcription factors could be the final target to bring about a physiological response such as tolerance or adaptation to low-K + condition. Putative factors of K + deficiency and adaptive responses are shown in boxes. Genes of related components are written below the boxes. doi:10.1371/journal.pone.0070321.g007 USA). Three biological and three technical replicates were taken for each treatment, and standard deviation and standard error were calculated. As an endogenous control, Actin was used for the normalization of Ct value obtained and the relative expression values were calculated by DDCt method. Figure S1 Comparison of functional category between rice and Arabidopsis in responses to K + deficiency.