An integrative piRNA analysis of mouse gametes and zygotes reveals new potential origins and gene regulatory roles

Piwi-interacting RNAs (piRNAs) are a subclass of the small non-coding RNAs (sncRNAs). Their main reported function was to exert control over transposable elements (TEs) in mammalian germline. In this study undertaking a deeper bioinformatics analysis of piRNAs present in mouse oocytes, sperm cells and zygotes, we first elaborated a new piRNA database based on sequences identified as piRNAs by immunoprecipitation with PIWI proteins. Our bioinformatics analysis revealed that, at least in gametes and zygotes, piRNAs could encompass multifunctional cell-dependent regulatory molecules. Indeed, genome analysis of the piRNA mapping density (reads/kb) evidenced in all samples an enrichment of intron-derived piRNAs. Further, piRNA population was classified into sequences not associated to TEs or repeats (NRapiRNAs) and associated to repetitive genome elements (RapiRNAs). In oocytes most of the NRapiRNAs mapped to the 5′UTRs of coding mRNAs, while higher proportion of NRapiRNAs was detected in sperm cells associated to the 3′UTRs of mRNAs. This piRNA complementarity to mRNA UTRs suggests key post-transcriptional regulatory roles over mRNAs such as those encoding MHC genes. In addition, a striking association of RapiRNA with long non-coding RNAs (lncRNAs) was identified. piRNAs associated with relevant lncRNAs such as: Rab26os and GAS5 and key mRNAs, were particularly assessed.