RECURRENT AND RECENT SELECTIVE SWEEPS IN THE piRNA PATHWAY

Uncontrolled transposable element (TE) insertions and excisions can cause chromosome breaks and mutations with dramatic deleterious effects. The PIWI interacting RNA (piRNA) pathway functions as an adaptive TE silencing system during germline development. Several essential piRNA pathway proteins appear to be rapidly evolving, suggesting that TEs and the silencing machinery may be engaged in a classical “evolutionary arms race.” Using a variety of molecular evolutionary and population genetic approaches, we find that the piRNA pathway genes rhino, krimper, and aubergine show patterns suggestive of extensive recurrent positive selection across Drosophila species. We speculate that selection on these proteins reflects crucial roles in silencing unfamiliar elements during vertical and horizontal transmission of TEs into naïve populations and species, respectively.     

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.     

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.     

Spindle-E is essential for gametogenesis in the silkworm,Bombyx mori

As a defense mechanism against transposable elements, the PIWI-interacting RNA (piRNA) pathway maintains genomic integrity and ensures proper gametogenesis in gonads. Numerous factors are orchestrated to ensure normal operation of the piRNA pathway. Spindle-E (Spn-E) gene was one of the first genes shown to participate in the piRNA pathway. In this study, we performed functional analysis of Spn-E in the model lepidopteran insect, Bombyx mori. Unlike the germline-specific expression pattern observed in Drosophila and mouse, BmSpn-E was ubiquitously expressed in all tissues tested, and it was highly expressed in gonads. Immunofluorescent staining showed that BmSpn-E was localized in both germ cells and somatic cells in ovary and was expressed in spermatocytes in testis. We used a binary transgenic CRISPR/Cas9 system to construct BmSpn-E mutants. Loss of BmSpn-E expression caused derepression of transposons in gonads. We also found that mutant gonads were much smaller than wild-type gonads and that the number of germ cells was considerably lower in mutant gonads. Quantitative real-time PCR analysis and TUNEL staining revealed that apoptosis was greatly enhanced in mutant gonads. Further, we found that the BmSpn-E mutation impacted gonadal development and gametogenesis at the early larval stage. In summary, our data provided the first evidence that BmSpn-E plays vital roles in gonadal development and gametogenesis in B. mori.     

Temporal expression and steroidal regulation of piRNA pathway genes (mael, piwi, vasa) during Silurana (Xenopus) tropicalis embryogenesis and early larval development

It has been extensively documented that exposure of amphibians and teleost fish to exogenous steroid hormones like estrogen, androgen, xenoestrogen or steroid biosynthesis inhibitors can impair their gonadal development or induce sex reversal against genotypic sex. However, the molecular pathways underlying sexual development and the effects of sex steroids or other exogenous hormones in these aquatic vertebrates remain elusive. Recently, a germ plasm-associated piRNA (piwi-interacting RNA) pathway has been shown to be a determinant in the development of animal gonadal germline cells. In the current study, we examined whether this piRNA pathway is involved in the regulation of sex steroid hormones in gonadal development. We firstly established developmental expression patterns of three key piRNA pathway genes (mael, piwi and vasa), during Silurana (Xenopus) tropicalis embryogenesis and early larval development. All three genes exhibit high expression at early developmental stages and have significantly decreased expression thereafter, indicating a very active involvement of piRNA pathway at the beginning of embryogenesis. We further examined gene expression changes of those genes in frog larvae exposed to two sex steroid biosynthesis inhibitors, fadrozole and finasteride, both of which are known to result in male-biased or female-biased phenotypes, respectively. We found that fadrozole and finasteride exposures increased the expression of piRNA pathway genes such as mael and vasa at the larval stage when the expression of piRNA pathway genes is programmed to be very low. Therefore, our results indicate that the piRNA pathway is likely a common pathway by which different sex steroid hormones regulate gonadal sex differentiation.     

Distinct functions for the Drosophila piRNA pathway in genome maintenance and telomere protection

Transposons and other selfish DNA elements can be found in all phyla, and mobilization of these elements can compromise genome integrity. The piRNA (PIWI-interacting RNA) pathway silences transposons in the germline, but it is unclear if this pathway has additional functions during development. Here we show that mutations in the Drosophila piRNA pathway genes, armi, aub, ago3, and rhi, lead to extensive fragmentation of the zygotic genome during the cleavage stage of embryonic divisions. Additionally, aub and armi show defects in telomere resolution during meiosis and the cleavage divisions; and mutations in lig-IV, which disrupt non-homologous end joining, suppress these fusions. By contrast, lig-IV mutations enhance chromosome fragmentation. Chromatin immunoprecipitation studies show that aub and armi mutations disrupt telomere binding of HOAP, which is a component of the telomere protection complex, and reduce expression of a subpopulation of 19- to 22-nt telomere-specific piRNAs. Mutations in rhi and ago3, by contrast, do not block HOAP binding or production of these piRNAs. These findings uncover genetically separable functions for the Drosophila piRNA pathway. The aub, armi, rhi, and ago3 genes silence transposons and maintain chromosome integrity during cleavage-stage embryonic divisions. However, the aub and armi genes have an additional function in assembly of the telomere protection complex.     

