Pitfalls of mapping high-throughput sequencing data to repetitive sequences: Piwi's genomic targets still not identified

Genome-wide distribution of 5-hydroxymethyluracil and chromatin accessibility in the Breviolum minutum genome

BACKGROUND

In dinoflagellates, a unique and extremely divergent genomic and nuclear organization has evolved. The highly unusual features of dinoflagellate nuclei and genomes include permanently condensed liquid crystalline chromosomes, primarily packaged by proteins other than histones, genes organized in very long unidirectional gene arrays, a general absence of transcriptional regulation, high abundance of the otherwise very rare DNA modification 5-hydroxymethyluracil (5-hmU), and many others. While most of these fascinating properties are originally identified in the 1970s and 1980s, they have not yet been investigated using modern genomic tools.

RESULTS

In this work, we address some of the outstanding questions regarding dinoflagellate genome organization by mapping the genome-wide distribution of 5-hmU (using both immunoprecipitation-based and basepair-resolution chemical mapping approaches) and of chromatin accessibility in the genome of the Symbiodiniaceae dinoflagellate Breviolum minutum. We find that the 5-hmU modification is preferentially enriched over certain classes of repetitive elements, often coincides with the boundaries between gene arrays, and is generally correlated with decreased chromatin accessibility, the latter otherwise being largely uniform along the genome. We discuss the potential roles of 5-hmU in the functional organization of dinoflagellate genomes and its relationship to the transcriptional landscape of gene arrays.

CONCLUSIONS

Our results provide the first window into the 5-hmU and chromatin accessibility landscapes in dinoflagellates.

AGO2 silences mobile transposons in the nucleus of quiescent cells

Argonaute 2 (AGO2) is a cytoplasmic component of the miRNA pathway, with essential roles in development and disease. Yet little is known about its regulation in vivo. Here we show that in quiescent mouse splenocytes, AGO2 localizes almost exclusively to the nucleus. AGO2 subcellular localization is modulated by the Pi3K-AKT-mTOR pathway, a well-established regulator of quiescence. Signaling through this pathway in proliferating cells promotes AGO2 cytoplasmic accumulation, at least in part by stimulating the expression of TNRC6, an essential AGO2 binding partner in the miRNA pathway. In quiescent cells in which mTOR signaling is low, AGO2 accumulates in the nucleus, where it binds to young mobile transposons co-transcriptionally to repress their expression via its catalytic domain. Our data point to an essential but previously unrecognized nuclear role for AGO2 during quiescence as part of a genome-defense system against young mobile elements and provide evidence of RNA interference in the soma of mammals.

The chromatin landscape of the euryarchaeon Haloferax volcanii

BACKGROUND

Archaea, together with Bacteria, represent the two main divisions of life on Earth, with many of the defining characteristics of the more complex eukaryotes tracing their origin to evolutionary innovations first made in their archaeal ancestors. One of the most notable such features is nucleosomal chromatin, although archaeal histones and chromatin differ significantly from those of eukaryotes, not all archaea possess histones and it is not clear if histones are a main packaging component for all that do. Despite increased interest in archaeal chromatin in recent years, its properties have been little studied using genomic tools.

RESULTS

Here, we adapt the ATAC-seq assay to archaea and use it to map the accessible landscape of the genome of the euryarchaeote Haloferax volcanii. We integrate the resulting datasets with genome-wide maps of active transcription and single-stranded DNA (ssDNA) and find that while H. volcanii promoters exist in a preferentially accessible state, unlike most eukaryotes, modulation of transcriptional activity is not associated with changes in promoter accessibility. Applying orthogonal single-molecule footprinting methods, we quantify the absolute levels of physical protection of H. volcanii and find that Haloferax chromatin is similarly or only slightly more accessible, in aggregate, than that of eukaryotes. We also evaluate the degree of coordination of transcription within archaeal operons and make the unexpected observation that some CRISPR arrays are associated with highly prevalent ssDNA structures.

