Mature piRNAs Research Articles

PiRNA Defense Against Endogenous Retroviruses

Infection by retroviruses and the mobilization of transposable elements cause DNA damage that can be catastrophic for a cell. If the cell survives, the mutations generated by retrotransposition may confer a selective advantage, although, more commonly, the effect of new integrants is neutral or detrimental. If retrotransposition occurs in gametes or in the early embryo, it introduces genetic modifications that can be transmitted to the progeny and may become fixed in the germline of that species. PIWI-interacting RNAs (piRNAs) are single-stranded, 21–35 nucleotide RNAs generated by the PIWI clade of Argonaute proteins that maintain the integrity of the animal germline by silencing transposons. The sequence specific manner by which piRNAs and germline-encoded PIWI proteins repress transposons is reminiscent of CRISPR, which retains memory for invading pathogen sequences. piRNAs are processed preferentially from the unspliced transcripts of piRNA clusters. Via complementary base pairing, mature antisense piRNAs guide the PIWI clade of Argonaute proteins to transposon RNAs for degradation. Moreover, these piRNA-loaded PIWI proteins are imported into the nucleus to modulate the co-transcriptional repression of transposons by initiating histone and DNA methylation. How retroviruses that invade germ cells are first recognized as foreign by the piRNA machinery, as well as how endogenous piRNA clusters targeting the sequences of invasive genetic elements are acquired, is not known. Currently, koalas (Phascolarctos cinereus) are going through an epidemic due to the horizontal and vertical transmission of the KoRV-A gammaretrovirus. This provides an unprecedented opportunity to study how an exogenous retrovirus becomes fixed in the genome of its host, and how piRNAs targeting this retrovirus are generated in germ cells of the infected animal. Initial experiments have shown that the unspliced transcript from KoRV-A proviruses in koala testes, but not the spliced KoRV-A transcript, is directly processed into sense-strand piRNAs. The cleavage of unspliced sense-strand transcripts is thought to serve as an initial innate defense until antisense piRNAs are generated and an adaptive KoRV-A-specific genome immune response is established. Further research is expected to determine how the piRNA machinery recognizes a new foreign genetic invader, how it distinguishes between spliced and unspliced transcripts, and how a mature genome immune response is established, with both sense and antisense piRNAs and the methylation of histones and DNA at the provirus promoter.

RNA helicase D1PAS1 resolves R-loops and forms a complex for mouse pachytene piRNA biogenesis required for male fertility.

During meiosis, RNA polymerase II transcribes pachytene piRNA precursors with unusually long and unspliced transcripts from discrete autosomal loci in the mouse genome. Despite the importance of piRNA for male fertility and a well-defined maturation process, the transcriptional machinery remains poorly understood. Here, we document that D1PAS1, an ATP-dependent RNA helicase, is critical for pachytene piRNA expression from multiple genomic loci and subsequent translocation into the cytoplasm to ensure mature piRNA biogenesis. Depletion of D1PAS1 in gene-edited mice results in the accumulation of R-loops in pachytene spermatocytes, leading to DNA-damage-induced apoptosis, disruption of piRNA biogenesis, spermatogenic arrest, and male infertility. Transcriptome, genome-wide R-loop profiling, and proteomic analyses document that D1PAS1 regulates pachytene piRNA transcript elongation and termination. D1PAS1 subsequently forms a complex with nuclear export components to ensure pachytene piRNA precursor translocation from the nucleus to the cytoplasm for processing into small non-coding RNAs. Thus, our study defines D1PAS1 as a specific transcription activator that promotes R-loop unwinding and is a critical factor in pachytene piRNA biogenesis.

An evolutionarily conserved stop codon enrichment at the 5\u2032 ends of mammalian piRNAs

