Controlled RISC loading efficiency of miR168 defined by miRNA duplex structure adjusts ARGONAUTE1 homeostasis

Micro RNAs (miRNAs) are processed from precursor RNA molecules with precisely defined secondary stem-loop structures. ARGONAUTE1 (AGO1) is the main executor component of miRNA pathway and its expression is controlled via the auto-regulatory feedback loop activity of miR168 in plants. Previously we have shown that AGO1 loading of miR168 is strongly restricted leading to abundant cytoplasmic accumulation of AGO-unbound miR168. Here, we report, that intrinsic RNA secondary structure of MIR168a precursor not only defines the processing of miR168, but also precisely adjusts AGO1 loading efficiency determining the biologically active subset of miR168 pool. Our results show, that modification of miRNA duplex structure of MIR168a precursor fragment or expression from artificial precursors can alter the finely adjusted loading efficiency of miR168. In dcl1-9 mutant where, except for miR168, production of most miRNAs is severely reduced this mechanism ensures the elimination of unloaded AGO1 proteins via enhanced AGO1 loading of miR168. Based on this data, we propose a new competitive loading mechanism model for miR168 action: the miR168 surplus functions as a molecular buffer for controlled AGO1 loading continuously adjusting the amount of AGO1 protein in accordance with the changing size of the cellular miRNA pool.     

AGO1 and AGO2 act redundantly in miR408-mediated Plantacyanin regulation

Background

In Arabidopsis, AGO1 and AGO2 associate with small RNAs that exhibit a Uridine and an Adenosine at their 5′ end, respectively. Because most plant miRNAs have a 5′U, AGO1 plays many essential roles in miRNA-mediated regulation of development and stress responses. In contrast, AGO2 has only been implicated in antibacterial defense in association with miR393*, which has a 5′A. AGO2 also participates in antiviral defense in association with viral siRNAs.

Principal Findings

This study reveals that miR408, which has a 5′A, regulates its target Plantacyanin through either AGO1 or AGO2. Indeed, neither ago1 nor ago2 single mutations abolish miR408-mediated regulation of Plantacyanin. Only an ago1 ago2 double mutant appears compromised in miR408-mediated regulation of Plantacyanin, suggesting that AGO1 and AGO2 have redundant roles in this regulation. Moreover, the nature of the 5′ nucleotide of miR408 does not appear essential for its regulatory role because both a wildtype 5′A-MIR408 and a mutant 5′U-MIR408 gene complement a mir408 mutant.

Conclusions/Significance

These results suggest that miR408 associates with both AGO1 and AGO2 based on criteria that differ from the 5′ end rule, reminiscent of miR390-AGO7 and miR165/166-AGO10 associations, which are not based on the nature of the 5′ nucleotide.     

Diversity,expression and mRNA targeting abilities of Argonaute-targeting miRNAs among selected vascular plants

Background

Micro (mi)RNAs are important regulators of plant development. Across plant lineages, Dicer-like 1 (DCL1) proteins process long ds-like structures to produce micro (mi) RNA duplexes in a stepwise manner. These miRNAs are incorporated into Argonaute (AGO) proteins and influence expression of RNAs that have sequence complementarity with miRNAs. Expression levels of AGOs are greatly regulated by plants in order to minimize unwarranted perturbations using miRNAs to target mRNAs coding for AGOs. AGOs may also have high promoter specificity-sometimes expression of AGO can be limited to just a few cells in a plant. Viral pathogens utilize various means to counter antiviral roles of AGOs including hijacking the host encoded miRNAs to target AGOs. Two host encoded miRNAs namely miR168 and miR403 that target AGOs have been described in the model plant Arabidopsis and such a mechanism is thought to be well conserved across plants because AGO sequences are well conserved.

Results

We show that the interaction between AGO mRNAs and miRNAs is species-specific due to the diversity in sequences of two miRNAs that target AGOs, sequence diversity among corresponding target regions in AGO mRNAs and variable expression levels of these miRNAs among vascular plants. We used miRNA sequences from 68 plant species representing 31 plant families for this analysis. Sequences of miR168 and miR403 are not conserved among plant lineages, but surprisingly they differ drastically in their sequence diversity and expression levels even among closely related plants. Variation in miR168 expression among plants correlates well with secondary structures/length of loop sequences of their precursors.

Conclusions

