Relaxed targeting rules help PIWI proteins silence transposons

TL;DR

PIWI proteins can silence transposons more effectively than AGO proteins by tolerating RNA mismatches, helping defend the genome without needing new RNA guides.

Contribution

PIWI proteins tolerate mismatches in RNA targets, making them more effective at silencing transposons compared to AGO proteins.

Findings

PIWI proteins can cleave transcripts with partial complementarity to their piRNA guides.

PIWI slicing does not require canonical seed pairing, unlike AGO proteins.

PIWI proteins avoid silencing host RNAs while efficiently targeting transposons.

Abstract

In eukaryotes, small RNA guides, such as small interfering RNAs and microRNAs, direct AGO-clade Argonaute proteins to regulate gene expression and defend the genome against external threats. Only animals make a second clade of Argonaute proteins: PIWI proteins. PIWI proteins use PIWI-interacting RNAs (piRNAs) to repress complementary transposon transcripts1,2. In theory, transposons could evade silencing through target site mutations that reduce piRNA complementarity. Here we report that, unlike AGO proteins, PIWI proteins efficiently cleave transcripts that are only partially paired to their piRNA guides. Examination of target binding and cleavage by mouse and sponge PIWI proteins revealed that PIWI slicing tolerates mismatches to any target nucleotide, including those flanking the scissile phosphate. Even canonical seed pairing is dispensable for PIWI binding or cleavage, unlike plant and animal AGOs, which require uninterrupted target pairing from the seed to the nucleotides past the scissile bond3,4. PIWI proteins are therefore better equipped than AGO proteins to target newly acquired or rapidly diverging endogenous transposons without recourse to new small RNA guides. Conversely, the minimum requirements for PIWI slicing are sufficient to avoid inadvertent silencing of host RNAs. Our results demonstrate the biological advantage of PIWI over AGO proteins in defending the genome against transposons and suggest an explanation for why the piRNA pathway was retained in animal evolution.

Of the two types of Argonaute proteins produced by animals, AGO and PIWI, PIWI proteins can bind RNAs with less complementarity, enabling efficient silencing of transposons without the need to produce new RNA guides.

Figures

PIWI proteins bind sites containing or lacking canonical seed pairing.a, Small RNA guides direct eukaryotic Argonaute proteins to complementary targets. nt, nucleotide. b, Binding affinities (*K*_d_ in pM) of MIWI, MILI, *Ef*Piwi and mouse AGO2 loaded with piRNA-1 for canonical and non-canonical target sites. c, Left, MIWI, MILI, *Ef*Piwi and mouse AGO2 binding affinities for targets contiguously paired from nucleotide g2. Right, relationship between binding energy Δ*G*^0^ calculated from *K*_d_ (mean of three independent trials) and predicted binding energy Δ*G*^0^. Goodness-of-fit for linear

PIWI slicing tolerates mismatches with any target nucleotide.a, Change in pre-steady-state cleavage rate for one or two mismatches (pink) between g2 and g20. For one mismatch, *n* = 456: all 19 possible positions × 3 geometries × 4 piRNAs × MILI and MIWI. For two mismatches, *n* = 1,368: all 171 possible combinations × 1 geometry × 4 piRNAs × MILI and MIWI. Box plots show the IQR and median. Statistical analysesare in Extended Data Fig. [3c](#Fig7). b, Change in pre-steady-state cleavage rate for one or two consecutive mismatches between g2 and g20 for contiguous g2–g21 or g2–g25 pairing for M

Mouse PIWI proteins cleave partially complementary targets in vivo.a, Schematic of the strategy used to identify 3′ cleavage products of piRNA-guided PIWI-catalysed slicing and to measure the fraction of targets cleaved by PIWI proteins in FACS-purified mouse primary spermatocytes. b, Fraction of cleaved MILI and MIWI targets in FACS-purified mouse primary spermatocytes for contiguous pairing from nucleotide g2. c, Fraction of cleaved targets in FACS-purified mouse primary spermatocytes for perfect matches (indicated in blue) and for pairing containing a single-nucleotide mismatch (indicated i

Determinants of PIWI slicing in vivo.a, Decision function coefficients for 400 logistic function fits (regression models) using around 3,500 distinct piRNA–target pairs detected in mouse primary spermatocytes. *n* = 16 permutations of 4 control (C57BL/6) and 4 *pi2*^*−/−*^*pi9*^*−/−*^*pi17*^*−/−*^ animals × 5-repeated × 5-fold cross validation. Box plots show IQR and median. b, Number of piRNAs and siRNAs predicted to cleave mutated versions of the *L1Md* *AI* transposon sequence. Data are median and IQR from 100 independent simulations. c, AGO and PIWI proteins use different rules to find and

c. Supplementary Data 1Source data for analyses of in vivo cleavage by MILI and MIWI (coordinates and identity of putative cleavage sites).

PIWI proteins bind sites containing or lacking canonical seed pairing.a, Small RNA guides direct eukaryotic Argonaute proteins to complementary targets. nt, nucleotide. b, Binding affinities (*K*_d_ in pM) of MIWI, MILI, *Ef*Piwi and mouse AGO2 loaded with piRNA-1 for canonical and non-canonical target sites. c, Left, MIWI, MILI, *Ef*Piwi and mouse AGO2 binding affinities for targets contiguously paired from nucleotide g2. Right, relationship between binding energy Δ*G*^0^ calculated from *K*_d_ (mean of three independent trials) and predicted binding energy Δ*G*^0^. Goodness-of-fit for linear

PIWI slicing tolerates mismatches with any target nucleotide.a, Change in pre-steady-state cleavage rate for one or two mismatches (pink) between g2 and g20. For one mismatch, *n* = 456: all 19 possible positions × 3 geometries × 4 piRNAs × MILI and MIWI. For two mismatches, *n* = 1,368: all 171 possible combinations × 1 geometry × 4 piRNAs × MILI and MIWI. Box plots show the IQR and median. Statistical analysesare in Extended Data Fig. [3c](#Fig7). b, Change in pre-steady-state cleavage rate for one or two consecutive mismatches between g2 and g20 for contiguous g2–g21 or g2–g25 pairing for M

Mouse PIWI proteins cleave partially complementary targets in vivo.a, Schematic of the strategy used to identify 3′ cleavage products of piRNA-guided PIWI-catalysed slicing and to measure the fraction of targets cleaved by PIWI proteins in FACS-purified mouse primary spermatocytes. b, Fraction of cleaved MILI and MIWI targets in FACS-purified mouse primary spermatocytes for contiguous pairing from nucleotide g2. c, Fraction of cleaved targets in FACS-purified mouse primary spermatocytes for perfect matches (indicated in blue) and for pairing containing a single-nucleotide mismatch (indicated i

Determinants of PIWI slicing in vivo.a, Decision function coefficients for 400 logistic function fits (regression models) using around 3,500 distinct piRNA–target pairs detected in mouse primary spermatocytes. *n* = 16 permutations of 4 control (C57BL/6) and 4 *pi2*^*−/−*^*pi9*^*−/−*^*pi17*^*−/−*^ animals × 5-repeated × 5-fold cross validation. Box plots show IQR and median. b, Number of piRNAs and siRNAs predicted to cleave mutated versions of the *L1Md* *AI* transposon sequence. Data are median and IQR from 100 independent simulations. c, AGO and PIWI proteins use different rules to find and

c. Supplementary Data 1Source data for analyses of in vivo cleavage by MILI and MIWI (coordinates and identity of putative cleavage sites).

Full text

Main

In prokaryotes and eukaryotes, small RNA or DNA guides direct Argonaute proteins to fight viruses, plasmids and transposons^5^, regulate gene expression^6–12^ or aid DNA replication^5^. Animals produce two distinct types of Argonaute protein: AGO and PIWI. AGO-clade proteins use small interfering RNAs (siRNAs; which are typically 21 nucleotides long) or microRNAs (miRNAs; which are most often 22 nucleotides long) to repress extensively or partially complementary transcripts^6^. AGO proteins initially find their targets through complementarity to a short 5′ region of their guide, the seed (nucleotides g2–g8; Fig. 1a). For miRNA-guided AGO proteins, seed complementarity is sufficient to repress the target RNA^6^. By contrast, PIWI-clade Argonaute proteins use piRNAs (which are 18–35 nucleotides long) as guides^1,2^. Although most eukaryotic genomes encode one or more AGO protein, only animals make PIWI proteins. With few exceptions, all animals use piRNAs to repress transposons^1,2^. The ancestral mechanism of piRNA-guided transposon silencing is PIWI-catalysed endonucleolytic cleavage (slicing) of complementary transposon RNAs in the cytoplasm^13–16^. Moreover, piRNA production itself requires piRNA-directed slicing of piRNA precursor transcripts^13,14,17–21^. In some animals, piRNAs also direct nuclear PIWI proteins to nascent transposon transcripts to silence transcription through repressive histone marks or DNA methylation^22–26^.Fig. 1PIWI proteins bind sites containing or lacking canonical seed pairing.a, Small RNA guides direct eukaryotic Argonaute proteins to complementary targets. nt, nucleotide. b, Binding affinities (Kd in pM) of MIWI, MILI, EfPiwi and mouse AGO2 loaded with piRNA-1 for canonical and non-canonical target sites. c, Left, MIWI, MILI, EfPiwi and mouse AGO2 binding affinities for targets contiguously paired from nucleotide g2. Right, relationship between binding energy ΔG^0^ calculated from Kd (mean of three independent trials) and predicted binding energy ΔG^0^. Goodness-of-fit for linear regression (r^2^) and P value for two-tailed permutation test for Pearson’s correlation are shown. All data are in Supplementary Fig. 2a. d, MIWI, MILI, EfPiwi and mouse AGO2 binding affinities for nine-nucleotide complementary stretches contiguously paired from all guide nucleotides. All data are in Supplementary Fig. 2b. Mean and standard deviation of data from three independent trials are shown (b,c (left), d).

siRNAs direct AGO proteins and piRNAs direct PIWI proteins to hydrolyse the phosphodiester bond that joins the target nucleotides opposite guide nucleotides g10 and g11 (that is, t10 and t11, respectively; Fig. 1a). Unlike siRNA-directed, AGO-catalysed transcript cleavage, efficient target slicing by PIWI proteins requires the auxiliary factor GTSF1 (ref. ^27^). GTSF1 accelerates the otherwise slow target cleavage by PIWIs by 10–100-fold, probably by stabilizing the catalytically competent conformation of PIWI proteins. Why PIWI slicing evolved to require an auxiliary protein is unknown.

Here we report that, compared with AGO proteins, the requirements for guide:target complementarity are relaxed for the PIWI proteins MILI (also known as PIWIL2) and MIWI (also known as PIWIL1) from mouse (Mus musculus) and Piwi from freshwater sponge (Ephydatia fluviatilis; hereafter denoted EfPiwi). PIWI proteins bind RNAs both with and without complementarity to the canonical 5′ seed of their guide. Both in vitro and in vivo, PIWI-catalysed slicing requires at least 15 contiguously paired nucleotides, and longer extents of complementarity tolerate guide:target mismatches at essentially any position. Unlike AGO proteins, guide pairing to any target nucleotide, including those that flank the scissile phosphate (t10 and t11), is dispensable for efficient slicing. Although pairing to at least four  piRNA 5′ terminal nucleotides facilitates target finding, in vitro and in vivo abundant piRNAs direct slicing of targets that lack 5′ complementarity. Notably, the minimum 15-nucleotide stretch of complementarity that licenses piRNA-guided target cleavage is sufficient to distinguish host transcripts from transposon RNAs. These findings suggest that the catalytic properties of PIWI proteins evolved to prevent transposons from escaping piRNA silencing through mutation while simultaneously retaining sufficient specificity to spare self-transcripts from inappropriate repression.

PIWI proteins bind without seed pairing

Mouse piRNAs guide MILI and MIWI to slice extensively complementary transposon transcripts^15,16,28^. piRNAs have also been proposed to direct MIWI to bind and regulate mRNA expression through the same mechanism by which miRNAs guide AGO proteins to their targets^29,30^. In this miRNA-like binding mode, base pairing to the canonical AGO seed sequence (guide nucleotides g2–g8; Fig. 1a) mediates the search for complementary sites and is sufficient to tether the Argonaute protein to its target RNA. We used the RNA Bind-’n-Seq method^31^ to measure the affinity of piRNA-guided PIWI Argonaute proteins for a library of 20-nucleotide-long random sequences (Extended Data Fig. 1a). We incubated the target RNA library with purified mouse MILI or MIWI or freshwater sponge EfPiwi loaded with a 5′ monophosphorylated synthetic RNA (26  nucleotides for MILI and EfPiwi, 30 nucleotides for MIWI; Extended Data Fig. 1b–d) and isolated and sequenced RNAs bound to the piRNA–PIWI protein complex (piRISC). The sequencing data were analysed using an approach that estimates the affinity (Kd) of piRISC for each binding-site type^32^.

Binding of MILI, MIWI or EfPiwi to RNAs with canonical AGO seed sites was weaker than for mammalian AGO2 proteins (Fig. 1b,c and Extended Data Fig. 2a). Compared with AGO proteins, PIWI protein affinity was 4–16-fold lower for an 8-mer (g2–g8•t1A), the canonical site type most effectively repressed by miRNA-guided AGO proteins^6^. In detail, for piRNA-1, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{AGO}}2,8{\rm{-mer}}}$$\end{document} = 10 ± 8 pM compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MILI}},8{\rm{-mer}}}$$\end{document} = 45 ± 7 pM, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MIWI}},8{\rm{-mer}}}$$\end{document} = 40 ± 20 pM and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{Ef{\rm{Piwi}},8{\rm{-mer}}}$$\end{document} = 160 ± 80 pM (Fig. 1b). The intracellular concentration of the most abundant piRNAs (19 nM) is less that of than the most abundant miRNAs in mouse primary spermatocytes (24 nM)^21^, which suggests that the weaker affinity of PIWI proteins for the canonical seed sites will result in lower occupancy of such targets. Our data therefore disfavour a model in which PIWI proteins find and productively regulate targets through seven-nucleotide canonical seed pairing.

