Emerg Infect DisEmerging Infect. DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention16318702336675105-051310.3201/eid1110.050513DispatchDispatchGenetic Characterization of Nipah Virus, Bangladesh, 2004Nipah Virus Genetic Characterization, BangladeshHarcourtBrian H.*LoweLuis*TaminAzaibi*LiuXin*BankampBettina*BowdenNadine*RollinPierre E.*ComerJames A.*KsiazekThomas G.*HossainMohammed JahangirGurleyEmily S.BreimanRobert F.*BelliniWilliam J.*RotaPaul A.*Centers for Disease Control and Prevention, Atlanta, Georgia, USAICDDR,B: Centre for Health and Population Research, Dhaka, BangladeshAddress for correspondence: Paul A. Rota, Centers for Disease Control and Prevention, MSC22, 1600 Clifton Rd, Atlanta, GA, 30333 USA; fax: 404-639-4187; email: prota@cdc.gov102005111015941597

Until 2004, identification of Nipah virus (NV)-like outbreaks in Bangladesh was based on serology. We describe the genetic characterization of a new strain of NV isolated during outbreaks in Bangladesh (NV-B) in 2004, which confirms that NV was the etiologic agent responsible for these outbreaks.

Keywords: Nipah virusencephalitishenipavirusdispatch

Nipah virus (NV) and Hendra virus (HV), the only members of the genus Henipavirus within the family Paramyxoviridae are different from most other paramyxoviruses because they have a broad host range in vivo and in vitro. Although HV and NV have genetic characteristics and replication strategies similar to those of other paramyxoviruses, the henipaviruses have several unique genetic features (1,24).

The first outbreak of NV occurred between late 1998 and early 1999 in peninsular Malaysia and Singapore and was associated with respiratory disease in swine and acute and febrile encephalitis in humans. Direct contact with sick pigs was the primary source of human infection (1). Of 265 human cases, 108 were fatal. Although NV is excreted in respiratory secretions and urine of patients (1,5), a survey of healthcare workers in Malaysia showed no evidence of human-to-human transmission (6). The reservoir of NV is presumed to be fruit bats, primarily of the genus Pteropus (7,8), and humans are infected through intermediate hosts such as pigs.

Recently, NV has been established as the cause of fatal, febrile encephalitis in human patients in Bangladesh during the winters of 2001, 2003, and 2004 (912). An NV-like virus was identified as the cause of the outbreaks in 2001 and 2003 on the basis of serologic testing (12). Two outbreaks consisting of 48 cases of NV were detected in 2004 in 2 adjacent districts (30 km apart) of central Bangladesh (Rajbari and Faridpur) with a case-fatality rate of nearly 75%. Because of heightened surveillance, other small clusters and isolated cases (n = 19) were identified during the same period in 7 other districts in central and northwest Bangladesh. Although antibodies to NV were detected in fruit bats from the affected areas in 2004, an intermediate animal host was not identified, which suggests that the virus was transmitted from bats to humans. Human-to-human transmission of NV was also documented during the Faridpur outbreak (10,11). We describe the genetic characteristics of 4 NV isolates from the outbreak in Bangladesh in 2004.

The Study

Virus isolation was performed in the BSL-4 laboratory at the Centers for Disease Control and Prevention in Atlanta. Vero E6 cells were inoculated and observed for characteristic cytopathic effect, syncytium formation (1,5). NV was isolated from 2 oropharyngeal swabs (SPB200401066, SPB200406506), 1 cerebrospinal fluid (SPB200401617), and 1 urine specimen (SPB200405758) from human patients, and isolation was confirmed by reverse transcription polymerase chain reaction (RT-PCR). Two isolates were from Rajbari, and 1 was from Faridpur; the fourth isolate, from the Rajshahi district (100 km from Rajbari), was not linked to the other 2 outbreaks. The complete genomic sequence of the first viral isolate (SPB200401066) from Rajbari was derived and submitted to GenBank (accession no. AY988601) as NV-Bangladesh (NV-B). The sequences of the open reading frame (ORF) coding for nucleoprotein (N) were obtained for the other 3 isolates. The methods used for RT-PCR, sequencing, cDNA cloning, rapid amplification of cDNA ends (RACE), and sequence analysis were previously described (2,3).

