Terminal RNA replication elements in human parechovirus 1

Evolutionary and network analysis of virus sequences from infants infected with an Australian recombinant strain of human parechovirus type 3

We present the near complete virus genome sequences with phylogenetic and network analyses of potential transmission networks of a total of 18 Australian cases of human parechovirus type 3 (HPeV3) infection in infants in the period from 2012–2015. Overall the results support our previous finding that the Australian outbreak strain/lineage is a result of a major recombination event that took place between March 2012 and November 2013 followed by further virus evolution and possibly recombination. While the nonstructural coding region of unknown provenance appears to evolve significantly both at the nucleotide and amino acid level, the capsid encoding region derived from the Yamagata 2011 lineage of HPeV3 appears to be very stable, particularly at the amino acid level. The phylogenetic and network analyses performed support a temporal evolution from the first Australian recombinant virus sequence from November 2013 to March/April 2014, onto the 2015 outbreak. The 2015 outbreak samples fall into two separate clusters with a possible common ancestor between March/April 2014 and September 2015, with each cluster further evolving in the period from September to November/December 2015.

An outbreak of severe infections among Australian infants caused by a novel recombinant strain of human parechovirus type 3

Human parechovirus types 1–16 (HPeV1–16) are positive strand RNA viruses in the family Picornaviridae. We investigated a 2015 outbreak of HPeV3 causing illness in infants in Victoria, Australia. Virus genome was extracted from clinical material and isolates and sequenced using a combination of next generation and Sanger sequencing. The HPeV3 outbreak genome was 98.7% similar to the HPeV3 Yamagata 2011 lineage for the region encoding the structural proteins up to nucleotide position 3115, but downstream of that the genome varied from known HPeV sequences with a similarity of 85% or less. Analysis indicated that recombination had occurred, may have involved multiple types of HPeV and that the recombination event/s occurred between March 2012 and November 2013. However the origin of the genome downstream of the recombination site is unknown. Overall, the capsid of this virus is highly conserved, but recombination provided a different non-structural protein coding region that may convey an evolutionary advantage. The indication that the capsid encoding region is highly conserved at the amino acid level may be helpful in directing energy towards the development of a preventive vaccine for expecting mothers or antibody treatment of young infants with severe disease.

Picornaviruses and nuclear functions: targeting a cellular compartment distinct from the replication site of a positive-strand RNA virus

The compartmentalization of DNA replication and gene transcription in the nucleus and protein production in the cytoplasm is a defining feature of eukaryotic cells. The nucleus functions to maintain the integrity of the nuclear genome of the cell and to control gene expression based on intracellular and environmental signals received through the cytoplasm. The spatial separation of the major processes that lead to the expression of protein-coding genes establishes the necessity of a transport network to allow biomolecules to translocate between these two regions of the cell. The nucleocytoplasmic transport network is therefore essential for regulating normal cellular functioning. The Picornaviridae virus family is one of many viral families that disrupt the nucleocytoplasmic trafficking of cells to promote viral replication. Picornaviruses contain positive-sense, single-stranded RNA genomes and replicate in the cytoplasm of infected cells. As a result of the limited coding capacity of these viruses, cellular proteins are required by these intracellular parasites for both translation and genomic RNA replication. Being of messenger RNA polarity, a picornavirus genome can immediately be translated upon entering the cell cytoplasm. However, the replication of viral RNA requires the activity of RNA-binding proteins, many of which function in host gene expression, and are consequently localized to the nucleus. As a result, picornaviruses disrupt nucleocytoplasmic trafficking to exploit protein functions normally localized to a different cellular compartment from which they translate their genome to facilitate efficient replication. Furthermore, picornavirus proteins are also known to enter the nucleus of infected cells to limit host-cell transcription and down-regulate innate antiviral responses. The interactions of picornavirus proteins and host-cell nuclei are extensive, required for a productive infection, and are the focus of this review.

Predicting pseudoknotted structures across two RNA sequences

Motivation: Laboratory RNA structure determination is demanding and costly and thus, computational structure prediction is an important task. Single sequence methods for RNA secondary structure prediction are limited by the accuracy of the underlying folding model, if a structure is supported by a family of evolutionarily related sequences, one can be more confident that the prediction is accurate. RNA pseudoknots are functional elements, which have highly conserved structures. However, few comparative structure prediction methods can handle pseudoknots due to the computational complexity.