Our data indicates a complex AGO targeting interaction among plant lineages due to miRNA sequence diversity and sequences of miRNA targeting regions among AGO mRNAs, thus leading to the assumption that the perturbations by viruses that use host miRNAs to target antiviral AGOs can only be species-specific. We also show that rapid evolution and likely loss of expression of miR168 isoforms in tobacco is related to the insertion of MITE-like transposons between miRNA and miRNA* sequences, a possible mechanism showing how miRNAs are lost in few plant lineages even though other close relatives have abundantly expressing miRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1049) contains supplementary material, which is available to authorized users.     

Distinct and Cooperative Activities of HESO1 and URT1 Nucleotidyl Transferases in MicroRNA Turnover in Arabidopsis

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.     

Intertwined pathways for Argonaute-mediated microRNA biogenesis in Drosophila

Although Dicer is essential for general microRNA (miRNA) biogenesis, vertebrate mir-451 is Dicer independent. Instead, its short pre-miRNA hairpin is ‘sliced’ by Ago2, then 3′-resected into mature miRNAs. Here, we show that Drosophila cells and animals generate functional small RNAs from mir-451-type precursors. However, their bulk maturation arrests as Ago-cleaved pre-miRNAs, which mostly associate with the RNAi effector AGO2. Routing of pre-mir-451 hairpins to the miRNA effector AGO1 was inhibited by Dicer-1 and its partner Loqs. Loss of these miRNA factors promoted association of pre-mir-451 with AGO1, which sliced them and permitted maturation into ∼23–26 nt products. The difference was due to the 3′ modification of single-stranded species in AGO2 by Hen1 methyltransferase, whose depletion permitted 3′ trimming of Ago-cleaved pre-miRNAs in AGO2. Surprisingly, Nibbler, a 3′–5′ exoribonuclease that trims ‘long’ mature miRNAs in AGO1, antagonized miR-451 processing. We used an in vitro reconstitution assay to identify a soluble, EDTA-sensitive activity that resects sliced pre-miRNAs in AGO1 complexes. Finally, we use deep sequencing to show that depletion of dicer-1 increases the diversity of small RNAs in AGO1, including some candidate mir-451-like loci. Altogether, we document unexpected aspects of miRNA biogenesis and Ago sorting, and provide insights into maturation of Argonaute-cleaved miRNA substrates.     

A non-canonical plant microRNA target site

Plant microRNAs (miRNAs) typically form near-perfect duplexes with their targets and mediate mRNA cleavage. Here, we describe an unconventional miRNA target of miR398 in Arabidopsis, an mRNA encoding the blue copper-binding protein (BCBP). BCBP mRNA carries an miR398 complementary site in its 5′-untranslated region (UTR) with a bulge of six nucleotides opposite to the 5′ region of the miRNA. Despite the disruption of a target site region thought to be especially critical for function, BCBP mRNAs are cleaved by ARGONAUTE1 between nucleotides 10th and 11th, opposite to the miRNA, like conventional plant target sites. Levels of BCBP mRNAs are inversely correlated to levels of miR398 in mutants lacking the miRNA, or transgenic plants overexpressing it. Introducing two mutations that disrupt the miRNA complementarity around the cleavage site renders the target cleavage-resistant. The BCBP site functions outside of the context of the BCBP mRNA and does not depend on 5′-UTR location. Reducing the bulge does not interfere with miR398-mediated regulation and completely removing it increases the efficiency of the slicing. Analysis of degradome data and target predictions revealed that the miR398-BCBP interaction seems to be rather unique. Nevertheless, our results imply that functional target sites with non-perfect pairings in the 5′ region of an ancient conserved miRNA exist in plants.     

Thermodynamic analysis of 5' and 3' single- and 3' double-nucleotide overhangs neighboring wobble terminal base pairs

Thermodynamic parameters are reported for duplex formation of 40 self-complementary RNA duplexes containing wobble terminal base pairs with all possible 3′ single and double-nucleotide overhangs, mimicking the structures of short interfering RNAs (siRNA) and microRNAs (miRNA). Based on nearest neighbor analysis, the addition of a single 3′ dangling nucleotide increases the stability of duplex formation up to 1 kcal/mol in a sequence-dependent manner. The addition of a second dangling nucleotide increases the stability of duplexes closed with wobble base pairs in an idiosyncratic manner. The results allow for the development of a nearest neighbor model, which improves the predication of free energy and melting temperature for duplexes closed by wobble base pairs with 3′ single or double-nucleotide overhangs. Phylogenetic analysis of naturally occurring miRNAs was performed. Selection of the effector miR strand of the mature miRNA duplex appears to be dependent on the orientation of the GU closing base pair rather than the identity of the 3′ double-nucleotide overhang. Thermodynamic parameters for the 5′ single terminal overhangs adjacent to wobble closing base pairs are also presented.     

