MicroRNAs visit the ER.

Little is known about where microRNAs (miRNAs) regulate their targets within the cell. In this issue, Li et al. identify a new player in the plant miRNA pathway that implicates the endoplasmic reticulum in miRNA-mediated gene silencing.

Specificity of ARGONAUTE7-miR390 interaction and dual functionality in TAS3 trans-acting siRNA formation.

Trans-acting siRNA form through a refined RNAi mechanism in plants. miRNA-guided cleavage triggers entry of precursor transcripts into an RNA-DEPENDENT RNA POLYMERASE6 pathway, and sets the register for phased tasiRNA formation by DICER-LIKE4. Here, we show that miR390-ARGONAUTE7 complexes function in distinct cleavage or noncleavage modes at two target sites in TAS3a transcripts. The AGO7 cleavage, but not the noncleavage, function could be provided by AGO1, the dominant miRNA-associated AGO, but only when AGO1 was guided to a modified target site through an alternate miRNA. AGO7 was highly selective for interaction with miR390, and miR390 in turn was excluded from association with AGO1 due entirely to an incompatible 5' adenosine. Analysis of AGO1, AGO2, and AGO7 revealed a potent 5' nucleotide discrimination function for some, although not all, ARGONAUTEs. miR390 and AGO7, therefore, evolved as a highly specific miRNA guide/effector protein pair to function at two distinct tasiRNA biogenesis steps.

Widespread roles for piRNAs and WAGO-class siRNAs in shaping the germline transcriptome of Caenorhabditis elegans.

Piwi-interacting RNAs (piRNAs) and small interfering RNAs (siRNAs) are distinct classes of small RNAs required for proper germline development. To identify the roles of piRNAs and siRNAs in regulating gene expression in Caenorhabditis elegans, we subjected small RNAs and mRNAs from the gonads of piRNA and siRNA defective mutants to high-throughput sequencing. We show that piRNAs and an abundant class of siRNAs known as WAGO-class 22G-RNAs are required for proper expression of spermatogenic and oogenic genes. WAGO-class 22G-RNAs are also broadly required for transposon silencing, whereas piRNAs are largely dispensable. piRNAs, however, have a critical role in controlling histone gene expression. In the absence of piRNAs, histone mRNAs are misrouted into the nuclear RNAi pathway involving the Argonaute HRDE-1, concurrent with a reduction in the expression of many histone mRNAs. We also show that high-level gene expression in the germline is correlated with high level 22G-RNA production. However, most highly expressed genes produce 22G-RNAs through a distinct pathway that presumably involves the Argonaute CSR-1. In contrast, genes targeted by the WAGO branch of the 22G-RNA pathway are typically poorly expressed and respond unpredictably to loss of 22G-RNAs. Our results point to broad roles for piRNAs and siRNAs in controlling gene expression in the C. elegans germline.

MUT-16 promotes formation of perinuclear mutator foci required for RNA silencing in the C. elegans germline.

RNA silencing can be initiated by endogenous or exogenously delivered siRNAs. In Caenorhabditis elegans, RNA silencing guided by primary siRNAs is inefficient and therefore requires an siRNA amplification step involving RNA-dependent RNA polymerases (RdRPs). Many factors involved in RNA silencing localize to protein- and RNA-rich nuclear pore-associated P granules in the germline, where they are thought to surveil mRNAs as they exit the nucleus. Mutator class genes are required for siRNA-mediated RNA silencing in both germline and somatic cells, but their specific roles and relationship to other siRNA factors are unclear. Here we show that each of the six mutator proteins localizes to punctate foci at the periphery of germline nuclei. The Mutator foci are adjacent to P granules but are not dependent on core P-granule components or other RNAi pathway factors for their formation or stability. The glutamine/asparagine (Q/N)-rich protein MUT-16 is specifically required for the formation of a protein complex containing the mutator proteins, and in its absence, Mutator foci fail to form at the nuclear periphery. The RdRP RRF-1 colocalizes with MUT-16 at Mutator foci, suggesting a role for Mutator foci in siRNA amplification. Furthermore, we demonstrate that genes that yield high levels of siRNAs, indicative of multiple rounds of siRNA amplification, are disproportionally affected in mut-16 mutants compared with genes that yield low levels of siRNAs. We propose that the mutator proteins and RRF-1 constitute an RNA processing compartment required for siRNA amplification and RNA silencing.

ALG-5 is a miRNA-associated Argonaute required for proper developmental timing in the Caenorhabditis elegans germline.

