Search

Noncoding and coding RNAs are key regulators of plant growth, development, and stress responses. To investigate the types of transcripts accumulated during the vegetative to reproductive transition and floral development in the Coffea arabica L., we sequenced small RNA libraries from eight developmental stages, up to anthesis. We combined these data with messenger RNA and PARE sequencing of two important development stages that marks the transition of an apparent latent to a rapid growth stage. In addition, we took advantage of multiple in silico tools to characterize genomic loci producing small RNAs such as phasiRNAs, miRNAs, and tRFs. Our differential and co-expression analysis showed that some types of small RNAs such as tRNAs, snoRNAs, snRNAs, and phasiRNAs preferentially accumulate in a stage-specific manner. Members of the miR482/miR2118 superfamily and their 21-nucleotide phasiRNAs originating from resistance genes show a robust co-expression pattern that is maintained across all the evaluated developmental stages. Finally, the majority of miRNAs accumulate in a family stage-specific manner, related to modulated hormonal responses and transcription factor expression., (© 2024 Society for Experimental Biology and John Wiley & Sons Ltd.)

Small RNAs (sRNAs), found extensively in plants, play an essential role in plant growth and development. Although various sRNA analysis tools have been developed for plants, the use of most of them depends on programming and command-line environments, which is a challenge for many wet-lab biologists. Furthermore, current sRNA analysis tools mostly focus on the analysis of certain type of sRNAs and are resource-intensive, normally demanding an immense amount of time and effort to learn the use of numerous tools or scripts and assemble them into a workable pipeline to get the final results. Here, we present sRNAminer, a powerful stand-alone toolkit with a user-friendly interface that integrates all common functions for the analysis of three major types of plant sRNAs: microRNAs (miRNAs), phased small interfering RNAs (phasiRNAs), and heterochromatic siRNAs (hc-siRNAs). We constructed a curated or "golden" set of MIRNA and PHAS loci, which was used to assess the performance of sRNAminer in comparison to other existing tools. The results showed that sRNAminer outperformed these tools in multiple aspects, highlighting its functionality. In addition, to enable an efficient evaluation of sRNA annotation results, we developed Integrative Genomics Viewer (IGV)-sRNA, a modified genome browser optimized from IGV and we incorporated it as a functional module in sRNAminer. IGV-sRNA can display a wealth of sRNA-specific features, enabling a more comprehensive understanding of sRNA data. sRNAminer and IGV-sRNA are both platform-independent software that can be run under all operating systems. They are now freely available at https://github.com/kli28/sRNAminer and https://gitee.com/CJchen/IGV-sRNA., Competing Interests: Conflict of interest The authors declare that they have no conflict of interest., (Copyright © 2023 Science China Press. Published by Elsevier B.V. All rights reserved.)

