Transposable element dynamics and PIWI regulation impacts lncRNA and gene expression diversity in Drosophila ovarian cell cultures

CRISPR-Cas9 Genome Editing and Rapid Selection of Cell Pools

The harnessing of the CRISPR-Cas9 system allows for quick and inexpensive genome editing in tissue culture models. Traditional CRISPR-Cas9 genome editing techniques rely on the ability of single progenitor cells to expand into new pools in a process known as clonal expansion. This is a significant technical challenge that is difficult to overcome for nontransformed cell culture models such as Drosophila ovarian somatic sheath cells (OSCs). OSCs are a unique ex vivo model for epigenetic regulation by PIWI-interacting RNAs (piRNAs) that establish restriction of mobile genetic elements in germ cells to protect genome integrity. Here, we provide a protocol to generate endogenously tagged proteins and gene knockouts without the need for clonal selection. We combine CRISPR-Cas genome editing and knockin of antibiotic selection markers to generate edited cell pools. At the example of Drosophila piwi in OSCs, we demonstrate a strategy that relies on the insertion of an artificial intron to accommodate a selection marker with minimal disturbance of the resulting mRNA. In brief, our donor cassette contains a peptide tag and an optimized intron that accommodates a selection marker driven by an independent promoter on the other genomic strand. The selection marker is transcribed as an independent mRNA, and the intron is efficiently removed from the mRNA encoding the endogenously tagged (endo-tagged) piwi gene. The endo-tagged Piwi protein is expressed at wild-type levels and appropriately localizes to the nucleus of OSCs. We also describe strategies for C-terminal tagging and generation of knockout alleles in OSCs and in human embryonic kidney cells, discuss different design strategies, and provide a plasmid toolkit (available at Addgene). Our protocol enables robust genome editing in OSCs for the first time and provides a simple and time-saving alternative for other cell culture systems. Published 2022. This article is a U.S. Government work and is in the public domain in the USA. Basic Protocol 1: Design and cloning of single-guide RNA plasmids Basic Protocol 2: Design and cloning of donor template plasmids for epitope tagging Alternate Protocol: Design and cloning of donor template plasmids for gene knockout Basic Protocol 3: Transfection and selection of edited cell pools.

Ongoing transposition in cell culture reveals the phylogeny of diverse Drosophila S2 sublines

Cultured cells are widely used in molecular biology despite poor understanding of how cell line genomes change in vitro over time. Previous work has shown that Drosophila cultured cells have a higher transposable element content than whole flies, but whether this increase in transposable element content resulted from an initial burst of transposition during cell line establishment or ongoing transposition in cell culture remains unclear. Here, we sequenced the genomes of 25 sublines of Drosophila S2 cells and show that transposable element insertions provide abundant markers for the phylogenetic reconstruction of diverse sublines in a model animal cell culture system. DNA copy number evolution across S2 sublines revealed dramatically different patterns of genome organization that support the overall evolutionary history reconstructed using transposable element insertions. Analysis of transposable element insertion site occupancy and ancestral states support a model of ongoing transposition dominated by episodic activity of a small number of retrotransposon families. Our work demonstrates that substantial genome evolution occurs during long-term Drosophila cell culture, which may impact the reproducibility of experiments that do not control for subline identity.

Transposable element landscapes in aging Drosophila

Genetic mechanisms that repress transposable elements (TEs) in young animals decline during aging, as reflected by increased TE expression in aged animals. Does increased TE expression during aging lead to more genomic TE copies in older animals? To address this question, we quantified TE Landscapes (TLs) via whole genome sequencing of young and aged Drosophila strains of wild-type and mutant backgrounds. We quantified TLs in whole flies and dissected brains and validated the feasibility of our approach in detecting new TE insertions in aging Drosophila genomes when small RNA and RNA interference (RNAi) pathways are compromised. We also describe improved sequencing methods to quantify extra-chromosomal DNA circles (eccDNAs) in Drosophila as an additional source of TE copies that accumulate during aging. Lastly, to combat the natural progression of aging-associated TE expression, we show that knocking down PAF1, a conserved transcription elongation factor that antagonizes RNAi pathways, may bolster suppression of TEs during aging and extend lifespan. Our study suggests that in addition to a possible influence by different genetic backgrounds, small RNA and RNAi mechanisms may mitigate genomic TL expansion despite the increase in TE transcripts during aging.

Transposable elements, also called transposons, are genetic parasites found in all animal genomes. Normally, transposons are compacted away in silent chromatin in young animals. But, as animals age and transposon-silencing defense mechanisms break down, transposon RNAs accumulate to significant levels in old animals like fruit flies. An open question is whether the increased levels of transposon RNAs in older animals also correspond to increased genomic copies of transposons. This study approached this question by sequencing the whole genomes of young and old wild-type and mutant flies lacking a functional RNA interference (RNAi) pathway, which naturally silences transposon RNAs. Although the wild-type flies with intact RNAi activity had little new accumulation of transposon copies, the sequencing approach was able to detect several transposon accumulation occurrences in some RNAi mutants. In addition, we found that some fly transposon families can also accumulate as extra-chromosomal circular DNA copies. Lastly, we showed that genetically augmenting the expression of RNAi factors can counteract the rising transposon RNA levels in aging and promote longevity. This study improves our understanding of the animal host genome relationship with transposons during natural aging processes.

Pathogenic variants in the human m6A reader YTHDC2 are associated with primary ovarian insufficiency

Primary ovarian insufficiency (POI) affects 1% of women and carries significant medical and psychosocial sequelae. Approximately 10% of POI has a defined genetic cause, with most implicated genes relating to biological processes involved in early fetal ovary development and function. Recently, Ythdc2, an RNA helicase and N6-methyladenosine (m6a) reader, has emerged as a novel regulator of meiosis in mice. Here, we describe homozygous pathogenic variants in YTHDC2 in three women with early-onset POI from two families: c. 2567C>G, p.P856R in the helicase-associated (HA2) domain; and c.1129G>T, p.E377*. We demonstrate that YTHDC2 is expressed in the developing human fetal ovary and is upregulated in meiotic germ cells, together with related meiosis-associated factors. The p.P856R variant results in a less flexible protein that likely disrupts downstream conformational kinetics of the HA2 domain, whereas the p.E377* variant truncates the helicase core. Taken together, our results reveal that YTHDC2 is a key new regulator of meiosis in humans and pathogenic variants within this gene are associated with POI.

