Evolution and Taxonomy of Positive-Strand RNA Viruses: Implications of Comparative Analysis of Amino Acid Sequences

Abstract

Despite the rapid mutational change that is typical of positive-strand RNA viruses, enzymes mediating the replication and expression of virus genomes contain arrays of conserved sequence motifs. Proteins with such motifs include RNA-dependent RNA polymerase, putative RNA helicase, chymotrypsin-like and papain-like proteases, and methyltransferases. The genes for these proteins form partially conserved modules in large subsets of viruses. A concept of the virus genome as a relatively evolutionarily stable “core” of housekeeping genes accompanied by a much more flexible “shell” consisting mostly of genes coding for virion components and various accessory proteins is discussed. Shuffling of the “shell” genes including genome reorganization and recombination between remote groups of viruses is considered to be one of the major factors of virus evolution.

Multiple alignments for the conserved viral proteins were constructed and used to generate the respective phylogenetic trees. Based primarily on the tentative phylogeny for the RNA-dependent RNA polymerase, which is the only universally conserved protein of positive-strand RNA viruses, three large classes of viruses, each consisting of distinct smaller divisions, were delineated. A strong correlation was observed between this grouping and the tentative phylogenies for the other conserved proteins as well as the arrangement of genes encoding these proteins in the virus genome. A comparable correlation with the polymerase phylogeny was not found for genes encoding virion components or for genome expression strategies. It is surmised that several types of arrangement of the “shell” genes as well as basic mechanisms of expression could have evolved independently in different evolutionary lineages.

The grouping revealed by phylogenetic analysis may provide the basis for revision of virus classification, and phylogenetic taxonomy of positive-strand RNA viruses is outlined. Some of the phylogenetically derived divisions of positive-strand RNA viruses also include double-stranded RNA viruses, indicating that in certain cases the type of genome nucleic acid may not be a reliable taxonomic criterion for viruses.

Hypothetical evolutionary scenarios for positive-strand RNA viruses are proposed. It is hypothesized that all positive-strand RNA viruses and some related double-stranded RNA viruses could have evolved from a common ancestor virus that contained genes for RNA-dependent RNA polymerase, a chymotrypsin-related protease that also functioned as the capsid protein, and possibly an RNA helicase.     

Molecular dissection of influenza virus nucleoprotein: deletion mapping of the RNA binding domain.

Influenza virus nucleoprotein (NP) is associated with the genome RNA, forming ribonucleoprotein cores. To identify the amino acid sequence involved in RNA binding, we performed Northwestern blot analysis with a set of N- and C-terminal deletion mutants of NP produced in Escherichia coli. The RNA binding region has been mapped between amino acid residues 91 and 188, a stretch of residues that contains a sequence that is not only highly conserved among NPs from A-, B-, and C-type influenza viruses but also similar to the RNA binding domain of a plant virus movement protein.     

Evolution of RNA-dependent RNA polymerases of positive riboviruses

A comparative analysis is presented of 24 known amino acid sequences of RNA-dependent RNA polymerases of positive strand RNA viruses infecting animals, plants and bacteria. Using a newly proposed methodology of group alignment for weakly similar sequences, evolutionary conserved fragments of all these proteins were unambiguously aligned. A unique pattern (consensus) of 7 invariant amino acid residues was revealed which is absent from the sequences of other RNA and DNA polymerases and is thought to unequivocally identify the RNA-dependent RNA polymerases of positive strand RNA viruses. Based on the obtained alignment a tentative phylogenetic tree of viral RNA polymerases was constructed for the first time. The RNA-dependent RNA polymerases of positive strand RNA viruses are concluded to comprise a distinct family of evolutionary related proteins.     

Relationships among the positive strand and double-strand RNA viruses as viewed through their RNA-dependent RNA polymerases.

The sequences of 50 RNA-dependent RNA polymerases (RDRPs) from 43 positive strand and 7 double strand RNA (dsRNA) viruses have been compared. The alignment permitted calculation of distances among the 50 viruses and a resultant dendrogram based on every amino acid, rather than just those amino acids in the conserved motifs. Remarkably, a large subgroup of these viruses, including vertebrate, plant, and insect viruses, forms a single cluster whose only common characteristic is exploitation of insect hosts or vectors. This similarity may be due to molecular constraints associated with a present and/or past ability to infect insects and/or to common descent from insect viruses. If common descent is important, as it appears to be, all the positive strand RNA viruses of eucaryotes except for the picornaviruses may have evolved from an ancestral dsRNA virus. Viral RDRPs appear to be inherited as modules rather than as portions of single RNA segments, implying that RNA recombination has played an important role in their dissemination.     

