MasterOfPores: A Workflow for the Analysis of Oxford Nanopore Direct RNA Sequencing Datasets

Abstract

The direct RNA sequencing platform offered by Oxford Nanopore Technologies allows for direct measurement of RNA molecules without the need of conversion to complementary DNA, fragmentation or amplification. As such, it is virtually capable of detecting any given RNA modification present in the molecule that is being sequenced, as well as provide polyA tail length estimations at the level of individual RNA molecules. Although this technology has been publicly available since 2017, the complexity of the raw Nanopore data, together with the lack of systematic and reproducible pipelines, have greatly hindered the access of this technology to the general user. Here we address this problem by providing a fully benchmarked workflow for the analysis of direct RNA sequencing reads, termed MasterOfPores. The pipeline starts with a pre-processing module, which converts raw current intensities into multiple types of processed data including FASTQ and BAM, providing metrics of the quality of the run, quality-filtering, demultiplexing, base-calling and mapping. In a second step, the pipeline performs downstream analyses of the mapped reads, including prediction of RNA modifications and estimation of polyA tail lengths. Four direct RNA MinION sequencing runs can be fully processed and analyzed in 10 h on 100 CPUs. The pipeline can also be executed in GPU locally or in the cloud, decreasing the run time fourfold. The software is written using the NextFlow framework for parallelization and portability, and relies on Linux containers such as Docker and Singularity for achieving better reproducibility. The MasterOfPores workflow can be executed on any Unix-compatible OS on a computer, cluster or cloud without the need of installing any additional software or dependencies, and is freely available in Github (https://github.com/biocorecrg/master_of_pores). This workflow simplifies direct RNA sequencing data analyses, facilitating the study of the (epi)transcriptome at single molecule resolution.

Classification of SARS-CoV-2 and non-SARS-CoV-2 using machine learning algorithms

“…Besides, other human, mammals, and birds’ coronaviruses [ 31 ], were incorporated into the non-SARS-CoV-2 group to assess the robustness and effectiveness of the proposed algorithms. We also downloaded a control sample from the Epitranscriptomics and RNA Dynamics Lab (Novoa Lab) and the Bioinformatics Core Facility (BioCore) at the Centre for Genome Regulation [ 32 ]. Table 1 depicts information on coronavirus species, sample size, and designated label for both SARS-CoV-2 and non-SARS-CoV-2 groups.…”

Updated insight into COVID-19 disease and health management to combat the pandemic

“…Nanopore sequencing is a third-generation sequencing technique that has been used for detailed viral genome analysis from various clinical specimens. 123 , 124 This technology also measures RNA molecules from the viral genomic sample without the necessity of generating cDNA for analysis. Such techniques offer complete genome information of the pathogen, thus providing novel research insights that help in understanding the nature of the pathogen.…”

Analysis of RNA Modifications by Second- and Third-Generation Deep Sequencing: 2020 Update

“…Another intermediate solution, a workflow for the analysis of direct RNA sequencing reads, termed MasterOfPores, converts raw current intensities into multiple types of processed data, including FASTQ and BAM, providing metrics of the quality of the run, quality filtering, demultiplexing, base calling, and mapping. In a second step, the pipeline performs downstream analyses of the mapped reads, including the prediction of RNA modifications and the estimation of polyA tail lengths [ 160 ].…”