Emerg Infect DisEmerging Infect. DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention15663869332336104-045210.3201/eid1012.040452DispatchDispatchIsolation and Molecular Identification of Nipah Virus from PigsMolecular Identification of Nipah Virus from PigsAbuBakarSazaly*ChangLi-Yen*AliA.R. MohdSharifahS.H.YusoffKhatijahZamrodZulkeflie§University of Malaya, Kuala Lumpur, Malaysia;Veterinary Research Institute, Ipoh, Perak, Malaysia;Universiti Putra Malaysia, Selangor, Malaysia;Universiti Kebangsaan Malaysia, Selangor, MalaysiaAddress for correspondence: Sazaly AbuBakar, Department of Medical Microbiology, Faculty of Medicine, University of Malaya, 50603, Kuala Lumpur, Malaysia; fax: 603-79675757; email: sazaly@um.edu.my122004101222282230

Nipah viruses from pigs from a Malaysian 1998 outbreak were isolated and sequenced. At least two different Nipah virus strains, including a previously unreported strain, were identified. The findings highlight the possibility that the Malaysia outbreaks had two origins of Nipah virus infections.

Key words: MalaysiaNipah virusPigsdispatch

An outbreak of febrile encephalitis with high death rates (≈56%) occurred among pig farmers in a pig-farming community in Tambun, Perak, Malaysia, in 1998 (1). The disease spread southward within 4 to 6 months to several other pig-farming communities near or around Seremban, a city ≈300 km from Tambun (2). By April 1999, at least 85 deaths were reported in Seremban; 15 were recorded in Tambun (3). All the cases from Tambun were recorded before the outbreak in Seremban. During the outbreaks, pigs in Seremban and Tambun manifested acute respiratory distress syndrome and encephalitis and subsequently died (2,3). Since the incidence of the disease in humans paralleled the occurrence of the disease in pigs, infected pigs were presumed to be the main and perhaps only source of infections in humans (4,5). Culling almost one million suspected infected pigs effectively curtailed the spread of the disease, except for a cluster of infections, the last to be reported, which occurred in the south, in Sungai Buloh, in 1999 (3).

Nipah virus (NV) was eventually isolated from patients manifesting the typical pig-farming–associated fatal encephalitis, confirming the viral source of the infection (4,5). The whole genome sequence of the virus was determined, and its close phylogenetic relationship to Hendra virus (HV) was shown (6,7). The viruses were designated as members of a new genus, Henipavirus, of the Paramyxoviridae family (8,9). Several other human isolates were also sequenced and found to share a high degree of sequence similarity to that of the first isolate (10). More recently, an NV was isolated from flying foxes of Tioman Island, located east of peninsular Malaysia (11). The virus had such high sequence similarity to all the human NV isolates that it was suggested as the potential source of NV (12). It was proposed that pigs in the north contracted the infection through contact with NV of flying foxes (from flying fox urine or leftover fruit) and that humans then acquired the infection by direct handling infected pigs (1214). Movement of infected pigs was responsible for subsequent foci of outbreaks in the south (2,3). Whether the index outbreak in Tambun was a result of a single event (transmission of a flying fox NV to pig), which resulted in clonal propagation and transmission of the virus to pigs and subsequently humans, was unclear. To date, no reports have indicated that NV found in humans shared similar genome sequences to those found in pigs, although NV has been detected and isolated from infected pig samples (15). We present results from analyses of the whole genome sequence of three representative NV isolates from pigs from the three outbreak clusters, one from the north, Tambun, and two from the south, Seremban and Sungai Buloh.

The Study

In this study, NV from pigs (NV/MY/99/VRI-0626, NV/MY/99/VRI-1413, and NV/MY/99/VRI-2794) were isolated by the Veterinary Research Institute, Malaysia. The isolate from a human patient sample from Seremban (NV/MY/99/UM-0128) was isolated at the Department of Medical Microbiology, Faculty of Medicine, University of Malaya, Malaysia. NV of pigs was isolated after infection of Vero (African green monkey kidney) cells with lung tissue samples of pigs that died from the acute respiratory distress syndrome and encephalitis. In Vero cell cultures, all the NV pig isolates manifested similar cytopathologic effects. RNA was extracted from NV-infected cells at passage three, and the genome of the four virus isolates was sequenced in its entirety and analyzed together with the genome sequences of all other NV isolates available in GenBank (Table). The genome sequence of the human NV isolate AF212302, designated as CDC, was used as the reference sequence. The sequences were aligned and manually edited. The phylogenetic trees were constructed and displayed as previously described (16). All new NV genome sequences were deposited in European Molecular Biology Laboratory (accession numbers AJ564621, AJ564622, AJ564623, and AJ627196).

