An RNA tertiary structure in the 3' untranslated region of enteroviruses is necessary for efficient replication

The 3' end of the coding region of senecavirus A contains a highly conserved sequence that potentially forms a stem-loop structure required for virus rescue

Senecavirus A (SVA) can cause a vesicular disease in swine. It is a positive-strand RNA virus belonging to the genus Senecavirus in the family Picornaviridae. Positive-strand RNA viruses possess positive-sense, single-stranded genomes whose untranslated regions (UTRs) have been reported to contain cis-acting RNA elements. In the present study, a total of 100 SVA isolates were comparatively analyzed at the genome level. A highly conserved fragment (HCF) was found to be located in the 3D sequence and to be close to the 3' UTR. The HCF was computationally predicted to form a stem-loop structure. Eight synonymous mutations can individually disrupt the formation of a single base pair within the stem region. We found that SVA itself was able to tolerate each of these mutations alone, as evidenced by the ability to rescue all eight single-site mutants from their individual cDNA clones, and all of them were genetically stable during serial passaging. However, the replication-competent SVA could not be rescued from another cDNA clone containing all eight mutations. The failure to recover SVA might be attributed to disruption of the predicted stem-loop structure, whereas introduction of a wild-type HCF into the cDNA clone with eight mutations still had no effect on virus recovery. These results suggest that the putative stem-loop structure at the 3' end of the 3D sequence is a cis-acting RNA element that is required for SVA growth.

A putative wild-type or wild-type-like hairpin structure is required within 3' untranslated region of Senecavirus A for virus replication

The 3' untranslated region (UTR) of Senecavirus A (SVA) was predicted to harbor two hairpin structures, hairpin-I and -II. The former is composed of two internal loops, one terminal loop and three stem regions; the latter comprises one internal loop, one terminal loop and two stem regions. In this study, we constructed a total of nine SVA cDNA clones, which contained different point mutations within a stem-formed motif in the hairpin-I or -II, for rescuing replication-competent viruses. Only three mutants were successfully rescued and moreover genetically stable during at least five serial passages. Computer-aided prediction showed these three mutants bearing either a wild-type or a wild-type-like hairpin-I in their individual 3' UTRs. Neither wild-type nor wild-type-like hairpin-I could be computationally predicted to exist in 3' UTRs of the other six unviable "viruses". The results suggested that the wild-type or wild-type-like hairpin-I was necessary in the 3' UTR for SVA replication.

Tolerance of Senecavirus A to Mutations in Its Kissing-Loop or Pseudoknot Structure Computationally Predicted in 3′ Untranslated Region

Senecavirus A (SVA) is an emerging virus that belongs to the genus Senecavirus in the family Picornaviridae. Its genome is a positive-sense and single-stranded RNA, containing two untranslated regions (UTRs). The 68-nt-long 3′ UTR is computationally predicted to possess two higher-order RNA structures: a kissing-loop interaction and an H-type-like pseudoknot, both of which, however, cannot coexist in the 3′ UTR. In this study, we constructed 17 full-length SVA cDNA clones (cD-1 to -17): the cD-1 to -7 contained different point mutations in a kissing-loop-forming motif (KLFM); the cD-8 to -17 harbored one single or multiple point mutations in a pseudoknot-forming motif (PFM). These 17 mutated cDNA clones were independently transfected into BSR-T7/5 cells for rescuing recombinant SVAs (rSVAs), named rSVA-1 to −17, corresponding to cD-1 to −17. The results showed that the rSVA-1, -2, -3, -4, -5, -6, -7, -9, -13, and -15 were successfully rescued from their individual cDNA clones. Moreover, all mutated motifs were genetically stable during 10 viral passages in vitro. This study unveiled viral abilities of tolerating mutations in the computationally predicted KLFM or PFMs. It can be concluded that the putative kissing-loop structure, even if present in the 3′ UTR, is unnecessary for SVA replication. Alternatively, if the pseudoknot formation potentially occurs in the 3′ UTR, its deformation would have a lethal effect on SVA propagation.

Oncolytic Activity of Targeted Picornaviruses Formulated as Synthetic Infectious RNA

Infectious nucleic acid has been proposed as a superior formulation for oncolytic virus therapy. Oncolytic picornaviruses can be formulated as infectious RNA (iRNA), and their unwanted tropisms eliminated by microRNA (miRNA) detargeting. However, genomic insertion of miRNA target sequences into coxsackievirus A21 (CVA21) iRNA compromised its specific infectivity, negating further development as a novel oncolytic virus formulation. To address this limitation, we substituted a muscle-specific miRNA response element for the spacer region downstream of the internal ribosomal entry site in the 5′ non-coding region of CVA21 iRNA, thereby preserving genome length while avoiding the disruption of known surrounding RNA structural elements. This new iRNA (R-CVA21) retained high specific infectivity, rapidly generating replicating miRNA-detargeted viruses following transfection in H1-HeLa cells. Further, in contrast with alternatively configured iRNAs that were tested in parallel, intratumoral administration of R-CVA21 generated a spreading oncolytic infection that was curative in treated animals without associated myotoxicity. Moreover, R-CVA21 also exhibited superior miRNA response element stability in vivo. This novel formulation is a promising agent for clinical translation.

The Roles of Picornavirus Untranslated Regions in Infection and Innate Immunity

Viral genomes have evolved to maximize their potential of overcoming host defense mechanisms and to induce a variety of disease syndromes. Structurally, a genome of a virus consists of coding and noncoding regions, and both have been shown to contribute to initiation and progression of disease. Accumulated work in picornaviruses has stressed out the importance of the noncoding RNAs, or untranslated 5′- and 3′-regions (UTRs), in both replication and translation of viral genomes. Unsurprisingly, defects in these processes have been reported to cause viral attenuation and affect viral pathogenicity. However, substantial evidence suggests that these untranslated RNAs may influence the outcome of the host innate immune response. This review discusses the involvement of 5′- and 3′-terminus UTRs in induction and regulation of host immunity and its consequences for viral life cycle and virulence.

Emergency Services of Viral RNAs: Repair and Remodeling

Reproduction of RNA viruses is typically error-prone due to the infidelity of their replicative machinery and the usual lack of proofreading mechanisms. The error rates may be close to those that kill the virus. Consequently, populations of RNA viruses are represented by heterogeneous sets of genomes with various levels of fitness. This is especially consequential when viruses encounter various bottlenecks and new infections are initiated by a single or few deviating genomes. Nevertheless, RNA viruses are able to maintain their identity by conservation of major functional elements. This conservatism stems from genetic robustness or mutational tolerance, which is largely due to the functional degeneracy of many protein and RNA elements as well as to negative selection. Another relevant mechanism is the capacity to restore fitness after genetic damages, also based on replicative infidelity. Conversely, error-prone replication is a major tool that ensures viral evolvability. The potential for changes in debilitated genomes is much higher in small populations, because in the absence of stronger competitors low-fit genomes have a choice of various trajectories to wander along fitness landscapes. Thus, low-fit populations are inherently unstable, and it may be said that to run ahead it is useful to stumble. In this report, focusing on picornaviruses and also considering data from other RNA viruses, we review the biological relevance and mechanisms of various alterations of viral RNA genomes as well as pathways and mechanisms of rehabilitation after loss of fitness. The relationships among mutational robustness, resilience, and evolvability of viral RNA genomes are discussed.

Complete genome characterization of a novel enterovirus type EV-B106 isolated in China, 2012

Human enterovirus B106 (EV-B106) is a recently identified member of enterovirus species B. In this study, we report the complete genomic characterization of an EV-B106 strain (148/YN/CHN/12) isolated from an acute flaccid paralysis patient in Yunnan Province, China. The new strain had 79.2–81.3% nucleotide and 89.1–94.8% amino acid similarity in the VP1 region with the other two EV-B106 strains from Bolivia and Pakistan. When compared with other EV serotypes, it had the highest (73.3%) VP1 nucleotide similarity with the EV-B77 prototype strain CF496-99. However, when aligned with all EV-B106 and EV-B77 sequences available from the GenBank database, two major frame shifts were observed in the VP1 coding region, which resulted in substantial (20.5%) VP1 amino acid divergence between the two serotypes. Phylogenetic analysis and similarity plot analysis revealed multiple recombination events in the genome of this strain. This is the first report of the complete genome of EV-B106.

