Detection of genome-scale ordered RNA structure (GORS) in genomes of positive-stranded RNA viruses: Implications for virus evolution and host persistence

Context-dependent structure formation of hairpin motifs in bacteriophage MS2 genomic RNA

Many functions of ribonucleic acid (RNA) rely on its ability to assume specific sequence-structure motifs. Packaging signals found in certain RNA viruses are one such prominent example of functional RNA motifs. These signals are short hairpin loops that interact with coat proteins and drive viral self-assembly. As they are found in different positions along the much longer genomic RNA, the formation of their correct structure occurs as a part of a larger context. Any changes to this context can consequently lead to changes in the structure of the motifs themselves. In fact, previous studies have shown that structure and function of RNA motifs can be highly context sensitive to the flanking sequence surrounding them. However, in what ways different flanking sequences influence the structure of an RNA motif they surround has yet to be studied in detail. We focus on a hairpin-rich region of the RNA genome of bacteriophage MS2-a well-studied RNA virus with a wide potential for use in biotechnology-and systematically examine context-dependent structural stability of 14 previously identified hairpin motifs, which include putative and confirmed packaging signals. Combining secondary and tertiary RNA structure prediction of the hairpin motifs placed in different contexts, ranging from the native genomic sequence to random RNA sequences and unstructured poly-U sequences, we determine different measures of motif structural stability. In this way, we show that while some motif structures can be stable in any context, others require specific context provided by the genome. Our results demonstrate the importance of context in RNA structure formation and how changes in the flanking sequence of an RNA motif sometimes lead to drastic changes in its structure. Structural stability of a motif in different contexts could provide additional insights into its functionality as well as assist in determining whether it remains functional when intentionally placed in other contexts.

Transformer model generated bacteriophage genomes are compositionally distinct from natural sequences

Novel applications of language models in genomics promise to have a large impact on the field. The megaDNA model is the first publicly available generative model for creating synthetic viral genomes. To evaluate megaDNA's ability to recapitulate the nonrandom genome composition of viruses and assess whether synthetic genomes can be algorithmically detected, compositional metrics for 4969 natural bacteriophage genomes and 1002 de novo synthetic bacteriophage genomes were compared. Transformer-generated sequences had varied but realistic genome lengths, and 58% were classified as viral by geNomad. However, the sequences demonstrated consistent differences in various compositional metrics when compared to natural bacteriophage genomes by rank-sum tests and principal component analyses. A simple neural network trained to detect transformer-generated sequences on global compositional metrics alone displayed a median sensitivity of 93.0% and specificity of 97.9% (n = 12 independent models). Overall, these results demonstrate that megaDNA does not yet generate bacteriophage genomes with realistic compositional biases and that genome composition is a reliable method for detecting sequences generated by this model. While the results are specific to the megaDNA model, the evaluated framework described here could be applied to any generative model for genomic sequences.

Molecular dynamic simulation reveals spider antimicrobial peptide Latarcin-1 and human eosinophil cationic protein as peptide inhibitors of SARS-CoV-2 variants

COVID-19 has rapidly proliferated around 180 countries, and new cases are reported frequently. No peptide medication has been developed that can reliably block SARS-CoV-2 infection. The investigation focuses on the crucial host receptors angiotensin-converting enzyme 2 (ACE2) , which can bind receptor-binding domain (RBD) on the SARS-CoV-2 spike protein (S). To investigate the inhibitory effects of human Eosinophil Cationic Protein (hECP) and Latarcin-1 (L1)on SARS-CoV-2 infection, we have selected them as research subjects. Further, we ran extensive molecular dynamics simulations to bring the docked peptide-ACE2 complex into its equilibrium state. The outcomes were then evaluated with g_MMPBSA and interaction analysis. We have also considered the Delta and Omicron variants to examine these peptides' inhibitory effects. The experimental findings revealed an enhanced capability of L1 and hECP as SARS-CoV-2 inhibitors, occupying hot spots and numerous key residues in ACE2. These include ASP30, ASP38, GLU35 and GLU75, which significantly inhibit the binding of RBD and ACE2 and are effective against two common variants in a similar manner. In addition, this study can serve as a springboard for future research on SARS-CoV-2 inhibitors.Communicated by Ramaswamy H. Sarma.

Base composition, adaptation, and evolution of goose astroviruses: codon-based investigation

