Genome-wide co-expression analysis predicts protein kinases as important regulators of phosphate deficiency-induced root hair remodeling in Arabidopsis
© Lan et al.; licensee BioMed Central Ltd. 2013
Received: 4 December 2012
Accepted: 20 March 2013
Published: 1 April 2013
Phosphorus (P) is one of the essential but often limiting elements for plants. Based on transcriptional profiling we reported previously that more than 3,000 genes are differentially expressed between phosphate (Pi)-deficient and Pi-sufficient Arabidopsis roots (MCP 11(11):1156–1166, 2012). The current study extends these findings by focusing on the analysis of genes that encode protein kinases (PK) and phosphatases (PP) by mining PK and PP genes that were differentially expressed in response to Pi deficiency.
Subsets of 1,118 and 205 annotated PK and PP genes were mined on the basis of the TAIR10 release of the Arabidopsis genome. Analysis of RNA-seq data showed that 92 PK and 19 PP genes were not detected in roots (zero reads in three biological repeats); 96 PK and 10 PP showed low abundance (≤ 10 reads). Gene ontology analysis revealed that the 188 PK genes with no or low expression level in Arabidopsis roots are mainly involved in pollen recognition, pollen tube growth or other processes not relevant for root hair formation. More than 50% of the cysteine-rich RLK (receptor-like protein kinase) subfamily genes belong to this group. Among the 29 PP genes with no or low expression level, purple acid phosphatases, haloacid dehalogenase-like hydrolases, and PP2C genes with functions in the dephosphorylation of RNA polymerase II C-terminal domain and mRNA capping were enriched. Subsets of 173 PK and 35 PP genes were differentially expressed under Pi-deficient conditions. Putative functional modules (clusters) of these PK and PP genes were constructed based on co-expression analysis using the MACCU toolbox. A co-expression network comprising 65 known or annotated PK and PP genes (60 PK and 5 PP genes, respectively) was subdivided into several highly co-expressed gene sub-clusters. The largest sub-cluster was composed of 22 genes, most of which have been assigned to the RLK superfamily and were associated with cell wall metabolism, pollen tube and/or root hair development and growth.
We here provide comprehensive ‘digital’ transcriptional information on PK and PP genes in Arabidopsis roots. The co-expression network derived from our data mining approach sets the stage for follow-up experimentation that helps to complete our understanding of the post-translational regulation of Pi deficiency-induced changes in root hair morphogenesis.
Phosphorous, mainly taken up as phosphate (Pi) by plants, is an essential micronutrient involved in signaling, metabolism and photosynthesis. The bioavailability of Pi is often very low due to its tendency to form complexes with soil cations. In agricultural systems, Pi deficiency is a major cause of severe yield losses in crops and poor quality of edible plant parts. Low Pi availability is often corrected by the application of large quantities of fertilizers, which is associated with environmental pollution and substantial costs. Understanding how plants adapt to low Pi availability is thus mandatory to develop Pi-efficient germplasms. To cope with low Pi availability, plants have evolved an array of adaptive processes aimed at improving Pi uptake and re-mobilization, comprising the acquisition and redistribution of Pi, alterations in developmental programs, and metabolic networks . Proteomic [2–4] and transcriptomic [5–11] profiling studies have uncovered several robustly changed processes in Pi-deficient plants, including the remodeling of lipid metabolism, changes in glycolytic carbon flux, alterations in root development, and rewired signaling pathways [12–14].
The mechanisms underlying the maintenance and recalibration of cellular Pi homeostasis are complex. The Myb-type transcription factor PHOSPHATE STARVATION RESPONSE1 (PHR1) is a central conserved regulator that controls a subset of Pi deficiency genes by binding to an imperfect palindromic sequence motif [15, 16]. Consistent with a critical regulatory role of PHR1 in Pi homeostasis, overexpression of PHR1 led to increased Pi accumulation . The activity of PHR1 is controlled by the SUMO E3 protein ligase SIZ1 , representing the most upstream component of the Pi deficiency signaling cascade identified so far. Another subset of Pi-responsive genes is regulated by the E2 ubiquitin conjugase PHOSPHATE2 (PHO2). A connection between these two central switches is established by MicroRNA399 (MiR399), which systemically controls PHO2 through transcript cleavage [19, 20]. MiR399 itself is strongly induced by Pi deficiency . The sensor for Pi remains to be discovered.
Besides the involvement of protein ubiquitination [18, 21], other posttranslational processes potentially involved in the Pi deficiency response have not been thoroughly investigated. An estimated one-third of all eukaryotic proteins undergoes reversible phosphorylation via protein kinases (PK) and phosphatases (PP), demonstrating the importance of this process. Modifications of protein with phosphate can affect protein structure, activity, localization, interaction, and stability, thereby regulating crucial processes such as metabolism and development. Several hundred genes encoding PKs and PPs were found to be differentially expressed upon Pi deficiency by transcriptional profiling of roots from Pi-deficient plants , suggesting that alterations in protein phosphorylation patterns induced by Pi deficiency are critical in the control of Pi homeostasis. For example, under Pi-limiting conditions the high-affinity phosphate transporterPHT1;1 was found to be induced and newly–synthesized PHT1;1 protein was phosphorylated by an yet unknown PK at the C-terminal 514 amino acid Ser, which is required for the precise localization of PHT1;1 to the plasma membrane .
Transcriptome analysis alone, however, is insufficient for defining potential roles of differentially expressed PKs and PPs genes in Pi homeostasis. Functional characterization of these genes by reverse genetic approaches such as increasing or decreasing their transcript level (by T-DNA insertion and/or RNAi) is required to elucidate their biological functions. However, individually assaying hundreds to thousands of differentially expressed genes without any selection filter would be extremely laborious. Systemic cluster analyses provide a means to filter and select genes of interest for the biological question addressed. Genes showing similar expression pattern under diverse conditions often have correlative functions , and the biological processes in which genes with unknown functions are involved can be predicted based on co-expression data (‘guilt by association’) [24, 25].
In the present study, the global expression of PK and PP genes in Arabidopsis roots was analyzed in order to gain insights into the regulation of the interplay of transcriptional and post-translational responses to Pi deficiency. By mining public databases, PK and PP genes that are differentially expressed upon Pi starvation were clustered into groups of closely correlated modules based on their co-expression under various sets of experimental conditions. Using this approach, we discovered several potentially critical regulatory PKs with roles in root hair development and growth under Pi deficiency.
Results and discussion
Digital information on the expression of protein kinase and phosphatase genes in Arabidopsis roots
Protein kinases and phosphoylases play key roles in regulating nearly all aspects of cellular processes. However, due to the use of microarray probe sets that have significant cross hybridization potential and are unable to distinguish highly similar genes of this subfamily, transcriptional information on the expression of PK and PP genes is incomplete in Arabidopsis. In addition, tissue-specific gene expression information of PKs and PPs is lacking. The RNA-seq technology has proven to provide precise digital information on gene expression, and is able to discriminate genes of high sequence identity . Using this technology, we previously examined global gene expression changes upon Pi deficiency in Arabidopsis roots . Focusing on the expression of PKs and PPs, we mined and re-analyzed our RNA-seq data set (NCBI: SRA050356.1). Based on the PK and PP gene families annotated in the TAIR10 release of Arabidopsis genome, a total of 1,118 PK (GO: 0004672, Additional file 1) and 205 PP (GO: 0004721, Additional file 2) genes was retrieved and compared with the RNA-seq data. We defined a gene as not being expressed if the unique read number was zero in all three biological repeats under normal (Pi-replete) conditions; low abundance of a transcript was defined by a unique read number ≤ 10 in either of the three biological repeats. A gene was defined as being highly expressed in Arabidopsis roots when the read number was higher than 2,000 in either of the three biological repeats. On the basis of these criteria, 92 PK genes were not detected in Arabidopsis roots, transcripts of 96 PK genes were low abundant; and 57 PK genes were highly expressed (Figure 1A and Additional file 3). For the 1,118 PK genes, 432 cognate proteins were identified with at least one uniquely matched peptide ( and Additional file 4). Generally, PK proteins were more likely detected when the associated transcript was highly abundant. However, proteins were also detected from about 20% of the genes with low abundant or absent transcripts (Figure 1C and Additional file 4) confirming the observation that gene expression is not always correlated with protein abundance .
Applying the same criteria used for the PK genes, subsets of 19 and 10 PP genes were not detected or lowly expressed in Arabidopsis roots, respectively (Additional file 7). The trend of PP protein expression was similar to that of PK proteins (Figure 1D and Additional file 4). Among these two subsets, purple acid phosphatases (PAP), haloacid dehalogenase-like hydrolases (HAD), and PP2C genes were enriched, with gene products mainly localized in nucleus and functioning in the biological processes ‘dephosphorylation of RNA polymerase II C-terminal domain’ and ‘mRNA capping’ (Additional file 8). Products of the nine highly expressed PP genes (Additional file 7) were mainly involved in the formation of PP-1 and PP2A complexes associated with the biological processes ‘embryonic root morphogenesis’ and ‘phosphate ion homeostasis’ (Additional file 9). The lack of a transcriptional response of PP genes to Pi deficiency suggest that under normal conditions in Arabidopsis roots cellular Pi homeostasis is regulated by PPs, probably by protein de-phosphorylation. Uncovering the substrates regulated by these PP genes would strongly facilitate our understanding of cellular Pi homeostasis.
Identification of Pi-responsive genes encoding protein kinase and phosphatase in Arabidopsis roots
Although approximately one thousand genes have been identified as being Pi-responsive by using the full-genome Affymetric ATH1 gene chip [5–11], the technical limitations of the microarray technology renders a precise estimation of the changes in the expression of Pi-responsive genes difficult. Using RNA-seq, we defined a total of 3,106 genes as differentially expressed between Pi-sufficient and Pi-deficient Arabidopsis roots . Subsets of 173 PK and 35 PP genes, comprising diverse subfamilies were differentially expressed under Pi-deficient conditions (Additional file 10), some of which have been listed as Pi-responsive genes in earlier microarray studies [5–11]. Among the differentially expressed PK genes, members of the RLK superfamily genes, especially from the LRR-RLK subfamily, and genes from the CDPK-SnRK superfamily were enriched, while members of the PAP and PP-2C families were predominant in the differentially expressed genes encoding PPs. In total, 205 PK and PP genes were finally mined for further analysis. Among them, three genes (At1g49580, At2g01830 and At2g20050) are annotated as harboring both PK and PP activity. Gene ontology analysis revealed that those genes whose products are localized on the plasma membrane, the cell surface, the ER, the cell wall, or are associated with the ubiquitin ligase and calcineurin complexes were highly enriched (P < 0.01; Additional file 11). Purple acid phosphatases, members of the largest class of plant acid phosphatases, are generally assumed to be involved in intra-and/or extra-cellular Pi scavenging and recycling of Pi under Pi-deficient conditions. In Arabidopsis, the PAP family is composed of 29 members sharing conserved domains. Most PAPs are induced by Pi deficiency, some of which in an organ-specific manner [31, 32]. Precise digital expression information of the 29 members of the PAP family in Pi-deficient Arabidopsis rootsis presented by RNA-seq (Additional file 12), which completes gene expression information in Arabidopsis roots so far uncovered by microarray analysis and classic molecular techniques. It is unclear, however, whether the members of the PAP family harbor PP catalytic activity.
Construction of Pi-responsive PK and PP co-expression networks
Co-expression networks of the differentially expressed PK and PP genes were constructed using the MACCU software . Co-expressed genes were selected with a Pearson correlation coefficient cutoff of 0.7. This cutoff also has been used in earlier studies [9, 33]. Co-expression networks constructed with this cutoff are well suited to guide follow-up experiments (i.e. networksare neither too big nor too small). It should be mentioned that the co-expression network constructed here is restricted to roots and the 300 public microarrays mined for generating co-expression relationship were root-related experiments . Because protein regulation by phosphorylation is reversible and requires both PKs and PPs, the network was constructed from both PK and PP genes. The 205 differentially expressed PK and PP genes were loaded as guide genes to calculate the correlations. Correlations between guide genes were visualized by Cytoscape (http://www.cytoscape.org). In the co-expression network, a node represents a gene and an edge represents the correlation between two genes. The network of PK and PP genes responsive to Pi deficiency consists of 65 nodes and 96 edges (Additional file 13). The 65 nodes contain 60 genes that encode PKs and five genes that encode PPs (Additional file 14). The network can further be divided into two larger and seven smaller modules. Genes with similar expression pattern under diverse conditions can have correlative functions and may form a functional module . Only 34% of the input genes were constituents of the co-expression network, suggesting that the majority of PKs and PPs responsive to Pi deficiency are functionally diverse and involved in a variety of biological processes and metabolic pathways. It is noteworthy that, compared to14% of the input PP genes associated with co-expression network, the group of PKs is represented by 34% of the genes, even though the number of differentially expressed PP gene is slightly higher than that of PK genes (35 out of 205 PP genes and 173 out of 1,118 PK genes). Some modules contain only a few or none PP genes. For instance, the largest module contains only one PP gene, PAP11. These observations suggest that the regulation of biological processes may require a cascaded and/or coordinated protein phosphorylation by different PKs to adapt to environmental stresses, while the removal of phosphate from a phosphorylated protein by a PP is less specific. Because of the importance of PK in signaling and metabolism, it is reasonable to speculate that protein phosphorylation is one of the proprietary processes for plant cell to use the limited Pi under Pi deficiency.
Genes involved in unidimensional cell growth form the major module
To gain insight into the function of the genes associated with the network, GO enrichment analyses were performed. The biological processes ‘glycolysis’, ‘negative regulation of anion channels activity by blue light’, ‘unidimensional cell growth’, and ‘glucosinolate biosynthetic process’ were most strongly enriched (Additional file 15). This analysis supports the assumption that protein phosphorylation is an important regulatory level for diverse processes associated with the recalibration of the cellular Pi homeostasis. It has been documented that Pi deficiency alters cellular metabolism, mainly by enhanced carbon flux via glycolysis for increased synthesis of organic acids . Our data support these findings and further suggest that under Pi deficiency protein phosphorylation is an important regulator of glycolytic flux.
Protein kinases involved in pollen tube development and growth are co-expressed in Pi-deficient roots
Genes comprising PKPP1
Mean ± SD (−Pi/+Pi; P < 0.01)
Protein kinase superfamily protein
1.93 ± 0.92
AtPERK12, IGI1, Protein kinase superfamily protein
1.60 ± 0.38
ATCIPK18, ATWL1, CIPK18, SnRK3.20, WL1, CBL-interacting protein kinase 18
1.50 ± 0.45
S-locus lectin protein kinase family protein
1.99 ± 0.51
Leucine-rich repeat protein kinase family protein
3.03 ± 0.79
ANX1, Malectin/receptor-like protein kinase family protein
1.36 ± 0.08
Protein kinase superfamily protein with octicosapeptide/Phox/Bem1p domain
0.65 ± 0.03
AGC (cAMP-dependent, cGMP-dependent and protein kinase C) kinase family protein
3.65 ± 0.62
Protein kinase superfamily protein
2.19 ± 0.47
ATPAP11, PAP11, purple acid phosphatase 11
∞ (De novo synthesis)
Leucine-rich receptor-like protein kinase family protein
0.87 ± 0.04
Protein kinase superfamily protein
1.39 ± 0.02
U-box domain-containing protein kinase family protein
1.33 ± 0.13
U-box domain-containing protein kinase family protein
3.52 ± 0.39
Leucine-rich repeat protein kinase family protein
1.75 ± 0.12
U-box domain-containing protein kinase family protein
6.13 ± 1.12
CPK25, calcium-dependent protein kinase 25
1.39 ± 0.08
RHS10, root hair specific 10
1.43 ± 0.18
AGC1.7, AGC kinase 1.7
1.98 ± 0.44
IRE, AGC (cAMP-dependent, cGMP-dependent and protein kinase C) kinase family protein
1.38 ± 0.25
Protein kinase superfamily protein
1.38 ± 0.13
RSH3, root hair specific 3
1.42 ± 0.17
Protein kinases from the AGC and PERK subfamilies are critical for root hair development and growth under Pi deficiency
Members of sub-module PKPP1B (Figure 3B and Table 1) may play key roles in root hair development and elongation under both normal and Pi-deficient conditions. Three genes in this submodule, including the two root hair-specific genes RHS3 and RHS10 and IRE1 (INCOMPLETE ROOT HAIR ELONGATION 1) , have been related to root hair development and elongation. RHS3 and IRE1 belong to the AGC family while RSH10 belongs to the PERK family, indicating that kinases in these families are particularly important for root hair development and growth. All three genes were induced by Pi deficiency at the transcriptional level. Reverse genetic studies would be of help to decipher their physiological functions under Pi deprivation. Several other genes in this sub-module, including At4g25160, At5g61550, At4g31250, and At2g41970, have been reported to be part of a gene regulatory network comprising 208 root epidermal ‘core’ genes in Arabidopsis . All genes in this sub-module were induced by Pi-deficiency, implicating their involvement in root hair development and growth under Pi starvation. AGC1.7, belonging to AGCVIII subfamily, and its homologue AGC1.5 have been reported to be critical for the polarized growth of pollen tubes , but have not yet been associated with root hair development. Our analysis revealed that AGC1.7 was co-expressed with RSH3 and CPK25 (calcium-independent CDPK ) under Pi deficiency, suggesting that AGC1.7 supports a possible function in root hair development under Pi limiting conditions. The only gene that was repressed by Pi deficiency in the sub-module PKPP1B is CLV1. By binding to a small protein ligand CLV3, CLV1restricts the proliferation and/or promotes the differentiation of stem cells in the shoot apical meristem . Interestingly, RNA-seq analysis revealed that CLV1 was also highly expressed in Arabidopsis roots and decreased in response to Pi deficiency (Table 1), suggesting that CLV1 might be negatively regulating root hair development in response to Pi starvation.
Root epidermial core gene associated with PKs and PPs under Pi deficiency
In summary, we here provide precise digital information on the transcription of protein kinase and phosphatase genes in Arabidopsis roots at a genome-wide level. A root-specific co-expression network of Pi-responsive genes encoding protein kinases and phosphatases has been generated and putative novel players in Pi deficiency-induced root hair formation have been uncovered. Combined with the previously published root hair core genes , a comprehensive inventory for the regulation of root hair development and metabolism was obtained. The approach applied here will be useful to direct further studies by reverse genetic methods and eventually decipher the mechanisms by which root epidermal cells are re-programmed to adapt to Pi deficiency.
Data collection and processing
Transcriptome data of roots from plants grown in the presence or absence of Pi from 13-d-old Arabidopsis seedlings by RNA-seq were downloaded from a public database (NCBI: SRA050356.1) and analyzed as described in . Microarray data of 2,671 ATH1 arrays from the NASCarray database (http://affymetrix.arabidopsis.info/) were downloaded and normalized using the RMA function of the Affy package of the Bioconductor software. Among the 2,671 arrays, 300 root-related arrays were manually identified as described in . PK and PP genes were retrieved on the basis of TAIR 10 release of Arabidopsis genome.
Generation of co-expression networks and modules of Pi-responsive PK and PP genes using the MACCU toolbox
To generate root-specific networks of Pi-responsive PK and PP genes, differentially expressed PK and PP genes in the Arabidopsis roots were obtained using a Student t-test at a P value <0.05. Gene networks were constructed based on 300 publicly available root-related arrays using the MACCU toolbox as described in , with a Pearson correlation threshold of 0.7. The generated co-expression networks were visualized by Cytoscape (http://www.cytoscape.org). If one cluster of genes did not have any connection (without any edges) to any other cluster in the co-expression network, we referred to such a cluster as a module.
Construction of root hair-specific networks of Pi-responsive PK and PP genes
To obtain a root hair cell-specific network of Pi-responsive PK and PP genes, the 208 root epidermal ‘core’ genes were first mined from the datasheet described in . Next, Pi-responsive PK and PP genes involved in the module of pollen tube/root hair development and growth were extracted (bait genes), combined with the genes from the root epidermal ‘core’ (preys), and used for generating co-expression network using the MACCU toolbox with a Pearson correlation threshold of 0.7. The resulting networks shows only those nodes (genes) and edges (relationships between genes) that were linked by at least one edge must with bait. Edges linked to two preys were excluded.
The study was partially supported by the starting career grant from the Institute of Soil Science, Chinese Academy Sciences (Y225070000). We thank Drs. Wen-Dar Lin and Jorge Rodríguez Celma for their help in using the MACCU software and two anonymous reviewers for their valuable comments.
- Chiou TJ, Lin SI: Signaling network in sensing phosphate availability in plants. Ann Rev Plant Biol. 2011, 62: 185-206. 10.1146/annurev-arplant-042110-103849.View ArticleGoogle Scholar
- Yao Y, Sun H, Xu F, Zhang X, Liu S: Comparative proteome analysis of metabolic changes by low phosphorus stress in two Brassica napus genotypes. Planta. 2011, 233 (3): 523-537. 10.1007/s00425-010-1311-x.View ArticlePubMedGoogle Scholar
- Li K, Xu C, Zhang K, Yang A, Zhang J: Proteomic analysis of roots growth and metabolic changes under phosphorus deficit in maize (Zea mays L.) plants. Proteomics. 2007, 7 (9): 1501-1512. 10.1002/pmic.200600960.View ArticlePubMedGoogle Scholar
- Torabi S, Wissuwa M, Heidari M, Naghavi MR, Gilany K, Hajirezaei MR, Omidi M, Yazdi-Samadi B, Ismail AM, Salekdeh GH: A comparative proteome approach to decipher the mechanism of rice adaptation to phosphorous deficiency. Proteomics. 2009, 9 (1): 159-170. 10.1002/pmic.200800350.View ArticlePubMedGoogle Scholar
- Muller R, Morant M, Jarmer H, Nilsson L, Nielsen TH: Genome-wide analysis of the Arabidopsis leaf transcriptome reveals interaction of phosphate and sugar metabolism. Plant Physiol. 2007, 143 (1): 156-171.PubMed CentralView ArticlePubMedGoogle Scholar
- Morcuende R, Bari R, Gibon Y, Zheng W, Pant BD, Blasing O, Usadel B, Czechowski T, Udvardi MK, Stitt M: Genome-wide reprogramming of metabolism and regulatory networks of Arabidopsis in response to phosphorus. Plant Cell Environ. 2007, 30 (1): 85-112. 10.1111/j.1365-3040.2006.01608.x.View ArticlePubMedGoogle Scholar
- Misson J, Raghothama KG, Jain A, Jouhet J, Block MA, Bligny R, Ortet P, Creff A, Somerville S, Rolland N: A genome-wide transcriptional analysis using Arabidopsis thaliana Affymetrix gene chips determined plant responses to phosphate deprivation. Proc Nat Acad Sci USA. 2005, 102 (33): 11934-11939. 10.1073/pnas.0505266102.PubMed CentralView ArticlePubMedGoogle Scholar
- Bari R, Datt Pant B, Stitt M, Scheible WR: PHO2, microRNA399, and PHR1 define a phosphate-signaling pathway in plants. Plant Physiol. 2006, 141 (3): 988-999. 10.1104/pp.106.079707.PubMed CentralView ArticlePubMedGoogle Scholar
- Lin WD, Liao YY, Yang TJ, Pan CY, Buckhout TJ, Schmidt W: Coexpression-based clustering of Arabidopsis root genes predicts functional modules in early phosphate deficiency signaling. Plant Physiol. 2011, 155 (3): 1383-1402. 10.1104/pp.110.166520.PubMed CentralView ArticlePubMedGoogle Scholar
- Hammond JP, Bennett MJ, Bowen HC, Broadley MR, Eastwood DC, May ST, Rahn C, Swarup R, Woolaway KE, White PJ: Changes in gene expression in Arabidopsis shoots during phosphate starvation and the potential for developing smart plants. Plant Physiol. 2003, 132 (2): 578-596. 10.1104/pp.103.020941.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu P, Ma L, Hou X, Wang M, Wu Y, Liu F, Deng XW: Phosphate starvation triggers distinct alterations of genome expression in Arabidopsis roots and leaves. Plant Physiol. 2003, 132 (3): 1260-1271. 10.1104/pp.103.021022.PubMed CentralView ArticlePubMedGoogle Scholar
- Plaxton WC, Tran HT: Metabolic adaptations of phosphate-starved plants. Plant Physiol. 2011, 156 (3): 1006-1015. 10.1104/pp.111.175281.PubMed CentralView ArticlePubMedGoogle Scholar
- Bayle V, Arrighi JF, Creff A, Nespoulous C, Vialaret J, Rossignol M, Gonzalez E, Paz-Ares J, Nussaume L: Arabidopsis thaliana high-affinity phosphate transporters exhibit multiple levels of posttranslational regulation. Plant Cell. 2011, 23 (4): 1523-1535. 10.1105/tpc.110.081067.PubMed CentralView ArticlePubMedGoogle Scholar
- Svistoonoff S, Creff A, Reymond M, Sigoillot-Claude C, Ricaud L, Blanchet A, Nussaume L, Desnos T: Root tip contact with low-phosphate media reprograms plant root architecture. Nature Gen. 2007, 39 (6): 792-796. 10.1038/ng2041.View ArticleGoogle Scholar
- Rubio V, Linhares F, Solano R, Martin AC, Iglesias J, Leyva A, Paz-Ares J: A conserved MYB transcription factor involved in phosphate starvation signaling both in vascular plants and in unicellular algae. Gen Dev. 2001, 15 (16): 2122-2133. 10.1101/gad.204401.View ArticleGoogle Scholar
- Bustos R, Castrillo G, Linhares F, Puga MI, Rubio V, Perez-Perez J, Solano R, Leyva A, Paz-Ares J: A central regulatory system largely controls transcriptional activation and repression responses to phosphate starvation in Arabidopsis. PLoS Genet. 2010, 6 (9): e1001102-10.1371/journal.pgen.1001102.PubMed CentralView ArticlePubMedGoogle Scholar
- Nilsson L, Muller R, Nielsen TH: Increased expression of the MYB-related transcription factor, PHR1, leads to enhanced phosphate uptake in Arabidopsis thaliana. Plant Cell Environ. 2007, 30 (12): 1499-1512. 10.1111/j.1365-3040.2007.01734.x.View ArticlePubMedGoogle Scholar
- Miura K, Rus A, Sharkhuu A, Yokoi S, Karthikeyan AS, Raghothama KG, Baek D, Koo YD, Jin JB, Bressan RA: The Arabidopsis SUMO E3 ligase SIZ1 controls phosphate deficiency responses. Proc Nat Acad Sci USA. 2005, 102 (21): 7760-7765. 10.1073/pnas.0500778102.PubMed CentralView ArticlePubMedGoogle Scholar
- Aung K, Lin SI, Wu CC, Huang YT, Su CL, Chiou TJ: pho2, a phosphate overaccumulator, is caused by a nonsense mutation in a microRNA399 target gene. Plant Physiol. 2006, 141 (3): 1000-1011. 10.1104/pp.106.078063.PubMed CentralView ArticlePubMedGoogle Scholar
- Pant BD, Buhtz A, Kehr J, Scheible WR: MicroRNA399 is a long-distance signal for the regulation of plant phosphate homeostasis. Plant J Cell Mol Biol. 2008, 53 (5): 731-738. 10.1111/j.1365-313X.2007.03363.x.View ArticleGoogle Scholar
- Li WF, Perry PJ, Prafulla NN, Schmidt W: Ubiquitin-specific protease 14 (UBP14) is involved in root responses to phosphate deficiency in Arabidopsis. Mol Plant. 2010, 3 (1): 212-223. 10.1093/mp/ssp086.View ArticlePubMedGoogle Scholar
- Lan P, Li W, Schmidt W: Complementary proteome and transcriptome profiling in phosphate-deficient Arabidopsis roots reveals multiple levels of gene regulation. MCP. 2012, 11 (11): 1156-1166.PubMed CentralPubMedGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Nat Acad Sci USA. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.PubMed CentralView ArticlePubMedGoogle Scholar
- Usadel B, Obayashi T, Mutwil M, Giorgi FM, Bassel GW, Tanimoto M, Chow A, Steinhauser D, Persson S, Provart NJ: Co-expression tools for plant biology: opportunities for hypothesis generation and caveats. Plant Cell Environ. 2009, 32 (12): 1633-1651. 10.1111/j.1365-3040.2009.02040.x.View ArticlePubMedGoogle Scholar
- Aoki K, Ogata Y, Shibata D: Approaches for extracting practical information from gene co-expression networks in plant biology. Plant Cell Physiol. 2007, 48 (3): 381-390. 10.1093/pcp/pcm013.View ArticlePubMedGoogle Scholar
- Ozsolak F, Milos PM: RNA sequencing: advances, challenges and opportunities. Nat Rev Gen. 2011, 12 (2): 87-98. 10.1038/nrg2934.View ArticleGoogle Scholar
- Tanaka H, Osakabe Y, Katsura S, Mizuno S, Maruyama K, Kusakabe K, Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K: Abiotic stress-inducible receptor-like kinases negatively control ABA signaling in Arabidopsis. Plant J Mol Biol. 2012, 70 (4): 599-613. 10.1111/j.1365-313X.2012.04901.x.View ArticleGoogle Scholar
- Wrzaczek M, Brosche M, Salojarvi J, Kangasjarvi S, Idanheimo N, Mersmann S, Robatzek S, Karpinski S, Karpinska B, Kangasjarvi J: Transcriptional regulation of the CRK/DUF26 group of receptor-like protein kinases by ozone and plant hormones in Arabidopsis. BMC Plant Biol. 2010, 10: 95-10.1186/1471-2229-10-95.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen K, Du L, Chen Z: Sensitization of defense responses and activation of programmed cell death by a pathogen-induced receptor-like protein kinase in Arabidopsis. Plant Mol Biol. 2003, 53 (1–2): 61-74.View ArticlePubMedGoogle Scholar
- Sassi M, Lu Y, Zhang Y, Wang J, Dhonukshe P, Blilou I, Dai M, Li J, Gong X, Jaillais Y: COP1 mediates the coordination of root and shoot growth by light through modulation of PIN1- and PIN2-dependent auxin transport in Arabidopsis. Development. 2012, 139 (18): 3402-3412. 10.1242/dev.078212.View ArticlePubMedGoogle Scholar
- Li D, Zhu H, Liu K, Liu X, Leggewie G, Udvardi M, Wang D: Purple acid phosphatases of Arabidopsis thaliana. Comparative analysis and differential regulation by phosphate deprivation. J Biol Chem. 2002, 277 (31): 27772-27781.PubMedGoogle Scholar
- Zhu H, Qian W, Lu X, Li D, Liu X, Liu K, Wang D: Expression patterns of purple acid phosphatase genes in Arabidopsis organs and functional analysis of AtPAP23 predominantly transcribed in flower. Plant Mol Biol. 2005, 59 (4): 581-594. 10.1007/s11103-005-0183-0.View ArticlePubMedGoogle Scholar
- Wang S, Yin Y, Ma Q, Tang X, Hao D, Xu Y: Genome-scale identification of cell-wall related genes in Arabidopsis based on co-expression network analysis. BMC Plant Biol. 2012, 12: 138-10.1186/1471-2229-12-138.PubMed CentralView ArticlePubMedGoogle Scholar
- Hammond JP, White PJ: Sucrose transport in the phloem: integrating root responses to phosphorus starvation. J Exp Bot. 2008, 59 (1): 93-109.View ArticlePubMedGoogle Scholar
- Oyama T, Shimura Y, Okada K: The IRE gene encodes a protein kinase homologue and modulates root hair growth in Arabidopsis. Plant J Mol Biol. 2002, 30 (3): 289-299. 10.1046/j.1365-313X.2002.01290.x.View ArticleGoogle Scholar
- Won SK, Lee YJ, Lee HY, Heo YK, Cho M, Cho HT: Cis-element- and transcriptome-based screening of root hair-specific genes and their functional characterization in Arabidopsis. Plant Physiol. 2009, 150 (3): 1459-1473. 10.1104/pp.109.140905.PubMed CentralView ArticlePubMedGoogle Scholar
- Nozawa A, Sawada Y, Akiyama T, Koizumi N, Sano H: Variable interactions between sucrose non-fermented 1-related protein kinases and regulatory proteins in higher plants. Biosci Biotechnol Biochem. 2003, 67 (12): 2533-2540. 10.1271/bbb.67.2533.View ArticlePubMedGoogle Scholar
- Lehti-Shiu MD, Shiu SH: Diversity, classification and function of the plant protein kinase superfamily. Phil Transac Royal Soc London Series B Biol Sci. 2012, 367 (1602): 2619-2639. 10.1098/rstb.2012.0003.View ArticleGoogle Scholar
- Lehti-Shiu MD, Zou C, Hanada K, Shiu SH: Evolutionary history and stress regulation of plant receptor-like kinase/pelle genes. Plant Physiol. 2009, 150 (1): 12-26. 10.1104/pp.108.134353.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 (5): 1220-1234. 10.1105/tpc.020834.PubMed CentralView ArticlePubMedGoogle Scholar
- Shiu SH, Bleecker AB: Expansion of the receptor-like kinase/Pelle gene family and receptor-like proteins in Arabidopsis. Plant Physiol. 2003, 132 (2): 530-543. 10.1104/pp.103.021964.View ArticlePubMedGoogle Scholar
- Shiu SH, Bleecker AB: Receptor-like kinases from Arabidopsis form a monophyletic gene family related to animal receptor kinases. Proc Nat Acad Sci USA. 2001, 98 (19): 10763-10768. 10.1073/pnas.181141598.PubMed CentralView ArticlePubMedGoogle Scholar
- Shiu SH, Bleecker AB: Plant receptor-like kinase gene family: diversity, function, and signaling. Sci STKE. 2001, 2001 (113): re22-PubMedGoogle Scholar
- Wang H, Chevalier D, Larue C, Ki Cho S, Walker JC: The Protein Phosphatases and Protein Kinases of Arabidopsis thaliana. Arabidopsis Book. 2007, 5: e0106-PubMed CentralPubMedGoogle Scholar
- Bruex A, Kainkaryam RM, Wieckowski Y, Kang YH, Bernhardt C, Xia Y, Zheng X, Wang JY, Lee MM, Benfey P: A gene regulatory network for root epidermis cell differentiation in Arabidopsis. PLoS Genet. 2012, 8 (1): e1002446-10.1371/journal.pgen.1002446.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim DW, Lee SH, Choi SB, Won SK, Heo YK, Cho M, Park YI, Cho HT: Functional conservation of a root hair cell-specific cis-element in angiosperms with different root hair distribution patterns. Plant Cell. 2006, 18 (11): 2958-2970. 10.1105/tpc.106.045229.PubMed CentralView ArticlePubMedGoogle Scholar
- Hwang I, Kim SY, Kim CS, Park Y, Tripathi GR, Kim SK, Cheong H: Over-expression of the IGI1 leading to altered shoot-branching development related to MAX pathway in Arabidopsis. Plant Mol Biol. 2010, 73 (6): 629-641. 10.1007/s11103-010-9645-0.PubMed CentralView ArticlePubMedGoogle Scholar
- Jain A, Poling MD, Karthikeyan AS, Blakeslee JJ, Peer WA, Titapiwatanakun B, Murphy AS, Raghothama KG: Differential effects of sucrose and auxin on localized phosphate deficiency-induced modulation of different traits of root system architecture in Arabidopsis. Plant Physiol. 2007, 144 (1): 232-247. 10.1104/pp.106.092130.PubMed CentralView ArticlePubMedGoogle Scholar
- Jones AR, Kramer EM, Knox K, Swarup R, Bennett MJ, Lazarus CM, Leyser HM, Grierson CS: Auxin transport through non-hair cells sustains root-hair development. Nat Cell Biol. 2009, 11 (1): 78-84. 10.1038/ncb1815.PubMed CentralView ArticlePubMedGoogle Scholar
- Bogre L, Okresz L, Henriques R, Anthony RG: Growth signalling pathways in Arabidopsis and the AGC protein kinases. Trends Plant Sci. 2003, 8 (9): 424-431. 10.1016/S1360-1385(03)00188-2.View ArticlePubMedGoogle Scholar
- Jacinto E, Lorberg A: TOR regulation of AGC kinases in yeast and mammals. Biochem J. 2008, 410 (1): 19-37. 10.1042/BJ20071518.View ArticlePubMedGoogle Scholar
- Pearce LR, Komander D, Alessi DR: The nuts and bolts of AGC protein kinases. Nat Rev Mol Cell Biol. 2010, 11 (1): 9-22.View ArticlePubMedGoogle Scholar
- Garcia AV, Al-Yousif M, Hirt H: Role of AGC kinases in plant growth and stress responses. CMLS. 2012, 69 (19): 3259-3267. 10.1007/s00018-012-1093-3.View ArticlePubMedGoogle Scholar
- Boisson-Dernier A, Roy S, Kritsas K, Grobei MA, Jaciubek M, Schroeder JI, Grossniklaus U: Disruption of the pollen-expressed FERONIA homologs ANXUR1 and ANXUR2 triggers pollen tube discharge. Dev (Cambridge, England). 2009, 136 (19): 3279-3288. 10.1242/dev.040071.View ArticleGoogle Scholar
- Miyazaki S, Murata T, Sakurai-Ozato N, Kubo M, Demura T, Fukuda H, Hasebe M: ANXUR1 and 2, sister genes to FERONIA/SIRENE, are male factors for coordinated fertilization. CB. 2009, 19 (15): 1327-1331.PubMedGoogle Scholar
- Zhang Y, He J, McCormick S: Two Arabidopsis AGC kinases are critical for the polarized growth of pollen tubes. Plant J Cell Mol Biol. 2009, 58 (3): 474-484. 10.1111/j.1365-313X.2009.03792.x.View ArticleGoogle Scholar
- Boudsocq M, Droillard MJ, Regad L, Lauriere C: Characterization of Arabidopsis calcium-dependent protein kinases: activated or not by calcium?. The. Biochem J. 2012, 447 (2): 291-299. 10.1042/BJ20112072.View ArticlePubMedGoogle Scholar
- Rojo E, Sharma VK, Kovaleva V, Raikhel NV, Fletcher JC: CLV3 is localized to the extracellular space, where it activates the Arabidopsis CLAVATA stem cell signaling pathway. Plant Cell. 2002, 14 (5): 969-977. 10.1105/tpc.002196.PubMed CentralView ArticlePubMedGoogle Scholar
- Ihmels J, Levy R, Barkai N: Principles of transcriptional control in the metabolic network of Saccharomyces cerevisiae. Nat Biotechnol. 2004, 22 (1): 86-92. 10.1038/nbt918.View ArticlePubMedGoogle Scholar
- Kharchenko P, Church GM, Vitkup D: Expression dynamics of a cellular metabolic network. Mol Sys Biol. 2005, 1: 2005.0016-Google Scholar
- Hirai MY, Sugiyama K, Sawada Y, Tohge T, Obayashi T, Suzuki A, Araki R, Sakurai N, Suzuki H, Aoki K: Omics-based identification of Arabidopsis Myb transcription factors regulating aliphatic glucosinolate biosynthesis. Proc Nat Acad Sci USA. 2007, 104 (15): 6478-6483. 10.1073/pnas.0611629104.PubMed CentralView ArticlePubMedGoogle Scholar
- Fu FF, Xue HW: Coexpression analysis identifies Rice Starch Regulator1, a rice AP2/EREBP family transcription factor, as a novel rice starch biosynthesis regulator. Plant Physiol. 2010, 154 (2): 927-938. 10.1104/pp.110.159517.PubMed CentralView ArticlePubMedGoogle Scholar
- Han X, Yin L, Xue H: Co-expression analysis identifies CRC and AP1 the regulator of Arabidopsis fatty acid biosynthesis. J Integr Plant Biol. 2012, 54 (7): 486-499. 10.1111/j.1744-7909.2012.01132.x.View ArticlePubMedGoogle Scholar
- Kulich I, Cole R, Drdova E, Cvrckova F, Soukup A, Fowler J, Zarsky V: Arabidopsis exocyst subunits SEC8 and EXO70A1 and exocyst interactor ROH1 are involved in the localized deposition of seed coat pectin. New Phytol. 2010, 188 (2): 615-625. 10.1111/j.1469-8137.2010.03372.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.