Email updates

Keep up to date with the latest news and content from Silence and BioMed Central.

This article is part of the series Non-coding RNA.

Open Access Research

A 5′-uridine amplifies miRNA/miRNA* asymmetry in Drosophila by promoting RNA-induced silencing complex formation

Hervé Seitz12*, Jogender S Tushir3 and Phillip D Zamore3*

Author Affiliations

1 Laboratoire de Biologie Moléculaire Eucaryote, 118 route de Narbonne, Université Toulouse III Paul Sabatier (UPS), F-31000 Toulouse, France

2 Laboratoire de Biologie Moléculaire Eucaryote, 118 route de Narbonne, Centre national de la recherche scientifique (CNRS), F-31000 Toulouse, France

3 Howard Hughes Medical Institute and Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, 364 Plantation street, Worcester, MA 01605, USA

For all author emails, please log on.

Silence 2011, 2:4  doi:10.1186/1758-907X-2-4


The electronic version of this article is the complete one and can be found online at: http://www.silencejournal.com/content/2/1/4


Received:26 January 2011
Accepted:7 June 2011
Published:7 June 2011

© 2011 Seitz et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

MicroRNA (miRNA) are diverse in sequence and have a single known sequence bias: they tend to start with uridine (U).

Results

Our analyses of fly, worm and mouse miRNA sequence data reveal that the 5′-U is recognized after miRNA production. Only one of the two strands can be assembled into Argonaute protein from a single miRNA/miRNA* molecule: in fly embryo lysate, a 5′-U promotes miRNA loading while decreasing the loading of the miRNA*.

Conclusion

We suggest that recognition of the 5′-U enhances Argonaute loading by a mechanism distinct from its contribution to weakening base pairing at the 5′-end of the prospective miRNA and, as recently proposed in Arabidopsis and in humans, that it improves miRNA precision by excluding incorrectly processed molecules bearing other 5′-nt.

Background

MicroRNA (miRNA) are approximately 22-nt regulatory RNA that direct members of the Argonaute protein family to their mRNA targets [1]. Together, miRNA guide and the Argonaute protein form the core of the RNA-induced silencing complex (RISC), which recognizes its mRNA targets primarily through its seed sequence, nt 2 through nt 7 [2].

The RNase III enzymes Drosha and Dicer excise most animal miRNA from long primary transcripts (pri-miRNA). Drosha cleaves pri-miRNA to release an approximately 65-nt pre-miRNA; Dicer cleaves the pre-miRNA to liberate a miRNA/miRNA* duplex. The duplex is then loaded into an Argonaute protein. The geometry of the miRNA/miRNA* duplex during the loading reaction determines the fate of each small RNA: the miRNA binds tightly to Argonaute, with its 5′-nt anchored in a positively charged pocket in the Mid domain of the protein [3,4]. The miRNA* assumes the same position as subsequent mRNA targets and is held to the complex predominantly by seed sequence base pairing. A seed sequence mismatch between the miRNA and its miRNA* is believed to promote miRNA* dissociation [5,6]. A subset of Argonaute proteins can cleave the miRNA* if it is extensively paired to the miRNA, triggering its destruction [7-10]. The orientation of the duplex during Argonaute loading is not random: the miRNA is usually the strand with the less stably paired 5′-end in the duplex [11,12]. Consequently, the duplex liberated by Dicer determines the identity of the miRNA.

miRNA sequences are diverse, and only one common sequence motif has been identified. Most miRNA begin with a 5′-uridine (5′-U). In plants, a 5′-U directs miRNA to AGO1, small RNA that begin with adenosine (A) load AGO2 and those that start with cytidine (C) load AGO5 [13-15]. Likewise, the 5′-nt of fly small RNA participates in sorting, with a 5′-U directing small RNA toward Ago1 and a 5′-C favoring Ago2 [16-19]. In mammals, the Mid domain of Ago2, the homolog of Drosophila Ago1, specifically recognizes a 5′-U or 5′-A [20], explaining why miRNA tend to start with those nucleotides, but fly and worm miRNA typically begin with 5′-U but not 5′-A. Moreover, small RNA sorting in flies and worms also reflects the secondary structure of the miRNA/miRNA* duplex, with centrally paired duplexes preferentially loaded into one Argonaute, - Ago2 in flies and RDE-1 in worms, - and duplexes bearing a central mismatch directed toward the major miRNA-binding Argonautes, - Ago1 in flies and the paralogous ALG-1/ALG-2 proteins in worms [5,6,17-19,21-23].

