Genome RNA terminus conservation and diversity among vesiculoviruses

Construction & establishment of two minigenome rescue systems for Chandipura virus driven by recombinant vaccinia virus expressing T7 polymerase

Background & objectives:

Chandipura virus (CHPV) is an emerging pathogenic rhabdovirus with a high case fatality rate. There are no reports of a minigenome system for CHPV, which could help its study without having to use the infectious agent. This study was, therefore, undertaken for the establishment of T7 polymerase-driven minigenome system for CHPV.

Methods:

The minigenome rescue system for CHPV consists of three helper plasmids expressing the nucleocapsid protein (N), phosphoprotein (P) and large protein (L) based on a recombinant vaccinia virus expressing bacteriophage T7 polymerase (vTF7-3). The minigenome construct is composed of a reporter gene, flanked by the non-coding regions of CHPV. Two minigenomes were constructed in an antigenome or complimentary sense, expressing luciferase or green fluorescent protein (GFP). The minigenome system was evaluated by co-transfection of the minigenome construct and three helper plasmids into CV-1 cells and analysis of the reporter gene activity.

Results:

All the helper proteins were expressed from the helper plasmids confirmed by Western blotting. Expression of reporter genes was observed from both the GFP and luciferase-based minigenomes. Green fluorescence could be visualized directly in live CV-1 cells. Luciferase activity was found to be significantly different from control.

Interpretation & conclusions:

The results showed that the helper plasmids provided all the necessary viral structural proteins required for the production of minigenome mRNA template, which in turn could rescue the expression of reporter genes. Thus, these minigenomes can be applied to mimic the manifestation of CHPV life cycle.

Complete genome sequence of Piry vesiculovirus

Piry virus (PIRYV) is a rhabdovirus (genus Vesiculovirus) and is described as a possible human pathogen, originally isolated from a Philander opossum trapped in Para State, Northern Brazil. This study describes the complete full coding sequence and the genetic characterization of PIRYV. The genome sequence reveals that PIRYV has a typical vesiculovirus-like organization, encoding the five genes typical of the genus. Phylogenetic analysis confirmed that PIRYV is most closely related to Perinet virus and clustered in the same clade as Chandipura and Isfahan vesiculoviruses.

Chandipura Virus: an emerging tropical pathogen

Chandipura Virus (CHPV), a member of Rhabdoviridae, is responsible for an explosive outbreak in rural areas of India. It affects mostly children and is characterized by influenza-like illness and neurologic dysfunctions. It is transmitted by vectors such as mosquitoes, ticks and sand flies. An effective real-time one step reverse-transcriptase PCR assay method is adopted for diagnosis of this virus. CHPV has a negative sense RNA genome encoding five different proteins (N, P, M, G, and L). P protein plays a vital role in the virus's life cycle, while M protein is lethal in nature. There is no specific treatment available to date, symptomatic treatment involves use of mannitol to reduce brain edema. A Vero cell based vaccine candidate against CHPV was evaluated efficiently as a preventive agent against it. Prevention is the best method to suppress CHPV infection. Containment of disease transmitting vectors, maintaining good nutrition, health, hygiene and awareness in rural areas will help in curbing the menace of CHPV. Thus, to control virus transmission some immense preventive measures need to be attempted until a good anti-CHPV agent is developed.

Whole genomes of Chandipura virus isolates and comparative analysis with other rhabdoviruses

The Chandipura virus (CHPV) belonging to the Vesiculovirus genus and Rhabdoviridae family, has recently been associated with a number of encephalitis epidemics, with high mortality in children, in different parts of India. No full length genome sequences of CHPV isolates were available in GenBank and little is known about the molecular markers for pathogenesis. In the present study, we provide the complete genomic sequences of four isolates from epidemics during 2003-2007. These sequences along with the deduced sequence of the prototype isolate of 1965 were analysed using phylogeny, motif search, homology modeling and epitope prediction methods. Comparison with other rhaboviruses was also done for functional extrapolations. All CHPV isolates clustered with the Isfahan virus and maintained several functional motifs of other rhabdoviruses. A notable difference with the prototype vesiculovirus, Vesicular Stomatitis Virus was in the L-domain flanking sequences of the M protein that are known to be crucial for interaction with host proteins. With respect to the prototype isolate, significant additional mutations were acquired in the 2003-2007 isolates. Several mutations in G mapped onto probable antigenic sites. A mutation in N mapped onto regions crucial for N-N interaction and a putative T-cell epitope. A mutation in the Casein kinase II phosphorylation site in P may attribute to increased rates of phosphorylation. Gene junction comparison revealed changes in the M-G junction of all the epidemic isolates that may have implications on read-through and gene transcription levels. The study can form the basis for further experimental verification and provide additional insights into the virulence determinants of the CHPV.