CONCLUSIONS

Our results provide the first comprehensive maps of chromatin accessibility and active transcription in Haloferax across conditions and thus a foundation for future functional studies of archaeal chromatin.

Genome-wide distribution of 5-hydroxymethyluracil and chromatin accessibility in the Breviolum minutum genome

In dinoflagellates, a unique and extremely divergent genomic and nuclear organization has evolved. The highly unusual features of dinoflagellate nuclei and genomes include permanently condensed liquid crystalline chromosomes, primarily packaged by proteins other than histones, genes organized in very long unidirectional gene arrays, a general absence of transcriptional regulation, high abundance of the otherwise very rare DNA modification 5-hydroxymethyluracil (5-hmU), and many others. While most of these fascinating properties were originally identified in the 1970s and 1980s, they have not yet been investigated using modern genomic tools. In this work, we address some of the outstanding questions regarding dinoflagellate genome organization by mapping the genome-wide distribution of 5-hmU (using both immunoprecipitation-based and basepair-resolution chemical mapping approaches) and of chromatin accessibility in the genome of the Symbiodiniaceae dinoflagellate Breviolum minutum. We find that the 5-hmU modification is preferentially enriched over certain classes of repetitive elements, often coincides with the boundaries between gene arrays, and is generally correlated with decreased chromatin accessibility, the latter otherwise being largely uniform along the genome. We discuss the potential roles of 5-hmU in the functional organization of dinoflagellate genomes and its relationship to the transcriptional landscape of gene arrays.

T3E: a tool for characterising the epigenetic profile of transposable elements using ChIP-seq data

BACKGROUND

Despite the advent of Chromatin Immunoprecipitation Sequencing (ChIP-seq) having revolutionised our understanding of the mammalian genome's regulatory landscape, many challenges remain. In particular, because of their repetitive nature, the sequencing reads derived from transposable elements (TEs) pose a real bioinformatics challenge, to the point that standard analysis pipelines typically ignore reads whose genomic origin cannot be unambiguously ascertained.

RESULTS

We show that discarding ambiguously mapping reads may lead to a systematic underestimation of the number of reads associated with young TE families/subfamilies. We also provide evidence suggesting that the strategy of randomly permuting the location of the read mappings (or the TEs) that is often used to compute the background for enrichment calculations at TE families/subfamilies can result in both false positive and negative enrichments. To address these problems, we present the Transposable Element Enrichment Estimator (T3E), a tool that makes use of ChIP-seq data to characterise the epigenetic profile of associated TE families/subfamilies. T3E weights the number of read mappings assigned to the individual TE copies of a family/subfamily by the overall number of genomic loci to which the corresponding reads map, and this is done at the single nucleotide level. In addition, T3E computes ChIP-seq enrichment relative to a background estimated based on the distribution of the read mappings in the input control DNA. We demonstrated the capabilities of T3E on 23 different ChIP-seq libraries. T3E identified enrichments that were consistent with previous studies. Furthermore, T3E detected context-specific enrichments that are likely to pinpoint unexplored TE families/subfamilies with individual TE copies that have been frequently exapted as cis-regulatory elements during the evolution of mammalian regulatory networks.

CONCLUSIONS

