ELAV links paused Pol II to alternative polyadenylation in the Drosophila nervous system

Dynamic regulation of alternative polyadenylation by PQBP1 during neurogenesis

Alternative polyadenylation (APA) is a critical post-transcriptional process that generates mRNA isoforms with distinct 3' untranslated regions (3' UTRs), thereby regulating mRNA localization, stability, and translational efficiency. Cell-type-specific APA extensively shapes the diversity of the cellular transcriptome, particularly during cell fate transition. Despite its recognized significance, the precise regulatory mechanisms governing cell-type-specific APA remain unclear. In this study, we uncover PQBP1 as an emerging APA regulator that actively maintains cell-specific APA profiles in neural progenitor cells (NPCs) and delicately manages the equilibrium between NPC proliferation and differentiation. Multi-omics analysis shows that PQBP1 directly interacts with the upstream UGUA elements, impeding the recruitment of the CFIm complex and influencing polyadenylation site selection within genes associated with the cell cycle. Our findings elucidate the molecular mechanism by which PQBP1 orchestrates dynamic APA changes during neurogenesis, providing valuable insights into the precise regulation of cell-type-specific APA and the underlying pathogenic mechanisms in neurodevelopmental disorders.

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

(Alternative) transcription start sites as regulators of RNA processing

Alternative transcription start site usage (ATSS) is a widespread regulatory strategy that enables genes to choose between multiple genomic loci for initiating transcription. This mechanism is tightly controlled during development and is often altered in disease states. In this review, we examine the growing evidence highlighting a role for transcription start sites (TSSs) in the regulation of mRNA isoform selection during and after transcription. We discuss how the choice of transcription initiation sites influences RNA processing and the importance of this crosstalk for cell identity and organism function. We also speculate on possible mechanisms underlying the integration of transcriptional and post-transcriptional processes.

Poly(A) tale: From A to A; RNA polyadenylation in prokaryotes and eukaryotes

Most eukaryotic mRNAs and different non-coding RNAs undergo a form of 3' end processing known as polyadenylation. Polyadenylation machinery is present in almost all organisms except few species. In bacteria, the machinery has evolved from PNPase, which adds heteropolymeric tails, to a poly(A)-specific polymerase. Differently, a complex machinery for accurate polyadenylation and several non-canonical poly(A) polymerases are developed in eukaryotes. The role of poly(A) tail has also evolved from serving as a degradative signal to a stabilizing modification that also regulates translation. In this review, we discuss poly(A) tail emergence in prokaryotes and its development into a stable, yet dynamic feature at the 3' end of mRNAs in eukaryotes. We also describe how appearance of novel poly(A) polymerases gives cells flexibility to shape poly(A) tail. We explain how poly(A) tail dynamics help regulate cognate RNA metabolism in a context-dependent manner, such as during oocyte maturation. Finally, we describe specific mRNAs in metazoans that bear stem-loops instead of poly(A) tails. We conclude with how recent discoveries about poly(A) tail can be applied to mRNA technology. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Processing > 3' End Processing RNA Turnover and Surveillance > Regulation of RNA Stability.

Alternative cleavage and polyadenylation of the Ccnb1 mRNA defines accumulation of cyclin protein during the meiotic cell cycle

Progression through the mitotic and meiotic cell cycle is driven by fluctuations in the levels of cyclins, the regulatory subunits controlling the localization and activity of CDK1 kinases. Cyclin levels are regulated through a precise balance of synthesis and degradation. Here we demonstrate that the synthesis of Cyclin B1 during the oocyte meiotic cell cycle is defined by the selective translation of mRNA variants generated through alternative cleavage and polyadenylation (APA). Using gene editing in mice, we introduced mutations into the proximal and distal polyadenylation elements of the 3' untranslated region (UTR) of the Ccnb1 mRNA. Through in vivo loss-of-function experiments, we demonstrate that the translation of mRNA with a short 3' UTR specifies Cyclin B1 protein levels that set the timing of meiotic re-entry. In contrast, translation directed by a long 3' UTR is necessary to direct Cyclin B1 protein accumulation during the MI/MII transition. These findings establish that the progression through the cell cycle is dependent on the selective translation of multiple mRNA variants generated by APA.

YTHDF1 facilitates PRC1-mediated H2AK119ub in human ES cells

Polycomb repressive complexes (PRCs) play critical roles in cell fate decisions during normal development as well as disease progression through mediating histone modifications such as H3K27me3 and H2AK119ub. How exactly PRCs recruited to chromatin remains to be fully illuminated. Here, we report that YTHDF1, the N6-methyladenine (m6 A) RNA reader that was previously known to be mainly cytoplasmic, associates with RNF2, a PRC1 protein that mediates H2AK119ub in human embryonic stem cells (hESCs). A portion of YTHDF1 localizes in the nuclei and associates with RNF2/H2AK119ub on a subset of gene loci related to neural development functions. Knock-down YTHDF1 attenuates H2AK119ub modification on these genes and promotes neural differentiation in hESCs. Our findings provide a noncanonical mechanism that YTHDF1 participates in PRC1 functions in hESCs.

Regulation of alternative splicing and polyadenylation in neurons

Cell-type-specific gene expression is a fundamental feature of multicellular organisms and is achieved by combinations of regulatory strategies. Although cell-restricted transcription is perhaps the most widely studied mechanism, co-transcriptional and post-transcriptional processes are also central to the spatiotemporal control of gene functions. One general category of expression control involves the generation of multiple transcript isoforms from an individual gene, whose balance and cell specificity are frequently tightly regulated via diverse strategies. The nervous system makes particularly extensive use of cell-specific isoforms, specializing the neural function of genes that are expressed more broadly. Here, we review regulatory strategies and RNA-binding proteins that direct neural-specific isoform processing. These include various classes of alternative splicing and alternative polyadenylation events, both of which broadly diversify the neural transcriptome. Importantly, global alterations of splicing and alternative polyadenylation are characteristic of many neural pathologies, and recent genetic studies demonstrate how misregulation of individual neural isoforms can directly cause mutant phenotypes.

stAPAminer: Mining Spatial Patterns of Alternative Polyadenylation for Spatially Resolved Transcriptomic Studies