Caenorhabditis elegans contains 25 Argonautes, of which, ALG-1 and ALG-2 are known to primarily interact with miRNAs. ALG-5 belongs to the AGO subfamily of Argonautes that includes ALG-1 and ALG-2, but its role in small RNA pathways is unknown. We analyzed by high-throughput sequencing the small RNAs associated with ALG-5, ALG-1 and ALG-2, as well as changes in mRNA expression in alg-5, alg-1 and alg-2 mutants. We show that ALG-5 defines a distinct branch of the miRNA pathway affecting the expression of genes involved in immunity, defense, and development. In contrast to ALG-1 and ALG-2, which associate with most miRNAs and have general roles throughout development, ALG-5 interacts with only a small subset of miRNAs and is specifically expressed in the germline where it localizes alongside the piRNA and siRNA machinery at P granules. alg-5 is required for optimal fertility and mutations in alg-5 lead to a precocious transition from spermatogenesis to oogenesis. Our results provide a near-comprehensive analysis of miRNA-Argonaute interactions in C. elegans and reveal a new role for miRNAs in the germline.

Pattern formation via small RNA mobility.

MicroRNAs and trans-acting siRNAs (ta-siRNAs) have important regulatory roles in development. Unlike other developmentally important regulatory molecules, small RNAs are not known to act as mobile signals during development. Here, we show that low-abundant, conserved ta-siRNAs, termed tasiR-ARFs, move intercellularly from their defined source of biogenesis on the upper (adaxial) side of leaves to the lower (abaxial) side to create a gradient of small RNAs that patterns the abaxial determinant AUXIN RESPONSE FACTOR3. Our observations have important ramifications for the function of small RNAs and suggest they can serve as mobile, instructive signals during development.

High-throughput sequencing reveals extraordinary fluidity of miRNA, piRNA, and siRNA pathways in nematodes.

The nematode Caenorhabditis elegans contains each of the broad classes of eukaryotic small RNAs, including microRNAs (miRNAs), endogenous small-interfering RNAs (endo-siRNAs), and piwi-interacting RNAs (piRNAs). To better understand the evolution of these regulatory RNAs, we deep-sequenced small RNAs from C. elegans and three closely related nematodes: C. briggsae, C. remanei, and C. brenneri. The results reveal a fluid landscape of small RNA pathways with essentially no conservation of individual sequences aside from a subset of miRNAs. We identified 54 miRNA families that are conserved in each of the four species, as well as numerous miRNAs that are species-specific or shared between only two or three species. Despite a lack of conservation of individual piRNAs and siRNAs, many of the features of each pathway are conserved between the different species. We show that the genomic distribution of 26G siRNAs and the tendency for piRNAs to cluster is conserved between C. briggsae and C. elegans. We also show that, in each species, 26G siRNAs trigger stage-specific secondary siRNA formation. piRNAs in each species also trigger secondary siRNA formation from targets containing up to three mismatches. Finally, we show that the production of male- and female-specific piRNAs is conserved in all four species, suggesting distinct roles for piRNAs in male and female germlines.

Functional analysis of three Arabidopsis ARGONAUTES using slicer-defective mutants.

In RNA-directed silencing pathways, ternary complexes result from small RNA-guided ARGONAUTE (AGO) associating with target transcripts. Target transcripts are often silenced through direct cleavage (slicing), destabilization through slicer-independent turnover mechanisms, and translational repression. Here, wild-type and active-site defective forms of several Arabidopsis thaliana AGO proteins involved in posttranscriptional silencing were used to examine several AGO functions, including small RNA binding, interaction with target RNA, slicing or destabilization of target RNA, secondary small interfering RNA formation, and antiviral activity. Complementation analyses in ago mutant plants revealed that the catalytic residues of AGO1, AGO2, and AGO7 are required to restore the defects of Arabidopsis ago1-25, ago2-1, and zip-1 (AGO7-defective) mutants, respectively. AGO2 had slicer activity in transient assays but could not trigger secondary small interfering RNA biogenesis, and catalytically active AGO2 was necessary for local and systemic antiviral activity against Turnip mosaic virus. Slicer-defective AGOs associated with miRNAs and stabilized AGO-miRNA-target RNA ternary complexes in individual target coimmunoprecipitation assays. In genome-wide AGO-miRNA-target RNA coimmunoprecipitation experiments, slicer-defective AGO1-miRNA associated with target RNA more effectively than did wild-type AGO1-miRNA. These data not only reveal functional roles for AGO1, AGO2, and AGO7 slicer activity, but also indicate an approach to capture ternary complexes more efficiently for genome-wide analyses.

