Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Uncovering Small RNA-Mediated Responses to Phosphate Deficiency in Arabidopsis by Deep Sequencing

Uncovering Small RNA-Mediated Responses to Phosphate Deficiency in Arabidopsis by Deep Sequencing Abstract Recent studies have demonstrated the important role of plant microRNAs (miRNAs) under nutrient deficiencies. In this study, deep sequencing of Arabidopsis (Arabidopsis thaliana) small RNAs was conducted to reveal miRNAs and other small RNAs that were differentially expressed in response to phosphate (Pi) deficiency. About 3.5 million sequence reads corresponding to 0.6 to 1.2 million unique sequence tags from each Pi-sufficient or Pi-deficient root or shoot sample were mapped to the Arabidopsis genome. We showed that upon Pi deprivation, the expression of miR156, miR399, miR778, miR827, and miR2111 was induced, whereas the expression of miR169, miR395, and miR398 was repressed. We found cross talk coordinated by these miRNAs under different nutrient deficiencies. In addition to miRNAs, we identified one Pi starvation-induced DICER-LIKE1-dependent small RNA derived from the long terminal repeat of a retrotransposon and a group of 19-nucleotide small RNAs corresponding to the 5′ end of tRNA and expressed at a high level in Pi-starved roots. Importantly, we observed an increased abundance of TAS4-derived trans-acting small interfering RNAs (ta-siRNAs) in Pi-deficient shoots and uncovered an autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(−) that regulates the biosynthesis of anthocyanin. This finding sheds light on the regulatory network between miRNA/ta-siRNA and its target gene. Of note, a substantial amount of miR399* accumulated under Pi deficiency. Like miR399, miR399* can move across the graft junction, implying a potential biological role for miR399*. This study represents a comprehensive expression profiling of Pi-responsive small RNAs and advances our understanding of the regulation of Pi homeostasis mediated by small RNAs. Plants acquire mineral nutrients that are channeled into diverse biochemical and metabolic pathways essential for the completion of their life cycle. Phosphorus (P), one of the essential macronutrients of plants, is often not readily accessible to plants because of the low availability of phosphate (Pi) in the soil, the major acquired form of P (Marschner, 1995). To overcome this limitation, plants have evolved a number of responses to conserve and mobilize internal Pi and to enhance the acquisition of external Pi (Raghothama, 1999; Poirier and Bucher, 2002; Lin et al., 2009). Recently, the molecular basis underlying these adaptive responses has been an intense subject of research. It is now recognized that small regulatory RNAs contribute to the complexity of gene regulation. MicroRNAs (miRNAs) and small interfering RNAs (siRNAs) of 21 to 24 nucleotides long are two major classes of small regulatory RNAs in plants, functioning as negative regulators of gene expression (Brodersen and Voinnet, 2006; Jones-Rhoades et al., 2006; Mallory and Vaucheret, 2006). miRNAs are processed from single-stranded RNA precursors capable of forming imperfectly complementary hairpin structures by the RNase III enzyme DICER-LIKE1 (DCL1) or DCL4. In most cases, plant miRNAs down-regulate gene expression by direct cleavage of their target transcripts, although translational repression has been proposed to be common as well (Brodersen et al., 2008). In contrast, siRNAs are generated from perfectly complementary long double-stranded RNA precursors, requiring the activity of RNA-dependent RNA polymerase. Endogenous siRNAs can be further classified into trans-acting siRNAs (ta-siRNAs), natural antisense transcript-derived siRNAs, and heterochromatin siRNAs according to their genomic origin and the interaction of distinct components (Vaucheret, 2006; Vazquez, 2006). ta-siRNAs are phased 21-nucleotide RNA molecules whose synthesis is triggered by the cleavage of miRNA on the TAS transcript (Peragine et al., 2004; Vazquez et al., 2004; Allen et al., 2005). They act in-trans to down-regulate the expression of unrelated loci from which they are produced. The biological functions of miRNA cover a broad range of physiological processes. In addition to their involvement in normal plant growth and development, the expression of many miRNAs is regulated by various biotic and abiotic stresses, which may contribute to the development of adaptive responses to overcome unfavorable conditions (Sunkar and Zhu, 2004; Navarro et al., 2006; Chiou, 2007; Li et al., 2008a; Lu et al., 2008). Recent findings have also demonstrated the important role of several miRNAs under nutrient deficiencies. miR395, miR398, and miR399 are up-regulated under sulfur (S), copper (Cu), and Pi deficiency, respectively. miR395 is involved in S assimilation and allocation by adjusting the activity of ATP sulfurylase (APS) and restricting the spatial expression of a sulfate transporter (AtSULTR2;1) responsible for root-to-shoot transport of sulfate (Jones-Rhoades and Bartel, 2004; Kawashima et al., 2009). Other than a role in alleviating oxidative stresses (Sunkar et al., 2006), miR398 modulates the expression of Cu/zinc (Zn) superoxide dismutase according to the availability of Cu (Yamasaki et al., 2007). Moreover, cell-specific expression of miR169 is crucial in the development of nitrogen (N)-fixing nodules in Medicago truncatula (Combier et al., 2006). miR399 regulates Pi homeostasis by controlling the expression of PHO2 encoding a ubiquitin-conjugating E2 enzyme, UBC24 (Fujii et al., 2005; Aung et al., 2006; Bari et al., 2006; Chiou et al., 2006). It acts as a positive regulator to activate the Pi uptake and translocation when Pi supply is insufficient, while PHO2 functions as a negative regulator to suppress these activities in order to prevent Pi excess when external Pi is ample. Notably, recent results from reciprocal grafting suggest that miR399 may serve as a long-distance signal from shoots to suppress PHO2 expression in roots as a way of communication between the Pi status of shoots and the Pi transport activity of roots (Lin et al., 2008; Pant et al., 2008). The long-distance movement of miR399 could be an early response to Pi deficiency. The shoot-derived miR399 is responsible for the degradation of PHO2 in the roots, where miR399 is not readily expressed (Lin et al., 2008). New technologies of high-throughput sequencing have become powerful tools to disclose the large inventory of small RNA species in plants. These deep sequencing strategies may discover new miRNAs or novel small RNA classes and provide quantitative profiling of small RNA expression (Lu et al., 2005; Rajagopalan et al., 2006; Nobuta et al., 2007; Heisel et al., 2008; Zhu et al., 2008). In this study, Solexa sequencing was employed to profile the changes of small RNA expression in response to Pi deficiency. Very recently, an independent study was reported in which P-responsive miRNAs were identified using quantitative reverse transcription (RT)-PCR of primary transcript of miRNAs followed by small RNA sequencing (Pant et al., 2009). In our study, we not only profiled the expression of small RNAs but also looked into their regulatory mechanisms and potential functions. Importantly, we found two unique classes of small RNAs up-regulated by Pi deficiency and uncovered an autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(−) that regulates the biosynthesis of anthocyanin. Changes in the abundance of specific groups of Pi-responsive small RNAs may mediate the activity or expression of target genes that are important in the development or regulation of adaptive responses, enabling plants to survive in a low-Pi environment. RESULTS AND DISCUSSION Overview of Small RNA Profiles in Response to Pi Deficiency Four small RNA libraries from Pi-sufficient roots (R+Pi), Pi-deficient roots (R−Pi), Pi-sufficient shoots (S+Pi), and Pi-deficient shoots (S−Pi) were constructed by Solexa high-throughput sequencing technology. After adaptor trimming, the numbers of sequence reads and unique sequence signatures from the raw data were calculated and then mapped to the Arabidopsis (Arabidopsis thaliana) genome (Supplemental Table S1). In summary, the sequencing yielded 4.8 to 5.1 million raw reads from each library, and approximately 90% of the reads remained after trimming the adaptor sequences. Approximately 3.5 million sequence reads (approximately 77%–83% of trimmed sequence reads) corresponding to 0.6 to 1.2 million unique sequence signatures could be mapped perfectly onto the genome (The Arabidopsis Information Resource 8). Only those perfectly mapped sequences were analyzed further. The sequence reads that mapped to single and multiple genomic locations were defined, respectively, as single-hit reads and multiple-hit reads. One-third of mapped sequences were single-hit reads in the root and one-half were single-hit reads in the shoot. The length of the mapped sequences ranged from 16 to 27 nucleotides. We first examined the correlation between the length of small RNAs and the proportion of total sequence reads or unique sequence signatures (Fig. 1 Figure 1. Open in new tabDownload slide Global profiling of small RNAs. Correlations between the small RNA length and the proportion of total sequence reads (A) or unique sequence signatures (B) are shown. nt, Nucleotides. Figure 1. Open in new tabDownload slide Global profiling of small RNAs. Correlations between the small RNA length and the proportion of total sequence reads (A) or unique sequence signatures (B) are shown. nt, Nucleotides. ). Similar distributions were observed in +Pi and −Pi libraries. Consistent with the earlier reports, there were two major peaks at 21 and 24 nucleotides in the total sequence reads of shoot libraries (Fig. 1A). Unexpectedly, an additional peak at 19 nucleotides was observed in the root libraries. This peak has not been noticed in the previous deep sequencing databases, probably because the pure root sample has not been examined. Further analysis revealed that the majority of these 19-nucleotide sequences originated from tRNAs (see below). When the unique sequence signatures were examined, the patterns of all four libraries were nearly identical (Fig. 1B). The 24-nucleotide small RNAs were dominant in either their reads or unique sequences, indicating that they are rich in sequence diversity. Overrepresentation of tRNA-Derived Small RNAs in the Pi-Starved Roots Mapped sequences were classified according to the annotation of the genome (Supplemental Table S2). Most sequences were mapped onto the nonannotated intergenic regions (30.12%–42.30%) and to a lesser extent onto MIRNA genes and protein-coding genes. Strikingly, a significantly higher proportion of small RNAs was mapped to tRNA genes in the root libraries (24.29%–34.32%) than in the shoot libraries (3.33%–5.58%). The percentage in the −Pi root library (34.32%) was even higher than that in the +Pi root library (24.29%). When the origin of these small RNAs was further analyzed (Supplemental Table S3), we did not find any correlation between the abundance of small RNAs from specific tRNA species and their codon usage. We were surprised to find that a 19-nucleotide sequence cleaved from the 5′ end of Gly-tRNATCC represented over 80% of tRNA-derived small RNAs in the roots and accounted for up to 18.44% and 27.70% of total sequence reads in the +Pi and −Pi root libraries, respectively, compared with only 1.00% to 1.79% in the shoot libraries (highlighted in Supplemental Table S3). For validation, we employed RNA gel analyses using the probes from three different segments corresponding to the 5′ and 3′ ends and the central anticodon loop of tRNA (Fig. 2B Figure 2. Open in new tabDownload slide Overrepresentation of tRNA-derived small RNAs in the Pi-starved roots. A, Small RNA gel analyses of tRNA-derived small RNAs in +Pi or −Pi roots or shoots. B, The tRNA regions corresponding to the 5′ and 3′ end probes are designated by the black line, and the anticodon probe is designated by the gray line. The black arrowheads in A and B indicate the 5′ cleaved products and corresponding cleavage site in tRNA, respectively. The vertical gray lines in A indicate the products of tRNA halves from the cleavage at the anticodon loop (gray arrowhead in B). nt, Nucleotides. Figure 2. Open in new tabDownload slide Overrepresentation of tRNA-derived small RNAs in the Pi-starved roots. A, Small RNA gel analyses of tRNA-derived small RNAs in +Pi or −Pi roots or shoots. B, The tRNA regions corresponding to the 5′ and 3′ end probes are designated by the black line, and the anticodon probe is designated by the gray line. The black arrowheads in A and B indicate the 5′ cleaved products and corresponding cleavage site in tRNA, respectively. The vertical gray lines in A indicate the products of tRNA halves from the cleavage at the anticodon loop (gray arrowhead in B). nt, Nucleotides. ). The small RNAs derived from Asp-tRNAGTC were also examined using the 5′ end probe. Consistent with the sequencing data, the 19-nucleotide RNAs from Gly-tRNATCC and Asp-tRNAGTC corresponding to the 5′ end of tRNA were highly accumulated in the roots but much less in the shoots (Fig. 2A). The −Pi roots showed the maximum accumulation. When probing with the 5′ or 3′ end probe, additional bands at 30 to 40 nucleotides likely representing tRNA halves cleaved at the anticodon loop were observed preferentially in the roots (Fig. 2A). These results revealed a spatial and temporal expression pattern of small RNAs derived from the specific cleavage on tRNA molecules rather than random degradation and may represent a novel processing of small RNAs. Studies in microorganisms or mammalian cells have suggested that the tRNA-derived small RNAs are accumulated in a developmentally regulated manner and become dominant in specific tissues or under specific stress conditions (Lee and Collins, 2005; Haiser et al., 2008; Jochl et al., 2008; Li et al., 2008b; Thompson et al., 2008; Pant et al., 2009; Thompson and Parker, 2009). These small RNAs could be associated with the quality control of protein synthesis or regulation of gene expression. Moreover, a recent report showed that tRNA half molecules present in the phloem sap of pumpkin (Cucurbita maxima) were able to inhibit translational activity in vitro (Zhang et al., 2009). The authors proposed that these phloem-delivered tRNA halves may be a long-distance signal to coordinate the metabolic status between source and sink tissues. The differential accumulation of tRNA-derived small RNAs between roots and shoots observed here may represent the consequence of such long-distance movement. Expression Profiling of miRNAs in Response to Pi Deprivation Forty-five and 55 differentially expressed MIRNA genes with greater than 1.5-fold relative change in sequence counts were identified in roots and shoots, respectively (Fig. 3 Figure 3. Open in new tabDownload slide Differentially expressed MIRNA genes in response to Pi deficiency in roots (A) and shoots (B). The significantly differentially expressed MIRNA genes with greater than 1.5-fold relative change are shown. The miRNAs with a relative change ratio greater than 3 are highlighted in boldface. Figure 3. Open in new tabDownload slide Differentially expressed MIRNA genes in response to Pi deficiency in roots (A) and shoots (B). The significantly differentially expressed MIRNA genes with greater than 1.5-fold relative change are shown. The miRNAs with a relative change ratio greater than 3 are highlighted in boldface. ; Supplemental Table S4). MIR399s, MIR778, MIR827, and MIR2111 were clearly up-regulated, whereas MIR169s, MIR395s, and MIR398s were down-regulated, in −Pi roots and/or shoots. The up-regulation of various species from the MIR399 family by −Pi was consistent with our previous observations (Lin et al., 2008), suggesting the reliability of our sequencing data. Among all MIRNA genes, the MIR399 family showed the highest degree of induction in both tissues. The expression of miRNAs with relative changes greater than a ratio of 3 (highlighted in boldface in Fig. 3) was validated by RNA gel analysis. To examine whether the response of these miRNAs was specific to −Pi or not, the RNA isolated from different nutrient deficiencies, including N, Pi, potassium (K), S, Cu, and iron (Fe), was inspected at once (Fig. 4 Figure 4. Open in new tabDownload slide Validation of differentially expressed miRNAs under different nutrient deficiencies. Small RNA gel analyses of Pi starvation-up-regulated miRNAs (A) and Pi starvation-down-regulated miRNAs (B) are shown. miR778 was detected using the probe corresponding to miR778-3p.1 and miR778-3p.2. Figure 4. Open in new tabDownload slide Validation of differentially expressed miRNAs under different nutrient deficiencies. Small RNA gel analyses of Pi starvation-up-regulated miRNAs (A) and Pi starvation-down-regulated miRNAs (B) are shown. miR778 was detected using the probe corresponding to miR778-3p.1 and miR778-3p.2. ). miR156 was highly expressed but unaffected by nutrient deficiencies in the shoots; however, it was up-regulated by −Pi, −N, or −K, with the highest induction under −N. miR156 regulates a group of SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE (SPL) transcription factors involved in juvenile-to-adult vegetative phase transition and other development processes (Wu and Poethig, 2006; Gandikota et al., 2007; Schwarz et al., 2008; Wang et al., 2008). The physiological role for the up-regulation of miR156 in the −Pi, −N, or −K root remains to be explored. On the other hand, miR399, miR778, and miR827 were specifically up-regulated by −Pi (Fig. 4A). It is interesting that several distinct small RNA sequences were generated by cleavage at different positions of the miR778 precursor (Supplemental Fig. S1). Because the corresponding miRNA* sequence could be detected for each cleaved RNA duplex, we suggest to rename miR778 as miR778-5p and the other two forms derived from the miR778 precursor as miR778-3p.1 and miR778-3p.2, respectively. Multiple processing of the miR778 precursor has not been reported previously, probably due to the low recovered sequence reads. All miR778 species were up-regulated by −Pi, with the highest sequence count of miR778-3p.1 (Supplemental Table S4). The target genes of miR778-5p encode SET domain-containing proteins, which are involved in the regulation of histone methylation (Ebbs and Bender, 2006). Although the target genes of miR778-3p.1 and miR778-3p.2 await identification, generation of multiple miRNAs from the same precursor has increased the complexity for the miR778 regulation networks. The target genes of miR827 encode SPX domain-containing proteins. In yeast cells, SPX domain-containing proteins participate in Pi transport or sensing (Lenburg and O'Shea, 1996). One of the target genes, At1g02860, encodes a nucleus-localized ubiquitin E3 ligase containing the RING and SPX domains, which was recently characterized as NITROGEN LIMITATION ADAPTATION (NLA) involved in the N recycle during N limitation (Peng et al., 2007) or as BENZOIC ACID HYPERSENSTIVE1 (BAH1) involved in immune responses (Yaeno and Iba, 2008). Here, we validated the targeting of miR827 to At1g02860 mRNA by RNA ligase-mediated (RLM) 5′-RACE analysis (Supplemental Fig. S2A) and observed the down-regulation of At1g02860 under −Pi (Supplemental Fig. S2B). At1g02860 may be regulated by many factors and have multiple functions in response to diverse external stimuli. The miR827-mediated regulation of At1g02860 under −Pi conditions may imply an interaction between N and P, although miR827 only responds to −Pi but not −N (Fig. 4A). The expression of several down-regulated miRNAs was also confirmed (Fig. 4B). miR169 was down-regulated in −Pi, −N, or −S roots but to a lesser extent in the shoots. Down-regulation of miR169 was much more severe under −N. miR395 was down-regulated in −Pi or −N roots and shoots. miR398 was down-regulated by −Pi, −N, −K, or −Fe in both roots and shoots. In contrast, miR395 and miR398 were up-regulated by −S and −Cu, respectively, which is in agreement with previous findings (Jones-Rhoades and Bartel, 2004; Yamasaki et al., 2007). These observations suggest that miRNAs can coordinate cross talk among different nutrient-deficient responses. Noticing the up-regulation of APS4 and SULTR2;1 under −Pi (Supplemental Fig. S2C), we speculate that the suppression of miR395 under −Pi may somewhat contribute to this effect. Up-regulation of APS4 and SULTR2;1 may result in the increase of sulfate assimilation and translocation for sulfolipid biosynthesis, a compensation for the reduced phospholipids during −Pi. miR169 and its target gene, HAP2/NFYA transcription factor, are critical in the development of nodules (Combier et al., 2006). Additionally, miR169a and NFYA5 have been shown to be involved in drought tolerance (Li et al., 2008a). NFYA5 functions upstream to activate several drought-responsive genes, including those involved in oxidative stress. The target genes of miR398 encoding Cu/Zn superoxide dismutase are also involved in the response to oxidative stress (Sunkar et al., 2006). Several nutrient deficiencies, such as N, Pi, K, and S, can trigger oxidative stresses (Shin et al., 2005; Schachtman and Shin, 2007). Changes of miR169 and miR398 expression could be involved in adjusting a plant's ability to survive the oxidative stress caused by nutrient deficiencies. Interestingly, the opposite expression of miR398 in response to Cu and Fe deficiencies may reflect an elegant regulation of the preferential expression of different superoxide dismutase, Cu/Zn-superoxide dismutase or Fe-superoxide dismutase, coupling with different metal cofactors. These Pi-responsive miRNAs and the potential functions of corresponding target genes relative to adaptive responses to −Pi are summarized in Table I Table I. Summary of validated Pi-responsive miRNAs and ta-siRNAs and the potential function involved under Pi deficiency miRNA/ta-siRNA . . . Target Genes . . . Familya . Speciesb . Expressionc . Gene Productd . Expressione . Potential Rolef . miR169 a, b, c, d, e, f, g, h, i, j, k, l, m, n Down HAP2 transcription factors (At1g17590, At1g54160, At1g72830, At3g05690, At3g20910, At5g06510, At5g12840) Up Involved in antioxidative stress miR395 b, c, f Down ATP sulfurylase (APS1, APS3, APS4; At3g22890, At4g14680, At5g43780) Up Increase sulfate translocation and utilization for sulfolipid biosynthesis Sulfate transporter (SULTR2;1; At5g10180) Up miR398 a, b, c Down Cu/Zn superoxide dismutase (CSD; CSD1, At1g08830; CSD2, At2g28190) Unchanged Scavenge reactive oxygen species Cytochrome c oxidase (At3g15640) Unchanged miR402 Downg DEMETER-LIKE protein 3 (DML3; At4g34060) Undetermined miR156 a, b, c, d, e, f, g Up (root) SPL transcription factors (At1g53160, At2g33810, At3g15270, At5g43270, At1g27360, At1g27370, At1g69170, At2g42200, At3g57920, At5g50570, At5g50670) Down miR399 a, b, c, d, f Up Ubiquitin conjugase E2 (UBC24/PHO2; At2g33770) Down Pi uptake and translocation miR778 5p Up SET domain-containing protein (At2g22740, At2g35160) Down 3p.1 Up 3p.2 Up miR827b Up Ubiquitin E3 ligase with RING and SPX domains (NLA/BAH1; At1g02860) SPX domain-containing protein (At1g63010) Down Unchanged Nutrient recycle miR2111 a-5p, b-5p Up Kelch repeat-containing F-box protein (At3g27150) Up a-3p, b-3p Up miR828 Up (shoot) TAS4 Anthocyanin biosynthesis TAS4-siRNAb D1(+), D2(+), D3(+) TAS4-siR81(−) Up (shoot) Up (shoot) MYB75, MYB90, and MYB113 transcription factors (At1g56650, At1g66390, At1g66370) Up Anthocyanin biosynthesis miRNA/ta-siRNA . . . Target Genes . . . Familya . Speciesb . Expressionc . Gene Productd . Expressione . Potential Rolef . miR169 a, b, c, d, e, f, g, h, i, j, k, l, m, n Down HAP2 transcription factors (At1g17590, At1g54160, At1g72830, At3g05690, At3g20910, At5g06510, At5g12840) Up Involved in antioxidative stress miR395 b, c, f Down ATP sulfurylase (APS1, APS3, APS4; At3g22890, At4g14680, At5g43780) Up Increase sulfate translocation and utilization for sulfolipid biosynthesis Sulfate transporter (SULTR2;1; At5g10180) Up miR398 a, b, c Down Cu/Zn superoxide dismutase (CSD; CSD1, At1g08830; CSD2, At2g28190) Unchanged Scavenge reactive oxygen species Cytochrome c oxidase (At3g15640) Unchanged miR402 Downg DEMETER-LIKE protein 3 (DML3; At4g34060) Undetermined miR156 a, b, c, d, e, f, g Up (root) SPL transcription factors (At1g53160, At2g33810, At3g15270, At5g43270, At1g27360, At1g27370, At1g69170, At2g42200, At3g57920, At5g50570, At5g50670) Down miR399 a, b, c, d, f Up Ubiquitin conjugase E2 (UBC24/PHO2; At2g33770) Down Pi uptake and translocation miR778 5p Up SET domain-containing protein (At2g22740, At2g35160) Down 3p.1 Up 3p.2 Up miR827b Up Ubiquitin E3 ligase with RING and SPX domains (NLA/BAH1; At1g02860) SPX domain-containing protein (At1g63010) Down Unchanged Nutrient recycle miR2111 a-5p, b-5p Up Kelch repeat-containing F-box protein (At3g27150) Up a-3p, b-3p Up miR828 Up (shoot) TAS4 Anthocyanin biosynthesis TAS4-siRNAb D1(+), D2(+), D3(+) TAS4-siR81(−) Up (shoot) Up (shoot) MYB75, MYB90, and MYB113 transcription factors (At1g56650, At1g66390, At1g66370) Up Anthocyanin biosynthesis a The Pi-responsive miRNAs reported by Pant et al. (2009) are underlined. b Only those species with statistically significant changes in expression are shown. Boldface letters indicate the existence of the P1BS cis-element (GNATATNC), a binding site for PHR1, in the 2-kb sequences upstream of miRNA hairpins. There are P1BS elements in the upstream sequences of MIR827 and TAS4 genes. c Change of small RNA expression due to Pi deficiency. d The target genes are based on the Arabidopsis Small RNA Project database (http://asrp.cgrb.oregonstate.edu/) except miR2111 (this study). The target genes with altered expression by Pi deficiency are indicated in boldface. e Change of target gene expression to Pi deficiency. The results were based on the microarry analysis from Genevestigator (https://www.genevestigator.ethz.ch/gv/index.jsp) except for the target genes of miR399 (Chiou et al., 2006), miR2111-5p and TAS4-siR81(−) (Figs. 5 and 7), miR395, and miR827 (Supplemental Fig. S2), whose expression was analyzed by quantitative RT-PCR. Changes larger than 2-fold are indicated. f Potential role of target gene under Pi deficiency. g Because of its low abundance, the expression of miR402 was validated by quantitative RT-PCR of primary transcript. Open in new tab Table I. Summary of validated Pi-responsive miRNAs and ta-siRNAs and the potential function involved under Pi deficiency miRNA/ta-siRNA . . . Target Genes . . . Familya . Speciesb . Expressionc . Gene Productd . Expressione . Potential Rolef . miR169 a, b, c, d, e, f, g, h, i, j, k, l, m, n Down HAP2 transcription factors (At1g17590, At1g54160, At1g72830, At3g05690, At3g20910, At5g06510, At5g12840) Up Involved in antioxidative stress miR395 b, c, f Down ATP sulfurylase (APS1, APS3, APS4; At3g22890, At4g14680, At5g43780) Up Increase sulfate translocation and utilization for sulfolipid biosynthesis Sulfate transporter (SULTR2;1; At5g10180) Up miR398 a, b, c Down Cu/Zn superoxide dismutase (CSD; CSD1, At1g08830; CSD2, At2g28190) Unchanged Scavenge reactive oxygen species Cytochrome c oxidase (At3g15640) Unchanged miR402 Downg DEMETER-LIKE protein 3 (DML3; At4g34060) Undetermined miR156 a, b, c, d, e, f, g Up (root) SPL transcription factors (At1g53160, At2g33810, At3g15270, At5g43270, At1g27360, At1g27370, At1g69170, At2g42200, At3g57920, At5g50570, At5g50670) Down miR399 a, b, c, d, f Up Ubiquitin conjugase E2 (UBC24/PHO2; At2g33770) Down Pi uptake and translocation miR778 5p Up SET domain-containing protein (At2g22740, At2g35160) Down 3p.1 Up 3p.2 Up miR827b Up Ubiquitin E3 ligase with RING and SPX domains (NLA/BAH1; At1g02860) SPX domain-containing protein (At1g63010) Down Unchanged Nutrient recycle miR2111 a-5p, b-5p Up Kelch repeat-containing F-box protein (At3g27150) Up a-3p, b-3p Up miR828 Up (shoot) TAS4 Anthocyanin biosynthesis TAS4-siRNAb D1(+), D2(+), D3(+) TAS4-siR81(−) Up (shoot) Up (shoot) MYB75, MYB90, and MYB113 transcription factors (At1g56650, At1g66390, At1g66370) Up Anthocyanin biosynthesis miRNA/ta-siRNA . . . Target Genes . . . Familya . Speciesb . Expressionc . Gene Productd . Expressione . Potential Rolef . miR169 a, b, c, d, e, f, g, h, i, j, k, l, m, n Down HAP2 transcription factors (At1g17590, At1g54160, At1g72830, At3g05690, At3g20910, At5g06510, At5g12840) Up Involved in antioxidative stress miR395 b, c, f Down ATP sulfurylase (APS1, APS3, APS4; At3g22890, At4g14680, At5g43780) Up Increase sulfate translocation and utilization for sulfolipid biosynthesis Sulfate transporter (SULTR2;1; At5g10180) Up miR398 a, b, c Down Cu/Zn superoxide dismutase (CSD; CSD1, At1g08830; CSD2, At2g28190) Unchanged Scavenge reactive oxygen species Cytochrome c oxidase (At3g15640) Unchanged miR402 Downg DEMETER-LIKE protein 3 (DML3; At4g34060) Undetermined miR156 a, b, c, d, e, f, g Up (root) SPL transcription factors (At1g53160, At2g33810, At3g15270, At5g43270, At1g27360, At1g27370, At1g69170, At2g42200, At3g57920, At5g50570, At5g50670) Down miR399 a, b, c, d, f Up Ubiquitin conjugase E2 (UBC24/PHO2; At2g33770) Down Pi uptake and translocation miR778 5p Up SET domain-containing protein (At2g22740, At2g35160) Down 3p.1 Up 3p.2 Up miR827b Up Ubiquitin E3 ligase with RING and SPX domains (NLA/BAH1; At1g02860) SPX domain-containing protein (At1g63010) Down Unchanged Nutrient recycle miR2111 a-5p, b-5p Up Kelch repeat-containing F-box protein (At3g27150) Up a-3p, b-3p Up miR828 Up (shoot) TAS4 Anthocyanin biosynthesis TAS4-siRNAb D1(+), D2(+), D3(+) TAS4-siR81(−) Up (shoot) Up (shoot) MYB75, MYB90, and MYB113 transcription factors (At1g56650, At1g66390, At1g66370) Up Anthocyanin biosynthesis a The Pi-responsive miRNAs reported by Pant et al. (2009) are underlined. b Only those species with statistically significant changes in expression are shown. Boldface letters indicate the existence of the P1BS cis-element (GNATATNC), a binding site for PHR1, in the 2-kb sequences upstream of miRNA hairpins. There are P1BS elements in the upstream sequences of MIR827 and TAS4 genes. c Change of small RNA expression due to Pi deficiency. d The target genes are based on the Arabidopsis Small RNA Project database (http://asrp.cgrb.oregonstate.edu/) except miR2111 (this study). The target genes with altered expression by Pi deficiency are indicated in boldface. e Change of target gene expression to Pi deficiency. The results were based on the microarry analysis from Genevestigator (https://www.genevestigator.ethz.ch/gv/index.jsp) except for the target genes of miR399 (Chiou et al., 2006), miR2111-5p and TAS4-siR81(−) (Figs. 5 and 7), miR395, and miR827 (Supplemental Fig. S2), whose expression was analyzed by quantitative RT-PCR. Changes larger than 2-fold are indicated. f Potential role of target gene under Pi deficiency. g Because of its low abundance, the expression of miR402 was validated by quantitative RT-PCR of primary transcript. Open in new tab . Up-regulation of miR399, miR778, and miR827 and down-regulation of miR398 by −Pi were also observed in a recent study in which a quantitative RT-PCR platform was established to profile the primary transcripts of miRNAs (Pant et al., 2009). The changes of these primary transcripts of miRNAs agreed with those of their mature miRNAs in response to −Pi, suggesting that they are mainly regulated at the transcriptional level. However, changes in miR156 and miR395 were not discovered in their system, probably because of differences in growth conditions or sample sources (Table I). To examine whether the regulation of these miRNAs depends on the PHOSPHATE STARVATION RESPONSE1 (PHR1) MYB transcription factor, which activates a group of Pi starvation-induced genes, including MIR399s (Rubio et al., 2001; Bari et al., 2006), we analyzed the existence of P1BS cis-element, the biding site of PHR1, within the 2-kb sequences upstream of miRNA hairpins. In addition to MIR399s, the P1BS element could be found in the upstream sequences of several species of MIR156, MIR169, MIR398c, and MIR827 genes (Table I), suggesting the possible regulation of these genes by PHR1. Characterization of the Pi-Responsive miRNA2111 Family We identified a miRNA gene family consisting of two genes (miR2111a and miR2111b) from the nonannotated sequence clusters, which were not reported as miRNAs during the time of our analysis. They were specifically up-regulated by −Pi but not by any of the other nutrient deficiencies tested (Fig. 5A Figure 5. Open in new tabDownload slide Characterization of miR2111 in response to Pi deficiency. A and B, Small RNA gel analyses of miR2111-5p and miR2111-3p under different nutrient-deficient conditions (A) and in small RNA biogenesis mutants (B). Col, Ecotype Columbia; Ler, ecotype Landsberg erecta; wt, wild type. C, Hairpin secondary structures of miR2111a and miR2111b precursors. D, RLM 5′-RACE analysis of miR2111-5p cleavage on At3g27150 mRNA. E, Quantitative RT-PCR analysis of At3g27150 mRNA in response to −Pi. One of two biological replicates is presented, and the error bars indicate the sd of two technical replicates. Figure 5. Open in new tabDownload slide Characterization of miR2111 in response to Pi deficiency. A and B, Small RNA gel analyses of miR2111-5p and miR2111-3p under different nutrient-deficient conditions (A) and in small RNA biogenesis mutants (B). Col, Ecotype Columbia; Ler, ecotype Landsberg erecta; wt, wild type. C, Hairpin secondary structures of miR2111a and miR2111b precursors. D, RLM 5′-RACE analysis of miR2111-5p cleavage on At3g27150 mRNA. E, Quantitative RT-PCR analysis of At3g27150 mRNA in response to −Pi. One of two biological replicates is presented, and the error bars indicate the sd of two technical replicates. ). They possessed a potential hairpin secondary structure and could produce a small RNA duplex with two-nucleotide overhang at the 3′ end (Fig. 5C), the primary criterion of miRNA definition (Meyers et al., 2008). Because a substantial amount of both complementary strands of the RNA duplex was detected, it was difficult to differentiate miRNA from miRNA*; thus, we named them miR2111a-5p, miR2111a-3p, miR2111b-5p, and miR2111b-3p. The sequences of miR2111a-5p and miR2111b-5p were identical, while those of miR2111a-3p and miR2111b-3p differed by four nucleotides. Although this family was also identified by other groups very recently (Fahlgren et al., 2009; Pant et al., 2009), its biogenesis and target genes have not been studied. The biogenesis of these miRNAs was analyzed using various mutants defective in the production of different classes of small RNAs (Fig. 5B). Here, we showed that the accumulation of these miRNAs was significantly reduced in the Pi-starved hen1 and dcl1 mutants. Nevertheless, some accumulation of miR2111 was still observed in the dcl1 mutant, suggesting that miR2111 can also be processed by another DCL, probably DCL4, in the absence of DCL1. The gene At3g27150, encoding a Kelch domain-containing F-box protein, was predicted to be a target of miR2111-5p. The cleavage of the mRNA of At3g27150 was confirmed by RLM 5′-RACE analysis (Fig. 5D). Interestingly, we did not observe the negative correlation between the expression of miR2111-5p and At3g27150. Instead, the steady-state level of At3g27150 mRNA was moderately induced by −Pi (Fig. 5E) regardless of the enhanced expression of miR2111-5p. Several reports have shown a positive temporal correlation between the expression of a miRNA and its target gene in spite of the negative regulatory role of miRNAs (Llave et al., 2002; Achard et al., 2004; Kawashima et al., 2009). The consequences of different miRNA regulatory circuits were classified into four categories, spatial restriction, mutual exclusion, dampening, and temporal regulation, in which the final expression level of the miRNA and its target gene may not be correlated (Voinnet, 2009). We hypothesize that miR2111-5p may regulate the spatial expression of At3g27150 by restricting its expression in specific cells and/or fine-tuning the expression of At3g27150 to maintain an optimal expression level. Alternatively, the possibility of translational repression cannot be ruled out. Nevertheless, regulation by miR2111 may be only one of many regulatory routes of At3g27150, because multiple regulatory pathways integrated from different levels are usually involved in gene regulation. It is important to note that the target genes of miR399, miR827, and miR2111-5p encode proteins involved in ubiquitin-mediated protein degradation, suggesting a crucial role of posttranslational regulation of protein abundance in the control of −Pi responses. The homologs of the miR2111 gene could be identified in many other dicotyledon plant species, such as Brassica rapa, M. truncatula, Lotus japonicus, and Vitis vinifera (Supplemental Fig. S3), but they could not be found in monocotyledon plants. This observation suggests that either the miR2111 gene had become lost in monocotyledons or, more likely, that it had emerged after the separation of the dicotyledons from the monocotyledons. The sequence of miR2111-5p is more conserved than that of miR2111-3p. The four miR2111 genes in M. truncatula are clustered in the same direction within a range of 40 kb. Interestingly, there are only 32 nucleotides between the middle two hairpins, which are likely derived from a polycistronic transcript. Conservation of the miR2111 species implies the importance of miR2111-mediated regulation for survival under the −Pi environment. Pi-Responsive Small RNAs Derived from LTR of Retrotransposon A cluster of small RNAs that originated from the long terminal repeat (LTR) of Copia95 retrotransposon was identified as showing a differential expression pattern in response to −Pi (Fig. 6A Figure 6. Open in new tabDownload slide Pi starvation-induced small RNAs derived from LTR of retrotransposon. A and B, Small RNA gel analysis of smRPi1LTR under different nutrient-deficient conditions (A) and in small RNA biogenesis mutants (B). Col, Ecotype Columbia; hetero, heterozygous plants; wt, wild type. C, Arrangement of Copia95 LTR in the intergenic region of At5g27990 and At5g28000 and the potential hairpin structure and location of smRPi1LTR. Figure 6. Open in new tabDownload slide Pi starvation-induced small RNAs derived from LTR of retrotransposon. A and B, Small RNA gel analysis of smRPi1LTR under different nutrient-deficient conditions (A) and in small RNA biogenesis mutants (B). Col, Ecotype Columbia; hetero, heterozygous plants; wt, wild type. C, Arrangement of Copia95 LTR in the intergenic region of At5g27990 and At5g28000 and the potential hairpin structure and location of smRPi1LTR. ). This small RNA cluster is located in the intergenic region of At5g27990 and At5g28000 containing a multiple array of truncated LTRs either as inverted or tandem repeats, suggesting the occurrence of several duplication events (Fig. 6C). The arrangement of inverted repeats allows the formation of a hairpin structure. The length of this hairpin is relatively long, at approximately 250 nucleotides. Unlike those miRNA-like small RNAs recently found in rice (Oryza sativa) grains (Heisel et al., 2008; Zhu et al., 2008), which are cleaved from a long hairpin, the small RNAs generated from this LTR hairpin are not phased. smRPi1LTR, the most abundant small RNA from the hairpin stem, accumulated specifically in −Pi roots but not in other nutrient-deficient tissues (Fig. 6A). DCL1 is essential for the biogenesis of smRPi1LTR, and its accumulation was not reduced in the dcl2,3,4 triple mutants or in the rdr1, rdr2, and rdr6 mutants (Fig. 6B). smRPi1LTR is likely generated from the cleavage of a single-stranded rather than a double-stranded RNA precursor, because the passenger strand of its duplex was identified from the hairpin stem with two mismatched base pairs. Because we were unable to detect smRPi1LTR in the Landsberg accession, we hypothesize that smRPi1LTR is a newly evolved small RNA due to rapid rearrangement of LTR and may represent an intermediate small RNA species from the transition of siRNAs to miRNAs. AtCopeg1 (for Copia evolved gene1; At2g04460), a protein-coding gene originating from a truncated AtCopia95 retrotransposon, was recently shown to be highly up-regulated by −Pi (Duan et al., 2008). It will be interesting to elucidate the relationship between AtCopeg1 and smRPi1LTR in terms of coevolution or cross-regulation. Up-Regulation of miR828- and TAS4-Derived siRNAs by Pi Deficiency and Autoregulation of PAP1/MYB75 The sequencing results revealed the up-regulation of ta-siRNAs derived from TAS4 and miR828, which initiates the cleavage of TAS4 RNA in −Pi shoots (Supplemental Table S4). These ta-siRNAs were up-regulated in the shoots by −Pi as well as −N (Fig. 7A Figure 7. Open in new tabDownload slide Autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(−). A, Small RNA gel analysis of TAS4-siRNAs under different nutrient-deficient conditions. B, Small RNA gel analysis of TAS4-siR81(−) and miR828 in the shoots of wild-type (wt) plants and tas4, mir828, and PAP1/MYB75 activation (pap1-D) mutants. C, Expression of PAP1/MYB75, PAP2/MYB90, and MYB113 by quantitative RT-PCR analysis in tas4 and mir828 mutants. One of two biological replicates is presented, and the error bars indicate the sd of two technical replicates. D, Alteration of anthocyanin accumulation in tas4, mir828, and pap1-D mutants. Error bars represent the sd (n = 5). FW, Fresh weight; OD, optical density. E, A proposed model of the autoregulatory pathway of PAP1/MYB75 via miR828 and TAS4-siR81(−). Figure 7. Open in new tabDownload slide Autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(−). A, Small RNA gel analysis of TAS4-siRNAs under different nutrient-deficient conditions. B, Small RNA gel analysis of TAS4-siR81(−) and miR828 in the shoots of wild-type (wt) plants and tas4, mir828, and PAP1/MYB75 activation (pap1-D) mutants. C, Expression of PAP1/MYB75, PAP2/MYB90, and MYB113 by quantitative RT-PCR analysis in tas4 and mir828 mutants. One of two biological replicates is presented, and the error bars indicate the sd of two technical replicates. D, Alteration of anthocyanin accumulation in tas4, mir828, and pap1-D mutants. Error bars represent the sd (n = 5). FW, Fresh weight; OD, optical density. E, A proposed model of the autoregulatory pathway of PAP1/MYB75 via miR828 and TAS4-siR81(−). ). TAS4-siR81(−) targets a group of MYB transcription factors, PAP1/MYB75, PAP2/MYB90, and MYB113, that are involved in the biosynthesis of anthocyanin (Rajagopalan et al., 2006). It is well documented that plants accumulate anthocyanin under −Pi, in which these MYB transcription factors are up-regulated (Misson et al., 2005). The negative regulatory role of siRNA seemed to be contradictory to the positive correlation between the expression of TAS4-siR81(−) and its target genes under the same conditions. This prompted us to carry out further investigation. TAS4 and MIR828 T-DNA knockout lines and a PAP1/MYB75 activation mutant, production of anthocyanin pigment1-Dominant (pap1-D; Borevitz et al., 2000), were examined. As predicted, there was no detectable TAS4-siR81(−) in the tas4 mutant, and neither miR828 nor TAS4-siR81(−) was detected in the mir828 mutant (Fig. 7B). The accumulation of miR828 and TAS4-siR81(−) was already observed in the +Pi shoots of pap1-D mutant and increased under −Pi (Fig. 7B). This suggests that PAP1/MYB75 can positively regulate MIR828 and/or TAS4 genes either by activating their expression directly or indirectly or by suppressing a negative regulator of MIR828 and/or TAS4 via a different feedback loop. Analyses of the 1-kb upstream sequences of MIR828 and TAS4 genes revealed multiple MYB binding sites (predicted by the PLACE database; http://www.dna.affrc.go.jp/PLACE/index.html) and one PAP1 cis-regulatory element, (C/T)(A/C)NCCAC(A/G/T)N(G/T) (Dare et al., 2008), at 305 and 747 bp upstream of the MIR828 hairpin precursor and TAS4, respectively. Thus, PAP1 likely acts directly on the MIR828 and TAS4 promoters to activate their transcription. On the other hand, the transcript level of PAP1/MYB75, PAP2/MYB90, and MYB113 was elevated in the tas4 and mir828 mutants under both +Pi and −Pi conditions (Fig. 7C). This could be explained by the relief of targeting of TAS4-siR81(−). Consistent with the increased expression of these MYB transcription factors, the anthocyanin accumulation was enhanced in the shoots of tas4 and mir828 mutants under both +Pi and −Pi conditions (Fig. 7D). These results uncovered an autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(−) in which an adequate expression level of PAP1/MYB75, PAP2/MYB90, and MYB113 and a proper accumulation of anthocyanin are maintained during −Pi (Fig. 7E). In this model, Pi deficiency results in up-regulation of those MYB transcription factors, which subsequently activate the biosynthesis of anthocyanin. On the other hand, the increased level of PAP1/MYB75 can also lead to production of TAS4-siR81(−) via the activation of miR828 and/or TAS4, which antagonizes its own expression and other MYB transcription factors. Whether the up-regulation of miR828 can be bypassed by the activation of PAP1/MYB75 remains to be studied. Because anthocyanin is a common stress pigment, we hypothesize that this autoregulation machinery may not be restricted to −Pi and could also be applied to other stress conditions in which anthocyanin is accumulated. miRNA399* and Its Movement In addition to miR399s, we observed a considerable amount of miR399* accumulated in −Pi tissues. The read number of different miR399* species was about half that of miR399 (e.g. miR399f*) or even higher than that of miR399 (e.g. miR399a* and miR399d*). Such high accumulation of miR399* could not be neglected because it was more abundant than many other miRNAs. RNA gel analyses confirmed the existence of miR399* with the high abundance in miR399d* and miR399f*, predominantly in the roots (Fig. 8A Figure 8. Open in new tabDownload slide Expression and long-distance movement of miR399*. A, Small RNA gel analysis of miR399*. B, Quantitative RT-PCR analysis of miR399f and miR399f* in the rootstocks of grafted plants. Combination of grafts is indicated as scion/rootstock. 399f, miR399f-overexpressing plants; wt, wild-type plants. The error bars represent sd (n = 4–5 independent grafted plants). Figure 8. Open in new tabDownload slide Expression and long-distance movement of miR399*. A, Small RNA gel analysis of miR399*. B, Quantitative RT-PCR analysis of miR399f and miR399f* in the rootstocks of grafted plants. Combination of grafts is indicated as scion/rootstock. 399f, miR399f-overexpressing plants; wt, wild-type plants. The error bars represent sd (n = 4–5 independent grafted plants). ), which was consistent with the high level of their primary transcripts and of mature miR399d and miR399f (Lin et al., 2008). We and others have previously shown that miR399 can move from shoots to roots and function as a long-distance signal in response to −Pi (Lin et al., 2008; Pant et al., 2008). We found that like miR399f, miR399f* could be detected in the wild-type rootstock grafted with the miR399f-overexpressing scions (Fig. 8B). This result indicates that both miR399f and miR399f* can move across the graft junction from scions to rootstocks. We hypothesize that miR399* may have certain physiological relevance, possibly functioning in concert with miR399. There may be several possible functions of miR399*. First, miR399* may be loaded into the RNA-induced silencing complex to target different genes. Many miRNA* species from Drosophila were reported to associate with the Argonaute protein and have inhibitory activity (Okamura et al., 2008). Several genes were predicted to be the putative targets of miR399* species; however, we were unable to validate them by RLM 5′-RACE analysis. Thus, whether these miR399* species possess gene-silencing activity remains uncertain. Second, given the ability of miR399 and miR399* to move across the graft junction, we suspect that miR399* may assist the long-distance movement of miR399 by forming a miR399/miR399* duplex in the phloem stream. Several miRNA* species were also detected in the phloem saps of Brassica napus (Buhtz et al., 2008; Pant et al., 2009). Although the authors concluded that they remain in the single-stranded form, based on the in vitro RNase digestion assay, they could not exclude the possibility of separation from a double-stranded RNA duplex during the RNA extraction. An in vivo nondestructive analysis needs to be developed to differentiate the double-stranded small RNA molecules from the single-stranded molecules. Alternatively, miR399* may regulate or buffer the targeting efficiency of miR399 on PHO2 mRNA via the formation of an RNA duplex. CONCLUSION In summary, we have presented an extensive survey of the small RNAs showing differential expression in response to −Pi. These Pi-responsive small RNAs and their target genes are likely involved in the development or regulation of adaptive responses to −Pi. Interestingly, while miR399, miR778, miR827, and miR2111 responded specifically to −Pi, miR156, miR169, miR395, and miR398 also responded to other nutrient stresses. These responses involve up-regulation or down-regulation of distinct miRNAs, which coordinate and balance the demands of different nutrients. We hypothesize that metabolic adjustment, carbon assimilation, or oxidative stresses may have an effect on the cross talk among different nutrients, as we observed here. Although there is still much more to learn about the biological significances of these Pi-responsive small RNAs, this work has opened a new avenue for the functional study of small RNA-mediated gene regulation in −Pi responses. MATERIALS AND METHODS Plant Materials and Growth Conditions Arabidopsis (Arabidopsis thaliana ecotype Columbia) wild-type plants and mutants (mir828 [SALK_097788], tas4 [SALK_066997], and pap1-D [CS3884]; Borevitz et al., 2000; Alonso et al., 2003) were grown hydroponically as described (Aung et al., 2006; Chiou et al., 2006). Seventeen-day-old seedlings were subjected to −Pi treatment by transferring them from +Pi (250 μm KH2PO4) into Pi-free half-strength Hoagland solution. Roots and shoots were collected separately after 5 d of treatment. Control samples were maintained in +Pi and collected at the same time. For treatments with different nutrient deficiency, the specific nutrient element was withdrawn from the nutrient solution for 5 d. To compensate for the ion originally associated with the omitted nutrient element, it was supplied in the form of chlorite salt at equivalent concentration except in the −Cu and −Fe treatments, in which they were omitted from the solution. RNA Preparation and Detection and High-Throughput Small RNA Sequencing Total RNA was isolated using Trizol reagent (Invitrogen). Four RNA samples from +Pi or −Pi roots or shoots of wild-type plants were used to generate small RNA libraries by the Solexa sequencing technology (Illumina). Small RNA gel-blot analysis was carried out as described (Chiou et al., 2006) with modification according to Pall et al. (2007). Quantitative RT-PCR analysis was performed as described (Lin et al., 2008). The sequences of probes and PCR primers are listed in Supplemental Table S5. Analysis of Sequencing Data The identical adapter-trimmed reads were grouped as unique sequences with associated counts of the individual sequence reads. Each unique sequence was aligned to the Arabidopsis genome (The Arabidopsis Information Resource 8; http://www.arabidopsis.org/), and only perfectly mapped sequences were retained and analyzed further (for details, see Supplemental Materials and Methods S1). We defined a sequence cluster as a group of same-strand sequences that were mapped to similar genomic locations, and the union of the mapped regions of the individual sequences of the cluster can form a consecutive nongapped region. For every cluster, the number of sequence reads from each of the four samples was counted separately and used for differential expression analysis. We used two approaches to perform the differential expression analysis, a cluster-based approach for finding novel differentially expressed gene models and a gene-model-based approach for analyzing known gene models, especially for miRNA genes and ta-siRNA genes. All sequence clusters or genes of interest were compared between two libraries for obtaining differentially expressed ones using a statistics test. For each of the differentially expressed sequence clusters that were not annotated, the sequence including the cluster union and 200-bp flanking genomic sequence extending from each side was extracted and folded by RNAfold (Hofacker, 2003). A sequence cluster was considered to be a novel miRNA if it met the following four criteria: (1) the structure forms a hairpin; (2) the sequence cluster is located in the duplex region of the hairpin structure; (3) there is another sequence cluster located on the other strand of the duplex; and (4) the most abundant sequences from each of the two clusters can form a duplex with two-nucleotide overhang at the 3′ end. The targets for the novel miRNAs were predicted using methods described elsewhere (Allen et al., 2005). Validation of the target gene was carried out by RLM 5′-RACE analysis using the GeneRacer kit (Invitrogen). Anthocyanin Analysis Anthocyanin analysis was performed according to Lange et al. (1971). The small RNA sequence data were submitted to the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/projects/geo/) under accession number GSE17741. Accession numbers of miRNAs are available at miRBase (http://microrna.sanger.ac.uk/sequences/index.shtml). Supplemental Data The following materials are available in the online version of this article. Supplemental Figure S1. Multiple processing of the miR778 precursor. Supplemental Figure S2. Expression of the target genes of miR395 and miR827. Supplemental Figure S3. Evolutionary analysis of MIR2111. Supplemental Table S1. Statistics of small RNA sequences. Supplemental Table S2. Categories of small RNA origins. Supplemental Table S3. tRNA-derived small RNAs. Supplemental Table S4. Expression of miRNA, ta-siRNA, and smRPi1LTR. Supplemental Table S5. Sequences of probes and primers. Supplemental Materials and Methods S1. ACKNOWLEDGMENTS We thank Shu-Hsing Wu and Ho-Ming Chen for helpful comments and Tzu-Yun Chang for assisting in the preparation of materials. LITERATURE CITED Achard P, Herr A, Baulcombe DC, Harberd NP ( 2004 ) Modulation of floral development by a gibberellin-regulated microRNA. Development 131 : 3357 – 3365 Crossref Search ADS PubMed Allen E, Xie Z, Gustafson AM, Carrington JC ( 2005 ) MicroRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell 121 : 207 – 221 Crossref Search ADS PubMed Alonso JM, Stepanova AN, Leisse TJ, Kim CJ, Chen H, Shinn P, Stevenson DK, Zimmerman J, Barajas P, Cheuk R, et al ( 2003 ) Genome-wide insertional mutagenesis of Arabidopsis thaliana. Science 301 : 653 – 657 Crossref Search ADS PubMed Aung K, Lin SI, Wu CC, Huang YT, Su CL, Chiou TJ ( 2006 ) pho2, a phosphate overaccumulator, is caused by a nonsense mutation in a microRNA399 target gene. Plant Physiol 141 : 1000 – 1011 Crossref Search ADS PubMed Bari R, Pant BD, Stitt M, Scheible WR ( 2006 ) PHO2, microRNA399, and PHR1 define a phosphate-signaling pathway in plants. Plant Physiol 141 : 988 – 999 Crossref Search ADS PubMed Borevitz JO, Xia Y, Blount J, Dixon RA, Lamb C ( 2000 ) Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell 12 : 2383 – 2394 Crossref Search ADS PubMed Brodersen P, Sakvarelidze-Achard L, Bruun-Rasmussen M, Dunoyer P, Yamamoto YY, Sieburth L, Voinnet O ( 2008 ) Widespread translational inhibition by plant miRNAs and siRNAs. Science 320 : 1185 – 1190 Crossref Search ADS PubMed Brodersen P, Voinnet O ( 2006 ) The diversity of RNA silencing pathways in plants. Trends Genet 22 : 268 – 280 Crossref Search ADS PubMed Buhtz A, Springer F, Chappell L, Baulcombe DC, Kehr J ( 2008 ) Identification and characterization of small RNAs from the phloem of Brassica napus. Plant J 53 : 739 – 749 Crossref Search ADS PubMed Chiou TJ ( 2007 ) The role of microRNAs in sensing nutrient stress. Plant Cell Environ 30 : 323 – 332 Crossref Search ADS PubMed Chiou TJ, Aung K, Lin SI, Wu CC, Chiang SF, Su Cl ( 2006 ) Regulation of phosphate homeostasis by microRNA in Arabidopsis. Plant Cell 18 : 412 – 421 Crossref Search ADS PubMed Combier JP, Frugier F, de Billy F, Boualem A, El-Yahyaoui F, Moreau S, Vernié T, Ott T, Gamas P, Crespi M, et al ( 2006 ) MtHAP2-1 is a key transcriptional regulator of symbiotic nodule development regulated by microRNA169 in Medicago truncatula. Genes Dev 20 : 3084 – 3088 Crossref Search ADS PubMed Dare AP, Schaffer RJ, Lin-Wang K, Allan AC, Hellens RP ( 2008 ) Identification of a cis-regulatory element by transient analysis of co-ordinately regulated genes. Plant Methods 4 : 17 Crossref Search ADS PubMed Duan K, Ding X, Zhang Q, Zhu H, Pan A, Huang J ( 2008 ) AtCopeg1, the unique gene originated from AtCopia95 retrotransposon family, is sensitive to external hormones and abiotic stresses. Plant Cell Rep 27 : 1065 – 1073 Crossref Search ADS PubMed Ebbs ML, Bender J ( 2006 ) Locus-specific control of DNA methylation by the Arabidopsis SUVH5 histone methyltransferase. Plant Cell 18 : 1166 – 1176 Crossref Search ADS PubMed Fahlgren N, Sullivan CM, Kasschau KD, Chapman EJ, Cumbie JS, Montgomery TA, Gilbert SD, Dasenko M, Backman TWH, Givan SA, et al ( 2009 ) Computational and analytical framework for small RNA profiling by high-throughput sequencing. RNA 15 : 992 – 1002 Crossref Search ADS PubMed Fujii H, Chiou TJ, Lin SI, Aung K, Zhu JK ( 2005 ) A miRNA involved in phosphate-starvation response in Arabidopsis. Curr Biol 15 : 2038 – 2043 Crossref Search ADS PubMed Gandikota M, Birkenbihl RP, Hohmann S, Cardon GH, Saedler H, Huijser P ( 2007 ) The miRNA156/157 recognition element in the 3′ UTR of the Arabidopsis SBP box gene SPL3 prevents early flowering by translational inhibition in seedlings. Plant J 49 : 683 – 693 Crossref Search ADS PubMed Haiser HJ, Karginov FV, Hannon GJ, Elliot MA ( 2008 ) Developmentally regulated cleavage of tRNAs in the bacterium Streptomyces coelicolor. Nucleic Acids Res 36 : 732 – 741 Crossref Search ADS PubMed Heisel SE, Zhang Y, Allen E, Guo L, Reynolds TL, Yang X, Kovalic D, Roberts JK ( 2008 ) Characterization of unique small RNA populations from rice grain. PLoS One 3 : e2871 Crossref Search ADS PubMed Hofacker IL ( 2003 ) Vienna RNA secondary structure server. Nucleic Acids Res 31 : 3429 – 3431 Crossref Search ADS PubMed Jochl C, Rederstorff M, Hertel J, Stadler PF, Hofacker IL, Schrettl M, Haas H, Huttenhofer A ( 2008 ) Small ncRNA transcriptome analysis from Aspergillus fumigatus suggests a novel mechanism for regulation of protein synthesis. Nucleic Acids Res 36 : 2677 – 2689 Crossref Search ADS PubMed Jones-Rhoades MW, Bartel DP ( 2004 ) Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell 14 : 787 – 799 Crossref Search ADS PubMed Jones-Rhoades MW, Bartel DP, Bartel B ( 2006 ) MicroRNAs and their regulatory roles in plants. Annu Rev Plant Biol 57 : 19 – 53 Crossref Search ADS PubMed Kawashima CG, Yoshimoto N, Maruyama-Nakashita A, Tsuchiya YN, Saito K, Takahashi H, Dalmay T ( 2009 ) Sulphur starvation induces the expression of microRNA-395 and one of its target genes but in different cell types. Plant J 57 : 313 – 321 Crossref Search ADS PubMed Lange H, Shropshire W Jr, Mohr H ( 1971 ) An analysis of phytochrome-mediated anthocyanin synthesis. Plant Physiol 47 : 649 – 655 Crossref Search ADS PubMed Lee SR, Collins K ( 2005 ) Starvation-induced cleavage of the tRNA anticodon loop in Tetrahymena thermophila. J Biol Chem 280 : 42744 – 42749 Crossref Search ADS PubMed Lenburg ME, O'Shea EK ( 1996 ) Signaling phosphate starvation. Trends Biochem Sci 21 : 383 – 387 Crossref Search ADS PubMed Li WX, Oono Y, Zhu J, He XJ, Wu JM, Iida K, Lu XY, Cui X, Jin H, Zhu JK ( 2008 a) The Arabidopsis NFYA5 transcription factor is regulated transcriptionally and posttranscriptionally to promote drought resistance. Plant Cell 20 : 2238 – 2251 Crossref Search ADS PubMed Li Y, Luo J, Zhou H, Liao JY, Ma LM, Chen YQ, Qu LH ( 2008 b) Stress-induced tRNA-derived RNAs: a novel class of small RNAs in the primitive eukaryote Giardia lamblia. Nucleic Acids Res 36 : 6048 – 6055 Crossref Search ADS PubMed Lin SI, Chiang SF, Lin WY, Chen JW, Tseng CY, Wu PC, Chiou TJ ( 2008 ) Regulatory network of microRNA399 and PHO2 by systemic signaling. Plant Physiol 147 : 732 – 746 Crossref Search ADS PubMed Lin WY, Lin SI, Chiou TJ ( 2009 ) Molecular regulators of phosphate homeostasis in plants. J Exp Bot 60 : 1427 – 1438 Crossref Search ADS PubMed Llave C, Xie Z, Kasschau KD, Carrington JC ( 2002 ) Cleavage of Scarecrow-like mRNA targets directed by a class of Arabidopsis miRNA. Science 297 : 2053 – 2056 Crossref Search ADS PubMed Lu C, Tej SS, Luo S, Haudenschild CD, Meyers BC, Green PJ ( 2005 ) Elucidation of the small RNA component of the transcriptome. Science 309 : 1567 – 1569 Crossref Search ADS PubMed Lu S, Sun YH, Chiang VL ( 2008 ) Stress-responsive microRNAs in Populus. Plant J 55 : 131 – 151 Crossref Search ADS PubMed Mallory AC, Vaucheret H ( 2006 ) Functions of microRNAs and related small RNAs in plants. Nat Genet 38 : S31 –S36 Crossref Search ADS PubMed Marschner H ( 1995 ). Mineral Nutrition of Higher Plants. Academic Press, London, pp 265–277, 537–594 Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, et al ( 2008 ) Criteria for annotation of plant microRNAs. Plant Cell 20 : 3186 – 3190 Crossref Search ADS PubMed Misson J, Raghothama KG, Jain A, Jouhet J, Block MA, Bligny R, Ortet P, Creff A, Somerville S, Rolland N, et al ( 2005 ) A genome-wide transcriptional analysis using Arabidopsis thaliana Affymetrix gene chips determined plant responses to phosphate deprivation. Proc Natl Acad Sci USA 102 : 11934 – 11939 Crossref Search ADS PubMed Navarro L, Dunoyer P, Jay F, Arnold B, Dharmasiri N, Estelle M, Voinnet O, Jones JDG ( 2006 ) A plant miRNA contributes to antibacterial resistance by repressing auxin signaling. Science 312 : 436 – 439 Crossref Search ADS PubMed Nobuta K, Venu RC, Lu C, Beló A, Vemaraju K, Kulkarni K, Wang W, Pillay M, Green PJ, Wang GL, et al ( 2007 ) An expression atlas of rice mRNAs and small RNAs. Nat Biotechnol 25 : 473 – 477 Crossref Search ADS PubMed Okamura K, Phillips MD, Tyler DM, Duan H, Chou YT, Lai EC ( 2008 ) The regulatory activity of microRNA* species has substantial influence on microRNA and 3′ UTR evolution. Nat Struct Mol Biol 15 : 354 – 363 Crossref Search ADS PubMed Pall GS, Codony-Servat C, Byrne J, Ritchie L, Hamilton A ( 2007 ) Carbodiimide-mediated cross-linking of RNA to nylon membranes improves the detection of siRNA, miRNA and piRNA by northern blot. Nucleic Acids Res 35 : e60 Crossref Search ADS PubMed Pant BD, Buhtz A, Kehr J, Scheible WR ( 2008 ) MicroRNA399 is a long-distance signal for the regulation of plant phosphate homeostasis. Plant J 53 : 731 – 738 Crossref Search ADS PubMed Pant BD, Musialak-Lange M, Nuc P, May P, Buhtz A, Kehr J, Walther D, Scheible WR ( 2009 ) Identification of nutrient-responsive Arabidopsis and rapeseed microRNAs by comprehensive real-time polymerase chain reaction profiling and small RNA sequencing. Plant Physiol 150 : 1541 – 1555 Crossref Search ADS PubMed Peng M, Hannam C, Gu H, Bi YM, Rothstein SJ ( 2007 ) A mutation in NLA, which encodes a RING-type ubiquitin ligase, disrupts the adaptability of Arabidopsis to nitrogen limitation. Plant J 50 : 320 – 337 Crossref Search ADS PubMed Peragine A, Yoshikawa M, Wu G, Albrecht HL, Poethig RS ( 2004 ) SGS3 and SGS2/SDE1/RDR6 are required for juvenile development and the production of trans-acting siRNAs in Arabidopsis. Genes Dev 18 : 2368 – 2379 Crossref Search ADS PubMed Poirier Y, Bucher M ( 2002 ) Phosphate transport and homeostasis in Arabidopsis. In CR Somerville, EM Meyerowitz, eds, The Arabidopsis Book. American Society of Plant Biologists, Rockville, MD, doi: 10.1199/tab.0024 , http://www.aspb.org/publications/arabidopsis/ Crossref Raghothama KG ( 1999 ) Phosphate acquisition. Annu Rev Plant Physiol Plant Mol Biol 50 : 665 – 693 Crossref Search ADS PubMed Rajagopalan R, Vaucheret H, Trejo J, Bartel DP ( 2006 ) A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev 20 : 3407 – 3425 Crossref Search ADS PubMed Rubio V, Linhares F, Solano R, Martin AC, Iglesias J, Leyva A, Paz-Ares J ( 2001 ) A conserved MYB transcription factor involved in phosphate starvation signaling both in vascular plants and in unicellular algae. Genes Dev 15 : 2122 – 2133 Crossref Search ADS PubMed Schachtman DP, Shin R ( 2007 ) Nutrient sensing and signaling: NPKS. Annu Rev Plant Biol 58 : 47 – 69 Crossref Search ADS PubMed Schwarz S, Grande A, Bujdoso N, Saedler H, Huijser P ( 2008 ) The microRNA regulated SBP-box genes SPL9 and SPL15 control shoot maturation in Arabidopsis. Plant Mol Biol 67 : 183 – 195 Crossref Search ADS PubMed Shin R, Berg RH, Schachtman DP ( 2005 ) Reactive oxygen species and root hairs in Arabidopsis root response to nitrogen, phosphorus and potassium deficiency. Plant Cell Physiol 46 : 1350 – 1357 Crossref Search ADS PubMed Sunkar R, Kapoor A, Zhu JK ( 2006 ) Posttranscriptional induction of two Cu/Zn superoxide dismutase genes in Arabidopsis is mediated by downregulation of miR398 and important for oxidative stress tolerance. Plant Cell 18 : 2051 – 2065 Crossref Search ADS PubMed Sunkar R, Zhu JK ( 2004 ) Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell 16 : 2001 – 2019 Crossref Search ADS PubMed Thompson DM, Lu C, Green PJ, Parker R ( 2008 ) tRNA cleavage is a conserved response to oxidative stress in eukaryotes. RNA 14 : 2095 – 2103 Crossref Search ADS PubMed Thompson DM, Parker R ( 2009 ) Stressing out over tRNA cleavage. Cell 138 : 215 – 219 Crossref Search ADS PubMed Vaucheret H ( 2006 ) Post-transcriptional small RNA pathways in plants: mechanisms and regulations. Genes Dev 20 : 759 – 771 Crossref Search ADS PubMed Vazquez F ( 2006 ) Arabidopsis endogenous small RNAs: highways and byways. Trends Plant Sci 11 : 460 – 468 Crossref Search ADS PubMed Vazquez F, Vaucheret H, Rajagopalan R, Lepers C, Gasciolli V, Mallory AC, Hilbert JL, Bartel DP, Crete P ( 2004 ) Endogenous trans-acting siRNAs regulate the accumulation of Arabidopsis mRNAs. Mol Cell 16 : 69 – 79 Crossref Search ADS PubMed Voinnet O ( 2009 ) Origin, biogenesis, and activity of plant microRNAs. Cell 136 : 669 – 687 Crossref Search ADS PubMed Wang JW, Schwab R, Czech B, Mica E, Weigel D ( 2008 ) Dual effects of miR156-targeted SPL genes and CYP78A5/KLUH on plastochron length and organ size in Arabidopsis thaliana. Plant Cell 20 : 1231 – 1243 Crossref Search ADS PubMed Wu G, Poethig RS ( 2006 ) Temporal regulation of shoot development in Arabidopsis thaliana by miR156 and its target SPL3. Development 133 : 3539 – 3547 Crossref Search ADS PubMed Yaeno T, Iba K ( 2008 ) BAH1/NLA, a RING-type ubiquitin E3 ligase, regulates the accumulation of salicylic acid and immune responses to Pseudomonas syringae DC3000. Plant Physiol 148 : 1032 – 1041 Crossref Search ADS PubMed Yamasaki H, Abdel-Ghany SE, Cohu CM, Kobayashi Y, Shikanai T, Pilon M ( 2007 ) Regulation of copper homeostasis by micro-RNA in Arabidopsis. J Biol Chem 282 : 16369 – 16378 Crossref Search ADS PubMed Zhang S, Sun L, Kragler F ( 2009 ) The phloem delivered RNA pool contains small non-coding RNAs and interferes with translation. Plant Physiol 150 : 378 – 387 Crossref Search ADS PubMed Zhu QH, Spriggs A, Matthew L, Fan L, Kennedy G, Gubler F, Helliwell C ( 2008 ) A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains. Genome Res 18 : 1456 – 1465 Crossref Search ADS PubMed Author notes 1 This work was supported by Academia Sinica, Taiwan (grant nos. AS97FPL20–1 and AS98CDAL11), and the National Science Council, Taiwan (grant no. 97–2321–B–001–018). 2 These authors contributed equally to the article. * Corresponding author; e-mail tjchiou@gate.sinica.edu.tw. The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: Tzyy-Jen Chiou (tjchiou@gate.sinica.edu.tw). [W] The online version of this article contains Web-only data. [OA] Open Access articles can be viewed online without a subscription. www.plantphysiol.org/cgi/doi/10.1104/pp.109.147280 © 2009 American Society of Plant Biologists © The Author(s) 2009. Published by Oxford University Press on behalf of American Society of Plant Biologists. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Plant Physiology Oxford University Press

Uncovering Small RNA-Mediated Responses to Phosphate Deficiency in Arabidopsis by Deep Sequencing

Loading next page...
 
/lp/oxford-university-press/uncovering-small-rna-mediated-responses-to-phosphate-deficiency-in-dtPtieIvl8

References (68)

Publisher
Oxford University Press
Copyright
Copyright © 2021 American Society of Plant Biologists
ISSN
0032-0889
eISSN
1532-2548
DOI
10.1104/pp.109.147280
Publisher site
See Article on Publisher Site

Abstract

Abstract Recent studies have demonstrated the important role of plant microRNAs (miRNAs) under nutrient deficiencies. In this study, deep sequencing of Arabidopsis (Arabidopsis thaliana) small RNAs was conducted to reveal miRNAs and other small RNAs that were differentially expressed in response to phosphate (Pi) deficiency. About 3.5 million sequence reads corresponding to 0.6 to 1.2 million unique sequence tags from each Pi-sufficient or Pi-deficient root or shoot sample were mapped to the Arabidopsis genome. We showed that upon Pi deprivation, the expression of miR156, miR399, miR778, miR827, and miR2111 was induced, whereas the expression of miR169, miR395, and miR398 was repressed. We found cross talk coordinated by these miRNAs under different nutrient deficiencies. In addition to miRNAs, we identified one Pi starvation-induced DICER-LIKE1-dependent small RNA derived from the long terminal repeat of a retrotransposon and a group of 19-nucleotide small RNAs corresponding to the 5′ end of tRNA and expressed at a high level in Pi-starved roots. Importantly, we observed an increased abundance of TAS4-derived trans-acting small interfering RNAs (ta-siRNAs) in Pi-deficient shoots and uncovered an autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(−) that regulates the biosynthesis of anthocyanin. This finding sheds light on the regulatory network between miRNA/ta-siRNA and its target gene. Of note, a substantial amount of miR399* accumulated under Pi deficiency. Like miR399, miR399* can move across the graft junction, implying a potential biological role for miR399*. This study represents a comprehensive expression profiling of Pi-responsive small RNAs and advances our understanding of the regulation of Pi homeostasis mediated by small RNAs. Plants acquire mineral nutrients that are channeled into diverse biochemical and metabolic pathways essential for the completion of their life cycle. Phosphorus (P), one of the essential macronutrients of plants, is often not readily accessible to plants because of the low availability of phosphate (Pi) in the soil, the major acquired form of P (Marschner, 1995). To overcome this limitation, plants have evolved a number of responses to conserve and mobilize internal Pi and to enhance the acquisition of external Pi (Raghothama, 1999; Poirier and Bucher, 2002; Lin et al., 2009). Recently, the molecular basis underlying these adaptive responses has been an intense subject of research. It is now recognized that small regulatory RNAs contribute to the complexity of gene regulation. MicroRNAs (miRNAs) and small interfering RNAs (siRNAs) of 21 to 24 nucleotides long are two major classes of small regulatory RNAs in plants, functioning as negative regulators of gene expression (Brodersen and Voinnet, 2006; Jones-Rhoades et al., 2006; Mallory and Vaucheret, 2006). miRNAs are processed from single-stranded RNA precursors capable of forming imperfectly complementary hairpin structures by the RNase III enzyme DICER-LIKE1 (DCL1) or DCL4. In most cases, plant miRNAs down-regulate gene expression by direct cleavage of their target transcripts, although translational repression has been proposed to be common as well (Brodersen et al., 2008). In contrast, siRNAs are generated from perfectly complementary long double-stranded RNA precursors, requiring the activity of RNA-dependent RNA polymerase. Endogenous siRNAs can be further classified into trans-acting siRNAs (ta-siRNAs), natural antisense transcript-derived siRNAs, and heterochromatin siRNAs according to their genomic origin and the interaction of distinct components (Vaucheret, 2006; Vazquez, 2006). ta-siRNAs are phased 21-nucleotide RNA molecules whose synthesis is triggered by the cleavage of miRNA on the TAS transcript (Peragine et al., 2004; Vazquez et al., 2004; Allen et al., 2005). They act in-trans to down-regulate the expression of unrelated loci from which they are produced. The biological functions of miRNA cover a broad range of physiological processes. In addition to their involvement in normal plant growth and development, the expression of many miRNAs is regulated by various biotic and abiotic stresses, which may contribute to the development of adaptive responses to overcome unfavorable conditions (Sunkar and Zhu, 2004; Navarro et al., 2006; Chiou, 2007; Li et al., 2008a; Lu et al., 2008). Recent findings have also demonstrated the important role of several miRNAs under nutrient deficiencies. miR395, miR398, and miR399 are up-regulated under sulfur (S), copper (Cu), and Pi deficiency, respectively. miR395 is involved in S assimilation and allocation by adjusting the activity of ATP sulfurylase (APS) and restricting the spatial expression of a sulfate transporter (AtSULTR2;1) responsible for root-to-shoot transport of sulfate (Jones-Rhoades and Bartel, 2004; Kawashima et al., 2009). Other than a role in alleviating oxidative stresses (Sunkar et al., 2006), miR398 modulates the expression of Cu/zinc (Zn) superoxide dismutase according to the availability of Cu (Yamasaki et al., 2007). Moreover, cell-specific expression of miR169 is crucial in the development of nitrogen (N)-fixing nodules in Medicago truncatula (Combier et al., 2006). miR399 regulates Pi homeostasis by controlling the expression of PHO2 encoding a ubiquitin-conjugating E2 enzyme, UBC24 (Fujii et al., 2005; Aung et al., 2006; Bari et al., 2006; Chiou et al., 2006). It acts as a positive regulator to activate the Pi uptake and translocation when Pi supply is insufficient, while PHO2 functions as a negative regulator to suppress these activities in order to prevent Pi excess when external Pi is ample. Notably, recent results from reciprocal grafting suggest that miR399 may serve as a long-distance signal from shoots to suppress PHO2 expression in roots as a way of communication between the Pi status of shoots and the Pi transport activity of roots (Lin et al., 2008; Pant et al., 2008). The long-distance movement of miR399 could be an early response to Pi deficiency. The shoot-derived miR399 is responsible for the degradation of PHO2 in the roots, where miR399 is not readily expressed (Lin et al., 2008). New technologies of high-throughput sequencing have become powerful tools to disclose the large inventory of small RNA species in plants. These deep sequencing strategies may discover new miRNAs or novel small RNA classes and provide quantitative profiling of small RNA expression (Lu et al., 2005; Rajagopalan et al., 2006; Nobuta et al., 2007; Heisel et al., 2008; Zhu et al., 2008). In this study, Solexa sequencing was employed to profile the changes of small RNA expression in response to Pi deficiency. Very recently, an independent study was reported in which P-responsive miRNAs were identified using quantitative reverse transcription (RT)-PCR of primary transcript of miRNAs followed by small RNA sequencing (Pant et al., 2009). In our study, we not only profiled the expression of small RNAs but also looked into their regulatory mechanisms and potential functions. Importantly, we found two unique classes of small RNAs up-regulated by Pi deficiency and uncovered an autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(−) that regulates the biosynthesis of anthocyanin. Changes in the abundance of specific groups of Pi-responsive small RNAs may mediate the activity or expression of target genes that are important in the development or regulation of adaptive responses, enabling plants to survive in a low-Pi environment. RESULTS AND DISCUSSION Overview of Small RNA Profiles in Response to Pi Deficiency Four small RNA libraries from Pi-sufficient roots (R+Pi), Pi-deficient roots (R−Pi), Pi-sufficient shoots (S+Pi), and Pi-deficient shoots (S−Pi) were constructed by Solexa high-throughput sequencing technology. After adaptor trimming, the numbers of sequence reads and unique sequence signatures from the raw data were calculated and then mapped to the Arabidopsis (Arabidopsis thaliana) genome (Supplemental Table S1). In summary, the sequencing yielded 4.8 to 5.1 million raw reads from each library, and approximately 90% of the reads remained after trimming the adaptor sequences. Approximately 3.5 million sequence reads (approximately 77%–83% of trimmed sequence reads) corresponding to 0.6 to 1.2 million unique sequence signatures could be mapped perfectly onto the genome (The Arabidopsis Information Resource 8). Only those perfectly mapped sequences were analyzed further. The sequence reads that mapped to single and multiple genomic locations were defined, respectively, as single-hit reads and multiple-hit reads. One-third of mapped sequences were single-hit reads in the root and one-half were single-hit reads in the shoot. The length of the mapped sequences ranged from 16 to 27 nucleotides. We first examined the correlation between the length of small RNAs and the proportion of total sequence reads or unique sequence signatures (Fig. 1 Figure 1. Open in new tabDownload slide Global profiling of small RNAs. Correlations between the small RNA length and the proportion of total sequence reads (A) or unique sequence signatures (B) are shown. nt, Nucleotides. Figure 1. Open in new tabDownload slide Global profiling of small RNAs. Correlations between the small RNA length and the proportion of total sequence reads (A) or unique sequence signatures (B) are shown. nt, Nucleotides. ). Similar distributions were observed in +Pi and −Pi libraries. Consistent with the earlier reports, there were two major peaks at 21 and 24 nucleotides in the total sequence reads of shoot libraries (Fig. 1A). Unexpectedly, an additional peak at 19 nucleotides was observed in the root libraries. This peak has not been noticed in the previous deep sequencing databases, probably because the pure root sample has not been examined. Further analysis revealed that the majority of these 19-nucleotide sequences originated from tRNAs (see below). When the unique sequence signatures were examined, the patterns of all four libraries were nearly identical (Fig. 1B). The 24-nucleotide small RNAs were dominant in either their reads or unique sequences, indicating that they are rich in sequence diversity. Overrepresentation of tRNA-Derived Small RNAs in the Pi-Starved Roots Mapped sequences were classified according to the annotation of the genome (Supplemental Table S2). Most sequences were mapped onto the nonannotated intergenic regions (30.12%–42.30%) and to a lesser extent onto MIRNA genes and protein-coding genes. Strikingly, a significantly higher proportion of small RNAs was mapped to tRNA genes in the root libraries (24.29%–34.32%) than in the shoot libraries (3.33%–5.58%). The percentage in the −Pi root library (34.32%) was even higher than that in the +Pi root library (24.29%). When the origin of these small RNAs was further analyzed (Supplemental Table S3), we did not find any correlation between the abundance of small RNAs from specific tRNA species and their codon usage. We were surprised to find that a 19-nucleotide sequence cleaved from the 5′ end of Gly-tRNATCC represented over 80% of tRNA-derived small RNAs in the roots and accounted for up to 18.44% and 27.70% of total sequence reads in the +Pi and −Pi root libraries, respectively, compared with only 1.00% to 1.79% in the shoot libraries (highlighted in Supplemental Table S3). For validation, we employed RNA gel analyses using the probes from three different segments corresponding to the 5′ and 3′ ends and the central anticodon loop of tRNA (Fig. 2B Figure 2. Open in new tabDownload slide Overrepresentation of tRNA-derived small RNAs in the Pi-starved roots. A, Small RNA gel analyses of tRNA-derived small RNAs in +Pi or −Pi roots or shoots. B, The tRNA regions corresponding to the 5′ and 3′ end probes are designated by the black line, and the anticodon probe is designated by the gray line. The black arrowheads in A and B indicate the 5′ cleaved products and corresponding cleavage site in tRNA, respectively. The vertical gray lines in A indicate the products of tRNA halves from the cleavage at the anticodon loop (gray arrowhead in B). nt, Nucleotides. Figure 2. Open in new tabDownload slide Overrepresentation of tRNA-derived small RNAs in the Pi-starved roots. A, Small RNA gel analyses of tRNA-derived small RNAs in +Pi or −Pi roots or shoots. B, The tRNA regions corresponding to the 5′ and 3′ end probes are designated by the black line, and the anticodon probe is designated by the gray line. The black arrowheads in A and B indicate the 5′ cleaved products and corresponding cleavage site in tRNA, respectively. The vertical gray lines in A indicate the products of tRNA halves from the cleavage at the anticodon loop (gray arrowhead in B). nt, Nucleotides. ). The small RNAs derived from Asp-tRNAGTC were also examined using the 5′ end probe. Consistent with the sequencing data, the 19-nucleotide RNAs from Gly-tRNATCC and Asp-tRNAGTC corresponding to the 5′ end of tRNA were highly accumulated in the roots but much less in the shoots (Fig. 2A). The −Pi roots showed the maximum accumulation. When probing with the 5′ or 3′ end probe, additional bands at 30 to 40 nucleotides likely representing tRNA halves cleaved at the anticodon loop were observed preferentially in the roots (Fig. 2A). These results revealed a spatial and temporal expression pattern of small RNAs derived from the specific cleavage on tRNA molecules rather than random degradation and may represent a novel processing of small RNAs. Studies in microorganisms or mammalian cells have suggested that the tRNA-derived small RNAs are accumulated in a developmentally regulated manner and become dominant in specific tissues or under specific stress conditions (Lee and Collins, 2005; Haiser et al., 2008; Jochl et al., 2008; Li et al., 2008b; Thompson et al., 2008; Pant et al., 2009; Thompson and Parker, 2009). These small RNAs could be associated with the quality control of protein synthesis or regulation of gene expression. Moreover, a recent report showed that tRNA half molecules present in the phloem sap of pumpkin (Cucurbita maxima) were able to inhibit translational activity in vitro (Zhang et al., 2009). The authors proposed that these phloem-delivered tRNA halves may be a long-distance signal to coordinate the metabolic status between source and sink tissues. The differential accumulation of tRNA-derived small RNAs between roots and shoots observed here may represent the consequence of such long-distance movement. Expression Profiling of miRNAs in Response to Pi Deprivation Forty-five and 55 differentially expressed MIRNA genes with greater than 1.5-fold relative change in sequence counts were identified in roots and shoots, respectively (Fig. 3 Figure 3. Open in new tabDownload slide Differentially expressed MIRNA genes in response to Pi deficiency in roots (A) and shoots (B). The significantly differentially expressed MIRNA genes with greater than 1.5-fold relative change are shown. The miRNAs with a relative change ratio greater than 3 are highlighted in boldface. Figure 3. Open in new tabDownload slide Differentially expressed MIRNA genes in response to Pi deficiency in roots (A) and shoots (B). The significantly differentially expressed MIRNA genes with greater than 1.5-fold relative change are shown. The miRNAs with a relative change ratio greater than 3 are highlighted in boldface. ; Supplemental Table S4). MIR399s, MIR778, MIR827, and MIR2111 were clearly up-regulated, whereas MIR169s, MIR395s, and MIR398s were down-regulated, in −Pi roots and/or shoots. The up-regulation of various species from the MIR399 family by −Pi was consistent with our previous observations (Lin et al., 2008), suggesting the reliability of our sequencing data. Among all MIRNA genes, the MIR399 family showed the highest degree of induction in both tissues. The expression of miRNAs with relative changes greater than a ratio of 3 (highlighted in boldface in Fig. 3) was validated by RNA gel analysis. To examine whether the response of these miRNAs was specific to −Pi or not, the RNA isolated from different nutrient deficiencies, including N, Pi, potassium (K), S, Cu, and iron (Fe), was inspected at once (Fig. 4 Figure 4. Open in new tabDownload slide Validation of differentially expressed miRNAs under different nutrient deficiencies. Small RNA gel analyses of Pi starvation-up-regulated miRNAs (A) and Pi starvation-down-regulated miRNAs (B) are shown. miR778 was detected using the probe corresponding to miR778-3p.1 and miR778-3p.2. Figure 4. Open in new tabDownload slide Validation of differentially expressed miRNAs under different nutrient deficiencies. Small RNA gel analyses of Pi starvation-up-regulated miRNAs (A) and Pi starvation-down-regulated miRNAs (B) are shown. miR778 was detected using the probe corresponding to miR778-3p.1 and miR778-3p.2. ). miR156 was highly expressed but unaffected by nutrient deficiencies in the shoots; however, it was up-regulated by −Pi, −N, or −K, with the highest induction under −N. miR156 regulates a group of SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE (SPL) transcription factors involved in juvenile-to-adult vegetative phase transition and other development processes (Wu and Poethig, 2006; Gandikota et al., 2007; Schwarz et al., 2008; Wang et al., 2008). The physiological role for the up-regulation of miR156 in the −Pi, −N, or −K root remains to be explored. On the other hand, miR399, miR778, and miR827 were specifically up-regulated by −Pi (Fig. 4A). It is interesting that several distinct small RNA sequences were generated by cleavage at different positions of the miR778 precursor (Supplemental Fig. S1). Because the corresponding miRNA* sequence could be detected for each cleaved RNA duplex, we suggest to rename miR778 as miR778-5p and the other two forms derived from the miR778 precursor as miR778-3p.1 and miR778-3p.2, respectively. Multiple processing of the miR778 precursor has not been reported previously, probably due to the low recovered sequence reads. All miR778 species were up-regulated by −Pi, with the highest sequence count of miR778-3p.1 (Supplemental Table S4). The target genes of miR778-5p encode SET domain-containing proteins, which are involved in the regulation of histone methylation (Ebbs and Bender, 2006). Although the target genes of miR778-3p.1 and miR778-3p.2 await identification, generation of multiple miRNAs from the same precursor has increased the complexity for the miR778 regulation networks. The target genes of miR827 encode SPX domain-containing proteins. In yeast cells, SPX domain-containing proteins participate in Pi transport or sensing (Lenburg and O'Shea, 1996). One of the target genes, At1g02860, encodes a nucleus-localized ubiquitin E3 ligase containing the RING and SPX domains, which was recently characterized as NITROGEN LIMITATION ADAPTATION (NLA) involved in the N recycle during N limitation (Peng et al., 2007) or as BENZOIC ACID HYPERSENSTIVE1 (BAH1) involved in immune responses (Yaeno and Iba, 2008). Here, we validated the targeting of miR827 to At1g02860 mRNA by RNA ligase-mediated (RLM) 5′-RACE analysis (Supplemental Fig. S2A) and observed the down-regulation of At1g02860 under −Pi (Supplemental Fig. S2B). At1g02860 may be regulated by many factors and have multiple functions in response to diverse external stimuli. The miR827-mediated regulation of At1g02860 under −Pi conditions may imply an interaction between N and P, although miR827 only responds to −Pi but not −N (Fig. 4A). The expression of several down-regulated miRNAs was also confirmed (Fig. 4B). miR169 was down-regulated in −Pi, −N, or −S roots but to a lesser extent in the shoots. Down-regulation of miR169 was much more severe under −N. miR395 was down-regulated in −Pi or −N roots and shoots. miR398 was down-regulated by −Pi, −N, −K, or −Fe in both roots and shoots. In contrast, miR395 and miR398 were up-regulated by −S and −Cu, respectively, which is in agreement with previous findings (Jones-Rhoades and Bartel, 2004; Yamasaki et al., 2007). These observations suggest that miRNAs can coordinate cross talk among different nutrient-deficient responses. Noticing the up-regulation of APS4 and SULTR2;1 under −Pi (Supplemental Fig. S2C), we speculate that the suppression of miR395 under −Pi may somewhat contribute to this effect. Up-regulation of APS4 and SULTR2;1 may result in the increase of sulfate assimilation and translocation for sulfolipid biosynthesis, a compensation for the reduced phospholipids during −Pi. miR169 and its target gene, HAP2/NFYA transcription factor, are critical in the development of nodules (Combier et al., 2006). Additionally, miR169a and NFYA5 have been shown to be involved in drought tolerance (Li et al., 2008a). NFYA5 functions upstream to activate several drought-responsive genes, including those involved in oxidative stress. The target genes of miR398 encoding Cu/Zn superoxide dismutase are also involved in the response to oxidative stress (Sunkar et al., 2006). Several nutrient deficiencies, such as N, Pi, K, and S, can trigger oxidative stresses (Shin et al., 2005; Schachtman and Shin, 2007). Changes of miR169 and miR398 expression could be involved in adjusting a plant's ability to survive the oxidative stress caused by nutrient deficiencies. Interestingly, the opposite expression of miR398 in response to Cu and Fe deficiencies may reflect an elegant regulation of the preferential expression of different superoxide dismutase, Cu/Zn-superoxide dismutase or Fe-superoxide dismutase, coupling with different metal cofactors. These Pi-responsive miRNAs and the potential functions of corresponding target genes relative to adaptive responses to −Pi are summarized in Table I Table I. Summary of validated Pi-responsive miRNAs and ta-siRNAs and the potential function involved under Pi deficiency miRNA/ta-siRNA . . . Target Genes . . . Familya . Speciesb . Expressionc . Gene Productd . Expressione . Potential Rolef . miR169 a, b, c, d, e, f, g, h, i, j, k, l, m, n Down HAP2 transcription factors (At1g17590, At1g54160, At1g72830, At3g05690, At3g20910, At5g06510, At5g12840) Up Involved in antioxidative stress miR395 b, c, f Down ATP sulfurylase (APS1, APS3, APS4; At3g22890, At4g14680, At5g43780) Up Increase sulfate translocation and utilization for sulfolipid biosynthesis Sulfate transporter (SULTR2;1; At5g10180) Up miR398 a, b, c Down Cu/Zn superoxide dismutase (CSD; CSD1, At1g08830; CSD2, At2g28190) Unchanged Scavenge reactive oxygen species Cytochrome c oxidase (At3g15640) Unchanged miR402 Downg DEMETER-LIKE protein 3 (DML3; At4g34060) Undetermined miR156 a, b, c, d, e, f, g Up (root) SPL transcription factors (At1g53160, At2g33810, At3g15270, At5g43270, At1g27360, At1g27370, At1g69170, At2g42200, At3g57920, At5g50570, At5g50670) Down miR399 a, b, c, d, f Up Ubiquitin conjugase E2 (UBC24/PHO2; At2g33770) Down Pi uptake and translocation miR778 5p Up SET domain-containing protein (At2g22740, At2g35160) Down 3p.1 Up 3p.2 Up miR827b Up Ubiquitin E3 ligase with RING and SPX domains (NLA/BAH1; At1g02860) SPX domain-containing protein (At1g63010) Down Unchanged Nutrient recycle miR2111 a-5p, b-5p Up Kelch repeat-containing F-box protein (At3g27150) Up a-3p, b-3p Up miR828 Up (shoot) TAS4 Anthocyanin biosynthesis TAS4-siRNAb D1(+), D2(+), D3(+) TAS4-siR81(−) Up (shoot) Up (shoot) MYB75, MYB90, and MYB113 transcription factors (At1g56650, At1g66390, At1g66370) Up Anthocyanin biosynthesis miRNA/ta-siRNA . . . Target Genes . . . Familya . Speciesb . Expressionc . Gene Productd . Expressione . Potential Rolef . miR169 a, b, c, d, e, f, g, h, i, j, k, l, m, n Down HAP2 transcription factors (At1g17590, At1g54160, At1g72830, At3g05690, At3g20910, At5g06510, At5g12840) Up Involved in antioxidative stress miR395 b, c, f Down ATP sulfurylase (APS1, APS3, APS4; At3g22890, At4g14680, At5g43780) Up Increase sulfate translocation and utilization for sulfolipid biosynthesis Sulfate transporter (SULTR2;1; At5g10180) Up miR398 a, b, c Down Cu/Zn superoxide dismutase (CSD; CSD1, At1g08830; CSD2, At2g28190) Unchanged Scavenge reactive oxygen species Cytochrome c oxidase (At3g15640) Unchanged miR402 Downg DEMETER-LIKE protein 3 (DML3; At4g34060) Undetermined miR156 a, b, c, d, e, f, g Up (root) SPL transcription factors (At1g53160, At2g33810, At3g15270, At5g43270, At1g27360, At1g27370, At1g69170, At2g42200, At3g57920, At5g50570, At5g50670) Down miR399 a, b, c, d, f Up Ubiquitin conjugase E2 (UBC24/PHO2; At2g33770) Down Pi uptake and translocation miR778 5p Up SET domain-containing protein (At2g22740, At2g35160) Down 3p.1 Up 3p.2 Up miR827b Up Ubiquitin E3 ligase with RING and SPX domains (NLA/BAH1; At1g02860) SPX domain-containing protein (At1g63010) Down Unchanged Nutrient recycle miR2111 a-5p, b-5p Up Kelch repeat-containing F-box protein (At3g27150) Up a-3p, b-3p Up miR828 Up (shoot) TAS4 Anthocyanin biosynthesis TAS4-siRNAb D1(+), D2(+), D3(+) TAS4-siR81(−) Up (shoot) Up (shoot) MYB75, MYB90, and MYB113 transcription factors (At1g56650, At1g66390, At1g66370) Up Anthocyanin biosynthesis a The Pi-responsive miRNAs reported by Pant et al. (2009) are underlined. b Only those species with statistically significant changes in expression are shown. Boldface letters indicate the existence of the P1BS cis-element (GNATATNC), a binding site for PHR1, in the 2-kb sequences upstream of miRNA hairpins. There are P1BS elements in the upstream sequences of MIR827 and TAS4 genes. c Change of small RNA expression due to Pi deficiency. d The target genes are based on the Arabidopsis Small RNA Project database (http://asrp.cgrb.oregonstate.edu/) except miR2111 (this study). The target genes with altered expression by Pi deficiency are indicated in boldface. e Change of target gene expression to Pi deficiency. The results were based on the microarry analysis from Genevestigator (https://www.genevestigator.ethz.ch/gv/index.jsp) except for the target genes of miR399 (Chiou et al., 2006), miR2111-5p and TAS4-siR81(−) (Figs. 5 and 7), miR395, and miR827 (Supplemental Fig. S2), whose expression was analyzed by quantitative RT-PCR. Changes larger than 2-fold are indicated. f Potential role of target gene under Pi deficiency. g Because of its low abundance, the expression of miR402 was validated by quantitative RT-PCR of primary transcript. Open in new tab Table I. Summary of validated Pi-responsive miRNAs and ta-siRNAs and the potential function involved under Pi deficiency miRNA/ta-siRNA . . . Target Genes . . . Familya . Speciesb . Expressionc . Gene Productd . Expressione . Potential Rolef . miR169 a, b, c, d, e, f, g, h, i, j, k, l, m, n Down HAP2 transcription factors (At1g17590, At1g54160, At1g72830, At3g05690, At3g20910, At5g06510, At5g12840) Up Involved in antioxidative stress miR395 b, c, f Down ATP sulfurylase (APS1, APS3, APS4; At3g22890, At4g14680, At5g43780) Up Increase sulfate translocation and utilization for sulfolipid biosynthesis Sulfate transporter (SULTR2;1; At5g10180) Up miR398 a, b, c Down Cu/Zn superoxide dismutase (CSD; CSD1, At1g08830; CSD2, At2g28190) Unchanged Scavenge reactive oxygen species Cytochrome c oxidase (At3g15640) Unchanged miR402 Downg DEMETER-LIKE protein 3 (DML3; At4g34060) Undetermined miR156 a, b, c, d, e, f, g Up (root) SPL transcription factors (At1g53160, At2g33810, At3g15270, At5g43270, At1g27360, At1g27370, At1g69170, At2g42200, At3g57920, At5g50570, At5g50670) Down miR399 a, b, c, d, f Up Ubiquitin conjugase E2 (UBC24/PHO2; At2g33770) Down Pi uptake and translocation miR778 5p Up SET domain-containing protein (At2g22740, At2g35160) Down 3p.1 Up 3p.2 Up miR827b Up Ubiquitin E3 ligase with RING and SPX domains (NLA/BAH1; At1g02860) SPX domain-containing protein (At1g63010) Down Unchanged Nutrient recycle miR2111 a-5p, b-5p Up Kelch repeat-containing F-box protein (At3g27150) Up a-3p, b-3p Up miR828 Up (shoot) TAS4 Anthocyanin biosynthesis TAS4-siRNAb D1(+), D2(+), D3(+) TAS4-siR81(−) Up (shoot) Up (shoot) MYB75, MYB90, and MYB113 transcription factors (At1g56650, At1g66390, At1g66370) Up Anthocyanin biosynthesis miRNA/ta-siRNA . . . Target Genes . . . Familya . Speciesb . Expressionc . Gene Productd . Expressione . Potential Rolef . miR169 a, b, c, d, e, f, g, h, i, j, k, l, m, n Down HAP2 transcription factors (At1g17590, At1g54160, At1g72830, At3g05690, At3g20910, At5g06510, At5g12840) Up Involved in antioxidative stress miR395 b, c, f Down ATP sulfurylase (APS1, APS3, APS4; At3g22890, At4g14680, At5g43780) Up Increase sulfate translocation and utilization for sulfolipid biosynthesis Sulfate transporter (SULTR2;1; At5g10180) Up miR398 a, b, c Down Cu/Zn superoxide dismutase (CSD; CSD1, At1g08830; CSD2, At2g28190) Unchanged Scavenge reactive oxygen species Cytochrome c oxidase (At3g15640) Unchanged miR402 Downg DEMETER-LIKE protein 3 (DML3; At4g34060) Undetermined miR156 a, b, c, d, e, f, g Up (root) SPL transcription factors (At1g53160, At2g33810, At3g15270, At5g43270, At1g27360, At1g27370, At1g69170, At2g42200, At3g57920, At5g50570, At5g50670) Down miR399 a, b, c, d, f Up Ubiquitin conjugase E2 (UBC24/PHO2; At2g33770) Down Pi uptake and translocation miR778 5p Up SET domain-containing protein (At2g22740, At2g35160) Down 3p.1 Up 3p.2 Up miR827b Up Ubiquitin E3 ligase with RING and SPX domains (NLA/BAH1; At1g02860) SPX domain-containing protein (At1g63010) Down Unchanged Nutrient recycle miR2111 a-5p, b-5p Up Kelch repeat-containing F-box protein (At3g27150) Up a-3p, b-3p Up miR828 Up (shoot) TAS4 Anthocyanin biosynthesis TAS4-siRNAb D1(+), D2(+), D3(+) TAS4-siR81(−) Up (shoot) Up (shoot) MYB75, MYB90, and MYB113 transcription factors (At1g56650, At1g66390, At1g66370) Up Anthocyanin biosynthesis a The Pi-responsive miRNAs reported by Pant et al. (2009) are underlined. b Only those species with statistically significant changes in expression are shown. Boldface letters indicate the existence of the P1BS cis-element (GNATATNC), a binding site for PHR1, in the 2-kb sequences upstream of miRNA hairpins. There are P1BS elements in the upstream sequences of MIR827 and TAS4 genes. c Change of small RNA expression due to Pi deficiency. d The target genes are based on the Arabidopsis Small RNA Project database (http://asrp.cgrb.oregonstate.edu/) except miR2111 (this study). The target genes with altered expression by Pi deficiency are indicated in boldface. e Change of target gene expression to Pi deficiency. The results were based on the microarry analysis from Genevestigator (https://www.genevestigator.ethz.ch/gv/index.jsp) except for the target genes of miR399 (Chiou et al., 2006), miR2111-5p and TAS4-siR81(−) (Figs. 5 and 7), miR395, and miR827 (Supplemental Fig. S2), whose expression was analyzed by quantitative RT-PCR. Changes larger than 2-fold are indicated. f Potential role of target gene under Pi deficiency. g Because of its low abundance, the expression of miR402 was validated by quantitative RT-PCR of primary transcript. Open in new tab . Up-regulation of miR399, miR778, and miR827 and down-regulation of miR398 by −Pi were also observed in a recent study in which a quantitative RT-PCR platform was established to profile the primary transcripts of miRNAs (Pant et al., 2009). The changes of these primary transcripts of miRNAs agreed with those of their mature miRNAs in response to −Pi, suggesting that they are mainly regulated at the transcriptional level. However, changes in miR156 and miR395 were not discovered in their system, probably because of differences in growth conditions or sample sources (Table I). To examine whether the regulation of these miRNAs depends on the PHOSPHATE STARVATION RESPONSE1 (PHR1) MYB transcription factor, which activates a group of Pi starvation-induced genes, including MIR399s (Rubio et al., 2001; Bari et al., 2006), we analyzed the existence of P1BS cis-element, the biding site of PHR1, within the 2-kb sequences upstream of miRNA hairpins. In addition to MIR399s, the P1BS element could be found in the upstream sequences of several species of MIR156, MIR169, MIR398c, and MIR827 genes (Table I), suggesting the possible regulation of these genes by PHR1. Characterization of the Pi-Responsive miRNA2111 Family We identified a miRNA gene family consisting of two genes (miR2111a and miR2111b) from the nonannotated sequence clusters, which were not reported as miRNAs during the time of our analysis. They were specifically up-regulated by −Pi but not by any of the other nutrient deficiencies tested (Fig. 5A Figure 5. Open in new tabDownload slide Characterization of miR2111 in response to Pi deficiency. A and B, Small RNA gel analyses of miR2111-5p and miR2111-3p under different nutrient-deficient conditions (A) and in small RNA biogenesis mutants (B). Col, Ecotype Columbia; Ler, ecotype Landsberg erecta; wt, wild type. C, Hairpin secondary structures of miR2111a and miR2111b precursors. D, RLM 5′-RACE analysis of miR2111-5p cleavage on At3g27150 mRNA. E, Quantitative RT-PCR analysis of At3g27150 mRNA in response to −Pi. One of two biological replicates is presented, and the error bars indicate the sd of two technical replicates. Figure 5. Open in new tabDownload slide Characterization of miR2111 in response to Pi deficiency. A and B, Small RNA gel analyses of miR2111-5p and miR2111-3p under different nutrient-deficient conditions (A) and in small RNA biogenesis mutants (B). Col, Ecotype Columbia; Ler, ecotype Landsberg erecta; wt, wild type. C, Hairpin secondary structures of miR2111a and miR2111b precursors. D, RLM 5′-RACE analysis of miR2111-5p cleavage on At3g27150 mRNA. E, Quantitative RT-PCR analysis of At3g27150 mRNA in response to −Pi. One of two biological replicates is presented, and the error bars indicate the sd of two technical replicates. ). They possessed a potential hairpin secondary structure and could produce a small RNA duplex with two-nucleotide overhang at the 3′ end (Fig. 5C), the primary criterion of miRNA definition (Meyers et al., 2008). Because a substantial amount of both complementary strands of the RNA duplex was detected, it was difficult to differentiate miRNA from miRNA*; thus, we named them miR2111a-5p, miR2111a-3p, miR2111b-5p, and miR2111b-3p. The sequences of miR2111a-5p and miR2111b-5p were identical, while those of miR2111a-3p and miR2111b-3p differed by four nucleotides. Although this family was also identified by other groups very recently (Fahlgren et al., 2009; Pant et al., 2009), its biogenesis and target genes have not been studied. The biogenesis of these miRNAs was analyzed using various mutants defective in the production of different classes of small RNAs (Fig. 5B). Here, we showed that the accumulation of these miRNAs was significantly reduced in the Pi-starved hen1 and dcl1 mutants. Nevertheless, some accumulation of miR2111 was still observed in the dcl1 mutant, suggesting that miR2111 can also be processed by another DCL, probably DCL4, in the absence of DCL1. The gene At3g27150, encoding a Kelch domain-containing F-box protein, was predicted to be a target of miR2111-5p. The cleavage of the mRNA of At3g27150 was confirmed by RLM 5′-RACE analysis (Fig. 5D). Interestingly, we did not observe the negative correlation between the expression of miR2111-5p and At3g27150. Instead, the steady-state level of At3g27150 mRNA was moderately induced by −Pi (Fig. 5E) regardless of the enhanced expression of miR2111-5p. Several reports have shown a positive temporal correlation between the expression of a miRNA and its target gene in spite of the negative regulatory role of miRNAs (Llave et al., 2002; Achard et al., 2004; Kawashima et al., 2009). The consequences of different miRNA regulatory circuits were classified into four categories, spatial restriction, mutual exclusion, dampening, and temporal regulation, in which the final expression level of the miRNA and its target gene may not be correlated (Voinnet, 2009). We hypothesize that miR2111-5p may regulate the spatial expression of At3g27150 by restricting its expression in specific cells and/or fine-tuning the expression of At3g27150 to maintain an optimal expression level. Alternatively, the possibility of translational repression cannot be ruled out. Nevertheless, regulation by miR2111 may be only one of many regulatory routes of At3g27150, because multiple regulatory pathways integrated from different levels are usually involved in gene regulation. It is important to note that the target genes of miR399, miR827, and miR2111-5p encode proteins involved in ubiquitin-mediated protein degradation, suggesting a crucial role of posttranslational regulation of protein abundance in the control of −Pi responses. The homologs of the miR2111 gene could be identified in many other dicotyledon plant species, such as Brassica rapa, M. truncatula, Lotus japonicus, and Vitis vinifera (Supplemental Fig. S3), but they could not be found in monocotyledon plants. This observation suggests that either the miR2111 gene had become lost in monocotyledons or, more likely, that it had emerged after the separation of the dicotyledons from the monocotyledons. The sequence of miR2111-5p is more conserved than that of miR2111-3p. The four miR2111 genes in M. truncatula are clustered in the same direction within a range of 40 kb. Interestingly, there are only 32 nucleotides between the middle two hairpins, which are likely derived from a polycistronic transcript. Conservation of the miR2111 species implies the importance of miR2111-mediated regulation for survival under the −Pi environment. Pi-Responsive Small RNAs Derived from LTR of Retrotransposon A cluster of small RNAs that originated from the long terminal repeat (LTR) of Copia95 retrotransposon was identified as showing a differential expression pattern in response to −Pi (Fig. 6A Figure 6. Open in new tabDownload slide Pi starvation-induced small RNAs derived from LTR of retrotransposon. A and B, Small RNA gel analysis of smRPi1LTR under different nutrient-deficient conditions (A) and in small RNA biogenesis mutants (B). Col, Ecotype Columbia; hetero, heterozygous plants; wt, wild type. C, Arrangement of Copia95 LTR in the intergenic region of At5g27990 and At5g28000 and the potential hairpin structure and location of smRPi1LTR. Figure 6. Open in new tabDownload slide Pi starvation-induced small RNAs derived from LTR of retrotransposon. A and B, Small RNA gel analysis of smRPi1LTR under different nutrient-deficient conditions (A) and in small RNA biogenesis mutants (B). Col, Ecotype Columbia; hetero, heterozygous plants; wt, wild type. C, Arrangement of Copia95 LTR in the intergenic region of At5g27990 and At5g28000 and the potential hairpin structure and location of smRPi1LTR. ). This small RNA cluster is located in the intergenic region of At5g27990 and At5g28000 containing a multiple array of truncated LTRs either as inverted or tandem repeats, suggesting the occurrence of several duplication events (Fig. 6C). The arrangement of inverted repeats allows the formation of a hairpin structure. The length of this hairpin is relatively long, at approximately 250 nucleotides. Unlike those miRNA-like small RNAs recently found in rice (Oryza sativa) grains (Heisel et al., 2008; Zhu et al., 2008), which are cleaved from a long hairpin, the small RNAs generated from this LTR hairpin are not phased. smRPi1LTR, the most abundant small RNA from the hairpin stem, accumulated specifically in −Pi roots but not in other nutrient-deficient tissues (Fig. 6A). DCL1 is essential for the biogenesis of smRPi1LTR, and its accumulation was not reduced in the dcl2,3,4 triple mutants or in the rdr1, rdr2, and rdr6 mutants (Fig. 6B). smRPi1LTR is likely generated from the cleavage of a single-stranded rather than a double-stranded RNA precursor, because the passenger strand of its duplex was identified from the hairpin stem with two mismatched base pairs. Because we were unable to detect smRPi1LTR in the Landsberg accession, we hypothesize that smRPi1LTR is a newly evolved small RNA due to rapid rearrangement of LTR and may represent an intermediate small RNA species from the transition of siRNAs to miRNAs. AtCopeg1 (for Copia evolved gene1; At2g04460), a protein-coding gene originating from a truncated AtCopia95 retrotransposon, was recently shown to be highly up-regulated by −Pi (Duan et al., 2008). It will be interesting to elucidate the relationship between AtCopeg1 and smRPi1LTR in terms of coevolution or cross-regulation. Up-Regulation of miR828- and TAS4-Derived siRNAs by Pi Deficiency and Autoregulation of PAP1/MYB75 The sequencing results revealed the up-regulation of ta-siRNAs derived from TAS4 and miR828, which initiates the cleavage of TAS4 RNA in −Pi shoots (Supplemental Table S4). These ta-siRNAs were up-regulated in the shoots by −Pi as well as −N (Fig. 7A Figure 7. Open in new tabDownload slide Autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(−). A, Small RNA gel analysis of TAS4-siRNAs under different nutrient-deficient conditions. B, Small RNA gel analysis of TAS4-siR81(−) and miR828 in the shoots of wild-type (wt) plants and tas4, mir828, and PAP1/MYB75 activation (pap1-D) mutants. C, Expression of PAP1/MYB75, PAP2/MYB90, and MYB113 by quantitative RT-PCR analysis in tas4 and mir828 mutants. One of two biological replicates is presented, and the error bars indicate the sd of two technical replicates. D, Alteration of anthocyanin accumulation in tas4, mir828, and pap1-D mutants. Error bars represent the sd (n = 5). FW, Fresh weight; OD, optical density. E, A proposed model of the autoregulatory pathway of PAP1/MYB75 via miR828 and TAS4-siR81(−). Figure 7. Open in new tabDownload slide Autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(−). A, Small RNA gel analysis of TAS4-siRNAs under different nutrient-deficient conditions. B, Small RNA gel analysis of TAS4-siR81(−) and miR828 in the shoots of wild-type (wt) plants and tas4, mir828, and PAP1/MYB75 activation (pap1-D) mutants. C, Expression of PAP1/MYB75, PAP2/MYB90, and MYB113 by quantitative RT-PCR analysis in tas4 and mir828 mutants. One of two biological replicates is presented, and the error bars indicate the sd of two technical replicates. D, Alteration of anthocyanin accumulation in tas4, mir828, and pap1-D mutants. Error bars represent the sd (n = 5). FW, Fresh weight; OD, optical density. E, A proposed model of the autoregulatory pathway of PAP1/MYB75 via miR828 and TAS4-siR81(−). ). TAS4-siR81(−) targets a group of MYB transcription factors, PAP1/MYB75, PAP2/MYB90, and MYB113, that are involved in the biosynthesis of anthocyanin (Rajagopalan et al., 2006). It is well documented that plants accumulate anthocyanin under −Pi, in which these MYB transcription factors are up-regulated (Misson et al., 2005). The negative regulatory role of siRNA seemed to be contradictory to the positive correlation between the expression of TAS4-siR81(−) and its target genes under the same conditions. This prompted us to carry out further investigation. TAS4 and MIR828 T-DNA knockout lines and a PAP1/MYB75 activation mutant, production of anthocyanin pigment1-Dominant (pap1-D; Borevitz et al., 2000), were examined. As predicted, there was no detectable TAS4-siR81(−) in the tas4 mutant, and neither miR828 nor TAS4-siR81(−) was detected in the mir828 mutant (Fig. 7B). The accumulation of miR828 and TAS4-siR81(−) was already observed in the +Pi shoots of pap1-D mutant and increased under −Pi (Fig. 7B). This suggests that PAP1/MYB75 can positively regulate MIR828 and/or TAS4 genes either by activating their expression directly or indirectly or by suppressing a negative regulator of MIR828 and/or TAS4 via a different feedback loop. Analyses of the 1-kb upstream sequences of MIR828 and TAS4 genes revealed multiple MYB binding sites (predicted by the PLACE database; http://www.dna.affrc.go.jp/PLACE/index.html) and one PAP1 cis-regulatory element, (C/T)(A/C)NCCAC(A/G/T)N(G/T) (Dare et al., 2008), at 305 and 747 bp upstream of the MIR828 hairpin precursor and TAS4, respectively. Thus, PAP1 likely acts directly on the MIR828 and TAS4 promoters to activate their transcription. On the other hand, the transcript level of PAP1/MYB75, PAP2/MYB90, and MYB113 was elevated in the tas4 and mir828 mutants under both +Pi and −Pi conditions (Fig. 7C). This could be explained by the relief of targeting of TAS4-siR81(−). Consistent with the increased expression of these MYB transcription factors, the anthocyanin accumulation was enhanced in the shoots of tas4 and mir828 mutants under both +Pi and −Pi conditions (Fig. 7D). These results uncovered an autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(−) in which an adequate expression level of PAP1/MYB75, PAP2/MYB90, and MYB113 and a proper accumulation of anthocyanin are maintained during −Pi (Fig. 7E). In this model, Pi deficiency results in up-regulation of those MYB transcription factors, which subsequently activate the biosynthesis of anthocyanin. On the other hand, the increased level of PAP1/MYB75 can also lead to production of TAS4-siR81(−) via the activation of miR828 and/or TAS4, which antagonizes its own expression and other MYB transcription factors. Whether the up-regulation of miR828 can be bypassed by the activation of PAP1/MYB75 remains to be studied. Because anthocyanin is a common stress pigment, we hypothesize that this autoregulation machinery may not be restricted to −Pi and could also be applied to other stress conditions in which anthocyanin is accumulated. miRNA399* and Its Movement In addition to miR399s, we observed a considerable amount of miR399* accumulated in −Pi tissues. The read number of different miR399* species was about half that of miR399 (e.g. miR399f*) or even higher than that of miR399 (e.g. miR399a* and miR399d*). Such high accumulation of miR399* could not be neglected because it was more abundant than many other miRNAs. RNA gel analyses confirmed the existence of miR399* with the high abundance in miR399d* and miR399f*, predominantly in the roots (Fig. 8A Figure 8. Open in new tabDownload slide Expression and long-distance movement of miR399*. A, Small RNA gel analysis of miR399*. B, Quantitative RT-PCR analysis of miR399f and miR399f* in the rootstocks of grafted plants. Combination of grafts is indicated as scion/rootstock. 399f, miR399f-overexpressing plants; wt, wild-type plants. The error bars represent sd (n = 4–5 independent grafted plants). Figure 8. Open in new tabDownload slide Expression and long-distance movement of miR399*. A, Small RNA gel analysis of miR399*. B, Quantitative RT-PCR analysis of miR399f and miR399f* in the rootstocks of grafted plants. Combination of grafts is indicated as scion/rootstock. 399f, miR399f-overexpressing plants; wt, wild-type plants. The error bars represent sd (n = 4–5 independent grafted plants). ), which was consistent with the high level of their primary transcripts and of mature miR399d and miR399f (Lin et al., 2008). We and others have previously shown that miR399 can move from shoots to roots and function as a long-distance signal in response to −Pi (Lin et al., 2008; Pant et al., 2008). We found that like miR399f, miR399f* could be detected in the wild-type rootstock grafted with the miR399f-overexpressing scions (Fig. 8B). This result indicates that both miR399f and miR399f* can move across the graft junction from scions to rootstocks. We hypothesize that miR399* may have certain physiological relevance, possibly functioning in concert with miR399. There may be several possible functions of miR399*. First, miR399* may be loaded into the RNA-induced silencing complex to target different genes. Many miRNA* species from Drosophila were reported to associate with the Argonaute protein and have inhibitory activity (Okamura et al., 2008). Several genes were predicted to be the putative targets of miR399* species; however, we were unable to validate them by RLM 5′-RACE analysis. Thus, whether these miR399* species possess gene-silencing activity remains uncertain. Second, given the ability of miR399 and miR399* to move across the graft junction, we suspect that miR399* may assist the long-distance movement of miR399 by forming a miR399/miR399* duplex in the phloem stream. Several miRNA* species were also detected in the phloem saps of Brassica napus (Buhtz et al., 2008; Pant et al., 2009). Although the authors concluded that they remain in the single-stranded form, based on the in vitro RNase digestion assay, they could not exclude the possibility of separation from a double-stranded RNA duplex during the RNA extraction. An in vivo nondestructive analysis needs to be developed to differentiate the double-stranded small RNA molecules from the single-stranded molecules. Alternatively, miR399* may regulate or buffer the targeting efficiency of miR399 on PHO2 mRNA via the formation of an RNA duplex. CONCLUSION In summary, we have presented an extensive survey of the small RNAs showing differential expression in response to −Pi. These Pi-responsive small RNAs and their target genes are likely involved in the development or regulation of adaptive responses to −Pi. Interestingly, while miR399, miR778, miR827, and miR2111 responded specifically to −Pi, miR156, miR169, miR395, and miR398 also responded to other nutrient stresses. These responses involve up-regulation or down-regulation of distinct miRNAs, which coordinate and balance the demands of different nutrients. We hypothesize that metabolic adjustment, carbon assimilation, or oxidative stresses may have an effect on the cross talk among different nutrients, as we observed here. Although there is still much more to learn about the biological significances of these Pi-responsive small RNAs, this work has opened a new avenue for the functional study of small RNA-mediated gene regulation in −Pi responses. MATERIALS AND METHODS Plant Materials and Growth Conditions Arabidopsis (Arabidopsis thaliana ecotype Columbia) wild-type plants and mutants (mir828 [SALK_097788], tas4 [SALK_066997], and pap1-D [CS3884]; Borevitz et al., 2000; Alonso et al., 2003) were grown hydroponically as described (Aung et al., 2006; Chiou et al., 2006). Seventeen-day-old seedlings were subjected to −Pi treatment by transferring them from +Pi (250 μm KH2PO4) into Pi-free half-strength Hoagland solution. Roots and shoots were collected separately after 5 d of treatment. Control samples were maintained in +Pi and collected at the same time. For treatments with different nutrient deficiency, the specific nutrient element was withdrawn from the nutrient solution for 5 d. To compensate for the ion originally associated with the omitted nutrient element, it was supplied in the form of chlorite salt at equivalent concentration except in the −Cu and −Fe treatments, in which they were omitted from the solution. RNA Preparation and Detection and High-Throughput Small RNA Sequencing Total RNA was isolated using Trizol reagent (Invitrogen). Four RNA samples from +Pi or −Pi roots or shoots of wild-type plants were used to generate small RNA libraries by the Solexa sequencing technology (Illumina). Small RNA gel-blot analysis was carried out as described (Chiou et al., 2006) with modification according to Pall et al. (2007). Quantitative RT-PCR analysis was performed as described (Lin et al., 2008). The sequences of probes and PCR primers are listed in Supplemental Table S5. Analysis of Sequencing Data The identical adapter-trimmed reads were grouped as unique sequences with associated counts of the individual sequence reads. Each unique sequence was aligned to the Arabidopsis genome (The Arabidopsis Information Resource 8; http://www.arabidopsis.org/), and only perfectly mapped sequences were retained and analyzed further (for details, see Supplemental Materials and Methods S1). We defined a sequence cluster as a group of same-strand sequences that were mapped to similar genomic locations, and the union of the mapped regions of the individual sequences of the cluster can form a consecutive nongapped region. For every cluster, the number of sequence reads from each of the four samples was counted separately and used for differential expression analysis. We used two approaches to perform the differential expression analysis, a cluster-based approach for finding novel differentially expressed gene models and a gene-model-based approach for analyzing known gene models, especially for miRNA genes and ta-siRNA genes. All sequence clusters or genes of interest were compared between two libraries for obtaining differentially expressed ones using a statistics test. For each of the differentially expressed sequence clusters that were not annotated, the sequence including the cluster union and 200-bp flanking genomic sequence extending from each side was extracted and folded by RNAfold (Hofacker, 2003). A sequence cluster was considered to be a novel miRNA if it met the following four criteria: (1) the structure forms a hairpin; (2) the sequence cluster is located in the duplex region of the hairpin structure; (3) there is another sequence cluster located on the other strand of the duplex; and (4) the most abundant sequences from each of the two clusters can form a duplex with two-nucleotide overhang at the 3′ end. The targets for the novel miRNAs were predicted using methods described elsewhere (Allen et al., 2005). Validation of the target gene was carried out by RLM 5′-RACE analysis using the GeneRacer kit (Invitrogen). Anthocyanin Analysis Anthocyanin analysis was performed according to Lange et al. (1971). The small RNA sequence data were submitted to the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/projects/geo/) under accession number GSE17741. Accession numbers of miRNAs are available at miRBase (http://microrna.sanger.ac.uk/sequences/index.shtml). Supplemental Data The following materials are available in the online version of this article. Supplemental Figure S1. Multiple processing of the miR778 precursor. Supplemental Figure S2. Expression of the target genes of miR395 and miR827. Supplemental Figure S3. Evolutionary analysis of MIR2111. Supplemental Table S1. Statistics of small RNA sequences. Supplemental Table S2. Categories of small RNA origins. Supplemental Table S3. tRNA-derived small RNAs. Supplemental Table S4. Expression of miRNA, ta-siRNA, and smRPi1LTR. Supplemental Table S5. Sequences of probes and primers. Supplemental Materials and Methods S1. ACKNOWLEDGMENTS We thank Shu-Hsing Wu and Ho-Ming Chen for helpful comments and Tzu-Yun Chang for assisting in the preparation of materials. LITERATURE CITED Achard P, Herr A, Baulcombe DC, Harberd NP ( 2004 ) Modulation of floral development by a gibberellin-regulated microRNA. Development 131 : 3357 – 3365 Crossref Search ADS PubMed Allen E, Xie Z, Gustafson AM, Carrington JC ( 2005 ) MicroRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell 121 : 207 – 221 Crossref Search ADS PubMed Alonso JM, Stepanova AN, Leisse TJ, Kim CJ, Chen H, Shinn P, Stevenson DK, Zimmerman J, Barajas P, Cheuk R, et al ( 2003 ) Genome-wide insertional mutagenesis of Arabidopsis thaliana. Science 301 : 653 – 657 Crossref Search ADS PubMed Aung K, Lin SI, Wu CC, Huang YT, Su CL, Chiou TJ ( 2006 ) pho2, a phosphate overaccumulator, is caused by a nonsense mutation in a microRNA399 target gene. Plant Physiol 141 : 1000 – 1011 Crossref Search ADS PubMed Bari R, Pant BD, Stitt M, Scheible WR ( 2006 ) PHO2, microRNA399, and PHR1 define a phosphate-signaling pathway in plants. Plant Physiol 141 : 988 – 999 Crossref Search ADS PubMed Borevitz JO, Xia Y, Blount J, Dixon RA, Lamb C ( 2000 ) Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell 12 : 2383 – 2394 Crossref Search ADS PubMed Brodersen P, Sakvarelidze-Achard L, Bruun-Rasmussen M, Dunoyer P, Yamamoto YY, Sieburth L, Voinnet O ( 2008 ) Widespread translational inhibition by plant miRNAs and siRNAs. Science 320 : 1185 – 1190 Crossref Search ADS PubMed Brodersen P, Voinnet O ( 2006 ) The diversity of RNA silencing pathways in plants. Trends Genet 22 : 268 – 280 Crossref Search ADS PubMed Buhtz A, Springer F, Chappell L, Baulcombe DC, Kehr J ( 2008 ) Identification and characterization of small RNAs from the phloem of Brassica napus. Plant J 53 : 739 – 749 Crossref Search ADS PubMed Chiou TJ ( 2007 ) The role of microRNAs in sensing nutrient stress. Plant Cell Environ 30 : 323 – 332 Crossref Search ADS PubMed Chiou TJ, Aung K, Lin SI, Wu CC, Chiang SF, Su Cl ( 2006 ) Regulation of phosphate homeostasis by microRNA in Arabidopsis. Plant Cell 18 : 412 – 421 Crossref Search ADS PubMed Combier JP, Frugier F, de Billy F, Boualem A, El-Yahyaoui F, Moreau S, Vernié T, Ott T, Gamas P, Crespi M, et al ( 2006 ) MtHAP2-1 is a key transcriptional regulator of symbiotic nodule development regulated by microRNA169 in Medicago truncatula. Genes Dev 20 : 3084 – 3088 Crossref Search ADS PubMed Dare AP, Schaffer RJ, Lin-Wang K, Allan AC, Hellens RP ( 2008 ) Identification of a cis-regulatory element by transient analysis of co-ordinately regulated genes. Plant Methods 4 : 17 Crossref Search ADS PubMed Duan K, Ding X, Zhang Q, Zhu H, Pan A, Huang J ( 2008 ) AtCopeg1, the unique gene originated from AtCopia95 retrotransposon family, is sensitive to external hormones and abiotic stresses. Plant Cell Rep 27 : 1065 – 1073 Crossref Search ADS PubMed Ebbs ML, Bender J ( 2006 ) Locus-specific control of DNA methylation by the Arabidopsis SUVH5 histone methyltransferase. Plant Cell 18 : 1166 – 1176 Crossref Search ADS PubMed Fahlgren N, Sullivan CM, Kasschau KD, Chapman EJ, Cumbie JS, Montgomery TA, Gilbert SD, Dasenko M, Backman TWH, Givan SA, et al ( 2009 ) Computational and analytical framework for small RNA profiling by high-throughput sequencing. RNA 15 : 992 – 1002 Crossref Search ADS PubMed Fujii H, Chiou TJ, Lin SI, Aung K, Zhu JK ( 2005 ) A miRNA involved in phosphate-starvation response in Arabidopsis. Curr Biol 15 : 2038 – 2043 Crossref Search ADS PubMed Gandikota M, Birkenbihl RP, Hohmann S, Cardon GH, Saedler H, Huijser P ( 2007 ) The miRNA156/157 recognition element in the 3′ UTR of the Arabidopsis SBP box gene SPL3 prevents early flowering by translational inhibition in seedlings. Plant J 49 : 683 – 693 Crossref Search ADS PubMed Haiser HJ, Karginov FV, Hannon GJ, Elliot MA ( 2008 ) Developmentally regulated cleavage of tRNAs in the bacterium Streptomyces coelicolor. Nucleic Acids Res 36 : 732 – 741 Crossref Search ADS PubMed Heisel SE, Zhang Y, Allen E, Guo L, Reynolds TL, Yang X, Kovalic D, Roberts JK ( 2008 ) Characterization of unique small RNA populations from rice grain. PLoS One 3 : e2871 Crossref Search ADS PubMed Hofacker IL ( 2003 ) Vienna RNA secondary structure server. Nucleic Acids Res 31 : 3429 – 3431 Crossref Search ADS PubMed Jochl C, Rederstorff M, Hertel J, Stadler PF, Hofacker IL, Schrettl M, Haas H, Huttenhofer A ( 2008 ) Small ncRNA transcriptome analysis from Aspergillus fumigatus suggests a novel mechanism for regulation of protein synthesis. Nucleic Acids Res 36 : 2677 – 2689 Crossref Search ADS PubMed Jones-Rhoades MW, Bartel DP ( 2004 ) Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell 14 : 787 – 799 Crossref Search ADS PubMed Jones-Rhoades MW, Bartel DP, Bartel B ( 2006 ) MicroRNAs and their regulatory roles in plants. Annu Rev Plant Biol 57 : 19 – 53 Crossref Search ADS PubMed Kawashima CG, Yoshimoto N, Maruyama-Nakashita A, Tsuchiya YN, Saito K, Takahashi H, Dalmay T ( 2009 ) Sulphur starvation induces the expression of microRNA-395 and one of its target genes but in different cell types. Plant J 57 : 313 – 321 Crossref Search ADS PubMed Lange H, Shropshire W Jr, Mohr H ( 1971 ) An analysis of phytochrome-mediated anthocyanin synthesis. Plant Physiol 47 : 649 – 655 Crossref Search ADS PubMed Lee SR, Collins K ( 2005 ) Starvation-induced cleavage of the tRNA anticodon loop in Tetrahymena thermophila. J Biol Chem 280 : 42744 – 42749 Crossref Search ADS PubMed Lenburg ME, O'Shea EK ( 1996 ) Signaling phosphate starvation. Trends Biochem Sci 21 : 383 – 387 Crossref Search ADS PubMed Li WX, Oono Y, Zhu J, He XJ, Wu JM, Iida K, Lu XY, Cui X, Jin H, Zhu JK ( 2008 a) The Arabidopsis NFYA5 transcription factor is regulated transcriptionally and posttranscriptionally to promote drought resistance. Plant Cell 20 : 2238 – 2251 Crossref Search ADS PubMed Li Y, Luo J, Zhou H, Liao JY, Ma LM, Chen YQ, Qu LH ( 2008 b) Stress-induced tRNA-derived RNAs: a novel class of small RNAs in the primitive eukaryote Giardia lamblia. Nucleic Acids Res 36 : 6048 – 6055 Crossref Search ADS PubMed Lin SI, Chiang SF, Lin WY, Chen JW, Tseng CY, Wu PC, Chiou TJ ( 2008 ) Regulatory network of microRNA399 and PHO2 by systemic signaling. Plant Physiol 147 : 732 – 746 Crossref Search ADS PubMed Lin WY, Lin SI, Chiou TJ ( 2009 ) Molecular regulators of phosphate homeostasis in plants. J Exp Bot 60 : 1427 – 1438 Crossref Search ADS PubMed Llave C, Xie Z, Kasschau KD, Carrington JC ( 2002 ) Cleavage of Scarecrow-like mRNA targets directed by a class of Arabidopsis miRNA. Science 297 : 2053 – 2056 Crossref Search ADS PubMed Lu C, Tej SS, Luo S, Haudenschild CD, Meyers BC, Green PJ ( 2005 ) Elucidation of the small RNA component of the transcriptome. Science 309 : 1567 – 1569 Crossref Search ADS PubMed Lu S, Sun YH, Chiang VL ( 2008 ) Stress-responsive microRNAs in Populus. Plant J 55 : 131 – 151 Crossref Search ADS PubMed Mallory AC, Vaucheret H ( 2006 ) Functions of microRNAs and related small RNAs in plants. Nat Genet 38 : S31 –S36 Crossref Search ADS PubMed Marschner H ( 1995 ). Mineral Nutrition of Higher Plants. Academic Press, London, pp 265–277, 537–594 Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, et al ( 2008 ) Criteria for annotation of plant microRNAs. Plant Cell 20 : 3186 – 3190 Crossref Search ADS PubMed Misson J, Raghothama KG, Jain A, Jouhet J, Block MA, Bligny R, Ortet P, Creff A, Somerville S, Rolland N, et al ( 2005 ) A genome-wide transcriptional analysis using Arabidopsis thaliana Affymetrix gene chips determined plant responses to phosphate deprivation. Proc Natl Acad Sci USA 102 : 11934 – 11939 Crossref Search ADS PubMed Navarro L, Dunoyer P, Jay F, Arnold B, Dharmasiri N, Estelle M, Voinnet O, Jones JDG ( 2006 ) A plant miRNA contributes to antibacterial resistance by repressing auxin signaling. Science 312 : 436 – 439 Crossref Search ADS PubMed Nobuta K, Venu RC, Lu C, Beló A, Vemaraju K, Kulkarni K, Wang W, Pillay M, Green PJ, Wang GL, et al ( 2007 ) An expression atlas of rice mRNAs and small RNAs. Nat Biotechnol 25 : 473 – 477 Crossref Search ADS PubMed Okamura K, Phillips MD, Tyler DM, Duan H, Chou YT, Lai EC ( 2008 ) The regulatory activity of microRNA* species has substantial influence on microRNA and 3′ UTR evolution. Nat Struct Mol Biol 15 : 354 – 363 Crossref Search ADS PubMed Pall GS, Codony-Servat C, Byrne J, Ritchie L, Hamilton A ( 2007 ) Carbodiimide-mediated cross-linking of RNA to nylon membranes improves the detection of siRNA, miRNA and piRNA by northern blot. Nucleic Acids Res 35 : e60 Crossref Search ADS PubMed Pant BD, Buhtz A, Kehr J, Scheible WR ( 2008 ) MicroRNA399 is a long-distance signal for the regulation of plant phosphate homeostasis. Plant J 53 : 731 – 738 Crossref Search ADS PubMed Pant BD, Musialak-Lange M, Nuc P, May P, Buhtz A, Kehr J, Walther D, Scheible WR ( 2009 ) Identification of nutrient-responsive Arabidopsis and rapeseed microRNAs by comprehensive real-time polymerase chain reaction profiling and small RNA sequencing. Plant Physiol 150 : 1541 – 1555 Crossref Search ADS PubMed Peng M, Hannam C, Gu H, Bi YM, Rothstein SJ ( 2007 ) A mutation in NLA, which encodes a RING-type ubiquitin ligase, disrupts the adaptability of Arabidopsis to nitrogen limitation. Plant J 50 : 320 – 337 Crossref Search ADS PubMed Peragine A, Yoshikawa M, Wu G, Albrecht HL, Poethig RS ( 2004 ) SGS3 and SGS2/SDE1/RDR6 are required for juvenile development and the production of trans-acting siRNAs in Arabidopsis. Genes Dev 18 : 2368 – 2379 Crossref Search ADS PubMed Poirier Y, Bucher M ( 2002 ) Phosphate transport and homeostasis in Arabidopsis. In CR Somerville, EM Meyerowitz, eds, The Arabidopsis Book. American Society of Plant Biologists, Rockville, MD, doi: 10.1199/tab.0024 , http://www.aspb.org/publications/arabidopsis/ Crossref Raghothama KG ( 1999 ) Phosphate acquisition. Annu Rev Plant Physiol Plant Mol Biol 50 : 665 – 693 Crossref Search ADS PubMed Rajagopalan R, Vaucheret H, Trejo J, Bartel DP ( 2006 ) A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev 20 : 3407 – 3425 Crossref Search ADS PubMed Rubio V, Linhares F, Solano R, Martin AC, Iglesias J, Leyva A, Paz-Ares J ( 2001 ) A conserved MYB transcription factor involved in phosphate starvation signaling both in vascular plants and in unicellular algae. Genes Dev 15 : 2122 – 2133 Crossref Search ADS PubMed Schachtman DP, Shin R ( 2007 ) Nutrient sensing and signaling: NPKS. Annu Rev Plant Biol 58 : 47 – 69 Crossref Search ADS PubMed Schwarz S, Grande A, Bujdoso N, Saedler H, Huijser P ( 2008 ) The microRNA regulated SBP-box genes SPL9 and SPL15 control shoot maturation in Arabidopsis. Plant Mol Biol 67 : 183 – 195 Crossref Search ADS PubMed Shin R, Berg RH, Schachtman DP ( 2005 ) Reactive oxygen species and root hairs in Arabidopsis root response to nitrogen, phosphorus and potassium deficiency. Plant Cell Physiol 46 : 1350 – 1357 Crossref Search ADS PubMed Sunkar R, Kapoor A, Zhu JK ( 2006 ) Posttranscriptional induction of two Cu/Zn superoxide dismutase genes in Arabidopsis is mediated by downregulation of miR398 and important for oxidative stress tolerance. Plant Cell 18 : 2051 – 2065 Crossref Search ADS PubMed Sunkar R, Zhu JK ( 2004 ) Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell 16 : 2001 – 2019 Crossref Search ADS PubMed Thompson DM, Lu C, Green PJ, Parker R ( 2008 ) tRNA cleavage is a conserved response to oxidative stress in eukaryotes. RNA 14 : 2095 – 2103 Crossref Search ADS PubMed Thompson DM, Parker R ( 2009 ) Stressing out over tRNA cleavage. Cell 138 : 215 – 219 Crossref Search ADS PubMed Vaucheret H ( 2006 ) Post-transcriptional small RNA pathways in plants: mechanisms and regulations. Genes Dev 20 : 759 – 771 Crossref Search ADS PubMed Vazquez F ( 2006 ) Arabidopsis endogenous small RNAs: highways and byways. Trends Plant Sci 11 : 460 – 468 Crossref Search ADS PubMed Vazquez F, Vaucheret H, Rajagopalan R, Lepers C, Gasciolli V, Mallory AC, Hilbert JL, Bartel DP, Crete P ( 2004 ) Endogenous trans-acting siRNAs regulate the accumulation of Arabidopsis mRNAs. Mol Cell 16 : 69 – 79 Crossref Search ADS PubMed Voinnet O ( 2009 ) Origin, biogenesis, and activity of plant microRNAs. Cell 136 : 669 – 687 Crossref Search ADS PubMed Wang JW, Schwab R, Czech B, Mica E, Weigel D ( 2008 ) Dual effects of miR156-targeted SPL genes and CYP78A5/KLUH on plastochron length and organ size in Arabidopsis thaliana. Plant Cell 20 : 1231 – 1243 Crossref Search ADS PubMed Wu G, Poethig RS ( 2006 ) Temporal regulation of shoot development in Arabidopsis thaliana by miR156 and its target SPL3. Development 133 : 3539 – 3547 Crossref Search ADS PubMed Yaeno T, Iba K ( 2008 ) BAH1/NLA, a RING-type ubiquitin E3 ligase, regulates the accumulation of salicylic acid and immune responses to Pseudomonas syringae DC3000. Plant Physiol 148 : 1032 – 1041 Crossref Search ADS PubMed Yamasaki H, Abdel-Ghany SE, Cohu CM, Kobayashi Y, Shikanai T, Pilon M ( 2007 ) Regulation of copper homeostasis by micro-RNA in Arabidopsis. J Biol Chem 282 : 16369 – 16378 Crossref Search ADS PubMed Zhang S, Sun L, Kragler F ( 2009 ) The phloem delivered RNA pool contains small non-coding RNAs and interferes with translation. Plant Physiol 150 : 378 – 387 Crossref Search ADS PubMed Zhu QH, Spriggs A, Matthew L, Fan L, Kennedy G, Gubler F, Helliwell C ( 2008 ) A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains. Genome Res 18 : 1456 – 1465 Crossref Search ADS PubMed Author notes 1 This work was supported by Academia Sinica, Taiwan (grant nos. AS97FPL20–1 and AS98CDAL11), and the National Science Council, Taiwan (grant no. 97–2321–B–001–018). 2 These authors contributed equally to the article. * Corresponding author; e-mail tjchiou@gate.sinica.edu.tw. The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: Tzyy-Jen Chiou (tjchiou@gate.sinica.edu.tw). [W] The online version of this article contains Web-only data. [OA] Open Access articles can be viewed online without a subscription. www.plantphysiol.org/cgi/doi/10.1104/pp.109.147280 © 2009 American Society of Plant Biologists © The Author(s) 2009. Published by Oxford University Press on behalf of American Society of Plant Biologists. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Journal

Plant PhysiologyOxford University Press

Published: Dec 4, 2009

There are no references for this article.