Results: A comparative pseudoknot prediction method called DotKnot-PW is introduced based on structural comparison of secondary structure elements and H-type pseudoknot candidates. DotKnot-PW outperforms other methods from the literature on a hand-curated test set of RNA structures with experimental support.

Availability: DotKnot-PW and the RNA structure test set are available at the web site http://dotknot.csse.uwa.edu.au/pw.

Contact:janaspe@csse.uwa.edu.au

Supplementary information:Supplementary data are available at Bioinformatics online.

Unusual loop-sequence flexibility of the proximal RNA replication element in EMCV

Picornaviruses contain stable RNA structures at the 5′ and 3′ ends of the RNA genome, OriL and OriR involved in viral RNA replication. The OriL RNA element found at the 5′ end of the enterovirus genome folds into a cloverleaf-like configuration. In vivo SELEX experiments revealed that functioning of the poliovirus cloverleaf depends on a specific structure in this RNA element. Little is known about the OriL of cardioviruses. Here, we investigated structural aspects and requirements of the apical loop of proximal stem-loop SL-A of mengovirus, a strain of EMCV. Using NMR spectroscopy, we showed that the mengovirus SL-A apical loop consists of an octaloop. In vivo SELEX experiments demonstrated that a large number of random sequences are tolerated in the apical octaloop that support virus replication. Mutants in which the SL-A loop size and the length of the upper part of the stem were varied showed that both stem-length and stability of the octaloop are important determinants for viral RNA replication and virus reproduction. Together, these data show that stem-loop A plays an important role in virus replication. The high degree of sequence flexibility and the lack of selective pressure on the octaloop argue against a role in sequence specific RNA-protein or RNA-RNA interactions in which octaloop nucleotides are involved.

DotKnot: pseudoknot prediction using the probability dot plot under a refined energy model

RNA pseudoknots are functional structure elements with key roles in viral and cellular processes. Prediction of a pseudoknotted minimum free energy structure is an NP-complete problem. Practical algorithms for RNA structure prediction including restricted classes of pseudoknots suffer from high runtime and poor accuracy for longer sequences. A heuristic approach is to search for promising pseudoknot candidates in a sequence and verify those. Afterwards, the detected pseudoknots can be further analysed using bioinformatics or laboratory techniques. We present a novel pseudoknot detection method called DotKnot that extracts stem regions from the secondary structure probability dot plot and assembles pseudoknot candidates in a constructive fashion. We evaluate pseudoknot free energies using novel parameters, which have recently become available. We show that the conventional probability dot plot makes a wide class of pseudoknots including those with bulged stems manageable in an explicit fashion. The energy parameters now become the limiting factor in pseudoknot prediction. DotKnot is an efficient method for long sequences, which finds pseudoknots with higher accuracy compared to other known prediction algorithms. DotKnot is accessible as a web server at http://dotknot.csse.uwa.edu.au.

Cis-acting RNA elements in human and animal plus-strand RNA viruses

The RNA genomes of plus-strand RNA viruses have the ability to form secondary and higher-order structures that contribute to their stability and to their participation in inter- and intramolecular interactions. Those structures that are functionally important are called cis-acting RNA elements because their functions cannot be complemented in trans. They can be involved not only in RNA/RNA interactions but also in binding of viral and cellular proteins during the complex processes of translation, RNA replication and encapsidation. Most viral cis-acting RNA elements are located in the highly structured 5'- and 3'-nontranslated regions of the genomes but sometimes they also extend into the adjacent coding sequences. In addition, some cis-acting RNA elements are embedded within the coding sequences far away from the genomic ends. Although the functional importance of many of these structures has been confirmed by genetic and biochemical analyses, their precise roles are not yet fully understood. In this review we have summarized what is known about cis-acting RNA elements in nine families of human and animal plus-strand RNA viruses with an emphasis on the most thoroughly characterized virus families, the Picornaviridae and Flaviviridae.

Interaction between polypeptide 3ABC and the 5'-terminal structural elements of the genome of Aichi virus: implication for negative-strand RNA synthesis

