The piRNA Response to Retroviral Invasion of the Koala Genome.

Antisense Piwi-interacting RNAs (piRNAs) guide silencing of established transposons during germline development, and sense piRNAs drive ping-pong amplification of the antisense pool, but how the germline responds to genome invasion is not understood. The KoRV-A gammaretrovirus infects the soma and germline and is sweeping through wild koalas by a combination of horizontal and vertical transfer, allowing direct analysis of retroviral invasion of the germline genome. Gammaretroviruses produce spliced Env mRNAs and unspliced transcripts encoding Gag, Pol, and the viral genome, but KoRV-A piRNAs are almost exclusively derived from unspliced genomic transcripts and are strongly sense-strand biased. Significantly, selective piRNA processing of unspliced proviral transcripts is conserved from insects to placental mammals. We speculate that bypassed splicing generates a conserved molecular pattern that directs proviral genomic transcripts to the piRNA biogenesis machinery and that this "innate" piRNA response suppresses transposition until antisense piRNAs are produced, establishing sequence-specific adaptive immunity.

The HP1 homolog rhino anchors a nuclear complex that suppresses piRNA precursor splicing.

piRNAs guide an adaptive genome defense system that silences transposons during germline development. The Drosophila HP1 homolog Rhino is required for germline piRNA production. We show that Rhino binds specifically to the heterochromatic clusters that produce piRNA precursors, and that binding directly correlates with piRNA production. Rhino colocalizes to germline nuclear foci with Rai1/DXO-related protein Cuff and the DEAD box protein UAP56, which are also required for germline piRNA production. RNA sequencing indicates that most cluster transcripts are not spliced and that rhino, cuff, and uap56 mutations increase expression of spliced cluster transcripts over 100-fold. LacI::Rhino fusion protein binding suppresses splicing of a reporter transgene and is sufficient to trigger piRNA production from a trans combination of sense and antisense reporters. We therefore propose that Rhino anchors a nuclear complex that suppresses cluster transcript splicing and speculate that stalled splicing differentiates piRNA precursors from mRNAs.

UAP56 couples piRNA clusters to the perinuclear transposon silencing machinery.

piRNAs silence transposons during germline development. In Drosophila, transcripts from heterochromatic clusters are processed into primary piRNAs in the perinuclear nuage. The nuclear DEAD box protein UAP56 has been previously implicated in mRNA splicing and export, whereas the DEAD box protein Vasa has an established role in piRNA production and localizes to nuage with the piRNA binding PIWI proteins Ago3 and Aub. We show that UAP56 colocalizes with the cluster-associated HP1 variant Rhino, that nuage granules containing Vasa localize directly across the nuclear envelope from cluster foci containing UAP56 and Rhino, and that cluster transcripts immunoprecipitate with both Vasa and UAP56. Significantly, a charge-substitution mutation that alters a conserved surface residue in UAP56 disrupts colocalization with Rhino, germline piRNA production, transposon silencing, and perinuclear localization of Vasa. We therefore propose that UAP56 and Vasa function in a piRNA-processing compartment that spans the nuclear envelope.

Adaptation to P element transposon invasion in Drosophila melanogaster.

Transposons evolve rapidly and can mobilize and trigger genetic instability. Piwi-interacting RNAs (piRNAs) silence these genome pathogens, but it is unclear how the piRNA pathway adapts to invasion of new transposons. In Drosophila, piRNAs are encoded by heterochromatic clusters and maternally deposited in the embryo. Paternally inherited P element transposons thus escape silencing and trigger a hybrid sterility syndrome termed P-M hybrid dysgenesis. We show that P-M hybrid dysgenesis activates both P elements and resident transposons and disrupts the piRNA biogenesis machinery. As dysgenic hybrids age, however, fertility is restored, P elements are silenced, and P element piRNAs are produced de novo. In addition, the piRNA biogenesis machinery assembles, and resident elements are silenced. Significantly, resident transposons insert into piRNA clusters, and these new insertions are transmitted to progeny, produce novel piRNAs, and are associated with reduced transposition. P element invasion thus triggers heritable changes in genome structure that appear to enhance transposon silencing.

