Search

Eukaryotic mRNA deadenylation is generally considered as a two-step process in which the PAN2-PAN3 complex initiates the poly(A) tail degradation while, in the second step, the CCR4-NOT complex completes deadenylation, leading to decapping and degradation of the mRNA body. However, the mechanism of the biphasic poly(A) tail deadenylation remains enigmatic in several points such as the timing of the switch between the two steps, the role of translation termination and the mRNAs population involved. Here, we have studied the deadenylation of endogenous mRNAs in human cells depleted in either PAN3 or translation termination factor eRF3. Among the mRNAs tested, we found that only the endogenous ATF4 mRNA meets the biphasic model for deadenylation and that eRF3 prevents the shortening of its poly(A) tail. For the other mRNAs, the poor effect of PAN3 depletion on their poly(A) tail shortening questions the mode of their deadenylation. It is possible that these mRNAs experience a single step deadenylation process. Alternatively, we propose that a very short initial deadenylation by PAN2-PAN3 is followed by a rapid transition to the second phase involving CCR4-NOT complex. These differences in the timing of the transition from one deadenylation step to the other could explain the difficulties encountered in the generalization of the biphasic deadenylation model., Competing Interests: Declaration of competing interest I have no conflict of interest., (Copyright © 2021 The Authors. Published by Elsevier B.V. All rights reserved.)

Dual RNA-seq experiments examining viral and bacterial pathogens are increasing, but vary considerably in their experimental designs, such as infection rates and RNA depletion methods. Here, we have applied dual RNA-seq to Chlamydia trachomatis infected epithelial cells to examine transcriptomic responses from both organisms. We compared two time points post infection (1 and 24 h), three multiplicity of infection (MOI) ratios (0.1, 1 and 10) and two RNA depletion methods (rRNA and polyA). Capture of bacterial-specific RNA were greatest when combining rRNA and polyA depletion, and when using a higher MOI. However, under these conditions, host RNA capture was negatively impacted. Although it is tempting to use high infection rates, the implications on host cell survival, the potential reduced length of infection cycles and real world applicability should be considered. This data highlights the delicate nature of balancing host-pathogen RNA capture and will assist future transcriptomic-based studies to achieve more specific and relevant infection-related biological insights.

The RNA-binding protein Ataxin-2 regulates translation and mRNA stability through cytoplasmic polyadenylation of the targets. Here we newly identified DDX6 as a positive regulator of the cytoplasmic polyadenylation. Analysis of Ataxin-2 interactome using LC-MS/MS revealed prominent interaction with the DEAD-box RNA helicase DDX6. DDX6 interacted with components of the Ataxin-2 polyadenylation machinery; Ataxin-2, PABPC1 and PAPD4. As in the case for Ataxin-2 downregulation, DDX6 downregulation led to an increase in Ataxin-2 target mRNAs with short poly(A) tails as well as a reduction in their protein expression. In contrast, Ataxin-2 target mRNAs with short poly(A) tails were decreased by the overexpression of Ataxin-2, which was compromised by the DDX6 downregulation. However, polyadenylation induced by Ataxin-2 tethering was not affected by the DDX6 downregulation. Taken together, these results suggest that DDX6 positively regulates Ataxin-2-induced cytoplasmic polyadenylation to maintain poly(A) tail length of the Ataxin-2 targets provably through accelerating binding of Ataxin-2 to the target mRNAs., Competing Interests: Declaration of competing interest The authors declare that they have no conflict of interest., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Translation of problematic mRNA sequences induces ribosome stalling, triggering quality-control events, including ribosome rescue and nascent polypeptide degradation. To define the timing and regulation of these processes, we developed a SunTag-based reporter to monitor translation of a problematic sequence (poly[A]) in real time on single mRNAs. Although poly(A)-containing mRNAs undergo continuous translation over the timescale of minutes to hours, ribosome load is increased by ∼3-fold compared to a control, reflecting long queues of ribosomes extending far upstream of the stall. We monitor the resolution of these queues in real time and find that ribosome rescue is very slow compared to both elongation and termination. Modulation of pause strength, collision frequency, and the collision sensor ZNF598 reveals how the dynamics of ribosome collisions and their recognition facilitate selective targeting for quality control. Our results establish that slow clearance of stalled ribosomes allows cells to distinguish between transient and deleterious stalls., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Many long non-coding RNAs (lncRNA) are highly dysregulated in cancer and are emerging as therapeutic targets. One example is NEAT1, which consists of two overlapping lncRNA isoforms, NEAT1_1 (3.7 kb) and NEAT1_2 (23 kb), that are functionally distinct. The longer NEAT1_2 is responsible for scaffolding gene-regulatory nuclear bodies termed paraspeckles, whereas NEAT1_1 is involved in paraspeckle-independent function. The NEAT1 isoform ratio is dependent on the efficient cleavage and polyadenylation of NEAT1_1 at the expense of NEAT1_2. Here, we developed a targeted antisense oligonucleotide (ASO) approach to sterically block NEAT1_1 polyadenylation processing, achieving upregulation of NEAT1_2 and abundant paraspeckles. We have applied these ASOs to cells of the heterogeneous infant cancer, neuroblastoma, as we found higher NEAT1_1:NEAT1_2 ratio and lack of paraspeckles in high-risk neuroblastoma cells. These ASOs decrease NEAT1_1 levels, increase NEAT1_2/paraspeckles and concomitantly reduce cell viability in high-risk neuroblastoma specifically. In contrast, overexpression of NEAT1_1 has the opposite effect, increasing cell proliferation. Transcriptomic analyses of high-risk neuroblastoma cells with altered NEAT1 ratios and increased paraspeckle abundance after ASO treatment showed an upregulation of differentiation pathways, as opposed to the usual aggressive neuroblastic phenotype. Thus, we have developed potential anti-cancer ASO drugs that can transiently increase growth-inhibiting NEAT1_2 RNA at the expense of growth-promoting NEAT1_1 RNA. These ASOs, unlike others that degrade lncRNAs, provide insights into the importance of altering lncRNA polyadenylation events to suppress tumorigenesis as a strategy to combat cancer.