PIWI-interacting RNAs (piRNAs) are small RNAs required to recognize and silence transposable elements. The 5’ ends of mature piRNAs are defined through cleavage of long precursor transcripts, primarily by Zucchini (Zuc). Zuc-dependent cleavage typically occurs immediately upstream of a uridine. However, Zuc lacks sequence preference in vitro, pointing towards additional unknown specificity factors. Here, we examine murine piRNAs and reveal a strong and specific enrichment of three sequences (UAA, UAG, UGA)—corresponding to stop codons—at piRNA 5’ ends. Stop codon sequences are also enriched immediately after piRNA processing intermediates, reflecting their Zuc-dependent tail-to-head arrangement. Further analyses reveal that a Zuc in vivo cleavage preference at four sequences (UAA, UAG, UGA, UAC) promotes 5’ end stop codons. This observation is conserved across mammals and possibly further. Our work provides new insights into Zuc-dependent cleavage and may point to a previously unrecognized connection between piRNA biogenesis and the translational machinery.

Delineating Characteristic Sequence and Structural Features of Precursor and Mature Piwi-interacting RNAs of Epithelial Ovarian Cancer

Background: Piwi-interacting RNAs (piRNAs) are an amazing class of small noncoding RNAs (sncRNAs) known for its promising role in germline and somatic cells. Myriad functional studies have been performed to unveil the true potential of this class of ncRNAs; however, global features encoded in their sequence and structure have not been explored. Objectives: We aim to identify the sequence and structural level characteristic features of piRNAs of normal ovary (NO), and two subtypes of epithelial ovarian cancer (EOCa)- endometrioid (ENOCa), and serous ovarian cancer (SOCa) that we had reported earlier and their precursors. Methods: We have performed sequence analysis of mature piRNAs and their upstream/downstream regions as well as structural analysis of precursor piRNAs (pre-piRNAs) by examining their minimal folding energy (MFE), adjusted minimal folding energy (AMFE) and minimal folding energy index (MFEI) etc using in silico approaches. Results: We observed enrichment of U at first position and G at several other positions of mature piRNAs, which might be associated with the processing of mature piRNAs similar to what is seen in the case of miRNAs and strong target binding, respectively. In addition, we found the richness of AU in and around 20 nts upstream and downstream of precursor piRNAs (pre-piRNAs). This characteristic feature of pre-piRNAs possibly contributes to lower MFE compared to random sequences and make its secondary structure less stable which decides biogenesis of piRNAs. We also noticed that MFE, AMFE and MFEI of pre-piRNAs are comparatively less than pre-miRNAs of metazoans, plants and viruses reported in other studies, which clearly discriminate pre-piRNAs from other RNA sequences including pre-miRNAs of other organisms. Conclusion: In summary, the present study reveals key characteristic features encoded within and around mature piRNAs as well as pre-piRNAs of NO and EOCa samples that distinguish piRNAs from miRNAs and other random RNA sequences. These findings might act as a cornerstone for better understanding of biogenesis and function of piRNAs as well as will aid in easier identification of new piRNAs from unknown stretch of sequences using the characteristic features.

The RNA polymerase II subunit RPB‐9 recruits the integrator complex to terminate Caenorhabditis elegans piRNA transcription

PIWI‐interacting RNAs (piRNAs) are genome‐encoded small RNAs that regulate germ cell development and maintain germline integrity in many animals. Mature piRNAs engage Piwi Argonaute proteins to silence complementary transcripts, including transposable elements and endogenous genes. piRNA biogenesis mechanisms are diverse and remain poorly understood. Here, we identify the RNA polymerase II (RNA Pol II) core subunit RPB‐9 as required for piRNA‐mediated silencing in the nematode Caenorhabditis elegans. We show that rpb‐9 initiates heritable piRNA‐mediated gene silencing at two DNA transposon families and at a subset of somatic genes in the germline. We provide genetic and biochemical evidence that RPB‐9 is required for piRNA biogenesis by recruiting the Integrator complex at piRNA genes, hence promoting transcriptional termination. We conclude that, as a part of its rapid evolution, the piRNA pathway has co‐opted an ancient machinery for high‐fidelity transcription.

Zucchini consensus motifs determine the mechanism of pre-piRNA production