Epigenetic and chromosomal features drive transposon insertion in Drosophila melanogaster.

Transposons are mobile genetic elements prevalent in the genomes of most species. The distribution of transposons within a genome reflects the actions of two opposing processes: initial insertion site selection, and selective pressure from the host. By analyzing whole-genome sequencing data from transposon-activated Drosophila melanogaster, we identified 43 316 de novo and 237 germline insertions from four long-terminal-repeat (LTR) transposons, one LINE transposon (I-element), and one DNA transposon (P-element). We found that all transposon types favored insertion into promoters de novo, but otherwise displayed distinct insertion patterns. De novo and germline P-element insertions preferred replication origins, often landing in a narrow region around transcription start sites and in regions of high chromatin accessibility. De novo LTR transposon insertions preferred regions with high H3K36me3, promoters and exons of active genes; within genes, LTR insertion frequency correlated with gene expression. De novo I-element insertion density increased with distance from the centromere. Germline I-element and LTR transposon insertions were depleted in promoters and exons, suggesting strong selective pressure to remove transposons from functional elements. Transposon movement is associated with genome evolution and disease; therefore, our results can improve our understanding of genome and disease biology.

The Drosophila HP1 homolog Rhino is required for transposon silencing and piRNA production by dual-strand clusters.

Piwi-interacting RNAs (piRNAs) silence transposons and maintain genome integrity during germline development. In Drosophila, transposon-rich heterochromatic clusters encode piRNAs either on both genomic strands (dual-strand clusters) or predominantly one genomic strand (uni-strand clusters). Primary piRNAs derived from these clusters are proposed to drive a ping-pong amplification cycle catalyzed by proteins that localize to the perinuclear nuage. We show that the HP1 homolog Rhino is required for nuage organization, transposon silencing, and ping-pong amplification of piRNAs. rhi mutations virtually eliminate piRNAs from the dual-strand clusters and block production of putative precursor RNAs from both strands of the major 42AB dual-strand cluster, but not of transcripts or piRNAs from the uni-strand clusters. Furthermore, Rhino protein associates with the 42AB dual-strand cluster,but does not bind to uni-strand cluster 2 or flamenco. Rhino thus appears to promote transcription of dual-strand clusters, leading to production of piRNAs that drive the ping-pong amplification cycle.

Collapse of germline piRNAs in the absence of Argonaute3 reveals somatic piRNAs in flies.

Piwi-interacting RNAs (piRNAs) silence transposons in animal germ cells. piRNAs are thought to derive from long transcripts spanning transposon-rich genomic loci and to direct an autoamplification loop in which an antisense piRNA, bound to Aubergine or Piwi protein, triggers production of a sense piRNA bound to the PIWI protein Argonaute3 (Ago3). In turn, the new piRNA is envisioned to produce a second antisense piRNA. Here, we describe strong loss-of-function mutations in ago3, allowing a direct genetic test of this model. We find that Ago3 acts to amplify piRNA pools and to enforce on them an antisense bias, increasing the number of piRNAs that can act to silence transposons. We also detect a second, Ago3-independent piRNA pathway centered on Piwi. Transposons targeted by this second pathway often reside in the flamenco locus, which is expressed in somatic ovarian follicle cells, suggesting a role for piRNAs beyond the germline.

A benchmark and an algorithm for detecting germline transposon insertions and measuring de novo transposon insertion frequencies.

Transposons are genomic parasites, and their new insertions can cause instability and spur the evolution of their host genomes. Rapid accumulation of short-read whole-genome sequencing data provides a great opportunity for studying new transposon insertions and their impacts on the host genome. Although many algorithms are available for detecting transposon insertions, the task remains challenging and existing tools are not designed for identifying de novo insertions. Here, we present a new benchmark fly dataset based on PacBio long-read sequencing and a new method TEMP2 for detecting germline insertions and measuring de novo 'singleton' insertion frequencies in eukaryotic genomes. TEMP2 achieves high sensitivity and precision for detecting germline insertions when compared with existing tools using both simulated data in fly and experimental data in fly and human. Furthermore, TEMP2 can accurately assess the frequencies of de novo transposon insertions even with high levels of chimeric reads in simulated datasets; such chimeric reads often occur during the construction of short-read sequencing libraries. By applying TEMP2 to published data on hybrid dysgenic flies inflicted by de-repressed P-elements, we confirmed the continuous new insertions of P-elements in dysgenic offspring before they regain piRNAs for P-element repression. TEMP2 is freely available at Github: https://github.com/weng-lab/TEMP2.