The complete genome sequence of a single-stranded RNA virus from the tarnished plant bug, Lygus lineolaris (Palisot de Beauvois)

The complete genome sequence of a single-stranded RNA virus infecting the tarnished plant bug, Lygus lineolaris (Palisot de Beauvois), was identified by sequencing cDNA prepared from insects collected from the Mississippi Delta. The 9655 nucleotide positive-sense single-stranded RNA genome of the L. lineolaris single-stranded RNA virus (LyLV-1) contained a single open reading frame of 8958 nucleotides encoding a 2986 amino acid genome polypeptide. The open reading frame was flanked by untranslated regions of 603 and 69 nucleotides at the 5'- and 3'- ends of the genome, respectively. Database searches and homology based modeling was used to identify four capsid proteins (VP1-VP4), helicase/AAA-ATPase, cysteine protease (C3P), protease 2A, and the RNA-directed RNA polymerase (RdRp). In addition, a region with weak similarity to the eukaryotic structural maintenance of chromosome (SMC) domain was identified near the amino-terminal of the polyprotein and adjacent to the VP1 domain. The amino acid sequence of LyLV-1 was approximately 44.4% similar to that of sacbrood virus (SBV) of the honey bee. The genomic organization of both viruses showed remarkable similarity with the exception of highly divergent amino acid regions flanking fairly conserved structural and non-structural polypeptide regions. High similarity to the SBV genome and similarities in the genome organization and amino acid sequence with the viruses of the family Iflaviridae suggested that LyLV-1 was a novel member of this family. Virus particles were 39 nm in diameter and appeared to transmit vertically via eggs. Although this virus may only cause covert infections under normal conditions, the potential for using this virus in biological control of L. lineolaris is discussed.     

副粘病毒Tianjin株NP、P、M及L蛋白的生物信息学分析

为了进一步明确副粘病毒Tianjin株的来源和种系进化地位,探讨其高致病性的机制.对Tianjin株NP、P、M及L蛋白进行了生物信息学分析.进化树显示:Tianjin株属于副粘病毒亚科呼吸道病毒属,且很可能为仙台病毒新的基因型.相似性比较表明,P蛋白变异最大.相似性仅为78.7%～91.9%;L蛋白相似性最高,为96.0%～98.0%.序列比对显示:NP蛋白氨基酸序列中存在15个独特的变异位点,P蛋白存在29个,M蛋白存在6个,L蛋白存在29个.这些独特变异位点的存在很可能是导致Tianjin株在宿主来源和致病特点等方面与已知仙台病毒株具有较大差异的原因.     

Structural properties of carnation mottle virus p7 movement protein and its RNA-binding domain

Plant viral movement proteins (MPs) participate actively in the intra- and intercellular movement of RNA plant viruses to such an extent that MP dysfunction impairs viral infection. However, the molecular mechanism(s) of their interaction with cognate nucleic acids are not well understood, partly due to the lack of structural information. In this work, a protein dissection approach was used to gain information on the structural and RNA-binding properties of this class of proteins, as exemplified by the 61-amino acid residue p7 MP from carnation mottle virus (CarMV). Circular dichroism spectroscopy showed that CarMV p7 is an alpha/beta RNA-binding soluble protein. Using synthetic peptides derived from the p7 sequence, we have identified three distinct putative domains within the protein. EMSA showed that the central region, from residue 17 to 35 (represented by peptide p7(17-35)), is responsible for the RNA binding properties of CarMV p7. This binding peptide populates a nascent alpha-helix in water solution that is further stabilized in the presence of either secondary structure inducers, such as trifluoroethanol and monomeric SDS, or RNA (which also changes its conformation upon binding to the peptide). Thus, the RNA recognition appears to occur via an "adaptive binding" mechanism. Interestingly, the amino acid sequence and structural properties of the RNA-binding domain of p7 seem to be conserved among carmoviruses and some other RNA-binding proteins and peptides. The low conserved N terminus of p7 (peptide p7(1-16)) is unstructured in solution. In contrast, the highly conserved C terminus motif (peptide p7(40-61)) adopts a beta-sheet conformation in aqueous solution. Alanine scanning mutagenesis of the RNA-binding motif showed how selected positive charged amino acids are more relevant than others in the RNA binding process and how hydrophobic amino acid side chains would participate in the stabilization of the protein-RNA complex.     

A new superfamily of putative NTP-binding domains encoded by genomes of small DNA and RNA viruses

