Open Access

Molecular characterization of China rabies virus vaccine strain

  • Wenqiang Jiao1,
  • Xiangping Yin1,
  • Zhiyong Li1,
  • Xi Lan1,
  • Xuerui Li1,
  • Xiaoting Tian1,
  • Baoyu Li1,
  • Bin Yang1,
  • Yun Zhang1 and
  • Jixing Liu1Email author
Virology Journal20118:521

https://doi.org/10.1186/1743-422X-8-521

Received: 20 September 2011

Accepted: 17 November 2011

Published: 17 November 2011

Abstract

Background

Rabies virus (RV), the agent of rabies, can cause a severe encephalomyelitis in several species of mammals, including humans. As a human rabies vaccine strain employed in China, the genetic knowledge of the aG strain has not been fully studied. The main goal of the present study is to amplify the whole genome of aG strain, and genetic relationships between other vaccine strains and wild strains were analyzed.

Results

The entire genome of human rabies virus vaccine strain aG employed in China was sequenced; this is the second rabies virus vaccine strain from China to be fully characterized. The overall organization and the length of the genome were similar to that of other lyssaviruses. The length of aG strain was 11925nt, comprising a leader sequence of 58nt, nucleoprotein (N) gene of 1353nt, phosphoprotein (P) gene of 894 nt, matrix protein (M) gene of 609nt, glycoprotein (G) gene of 1575nt, RNA-dependent RNA polymerase (RdRp,L) gene of 6384nt, and a trailer region of 70 nt. There was TGAAAAAAA (TGA7) consensus sequence in the end of each gene, except AGA7 at the end of G gene. There was AACAYYYCT consensus start signal at the beginning of each gene.

Conclusions

In this report, we analyzed the full genome of China human rabies vaccine strain aG. Our studies indicated that the genome of aG retained the basic characteristics of RV. At gene level, N was the most conserved among the five coding genes, indicating this gene is the most appropriate for quantitative genotype definition. The phylogenetic analysis of the N indicated the aG strain clustered most closely with Japanese and Russian rabies vaccine strains, suggesting that they may share the same ancestor; also, the aG strain did not share high homology with wild strains isolated from China, making it may not be the best vaccine strain, more research is needed to elucidate the genetic relationship among the RV circulating in China.

Keywords

Rabies virusaGChinafull-length genome

1. Background

Rabies is a widespread neurological zoometric disease which affects almost all kinds of mammals, including humans. The morality is almost 100%. It is reported that nearly 55,000 human die of rabies all over the world each year, most occurred in Asia and Africa [1]. Following India, China has the second highest number of human cases in the world [2]. The causal agent of the disease is rabies virus (RV), a member of Lyssavirus genus, Rhabdoviridae family. The RV genome is almost 12 kb in size consisting of five monocistronic RNAs, which encodes the nucleocapsid protein (N), phosphoprotein (P), matrix protein (M), glycoprotein (G) and RNA-dependent RNA polymerase (large protein). Between the monocistronic there are four intergenic regions (IGRS) with different length [3]. Besides, at the end of the genome there are two underinflated regions (UTRS), which play an important role in virus replication.

So far, vaccination is the most effective method to prevent rabies. The vaccine strains employed in China are CTN, aG, PM and PV. The complete genome sequence of CTN strain has been sequenced already [4]. Although novel vaccines, such as oral vaccine and live attenuated vaccine have been developed in recent years, the virus strain used for the production of vaccine are still pathogenic for laboratory and wildlife spices, what is more important, these rabies vaccines have the potential of causing vaccine-induced rabies and appear to have limited efficacy in a broad spectrum of animal species [5]. Therefore, genetic knowledge of vaccine strain is essential. Also, comparison of sequence of vaccine strain and field strain circulating in the country would prove how effective of the vaccine against the field virus [6].

The aG strain, also known as pG strain, was obtained from a rabid dog in 1931, Beijing. After 50 passages in the brain of rabbit, the strain was named aG strain. Then, the aG strain has been adapted to VERO cell cultures after 68 passages. Much higher antibody was obtained in immunized dogs, and the morality was lower. Since then, the aG strain has been chosen for the production of human rabies vaccine. Before this study, P, M, G and L genes of aG strain have been sequenced; however, this is the first time that the complete genome has been sequenced. In the present paper, the complete sequence of the rabies virus vaccine strain aG employed in China was sequenced. The sequence obtained indicates that aG strain genome is similar to the sequence of other complete lyssavirus genome available in Genbank. Sequence comparison with other rabies virus strain and phylogenetic analysis according to the N gene also falls within the scope of this study.

2. Materials and methods

2.1. Virus strain

The aG strain was prepared and stored by State Key Laboratory of Veterinary Etiological Biology, Lanzhou Veterinary Research Institite, Chinese Academy of Agricultural Science (LVRI, CAAS). The parental virus was isolated from a rabid dog's brain from Beijing, China, of a patient died of rabies from Beijing 1931, China. After 31 passages in primary hamster kidney (PHKC) cell, it was chosen for vaccine production in China [4].

