Terminator oligo blocking efficiently eliminates rRNA from Drosophila small RNA sequencing libraries

DNMT3L interacts with Piwi and modulates the expression of piRNAs in transgenic Drosophila

Aim: To explore the role of Piwi protein and piRNAs in DNMT3L-mediated epigenetic inheritance. Materials & methods: Transgenic Drosophila were used to examine the effect of ectopically expressed DNMT3L on the profile of piRNAs by sequencing of small RNAs. Results & conclusion: Our previous work showed accumulation and inheritance of epimutations across multiple generations in transgenic DNMT3L Drosophila. Here, we show interaction of DNMT3L with Piwi and a significant alteration in the piRNA profile across multiple generations in transgenic Drosophila. In the light of its interaction with histone H1, we propose that in addition to its role of modulating core histone modifications, DNMT3L allows for inheritance of epigenetic information through its collaboration with Piwi, piRNAs and histone H1.

Transcriptional programs are activated and microRNAs are repressed within minutes after mating in the Drosophila melanogaster female reproductive tract

BACKGROUND

The female reproductive tract is exposed directly to the male's ejaculate, making it a hotspot for mating-induced responses. In Drosophila melanogaster, changes in the reproductive tract are essential to optimize fertility. Many changes occur within minutes after mating, but such early timepoints are absent from published RNA-seq studies. We measured transcript abundances using RNA-seq and microRNA-seq of reproductive tracts of unmated and mated females collected at 10-15 min post-mating. We further investigated whether early transcriptome changes in the female reproductive tract are influenced by inhibiting BMPs in secondary cells, a condition that depletes exosomes from the male's ejaculate.

RESULTS

We identified 327 differentially expressed genes. These were mostly upregulated post-mating and have roles in tissue morphogenesis, wound healing, and metabolism. Differentially abundant microRNAs were mostly downregulated post-mating. We identified 130 predicted targets of these microRNAs among the differentially expressed genes. We saw no detectable effect of BMP inhibition in secondary cells on transcript levels in the female reproductive tract.

CONCLUSIONS

Our results indicate that mating induces early changes in the female reproductive tract primarily through upregulation of target genes, rather than repression. The upregulation of certain target genes might be mediated by the mating-induced downregulation of microRNAs. Male-derived exosomes and other BMP-dependent products were not uniquely essential for this process. Differentially expressed genes and microRNAs provide candidates that can be further examined for their participation in the earliest alterations of the reproductive tract microenvironment.

P-element invasion fuels molecular adaptation in laboratory populations of Drosophila melanogaster

Transposable elements (TEs) are mobile genetic parasites that frequently invade new host genomes through horizontal transfer. Invading TEs often exhibit a burst of transposition, followed by reduced transposition rates as repression evolves in the host. We recreated the horizontal transfer of P-element DNA transposons into a Drosophila melanogaster host and followed the expansion of TE copies and evolution of host repression in replicate laboratory populations reared at different temperatures. We observed that while populations maintained at high temperatures rapidly go extinct after TE invasion, those maintained at lower temperatures persist, allowing for TE spread and the evolution of host repression. We also surprisingly discovered that invaded populations experienced recurrent insertion of P-elements into a specific long non-coding RNA, lncRNA:CR43651, and that these insertion alleles are segregating at unusually high frequency in experimental populations, indicative of positive selection. We propose that, in addition to driving the evolution of repression, transpositional bursts of invading TEs can drive molecular adaptation.

Stress-sensitive dynamics of miRNAs and Elba1 in Drosophila embryogenesis

Early-life stress can result in life-long effects that impact adult health and disease risk, but little is known about how such programming is established and maintained. Here, we show that such epigenetic memories can be initiated in the Drosophila embryo before the major wave of zygotic transcription, and higher-order chromatin structures are established. An early short heat shock results in elevated levels of maternal miRNA and reduced levels of a subgroup of zygotic genes in stage 5 embryos. Using a Dicer-1 mutant, we show that the stress-induced decrease in one of these genes, the insulator-binding factor Elba1, is dependent on functional miRNA biogenesis. Reduction in Elba1 correlates with the upregulation of early developmental genes and promotes a sustained weakening of heterochromatin in the adult fly as indicated by an increased expression of the PEV w^m4h reporter. We propose that maternal miRNAs, retained in response to an early embryonic heat shock, shape the subsequent de novo heterochromatin establishment that occurs during early development via direct or indirect regulation of some of the earliest expressed genes, including Elba1.

Off-target piRNA gene silencing in Drosophila melanogaster rescued by a transposable element insertion