Herein the nucleic acid aptamers were attached to the polydeoxyadenylic acid (poly(dA)) tail for improving the tumor-targetability and cellular internalization of s-LNT/poly(dA) composite composed of two single chains of triple helical β-glucan lentinan (s-LNT) and one poly(dA) chain. The in vitro results demonstrate that the cellular uptake of s-LNT/poly(dA) composites in MCF-7 cancer cells was enhanced effectively after attaching the aptamer. The as-prepared fluorescin isothiocyanate (FITC)-labelled LNT (LNT-FITC) through grafting was used for tracing the enhanced tumor-targetability of the composites. As a result, the cellular internalization of the LNT-FITC into MCF-7 and 4T1 cancer cells was further increased by the aptamer conjugated to poly(dA). Meanwhile, the in vivo experiments further demonstrate more s-LNT/poly(dA)-aptamer composites were effectively accumulated at the tumor site compared with s-LNT alone. This work provides a novel strategy for fabricating triplex β-glucan as delivery vectors with active tumor-targetability., (Copyright © 2020. Published by Elsevier Ltd.)

The termination of transcription is a complex process that substantially contributes to gene regulation in eukaryotes. Previously, it was noted that a single cytosine deletion at the position + 32 bp relative to the single polyadenylation signal AAUAAA (hereafter the dC mutation) causes a 2-fold increase in the transcription level of the upstream eGFP reporter in mouse embryonic stem cells. Here, we analyzed the conservation of this phenomenon in immortalized mouse, human and drosophila cell lines and the influence of the dC mutation on the choice of the pre-mRNA cleavage sites. We have constructed dual-reporter plasmids to accurately measure the effect of the dC and other nearby located mutations on eGFP mRNA level by RT-qPCR. In this way, we found that the dC mutation leads to a 2-fold increase in the expression level of the upstream eGFP reporter gene in cultured mouse and human, but not in drosophila cells. In addition, 3' RACE analysis demonstrated that eGFP pre-mRNAs are cut at multiple positions between + 14 to + 31, and that the most proximal cleavage site becomes almost exclusively utilized in the presence of the dC mutation. We also identified new short sequence variations located within positions + 25.. + 40 and + 33.. + 48 that increase eGFP expression up to ~2-4-fold. Altogether, the positive effect of the dC mutation seems to be conserved in mouse embryonic stem cells, mouse embryonic 3T3 fibroblasts and human HEK293T cells. In the latter cells, the dC mutation appears to be involved in regulating pre-mRNA cleavage site selection. Finally, a multiplexed approach is proposed to identify motifs located downstream of cleavage site(s) that are essential for transcription termination.