PIWI associated siRNAs and piRNAs specifically require the Caenorhabditis elegans HEN1 ortholog henn-1.

Small RNAs--including piRNAs, miRNAs, and endogenous siRNAs--bind Argonaute proteins to form RNA silencing complexes that target coding genes, transposons, and aberrant RNAs. To assess the requirements for endogenous siRNA formation and activity in Caenorhabditis elegans, we developed a GFP-based sensor for the endogenous siRNA 22G siR-1, one of a set of abundant siRNAs processed from a precursor RNA mapping to the X chromosome, the X-cluster. Silencing of the sensor is also dependent on the partially complementary, unlinked 26G siR-O7 siRNA. We show that 26G siR-O7 acts in trans to initiate 22G siRNA formation from the X-cluster. The presence of several mispairs between 26G siR-O7 and the X-cluster mRNA, as well as mutagenesis of the siRNA sensor, indicates that siRNA target recognition is permissive to a degree of mispairing. From a candidate reverse genetic screen, we identified several factors required for 22G siR-1 activity, including the chromatin factors mes-4 and gfl-1, the Argonaute ergo-1, and the 3' methyltransferase henn-1. Quantitative RT-PCR of small RNAs in a henn-1 mutant and deep sequencing of methylated small RNAs indicate that siRNAs and piRNAs that associate with PIWI clade Argonautes are methylated by HENN-1, while siRNAs and miRNAs that associate with non-PIWI clade Argonautes are not. Thus, PIWI-class Argonaute proteins are specifically adapted to associate with methylated small RNAs in C. elegans.

mut-16 and other mutator class genes modulate 22G and 26G siRNA pathways in Caenorhabditis elegans.

Argonaute-associated siRNAs and Piwi-associated piRNAs have overlapping roles in silencing mobile genetic elements in animals. In Caenorhabditis elegans, mutator (mut) class genes mediate siRNA-guided repression of transposons as well as exogenous RNAi, but their roles in endogenous RNA silencing pathways are not well-understood. To characterize the endogenous small RNAs dependent on mut class genes, small RNA populations from a null allele of mut-16 as well as a regulatory mut-16(mg461) allele that disables only somatic RNAi were subjected to deep sequencing. Additionally, each of the mut class genes was tested for a requirement in 26G siRNA pathways. The results indicate that mut-16 is an essential factor in multiple endogenous germline and somatic siRNA pathways involving several distinct Argonautes and RNA-dependent RNA polymerases. The results also reveal essential roles for mut-2 and mut-7 in the ERGO-1 class 26G siRNA pathway and less critical roles for mut-8, mut-14, and mut-15. We show that transposons are hypersusceptible to mut-16-dependent silencing and identify a requirement for the siRNA machinery in piRNA biogenesis from Tc1 transposons. We also show that the soma-specific mut-16(mg461) mutant allele is present in multiple C. elegans laboratory strains.

The ERI-6/7 helicase acts at the first stage of an siRNA amplification pathway that targets recent gene duplications.

Endogenous small interfering RNAs (siRNAs) are a class of naturally occuring regulatory RNAs found in fungi, plants, and animals. Some endogenous siRNAs are required to silence transposons or function in chromosome segregation; however, the specific roles of most endogenous siRNAs are unclear. The helicase gene eri-6/7 was identified in the nematode Caenorhabditis elegans by the enhanced response to exogenous double-stranded RNAs (dsRNAs) of the null mutant. eri-6/7 encodes a helicase homologous to small RNA factors Armitage in Drosophila, SDE3 in Arabidopsis, and Mov10 in humans. Here we show that eri-6/7 mutations cause the loss of 26-nucleotide (nt) endogenous siRNAs derived from genes and pseudogenes in oocytes and embryos, as well as deficiencies in somatic 22-nucleotide secondary siRNAs corresponding to the same loci. About 80 genes are eri-6/7 targets that generate the embryonic endogenous siRNAs that silence the corresponding mRNAs. These 80 genes share extensive nucleotide sequence homology and are poorly conserved, suggesting a role for these endogenous siRNAs in silencing of and thereby directing the fate of recently acquired, duplicated genes. Unlike most endogenous siRNAs in C. elegans, eri-6/7-dependent siRNAs require Dicer. We identify that the eri-6/7-dependent siRNAs have a passenger strand that is ∼19 nt and is inset by ∼3-4 nts from both ends of the 26 nt guide siRNA, suggesting non-canonical Dicer processing. Mutations in the Argonaute ERGO-1, which associates with eri-6/7-dependent 26 nt siRNAs, cause passenger strand stabilization, indicating that ERGO-1 is required to separate the siRNA duplex, presumably through endonucleolytic cleavage of the passenger strand. Thus, like several other siRNA-associated Argonautes with a conserved RNaseH motif, ERGO-1 appears to be required for siRNA maturation.

