The influence of viral RNA secondary structure on interactions with innate host cell defences

Raman Fingerprints of SARS-CoV-2 Omicron Subvariants: Molecular Roots of Virological Characteristics and Evolutionary Directions

The latest RNA genomic mutation of SARS-CoV-2 virus, termed the Omicron variant, has generated a stream of highly contagious and antibody-resistant strains, which in turn led to classifying Omicron as a variant of concern. We systematically collected Raman spectra from six Omicron subvariants available in Japan (i.e., BA.1.18, BA.2, BA.4, BA.5, XE, and BA.2.75) and applied machine-learning algorithms to decrypt their structural characteristics at the molecular scale. Unique Raman fingerprints of sulfur-containing amino acid rotamers, RNA purines and pyrimidines, tyrosine phenol ring configurations, and secondary protein structures clearly differentiated the six Omicron subvariants. These spectral characteristics, which were linked to infectiousness, transmissibility, and propensity for immune evasion, revealed evolutionary motifs to be compared with the outputs of genomic studies. The availability of a Raman "metabolomic snapshot", which was then translated into a barcode to enable a prompt subvariant identification, opened the way to rationalize in real-time SARS-CoV-2 activity and variability. As a proof of concept, we applied the Raman barcode procedure to a nasal swab sample retrieved from a SARS-CoV-2 patient and identified its Omicron subvariant by coupling a commercially available magnetic bead technology with our newly developed Raman analyses.

Transcriptome Profiling of Vero E6 Cells during Original Parental or Cell-Attenuated Porcine Epidemic Diarrhea Virus Infection

Porcine epidemic diarrhea virus (PEDV) has led to significant economic losses in the global porcine industry since the emergence of variant strains in 2010. The high mutability of coronaviruses endows PEDV with the ability to evade the host immune response, which impairs the effectiveness of vaccines. In our previous study, we generated a highly cell-passaged PEDV strain, CT-P120, which showed promise as a live attenuated vaccine candidate by providing satisfactory protection against variant PEDV infection in piglets. However, the mechanism by which the attenuated CT-P120 adapts to cells during passage, resulting in increased replication efficiency, remains unclear. To address this question, we conducted a comparative transcriptomic analysis of Vero E6 cells infected with either the original parental strain (CT-P10) or the cell-attenuated strain (CT-P120) of PEDV at 6, 12, and 24 h post-infection. Compared to CT-P10, CT-P120 infection resulted in a significant decrease in the number of differentially expressed genes (DEGs) at each time point. Functional enrichment analysis of genes revealed the activation of various innate immune-related pathways by CT-P10, notably attenuated during CT-P120 infection. To validate these results, we selected eight genes (TRAF3, IRF3, IFNL1, ISG15, NFKB1, MAP2K3, IL1A, and CCL2) involved in antiviral processes and confirmed their mRNA expression patterns using RT-qPCR, in line with the transcriptomic data. Subsequent protein-level analysis of selected genes via Western blotting and enzyme-linked immunosorbent assay corroborated these results, reinforcing the robustness of our findings. Collectively, our research elucidates the strategies underpinning PEDV attenuation and immune evasion, providing invaluable insights for the development of effective PEDV vaccines.

Epistasis-Driven Evolution of the SARS-CoV-2 Secondary Structure

Epistasis is an evolutionary phenomenon whereby the fitness effect of a mutation depends on the genetic background in which it arises. A key source of epistasis in an RNA molecule is its secondary structure, which contains functionally important topological motifs held together by hydrogen bonds between Watson–Crick (WC) base pairs. Here we study epistasis in the secondary structure of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by examining properties of derived alleles arising from substitution mutations at ancestral WC base-paired and unpaired (UP) sites in 15 conserved topological motifs across the genome. We uncover fewer derived alleles and lower derived allele frequencies at WC than at UP sites, supporting the hypothesis that modifications to the secondary structure are often deleterious. At WC sites, we also find lower derived allele frequencies for mutations that abolish base pairing than for those that yield G·U “wobbles,” illustrating that weak base pairing can partially preserve the integrity of the secondary structure. Last, we show that WC sites under the strongest epistatic constraint reside in a three-stemmed pseudoknot motif that plays an essential role in programmed ribosomal frameshifting, whereas those under the weakest epistatic constraint are located in 3’ UTR motifs that regulate viral replication and pathogenicity. Our findings demonstrate the importance of epistasis in the evolution of the SARS-CoV-2 secondary structure, as well as highlight putative structural and functional targets of different forms of natural selection.

Lessons Learned and Yet-to-Be Learned on the Importance of RNA Structure in SARS-CoV-2 Replication

SARS-CoV-2, the etiological agent responsible for the COVID-19 pandemic, is a member of the virus family Coronaviridae, known for relatively extensive (~30-kb) RNA genomes that not only encode for numerous proteins but are also capable of forming elaborate structures. As highlighted in this review, these structures perform critical functions in various steps of the viral life cycle, ultimately impacting pathogenesis and transmissibility. We examine these elements in the context of coronavirus evolutionary history and future directions for curbing the spread of SARS-CoV-2 and other potential human coronaviruses. While we focus on structures supported by a variety of biochemical, biophysical, and/or computational methods, we also touch here on recent evidence for novel structures in both protein-coding and noncoding regions of the genome, including an assessment of the potential role for RNA structure in the controversial finding of SARS-CoV-2 integration in “long COVID” patients. This review aims to serve as a consolidation of previous works on coronavirus and more recent investigation of SARS-CoV-2, emphasizing the need for improved understanding of the role of RNA structure in the evolution and adaptation of these human viruses.

The In Vivo and In Vitro Architecture of the Hepatitis C Virus RNA Genome Uncovers Functional RNA Secondary and Tertiary Structures