We investigated the function of 5′-U in animal miRNA. Our statistical analyses of sequencing data from flies, worms and mice reveal that 5′-U is recognized after miRNA/miRNA* production by Dicer cleavage of the pre-miRNA. Our experimental results show that 5′-U facilitates loading of miRNA while decreasing loading of miRNA*, consistent with the view that only one of the two strands can be assembled from a single miRNA/miRNA* molecule. Our data support the view that 5′-U enhances RISC assembly by a mechanism distinct from its contribution to destabilizing base pairing at the 5′-end of miRNA. Similarly to what has been proposed in Arabidopsis thaliana and in Homo sapiens [13,20], our data also suggest that recognition of the first miRNA nucleotide during loading may select against incorrectly processed molecules bearing 5′-nt other than 5′-U.

Results and discussion

5′-U acts after miRNA processing

We used high-throughput sequencing data to examine the 5′-sequence bias of miRNA and miRNA*. miRNA are far more likely to begin with U in flies (P value <10-15), worms (P value <10-15) or mice (P value = 1.1 × 10-14) than would be expected from their general nucleotide composition (Figure 1, Additional file 1, Figure S1, and Additional file 2, Figure S2). Conversely, miRNA* were less likely than expected to begin with U in flies (P value = 0.0029), worms (P value = 0.017) or mice (P value = 0.0020).

thumbnailFigure 1. Fly miRNA tend to start with U. Each miRNA or miRNA* isoform derived from a common pre-miRNA was weighted according to its abundance in the pooled deep-sequencing libraries, and the sequence composition analyses for all small RNA from different pre-miRNA that were read at least 100 times in the pooled libraries were weighted equally. Gray, nucleotide frequency at position 1; white, 100 sets of nucleotides randomly selected from nt 1-18 of the miRNA and miRNA* species to assess the overall nucleotide composition of miRNA and miRNA*. Each random set had the same size as the corresponding set of miRNA or miRNA* 5′-nt. P values measure the probability of picking a random set from nt 1-18 with the same nucleotide frequency as the actual set of 5′-nt.

Additional file 1. Figure S1. Caenorhabditis elegans miRNA tend to start with a uridine. Gray, nucleotide frequency at position 1; white, nucleotide frequency at random positions in the miRNA or miRNA* sequence (means ± standard deviation (SD)).

Format: PDF Size: 170KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 2. Figure S2. Mouse miRNA tend to start with a uridine. Gray, nucleotide frequency at position 1; white, nucleotide frequency at random positions in the miRNA or miRNA* sequence (means ± SD).

Format: PDF Size: 172KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

In theory, a 5′-U might facilitate Drosha cleavage of the pri-miRNA or pre-miRNA export from the nucleus. Such a role for a 5′-U would be reflected in a greater likelihood of both miRNA and miRNA* derived from the 5′-arm of the pre-miRNA stem to begin with U compared to those residing in the 3′ arm. We compared the approximately 40% of fly, 35% of worm and 50% of mouse miRNA that reside in the 5′-arm of their pre-miRNA to their 3′ counterparts. Our analysis argues against a role for a 5′-U in Drosha processing or nuclear export. miRNA tend to start with a U, regardless of their position in the pre-miRNA (Figure 1, Additional file 1, Figure S1, and Additional file 2, Figure S2). Moreover, miRNA* sequences tend not to begin with U, even when they derive from the pre-miRNA 5′-arm. Our data similarly exclude a role for a 5′-U in cleavage of the pre-miRNA by Dicer, which would favor a 5′-U for miRNA and miRNA* derived from the 3′-arm.

