A cis-acting replication element in the sequence encoding the NS5B RNA-dependent RNA polymerase is required for hepatitis C virus RNA replication

Conformational dynamics of the hepatitis C virus 3'X RNA

The 3' end of the hepatitis C virus genome is terminated by a highly conserved, 98 nt sequence called 3'X. This untranslated structural element is thought to regulate several essential RNA-dependent processes associated with infection. 3'X has two proposed conformations comprised of either three or two stem-loop structures that result from the different base-pairing interactions within the first 55 nt. Here, we used single-molecule Förster resonance energy transfer spectroscopy to monitor the conformational status of fluorescently labeled constructs that isolate this region of the RNA (3'X55). We observed that 3'X55 can adopt both proposed conformations and the relative abundance of them can be modulated by either solution conditions or nucleotide deletions. Furthermore, interconversion between the two conformations takes place over the course of several hours. The simultaneous existence of two slowly interconverting conformations may help prime individual copies of the viral genome for either viral protein or RNA synthesis, thereby minimizing conflicts between these two competing processes.

Comprehensive survey of conserved RNA secondary structures in full-genome alignment of Hepatitis C virus

Hepatitis C virus (HCV) is a plus-stranded RNA virus that often chronically infects liver hepatocytes and causes liver cirrhosis and cancer. These viruses replicate their genomes employing error-prone replicases. Thereby, they routinely generate a large 'cloud' of RNA genomes (quasispecies) which-by trial and error-comprehensively explore the sequence space available for functional RNA genomes that maintain the ability for efficient replication and immune escape. In this context, it is important to identify which RNA secondary structures in the sequence space of the HCV genome are conserved, likely due to functional requirements. Here, we provide the first genome-wide multiple sequence alignment (MSA) with the prediction of RNA secondary structures throughout all representative full-length HCV genomes. We selected 57 representative genomes by clustering all complete HCV genomes from the BV-BRC database based on k-mer distributions and dimension reduction and adding RefSeq sequences. We include annotations of previously recognized features for easy comparison to other studies. Our results indicate that mainly the core coding region, the C-terminal NS5A region, and the NS5B region contain secondary structure elements that are conserved beyond coding sequence requirements, indicating functionality on the RNA level. In contrast, the genome regions in between contain less highly conserved structures. The results provide a complete description of all conserved RNA secondary structures and make clear that functionally important RNA secondary structures are present in certain HCV genome regions but are largely absent from other regions. Full-genome alignments of all branches of Hepacivirus C are provided in the supplement.

Hepatitis C virus NS5B triggers an MDA5-mediated innate immune response by producing dsRNA without the replication of viral genomes

During the replication of viral genomes, RNA viruses produce double-stranded RNA (dsRNA), through the activity of their RNA-dependent RNA polymerases (RdRps) as viral replication intermediates. Recognition of viral dsRNA by host pattern recognition receptors - such as retinoic acid-induced gene-I (RIG-I)-like receptors and Toll-like receptor 3 - triggers the production of interferon (IFN)-β via the activation of IFN regulatory factor (IRF)-3. It has been proposed that, during the replication of viral genomes, each of RIG-I and melanoma differentiation-associated gene 5 (MDA5) form homodimers for the efficient activation of a downstream signalling pathway in host cells. We previously reported that, in the non-neoplastic human hepatocyte line PH5CH8, the RdRp NS5B derived from hepatitis C virus (HCV) could induce IFN-β expression by its RdRp activity without the actual replication of viral genomes. However, the exact mechanism by which HCV NS5B produced IFN-β remained unknown. In the present study, we first showed that NS5B derived from another Flaviviridae family member, GB virus B (GBV-B), also possessed the ability to induce IFN-β in PH5CH8 cells. Similarly, HCV NS5B, but not its G317V mutant, which lacks RdRp activity, induced the dimerization of MDA5 and subsequently the activation of IRF-3. Interestingly, immunofluorescence analysis showed that HCV NS5B produced dsRNA. Like HCV NS5B, GBV-B NS5B also triggered the production of dsRNA and subsequently the dimerization of MDA5. Taken together, our results show that HCV NS5B triggers an MDA5-mediated innate immune response by producing dsRNA without the replication of viral genomes in human hepatocytes.

Dimerization of an umbravirus RNA genome activates subgenomic mRNA transcription

Many eukaryotic RNA viruses transcribe subgenomic (sg) mRNAs during infections to control expression of a subset of viral genes. Such transcriptional events are commonly regulated by local or long-range intragenomic interactions that form higher-order RNA structures within these viral genomes. In contrast, here we report that an umbravirus activates sg mRNA transcription via base pair-mediated dimerization of its plus-strand RNA genome. Compelling in vivo and in vitro evidence demonstrate that this viral genome dimerizes via a kissing-loop interaction involving an RNA stem-loop structure located just upstream from its transcriptional initiation site. Both specific and non-specific features of the palindromic kissing-loop complex were found to contribute to transcriptional activation. Structural and mechanistic aspects of the process in umbraviruses are discussed and compared with genome dimerization events in other RNA viruses. Notably, probable dimer-promoting RNA stem-loop structures were also identified in a diverse group of umbra-like viruses, suggesting broader utilization of this unconventional transcriptional strategy.

IndelsRNAmute: predicting deleterious multiple point substitutions and indels mutations

Background

RNA deleterious point mutation prediction was previously addressed with programs such as RNAmute and MultiRNAmute. The purpose of these programs is to predict a global conformational rearrangement of the secondary structure of a functional RNA molecule, thereby disrupting its function. RNAmute was designed to deal with only single point mutations in a brute force manner, while in MultiRNAmute an efficient approach to deal with multiple point mutations was developed. The approach used in MultiRNAmute is based on the stabilization of the suboptimal RNA folding prediction solutions and/or destabilization of the optimal folding prediction solution of the wild type RNA molecule. The MultiRNAmute algorithm is significantly more efficient than the brute force approach in RNAmute, but in the case of long sequences and large m-point mutation sets the MultiRNAmute becomes exponential in examining all possible stabilizing and destabilizing mutations.

Results

An inherent limitation in the RNAmute and MultiRNAmute programs is their ability to predict only substitution mutations, as these programs were not designed to work with deletion or insertion mutations. To address this limitation we herein develop a very fast algorithm, based on suboptimal folding solutions, to predict a predefined number of multiple point deleterious mutations as specified by the user. Depending on the user’s choice, each such set of mutations may contain combinations of deletions, insertions and substitution mutations. Additionally, we prove the hardness of predicting the most deleterious set of point mutations in structural RNAs.

Conclusions