Comparative genomic analysis of coxsackievirus A6 strains of different clinical disease entities

Background

Studies regarding coxsackievirus A6 (CVA6) infection were limited. In Taiwan, outbreaks of CVA6 occurred in 2009 and 2010, respectively, but the clinical manifestations were markedly different. We conducted a study to compare the clinical features and genomic sequence between the two years.

Methodology/Principal Findings

In 2009 and 2010, 205 patients with coxsackievirus A6 (CVA6) infection were treated at Chang Gung Memorial Hospital. Detailed clinical features were obtained from 126 inpatients, 62 in 2009 and 64 in 2010. Between the inpatients in 2009 and 2010, no statistically significant difference was noted in terms of demographics, length of hospital stay and laboratory data. Significantly more patients in 2009 presented with herpangina (82%) while more patients in 2010 presented with hand-foot-mouth disease (HFMD; 67%) and skin rash beyond the typical sites for HFMD. Complete genomic sequences were determined and compared for three isolates from patients with herpangina in 2009 and three isolates from patients with HFMD in 2010. The complete sequences showed that 2009 and 2010 CVA6 isolates were indistinguishable by partial VP1 genes, but there were 5 unique nucleotide changes in 3′ UTR, and 23 out of 2201 (1%) amino acids were different. 2010 viruses underwent the largest number of amino acid changes in 3CD protein, which is the precursor of both 3C protease and 3D polymerase.

Conclusions

Since 2008 in Finland, outbreaks of HFMD due to CVA6 were noted internationally. CVA6 of different genetic background may cause different clinical manifestations such as herpangina and HFMD.

Specific temperature-induced perturbations of secondary mRNA structures are associated with the cold-adapted temperature-sensitive phenotype of influenza A virus

For decades, cold-adapted, temperature-sensitive (ca/ts) strains of influenza A virus have been used as live attenuated vaccines. Due to their great public health importance it is crucial to understand the molecular mechanism(s) of cold adaptation and temperature sensitivity that are currently unknown. For instance, secondary RNA structures play important roles in influenza biology. Thus, we hypothesized that a relatively minor change in temperature (32-39°C) can lead to perturbations in influenza RNA structures and, that these structural perturbations may be different for mRNAs of the wild type (wt) and ca/ts strains. To test this hypothesis, we developed a novel in silico method that enables assessing whether two related RNA molecules would undergo (dis)similar structural perturbations upon temperature change. The proposed method allows identifying those areas within an RNA chain where dissimilarities of RNA secondary structures at two different temperatures are particularly pronounced, without knowing particular RNA shapes at either temperature. We identified such areas in the NS2, PA, PB2 and NP mRNAs. However, these areas are not identical for the wt and ca/ts mutants. Differences in temperature-induced structural changes of wt and ca/ts mRNA structures may constitute a yet unappreciated molecular mechanism of the cold adaptation/temperature sensitivity phenomena.

Non-template functions of viral RNA in picornavirus replication

The genomic RNA of poliovirus and closely related picornaviruses perform template and non-template functions during viral RNA replication. The non-template functions are mediated by cis-active RNA sequences that bind viral and cellular proteins to form RNP complexes. The RNP complexes mediate temporally dynamic, long-range interactions in the viral genome and ensure the specificity of replication. The 5' cloverleaf (5' CL)-RNP complex serves as a key cis-active element in all of the non-template functions of viral RNA. The 5'CL-RNP complex is proposed to interact with the cre-RNP complex during VPgpUpU synthesis, the 3'NTR-poly(A) RNP complex during negative-strand initiation and the 30 end negative-strand-RNP complex during positive-strand initiation. Co-ordinating these long-range interactions is important in regulating each step in the replication cycle.

Genetic characterization of enterovirus 71 isolated from patients with severe disease by comparative analysis of complete genomes

Enterovirus 71 (EV71) which causes mild illness in children is also associated with severe neurological complications. This study analyzed the complete genomes of EV71 strains derived from mild and severe diseases in order to determine whether the differences of EV71 genomes were responsible for different clinical presentations. Compared to complete genomes of EV71 strains derived from mild cases (less virulent strains), nucleotide differences in EV71 strains isolated from severe cases (more virulent strains) were observed primarily in the internal ribosomal entry site (IRES) of the 5'-untranslated region (UTR), which is vital for the cap-independent translation of viral proteins. In the protein-coding region, an E-Q substitution at amino acid position 145 of structural protein VP1 that occurred in more than one of more virulent strains was observed. This site is known to be related functionally to receptor binding and virulence in mice. Overall, strains (Group III) isolated from patients with fatal or severe sequelae outcomes had greater sequence substitutions in the 5'-UTR and/or protein-coding region and exhibited a relatively low-average homology to less virulent strains across the entire genome, indicating the possibility of significant genomic diversity in the most virulent EV71 strains. Further studies of EV71 pathogenesis should examine the significance of genomic diversity and the effects of multiple mutations in a viral population.

Evolution of newly described enteroviruses

Several new enterovirus serotypes and a new human rhinovirus species have been characterized in the Enterovirus genus recently, raising a question about the origin of the new viruses. In this article we attempt to outline the general patterns of enterovirus evolution, ultimately leading to the emergence of new serotypes or species. Different evolutionary and epidemiological patterns can be deduced between different enterovirus species, between entero- and rhino-viruses and between different serotypes within a species. This article presents a hypothesis that the divergent evolution leading to a new serotype is likely to involve adaptation to a new ecological niche either within a single host species or due to interspecies transmission. By contrast, evolution within a serotype appears to occur primarily by genetic drift.

Heuristic RNA pseudoknot prediction including intramolecular kissing hairpins

Pseudoknots are an essential feature of RNA tertiary structures. Simple H-type pseudoknots have been studied extensively in terms of biological functions, computational prediction, and energy models. Intramolecular kissing hairpins are a more complex and biologically important type of pseudoknot in which two hairpin loops form base pairs. They are hard to predict using free energy minimization due to high computational requirements. Heuristic methods that allow arbitrary pseudoknots strongly depend on the quality of energy parameters, which are not yet available for complex pseudoknots. We present an extension of the heuristic pseudoknot prediction algorithm DotKnot, which covers H-type pseudoknots and intramolecular kissing hairpins. Our framework allows for easy integration of advanced H-type pseudoknot energy models. For a test set of RNA sequences containing kissing hairpins and other types of pseudoknot structures, DotKnot outperforms competing methods from the literature. DotKnot is available as a web server under http://dotknot.csse.uwa.edu.au.

Genomic characterization of an enterovirus 97 strain isolated in Shandong, China

The genomic characterization of human enterovirus 97 (EV97) strain isolated from an acute flaccid paralysis case in Shandong province, China in 1999, is described. The strain, designated as 99188/SD/CHN/1999/EV97 (abbreviated as 99188), had a genome of 7394 nucleotides. Compared with other EV97 strains, it had 81.3-83.3% nucleotide similarity and 94.0-95.4% amino acid similarity in VP1 coding region, and it had 81.4% complete genomic similarity with prototype strain BAN99-10355. The most striking feature was the deletion of 18 nucleotides in the 3' end of VP1 coding region, combined with two deletions and one insertion in 5' and 3' untranslated regions. All these findings demonstrated the strain 99188 had a distant genetic relationship with other EV97 strains. In the phylogenetic trees generated from VP1 and 3D sequences of human enterovirus species B (HEV-B), the lineages of strain 99188 were not congruent, suggesting the event of recombination. Similarity plot analysis further provided the evidence of recombination with other strains of HEV-B in P2 and P3 coding region. This is the first finding of EV97 in China and the third genomic sequence of EV97 reported.