Transposable elements (TE) are selfish genetic elements that can cause harmful mutations. In Drosophila, it has been estimated that half of all spontaneous visible marker phenotypes are mutations caused by TE insertions. Several factors likely limit the accumulation of exponentially amplifying TEs within genomes. First, synergistic interactions between TEs that amplify their harm with increasing copy number are proposed to limit TE copy number. However, the nature of this synergy is poorly understood. Second, because of the harm posed by TEs, eukaryotes have evolved systems of small RNA-based genome defense to limit transposition. However, as in all immune systems, there is a cost of autoimmunity and small RNA-based systems that silence TEs can inadvertently silence genes flanking TE insertions. In a screen for essential meiotic genes in Drosophila melanogaster, a truncated Doc retrotransposon within a neighboring gene was found to trigger the germline silencing of ald, the Drosophila Mps1 homolog, a gene essential for proper chromosome segregation in meiosis. A subsequent screen for suppressors of this silencing identified a new insertion of a Hobo DNA transposon in the same neighboring gene. Here we describe how the original Doc insertion triggers flanking piRNA biogenesis and local gene silencing. We show that this local gene silencing occurs in cis and is dependent on deadlock, a component of the Rhino-Deadlock-Cutoff (RDC) complex, to trigger dual-strand piRNA biogenesis at TE insertions. We further show how the additional Hobo insertion leads to de-silencing by reducing flanking piRNA biogenesis triggered by the original Doc insertion. These results support a model of TE-mediated gene silencing by piRNA biogenesis in cis that depends on local determinants of transcription. This may explain complex patterns of off-target gene silencing triggered by TEs within populations and in the laboratory. It also provides a mechanism of sign epistasis among TE insertions, illuminates the complex nature of their interactions and supports a model in which off-target gene silencing shapes the evolution of the RDC complex.

Genetic variation in P-element dysgenic sterility is associated with double-strand break repair and alternative splicing of TE transcripts

The germline mobilization of transposable elements (TEs) by small RNA mediated silencing pathways is conserved across eukaryotes and critical for ensuring the integrity of gamete genomes. However, genomes are recurrently invaded by novel TEs through horizontal transfer. These invading TEs are not targeted by host small RNAs, and their unregulated activity can cause DNA damage in germline cells and ultimately lead to sterility. Here we use hybrid dysgenesis-a sterility syndrome of Drosophila caused by transposition of invading P-element DNA transposons-to uncover host genetic variants that modulate dysgenic sterility. Using a panel of highly recombinant inbred lines of Drosophila melanogaster, we identified two linked quantitative trait loci (QTL) that determine the severity of dysgenic sterility in young and old females, respectively. We show that ovaries of fertile genotypes exhibit increased expression of splicing factors that suppress the production of transposase encoding transcripts, which likely reduces the transposition rate and associated DNA damage. We also show that fertile alleles are associated with decreased sensitivity to double-stranded breaks and enhanced DNA repair, explaining their ability to withstand high germline transposition rates. Together, our work reveals a diversity of mechanisms whereby host genotype modulates the cost of an invading TE, and points to genetic variants that were likely beneficial during the P-element invasion.

Invading viral DNA triggers dsRNA synthesis by RNA polymerase II to activate antiviral RNA interference in Drosophila

dsRNA sensing triggers antiviral responses against RNA and DNA viruses in diverse eukaryotes. In Drosophila, Invertebrate iridescent virus 6 (IIV-6), a large DNA virus, triggers production of small interfering RNAs (siRNAs) by the dsRNA sensor Dicer-2. Here, we show that host RNA polymerase II (RNAPII) bidirectionally transcribes specific AT-rich regions of the IIV-6 DNA genome to generate dsRNA. Both replicative and naked IIV-6 genomes trigger production of dsRNA in Drosophila cells, implying direct sensing of invading DNA. Loquacious-PD, a Dicer-2 co-factor essential for the biogenesis of endogenous siRNAs, is dispensable for processing of IIV-6-derived dsRNAs, which suggests that they are distinct. Consistent with this finding, inhibition of the RNAPII co-factor P-TEFb affects the synthesis of endogenous, but not virus-derived, dsRNA. Altogether, our results suggest that a non-canonical RNAPII complex recognizes invading viral DNA to synthesize virus-derived dsRNA, which activates the antiviral siRNA pathway in Drosophila.

Cloning, Sequencing, and Linkage Analysis of piRNAs

