# Coupled protein synthesis and ribosome-guided piRNA processing on mRNAs

Yu Sun, Ruoqiao Huiyi Wang, Khai Du, Jiang Zhu, Jihong Zheng, Li Huitong

Xie, Amanda Pereira, Chao Zhang, Emiliano Ricci, Xin Zhiguo Li

► To cite this version:

Yu Sun, Ruoqiao Huiyi Wang, Khai Du, Jiang Zhu, Jihong Zheng, et al.. Coupled protein synthesis and ribosome-guided piRNA processing on mRNAs. Nature Communications, 2021, 12 (1), pp.5970. 10.1038/s41467-021-26233-8 . hal-04892075

HAL Id: hal-04892075

https://hal.science/hal-04892075v1

Submitted on 16 Jan 2025

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
ARTICLE

https://doi.org/10.1038/s41467-021-26233-8

OPEN

<figure><img src="image_2.png" /></figure>

Check for updates

# Coupled protein synthesis and ribosome-guided piRNA processing on mRNAs

Yu H. Sun<fnref n="1" />, Ruoqiao Huiyi Wang<fnref n="1" />, Khai Du<fnref n="1" />, Jiang Zhu<fnref n="1" />, Jihong Zheng², Li Huitong Xie<fnref n="1" />, Amanda A. Pereira<fnref n="1" />, Chao Zhang², Emiliano P. Ricci³ & Xin Zhiguo Li<fnref n="1" />✉

PIWI-interacting small RNAs (piRNAs) protect the germline genome and are essential for fertility. piRNAs originate from transposable element (TE) RNAs, long non-coding RNAs, or 3' untranslated regions (3'UTRs) of protein-coding messenger genes, with the last being the least characterized of the three piRNA classes. Here, we demonstrate that the precursors of 3'UTR piRNAs are full-length mRNAs and that post-termination 80S ribosomes guide piRNA production on 3'UTRs in mice and chickens. At the pachytene stage, when other co-translational RNA surveillance pathways are sequestered, piRNA biogenesis degrades mRNAs right after pioneer rounds of translation and fine-tunes protein production from mRNAs. Although 3'UTR piRNA precursor mRNAs code for distinct proteins in mice and chickens, they all harbor embedded TEs and produce piRNAs that cleave TEs. Altogether, we discover a function of the piRNA pathway in fine-tuning protein production and reveal a conserved piRNA biogenesis mechanism that recognizes translating RNAs in amniotes.

<footnote n="1">¹Center for RNA Biology: From Genome to Therapeutics, Department of Biochemistry and Biophysics, University of Rochester Medical Center, Rochester, NY 14642, USA. ²School of Life Sciences and Technology, Tongji University, Shanghai 200092, China. ³Université de Lyon, ENSL, UCBL, INSERM, CNRS, LBMC, 46 Allée d'Italie, Lyon 69007, France. ✉email: Xin_Li@urmc.rochester.edu</footnote>
PIWI-interacting RNAs (piRNAs) are a class of small silencing RNAs that are essential for animal fertility. piRNAs guide a specialized class of Argonaute proteins primarily found in germ cells, known as PIWI proteins¹⁻⁶, to target RNAs via base-pair complementarity. Unlike microRNAs (miRNAs) and small interfering RNAs (siRNAs) that are derived from precursors with double-stranded RNA structure, most piRNAs are derived from long single-stranded transcripts⁷,⁸. The 5'-monophosphate ends (5'P) of piRNAs are formed when piRNA precursors are fragmented into pre-piRNAs that are loaded onto PIWI proteins. This fragmentation process requires MOV10L1, a UPF1-like RNA helicase, and PLD6, an RNA endonuclease on mitochondrial outer membranes⁹⁻¹⁶. The PIWI-loaded pre-piRNAs are further trimmed by a 3'-to-5' exonuclease (PNLDC1 in mice), resulting in a length distribution characteristic of each particular PIWI protein¹⁷⁻²⁶. Finally, HEN1 binds PIWI proteins and adds a 2'-O-methyl group to the 3'-end of mature piRNA²⁷,²⁸. The resulting primary piRNAs can be further amplified via the ping-pong pathway⁷,²⁹, yielding secondary piRNAs.

Compared to our understanding of the piRNA pathway in fruit flies, it remains unclear what marks a transcript for piRNA processing in mammals. In fruit flies, piRNA loci in germ cells are epigenetically marked by heterochromatin-bound factor Rhino³⁰⁻³². piRNA precursors are derived from promoter-independent, unspliced cryptic transcripts in highly repetitive regions that harbor diverse transposable elements (TEs)³³,³⁴. piRNA precursors are then directly channeled for piRNA processing on perinuclear RNA granules that are proximal to nuclear pores³⁵. Maternally deposited piRNAs further recognize the piRNA precursors post-transcriptionally to trigger a cascade of piRNA processing³⁶,³⁷. Unlike germ cells in fruit flies predetermined by maternal deposition, mouse germ cells are induced from somatic cells³⁸,³⁹, leaving parental piRNA deposition unlikely. The majority of mammalian piRNAs come from 5' capped and 3' polyadenylated, long continuous non-coding RNAs (lncRNAs) that are depleted of TEs⁸. The transcription factor A-MYB directs the expression of these lncRNAs from euchromatin regions without identifiable epigenetic markers distinguishing mouse piRNA loci, along with mRNAs that do not produce piRNAs, during pachynema (the pachytene stage of meiosis)⁸,⁴⁰. Despite these piRNA precursors being annotated as lncRNAs, rather than directly channeling to RNA granules for processing, we recently demonstrated that ribosomes translate their upstream open reading frames (uORFs). Thus, transcription and translation of mammalian piRNA precursors utilize conventional machineries, therefore leaving the features that distinguish them from other lncRNAs and mRNAs elusive.

In eukaryotes, translation-dependent RNA quality-control mechanisms identify and degrade aberrant transcripts with premature termination codons (nonsense-mediated decay (NMD)), lacking in-frame stop codons (non-stop decay), or containing elongation inhibitory structures (no-go decay). These quality-control pathways are coupled with translation and facilitate the recycling of ribosomes trapped on faulty transcripts⁴¹⁻⁴⁴. The finding of ribosome-guided piRNA biogenesis downstream of uORFs on lncRNAs⁴⁵ introduces the possibility that aberrant translation intermediates could act as substrates for piRNA biogenesis. It nevertheless remains unclear how 80S ribosomes assemble on the non-coding regions of piRNA precursors, which contain frequent stop codons, and how these ribosomes escape the translation-dependent surveillance mechanisms that recycle ribosomes and degrade RNAs⁴⁴. The detection of ribosomes downstream of uORFs on lncRNA piRNA precursors is reminiscent of the re-initiation of ribosomes downstream of uORFs present in the 5' untranslated regions (5'UTRs) of canonical

mRNAs when initiation factors have not yet dissociated from the 40S subunit⁴⁶. Given that short uORF length is an important factor for translation of the main ORFs on mRNAs⁴⁶, which explains why eukaryotic mRNAs are generally monocistronic, the presence of long ORF on piRNA precursors should inhibit ribosome-guided piRNA biogenesis. It is thus mechanistically essential for both piRNA biogenesis and translational regulation to test whether ribosome-guided biogenesis can only occur after the translation of a short ORF.

