Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data

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

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

Dissecting the cis-regulatory syntax of transcription initiation with deep learning

Despite extensive characterization of mammalian Pol II transcription, the DNA sequence determinants of transcription initiation at a third of human promoters and most enhancers remain poorly understood. Hence, we trained and interpreted a neural network called ProCapNet that accurately models base-resolution initiation profiles from PRO-cap experiments using local DNA sequence. ProCapNet learns sequence motifs with distinct effects on initiation rates and TSS positioning and uncovers context-specific cryptic initiator elements intertwined within other TF motifs. ProCapNet annotates predictive motifs in nearly all actively transcribed regulatory elements across multiple cell-lines, revealing a shared cis-regulatory logic across promoters and enhancers mediated by a highly epistatic sequence syntax of cooperative and competitive motif interactions. ProCapNet models of RAMPAGE profiles measuring steady-state RNA abundance at TSSs distill initiation signals on par with models trained directly on PRO-cap profiles. ProCapNet learns a largely cell-type-agnostic cis-regulatory code of initiation complementing sequence drivers of cell-type-specific chromatin state critical for accurate prediction of cell-type-specific transcription initiation.

Characterizing nascent transcription patterns of PROMPTs, eRNAs, and readthrough transcripts in the ENCODE4 deeply profiled cell lines

Arising as co-products of canonical gene expression, transcription-associated lincRNAs, such as promoter upstream transcripts (PROMPTs), enhancer RNAs (eRNAs), and readthrough (RT) transcripts, are often regarded as byproducts of transcription, although they may be important for the expression of nearby genes. We identified regions of nascent expression of these lincRNA in 16 human cell lines using Bru-seq techniques, and found distinctly regulated patterns of PROMPT, eRNA, and RT transcription using the diverse biochemical approaches in the ENCODE4 deeply profiled cell lines collection. Transcription of these lincRNAs was influenced by sequence-specific features and the local or 3D chromatin landscape. However, these sequence and chromatin features do not describe the full spectrum of lincRNA expression variability we identify, highlighting the complexity of their regulation. This may suggest that transcription-associated lincRNAs are not merely byproducts, but rather that the transcript itself, or the act of its transcription, is important for genomic function.

A transcription factor (TF) inference method that broadly measures TF activity and identifies mechanistically distinct TF networks

TF profiler is a method of inferring transcription factor regulatory activity, i.e. when a TF is present and actively regulating transcription, directly directly from nascent sequencing assays such as PRO-seq and GRO-seq. Transcription factors orchestrate transcription and play a critical role in cellular maintenance, identity and response to external stimuli. While ChIP assays have measured DNA localization, they fall short of identifying when and where transcription factors are actively regulating transcription. Our method, on the other hand, uses RNA polymerase activity to infer TF activity across hundreds of data sets and transcription factors. Based on these classifications we identify three distinct classes of transcription factors: ubiquitous factors that play roles in cellular homeostasis, driving basal gene programs across tissues and cell types, tissue specific factors that act almost exclusively at enhancers and are themselves regulated at transcription, and stimulus responsive TFs which are regulated post-transcriptionally but act predominantly at enhancers. TF profiler is broadly applicable, providing regulatory insights on any PRO-seq sample for any transcription factor with a known binding motif.

Tdrd3-null mice show post-transcriptional and behavioral impairments associated with neurogenesis and synaptic plasticity

The Topoisomerase 3B (Top3b) - Tudor domain containing 3 (Tdrd3) protein complex is the only dual-activity topoisomerase complex that can alter both DNA and RNA topology in animals. TOP3B mutations in humans are associated with schizophrenia, autism and cognitive disorders; and Top3b-null mice exhibit several phenotypes observed in animal models of psychiatric and cognitive disorders, including impaired cognitive and emotional behaviors, aberrant neurogenesis and synaptic plasticity, and transcriptional defects. Similarly, human TDRD3 genomic variants have been associated with schizophrenia, verbal short-term memory and educational attainment. However, the importance of Tdrd3 in normal brain function has not been examined in animal models. Here we generated a Tdrd3-null mouse strain and demonstrate that these mice display both shared and unique defects when compared to Top3b-null mice. Shared defects were observed in cognitive behaviors, synaptic plasticity, adult neurogenesis, newborn neuron morphology, and neuronal activity-dependent transcription; whereas defects unique to Tdrd3-deficient mice include hyperactivity, changes in anxiety-like behaviors, olfaction, increased new neuron complexity, and reduced myelination. Interestingly, multiple genes critical for neurodevelopment and cognitive function exhibit reduced levels in mature but not nascent transcripts. We infer that the entire Top3b-Tdrd3 complex is essential for normal brain function, and that defective post-transcriptional regulation could contribute to cognitive and psychiatric disorders.