2.2 Primer design

According to the conserved regions of the genome of rabies virus published in GenBank, particular the RC-HL(Genbank number: AB009663), RV-97(Genbank number EF542830), Nishigahara(Genbank number: AB044824) and NI-CE(Genbank number: Ab128149), 12 pairs of primers were designed to amplify the unknown regions of aG strain by primer premier 5.0 software(Table 1). All primers were synthesized by Shanghai Sangon Biological Engineering Technology & Service Co., Ltd.
Table 1

Primers employed for the amplication of unknown genomic RNA of rabies virus aG strain

Fragments

Sequence(5'-3')

Lengths(bases)

Base position

F1aFb

ACGCTTAACAACAAGATCA

19

1-1581

F1R

CTTCAGCCATCTCAAGAT

18

 

F2F

TATGTCTCAGTCAGTTCC

18

1340-2933

F2R

TCATCCCAAGTGATCTCC

18

 

F3F

ATACTCTGGGAATCATAGGAT

21

2742-4280

F3R

TGAGACGTCTGAAACTCACTG

21

 

F4F

AACATGGGTCGCGATGCA

18

4090-5453

F4R

GGTCATCATAGACCTCTC

18

 

F5F

CATGTAGATTCTCATAAT

18

5063-6410

F5R

TAAATACAGGAAAGTCTC

18

 

F6F

TGTGGAAACTCCGGCTAT

18

6290-7430

F6R

CGATGAGGTCTGATCTGTCTG

21

 

F7F

TGTTTGGATTGAAGAGAGTGT

21

7344-8513

F7R

GCTCACTGAGAAATCGAG

18

 

F8F

TTCAGAGTTTCGAGAGGCAAT

21

8410-9648

F8R

ATGTTGACAGGGAAGATGGTT

21

 

F9F

ATGTTCCAGCCATTGATGCTT

21

9272-10598

F9R

TGAACACAAGCTTGGCATCAG

21

 

F10F

TGCTCTGCTCAACAGGTT

18

10331-11932

F10R

ACGCTTAACAAATAAACAACA

21

 

F11F

GGTCTGGCGACACCCCGGTGTTCA

24

 

F11R

GAGTCTGTCATCTCACTGGATCA

23

 

F12F

TGCTGCTGCCAAGTAGGAACAT

22

 

F12R

AGTTTGGCGGCATCCATGCCTG

22

 

a The 12 fragments are labeled from "F1" to "F12"

b "F" and "R" strand for forward and reverse direction of the primer extension

2.3 Reverse transcription-PCR and direct sequencing

Total RNA were extracted from lyophilized brains of mice infected with aG strain by using an Rneasy Mini Plus Total RNA extraction kit (Qiagen) according to the manufacture's instruction. Single cDNA was synthesized by using First Strand cDNA Synthesis Kit (ToYOBO), and the cDNA were amplified using KOD Plus Neo (TOYOBO). Amplified DNA products were examined via electrophoresis using 0.8% ultrapure TM Agarose (Invitrogen) gels. DNA products of expected size were purified using Agarose Gel DNA Purification Kit (TAKARA). The purified products were added to A at the end of 3'end using DNA A-tailing Kit (TAKARA) and then ligated with pMD18-T (TAKARA), the ligated products were transformed into competent cell JM109 (TianGen), at least 3 clones carrying the PCR products of expected size were identified by PCR and the positive clones were sequenced commercially (Shanghai Sangon biological Engineering Technology &Services Co., Ltd., China).

2.4 Amplification and sequencing of the terminal ends

The 5'and 3'terminal of the genome was confirmed using the approach described by Kuzmin et al [7]. Briefly, Total RNA was subjected to ligation by T4 RNA ligase (Invitrogen) to circularize the genome. Thereafter the ligated RNA was subjected to nested PCR with the primer marked with F 11F, 11R, 12F and 12 R located within the 5'and 3'end of the genome. Products of the expected size were dealt with above.

2.4 Sequence alignment and analysis

Nucleotide sequences obtained were edited and assembled manually using EditSeq of DNAStar program 7.0 (DNAStar). Phylogenetic analysis of N gene of Lyssavirus, including the construction of phylogenetic tree were conducted using Molecular Evolutionary Genetics Analysis (MEGA) Version 5.0 with the methods of neighbor-joining (NJ) algorithm with the Kimura two-parameter model. The reference sequences have been listed in Table 2. The reliability of the phylogeny groupings was evaluated using bootstrapping with 1000 replicates. Bootstrap values of 70% or greater were viewed significant.
Table 2

Complete genome of rabies virus used in this study from Genbank

Strain

Origins

Country

Accession no

Length

CTN181

Attenuated vaccine strain, derived from the brain tissue of a patient in Shandong, China