The genome of NV-B is 18,252 nt in length, 6 nt longer than NV-Malaysia (NV-M), the prototype strain of NV (SPB199901924). The additional 6 nt map to the 5´ nontranslated region of the fusion protein (F) gene. The length of the NV-B genome is evenly divisible by 6, suggesting that NV-B follows the "rule of six" (3). The gene order and sizes of all the ORFs except V are conserved between NV-B and NV-M (Table, Figure, A). The overall nucleotide homology between the genomes of NV-B and NV-M is 91.8%, but the changes are not uniformly distributed throughout the genome. Nucleotide homologies are higher in the protein coding regions than in the noncoding regions, although the sizes of the nontranslated regions remain highly conserved (Table). The predicted amino acid homologies between the proteins expressed by NV-M and NV-B are all >92% (Table).

Comparison of gene sequences NV-Malaysia and NV-Bangladesh*
GeneVirusOpen reading frame
5´ nontranslated
3´ nontranslated
Length†% amino acid identity‡% nucleotide homology‡Length§% nucleotide homology‡Length§% nucleotide homology‡
NMalay53298.394.35710058690.8
Bang53257586
PMalay70992.092.010591.446988.1
Bang709105469
VMalay5292.595.7
Bang55
WMalay4710098.5
Bang47
CMalay16695.297.6
Bang166
MMalay35298.993.410086.020083.5
Bang352100200
FMalay54698.493.428483.141279.4
Bang546290412
GMalay60295.593.023375.550480.8
Bang602233504
LMalay2,24498.293.415382.46780.6
Bang2,24415367

*NV, Nipah virus.
†Length in amino acids.
‡Percentage amino acid identity or nucleotide homology after sequences were aligned by using GAP from GCG.
§Length in nucleotides. NV-Malaysia gene lengths and 5´ and 3´ nontranslated sequences and gene lengths were obtained from GenBank accession no. AF212302. The percentage amino acid identity does not take into consideration conserved amino acid changes.

A) Schematic representation of the genome of Nipah virus (NV). Negative-sense genomic RNA is shown in 3´ to 5´ orientation. Open reading frames (ORFs) are indicated by shaded boxes: N, nucleoprotein; P, phosphoprotein; M, matrix protein; F, fusion protein; G, attachment protein; L, polymerase protein. B) Phylogenetic analysis of the N ORFs from members of the subfamily Paramyxovirinae. Arrows identify the 5 genera. A phenogram of the N ORFs of members of this subfamily was created by using maximum parsimony analysis with PAUP 4.02 (Sinauer Associates, Sunderland, MA, USA). Abbreviations and accession numbers: HPIV-1, human parainfluenza virus, D01070; Sendai, X00087; HPIV-3, D10025; CDV, canine distemper virus, AF014953; PDV, phocine distemper virus, X75717; RPV, Rinderpest virus, X68311; MV, K01711; DMV, dolphin Morbillivirus, X75961; NDV, Newcastle disease virus, AF064091; GP, goose paramyxovirus, AF473851; HPIV-4b, M32983; HPIV-4a, M32982; Tioman, AF298895; Menangle, AF326114; HPIV-2, M55320; Simian virus 5 (SV5), M81442; Mumps, D86172; NV-UMCC1, AY029767; NV-Malaysia, AF212302; NV-P. hypomelanus,* AF376747; NV-Bangladesh; Hendra virus, AF017149; Tupaia paramyxovirus, AF079780; Mossman virus, AY286409; and Salem virus, AF237881. C) The phylogenetic relationship between the N gene sequences of the 4 human NV isolates from the Bangladesh outbreak in 2004 and the N gene sequences from pig and human NV isolates from Malaysia. Accession numbers for the pig isolates of NV (13) are AJ627196, AJ564622, and AJ564621. *NV-P. hypomelanus is sequence from a virus isolated from Pteropus hypomelanus, the Island Flying Fox (8).

Overall, the predicted amino acid homologies of the surface glycoproteins, F and G, of NV-B and NV-M are high (Table). In the F protein, the predicted cleavage site, F1 amino-terminal domain, transmembrane domain, and predicted N-glycosylation sites are identical in NV-B and NV-M (2). Four of the 9 predicted amino acid changes occur in the first 11 amino acids (aa) of the precursor of the F protein, F0, which fall within the predicted signal peptide and would be cleaved from the mature protein. Within the G proteins, the predicted transmembrane domains, and the positions of all 17 cysteine residues are conserved between NV-B and NV-M. Of 8 predicted N-linked glycosylation sites in the G protein of NV-M, 6 are conserved in NV-B and in HV (2).

