TRIBE: Hijacking an RNA-Editing Enzyme to Identify Cell-Specific Targets of RNA-Binding Proteins.

RNA transcripts are bound and regulated by RNA-binding proteins (RBPs). Current methods for identifying in vivo targets of an RBP are imperfect and not amenable to examining small numbers of cells. To address these issues, we developed TRIBE (targets of RNA-binding proteins identified by editing), a technique that couples an RBP to the catalytic domain of the Drosophila RNA-editing enzyme ADAR and expresses the fusion protein in vivo. RBP targets are marked with novel RNA editing events and identified by sequencing RNA. We have used TRIBE to identify the targets of three RBPs (Hrp48, dFMR1, and NonA). TRIBE compares favorably to other methods, including CLIP, and we have identified RBP targets from as little as 150 specific fly neurons. TRIBE can be performed without an antibody and in small numbers of specific cells.

Transposable element landscapes in aging Drosophila.

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

TDP-43 dysfunction restricts dendritic complexity by inhibiting CREB activation and altering gene expression.

Amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) are two related neurodegenerative diseases that present with similar TDP-43 pathology in patient tissue. TDP-43 is an RNA-binding protein which forms aggregates in neurons of ALS and FTD patients as well as in a subset of patients diagnosed with other neurodegenerative diseases. Despite our understanding that TDP-43 is essential for many aspects of RNA metabolism, it remains obscure how TDP-43 dysfunction contributes to neurodegeneration. Interestingly, altered neuronal dendritic morphology is a common theme among several neurological disorders and is thought to precede neurodegeneration. We previously found that both TDP-43 overexpression (OE) and knockdown (KD) result in reduced dendritic branching of cortical neurons. In this study, we used TRIBE (targets of RNA-binding proteins identified by editing) as an approach to identify signaling pathways that regulate dendritic branching downstream of TDP-43. We found that TDP-43 RNA targets are enriched for pathways that signal to the CREB transcription factor. We further found that TDP-43 dysfunction inhibits CREB activation and CREB transcriptional output, and restoring CREB signaling rescues defects in dendritic branching. Finally, we demonstrate, using RNA sequencing, that TDP-43 OE and KD cause similar changes in the abundance of specific messenger RNAs, consistent with their ability to produce similar morphological defects. Our data therefore provide a mechanism by which TDP-43 dysfunction interferes with dendritic branching, and may define pathways for therapeutic intervention in neurodegenerative diseases.

Mechanistic implications of enhanced editing by a HyperTRIBE RNA-binding protein.

We previously developed TRIBE, a method for the identification of cell-specific RNA-binding protein targets. TRIBE expresses an RBP of interest fused to the catalytic domain (cd) of the RNA-editing enzyme ADAR and performs adenosine-to-inosine editing on RNA targets of the RBP. However, target identification is limited by the low editing efficiency of the ADARcd. Here we describe HyperTRIBE, which carries a previously characterized hyperactive mutation (E488Q) of the ADARcd. HyperTRIBE identifies dramatically more editing sites, many of which are also edited by TRIBE but at a much lower editing frequency. HyperTRIBE therefore more faithfully recapitulates the known binding specificity of its RBP than TRIBE. In addition, separating RNA binding from the enhanced editing activity of the HyperTRIBE ADAR catalytic domain sheds light on the mechanism of ADARcd editing as well as the enhanced activity of the HyperADARcd.

RNA-seq analysis of Drosophila clock and non-clock neurons reveals neuron-specific cycling and novel candidate neuropeptides.

Locomotor activity rhythms are controlled by a network of ~150 circadian neurons within the adult Drosophila brain. They are subdivided based on their anatomical locations and properties. We profiled transcripts "around the clock" from three key groups of circadian neurons with different functions. We also profiled a non-circadian outgroup, dopaminergic (TH) neurons. They have cycling transcripts but fewer than clock neurons as well as low expression and poor cycling of clock gene transcripts. This suggests that TH neurons do not have a canonical circadian clock and that their gene expression cycling is driven by brain systemic cues. The three circadian groups are surprisingly diverse in their cycling transcripts and overall gene expression patterns, which include known and putative novel neuropeptides. Even the overall phase distributions of cycling transcripts are distinct, indicating that different regulatory principles govern transcript oscillations. This surprising cell-type diversity parallels the functional heterogeneity of the different neurons.

Conserved piRNA Expression from a Distinct Set of piRNA Cluster Loci in Eutherian Mammals.

