The 5' and 3' downstream AUG region elements are required for mosquito-borne flavivirus RNA replication

Flaviviruses require complementarity between the 5' and 3' ends of the genome for RNA replication. For mosquito-borne flaviviruses, the cyclization sequences (CS) and upstream of AUG region (UAR) elements at the genomic termini are necessary for viral RNA replication, and a third motif, the downstream of AUG region (DAR), was recently designated for dengue virus. The 3' DAR sequence is also part of a hairpin (HP-3'SL), and both complementarity between 5' and 3' DAR motifs and formation of the HP-3'SL in the absence of the 5' end are conserved among mosquito-borne flaviviruses. Using West Nile virus as a model, we demonstrate that 5'-3' DAR complementarity and HP-3'SL formation are essential for viral RNA replication.     

Replication enhancer elements within the open reading frame of tick-borne encephalitis virus and their evolution within the Flavivirus genus

We provide experimental evidence of a replication enhancer element (REE) within the capsid gene of tick-borne encephalitis virus (TBEV, genus Flavivirus). Thermodynamic and phylogenetic analyses predicted that the REE folds as a long stable stem-loop (designated SL6), conserved among all tick-borne flaviviruses (TBFV). Homologous sequences and potential base pairing were found in the corresponding regions of mosquito-borne flaviviruses, but not in more genetically distant flaviviruses. To investigate the role of SL6, nucleotide substitutions were introduced which changed a conserved hexanucleotide motif, the conformation of the terminal loop and the base-paired dsRNA stacking. Substitutions were made within a TBEV reverse genetic system and recovered mutants were compared for plaque morphology, single-step replication kinetics and cytopathic effect. The greatest phenotypic changes were observed in mutants with a destabilized stem. Point mutations in the conserved hexanucleotide motif of the terminal loop caused moderate virus attenuation. However, all mutants eventually reached the titre of wild-type virus late post-infection. Thus, although not essential for growth in tissue culture, the SL6 REE acts to up-regulate virus replication. We hypothesize that this modulatory role may be important for TBEV survival in nature, where the virus circulates by non-viraemic transmission between infected and non-infected ticks, during co-feeding on local rodents.     

Kissing-loop interaction between 5′ and 3′ ends of tick-borne Langat virus genome ‘bridges the gap’ between mosquito- and tick-borne flaviviruses in mechanisms of viral RNA cyclization: applications for virus attenuation and vaccine development

Insertion of microRNA target sequences into the flavivirus genome results in selective tissue-specific attenuation and host-range restriction of live attenuated vaccine viruses. However, previous strategies for miRNA-targeting did not incorporate a mechanism to prevent target elimination under miRNA-mediated selective pressure, restricting their use in vaccine development. To overcome this limitation, we developed a new approach for miRNA-targeting of tick-borne flavivirus (Langat virus, LGTV) in the duplicated capsid gene region (DCGR). Genetic stability of viruses with DCGR was ensured by the presence of multiple cis-acting elements within the N-terminal capsid coding region, including the stem-loop structure (5′SL6) at the 3′ end of the promoter. We found that the 5′SL6 functions as a structural scaffold for the conserved hexanucleotide motif at its tip and engages in a complementary interaction with the region present in the 3′ NCR to enhance viral RNA replication. The resulting kissing-loop interaction, common in tick-borne flaviviruses, supports a single pair of cyclization elements (CYC) and functions as a homolog of the second pair of CYC that is present in the majority of mosquito-borne flaviviruses. Placing miRNA targets into the DCGR results in superior attenuation of LGTV in the CNS and does not interfere with development of protective immunity in immunized mice.     

Long-range RNA-RNA interactions circularize the dengue virus genome

