The RNA exosome shapes the expression of key protein-coding genes

True length of diverse capped RNA sequencing (TLDR-seq): 5'-3'-end sequencing of capped RNAs regardless of 3'-end status

Analysis of transcript function is greatly aided by knowledge of the full-length RNA sequence. New long-read sequencing enabled by Oxford Nanopore and PacBio devices have the potential to provide full-length transcript information; however, standard methods still lack the ability to capture true RNA 5' ends and select for polyadenylated (pA+) transcripts only. Here, we present a method that, by utilizing cap trapping and 3'-end adapter ligation, sequences transcripts between their exact 5' and 3' ends regardless of polyadenylation status and without the need for ribosomal RNA depletion, with the ability to characterize polyadenylation length of RNAs, if any. The method shows high reproducibility, can faithfully detect 5' ends, 3' ends and splice junctions, and produces gene-expression estimates that are highly correlated to those of short-read sequencing techniques. We also demonstrate that the method can detect and sequence full-length nonadenylated (pA-) RNAs, including long noncoding RNAs, promoter upstream transcripts, and enhancer RNAs, and present cases where pA+ and pA- RNAs show preferences for different but closely located transcription start sites. Our method is therefore useful for the characterization of diverse capped RNA species and analysis of relationships between transcription initiation, termination, and RNA processing.

Nuclear sorting of short RNA polymerase II transcripts

Mammalian genomes produce an abundance of short RNA. This is, to a large extent, due to the genome-wide and spurious activity of RNA polymerase II (RNAPII). However, it is also because the vast majority of initiating RNAPII, regardless of the transcribed DNA unit, terminates within a ∼3-kb early "pausing zone." Given that the resultant RNAs constitute both functional and non-functional species, their proper sorting is critical. One way to think about such quality control (QC) is that transcripts, from their first emergence, are relentlessly targeted by decay factors, which may only be avoided by engaging protective processing pathways. In a molecular materialization of this concept, recent progress has found that both "destructive" and "productive" RNA effectors assemble at the 5' end of capped RNA, orchestrated by the essential arsenite resistance protein 2 (ARS2) protein. Based on this principle, we here discuss early QC mechanisms and how these might sort short RNAs to their final fates.

Nuclear mRNA decay: regulatory networks that control gene expression

Proper regulation of mRNA production in the nucleus is critical for the maintenance of cellular homoeostasis during adaptation to internal and environmental cues. Over the past 25 years, it has become clear that the nuclear machineries governing gene transcription, pre-mRNA processing, pre-mRNA and mRNA decay, and mRNA export to the cytoplasm are inextricably linked to control the quality and quantity of mRNAs available for translation. More recently, an ever-expanding diversity of new mechanisms by which nuclear RNA decay factors finely tune the expression of protein-encoding genes have been uncovered. Here, we review the current understanding of how mammalian cells shape their protein-encoding potential by regulating the decay of pre-mRNAs and mRNAs in the nucleus.

Autism candidate gene rbm-26 (RBM26/27) regulates MALS-1 to protect against mitochondrial dysfunction and axon degeneration during neurodevelopment

Mitochondrial dysfunction is thought to be a key component of neurodevelopmental disorders such as autism, intellectual disability, and ADHD. However, little is known about the molecular mechanisms that protect against mitochondrial dysfunction during neurodevelopment. Here, we address this question through the investigation of rbm-26, the C. elegans ortholog of the RBM27 autism candidate gene, which encodes an RNA-binding protein whose role in neurons is unknown. We report that RBM-26 (RBM26/27) protects against axonal defects by negatively regulating expression of the MALS-1 (MALSU1) mitoribosomal assembly factor. Autism-associated missense variants in RBM-26 cause a sharp decrease in RBM-26 protein expression along with defects in in axon overlap and axon degeneration that occurs during larval development. Using a biochemical screen, we identified the mRNA for the MALS-1 mitoribosomal assembly factor as a binding partner for RBM-26. Loss of RBM-26 function causes a dramatic overexpression of mals-1 mRNA and MALS-1 protein. Moreover, genetic analysis indicates that this overexpression of MALS-1 is responsible for the mitochondrial and axon degeneration defects in rbm-26 mutants. These observations reveal a mechanism that regulates expression of a mitoribosomal assembly factor to protect against axon degeneration during neurodevelopment.

Genome-wide quantification of RNA flow across subcellular compartments reveals determinants of the mammalian transcript life cycle

Dissecting the regulatory mechanisms controlling mammalian transcripts from production to degradation requires quantitative measurements of mRNA flow across the cell. We developed subcellular TimeLapse-seq to measure the rates at which RNAs are released from chromatin, exported from the nucleus, loaded onto polysomes, and degraded within the nucleus and cytoplasm in human and mouse cells. These rates varied substantially, yet transcripts from genes with related functions or targeted by the same transcription factors and RNA-binding proteins flowed across subcellular compartments with similar kinetics. Verifying these associations uncovered a link between DDX3X and nuclear export. For hundreds of RNA metabolism genes, most transcripts with retained introns were degraded by the nuclear exosome, while the remaining molecules were exported with stable cytoplasmic lifespans. Transcripts residing on chromatin for longer had extended poly(A) tails, whereas the reverse was observed for cytoplasmic mRNAs. Finally, machine learning identified molecular features that predicted the diverse life cycles of mRNAs.