Phosphorylated H2A.X is a critical chromatin marker of DNA damage repair (DDR) in higher eukaryotes. However, H2A.X gene expression remains relatively uncharacterised. Replication-dependent (RD) histone genes generate poly(A)- mRNA encoding new histones to package DNA during replication. In contrast, replication-independent (RI) histone genes synthesise poly(A)+ mRNA throughout the cell cycle, translated into histone variants that confer specific epigenetic patterns on chromatin. Remarkably H2AFX, encoding H2A.X, is a hybrid histone gene, generating both poly(A)+ and poly(A)- mRNA isoforms. Here we report that the selective removal of either mRNA isoform reveals different effects in different cell types. In some cells, RD H2A.X poly(A)- mRNA generates sufficient histone for deposition onto DDR associated chromatin. In contrast, cells making predominantly poly(A)+ mRNA require this isoform for de novo H2A.X synthesis, required for efficient DDR. This highlights the importance of differential H2A.X mRNA 3'-end processing in the maintenance of effective DDR.

RNA sequencing (RNA-seq) is extensively used to quantify gene expression transcriptome-wide. Although often paired with polyadenylate (poly(A)) selection to enrich for messenger RNA (mRNA), many applications require alternate approaches to counteract the high proportion of ribosomal RNA (rRNA) in total RNA. Recently, digestion using RNaseH and antisense DNA oligomers tiling target rRNAs has emerged as an alternative to commercial rRNA depletion kits. Here, we present a streamlined, more economical RNaseH-mediated rRNA depletion with substantially lower up-front costs, using shorter antisense oligos only sparsely tiled along the target RNA in a 5-min digestion reaction. We introduce a novel Web tool, Oligo-ASST, that simplifies oligo design to target regions with optimal thermodynamic properties, and additionally can generate compact, common oligo pools that simultaneously target divergent RNAs, e.g. across different species. We demonstrate the efficacy of these strategies by generating rRNA-depletion oligos for Xenopus laevis and for zebrafish, which expresses two distinct versions of rRNAs during embryogenesis. The resulting RNA-seq libraries reduce rRNA to <5% of aligned reads, on par with poly(A) selection, and also reveal expression of many non-adenylated RNA species. Oligo-ASST is freely available at https://mtleelab.pitt.edu/oligo to design antisense oligos for any taxon or to target any abundant RNA for depletion., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (Copyright © 2021 Elsevier Inc. All rights reserved.)

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.

Alternative polyadenylation (APA) is an essential regulatory mechanism for gene expression. The next generation sequencing provides ample opportunity to precisely delineate APA sites genome-wide. Various methods for profiling transcriptome-wide poly(A) sites were developed. By comparing available methods, the ways for adding sequencing adaptors to fit with the Illumina sequencing platform are different. These methods have identified more than 50% genes that undergo APA in eukaryotes. However, due to the unbalanced PCR during library preparation, accurate quantification of poly(A) sites is still a challenge. Here, we describe an updated poly(A) tag sequencing method that incorporates unique molecular identifier (UMI) into the adaptor for removing quantification bias induced by PCR duplicates. Hence, quantification of poly(A) site usages can be achieved by counting UMIs. This protocol, quantifying poly(A) tag sequencing (QPAT-seq), can be finished in 1 day with reduced cost, and is particularly useful for application with a large number of samples., (Copyright © 2021 Elsevier Inc. All rights reserved.)

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.

Polyadenylation of the 3' end of mRNAs is an important mechanism for regulating their stability and translation. We developed a nucleotide-resolution, transcriptome-wide, single-molecule SM-PAT-Seq method to accurately measure the polyA tail lengths of individual transcripts using long-read sequencing. The method generates cDNA using a double stranded splint adaptor targeting the far 3' end of the polyA tail for first strand synthesis along with random hexamers for second strand synthesis. This straight-forward method yields accurate polyA tail sequence lengths, can identify non-A residues in those tails, and quantitate transcript abundance., (2021 Published by Elsevier Inc.)

