Web-based platform for analysis of RNA folding from high throughput chemical probing data

RNA structures play critical roles in regulating gene expression across all domains of life and viruses. Chemical probing methods coupled with massively parallel sequencing have revolutionized the RNA structure field by enabling the assessment of many structures in their native, physiological context. Previously, we developed Dimethyl-Sulfate-based Mutational Profiling and Sequencing (DMS-MaPseq), which uses DMS to label the Watson-Crick face of open and accessible adenine and cytosine bases in the RNA. We used this approach to determine the genome-wide structures of HIV-1 and SARS-CoV-2 in infected cells, which permitted uncovering new biology and identifying therapeutic targets. Due to the simplicity and ease of the experimental procedure, DMS-MaPseq has been adopted by labs worldwide. However, bioinformatic analysis remains a substantial hurdle for labs that often lack the necessary infrastructure and computational expertise. Here we present a modern web-based interface that automates the analysis of chemical probing profiles from raw sequencing files (http://rnadreem.org). The availability of a simple web-based platform for DMSMaPseq analysis will dramatically expand studies of RNA structure and aid in the design of structurebased therapeutics

Genome-wide probing of RNA structure reveals active unfolding of mRNA structures in vivo

RNA has a dual role as an informational molecule and a direct effector of biological tasks. The latter function is enabled by RNA’s ability to adopt complex secondary and tertiary folds and thus has motivated extensive computational and experimental efforts for determining RNA structures. Existing approaches for evaluating RNA structure have been largely limited to in vitro systems, yet the thermodynamic forces which drive RNA folding in vitro may not be sufficient to predict stable RNA structures in vivo. Indeed, the presence of RNA-binding proteins and ATP-dependent helicases can influence which structures are present inside cells. Here we present an approach for globally monitoring RNA structure in native conditions in vivo with single-nucleotide precision. This method is based on in vivo modification with dimethyl sulphate (DMS), which reacts with unpaired adenine and cytosine residues, followed by deep sequencing to monitor modifications. Our data from yeast and mammalian cells are in excellent agreement with known messenger RNA structures and with the high-resolution crystal structure of the Saccharomyces cerevisiae ribosome. Comparison between in vivo and in vitro data reveals that in rapidly dividing cells there are vastly fewer structured mRNA regions in vivo than in vitro. Even thermostable RNA structures are often denatured in cells, highlighting the importance of cellular processes in regulating RNA structure. Indeed, analysis of mRNA structure under ATP-depleted conditions in yeast shows that energy-dependent processes strongly contribute to the predominantly unfolded state of mRNAs inside cells. Our studies broadly enable the functional analysis of physiological RNA structures and reveal that, in contrast to the Anfinsen view of protein folding whereby the structure formed is the most thermodynamically favourable, thermodynamics have an incomplete role in determining mRNA structure in vivo

Microarray detection of human parainfluenzavirus 4 infection associated with respiratory failure in an immunocompetent adult.

A pan-viral DNA microarray, the Virochip (University of California, San Francisco), was used to detect human parainfluenzavirus 4 (HPIV-4) infection in an immunocompetent adult presenting with a life-threatening acute respiratory illness. The virus was identified in an endotracheal aspirate specimen, and the microarray results were confirmed by specific polymerase chain reaction and serological analysis for HPIV-4. Conventional clinical laboratory testing using an extensive panel of microbiological tests failed to yield a diagnosis. This case suggests that the potential severity of disease caused by HPIV-4 in adults may be greater than previously appreciated and illustrates the clinical utility of a microarray for broad-based viral pathogen screening

Operon mRNAs are organized into ORF-centric structures that predict translation efficiency

Bacterial mRNAs are organized into operons consisting of discrete open reading frames (ORFs) in a single polycistronic mRNA. Individual ORFs on the mRNA are differentially translated, with rates varying as much as 100-fold. The signals controlling differential translation are poorly understood. Our genome-wide mRNA secondary structure analysis indicated that operonic mRNAs are comprised of ORF-wide units of secondary structure that vary across ORF boundaries such that adjacent ORFs on the same mRNA molecule are structurally distinct. ORF translation rate is strongly correlated with its mRNA structure in vivo, and correlation persists, albeit in a reduced form, with its structure when translation is inhibited and with that of in vitro refolded mRNA. These data suggest that intrinsic ORF mRNA structure encodes a rough blueprint for translation efficiency. This structure is then amplified by translation, in a self-reinforcing loop, to provide the structure that ultimately specifies the translation of each ORF

An intranasal ASO therapeutic targeting SARS-CoV-2

The COVID-19 pandemic is exacting an increasing toll worldwide, with new SARS-CoV-2 variants emerging that exhibit higher infectivity rates and that may partially evade vaccine and antibody immunity. Rapid deployment of non-invasive therapeutic avenues capable of preventing infection by all SARS-CoV-2 variants could complement current vaccination efforts and help turn the tide on the COVID-19 pandemic. Here, we describe a novel therapeutic strategy targeting the SARS-CoV-2 RNA using locked nucleic acid antisense oligonucleotides (LNA ASOs). We identify an LNA ASO binding to the 5′ leader sequence of SARS-CoV-2 that disrupts a highly conserved stem-loop structure with nanomolar efficacy in preventing viral replication in human cells. Daily intranasal administration of this LNA ASO in the COVID-19 mouse model potently suppresses viral replication (>80-fold) in the lungs of infected mice. We find that the LNA ASO is efficacious in countering all SARS-CoV-2 “variants of concern” tested both in vitro and in vivo. Hence, inhaled LNA ASOs targeting SARS-CoV-2 represents a promising therapeutic approach to reduce or prevent transmission and decrease severity of COVID-19 in infected individuals. LNA ASOs are chemically stable and can be flexibly modified to target different viral RNA sequences and could be stockpiled for future coronavirus pandemics

Genome-wide diversity and selective pressure in the human rhinovirus

BACKGROUND: The human rhinoviruses (HRV) are one of the most common and diverse respiratory pathogens of humans. Over 100 distinct HRV serotypes are known, yet only 6 genomes are available. Due to the paucity of HRV genome sequence, little is known about the genetic diversity within HRV or the forces driving this diversity. Previous comparative genome sequence analyses indicate that recombination drives diversification in multiple genera of the picornavirus family, yet it remains unclear if this holds for HRV. RESULTS: To resolve this and gain insight into the forces driving diversification in HRV, we generated a representative set of 34 fully sequenced HRVs. Analysis of these genomes shows consistent phylogenies across the genome, conserved non-coding elements, and only limited recombination. However, spikes of genetic diversity at both the nucleotide and amino acid level are detectable within every locus of the genome. Despite this, the HRV genome as a whole is under purifying selective pressure, with islands of diversifying pressure in the VP1, VP2, and VP3 structural genes and two non-structural genes, the 3C protease and 3D polymerase. Mapping diversifying residues in these factors onto available 3-dimensional structures revealed the diversifying capsid residues partition to the external surface of the viral particle in statistically significant proximity to antigenic sites. Diversifying pressure in the pleconaril binding site is confined to a single residue known to confer drug resistance (VP1 191). In contrast, diversifying pressure in the non-structural genes is less clear, mapping both nearby and beyond characterized functional domains of these factors. CONCLUSION: This work provides a foundation for understanding HRV genetic diversity and insight into the underlying biology driving evolution in HRV. It expands our knowledge of the genome sequence space that HRV reference serotypes occupy and how the pattern of genetic diversity across HRV genomes differs from other picornaviruses. It also reveals evidence of diversifying selective pressure in both structural genes known to interact with the host immune system and in domains of unassigned function in the non-structural 3C and 3D genes, raising the possibility that diversification of undiscovered functions in these essential factors may influence HRV fitness and evolution

COMRADES determines in vivo RNA structures and interactions.

The structural flexibility of RNA underlies fundamental biological processes, but there are no methods for exploring the multiple conformations adopted by RNAs in vivo. We developed cross-linking of matched RNAs and deep sequencing (COMRADES) for in-depth RNA conformation capture, and a pipeline for the retrieval of RNA structural ensembles. Using COMRADES, we determined the architecture of the Zika virus RNA genome inside cells, and identified multiple site-specific interactions with human noncoding RNAs.This work was supported by Cancer Research UK (C13474/A18583, C6946/A14492) and the Wellcome Trust (104640/Z/14/Z, 092096/Z/10/Z) to E.A.M. O.Z. was supported by the Human Frontier Science Program (HFSP, LT000558/2015), the European Molecular Biology Organization (EMBO, ALTF1622-2014), and the Blavatnik Family Foundation postdoctoral fellowship. G.K. and M.G. were supported by Wellcome Trust grant 207507 and UK Medical Research Council. A.T.L.L. and J.C.M. were supported by core funding from Cancer Research UK (award no. 17197 to JCM). J.C.M was also supported by core funding from EMBL. I.G. and L.W.M. were supported by the Wellcome Trust Senior Fellowship in Basic Biomedical Science to I.G. (207498/Z/17/Z). I.J.M., L.F.G. and J.S.-G. were supported by grants R01GM104475 and R01GM115649 from NIGMS. C.K.K was supported by City University of Hong Kong Projects 9610363 and 7200520, Croucher Foundation Project 9500030 and Hong Kong RGC Projects 9048103 and 9054020. C.-F.Q. was supported by the NSFC Excellent Young Scientist Fund 81522025 and the Newton Advanced Fellowship from the Academy of Medical Sciences, UK

Deciphering the universe of RNA structures and trans RNA-RNA interactions of transcriptomes in vivo: from experimental protocols to computational analyses

The last few years have seen an explosion of experimental and computational methods for investigating RNA structures of entire transcriptomes in vivo. Very recent experimental protocols now also allow trans RNA–RNA interactions to be probed in a transcriptome-wide manner. All of the experimental strategies require comprehensive computational pipelines for analysing the raw data and converting it back into actual RNA structure features or trans RNA–RNA interactions. The overall performance of these methods thus strongly depends on the experimental and the computational protocols employed. In order to get the best out of both worlds, both aspects need to be optimised simultaneously. This review introduced the methods and proposes ideas how they could be further improved