Trailer Region Research Articles

Construction and Rescue of a Minigenome System for Lentogenic D58 Strain of Avian orthoavulavirus 1 (Newcastle Disease Virus)

Background: In reverse genetics studies, minigenomes of negative sense RNA viruses are used to model the process of virus replication and evaluate the virus rescue system in the course of construction of recombinant virus. In the current study, a minigenome (MG) rescue system for lentogenic D58 vaccine strain of Avian orthoavulavirus 1 (AOAV-1) (Newcastle disease virus) was developed to establish a reverse genetic system for this strain. Methods: The minigenome was constructed by assembling the trailer and leader regions of AOAV-1 D58 strain flanking the reporter gene, Aequorea coerulescens green fluorescence protein (AcGFP) in a modified transcription vector. Construction of this Minigenome was done using a simple restriction free technique-Splicing by Overlap Extension PCR. The constructed minigenome was evaluated for the reporter gene expression by transfection in T7 RNA polymerase expressing BSR/T7 cell line. Further, the replication and packaging ability of the constructed minigenome were also analyzed. Result: The reporter gene expression was confirmed in BSR/T7 cells which simulates the virus replication process and verifies the virus rescue system. In addition, this proves the applicability of restriction free cloning by Splicing by Overlap Extension (SOE) PCR for precise construction of the minigenome. The developed minigenome system can be used as a potential tool in the reverse genetic rescue of AOAV-1 D58 strain for the development of marker vaccine.

Near-Complete Genome Sequence of an Avian Orthoavulavirus Type 13 Strain Isolated in South Korea in 2020.

ABSTRACTWe report the near-complete genome sequence of an avian orthoavulavirus 13 (AOAV-13) strain isolated from a wild goose fecal sample collected in South Korea in early 2020. The AOAV-13 sequence had a unique 3′ trailer region, including an 84-nucleotide (nt) deletion and a 24-nt insertion, compared to the most closely related Chinese genome sequence from 2015.

Novel avian orthoavulavirus 13 in wild migratory waterfowl: biological and genetic considerations.

Avian orthoavulavirus 13 (AOAV-13), formerly known as Avian paramyxovirus 13 (APMV-13), is found scatteredly in wild birds around the world. Although four complete genome sequences of AOAV-13 had been identified since the first discovery in Japan in 2003, the information available on the genetic variation and biological characteristics of AOAV-13 is still limited. In the present study, we isolated six AOAV-13 strains from fecal samples of wild migratory waterfowls during annual (2014-2018) viral surveillance of wild bird populations from wetland and domestic poultry of live bird markets (LBMs) in China. The phylogenetic analyses based on the HN and F genes showed that they had very close relationship and the molecular clock estimations showed a low evolutionary rate of AOAV-13. However, Bean goose/Hubei/V97-1/2015 is 1953nt in size (ORF, 1, 776nt), which is a unique size and longer than other reported AOAV-13 strains. Additionally, four repeats of conserved sequences "AAAAAT" was presented in the 5'-end trailer region of Swan goose/Hubei/VI49-1/2016, which is unprecedented in the AOAV-13. These findings highlight the importance of continuous monitoring the specific species of APMVs.

Joá yellow blotch-associated virus, a new alphanucleorhabdovirus from a wild solanaceous plant in Brazil.

We identified a novel plant rhabdovirus infecting native joá (Solanum aculeatissimum) plants in Brazil. Infected plants showed yellow blotches on the leaves, and typical enveloped bacilliform rhabdovirus particles associated with the nucleus were seen in thin sections by electron microscopy. The virus could be graft-transmitted to healthy joá and tomato plants but was not mechanically transmissible. RT-PCR using degenerate plant rhabdovirus L gene primers yielded an amplicon from extracted total RNA, the sequence of which was similar to those of alphanucleorhabdoviruses. Based on close sequence matches, especially with the type member potato yellow dwarf virus (PYDV), we adopted a degenerate-primer-walking strategy towards both genome ends. The complete genome of joá yellow blotch-associated virus (JYBaV) is comprised of 12,965 nucleotides, is less than 75% identical to that of its closest relative PYDV, and clusters with PYDV and other alphanucleorhabdoviruses in L protein phylogenetic trees, suggesting that it should be taxonomically classified in a new species in the genus Alphanucleorhabdovirus, family Rhabdoviridae. The genome organization of JYBaV is typical of the 'PYDV-like' subgroup of alphanucleorhabdoviruses, with seven genes (N-X-P-Y-M-G-L) separated by conserved intergenic regions and flanked by partly complementary 3' leader and 5' trailer regions.