Genetic and antigenic relationships of vesicular stomatitis viruses from South America

Vesicular stomatitis (VS) viruses have been classified into two serotypes: New Jersey (VSNJV) and Indiana (VSIV). Here, we have characterized field isolates causing vesicular stomatitis in Brazil and Argentina over a 35-year span. Cluster analysis based on either serological relatedness, as inferred from virus neutralization and complement fixation assays, or nucleotide sequences of two separate genes (phosphoprotein or glycoprotein) grouped the field isolates into two distinct monophyletic groups within the Indiana serogroup. One group included seven viruses from Brazil and Argentina that were serologically classified as Indiana-2 and Cocal virus (COCV). The other group contained three viruses from Brazil that were serologically classified as Indiana-3 and the prototype of this group, Alagoas virus (VSAV). Interestingly, two vesiculoviruses that were isolated from insects but do not cause disease in animals, one from Brazil (Maraba virus; MARAV) and the other from Colombia (CoAr 171638), grouped into two separate genetic lineages within the Indiana serotype. Our data provide support for the classification of viruses causing clinical VS in livestock in Brazil and Argentina into two distinct groups: Indiana-2 (VSIV-2) and Indiana-3 (VSIV-3). We suggest using nomenclature for these viruses that includes the serotype, year and place of occurrence, and affected host. This nomenclature is consistent with that currently utilized to describe field isolates of VSNJV or VSIV in scientific literature.

Genome of turbot rhabdovirus exhibits unusual non-coding regions and an additional ORF that could be expressed in fish cell

Genomic sequence of Scophthalmus maximus rhabdovirus (SMRV) isolated from diseased turbot has been characterized. The complete genome of SMRV comprises 11,492 nucleotides and encodes five typical rhabdovirus genes N, P, M, G and L. In addition, two open reading frames (ORF) are predicted overlapping with P gene, one upstream of P and smaller than P (temporarily called Ps), and another in P gene which may encodes a protein similar to the vesicular stomatitis virus C protein. The C ORF is contained within the P ORF. The five typical proteins share the highest sequence identities (48.9%) with the corresponding proteins of rhabdoviruses in genus Vesiculovirus. Phylogenetic analysis of partial L protein sequence indicates that SMRV is close to genus Vesiculovirus. The first 13 nucleotides at the ends of the SMRV genome are absolutely inverse complementarity. The gene junctions between the five genes show conserved polyadenylation signal (CATGA(7)) and intergenic dinucleotide (CT) followed by putative transcription initiation sequence A(A/G)(C/G)A(A/G/T), which are different from known rhabdoviruses. The entire Ps ORF was cloned and expressed, and used to generate polyclonal antibody in mice. One obvious band could be detected in SMRV-infected carp leucocyte cells (CLCs) by anti-Ps/C serum via Western blot, and the subcellular localization of Ps-GFP fusion protein exhibited cytoplasm distribution as multiple punctuate or doughnut shaped foci of uneven size.

Reviewing Chandipura: a vesiculovirus in human epidemics

Chandipura virus, a member of the rhabdoviridae family and vesiculovirus genera, has recently emerged as human pathogen that is associated with a number of outbreaks in different parts of India. Although, the virus closely resembles with the prototype vesiculovirus, Vesicular Stomatitis Virus, it could be readily distinguished by its ability to infect humans. Studies on Chandipura virus while shed light into distinct stages of viral infection; it may also allow us to identify potential drug targets for antiviral therapy. In this review, we have summarized our current understanding of Chandipura virus life cycle at the molecular detail with particular interest in viral RNA metabolisms, namely transcription, replication and packaging of viral RNA into nucleocapsid structure. Contemporary research on otherwise extensively studied family member Vesicular Stomatitis Virus has also been addressed to present a more comprehensive picture of vesiculovirus life cycle. Finally, we reveal examples of protein economy in Chandipura virus life-cycle whereby each viral protein has evolved complexity to perform multiple tasks.