Precision run-on sequencing (PRO-seq) for microbiome transcriptomics

Bacteria respond to environmental stimuli through precise regulation of transcription initiation and elongation. Bulk RNA sequencing primarily characterizes mature transcripts, so to identify actively transcribed loci we need to capture RNA polymerase (RNAP) complexed with nascent RNA. However, such capture methods have only previously been applied to culturable, genetically tractable organisms such as E. coli and B. subtilis. Here we apply precision run-on sequencing (PRO-seq) to profile nascent transcription in cultured E. coli and diverse uncultured bacteria. We demonstrate that PRO-seq can characterize the transcription of small, structured, or post-transcriptionally modified RNAs, which are often absent from bulk RNA-seq libraries. Applying PRO-seq to the human microbiome highlights taxon-specific RNAP pause motifs and pause-site distributions across non-coding RNA loci that reflect structure-coincident pausing. We also uncover concurrent transcription and cleavage of CRISPR guide RNAs and transfer RNAs. We demonstrate the utility of PRO-seq for exploring transcriptional dynamics in diverse microbial communities.

Biodistribution of RNA Vaccines and of Their Products: Evidence from Human and Animal Studies

Explosive developments in mRNA vaccine technology in the last decade have made it possible to achieve great success in clinical trials of mRNA vaccines to prevent infectious diseases and develop cancer treatments and mRNA-based gene therapy products. The approval of the mRNA-1273 and BNT162b2 mRNA vaccines against SARS-CoV-2 by the U.S. Food and Drug Administration has led to mass vaccination (with mRNA vaccines) of several hundred million people around the world, including children. Despite its effectiveness in the fight against COVID-19, rare adverse effects of the vaccination have been shown in some studies, including vascular microcirculation disorders and autoimmune and allergic reactions. The biodistribution of mRNA vaccines remains one of the most poorly investigated topics. This mini-review discussed the results of recent experimental studies on humans and rodents regarding the biodistribution of mRNA vaccines, their constituents (mRNA and lipid nanoparticles), and their encoded antigens. We focused on the dynamics of the biodistribution of mRNA vaccine products and on the possibility of crossing the blood-brain and blood-placental barriers as well as transmission to infants through breast milk. In addition, we critically assessed the strengths and weaknesses of the detection methods that have been applied in these articles, whose results' reliability is becoming a subject of debate.

Subcytoplasmic location of translation controls protein output

The cytoplasm is highly compartmentalized, but the extent and consequences of subcytoplasmic mRNA localization in non-polarized cells are largely unknown. We determined mRNA enrichment in TIS granules (TGs) and the rough endoplasmic reticulum (ER) through particle sorting and isolated cytosolic mRNAs by digitonin extraction. When focusing on genes that encode non-membrane proteins, we observed that 52% have transcripts enriched in specific compartments. Compartment enrichment correlates with a combinatorial code based on mRNA length, exon length, and 3' UTR-bound RNA-binding proteins. Compartment-biased mRNAs differ in the functional classes of their encoded proteins: TG-enriched mRNAs encode low-abundance proteins with strong enrichment of transcription factors, whereas ER-enriched mRNAs encode large and highly expressed proteins. Compartment localization is an important determinant of mRNA and protein abundance, which is supported by reporter experiments showing that redirecting cytosolic mRNAs to the ER increases their protein expression. In summary, the cytoplasm is functionally compartmentalized by local translation environments.

Model-based characterization of the equilibrium dynamics of transcription initiation and promoter-proximal pausing in human cells

