Full-Length Genomic Analysis of Korean Porcine Sapelovirus Strains

Porcine sapelovirus (PSV), a species of the genus Sapelovirus within the family Picornaviridae, is associated with diarrhea, pneumonia, severe neurological disorders, and reproductive failure in pigs. However, the structural features of the complete PSV genome remain largely unknown. To analyze the structural features of PSV genomes, the full-length nucleotide sequences of three Korean PSV strains were determined and analyzed using bioinformatic techniques in comparison with other known PSV strains. The Korean PSV genomes ranged from 7,542 to 7,566 nucleotides excluding the 3′ poly(A) tail, and showed the typical picornavirus genome organization; 5′untranslated region (UTR)-L-VP4-VP2-VP3-VP1-2A-2B-2C-3A-3B-3C-3D-3′UTR. Three distinct cis-active RNA elements, the internal ribosome entry site (IRES) in the 5′UTR, a cis-replication element (CRE) in the 2C coding region and 3′UTR were identified and their structures were predicted. Interestingly, the structural features of the CRE and 3′UTR were different between PSV strains. The availability of these first complete genome sequences for PSV strains will facilitate future investigations of the molecular pathogenesis and evolutionary characteristics of PSV.     

Effects of Naturally Occurring Six- and Twelve-Nucleotide Inserts on Newcastle Disease Virus Replication and Pathogenesis

Newcastle disease virus (NDV) isolates contain genomes of 15,186, 15,192 or 15,198 nucleotides (nt). The length differences reflect a 6-nt insert in the 5′ (downstream) non-translated region (NTR) of the N gene (15,192-nt genome) or a 12-nt insert in the ORF encoding the P and V proteins (causing a 4-amino acid insert; 15,198-nt genome). We evaluated the role of these inserts in the N and P genes on viral replication and pathogenicity by inserting them into genomes of two NDV strains that have natural genome lengths of 15,186 nt and represent two different pathotypes, namely the mesogenic strain Beaudette C (BC) and the velogenic strain GB Texas (GBT). Our results showed that the 6-nt and 12-nt inserts did not detectably affect N gene expression or P protein function. The inserts had no effect on the replication or virulence of the highly virulent GBT strain but showed modest degree of attenuation in mesogenic strain BC. We also deleted a naturally-occurring 6-nt insertion in the N gene from a highly virulent 15,192-nt genome-length virus, strain Banjarmasin. This resulted in reduced replication in vitro and reduced virulence in vivo. Thus, although these inserts had no evident effect on gene expression, protein function, or replication in vivo, they did affect virulence in two of the three tested strains.     

Tombusvirus recruitment of host translational machinery via the 3′ UTR

RNA viruses recruit the host translational machinery by different mechanisms that depend partly on the structure of their genomes. In this regard, the plus-strand RNA genomes of several different pathogenic plant viruses do not contain traditional translation-stimulating elements, i.e., a 5′-cap structure and a 3′-poly(A) tail, and instead rely on a 3′-cap-independent translational enhancer (3′CITE) located in their 3′ untranslated regions (UTRs) for efficient synthesis of viral proteins. We investigated the structure and function of the I-shaped class of 3′CITE in tombusviruses—also present in aureusviruses and carmoviruses—using biochemical and molecular approaches and we determined that it adopts a complex higher-order RNA structure that facilitates translation by binding simultaneously to both eukaryotic initiation factor (eIF) 4F and the 5′ UTR of the viral genome. The specificity of 3′CITE binding to eIF4F is mediated, at least in part, through a direct interaction with its eIF4E subunit, whereas its association with the viral 5′ UTR relies on complementary RNA–RNA base-pairing. We show for the first time that this tripartite 5′ UTR/3′CITE/eIF4F complex forms in vitro in a translationally relevant environment and is required for recruitment of ribosomes to the 5′ end of the viral RNA genome by a mechanism that shares some fundamental features with cap-dependent translation. Notably, our results demonstrate that the 3′CITE facilitates the initiation step of translation and validate a molecular model that has been proposed to explain how several different classes of 3′CITE function. Moreover, the virus–host interplay defined in this study provides insights into natural host resistance mechanisms that have been linked to 3′CITE activity.     