High-resolution analysis of differential gene expression during skeletal muscle atrophy and programmed cell death.

Skeletal muscles can undergo atrophy and/or programmed cell death (PCD) during development or in response to a wide range of insults, including immobility, cachexia, and spinal cord injury. However, the protracted nature of atrophy and the presence of multiple cell types within the tissue complicate molecular analyses. One model that does not suffer from these limitations is the intersegmental muscle (ISM) of the tobacco hawkmoth Manduca sexta. Three days before the adult eclosion (emergence) at the end of metamorphosis, the ISMs initiate a nonpathological program of atrophy that results in a 40% loss of mass. The ISMs then generate the eclosion behavior and initiate a nonapoptotic PCD during the next 30 h. We have performed a comprehensive transcriptomics analysis of all mRNAs and microRNAs throughout ISM development to better understand the molecular mechanisms that mediate atrophy and death. Atrophy involves enhanced protein catabolism and reduced expression of the genes involved in respiration, adhesion, and the contractile apparatus. In contrast, PCD involves the induction of numerous proteases, DNA methylases, membrane transporters, ribosomes, and anaerobic metabolism. These changes in gene expression are largely repressed when insects are injected with the insect steroid hormone 20-hydroxyecdysone, which delays death. The expression of the death-associated proteins may be greatly enhanced by reductions in specific microRNAs that function to repress translation. This study not only provides fundamental new insights into basic developmental processes, it may also represent a powerful resource for identifying potential diagnostic markers and molecular targets for therapeutic intervention.

Structural insights into Rhino-Deadlock complex for germline piRNA cluster specification.

PIWI-interacting RNAs (piRNAs) silence transposons in germ cells to maintain genome stability and animal fertility. Rhino, a rapidly evolving heterochromatin protein 1 (HP1) family protein, binds Deadlock in a species-specific manner and so defines the piRNA-producing loci in the Drosophila genome. Here, we determine the crystal structures of Rhino-Deadlock complex in Drosophila melanogaster and simulans In both species, one Rhino binds the N-terminal helix-hairpin-helix motif of one Deadlock protein through a novel interface formed by the beta-sheet in the Rhino chromoshadow domain. Disrupting the interface leads to infertility and transposon hyperactivation in flies. Our structural and functional experiments indicate that electrostatic repulsion at the interaction interface causes cross-species incompatibility between the sibling species. By determining the molecular architecture of this piRNA-producing machinery, we discover a novel HP1-partner interacting mode that is crucial to piRNA biogenesis and transposon silencing. We thus explain the cross-species incompatibility of two sibling species at the molecular level.

A systems level approach to temporal expression dynamics in Drosophila reveals clusters of long term memory genes.

The ability to integrate experiential information and recall it in the form of memory is observed in a wide range of taxa, and is a hallmark of highly derived nervous systems. Storage of past experiences is critical for adaptive behaviors that anticipate both adverse and positive environmental factors. The process of memory formation and consolidation involve many synchronized biological events including gene transcription, protein modification, and intracellular trafficking: However, many of these molecular mechanisms remain illusive. With Drosophila as a model system we use a nonassociative memory paradigm and a systems level approach to uncover novel transcriptional patterns. RNA sequencing of Drosophila heads during and after memory formation identified a number of novel memory genes. Tracking the dynamic expression of these genes over time revealed complex gene networks involved in long term memory. In particular, this study focuses on two functional gene clusters of signal peptides and proteases. Bioinformatics network analysis and prediction in combination with high-throughput RNA sequencing identified previously unknown memory genes, which when genetically knocked down resulted in behaviorally validated memory defects.