China

EF564174

11923 bp

Nishigahara

Seed strain of animal rabies vaccine strain used in Japan

Japan

AB044824

11926 bp

SAD-B19

A highly attenuated mutant of SAD strain

America

AB044824

11928 bp

RC-HL

Animal vaccine strain used in Japan

Japan

AB009663

11926 bp

RV-97

Vaccine strain employed in Russia

Russia

EF542830

11932 bp

PV

Pasteur vaccine strain

France

M13215

11932 bp

ERA

Attenuated rabies vaccine strain

America

EF206707

11931 bp

SRV9

Avirulent vaccine strain

China

AF499686

11928 bp

NI-CE

Developed from Nishigahara

Japan

AB128149

11926 bp

HEP-Flury

Derived from a human isolate in Georgia, USA, 1939

America

AB085828

11615 bp

DRV

Isolated from a deer in Ji lin, China

China

DQ875051.

11863 bp

MRV

Isolated from a mouse in China

China

DQ875050.

11869 bp

HN10

Isolated fron a rabies patient in Hunan, China

China

EU643590

11923 bp

SHBRV-18

Isolated from a silver-haired bat

America

AY705272

11923 bp

BD06

Isolated in China

China

EU643590

11924 bp

NNV-RAB-H

Isolated from the brain of human

India

EF437215

11928 bp

SAG2

Vaccine strain Derived from SAD

America

EF206719

11928 bp

8743THA

A street strain

Thailand

EU293121

11923 bp

8764THA

A street strain

Thailand

EU293111

11925 bp

H-08-1320

Isolated from the brain of human

Sri Lanka

AB569299

11926 bp

RRV ON-99-2

Isolated from a raccoon

Canada

EU311738

11923 bp

BR-Pfx1

Isolated from Dusicyon sp

Brazil

AB362483

11924 bp

9704ARG

Isolated from tadarida brasiliensis

Argentina

EU293116

11923 bp

9147FRA

Isolated from fox

France

EU293115

11923 bp

9001FRA

Isolated from dog

Guyana

EU293113

11922 bp

FJ08

Isolated from dog

China

FJ866835

11924 bp

FJ09

Isolated from dog

China

FJ866836

11924 bp

D01

Isolated from the brain of dog

China

FJ12193

11925 bp

D02

Isolated form the brain of dog

China

FJ12194

11925 bp

MOKOLA

Isolated from Mokola Virus

Africa

NC_006429

11940 bp

KRC5-04

Isolated from dog

South Korea

AY730597

1353 bp

HNDB28

Isolated from dog

China

EU008922

1353 bp

HNDB12

Isolated from dog

China

EU008920

1353 bp

GXPL

Isolated from Canis lupus familiaris

China

GQ472474

1353 bp

Hebei0(H)

Isolated from homo sapiens

China

EU267777

1353 bp

Zhejiang Wz1(H)

Isolated from homo sapiens

China

EU700032

1353 bp

West Caucasian bat virus

Isolated from West Caucasian bat

Russia

EF614258

12278 bp

Flury-LEP

Derived from a human isolate in Georgia, USA, 1939

America

DQ099524

11711 bp

JS08-45

Isolated from Chinese ferret badger

China

GU647092

11922 bp

3. Results

3.1 Genome organization of the aG strain

Using a total of 24 primers (as shown in table 1) the strain aG genome was obtained as 12 separate overlapping PCR products. The 5'- and 3'-terminal were obtained as mentioned above. The length of the genome is 11925 nucleotides(nt), and the genome organization of aG strain, which follows typical rabies virus organization, is summarized as follows: a 3' leader region of 58 nt(1-58), the N gene (59-1482), P gene (1485-2476), M gene (2482-3284), G gene (3290-5356), L gene (5381-11855), and the 5'trailer region of 70 nt (11855-11925). The coding sequence (CDS) of each of the five structural genes are as follows:1353-nt nucleotide (71-1423), 894-nt phosphoprotein(1515-2408), 609-nt matrix protein (2497-3105), 1575-nt glycoprotein (3317-4891) and 6384-nt large protein(5411-11794). The G-L non-coding region is 455nt. The sequences of the transcription initiation sites (TIS) and the termination signals of the structural genes and the intergenic regions (IGRS) among these genes are shown in table 3.
Table 3

Transcription initiation and termination signals for all strain aG

Gene or region

Termination signal

IGRs

Initiation signal

3'leader/N

------

-------

AACACCTCT

N/P

TGAAAAAAA

CT

AACACCCCT

P/M

TGAAAAAAA

CGAGC

AACACCACT

M/G

TGAAAAAAA

CTATT

AACATCCCT

G/L

AGAAAAAAA

CTGTAGATCGAAAGAGCAACTGGC

AACACCTCT

L/5'trail

CGAAAAAAA

-------------------------------------------------

-------

3.1.1. Nucleotide feathers of strain aG