PIWI-interacting RNAs (piRNAs) of between approximately 24 and 31 nucleotides in length guide PIWI proteins to silence transposons in animal gonads, thereby ensuring fertility1. In the biogenesis of piRNAs, PIWI proteins are first loaded with 5'-monophosphorylated RNA fragments called pre-pre-piRNAs, which then undergo endonucleolytic cleavage to produce pre-piRNAs1,2. Subsequently, the 3'-ends of pre-piRNAs are trimmed by the exonuclease Trimmer (PNLDC1 in mouse)3-6 and 2'-O-methylated by the methyltransferase Hen1 (HENMT1 in mouse)7-9, generating mature piRNAs. It is assumed that the endonuclease Zucchini (MitoPLD in mouse) is a major enzyme catalysing the cleavage of pre-pre-piRNAs into pre-piRNAs10-13. However, direct evidence for this model is lacking, and how pre-piRNAs are generated remains unclear. Here, to analyse pre-piRNA production, we established a Trimmer-knockout silkworm cell line and derived a cell-free system that faithfully recapitulates Zucchini-mediated cleavage of PIWI-loaded pre-pre-piRNAs. We found that pre-piRNAs are generated by parallel Zucchini-dependent and -independent mechanisms. Cleavage by Zucchini occurs at previously unrecognized consensus motifs on pre-pre-piRNAs, requires the RNA helicase Armitage, and is accompanied by 2'-O-methylation of pre-piRNAs. By contrast, slicing of pre-pre-piRNAs with weak Zucchini motifs is achieved by downstream complementary piRNAs, producing pre-piRNAs without 2'-O-methylation. Regardless of the endonucleolytic mechanism, pre-piRNAs are matured by Trimmer and Hen1. Our findings highlight multiplexed processing of piRNA precursors that supports robust and flexible piRNA biogenesis.

The piRNA pathway in Drosophila ovarian germ and somatic cells.

RNA silencing refers to gene silencing pathways mediated by small non-coding RNAs, including microRNAs. Piwi-interacting RNAs (piRNAs) constitute the largest class of small non-coding RNAs in animal gonads, which repress transposons to protect the germline genome from the selfish invasion of transposons. Deterioration of the system causes DNA damage, leading to severe defects in gametogenesis and infertility. Studies using Drosophila ovaries show that piRNAs originate from specific genomic loci, termed piRNA clusters, and that in piRNA biogenesis, cluster transcripts are processed into mature piRNAs via three distinct pathways: initiator or responder for ping-pong piRNAs and trailing for phased piRNAs. piRNAs then assemble with PIWI members of the Argonaute family of proteins to form piRNA-induced RNA silencing complexes (piRISCs), the core engine of the piRNA-mediated silencing pathway. Upon piRISC assembly, the PIWI member, Piwi, is translocated to the nucleus and represses transposons co-transcriptionally by inducing local heterochromatin formation at target transposon loci.

Somatic expression of piRNA and associated machinery in the mouse identifies short, tissue-specific piRNA

ABSTRACTPiwi-interacting RNAs (piRNAs) are small non-coding RNAs that associate with PIWI proteins for transposon silencing via DNA methylation and are highly expressed and extensively studied in the germline. Mature germline piRNAs typically consist of 24–32 nucleotides, with a strong preference for a 5ʹ uridine signature, an adenosine signature at position 10, and a 2ʹ-O-methylation signature at the 3ʹ end. piRNA presence in somatic tissues, however, is not well characterized and requires further systematic evaluation. In the current study, we identified piRNAs and associated machinery from mouse somatic tissues representing the three germ layers. piRNA specificity was improved by combining small RNA size selection, sodium periodate treatment enrichment for piRNA over other small RNA, and small RNA next-generation sequencing. We identify PIWIL1, PIWIL2, and PIWIL4 expression in brain, liver, kidney, and heart. Of note, somatic piRNAs are shorter in length and tissue-specific, with increased occurrence of unique piRNAs in hippocampus and liver, compared to the germline. Hippocampus contains 5,494 piRNA-like peaks, the highest expression among all tested somatic tissues, followed by cortex (1,963), kidney (580), and liver (406). The study identifies 26 piRNA sequence species and 40 piRNA locations exclusive to all examined somatic tissues. Although piRNA expression has long been considered exclusive to the germline, our results support that piRNAs are expressed in several somatic tissues that may influence piRNA functions in the soma. Once confirmed, the PIWI/piRNA system may serve as a potential tool for future research in epigenome editing to improve human health by manipulating DNA methylation.