In metazoans, both transcription initiation and the escape of RNA polymerase (RNAP) from promoter-proximal pausing are key rate-limiting steps in gene expression. These processes play out at physically proximal sites on the DNA template and appear to influence one another through steric interactions. Here, we examine the dynamics of these processes using a combination of statistical modeling, simulation, and analysis of real nascent RNA sequencing data. We develop a simple probabilistic model that jointly describes the kinetics of transcription initiation, pause-escape, and elongation, and the generation of nascent RNA sequencing read counts under steady-state conditions. We then extend this initial model to allow for variability across cells in promoter-proximal pause site locations and steric hindrance of transcription initiation from paused RNAPs. In an extensive series of simulations, we show that this model enables accurate estimation of initiation and pause-escape rates. Furthermore, we show by simulation and analysis of real data that pause-escape is often strongly rate-limiting and that steric hindrance can dramatically reduce initiation rates. Our modeling framework is applicable to a variety of inference problems, and our software for estimation and simulation is freely available.

Hyperbranched polyamidoamine-RGD peptide/si-circICA1 in the treatment of invasive thyroid cancer through targeting of the miR-486-3p/SERPINA1 axis

Aim: This study aimed to identify molecular markers associated with papillary thyroid cancer (PTC) and investigate the therapeutic potential of targeted nanoscale drugs. Materials & methods: We analyzed the effects of circICA1 and miR-486-3p on B-CPAP cells' proliferation, apoptosis, migration and invasion. The regulation of the miR-486-3p/SERPINA1 axis was explored using quantitative real-time reverse transcription PCR and western blot analyses for metastasis. In vivo, we evaluated the effects of hyperbranched polyamidoamine-RGD peptide/si-circICA1 on PTC growth and metastasis. Results: Enhanced miR-486-3p expression inhibits B-CPAP cells' proliferation and invasion. si-circICA1 delivered via hyperbranched polyamidoamine-RGD peptide nanoparticles shows potential for treating metastasis in PTC. Conclusion: This study identifies key molecular mechanisms underlying PTC invasiveness and suggests a promising therapeutic strategy for PTC using targeted nanoscale drugs.

Functional analysis of a random-sequence chromosome reveals a high level and the molecular nature of transcriptional noise in yeast cells

We measure transcriptional noise in yeast by analyzing chromatin structure and transcription of an 18-kb region of DNA whose sequence was randomly generated. Nucleosomes fully occupy random-sequence DNA, but nucleosome-depleted regions (NDRs) are much less frequent, and there are fewer well-positioned nucleosomes and shorter nucleosome arrays. Steady-state levels of random-sequence RNAs are comparable to yeast mRNAs, although transcription and decay rates are higher. Transcriptional initiation from random-sequence DNA occurs at numerous sites, indicating very low intrinsic specificity of the RNA Pol II machinery. In contrast, poly(A) profiles of random-sequence RNAs are roughly comparable to those of yeast mRNAs, suggesting limited evolutionary restraints on poly(A) site choice. Random-sequence RNAs show higher cell-to-cell variability than yeast mRNAs, suggesting that functional elements limit variability. These observations indicate that transcriptional noise occurs at high levels in yeast, and they provide insight into how chromatin and transcription patterns arise from the evolved yeast genome.

Low-nutrient diet in Drosophila larvae stage causes enhancement in dopamine modulation in adult brain due epigenetic imprinting

Nutrient scarcity is a frequent adverse condition that organisms face during their development. This condition may lead to long-lasting effects on the metabolism and behaviour of adults due to developmental epigenetic modifications. Here, we show that reducing nutrient availability during larval development affects adult spontaneous activity and sleep behaviour, together with changes in gene expression and epigenetic marks in the mushroom bodies (MBs). We found that open chromatin regions map to 100 of 241 transcriptionally upregulated genes in the adult MBs, these new opening zones are preferentially located in regulatory zones such as promoter-TSS and introns. Importantly, opened chromatin at the Dopamine 1-like receptor 2 regulatory zones correlate with increased expression. In consequence, adult administration of a dopamine antagonist reverses increased spontaneous activity and diminished sleep time observed in response to early-life nutrient restriction. In comparison, reducing dop1R2 expression in MBs also ameliorates these effects, albeit to a lesser degree. These results lead to the conclusion that increased dopamine signalling in the MBs of flies reared in a poor nutritional environment underlies the behavioural changes observed due to this condition during development.

U1 snRNP increases RNA Pol II elongation rate to enable synthesis of long genes