Several protein families participate in the biogenesis and function of small RNAs (sRNAs) in plants. Those with primary roles include Dicer-like (DCL), RNA-dependent RNA polymerase (RDR), and Argonaute (AGO) proteins. Protein families such as double-stranded RNA-binding (DRB), SERRATE (SE), and SUPPRESSION OF SILENCING 3 (SGS3) act as partners of DCL or RDR proteins. Here, we present curated annotations and phylogenetic analyses of seven sRNA pathway protein families performed on 196 species in the Viridiplantae (aka green plants) lineage. Our results suggest that the RDR3 proteins emerged earlier than RDR1/2/6. RDR6 is found in filamentous green algae and all land plants, suggesting that the evolution of RDR6 proteins coincides with the evolution of phased small interfering RNAs (siRNAs). We traced the origin of the 24-nt reproductive phased siRNA-associated DCL5 protein back to the American sweet flag (Acorus americanus), the earliest diverged, extant monocot species. Our analyses of AGOs identified multiple duplication events of AGO genes that were lost, retained, or further duplicated in subgroups, indicating that the evolution of AGOs is complex in monocots. The results also refine the evolution of several clades of AGO proteins, such as AGO4, AGO6, AGO17, and AGO18. Analyses of nuclear localization signal sequences and catalytic triads of AGO proteins shed light on the regulatory roles of diverse AGOs. Collectively, this work generates a curated and evolutionarily coherent annotation for gene families involved in plant sRNA biogenesis/function and provides insights into the evolution of major sRNA pathways., Competing Interests: Conflict of interest statement. None declared., (© American Society of Plant Biologists 2023. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Plant cells accumulate small RNA molecules that regulate plant development, genome stability, and environmental responses. These small RNAs fall into three major classes based on their function and mechanisms of biogenesis-microRNAs, heterochromatic small interfering RNAs, and secondary small interfering RNAs-plus several other less well-characterized categories. Biogenesis of each small RNA class requires a pathway of factors, some specific to each pathway and others involved in multiple pathways. Diverse sequenced plant genomes, along with rapid developments in sequencing, imaging, and genetic transformation techniques, have enabled significant progress in understanding the biogenesis, functions, and evolution of plant small RNAs, including those that had been poorly characterized because they were absent or had low representation in Arabidopsis ( Arabidopsis thaliana ). Here, we review recent findings about plant small RNAs and discuss our current understanding of their biogenesis mechanisms, targets, modes of action, mobility, and functions in Arabidopsis and other plant species, including economically important crops.

Plant microRNAs (miRNAs) are short, noncoding RNA molecules that restrict gene expression via posttranscriptional regulation and function in several essential pathways, including development, growth, and stress responses. Accurately identifying miRNAs in populations of small RNA sequencing libraries is a computationally intensive process that has resulted in the misidentification of inaccurately annotated miRNA sequences. In recent years, criteria for miRNA annotation have been refined with the aim to reduce these misannotations. Here, we describe miRador, a miRNA identification tool that utilizes the most up-to-date, community-established criteria for accurate identification of miRNAs in plants. We combined target prediction and Parallel Analysis of RNA Ends (PARE) data to assess the precision of the miRNAs identified by miRador. We compared miRador to other commonly used miRNA prediction tools and found that miRador is at least as precise as other prediction tools while being substantially faster than other tools. miRador should be broadly useful for the plant community to identify and annotate miRNAs in plant genomes., Competing Interests: Conflict of interest statement. None declared., (© American Society of Plant Biologists 2022. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

The anther-enriched phased, small interfering RNAs (phasiRNAs) play vital roles in sustaining male fertility in grass species. Their long non-coding precursors are synthesized by RNA polymerase II and are likely regulated by transcription factors (TFs). A few putative transcriptional regulators of the 21- or 24-nucleotide phasiRNA loci (referred to as 21- or 24-PHAS loci) have been identified in maize (Zea mays), but whether any of the individual TFs or TF combinations suffice to activate any PHAS locus is unclear. Here, we identified the temporal gene coexpression networks (modules) associated with maize anther development, including two modules highly enriched for the 21- or 24-PHAS loci. Comparisons of these coexpression modules and gene sets dysregulated in several reported male sterile TF mutants provided insights into TF timing with regard to phasiRNA biogenesis, including antagonistic roles for OUTER CELL LAYER4 and MALE STERILE23. Trans-activation assays in maize protoplasts of individual TFs using bulk-protoplast RNA-sequencing showed that two of the TFs coexpressed with 21-PHAS loci could activate several 21-nucleotide phasiRNA pathway genes but not transcription of 21-PHAS loci. Screens for combinatorial activities of these TFs and, separately, the recently reported putative transcriptional regulators of 24-PHAS loci using single-cell (protoplast) RNA-sequencing, did not detect reproducible activation of either 21-PHAS or 24-PHAS loci. Collectively, our results suggest that the endogenous transcriptional machineries and/or chromatin states in the anthers are necessary to activate reproductive PHAS loci., (© 2022 Society for Experimental Biology and John Wiley & Sons Ltd.)

Phased secondary siRNAs (phasiRNAs) are broadly present in the reproductive tissues of flowering plants, with spatial-temporal specificity. However, the ARGONAUTE (AGO) proteins associated with phasiRNAs and their miRNA triggers remain elusive. Here, through histological and high-throughput sequencing analyses, we show that rice AGO1d, which is specifically expressed in anther wall cells before and during meiosis, associates with both miR2118 and miR2275 to mediate phasiRNA biogenesis. AGO1d preferentially binds to miR2118-triggered 21-nucleotide (nt) phasiRNAs with a 5'-terminal uridine, suggesting a dual role in phasiRNA biogenesis and function. Depletion of AGO1d causes a reduction of 21- and 24-nt phasiRNAs and temperature-sensitive male sterility. At lower temperatures, anthers of the ago1d mutant predominantly show excessive tapetal cells with little starch accumulation during pollen formation, possibly caused by the dysregulation of cell metabolism. These results uncover an essential role of AGO1d in rice anther development at lower temperatures and demonstrate coordinative roles of AGO proteins during reproductive phasiRNA biogenesis and function., (© 2022 The Authors. New Phytologist © 2022 New Phytologist Foundation.)

The Solanaceae or "nightshade" family is an economically important group with remarkable diversity. To gain a better understanding of how the unique biology of the Solanaceae relates to the family's small RNA (sRNA) genomic landscape, we downloaded over 255 publicly available sRNA data sets that comprise over 2.6 billion reads of sequence data. We applied a suite of computational tools to predict and annotate two major sRNA classes: (1) microRNAs (miRNAs), typically 20- to 22-nucleotide (nt) RNAs generated from a hairpin precursor and functioning in gene silencing and (2) short interfering RNAs (siRNAs), including 24-nt heterochromatic siRNAs typically functioning to repress repetitive regions of the genome via RNA-directed DNA methylation, as well as secondary phased siRNAs and trans-acting siRNAs generated via miRNA-directed cleavage of a polymerase II-derived RNA precursor. Our analyses described thousands of sRNA loci, including poorly understood clusters of 22-nt siRNAs that accumulate during viral infection. The birth, death, expansion, and contraction of these sRNA loci are dynamic evolutionary processes that characterize the Solanaceae family. These analyses indicate that individuals within the same genus share similar sRNA landscapes, whereas comparisons between distinct genera within the Solanaceae reveal relatively few commonalities., (© The Author(s) 2022. Published by Oxford University Press on behalf of American Society of Plant Biologists.)

RNA silencing is a conserved mechanism in eukaryotes involved in development and defense against viruses. In plants, ARGONAUTE1 (AGO1) protein plays a central role in both microRNA- and small interfering RNA-directed silencing, and its expression is regulated at multiple levels. Here, we report that the F-box protein FBW2 assembles an SCF complex that selectively targets for proteolysis AGO1 when it is unloaded and mutated. Although FBW2 loss of function does not lead to strong growth or developmental defects, it significantly increases RNA-silencing activity. Interestingly, under conditions in which small-RNA accumulation is affected, the failure to degrade AGO1 in fbw2 mutants becomes more deleterious for the plant. Accordingly, the non-degradable AGO1 protein assembles high-molecular-weight complexes and binds illegitimate small RNA, leading to off-target cleavage. Therefore, control of AGO1 homeostasis by FBW2 plays an important role in quality control of RNA silencing., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

Small RNAs play important roles in plant growth and development by modulating expression of genes and transposons. In many flowering plant species, male reproductive organs, the anthers, produce abundant phased small interfering RNAs (phasiRNAs). Two classes of reproductive phasiRNAs are generally known, mostly from monocots: (i) pre-meiotic 21-nucleotide (nt) phasiRNAs triggered by miR2118 and (ii) meiotic 24-nt phasiRNAs triggered by miR2275. Here, we describe conserved and non-conserved triggers of 24-nt phasiRNAs in several eudicots. We found that the abundant 24-nt phasiRNAs in the basal eudicot columbine (Aquilegia coerulea) are produced by the canonical trigger miR2275, as well as by other non-canonical triggers, miR482/2118 and miR14051. These triggering microRNAs (miRNAs) are localized in microspore mother cells and tapetal cells of meiotic and post-meiotic stage anthers. Furthermore, we identified a lineage-specific trigger (miR11308) of 24-nt phasiRNAs and an expanded number of 24-PHAS loci in wild strawberry (Fragaria vesca). We validated the presence of the miR2275-derived 24-nt phasiRNA pathway in rose (Rosa chinensis). Finally, we evaluated all eudicots that have been validated for the presence of 24-nt phasiRNAs as possible model systems in which to study the biogenesis and function of 24-nt phasiRNAs. We conclude that columbine (Aquilegia coerulea) would be a strong model because of its extensive number of 24-PHAS loci and its diversity of trigger miRNAs., (© 2021 Society for Experimental Biology and John Wiley & Sons Ltd.)

In grasses, two types of phased, small interfering RNAs (phasiRNAs) are expressed largely in young, developing anthers. They are 21 or 24 nucleotides (nt) in length and are triggered by miR2118 or miR2275, respectively. However, most of their functions and activities are not fully understood. We performed comparative genomic analysis of their source loci (PHAS) in five Oryza genomes and combined this with analysis of high-throughput sRNA and degradome datasets. In total, we identified 8216 21-PHAS and 626 24-PHAS loci. Local tandem and segmental duplications mainly contributed to the expansion and supercluster distribution of the 21-PHAS loci. Despite their relatively conserved genomic positions, PHAS sequences diverged rapidly, except for the miR2118/2275 target sites, which were under strong selection for conservation. We found that 21-nt phasiRNAs with a 5'-terminal uridine (U) demonstrated cis-cleavage at PHAS precursors, and these cis-acting sites were also variable among close species. miR2118 could trigger phasiRNA production from its own antisense transcript and the derived phasiRNAs might reversibly regulate miR2118 precursors. We hypothesised that successful initiation of phasiRNA biogenesis is conservatively maintained, while phasiRNA products diverged quickly and are not individually conserved. In particular, phasiRNA production is under the control of multiple reciprocal regulation mechanisms., (© 2020 The Authors New Phytologist © 2020 New Phytologist Foundation.)

MicroRNAs (miRNAs) are 20-24 nucleotides (nt) small RNAs functioning in eukaryotes. The length and sequence of miRNAs are not only related to the biogenesis of miRNAs but are also important for downstream physiological processes like ta-siRNA production. To investigate these roles, it is informative to create small mutations within mature miRNA sequences. We used both TALENs (transcription activator-like effector nucleases) and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) to introduce heritable base pair mutations in mature miRNA sequences. For rice, TALEN constructs were built targeting five different mature miRNA sequences and yielding heritable mutations. Among the resulting mutants, mir390 mutant showed a severe defect in the shoot apical meristem (SAM), a shootless phenotype, which could be rescued by the wild-type MIR390. Small RNA sequencing showed the two base pair deletion in mir390 substantially interfered with miR390 biogenesis. In Arabidopsis, CRISPR/Cas9-mediated editing of the miR160* strand confirmed that the asymmetric structure of miRNA is not a necessary determinant for secondary siRNA production. CRISPR/Cas9 with double-guide RNAs successfully generated mir160a null mutants with fragment deletions, at a higher efficiency than a single-guide RNA. The difference between the phenotypic severity of miR160a mutants in Col-0 versus Ler backgrounds highlights a diverged role for miR160a in different ecotypes. Overall, we demonstrated that TALENs and CRISPR/Cas9 are both effective in modifying miRNA precursor structure, disrupting miRNA processing and generating miRNA null mutant plants., (© 2019 The Authors. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

Maize somatic embryogenesis (SE) requires the induction of embryogenic callus and establishment of proliferation before plant regeneration. The molecular mechanisms underlying callus embryogenic potential are not well understood. Here we explored the role of small RNAs (sRNAs) and the accumulation of their target transcripts in maize SE at the dedifferentiation step using VS-535 zygotic embryos collected at distinct developmental stages and displaying contrasting in vitro embryogenic potential and morphology. MicroRNAs (miRNAs), trans-acting siRNAs (tasiRNAs), heterochromatic siRNAs (hc-siRNAs) populations and their RNA targets were analyzed by high-throughput sequencing. Abundances of specific miRNAs, tasiRNAs and targets were validated by qRT-PCR. Unique accumulation patterns were found for immature embryo at 15 Days After Pollination (DAP) and for the callus induction from this explant, as compared to 23 DAP and mature embryos. miR156, miR164, miR166, tasiARFs and the 24 nt hc-siRNAs displayed the most strikingly different patterns between explants and during dedifferentiation. According to their role in auxin responses and developmental cues, we conclude that sRNA-target regulation operating within the 15 DAP immature embryo explant provides key molecular hints as to why this stage is relevant for callus induction with successful proliferation and plant regeneration.

The Lemnaceae family comprises aquatic plants of angiosperms gaining attention due to their utility in wastewater treatment, and rapid production of biomass that can be used as feed, fuel, or food. Moreover, it can serve as a model species for neotenous growth and environmental adaptation. The latter properties are subject to post-transcriptional regulation of gene expression, meriting investigation of how miRNAs in Spirodela polyrhiza, the most basal and most thoroughly sequenced member of the family, are expressed under different growth conditions. To further scientific understanding of its capacity to adapt to environmental cues, we measured miRNA expression and processing of their target sequences under different temperatures, and in the presence of abscisic acid, copper, kinetin, nitrate, and sucrose. Using two small RNA sequencing experiments and one degradome sequencing experiment, we provide evidence for 108 miRNAs. Sequencing cleaved mRNAs validated 42 conserved miRNAs with 83 targets and 24 novel miRNAs regulating 66 targets and created a list of 575 predicted and verified targets. These analyses revealed condition-induced changes in miRNA expression and cleavage activity, and resulted in the addition of stringently reviewed miRNAs to miRBase. This combination of small RNA and degradome sequencing provided not only high confidence predictions of conserved and novel miRNAs and targets, but also a view of the post-transcriptional regulation of adaptations. A unique aspect is the role of miR156 and miR172 expression and activity in its clonal propagation and neoteny. Additionally, low levels of 24 nt sRNAs were observed, despite the lack of recent retrotransposition., (© 2019 The Authors The Plant Journal © 2019 John Wiley & Sons Ltd.)

Many evolutionarily conserved microRNAs (miRNAs) in plants regulate transcription factors with key functions in development. Hence, mutations in the core components of the miRNA biogenesis machinery cause strong growth defects. An essential aspect of miRNA biogenesis is the precise excision of the small RNA from its precursor. In plants, miRNA precursors are largely variable in size and shape and can be processed by different modes. Here, we optimized an approach to detect processing intermediates during miRNA biogenesis. We characterized a miRNA whose processing is triggered by a terminal branched loop. Plant miRNA processing can be initiated by internal bubbles, small terminal loops or branched loops followed by dsRNA segments of 15-17 bp. Interestingly, precision and efficiency vary with the processing modes. Despite the various potential structural determinants present in a single a miRNA precursor, DCL1 is mostly guided by a predominant structural region in each precursor in wild-type plants. However, our studies in fiery1, hyl1 and se mutants revealed the existence of cleavage signatures consistent with the recognition of alternative processing determinants. The results provide a general view of the mechanisms underlying the specificity of miRNA biogenesis in plants.

Phased, secondary siRNAs (phasiRNAs) represent a class of small RNAs in plants generated via distinct biogenesis pathways, predominantly dependent on the activity of 22-nt miRNAs. Most 22-nt miRNAs are processed by DCL1 from miRNA precursors containing an asymmetric bulge, yielding a 22/21-nt miRNA/miRNA* duplex. Here we show that miR1510, a soybean miRNA capable of triggering phasiRNA production from numerous nucleotide-binding leucine-rich repeat ( NB-LRR s), previously described as 21 nt in its mature form, primarily accumulates as a 22-nt isoform via monouridylation. We demonstrate that, in Arabidopsis , this uridylation is performed by HESO1. Biochemical experiments showed that the 3' terminus of miR1510 is only partially 2'- O -methylated because of the terminal mispairing in the miR1510/miR1510* duplex that inhibits HEN1 activity in soybean. miR1510 emerged in the Phaseoleae ∼41-42 million years ago with a conserved precursor structure yielding a 22-nt monouridylated form, yet a variant in mung bean is processed directly in a 22-nt mature form. This analysis of miR1510 yields two observations: ( i ) plants can utilize postprocessing modification to generate abundant 22-nt miRNA isoforms to more efficiently regulate target mRNA abundances; and ( ii ) comparative analysis demonstrates an example of selective optimization of precursor processing of a young plant miRNA., Competing Interests: The authors declare no conflict of interest.

Vegetative phase change is regulated by a decrease in the abundance of the miRNAs, miR156 and miR157, and the resulting increase in the expression of their targets, SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors. To determine how miR156/miR157 specify the quantitative and qualitative changes in leaf morphology that occur during vegetative phase change, we measured their abundance in successive leaves and characterized the phenotype of mutations in different MIR156 and MIR157 genes. miR156/miR157 decline rapidly between leaf 1&2 and leaf 3 and decrease more slowly after this point. The amount of miR156/miR157 in leaves 1&2 greatly exceeds the threshold required to specify their identity. Subsequent leaves have relatively low levels of miR156/miR157 and are sensitive to small changes in their abundance. In these later-formed leaves, the amount of miR156/miR157 is close to the threshold required to specify juvenile vs. adult identity; a relatively small decrease in the abundance of miR156/157 in these leaves produces a disproportionately large increase in SPL proteins and a significant change in leaf morphology. miR157 is more abundant than miR156 but has a smaller effect on shoot morphology and SPL gene expression than miR156. This may be attributable to the inefficiency with which miR157 is loaded onto AGO1, as well as to the presence of an extra nucleotide at the 5' end of miR157 that is mis-paired in the miR157:SPL13 duplex. miR156 represses different targets by different mechanisms: it regulates SPL9 by a combination of transcript cleavage and translational repression and regulates SPL13 primarily by translational repression. Our results offer a molecular explanation for the changes in leaf morphology that occur during shoot development in Arabidopsis and provide new insights into the mechanism by which miR156 and miR157 regulate gene expression.

MicroRNAs (miRNAs) are key regulators of gene expression. A handful of miRNAs are broadly conserved in land plants, while the majority are lineage specific; this review describes the processes by which new miRNAs are hypothesized to have emerged. Two major models describe miRNA origins, firstly, de novo emergence via inverted duplication of target gene fragments, and secondly, the expansion and neofunctionalization of existing miRNA families. The occasional acquisition of target sites by previously un-targeted genes adds further dynamism to the process by which miRNAs may shift roles during evolution. Additional factors guiding miRNA evolution include functional constraints on their length and the importance of precursor conservation that is observed in regions above or below the mature miRNA duplex; these regions represent recognition sites for components of biogenesis machinery and direct precursor processing. Insights into the mechanisms of miRNA emergence and divergence are important for understanding plant genome evolution and the impact of miRNA regulatory networks., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

MicroRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs) play vital regulatory roles in plant growth and development. Little is known about these small RNAs in litchi (Litchi chinensis), an economically important fruit crop widely cultivated in Southeast Asia. We profiled the litchi small RNA population with various deep-sequencing techniques and in-depth bioinformatic analyses. The genome-wide identification of miRNAs, their target genes, and phasiRNA-generating (PHAS) genes/loci showed that the function of miR482/2118 has expanded, relative to its canonical function. We also discovered that, for 29 PHAS loci, miRNA-mediated phasiRNA production was coupled with alternative splicing (AS) and alternative polyadenylation (APA). Most of these loci encoded long noncoding RNAs. An miR482/2118 targeted locus gave rise to four main transcript isoforms through AS/APA, and diverse phasiRNAs generated from these isoforms appeared to target long terminal repeat (LTR) retrotransposons and other unrelated genes. This coupling enables phasiRNA production from different exons of noncoding PHAS genes and yields diverse phasiRNA populations, both broadening and altering the range of downstream phasiRNA-regulated genes. Our results reveal the diversity of miRNA and phasiRNA in litchi, and demonstrate AS/APA as a new layer of regulation in small RNA-mediated gene silencing., (© 2017 The Authors. New Phytologist © 2017 New Phytologist Trust.)

MicroRNAs (miRNAs) are ∼21-nucleotide-long regulatory RNAs that arise from endonucleolytic processing of hairpin precursors. Many function as essential posttranscriptional regulators of target mRNAs and long noncoding RNAs. Alongside miRNAs, plants also produce large numbers of short interfering RNAs (siRNAs), which are distinguished from miRNAs primarily by their biogenesis (typically processed from long double-stranded RNA instead of single-stranded hairpins) and functions (typically via roles in transcriptional regulation instead of posttranscriptional regulation). Next-generation DNA sequencing methods have yielded extensive data sets of plant small RNAs, resulting in many miRNA annotations. However, it has become clear that many miRNA annotations are questionable. The sheer number of endogenous siRNAs compared with miRNAs has been a major factor in the erroneous annotation of siRNAs as miRNAs. Here, we provide updated criteria for the confident annotation of plant miRNAs, suitable for the era of "big data" from DNA sequencing. The updated criteria emphasize replication and the minimization of false positives, and they require next-generation sequencing of small RNAs. We argue that improved annotation systems are needed for miRNAs and all other classes of plant small RNAs. Finally, to illustrate the complexities of miRNA and siRNA annotation, we review the evolution and functions of miRNAs and siRNAs in plants., (© American Society of Plant Biologists. All rights reserved.)

The altered carbon assimilation pathway of crassulacean acid metabolism (CAM) photosynthesis results in an up to 80% higher water-use efficiency than C 3 photosynthesis in plants making it a potentially useful pathway for engineering crop plants with improved drought tolerance. Here we surveyed detailed temporal (diel time course) and spatial (across a leaf gradient) gene and microRNA (miRNA) expression patterns in the obligate CAM plant pineapple [Ananas comosus (L.) Merr.]. The high-resolution transcriptome atlas allowed us to distinguish between CAM-related and non-CAM gene copies. A differential gene co-expression network across green and white leaf diel datasets identified genes with circadian oscillation, CAM-related functions, and source-sink relations. Gene co-expression clusters containing CAM pathway genes are enriched with clock-associated cis-elements, suggesting circadian regulation of CAM. About 20% of pineapple microRNAs have diel expression patterns, with several that target key CAM-related genes. Expression and physiology data provide a model for CAM-specific carbohydrate flux and long-distance hexose transport. Together these resources provide a list of candidate genes for targeted engineering of CAM into C 3 photosynthesis crop species., (© 2017 The Authors The Plant Journal © 2017 John Wiley & Sons Ltd.)

Small RNAs, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), play a variety of important regulatory roles in many eukaryotes. Their small size has made it challenging to study them directly in live cells. Here we describe an RNA-based fluorescent sensor for small RNA detection both in vitro and in vivo, adaptable for any small RNA. It utilizes an sxRNA switch for detection of miRNA-mRNA interactions combined with a fluorophore-binding sequence 'Spinach', a GFP-like RNA aptamer for which the RNA-fluorophore complex exhibits strong and consistent fluorescence under an excitation wavelength. Two example sensors, FASTmiR171 and FASTmiR122, can rapidly detect and quantify the levels of miR171 and miR122 in vitro. The sensors can determine relative levels of miRNAs in total RNA extracts with sensitivity similar to small RNA sequencing and northern blots. FASTmiR sensors were also used to estimate the copy number range of miRNAs in total RNA extracts. To localize and analyze the spatial distribution of small RNAs in live, single cells, tandem copies of FASTmiR122 were expressed in different cell lines. FASTmiR122 was able to quantitatively detect the differences in miR122 levels in Huh7 and HEK293T cells demonstrating its potential for tracking miRNA expression and localization in vivo., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