Heterotypic piRNA Ping-Pong requires qin, a protein with both E3 ligase and Tudor domains.

piRNAs guide PIWI proteins to silence transposons in animal germ cells. Reciprocal cycles of piRNA-directed RNA cleavage--catalyzed by the PIWI proteins Aubergine (Aub) and Argonaute3 (Ago3) in Drosophila melanogaster--expand the population of antisense piRNAs in response to transposon expression, a process called the Ping-Pong cycle. Heterotypic Ping-Pong between Aub and Ago3 ensures that antisense piRNAs predominate. We show that qin, a piRNA pathway gene whose protein product contains both E3 ligase and Tudor domains, colocalizes with Aub and Ago3 in nuage, a perinuclear structure implicated in transposon silencing. In qin mutants, less Ago3 binds Aub, futile Aub:Aub homotypic Ping-Pong prevails, antisense piRNAs decrease, many families of mobile genetic elements are reactivated, and DNA damage accumulates in nurse cells and oocytes. We propose that Qin enforces heterotypic Ping-Pong between Aub and Ago3, ensuring that transposons are silenced and maintaining the integrity of the germline genome.

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

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

Aub, Vasa and Armi localization to phase separated nuage is dispensable for piRNA biogenesis and transposon silencing in Drosophila.

From nematodes to placental mammals, key components of the germline transposon silencing piRNAs pathway localize to phase separated perinuclear granules. In Drosophila, the PIWI protein Aub, DEAD box protein Vasa and helicase Armi localize to nuage granules and are required for ping-pong piRNA amplification and phased piRNA processing. Drosophila piRNA mutants lead to genome instability and Chk2 kinase DNA damage signaling. By systematically analyzing piRNA pathway organization, small RNA production, and long RNA expression in single piRNA mutants and corresponding chk2/mnk double mutants, we show that Chk2 activation disrupts nuage localization of Aub and Vasa, and that the HP1 homolog Rhino, which drives piRNA precursor transcription, is required for Aub, Vasa, and Armi localization to nuage. However, these studies also show that ping-pong amplification and phased piRNA biogenesis are independent of nuage localization of Vasa, Aub and Armi. Dispersed cytoplasmic proteins thus appear to mediate these essential piRNA pathway functions.

Drosophila rasiRNA pathway mutations disrupt embryonic axis specification through activation of an ATR/Chk2 DNA damage response.

Small repeat-associated siRNAs (rasiRNAs) mediate silencing of retrotransposons and the Stellate locus. Mutations in the Drosophila rasiRNA pathway genes armitage and aubergine disrupt embryonic axis specification, triggering defects in microtubule polarization as well as asymmetric localization of mRNA and protein determinants in the developing oocyte. Mutations in the ATR/Chk2 DNA damage signal transduction pathway dramatically suppress these axis specification defects, but do not restore retrotransposon or Stellate silencing. Furthermore, rasiRNA pathway mutations lead to germline-specific accumulation of gamma-H2Av foci characteristic of DNA damage. We conclude that rasiRNA-based gene silencing is not required for axis specification, and that the critical developmental function for this pathway is to suppress DNA damage signaling in the germline.

Kinesin I-dependent cortical exclusion restricts pole plasm to the oocyte posterior.

Microtubules and the plus-end-directed microtubule motor Kinesin I are required for the selective accumulation of oskar mRNA at the posterior cortex of the Drosophila melanogaster oocyte, which is essential to posterior patterning and pole plasm assembly. We present evidence that microtubule minus ends associate with the entire cortex, and that Kinesin and microtubules are not required for oskar mRNA association with the posterior pole, but prevent ectopic localization of this transcript and the pole plasm proteins Oskar and Vasa to other cortical regions. Cortical binding of oskar mRNA seems to be dependent on the actin cytoskeleton. We conclude that most of the actin-rich oocyte cortex can support pole plasm assembly, and propose that Kinesin restricts pole plasm formation to the posterior by moving oskar mRNA away from microtubule-rich lateral and anterior cortical regions.