Regulation of Microprocessor assembly and localization via Pasha's WW domain in C. elegans.

Primary microRNA (pri-miRNA) transcripts are processed by the Microprocessor, a protein complex that includes the ribonuclease Drosha and its RNA binding partner DGCR8/Pasha. We developed a live, whole animal, fluorescence-based sensor that reliably monitors pri-miRNA processing with high sensitivity in C. elegans. Through a forward genetic selection for alleles that desilence the sensor, we identified a mutation in the conserved G residue adjacent to the namesake W residue of Pasha's WW domain. Using genome editing we also mutated the W residue and reveal that both the G and W residue are required for dimerization of Pasha and proper assembly of the Microprocessor. Surprisingly, we find that the WW domain also facilitates nuclear localization of Pasha, which in turn promotes nuclear import or retention of Drosha. Furthermore, depletion of Pasha or Drosha causes both components of the Microprocessor to mislocalize to the cytoplasm. Thus, Pasha and Drosha mutually regulate each other's spatial expression in C. elegans.

Identification of MIR390a precursor processing-defective mutants in Arabidopsis by direct genome sequencing.

Transacting siRNA (tasiRNA) biogenesis in Arabidopsis is initiated by microRNA (miRNA) -guided cleavage of primary transcripts. In the case of TAS3 tasiRNA formation, ARGONAUTE7 (AGO7)-miR390 complexes interact with primary transcripts at two sites, resulting in recruitment of RNA-DEPENDENT RNA POLYMERASE6 for dsRNA biosynthesis. An extensive screen for Arabidopsis mutants with specific defects in TAS3 tasiRNA biogenesis or function was done. This yielded numerous ago7 mutants, one dcl4 mutant, and two mutants that accumulated low levels of miR390. A direct genome sequencing-based approach to both map and rapidly identify one of the latter mutant alleles was developed. This revealed a G-to-A point mutation (mir390a-1) that was calculated to stabilize a relatively nonpaired region near the base of the MIR390a foldback, resulting in misprocessing of the miR390/miR390* duplex and subsequent reduced TAS3 tasiRNA levels. Directed substitutions, as well as analysis of variation at paralogous miR390-generating loci (MIR390a and MIR390b), indicated that base pair properties and nucleotide identity within a region 4-6 bases below the miR390/miR390* duplex region contributed to the efficiency and accuracy of precursor processing.

tiny-count: a counting tool for hierarchical classification and quantification of small RNA-seq reads with single-nucleotide precision.

Summary: tiny-count is a highly flexible counting tool that allows for hierarchical classification and quantification of small RNA reads from high-throughput sequencing data. Selection rules can be used to filter reads by 5' nucleotide, length, position of alignments in relation to reference features, and by the number of mismatches to reference sequences. tiny-count can quantify reads aligned to a genome or directly to small RNA or transcript sequences. With tiny-count, users can quantify a single class of small RNAs or multiple classes in parallel. tiny-count can resolve distinct classes of small RNAs, for example, piRNAs and siRNAs, produced from the same locus. It can distinguish small RNA variants, such as miRNAs and isomiRs, with single-nucleotide precision. tRNA, rRNA, and other RNA fragments can also be quantified. tiny-count can be run alone or as part of tinyRNA, a workflow that provides a basic all-in-one command line-based solution for small RNA-seq data analysis, with documentation and statistics generated at each step for accurate and reproducible results. Availability and implementation: tiny-count and other tinyRNA tools are implemented in Python, C++, Cython, and R, and the workflow is coordinated with CWL. tiny-count and tinyRNA are free and open-source software distributed under the GPLv3 license. tiny-count can be installed via Bioconda (https://anaconda.org/bioconda/tiny-count) and both tiny-count and tinyRNA documentation and software downloads are available at https://github.com/MontgomeryLab/tinyRNA. Reference data, including genome and feature information, for certain species can be found at https://www.MontgomeryLab.org.