We developed a method that extends our previous MultiRNAmute method to predict insertion and deletion mutations in addition to substitutions. The additional advantage of the new method is its efficiency to find a predefined number of deleterious mutations. Our new method may be exploited by biologists and virologists prior to site-directed mutagenesis experiments, which involve indel mutations along with substitutions. For example, our method may help to investigate the change of function in an RNA virus via mutations that disrupt important motifs in its secondary structure.

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.

Experimental demonstration and pan-structurome prediction of climate-associated riboSNitches in Arabidopsis

Background

Genome-wide association studies (GWAS) aim to correlate phenotypic changes with genotypic variation. Upon transcription, single nucleotide variants (SNVs) may alter mRNA structure, with potential impacts on transcript stability, macromolecular interactions, and translation. However, plant genomes have not been assessed for the presence of these structure-altering polymorphisms or “riboSNitches.”

Results

We experimentally demonstrate the presence of riboSNitches in transcripts of two Arabidopsis genes, ZINC RIBBON 3 (ZR3) and COTTON GOLGI-RELATED 3 (CGR3), which are associated with continentality and temperature variation in the natural environment. These riboSNitches are also associated with differences in the abundance of their respective transcripts, implying a role in regulating the gene's expression in adaptation to local climate conditions. We then computationally predict riboSNitches transcriptome-wide in mRNAs of 879 naturally inbred Arabidopsis accessions. We characterize correlations between SNPs/riboSNitches in these accessions and 434 climate descriptors of their local environments, suggesting a role of these variants in local adaptation. We integrate this information in CLIMtools V2.0 and provide a new web resource, T-CLIM, that reveals associations between transcript abundance variation and local environmental variation.

Conclusion

We functionally validate two plant riboSNitches and, for the first time, demonstrate riboSNitch conditionality dependent on temperature, coining the term “conditional riboSNitch.” We provide the first pan-genome-wide prediction of riboSNitches in plants. We expand our previous CLIMtools web resource with riboSNitch information and with 1868 additional Arabidopsis genomes and 269 additional climate conditions, which will greatly facilitate in silico studies of natural genetic variation, its phenotypic consequences, and its role in local adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02656-4.

OUP accepted manuscript

Subdomain 5BSL3.2 of hepatitis C virus RNA lies at the core of a network of distal RNA–RNA contacts that connect the 5′ and 3′ regions of the viral genome and regulate the translation and replication stages of the viral cycle. Using small-angle X-ray scattering and NMR spectroscopy experiments, we have determined at low resolution the structural models of this subdomain and its distal complex with domain 3′X, located at the 3′-terminus of the viral RNA chain. 5BSL3.2 adopts a characteristic ‘L’ shape in solution, whereas the 5BSL3.2–3′X distal complex forms a highly unusual ‘Y’-shaped kissing junction that blocks the dimer linkage sequence of domain 3′X and promotes translation. The structure of this complex may impede an effective association of the viral polymerase with 5BSL3.2 and 3′X to start negative-strand RNA synthesis, contributing to explain the likely mechanism used by these sequences to regulate viral replication and translation. In addition, sequence and shape features of 5BSL3.2 are present in functional RNA motifs of flaviviruses, suggesting conserved regulatory processes within the Flaviviridae family.

Development of full-length cell-culture infectious clone and subgenomic replicon for a genotype 3a isolate of hepatitis C virus

Hepatitis C virus (HCV) genotype 3 is widely distributed, and genotype 3-infected patients achieve a lower cure rate in direct-acting antiviral (DAA) therapy and are associated with a higher risk of hepatic steatosis than patients with other genotypes. Thus, the study of the virology and pathogenesis of genotype 3 HCV is increasingly relevant. Here, we developed a full-length infectious clone and a subgenomic replicon for the genotype 3a isolate, CH3a. From an infected serum, we constructed a full-length CH3a clone, however, it was nonviable in Huh7.5.1 cells. Next, we systematically adapted several intergenotypic recombinants containing Core-NS2 and 5′UTR-NS5A from CH3a, and other sequences from a replication-competent genotype 2 a clone JFH1. Adaptive mutations were identified, of which several combinations facilitated the replication of CH3a-JFH1 recombinants; however, they failed to adapt to the full-length CH3a and the recombinants containing CH3a NS5B. Thus, we attempted to separately adapt CH3a NS5B-3′UTR by constructing an intragenotypic recombinant using 5′UTR-NS5A from an infectious genotype 3a clone, DBN3acc, from which L3004P/M in NS5B and a deletion of 11 nucleotides (Δ11nt) downstream of the polyU/UC tract of the 3′UTR were identified and demonstrated to efficiently improve virus production. Finally, we combined functional 5′UTR-NS5A and NS5B-3′UTR sequences that carried the selected mutations to generate full-length CH3a with 26 or 27 substitutions (CH3acc), and both revealed efficient replication and virus spread in transfected and infected cells, releasing HCV of 104.2 f.f.u. ml−1. CH3acc was inhibited by DAAs targeting NS3/4A, NS5A and NS5B in a dose-dependent manner. The selected mutations permitted the development of subgenomic replicon CH3a-SGRep, by which L3004P, L3004M and Δ11nt were proven, together with a single-cycle virus production assay, to facilitate virus assembly, release, and RNA replication. CH3acc clones and CH3a-SGRep replicon provide new tools for the study of HCV genotype 3.

RNA Structures and Their Role in Selective Genome Packaging

To generate infectious viral particles, viruses must specifically select their genomic RNA from milieu that contains a complex mixture of cellular or non-genomic viral RNAs. In this review, we focus on the role of viral encoded RNA structures in genome packaging. We first discuss how packaging signals are constructed from local and long-range base pairings within viral genomes, as well as inter-molecular interactions between viral and host RNAs. Then, how genome packaging is regulated by the biophysical properties of RNA. Finally, we examine the impact of RNA packaging signals on viral evolution.

A Mathematical Analysis of HDV Genotypes: From Molecules to Cells

Hepatitis D virus (HDV) is classified according to eight genotypes. The various genotypes are included in the HDVdb database, where each HDV sequence is specified by its genotype. In this contribution, a mathematical analysis is performed on RNA sequences in HDVdb. The RNA folding predicted structures of the Genbank HDV genome sequences in HDVdb are classified according to their coarse-grain tree-graph representation. The analysis allows discarding in a simple and efficient way the vast majority of the sequences that exhibit a rod-like structure, which is important for the virus replication, to attempt to discover other biological functions by structure consideration. After the filtering, there remain only a small number of sequences that can be checked for their additional stem-loops besides the main one that is known to be responsible for virus replication. It is found that a few sequences contain an additional stem-loop that is responsible for RNA editing or other possible functions. These few sequences are grouped into two main classes, one that is well-known experimentally belonging to genotype 3 for patients from South America associated with RNA editing, and the other that is not known at present belonging to genotype 7 for patients from Cameroon. The possibility that another function besides virus replication reminiscent of the editing mechanism in HDV genotype 3 exists in HDV genotype 7 has not been explored before and is predicted by eigenvalue analysis. Finally, when comparing native and shuffled sequences, it is shown that HDV sequences belonging to all genotypes are accentuated in their mutational robustness and thermodynamic stability as compared to other viruses that were subjected to such an analysis.