The length and sequence of leader region is highly conserved in all lyssavirus. The length of the leader region is 58nt. The first 12 nucleotides, especially the first 9 are the same in all the lyssavirus. However, the conservation is not so stringent after residue 25. Following most lyssavirus, the trailer regions of aG strain is 70nts.

The complementarity of the 3'and 5' ends of the genome is another typical feather of the Mononegavirals, with up to 16 nts are complemented. The first 9 nts of the genome leader regions are absolutely complementary to the corresponding regions of the trailer regions [7] (Figure 1).
Figure 1

Comparison of the 3' and 5' genomic termini of the antigenomic (+) sense RNA of aG strain.

The transcription initiation signals (TISs) and transcription termination signals (TTS) are conserved as AACAYYHCT and G (A) 7, but the TTS between G-L was AGAAAAAAA at the N-P junction. The intergenic regions between the cistrons CT, CCGAA, ACTATT and CTGTAGATCGAAAGAGCAACTGGC (Table 3). Pairwise comparisons of nucleotide and predicted amino acid sequence identities of the aG with the vaccine strain of lyssaviruses are listed in Table 4.
Table 4

Pairwise comparisons of nucleotide and predicted amino acid sequence identities of the aG with the vaccine strain of lyssaviruses

Strain

 

N (%)

P (%)

M (%)

G (%)

L (%)

SRV9

Nucleotide

92.5

89.6

91.1

91

92.2

 

Amino acid

96.2

88.9

87.7

89

95.3

CTN181

Nucleotide

87.3

82.9

83.9

84.1

25.1

 

Amino acid

96

87.2

93.1

88.8

5.9

ERA

Nucleotide

92.6

89.7

92.3

91.4

92.4

 

Amino acid

96.2

89.3

90.6

89.9

95.9

Flury-LEP

Nucleotide

92.1

88.6

92.6

90.5

91

 

Amino acid

96.9

87.9

93.6

91.2

96.5

Flury-HEP

Nucleotide

92.1

87.9

91.5

89.8

91.5

 

Amino acid

96.9

87.2

92.1

89.9

96.4

NI-CE

Nucleotide

97.9

95.6

96.9

97.1

98.6

 

Amino acid

98.4

93.3

96.1

96.6

98.6

Nishigahara

Nucleotide

98.1

96.8

97.2

97.1

98.2

 

Amino acid

99.1

95

97

96.4

98.7

PV

Nucleotide

92.3

89.6

92

91.6

92.3

 

Amino acid

95.8

89.9

90.6

90.5

95.5

RV97

Nucleotide

95.9

92.5

92.3

93

94.9

 

Amino acid

96

87.9

90.6

89.7

96.1

SAD-B19

Nucleotide

92.5

89.6

91.1

91.2

92.2

 

Amino acid

96.2

88.9

87.7

89.7

95.3

SHBRV-18

Nucleotide

84.2

79.6

80.8

80

82

 

Amino acid

94.7

83.9

88.7

84.2

5.8

3.2 Structural feathers of the proteins of aG strain

Consistent with previous studies, at the nucleotide level N was the most conserved gene [8]. The N protein of aG strain is 450 amino acid residues. Only 10Ala amino acid substitution is unique to the strain aG. Antigenic site I (residues 358-367), antigenic site IV (375-383, 359-366) and RNA-binding domain (residues 298-352) were found to be conserved in all the isolates analyzed. Ser389, which was considered to be related to casein-type phosphorylation site and regulation of viral RNA transcription and replication, was highly conserved [9] (Figure 2).
Figure 2

Alignment of Antigenic sites I and IV, RNA-binding site of N. Dots represent identity to the majority, shade (with solid black) residues differ from consensus. A. Alignment of Antigenic sites I and IV. Dots represent identity to the majority. B. Alignment of RNA-binding site of N, shade (with solid black) residues differ from consensus.

P gene was considered to be the most divergent protein among the five structural proteins. Six amino acid substitutions, Lys68, Tyr95, Thr174, Thr184, Glu254 and Glu295, unique to the aG were observed. The interaction of the motif [(K/R)XTQT] between residues 145-149 with cytoplasmic dynein light chain (LC8), which plays an important role in viral nucleocapsid axoplasmic transport, was encoded as KSTQT [10]. Compared with vaccine strain, the N protein binding site (69-177 and 268-297) of aG strain had one unique amino acid residues substitution, Tyr95. The L binding site in P, the first 19 amino acid residues, was absolutely conserved in aG. Ser162, Ser210, Ser271, which have been shown to be involved in the phosphorylation of P were well conserved in the aG strain [11] (Figure 3).
Figure 3

Alignment of N binding sites (69-177), (268-297)and [(K/R) XTQT], dots represent identity to consensus. A. Alignment of N binding sites (69-177) and [(K/R) XTQT], dots represent identity to consensus. B. Alignment of N binding site (268-297), dots represent identity to consensus.