Arabidopsis Argonaute MID domains use their nucleotide specificity loop to sort small RNAs

The 5'-nucleotide of small RNAs associates directly with the MID domain of Argonaute (AGO) proteins. In humans, the identity of the 5'-base is sensed by the MID domain nucleotide specificity loop and regulates the integrity of miRNAs. In Arabidopsis thaliana, the 5'-nucleotide also controls sorting of small RNAs into the appropriate member of the AGO family; however, the structural basis for this mechanism is unknown. Here, we present crystal structures of the MID domain from three Arabidopsis AGOs, AtAGO1, AtAGO2 and AtAGO5, and characterize their interactions with nucleoside monophosphates (NMPs). In AtAGOs, the nucleotide specificity loop also senses the identity of the 5'-nucleotide but uses more diverse modes of recognition owing to the greater complexity of small RNAs found in plants. Binding analyses of these interactions reveal a strong correlation between their affinities and evolutionary conservation.     

HEN1 recognizes 21-24 nt small RNA duplexes and deposits a methyl group onto the 2' OH of the 3' terminal nucleotide

microRNAs (miRNAs) and small interfering RNAs (siRNAs) in plants bear a methyl group on the ribose of the 3′ terminal nucleotide. We showed previously that the methylation of miRNAs and siRNAs requires the protein HEN1 in vivo and that purified HEN1 protein methylates miRNA/miRNA* duplexes in vitro. In this study, we show that HEN1 methylates both miRNA/miRNA* and siRNA/siRNA* duplexes in vitro with a preference for 21–24 nt RNA duplexes with 2 nt overhangs. We also demonstrate that HEN1 deposits the methyl group on to the 2′ OH of the 3′ terminal nucleotide. Among various modifications that can occur on the ribose of the terminal nucleotide, such as 2′-deoxy, 3′-deoxy, 2′-O-methyl and 3′-O-methyl, only 2′-O-methyl on a small RNA inhibits the activity of yeast poly(A) polymerase (PAP). These findings indicate that HEN1 specifically methylates miRNAs and siRNAs and implicate the importance of the 2′-O-methyl group in the biology of RNA silencing.     

Sorting of Drosophila small silencing RNAs partitions microRNA* strands into the RNA interference pathway

In flies, small silencing RNAs are sorted between Argonaute1 (Ago1), the central protein component of the microRNA (miRNA) pathway, and Argonaute2 (Ago2), which mediates RNA interference. Extensive double-stranded character—as is found in small interfering RNAs (siRNAs)—directs duplexes into Ago2, whereas central mismatches, like those found in miRNA/miRNA* duplexes, direct duplexes into Ago1. Central to this sorting decision is the affinity of the small RNA duplex for the Dcr-2/R2D2 heterodimer, which loads small RNAs into Ago2. Here, we show that while most Drosophila miRNAs are bound to Ago1, miRNA* strands accumulate bound to Ago2. Like siRNA loading, efficient loading of miRNA* strands in Ago2 favors duplexes with a paired central region and requires both Dcr-2 and R2D2. Those miRNA and miRNA* sequences bound to Ago2, like siRNAs diced in vivo from long double-stranded RNA, typically begin with cytidine, whereas Ago1-bound miRNA and miRNA* disproportionately begin with uridine. Consequently, some pre-miRNA generate two or more isoforms from the same side of the stem that differentially partition between Ago1 and Ago2. Our findings provide the first genome-wide test for the idea that Drosophila small RNAs are sorted between Ago1 and Ago2 according to their duplex structure and the identity of their first nucleotide.     

Differential effects of viral silencing suppressors on siRNA and miRNA loading support the existence of two distinct cellular pools of ARGONAUTE1