Statistically significant similarity was revealed between amino acid sequences of NTP-binding pattern-containing domains which are among the most conserved protein segments in dissimilar groups of ss and dsDNA viruses (papova-, parvo-, geminiviruses and P4 bacteriophage), and RNA viruses (picorna-, como- and nepoviruses) with small genomes. Within the aligned domains of 100-120 amino acid residues, three highly conserved sequence segments have been identified, i.e. 'A' and 'B' motifs of the NTP-binding pattern, and a third, C-terminal motif 'C', not described previously. The sequence of the 'B' motif in the proteins of the new superfamily is unusually variable, with substitutions, in some of the members, of the Asp residue conserved in other NTP-binding proteins. The 'C' motif is characterized by an invariant Asn residue preceded by a stretch of hydrophobic residues. As the new superfamily included a well studied DNA and RNA helicase, T antigen of SV40, helicase function could be tentatively assigned also to the other related viral putative NTP-binding proteins. On the other hand, the possibility of different and/or multiple functions for some of these proteins is discussed.     

Manganese-dependent polioviruses caused by mutations within the viral polymerase

Viral RNA-dependent RNA polymerases exhibit great sequence diversity. Only six core amino acids are conserved across all polymerases of positive-strand RNA viruses of eukaryotes. While exploring the function of one of these completely conserved residues, asparagine 297 in the prototypic poliovirus polymerase 3D(pol), we identified three viable mutants with noncanonical amino acids at this conserved position. Although asparagine 297 could be replaced by glycine or alanine in these mutants, the viruses exhibited Mn(2+)-dependent RNA replication and viral growth. All known RNA polymerases and replicative polymerases of bacterial, eukaryotic, and viral organisms are thought to be magnesium dependent in vivo, and therefore these mutant polioviruses may represent the first viruses with a requirement for an alternative polymerase cation. These results demonstrate the extreme functional flexibility of viral RNA-dependent RNA polymerases. Furthermore, the finding that strictly conserved residues in the nucleotide binding pocket of the polymerase can be altered in a manner that supports virus production suggests that drugs targeting this region of the enzyme will still be susceptible to the problem of drug-resistant escape mutants.     

The protein family of pyruvate:quinone oxidoreductases: Amino acid sequence conservation and taxonomic distribution

Pyruvate:quinone oxidoreductases (PQOs) catalyse the oxidative decarboxylation of pyruvate to acetate and concomitant reduction of quinone to quinol with the release of CO₂. They are thiamine pyrophosphate (TPP) and flavin-adenine dinucleotide (FAD) containing enzymes, which interact with the membrane in a monotopic way. PQOs are considered as part of alternatives to most recognized pyruvate catabolizing pathways, and little is known about their taxonomic distribution and structural/functional relationship.In this bioinformatics work we tackled these gaps in PQO knowledge. We used the KEGG database to identify PQO coding genes, performed a multiple sequence analysis which allowed us to study the amino acid conservation on these enzymes, and looked at their possible cellular function. We observed that PQOS are enzymes exclusively present in prokaryotes with most of the sequences identified in bacteria. Regarding the amino acid sequence conservation, we found that 75 amino acid residues (out of 570, on average) have a conservation over 90 %, and that the most conserved regions in the protein are observed around the TPP and FAD binding sites. We systematized the presence of conserved features involved in Mg²⁺, TPP and FAD binding, as well as residues directly linked to the catalytic mechanism. We also established the presence of a new motif named “HEH lock”, possibly involved in the dimerization process. The results here obtained for the PQO protein family contribute to a better understanding of the biochemistry of these respiratory enzymes.     

Structural homology among four nodaviruses as deduced by sequencing and X-ray crystallography

The genomic RNA2s of nodaviruses encode a single gene, that of protein alpha, the precursor of virion proteins beta and gamma. We compared the sequences of the RNA2s of the nodaviruses, black beetle virus (BBV), flock house virus, boolarra virus and nodamura virus, with the objective of identifying homologies in the primary and secondary structure of these RNAs and in the structure of their encoded protein. The sequences of the four RNAs were found to be similar, so that homologous regions relating to translation and RNA replication were readily identified. However, the overall, secondary structures in solution, deduced from calculations of optimal Watson-Crick base-pairing configurations, were very different for the four RNAs. We conclude that a particular, overall, secondary structure in solution within host cells is not required for virus viability. The partially refined X-ray structure of BBV (R = 26.4% for the current model) was used as a framework for comparing the structure of the encoded proteins of the four viruses. Mapping of the four protein sequences onto the BBV capsid showed many amino acid differences on the outer surface, indicating that the exteriors of the four virions are substantially different. Mapping in the beta-barrel region showed an intermediate level of differences, indicating that some freedom in choice of amino acid residues is possible there although the basic framework of the capsids is evidently conserved. Mapping onto the interior surface of the BBV capsid showed a high degree of conservation of amino acid residues, particularly near the protein cleavage site, implying that that region is nearly identical in all four virions and has an essential role in virion maturation, and also suggests that all four capsid interior surfaces have similar surfaces exposed to the viral RNA. Apart from a small portion of the C promoter, the amino terminus of the BBV protein (residues 1 to 60) is crystallographically disordered and the amino acid residues in that region are not well conserved. The disordered portion of the BBV protein clearly projects from the capsid inner surface into the interior of the virion, the region occupied by the viral RNA. In all four viruses, residues 1 to 60 had a high proportion of basic residues, suggesting a virus-specific interaction of the amino terminus with the virion RNA.     

