Genome-wide assessment of post-transcriptional control in the fly brain

Post-transcriptional control of gene expression has central importance during development and adulthood and in physiology in general. However, little is known about the extent of post-transcriptional control of gene expression in the brain. Most post-transcriptional regulatory effectors (e.g., miRNAs) destabilize target mRNAs by shortening their polyA tails. Hence, the fraction of a given mRNA that it is fully polyadenylated should correlate with its stability and serves as a good measure of post-transcriptional control. Here, we compared RNA-seq datasets from fly brains that were generated either from total (rRNA-depleted) or polyA-selected RNA. By doing this comparison we were able to compute a coefficient that measures the extent of post-transcriptional control for each brain-expressed mRNA. In agreement with current knowledge, we found that mRNAs encoding ribosomal proteins, metabolic enzymes, and housekeeping genes are among the transcripts with least post-transcriptional control, whereas mRNAs that are known to be highly unstable, like circadian mRNAs and mRNAs expressing synaptic proteins and proteins with neuronal functions, are under strong post-transcriptional control. Surprisingly, the latter group included many specific groups of genes relevant to brain function and behavior. In order to determine the importance of miRNAs in this regulation, we profiled miRNAs from fly brains using oligonucleotide microarrays. Surprisingly, we did not find a strong correlation between the expression levels of miRNAs in the brain and the stability of their target mRNAs; however, genes identified as highly regulated post-transcriptionally were strongly enriched for miRNA targets. This demonstrates a central role of miRNAs for modulating the levels and turnover of brain-specific mRNAs in the fly

Drosophila PSI controls circadian period and the phase of circadian behavior under temperature cycle via tim splicing

The Drosophila circadian pacemaker consists of transcriptional feedback loops subjected to post-transcriptional and post-translational regulation. While post-translational regulatory mechanisms have been studied in detail, much less is known about circadian post-transcriptional control. Thus, we targeted 364 RNA binding and RNA associated proteins with RNA interference. Among the 43 hits we identified was the alternative splicing regulator P-element somatic inhibitor (PSI). PSI regulates the thermosensitive alternative splicing of timeless (tim), promoting splicing events favored at warm temperature over those increased at cold temperature. Psi downregulation shortens the period of circadian rhythms and advances the phase of circadian behavior under temperature cycle. Interestingly, both phenotypes were suppressed in flies that could produce TIM proteins only from a transgene that cannot form the thermosensitive splicing isoforms. Therefore, we conclude that PSI regulates the period of Drosophila circadian rhythms and circadian behavior phase during temperature cycling through its modulation of the tim splicing pattern

Defining the 5 and 3 landscape of the Drosophila transcriptome with Exo-seq and RNaseH-seq

Cells regulate biological responses in part through changes in transcription start sites (TSS) or cleavage and polyadenylation sites (PAS). To fully understand gene regulatory networks, it is therefore critical to accurately annotate cell type-specific TSS and PAS. Here we present a simple and straightforward approach for genome-wide annotation of 5- and 3-RNA ends. Our approach reliably discerns bona fide PAS from false PAS that arise due to internal poly(A) tracts, a common problem with current PAS annotation methods. We applied our methodology to study the impact of temperature on the Drosophila melanogaster head transcriptome. We found hundreds of previously unidentified TSS and PAS which revealed two interesting phenomena: first, genes with multiple PASs tend to harbor a motif near the most proximal PAS, which likely represents a new cleavage and polyadenylation signal. Second, motif analysis of promoters of genes affected by temperature suggested that boundary element association factor of 32 kDa (BEAF-32) and DREF mediates a transcriptional program at warm temperatures, a result we validated in a fly line where beaf-32 is downregulated. These results demonstrate the utility of a high-throughput platform for complete experimental and computational analysis of mRNA-ends to improve gene annotation

An in vivo strategy for knockdown of circular RNAs

Exonic circular RNAs (circRNAs) are highly abundant RNAs generated mostly from exons of protein-coding genes. Assaying the functions of circRNAs is not straightforward as common approaches for circRNA depletion tend to also alter the levels of mRNAs generated from the hosting gene. Here we describe a methodology for specific knockdown of circRNAs in vivo with tissue and cell resolution. We also describe an experimental and computational platform for determining the potential off-target effects as well as for verifying the obtained phenotypes. Briefly, we utilize shRNAs targeted to the circRNA-specific back-splice junction to specifically downregulate the circRNA. We utilized this methodology to downregulate five circRNAs that are highly expressed in Drosophila. There were no effects on the levels of their linear counterparts or any RNA with complementarity to the expressed shRNA. Interestingly, downregulation of circCtrip resulted in developmental lethality that was recapitulated with a second shRNA. Moreover, downregulation of individual circRNAs caused specific changes in the fly head transcriptome, suggesting roles for these circRNAs in the fly nervous system. Together, our results provide a methodological approach that enables the comprehensive study of circRNAs at the organismal and cellular levels and generated for the first time flies in which specific circRNAs are downregulated

End Sequence Analysis Toolkit (ESAT) expands the extractable information from single-cell RNA-seq data