Hepatitis C virus (HCV) is a positive-strand RNA virus that remains one of the main contributors to chronic liver disease worldwide. Studies over the last 30 years have demonstrated that HCV contains a highly structured RNA genome and many of these structures play essential roles in the HCV life cycle. Despite the importance of riboregulation in this virus, most of the HCV RNA genome remains functionally unstudied. Here, we report a complete secondary structure map of the HCV RNA genome in vivo, which was studied in parallel with the secondary structure of the same RNA obtained in vitro. Our results show that HCV is folded extensively in the cellular context. By performing comprehensive structural analyses on both in vivo data and in vitro data, we identify compact and conserved secondary and tertiary structures throughout the genome. Genetic and evolutionary functional analyses demonstrate that many of these elements play important roles in the virus life cycle. In addition to providing a comprehensive map of RNA structures and riboregulatory elements in HCV, this work provides a resource for future studies aimed at identifying therapeutic targets and conducting further mechanistic studies on this important human pathogen. IMPORTANCE HCV has one of the most highly structured RNA genomes studied to date, and it is a valuable model system for studying the role of RNA structure in protein-coding genes. While previous studies have identified individual cases of regulatory RNA structures within the HCV genome, the full-length structure of the HCV genome has not been determined in vivo. Here, we present the complete secondary structure map of HCV determined both in cells and from corresponding transcripts generated in vitro. In addition to providing a comprehensive atlas of functional secondary structural elements throughout the genomic RNA, we identified a novel set of tertiary interactions and demonstrated their functional importance. In terms of broader implications, the pipeline developed in this study can be applied to other long RNAs, such as long noncoding RNAs. In addition, the RNA structural motifs characterized in this study broaden the repertoire of known riboregulatory elements.

Conformational dynamics and energetics of viral RNA recognition by lab-evolved proteins

The conserved and structured elements in viral RNA genomes interact with proteins to regulate various events in the viral life cycle and have become key targets for developing novel therapeutic approaches. We probe physical interactions between lab-evolved proteins and a viral RNA element from the HIV-1 genome. Specifically, we study the role of an arginine-rich loop in recognition of designed proteins by the viral RNA element. We report free energy calculations to quantitatively estimate the protein/RNA binding energetics, focusing on the mutations of arginine residues involved in recognition of the major groove of RNA by proteins.

G-quadruplexes in genomes of viruses infecting eukaryotes or prokaryotes are under different selection pressures from hosts

