Enhanced detection of RNA modifications and read mapping with high-accuracy nanopore RNA basecalling models.

In recent years, nanopore direct RNA sequencing (DRS) became a valuable tool for studying the epitranscriptome, due to its ability to detect multiple modifications within the same full-length native RNA molecules. While RNA modifications can be identified in the form of systematic basecalling 'errors' in DRS datasets, N6-methyladenosine (m6A) modifications produce relatively low 'errors' compared to other RNA modifications, limiting the applicability of this approach to m6A sites that are modified at high stoichiometries. Here, we demonstrate that the use of alternative RNA basecalling models, trained with fully unmodified sequences, increases the 'error'signal of m6A, leading to enhanced detection and improved sensitivity even at low stoichiometries. Moreover, we find that high-accuracy alternative RNA basecalling models can show up to 97% median basecalling accuracy, outperforming currently available RNA basecalling models, which show 91% median basecalling accuracy. Notably, the use of high-accuracy basecalling models is accompanied by a significant increase in the number of mapped reads -especially in shorter RNA fractions- and increased basecalling error signatures at pseudouridine (Ψ) and N1-methylpseudouridine (m1Ψ) modified sites. Overall, our work demonstrates that alternative RNA basecalling models can be used to improve the detection of RNA modifications, read mappability, and basecalling accuracy in nanopore DRS datasets.

RNA Pol II-dependent transcription efficiency fine-tunes A-to-I editing levels.

A-to-I RNA editing is a widespread epitranscriptomic phenomenon leading to the conversion of adenosines to inosines, which are primarily interpreted as guanosines by cellular machines. Consequently, A-to-I editing can alter splicing or lead to recoding of transcripts. As misregulation of editing can cause a variety of human diseases, A-to-I editing requires tight regulation of the extent of deamination, particularly in protein-coding regions. The bulk of A-to-I editing occurs cotranscriptionally. Thus, we studied A-to-I editing regulation in the context of transcription and pre-mRNA processing. We show that stimulation of transcription impacts editing levels. Activation of the transcription factor MYC leads to an up-regulation of A-to-I editing, particularly in transcripts that are suppressed upon MYC activation. Moreover, low pre-mRNA synthesis rates and low pre-mRNA expression levels support high levels of editing. We also show that editing levels greatly differ between nascent pre-mRNA and mRNA in a cellular system, as well as in mouse tissues. Editing levels can increase or decrease from pre-mRNA to mRNA and can vary across editing targets and across tissues, showing that pre-mRNA processing is an important layer of editing regulation. Several lines of evidence suggest that the differences emerge during pre-mRNA splicing. Moreover, actinomycin D treatment of primary neuronal cells and editing level analysis suggests that regulation of editing levels also depends on transcription.

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

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