Imaging translational control by Argonaute with single-molecule resolution in live cells.

A major challenge to our understanding of translational control has been deconvolving the individual impact specific regulatory factors have on the complex dynamics of mRNA translation. MicroRNAs (miRNAs), for example, guide Argonaute and associated proteins to target mRNAs, where they direct gene silencing in multiple ways that are not well understood. To better deconvolve these dynamics, we have developed technology to directly visualize and quantify the impact of human Argonaute2 (Ago2) on the translation and subcellular localization of individual reporter mRNAs in living cells. We show that our combined translation and Ago2 tethering sensor reflects endogenous miRNA-mediated gene silencing. Using the sensor, we find that Ago2 association leads to progressive silencing of translation at individual mRNA. Silencing was occasionally interrupted by brief bursts of translational activity and took 3-4 times longer than a single round of translation, consistent with a gradual increase in the inhibition of translation initiation. At later time points, Ago2-tethered mRNAs cluster and coalesce with P-bodies, where a translationally silent state is maintained. These results provide a framework for exploring miRNA-mediated gene regulation in live cells at the single-molecule level. Furthermore, our tethering-based, single-molecule reporter system will likely have wide-ranging application in studying RNA-protein interactions.

Computational and analytical framework for small RNA profiling by high-throughput sequencing.

The advent of high-throughput sequencing (HTS) methods has enabled direct approaches to quantitatively profile small RNA populations. However, these methods have been limited by several factors, including representational artifacts and lack of established statistical methods of analysis. Furthermore, massive HTS data sets present new problems related to data processing and mapping to a reference genome. Here, we show that cluster-based sequencing-by-synthesis technology is highly reproducible as a quantitative profiling tool for several classes of small RNA from Arabidopsis thaliana. We introduce the use of synthetic RNA oligoribonucleotide standards to facilitate objective normalization between HTS data sets, and adapt microarray-type methods for statistical analysis of multiple samples. These methods were tested successfully using mutants with small RNA biogenesis (miRNA-defective dcl1 mutant and siRNA-defective dcl2 dcl3 dcl4 triple mutant) or effector protein (ago1 mutant) deficiencies. Computational methods were also developed to rapidly and accurately parse, quantify, and map small RNA data.

AGO1-miR173 complex initiates phased siRNA formation in plants.

MicroRNA (miRNA)-guided cleavage initiates entry of primary transcripts into the transacting siRNA (tasiRNA) biogenesis pathway involving RNA-DEPENDENT RNA POLYMERASE6, DICER-LIKE4, and SUPPRESSOR OF GENE SILENCING3. Arabidopsis thaliana TAS1 and TAS2 families yield tasiRNA that form through miR173-guided initiation-cleavage of primary transcripts and target several transcripts encoding pentatricopeptide repeat proteins and proteins of unknown function. Here, the TAS1c locus was modified to produce synthetic (syn) tasiRNA to target an endogenous transcript encoding PHYTOENE DESATURASE and used to analyze the role of miR173 in routing of transcripts through the tasiRNA pathway. miR173 was unique from other miRNAs in its ability to initiate TAS1c-based syn-tasiRNA formation. A single miR173 target site was sufficient to route non-TAS transcripts into the pathway to yield phased siRNA. We also show that miR173 functions in association with ARGONAUTE 1 (AGO1) during TAS1 and TAS2 tasiRNA formation, and we provide data indicating that the miR173-AGO1 complex possesses unique functionality that many other miRNA-AGO1 complexes lack.

Dual roles for piRNAs in promoting and preventing gene silencing in C. elegans.

Piwi-interacting RNAs (piRNAs) regulate many biological processes through mechanisms that are not fully understood. In Caenorhabditis elegans, piRNAs intersect the endogenous RNA interference (RNAi) pathway, involving a distinct class of small RNAs called 22G-RNAs, to regulate gene expression in the germline. In the absence of piRNAs, 22G-RNA production from many genes is reduced, pointing to a role for piRNAs in facilitating endogenous RNAi. Here, however, we show that many genes gain, rather than lose, 22G-RNAs in the absence of piRNAs, which is in some instances coincident with RNA silencing. Aberrant 22G-RNA production is somewhat stochastic but once established can occur within a population for at least 50 generations. Thus, piRNAs both promote and suppress 22G-RNA production and gene silencing. rRNAs and histones are hypersusceptible to aberrant silencing, but we do not find evidence that their misexpression is the primary cause of the transgenerational sterility observed in piRNA-defective mutants.