G-quadruplexes in viral genomes can be applied as the targets of antiviral therapies, which has attracted wide interest. However, it is still not clear whether the pervasive number of such elements in the viral world is the result of natural selection for functionality. In this study, we identified putative quadruplex-forming sequences (PQSs) across the known viral genomes and analyzed the abundance, structural stability, and conservation of viral PQSs. A Viral Putative G-quadruplex Database (ViPGD,http://jsjds.hzau.edu.cn/MBPC/ViPGD/index.php/home/index) was constructed to collect the details of each viral PQS, which provides guidance for selecting the desirable PQS. The PQS with two putative G-tetrads (G₂-PQS) was significantly enriched in both eukaryotic viruses and prokaryotic viruses, while the PQSs with three putative G-tetrads (G₃-PQS) were only enriched in eukaryotic viruses and depleted in prokaryotic viruses. The structural stability of PQSs in prokaryotic viruses was significantly lower than that in eukaryotic viruses. Conservation analysis showed that the G₂-PQS, instead of G₃-PQS, was highly conserved within the genus. This suggested that the G₂-quadruplex might play an important role in viral biology, and the difference in the occurrence of G-quadruplex between eukaryotic viruses and prokaryotic viruses may result from the different selection pressures from hosts.

Mutagenesis Mapping of RNA Structures within the Foot-and-Mouth Disease Virus Genome Reveals Functional Elements Localized in the Polymerase (3Dpol)-Encoding Region

RNA structures can form functional elements that play crucial roles in the replication of positive-sense RNA viruses. While RNA structures in the untranslated regions (UTRs) of several picornaviruses have been functionally characterized, the roles of putative RNA structures predicted for protein coding sequences (or open reading frames [ORFs]) remain largely undefined. Here, we have undertaken a bioinformatic analysis of the foot-and-mouth disease virus (FMDV) genome to predict 53 conserved RNA structures within the ORF. Forty-six of these structures were located in the regions encoding the nonstructural proteins (nsps). To investigate whether structures located in the regions encoding the nsps are required for FMDV replication, we used a mutagenesis method, CDLR mapping, where sequential coding segments were shuffled to minimize RNA secondary structures while preserving protein coding, native dinucleotide frequencies, and codon usage. To examine the impact of these changes on replicative fitness, mutated sequences were inserted into an FMDV subgenomic replicon. We found that three of the RNA structures, all at the 3' termini of the FMDV ORF, were critical for replicon replication. In contrast, disruption of the other 43 conserved RNA structures that lie within the regions encoding the nsps had no effect on replicon replication, suggesting that these structures are not required for initiating translation or replication of viral RNA. Conserved RNA structures that are not essential for virus replication could provide ideal targets for the rational attenuation of a wide range of FMDV strains. IMPORTANCE Some RNA structures formed by the genomes of RNA viruses are critical for viral replication. Our study shows that of 46 conserved RNA structures located within the regions of the foot-and-mouth disease virus (FMDV) genome that encode the nonstructural proteins, only three are essential for replication of an FMDV subgenomic replicon. Replicon replication is dependent on RNA translation and synthesis; thus, our results suggest that the three RNA structures are critical for either initiation of viral RNA translation and/or viral RNA synthesis. Although further studies are required to identify whether the remaining 43 RNA structures have other roles in virus replication, they may provide targets for the rational large-scale attenuation of a wide range of FMDV strains. FMDV causes a highly contagious disease, posing a constant threat to global livestock industries. Such weakened FMDV strains could be investigated as live-attenuated vaccines or could enhance biosecurity of conventional inactivated vaccine production.

Comparative Genomics Reveals Early Emergence and Biased Spatiotemporal Distribution of SARS-CoV-2

Effective systems for the analysis of molecular data are fundamental for monitoring the spread of infectious diseases and studying pathogen evolution. The rapid identification of emerging viral strains, and/or genetic variants potentially associated with novel phenotypic features is one of the most important objectives of genomic surveillance of human pathogens and represents one of the first lines of defense for the control of their spread. During the COVID 19 pandemic, several taxonomic frameworks have been proposed for the classification of SARS-Cov-2 isolates. These systems, which are typically based on phylogenetic approaches, represent essential tools for epidemiological studies as well as contributing to the study of the origin of the outbreak. Here, we propose an alternative, reproducible, and transparent phenetic method to study changes in SARS-CoV-2 genomic diversity over time. We suggest that our approach can complement other systems and facilitate the identification of biologically relevant variants in the viral genome. To demonstrate the validity of our approach, we present comparative genomic analyses of more than 175,000 genomes. Our method delineates 22 distinct SARS-CoV-2 haplogroups, which, based on the distribution of high-frequency genetic variants, fall into four major macrohaplogroups. We highlight biased spatiotemporal distributions of SARS-CoV-2 genetic profiles and show that seven of the 22 haplogroups (and of all of the four haplogroup clusters) showed a broad geographic distribution within China by the time the outbreak was widely recognized—suggesting early emergence and widespread cryptic circulation of the virus well before its isolation in January 2020. General patterns of genomic variability are remarkably similar within all major SARS-CoV-2 haplogroups, with UTRs consistently exhibiting the greatest variability, with s2m, a conserved secondary structure element of unknown function in the 3′-UTR of the viral genome showing evidence of a functional shift. Although several polymorphic sites that are specific to one or more haplogroups were predicted to be under positive or negative selection, overall our analyses suggest that the emergence of novel types is unlikely to be driven by convergent evolution and independent fixation of advantageous substitutions, or by selection of recombined strains. In the absence of extensive clinical metadata for most available genome sequences, and in the context of extensive geographic and temporal biases in the sampling, many questions regarding the evolution and clinical characteristics of SARS-CoV-2 isolates remain open. However, our data indicate that the approach outlined here can be usefully employed in the identification of candidate SARS-CoV-2 genetic variants of clinical and epidemiological importance.

K160 in the RNA-binding domain of the orf virus virulence factor OV20.0 is critical for its functions in counteracting host antiviral defense

The OV20.0 virulence factor of orf virus antagonizes host antiviral responses. One mechanism through which it functions is by inhibiting activation of the dsRNA-activated protein kinase R (PKR) by sequestering dsRNA and by physically interacting with PKR. Sequence alignment indicated that several key residues critical for dsRNA binding were conserved in OV20.0, and their contribution to OV20.O function was investigated in this study. We found that residues F141, K160, and R164 were responsible for the dsRNA-binding ability of OV20.0. Interestingly, mutation at K160 (K160A) diminished the OV20.0-PKR interaction and further reduced the inhibitory effect of OV20.0 on PKR activation. Nevertheless, OV20.0 homodimerization was not influenced by K160A. The contribution of the dsRNA-binding domain and K160 to the suppression of RNA interference by OV20.0 was further demonstrated in plants. In summary, K160 is essential for the function of OV20.0, particularly its interaction with dsRNA and PKR that ultimately contributes to the suppression of PKR activation.

Structural landscape of the complete genomes of dengue virus serotypes and other viral hemorrhagic fevers

BACKGROUND

With more than 300 million potentially infected people every year, and with the expanded habitat of mosquitoes due to climate change, Dengue virus (DENV) cannot be considered anymore only a tropical disease. The RNA secondary structure is a functional characteristic of RNA viruses, and together with the accumulated high-throughput sequencing data could provide general insights towards understanding virus biology. Here, we profiled the RNA secondary structure of > 7000 complete viral genomes from 11 different species focusing on viral hemorrhagic fevers, including DENV serotypes, EBOV, and YFV.

RESULTS

In our work we demonstrated that the secondary structure and presence of protein-binding domains in the genomes can be used as intrinsic signature to further classify the viruses. With our predictive approach, we achieved high prediction scores of the secondary structure (AUC up to 0.85 with experimental data), and computed consensus secondary structure profiles using hundreds of in silico models. We observed that viruses show different structural patterns, where e.g., DENV-2 and Ebola virus tend to be less structured than the other viruses. Furthermore, we observed virus-specific correlations between secondary structure and the number of interaction sites with human proteins, reaching a correlation of 0.89 in the case of Zika virus. We also identified that helicases-encoding regions are more structured in several flaviviruses, while the regions encoding for the contact proteins exhibit virus-specific clusters in terms of RNA structure and potential protein-RNA interactions. We also used structural data to study the geographical distribution of DENV, finding a significant difference between DENV-3 from Asia and South-America, where the structure is also driving the clustering more than sequence identity, which could imply different evolutionary routes of this subtype.

CONCLUSIONS

Our massive computational analysis provided novel results regarding the secondary structure and the interaction with human proteins, not only for DENV serotypes, but also for other flaviviruses and viral hemorrhagic fevers-associated viruses. We showed how the RNA secondary structure can be used to categorise viruses, and even to further classify them based on the interaction with proteins. We envision that these approaches can be used to further classify and characterise these complex viruses.

Neutralism versus selectionism: Chargaff's second parity rule, revisited

Of Chargaff's four "rules" on DNA base frequencies, the functional interpretation of his second parity rule (PR2) is the most contentious. Thermophile base compositions (GC%) were taken by Galtier and Lobry (1997) as favoring Sueoka's neutral PR2 hypothesis over Forsdyke's selective PR2 hypothesis, namely that mutations improving local within-species recombination efficiency had generated a genome-wide potential for the strands of duplex DNA to separate and initiate recombination through the "kissing" of the tips of stem-loops. However, following Chargaff's GC rule, base composition mainly reflects a species-specific, genome-wide, evolutionary pressure. GC% could not have consistently followed the dictates of temperature, since it plays fundamental roles in both sustaining species integrity and, through primarily neutral genome-wide mutation, fostering speciation. Evidence for a local within-species recombination-initiating role of base order was obtained with a novel technology that masked the contribution of base composition to nucleic acid folding energy. Forsdyke's results were consistent with his PR2 hypothesis, appeared to resolve some root problems in biology and provided a theoretical underpinning for alignment-free taxonomic analyses using relative oligonucleotide frequencies (k-mer analysis). Moreover, consistent with Chargaff's cluster rule, discovery of the thermoadaptive role of the "purine-loading" of open reading frames made less tenable the Galtier-Lobry anti-selectionist arguments.

Insect-Specific Flavivirus Replication in Mammalian Cells Is Inhibited by Physiological Temperature and the Zinc-Finger Antiviral Protein

The genus Flavivirus contains pathogenic vertebrate-infecting flaviviruses (VIFs) and insect-specific flaviviruses (ISF). ISF transmission to vertebrates is inhibited at multiple stages of the cellular infection cycle, via yet to be elucidated specific antiviral responses. The zinc-finger antiviral protein (ZAP) in vertebrate cells can bind CpG dinucleotides in viral RNA, limiting virus replication. Interestingly, the genomes of ISFs contain more CpG dinucleotides compared to VIFs. In this study, we investigated whether ZAP prevents two recently discovered lineage II ISFs, Binjari (BinJV) and Hidden Valley viruses (HVV) from replicating in vertebrate cells. BinJV protein and dsRNA replication intermediates were readily observed in human ZAP knockout cells when cultured at 34 °C. In ZAP-expressing cells, inhibition of the interferon response via interferon response factors 3/7 did not improve BinJV protein expression, whereas treatment with kinase inhibitor C16, known to reduce ZAP’s antiviral function, did. Importantly, at 34 °C, both BinJV and HVV successfully completed the infection cycle in human ZAP knockout cells evident from infectious progeny virus in the cell culture supernatant. Therefore, we identify vertebrate ZAP as an important barrier that protects vertebrate cells from ISF infection. This provides new insights into flavivirus evolution and the mechanisms associated with host switching.

The global and local distribution of RNA structure throughout the SARS-CoV-2 genome

SARS-CoV-2 is the causative viral agent of COVID-19, the disease at the center of the current global pandemic. While knowledge of highly structured regions is integral for mechanistic insights into the viral infection cycle, very little is known about the location and folding stability of functional elements within the massive, ∼30kb SARS-CoV-2 RNA genome. In this study, we analyze the folding stability of this RNA genome relative to the structural landscape of other well-known viral RNAs. We present an in-silico pipeline to predict regions of high base pair content across long genomes and to pinpoint hotspots of well-defined RNA structures, a method that allows for direct comparisons of RNA structural complexity within the several domains in SARS-CoV-2 genome. We report that the SARS-CoV-2 genomic propensity for stable RNA folding is exceptional among RNA viruses, superseding even that of HCV, one of the most structured viral RNAs in nature. Furthermore, our analysis suggests varying levels of RNA structure across genomic functional regions, with accessory and structural ORFs containing the highest structural density in the viral genome. Finally, we take a step further to examine how individual RNA structures formed by these ORFs are affected by the differences in genomic and subgenomic contexts, which given the technical difficulty of experimentally separating cellular mixtures of sgRNA from gRNA, is a unique advantage of our in-silico pipeline. The resulting findings provide a useful roadmap for planning focused empirical studies of SARS-CoV-2 RNA biology, and a preliminary guide for exploring potential SARS-CoV-2 RNA drug targets.Importance The RNA genome of SARS-CoV-2 is among the largest and most complex viral genomes, and yet its RNA structural features remain relatively unexplored. Since RNA elements guide function in most RNA viruses, and they represent potential drug targets, it is essential to chart the architectural features of SARS-CoV-2 and pinpoint regions that merit focused study. Here we show that RNA folding stability of SARS-CoV-2 genome is exceptional among viral genomes and we develop a method to directly compare levels of predicted secondary structure across SARS-CoV-2 domains. Remarkably, we find that coding regions display the highest structural propensity in the genome, forming motifs that differ between the genomic and subgenomic contexts. Our approach provides an attractive strategy to rapidly screen for candidate structured regions based on base pairing potential and provides a readily interpretable roadmap to guide functional studies of RNA viruses and other pharmacologically relevant RNA transcripts.

Metallohelices emulate the properties of short cationic α-helical peptides

Naturally occurring peptides in many living systems perform antimicrobial and anticancer host defence roles, but their potential for clinical application is limited by low metabolic stability and relatively high costs of goods. Self-assembled helical metal complexes provide an attractive synthetic platform for non-peptidic architectures that can emulate some of the properties of short cationic α-helical peptides, with tuneable charge, shape, size and amphipathicity. Correspondingly there is a growing body of evidence demonstrating that these supramolecular architectures exhibit bioactivity that emulates that of the natural systems. We review that evidence in the context of synthetic advances in the area, driven by the potential for biomedical applications. We note some design considerations for new biologically-relevant metallohelices, and give our outlook on the future of these compounds as therapeutic peptidomimetics.

New RNA Structural Elements Identified in the Coding Region of the Coxsackie B3 Virus Genome

Here we present a set of new structural elements formed within the open reading frame of the virus, which are highly probable, evolutionarily conserved and may interact with host proteins. This work focused on the coding regions of the CVB3 genome (particularly the V4-, V1-, 2C-, and 3D-coding regions), which, with the exception of the cis-acting replication element (CRE), have not yet been subjected to experimental analysis of their structures. The SHAPE technique, chemical modification with DMS and RNA cleavage with Pb²⁺, were performed in order to characterize the RNA structure. The experimental results were used to improve the computer prediction of the structural models, whereas a phylogenetic analysis was performed to check universality of the newly identified structural elements for twenty CVB3 genomes and 11 other enteroviruses. Some of the RNA motifs turned out to be conserved among different enteroviruses. We also observed that the 3'-terminal region of the genome tends to dimerize in a magnesium concentration-dependent manner. RNA affinity chromatography was used to confirm RNA-protein interactions hypothesized by database searches, leading to the discovery of several interactions, which may be important for virus propagation.

Norovirus infection results in eIF2α independent host translation shut-off and remodels the G3BP1 interactome evading stress granule formation

Viral infections impose major stress on the host cell. In response, stress pathways can rapidly deploy defence mechanisms by shutting off the protein synthesis machinery and triggering the accumulation of mRNAs into stress granules to limit the use of energy and nutrients. Because this threatens viral gene expression, viruses need to evade these pathways to propagate. Human norovirus is responsible for gastroenteritis outbreaks worldwide. Here we examined how norovirus interacts with the eIF2α signaling axis controlling translation and stress granules. While norovirus infection represses host cell translation, our mechanistic analyses revealed that eIF2α signaling mediated by the stress kinase GCN2 is uncoupled from translational stalling. Moreover, infection results in a redistribution of the RNA-binding protein G3BP1 to replication complexes and remodelling of its interacting partners, allowing the avoidance from canonical stress granules. These results define novel strategies by which norovirus undergo efficient replication whilst avoiding the host stress response and manipulating the G3BP1 interactome.

Viruses have evolved elegant strategies to evade host responses that restrict viral propagation by targeting the protein synthesis machinery and stress granules, which are membrane-less RNA granules with antiviral properties. Previous studies have unravelled how viruses, including norovirus the leading cause of gastroenteritis, regulate the activity of translation factors to affect the antiviral response. Furthermore, stress granules evasion strategies have been linked to targeting the scaffolding protein G3BP1. Here we dissect how murine norovirus, the main model for norovirus, evades the cellular stress responses. Our work challenges the dogma that translational control during infection is mainly mediated by eIF2α and demonstrate that norovirus evades this stress pathway. We further show that norovirus evades the stress granule response in a novel way by isolating and characterising the G3BP1 interactome for the first time in the context of a viral infection. We conclude that norovirus infection results in a redistribution of G3BP1 and its cellular partners to replication complexes, thereby preventing the assembly of stress granules. Overall, we define a novel evasion strategy by which norovirus escapes stress granule formation by rewiring the G3BP1 interactome.

Inhibition of Microprocessor Function during the Activation of the Type I Interferon Response

Type I interferons (IFNs) are central components of the antiviral response. Most cell types respond to viral infections by secreting IFNs, but the mechanisms that regulate correct expression of these cytokines are not completely understood. Here, we show that activation of the type I IFN response regulates the expression of miRNAs in a post-transcriptional manner. Activation of IFN expression alters the binding of the Microprocessor complex to pri-miRNAs, reducing its processing rate and thus leading to decreased levels of a subset of mature miRNAs in an IRF3-dependent manner. The rescue of Microprocessor function during the antiviral response downregulates the levels of IFN-β and IFN-stimulated genes. All these findings support a model by which the inhibition of Microprocessor activity is an essential step to induce a robust type I IFN response in mammalian cells.

Conserved structural RNA domains in regions coding for cleavage site motifs in hemagglutinin genes of influenza viruses

The acquisition of a multibasic cleavage site (MBCS) in the hemagglutinin (HA) glycoprotein is the main determinant of the conversion of low pathogenic avian influenza viruses into highly pathogenic strains, facilitating HA cleavage and virus replication in a broader range of host cells. In nature, substitutions or insertions in HA RNA genomic segments that code for multiple basic amino acids have been observed only in the HA genes of two out of sixteen subtypes circulating in birds, H5 and H7. Given the compatibility of MBCS motifs with HA proteins of numerous subtypes, this selectivity was hypothesized to be determined by the existence of specific motifs in HA RNA, in particular structured domains. In H5 and H7 HA RNAs, predictions of such domains have yielded alternative conserved stem-loop structures with the cleavage site codons in the hairpin loops. Here, potential RNA secondary structures were analyzed in the cleavage site regions of HA segments of influenza viruses of different types and subtypes. H5- and H7-like stem-loop structures were found in all known influenza A virus subtypes and in influenza B and C viruses with homology modeling. Nucleotide covariations supported this conservation to be determined by RNA structural constraints that are stronger in the domain-closing bottom stems as compared to apical parts. The structured character of this region in (sub-)types other than H5 and H7 indicates its functional importance beyond the ability to evolve toward an MBCS responsible for a highly pathogenic phenotype.

Harnessing host-virus evolution in antiviral therapy and immunotherapy

Pathogen resistance and development costs are major challenges in current approaches to antiviral therapy. The high error rate of RNA synthesis and reverse‐transcription confers genome plasticity, enabling the remarkable adaptability of RNA viruses to antiviral intervention. However, this property is coupled to fundamental constraints including limits on the size of information available to manipulate complex hosts into supporting viral replication. Accordingly, RNA viruses employ various means to extract maximum utility from their informationally limited genomes that, correspondingly, may be leveraged for effective host‐oriented therapies. Host‐oriented approaches are becoming increasingly feasible because of increased availability of bioactive compounds and recent advances in immunotherapy and precision medicine, particularly genome editing, targeted delivery methods and RNAi. In turn, one driving force behind these innovations is the increasingly detailed understanding of evolutionarily diverse host–virus interactions, which is the key concern of an emerging field, neo‐virology. This review examines biotechnological solutions to disease and other sustainability issues of our time that leverage the properties of RNA and DNA viruses as developed through co‐evolution with their hosts.

A novel method for the capture-based purification of whole viral native RNA genomes

Current technologies for targeted characterization and manipulation of viral RNA primarily involve amplification or ultracentrifugation with isopycnic gradients of viral particles to decrease host RNA background. The former strategy is non-compatible for characterizing properties innate to RNA strands such as secondary structure, RNA-RNA interactions, and also for nanopore direct RNA sequencing involving the sequencing of native RNA strands. The latter strategy, ultracentrifugation, causes loss in genomic information due to its inability to retrieve unassembled viral RNA. To address this, we developed a novel application of current nucleic acid hybridization technologies for direct characterization of RNA. In particular, we modified a current enrichment protocol to capture whole viral native RNA genomes for downstream RNA assays to circumvent the abovementioned problems. This technique involves hybridization of biotinylated baits at 500 nucleotides (nt) intervals, stringent washes and release of free native RNA strands using DNase I treatment, with a turnaround time of about 6 h 15 min. RT-qPCR was used as the primary proof of concept that capture-based purification indeed removes host background. Subsequently, capture-based purification was applied to direct RNA sequencing as proof of concept that capture-based purification can be coupled with downstream RNA assays. We report that this protocol was able to successfully purify viral RNA by 561- to 791-fold. We also report that application of this protocol to direct RNA sequencing yielded a reduction in human host RNA background by 1580-fold, a 99.91% recovery of viral genome with at least 15× coverage, and a mean coverage across the genome of 120×. This report is, to the best of our knowledge, the first description of a capture-based purification method for assays that involve direct manipulation or characterisation of native RNA. This report also describes a successful application of capture-based purification as a direct RNA sequencing strategy that addresses certain limitations of current strategies in sequencing RNA viral genomes.

Metallo supramolecular cylinders inhibit HIV-1 TAR-TAT complex formation and viral replication in cellulo

Shape-selective recognition of nucleic acid structures by supramolecular drugs offers the potential to treat disease. The Trans Activation Response (TAR) region is a region of high secondary structure within the human immunodeficiency virus-1 (HIV-1) RNA that complexes with the virus-encoded Transactivator protein (TAT) and regulates viral transcription. Herein, we explore different metallo-supramolecular triple stranded helicates (cylinders) that target the TAR bulge motif and inhibit the formation of TAR-TAT complexes and HIV infection. Cylinders that incorporate Ni(II) and Ru(II) showed the most potent anti-viral activity with limited evidence of cellular cytotoxicity. These metallo-supramolecular compounds provide an exciting avenue for developing a new class of anti-viral agents.

Small RNA-based interactions between rice and the viruses which cause the tungro disease

Rice tungro disease is caused by a complex of two viruses, Rice tungro bacilliform virus (RTBV) and Rice tungro spherical virus (RTSV). To examine the RNAi-based defence response in rice during tungro disease, we characterized the virus-derived small RNAs and miRNAs by Deep Sequencing. We found that, while 21 nt/22 nt (nucleotide) siRNAs are predominantly produced in a continuous, overlapping and asymmetrical manner from RTBV, siRNA accumulation from RTSV were negligible. Additionally, 54 previously known miRNAs from rice, predicted to be regulating genes involved in plant defence, hormone signaling and developmental pathways were differentially expressed in the infected samples, compared to the healthy ones. This is the first study of sRNA profile of tungro virus complex from infected rice plants. The biased response of the host antiviral machinery against the two viruses and the differentially-expressed miRNAs are novel observations, which entail further studies.

The 5BSL3.2 Functional RNA Domain Connects Distant Regions in the Hepatitis C Virus Genome

Viral genomes are complexly folded entities that carry all the information required for the infective cycle. The nucleotide sequence of the RNA virus genome encodes proteins and functional information contained in discrete, highly conserved structural units. These so-called functional RNA domains play essential roles in the progression of infection, which requires their preservation from one generation to the next. Numerous functional RNA domains exist in the genome of the hepatitis C virus (HCV). Among them, the 5BSL3.2 domain in the cis-acting replication element (CRE) at the 3' end of the viral open reading frame has become of particular interest given its role in HCV RNA replication and as a regulator of viral protein synthesis. These functionalities are achieved via the establishment of a complex network of long-distance RNA-RNA contacts involving (at least as known to date) the highly conserved 3'X tail, the apical loop of domain IIId in the internal ribosome entry site, and/or the so-called Alt region upstream of the CRE. Changing contacts promotes the execution of different stages of the viral cycle. The 5BSL3.2 domain thus operates at the core of a system that governs the progression of HCV infection. This review summarizes our knowledge of the long-range RNA-RNA interaction network in the HCV genome, with special attention paid to the structural and functional consequences derived from the establishment of different contacts. The potential implications of such interactions in switching between the different stages of the viral cycle are discussed.

Quasispecies evolution of the prototypical genotype 1 porcine reproductive and respiratory syndrome virus early during in vivo infection is rapid and tissue specific

Porcine reproductive and respiratory syndrome virus (PRRSV) is a major infectious threat to the pig industry worldwide. Increasing evidence suggests that microevolution within a quasispecies population can give rise to high sequence heterogeneity in PRRSV; potentially impacting the pathogenicity of the virus. Here, we report on micro-evolutionary events taking place within the viral quasispecies population in lung and lymph node 3 days post infection (dpi) following experimental in vivo infection with the prototypical Lelystad PRRSV (LV). Sequence analysis revealed 16 high frequency single nucleotide variants (SNV) or differences from the reference LV genome which are assumed to be representative of the consensus inoculum genome. Additionally, 49 other low frequency SNVs were also found in the inoculum population. At 3 dpi, a total of 9 and 10 SNVs of varying frequencies could already be detected in the LV population infecting the lung and lymph nodes, respectively. Interestingly, of these, three and four novel SNVs emerged independently in the two respective tissues when compared to the inoculum. The remaining variants, though already present at lower frequencies in the inoculum, were positively selected and their frequency increased within the quasispecies population. Hence, we were able to determine directly from tissues infected with PRRSV the repertoire of genetic variants within the viral quasispecies population. Our data also suggest that microevolution of these variants is rapid and some may be tissue-specific.

Synonymous Co-Variation across the E1/E2 Gene Junction of Hepatitis C Virus Defines Virion Fitness

Hepatitis C virus is a positive-sense single-stranded RNA virus. The gene junction partitioning the viral glycoproteins E1 and E2 displays concurrent sequence evolution with the 3'-end of E1 highly conserved and the 5'-end of E2 highly heterogeneous. This gene junction is also believed to contain structured RNA elements, with a growing body of evidence suggesting that such structures can act as an additional level of viral replication and transcriptional control. We have previously used ultradeep pyrosequencing to analyze an amplicon library spanning the E1/E2 gene junction from a treatment naïve patient where samples were collected over 10 years of chronic HCV infection. During this timeframe maintenance of an in-frame insertion, recombination and humoral immune targeting of discrete virus sub-populations was reported. In the current study, we present evidence of epistatic evolution across the E1/E2 gene junction and observe the development of co-varying networks of codons set against a background of a complex virome with periodic shifts in population dominance. Overtime, the number of codons actively mutating decreases for all virus groupings. We identify strong synonymous co-variation between codon sites in a group of sequences harbouring a 3 bp in-frame insertion and propose that synonymous mutation acts to stabilize the RNA structural backbone.

Sequence-Specific Sensing of Nucleic Acids

Innate immune cells are endowed with many nucleic acid receptors, but the role of sequence in the detection of foreign organisms remains unclear. Can sequence patterns influence recognition? In addition, how can we infer those patterns from sequence data? Here, we detail recent computational and experimental evidence associated with sequence-specific sensing. We review the mechanisms underlying the detection and discrimination of foreign sequences from self. We also describe quantitative approaches used to infer the stimulatory capacity of a given pathogen nucleic acid species, and the influence of sequence-specific sensing on host-pathogen coevolution, including endogenous sequences of foreign origin. Finally, we speculate how further studies of sequence-specific sensing will be useful to improve vaccine design, gene therapy and cancer treatment.

Movements of HIV-1 genomic RNA-APOBEC3F complexes and PKR reveal cytoplasmic and nuclear PKR defenses and HIV-1 evasion strategies

APOBEC3 cytidine deaminases and viral genomic RNA (gRNA) occur in virions, polysomes, and cytoplasmic granules, but have not been tracked together. Moreover, gRNA traffic is important, but the factors that move it into granules are unknown. Using in situ hybridization of transfected cells and protein synthesis inhibitors that drive mRNAs between locales, we observed APOBEC3F cotrafficking with gRNA without altering its movements. Whereas cells with little cytoplasmic gRNA were translationally active and accumulated Gag, suprathreshold amounts induced autophosphorylation of the cytoplasmic double-stranded RNA (dsRNA)-dependent protein kinase (PKR), causing eIF2α phosphorylation, protein synthesis suppression, and gRNA sequestration in stress granules. Additionally, we confirmed recent evidence that PKR is activated by chromosome-associated cellular dsRNAs after nuclear membranes disperse in prophase. By arresting cells in G2, HIV-1 blocks this mechanism for PKR activation and eIF2α phosphorylation. However, cytopathic membrane damage in CD4- and coreceptor-positive cultures infected with laboratory-adapted fusogenic HIV-1LAI eventually enabled PKR entry and activation in interphase nuclei. These results reveal multiple stages in the PKR-HIV-1 battleground that culminate in cell death. We discuss evidence suggesting that HIV-1s evolve in vivo to prevent or delay PKR activation by all these mechanisms.

Innate immune system activation by viral RNA: How to predict it?

The immune system is able to identify foreign pathogens via different pathways. In the case of viral infection, recognition of the viral RNA is a crucial step, and many efforts have been made to understand which features of viral RNA are detected by the immune system. The biased viral RNA composition, measured as host-virus nucleotidic divergence, or CpG enrichment, has been proposed as salient signal. Peculiar structural features of these RNA could also be related to the immune system activation. Here, we gather multiple datasets and proceed to a meta-analysis to uncover the best predictors of immune system activation by viral RNA. "A" nucleotide content and Minimum Folding Energy are good predictors, and are more easily generalized than more complex indicators suggested previously. As RNA composition and structure are highly correlated, we suggest further experiments on synthetic sequences to identify the viral RNA sensing mechanisms by immune system receptors.

Requirement for Host RNA-Silencing Components and the Virus-Silencing Suppressor when Second-Site Mutations Compensate for Structural Defects in the 3' Untranslated Region

Turnip crinkle virus (TCV) contains a structured 3' region with hairpins and pseudoknots that form a complex network of noncanonical RNA:RNA interactions supporting higher-order structure critical for translation and replication. We investigated several second-site mutations in the p38 coat protein open reading frame (ORF) that arose in response to a mutation in the asymmetric loop of a critical 3' untranslated region (UTR) hairpin that disrupts local higher-order structure. All tested second-site mutations improved accumulation of TCV in conjunction with a partial reversion of the primary mutation (TCV-rev1) but had neutral or a negative effect on wild-type (wt) TCV or TCV with the primary mutation. SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) structure probing indicated that these second-site mutations reside in an RNA domain that includes most of p38 (domain 2), and evidence for RNA:RNA interactions between domain 2 and 3'UTR-containing domain 1 was found. However, second-site mutations were not compensatory in the absence of p38, which is also the TCV silencing suppressor, or in dcl-2/dcl4 or ago1/ago2 backgrounds. One second-site mutation reduced silencing suppressor activity of p38 by altering one of two GW motifs that are required for p38 binding to double-stranded RNAs (dsRNAs) and interaction with RNA-induced silencing complex (RISC)-associated AGO1/AGO2. Another second-site mutation substantially reduced accumulation of TCV-rev1 in the absence of p38 or DCL2/DCL4. We suggest that the second-site mutations in the p38 ORF exert positive effects through a similar downstream mechanism, either by enhancing accumulation of beneficial DCL-produced viral small RNAs that positively regulate the accumulation of TCV-rev1 or by affecting the susceptibility of TCV-rev1 to RISC loaded with viral small RNAs.