Piwi-interacting RNAs (piRNAs) are 25- to 32-nucleotide-long small RNAs that silence transposable elements (transposons) in animal gonads. piRNAs have a large sequence diversity (over one million different sequences per organism) to target a variety of transposon sequences. This is achieved by flexible and distinct biogenesis pathways that are evolutionarily conserved. In this chapter, I describe a detailed method of purifying and cloning piRNAs from freshly dissected tissue samples, such as fruit fly ovaries, for the high-throughput sequencing. I also describe how to computationally process the sequencing data and interrogate the characteristic pattern of piRNA biogenesis, including ping-pong amplification and head-to-tail phasing.

5´XP sRNA-seq: efficient identification of transcripts with and without 5´ phosphorylation reveals evolutionary conserved small RNA

Small RNA (sRNA) sequencing has been critical for our understanding of many cellular processes, including gene regulation. Nonetheless, the varying biochemical properties of sRNA, such as 5´ nucleotide modifications, make many sRNA subspecies incompatible with common protocols for sRNA sequencing. Here we describe 5XP-seq that outlines a novel strategy that captures a more complete picture of sRNA. By tagging 5´P sRNA during library preparation, 5XP-seq combines an open approach that includes all types of 5'-terminal modifications (5´X), with a selective approach for 5-phosphorylated sRNA (5´P). We show that 5XP-seq not only enriches phosphorylated miRNA and piRNA but successfully discriminates these sRNA from all other sRNA species. We further demonstrate the importance of this strategy by successful inter-species validation of sRNAs that would have otherwise failed, including human to insect translation of several tRNA (tRFs) and rRNA (rRFs) fragments. By combining 5´ insensitive library strategies with 5´ sensitive tagging, we have successfully tackled an intrinsic bias in modern sRNA sequencing that will help us reveal the true complexity and the evolutionary significance of the sRNA world.

Heterochromatin-dependent transcription of satellite DNAs in the Drosophila melanogaster female germline

Large blocks of tandemly repeated DNAs-satellite DNAs (satDNAs)-play important roles in heterochromatin formation and chromosome segregation. We know little about how satDNAs are regulated; however, their misregulation is associated with genomic instability and human diseases. We use the Drosophila melanogaster germline as a model to study the regulation of satDNA transcription and chromatin. Here we show that complex satDNAs (>100-bp repeat units) are transcribed into long noncoding RNAs and processed into piRNAs (PIWI interacting RNAs). This satDNA piRNA production depends on the Rhino-Deadlock-Cutoff complex and the transcription factor Moonshiner-a previously described non-canonical pathway that licenses heterochromatin-dependent transcription of dual-strand piRNA clusters. We show that this pathway is important for establishing heterochromatin at satDNAs. Therefore, satDNAs are regulated by piRNAs originating from their own genomic loci. This novel mechanism of satDNA regulation provides insight into the role of piRNA pathways in heterochromatin formation and genome stability.

Small RNA-Sequencing: Approaches and Considerations for miRNA Analysis

MicroRNAs (miRNAs) are a class of small RNA molecules that have an important regulatory role in multiple physiological and pathological processes. Their disease-specific profiles and presence in biofluids are properties that enable miRNAs to be employed as non-invasive biomarkers. In the past decades, several methods have been developed for miRNA analysis, including small RNA sequencing (RNA-seq). Small RNA-seq enables genome-wide profiling and analysis of known, as well as novel, miRNA variants. Moreover, its high sensitivity allows for profiling of low input samples such as liquid biopsies, which have now found applications in diagnostics and prognostics. Still, due to technical bias and the limited ability to capture the true miRNA representation, its potential remains unfulfilled. The introduction of many new small RNA-seq approaches that tried to minimize this bias, has led to the existence of the many small RNA-seq protocols seen today. Here, we review all current approaches to cDNA library construction used during the small RNA-seq workflow, with particular focus on their implementation in commercially available protocols. We provide an overview of each protocol and discuss their applicability. We also review recent benchmarking studies comparing each protocol's performance and summarize the major conclusions that can be gathered from their usage. The result documents variable performance of the protocols and highlights their different applications in miRNA research. Taken together, our review provides a comprehensive overview of all the current small RNA-seq approaches, summarizes their strengths and weaknesses, and provides guidelines for their applications in miRNA research.

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.

A kinesin Klp10A mediates cell cycle-dependent shuttling of Piwi between nucleus and nuage