For AGO proteins, extending pairing beyond the canonical seed did not increase the affinity of RISC for a target. For piRNA-1, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{AGO}}2,{\rm{g}}2{\rm{-}}{\rm{g}}10}$$\end{document} = 10 ± 5 pM compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{AGO}}2,8{\rm{-mer}}}$$\end{document} = 10 ± 8 pM (Fig. 1b,c and Extended Data Fig. 2a). By contrast, extending guide:target complementarity from g2–g8 to g2–g10 increased MILI and MIWI target affinity by 4–6-fold. For piRNA-1, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MILI}},{\rm{g}}2{\rm{-}}{\rm{g}}10}$$\end{document} = 7 ± 3 pM compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MILI}},8{\rm{-mer}}}$$\end{document} = 45 ± 7 pM, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MIWI}},{\rm{g}}2{\rm{-}}{\rm{g}}10}$$\end{document} = 11 ± 5 pM compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MIWI}},8{\rm{-mer}}}$$\end{document} = 40 ± 20 pM (Fig. 1b,c). MILI and MIWI therefore bound to g2–g10 complementary sites with an affinity indistinguishable from that of mouse AGO2 for 8-mer sites (Fig. 1b,c). Compared to MILI and MIWI, EfPiwi required longer guide:target pairing (g2–g12) to achieve the affinity of AGO2 for a canonical 8-mer seed site, with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{Ef{\rm{Piwi}},{\rm{g}}2{\rm{-}}{\rm{g}}12}$$\end{document}  = 10 ± 4 pM (Fig. 1b,c). For both a piRNA of synthetic sequence (piRNA-1) and a piRNA found in vivo (L1MC piRNA, antisense to the mouse L1MC retrotransposon), the binding affinity of PIWI proteins increased linearly with increasing predicted base pairing energy (ΔG^0^; Fig. 1c). Thus, MILI, MIWI and EfPiwi^33^ require a longer extent of guide:target base pairing to approach the affinity of AGO2 RISC for a canonical seed match.

Unlike AGO proteins, PIWI proteins bound with similar affinity to sites both with and without full pairing to the canonical seed (Fig. 1b,d and Extended Data Fig. 2b). For instance, MILI and MIWI bound to targets with nine-nucleotide uninterrupted complementarity to the guide starting at g2 or g5 with nearly indistinguishable affinities. For piRNA-1, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MILI}},{\rm{g}}2{\rm{-}}{\rm{g}}10}$$\end{document} = 7 ± 3 pM compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MILI}},{\rm{g}}5{\rm{-}}{\rm{g}}13}$$\end{document} = 12 ± 4 pM, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MIWI}},{\rm{g}}2{\rm{-}}{\rm{g}}10}$$\end{document} = 11 ± 5 pM compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MIWI}},{\rm{g}}5{\rm{-}}{\rm{g}}13}$$\end{document} = 17 ± 7 pM (Fig. 1b,d). Even when guide:target pairing did not start until g8, MILI and MIWI binding was only fivefold weaker compared with targets paired from g2. For piRNA-1, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MIWI}},{\rm{g}}2{\rm{-}}{\rm{g}}10}$$\end{document}  = 11 ± 5 pM compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MIWI}},{\rm{g}}8{\rm{-}}{\rm{g}}16}$$\end{document}  = 60 ± 20 pM (Fig. 1b,d). These RNA Bind-’n-Seq data agreed well with direct measurement of binding affinity for individual RNA targets. For piRNA-1, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MIWI}},{\rm{g}}2{\rm{-}}{\rm{g}}10}$$\end{document}  = 14 ± 9 pM compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{MIWI}},{\rm{g}}8{\rm{-}}{\rm{g}}16}$$\end{document}  = 80 ± 20 pM (Extended Data Fig. 2c).

EfPiwi also bound to nine-nucleotide sites starting at g6, g7 or g8 only 3–4-fold less tightly than when complementarity started at g2. For piRNA-1, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{Ef{\rm{Piwi}},{\rm{g}}2{\rm{-}}{\rm{g}}10}$$\end{document}  = 70 ± 30 pM compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{Ef{\rm{Piwi}},{\rm{g}}8{\rm{-}}{\rm{g}}16}$$\end{document}  = 200 ± 100 pM (Fig. 1b,d). By contrast, AGO2 bound sites lacking seed pairing 10–100-fold more weakly than those containing a seed match (Fig. 1b,d and Extended Data Fig. 2b). For piRNA-1, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{AGO}}2,{\rm{g}}2{\rm{-}}{\rm{g}}10}$$\end{document}  = 10 ± 5 pM compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{{\rm{d}}}^{{\rm{AGO}}2,{\rm{g}}8{\rm{-}}{\rm{g}}16}$$\end{document}  = 500 ± 100 pM. Compared to AGO2, EfPiwi had 3–15-fold higher affinity for longer complementary sites (≥11 nucleotides), with pairing starting at g2, g3, g4, g5 and g6 (Extended Data Fig. 2d). We conclude that PIWI proteins are more flexible than AGO proteins in the types of sites they can bind but require longer complementarity for high-affinity binding.

Slicing of partially complementary RNA

The modal length of piRNAs is 5–10 nucleotides longer than that of siRNAs (26–31 nucleotides compared with 21 nucleotides)^1^, yet MILI, MIWI and EfPiwi do not require pairing to these additional 3′ nucleotides to cleave a target RNA^15,27,33^. In vitro, 16–23-nucleotide-long contiguous complementarity was sufficient for MILI, MIWI and EfPiwi to reach their maximum endonuclease rate (Extended Data Fig. 3a). Moreover, MILI and MIWI, directed by 21-nucleotide-long guides, cleave targets as efficiently as when loaded with full-length 26-nucleotide or 30-nucleotide piRNAs^27^. Nonetheless, we found that extending pairing beyond piRNA nucleotide g20 enabled MILI and MIWI to tolerate guide:target mismatches.

We used the high-throughput Cleave-’n-Seq approach^3^ to determine the rates of cleavage for thousands of target variants (Extended Data Fig. 3b). We incubated purified MILI, MIWI or EfPiwi piRISC complexes (1 nM) and the PIWI auxiliary factor GTSF1 (500 nM) with a library containing 7,700–10,400 30-nucleotide-long target RNAs for different lengths of time (60 s to 16 h). Uncleaved RNAs were reverse-transcribed and sequenced, and their abundance at each time point was used to determine their pre-steady-state cleavage rate, k (Extended Data Fig. 3b and Supplementary Table 1).

Consistent with the idea that additional complementarity to piRNA 3′ nucleotides accelerates cleavage by MILI or MIWI of imperfectly paired targets, a mismatch between g2 and g20 decreased the median k by 3.6-fold when all nucleotides after g20 were unpaired, but by only 1.4-fold when g21–g25 were also base paired (Fig. 2a and Extended Data Fig. 3c). Two mismatches between g2 and g20 caused a 30-fold median reduction in k for targets with no pairing beyond g20, but a reduction of only 3.4-fold when g21–g25 were paired (Fig. 2a and Extended Data Fig. 3c). Thus, endonucleolytic cleavage by MILI or MIWI does not require target pairing to piRNA 3′ sequences, but such extended complementarity readily compensates for guide:target mismatches. We did not observe the same compensatory effect for EfPiwi. Compared to MILI and MIWI, slicing by EfPiwi was generally slower (Extended Data Fig. 3a), probably because Gtsf1 from Ephydatia muelleri (EmGtsf1) was used to stimulate EfPiwi-catalysed target cleavage (Methods).Fig. 2PIWI slicing tolerates mismatches with any target nucleotide.a, Change in pre-steady-state cleavage rate for one or two mismatches (pink) between g2 and g20. For one mismatch, n = 456: all 19 possible positions × 3 geometries × 4 piRNAs × MILI and MIWI. For two mismatches, n = 1,368: all 171 possible combinations × 1 geometry × 4 piRNAs × MILI and MIWI. Box plots show the IQR and median. Statistical analysesare in Extended Data Fig. 3c. b, Change in pre-steady-state cleavage rate for one or two consecutive mismatches between g2 and g20 for contiguous g2–g21 or g2–g25 pairing for MILI, MIWI, EfPiwi and mouse AGO2. Median and IQR are shown. For one mismatch, n = 24 (3 geometries × 4 piRNAs × MILI and MIWI); n = 6 for EfPiwi (3 geometries × 2 piRNAs); n = 21 for AGO2 (3 geometries for L1MC guide and 3 geometries × 3 contexts for let-7a and miR-21 guides). For two consecutive mismatches, n = 8 (1 geometry × 4 piRNAs × MILI and MIWI); n = 2 for EfPiwi (1 geometry × 2 piRNAs); n = 19 for AGO2 (1 geometry for L1MC RISC and 9 geometries for let-7a and miR-21 RISCs). All data and statistical analyses are in Extended Data Fig. 4. ND, not detected. c, MIWI, MILI and EfPiwi pre-steady-state cleavage rates (k) for targets of L1MC piRNA containing a single unpaired nucleotide. Position and identity of mononucleotide mismatch in targets (indicated in blue) of L1MC piRNA (indicated in red) are on the top of the chart.

AGO-clade Argonaute proteins secure target nucleotide t1 in a pocket that often displays specific nucleotide preferences. Although MILI, MIWI and EfPiwi showed a slight binding preference for t1A and t1U (≤3-fold stronger affinity compared with t1S; Extended Data Fig. 3d), target cleavage showed no t1 preference (Extended Data Fig. 3a).

PIWI proteins tolerate mismatch at any position

Efficient target RNA cleavage by animal and plant AGO proteins requires uninterrupted base pairing to siRNA nucleotides g9–g13 (refs. ^3,4,34^). By contrast, MILI, MIWI and EfPiwi slicing tolerated mismatches at any position within the region of complementarity (Fig. 2b,c, Extended Data Fig. 4 and Supplementary Figs. 3 and 4).

For g2–g21-paired targets of AGO2, a single mononucleotide mismatch at g9, g10, g11 or g13 decreased the median k by 10–200-fold^3^ (Fig. 2b and Extended Data Fig. 4). For the same extent of pairing, the median reduction in k was ≤5-fold for MILI and MIWI and ≤7-fold for EfPiwi for a mononucleotide mismatch at any position between g2 and g20 (Fig. 2b and Extended Data Fig. 4). For MILI and MIWI, guide:target complementarity from g2 to g25 reduced the median effect of a mononucleotide mismatch to ≤3-fold at any position between g2 and g20 (Fig. 2b and Extended Data Fig. 4). Among mismatch types, GU wobbles had the smallest impact on endonucleolytic rate, decreasing k by 1.2-fold (inter-quartile range (IQR) of 1–1.65; Extended Data Fig. 5a).

Unlike AGO2 RISC, pairing to target nucleotides adjacent to the scissile phosphate was dispensable for target slicing by piRISC (Fig. 2b,c, Extended Data Fig. 4 and Supplementary Figs. 3 and 4). Loss of pairing to either t10 or t11 in a g2–g21 match decreased the median cleavage rate by 4-fold for MILI and MIWI (IQR of 2.5–19 for t10, and IQR of 2.8–5.7 for t11), and by 5–6.4-fold for EfPiwi (IQR of 4–18.9 for t10, and IQR of 3.2–5.3 for t11; Fig. 2b and Extended Data Fig. 4). Purine–purine mismatches at these positions appeared to be the least tolerated by piRISC, perhaps because their greater bulk is poorly accommodated in the PIWI catalytic centre (Extended Data Fig. 5b). piRISC-catalysed slicing was detectable even when both t10 and t11 were unpaired. The median decrease in k was 60-fold for MILI and MIWI (IQR of 40–70) and 67-fold for EfPiwi (IQR of 54–88; Fig. 2b and Extended Data Fig. 4). For the same t10–t11 dinucleotide mismatch, target cleavage was undetectable with AGO2 RISC (Fig. 2b and Extended Data Fig. 4). We conclude that MILI, MIWI and EfPiwi, unlike AGO2, can efficiently cleave partially paired RNAs with mismatches anywhere in a target site.

