Bioinformatic and physical characterizations of genome-scale ordered RNA structure in mammalian RNA viruses

Determining the efficacy of disinfectants at nucleic acid degradation

AIMS

Nucleic acids, particularly antibiotic resistance genes, are commonly found on surfaces within healthcare environments, with levels not reducing following cleaning. Within the UK, there are no regulations for testing disinfectants against nucleic acids.

METHODS AND RESULTS

A series of commonplace in vitro methods were used to determine disinfectant-induced physical and functional damage to various nucleic acids; RNA (10 μg), genomic DNA (2 μg), and plasmids (1 μg). Using these methods, the optimal residence time (10 minutes) and working concentration (10%) were determined for a new disinfectant. Furthermore, comparison of disinfectants with different active ingredients including lactic acid (LA), sodium hydroxide (NaOH), chloroxylenol (PCMX), and quaternary ammonium compounds (QACs), were compared to controls. All disinfectants showed greater degradation by gel electrophoresis of genomic DNA and RNA than of purified plasmids. Functional analysis using quantitative polymerase chain reaction (qPCR) and polymerase chain reaction (PCR) demonstrated that no disinfectant tested (apart from control) could damage DNA to the level where PCR amplification was not possible, and only the NaOH reagent could achieve this for RNA.

CONCLUSIONS

The set of methods described herein provides a platform for future standardization and potential regulation regarding monitoring cleaning solutions for their activity against nucleic acids.

Let's phase it: viruses are master architects of biomolecular condensates

Viruses compartmentalize their replication and assembly machinery to both evade detection and concentrate the viral proteins and nucleic acids necessary for genome replication and virion production. Accumulating evidence suggests that diverse RNA and DNA viruses form replication organelles and nucleocapsid assembly sites using phase separation. In general, the biogenesis of these compartments is regulated by two types of viral protein, collectively known as antiterminators and nucleocapsid proteins, respectively. Herein, we discuss how RNA viruses establish replication organelles and nucleocapsid assembly sites, and the evidence that these compartments form through phase separation. While this review focuses on RNA viruses, accumulating evidence suggests that all viruses rely on phase separation and form biomolecular condensates important for completing the infectious cycle.

Generated Randomly and Selected Functionally? The Nature of Enterovirus Recombination

Genetic recombination in RNA viruses is an important evolutionary mechanism. It contributes to population diversity, host/tissue adaptation, and compromises vaccine efficacy. Both the molecular mechanism and initial products of recombination are relatively poorly understood. We used an established poliovirus-based in vitro recombination assay to investigate the roles of sequence identity and RNA structure, implicated or inferred from an analysis of circulating recombinant viruses, in the process. In addition, we used next-generation sequencing to investigate the early products of recombination after cellular coinfection with different poliovirus serotypes. In independent studies, we find no evidence for a role for RNA identity or structure in determining recombination junctions location. Instead, genome function and fitness are of greater importance in determining the identity of recombinant progeny. These studies provide further insights into this important evolutionary mechanism and emphasize the critical nature of the selection process on a mixed virus population.

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.

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.

A map of the SARS-CoV-2 RNA structurome