Molecular characterization of a Chinese vaccine strain of transmissible gastroenteritis virus: mutations that may contribute to attenuation

Transmissible gastroenteritis virus (TGEV), the etiological agent of transmissible gastroenteritis (TGE), is the major cause of viral enteritis and fetal diarrhea in swine neonates, resulting in significant economic losses to the swine industry. The Chinese vaccine strain H165 of TGEV was derived from a virulent field strain H16 by serial passage in vitro. Strain H165 has been proven to be safe in piglets and pregnant sows and displays efficacy against TGEV infection. In this study, we report the complete genome sequences of strains H165 and H16, obtained by sequencing several overlapping fragments amplified from viral RNA and our findings from sequence and phylogenetic analyses. The genomes were 28,569 nucleotides in length, including the poly (A) tail. No deletions or insertions were detected in the H16 genome sequence after continuous passage in vitro; however, we found 27 nucleotide mutations in strain H165 compared with strain H16, resulting in 16 amino acid changes distributed among the genes 1, S, 3, and sM. An A to G nucleotide mutation was found in the intergenic region between the 3a and 3b genes. Furthermore, six unique nucleotides identified in the genome sequence of H165 could be used as makers to differentiate the H165 vaccine strain from wild-type TGEV strains. Our findings from phylogenetic analysis may enhance our understanding of the evolution of TGEV, as well as the other coronaviruses.

The structure-function relationship of the enterovirus 3'-UTR

Essential processes in living cells are carried out by large complex assemblies, which typically consist of a large number of proteins and frequently also contain nucleic acids, mostly RNA [Alberts, B., 1998. The cell as a collection of protein machines: preparing the next generation of molecular biologists. Cell 92, 291-294]. These large biomolecular complexes carry out biological processes in highly sophisticated ways: molecules do not move around randomly in the cell and interact by chance, but are guided to these "macromolecular machines", in which the number of possible collisions is restricted to a few possibilities, based, e.g., on the specificity of protein-RNA recognition. While the coding capacity of RNA lies within its sequence, the shape of an RNA molecule determines other functionalities such as stability, intra- and intermolecular interactions, catalytic activity, regulation of cellular processes, etc. [Doudna, J.A., 2000. Structural genomics of RNA. Nat. Struct. Biol. 7, 954-956; Cech, T.R. 2000. Structural biology. The ribosome is a ribozyme. Science 289, 878-879]. RNA structures in macromolecular machines are important features in assembly, target recognition and activity. Viral RNA molecules contain cis- and/or trans-acting control elements that, as exemplified by internal ribosomal entry sites and origins of genome replication, consist of complex multidomain structures [Andino, R., Rieckhof, G.E., Achacoso, P.L., Baltimore D., 1993. Poliovirus RNA synthesis utilizes an RNP complex formed around the 5'-end of viral RNA. EMBO J. 12, 3587-3598; Melchers, W.J.G., Hoenderop, J.G.J., Bruins Slot, H.J., Pleij, C.W.A., Pilipenko, E.V., Agol, V.I., Galama, J.M.D., 1997. Kissing of the two predominant hairpin loops in the coxsackie B virus 3' untranslated region is the essential structural feature of the origin of replication required for negative-strand RNA synthesis. J. Virol. 71, 686-696]. The formation of these structures is involved in the specific recognition of ligands or serves to support the structural integrity of the whole element. The replication of the enterovirus RNA is carried out by a large biomolecular complex formed by cis-acting RNA elements found in the 5'- and 3'-UTR of the virus genome and several cellular and viral proteins. This review will focus on RNA elements in the 3'-UTR of enteroviruses.

Complete genome sequence analysis of Seneca Valley virus-001, a novel oncolytic picornavirus

The complete genome sequence of Seneca Valley virus-001 (SVV-001), a small RNA virus, was determined and was shown to have typical picornavirus features. The 7280 nt long genome was predicted to contain a 5' untranslated region (UTR) of 666 nt, followed by a single long open reading frame consisting of 6543 nt, which encodes a 2181 aa polyprotein. This polyprotein could potentially be cleaved into 12 polypeptides in the standard picornavirus L-4-3-4 layout. A 3' UTR of 71 nt was followed by a poly(A) tail of unknown length. Comparisons with other picornaviruses showed that the P1, 2C, 3C and 3D polypeptides of SVV-001 were related most closely to those of the cardioviruses, although they were not related as closely to those of encephalomyocarditis virus and Theiler's murine encephalomyelitis virus as the latter were to each other. Most other regions of the polyprotein differed considerably from those of all other known picornaviruses. SVV-001 contains elements of an internal ribosome entry site reminiscent of that found in hepatitis C virus and a number of genetically diverse picornaviruses. SVV-001 is a novel picornavirus and it is proposed that it be classified as the prototype species in a novel genus named 'Senecavirus'.

Bioinformatic and functional analysis of RNA secondary structure elements among different genera of human and animal caliciviruses

The mechanism and role of RNA structure elements in the replication and translation of Caliciviridae remains poorly understood. Several algorithmically independent methods were used to predict secondary structures within the Norovirus, Sapovirus, Vesivirus and Lagovirus genera. All showed profound suppression of synonymous site variability (SSSV) at genomic 5' ends and the start of the sub-genomic (sg) transcript, consistent with evolutionary constraints from underlying RNA structure. A newly developed thermodynamic scanning method predicted RNA folding mapping precisely to regions of SSSV and at the genomic 3' end. These regions contained several evolutionarily conserved RNA secondary structures, of variable size and positions. However, all caliciviruses contained 3' terminal hairpins, and stem-loops in the anti-genomic strand invariably six bases upstream of the sg transcript, indicating putative roles as sg promoters. Using the murine norovirus (MNV) reverse-genetics system, disruption of 5' end stem-loops produced approximately 15- to 20-fold infectivity reductions, while disruption of the RNA structure in the sg promoter region and at the 3' end entirely destroyed replication ability. Restoration of infectivity by repair mutations in the sg promoter region confirmed a functional role for the RNA secondary structure, not the sequence. This study provides comprehensive bioinformatic resources for future functional studies of MNV and other caliciviruses.

Identification of a contemporary human parechovirus type 1 by VIDISCA and characterisation of its full genome

Background

Enteritis is caused by a spectrum of viruses that is most likely not fully characterised. When testing stool samples by cell culture, virus isolates are sometimes obtained which cannot be typed by current methods. In this study we used VIDISCA, a virus identification method which has not yet been widely applied, on such an untyped virus isolate.

Results

We found a human parechovirus (HPeV) type 1 (strain designation: BNI-788st). Because genomes of contemporary HPeV1 were not available, we determined its complete genome sequence. We found that the novel strain was likely the result of recombination between structural protein genes of an ancestor of contemporary HPeV1 strains and nonstructural protein genes from an unknown ancestor, most closely related to HPeV3. In contrast to the non-structural protein genes of other HPeV prototype strains, the non-structural protein genes of BNI-788st and HPeV3 prototype strains did not co-segregate in bootscan analysis with that of other prototype strains.

Conclusion

HPeV3 nonstructural protein genes may form a distinct element in a pool of circulating HPeV non-structural protein genes. More research into the complex HPeV evolution is required to connect virus ecology with disease patterns in humans.

Complete genome sequences for nine simian enteroviruses