Secondary and tertiary RNA structures present in viral RNA genomes play essential regulatory roles during translation, RNA replication, and assembly of new viral particles. In the case of flaviviruses, RNA-RNA interactions between the 5' and 3' ends of the genome have been proposed to be required for RNA replication. We found that two RNA elements present at the ends of the dengue virus genome interact in vitro with high affinity. Visualization of individual molecules by atomic force microscopy revealed that physical interaction between these RNA elements results in cyclization of the viral RNA. Using RNA binding assays, we found that the putative cyclization sequences, known as 5' and 3' CS, present in all mosquito-borne flaviviruses, were necessary but not sufficient for RNA-RNA interaction. Additional sequences present at the 5' and 3' untranslated regions of the viral RNA were also required for RNA-RNA complex formation. We named these sequences 5' and 3' UAR (upstream AUG region). In order to investigate the functional role of 5'-3' UAR complementarity, these sequences were mutated either separately, to destroy base pairing, or simultaneously, to restore complementarity in the context of full-length dengue virus RNA. Nonviable viruses were recovered after transfection of dengue virus RNA carrying mutations either at the 5' or 3' UAR, while the RNA containing the compensatory mutations was able to replicate. Since sequence complementarity between the ends of the genome is required for dengue virus viability, we propose that cyclization of the RNA is a required conformation for viral replication.     

Essential Role of Cyclization Sequences in Flavivirus RNA Replication

A possible role in RNA replication for interactions between conserved complementary (cyclization) sequences in the 5'- and 3'-terminal regions of Flavivirus RNA was previously suggested but never tested in vivo. Using the M-fold program for RNA secondary-structure predictions, we examined for the first time the base-pairing interactions between the covalently linked 5' genomic region (first ~160 nucleotides) and the 3' untranslated region (last ~115 nucleotides) for a range of mosquito-borne Flavivirus species. Base-pairing occurred as predicted for the previously proposed conserved cyclization sequences. In order to obtain experimental evidence of the predicted interactions, the putative cyclization sequences (5' or 3') in the replicon RNA of the mosquito-borne Kunjin virus were mutated either separately, to destroy base-pairing, or simultaneously, to restore the complementarity. None of the RNAs with separate mutations in only the 5' or only the 3' cyclization sequences was able to replicate after transfection into BHK cells, while replicon RNA with simultaneous compensatory mutations in both cyclization sequences was replication competent. This was detected by immunofluorescence for expression of the major nonstructural protein NS3 and by Northern blot analysis for amplification and accumulation of replicon RNA. We then used the M-fold program to analyze RNA secondary structure of the covalently linked 5'- and 3'-terminal regions of three tick-borne virus species and identified a previously undescribed additional pair of conserved complementary sequences in locations similar to those of the mosquito-borne species. They base-paired with DeltaG values of approximately -20 kcal, equivalent or greater in stability than those calculated for the originally proposed cyclization sequences. The results show that the base-pairing between 5' and 3' complementary sequences, rather than the nucleotide sequence per se, is essential for the replication of mosquito-borne Kunjin virus RNA and that more than one pair of cyclization sequences might be involved in the replication of the tick-borne Flavivirus species.     

Nucleotide sequence of the genome and complete amino acid sequence of a polyprotein of the tick-borne encephalitis virus

We have cloned and sequenced RNA encoding all virion and nonstructural proteins of tick-borne encephalitis virus (TBEV). Its length is 10,477 bases with a single open reading frame (nucleotides 127-10,363) encoding 3412 amino acids. The 5'- and 3'-noncoding regions have stem- and- loop structure. The polyprotein precursor is proteolytically cleaved, apparently, by a mechanism resembling that proposed for the expression of polyproteins of other flaviviruses, such as yellow fever, West Nile and Kunjin viruses. The deduced TBEV gene order is 5'-C-preM (M)-E-NS1-NS2A-NS2B-NS3-NS4A-NS4B-NS5++ +-3'. The genome and the polyprotein of TBEV and other flaviviruses appears to be structurally similar, although these flaviviruses are transmitted to and from their vertebrate hosts by different carriers, such as ticks or mosquitoes. Analysis of sequence homologies of polyproteins of flaviviruses suggests that TBEV is more closely related to yellow fever virus than to other serological subgroups of flaviviruses (West Nile or Dengue viruses). The hydrophobic profiles of the flaviviruses are highly conservative. Nonstructural proteins NS2A, NS2B, NS4A and NS4B are extremely hydrophobic, suggesting that they are likely to be associated with cellular membranes. Proteins E, NS1, NS3 and NS5 are the most conservative and may be involved in general enzymatic activities related to viral replication and virion assembly.     

Structural and Functional Studies of the Promoter Element for Dengue Virus RNA Replication