四株不同宿主来源的Class I基因3型新城疫病毒全基因组序列测定及比较分析

为了解析基因型相同但宿主来源不同的新城疫病毒的全基因组差异。本文采用RT-PCR方法分别获得4株（JS/3/09/Ch,ZJ/3/10/Ch,AH/2/10/Du,JS/9/08/Go）Class I 基因3型病毒的全基因组核苷酸序列,并与GenBank中已公布的Class I基因3型病毒全基因组序列进行比对分析。本实验4株病毒的基因组长度均为15198bp,在基因组1607～1608位有6碱基的缺失,在2381～2382位有12碱基的插入,裂解位点为112EQ/RQE/GRL117是标准弱毒特征。5株Class I基因3型病毒之间全基因组同源性超过93%;而与Class II弱毒株同源性最低只有72.2%;比较6个结构蛋白基因的同源性,NP基因的同源性最高(98.3%～96.4%),而P基因最低(96.1%～91.9%)。结果表明不同宿主来源的Class I基因3型新城疫病毒在遗传信息方面差异不大,但NP/F/L基因的变异幅度较P/M/HN基因明显。     

Discovery of functional genomic motifs in viruses with ViReMa–a Virus Recombination Mapper–for analysis of next-generation sequencing data

We developed an algorithm named ViReMa (Viral-Recombination-Mapper) to provide a versatile platform for rapid, sensitive and nucleotide-resolution detection of recombination junctions in viral genomes using next-generation sequencing data. Rather than mapping read segments of pre-defined lengths and positions, ViReMa dynamically generates moving read segments. ViReMa initially attempts to align the 5′ end of a read to the reference genome(s) with the Bowtie seed-based alignment. A new read segment is then made by either extracting any unaligned nucleotides at the 3′ end of the read or by trimming the first nucleotide from the read. This continues iteratively until all portions of the read are either mapped or trimmed. With multiple reference genomes, it is possible to detect virus-to-host or inter-virus recombination. ViReMa is also capable of detecting insertion and substitution events and multiple recombination junctions within a single read. By mapping the distribution of recombination events in the genome of flock house virus, we demonstrate that this information can be used to discover de novo functional motifs located in conserved regions of the viral genome.     

Identification of cis-Acting Nucleotides and a Structural Feature in West Nile Virus 3′-Terminus RNA That Facilitate Viral Minus Strand RNA Synthesis

The 3′-terminal nucleotides (nt) of West Nile virus (WNV) genomic RNA form a penultimate 16-nt small stem-loop (SSL) and an 80-nt terminal stem-loop (SL). These RNA structures are conserved in divergent flavivirus genomes. A previous in vitro study using truncated WNV 3′ RNA structures predicted a putative tertiary interaction between the 5′ side of the 3′-terminal SL and the loop of the SSL. Although substitution or deletion of the 3′ G (nt 87) within the SSL loop, which forms the only G-C pair in the predicted tertiary interaction, in a WNV infectious clone was lethal, a finding consistent with the involvement in a functionally relevant pseudoknot interaction, extensive mutagenesis of nucleotides in the terminal SL did not identify a cis-acting pairing partner for this SSL 3′ G. However, both the sequence and the structural context of two adjacent base pairs flanked by symmetrical internal loops in the 3′-terminal SL were shown to be required for efficient viral RNA replication. Nuclear magnetic resonance analysis confirmed the predicted SSL and SL structures but not the tertiary interaction. The SSL was previously reported to contain one of three eEF1A binding sites, and G87 in the SSL loop was shown to be involved in eEF1A binding. The nucleotides at the bottom part of the 3′-terminal SL switch between 3′ RNA-RNA and 3′-5′ RNA-RNA interactions. The data suggest that interaction of the 3′ SL RNA with eEF1A at three sites and a unique metastable structural feature may participate in regulating structural changes in the 3′-terminal SL.     

