Sheng-Li MENG, Ge-Lin XU, Jia-Xin YAN, Ping-Gang MING, Jie WU, Xiao-Ming YANG, He-Tian MING, Feng-Cai ZHU, Dun-Jin ZHOU, QI-You XIAO and Guan-Mu DONG. Molecular Epidemiology and Sequencing of the G-L Intergenic Region of RabiesViruses Isolated in China[J]. Virologica Sinica, 2007, 22(1): 26-33.
Citation: Sheng-Li MENG, Ge-Lin XU, Jia-Xin YAN, Ping-Gang MING, Jie WU, Xiao-Ming YANG, He-Tian MING, Feng-Cai ZHU, Dun-Jin ZHOU, QI-You XIAO, Guan-Mu DONG. Molecular Epidemiology and Sequencing of the G-L Intergenic Region of RabiesViruses Isolated in China .VIROLOGICA SINICA, 2007, 22(1) : 26-33.

中国分离的狂犬病毒G-L基因间隔区的序列测定与分子流行病学

  • 通讯作者: 严家新, yanjx45@163.com
  • 收稿日期: 2006-07-31
    录用日期: 2006-08-31
  • 摘要:在2004-2006年期间,从中国不同地区共收集到25株狂犬病毒,其中24株是从犬脑分离,1株是从人脑分离。对这25株街毒株和7个代次的CTN疫苗株的G-L基因间隔区进行了序列测定和系统进化关系分析。结果显示:所有中国毒株的G-L基因间隔区都由519个碱基构成,这25株街毒株之间核苷酸序列同源性为95.5%-100%。本文所分析的所有中国毒株都属于基因 1型,其详细分布情况与地理上的来源相关。这些毒株共可分为密切相关的两个类群,所有街毒株为一个类群,不同代次的CTN疫苗株为另一个类群。此项研究提供了狂犬病毒基于G-L基因间隔区序列的分子流行病学详细资料。