Transposable element profiles reveal cell line identity and loss of heterozygosity in Drosophila cell culture

Cell culture systems allow key insights into biological mechanisms yet suffer from irreproducible outcomes in part because of cross-contamination or mislabeling of cell lines. Cell line misidentification can be mitigated by the use of genotyping protocols, which have been developed for human cell lines but are lacking for many important model species. Here, we leverage the classical observation that transposable elements (TEs) proliferate in cultured Drosophila cells to demonstrate that genome-wide TE insertion profiles can reveal the identity and provenance of Drosophila cell lines. We identify multiple cases where TE profiles clarify the origin of Drosophila cell lines (Sg4, mbn2, and OSS_E) relative to published reports, and also provide evidence that insertions from only a subset of long-terminal repeat retrotransposon families are necessary to mark Drosophila cell line identity. We also develop a new bioinformatics approach to detect TE insertions and estimate intra-sample allele frequencies in legacy whole-genome sequencing data (called ngs_te_mapper2), which revealed loss of heterozygosity as a mechanism shaping the unique TE profiles that identify Drosophila cell lines. Our work contributes to the general understanding of the forces impacting metazoan genomes as they evolve in cell culture and paves the way for high-throughput protocols that use TE insertions to authenticate cell lines in Drosophila and other organisms.

Genome-wide mapping of Piwi association with specific loci in Drosophila ovaries

Small noncoding RNA pathways have been implicated in diverse mechanisms of gene regulation. In Drosophila ovaries, Piwi binds to Piwi-interacting RNAs (piRNAs) of mostly 24–28 nucleotides (nt) and plays an important role in germline stem cell maintenance, transposon repression, and epigenetic regulation. To understand the mechanism underlying these functions, we report the application of the DamID-seq method to identify genome-wide binding sites of Piwi in Drosophila ovaries. Piwi localizes to at least 4535 euchromatic regions that are enriched with piRNA target sites. Surprisingly, the density of Piwi binding to euchromatin is much higher than in heterochromatin. Disrupting the piRNA binding of Piwi results in an overall change of the genomic binding profile, which indicates the role of piRNAs in directing Piwi to specific genomic sites. Most Piwi binding sites were either within or in the vicinity of protein-coding genes, particularly enriched near the transcriptional start and termination sites. The methylation signal near the transcriptional termination sites is significantly reduced when Piwi was mutated to become defective in piRNA binding. These observations indicate that Piwi might directly regulate the expression of many protein-coding genes, especially through regulating the 3' ends of targeted transcripts.

Maternal Piwi Regulates Primordial Germ Cell Development to Ensure the Fertility of Female Progeny in Drosophila

In many animals, germline development is initiated by proteins and RNAs that are expressed maternally. PIWI proteins and their associated small noncoding PIWI-interacting RNAs (piRNAs), which guide PIWI to target RNAs by base-pairing, are among the maternal components deposited into the germline of the Drosophila early embryo. Piwi has been extensively studied in the adult ovary and testis, where it is required for transposon suppression, germline stem cell self-renewal, and fertility. Consequently, loss of Piwi in the adult ovary using piwi-null alleles or knockdown from early oogenesis results in complete sterility, limiting investigation into possible embryonic functions of maternal Piwi. In this study, we show that the maternal Piwi protein persists in the embryonic germline through gonad coalescence, suggesting that maternal Piwi can regulate germline development beyond early embryogenesis. Using a maternal knockdown strategy, we find that maternal Piwi is required for the fertility and normal gonad morphology of female, but not male, progeny. Following maternal piwi knockdown, transposons were mildly derepressed in the early embryo but were fully repressed in the ovaries of adult progeny. Furthermore, the maternal piRNA pool was diminished, reducing the capacity of the PIWI/piRNA complex to target zygotic genes during embryogenesis. Examination of embryonic germ cell proliferation and ovarian gene expression showed that the germline of female progeny was partially masculinized by maternal piwi knockdown. Our study reveals a novel role for maternal Piwi in the germline development of female progeny and suggests that the PIWI/piRNA pathway is involved in germline sex determination in Drosophila.

Roles of piRNAs in transposon and pseudogene regulation of germline mRNAs and lncRNAs

PIWI proteins, a subfamily of PAZ/PIWI Domain family RNA-binding proteins, are best known for their function in silencing transposons and germline development by partnering with small noncoding RNAs called PIWI-interacting RNAs (piRNAs). However, recent studies have revealed multifaceted roles of the PIWI-piRNA pathway in regulating the expression of other major classes of RNAs in germ cells. In this review, we summarize how PIWI proteins and piRNAs regulate the expression of many disparate RNAs, describing a highly complex global genomic regulatory relationship at the RNA level through which piRNAs functionally connect all major constituents of the genome in the germline.

Adaptive evolution among cytoplasmic piRNA proteins leads to decreased genomic auto-immunity

In metazoan germlines, the piRNA pathway acts as a genomic immune system, employing small RNA-mediated silencing to defend host DNA from the harmful effects of transposable elements (TEs). Expression of genomic TEs is proposed to initiate self regulation by increasing the production of repressive piRNAs, thereby “adapting” piRNA-mediated control to the most active TE families. Surprisingly, however, piRNA pathway proteins, which execute piRNA biogenesis and enforce silencing of targeted sequences, evolve rapidly and adaptively in animals. If TE silencing is ensured through piRNA biogenesis, what necessitates changes in piRNA pathway proteins? Here we used interspecific complementation to test for functional differences between Drosophila melanogaster and D. simulans alleles of three adaptively evolving piRNA pathway proteins: Armitage, Aubergine and Spindle-E. In contrast to piRNA-mediated transcriptional regulators examined in previous studies, these three proteins have cytoplasmic functions in piRNA maturation and post-transcriptional silencing. Across all three proteins we observed interspecific divergence in the regulation of only a handful of TE families, which were more robustly silenced by the heterospecific piRNA pathway protein. This unexpected result suggests that unlike transcriptional regulators, positive selection has not acted on cytoplasmic piRNA effector proteins to enhance their function in TE repression. Rather, TEs may evolve to “escape” silencing by host proteins. We further discovered that D. simulans alleles of aub and armi exhibit enhanced off-target effects on host transcripts in a D. melanogaster background, as well as modest reductions in the efficiency of piRNA biogenesis, suggesting that promiscuous binding of D. simulans Aub and Armi proteins to host transcripts reduces their participation in piRNA production. Avoidance of genomic auto-immunity may therefore be a critical target of selection. Our observations suggest that piRNA effector proteins are subject to an evolutionary trade-off between defending the host genome from the harmful effect of TEs while also minimizing collateral damage to host genes.