M protein connects the plasma membranes, RNPS and G protein together, and is deeply involved in the budding of rabies virus. Thr100, Ser111 and Ser174 were the three amino acid substitution to the strain aG. The motif proline-rich (PPxY) between residues 35-38 is thought to be associated with the interaction with WW domains of cellular components, and was encoded PPEY, which is very conserved [12]. Residue 58, which is crucial for regulation of RV RNA synthesis [13] (Figure 4), was encoded Glu in aG strain and Gly in PV strain.
Figure 4

Alignment of proline-rich motif and residue 58. Dots represent identity to the consensus.

The G protein of 514 amino aicd residues is the only viral protein exposed on the surface of the virus and responsible for the immune responses of the host. Compared with vaccine strain analysed in this study, there are eight amino acid substitutions unique to aG strain: His69, Pro184, Pro250, Gly427, Ile431, Ile477, Lys481, Asn160. Antigenic site I (residue 231), antigenic site II(residues 34-42,198-200), antigenic site III (residues 330-338), antigenic site IV(residue 264) and antigenic site a (residue 342) were conserved. Ala242, sp255, Ile268 and Arg333, which are considered to be associated with viral pathogenicity, were replaced with Cys, Met, Thr and Ala respectly (Figure 5).
Figure 5

Alignment of antigenic antigenic siteI (residue 231), site II (residues 34-42), (residues 198-200), shade (with solid black) residues differ from the consensus. A. Alignment of antigenic site II (residues 34-42), shade (with solid black) residues differ from the consensus. B. Alignment of antigenic site I (residue 231) and II (residues 198-200), shade (with solid black) residues differ from consensus.

L protein is an important component of RNP, although the N and P protein are also needed to form RNP. The 2,127 amino acid residues of L protein are distributed into six conserved motifs.

Three potential N-glycosylation sites, which are located at 37-39, 158-160 and 319-321 were conserved in aG. The L gene of aG encoded a protein of 2,127-amino acid, and thirteen amino acid substitutions unique to strain aG were observed: Gly48, Gly58, Gln189, Tyr353, Gly489, Gly496, Ser606, Lys995, Val1043, Ser1133, Tyr1588, Leu1658, His1801. The motifs among 544-563(A), 728-732(B) and 1705-1710(C) have been regarded as functional motifs. These motifs are involved in RNA binding, active site of polymerase, polyadenylation or a core of ATP binding site. All these amino acid residues were conserved in strain aG (Figure 6). In the absence of clear functions for these residues, the importance of any of these changes is presently unknown.
Figure 6

Alignment of motif A, B, C. shade (with solid black) residues differ from consensus. A Alignment of motif A. shade (with solid black) residues differ from consensus. B Alignment motif B. shade (with solid black) residues differ from consensus. C. Alignment of motif C. shade (with solid black) residues differ from consensus.

3.2 Phylogenetic analysis of N

As shown by the phylogenetic analysis of the nucleoprotein gene (Figure 7), the three vaccine strain Nishigahara, NI-CE and RC-HL were grouped with aG most closely.
Figure 7

Phylogenetic analysis of N gene sequence of aG with nucleoprotein sequence available in Genbank by Neighbor-Joining program in MEGA 5 software. The number below the brace is bootstrap values for 100 replicates.

4. Discussion

In the present study, the entire aG genomic sequence, which is composed of 11925 bp, was determined by 12 paris overlapping fragments covering the whole genome to improve our knowledge of the genetic relationship between aG strain and other strain of rabies virus available in Genbank. In contrast to the previous report [14], the length of the whole genomes in GT 1 lyssaviruses varies from 11711-nt to 12278-nt, with both even and odd numbers of nucleotides being represented [15, 8]. In contrast to previous studies, the laboratory adapted vaccine strain aG has a odd nucleotide, while other laboratory isolates, such as Nishigahara, PV and SAD-B19 strain, have genomes of even nucleotides numbers [16]. It seems that genomes of odd numbers of nucleotides are more common in wild type strains than fixed strain but both odd and even numbers of genomes occurred in wild type strains and fixed strain. It remains to be determined that whether the isolates with odd or even numbers of the genome has any biological significance. However, the number of the sequence does not follow the 'rule of six' [17].

On average, the G+C content of the lyssaviruses genome is 44.57%, the G+C content was 45.18 mol% for aG strain. This is in agreement with the idea that the G+C content of the negative-stranded RNA viruses is lower than that of positive stranded RNA viruses, which has been considered to be associated with host cell RNA editing [18]. In agreement with previous studies, the genome of aG strain is similar to that of other lyssaviruses with few variations. The complete nucleotide length of aG strain is 11,925nt, which was one, seven and three nucleotides shorter than that of Nishigahara, PV and SAD-B19 strain. Transcriptional initiation and termination signals of aG strain at the beginning and end of the each monocistron, including the motif AACAYYHCT initiatite transcription and a nine nucleotide motif WG(A)7, except the AGAAAAAAA TTP at the N-P junction, which is involved in transcription termination and polyadenylation (TTP), are precisely the same as described by Marston. All the intergenic sequences (IGS) in the aG agree precisely with what was described in that report. Unlike PV and SAD strains, the G-L intergenic region of aG strain has only one polyadenylation site. The 3' and 5' of the genome not only play an important role in transcription, replication and the switching between these two functions, but also in the initiation of encapsidation. The 11 nts at the leader and trailer region of the genome were complementary conserved. These regions showed that the basic structure of rabies virus have been retained in strain aG [19].