Development of cell culture infectious clones for hepatitis C virus genotype 1b and transcription analysis of 1b-infected hepatoma cells

Globally, hepatitis C virus (HCV) genotype 1b is the most prevalent, and its infection has been found to associate with a higher risk of hepatocellular carcinoma (HCC) than other genotype viruses. However, an efficient infectious HCV genotype 1b culture system is unavailable, which has largely hampered the study of this important genotype virus. In this study, by using a systematic approach combining the sequences of infectious 1a TNcc clone and adaptive mutations, we succeeded in culture adaption of two full-length 1b clones for the reference strain Con1 and a clinical isolate A6, and designated as Con1cc and A6cc, respectively. Con1cc and A6cc replicated efficiently in hepatoma Huh7.5.1 cells, released HCV infectivity titers of 10^4.1 and 10^3.72 focus forming units per milliliter, respectively, and maintained the engineered mutations after passages. Both viruses responded to sofosbuvir and velpatasvir in a dose-dependent manner. With culture infectious 1b clones, we characterized the transcriptomes of 1b Con1cc-infected cells, in comparison with 2a-infected and uninfected cells. In conclusion, we have developed two infectious clones for genotype 1b and shown a novel strategy for culture adaptation of HCV isolates by using a genetically close backbone sequence. Furthermore, this study provides transcriptional landscape of HCV 1b-infected hepatoma cells facilitating the study of genotype 1b infection.

Dengue virus 2 capsid protein chaperones the strand displacement of 5'-3' cyclization sequences

By virtue of its chaperone activity, the capsid protein of dengue virus strain 2 (DENV2C) promotes nucleic acid structural rearrangements. However, the role of DENV2C during the interaction of RNA elements involved in stabilizing the 5′-3′ panhandle structure of DENV RNA is still unclear. Therefore, we determined how DENV2C affects structural functionality of the capsid-coding region hairpin element (cHP) during annealing and strand displacement of the 9-nt cyclization sequence (5CS) and its complementary 3CS. cHP has two distinct functions: a role in translation start codon selection and a role in RNA synthesis. Our results showed that cHP impedes annealing between 5CS and 3CS. Although DENV2C does not modulate structural functionality of cHP, it accelerates annealing and specifically promotes strand displacement of 3CS during 5′-3′ panhandle formation. Furthermore, DENV2C exerts its chaperone activity by favouring one of the active conformations of cHP. Based on our results, we propose mechanisms for annealing and strand displacement involving cHP. Thus, our results provide mechanistic insights into how DENV2C regulates RNA synthesis by modulating essential RNA elements in the capsid-coding region, that in turn allow for DENV replication.

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.

Identification of Potential HCV Inhibitors Based on the Interaction of Epigallocatechin-3-Gallate with Viral Envelope Proteins

Hepatitis C is affecting millions of people around the globe annually, which leads to death in very high numbers. After many years of research, hepatitis C virus (HCV) remains a serious threat to the human population and needs proper management. The in silico approach in the drug discovery process is an efficient method in identifying inhibitors for various diseases. In our study, the interaction between Epigallocatechin-3-gallate, a component of green tea, and envelope glycoprotein E2 of HCV is evaluated. Epigallocatechin-3-gallate is the most promising polyphenol approved through cell culture analysis that can inhibit the entry of HCV. Therefore, various in silico techniques have been employed to find out other potential inhibitors that can behave as EGCG. Thus, the homology modelling of E2 protein was performed. The potential lead molecules were predicted using ligand-based as well as structure-based virtual screening methods. The compounds obtained were then screened through PyRx. The drugs obtained were ranked based on their binding affinities. Furthermore, the docking of the topmost drugs was performed by AutoDock Vina, while its 2D interactions were plotted in LigPlot+. The lead compound mms02387687 (2-[[5-[(4-ethylphenoxy) methyl]-4-prop-2-enyl-1,2,4-triazol-3-yl] sulfanyl]-N-[3(trifluoromethyl) phenyl] acetamide) was ranked on top, and we believe it can serve as a drug against HCV in the future, owing to experimental validation.

Comprehensive in vivo secondary structure of the SARS-CoV-2 genome reveals novel regulatory motifs and mechanisms

Severe-acute-respiratory-syndrome-related coronavirus 2 (SARS-CoV-2) is the positive-sense RNA virus that causes coronavirus disease 2019 (COVID-19). The genome of SARS-CoV-2 is unique among viral RNAs in its vast potential to form RNA structures, yet as much as 97% of its 30 kilobases have not been structurally explored. Here, we apply a novel long amplicon strategy to determine the secondary structure of the SARS-CoV-2 RNA genome at single-nucleotide resolution in infected cells. Our in-depth structural analysis reveals networks of well-folded RNA structures throughout Orf1ab and reveals aspects of SARS-CoV-2 genome architecture that distinguish it from other RNA viruses. Evolutionary analysis shows that several features of the SARS-CoV-2 genomic structure are conserved across β-coronaviruses, and we pinpoint regions of well-folded RNA structure that merit downstream functional analysis. The native, secondary structure of SARS-CoV-2 presented here is a roadmap that will facilitate focused studies on the viral life cycle, facilitate primer design, and guide the identification of RNA drug targets against COVID-19.

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.

Ribosome Pausing at Inefficient Codons at the End of the Replicase Coding Region Is Important for Hepatitis C Virus Genome Replication

Hepatitis C virus (HCV) infects liver cells and often causes chronic infection, also leading to liver cirrhosis and cancer. In the cytoplasm, the viral structural and non-structural (NS) proteins are directly translated from the plus strand HCV RNA genome. The viral proteins NS3 to NS5B proteins constitute the replication complex that is required for RNA genome replication via a minus strand antigenome. The most C-terminal protein in the genome is the NS5B replicase, which needs to initiate antigenome RNA synthesis at the very 3′-end of the plus strand. Using ribosome profiling of cells replicating full-length infectious HCV genomes, we uncovered that ribosomes accumulate at the HCV stop codon and about 30 nucleotides upstream of it. This pausing is due to the presence of conserved rare, inefficient Wobble codons upstream of the termination site. Synonymous substitution of these inefficient codons to efficient codons has negative consequences for viral RNA replication but not for viral protein synthesis. This pausing may allow the enzymatically active replicase core to find its genuine RNA template in cis, while the protein is still held in place by being stuck with its C-terminus in the exit tunnel of the paused ribosome.