The piRNA pathway protects germline genomes from selfish genetic elements (e.g. transposons) through their transcript cleavage in the cytoplasm and/or their transcriptional silencing in the nucleus. Here, we describe a mechanism by which the nuclear and cytoplasmic arms of the piRNA pathway are linked. We find that during mitosis of Drosophila spermatogonia, nuclear Piwi interacts with nuage, the compartment that mediates the cytoplasmic arm of the piRNA pathway. At the end of mitosis, Piwi leaves nuage to return to the nucleus. Dissociation of Piwi from nuage occurs at the depolymerizing microtubules of the central spindle, mediated by a microtubule-depolymerizing kinesin, Klp10A. Depletion of klp10A delays the return of Piwi to the nucleus and affects piRNA production, suggesting the role of nuclear-cytoplasmic communication in piRNA biogenesis. We propose that cell cycle-dependent communication between the nuclear and cytoplasmic arms of the piRNA pathway may play a previously unappreciated role in piRNA regulation.

The piRNA pathway that defends germline from selfish elements operates in two subpathways, one mediated by Piwi in Drosophila to silence transcription of targets in the nucleus and the other mediated by Aub and Ago3 to cleave transcripts of targets in the cytoplasm. How these two subpathways might coordinate with each other, particularly at the cell biological level, remains elusive. This study shows that Piwi interacts with Aub/Ago3 specifically in mitosis in nuage, the organelle that serves as the platform for piRNA cytoplasmic subpathway. Piwi returns to the nucleus at the end of mitosis, and our study suggests that dissociation of Piwi from nuage is facilitated by microtubule depolymerization by a kinesin Klp10A at the central spindle. We propose that cell-cycle-dependent interaction of two piRNA subpathways may play an important role in piRNA production.

Depletion of erythropoietic miR-486-5p and miR-451a improves detectability of rare microRNAs in peripheral blood-derived small RNA sequencing libraries

Erythroid-specific miR-451a and miR-486-5p are two of the most dominant microRNAs (miRNAs) in human peripheral blood. In small RNA sequencing libraries, their overabundance reduces diversity as well as complexity and consequently causes negative effects such as missing detectability and inaccurate quantification of low abundant miRNAs. Here we present a simple, cost-effective and easy to implement hybridization-based method to deplete these two erythropoietic miRNAs from blood-derived RNA samples. By utilization of blocking oligonucleotides, this method provides a highly efficient and specific depletion of miR-486-5p and miR-451a, which leads to a considerable increase of measured expression as well as detectability of low abundant miRNA species. The blocking oligos are compatible with common 5′ ligation-dependent small RNA library preparation protocols, including commercially available kits, such as Illumina TruSeq and Perkin Elmer NEXTflex. Furthermore, the here described method and oligo design principle can be easily adapted to target many other miRNA molecules, depending on context and research question.

CRISPR/Cas9-targeted removal of unwanted sequences from small-RNA sequencing libraries

In small RNA (smRNA) sequencing studies, highly abundant molecules such as adapter dimer products and tissue-specific microRNAs (miRNAs) inhibit accurate quantification of lowly expressed species. We previously developed a method to selectively deplete highly abundant miRNAs. However, this method does not deplete adapter dimer ligation products that, unless removed by gel-separation, comprise most of the library. Here, we have adapted and modified recently described methods for CRISPR/Cas9-based Depletion of Abundant Species by Hybridization ('DASH') to smRNA-seq, which we have termed miRNA and Adapter Dimer-DASH (MAD-DASH). In MAD-DASH, Cas9 is complexed with single guide RNAs (sgRNAs) targeting adapter dimer ligation products, alongside highly expressed tissue-specific smRNAs, for cleavage in vitro. This process dramatically reduces adapter dimer and targeted smRNA sequences, can be multiplexed, shows minimal off-target effects, improves the quantification of lowly expressed miRNAs from human plasma and tissue derived RNA, and obviates the need for gel-separation, greatly increasing sample throughput. Additionally, the method is fully customizable to other smRNA-seq preparation methods. Like depletion of ribosomal RNA for mRNA-seq and mitochondrial DNA for ATAC-seq, our method allows for greater proportional read-depth of non-targeted sequences.

The RNA-Binding ATPase, Armitage, Couples piRNA Amplification in Nuage to Phased piRNA Production on Mitochondria

PIWI-interacting RNAs (piRNAs) silence transposons in Drosophila ovaries, ensuring female fertility. Two coupled pathways generate germline piRNAs: the ping-pong cycle, in which the PIWI proteins Aubergine and Ago3 increase the abundance of pre-existing piRNAs, and the phased piRNA pathway, which generates strings of tail-to-head piRNAs, one after another. Proteins acting in the ping-pong cycle localize to nuage, whereas phased piRNA production requires Zucchini, an endonuclease on the mitochondrial surface. Here, we report that Armitage (Armi), an RNA-binding ATPase localized to both nuage and mitochondria, links the ping-pong cycle to the phased piRNA pathway. Mutations that block phased piRNA production deplete Armi from nuage. Armi ATPase mutants cannot support phased piRNA production and inappropriately bind mRNA instead of piRNA precursors. We propose that Armi shuttles between nuage and mitochondria, feeding precursor piRNAs generated by Ago3 cleavage into the Zucchini-dependent production of Aubergine- and Piwi-bound piRNAs on the mitochondrial surface.