Diverse RNA interference strategies in early-branching metazoans

BackgroundMicro RNAs (miRNAs) and piwi interacting RNAs (piRNAs), along with the more ancient eukaryotic endogenous small interfering RNAs (endo-siRNAs) constitute the principal components of the RNA interference (RNAi) repertoire of most animals. RNAi in non-bilaterians – sponges, ctenophores, placozoans and cnidarians - appears to be more diverse than that of bilaterians, and includes structurally variable miRNAs in sponges, an enormous number of piRNAs in cnidarians and the absence of miRNAs in ctenophores and placozoans.ResultsHere we identify thousands of endo-siRNAs and piRNAs from the sponge Amphimedon queenslandica, the ctenophore Mnemiopsis leidyi and the cnidarian Nematostella vectensis using a computational approach that clusters mapped small RNA sequences and annotates each cluster based on the read length and relative abundance of the constituent reads. This approach was validated on 11 small RNA libraries in Drosophila melanogaster, demonstrating the successful annotation of RNAi-associated loci with properties consistent with previous reports. In the non-bilaterians we uncover seven new miRNAs from Amphimedon and four from Nematostella as well as sub-populations of candidate cis-natural antisense transcript (cis-NAT) endo-siRNAs. We confirmed the absence of miRNAs in Mnemiopsis but detected an abundance of endo-siRNAs in this ctenophore. Analysis of putative piRNA structure suggests that conserved localised secondary structures in primary transcripts may be important for the production of mature piRNAs in Amphimedon and Nematostella, as is also the case for endo-siRNAs.ConclusionTogether, these findings suggest that the last common ancestor of extant animals did not have the entrained RNAi system that typifies bilaterians. Instead it appears that bilaterians, cnidarians, ctenophores and sponges express unique repertoires and combinations of miRNAs, piRNAs and endo-siRNAs.

Editorial on “PNLDC1 is essential for piRNA 3' end trimming and transposon silencing during spermatogenesis in mice”

Piwi-interacting RNA (piRNA) are a class of small non-coding RNAs expressed in the germline. These piRNAs are processed to a mature piRNA by endonucleases and exonucleases that are bound to PIWI family proteins (1-3). In mouse, three Piwi proteins have been shown to be critical for spermatogenesis, namely Miwi, Mili and Miwi2 (4-6). Mili and Miwi show different temporal expression during testes development.

Hierarchical roles of mitochondrial Papi and Zucchini in Bombyx germline piRNA biogenesis.

PIWI-interacting RNAs (piRNAs) are small regulatory RNAs that bind to PIWI proteins to control transposons and maintain genome integrity in animal germ lines. piRNA 3' end formation in the silkworm Bombyx mori has been shown to be mediated by the 3'-to-5' exonuclease Trimmer (Trim; known as PNLDC1 in mammals), and piRNA intermediates are bound with PIWI anchored onto mitochondrial Tudor domain protein Papi. However, it remains unclear whether the Zucchini (Zuc) endonuclease and Nibbler (Nbr) 3'-to-5' exonuclease, both of which have pivotal roles in piRNA biogenesis in Drosophila, are required for piRNA processing in other species. Here we show that the loss of Zuc in Bombyx had no effect on the levels of Trim and Nbr, but resulted in the aberrant accumulation of piRNA intermediates within the Papi complex, and that these were processed to form mature piRNAs by recombinant Zuc. Papi exerted its RNA-binding activity only when bound with PIWI and phosphorylated, suggesting that complex assembly involves a hierarchical process. Both the 5' and 3' ends of piRNA intermediates within the Papi complex showed hallmarks of PIWI 'slicer' activity, yet no phasing pattern was observed in mature piRNAs. The loss of Zuc did not affect the 5'- and 3'-end formation of the intermediates, strongly supporting the idea that the 5' end of Bombyx piRNA is formed by PIWI slicer activity, but independently of Zuc, whereas the 3' end is formed by the Zuc endonuclease. The Bombyx piRNA biogenesis machinery is simpler than that of Drosophila, because Bombyx has no transcriptional silencing machinery that relies on phased piRNAs.