Utilization of homotypic and heterotypic proteins of vesicular stomatitis virus by defective interfering particle genomes for RNA replication and virion assembly: implications for the mechanism of homologous viral interference

Defective interfering (DI) particles of Indiana serotype of vesicular stomatitis virus (VSV(Ind)) are capable of interfering with the replication of both homotypic VSV(Ind) and heterotypic New Jersey serotype (VSV(NJ)) standard virus. In contrast, DI particles from VSV(NJ) do not interfere with the replication of VSV(Ind) standard virus but do interfere with VSV(NJ) replication. The differences in the interfering activities of VSV(Ind) DI particles and VSV(NJ) DI particles against heterotypic standard virus were investigated. We examined the utilization of homotypic and heterotypic VSV proteins by DI particle genomic RNAs for replication and maturation into infectious DI particles. Here we show that the RNA-nucleocapsid protein (N) complex of one serotype does not utilize the polymerase complex (P and L) of the other serotype for RNA synthesis, while DI particle genomic RNAs of both serotypes can utilize the N, P, and L proteins of either serotype without serotypic restriction but with differing efficiencies as long as all three proteins are derived from the same serotype. The genomic RNAs of VSV(Ind) DI particles assembled and matured into DI particles by using either homotypic or heterotypic viral proteins. In contrast, VSV(NJ) DI particles could assemble only with homotypic VSV(NJ) viral proteins, although the genomic RNAs of VSV(NJ) DI particles could be replicated by using heterotypic VSV(Ind) N, P, and L proteins. Thus, we concluded that both efficient RNA replication and assembly of DI particles are required for the heterotypic interference by VSV DI particles.

Monomer and dimer of Chandipura virus unphosphorylated P-protein binds leader RNA differently: implications for viral RNA synthesis

Interaction of the leader RNA with the unphosphorylated P-protein has been proposed to play a key role in the transcription-replication transition of Chandipura virus, a model rhabdovirus. Electrophoretic mobility shift assay with the leader RNA and the unphosphorylated P-protein demonstrated existence of two distinct complexes in vitro. Measurements of stoichiometry indicate the protein monomer/RNA ratio to be 1:1 and 2:1 for faster and slower migrating bands, respectively. We have also observed a concentration-dependent oligomerization of the unphosphorylated P-protein, in sub-micromolar to low micromolar range. Sedimentation velocity, dynamic light scattering and large zone gel filtration experiments suggest a monomer-dimer-tetramer model of association. RNA binding experiments suggest that the two complexes assembled from one molecule of the leader RNA binding to either a protein monomer or a dimer. A truncated RNA consisting of a 3' region of the leader transcript exclusively formed the 1:1 complex, whereas a RNA consisting of only the 5' region forms the 2:1 complex exclusively. RNA binding experiments at different protein concentrations suggest that binding of the RNA comprising the 3' region weakens significantly at higher P(0) concentrations, whereas in contrast the binding of the RNA comprising the 5' region becomes modestly tighter. Implications of two different types of leader RNA-P-protein complexes in viral RNA synthesis are discussed.

Leader RNA binding ability of Chandipura virus P protein is regulated by its phosphorylation status: a possible role in genome transcription-replication switch

The molecular events associated with the transcriptive and replicative cycle of negative-stranded RNA viruses are still an enigma. We took Chandipura virus, a member of the Rhabdoviridae family, as our model system to demonstrate that Phosphoprotein P, besides Nucleocapsid protein N, also acts as a leader RNA-binding protein in its unphosphorylated form, whereas CKII-mediated phosphorylation totally abrogates its RNA-binding ability. However, interaction between P protein and leader RNA can be distinguished from N-mediated encapsidation of viral sequences. Furthermore, P protein bound to leader chain can successively recruit N protein on RNA while itself being replaced. We also observed that the accumulation of phosphorylation null mutant of P protein in cells results in enhanced genome RNA replication with concurrent increase in the viral yield. All these results led us to propose a model explaining viral transcription-replication switch where Phosphoprotein P acts as a modulator of genome transcription and replication by its ability to bind to the nascent leader RNA in its unphosphorylated form, promoting read-through of the transcription termination signals and initiating nucleocapsid assembly on the nascent RNA chain.

