Processing of Drosophila endo-siRNAs depends on a specific Loquacious isoform

Drosophila melanogaster expresses three classes of small RNAs, which are classified according to their mechanisms of biogenesis. MicroRNAs are ∼22–23 nucleotides (nt), ubiquitously expressed small RNAs that are sequentially processed from hairpin-like precursors by Drosha/Pasha and Dcr-1/Loquacious complexes. MicroRNAs usually associate with AGO1 and regulate the expression of protein-coding genes. Piwi-interacting RNAs (piRNAs) of ∼24–28 nt associate with Piwi-family proteins and can arise from single-stranded precursors. piRNAs function in transposon silencing and are mainly restricted to gonadal tissues. Endo-siRNAs are found in both germline and somatic tissues. These ∼21-nt RNAs are produced by a distinct Dicer, Dcr-2, and do not depend on Drosha/Pasha complexes. They predominantly bind to AGO2 and target both mobile elements and protein-coding genes. Surprisingly, a subset of endo-siRNAs strongly depend for their production on the dsRNA-binding protein Loquacious (Loqs), thought generally to be a partner for Dcr-1 and a cofactor for miRNA biogenesis. Endo-siRNA production depends on a specific Loqs isoform, Loqs-PD, which is distinct from the one, Loqs-PB, required for the production of microRNAs. Paralleling their roles in the biogenesis of distinct small RNA classes, Loqs-PD and Loqs-PB bind to different Dicer proteins, with Dcr-1/Loqs-PB complexes and Dcr-2/Loqs-PD complexes driving microRNA and endo-siRNA biogenesis, respectively.     

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.     

The hub protein loquacious connects the microRNA and short interfering RNA pathways in mosquitoes

Aedes aegypti mosquitoes vector several arboviruses of global health significance, including dengue viruses and chikungunya virus. RNA interference (RNAi) plays an important role in antiviral immunity, gene regulation and protection from transposable elements. Double-stranded RNA binding proteins (dsRBPs) are important for efficient RNAi; in Drosophila functional specialization of the miRNA, endo-siRNA and exo-siRNA pathway is aided by the dsRBPs Loquacious (Loqs-PB, Loqs-PD) and R2D2, respectively. However, this functional specialization has not been investigated in other dipterans. We were unable to detect Loqs-PD in Ae. aegypti; analysis of other dipteran genomes demonstrated that this isoform is not conserved outside of Drosophila. Overexpression experiments and small RNA sequencing following depletion of each dsRBP revealed that R2D2 and Loqs-PA cooperate non-redundantly in siRNA production, and that these proteins exhibit an inhibitory effect on miRNA levels. Conversely, Loqs-PB alone interacted with mosquito dicer-1 and was essential for full miRNA production. Mosquito Loqs interacted with both argonaute 1 and 2 in a manner independent of its interactions with dicer. We conclude that the functional specialization of Loqs-PD in Drosophila is a recently derived trait, and that in other dipterans, including the medically important mosquitoes, Loqs-PA participates in both the miRNA and endo-siRNA based pathways.     

Transposon Defense by Endo-siRNAs,piRNAs and Somatic pilRNAs in Drosophila: Contributions of Loqs-PD and R2D2

Transposable elements are a serious threat for genome integrity and their control via small RNA mediated silencing pathways is an ancient strategy. The fruit fly Drosophila melanogaster has two silencing activities that target transposons: endogenous siRNAs (esiRNAs or endo-siRNAs) and Piwi-interacting small RNAs (piRNAs). The biogenesis of endo-siRNAs involves the Dicer-2 co-factors Loqs-PD, which acts predominantly during processing of dsRNA by Dcr-2, and R2D2, which primarily helps to direct siRNAs into the RNA interference effector Ago2. Nonetheless, loss of either protein is not sufficient to produce a phenotype comparable with a dcr-2 mutation. We provide further deep sequencing evidence supporting the notion that R2D2 and Loqs-PD have partially overlapping function. Certain transposons display a preference for either dsRBD-protein during production or loading; this appeared to correlate neither with overall abundance, classification of the transposon or a specific site of genomic origin. The endo-siRNA biogenesis pathway in germline operates according to the same principles as the existing model for the soma, and its impairment does not significantly affect piRNAs. Expanding the analysis, we confirmed the occurrence of somatic piRNA-like RNAs (pilRNAs) that show a ping-pong signature. We detected expression of the Piwi-family protein mRNAs only barely above background, indicating that the somatic pilRNAs may arise from a small sub-population of somatic cells that express a functional piRNA pathway.     

The Bombyx ovary-derived cell line endogenously expresses PIWI/PIWI-interacting RNA complexes