In vitro trackable assembly of RNA-specific nucleocapsids of the respiratory syncytial virus

In vitro trackable assembly of RNA-specific nucleocapsids of the respiratory syncytial virus

The Emergence of Avian Orthoavulavirus 13 in Wild Migratory Waterfowl in China Revealed the Existence of Diversified Trailer Region Sequences and HN Gene Lengths within this Serotype.

Avian orthoavulavirus 13 (AOAV-13), also named avian paramyxovirus 13 (APMV-13), has been found sporadically in wild birds around the world ever since the discovery of AOAV-13 (AOAV-13/wild goose/Shimane/67/2000) in a wild goose from Japan in 2000. However, there are no reports of AOAV-13 in China. In the present study, a novel AOAV-13 virus (AOAV-13/wild goose/China/Hubei/V93-1/2015), isolated from a wild migratory waterfowl in a wetland of Hubei province of China, during active surveillance from 2013 to 2018, was biologically and genetically characterized. Phylogenetic analyses demonstrated a very close genetic relationship among all AOAV-13 strains, as revealed by very few genetic variations. Moreover, pathogenicity tests indicated that the V93-1 strain is a low virulent virus for chickens. However, the genome of the V93-1 virus was found to be 16,158 nucleotides (nt) in length, which is 12 nt or 162 nt longer than the other AOAV-13 strains that have been reported to date. The length difference of 12 nt in strain V93-1 is due to the existence of three repeats of the conserved sequence, “AAAAAT”, in the 5′-end trailer of the genome. Moreover, the HN gene of the V93-1 virus is 2070 nt in size, encoding 610 aa, which is the same size as the AOAV-13 strain from Japan, whereas that of two strains from Ukraine and Kazakhstan are 2080 nt in length, encoding 579 aa. We describe a novel AOAV-13 in migratory waterfowl in China, which suggests that diversified trailer region sequences and HN gene lengths exist within serotype AOAV-13, and highlight the need for its constant surveillance in poultry from live animal markets, and especially migratory birds.

Genetic characterization and phylogenetic analysis of Newcastle disease virus from China

The avian infectious disease, Newcastle disease (ND), caused by Newcastle disease virus (NDV) can cause severe economic losses to poultry whether vaccinated or not in many countries. In this study, a strain of NDV isolated from an outbreak in China was subjected to biological, phylogenetic and genetic characterization. The results showed that the mean death time (MDT) was 52.4 h and the intracerebral pathogenicity indices (ICPI) value was 1.95. In addition, amino acid sequencing result showed that it had a sequence 112R-R-Q-R-R↓F117 at fusion protein cleaving site (FPCS) indicating a velogenic strain. And its genome length is 15,192 nucleotide (nt) with the conserved complementary 3′ leader and 5′ trailer regions encoding six genes, 3′-NP-P-M-F-HN-L-5′. Based on phylogenetic analyses for hyper-variable region and complete genome of F gene, the strain studied here can be clustered into genotype IX, Class II, which has little evolution distance with strains of genotype III, being considered as a transitional strain in the evolution history of NDV. The rescue of infectious cDNA is proceeded in 9-day-old embryonated SPF chicken eggs. Despite the death of the first generation, the allantoic fluid harvested from the first generation lost its pathogenicity after passage. And we found the phenomenon happened due to the antibody appearing in the allantoic fluid. These findings offer our understanding of circulating strains of NDV in China and lay scientific foundations for making more efficient vaccines for Newcastle disease.

The Sub-Chromosomic Macronuclear Pheromone Genes of the Ciliate Euplotes raikovi: Comparative Structural Analysis and Insights into the Mechanism of Expression.