Lint‐O cooperates with L(3)mbt in target gene suppression to maintain homeostasis in fly ovary and brain

Loss‐of‐function mutations in Drosophila lethal(3)malignant brain tumor [l(3)mbt] cause ectopic expression of germline genes and brain tumors. Loss of L(3)mbt function in ovarian somatic cells (OSCs) aberrantly activates germ‐specific piRNA amplification and leads to infertility. However, the underlying mechanism remains unclear. Here, ChIP‐seq for L(3)mbt in cultured OSCs and RNA‐seq before and after L(3)mbt depletion shows that L(3)mbt genomic binding is not necessarily linked to gene regulation and that L(3)mbt controls piRNA pathway genes in multiple ways. Lack of known L(3)mbt co‐repressors, such as Lint‐1, has little effect on the levels of piRNA amplifiers. Identification of L(3)mbt interactors in OSCs and subsequent analysis reveals CG2662 as a novel co‐regulator of L(3)mbt, termed “L(3)mbt interactor in OSCs” (Lint‐O). Most of the L(3)mbt‐bound piRNA amplifier genes are also bound by Lint‐O in a similar fashion. Loss of Lint‐O impacts the levels of piRNA amplifiers, similar to the lack of L(3)mbt. The lint‐O‐deficient flies exhibit female sterility and tumorous brains. Thus, L(3)mbt and its novel co‐suppressor Lint‐O cooperate in suppressing target genes to maintain homeostasis in the ovary and brain.     

The piRNA response to BmNPV infection in the silkworm fat body and midgut

Bombyx mori nucleopolyhedrovirus (BmNPV) is a DNA virus that causes huge losses to the silkworm industry but the piRNA responses during BmNPV infection in the silkworm remain uninvestigated. Here, silkworm piRNA profiles of uninfected and BmNPV-infected fat body and midgut were determined by high-through sequencing in the early stages of BmNPV infection. A total of 2675 and 3396 genome-derived piRNAs were identified from fat body and midgut, respectively. These genome-derived piRNAs mainly originated from unannotated instead of transposon regions in the silkworm genome. In total, 572 piRNAs were associated with 280 putative target genes in fat body and 805 piRNAs with 380 target genes in midgut. Compared to uninfected tissues, 322 and 129 piRNAs were significantly upregulated in BmNPV-infected fat body and midgut, respectively. In addition, 276 and 117 piRNAs were significantly downregulated. Moreover, differentially expressed (DE) piRNAs during BmNPV infection differed significantly between fat body and midgut. Putative DE piRNA–targeted genes were associated with “response to stimulus” and “environmental information processing” in fat body after infection with BmNPV, which may indicate an active piRNA response to BmNPV infection in fat body. This study may lay the foundation for future research of the potential roles of the piRNA pathway and specific piRNAs in BmNPV pathogenesis.     

Interplay of transposon-silencing genes in the germline of <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

Complexes of Piwi family proteins with short piRNAs (Piwi-interacting RNAs) are responsible for silencing transposable elements in animal reproductive organs. In Drosophila melanogaster, three proteins (Piwi, Aub, and Ago3) are members of the Piwi family. Piwi is the nuclear protein of somatic and germinal ovarian cells, whereas Aub and Ago3 are cytoplasmic proteins involved in piRNA amplification in perinuclear granules that constitute special organelles of germinal cells called nuage. Mutations in genes of the piRNA system are known to cause derepression of several transposable elements. In this study, we compared quantitatively changes in expression of a larger number of elements in the case of mutations in the piwi gene, genes aub, mael, and spn-E, which encode proteins of nuage granules, and armi gene coding an RNA helicase, the lack of which does not interfere with cytoplasmic piRNA amplification but disturbs nuclear localization of Piwi protein. We found that the genes piwi, armi, aub, spn-E, and mael interact to induce silencing of some retrotransposons (HMS-Beagle, Gate and HeT-A); the same genes, except piwi, are involved in repression of I and G elements. We propose that Armi is involved in control of not only nuclear Piwi localization. Our data suggest the relation of nuage proteins Aub, Spn-E, and Mael to Piwi-mediated silencing of retrotransposons Gate and HMS-Beagle in the nucleus. In general, our results corroborate the idea of genome stabilization by means of various silencing strategies specific to different transposable elements. At the same time, our data suggest the existence of yet unknown mechanisms of interplay between nuclear and cytoplasmic components of the piRNA machinery in germinal cells.