Small silencing RNAs, including microRNAs, endogenous small interfering RNAs (endo-siRNAs) and Piwi-interacting RNAs (piRNAs), have been shown to play important roles in fine-tuning gene expression, defending virus and controlling transposons. Loss of small silencing RNAs or components in their pathways often leads to severe developmental defects, including lethality and sterility. Recently, non-templated addition of nucleotides to the 3' end, namely tailing, was found to associate with the processing and stability of small silencing RNAs. Next Generation Sequencing has made it possible to detect such modifications at nucleotide resolution in an unprecedented throughput. Unfortunately, detecting such events from millions of short reads confounded by sequencing errors and RNA editing is still a tricky problem. Here, we developed a computational framework, Tailor, driven by an efficient and accurate aligner specifically designed for capturing the tailing events directly from the alignments without extensive post-processing. The performance of Tailor was fully tested and compared favorably with other general-purpose aligners using both simulated and real datasets for tailing analysis. Moreover, to show the broad utility of Tailor, we used Tailor to reanalyze published datasets and revealed novel findings worth further experimental validation. The source code and the executable binaries are freely available at https://github.com/jhhung/Tailor.

PIWI-interacting RNAs (piRNAs) silence transposons in animal germ cells. PIWI proteins bind and amplify piRNAs via the "Ping-Pong" pathway. Because PIWI proteins cleave RNAs between target nucleotides t10 and t11-the nucleotides paired to piRNA guide positions g10 and g11-the first ten nucleotides of piRNAs participating in the Ping-Pong amplification cycle are complementary. Drosophila piRNAs bound to the PIWI protein Aubergine typically begin with uridine (1U), while piRNAs bound to Argonaute3, which are produced by Ping-Pong amplification, often have adenine at position 10 (10A). The Ping-Pong model proposes that the 10A is a consequence of 1U. We find that 10A is not caused by 1U. Instead, fly Aubergine as well as its homologs, Siwi in silkmoth and MILI in mice, have an intrinsic preference for adenine at the t1 position of their target RNAs; during Ping-Pong amplification, this t1A subsequently becomes the g10A of a piRNA bound to Argonaute3.

Adoption of a streamlined version of the bacterial clustered regular interspersed short palindromic repeat (CRISPR)/Cas9 defense system has accelerated targeted genome engineering. The Streptococcus pyogenes Cas9 protein, directed by a simplified, CRISPR-like single-guide RNA, catalyzes a double-stranded DNA break at a specific genomic site; subsequent repair by end joining can introduce mutagenic insertions or deletions, while repair by homologous recombination using an exogenous DNA template can incorporate new sequences at the target locus. However, the efficiency of Cas9-directed mutagenesis is low in Drosophila melanogaster Here, we describe a strategy that reduces the time and effort required to identify flies with targeted genomic changes. The strategy uses editing of the white gene, evidenced by altered eye color, to predict successful editing of an unrelated gene-of-interest. The red eyes of wild-type flies are readily distinguished from white-eyed (end-joining-mediated loss of White function) or brown-eyed (recombination-mediated conversion to the whitecoffee allele) mutant flies. When single injected G0 flies produce individual G1 broods, flies carrying edits at a gene-of-interest were readily found in broods in which all G1 offspring carried white mutations. Thus, visual assessment of eye color substitutes for wholesale PCR screening of large numbers of G1 offspring. We find that end-joining-mediated mutations often show signatures of microhomology-mediated repair and that recombination-based mutations frequently involve donor plasmid integration at the target locus. Finally, we show that gap repair induced by two guide RNAs more reliably converts the intervening target sequence, whereas the use of Lig4169 mutants to suppress end joining does not improve recombination efficacy.

We report a genome-editing strategy to correct compound heterozygous mutations, a common genotype in patients with recessive genetic disorders. Adeno-associated viral vector delivery of Cas9 and guide RNA induces allelic exchange and rescues the disease phenotype in mouse models of hereditary tyrosinemia type I and mucopolysaccharidosis type I. This approach recombines non-mutated genetic information present in two heterozygous alleles into one functional allele without using donor DNA templates.