A repertoire of T cells with diverse antigen receptors is selected in the thymus. However, detailed mechanisms underlying this thymic positive selection are not clear. Here we show that the CCR4-NOT complex limits expression of specific genes through deadenylation of mRNA poly(A) tails, enabling positive selection. Specifically, the CCR4-NOT complex is up-regulated in thymocytes before initiation of positive selection, where in turn, it inhibits up-regulation of pro-apoptotic Bbc3 and Dab2ip. Elimination of the CCR4-NOT complex permits up-regulation of Bbc3 during a later stage of positive selection, inducing thymocyte apoptosis. In addition, CCR4-NOT elimination up-regulates Dab2ip at an early stage of positive selection. Thus, CCR4-NOT might control thymocyte survival during two-distinct stages of positive selection by suppressing expression levels of pro-apoptotic molecules. Taken together, we propose a link between CCR4-NOT-mediated mRNA decay and T cell selection in the thymus.

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-BH 3 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., (© 2020 Strzelecka et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Unstable repeat expansions and insertions cause more than 30 neurodegenerative and neuromuscular diseases. Remarkably, bidirectional transcription of repeat expansions has been identified in at least 14 of these diseases. More remarkably, a growing number of studies has been showing that both sense and antisense repeat RNAs are able to dysregulate important cellular pathways, contributing together to the observed clinical phenotype. Notably, antisense repeat RNAs from spinocerebellar ataxia type 7, myotonic dystrophy type 1, Huntington's disease and frontotemporal dementia/amyotrophic lateral sclerosis associated genes have been implicated in transcriptional regulation of sense gene expression, acting either at a transcriptional or posttranscriptional level. The recent evidence that antisense repeat RNAs could modulate gene expression broadens our understanding of the pathogenic pathways and adds more complexity to the development of therapeutic strategies for these disorders. In this review, we cover the amazing progress made in the understanding of the pathogenic mechanisms associated with repeat expansion neurodegenerative and neuromuscular diseases with a focus on the impact of antisense repeat transcription in the development of efficient therapies.

Embryonic stem (ES) cells are unique in their ability to self-renew indefinitely while maintaining pluripotency. Krüppel-like factor (Klf) 4 is an important member of the Klf family that is known to play a key role in pluripotency and somatic cell reprogramming. However, the identification and functional comparison of Klf4 splicing isoforms in mouse ESCs (mESCs) remains to be elucidated. Here, we identified three novel alternative splicing variants of Klf4 in mESCs-mKlf4-108, mKlf4-375 and mKlf4-1482-that are distinct from the previously known mKlf4-1449. mKlf4-1449 and mKlf4-1482 may stimulate the growth of ESCs, while mKlf4-108 can only promote the growth of ESCs in LIF low /serum conditions. In addition, both mKlf4-1449 and mKlf4-1482 can inhibit the differentiation of mESCs. However, the ability of mKlf4-1482 to promote self-renewal and inhibit differentiation is not as strong as that of mKlf4-1449. In contrast, both mKlf4-108 and mKlf4-375 may have the ability to induce endodermal differentiation. Taken together, we have identified for the first time the existence of alternative splicing variants of mKlf4 and have revealed their different roles, which provide new insights into the contribution of Klf4 to the self-renewal and pluripotency of mouse ESCs., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier Inc. All rights reserved.)

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.