Comprehensive in-vivo secondary structure of the SARS-CoV-2 genome reveals novel regulatory motifs and mechanisms

SARS-CoV-2 is the positive-sense RNA virus that causes COVID-19, a disease that has triggered a major human health and economic crisis. The genome of SARS-CoV-2 is unique among viral RNAs in its vast potential to form stable RNA structures and yet, as much as 97% of its 30 kilobases have not been structurally explored in the context of a viral infection. Our limited knowledge of SARS-CoV-2 genomic architecture is a fundamental limitation to both our mechanistic understanding of coronavirus life cycle and the development of COVID-19 RNA-based therapeutics. Here, we apply a novel long amplicon strategy to determine for the first time the secondary structure of the SARS-CoV-2 RNA genome probed in infected cells. In addition to the conserved structural motifs at the viral termini, we report new structural features like a conformationally flexible programmed ribosomal frameshifting pseudoknot, and a host of novel RNA structures, each of which highlights the importance of studying viral structures in their native genomic context. Our in-depth structural analysis reveals extensive networks of well-folded RNA structures throughout Orf1ab and reveals new aspects of SARS-CoV-2 genome architecture that distinguish it from other single-stranded, positive-sense RNA viruses. Evolutionary analysis of RNA structures in SARS-CoV-2 shows that several features of its genomic structure are conserved across beta coronaviruses and we pinpoint individual regions of well-folded RNA structure that merit downstream functional analysis. The native, complete secondary structure of SAR-CoV-2 presented here is a roadmap that will facilitate focused studies on mechanisms of replication, translation and packaging, and guide the identification of new RNA drug targets against COVID-19.

Hepatitis C Virus Translation Regulation

Translation of the hepatitis C virus (HCV) RNA genome is regulated by the internal ribosome entry site (IRES), located in the 5’-untranslated region (5′UTR) and part of the core protein coding sequence, and by the 3′UTR. The 5′UTR has some highly conserved structural regions, while others can assume different conformations. The IRES can bind to the ribosomal 40S subunit with high affinity without any other factors. Nevertheless, IRES activity is modulated by additional cis sequences in the viral genome, including the 3′UTR and the cis-acting replication element (CRE). Canonical translation initiation factors (eIFs) are involved in HCV translation initiation, including eIF3, eIF2, eIF1A, eIF5, and eIF5B. Alternatively, under stress conditions and limited eIF2-Met-tRNA_i^Met availability, alternative initiation factors such as eIF2D, eIF2A, and eIF5B can substitute for eIF2 to allow HCV translation even when cellular mRNA translation is downregulated. In addition, several IRES trans-acting factors (ITAFs) modulate IRES activity by building large networks of RNA-protein and protein–protein interactions, also connecting 5′- and 3′-ends of the viral RNA. Moreover, some ITAFs can act as RNA chaperones that help to position the viral AUG start codon in the ribosomal 40S subunit entry channel. Finally, the liver-specific microRNA-122 (miR-122) stimulates HCV IRES-dependent translation, most likely by stabilizing a certain structure of the IRES that is required for initiation.

Hepatitis C Virus Replication

Replication and amplification of the viral genome is a key process for all viruses. For hepatitis C virus (HCV), a positive-strand RNA virus, amplification of the viral genome requires the synthesis of a negative-sense RNA template, which is in turn used for the production of new genomic RNA. This process is governed by numerous proteins, both host and viral, as well as distinct lipids and specific RNA elements within the positive- and negative-strand RNAs. Moreover, this process requires specific changes to host cell ultrastructure to create microenvironments conducive to viral replication. This review will focus on describing the processes and factors involved in facilitating or regulating HCV genome replication.

The Role of the RNA-RNA Interactome in the Hepatitis C Virus Life Cycle

RNA virus genomes are multifunctional entities endowed with conserved structural elements that control translation, replication and encapsidation, among other processes. The preservation of these structural RNA elements constraints the genomic sequence variability. The hepatitis C virus (HCV) genome is a positive, single-stranded RNA molecule with numerous conserved structural elements that manage different steps during the infection cycle. Their function is ensured by the association of protein factors, but also by the establishment of complex, active, long-range RNA-RNA interaction networks-the so-called HCV RNA interactome. This review describes the RNA genome functions mediated via RNA-RNA contacts, and revisits some canonical ideas regarding the role of functional high-order structures during the HCV infective cycle. By outlining the roles of long-range RNA-RNA interactions from translation to virion budding, and the functional domains involved, this work provides an overview of the HCV genome as a dynamic device that manages the course of viral infection.

Putative Auxin and Light Responsive Promoter Elements From the Tomato spotted wilt tospovirus Genome, When Expressed as cDNA, Are Functional in Arabidopsis

Members of the virus order Bunyavirales cause serious diseases in animals, humans and plants. Family Tospoviridae in this order contains only one genus Orthotospovirus, and members in this genus exclusively infect plants. Tomato spotted wilt tospovirus (TSWV) is considered one of the most economically important plants viruses. Little is known about the regulatory elements in the TSWV genome. Here we show that, when in the cDNA form, the 5'-upstream region of the TSWV-coded G_N/G_C gene (pG_N/G_C) possesses putative cis-regulatory elements, including an auxin responsive element (AuxRE) for binding of auxin response factors (ARFs), as well as a circadian clock-associated 1 (CCA1) protein binding site (CBS). Due to the lack of a reverse genetics system, we verified the functionality of these elements in Arabidopsis. pG_N/G_C showed light-suppressive promoter activity in transgenic Arabidopsis, and mutation in the CBS was sufficient to switch the activity to light inducible. Additionally, exogenous auxin treatments repressed the promoter activity of both wild type and CBS-mutated pG_N/G_C. Mutation in AuxRE in both promoters abolished their sensitivity to auxin. As transcriptional repressors, both CCA1 and ARF2 were able to bind to pG_N/G_C directly. To our knowledge, this is the first report that a 5'-terminal sequence of an RNA virus has light-and hormone-responsive promoter activities when expressed as cDNA in host plant's nuclear background. Our findings suggest new clues on the possible origin, evolution and function of the TSWV genomic sequence and its non-coding regions.

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).

Activation of protein kinase R by hepatitis C virus RNA-dependent RNA polymerase