The expansion of introns within mammalian genomes poses a challenge for the production of full-length messenger RNAs (mRNAs), with increasing evidence that these long AT-rich sequences present obstacles to transcription. Here, we investigate RNA polymerase II (RNAPII) elongation at high resolution in mammalian cells and demonstrate that RNAPII transcribes faster across introns. Moreover, we find that this acceleration requires the association of U1 snRNP (U1) with the elongation complex at 5' splice sites. The role of U1 to stimulate elongation rate through introns reduces the frequency of both premature termination and transcriptional arrest, thereby dramatically increasing RNA production. We further show that changes in RNAPII elongation rate due to AT content and U1 binding explain previous reports of pausing or termination at splice junctions and the edge of CpG islands. We propose that U1-mediated acceleration of elongation has evolved to mitigate the risks that long AT-rich introns pose to transcript completion.

Spurious intragenic transcription is a feature of mammalian cellular senescence and tissue aging

Mammalian aging is characterized by the progressive loss of tissue function and increased risk for disease. Accumulation of senescent cells in aging tissues partly contributes to this decline, and targeted depletion of senescent cells in vivo ameliorates many age-related phenotypes. The fundamental molecular mechanisms responsible for the decline of cellular health and fitness during senescence and aging are largely unknown. In this study, we investigated whether chromatin-mediated loss of transcriptional fidelity, known to contribute to fitness and survival in yeast and worms, also occurs during human cellular senescence and mouse aging. Our findings reveal aberrant transcription initiation inside genes during senescence and aging that co-occurs with changes in the chromatin landscape. Interventions that alter these spurious transcripts have profound consequences on cellular health, primarily affecting intracellular signal transduction pathways. We propose that age-related spurious transcription promotes a noisy transcriptome and degradation of coherent transcriptional networks.

Tdrd3-null mice show post-transcriptional and behavioral impairments associated with neurogenesis and synaptic plasticity

The Topoisomerase 3B (Top3b) - Tudor domain containing 3 (Tdrd3) protein complex is the only dual-activity topoisomerase complex in animals that can alter the topology of both DNA and RNA. TOP3B mutations in humans are associated with schizophrenia, autism and cognitive disorders; and Top3b-null mice exhibit several phenotypes observed in animal models of psychiatric and cognitive disorders, including impairments in cognitive and emotional behaviors, aberrant neurogenesis and synaptic plasticity, and transcriptional defects. Similarly, human TDRD3 genomic variants have been associated with schizophrenia, verbal shorten-memory and learning, and educational attainment. However, the importance of Tdrd3 in normal brain function has not been examined in animal models. Here we built a Tdrd3-null mouse strain and demonstrate that these mice display both shared and unique defects when compared to Top3b-null mice. Shared defects were observed in cognitive behaviors, synaptic plasticity, adult neurogenesis, newborn neuron morphology, and neuronal activity-dependent transcription; whereas defects unique to Tdrd3-deficient mice include hyperactivity, changes in anxiety-like behaviors, increased new neuron complexity, and reduced myelination. Interestingly, multiple genes critical for neurodevelopment and cognitive function exhibit reduced levels in mature but not nascent transcripts. We infer that the entire Top3b-Tdrd3 complex is essential for normal brain function, and that defective post-transcriptional regulation could contribute to cognitive impairment and psychiatric disorders.

The genetic and biochemical determinants of mRNA degradation rates in mammals

BACKGROUND

Degradation rate is a fundamental aspect of mRNA metabolism, and the factors governing it remain poorly characterized. Understanding the genetic and biochemical determinants of mRNA half-life would enable more precise identification of variants that perturb gene expression through post-transcriptional gene regulatory mechanisms.

RESULTS

We establish a compendium of 39 human and 27 mouse transcriptome-wide mRNA decay rate datasets. A meta-analysis of these data identified a prevalence of technical noise and measurement bias, induced partially by the underlying experimental strategy. Correcting for these biases allowed us to derive more precise, consensus measurements of half-life which exhibit enhanced consistency between species. We trained substantially improved statistical models based upon genetic and biochemical features to better predict half-life and characterize the factors molding it. Our state-of-the-art model, Saluki, is a hybrid convolutional and recurrent deep neural network which relies only upon an mRNA sequence annotated with coding frame and splice sites to predict half-life (r=0.77). The key novel principle learned by Saluki is that the spatial positioning of splice sites, codons, and RNA-binding motifs within an mRNA is strongly associated with mRNA half-life. Saluki predicts the impact of RNA sequences and genetic mutations therein on mRNA stability, in agreement with functional measurements derived from massively parallel reporter assays.