The 5′ untranslated region (5′UTR) of the dengue virus (DENV) genome contains two defined elements essential for viral replication. At the 5′ end, a large stem-loop (SLA) structure functions as the promoter for viral polymerase activity. Next to the SLA, there is a short stem-loop that contains a cyclization sequence known as the 5′ upstream AUG region (5′UAR). Here, we analyzed the secondary structure of the SLA in solution and the structural requirements of this element for viral replication. Using infectious DENV clones, viral replicons, and in vitro polymerase assays, we defined two helical regions, a side stem-loop, a top loop, and a U bulge within SLA as crucial elements for viral replication. The determinants for SLA-polymerase recognition were found to be common in different DENV serotypes. In addition, structural elements within the SLA required for DENV RNA replication were also conserved among different mosquito- and tick-borne flavivirus genomes, suggesting possible common strategies for polymerase-promoter recognition in flaviviruses. Furthermore, a conserved oligo(U) track present downstream of the SLA was found to modulate RNA synthesis in transfected cells. In vitro polymerase assays indicated that a sequence of at least 10 residues following the SLA, upstream of the 5′UAR, was necessary for efficient RNA synthesis using the viral 3′UTR as template.     

Fine mapping of a cis-acting sequence element in yellow fever virus RNA that is required for RNA replication and cyclization

We present fine mapping of a cis-acting nucleotide sequence found in the 5' region of yellow fever virus genomic RNA that is required for RNA replication. There is evidence that this sequence interacts with a complementary sequence in the 3' region of the genome to cyclize the RNA. Replicons were constructed that had various deletions in the 5' region encoding the capsid protein and were tested for their ability to replicate. We found that a sequence of 18 nucleotides (residues 146 to 163 of the yellow fever virus genome, which encode amino acids 9 to 14 of the capsid protein) is essential for replication of the yellow fever virus replicon and that a slightly longer sequence of 21 nucleotides (residues 146 to 166, encoding amino acids 9 to 15) is required for full replication. This region is larger than the core sequence of 8 nucleotides conserved among all mosquito-borne flaviviruses and contains instead the entire sequence previously proposed to be involved in cyclization of yellow fever virus RNA.     

Identification of mutated cyclization sequences that permit efficient replication of West Nile virus genomes: use in safer propagation of a novel vaccine candidate

Existing live-attenuated flavivirus vaccines (LAV) could be improved by reducing their potential to recombine with naturally circulating viruses in the field. Since the highly conserved cyclization sequences (CS) found in the termini of flavivirus genomes must be complementary to each other to support genome replication, we set out to identify paired mutant CS that could support the efficient replication of LAV but would be unable to support replication in recombinant viruses harboring one wild-type (WT) CS. By systematic evaluation of paired mutated CS encoded in West Nile virus (WNV) replicons, we identified variants having single and double mutations in the 5'- and 3'-CS components that could support genome replication at WT levels. Replicons containing only the double-mutated CS in the 5' or the 3' ends of the genome were incapable of replication, indicating that mutated CS could be useful for constructing safer LAV. Despite the identity of the central portion of the CS in all mosquito-borne flaviviruses, viruses carrying complementary the double mutations in both the 5'- and the 3'-CS were indistinguishable from WT WNV in their replication in insect and mammalian cell lines. In addition to the utility of our novel CS pair in constructing safer LAV, we demonstrated that introduction of these mutated CS into one component of a recently described two-component genome system (A. V. Shustov, P. W. Mason, and I. Frolov, J. Virol. 81:11737-11748, 2007) enabled us to engineer a safer single-cycle WNV vaccine candidate with reduced potential for recombination during its propagation.     

Sequence comparison and secondary structure analysis of the 3' noncoding region of flavivirus genomes reveals multiple pseudoknots.