A Transgenic Flock House Virus Replicon Reveals an RNAi Independent Antiviral Mechanism Acting in Drosophila Follicular Somatic Cells

The small interfering RNA (siRNA) pathway is the main and best studied invertebrate antiviral response. Other poorly characterized protein based antiviral mechanisms also contribute to the control of viral replication in insects. In addition, it remains unclear whether tissue specific factors contribute to RNA and protein-based antiviral immunity mechanisms. In vivo screens to identify such factors are challenging and time consuming. In addition, the scored phenotype is usually limited to survival and/or viral load. Transgenic viral replicons are valuable tools to overcome these limitations and screen for novel antiviral factors. Here we describe transgenic Drosophila melanogaster lines encoding a Flock House Virus-derived replicon (FHV∆B2eGFP), expressing GFP as a reporter of viral replication. This replicon is efficiently controlled by the siRNA pathway in most somatic tissues, with GFP fluorescence providing a reliable marker for the activity of antiviral RNAi. Interestingly, in follicular somatic cells (FSC) of ovaries, this replicon is still partially repressed in an siRNA independent manner. We did not detect replicon derived Piwi-interacting RNAs in FSCs and identified 31 differentially expressed genes between restrictive and permissive FSCs. Altogether, our results uncovered a yet unidentified RNAi-independent mechanism controlling FHV replication in FSCs of ovaries and validate the FHV∆B2eGFP replicon as a tool to screen for novel tissue specific antiviral mechanisms.

Maelstrom Represses Canonical Polymerase II Transcription within Bi-directional piRNA Clusters in Drosophila melanogaster

In Drosophila, 23-30 nt long PIWI-interacting RNAs (piRNAs) direct the protein Piwi to silence germline transposon transcription. Most germline piRNAs derive from dual-strand piRNA clusters, heterochromatic transposon graveyards that are transcribed from both genomic strands. These piRNA sources are marked by the heterochromatin protein 1 homolog Rhino (Rhi), which facilitates their promoter-independent transcription, suppresses splicing, and inhibits transcriptional termination. Here, we report that the protein Maelstrom (Mael) represses canonical, promoter-dependent transcription in dual-strand clusters, allowing Rhi to initiate piRNA precursor transcription. Mael also represses promoter-dependent transcription at sites outside clusters. At some loci, Mael repression requires the piRNA pathway, while at others, piRNAs play no role. We propose that by repressing canonical transcription of individual transposon mRNAs, Mael helps Rhi drive non-canonical transcription of piRNA precursors without generating mRNAs encoding transposon proteins.

Pathogenic tau-induced piRNA depletion promotes neuronal death through transposable element dysregulation in neurodegenerative tauopathies

Transposable elements, known colloquially as “jumping genes,” constitute approximately 45% of the human genome. Cells utilize epigenetic defenses to limit transposable element jumping, including formation of silencing heterochromatin and generation of piwi-interacting RNAs (piRNAs), small RNAs that facilitate clearance of transposable element transcripts. Here we identify transposable element dysregulation as a key mediator of neuronal death in tauopathies, a group of neurodegenerative disorders that are pathologically characterized by deposits of tau protein in the brain. Mechanistically, we find that heterochromatin decondensation and reduction of piwi/piRNAs drive transposable element dysregulation in tauopathy. We further report a significant increase in transcripts of the endogenous retrovirus class of transposable elements in human Alzheimer’s disease and progressive supranuclear palsy, suggesting that transposable element dysregulation is conserved in human tauopathy. Taken together, our data identify heterochromatin decondensation, piwi/piRNA depletion and consequent transposable element dysregulation as a novel, pharmacologically targetable, mechanistic driver of neurodegeneration in tauopathy.

Small RNA populations revealed by blocking rRNA fragments in Drosophila melanogaster reproductive tissues