SARS-CoV-2 has exploded throughout the human population. To facilitate efforts to gain insights into SARS-CoV-2 biology and to target the virus therapeutically, it is essential to have a roadmap of likely functional regions embedded in its RNA genome. In this report, we used a bioinformatics approach, ScanFold, to deduce the local RNA structural landscape of the SARS-CoV-2 genome with the highest likelihood of being functional. We recapitulate previously-known elements of RNA structure and provide a model for the folding of an essential frameshift signal. Our results find that SARS-CoV-2 is greatly enriched in unusually stable and likely evolutionarily ordered RNA structure, which provides a large reservoir of potential drug targets for RNA-binding small molecules. Results are enhanced via the re-analyses of publicly-available genome-wide biochemical structure probing datasets that are broadly in agreement with our models. Additionally, ScanFold was updated to incorporate experimental data as constraints in the analysis to facilitate comparisons between ScanFold and other RNA modelling approaches. Ultimately, ScanFold was able to identify eight highly structured/conserved motifs in SARS-CoV-2 that agree with experimental data, without explicitly using these data. All results are made available via a public database (the RNAStructuromeDB: https://structurome.bb.iastate.edu/sars-cov-2) and model comparisons are readily viewable at https://structurome.bb.iastate.edu/sars-cov-2-global-model-comparisons.

Extensive C->U transition biases in the genomes of a wide range of mammalian RNA viruses; potential associations with transcriptional mutations, damage- or host-mediated editing of viral RNA

The rapid evolution of RNA viruses has been long considered to result from a combination of high copying error frequencies during RNA replication, short generation times and the consequent extensive fixation of neutral or adaptive changes over short periods. While both the identities and sites of mutations are typically modelled as being random, recent investigations of sequence diversity of SARS coronavirus 2 (SARS-CoV-2) have identified a preponderance of C->U transitions, proposed to be driven by an APOBEC-like RNA editing process. The current study investigated whether this phenomenon could be observed in datasets of other RNA viruses. Using a 5% divergence filter to infer directionality, 18 from 36 datasets of aligned coding region sequences from a diverse range of mammalian RNA viruses (including Picornaviridae, Flaviviridae, Matonaviridae, Caliciviridae and Coronaviridae) showed a >2-fold base composition normalised excess of C->U transitions compared to U->C (range 2.1x-7.5x), with a consistently observed favoured 5' U upstream context. The presence of genome scale RNA secondary structure (GORS) was the only other genomic or structural parameter significantly associated with C->U/U->C transition asymmetries by multivariable analysis (ANOVA), potentially reflecting RNA structure dependence of sites targeted for C->U mutations. Using the association index metric, C->U changes were specifically over-represented at phylogenetically uninformative sites, potentially paralleling extensive homoplasy of this transition reported in SARS-CoV-2. Although mechanisms remain to be functionally characterised, excess C->U substitutions accounted for 11-14% of standing sequence variability of structured viruses and may therefore represent a potent driver of their sequence diversification and longer-term evolution.

Whole genome analysis of Aichivirus A, isolated from a child, suffering from gastroenteritis, in Pakistan

Viruses are the primary cause of acute gastroenteritis in children all over the world. Understanding the emergence and genetic variation of these viruses may help to prevent infections. Aichivirus (AiV) is a member of the Kobuvirus genus, which currently contains six officially recognized species: Aichivirus A-F. The species AiV A contains six types including Aichivirus 1 (AiV 1) and eventually, three genotypes have been identified in the human AiV 1 (named A to C). The present study describes the identification and sequencing of the polyprotein gene of a human AiV 1 strain PAK419 via NGS in Pakistani children with acute gastroenteritis. Our study strain PAK419 was classified as AiV 1 genotype A, most commonly found in Japan and Europe, and closely related to non-Japanese and European strains on the phylogenetic tree. PAK419 showed 95-98 % nucleotide sequence identity with strains isolated from Ethiopia (ETH/2016/P4), Australia (FSS693) and China (Chshc7). On phylogenetic observation PAK419 formed a distinct cluster in the AiV 1 genotype A with the above mentioned and other human AiV strains detected around the world (Germany, Brazil, Japan, Thailand, Korea and Vietnam). The data clearly showed that Pakistani AiV strains and human strains identified from all over the world are distinct from Aichivirus strains found in bovine, swine, canine, feline, caprine, ferret, bat, and environmental samples. The distinguishing characteristics of the AiV genome showed a lower probability of inter-genotypic recombination events, which may support the lack of AiV serotypes. PAK419 also had a high content of C nucleotide (37.4 %), as found in previous studies, which could also restrict the possible genetic variation of AiV. This study demonstrate the power of NGS in uncovering unknown gastroenteric etiological agents circulating in the population.

Understanding the outcomes of COVID-19 - does the current model of an acute respiratory infection really fit?

Although coronavirus disease 2019 (COVID-19) is regarded as an acute, resolving infection followed by the development of protective immunity, recent systematic literature review documents evidence for often highly prolonged shedding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in respiratory and faecal samples, periodic recurrence of PCR positivity in a substantial proportion of individuals and increasingly documented instances of reinfection associated with a lack of protective immunity. This pattern of infection is quite distinct from the acute/resolving nature of other human pathogenic respiratory viruses, such as influenza A virus and respiratory syncytial virus. Prolonged shedding of SARS-CoV-2 furthermore occurs irrespective of disease severity or development of virus-neutralizing antibodies. SARS-CoV-2 possesses an intensely structured RNA genome, an attribute shared with other human and veterinary coronaviruses and with other mammalian RNA viruses such as hepatitis C virus. These are capable of long-term persistence, possibly through poorly understood RNA structure-mediated effects on innate and adaptive host immune responses. The assumption that resolution of COVID-19 and the appearance of anti-SARS-CoV-2 IgG antibodies represents virus clearance and protection from reinfection, implicit for example in the susceptible-infected-recovered (SIR) model used for epidemic prediction, should be rigorously re-evaluated.

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.

Impact of virus subtype and host IFNL4 genotype on large-scale RNA structure formation in the genome of hepatitis C virus

Mechanisms underlying the ability of hepatitis C virus (HCV) to establish persistent infections and induce progressive liver disease remain poorly understood. HCV is one of several positive-stranded RNA viruses capable of establishing persistence in their immunocompetent vertebrate hosts, an attribute previously associated with formation of large-scale RNA structure in their genomic RNA. We developed novel methods to analyze and visualize genome-scale ordered RNA structure (GORS) predicted from the increasingly large data sets of complete genome sequences of HCV. Structurally conserved RNA secondary structure in coding regions of HCV localized exclusively to polyprotein ends (core, NS5B). Coding regions elsewhere were also intensely structured based on elevated minimum folding energy difference (MFED) values, but the actual stem-loop elements involved in genome folding were structurally poorly conserved, even between subtypes 1a and 1b. Dynamic remodeling was further evident from comparison of HCV strains in different host genetic backgrounds. Significantly higher MFED values, greater suppression of UpA dinucleotide frequencies, and restricted diversification were found in subjects with the TT genotype of the rs12979860 SNP in the IFNL4 gene compared to the CC (nonexpressing) allele. These structural and compositional associations with expression of interferon-λ4 were recapitulated on a larger scale by higher MFED values and greater UpA suppression of genotype 1 compared to genotype 3a, associated with previously reported HCV genotype-associated differences in hepatic interferon-stimulated gene induction. Associations between innate cellular responses with HCV structure and further evolutionary constraints represent an important new element in RNA virus evolution and the adaptive interplay between virus and host.

Pervasive RNA Secondary Structure in the Genomes of SARS-CoV-2 and Other Coronaviruses

The ultimate outcome of the coronavirus disease 2019 (COVID-19) pandemic is unknown and is dependent on a complex interplay of its pathogenicity, transmissibility, and population immunity. In the current study, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was investigated for the presence of large-scale internal RNA base pairing in its genome. This property, termed genome-scale ordered RNA structure (GORS) has been previously associated with host persistence in other positive-strand RNA viruses, potentially through its shielding effect on viral RNA recognition in the cell. Genomes of SARS-CoV-2 were remarkably structured, with minimum folding energy differences (MFEDs) of 15%, substantially greater than previously examined viruses such as hepatitis C virus (HCV) (MFED of 7 to 9%). High MFED values were shared with all coronavirus genomes analyzed and created by several hundred consecutive energetically favored stem-loops throughout the genome. In contrast to replication-associated RNA structure, GORS was poorly conserved in the positions and identities of base pairing with other sarbecoviruses-even similarly positioned stem-loops in SARS-CoV-2 and SARS-CoV rarely shared homologous pairings, indicative of more rapid evolutionary change in RNA structure than in the underlying coding sequences. Sites predicted to be base paired in SARS-CoV-2 showed less sequence diversity than unpaired sites, suggesting that disruption of RNA structure by mutation imposes a fitness cost on the virus that is potentially restrictive to its longer evolution. Although functionally uncharacterized, GORS in SARS-CoV-2 and other coronaviruses represents important elements in their cellular interactions that may contribute to their persistence and transmissibility.IMPORTANCE The detection and characterization of large-scale RNA secondary structure in the genome of SARS-CoV-2 indicate an extraordinary and unsuspected degree of genome structural organization; this could be effectively visualized through a newly developed contour plotting method that displays positions, structural features, and conservation of RNA secondary structure between related viruses. Such RNA structure imposes a substantial evolutionary cost; paired sites showed greater restriction in diversity and represent a substantial additional constraint in reconstructing its molecular epidemiology. Its biological relevance arises from previously documented associations between possession of structured genomes and persistence, as documented for HCV and several other RNA viruses infecting humans and mammals. Shared properties potentially conferred by large-scale structure in SARS-CoV-2 include increasing evidence for prolonged infections and induced immune dysfunction that prevents development of protective immunity. The findings provide an additional element to cellular interactions that potentially influences the natural history of SARS-CoV-2, its pathogenicity, and its transmission.

The First Nonmammalian Pegivirus Demonstrates Efficient In Vitro Replication and High Lymphotropism

Members of the Pegivirus genus, family Flaviviridae, widely infect humans and other mammals, including nonhuman primates, bats, horses, pigs, and rodents, but are not associated with disease. Here, we report a new, genetically distinct pegivirus in goose (Anser cygnoides), the first identified in a nonmammalian host species. Goose pegivirus (GPgV) can be propagated in goslings, embryonated goose eggs, and primary goose embryo fibroblasts, and is thus the first pegivirus that can be efficiently cultured in vitro Experimental infection of GPgV in goslings via intravenous injection revealed robust replication and high lymphotropism. Analysis of the tissue tropism of GPgV revealed that the spleen and thymus were the organs bearing the highest viral loads. Importantly, GPgV could promote clinical manifestations of goose parvovirus infection, including reduced weight gain and 7% mortality. This finding contrasts with the lack of pathogenicity that is characteristic of previously reported pegiviruses.IMPORTANCE Members of the Pegivirus genus, family Flaviviridae, widely infect humans and other mammals, but are described as causing persistent infection and lacking pathogenicity. The efficiency of in vitro replication systems for pegivirus is poor, thus limiting investigation into viral replication steps. Because of that, the pathogenesis, cellular tropism, route of transmission, biology, and epidemiology of pegiviruses remain largely uncovered. Here, we report a phylogenetically distinct goose pegivirus (GPgV) that should be classified as a new species. GPgV proliferated in cell culture in a species- and cell-type-specific manner. Animal experiments show GPgV lymphotropism and promote goose parvovirus clinical manifestations. This study provides the first cell culture model for pegivirus, opening new possibilities for studies of pegivirus molecular biology. More importantly, our findings stand in contrast to the lack of identified pathogenicity of previously reported pegiviruses, which sheds lights on the pathobiology of pegivirus.

An in silico map of the SARS-CoV-2 RNA Structurome

SARS-CoV-2 is a positive-sense single-stranded RNA virus that has exploded throughout the global human population. This pandemic coronavirus strain has taken scientists and public health researchers by surprise and knowledge of its basic biology (e.g. structure/function relationships in its genomic, messenger and template RNAs) and modes for therapeutic intervention lag behind that of other human pathogens. In this report we used a recently-developed bioinformatics approach, ScanFold, to deduce the RNA structural landscape of the SARS-CoV-2 transcriptome. We recapitulate known elements of RNA structure and provide a model for the folding of an essential frameshift signal. Our results find that the SARS-CoV-2 is greatly enriched in unusually stable and likely evolutionarily ordered RNA structure, which provides a huge reservoir of potential drug targets for RNA-binding small molecules. Our results also predict regions that are accessible for intermolecular interactions, which can aid in the design of antisense therapeutics. All results are made available via a public database (the RNAStructuromeDB) where they may hopefully drive drug discovery efforts to inhibit SARS-CoV-2 pathogenesis.

Translation of the intrinsically disordered protein α-synuclein is inhibited by a small molecule targeting its structured mRNA

Many proteins are refractory to targeting because they lack small-molecule binding pockets. An alternative to drugging these proteins directly is to target the messenger (m)RNA that encodes them, thereby reducing protein levels. We describe such an approach for the difficult-to-target protein α-synuclein encoded by the SNCA gene. Multiplication of the SNCA gene locus causes dominantly inherited Parkinson's disease (PD), and α-synuclein protein aggregates in Lewy bodies and Lewy neurites in sporadic PD. Thus, reducing the expression of α-synuclein protein is expected to have therapeutic value. Fortuitously, the SNCA mRNA has a structured iron-responsive element (IRE) in its 5' untranslated region (5' UTR) that controls its translation. Using sequence-based design, we discovered small molecules that target the IRE structure and inhibit SNCA translation in cells, the most potent of which is named Synucleozid. Both in vitro and cellular profiling studies showed Synucleozid directly targets the α-synuclein mRNA 5' UTR at the designed site. Mechanistic studies revealed that Synucleozid reduces α-synuclein protein levels by decreasing the amount of SNCA mRNA loaded into polysomes, mechanistically providing a cytoprotective effect in cells. Proteome- and transcriptome-wide studies showed that the compound's selectivity makes Synucleozid suitable for further development. Importantly, transcriptome-wide analysis of mRNAs that encode intrinsically disordered proteins revealed that each has structured regions that could be targeted with small molecules. These findings demonstrate the potential for targeting undruggable proteins at the level of their coding mRNAs. This approach, as applied to SNCA, is a promising disease-modifying therapeutic strategy for PD and other α-synucleinopathies.

The role of ZAP and OAS3/RNAseL pathways in the attenuation of an RNA virus with elevated frequencies of CpG and UpA dinucleotides

Zinc finger antiviral protein (ZAP) is a powerful restriction factor for viruses with elevated CpG dinucleotide frequencies. We report that ZAP similarly mediates antiviral restriction against echovirus 7 (E7) mutants with elevated frequencies of UpA dinucleotides. Attenuation of both CpG- and UpA-high viruses and replicon mutants was reversed in ZAP k/o cell lines, and restored by plasmid-derived reconstitution of expression in k/o cells. In pull-down assays, ZAP bound to viral RNA transcripts with either CpG- and UpA-high sequences inserted in the R2 region. We found no evidence that attenuation of CpG- or UpA-high mutants was mediated through either translation inhibition or accelerated RNA degradation. Reversal of the attenuation of CpG-high, and UpA-high E7 viruses and replicons was also achieved through knockout of RNAseL and oligodenylate synthetase 3 (OAS3), but not OAS1. WT levels of replication of CpG- and UpA-high mutants were observed in OAS3 k/o cells despite abundant expression of ZAP, indicative of synergy or complementation of these hitherto unconnected pathways. The dependence on expression of ZAP, OAS3 and RNAseL for CpG/UpA-mediated attenuation and the variable and often low level expression of these pathway proteins in certain cell types, such as those of the central nervous system, has implications for the use of CpG-elevated mutants as attenuated live vaccines against neurotropic viruses.

Breast Cancer Resistance Protein (BCRP/ABCG2) Inhibits Extra Villous Trophoblast Migration: The Impact of Bacterial and Viral Infection

Extravillous trophoblasts (EVT) migration into the decidua is critical for establishing placental perfusion and when dysregulated, may lead to pre-eclampsia (PE) and intrauterine growth restriction (IUGR). The breast cancer resistance protein (BCRP; encoded by ABCG2) regulates the fusion of cytotrophoblasts into syncytiotrophoblasts and protects the fetus from maternally derived xenobiotics. Information about BCRP function in EVTs is limited, however placental exposure to bacterial/viral infection leads to BCRP downregulation in syncitiotrophoblasts. We hypothesized that BCRP is involved in the regulation of EVT function and is modulated by infection/inflammation. We report that besides syncitiotrophoblasts and cytotrophoblasts, BCRP is also expressed in EVTs. BCRP inhibits EVT cell migration in HTR8/SVneo (human EVT-like) cells and in human EVT explant cultures, while not affecting cell proliferation. We have also shown that bacterial-lipopolysaccharide (LPS)-and viral antigens-single stranded RNA (ssRNA)-have a profound effect in downregulating ABCG2 and BCRP levels, whilst simultaneously increasing the migration potential of EVT-like cells. Our study reports a novel function of BCRP in early placentation and suggests that exposure of EVTs to maternal infection/inflammation could disrupt their migration potential via the downregulation of BCRP. This could negatively influence placental development/function, contribute to existing obstetric pathologies, and negatively impact pregnancy outcomes and maternal/neonatal health.

Physical and Functional Analysis of Viral RNA Genomes by SHAPE

RNA viruses encode the information required to usurp cellular metabolism and gene regulation and to enable their own replication in two ways: in the linear sequence of their RNA genomes and in higher-order structures that form when the genomic RNA strand folds back on itself. Application of high-resolution SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) structure probing to viral RNA genomes has identified numerous new regulatory elements, defined new principles by which viral RNAs interact with the cellular host and evade host immune responses, and revealed relationships between virus evolution and RNA structure. This review summarizes our current understanding of genome structure-function interrelationships for RNA viruses, as informed by SHAPE structure probing, and outlines opportunities for future studies.

The Effect of RNA Substitution Models on Viroid and RNA Virus Phylogenies

Many viroids and RNA viruses have genomes that exhibit secondary structure, with paired nucleotides forming stems and loops. Such structures violate a key assumption of most methods of phylogenetic reconstruction, that sequence change is independent among sites. However, phylogenetic analyses of these transmissible agents rarely use evolutionary models that account for RNA secondary structure. Here, we assess the effect of using RNA-specific nucleotide substitution models on the phylogenetic inference of viroids and RNA viruses. We obtained data sets comprising full-genome nucleotide sequences from six viroid and ten single-stranded RNA virus species. For each alignment, we inferred consensus RNA secondary structures, then evaluated different DNA and RNA substitution models. We used model selection to choose the best-fitting model and evaluate estimated Bayesian phylogenies. Further, for each data set we generated and compared Robinson–Foulds (RF) statistics in order to test whether the distributions of trees generated under alternative models are notably different to each other. In all alignments, the best-fitting model was one that considers RNA secondary structure: RNA models that allow a nonzero rate of double substitution (RNA16A and RNA16C) fitted best for both viral and viroid data sets. In 14 of 16 data sets, the use of an RNA-specific model led to significantly longer tree lengths, but only in three cases did it have a significant effect on RFs. In conclusion, using RNA model when undertaking phylogenetic inference of viroids and RNA viruses can provide a better model fit than standard approaches and model choice can significantly affect branch length estimates.

Potential of the Other Genetic Information Coded by the Viral RNA Genomes as Antiviral Target

In addition to the protein coding information, viral RNA genomes code functional information in structurally conserved units termed functional RNA domains. These RNA domains play essential roles in the viral cycle (e.g., replication and translation). Understanding the molecular mechanisms behind their function is essential to understanding the viral infective cycle. Further, interfering with the function of the genomic RNA domains offers a potential means of developing antiviral strategies. Aptamers are good candidates for targeting structural RNA domains. Besides its potential as therapeutics, aptamers also provide an excellent tool for investigating the functionality of RNA domains in viral genomes. This review briefly summarizes the work carried out in our laboratory aimed at the structural and functional characterization of the hepatitis C virus (HCV) genomic RNA domains. It also describes the efforts we carried out for the development of antiviral aptamers targeting specific genomic domains of the HCV and the human immunodeficiency virus type-1 (HIV-1).

Genomic-Scale Interaction Involving Complementary Sequences in the Hepatitis C Virus 5'UTR Domain IIa and the RNA-Dependent RNA Polymerase Coding Region Promotes Efficient Virus Replication

The hepatitis C virus (HCV) genome contains structured elements thought to play important regulatory roles in viral RNA translation and replication processes. We used in vitro RNA binding assays to map interactions involving the HCV 5′UTR and distal sequences in NS5B to examine their impact on viral RNA replication. The data revealed that 5′UTR nucleotides (nt) 95–110 in the internal ribosome entry site (IRES) domain IIa and matching nt sequence 8528–8543 located in the RNA-dependent RNA polymerase coding region NS5B, form a high-affinity RNA-RNA complex in vitro. This duplex is composed of both wobble and Watson-Crick base-pairings, with the latter shown to be essential to the formation of the high-affinity duplex. HCV genomic RNA constructs containing mutations in domain IIa nt 95–110 or within the genomic RNA location comprising nt 8528–8543 displayed, on average, 5-fold less intracellular HCV RNA and 6-fold less infectious progeny virus. HCV genomic constructs containing complementary mutations for IRES domain IIa nt 95–110 and NS5B nt 8528–8543 restored intracellular HCV RNA and progeny virus titers to levels obtained for parental virus RNA. We conclude that this long-range duplex interaction between the IRES domain IIa and NS5B nt 8528–8543 is essential for optimal virus replication.

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.

Inhibition of dengue virus infection by small interfering RNAs that target highly conserved sequences in the NS4B or NS5 coding regions

Dengue fever is one of the most common viral infections in the world. Although a vaccine against dengue virus (DENV) has been approved in several countries, this disease is still considered a public health priority worldwide. The ability of three small interfering RNAs (FG-siRNAs) targeting conserved sequences in the NS4B and NS5 regions of the DENV genome to inhibit DENV replication was tested in vitro in both Vero and C6/36 cells. The FG-siRNAs were effective against DENV-1, -3, and -4, but not DENV-2. A fourth siRNA specifically targeting the NS5 region of the DENV-2 genome (SG-siRNA) was designed and tested against two different DENV-2 strains, showing high levels of inhibition in both mammalian and insect cells.

An ultraprocessive, accurate reverse transcriptase encoded by a metazoan group II intron

Group II introns and non-LTR retrotransposons encode a phylogenetically related family of highly processive reverse transcriptases (RTs) that are essential for mobility and persistence of these retroelements. Recent crystallographic studies on members of this RT family have revealed that they are structurally distinct from the retroviral RTs that are typically used in biotechnology. However, quantitative, structure-guided analysis of processivity, efficiency, and accuracy of this alternate RT family has been lacking. Here, we characterize the processivity of a group II intron maturase RT from Eubacterium rectale (E.r), for which high-resolution structural information is available. We find that the E.r. maturase RT (MarathonRT) efficiently copies transcripts at least 10 kb in length and displays superior intrinsic RT processivity compared to commercial enzymes such as Superscript IV (SSIV). The elevated processivity of MarathonRT is at least partly mediated by a loop structure in the finger subdomain that acts as a steric guard (the α-loop). Additionally, we find that a positively charged secondary RNA binding site on the surface of the RT diminishes the primer utilization efficiency of the enzyme, and that reengineering of this surface enhances capabilities of the MarathonRT. Finally, using single-molecule sequencing, we show that the error frequency of MarathonRT is comparable to that of other high-performance RTs, such as SSIV, which were tested in parallel. Our results provide a structural framework for understanding the enhanced processivity of retroelement RTs, and they demonstrate the potential for engineering a powerful new generation of RT tools for application in biotechnology and research.

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.

The chaperone-like activity of the hepatitis C virus IRES and CRE elements regulates genome dimerization

The RNA genome of the hepatitis C virus (HCV) establishes a network of long-distance RNA-RNA interactions that direct the progression of the infective cycle. This work shows that the dimerization of the viral genome, which is initiated at the dimer linkage sequence (DLS) within the 3′UTR, is promoted by the CRE region, while the IRES is a negative regulatory partner. Using differential 2′-acylation probing (SHAPE-dif) and molecular interference (HMX) technologies, the CRE activity was found to mainly lie in the critical 5BSL3.2 domain, while the IRES-mediated effect is dependent upon conserved residues within the essential structural elements JIIIabc, JIIIef and PK2. These findings support the idea that, along with the DLS motif, the IRES and CRE are needed to control HCV genome dimerization. They also provide evidences of a novel function for these elements as chaperone-like partners that fine-tune the architecture of distant RNA domains within the HCV genome.

Evaluation of Viremia Frequencies of a Novel Human Pegivirus by Using Bioinformatic Screening and PCR

Bioinformatic screening and PCR-based approaches detected active infection with human hepegivirus-1 in exposed populations.

Next-generation sequencing has critical applications in virus discovery, diagnostics, and environmental surveillance. We used metagenomic sequence libraries for retrospective screening of plasma samples for the recently discovered human hepegivirus 1 (HHpgV-1). From a cohort of 150 hepatitis C virus (HCV)–positive case-patients, we identified 2 persons with HHpgV-1 viremia and a high frequency of human pegivirus (HPgV) viremia (14%). Detection of HHpgV-1 and HPgV was concordant with parallel PCR-based screening using conserved primers matching groups 1 (HPgV) and 2 (HHPgV-1) nonstructural 3 region sequences. PCR identified 1 HHPgV-1–positive person with viremia from a group of 195 persons with hemophilia who had been exposed to nonvirally inactivated factor VII/IX; 18 (9%) were HPgV-positive. Relative to HCV and HPgV, active infections with HHpgV-1 were infrequently detected in blood, even in groups that had substantial parenteral exposure. Our findings are consistent with lower transmissibility or higher rates of virus clearance for HHpgV-1 than for other bloodborne human flaviviruses.

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.

Impact of the biased nucleotide composition of viral RNA genomes on RNA structure and codon usage

We are interested in the influence of nucleotide composition on the fundamental characteristics of the virus RNA genome. Most RNA viruses have genomes with a distinct nucleotide composition, e.g. ranging from minimally 12.9 % to maximally 40.3 % (C- and U-count, respectively, in coronavirus HKU). We present a global analysis of diverse virus types, including plus-strand, minus-strand and double-strand RNA viruses, for the impact of this nucleotide preference on the predicted structure of the RNA genome that is packaged in virion particles and on the codon usage in the viral open reading frames. Several virus-specific features will be described, but also some general conclusions were drawn. Without exception, the virus-specific nucleotide bias was enriched in the unpaired, single-stranded regions of the RNA genome, thus creating an even more striking virus-specific signature. We present a simple mechanism that is based on elementary aspects of RNA structure folding to explain this general trend. In general, the nucleotide bias was the major determinant of the virus-specific codon usages, thus limiting a role for codon selection and translational control. We will discuss molecular and evolutionary scenarios that may be responsible for the diverse nucleotide biases of RNA viruses.

The Coding Region of the HCV Genome Contains a Network of Regulatory RNA Structures

RNA is a versatile macromolecule that accommodates functional information in primary sequence and secondary and tertiary structure. We use a combination of chemical probing, RNA structure modeling, comparative sequence analysis, and functional assays to examine the role of RNA structure in the hepatitis C virus (HCV) genome. We describe a set of conserved but functionally diverse structural RNA motifs that occur in multiple coding regions of the HCV genome, and we demonstrate that conformational changes in these motifs influence specific stages in the virus' life cycle. Our study shows that these types of structures can pervade a genome, where they play specific mechanistic and regulatory roles, constituting a "code within the code" for controlling biological processes.

Low-Fidelity Polymerases of Alphaviruses Recombine at Higher Rates To Overproduce Defective Interfering Particles

Low-fidelity RNA-dependent RNA polymerases for many RNA virus mutators have been shown to confer attenuated phenotypes, presumably due to increased mutation rates. Additionally, for many RNA viruses, replication to high titers results in the production of defective interfering particles (DIs) that also attenuate infection. We hypothesized that fidelity, recombination, and DI production are tightly linked. We show that a Sindbis virus mutator replicating at a high multiplicity of infection manifests an earlier and greater accumulation of DIs than its wild-type counterpart. The isolated DIs interfere with the replication of full-length virus in a dose-dependent manner. Importantly, the ability of the mutator virus to overproduce DIs could be linked to an increased recombination frequency. These data confirm that RNA-dependent RNA polymerase fidelity and recombination are inversely correlated for this mutator. Our findings suggest that defective interference resulting from higher recombination rates may be more detrimental to RNA virus mutators than the increase in mutational burden.

IMPORTANCE Replication, adaptation, and evolution of RNA viruses rely in large part on their low-fidelity RNA-dependent RNA polymerase. Viruses artificially modified in their polymerases to decrease fidelity (mutator viruses) are attenuated in vivo, demonstrating the important role of fidelity in viral fitness. However, attenuation was attributed solely to the modification of the viral mutation rate and the accumulation of detrimental point mutations. In this work, we described an additional phenotype of mutator viruses: an increased recombination rate leading to defective interfering particle (DI) overproduction. Because DIs are known for their inhibitory effect on viral replication, our work suggests that fidelity variants may be attenuated in vivo via several mechanisms. This has important implications in the development of fidelity variants as live attenuated vaccine strains.

Large-scale screening and characterization of enteroviruses and kobuviruses infecting pigs in Vietnam

A recent survey of pigs in Dong Thap province, Vietnam identified a high frequency of enterovirus species G (EV-G) infection (144/198; 72.7%). Amongst these was a plethora of EV-G types (EV-G1, EV-G6 and four new types EV-G8-EV-G11). To better characterize the genetic diversity of EV-G and investigate the possible existence of further circulating types, we performed a larger-scale study on 484 pig and 45 farm-bred boar faecal samples collected in 2012 and 2014, respectively. All samples from the previous and current studies were also screened for kobuviruses. The overall EV infection frequency remained extremely high (395/484; 81.6%), but with comparable detection rates and viral loads between healthy and diarrhoeic pigs; this contrasted with less frequent detection of EV-G in boars (4/45; 8.9%). EV was most frequently detected in pigs ≤ 14 weeks old (∼ 95%) and declined in older pigs. Infections with EV-G1 and EV-G6 were most frequent, whilst less commonly detected types included EV-G3, EV-G4 and EV-G8-EV-G11, and five new types (EV-G12-EV-G16). In contrast, kobuvirus infection frequency was significantly higher in diarrhoeic pigs (40.9 versus 27.6%; P = 0.01). Kobuviruses also showed contrasting epizootiologies and age associations; a higher prevalence was found in boars (42%) compared with domestic pigs (29%), with the highest infection frequency amongst pigs >52 weeks old. Although genetically diverse, all kobuviruses identified belonged to the species Aichivirus C. In summary, this study confirms infection with EV-G was endemic in Vietnamese domestic pigs and exhibits high genetic diversity and extensive inter-type recombination.

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.

Full-length RNA structure prediction of the HIV-1 genome reveals a conserved core domain

A distance constrained secondary structural model of the ≈10 kb RNA genome of the HIV-1 has been predicted but higher-order structures, involving long distance interactions, are currently unknown. We present the first global RNA secondary structure model for the HIV-1 genome, which integrates both comparative structure analysis and information from experimental data in a full-length prediction without distance constraints. Besides recovering known structural elements, we predict several novel structural elements that are conserved in HIV-1 evolution. Our results also indicate that the structure of the HIV-1 genome is highly variable in most regions, with a limited number of stable and conserved RNA secondary structures. Most interesting, a set of long distance interactions form a core organizing structure (COS) that organize the genome into three major structural domains. Despite overlapping protein-coding regions the COS is supported by a particular high frequency of compensatory base changes, suggesting functional importance for this element. This new structural element potentially organizes the whole genome into three major domains protruding from a conserved core structure with potential roles in replication and evolution for the virus.

Virome Analysis of Transfusion Recipients Reveals a Novel Human Virus That Shares Genomic Features with Hepaciviruses and Pegiviruses

To investigate the transmission of novel infectious agents by blood transfusion, we studied changes in the virome composition of blood transfusion recipients pre- and posttransfusion. Using this approach, we detected and genetically characterized a novel human virus, human hepegivirus 1 (HHpgV-1), that shares features with hepatitis C virus (HCV) and human pegivirus (HPgV; formerly called GB virus C or hepatitis G virus). HCV and HPgV belong to the genera Hepacivirus and Pegivirus of the family Flaviviridae. HHpgV-1 was found in serum samples from two blood transfusion recipients and two hemophilia patients who had received plasma-derived clotting factor concentrates. In the former, the virus was detected only in the posttransfusion samples, indicating blood-borne transmission. Both hemophiliacs were persistently viremic over periods of at least 201 and 1,981 days. The 5′ untranslated region (UTR) of HHpgV-1 contained a type IV internal ribosome entry site (IRES), structurally similar to although highly divergent in sequence from that of HCV and other hepaciviruses. However, phylogenetic analysis of nonstructural genes (NS3 and NS5B) showed that HHpgV-1 forms a branch within the pegivirus clade distinct from HPgV and homologs infecting other mammalian species. In common with some pegivirus variants infecting rodents and bats, the HHpgV-1 genome encodes a short, highly basic protein upstream of E1, potentially possessing a core-like function in packaging RNA during assembly. Identification of this new human virus, HHpgV-1, expands our knowledge of the range of genome configurations of these viruses and may lead to a reevaluation of the original criteria by which the genera Hepacivirus and Pegivirus are defined.

More than 30 million blood components are transfused annually in the United States alone. Surveillance for infectious agents in the blood supply is key to ensuring the safety of this critical resource for medicine and public health. Here, we report the identification of a new and highly diverse HCV/GB virus (GBV)-like virus from human serum samples. This new virus, human hepegivirus 1 (HHpgV-1), was found in serum samples from blood transfusion recipients, indicating its potential for transmission via transfusion products. We also found persistent long-term HHpgV-1 viremia in two hemophilia patients. HHpgV-1 is unique because it shares genetic similarity with both highly pathogenic HCV and the apparently nonpathogenic HPgV (GBV-C). Our results add to the list of human viruses and provide data to develop reagents to study virus transmission and disease association and for interrupting virus transmission and new human infections.

Synonymous mutations reduce genome compactness in icosahedral ssRNA viruses

Recent studies have shown that single-stranded (ss) viral RNAs fold into more compact structures than random RNA sequences with similar chemical composition and identical length. Based on this comparison, it has been suggested that wild-type viral RNA may have evolved to be atypically compact so as to aid its encapsidation and assist the viral assembly process. To further explore the compactness selection hypothesis, we systematically compare the predicted sizes of >100 wild-type viral sequences with those of their mutants, which are evolved in silico and subject to a number of known evolutionary constraints. In particular, we enforce mutation synonynimity, preserve the codon-bias, and leave untranslated regions intact. It is found that progressive accumulation of these restricted mutations still suffices to completely erase the characteristic compactness imprint of the viral RNA genomes, making them in this respect physically indistinguishable from randomly shuffled RNAs. This shows that maintaining the physical compactness of the genome is indeed a primary factor among ssRNA viruses' evolutionary constraints, contributing also to the evidence that synonymous mutations in viral ssRNA genomes are not strictly neutral.

In silico discovery and modeling of non-coding RNA structure in viruses

This review covers several computational methods for discovering structured non-coding RNAs in viruses and modeling their putative secondary structures. Here we will use examples from two target viruses to highlight these approaches: influenza A virus-a relatively small, segmented RNA virus; and Epstein-Barr virus-a relatively large DNA virus with a complex transcriptome. Each system has unique challenges to overcome and unique characteristics to exploit. From these particular cases, generically useful approaches can be derived for the study of additional viral targets.

Visualizing the global secondary structure of a viral RNA genome with cryo-electron microscopy

The lifecycle, and therefore the virulence, of single-stranded (ss)-RNA viruses is regulated not only by their particular protein gene products, but also by the secondary and tertiary structure of their genomes. The secondary structure of the entire genomic RNA of satellite tobacco mosaic virus (STMV) was recently determined by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE). The SHAPE analysis suggested a single highly extended secondary structure with much less branching than occurs in the ensemble of structures predicted by purely thermodynamic algorithms. Here we examine the solution-equilibrated STMV genome by direct visualization with cryo-electron microscopy (cryo-EM), using an RNA of similar length transcribed from the yeast genome as a control. The cryo-EM data reveal an ensemble of branching patterns that are collectively consistent with the SHAPE-derived secondary structure model. Thus, our results both elucidate the statistical nature of the secondary structure of large ss-RNAs and give visual support for modern RNA structure determination methods. Additionally, this work introduces cryo-EM as a means to distinguish between competing secondary structure models if the models differ significantly in terms of the number and/or length of branches. Furthermore, with the latest advances in cryo-EM technology, we suggest the possibility of developing methods that incorporate restraints from cryo-EM into the next generation of algorithms for the determination of RNA secondary and tertiary structures.