Determination of the complete genomic sequence and analysis of the gene products of the virus of Spring Viremia of Carp, a fish rhabdovirus

The complete genome of spring viremia of carp virus (SVCV) was cloned and the sequence of 11019 nucleotides was determined. It contains five open reading frames (ORF's) encoding for the nucleoprotein N; phosphoprotein P; matrix protein M; glycoprotein G; and the viral RNA dependent RNA polymerase L. Genes are organised in the order typical for rhabdoviruses: 3'-N-P-M-G-L-5'. The short leader and trailer regions of SVCV exhibit inverse complementarity and are similar to the respective 3' and 5' ends of the genome of vesicular stomatitis virus. To verify the predicted open reading frames proteins were expressed in bacteria and analysed with a polyclonal anti-SVCV serum. Furthermore, monospecific antisera against the distinct viral proteins were generated. Comparison of genome and protein confirm the assignment of SVCV to the genus Vesiculovirus.

Matrix protein of Chandipura virus inhibits transcription from an RNA polymerase II promoter

Chandipura virus (CHPV) is a Vesiculovirus, related to, but phylogenetically distinct from, vesicular stomatitis virus (VSV). The matrix protein of VSV, as well as its role in virus assembly, inhibits the transcription from promoters for host RNA polymerases I and II. Cloning and expression of the matrix protein of CHPV in human cells showed that this protein is also functional in its inhibitory effect on transcription of a reporter gene from the cytomegalovirus immediate-early promoter, despite sharing only 28% amino acid sequence identity with the matrix protein of VSV.

Carboxy-terminal five amino acids of the nucleocapsid protein of vesicular stomatitis virus are required for encapsidation and replication of genome RNA

The encapsidation of vesicular stomatitis virus (VSV) genome RNA, a prerequisite step to the replication process by the nucleocapsid protein (N) was studied by its ability to package VSV leader RNA in vitro in a RNase-resistant form. The VSV leader RNA was derived from the SP6 transcription vector while the N protein was made in rabbit reticulocyte lysate. The in vitro encapsidation was carried out by translating N mRNA in the presence of 32P-labeled presynthesized leader RNA. The RNA encapsidation property of the N protein was completely abrogated when the C-terminal five amino acids (VEFDK-COOH) were deleted. Systematic mutational analyses within the C-terminal five amino acid regions reveal that the RNA encapsidation activity was lost in all mutants except K --> A and K --> R, indicating that C-terminal five amino acids, in particular the lysine residue play critical role in genome RNA encapsidation. To correlate the in vitro encapsidation abilities of these mutant N proteins with genome RNA replication, we have used a full-length cDNA clone of VSV genome RNA to rescue infectious virions from cells expressing L, P, and wt or mutant N proteins and measured the recovery of plaque forming units. The results indicate that the N mutants that are defective in in vitro encapsidation of leader RNA do not support replication, establishing the requirement of C-terminal five amino acids of the N protein in viral replication.

The complete genome structure and phylogenetic relationship of infectious hematopoietic necrosis virus

Infectious hematopoietic necrosis virus (IHNV), a member of the family Rhabdoviridae, causes a severe disease with high mortality in salmonid fish. The nucleotide sequence (11,131 bases) of the entire genome was determined for the pathogenic WRAC strain of IHNV from southern Idaho. This allowed detailed analysis of all 6 genes, the deduced amino acid sequences of their encoded proteins, and important control motifs including leader, trailer and gene junction regions. Sequence analysis revealed that the 6 virus genes are located along the genome in the 3' to 5' order: nucleocapsid (N), polymerase-associated phosphoprotein (P or M1), matrix protein (M or M2), surface glycoprotein (G), a unique non-virion protein (NV) and virus polymerase (L). The IHNV genome RNA was found to have highly complementary termini (15 of 16 nucleotides). The gene junction regions display the highly conserved sequence UCURUC(U)7RCCGUG(N)4CACR (in the vRNA sense), which includes the typical rhabdovirus transcription termination/polyadenylation signal and a novel putative transcription initiation signal. Phylogenetic analysis of M, G and L protein sequences allowed insights into the evolutionary and taxonomic relationship of rhabdoviruses of fish relative to those of insects or mammals, and a broader sense of the relationship of non-segmented negative-strand RNA viruses. Based on these data, a new genus, piscivirus, is proposed which will initially contain IHNV, viral hemorrhagic septicemia virus and Hirame rhabdovirus.

Assembly and transcription of synthetic vesicular stomatitis virus nucleocapsids

The functional template for transcription of vesicular stomatitis virus (VSV) RNA is a ribonucleoprotein particle (nucleocapsid) consisting of the negative-strand sense genomic RNA completely encapsidated by the viral nucleocapsid (N) protein. As an approach to create nucleocapsids in vitro, we demonstrate here the specific encapsidation by purified N protein of in vitro-synthesized RNA sequences representing the 5' end of both the negative- and positive-strand VSV genome-length RNAs. As few as 19 nucleotides from the 5'-end of positive-strand RNA allowed maximal encapsidation, although the 5' terminal 10 nucleotides would allow partial (50%) encapsidation. Sequences downstream of the binding site can be of any origin. Specific encapsidation of VSV sequences was dependent on the presence of uninfected cell cytoplasmic extracts or poly(A). The synthetic nucleocapsids have the properties of RNase resistance and a buoyant density typical of wild-type VSV nucleocapsids. We have encapsidated a synthetic virionlike RNA species which contained just the terminal sequences of the virion RNA: the N encapsidation signal from the 5' end and the leader gene from the 3' end. This assembled nucleocapsid could function in vitro as a transcription template for the VSV RNA polymerase.

Polymerase errors accumulating during natural evolution of the glycoprotein gene of vesicular stomatitis virus Indiana serotype isolates

We report the entire glycoprotein (G) gene nucleotide sequences of 26 vesicular stomatitis virus Indiana serotype (VSV IND) type 1 isolates from North and Central America. These sequences are also compared with partial G gene sequences of VSV IND type 2 (Cocal) and type 3 (Alagoas) viruses and the complete G gene sequences of the more distantly related VSV New Jersey (NJ) and Chandipura viruses. Phylogenetic analysis of the G gene sequences by maximum parsimony revealed four major lineages or subtypes within the classical VSV IND (type 1) viruses, each with a distinct geographic distribution. A high degree of VSV genetic diversity was found in Central America, with several virus subtypes of both VSV IND and NJ serotypes existing in this mainly enzootic disease region. Nineteen percent sequence variation but no deletions or insertions were evident within the 5' noncoding and the coding regions of the VSV IND type 1 G genes. In addition to numerous base substitutions, the 3' noncoding regions of these viruses also contained numerous base insertions and deletions. This resulted in striking variation in G gene sizes, with gene lengths ranging from 1,652 to 1,868 nucleotides. As the VSV IND type 1 subtypes have diverged from the common ancestor with the NJ subtypes, their G mRNAs have accumulated more 3' noncoding sequence inserts, ranging up to 303 nucleotides in length. These primarily consist of an imprecise reiteration of the sequence UUUUUAA, apparently generated by a unique polymerase stuttering error. Analysis of the deduced amino acid sequence differences among VSV IND type 1 viruses revealed numerous substitutions within defined antigenic epitopes, suggesting that immune selection may play a role in the evolution of these viruses.

From rabies to rabies-related viruses

Antigenic differences between rabies virus strains characterized with monoclonal antibodies presently define at least four serotypes within the Lyssavirus genus of the Rhabdoviridae family: classical rabies virus strains (serotype 1), Lagos bat virus (serotype 2), Mokola virus (serotype 3) and Duvenhage virus (serotype 4). The wide distribution of rabies-related virus strains (serotypes 2, 3 and 4) and above all, the weak protection conferred by rabies vaccines against some of them (principally Mokola virus) necessitates the development of new specific vaccines. We first determined the complete nucleotide sequence of a rabies virus strain of serotype 1 (Pasteur virus) and characterized the structure of the viral genes and their regulatory sequences. We then extended this study to the Mokola virus genome. Five non-overlapping open reading frames were found in both viruses and had similar sizes and positions in both. Similarities were also found in the mRNA start and stop sequences and at the genomic extremities. Comparison of both genomes helps to analyze the basis of the particular antigenicity of these two serotypes. The sequence homology in the region coding for the viral glycoprotein was only 58% between the two viruses, compared with 94% between different rabies virus strains within serotype 1. This comparison, extended to other unsegmented negative strand RNA viruses, gives new insight into the understanding of rhabdoviruses and paramyxoviruses. Furthermore, molecular cloning provides a rationale for the genetic engineering of a future vaccine.

High nucleotide substitution error frequencies in clonal pools of vesicular stomatitis virus

Nucleotide substitution error frequencies were determined for several specific guanine base positions in the genomes of cloned vesicular stomatitis virus populations. Predetermined sites were examined in coding regions for the N, M, and L proteins and at a site in the genome 5'-end regulatory region. Misincorporation frequencies were estimated to be on the order of 10(-3) to 10(-4) at all positions analyzed. Isolates taken from virus populations after disruption of equilibrium conditions displayed replicase fidelity similar to that of cloned wild-type vesicular stomatitis virus. These mutation frequencies apply to all virus genomes present, including viruses rendered nonviable by lethal mutations. At one selected site in the N gene, two of three G----N base substitutions generated lethal nonsense mutations, yet their frequency was also very high. Biological implications for rapid virus evolution are discussed.

Completion of the rabies virus genome sequence determination: highly conserved domains among the L (polymerase) proteins of unsegmented negative-strand RNA viruses

We have now completed the rabies genome structure by the cloning and the sequencing of the entire L gene and the 5' untranscribed region. The L gene encodes a single open reading frame 2142 amino acids in length (244,206 Da) that corresponds to the viral RNA-dependent RNA polymerase. In contrast with other isofunctional proteins, the rabies polymerase exhibits a high degree of homology with the vesicular stomatitis virus polymerase, and a lesser degree, although significant, with those of Sendai virus and Newcastle disease virus, which suggests a differential evolution of the different cistrons. We have observed several strongly conserved stretches which may designate the independent functional domains of this multifunctional protein. In addition to the conservation of related transcription signals (N. Tordo et al. (1986) Proc. Natl. Acad. Sci. USA 83, 3914-3918.), this highlights the striking selective pressure on elements involved in transcription and replication mechanisms, and provides further evidence for a common ancestry of Rhabdoviridae and Paramyxoviridae families. The terminal complementarity observed in the rabies genome suggests the conservation of important genomic signals.

The 5'-terminal sequence of VSV(NJ) (Ogden): is the interaction of the NS protein with the NS binding site responsible for heterotypic interference activity?

The 5'-terminal sequence of VSV(NJ) (Ogden) and VSV(NJ) (Hazelhurst) was compared in an attempt to understand why the defective interfering particle, DI-LT, heterotypically interferes with VSV(NJ) (Ogden) but not with VSV(NJ) (Hazelhurst). The 5'-terminal sequence of VSV(NJ) (Ogden) genomic RNA was determined by direct RNA sequencing and by DNA sequencing of cDNA clones of the 3'-terminal sequence of VSV(NJ) (Ogden) DI particle genome. Primer extension analysis of the 5'-terminus of VSV(NJ) (Ogden) standard genomic RNA confirmed these data. Within the last 47 nucleotides, equivalent to the negative-strand leader RNA, the only nucleotide changes between VSV(NJ) (Ogden) and VSV(NJ) (Hazelhurst) occur between nucleotides 19 and 26, representing part of the putative NS binding region described by Isaac and Keene (J. Virol. 43, 241-249 (1982] for VSV(IND) DI particles. The spacer (S) region, located between the polyadenylation signal of the L gene and the 47th nucleotide of the leader RNA, contains more differences. The polyadenylation signal of the L gene is fully conserved, but the remainder of the L gene region (177 nucleotides) has highly diverged between VSV(NJ) (Ogden) and VSV(NJ) (Hazelhurst). The changes in the NS binding region of the negative-strand leader RNA provide further evidence for the divergent evolution of VSV(NJ) (Ogden) and VSV(NJ) (Hazelhurst). The NS binding region has been implicated as a crucial site for the initiation of RNA transcription and replication. The interaction of the NS protein with this site may determine the ability of DI particles to interfere heterotypically.