Alternative polyadenylation (APA) contributes to transcriptome complexity and gene expression regulation and has been implicated in various cellular processes and diseases. Single-cell RNA sequencing (scRNA-seq) has enabled the profiling of APA at the single-cell level; however, the spatial information of cells is not preserved in scRNA-seq. Alternatively, spatial transcriptomics (ST) technologies provide opportunities to decipher the spatial context of the transcriptomic landscape. Pioneering studies have revealed potential spatially variable genes and/or splice isoforms; however, the pattern of APA usage in spatial contexts remains unappreciated. In this study, we developed a toolkit called stAPAminer for mining spatial patterns of APA from spatially barcoded ST data. APA sites were identified and quantified from the ST data. In particular, an imputation model based on the k-nearest neighbors algorithm was designed to recover APA signals, and then APA genes with spatial patterns of APA usage variation were identified. By analyzing well-established ST data of the mouse olfactory bulb (MOB), we presented a detailed view of spatial APA usage across morphological layers of the MOB. We compiled a comprehensive list of genes with spatial APA dynamics and obtained several major spatial expression patterns that represent spatial APA dynamics in different morphological layers. By extending this analysis to two additional replicates of the MOB ST data, we observed that the spatial APA patterns of several genes were reproducible among replicates. stAPAminer employs the power of ST to explore the transcriptional atlas of spatial APA patterns with spatial resolution. This toolkit is available at https://github.com/BMILAB/stAPAminer and https://ngdc.cncb.ac.cn/biocode/tools/BT007320.

Sites of transcription initiation drive mRNA isoform selection

The generation of distinct messenger RNA isoforms through alternative RNA processing modulates the expression and function of genes, often in a cell-type-specific manner. Here, we assess the regulatory relationships between transcription initiation, alternative splicing, and 3' end site selection. Applying long-read sequencing to accurately represent even the longest transcripts from end to end, we quantify mRNA isoforms in Drosophila tissues, including the transcriptionally complex nervous system. We find that in Drosophila heads, as well as in human cerebral organoids, 3' end site choice is globally influenced by the site of transcription initiation (TSS). "Dominant promoters," characterized by specific epigenetic signatures including p300/CBP binding, impose a transcriptional constraint to define splice and polyadenylation variants. In vivo deletion or overexpression of dominant promoters as well as p300/CBP loss disrupted the 3' end expression landscape. Our study demonstrates the crucial impact of TSS choice on the regulation of transcript diversity and tissue identity.

Regulation of neuronal RNA signatures by ELAV/Hu proteins

The RNA-binding proteins encoded by the highly conserved elav/Hu gene family, found in all metazoans, regulate the expression of a wide range of genes, at both the co-transcriptional and posttranscriptional level. Nervous-system-specific ELAV/Hu proteins are prominent for their essential role in neuron differentiation, and mutations have been associated with human neurodevelopmental and neurodegenerative diseases. Drosophila ELAV, the founding member of the protein family, mediates the synthesis of neuronal RNA signatures by promoting alternative splicing and alternative polyadenylation of hundreds of genes. The recent identification of ELAV's direct RNA targets revealed the protein's central role in shaping the neuronal transcriptome, and highlighted the importance of neuronal transcript signatures for neuron maintenance and organism survival. Animals have evolved multiple cellular mechanisms to ensure robustness of ELAV/Hu function. In Drosophila, elav autoregulates in a 3'UTR-dependent manner to maintain optimal protein levels. A complete absence of ELAV causes the activation and nuclear localization of the normally cytoplasmic paralogue FNE, in a process termed EXon-Activated functional Rescue (EXAR). Other species, including mammals, seem to utilize different strategies, such as protein redundancy, to maintain ELAV protein function and effectively safeguard the identity of the neuronal transcriptome. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Development RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

Regulation and function of alternative polyadenylation in development and differentiation

Alternative processing of nascent mRNAs is widespread in eukaryotic organisms and greatly impacts the output of gene expression. Specifically, alternative cleavage and polyadenylation (APA) is a co-transcriptional molecular process that switches the polyadenylation site (PAS) at which a nascent mRNA is cleaved, resulting in mRNA isoforms with different 3'UTR length and content. APA can potentially affect mRNA translation efficiency, localization, stability, and mRNA seeded protein-protein interactions. APA naturally occurs during development and cellular differentiation, with around 70% of human genes displaying APA in particular tissues and cell types. For example, neurons tend to express mRNAs with long 3'UTRs due to preferential processing at PASs more distal than other PASs used in other cell types. In addition, changes in APA mark a variety of pathological states, including many types of cancer, in which mRNAs are preferentially cleaved at more proximal PASs, causing expression of mRNA isoforms with short 3'UTRs. Although APA has been widely reported, both the function of APA in development and the mechanisms that regulate the choice of 3'end cut sites in normal and pathogenic conditions are still poorly understood. In this review, we summarize current understanding of how APA is regulated during development and cellular differentiation and how the resulting change in 3'UTR content affects multiple aspects of gene expression. With APA being a widespread phenomenon, the advent of cutting-edge scientific techniques and the pressing need for in-vivo studies, there has never been a better time to delve into the intricate mechanisms of alternative cleavage and polyadenylation.

Context-specific regulation and function of mRNA alternative polyadenylation

Alternative cleavage and polyadenylation (APA) is a widespread mechanism to generate mRNA isoforms with alternative 3' untranslated regions (UTRs). The expression of alternative 3' UTR isoforms is highly cell type specific and is further controlled in a gene-specific manner by environmental cues. In this Review, we discuss how the dynamic, fine-grained regulation of APA is accomplished by several mechanisms, including cis-regulatory elements in RNA and DNA and factors that control transcription, pre-mRNA cleavage and post-transcriptional processes. Furthermore, signalling pathways modulate the activity of these factors and integrate APA into gene regulatory programmes. Dysregulation of APA can reprogramme the outcome of signalling pathways and thus can control cellular responses to environmental changes. In addition to the regulation of protein abundance, APA has emerged as a major regulator of mRNA localization and the spatial organization of protein synthesis. This role enables the regulation of protein function through the addition of post-translational modifications or the formation of protein-protein interactions. We further discuss recent transformative advances in single-cell RNA sequencing and CRISPR-Cas technologies, which enable the mapping and functional characterization of alternative 3' UTRs in any biological context. Finally, we discuss new APA-based RNA therapeutics, including compounds that target APA in cancer and therapeutic genome editing of degenerative diseases.

RNA-Binding Proteins as Epigenetic Regulators of Brain Functions and Their Involvement in Neurodegeneration