CONCLUSIONS

Our work produces a more robust ground truth for transcriptome-wide mRNA half-lives in mammalian cells. Using these revised measurements, we trained Saluki, a model that is over 50% more accurate in predicting half-life from sequence than existing models. Saluki succinctly captures many of the known determinants of mRNA half-life and can be rapidly deployed to predict the functional consequences of arbitrary mutations in the transcriptome.

Microfluidic live tracking and transcriptomics of cancer-immune cell doublets link intercellular proximity and gene regulation

Cell–cell communication and physical interactions play a vital role in cancer initiation, homeostasis, progression, and immune response. Here, we report a system that combines live capture of different cell types, co-incubation, time-lapse imaging, and gene expression profiling of doublets using a microfluidic integrated fluidic circuit that enables measurement of physical distances between cells and the associated transcriptional profiles due to cell–cell interactions. We track the temporal variations in natural killer—triple-negative breast cancer cell distances and compare them with terminal cellular transcriptome profiles. The results show the time-bound activities of regulatory modules and allude to the existence of transcriptional memory. Our experimental and bioinformatic approaches serve as a proof of concept for interrogating live-cell interactions at doublet resolution. Together, our findings highlight the use of our approach across different cancers and cell types.

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.

Genetic dissection of the RNA polymerase II transcription cycle

How DNA sequence affects the dynamics and position of RNA Polymerase II (Pol II) during transcription remains poorly understood. Here, we used naturally occurring genetic variation in F1 hybrid mice to explore how DNA sequence differences affect the genome-wide distribution of Pol II. We measured the position and orientation of Pol II in eight organs collected from heterozygous F1 hybrid mice using ChRO-seq. Our data revealed a strong genetic basis for the precise coordinates of transcription initiation and promoter proximal pause, allowing us to redefine molecular models of core transcriptional processes. Our results implicate DNA sequence, including both known and novel DNA sequence motifs, as key determinants of the position of Pol II initiation and pause. We report evidence that initiation site selection follows a stochastic process similar to Brownian motion along the DNA template. We found widespread differences in the position of transcription termination, which impact the primary structure and stability of mature mRNA. Finally, we report evidence that allelic changes in transcription often affect mRNA and ncRNA expression across broad genomic domains. Collectively, we reveal how DNA sequences shape core transcriptional processes at single nucleotide resolution in mammals.

A dual-activity topoisomerase complex regulates mRNA translation and turnover

Topoisomerase 3β (TOP3B) and TDRD3 form a dual-activity topoisomerase complex that interacts with FMRP and can change the topology of both DNA and RNA. Here, we investigated the post-transcriptional influence of TOP3B and associated proteins on mRNA translation and turnover. First, we discovered that in human HCT116 colon cancer cells, knock-out (KO) of TOP3B had similar effects on mRNA turnover and translation as did TDRD3-KO, while FMRP-KO resulted in rather distinct effects, indicating that TOP3B had stronger coordination with TDRD3 than FMRP in mRNA regulation. Second, we identified TOP3B-bound mRNAs in HCT116 cells; we found that while TOP3B did not directly influence the stability or translation of most TOP3B target mRNAs, it stabilized a subset of target mRNAs but had a more complex effect on translation-enhancing for some mRNAs whereas reducing for others. Interestingly, a point mutation that specifically disrupted TOP3B catalytic activity only partially recapitulated the effects of TOP3B-KO on mRNA stability and translation, suggesting that the impact of TOP3B on target mRNAs is partly linked to its ability to change topology of mRNAs. Collectively, our data suggest that TOP3B-TDRD3 can regulate mRNA translation and turnover by mechanisms that are dependent and independent of topoisomerase activity.

Candidate master microRNA regulator of arsenic-induced pancreatic beta cell impairment revealed by multi-omics analysis

