Methyl-5-uridine (m5U) is one the most abundant non-canonical bases present in cellular RNA, and in yeast is found at position U54 of tRNAs where modification is catalysed by the methyltransferase Trm2. Although the mammalian enzymes that catalyse m5U formation are yet to be identified via experimental evidence, based on sequence homology to Trm2, two candidates currently exist, TRMT2A and TRMT2B. Here we developed a genome-wide single-nucleotide resolution mapping method, Fluorouracil-Induced-Catalytic-Crosslinking-Sequencing (FICC-Seq), in order to identify the relevant enzymatic targets. We demonstrate that TRMT2A is responsible for the majority of m5U present in human RNA, and that it commonly targets U54 of cytosolic tRNAs. By comparison to current methods, we show that FICC-Seq is a particularly robust method for accurate and reliable detection of relevant enzymatic target sites. Our associated finding of extensive irreversible TRMT2A-tRNA crosslinking in vivo following 5-Fluorouracil exposure is also intriguing, as it suggests a tangible mechanism for a previously suspected RNA-dependent route of Fluorouracil-mediated cytotoxicity.

Recent studies have revealed multiple mechanisms that can lead to heterogeneity in ribosomal composition. This heterogeneity can lead to preferential translation of specific panels of mRNAs, and is defined in large part by the ribosomal protein (RP) content, amongst other things. However, it is currently unknown to what extent ribosomal composition is heterogeneous across tissues, which is compounded by a lack of tools available to study it. Here we present dripARF, a method for detecting differential RP incorporation into the ribosome using Ribosome Profiling (Ribo-seq) data. We combine the 'waste' rRNA fragment data generated in Ribo-seq with the known 3D structure of the human ribosome to predict differences in the composition of ribosomes in the material being studied. We have validated this approach using publicly available data, and have revealed a potential role for eS25/RPS25 in development. Our results indicate that ribosome heterogeneity can be detected in Ribo-seq data, providing a new method to study this phenomenon. Furthermore, with dripARF, previously published Ribo-seq data provides a wealth of new information, allowing the identification of RPs of interest in many disease and normal contexts. dripARF is available as part of the ARF R package and can be accessed through https://github.com/fallerlab/ARF.

The structure of mRNA molecules plays an important role in its interactions with trans-acting factors, notably RNA binding proteins (RBPs), thus contributing to the functional consequences of this interplay. However, current transcriptome-wide experimental methods to chart these interactions are limited by their poor sensitivity. Here we extend the hiCLIP atlas of duplexes bound by Staufen1 (STAU1) ∼10-fold, through careful consideration of experimental assumptions, and the development of bespoke computational methods which we apply to existing data. We present Tosca, a Nextflow computational pipeline for the processing, analysis and visualisation of proximity ligation sequencing data generally. We use our extended duplex atlas to discover insights into the RNA selectivity of STAU1, revealing the importance of structural symmetry and duplex-span-dependent nucleotide composition. Furthermore, we identify heterogeneity in the relationship between transcripts with STAU1-bound 3' UTR duplexes and metabolism of the associated RNAs that we relate to RNA structure: transcripts with short-range proximal 3' UTR duplexes have high degradation rates, but those with long-range duplexes have low rates. Overall, our work enables the integrative analysis of proximity ligation data delivering insights into specific features and effects of RBP-RNA structure interactions.

Schlafen-5 (SLFN5) is an interferon-induced protein of the Schlafen family, which are involved in immune responses and oncogenesis. To date, little is known regarding its anti-HIV-1 function. Here, the authors report that overexpression of SLFN5 inhibits HIV-1 replication and reduces viral mRNA levels, whereas depletion of endogenous SLFN5 promotes HIV-1 replication. Moreover, they show that SLFN5 markedly decreases the transcriptional activity of HIV-1 long terminal repeat (LTR) via binding to two sequences in the U5-R region, which consequently represses the recruitment of RNA polymerase II to the transcription initiation site. Mutagenesis studies show the importance of nuclear localization and the N-terminal 1-570 amino acids fragment in the inhibition of HIV-1. Further mechanistic studies demonstrate that SLFN5 interacts with components of the PRC2 complex, G9a and Histone H3, thereby promoting H3K27me2 and H3K27me3 modification leading to silencing HIV-1 transcription. In concert with this, they find that SLFN5 blocks the activation of latent HIV-1. Altogether, their findings demonstrate that SLFN5 is a transcriptional repressor of HIV-1 through epigenetic modulation and a potential determinant of HIV-1 latency.

The 3' untranslated regions (3' UTRs) of transcripts serve as important hubs for posttranscriptional gene expression regulation. Here, we find that the exonisation of intergenic Alu elements introduced new terminal exons and polyadenylation sites during human genome evolution. While Alu exonisation from introns has been described previously, we shed light on a novel mechanism to create alternative 3' UTRs, thereby opening opportunities for differential posttranscriptional regulation. On the mechanistic level, we show that intergenic Alu exonisation can compete both with alternative splicing and polyadenylation in the upstream gene. Notably, the Alu-derived isoforms are often expressed in a tissue-specific manner, and the Alu-derived 3' UTRs can alter mRNA stability. In summary, we demonstrate that intergenic elements can affect processing of preceding genes, and elucidate how intergenic Alu exonisation can contribute to tissue-specific posttranscriptional regulation by expanding the repertoire of 3' UTRs.