Argonaute proteins loaded with microRNAs (miRNAs) or small interfering RNAs (siRNAs) form the RNA-induced silencing complex (RISC), which represses target RNA expression. Predicting the biological targets, specificity, and efficiency of both miRNAs and siRNAs has been hamstrung by an incomplete understanding of the sequence determinants of RISC binding and cleavage. We applied high-throughput methods to measure the association kinetics, equilibrium binding energies, and single-turnover cleavage rates of mouse AGO2 RISC. We find that RISC readily tolerates insertions of up to 7 nt in its target opposite the central region of the guide. Our data uncover specific guide:target mismatches that enhance the rate of target cleavage, suggesting novel siRNA design strategies. Using these data, we derive quantitative models for RISC binding and target cleavage and show that our in vitro measurements and models predict knockdown in an engineered cellular system.

In animals, small RNA molecules termed PIWI-interacting RNAs (piRNAs) silence transposable elements (TEs), protecting the germline from genomic instability and mutation. piRNAs have been detected in the soma in a few animals, but these are believed to be specific adaptations of individual species. Here, we report that somatic piRNAs were probably present in the ancestral arthropod more than 500 million years ago. Analysis of 20 species across the arthropod phylum suggests that somatic piRNAs targeting TEs and messenger RNAs are common among arthropods. The presence of an RNA-dependent RNA polymerase in chelicerates (horseshoe crabs, spiders and scorpions) suggests that arthropods originally used a plant-like RNA interference mechanism to silence TEs. Our results call into question the view that the ancestral role of the piRNA pathway was to protect the germline and demonstrate that small RNA silencing pathways have been repurposed for both somatic and germline functions throughout arthropod evolution.

In Drosophila, 23-30 nt long PIWI-interacting RNAs (piRNAs) direct the protein Piwi to silence germline transposon transcription. Most germline piRNAs derive from dual-strand piRNA clusters, heterochromatic transposon graveyards that are transcribed from both genomic strands. These piRNA sources are marked by the heterochromatin protein 1 homolog Rhino (Rhi), which facilitates their promoter-independent transcription, suppresses splicing, and inhibits transcriptional termination. Here, we report that the protein Maelstrom (Mael) represses canonical, promoter-dependent transcription in dual-strand clusters, allowing Rhi to initiate piRNA precursor transcription. Mael also represses promoter-dependent transcription at sites outside clusters. At some loci, Mael repression requires the piRNA pathway, while at others, piRNAs play no role. We propose that by repressing canonical transcription of individual transposon mRNAs, Mael helps Rhi drive non-canonical transcription of piRNA precursors without generating mRNAs encoding transposon proteins.

Single-molecule binding assays enable the study of how molecular machines assemble and function. Current algorithms can identify and locate individual molecules, but require tedious manual validation of each spot. Moreover, no solution for high-throughput analysis of single-molecule binding data exists. Here, we describe an automated pipeline to analyze single-molecule data over a wide range of experimental conditions. In addition, our method enables state estimation on multivariate Gaussian signals. We validate our approach using simulated data, and benchmark the pipeline by measuring the binding properties of the well-studied, DNA-guided DNA endonuclease, TtAgo, an Argonaute protein from the Eubacterium Thermus thermophilus. We also use the pipeline to extend our understanding of TtAgo by measuring the protein's binding kinetics at physiological temperatures and for target DNAs containing multiple, adjacent binding sites.

Small interfering RNA (siRNA)-based drugs require chemical modifications or formulation to promote stability, minimize innate immunity, and enable delivery to target tissues. Partially modified siRNAs (up to 70% of the nucleotides) provide significant stabilization in vitro and are commercially available; thus are commonly used to evaluate efficacy of bio-conjugates for in vivo delivery. In contrast, most clinically-advanced non-formulated compounds, using conjugation as a delivery strategy, are fully chemically modified (100% of nucleotides). Here, we compare partially and fully chemically modified siRNAs in conjugate mediated delivery. We show that fully modified siRNAs are retained at 100x greater levels in various tissues, independently of the nature of the conjugate or siRNA sequence, and support productive mRNA silencing. Thus, fully chemically stabilized siRNAs may provide a better platform to identify novel moieties (peptides, aptamers, small molecules) for targeted RNAi delivery.