Sequences of 191 flavivirus RNAs belonging to four sero-groups were used to predict the secondary structure of the 3' noncoding region (3' NCR) directly upstream of the conserved terminal hairpin. In mosquito-borne flavivirus RNAs (n = 164) a characteristic structure element was identified that includes a phylogenetically well-supported pseudoknot. This element is repeated in the dengue and Japanese encephalitis RNAs and centers around the conserved sequences CS2 and RCS2. In yellow fever virus RNAs that contain one CS2 motif, only one copy of this pseudoknotted structure was found. The conserved pseudoknotted element is absent from the 3' NCR of tick-borne virus RNAs, which altogether adopt a secondary structure that is very different from that of mosquito-borne virus RNAs. The strong conservation of the pseudoknot in mosquito-borne flavivirus RNAs implies a stronger relationship between these viruses than concluded from previous secondary structure analyses. The role of the (tandem) pseudoknots in flavivirus replication is discussed.     

Novel cis-Acting Element within the Capsid-Coding Region Enhances Flavivirus Viral-RNA Replication by Regulating Genome Cyclization

cis-Acting elements in the viral genome RNA (vRNA) are essential for the translation, replication, and/or encapsidation of RNA viruses. In this study, a novel conserved cis-acting element was identified in the capsid-coding region of mosquito-borne flavivirus. The downstream of 5′ cyclization sequence (5′CS) pseudoknot (DCS-PK) element has a three-stem pseudoknot structure, as demonstrated by structure prediction and biochemical analysis. Using dengue virus as a model, we show that DCS-PK enhances vRNA replication and that its function depends on its secondary structure and specific primary sequence. Mutagenesis revealed that the highly conserved stem 1 and loop 2, which are involved in potential loop-helix interactions, are crucial for DCS-PK function. A predicted loop 1-stem 3 base triple interaction is important for the structural stability and function of DCS-PK. Moreover, the function of DCS-PK depends on its position relative to the 5′CS, and the presence of DCS-PK facilitates the formation of 5′-3′ RNA complexes. Taken together, our results reveal that the cis-acting element DCS-PK enhances vRNA replication by regulating genome cyclization, and DCS-PK might interplay with other cis-acting elements to form a functional vRNA cyclization domain, thus playing critical roles during the flavivirus life cycle and evolution.     

TRIM79α, an interferon-stimulated gene product, restricts tick-borne encephalitis virus replication by degrading the viral RNA polymerase

In response to virus infection, type I interferons (IFNs) induce several genes, most of whose functions are largely unknown. Here, we show that the tripartite motif (TRIM) protein, TRIM79α, is an IFN-stimulated gene (ISG) product that specifically targets tick-borne encephalitis virus (TBEV), a Flavivirus that causes encephalitides in humans. TRIM79α restricts TBEV replication by mediating lysosome-dependent degradation of the flavivirus NS5 protein, an RNA-dependent RNA polymerase essential for virus replication. NS5 degradation was specific to tick-borne flaviviruses, as TRIM79α did not recognize NS5 from West Nile virus (WNV) or inhibit WNV replication. In the absence of TRIM79α, IFN-β was less effective in inhibiting tick-borne flavivirus infection of mouse macrophages, highlighting the importance of a single virus-specific ISG in establishing an antiviral state. The specificity of TRIM79α for TBEV reveals?a remarkable ability of the innate IFN response to discriminate between closely related flaviviruses.     

Molecular evolution of the tick-borne encephalitis and Powassan viruses

The problem of emerging viruses, their genetic diversity and viral evolution in nature are attracting more attention. The phylogenetic analysis and evaluationary rate estimation were made for pathogenic flaviviruses such as tick-borne encephalitis virus (TBEV) and Powassan (PV) circulated in natural foci in Russia. 47 nucleotide sequences of encoded protein E of the TBEV and 17 sequences of NS5 genome region of the PV have been used. It was found that the rate of accumulation of nucleotide substitutions for E genome region of TBEV was approximately 1.4 x 10(-4) and 5.4 x 10(-5) substitutions per site per year for NS5 genome region of PV. The ratio of non-synonymous nucleotide substitutions to synonymous substitution (dN/dS) for viral sequences were estimated of 0.049 for TBEV and 0.098 for PV. Maximum value dN/dS was 0.201-0.220 for sub-cluster of Russian and Canadian strains of PV and the minimum - 0.024 for cluster of Russian and Chinese strains of Far Eastern genotype TBEV. Evaluation of time intervals of evolutionary events associated with these viruses showed that European subtype TBEV are diverged from all-TBEV ancestor within approximately 2750 years and the Siberian and Far Eastern subtypes are emerged about 2250 years ago. The PV was introduced into natural foci of the Primorsky Krai of Russia only about 70 years ago and PV is a very close to Canadian strains of PV. Evolutionary picture for PV in North America is similar to evolution of Siberian and Far Eastern subtypes TBEV in Asia. The divergence time for main genetic groups of TBEV and PV are correlated with historical periods of warming and cooling. These allow to propose a hypothesis that climate changes were essential to the evolution of the flaviviruses in the past millenniums.     