IMPORTANCE

Genomes of positive-strand RNA viruses fold into high-order RNA structures. Viruses with mutations in regions critical for translation and replication often acquire second-site mutations that exert a positive compensatory effect through reestablishment of canonical base pairing with the altered region. In this study, two distal second-site mutations that individually arose in response to a primary mutation in a critical 3' UTR hairpin in the genomic RNA of turnip crinkle virus did not directly interact with the primary mutation. Although different second-site changes had different attributes, compensation was dependent on the production of the viral p38 silencing suppressor and on the presence of silencing-required DCL and AGO proteins. Our results provide an unexpected connection between a 3' UTR primary-site mutation proposed to disrupt higher-order structure and the RNA-silencing machinery.

Identification of a Natural Viral RNA Motif That Optimizes Sensing of Viral RNA by RIG-I

Stimulation of the antiviral response depends on the sensing of viral pathogen-associated molecular patterns (PAMPs) by specialized cellular proteins. During infection with RNA viruses, 5'-di- or -triphosphates accompanying specific single or double-stranded RNA motifs trigger signaling of intracellular RIG-I-like receptors (RLRs) and initiate the antiviral response. Although these molecular signatures are present during the replication of many viruses, it is unknown whether they are sufficient for strong activation of RLRs during infection. Immunostimulatory defective viral genomes (iDVGs) from Sendai virus (SeV) are among the most potent natural viral triggers of antiviral immunity. Here we describe an RNA motif (DVG(70-114)) that is essential for the potent immunostimulatory activity of 5'-triphosphate-containing SeV iDVGs. DVG(70-114) enhances viral sensing by the host cell independently of the long stretches of complementary RNA flanking the iDVGs, and it retains its stimulatory potential when transferred to otherwise inert viral RNA. In vitro analysis showed that DVG(70-114) augments the binding of RIG-I to viral RNA and promotes enhanced RIG-I polymerization, thereby facilitating the onset of the antiviral response. Together, our results define a new natural viral PAMP enhancer motif that promotes viral recognition by RLRs and confers potent immunostimulatory activity to viral RNA.