miRNA asymmetry correlates with first nucleotide identity

To test whether 5′-U plays a role in assembling a miRNA into RISC, we separately evaluated the 5′-nt frequencies in flies of highly asymmetric duplexes (miRNA/miRNA* ≥10; 79 duplexes), moderately asymmetric duplexes (2 < miRNA/miRNA* < 10; 33 duplexes) and quasisymmetric duplexes (miRNA/miRNA* < 2; 10 duplexes). If the identity of the 5′-nt affects miRNA loading, then the most asymmetric miRNA should exhibit a higher 5′-U bias than the least asymmetric miRNA. Indeed, the most highly asymmetric miRNA have a higher frequency of 5′-U (79%) than moderately asymmetric miRNA (61%) or quasisymmetric miRNA and miRNA* (32%) (Figure 2), which is in line with the previously published observation that the most asymmetric human miRNA tend to be richer in 5′-U [24]. Moreover, miRNA* strands from highly asymmetric duplexes have a significantly lower frequency of 5′-U (16.5%) than those from moderately asymmetric or quasisymmetric duplexes. In fact, miRNA* strands have a significantly lower frequency of U at their 5′-ends than across their entire sequence, while the frequency of an initial U was indistinguishable from the overall U frequency in miRNA* from moderately asymmetric or quasisymmetric duplexes.

thumbnailFigure 2. Fly miRNA asymmetry correlates with the identity of the first nucleotide of the small RNA. miRNA/miRNA* duplexes were binned according to their asymmetry: highly asymmetric, miRNA/miRNA* ≥10 in the pooled deep-sequencing libraries; moderately asymmetric, 10 > miRNA/miRNA* ≥ 2; quasisymmetric: miRNA/miRNA* <2; and analyzed as in Figure 1.

Strikingly, the most asymmetric miRNA also exhibit a lower than expected frequency of 5′-A (Figure 2, top left), whereas the thermodynamic stability rule would have predicted a high frequency of both U and A. This observation suggests that 5′-nt identity, not just thermodynamic asymmetry, contributes to the differential loading of miRNA and miRNA* in vivo.

Initial nucleotide identity influences miRNA loading in vitro

Several studies have proposed that a U at the 5′-end of a small RNA directly promotes its loading into Ago1 in flies [18,19,25,26]. We measured the effect of initial nucleotide identity on the efficiency of loading of the miR-2a/miR-2a-1* duplex in Drosophila embryo lysate. To avoid altering the thermodynamic stability of the 5′-ends of the duplex, we designed them so that changing the 5′-nt preserved the pattern and strength of base pairing. To measure the association of miR-2a and miR-2a-1* with mature RISC, we assembled RISC in Drosophila embryo lysate using a duplex in which one strand was 5′-32P-radiolabeled, then captured the radiolabeled strand using a complementary 2′-O-methyl oligonucleotide tethered to a magnetic bead (Figure 3). Labeling either the miRNA or the miRNA* strand (always capturing RISC with an oligonucleotide complementary to the labeled strand), we were able to quantify precisely both miRNA and miRNA* loading by scintillation counting. Ultraviolet cross-linking and RISC capture control experiments demonstrated that the amount of radioactivity captured minus the amount recovered when the duplex was incubated with N-ethylmaleimide (NEM)-inactivated lysate reflected the amount of single-stranded miRNA or miRNA* produced by assembly of Ago1 RISC (Additional file 3, Figure S3, and Additional file 4, Figure S4).

thumbnailFigure 3. Identity of the first miRNA nucleotide affects duplex asymmetry. RNA-induced silencing complex (RISC) loading of miRNA and miRNA* strands was measured after in vitro assembly. Data are reported as means ± standard deviation for three independent replicate experiments. (A) Swapping the terminal uridine:adenosine (U:A) pair of the miR-2a/miR-2a-1* duplex decreased miRNA loading and increased miRNA* loading. (B) The effect of the terminal A:U pair mainly reflects the identity of the first miRNA nucleotide, which affects both miRNA and miRNA* loading (C), whereas the identity of the facing miRNA* nucleotide has no detectable effect (D).