Hepatitis C virus (HCV) was shown to activate protein kinase R (PKR), which inhibits expression of interferon (IFN) and IFN-stimulated genes by controlling the translation of newly transcribed mRNAs. However, it is unknown exactly how HCV activates PKR. To address the molecular mechanism(s) of PKR activation mediated by HCV infection, we examined the effects of viral proteins on PKR activation. Here, we show that expression of HCV NS5B strongly induced PKR and eIF2α phosphorylation, and attenuated MHC class I expression. In contrast, expression of Japanese encephalitis virus RNA-dependent RNA polymerase did not induce phosphorylation of PKR. Co-immunoprecipitation analyses showed that HCV NS5B interacted with PKR. Furthermore, expression of NS5B with polymerase activity-deficient mutation failed to phosphorylate PKR, suggesting that RNA polymerase activity is required for PKR activation. These results suggest that HCV activates PKR by association with NS5B, resulting in translational suppression of MHC class I to establish chronic infection.

The HCV genome domains 5BSL3.1 and 5BSL3.3 act as managers of translation

The RNA genome of the hepatitis C virus (HCV) encodes a single open reading frame (ORF) containing numerous functional elements. Among these, the cis-acting replication element (CRE) at the 3′ end of the viral ORF, has become of increasing interest given its dual role as a viral translation repressor and replication enhancer. Long-range RNA-RNA contacts mediated by the CRE build the structural scaffold required for its proper functioning. The recruitment of different cellular factors, many related to the functioning of the translation machinery, might aid in the CRE-exerted downregulation of viral translation. The present data show that the CRE promotes a defect in polysome production, and hinders the assembly of the 80S complex, likely through the direct, high affinity recruitment of the 40S ribosomal subunit. This interaction involves the highly conserved 5BSL3.1 and 5BSL3.3 domains of the CRE, and is strictly dependent on RNA-protein contacts, particularly with the ribosomal proteins RPSA and RPS29. These observations support a model in which the CRE-mediated inhibition of viral translation is a multifactorial process defined by the establishment of long-range RNA-RNA interactions between the 5′ and 3′ ends of the viral genome, the sequestration of the 40S subunit by the CRE, and the subsequent stalling of polysome elongation at the 3′ end of the ORF, all governed by the highly stable hairpin domains 5BSL3.1 and 5BSL3.3. The present data thus suggest a new managerial role in HCV translation for these 5BSL3.1 and 5BSL3.3 domains.

Form confers function: Case of the 3’X region of the hepatitis C virus genome

At the 3’ end of genomic hepatitis C virus (HCV) RNA there is a highly conserved untranslated region, the 3’X-tail, which forms part of the 3’UTR. This region plays key functions in regulation of critical processes of the viral life cycle. The 3’X region is essential for viral replication and infectivity. It is also responsible for regulation of switching between translation and transcription of the viral RNA. There is some evidence indicating the contribution of the 3’X region to the translation efficiency of the viral polyprotein and to the encapsidation process. Several different secondary structure models of the 3’X region, based on computer predictions and experimental structure probing, have been proposed. It is likely that the 3’X region adopts more than one structural form in infected cells and that a specific equilibrium between the various forms regulates several aspects of the viral life cycle. The most intriguing explanations of the structural heterogeneity problem of the 3’X region came with the discovery of its involvement in long-range RNA-RNA interactions and the potential for homodimer formation. This article summarizes current knowledge on the structure and function of the 3’X region of hepatitis C genomic RNA, reviews previous opinions, presents new hypotheses and summarizes the questions that still remain unanswered.

Signals Involved in Regulation of Hepatitis C Virus RNA Genome Translation and Replication

Hepatitis C virus (HCV) preferentially replicates in the human liver and frequently causes chronic infection, often leading to cirrhosis and liver cancer. HCV is an enveloped virus classified in the genus Hepacivirus in the family Flaviviridae and has a single-stranded RNA genome of positive orientation. The HCV RNA genome is translated and replicated in the cytoplasm. Translation is controlled by the Internal Ribosome Entry Site (IRES) in the 5' untranslated region (5' UTR), while also downstream elements like the cis-replication element (CRE) in the coding region and the 3' UTR are involved in translation regulation. The cis-elements controlling replication of the viral RNA genome are located mainly in the 5'- and 3'-UTRs at the genome ends but also in the protein coding region, and in part these signals overlap with the signals controlling RNA translation. Many long-range RNA-RNA interactions (LRIs) are predicted between different regions of the HCV RNA genome, and several such LRIs are actually involved in HCV translation and replication regulation. A number of RNA cis-elements recruit cellular RNA-binding proteins that are involved in the regulation of HCV translation and replication. In addition, the liver-specific microRNA-122 (miR-122) binds to two target sites at the 5' end of the viral RNA genome as well as to at least three additional target sites in the coding region and the 3' UTR. It is involved in the regulation of HCV RNA stability, translation and replication, thereby largely contributing to the hepatotropism of HCV. However, we are still far from completely understanding all interactions that regulate HCV RNA genome translation, stability, replication and encapsidation. In particular, many conclusions on the function of cis-elements in HCV replication have been obtained using full-length HCV genomes or near-full-length replicon systems. These include both genome ends, making it difficult to decide if a cis-element in question acts on HCV replication when physically present in the plus strand genome or in the minus strand antigenome. Therefore, it may be required to use reduced systems that selectively focus on the analysis of HCV minus strand initiation and/or plus strand initiation.

Molecular Basis of Encapsidation of Hepatitis C Virus Genome

Hepatitis C virus (HCV), a major etiologic agent of human liver diseases, is a positive-sense single-stranded RNA virus and is classified in the Flaviviridae family. Although research findings for the assembly of HCV particles are accumulating due to development of HCV cell culture system, the mechanism(s) by which the HCV genome becomes encapsidated remains largely unclear. In general, viral RNA represents only a small fraction of the RNA molecules in the cells infected with RNA viruses, but the viral genomic RNA is considered to selectively packaged into virions. It was recently demonstrated that HCV RNAs containing 3' end of the genome are selectively incorporated into virus particles during the assembly process and the 3' untranslated region functions as a cis-acting element for RNA packaging. Here, we discuss the molecular basis of RNA encapsidation of HCV and classical flaviviruses, contrast with the packaging mechanism of HIV-1.

Characterization of Recombinant Flaviviridae Viruses Possessing a Small Reporter Tag