Exploring the architecture of viral RNA genomes

The genomes of RNA viruses contain local structural elements and long-range interactions that control various steps in virus replication. While many individual RNA elements have been characterized, it remains less clear how the structure and activity of such elements are integrated and regulated within the complex context of complete viral genomes. Recent technical advances, particularly the development of high-throughput solution structure mapping methods, have made secondary structural analysis of entire viral RNA genomes feasible. As a consequence, whole-genome structural models have been deduced for a number of plus-strand RNA viruses and retroviruses and these structures have provided intriguing functional and evolutionary insights into global genome architecture.

Functionally conserved architecture of hepatitis C virus RNA genomes

Hepatitis C virus (HCV) infects over 170 million people worldwide and is a leading cause of liver disease and cancer. The virus has a 9,650-nt, single-stranded, messenger-sense RNA genome that is infectious as an independent entity. The RNA genome has evolved in response to complex selection pressures, including the need to maintain structures that facilitate replication and to avoid clearance by cell-intrinsic immune processes. Here we used high-throughput, single-nucleotide resolution information to generate and functionally test data-driven structural models for three diverse HCV RNA genomes. We identified, de novo, multiple regions of conserved RNA structure, including all previously characterized cis-acting regulatory elements and also multiple novel structures required for optimal viral fitness. Well-defined RNA structures in the central regions of HCV genomes appear to facilitate persistent infection by masking the genome from RNase L and double-stranded RNA-induced innate immune sensors. This work shows how structure-first comparative analysis of entire genomes of a pathogenic RNA virus enables comprehensive and concise identification of regulatory elements and emphasizes the extensive interrelationships among RNA genome structure, viral biology, and innate immune responses.

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.