The Piwi pathway is deeply conserved amongst animals because one of its essential functions is to repress transposons. However, many Piwi-interacting RNAs (piRNAs) do not base-pair to transposons and remain mysterious in their targeting function. The sheer number of piRNA cluster (piC) loci in animal genomes and infrequent piRNA sequence conservation also present challenges in determining which piC loci are most important for development. To address this question, we determined the piRNA expression patterns of piC loci across a wide phylogenetic spectrum of animals, and reveal that most genic and intergenic piC loci evolve rapidly in their capacity to generate piRNAs, regardless of known transposon silencing function. Surprisingly, we also uncovered a distinct set of piC loci with piRNA expression conserved deeply in Eutherian mammals. We name these loci Eutherian-Conserved piRNA cluster (ECpiC) loci. Supporting the hypothesis that conservation of piRNA expression across ~100 million years of Eutherian evolution implies function, we determined that one ECpiC locus generates abundant piRNAs antisense to the STOX1 transcript, a gene clinically associated with preeclampsia. Furthermore, we confirmed reduced piRNAs in existing mouse mutations at ECpiC-Asb1 and -Cbl, which also display spermatogenic defects. The Asb1 mutant testes with strongly reduced Asb1 piRNAs also exhibit up-regulated gene expression profiles. These data indicate ECpiC loci may be specially adapted to support Eutherian reproduction.

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

Piwi proteins and Piwi-interacting RNAs (piRNAs) repress transposable elements (TEs) from mobilizing in gonadal cells. To determine the spectrum of piRNA-regulated targets that may extend beyond TEs, we conducted a genome-wide survey for transcripts associated with PIWI and for transcripts affected by PIWI knockdown in Drosophila ovarian somatic sheet (OSS) cells, a follicle cell line expressing the Piwi pathway. Despite the immense sequence diversity among OSS cell piRNAs, our analysis indicates that TE transcripts are the major transcripts associated with and directly regulated by PIWI. However, several coding genes were indirectly regulated by PIWI via an adjacent de novo TE insertion that generated a nascent TE transcript. Interestingly, we noticed that PIWI-regulated genes in OSS cells greatly differed from genes affected in a related follicle cell culture, ovarian somatic cells (OSCs). Therefore, we characterized the distinct genomic TE insertions across four OSS and OSC lines and discovered dynamic TE landscapes in gonadal cultures that were defined by a subset of active TEs. Particular de novo TEs appeared to stimulate the expression of novel candidate long noncoding RNAs (lncRNAs) in a cell lineage-specific manner, and some of these TE-associated lncRNAs were associated with PIWI and overlapped PIWI-regulated genes. Our analyses of OSCs and OSS cells demonstrate that despite having a Piwi pathway to suppress endogenous mobile elements, gonadal cell TE landscapes can still dramatically change and create transcriptome diversity.

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.

Identifying complete RNA structural ensembles including pseudoknots.

The close relationship between RNA structure and function underlines the significance of accurately predicting RNA structures from sequence information. Structural topologies such as pseudoknots are of particular interest due to their ubiquity and direct involvement in RNA function, but identifying pseudoknots is a computationally challenging problem and existing heuristic approaches usually perform poorly for RNA sequences of even a few hundred bases. We survey the performance of pseudoknot prediction methods on a data set of full-length RNA sequences representing varied sequence lengths, and biological RNA classes such as RNase P RNA, Group I Intron, tmRNA and tRNA. Pseudoknot prediction methods are compared with minimum free energy and suboptimal secondary structure prediction methods in terms of correct base-pairs, stems and pseudoknots and we find that the ensemble of suboptimal structure predictions succeeds in identifying correct structural elements in RNA that are usually missed in MFE and pseudoknot predictions. We propose a strategy to identify a comprehensive set of non-redundant stems in the suboptimal structure space of a RNA molecule by applying heuristics that reduce the structural redundancy of the predicted suboptimal structures by merging slightly varying stems that are predicted to form in local sequence regions. This reduced-redundancy set of structural elements consistently outperforms more specialized approaches.in data sets. Thus, the suboptimal folding space can be used to represent the structural diversity of an RNA molecule more comprehensively than optimal structure prediction approaches alone.

Protocol for using TRIBE to study RNA-protein interactions and nuclear organization in mammalian cells.

Targets of RNA-binding proteins discovered by editing (TRIBE) determines RNA-proteins interactions and nuclear organization with minimal false positives. We detail necessary steps for performing mammalian cell RBP-TRIBE to determine the targets of RNA-binding proteins and MS2-TRIBE to determine RNA-RNA interactions within the nucleus. Necessary steps for performing a TRIBE experiment are detailed, starting with plasmid/cell line generation, cellular transfection, and RNA sequencing library preparation and concluding with bioinformatics analysis of RNA editing sites and identification of target RNAs. For complete details on the use and execution of this protocol, please refer to Biswas et al. (2020).