RNA interference (RNAi) is a complex and highly conserved regulatory mechanism mediated via small RNAs (sRNAs). Recent technical advances in high throughput sequencing have enabled an increasingly detailed analysis of sRNA abundances and profiles in specific body parts and tissues. This enables investigations of the localized roles of microRNAs (miRNAs) and small interfering RNAs (siRNAs). However, variation in the proportions of non-coding RNAs in the samples being compared can hinder these analyses. Specific tissues may vary significantly in the proportions of fragments of longer non-coding RNAs (such as ribosomal RNA or transfer RNA) present, potentially reflecting tissue-specific differences in biological functions. For example, in Drosophila, some tissues contain a highly abundant 30nt rRNA fragment (the 2S rRNA) as well as abundant 5’ and 3’ terminal rRNA fragments. These can pose difficulties for the construction of sRNA libraries as they can swamp the sequencing space and obscure sRNA abundances. Here we addressed this problem and present a modified “rRNA blocking” protocol for the construction of high-definition (HD) adapter sRNA libraries, in D. melanogaster reproductive tissues. The results showed that 2S rRNAs targeted by blocking oligos were reduced from >80% to < 0.01% total reads. In addition, the use of multiple rRNA blocking oligos to bind the most abundant rRNA fragments allowed us to reveal the underlying sRNA populations at increased resolution. Side-by-side comparisons of sequencing libraries of blocked and non-blocked samples revealed that rRNA blocking did not change the miRNA populations present, but instead enhanced their abundances. We suggest that this rRNA blocking procedure offers the potential to improve the in-depth analysis of differentially expressed sRNAs within and across different tissues.

Pan-arthropod analysis reveals somatic piRNAs as an ancestral defence against transposable elements

In animals, small RNA molecules termed PIWI-interacting RNAs (piRNAs) silence transposable elements (TEs), protecting the germline from genomic instability and mutation. piRNAs have been detected in the soma in a few animals, but these are believed to be specific adaptations of individual species. Here, we report that somatic piRNAs were probably present in the ancestral arthropod more than 500 million years ago. Analysis of 20 species across the arthropod phylum suggests that somatic piRNAs targeting TEs and messenger RNAs are common among arthropods. The presence of an RNA-dependent RNA polymerase in chelicerates (horseshoe crabs, spiders and scorpions) suggests that arthropods originally used a plant-like RNA interference mechanism to silence TEs. Our results call into question the view that the ancestral role of the piRNA pathway was to protect the germline and demonstrate that small RNA silencing pathways have been repurposed for both somatic and germline functions throughout arthropod evolution.

Circulating Immune Cells Mediate a Systemic RNAi-Based Adaptive Antiviral Response in Drosophila

Effective antiviral protection in multicellular organisms relies on both cell-autonomous and systemic immunity. Systemic immunity mediates the spread of antiviral signals from infection sites to distant uninfected tissues. In arthropods, RNA interference (RNAi) is responsible for antiviral defense. Here, we show that flies have a sophisticated systemic RNAi-based immunity mediated by macrophage-like haemocytes. Haemocytes take up dsRNA from infected cells and, through endogenous transposon reverse transcriptases, produce virus-derived complementary DNAs (vDNA). These vDNAs template de novo synthesis of secondary viral siRNAs (vsRNA), which are secreted in exosome-like vesicles. Strikingly, exosomes containing vsRNAs, purified from haemolymph of infected flies, confer passive protection against virus challenge in naive animals. Thus, similar to vertebrates, insects use immune cells to generate immunological memory in the form of stable vDNAs that generate systemic immunity, which is mediated by the vsRNA-containing exosomes.

Co-chaperone Hsp70/Hsp90-organizing protein (Hop) is required for transposon silencing and Piwi-interacting RNA (piRNA) biogenesis

Piwi-interacting RNAs (piRNAs) are 26-30-nucleotide germ line-specific small non-coding RNAs that have evolutionarily conserved function in mobile genetic element (transposons) silencing and maintenance of genome integrity. Drosophila Hsp70/90-organizing protein homolog (Hop), a co-chaperone, interacts with piRNA-binding protein Piwi and mediates silencing of phenotypic variations. However, it is not known whether Hop has a direct role in piRNA biogenesis and transposon silencing. Here, we show that knockdown of Hop in the germ line nurse cells (GLKD) of Drosophila ovaries leads to activation of transposons. Hop GLKD females can lay eggs at the same rate as wild-type counterparts, but the eggs do not hatch into larvae. Hop GLKD leads to the accumulation of γ-H2Av foci in the germ line, indicating increased DNA damage in the ovary. We also show that Hop GLKD-induced transposon up-regulation is due to inefficient piRNA biogenesis. Based on these results, we conclude that Hop is a critical component of the piRNA pathway and that it maintains genome integrity by silencing transposons.