In Euplotes raikovi, we have determined the full-length sequences of a family of macronuclear genes that are the transcriptionally active versions of codominant alleles inherited at the mating-type (mat) locus of the micronuclear genome, and encode cell type-distinctive signaling pheromones. These genes include a 225-231-bp coding region flanked by a conserved 544-bp 5'-leader region and a more variable 3'-trailer region. Two transcription initiation start sites and two polyadenylation sites associated with nonconventional signals cooperate with a splicing phenomenon of a 326-bp intron residing in the 5'-leader region in the generation of multiple transcripts from the same gene. In two of them, the synthesis of functional products depends on the reassignment to a sense codon, or readthrough of a strictly conserved leaky UAG stop codon. That this reassignment may take place is suggested by the position this codon occupies in the transcripts, close to the transcript extremity and far from the poly(A) tail. In such a case, one product is a 69-amino acid protein in search of function and the second product is a 126-amino acid protein that represents a membrane-bound pheromone isoform candidate to function as a cell type-specific binding site (receptor) of the soluble pheromones.

A Rapid Screen for Host-Encoded miRNAs with Inhibitory Effects against Ebola Virus Using a Transcription- and Replication-Competent Virus-Like Particle System.

MicroRNAs (miRNAs) may become efficient antiviral agents against the Ebola virus (EBOV) targeting viral genomic RNAs or transcripts. We previously conducted a genome-wide search for differentially expressed miRNAs during viral replication and transcription. In this study, we established a rapid screen for miRNAs with inhibitory effects against EBOV using a tetracistronic transcription- and replication-competent virus-like particle (trVLP) system. This system uses a minigenome comprising an EBOV leader region, luciferase reporter, VP40, GP, VP24, EBOV trailer region, and three noncoding regions from the EBOV genome and can be used to model the life cycle of EBOV under biosafety level (BSL) 2 conditions. Informatic analysis was performed to select up-regulated miRNAs targeting the coding regions of the minigenome with the highest binding energy to perform inhibitory effect screening. Among these miRNAs, miR-150-3p had the most significant inhibitory effect. Reverse transcription polymerase chain reaction (RT-PCR), Western blot, and double fluorescence reporter experiments demonstrated that miR-150-3p inhibited the reproduction of trVLPs via the regulation of GP and VP40 expression by directly targeting the coding regions of GP and VP40. This novel, rapid, and convenient screening method will efficiently facilitate the exploration of miRNAs against EBOV under BSL-2 conditions.

An RNA polymerase I-driven human respiratory syncytial virus minigenome as a tool for quantifying virus titers and screening antiviral drug

Human respiratory syncytial virus (RSV) is an important pediatric pathogen of lower respiratory tract worldwide. No vaccines and antiviral drugs are available. Herein the use of an RNA polymerase I-driven RSV minigenome for analyzing RSV replication and screening anti-RSV drugs was investigated. The RNA polymerase I (Pol I) was used to transcribe RSV minigenome from the constructed plasmid, designated pHM-RSV-Gluc, of minigenome cDNA which comprised trailer region, gene start sequence (GS), reverse complementary copy of Gaussia luciferase (Gluc) gene, gene end sequence (GE), and leader region in the direction of 5′–3′ end and was flanked by promoter and terminator of Pol I. The expression of Gluc was confirmed in pHM-RSV-Gluc transfected HEp-2 cells following RSV infection and had the characteristics of dose-dependent, which provided a rapid, sensitive, and quantitative method for quantifying virus titers and screening antiviral drugs.

Sequencing and molecular characterization of CTNCEC25, a China fixed rabies virus vaccine strain CTN-1 adapted to primary chicken embryo cells.