Goose astroviruses (GoAstVs) are causative agents that account for fatal infection of goslings characterized by visceral urate deposition, resulting in severe economic losses in major goose-producing regions in China since 2017. In this study, we sought to unravel the intrinsic properties associated with adaptation and evolution in the host environment of GoAstVs. Consistent results from phylogenetic analysis and correspondence analysis performed on the codon usage patterns (CUPs) reveal 2 clusters of GoAstVs, namely, GoAstV-1 and GoAstV-2. However, multiple similar compositional characteristics were found, despite the high divergence between GoAstV-1 and GoAstV-2. Studies on the base composition of GoAstVs reveal an A/U bias, indicating a compositional constraint, while natural selection prevailed in determining the CUPs in the virus genome based on our neutrality plot analysis, reflecting high adaptive pressure to fit the host environment. Codon adaptation index (CAI) analysis revealed a higher degree of fitness to the CUPs of the corresponding host for GoAstVs than avian influenza virus and betacoronaviruses, which may be a favorable factor contributing to the high pathogenicity and wide distribution of GoAstVs in goslings. In addition, GoAstVs were less adapted to ducks and chickens, with significantly lower CAI values than to geese, which may be a reason for the different prevalence of GoAstVs among these species. Extensive investigations on dinucleotide distribution revealed a significant suppression of the CpG and UpA motifs in the virus genome, which may facilitate adaptation to the host's innate immune system by evading surveillance. In addition, our study reported the trends of increasing fitness to the host's microenvironment for GoAstVs through increasing adaptation to host CUPs and ongoing reduction of CpG motifs in the virus genome. The present analysis deepens our understanding of the basic biology, pathogenesis, adaptation and evolutionary pattern of GoAstVs, and contributes to the development of novel antiviral strategies.

Enterovirus A Shows Unique Patterns of Codon Usage Bias in Conventional Versus Unconventional Clade

Enterovirus A (EV-A) species cause hand, foot and mouth disease (HFMD), threatening the health of young children. Understanding the mutual codon usage pattern of the virus and its host(s) has fundamental and applied values. Here, through examining multiple codon usage parameters, we found that the codon usage bias among EV-A strains varies and is clade-specific. EVA76, EVA89, EVA90, EVA91 and EVA92, the unconventional clade of EV-A strains, show unique codon usage pattern relative to the two conventional clades, including EVA71, CVA16, CVA6 and CVA10, etc. Analyses of Effective Number of Codon (ENC), Correspondence Analysis (COA) and Parity Rule 2 (PR2), etc., revealed that the codon usage patterns of EV-A strains are shaped by mutation pressure and natural selection. Based on the neutrality analysis, we determined the dominant role of natural selection in the formation of the codon usage bias of EV-A. In addition, we have determined the codon usage compatibility of potential hosts for EV-A strains using codon adaptation index (CAI), relative codon deoptimization index (RCDI) and similarity index (SiD) analyses, and found that EV-A showed host-specific codon adaptation patterns in different clades. Finally, we confirmed that the unique codon usage pattern of the unconventional clade affected protein expression level in human cell lines. In conclusion, we identified novel characteristics of codon usage bias in distinct EV-A clades associated with their host range, transmission and pathogenicity.

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.

Transcription, mRNA export and immune evasion shape the codon usage of viruses

The nucleotide composition, dinucleotide composition, and codon usage of many viruses differs from their hosts. These differences arise because viruses are subject to unique mutation and selection pressures that do not apply to host genomes; however, the molecular mechanisms that underlie these evolutionary forces are unclear. Here, we analysed the patterns of codon usage in 1,520 vertebrate-infecting viruses, focusing on parameters known to be under selection and associated with gene regulation. We find that GC content, dinucleotide content, and splicing and m6A modification-related sequence motifs are associated with the type of genetic material (DNA or RNA), strandedness, and replication compartment of viruses. In an experimental follow-up, we find that the effects of GC content on gene expression depend on whether the genetic material is delivered to the cell as DNA or mRNA, whether it is transcribed by endogenous or exogenous RNA polymerase, and whether transcription takes place in the nucleus or cytoplasm. Our results suggest that viral codon usage cannot be explained by a simple adaptation to the codon usage of the host - instead, it reflects the combination of multiple selective and mutational pressures, including the need for efficient transcription, export, and immune evasion.

Insights into the evolutionary forces that shape the codon usage in the viral genome segments encoding intrinsically disordered protein regions

Intrinsically disordered regions/proteins (IDRs) are abundant across all the domains of life, where they perform important regulatory roles and supplement the biological functions of structured proteins/regions (SRs). Despite the multifunctionality features of IDRs, several interrogations on the evolution of viral genomic regions encoding IDRs in diverse viral proteins remain unreciprocated. To fill this gap, we benchmarked the findings of two most widely used and reliable intrinsic disorder prediction algorithms (IUPred2A and ESpritz) to a dataset of 6108 reference viral proteomes to unravel the multifaceted evolutionary forces that shape the codon usage in the viral genomic regions encoding for IDRs and SRs. We found persuasive evidence that the natural selection predominantly governs the evolution of codon usage in regions encoding IDRs by most of the viruses. In addition, we confirm not only that codon usage in regions encoding IDRs is less optimized for the protein synthesis machinery (transfer RNAs pool) of their host than for those encoding SRs, but also that the selective constraints imposed by codon bias sustain this reduced optimization in IDRs. Our analysis also establishes that IDRs in viruses are likely to tolerate more translational errors than SRs. All these findings hold true, irrespective of the disorder prediction algorithms used to classify IDRs. In conclusion, our study offers a novel perspective on the evolution of viral IDRs and the evolutionary adaptability to multiple taxonomically divergent hosts.