Molecular Epidemiology and Sequencing of the G-L Intergenic Region of RabiesViruses Isolated in China

  • Corresponding author: Jia-Xin YAN, yanjx45@163.com
  • Received Date: 31 July 2006
    Accepted Date: 31 August 2006

    Fund Project: The Key Technologies R & D Programme of the 10th National five-year-plan 2004BA718 b03

  • A group of 25 rabies viruses (RABVs), recovered from 24 dogs and one human case, were collected from various areas in China between 2004 and 2006. Genetic and phylogenetic analyses of the G-L intergenic region were carried out in 25 street RABV isolates and CTN vaccine strains of 7 generations. The study was based on the comparison of a 519 bp nucleotide sequence, encompassing the G-L intergenic region. The nucleotide sequence homologies of Chinese street strains were from 95.5% to 100%. The phylogenetic analysis showed that all Chinese isolates clearly supported the placement of all Chinese viruses in Lyssavirus genotype 1 and they were distributed according to their geographical origins. All of the Chinese strains were closely related but they could still be divided into two groups: group of street strains and group of CTN strains. This study presents details about the molecular epidemiology of rabies viruses based on the sequences of the G-L Intergenic region.

  • 加载中
    1. Badrane H,Bahloul C,Perrin P,et al. 2001. Evidence of two Lyssavirus phylogroups with distinct pathogenicity and immunogenicity[J]. J Virol,75(7): 3268-3276.
        doi: 10.1128/JVI.75.7.3268-3276.2001

    2. Dodet B,Meslin F-X,Aubert. 2001. In: Rabies Control in Asia[M]. Paris. John Libbey Eurotext,191-196.

    3. Hillis D M,Bull J J. 1993. An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis[J]. Systematic Biol,42: 182-192
        doi: 10.1093/sysbio/42.2.182

    4. Kimura M. 2004. A simple method for estimating evolutionary rates of base substitutions through compara-tive studies of nucleotide sequences[J]. Briefings in Bioinformatics,5: 150-163
        doi: 10.1093/bib/5.2.150

    5. Kumar S,Tamura K,Nei M. 2004. MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment[J]. Brief Bioinform,5(2): 150-163.
        doi: 10.1093/bib/5.2.150

    6. Kuzmin I V,Hughes G J,Botvinkin A D,et al. 2005. Phylogenetic relationships of Irkut and West Caucasian bat viruses within the Lyssavirus genus and suggested quantitative criteria based on the N gene sequence for lyssavirus genotype definition[J]. Virus Res,111(1): 28-43.
        doi: 10.1016/j.virusres.2005.03.008

    7. Meslin F X,Kaplan M M,Koprowski H,et al. 1996. Labo-ratory techniques in rabies[M]. 4th ed. Geneva. World Health Organization

    8. Nadin-Davis S A,Abdel-Malik M,Armstrong J,et al. 2002. Lyssavirus P gene characterisation provides insights into the phylogeny of the genus and identifies structural similarities and diversity within the encoded phosphoprotein[J]. Virology,298(2): 286-305.
        doi: 10.1006/viro.2002.1492

    9. Nel L H,Sabeta C T,von Teichman B,et al. 2005. Mongoose rabies in southern Africa: a re-evaluation based on mo-lecular epidemiology[J]. Virus Res,109(2): 165-173.
        doi: 10.1016/j.virusres.2004.12.003

    10. Ravkov E V,Smith J S,Nichol S T. 1995. Rabies virus glycoprotein gene contains a long 3' noncoding region which lacks pseudogene properties[J]. Virology,206(1): 718-723.
        doi: 10.1016/S0042-6822(95)80095-6

    11. Sato G,Itou T,Shoji Y,et al. 2004. Genetic and phylogenetic analysis of glycoprotein of rabies virus isolated from several species in Brazil[J]. J Vet Med Sci,66(7): 747-753.
        doi: 10.1292/jvms.66.747

    12. Thompson J D,Higgins D G,Gibson T. 1994. CLUSTALW: improving the sensitivity of progressive multiple sequence alignment through sequence weighting,position-specific gap penalties and weight matrix choice [J]. Nucleic Acids Res,22: 4673-4680
        doi: 10.1093/nar/22.22.4673

    13. Velazquez-Monroy O,Vargas-Pino F,Gutierrez-Cedillo V,et al. 2003. Advances in canine rabies control in Mexico[D]. The XVI international conference-Rabies in the Americas Philadelphia: Thomas Jefferson University. p78.

    14. WHO. 2005. Expert Consultation on rabies[J]. World Health Organ Tech Rep Ser,931: 1-88.

    15. Xu G L,Ku L,Wu J,et al. 2002. Sequence analysis of N Gene among 19 rabies virus street strains from China[J]. Vi-rologica Sinica,18: 48-51.(in Chinese)

    16. Zhang Y Z,Xiong C L,Zou Y,et al. 2006. Molecular cha-racterization of rabies virus isolates in China during 2004[J]. Virus Res,121(2): 179-188.
        doi: 10.1016/j.virusres.2006.05.010

  • 加载中

Figures(1) / Tables(4)

Article Metrics

Article views(3529) PDF downloads(20) Cited by(0)