Analysis of the VP1 capsid-coding sequences of the simian picornaviruses has suggested that baboon enterovirus (BaEV), SV19, SV43 and SV46 belong to the species Human enterovirus A (HEV-A) and SA5 belongs to HEV-B, whereas SV4/A2 plaque virus (two isolates of a single serotype), SV6 and N125/N203 (two isolates of a single serotype) appear to represent new species in the genus. We have further characterized by complete genomic sequencing the genetic relationships among the simian enteroviruses serotypes (BaEV, N125/N203, SA5, SV4/A2 plaque virus, SV6, SV19, SV43 and SV46) and to other enteroviruses. Phylogenetic and pairwise sequence relationships for the P1 region paralleled those of VP1 alone, and confirmed that SV4/A-2 plaque virus, SV6 and N125/N203 represent unique genetic clusters that probably correspond to three new species. However, sequence relationships in the P2 and P3 regions were quite different. In 2C, SV19, SV43 and SV46 remain clustered with the human viruses of HEV-A, but BaEV, SV6 and N125/N203 cluster together; in 3CD, SA5 (HEV-B) also joined this cluster. The 3'-non-translated region (NTR) sequences are highly conserved within each of the four human enterovirus species, but the 3'-NTRs of the simian enteroviruses are distinct from those of all human enteroviruses and generally distinct from one another. These results suggest that host species may have a significant influence on the evolution of enterovirus non-capsid sequences.

An RNA pseudoknot in the 3' end of the arterivirus genome has a critical role in regulating viral RNA synthesis

In the life cycle of plus-strand RNA viruses, the genome initially serves as the template for both translation of the viral replicase gene and synthesis of minus-strand RNA and is ultimately packaged into progeny virions. These various processes must be properly balanced to ensure efficient viral proliferation. To achieve this, higher-order RNA structures near the termini of a variety of RNA virus genomes are thought to play a key role in regulating the specificity and efficiency of viral RNA synthesis. In this study, we have analyzed the signals for minus-strand RNA synthesis in the prototype of the arterivirus family, equine arteritis virus (EAV). Using site-directed mutagenesis and an EAV reverse genetics system, we have demonstrated that a stem-loop structure near the 3' terminus of the EAV genome is required for RNA synthesis. We have also obtained evidence for an essential pseudoknot interaction between the loop region of this stem-loop structure and an upstream hairpin residing in the gene encoding the nucleocapsid protein. We propose that the formation of this pseudoknot interaction may constitute a molecular switch that could regulate the specificity or timing of viral RNA synthesis. This hypothesis is supported by the fact that phylogenetic analysis predicted the formation of similar pseudoknot interactions near the 3' end of all known arterivirus genomes, suggesting that this interaction has been conserved in evolution.

Specific interaction of HeLa cell proteins with coxsackievirus B3 3'UTR: La autoantigen binds the 3' and 5'UTR independently of the poly(A) tail

Coxsackievirus B3 (CVB3) is a positive, single-stranded RNA virus. The secondary structure of the 3' untranslated region (3'UTR) of CVB3 RNA consists of three stem-loops and is followed by a poly(A) tail sequence. These stem-loop structures have been suggested to participate in the regulation of viral replication through interaction with cellular proteins that are yet to be identified. In this study, by competitive UV cross-linking using mutated 3'UTR probes we have demonstrated that the poly(A) tail is essential for promoting HeLa cell protein interactions with the 3'UTR because deletion of this sequence abolished most of the protein interactions. Unexpectedly, mutations that disrupted the tertiary loop-loop interactions without affecting the stem-loops did not apparently affect these protein interactions, indicating that secondary structure rather than the high-order structure may play a major role in recruiting these RNA binding proteins. Among the observed 3'UTR RNA binding proteins, we have confirmed a 52 kDa protein as the human La autoantigen by using purified recombinant protein and a polyclonal La antibody. This protein can interact with both the 3' and 5'UTRs independently of the poly(A) tail. Further analysis by two-stage UV cross-linking, we found that the 3' and 5'UTR sequences may share the same binding site on the La protein.

Enterovirus 94, a proposed new serotype in human enterovirus species D

The genus Enterovirus (family Picornaviridae) contains five species with strains isolated from humans: Human enterovirus A (HEV-A), HEV-B, HEV-C, HEV-D and Poliovirus. In this study, a proposed new serotype of HEV-D was characterized. Four virus strains were isolated from sewage in Egypt and one strain from acute flaccid paralysis cases in the Democratic Republic of the Congo. The complete genome of one environmental isolate, the complete coding sequence of one clinical isolate and complete VP1 regions from the other isolates were sequenced. These isolates had 66.6–69.4 % nucleotide similarity and 74.7–76.6 % amino acid sequence similarity in the VP1 region with the closest enterovirus serotype, enterovirus 70 (EV70), suggesting that the isolates form a new enterovirus type, tentatively designated enterovirus 94 (EV94). Phylogenetic analyses including sequences of the 5′ UTR, VP1 and 3D regions demonstrated that EV94 isolates formed a monophyletic group within the species HEV-D. No evidence of recombination was found between EV94 and the other HEV-D serotypes, EV68 and EV70. Further biological characterization showed that EV94 was acid stable and had a wide cell tropism in vitro. Attempts to prevent replication with protective antibodies to known enterovirus receptors (poliovirus receptor, vitronectin α v β 3 receptor and decay accelerating factor) were not successful. Seroprevalence studies in the Finnish population revealed a high prevalence of this virus over the past two decades.

Analysis of the genetic and the corresponding antigenic variability of the VP1 3' end of ECHO virus type 11 and ECHO virus type 30

The enteroviruses (EV), RNA viruses belonging to the Picornaviridae family, have a high genetic variability due to the absence of the efficient proofreading and post replicative repair activities associated with the RNA polymerase. In the present work, we studied the genetic and the antigenic variability of ECHO virus types 11 (E11) and 30 (E30), which are the most isolated echoviruses serotypes in clinical and environmental samples. We established on the 3' end of the VP1 gene, consensus sequences of E11 and E30 by alignment of 67 E11 and 247 E30 published sequences in GenBank. Our results of sequences comparison showed that the majority of the mutational sites are situated on the third nucleotide of the codon. These mutations were without consequence on the antigenic sequences of the VP1 protein. Thus, E11 and E30 have a high genetic variability (1/3 of the nucleotides are variable), but a relative antigenic conservation. The analysis of the intertypic antigenic variability between E11 and E30 was obtained by the alignment of the corresponding amino acids sequences relative to the N-terminal part of the VP1 protein. Two discriminating parts were highlighted, probably representing antigenic sites for neutralisation antibodies.

Analysis of a new human parechovirus allows the definition of parechovirus types and the identification of RNA structural domains

Human parechoviruses (HPeV), members of the Parechovirus genus of Picornaviridae, are frequent pathogens but have been comparatively poorly studied, and little is known of their diversity, evolution, and molecular biology. To increase the amount of information available, we have analyzed 7 HPeV strains isolated in California between 1973 and 1992. We found that, on the basis of VP1 sequences, these fall into two genetic groups, one of which has not been previously observed, bringing the number of known groups to five. While these correlate partly with the three known serotypes, two members of the HPeV2 serotype belong to different genetic groups. In view of the growing importance of molecular techniques in diagnosis, we suggest that genotype is an important criterion for identifying viruses, and we propose that the genetic groups we have defined should be termed human parechovirus types 1 to 5. Complete nucleotide sequence analysis of two of the Californian isolates, representing two types, confirmed the identification of a new genetic group and suggested a role for recombination in parechovirus evolution. It also allowed the identification of a putative HPeV1 cis-acting replication element, which is located in the VP0 coding region, as well as the refinement of previously predicted 5' and 3' untranslated region structures. Thus, the results have significantly improved our understanding of these common pathogens.

Polyadenylation of genomic RNA and initiation of antigenomic RNA in a positive-strand RNA virus are controlled by the same cis-element