Enhancers, as specialized genomic cis-regulatory elements, activate transcription of their target genes and play an important role in pathogenesis of many human complex diseases. Despite recent systematic identification of them in the human genome, currently there is an urgent need for comprehensive annotation databases of human enhancers with a focus on their disease connections. In response, we built the Human Enhancer Disease Database (HEDD) to facilitate studies of enhancers and their potential roles in human complex diseases. HEDD currently provides comprehensive genomic information for ∼2.8 million human enhancers identified by ENCODE, FANTOM5 and RoadMap with disease association scores based on enhancer-gene and gene-disease connections. It also provides Web-based analytical tools to visualize enhancer networks and score enhancers given a set of selected genes in a specific gene network. HEDD is freely accessible at http://zdzlab.einstein.yu.edu/1/hedd.php.

Alternative splicing (AS) is a key component of gene expression programs that drive cellular differentiation. Smooth muscle cells (SMCs) are important in the function of a number of physiological systems; however, investigation of SMC AS has been restricted to a handful of events. We profiled transcriptome changes in mouse de-differentiating SMCs and observed changes in hundreds of AS events. Exons included in differentiated cells were characterized by particularly weak splice sites and by upstream binding sites for Polypyrimidine Tract Binding protein (PTBP1). Consistent with this, knockdown experiments showed that that PTBP1 represses many smooth muscle specific exons. We also observed coordinated splicing changes predicted to downregulate the expression of core components of U1 and U2 snRNPs, splicing regulators and other post-transcriptional factors in differentiated cells. The levels of cognate proteins were lower or similar in differentiated compared to undifferentiated cells. However, levels of snRNAs did not follow the expression of splicing proteins, and in the case of U1 snRNP we saw reciprocal changes in the levels of U1 snRNA and U1 snRNP proteins. Our results suggest that the AS program in differentiated SMCs is orchestrated by the combined influence of auxiliary RNA binding proteins, such as PTBP1, along with altered activity and stoichiometry of the core splicing machinery.

RNA-binding proteins (RBPs) play diverse roles in regulating co-transcriptional RNA-processing and chromatin functions, but our knowledge of the repertoire of chromatin-associated RBPs (caRBPs) and their interactions with chromatin remains limited. Here, we developed SPACE (Silica Particle Assisted Chromatin Enrichment) to isolate global and regional chromatin components with high specificity and sensitivity, and SPACEmap to identify the chromatin-contact regions in proteins. Applied to mouse embryonic stem cells, SPACE identified 1459 chromatin-associated proteins, ∼48% of which are annotated as RBPs, indicating their dual roles in chromatin and RNA-binding. Additionally, SPACEmap stringently verified chromatin-binding of 403 RBPs and identified their chromatin-contact regions. Notably, SPACEmap showed that about 40% of the caRBPs bind chromatin by intrinsically disordered regions (IDRs). Studying SPACE and total proteome dynamics from mES cells grown in 2iL and serum medium indicates significant correlation (R = 0.62). One of the most dynamic caRBPs is Dazl, which we find co-localized with PRC2 at transcription start sites of genes that are distinct from Dazl mRNA binding. Dazl and other PRC2-colocalised caRBPs are rich in intrinsically disordered regions (IDRs), which could contribute to the formation and regulation of phase-separated PRC condensates. Together, our approach provides an unprecedented insight into IDR-mediated interactions and caRBPs with moonlighting functions in native chromatin.

Extracellular vesicles (EVs) are small membranous vesicles that contain an abundant cargo of different RNA species with specialized functions and clinical implications. Here, we introduce an updated online database (http://www.exoRBase.org), exoRBase 2.0, which is a repository of EV long RNAs (termed exLRs) derived from RNA-seq data analyses of diverse human body fluids. In exoRBase 2.0, the number of exLRs has increased to 19 643 messenger RNAs (mRNAs), 15 645 long non-coding RNAs (lncRNAs) and 79 084 circular RNAs (circRNAs) obtained from ∼1000 human blood, urine, cerebrospinal fluid (CSF) and bile samples. Importantly, exoRBase 2.0 not only integrates and compares exLR expression profiles but also visualizes the pathway-level functional changes and the heterogeneity of origins of circulating EVs in the context of different physiological and pathological conditions. Our database provides an attractive platform for the identification of novel exLR signatures from human biofluids that will aid in the discovery of new circulating biomarkers to improve disease diagnosis and therapy.

The non-canonical initiation factor DENR promotes translation reinitiation on mRNAs harbouring upstream open reading frames (uORFs). Moreover, DENR depletion shortens circadian period in mouse fibroblasts, suggesting involvement of uORF usage and reinitiation in clock regulation. To identify DENR-regulated translation events transcriptome-wide and, in particular, specific core clock transcripts affected by this mechanism, we have used ribosome profiling in DENR-deficient NIH3T3 cells. We uncovered 240 transcripts with altered translation rate, and used linear regression analysis to extract 5' UTR features predictive of DENR dependence. Among core clock genes, we identified Clock as a DENR target. Using Clock 5' UTR mutants, we mapped the specific uORF through which DENR acts to regulate CLOCK protein biosynthesis. Notably, these experiments revealed an alternative downstream start codon, likely representing the bona fide CLOCK N-terminus. Our findings provide insights into uORF-mediated translational regulation that can regulate the mammalian circadian clock and gene expression at large.