Sequencing the 3′ products of piRNA-directed, MIWI-catalysed slicing showed that piRISC invariably hydrolysed the RNA at the canonical scissile phosphodiester bond, between target nucleotides t10 and t11, even when both g10 and g11 were unpaired or when contiguous pairing did not start until g11 (Extended Data Fig. 5c). These data suggest that piRNA–target base pairing near the cleavage site has little if any role in positioning the scissile phosphate within the MIWI catalytic centre.

Mismatch tolerance is intrinsic to PIWI proteins

Unlike AGO-clade Argonaute proteins, PIWI proteins require the auxiliary factor GTSF1 to achieve their maximal catalytic rate^27^ (Extended Data Fig. 6a). Without GTSF1, PIWI Argonaute proteins inefficiently slice even perfectly complementary RNAs^27^.

Although GTSF1 potentiates target cleavage by PIWI proteins, it is not required for PIWI tolerance of guide:target mismatches. GTSF1 accelerated slicing of fully and partially complementary RNAs to a similar extent (Extended Data Fig. 6b). In the presence of GTSF1, the median increase in MILI cleavage rate was 25-fold (95% confidence interval (CI) of [19, 29]) for perfectly complementary RNAs and 17–50-fold for mismatched targets. GTSF1 enhanced MIWI-catalysed slicing of fully (median increase, 14-fold; CI of [12.8, 15.4]) and partially (11–35-fold) complementary targets with similar efficacy. EmGtsf1 also accelerated EfPiwi slicing for perfectly paired (15-fold, CI of [8, 28]) and mismatched RNAs (4–19-fold) to a comparable degree (Extended Data Fig. 6b). These data suggest that, unlike AGOs, PIWI inherently accommodates unpaired target nucleotides and that GTSF1 only accelerates cleavage.

Relaxed rules for slicing apply in vivo

In mice, piRNAs direct MILI and MIWI to slice complementary transposon transcripts, mRNAs or long non-coding transcripts (lncRNAs)^7,8,10,11,15,16^. As we observed for purified piRISC, piRNAs directed MILI and MIWI to cleave targets with as few as 15–19-nucleotide complementary nucleotides in vivo in mouse primary spermatocytes.

Mouse primary spermatocytes produce a class of MILI-loaded and MIWI-loaded piRNAs called pachytene piRNAs, which first appear at the pachytene stage of meiosis^1^. Mouse GTSF1 is present in all meiotic male germ cells^27^. Because endonucleolytic cleavage by Argonaute proteins leaves a 5′-monophosphate^35,36^, we sequenced 5′-monophosphorylated RNAs from mouse primary spermatocytes purified by fluorescence activated cell sorting (FACS) to identify potential 3′ cleavage products generated in vivo by MILI and MIWI (Fig. 3a). Restricting our analysis to pachytene piRNAs, >80% of which derive from non-repetitive sequences, ensured unambiguous assignment of piRNAs to candidate cleavage products. To identify those RNAs corresponding to 3′ cleavage products generated by piRNA-directed slicing, we searched for 5′-monophosphate-bearing RNAs present in control C57BL/6 mice, but the abundance of which was reduced by ≥8-fold in a triple mutant lacking all piRNAs from three major pachytene piRNA-producing loci on chromosomes 2, 9 and 17: 2-qE1-35981(+); 9-qC-31469(−),10667(+); 17-qA3.3-27363(−),26735(+) (refs. ^10,37^). For simplicity, we refer to these loci as pi2, pi9 and pi17, respectively. The pi2^−/−^pi9^−/−^pi17^−/−^ triple mutation removes approximately 22% of all pachytene piRNAs (Fig. 3a and Extended Data Fig. 7).Fig. 3. Mouse PIWI proteins cleave partially complementary targets in vivo.a, Schematic of the strategy used to identify 3′ cleavage products of piRNA-guided PIWI-catalysed slicing and to measure the fraction of targets cleaved by PIWI proteins in FACS-purified mouse primary spermatocytes. b, Fraction of cleaved MILI and MIWI targets in FACS-purified mouse primary spermatocytes for contiguous pairing from nucleotide g2. c, Fraction of cleaved targets in FACS-purified mouse primary spermatocytes for perfect matches (indicated in blue) and for pairing containing a single-nucleotide mismatch (indicated in pink). Horizontal dotted lines indicate the medians for perfect matches. d, MIWI, MILI, and EfPiwi pre-steady-state cleavage rates in vitro for all possible stretches of ≥6-nucleotide contiguous pairing starting from nucleotides g2–g15 of L1MC piRNA. e, Fraction of cleaved targets in FACS-purified mouse primary spermatocytes for 14-nucleotide contiguous pairing starting from nucleotides g2 to g5. Data are binned by piRNA intracellular concentration (<30, 30–50, 50–100, 100–500, >500 pM). For b,c and e, box plots show IQR and median; 95% CI was calculated with 10,000 bootstrapping iterations; n = 16 permutations of 4 control (C57BL/6) and 4 pi2^−/−^pi9^−/−^pi17^−/−^ animals.

Among the 5′-monophosphorylated RNAs detected in the control C57BL/6 mice, we selected candidate cleavage products for which production could be explained by a pi2, pi9 or pi17 piRNA directing cleavage between nucleotides t10 and t11. For each pattern of piRNA–target complementarity, for example, g2–gX pairing with t2–tX, we calculated the fraction of cleavage product candidates for which abundance was reduced by ≥8-fold in pi2^−/−^pi9^−/−^pi17^−/−^ mutant primary spermatocytes. We denote this fraction as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${f}_{{\rm{d}}{\rm{e}}{\rm{c}}{\rm{r}}{\rm{e}}{\rm{a}}{\rm{s}}{\rm{e}}{\rm{d}}\,{\rm{i}}{\rm{n}}\,{\rm{m}}{\rm{u}}{\rm{t}}{\rm{a}}{\rm{n}}{\rm{t}}}^{{p}{i}{2}{,}{p}{i}{9}{,}{p}{i}{17}}$$\end{document} . A pairing configuration that can support target slicing is predicted to have a high fraction of such candidate 3′ cleavage products present in the control mice but reduced or absent in the triple mutant mice (Fig. 3a).

To account for sampling error arising from the short-lived nature of 5′-monophosphorylated fragments in vivo, we identified all 5′-monophosphorylated RNAs in C57BL/6 explained by piRNAs not removed in pi2^−/−^pi9^−/−^pi17^−/−^ mice (control piRNAs), and then calculated the fraction of these RNAs reduced by ≥8-fold in pi2^−/−^pi9^−/−^pi17^−/−^ animals. The fraction of cleaved targets for each pairing arrangement g2–gX was then calculated as the observed signal minus the sampling error (Fig. 3a) as follows:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${f}_{{\rm{c}}{\rm{l}}{\rm{e}}{\rm{a}}{\rm{v}}{\rm{e}}{\rm{d}}}^{{\rm{g}}2-{\rm{g}}X}={f}_{{\rm{d}}{\rm{e}}{\rm{c}}{\rm{r}}{\rm{e}}{\rm{a}}{\rm{s}}{\rm{e}}{\rm{d}}\,{\rm{i}}{\rm{n}}\,{\rm{m}}{\rm{u}}{\rm{t}}{\rm{a}}{\rm{n}}{\rm{t}}}^{{p}{i}2,{p}{i}9,{p}{i}17}-{f}_{{\rm{d}}{\rm{e}}{\rm{c}}{\rm{r}}{\rm{e}}{\rm{a}}{\rm{s}}{\rm{e}}{\rm{d}}\,{\rm{i}}{\rm{n}}\,{\rm{m}}{\rm{u}}{\rm{t}}{\rm{a}}{\rm{n}}{\rm{t}}}^{{\rm{c}}{\rm{o}}{\rm{n}}{\rm{t}}{\rm{r}}{\rm{o}}{\rm{l}}}$$\end{document}

Cleavage by MILI or MIWI is indistinguishable in our data, thus \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${f}_{{\rm{c}}{\rm{l}}{\rm{e}}{\rm{a}}{\rm{v}}{\rm{e}}{\rm{d}}}^{{\rm{g}}2-{\rm{g}}\,X}$$\end{document} corresponds to the sum of targets sliced in mouse primary spermatocytes by both PIWI proteins.

These analyses showed that PIWI-catalysed cleavage was detected in primary spermatocytes for piRNA–target base pairing as short as g2–g16, with median \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${f}_{{\rm{cleaved}}}^{{\rm{g}}2{\rm{-}}{\rm{g}}16}$$\end{document}  = 0.18, CI of [0.08, 0.23] (Fig. 3b, Supplementary Table 2 and Supplementary Data 1). The efficiency of piRNA-directed target cleavage increased with longer complementarity. For example, median \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${f}_{{\rm{cleaved}}}^{{\rm{g}}2{\rm{-}}{\rm{g}}20}$$\end{document}  = 0.52 (CI of [0.37, 0.65]; Fig. 3b and Supplementary Table 2). Note that pairing longer than g2–g20 contained too few data points to measure the corresponding \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${f}_{{\rm{c}}{\rm{l}}{\rm{e}}{\rm{a}}{\rm{v}}{\rm{e}}{\rm{d}}}^{{\rm{g}}2-{\rm{g}}\,X}$$\end{document} .

When the piRNA–target complementarity was ≥17 nucleotides long, single mismatches were tolerated at most positions (Fig. 3c, Extended Data Fig. 8a, Supplementary Table 2 and Supplementary Data 1). For example, for g2–g18 complementarity, the median fcleaved for perfect pairing (0.18, CI of [0.14, 0.23]) was similar to g2–g18 matches bearing a single nucleotide mismatch at positions g2, g5, g6, g9–g14 or g17–g18 (0.15–0.25; Fig. 3c, Supplementary Table 2 and Supplementary Data 1). For targets complementary to piRNA nucleotides g2–g18 or g2–g19, pairing to t10 and t11, the target nucleotides flanking the scissile phosphate, was dispensable for slicing (Fig. 3c, Supplementary Table 2 and Supplementary Data 1). The lower median fcleaved values for targets mismatched to piRNA 5′ sequences compared with other piRNA regions may reflect slower on-rates for piRNAs of low intracellular concentration (see also the next section). For g2–g18 targets, the median \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${f}_{{\rm{cleaved}}}^{{\rm{mismatches}}\,{\rm{to}}\,{\rm{g3}}{\rm{-g8}}}$$\end{document}  = 0.10 compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${f}_{{\rm{cleaved}}}^{{\rm{mismatches}}\,{\rm{to}}\,{\rm{g9}}{\rm{-g18}}}$$\end{document}  = 0.18, whereas for g2–g19 targets, the median \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${f}_{{\rm{cleaved}}}^{{\rm{mismatches}}\,{\rm{to}}\,{\rm{g3}}{\rm{-g8}}}$$\end{document}  = 0.06 compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${f}_{{\rm{cleaved}}}^{{\rm{mismatches}}\,{\rm{to}}\,{\rm{g9}}{\rm{-g19}}}$$\end{document}  = 0.21.

Together, these in vivo data corroborate our in vitro target cleavage experiments, which found that the median decrease in cleavage rate was fourfold or less for a mononucleotide mismatch at any position between g2 and g20 (Fig. 2b and Extended Data Fig. 4). Thus, a wide variety of piRNA–target pairing patterns can efficiently direct MILI and MIWI to cleave targets, unlike the relatively limited pairing configurations tolerated by AGO-clade Argonaute proteins.

Abundant piRNAs slice without seed match

AGO2-catalysed slicing of sites contiguously paired from g4 or g5 (that is, without canonical seed pairing) has been observed in vitro but not detected in vivo^3,38^. By contrast, our data identified piRNA-directed cleavage in vivo in mouse primary spermatocytes of targets for which pairing to the guide starts at g3, g4 or g5.

In vitro, MILI and MIWI did not require target pairing to piRNA 5′ terminal nucleotides for binding (Fig. 1b,d) or slicing (Fig. 3d and Extended Data Fig. 8b). Similarly, we detected in vivo piRNA-directed cleavage of targets that lack complementarity to nucleotides g2–g4 (Fig. 3e, Supplementary Table 2 and Supplementary Data 1). Pre-organization of seed nucleotides g2–g6 in an A-form-like helix accelerates target finding by AGO proteins (for the let-7a 8-mer, target-finding rate constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_{{\rm{on}}}^{{\rm{guide}}\,{\rm{alone}}}$$\end{document} of 5 × 10^6^ M^−1 ^s^−1^ compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_{{\rm{on}}}^{{\rm{guide}}\,{\rm{in}}\,{\rm{AGO}}2}$$\end{document} of 2.4 × 10^8^ M^−1^ s^−1^)^39^. By contrast, PIWI proteins pre-organize only nucleotides g2–g4 (refs. ^33,40,41^), which suggests that there is slower target finding when piRNA–target complementarity begins after nucleotide g4. In theory, high piRISC concentration could compensate for a slower kon value. In vivo, piRNA concentrations vary widely, and about 1,500 piRNAs are present in mouse primary spermatocytes at ≥500 pM (Extended Data Fig. 8c). We observed piRNA-directed cleavage of targets with contiguous 14-nucleotide pairing beginning at g4 (median \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${f}_{{\rm{cleaved}}}^{{\rm{g}}4{\rm{-}}{\rm{g}}17}$$\end{document}  = 0.43, CI of [0.29–0.56]) or g5 (median \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${f}_{{\rm{cleaved}}}^{{\rm{g}}5{\rm{-}}{\rm{g}}18}$$\end{document}  = 0.33, CI of [0.2–0.4]) only for highly abundant piRNAs (≥500 pM). Cleavage of targets paired from g3 was detectable for piRNAs for which the in vivo concentration was ≥100 pM, and targets paired from g2 were sliced by piRNAs present at ≥50 pM (Fig. 3e, Supplementary Table 2 and Supplementary Data 1). These in vivo data were consistent with equilibrium binding measurements showing that the affinity (Kd) of MILI or MIWI for target sites with pairing starting at g2, g3, g4 or g5 was <500 pM (Fig. 1d).