Secondary structural elements at the 5' end of picornavirus genomic RNA function as cis-acting replication elements and are known to interact specifically with viral P3 proteins in several picornaviruses. In poliovirus, ribonucleoprotein complex formation at the 5' end of the genome is required for negative-strand synthesis. We have previously shown that the 5'-end 115 nucleotides of the Aichi virus genome, which are predicted to fold into two stem-loops (SL-A and SL-C) and one pseudoknot (PK-B), act as a cis-acting replication element and that correct folding of these structures is required for negative-strand synthesis. In this study, we investigated the interaction between the 5'-terminal 120 nucleotides of the genome and the P3 proteins, 3AB, 3ABC, 3C, and 3CD, by gel shift assay and Northwestern analysis. The results showed that 3ABC and 3CD bound to the 5'-terminal region specifically. The binding of 3ABC was observed on both assays, while that of 3CD was detected only on Northwestern analysis. No binding of 3AB or 3C was observed. Binding assays using mutant RNAs demonstrated that disruption of the base pairings of the stem of SL-A and one of the two stem segments of PK-B (stem-B1) abolished the 3ABC binding. In addition, the specific nucleotide sequence of stem-B1 was responsible for the efficient 3ABC binding. These results suggest that the interaction of 3ABC with the 5'-terminal region of the genome is involved in negative-strand synthesis. On the other hand, the ability of 3CD to interact with the 5'-terminal region did not correlate with the RNA replication ability.

KnotSeeker: heuristic pseudoknot detection in long RNA sequences

Pseudoknots are folded structures in RNA molecules that perform essential functions as part of cellular transcription machinery and regulatory processes. The prediction of these structures in RNA molecules has important implications in antiviral drug design. It has been shown that the prediction of pseudoknots is an NP-complete problem. Practical structure prediction algorithms based on free energy minimization employ a restricted problem class and dynamic programming. However, these algorithms are computationally very expensive, and their accuracy deteriorates if the input sequence containing the pseudoknot is too long. Heuristic methods can be more efficient, but do not guarantee an optimal solution in regards to the minimum free energy model. We present KnotSeeker, a new heuristic algorithm for the detection of pseudoknots in RNA sequences as a preliminary step for structure prediction. Our method uses a hybrid sequence matching and free energy minimization approach to perform a screening of the primary sequence. We select short sequence fragments as possible candidates that may contain pseudoknots and verify them by using an existing dynamic programming algorithm and a minimum weight independent set calculation. KnotSeeker is significantly more accurate in detecting pseudoknots compared to other common methods as reported in the literature. It is very efficient and therefore a practical tool, especially for long sequences. The algorithm has been implemented in Python and it also uses C/C++ code from several other known techniques. The code is available from http://www.csse.uwa.edu.au/~datta/pseudoknot.

Replication of Ljungan virus in cell culture: The genomic 5′-end, infectious cDNA clones and host cell response to viral infections

Ljungan virus (LV) is a picornavirus recently isolated from bank voles (Clethrionomys glareolus). The previously uncharacterised 5'-end sequence of the LV genome was determined. Infectious cDNA clones were constructed of the wild type LV prototype strain 87-012 and of the cytolytically replicating cell culture adapted variant 87-012G. Virus generated from cDNA clones showed identical growth characteristics as uncloned virus stocks. Cell culture adapted LV, 87-012G, showed a clear cytopathic effect (CPE) at 3-4 days post-infection (p.i.). Virus titers, determined by plaque titration, increased however only within the first 18h p.i. Replication of LV (+) strand RNA was determined by real-time PCR and corresponded in time with increasing titers. In contrast, the amounts of the replication intermediate, the (-) strand, continued to increase until the cells showed CPE. This indicates separate controlling mechanisms for replication of LV (+) and (-) genome strands. Replication was also monitored by immunofluorescence (IF) staining. IF staining of both prototype 87-012 and the CPE causing 87-012G showed groups of 5-25 infected cells at 48h p.i., suggesting a, for picornaviruses, not previously described direct cell-to-cell transmission.

A highly divergent picornavirus in a marine mammal