Genomes and antigenomes of many positive-strand RNA viruses contain 3'-poly(A) and 5'-poly(U) tracts, respectively, serving as mutual templates. Mechanism(s) controlling the length of these homopolymeric stretches are not well understood. Here, we show that in coxsackievirus B3 (CVB3) and three other enteroviruses the poly(A) tract is approximately 80-90 and the poly(U) tract is approximately 20 nt-long. Mutagenesis analysis indicate that the length of the CVB3 3'-poly(A) is determined by the oriR, a cis-element in the 3'-noncoding region of viral RNA. In contrast, while mutations of the oriR inhibit initiation of (-) RNA synthesis, they do not affect the 5'-poly(U) length. Poly(A)-lacking genomes are able to acquire genetically unstable AU-rich poly(A)-terminated 3'-tails, which may be generated by a mechanism distinct from the cognate viral RNA polyadenylation. The aberrant tails ensure only inefficient replication. The possibility of RNA replication independent of oriR and poly(A) demonstrate that highly debilitated viruses are able to survive by utilizing 'emergence', perhaps atavistic, mechanisms.

Molecular analysis of an echovirus 3 strain isolated from an individual concurrently with appearance of islet cell and IA-2 autoantibodies

Growing evidence has implicated members of the genus Enterovirus of the family Picornaviridae in the etiology of some cases of type 1 diabetes (T1D). To contribute to an understanding of the molecular determinants underlying this association, we determined the complete nucleotide sequence of a strain of echovirus 3 (E3), Human enterovirus B (HEV-B) species, isolated from an individual who soon after virus isolation developed autoantibodies characteristic of T1D. The individual has remained positive for over 6 years for tyrosine phosphatase-related IA-2 protein autoantibodies and islet cell autoantibodies, indicating an ongoing autoimmune process, although he has not yet developed clinical T1D. The sequence obtained adds weight to the observation that recent enterovirus isolates differ significantly from prototype strains and provides further evidence of a role for recombination in enterovirus evolution. In common with most HEV-B species members, the isolate exhibits 2C and VP1 sequences suggested as triggers of autoimmunity through molecular mimicry. However, comparisons with the E3 prototype strain and previously reported diabetogenic and nondiabetogenic HEV-B strains do not reveal clear candidates for sequence features of PicoBank/DM1/E3 that could be associated with autoantibody appearance. This is the first time a virus strain isolated at the time of commencement of beta-cell damage has been analyzed and is an invaluable addition to enterovirus strains isolated previously at the onset of T1D in the search for specific molecular features which could be associated with diabetes induction.

Structural and functional integrity of the coxsackievirus B3 oriR: spacing between coaxial RNA helices

The enterovirus oriR is composed of two helices, X and Y, anchored by a kissing (K) interaction. For proper oriR function, certain areas of these helices should be specifically oriented towards each other. It was hypothesized that the single-stranded nucleotides bridging the coaxial helices (Y-X and K-Y linkers) are important to determine this orientation. Spatial changes were introduced by altering the linker length between the helices of the coxsackievirus B3 oriR. Changing the linker lengths resulted in defective RNA replication, probably because of an altered oriR geometry. The identity of the linker residues also played a role, possibly because of sequence-specific ligand recognition. Although each point mutation altering the primary sequence of the Y-X spacer resulted in defective growth at 36 degrees C, the mutations had a wild-type phenotype at 39 degrees C, indicating a cold-sensitive phenotype. The results show that the intrinsic connection between oriR structure and function is fine-tuned by the spacing between the coaxial RNA helices.

Structural and functional characterization of the coxsackievirus B3 CRE(2C): role of CRE(2C) in negative- and positive-strand RNA synthesis

A stem-loop element located within the 2C-coding region of the coxsackievirus B3 (CVB3) genome has been proposed to function as a cis-acting replication element (CRE). It is shown here that disruption of this structure indeed interfered with viral RNA replication in vivo and abolished uridylylation of VPg in vitro. Site-directed mutagenesis demonstrated that the previously proposed enteroviral CRE consensus loop sequence, R(1)NNNAAR(2)NNNNNNR(3), is also applicable to CVB3 CRE(2C) and that a positive correlation exists between the ability of CRE(2C) mutants to serve as template in the uridylylation reaction and the capacity of these mutants to support viral RNA replication. To further investigate the effects of the mutations on negative-strand RNA synthesis, an in vitro translation/replication system containing HeLa S10 cell extracts was used. Similar to the results observed for poliovirus and rhinovirus, it was found that a complete disruption of the CRE(2C) structure interfered with positive-strand RNA synthesis, but not with negative-strand synthesis. All CRE(2C) point mutants affecting the enteroviral CRE consensus loop, however, showed a marked decrease in efficiency to induce negative-strand synthesis. Moreover, a transition (A(5)G) regarding the first templating adenosine residue in the loop was even unable to initiate complementary negative-strand synthesis above detectable levels. Taken together, these results indicate that the CVB3 CRE(2C) is not only required for the initiation of positive-strand RNA synthesis, but also plays an essential role in the efficient initiation of negative-strand RNA synthesis, a conclusion that has not been reached previously by using the cell-free system.

Species-specific RT-PCR amplification of human enteroviruses: a tool for rapid species identification of uncharacterized enteroviruses

The 65 serotypes of human enteroviruses are classified into four species, Human enterovirus (HEV) A to D, based largely on phylogenetic relationships in multiple genome regions. The 3'-non-translated region of enteroviruses is highly conserved within a species but highly divergent between species. From this information, species-specific RT-PCR primers were developed that can be used to rapidly screen collections of enterovirus isolates to identify species of interest. The four primer pairs were 100 % specific when tested against enterovirus prototype strains and panels of isolates of known serotype (a total of 193 isolates). For evaluation in a typical application, the species-specific primers were used to screen 186 previously uncharacterized non-polio enterovirus isolates. The HEV-B primers amplified 68.3 % of isolates, while the HEV-A and HEV-C primers accounted for 9.7 and 11.3 % of isolates, respectively; no isolates were amplified with the HEV-D primers. Twelve isolates (6.5 %) were amplified by more than one primer set and eight isolates (4.3 %) were not amplified by any of the four primer pairs. Serotypes were identified by partial sequencing of the VP1 capsid gene, and in every case sequencing confirmed that the species-specific PCR result was correct; the isolates that were amplified by more than one species-specific primer pair were mixtures of two (11 isolates) or three (one isolate) species of viruses. The eight isolates that were not amplified by the species-specific primers comprised four new serotypes (EV76, EV89, EV90 and EV91) that appear to be unique members of HEV-A based on VP1, 3D and 3'-non-translated region sequences.

Molecular characterization of a Canadian human parechovirus (HPeV)-3 isolate and its relationship to other HPeVs

A new member of the Parechovirus genus (HPeV-3) has been recently isolated from a stool specimen of a young Japanese child with transient paralysis (isolate A308/99). Three additional HPeV-3 Canadian cases associated with severe neonatal infections were subsequently described by our group. At this time, limited information is available on the molecular characterization of the HPeV-3 genome. In this study, we report the complete genome sequence of a Canadian HPeV-3 isolate (Can82853-01) recovered in 2001 from the nasopharyngeal aspirate of a neonate with a sepsis-like syndrome. The Can82853-01 genome was 7,322 nucleotides (nt) in length and had 96.4%, 77.7%, and 77.5% nt identity with the A308/99 (HPeV-3), HPeV-1, and HPeV-2 reference strains, respectively. Furthermore, the Can82853-01 polyprotein had 98.2%, 86.9%, and 84.7% amino acids identity compared to the same viruses. Most differences between the polyproteins of HPeV-3 isolates and those of HPeV-1 and -2 were located in the VP0-VP3-VP1 capsid region whereas the non-structural proteins were relatively conserved. The primary and predicted secondary structures of 5'UTR and 3'UTR of the HPeV-3 isolates were also different than those of HPeV-1 and -2 reference strains. Phylogenetic analysis confirmed the relatedness between Can82853-01 and A308/99 and their divergence from other picornaviruses.

Relationship between poliovirus negative-strand RNA synthesis and the length of the 3' poly(A) tail