Transposable elements are mobile fragments of selfish DNA that burden host genomes with deleterious mutations and incite genome instability. Host cells employ a specialized small-RNA mediated silencing pathway, the piRNA pathway, to act as a genomic immune system suppressing the mobilization of TEs. Changes in genomic TE content are met with rapid changes in the piRNA pool, thereby maintaining host control over transposition. However, piRNA pathway proteins—which enact piRNA biogenesis and silence target TEs—also evolve adaptively. To isolate forces that underlie this adaptive evolution, we examined functional divergence between two Drosophila species for three adaptively evolving piRNA pathway proteins. To our surprise, we found very few differences in TE regulation, suggesting that evolution has not generally acted to enhance control of TE parasites. Rather, we discovered interspecific differences in the regulation of host mRNAs for two proteins, which suggested that proteins evolve to avoid off-target silencing of host transcripts. We propose that the avoidance of such “genomic autoimmunity” is an important and underappreciated force driving the adaptive evolution of piRNA proteins.

Har-P, a short P-element variant, weaponizes P-transposase to severely impair Drosophila development

Without transposon-silencing Piwi-interacting RNAs (piRNAs), transposition causes an ovarian atrophy syndrome in Drosophila called gonadal dysgenesis (GD). Harwich (Har) strains with P-elements cause severe GD in F1 daughters when Har fathers mate with mothers lacking P-element-piRNAs (i.e. ISO1 strain). To address the mystery of why Har induces severe GD, we bred hybrid Drosophila with Har genomic fragments into the ISO1 background to create HISR-D or HISR-N lines that still cause Dysgenesis or are Non-dysgenic, respectively. In these lines, we discovered a highly truncated P-element variant we named ‘Har-P’ as the most frequent de novo insertion. Although HISR-D lines still contain full-length P-elements, HISR-N lines lost functional P-transposase but retained Har-P’s that when crossed back to P-transposase restores GD induction. Finally, we uncovered P-element-piRNA-directed repression on Har-P’s transmitted paternally to suppress somatic transposition. The Drosophila short Har-P’s and full-length P-elements relationship parallels the MITEs/DNA-transposase in plants and SINEs/LINEs in mammals.

DNA provides the instructions needed for life, a role that relies on it being a very stable and organized molecule. However, some sections of DNA are able to move from one place in the genome to another. When these “mobile genetic elements” move they may disrupt other genes and cause disease. For example, a mobile section of DNA known as the P-element causes a condition called gonadal dysgenesis in female fruit flies, leading to infertility.

Only certain strains of fruit flies carry P-elements and the severity of gonadal dysgenesis in their daughters varies. For example, when male fruit flies of a strain known as Harwich (or Har for short) is crossed with female fruit flies that do not contain P-elements, all of their daughters develop severe gonadal dysgenesis and are infertile. However, if the cross is done the other way around, and female Har flies mate with males that do not contain P-elements, the daughters are fertile because the Har mothers provide their daughters with protective molecules that silence the P-elements. But it was a mystery as to why the P-elements from the Har fathers always caused such severe gonadal dysgenesis in all the daughters.

Here, Srivastav et al. bred fruit flies to create offspring that had different pieces of Har DNA in a genetic background that was normally free from P-elements; they then analyzed the ‘hybrid’ offspring to identify which pieces of the Har genome caused gonadal dysgenesis in the daughter flies. These experiments showed that Har flies possess a very short variant of the P-element (named “Har-P”) that is more mobile than other variants. However, the Har-P variants still depended on an enzyme known as P-transposase encoded by the full-length P-elements to move around the genome. Further experiments showed that other strains of fruit flies that cause severe gonadal dysgenesis also had very short P-element variants that were almost identical to Har-P.

These findings may explain why Har and some other strains of fruit flies drive severe gonadal dysgenesis. In the future, it may be possible to transfer P-transposase and Har-P into mosquitoes, ticks and other biting insects to make them infertile and help reduce the spread of certain diseases in humans.

Sex and the TEs: transposable elements in sexual development and function in animals

Transposable elements are endogenous DNA sequences able to integrate into and multiply within genomes. They constitute a major source of genetic innovations, as they can not only rearrange genomes but also spread ready-to-use regulatory sequences able to modify host gene expression, and even can give birth to new host genes. As their evolutionary success depends on their vertical transmission, transposable elements are intrinsically linked to reproduction. In organisms with sexual reproduction, this implies that transposable elements have to manifest their transpositional activity in germ cells or their progenitors. The control of sexual development and function can be very versatile, and several studies have demonstrated the implication of transposable elements in the evolution of sex. In this review, we report the functional and evolutionary relationships between transposable elements and sexual reproduction in animals. In particular, we highlight how transposable elements can influence expression of sexual development genes, and how, reciprocally, they are tightly controlled in gonads. We also review how transposable elements contribute to the organization, expression and evolution of sexual development genes and sex chromosomes. This underscores the intricate co-evolution between host functions and transposable elements, which regularly shift from a parasitic to a domesticated status useful to the host.

Long-read direct RNA sequencing by 5'-Cap capturing reveals the impact of Piwi on the widespread exonization of transposable elements in locusts

The large genome of the migratory locust (Locusta migratoria) genome accumulates massive amount of accumulated transposable elements (TEs), which show intrinsic transcriptional activities. Hampering the ability to precisely determine full-length RNA transcript sequences are exonized TEs, which produce numerous highly similar fragments that are difficult to resolve using short-read sequencing technology. Here, we applied a 5'-Cap capturing method using Nanopore long-read direct RNA sequencing to characterize full-length transcripts in their native RNA form and to analyze the TE exonization pattern in the locust transcriptome. Our results revealed the widespread establishment of TE exonization and a substantial contribution of TEs to RNA splicing in the locust transcriptome. The results of the transcriptomic spectrum influenced by Piwi expression indicated that TE-derived sequences were the main targets of Piwi-mediated repression. Furthermore, our study showed that Piwi expression regulates the length of RNA transcripts containing TE-derived sequences, creating an alternative UTR usage. Overall, our results reveal the transcriptomic characteristics of TE exonization in the species characterized by large and repetitive genomes.