It is reported that the 3'NCR of G mRNA to be the pseudogene, because there were two TTP motifs in PV. However, in aG strain, there is only one TTP motif, which holds true for many other isolates, including laboratory fixed strains and street rabies virus strains. Thus, a conclusion can be drawn that the psdudogene is only part of non-translated region of G mRNA.

It is assumed that the aG strain in China is more closely with rabies virus in northern and northeast part of China, since these strains are closely related geographically and ethnically. The phylogenetic tree based on the N gene showed that aG strain clustered with the Japanese vaccine strain, supporting the idea that rabies virus from the same region tend to cluster together.

Consistent with previous observations that viruses from the same geographical area tend to group together, the aG strain tends to cluster most closely with DRV, which was isolated from deer in Jilin province, China. Also, the aG strain grouped closely to other vaccine strain, especially the vaccine strain from Japan [20]. However, the aG strain did not share high nucleotide homology with wild strain, such as FJ008, FJ009 and HN10, making aG not the best vaccine in China [21], although the strain aG has been employed for human vaccine strain in China for a long time. Of course, more research is needed to confirm the hypothesis.

To date, several methods, such as rapid amplification of cDNA (RACE) [2224] and gene-walking approach have been employed to obtain the whole genome sequence, including the 5' and 3' untranslated regions (UTR). In the present study, a simple method of amplification and sequencing of full length of rabies virus genome was performed, according to reference 7. Before this study, similar method has been used to sequence the genome of different viruses [25]. Compared with the method previously, the method described here is easy to perform, does not need complex principle and does not need expensive device; thus it is recommended for other virus.

Comparison of nucleotide sequences of rabies virus available in Genbank with aG strain showed extensive divergence excluding the transcription regulatory signals and limited stretches of the nucleotide. However, according to the phylogenetic analysis of the N gene of aG strain, the strain probably emerged from the DRV, or they may share the same ancestor, which was isolated deer in Jilin province, China. More researches are needed to determine whether the aG strain came originally as a result of the DRV adapting to the dog host or not. The Ser389, the putative casein-type phosphorylation site and regulation of viral RNA transcription and replication was conserved in strain aG. The antigenic sites in strain aG were conserved, indicating the antigenicity of aG was well retained, which was suitable for vaccine production.

P has been considered to be the most divergent of the five coding protein. The LC8-interacting motif was conserved in aG strain, indicating that the motif serve to transport the virus through neurons. Also, the L-binding region of P, as well as the N-binding site of P was also conserved, except Thr174, Glu295. These results, together with the five Ser and four Met, indicating that the same function which have been reported in fixed strains were also retained in aG strain. Mebatsion and Rasalingam questioned the function of the interaction between P and LC8 by constructing of recombinant virus in the absence of the LC8 binding site. The results indicated that the virus lacking the binding site was as pathogenic as wild virus, suggesting the interaction is not indispensible to the spread of RV.

M protein has been considered to be important in membrane budding and interacting with cellular components. Finke found that substitution of Arg58 to Gly58 reduced the expression of full-length RNA. The aG strain retains the Glu58, suggesting that it may serve to RNA regulation.

The G protein is related to cell attachment, responsible for the induction of neutralizing antibodies and cell-mediated immune response [26]. Ala242, Asp255, Ile268, and Arg333, which are correlated with pathogenicity of rabies virus, were encoded Cys, Met, Thr and Ala respectively. The antigenic sites were all retained in aG. These results suggest the antigenicity of G protein was retained. However, the strain aG was strongly attenuated, although other factors may be responsible for the pathogenicity of rabies virus.

The L protein, together with P compose RNA polymerase complex, which is involved in enzymic activities of polymerase. The RNA-binding region, active site of polymerase, the Pro residues rich motif GXGXG, which is involved in polyadenylation or protein kinase activities, were conserved; These motifs, together with GHP were conserved in strain aG. All These results suggest the L protein possess the same function as other strain. However, little is known about the mechanism of the function of L protein, therefore, more researches are needed and necessary.

Conclusion