The precise relationship between the length of the 3' poly(A) tail and the replication and infectivity of poliovirus RNA was examined in this study. With both poly(A)(11) and poly(A)(12) RNAs, negative-strand synthesis was 1-3% of the level observed with poly(A)(80) RNA. In contrast, increasing the length of the poly(A) tail from (A)(12) to (A)(13) resulted in about a ten-fold increase in negative-strand synthesis. This increase continued with each successive increase in poly(A) tail length. With poly(A)(20) RNA, RNA synthesis approached the level observed with poly(A)(80) RNA. A similar relationship was observed between poly(A) tail length and the infectivity of the viral RNA. A replication model is described which suggests that viral RNA replication is dependent on a poly(A) tail that is long enough to bind poly(A) binding protein and to act as a template for VPg uridylylation and negative-strand initiation.

An authentic 3' noncoding region is necessary for efficient poliovirus replication

Picornavirus RNA replication involves the specific synthesis of negative-strand intermediates followed by an accumulation of positive-strand viral RNA in the presence of a multitude of cellular mRNAs. Previously, in an effort to identify cis-acting elements required for initiation of negative-strand RNA synthesis, we deleted the entire 3' noncoding regions from human rhinovirus and poliovirus genomic RNAs. These deletion mutation transcripts displayed a severe delay in RNA accumulation following transfection of HeLa cells. Interestingly, in subsequent infection of HeLa cells, the deletion-mutant poliovirus displayed only a moderate deficiency in RNA synthesis. These data suggested that the delay in the production of cytopathic effects after transfection may have been due to an RNA replication defect overcome by the accumulation of a compensatory mutation(s) generated during initial rounds of RNA synthesis. In this study, we have sequenced the entire genome of the deletion-mutant virus and found only two nucleotide changes from the parental clone. Transfection analysis of these sequence variants revealed that the sequence changes did not provide compensatory functions for the 3' noncoding region deletion mutation replication defect. Further examination of the deletion mutant phenotype revealed that the severe replication defect following RNA transfection is due, in part, to nonviral terminal sequences present in the in vitro-derived deletion mutation transcripts. Our data suggest that poliovirus RNA harboring a complete 3' noncoding region deletion mutation is infectious (not merely quasi-infectious).

Human enterovirus 71 subgenotype B3 lacks coxsackievirus A16-like neurovirulence in mice infection

Background

At least three different EV-71 subgenotypes were identified from an outbreak in Malaysia in 1998. The subgenotypes C2 and B4 were associated with the severe and fatal infections, whereas the B3 virus was associated with mild to subclinical infections. The B3 virus genome sequences had ≥85% similarity at the 3' end to CV-A16. This offers opportunities to examine if there are characteristic similarities and differences in virulence between CV-A16, EV-71 B3 and EV-71 B4 and to determine if the presence of the CV-A16-liked genes in EV-71 B3 would also confer the virus with a CV-A16-liked neurovirulence in mice model infection.

Results

Analysis of human enterovirus 71 (EV-71) subgenotype B3 genome sequences revealed that the 3D RNA polymerase and domain Z of the 3'-untranslating region RNA secondary structure had high similarity to CV-A16. Intracerebral inoculation of one-day old mice with the virus resulted in 16% of the mice showing swollen hind limbs and significantly lower weight gain in comparison to EV-71 B4-infected mice. None of the mice presented with hind leg paralysis typical in all the CV-A16 infected mice. CV-A16 genome sequences were amplified from the CV-A16-infected mice brain but no amplification was obtained from all the EV-71-inoculated mice suggesting that no replication had taken place in the suckling mice brain.

Conclusion

The findings presented here suggest that EV-71 B3 viruses had CV-A16-liked non-structural gene features at the 3'-end of the genome. Their presence could have affected virulence by affecting the mice general health but was insufficient to confer the EV-71 B3 virus a CV-A16-liked neurovirulence in mice model infection.

Feline calicivirus VP2 is essential for the production of infectious virions

The third open reading frame (ORF3) located at the 3' end of the genomic RNA of feline calicivirus (FCV) encodes a small (12.2-kDa) minor structural protein of 106 amino acids designated VP2. Point mutations and deletions were introduced into an infectious FCV cDNA clone in order to evaluate the functional importance of ORF3 and its encoded protein, VP2. Deletion of the entire ORF3 sequence was lethal for the virus, and evidence was found for strong selective pressure to produce the VP2 protein. Extended deletions in the 5' end and small deletions in the 3' end of ORF3, as well as the introduction of stop codons into the ORF3 sequence, were tolerated by the viral replication machinery, but infectious virus could not be recovered. Infectious virus particles could be rescued from a full-length FCV cDNA clone encoding a nonfunctional VP2 when VP2 was provided in trans from a eukaryotic expression plasmid. Our data indicate that VP2, a protein apparently unique to the caliciviruses, is essential for productive replication that results in the synthesis and maturation of infectious virions and that the ORF3 nucleotide sequence itself overlaps a cis-acting RNA signal at the genomic 3' end.

3'-Terminal sequence in poliovirus negative-strand templates is the primary cis-acting element required for VPgpUpU-primed positive-strand initiation

The 5' cloverleaf in poliovirus RNA has a direct role in regulating the stability, translation, and replication of viral RNA. In this study, we investigated the role of stem a in the 5' cloverleaf in regulating the stability and replication of poliovirus RNA in HeLa S10 translation-replication reactions. Our results showed that disrupting the duplex structure of stem a destabilized viral RNA and inhibited efficient negative-strand synthesis. Surprisingly, the duplex structure of stem a was not required for positive-strand synthesis. In contrast, altering the primary sequence at the 5'-terminal end of stem a had little or no effect on negative-strand synthesis but dramatically reduced positive-strand initiation and the formation of infectious virus. The inhibition of positive-strand synthesis observed in these reactions was most likely a consequence of nucleotide alterations in the conserved sequence at the 3' ends of negative-strand RNA templates. Previous studies suggested that VPgpUpU synthesized on the cre(2C) hairpin was required for positive-strand synthesis. Therefore, these results are consistent with a model in which preformed VPgpUpU serves as the primer for positive-strand initiation on the 3'AAUUUUGUC5' sequence at the 3' ends of negative-strand templates. Our results suggest that this sequence is the primary cis-acting element that is required for efficient VPgpUpU-primed positive-strand initiation.

A phosphorothioate antisense oligodeoxynucleotide specifically inhibits coxsackievirus B3 replication in cardiomyocytes and mouse hearts

Antisense oligodeoxynucleotides (AS-ODNs) are promising therapeutic agents for the treatment of virus-induced diseases. We previously reported that coxsackievirus B3 (CVB3) infectivity could be inhibited effectively in HeLa cells by phosphorothioate AS-ODNs complementary to different regions of the 5' and 3' untranslated regions of CVB3 RNA. The most effective target is the proximal terminus of the 3' untranslated region. To further investigate the potential antiviral role of the AS-ODN targeting this site in cardiomyocytes (HL-1 cell line), corresponding AS-ODN (AS-7) was transfected into the HL-1 cells and followed by CVB3 infection. Analyses by RT-PCR, Western blotting and plaque assay demonstrated that AS-7 strongly inhibits viral RNA and viral protein synthesis as compared to scrambled AS-ODNs. The percent inhibitions of viral RNA transcription and capsid protein VP1 synthesis were 87.6 and 40.1, respectively. Moreover, AS-7 could inhibit ongoing CVB3 infection when it was given after virus infection. The antiviral activity was further evaluated in a CVB3 myocarditis mouse model. Adolescent A/J mice were intravenously administrated with AS-7 or scrambled AS-ODNs prior to and after CVB3 infection. Following a 4-day therapy, the myocardium CVB3 RNA replication decreased by 68% and the viral titers decreased by 0.5 log(10) in the AS-7-treated group as compared to the group treated with the scrambled AS-ODNs as determined by RT-PCR, in situ hybridization and viral plaque assay. Taken together, our results demonstrated a great potential for AS-7 to be further developed into an effective treatment towards viral myocarditis as well as other diseases caused by CVB3 infection.