A central aspect of nervous system development and function is the post-transcriptional regulation of mRNA fate, which implies time- and site-dependent translation, in response to cues originating from cell-to-cell crosstalk. Such events are fundamental for the establishment of brain cell asymmetry, as well as of long-lasting modifications of synapses (long-term potentiation: LTP), responsible for learning, memory, and higher cognitive functions. Post-transcriptional regulation is in turn dependent on RNA-binding proteins that, by recognizing and binding brief RNA sequences, base modifications, or secondary/tertiary structures, are able to control maturation, localization, stability, and translation of the transcripts. Notably, most RBPs contain intrinsically disordered regions (IDRs) that are thought to be involved in the formation of membrane-less structures, probably due to liquid-liquid phase separation (LLPS). Such structures are evidenced as a variety of granules that contain proteins and different classes of RNAs. The other side of the peculiar properties of IDRs is, however, that, under altered cellular conditions, they are also prone to form aggregates, as observed in neurodegeneration. Interestingly, RBPs, as part of both normal and aggregated complexes, are also able to enter extracellular vesicles (EVs), and in doing so, they can also reach cells other than those that produced them.

A Survey on Methods for Predicting Polyadenylation Sites from DNA Sequences, Bulk RNA-seq, and Single-cell RNA-seq

Alternative polyadenylation (APA) plays important roles in modulating mRNA stability, translation, and subcellular localization, and contributes extensively to shaping eukaryotic transcriptome complexity and proteome diversity. Identification of poly(A) sites (pAs) on a genome-wide scale is a critical step toward understanding the underlying mechanism of APA-mediated gene regulation. A number of established computational tools have been proposed to predict pAs from diverse genomic data. Here we provided an exhaustive overview of computational approaches for predicting pAs from DNA sequences, bulk RNA sequencing (RNA-seq) data, and single-cell RNA sequencing (scRNA-seq) data. Particularly, we examined several representative tools using bulk RNA-seq and scRNA-seq data from peripheral blood mononuclear cells and put forward operable suggestions on how to assess the reliability of pAs predicted by different tools. We also proposed practical guidelines on choosing appropriate methods applicable to diverse scenarios. Moreover, we discussed in depth the challenges in improving the performance of pA prediction and benchmarking different methods. Additionally, we highlighted outstanding challenges and opportunities using new machine learning and integrative multi-omics techniques, and provided our perspective on how computational methodologies might evolve in the future for non-3' untranslated region, tissue-specific, cross-species, and single-cell pA prediction.

Enhancers regulate 3' end processing activity to control expression of alternative 3'UTR isoforms

Multi-UTR genes are widely transcribed and express their alternative 3'UTR isoforms in a cell type-specific manner. As transcriptional enhancers regulate mRNA expression, we investigated if they also regulate 3'UTR isoform expression. Endogenous enhancer deletion of the multi-UTR gene PTEN did not impair transcript production but prevented 3'UTR isoform switching which was recapitulated by silencing of an enhancer-bound transcription factor. In reporter assays, enhancers increase transcript production when paired with single-UTR gene promoters. However, when combined with multi-UTR gene promoters, they change 3'UTR isoform expression by increasing 3' end processing activity of polyadenylation sites. Processing activity of polyadenylation sites is affected by transcription factors, including NF-κB and MYC, transcription elongation factors, chromatin remodelers, and histone acetyltransferases. As endogenous cell type-specific enhancers are associated with genes that increase their short 3'UTRs in a cell type-specific manner, our data suggest that transcriptional enhancers integrate cellular signals to regulate cell type-and condition-specific 3'UTR isoform expression.

Regulation of the Alternative Neural Transcriptome by ELAV/Hu RNA Binding Proteins

The process of alternative polyadenylation (APA) generates multiple 3' UTR isoforms for a given locus, which can alter regulatory capacity and on occasion change coding potential. APA was initially characterized for a few genes, but in the past decade, has been found to be the rule for metazoan genes. While numerous differences in APA profiles have been catalogued across genetic conditions, perturbations, and diseases, our knowledge of APA mechanisms and biology is far from complete. In this review, we highlight recent findings regarding the role of the conserved ELAV/Hu family of RNA binding proteins (RBPs) in generating the broad landscape of lengthened 3' UTRs that is characteristic of neurons. We relate this to their established roles in alternative splicing, and summarize ongoing directions that will further elucidate the molecular strategies for neural APA, the in vivo functions of ELAV/Hu RBPs, and the phenotypic consequences of these regulatory paradigms in neurons.

The Role of Alternative Polyadenylation in the Regulation of Subcellular RNA Localization

Alternative polyadenylation (APA) is a widespread and conserved regulatory mechanism that generates diverse 3' ends on mRNA. APA patterns are often tissue specific and play an important role in cellular processes such as cell proliferation, differentiation, and response to stress. Many APA sites are found in 3' UTRs, generating mRNA isoforms with different 3' UTR contents. These alternate 3' UTR isoforms can change how the transcript is regulated, affecting its stability and translation. Since the subcellular localization of a transcript is often regulated by 3' UTR sequences, this implies that APA can also change transcript location. However, this connection between APA and RNA localization has only recently been explored. In this review, we discuss the role of APA in mRNA localization across distinct subcellular compartments. We also discuss current challenges and future advancements that will aid our understanding of how APA affects RNA localization and molecular mechanisms that drive these processes.

A working model for condensate RNA-binding proteins as matchmakers for protein complex assembly

Most cellular processes are carried out by protein complexes, but it is still largely unknown how the subunits of lowly expressed complexes find each other in the crowded cellular environment. Here, we will describe a working model where RNA-binding proteins in cytoplasmic condensates act as matchmakers between their bound proteins (called protein targets) and newly translated proteins of their RNA targets to promote their assembly into complexes. Different RNA-binding proteins act as scaffolds for various cytoplasmic condensates with several of them supporting translation. mRNAs and proteins are recruited into the cytoplasmic condensates through binding to specific domains in the RNA-binding proteins. Scaffold RNA-binding proteins have a high valency. In our model, they use homotypic interactions to assemble condensates and they use heterotypic interactions to recruit protein targets into the condensates. We propose that unoccupied binding sites in the scaffold RNA-binding proteins transiently retain recruited and newly translated proteins in the condensates, thus promoting their assembly into complexes. Taken together, we propose that lowly expressed subunits of protein complexes combine information in their mRNAs and proteins to colocalize in the cytoplasm. The efficiency of protein complex assembly is increased by transient entrapment accomplished by multivalent RNA-binding proteins within cytoplasmic condensates.

The Fine Art of Writing a Message: RNA Metabolism in the Shaping and Remodeling of the Nervous System