A data-driven genome annotation approach for cassava

Genome annotation files play a critical role in dictating the quality of downstream analyses by providing essential predictions for gene positions and structures. These files are pivotal in decoding the complex information encoded within DNA sequences. Here, we generated experimental data resolving RNA 5'- and 3'-ends as well as full-length RNAs for cassava TME12 sticklings in ambient temperature and cold. We used these data to generate genome annotation files using the TranscriptomeReconstructoR (TR) tool. A careful comparison to high-quality genome annotations suggests that our new TR genome annotations identified additional genes, resolved the transcript boundaries more accurately and identified additional RNA isoforms. We enhanced existing cassava genome annotation files with the information from TR that maintained the different transcript models as RNA isoforms. The resultant merged annotation was subsequently utilized for comprehensive analysis. To examine the effects of genome annotation files on gene expression studies, we compared the detection of differentially expressed genes during cold using the same RNA-seq data but alternative genome annotation files. We found that our merged genome annotation that included cold-specific TR gene models identified about twice as many cold-induced genes. These data indicate that environmentally induced genes may be missing in off-the-shelf genome annotation files. In conclusion, TR offers the opportunity to enhance crop genome annotations with implications for the discovery of differentially expressed candidate genes during plant-environment interactions.

Dysregulation of RNA-Exosome machinery is directly linked to major cancer hallmarks in prostate cancer: Oncogenic role of PABPN1

Novel biomarkers and therapeutic strategies for prostate-cancer (PCa) are required to overcome its lethal progression. The dysregulation/implication of the RNA-Exosome-complex (REC; cellular machinery controlling the 3'-5'processing/degradation of most RNAs) in different cancer-types, including PCa, is poorly known. Herein, different cellular/molecular/preclinical approaches with human PCa-samples (tissues and/or plasma of 7 independent cohorts), and in-vitro/in-vivo PCa-models were used to comprehensively characterize the REC-profile and explore its role in PCa. Moreover, isoginkgetin (REC-inhibitor) effects were evaluated on PCa-cells. We demonstrated a specific dysregulation of the REC-components in PCa-tissues, identifying the Poly(A)-Binding-Protein-Nuclear 1 (PABPN1) factor as a critical regulator of major cancer hallmarks. PABPN1 is consistently overexpressed in different human PCa-cohorts and associated with poor-progression, invasion and metastasis. PABPN1 silencing decreased relevant cancer hallmarks in multiple PCa-models (proliferation/migration/tumourspheres/colonies, etc.) through the modulation of key cancer-related lncRNAs (PCA3/FALEC/DLEU2) and mRNAs (CDK2/CDK6/CDKN1A). Plasma PABPN1 levels were altered in patients with metastatic and tumour-relapse. Finally, pharmacological inhibition of REC-activity drastically inhibited PCa-cell aggressiveness. Altogether, the REC is drastically dysregulated in PCa, wherein this novel molecular event/mechanism, especially PABPN1 alteration, may be potentially exploited as a novel prognostic and therapeutic tool for PCa.

Temporal-iCLIP captures co-transcriptional RNA-protein interactions

Dynamic RNA-protein interactions govern the co-transcriptional packaging of RNA polymerase II (RNAPII)-derived transcripts. Yet, our current understanding of this process in vivo primarily stems from steady state analysis. To remedy this, we here conduct temporal-iCLIP (tiCLIP), combining RNAPII transcriptional synchronisation with UV cross-linking of RNA-protein complexes at serial timepoints. We apply tiCLIP to the RNA export adaptor, ALYREF; a component of the Nuclear Exosome Targeting (NEXT) complex, RBM7; and the nuclear cap binding complex (CBC). Regardless of function, all tested factors interact with nascent RNA as it exits RNAPII. Moreover, we demonstrate that the two transesterification steps of pre-mRNA splicing temporally separate ALYREF and RBM7 binding to splicing intermediates, and that exon-exon junction density drives RNA 5'end binding of ALYREF. Finally, we identify underappreciated steps in snoRNA 3'end processing performed by RBM7. Altogether, our data provide a temporal view of RNA-protein interactions during the early phases of transcription.

RNA-regulatory exosome complex suppresses an apoptotic program to confer erythroid progenitor cell survival in vivo