BackgroundRabies virus is the main etiologic agent of the widespread neurological disease rabies. Recently, the China rabies virus vaccine strain CTN-1 adapted to chicken embryo cells, which has been designated as CTNCEC25, was obtained and demonstrated to have high immunogenicity. However, the full genome sequence of CTNCEC25 and its phylogenetic relationship with other rabies virus street and vaccine strains have not been characterized.ResultsThe complete genome of CTNCEC25 was sequenced and analyzed. The length of CTNCEC25 genome is 11,924 nucleotides (nt), comprising a 3′ leader sequence of 59 nt, nucleoprotein (N) gene of 1,425 nt, phosphoprotein (P) gene of 989 nt, matrix protein (M) gene of 803 nt, glycoprotein (G) gene of 2,067 nt, RNA-dependent RNA polymerase gene (L) of 6,474 nt and a 5′ trailer region of 71 nt. A comparison of the entire genomes of CTN-1 and CTNCEC25 identified 16 nt substitutions and 1 deletion, resulting in 8 amino acid (aa) changes in the five structural proteins with one in L (aa 1602), two in M (aa 99 and 191) and six in mature G (aa 147, 333, 389, 421 and 485). The percentage homology of the CTNCEC25 genomic sequence with other fully sequenced rabies virus strains ranged from 81.4% to 99.9%. Phylogenetic analysis indicated that CTNCEC25 was more closely related with those recently isolated China street strains than other vaccine strains. Virus growth analysis showed that CTNCEC25 achieved high rate of propagation in cultured cells.ConclusionsIn this study, the complete genome of CTNCEC25 was sequenced and characterized. Our results showed that CTNCEC25 was more closely related to wild street strains circulating in China than other vaccine strains. Sequence analysis showed that the G protein ectodomain amino acid sequence identity between CTNCEC25 and other rabies virus strains was at least 90% identical. Furthermore, CTNCEC25 achieved high virus titers in cultured cells. Given that CTNCEC25 has high immunogenicity and induced strong protective immune response in animals, these results collectively demonstrated that CTNCEC25 is an ideal vaccine strain candidate for producing human vaccine with high quality and safety in China.

Genotype characterization of commonly used Newcastle disease virus vaccine strains of India.

Newcastle disease is an avian pathogen causing severe economic losses to the Indian poultry industry due to recurring outbreaks in vaccinated and unvaccinated flocks. India being an endemic country, advocates vaccination against the virus using lentogenic and mesogenic strains. Two virus strains which are commonly used for vaccination are strain F (a lentogenic virus) and strain R2B (a mesogenic virus). Strain F is given to 0–7 days old chicks and R2B is given to older birds which are around 6–8 weeks old. To understand the genetic makeup of these two strains, a complete genome study and phylogenetic analysis of the F, HN genes of these vaccine strains were carried out. Both the viral strains had a genome length of 15,186 nucleotides and consisted of six genes with conserved complimentary 3' leader and 5' trailer regions. The fusion protein cleavage site of strain F is GGRQGRL and strain R2B is RRQKRF. Although both the viral strains had different virulence attributes, the length of the HN protein was similar with 577 amino acids. Phylogenetic analysis of F, HN and complete genome sequences grouped these two strains in genotype II category which are considered as early genotypes and corroborated with their years of isolation.

Evaluation of the contributions of individual viral genes to newcastle disease virus virulence and pathogenesis.

Naturally occurring Newcastle disease virus (NDV) strains vary greatly in virulence. The presence of multibasic residues at the proteolytic cleavage site of the fusion (F) protein has been shown to be a primary determinant differentiating virulent versus avirulent strains. However, there is wide variation in virulence among virulent strains. There also are examples of incongruity between cleavage site sequence and virulence. These observations suggest that additional viral factors contribute to virulence. In this study, we evaluated the contribution of each viral gene to virulence individually and in different combinations by exchanging genes between velogenic (highly virulent) strain GB Texas (GBT) and mesogenic (moderately virulent) strain Beaudette C (BC). These two strains are phylogenetically closely related, and their F proteins contain identical cleavage site sequences, (112)RRQKR↓F(117). A total of 20 chimeric viruses were constructed and evaluated in vitro, in 1-day-old chicks, and in 2-week-old chickens. The results showed that both the envelope-associated and polymerase-associated proteins contribute to the difference in virulence between rBC and rGBT, with the envelope-associated proteins playing the greater role. The F protein was the major individual contributor and was sometimes augmented by the homologous M and HN proteins. The dramatic effect of F was independent of its cleavage site sequence since that was identical in the two strains. The polymerase L protein was the next major individual contributor and was sometimes augmented by the homologous N and P proteins. The leader and trailer regions did not appear to contribute to the difference in virulence between BC and GBT. This study is the first comprehensive and systematic study of NDV virulence and pathogenesis. Genetic exchanges between a mesogenic and a velogenic strain revealed that the fusion glycoprotein is the major virulence determinant regardless of the identical virulence protease cleavage site sequence present in both strains. The contribution of the large polymerase protein to NDV virulence is second only to that of the fusion glycoprotein. The identification of virulence determinants is of considerable importance, because of the potential to generate better live attenuated NDV vaccines. It may also be possible to apply these findings to other paramyxoviruses.