Genetic evidence for an interaction between a picornaviral cis-acting RNA replication element and 3CD protein

Internally located, cis-acting RNA replication elements, termed cres, are essential for replication of the genomes of picornaviruses such as human rhinovirus 14 (HRV-14) and poliovirus because they template uridylylation of the protein primer, VPg, by the polymerase 3D(pol). These cres form stem-loop structures sharing a common loop motif, and the HRV-14 cre can substitute functionally for the poliovirus cre in both uridylylation in vitro and RNA replication in vivo. We show, however, that the poliovirus cre is unable to support HRV-14 RNA replication. This lack of complementation maps to the stem of the poliovirus cre and was reversed by single nucleotide substitutions in the stem as well as the base of the loop. Replication-competent, revertant viruses rescued from dicistronic HRV-14 RNAs containing the poliovirus cre, or a chimeric cre containing the poliovirus stem, contained adaptive amino acid substitutions. These mapped to the surface of both the polymerase 3D(pol), at the tip of the "thumb" domain, and the protease 3C(pro), on the side opposing the active site and near the end of an extended strand segment implicated previously in RNA binding. These mutations substantially enhanced replication competence when introduced into HRV-14 RNAs containing the poliovirus cre, and they were additive in their effects. The data support a model in which 3CD or its derivatives 3C(pro) and 3D(pol) interact directly with the stem of the cre during uridylylation of VPg.

Poliovirus cre(2C)-dependent synthesis of VPgpUpU is required for positive- but not negative-strand RNA synthesis

The cre(2C) hairpin is a cis-acting replication element in poliovirus RNA and serves as a template for the synthesis of VPgpUpU. We investigated the role of the cre(2C) hairpin on VPgpUpU synthesis and viral RNA replication in preinitiation RNA replication complexes isolated from HeLa S10 translation-RNA replication reactions. cre(2C) hairpin mutations that block VPgpUpU synthesis in reconstituted assays with purified VPg and poliovirus polymerase were also found to completely inhibit VPgpUpU synthesis in preinitiation replication complexes. Surprisingly, blocking VPgpUpU synthesis by mutating the cre(2C) hairpin had no significant effect on negative-strand synthesis but completely inhibited positive-strand synthesis. Negative-strand RNA synthesized in these reactions immunoprecipitated with anti-VPg antibody and demonstrated that it was covalently linked to VPg. This indicated that VPg was used to initiate negative-strand RNA synthesis, although the cre(2C)-dependent synthesis of VPgpUpU was inhibited. Based on these results, we concluded that the cre(2C)-dependent synthesis of VPgpUpU was required for positive- but not negative-strand RNA synthesis. These findings suggest a replication model in which negative-strand synthesis initiates with VPg uridylylated in the 3' poly(A) tail in virion RNA and positive-strand synthesis initiates with VPgpUpU synthesized on the cre(2C) hairpin. The pool of excess VPgpUpU synthesized on the cre(2C) hairpin should support high levels of positive-strand synthesis and thereby promote the asymmetric replication of poliovirus RNA.

Poliovirus CRE-dependent VPg uridylylation is required for positive-strand RNA synthesis but not for negative-strand RNA synthesis

The cis-acting replication element (CRE) is a 61-nucleotide stem-loop RNA structure found within the coding sequence of poliovirus protein 2C. Although the CRE is required for viral RNA replication, its precise role(s) in negative- and positive-strand RNA synthesis has not been defined. Adenosine in the loop of the CRE RNA structure functions as the template for the uridylylation of the viral protein VPg. VPgpUpU(OH), the predominant product of CRE-dependent VPg uridylylation, is a putative primer for the poliovirus RNA-dependent RNA polymerase. By examining the sequential synthesis of negative- and positive-strand RNAs within preinitiation RNA replication complexes, we found that mutations that disrupt the structure of the CRE prevent VPg uridylylation and positive-strand RNA synthesis. The CRE mutations that inhibited the synthesis of VPgpUpU(OH), however, did not inhibit negative-strand RNA synthesis. A Y3F mutation in VPg inhibited both VPgpUpU(OH) synthesis and negative-strand RNA synthesis, confirming the critical role of the tyrosine hydroxyl of VPg in VPg uridylylation and negative-strand RNA synthesis. trans-replication experiments demonstrated that the CRE and VPgpUpU(OH) were not required in cis or in trans for poliovirus negative-strand RNA synthesis. Because these results are inconsistent with existing models of poliovirus RNA replication, we propose a new four-step model that explains the roles of VPg, the CRE, and VPgpUpU(OH) in the asymmetric replication of poliovirus RNA.

Complete sequence analyses of enterovirus 71 strains from fatal and non-fatal cases of the hand, foot and mouth disease outbreak in Singapore (2000)

Enterovirus 71 (EV71) is a major aetiological agent of hand, foot and mouth disease (HFMD). In recent years, several outbreaks in East Asia were associated with neurological complications and numerous deaths. An outbreak in Singapore in October 2000 afflicted thousands of children, resulting in four fatal cases from three of whom EV71 was isolated. The genomes of two representative EV71 strains isolated from a fatal case and a surviving patient were completely sequenced, and their nucleotide and amino acid sequences compared with known EV71 strains. The two outbreak strains were classified under genogroup B, together with those previously isolated in Singapore, Malaysia and Japan. Comparative sequence analysis of the two Singapore strains revealed 99% nucleotide similarity, while their deduced amino acid sequences were almost identical except for residue 1506 in the 3A non-structural region. Given that the outbreak involved closely related genetic variants of EV71, the broad spectrum of disease severity may be attributed to critical factors such as varying viral inoculation doses or differing host immune responses following infection, but is less likely to be due to the emergence of EV71 strains with heightened virulence.

Interaction of glyceraldehyde-3-phosphate dehydrogenase with secondary and tertiary RNA structural elements of the hepatitis A virus 3' translated and non-translated regions

Proteins interacting with RNA structures at the 3' non-translated region (3'NTR) of picornaviruses are probably important during viral RNA replication. We have shown previously that a dominant cellular cytoplasmic protein of 38 kDa (p38) interacts with the 3'NTR and upstream regions of the hepatitis A virus (HAV) RNA (Kusov et al., J Virol 70, 1890-1897, 1996). Immunological and biochemical analyses of p38 have indicated that it is identical to GAPDH, which has previously been described as modulating translational regulation of the HAV RNA by interacting with the 5'NTR (Schultz et al., J Biol Chem 271, 14134-14142, 1996). Three separate binding regions for GAPDH in the 3'NTR and in the upstream 3D polymerase-coding region were identified. Structural analysis of these RNA regions by computer modelling and direct enzymatic cleavage suggested the presence of several AU-rich stem-loop structures having the potential for tertiary interactions. Binding of GAPDH to these structures was confirmed by RNA footprint analysis and resulted in the loss of double-stranded RNA regions. A different panel of RNA binding proteins (p28, p41 and p65) was detected in the ribosomal fractions of several cell lines (BSC-1, FRhK-4 and HeLa), whereas RNA binding of the GAPDH that was also present in these fractions was only marginal or absent.

Initiation of poliovirus negative-strand RNA synthesis requires precursor forms of p2 proteins