Interplay between small RNA pathways shapes chromatin landscapes in C. elegans

The nematode Caenorhabditis elegans contains several types of endogenous small interfering RNAs (endo-siRNAs) produced by RNA-dependent RNA polymerase (RdRP) complexes. Both 'silencing' siRNAs bound by Worm-specific Argonautes (WAGO) and 'activating' siRNAs bound by the CSR-1 Argonaute require the DRH-3 helicase, an RdRP component. Here, we show that, in the drh-3(ne4253) mutant deficient in RdRP-produced secondary endo-siRNAs, the silencing histone mark H3K9me3 is largely depleted, whereas in the csr-1 partially rescued null mutant strain (WM193), this mark is ectopically deposited on CSR-1 target genes. Moreover, we observe ectopic H3K9me3 at enhancer elements and an increased number of small RNAs that match enhancers in both drh-3 and csr-1 mutants. Finally, we detect accumulation of H3K27me3 at highly expressed genes in the drh-3(ne4253) mutant, which correlates with their reduced transcription. Our study shows that when abundant RdRP-produced siRNAs are depleted, there is ectopic elevation of noncoding RNAs linked to sites with increased silencing chromatin marks. Moreover, our results suggest that enhancer small RNAs may guide local H3K9 methylation.

A unique cluster of roo insertions in the promoter region of a stress response gene in Drosophila melanogaster

Transposable elements (TEs) are not randomly distributed in the genome. A genome-wide analysis of the D. melanogaster genome found that differences in TE density across 50 kb genomic regions was due both to transposition and duplication. At smaller genomic scales, promoter regions of hsp genes and the promoter region of CG18446 have been shown to accumulate TE insertions. In this work, we have further analyzed the promoter region of CG18446. We screened 218 strains collected in 15 natural populations, and we found that the CG18446 promoter region contains 20 independent roo insertions. Based on phylogenetic analysis, we suggest that the presence of multiple roo insertions in this region is likely to be the result of several bursts of transposition. Moreover, we found that the roo insertional cluster in the CG18446 promoter region is unique: no other promoter region in the genome contains a similar number of roo insertions. We found that, similar to hsp gene promoters, chromatin accessibility could be one of the factors explaining the recurrent insertions of roo elements in CG18446 promoter region.

Genome-wide identification and differentially expression analysis of lncRNAs in tilapia

Background

Long noncoding RNAs (LncRNAs) play important roles in fundamental biological processes. However, knowledge about the genome-wide distribution and stress-related expression of lncRNAs in tilapia is still limited.

Results

Genome-wide identification of lncRNAs in the tilapia genome was carried out in this study using bioinformatics tools. 103 RNAseq datasets that generated in our laboratory or collected from NCBI database were analyzed. In total, 72,276 high-confidence lncRNAs were identified. The averaged positive correlation coefficient (r_mean = 0.286) between overlapped lncRNA and mRNA pairs showed significant differences with the values for all lncRNA-mRNA pairs (r_mean = 0.176, z statistics = − 2.45, p value = 0.00071) and mRNA-mRNA pairs (r_mean = 0.186, z statistics = − 2.23, p value = 0.0129). Weighted correlation network analysis of the lncRNA and mRNA datasets from 12 tissues identified 21 modules and many interesting mRNA genes that clustered with lncRNAs. Overrepresentation test indicated that these mRNAs enriched in many biological processes, such as meiosis (p = 0.00164), DNA replication (p = 0.00246), metabolic process (p = 0.000838) and in molecular function, e.g., helicase activity (p = 0.000102) and catalytic activity (p = 0.0000612). Differential expression (DE) analysis identified 99 stress-related lncRNA genes and 1955 tissue-specific DE lncRNA genes. MiRNA-lncRNA interaction analysis detected 72,267 lncRNAs containing motifs with sequence complementary to 458 miRNAs.

Conclusions

This study provides an invaluable resource for further studies on molecular bases of lncRNAs in tilapia genomes. Further function analysis of the lncRNAs will help to elucidate their roles in regulating stress-related adaptation in tilapia.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5115-x) contains supplementary material, which is available to authorized users.

Nuclear lamina dysfunction triggers a germline stem cell checkpoint

LEM domain (LEM-D) proteins are conserved components of the nuclear lamina (NL) that contribute to stem cell maintenance through poorly understood mechanisms. The Drosophila emerin homolog Otefin (Ote) is required for maintenance of germline stem cells (GSCs) and gametogenesis. Here, we show that ote mutants carry germ cell-specific changes in nuclear architecture that are linked to GSC loss. Strikingly, we found that both GSC death and gametogenesis are rescued by inactivation of the DNA damage response (DDR) kinases, ATR and Chk2. Whereas the germline checkpoint draws from components of the DDR pathway, genetic and cytological features of the GSC checkpoint differ from the canonical pathway. Instead, structural deformation of the NL correlates with checkpoint activation. Despite remarkably normal oogenesis, rescued oocytes do not support embryogenesis. Taken together, these data suggest that NL dysfunction caused by Otefin loss triggers a GSC-specific checkpoint that contributes to maintenance of gamete quality.

Otefin is a nuclear lamina protein required for survival of Drosophila germ stem cells. Here the authors show that nuclear lamina dysfunction resulting from loss of Otefin activates a DNA damage-independent germ stem cell-specific checkpoint, mediated by the ATR and Chk2 kinases, which ensures that healthy gametes are passed on to the next generation.

MIWI2 targets RNAs transcribed from piRNA-dependent regions to drive DNA methylation in mouse prospermatogonia

Argonaute/Piwi proteins can regulate gene expression via RNA degradation and translational regulation using small RNAs as guides. They also promote the establishment of suppressive epigenetic marks on repeat sequences in diverse organisms. In mice, the nuclear Piwi protein MIWI2 and Piwi-interacting RNAs (piRNAs) are required for DNA methylation of retrotransposon sequences and some other sequences. However, its underlying molecular mechanisms remain unclear. Here, we show that piRNA-dependent regions are transcribed at the stage when piRNA-mediated DNA methylation takes place. MIWI2 specifically interacts with RNAs from these regions. In addition, we generated mice with deletion of a retrotransposon sequence either in a representative piRNA-dependent region or in a piRNA cluster. Both deleted regions were required for the establishment of DNA methylation of the piRNA-dependent region, indicating that piRNAs determine the target specificity of MIWI2-mediated DNA methylation. Our results indicate that MIWI2 affects the chromatin state through base-pairing between piRNAs and nascent RNAs, as observed in other organisms possessing small RNA-mediated epigenetic regulation.