In plants, miR390 directs the production of tasiRNAs from TRANS-ACTING SIRNA3 ( TAS3 ) transcripts to regulate AUXIN RESPONSIVE FACTOR ( ARF ) genes, critical for auxin signaling; these tasiRNAs are known as tasiARFs. To understand the evolution of this miR390- TAS3 - ARF pathway, we characterized homologs of these three genes from thousands of plant species, from bryophytes to angiosperms. We found the lower-stem region of MIR390 genes, critical for accurate DICER-LIKE1 processing, is conserved in sequence in seed plants. We propose a model for the transition of functional tasiRNA sequences in TAS3 genes occurred at the emergence of vascular plants, in which the two miR390 target sites of TAS3 genes showed distinct pairing patterns. Based on the cleavability of miR390 target sites and the distance between target site and tasiARF, we inferred a potential bidirectional processing mechanism exists for some TAS3 genes. We also demonstrated a tight mutual selection between tasiARF and its target genes and that ARGONAUTE7, the partner of miR390, was specified later than other factors in the pathway. All these data illuminate the evolutionary path of the miR390- TAS3 - ARF pathway in land plants and demonstrate the significant variation that occurs in this functionally important and archetypal regulatory circuit., (© 2017 American Society of Plant Biologists. All rights reserved.)