The dinucleotide composition of the Zika virus genome is shaped by conflicting evolutionary pressures in mammalian hosts and mosquito vectors

Most vertebrate RNA viruses show pervasive suppression of CpG and UpA dinucleotides, closely resembling the dinucleotide composition of host cell transcriptomes. In contrast, CpG suppression is absent in both invertebrate mRNA and RNA viruses that exclusively infect arthropods. Arthropod-borne (arbo) viruses are transmitted between vertebrate hosts by invertebrate vectors and thus encounter potentially conflicting evolutionary pressures in the different cytoplasmic environments. Using a newly developed Zika virus (ZIKV) model, we have investigated how demands for CpG suppression in vertebrate cells can be reconciled with potentially quite different compositional requirements in invertebrates and how this affects ZIKV replication and transmission. Mutant viruses with synonymously elevated CpG or UpA dinucleotide frequencies showed attenuated replication in vertebrate cell lines, which was rescued by knockout of the zinc-finger antiviral protein (ZAP). Conversely, in mosquito cells, ZIKV mutants with elevated CpG dinucleotide frequencies showed substantially enhanced replication compared to wild type. Host-driven effects on virus replication attenuation and enhancement were even more apparent in mouse and mosquito models. Infections with CpG- or UpA-high ZIKV mutants in mice did not cause typical ZIKV-induced tissue damage and completely protected mice during subsequent challenge with wild-type virus, which demonstrates their potential as live-attenuated vaccines. In contrast, the CpG-high mutants displayed enhanced replication in Aedes aegypti mosquitoes and a larger proportion of mosquitoes carried infectious virus in their saliva. These findings show that mosquito cells are also capable of discriminating RNA based on dinucleotide composition. However, the evolutionary pressure on the CpG dinucleotides of viral genomes in arthropod vectors directly opposes the pressure present in vertebrate host cells, which provides evidence that an adaptive compromise is required for arbovirus transmission. This suggests that the genome composition of arbo flaviviruses is crucial to maintain the balance between high-level replication in the vertebrate host and persistent replication in the mosquito vector.

The genome of the flavivirus Zika virus is stuck in a tug-of-war between two directly opposing evolutionary pressures that are present in the cells of mammalian host organisms and mosquito vectors; this results in an adaptive compromise, as manifested in the virus’s genome dinucleotide composition.

Covid-19 and the need for more history and philosophy of RNA

RNA is central to the COVID-19 pandemic-it shapes how the SARS Coronavirus 2 (SARS-CoV-2) behaves, and how researchers investigate and fight it. However, RNA has received relatively little attention in the history and philosophy of the life sciences. By analysing RNA biology in more detail, philosophers and historians of science could gain new and powerful tools to assess the current pandemic, and the biological sciences more generally.

A comprehensive profile of genomic variations in the SARS-CoV-2 isolates from the state of Telangana, India

The novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causing COVID-19 has rapidly turned into a pandemic, infecting millions and causing 1 157 509 (as of 27 October 2020) deaths across the globe. In addition to studying the mode of transmission and evasion of host immune system, analysing the viral mutational landscape constitutes an area under active research. The latter is expected to impart knowledge on the emergence of different clades, subclades, viral protein functions and protein-protein and protein-RNA interactions during replication/transcription cycle of virus and response to host immune checkpoints. In this study, we have attempted to bring forth the viral genomic variants defining the major clade(s) as identified from samples collected from the state of Telangana, India. We further report a comprehensive draft of all genomic variations (including unique mutations) present in SARS-CoV-2 strain in the state of Telangana. Our results reveal the presence of two mutually exclusive subgroups defined by specific variants within the dominant clade present in the population. This work attempts to bridge the critical gap regarding the genomic landscape and associate mutations in SARS-CoV-2 from a highly infected southern region of India, which was lacking to date.

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.

Cross-order host switches of hepatitis C-related viruses illustrated by a novel hepacivirus from sloths