iCLIP of the PIWI Protein Aubergine in Drosophila Embryos

Piwi-interacting RNAs (piRNAs) are a class of small noncoding RNAs bound to specific Argonaute proteins, the PIWI proteins. piRNAs target mRNAs by complementarity to silence them; they play an important role in the repression of transposable elements in the germ line of many species. piRNAs and PIWI proteins are also involved in diverse biological processes through their role in the regulation of cellular mRNAs. In the Drosophila embryo, they contribute to the maternal mRNA decay occurring during the maternal-to-zygotic transition. CLIP (UV cross-linking and immunoprecipitation) techniques have been used to identify target mRNAs of Argonaute proteins. Here we describe the iCLIP (individual-nucleotide resolution CLIP) protocol that we have adapted for the PIWI protein Aubergine in Drosophila embryos.

Identification of transposons near predicted lncRNA and mRNA pools of Prunus mume using an integrative transposable element database constructed from Rosaceae plant genomes

This study focused on the construction of a database of transposable elements (TEs) from Rosaceae plants, the third most economically important plant family in temperate regions, and its transcriptomics applications. The evolutionary effects of TEs on gene regulation have been explored, and TE insertions can be the molecular bases of changes in gene structure and function. However, a specific Rosaceae plant TE database (RPTEdb) is lacking. The genomes of several Rosaceae plants have been sequenced, providing the opportunity to mine TE data at a whole-genome level. Therefore, we constructed the RPTEdb, a collective and comprehensive database of 19,596 annotated TEs in the genomes of Rosaceae plants using previously described identification and annotation methods and published genome sequences. The user-friendly web-based database provides access to research tools through hyperlinks, including Browse, TE tree, tools, JBrowse, and search sections, and through the inputting of sequences on the main webpage. Next, we performed one advanced application in which TEs near predicted long non-coding RNA (lncRNA) and mRNA domains within white and red petal-tissue transcriptomes of Prunus mume 'Fuban Tiaozhi' were identified, revealing 16 TEs that overlapped or were near 16 differentially expressed lncRNA domains, and 54 TEs that overlapped or were near 54 differentially expressed mRNA domains, and the TEs' possible functions were also discussed. We believe that the RPTEdb will contribute to the understanding of TE roles in the structural, functional and evolutionary dynamics of Rosaceae plant genomes.

Interactions between short and long noncoding RNAs

It is now evident that noncoding RNAs play key roles in regulatory networks determining cell fate and behavior, in a myriad of different conditions, and across all species. Among these noncoding RNAs are short RNAs, such as MicroRNAs, snoRNAs, and Piwi-interacting RNAs, and the functions of those are relatively well understood. Other noncoding RNAs are longer, and their modes of action and functions are also increasingly explored and deciphered. Short RNAs and long noncoding RNAs (lncRNAs) interact with each other with reciprocal consequences for their fates and functions. LncRNAs serve as precursors for many types of small RNAs and, therefore, the pathways for small RNA biogenesis can impinge upon the fate of lncRNAs. In addition, lncRNA expression can be repressed by small RNAs, and lncRNAs can affect small RNA activity and abundance through competition for binding or by triggering small RNA degradation. Here, I review the known types of interactions between small and long RNAs, discuss their outcomes, and bring representative examples from studies in mammals.

Piwi interacts with chromatin at nuclear pores and promiscuously binds nuclear transcripts in Drosophila ovarian somatic cells

Piwi in a complex with Piwi-interacting RNAs (piRNAs) triggers transcriptional silencing of transposable elements (TEs) in Drosophila ovaries, thus ensuring genome stability. To do this, Piwi must scan the nascent transcripts of genes and TEs for complementarity to piRNAs. The mechanism of this scanning is currently unknown. Here we report the DamID-seq mapping of multiple Piwi-interacting chromosomal domains in somatic cells of Drosophila ovaries. These domains significantly overlap with genomic regions tethered to Nuclear Pore Complexes (NPCs). Accordingly, Piwi was coimmunoprecipitated with the component of NPCs Elys and with the Xmas-2 subunit of RNA transcription and export complex, known to interact with NPCs. However, only a small Piwi fraction has transient access to DNA at nuclear pores. Importantly, although 36% of the protein-coding genes overlap with Piwi-interacting domains and RNA-immunoprecipitation results demonstrate promiscuous Piwi binding to numerous genic and TE nuclear transcripts, according to available data Piwi does not silence these genes, likely due to the absence of perfect base-pairing between piRNAs and their transcripts.

Drosophila PAF1 Modulates PIWI/piRNA Silencing Capacity

To test the directness of factors in initiating PIWI-directed gene silencing, we employed a Piwi-interacting RNA (piRNA)-targeted reporter assay in Drosophila ovary somatic sheet (OSS) cells [1]. This assay confirmed direct silencing roles for piRNA biogenesis factors and PIWI-associated factors [2-12] but suggested that chromatin-modifying proteins may act downstream of the initial silencing event. Our data also revealed that RNA-polymerase-II-associated proteins like PAF1 and RTF1 antagonize PIWI-directed silencing. PAF1 knockdown enhances PIWI silencing of reporters when piRNAs target the transcript region proximal to the promoter. Loss of PAF1 suppresses endogenous transposable element (TE) transcript maturation, whereas a subset of gene transcripts and long-non-coding RNAs adjacent to TE insertions are affected by PAF1 knockdown in a similar fashion to piRNA-targeted reporters. Additionally, transcription activation at specific TEs and TE-adjacent loci during PIWI knockdown is suppressed when PIWI and PAF1 levels are both reduced. Our study suggests a mechanistic conservation between fission yeast PAF1 repressing AGO1/small interfering RNA (siRNA)-directed silencing [13, 14] and Drosophila PAF1 opposing PIWI/piRNA-directed silencing.

Xenopus Piwi proteins interact with a broad proportion of the oocyte transcriptome