The coding strategy of the P gene is identical in NV-B and NV-M. In these viruses, the P gene contains the C, V, and W ORFs in addition to the P ORF. Like most other paramyxoviruses, the henipaviruses have a conserved AG-rich region that acts as an editing site to facilitate the addition of nontemplated G residues into the transcripts of the P gene. The edited transcripts encode 2 proteins, V and W, which are co-amino-terminal with P but have unique carboxy termini (2). The addition of 1 G residue generates the mRNA for the V protein, and the addition of 2 G residues produces the mRNA for the W protein. Sequence analysis of multiple cDNA clones containing the editing site of NV-B identified edited transcripts that encoded both the V and W proteins (data not shown). The conserved 20-nt region encompassing the editing site is identical in NV-M, HV, and measles virus (MV) (2); however, the editing of site of NV-M (UGGGUAAUUUUUCCCGUGUC) differs from NV-B (GGGAUAAUUUUUCCCGUGUC) at 2 nt positions (underlined). The functional significance of these 2 substitutions is under investigation. All of the cysteine residues are conserved in the V proteins of NV-B and NV-M; however, the unique portion of the V protein of NV-B is predicted to be 55 aa, 3 aa longer than the V protein of NV-M. The predicted W protein of NV-B is identical in size and sequence to the W protein of NV-M. Recently, aa 100–160 and 230–237 of the V protein of NV-M have been identified as necessary for inhibition of interferon signaling (14). These regions are highly conserved in the V protein of NV-B, which has 4 predicted amino acids substitutions (1 conservative) between positions 100–160 and no predicted substitutions between positions 230–237. In addition, the V protein of NV-B has 3 predicted amino acid substitutions in the CRM1-dependent nuclear export signal that was identified between aa 174–193 in the V protein of NV-M (14). The biologic effects of these amino acid substitutions are under investigation.

The L proteins of NV-B and NV-M had a high level of predicted amino acid conservation (Table). The 6 highly conserved domains of viral polymerases, originally described by Poch et al. (15) and delineated for NV-M by Harcourt et al. (3), remain largely unchanged between NV-B and NV-M. Domains 1, 2, and 5 have 1 conservative amino acid change each and domain 3, which is considered the most conserved domain within the L proteins of paramyxoviruses, has 2 aa changes. The 4 motifs identified in domain 3, including the QGDNE motif, which is assumed to be the active site of the polymerase, are identical between the NV-B and NV-M, as is the predicted nucleotide-binding motif in domain 6.

The cis-acting control sequences are highly conserved in the genomes of NV-B and NV-M. As in NV-M, the intergenic sequences in NV-B are GAA, with the exception of the sequence between the G and L genes, UAA, which is unique among the henipaviruses. However, the intergenic sequence between the G and L genes is GAA in the second isolate from Bangladesh. The transcriptional start and stop signals of each gene of NV-B are highly conserved in relation to the other henipaviruses. The 3´ leader sequence of NV-B is identical in length to those of all other paramyxoviruses and has nucleotide changes at positions 14 and 47 compared to NV-M. The 5´ trailer of NV-B is identical in length and sequence to NV-M.

Phylogenetic analysis was used to compare the sequence of the N ORF of NV-B to the sequences of the N ORFs from other members of the subfamily Paramyxovirinae. The results confirmed the results of the sequence comparisons, which show that NV-B is most closely related to the henipavirus NV-M, and support the conclusion that NV-B should be regarded as new strain of NV (Figure, B). Phylogenetic analyses conducted with the sequences of the other genes produced similar results (data not shown). The sequences of the N ORFs of 4 NV isolates from Bangladesh share 99.1% nt homology (Figure, C) but exhibited more interstrain nucleotide heterogeneity than the sequences of the human isolates in Malaysia, which were nearly identical (1,8,13). These varying amounts of genetic variability may reflect differences in the mode of transmission of NV in the 2 countries. In Malaysia, molecular evidence suggests that at least 2 introductions of NV into pigs occurred (Figure, C). However, the nearly identical sequences of human and pig isolates from the later phase of the outbreak suggest that only 1 of the variants spread rapidly in pigs and was associated with most human cases (13). In contrast, the sequence heterogeneity observed in Bangladesh may be the result of multiple introductions of NV into humans from different colonies of fruit bats.

Conclusions

This first look at strain variation in NV indicates that viruses circulating in different areas have unique genetic signatures and suggests that these strains may have coevolved within the local natural reservoirs. Until 2004, identification of NV outbreaks in Bangladesh had been based only on serologic testing. The isolation and genetic characterization of NV-B confirm that NV was the etiologic agent responsible for these outbreaks.