cis-Acting RNA elements in the hepatitis C virus RNA genome

Hepatitis C virus (HCV) infection is a rapidly increasing global health problem with an estimated 170 million people infected worldwide. HCV is a hepatotropic, positive-sense RNA virus of the family Flaviviridae. As a positive-sense RNA virus, the HCV genome itself must serve as a template for translation, replication and packaging. The viral RNA must therefore be a dynamic structure that is able to readily accommodate structural changes to expose different regions of the genome to viral and cellular proteins to carry out the HCV life cycle. The ∼ 9600 nucleotide viral genome contains a single long open reading frame flanked by 5' and 3' non-coding regions that contain cis-acting RNA elements important for viral translation, replication and stability. Additional cis-acting RNA elements have also been identified in the coding sequences as well as in the 3' end of the negative-strand replicative intermediate. Herein, we provide an overview of the importance of these cis-acting RNA elements in the HCV life cycle.

Computational tools and algorithms for designing customized synthetic genes

Advances in DNA synthesis have enabled the construction of artificial genes, gene circuits, and genomes of bacterial scale. Freedom in de novo design of synthetic constructs provides significant power in studying the impact of mutations in sequence features, and verifying hypotheses on the functional information that is encoded in nucleic and amino acids. To aid this goal, a large number of software tools of variable sophistication have been implemented, enabling the design of synthetic genes for sequence optimization based on rationally defined properties. The first generation of tools dealt predominantly with singular objectives such as codon usage optimization and unique restriction site incorporation. Recent years have seen the emergence of sequence design tools that aim to evolve sequences toward combinations of objectives. The design of optimal protein-coding sequences adhering to multiple objectives is computationally hard, and most tools rely on heuristics to sample the vast sequence design space. In this review, we study some of the algorithmic issues behind gene optimization and the approaches that different tools have adopted to redesign genes and optimize desired coding features. We utilize test cases to demonstrate the efficiency of each approach, as well as identify their strengths and limitations.