A somatic piRNA pathway in the Drosophila fat body ensures metabolic homeostasis and normal lifespan

In gonadal tissues, the Piwi-interacting (piRNA) pathway preserves genomic integrity by employing 23–29 nucleotide (nt) small RNAs complexed with argonaute proteins to suppress parasitic mobile sequences of DNA called transposable elements (TEs). Although recent evidence suggests that the piRNA pathway may be present in select somatic cells outside the gonads, the role of a non-gonadal somatic piRNA pathway is not well characterized. Here we report a functional somatic piRNA pathway in the adult Drosophila fat body including the presence of the piRNA effector protein Piwi and canonical 23–29 nt long TE-mapping piRNAs. The piwi mutants exhibit depletion of fat body piRNAs, increased TE mobilization, increased levels of DNA damage and reduced lipid stores. These mutants are starvation sensitive, immunologically compromised and short-lived, all phenotypes associated with compromised fat body function. These findings demonstrate the presence of a functional non-gonadal somatic piRNA pathway in the adult fat body that affects normal metabolism and overall organismal health.

The Piwi-interacting RNA (piRNA) pathway is known to suppress transposable elements in gonadal tissues. Here the authors provide evidence for a functional piRNA pathway in the somatic cells of the Drosophila fat body with roles in metabolism, immunological function and overall health.

Wolbachia Blocks Viral Genome Replication Early in Infection without a Transcriptional Response by the Endosymbiont or Host Small RNA Pathways

The intracellular endosymbiotic bacterium Wolbachia can protect insects against viral infection, and is being introduced into mosquito populations in the wild to block the transmission of arboviruses that infect humans and are a major public health concern. To investigate the mechanisms underlying this antiviral protection, we have developed a new model system combining Wolbachia-infected Drosophila melanogaster cell culture with the model mosquito-borne Semliki Forest virus (SFV; Togaviridae, Alphavirus). Wolbachia provides strong antiviral protection rapidly after infection, suggesting that an early stage post-infection is being blocked. Wolbachia does appear to have major effects on events distinct from entry, assembly or exit as it inhibits the replication of an SFV replicon transfected into the cells. Furthermore, it causes a far greater reduction in the expression of proteins from the 3´ open reading frame than the 5´ non-structural protein open reading frame, indicating that it is blocking the replication of viral RNA. Further to this separation of the replicase proteins and viral RNA in transreplication assays shows that uncoupling of viral RNA and replicase proteins does not overcome Wolbachia’s antiviral activity. This further suggests that replicative processes are disrupted, such as translation or replication, by Wolbachia infection. This may occur by Wolbachia mounting an active antiviral response, but the virus did not cause any transcriptional response by the bacterium, suggesting that this is not the case. Host microRNAs (miRNAs) have been implicated in protection, but again we found that host cell miRNA expression was unaffected by the bacterium and neither do our findings suggest any involvement of the antiviral siRNA pathway. We conclude that Wolbachia may directly interfere with early events in virus replication such as translation of incoming viral RNA or RNA transcription, and this likely involves an intrinsic (as opposed to an induced) mechanism.

The intracellular endosymbiotic bacterium Wolbachia can protect insects against viral infection. However, the mechanisms underlying this antiviral activity are poorly understood. We have developed a new model system combining Wolbachia-infected Drosophila melanogaster cell culture and the model mosquito-borne virus, Semliki Forest virus. Wolbachia confers strong antiviral activity against SFV. Our study indicates that viral replication appears to be inhibited at a very early stage, such as initial translation or replication. Results indicate that Wolbachia does not mount a transcriptional response to SFV infection and that host small RNA pathways are not involved in Wolbachia mediated antiviral activity in our system. We conclude that Wolbachia may directly interfere with early events in virus replication such as translation of incoming viral RNA or RNA transcription, and this likely involves an intrinsic (as opposed to an induced) mechanism.

RNA helicase Belle/DDX3 regulates transgene expression in Drosophila