The genealogy of the hepatitis C virus (HCV) and the genus Hepacivirus remains elusive despite numerous recently discovered animal hepaciviruses (HVs). Viruses from evolutionarily ancient mammals might elucidate the HV macro-evolutionary patterns. Here, we investigated sixty-seven two-toed and nine three-toed sloths from Costa Rica for HVs using molecular and serological tools. A novel sloth HV was detected by reverse transcription polymerase chain reaction (RT-PCR) in three-toed sloths (2/9, 22.2%; 95% confidence interval (CI), 5.3-55.7). Genomic characterization revealed typical HV features including overall polyprotein gene structure, a type 4 internal ribosomal entry site in the viral 5'-genome terminus, an A-U-rich region and X-tail structure in the viral 3'-genome terminus. Different from other animal HVs, HV seropositivity in two-toed sloths was low at 4.5 per cent (3/67; CI, 1.0-12.9), whereas the RT-PCR-positive three-toed sloths were seronegative. Limited cross-reactivity of the serological assay implied exposure of seropositive two-toed sloths to HVs of unknown origin and recent infections in RT-PCR-positive animals preceding seroconversion. Recent infections were consistent with only 9 nucleotide exchanges between the two sloth HVs, located predominantly within the E1/E2 encoding regions. Translated sequence distances of NS3 and NS5 proteins and host comparisons suggested that the sloth HV represents a novel HV species. Event- and sequence distance-based reconciliations of phylogenies of HVs and of their hosts revealed complex macro-evolutionary patterns, including both long-term evolutionary associations and host switches, most strikingly from rodents into sloths. Ancestral state reconstructions corroborated rodents as predominant sources of HV host switches during the genealogy of extant HVs. Sequence distance comparisons, partial conservation of critical amino acid residues associated with HV entry and selection pressure signatures of host genes encoding entry and antiviral protein orthologs were consistent with HV host switches between genetically divergent mammals, including the projected host switch from rodents into sloths. Structural comparison of HCV and sloth HV E2 proteins suggested conserved modes of hepaciviral entry. Our data corroborate complex macro-evolutionary patterns shaping the genus Hepacivirus, highlight that host switches are possible across highly diverse host taxa, and elucidate a prominent role of rodent hosts during the Hepacivirus genealogy.

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.

Evolutionary Analyses of Base-Pairing Interactions in DNA and RNA Secondary Structures

Pairs of nucleotides within functional nucleic acid secondary structures often display evidence of coevolution that is consistent with the maintenance of base-pairing. Here, we introduce a sequence evolution model, MESSI (Modeling the Evolution of Secondary Structure Interactions), that infers coevolution associated with base-paired sites in DNA or RNA sequence alignments. MESSI can estimate coevolution while accounting for an unknown secondary structure. MESSI can also use graphics processing unit parallelism to increase computational speed. We used MESSI to infer coevolution associated with GC, AU (AT in DNA), GU (GT in DNA) pairs in noncoding RNA alignments, and in single-stranded RNA and DNA virus alignments. Estimates of GU pair coevolution were found to be higher at base-paired sites in single-stranded RNA viruses and noncoding RNAs than estimates of GT pair coevolution in single-stranded DNA viruses. A potential biophysical explanation is that GT pairs do not stabilize DNA secondary structures to the same extent that GU pairs do in RNA. Additionally, MESSI estimates the degrees of coevolution at individual base-paired sites in an alignment. These estimates were computed for a SHAPE-MaP-determined HIV-1 NL4-3 RNA secondary structure. We found that estimates of coevolution were more strongly correlated with experimentally determined SHAPE-MaP pairing scores than three nonevolutionary measures of base-pairing covariation. To assist researchers in prioritizing substructures with potential functionality, MESSI automatically ranks substructures by degrees of coevolution at base-paired sites within them. Such a ranking was created for an HIV-1 subtype B alignment, revealing an excess of top-ranking substructures that have been previously identified as having structure-related functional importance, among several uncharacterized top-ranking substructures.

Interaction of virus populations with their hosts

Viral population numbers are extremely large compared with those of their host species. Population bottlenecks are frequent during the life cycle of viruses and can reduce viral populations transiently to very few individuals. Viruses have to confront several types of constraints that can be divided into basal, cell-dependent, and organism-dependent constraints. Viruses overcome them exploiting a number of molecular mechanisms, with an important contribution of population numbers and genome variation. The adaptive potential of viruses is reflected in modifications of cell tropism and host range, escape to components of the host immune response, and capacity to alternate among different host species, among other phenotypic changes. Despite a fitness cost of most mutations required to overcome a selective constraint, viruses can find evolutionary pathways that ensure their survival in equilibrium with their hosts.

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.

Human pegivirus 2 exhibits minimal geographic and temporal genetic diversity