The degradation of small RNAs in plants and animals is associated with small RNA 3' truncation and 3' uridylation and thus relies on exonucleases and nucleotidyl transferases. ARGONAUTE (AGO) proteins associate with small RNAs in vivo and are essential for not only the activities but also the stability of small RNAs. AGO1 is the microRNA (miRNA) effector in Arabidopsis, and its closest homolog, AGO10, maintains stem cell homeostasis in meristems by sequestration of miR165/6, a conserved miRNA acting through AGO1. Here, we show that SMALL RNA DEGRADING NUCLEASES (SDNs) initiate miRNA degradation by acting on AGO1-bound miRNAs to cause their 3' truncation, and the truncated species are uridylated and degraded. We report that AGO10 reduces miR165/6 accumulation by enhancing its degradation by SDN1 and SDN2 in vivo. In vitro, AGO10-bound miR165/6 is more susceptible to SDN1-mediated 3' truncation than AGO1-bound miR165/6. Thus, AGO10 promotes the degradation of miR165/6, which is contrary to the stabilizing effect of AGO1. Our work identifies a class of exonucleases responsible for miRNA 3' truncation in vivo and uncovers a mechanism of specificity determination in miRNA turnover. This work, together with previous studies on AGO10, suggests that spatially regulated miRNA degradation underlies stem cell maintenance in plants.