Additional file 3. Figure S3. Faithful in vitro reconstitution of miRNA loading. (A and B) miR-2a, let-7, miR-2c and miR-184 are correctly loaded into Ago1 in fly embryo lysate, and a let-7/anti-let-7 small interfering RNA is correctly loaded into Ago2. Left: lysate prepared from embryos from dcr-2L811fsX- and ago2414-mutant mothers; Ago1-depleted and HA-depleted, wild-type embryo lysate-immunodepleted using anti-Ago1 or anti-HA (hemagglutinin epitope) antibody. Lysate was incubated for one hour with 5′-32P-radiolabeled Drosophila melanogaster miR-2a paired with 5′ phosphorylated miR-2a-1*. Single-stranded (ssRNA), 5′-32P-radiolabeled miR-2a was incubated for one hour in embryo lysate. Each sample was cross-linked using 254 nm ultraviolet light. (C and D) miRNA/miRNA* asymmetry is recapitulated in vitro. (E) In S2 cells, pre-miR-2a-1 liberates both miRNA and miRNA*, and both strands are efficiently loaded into RNA-induced silencing complex as we observed in vitro (Figure 3).

Format: PDF Size: 1.6MB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 4. Figure S4. The Ago2 loading machinery has a moderate effect on miR-184 loading preferences, while it strongly affects miR-184* loading preferences. Left: miR-184 and miR-184* capture assay in dcr-2L811fsX-null mutant embryo lysate. Right: miR-184 and miR-184* capture assay in wild-type lysate.

Format: PDF Size: 159KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Both authentic miR-2a and miR-2a-1* begin with U; the 5′-U of miR-2a is paired to A19 of miR-2a-1*. Inverting this U:A base pair so that miR-2a began with A nearly halved the amount of miRNA assembled into RISC and more than doubled the amount of miR-2a-1* (Figure 3A). Thus, a change in the identity of the first nucleotide of the miRNA decreased the efficiency of assembly of the miRNA into RISC and increased assembly of the miRNA* while preserving the relative thermodynamic asymmetry of the duplex.

When the initial U:A base pair of miR-2a/miR-2a-1* was altered, UU assembled more miRNA into RISC than did AA (Figure 3B). Notably, an AA mismatch at the 5′-end of the miRNA more than doubled the amount of miRNA* incorporated into RISC. Next, we examined a series of miR-2a/miR-2a* derivatives in which the 19th base of miR-2a* was always C, ensuring that duplex stability was the same when the miRNA began with U or A. Again, a 5′-U favored miRNA loading and disfavored miRNA* loading (Figure 3C). When the 5′-U was replaced with inosine, which can pair to the miRNA* C at position 19, only slightly less miRNA was assembled into RISC than that observed for an A/C mismatch. We conclude that the identity of the first miRNA nucleotide contributes more to the loading of miR-2a than do differences in the stability of the duplex termini. Reciprocally, when the first miRNA nucleotide was C, the identity of miRNA* nt 19 did not have any significant effect on miRNA or miRNA* loading (Figure 3D), demonstrating that the effect shown in Figure 3A reflects a mutation of the first miRNA nucleotide, not the change in miRNA* nt 19. Experiments using miR-14 and miR-184 gave similar results (Additional file 5, Figure S5).

Additional file 5. Figure S5. Identity of miRNA nucleotide 1 affects duplex asymmetry. (A) Changing the 5′-uridine (5′-U) of miR-14 to 5′-adenosine (5′-A) decreased miRNA loading. (B) Changing the 5′-U of miR-184 into 5′-cytidine (5′-C) decreased miRNA loading; mutating it to ribothymidine increased miRNA loading. Changing the 5′-nt of miR-2a (C) or miR-184 (D) into 5′-guanidine (5′-G) decreases miRNA loading (relatively to a 5′-A).

