Distinct sets of PIWI proteins produce arbovirus and transposon-derived piRNAs in Aedes aegypti mosquito cells

The PIWI-interacting RNA (piRNA) pathway is essential for transposon silencing in many model organisms. Its remarkable efficiency relies on a sophisticated amplification mechanism known as the ping-pong loop. In Alphavirus-infected Aedes mosquitoes, piRNAs with sequence features that suggest ping-pong-dependent biogenesis are produced from viral RNA. The PIWI family in Aedes mosquitoes is expanded when compared to other model organisms, raising the possibility that individual PIWI proteins have functionally diversified in these insects. Here, we show that Piwi5 and Ago3, but none of the other PIWI family members, are essential for piRNA biogenesis from Sindbis virus RNA in infected Aedes aegypti cells. In contrast, the production of piRNAs from transposons relies on a more versatile set of PIWI proteins, some of which do not contribute to viral piRNA biogenesis. These results indicate that functional specialization allows distinct mosquito PIWI proteins to process RNA from different endogenous and exogenous sources.     

An in vivo RNAi assay identifies major genetic and cellular requirements for primary piRNA biogenesis in Drosophila

In Drosophila, PIWI proteins and bound PIWI‐interacting RNAs (piRNAs) form the core of a small RNA‐mediated defense system against selfish genetic elements. Within germline cells, piRNAs are processed from piRNA clusters and transposons to be loaded into Piwi/Aubergine/AGO3 and a subset of piRNAs undergoes target‐dependent amplification. In contrast, gonadal somatic support cells express only Piwi, lack signs of piRNA amplification and exhibit primary piRNA biogenesis from piRNA clusters. Neither piRNA processing/loading nor Piwi‐mediated target silencing is understood at the genetic, cellular or molecular level. We developed an in vivo RNAi assay for the somatic piRNA pathway and identified the RNA helicase Armitage, the Tudor domain containing RNA helicase Yb and the putative nuclease Zucchini as essential factors for primary piRNA biogenesis. Lack of any of these proteins leads to transposon de‐silencing, to a collapse in piRNA levels and to a failure in Piwi‐nuclear accumulation. We show that Armitage and Yb interact physically and co‐localize in cytoplasmic Yb bodies, which flank P bodies. Loss of Zucchini leads to an accumulation of Piwi and Armitage in Yb bodies, indicating that Yb bodies are sites of primary piRNA biogenesis.     

The Bombyx ovary-derived cell line endogenously expresses PIWI/PIWI-interacting RNA complexes

Genetic studies and large-scale sequencing experiments have revealed that the PIWI subfamily proteins and PIWI-interacting RNAs (piRNAs) play an important role in germ line development and transposon control. Biochemical studies in vitro have greatly contributed to the understanding of small interfering RNA (siRNA) and microRNA (miRNA) pathways. However, in vitro analyses of the piRNA pathway have been thus far quite challenging, because their expression is largely restricted to the germ line. Here we report that Bombyx mori ovary-derived cultured cell line, BmN4, endogenously expresses two PIWI subfamily proteins, silkworm Piwi (Siwi) and Ago3 (BmAgo3), and piRNAs associated with them. Siwi-bound piRNAs have a strong bias for uridine at their 5′ end and BmAgo3-bound piRNAs are enriched for adenine at position 10. In addition, Siwi preferentially binds antisense piRNAs, whereas BmAgo3 binds sense piRNAs. Moreover, we identified many pairs in which Siwi-bound antisense and BmAgo3-bound sense piRNAs are overlapped by precisely 10 nt at their 5′ ends. These signatures are known to be important for secondary piRNA biogenesis in other organisms. Taken together, BmN4 is a unique cell line in which both primary and secondary steps of piRNA biogenesis pathways are active. This cell line would provide useful tools for analysis of piRNA biogenesis and function.     

Hsp90 facilitates accurate loading of precursor piRNAs into PIWI proteins