IMPORTANCE

A discrete group of molecular motifs, including 5'-triphosphates associated with double-stranded RNA, have been identified as essential for the triggering of antiviral immunity. Most RNA viruses expose these motifs during their replication; however, successful viruses normally evade immune recognition and replicate to high levels before detection, indicating that unknown factors drive antiviral immunity. DVGs from SeV are among the most potent natural viral stimuli of the antiviral response known to date. These studies define a new natural viral motif present in DVGs that maximizes viral recognition by the intracellular sensor RIG-I, allowing fast and strong antiviral responses even in the presence of viral-encoded immune antagonists. This motif can be harnessed to increase the immunostimulatory potential of otherwise inert viral RNAs and represents a novel immunostimulatory enhancer that could be used in the development of vaccine adjuvants and antivirals.

Inhibition of HCV translation by disrupting the structure and interactions of the viral CRE and 3' X-tail

A phylogenetically conserved RNA structure within the NS5B coding region of hepatitis C virus functions as a cis-replicating element (CRE). Integrity of this CRE, designated SL9266 (alternatively 5BSL3.2), is critical for genome replication. SL9266 forms the core of an extended pseudoknot, designated SL9266/PK, involving long distance RNA–RNA interactions between unpaired loops of SL9266 and distal regions of the genome. Previous studies demonstrated that SL9266/PK is dynamic, with ‘open’ and ‘closed’ conformations predicted to have distinct functions during virus replication. Using a combination of site-directed mutagenesis and locked nucleic acids (LNA) complementary to defined domains of SL9266 and its interacting regions, we have explored the influence of this structure on genome translation and replication. We demonstrate that LNAs which block formation of the closed conformation inhibit genome translation. Inhibition was at least partly independent of the initiation mechanism, whether driven by homologous or heterologous internal ribosome entry sites or from a capped message. Provision of SL9266/PK in trans relieved translational inhibition, and mutational analysis implied a mechanism in which the closed conformation recruits a cellular factor that would otherwise suppresses translation. We propose that SL9266/PK functions as a temporal switch, modulating the mutually incompatible processes of translation and replication.