Appearances can be deceptive: revealing a hidden viral infection with deep sequencing in a plant quarantine context

Comprehensive inventories of plant viral diversity are essential for effective quarantine and sanitation efforts. The safety of regulated plant material exchanges presently relies heavily on techniques such as PCR or nucleic acid hybridisation, which are only suited to the detection and characterisation of specific, well characterised pathogens. Here, we demonstrate the utility of sequence-independent next generation sequencing (NGS) of both virus-derived small interfering RNAs (siRNAs) and virion-associated nucleic acids (VANA) for the detailed identification and characterisation of viruses infecting two quarantined sugarcane plants. Both plants originated from Egypt and were known to be infected with Sugarcane streak Egypt Virus (SSEV; Genus Mastrevirus, Family Geminiviridae), but were revealed by the NGS approaches to also be infected by a second highly divergent mastrevirus, here named Sugarcane white streak Virus (SWSV). This novel virus had escaped detection by all routine quarantine detection assays and was found to also be present in sugarcane plants originating from Sudan. Complete SWSV genomes were cloned and sequenced from six plants and all were found to share >91% genome-wide identity. With the exception of two SWSV variants, which potentially express unusually large RepA proteins, the SWSV isolates display genome characteristics very typical to those of all other previously described mastreviruses. An analysis of virus-derived siRNAs for SWSV and SSEV showed them to be strongly influenced by secondary structures within both genomic single stranded DNA and mRNA transcripts. In addition, the distribution of siRNA size frequencies indicates that these mastreviruses are likely subject to both transcriptional and post-transcriptional gene silencing. Our study stresses the potential advantages of NGS-based virus metagenomic screening in a plant quarantine setting and indicates that such techniques could dramatically reduce the numbers of non-intercepted virus pathogens passing through plant quarantine stations.