Zucchini-dependent piRNA processing is triggered by recruitment to the cytoplasmic processing machinery.

The piRNA pathway represses transposable elements in the gonads and thereby plays a vital role in protecting the integrity of germline genomes of animals. Mature piRNAs are processed from longer transcripts, piRNA precursors (pre-piRNAs). In Drosophila, processing of pre-piRNAs is initiated by piRNA-guided Slicer cleavage or the endonuclease Zucchini (Zuc). As Zuc does not have any sequence or structure preferences in vitro, it is not known how piRNA precursors are selected and channeled into the Zuc-dependent processing pathway. We show that a heterologous RNA that lacks complementary piRNAs is processed into piRNAs upon recruitment of several piRNA pathway factors. This processing requires Zuc and the helicase Armitage (Armi). Aubergine (Aub), Argonaute 3 (Ago3), and components of the nuclear RDC complex, which are required for normal piRNA biogenesis in germ cells, are dispensable. Our approach allows discrimination of proteins involved in the transcription and export of piRNA precursors from components required for the cytoplasmic processing steps. piRNA processing correlates with localization of the substrate RNA to nuage, a distinct membraneless cytoplasmic compartment, which surrounds the nucleus of germ cells, suggesting that sequestration of RNA to this subcellular compartment is both necessary and sufficient for selecting piRNA biogenesis substrates.

Increasing cell density globally enhances the biogenesis of Piwi-interacting RNAs in Bombyx mori germ cells

Piwi proteins and their bound Piwi-interacting RNAs (piRNAs) are predominantly expressed in the germline and play crucial roles in germline development by silencing transposons and other targets. Bombyx mori BmN4 cells are culturable germ cells that equip the piRNA pathway. Because of the scarcity of piRNA-expressing culturable cells, BmN4 cells are being utilized for the analyses of piRNA biogenesis. We here report that the piRNA biogenesis in BmN4 cells is regulated by cell density. As cell density increased, the abundance of Piwi proteins and piRNA biogenesis factors was commonly upregulated, resulting in an increased number of perinuclear nuage-like granules where Piwi proteins localize. Along with these phenomena, the abundance of mature piRNAs also globally increased, whereas levels of long piRNA precursor and transposons decreased, suggesting that increasing cell density promotes piRNA biogenesis pathway and that the resultant accumulation of mature piRNAs is functionally significant for transposon silencing. Our study reveals a previously uncharacterized link between cell density and piRNA biogenesis, designates cell density as a critical variable in piRNA studies using BmN4 cell system, and suggests the alteration of cell density as a useful tool to monitor piRNA biogenesis and function.

Roles of MIWI, MILI and PLD6 in small RNA regulation in mouse growing oocytes.

The mouse PIWI-interacting RNA (piRNA) pathway produces a class of 26–30-nucleotide (nt) small RNAs and is essential for spermatogenesis and retrotransposon repression. In oocytes, however, its regulation and function are poorly understood. In the present study, we investigated the consequences of loss of piRNA-pathway components in growing oocytes. When MILI (or PIWIL2), a PIWI family member, was depleted by gene knockout, almost all piRNAs disappeared. This severe loss of piRNA was accompanied by an increase in transcripts derived from specific retrotransposons, especially IAPs. MIWI (or PIWIL1) depletion had a smaller effect. In oocytes lacking PLD6 (or ZUCCHINI or MITOPLD), a mitochondrial nuclease/phospholipase involved in piRNA biogenesis in male germ cells, the piRNA level was decreased to 50% compared to wild-type, a phenotype much milder than that in males. Since PLD6 is essential for the creation of the 5΄ ends of primary piRNAs in males, the presence of mature piRNA in PLD6-depleted oocytes suggests the presence of compensating enzymes. Furthermore, we identified novel 21–23-nt small RNAs, termed spiRNAs, possessing a 10-nt complementarity with piRNAs, which were produced dependent on MILI and independent of DICER. Our study revealed the differences in the biogenesis and function of the piRNA pathway between sexes.