Complete Genome of a Porcine Calicivirus Strain in Anhui Province,China, Is Significantly Shorter than That of the Other Chinese Strain

Caliciviruses that cause diarrhea have been reported in both industrial and developing countries, including China, in recent years. Here, we report the complete genome of a porcine calicivirus strain, Ah-1, which is prevalent in swine groups in Anhui Province. This viral genome is 7,342 nucleotides (nt) long, excluding the poly(A) of the 3′ end, which is 202 nt shorter in the 3′ untranslated region (UTR) than that of the other Chinese porcine calicivirus strain (Ch-sw-sav1; GenBank accession number {"type":"entrez-nucleotide","attrs":{"text":"FJ387164","term_id":"227976393","term_text":"FJ387164"}}FJ387164), previously isolated in the Shanghai area, China, though they shared 98.8% sequence identity over the whole genome excluding the 202-nt-shorter region.     

Transmissible Gastroenteritis Coronavirus Genome Packaging Signal Is Located at the 5′ End of the Genome and Promotes Viral RNA Incorporation into Virions in a Replication-Independent Process

Preferential RNA packaging in coronaviruses involves the recognition of viral genomic RNA, a crucial process for viral particle morphogenesis mediated by RNA-specific sequences, known as packaging signals. An essential packaging signal component of transmissible gastroenteritis coronavirus (TGEV) has been further delimited to the first 598 nucleotides (nt) from the 5′ end of its RNA genome, by using recombinant viruses transcribing subgenomic mRNA that included potential packaging signals. The integrity of the entire sequence domain was necessary because deletion of any of the five structural motifs defined within this region abrogated specific packaging of this viral RNA. One of these RNA motifs was the stem-loop SL5, a highly conserved motif in coronaviruses located at nucleotide positions 106 to 136. Partial deletion or point mutations within this motif also abrogated packaging. Using TGEV-derived defective minigenomes replicated in trans by a helper virus, we have shown that TGEV RNA packaging is a replication-independent process. Furthermore, the last 494 nt of the genomic 3′ end were not essential for packaging, although this region increased packaging efficiency. TGEV RNA sequences identified as necessary for viral genome packaging were not sufficient to direct packaging of a heterologous sequence derived from the green fluorescent protein gene. These results indicated that TGEV genome packaging is a complex process involving many factors in addition to the identified RNA packaging signal. The identification of well-defined RNA motifs within the TGEV RNA genome that are essential for packaging will be useful for designing packaging-deficient biosafe coronavirus-derived vectors and providing new targets for antiviral therapies.     

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

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

Sequences within Both the 5′ UTR and Gag Are Required for Optimal In Vivo Packaging and Propagation of Mouse Mammary Tumor Virus (MMTV) Genomic RNA

Background

This study mapped regions of genomic RNA (gRNA) important for packaging and propagation of mouse mammary tumor virus (MMTV). MMTV is a type B betaretrovirus which preassembles intracellularly, a phenomenon distinct from retroviruses that assemble the progeny virion at cell surface just before budding such as the type C human and feline immunodeficiency viruses (HIV and FIV). Studies of FIV and Mason-Pfizer monkey virus (MPMV), a type D betaretrovirus with similar intracellular virion assembly processes as MMTV, have shown that the 5′ untranslated region (5′ UTR) and 5′ end of gag constitute important packaging determinants for gRNA.

Methodology

Three series of MMTV transfer vectors containing incremental amounts of gag or 5′ UTR sequences, or incremental amounts of 5′ UTR in the presence of 400 nucleotides (nt) of gag were constructed to delineate the extent of 5′ sequences that may be involved in MMTV gRNA packaging. Real time PCR measured the packaging efficiency of these vector RNAs into MMTV particles generated by co-transfection of MMTV Gag/Pol, vesicular stomatitis virus envelope glycoprotein (VSV-G Env), and individual transfer vectors into human 293T cells. Transfer vector RNA propagation was monitored by measuring transduction of target HeLaT4 cells following infection with viral particles containing a hygromycin resistance gene expression cassette on the packaged RNA.

