Control of poly(A)-tail length and translation in vertebrate oocytes and early embryos

During oocyte maturation and early embryogenesis, changes in mRNA poly(A)-tail lengths strongly influence translation, but how these tail-length changes are orchestrated has been unclear. Here, we performed tail-length and translational profiling of mRNA reporter libraries (each with millions of 3' UTR sequence variants) in frog oocytes and embryos and in fish embryos. Contrasting to previously proposed cytoplasmic polyadenylation elements (CPEs), we found that a shorter element, UUUUA, together with the polyadenylation signal (PAS), specify cytoplasmic polyadenylation, and we identified contextual features that modulate the activity of both elements. In maturing oocytes, this tail lengthening occurs against a backdrop of global deadenylation and the action of C-rich elements that specify tail-length-independent translational repression. In embryos, cytoplasmic polyadenylation becomes more permissive, and additional elements specify waves of stage-specific deadenylation. Together, these findings largely explain the complex tapestry of tail-length changes observed in early frog and fish development, with strong evidence of conservation in both mice and humans.

Rational design of an artificial tethered enzyme for non-templated post-transcriptional mRNA polyadenylation by the second generation of the C3P3 system

We have previously introduced the first generation of C3P3, an artificial system that allows the autonomous in-vivo production of mRNA with m7GpppN-cap. While C3P3-G1 synthesized much larger amounts of capped mRNA in human cells than conventional nuclear expression systems, it produced a proportionately much smaller amount of the corresponding proteins, indicating a clear defect of mRNA translatability. A possible mechanism for this poor translatability could be the rudimentary polyadenylation of the mRNA produced by the C3P3-G1 system. We therefore sought to develop the C3P3-G2 system using an artificial enzyme to post-transcriptionally lengthen the poly(A) tail. This system is based on the mutant mouse poly(A) polymerase alpha fused at its N terminus with an N peptide from the λ virus, which binds to BoxBr sequences placed in the 3'UTR region of the mRNA of interest. The resulting system selectively brings mPAPαm7 to the target mRNA to elongate its poly(A)-tail to a length of few hundred adenosine. Such elongation of the poly(A) tail leads to an increase in protein expression levels of about 2.5-3 times in cultured human cells compared to the C3P3-G1 system. Finally, the coding sequence of the tethered mutant poly(A) polymerase can be efficiently fused to that of the C3P3-G1 enzyme via an F2A sequence, thus constituting the single-ORF C3P3-G2 enzyme. These technical developments constitute an important milestone in improving the performance of the C3P3 system, paving the way for its applications in bioproduction and non-viral human gene therapy.

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.

Differential Poly(A) Tail Length Analysis Using Nanopore Sequencing

The poly(A) tail is a sequence of several adenosine nucleotides added to the 3' end of RNA molecules transcribed by polymerase II. The dynamics of poly(A) tail length play a significant role in regulating post-transcriptional gene expression by regulating the stability, translation, and decay of messenger RNAs. As a result, an accurate measurement of poly(A) tail length changes is important for understanding its regulatory function in different cellular contexts. Here, we outline a method for using nanopore sequencing and linear mixed models to analyze differences in poly(A) tail length across conditions.

Sequencing of Transcriptome-Wide Poly(A) Tails by PAIso-seq

Poly(A) tails are added to most eukaryotic mRNA and have essential regulatory functions. However, due to its homopolymeric nature, the sequence information in poly(A) tails is challenging to obtain in transcriptome measurement studies. In this chapter, we describe the detailed procedures of poly(A) inclusive full-length RNA isoform-sequencing (PAIso-seq), a method that can measure transcriptome-wide poly(A) tails from as low as nanogram level of total RNA based on the PacBio HiFi sequencing platform. The accurate length and base composition of poly(A) tails can be obtained along with the full-length cDNA.

Analysis of Human Endogenous mRNA Deadenylation Complexes by High-Resolution Gel Electrophoresis

In eukaryotic cells, poly(A) tails stabilize mRNA molecules and play a pivotal role in enhancing translational efficiency. Consequently, the enzymatic shortening of these poly(A) tails by deadenylase enzymes has a critical role in post-transcriptional gene regulation. However, deadenylases are usually large, multisubunit, and multifunctional complexes, which complicates their biochemical analysis. This chapter presents a methodology for isolating human deadenylation complexes from endogenous sources and conducting an in vitro deadenylation assay to examine their enzymatic activity. The reactions involving fluorescently labeled synthetic polyadenylated RNAs are subsequently analyzed using high-resolution denaturing polyacrylamide gel electrophoresis.

An RNA-Ligation-Based RACE-PAT Assay to Monitor Poly(A) Tail Length of mRNAs of Interest

In eukaryotes, a non-templated poly-adenosine (poly(A)) tail is added co-transcriptionally to almost every messenger RNA (mRNA). The length of this poly(A) tail changes during the lifetime of mRNAs and has been shown in many circumstances to be an important factor controlling transcript fates. Yet, the measure of the length of this homogenous nucleotide sequence is technically challenging, making it difficult to assess its dynamic variation. In this chapter, we describe an RNA-ligation-based RACE-PAT (Rapid Amplification of cDNA End-Poly(A) Tail) assay to monitor the poly(A) tail length of mRNAs. In the first step, an RNA oligonucleotide is ligated to mRNA 3' ends providing an anchoring site to prime cDNA synthesis, avoiding the bias introduced by oligo(dT)-derived primers. Afterward, reverse transcription is performed with an anchor primer with a unique 5' extension. The choice of the oligonucleotide 3' end at this step allows further flexibility to amplify modified tails, for example, by uridylation. Next, short DNA fragments encompassing the poly(A) tails are amplified by Polymerase Chain Reaction (PCR) using as forward primer, a transcript-specific primer hybridizing close to the transcript polyadenylation signal, and as reverse primer, an oligonucleotide corresponding to the 5' extension of the primer used for cDNA synthesis, ensuring that only cDNAs are amplified. The resulting DNA fragments are then visualized after size fractionation by electrophoresis. This method does not provide exact nucleotide count and composition but has the advantage of allowing the processing of many samples in parallel at a low cost.

Quantification of Poly(A) Tail Length and Terminal Modifications Using Direct RNA Sequencing

Poly(A) tail metabolism is critical for various biological processes, including early embryogenesis and cell differentiation. While traditional biochemical methods to measure poly(A) tail length allow for the study of selected transcripts, the advent of long-read sequencing technologies enabled the development of simple and robust protocols to measure poly(A) tail length at the transcriptome level. Here, we describe a direct RNA sequencing protocol to capture poly(A) tail terminal additions based on the splint ligation of barcoded oligos compatible with terminal guanylation and uridylation. We cover how to prepare the libraries and perform the bioinformatics analysis to simultaneously determine the length of the transcripts' poly(A) tails and detect the presence of terminal guanylation and uridylation.

Measuring Poly-Adenosine Tail Length of RNAs by High-Resolution Northern Blotting Coupled with RNase H Cleavage

The poly-adenosine, or poly(A) tail, plays key roles in controlling the stability and translation of messenger RNAs in all eukaryotes, and, as such, facile assays that can measure poly(A) length are needed. This chapter describes an approach that couples RNase H-mediated cleavage of an RNA of interest with high-resolution denaturing gel electrophoresis and northern blot-based detection. Major advantages of this method include the ability to directly measure the abundance of any RNA and the length of its poly(A) tail without amplification steps. The assay provides high specificity, sensitivity, and reproducibility for accurate quantitation using standard molecular biology equipment and reagents. Overall, the high-resolution northern blotting approach offers a cost-effective means of poly(A) RNA analysis that is especially useful for small numbers of transcripts and comparisons between experimental conditions or time points.

Nano3'RACE: A Method to Analyze Poly(A) Tail Length and Nucleotide Additions at the 3' Extremity of Selected mRNAs Using Nanopore Sequencing

Deadenylation is a major process that regulates gene expression by shaping the length of mRNA poly(A) tails. Deadenylation is controlled by factors in trans that recruit or impede deadenylases, by the incorporation of non-adenosines during poly(A) tail synthesis, and by the posttranscriptional addition of 3' nucleotides to poly(A) tails. Deciphering the regulation of poly(A) tail shortening requires both transcriptome-wide approaches and more targeted methodologies, allowing deep analyses of specific mRNAs. In this chapter, we present Nano3'RACE, a nanopore-based cDNA sequencing method that allows in-depth analysis to precisely measure poly(A) tail length and detect 3' terminal nucleotide addition, such as uridylation, for mRNAs of interest.

Single-Molecule Poly(A) Tail Sequencing (SM-PATseq) Using the PacBio Platform

The polyadenylation of the 3' ends of messenger RNAs is an important regulator of stability and translation. We developed the single-molecule poly(A) tail sequencing method, SM-PATseq, to assay tail lengths of the whole transcriptome at nucleotide resolution using long-read sequencing. This method generates cDNA using an oligo-dT 3' splint adaptor ligation to prime first-strand cDNA synthesis, followed by random hexamer priming for second-strand synthesis. By directly sequencing the cDNA on long-read platforms, we can resolve tail lengths at nucleotide resolution, identify non-A bases within the tail, and quantify transcript abundance analogous to traditional RNAseq methods. Here, we discuss the method for generating, sequencing, and primary analysis of poly(A) tail data from total RNA using the Pacific Biosciences Sequel platform.

Profiling the polyadenylated transcriptome of extracellular vesicles with long-read nanopore sequencing

BACKGROUND

While numerous studies have described the transcriptomes of extracellular vesicles (EVs) in different cellular contexts, these efforts have typically relied on sequencing methods requiring RNA fragmentation, which limits interpretations on the integrity and isoform diversity of EV-targeted RNA populations. It has been assumed that mRNA signatures in EVs are likely to be fragmentation products of the cellular mRNA material, and the extent to which full-length mRNAs are present within EVs remains to be clarified.

RESULTS

Using long-read nanopore RNA sequencing, we sought to characterize the full-length polyadenylated (poly-A) transcriptome of EVs released by human chronic myelogenous leukemia K562 cells. We detected 443 and 280 RNAs that were respectively enriched or depleted in EVs. EV-enriched poly-A transcripts consist of a variety of biotypes, including mRNAs, long non-coding RNAs, and pseudogenes. Our analysis revealed that 10.58% of all EV reads, and 18.67% of all cellular (WC) reads, corresponded to known full-length transcripts, with mRNAs representing the largest biotype for each group (EV = 58.13%, WC = 43.93%). We also observed that for many well-represented coding and non-coding genes, diverse full-length transcript isoforms were present in EV specimens, and these isoforms were reflective-of but often in different ratio compared to cellular samples.

CONCLUSION

This work provides novel insights into the compositional diversity of poly-A transcript isoforms enriched within EVs, while also underscoring the potential usefulness of nanopore sequencing to interrogate secreted RNA transcriptomes.

Only fourteen 3'-end poly(A)s sufficient for rescuing Senecavirus A from its cDNA clone, but inadequate to meet requirement of viral replication

Senecavirus A (SVA) belongs to the genus Senecavirus in the family Picornaviridae. Its genome is a positive-sense, single-strand RNA that has 5' and 3' untranslated regions. There is a poly(A) tail at the 3' end of viral genome. Although the number of poly(A)s is variable, the length of poly(A) tail generally has the minimum nucleotide limit for picornaviral replication. To identify a range limit of poly(A)s for SVA recovery, five SVA cDNA clones, separately containing 25, 20, 15, 10 and 5 poly(A)s, were constructed for rescuing viruses. Replication-competent SVAs could be rescued from the first three cDNA clones, implying the range limit of poly(A)s was (A)₁₅ to (A)₁₀. To recognize the precise limit, four extra cDNA clones, separately containing 14, 13, 12 and 11 poly(A)s, were constructed to rescue SVAs further. The replication-competent SVA was rescued only from the poly(A)₁₄-containing plasmid, indicating that the precise limit was poly(A)₁₄ at the 3' end of cDNA clone for SVA recovery. The rescued SVA was serially passaged in cells. The passage-5 and -10 progenies were independently subjected to the analysis of 3'-rapid amplification of cDNA ends. Both progenies showed their own poly(A) tails far more than 14 (A)s, implying extra (A)s added to the poly(A)₁₄ sequence during viral passaging. It can be concluded that fourteen (A)s are sufficient for rescuing a replication-competent SVA from its cDNA clone, but inadequate for maintaining viral propagation in cells.

Nano3P-seq: transcriptome-wide analysis of gene expression and tail dynamics using end-capture nanopore cDNA sequencing

RNA polyadenylation plays a central role in RNA maturation, fate, and stability. In response to developmental cues, polyA tail lengths can vary, affecting the translation efficiency and stability of mRNAs. Here we develop Nanopore 3' end-capture sequencing (Nano3P-seq), a method that relies on nanopore cDNA sequencing to simultaneously quantify RNA abundance, tail composition, and tail length dynamics at per-read resolution. By employing a template-switching-based sequencing protocol, Nano3P-seq can sequence RNA molecule from its 3' end, regardless of its polyadenylation status, without the need for PCR amplification or ligation of RNA adapters. We demonstrate that Nano3P-seq provides quantitative estimates of RNA abundance and tail lengths, and captures a wide diversity of RNA biotypes. We find that, in addition to mRNA and long non-coding RNA, polyA tails can be identified in 16S mitochondrial ribosomal RNA in both mouse and zebrafish models. Moreover, we show that mRNA tail lengths are dynamically regulated during vertebrate embryogenesis at an isoform-specific level, correlating with mRNA decay. Finally, we demonstrate the ability of Nano3P-seq in capturing non-A bases within polyA tails of various lengths, and reveal their distribution during vertebrate embryogenesis. Overall, Nano3P-seq is a simple and robust method for accurately estimating transcript levels, tail lengths, and tail composition heterogeneity in individual reads, with minimal library preparation biases, both in the coding and non-coding transcriptome.

Review: Mechanisms underlying alternative polyadenylation in plants - looking in the right places

Recent years have seen an explosion of interest in the subject of alternative polyadenylation in plants. Connections between the polyadenylation complex and numerous developmental and stress responses are well-established. However, those that link stimuli with the functioning of the polyadenylation complex are less well understood. To this end, it is imperative to clearly delineate the roles of the polyadenylation complex in both plant growth AND alternative polyadenylation. It is also necessary to understand the ways by which other molecular processes may contribute to alternative polyadenylation. This review discusses these issues, with a focus on instances that reveal mechanisms by which mRNA polyadenylation may be regulated. Insights from from characterizations of mutants affected in the polyadenylation complex are discussed, as are the limitations of such characterizations when it comes to teasing out cause and effect. These limitations encourage explorations to other processes that are beyond the core polyadenylation complex. Two such processes that sculpt the plant transcriptome - transcription termination and the epigenetic control of transposon activity - also contribute to regulated poly(A) site choice. These subjects define "the right places" - molecular mechanisms that contribute to the wide-ranging control of gene expression via mRNA polyadenylation.

PABPN1 functions as a hub in the assembly of nuclear poly(A) domains that are essential for mouse oocyte development

Growing oocytes store a large amount of maternal mRNA to support the subsequent "maternal-zygotic transition" process. At present, it is not clear how the growing oocytes store and process the newly transcribed mRNA under physiological conditions. In this study, we report non-membrane-bound compartments, nuclear poly(A) domains (NPADs), as the hub for newly transcribed mRNA, in developing mouse oocytes. The RNA binding protein PABPN1 promotes the formation of NPAD through its N-terminal disordered domain and RNA-recognized motif by means of liquid phase separation. Pabpn1-null growing oocytes cannot form NPAD normally in vivo and have defects in stability of oocyte growing-related transcripts and formation of long 3' untranslated region isoform transcripts. Ultimately, Pabpn1^fl/fl;Gdf9-Cre mice are completely sterile with primary ovarian insufficiency. These results demonstrate that NPAD formed by the phase separation properties of PABPN1-mRNA are the hub of the newly transcribed mRNA and essential for the development of oocytes and female reproduction.

3' Untranslated Regions Are Modular Entities That Determine Polyadenylation Profiles

The 3' ends of eukaryotic mRNAs are generated by cleavage of nascent transcripts followed by polyadenylation, which occurs at numerous sites within 3' untranslated regions (3' UTRs) but rarely within coding regions. An individual gene can yield many 3'-mRNA isoforms with distinct half-lives. We dissect the relative contributions of protein-coding sequences (open reading frames [ORFs]) and 3' UTRs to polyadenylation profiles in yeast. ORF-deleted derivatives often display strongly decreased mRNA levels, indicating that ORFs contribute to overall mRNA stability. Poly(A) profiles, and hence relative isoform half-lives, of most (9 of 10) ORF-deleted derivatives are very similar to their wild-type counterparts. Similarly, in-frame insertion of a large protein-coding fragment between the ORF and 3' UTR has minimal effect on the poly(A) profile in all 15 cases tested. Last, reciprocal ORF/3'-UTR chimeric genes indicate that the poly(A) profile is determined by the 3' UTR. Thus, 3' UTRs are self-contained modular entities sufficient to determine poly(A) profiles and relative 3'-isoform half-lives. In the one atypical instance, ORF deletion causes an upstream shift of poly(A) sites, likely because juxtaposition of an unusually high AT-rich stretch directs polyadenylation closely downstream. This suggests that long AT-rich stretches, which are not encountered until after coding regions, are important for restricting polyadenylation to 3' UTRs.

An atlas of plant full-length RNA reveals tissue-specific and monocots-dicots conserved regulation of poly(A) tail length

Poly(A) tail is a hallmark of eukaryotic messenger RNA and its length plays an essential role in regulating mRNA metabolism. However, a comprehensive resource for plant poly(A) tail length has yet to be established. Here, we applied a poly(A)-enrichment-free, nanopore-based method to profile full-length RNA with poly(A) tail information in plants. Our atlas contains over 120 million polyadenylated mRNA molecules from seven different tissues of Arabidopsis, as well as the shoot tissue of maize, soybean and rice. In most tissues, the size of plant poly(A) tails shows peaks at approximately 20 and 45 nucleotides, while the poly(A) tails in pollen exhibit a distinct pattern with strong peaks centred at 55 and 80 nucleotides. Moreover, poly(A) tail length is regulated in a gene-specific manner-mRNAs with short half-lives in general have long poly(A) tails, while mRNAs with long half-lives are featured with relatively short poly(A) tails that peak at ~45 nucleotides. Across species, poly(A) tails in the nucleus are almost twice as long as in the cytoplasm. Our comprehensive dataset lays the groundwork for future functional and evolutionary studies on poly(A) tail length regulation in plants.

Dynamic and widespread control of poly(A) tail length during macrophage activation

The poly(A) tail enhances translation and transcript stability, and tail length is under dynamic control during cell state transitions. Tail regulation plays essential roles in translational timing and fertilization in early development, but poly(A) tail dynamics have not been fully explored in post-embryonic systems. Here, we examined the landscape and impact of tail length control during macrophage activation. Upon activation, more than 1500 mRNAs, including proinflammatory genes, underwent distinctive changes in tail lengths. Increases in tail length correlated with mRNA levels regardless of transcriptional activity, and many mRNAs that underwent tail extension encode proteins necessary for immune function and post-transcriptional regulation. Strikingly, we found that ZFP36, whose protein product destabilizes target transcripts, undergoes tail extension. Our analyses indicate that many mRNAs undergoing tail lengthening are, in turn, degraded by elevated levels of ZFP36, constituting a post-transcriptional feedback loop that ensures transient regulation of transcripts integral to macrophage activation. Taken together, this study establishes the complexity, relevance, and widespread nature of poly(A) tail dynamics, and the resulting post-transcriptional regulation during macrophage activation.

The interplay between translational efficiency, poly(A) tails, microRNAs, and neuronal activation

Neurons provide a rich setting for studying post-transcriptional control. Here, we investigate the landscape of translational control in neurons and search for mRNA features that explain differences in translational efficiency (TE), considering the interplay between TE, mRNA poly(A)-tail lengths, microRNAs, and neuronal activation. In neurons and brain tissues, TE correlates with tail length, and a few dozen mRNAs appear to undergo cytoplasmic polyadenylation upon light or chemical stimulation. However, the correlation between TE and tail length is modest, explaining <5% of TE variance, and even this modest relationship diminishes when accounting for other mRNA features. Thus, tail length appears to affect TE only minimally. Accordingly, miRNAs, which accelerate deadenylation of their mRNA targets, primarily influence target mRNA levels, with no detectable effect on either steady-state tail lengths or TE. Larger correlates with TE include codon composition and predicted mRNA folding energy. When combined in a model, the identified correlates explain 38%-45% of TE variance. These results provide a framework for considering the relative impact of factors that contribute to translational control in neurons. They indicate that when examined in bulk, translational control in neurons largely resembles that of other types of post-embryonic cells. Thus, detection of more specialized control might require analyses that can distinguish translation occurring in neuronal processes from that occurring in cell bodies.

Oppositional poly(A) tail length regulation by FMRP and CPEB1

Poly(A) tail length is regulated in both the nucleus and cytoplasm. One factor that controls polyadenylation in the cytoplasm is CPEB1, an RNA binding protein that associates with specific mRNA 3'UTR sequences to tether enzymes that add and remove poly(A). Two of these enzymes, the noncanonical poly(A) polymerases GLD2 (TENT2, PAPD4, Wispy) and GLD4 (TENT4B, PAPD5, TRF4, TUT3), interact with CPEB1 to extend poly(A). To identify additional RNA binding proteins that might anchor GLD4 to RNA, we expressed double tagged GLD4 in U87MG cells, which was used for sequential immunoprecipitation and elution followed by mass spectrometry. We identified several RNA binding proteins that coprecipitated with GLD4, among which was FMRP. To assess whether FMRP regulates polyadenylation, we performed TAIL-seq from WT and FMRP-deficient HEK293 cells. Surprisingly, loss of FMRP resulted in an overall increase in poly(A), which was also observed for several specific mRNAs. Conversely, loss of CPEB1 elicited an expected decrease in poly(A), which was examined in cultured neurons. We also examined polyadenylation in wild type (WT) and FMRP-deficient mouse brain cortex by direct RNA nanopore sequencing, which identified RNAs with both increased and decreased poly(A). Our data show that FMRP has a role in mediating poly(A) tail length, which adds to its repertoire of RNA regulation.

Precise gene models using long-read sequencing reveal a unique poly(A) signal in Giardia lamblia

During pre-mRNA processing, the poly(A) signal is recognized by a protein complex that ensures precise cleavage and polyadenylation of the nascent transcript. The location of this cleavage event establishes the length and sequence of the 3' UTR of an mRNA, thus determining much of its post-transcriptional fate. Using long-read sequencing, we characterize the polyadenylation signal and related sequences surrounding Giardia lamblia cleavage sites for over 2600 genes. We find that G. lamblia uses an AGURAA poly(A) signal, which differs from the mammalian AAUAAA. We also describe how G. lamblia lacks common auxiliary elements found in other eukaryotes, along with the proteins that recognize them. Further, we identify 133 genes with evidence of alternative polyadenylation. These results suggest that despite pared-down cleavage and polyadenylation machinery, 3' end formation still appears to be an important regulatory step for gene expression in G. lamblia.

Genome-Wide Identification of Polyadenylation Dynamics with TED-Seq

Polyadenylation and deadenylation of mRNA are major RNA modifications associated with nucleus-to-cytoplasm translocation, mRNA stability, translation efficiency, and mRNA decay pathways. Our current knowledge of polyadenylation and deadenylation has been expanded due to recent advances in transcriptome-wide poly(A) tail length assays. Whereas these methods measure poly(A) length by quantifying the adenine (A) base stretch at the 3' end of mRNA, we developed a more cost-efficient technique that does not rely on A-base counting, called tail-end-displacement sequencing (TED-seq). Through sequencing highly size-selected 3' RNA fragments including the poly(A) tail pieces, TED-seq provides accurate measure of transcriptome-wide poly(A)-tail lengths in high resolution, economically suitable for larger scale analysis under various biologically transitional contexts.

Dynamic bulge nucleotides in the KSHV PAN ENE triple helix provide a unique binding platform for small molecule ligands

Cellular and virus-coded long non-coding (lnc) RNAs support multiple roles related to biological and pathological processes. Several lncRNAs sequester their 3' termini to evade cellular degradation machinery, thereby supporting disease progression. An intramolecular triplex involving the lncRNA 3' terminus, the element for nuclear expression (ENE), stabilizes RNA transcripts and promotes persistent function. Therefore, such ENE triplexes, as presented here in Kaposi's sarcoma-associated herpesvirus (KSHV) polyadenylated nuclear (PAN) lncRNA, represent targets for therapeutic development. Towards identifying novel ligands targeting the PAN ENE triplex, we screened a library of immobilized small molecules and identified several triplex-binding chemotypes, the tightest of which exhibits micromolar binding affinity. Combined biophysical, biochemical, and computational strategies localized ligand binding to a platform created near a dinucleotide bulge at the base of the triplex. Crystal structures of apo (3.3 Å) and ligand-soaked (2.5 Å) ENE triplexes, which include a stabilizing basal duplex, indicate significant local structural rearrangements within this dinucleotide bulge. MD simulations and a modified nucleoside analog interference technique corroborate the role of the bulge and the base of the triplex in ligand binding. Together with recently discovered small molecules that reduce nuclear MALAT1 lncRNA levels by engaging its ENE triplex, our data supports the potential of targeting RNA triplexes with small molecules.

PARN-like Proteins Regulate Gene Expression in Land Plant Mitochondria by Modulating mRNA Polyadenylation

Mitochondria have their own double-stranded DNA genomes and systems to regulate transcription, mRNA processing, and translation. These systems differ from those operating in the host cell, and among eukaryotes. In recent decades, studies have revealed several plant-specific features of mitochondrial gene regulation. The polyadenylation status of mRNA is critical for its stability and translation in mitochondria. In this short review, I focus on recent advances in understanding the mechanisms regulating mRNA polyadenylation in plant mitochondria, including the role of poly(A)-specific ribonuclease-like proteins (PARNs). Accumulating evidence suggests that plant mitochondria have unique regulatory systems for mRNA poly(A) status and that PARNs play pivotal roles in these systems.

A comprehensive map of alternative polyadenylation in African American and European American lung cancer patients

Deciphering the post-transcriptional mechanisms (PTM) regulating gene expression is critical to understand the dynamics underlying transcriptomic regulation in cancer. Alternative polyadenylation (APA)-regulation of mRNA 3'UTR length by alternating poly(A) site usage-is a key PTM mechanism whose comprehensive analysis in cancer remains an important open challenge. Here we use a method and analysis pipeline that sequences 3'end-enriched RNA directly to overcome the saturation limitation of traditional 5'-3' based sequencing. We comprehensively map the APA landscape in lung cancer in a cohort of 98 tumor/non-involved tissues derived from European American and African American patients. We identify a global shortening of 3'UTR transcripts in lung cancer, with notable functional implications on the expression of both coding and noncoding genes. We find that APA of non-coding RNA transcripts (long non-coding RNAs and microRNAs) is a recurrent event in lung cancer and discover that the selection of alternative polyA sites is a form of non-coding RNA expression control. Our results indicate that mRNA transcripts from EAs are two times more likely than AAs to undergo APA in lung cancer. Taken together, our findings comprehensively map and identify the important functional role of alternative polyadenylation in determining transcriptomic heterogeneity in lung cancer.

Detecting and Profiling Endogenous RNA G-Quadruplexes in the Human Transcriptome

G-quadruplexes are the non-canonical nucleic acid structures that are preferentially formed in G-rich regions. This structure has been shown to be associated with many biological functions. Regardless of the broad efforts on DNA G-quadruplexes, we still have limited knowledge on RNA G-quadruplexes, especially in a transcriptome-wide manner. Herein, by integrating the DMS-seq and the bioinformatics pipeline, we profiled and depicted the RNA G-quadruplexes in the human transcriptome. The genes that contain RNA G-quadruplexes in their specific regions are significantly related to immune pathways and the COVID-19-related gene sets. Bioinformatics analysis reveals the potential regulatory functions of G-quadruplexes on miRNA targeting at the scale of the whole transcriptome. In addition, the G-quadruplexes are depleted in the putative, not the real, PAS-strong poly(A) sites, which may weaken the possibility of such sites being the real cleaved sites. In brief, our study provides insight into the potential function of RNA G-quadruplexes in post-transcription.

A computational pipeline to infer alternative poly-adenylation from 3' sequencing data

An increasing number of investigations have established alternative polyadenylation (APA) as a key mechanism of gene regulation through altering the length of 3' untranslated region (UTR) and generating distinct mRNA termini. Further, appreciation for the significance of APA in disease contexts propelled the development of several 3' sequencing techniques. While these RNA sequencing technologies have advanced APA analysis, the intrinsic limitation of 3' read coverage and lack of appropriate computational tools constrain precise mapping and quantification of polyadenylation sites. Notably, Poly(A)-ClickSeq (PAC-seq) overcomes limiting factors such as poly(A) enrichment and 3' linker ligation steps using click-chemistry. Here we provide an updated PolyA-miner protocol, a computational approach to analyze PAC-seq or other 3'-Seq datasets. As a key practical constraint, we also provide a detailed account on the impact of sequencing depth on the number of detected polyadenylation sites and APA changes. This protocol is also updated to handle unique molecular identifiers used to address PCR duplication potentially observed in PAC-seq.

Different roles for the adjoining and structurally similar A-rich and poly(A) domains of oskar mRNA: Only the A-rich domain is required for oskar noncoding RNA function, which includes MTOC positioning

Drosophila oskar (osk) mRNA has both coding and noncoding functions, with the latter required for progression through oogenesis. Noncoding activity is mediated by the osk 3' UTR. Three types of cis elements act most directly and are clustered within the final ~120 nucleotides of the 3' UTR: multiple binding sites for the Bru1 protein, a short highly conserved region, and A-rich sequences abutting the poly(A) tail. Here we extend the characterization of these elements and their functions, providing new insights into osk noncoding RNA function and the makeup of the cis elements. We show that all three elements are required for correct positioning of the microtubule organizing center (MTOC), a defect not previously reported for any osk mutant. Normally, the MTOC is located at the posterior of the oocyte during previtellogenic stages of oogenesis, and this distribution underlies the strong posterior enrichment of many mRNAs transported into the oocyte from the nurse cells. When osk noncoding function was disrupted the MTOC was dispersed in the oocyte and osk mRNA failed to be enriched at the posterior, although transport to the oocyte was not affected. A previous study did not detect loss of posterior enrichment for certain osk mutants lacking noncoding activity (Kanke et al., 2015). This discrepancy may be due to use of imaging aimed at monitoring transport to the oocyte rather than posterior enrichment. Involvement in MTOC positioning suggests that the osk noncoding function may act in conjunction with genes whose loss has similar effects, and that osk function may extend to other processes requiring those genes. Further characterization of the cis elements required for osk noncoding function included completion of saturation mutagenesis of the most highly conserved region, providing critical information for evaluating the possible contribution of candidate binding factors. The 3'-most cis element is a cluster of A-rich sequences, the ARS. The close juxtaposition and structural similarity of the ARS and poly(A) tail raised the possibility that they comprise an extended A-rich element required for osk noncoding function. We found that absence of the poly(A) tail did not mimic the effects of mutation of the ARS, causing neither arrest of oogenesis nor mispositioning of osk mRNA in previtellogenic stage oocytes. Thus, the ARS and the poly(A) tail are not interchangeable for osk noncoding RNA function, suggesting that the role of the ARS is not in recruitment of Poly(A) binding protein (PABP), the protein that binds the poly(A) tail. Furthermore, although PABP has been implicated in transport of osk mRNA from the nurse cells to the oocyte, mutation of the ARS in combination with loss of the poly(A) tail did not disrupt transport of osk mRNA into the oocyte. We conclude that PABP acts indirectly in osk mRNA transport, or is associated with osk mRNA independent of an A-rich binding site. Although the poly(A) tail was not required for osk mRNA transport into the oocyte, its absence was associated with a novel osk mRNA localization defect later in oogenesis, potentially revealing a previously unrecognized step in the localization process.

Transcript Isoform-Specific Estimation of Poly(A) Tail Length by Nanopore Sequencing of Native RNA

The poly(A) tail is a homopolymeric stretch of adenosine at the 3'-end of mature RNA transcripts and its length plays an important role in nuclear export, stability, and translational regulation of mRNA. Existing techniques for genome-wide estimation of poly(A) tail length are based on short-read sequencing. These methods are limited because they sequence a synthetic DNA copy of mRNA instead of the native transcripts. Furthermore, they can identify only a short segment of the transcript proximal to the poly(A) tail which makes it difficult to assign the measured poly(A) length uniquely to a single transcript isoform. With the introduction of native RNA sequencing by Oxford Nanopore Technologies, it is now possible to sequence full-length native RNA. A single long read contains both the transcript and the associated poly(A) tail, thereby making transcriptome-wide isoform-specific poly(A) tail length assessment feasible. We developed tailfindr-an R-based package for estimating poly(A) tail length from Oxford Nanopore sequencing data. In this chapter, we describe in detail the pipeline for transcript isoform-specific poly(A) tail profiling based on native RNA Nanopore sequencing-from library preparation to downstream data analysis with tailfindr.

Exploring Translational Control of Maternal mRNAs in Zebrafish

The study of translational regulation requires reliable measurement of both mRNA levels and protein synthesis. Cytoplasmic polyadenylation is a prevalent mode of translational regulation during oogenesis and early embryogenesis. Here the length of the poly(A) tail of an mRNA is coupled to its translatability. We describe a protocol to identify translationally regulated genes and measure their translation rate in the early zebrafish embryo using genome-wide polysome profiling. This protocol relies on the isolation of mRNA by means of an rRNA depletion strategy, which avoids capture bias due to short poly(A) tail that can occur when using conventional oligo(dT)-based methods. We also present a simple PCR-based method to measure the poly(A) tail length of selected mRNAs.

Phosphodiester modifications in mRNA poly(A) tail prevent deadenylation without compromising protein expression

Chemical modifications enable preparation of mRNAs with augmented stability and translational activity. In this study, we explored how chemical modifications of 5',3'-phosphodiester bonds in the mRNA body and poly(A) tail influence the biological properties of eukaryotic mRNA. To obtain modified and unmodified in vitro transcribed mRNAs, we used ATP and ATP analogs modified at the α-phosphate (containing either O-to-S or O-to-BH3 substitutions) and three different RNA polymerases-SP6, T7, and poly(A) polymerase. To verify the efficiency of incorporation of ATP analogs in the presence of ATP, we developed a liquid chromatography-tandem mass spectrometry (LC-MS/MS) method for quantitative assessment of modification frequency based on exhaustive degradation of the transcripts to 5'-mononucleotides. The method also estimated the average poly(A) tail lengths, thereby providing a versatile tool for establishing a structure-biological property relationship for mRNA. We found that mRNAs containing phosphorothioate groups within the poly(A) tail were substantially less susceptible to degradation by 3'-deadenylase than unmodified mRNA and were efficiently expressed in cultured cells, which makes them useful research tools and potential candidates for future development of mRNA-based therapeutics.

Paired rRNA-depleted and polyA-selected RNA sequencing data and supporting multi-omics data from human T cells

Both poly(A) enrichment and ribosomal RNA depletion are commonly used for RNA sequencing. Either has its advantages and disadvantages that may lead to biases in the downstream analyses. To better access these effects, we carried out both ribosomal RNA-depleted and poly(A)-selected RNA-seq for CD4+ T naive cells isolated from 40 healthy individuals from the Blueprint Project. For these 40 individuals, the genomic and epigenetic data were also available. This dataset offers a unique opportunity to understand how library construction influences differential gene expression, alternative splicing and molecular QTL (quantitative loci) analyses for human primary cells.

Tob2 phosphorylation regulates global mRNA turnover to reshape transcriptome and impact cell proliferation

Tob2, an anti-proliferative protein, promotes deadenylation through recruiting Caf1 deadenylase to the mRNA poly(A) tail by simultaneously interacting with both Caf1 and poly(A)-binding protein (PABP). Previously, we found that changes in Tob2 phosphorylation can alter its PABP-binding ability and deadenylation-promoting function. However, it remained unknown regarding the relevant kinase(s). Moreover, it was unclear whether Tob2 phosphorylation modulates the transcriptome and whether the phosphorylation is linked to Tob2's anti-proliferative function. In this study, we found that c-Jun amino-terminal kinase (JNK) increases phosphorylation of Tob2 at many Ser/Thr sites in the intrinsically disordered region (IDR) that contains two separate PABP-interacting PAM2 motifs. JNK-induced phosphorylation or phosphomimetic mutations at these sites weaken the Tob2-PABP interaction. In contrast, JNK-independent phosphorylation of Tob2 at serine 254 (S254) greatly enhances Tob2 interaction with PABP and its ability to promote deadenylation. We discovered that both PAM2 motifs are required for Tob2 to display these features. Combining mass spectrometry analysis, poly(A) size-distribution profiling, transcriptome-wide mRNA turnover analyses, and cell proliferation assays, we found that the phosphomimetic mutation at S254 (S254D) enhances Tob2's association with PABP, leading to accelerated deadenylation and decay of mRNAs globally. Moreover, the Tob2-S254D mutant accelerates the decay of many transcripts coding for cell cycle related proteins and enhances anti-proliferation function. Our findings reveal a novel mechanism by which Ccr4-Not complex is recruited by Tob2 to the mRNA 3' poly(A)-PABP complex in a phosphorylation dependent manner to promote rapid deadenylation and decay across the transcriptome, eliciting transcriptome reprogramming and suppressed cell proliferation.

Parallel poly(A) homo- and hetero-duplex formation detection with an adapted DNA nanoswitch technique

Polyriboadenylic [poly(rA)] strands of sufficient length form parallel double helices in acidic and/or ammonium-containing conditions. Poly(rA) duplexes in acidic conditions are held together by A⁺-A⁺ base-pairing also involving base interactions with the phosphate backbone. Traditional UV-melting studies of parallel poly(A) duplexes have typically examined homo-duplex formation of a single nucleic acid species in solution. We have adapted a technique utilizing a DNA nanoswitch that detects interaction of two different strands either with similar or differing lengths or modifications. Our method detected parallel duplex formation as a function of length, chemical modifications, and pH, and at a sensitivity that required over 100-fold less concentration of sample than prior UV-melting methods. While parallel polyriboadenylic acid and poly-2'-O-methyl-adenylic acid homo-duplexes formed, we did not detect homo-duplexes of polydeoxyriboadenylic acid strands or poly-locked nucleic acid (LNA)-adenylic strands. Importantly however, a poly-locked nucleic acid (LNA)-adenylic strand, as well as a poly-2'-O-methyl-adenylic strand, formed a hetero-duplex with a polyriboadenylic strand. Overall, our work validates a new tool for studying parallel duplexes and reveals fundamental properties of poly(A) parallel duplex formation. Parallel duplexes may find use in DNA nanotechnology and in molecular biology applications such as a potential poly(rA) tail capture tool as an alternative to traditional oligo(dT) based purification.

One of the Triple Poly(A) Signals in the M112-113 Gene Is Important and Sufficient for Stabilizing the M112-113 mRNA and the Replication of Murine Cytomegalovirus

The M112-113 gene is the first early gene of the murine cytomegalovirus (MCMV), and its expression is activated by the immediate-early 3 (IE3) protein during MCMV infection in permissive cells. At its 5' terminus, a 10-bp motif, upstream of the TATA box of the M112-113 gene, was identified to bind to IE3, and it is necessary for IE3 to activate M112-113 gene expression (Perez KJ et al. 2013 JVI). At the 3' terminus of the M112-113 gene, three poly(A) signals (PASs) are arranged closely, forming a PAS cluster. We asked whether it is necessary to have the PAS cluster for the M112-113 gene and wondered which PAS is required or important for M112-113 gene expression. In this study, we mutated one, two, or all three PASs in expressing plasmids. Then, we applied bacterial artificial chromosome (BAC) techniques to mutate PASs in viruses. Gene expression and viral replication were analyzed. We found that not all three PASs are needed for M112-113 gene expression. Moreover, we revealed that just one of the three poly(A)s is enough for MCMV replication. However, the deletion of all three PASs did not kill MCMV, although it significantly attenuated viral replication. Finally, an mRNA stability assay was performed and demonstrated that PASs are important to stabilize M112-113 mRNA. Therefore, we conclude that just one of the PASs of the M112-113 gene is sufficient and important for MCMV replication through the stabilization of M112-113 mRNA.

Requirement for cleavage factor IIm in the control of alternative polyadenylation in breast cancer cells

Alternative polyadenylation (APA) determines stability, localization and translation potential of the majority of mRNA in eukaryotic cells. The heterodimeric mammalian cleavage factor II (CF II_m) is required for pre-mRNA 3' end cleavage and is composed of the RNA kinase hClp1 and the termination factor hPcf11; the latter protein binds to RNA and the RNA polymerase II carboxy-terminal domain. Here, we used siRNA mediated knockdown and poly(A) targeted RNA sequencing to analyze the role of CF II_m in gene expression and APA in estrogen receptor positive MCF7 breast cancer cells. Identified gene ontology terms link CF II_m function to regulation of growth factor activity, protein heterodimerization and the cell cycle. An overlapping requirement for hClp1 and hPcf11 suggested that CF II_m protein complex was involved in the selection of proximal poly(A) sites. In addition to APA shifts within 3' untranslated regions (3'-UTRs), we observed shifts from promoter proximal regions to the 3'-UTR facilitating synthesis of full-length mRNAs. Moreover, we show that several truncated mRNAs that resulted from APA within introns in MCF7 cells cosedimented with ribosomal components in an EDTA sensitive manner suggesting that those are translated into protein. We propose that CF II_m contributes to the regulation of mRNA function in breast cancer.

The changing mouse embryo transcriptome at whole tissue and single-cell resolution

During mammalian embryogenesis, differential gene expression gradually builds the identity and complexity of each tissue and organ system1. Here we systematically quantified mouse polyA-RNA from day 10.5 of embryonic development to birth, sampling 17 tissues and organs. The resulting developmental transcriptome is globally structured by dynamic cytodifferentiation, body-axis and cell-proliferation gene sets that were further characterized by the transcription factor motif codes of their promoters. We decomposed the tissue-level transcriptome using single-cell RNA-seq (sequencing of RNA reverse transcribed into cDNA) and found that neurogenesis and haematopoiesis dominate at both the gene and cellular levels, jointly accounting for one-third of differential gene expression and more than 40% of identified cell types. By integrating promoter sequence motifs with companion ENCODE epigenomic profiles, we identified a prominent promoter de-repression mechanism in neuronal expression clusters that was attributable to known and novel repressors. Focusing on the developing limb, single-cell RNA data identified 25 candidate cell types that included progenitor and differentiating states with computationally inferred lineage relationships. We extracted cell-type transcription factor networks and complementary sets of candidate enhancer elements by using single-cell RNA-seq to decompose integrative cis-element (IDEAS) models that were derived from whole-tissue epigenome chromatin data. These ENCODE reference data, computed network components and IDEAS chromatin segmentations are companion resources to the matching epigenomic developmental matrix, and are available for researchers to further mine and integrate.

Detection and Degradation of Stalled Nascent Chains via Ribosome-Associated Quality Control

Stalled protein synthesis produces defective nascent chains that can harm cells. In response, cells degrade these nascent chains via a process called ribosome-associated quality control (RQC). Here, we review the irregularities in the translation process that cause ribosomes to stall as well as how cells use RQC to detect stalled ribosomes, ubiquitylate their tethered nascent chains, and deliver the ubiquitylated nascent chains to the proteasome. We additionally summarize how cells respond to RQC failure.

Poly(A)-seq: A method for direct sequencing and analysis of the transcriptomic poly(A)-tails

Poly(A) tails at the 3' end of eukaryotic messenger RNAs control mRNA stability and translation efficiency. Facilitated by various NGS methods, alternative polyadenylation sites determining the 3'-UTR length of gene transcripts have been extensively studied. However, poly(A) lengths demonstrating dynamic and developmental regulation remain largely unexplored. The recently developed NGS-based methods for genome-wide poly(A) profiling have promoted the study of genom-wide poly(A) dynamics. Here we present a straight forward NGS-method for poly(A) profiling, which applies a direct 3'-end adaptor ligation and the template switching for 5'-end adaptor ligation for cDNA library construction. Poly(A) lengths are directly calculated from base call data using a self-developed pipeline pA-finder. The libraries were directly sequenced from the 3'-UTR regions into the followed poly(A) tails, firstly on NextSeq 500 to produce single-end 300-nt reads, demonstrating the method feasibility and that optimization of the fragmented RNA size for cDNA library construction could detecting longer poly (A) tails. We next applied Poly(A)-seq cDNA libraries containing 40-nt and 120-nt poly(A) tail spike-in RNAs on HiSeq X-ten and NovaSeq 6000 to obtain 150-nt and 250-nt pair-end reads. The sequencing profiles of the spike-in RNAs demonstrated both high accuracy and high quality score in reading poly(A) tails. The poly(A) signal bleeding into the 3' adaptor sequence and a sharp decreased quality score at the junction were observed, allowing the modification of pA-finder to remove homopolymeric signal bleeding. We hope that wide applications of Poly(A)-seq help facilitate the study of the development- and disease-related poly(A) dynamics and regulation, and of the recent emerging mixed tailing regulation.

Establishment of 5'-3' interactions in mRNA independent of a continuous ribose-phosphate backbone

Functions of eukaryotic mRNAs are characterized by intramolecular interactions between their ends. We have addressed the question whether 5' and 3' ends meet by diffusion-controlled encounter "through solution" or by a mechanism involving the RNA backbone. For this purpose, we used a translation system derived from Drosophila embryos that displays two types of 5'-3' interactions: Cap-dependent translation initiation is stimulated by the poly(A) tail and inhibited by Smaug recognition elements (SREs) in the 3' UTR. Chimeric RNAs were made consisting of one RNA molecule carrying a luciferase coding sequence and a second molecule containing SREs and a poly(A) tail; the two were connected via a protein linker. The poly(A) tail stimulated translation of such chimeras even when disruption of the RNA backbone was combined with an inversion of the 5'-3' polarity between the open reading frame and poly(A) segment. Stimulation by the poly(A) tail also decreased with increasing RNA length. Both observations suggest that contacts between the poly(A) tail and the 5' end are established through solution, independently of the RNA backbone. In the same chimeric constructs, SRE-dependent inhibition of translation was also insensitive to disruption of the RNA backbone. Thus, tracking of the backbone is not involved in the repression of cap-dependent initiation. However, SRE-dependent repression was insensitive to mRNA length, suggesting that the contact between the SREs in the 3' UTR and the 5' end of the RNA might be established in a manner that differs from the contact between the poly(A) tail and the cap.

Critical role of deadenylation in regulating poly(A) rhythms and circadian gene expression

The mammalian circadian clock is deeply rooted in rhythmic regulation of gene expression. Rhythmic transcriptional control mediated by the circadian transcription factors is thought to be the main driver of mammalian circadian gene expression. However, mounting evidence has demonstrated the importance of rhythmic post-transcriptional controls, and it remains unclear how the transcriptional and post-transcriptional mechanisms collectively control rhythmic gene expression. In mouse liver, hundreds of genes were found to exhibit rhythmicity in poly(A) tail length, and the poly(A) rhythms are strongly correlated with the protein expression rhythms. To understand the role of rhythmic poly(A) regulation in circadian gene expression, we constructed a parsimonious model that depicts rhythmic control imposed upon basic mRNA expression and poly(A) regulation processes, including transcription, deadenylation, polyadenylation, and degradation. The model results reveal the rhythmicity in deadenylation as the strongest contributor to the rhythmicity in poly(A) tail length and the rhythmicity in the abundance of the mRNA subpopulation with long poly(A) tails (a rough proxy for mRNA translatability). In line with this finding, the model further shows that the experimentally observed distinct peak phases in the expression of deadenylases, regardless of other rhythmic controls, can robustly cluster the rhythmic mRNAs by their peak phases in poly(A) tail length and abundance of the long-tailed subpopulation. This provides a potential mechanism to synchronize the phases of target gene expression regulated by the same deadenylases. Our findings highlight the critical role of rhythmic deadenylation in regulating poly(A) rhythms and circadian gene expression.

Robust mapping of polyadenylated and non-polyadenylated RNA 3' ends at nucleotide resolution by 3'-end sequencing

3'-end poly(A)+ sequencing is an efficient and economical method for global measurement of mRNA levels and alternative poly(A) site usage. A common method involves oligo(dT)19V reverse-transcription (RT)-based library preparation and high-throughput sequencing with a custom primer ending in (dT)19. While the majority of library products have the first sequenced nucleotide reflect the bona fide poly(A) site (pA), a substantial fraction of sequencing reads arise from various mis-priming events. These can result in incorrect pA site calls anywhere from several nucleotides downstream to several kilobases upstream from the bona fide pA site. While these mis-priming events can be mitigated by increasing annealing stringency (e.g. increasing temperature from 37 °C to 42 °C), they still persist at an appreciable level (∼10%) and computational methods must be used to prevent artifactual calls. Here we present a bioinformatics workflow for precise mapping of poly(A)+ 3' ends and handling of artifacts due to oligo(dT) mis-priming and sample polymorphisms. We test pA site calling with three different read mapping programs (STAR, BWA, and BBMap), and show that the way in which each handles terminal mismatches and soft clipping has a substantial impact on identifying correct pA sites, with BWA requiring the least post-processing to correct artifacts. We demonstrate the use of this pipeline for mapping pA sites in the model eukaryote S. cerevisiae, and further apply this technology to non-polyadenylated transcripts by employing in vitro polyadenylation prior to library prep (IVP-seq). As proof of principle, we show that a fraction of tRNAs harbor CCU 3' tails instead of the canonical CCA tail, and globally identify 3' ends of splicing intermediates arising from inefficiently spliced transcripts.

Functional Cyclization of Eukaryotic mRNAs

The closed-loop model of eukaryotic translation states that mRNA is circularized by a chain of the cap-eIF4E-eIF4G-poly(A)-binding protein (PABP)-poly(A) interactions that brings 5' and 3' ends together. This circularization is thought to promote the engagement of terminating ribosomes to a new round of translation at the same mRNA molecule, thus enhancing protein synthesis. Despite the general acceptance and the elegance of the hypothesis, it has never been proved experimentally. Using continuous in situ monitoring of luciferase synthesis in a mammalian in vitro system, we show here that the rate of translation initiation at capped and polyadenylated reporter mRNAs increases after the time required for the first ribosomes to complete mRNA translation. Such acceleration strictly requires the presence of a poly(A)-tail and is abrogated by the addition of poly(A) RNA fragments or m⁷GpppG cap analog to the translation reaction. The optimal functional interaction of mRNA termini requires 5' untranslated region (UTR) and 3' UTR of moderate lengths and provides stronger acceleration, thus a longer poly(A)-tail. Besides, we revealed that the inhibitory effect of the dominant negative R362Q mutant of initiation factor eIF4A diminishes in the course of translation reaction, suggesting a relaxed requirement for ATP. Taken together, our results imply that, upon the functional looping of an mRNA, the recycled ribosomes can be recruited to the start codon of the same mRNA molecule in an eIF4A-independent fashion. This non-canonical closed-loop assisted reinitiation (CLAR) mode provides efficient translation of the functionally circularized mRNAs.

Nanopore direct RNA sequencing maps the complexity of Arabidopsis mRNA processing and m6A modification

Understanding genome organization and gene regulation requires insight into RNA transcription, processing and modification. We adapted nanopore direct RNA sequencing to examine RNA from a wild-type accession of the model plant Arabidopsis thaliana and a mutant defective in mRNA methylation (m6A). Here we show that m6A can be mapped in full-length mRNAs transcriptome-wide and reveal the combinatorial diversity of cap-associated transcription start sites, splicing events, poly(A) site choice and poly(A) tail length. Loss of m6A from 3' untranslated regions is associated with decreased relative transcript abundance and defective RNA 3' end formation. A functional consequence of disrupted m6A is a lengthening of the circadian period. We conclude that nanopore direct RNA sequencing can reveal the complexity of mRNA processing and modification in full-length single molecule reads. These findings can refine Arabidopsis genome annotation. Further, applying this approach to less well-studied species could transform our understanding of what their genomes encode.

PlantAPAdb: A Comprehensive Database for Alternative Polyadenylation Sites in Plants

Alternative cleavage and polyadenylation (APA) is increasingly recognized as an important regulatory mechanism in eukaryotic gene expression and is dynamically modulated in a developmental, tissue-specific, or environmentally responsive manner. Given the functional importance of APA and the rapid accumulation of APA sites in plants, a comprehensive and easily accessible APA site database is necessary for improved understanding of APA-mediated gene expression regulation. We present a database called PlantAPAdb that catalogs the most comprehensive APA site data derived from sequences from diverse 3' sequencing protocols and biological samples in plants. Currently, PlantAPAdb contains APA sites in six species, Oryza sativa (japonica and indica), Arabidopsis (Arabidopsis thaliana), Medicago truncatula, Trifolium pratense, Phyllostachys edulis, and Chlamydomonas reinhardtii APA sites in PlantAPAdb are available for bulk download and can be queried in a Google-like manner. PlantAPAdb provides rich information of the whole-genome APA sites, including genomic locations, heterogeneous cleavage sites, expression levels, and sample information. It also provides comprehensive poly(A) signals for APA sites in different genomic regions according to distinct profiles of cis-elements in plants. In addition, PlantAPAdb contains events of 3' untranslated region shortening/lengthening resulting from APA, which helps to understand the mechanisms underlying systematic changes in 3' untranslated region lengths. Additional information about conservation of APA sites in plants is also available, providing insights into the evolutionary polyadenylation configuration across species. As a user-friendly database, PlantAPAdb is a large and extendable resource for elucidating APA mechanisms, APA conservation, and gene expression regulation.

Structural and functional characterisation of human RNA helicase DHX8 provides insights into the mechanism of RNA-stimulated ADP release

DHX8 is a crucial DEAH-box RNA helicase involved in splicing and required for the release of mature mRNA from the spliceosome. Here, we report the biochemical characterisation of full-length human DHX8 and the catalytically active helicase core DHX8Δ547, alongside crystal structures of DHX8Δ547 bound to ADP and a structure of DHX8Δ547 bound to poly(A)6 single-strand RNA. Our results reveal that DHX8 has an in vitro binding preference for adenine-rich RNA and that RNA binding triggers the release of ADP through significant conformational flexibility in the conserved DEAH-, P-loop and hook-turn motifs. We demonstrate the importance of R620 and both the hook-turn and hook-loop regions for DHX8 helicase activity and propose that the hook-turn acts as a gatekeeper to regulate the directional movement of the 3' end of RNA through the RNA-binding channel. This study provides an in-depth understanding of the activity of DHX8 and contributes insights into the RNA-unwinding mechanisms of the DEAH-box helicase family.

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 comprehensive examination of Nanopore native RNA sequencing for characterization of complex transcriptomes

A platform for highly parallel direct sequencing of native RNA strands was recently described by Oxford Nanopore Technologies, but despite initial efforts it remains crucial to further investigate the technology for quantification of complex transcriptomes. Here we undertake native RNA sequencing of polyA + RNA from two human cell lines, analysing ~5.2 million aligned native RNA reads. To enable informative comparisons, we also perform relevant ONT direct cDNA- and Illumina-sequencing. We find that while native RNA sequencing does enable some of the anticipated advantages, key unexpected aspects currently hamper its performance, most notably the quite frequent inability to obtain full-length transcripts from single reads, as well as difficulties to unambiguously infer their true transcript of origin. While characterising issues that need to be addressed when investigating more complex transcriptomes, our study highlights that with some defined improvements, native RNA sequencing could be an important addition to the mammalian transcriptomics toolbox.