Molecular evolution of the tick-borne encephalitis and Powassan viruses

Issues associated with newly emerging viruses, their genetic diversity, and viral evolution in modern environments are currently attracting growing attention. In this study, a phylogenetic analysis was performed and the evolution rate was evaluated for such pathogenic flaviviruses endemic to Russia as tick-borne encephalitis virus (TBEV) and Powassan virus (PV). The analysis involved 47 nucleotide sequences of the TBEV genome region encoding protein E and 17 sequences of the PV NS5-encoding region. The nucleotide substitution rate was estimated as 1.4 × 10⁻⁴ and 5.4 × 10⁻⁵ substitutions per site per year for the E protein-encoding region of the TBEV genome and for the NS5 genome region of PV, respectively. The ratio of non-synonymous to synonymous nucleotide substitutions (dN/dS) in viral sequences was calculated as 0.049 for TBEV and 0.098 for PV. The highest dN/dS values of 0.201–0.220 were found in the subcluster of Russian and Canadian PV strains, and the lowest value of 0.024 was observed in the cluster of Russian and Chinese strains of the Far Eastern TBEV genotype. Evaluation of time intervals between the events of viral evolution showed that the European subtype of TBEV diverged from the common TBEV ancestor approximately 2750 years ago, while the Siberian and Far Eastern subtypes emerged approximately 2250 years ago. The PV was introduced into its natural foci of the Russian Primorskii krai only approximately 70 years ago; these strains were very close to Canadian PV strains. The pattern of PV evolution in North America was similar to the evolution of the Siberian and Far Eastern TBEV subtypes in Asia. The moments of divergence between major genetic groups of TBEV and PV coincide with historical periods of climate warming and cooling, suggesting that climate change was a key factor in the evolution of flaviviruses in past millennia.     

A tick-borne Langat virus mutant that is temperature sensitive and host range restricted in neuroblastoma cells and lacks neuroinvasiveness for immunodeficient mice

Langat virus (LGT), the naturally attenuated member of the tick-borne encephalitis virus (TBEV) complex, was tested extensively in clinical trials as a live TBEV vaccine and was found to induce a protective, durable immune response; however, it retained a low residual neuroinvasiveness in mice and humans. In order to ablate or reduce this property, LGT mutants that produced a small plaque size or temperature-sensitive (ts) phenotype in Vero cells were generated using 5-fluorouracil. One of these ts mutants, clone E5-104, exhibited a more than 10(3)-fold reduction in replication at the permissive temperature in both mouse and human neuroblastoma cells and lacked detectable neuroinvasiveness for highly sensitive immunodeficient mice. The E5-104 mutant possessed five amino acid substitutions in the structural protein E and one change in each of the nonstructural proteins NS3 and NS5. Using reverse genetics, we demonstrated that a Lys(46)-->Glu substitution in NS3 as well as a single Lys(315)-->Glu change in E significantly impaired the growth of LGT in neuroblastoma cells and reduced its peripheral neurovirulence for SCID mice. This study and our previous experience with chimeric flaviviruses indicated that a decrease in viral replication in neuroblastoma cells might serve as a predictor of in vivo attenuation of the neurotropic flaviviruses. The combination of seven mutations identified in the nonneuroinvasive E5-104 mutant provided a useful foundation for further development of a live attenuated TBEV vaccine. An evaluation of the complete sequence of virus recovered from brain of SCID mice inoculated with LGT mutants identified sites in the LGT genome that promoted neurovirulence/neuroinvasiveness.     

Induction and suppression of tick cell antiviral RNAi responses by tick-borne flaviviruses