Related
Proportional views
    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Molecular Epidemiology and Sequencing of the G-L Intergenic Region of RabiesViruses Isolated in China

      Corresponding author: Jia-Xin YAN, yanjx45@163.com
    • 1. Wuhan Institute of Biological Products,Wuhan 430060,China
    • 2. Fuyang Center for Disease Control and Prevention,Fuyang 236000,China
    • 3. Jiangsu Center for Disease Control and Prevention,Nanjing 210009 China
    • 4. Wuhan Center for Disease Control and Prevention,Wuhan 430015,China
    • 5. Hunan Center for Disease Control and Prevention,Changsha 410005,China
    • 6. National Institute for the Control of Pharmaceutical and Biological Products,Beijing 100050,China
    Fund Project:  The Key Technologies R & D Programme of the 10th National five-year-plan 2004BA718 b03

    Abstract: A group of 25 rabies viruses (RABVs), recovered from 24 dogs and one human case, were collected from various areas in China between 2004 and 2006. Genetic and phylogenetic analyses of the G-L intergenic region were carried out in 25 street RABV isolates and CTN vaccine strains of 7 generations. The study was based on the comparison of a 519 bp nucleotide sequence, encompassing the G-L intergenic region. The nucleotide sequence homologies of Chinese street strains were from 95.5% to 100%. The phylogenetic analysis showed that all Chinese isolates clearly supported the placement of all Chinese viruses in Lyssavirus genotype 1 and they were distributed according to their geographical origins. All of the Chinese strains were closely related but they could still be divided into two groups: group of street strains and group of CTN strains. This study presents details about the molecular epidemiology of rabies viruses based on the sequences of the G-L Intergenic region.

    • Rabies is an acute,progressive,incurable viral encephalitis,caused by a single stranded RNA virus belonging to the genus Lyssavirus of the family Rhabdoviridae. Human mortality from endemic canine rabies is estimated to be 55 000 deaths per year with 56 % of the deaths estimated to occur in Asia and 44% in Africa (14). While many distinct RABV variants are harbored by a variety of mammalian host species,the main transmission route for human rabies in the developing world is a rabid dog bite: between 94% and 98% of human rabies deaths are due to canine rabies (2).

      In China,there were 110 983 human rabies cases and three major epidemic peaks between 1950 and 2005. These human cases occurred mainly in poor,rural areas and were transmitted by dogs. The first epidemic occurred in the mid 1950s when cases rose to a peak of about 2 000 annually. After a decline in the 1960s,cases again started to increase in the early 1970s reaching a peak in 1982 but remaining at levels of 5-6 000 cases per year till the end of the decade (16). Since the vast majority of cases have been due to canine rabies,an extensive dog vaccination and culling program was initiated in 1987 and between the years 1990 and 1996 canine rabies decreased. It is important to note that 98% of the people who died of rabies in China had not received any anti-rabies vaccine treatment and 70% of the people who died of rabies between 1996 and 2002 lived in Guangxi,Jiangxi,Hunan,Guangdong and Jiangsu provinces. There were 2 651 and 2 571 cases in 2004 and 2005,respectively,which suggest the new epidemic peak is coming. All warm-blooded mammals are susceptible to RABV infection,but dogs represent the major host of rabies virus in China. The main transmission route for rabies is canine bites. 95% rabies deaths are due to canine bites and rabid cats are responsible for 4% human cases in China. Increased dog populations and limited vacci-nation against rabies together with improved transportation may have resulted in this recent increase in rabies cases.

      The lyssavirus particle has a bullet-shaped form,100-300 nm in length and 75 nm in diameter (7). The RABV genome consists of a single-stranded,unsegmented,negative-sense RNA of about 12 kb,which encodes five viral proteins (3'-N-P-M-G-L-5') and is contained in a bullet-shaped and bilayered envelope. The RNA polymerase(L) and phos-phoprotein(P) complex with the nucleoprotein (N) form the nucleocapsid (NC),and the matrix protein(M) and the glycoprotein (G) form the inner and outer layers of the envelope,respectively.

      Evolutionary studies of lyssaviruses have tended to focus on the N protein,a well conserved structural protein,and the G protein which forms the envelope and capsid of the virus. The G protein contains domains responsible for host cell receptor recognition and membrane fusion and is the major target for the host neutralizing-antibody response (1) ,and the G-L intergenic region or pseudogene (Ψ) which is a noncoding region,highly variable,not subject to immunological selective pressure,and has therefore been used in studies of molecular epidemiology of the rabies virus (11).

      Sequence comparison of several regions of the genome,including the N,P,and G genes,has provided consistent evidence for division of the Lyssavirus genus into two major phylogroups and seven established genotypes as follows: Phylogroup Ⅰ comprises Rabies virus(RABV; genotype 1),Duvenhage virus(DUVV; genotype 4),European bat lyssavirus 1 (EBLV-1; genot-ype 5),European bat lyssavirus 2 (EBLV-2; genotype 6),Australian bat lyssavirus(ABLV; genotype 7); Phylogroup Ⅱ contains Lagos bat virus(LBV; genotype 2),Mokola virus,(MOKV; genotype 3) (1, 8). Four recent lyssavirus isolates from bats of Eurasia,designated Aravan virus(ARAV),Khujand virus(KHUV),Irkut virus(IRKV) and West Caucasian bat virus(WCBV),have been described (6). While these viruses await formal classification,their genetic diversity between each other and currently classified lyssaviruses would suggest that they represent several additional genotypes.

      By analyses of the nucleotide sequence data of the non-coding G-L intergenic region,this study provides a new description of the molecular epidemiology of rabies in China and a regional comparison of virus relationships.

    • Between 2004 and 2006,brain samples from 24 dogs and one human were collected from 5 provinces in China: Anhui,Hubei,Jiangsu,Hunan,and Jiangxi provinces. The CTN vaccine strains of 7 generations from 7 to 35 passages were also included in the analysis. The species origin,year of isolation and geographical distribution of the rabies virus samples used in this study are described in and Table 1.

      Table 1.  List of rabies virus Chinese strains employed for evolutionary analysis

    • All the brain samples were tested for rabies virus by direct immunofluorescence antibody(IFA) test and enzyme-linked immunosorbent assay (ELISA),which are the standard diagnostic tests currently used in China. The anti-rabies nucleoprotein monoclonal antibodies and FITC-labelled monoclonal antibody were manufactured in Wuhan Institute of Biological Products,China.

    • Total RNA was extracted by TRIzol(Invitrogen,Carlsbad,CA,USA),from brain tissue of naturally infected animals for PCR assay according to the manufacturer's instructions.

    • Reverse transcription(RT) of viral RNA was Performed in a 100µL reaction volume using 1 µL extracted RNA,10units AMV Reverse Transcriptase (Promega),10 units RNasin ribounclease inhibitor,2 µL RT reaction buffer (250 mmol/L Tris-HCl [pH 8.3,25℃],250 mmol/L KCl,50 mmol/L MgCl2,2.5mmol/L spermidine and 50mmol/L ditheothreitol),250mmol/L for each of the four deoxynucleotides,and 30 pmol sense primerGH3.3(5'-GAYTACACCATC-TGGATGCC-3') corresponding to nucleotides 3894-3913 of the PV genome. RT reactions were incubated at 42℃ for 90 min before heat-inactivating AMV Reverse Transcriptase at 100℃ for 3min. The cDNA product was stored at -20℃ or PCR amplification was done at once.

    • In a 100µL eppendorf tube,add the following components in the indicated order: 37µL dd H2O,1µL sequence-specific sense primer MSL(5'-TGGATTTGTGGATKAAAGAGGC-3'),and 1µL anti-sense primer L1(5'-GAGTTNAGRTT GTART-CAGAG-3') (30.3 pmol/µL),corresp-onding to bases 3995-4016 and 5516-5536 of the PV genome,respectively,2 µL template cDNA,5µL 10×PCR buffer(100mmol/L Tris-HCl[pH9.0],500 mmol/L KCl,0.5% Tween-20),3µL 20mM MgCl2,Mix gently and collect drops by brief centrifugation,then incubate the mixture at 100℃ for 3min. For collecting drops,chill on ice and spin down by brief centrifugation,then place the tube on ice and add 0.5µL Taq DNA polymerase(Generay Biotech Shanghai,5U/µL). Collect drops by brief centrifugation,then PCR amplify. Thermal cycling conditions were: 1 cycle at 94℃ for 1min; 40 cycles at 94℃ for 1min,at 58℃ for 50s and at 72℃ for 90s,followed by 1 final incubation at 72℃ for 10 min; holding at 4℃. The PCR product was analyzed on a 1% agarose gel containing ethidium bromide.

    • PCR products were purified using the quick-spin PCR Purification Kit and sequenced using the Taq Big Dye Terminator Cycle Sequencing Ready Reaction Kit on an Applied Biosystems 3770 DNA automated sequencer(Applied Biosystems Inc. Foster City,CA,USA). A multiple alignment of the nucleotide sequence of G-L intergenic region (519 bp) was generated with ClustalW1.8 (12). Distance calcula-tions were done using the Kimura 2-parameter model (4). At the same time,for construction of the phylogenetic trees,the Neighbour Joining(NJ) and Maximum Likelihood (ML) method were conducted using MEGA version 3.1 (3, 5). The statistical significance of the phylogenies constructed was estimated by bootstrap analysis with 1000 replicates,and bootstrap values above 70% were considered significant. The new nucleotide sequences obtai-ned in this study have been submitted to GenBank and their accession numbers are listed in Table 1.

    • Between 2004 and 2006,from the total of 381 dogs,80 cats,100 bats and a human brain samples examined,the rabies infection were confirmed in 24 dogs and one human brain samples (Table 1) by FAT,ELISA and RT-PCR.

      The G-L intergenic region of 25 street rabies virus isolates from China and CTN vaccine strains of 7 gene-rations were amplified by the RT-PCR and sequenced. A 519 nucleotides sequence of the non-coding G-L intergenic region was analyzed with the related sequences of other rabies viruses obtained from GenBank (Table 2). The sequence homologies of the 25 Chinese strains against other non-Chinese rabies virus strains were at least 67.5%,and the pairwise distances between strains were calculated using the Kimura 2-parameter model in MEGA 3.1 (data not shown). Analysis of the sequences within the G-L intergenic region showed that street rabies viruses isolated from China had the closest relationship to the Thailand strain (THA1-HM) and the Malaysia strain (MAL1-HM). The homologies were at least 82% and 82.3%,respectively. However,there were higher nucleotide homologies of the sequences of the G-L intergenic region among all rabies virus isolates from China,ranging from 95.5% to 100%. In addition,alignment of the G-L intergenic region of Chinese street strains showed that nucleotide homologies were ranging from 85.5% to 86.4%,68.5% to 69.8%,68.1% to 69.8%,70.7% to 73.2% and 69.9% to 72.8% as compared with 4tive results of nucleotide homologies are shown in Table 3.

      Table 2.  Sequences of rabies virus strains obtained from GenBank in this G-L intergenic region analysis

      Table 3.  Nucleotide homology (%) and Std. Err. of the G-L intergenic regions of some rabies viruses

      For a better understanding of phylogeny of Chinese strains,pairwise comparisons of the G-L intergenic region were performed. One consensus tree (not shown) was constructed with the NJ method in rabies viruses isolated from China,showing the Chinese strains segregated into two groups (ⅠandⅡ) with 100% bootstrap support. The CTN vaccine strains formed GroupⅡand the Chinese street strains composed of GroupⅠ,which was further divided into four lineages (Ⅰa -Ⅰd) with at least 80% bootstrap supports and contained rabies viruses originating from different areas of China. The Ⅰa Ⅰd lineage originated from Anhui,Hubei,Jiangsu and Hunan provinces,respectively. As in the Neighbor joining phylogenetic (NJ) tree,the NC specimens lay on a separate branch within groupⅠ. The nucleotide sequence differences of Chinese street strains are shown in Table 4. The 3 regions(8bp to 49bp,254bp to 3 08bp,and 505bp to 514bp) are highly diverse. A phylogenetic NJ tree of interior branch test based on alignment of the G-L intergenic region nucleotide sequences and the neighbor-joining method is shown in Fig. 1. A similar ML tree was obtained by the maximum likelihood method (data not shown).

      Table 4.  The comparison of the G-L intergenic region nucleotide sequences of Chinese street strains

      Figure 1.  Neighbor joining phylogenetic tree based on the G-L intergenic region nucleotide sequences of rabies virus by Kimura two parameters (MEGA version 3.1).The length of the horizontal branches reflects phylogenetic distance relationship. Percentage bootstrap values above 70% are shown at the branch nodes.

    • Since the SARS epidemic of 2003,the Chinese government has set up a systematic surveillance net for zoonotic diseases. Increased surveillance together with increased dog populations may be the principal factors explaining the increasing number of human rabies cases reported in China in the last few years. A mass vacci-nation of dogs in the epidemic regions,as has been successfully undertaken in certain countries of Latin Ame-rica (13) ,could virtually eliminate rabies,but the high cost of vaccines currently prevents this undertaking. The data reported here represent an early step in understanding the epidemiology of dog rabies in China. This knowledge will be important to future control efforts.

      According to the phylogenetic analysis of the G-L intergenic region in different rabies virus strains,all of the Chinese street isolates and the CTN vaccine strains of 7 generations were closely related but could be divided into two groups. The Chinese street isolates formed Group Ⅰ,which was further divided into four lineages related to origins of distinct geographic regions. Inte-restingly,we recovered 15 RABV isolates from Fuyang and its surrounding area,and although the homology among these 15 isolates,which comprise lineage 1a,is very high,these isolates could still be further subdivided into several variants. Group Ⅱ included CTN vaccine strains of different generations,which were originally isolated from Shangdong province in 1956 and became a vaccine strain through serial passages in mouse brain and human diploid cell(KMB-17). Group Ⅰ and Group Ⅱ were closely related. A previous study on the N genes of Chinese strains has obtained similar phyloge-netic patterns (15).

      The G-L intergenic region in different rabies virus strains was highly variable and even the nucleotide numbers of this region were different. In Chinese street strains,this region had 519 nucleotides and only one transcription termination and one polyadenyla-tion(PPT) motif. However,there were two transcrip-tion termination and two polyadenylation(PPT) motifs found in this region in the Pasteur strain derived from an 1882 French strain and the SAD strain derived from a 1935 U.S. rabies isolate (10).

      The G-L intergenic region is a non-coding region,which is highly susceptible to random mutations,unrestricted by structure and function requirements or by immunological pressure. Being the most divergent region of the rabies genome,it might be more sensi-tive in demonstrating recent evolutionary events (9).

      The CTN,PV and PM strains are the human vaccine strains currently used in China. The nucleotide homology between the CTN strains and the Chinese street strains was much higher(at least 12.3%) than that of any other vaccine strains. When the N and G Genes of the same vaccine strains and the same Chinese street strains were sequenced and analyzed,we obtained very similar phylogenetic trees(the data are not showed here and they will be published later). So it is possible that vaccines produced with the CTN vaccine strain might be more efficient than that produced with any other vaccine strains in China. However,vaccine trials should be undertaken to confirm this hypothesis.

      At the same time,we find that SAD,PV,HEP-Flury and PM vaccine strains have very high nucleotide homologies with South Korean,Indian,Canadian and Brazilian isolates. All these vaccine strains originated from America or Europe. For nearly half a century there was virtually no connection between China and South Korea and there is natural barrier between India and Tibet. These might explain the difference in the strains between China and South Korea or India.

      In conclusion,the findings from this study confirmed the fact that rabies viruses are genetically variable and the evolutions of different groups are apparently dominated by geographical influences. The nucleotide homology between the CTN vaccine strain and the Chinese street strains is higher than for other vaccine strains. Our study provides new data to support the national rabies control programme.

    Figure (1)  Table (4) Reference (16) Relative (20)

    目录

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return