PIWI-interacting RNAs (piRNAs) defend the genome against transposon activity in animal gonads. The Hsp90 chaperone machinery has been implicated in the piRNA pathway, but its exact role remains obscure. Here, we examined the effect of 17-N-allylamino-17-demethoxygeldanamycin (17-AAG), an Hsp90-specific inhibitor, on the piRNA pathway. In the silkworm ovary-derived BmN4 cells, 17-AAG treatment reduced the level of piRNAs and PIWI proteins. In vitro, the 5′-nucleotide preference upon precursor piRNA loading was compromised by 17-AAG, whereas 3′-end trimming and 2′-O-methylation were unaffected. Our data highlight a role of Hsp90 in accurate loading of precursor piRNAs into PIWI proteins.     

Mitochondrial protein BmPAPI modulates the length of mature piRNAs

PIWI proteins and their associated PIWI-interacting RNAs (piRNAs) protect genome integrity by silencing transposons in animal germlines. The molecular mechanisms and components responsible for piRNA biogenesis remain elusive. PIWI proteins contain conserved symmetrical dimethylarginines (sDMAs) that are specifically targeted by TUDOR domain-containing proteins. Here we report that the sDMAs of PIWI proteins play crucial roles in PIWI localization and piRNA biogenesis in Bombyx mori-derived BmN4 cells, which harbor fully functional piRNA biogenesis machinery. Moreover, RNAi screenings for Bombyx genes encoding TUDOR domain-containing proteins identified BmPAPI, a Bombyx homolog of Drosophila PAPI, as a factor modulating the length of mature piRNAs. BmPAPI specifically recognized sDMAs and interacted with PIWI proteins at the surface of the mitochondrial outer membrane. BmPAPI depletion resulted in 3′-terminal extensions of mature piRNAs without affecting the piRNA quantity. These results reveal the BmPAPI-involved piRNA precursor processing mechanism on mitochondrial outer membrane scaffolds.     

The MID-PIWI module of Piwi proteins specifies nucleotide- and strand-biases of piRNAs

Piwi-interacting RNAs (piRNAs) guide Piwi Argonautes to suppress transposon activity in animal gonads. Known piRNA populations are extremely complex, with millions of individual sequences present in a single organism. Despite this complexity, specific Piwi proteins incorporate piRNAs with distinct nucleotide- and transposon strand-biases (antisense or sense) of unknown origin. Here, we examined the contribution of structural domains in Piwi proteins toward defining these biases. We report the first crystal structure of the MID domain from a Piwi Argonaute and use docking experiments to show its ability to specify recognition of 5′ uridine (1U-bias) of piRNAs. Mutational analyses reveal the importance of 5′ end-recognition within the MID domain for piRNA biogenesis in vivo. Finally, domain-swapping experiments uncover an unexpected role for the MID-PIWI module of a Piwi protein in dictating the transposon strand-orientation of its bound piRNAs. Our work identifies structural features that allow distinguishing individual Piwi members during piRNA biogenesis.     

piRNA clusters as a main source of small RNAs in the animal germline

PIWI subfamily Argonaute proteins and small RNAs bound to them (PIWI interacting RNA, piRNA) control mobilization of transposable elements (TE) in the animal germline. piRNAs are generated by distinct genomic regions termed piRNA clusters. piRNA clusters are often extensive loci enriched in damaged fragments of TEs. New TE integration into piRNA clusters causes production of TE-specific piRNAs and repression of cognate sequences. piRNAs are thought to be generated from long single-stranded precursors encoded by piRNA clusters. Special chromatin structures might be essential to distinguish these genomic loci as a source for piRNAs. In this review, we present recent findings on the structural organization of piRNA clusters and piRNA biogenesis in Drosophila and other organisms, which are important for understanding a key epigenetic mechanism that provides defense against TE expansion.     

Paramutation in Drosophila Requires Both Nuclear and Cytoplasmic Actors of the piRNA Pathway and Induces Cis-spreading of piRNA Production

