Nanopore direct RNA sequencing (dRNA-seq) enables unique insights into RNA biology. However, applications are currently limited by the lack of accurate and cost-effective sample multiplexing. Here we introduce WarpDemuX, an ultra-fast and highly accurate adapter-barcoding and demultiplexing approach for dRNA-seq with SQK-RNA002 and SQK-RNA004 chemistries. WarpDemuX enhances speed and accuracy by fast processing of the raw nanopore signal, use of a light-weight machine-learning algorithm and design of optimized barcode sets. We demonstrate its utility by performing rapid phenotypic profiling of different SARS-CoV-2 viruses through multiplexed sequencing of longitudinal samples on a single flowcell, identifying systematic differences in transcript abundance and poly(A) tail lengths during infection. Additionally, integrating WarpDemuX into sequencing control software enables real-time enrichment of target molecules through barcode-specific adaptive sampling, which we demonstrate by enriching low abundance viral RNA. In summary, WarpDemuX represents a broadly applicable, high-performance, economical multiplexing solution for dRNA-seq, facilitating advanced (epi-) transcriptomic research.
Human immunodeficiency virus-1 (HIV-1) uses a number of strategies to modulate viral and host gene expression during its life cycle. To characterize the transcriptional and translational landscape of HIV-1 infected cells, we used a combination of ribosome profiling, disome sequencing and RNA sequencing. We show that HIV-1 messenger RNAs are efficiently translated at all stages of infection, despite evidence for a substantial decrease in the translational efficiency of host genes that are implicated in host cell translation. Our data identify upstream open reading frames in the HIV-1 5′-untranslated region as well as internal open reading frames in the Vif and Pol coding domains. We also observed ribosomal collisions in Gag-Pol upstream of the ribosome frameshift site that we attributed to an RNA structural fold using RNA structural probing and functional analysis. Antisense oligonucleotides designed to alter the base of this structure decreased frameshift efficiency. Overall, our data highlight the complexity of HIV-1 gene regulation and provide a key resource for decoding of host–pathogen interactions upon HIV-1 infection. Furthermore, we provide evidence for a RNA structural fold including the frameshift site that could serve as a target for antiviral therapy. A study by Kibe et al. of the translational landscape of HIV-1 reveals alternative translation events, including use of upstream and internal open reading frames, an RNA fold regulating gag-pol frameshifting and ribosome collisions, providing new targets for antivirals.
Targeting of diseased cells is one of the most urgently needed prerequisites for a next generation of potent pharmaceuticals. Different approaches pursued fail mainly due to a lack of specific surface markers. Developing an RNA-based methodology, we can now ensure precise cell targeting combined with selective expression of effector proteins for therapy, diagnostics or cell steering. The specific combination of the molecular properties of antisense technology and mRNA therapy with functional RNA secondary structures allowed us to develop selectively expressed RNA molecules for medical applications. These seRNAs remain inactive in non-target cells and induce translation by partial degradation only in preselected cell types of interest. Cell specificity and type of functionalization are easily adaptable based on a modular system. In proof-of-concept studies we use seRNAs as platform technology for highly selective cell targeting. We effectively treat breast tumor cell clusters in mixed cell systems and shrink early U87 glioblastoma cell clusters in the brain of male mice without detectable side effects. Our data open up potential avenues for various therapeutic applications.
Background Direct RNA sequencing (dRNA-seq) on the Oxford Nanopore Technologies (ONT) platforms can produce reads covering up to full-length gene transcripts, while containing decipherable information about RNA base modifications and poly-A tail lengths. Although many published studies have been expanding the potential of dRNA-seq, its sequencing accuracy and error patterns remain understudied. Results We present the first comprehensive evaluation of sequencing accuracy and characterisation of systematic errors in dRNA-seq data from diverse organisms and synthetic in vitro transcribed RNAs. We found that for sequencing kits SQK-RNA001 and SQK-RNA002, the median read accuracy ranged from 87% to 92% across species, and deletions significantly outnumbered mismatches and insertions. Due to their high abundance in the transcriptome, heteropolymers and short homopolymers were the major contributors to the overall sequencing errors. We also observed systematic biases across all species at the levels of single nucleotides and motifs. In general, cytosine/uracil-rich regions were more likely to be erroneous than guanines and adenines. By examining raw signal data, we identified the underlying signal-level features potentially associated with the error patterns and their dependency on sequence contexts. While read quality scores can be used to approximate error rates at base and read levels, failure to detect DNA adapters may be a source of errors and data loss. By comparing distinct basecallers, we reason that some sequencing errors are attributable to signal insufficiency rather than algorithmic (basecalling) artefacts. Lastly, we generated dRNA-seq data using the latest SQK-RNA004 sequencing kit released at the end of 2023 and found that although the overall read accuracy increased, the systematic errors remain largely identical compared to the previous kits. Conclusions As the first systematic investigation of dRNA-seq errors, this study offers a comprehensive overview of reproducible error patterns across diverse datasets, identifies potential signal-level insufficiency, and lays the foundation for error correction methods.
To better understand viral pathogenesis, host-virus interactions, and potential therapeutic interventions, the development of robust reverse genetics systems for SARS-CoV-2 is crucial. Here, we present a reverse genetics platform that enables the efficient manipulation, assembly, and rescue of recombinant SARS-CoV-2. The versatility of our reverse genetics system was demonstrated by generating recombinant SARS-CoV-2 viruses. We used this system to generate N501Y and Y453F spike protein mutants. Characterization studies revealed distinct phenotypic effects, impact on viral fitness, cell binding, and replication kinetics. We also investigated a recently discovered priming site for NSP9, which is postulated to produce a short RNA antisense leader sequence. By introducing the U76G mutation into the 5’UTR, we show that this priming site is necessary for the correct production of genomic and subgenomic RNAs, and also for efficient viral replication. In conclusion, our developed reverse genetics system provides a robust and adaptable platform for the efficient generation of recombinant SARS-CoV-2 viruses for their comprehensive characterization.Significance statement In this study, we present a versatile reverse genetics platform facilitating the efficient manipulation, assembly, and rescue of recombinant SARS-CoV-2. Demonstrating its adaptability, we successfully engineered N501Y and Y453F spike protein mutants, each exhibiting distinct phenotypic effects on viral fitness, cell binding, and replication kinetics. We also investigated a novel negative sense priming site for NSP9, demonstrating a role in RNA production and viral replication. This straightforward reverse genetic system is therefore a powerful tool to generate recombinant viruses for advancing our understanding of SARS-CoV-2 biology.### Competing Interest StatementThe authors have declared no competing interest.
How do segmented RNA viruses correctly recruit their genome has yet to be clarified. Bluetongue virus is a double-stranded RNA virus with 10 segments of different sizes, but it assembles its genome in single-stranded form through a series of specific RNA-RNA interactions prior to packaging. In this study, we determined the structure of each BTV transcript, individually and in different combinations, using 2'-hydroxyl acylation analysed by primer extension and mutational profiling (SHAPE-MaP). SHAPE-MaP identified RNA structural changes during complex formation and putative RNA-RNA interaction sites. Our data also revealed a core RNA-complex of smaller segments which serves as the foundation ('anchor') for the assembly of a complete network composed of ten ssRNA segments. The same order of core RNA complex formation was identified in cells transfected with viral RNAs. No viral protein was required for these assembly reactions. Further, substitution mutations in the interacting bases within the core assemblies, altered subsequent segment addition and affected virus replication. These data identify a wholly RNA driven reaction that may offer novel opportunities for designed attenuation or antiviral therapeutics.
RNA structure determination is essential to understand how RNA carries out its diverse biological functions. In cells, RNA isoforms are readily expressed with partial variations within their sequences due, for example, to alternative splicing, heterogeneity in the transcription start site, RNA processing or differential termination/polyadenylation. Nanopore dimethyl sulfate mutational profiling (Nano-DMS-MaP) is a method for in situ isoform-specific RNA structure determination. Unlike similar methods that rely on short sequencing reads, Nano-DMS-MaP employs nanopore sequencing to resolve the structures of long and highly similar RNA molecules to reveal their previously hidden structural differences. This Protocol describes the development and applications of Nano-DMS-MaP and outlines the main considerations for designing and implementing a successful experiment: from bench to data analysis. In cell probing experiments can be carried out by an experienced molecular biologist in 3–4 d. Data analysis requires good knowledge of command line tools and Python scripts and requires a further 3–5 d.
A limited coding capacity means that RNA viruses need strategies to diversify their proteome. The nucleoprotein (NP) of the highly pathogenic arenavirus Junín virus (JUNV) produces three N-terminally truncated isoforms: two (NP 47kD and NP 40kD ) are known to be produced by caspase cleavage, while, here, we show that NP 53kD is produced by alternative translation initiation. Recombinant JUNVs lacking individual NP isoforms revealed that all three isoforms contribute to inhibiting caspase activation during infection, but cleavage to generate NP 40kD makes the biggest contribution. Importantly, all three isoforms retain their ability to digest double-stranded (ds)RNA and inhibit interferon promoter activation but have a diffuse cytoplasmic distribution. Given the cytoplasmic localization of both aberrant viral dsRNAs, as well as dsRNA sensors and many other cellular components of innate immune activation pathways, we suggest that the generation of NP isoforms not only contributes to evasion of apoptosis but also robust control of the antiviral response.
Long noncoding (lnc)RNAs emerge as regulators of genome stability. The nuclear-enriched abundant transcript 1 (NEAT1) is overexpressed in many tumors and is responsive to genotoxic stress. However, the mechanism that links NEAT1 to DNA damage response (DDR) is unclear. Here, we investigate the expression, modification, localization, and structure of NEAT1 in response to DNA double-strand breaks (DSBs). DNA damage increases the levels and N6-methyladenosine (m6A) marks on NEAT1, which promotes alterations in NEAT1 structure, accumulation of hypermethylated NEAT1 at promoter-associated DSBs, and DSB signaling. The depletion of NEAT1 impairs DSB focus formation and elevates DNA damage. The genome-protective role of NEAT1 is mediated by the RNA methyltransferase 3 (METTL3) and involves the release of the chromodomain helicase DNA binding protein 4 (CHD4) from NEAT1 to fine-tune histone acetylation at DSBs. Our data suggest a direct role for NEAT1 in DDR.
Genome-wide measurements of RNA structure can be obtained using reagents that react with unpaired bases, leading to adducts that can be identified by mutational profiling on next-generation sequencing machines. One drawback of these experiments is that short sequencing reads can rarely be mapped to specific transcript isoforms. Consequently, information is acquired as a population average in regions that are shared between transcripts, thus blurring the underlying structural landscape. Here, we present nanopore dimethylsulfate mutational profiling (Nano-DMS-MaP)—a method that exploits long-read sequencing to provide isoform-resolved structural information of highly similar RNA molecules. We demonstrate the value of Nano-DMS-MaP by resolving the complex structural landscape of human immunodeficiency virus-1 transcripts in infected cells. We show that unspliced and spliced transcripts have distinct structures at the packaging site within the common 5′ untranslated region, likely explaining why spliced viral RNAs are excluded from viral particles. Thus, Nano-DMS-MaP is a straightforward method to resolve biologically important transcript-specific RNA structures that were previously hidden in short-read ensemble analyses.
Nucleic acid sequencing is the process of identifying the sequence of DNA or RNA, with DNA used for genomes and RNA for transcriptomes. Deciphering this information has the potential to greatly advance our understanding of genomic features and cellular functions. In comparison to other available sequencing methods, nanopore sequencing stands out due to its unique advantages of processing long nucleic acid strands in real time, within a small portable device, enabling the rapid analysis of samples in diverse settings. Evolving over the past decade, nanopore sequencing remains in a state of ongoing development and refinement, resulting in persistent challenges in protocols and technology. This article employs an interdisciplinary approach, evaluating experimental and computational methods to address critical gaps in our understanding in order to maximise the information gain from this advancing technology. We present a robust analysis of all aspects of nanopore sequencing by providing statistically supported insights, thus aiming to provide comprehensive guidelines for the diverse challenges that frequently impede optimal experimental outcomes.Here we present a robust analysis, bridging the gap by providing statistically supported insights into genomic and transcriptomic studies, providing fresh perspectives on sequencing.### Competing Interest StatementDM, CD and BS are employees of oncgnostics GmbH, a company that aims to commercialize DNA methylation markers.
RNA dimerization is the noncovalent association of two human immunodeficiency virus-1 (HIV-1) genomes. It is a conserved step in the HIV-1 life cycle and assumed to be a prerequisite for binding to the viral structural protein Pr55Gag during genome packaging. Here, we developed functional analysis of RNA structure-sequencing (FARS-seq) to comprehensively identify sequences and structures within the HIV-1 5' untranslated region (UTR) that regulate this critical step. Using FARS-seq, we found nucleotides important for dimerization throughout the HIV-1 5' UTR and identified distinct structural conformations in monomeric and dimeric RNA. In the dimeric RNA, key functional domains, such as stem-loop 1 (SL1), polyadenylation signal (polyA) and primer binding site (PBS), folded into independent structural motifs. In the monomeric RNA, SL1 was reconfigured into long- and short-range base pairings with polyA and PBS, respectively. We show that these interactions disrupt genome packaging, and additionally show that the PBS-SL1 interaction unexpectedly couples the PBS with dimerization and Pr55Gag binding. Altogether, our data provide insights into late stages of HIV-1 life cycle and a mechanistic explanation for the link between RNA dimerization and packaging.
The arenavirus nucleoprotein (NP) plays an important role in the virus’ ability to block interferon (IFN) production, and its exonuclease function appears to contribute to this activity. However, efforts to analyze this contribution are complicated by the functional overlap between the exonuclease active site and a neighboring region involved in IKKε-binding and subsequent inhibition of IRF3 activation, which also plays an important role in IFN production. To circumvent this issue, we mutated a residue located away from the active site that is involved in binding of the dsRNA substrate being targeted for exonuclease digestion, i.e. H426A. We found that expression of Tacaribe virus (TCRV) NP containing this RNA-binding H426A mutation was still able to efficiently block IFN-β promoter activity in response to Sendai virus infection, despite being strongly impaired in its exonuclease activity. This was in contrast to a conventional exonuclease active site mutant (E388A), which was impaired with respect to both exonuclease activity and IFN antagonism. Importantly, growth of a recombinant virus encoding the RNA-binding mutation (rTCRV-H426A) was similar to wild-type in IFN-deficient cells, unlike the active site mutant (rTCRV-E388A), which was already markedly impaired in these cells. Further, in IFN-competent cells, the TCRV-H426A RNA-binding mutant showed more robust growth and delayed IFN-β mRNA upregulation compared to the TCRV-E388A active site mutant. Taken together, this novel mutational approach, which allows us to now dissect the different contributions of the NP exonuclease activity and IKKε-binding/IRF3 inhibition to IFN antagonism, clearly suggests that conventional exonuclease mutants targeting the active site overestimate the contribution of the exonuclease function, and that rather other IFN antagonistic functions of NP play the dominant role in IFN-antagonism.
The RNA genome of SARS-CoV-2 contains a frameshift stimulatory element (FSE) that allows access to an alternative reading frame through -1 programmed ribosomal frameshifting (PRF). -1PRF in the 1a/1b gene is essential for efficient viral replication and transcription of the viral genome. -1PRF efficiency relies on the presence of conserved RNA elements within the FSE. One of these elements is a three-stemmed pseudoknot, although alternative folds of the frameshift site might have functional roles as well. Here, by complementing ensemble and single-molecule structural analysis of SARS-CoV-2 frameshift RNA variants with functional data, we reveal a conformational interplay of the 5 ' and 3 ' immediate regions with the FSE and show that the extended FSE exists in multiple conformations. Furthermore, limiting the base pairing of the FSE with neighboring nucleotides can favor or impair the formation of the alternative folds, including the pseudoknot. Our results demonstrate that co-existing RNA structures can function together to fine-tune SARS-CoV-2 gene expression, which will aid efforts to design specific inhibitors of viral frameshifting.
To generate infectious viral particles, viruses must specifically select their genomic RNA from milieu that contains a complex mixture of cellular or non-genomic viral RNAs. In this review, we focus on the role of viral encoded RNA structures in genome packaging. We first discuss how packaging signals are constructed from local and long-range base pairings within viral genomes, as well as inter-molecular interactions between viral and host RNAs. Then, how genome packaging is regulated by the biophysical properties of RNA. Finally, we examine the impact of RNA packaging signals on viral evolution.
Macrophages have important protective functions during infection with herpes simplex virus type 1 (HSV-1). However, molecular mechanisms that restrict viral propagation and protect from severe disease are unclear. Here we show that macrophages take up HSV-1 via endocytosis and transport the virions into multivesicular bodies (MVBs). In MVBs, acid ceramidase (aCDase) converts ceramide into sphingosine and increases the formation of sphingosine-rich intraluminal vesicles (ILVs). Once HSV-1 particles reach MVBs, sphingosine-rich ILVs bind to HSV-1 particles, which restricts fusion with the limiting endosomal membrane and prevents cellular infection. Lack of aCDase in macrophage cultures or in Asah1 −/− mice results in replication of HSV-1 and Asah1 −/− mice die soon after systemic or intravaginal inoculation. The treatment of macrophages with sphingosine enhancing compounds blocks HSV-1 propagation, suggesting a therapeutic potential of this pathway. In conclusion, aCDase loads ILVs with sphingosine, which prevents HSV-1 capsids from penetrating into the cytosol.
Arenaviruses cause several viral hemorrhagic fevers endemic to Africa and South America. The respective causative agents are classified as biosafety level (BSL) 4 pathogens. Unlike for most other BSL4 agents, for the New World arenavirus Junín virus (JUNV) both a highly effective vaccination (Candid#1) and a post-exposure treatment, based on convalescent plasma transfer, are available. In particular, neutralizing antibodies (nAbs) represent a key protective determinant in JUNV infection, which is supported by the correlation between successful passive antibody therapy and the levels of nAbs administered. Unfortunately, comparable resources for the management of other closely related arenavirus infections are not available. Given the significant challenges inherent in studying BSL4 pathogens, our goal was to first assess the suitability of a JUNV transcription and replication-competent virus-like particle (trVLP) system for measuring virus neutralization under BSL1/2 conditions. Indeed, we could show that infection with JUNV trVLPs is glycoprotein (GP) dependent, that trVLP input has a direct correlation to reporter readout, and that these trVLPs can be neutralized by human serum with kinetics similar to those obtained using authentic virus. These properties make trVLPs suitable for use as a proxy for virus in neutralization assays. Using this platform we then evaluated the potential of JUNV nAbs to cross-neutralize entry mediated by GPs from other arenaviruses using JUNV (strain Romero)-based trVLPs bearing GPs either from other JUNV strains, other closely related New World arenaviruses (e.g. Tacaribe, Machupo, Sabiá), or the distantly related Lassa virus. While nAbs against the JUNV vaccine strain are also active against a range of other JUNV strains, they appear to have little or no capacity to neutralize other arenavirus species, suggesting that therapy with whole plasma directed against another species is unlikely to be successful and that the targeted development of cross-specific monoclonal antibody-based resources is likely needed. Such efforts will be supported by the availability of this BSL1/2 screening platform which provides a rapid and easy means to characterize the potency and reactivity of anti-arenavirus neutralizing antibodies against a range of arenavirus species.