The family Flaviviridae consists of four genera, Flavivirus, Pestivirus, Pegivirus, and Hepacivirus, and comprises important pathogens of human and animals. Although the construction of recombinant viruses carrying reporter genes encoding fluorescent and bioluminescent proteins has been reported, the stable insertion of foreign genes into viral genomes retaining infectivity remains difficult. Here, we applied the 11-amino-acid subunit derived from NanoLuc luciferase to the engineering of the Flaviviridae viruses and then examined the biological characteristics of the viruses. We successfully generated recombinant viruses carrying the split-luciferase gene, including dengue virus, Japanese encephalitis virus, hepatitis C virus (HCV), and bovine viral diarrhea virus. The stability of the viruses was confirmed by five rounds of serial passages in the respective susceptible cell lines. The propagation of the recombinant luciferase viruses in each cell line was comparable to that of the parental viruses. By using a purified counterpart luciferase protein, this split-luciferase assay can be applicable in various cell lines, even when it is difficult to transduce the counterpart gene. The efficacy of antiviral reagents against the recombinant viruses could be monitored by the reduction of luciferase expression, which was correlated with that of viral RNA, and the recombinant HCV was also useful to examine viral dynamics in vivo Taken together, our findings indicate that the recombinant Flaviviridae viruses possessing the split NanoLuc luciferase gene generated here provide powerful tools to understand viral life cycle and pathogenesis and a robust platform to develop novel antivirals against Flaviviridae viruses.IMPORTANCE The construction of reporter viruses possessing a stable transgene capable of expressing specific signals is crucial to investigations of viral life cycle and pathogenesis and the development of antivirals. However, it is difficult to maintain the stability of a large foreign gene, such as those for fluorescence and bioluminescence, after insertion into a viral genome. Here, we successfully generated recombinant Flaviviridae viruses carrying the 11-amino-acid subunit derived from NanoLuc luciferase and demonstrated that these viruses are applicable to in vitro and in vivo experiments, suggesting that these recombinant Flaviviridae viruses are powerful tools for increasing our understanding of viral life cycle and pathogenesis and that these recombinant viruses will provide a robust platform to develop antivirals against Flaviviridae viruses.

NF90-NF45 is a selective RNA chaperone that rearranges viral and cellular riboswitches: biochemical analysis of a virus host factor activity

The heterodimer NF90-NF45 is an RNA-binding protein complex that modulates the expression of various cellular mRNAs on the post-transcriptional level. Furthermore, it acts as a host factor that supports the replication of several RNA viruses. The molecular mechanisms underlying these activities have yet to be elucidated. Recently, we showed that the RNA-binding capabilities and binding specificity of NF90 considerably improves when it forms a complex with NF45. Here, we demonstrate that NF90 has a substrate-selective RNA chaperone activity (RCA) involving RNA annealing and strand displacement activities. The mechanism of the NF90-catalyzed RNA annealing was elucidated to comprise a combination of 'matchmaking' and compensation of repulsive charges, which finally results in the population of dsRNA products. Heterodimer formation with NF45 enhances 'matchmaking' of complementary ssRNAs and substantially increases the efficiency of NF90's RCA. During investigations of the relevance of the NF90-NF45 RCA, the complex was shown to stimulate the first step in the RNA replication process of hepatitis C virus (HCV) in vitro and to stabilize a regulatory element within the mRNA of vascular endothelial growth factor (VEGF) by protein-guided changes of the RNAs' structures. Thus, our study reveals how the intrinsic properties of an RNA-binding protein determine its biological activities.

Theory Meets Experiment: Metal Ion Effects in HCV Genomic RNA Kissing Complex Formation

The long-range base pairing between the 5BSL3. 2 and 3′X domains in hepatitis C virus (HCV) genomic RNA is essential for viral replication. Experimental evidence points to the critical role of metal ions, especially Mg²⁺ ions, in the formation of the 5BSL3.2:3′X kissing complex. Furthermore, NMR studies suggested an important ion-dependent conformational switch in the kissing process. However, for a long time, mechanistic understanding of the ion effects for the process has been unclear. Recently, computational modeling based on the Vfold RNA folding model and the partial charge-based tightly bound ion (PCTBI) model, in combination with the NMR data, revealed novel physical insights into the role of metal ions in the 5BSL3.2-3′X system. The use of the PCTBI model, which accounts for the ion correlation and fluctuation, gives reliable predictions for the ion-dependent electrostatic free energy landscape and ion-induced population shift of the 5BSL3.2:3′X kissing complex. Furthermore, the predicted ion binding sites offer insights about how ion-RNA interactions shift the conformational equilibrium. The integrated theory-experiment study shows that Mg²⁺ ions may be essential for HCV viral replication. Moreover, the observed Mg²⁺-dependent conformational equilibrium may be an adaptive property of the HCV genomic RNA such that the equilibrium is optimized to the intracellular Mg²⁺ concentration in liver cells for efficient viral replication.

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.

Nuclear Magnetic Resonance Study of RNA Structures at the 3'-End of the Hepatitis C Virus Genome