To test whether a short uORF is required for ribosome-guided piRNA biogenesis, we studied a specific class of piRNAs, namely, 3'UTR piRNAs. In mice, piRNAs are divided into three major classes based on their origin⁴⁷,⁴⁸: (i) piRNAs from TEs (TE piRNAs), (ii) intergenic piRNAs derived from lncRNAs (pachytene piRNAs), and (iii) genic piRNAs that map to 3'UTRs of protein-coding regions in the sense orientation (3'UTR piRNAs). TE piRNAs protect the animal germline genome from TEs and are essential for animal fertility, a function that is evolutionarily conserved in bilateral animals¹⁻³,⁵,⁴⁹,⁵⁰. Pachytene piRNAs have only been reported in mammals thus far and have been shown to trigger the decay of mRNAs required for sperm formation⁵¹⁻⁵⁴. 3' UTR piRNAs have been detected in fruit flies, frogs, and diverse mammalian species⁵⁵⁻⁵⁷, but we currently do not know their function(s), how their production is regulated, nor whether their precursors are full-length mRNAs or cryptic transcripts corresponding exclusively to the 3'UTR⁵⁶. The lack of understanding of 3'UTR piRNAs hindered our efforts to identify either common features that mark such transcripts for piRNA biogenesis and/or machinery that sorts diverse RNAs for piRNA biogenesis.

Here we characterize the biogenesis of 3'UTR piRNAs in mice, demonstrating that their precursors are full-length protein-coding mRNAs. We further show that piRNA biogenesis from these precursors is coupled with efficient translation and that ribosomes guide piRNA precursor fragmentation on mRNA 3'UTRs. We demonstrate that this tight coupling of ribosome binding and piRNA biogenesis fine-tunes protein synthesis from mRNAs. Ribosome-guided piRNA processing occurs at the meiotic stage when ribosome recycling factors and NMD are temporally inhibited. Lastly, we demonstrate that 3'UTR ribosome-guided piRNA processing also occurs in chickens. Although 3'UTR piRNAs are derived from distinct sets of genes in mice and chickens, we found the presence of TE sequences to be a shared feature that serves to produce anti-sense TE piRNAs that cleave TEs post-transcriptionally, indicating that TE suppression is a conserved evolutionary force driving 3'UTR piRNA production. Taken together with our previous studies, we find that a general and conserved piRNA biogenesis pathway recognizes translating RNAs regardless of their ORF length.

## Results

**3'UTR piRNAs are produced from full-length mRNAs.** To test whether 3'UTR piRNAs are derived from full-length mRNAs or cryptic transcripts with alternative transcription start sites (TSSs), we set out to identify the structure of the precursor RNAs by blocking their transcription and inhibiting the post-transcriptional processing of 3'UTR piRNAs. We previously defined a group of 3'UTR piRNAs increasing from 12.5 days postpartum (dpp) to 17.5 dpp in mice⁸, when spermatocytes enter pachynema (Supplementary Fig. 1a). Thus, we hypothesized that these 3'UTR piRNA precursors may be controlled by the transcription factor A-MYB, which also promotes the synthesis of pachytene piRNA precursors⁸. Chromatin immunoprecipitation sequencing (ChIP-seq) revealed that A-MYB-binding sites are far from the 3'UTRs of mRNAs that come from 3'UTR piRNA-
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

**Fig. 1** mRNA precursors of 3'UTR piRNA are regulated by A-MYB. **a** The distance from the annotated transcription start site of each mRNA that comes from 3'UTR piRNA-producing loci (uppl) (left, *n* = 30) and from the 5'-ends of uppl mRNA 3'UTRs (right, *n* = 30) to the nearest A-MYB ChIP-seq peak. **bp:** base pair, **kbp:** kilobase pair. **b** Top, the MYB motif from the mouse UniPROBE database. Bottom, MEME-identified sequence motif in the A-MYB ChIP-seq peaks near the uppl mRNA transcription state sites. **E-value** computed by MEME measures the statistical significance of the motif. The information content is measured in bits and, in the case of DNA sequences, ranges from 0 to 2 bits. A position in the motif at which all nucleotides occur with equal probability has an information content of 0 bits, while a position at which only a single nucleotide can occur has an information content of 2 bits. **c** Boxplots showing the abundance of uppl mRNA (left) and 3'UTR piRNA (right) per uppl gene in A-Myb mutants (red) and their heterozygous littermates (blue) in testes at 14.5 and 17.5 dpp. **Dpp:** days postpartum, **tpm:** transcript per million, **rpkm:** reads per kilobase million. **uppl n = 30.** Box plots show the 25th and 75th percentiles, whiskers represent the 5th and 95th percentiles, and midlines show median values. **d** Aggregated data for 3'UTR piRNA abundance (upper) and uppl mRNA abundance (bottom) on uppl mRNAs (10% trimmed mean) from adult *Mov10l1*<sup>CKO/Δ</sup> testes (left) and *Mov10l1*<sup>CKO/Δ</sup> Neurog3-*cre* (right). Signals are aligned to the transcriptional start site (5'cap) and site of polyadenylation (3'PolyA) and further aligned to the ORF regions. Dotted lines show the translation start codon (ORF Start) and stop codon (ORF End). The x-axis shows the median length of these regions. **Ppm:** parts per million.

producing loci (uppl) (median distance >14 kb, Fig. 1a) but close to the TSSs of uppl mRNAs (median distance 81.5 nucleotides, nts; Fig. 1a). These sites are also enriched for MYB consensus sequences (Fig. 1b). To verify that A-MYB regulates uppl mRNAs, we studied mice carrying an *A-Myb* hypomorphic mutation<sup>58</sup>, which displayed significantly decreased levels of uppl mRNAs based on RNA sequencing (RNA-seq) data from 14.5 dpp (*p* = 4.7 × 10<sup>-2</sup>) and 17.5 dpp testes (*p* = 1.0 × 10<sup>-3</sup>, Fig. 1c, left, and Supplementary Fig. 1b). Therefore, as it does for pachytene piRNA precursors, A-MYB directly regulates the transcription of uppl mRNAs.

To test whether uppl mRNAs are 3'UTR piRNA precursors, we examined the effect of the loss of transcription of uppl mRNAs on 3'UTR piRNA biogenesis in *A-Myb* mutant mice. Like uppl mRNAs, 3'UTR piRNA abundance significantly decreased in mutant testes (Fig. 1c, right, *p* ≤ 4.8 × 10<sup>-10</sup>, and Supplementary Fig. 1b). 3'UTR piRNA depletion associated with uppl mRNA transcription loss could reflect an indirect effect of the meiotic arrest caused by the *A-Myb* mutant<sup>58</sup> or the lack of piRNA biogenic factors, the transcriptional activation of which requires A-MYB<sup>8</sup>. To test these possibilities, we blocked piRNA processing but not meiosis by conditionally knocking out (CKO) *Mov10l1* in
spermatocytes (Supplementary Fig. 1c)¹¹. In the following analyses, we defined a group of non-piRNA-producing mRNAs that display similar expression dynamics in mouse testes as uppl mRNAs to be the developmentally matched-control mRNAs for comparison (Supplementary Data 1). These control mRNAs remained unchanged in *Mov10l1* CKO mutants compared to control littermates (Supplementary Fig. 1d, e). We found that 3' UTR piRNAs are depleted in *Mov10l1* CKO mutants (Fig. 1d, upper), whereas uppl mRNAs showed significant accumulation in *Mov10l1* CKO mutants (Fig. 1d and Supplementary Fig. 1e, *p* = 4.1 × 10⁻¹¹). We ruled out the possibility that the increased steady-state of uppl mRNAs is due to their transcriptional activation in *Mov10l1* CKO mutants (Supplementary Fig. 1f, *p* = 0.13). Together with the cross-linking immunoprecipitation data (CLIP)¹³ that indicate MOV10L1 specifically binds to the entire length of uppl mRNAs (Supplementary Fig. 1g), our results suggest that the increased steady-state abundance of uppl mRNAs is due to the lack of 3'UTR piRNA biogenesis. Overall, our results indicate that 3'UTR piRNAs come from full-length uppl mRNAs, not from isoforms derived from uppl 3'UTRs.

**3'UTR piRNAs are derived from processed mRNAs.** A recent study suggests that unspliced transcripts are substrates for piRNA processing⁵⁹. However, our data indicated that 3'UTR piRNAs were disproportionately produced after intron removal: >99% of piRNAs mapped to exons, compared to <0.1% that mapped to introns (17,200 ppm unique mapping 3'UTR piRNA reads from adult testes). After correcting for the length of exons and introns in piRNA-producing primary transcripts, exon-derived piRNAs were enriched by 630-fold compared to intron-derived piRNAs (Supplementary Fig. 2a). We detected piRNAs that failed to map to the genome but instead mapped to exon-exon junctions (Fig. 2a), further indicating that piRNAs were produced after intron removal and exon-exon joining. In uppl mRNA genes, 98% of the introns (207 out of 211) contained canonical GT-AG splice sites, not significantly different from the introns in control mRNA genes (326 out of 329 mRNA loci, χ² test, *p* = 0.55), suggesting that uppl mRNAs were spliced conventionally. Furthermore, the density of piRNAs falls off sharply after the 3'end of the transcript, i.e., the site of polyadenylation (Supplementary Fig. 2a). Taken together with our previous piRNA precursor studies based on a combination of cap analysis of gene expression (CAGE) and the polyadenylation site sequencing (PAS-Seq)⁸, we conclude that 3' UTR piRNAs were produced after their precursor transcripts are fully processed—capped, spliced, and polyadenylated.

**uppl mRNA precursors harbor extensive 3'UTRs.** We used our recent reconstruction of the mouse testis transcriptome⁶⁰, combining CAGE, PAS-seq, and single-molecule long-read sequencing, to characterize the transcript structure of uppl mRNAs in comparison to that of developmentally matched and non-piRNA-producing control mRNAs. On average, mature uppl transcripts (*p* < 2.2 × 10⁻¹⁶), but not primary transcripts (*p* = 0.18), were longer than control mRNAs (Fig. 2b and Supplementary Fig. 2b). Given similar numbers of splicing events (no significant difference in the number of introns, *p* = 0.23, Supplementary Fig. 2c), uppl mRNAs harbored longer exons with a median length of 153 nts (Supplementary Fig. 2d; *p* = 3.1 × 10⁻¹⁴) versus the control mRNAs with a median exon length of 109 nts, similar to that of typical vertebrate mRNAs⁶¹. Since exon length distribution is associated with the relative position in the transcripts, we further separated exons into first, middle, and last exons. The last exons for the control mRNAs and uppl mRNAs are the longest among the three categories, whereas the middle exons are the shortest (Fig. 2c). While there is no significant difference in the length

distribution of the middle exons (*p* = 0.65), both the first and the last exons of uppl mRNAs are significantly longer than those of the control mRNAs (*p* ≤ 4.2 × 10⁻⁸, Fig. 2c). Long first exons were recently reported to be a conserved feature for pachytene piRNA precursors⁴⁰, suggesting that the unique feature of piRNA precursors can be traced back to exon-intron structures and the selection of poly-A cleavage sites.

Since 5'UTRs and 3'UTRs are typically localized to the first and last exons, respectively, we further separated the spliced transcripts according to their longest ORFs. Compared to mRNAs in other tissues, testis mRNAs are reported to have shorter 3'UTRs⁶²⁻⁶⁵. Consistent with these reports, the 3'UTRs of the control mRNAs have a median length of 392 nts. However, against the trend in testes to produce mRNAs with proximal polyadenylation sites, the 3'UTRs of the uppl mRNAs are significantly longer than those of control mRNAs, with a median length of 4583 nts (Fig. 2d; *p* = 1.5 × 10⁻¹¹). While the lengths of the ORFs are similar (Fig. 2d, 1315.5 nts versus 945 nts; *p* > 0.01), the 5'UTRs of uppl mRNAs are also significantly longer than those of the control mRNAs (Fig. 2d, 346 nts versus 77 nts; *p* = 3.3 × 10⁻⁷). Therefore, regulatory sequences (UTRs), rather than ORFs, contribute to the majority of the length difference between the spliced transcripts of uppl mRNAs and control mRNAs (Fig. 2b) and render more diverse translational regulations on uppl mRNAs.

**3'UTR ribosomes guide piRNA formation.** The biogenesis of piRNAs from mRNA 3'UTRs suggests that cellular translation machinery may be used by piRNA biogenesis machinery to distinguish between 3'UTRs and ORFs. To test whether uppl mRNAs are translated, we performed ribosome profiling (Ribo-seq) in which RNA fragments protected from RNase A and T1 digestion are isolated from 80S fractions and sequenced⁴⁵. We found that ribosomeprotected fragments (RPFs) from uppl mRNA ORFs displayed a three-nucleotide periodicity (Fig. 2e), indicating that elongating ribosomes translate these ORFs. We calculated their translational efficiency (RPF reads normalized to transcript abundance) and found that uppl mRNAs have a slightly lower translational efficiency compared to control mRNAs (Fig. 2f, 59% of the median translational efficiency of control mRNAs, *p* = 0.037), as expected given that their longer 5' UTRs (Fig. 2c) will take longer time to scan through before the next 40S subunits start⁶⁶ and may also contain more upstream initiation sites and/or have more complex secondary structures that can reduce translational efficiency⁶⁷. We tested codon usage, another key determination factor for translational efficiency⁶⁸, and found that uppl mRNAs and control mRNAs have a similar codon usage (Supplementary Fig. 2e, f)⁶⁹⁻⁷¹. Overall, uppl mRNAs are efficiently translated on their main ORFs without obvious signs of aberrant initiation or elongation distinguishing them from control mRNAs. Taken together, uppl mRNAs function to produce both piRNAs and proteins.

To test whether ribosomes also guide piRNA biogenesis from uppl mRNAs as they do from lncRNAs⁴⁵, we identified RPF signals coming from uppl mRNA 3'UTRs (Fig. 3a). 3'UTR RPF signals were not seen in the control mRNAs (Supplementary Fig. 3a, b, *p* = 1.1 × 10⁻⁸), indicating that the distribution of RPFs on uppl mRNA 3'UTRs is RNA-specific, rather than the characteristic of typical testicular mRNAs. The in vitro RNase digestion used to obtain RPFs removed >98% of the mature 3' UTR piRNAs in testis lysates (Supplementary Fig. 3c). We performed Ribo-seq on adult *RiboTag* mice after activating the expression of hemagglutinin (HA)-tagged RPL22 (a ribosomal large subunit protein) in germ cells⁴⁵ and affinity-purified RPFs before sequencing. After IP, sequences from mitochondrial
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

**Fig. 2** Translation and piRNA processing are coupled on uppl mRNAs. **a** Top, aggregation plots of piRNA reads surrounding the 5' and 3' splice sites of uppl mRNAs from the adult mouse testes. Bottom, the signal was calculated for non-genome matching piRNA reads mapping to exon-exon junction sequences. **b** Boxplots showing spliced RNA transcript length distributions. Control mRNA n = 43, uppl mRNA n = 30. **c** Boxplots showing exon length distributions. Control mRNA first or last exon n = 43, uppl mRNA first or last exon n = 30. Control mRNA middle exon n = 145, uppl mRNA middle exon n = 69. **d** Boxplots showing ORF and UTR length distributions. Control mRNA n = 43, uppl mRNA n = 30. **e** Discrete Fourier transformation of the distance spectrum of 5'-ends of RPFs across ORFs (gold) and 3'UTRs (green) of uppl mRNAs in adult wild-type testes. **f** Boxplots showing translational efficiency (RPF abundance divided by mRNA abundance) in adult wild-type testes. Control mRNA n = 43, uppl mRNA n = 30. Box plots in (b-d and f) show the 25th and 75th percentiles, whiskers represent the 5th and 95th percentiles, and midlines show median values.

coding transcripts were also depleted (Supplementary Fig. 3d) given that they are translated by untagged 55S mitochondrial ribosomes. In contrast, uppl mRNA ORF and 3'UTR sequences were retained, similar to RPFs from control mRNAs (Supplementary Fig. 3d). Overall, these results indicate that the 3'UTR RPF signals are bona fide ribosomal footprints.

To test whether ribosome-bound 3'UTRs are recognized as precursors for 3'UTR piRNA production, we performed partial correlation analyses⁷² between RPF abundance and piRNA abundance in uppl 3'UTRs while controlling for the abundance of uppl mRNAs, as measured by RNA-seq. These analyses can distinguish a biogenic relationship between RPFs and piRNAs or the independent correlations of RPFs and piRNAs with their uppl mRNA precursors. We found that the abundance of RPFs and piRNAs at 3'UTRs are directly correlated with each other (Fig. 3b, right, r = 0.69, *p* = 2.9 × 10⁻⁵). Thus, uppl 3'UTRs bound by ribosomes are processed into piRNAs.

To test whether ribosomes guide the degradation of 3'UTRs for piRNA production, we analyzed the position of the RPFs on uppl mRNA 3'UTRs. We aligned piRNAs to RPFs and plotted the 5'-ends of RPFs that overlapped with piRNAs (Fig. 3c). For RPFs, their first nucleotide overlapped with the first nucleotide of

piRNAs significantly more than with nucleotides residing 50 nts upstream or downstream (Fig. 3c, right, Z score = 36 ± 1, Z scores indicate how many standard deviations an element is from the mean; Z score > 3.3 corresponds to *p* < 0.01). We ruled out the possibility that this overlap could occur by chance or be due to ligation bias (Supplementary Fig. 3e, f). Consistent with their 5'-overlap, we found that the 5'-ends of RPFs from 3'UTRs displayed a uridine bias at the 5'-most position (1U) (Fig. 4a), reminiscent of the 1U bias in piRNAs. These results indicate that 3'UTR ribosomes dwell at the sites that represent future piRNAs.

Given that in vitro RNA digestion (RNase T1&A) used to obtain RPFs does not yield a 1U bias, we tested the possibility that these RPFs are processed in vivo by piRNA processing machinery. We modified the conventional Ribo-seq procedure (which detects both 5'P and 5'OH RPFs) to specifically capture 5' P RPF and 5'OH RPF species separately (Supplementary Fig. 4a). 5'P species are products of in vivo enzymatic cleavage, whereas the 5'-hydroxyl (5'OH) species mainly arise from in vitro RNase treatment (RNase T1&A in our procedure). To prevent mature piRNA contamination, we used the RPFs from affinity-purified 80S ribosomes for library construction (Supplementary Fig. 4a). The 5'OH RPFs showed the expected in vitro digestion signature,
<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

**Fig. 3 3'UTR ribosomes guide piRNA 5'-end formation.** a Normalized reads of anti-HA-immunoprecipitated RPFs and piRNAs mapping to four representative uppl genes in adult testes. b Residual plot of partial correlation between piRNA and RPF abundance conditional on RNA-seq abundance at the uppl ORFs (left) and 3'UTRs (right) in adult testes. Linear regression was performed between piRNA and RNA-seq abundance and between RPF and RNA-seq abundance separately, and residuals were plotted with smooth parameter "method = lm". Data are presented as mean values ± standard deviation shown as error bands. c Distance spectrum of 5'-ends of RPFs (26-32 nts) from ORFs (left) and from 3'UTRs (right) that overlap piRNAs in adult testes (sample size n = 3 independent biological replicates). Data are mean ± standard deviation.

even at uppl 3'UTRs (Supplementary Fig. 4b). However, the 5'P RPFs are >32-fold enriched relative to 5'OH RPFs at uppl 3'UTRs in comparison to uppl ORFs (Fig. 4b, *p* = 1.0 × 10⁻⁶). This indicates that RPFs predominantly present with 5'P ends in 3' UTRs at steady state, suggesting that in vivo cleavage occurs efficiently on ribosome-bound 3'UTRs.

Consistent with these sites representing hot spots for efficient in vivo cleavage, 3'UTR ribosomes are significantly enriched in the monosome fractions relative to polysome fractions in comparison to ORFs of the control mRNAs (Fig. 4c, *p* = 2.5 × 10⁻¹⁰), as measured by Ribo-seq performed on purified monosome and polysome fractions as we reported previously⁴⁵. The 1U signature of 5'P RPFs occurs specifically at the uppl 3' UTRs but not at the ORFs of uppl RNA or the ORFs of control mRNAs (Supplementary Fig. 4b), indicating that 3'UTR RPFs are cleaved by the piRNA processing machinery. Altogether, piRNA

processing machinery generates the 5'P in vivo cleavage products with a ribosome bound at their 5' extremities, and these 5'P ends become the 5'-ends of future piRNAs.

**3'UTR piRNA biogenesis is coupled with upstream translation.**

The translation of uppl mRNAs could either be coupled or uncoupled with piRNA biogenesis. If coupled, there would be a limited time window for translation to occur before the mRNA can be cleaved to generate piRNAs. If uncoupled, uppl mRNAs should undergo multiple rounds of translation. In mammals, two types of cap-binding proteins (CBPs) participate in protein synthesis⁷³. The largely nuclear CBPs 80 and 20 (CBP80/20), which mediate the pioneer round of translation once newly synthesized mRNAs reach the cytoplasm⁷⁴,⁷⁵, are ultimately replaced by the major cytoplasmic CBP eukaryotic translation
<figure><img src="image_18.png" /><figcaption>**Fig. 4** Biphasic piRNA biogenesis. **a** Sequence logos depicting nucleotide (nt) bias at 5'-ends and 1 nt upstream of 5'-ends. Top to bottom: RPF species, anti-HA-immunoprecipitated 5'P RNA species, and piRNAs in adult testes, which map to mRNAs, uppl mRNA ORFs, and 3'UTRs, respectively. **b** Boxplots of the ratios of 5'P RPF versus 5'OH RPFs in adult wild-type testes. Control mRNA *n* = 43, uppl mRNA *n* = 30. **c** Boxplots of the ratios of RPF abundance in monosome fractions versus those in polysome fractions in adult wild-type testes. Control mRNA *n* = 43, uppl mRNA *n* = 30. **d** Boxplots of the ratios of RNA abundance in CBP80 IP versus eIF4E IP in adult wild-type (upper) and *Mov10l1* CKO mutant (lower) testes. Control mRNA *n* = 43, uppl mRNA *n* = 30. **e** Aggregated data for RPF abundance from untreated adult testes (top) and from harringtonine-treated adult testes (bottom) across 5'UTRs, ORFs, and 3' UTRs of the uppl mRNAs (left) and control mRNAs (right). The x-axis shows the median length of these regions, and the y-axis represents the 10% trimmed mean of relative abundance. **f** piRNA abundance on uppl mRNAs (10% trimmed mean) in adult testes. Control (left) and *Tdrd5* CKO (right). Box plots in (b–d) show the 25th and 75th percentiles, whiskers represent the 5th and 95th percentiles, and midlines show median values.</figcaption></figure>

initiation factor 4E (eIF4E)⁷⁶. To detect the nature of these CBPs during the translation of mRNAs that generate piRNAs, we immunoprecipitated (IP) CBP80-bound RNAs and eIF4E-bound RNAs and subjected them to RNA-seq. We calculated RNA abundance after CBP80 IP versus eIF4E IP and found that uppl mRNAs are significantly enriched in the CBP80 IP in comparison to control mRNAs (Fig. 4d, upper, *p* = 1.3 × 10⁻⁵). To test whether the enrichment of uppl mRNA with CBP80 is due to piRNA processing, we repeated the CBP80 IPs and eIF4E IPs using *Mov10l1* CKO mutants. The enrichment of uppl mRNAs in CBP80 IPs was significantly reduced in *Mov10l1* CKO mutants compared to the wild type (Fig. 4d, *p* = 3.7 × 10⁻⁴), while the distribution of control mRNAs was unaffected (*p* = 0.32). Our results indicate that our failure to detect appreciable translation of eIF4E-bound uppl mRNAs is due to their processing into piRNAs while they are CBP80-bound, suggesting that uppl mRNA

translation and piRNA biogenesis are temporally coupled and carried out using the same uppl RNA molecule.

The coupling of piRNA processing and translation indicates that the uppl mRNAs are targeted by the piRNA processing machinery immediately after or while they are being translated. To investigate how this coupling is achieved, we tested whether 3' UTR ribosome binding requires upstream initiation near the 5'-cap. We treated testes with harringtonine⁴⁵, which blocks early elongation by binding to the A site of newly assembled ribosomes⁷⁷. Given that uppl 3'UTRs do not harbor long ORFs (median putative ORF length is 42 nts, Supplementary Fig. 4c), 3' UTR ribosomes should accumulate and be resistant to harringtonine treatment if there is internal initiation. After a 2-h treatment, the control mRNAs showed a significant reduction in elongating ribosomes on ORFs (*p* = 8.7 × 10⁻¹², Fig. 4e and Supplementary Fig. 4d), compared to the initiating ribosomes accumulated
around the start codons. On uppl mRNAs, ribosomes accumulated at the translation start sites, similar to control mRNA ($p = 0.17$). In comparison to the ribosomes at translation start sites, ribosomes at ORFs ($p = 4.4 \times 10^{-11}$) and 3'UTRs ($p = 1.2 \times 10^{-11}$) were substantially reduced (Fig. 4e and Supplementary Fig. 4d). Thus, 3'UTR ribosomes migrate from upstream long ORFs rather than loading internally. Unlike the conventional initiation mechanisms on mRNA ORFs after uORF that require a short uORF length, our study suggests compromised post-termination recycling, which underlies the coupling between translation at ORFs and piRNA processing at 3'UTRs.

**Biphasic biogenesis before and after the stop codon.** uppl ORFs also produce authentic piRNAs with a 1U bias (Fig. 4a), although piRNAs derived from 3'UTRs are >31-fold more abundant than they are in ORFs, which equates to an 8-fold difference when normalized to 3'UTR and ORF length, respectively (Supplementary Fig. 4e). To understand why significantly different amounts of piRNAs are detected from ORFs and 3'UTRs ($p = 1.4 \times 10^{-7}$, Supplementary Fig. 4e) and how this is related to ribosomes, we analyzed uppl ORF RPFs. Unlike ORF piRNAs with a 1U bias, ORF RPFs exhibit an in vitro digestion signature without a 1U bias (Fig. 4a). The ratios of 5'P RPFs to 5'OH RPFs from uppl ORFs are comparable to those of RPFs from the control mRNA ORFs (Fig. 4b, $p = 0.71$), suggesting that RPFs present at uppl ORFs predominantly harbor 5'OH ends at the steady state. RPFs from uppl ORFs have a significantly less pronounced 5' overlap with mature piRNAs than 3'UTR RPFs (Fig. 3c, Student's t test, $p = 1.4 \times 10^{-4}$). Consistent with the lack of in vivo processing of ORF RPFs, we found that ORF RPFs have a similar distribution between monosome and polysome fractions when compared to RPFs from the control mRNA ORFs (Fig. 4c, $p = 0.45$). Furthermore, the abundance of RPFs and piRNAs at ORFs do not correlate with each other when controlling for the abundance of uppl mRNAs (Fig. 3b, left, $r = -0.3$, $p = 0.1$), indicating a lack of a biogenic relationship between ribosome-bound ORFs and piRNAs. Therefore, our results indicate that ribosome-guided piRNA processing occurs at uppl 3'UTRs but not at uppl ORFs.

The lack of a biogenic relationship between RPFs and piRNAs at uppl ORFs could be due to: (1) a different mechanism of piRNA processing where only a translationally suppressed subpopulation of uppl mRNA ORFs is processed into piRNAs, with the majority protected from processing; or (2) the mRNA ORFs are processed but do not generate piRNAs efficiently. To distinguish between these two possibilities, we monitored the fragmentation process on ribosomes. Using degradome sequencing (degradome-seq) of long RNAs (>200 nt) from affinity-purified ribosomes, we captured ribosome-bound 5'P RNAs. We found that the abundance of ribosome-bound 5'P RNAs from 3' UTRs was not higher than that from ORFs ($p = 0.3$, Supplementary Fig. 4e), arguing against the possibility that the majority of ORFs were protected from processing. We then tested whether these 5'P RNAs have their 5'-ends aligned with the 5'-ends of piRNAs, which would suggest that the cleavage process forms the piRNA 5'-ends. We found that the 5'P RNAs from uppl ORFs had a significantly lower 5' overlap (Z score = 3.2 ± 0.4) with the 5'-ends of uppl ORF piRNAs compared to the 5' overlap between 5'P RNAs and piRNAs at uppl 3'UTRs (Z score = 18 ± 1, Student's t test, $p = 1.8 \times 10^{-5}$, Supplementary Fig. 4f). Consistent with this, neither ribosome-bound 5'P RNAs nor RPFs from uppl ORFs have the 1U bias, unlike that of the 5'P RNAs and RPFs from uppl 3'UTRs (Fig. 4a). Thus, although uppl ORFs are cleaved, the lack of 5' overlap between the cleavage products and uppl ORF piRNAs indicates that the cleavage products are inefficiently processed into piRNAs.

To determine whether auxiliary factors facilitate efficient processing of piRNAs at uppl 3'UTRs, we tested for TDRD5 (tudor domain containing 5) function, the disruption of which impacts pachytene piRNA production$^{45,78}$. We found that piRNAs derived from uppl ORFs were still produced, but 3' UTR piRNAs were depleted in *Tdrd5* mutants (Fig. 4f and Supplementary Fig. 4f), indicating that uppl ORF piRNAs do not require TDRD5 for their production, whereas 3'UTR piRNAs do. Therefore, 3'UTR piRNAs derived from uppl 3'UTRs and ORFs have linked but distinct biogenic requirements. Taken together, our results indicate that efficient piRNA processing at uppl 3' UTRs requires both ribosomes and TDRD5 and that post-transcriptional processing differences at 3'UTRs and ORFs explain why even though the entire length of uppl mRNAs generates piRNAs, >96% of uppl piRNAs are derived from uppl 3' UTRs. These results are similar to the biphasic piRNA biogenic mechanism we identified previously in pachytene piRNAs from lncRNA piRNA precursors$^{45}$. Therefore, combined with our previous work on lncRNAs$^{45}$, our study reveals a general ribosome-guided mechanism by which piRNA precursors, regardless of their source, are converted into piRNA sequences (Supplementary Fig. 4g).

**Inhibition of co-translational surveillance pathways.** Ribosomes that failed to recycle at the 3'UTRs of mRNAs should be rescued/ resolved by ribosome recycling factors such as PELOTA (the mouse homolog of yeast DOM34)$^{79-84}$. To understand how 3' UTR ribosomes avoid recycling by PELOTA, we performed immunostaining on squashed spermatocytes and spermatids. We found that PELOTA localized to the nuclei of pachytene spermatocytes, but not to the nuclei of round spermatids (Fig. 5a). Similar nuclear localization, specifically at the pachytene stage, has been reported for other translational control proteins$^{85}$. Considering that PELOTA is the major player in no-go decay$^{86}$ and no-stop decay$^{87}$, the sequestration of ribosome recycling factors in the nucleus suggests that their associated mRNA decay pathways are inhibited at pachynema.

Although uppl mRNAs that have extensive 3'UTRs (Fig. 2d) should represent conventional substrates to elicit NMD$^{88}$, the spreading of ribosomes on 3'UTRs should inhibit the NMD as demonstrated by polycistronic viruses or inhibition of ribosome recycling factors$^{89,90}$. To test this idea, we performed IP in testis lysates using anti-UPF1 and anti-phosphorylated UPF1, followed by RNA-seq. UPF1 plays a role in NMD target recognition and elimination, and the phosphorylation of UPF1 is required to activate NMD$^{74,91}$. We found that uppl mRNAs are enriched in UPF1 IP (Fig. 5b, $p = 3.9 \times 10^{-8}$), but this uppl mRNA-associated UPF1 did not show increased levels of phosphorylation compared to control mRNAs (Fig. 4b, lower, $p = 0.15$), indicating the NMD is not activated on the uppl mRNAs given that they appear to be NMD substrates. The NMD pathway is further repressed globally at the pachynema in testes$^{92}$. The lack of NMD activity, in conjunction with the nuclear localization of PELOTA/DOM34 at the pachytene stage, may be a prerequisite for piRNA biogenesis during normal development. Taken together, ribosomes binding to the piRNA precursor uppl mRNAs temporally stagger with other translation-dependent mRNA decay pathways, allowing massive piRNA production from mRNAs and lncRNA precursors during the pachytene stage with an abundance of 3.8–8.4 million piRNA molecules in each spermatocyte$^{25}$ constituting >95% of the total piRNAs in adult mouse testes$^{8}$.

**piRNA biogenesis fine-tunes protein production.** Given that other translation-dependent mRNA decay targets endogenous mRNAs to fine-tune protein production$^{93}$, we tested whether
<figure><img src="image_19.png" /><figcaption>**Fig. 5** 3′UTR piRNA biogenesis as a special co-translational mRNA decay pathway. **a** Immunolabeling of squashed pachytene spermatocytes (left) and round spermatids (right). Three examples of each stage were shown using anti-PELOTA, DAPI, and merged. Experiments were repeated at least three times independently with similar results. Scale bar, 5 µm. **b** Boxplots of the ratios of RNA abundance in UPF1 IP (upper) and phosphorylated UPF1 IP mutant (lower) versus their corresponding input in adult wild-type testes. Control mRNA *n* = 43, uppl mRNA *n* = 30. **c** Boxplots of the changes of protein abundance coded per gene in *Mov10l1* mutants (sample size *n* = 5 independent biological replicates) compared to littermate controls in adult testes. Control mRNA *n* = 43, uppl mRNA *n* = 30. **d** AGO2 western blot (top) and TUBULIN western blot (bottom) of adult testis lysates from serial dilutions of adult wild-type lysates (left) and four biological replicates of *Mov10l1*CKO/Δ Neurog3-cre (middle) and *Mov10l1*CKO/Δ (right). Abundance of AGO2 and TUBULIN was quantified using ImageJ¹⁵⁸ and then fitted to the standard curve of the wild-type testis lysates. Relative abundance of AGO2 were calculated as the abundance of AGO2 normalized by the abundance of TUBULIN (sample size *n* = 4 independent biological samples). Data are mean ± standard deviation; \* *p* = 0.048, one-side Student's *t* test. Each data point was overlaid as dot plots. **e** Boxplots of the 5'P decay intermediates per mRNA in adult wild-type testes. AGO2 target transcripts *n* = 41. Box plots in (b, c, e) show the 25th and 75th percentiles, whiskers represent the 5th and 95th percentiles, and midlines show median values.</figcaption></figure>

the piRNA biogenesis pathway impacts protein production from uppl mRNAs. We performed quantitative mass spectroscopy using *Mov10l1* CKO mutants and their littermate controls. We observed a significant increase in steady-state protein levels from uppl mRNAs, with a median increase of 13%, in comparison to the control mRNAs (*p* = 3.2 × 10⁻², Fig. 5c). These results indicate that piRNA biogenesis regulates the abundance of a subset of proteins, including AGO2 (encoded by *pi-Ago2/Eif2c2.1*), which plays a central role in small RNA function, DDX19B (encoded by *pi-Ddx19b.1*), the homolog of yeast *Dbp5*), which controls mRNA export⁹⁴,⁹⁵, and the translational repressor 4E-BP1 (encoded by *pi-Eif4ebp2.1*)⁹⁶. Among the 16 out of 30 uppl genes with reported mutant phenotypes⁹⁷, half of them have been demonstrated to be essential, 5 of which (*Ddx19b*, *Strbp*, *Asb1*, *Ip6k1*, and *Nr2c2*) exhibit impaired male fertility⁹⁸⁻¹⁰² and 3 of which (*Ago2*,

*Exoc8*, and *Ipmk*) display embryonic lethality⁹⁸,¹⁰³,¹⁰⁴. This is significantly different (*χ*² test, *p* = 1.3 × 10⁻⁴) from the reported phenotypes in all mouse mutants wherein 7% and 30% of mutants display fertility defects and embryonic lethality, respectively¹⁰⁵. Thus, 3′UTR piRNA processing affects the level of proteins from essential genes.

To test the impact of protein amount changes due to 3′UTR piRNA biogenesis, we chose to focus on the miRNA function of AGO2, whose abundance significantly increased by 115% according to western blot quantification (Student's *t* test, *n* = 4, *p* < 0.05, Fig. 5d) in total testes of *Mov10l1* CKO mutants where piRNA biogenesis is blocked. We focused on miRNAs that are regulated by A-MYB to ensure that the targeted events occur at the same developmental stage when 3′UTR piRNAs are produced. We identified six miRNAs (mmu-mir-449a,c, mmu-mir-34b,c, mmu-mir-184, and mmu-mir-191) whose expression are
significantly reduced in the *A-Myb* mutant (*q* value < 0.05) and have an A-MYB binding peak nearby (within 1 kb) (Supplementary Fig. 5a, b). Given that miRNAs can be loaded into all four AGO proteins (AGO1–4), to distinguish the function of AGO2-bound miRNAs, we took advantage of the fact that only AGO2 has slicer activity¹⁰⁶ and searched for the targets of miRNA-guided AGO2-mediated cleavage with the following criteria: (1) the cleavage products are detected in a degradome-seq library that captures the 5'P species in adult testis; and (2) the target RNAs are directly regulated by A-MYB to ensure that the targeted events occur at the same cell type and same developmental stage when 3'UTR piRNAs are produced. We detected 41 target transcripts meeting these criteria. The number of cleavage events in these transcripts, as measured by degradome-seq, significantly increased with a median increase of 146% in *Mov10l1* CKO mutants (Fig. 5e, *p* = 2.5 × 10⁻⁴). In sum, the increased protein level of AGO2 in *Mov10l1* CKO mutants leads to increased activity of miRNAs expressed during pachynema.

Among the 17 out of 41 target transcripts of miRNA-guided AGO2-mediated cleavage with reported mutant phenotypes⁹⁷, 14 (>80%) have been demonstrated to be essential, 4 (*Atp8b3*, *Btrc*, *Cfap206*, and *Ptdss2*) exhibit impaired male fertility¹⁰⁷⁻¹¹⁰, and 10 (*Dnaaf1*, *Eif3e*, *Ip011*, *Ndufa7*, *Rpa1*, *Slc2a3*, *Smarcb1*, *Smc5*, *Ssrp1*, and *Tmem258*) display embryonic lethality⁹⁸,¹¹¹⁻¹¹⁵. *Rpa1* (replication protein A1) and *Slc2a3* (solute carrier family 2, member 3, also known as *Glut3*) are haploinsufficient¹¹⁶,¹¹⁷, indicating the essence of sufficient dosage of gene products for their normal functions. As expected with increased miRNA-guided AGO2-mediated cleavage, we detected significantly decreased RPF abundance from these target RNAs in comparison to the control mRNAs with a median decrease of 29% (Supplementary Fig. 5c, *p* = 1.8 × 10⁻⁵). Given the sensitivity to gene dosage and the essential functions of these target genes, increased miRNA-guided AGO2-mediated cleavage of their transcripts and decreased protein synthesis may contribute to the infertility of *Mov10l1* CKO mutants. Considering that AGO2 is just one of the uppl mRNA protein products, our data support the biological significance of 3'UTR piRNA biogenesis in fine-tuning protein abundance during normal development.

**The biogenesis of 3'UTR piRNAs is evolutionarily conserved.**
To test whether the biogenic mechanisms for 3'UTR piRNAs are also seen for 3'UTR piRNAs found in other amniotes, we identified and annotated 3'UTR piRNAs in roosters (*Gallus gallus*). We used RNA-seq from roosters to assemble the testis-specific transcriptome and then aligned piRNAs to annotated mRNAs. To identify precursor transcripts, we required a piRNA abundance of >100 ppm and ≥90% of piRNAs mapping to 3'UTRs (a median percentage of 96.3% of mouse 3'UTR piRNAs derived from 3' UTRs). To ensure that mRNAs are translated in rooster testes, we also required an RPF abundance ≥1 ppm from their ORFs. Using these criteria, we detected, in total, 23 transcripts that both produce piRNAs (Fig. 6a) and code for proteins (as shown by 3-nucleotide periodicity, Fig. 6b), thus representing uppl mRNAs. The TSSs of these transcripts (but not their 3'UTRs) have a nearby H3K4me3 ChIP-seq peak, a signature of RNA Pol II TSSs¹¹⁸ (Fig. 6a), arguing against the existence of 3'UTR-specific isoforms, and 21 out of 23 H3K4me3 ChIP-seq peaks completely overlap with A-MYB ChIP-seq peaks (Supplementary Fig. 6a), indicating that 3'UTR piRNAs also exist in chickens and are derived from full-length mRNA precursors.

The 23 chicken uppl mRNAs include genes such as *pi-CRTC1* (CREB regulated transcription coactivator 1), *pi-DOT1L* (DOR1-like, histone H3 methyltransferase), and *pi-USP53* (ubiquitin specific peptidase 53). All 23 genes have mouse homologs, and

mouse mutants for 17 of the 23 have been reported⁹⁷. Out of the 17 genes, 10 (59%) are essential in mice, with 3 impairing fertility and the rest causing embryonic lethality. Among the 30 mouse uppl genes, 23 have a homolog in chickens, none of which robustly produce piRNAs in rooster testes (Supplementary Fig. 6b). The 23 chicken uppl genes have no overlap with the mouse uppl genes, and their mouse homologs produce few piRNAs in adult mouse testes (Supplementary Fig. 6c). We used these 23 chicken transcripts as non-piRNA-producing control mRNAs in the following analyses in chickens. While none of the mouse uppl mRNAs come from sex chromosomes (X or Y chromosomes), explained by meiotic sex chromosome inactivation, 8 out of 23 (35%) of the chicken uppl genes mapped to the Z chromosomes (the bird sex chromosomes that are not inactivated because rooster is the homogametic sex), which are significantly enriched compared to all the mRNA-encoding genes we assembled for rooster testes (5.7%, χ² test, *p* = 4.0 × 10⁻⁶). Thus, mouse and chicken 3'UTR piRNAs are derived from diverse, non-overlapping genes, the majority of which are essential for viability and fertility.

We found that RPFs also extended into the 3'UTRs of chicken uppl mRNAs (Fig. 6a). After controlling for the mRNA levels measured by RNA-seq, we found a significant partial correlation between piRNA abundance and RPF abundance from each chicken uppl mRNA 3'UTR (*r* = 0.96, *p* = 1.5 × 10⁻¹², Supplementary Fig. 6d, left), indicating that the ribosome-bound uppl 3' UTRs are processed into piRNAs in roosters. The 5'-ends of RPFs from chicken uppl 3'UTRs significantly overlapped with the 5'-ends of piRNAs (Fig. 6c, right, and Supplementary Fig. 6e). The 5'-ends of 3'UTR RPFs also displayed a 1U bias (Fig. 6d). Authentic piRNAs with 1U bias are produced from uppl ORFs (Fig. 6d). Unlike 3'UTR RPFs, the abundance of uppl ORF RPFs does not correlate with piRNA abundance (*r* = −0.09, *p* = 0.69, Supplementary Fig. 6d, right), and the ORF RPFs did not display a signature of in vivo cleavage (Fig. 6d) nor correspond to future piRNA sites (Fig. 6c, left). Thus, chicken uppl ribosomes also guide endonucleolytic cleavages that generate piRNA 5'-ends in 3' UTRs, but do not do so in ORFs. In sum, although the mRNAs that produce piRNAs do not overlap between mice and chickens, the existence of ribosome-guided piRNA biogenesis from mRNA 3'UTRs in both mice and chickens suggests that an evolutionary conserved biogenic mechanism predates the divergence of mice and chickens approximately 330 million years ago¹¹⁹.

**Transposon fragments are embedded in precursor mRNAs.** To determine the common features of 3'UTR piRNA precursors in chicken and mice, we revisited the debate over whether 3'UTR piRNA precursors harbor TE sequences⁵⁶,¹²⁰,¹²¹. We tested whether mouse uppl mRNAs are enriched in TE sequences or are produced from regions with TEs nearby when compared to the control mRNAs. We determined the distance from the start site of each transcript to the nearest TE in the genome and found that uppl mRNAs are not significantly closer to TEs than the control mRNAs (Supplementary Fig. 7a, *p* ≥ 0.20). We computed the percentages of exonic and intronic nucleotides that are annotated as part of a TE for each locus and found that a significantly higher fraction of uppl mRNA exons harbor short interspersed nuclear elements (SINEs; *p* ≤ 2.7 × 10⁻⁸) in comparison to the control mRNA genes that contained no exonic nucleotides corresponding to TEs (Fig. 7a). The intronic regions of uppl mRNA genes and the control mRNA genes contained a similar fraction of TEs (*p* ≥ 0.02, Supplementary Fig. 7b). These SINE-containing exonic regions almost exclusively correspond to the 3'UTRs (Supplementary Fig. 7c). Thus, TEs are enriched in spliced uppl mRNAs in mice.
<figure><img src="image_20.png" /><figcaption>**Fig. 6** Ribosome-guided 3′UTR piRNA formation is evolutionarily conserved. **a** Normalized reads of RPFs and piRNAs mapping to four representative chicken uppl genes in adult rooster testes. In roosters, both strands of uppl mRNAs code for piRNAs; blue represents Watson strand mapping reads; red represents Crick strand mapping reads. **b** Discrete Fourier transformation of the distance spectrum of 5'-ends of RPFs across ORFs (gold) and 3'UTRs (green) of chicken uppl mRNAs in adult rooster testes. **c** Distance spectrum of 5'-ends of RPFs from uppl ORFs (left) and from 3'UTRs (right) that overlap piRNAs in adult rooster testes (sample size n = 3 independent biological replicates). Data are mean ± standard deviation. **d** Sequence logos depicting nucleotide bias at 5'-ends and 1 nt upstream of 5'-ends of the following species from adult rooster testes. Top to bottom: RPF species and piRNAs, which map to control mRNAs, uppl ORFs, and 3'UTRs, respectively.</figcaption></figure>

In adult mouse testes, we detected TE piRNAs produced from uppl mRNAs (Fig. 7b) that uniquely map to SINEs embedded in uppl mRNA 3′UTRs. Both sense and antisense piRNAs (according to their orientation with SINEs) exhibited a 1U bias not a 10A bias (Fig. 7b and Supplementary Fig. 7d), displaying a signature of primary piRNAs rather than secondary piRNAs. To test whether these piRNAs could trigger ping-pong amplification cycles in trans, we searched for secondary piRNAs from total SINE piRNAs in adult testis (piRNAs mapping to consensus SINE sequences with up to three mismatches allowed throughout the piRNA sequences). We found that antisense SINE-piRNAs, but not sense SINE-piRNAs, produced from uppl mRNAs display a significant ping-pong signature with total SINE-piRNAs (Fig. 7c

and Supplementary Fig. 7e), indicating that these 3′UTR piRNAs can trigger post-transcriptional cleavage of SINE transcripts. To test the biological significance of why only SINEs but not other TE superfamilies are embedded in mouse uppl mRNAs, which cannot be explained by their percentage in the mouse genome (8.22% is SINE, 19.20% is long interspersed nuclear element (LINE), 9.87% is long terminal repeat (LTR), and 0.88% is DNA transposons)¹²², we compared their expression around the developmental stage when uppl piRNAs are produced. We found that SINEs are more highly expressed compared to DNA, LINE, and LTR transposons (Fig. 7d). Therefore, our data indicate that a subset of piRNAs produced from mouse uppl mRNAs post-transcriptionally silence TEs.
<figure><img src="image_21.png" /><figcaption>**Fig. 7** Spliced piRNA precursor mRNAs contain transposon fragments. **a** Boxplots showing the fraction of transcript exon sequence correspondence to sense (blue) and antisense (red) transposon sequences in mouse genome. Control mRNA *n* = 43, uppl mRNA *n* = 30. **b** Sequence logo showing the nucleotide composition of antisense SINE-piRNA species that uniquely map to mouse uppl mRNAs. **c** The 5'-5' overlap between piRNAs from opposite strands of SINE consensus sequences was analyzed to determine whether antisense SINE piRNAs from mouse uppl mRNAs display Ping-Pong amplification in trans. The number of pairs of piRNA reads at each position is reported (sample size *n* = 3 independent biological replicates). Data are mean ± standard deviation. The Z score indicates that a significant ten-nucleotide overlap (Ping-Pong) was detected. Z score = 1.96 corresponds to *p* value = 0.025. **d** The RNA abundance of each TE superfamily in testes at the five developmental time points. Dpp: days postpartum, TPM: transcript per million. **e** Boxplots showing the fraction of transcript exon sequence correspondence to sense (blue) and antisense (red) transposon sequences in chicken genome. Control mRNA *n* = 23, uppl mRNA *n* = 23. **f** Sequence logo showing the nucleotide composition of antisense LINE-piRNA species that uniquely map to chicken uppl mRNAs. **g** The 5'-5' overlap between piRNAs from opposite strands of LINE consensus sequences was analyzed to determine whether antisense LINE-piRNAs from chicken uppl mRNAs display Ping-Pong amplification in trans (sample size *n* = 3 independent biological replicates). Data are mean ± standard deviation. Z score = 1.96 corresponds to *p* value = 0.025. Box plots in **a**, **e** show the 25th and 75th percentiles, whiskers represent the 5th and 95th percentiles, and midlines show median values.</figcaption></figure>

To test whether the function of producing antisense TE piRNAs is biologically significant during evolution, we performed a similar analysis in chickens and found that the processing of chicken 3'UTR piRNA precursors also generates antisense TE piRNAs. Chicken piRNA precursor mRNAs harbor TEs in their

exonic regions that correspond to 3'UTRs (Fig. 7e and Supplementary Fig. 7f). While SINEs are largely absent from the chicken genome¹²³, we found significantly more LINEs embedded in chicken uppl mRNA exons (*p* = 2.5 × 10⁻⁴), whereas control mRNA genes contained no exonic nucleotides
corresponding to TEs (Fig. 7e). Unlike mice, where both sense and antisense TEs are embedded in uppl mRNAs, an antisense bias of TEs was detected in chickens (sense median = 0.0%, antisense median = 2.1%). LINE piRNAs produced in chicken adult testes manifest a 1U bias but not a 10A bias (Fig. 7f). These piRNAs, uniquely mapping to uppl mRNAs, can also target TEs in trans, post-transcriptionally generating secondary piRNAs (Fig. 7g). Although neither uppl genes nor TE families (Fig. 7a, e) are conserved between mice and chickens, which may accommodate rapidly changing populations of TEs, our data show that the use of mRNAs embedded with TEs to produce antisense TE piRNAs that cleave TEs post-transcriptionally is a common strategy in amniotes.

## Discussion

Here we systematically characterized 3'UTR piRNAs (which should be called genic piRNAs, as ORF regions also produce piRNAs). By defining the transcription factors associated with piRNA biogenesis and characterizing mutants with transcriptional and post-transcriptional processing defects, we demonstrate that full-length mRNAs are the precursors of 3'UTR piRNAs. These mRNAs undergo pioneer rounds of translation that are followed by the production of piRNAs. This coupling is mediated by post-termination 80S ribosomes on 3'UTRs that guide endonucleolytic cleavages and generate the 5'-ends of piRNAs. Together with our previous studies on piRNA biogenesis from lncRNA precursors⁴⁵, we find that ribosomes guide piRNA biogenesis downstream of ORFs regardless of ORF length. Similar to other co-translational mRNA quality-control pathways, piRNA processing from mRNAs fine-tunes the protein products from these mRNAs. This co-translational processing of piRNAs is found in both mice and chickens. Therefore, we reveal a general and conserved mechanism by which post-termination ribosomes guide piRNA 5'-end formation from non-protein-coding regions of RNAs in amniotes.

TE silencing in mice is thought to be carried out principally by prenatal piRNAs¹²⁴,¹²⁵. However, TE surveillance is still required throughout spermatogenesis. The murine PIWI protein MILI is not detected after prenatal piRNA expression ends and before pachytene piRNA expression begins¹²⁶, resulting in a loss of production of prenatal piRNAs and potentially leaving TEs unrestrained after birth and before pachynema. Indeed, we observed a very significant increase in the expression load of SINEs from 10.5 dpp to 20.5 dpp. The expression of 3'UTR piRNAs throughout spermatogenesis (Supplementary Fig. 1a) may be critical for germline genome protection against activated TEs when TE piRNAs produced from other sources are low. Since TE piRNAs can also cleave other mRNAs that harbor TEs, as previously reported¹²⁷, the TE piRNAs produced from uppl mRNAs may lead to a cascade of events to regulate mRNA stability. Given the rapid changes in TE families in each species over evolutionary time, the shared mechanisms for producing TE piRNAs from distinct gene subsets in divergent lineages suggest the biological significance of this strategy.

Furthermore, we discovered that fine-tuning of essential protein levels is a function of 3'UTR piRNA processing. One target of co-translational piRNA processing in mice spermatocytes is Ago2 mRNAs. Consistent with the idea that piRNA processing attenuates AGO2 function, defective piRNA processing in *Mov10l1* mutants leads to increased AGO2 protein levels and decreased protein synthesis of AGO2-miRNA targets. Moreover, most AGO2-miRNA targets in spermatocytes harbor essential functions, allowing the initial impact on AGO2 protein levels to be further amplified. Although it is technically challenging, if even possible, to single out the impact of each uppl protein change on

spermatogenesis, the combination of the dysregulation of these genes is likely to be a non-negligible contributor to the *Mov10l1* mutant phenotype. While the set of protein-coding genes that produce piRNAs is not evolutionarily conserved between mice and chickens, the functions of uppl mRNA genes are in general critical for spermatogenesis and viability in mice and chickens. The high enrichment of chicken uppl mRNA genes from the Z chromosome might reflect the requirements for dosage compensation from the two Z chromosomes in males, given that in chickens females are ZW with only one Z chromosome and the piRNA pathway is not active in the chicken ovary¹²⁸. Taken together, although fine-tuning protein is unlikely to represent the primary selective force that marks a subset of mRNAs as piRNA precursors, we expect the better fitness gained by the rapid turnover of these essential mRNAs to further reinforce their piRNA precursor identity.

The ribosome mitigation downstream of ORFs regardless of ORF length suggests compromised ribosome recycling. During conventional termination, ABCE1 binds to post-termination ribosomes and splits the large and small subunits¹²⁹,¹³⁰. After 60S dissociation, a complex including eIF2D/Tma64 (the homolog of eIF2D in yeast), MCST1/Tma20, and DENR/Tma22 is required for 40S ribosome recycling in yeast. In general, translation termination is a highly controlled process, with the binding of 80S ribosomes to 3'UTRs only being demonstrated in either terminally differentiated erythrocytes or in the context of mutated translational machinery⁷⁹⁻⁸⁴,¹³¹⁻¹³⁴. In our study, the regulation of ribosome recycling occurs during meiosis prophase in cells that are not terminally differentiated, as spermatocytes will undergo two more rounds of cell division, weeks of further differentiation, and continued protein synthesis. This ribosome recycling regulation is specific to uppl mRNAs and may be due to the localization of uppl mRNAs in proximity to the mitochondria where piRNAs are processed. The process is enabled, at least partially, by the reprogramming of RNA metabolism at pachynema, including the inhibition of ribosome recycling and the NMD pathway, allowing for robust piRNA production in a short but critical time window during spermatogenesis. Our study indicates that reprogramming of ribosome recycling can occur locally and stage-specifically to enable biologically significant processes.

Our systematic characterization of 3'UTR piRNAs provides insights into the fundamental question: why are some transcripts processed to piRNAs while others are not? piRNA precursors could be either marked epigenetically during transcription or are recognized post-transcriptionally in the cytosol. Unlike *Drosophila* germline piRNA loci marked with a chromatin-bound protein Rhino³⁰, currently no epigenetic factors have been identified to specially bind to piRNA loci in mammals. Although splicing signature is unlikely to be the hallmark for piRNA precursors in mammals as proposed for *Drosophila* piRNA biogenesis⁵⁹, both lncRNA piRNA precursors and mRNA piRNA precursors exhibit longer first exons, suggesting that unique exon-intron structure could be one unique feature. Furthermore, given the translation of short upstream ORFs is insufficient to mark a transcript for piRNA biogenesis and the uppl mRNAs do not exhibit faulty translation on their main ORFs, the translation intermediates with post-termination ribosomes on a long 3'UTR is a prime candidate for further testing. Last but not least, if the TE-rich prenatal piRNAs target and initiate the processing of the 3'UTR piRNA precursors, the embedding of TE elements could also serve as a determining feature of a transcript for piRNA processing. Thus, the study of 3'UTR piRNAs allows a more comprehensive investigation into the unique features defining piRNA precursors.

The linked but distinct biogenic requirements before and after the stop codon of uppl main ORFs suggest that post-
transcriptional processing of piRNAs can be further divided into two phases: substrate recognition and efficient processing. Since uppl ORFs are still processed into piRNAs distinct from other cellular mRNAs, this suggests that uppl mRNAs are recognized as substrates for piRNA processing as entire transcripts before endonucleolytic cleavages occur. Otherwise, if the substrate recognition couples with piRNA processing, as soon as the initial cleavage starts on 3′UTRs, the ORF portion of the transcripts would display no unique feature compared to other translating mRNAs. Our finding that the ORF portion of transcripts still enter piRNA processing but in an inefficient manner suggests that either post-termination ribosomes streamline piRNA production and/or the translating ribosomes inhibit piRNA processing machinery from accessing the RNAs. The requirement of TDRD5 for 3′UTR piRNA biogenesis suggests that TDRD5 may function by coordinating post-termination ribosomes with piRNA processing machinery. TDRD5 detected in amniotes (which harbor a homolog with >50% protein sequence identity with the mouse TDRD5) but not in fish (<10% of potential homologs)¹³⁵ may have co-evolved with appearance of ribosome-guided piRNA biogenesis. Furthermore, although 3′UTR piRNAs have been characterized in fruit flies and the cis RNA elements sufficient to promote piRNA biogenesis in somatic cells have been identified¹³⁶,¹³⁷, their biogenesis remains when upstream translation is inhibited¹³⁷, and the region downstream of the cis RNA element produce fewer piRNAs than the element region itself¹³⁶. Thus, consistent with the recent appearance of TDRD5, ribosome-guided piRNA biogenesis is unlikely to be conserved in invertebrates and further studies are required to trace its evolutionary origins.

In summary, we reveal a conserved and general piRNA biogenesis mechanism that recognizes translating RNAs regardless of whether they harbor long ORFs or not. The assembly of 80S ribosomes on non-coding regions of RNA is not restricted by the length of the upstream ORFs and is temporally staggered with translation-dependent RNA quality-control pathways, suggesting compromised ribosome recycling. The coupling of piRNA biogenesis with translation fine-tunes the abundance of proteins that are critical for spermatogenesis in both mice and chickens.

## Methods

**Animals.** Mice were maintained and used according to guidelines for animal care of the NIH and the University Committee on Animal Resources at the University of Rochester. Mice of the following strains C57BL/6J (Jackson Labs, Bar Harbor, ME, USA; stock number 664); Rpl22tm1.1.Psam on a C57BL/6J background (Jackson Labs; stock number 011029)¹³⁸; Mov10l1tm1.1/w on a mixed 129×1/SvJ × C57BL/6J background¹¹; and Tg(Neurog3-cre)C1Able/J on a B6.FVB(Cg) background (Jackson Labs; stock number 006333)¹³⁹ were genotyped as described. Comparisons of compound mutants and controls involving Mov10l1 CKO mutation were performed using siblings from individual litters. White Leghorn testes of the Cornell Special C strain from 1-year-old roosters were used according to guidelines for animal care of the NIH and the University Committee on Animal Resources at the University of Rochester.

**Small RNA sequencing library construction.** Small RNA libraries were constructed and sequenced, as previously described⁴⁵, using oxidation to enrich for piRNAs by virtue of their 2'-O-methyl-modified 3' termini. The oxidation procedure selects against in vitro-digested products with a 3' phosphate. A 25-mer RNA with 2'-O-methyl-modified 3' termini (Supplementary Table 1, Spike-in RNA) was used as a spike-in control.

**Western blotting.** For protein detection, testis lysates were resolved by electrophoresis on 10% sodium dodecyl sulfate (SDS)-polyacrylamide gels. The proteins were transferred to a 0.45 µm polyvinylidene difluoride membrane (EMD Millipore, Billerica, MA, USA), and the blot was probed with anti-AGO2 mouse monoclonal antibody (Wako Pure Chemical Corporation, 018-22021, 1:500), anti-TUBULIN rabbit antibody (Bimake, Houston, TX, USA, A5105, 1:1000), and then detected with sheep anti-mouse immunoglobulin G-horseradish peroxidase (IgG-HRP; GE Healthcare, Little Chalfont, UK; NA931V, 1:5000), and donkey anti-

rabbit IgG-HRP (GE Healthcare, NA934V, 1:5000). Western blotting images were taken using Azure c300 imaging system, with the cSeries capture software (v1.6).

**Polysome gradient.** Fresh testes were lysed in 1 ml lysis buffer (10 mM Tris-HCl, pH 7.5, 5 mM MgCl₂, 100 mM KCl, 1% Triton X-100, 2 mM dithiothreitol (DTT), 100 µg/ml cycloheximide, and 1× protease-inhibitor cocktail) as previously described⁴⁵. Five A₂₆₀ absorbance units were loaded on a 10–50% (w/v) linear sucrose gradient prepared in buffer (20 mM HEPES-KOH, pH 7.4, 5 mM MgCl₂, 100 mM KCl, 2 mM DTT, 100 µg/ml cycloheximide) and centrifuged in a SW-40ti rotor at 154,348 × g for 2 h 40 min at 4 °C. Samples were collected from the top of the gradient using a gradient Fractionation system (Brandel, Boca Raton, FL, USA; BR-188) while monitoring absorbance at 254 nm. Synthetic spike-in RNAs were added to each collected fraction before RNA purification. RNA was purified by Trizol (Ambion, Waltham, MA, USA) using the Direct-zol™ RNA MiniPrep plus Kit (Zymo Research, Irvine, CA, USA).

**Ribosome profiling.** Ribo-seq was performed as previously described⁴⁵. Cleared testis lysates were incubated with 60 units of RNase T1 (Fermentas, Waltham, MA, USA) and 100 ng of RNase A (Ambion) per A₂₆₀ unit for 30 min at room temperature. Samples were loaded on sucrose gradients, and after centrifugation, the fractions corresponding to 80S monosomes were recovered for library construction.

The 5'P and 5'OH Ribo-seq library was prepared as follows (Supplementary Fig. 4a). (i) Size-selected and rRNA-depleted RPFs were ligated to a 5' adapter. The ligation products were resolved on a 15% denaturing gel, and ligated products and unligated RNAs were purified separately. (ii) The 3'-ends of recovered RNAs were dephosphorylated and ligated to 3' adapters (Supplementary Table 1, 3' adapter). (iii) The unligated RNAs proceed to the conventional Ribo-seq library construction. (iv) The ligated RNAs directly proceed to the reverse transcription steps for library construction.

**RNA sequencing.** Strand-specific RNA-seq libraries were constructed following the TruSeq RNA sample preparation protocol, as previously described⁴⁵. rRNAs were depleted from total RNAs with complementary DNA oligomers (IDT) and RNase H (Invitrogen, Waltham, MA, USA)¹⁴⁰,¹⁴¹. RNA-seq data were generated on HiSeq 2000 instrument (Illumina, San Diego, CA, USA).

**Degradome-seq library construction.** Degradome-seq library construction was performed as described previously⁴⁵. The RNAs were first oxidized at room temperature for 30 min with sodium periodate (Sigma, St. Louis, MO, USA) to block the 3'-ends from ligation and were then size-selected to isolate RNA ≥200 nts (DNA Clean & Concentrator™-5, Zymo Research). 5' Adapters were attached using T4 RNA ligase (Ambion) at 20 °C for 3 h. The ligated products were subjected to rRNA depletion with complementary DNA oligomers (IDT) and RNase H (Invitrogen)¹⁴⁰,¹⁴¹. The rRNA-depleted ligation products were reverse transcribed using a degenerate primer (Supplementary Table 1, Degenerate primer). cDNA was amplified by PCR using KAPA HIFI Hotstart polymerase (Kapa Biosystems, Wilmington, MA, USA), and 250–350 nts double-stranded DNA was isolated on 8% native PAGE gels.

**Histology and immunostaining.** For histologic analysis, testes were fixed in Bouin's solution overnight, embedded in paraffin, and sectioned at 4 µm. Following standard protocols, sections were deparaffinized, rehydrated, and then stained with hematoxylin and eosin.

**PELOTA/DOM34 immunostaining** and mouse monoclonal anti-HA antibody (ascites fluid, 1:2000 dilution; Covance, Princeton, NJ, USA; MMS-101P) were performed on squashed spermatocytes and spermatids as previously described¹⁴². Seminiferous tubules were fixed in 2% paraformaldehyde containing 0.1% Triton X-100 for 10 min at room temperature, placed on a slide coated with 1 mg/ml Poly-L-lysine (Sigma) with a small drop of fixative, gently minced with tweezers, and squashed. The coverslip was removed after freezing in liquid nitrogen. The slides were later rinsed three times for 5 min in phosphate-buffered saline (PBS) and incubated for 12 h at 4 °C with rabbit anti-PELOTA antibody (1:50 dilution; Thermo, PA5-31697). Secondary antibodies conjugated with Alexa Fluor 488 (Molecular Probes, Eugene, OR, USA) were used at a dilution of 1:500. Histology and immunostaining images were taken using Leica DM4000 B LED microscope system with the Leica software: Leica Application Suite X v1.1.0.12420.

**IP from mouse testes.** For CBP80 and eIF4E IPs, anti-CBP80 rabbit antibody (Bethyl Laboratories, A301-794A) or anti-eIF4E rabbit antibody (Bethyl Laboratories, A301-153A) from mouse testis lysates were performed as described¹⁴³, except Protein A magnetic Beads (161-4013, Bio-Rad, Hercules, CA, USA) were used.

**Intra-testicular injection.** The mice were anesthetized with ketamine/xylazine mixture (ketamine 100 mg/kg; xylazine 25 mg/kg) via intraperitoneal injection. After complete anesthesia, testes were exteriorized with a longitudinal incision around 1 cm at the center of abdomen. The tunica albuginea was penetrated using
a sharp 26 G needle (BD, Franklin Lakes, NJ, USA, 30511) 1 mm from the vascular pedicle, and the needle was withdrawn to generate a path for introducing a blunt end Hamilton needle (Hamilton, Reno, NV, USA, 7786-02). PBS containing 0.02% Fast Green FCF (Thermo Fisher Scientific, Waltham, MA, USA) with harringtonine (LKT labs, 0.5 µg/µl in a total volume of 10 µl) or with Okad acid (LX Laboratories, Woburn, MA, USA, O-2220, 16 nM per testis in a total volume of 10 µl) was slowly injected using a Hamilton microsyringe (1705RN) into one testis, and a vehicle control without the drug was injected into the other testis of the same animal. The needle was held in place for 30 s before removal to prevent backflow of the solution. Successful completion of injection was indicated by testis filled with green solution. The testes were returned to the abdominal cavity after injection. The incisions were sutured. At the end point, the mice were euthanized by cervical dislocation, and the testes were collected.

**General bioinformatics analyses.** Analyses were performed using piPipes v1.4¹⁴⁴. All data from the small RNA sequencing, Ribo-seq, 40S footprinting, 80S foot-printing, RNA-seq, degradome-seq, CLIP sequencing, and ChIP-seq were analyzed using the latest mouse genome release mm10 (GCA_000001635.7) and chicken genome release galGal6 (GCA_000002315.5). Generally, one mismatch is allowed for genome mapping. For mouse transcriptome annotation, 30 upl mRNAs defined in our previous studies with mm9⁸ were converted to mm10 coordinates with using *liftOver*¹⁴⁵ with minor manual correction (Supplementary Data 1). We selected 43 control mRNAs (Supplementary Data 1) from our recently reconstructed mouse testis transcriptome⁶⁰ by their similar expression dynamics as upl mRNAs from 10.5 dpp to 20.5 dpp and by their lack of any piRNA production. For genes with alternative transcripts, the most abundant transcript involving that gene was selected. We reassembled mRNAs using RNA-seq data from rooster testes and have defined 3'UTR piRNA precursor mRNAs as described in the text. Statistics pertaining to the high-throughput sequencing libraries constructed for this study are provided in Supplementary Data 1.

For small RNA-seq, libraries were normalized to the sum of total miRNA reads; spike-in RNA was used to normalize the libraries from each fraction of polysome profiling. Uniquely mapping reads >23 nts were selected for further piRNA analysis. We analyzed previously published small RNA libraries from adult wild-type mouse testis (GSM1096604)⁸, mouse *Tdrd5* CKO, control testis (SRP093845)⁷⁸, *A-Myb* mutant, control testis at 14.5 and 17.5 dpp (GSM1096588, GSM1096589, GSM1096590, GSM1096591, GSM1096605, GSM1096606, GSM1096607, GSM1096608), RNase treated (GSM4160780) and untreated adult testis (GSM4160781), *Mov10l1* CKO mutants (GSM4160774, GSM4160775, GSM4160776, GSM4160777, GSM4160778, and GSM4160779), and littermate controls (GSM4160768, GSM4160769, GSM4160770, GSM4160771, GSM4160772, and GSM4160773)⁴⁵. We report piRNA abundance either as parts per million reads mapped to the genome (ppm), or as reads per kilobase pair per million reads mapped to the genome (rpkm) using a pseudo count of 0.001. To ensure precision of the mapping, a piRNA is counted only when the 5'-end of the piRNA maps to the ORF or 3'UTR of a transcript.

For RNA-seq reads, the expression per transcript was normalized to the top quartile of expressed transcripts per library calculated by Cufflinks v2.2.1⁴⁶, and the transcripts per million value was quantified using the Salmon v0.8.2 algorithm¹⁴⁷. We analyzed previously published RNA-seq libraries from mouse *A-Myb* mutant, control testis at 14.5 and 17.5 dpp (GSM1088426, GSM1088427, GSM1088428, GSM1088429), wild-type mouse testis at 10.5 dpp, 12.5 dpp, 14.5 dpp, 17.5 dpp, 20.5 dpp, and adult stage (GSM1088421, GSM1088422, GSM1088423, GSM1088424, GSM1088425, and GSM1088420)⁸, *Mov10l1* CKO mutants (GSM4160761, GSM4160762, GSM4160753), and littermate controls (GSM4160758, GSM4160759, and GSM4160760)⁴⁵.

Ribo-seq analysis followed the modified small RNA pipeline including the junction mapping reads as previously described⁴⁵. Uniquely mapping reads between 26 and 32 nts were selected for further analysis except for the analysis on the 40S footprints where 18–80 nts were selected. RPFs and 80S footprints from different developmental stages were normalized to the sum of reads mapping to mRNA protein-coding regions, assuming that mRNA translation was largely unchanged during spermatogenesis. Libraries from harringtonine treatment were further normalized to the sum of reads mapping to mitochondrial coding sequences (CDSs) as previously described¹⁴⁸. 40S footprints were normalized to the sum of reads mapping to mRNA 5'UTRs as it has been shown that the 40S binds 5' UTRs in a cap-tethered fashion, thus 40S ribosomes do not accumulate upon harringtonine treatment⁶⁶. We analyzed published Ribo-seq libraries from *Mov10l1* CKO mutants (GSM4160728, GSM4160729, and GSM4160730) and litter mate controls (GSM4160725, GSM4160726, and GSM4160727), anti-HA IP RPF (GSM4160731, GSM4160732, and GSM4160733) and input (GSM4160734, GSM4160735, and GSM4160736) from *RiboTag* mice, wild-type adult testis injected with (GSM4160744 and GSM4160745) and without harringtonine (GSM4160746), in vitro sheared RNA fragments (GSM4160749), monosome (GSM4160752), polysome (GSM4160753), adult testis (GSM4160756 and GSM4160757)⁴⁵, and adult rooster wild-type testis (GSM2454692)¹²⁸.

ChIP-seq reads were analyzed as previously described⁸. Multiple mapping reads were apportioned randomly to each location (-k 1 switch) and one mismatch was allowed (-v 1). ChIP peaks were identified using MACS2 v2.1.1.2063091⁴⁹ using default arguments and input as control. We analyzed published A-MYB ChIP-seq

(GSM1087281) and input DNA libraries (GSM1087284) from mouse wild-type adult testis and A-MYB ChIP-seq (GSM1087285), H3K4me3 ChIP-seq (GSM1087286), and input DNA libraries (GSM1087287) from rooster wild-type adult testes⁸.

Degradome reads and CLIP sequencing reads were aligned to the genome using TopHat 2.0.12¹⁵⁰. Reads were mapped uniquely using the “-g 1” flag. Uniquely mapping reads were selected for further analysis. Libraries were normalized to the sum of reads mapping to mRNA protein-coding regions, assuming that mRNA cleavage was largely unchanged during spermatogenesis. We analyzed a published anti-HA IP degradome library from adult wild-type testis (GSM4160721)¹²⁸ and MOV10L1 CLIP library mouse wild-type testis (PRJNA230507)¹³.

Statistical analyses were performed in R v3.5.0¹⁵¹. The significance of the differences was calculated by Wilcoxon rank-sum test except as indicated in the text. The significance of correlation was performed using partial correlation analysis in addition to simple correlations⁷².

**Chicken transcriptome assembly and annotation.** We used StringTie v1.3.3b¹⁵² with parameter “-m 50 -p 8 -G RefSeq -rf -l prefix” to assemble transcripts from RNA-seq reads from chicken testes. After assembling each sample, StringTie’s merge function was used to merge each assembly. We required the transcript abundance to be at least 8 fpkm and filtered 406,332 assembled transcripts to 39,940 transcripts. We further added transcripts that are not covered in our own assembly but reported in the latest RefSeq (downloaded on Feb 19, 2020), resulting in 81,079 total transcripts. Finally, using TransDecoder v5.5.0¹⁵³ with the BlastP (v2.10.0 +) and Hmmer (v3.3) search, we identified 44,856 mRNA transcripts and considered the remaining 36,223 transcripts as lncRNAs.

**5'-end overlap analysis and phasing analysis.** Reads were mapped to the transcriptome, and their relative positions on transcripts were reported. For alternative transcription with overlapping annotations, we chose the longest transcript. We calculated the distance spectrum of 5'-ends of Set A (RPFs or degradome reads) that overlapped with Set B (piRNAs or simulated sequences) as follows: for each read *b* in Set B, we identified all the reads in Set A whose 5'-ends overlapped within the 50-nT region upstream of *b* including the 5'-end of *b* and 50 nts downstream of *b* (200-nT window of *b* reads). We assigned the *b* spectrum as the fractions of 5'-ends of *a* reads distributed across the 100-nT window of *b* reads. The height of the *b* spectrum at each nucleotide position in the 100-nT window of *b* reads was based on the number of *a* reads whose 5'-ends overlapped at this position divided by the total number of *a* reads whose 5'-ends overlapped with the 100-nT window of *b* reads. The sum of all *b* spectra was then divided by the total number of reads in Set B. We defined this average fraction of 5'-ends of *a* reads that overlapped with the 100-nT window of *b* reads as the distance spectrum of 5'-ends of Set A that overlap Set B. The Z score for overlap at the 5'-end position was calculated using the spectral value at positions -50–-1 and 2–50 as background.

**Nucleotide periodicity.** Nucleotide periodicity was computed as previously described⁴⁵. We first aligned the RPFs to each other using 5'-end overlap analysis and reported the distance spectrum. An annotated ORF was not a prerequisite for this analysis as the distance spectrum of RPFs from mRNAs already showed a 3-n periodicity pattern. We then transformed the distance spectrum using the “per iodogram” function from the *GeneCycle* v1.1.4 package¹⁵⁴ with the “clone” method. The relative spectral density was calculated by normalizing to the value at the first position.

**Generating simulated sequences as negative controls.** We generated a random pool of 28-mer sequences using a sliding window of 1 nucleotide from 5' to 3' of the piRNA precursors. We then sampled from this 28-mer pool to match the first nucleotide composition of the real reads. These simulated sequences from piRNA precursors were used as random controls for piRNAs (source code available upon request).

**Mass spectrometry sample preparation and analysis.** After testes lysis, protein concentration was determined by BCA (Thermo Scientific). Samples were then diluted to 1 mg/ml in 5% SDS, 50 mM triethylammonium bicarbonate (TEAB). In all, 25 µg of protein from each sample was reduced with DTT to 2 mM, followed by incubation at 55 °C for 60 min. Iodoacetamide was added to 10 mM and incubated in the dark at room temperature for 30 min to alkylate the proteins. Phosphoric acid was added to 1.2%, followed by six volumes of 90% methanol and 100 mM TEAB. The resulting solution was added to S-Trap micros (Protifi) and centrifuged at 4000 × g for 1 min. The S-Traps containing trapped proteins were washed twice by centrifuging through 90% methanol and 100 mM TEAB. One microgram of trypsin was brought up in 20 µl of 100 mM TEAB and added to the S-Trap, followed by an additional 20 µl of TEAB to ensure the sample did not dry out. The cap to the S-Trap was loosely screwed on but not tightened to ensure the solution was not pushed out of the S-Trap during digestion. Samples were placed in a humidity chamber at 37 °C overnight. The next morning, the S-Trap was centrifuged at 4000 × g for 1 min to collect the digested peptides. Sequential additions of 0.1% trifluoroacetic acid (TFA) in acetonitrile and 0.1% TFA in 50% acetonitrile
were added to the S-trap, centrifuged, and pooled. Samples were frozen and dried down in a Speed Vac (Labconco), then re-suspended in 0.1% TFA prior to analysis.

Peptides were loaded onto a 100 µm x 30 cm C18 nano-column packed with 1.8 µm beads (Sepax), using an Easy nLC-1200 HPLC (Thermo Fisher) connected to a Orbitrap Fusion Lumos mass spectrometer (Thermo Fisher). Solvent A was 0.1% formic acid in water, and solvent B was 0.1% formic acid in 80% acetonitrile. Ions were delivered to the mass spectrometer using a Nanospray Flex source operating at 2 kV. Peptides were eluted off the column using a multi-step gradient, which started at 3% B and held for 2 min, quickly ramped to 10% B over 7 min, increased to 38% B over 152 min, then ramped up to 90% B in 6 min and held there for 4 min to wash the column before returning to starting conditions in 2 min. The column was re-equilibrated for 7 min for a total run time of 180 min. The flow rate was 300 nl/min. The Fusion Lumos was operated in data-dependent mode, performing a full scan followed by as many MS2 scans as possible in 3 s. The full scan was done over a range of 375–1400 m/z, with a resolution of 120,000 at m/z of 200, an AGC target of 4e5, and a maximum injection time of 50 ms. Peptides with a charge state between 2 and 5 were selected for fragmentation. Precursor ions were fragmented by collision-induced dissociation using a collision energy of 30 and an isolation width of 1.1 m/z. MS2 scans were collected in the ion trap with the scan rate set to rapid, a maximum injection time of 35 ms, and an AGC setting of 1e4. Dynamic exclusion was set to 45 s.

Raw data were searched using SEQUEST within the Proteome Discoverer software platform, v2.2 (Thermo Fisher) employing the SwissProt mouse database, along with a custom fasta database that included both test and control proteins. Trypsin was selected as the enzyme allowing up to 2 missed cleavages, with an MS1 mass tolerance of 10 ppm, and an MS2 mass tolerance of 0.6 Da. Carbamidomethyl on cysteine was selected as a fixed modification. Oxidation of methionine was set as a variable modification. A percolator was used as the false discovery rate calculator, filtering out peptides with a q value > 0.01. Label-free quantitation was performed using the Minora Feature Detector node with a minimum trace length of 5. The Precursor Ions Quantifier node was then used to calculate protein abundance ratios using only unique and razor peptides. The summed abundance-based method was employed, which sums the peak areas for all the peptides for a given protein to determine protein ratios.

**Codon Adaptation Index (CAI) analysis.** CAI was calculated for coding region of uppl mRNAs and control mRNAs using two different software packages: DAMBE v7.2.1¹⁵⁵ and CAIcal v1.4¹⁵⁶. Both DAMBE software package by Xia et al. and CAIcal from Puigbo's group use alternative implementation⁷⁰ of Sharp's formula⁶⁹ to calculate CAI. CAI calculates normalized CAI as quotient between the CAI of the query sequence and expected CAI (eCAI) of 1000 randomly generated sequences with G + C and amino acid content similar to that of the query sequence⁷¹. Random sequences were generated using the Markov method, and eCAI was estimated at 99% level of confidence and 99% coverage. CDSs of housekeeping genes were used as the reference set for all the calculations.

**miRNA target search.** miRNA target search was performed by miRanda v3.3¹⁵⁷ using the -strict parameter.

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

## Data availability

The data supporting the findings of this study are available from the corresponding author upon reasonable request. Next-generation sequencing data used in this study have been deposited at the NCBI Gene Expression Omnibus under the accession number GSE155350. Mass spectrometry data have been uploaded to the ProteomeXchange Consortium via the PRIDE database under accession number PXD027489. Source data are provided with this paper.

Received: 6 July 2021; Accepted: 17 September 2021;
Published online: 13 October 2021

## References

1. Kumar, M. & Carmichael, G. G. Antisense RNA: function and fate of duplex RNA in cells of higher eukaryotes. *Microbiol. Mol. Biol. Rev.* **62**, 1415–1434 (1998).

2. Aravin, A. Å. & Hannon, G. J. Small RNA silencing pathways in germ and stem cells. *Cold Spring Harb. Symp. Quant. Biol.* **73**, 283–290 (2008).

3. Farazi, T. A., Juranek, S. A. & Tuschl, T. The growing catalog of small RNAs and their association with distinct Argonaute/Piwi family members. *Development* **135**, 1201–1214 (2008).

4. Kim, V. N., Han, J. & Siomi, M. C. Biogenesis of small RNAs in animals. *Nat. Rev. Mol. Cell Biol.* **10**, 126–139 (2009).

5. Thomson, T. & Lin, H. The biogenesis and function of PIWI proteins and piRNAs: progress and prospect. *Annu. Rev. Cell Dev. Biol.* **25**, 355–376 (2009).

6. Cenik, E. S. & Zamore, P. D. Argonaute proteins. *Curr. Biol.* **21**, R446–R449 (2011).

7. Brennecke, J. et al. Discrete small RNA-generating loci as master regulators of transposon activity in Drosophila. *Cell* **128**, 1089–1103 (2007).

8. Li, X. Z. et al. An ancient transcription factor initiates the burst of piRNA production during early meiosis in mouse testes. *Mol. Cell* **50**, 67–81 (2013).

9. Frost, R. J. et al. MOV10L1 is necessary for protection of spermatocytes against retrotransposons by Piwi-interacting RNAs. *Proc. Natl Acad. Sci. USA* **107**, 11847–11852 (2010).

10. Zheng, K. et al. Mouse MOV10L1 associates with Piwi proteins and is an essential component of the Piwi-interacting RNA (piRNA) pathway. *Proc. Natl Acad. Sci. USA* **107**, 11841–11846 (2010).

11. Zheng, K. & Wang, P. J. Blockade of pachytene piRNA biogenesis reveals a novel requirement for maintaining post-meiotic germline genome integrity. *PLoS Genet.* **8**, e1003038 (2012).

12. Vourekas, A. et al. Mili and Miwi target RNA repertoire reveals piRNA biogenesis and function of Miwi in spermiogenesis. *Nat. Struct. Mol. Biol.* **19**, 773–781 (2012).

13. Vourekas, A. et al. The RNA helicase MOV10L1 binds piRNA precursors to initiate piRNA processing. *Genes Dev.* **29**, 617–629 (2015).

14. Ipsaro, J. J., Haase, A. D., Knott, S. R., Joshua-Tor, L. & Hannon, G. J. The structural biochemistry of Zucchini implicates it as a nuclease in piRNA biogenesis. *Nature* **491**, 279–283 (2012).

15. Nishimasu, H. et al. Structure and function of Zucchini endoribonuclease in piRNA biogenesis. *Nature* **491**, 284–287 (2012).

16. Voigt, F. et al. Crystal structure of the primary piRNA biogenesis factor Zucchini reveals similarity to the bacterial PLD endonuclease Nuc. *RNA* **18**, 2128–2134 (2012).

17. Vagin, V. V. et al. A distinct small RNA pathway silences selfish genetic elements in the germline. *Science* **313**, 320–324 (2006).

18. Kawaoka, S. et al. Zygotic amplification of secondary piRNAs during silkworm embryogenesis. *RNA* **17**, 1401–1407 (2011).

19. Tang, W., Tu, S., Lee, H. C., Weng, Z. & Mello, C. C. The RNase PARN-1 trims piRNA 3' ends to promote transcriptome surveillance in *C. elegans*. *Cell* **164**, 974–984 (2016).

20. Izumi, N. et al. Identification and functional analysis of the pre-piRNA 3' trimmer in silkworms. *Cell* **164**, 962–973 (2016).

21. Lim, S. L. et al. HENMT1 and piRNA stability are required for adult male germ cell transposon repression and to define the spermatogenic program in the mouse. *PLoS Genet.* **11**, e1005620 (2015).

22. Ding, D. et al. PNLDC1 is essential for piRNA 3' end trimming and transposon silencing during spermatogenesis in mice. *Nat. Commun.* **8**, 819 (2017).

23. Nishimura, T. et al. PNLDC1, mouse pre-piRNA Trimmer, is required for meiotic and post-meiotic male germ cell development. *EMBO Rep.* **19**, e44957 (2018).

24. Zhang, Y. et al. An essential role for PNLDC1 in piRNA 3' end trimming and male fertility in mice. *Cell Res.* **27**, 1392–1396 (2017).

25. Gainetdinov, I., Colpan, C., Arif, A., Cecchini, K. & Zamore, P. D. A single mechanism of biogenesis, initiated and directed by PIWI proteins, explains piRNA production in most animals. *Mol. Cell* **71**, 775.e5–790.e5 (2018).

26. Saxe, J. P., Chen, M., Zhao, H. & Lin, H. Tdrkh is essential for spermatogenesis and participates in primary piRNA biogenesis in the germline. *EMBO J.* **32**, 1869–1885 (2013).

27. Horwich, M. D. et al. The Drosophila RNA methyltransferase, DmHen1, modifies germline piRNAs and single-stranded siRNAs in RISC. *Curr. Biol.* **17**, 1265–1272 (2007).

28. Saito, K. et al. Pimet, the Drosophila homolog of HEN1, mediates 2'-O-methylation of Piwi- interacting RNAs at their 3' ends. *Genes Dev.* **21**, 1603–1608 (2007).

29. Gunawardane, L. S. et al. A slicer-mediated mechanism for repeat-associated siRNA 5' end formation in Drosophila. *Science* **315**, 1587–1590 (2007).

30. Klattenhoff, C. et al. The Drosophila HP1 homolog Rhino is required for transposon silencing and piRNA production by dual-strand clusters. *Cell* **138**, 1137–1149 (2009).

31. Le Thomas, A. et al. Transgenerationally inherited piRNAs trigger piRNA biogenesis by changing the chromatin of piRNA clusters and inducing precursor processing. *Genes Dev.* **28**, 1667–1680 (2014).

32. Mohn, F., Sienski, G., Handler, D. & Brennecke, J. The rhino-deadlock-cutoff complex licenses noncanonical transcription of dual-strand piRNA clusters in Drosophila. *Cell* **157**, 1364–1379 (2014).

33. Andersen, P. R., Tirian, L., Vunjak, M. & Brennecke, J. A heterochromatin-dependent transcription machinery drives piRNA expression. *Nature* **549**, 54–59 (2017).
34. Zhang, Z. et al. The HP1 homolog rhino anchors a nuclear complex that suppresses piRNA precursor splicing. *Cell* **157**, 1353–1363 (2014).

35. Zhang, F. et al. UAP56 couples piRNA clusters to the perinuclear transposon silencing machinery. *Cell* **151**, 871–884 (2012).

36. Mohn, F., Handler, D. & Brennecke, J. Noncoding RNA. piRNA-guided slicing specifies transcripts for Zucchini-dependent, phased piRNA biogenesis. *Science* **348**, 812–817 (2015).

37. Han, B. W., Wang, W., Li, C., Weng, Z. & Zamore, P. D. Noncoding RNA. piRNA-guided transposon cleavage initiates Zucchini-dependent, phased piRNA production. *Science* **348**, 817–821 (2015).

38. Chiquoine, A. D. The identification, origin, and migration of the primordial germ cells in the mouse embryo. *Anat. Rec.* **118**, 135–146 (1954).

39. Ohinata, Y. et al. A signaling principle for the specification of the germ cell lineage in mice. *Cell* **137**, 571–584 (2009).

40. Yu, T. et al. Long first exons and epigenetic marks distinguish conserved pachytene piRNA clusters from other mammalian genes. *Nat. Commun.* **12**, 73 (2021).

41. Schoenberg, D. R. & Maquat, L. E. Regulation of cytoplasmic mRNA decay. *Nat. Rev. Genet.* **13**, 246–259 (2012).

42. Kervestin, S. & Jacobson, A. NMD: a multifaceted response to premature translational termination. *Nat. Rev. Mol. Cell Biol.* **13**, 700–712 (2012).

43. Graille, M. & Séraphin, B. Surveillance pathways rescuing eukaryotic ribosomes lost in translation. *Nat. Rev. Mol. Cell Biol.* **13**, 727–735 (2012).

44. Shoemaker, C. J. & Green, R. Translation drives mRNA quality control. *Nat. Struct. Mol. Biol.* **19**, 594–601 (2012).

45. Sun, Y. H. et al. Ribosomes guide pachytene piRNA formation on long intergenic piRNA precursors. *Nat. Cell Biol.* **22**, 200–212 (2020).

46. Pöyry, T. A., Kaminski, A. & Jackson, R. J. What determines whether mammalian ribosomes resume scanning after translation of a short upstream open reading frame. *Genes Dev.* **18**, 62–75 (2004).

47. Han, B. W. & Zamore, P. D. piRNAs. *Curr. Biol.* **24**, R730–R733 (2014).

48. Wu, P. H. & Zamore, P. D. Defining the functions of PIWI-interacting RNAs. *Nat. Rev. Mol. Cell Biol.* **22**, 239–240 (2021).

49. Kim, V. N. Sorting out small RNAs. *Cell* **133**, 25–26 (2008).

50. Grimson, A. et al. Early origins and evolution of microRNAs and Piwi-interacting RNAs in animals. *Nature* **455**, 1193–1197 (2008).

51. Goh, W. S. et al. piRNA-directed cleavage of meiotic transcripts regulates spermatogenesis. *Genes Dev.* **29**, 1032–1044 (2015).

52. Gou, L. T. et al. Pachytene piRNAs instruct massive mRNA elimination during late spermiogenesis. *Cell Res.* **24**, 680–700 (2014).

53. Zhang, P. et al. MIWI and piRNA-mediated cleavage of messenger RNAs in mouse testes. *Cell Res.* **25**, 193–207 (2015).

54. Wu, P. H. et al. The evolutionarily conserved piRNA-producing locus pi6 is required for male mouse fertility. *Nat. Genet.* **52**, 728–739 (2020).

55. Saito, K. et al. A regulatory circuit for piwi by the large Maf gene traffic jam in Drosophila. *Nature* **461**, 1296–1299 (2009).

56. Robine, N. et al. A broadly conserved pathway generates 3'UTR-directed primary piRNAs. *Curr. Biol.* **19**, 2066–2076 (2009).

57. Chirn, G. W. et al. Conserved piRNA expression from a distinct set of piRNA cluster loci in eutherian mammals. *PLoS Genet.* **11**, e1005652 (2015).

58. Bolcun-Filas, E. et al. A-MYB (MYBL1) transcription factor is a master regulator of male meiosis. *Development* **138**, 3319–3330 (2011).

59. Yu, T. et al. The piRNA response to retroviral invasion of the koala genome. *Cell* **179**, 632.e12–643.e12 (2019).

60. Sun, Y.H., Wang, A., Song, C. et al. Single-molecule long-read sequencing reveals a conserved intact long RNA profile in sperm. *Nat. Commun.* **12**, 1361 (2021).

61. Gudlaugsdottir, S., Boswell, D. R., Wood, G. R. & Ma, J. Exon size distribution and the origin of introns. *Genetica* **131**, 299–306 (2007).

62. Li, W. et al. Alternative cleavage and polyadenylation in spermatogenesis connects chromatin regulation with post-transcriptional control. *BMC Biol.* **14**, 6 (2016).

63. Bao, J. et al. UPF2-dependent nonsense-mediated mRNA decay pathway is essential for spermatogenesis by selectively eliminating longer 3'UTR transcripts. *PLoS Genet.* **12**, e1005863 (2016).

64. Zhang, Y. et al. MicroRNAs control mRNA fate by compartmentalization based on 3' UTR length in male germ cells. *Genome Biol.* **18**, 105 (2017).

65. Tang, C. et al. ALKBH5-dependent m6A demethylation controls splicing and stability of long 3'-UTR mRNAs in male germ cells. *Proc. Natl Acad. Sci. USA* **115**, E325–E333 (2018).

66. Bohlen, J., Fenzl, K., Kramer, G., Bukau, B. & Teleman, A. A. Selective 40S footprinting reveals cap-tethered ribosome scanning in human cells. *Mol. Cell* **79**, 561.e5–574.e5 (2020).

67. Leppek, K., Das, R. & Barna, M. Functional 5' UTR mRNA structures in eukaryotic translation regulation and how to find them. *Nat. Rev. Mol. Cell Biol.* **19**, 158–174 (2018).

68. Tuller, T., Waldman, Y. Y., Kupiec, M. & Ruppin, E. Translation efficiency is determined by both codon bias and folding energy. *Proc. Natl Acad. Sci. USA* **107**, 3645–3650 (2010).

69. Sharp, P. M. & Li, W. H. The codon Adaptation Index—a measure of directional synonymous codon usage bias, and its potential applications. *Nucleic Acids Res.* **15**, 1281–1295 (1987).

70. Xia, X. An improved implementation of codon adaptation index. *Evol. Bioinform. Online* **3**, 53–58 (2007).

71. Puigbó, P., Bravo, I. G. & Garcia-Vallvé, S. E-CAI: a novel server to estimate an expected value of Codon Adaptation Index (eCAI). *BMC Bioinformatics* **9**, 65 (2008).

72. Kim, S. ppcor: An R package for a fast calculation to semi-partial correlation coefficients. *Commun. Stat. Appl. Methods* **22**, 665–685 (2015).

73. Maquat, L. E., Tarn, W. Y. & Isken, O. The pioneer round of translation: features and functions. *Cell* **142**, 368–374 (2010).

74. Ishigaki, Y., Li, X., Serin, G. & Maquat, L. E. Evidence for a pioneer round of mRNA translation: mRNAs subject to nonsense-mediated decay in mammalian cells are bound by CBP80 and CBP20. *Cell* **106**, 607–617 (2001).

75. Trcek, T., Sato, H., Singer, R. H. & Maquat, L. E. Temporal and spatial characterization of nonsense-mediated mRNA decay. *Genes Dev.* **27**, 541–551 (2013).

76. Lejeune, F., Ishigaki, Y., Li, X. & Maquat, L. E. The exon junction complex is detected on CBP80-bound but not eIF4E-bound mRNA in mammalian cells: dynamics of mRNP remodeling. *EMBO J.* **21**, 3536–3545 (2002).

77. Fresno, M., Jiménez, A. & Vázquez, D. Inhibition of translation in eukaryotic systems by harringtonine. *Eur. J. Biochem.* **72**, 323–330 (1977).

78. Ding, D. et al. TDRD5 binds piRNA precursors and selectively enhances pachytene piRNA processing in mice. *Nat. Commun.* **9**, 127 (2018).

79. Skabkin, M. A., Skabkina, O. V., Hellen, C. U. & Pestova, T. V. Reinitiation and other unconventional posttermination events during eukaryotic translation. *Mol. Cell* **51**, 249–264 (2013).

80. Guydosh, N. R. & Green, R. Dom34 rescues ribosomes in 3' untranslated regions. *Cell* **156**, 950–962 (2014).

81. Miettinen, T. P. & Björklund, M. Modified ribosome profiling reveals high abundance of ribosome protected mRNA fragments derived from 3' untranslated regions. *Nucleic Acids Res.* **43**, 1019–1034 (2015).

82. Young, D. J., Guydosh, N. R., Zhang, F., Hinnebusch, A. G. & Green, R. Rli1/ABCE1 recycles terminating ribosomes and controls translation reinitiation in 3'UTRs in vivo. *Cell* **162**, 872–884 (2015).

83. Mills, E. W., Wangen, J., Green, R. & Ingolia, N. T. Dynamic regulation of a ribosome rescue pathway in erythroid cells and platelets. *Cell Rep.* **17**, 1–10 (2016).

84. Schuller, A. P., Wu, C. C., Dever, T. E., Buskirk, A. R. & Green, R. eIF5A functions globally in translation elongation and termination. *Mol. Cell* **66**, 194.e5–205.e5 (2017).

85. Hu, J., Sun, F. & Handel, M. A. Nuclear localization of EIF4G3 suggests a role for the XY body in translational regulation during spermatogenesis in mice. *Biol. Reprod.* **98**, 102–114 (2018).

86. Passos, D. O. et al. Analysis of Dom34 and its function in no-go decay. *Mol. Biol. Cell* **20**, 3025–3032 (2009).

87. Saito, S., Hosoda, N. & Hoshino, S. The Hbs1-Dom34 protein complex functions in non-stop mRNA decay in mammalian cells. *J. Biol. Chem.* **288**, 17832–17843 (2013).

88. Hogg, J. R. & Goff, S. P. Upf1 senses 3'UTR length to potentiate mRNA decay. *Cell* **143**, 379–389 (2010).

89. Mocquet, V., Durand, S. & Jalinot, P. How retroviruses escape the nonsense-mediated mRNA decay. *AIDS Res. Hum. Retroviruses* **31**, 948–958 (2015).

90. Annibaldis, G. et al. Readthrough of stop codons under limiting ABCE1 concentration involves frameshifting and inhibits nonsense-mediated mRNA decay. *Nucleic Acids Res.* **48**, 10259–10279 (2020).

91. Schweingruber, C., Rufener, S. C., Zünd, D., Yamashita, A. & Mühlemann, O. Nonsense-mediated mRNA decay - mechanisms of substrate mRNA recognition and degradation in mammalian cells. *Biochim. Biophys. Acta* **1829**, 612–623 (2013).

92. Shum, E. Y. et al. The antagonistic gene paralogs Upf3a and Upf3b govern nonsense-mediated RNA decay. *Cell* **165**, 382–395 (2016).

93. Kurosaki, T., Popp, M. W. & Maquat, L. E. Quality and quantity control of gene expression by nonsense-mediated mRNA decay. *Nat. Rev. Mol. Cell Biol.* **20**, 406–420 (2019).

94. Lund, M. K. & Guthrie, C. The DEAD-box protein Dbp5p is required to dissociate Mex67p from exported mRNPs at the nuclear rim. *Mol. Cell* **20**, 645–651 (2005).

95. Tran, E. J., Zhou, Y., Corbett, A. H. & Wente, S. R. The DEAD-box protein Dbp5 controls mRNA export by triggering specific RNA:protein remodeling events. *Mol. Cell* **28**, 850–859 (2007).

96. Burnett, P. E., Barrow, R. K., Cohen, N. A., Snyder, S. H. & Sabatini, D. M. RAFT1 phosphorylation of the translational regulators p70 S6 kinase and 4E-BP1. *Proc. Natl Acad. Sci. USA* **95**, 1432–1437 (1998).
97. Smith, C. L. & Eppig, J. T. The mammalian phenotype ontology: enabling robust annotation and comparative analysis. *Wiley Interdiscip. Rev. Syst. Biol. Med.* **1**, 390–399 (2009).

98. Dickinson, M. E. et al. High-throughput discovery of novel developmental phenotypes. *Nature* **537**, 508–514 (2016).

99. Pires-daSilva, A. et al. Mice deficient for spermatid perinuclear RNA-binding protein show neurologic, spermatogenic, and sperm morphological abnormalities. *Dev. Biol.* **233**, 319–328 (2001).

100. Kile, B. T. et al. Functional analysis of Asb-1 using genetic modification in mice. *Mol. Cell Biol.* **21**, 6189–6197 (2001).

101. Bhandari, R., Juluri, K. R., Resnick, A. C. & Snyder, S. H. Gene deletion of inositol hexakisphosphate kinase 1 reveals inositol pyrophosphate regulation of insulin secretion, growth, and spermiogenesis. *Proc. Natl Acad. Sci. USA* **105**, 2349–2353 (2008).

102. Mu, X. et al. Targeted inactivation of testicular nuclear orphan receptor 4 delays and disrupts late meiotic prophase and subsequent meiotic divisions of spermatogenesis. *Mol. Cell Biol.* **24**, 5887–5899 (2004).

103. Morita, S. et al. One Argonaute family member, Eif2c2 (Ago2), is essential for development and appears not to be involved in DNA methylation. *Genomics* **89**, 687–696 (2007).

104. Frederick, J. P. et al. An essential role for an inositol polyphosphate multikinase, Ipk2, in mouse embryogenesis and second messenger production. *Proc. Natl Acad. Sci. USA* **102**, 8454–8459 (2005).

105. Ayadi, A. et al. Mouse large-scale phenotyping initiatives: overview of the European Mouse Disease Clinic (EUMODIC) and of the Wellcome Trust Sanger Institute Mouse Genetics Project. *Mamm. Genome* **23**, 600–610 (2012).

106. Liu, J. et al. Argonaute2 is the catalytic engine of mammalian RNAi. *Science* **305**, 1437–1441 (2004).

107. Wang, L., Beserra, C. & Garbers, D. L. A novel aminophospholipid transporter exclusively expressed in spermatozoa is required for membrane lipid asymmetry and normal fertilization. *Dev. Biol.* **267**, 203–215 (2004).

108. Guardavaccaro, D. et al. Control of meiotic and mitotic progression by the F box protein beta-Trcp1 in vivo. *Dev. Cell* **4**, 799–812 (2003).

109. Beckers, A. et al. The FOXJ1 target Ccap206 is required for sperm motility, mucociliary clearance of the airways and brain development. *Development* **147**, dev188052 (2020).

110. Bergo, M. O. et al. Defining the importance of phosphatidylylserine synthase 2 in mice. *J. Biol. Chem.* **277**, 47701–47708 (2002).

111. Ha, S., Lindsay, A. M., Timms, A. E. & Beier, D. R. Mutations in Dnaa1f and Lrrc48 cause hydrocephalus, laterality defects, and sinusitis in mice. *G3* **6**, 2479–2487 (2016).

112. Sadato, D. et al. Eukaryotic translation initiation factor 3 (eIF3) subunit e is essential for embryonic development and cell proliferation. *FEBS Open Bio* **8**, 1188–1201 (2018).

113. Zurita Rendón, O., Silva Neiva, L., Sasarman, F. & Shoubridge, E. A. The arginine methyltransferase NDUFAF7 is essential for complex I assembly and early vertebrate embryogenesis. *Hum. Mol. Genet.* **23**, 5159–5170 (2014).

114. Cao, S. et al. The high-mobility-group box protein SSRP1/T160 is essential for cell viability in day 3.5 mouse embryos. *Mol. Cell Biol.* **23**, 5301–5307 (2003).

115. Graham, D. B. et al. TMEM258 is a component of the oligosaccharyltransferase complex controlling ER stress and intestinal inflammation. *Cell Rep.* **17**, 2955–2965 (2016).

116. Wang, Y. et al. Mutation in Rpa1 results in defective DNA double-strand break repair, chromosomal instability and cancer in mice. *Nat. Genet.* **37**, 750–755 (2005).

117. Hildebrand, M. S. et al. Glucose metabolism transporters and epilepsy: only GLUT1 has an established role. *Epilepsia* **55**, e18–e21 (2014).

118. Guenther, M. G., Levine, S. S., Boyer, L. A., Jaenisch, R. & Young, R. A. A chromatin landmark and transcription initiation at most promoters in human cells. *Cell* **130**, 77–88 (2007).

119. Benton, M. J. & Donoghue, P. C. Paleontological evidence to date the tree of life. *Mol. Biol. Evol.* **24**, 26–53 (2007).

120. Aravin, A. A., Sachidanandam, R., Girard, A., Fejes-Toth, K. & Hannon, G. J. Developmentally regulated piRNA clusters implicate MILI in transposon control. *Science* **316**, 744–747 (2007).

121. Aravin, A. A. et al. A piRNA pathway primed by individual transposons is linked to de novo DNA methylation in mice. *Mol. Cell* **31**, 785–799 (2008).

122. Mouse, G. S. C. et al. Initial sequencing and comparative analysis of the mouse genome. *Nature* **420**, 520–562 (2002).

123. International Chicken Genome Sequencing Consortium. Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. *Nature* **432**, 695–716 (2004).

124. Ernst, C., Odom, D. T. & Kutter, C. The emergence of piRNAs against transposon invasion to preserve mammalian genome integrity. *Nat. Commun.* **8**, 1411 (2017).

125. Wang, C. & Lin, H. Roles of piRNAs in transposon and pseudogene regulation of germline mRNAs and lncRNAs. *Genome Biol.* **22**, 27 (2021).

126. Di Giacomo, M. et al. Multiple epigenetic mechanisms and the piRNA pathway enforce LINE1 silencing during adult spermatogenesis. *Mol. Cell* **50**, 601–608 (2013).

127. Watanabe, T., Cheng, E. C., Zhong, M. & Lin, H. Retrotransposons and pseudogenes regulate mRNAs and lncRNAs via the piRNA pathway in the germline. *Genome Res.* **25**, 368–380 (2015).

128. Sun, Y. H. et al. Domestic chickens activate a piRNA defense against avian leukosis virus. *Elife* **6**, e24695 (2017).

129. Schuller, A. P. & Green, R. Roadblocks and resolutions in eukaryotic translation. *Nat. Rev. Mol. Cell Biol.* **19**, 526–541 (2018).

130. Dever, T. E. & Green, R. The elongation, termination, and recycling phases of translation in eukaryotes. *Cold Spring Harb. Perspect. Biol.* **4**, a013706 (2012).

131. Zhu, X., Zhang, H. & Mendell, J. T. Ribosome recycling by ABCE1 links lysosomal function and iron homeostasis to 3' UTR-directed regulation and nonsense-mediated decay. *Cell Rep.* **32**, 107895 (2020).

132. Young, D. J. et al. Tma64/eIF2D, Tma20/MCT-1, and Tma22/DENR recycle post-termination 40S subunits in vivo. *Mol. Cell* **71**, 761.e5–774.e5 (2018).

133. Bohlen, J. et al. DENR promotes translation reinitiation via ribosome recycling to drive expression of oncogenes including ATF4. *Nat. Commun.* **11**, 4676 (2020).

134. Young, D. J., Meydan, S. & Guydosh, N. R. 40S ribosome profiling reveals distinct roles for Tma20/Tma22 (MCT-1/DENR) and Tma64 (eIF2D) in 40S subunit recycling. *Nat. Commun.* **12**, 2976 (2021).

135. Howe, K. L. et al. Ensembl 2021. *Nucleic Acids Res.* **49**, D884–D891 (2021).

136. Homolka, D. et al. PIWI slicing and RNA elements in precursors instruct directional primary piRNA biogenesis. *Cell Rep.* **12**, 418–428 (2015).

137. Ishizu, H. et al. Somatic primary piRNA biogenesis driven by cis-acting RNA elements and trans-acting Yb. *Cell Rep.* **12**, 429–440 (2015).

138. Sanz, E. et al. Cell-type-specific isolation of ribosome-associated mRNA from complex tissues. *Proc. Natl Acad. Sci. USA* **106**, 13939–13944 (2009).

139. Schonhoff, S. E., Giel-Moloney, M. & Leiter, A. B. Neurogenin 3-expressing progenitor cells in the gastrointestinal tract differentiate into both endocrine and non-endocrine cell types. *Dev. Biol.* **270**, 443–454 (2004).

140. Morlan, J. D., Qu, K. & Sinicropi, D. V. Selective depletion of rRNA enables whole transcriptome profiling of archival fixed tissue. *PLoS ONE* **7**, e42882 (2012).

141. Adiconis, X. et al. Comparative analysis of RNA sequencing methods for degraded or low-input samples. *Nat. Methods* **10**, 623–629 (2013).

142. Page, J., Suja, J. A., Santos, J. L. & Rufas, J. S. Squash procedure for protein immunolocalization in meiotic cells. *Chromosome Res.* **6**, 639–642 (1998).

143. Aravin, A. et al. A novel class of small RNAs bind to MILI protein in mouse testes. *Nature* **442**, 203–207 (2006).

144. Han, B. W., Wang, W., Zamore, P. D. & Weng, Z. PiPipes: a set of pipelines for piRNA and transposon analysis via small RNA-seq, RNA-seq, degradome- and CAGE-seq, ChIP-seq and genomic DNA sequencing. *Bioinformatics* **31**, 593–595 (2015).

145. Hinrichs, A. S. et al. The UCSC Genome Browser Database: update 2006. *Nucleic Acids Res.* **34**, D590–D598 (2006).

146. Trapnell, C. et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. *Nat. Protoc.* **7**, 562–578 (2012).

147. Patro, R., Duggal, G., Love, M. I., Irizarry, R. A. & Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. *Nat. Methods* **14**, 417–419 (2017).

148. Ingolia, N. T. Ribosome footprint profiling of translation throughout the genome. *Cell* **165**, 22–33 (2016).

149. Zhang, Y. et al. Model-based analysis of ChIP-Seq (MACS). *Genome Biol.* **9**, R137 (2008).

150. Trapnell, C., Pachter, L. & Salzberg, S. L. TopHat: discovering splice junctions with RNA-Seq. *Bioinformatics* **25**, 1105–1111 (2009).

151. R Core Team. *R: A Language and Environment for Statistical Computing* (R Foundation for Statistical Computing, 2014).

152. Pertea, M. et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. *Nat. Biotechnol.* **33**, 290–295 (2015).

153. Haas, B. J. et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. *Nat. Protoc.* **8**, 1494–1512 (2013).

154. Wichert, S., Fokianos, K. & Strimmer, K. Identifying periodically expressed transcripts in microarray time series data. *Bioinformatics* **20**, 5–20 (2004).

155. Xia, X. DAMBE7: new and improved tools for data analysis in molecular biology and evolution. *Mol. Biol. Evol.* **35**, 1550–1552 (2018).

156. Puigbò, P., Bravo, I. G. & Garcia-Vallve, S. CAIcal: a combined set of tools to assess codon usage adaptation. *Biol. Direct* **3**, 38 (2008).

157. Betel, D., Koppal, A., Agius, P., Sander, C. & Leslie, C. Comprehensive modeling of microRNA targets predicts functional non-conserved and non-canonical sites. *Genome Biol.* **11**, R90 (2010).

158. Rueden, C. T. et al. ImageJ2: ImageJ for the next generation of scientific image data. *BMC Bioinformatics* **18**, 529 (2017).
## Acknowledgements

We thank L. Maquat, T. Kurosaki, D. Ermolenko, J. Lin, X. Zhuo, and members of the Li laboratory for discussions; G. Riddihough and J. Lozada for help with editing the manuscript; and J. P. Wang, UR pathology core, URMC Mass Spectrometry Resource Facility, and UR Genomics Research Center for help with the experiments. We thank N. Chen for creating the cartoon pictures. This work was supported by National Institutes of Health grant R35GM128782 and Agriculture and Food Research Initiative Competitive Grant no. 2018-67015-27615 from the USDA National Institute of Food and Agriculture to X.Z.L.

## Author contributions

Y.H.S. and J. Zheng analyzed the data with input from C.Z., E.P.R., and X.Z.L.; R.H.W, K.D., J. Zhu, L.H.X., and A.A.P. performed the experiments with input from E.P.R. and X.Z.L. X.Z.L. designed the study and drafted the manuscript, and all authors contributed to the preparation of the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-021-26233-8.

**Correspondence** and requests for materials should be addressed to Xin Zhiguo Li.

**Peer review information** *Nature Communications* thanks Haruhiko Siomi and Shu-Bing Qian for their contributions to the peer review of this work. Peer review reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_22.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021