T3E is a novel open-source computational tool (available for use at: https://github.com/michelleapaz/T3E ) that overcomes some of the pitfalls associated with the analysis of ChIP-seq data arising from the repetitive mammalian genome and provides a framework to shed light on the epigenetics of entire TE families/subfamilies.

Chromatin structure can introduce systematic biases in genome-wide analyses of Plasmodium falciparum

Background: The maintenance, regulation, and dynamics of heterochromatin in the human malaria parasite, Plasmodium falciparum, has drawn increasing attention due to its regulatory role in mutually exclusive virulence gene expression and the silencing of key developmental regulators. The advent of genome-wide analyses such as chromatin-immunoprecipitation followed by sequencing (ChIP-seq) has been instrumental in understanding chromatin composition; however, even in model organisms, ChIP-seq experiments are susceptible to intrinsic experimental biases arising from underlying chromatin structure. Methods: We performed a control ChIP-seq experiment, re-analyzed previously published ChIP-seq datasets and compared different analysis approaches to characterize biases of genome-wide analyses in P. falciparum. Results: We found that heterochromatic regions in input control samples used for ChIP-seq normalization are systematically underrepresented in regard to sequencing coverage across the P. falciparum genome. This underrepresentation, in combination with a non-specific or inefficient immunoprecipitation, can lead to the identification of false enrichment and peaks across these regions. We observed that such biases can also be seen at background levels in specific and efficient ChIP-seq experiments. We further report on how different read mapping approaches can also skew sequencing coverage within highly similar subtelomeric regions and virulence gene families. To ameliorate these issues, we discuss orthogonal methods that can be used to characterize bona fide chromatin-associated proteins. Conclusions: Our results highlight the impact of chromatin structure on genome-wide analyses in the parasite and the need for caution when characterizing chromatin-associated proteins and features.

Stonewall prevents expression of ectopic genes in the ovary and accumulates at insulator elements in D. melanogaster

Germline stem cells (GSCs) are the progenitor cells of the germline for the lifetime of an animal. In Drosophila, these cells reside in a cellular niche that is required for both their maintenance (self-renewal) and differentiation (asymmetric division resulting in a daughter cell that differs from the GSC). The stem cell—daughter cell transition is tightly regulated by a number of processes, including an array of proteins required for genome stability. The germline stem-cell maintenance factor Stonewall (Stwl) associates with heterochromatin, but its molecular function is poorly understood. We performed RNA-Seq on stwl mutant ovaries and found significant derepression of many transposon families but not heterochromatic genes. We also discovered inappropriate expression of multiple classes of genes. Most prominent are testis-enriched genes, including the male germline sex-determination switch Phf7, the differentiation factor bgcn, and a large testis-specific gene cluster on chromosome 2, all of which are upregulated or ectopically expressed in stwl mutant ovaries. Surprisingly, we also found that RNAi knockdown of stwl in somatic S2 cells results in ectopic expression of these testis genes. Using parallel ChIP-Seq and RNA-Seq experiments in S2 cells, we discovered that Stwl localizes upstream of transcription start sites and at heterochromatic sequences including repetitive sequences associated with telomeres. Stwl is also enriched at bgcn, suggesting that it directly regulates this essential differentiation factor. Finally, we identify Stwl binding motifs that are shared with known insulator binding proteins. We propose that Stwl affects gene regulation, including repression of male transcripts in the female germline, by binding insulators and establishing chromatin boundaries.

Stem cells are defined by their ability to divide asymmetrically, resulting in a differentiated cell and a stem cell daughter. In fruit flies, sperm and egg production begins with germline stem cells (GSCs). The ability of a GSC to differentiate or self-renew is tightly regulated by a myriad of factors. Some of these are transcription factors, which are responsible for activating or suppressing other genes to promote one state in favor of another. Stonewall is an ovarian nuclear protein required for GSC self-renewal, whose molecular function is poorly understood. Here we show that Stonewall is responsible for preventing the activation of “male” molecular programming in the fruit fly ovary. When Stonewall is absent from the ovary, egg production is terminated and testis-specific genes become highly expressed, including the male transcript of Phf7, which induces male sexual identity in female germ cells. We also show that Stonewall is likely localizing to genomic insulators, which are regions of the genome that shield genes from nearby regulators. Our findings suggest that Stonewall helps to organize the genome in ovarian germ cells and prevent expression of male genes.

The chromatin organization of a chlorarachniophyte nucleomorph genome

BACKGROUND

Nucleomorphs are remnants of secondary endosymbiotic events between two eukaryote cells wherein the endosymbiont has retained its eukaryotic nucleus. Nucleomorphs have evolved at least twice independently, in chlorarachniophytes and cryptophytes, yet they have converged on a remarkably similar genomic architecture, characterized by the most extreme compression and miniaturization among all known eukaryotic genomes. Previous computational studies have suggested that nucleomorph chromatin likely exhibits a number of divergent features.

RESULTS

In this work, we provide the first maps of open chromatin, active transcription, and three-dimensional organization for the nucleomorph genome of the chlorarachniophyte Bigelowiella natans. We find that the B. natans nucleomorph genome exists in a highly accessible state, akin to that of ribosomal DNA in some other eukaryotes, and that it is highly transcribed over its entire length, with few signs of polymerase pausing at transcription start sites (TSSs). At the same time, most nucleomorph TSSs show very strong nucleosome positioning. Chromosome conformation (Hi-C) maps reveal that nucleomorph chromosomes interact with one other at their telomeric regions and show the relative contact frequencies between the multiple genomic compartments of distinct origin that B. natans cells contain.

CONCLUSIONS

We provide the first study of a nucleomorph genome using modern functional genomic tools, and derive numerous novel insights into the physical and functional organization of these unique genomes.

Transposable elements and Piwi‑interacting RNAs in hemato‑oncology with a focus on myelodysplastic syndrome (Review)

Our current understanding of hematopoietic stem cell differentiation and the abnormalities that lead to leukemogenesis originates from the accumulation of knowledge regarding protein‑coding genes. However, the possible impact of transposable element (TE) mobilization and the expression of P‑element‑induced WImpy testis‑interacting RNAs (piRNAs) on leukemogenesis has been beyond the scope of scientific interest to date. The expression profiles of these molecules and their importance for human health have only been characterized recently due to the rapid progress of high‑throughput sequencing technology development. In the present review, current knowledge on the expression profile and function of TEs and piRNAs was summarized, with specific focus on their reported involvement in leukemogenesis and pathogenesis of myelodysplastic syndrome.

piRNA- and siRNA-mediated transcriptional repression in Drosophila, mice, and yeast: new insights and biodiversity

The PIWI‐interacting RNA (piRNA) pathway acts as a self‐defense mechanism against transposons to maintain germline genome integrity. Failures in the piRNA pathway cause DNA damage in the germline genome, disturbing inheritance of “correct” genetic information by the next generations and leading to infertility. piRNAs execute transposon repression in two ways: degrading their RNA transcripts and compacting the genomic loci via heterochromatinization. The former event is mechanistically similar to siRNA‐mediated RNA cleavage that occurs in the cytoplasm and has been investigated in many species including nematodes, fruit flies, and mammals. The latter event seems to be mechanistically parallel to siRNA‐centered kinetochore assembly and subsequent chromosome segregation, which has so far been studied particularly in fission yeast. Despite the interspecies conservations, the overall schemes of the nuclear events show clear biodiversity across species. In this review, we summarize the recent progress regarding piRNA‐mediated transcriptional silencing in Drosophila and discuss the biodiversity by comparing it with the equivalent piRNA‐mediated system in mice and the siRNA‐mediated system in fission yeast.

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.

A dual role of dLsd1 in oogenesis: regulating developmental genes and repressing transposons

The histone demethylase LSD1 is a key chromatin regulator that is often deregulated in cancer. Its ortholog, dLsd1 plays a crucial role in Drosophila oogenesis; however, our knowledge of dLsd1 function is insufficient to explain its role in the ovary. Here, we have performed genome-wide analysis of dLsd1 binding in the ovary, and we document that dLsd1 is preferentially associated to the transcription start site of developmental genes. We uncovered an unanticipated interplay between dLsd1 and the GATA transcription factor Serpent and we report an unexpected role for Serpent in oogenesis. Besides, our transcriptomic data show that reducing dLsd1 levels results in ectopic transposable elements (TE) expression correlated with changes in H3K4me2 and H3K9me2 at TE loci. In addition, our results suggest that dLsd1 is required for Piwi dependent TE silencing. Hence, we propose that dLsd1 plays crucial roles in establishing specific gene expression programs and in repressing transposons during oogenesis.

Tools and best practices for retrotransposon analysis using high-throughput sequencing data

Background

Sequencing technologies give access to a precise picture of the molecular mechanisms acting upon genome regulation. One of the biggest technical challenges with sequencing data is to map millions of reads to a reference genome. This problem is exacerbated when dealing with repetitive sequences such as transposable elements that occupy half of the mammalian genome mass. Sequenced reads coming from these regions introduce ambiguities in the mapping step. Therefore, applying dedicated parameters and algorithms has to be taken into consideration when transposable elements regulation is investigated with sequencing datasets.

Results

Here, we used simulated reads on the mouse and human genomes to define the best parameters for aligning transposable element-derived reads on a reference genome. The efficiency of the most commonly used aligners was compared and we further evaluated how transposable element representation should be estimated using available methods. The mappability of the different transposon families in the mouse and the human genomes was calculated giving an overview into their evolution.

Conclusions

Based on simulated data, we provided recommendations on the alignment and the quantification steps to be performed when transposon expression or regulation is studied, and identified the limits in detecting specific young transposon families of the mouse and human genomes. These principles may help the community to adopt standard procedures and raise awareness of the difficulties encountered in the study of transposable elements.

piRNA-guided co-transcriptional silencing coopts nuclear export factors

The PIWI-interacting RNA (piRNA) pathway is a small RNA-based immune system that controls the expression of transposons and maintains genome integrity in animal gonads. In Drosophila, piRNA-guided silencing is achieved, in part, via co-transcriptional repression of transposons by Piwi. This depends on Panoramix (Panx); however, precisely how an RNA binding event silences transcription remains to be determined. Here we show that Nuclear Export Factor 2 (Nxf2) and its co-factor, Nxt1, form a complex with Panx and are required for co-transcriptional silencing of transposons in somatic and germline cells of the ovary. Tethering of Nxf2 or Nxt1 to RNA results in silencing of target loci and the concomitant accumulation of repressive chromatin marks. Nxf2 and Panx proteins are mutually required for proper localization and stability. We mapped the protein domains crucial for the Nxf2/Panx complex formation and show that the amino-terminal portion of Panx is sufficient to induce transcriptional silencing.

Ro60 and Y RNAs: structure, functions, and roles in autoimmunity

Ro60, also known as SS-A or TROVE2, is an evolutionarily conserved RNA-binding protein that is found in most animal cells, approximately 5% of sequenced prokaryotic genomes and some archaea. Ro60 is present in cells as both a free protein and as a component of a ribonucleoprotein complex, where its best-known partners are members of a class of noncoding RNAs called Y RNAs. Structural and biochemical analyses have revealed that Ro60 is a ring-shaped protein that binds Y RNAs on its outer surface. In addition to Y RNAs, Ro60 binds misfolded and aberrant noncoding RNAs in some animal cell nuclei. Although the fate of these defective Ro60-bound noncoding RNAs in animal cells is not well-defined, a bacterial Ro60 ortholog functions with 3' to 5' exoribonucleases to assist structured RNA degradation. Studies of Y RNAs have revealed that these RNAs regulate the subcellular localization of Ro60, tether Ro60 to effector proteins and regulate the access of other RNAs to its central cavity. As both mammalian cells and bacteria lacking Ro60 are sensitized to ultraviolet irradiation, Ro60 function may be important during exposure to some environmental stressors. Here we summarize the current knowledge regarding the functions of Ro60 and Y RNAs in animal cells and bacteria. Because the Ro60 RNP is a clinically important target of autoantibodies in patients with rheumatic diseases such as Sjogren's syndrome, systemic lupus erythematosus, and neonatal lupus, we also discuss potential roles for Ro60 RNPs in the initiation and pathogenesis of systemic autoimmune rheumatic disease.

Ensemble learning method for the prediction of new bioactive molecules

Pharmacologically active molecules can provide remedies for a range of different illnesses and infections. Therefore, the search for such bioactive molecules has been an enduring mission. As such, there is a need to employ a more suitable, reliable, and robust classification method for enhancing the prediction of the existence of new bioactive molecules. In this paper, we adopt a recently developed combination of different boosting methods (Adaboost) for the prediction of new bioactive molecules. We conducted the research experiments utilizing the widely used MDL Drug Data Report (MDDR) database. The proposed boosting method generated better results than other machine learning methods. This finding suggests that the method is suitable for inclusion among the in silico tools for use in cheminformatics, computational chemistry and molecular biology.

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.

Alternative Computational Analysis Shows No Evidence for Nucleosome Enrichment at Repetitive Sequences in Mammalian Spermatozoa

Samans et al. (2014) reported the enrichment of nucleosomes in human and bovine spermatozoa at centromere repeats and retrotransposon sequences such as LINE-1 and SINE. We demonstrate here that nucleosomal enrichments at repetitive sequences as reported result from bioinformatic analyses that make redundant use of sequencing reads that map to multiple locations in the genome. To illustrate that this computational approach is flawed, we observed comparable artificial enrichments at repetitive sequences when aligning control genomic DNA or simulated reads of uniform genome coverage. These results imply that the main conclusions of the article by Samans et al. (2014) are confounded by an inappropriate computational methodology used to analyze the primary data.

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.

Selectively Constrained RNA Editing Regulation Crosstalks with piRNA Biogenesis in Primates

Although millions of RNA editing events have been reported to modify hereditary information across the primate transcriptome, evidence for their functional significance remains largely elusive, particularly for the vast majority of editing sites in noncoding regions. Here, we report a new mechanism for the functionality of RNA editing—a crosstalk with PIWI-interacting RNA (piRNA) biogenesis. Exploiting rhesus macaque as an emerging model organism closely related to human, in combination with extensive genome and transcriptome sequencing in seven tissues of the same animal, we deciphered accurate RNA editome across both long transcripts and the piRNA species. Superimposing and comparing these two distinct RNA editome profiles revealed 4,170 editing-bearing piRNA variants, or epiRNAs, that primarily derived from edited long transcripts. These epiRNAs represent distinct entities that evidence an intersection between RNA editing regulations and piRNA biogenesis. Population genetics analyses in a macaque population of 31 independent animals further demonstrated that the epiRNA-associated RNA editing is maintained by purifying selection, lending support to the functional significance of this crosstalk in rhesus macaque. Correspondingly, these findings are consistent in human, supporting the conservation of this mechanism during the primate evolution. Overall, our study reports the earliest lines of evidence for a crosstalk between selectively constrained RNA editing regulation and piRNA biogenesis, and further illustrates that such an interaction may contribute substantially to the diversification of the piRNA repertoire in primates.

Reassessment of Piwi binding to the genome and Piwi impact on RNA polymerase II distribution

Drosophila Piwi was reported by Huang et al. (2013) to be guided by piRNAs to piRNA-complementary sites in the genome, which then recruits heterochromatin protein 1a and histone methyltransferase Su(Var)3-9 to the sites. Among additional findings, Huang et al. (2013) also reported Piwi binding sites in the genome and the reduction of RNA polymerase II in euchromatin but its increase in pericentric regions in piwi mutants. Marinov et al. (2015) disputed the validity of the Huang et al. bioinformatic pipeline that led to the last two claims. Here we report our independent reanalysis of the data using current bioinformatic methods. Our reanalysis agrees with Marinov et al. (2015) that Piwi's genomic targets still remain to be identified but confirms the Huang et al. claim that Piwi influences RNA polymerase II distribution in the genome. This Matters Arising Response addresses the Marinov et al. (2015) Matters Arising, published concurrently in this issue of Developmental Cell.