The 3'-end of the genomic RNA of the hepatitis C virus (HCV) embeds conserved elements that regulate viral RNA synthesis and protein translation by mechanisms that have yet to be elucidated. Previous studies with oligo-RNA fragments have led to multiple, mutually exclusive secondary structure predictions, indicating that HCV RNA structure may be context-dependent. Here we employed a nuclear magnetic resonance (NMR) approach that involves long-range adenosine interaction detection, coupled with site-specific ²H labeling, to probe the structure of the intact 3'-end of the HCV genome (385 nucleotides). Our data reveal that the 3'-end exists as an equilibrium mixture of two conformations: an open conformation in which the 98 nucleotides of the 3'-tail (3'X) form a two-stem-loop structure with the kissing-loop residues sequestered and a closed conformation in which the 3'X rearranges its structure and forms a long-range kissing-loop interaction with an upstream cis-acting element 5BSL3.2. The long-range kissing species is favored under high-Mg²⁺ conditions, and the intervening sequences do not affect the equilibrium as their secondary structures remain unchanged. The open and closed conformations are consistent with the reported function regulation of viral RNA synthesis and protein translation, respectively. Our NMR detection of these RNA conformations and the structural equilibrium in the 3'-end of the HCV genome support its roles in coordinating various steps of HCV replication.

Profiling the Mismatch Tolerance of Argonaute 2 through Deep Sequencing of Sliced Polymorphic Viral RNAs

Low allelic and clonal variability among endogenous RNAi targets has focused mismatch tolerance studies to RNAi-active guide strands. However, the inherent genomic instability of RNA viruses such as hepatitis C virus (HCV) gives rise to quasi-species mutants within discrete clones: this facilitates mismatch tolerance studies from a target perspective. We recently quantified the slicing imprecision of Argonaute 2 using small interfering RNA (siRNA) analogs of the DNA-directed RNAi drug TT-034 and next-generation sequencing of 5' RNA ligase-mediated rapid amplification of cDNA ends (RACE-SEQ). Here, we present an open-source, customizable, and computationally light RACE-SEQ bioinformatic pipeline, describing adaptations that semiquantitatively report the impact of RNAi hybridization site mismatches from the target perspective. The analysis shows that Argonaute 2 has a substitution-specific, 3- to 5-log activity window between fully complementary targets and targets with mismatches across positions 10-11. It further focuses the endonucleotic Slicer imprecision around positions 13-17, demonstrating its dependence on guide strand central region complementarity, and potentiation by even a single mismatch. We further propose pharmacogenomics value in testing endogenous targets using recombinant replicon systems and RACE-SEQ to report the pharmacodynamics of sequence-specific oligonucleotide therapeutics against all possible polymorphisms in a population, in a minimally biased, patient-free manner.

Predicting Ion Effects in an RNA Conformational Equilibrium

We develop a partial charge-based tightly bound ion (PCTBI) model for the ion effects in RNA folding. On the basis of the Monte Carlo tightly bound ion (MCTBI) approach, the model can account for ion fluctuation and correlation effects, and can predict the ion distribution around the RNA. Furthermore, unlike the previous coarse-grained RNA charge models, where negative charges are placed on the phosphates only, the current new model considers the detailed all-atom partial charge distribution on the RNA. Thus, the model not only keeps the advantage of the MCTBI model, but also has the potential to provide important detailed information unattainable by the previous MCTBI models. For example, the model predicts the reduction in ion binding upon protein binding and ion-induced conformational switches. For hepatitis C virus genomic RNA, the model predicts a Mg²⁺-induced stabilization of a kissing motif for a cis-acting regulatory element in the genomic RNA. Extensive theory-experiment comparisons support the reliability of the theoretical predictions. Therefore, the model may serve as a robust starting point for further development of an accurate method for ion effects in an RNA conformational equilibrium and RNA-cofactor interactions.

In silico methods for co-transcriptional RNA secondary structure prediction and for investigating alternative RNA structure expression

RNA transcripts are the primary products of active genes in any living organism, including many viruses. Their cellular destiny not only depends on primary sequence signals, but can also be determined by RNA structure. Recent experimental evidence shows that many transcripts can be assigned more than a single functional RNA structure throughout their cellular life and that structure formation happens co-transcriptionally, i.e. as the transcript is synthesised in the cell. Moreover, functional RNA structures are not limited to non-coding transcripts, but can also feature in coding transcripts. The picture that now emerges is that RNA structures constitute an additional layer of information that can be encoded in any RNA transcript (and on top of other layers of information such as protein-context) in order to exert a wide range of functional roles. Moreover, different encoded RNA structures can be expressed at different stages of a transcript's life in order to alter the transcript's behaviour depending on its actual cellular context. Similar to the concept of alternative splicing for protein-coding genes, where a single transcript can yield different proteins depending on cellular context, it is thus appropriate to propose the notion of alternative RNA structure expression for any given transcript. This review introduces several computational strategies that my group developed to detect different aspects of RNA structure expression in vivo. Two aspects are of particular interest to us: (1) RNA secondary structure features that emerge during co-transcriptional folding and (2) functional RNA structure features that are expressed at different times of a transcript's life and potentially mutually exclusive.

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.

The HCV Replicase Complex and Viral RNA Synthesis

Replication of hepatitis C virus (HCV) is tightly linked to membrane alterations designated the membranous web, harboring the viral replicase complex. In this chapter we describe the morphology and 3D architecture of the HCV-induced replication organelles, mainly consisting of double membrane vesicles, which are generated by a concerted action of the nonstructural proteins NS3 to NS5B. Recent studies have furthermore identified a number of host cell proteins and lipids contributing to the biogenesis of the membranous web, which are discussed in this chapter. Viral RNA synthesis is tightly associated with these membrane alterations and mainly driven by the viral RNA dependent RNA polymerase NS5B. We summarize our current knowledge of the structure and function of NS5B, the role of cis-acting replication elements at the termini of the genome in regulating RNA synthesis and the contribution of additional viral and host factors to viral RNA synthesis, which is still ill defined.

Affinity Purification of the Hepatitis C Virus Replicase Identifies Valosin-Containing Protein, a Member of the ATPases Associated with Diverse Cellular Activities Family, as an Active Virus Replication Modulator

Like almost all of the positive-strand RNA viruses, hepatitis C virus (HCV) induces host intracellular membrane modification to form the membrane-bound viral replication complex (RC), within which viral replicases amplify the viral RNA genome. Despite accumulated information about how HCV co-opts host factors for viral replication, our knowledge of the molecular mechanisms by which viral proteins hijack host factors for replicase assembly has only begun to emerge. Purification of the viral replicase and identification of the replicase-associated host factors to dissect their roles in RC biogenesis will shed light on the molecular mechanisms of RC assembly. To purify the viral replicase in the context of genuine viral replication, we developed an HCV subgenomic replicon system in which two different affinity tags were simultaneously inserted in frame into HCV NS5A and NS5B. After solubilizing the replicon cells, we purified the viral replicase by two-step affinity purification and identified the associated host factors by mass spectrometry. We identified valosin-containing protein (VCP), a member of the ATPases associated with diverse cellular activities (AAA+ATPase) family, as an active viral replication modulator whose ATPase activity is required for viral replication. A transient replication assay indicated that VCP is involved mainly in viral genome amplification. VCP associated with viral replicase and colocalized with a viral RC marker. Further, in an HCV replicase formation surrogate system, abolishing VCP function resulted in aberrant distribution of HCV NS5A. We propose that HCV may co-opt a host AAA+ATPase for its replicase assembly.

IMPORTANCE

Almost all of the positive-strand RNA viruses share a replication strategy in which viral proteins modify host membranes to form the membrane-associated viral replicase. Viruses hijack host factors to facilitate this energy-unfavorable process. Understanding of this fundamental process is hampered by the challenges of purifying the replicase because of the technical difficulties involved. In this study, we developed an HCV subgenomic replicon system in which two different affinity tags were simultaneously inserted in frame into two replicase components. Using this dual-affinity-tagged replicon system, we purified the viral replicase and identified valosin-containing protein (VCP) AAA+ATPase as a pivotal viral replicase-associated host factor that is required for viral genome replication. Abolishing VCP function resulted in aberrant viral protein distribution. We propose that HCV hijacks a host AAA+ATPase for its replicase assembly. Understanding the molecular mechanism of VCP regulates viral replicase assembly may lead to novel antiviral strategies targeting the most conserved viral replication step.

Future landscape of hepatitis C research - Basic, translational and clinical perspectives

With the latest all-oral interferon- and ribavirin-free regimens based on direct acting antivirals against the hepatitis C virus (HCV), sustained virological response rates of >90% are achieved, which is equivalent to cure. This has become possible for all genotypes and all subgroups of patients, including many of the most difficult-to-treat populations so far. Since a prophylactic HCV vaccine is not yet available, control of HCV infection will for the time being have to rely on the use of effective and safe antiviral treatments as well as their accessibility and affordability. Different approaches may apply to different parts of the world, eradication of HCV representing a major long-term goal. Whether hepatitis C becomes the first chronic viral infection to be eradicated without a prophylactic vaccine remains to be shown. Here, we briefly summarize advances in the molecular virology of hepatitis C, highlight lessons of biological relevance that were learned through the study of HCV, and its translational and clinical implications. We have also listed selected unsolved challenges, emphasizing that HCV is a unique model and that advances in this direction may yield knowledge of broad biological significance, novel technologies and insights into related important human pathogens.

Mapping the genomic diversity of HCV subtypes 1a and 1b: Implications of structural and immunological constraints for vaccine and drug development

Despite significant progress in hepatitis C (HCV) treatment, global viral eradication remains a challenge. An in-depth map of its genome diversity within the context of structural and immunological constraints could contribute to the design of pan-genotypic antivirals and preventive vaccines. For such analyses, extensive information is only available for the highly prevalent HCV genotypes (GT) 1a and 1b. Using 647 GT1a and 408 GT1b full-genome sequences obtained from the Los Alamos database, we found that respectively 3 per cent and 82 per cent of all codon positions are under positive and negative selective pressure, suggesting variation mainly accumulates due to random genetic drift. An association between conservation and both structured RNA and secondary protein structures confirmed the important role of structural elements at nucleotide and at amino acid level. Remarkably, CD8⁺ T-cell epitopes in HCV GT1a were significantly more conserved, while at the same time containing more sites under positive selection. Similarly, CD4⁺ T-cell epitopes were significantly more conserved in both HCV subtypes, but under less positive selective pressure in GT1b and more negative selective pressure in GT1a. In contrast, B-cell epitopes in both subtypes were less conserved and under less stringent negative selection. These findings argue against immune selective pressure as the main force of between-host diversifying evolution. Despite its high variability, HCV is under strict evolutionary constraints, most probably to keep its genes and proteins functional during the replication cycle. These are encouraging findings for vaccine and drug design, which could consider these newly established genetic diversity profiles.

Flexible and rapid construction of viral chimeras applied to hepatitis C virus

A novel and broadly applicable strategy combining site-directed mutagenesis and DNA assembly for constructing seamless viral chimeras is described using hepatitis C virus (HCV) as an exemplar. Full-length HCV genomic cloning cassettes, which contained flexibly situated restriction endonuclease sites, were prepared via a single, site-directed mutagenesis reaction and digested to receive PCR-amplified virus envelope genes by In-Fusion cloning. Using this method, we were able to construct gene-shuttle cassettes for generation of cell culture-infectious JFH-1-based chimeras containing genotype 1–3 E1E2 genes. Importantly, using this method we also show that E1E2 clones that were not able to support cell entry in the HCV pseudoparticle assay did confer entry when shuttled into the chimeric cell culture chimera system. This method can be easily applied to other genes of study and other viruses and, as such, will greatly simplify reverse genetics studies of variable viruses.

The cis-acting replication element of the Hepatitis C virus genome recruits host factors that influence viral replication and translation

The cis-acting replication element (CRE) of the hepatitis C virus (HCV) RNA genome is a region of conserved sequence and structure at the 3′ end of the open reading frame. It participates in a complex and dynamic RNA-RNA interaction network involving, among others, essential functional domains of the 3′ untranslated region and the internal ribosome entry site located at the 5′ terminus of the viral genome. A proper balance between all these contacts is critical for the control of viral replication and translation, and is likely dependent on host factors. Proteomic analyses identified a collection of proteins from a hepatoma cell line as CRE-interacting candidates. A large fraction of these were RNA-binding proteins sharing highly conserved RNA recognition motifs. The vast majority of these proteins were validated by bioinformatics tools that consider RNA-protein secondary structure. Further characterization of representative proteins indicated that hnRNPA1 and HMGB1 exerted negative effects on viral replication in a subgenomic HCV replication system. Furthermore DDX5 and PARP1 knockdown reduced the HCV IRES activity, suggesting an involvement of these proteins in HCV translation. The identification of all these host factors provides new clues regarding the function of the CRE during viral cycle progression.

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.

Polyprotein-Driven Formation of Two Interdependent Sets of Complexes Supporting Hepatitis C Virus Genome Replication

Hepatitis C virus (HCV) requires proteins from the NS3-NS5B polyprotein to create a replicase unit for replication of its genome. The replicase proteins form membranous compartments in cells to facilitate replication, but little is known about their functional organization within these structures. We recently reported on intragenomic replicons, bicistronic viral transcripts expressing an authentic replicase from open reading frame 2 (ORF2) and a second duplicate nonstructural (NS) polyprotein from ORF1. Using these constructs and other methods, we have assessed the polyprotein requirements for rescue of different lethal point mutations across NS3-5B. Mutations readily tractable to rescue broadly fell into two groupings: those requiring expression of a minimum NS3-5A and those requiring expression of a minimum NS3-5B polyprotein. A cis-acting mutation that blocked NS3 helicase activity, T1299A, was tolerated when introduced into either ORF within the intragenomic replicon, but unlike many other mutations required the other ORF to express a functional NS3-5B. Three mutations were identified as more refractile to rescue: one that blocked cleavage of the NS4B5A boundary (S1977P), another in the NS3 helicase (K1240N), and a third in NS4A (V1665G). Introduced into ORF1, these exhibited a dominant negative phenotype, but with K1240N inhibiting replication as a minimum NS3-5A polyprotein whereas V1665G and S1977P only impaired replication as a NS3-5B polyprotein. Furthermore, an S1977P-mutated NS3-5A polyprotein complemented other defects shown to be dependent on NS3-5A for rescue. Overall, our findings suggest the existence of two interdependent sets of protein complexes supporting RNA replication, distinguishable by the minimum polyprotein requirement needed for their formation.

IMPORTANCE Positive-strand RNA viruses reshape the intracellular membranes of cells to form a compartment within which to replicate their genome, but little is known about the functional organization of viral proteins within this structure. We have complemented protein-encoded defects in HCV by constructing subgenomic HCV transcripts capable of simultaneously expressing both a mutated and functional polyprotein precursor needed for RNA genome replication (intragenomic replicons). Our results reveal that HCV relies on two interdependent sets of protein complexes to support viral replication. They also show that the intragenomic replicon offers a unique way to study replication complex assembly, as it enables improved composite polyprotein complex formation compared to traditional trans-complementation systems. Finally, the differential behavior of distinct NS3 helicase knockout mutations hints that certain conformations of this enzyme might be particularly deleterious for replication.

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.