In this report, we analyzed the full genome of China human rabies vaccine strain aG. Our studies indicated that the genome of aG retained the basic characteristics of RV, including the organization of the genome, the functional domains of the each ORF. At gene level, N was the most conserved among the five coding genes, indicating We this gene is the most appropriate for quantitative genotype definition. The phylogenetic analysis of the N indicated the aG strain clustered most closely with Japanese and Russian rabies vaccine strains, suggesting that they may share the same ancestor; also, the aG strain did not share high homology with wild strains isolated from China, making it may not be the best vaccine strain, more research is needed to elucidate the genetic relationship among the RV circulating in China.

List of abbreviations

RACE: 

rapid amplification of cDNA ends

RV: 

rabies virus

RT-PCR: 

Reverse transcription polymerase chain reaction

RNP: 

ribonucleoprotein

PHKC: 

primary hamster kidney cell.

Declarations

Acknowledgements

This work was financially supported by the Lanzhou biological medicine foundation (2011-1-61).

Author Details

State Key Laboratory of Veterinary Etiological biology, Key Laboratory of Epizootic disease of Grazing Animal of Ministry of Agriculture, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science (LVRI, CAAS), Lan Zhou, 730046, China.

Authors’ Affiliations

(1)
State Key Laboratory of Veterinary Etiological biology, Key Laboratory of Epizootic disease of Grazing Animal of Ministry of Agriculture, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science(LVRI, CAAS), Xujia ping1, Yanchang bu