The RNA-regulatory exosome complex (EC) posttranscriptionally and cotranscriptionally processes and degrades RNAs in a context-dependent manner. Although the EC functions in diverse cell types, its contributions to stem and progenitor cell development are not well understood. Previously, we demonstrated that the transcriptional regulator of erythrocyte development, GATA1, represses EC subunit genes, and the EC maintains erythroid progenitors in vitro. To determine if this mechanism operates in vivo, we used the hematopoietic-specific Vav1-Cre and "conditional by inversion" mouse system to ablate Exosc3, encoding an EC structural subunit. Although Exosc3C/C Cre+ embryos developed normally until embryonic day 14.5, Exosc3 ablation was embryonic lethal and severely reduced erythromyeloid progenitor activity. RNA sequencing analysis of Exosc3-ablated burst-forming unit-erythroid revealed elevated transcripts encoding multiple proapoptotic factors, and the mutant erythroid progenitors exhibited increased apoptosis. We propose that the EC controls an ensemble of apoptosis-regulatory RNAs, thereby promoting erythroid progenitor survival and developmental erythropoiesis in vivo.

Computational identification of signals predictive for nuclear RNA exosome degradation pathway targeting

The RNA exosome degrades transcripts in the nucleoplasm of mammalian cells. Its substrate specificity is mediated by two adaptors: the ‘nuclear exosome targeting (NEXT)’ complex and the ‘poly(A) exosome targeting (PAXT)’ connection. Previous studies have revealed some DNA/RNA elements that differ between the two pathways, but how informative these features are for distinguishing pathway targeting, or whether additional genomic features that are informative for such classifications exist, is unknown. Here, we leverage the wealth of available genomic data and develop machine learning models that predict exosome targets and subsequently rank the features the models use by their predictive power. As expected, features around transcript end sites were most predictive; specifically, the lack of canonical 3′ end processing was highly predictive of NEXT targets. Other associated features, such as promoter-proximal G/C content and 5′ splice sites, were informative, but only for distinguishing NEXT and not PAXT targets. Finally, we discovered predictive features not previously associated with exosome targeting, in particular RNA helicase DDX3X binding sites. Overall, our results demonstrate that nucleoplasmic exosome targeting is to a large degree predictable, and our approach can assess the predictive power of previously known and new features in an unbiased way.

Chromatin modifier HUSH co-operates with RNA decay factor NEXT to restrict transposable element expression

Transposable elements (TEs) are widespread genetic parasites known to be kept under tight transcriptional control. Here, we describe a functional connection between the mouse-orthologous “nuclear exosome targeting” (NEXT) and “human silencing hub” (HUSH) complexes, involved in nuclear RNA decay and the epigenetic silencing of TEs, respectively. Knocking out the NEXT component ZCCHC8 in embryonic stem cells results in elevated TE RNA levels. We identify a physical interaction between ZCCHC8 and the MPP8 protein of HUSH and establish that HUSH recruits NEXT to chromatin at MPP8-bound TE loci. However, while NEXT and HUSH both dampen TE RNA expression, their activities predominantly affect shorter non-polyadenylated and full-length polyadenylated transcripts, respectively. Indeed, our data suggest that the repressive action of HUSH promotes a condition favoring NEXT RNA decay activity. In this way, transcriptional and post-transcriptional machineries synergize to suppress the genotoxic potential of TE RNAs.

Garland et al. report a physical and functional connection between the NEXT complex, involved in RNA decay, and the HUSH complex, involved in chromatin regulation. Together, NEXT and HUSH cooperate to control transposable element (TE) RNA expression in embryonic stem cells, suppressing pA⁻ and pA⁺ transcripts, respectively.

Rapid factor depletion highlights intricacies of nucleoplasmic RNA degradation

Turnover of nucleoplasmic transcripts by the mammalian multi-subunit RNA exosome is mediated by two adaptors: the Nuclear EXosome Targeting (NEXT) complex and the Poly(A) tail eXosome Targeting (PAXT) connection. Functional analyses of NEXT and PAXT have largely utilized long-term factor depletion strategies, facilitating the appearance of indirect phenotypes. Here, we rapidly deplete NEXT, PAXT and core exosome components, uncovering the direct consequences of their acute losses. Generally, proteome changes are sparse and largely dominated by co-depletion of other exosome and adaptor subunits, reflecting possible subcomplex compositions. While parallel high-resolution 3′ end sequencing of newly synthesized RNA confirms previously established factor specificities, it concomitantly demonstrates an inflation of long-term depletion datasets by secondary effects. Most strikingly, a general intron degradation phenotype, observed in long-term NEXT depletion samples, is undetectable upon short-term depletion, which instead emphasizes NEXT targeting of snoRNA-hosting introns. Further analysis of these introns uncovers an unusual mode of core exosome-independent RNA decay. Our study highlights the accumulation of RNAs as an indirect result of long-term decay factor depletion, which we speculate is, at least partly, due to the exhaustion of alternative RNA decay pathways.

Using TIF-Seq2 to investigate association between 5´ and 3´mRNA ends

The development of high-throughput technologies has revealed pervasive transcription in all genomes that have been investigated so far. This has uncovered a highly interleaved transcriptome organization involving thousands of overlapping coding and non-coding RNA isoforms that challenge our traditional definitions of genes and functional regions of the genome. In this chapter, we discuss the application of an improved Transcript Isoform Sequencing approach (TIF-Seq2) able to concurrently determine the start and end sites of individual RNA molecules. We exemplify its use for the investigation of the human transcriptome and show how it is especially well suited to discriminate between overlapping molecules and accurately define their boundaries.