Genetic and mechanistic diversity of piRNA 3'-end formation.

Small regulatory RNAs guide Argonaute (Ago) proteins in a sequence-specific manner to their targets and thereby play important roles in eukaryotic gene silencing1. Of the three small RNA classes, microRNAs and siRNAs are processed from double-stranded precursors into defined 21- to 23-mers by Dicer, an endoribonuclease with intrinsic ruler function. piRNAs—the 22-30 nt long guides for PIWI-clade Ago proteins that silence transposons in animal gonads—are generated Dicer-independently from single-stranded precursors2,3. piRNA 5' ends are defined either by Zucchini, a mitochondria-anchored endonuclease4,5, or by piRNA-guided target cleavage6,7. Formation of piRNA 3' ends is poorly understood. Here, we find that two genetically and mechanistically distinct pathways generate piRNA 3' ends in Drosophila. The initiating nucleases are either Zucchini or the PIWI-clade proteins Aubergine (Aub)/Ago3. While Zucchini-mediated cleavages directly define mature piRNA 3' ends8,9, Aub/Ago3-mediated cleavages liberate pre-piRNAs that require extensive resection by the 3'-to-5' exoribonuclease Nibbler/Mut-710–13. The relative activity of these two pathways dictates the extent to which piRNAs are fueled into cytoplasmic or nuclear PIWI-clade proteins and thereby sets the balance between post-transcriptional and transcriptional silencing. Strikingly, loss of both Zucchini and Nibbler reveals a minimal, Argonaute-driven small RNA biogenesis pathway where piRNA 5' and 3' ends are directly produced by closely spaced Aub/Ago3-mediated cleavage events. Our data establish a coherent blueprint for piRNA biogenesis, and set the stage for the mechanistic dissection of the processes that govern piRNA 3' end formation.

A novel model of THO/TREX loading onto target RNAs in metazoan gene expression.

The THO/TREX complex consists of several conserved subunits and is required for mRNA export. In metazoans, THO/TREX binds a subset of mRNAs during RNA splicing, and facilitates their nuclear export. How THO/TREX selects RNA targets is, however, incompletely understood. In our recent study, we reported that THO is loaded onto Piwi-interacting RNA (piRNA) precursor transcripts independent of splicing, and facilitates convergent transcription in Drosophila ovary. The precursors are later processed into mature piRNAs, small noncoding RNAs that silence transposable elements (TEs). We observed that piRNAs originating from dual-strand clusters, where precursors are transcribed from both strands, were specifically affected by THO mutation. Analysis of THO-bound RNAs showed enrichment of dual-strand cluster transcripts. Interestingly, THO loading onto piRNA precursors was dependent on Cutoff (Cuff), which comprises the Rhino-Deadlock-Cutoff (RDC) complex that is recruited to dual-strand clusters by recognizing H3K9me3 and licenses convergent transcription from he cluster. We also found that THO mutation affected transcription from dual-strand clusters. Therefore, we concluded that THO/TREX is recruited to dual-strand piRNA clusters, independent of splicing events, via multi-protein interactions with chromatin structure. Then, it facilitates transcription likely by suppressing premature termination to ensure adequate expression of piRNA precursors. [BMB Reports 2016; 49(7): 355-356]

Identification and Functional Analysis of the Pre-piRNA 3′ Trimmer in Silkworms

