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.     

Novel quality metrics allow identifying and generating high-quality assemblies of piRNA clusters

In most animals, it is thought that the proliferation of a transposable element (TE) is stopped when the TE jumps into a piRNA cluster. Despite this central importance, little is known about the composition and the evolutionary dynamics of piRNA clusters. This is largely because piRNA clusters are notoriously difficult to assemble as they are frequently composed of highly repetitive DNA. With long reads, we may finally be able to obtain reliable assemblies of piRNA clusters. Unfortunately, it is unclear how to generate and identify the best assemblies, as many assembly strategies exist and standard quality metrics are ignorant of TEs. To address these problems, we introduce several novel quality metrics that assess: (a) the fraction of completely assembled piRNA clusters, (b) the quality of the assembled clusters and (c) whether an assembly captures the overall TE landscape of an organisms (i.e. the abundance, the number of SNPs and internal deletions of all TE families). The requirements for computing these metrics vary, ranging from annotations of piRNA clusters to consensus sequences of TEs and genomic sequencing data. Using these novel metrics, we evaluate the effect of assembly algorithm, polishing, read length, coverage, residual polymorphisms and finally identify strategies that yield reliable assemblies of piRNA clusters. Based on an optimized approach, we provide assemblies for the two Drosophila melanogaster strains Canton-S and Pi2. About 80% of known piRNA clusters were assembled in both strains. Finally, we demonstrate the generality of our approach by extending our metrics to humans and Arabidopsis thaliana.     

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.     

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.     

Drosophila Interspecific Hybrids Phenocopy piRNA-Pathway Mutants

The Piwi-interacting RNA (piRNA) pathway defends the germline of animals from the deleterious activity of selfish transposable elements (TEs) through small-RNA mediated silencing. Adaptation to novel invasive TEs is proposed to occur by incorporating their sequences into the piRNA pool that females produce and deposit into their eggs, which then propagates immunity against specific TEs to future generations. In support of this model, the F1 offspring of crosses between strains of the same Drosophila species sometimes suffer from germline derepression of paternally inherited TE families, caused by a failure of the maternal strain to produce the piRNAs necessary for their regulation. However, many protein components of the Drosophila piRNA pathway exhibit signatures of positive selection, suggesting that they also contribute to the evolution of host genome defense. Here we investigate piRNA pathway function and TE regulation in the F1 hybrids of interspecific crosses between D. melanogaster and D. simulans and compare them with intraspecific control crosses of D. melanogaster. We confirm previous reports showing that intraspecific crosses are characterized by derepression of paternally inherited TE families that are rare or absent from the maternal genome and piRNA pool, consistent with the role of maternally deposited piRNAs in shaping TE silencing. In contrast to the intraspecific cross, we discover that interspecific hybrids are characterized by widespread derepression of both maternally and paternally inherited TE families. Furthermore, the pattern of derepression of TE families in interspecific hybrids cannot be attributed to their paucity or absence from the piRNA pool of the maternal species. Rather, we demonstrate that interspecific hybrids closely resemble piRNA effector-protein mutants in both TE misregulation and aberrant piRNA production. We suggest that TE derepression in interspecific hybrids largely reflects adaptive divergence of piRNA pathway genes rather than species-specific differences in TE-derived piRNAs.     

A comprehensive analysis of piRNAs from adult human testis and their relationship with genes and mobile elements

Background

Piwi-interacting RNAs (piRNAs) are a recently discovered class of small non-coding RNAs whose best-understood function is to repress mobile element (ME) activity in animal germline. To date, nearly all piRNA studies have been conducted in model organisms and little is known about piRNA diversity, target specificity and biological function in human.

Results

Here we performed high-throughput sequencing of piRNAs from three human adult testis samples. We found that more than 81% of the ~17 million putative piRNAs mapped to ~6,000 piRNA-producing genomic clusters using a relaxed definition of clusters. A set of human protein-coding genes produces a relatively large amount of putative piRNAs from their 3’UTRs, and are significantly enriched for certain biological processes, suggestive of non-random sampling by the piRNA biogenesis machinery. Up to 16% of putative piRNAs mapped to a few hundred annotated long non-coding RNA (lncRNA) genes, suggesting that some lncRNA genes can act as piRNA precursors. Among major ME families, young families of LTR and endogenous retroviruses have a greater association with putative piRNAs than other MEs. In addition, piRNAs preferentially mapped to specific regions in the consensus sequences of several ME (sub)families and some piRNA mapping peaks showed patterns consistent with the “ping-pong” cycle of piRNA targeting and amplification.

Conclusions

Overall our data provide a comprehensive analysis and improved annotation of human piRNAs in adult human testes and shed new light into the relationship of piRNAs with protein-coding genes, lncRNAs, and mobile genetic elements in human.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-545) contains supplementary material, which is available to authorized users.     