Neuronal morphogenesis, integration into circuits, and remodeling of synaptic connections occur in temporally and spatially defined steps. Accordingly, the expression of proteins and specific protein isoforms that contribute to these processes must be controlled quantitatively in time and space. A wide variety of post-transcriptional regulatory mechanisms, which act on pre-mRNA and mRNA molecules contribute to this control. They are thereby critically involved in physiological and pathophysiological nervous system development, function, and maintenance. Here, we review recent findings on how mRNA metabolism contributes to neuronal development, from neural stem cell maintenance to synapse specification, with a particular focus on axon growth, guidance, branching, and synapse formation. We emphasize the role of RNA-binding proteins, and highlight their emerging roles in the poorly understood molecular processes of RNA editing, alternative polyadenylation, and temporal control of splicing, while also discussing alternative splicing, RNA localization, and local translation. We illustrate with the example of the evolutionary conserved Musashi protein family how individual RNA-binding proteins are, on the one hand, acting in different processes of RNA metabolism, and, on the other hand, impacting multiple steps in neuronal development and circuit formation. Finally, we provide links to diseases that have been associated with the malfunction of RNA-binding proteins and disrupted post-transcriptional regulation.

Chromatin-contact atlas reveals disorder-mediated protein interactions and moonlighting chromatin-associated RBPs

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.

Repression of the Hox gene abd-A by ELAV-mediated Transcriptional Interference

Intergenic transcription is a common feature of eukaryotic genomes and performs important and diverse cellular functions. Here, we investigate the iab-8 ncRNA from the Drosophila Bithorax Complex and show that this RNA is able to repress the transcription of genes located at its 3’ end by a sequence-independent, transcriptional interference mechanism. Although this RNA is expressed in the early epidermis and CNS, we find that its repressive activity is limited to the CNS, where, in wild-type embryos, it acts on the Hox gene, abd-A, located immediately downstream of it. The CNS specificity is achieved through a 3’ extension of the transcript, mediated by the neuronal-specific, RNA-binding protein, ELAV. Loss of ELAV activity eliminates the 3’ extension and results in the ectopic activation of abd-A. Thus, a tissue-specific change in the length of a ncRNA is used to generate a precise pattern of gene expression in a higher eukaryote.

Although all of the cells making up complex organisms contain the same genetic material, they are nevertheless able to create the diverse tissues of the body. They do this by changing the genes they express. Thus, understanding how genes are controlled in a tissue-specific fashion is one of the primary interests of molecular genetics. Within the bithorax homeotic complex of the fruit fly Drosophila melanogaster, we, and others, previously showed that a >92 kb-long non-coding RNA, called the iab-8 ncRNA, downregulates many important developmental genes, including its genomic downstream neighbor, the homeotic gene abd-A. This downregulation is important as its loss is linked to female sterility. Interestingly, we find that the iab-8 ncRNA regulates abd-A through a mechanism called transcriptional interference, where one gene downregulates a target gene by transcribing over it. In the case of iab-8, this process is limited to the posterior central nervous system, where the iab-8 ncRNA is specifically extended into the abd-A gene by the action of the neuronal-specific RNA binding protein, ELAV. Overall, our work highlights a largely unexplored mechanism by which tissue-specific gene regulation is achieved.

Repression of the Hox gene abd-A by ELAV-mediated Transcriptional Interference.. [DOI: 10.1101/2021.09.29.462302]

Intergenic transcription is a common feature of eukaryotic genomes and performs important and diverse cellular functions. Here, we investigate the iab-8 ncRNA from the Drosophila Bithorax Complex and show that this RNA is able to repress the transcription of genes located at its 3’ end by a sequence-independent, transcriptional interference mechanism. Although this RNA is expressed in the early epidermis and CNS, we find that its repressive activity is limited to the CNS, where in wild-type embryos, it acts on the Hox gene, abd-A located immediately downstream of it. The CNS specificity is achieved through a 3’ extension of the transcript, mediated by the neuronal-specific, RNA-binding protein, ELAV. Loss of ELAV activity eliminates the 3’ extension and results in the ectopic activation of abd-A. Thus, a tissue-specific change in the length of a ncRNA is used to generate a precise pattern of gene expression in a higher eukaryote.

The landscape of alternative polyadenylation in single cells of the developing mouse embryo

3′ untranslated regions (3′ UTRs) post-transcriptionally regulate mRNA stability, localization, and translation rate. While 3′-UTR isoforms have been globally quantified in limited cell types using bulk measurements, their differential usage among cell types during mammalian development remains poorly characterized. In this study, we examine a dataset comprising ~2 million nuclei spanning E9.5–E13.5 of mouse embryonic development to quantify transcriptome-wide changes in alternative polyadenylation (APA). We observe a global lengthening of 3′ UTRs across embryonic stages in all cell types, although we detect shorter 3′ UTRs in hematopoietic lineages and longer 3′ UTRs in neuronal cell types within each stage. An analysis of RNA-binding protein (RBP) dynamics identifies ELAV-like family members, which are concomitantly induced in neuronal lineages and developmental stages experiencing 3′-UTR lengthening, as putative regulators of APA. By measuring 3′-UTR isoforms in an expansive single cell dataset, our work provides a transcriptome-wide and organism-wide map of the dynamic landscape of alternative polyadenylation during mammalian organogenesis.

Alternative polyadenylation regulates localization, half-life and translation of mRNA isoforms. Here the authors investigate alternative polyadenylation using single cell RNA sequencing data from mouse embryos and identify 3’-UTR isoforms that are regulated across cell types and developmental time.

Human retroviral antisense mRNAs are retained in the nuclei of infected cells for viral persistence

Human retroviruses, including human T cell leukemia virus type 1 (HTLV-1) and HIV type 1 (HIV-1), encode an antisense gene in the negative strand of the provirus. Besides coding for proteins, the messenger RNAs (mRNAs) of retroviral antisense genes have also been found to regulate transcription directly. Thus, it has been proposed that retroviruses likely localize their antisense mRNAs to the nucleus in order to regulate nuclear events; however, this opposes the coding function of retroviral antisense mRNAs that requires a cytoplasmic localization for protein translation. Here, we provide direct evidence that retroviral antisense mRNAs are localized predominantly in the nuclei of infected cells. The retroviral 3' LTR induces inefficient polyadenylation and nuclear retention of antisense mRNA. We further reveal that retroviral antisense RNAs retained in the nucleus associate with chromatin and have transcriptional regulatory function. While HTLV-1 antisense mRNA is recruited to the promoter of C-C chemokine receptor type 4 (CCR4) and enhances transcription from it to support the proliferation of HTLV-1-infected cells, HIV-1 antisense mRNA is recruited to the viral LTR and inhibits sense mRNA expression to maintain the latency of HIV-1 infection. In summary, retroviral antisense mRNAs are retained in nucleus, act like long noncoding RNAs instead of mRNAs, and contribute to viral persistence.