We report a draft assembly of the genome of Hi5 cells from the lepidopteran insect pest, Trichoplusia ni, assigning 90.6% of bases to one of 28 chromosomes and predicting 14,037 protein-coding genes. Chemoreception and detoxification gene families reveal T. ni-specific gene expansions that may explain its widespread distribution and rapid adaptation to insecticides. Transcriptome and small RNA data from thorax, ovary, testis, and the germline-derived Hi5 cell line show distinct expression profiles for 295 microRNA- and >393 piRNA-producing loci, as well as 39 genes encoding small RNA pathway proteins. Nearly all of the W chromosome is devoted to piRNA production, and T. ni siRNAs are not 2´-O-methylated. To enable use of Hi5 cells as a model system, we have established genome editing and single-cell cloning protocols. The T. ni genome provides insights into pest control and allows Hi5 cells to become a new tool for studying small RNAs ex vivo.

In animals, PIWI-interacting RNAs (piRNAs) guide PIWI proteins to silence transposons and regulate gene expression. The mechanisms for making piRNAs have been proposed to differ among cell types, tissues, and animals. Our data instead suggest a single model that explains piRNA production in most animals. piRNAs initiate piRNA production by guiding PIWI proteins to slice precursor transcripts. Next, PIWI proteins direct the stepwise fragmentation of the sliced precursor transcripts, yielding tail-to-head strings of phased precursor piRNAs (pre-piRNAs). Our analyses detect evidence for this piRNA biogenesis strategy across an evolutionarily broad range of animals, including humans. Thus, PIWI proteins initiate and sustain piRNA biogenesis by the same mechanism in species whose last common ancestor predates the branching of most animal lineages. The unified model places PIWI-clade Argonautes at the center of piRNA biology and suggests that the ancestral animal-the Urmetazoan-used PIWI proteins both to generate piRNA guides and to execute piRNA function.

Recombinant adeno-associated viruses (rAAVs) that can cross the blood-brain-barrier and achieve efficient and stable transvascular gene transfer to the central nervous system (CNS) hold significant promise for treating CNS disorders. However, following intravascular delivery, these vectors also target liver, heart, skeletal muscle, and other tissues, which may cause untoward effects. To circumvent this, we used tissue-specific, endogenous microRNAs (miRNAs) to repress rAAV expression outside the CNS, by engineering perfectly complementary miRNA-binding sites into the rAAV9 genome. This approach allowed simultaneous multi-tissue regulation and CNS-directed stable transgene expression without detectably perturbing the endogenous miRNA pathway. Regulation of rAAV expression by miRNA was primarily via site-specific cleavage of the transgene mRNA, generating specific 5' and 3' mRNA fragments. Our findings promise to facilitate the development of miRNA-regulated rAAV for CNS-targeted gene delivery and other applications.

MicroRNAs (miRNAs) and small interfering RNAs (siRNAs) guide Argonaute proteins to silence mRNA expression. Argonaute binding alters the properties of an RNA guide, creating functional domains. We show that the domains established by Argonaute-the anchor, seed, central, 3' supplementary, and tail regions-have distinct biochemical properties that explain the differences between how animal miRNAs and siRNAs bind their targets. Extensive complementarity between an siRNA and its target slows the rate at which fly Argonaute2 (Ago2) binds to and dissociates from the target. Highlighting its role in antiviral defense, fly Ago2 dissociates so slowly from extensively complementary target RNAs that essentially every fully paired target is cleaved. Conversely, mouse AGO2, which mainly mediates miRNA-directed repression, dissociates rapidly and with similar rates for fully paired and seed-matched targets. Our data narrow the range of biochemically reasonable models for how Argonaute-bound siRNAs and miRNAs find, bind, and regulate their targets.

Hundreds of microRNAs (miRNAs) and endogenous small interfering RNAs (siRNAs) have been identified from both plants and animals, yet little is known about their biochemical modes of action or biological functions. Here we report that 2'-O-methyl oligonucleotides can act as irreversible, stoichiometric inhibitors of small RNA function. We show that a 2'-O-methyl oligonucleotide complementary to an siRNA can block mRNA cleavage in Drosophila embryo lysates and HeLa cell S100 extracts and in cultured human HeLa cells. In Caenorhabditis elegans, injection of the 2'-O-methyl oligonucleotide complementary to the miRNA let-7 can induce a let-7 loss-of-function phenocopy. Using an immobilized 2'-O-methyl oligonucleotide, we show that the C. elegans Argonaute proteins ALG-1 and ALG-2, which were previously implicated in let-7 function through genetic studies, are constituents of a let-7-containing protein-RNA complex. Thus, we demonstrate that 2'-O-methyl RNA oligonucleotides can provide an efficient and straightforward way to block small RNA function in vivo and furthermore can be used to identify small RNA-associated proteins that mediate RNA silencing pathways.