Plant viruses encode RNA silencing suppressors (VSRs) to counteract the antiviral RNA silencing response. Based on in-vitro studies, several VSRs were proposed to suppress silencing through direct binding of short-interfering RNAs (siRNAs). Because their expression also frequently hinders endogenous miRNA-mediated regulation and stabilizes labile miRNA* strands, VSRs have been assumed to prevent both siRNA and miRNA loading into their common effector protein, AGO1, through sequestration of small RNA (sRNA) duplexes in vivo. These assumptions, however, have not been formally tested experimentally. Here, we present a systematic in planta analysis comparing the effects of four distinct VSRs in Arabidopsis. While all of the VSRs tested compromised loading of siRNAs into AGO1, only P19 was found to concurrently prevent miRNA loading, consistent with a VSR strategy primarily based on sRNA sequestration. By contrast, we provide multiple lines of evidence that the action of the other VSRs tested is unlikely to entail siRNA sequestration, indicating that in-vitro binding assays and in-vivo miRNA* stabilization are not reliable indicator of VSR action. The contrasted effects of VSRs on siRNA versus miRNA loading into AGO1 also imply the existence of two distinct pools of cellular AGO1 that are specifically loaded by each class of sRNAs. These findings have important implications for our current understanding of RNA silencing and of its suppression in plants.     

Genome-Wide Identification of Reverse Complementary microRNA Genes in Plants

MicroRNAs (miRNAs) are ∼21-nucleotide small RNAs (sRNAs) with essential regulatory roles in plants. They are generated from stem-loop-structured precursors through two sequential Dicer-like 1 (DCL1)-mediated cleavages. To date, hundreds of plant miRNAs have been uncovered. However, the question, whether the sequences reverse complementary (RC) to the miRNA precursors could form hairpin-like structures and produce sRNA duplexes similar to the miRNA/miRNA* pairs has not been solved yet. Here, we interrogated this possibility in 16 plant species based on sRNA high-throughput sequencing data and secondary structure prediction. A total of 59 RC sequences with great potential to form stem-loop structures and generate miRNA/miRNA*-like duplexes were identified in ten plants, which were named as RC-miRNA precursors. Unlike the canonical miRNAs, only a few cleavage targets of the RC-miRNAs were identified in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), and none in Soybean (Glycine max) based on degradome data. Surprisingly, the genomic regions surrounding some of the RC-miRNA target recognition sites were observed to be specifically methylated in both Arabidopsis and rice. Taken together, we reported a new class of miRNAs, called RC-miRNAs, which were generated from the antisense strands of the miRNA precursors. Based on the results, we speculated that the mature RC-miRNAs might have subtle regulatory activity through target cleavages, but might possess short interfering RNA-like activity by guiding sequence-specific DNA methylation.     

Transposable element-associated microRNA hairpins produce 21-nt sRNAs integrated into typical microRNA pathways in rice

microRNAs (miRNAs) are a class of small RNAs (sRNAs) of ～21 nucleotides (nt) in length processed from foldback hairpins by DICER-LIKE1 (DCL1) or DCL4. They regulate the expression of target mRNAs by base pairing through RNA-induced silencing complex (RISC). In the RISC, ARGONAUTE1 (AGO1) is the key protein that cleaves miRNA targets at position ten of a miRNA:target duplex. The authenticity of many annotated rice miRNA hairpins is under debate because of their homology to repeat sequences. Some of them, like miR1884b, have been removed from the current release of miRBase based on incomplete information. In this study, we investigated the association of transposable element (TE)-derived miRNAs with typical miRNA pathways (DCL1/4- and AGO1-dependent) using publicly available deep sequencing datasets. Seven miRNA hairpins with 13 unique sRNAs were specifically enriched in AGO1 immunoprecipitation samples and relatively reduced in DCL1/4 knockdown genotypes. Interestingly, these species are ～21-nt long, instead of 24-nt as annotated in miRBase and the literature. Their expression profiles meet current criteria for functional annotation of miRNAs. In addition, diagnostic cleavage tags were found in degradome datasets for predicted target mRNAs. Most of these miRNA hairpins share significant homology with miniature inverted-repeat transposable elements, one type of abundant DNA transposons in rice. Finally, the root-specific production of a 24-nt miRNA-like sRNA was confirmed by RNA blot for a novel EST that maps to the 3′-UTR of a candidate pseudogene showing extensive sequence homology to miR1884b hairpin. Our data are consistent with the hypothesis that TEs can serve as a driving force for the evolution of some MIRNAs, where co-opting of DICER-LIKE1/4 processing and integration into AGO1 could exapt transcribed TE-associated hairpins into typical miRNA pathways.