Dissection of the genetic pathways and mechanisms by which anther development occurs in grasses is crucial for both a basic understanding of plant development and for examining traits of agronomic importance such as male sterility. In rice, MULTIPLE SPOROCYTES1 (MSP1), a leucine-rich-repeat receptor kinase, plays an important role in anther development by limiting the number of sporocytes. OsTDL1a (a TPD1-like gene in rice) encodes a small protein that acts as a cofactor of MSP1 in the same regulatory pathway. In this study, we analyzed small RNA and mRNA changes in different stages of spikelets from wild-type rice, and from msp1 and ostdl1a mutants. Analysis of the small RNA data identified miRNAs demonstrating differential abundances. miR2275 was depleted in the two rice mutants; this miRNA is specifically enriched in anthers and functions to trigger the production of 24-nt phased secondary siRNAs (phasiRNAs) from PHAS loci. We observed that the 24-nt phasiRNAs as well as their precursor PHAS mRNAs were also depleted in the two mutants. An analysis of co-expression identified three Argonaute-encoding genes (OsAGO1d, OsAGO2b, and OsAGO18) that accumulate transcripts coordinately with phasiRNAs, suggesting a functional relationship. By mRNA in situ analysis, we demonstrated a strong correlation between the spatiotemporal pattern of these OsAGO transcripts and phasiRNA accumulations., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Experimental Biology.)

High expression of plant nucleotide binding site leucine-rich repeat (NBS-LRR) defense genes is often lethal to plant cells, a phenotype perhaps associated with fitness costs. Plants implement several mechanisms to control the transcript level of NBS-LRR defense genes. As negative transcriptional regulators, diverse miRNAs target NBS-LRRs in eudicots and gymnosperms. To understand the evolutionary benefits of this miRNA-NBS-LRR regulatory system, we investigated the NBS-LRRs of 70 land plants, coupling this analysis with extensive small RNA data. A tight association between the diversity of NBS-LRRs and miRNAs was found. The miRNAs typically target highly duplicated NBS-LRRs In comparison, families of heterogeneous NBS-LRRs were rarely targeted by miRNAs in Poaceae and Brassicaceae genomes. We observed that duplicated NBS-LRRs from different gene families periodically gave birth to new miRNAs. Most of these newly emerged miRNAs target the same conserved, encoded protein motif of NBS-LRRs, consistent with a model of convergent evolution for these miRNAs. By assessing the interactions between miRNAs and NBS-LRRs, we found nucleotide diversity in the wobble position of the codons in the target site drives the diversification of miRNAs. Taken together, we propose a co-evolutionary model of plant NBS-LRRs and miRNAs hypothesizing how plants balance the benefits and costs of NBS-LRR defense genes., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Plant small RNAs play important roles in transcriptional and posttranscriptional regulation, with ongoing work demonstrating their functions in diverse pathways. Their roles in defense responses are a topic of active investigation, particularly the rich set of micro (mi)RNAs that target disease resistance genes such as nucleotide binding/leucine-rich repeat (NB-LRR) genes. The miRNA-NB-LRR interactions result in the production of phased, secondary small interfering (phasi)RNAs, and phasiRNAs function in both cis and trans to propagate negative regulatory effects across additional members of the target gene family. Yet, while phasiRNAs have the capacity to trigger targeted decay of specific targets, both in cis and trans, their functional relevance in NB-LRR regulation remains largely a matter of speculation.

Summary: We describe miTRATA, the first web-based tool for microRNA Truncation and Tailing Analysis--the analysis of 3' modifications of microRNAs including the loss or gain of nucleotides relative to the canonical sequence. miTRATA is implemented in Python (version 3) and employs parallel processing modules to enhance its scalability when analyzing multiple small RNA (sRNA) sequencing datasets. It utilizes miRBase, currently version 21, as a source of known microRNAs for analysis. miTRATA notifies user(s) via email to download as well as visualize the results online. miTRATA's strengths lie in (i) its biologist-focused web interface, (ii) improved scalability via parallel processing and (iii) its uniqueness as a webtool to perform microRNA truncation and tailing analysis., Availability and Implementation: miTRATA is developed in Python and PHP. It is available as a web-based application from https://wasabi.dbi.udel.edu/∼apps/ta/., Contact: meyers@dbi.udel.edu, Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

Background: Small RNAs (sRNAs) are endogenous sRNAs that play regulatory roles in plant growth, development, and biotic and abiotic stress responses. In plants, one subset of sRNAs, microRNAs (miRNAs) exhibit tissue-differential expression and regulate gene expression mainly through direct cleavage of mRNA or indirectly via production of secondary phased siRNAs (phasiRNAs) that silence cognate target transcripts in trans., Results: Here, we have identified cassava (Manihot esculenta Crantz) miRNAs using high resolution sequencing of sRNA libraries from leaf, stem, callus, male and female flower tissues. To analyze the data, we built a cassava genome database and, via sequence analysis and secondary structure prediction, 38 miRNAs not previously reported in cassava were identified. These new cassava miRNAs included two miRNAs not previously been reported in any plant species. The miRNAs exhibited tissue-differential accumulation as confirmed by quantitative RT-PCR and Northern blot analysis, largely reflecting levels observed in sequencing data. Some of the miRNAs identified were predicted to trigger production of secondary phased siRNAs (phasiRNAs) from 80 PHAS loci., Conclusions: Cassava is a woody perennial shrub, grown principally for its starch-rich storage roots, which are rich in calories. In this study, new miRNAs were identified and their expression was validated using qRT-PCR of RNA from five different tissues. The data obtained expand the list of annotated miRNAs and provide additional new resources for cassava improvement research.

Three soybean [Glycine max (L) Merr.] small RNA libraries were generated and sequenced using the Illumina platform to examine the role of miRNAs during soybean nodulation. The small RNA libraries were derived from root hairs inoculated with Bradyrhizobium japonicum (In_RH) or mock-inoculated with water (Un_RH), as well as from the comparable inoculated stripped root samples (i.e. inoculated roots with the root hairs removed). Sequencing of these libraries identified a total of 114 miRNAs, including 22 novel miRNAs. A comparison of miRNA abundance among the 114 miRNAs identified 66 miRNAs that were differentially expressed between root hairs and stripped roots, and 48 miRNAs that were differentially regulated in infected root hairs in response to B. japonicum when compared to uninfected root hairs (P ≤ 0.05). A parallel analysis of RNA ends (PARE) library was constructed and sequenced to reveal a total of 405 soybean miRNA targets, with most predicted to encode transcription factors or proteins involved in protein modification, protein degradation and hormone pathways. The roles of gma-miR4416 and gma-miR2606b during nodulation were further analysed. Ectopic expression of these two miRNAs in soybean roots resulted in significant changes in nodule numbers. miRNA target information suggested that gma-miR2606b regulates a Mannosyl-oligosaccharide 1, 2-alpha-mannosidase gene, while gma-miR4416 regulates the expression of a rhizobium-induced peroxidase 1 (RIP1)-like peroxidase gene, GmRIP1, during nodulation., (© 2015 Society for Experimental Biology, Association of Applied Biologists and John Wiley & Sons Ltd.)

Small nonprotein-coding microRNAs (miRNAs) are present in most eukaryotes and are central effectors of RNA silencing-mediated mechanisms for gene expression regulation. In plants, DICER-LIKE1 (DCL1) is the founding member of a highly conserved family of RNase III-like endonucleases that function as core machinery proteins to process hairpin-like precursor transcripts into mature miRNAs, small regulatory RNAs, 21-22 nucleotides in length. Zinc finger nucleases (ZFNs) were used to generate single and double-mutants of putative soybean DCL1 homologs, DCL1a and DCL1b, to confirm their functional role(s) in the soybean miRNA pathway. Neither DCL1 single mutant, dcl1a or dcl1b plants, exhibited a pronounced morphological or molecular phenotype. However, the dcl1a/dcl1b double mutant expressed a strong morphological phenotype, characterized by reduced seed size and aborted seedling development, in addition to defective miRNA precursor transcript processing efficiency and deregulated miRNA target gene expression. Together, these findings indicate that the two soybean DCL1 paralogs, DCL1a and DCL1b, largely play functionally redundant roles in the miRNA pathway and are essential for normal plant development., (Copyright © 2016 Curtin et al.)