Endogenous siRNAs and piRNAs derived from transposable elements and genes in the malaria vector mosquito Anopheles gambiae

Background

The siRNA and piRNA pathways have been shown in insects to be essential for regulation of gene expression and defence against exogenous and endogenous genetic elements (viruses and transposable elements). The vast majority of endogenous small RNAs produced by the siRNA and piRNA pathways originate from repetitive or transposable elements (TE). In D. melanogaster, TE-derived endogenous siRNAs and piRNAs are involved in genome surveillance and maintenance of genome integrity. In the medically relevant malaria mosquito Anopheles gambiae TEs constitute 12-16% of the genome size. Genetic variations induced by TE activities are known to shape the genome landscape and to alter the fitness in An. gambiae.

Results

Here, using bioinformatics approaches we analyzed the small RNA data sets from 6 libraries formally reported in a previous study and examined the expression of the mixed germline/somatic siRNAs and piRNAs produced in adult An. gambiae females. We characterized a large population of TE-derived endogenous siRNAs and piRNAs, which constitutes 56-60% of the total siRNA and piRNA reads in the analysed libraries. Moreover, we identified a number of protein coding genes producing gene-specific siRNAs and piRNAs that were generally expressed at much lower levels than the TE-associated small RNAs. Detailed sequence analysis revealed that An. gambiae piRNAs were produced by both “ping-pong” dependent (TE-associated piRNAs) and independent mechanisms (genic piRNAs). Similarly to D. melanogaster, more than 90% of the detected piRNAs were produced from TE-associated clusters in An. gambiae. We also found that biotic stress as blood feeding and infection with Plasmodium parasite, the etiological agent of malaria, modulated the expression levels of the endogenous siRNAs and piRNAs in An. gambiae.

Conclusions

We identified a large and diverse set of the endogenously derived siRNAs and piRNAs that share common and distinct aspects of small RNA expression across insect species, and inferred their impact on TE and gene activity in An. gambiae. The TE-specific small RNAs produced by both the siRNA and piRNA pathways represent an important aspect of genome stability and genetic variation, which might have a strong impact on the evolution of the genome and vector competence in the malaria mosquitoes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1436-1) contains supplementary material, which is available to authorized users.     

Small RNA profiling and characterization of piRNA clusters in the adult testes of the common marmoset,a model primate

Small RNAs mediate gene silencing by binding Argonaute/Piwi proteins to regulate target RNAs. Here, we describe small RNA profiling of the adult testes of Callithrix jacchus, the common marmoset. The most abundant class of small RNAs in the adult testis was piRNAs, although 353 novel miRNAs but few endo-siRNAs were also identified. MARWI, a marmoset homolog of mouse MIWI and a very abundant PIWI in adult testes, associates with piRNAs that show characteristics of mouse pachytene piRNAs. As in other mammals, most marmoset piRNAs are derived from conserved clustered regions in the genome, which are annotated as intergenic regions. However, unlike in mice, marmoset piRNA clusters are also found on the X chromosome, suggesting escape from meiotic sex chromosome inactivation by the X-linked clusters. Some of the piRNA clusters identified contain antisense-orientated pseudogenes, suggesting the possibility that pseudogene-derived piRNAs may regulate parental functional protein-coding genes. More piRNAs map to transposable element (TE) subfamilies when they have copies in piRNA clusters. In addition, the strand bias observed for piRNAs mapped to each TE subfamily correlates with the polarity of copies inserted in clusters. These findings suggest that pachytene piRNA clusters determine the abundance and strand-bias of TE-derived piRNAs, may regulate protein-coding genes via pseudogene-derived piRNAs, and may even play roles in meiosis in the adult marmoset testis.     

Quantitatively Increased Somatic Transposition of Transposable Elements in Drosophila Strains Compromised for RNAi

In Drosophila melanogaster, small RNAs homologous to transposable elements (TEs) are of two types: piRNA (piwi-interacting RNA) with size 23-29nt and siRNA (small interfering RNA) with size 19-22nt. The siRNA pathway is suggested to silence TE activities in somatic tissues based on TE expression profiles, but direct evidence of transposition is lacking. Here we developed an efficient FISH (fluorescence in Situ hybridization) based method for polytene chromosomes from larval salivary glands to reveal new TE insertions. Analysis of the LTR-retrotransposon 297 and the non-LTR retroposon DOC shows that in the argonaut 2 (Ago2) and Dicer 2 (Dcr2) mutant strains, new transposition events are much more frequent than in heterozygous strains or wild type strains. The data demonstrate that the siRNA pathway represses TE transposition in somatic cells. Nevertheless, we found that loss of one functional copy of Ago2 or Dcr2 increases somatic transpositions of the elements at a lower level depending on the genetic background, suggesting a quantitative role for RNAi core components on mutation frequency.     

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.