Genetic studies and large-scale sequencing experiments have revealed that the PIWI subfamily proteins and PIWI-interacting RNAs (piRNAs) play an important role in germ line development and transposon control. Biochemical studies in vitro have greatly contributed to the understanding of small interfering RNA (siRNA) and microRNA (miRNA) pathways. However, in vitro analyses of the piRNA pathway have been thus far quite challenging, because their expression is largely restricted to the germ line. Here we report that Bombyx mori ovary-derived cultured cell line, BmN4, endogenously expresses two PIWI subfamily proteins, silkworm Piwi (Siwi) and Ago3 (BmAgo3), and piRNAs associated with them. Siwi-bound piRNAs have a strong bias for uridine at their 5′ end and BmAgo3-bound piRNAs are enriched for adenine at position 10. In addition, Siwi preferentially binds antisense piRNAs, whereas BmAgo3 binds sense piRNAs. Moreover, we identified many pairs in which Siwi-bound antisense and BmAgo3-bound sense piRNAs are overlapped by precisely 10 nt at their 5′ ends. These signatures are known to be important for secondary piRNA biogenesis in other organisms. Taken together, BmN4 is a unique cell line in which both primary and secondary steps of piRNA biogenesis pathways are active. This cell line would provide useful tools for analysis of piRNA biogenesis and function.     

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.     

Dicer-1, but not Loquacious, is critical for assembly of miRNA-induced silencing complexes

Double-stranded RNA-binding proteins (dsRBPs), such as R2D2 and Loquacious (Loqs), function in tandem with Dicer (Dcr) enzymes in RNA interference (RNAi). In Drosophila, Dcr-1/Loqs and Dcr-2/R2D2 complexes generate microRNAs (miRNAs) and small interfering RNAs (siRNAs), respectively. Although R2D2 does not regulate siRNA production, R2D2 and Dcr-2 coordinately bind siRNAs to promote assembly of the siRNA-induced silencing (siRISC) complexes. Conversely, Loqs enhances miRNA production. It is uncertain if Dcr-1 and Loqs facilitate miRNA loading onto the miRISC complexes. Here we used loqs knockout (KO) flies to characterize the physiological functions of Loqs in the miRNA pathway. Northern analysis revealed consistent accumulation of precursor (pre)-miRNAs in loqs(KO) flies. However, the lack of Loqs had differential effects on mature miRNAs: some are diminished, whereas others maintain wild-type levels. Importantly, the data suggest that miRNA production is not the rate-limiting step of the miRNA pathway. We show that Dcr-1, but not Loqs, is critical for assembly of miRISCs by using dcr-1 or loqs null egg extract. Consistent with this, recombinant Dcr-1 could efficiently interact with miRNA duplex in the absence of Loqs. Together, our results indicate that Loqs plays a prominent role in miRNA biogenesis, but is largely dispensable for miRISC assembly. Thus, Loqs and R2D2 represent two distinct functional modes for dsRBPs in the RNAi pathways.     

The impact of age,biogenesis, and genomic clustering on Drosophila microRNA evolution

The molecular evolutionary signatures of miRNAs inform our understanding of their emergence, biogenesis, and function. The known signatures of miRNA evolution have derived mostly from the analysis of deeply conserved, canonical loci. In this study, we examine the impact of age, biogenesis pathway, and genomic arrangement on the evolutionary properties of Drosophila miRNAs. Crucial to the accuracy of our results was our curation of high-quality miRNA alignments, which included nearly 150 corrections to ortholog calls and nucleotide sequences of the global 12-way Drosophilid alignments currently available. Using these data, we studied primary sequence conservation, normalized free-energy values, and types of structure-preserving substitutions. We expand upon common miRNA evolutionary patterns that reflect fundamental features of miRNAs that are under functional selection. We observe that melanogaster-subgroup-specific miRNAs, although recently emerged and rapidly evolving, nonetheless exhibit evolutionary signatures that are similar to well-conserved miRNAs and distinct from other structured noncoding RNAs and bulk conserved non-miRNA hairpins. This provides evidence that even young miRNAs may be selected for regulatory activities. More strikingly, we observe that mirtrons and clustered miRNAs both exhibit distinct evolutionary properties relative to solo, well-conserved miRNAs, even after controlling for sequence depth. These studies highlight the previously unappreciated impact of biogenesis strategy and genomic location on the evolutionary dynamics of miRNAs, and affirm that miRNAs do not evolve as a unitary class.     

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.     

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.     

Structural and functional insights into the fly microRNA biogenesis factor Loquacious

In the microRNA (miRNA) pathway, Dicer processes precursors to mature miRNAs. For efficient processing, double-stranded RNA-binding proteins support Dicer proteins. In flies, Loquacious (Loqs) interacts with Dicer1 (dmDcr1) to facilitate miRNA processing. Here, we have solved the structure of the third double-stranded RNA-binding domain (dsRBD) of Loqs and define specific structural elements that interact with dmDcr1. In addition, we show that the linker preceding dsRBD3 contributes significantly to dmDcr1 binding. Furthermore, our structural work demonstrates that the third dsRBD of Loqs forms homodimers. Mutations in the dimerization interface abrogate dmDcr1 interaction. Loqs, however, binds to dmDcr1 as a monomer using the identified dimerization surface, which suggests that Loqs might form dimers under conditions where dmDcr1 is absent or not accessible. Since critical sequence elements are conserved, we suggest that dimerization might be a general feature of dsRBD proteins in gene silencing.