Format: PDF Size: 195KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Strikingly, the order of preference for nt 1 was not the same across the three tested miRNA: miR-2a preferred U > A > C (Figure 3), miR-14 preferred U ~ C > A and miR-184 preferred U ~ A > C (Additional file 6, Figure S6). Hence additional features in the miRNA/miRNA* duplex must influence the order of preference for miRNA nt 1. Mutating the overhanging nucleotide in miR-184* did not alter the efficiency of loading miR-184 (Additional file 7, Figure S7), excluding a role for base pairing between nt 1 and the 3′ overhang of the miRNA*.

Additional file 6. Figure S6. Duplex-specific order of preference on the identity of the first nucleotide. (A) Regardless of the identity of the facing (p19*) nucleotide, miR-2a is better loaded if it starts with a U than if it starts with an A than if it starts with a C (U > A > C). (B) miR-14 is better loaded if it starts with a U or a C than if it starts with an A (U ~ C > A). (C) miR-184 is better loaded if it starts with a ribothymidine than if it starts with a U or an A than if it starts with a C (T > U ~ A > C).

Format: PDF Size: 157KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 7. Figure S7. The sequence of the miRNA* 3′ overhang is not responsible for miRNA-specific preferences for nt 1. The modest effect of the identity of the miRNA* 3′-most nucleotide does not correlate with the base-pairing ability of miRNA nt 1 to the miRNA* 3′ terminus.

Format: PDF Size: 125KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Covarying features in miRNA/miRNA* duplexes suggest that the identity of nt 2 affects the order of preference for miRNA nt 1

If a sequence or structural feature affects the order of preference for nt 1, then these two features should evolve together. We searched for significant covariation between nt 1 identity and other sequence or structural motifs in miRNA/miRNA* duplexes. For Drosophila miRNA/miRNA*, the identity of miRNA nt 1 covaries with the identity of the facing nucleotide on the miRNA* strand, the identity of the second nucleotide of the miRNA strand and the base-pairing status of the 15th nucleotide of the miRNA strand (Figure 4A). Mutating miRNA nt 2 in miR-2a and miR-184 influenced the order of preference for nt 1 in flies (Figures 4B and 4C).