We conclude that because pairing to piRNA 5′ terminal nucleotides is dispensable for both target finding and slicing, piRISC efficiently cleaves targets that lack full complementarity to the canonical 5′ seed.

Insertions and deletions thwart piRNA-directed slicing

As observed for AGO2 (ref. ^3^), mononucleotide target insertions between t9 and t15 slowed cleavage catalysed by MILI, MIWI or EfPiwi in vitro by ≥10-fold (Extended Data Fig. 9a). For both AGO2 (refs. ^3,42^) and PIWI proteins^33,40,41^, target nucleotides t9–t15 face the protein surface, which makes insertions likely to distort the catalytic centre.

Single-nucleotide target deletions between t6 and t15 were also poorly tolerated (≥5-fold lower k in vitro for MILI, MIWI and EfPiwi relative to a fully complementary target; Extended Data Fig. 9b). Such target sequence deletions result in mononucleotide bulges in the piRNA guide. Similar to mammalian AGO2 (ref. ^42^), EfPiwi restricts piRNA nucleotides g6–g10 to the central cleft of the protein^33^, and a mononucleotide piRNA bulge between t6 and t10 is unlikely to fit in this narrow furrow, which potentially explains why PIWI proteins do not tolerate such target deletions. By contrast, single-nucleotide deletions between t11 and t15 create solvent-facing, mononucleotide loops of guide nucleotides g11–g15 that are predicted to be accommodated. Indeed, AGO2 tolerates target deletions between t11 and t15 (ref. ^3^). Notably, PIWI proteins did not tolerate t11–t15 target deletions (Extended Data Fig. 9b). We speculate that piRNA guide bulges between t11 and t15 specifically affect PIWI proteins because they impair interactions with GTSF1.

Target insertions or deletions near the centre of the piRNA–target duplex were also not tolerated in vivo in mouse primary spermatocytes (Extended Data Fig. 9c, Supplementary Table 2 and Supplementary Data 1). We note that insertions and deletions occur in mammalian genomes 30-fold less frequently than single-nucleotide polymorphisms^43^.

Predicting in vivo target cleavage sites