A high sequence similarity (>99%) between NV sequences of pigs and humans was noted after aligning the whole genome sequences. A maximum likelihood phylogenetic tree, constructed by using pig NV genome sequences, showed that NV isolates of pigs clustered tightly together with all known human and flying fox Nipah viruses (Figure A). As expected, all NV sequences clustered with HV to form a distinct group from other established genera within the Paramyxovirinae subfamily. The NV-Tambun showed the most divergent genome sequence from all other known Nipah viruses, in a phylogram constructed by using cumulative nucleotide sequence differences (Figure B). The overall divergence value was nonetheless small (<1%); hence, the isolate remained within the NV cluster.

Phylogenetic trees illustrating the relationships of the pig Nipah virus isolates to all other known Nipah viruses and the related members of subfamily Paramyxovirinae. A) The maximum likelihood tree was drawn by using alignments of the full genome sequences. All the new isolates described in the study are shown in boxes. Abbreviations used and accession numbers not described elsewhere in the text are in parenthesis: Tupaia paramyxovirus (Tupaia-A) (AF079780); (Tupaia-N) (NC_002199), and Hendra virus (HV-A) (AF017149); (HV-N) (NC_001906). B) Unrooted maximum likelihood plot was constructed by using alignments of all the nucleotide differences in the Nipah virus gene coding regions (N, P, M, F, and G) shown in Table A1, Table A2, and Table A3 by artificially treating all the differences as a single stretch of nucleotide sequence.

Subsequent analysis of the deduced amino acid sequences of the NV-Seremban pig isolate showed that they were identical to the human NV isolates CDC and UMMC2 (Table A1, Table A2, Table A3). By contrast, the NV-Sungai Buloh pig isolate–deduced amino acid sequences were identical to the human NV isolates UMMC1 and UM-0128. NV-Seremban differed from NV-Sungai Buloh at only one amino acid position (1645) within the polymerase protein (L). Both isolates, however, differed from the flying fox NV isolate, NV-Flying Fox, at three amino acid positions, residues 30, 206, and 348 in the coding regions of nucleoprotein (N), phosphoprotein (P), and fusion protein (F), respectively. In contrast, the NV-Tambun pig isolate had a distinct signature sequence in comparison to all other NV. NV-Tambun differed from all known NV at 47 nucleotide positions; 28 of these differences occurred within the virus coding regions. The nucleotide differences were translated into amino acid changes at 11 positions; residues 274, 304, and 378 of the P protein, residues 147 and 250 of the matrix protein (M) and F protein, respectively, residues 20 and 272 of the glycoprotein (G), and residues 223, 1645, 1753, and 2039 of the L protein. Amino acid changes noted in the highly phosphorylated P protein at positions 274 and 304 resulted in residue changes from serine to arginine and threonine to alanine, respectively. These changes may reduce the potential phosphorylation sites in the P protein since serine, threonine, tyrosine, and histidine residues are the common targets for protein phosphorylation. Substitution of amino acids at positions 223 (threonine → asparagine), 1645 (serine → phenylalanine) and 2039 (histidine → asparagine) in L protein may also reduce the number of predicted potential phosphorylation sites on L. A substitution of amino acid isoleucine for asparagine at position 20 of the G protein added a potential glycosylation site apart from the eight identified N-linked glycosylation sites (6). However, the additional potential glycosylation site is located at the cytoplasmic tail of the protein and might make the addition of an additional glycan extremely unlikely, if not impossible. No nucleotide and consequently deduced amino acid variations, however, were observed in the N protein.

Conclusions

