Abstract:Abstract Pseudouridine is one of the most abundant RNA modifications, occurring when uridines are catalyzed by Pseudouridine synthase proteins. It plays an important role in many biological processes and also has an importance in drug development. Recently, the single-molecule sequencing techniques such as the direct RNA sequencing platform offered by Oxford Nanopore technologies enable direct detection of RNA modifications on the molecule that is being sequenced, but to our knowledge this technology has not been used to identify RNA Pseudouridine sites. To this end, in this paper, we address this limitation by introducing a tool called Penguin that integrates several developed machine learning (ML) models (i.e., predictors) to identify RNA Pseudouridine sites in Nanopore direct RNA sequencing reads. Penguin extracts a set of features from the raw signal measured by the Oxford Nanopore and the corresponding basecalled k-mer. Those features are used to train the predictors included in Penguin, which in turn, is able to predict whether the signal is modified by the presence of Pseudouridine sites. We have included various predictors in Penguin including Support vector machine (SVM), Random Forest (RF), and Neural network (NN). The results on the two benchmark data sets show that Penguin is able to identify Pseudouridine sites with a high accuracy of 93.38% and 92.61% using SVM in random split testing and independent validation testing respectively. Thus, Penguin outperforms the existing Pseudouridine predictors in the literature that achieved an accuracy of 76.0 at most with an independent validation testing. A GitHub of the tool is accessible at https://github.com/Janga-Lab/Penguin . HIGHLIGHTS Penguin integrates several developed ML learning models (i.e., predictors) to identify RNA Ψ sites in Nanopore direct RNA sequencing reads. The pipeline of penguin automates the data preprocessing including Nanopore direct RNA reads alignment using Minimap2, and Nanopore signal extraction using Nanopolish, feature extraction from raw Nanopore signal for training ML predictors integrated in its platform, and the prediction of RNA Ψ sites with those predictors. Penguin can predict Ψ sites with a performance that outperforms the performance of the state-of-the-art research methods existing in the literature. Penguin platform can be adopted to be used for predicting other/various types of RNA modification. There are 6137606 U-mers samples predicted by penguin best ML model (SVM) as Ψ ones from a total of 67491289 U-mers samples in the complete Hek293 cell line with 556813 unique genomic location of Ψ. There are 1193192 U-mers samples predicted by penguin best ML model (SVM) as Ψ ones from a total of 229637931 U-mers samples in the complete Hela cell line with 39384 unique genomic locations of Ψ. There is a small fraction of 0.01% (6482 unique genomic locations) of Ψ that are common (overlapped) between both Hek293 and Hela cell lines. The extend of Ψ modification (the number of U-mers samples predicted as Ψ samples to the total number of U-mer samples in the complete RNA sequence of the cell line) in RNA sequence of Hek293 cell line is much greater than its counterpart for Hela cell line (9% for Hek293 versus 0. 5 % for Hela cell line).

Interferon inducible pseudouridine modification in human mRNA by quantitative nanopore profiling

Quantitative profiling of pseudouridylation dynamics in native RNAs with nanopore sequencing

Simultaneous nanopore profiling of mRNA m 6 A and pseudouridine reveals translation coordination

Integrative analysis of nanopore direct RNA sequencing data reveals a role of PUS7-dependent pseudouridylation in regulation of m6A and m5C modifications

Integrative analysis of nanopore direct RNA sequencing data reveals a role of PUS7-dependent pseudouridylation in regulation of m 6 A and m 5 C modifications

Nanopore signal deviations from pseudouridine modifications in RNA are sequence-specific: quantification requires dedicated synthetic controls

Quantitative Profiling of Pseudouridylation Landscape in the Human Transcriptome

Quantitative sequencing using BID-seq uncovers abundant pseudouridines in mammalian mRNA at base resolution

NanoMUD: Profiling of pseudouridine and N1-methylpseudouridine using Oxford Nanopore direct RNA sequencing

Absolute quantitative and base-resolution sequencing reveals comprehensive landscape of pseudouridine across the human transcriptome

mRNA psi profiling using nanopore DRS reveals cell-type-specific pseudouridylation

BID-seq for transcriptome-wide quantitative sequencing of mRNA pseudouridine at base resolution

Bisulfite and Nanopore Sequencing for Pseudouridine in RNA

Quantitative detection of pseudouridine in RNA by mass spectrometry

Penguin: A Tool for Predicting Pseudouridine Sites in Direct RNA Nanopore Sequencing Data

A comprehensive survey of RNA modifications in a human transcriptome

A bisulfite-assisted and ligation-based qPCR amplification technology for locus-specific pseudouridine detection at base resolution

Identification of nucleoside monophosphates and their epigenetic modifications using an engineered nanopore

Quantification of mRNA Expression Using Single-Molecule Nanopore Sensing

Quantitative analysis of tRNA abundance and modifications by nanopore RNA sequencing