On the function and relevance of alternative 3'-UTRs in gene expression regulation

Messanger RNA (mRNA) isoforms with alternative 3'-untranslated regions (3'-UTRs) are produced by alternative polyadenylation (APA), which occurs during transcription in most eukaryotic genes. APA fine-tunes gene expression in a cell-type- and cellular state-dependent manner. Selection of an APA site entails the binding of core cleavage and polyadenylation factors to a particular polyadenylation site localized in the pre-mRNA and is controlled by multiple regulatory determinants, including transcription, pre-mRNA cis-regulatory sequences, and protein factors. Alternative 3'-UTRs serve as platforms for specific RNA binding proteins and microRNAs, which regulate gene expression in a coordinated manner by controlling mRNA fate and function in the cell. Genome-wide studies illustrated the full extent of APA prevalence and revealed that specific 3'-UTR profiles are associated with particular cellular states and diseases. Generally, short 3'-UTRs are associated with proliferative and cancer cells, and long 3'-UTRs are mostly found in polarized and differentiated cells. Fundamental new insights on the physiological consequences of this widespread event and the molecular mechanisms involved have been revealed through single-cell studies. Publicly available comprehensive databases that cover all APA mRNA isoforms identified in many cellular states and diseases reveal specific APA signatures. Therapies tackling APA mRNA isoforms or APA regulators may be regarded as innovative and attractive tools for diagnostics or treatment of several pathologies. We highlight the function of APA and alternative 3'-UTRs in gene expression regulation, the control of these mechanisms, their physiological consequences, and their potential use as new biomarkers and therapeutic tools. This article is categorized under: RNA Processing > 3' End Processing RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA in Disease and Development > RNA in Disease.

Cis-acting variation is common across regulatory layers but is often buffered during embryonic development

Precise patterns of gene expression are driven by interactions between transcription factors, regulatory DNA sequences, and chromatin. How DNA mutations affecting any one of these regulatory "layers" are buffered or propagated to gene expression remains unclear. To address this, we quantified allele-specific changes in chromatin accessibility, histone modifications, and gene expression in F1 embryos generated from eight Drosophila crosses at three embryonic stages, yielding a comprehensive data set of 240 samples spanning multiple regulatory layers. Genetic variation (allelic imbalance) impacts gene expression more frequently than chromatin features, with metabolic and environmental response genes being most often affected. Allelic imbalance in cis-regulatory elements (enhancers) is common and highly heritable, yet its functional impact does not generally propagate to gene expression. When it does, genetic variation impacts RNA levels through two alternative mechanisms involving either H3K4me3 or chromatin accessibility and H3K27ac. Changes in RNA are more predictive of variation in H3K4me3 than vice versa, suggesting a role for H3K4me3 downstream from transcription. The impact of a substantial proportion of genetic variation is consistent across embryonic stages, with 50% of allelic imbalanced features at one stage being also imbalanced at subsequent developmental stages. Crucially, buffering, as well as the magnitude and evolutionary impact of genetic variants, is influenced by regulatory complexity (i.e., number of enhancers regulating a gene), with transcription factors being most robust to cis-acting, but most influenced by trans-acting, variation.

Alternative polyadenylation: methods, mechanism, function, and role in cancer

Occurring in over 60% of human genes, alternative polyadenylation (APA) results in numerous transcripts with differing 3'ends, thus greatly expanding the diversity of mRNAs and of proteins derived from a single gene. As a key molecular mechanism, APA is involved in various gene regulation steps including mRNA maturation, mRNA stability, cellular RNA decay, and protein diversification. APA is frequently dysregulated in cancers leading to changes in oncogenes and tumor suppressor gene expressions. Recent studies have revealed various APA regulatory mechanisms that promote the development and progression of a number of human diseases, including cancer. Here, we provide an overview of four types of APA and their impacts on gene regulation. We focus particularly on the interaction of APA with microRNAs, RNA binding proteins and other related factors, the core pre-mRNA 3'end processing complex, and 3'UTR length change. We also describe next-generation sequencing methods and computational tools for use in poly(A) signal detection and APA repositories and databases. Finally, we summarize the current understanding of APA in cancer and provide our vision for future APA related research.

Crosstalk of promoter and terminator during RNA polymerase II transcription cycle

The transcription cycle of RNAPII is comprised of three consecutive steps; initiation, elongation and termination. It has been assumed that the initiation and termination steps occur in spatial isolation, essentially as independent events. A growing body of evidence, however, has challenged this dogma. First, factors involved in initiation and termination exhibit both a genetic and a physical interaction during transcription. Second, the initiation and termination factors have been found to occupy both ends of a transcribing gene. Third, physical interaction of initiation and termination factors occupying distal ends of a gene sometime results in the entire terminator region of a genes looping back and contact its cognate promoter, thereby forming a looped gene architecture during transcription. A logical interpretation of these findings is that the initiation and termination steps of transcription do not occur in isolation. There is extensive communication of factors occupying promoter and terminator ends of a gene during transcription cycle. This review entails a discussion of the promoter-terminator crosstalk and its implication in the context of transcription.

Toxin-like neuropeptides in the sea anemone Nematostella unravel recruitment from the nervous system to venom

The sea anemone Nematostella vectensis (Anthozoa, Cnidaria) is a powerful model for characterizing the evolution of genes functioning in venom and nervous systems. Although venom has evolved independently numerous times in animals, the evolutionary origin of many toxins remains unknown. In this work, we pinpoint an ancestral gene giving rise to a new toxin and functionally characterize both genes in the same species. Thus, we report a case of protein recruitment from the cnidarian nervous to venom system. The ShK-like1 peptide has a ShKT cysteine motif, is lethal for fish larvae and packaged into nematocysts, the cnidarian venom-producing stinging capsules. Thus, ShK-like1 is a toxic venom component. Its paralog, ShK-like2, is a neuropeptide localized to neurons and is involved in development. Both peptides exhibit similarities in their functional activities: They provoke contraction in Nematostella polyps and are toxic to fish. Because ShK-like2 but not ShK-like1 is conserved throughout sea anemone phylogeny, we conclude that the two paralogs originated due to a Nematostella-specific duplication of a ShK-like2 ancestor, a neuropeptide-encoding gene, followed by diversification and partial functional specialization. ShK-like2 is represented by two gene isoforms controlled by alternative promoters conferring regulatory flexibility throughout development. Additionally, we characterized the expression patterns of four other peptides with structural similarities to studied venom components and revealed their unexpected neuronal localization. Thus, we employed genomics, transcriptomics, and functional approaches to reveal one venom component, five neuropeptides with two different cysteine motifs, and an evolutionary pathway from nervous to venom system in Cnidaria.