Arsenic is a pervasive environmental toxin that is listed as the top priority for investigation by the Agency for Toxic Substance and Disease Registry. While chronic exposure to arsenic is associated with type 2 diabetes (T2D), the underlying mechanisms are largely unknown. We have recently demonstrated that arsenic treatment of INS-1 832/13 pancreatic beta cells impairs glucose-stimulated insulin secretion (GSIS), a T2D hallmark. We have also shown that arsenic alters the microRNA profile of beta cells. MicroRNAs have a well-established post-transcriptional regulatory role in both normal beta cell function and T2D pathogenesis. We hypothesized that there are microRNA master regulators that shape beta cell gene expression in pathways pertinent to GSIS after exposure to arsenicals. To test this hypothesis, we first treated INS-1 832/13 beta cells with either inorganic arsenic (iAsIII) or monomethylarsenite (MAsIII) and confirmed GSIS impairment. We then performed multi-omic analysis using chromatin run-on sequencing, RNA-sequencing, and small RNA-sequencing to define profiles of transcription, gene expression, and microRNAs, respectively. Integrating across these data sets, we first showed that genes downregulated by iAsIII treatment are enriched in insulin secretion and T2D pathways, whereas genes downregulated by MAsIII treatment are enriched in cell cycle and critical beta cell maintenance factors. We also defined the genes that are subject primarily to post-transcriptional control in response to arsenicals and demonstrated that miR-29a is the top candidate master regulator of these genes. Our results highlight the importance of microRNAs in arsenical-induced beta cell dysfunction and reveal both shared and unique mechanisms between iAsIII and MAsIII.

A compensatory link between cleavage/polyadenylation and mRNA turnover regulates steady-state mRNA levels in yeast

Cells have compensatory mechanisms to coordinate the rates of major biological processes, thereby permitting growth in a wide variety of conditions. Here, we uncover a compensatory link between cleavage/polyadenylation in the nucleus and messenger RNA (mRNA) turnover in the cytoplasm. On a global basis, same-gene 3' mRNA isoforms with twofold or greater differences in half-lives have steady-state mRNA levels that differ by significantly less than a factor of 2. In addition, increased efficiency of cleavage/polyadenylation at a specific site is associated with reduced stability of the corresponding 3' mRNA isoform. This inverse relationship between cleavage/polyadenylation and mRNA isoform half-life reduces the variability in the steady-state levels of mRNA isoforms, and it occurs in all four growth conditions tested. These observations suggest that during cleavage/polyadenylation in the nucleus, mRNA isoforms are marked in a manner that persists upon translocation to the cytoplasm and affects the activity of mRNA degradation machinery, thus influencing mRNA stability.

Deconvolution of Expression for Nascent RNA sequencing data (DENR) highlights pre-RNA isoform diversity in human cells

Motivation

Quantification of isoform abundance has been extensively studied at the mature RNA level using RNA-seq but not at the level of precursor RNAs using nascent RNA sequencing.

Results

We address this problem with a new computational method called Deconvolution of Expression for Nascent RNA-sequencing data (DENR), which models nascent RNA-sequencing read-counts as a mixture of user-provided isoforms. The baseline algorithm is enhanced by machine-learning predictions of active transcription start sites and an adjustment for the typical ‘shape profile’ of read-counts along a transcription unit. We show that DENR outperforms simple read-count-based methods for estimating gene and isoform abundances, and that transcription of multiple pre-RNA isoforms per gene is widespread, with frequent differences between cell types. In addition, we provide evidence that a majority of human isoform diversity derives from primary transcription rather than from post-transcriptional processes.

Availability and implementation

DENR and nascentRNASim are freely available at https://github.com/CshlSiepelLab/DENR (version v1.0.0) and https://github.com/CshlSiepelLab/nascentRNASim (version v0.3.0).

Supplementary information

Supplementary data are available at Bioinformatics online.

One-shot analysis of translated mammalian lncRNAs with AHARIBO

A vast portion of the mammalian genome is transcribed as long non-coding RNAs (lncRNAs) acting in the cytoplasm with largely unknown functions. Surprisingly, lncRNAs have been shown to interact with ribosomes, encode peptides, or act as ribosome sponges. These functions still remain mostly undetected and understudied owing to the lack of efficient tools for genome-wide simultaneous identification of ribosome-associated and peptide-producing lncRNAs. Here, we present AHA-mediated RIBOsome isolation (AHARIBO), a method for the detection of lncRNAs either untranslated, but associated with ribosomes, or encoding small peptides. Using AHARIBO in mouse embryonic stem cells during neuronal differentiation, we isolated ribosome-protected RNA fragments, translated RNAs, and corresponding de novo synthesized peptides. Besides identifying mRNAs under active translation and associated ribosomes, we found and distinguished lncRNAs acting as ribosome sponges or encoding micropeptides, laying the ground for a better functional understanding of hundreds of lncRNAs.