In eudicot plants, the miR482/miR2118 superfamily regulates and instigates the production of phased secondary small interfering RNAs (siRNAs) from NB-LRR (nucleotide binding leucine-rich repeat) genes that encode disease resistance proteins. In grasses, this miRNA family triggers siRNA production specifically in reproductive tissues from long noncoding RNAs. To understand this functional divergence, we examined the small RNA population in the ancient gymnosperm Norway spruce (Picea abies). As many as 41 miRNA families in spruce were found to trigger phasiRNA (phased, secondary siRNAs) production from diverse PHAS loci, with a remarkable 19 miRNA families capable of targeting over 750 NB-LRR genes to generate phasiRNAs. miR482/miR2118, encoded in spruce by at least 24 precursor loci, targets not only NB-LRR genes to trigger phasiRNA production (as in eudicots) but also noncoding PHAS loci, generating phasiRNAs preferentially in male or female cones, reminiscent of its role in the grasses. These data suggest a dual function of miR482/miR2118 present in gymnosperms that was selectively yet divergently retained in flowering plants. A few MIR482/MIR2118 precursors possess an extremely long stem-loop structure, one arm of which shows significant sequence similarity to spruce NB-LRR genes, suggestive of an evolutionary origin from NB-LRR genes through gene duplication. We also characterized an expanded miR390-TAS3 (TRANS-ACTING SIRNA GENE 3)-ARF (AUXIN RESPONSIVE FACTOR) pathway, comprising 18 TAS3 genes of diverse features. Finally, we annotated spruce miRNAs and their targets. Taken together, these data expand our understanding of phasiRNA network in plants and the evolution of plant miRNAs, particularly miR482/miR2118 and its functional diversification., (© The Author 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The wild strawberry (Fragaria vesca) has recently emerged as an excellent model for cultivated strawberry (Fragaria × ananassa) as well as other Rosaceae fruit crops due to its short seed-to-fruit cycle, diploidy, and sequenced genome. Deep sequencing and parallel analysis of RNA ends were used to identify F. vesca microRNAs (miRNAs) and their target genes, respectively. Thirty-eight novel and 31 known miRNAs were identified. Many known miRNAs targeted not only conserved mRNA targets but also developed new target genes in F. vesca. Significantly, two new clusters of miRNAs were found to collectively target 94 F-BOX (FBX) genes. One of the miRNAs in the new cluster is 22 nucleotides and triggers phased small interfering RNA production from six FBX genes, which amplifies the silencing to additional FBX genes. Comparative genomics revealed that the main novel miRNA cluster evolved from duplications of FBX genes. Finally, conserved trans-acting siRNA pathways were characterized and confirmed with distinct features. Our work identified novel miRNA-FBX networks in F. vesca and shed light on the evolution of miRNAs/phased small interfering RNA networks that regulate large gene families in higher plants., (© 2015 American Society of Plant Biologists. All Rights Reserved.)

Plant microRNAs (miRNAs) play important regulatory roles in a number of developmental processes. The present work investigated the roles of miRNAs during nodule development in the crop legume soybean (Glycine max). Fifteen soybean small RNA libraries were sequenced from different stages of nodule development, including young nodules, mature nodules and senescent nodules. In order to identify the regulatory targets of the miRNAs, five parallel analysis of RNA ends (PARE) libraries were also sequenced from the same stages of nodule development. Sequencing identified 284 miRNAs, including 178 novel soybean miRNAs. Analysis of miRNA abundance identified 139 miRNAs whose expression was significantly regulated during nodule development, including 12 miRNAs whose expression changed > 10-fold. Analysis of the PARE libraries identified 533 miRNA targets, including three nodulation-related genes and eight nodule-specific genes. miR393j-3p was selected for detailed analysis as its expression was significantly regulated during nodule formation, and it targeted a nodulin gene, Early Nodulin 93 (ENOD93). Strong, ectopic expression of miR393j-3p, as well as RNAi silencing of ENOD93 expression, significantly reduced nodule formation. The data indicate that miR393j-3p regulation of ENOD93 mRNA abundance is a key control point for soybean nodule formation., (© 2015 The Authors. New Phytologist © 2015 New Phytologist Trust.)

Small RNAs are a class of non-coding RNAs that are of great importance in gene expression regulatory networks. Different families of small RNAs are generated via distinct biogenesis pathways. One such family specific to plants is that of phased, secondary siRNAs (phasiRNAs); these require RDR6, DCL4, and (typically) a microRNA (miRNA) trigger for their biogenesis. Protein-encoding genes are an important source of phasi-RNAs. The model legume Medicago truncatula generates phasiRNAs from many PHAS loci, and we aimed to investigate their biogenesis and mechanism by which miRNAs trigger these molecules. We modulated miRNA abundances in transgenic tissues showing that the abundance of phasiRNAs correlates with the levels of both miRNA triggers and the target, precursor transcripts. We identified sets of phasiRNAs or PHAS loci that predominantly and substantially increase in response to miRNA overexpression. In the process of validating targets from miRNA overexpression tissues, we found that in the miRNA-mRNA target pairing, the 3' terminal nucleotide (the 22nd position), but not the 10th position, is important for phasiRNA production. Mutating the single 3' terminal nucleotide dramatically diminishes phasiRNA production. Ectopic expression of Medicago NB-LRR-targeting miRNAs in Arabidopsis showed that only a few NB-LRRs are capable of phasiRNA production; our data indicate that this might be due to target inaccessibility determined by sequences flanking target sites. Our results suggest that target accessibility is an important component in miRNA-target interactions that could be utilized in target prediction, and the evolution of mRNA sequences flanking miRNA-target sites may be impacted., (© 2015 The Authors The Plant Journal © 2015 John Wiley & Sons Ltd.)

Background: miRNAs are the most abundant class of small non-coding RNAs, and they are involved in post-transcriptional regulations, playing a crucial role in the refinement of genetic programming during plant development. Here we present a comprehensive picture of miRNA regulation in Vitis vinifera L. plant during its complete life cycle. Furthering our knowledge about the post-transcriptional regulation of plant development is fundamental to understand the biology of such an important crop., Results: We analyzed 70 small RNA libraries, prepared from berries, inflorescences, tendrils, buds, carpels, stamens and other samples at different developmental stages. One-hundred and ten known and 175 novel miRNAs have been identified and a wide grapevine expression atlas has been described. The distribution of miRNA abundance reveals that 22 novel miRNAs are specific to stamen, and two of them are, interestingly, involved in ethylene biosynthesis, while only few miRNAs are highly specific to other organs. Thirty-eight miRNAs are present in all our samples, suggesting a role in key regulatory circuit. On the basis of miRNAs abundance and distribution across samples and on the estimated correlation, we suggest that miRNA expression define organ identity. We performed target prediction analysis and focused on miRNA expression analysis in berries and inflorescence during their development, providing an initial functional description of the identified miRNAs., Conclusions: Our findings represent a very extensive miRNA expression atlas in grapevine, allowing the definition of how the spatio-temporal distribution of miRNAs defines organ identity. We describe miRNAs abundance in specific tissues not previously described in grapevine and contribute to future targeted functional analyses. Finally, we present a deep characterization of miRNA involvement in berry and inflorescence development, suggesting a role for miRNA-driven hormonal regulation.

3' uridylation is increasingly recognized as a conserved RNA modification process associated with RNA turnover in eukaryotes. 2'-O-methylation on the 3' terminal ribose protects micro(mi)RNAs from 3' truncation and 3' uridylation in Arabidopsis. Previously, we identified HESO1 as the nucleotidyl transferase that uridylates most unmethylated miRNAs in vivo, but substantial 3' tailing of miRNAs still remains in heso1 loss-of-function mutants. In this study, we found that among nine other potential nucleotidyl transferases, UTP:RNA uridylyltransferase 1 (URT1) is the single most predominant nucleotidyl transferase that tails miRNAs. URT1 and HESO1 prefer substrates with different 3' end nucleotides in vitro and act cooperatively to tail different forms of the same miRNAs in vivo. Moreover, both HESO1 and URT1 exhibit nucleotidyl transferase activity on AGO1-bound miRNAs. Although these enzymes are able to add long tails to AGO1-bound miRNAs, the tailed miRNAs remain associated with AGO1. Moreover, tailing of AGO1-bound miRNA165/6 drastically reduces the slicing activity of AGO1-miR165/6, suggesting that tailing reduces miRNA activity. However, monouridylation of miR171a by URT1 endows the miRNA the ability to trigger the biogenesis of secondary siRNAs. Therefore, 3' tailing could affect the activities of miRNAs in addition to leading to miRNA degradation.

The evolutionary dynamics of duplicated protein-encoding genes (PEGs) is well documented. However, the evolutionary patterns and consequences of duplicated MIRNAs and the potential influence on the evolution of their PEG targets are poorly understood. Here, we demonstrate the evolution of plant MIRNAs subsequent to a recent whole-genome duplication. Overall, the retention of MIRNA duplicates was correlated to the retention of adjacent PEG duplicates, and the retained MIRNA duplicates exhibited a higher level of interspecific preservation of orthologs than singletons, suggesting that the retention of MIRNA duplicates is related to their functional constraints and local genomic stability. Nevertheless, duplication status, rather than local genic collinearity, was the primary determinant of levels of nucleotide divergence of MIRNAs. In addition, the retention of duplicated MIRNAs appears to be associated with the retention of their corresponding duplicated PEG targets. Furthermore, we characterized the evolutionary novelty of a legume-specific microRNA (miRNA) family, which resulted from rounds of genomic duplication, and consequent dynamic evolution of its NB-LRR targets, an important gene family with primary roles in plant-pathogen interactions. Together, these observations depict evolutionary patterns and novelty of MIRNAs in the context of genomic duplication and evolutionary interplay between MIRNAs and their PEG targets mediated by miRNAs., (© 2015 American Society of Plant Biologists. All rights reserved.)

Small RNAs are ubiquitous, versatile repressors and include (1) microRNAs (miRNAs), processed from mRNA forming stem-loops; and (2) small interfering RNAs (siRNAs), the latter derived in plants by a process typically requiring an RNA-dependent RNA polymerase. We constructed and analyzed an expression atlas of soybean (Glycine max) small RNAs, identifying over 500 loci generating 21-nucleotide phased siRNAs (phasiRNAs; from PHAS loci), of which 483 overlapped annotated protein-coding genes. Via the integration of miRNAs with parallel analysis of RNA end (PARE) data, 20 miRNA triggers of 127 PHAS loci were detected. The primary class of PHAS loci (208 or 41% of the total) corresponded to NB-LRR genes; some of these small RNAs preferentially accumulate in nodules. Among the PHAS loci, novel representatives of TAS3 and noncanonical phasing patterns were also observed. A noncoding PHAS locus, triggered by miR4392, accumulated preferentially in anthers; the phasiRNAs are predicted to target transposable elements, with their peak abundance during soybean reproductive development. Thus, phasiRNAs show tremendous diversity in dicots. We identified novel miRNAs and assessed the veracity of soybean miRNAs registered in miRBase, substantially improving the soybean miRNA annotation, facilitating an improvement of miRBase annotations and identifying at high stringency novel miRNAs and their targets., (© 2014 American Society of Plant Biologists. All rights reserved.)

Parallel analysis of RNA ends (PARE) is a technique utilizing high-throughput sequencing to profile uncapped, mRNA cleavage or decay products on a genome-wide basis. Tools currently available to validate miRNA targets using PARE data employ only annotated genes, whereas important targets may be found in unannotated genomic regions. To handle such cases and to scale to the growing availability of PARE data and genomes, we developed a new tool, 'sPARTA' (small RNA-PARE target analyzer) that utilizes a built-in, plant-focused target prediction module (aka 'miRferno'). sPARTA not only exhibits an unprecedented gain in speed but also it shows greater predictive power by validating more targets, compared to a popular alternative. In addition, the novel 'seed-free' mode, optimized to find targets irrespective of complementarity in the seed-region, identifies novel intergenic targets. To fully capitalize on the novelty and strengths of sPARTA, we developed a web resource, 'comPARE', for plant miRNA target analysis; this facilitates the systematic identification and analysis of miRNA-target interactions across multiple species, integrated with visualization tools. This collation of high-throughput small RNA and PARE datasets from different genomes further facilitates re-evaluation of existing miRNA annotations, resulting in a 'cleaner' set of microRNAs., (© The Author(s) 2014. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The genus Phytophthora consists of many notorious pathogens of crops and forestry trees. At present, battling Phytophthora diseases is challenging due to a lack of understanding of their pathogenesis. We investigated the role of small RNAs in regulating soybean defense in response to infection by Phytophthora sojae, the second most destructive pathogen of soybean. Small RNAs, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), are universal regulators that repress target gene expression in eukaryotes. We identified known and novel small RNAs that differentially accumulated during P. sojae infection in soybean roots. Among them, miR393 and miR166 were induced by heat-inactivated P. sojae hyphae, indicating that they may be involved in soybean basal defense. Indeed, knocking down the level of mature miR393 led to enhanced susceptibility of soybean to P. sojae; furthermore, the expression of isoflavonoid biosynthetic genes was drastically reduced in miR393 knockdown roots. These data suggest that miR393 promotes soybean defense against P. sojae. In addition to miRNAs, P. sojae infection also resulted in increased accumulation of phased siRNAs (phasiRNAs) that are predominantly generated from canonical resistance genes encoding nucleotide binding-leucine rich repeat proteins and genes encoding pentatricopeptide repeat-containing proteins. This work identifies specific miRNAs and phasiRNAs that regulate defense-associated genes in soybean during Phytophthora infection., (© 2014 The Authors The Plant Journal © 2014 John Wiley & Sons Ltd.)

Small RNAs are pivotal regulators of gene expression that guide transcriptional and post-transcriptional silencing mechanisms in eukaryotes, including plants. Here we report a comprehensive atlas of sRNA and miRNA from 3 species of algae and 31 representative species across vascular plants, including non-model plants. We sequence and quantify sRNAs from 99 different tissues or treatments across species, resulting in a data set of over 132 million distinct sequences. Using miRBase mature sequences as a reference, we identify the miRNA sequences present in these libraries. We apply diverse profiling methods to examine critical sRNA and miRNA features, such as size distribution, tissue-specific regulation and sequence conservation between species, as well as to predict putative new miRNA sequences. We also develop database resources, computational analysis tools and a dedicated website, http://smallrna.udel.edu/. This study provides new insights on plant sRNAs and miRNAs, and a foundation for future studies.

In plants, post-transcriptional gene silencing (PTGS) is mediated by DICER-LIKE 1 (DCL1)-dependent microRNAs (miRNAs), which also trigger 21-nucleotide secondary short interfering RNAs (siRNAs) via RNA-DEPENDENT RNA POLYMERASE 6 (RDR6), DCL4 and ARGONAUTE 1 (AGO1), whereas transcriptional gene silencing (TGS) of transposons is mediated by 24-nucleotide heterochromatic (het)siRNAs, RDR2, DCL3 and AGO4 (ref. 4). Transposons can also give rise to abundant 21-nucleotide 'epigenetically activated' small interfering RNAs (easiRNAs) in DECREASED DNA METHYLATION 1 (ddm1) and DNA METHYLTRANSFERASE 1 (met1) mutants, as well as in the vegetative nucleus of pollen grains and in dedifferentiated plant cell cultures. Here we show that easiRNAs in Arabidopsis thaliana resemble secondary siRNAs, in that thousands of transposon transcripts are specifically targeted by more than 50 miRNAs for cleavage and processing by RDR6. Loss of RDR6, DCL4 or DCL1 in a ddm1 background results in loss of 21-nucleotide easiRNAs and severe infertility, but 24-nucleotide hetsiRNAs are partially restored, supporting an antagonistic relationship between PTGS and TGS. Thus miRNA-directed easiRNA biogenesis is a latent mechanism that specifically targets transposon transcripts, but only when they are epigenetically reactivated during reprogramming of the germ line. This ancient recognition mechanism may have been retained both by transposons to evade long-term heterochromatic silencing and by their hosts for genome defence.

Background: The wild grass Brachypodium distachyon has emerged as a model system for temperate grasses and biofuel plants. However, the global analysis of miRNAs, molecules known to be key for eukaryotic gene regulation, has been limited in B. distachyon to studies examining a few samples or that rely on computational predictions. Similarly an in-depth global analysis of miRNA-mediated target cleavage using parallel analysis of RNA ends (PARE) data is lacking in B. distachyon., Results: B. distachyon small RNAs were cloned and deeply sequenced from 17 libraries that represent different tissues and stresses. Using a computational pipeline, we identified 116 miRNAs including not only conserved miRNAs that have not been reported in B. distachyon, but also non-conserved miRNAs that were not found in other plants. To investigate miRNA-mediated cleavage function, four PARE libraries were constructed from key tissues and sequenced to a total depth of approximately 70 million sequences. The roughly 5 million distinct genome-matched sequences that resulted represent an extensive dataset for analyzing small RNA-guided cleavage events. Analysis of the PARE and miRNA data provided experimental evidence for miRNA-mediated cleavage of 264 sites in predicted miRNA targets. In addition, PARE analysis revealed that differentially expressed miRNAs in the same family guide specific target RNA cleavage in a correspondingly tissue-preferential manner., Conclusions: B. distachyon miRNAs and target RNAs were experimentally identified and analyzed. Knowledge gained from this study should provide insights into the roles of miRNAs and the regulation of their targets in B. distachyon and related plants.

MicroRNAs (miRNAs) are small RNAs that derive from endogenous precursors harboring foldback structures. Plant miRNA precursors are quite variable in their size and shape. Still, the miRNA processing machinery, consisting of DICER-LIKE1 (DCL1) and accessory proteins recognize structural features on the precursors to cleave them at specific places releasing the mature miRNAs. The identification of miRNA processing intermediates in plants has mostly relied on a modified 5' RACE method, designed to detect the 5' end of uncapped RNAs. However, this method is time consuming and is, therefore, only practical for the analysis of a handful miRNAs. Here, we present a modification of this approach in order to perform genome-wide analysis of miRNA processing intermediates. Briefly, a reverse transcription is performed with a mixture of specific primers designed against all known miRNA precursors. miRNA processing intermediates are then specifically amplified to generate a library and subjected to deep sequencing. This method, called SPARE (Specific Parallel Amplification of 5' RNA Ends) allows the identification of processing intermediates for most of the Arabidopsis miRNAs. The results enable the determination of the DCL1 processing direction and the cleavage sites introduced by miRNA processing machinery in the precursors. The SPARE method can be easily adapted to detect miRNA-processing intermediates in other systems., (Copyright © 2013 Elsevier Inc. All rights reserved.)

MicroRNAs (miRNAs) derive from longer precursors with fold-back structures. While animal miRNA precursors have homogenous structures, plant precursors comprise a collection of fold-backs with variable size and shape. Here, we design an approach to systematically analyze miRNA processing intermediates and characterize the biogenesis of most of the evolutionarily conserved miRNAs present in Arabidopsis thaliana. We found that plant miRNAs are processed by four mechanisms, depending on the sequential direction of the processing machinery and the number of cuts required to release the miRNA. Classification of the precursors according to their processing mechanism revealed specific structural determinants for each group. We found that the complexity of the miRNA processing pathways occurs in both ancient and evolutionarily young sequences and that members of the same family can be processed in different ways. We observed that different structural determinants compete for the processing machinery and that alternative miRNAs can be generated from a single precursor. The results provide an explanation for the structural diversity of miRNA precursors in plants and new insights toward the understanding of the biogenesis of small RNAs.

Plant small RNAs are 3' methylated by the methyltransferase HUA1 ENHANCER1 (HEN1). In plant hen1 mutants, 3' modifications of small RNAs, including oligo-uridylation (tailing), are associated with accelerated degradation of microRNAs (miRNAs). By sequencing small RNAs of the wild type and hen1 mutants from Arabidopsis thaliana, rice (Oryza sativa), and maize (Zea mays), we found 3' truncation prior to tailing is widespread in these mutants. Moreover, the patterns of miRNA truncation and tailing differ substantially among miRNA families but are conserved across species. The same patterns are also observable in wild-type libraries from a broad range of species, only at lower abundances. ARGONAUTE (AGO1), even with defective slicer activity, can bind these truncated and tailed variants of miRNAs. An ago1 mutation in hen1 suppressed such 3' modifications, indicating that they occur while miRNAs are in association with AGO1, either during or after RNA-induced silencing complex assembly. Our results showed AGO1-bound miRNAs are actively 3' truncated and tailed, possibly reflecting the activity of cofactors acting in conserved patterns in miRNA degradation.

MicroRNAs (miRNAs) are a class of small RNAs that typically function by guiding the cleavage of target messenger RNAs. They have been shown to play major roles in a variety of plant processes, including development, and responses to pathogens and environmental stresses. To identify new miRNAs and regulation in Arabidopsis (Arabidopsis thaliana), 27 small RNA libraries were constructed and sequenced from various tissues, stresses, and small RNA biogenesis mutants, resulting in 95 million genome-matched sequences. The use of rdr2 to enrich the miRNA population greatly enhanced this analysis and led to the discovery of new miRNAs arising from both known and new precursors, increasing the total number of Arabidopsis miRNAs by about 10%. Parallel Analysis of RNA Ends data provide evidence that the majority guide target cleavage. Many libraries represented novel stress/tissue conditions, such as submergence-stressed flowers, which enabled the identification of new stress regulation of both miRNAs and their targets, all of which were validated in wild-type plants. By combining small RNA expression analysis with ARGONAUTE immunoprecipitation data and global target cleavage data from Parallel Analysis of RNA Ends, a much more complete picture of Arabidopsis miRNAs was obtained. In particular, the discovery of ARGONAUTE loading and target cleavage biases gave important insights into tissue-specific expression patterns, pathogen responses, and the role of sequence variation among closely related miRNA family members that would not be evident without this combinatorial approach.

Non-coding RNAs, especially small RNAs, play important roles in many biological processes. Several small RNA types, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), are well-described in rice (Oryza sativa), although much remains to be learned about their function. Many small RNAs along with their targets have been characterized with deep sequencing technologies. Some special classes of these small RNAs have been found to be unique to rice or within the larger group of grasses. The functional and biological roles of numerous plants small RNAs have been described in detail, including functions as varied as the regulation of tissue development, phase transition, or abiotic and biotic stress resistance. Mutant analysis has proven useful in the genetic identification of components involved in small RNA biogenesis and also in identification of regulatory functions of small RNAs. Although many small RNAs have been identified by deep sequencing in rice, their precise regulatory functions and cell-type specificity are in many cases still unknown., (Copyright © 2013 Elsevier Ltd. All rights reserved.)