Overlapping Activities of ELAV/Hu Family RNA Binding Proteins Specify the Extended Neuronal 3' UTR Landscape in Drosophila

The tissue-specific deployment of highly extended neural 3' UTR isoforms, generated by alternative polyadenylation (APA), is a broad and conserved feature of metazoan genomes. However, the factors and mechanisms that control neural APA isoforms are not well understood. Here, we show that three ELAV/Hu RNA binding proteins (Elav, Rbp9, and Fne) have similar capacities to induce a lengthened 3' UTR landscape in an ectopic setting. These factors promote accumulation of chromatin-associated, 3' UTR-extended, nascent transcripts, through inhibition of proximal polyadenylation site (PAS) usage. Notably, Elav represses an unannotated splice isoform of fne, switching the normally cytoplasmic Fne toward the nucleus in elav mutants. We use genomic profiling to reveal strong and broad loss of neural APA in elav/fne double mutant CNS, the first genetic background to largely abrogate this distinct APA signature. Overall, we demonstrate how regulatory interplay and functionally overlapping activities of neural ELAV/Hu RBPs drives the neural APA landscape.

ELAV and FNE Determine Neuronal Transcript Signatures through EXon-Activated Rescue

The production of alternative RNA variants contributes to the tissue-specific regulation of gene expression. In the animal nervous system, a systematic shift toward distal sites of transcription termination produces transcript signatures that are crucial for neuron development and function. Here, we report that, in Drosophila, the highly conserved protein ELAV globally regulates all sites of neuronal 3' end processing and directly binds to proximal polyadenylation sites of target mRNAs in vivo. We uncover an endogenous strategy of functional gene rescue that safeguards neuronal RNA signatures in an ELAV loss-of-function context. When not directly repressed by ELAV, the transcript encoding the ELAV paralog FNE acquires a mini-exon, generating a new protein able to translocate to the nucleus and rescue ELAV-mediated alternative polyadenylation and alternative splicing. We propose that exon-activated functional rescue is a more widespread mechanism that ensures robustness of processes regulated by a hierarchy, rather than redundancy, of effectors.

So close, no matter how far: multiple paths connecting transcription to mRNA translation in eukaryotes

Transcription of DNA into mRNA and translation of mRNA into proteins are two major processes underlying gene expression. Due to the distinct molecular mechanisms, timings, and locales of action, these processes are mainly considered to be independent. During the last two decades, however, multiple factors and elements were shown to coordinate transcription and translation, suggesting an intricate level of synchronization. This review discusses the molecular mechanisms that impact both processes in eukaryotic cells of different origins. The emerging global picture suggests evolutionarily conserved regulation and coordination between transcription and mRNA translation, indicating the importance of this phenomenon for the fine-tuning of gene expression and the adjustment to constantly changing conditions.

Elav-Mediated Exon Skipping and Alternative Polyadenylation of the Dscam1 Gene Are Required for Axon Outgrowth

Many metazoan genes express alternative long 3′ UTR isoforms in the nervous system, but their functions remain largely unclear. In Drosophila melanogaster, the Dscam1 gene generates short and long (Dscam1-L) 3′ UTR isoforms because of alternative polyadenylation (APA). Here, we found that the RNA-binding protein Embryonic Lethal Abnormal Visual System (Elav) impacts Dscam1 biogenesis at two levels, including regulation of long 3′ UTR biogenesis and skipping of an upstream exon (exon 19). MinION long-read sequencing confirmed the connectivity of this alternative splicing event to the long 3′ UTR. Knockdown or CRISPR deletion of Dscam1-L impaired axon outgrowth in Drosophila. The Dscam1 long 3′ UTR was found to be required for correct Elav-mediated skipping of exon 19. Elav thus co-regulates APA and alternative splicing to generate specific Dscam1 transcripts that are essential for neural development. This coupling of APA to alternative splicing might represent a new class of regulated RNA processing.

Like most metazoan genes, Dscam1 expresses alternative short and long 3′ UTR mRNAs. Zhang et al. find that loss of Dscam1 long 3′ UTR transcripts impairs axon outgrowth in Drosophila. Long-read sequencing reveals that these long 3′ UTR mRNAs preferentially skip an upstream exon, altering Dscam1 amino acid sequence.

Alternative polyadenylation mediates genetic regulation of gene expression

Little is known about co-transcriptional or post-transcriptional regulatory mechanisms linking noncoding variation to variation in organismal traits. To begin addressing this gap, we used 3’ Seq to study the impact of genetic variation on alternative polyadenylation (APA) in the nuclear and total mRNA fractions of 52 HapMap Yoruba human lymphoblastoid cell lines. We mapped 602 APA quantitative trait loci (apaQTLs) at 10% FDR, of which 152 were nuclear specific. Effect sizes at intronic apaQTLs are negatively correlated with eQTL effect sizes. These observations suggest genetic variants can decrease mRNA expression levels by increasing usage of intronic PAS. We also identified 24 apaQTLs associated with protein levels, but not mRNA expression. Finally, we found that 19% of apaQTLs can be associated with disease. Thus, our work demonstrates that APA links genetic variation to variation in gene expression, protein expression, and disease risk, and reveals uncharted modes of genetic regulation.

3'UTRs Regulate Mouse Ntrk2 mRNA Distribution in Cortical Neurons

There are two major isoforms of NTRK2 (neurotrophic receptor tyrosine kinase 2, or TrkB), full-length isoform with tyrosine kinase (TK) domain intact (+) and spliced isoform without tyrosine kinase domain (TK(−)). Within each isoform, there exist subtypes with minor modifications of the protein sequences. In human, the NTRK2 mRNA transcripts encoding TK(+) have same 3′UTRs, while the transcripts encoding subtypes of NTRK2 TK(−) have two completely different 3′UTRs. In mouse, the mRNA transcripts encoding same NTRK2 protein sequence for either TK(+) or TK(−) have long or short 3′UTRs, respectively. The physiological functions of these different 3′UTRs are still unknown. Pilocarpine stimulation increased Ntrk2 mRNA levels in soma, while the increase in synaptosome was smaller. FISH results further showed that mouse Ntrk2 transcripts with different 3′UTRs were distributed differently in cultured cortical neurons. The transcripts with long 3′UTR were distributed more in apical dendrites compared with transcripts with short 3′UTR. Our results provide evidence of non-coding 3′UTR function in regulating mRNA distribution in neurons.

Emerging Roles for 3' UTRs in Neurons

