The characterisation of piRNA-related 19mers in the mouse
© Oey et al; licensee BioMed Central Ltd. 2011
Received: 25 February 2011
Accepted: 15 June 2011
Published: 15 June 2011
Piwi interacting RNA, or piRNA, is a class of small RNA almost exclusively expressed in the germline where they serve essential roles in retrotransposon silencing. There are two types, primary and secondary piRNA, and the latter is a product of enzymatic cleavage of retrotransposons' transcripts directed by the former. Recently, a new class of 19nt long RNA was discovered that is specific to testis and appears to be linked to secondary piRNA biogenesis.
We locate clusters of the testis-specific 19mers, which we call piRNA-related 19mers (pr19RNA), and characterise the transcripts from which they are derived. Most pr19RNA clusters were associated with retrotransposons and unannotated antisense transcripts overlapping piRNA clusters. At these loci the abundance of 19mers was found to be greater than that of secondary piRNAs.
We find that pr19RNAs are distinguished from other RNA populations by their length and flanking sequence, allowing their identification without requiring overlapping piRNAs. Using such sequence features allows identification of the source transcripts, and we suggest that these likely represent the substrates of primary piRNA-guided RNA cleavage events. While pr19RNAs appear not to bind directly to Miwi or Mili, their abundance relative to secondary piRNAs, in combination with their precise length, suggests they may be more than by-products of secondary piRNA biogenesis.
Piwi-interacting small RNAs (piRNAs) are almost exclusively expressed in the germline where they serve essential roles in retrotransposon silencing [1–5]. There are two types, the primary piRNAs (p-piRNAs) and secondary piRNAs (s-piRNAs), and they associate with three developmentally regulated Ago-related proteins, Miwi, Mili and Miwi2, in mice . The p-piRNAs are identified by their length, spanning ~24-30 ribonucleotides (nt) and the presence of a 5' U. They are mostly derived from long non-coding RNAs transcribed from distinct genomic loci known as piRNA clusters [5, 7].
Some p-piRNAs target the transcripts of active retrotransposons and retroviruses, such as the long interspersed nuclear elements (LINEs) and the intracisternal A particles (IAPs). The retrotransposon transcripts are then cleaved enzymatically resulting in the release of s-piRNAs. The initial cleavage site is located within the base-paired region 10 nt downstream of the 5' terminal U of the p-piRNA and the resulting s-piRNAs are therefore distinguished by an A at position 10 [8–10]. Some s-piRNAs are expected to be reverse complementary to the original p-piRNA precursor transcript and may themselves be able to direct cleavage of these to recreate the original p-piRNA. The formation of s-piRNAs directed by p-piRNAs and p-piRNAs directed by s-piRNAs has been been called the ping pong amplification cycle . The ping-pong amplification, identified by the presence of secondary piRNAs, is particularly active at the pre-pachytene stage. At the pachytene stage the population of piRNAs shifts towards p-piRNAs [10, 11]. The precise nature of these events in mice is still being investigated and appears to vary according to the complement of Piwi-related proteins that are expressed .
In this report we study a class of 19 nt long piRNA-related RNAs (pr19RNA) in mice testes that appears to result from piRNA directed cleavage of transcripts. While this manuscript was in preparation these RNAs were also reported by Berninger et al. . However, our analyses differ in that we investigate the genomic origin of most of these 19mers, including those derived from the repetitive IAP and LINE elements, rather than limiting the analyses to non-repetitive RNA. We find that the 19mers are readily distinguished from other RNAs by the composition of their downstream flanking sequence and use this feature to identify likely substrates of piRNA-directed RNA cleavage.
Deep sequencing datasets and initial data processing
Three publicly available spermatogenic tubule small RNA deep sequencing datasets published by Robine et al . were obtained from the Gene Expression Omnibus (GEO)  under the accession number GSM475279, GSM475280 and GSM475281, and used for the primary analyses. Additionally, small RNA sequencing data from several tissues published by Chiang et al . were also obtained under the GEO Series accession GSE20384 and used for some analyses, as specified in the results. Prior to mapping, all reads containing homopolymers > 8 nt were removed and reads matching the partial sequence CCGGGTTTCGGCACC, identified as tRNA-derived, were also removed. In the case of reads mapping to LINE elements, reads matching the patterns (GGA)4, (CCT)4, (GAA)4 and (CTT)4 were also removed as these were found to cause mapping artefacts against common simple repeats inside many LINE elements. All other data processing and sequence analyses were carried out using scripts available as part of the Biopieces (http://www.biopieces.org).
Mapping of reads to LINE and IAP
Reads likely to originate from IAP or LINE elements were identified by mapping the reads to a database containing all mouse IAP elements and one containing all LINE elements annotated in the UCSC Genome Browser RepeatMasker track (accessed 01.09.2010) . To obtain reads likely to originate from IAP or LINE repeats the program Bowtie  was used to map reads directly to the IAP and LINE sequences using the options "-v 0 -m 50000" to get high-confidence IAP and LINE small RNA datasets. These reads were then remapped to IAP or LINE consensus sequences with less stringent parameters allowing for two mismatches (Bowtie option -v 2). The consensus sequences that were used were obtained from RepBase  under the identifiers "IAPLTR1_MM_LTR" and "IAPLTR1_MM" for IAP and "L1_MM" for LINE. The two IAP reference sequences were joined such that the IAPLTR1_MM sequence was flanked on both sides by a copy of the IAPLTR1_Mm_LTR sequence to obtain an IAP reference sequence that includes LTRs. The sequence composition of pr19RNA and piRNA reads were visualised in the form of sequence logos that were plotted using the "plot_seqlogo" script, which is available as part of the Biopieces (http://www.biopieces.org). The script calculates logos using Shannons general formula for uncertainty [18, 19].
Genome-wide annotation of 19mers
All 19mers from GSM475281 were mapped to the mouse genome (NCBI37/mm9) using Bowtie  allowing reads to map to < 6 loci. Reads mapping to multiple loci were randomly assigned to single loci. The sequence of all the reads plus 10nt of downstream flanking sequence was obtained, and the reads were binned according to the identity of the 10th downstream base. The origin of the reads were then found by intersection against UCSC Genome Browser RepeatMasker annotation tracks followed by intersection with piRNA clusters, using 94 piRNA clusters annotated by Lau et al. , and finally with RefSeq genes, which were also obtained through the UCSC Genome Browser .
Positional relationship of IAP-derived p-piRNAs and s-piRNAs
The tendency of p-piRNAs and s-piRNAs to be positionally fixed to pr19RNAs was investigated by first identifying high-confidence pr19RNAs as well as p-piRNAs and s-piRNAs. The pr19RNAs were defined as reads in same orientation as the IAP reference sequence with the sequence motif N19VN8A, where N19 represents the pr19RNA and V represents any nucleotide except U. The p-piRNAs were reads in the opposite sense of the IAP reference sequence, 24-30 nt in length and with the sequence motif UN8BNn, where B represents any nucleotide except A. The s-piRNAs were any read in the same orientation as the IAP reference sequence, 24-30 nt in length and with the sequence motif VN8ANn, where V represents any nucleotide except U. The 5' ends of the s-piRNA 5' ends and p-piRNA 5' ends were then tallied at incrementally increasing distances from pr19RNAs. The positional relationship between these were elucidated by plotting the density of piRNA 5' ends relative to pr19RNA 5' ends, similar to the method employed by Aravin et al. .
Distribution of prRNA relative to s-piRNAs across the genome
The pr19RNA, s-piRNAs and p-piRNA sequence patterns were also used to define regions in the mouse genome enriched for pr19RNA and s-piRNAs by selecting those regions where uniquely mapping pr19RNAs outnumber other 19mers, and where uniquely mapping s-piRNAs outnumber p-piRNAs. Regions enriched for both were used to estimate the typical ratio of pr19RNAs relative to s-piRNAs genome-wide.
Miwi and Mili IP analyses
The IP datasets were mapped to the LINE and IAP consensus sequences as well as miRNAs (miRBase v.16) , tRNAs , snoRNAs [23, 24] and the mitochondrial genome (NCBI37/mm9). The number of mapped reads were then normalised to the number of genomically mappable reads in each library as reads per million.
Strand-specific real-time quantitative reverse transcription PCR
RNA was extracted from 4-month old C57BL/6 mouse testes with Trizol (Invitrogen), 2 μg was treated with DNase I (New England Biolabs) and used for reverse transcription (RT). RT was carried out using Superscript III (Invitrogen) according to manufacturer's instructions, apart from the following changes; MgCl2 was added to a final concentration of 5 mM and the RT step was performed at 50°C. Two strand-specific RT reactions were carried out against each strand of the piRNA cluster using the forward and reverse primers to prime the first-strand synthesis (Target1-forward: GCCTATTGCTCCTTGGACTG, Target1-reverse: AGCTAGCTGAAATCGGATGG, Target2-forward: AGGAAAAAGGTTGCCGATCT, Target2-reverse: CTCTTGAAAGGGCTTCTTGC. GAPDH was reverse transcribed in parallel and used to normalise between the reactions (GAPDH-forward: TCCAGAGACAGCCGCATCT, GAPDH-reverse: ACACCGACCTTCACCATTTTG). Real-time quantitative RT PCR (qRT-PCR) was performed in triplicate using Platinum SYBR Green qPCR mix (Invitrogen) on a Corbett Rotorgene qPCR machine at the default settings. Differences in expression between strands were estimated using the 2-ΔΔCT method .
Small RNAs derived from IAPs in testis
Because s-piRNAs have a known preference for A at position 10, it is likely that the 5' ends of the 19mers were derived from a cleavage event 19 nt upstream of the 5' termini of these s-piRNAs. This suggested that many of the 19mers would be flanked by s-piRNAs at their 3' end and overlap with opposite sense p-piRNAs. We tested this hypothesis by tallying the 5' termini of same-sense s-piRNAs and opposite-sense p-piRNAs relative to the 5' ends of the 19mers. As predicted, we found this to be the case (Figure 1D and Additional file 1 Figure S1A). This positional relationship between the 19mers and piRNAs was also described by Berninger et al.  using non-repetitive reads and it mirrors the relationship between p-piRNAs and s-piRNAs previously described by others [9, 11]. The correlation in the distribution of 19mers relative to p-piRNAs and and s-piRNAs is also evident when plotting the density of these along the IAP reference sequence where the 19mers are often found to overlap with same-sense s-piRNAs and opposite sense p-piRNAs (Figure 1C). Out of the 2,743 reads identified as s-piRNAs (RNAs 24-30nt in length with A in position 10 and absence of U in position 1) mapped to the sense-strand of the IAP consensus sequence, 1,931 (70%) were found immediately adjacent to a 19mer (configuration illustrated in Figure 1D). Similarly, out of 4,985 p-piRNAs (RNAs 24-30nt in length with U in position 1 and absence of A in position 10) mapped in antisense to the IAP consensus sequence, 4,434 (89%) were found to have 5' termini located precisely 10 bases downstream of the 3' terminus of a 19mer mapped to IAP in sense (configuration illustrated in Figure 1D).
It did not escape our notice that some 19mers had a strong bias towards U immediately downstream of their 3' ends (Figure 1B). We found that many of these 19mers were positionally locked to the 5' end of same-sense p-piRNAs while they overlapped with opposite sense s-piRNA (Additional file 1 Figure S1A). This suggests that this latter group is formed from enzymatic cleavage of IAP guided by s-piRNAs, thus causing the release of a 19mer with a nucleotide bias characterised by a U immediately downstream of the 3' terminus, instead of an A 10nt downstream. We tallied the number of 19mers that were either immediately adjacent to a piRNA-sized read, or antisense to one, in the two possible configurations indicated in Figure 1D and Additional file 1 Figure S1A, but disregarding whether the piRNAs where p-piRNAs or s-piRNAs. We found that out of the 9,648 19mers mapped to the sense strand, 8,861 (91%) were found with at least one piRNA in such a configuration, while out of 8,139 19mers mapped in antisense, 7,238 (89%) were arranged in this manner. The remaining IAP-derived 19mers have minimal downstream sequence bias, and likely represent RNA degradation products. As the 19mers seem likely to be formed via a piRNA-dependent mechanism, we call them piRNA-related 19mers, or pr19RNAs.
We also investigated the full complement of small RNAs that mapped to any mouse genomic IAP element in our database, as many were not mappable to the consensus IAP sequence. This increased the number of IAP-derived reads two-fold but had only a minor effect on the distribution of 19mers relative to piRNAs (Additional file 1 Figure S1B and S1C). The number of s-piRNAs in the database was 7,700 (defined as reads 24-30 nt in length, the absence of 5'U and the presence of A in position 10), while that of pr19RNAs was 8,839 (defined as 19mers with an A 10nt downstream and absence of U immediately downstream of the 3' termini).
Small RNAs mapped to LINE elements
Most 19mers in mouse spermatogenic tubules are pr19RNAs
Origin of non-repeat associated 19mers
Little evidence for binding of pr19RNA to Miwi or Mili
pr19RNAs in somatic tissues
In this report we characterise novel members of a recently discovered class of 19nt long RNAs that originate from piRNA-directed RNA cleavage. In a previous paper Berninger et al.  elegantly show that these small RNAs can be distinguished by their spatial relationship to p-piRNAs and s-piRNAs. As expected, we find that this class of RNA is readily distinguished from other RNAs by the composition of the downstream flanking sequence. By making use of sequence composition rather than overlap with piRNAs a greater number of pr19RNAs can be identified. Their precise length suggests they are formed from two separate enzymatic cleavage events with the 3' end likely to be linked to the formation of s-piRNA 5' ends and the 5' end likely to be representing a novel event. The biogenesis of s-piRNAs has been hypothesised to require only a single enzymatic cleavage event at the 5' end, with the 3' end defined by a less specific mechanism resulting in their variable lengths .
Curiously, the IAP and LINE elements, as well as the reverse strand of many piRNA clusters, were found to have a greater number of associated pr19RNAs than s-piRNAs. This difference is puzzling as their biogenesis appears linked and the latter is expected to accumulate more readily due to their association with Miwi and Mili. It may suggest that subsequent steps in the formation of these RNAs are not synchronised. Alternatively, pr19RNAs may themselves be protected from degradation via a mechanism independent of Miwi or Mili as neither appear to associate strongly with pr19RNAs. The abundance of pr19RNAs mappable to the LINE consensus sequence is particularly interesting as they also outnumber p-piRNAs. It has been shown that miRNAs are incorporated into multiple-turnover enzymes such that a single microRNA can guide several consecutive cleavage events . If p-piRNAs function in a similar manner, then the number of pr19RNAs would depend on the abundance of substrate (e.g. LINE transcripts) rather than the abundance of p-piRNAs.
By locating regions enriched for 19mers with an A 10nt downstream we find several genes with a reciprocal pattern in the distribution of 19mers and p-piRNAs. The transcripts from some of these genes appear to be the the target of piRNA-guided cleavage resulting in release of a 19mer. A similar mechanism where p-piRNAs target the transcript of a protein-coding gene was recently proposed for piRNAs in Drosophila . This process appears analogous to the endogenous short interfering RNAs that operate in some mammalian tissues [28, 32, 33]. However, most pr19RNAs co-localising with genes were found to occur in the opposite sense. These 19mers likely arise when genic piRNAs, which are abundant in testis , direct the cleavage of overlapping antisense transcripts. The purpose of such events is not immediately clear.
Based on the sequence composition of the pr19RNA flanking regions and the positional relationship of pr19RNAs to p-piRNAs and s-piRNAs, we conclude that pr19RNAs are produced by enzymatic cleavage events guided mainly by p-piRNAs in adult mouse spermatogenic tubules. This was also the conclusion reached by Berninger et al.  in their analyses of overlapping piRNAs and pr19RNAs. Like many newly described classes of small RNAs, a function for pr19RNAs has not yet been ascertained. However, the 5' cleavage site, precise length and abundance of these RNAs suggest the underlying mechanism serves an important role in piRNA-directed post-transcriptional gene regulation. The characterisation of this novel class of small RNA offers new insight into the biogenesis and function of piRNAs in the germline and is likely to facilitate new methods for elucidating the targets of these enigmatic small RNAs.
We wish to thank Martin Asser Hansen for helpful bioinformatics advice and comments on the manuscript, and Lutz Krause and Lucia Clemens-Daxinger for providing helpful comments. This work was supported by NHMRC Australia.
- Aravin A, Gaidatzis D, Pfeffer S, Lagos-Quintana M, Landgraf P, Iovino N, Morris P, Brownstein MJ, Kuramochi-Miyagawa S, Nakano T, Chien M, Russo JJ, Ju J, Sheridan R, Sander C, Zavolan M, Tuschl T: A novel class of small RNAs bind to MILI protein in mouse testes. Nature. 2006, 442: 203-207.PubMedGoogle Scholar
- Grivna ST, Beyret E, Wang Z, Lin H: A novel class of small RNAs in mouse spermatogenic cells. Genes Dev. 2006, 20: 1709-1714. 10.1101/gad.1434406.View ArticlePubMedPubMed CentralGoogle Scholar
- Girard A, Sachidanandam R, Hannon GJ, Carmell MA: A germline-specific class of small RNAs binds mammalian Piwi proteins. Nature. 2006, 442: 199-202.PubMedGoogle Scholar
- Watanabe T, Takeda A, Tsukiyama T, Mise K, Okuno T, Sasaki H, Minami N, Imai H: Identification and characterization of two novel classes of small RNAs in the mouse germline: retrotransposon-derived siRNAs in oocytes and germline small RNAs in testes. Genes Dev. 2006, 20: 1732-1743. 10.1101/gad.1425706.View ArticlePubMedPubMed CentralGoogle Scholar
- Ro S, Park C, Song R, Nguyen D, Jin J, Sanders KM, McCarrey JR, Yan W: Cloning and expression profiling of testis-expressed piRNA-like RNAs. RNA. 2007, 13: 1693-1702. 10.1261/rna.640307.View ArticlePubMedPubMed CentralGoogle Scholar
- Thomson T, Lin H: The biogenesis and function of PIWI proteins and piRNAs: progress and prospect. Annu Rev Cell Dev Biol. 2009, 25: 355-376. 10.1146/annurev.cellbio.24.110707.175327.View ArticlePubMedPubMed CentralGoogle Scholar
- Robine N, Lau NC, Balla S, Jin Z, Okamura K, Kuramochi-Miyagawa S, Blower MD, Lai EC: A broadly conserved pathway generates 3'UTR-directed primary piRNAs. Curr Biol. 2009, 19: 2066-2076. 10.1016/j.cub.2009.11.064.View ArticlePubMedPubMed CentralGoogle Scholar
- Gunawardane LS, Saito K, Nishida KM, Miyoshi K, Kawamura Y, Nagami T, Siomi H, Siomi MC: A slicer-mediated mechanism for repeat-associated siRNA 5' end formation in Drosophila. Science (80-). 2007, 315: 1587-1590.View ArticleGoogle Scholar
- Brennecke J, Aravin AA, Stark A, Dus M, Kellis M, Sachidanandam R, Hannon GJ: Discrete small RNA-generating loci as master regulators of transposon activity in Drosophila. Cell. 2007, 128: 1089-1103. 10.1016/j.cell.2007.01.043.View ArticlePubMedGoogle Scholar
- Aravin AA, Sachidanandam R, Bourc'his D, Schaefer C, Pezic D, Toth KF, Bestor T, Hannon GJ: A piRNA pathway primed by individual transposons is linked to de novo DNA methylation in mice. Mol Cell. 2008, 31: 785-799. 10.1016/j.molcel.2008.09.003.View ArticlePubMedPubMed CentralGoogle Scholar
- Aravin AA, Sachidanandam R, Girard A, Fejes-Toth K, Hannon GJ: Developmentally regulated piRNA clusters implicate MILI in transposon control. Science (80-). 2007, 316: 744-747.View ArticleGoogle Scholar
- Berninger P, Jaskiewicz L, Khorshid M, Zavolan M: Conserved generation of short products at piRNA loci. BMC Genomics. 2011, 12: 46-10.1186/1471-2164-12-46.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.View ArticlePubMedPubMed CentralGoogle Scholar
- Chiang HR, Schoenfeld LW, Ruby JG, Auyeung VC, Spies N, Baek D, Johnston WK, Russ C, Luo S, Babiarz JE, Blelloch R, Schroth GP, Nusbaum C, Bartel DP: Mammalian microRNAs: experimental evaluation of novel and previously annotated genes. Genes Dev. 2010, 24: 992-1009. 10.1101/gad.1884710.View ArticlePubMedPubMed CentralGoogle Scholar
- Rhead B, Karolchik D, Kuhn RM, Hinrichs AS, Zweig AS, Fujita PA, Diekhans M, Smith KE, Rosenbloom KR, Raney BJ, Pohl A, Pheasant M, Meyer LR, Learned K, Hsu F, Hillman-Jackson J, Harte RA, Giardine B, Dreszer TR, Clawson H, Barber GP, Haussler D, Kent WJ: The UCSC Genome Browser database: update 2010. Nucleic Acids Res. 2010, 38: D613-9. 10.1093/nar/gkp939.View ArticlePubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.View ArticlePubMedPubMed CentralGoogle Scholar
- Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110: 462-467. 10.1159/000084979.View ArticlePubMedGoogle Scholar
- Schneider TD, Stephens RM: Sequence logos: a new way to display consensus sequences. Nucleic Acids Res. 1990, 18: 6097-6100. 10.1093/nar/18.20.6097.View ArticlePubMedPubMed CentralGoogle Scholar
- Shannon CE: A mathematical theory of communication. Bell System Technical Journal. 1948, 27: 379-423. 623-656View ArticleGoogle Scholar
- Lau NC, Seto AG, Kim J, Kuramochi-Miyagawa S, Nakano T, Bartel DP, Kingston RE: Characterization of the piRNA complex from rat testes. Science (80-). 2006, 313: 363-367.View ArticleGoogle Scholar
- Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, 34: D140-4. 10.1093/nar/gkj112.View ArticlePubMedGoogle Scholar
- Chan PP, Lowe TM: GtRNAdb: a database of transfer RNA genes detected in genomic sequence. Nucleic Acids Res. 2009, 37: D93-7. 10.1093/nar/gkn787.View ArticlePubMedGoogle Scholar
- Lestrade L, Weber MJ: snoRNA-LBME-db, a comprehensive database of human H/ACA and C/D box snoRNAs. Nucleic Acids Res. 2006, 34: D158-62. 10.1093/nar/gkj002.View ArticlePubMedGoogle Scholar
- Yang J, Shao P, Zhou H, Chen Y, Qu L: deepBase: a database for deeply annotating and mining deep sequencing data. Nucleic Acids Res. 2010, 38: D123-30. 10.1093/nar/gkp943.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
- Svoboda P, Stein P, Anger M, Bernstein E, Hannon GJ, Schultz RM: RNAi and expression of retrotransposons MuERV-L and IAP in preimplantation mouse embryos. Dev Biol. 2004, 269: 276-285. 10.1016/j.ydbio.2004.01.028.View ArticlePubMedGoogle Scholar
- Betel D, Sheridan R, Marks DS, Sander C: Computational analysis of mouse piRNA sequence and biogenesis. PLoS Comput Biol. 2007, 3: e222-10.1371/journal.pcbi.0030222.View ArticlePubMedPubMed CentralGoogle Scholar
- Tam OH, Aravin AA, Stein P, Girard A, Murchison EP, Cheloufi S, Hodges E, Anger M, Sachidanandam R, Schultz RM, Hannon GJ: Pseudogene-derived small interfering RNAs regulate gene expression in mouse oocytes. Nature. 2008, 453: 534-538. 10.1038/nature06904.View ArticlePubMedPubMed CentralGoogle Scholar
- Zheng K, Xiol J, Reuter M, Eckardt S, Leu NA, McLaughlin KJ, Stark A, Sachidanandam R, Pillai RS, Wang PJ: Mouse MOV10L1 associates with Piwi proteins and is an essential component of the Piwi-interacting RNA (piRNA) pathway. Proc Natl Acad Sci USA. 2010, 107: 11841-11846. 10.1073/pnas.1003953107.View ArticlePubMedPubMed CentralGoogle Scholar
- Hutvágner G, Zamore PD: A microRNA in a multiple-turnover RNAi enzyme complex. Science (80-). 2002, 297: 2056-2060.View ArticleGoogle Scholar
- Rouget C, Papin C, Boureux A, Meunier A, Franco B, Robine N, Lai EC, Pelisson A, Simonelig M: Maternal mRNA deadenylation and decay by the piRNA pathway in the early Drosophila embryo. Nature. 2010, 467: 1128-1132. 10.1038/nature09465.View ArticlePubMedPubMed CentralGoogle Scholar
- Lau NC, Robine N, Martin R, Chung W, Niki Y, Berezikov E, Lai EC: Abundant primary piRNAs, endo-siRNAs, and microRNAs in a Drosophila ovary cell line. Genome Res. 2009, 19: 1776-1785. 10.1101/gr.094896.109.View ArticlePubMedPubMed CentralGoogle Scholar
- Watanabe T, Totoki Y, Toyoda A, Kaneda M, Kuramochi-Miyagawa S, Obata Y, Chiba H, Kohara Y, Kono T, Nakano T, Surani MA, Sakaki Y, Sasaki H: Endogenous siRNAs from naturally formed dsRNAs regulate transcripts in mouse oocytes. Nature. 2008, 453: 539-543. 10.1038/nature06908.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.