Principal Findings

MMTV requires the entire 5′ UTR and a minimum of ∼120 nucleotide (nt) at the 5′ end of gag for not only efficient gRNA packaging but also propagation of MMTV-based transfer vector RNAs. Vector RNAs without the entire 5′ UTR were defective for both efficient packaging and propagation into target cells.

Conclusions/Significance

These results reveal that the 5′ end of MMTV genome is critical for both gRNA packaging and propagation, unlike the recently delineated FIV and MPMV packaging determinants that have been shown to be of bipartite nature.     

Challenging the Roles of NSP3 and Untranslated Regions in Rotavirus mRNA Translation

Rotavirus NSP3 is a translational surrogate of the PABP-poly(A) complex for rotavirus mRNAs. To further explore the effects of NSP3 and untranslated regions (UTRs) on rotavirus mRNAs translation, we used a quantitative in vivo assay with simultaneous cytoplasmic NSP3 expression (wild-type or deletion mutant) and electroporated rotavirus-like and standard synthetic mRNAs. This assay shows that the last four GACC nucleotides of viral mRNA are essential for efficient translation and that both the NSP3 eIF4G- and RNA-binding domains are required. We also show efficient translation of rotavirus-like mRNAs even with a 5’UTR as short as 5 nucleotides, while more than eleven nucleotides are required for the 3’UTR. Despite the weak requirement for a long 5’UTR, a good AUG environment remains a requirement for rotavirus mRNAs translation.     

RNA Polymerase Activity and Specific RNA Structure Are Required for Efficient HCV Replication in Cultured Cells

We have previously reported that the NS3 helicase (N3H) and NS5B-to-3′X (N5BX) regions are important for the efficient replication of hepatitis C virus (HCV) strain JFH-1 and viral production in HuH-7 cells. In the current study, we investigated the relationships between HCV genome replication, virus production, and the structure of N5BX. We found that the Q377R, A450S, S455N, R517K, and Y561F mutations in the NS5B region resulted in up-regulation of J6CF NS5B polymerase activity in vitro. However, the activation effects of these mutations on viral RNA replication and virus production with JFH-1 N3H appeared to differ. In the presence of the N3H region and 3′ untranslated region (UTR) of JFH-1, A450S, R517K, and Y561F together were sufficient to confer HCV genome replication activity and virus production ability to J6CF in cultured cells. Y561F was also involved in the kissing-loop interaction between SL3.2 in the NS5B region and SL2 in the 3′X region. We next analyzed the 3′ structure of HCV genome RNA. The shorter polyU/UC tracts of JFH-1 resulted in more efficient RNA replication than J6CF. Furthermore, 9458G in the JFH-1 variable region (VR) was responsible for RNA replication activity because of its RNA structures. In conclusion, N3H, high polymerase activity, enhanced kissing-loop interactions, and optimal viral RNA structure in the 3′UTR were required for J6CF replication in cultured cells.     

The Footprint of Genome Architecture in the Largest Genome Expansion in RNA Viruses

