ABSTRACT
N6-methyladenosine (m6A) is the most abundant internal eukaryotic mRNA modification, and is involved in the regulation of various biological processes. Direct Nanopore sequencing of native RNA (dRNA-seq) emerged as a leading approach for its identification. Several software were published for m6A detection and there is a strong need for independent studies benchmarking their performance on data from different species, and against various reference datasets. Moreover, a computational workflow is needed to streamline the execution of tools whose installation and execution remains complicated. We developed NanOlympicsMod, a Nextflow pipeline exploiting containerized technology for comparing 14 tools for m6A detection on dRNA-seq data. NanOlympicsMod was tested on dRNA-seq data generated from in vitro (un)modified synthetic oligos. The m6A hits returned by each tool were compared to the m6A position known by design of the oligos. In addition, NanOlympicsMod was used on dRNA-seq datasets from wild-type and m6A-depleted yeast, mouse and human, and each tool's hits were compared to reference m6A sets generated by leading orthogonal methods. The performance of the tools markedly differed across datasets, and methods adopting different approaches showed different preferences in terms of precision and recall. Changing the stringency cut-offs allowed for tuning the precision-recall trade-off towards user preferences. Finally, we determined that precision and recall of tools are markedly influenced by sequencing depth, and that additional sequencing would likely reveal additional m6A sites. Thanks to the possibility of including novel tools, NanOlympicsMod will streamline the benchmarking of m6A detection tools on dRNA-seq data, improving future RNA modification characterization.
Subject(s)
Adenine/analogs & derivatives , Nanopore Sequencing , Nanopores , Humans , Animals , Mice , RNA/genetics , Benchmarking , Sequence Analysis, RNA/methodsABSTRACT
The SARS-CoV-2 virus has a complex transcriptome characterised by multiple, nested subgenomic RNAsused to express structural and accessory proteins. Long-read sequencing technologies such as nanopore direct RNA sequencing can recover full-length transcripts, greatly simplifying the assembly of structurally complex RNAs. However, these techniques do not detect the 5' cap, thus preventing reliable identification and quantification of full-length, coding transcript models. Here we used Nanopore ReCappable Sequencing (NRCeq), a new technique that can identify capped full-length RNAs, to assemble a complete annotation of SARS-CoV-2 sgRNAs and annotate the location of capping sites across the viral genome. We obtained robust estimates of sgRNA expression across cell lines and viral isolates and identified novel canonical and non-canonical sgRNAs, including one that uses a previously un-annotated leader-to-body junction site. The data generated in this work constitute a useful resource for the scientific community and provide important insights into the mechanisms that regulate the transcription of SARS-CoV-2 sgRNAs.
Subject(s)
COVID-19 , Nanopores , RNA, Guide, Kinetoplastida/chemistry , COVID-19/genetics , Genome, Viral/genetics , Humans , RNA Caps , RNA, Viral/genetics , RNA, Viral/metabolism , SARS-CoV-2/geneticsABSTRACT
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a positive single-stranded RNA virus, engages in complex interactions with host cell proteins throughout its life cycle. While these interactions enable the host to recognize and inhibit viral replication, they also facilitate essential viral processes such as transcription, translation, and replication. Many aspects of these virus-host interactions remain poorly understood. Here, we employed the catRAPID algorithm and utilized the RNA-protein interaction detection coupled with mass spectrometry technology to predict and validate the host proteins that specifically bind to the highly structured 5' and 3' terminal regions of the SARS-CoV-2 RNA. Among the interactions identified, we prioritized pseudouridine synthase PUS7, which binds to both ends of the viral RNA. Using nanopore direct RNA sequencing, we discovered that the viral RNA undergoes extensive post-transcriptional modifications. Modified consensus regions for PUS7 were identified at both terminal regions of the SARS-CoV-2 RNA, including one in the viral transcription regulatory sequence leader. Collectively, our findings offer insights into host protein interactions with the SARS-CoV-2 UTRs and highlight the likely significance of pseudouridine synthases and other post-transcriptional modifications in the viral life cycle. This new knowledge enhances our understanding of virus-host dynamics and could inform the development of targeted therapeutic strategies.