Thermodynamic Modeling of Variations in the Rate of RNA Chain Elongation of E. coli rrn Operons

Previous electron-microscopic imaging has shown high RNA polymerase occupation densities in the 16S and 23S encoding regions and low occupation densities in the noncoding leader, spacer, and trailer regions of the rRNA (rrn) operons in E. coli. This indicates slower transcript elongation within the coding regions and faster elongation within the noncoding regions of the operon. Inactivation of four of the seven rrn operons increases the transcript initiation frequency at the promoters of the three intact operons and reduces the time for RNA polymerase to traverse the operon. We have used the DNA sequence-dependent standard free energy variation of the transcription complex to model the experimentally observed changes in the elongation rate along the rrnB operon. We also model the stimulation of the average transcription rate over the whole operon by increasing rate of transcript initiation. Monte Carlo simulations, taking into account initiation of transcription, translocation, and backward and forward tracking of RNA polymerase, partially reproduce the observed transcript elongation rate variations along the rrn operon and fully account for the increased average rate in response to increased frequency of transcript initiation.

Complete genome sequence of rabies virus CVS-24 from China

The entire genome of the mouse-adapted rabies virus strain CVS-24 (challenge virus standard 24), was sequenced. The overall length of the genome was 11,927 nucleotide (nt), comprising a leader sequence of 58nt, a nucleoprotein (N) gene of 1353nt, phosphoprotein (P) gene of 894nt, a matrix protein (M) gene of 609nt, a glycoprotein (G) gene of 1575nt, an RNA-dependent RNA polymerase (RdRp, L) gene of 6384nt and a trailer region of 70nt. There was a TGAAAAAAA (TG7) consensus sequence at the end of each gene, except the G gene which had an AGAAAAAAA sequence at the end, and the L/trailer region had the sequence CGAAAAAAA. Three were AACAYYYCT consensus start signals close to TG7. The five cistrons were separated by intergenic regions (IGRS) of 2, 5, 5, 24nt, respectively. Residue 333 of the mature G protein, which is considered to be associated with pathogenicity, was Ala in CVS-24. The topology of the phylogenetic trees generated using N protein sequences suggested that CVS-11 and CVS-N2C have a close relationship to CVS-24.

Characterization of Malaysian velogenic NDV strain AF2240-I genomic sequence: a comparative study

Newcastle disease virus (NDV) strain AF2240 is a viscerotropic velogenic strain that is used as a vaccine challenge virus in Malaysia. The identification of the full length genome will be a crucial platform for further studies of this isolate. In this study, we fully sequenced the genome of a derivative of this strain named AF2240-I. The 15,192 nt long genome contains a 55-nt leader sequence at the 3' whereas the trailer region consists of 114 nt at the 5'. The intergenic sequences between the NP-P, P-M, M-F, F-HN, and HN-L genes comprise 1, 1, 1, 31, and 47 nt, respectively. The acknowledged cleavage site of fusion protein showed amino acid sequence of 112-R-R-Q-K-R-F-117, which corresponds to those of virulent NDV strains. Phylogenetic analysis of the whole virus genome shows that the strain AF2240-I belongs to genotype VIII and is more closely related to velogenic strains QH1, QH4, Fontana, Largo, and Italienas compared to other strains of NDV. Differences are noticed in the hemagglutinin-neuraminidase (HN) and matrix (M) gene between AF2240 and its derivative AF2240-I. This is the first report of a complete genome sequence of an NDV strain isolated in Malaysia.