The small size of RNA virus genomes (2-to-32 kb) has been attributed to high mutation rates during replication, which is thought to lack proof-reading. This paradigm is being revisited owing to the discovery of a 3′-to-5′ exoribonuclease (ExoN) in nidoviruses, a monophyletic group of positive-stranded RNA viruses with a conserved genome architecture. ExoN, a homolog of canonical DNA proof-reading enzymes, is exclusively encoded by nidoviruses with genomes larger than 20 kb. All other known non-segmented RNA viruses have smaller genomes. Here we use evolutionary analyses to show that the two- to three-fold expansion of the nidovirus genome was accompanied by a large number of replacements in conserved proteins at a scale comparable to that in the Tree of Life. To unravel common evolutionary patterns in such genetically diverse viruses, we established the relation between genomic regions in nidoviruses in a sequence alignment-free manner. We exploited the conservation of the genome architecture to partition each genome into five non-overlapping regions: 5′ untranslated region (UTR), open reading frame (ORF) 1a, ORF1b, 3′ORFs (encompassing the 3′-proximal ORFs), and 3′ UTR. Each region was analyzed for its contribution to genome size change under different models. The non-linear model statistically outperformed the linear one and captured >92% of data variation. Accordingly, nidovirus genomes were concluded to have reached different points on an expansion trajectory dominated by consecutive increases of ORF1b, ORF1a, and 3′ORFs. Our findings indicate a unidirectional hierarchical relation between these genome regions, which are distinguished by their expression mechanism. In contrast, these regions cooperate bi-directionally on a functional level in the virus life cycle, in which they predominantly control genome replication, genome expression, and virus dissemination, respectively. Collectively, our findings suggest that genome architecture and the associated region-specific division of labor leave a footprint on genome expansion and may limit RNA genome size.     

Major alteration in coxsackievirus B3 genomic RNA structure distinguishes a virulent strain from an avirulent strain

Coxsackievirus B3 (CV-B3) is a cardiovirulent enterovirus that utilizes a 5′ untranslated region (5′UTR) to complete critical viral processes. Here, we directly compared the structure of a 5′UTR from a virulent strain with that of a naturally occurring avirulent strain. Using chemical probing analysis, we identified a structural difference between the two 5′UTRs in the highly substituted stem-loop II region (SLII). For the remainder of the 5′UTR, we observed conserved structure. Comparative sequence analysis of 170 closely related enteroviruses revealed that the SLII region lacks conservation. To investigate independent folding and function, two chimeric CV-B3 strains were created by exchanging nucleotides 104–184 and repeating the 5′UTR structural analysis. Neither the parent SLII nor the remaining domains of the background 5′UTR were structurally altered by the exchange, supporting an independent mechanism of folding and function. We show that the attenuated 5′UTR lacks structure in the SLII cardiovirulence determinant.     

Genome-wide analysis of single non-templated nucleotides in plant endogenous siRNAs and miRNAs

Plant small RNAs are subject to various modifications. Previous reports revealed widespread 3′ modifications (truncations and non-templated tailing) of plant miRNAs when the 2′-O-methyltransferase HEN1 is absent. However, non-templated nucleotides in plant heterochromatic siRNAs have not been deeply studied, especially in wild-type plants. We systematically studied non-templated nucleotide patterns in plant small RNAs by analyzing small RNA sequencing libraries from Arabidopsis, tomato, Medicago, rice, maize and Physcomitrella. Elevated rates of non-templated nucleotides were observed at the 3′ ends of both miRNAs and endogenous siRNAs from wild-type specimens of all species. ‘Off-sized’ small RNAs, such as 25 and 23 nt siRNAs arising from loci dominated by 24 nt siRNAs, often had very high rates of 3′-non-templated nucleotides. The same pattern was observed in all species that we studied. Further analysis of 24 nt siRNA clusters in Arabidopsis revealed distinct patterns of 3′-non-templated nucleotides of 23 nt siRNAs arising from heterochromatic siRNA loci. This pattern of non-templated 3′ nucleotides on 23 nt siRNAs is not affected by loss of known small RNA 3′-end modifying enzymes, and may result from modifications added to longer heterochromatic siRNA precursors.     

Genome 3'-end repair in dengue virus type 2