MicroRNAs (miRNAs) are short regulatory RNAs that direct repression of their mRNA targets. The miRNA "seed"-nucleotides 2-7-establishes target specificity by mediating target binding. Accurate processing of the miRNA 5' end is thought to be under strong selective pressure because a shift by just one nucleotide in the 5' end of a miRNA alters its seed sequence, redefining its repertoire of targets (Figure 1). Animal miRNAs are produced by the sequential cleavage of partially double-stranded precursors by the RNase III endonucleases Drosha and Dicer, thereby generating a transitory double-stranded intermediate comprising the miRNA paired to its partially complementary miRNA* strand. Here, we report that in flies, the 5' ends of miRNAs and miRNA* strands are typically more precisely defined than their 3' ends. Surprisingly, the precision of the 5' ends of both miRNA and miRNA* sequences increases after Argonaute2 (Ago2) loading. Our data imply that either many miRNA* sequences are under evolutionary pressure to maintain their seed sequences-that is, they have targets-or that secondary constraints, such as the sequence requirements for loading small RNAs into functional Argonaute complexes, narrow the range of miRNA and miRNA 5' ends that accumulate in flies.

In animals, PIWI-interacting RNAs (piRNAs) silence transposons, fight viral infections, and regulate gene expression. piRNA biogenesis concludes with 3' terminal trimming and 2'-O-methylation. Both trimming and methylation influence piRNA stability. Our biochemical data show that multiple mechanisms destabilize unmethylated mouse piRNAs, depending on whether the piRNA 5' or 3' sequence is complementary to a trigger RNA. Unlike target-directed degradation of microRNAs, complementarity-dependent destabilization of piRNAs in mice and flies is blocked by 3' terminal 2'-O-methylation and does not require base pairing to both the piRNA seed and the 3' sequence. In flies, 2'-O-methylation also protects small interfering RNAs (siRNAs) from complementarity-dependent destruction. By contrast, pre-piRNA trimming protects mouse piRNAs from a degradation pathway unaffected by trigger complementarity. In testis lysate and in vivo, internal or 3' terminal uridine- or guanine-rich tracts accelerate pre-piRNA decay. Loss of both trimming and 2'-O-methylation causes the mouse piRNA pathway to collapse, demonstrating that these modifications collaborate to stabilize piRNAs.

Many RNAs, including pre-mRNAs and long non-coding RNAs, can be thousands of nucleotides long and undergo complex post-transcriptional processing. Multiple sites of alternative splicing within a single gene exponentially increase the number of possible spliced isoforms, with most human genes currently estimated to express at least ten. To understand the mechanisms underlying these complex isoform expression patterns, methods are needed that faithfully maintain long-range exon connectivity information in individual RNA molecules. In this study, we describe SeqZip, a methodology that uses RNA-templated DNA-DNA ligation to retain and compress connectivity between distant sequences within single RNA molecules. Using this assay, we test proposed coordination between distant sites of alternative exon utilization in mouse Fn1, and we characterize the extraordinary exon diversity of Drosophila melanogaster Dscam1.

PIWI-interacting RNAs (piRNAs) protect the animal germ line by silencing transposons. Primary piRNAs, generated from transcripts of genomic transposon "junkyards" (piRNA clusters), are amplified by the "ping-pong" pathway, yielding secondary piRNAs. We report that secondary piRNAs, bound to the PIWI protein Ago3, can initiate primary piRNA production from cleaved transposon RNAs. The first ~26 nucleotides (nt) of each cleaved RNA becomes a secondary piRNA, but the subsequent ~26 nt become the first in a series of phased primary piRNAs that bind Piwi, allowing piRNAs to spread beyond the site of RNA cleavage. The ping-pong pathway increases only the abundance of piRNAs, whereas production of phased primary piRNAs from cleaved transposon RNAs adds sequence diversity to the piRNA pool, allowing adaptation to changes in transposon sequence.

Huang et al. (2013) recently reported that chromatin immunoprecipitation sequencing (ChIP-seq) reveals the genome-wide sites of occupancy by Piwi, a piRNA-guided Argonaute protein central to transposon silencing in Drosophila. Their study also reported that loss of Piwi causes widespread rewiring of transcriptional patterns, as evidenced by changes in RNA polymerase II occupancy across the genome. Here we reanalyze their data and report that the underlying deep-sequencing dataset does not support the authors' genome-wide conclusions.

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