We applied an NGS based target capture approach to amplify HPgV-2 sequences from metagenomic libraries and enable full genome characterization. Despite expanded geographical sampling, sequence variability remains low, with diversity concentrated in approximately 3.3% of all amino acids. Serial samples from one HPgV-2 positive individual co-infected with comparable titers of HIV, HCV, and GBV-C showed that HPgV-2 remains highly stable over several weeks compared to other RNA viruses, despite a similarly error-prone polymerase. The consistent epidemiological association with and structural similarities to HCV, and the weak positive correlation of HCV and HPgV-2 titers shown here, suggests it may benefit from co-infection. While minimal selective pressure on HPgV-2 to evolve could suggest fitness, the rarity of HPgV-2 and the tight phylogenetic clustering of global strains likely indicates origination from a common source and a virus that is ill-suited to its host. Sporadic infections may explain the limited genetic diversity observed worldwide.

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.

Conserved Secondary Structures in Viral mRNAs

RNA secondary structure in untranslated and protein coding regions has been shown to play an important role in regulatory processes and the viral replication cycle. While structures in non-coding regions have been investigated extensively, a thorough overview of the structural repertoire of protein coding mRNAs, especially for viruses, is lacking. Secondary structure prediction of large molecules, such as long mRNAs remains a challenging task, as the contingent of structures a sequence can theoretically fold into grows exponentially with sequence length. We applied a structure prediction pipeline to Viral Orthologous Groups that first identifies the local boundaries of potentially structured regions and subsequently predicts their functional importance. Using this procedure, the orthologous groups were split into structurally homogenous subgroups, which we call subVOGs. This is the first compilation of potentially functional conserved RNA structures in viral coding regions, covering the complete RefSeq viral database. We were able to recover structural elements from previous studies and discovered a variety of novel structured regions. The subVOGs are available through our web resource RNASIV (RNA structure in viruses; http://rnasiv.bio.wzw.tum.de).

Differential interaction strategies of hepatitis c virus genotypes during entry - An in silico investigation of envelope glycoprotein E2 - CD81 interaction

Hepatitis C Virus is a blood borne pathogen responsible for chronic hepatitis in more than 71 million people. Wide variations across strains and genotypes are one of the major hurdles in therapeutic development. While genotype 1 remains the most extensively studied and abundant strain, genotype 3 is more virulent and second most prevalent. This study aimed to compare differences in the glycoprotein E2 across HCV genotypes at nucleotide, protein and structural levels. Nucleotide sequences of E2 from 29 strains across genotypes 1a, 1b, 3a and 3b revealed a stark preference for C-richness which was attributed to a distinct bias for C-rich codons in genotype 1. Genotype 3 exhibited a similar preference to a lesser extent. Amino acid level comparison revealed majority of the changes at the C-terminal half of the proteins leaving the N-terminal region conspicuously conserved apart from the two hyper variable regions. Amino acid changes across genotypes were mostly polar-nonpolar alterations. In silico models of E2 glycoproteins and docking analysis with the energy minimized PDB-CD81 model revealed unique interacting residues in both E2 and CD81. While several CD81 binding residues were common for all four genotypes, number and composition of interacting residues varied. The interacting residues of E2 were however unique for each genotype. E2 of genotype 3a and CD81 had the strongest interaction. In conclusion this is the first comprehensive study comparing E2 sequences across genotypes 1a, 1b, 3a and 3b revealing stark genotype-specific differences which requires more extensive investigation.

Evolutionary effects of the AID/APOBEC family of mutagenic enzymes on human gamma-herpesviruses

The human gamma-herpesviruses, Epstein-Barr virus and Kaposi's sarcoma-associated herpesvirus, establish lifelong latency in B cells and are associated with multiple malignancies. Virus-host coevolution often drive changes in both host immunity and in the viral genome. We consider one host immune mechanism, the activation-induced deaminase (AID)/APOBEC family of cytidine deaminases, that induces mutations in viral DNA. AID, the ancestral gene in the family has a conserved role in somatic hypermutation, a key step in antibody affinity maturation. The APOBEC3 subfamily, of which there are seven genes in human, have evolved antiviral functions and have diversified in terms of their expression pattern, subcellular localization, and DNA mutation motifs (hotspots). In this study, we investigated how the human gamma-herpesviruses have evolved to avoid the action of the AID/APOBEC enzymes and determine if these enzymes are contributing to the ongoing evolution of the viruses. We used computational methods to evaluate observed versus expected frequency of AID/APOBEC hotspots in viral genomes and found that the viruses have evolved to limit the representation of AID and certain APOBEC3 motifs. At the same time, the remaining hotspots were highly likely to cause amino acid changes, suggesting prolonged evolutionary pressure of the enzymes on the viruses. To study current hypermutation, as opposed to historical mutation processes, we also analyzed putative mutations derived from alignments of published viral genomes and found again that AID and APOBEC3 appear to target the genome most frequently. New protein variants resulting from AID/APOBEC activity may have important consequences in health, including vaccine development (epitope evolution) and host immune evasion.

Structural and functional analysis of the roles of the HCV 5' NCR miR122-dependent long-range association and SLVI in genome translation and replication

The hepatitis C virus RNA genome possesses a variety of conserved structural elements, in both coding and non-coding regions, that are important for viral replication. These elements are known or predicted to modulate key life cycle events, such as translation and genome replication, some involving conformational changes induced by long-range RNA–RNA interactions. One such element is SLVI, a stem-loop (SL) structure located towards the 5′ end of the core protein-coding region. This element forms an alternative RNA–RNA interaction with complementary sequences in the 5′ untranslated regions that are independently involved in the binding of the cellular microRNA 122 (miR122). The switch between ‘open’ and ‘closed’ structures involving SLVI has previously been proposed to modulate translation, with lower translation efficiency associated with the ‘closed’ conformation. In the current study, we have used selective 2′-hydroxyl acylation analysed by primer extension to validate this RNA–RNA interaction in the absence and presence of miR122. We show that the long-range association (LRA) only forms in the absence of miR122, or otherwise requires the blocking of miR122 binding combined with substantial disruption of SLVI. Using site-directed mutations introduced to promote open or closed conformations of the LRA we demonstrate no correlation between the conformation and the translation phenotype. In addition, we observed no influence on virus replication compared to unmodified genomes. The presence of SLVI is well-documented to suppress translation, but these studies demonstrate that this is not due to its contribution to the LRA. We conclude that, although there are roles for SLVI in translation, the LRA is not a riboswitch regulating the translation and replication phenotypes of the virus.

Evolution of Codon Usage Bias in Henipaviruses Is Governed by Natural Selection and Is Host-Specific

Hendra virus (HeV) and Nipah virus (NiV) are among a group of emerging bat-borne paramyxoviruses that have crossed their species-barrier several times by infecting several hosts with a high fatality rate in human beings. Despite the fatal nature of their infection, a comprehensive study to explore their evolution and adaptation in different hosts is lacking. A study of codon usage patterns in henipaviruses may provide some fruitful insight into their evolutionary processes of synonymous codon usage and host-adapted evolution. Here, we performed a systematic evolutionary and codon usage bias analysis of henipaviruses. We found a low codon usage bias in the coding sequences of henipaviruses and that natural selection, mutation pressure, and nucleotide compositions shapes the codon usage patterns of henipaviruses, with natural selection being more important than the others. Also, henipaviruses showed the highest level of adaptation to bats of the genus Pteropus in the codon adaptation index (CAI), relative to the codon de-optimization index (RCDI), and similarity index (SiD) analyses. Furthermore, a comparison to recently identified henipa-like viruses indicated a high tRNA adaptation index of henipaviruses for human beings, mainly due to F, G and L proteins. Consequently, the study concedes the substantial emergence of henipaviruses in human beings, particularly when paired with frequent exposure to direct/indirect bat excretions.

Mechanisms and consequences of positive-strand RNA virus recombination

Genetic recombination in positive-strand RNA viruses is a significant evolutionary mechanism that drives the creation of viral diversity by the formation of novel chimaeric genomes. The process and its consequences, for example the generation of viruses with novel phenotypes, has historically been studied by analysis of the end products. More recently, with an appreciation that there are both replicative and non-replicative mechanisms at work, and with new approaches and techniques to analyse intermediate products, the viral and cellular factors that influence the process are becoming understood. The major influence on replicative recombination is the fidelity of viral polymerase, although RNA structures and sequences may also have an impact. In replicative recombination the viral polymerase is necessary and sufficient, although roles for other viral or cellular proteins may exist. In contrast, non-replicative recombination appears to be mediated solely by cellular components. Despite these insights, the relative importance of replicative and non-replicative mechanisms is not clear. Using single-stranded positive-sense RNA viruses as exemplars, we review the current state of understanding of the processes and consequences of recombination.

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.

Probing of RNA structures in a positive sense RNA virus reveals selection pressures for structural elements

In single stranded (+)-sense RNA viruses, RNA structural elements (SEs) play essential roles in the infection process from replication to encapsidation. Using selective 2'-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq) and covariation analysis, we explore the structural features of the third genome segment of cucumber mosaic virus (CMV), RNA3 (2216 nt), both in vitro and in plant cell lysates. Comparing SHAPE-Seq and covariation analysis results revealed multiple SEs in the coat protein open reading frame and 3' untranslated region. Four of these SEs were mutated and serially passaged in Nicotiana tabacum plants to identify biologically selected changes to the original mutated sequences. After passaging, loop mutants showed partial reversion to their wild-type sequence and SEs that were structurally disrupted by mutations were restored to wild-type-like structures via synonymous mutations in planta. These results support the existence and selection of virus open reading frame SEs in the host organism and provide a framework for further studies on the role of RNA structure in viral infection. Additionally, this work demonstrates the applicability of high-throughput chemical probing in plant cell lysates and presents a new method for calculating SHAPE reactivities from overlapping reverse transcriptase priming sites.

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.

Resistance of Echovirus 11 to ClO2 Is Associated with Enhanced Host Receptor Use, Altered Entry Routes, and High Fitness

Waterborne viruses can exhibit resistance to common water disinfectants, yet the mechanisms that allow them to tolerate disinfection are poorly understood. Here, we generated echovirus 11 (E11) with resistance to chlorine dioxide (ClO₂) by experimental evolution, and we assessed the associated genotypic and phenotypic traits. ClO₂ resistance emerged after E11 populations were repeatedly reduced (either by ClO₂-exposure or by dilution) and then regrown in cell culture. The resistance was linked to an improved capacity of E11 to bind to its host cells, which was further attributed to two potential causes: first, the resistant E11 populations possessed mutations that caused amino acid substitutions from ClO₂-labile to ClO₂-stable residues in the viral proteins, which likely increased the chemical stability of the capsid toward ClO₂. Second, resistant E11 mutants exhibited the capacity to utilize alternative cell receptors for host binding. Interestingly, the emergence of ClO₂ resistance resulted in an enhanced replicative fitness compared to the less resistant starting population. Overall this study contributes to a better understanding of the mechanism underlying disinfection resistance in waterborne viruses, and processes that drive resistance development.

Differential Infection Patterns and Recent Evolutionary Origins of Equine Hepaciviruses in Donkeys

The hepatitis C virus (HCV) is a major human pathogen. Genetically related viruses in animals suggest a zoonotic origin of HCV. The closest relative of HCV is found in horses (termed equine hepacivirus [EqHV]). However, low EqHV genetic diversity implies relatively recent acquisition of EqHV by horses, making a derivation of HCV from EqHV unlikely. To unravel the EqHV evolutionary history within equid sister species, we analyzed 829 donkeys and 53 mules sampled in nine European, Asian, African, and American countries by molecular and serologic tools for EqHV infection. Antibodies were found in 278 animals (31.5%), and viral RNA was found in 3 animals (0.3%), all of which were simultaneously seropositive. A low RNA prevalence in spite of high seroprevalence suggests a predominance of acute infection, a possible difference from the mostly chronic hepacivirus infection pattern seen in horses and humans. Limitation of transmission due to short courses of infection may explain the existence of entirely seronegative groups of animals. Donkey and horse EqHV strains were paraphyletic and 97.5 to 98.2% identical in their translated polyprotein sequences, making virus/host cospeciation unlikely. Evolutionary reconstructions supported host switches of EqHV between horses and donkeys without the involvement of adaptive evolution. Global admixture of donkey and horse hepaciviruses was compatible with anthropogenic alterations of EqHV ecology. In summary, our findings do not support EqHV as the origin of the significantly more diversified HCV. Identification of a host system with predominantly acute hepacivirus infection may enable new insights into the chronic infection pattern associated with HCV.

IMPORTANCE

The evolutionary origins of the human hepatitis C virus (HCV) are unclear. The closest animal-associated relative of HCV occurs in horses (equine hepacivirus [EqHV]). The low EqHV genetic diversity implies a relatively recent acquisition of EqHV by horses, limiting the time span for potential horse-to-human infections in the past. Horses are genetically related to donkeys, and EqHV may have cospeciated with these host species. Here, we investigated a large panel of donkeys from various countries using serologic and molecular tools. We found EqHV to be globally widespread in donkeys and identify potential differences in EqHV infection patterns, with donkeys potentially showing enhanced EqHV clearance compared to horses. We provide strong evidence against EqHV cospeciation and for its capability to switch hosts among equines. Differential hepacivirus infection patterns in horses and donkeys may enable new insights into the chronic infection pattern associated with HCV.

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.

The influenza A virus NS genome segment displays lineage-specific patterns in predicted RNA secondary structure

BACKGROUND

Influenza A virus (IAV) is a segmented negative-sense RNA virus that causes seasonal epidemics and periodic pandemics in humans. Two regions (nucleotide positions 82-148 and 497-564) in the positive-sense RNA of the NS segment fold into a multi-branch loop or hairpin structures.

RESULTS

We studied 25,384 NS segment positive-sense RNA unique sequences of human and non-human IAVs in order to predict secondary RNA structures of the 82-148 and 497-564 regions using RNAfold software, and determined their host- and lineage-specific distributions. Hairpins prevailed in avian and avian-origin human IAVs, including H1N1pdm1918 and H5N1. In human and swine IAV hairpins distribution varied between evolutionary lineages.

CONCLUSIONS

These results suggest a possible functional role for these RNA secondary structures and the need for experimental evaluation of these structures in the influenza life cycle.

Identification of novel RNA secondary structures within the hepatitis C virus genome reveals a cooperative involvement in genome packaging

The specific packaging of the hepatitis C virus (HCV) genome is hypothesised to be driven by Core-RNA interactions. To identify the regions of the viral genome involved in this process, we used SELEX (systematic evolution of ligands by exponential enrichment) to identify RNA aptamers which bind specifically to Core in vitro. Comparison of these aptamers to multiple HCV genomes revealed the presence of a conserved terminal loop motif within short RNA stem-loop structures. We postulated that interactions of these motifs, as well as sub-motifs which were present in HCV genomes at statistically significant levels, with the Core protein may drive virion assembly. We mutated 8 of these predicted motifs within the HCV infectious molecular clone JFH-1, thereby producing a range of mutant viruses predicted to possess altered RNA secondary structures. RNA replication and viral titre were unaltered in viruses possessing only one mutated structure. However, infectivity titres were decreased in viruses possessing a higher number of mutated regions. This work thus identified multiple novel RNA motifs which appear to contribute to genome packaging. We suggest that these structures act as cooperative packaging signals to drive specific RNA encapsidation during HCV assembly.

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.

Interaction of Virus Populations with Their Hosts

Viral population numbers are extremely large compared with those of their host species. Population bottlenecks are frequent during the life cycle of viruses and can reduce viral populations transiently to very few individuals. Viruses have to confront several types of constraints that can be divided in basal, cell-dependent, and organism-dependent constraints. Viruses overcome them exploiting a number of molecular mechanisms, with an important contribution of population numbers and genome variation. The adaptive potential of viruses is reflected in modifications of cell tropism and host range, escape to components of the host immune response, and capacity to alternate among different host species, among other phenotypic changes. Despite a fitness cost of most mutations required to overcome a selective constraint, viruses can find evolutionary pathways that ensure their survival in equilibrium with their hosts.

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.

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.

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.

Widespread signatures of local mRNA folding structure selection in four Dengue virus serotypes

BACKGROUND

It is known that mRNA folding can affect and regulate various gene expression steps both in living organisms and in viruses. Previous studies have recognized functional RNA structures in the genome of the Dengue virus. However, these studies usually focused either on the viral untranslated regions or on very specific and limited regions at the beginning of the coding sequences, in a limited number of strains, and without considering evolutionary selection.

RESULTS

Here we performed the first large scale comprehensive genomics analysis of selection for local mRNA folding strength in the Dengue virus coding sequences, based on a total of 1,670 genomes and 4 serotypes. Our analysis identified clusters of positions along the coding regions that may undergo a conserved evolutionary selection for strong or weak local folding maintained across different viral variants. Specifically, 53-66 clusters for strong folding and 49-73 clusters for weak folding (depending on serotype) aggregated of positions with a significant conservation of folding energy signals (related to partially overlapping local genomic regions) were recognized. In addition, up to 7% of these positions were found to be conserved in more than 90% of the viral genomes. Although some of the identified positions undergo frequent synonymous / non-synonymous substitutions, the selection for folding strength therein is preserved, and thus cannot be trivially explained based on sequence conservation alone.

CONCLUSIONS

The fact that many of the positions with significant folding related signals are conserved among different Dengue variants suggests that a better understanding of the mRNA structures in the corresponding regions may promote the development of prospective anti- Dengue vaccination strategies. The comparative genomics approach described here can be employed in the future for detecting functional regions in other pathogens with very high mutations rates.

CodonShuffle: a tool for generating and analyzing synonymously mutated sequences

Because synonymous mutations do not change the amino acid sequence of a protein, they are generally considered to be selectively neutral. Empiric data suggest, however, that a significant fraction of viral mutational fitness effects may be attributable to synonymous mutation. Bias in synonymous codon usage in viruses may result from selection for translational efficiency, mutational bias, base pairing requirements in RNA structures, or even selection against specific dinucleotides by innate immune effectors. Experimental analyses of codon usage and genome evolution have been facilitated by advances in synthetic biology, which now make it feasible to generate viral genomes that contain large numbers of synonymous mutations. The generally pleiotropic effects of synonymous mutation on viral fitness have, at times, made it difficult to define the mechanistic basis for the observed attenuation of these heavily mutated viruses. We have addressed this problem by developing a bioinformatic tool for the generation and analysis of viral sequences with large-scale synonymous mutation. A variety of permutation strategies are applied to shuffle codons within an open reading frame. After measuring the dinucleotide frequency, codon usage, codon pair bias, and free energy of RNA folding for each permuted genome, we used z-score normalization and a least squares regression model to quantify their overall distance from the starting sequence. Using this approach, the user can easily identify a large number of synonymously mutated sequences with varying similarity to a wild-type genome across a range of nucleic-acid-based determinants of viral fitness. We believe that this tool will be useful in designing genomes for subsequent experimental studies of the fitness impacts of synonymous mutation.