Note: After this article was accepted for publication, Nipah virus was isolated from Pteropus lylei in Cambodia (16). Phylogenetic analysis of the N gene sequences demonstrated that this virus is more closely related to Nipah-Malaysia than to Nipah-Bangladesh and represented another lineage of Nipah virus.

Suggested citation for this article: Harcourt BH, Lowe L, Tamin A, Liu X, Bankamp B, Bowden N, et al. Genetic characterization of Nipah virus, Bangladesh, 2004. Emerg Infect Dis [serial on the Internet]. 2005 Oct [date cited]. http://dx.doi.org/10.3201/eid1110.050513

Dr. Harcourt is a research associate in the Division of Viral and Rickettsial Diseases, National Center for Infectious Diseases, Centers for Disease Control and Prevention. His research interests include viral pathogenesis and emerging viral diseases.

ReferencesChua KB, Bellini WJ, Rota PA, Harcourt BH, Tamin A, Lam SK, Nipah virus: a recently emergent deadly paramyxovirus. Science. 2000;288:14325 10.1126/science.288.5470.143210827955Harcourt BH, Tamin A, Ksiazek TG, Rollin PE, Anderson LJ, Bellini WJ, Molecular characterization of Nipah virus, a newly emergent paramyxovirus. Virology. 2000;271:33449 10.1006/viro.2000.034010860887Harcourt BH, Tamin A, Halpin K, Ksiazek TG, Rollin PE, Bellini WJ, Molecular characterization of the polymerase gene and genomic termini of Nipah virus. Virology. 2001;287:192201 10.1006/viro.2001.102611504554Wang LF, Harcourt BH, Yu M, Tamin A, Rota PA, Bellini WJ, Molecular biology of Hendra and Nipah viruses. Microbes Infect. 2001;3:27987 10.1016/S1286-4579(01)01381-811334745Chua KB, Lam SK, Goh KJ, Hooi PS, Ksiazek TG, Kamarulzaman A, The presence of Nipah virus in respiratory secretions and urine of patients during an outbreak of Nipah virus encephalitis in Malaysia. J Infect. 2001;42:403 10.1053/jinf.2000.078211243752Mounts AW, Kaur H, Parashar UD, Ksiazek TG, Cannon D, Arokiasamy JT, ; Nipah Virus Nosocomial Study Group A cohort study of health care workers to assess nosocomial transmissibility of Nipah virus. J Infect Dis. 2001;183:8103 10.1086/31882211181159Yob JM, Field H, Rashdi AM, Morrissy C, van der Heide B, Rota P, Nipah virus infection in bats (order Chiroptera) in peninsular Malaysia. Emerg Infect Dis. 2001;7:4394111384522Chua KB, Koh CL, Hooi PS, Wee KF, Khong JH, Chua BH, Isolation of Nipah virus from Malaysian island flying-foxes. Microbes Infect. 2002;4:14551 10.1016/S1286-4579(01)01522-211880045ICDDRB Nipah encephalitis outbreak over wide area of western Bangladesh. Health Science Bulletin. 2004;2:711ICDDRB Person-to-person transmission of Nipah virus during outbreak in Faridpur District. Health Science Bulletin. 2004;2:59World Health Organization Nipah virus outbreak(s) in Bangladesh, January–April 2004. Wkly Epidemiol Rec. 2004;79:1687115132054Hsu VP, Hossain MJ, Parashar UD, Ali MM, Ksiazek TG, Kuzmin I, Nipah virus encephalitis reemergence, Bangladesh. Emerg Infect Dis. 2004;10:2082715663842AbuBakar S, Chang LY, Ali AR, Sharifah SH, Yusoff K, Zamrod Z Isolation and molecular identification of Nipah virus from pigs. Emerg Infect Dis. 2004;10:22283015663869Rodriguez JJ, Cruz CD, Horvath CM Identification of the nuclear export signal and STAT-binding domains of the Nipah virus V protein reveals mechanisms underlying interferon evasion. J Virol. 2004;78:535867 10.1128/JVI.78.10.5358-5367.200415113915Poch O, Blumberg BM, Bougueleret L, Tordo N Sequence comparison of five polymerases (L proteins) of unsegmented negative-strand RNA viruses: theoretical assignment of functional domains. J Gen Virol. 1990;71:115362 10.1099/0022-1317-71-5-11532161049Reynes JM, Counor D, Ong S, Faure C, Seng V, Molia S, Nipah virus in Lyle's Flying Foxes, Cambodia. Emerg Infect Dis. 2005;11:1042716022778