PIWI-interacting RNAs (piRNAs) play a crucial role in transposon silencing in animal germ cells. In piRNA biogenesis, single-stranded piRNA intermediates are loaded into PIWI-clade proteins and cleaved by Zucchini/MitoPLD, yielding precursor piRNAs (pre-piRNAs). Pre-piRNAs that are longer than the mature piRNA length are then trimmed at their 3' ends. Although recent studies implicated the Tudor domain protein Papi/Tdrkh in pre-piRNA trimming, the identity of Trimmer and its relationship with Papi/Tdrkh remain unknown. Here, we identified PNLDC1, an uncharacterized 3'-5' exonuclease, as Trimmer in silkworms. Trimmer is enriched in the mitochondrial fraction and binds to Papi/Tdrkh. Depletion of Trimmer and Papi/Tdrkh additively inhibits trimming, causing accumulation of ∼35-40-nt pre-piRNAs that are impaired for target cleavage and prone to degradation. Our results highlight the cooperative action of Trimmer and Papi/Tdrkh in piRNA maturation.

PIWI-Interacting RNA: Its Biogenesis and Functions.

PIWI-interacting RNAs (piRNAs) are a class of small RNAs that are 24-31 nucleotides in length. They associate with PIWI proteins, which constitute a germline-specific subclade of the Argonaute family, to form effector complexes known as piRNA-induced silencing complexes, which repress transposons via transcriptional or posttranscriptional mechanisms and maintain germline genome integrity. In addition to having a role in transposon silencing, piRNAs in diverse organisms function in the regulation of cellular genes. In some cases, piRNAs have shown transgenerational inheritance to pass on the memory of "self" and "nonself," suggesting a contribution to various cellular processes over generations. Many piRNA factors have been identified; however, both the molecular mechanisms leading to the production of mature piRNAs and the effector phases of gene silencing are still enigmatic. Here, we summarize the current state of our knowledge on the biogenesis of piRNA, its biological functions, and the underlying mechanisms.

Prediction of mature microRNA and piwi-interacting RNA without a genome reference or precursors.

The discovery of novel microRNA (miRNA) and piwi-interacting RNA (piRNA) is an important task for the understanding of many biological processes. Most of the available miRNA and piRNA identification methods are dependent on the availability of the organism’s genome sequence and the quality of its annotation. Therefore, an efficient prediction method based solely on the short RNA reads and requiring no genomic information is highly desirable. In this study, we propose an approach that relies primarily on the nucleotide composition of the read and does not require reference genomes of related species for prediction. Using an empirical Bayesian kernel method and the error correcting output codes framework, compact models suitable for large-scale analyses are built on databases of known mature miRNAs and piRNAs. We found that the usage of an L1-based Gaussian kernel can double the true positive rate compared to the standard L2-based Gaussian kernel. Our approach can increase the true positive rate by at most 60% compared to the existing piRNA predictor based on the analysis of a hold-out test set. Using experimental data, we also show that our approach can detect about an order of magnitude or more known miRNAs than the mature miRNA predictor, miRPlex.

Guardian small RNAs and sex determination

The W chromosome of the silkworm Bombyx mori has been known to determine femaleness for more than 80 years. However, the feminizing gene has not been molecularly identified, because the B. mori W chromosome is almost fully occupied by a large number of transposable elements. The W chromosome-derived feminizing factor of B. mori was recently shown to be a female-specific PIWI-interacting RNA (piRNA). piRNAs are small RNAs that potentially repress invading “non-self” elements (e.g., transposons and virus-like elements) by associating with PIWI proteins. Our results revealed that female-specific piRNA precursors, which we named Fem, are transcribed from the sex-determining region of the W chromosome at the early embryonic stage and are processed into a single mature piRNA (Fem piRNA). Fem piRNA forms a complex with Siwi (silkworm Piwi), which cleaves a protein-coding mRNA transcribed from the Z chromosome. RNA interference of this Z-linked gene, which we named Masc, revealed that this gene encodes a protein required for masculinization and dosage compensation. Fem and Masc both participate in the ping-pong cycle of the piRNA amplification loop by associating with the 2 B. mori PIWI proteins Siwi and BmAgo3 (silkworm Ago3), respectively, indicating that the piRNA-mediated interaction between the 2 sex chromosomes is the primary signal for the B. mori sex determination cascade. Fem is a non-transposable repetitive sequence on the W chromosome, whereas Masc is a single-copy protein-coding gene. It is of great interest how the piRNA system recognizes “self ”Masc mRNA as “non-self” RNA.