The replication proteins encoded in the P2 region of the poliovirus genome induce extensive rearrangement of cellular membranes into vesicles and are a required component of viral RNA replication complexes. To identify distinct viral protein(s) from the P2 region of the genome that were required to form functional RNA replication complexes, the P2 proteins were expressed in addition to P3 in HeLa S10 translation-RNA replication reactions. Membrane-associated preinitiation replication complexes were isolated from these reactions and used to measure negative-strand synthesis. The formation of replication complexes capable of initiating negative-strand synthesis was observed when either P23 or when P2 and P3 were expressed in the HeLa S10 translation-replication reactions. The amount of negative-strand RNA synthesized with P2 and P3 was approximately 50% of that observed with P23. Negative-strand synthesis was not observed when the processed forms of the P2 proteins (e.g., 2A, 2B, 2C, 2AB, and 2BC) were used in various combinations in place of P2. In contrast, the expression of 2A and 2BCP3 supported negative-strand synthesis at the same level observed with P23. Therefore, functional replication complexes were formed in reaction mixtures that contained either 2A and 2BCP3 or P2 and P3. Genetic complementation analysis of P23 RNA that contained a lethal mutation in 2C confirmed these results. The expression of 2BCP3 in trans restored the replication of P23-2C(P131N) RNA to wild-type levels. The expression of P2 and P3 also complemented the replication of this mutant RNA, although very inefficiently. Complementation was not observed in reactions that contained P2 alone, 2BC, or 2C. Based on these results, we propose that RNA replication complexes are initially formed with the primary cleavage products of P23 (i.e., P2 and P3 or 2A and 2BCP3), and that 2A and 2BCP3 are preferentially used in this process.

Terminal RNA replication elements in human parechovirus 1

To define structural elements critical for RNA replication in human parechovirus 1 (HPeV1), a replicon with chloramphenicol acetyltransferase as a reporter gene and an infectious virus cDNA clone have been used. It was observed that there are cis-acting signals required for HPeV1 replication located within the 5'-terminal 112 nucleotides of the genome and that these include two terminal stem-loops, SL-A and SL-B, together with a pseudoknot element. Significant disruption of any of these structures impaired both RNA replication and virus growth. In view of the similarity in terminal structures to several picornaviruses, such as cardioviruses and hepatoviruses, the insights generated in this work are of wider significance for understanding picornavirus replication.

Biological significance of a human enterovirus B-specific RNA element in the 3' nontranslated region

The secondary structures predicted for the enteroviral 3' nontranslated region (3'NTR) all seem to indicate a conformation consisting of two (X and Y) hairpin structures. The higher-order RNA structure of the 3'NTR appears to exist as an intramolecular kissing interaction between the loops of these two hairpin structures. The enterovirus B-like subgroup possesses an additional stem-loop structure, domain Z, which is not present in the poliovirus-like enteroviruses. It has been suggested that the Z domain originated from a burst of short sequence repetitions (E. V. Pilipenko, S. V. Maslova, A. N. Sinyakov, and V. I. Agol, Nucleic Acids Res. 20:1739-1745, 1992). However, no functional features have yet been ascribed to this enterovirus B-like-specific RNA element in the 3'NTR. In this study, we tested the functional characteristics and biological significance of domain Z. A mutant of the cardiovirulent coxsackievirus group B3 strain Nancy which completely lacked the Z domain and which therefore acquired enterovirus C-like secondary structures exhibited a wild-type growth phenotype, as determined by single-cycle growth analysis with BGM cells. This result proves that the Z domain is virtually dispensable for viral growth in tissue cultures. Partial distortion of the Z domain structure resulted in a disabled virus with reduced growth kinetics, probably due to alternative conformations of the overall structure of the domain. Infection of mice showed that the recombinant coxsackievirus group B3 mutant which completely lacked the Z domain was less virulent. Pancreatic tissues from mice infected with wild-type virus and recombinant virus were equally affected. However, the heart tissue from mice infected with the recombinant virus showed only slight signs of myocarditis. These results suggest that the enterovirus B-like-specific Z domain plays a role in coxsackievirus-induced pathogenesis.

Sequence requirements for viral RNA replication and VPg uridylylation directed by the internal cis-acting replication element (cre) of human rhinovirus type 14

Until recently, the cis-acting signals required for replication of picornaviral RNAs were believed to be restricted to the 5' and 3' noncoding regions of the genome. However, an RNA stem-loop in the VP1-coding sequence of human rhinovirus type 14 (HRV-14) is essential for viral minus-strand RNA synthesis (K. L. McKnight and S. M. Lemon, RNA 4:1569-1584, 1998). The nucleotide sequence of the apical loop of this internal cis-acting replication element (cre) was critical for RNA synthesis, while secondary RNA structure, but not primary sequence, was shown to be important within the duplex stem. Similar cres have since been identified in other picornaviral genomes. These RNA segments appear to serve as template for the uridylylation of the genome-linked protein, VPg, providing the VPg-pUpU primer required for viral RNA transcription (A. V. Paul et al., J. Virol. 74:10359-10370, 2000). Here, we show that the minimal functional HRV-14 cre resides within a 33-nucleotide (nt) RNA segment that is predicted to form a simple stem-loop with a 14-nt loop sequence. An extensive mutational analysis involving every possible base substitution at each position within the loop segment defined the sequence that is required within this loop for efficient replication of subgenomic HRV-14 replicon RNAs. These results indicate that three consecutive adenosine residues (nt 2367 to 2369) within the 5' half of this loop are critically important for cre function and suggest that a common RNNNAARNNNNNNR loop motif exists among the cre sequences of enteroviruses and rhinoviruses. We found a direct, positive correlation between the capacity of mutated cres to support RNA replication and their ability to function as template in an in vitro VPg uridylylation reaction, suggesting that these functions are intimately linked. These data thus define more precisely the sequence and structural requirements of the HRV-14 cre and provide additional support for a model in which the role of the cre in RNA replication is to act as template for VPg uridylylation.

Sequencing of porcine enterovirus groups II and III reveals unique features of both virus groups

The molecular classification of the porcine enterovirus (PEV) groups II and III was investigated. The sequence of the almost complete PEV-8 (group II) genome reveals that this virus has unique L and 2A gene regions. A reclassification of this group into a new picornavirus genus is suggested. PEV group III viruses are typical enteroviruses. They differ from other enteroviruses by a prolonged stem-loop D of the 5'-cloverleaf structure.

The 3prime prime or minute-terminal structure required for replication of Barley yellow dwarf virus RNA contains an embedded 3prime prime or minute end

We determined the 3prime prime or minute-terminal primary and secondary structures required for replication of Barley yellow dwarf virus (BYDV) RNA in oat protoplasts. Computer predictions, nuclease probing, phylogenetic comparisons, and replication assays of specific mutants and chimeras revealed that the 3prime prime or minute-terminal 109 nucleotides (nt) form a structure with three to four stem-loops followed by a coaxially stacked helix incorporating the last four nt [(A/U)CCC]. Sequences upstream of the 109-nt region also contributed to RNA accumulation. The base-pairing but not the sequences or bulges in the stems were essential for replication, but any changes to the 3prime prime or minute-terminal helix destroyed replication. The two 3prime prime or minute-proximal tetraloops tolerated all changes, but the two 3prime prime or minute-distal tetraloops gave most efficient replication if they fit the GNRA consensus. A mutant lacking the 3prime prime or minute-proximal stem-loop produced elevated levels of less-than-full-length minus strands, and no (+) strand. We propose that a "pocket" structure is the origin of (minus sign)-strand synthesis, which is negatively regulated by the inaccessible conformation of the 3prime prime or minute terminus, thus favoring a high (+)/(minus sign) ratio. This 3prime prime or minute structure and the polymerase homologies suggest that genus Luteovirus is more closely related to the Tombusviridae family than to other Luteoviridae genera.

Conserved RNA secondary structures in Picornaviridae genomes

The family Picornaviridae contains important pathogens including, for example, hepatitis A virus and foot-and-mouth disease virus. The genome of these viruses is a single messenger-active (+)-RNA of 7200-8500 nt. Besides coding for the viral proteins, it also contains functionally important RNA secondary structures, among them an internal ribosomal entry site (IRES) region towards the 5'-end. This contribution provides a comprehensive computational survey of the complete genomic RNAs and a detailed comparative analysis of the conserved structural elements in seven of the currently nine genera in the family PICORNAVIRIDAE: Compared with previous studies we find: (i) that only smaller sections of the IRES region than previously reported are conserved at single base-pair resolution and (ii) that there is a number of significant structural elements in the coding region. Furthermore, we identify potential cis-acting replication elements in four genera where this feature has not been reported so far.