The 3′ untranslated regions (3′ UTRs) of mRNAs serve as hubs for post-transcriptional control as the targets of microRNAs (miRNAs) and RNA-binding proteins (RBPs). Sequences in 3′ UTRs confer alterations in mRNA stability, direct mRNA localization to subcellular regions, and impart translational control. Thousands of mRNAs are localized to subcellular compartments in neurons—including axons, dendrites, and synapses—where they are thought to undergo local translation. Despite an established role for 3′ UTR sequences in imparting mRNA localization in neurons, the specific RNA sequences and structural features at play remain poorly understood. The nervous system selectively expresses longer 3′ UTR isoforms via alternative polyadenylation (APA). The regulation of APA in neurons and the neuronal functions of longer 3′ UTR mRNA isoforms are starting to be uncovered. Surprising roles for 3′ UTRs are emerging beyond the regulation of protein synthesis and include roles as RBP delivery scaffolds and regulators of alternative splicing. Evidence is also emerging that 3′ UTRs can be cleaved, leading to stable, isolated 3′ UTR fragments which are of unknown function. Mutations in 3′ UTRs are implicated in several neurological disorders—more studies are needed to uncover how these mutations impact gene regulation and what is their relationship to disease severity.

Neuronal upregulation of Prospero protein is driven by alternative mRNA polyadenylation and Syncrip-mediated mRNA stabilisation

During Drosophila and vertebrate brain development, the conserved transcription factor Prospero/Prox1 is an important regulator of the transition between proliferation and differentiation. Prospero level is low in neural stem cells and their immediate progeny, but is upregulated in larval neurons and it is unknown how this process is controlled. Here, we use single molecule fluorescent in situ hybridisation to show that larval neurons selectively transcribe a long prospero mRNA isoform containing a 15 kb 3' untranslated region, which is bound in the brain by the conserved RNA-binding protein Syncrip/hnRNPQ. Syncrip binding increases the stability of the long prospero mRNA isoform, which allows an upregulation of Prospero protein production. Adult flies selectively lacking the long prospero isoform show abnormal behaviour that could result from impaired locomotor or neurological activity. Our findings highlight a regulatory strategy involving alternative polyadenylation followed by differential post-transcriptional regulation.This article has an associated First Person interview with the first author of the paper.

Full-length transcriptome reconstruction reveals a large diversity of RNA and protein isoforms in rat hippocampus

Gene annotation is a critical resource in genomics research. Many computational approaches have been developed to assemble transcriptomes based on high-throughput short-read sequencing, however, only with limited accuracy. Here, we combine next-generation and third-generation sequencing to reconstruct a full-length transcriptome in the rat hippocampus, which is further validated using independent 5´ and 3´-end profiling approaches. In total, we detect 28,268 full-length transcripts (FLTs), covering 6,380 RefSeq genes and 849 unannotated loci. Based on these FLTs, we discover co-occurring alternative RNA processing events. Integrating with polysome profiling and ribosome footprinting data, we predict isoform-specific translational status and reconstruct an open reading frame (ORF)-eome. Notably, a high proportion of the predicted ORFs are validated by mass spectrometry-based proteomics. Moreover, we identify isoforms with subcellular localization pattern in neurons. Collectively, our data advance our knowledge of RNA and protein isoform diversity in the rat brain and provide a rich resource for functional studies.

Alternative polyadenylation of mRNA and its role in cancer

Alternative polyadenylation (APA) is a molecular process that generates diversity at the 3′ end of RNA polymerase II transcripts from over 60% of human genes. APA is derived from the existence of multiple polyadenylation signals (PAS) within the same transcript, and results in the differential inclusion of sequence information at the 3′ end. While APA can occur between two PASs allowing for generation of transcripts with distinct coding potential from a single gene, most APA occurs within the untranslated region (3′UTR) and changes the length and content of these non-coding sequences. APA within the 3′UTR can have tremendous impact on its regulatory potential of the mRNA through a variety of mechanisms, and indeed this layer of gene expression regulation has profound impact on processes vital to cell growth and development. Recent studies have particularly highlighted the importance of APA dysregulation in cancer onset and progression. Here, we review the current knowledge of APA and its impacts on mRNA stability, translation, localization and protein localization. We also discuss the implications of APA dysregulation in cancer research and therapy.

What Are 3' UTRs Doing?

3' untranslated regions (3' UTRs) of messenger RNAs (mRNAs) are best known to regulate mRNA-based processes, such as mRNA localization, mRNA stability, and translation. In addition, 3' UTRs can establish 3' UTR-mediated protein-protein interactions (PPIs), and thus can transmit genetic information encoded in 3' UTRs to proteins. This function has been shown to regulate diverse protein features, including protein complex formation or posttranslational modifications, but is also expected to alter protein conformations. Therefore, 3' UTR-mediated information transfer can regulate protein features that are not encoded in the amino acid sequence. This review summarizes both 3' UTR functions-the regulation of mRNA and protein-based processes-and highlights how each 3' UTR function was discovered with a focus on experimental approaches used and the concepts that were learned. This review also discusses novel approaches to study 3' UTR functions in the future by taking advantage of recent advances in technology.

Targeting mRNA processing as an anticancer strategy

Discoveries in the past decade have highlighted the potential of mRNA as a therapeutic target for cancer. Specifically, RNA sequencing revealed that, in addition to gene mutations, alterations in mRNA can contribute to the initiation and progression of cancer. Indeed, precursor mRNA processing, which includes the removal of introns by splicing and the formation of 3' ends by cleavage and polyadenylation, is frequently altered in tumours. These alterations result in numerous cancer-specific mRNAs that generate altered levels of normal proteins or proteins with new functions, leading to the activation of oncogenes or the inactivation of tumour-suppressor genes. Abnormally spliced and polyadenylated mRNAs are also associated with resistance to cancer treatment and, unexpectedly, certain cancers are highly sensitive to the pharmacological inhibition of splicing. This Review summarizes recent progress in our understanding of how splicing and polyadenylation are altered in cancer and highlights how this knowledge has been translated for drug discovery, resulting in the production of small molecules and oligonucleotides that modulate the spliceosome and are in clinical trials for the treatment of cancer.

FLAM-seq: full-length mRNA sequencing reveals principles of poly(A) tail length control