RNA-seq protocols that focus on transcript termini are well suited for applications in which template quantity is limiting. Here we show that, when applied to end-sequencing data, analytical methods designed for global RNA-seq produce computational artifacts. To remedy this, we created the End Sequence Analysis Toolkit (ESAT). As a test, we first compared end-sequencing and bulk RNA-seq using RNA from dendritic cells stimulated with lipopolysaccharide (LPS). As predicted by the telescripting model for transcriptional bursts, ESAT detected an LPS-stimulated shift to shorter 3\u27-isoforms that was not evident by conventional computational methods. Then, droplet-based microfluidics was used to generate 1000 cDNA libraries, each from an individual pancreatic islet cell. ESAT identified nine distinct cell types, three distinct beta-cell types, and a complex interplay between hormone secretion and vascularization. ESAT, then, offers a much-needed and generally applicable computational pipeline for either bulk or single-cell RNA end-sequencing

Translation of circRNAs

Circular RNAs (circRNAs) are abundant and evolutionarily conserved RNAs of largely unknown function. Here, we show that a subset of circRNAs is translated in vivo. By performing ribosome footprinting from fly heads, we demonstrate that a group of circRNAs is associated with translating ribosomes. Many of these ribo-circRNAs use the start codon of the hosting mRNA, are bound by membrane-associated ribosomes, and have evolutionarily conserved termination codons. In addition, we found that a circRNA generated from the muscleblind locus encodes a protein, which we detected in fly head extracts by mass spectrometry. Next, by performing in vivo and in vitro translation assays, we show that UTRs of ribo-circRNAs (cUTRs) allow cap-independent translation. Moreover, we found that starvation and FOXO likely regulate the translation of a circMbl isoform. Altogether, our study provides strong evidence for translation of circRNAs, revealing the existence of an unexplored layer of gene activity

The Molecular Clockwork of the Fire Ant Solenopsis invicta

This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication

Regulatory mechanisms for 3′-end alternative splicing and polyadenylation of the Glial Fibrillary Acidic Protein, GFAP, transcript

The glial fibrillary acidic protein, GFAP, forms the intermediate cytoskeleton in cells of the glial lineage. Besides the common GFAPα transcript, the GFAPε and GFAPκ transcripts are generated by alternative mRNA 3′-end processing. Here we use a GFAP minigene to characterize molecular mechanisms participating in alternative GFAP expression. Usage of a polyadenylation signal within the alternatively spliced exon 7a is essential to generate the GFAPκ and GFAPκ transcripts. The GFAPκ mRNA is distinct from GFAPε mRNA given that it also includes intron 7a. Polyadenylation at the exon 7a site is stimulated by the upstream splice site. Moreover, exon 7a splice enhancer motifs supported both exon 7a splicing and polyadenylation. SR proteins increased the usage of the exon 7a polyadenylation signal but not the exon 7a splicing, whereas the polypyrimidine tract binding (PTB) protein enhanced both exon 7a polyadenylation and exon 7a splicing. Finally, increasing transcription by the VP16 trans-activator did not affect the frequency of use of the exon 7a polyadenylation signal whereas the exon 7a splicing frequency was decreased. Our data suggest a model with the selection of the exon 7a polyadenylation site being the essential and primary event for regulating GFAP alternative processing

Dynamics of the Drosophila Circadian Clock: Theoretical Anti-Jitter Network and Controlled Chaos

Background: Electronic clocks exhibit undesirable jitter or time variations in periodic signals. The circadian clocks of humans, some animals, and plants consist of oscillating molecular networks with peak-to-peak time of approximately 24 hours. Clockwork orange (CWO) is a transcriptional repressor of Drosophila direct target genes. Methodology/Principal Findings: Theory and data from a model of the Drosophila circadian clock support the idea that CWO controls anti-jitter negative circuits that stabilize peak-to-peak time in light-dark cycles (LD). The orbit is confined to chaotic attractors in both LD and dark cycles and is almost periodic in LD; furthermore, CWO diminishes the Euclidean dimension of the chaotic attractor in LD. Light resets the clock each day by restricting each molecular peak to the proximity of a prescribed time. Conclusions/Significance: The theoretical results suggest that chaos plays a central role in the dynamics of the Drosophila circadian clock and that a single molecule, CWO, may sense jitter and repress it by its negative loops

Circadian Transcription Contributes to Core Period Determination in Drosophila

The Clock–Cycle (CLK–CYC) heterodimer constitutes a key circadian transcription complex in Drosophila. CYC has a DNA-binding domain but lacks an activation domain. Previous experiments also indicate that most of the transcriptional activity of CLK–CYC derives from the glutamine-rich region of its partner CLK. To address the role of transcription in core circadian timekeeping, we have analyzed the effects of a CYC–viral protein 16 (VP16) fusion protein in the Drosophila system. The addition of this potent and well-studied viral transcriptional activator (VP16) to CYC imparts to the CLK–CYC-VP16 complex strongly enhanced transcriptional activity relative to that of CLK–CYC. This increase is manifested in flies expressing CYC-VP16 as well as in S2 cells. These flies also have increased levels of CLK–CYC direct target gene mRNAs as well as a short period, implicating circadian transcription in period determination. A more detailed examination of reporter gene expression in CYC-VP16–expressing flies suggests that the short period is due at least in part to a more rapid transcriptional phase. Importantly, the behavioral effects require a period (per) promoter and are therefore unlikely to be merely a consequence of generally higher PER levels. This indicates that the CLK–CYC-VP16 behavioral effects are a consequence of increased per transcription. All of this also suggests that the timing of transcriptional activation and not the activation itself is the key event responsible for the behavioral effects observed in CYC-VP16-expressing flies. The results taken together indicate that circadian transcription contributes to core circadian function in Drosophila