Recombination in enteroviruses is a biphasic replicative process involving the generation of greater-than genome length 'imprecise' intermediates

Recombination in enteroviruses provides an evolutionary mechanism for acquiring extensive regions of novel sequence, is suggested to have a role in genotype diversity and is known to have been key to the emergence of novel neuropathogenic variants of poliovirus. Despite the importance of this evolutionary mechanism, the recombination process remains relatively poorly understood. We investigated heterologous recombination using a novel reverse genetic approach that resulted in the isolation of intermediate chimeric intertypic polioviruses bearing genomes with extensive duplicated sequences at the recombination junction. Serial passage of viruses exhibiting such imprecise junctions yielded progeny with increased fitness which had lost the duplicated sequences. Mutations or inhibitors that changed polymerase fidelity or the coalescence of replication complexes markedly altered the yield of recombinants (but did not influence non-replicative recombination) indicating both that the process is replicative and that it may be possible to enhance or reduce recombination-mediated viral evolution if required. We propose that extant recombinants result from a biphasic process in which an initial recombination event is followed by a process of resolution, deleting extraneous sequences and optimizing viral fitness. This process has implications for our wider understanding of ‘evolution by duplication’ in the positive-strand RNA viruses.

The rapid evolution of most positive-sense RNA viruses enables them to escape immune surveillance and adapt to new hosts. Genetic variation arises due to their error-prone RNA polymerases and by recombination of viral genomes in co-infected cells. We have developed a novel approach to analyse the poorly understood mechanism of recombination using a poliovirus model system. We characterised the initial viable recombinants and demonstrate the majority are longer than genome length due to an imprecise crossover event that duplicates part of the genome. These viruses are unfit, but rapidly lose the duplicated material and regain full fitness upon serial passage, a process we term resolution. We show this is a replicative recombination process by modifying the fidelity of the viral polymerase, or replication complex coalescence, using methods that have no influence on a previously reported, less efficient, non-replicative recombination mechanism. We conclude that recombination is a biphasic process involving separate generation and resolution events. These new insights into an important evolutionary mechanism have implications for our understanding of virus evolution through partial genome duplication, they suggest ways in which recombination might be modified and provides an approach that may be exploited to analyse recombination in other RNA viruses.