Molecular characterization of KGH, the first human isolate of rabies virus in Korea

The complete genome sequence of the KGH strain of the first human rabies virus, which was isolated from a skin biopsy of a patient with rabies, whose symptoms developed due to bites from a raccoon dog in 2001. The size of the KGH strain genome was determined to be 11,928 nucleotides (nt) with a leader sequence of 58 nt, nucleoprotein gene of 1,353 nt, phosphoprotein gene of 894 nt, matrix protein gene of 609 nt, glycoprotein gene of 1,575 nt, RNA-dependent RNA polymerase gene of 6,384 nt, and trailer region of 69 nt. Sequence similarity was compared with 39 fully sequenced rabies virus genomes currently available, and the result showed 70.6-91.6 % at the nucleotide level, and 82.8-97.9 % at the amino acid level. The deduced amino acids in the viral protein were compared with those of other rabies viruses, and various functional regions were investigated. As a result, we found that the KGH strain only had a unique amino acid substitution that was identified to be associated either with host immune response and pathogenicity in the N protein, or with a related region regulating STAT1 in the P protein, and related to pathogenicity in G protein. Based on phylogenetic analyses using the complete genome of 39 rabies viruses, the KGH strain was determined to be closely related with the NNV-RAB-H strain and transplant rabies virus serotype 1, which are Indian isolates, and was confirmed to belong to the Arctic-like 2 clade. The KGH strain was most closely related to the SKRRD0204HC and SKRRD0205HC strain when compared with Korean animal isolates, which was separated around the same time and place, and belonged to the Gangwon III subgroup.

Phylogenetic analysis of the bovine parainfluenza virus type 3 from cattle herds revealing the existence of a genotype A strain in China

In 2009, a bovine parainfluenza virus (BPIV3), named as NM09, was isolated using MDBK cell culture from the nasal swabs of normal cattle in China. The NM09 isolate was characterized by RT-PCR and nucleotide sequence analysis. Its complete genome was 15,456 nucleotides in length. Similar to other sequenced PIV strains, the NM09 virus consisted of six non-overlapping genes, which were predicted to encode nine proteins with conserved and complementary 3' leader and 5' trailer regions, conserved gene starts, gene stops, and trinucleotide intergenic sequences. Nucleotide phylogenetic analysis of matrix and hemagglutinin-neuraminidase gene demonstrated that this NM09 isolate belonged to BPIV3 genotype A instead of the previously reported BPIV3 genotype C in China. It is implicated that the different genotypes A and C might coexist infection for a long time in China.

Genotypic and Pathotypic Characterization of Newcastle Disease Viruses from India

Newcastle disease virus (NDV) is an avian paramyxovirus that causes significant economic losses to the poultry industry in most parts of the world. The susceptibility of a wide variety of avian species coupled with synanthropic bird reservoirs has contributed to the vast genomic diversity of this virus as well as diagnostic failures. Since the first panzootic in 1926, Newcastle disease (ND) became enzootic in India with recurrent outbreaks in multiple avian species. The genetic characteristics of circulating strains in India, however, are largely unknown. To understand the nature of NDV genotypes in India, we characterized two representative strains isolated 13 years apart from a chicken and a pigeon by complete genome sequence analysis and pathotyping. The viruses were characterized as velogenic by pathogenicity indices devised to distinguish these strains. The genome length was 15,186 nucleotides (nt) and consisted of six non-overlapping genes, with conserved and complementary 3′ leader and 5′ trailer regions, conserved gene starts, gene stops, and intergenic sequences similar to those in avian paramyxovirus 1 (APMV-1) strains. Matrix gene sequence analysis grouped the pigeon isolate with APMV-1 strains. Phylogeny based on the fusion (F), and hemagglutinin (HN) genes and complete genome sequence grouped these viruses into genotype IV. Genotype IV strains are considered to have “died out” after the first panzootic (1926–1960) of ND. But, our results suggest that there is persistence of genotype IV strains in India.

Molecular characterization of China rabies virus vaccine strain

BackgroundRabies 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.ResultsThe 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.ConclusionsIn 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.