Although messenger RNAs are key molecules for understanding life, until now, no method has existed to determine the full-length sequence of endogenous mRNAs including their poly(A) tails. Moreover, although non-A nucleotides can be incorporated in poly(A) tails, there also exists no method to accurately sequence them. Here, we present full-length poly(A) and mRNA sequencing (FLAM-seq), a rapid and simple method for high-quality sequencing of entire mRNAs. We report a complementary DNA library preparation method coupled to single-molecule sequencing to perform FLAM-seq. Using human cell lines, brain organoids and Caenorhabditis elegans we show that FLAM-seq delivers high-quality full-length mRNA sequences for thousands of different genes per sample. We find that 3' untranslated region length is correlated with poly(A) tail length, that alternative polyadenylation sites and alternative promoters for the same gene are linked to different tail lengths, and that tails contain a substantial number of cytosines.

The Length of the Expressed 3' UTR Is an Intermediate Molecular Phenotype Linking Genetic Variants to Complex Diseases

In the last decades, genome-wide association studies (GWAS) have uncovered tens of thousands of associations between common genetic variants and complex diseases. However, these statistical associations can rarely be interpreted functionally and mechanistically. As the majority of the disease-associated variants are located far from coding sequences, even the relevant gene is often unclear. A way to gain insight into the relevant mechanisms is to study the genetic determinants of intermediate molecular phenotypes, such as gene expression and transcript structure. We propose a computational strategy to discover genetic variants affecting the relative expression of alternative 3′ untranslated region (UTR) isoforms, generated through alternative polyadenylation, a widespread posttranscriptional regulatory mechanism known to have relevant functional consequences. When applied to a large dataset in which whole genome and RNA sequencing data are available for 373 European individuals, 2,530 genes with alternative polyadenylation quantitative trait loci (apaQTL) were identified. We analyze and discuss possible mechanisms of action of these variants, and we show that they are significantly enriched in GWAS hits, in particular those concerning immune-related and neurological disorders. Our results point to an important role for genetically determined alternative polyadenylation in affecting predisposition to complex diseases, and suggest new ways to extract functional information from GWAS data.

FLAM-seq: full-length mRNA sequencing reveals principles of poly(A) tail length control

Although messenger RNAs are key molecules for understanding life, until now, no method has existed to determine the full-length sequence of endogenous mRNAs including their poly(A) tails. Moreover, although non-A nucleotides can be incorporated in poly(A) tails, there also exists no method to accurately sequence them. Here, we present full-length poly(A) and mRNA sequencing (FLAM-seq), a rapid and simple method for high-quality sequencing of entire mRNAs. We report a complementary DNA library preparation method coupled to single-molecule sequencing to perform FLAM-seq. Using human cell lines, brain organoids and Caenorhabditis elegans we show that FLAM-seq delivers high-quality full-length mRNA sequences for thousands of different genes per sample. We find that 3' untranslated region length is correlated with poly(A) tail length, that alternative polyadenylation sites and alternative promoters for the same gene are linked to different tail lengths, and that tails contain a substantial number of cytosines.

A survey on identification and quantification of alternative polyadenylation sites from RNA-seq data

Alternative polyadenylation (APA) has been implicated to play an important role in post-transcriptional regulation by regulating mRNA abundance, stability, localization and translation, which contributes considerably to transcriptome diversity and gene expression regulation. RNA-seq has become a routine approach for transcriptome profiling, generating unprecedented data that could be used to identify and quantify APA site usage. A number of computational approaches for identifying APA sites and/or dynamic APA events from RNA-seq data have emerged in the literature, which provide valuable yet preliminary results that should be refined to yield credible guidelines for the scientific community. In this review, we provided a comprehensive overview of the status of currently available computational approaches. We also conducted objective benchmarking analysis using RNA-seq data sets from different species (human, mouse and Arabidopsis) and simulated data sets to present a systematic evaluation of 11 representative methods. Our benchmarking study showed that the overall performance of all tools investigated is moderate, reflecting that there is still lot of scope to improve the prediction of APA site or dynamic APA events from RNA-seq data. Particularly, prediction results from individual tools differ considerably, and only a limited number of predicted APA sites or genes are common among different tools. Accordingly, we attempted to give some advice on how to assess the reliability of the obtained results. We also proposed practical recommendations on the appropriate method applicable to diverse scenarios and discussed implications and future directions relevant to profiling APA from RNA-seq data.

Pervasive Chromatin-RNA Binding Protein Interactions Enable RNA-Based Regulation of Transcription

Increasing evidence suggests that transcriptional control and chromatin activities at large involve regulatory RNAs, which likely enlist specific RNA-binding proteins (RBPs). Although multiple RBPs have been implicated in transcription control, it has remained unclear how extensively RBPs directly act on chromatin. We embarked on a large-scale RBP ChIP-seq analysis, revealing widespread RBP presence in active chromatin regions in the human genome. Like transcription factors (TFs), RBPs also show strong preference for hotspots in the genome, particularly gene promoters, where their association is frequently linked to transcriptional output. Unsupervised clustering reveals extensive co-association between TFs and RBPs, as exemplified by YY1, a known RNA-dependent TF, and RBM25, an RBP involved in splicing regulation. Remarkably, RBM25 depletion attenuates all YY1-dependent activities, including chromatin binding, DNA looping, and transcription. We propose that various RBPs may enhance network interaction through harnessing regulatory RNAs to control transcription.

APAtrap: identification and quantification of alternative polyadenylation sites from RNA-seq data

Motivation

Alternative polyadenylation (APA) has been increasingly recognized as a crucial mechanism that contributes to transcriptome diversity and gene expression regulation. As RNA-seq has become a routine protocol for transcriptome analysis, it is of great interest to leverage such unprecedented collection of RNA-seq data by new computational methods to extract and quantify APA dynamics in these transcriptomes. However, research progress in this area has been relatively limited. Conventional methods rely on either transcript assembly to determine transcript 3' ends or annotated poly(A) sites. Moreover, they can neither identify more than two poly(A) sites in a gene nor detect dynamic APA site usage considering more than two poly(A) sites.

Results

We developed an approach called APAtrap based on the mean squared error model to identify and quantify APA sites from RNA-seq data. APAtrap is capable of identifying novel 3' UTRs and 3' UTR extensions, which contributes to locating potential poly(A) sites in previously overlooked regions and improving genome annotations. APAtrap also aims to tally all potential poly(A) sites and detect genes with differential APA site usages between conditions. Extensive comparisons of APAtrap with two other latest methods, ChangePoint and DaPars, using various RNA-seq datasets from simulation studies, human and Arabidopsis demonstrate the efficacy and flexibility of APAtrap for any organisms with an annotated genome.

Availability and implementation

Freely available for download at https://apatrap.sourceforge.io.

Contact

liqq@xmu.edu.cn or xhuister@xmu.edu.cn.

Supplementary information

Supplementary data are available at Bioinformatics online.