MS2-TRIBE Evaluates Both Protein-RNA Interactions and Nuclear Organization of Transcription by RNA Editing.

Both UV-cross-linking and immunoprecipitation (CLIP) and RNA editing (TRIBE) can identify the targets of RNA-binding proteins. To evaluate false-positives of CLIP and TRIBE, endogenous ß-actin mRNA was tagged with MS2 stem loops, making it the only bona fide target mRNA for the MS2 capsid protein (MCP). CLIP and TRIBE detected ß-actin, albeit with false-positives. False-positive CLIP signals were attributed to nonspecific antibody interactions. In contrast, putative false-positive TRIBE targets were genes spatially proximal to the ß-actin gene. MCP-ADAR edited nearby nascent transcripts consistent with interchromosomal contacts observed in Hi-C. The identification of nascent contacts implies RNA regulatory proteins (e.g., splicing factors) associated with multiple nascent transcripts, forming domains of post-transcriptional activity. Repeating these results with an integrated inducible MS2 reporter indicated that MS2-TRIBE can be applied to a broad array of cells and transcripts to study spatial organization and nuclear RNA regulation.

TRIBE editing reveals specific mRNA targets of eIF4E-BP in Drosophila and in mammals.

4E-BP (eIF4E-BP) represses translation initiation by binding to the 5' cap-binding protein eIF4E and inhibiting its activity. Although 4E-BP has been shown to be important in growth control, stress response, cancer, neuronal activity, and mammalian circadian rhythms, it is not understood how it preferentially represses a subset of mRNAs. We successfully used HyperTRIBE (targets of RNA binding proteins identified by editing) to identify in vivo 4E-BP mRNA targets in both Drosophila and mammals under conditions known to activate 4E-BP. The protein associates with specific mRNAs, and ribosome profiling data show that mTOR inhibition changes the translational efficiency of 4E-BP TRIBE targets more substantially compared to nontargets. In both systems, these targets have specific motifs and are enriched in translation-related pathways, which correlate well with the known activity of 4E-BP and suggest that it modulates the binding specificity of eIF4E and contributes to mTOR translational specificity.

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

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

Identification of RNA-binding protein targets with HyperTRIBE.

RNA-binding proteins (RBPs) accompany RNA from birth to death, affecting RNA biogenesis and functions. Identifying RBP-RNA interactions is essential to understanding their complex roles in different cellular processes. However, detecting in vivo RNA targets of RBPs, especially in a small number of discrete cells, has been a technically challenging task. We previously developed a novel technique called TRIBE (targets of RNA-binding proteins identified by editing) to overcome this problem. TRIBE expresses a fusion protein consisting of a queried RBP and the catalytic domain of the RNA-editing enzyme ADAR (adenosine deaminase acting on RNA) (ADARcd), which marks target RNA transcripts by converting adenosine to inosine near the RBP binding sites. These marks can be subsequently identified via high-throughput sequencing. In spite of its usefulness, TRIBE is constrained by a low editing efficiency and editing-sequence bias from the ADARcd. Therefore, we developed HyperTRIBE by incorporating a previously characterized hyperactive mutation, E488Q, into the ADARcd. This strategy increases the editing efficiency and reduces sequence bias, which markedly increases the sensitivity of this technique without sacrificing specificity. HyperTRIBE provides a more powerful strategy for identifying RNA targets of RBPs with an easy experimental and computational protocol at low cost, that can be performed not only in flies, but also in mammals. The HyperTRIBE experimental protocol described below can be carried out in cultured Drosophila S2 cells in 1 week, using tools available in a common molecular biology laboratory; the computational analysis requires 3 more days.

Drosophila PAF1 Modulates PIWI/piRNA Silencing Capacity.

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

Genome-wide identification of neuronal activity-regulated genes in Drosophila.

Activity-regulated genes (ARGs) are important for neuronal functions like long-term memory and are well-characterized in mammals but poorly studied in other model organisms like Drosophila. Here we stimulated fly neurons with different paradigms and identified ARGs using high-throughput sequencing from brains as well as from sorted neurons: they included a narrow set of circadian neurons as well as dopaminergic neurons. Surprisingly, many ARGs are specific to the stimulation paradigm and very specific to neuron type. In addition and unlike mammalian immediate early genes (IEGs), fly ARGs do not have short gene lengths and are less enriched for transcription factor function. Chromatin assays using ATAC-sequencing show that the transcription start sites (TSS) of ARGs do not change with neural firing but are already accessible prior to stimulation. Lastly based on binding site enrichment in ARGs, we identified transcription factor mediators of firing and created neuronal activity reporters.