References

  1. Knobel DL, Cleaveland S, Coleman PG, Fevre EM, Meltzer Ml, Miranda Me, Shaw A, Zinsstag J, Meslin FX: Re-evaluating the burden od rabies virus in Africa and Asia. Bull World Health Organization. 2005, 83 (5): 360-368.Google Scholar
  2. Tang Q, Li H: Epidemic situation and related factors analysis of rabies in China. Chin J Epidemiol. 2005, 26: 223-224.Google Scholar
  3. Mita T, Shimizu K, Ito N, Yamada K, Ito Y, Sugiyama M, Minamoto N: Amino acid at position 95 of the matrix protein is a cytopathic determinant of rabies virus. Virus Res. 2008, 137 (1): 33-39. 10.1016/j.virusres.2008.05.011.View ArticlePubMedGoogle Scholar
  4. Du J, Zhang Q, Tang Q, Li H, Tao X, Morimoto K, Nadin-Davis SA, Liang G: Characterization of human rabies virus vaccine strain in China. Virus Res. 2008, 135 (2): 260-266. 10.1016/j.virusres.2008.04.002.View ArticlePubMedGoogle Scholar
  5. Blanton JD, Self J, Niezgoda M, Faber ML, Dietzschold B, Rupprecht C: Oral vaccination of raccoons (Procyon lotor) with genetically modified rabies virus vaccines. Vaccine. 2007, 25 (42): 7296-7300. 10.1016/j.vaccine.2007.08.004.PubMed CentralView ArticlePubMedGoogle Scholar
  6. Metlin A, Paulin L, Suomalainen S, Neuvonen E, Rybakov S, Mikhalishin V, Huovilainen A: Characterization of Russian rabies virus vaccine strain RV-97. Virus Res. 2008, 132 (1-2): 242-247. 10.1016/j.virusres.2007.11.016.View ArticlePubMedGoogle Scholar
  7. Kuzmin IV, Wu X, Tordo N, Rupprecht CE: Complete genomes of Aravan, Khujand, Irkut and West Caucasian bat viruses, with special attention to the polymerase genen and non-coding regions. Virus Res. 2008, 136 (1-2): 81-90. 10.1016/j.virusres.2008.04.021.View ArticlePubMedGoogle Scholar
  8. Szanto AG, Nadin-Davis SA, White BN: Complete genome sequence of a raccoon rabies virus isolate. Virus Res. 2008, 136 (1-2): 130-9. 10.1016/j.virusres.2008.04.029.View ArticlePubMedGoogle Scholar
  9. Yang J, Koprowski H, Dietzschold B, Fu ZF: Phosphorylation of rabies virus nucleoprotein regulates viral RNA transcription and replication by modulating leader RNA encapsidation. J Virol. 1999, 73 (2): 1661-1664.PubMed CentralPubMedGoogle Scholar
  10. Lo Kw, Naisbitt S, Fan JS, Sheng M, Zhang M: The 8-Kda dynein light chain binds to its targets via a conserved (K/R) XTQT motif. J Biol Chem. 2001, 276 (17): 14059-14066.PubMedGoogle Scholar
  11. Gupta AK, Blondel D, Choudhary S, Banerjee AK: The phosphoprotein of rabies virus is phosphorylated by a unique cellular protein kinase and specific isomers of protein kinase C. J Virol. 2000, 74 (1): 91-98. 10.1128/JVI.74.1.91-98.2000.PubMed CentralView ArticlePubMedGoogle Scholar
  12. Jayakar HR, Murti KG, Whitt MA: Mutations in the PPPY motif of vescular stomatitis virus matrix protein reduce virus budding by inhibiting a late step in virion replace. J Virol. 2000, 74 (21): 9818-9827. 10.1128/JVI.74.21.9818-9827.2000.PubMed CentralView ArticlePubMedGoogle Scholar
  13. Finke S, Conzelmann KK: Dissociation of rabies virus matrix protein functions in regulation of viral RNA synthesis and virus assembly. J Virol. 2003, 77 (22): 12074-12082. 10.1128/JVI.77.22.12074-12082.2003.PubMed CentralView ArticlePubMedGoogle Scholar
  14. Delmas O, Holmes EC, Talbi C, Larrou F, Dacheux L, Bouchier C, Bourhy H: Genomic diversity and evolution of the Lyssavirus. Plos ONE. 2008, 3 (4): e2057-10.1371/journal.pone.0002057.PubMed CentralView ArticlePubMedGoogle Scholar
  15. Warrilow D, Smith IL, Harrower B, Smith GA: Sequence analysis of an isolate from a fatal human infection of Australian bat Lyssavirus. Virology. 2002, 297 (1): 109-119. 10.1006/viro.2002.1417.View ArticlePubMedGoogle Scholar
  16. Warrilow D, Smith IL, Harrower B, Smith GA: Sequence analysis of an isolate from a fatal human infection of Australian bat lyssavirus. Virology. 2002, 297 (1): 109-119. 10.1006/viro.2002.1417.View ArticlePubMedGoogle Scholar
  17. Calain P, Roux L: The rule of six, a basic feather for efficient replication of Sendai Virus Defective Interfering RNA. J Virol. 1993, 67 (8): 4822-4830.PubMed CentralPubMedGoogle Scholar
  18. Auewarakul P: Composition bias and genome polarity of RNA viruses. Virus Res. 2005, 109 (1): 33-37. 10.1016/j.virusres.2004.10.004.View ArticlePubMedGoogle Scholar
  19. Marston DA, Mcelhinney LM, Johnson N, Muller T, Conzelmann KK, Tordo N, Fooks AR: Comparative analysis of the full genome sequence of European bat lyssavirus type 1 and type 2 with other lyssaviruses and evidence for a conserved transcription termination and polyadenylation motif in the G-L 3'non-translated region. J Gen Virol. 2007, 88 (pt4): 1302-1314.View ArticlePubMedGoogle Scholar
  20. Metlin AE, Rybakov S, Gruzdev K, Neuvonen E, Huovilainen A: Genetic heterogeneity of Russian, Estonian and Finnish field rabies viruses. Arch virol. 2007, 152 (9): 1645-54. 10.1007/s00705-007-1001-6.View ArticlePubMedGoogle Scholar
  21. Faber M, Pulmanausahakul R, Nagao K, Prosniak M, Rice AB, Koprowski H, Schnell MJ, Dietzschold B: Identification of viral genomic elements responsible for rabies virus neuroinvasiveness. Proc Natl Acad Sci USA. 2004, 101 (46): 16328-32. 10.1073/pnas.0407289101.PubMed CentralView ArticlePubMedGoogle Scholar
  22. Geue L, Schares S, Schnick C, Kliemt J, Bekert A, Freuling C, Hoffmann B, Zanoni R, Martston D, Mcelhinney L, Johnson N, Fooks AR, Tordo N, Muller T: Genetic characterization of attenuated SAD rabies virus strains used for oral vaccination of wildlife. Vaccine. 2008, 26 (26): 3227-35. 10.1016/j.vaccine.2008.04.007.View ArticlePubMedGoogle Scholar
  23. Nagaraja T, Madhusudana S, Desai A: Molecular characterization of the full-length genome of a rabies virus isolate from India. Virus Genes. 2008, 36 (3): 449-59. 10.1007/s11262-008-0223-3.View ArticlePubMedGoogle Scholar
  24. Matsumoto T, Ahmed K, Wimalaratne O, Yamada K, Nanayakkara S, Perera D, Karunanayake D, Nishizono A: Whole-genome analysis of a human rabies virus from Sri Lanka. Arch Virol. 2011, 156 (4): 659-69. 10.1007/s00705-010-0905-8.View ArticlePubMedGoogle Scholar
  25. Morzunov SP, Winton JR, Nichol ST: The complete genome structure and phylogenetic relationship of infectious hematopoietic necrosis virus. Virus Res. 1995, 38 (2-3): 175-192. 10.1016/0168-1702(95)00056-V.View ArticlePubMedGoogle Scholar
  26. Takayama-Ito M, Inoue K, Shoji Y, Inoue S, Iijima T, Sakai T, Kurane I, Morimoto K: A highly attenuated rabies virus HEP-Flury strain reverts to virulent by single amino acid substitution to arginine at position 333 in glycoprotein. Virus Res. 2006, 119 (2): 208-15. 10.1016/j.virusres.2006.01.014.View ArticlePubMedGoogle Scholar

Copyright

© Jiao et al; licensee BioMed Central Ltd. 2011

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.