Belle (Bel), the Drosophila homolog of the yeast DEAD-box RNA helicase DED1 and human DDX3, has been shown to be required for oogenesis and female fertility. Here we report a novel role of Bel in regulating the expression of transgenes. Abrogation of Bel by mutations or RNAi induces silencing of a variety of P-element-derived transgenes. This silencing effect depends on downregulation of their RNA levels. Our genetic studies have revealed that the RNA helicase Spindle-E (Spn-E), a nuage RNA helicase that plays a crucial role in regulating RNA processing and PIWI-interacting RNA (piRNA) biogenesis in germline cells, is required for loss-of-bel-induced transgene silencing. Conversely, Bel abrogation alleviates the nuage-protein mislocalization phenotype in spn-E mutants, suggesting a competitive relationship between these two RNA helicases. Additionally, disruption of the chromatin remodeling factor Mod(mdg4) or the microRNA biogenesis enzyme Dicer-1 (Dcr-1) also alleviates the transgene-silencing phenotypes in bel mutants, suggesting the involvement of chromatin remodeling and microRNA biogenesis in loss-of-bel-induced transgene silencing. Finally we show that genetic inhibition of Bel function leads to de novo generation of piRNAs from the transgene region inserted in the genome, suggesting a potential piRNA-dependent mechanism that may mediate transgene silencing as Bel function is inhibited.

The use of duplex-specific nuclease in ribosome profiling and a user-friendly software package for Ribo-seq data analysis

Ribosome profiling is a technique that permits genome-wide, quantitative analysis of translation and has found broad application in recent years. Here we describe a modified profiling protocol and software package designed to benefit more broadly the translation community in terms of simplicity and utility. The protocol, applicable to diverse organisms, including organelles, is based largely on previously published profiling methodologies, but uses duplex-specific nuclease (DSN) as a convenient, species-independent way to reduce rRNA contamination. We show that DSN-based depletion compares favorably with other commonly used rRNA depletion strategies and introduces little bias. The profiling protocol typically produces high levels of triplet periodicity, facilitating the detection of coding sequences, including upstream, downstream, and overlapping open reading frames (ORFs) and an alternative ribosome conformation evident during termination of protein synthesis. In addition, we provide a software package that presents a set of methods for parsing ribosomal profiling data from multiple samples, aligning reads to coding sequences, inferring alternative ORFs, and plotting average and transcript-specific aspects of the data. Methods are also provided for extracting the data in a form suitable for differential analysis of translation and translational efficiency.

Noncoding RNA. piRNA-guided transposon cleavage initiates Zucchini-dependent, phased piRNA production

PIWI-interacting RNAs (piRNAs) protect the animal germ line by silencing transposons. Primary piRNAs, generated from transcripts of genomic transposon "junkyards" (piRNA clusters), are amplified by the "ping-pong" pathway, yielding secondary piRNAs. We report that secondary piRNAs, bound to the PIWI protein Ago3, can initiate primary piRNA production from cleaved transposon RNAs. The first ~26 nucleotides (nt) of each cleaved RNA becomes a secondary piRNA, but the subsequent ~26 nt become the first in a series of phased primary piRNAs that bind Piwi, allowing piRNAs to spread beyond the site of RNA cleavage. The ping-pong pathway increases only the abundance of piRNAs, whereas production of phased primary piRNAs from cleaved transposon RNAs adds sequence diversity to the piRNA pool, allowing adaptation to changes in transposon sequence.

Probe-Directed Degradation (PDD) for Flexible Removal of Unwanted cDNA Sequences from RNA-Seq Libraries

Most applications for RNA-seq require the depletion of abundant transcripts to gain greater coverage of the underlying transcriptome. The sequences to be targeted for depletion depend on application and species and in many cases may not be supported by commercial depletion kits. This unit describes a method for generating RNA-seq libraries that incorporates probe-directed degradation (PDD), which can deplete any unwanted sequence set, with the low-bias split-adapter method of library generation (although many other library generation methods are in principle compatible). The overall strategy is suitable for applications requiring customized sequence depletion or where faithful representation of fragment ends and lack of sequence bias is paramount. We provide guidelines to rapidly design specific probes against the target sequence, and a detailed protocol for library generation using the split-adapter method including several strategies for streamlining the technique and reducing adapter dimer content.

Variation in piRNA and transposable element content in strains of Drosophila melanogaster

Transposable elements (TEs) are one of the most important features of genome architecture, so their evolution and relationship with host defense mechanisms have been topics of intense study, especially in model systems such as Drosophila melanogaster. Recently, a novel small RNA-based defense mechanism in animals called the Piwi-interacting RNA (piRNA) pathway was discovered to form an adaptive defense mechanism against TEs. To investigate the relationship between piRNA and TE content between strains of a species, we sequenced piRNAs from 16 inbred lines of D. melanogaster from the Drosophila Genetic Reference Panel. Instead of a global correlation of piRNA expression and TE content, we found evidence for a host response through de novo piRNA production from novel TE insertions. Although approximately 20% of novel TE insertions induced de novo piRNA production, the abundance of de novo piRNAs was low and did not markedly affect the global pool of ovarian piRNAs. Our results provide new insights into the evolution of TEs and the piRNA system in an important model organism.