Nucleic acids from an unidentified virus from ringed seals (Phoca hispida) were amplified using sequence-independent PCR, subcloned, and then sequenced. The full genome of a novel RNA virus was derived, identifying the first sequence-confirmed picornavirus in a marine mammal. The phylogenetic position of the tentatively named seal picornavirus 1 (SePV-1) as an outlier to the grouping of parechoviruses was found consistently in alignable regions of the genome. A mean protein sequence identity of only 19.3 to 30.0% was found between the 3D polymerase gene sequence of SePV-1 and those of other picornaviruses. The predicted secondary structure of the short 506-base 5'-untranslated region showed some attributes of a type IVB internal ribosome entry site, and the polyprotein lacked an apparent L peptide, both properties associated with the Parechovirus genus. The presence of two SePV-1 2A genes and of the canonical sequence required for cotranslational cleavage resembled the genetic organization of Ljungan virus. Minor genetic variants were detected in culture supernatants derived from 8 of 108 (7.4%) seals collected in 2000 to 2002, indicating a high prevalence of SePV-1 in this hunted seal population. The high level of genetic divergence of SePV-1 compared to other picornaviruses and its mix of characteristics relative to its closest relatives support the provisional classification of SePV-1 as the prototype for a new genus in the family Picornaviridae.

[Analysis of Aichi virus replication]

Aichi virus is a member of the Family Picornaviridae. This virus was first isolated in 1989 from a stool specimen from a patient with oyster-associated gastroenteritis in Aichi, Japan. We analyzed the function of the 5' terminal region of the genome and the leader protein in virus replication. The results indicate that both the 5' terminal region of the genome and the leader protein are involved in viral RNA replication and encapsidation.

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.

A heuristic approach for detecting RNA H-type pseudoknots

MOTIVATION

RNA H-type pseudoknots are ubiquitous pseudoknots that are found in almost all classes of RNA and thought to play very important roles in a variety of biological processes. Detection of these RNA H-type pseudoknots can improve our understanding of RNA structures and their associated functions. However, the currently existing programs for detecting such RNA H-type pseudoknots are still time consuming and sometimes even ineffective. Therefore, efficient and effective tools for detecting the RNA H-type pseudoknots are needed.

RESULTS

In this paper, we have adopted a heuristic approach to develop a novel tool, called HPknotter, for efficiently and accurately detecting H-type pseudoknots in an RNA sequence. In addition, we have demonstrated the applicability and effectiveness of HPknotter by testing on some sequences with known H-type pseudoknots. Our approach can be easily extended and applied to other classes of more general pseudoknots.

AVAILABILITY

The web server of our HPknotter is available for online analysis at http://bioalgorithm.life.nctu.edu.tw/HPKNOTTER/ CONTACT: cllu@mail.nctu.edu.tw, chiu@cc.nctu.edu.tw

The 5'-terminal region of the Aichi virus genome encodes cis-acting replication elements required for positive- and negative-strand RNA synthesis

Aichi virus is a member of the family Picornaviridae. It has already been shown that three stem-loop structures (SL-A, SL-B, and SL-C, from the 5' end) formed at the 5' end of the genome are critical elements for viral RNA replication. In this study, we further characterized the 5'-terminal cis-acting replication elements. We found that an additional structural element, a pseudoknot structure, is formed through base-pairing interaction between the loop segment of SL-B (nucleotides [nt] 57 to 60) and a sequence downstream of SL-C (nt 112 to 115) and showed that the formation of this pseudoknot is critical for viral RNA replication. Mapping of the 5'-terminal sequence of the Aichi virus genome required for RNA replication using a series of Aichi virus-encephalomyocarditis virus chimera replicons indicated that the 5'-end 115 nucleotides including the pseudoknot structure are the minimum requirement for RNA replication. Using the cell-free translation-replication system, we examined the abilities of viral RNAs with a lethal mutation in the 5'-terminal structural elements to synthesize negative- and positive-strand RNAs. The results showed that the formation of three stem-loops and the pseudoknot structure at the 5' end of the genome is required for negative-strand RNA synthesis. In addition, specific nucleotide sequences in the stem of SL-A or its complementary sequences at the 3' end of the negative-strand were shown to be critical for the initiation of positive-strand RNA synthesis but not for that of negative-strand synthesis. Thus, the 5' end of the Aichi virus genome encodes elements important for not only negative-strand synthesis but also positive-strand synthesis.