thumbnailFigure 4. In Drosophila, the identity of the second miRNA nucleotide influences the effect of the first nucleotide. (A) The y-axis shows the significance (Fisher's exact test) of observed covariation between the identity of miRNA nt 1 and nucleotide identity (black) or base-pairing status (gray) along the miRNA and the miRNA* strand. The low values for position 1 in the black curve occur because the identity of miRNA nt 1 correlates strongly with itself. (B and C) Mutating nt 2 in miR-2a (B) or miR-184 (C) changed the order of nt 1 preference for miRNA loading.

Strikingly, the influence of nt 2 on nt 1 seems to be specific for flies. Neither worm nor mouse miRNA/miRNA* show such covariation (Additional file 8, Figure S8). Caenorhabditis elegans miRNA nt 1 covaries mostly with the base-pairing status of miRNA nt 18 and the identity of the miRNA* nucleotide facing miRNA nt 3. In mouse, nt 1 covaries with the identity of miRNA nt 12 as well as several positions at the 3′ end of the miRNA strand. The sequence composition of miRNA differs greatly between flies and humans [24], suggesting that the nucleotide preference of the miRNA loading machinery has evolved since the divergence of protostomes and deuterostomes, with only the overall tendency for miRNA to start with U remaining conserved.

Additional file 8. Figure S8. Structure and sequence features covarying with the identity of miRNA nt 1 in C. elegans and mouse. See Figure 4A legend for details. (A) Covariation in Mus musculus miRNA/miRNA* duplexes. (B) Covariation in C. elegans miRNA/miRNA* duplexes.

Format: PDF Size: 95KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Conclusions

Our data support the view that a U at the 5′-end of a miRNA favors RISC loading in flies and, given both our informatics data and the broad phylogenetic conservation of the 5′-U bias among miRNA in worms and mice, likely in animals generally.

The Drosophila Ago1 loading machinery remains to be identified, although chaperones have been implicated in assembling miRNA into RISC [6,27,28]. It is tempting to speculate that the requirement for the miRNA 5′-end to be the less thermodynamically stable in a miRNA/miRNA* duplex reflects the need for the first nucleotide to be single-stranded to present it to components of the RISC loading machinery or to Ago1 itself.

Why has the miRNA pathway evolved to prefer a 5′-U? The likely answer is that preferential loading of miRNA starting with U improves the precision of the miRNA 5′-end [13]. Drosha and Dicer generate pools of miRNA/miRNA* duplexes with alternative 5′- and 3′-ends; loading of these duplexes into Drosophila Ago2, -which prefers 5′-C, - has been shown to purify this population of miRNA [29], loading preferentially the miRNA isoforms bearing a 5′-C [19,25]. The preference of the Ago1 loading machinery or of Ago1 itself for 5′-U could similarly restrict entry into the Ago1 pathway by loading only miRNA isoforms that begin with U. Consistent with this idea, the pre-miRNA nucleotides flanking miRNA nt 1 tend to be depleted in U (Additional file 9, Figure S9). Such a purifying selection could ensure that most mature miRNA have the correct 5′-end and therefore the correct seed sequence, ensuring that they regulate the appropriate mRNA targets.

Additional file 9. Figure S9. Pre-miRNA nucleotides flanking miRNA nt 1 are depleted of U. U frequency was measured in pre-miRNA covered by at least 100 reads in the pooled deep-sequencing libraries. The 5′-most nucleotide of mature miRNA is enriched in U (position 0 on the x-axis), while its flanking nucleotides are depleted. The horizontal line indicates the mean U frequency in 100 random sets of nucleotides picked from the corresponding 21-nt segment in the analyzed pre-miRNA. Dashed lines indicate 95% confidence intervals.

Format: PDF Size: 982KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Methods

In vitro reconstitution of miRNA/miRNA* loading 5′ phosphorylated miRNA/miRNA* (approximately 20 nM; the strand measured was 32P-radiolabeled) was incubated with zero- to two-hour fly embryo lysate for one hour at 25°C [30]. Assembly was stopped with NEM [7]. Two-thirds of each assembly reaction were incubated with biotinylated 2′-O-methyl capture oligonucleotide (Table 1) tethered to streptavidin-coated magnetic beads (MyOne Streptavidin C1 DYNAL Magnetic Beads; Invitrogen Corp., Carlsbad, CA, USA) for one hour at 25°C. The radioactivity in the remaining one-third of each reaction was measured by scintillation counting to allow data normalization. Typical replicate-to-replicate variability (standard deviation/mean) was approximately 5%. P values were calculated using Student's t-test assuming equal variances, and distribution normality and homogeneity of variances were assessed using the Shapiro-Wilk test and Levene's test.

Table 1. Synthetic oligonucleotides used in this studya

Covariation analysis

miRNA with ≥100 reads in the pooled deep-sequencing libraries were selected (see Table 2 for the list of analyzed deep-sequencing libraries). The most abundant isoform of each strand was retained. We evaluated the identity and base-pairing status (using RNAcofold, part of the Vienna RNA Secondary Structure Package; available at http://www.tbi.univie.ac.at/RNA/ webcite) of each of the first 18 nt. If the pairing probability of a nucleotide was >0.5, it was called paired. The analysis defined 18 nt identities, starting from either the 5′- or the 3′-end, and 18 base-pairing statuses, starting from either the 5′- or the 3′-end, with a total of 144 features per miRNA/miRNA* duplex. Fisher's exact test was used to evaluate the significance of covariation between these 144 features and the identity of the first miRNA nucleotide using the R Project for Statistical Computing statistical package (http://www.r-project.org/ webcite).

Table 2. High throughput sequencing data used in this studya

Competing interests

PDZ is a cofounder and member of the scientific advisory board of Alnylam Pharmaceuticals, Inc., and a member of the scientific advisory board of Regulus Therapeutics, L.L.C.

Authors' contributions

HS and PDZ planned the experiments and wrote the manuscript. HS and JST performed the experiments. All authors read and approved the final manuscript.

Acknowledgements

We thank Yukihide Tomari, Tomoko Kawamata and Miyuki Mitomi for dcr-2 embryo lysate and sharing unpublished data, as well as members of the Zamore laboratory for critical comments on the manuscript. This work was supported in part by National Institutes of Health grants GM62862 and GM65236 (to PDZ) and European Molecular Biology Organization long-term (ALTF 9102004) and Human Frontier Science Program (LT00575/2005-L) fellowships (to HS).

References

  1. Ghildiyal M, Zamore PD: Small silencing RNAs: an expanding universe.

    Nat Rev Genet 2009, 10:94-108. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Bartel DP: MicroRNAs: target recognition and regulatory functions.

    Cell 2009, 136:215-233. PubMed Abstract | Publisher Full Text OpenURL

  3. Ma JB, Yuan YR, Meister G, Pei Y, Tuschl T, Patel DJ: Structural basis for 5′-end-specific recognition of guide RNA by the A. fulgidus Piwi protein.

    Nature 2005, 434:666-670. PubMed Abstract | Publisher Full Text OpenURL

  4. Parker JS, Roe SM, Barford D: Structural insights into mRNA recognition from a PIWI domain-siRNA guide complex.

    Nature 2005, 434:663-666. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Tomari Y, Du T, Zamore PD: Sorting of Drosophila small silencing RNAs.

    Cell 2007, 130:299-308. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Kawamata T, Seitz H, Tomari Y: Structural determinants of miRNAs for RISC loading and slicer-independent unwinding.

    Nat Struct Mol Biol 2009, 16:953-960. PubMed Abstract | Publisher Full Text OpenURL

  7. Matranga C, Tomari Y, Shin C, Bartel DP, Zamore PD: Passenger-strand cleavage facilitates assembly of siRNA into Ago2-containing RNAi enzyme complexes.

    Cell 2005, 123:607-620. PubMed Abstract | Publisher Full Text OpenURL

  8. Miyoshi K, Tsukumo H, Nagami T, Siomi H, Siomi MC: Slicer function of Drosophila Argonautes and its involvement in RISC formation.

    Genes Dev 2005, 19:2837-2848. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Rand TA, Petersen S, Du F, Wang X: Argonaute2 cleaves the anti-guide strand of siRNA during RISC activation.

    Cell 2005, 123:621-629. PubMed Abstract | Publisher Full Text OpenURL

  10. Leuschner PJ, Ameres SL, Kueng S, Martinez J: Cleavage of the siRNA passenger strand during RISC assembly in human cells.

    EMBO Rep 2006, 7:314-320. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Khvorova A, Reynolds A, Jayasena SD: Functional siRNAs and miRNAs exhibit strand bias.

    Cell 2003, 115:209-216. PubMed Abstract | Publisher Full Text OpenURL

  12. Schwarz DS, Hutvágner G, Du T, Xu Z, Aronin N, Zamore PD: Asymmetry in the assembly of a the RNAi enzyme complex.

    Cell 2003, 115:199-208. PubMed Abstract | Publisher Full Text OpenURL

  13. Mi S, Cai T, Hu Y, Chen Y, Hodges E, Ni F, Wu L, Li S, Zhou H, Long C, Chen S, Hannon GJ, Qi Y: Sorting of small RNAs into Arabidopsis Argonaute complexes is directed by the 5′ terminal nucleotide.

    Cell 2008, 133:116-127. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Montgomery TA, Howell MD, Cuperus JT, Li D, Hansen JE, Alexander AL, Chapman EJ, Fahlgren N, Allen E, Carrington JC: Specificity of ARGONAUTE7-miR390 interaction and dual functionality in TAS3 trans-acting siRNA formation.

    Cell 2008, 133:128-141. PubMed Abstract | Publisher Full Text OpenURL

  15. Takeda A, Iwasaki S, Watanabe T, Utsumi M, Watanabe Y: The mechanism selecting the guide strand from small RNA duplexes is different among Argonaute proteins.

    Plant Cell Physiol 2008, 49:493-500. PubMed Abstract | Publisher Full Text OpenURL

  16. Ghildiyal M, Seitz H, Horwich MD, Li C, Du T, Lee S, Xu J, Kittler EL, Zapp ML, Weng Z, Zamore PD: Endogenous siRNAs derived from transposons and mRNAs in Drosophila somatic cells.

    Science 2008, 320:1077-1081. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Czech B, Zhou R, Erlich Y, Brennecke J, Binari R, Villalta C, Gordon A, Perrimon N, Hannon GJ: Hierarchical rules for Argonaute loading in Drosophila.

    Mol Cell 2009, 36:445-456. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Okamura K, Liu N, Lai EC: Distinct mechanisms for microRNA strand selection by Drosophila Argonautes.

    Mol Cell 2009, 36:431-444. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Ghildiyal M, Xu J, Seitz H, Weng Z, Zamore PD: Sorting of Drosophila small silencing RNAs partitions microRNA* strands into the RNA interference pathway.

    RNA 2010, 16:43-56. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Frank F, Sonenberg N, Nagar B: Structural basis for 5′-nucleotide base-specific recognition of guide RNA by human AGO2.

    Nature 2010, 465:818-822. PubMed Abstract | Publisher Full Text OpenURL

  21. Förstemann K, Horwich MD, Wee L, Tomari Y, Zamore PD: Drosophila microRNAs are sorted into functionally distinct Argonaute complexes after production by Dicer-1.

    Cell 2007, 130:287-297. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Steiner FA, Hoogstrate SW, Okihara KL, Thijssen KL, Ketting RF, Plasterk RH, Sijen T: Structural features of small RNA precursors determine Argonaute loading in Caenorhabditis elegans.

    Nat Struct Mol Biol 2007, 14:927-933. PubMed Abstract | Publisher Full Text OpenURL

  23. Jannot G, Boisvert ME, Banville IH, Simard MJ: Two molecular features contribute to the Argonaute specificity for the microRNA and RNAi pathways in C. elegans.

    RNA 2008, 14:829-835. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Hu HY, Yan Z, Xu Y, Hu H, Menzel C, Zhou YH, Chen W, Khaitovich P: Sequence features associated with microRNA strand selection in humans and flies.

    BMC Genomics 2009, 10:413. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  25. Ameres SL, Hung JH, Xu J, Weng Z, Zamore PD: Target RNA-directed tailing and trimming purifies the sorting of endo-siRNAs between the two Drosophila Argonaute proteins.

    RNA 2011, 17:54-63. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Czech B, Hannon GJ: Small RNA sorting: matchmaking for Argonautes.

    Nat Rev Genet 2011, 12:19-31. PubMed Abstract | Publisher Full Text OpenURL

  27. Iwasaki S, Kobayashi M, Yoda M, Sakaguchi Y, Katsuma S, Suzuki T, Tomari Y: Hsc70/Hsp90 chaperone machinery mediates ATP-dependent RISC loading of small RNA duplexes.

    Mol Cell 2010, 39:292-299. PubMed Abstract | Publisher Full Text OpenURL

  28. Johnston M, Geoffroy MC, Sobala A, Hay R, Hutvágner G: HSP90 protein stabilizes unloaded Argonaute complexes and microscopic P-bodies in human cells.

    Mol Biol Cell 2010, 21:1462-1469. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Seitz H, Ghildiyal M, Zamore PD: Argonaute loading improves the 5′ precision of both microRNAs and their miRNA* strands in flies.

    Curr Biol 2008, 18:147-151. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Haley B, Tang G, Zamore PD: In vitro analysis of RNA interference in Drosophila melanogaster.

    Methods 2003, 30:330-336. PubMed Abstract | Publisher Full Text OpenURL