Piwi proteins utilize small RNAs (piRNAs) to recognize target transcripts such as transposable elements (TE). However, extensive piRNA sequence diversity also suggests that Piwi/piRNA complexes interact with many transcripts beyond TEs. To determine Piwi target RNAs, we used ribonucleoprotein-immunoprecipitation (RIP) and cross-linking and immunoprecipitation (CLIP) to identify thousands of transcripts associated with the Piwi proteins XIWI and XILI (Piwi-protein-associated transcripts, PATs) from early stage oocytes of X. laevis and X. tropicalis Most PATs associate with both XIWI and XILI and include transcripts of developmentally important proteins in oogenesis and embryogenesis. Only a minor fraction of PATs in both frog species displayed near perfect matches to piRNAs. Since predicting imperfect pairing between all piRNAs and target RNAs remains intractable, we instead determined that PAT read counts correlate well with the lengths and expression levels of transcripts, features that have also been observed for oocyte mRNAs associated with Drosophila Piwi proteins. We used an in vitro assay with exogenous RNA to confirm that XIWI associates with RNAs in a length- and concentration-dependent manner. In this assay, noncoding transcripts with many perfectly matched antisense piRNAs were unstable, whereas coding transcripts with matching piRNAs were stable, consistent with emerging evidence that Piwi proteins both promote the turnover of TEs and other RNAs, and may also regulate mRNA localization and translation. Our study suggests that Piwi proteins play multiple roles in germ cells and establishes a tractable vertebrate system to study the role of Piwi proteins in transcript regulation.

Transposable Element Bm1645 is a Source of BmAGO2-associated Small RNAs that affect its expression in Bombyx mori

Background

A transposable element (TE) is a DNA fragment that can change its position within a genome. Transposable elements play important roles in maintaining the stability and diversity of organisms by transposition. Recent studies have shown that approximately half of the genes in Bombyx mori are TEs.

Results

We systematically identified and analyzed the BmAGO2-associated TEs, which exceed 100 in the B. mori genome. Additionally, we also mapped the small RNAs associated with BmAGO2 in B.mori. The transposon Bm1645 is the most abundant TE associated with BmAGO2, and Bm1645-derived small RNAs represent a small RNA pool. We determined the expression patterns of several Bm1645-derived small RNAs by northern blotting, and the results showed there was differential expression of multiple small RNAs in normal and BmNPV-infected BmN cells and silkworms from various developmental stages. We confirmed that four TE-siRNAs could bind to BmAGO2 using EMSA and also validated the recognition sites of these four TE-siRNAs in Bm1645 by dual-luciferase reporter assays. Furthermore, qRT-PCR analysis revealed the overexpression of the four TE-siRNAs could downregulate the expression of Bm1645 in BmN cells, and the transcription of Bm1645 was upregulated by the downregulation of BmAGO2.

Conclusions

Our results suggest Bm1645 functions as a source of small RNAs pool and this pool can produce many BmAGO2-associated small RNAs that regulate TE’s expression.

Electronic supplementary material

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

Telomeric Retrotransposon HeT-A Contains a Bidirectional Promoter that Initiates Divergent Transcription of piRNA Precursors in Drosophila Germline

PIWI-interacting RNAs (piRNAs) provide the silencing of transposable elements in the germline. Drosophila telomeres are maintained by transpositions of specialized telomeric retroelements. piRNAs generated from sense and antisense transcripts of telomeric elements provide telomere length control in the germline. Previously, we have found that antisense transcription of the major telomeric retroelement HeT-A is initiated upstream of the HeT-A sense transcription start site. Here, we performed a deletion analysis of the HeT-A promoter and show that common regulatory elements are shared by sense and antisense promoters of HeT-A. Therefore, the HeT-A promoter is a bidirectional promoter capable of processive sense and antisense transcription. Ovarian small RNA data show that a solo HeT-A promoter within an euchromatic transgene initiates the divergent transcription of transgenic reporter genes and subsequent processing of these transcripts into piRNAs. These events lead to the formation of a divergent unistrand piRNA cluster at solo HeT-A promoters, in contrast to endogenous telomeres that represent strong dual-strand piRNA clusters. Solo HeT-A promoters are not immunoprecipitated with heterochromatin protein 1 (HP1) homolog Rhino, a marker of the dual-strand piRNA clusters, but are associated with HP1 itself, which provides piRNA-mediated transcriptional repression of the reporter genes. Unlike endogenous dual-strand piRNA clusters, the solo HeT-A promoter does not produce overlapping transcripts. In a telomeric context, however, bidirectional promoters of tandem HeT-A repeats provide a read-through transcription of both genomic strands, followed by Rhi binding. These data indicate that Drosophila telomeres share properties of unistrand and dual-strand piRNA clusters.

Small RNAs from a Big Genome: The piRNA Pathway and Transposable Elements in the Salamander Species Desmognathus fuscus

Most of the largest vertebrate genomes are found in salamanders, a clade of amphibians that includes 686 species. Salamander genomes range in size from 14 to 120 Gb, reflecting the accumulation of large numbers of transposable element (TE) sequences from all three TE classes. Although DNA loss rates are slow in salamanders relative to other vertebrates, high levels of TE insertion are also likely required to explain such high TE loads. Across the Tree of Life, novel TE insertions are suppressed by several pathways involving small RNA molecules. In most known animals, TE activity in the germline is primarily regulated by the Piwi-interacting RNA (piRNA) pathway. In this study, we test the hypothesis that salamanders' unusually high TE loads reflect the loss of the ancestral piRNA-mediated TE-silencing machinery. We characterized the small RNA pool in the female and male adult gonads, testing for the presence of small RNA molecules that bear the characteristics of TE-targeting piRNAs. We also analyzed the amino acid sequences of piRNA pathway proteins from salamanders and other vertebrates, testing whether the overall patterns of sequence divergence are consistent with conserved pathway function across the vertebrate clade. Our results do not support the hypothesis of piRNA pathway loss; instead, they suggest that the piRNA pathway is expressed in salamanders. Given these results, we propose hypotheses to explain how the extraordinary TE loads in salamander genomes could have accumulated, despite the expression of TE-silencing machinery.

The Role of PIWIL4, an Argonaute Family Protein, in Breast Cancer

P-element-induced wimpy testis (PIWI) proteins bind to PIWI-interacting RNAs and play key roles in the biogenesis and functions of PIWI-interacting RNAs. It has been reported that PIWI proteins are essential for stem cell self-renewal and germline development in diverse organisms and that they are ectopically expressed in multiple forms of cancer. However, the role of PIWI in cancer remains elusive. Here we report that one of the four PIWI proteins in humans, PIWIL4, is highly expressed in both breast cancer tissues and the cytoplasm of MDA-MB-231 cells derived from breast cancer. Reducing PIWIL4 expression drastically impairs the migration ability of MDA-MB-231 cells, significantly increases their apoptosis, and mildly affects their proliferation. Our transcriptome and proteome analysis reveal that these functions are at least partially achieved via the PIWIL4 regulation of TGF-β and FGF signaling pathways and MHC class II proteins. These findings suggest that PIWIL4 may serve as a potential therapeutic target for breast cancer.