Overlapping open reading frames revealed by complete nucleotide sequencing of turnip yellow mosaic virus genomic RNA.

The complete nucleotide sequence of turnip yellow mosaic virus (TYMV) genomic RNA has been determined on a set of overlapping cDNA clones using a sequential sequencing strategy. The RNA is 6318 nucleotides long, excluding the cap structure. The genome organization deduced from the sequence confirms previous results of in vitro translation. A novel open reading frame (ORF) putatively encoding a Pro-rich and very basic 69K (K = kilodalton) protein is detected at the 5' end of the genome. It is initiated at the first AUG codon on the RNA and overlaps the major ORF that encodes the non structural 206K (previously referred to as 195K) protein of TYMV; its function is unknown. Several amino acid consensus sequences already described among plant and animal viruses are also found in the TYMV-encoded polypeptides. A comparison with other viruses whose RNA sequence is known leads to the conclusion that TYMV belongs to the "Sindbis-like" supergroup of viruses and could be related to Semliki forest virus.     

Identification and classification of a resistance breaking strain of barley yellow mosaic virus

A “resistance breaking” isolate of barley yellow mosaic virus-2 (BaYMV-2) was cloned as a cDNA and characterised. Restriction mapping and comparison with a German and a Japanese isolate of BaYMV (BaYMV-G and BaYMV-J) revealed a high level of restriction site conservation for RNAl and the majority of RNA2. However, in a box of approximately 600 nucleotides located on RNA2, striking differences in the restriction pattern could be identified. The nucleotide sequence of this box, as well as of the 3‘-terminal region of RNA1 including the coat protein coding region, and the deduced amino acid sequences were analysed. Identity at the amino acid level was between 99.3% and 92.3% in comparison with the corresponding sequences of BaYMV-G and BaYMV-J, suggesting that BaYMV-2 is closely related to BaYMV. Consequently, the classification of BaYMV-2 as a “resistance breaking” strain of BaYMV is justified.     

利用小RNA深度测序技术检测灰飞虱病毒种类

灰飞虱是一种重要农业害虫,作为病毒介体,可以传播多种植物病毒引起水稻、小麦和玉米等粮食作物病毒病害。目前对于灰飞虱体内昆虫病毒种类尚缺乏系统认识,难以开展利用昆虫病毒防治灰飞虱相关研究工作。为挖掘灰飞虱体内昆虫病毒资源,本文通过小RNA深度测序技术对灰飞虱所携带的病毒种类进行分析鉴定。结果显示,测序数据比对得到13种病毒,涉及8个病毒科和2种未分类病毒。除占据优势的水稻条纹病毒外,其余均为专性寄生的昆虫病毒,包括5种正链RNA病毒、2种单链DNA病毒和5种双链DNA病毒。在这些病毒中,发现了一种与果蝇A病毒较相似的新病毒,克隆出其依赖RNA的RNA聚合酶(RdRP)基因1-1 932位核酸序列,经Blast比对和系统进化分析,在氨基酸序列中发现RdRP掌型亚结构域保守区呈现复制酶置换四体病毒科(Permutotetraviridae)病毒所具有的“C-A-B”排列样式,确定该病毒是一种新的类复制酶置换四体病毒,暂命名为Laodelphax striatellus permutotetra-like virus(LsPLV)。这是首次在半翅目昆虫中发现类复制酶置换四体病毒。本研究表明灰飞...     

Correlation of co-ordinated amino acid substitutions with function in viruses related to tobacco mosaic virus

Sequence data are available for the coat proteins of seven tobamoviruses, with homologies ranging from at least 26% to 82%, and atomic co-ordinates are known for tobacco mosaic virus (TMV) vulgare. A significant spatial relationship has been found between groups of residues with identical amino acid substitution patterns. This strongly suggest that their location is linked to a particular function, at least in viruses identical with the wild-type for these residues. The most conserved feature of TMV is the RNA binding region. Core residues are conserved in all viruses or show mutations complementary in volume. The specificity of inter-subunit contacts is achieved in different ways in the three more distantly related viruses.