Template-independent terminal ribonucleotide transferases (TENTs) catalyze the addition of nucleotide monophosphates to the 3'-end of RNA molecules regulating their fate. TENTs include poly(U) polymerases (PUPs) with a subgroup of 3' CUCU-tagging enzymes, such as CutA in Aspergillus nidulans. CutA preferentially incorporates cytosines, processively polymerizes only adenosines and does not incorporate or extend guanosines. The basis of this peculiar specificity remains to be established. Here, we describe crystal structures of the catalytic core of CutA in complex with an incoming non-hydrolyzable CTP analog and an RNA with three adenosines, along with biochemical characterization of the enzyme. The binding of GTP or a primer with terminal guanosine is predicted to induce clashes between 2-NH2 of the guanine and protein, which would explain why CutA is unable to use these ligands as substrates. Processive adenosine polymerization likely results from the preferential binding of a primer ending with at least two adenosines. Intriguingly, we found that the affinities of CutA for the CTP and UTP are very similar and the structures did not reveal any apparent elements for specific NTP binding. Thus, the properties of CutA likely result from an interplay between several factors, which may include a conformational dynamic process of NTP recognition., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The stability and processing of cellular RNA transcripts are efficiently controlled via non-templated addition of single or multiple nucleotides, which is catalyzed by various nucleotidyltransferases including poly(A) polymerases (PAPs). Germline development defective 2 (GLD-2) is among the first reported cytoplasmic non-canonical PAPs that promotes the translation of germline-specific mRNAs by extending their short poly(A) tails in metazoan, such as Caenorhabditis elegans and Xenopus. On the other hand, the function of mammalian GLD-2 seems more diverse, which includes monoadenylation of certain microRNAs. To understand the structural basis that underlies the difference between mammalian and non-mammalian GLD-2 proteins, we determine crystal structures of two rodent GLD-2s. Different from C. elegans GLD-2, mammalian GLD-2 is an intrinsically robust PAP with an extensively positively charged surface. Rodent and C. elegans GLD-2s have a topological difference in the β-sheet region of the central domain. Whereas C. elegans GLD-2 prefers adenosine-rich RNA substrates, mammalian GLD-2 can work on RNA oligos with various sequences. Coincident with its activity on microRNAs, mammalian GLD-2 structurally resembles the mRNA and miRNA processor terminal uridylyltransferase 7 (TUT7). Our study reveals how GLD-2 structurally evolves to a more versatile nucleotidyltransferase, and provides important clues in understanding its biological function in mammals., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (© 2020 Chen et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

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., (© 2020 Pickard et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

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.

Yeast cells undergoing the diauxic response show a striking upstream shift in poly(A) site utilization, with increased use of ORF-proximal poly(A) sites resulting in shorter 3' mRNA isoforms for most genes. This altered poly(A) pattern is extremely similar to that observed in cells containing Pol II derivatives with slow elongation rates. Conversely, cells containing derivatives with fast elongation rates show a subtle downstream shift in poly(A) sites. Polyadenylation patterns of many genes are sensitive to both fast and slow elongation rates, and a global shift of poly(A) utilization is strongly linked to increased purine content of sequences flanking poly(A) sites. Pol II processivity is impaired in diauxic cells, but strains with reduced processivity and normal Pol II elongation rates have normal polyadenylation profiles. Thus, Pol II elongation speed is important for poly(A) site selection and for regulating poly(A) patterns in response to environmental conditions., Competing Interests: JG, ZM No competing interests declared, KS Senior editor, eLife, (© 2020, Geisberg et al.)

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., (© 2020 Turner et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

During mammalian embryogenesis, differential gene expression gradually builds the identity and complexity of each tissue and organ system 1 . 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.

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) 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., Competing Interests: This study was supported by ABLife Inc.. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Gene expression is regulated by the rates of synthesis and degradation of mRNAs, but how these processes are coordinated is poorly understood. Here, we show that reduced transcription dynamics of specific genes leads to enhanced m 6 A deposition, preferential activity of the CCR4-Not complex, shortened poly(A) tails, and reduced stability of the respective mRNAs. These effects are also exerted by internal ribosome entry site (IRES) elements, which we found to be transcriptional pause sites. However, when transcription dynamics, and subsequently poly(A) tails, are globally altered, cells buffer mRNA levels by adjusting the expression of mRNA degradation machinery. Stress-provoked global impediment of transcription elongation leads to a dramatic inhibition of the mRNA degradation machinery and massive mRNA stabilization. Accordingly, globally enhanced transcription, such as following B cell activation or glucose stimulation, has the opposite effects. This study uncovers two molecular pathways that maintain balanced gene expression in mammalian cells by linking transcription to mRNA stability., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 Elsevier Inc. All rights reserved.)

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., (© 2020 Kluge et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

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., Competing Interests: The authors have declared that no competing interests exist.

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) 19 V 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., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

The biological significance of deadenylation in global gene expression is not fully understood. Here, we show that the CCR4-NOT deadenylase complex maintains expression of mRNAs, such as those encoding transcription factors, cell cycle regulators, DNA damage response-related proteins, and metabolic enzymes, at appropriate levels in the liver. Liver-specific disruption of Cnot1 , encoding a scaffold subunit of the CCR4-NOT complex, leads to increased levels of mRNAs for transcription factors, cell cycle regulators, and DNA damage response-related proteins because of reduced deadenylation and stabilization of these mRNAs. CNOT1 suppression also results in an increase of immature, unspliced mRNAs (pre-mRNAs) for apoptosis-related and inflammation-related genes and promotes RNA polymerase II loading on their promoter regions. In contrast, mRNAs encoding metabolic enzymes become less abundant, concomitant with decreased levels of these pre-mRNAs. Lethal hepatitis develops concomitantly with abnormal mRNA expression. Mechanistically, the CCR4-NOT complex targets and destabilizes mRNAs mainly through its association with Argonaute 2 (AGO2) and butyrate response factor 1 (BRF1) in the liver. Therefore, the CCR4-NOT complex contributes to liver homeostasis by modulating the liver transcriptome through mRNA deadenylation., (© 2020 Takahashi et al.)

Introns in different genes, or even different introns within the same gene, often have different splice sites and differ in splicing efficiency (SE). One expects mass-transcribed genes to have introns with higher SE than weakly transcribed genes. However, such a simple expectation cannot be tested directly because variable SE for these genes is often not measured. Mechanistically, SE should depend on signal strength at key splice sites (SS) such as 5'SS, 3'SS and branchpoint site (BPS), i.e., SE = F(5'SS, 3'SS, BPS). However, without SE, we again cannot model how these splice sites contribute to SE. Here I present an RNA-Seq approach to quantify SE for each of the 304 introns in yeast (Saccharomyces cerevisiae) genes, including 24 in the 5'UTR, by measuring 1) number of reads mapped to exon-exon junctions (N EE ) as a proxy for the abundance of spliced form, and 2) number of reads mapped to exon-intron junction (N EI5 and N EI3 at 5' and 3' ends of intron) as a proxy for the abundance of unspliced form. The total mRNA is N Total  = N EE  + p * N EI5  + (1-p) * N EI3 , with the simplest p = 0.5 but statistical methods were presented to estimate p from data. An estimated p is needed because N EI5 is expected to be smaller than N EI3 due to 1) step 1 splicing occurs before step 2 so EI5 is broken before EI3, 2) enrichment of poly(A) mRNA by oligo-dT, and 3) 5' degradation. SE is defined as the proportion (N EE /N Total ). Application of the method shows that ribosomal protein messages are efficiently and mostly cotranscriptionally spliced. Yeast genes with long introns are also spliced efficiently. HAC1/YFL031W is poorly spliced partly because its splicing involves a nonspliceosome mechanism and partly because Ire1p, which participate in splicing HAC1, is hardly expressed. Many putative yeast genes have low SE, and some splice sites are incorrectly annotated., (Copyright © 2019 Elsevier Inc. All rights reserved.)

To identify non-protein coding as well as truncated or premature RNA sequences expressed and obtain more complete transcriptome information, we combined the MinION direct RNA-sequencing of a conventional poly(A) RNA purification method with poly(A)-tagging of the non-coding RNA (ncRNA) fraction. This approach was applied to transcriptome sequencing of the dicyemid mesozoan, Dicyema misakiense, which has minicircular mitochondrial DNA molecules where each molecule encodes a single gene, as well as the host. Using informatics analysis, we distinguished dicyemid RNAs from those of the host squid. The poly(A) RNAs were assigned to host mitochondrial genes, host nuclear protein-coding genes, Dicyema nuclear protein-coding genes, and Dicyema mitochondrial genes in the decreasing order. Our poly(A)-tailing method recovered significantly more ncRNAs from the host compared with the sequencing of poly(A) RNAs. Furthermore, our method captured various lengths of squid mitochondrial DNA (mtDNA) transcripts at different steps of maturation including a read of 3,500 bp, which covers 21% of the squid mitochondrial genome, possibly a premature host RNA product. In contrast, shorter and less abundant reads were recovered from the dicyemid mitochondrial RNAs (mtRNAs). Even the longest read was 307 bp covering only a part of a minicircle. This study revealed significantly different modes of the mitochondrial transcription between a mesozoan and the host. Our approach to perform direct RNA-sequencing combined with the poly(A)-tailing reaction can be an effective method to fully capture non-poly(A) transcripts in a wide range of organisms., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Shortening of mRNA poly(A) tails (deadenylation) to trigger their decay is mediated mainly by the CCR4-NOT deadenylase complex. While four catalytic subunits (CNOT6, 6L 7, and 8) have been identified in the mammalian CCR4-NOT complex, their individual biological roles are not fully understood. In this study, we addressed the contribution of CNOT7/8 to viability of primary mouse embryonic fibroblasts (MEFs). We found that MEFs lacking CNOT7/8 expression [ Cnot7/8 -double knockout (dKO) MEFs] undergo cell death, whereas MEFs lacking CNOT6/6L expression ( Cnot6/6l -dKO MEFs) remain viable. Co-immunoprecipitation analyses showed that CNOT6/6L are also absent from the CCR4-NOT complex in Cnot7/8 -dKO MEFs. In contrast, either CNOT7 or CNOT8 still interacts with other subunits in the CCR4-NOT complex in Cnot6/6l -dKO MEFs. Exogenous expression of a CNOT7 mutant lacking catalytic activity in Cnot7/8 -dKO MEFs cannot recover cell viability, even though CNOT6/6L exists to some extent in the CCR4-NOT complex, confirming that CNOT7/8 is essential for viability. Bulk poly(A) tail analysis revealed that mRNAs with longer poly(A) tails are more numerous in Cnot7/8 -dKO MEFs than in Cnot6/6l -dKO MEFs. Consistent with elongated poly(A) tails, more mRNAs are upregulated and stabilized in Cnot7/8 -dKO MEFs than in Cnot6/6l -dKO MEFs. Importantly, Cnot6/6l -dKO mice are viable and grow normally to adulthood. Taken together, the CNOT7/8 catalytic subunits are essential for deadenylation, which is necessary to maintain cell viability, whereas CNOT6/6L are not.

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 7 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., Competing Interests: The authors declare no conflict of interest.

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 (m 6 A). Here we show that m 6 A 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 m 6 A from 3' untranslated regions is associated with decreased relative transcript abundance and defective RNA 3' end formation. A functional consequence of disrupted m 6 A 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., Competing Interests: MP, KK, AS, NS, KM, PG, AH, GB, GS No competing interests declared, (© 2020, Parker et al.)

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., (© 2020 The authors. All Rights Reserved.)

Sex hormones are best known for their influences on reproduction, but they also have profound influences on the immune response. Examples of sex-specific differences include: (i) the relatively poor control of influenza virus infections in males compared to females, (ii) allergic asthma, an IgE-associated hypersensitivity reaction that is exacerbated in adolescent females compared to males, and (iii) systemic lupus erythematosus, a life-threatening autoimmune disease with a 9:1 female:male bias. Here we consider how estrogen and estrogen receptor α (ERα) may influence the immune response by modifying class switch recombination (CSR) and immunoglobulin expression patterns. We focus on ERα binding to enhancers (Eμ and the 3' regulatory region) and switch sites (Sµ and Sε) in the immunoglobulin heavy chain locus. Our preliminary data from ChIP-seq analyses of purified, activated B cells show estrogen-mediated changes in the positioning of ERα binding within and near Sµ and Sε. In the presence of estrogen, ERα is bound not only to estrogen response elements (ERE), but also to adenosine-cytidine (AC)-repeats and poly adenosine (poly A) sequences, in some cases within constant region gene introns. We propose that by binding these sites, estrogen and ERα directly participate in the DNA loop formation required for CSR. We further suggest that estrogen regulates immunoglobulin expression patterns and can thereby influence life-and-death outcomes of infection, hypersensitivity, and autoimmune disease., (Copyright © 2019 Elsevier Inc. All rights reserved.)

High-throughput complementary DNA sequencing technologies have advanced our understanding of transcriptome complexity and regulation. However, these methods lose information contained in biological RNA because the copied reads are often short and modifications are not retained. We address these limitations using a native poly(A) RNA sequencing strategy developed by Oxford Nanopore Technologies. Our study generated 9.9 million aligned sequence reads for the human cell line GM12878, using thirty MinION flow cells at six institutions. These native RNA reads had a median length of 771 bases, and a maximum aligned length of over 21,000 bases. Mitochondrial poly(A) reads provided an internal measure of read-length quality. We combined these long nanopore reads with higher accuracy short-reads and annotated GM12878 promoter regions to identify 33,984 plausible RNA isoforms. We describe strategies for assessing 3' poly(A) tail length, base modifications and transcript haplotypes.

Small cytoplasmic mRNP granules are implicated in mRNA transport, translational control, and decay. Using super-resolution microscopy and fluorescence correlation spectroscopy, we analyzed the molecular composition and dynamics of single cytoplasmic YBX1&#95;IMP1 mRNP granules in live cells. Granules appeared elongated and branched, with patches of IMP1 and YBX1 distributed along mRNA, reflecting the attachment of the two RNA-binding proteins in cis. Particles form at the nuclear pore and do not associate with translating ribosomes, so the mRNP is a repository for mRNAs awaiting translation. In agreement with the average number of mRNA-binding sites derived from crosslinked immunoprecipitation (CLIP) analyses, individual mRNPs contain 5-15 molecules of YBX1 and IMP1 and a single poly(A) tail identified by PABPC1. Taken together, we conclude that small cytoplasmic mRNP granules are mRNA singletons, thus depicting the cellular transcriptome. Consequently, expression of functionally related mRNAs in RNA regulons is unlikely to result from coordinated assembly., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

This is a general protocol for the isolation of mRNA from total RNA using oligo(dT) coupled to magnetic beads. First, total RNA is dissolved in a high-salt buffer and heated briefly to 65°C-70°C, followed by immediate cooling on ice to disrupt secondary structures. The RNA is subsequently annealed to the oligo(dT)-magnetic beads at room temperature; the high-salt binding buffer stabilizes the poly(A)-oligo(dT) complexes. A high-salt washing buffer is then used to wash away unbound RNAs while retaining oligo(dT)-bound poly(A) + mRNAs. To elute the poly(A) + mRNAs from the beads, a low-salt buffer (or water) is used to destabilize the poly(A)-oligo(dT) complexes. Alternatively, poly(A) + mRNAs can be retained on the beads for downstream applications (e.g., solid-phase cDNA synthesis)., (© 2019 Cold Spring Harbor Laboratory Press.)

Single cell RNA sequencing methods have been increasingly used to understand cellular heterogeneity. Nevertheless, most of these methods suffer from one or more limitations, such as focusing only on polyadenylated RNA, sequencing of only the 3' end of the transcript, an exuberant fraction of reads mapping to ribosomal RNA, and the unstranded nature of the sequencing data. Here, we developed a novel single cell strand-specific total RNA library preparation method addressing all the aforementioned shortcomings. Our method was validated on a microfluidics system using three different cancer cell lines undergoing a chemical or genetic perturbation and on two other cancer cell lines sorted in microplates. We demonstrate that our total RNA-seq method detects an equal or higher number of genes compared to classic polyA[+] RNA-seq, including novel and non-polyadenylated genes. The obtained RNA expression patterns also recapitulate the expected biological signal. Inherent to total RNA-seq, our method is also able to detect circular RNAs. Taken together, SMARTer single cell total RNA sequencing is very well suited for any single cell sequencing experiment in which transcript level information is needed beyond polyadenylated genes., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (© 2019 The Author(s).)

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.

Alternative cleavage and polyadenylation (APA) can occur at more than half of all human genes, greatly enhancing the cellular repertoire of mRNA isoforms. As these isoforms can have altered stability, localisation and coding potential, deregulation of APA can disrupt gene expression and this has been linked to many diseases including cancer progression. How APA generates cancer-specific isoform profiles and what their physiological consequences are, however, is largely unclear. Here we use a subcellular fractionation approach to determine the nuclear and cytoplasmic APA profiles of successive stages of colon cancer using a cell line-based model. Using this approach, we show that during cancer progression specific APA profiles are established. We identify that overexpression of hnRNPC has a critical role in the establishment of APA profiles characteristic for metastatic colon cancer cells, by regulating poly(A) site selection in a subset of genes that have been implicated in cancer progression including MTHFD1L., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Polyadenylation signals (PAS) are found in most protein-coding and some non-coding genes in eukaryotes. Their accurate recognition improves understanding gene regulation mechanisms and recognition of the 3'-end of transcribed gene regions where premature or alternate transcription ends may lead to various diseases. Although different methods and tools for in-silico prediction of genomic signals have been proposed, the correct identification of PAS in genomic DNA remains challenging due to a vast number of non-relevant hexamers identical to PAS hexamers. In this study, we developed a novel method for PAS recognition. The method is implemented in a hybrid PAS recognition model (HybPAS), which is based on deep neural networks (DNNs) and logistic regression models (LRMs). One of such models is developed for each of the 12 most frequent human PAS hexamers. DNN models appeared the best for eight PAS types (including the two most frequent PAS hexamers), while LRM appeared best for the remaining four PAS types. The new models use different combinations of signal processing-based, statistical, and sequence-based features as input. The results obtained on human genomic data show that HybPAS outperforms the well-tuned state-of-the-art Omni-PolyA models, reducing the classification error for different PAS hexamers by up to 57.35% for 10 out of 12 PAS types, with Omni-PolyA models being better for two PAS types. For the most frequent PAS types, 'AATAAA' and 'ATTAAA', HybPAS reduced the error rate by 35.14% and 34.48%, respectively. On average, HybPAS reduces the error by 30.29%. HybPAS is implemented partly in Python and in MATLAB available at https://github.com/EMANG-KAUST/PolyA&#95;Prediction&#95;LRM&#95;DNN., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

Synthesis of the negative-strand ((-)-strand) counterpart is the first step of coronavirus (CoV) replication; however, the detailed mechanism of the early event and the factors involved remain to be determined. Here, using bovine coronavirus (BCoV)-defective interfering (DI) RNA, we showed that (a) a poly(A) tail with a length of 15 nucleotides (nt) was sufficient to initiate efficient (-)-strand RNA synthesis and (b) substitution of the poly(A) tail with poly(U), (C) or (G) only slightly decreased the efficiency of (-)-strand synthesis. The findings indicate that in addition to the poly(A) tail, other factors acting in trans may also participate in (-)-strand synthesis. The BCoV nucleocapsid (N) protein, an RNA-binding protein, was therefore tested as a candidate. Based on dissociation constant (K d ) values, it was found that the binding affinity between N protein, but not poly(A)-binding protein, and the 3'-terminal 55 nt plus a poly(A), poly(U), poly(C) or poly(G) tail correlates with the efficiency of (-)-strand synthesis. Such an association was also evidenced by the binding affinity between the N protein and 5'- and 3'-terminal cis-acting elements important for (-)-strand synthesis. Further analysis demonstrated that N protein can act as a bridge to facilitate interaction between the 5'- and 3'-ends of the CoV genome, leading to circularization of the genome. Together, the current study extends our understanding of the mechanism of CoV (-)-strand RNA synthesis through involvement of N protein and genome circularization and thus may explain why the addition of N protein in trans is required for efficient CoV replication., (© 2019 Federation of European Biochemical Societies.)