Transposable element activity is repressed in the germline in animals by PIWI-interacting RNAs (piRNAs), a class of small RNAs produced by genomic loci mostly composed of TE sequences. The mechanism of induction of piRNA production by these loci is still enigmatic. We have shown that, in Drosophila melanogaster, a cluster of tandemly repeated P-lacZ-white transgenes can be activated for piRNA production by maternal inheritance of a cytoplasm containing homologous piRNAs. This activated state is stably transmitted over generations and allows trans-silencing of a homologous transgenic target in the female germline. Such an epigenetic conversion displays the functional characteristics of a paramutation, i.e., a heritable epigenetic modification of one allele by the other. We report here that piRNA production and trans-silencing capacities of the paramutated cluster depend on the function of the rhino, cutoff, and zucchini genes involved in primary piRNA biogenesis in the germline, as well as on that of the aubergine gene implicated in the ping-pong piRNA amplification step. The 21-nt RNAs, which are produced by the paramutated cluster, in addition to 23- to 28-nt piRNAs are not necessary for paramutation to occur. Production of these 21-nt RNAs requires Dicer-2 but also all the piRNA genes tested. Moreover, cytoplasmic transmission of piRNAs homologous to only a subregion of the transgenic locus can generate a strong paramutated locus that produces piRNAs along the whole length of the transgenes. Finally, we observed that maternally inherited transgenic small RNAs can also impact transgene expression in the soma. In conclusion, paramutation involves both nuclear (Rhino, Cutoff) and cytoplasmic (Aubergine, Zucchini) actors of the piRNA pathway. In addition, since it is observed between nonfully homologous loci located on different chromosomes, paramutation may play a crucial role in epigenome shaping in Drosophila natural populations.     

Rapid Intraspecific Evolution of miRNA and siRNA Genes in the Mosquito Aedes aegypti

RNA silencing, or RNA interference (RNAi) in metazoans mediates development, reduces viral infection and limits transposon mobility. RNA silencing involves 21–30 nucleotide RNAs classified into microRNA (miRNA), exogenous and endogenous small interfering RNAs (siRNA), and Piwi-interacting RNA (piRNA). Knock-out, silencing and mutagenesis of genes in the exogenous siRNA (exo-siRNA) regulatory network demonstrate the importance of this RNAi pathway in antiviral immunity in Drosophila and mosquitoes. In Drosophila, genes encoding components for processing exo-siRNAs are among the fastest evolving 3% of all genes, suggesting that infection with pathogenic RNA viruses may drive diversifying selection in their host. In contrast, paralogous miRNA pathway genes do not evolve more rapidly than the genome average. Silencing of exo-siRNA pathway genes in mosquitoes orally infected with arboviruses leads to increased viral replication, but little is known about the comparative patterns of molecular evolution among the exo-siRNA and miRNA pathways genes in mosquitoes. We generated nearly complete sequences of all exons of major miRNA and siRNA pathway genes dicer-1 and dicer-2, argonaute-1 and argonaute-2, and r3d1 and r2d2 in 104 Aedes aegypti mosquitoes collected from six distinct geographic populations and analyzed their genetic diversity. The ratio of replacement to silent amino acid substitutions was 1.4 fold higher in dicer-2 than in dicer-1, 27.4 fold higher in argonaute-2 than in argonaute-1 and similar in r2d2 and r3d1. Positive selection was supported in 32% of non-synonymous sites in dicer-1, in 47% of sites in dicer-2, in 30% of sites in argonaute-1, in all sites in argonaute-2, in 22% of sites in r3d1 and in 55% of sites in r2d2. Unlike Drosophila, in Ae. aegypti, both exo-siRNA and miRNA pathway genes appear to be undergoing rapid, positive, diversifying selection. Furthermore, refractoriness of mosquitoes to infection with dengue virus was significantly positively correlated for nucleotide diversity indices in dicer-2.