Structure-function relationship in viral RNA genomes: The case of hepatitis C virus

The acquisition of a storage information system beyond the nucleotide sequence has been a crucial issue for the propagation and dispersion of RNA viruses. This system is composed by highly conserved, complex structural units in the genomic RNA, termed functional RNA domains. These elements interact with other regions of the viral genome and/or proteins to direct viral translation, replication and encapsidation. The genomic RNA of the hepatitis C virus (HCV) is a good model for investigating about conserved structural units. It contains functional domains, defined by highly conserved structural RNA motifs, mostly located in the 5’-untranslatable regions (5’UTRs) and 3’UTR, but also occupying long stretches of the coding sequence. Viral translation initiation is mediated by an internal ribosome entry site located at the 5’ terminus of the viral genome and regulated by distal functional RNA domains placed at the 3’ end. Subsequent RNA replication strongly depends on the 3’UTR folding and is also influenced by the 5’ end of the HCV RNA. Further increase in the genome copy number unleashes the formation of homodimers by direct interaction of two genomic RNA molecules, which are finally packed and released to the extracellular medium. All these processes, as well as transitions between them, are controlled by structural RNA elements that establish a complex, direct and long-distance RNA-RNA interaction network. This review summarizes current knowledge about functional RNA domains within the HCV RNA genome and provides an overview of the control exerted by direct, long-range RNA-RNA contacts for the execution of the viral cycle.