Diverse roles and interactions of RNA structures during the replication of positive-stranded RNA viruses of humans and animals

Positive-stranded RNA viruses include important human, animal and plant pathogens. Their genomes are able to fold into complex structures stabilized by base pairing between individual nucleotides, many of which are highly conserved and have essential functions during virus replication. With new studies and technological advances the diversity of roles, mechanisms and interactions in which such structured viral RNA functions is becoming increasingly clear. It is also evident that many RNA structures do not function as discrete elements but through mechanisms involving multiple, long-range and often dynamic RNARNA interactions. Through a range of examples and recent advances, this review illustrates the diverse roles and mechanisms of structured viral RNA during the replication of positive-stranded RNA viruses infecting humans and animals.

Influenza A viruses with different amino acid residues at PB2-627 display distinct replication properties in vitro and in vivo: revealing the sequence plasticity of PB2-627 position

Sequence analyses of influenza PB2 sequences indicate that the 627 position almost exclusively contains either lysine (K) or glutamic acid (E), suggesting a high sequence constraint at this genetic marker. Here, we used a site-directed random mutagenesis method to demonstrate that PB2-627 position has a high sequence plasticity. Recombinant viruses carrying various amino acid residues at this position are viable in cell cultures. These PB2-627 mutants showed various polymerase activities and replication kinetics in mammalian and avian cells as well as pathogenicity in mice. Serially passaging these mutants in MDCK cells generated some compensatory PB2 mutations that can restore polymerase activities of the PB2-627 mutants. Of these, PB2-D309N was identified as a novel one. Besides showing that influenza virus can tolerate a wide range of amino acid residues at the PB2-627 position, this study also demonstrates a potential strategy to identify novel mutations that can enhance viral polymerase.

Sensing viral invasion by RIG-I like receptors

Cellular responses to pathogen invasion are crucial for maintaining cell homeostasis and survival. The interferon (IFN) system is one of the most effective cellular responses to viral intrusion in mammals. Viral recognition by innate immune sensors activates the antiviral IFN system. Retinoic acid-inducible gene I (RIG-I) like receptors (RLRs) are DExD/H box RNA helicases that sense viral invasion. RLRs recognize cytoplasmic viral RNAs and trigger antiviral responses, resulting in production of type I IFN and inflammatory cytokines. Unique and common sensing mechanisms among RLRs have been reported. In this review, recent progress in the understanding of antiviral responses by RLRs is summarized and discussed.