Zfrp8/PDCD2 Interacts with RpS2 Connecting Ribosome Maturation and Gene-Specific Translation

Zfrp8/PDCD2 is a highly conserved protein essential for stem cell maintenance in both flies and mammals. It is also required in fast proliferating cells such as cancer cells. Our previous studies suggested that Zfrp8 functions in the formation of mRNP (mRNA ribonucleoprotein) complexes and also controls RNA of select Transposable Elements (TEs). Here we show that in Zfrp8/PDCD2 knock down (KD) ovaries, specific mRNAs and TE transcripts show increased nuclear accumulation. We also show that Zfrp8/PDCD2 interacts with the (40S) small ribosomal subunit through direct interaction with RpS2 (uS5). By studying the distribution of endogenous and transgenic fluorescently tagged ribosomal proteins we demonstrate that Zfrp8/PDCD2 regulates the cytoplasmic levels of components of the small (40S) ribosomal subunit, but does not control nuclear/nucleolar localization of ribosomal proteins. Our results suggest that Zfrp8/PDCD2 functions at late stages of ribosome assembly and may regulate the binding of specific mRNA-RNPs to the small ribosomal subunit ultimately controlling their cytoplasmic localization and translation.

Piwi maintains germline stem cells and oogenesis in Drosophila through negative regulation of Polycomb group proteins

The Drosophila Piwi protein regulates both niche and intrinsic mechanisms to maintain germline stem cells, but its underlying mechanism remains unclear. Here we report that Piwi cooperates with Polycomb Group complexes PRC1 and PRC2 in niche and germline cells to regulate ovarian germline stem cells and oogenesis. Piwi physically interacts with PRC2 subunits Su(z)12 and Esc in the ovary and in vitro. Chromatin co-immunoprecipitation of Piwi, the PRC2 enzymatic subunit E(z), lysine-27-tri-methylated histone 3 (H3K27m3), and RNA polymerase II in wild-type and piwi mutant ovaries reveals that Piwi binds a conserved DNA motif at ~72 genomic sites, and inhibits PRC2 binding to many non-Piwi-binding genomic targets and H3K27 tri-methylation. Moreover, Piwi influences RNA Polymerase II activities in Drosophila ovaries likely via inhibiting PRC2. We hypothesize that Piwi negatively regulates PRC2 binding by sequestering PRC2 in the nucleoplasm, thus reducing PRC2 binding to many targets and influences transcription during oogenesis.

Unique transposon landscapes are pervasive across Drosophila melanogaster genomes

To understand how transposon landscapes (TLs) vary across animal genomes, we describe a new method called the Transposon Insertion and Depletion AnaLyzer (TIDAL) and a database of >300 TLs in Drosophila melanogaster (TIDAL-Fly). Our analysis reveals pervasive TL diversity across cell lines and fly strains, even for identically named sub-strains from different laboratories such as the ISO1 strain used for the reference genome sequence. On average, >500 novel insertions exist in every lab strain, inbred strains of the Drosophila Genetic Reference Panel (DGRP), and fly isolates in the Drosophila Genome Nexus (DGN). A minority (<25%) of transposon families comprise the majority (>70%) of TL diversity across fly strains. A sharp contrast between insertion and depletion patterns indicates that many transposons are unique to the ISO1 reference genome sequence. Although TL diversity from fly strains reaches asymptotic limits with increasing sequencing depth, rampant TL diversity causes unsaturated detection of TLs in pools of flies. Finally, we show novel transposon insertions negatively correlate with Piwi-interacting RNA (piRNA) levels for most transposon families, except for the highly-abundant roo retrotransposon. Our study provides a useful resource for Drosophila geneticists to understand how transposons create extensive genomic diversity in fly cell lines and strains.

Panoramix enforces piRNA-dependent cotranscriptional silencing

The Piwi-interacting RNA (piRNA) pathway is a small RNA-based innate immune system that defends germ cell genomes against transposons. In Drosophila ovaries, the nuclear Piwi protein is required for transcriptional silencing of transposons, though the precise mechanisms by which this occurs are unknown. Here we show that the CG9754 protein is a component of Piwi complexes that functions downstream of Piwi and its binding partner, Asterix, in transcriptional silencing. Enforced tethering of CG9754 to nascent messenger RNA transcripts causes cotranscriptional silencing of the source locus and the deposition of repressive chromatin marks. We have named CG9754 "Panoramix," and we propose that this protein could act as an adaptor, scaffolding interactions between the piRNA pathway and the general silencing machinery that it recruits to enforce transcriptional repression.

Genome-wide identification of long noncoding RNA genes and their potential association with fecundity and virulence in rice brown planthopper, Nilaparvata lugens

Background

The functional repertoire of long noncoding RNA (lncRNA) has been characterized in several model organisms, demonstrating that lncRNA plays important roles in fundamental biological processes. However, they remain largely unidentified in most species. Understanding the characteristics and functions of lncRNA in insects would be useful for insect resources utilization and sustainable pest control.

Methods

A computational pipeline was developed to identify lncRNA genes in the rice brown planthopper, Nilaparvata lugens, a destructive rice pest causing huge yield losses. Strand specific RT-PCR were used to determine the transcription orientation of lncRNAs.

Results

In total, 2,439 lncRNA transcripts corresponding to 1,882 loci were detected from 12 whole transcriptomes (RNA-seq) datasets, including samples from high fecundity (HFP), low fecundity (LFP), I87i and C89i populations, in addition Mudgo and TN1 virulence strains. The identified N. lugens lncRNAs had low sequence similarities with other known lncRNAs. However, their structural features were similar with mammalian counterparts. N. lugens lncRNAs had shorter transcripts than protein-coding genes due to the lower exon number though their exons and introns were longer. Only 19.9% of N. lugens lncRNAs had multiple alternatively spliced isoforms. We observed biases in the genome location of N. lugens lncRNAs. More than 30% of the lncRNAs overlapped with known protein-coding genes. These lncRNAs tend to be co-expressed with their neighboring genes (Pearson correlation, p < 0.01, T-test) and might interact with adjacent protein-coding genes. In total, 19-148 lncRNAs were specifically-expressed in the samples of HFP, LFP, Mudgo, TN1, I87i and C89i populations. Three lncRNAs specifically expressed in HFP and LFP populations overlapped with reproductive-associated genes.