Replication of many human viruses is refractory to inhibition by endogenous cellular microRNAs

The issue of whether viruses are subject to restriction by endogenous microRNAs (miRNAs) and/or by virus-induced small interfering RNAs (siRNAs) in infected human somatic cells has been controversial. Here, we address this question in two ways. First, using deep sequencing, we demonstrate that infection of human cells by the RNA virus dengue virus (DENV) or West Nile virus (WNV) does not result in the production of any virus-derived siRNAs or viral miRNAs. Second, to more globally assess the potential of small regulatory RNAs to inhibit virus replication, we used gene editing to derive human cell lines that lack a functional Dicer enzyme and that therefore are unable to produce miRNAs or siRNAs. Infection of these cells with a wide range of viruses, including DENV, WNV, yellow fever virus, Sindbis virus, Venezuelan equine encephalitis virus, measles virus, influenza A virus, reovirus, vesicular stomatitis virus, human immunodeficiency virus type 1, or herpes simplex virus 1 (HSV-1), failed to reveal any enhancement in the replication of any of these viruses, although HSV-1, which encodes at least eight Dicer-dependent viral miRNAs, did replicate somewhat more slowly in the absence of Dicer. We conclude that most, and perhaps all, human viruses have evolved to be resistant to inhibition by endogenous human miRNAs during productive replication and that dependence on a cellular miRNA, as seen with hepatitis C virus, is rare. How viruses have evolved to avoid inhibition by endogenous cellular miRNAs, which are generally highly conserved during metazoan evolution, remains to be determined. Importance: Eukaryotic cells express a wide range of small regulatory RNAs, including miRNAs, that have the potential to inhibit the expression of mRNAs that show sequence complementarity. Indeed, previous work has suggested that endogenous miRNAs have the potential to inhibit viral gene expression and replication. Here, we demonstrate that the replication of a wide range of pathogenic viruses is not enhanced in human cells engineered to be unable to produce miRNAs, indicating that viruses have evolved to be resistant to inhibition by miRNAs. This result is important, as it implies that manipulation of miRNA levels is not likely to prove useful in inhibiting virus replication. It also focuses attention on the question of how viruses have evolved to resist inhibition by miRNAs and whether virus mutants that have lost this resistance might prove useful, for example, in the development of attenuated virus vaccines.

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

RNA viruses infecting vertebrates differ fundamentally in their ability to establish persistent infections with markedly different patterns of transmission, disease mechanisms and evolutionary relationships with their hosts. Although interactions with host innate and adaptive responses are complex and persistence mechanisms likely multi-factorial, we previously observed associations between bioinformatically predicted RNA secondary formation in genomes of positive-stranded RNA viruses with their in vivo fitness and persistence. To analyse this interactions functionally, we transfected fibroblasts with non-replicating, non-translated RNA transcripts from RNA viral genomes with differing degrees of genome-scale ordered RNA structure (GORS). Single-stranded RNA transcripts induced interferon-β mediated though RIG-I and PKR activation, the latter associated with rapid induction of antiviral stress granules. A striking inverse correlation was observed between induction of both cellular responses with transcript RNA structure formation that was independent of both nucleotide composition and sequence length. The consistent inability of cells to recognize RNA transcripts possessing GORS extended to downstream differences from unstructured transcripts in expression of TNF-α, other interferon-stimulated genes and induction of apoptosis. This functional association provides novel insights into interactions between virus and host early after infection and provides evidence for a novel mechanism for evading intrinsic and innate immune responses.

Ribonuclease L and metal-ion-independent endoribonuclease cleavage sites in host and viral RNAs

Ribonuclease L (RNase L) is a metal-ion–independent endoribonuclease associated with antiviral and antibacterial defense, cancer and lifespan. Despite the biological significance of RNase L, the RNAs cleaved by this enzyme are poorly defined. In this study, we used deep sequencing methods to reveal the frequency and location of RNase L cleavage sites within host and viral RNAs. To make cDNA libraries, we exploited the 2′, 3′-cyclic phosphate at the end of RNA fragments produced by RNase L and other metal-ion–independent endoribonucleases. We optimized and validated 2′, 3′-cyclic phosphate cDNA synthesis and Illumina sequencing methods using viral RNAs cleaved with purified RNase L, viral RNAs cleaved with purified RNase A and RNA from uninfected and poliovirus-infected HeLa cells. Using these methods, we identified (i) discrete regions of hepatitis C virus and poliovirus RNA genomes that were profoundly susceptible to RNase L and other single-strand specific endoribonucleases, (ii) RNase L-dependent and RNase L-independent cleavage sites within ribosomal RNAs (rRNAs) and (iii) 2′, 3′-cyclic phosphates at the ends of 5S rRNA and U6 snRNA. Monitoring the frequency and location of metal-ion–independent endoribonuclease cleavage sites within host and viral RNAs reveals, in part, how these enzymes contribute to health and disease.