We used a logistic regression classifier approach to identify factors that predict effective piRNA slicing in vivo. Cleavage data from pachytene spermatocytes were used to fit a logistic function representing the probability of piRNA-guided cleavage, P(cleaved), determined by 35 variables (x1, x2, … x35): the presence or absence of pairing with each guide nucleotide between g2 and g25; the total number of paired nucleotides; the predicted binding energy; the piRNA abundance; the target site location in the transcript (5′ untranslated region (UTR), open reading frame (ORF), 3′ UTR or lncRNA; and the identity of target nucleotide t1 (Fig. 4a). The coefficient for each variable in the fitted logistic decision function (β1, β2, … β35) estimates the importance of each feature as follows:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P({\rm{cleaved}})=\frac{1}{1+{e}^{-({\beta }_{0}+{\beta }_{1}{x}_{1}+{\beta }_{2}{x}_{2}+\ldots +{\beta }_{35}{x}_{35})}}$$\end{document}

Fig. 4. Determinants of PIWI slicing in vivo.a, Decision function coefficients for 400 logistic function fits (regression models) using around 3,500 distinct piRNA–target pairs detected in mouse primary spermatocytes. n = 16 permutations of 4 control (C57BL/6) and 4 pi2^−/−^pi9^−/−^pi17^−/−^ animals × 5-repeated × 5-fold cross validation. Box plots show IQR and median. b, Number of piRNAs and siRNAs predicted to cleave mutated versions of the L1Md AI transposon sequence. Data are median and IQR from 100 independent simulations. c, AGO and PIWI proteins use different rules to find and slice targets.

We selected about 3,500 distinct pairs of pi2, pi9 and p17 piRNAs and target sites for which 5′-monophosphorylated 3′ cleavage products both were detected in the control C57BL/6 mice and had ≥19 nucleotides paired between g2 and g25. Target sites were considered cleaved if the abundance of the 3′ cleavage products decreased by ≥8-fold in the pi2^−/−^pi9^−/−^pi17^−/−^ mutant mice compared with control mice; all other target sites were assigned as not cleaved. To test whether the logistic regression models created with the pi2, pi9 and pi17 piRNA data could predict cleavage by non-pi2, non-pi9 and non-pi17 piRNAs, we generated independent datasets from mice with a mutation disrupting a pachytene piRNA-producing locus on chromosome 7, pi7 (7-qD2-24830(−)11976(+); Extended Data Fig. 7). Note that >99% of piRNAs eliminated in pi7^−/−^ mice are not found in pi2, pi9 or pi17. The performance of the logistic function fitted to pi2^−/−^pi9^−/−^pi17^−/−^ data was similar when tested with either pi2^−/−^pi9^−/−^pi17^−/−^ data (Extended Data Fig. 10a) or pi7^−/−^ data (Extended Data Fig. 10a).

The most predictive features, that is, highest median coefficients in the logistic decision function, were in vivo piRNA concentration (+2.11), predicted energy of piRNA–target base pairing (+1.82) and the total number of paired nucleotides (+1.20; Fig. 4a). These results show that in vivo, piRISC behaves as a conventional enzyme. That is, its concentration and substrate-binding strength determine the efficacy of target cleavage. A high aggregate number of targeting piRNAs was also recently shown to be required for potent transcriptional silencing^44^.

Consistent with the idea that PIWI slicing does not rely on complementarity to specific target nucleotides (Figs. 2b and 3c), the median decision function coefficients were ≤+0.6 for guide:target pairing at any individual position (Fig. 4a). Extensive complementarity anywhere in the piRNA 5′ half seems to initiate target binding (Fig. 1b,d), and highly abundant piRNAs (≥500 pM) even direct slicing of targets without pairing to positions g2–g4 (Fig. 3e). In agreement with these data, logistic function coefficients for matches to g2–g10 were higher than those for pairing to other nucleotides ([+0.25, +0.6] compared with [−0.12, +0.25]; Fig. 4a).

The features with the lowest median coefficients were the identity of nucleotide t1 and the location of the target site within a transcript (−0.18 to +0.13; Fig. 4a). This result suggests that these factors are not rate-determining for piRNA-guided cleavage in vivo.

Together, these analyses show that simple biochemical principles are sufficient to predict efficient piRNA-directed cleavage in vivo. First, piRNA concentration determines how frequently a target encounters piRISC and therefore the concentration of the piRISC–target complex. Second, tighter guide:target base pairing (binding energy) extends the lifetime of the piRISC–target complex, which increases the likelihood of cleavage.

Implications of piRNA targeting rules

Because PIWI-catalysed slicing does not depend on pairing to a specific piRNA nucleotide position, target mutations are predicted to be better tolerated by PIWI proteins than by AGO proteins. Our computational simulations estimated that when a transposon sequence mutates but the small RNA repertoire does not change, the number of guides capable of directing target slicing decreases by fourfold more slowly for PIWI than AGO proteins (Fig. 4b and Extended Data Fig. 10b). We simulated 1,000 rounds of single-nucleotide substitutions in the LINE1 consensus sequences, excluding non-synonymous mutations in ORFs. At each round, we recorded the number of embryonic testicular piRNAs or siRNAs (simulated using piRNA 21 nucleotide prefixes) expected to productively slice the mutated transposons. The number of MILI-loaded piRNAs capable of cleavage decreased at 0.01% of guides per single-nucleotide substitution in LINE1 elements, whereas the number of simulated AGO2-loaded siRNAs decreased at 0.04% of guides per mutation in the transposon sequence (Fig. 4b and Extended Data Fig. 10b).

Discussion

Our data highlight several distinct features of PIWI proteins that set them apart from AGO-clade Argonaute proteins. First, PIWI proteins do not limit target finding to the seven-nucleotide, 5′ canonical seed. The central cleft that cradles the guide RNA is wider in PIWI than in AGO proteins^33,40,41^, which perhaps enables PIWI proteins to productively use piRNA nucleotides 3′ to g8 to initiate pairing with targets. At high concentrations, piRISC efficiently binds and slices RNAs unpaired to nucleotides g2–g4, so targeting capacity is similar for piRNAs for which 5′ ends are several nucleotides apart along a piRNA precursor transcript. This observation explains why piRNAs can tolerate a high degree of 5′ heterogeneity, which is an intrinsic feature of phased (trailing) piRNA biogenesis^17–21^. By contrast, miRNA 5′-isoforms have distinct target repertoires^6^, and any change in the 5′ position of a siRNA duplex can invert which strand becomes a guide for an AGO protein.

Second, piRNA-directed slicing tolerates mismatches to any nucleotide of the piRNA guide. Despite more than 900 million years of independent evolution, both mammalian and poriferan PIWI proteins can slice targets even with mismatches adjacent to the scissile bond. In fact, both a perfectly paired and a t10 or t11 mismatched target site were cleaved with similar efficiency when introduced into the 3′ UTR of mouse Ythdc2 mRNA^45^. Mismatches near the scissile phosphate bond not only block target cleavage by AGO proteins but also promote AGO protein degradation^46,47^. We propose that because piRNAs are generally longer than siRNAs, PIWI proteins can extend the lifetime of the piRISC–target complex through compensatory pairing to piRNA 3′ sequences, which enables piRISC to tolerate multiple guide:target mismatches. Consistent with this model, piRNA 3′ ends are 2′-O-methylated to protect them from decay triggered by extensive pairing to targets^48^. This hypothesis also predicts that cleavage by the ovary-specific mammalian PIWIL3 protein, which uses unusually short, 19-nucleotide-long piRNA guides^49,50^, will require full complementarity between their guides and targets.

The catalytic centres of AGO-clade and PIWI-clade Argonaute proteins probably have distinct requirements for effective catalysis of target cleavage. For AGO proteins, perfect pairing between a siRNA and its target moves target nucleotides t10 and t11 into the endonuclease active site^4^. Our data suggest that the catalytically competent geometry of PIWI proteins does not intrinsically rely on perfect complementarity between the target and piRNA near the cleavage site.

Third, the ability of piRNAs to direct cleavage of imperfectly complementary RNAs, but only when the extent of complementarity is more than 15 contiguous base pairs, may explain the rapid evolution of the piRNA target repertoire (Fig. 3b and Extended Data Fig. 3a). Our analyses suggest that a 15-nucleotide contiguous match is sufficient to prevent inappropriate targeting of mRNAs and other self-transcripts. We determined the fraction of all mouse mRNAs and lncRNAs that contain at least one k-mer from the consensus sequences of transpositionally active families of mouse LTR or LINE transposons. For k ≥ 15, less than 5% of mRNAs and lncRNAs shared at least one *k-*mer with transposon sequence (Extended Data Fig. 10c, top). We obtained similar results when these analyses were conducted for intron sequences from the same mouse transcripts (Extended Data Fig. 10c, bottom). This result suggests that strong negative selection against transposon-derived ≥15-mers in mRNAs and lncRNAs is unlikely. Together, our findings (Fig. 4c) provide a plausible explanation for why the piRNA pathway, not RNA interference, was favoured by evolution as the mainstay of transposon defence in animals.

Methods

Mouse strains and mutants

Mice (Supplementary Table 3; C57BL/6J, International Mouse Strain Resource JAX:000664) were housed in an Association for Assessment and Accreditation of Laboratory Animal Care International -accredited barrier facility at controlled temperature (22 ± 2 °C), relative humidity (40 ± 15%) and a 12-h day–light cycle. All experimental animals were 2–6 months old. All procedures were reviewed and performed in compliance with the guidelines of the Institutional Animal Care and Use Committee (IACUC) of the University of Massachusetts Chan Medical School (IACUC protocol number A201900331). Single-guide RNAs (sgRNAs; Supplementary Table 3) were designed using a CRISPR design tool (https://portals.broadinstitute.org/gppx/crispick/public). sgRNAs were transcribed with T7 RNA polymerase and then purified by electrophoresis on a 10% denaturing polyacrylamide gel. gRNA (20 ng µl^–1^) and Cas9 mRNA (50 ng µl^–1^, TriLink Biotechnologies, L-7206) were injected together into the pronucleus of one-cell C57BL/6 zygotes in M2 medium (Sigma, M7167). After injection, the zygotes were cultured in EmbryoMax Advanced KSOM medium (Sigma, MR-106-D) at 37 °C under 5% CO_2_ until the blastocyst stage (3.5 days), then transferred into the uterus of pseudopregnant ICR females 2.5 days post coitum. To screen for mutant founders, gDNA extracted from tail tissues was analysed by PCR using the primers listed in Supplementary Table 3.

No statistical method was used to determine the sample size. For biological samples, the maximum possible sample size (n = 4–12) was used for each type of data, which ensured that variability arising from all accountable sources was incorporated in the analyses (animal, day of data collection, reagent lots). No data were excluded from the analyses. Randomization is not relevant to this study because it did not involve treatment or exposure of animals to any agent. Instead, untreated wild-type mice were compared with untreated mutant mice lacking piRNAs from four genomic loci. Blinding is not relevant to this study because during analyses, wild-type control and mutant datasets were easily identified. Blinding was not performed during data acquisition and/or analysis.

piRNA loading and recombinant piRISC purification for MILI and MIWI

Synthetic piRNA guides (IDT) were purified by electrophoresis through a 15% denaturing polyacrylamide gel. HEK293T cells (American Type Culture Collection) expressing SNAP-tagged, 3×Flag-tagged MILI or MIWI were generated as previously described^27^. Cells were collected at 70% confluency using a TC cell scraper (ThermoFisher, 50809263) into ice-cold PBS and collected by centrifugation at 500g. Supernatant was removed, and the pellet was stored at −80 °C until lysed in 10 ml of 30 mM HEPES-KOH, pH 7.5, 100 mM potassium acetate, 3.5 mM magnesium acetate, 2 mM DTT, 0.1% (v/v) Triton X-100, 15% (v/v) glycerol and 1× protease inhibitor cocktail (1 mM 4-(2-aminoethyl)benzenesulfonyl fluoride hydrochloride (Sigma, A8456), 0.3 μM aprotinin, 40 μM betanin hydrochloride, 10 μM E-64 (Sigma, E3132) and 10 μM leupeptin hemisulfate) per g frozen cells. Cell lysis was monitored by staining with trypan blue. Crude cytoplasmic lysate was clarified at 20,000g, flash frozen in liquid nitrogen and stored at −80 °C.

To capture MILI or MIWI, 1 ml of clarified lysate was incubated with 20 μl anti-Flag M2 paramagnetic beads (Sigma, M8823) for 4 h to overnight rotating at 4 °C. Beads were washed four times with extract buffer (30 mM HEPES-KOH, pH 7.5, 3.5 mM magnesium acetate, 2 mM DTT, 15% (v/v) glycerol and 0.01% (v/v) Triton X-100) containing 2 M potassium acetate and four times with extract buffer containing 100 mM potassium acetate. To assemble MILI or MIWI piRISC, beads were resuspended in extract buffer containing 100 mM potassium acetate and 100 nM synthetic piRNA guide (Supplementary Table 4) and incubated with rotation for 30 min at 37 °C or room temperature. After five washes in 2 M potassium acetate extract buffer and five washes in 100 mM potassium acetate extract buffer, MILI or MIWI piRISC was eluted from the beads twice with 200 ng µl^–1^ 3×Flag peptide in 100 µl of 100 mM potassium extract buffer with rotation for 1 h at room temperature. The combined 200 µl eluate was used immediately for capture oligonucleotide affinity purification or flash frozen in liquid nitrogen and stored at −80 °C.

To purify MILI or MIWI piRISC loaded with a single synthetic RNA, 200 µl Dynabeads MyOne Streptavidin T1 paramagnetic beads (ThermoFisher, 65601) was washed and incubated with 800 pmol 5′ biotinylated, 2′-O-methyl capture oligonucleotide (Supplementary Table 4) according to the manufacturer’s instructions, then resuspended in 100 µl of 100 mM potassium extract buffer. The 200 µl eluate with piRISC from the previous step was added to the capture oligonucleotide-conjugated beads in the extract buffer and incubated with rotation for 1 h at room temperature. The supernatant was removed, and then the beads were washed five times with 100 mM potassium extract, followed by five washes with 2 M potassium extract buffer. piRISC was eluted by rotating the beads for 2 h at room temperature in 200 µl of 100 mM potassium extract buffer containing 1,200 pmol of 5′ biotinylated competitor DNA oligonucleotide (Supplementary Table 4) and S20 testis lysate (total protein 200 ng µl^–1^ final concentration (f.c., see below)). The supernatant containing eluted piRISC was then incubated for 30 min at room temperature with 300 µl Dynabeads MyOne Streptavidin T1 paramagnetic beads (prewashed according to the manufacturer’s instructions followed by two washes in 100 mM potassium extract buffer) to remove excess competitor DNA oligonucleotide. After removing streptavidin beads, 20 μl anti-Flag M2 paramagnetic beads (Sigma, M8823) was added and incubated with the supernatant for 4 h rotating at 4 °C to isolate piRISC from testis lysate. Beads were then washed four times with 2 M potassium acetate extract buffer and four times with 100 mM potassium acetate extract buffer. piRISC was eluted from the beads twice with 200 ng µl^–1^ 3×Flag peptide in 100 µl of 100 mM potassium extract buffer with rotation for 1 h at room temperature. The combined 200 µl eluate was aliquoted, flash frozen in liquid nitrogen and stored at −80 °C.

Testis lysate for eluting MILI and MIWI piRISC from capture oligonucleotide

Dissected animal tissue samples were homogenized at 4 °C in 5 volumes of 30 mM HEPES-KOH, pH 7.5, 100 mM potassium acetate, 3.5 mM magnesium acetate, 1 mM DTT and 15% (v/v) glycerol in a dounce homogenizer using 10 strokes of the loose-fitting pestle A, followed by 20 strokes of tight-fitting pestle B to generate crude lysate. S20 was prepared by clarifying the crude lysate at 20,000g. The protein concentration was estimated using a BCA assay (ThermoFisher, 23200). Crude and fractionated testis lysate were flash frozen in liquid nitrogen and stored at −80 °C.

piRNA loading and recombinant piRISC purification for EfPiwi

Synthetic piRNA guides (IDT) were purified by electrophoresis through a 15% denaturing polyacrylamide gel. EfPiwi protein was expressed as a His_6_-TEV-EfPiwi construct using the Bac-to-Bac Baculovirus Expression System (ThermoFisher, 10359016) and Sf9 cells (American Type Culture Collection). Sf9 infection with EfPiwi-expressing baculovirus^33^ was performed in 750 ml cultures of 1,275 × 10^6^ cells for 72 h at 27 °C. Each 750 ml culture of Sf9 cells was pelleted and resuspended in 25 ml lysis buffer (50 mM Tris pH 8.0, 300 mM NaCl and 0.5 mM TCEP) and lysed using a high-pressure (18,000 p.s.i.) microfluidizer (Microfluidics M100P). Debris was pelleted by centrifugation, and the clarified lysate was incubated with 1 ml Ni-NTA resin (Qiagen, 30210) per 750 ml culture for 1 h at 4 °C, followed by washing twice in nickel wash buffer (50 mM Tris pH 8.0, 300 mM NaCl, 20 mM imidazole and 0.5 mM TCEP). The resin was then washed once with wash buffer supplemented with 5 mM CaCl_2_ in preparation for micrococcal nuclease treatment to degrade co-purifying cellular RNAs. The washed resin was resuspended in nickel wash buffer supplemented with 5 mM CaCl_2_ (final volume of 20 ml). Next, 100 U micrococcal nuclease (Takara Bio, 2910A) was added per 750 ml culture and incubated at room temperature for 1  h, inverting gently every 15 min to resuspend the resin. After three washes with nickel wash buffer without CaCl_2_, protein was eluted with 6× column volumes of nickel elution buffer (wash buffer supplemented with 300 mM imidazole). Eluted protein was supplemented with 5 mM EGTA to chelate any remaining calcium and dialysed (10,000 MWCO) against 50 mM Tris pH 8.0, 300 mM NaCl, 0.5 mM TCEP buffer overnight at 4 °C.

For each loading procedure, an aliquot of EfPiwi (1/50 of the protein yield from 750 ml of Sf9 culture) was incubated with a synthetic piRNA guide (15 µM f.c.) for 15 min at room temperature and then dialysed into 50 mM Tris pH 8.0, 300 mM NaCl, 0.5 mM TCEP, 0.02% CHAPS buffer overnight at 4 °C (12,000 MWCO). To prepare for capturing guide-loaded EfPiwi, 2.5 nmol of biotinylated capture oligonucleotide was incubated with 40 µl high capacity neutravidin resin (ThermoFisher, 29204) in 1 ml wash A buffer (30 mM Tris pH 8.0, 0.1 M potassium acetate, 2 mM magnesium acetate, 0.02% CHAPS and 0.5 mM TCEP) for 30 min at 4 °C, followed by two washes with 2 ml wash A buffer. EfPiwi–guide complex was captured by incubating with the capture oligonucleotide-conjugated neutravidin resin at room temperature for 1.5 h with rotation. The resin was then washed three times with 2 ml wash A buffer, four times with 2 ml wash B buffer (30 mM Tris pH 8.0, 2 M potassium acetate, 2 mM magnesium acetate, 0.02% CHAPS and 0.5 mM TCEP), and three times with 2 ml wash C buffer (30 mM Tris pH 8.0, 1 M potassium acetate, 2 mM magnesium acetate, 0.02% CHAPS and 0.5 mM TCEP) at 4 °C. The resin was then resuspended in 250 µl wash C buffer containing biotinylated competitor oligonucleotide (50 µM f.c.) and incubated with rotation at room temperature for 3 h. To remove excess competitor oligonucleotide, the supernatant was incubated for 30 min at 4 °C with 60 µl fresh neutravidin resin (prewashed twice in wash C buffer), and the supernatant was dialysed overnight at 4 °C into extract buffer (30 mM HEPES-KOH, pH 7.5, 3.5 mM magnesium acetate, 2 mM DTT, 15% (v/v) glycerol and 0.01% (v/v) Triton X-100). The dialysed EfPiwi–guide RNA complex was aliquoted, flash frozen in liquid nitrogen and stored at −80 °C.

Recombinant mouse GTSF1 purification

pCold-GST GTSF-expression vectors were transformed into Rosetta-Gami 2 competent cells (Sigma, 71351). Cells were grown to an OD_600_ of 0.6–0.8 in the presence of 1 μM ZnSO_4_ at 37 °C, then chilled on ice for 30 min to initiate cold shock. Protein expression was induced with 0.5 mM IPTG for 18 h at 15 °C. Cells were collected by centrifugation, washed twice with PBS and cell pellets were flash frozen and stored at −80 °C. Cell pellets were resuspended in lysis/GST column buffer containing 20 mM Tris-HCl pH 7.5, 500 mM NaCl, 1 mM DTT, 5% (v/v) glycerol and 1× protease inhibitor cocktail (1 mM 4-(2-aminoethyl)benzenesulfonyl fluoride hydrochloride (Sigma, A8456), 0.3 μM aprotinin, 40 μM betanin hydrochloride, 10 μM E-64 (Sigma, E3132) and 10 μM leupeptin hemisulfate). Cells were lysed by a single pass at 18,000 p.s.i. through a high-pressure microfluidizer (Microfluidics, M110P), and the resulting lysate clarified at 30,000g for 1 h at 4 °C. Clarified lysate was filtered through a 0.22 µm Millex Durapore low-protein-binding syringe filter (EMD Millipore) and applied to glutathione Sepharose 4b resin (Cytiva, 17075604) equilibrated with GST column buffer. After draining the flow through, the resin was washed with 50 column-volumes of GST column buffer. To elute the bound protein and cleave the GST tag in a single step, 50 U HRV3C protease (Millipore, 71493) in 2.5 ml 20 mM Tris-HCl, pH 7.5, 50 mM NaCl, 1 mM DTT and 5% (v/v) glycerol was added to the column, and the column sealed and incubated for 3 h at 4 °C. Next, the column was drained to collect the cleaved protein. The eluate was diluted to 50 mM NaCl and further purified using a HiTrap Q (Cytiva, 29051325) anion exchange column equilibrated with 20 mM Tris-HCl, pH 7.5, 50 mM NaCl, 1 mM DTT and 5% (v/v) glycerol. The bound protein was eluted using a 100–500 mM NaCl gradient in the same buffer. Peak fractions were analysed for purity by SDS–PAGE and the purest were pooled and dialysed into storage buffer containing 30 mM HEPES-KOH, pH 7.5, 100 mM potassium acetate, 3.5 mM magnesium acetate, 1 mM DTT and 20% (v/v) glycerol. Aliquots of the pooled fractions were flash frozen in liquid nitrogen and stored at −80 °C.

Recombinant EmGtsf1 purification

The high-quality draft genome of E. muelleri was used to design the expression construct of Ephydatia sp. Gtsf1 orthologue. The EmGtsf1 expression vector was transformed into BL21(DE3) cells (NEB, C2527H). Transformed cells were grown in LB medium supplemented with 1  μM ZnSO_4_ at 37 °C until an OD_600_ of 0.6–0.8. The incubation temperature was lowered to 16 °C and protein expression was induced by the addition of 1 mM IPTG for 16 h. Cells were collected by centrifugation and cell pellets flash frozen in liquid nitrogen and stored at −80 °C. Thawed cell pellets were resuspended in lysis buffer (50 mM Tris, pH 8, 300 mM NaCl and 0.5 mM TCEP) and passed through a high-pressure (18,000 p.s.i.) microfluidizer (Microfluidics, M110P) to induce cell lysis. The lysate was clarified by centrifugation at 30,000g for 20 min at 4 °C. Clarified lysate was applied to Ni-NTA resin (Qiagen) and incubated for 1 h. The resin was washed with nickel wash buffer (300 mM NaCl, 20 mM imidazole, 0.5 mM TCEP and 50 mM Tris, pH 8.0). Protein was eluted in four column volumes of nickel elution buffer (300 mM NaCl, 300 mM imidazole, 0.5 mM TCEP and 50 mM Tris, pH 8.0). TEV protease was added to the eluted protein to remove the amino-terminal His_6_ and MBP tags. The resulting mixture was dialysed against HiTrap dialysis buffer (300 mM NaCl, 20 mM imidazole, 0.5 mM TCEP and 50 mM Tris, pH 8.0) at 4 °C overnight. The dialysed protein was then passed through a 5 ml HiTrap chelating column (Cytiva) and the unbound material collected. Unbound material was concentrated and further purified by size-exclusion chromatography using a Superdex 75 Increase 10/300 column (Cytiva) equilibrated in 50 mM Tris, pH 8.0, 300 mM NaCl and 0.5 mM TCEP. Peak fractions were analysed for purity by SDS–PAGE, and the purest were pooled, concentrated to 100 µM, aliquoted and stored at −80 °C.

Determination of the active fraction of piRISC

In vitro cleavage assays were used to determine the fraction of active piRISC. Target RNA substrates for cleavage assays were prepared as previously described^34^. Fully complementary piRNA target site-containing templates were PCR amplified from pGL2 (primers listed in Supplementary Table 4), in vitro transcribed with T7 RNA polymerase, purified using a 7% denaturing polyacrylamide gel and capped using α-[^32^P]GTP (Perkin Elmer) and a Vaccinia Capping System (NEB, M2080S). Unincorporated α-[^32^P]GTP was removed using a G-25 spin column (Cytiva, 27532501), target RNA was purified using a 7% denaturing polyacrylamide gel, eluted overnight with rotation in 0.4 M NaCl at 4 °C and collected by ethanol precipitation. Radiolabelled target (10 nM f.c.) was added to a mix of purified piRISC and GTSF1 (500 nM f.c.) to assemble a 30 μl cleavage reaction. At 0, 5, 15, 30 and 60 min, a 5 μl sample was quenched in 280 μl 50 mM Tris-HCl, pH 7.5, 100 mM NaCl, 25 mM EDTA and 1% (w/v) SDS, then proteinase K (1 mg ml^–1^ f.c.) was added and the mix incubated at 45 °C for 15 min, followed by extraction with phenol–chloroform–isoamyl alcohol (25:24:1, pH 6.7) and ethanol precipitation. RNA was resuspended in 10 μl 95% (v/v) formamide, 5 mM EDTA, 0.025% (w/v) bromophenol blue and 0.025% (w/v) xylene cyanol, heated at 95 °C for 2 min, and resolved on a 7% denaturing polyacrylamide gel. Gels were dried, exposed to a storage phosphor screen and imaged on a Typhoon FLA 7000 (GE). The raw image file was used to quantify the substrate and product bands, corrected for background. Data were fit to the burst-and-steady-state equation to determine the concentration of active piRISC (see equation and fitting procedure in the section ‘Analysis of CNS data’).

RNA bind-’n-Seq for Kd measurements

RNA bind-’n-Seq (RBNS) was performed as previously described^32^ with modifications. A library of RNA oligonucleotides containing a central region of 20 random-sequence positions (Extended Data Fig. 1a) was obtained from IDT, 5′[^32^P]-radiolabelled with α-[^32^P]GTP (Perkin Elmer) and T4 PNK (NEB, M0201) and purified using a 15% denaturing polyacrylamide gel, extracted with phenol–chloroform–isoamyl alcohol (25:24:1, pH 6.7) and collected by ethanol precipitation. To sequence the input library, RNA was denatured at 90 °C for 1 min, annealed to a RT primer (Supplementary Table 4) and reverse transcribed with SuperScript III. RNA was degraded by alkaline hydrolysis in 0.4 M NaOH for 1 h at 55 °C, and cDNA was recovered by ethanol precipitation. The sample was then amplified in 25 µl using AccuPrime Pfx DNA polymerase (ThermoFisher, 12344024; 95 °C for 2 min, 15 cycles of 95 °C for 15 s, 65 °C for 30 s, 68 °C for 15 s; primers listed in Supplementary Table 4). PCR products were purified with a 2% agarose gel and sequenced on a NextSeq 550 (Illumina) to obtain 79-nucleotide, single-end reads.

For RBNS, DNA-blocking oligonucleotides (Supplementary Table 4) were annealed to the RNA library in 30 mM HEPES-KOH, pH 7.5, 120 mM potassium acetate and 3.5 mM magnesium acetate using a 1:1.2 molar ratio of RNA pool to DNA blockers by first incubating at 95 °C for 1 min, then at 65 °C for 10 min and finally cooled to room temperature. For each trial, the final piRISC concentrations in the six RBNS reactions were 0.003 nM, 0.01 nM, 0.032 nM, 0.1 nM, 0.316 nM and 1 nM active piRISC. Each trial also included a control in which protein storage buffer replaced piRISC. Binding for each piRISC concentration was performed in 20 μl 25 mM HEPES-KOH, pH 7.9, 110 mM potassium acetate, 3.5 mM magnesium acetate, 0.01% (w/v) Triton X-100, 2 mM DTT, 10% (w/v) glycerol and 100 nM (f.c.) RNA library. To reduce non-specific binding, each reaction also included 2.5 mg ml^−1^ BSA and 0.5 mg ml^−1^ yeast tRNA. Reactions were incubated for 2 h at 33 °C (ref. ^51^) and then filtered through a Whatman Protran nitrocellulose membrane (Sigma, WHA10402506) on top of an Amersham Hybond-XL (Cytiva, RPN2222S) nylon membrane in a Bio-Dot apparatus (Bio-Rad, 1706545). To reduce the retention of free single-stranded RNA, nitrocellulose and nylon membranes were pre-conditioned as previously described^52,53^. In brief, the nitrocellulose filter was pre-soaked in 0.4 M KOH for 10 min, and the nylon filter was incubated in 0.1 M EDTA, pH 8.2 for 10 min, washed three times in 1 M NaCl for 10 min each followed by an 15 s rinse in 0.5 M NaOH. Nitrocellulose and nylon filters were then rinsed in water until the pH returned to neutral and then equilibrated in wash buffer (20 mM HEPES-KOH, pH 7.9, 100 mM potassium acetate, 3.5 mM magnesium acetate and 1 mM DTT) for at least 1 h at 37 °C. After applying seven samples (no-piRISC and six piRISC concentrations; Extended Data Fig. 1a) onto the nitrocellulose and nylon membrane under vacuum, the two membranes were washed with 100 ml wash buffer for 3 s. Membranes were air-dried and signals detected by phosphorimaging to monitor binding. The nitrocellulose membrane areas containing piRISC-bound RNA were excised and incubated with 1 mg ml^–1^ proteinase K in 100 mM Tris-HCl, pH 7.5, 10 mM EDTA, 150 mM sodium chloride and 1% (w/v) SDS for 1 h at 45 °C shaking at 300 r.p.m. After phenol–chloroform extraction and ethanol precipitation, RNA was reverse transcribed, amplified and sequenced as described above for the input RNA pool.

Determining equilibrium dissociation constants by double-filter binding assays

Binding assays were performed as previously described^34^ in 5 μl in 30 mM HEPES-KOH, pH 7.5, 120 mM potassium acetate, 3.5 mM magnesium acetate, 2 mM DTT and 0.01% (w/v) Triton X-100. The 5′[^32^P]-radiolabelled RNA targets (0.1 nM; listed in Supplementary Table 4) were incubated with 0.001–0.8 nM piRISC. The assay also included a control reaction using piRISC storage buffer. Binding reactions were incubated at 33 °C for 2 h. RNA binding was measured by capturing protein–RNA complexes on Protran nitrocellulose (GE, GE10600002) and unbound RNA on a Hybond-XL (Cytiva, 45-001-151) in a Bio-Dot apparatus (Bio-Rad). After applying the sample under vacuum, membranes were washed with 10 μl equilibration buffer (30 mM HEPES-KOH, pH 7.5, 120 mM potassium acetate, 3.5 mM magnesium acetate and 2 mM DTT). Membranes were air-dried and signals detected by phosphorimaging. Because Kd < [RNA target], all binding data were fit to the following equation using IgorPro 6.11 (WaveMetrics):

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f=\frac{\left(\left[{E}_{{\rm{T}}}\right]+\left[{S}_{{\rm{T}}}\right]+{K}_{{\rm{d}}}\right)-\sqrt{{\left(\left[{E}_{{\rm{T}}}\right]+\left[{S}_{{\rm{T}}}\right]+{K}_{{\rm{d}}}\right)}^{2}-4\left[{E}_{{\rm{T}}}\right]\left[{S}_{{\rm{T}}}\right]}}{2\left[{S}_{{\rm{T}}}\right]},$$\end{document}

where f is the fraction target bound, [ET] is the total piRISC concentration, [ST] is the total RNA target concentration and Kd is the apparent equilibrium dissociation constant.

Cleave-’n-Seq to determine target cleavage rates

To the ssDNA oligonucleotide pool of Cleave-’n-Seq (CNS) targets (Extended Data Fig. 3b) obtained from TWIST Bioscience, a T7 promoter was added by PCR (primers listed in Supplementary Table 4). The PCR products were in vitro transcribed with T7 RNA polymerase, then treated with TURBO DNase (ThermoFisher, AM2238), and the CNS target RNA library was purified using a 7% denaturing polyacrylamide gel.

DNA-blocking oligonucleotides (Supplementary Table 4) in 1.2-fold excess were annealed to 100 nM CNS target RNA library in 10 mM Tris-HCl (pH 7.4) and 20 mM NaCl by heating the mixture to 95 °C, cooling it at −0.1 °C s^–1^ to 22 °C and incubating at 22 °C for 5 min. The 100 nM target RNA library was diluted with water to 1 nM, aliquoted and stored at −80 °C. Cleavage assays were performed in 20 mM HEPES-KOH (pH 7.5), 80 mM potassium acetate, 3.5 mM magnesium acetate, 4 mM DTT, 10% (v/v) glycerol and 0.01% (v/v) Triton X-100. Each reaction contained 0.1 nM RNA library, 500 nM GTSF1 and 1 nM active piRISC (single turnover conditions). Reactions were conducted by prewarming components to 33 °C, first mixing piRISC with GTSF1 and then adding target libraries, immediately before incubating at 33 °C (ref. ^51^) for 0, 1, 2, 4 and 8 min or 0, 20, 60, 120, 240, 480 and 960 min. At each time point, a 5 μl sample was quenched in 280 μl 50 mM Tris-HCl, pH 7.5, 100 mM NaCl, 25 mM EDTA and 1% (w/v) SDS, then proteinase K (1 mg ml^–1^ f.c.) added, and the mix incubated at 45 °C for 15 min followed by extraction with phenol–chloroform–isoamyl alcohol (25:24:1, pH 6.7). The RNA was collected by ethanol precipitation, resuspended in 10 μl water, denatured at 90 °C for 1 min, annealed with 1 μl 50 μM RT primer (Supplementary Table 4) at 65 °C for 5 min, and reverse transcribed with SuperScript III (ThermoFisher, 18080093). cDNA was recovered by ethanol precipitation, and the sample was then amplified in 25 µl using AccuPrime Pfx DNA polymerase (ThermoFisher, 12344024; 95 °C for 2 min, 15 cycles of 95 °C for 15 s, 65 °C for 30 s and 68 °C for 15 s; primers listed in Supplementary Table 4). PCR products were purified using a 2% agarose gel and sequenced on a NextSeq 550 (Illumina) to obtain 60-nucleotide, single-end reads. All time points in each trial (that is, both the 0–8 min and the 0–960 min subsets) were sequenced in the same NextSeq 550 run. Data from three trials of each of the 0–8 min and 0–960 min subsets were combined to estimate pre-steady-state cleavage rates.

Cloning and sequencing 3′ cleavage products from CNS reactions

For Extended Data Fig. 5c, a modified DNA library of CNS targets with 8 nucleotide barcodes (each unique to a target variant) was obtained from TWIST Bioscience (Supplementary Table 5). The procedure was identical to the CNS protocol described in the previous section except for the addition of the 5′ adapter ligation step after annealing the RT primer and before reverse transcription: 3′ cleavage products were ligated to a mixed pool of equimolar amount of two 5′ RNA adapters (to increase nucleotide diversity at the 5′ end of the sequencing read; Supplementary Table 4) in 20 µl 50 mM Tris-HCl (pH 7.8), 10 mM MgCl_2_, 10 mM DTT, 1 mM ATP with 60 U high concentration T4 RNA ligase (NEB, M0437M) at 16 °C overnight. Ligation was followed by ethanol precipitation. Cleavage reactions were for performed for 2 h at 33 °C. To account for 5′-to-3′ exonucleolytic digestion or addition of non-templated nucleotides to RNA 5′ ends, a set of five synthetic 5′ monophosphorylated oligonucleotides (Supplementary Table 4) was added to the CNS library before starting the cleavage reaction.

FACS isolation and immunostaining of mouse germ cells

Testes of 2–6-month-old mice were isolated, decapsulated and incubated for 15 min at 33 °C in 1× Gey′s balanced salt solution (GBSS, Sigma, G9779) containing 0.4 mg ml^–1^ collagenase type 4 (Worthington, LS004188) rotating at 150 r.p.m. Seminiferous tubules were then washed twice with 1× GBSS and incubated for 15 min at 33 °C in 1× GBSS with 0.5 mg ml^–1^ trypsin and 1 µg ml^–1^ DNase I, rotating at 150 r.p.m. Next, tubules were homogenized by pipetting through a glass Pasteur pipette for 3 min at 4 °C. FBS (7.5% f.c., v/v) was added to inactivate trypsin, and the cell suspension was then strained through a pre-wetted 70 µm cell strainer (ThermoFisher, 22363548). Cells were collected by centrifugation at 300g for 10 min. The supernatant was removed, cells resuspended in 1× GBSS containing 5% (v/v) FBS, 1 µg ml^–1^ DNase I and 5 μg ml^–1^ Hoechst 33342 (ThermoFisher, 62249) and rotated at 150 r.p.m. for 45 min at 33 °C. Propidium iodide (0.2 μg ml^–1^, f.c.; ThermoFisher, P3566) was added, and cells strained through a pre-wetted 40 µm cell strainer (ThermoFisher, 22363547). Spermatogonia, primary spermatocytes, secondary spermatocytes and round spermatids were purified^48,54^ (Supplementary Fig. 5) using a FACSAria II Cell Sorter (BD Biosciences; University of Massachusetts Medical School FACS Core). The 355 nm laser was used to excite Hoechst 33342; the 488 nm laser was used to record forward and side scatter and to excite propidium iodide. Propidium iodide emission was detected using a 610/20 bandpass filter. Hoechst 33342 emission was recorded using 450/50 and 670/50 band pass filters.

Germ cell stages in the unsorted population and the purity of sorted fractions were assessed by immunostaining aliquots of cells. Cells were incubated for 20 min in 25 mM sucrose and then fixed on a slide with 1% (w/v) paraformaldehyde containing 0.15% (v/v) Triton X-100 for 2 h at room temperature in a humidifying chamber. Slides were washed sequentially for 10 min in the following solutions: (1) PBS containing 0.4% (v/v) Photo-Flo 200 (Kodak, 1464510); (2) PBS containing 0.1% (v/v) Triton X-100; and (3) PBS containing 0.3% (w/v) BSA, 1% (v/v) donkey serum (Sigma, D9663), and 0.05% (v/v) Triton X-100. After washing, slides were incubated with primary antibodies in PBS containing 3% (w/v) BSA, 10% (v/v) donkey serum and 0.5% (v/v) Triton X-100 overnight at room temperature in a humidified chamber. Rabbit polyclonal anti-SYCP3 (Abcam, ab15093, RRID:AB_301639; 1:1,000 dilution) and mouse monoclonal anti-γH2AX (Millipore, 05-636, RRID:AB_309864; 1:1,000 dilution) were used as primary antibodies. Slides were washed again as described above and then incubated with secondary donkey anti-mouse IgG (H+L) Alexa Fluor 594 (ThermoFisher, A-21203, RRID:AB_2535789; 1:2,000 dilution) or donkey anti-rabbit IgG (H+L) Alexa Fluor 488 (ThermoFisher, A-21206, RRID:AB_2535792; 1:2,000 dilution) for 1 h at room temperature in a humidified chamber. After incubation, slides were washed three times (10 min each) in PBS containing 0.4% (v/v) Photo-Flo 200 and once for 10 min in 0.4% (v/v) Photo-Flo 200. Finally, slides were dried and mounted in ProLong Gold antifade mountant with DAPI (ThermoFisher, P36931). To assess the purity of sorted fractions, 50–100 cells were staged by DNA, γH2AX and SYCP3 staining^54^. All samples used in this study met the following criteria: spermatogonia, 95–100% pure with ≤5% pre-leptotene spermatocytes; primary spermatocytes, 10–15% leptotene/zygotene spermatocytes, 45–50% pachytene spermatocytes, 35–40% diplotene spermatocytes; secondary spermatocytes, 100%; round spermatids, 95–100%, ≤5% elongated spermatids.

Small RNA sequencing library preparation

Total RNA from sorted mouse germ cells was extracted using a mirVana miRNA isolation kit (ThermoFisher, AM1560). Small RNA libraries were constructed as previously described^48^ with modifications. Before library preparation, an equimolar mix of nine synthetic spike-in RNA oligonucleotides (Supplementary Table 4) was added to each RNA sample to enable absolute quantification of small RNAs (Supplementary Table 6); the median cell volume from ref. ^21^ was used to calculate the intracellular concentration. To reduce ligation bias and to eliminate PCR duplicates, the 3′ and 5′ adaptors both contained nine random nucleotides at their 5′ and 3′ ends, respectively^55^ (Supplementary Table 4) and 3′ adaptor ligation reactions contained 25% (w/v) PEG-8000 (f.c.): 500–1,000 ng total RNA was first ligated to 25 pmol of 3′ DNA adapter (Supplementary Table 4) with adenylated 5′ and dideoxycytosine-blocked 3′ ends in 30 µl of 50 mM Tris-HCl (pH 7.5), 10 mM MgCl_2_, 10 mM DTT and 25% (w/v) PEG-8000 (NEB) with 600 U of homemade T4 Rnl2tr K227Q at 16 °C overnight. After ethanol precipitation, the 50–90 nucleotide (14–54 nucleotide small RNA + 36 nucleotide 3′ adapter containing unique molecular identifiers) 3′ ligated product was purified from a 15% denaturing urea–polyacrylamide gel (National Diagnostics). After overnight elution in 0.4 M NaCl followed by ethanol precipitation, the 3′ ligated product was denatured in 14 µl water at 90 °C for 60 s, 1 µl of 50 µM RT primer (Supplementary Table 4) was added and annealed at 65 °C for 5 min to suppress the formation of 5′-adapter–3′-adapter dimers during the next step. The resulting mix was then ligated to a mixed pool of equimolar amount of two 5′ RNA adapters (to increase the nucleotide diversity at the 5′ end of the sequencing read; Supplementary Table 4) in 20 µl of 50 mM Tris-HCl (pH 7.8), 10 mM MgCl_2_, 10 mM DTT and 1 mM ATP with 20 U of T4 RNA ligase (ThermoFisher, EL0021) at 25 °C for 2 h. The ligated product was precipitated with ethanol, cDNA synthesis was performed in 20 µl at 42 °C for 1 h using AMV reverse transcriptase (NEB, M0277), and 5 µl of the RT reaction was amplified in 25 µl using AccuPrime Pfx DNA polymerase (ThermoFisher, 12344024; 95 °C for 2 min, 15 cycles of 95 °C for 15 s, 65 °C for 30 s and 68 °C for 15 s; primers listed in Supplementary Table 4). Finally, the PCR product was purified in a 2% agarose gel. Small RNA sequencing (RNA-seq) libraries samples were sequenced using a NextSeq 550 (Illumina) to obtain 79 nucleotide, single-end reads.

RNA-seq library preparation

Total RNA from sorted germ cells was extracted using a mirVana miRNA isolation kit (ThermoFisher, AM1560) and used for library preparation as previously described^56^ with modifications, including the addition of the ERCC spike-in mix to enable absolute quantification of RNAs and the use of unique molecular identifiers in adapters (Supplementary Table 4) to eliminate PCR duplicates^55^. Before library preparation, 1 µl of 1:100 diluted ERCC spike-in mix 1 (ThermoFisher, 4456740) was added to 1 µg total RNA. To remove rRNA, 1 µg total RNA was hybridized in 10 µl to a pool of 186 rRNA antisense oligos (0.05 µM f.c. each) in 10 mM Tris-HCl (pH 7.4), 20 mM NaCl by heating the mixture to 95 °C, cooling at −0.1 °C s^–1^ to 22 °C, and incubating at 22 °C for 5 min. RNase H (10 U; Lucigen, H39500) was added and the mixture incubated at 45 °C for 30 min in 20 µl containing 50 mM Tris-HCl (pH 7.4), 100 mM NaCl and 20 mM MgCl_2_. The reaction volume was adjusted to 50 µl with 1× Turbo DNase buffer (ThermoFisher, AM2238) and then incubated with 4 U Turbo DNase (ThermoFisher, AM2238) for 20 min at 37 °C. Next, RNA was purified using RNA Clean & Concentrator-5 (Zymo Research, R1016) to retain ≥200 nucleotide RNAs, followed by the stranded, dUTP-based RNA-seq protocol as previously described^56^. RNA-seq libraries were sequenced using a NextSeq 550 (Illumina) to obtain 79+79 nucleotide, paired-end reads.

Sequencing of 5′-monophosphorylated long RNAs

Total RNA from sorted mouse germ cells was extracted using a mirVana miRNA isolation kit (ThermoFisher, AM1560) and used to prepare a library of 5′-monophosphorylated long RNAs as previously described^21,36^ with modifications. rRNA was depleted as described above for RNA-seq libraries. RNA was ligated to a mixed pool of equimolar amount of two 5′ RNA adapters (to increase the nucleotide diversity at the 5′ end of the sequencing read; Supplementary Table 4) in 20 µl of 50 mM Tris-HCl (pH 7.8), 10 mM MgCl_2_, 10 mM DTT and 1 mM ATP with 60 U of high concentration T4 RNA ligase (NEB, M0437M) at 16 °C overnight. The ligated product was isolated using RNA Clean & Concentrator-5 (Zymo Research, R1016) to retain ≥200 nucleotide RNAs and reverse transcribed in 25 µl with 50 pmol RT primer (Supplementary Table 4) using SuperScript III (ThermoFisher, 18080093). After purification with 50 µl Ampure XP beads (Beckman Coulter, A63880), cDNA was PCR amplified using NEBNext High-Fidelity (NEB, M0541; 98 °C for 30 s; 4 cycles of 98 °C for 10 s, 59 °C for 30 s, 72 °C for 12 s; 6 cycles of 98 °C for 10 s, 68 °C for 10 s, 72 °C for 12 s; and 72 °C for 3 min; primers listed in Supplementary Table 4). PCR products between 200 and 400 bp were isolated from a 1% agarose gel, purified using a QIAquick Gel Extraction kit (Qiagen, 28706), and amplified again with NEBNext High-Fidelity (NEB, M0541; 98 °C for 30 s; 3 cycles of 98 °C for 10 s, 68 °C for 30 s, 72 °C for 14 s; 6 cycles of 98 °C for 10 s, 72 °C for 14 s; and 72 °C for 3 min; primers listed in Supplementary Table 4). The PCR product was purified from a 1% agarose gel and sequenced using a NextSeq 550 or NovaSeq 6000 (Illumina) to obtain 79+79 nucleotide or 150+150 nucleotide, paired-end reads.

Analysis of RBNS data

To analyse RBNS^31^ data, the sequence of the 3′ adapter (5′-TGGAATTCTCGGGTGCCAAGG-3′) was removed using fastx toolkit (v.0.0.14), then each sequencing read in the RNA input library and piRISC-bound libraries was interrogated for the presence of all binding sites of interest. The entire single-stranded 20 nucleotide random-sequence region flanked by four nucleotides of constant primer-binding sequence on either side (GATCNNNNNNNNNNNNNNNNNNNNTGGA) was searched for the presence of piRISC-binding sites. The sequencing depth of the input library (about 50 × 10^6^ reads) allowed measurement of input frequencies for ≤12 nucleotide motifs. To interrogate non-overlapping target sets, each ≤10 nucleotide contiguous binding site was required to be flanked by nucleotides that not complementary to the guide: for example, a g4–g12 contiguous target site did not pair to guide positions g3 and g13. Each 11-nucleotide-long contiguous complementary site was required to be flanked by a non-matching nucleotide only at its 5′ end: for example, a g4–g14 contiguous target site did not pair to guide position g3. To eliminate interference from potential piRISC cleavage activity, GTSF1 was omitted from binding reactions; we also relied on the fact that, in our analyses, we do not interrogate sites that are long enough (≥15 nucleotides) to be cleaved by piRISC.

A read was assigned to a site category if it contained one single binding motif. Reads containing multiple instances of binding sites (from the same or a different site category) and reads containing partially overlapping sites were not included in the analysis. Reads that did not have any of the binding motifs of interest were classified as reads with no site. Fitting of the binding model from a previously described method^32^ to estimate Kd values for binding sites was performed using a Python-based implementation (MLE_KD.py from https://figshare.com/articles/software/MicroRNA-binding_thermodynamics_and_kinetics_by_RNA_Bind-n-Seq/19180952) on each of the 49 different combinations of 7 initial guesses of piRISC concentration (0.1, 0.2, 0.5, 1, 2, 5 or 10 nM) and 7 initial guesses of Kd for RNA with no enriched site (0.1, 0.2, 0.5, 1, 2, 5 or 10 nM). For each trial, the median of the 49 estimates was reported. Two trials of mouse AGO2 RBNS data for let-7a (piRNA-1 in Fig. 1) are from a previous study^32^; the third trial was conducted separately for this study. All other mouse AGO2 RBNS data are from a previous study^32^. All human AGO2 RBNS data are from a previous study^57^. AGO2 RBNS data were downloaded from the National Center for Biotechnology Information and analysed using the same binding model as previously described^32^. Predicted binding energy, ΔG^0^, was estimated using the RNAplex nearest neighbour algorithm^58^.

Analysis of CNS data

After the sequence of the 3′ adapter (5′-TGGAATTCTCGGGTGCCAAGG-3′) was removed using fastx toolkit (v.0.0.14), CNS library target sites (Supplementary Table 1) were identified without allowing mismatches or insertions or deletions. The 8 nucleotide barcodes were used when 3′ cleavage products were cloned and sequenced (Supplementary Table 5). Sequencing data (representing the abundance of uncleaved targets) were first normalized to the sequencing depth (parts per million (ppm)). To adjust for the decrease in total abundance of the library over the course of cleavage reaction, each ppm value was divided by the sum of ppm values of targets that contained ≤7 nucleotide complementarity to the piRISC piRNA guide. Next, the relative abundance of cleaved product at non-zero time points was inferred as follows: [Prelative] = (ppm_0 min_ − ppm_X min_)/ppm_0 min_. [Prelative] ranged from 0 (no cleaved product) to 1 (all substrate cleaved). The combined [Prelative] data from three independent trials of each 0–8 min and 0–960 min subsets (that is, three trials of each 1, 2, 4, 8, 20, 60, 120, 240, 480 and 960 min) were used to fit the burst-and-steady-state scheme \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E+S\mathop{\mathop{{\rm{\rightleftharpoons }}}\limits^{{k}_{1}}}\limits_{{k}_{-1}}ES\mathop{\to }\limits^{{k}_{2}}EP\mathop{\to }\limits^{{k}_{3}}E+P$$\end{document} , using equation:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}[{P}_{{\rm{relative}}}]=f({\rm{t}})=[{E}_{{\rm{relative}}}]({[{k}_{2}/({k}_{2}+{k}_{3})]}^{2}\times (1-{e}^{-[k2+k3]{\rm{t}}})\\ \,\,\,+\,[{k}_{2}{k}_{3}/({k}_{2}+{k}_{3}){\rm{t}}])([{E}_{{\rm{relative}}}]).\end{array}$$\end{document}

The fit was performed using the Trust Region Reflective algorithm implemented in the optimize.curve_fit function from Python module scipy (v.1.8.1)^59^ for the maximum number of 10,000 function evaluations before the termination. The following physically meaningful constraints on the parameters were used: 0.5 ≤ [Erelative] ≤ 1; 0 ≤ k2 ≤ 100 min^−1^; and for the single turnover experiment setup, 0 ≤ k3 ≤ 0.0001 min^−1^. For each fitting procedure, the mean and the standard deviation of the estimate for each parameter are reported in a Supplementary Table 1. The resulting (k2 + k3) was reported as the pre-steady-state cleavage rate (k).

Mouse AGO2 CNS data for let-7a and miR-21 RISCs are from a previous study^3^, and mouse AGO2 CNS data for L1MC RISC was generated for this study.

Analysis of small RNA datasets

The 3′ adapter (5′-TGGAATTCTCGGGTGCCAAGG-3′) was removed using fastx toolkit (v.0.0.14), and PCR duplicates were eliminated as previously described^55^. rRNA matching reads were removed using bowtie (parameter -v 1; v.1.0.0)^60^ against the M. musculus set in the SILVA rRNA database^61^. Deduplicated and filtered data were analysed using Tailor^62^ to account for non-templated tailing of small RNAs. Sequences of synthetic spike-in oligonucleotides (Supplementary Table 4) were identified, allowing no mismatches with bowtie (parameter -v 0; v.1.0.0)^60^, and the absolute abundance of small RNAs calculated (Supplementary Table 6). Because piRNA 3′ trimming by PNLDC1 results in piRNA 3′ end heterogeneity, sequencing reads were next grouped by their 5′, 25 nucleotide prefix. For further analyses, we kept only prefix groups that met three criteria. First, the 25 nucleotide prefix unambiguously mapped to a single genomic position (>80% of the 25 nucleotide piRNA prefixes met this criterion). Second, the prefix group total abundance was ≥1 ppm (that is, ≥10 piRNAs/mouse primary spermatocyte), ensuring that, assuming a Poisson or a negative binomial distribution for piRNA concentration in different cells, ≥99.99% of primary spermatocytes contained at least 1 molecule of the piRNA 25 nucleotide prefix. Third, the prefix group total abundance was ≥1 ppm in all 12 replicates of control C57BL/6 samples (Supplementary Table 6). piRNAs were considered undetectable in pi2^−/−^pi9^−/−^pi17^−/−^ mutants if their mean abundance in mutants (n = 9) was ≤0.1 ppm.

RNA-seq analysis

RNA-seq analysis was performed using piPipes for genomic alignment^63^. Before starting piPipes, sequences were reformatted to extract unique molecular identifiers^55^. The reformatted reads were then aligned to rRNA using bowtie2 (v.2.2.0)^64^. Unaligned reads were mapped to mouse genome mm10 using STAR (v.2.3.1)^65^ and PCR duplicates were removed^55^. Transcript abundance was calculated using StringTie (v.1.3.4)^66^. Differential expression analysis was performed using DESeq2 (v.1.18.1)^67^. In parallel, reformatted reads were aligned to an index of ERCC spike-in transcripts (ThermoFisher, 4456740) using bowtie (v.1.0.0)^60^, PCR duplicates were removed as previously described^55^ and the absolute quantity of transcripts calculated (Supplementary Table 7).

Analysis of 5′-monophosphorylated long RNA-sequencing data

Sequencing data for 5′-monophosphorylated long RNAs was aligned to the mouse genome using piPipes^63^. Before starting piPipes, the degenerate portions of the 5′ adapter sequences were removed (nucleotides 1–15 of read 1). Because each library was sequenced at least twice to increase the sequencing depth, to harmonize the length of paired-end reads from different runs, sequences were trimmed to 64 nucleotide (read 1) + 79 nucleotide (read 2) paired reads. The trimmed reads were then aligned to rRNA using bowtie2 (v.2.2.0)^64^. Unaligned reads were mapped to mouse genome mm10 using STAR (v.2.3.1)^65^, alignments with soft clipping of ends were removed using SAMtools (v.1.0.0)^68^ and reads with the same 5′ end were merged to represent a single 5′-monophosphorylated RNA species. For further analyses, only unambiguously mapping 5′-monophosphorylated RNA species for which abundance was ≥0.04 ppm were used. For 5′-monophosphorylated RNAs mapped in annotated transcripts, the nucleotide sequence of the corresponding transcript was used to find piRNAs potentially explaining the cleavage, and we used the genomic sequence for 5′-monophosphorylated RNAs mapped outside any annotated transcript. To ensure that piRNA–target combinations for all pairing configurations did not overlap, the piRNA nucleotide immediately after the paired region was required to be unpaired with the target: for example, for g2–g10, g11 was unpaired and thus did not overlap with g2–g11, g2–g12, an so on. Calculating of the fraction of cleaved sites was performed for a collapsed, non-redundant list of cleavage sites, that is, even if a cleavage site was explained by several piRNAs, it was counted only once. Cumulative abundance of all piRNAs explaining each site was used to assess the effect of piRNA concentration.

Logistic regression classifier implementation

For each of the 16 permutations of 4 C57BL/6 control and 4 pi2^−/−^pi9^−/−^pi17^−/−^ mutant primary spermatocyte datasets, we identified 3,150–3,750 5′-monophosphorylated RNAs (that is, potential 3′ cleavage products of piRNA-guided slicing) for which abundance was ≥0.1 ppm and that were explained by ≥19 paired nucleotides between g2 and g25 of pi2, pi9 and pi17 piRNAs (target insertions or deletions were not allowed). Note that although abundance and binding energy remained the best predictive features regardless of the minimum number of paired nucleotides used as a threshold, requiring <19 paired nucleotides produced too few piRNA–target data points to inform the importance of pairing to piRNA 5′ terminal nucleotides. A target site was considered cleaved, that is, P(cleaved) = 1, if the abundance of the 5′-monophosphorylated RNAs decreased by ≥8-fold in pi2^−/−^pi9^−/−^pi17^−/−^ mutants compared with C57BL/6 controls. All other sites were assigned as uncleaved, that is, P(cleaved) = 0.

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P({\rm{c}}{\rm{l}}{\rm{e}}{\rm{a}}{\rm{v}}{\rm{e}}{\rm{d}})=\frac{1}{1+{{\rm{e}}}^{-({\beta }_{0}+{\beta }_{1}{x}_{1}+{\beta }_{2}{x}_{2}+\ldots +{\beta }_{35}{x}_{35})}}$$\end{document}

The logistic function representing the probability of target site cleavage, P(cleaved), contained 35 independent variables as follows: x1–x24: absence (0) or presence (1) of pairing to g2–g25; x25: total number of paired nucleotides between g2–g25, rescaled to [0,1]; x26: piRNA abundance, that is, the total abundance of all piRNAs with the same 25 nucleotide, 5′ prefix (see the section ‘Analysis of small RNA datasets’), rescaled to [0,1]; x27: negative of the predicted energy of piRNA–target pairing ΔG^0^ estimated with RNAplex^58^, rescaled to [0,1] (use of the negative of ΔG^0^ creates a positive relationship between strength of binding and probability of cleavage). Moreover, x28: equals 1 if t1A, 0 if t1B; x29: equals 1 if t1U, 0 if t1V; x30: equals 1 if t1C, 0 if t1D; x31: equals 1 if t1G, 0 if t1H; x32: equals 1 target site is in the 5′ UTR, 0 if outside the 5′ UTR; x33: equals 1 if the target site is in the ORF, 0 if outside the ORF; x34: equals 1 if the target site is in the 3′ UTR, 0 if outside the 3′ UTR; x35: equals 1 if the target site is in lncRNA, 0 if in mRNA.

The logistic function was fit using the Limited-memory Broyden–Fletcher–Goldfarb–Shanno algorithm (L-BFGS) implemented in LogisticRegression from the Python module scikit-learn^69^ using L2-regularization (λ = 1) with default parameters on acceptance of convergence and the maximum number of iterations set at 1,000. To balance cleaved and uncleaved classes, weights inversely proportional to class frequencies were used. RepeatedStratifiedKFold and cross_validate from scikit-learn were used to perform 5× repeated 5-fold cross validation, resulting in 5 × 5 = 25 logistic function fits for each of the 16 permutations of 4 control 4 mutant datasets, generating the total of 25 × 16 = 400 logistic regression models. To assess model performance, area under the precision-recall curve (AUC) for each of the 400 logistic functions was calculated either with the corresponding pi2^−/−^pi9^−/−^pi17^−/−^ dataset (400 AUC values total) or with each of the 16 permutations of 4 C57BL/6 and 4 pi7^−/−^ mutant datasets (6,400 AUC values total).

Simulation of transposon sequence mutagenesis

The consensus sequences of active mouse LINE transposons^70,71^ were mutagenized by adding 1,000 random single-nucleotide substitutions at ratios that reflect the established mouse germline mutation rates^72^. Only synonymous substitutions were allowed in LINE ORFs, and 100 independent simulations were performed for each consensus sequence. piRNAs from fetal mouse testis (embryonic day 16.5) were sequenced and used for the analyses, and 21 nucleotide siRNAs were simulated using piRNA 5′ prefixes. piRNA and siRNA guides were predicted to cleave the mutated transposon sequence using the following rules: ≤6 total mismatches at any position were allowed for 26 nucleotide piRNAs; ≤5 total mismatches, ≤1 mismatch between g2 and g8, and no mismatches at g9, g10, g11 and g13 were allowed for 21 nucleotide siRNAs.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Online content

Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at 10.1038/s41586-023-06257-4.

Supplementary information

Supplementary FiguresSupplementary Figs. 1–5. Uncropped gel source data. Binding affinities and cleavage rates for MILI, MIWI and EfPiwi. FACS gating strategy to purify mouse primary germ cells. Reporting Summary Peer Review File Supplementary Table 1. Cleavage rates reported in this study. Supplementary Tables Supplementary Tables 2–4, 6 and 7. Source data for analyses of in vivo cleavage by MILI and MIWI (number of putative cleavage sites). Mouse strains and sequences of oligonucleotides used in this study. Amount of spike-in RNAs used to prepare small RNA-seq and RNA-seq libraries in this study. Supplementary Table 5. Target, expected 3′ cleavage product and the corresponding 8 nucleotide barcode for data in Extended Data Fig. 5c. Supplementary Data 1Source data for analyses of in vivo cleavage by MILI and MIWI (coordinates and identity of putative cleavage sites).