piRNA-independent transposon silencing by the Drosophila THO complex.

piRNAs guide Piwi/Panoramix-dependent H3K9me3 chromatin modification and transposon silencing during Drosophila germline development. The THO RNA export complex is composed of Hpr1, Tho2, and Thoc5-7. Null thoc7 mutations, which displace Thoc5 and Thoc6 from a Tho2-Hpr1 subcomplex, reduce expression of a subset of germline piRNAs and increase transposon expression, suggesting that THO silences transposons by promoting piRNA biogenesis. Here, we show that the thoc7-null mutant combination increases transposon transcription but does not reduce anti-sense piRNAs targeting half of the transcriptionally activated transposon families. These mutations also fail to reduce piRNA-guided H3K9me3 chromatin modification or block Panoramix-dependent silencing of a reporter transgene, and unspliced transposon transcripts co-precipitate with THO through a Piwi- and Panoramix-independent mechanism. Mutations in piwi also dominantly enhance germline defects associated with thoc7-null alleles. THO thus functions in a piRNA-independent transposon-silencing pathway, which acts cooperatively with Piwi to support germline development.

Long first exons and epigenetic marks distinguish conserved pachytene piRNA clusters from other mammalian genes.

In the male germ cells of placental mammals, 26-30-nt-long PIWI-interacting RNAs (piRNAs) emerge when spermatocytes enter the pachytene phase of meiosis. In mice, pachytene piRNAs derive from ~100 discrete autosomal loci that produce canonical RNA polymerase II transcripts. These piRNA clusters bear 5' caps and 3' poly(A) tails, and often contain introns that are removed before nuclear export and processing into piRNAs. What marks pachytene piRNA clusters to produce piRNAs, and what confines their expression to the germline? We report that an unusually long first exon (≥ 10 kb) or a long, unspliced transcript correlates with germline-specific transcription and piRNA production. Our integrative analysis of transcriptome, piRNA, and epigenome datasets across multiple species reveals that a long first exon is an evolutionarily conserved feature of pachytene piRNA clusters. Furthermore, a highly methylated promoter, often containing a low or intermediate level of CG dinucleotides, correlates with germline expression and somatic silencing of pachytene piRNA clusters. Pachytene piRNA precursor transcripts bind THOC1 and THOC2, THO complex subunits known to promote transcriptional elongation and mRNA nuclear export. Together, these features may explain why the major sources of pachytene piRNA clusters specifically generate these unique small RNAs in the male germline of placental mammals.

Somatic piRNAs and Transposons are Differentially Expressed Coincident with Skeletal Muscle Atrophy and Programmed Cell Death.

PIWI-interacting RNAs (piRNAs) are small single-stranded RNAs that can repress transposon expression via epigenetic silencing and transcript degradation. They have been identified predominantly in the ovary and testis, where they serve essential roles in transposon silencing in order to protect the integrity of the genome in the germline. The potential expression of piRNAs in somatic cells has been controversial. In the present study we demonstrate the expression of piRNAs derived from both genic and transposon RNAs in the intersegmental muscles (ISMs) from the tobacco hawkmoth Manduca sexta. These piRNAs are abundantly expressed, ∼27 nt long, map antisense to transposons, are oxidation resistant, exhibit a 5' uridine bias, and amplify via the canonical ping-pong pathway. An RNA-seq analysis demonstrated that 19 piRNA pathway genes are expressed in the ISMs and are developmentally regulated. The abundance of piRNAs does not change when the muscles initiate developmentally-regulated atrophy, but are repressed coincident with the commitment of the muscles undergo programmed cell death at the end of metamorphosis. This change in piRNA expression is correlated with the repression of several retrotransposons and the induction of specific DNA transposons. The developmentally-regulated changes in the expression of piRNAs, piRNA pathway genes, and transposons are all regulated by 20-hydroxyecdysone, the steroid hormone that controls the timing of ISM death. Taken together, these data provide compelling evidence for the existence of piRNA in somatic tissues and suggest that they may play roles in developmental processes such as programmed cell death.