Findings reported here present for the first time molecular evidence that at least two major strains of NV of pigs were circulating during the 1998 NV outbreak in Malaysia, one strain from the initial outbreak in the north (NV-Tambun) and the other strain from the subsequent outbreak approximately 4 months later in the south (NV-Seremban and NV-Sungai Buloh). The NV-Seremban and NV-Sungai Buloh pig isolates had identical sequences to those reported from human infections, which confirmed that the infections in humans during the southern outbreak originated from infected pigs. No record of isolation of the NV-Tambun is available from patients from the initial outbreak in Tambun or from the subsequent outbreaks. Isolation of NV-Seremban and NV-Sungai Buloh was not reported from the Tambun outbreak. Hence, ascertaining if the two major strains originated from the same initial focus of infection, Tambun, is not possible. Alternatively, the NV-Tambun could be the basal ancestral strain from which the later southern strain evolved. Two findings suggested this hypothesis: the Tambun outbreak occurred at least 4 months before the Seremban outbreak, and the sequence differences between NV-Seremban and NV-Sungai Buloh occurred as a result of genetic drift, a phenomenon not uncommon amongst RNA viruses. Then again, this occurrence is unlikely considering that the genome sequence of NV-Tambun diverges from the NV-Flying Fox, purportedly the initial source of NV infections. In addition, given that the NV sequences of both humans and pigs (UM-0128, NV-Seremban and NV-Sungai Buloh), sequenced independently in different laboratories, were practically identical, the sequence differences were not likely caused by inherent polymerase chain reaction errors or adaptation to tissue culture conditions. Therefore, the NV-Tambun strain is the more likely causal agent for the initial outbreak among pigs in Tambun, resulting from an infection acquired from a yet-to-be-identified source. By contrast, the subsequent outbreaks in the south were due to the pig NV isolates with higher sequence similarities to the NV-Flying Fox of Tioman Island. This finding implied that the 1998 Malaysia NV outbreak is unlikely to be due to a single transmission of NV from flying foxes of Tioman Island to pigs, but it points to the possibility of at least two different origins of NV infections.

Suggested citation for this article: AbuBakar S, Chang LY, Mohd Ali AR, Sharifah SH, Yusoff K, Zamrod Z. Isolation and molecular identification of Nipah virus from pigs. Emerg Infect Dis [serial on the Internet]. 2004 Dec [date cited]. http://dx.doi.org/10.3201/eid1012.040452

This project was supported by the Ministry of Science, Technology and Innovation, Malaysia research grant #01-02-03-004BTK/ER/28. The Department of Veterinary Services, Veterinary Research Institute, Ipoh, Perak, Malaysia, and the Department of Medical Microbiology, Faculty of Medicine, University of Malaya permitted us to use their facilities and provided technical and laboratory assistance.

Comparisons of the genome sequences of Nipah virus isolates<sup>a</sup>
Gene:
N
P
M
F
G
L
NV isolateNucleotide position
896178208289101132212431243949552356774995110171043592643064505438149117153668222625412704365747144749493452575268611564866609
CDCccacagtaatcgcgcttgcgaatggcgccacccaccag
UMMC1................................t.....
UMMC2...............................t......
NV-Flying
Foxtt.............c..ta..----------------
UM-0128...................a..........t.t.....
NV-Tambun..ctgacggctata..aa..ggca.aat.g..tgtaga
NV-Seremban............................t.........
NV-Sungai
Buloh..............t....a............t.....

aResults shown are a summary of all the nucleotides (a) and its corresponding deduced amino acids (b) differences in the NV gene coding regions and nucleotides differences in the NV noncoding regions (c) of the NV-Tambun, NV-Seremban, NV-Sungai Buloh, and UM-0128 in comparison to other known NV isolates.

. Comparisons of the genome sequences of Nipah virus isolates<sup>a</sup>
Gene:
N
P
M
F
G
L
NV isolateAmino acid position
30206274276304378708104147165175189250317339348208810215018127233951223742847902121915721583164517531756203921622203
CDCTPSGTENQSLLKTLIMIQGLLTDQRTLPLPLHSMDHKK
UMMC1................................F.....
UMMC2......................................
NV-Flying
FoxIL.............T......----------------
UM-0128................................F.....
NV-Tambun..R.AK..G...I...N....A...N......FV.N..
NV-Seremban......................................
NV-Sungai
Buloh................................F.....