Arboviruses are transmitted by distantly related arthropod vectors such as mosquitoes (class Insecta) and ticks (class Arachnida). RNA interference (RNAi) is the major antiviral mechanism in arthropods against arboviruses. Unlike in mosquitoes, tick antiviral RNAi is not understood, although this information is important to compare arbovirus/host interactions in different classes of arbovirus vectos. Using an Ixodes scapularis-derived cell line, key Argonaute proteins involved in RNAi and the response against tick-borne Langat virus (Flaviviridae) replication were identified and phylogenetic relationships characterized. Analysis of small RNAs in infected cells showed the production of virus-derived small interfering RNAs (viRNAs), which are key molecules of the antiviral RNAi response. Importantly, viRNAs were longer (22 nucleotides) than those from other arbovirus vectors and mapped at highest frequency to the termini of the viral genome, as opposed to mosquito-borne flaviviruses. Moreover, tick-borne flaviviruses expressed subgenomic flavivirus RNAs that interfere with tick RNAi. Our results characterize the antiviral RNAi response in tick cells including phylogenetic analysis of genes encoding antiviral proteins, and viral interference with this pathway. This shows important differences in antiviral RNAi between the two major classes of arbovirus vectors, and our data broadens our understanding of arthropod antiviral RNAi.     

Solution structure and structural dynamics of envelope protein domain III of mosquito- and tick-borne flaviviruses

The mosquito-borne West Nile (WNV) and dengue 2 (DEN2V) viruses and tick-borne Langat (LGTV) and Omsk hemorrhagic fever (OHFV) viruses are arthropod-borne flaviviruses (family Flaviviridae, genus Flavivirus). These viruses are quite similar at both the nucleotide and amino acid level, yet they are very divergent in their biological properties and in the diseases they cause. The objective of this study was to examine the putative receptor-binding domains of the flaviviruses, the envelope (E) protein domain III (D3), which assume very similar structures either as part of the whole envelope protein or as individual entities, and to define the biophysical properties that distinguish among these viruses. Circular dichroism and Fourier transform infrared spectroscopy were employed to monitor the solution structure of these proteins. While the spectroscopic results found that the D3 from each of these viruses is composed of either beta-sheets or beta-turns, which is consistent with X-ray crystal data for tick-borne encephalitis and dengue viruses, these results reveal that recombinant D3s (rED3s) derived from tick-borne flaviviruses (LGT-rED3 and OHF-rED3) were similar to each other, while those from mosquito-borne flaviviruses (WN-rED3 and DEN-rED3) were similar to each other yet distinct from rED3 of the tick-borne viruses. Protein dynamic studies probed by fluorescence quenching and hydrogen/deuterium exchange found that the rED3s are dynamic entities. The tick-borne proteins again exhibit very similar dynamic properties, which are different from the mosquito-borne proteins. The WN-rED3 is significantly less stable than the other three rED3s. Overall, these differences in biophysical properties correlate with biological properties of these viruses that tick-borne flaviviruses are more stable than mosquito-borne flaviviruses.     

RNA secondary structure in the coding region of dengue virus type 2 directs translation start codon selection and is required for viral replication

Dengue virus is a positive-strand RNA virus and a member of the genus Flavivirus, which includes West Nile, yellow fever, and tick-borne encephalitis viruses. Flavivirus genomes are translated as a single polyprotein that is subsequently cleaved into 10 proteins, the first of which is the viral capsid (C) protein. Dengue virus type 2 (DENV2) and other mosquito-borne flaviviruses initiate translation of C from a start codon in a suboptimal context and have multiple in-frame AUGs downstream. Here, we show that an RNA hairpin structure in the capsid coding region (cHP) directs translation start site selection in human and mosquito cells. The ability of the cHP to direct initiation from the first start codon is proportional to its thermodynamic stability, is position dependent, and is sequence independent, consistent with a mechanism in which the scanning initiation complex stalls momentarily over the first AUG as it begins to unwind the cHP. The cHP of tick-borne flaviviruses is not maintained in a position to influence start codon selection, which suggests that this coding region cis element may serve another function in the flavivirus life cycle. Here, we demonstrate that the DENV2 cHP and both the first and second AUGs of C are necessary for efficient viral replication in human and mosquito cells. While numerous regulatory elements have been identified in the untranslated regions of RNA viral genomes, we show that the cHP is a coding-region RNA element that directs start codon selection and is required for viral replication.