Genomes of RNA viruses encounter a continual threat from host cellular ribonucleases. Therefore, viruses have evolved mechanisms to protect the integrity of their genomes. To study the mechanism of 3′-end repair in dengue virus-2 in mammalian cells, a series of 3′-end deletions in the genome were evaluated for virus replication by detection of viral antigen NS1 and by sequence analysis. Limited deletions did not cause any delay in the detection of NS1 within 5 d. However, deletions of 7–10 nucleotides caused a delay of 9 d in the detection of NS1. Sequence analysis of RNAs from recovered viruses showed that at early times, virus progenies evolved through RNA molecules of heterogeneous lengths and nucleotide sequences at the 3′ end, suggesting a possible role for terminal nucleotidyl transferase activity of the viral polymerase (NS5). However, this diversity gradually diminished and consensus sequences emerged. Template activities of 3′-end mutants in the synthesis of negative-strand RNA in vitro by purified NS5 correlate well with the abilities of mutant RNAs to repair and produce virus progenies. Using the Mfold program for RNA structure prediction, we show that if the 3′ stem–loop (3′ SL) structure was abrogated by mutations, viruses eventually restored the 3′ SL structure. Taken together, these results favor a two-step repair process: non-template-based nucleotide addition followed by evolutionary selection of 3′-end sequences based on the best-fit RNA structure that can support viral replication.     

Improved Annotation of 3′ Untranslated Regions and Complex Loci by Combination of Strand-Specific Direct RNA Sequencing,RNA-Seq and ESTs

The reference annotations made for a genome sequence provide the framework for all subsequent analyses of the genome. Correct and complete annotation in addition to the underlying genomic sequence is particularly important when interpreting the results of RNA-seq experiments where short sequence reads are mapped against the genome and assigned to genes according to the annotation. Inconsistencies in annotations between the reference and the experimental system can lead to incorrect interpretation of the effect on RNA expression of an experimental treatment or mutation in the system under study. Until recently, the genome-wide annotation of 3′ untranslated regions received less attention than coding regions and the delineation of intron/exon boundaries. In this paper, data produced for samples in Human, Chicken and A. thaliana by the novel single-molecule, strand-specific, Direct RNA Sequencing technology from Helicos Biosciences which locates 3′ polyadenylation sites to within +/− 2 nt, were combined with archival EST and RNA-Seq data. Nine examples are illustrated where this combination of data allowed: (1) gene and 3′ UTR re-annotation (including extension of one 3′ UTR by 5.9 kb); (2) disentangling of gene expression in complex regions; (3) clearer interpretation of small RNA expression and (4) identification of novel genes. While the specific examples displayed here may become obsolete as genome sequences and their annotations are refined, the principles laid out in this paper will be of general use both to those annotating genomes and those seeking to interpret existing publically available annotations in the context of their own experimental data.     

The 3' proximal translational enhancer of Turnip crinkle virus binds to 60S ribosomal subunits

During cap-dependent translation of eukaryotic mRNAs, initiation factors interact with the 5′ cap to attract ribosomes. When animal viruses translate in a cap-independent fashion, ribosomes assemble upstream of initiation codons at internal ribosome entry sites (IRES). In contrast, many plant viral genomes do not contain 5′ ends with substantial IRES activity but instead have 3′ translational enhancers that function by an unknown mechanism. A 393-nucleotide (nt) region that includes the entire 3′ UTR of the Turnip crinkle virus (TCV) synergistically enhances translation of a reporter gene when associated with the TCV 5′ UTR. The major enhancer activity was mapped to an internal region of ~140 nt that partially overlaps with a 100-nt structural domain previously predicted to adopt a form with some resemblance to a tRNA, according to a recent study by J.C. McCormack and colleagues. The T-shaped structure binds to 80S ribosomes and 60S ribosomal subunits, and binding is more efficient in the absence of surrounding sequences and in the presence of a pseudoknot that mimics the tRNA-acceptor stem. Untranslated TCV satellite RNA satC, which contains the TCV 3′ end and 6-nt differences in the region corresponding to the T-shaped element, does not detectably bind to 80S ribosomes and is not predicted to form a comparable structure. Binding of the TCV T-shaped element by 80S ribosomes was unaffected by salt-washing, reduced in the presence of AcPhe-tRNA, which binds to the P-site, and enhanced binding of Phe-tRNA to the ribosome A site. Mutations that reduced translation in vivo had similar effects on ribosome binding in vitro. This strong correlation suggests that ribosome entry in the 3′ UTR is a key function of the 3′ translational enhancer of TCV and that the T-shaped element contains some tRNA-like properties.