aResults shown are a summary of all the nucleotides (a) and its corresponding deduced amino acids (b) differences in the NV gene coding regions and nucleotides differences in the NV noncoding regions (c) of the NV-Tambun, NV-Seremban, NV-Sungai Buloh, and UM-0128 in comparison to other known NV isolates.

. Comparisons of the genome sequences of Nipah virus isolates<sup>a</sup>
Genome:
Nucleotide position
NV isolate20129552052911471669178020352097221523814952624164618847886488721112411309113151817318180
CDCctacgagacgcttaatgtagtt
UMMC1......................
UMMC2...................a..
NV-Flying
Foxt.........t.......----
UM-0128......................
NV-Tambun.cgaagtcaa.accgcact.cg
NV-Seremban......................
NV-Sungai
Buloh......................

aResults shown are a summary of all the nucleotides (a) and its corresponding deduced amino acids (b) differences in the NV gene coding regions and nucleotides differences in the NV noncoding regions (c) of the NV-Tambun, NV-Seremban, NV-Sungai Buloh, and UM-0128 in comparison to other known NV isolates.

Nipah virus isolates used in the study
IsolateHostAccession no.
CDCHumanAF212302
UMMC1HumanAY029767
UMMC2HumanAY029768
NV-Flying FoxFlying foxesAF376747
UM-0128 (NV/MY/99/UM-0128)aHumanAJ564623
NV-Tambun (NV/MY/99/VRI-0626)aPigAJ627196
NV-Seremban (NV/MY/99/VRI-1413)aPigAJ564622
NV-Sungai Buloh (NV/MY/99/VRI-2794)aPigAJ564621

aNew isolates described in the study.

Dr. AbuBakar is a virologist with the Faculty of Medicine, University of Malaya. His research interest is in viral diseases caused by dengue virus, human enterovirus 71, and Nipah virus.

ReferencesWong KT, Shieh WJ, Zaki SR, Tan CT Nipah virus infection, an emerging paramyxoviral zoonosis. Springer Semin Immunopathol. 2002;24:21528 10.1007/s00281-002-0106-y12503066Centers for Disease Control and Prevention Outbreak of Hendra-like virus—Malaysia and Singapore, 1998–1999. MMWR Morb Mortal Wkly Rep. 1999;48:265910227800Centers for Disease Control and Prevention Update: outbreak of Nipah virus—Malaysia and Singapore, 1999. MMWR Morb Mortal Wkly Rep. 1999;48:335710366143Chua KB, Goh KJ, Wong KT, Kamarulzaman A, Tan PSK, Ksiazek TG, Fatal encephalitis due to Nipah virus among pig-farmers in Malaysia. Lancet. 1999;354:12579 10.1016/S0140-6736(99)04299-310520635Chua 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.102611504554Mayo MA Virus taxonomy—Houston 2002. Arch Virol. 2002;147:10716 10.1007/s00705020003612021875Mayo MA A summary of taxonomic changes recently approved by ICTV. Arch Virol. 2002;147:16556 10.1007/s00705020003912181683Chan YP, Chua KB, Koh CL, Lim ME, Lam SK Complete nucleotide sequences of Nipah virus isolates from Malaysia. J Gen Virol. 2001;82:2151511514724Chua KB, Koh CL, Hooi PS, Wee KF, Khong JH, Chua BH, Isolation of Nipah virus from Malaysia Island flying-foxes. Microbes Infect. 2002;4:14551 10.1016/S1286-4579(01)01522-211880045Chua KB Nipah virus outbreak in Malaysia. J Clin Virol. 2003;26:26575 10.1016/S1386-6532(02)00268-812637075Chua KB, Chua BH, Wang CW Anthropogenic deforestation, El Niño and the emergence of Nipah virus in Malaysia. Malays J Pathol. 2002;24:152116329551Lam SK Nipah virus—a potential agent of bioterrorism? Antiviral Res. 2003;57:1139 10.1016/S0166-3542(02)00204-812615307Maizan M, Mohd Ali AR, Sharifah SH The identification and distinction between Nipah virus and Hendra virus by using RT-PCR, sequencing and restriction enzyme analysis. Asia Pac J Mol Biol Biotechnol. 2000;8:1016AbuBakar S, Wong PF, Chan YF Emergence of dengue virus type 4 genotype IIA in Malaysia. J Gen Virol. 2002;83:24374212237425