Discussion

The structural features of N. lugens lncRNAs are similar to mammalian counterparts. Coexpression and function analysis suggeste that N. lugens lncRNAs might have important functions in high fecundity and virulence adaptability.

Conclusions

This study provided the first catalog of lncRNA genes in rice brown planthopper. Gene expression and genome location analysis indicated that lncRNAs might play important roles in high fecundity and virulence adaptation in N. lugens.

Electronic supplementary material

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

Long Noncoding RNAs as New Architects in Cancer Epigenetics, Prognostic Biomarkers, and Potential Therapeutic Targets

Recent advances in genome-wide analysis have revealed that 66% of the genome is actively transcribed into noncoding RNAs (ncRNAs) while less than 2% of the sequences encode proteins. Among ncRNAs, high-resolution microarray and massively parallel sequencing technologies have identified long ncRNAs (>200 nucleotides) that lack coding protein function. LncRNAs abundance, nuclear location, and diversity allow them to create in association with protein interactome, a complex regulatory network orchestrating cellular phenotypic plasticity via modulation of all levels of protein-coding gene expression. Whereas lncRNAs biological functions and mechanisms of action are still not fully understood, accumulating data suggest that lncRNAs deregulation is pivotal in cancer initiation and progression and metastatic spread through various mechanisms, including epigenetic effectors, alternative splicing, and microRNA-like molecules. Mounting data suggest that several lncRNAs expression profiles in malignant tumors are associated with prognosis and they can be detected in biological fluids. In this review, we will briefly discuss characteristics and functions of lncRNAs, their role in carcinogenesis, and their potential usefulness as diagnosis and prognosis biomarkers and novel therapeutic targets.

The Discovery, Distribution, and Evolution of Viruses Associated with Drosophila melanogaster

Drosophila melanogaster is a valuable invertebrate model for viral infection and antiviral immunity, and is a focus for studies of insect-virus coevolution. Here we use a metagenomic approach to identify more than 20 previously undetected RNA viruses and a DNA virus associated with wild D. melanogaster. These viruses not only include distant relatives of known insect pathogens but also novel groups of insect-infecting viruses. By sequencing virus-derived small RNAs, we show that the viruses represent active infections of Drosophila. We find that the RNA viruses differ in the number and properties of their small RNAs, and we detect both siRNAs and a novel miRNA from the DNA virus. Analysis of small RNAs also allows us to identify putative viral sequences that lack detectable sequence similarity to known viruses. By surveying >2,000 individually collected wild adult Drosophila we show that more than 30% of D. melanogaster carry a detectable virus, and more than 6% carry multiple viruses. However, despite a high prevalence of the Wolbachia endosymbiont--which is known to be protective against virus infections in Drosophila--we were unable to detect any relationship between the presence of Wolbachia and the presence of any virus. Using publicly available RNA-seq datasets, we show that the community of viruses in Drosophila laboratories is very different from that seen in the wild, but that some of the newly discovered viruses are nevertheless widespread in laboratory lines and are ubiquitous in cell culture. By sequencing viruses from individual wild-collected flies we show that some viruses are shared between D. melanogaster and D. simulans. Our results provide an essential evolutionary and ecological context for host-virus interaction in Drosophila, and the newly reported viral sequences will help develop D. melanogaster further as a model for molecular and evolutionary virus research.

The capacity of target silencing by Drosophila PIWI and piRNAs

Although Piwi proteins and Piwi-interacting RNAs (piRNAs) genetically repress transposable elements (TEs), it is unclear how the highly diverse piRNA populations direct Piwi proteins to silence TE targets without silencing the entire transcriptome. To determine the capacity of piRNA-mediated silencing, we introduced reporter genes into Drosophila OSS cells, which express microRNAs (miRNAs) and piRNAs, and compared the Piwi pathway to the Argonaute pathway in gene regulation. Reporter constructs containing several target sites that were robustly silenced by miRNAs were not silenced to the same degrees by piRNAs. However, another set of reporters we designed to enable a large number of both TE-directed and genic piRNAs to bind were robustly silenced by the PIWI/piRNA complex in OSS cells. These reporters show that a bulk of piRNAs are required to pair to the reporter's transcripts and not the reporter's DNA sequence to engage PIWI-mediated silencing. Following our genome-wide study of PIWI-regulated targets in OSS cells, we assessed candidate gene elements with our reporter platform. These results suggest TE sequences are the most direct of PIWI regulatory targets while coding genes are less directly affected by PIWI targeting. Finally, our study suggests that the PIWI transcriptional silencing mechanism triggers robust chromatin changes on targets with sufficient piRNA binding, and preferentially regulates TE transcripts because protein-coding transcripts lack a threshold of targeting by piRNA populations. This reporter platform will facilitate future dissections of the PIWI-targeting mechanism.

Piwi Proteins and piRNAs step onto the systems biology stage

Animal germ cells are totipotent because they maintain a highly unique and specialized epigenetic state for its genome. To accomplish this, germ cells express a rich repertoire of specialized RNA-binding protein complexes such as the Piwi proteins and Piwi-interacting RNAs (piRNAs): a germ-cell branch of the RNA interference (RNAi) phenomenon which includes microRNA and endogenous small interfering RNA pathways. Piwi proteins and piRNAs are deeply conserved in animal evolution and play essential roles in fertility and regeneration. Molecular mechanisms for how these ribonucleoproteins act upon the transcriptome and genome are only now coming to light with the application of systems-wide approaches in both invertebrates and vertebrates. Systems biology studies on invertebrates have revealed that transcriptional and heritable silencing is a main mechanism driven by Piwi proteins and piRNA complexes. In vertebrates, Piwi-targeting mechanisms and piRNA biogenesis have progressed, while the discovery that the nuclease activity of Piwi protein is essential for vertebrate germ cell development but not completely required in invertebrates highlights the many complexities of this pathway in different animals. This review recounts how recent systems-wide approaches have rapidly accelerated our appreciation for the broad reach of the Piwi pathway on germline genome regulation and what questions facing the field await to be unraveled.