- Short report
- Open Access
Molecular characterization of Nipah virus from Pteropus hypomelanus in Southern Thailand
Virology Journal volume 13, Article number: 53 (2016)
Nipah virus (NiV) first emerged in Malaysia in 1998, with two bat species (Pteropus hypomelanus and P. vampyrus) as the putative natural reservoirs. In 2002, NiV IgG antibodies were detected in these species from Thailand, but viral RNA could not be detected for strain characterization. Two strains of NiV (Malaysia and Bangladesh) have been found in P. lylei in central Thailand, although Bangladesh strain, the causative strain for the outbreak in Bangladesh since 2001, was dominant. To understand the diversity of NiV in Thailand, this study identified NiV strain, using molecular characterizations, from P. hypomelanus in southern Thailand.
Pooled bat urine specimens were collected from plastic sheet underneath bat roosts in April 2010, and then monthly from December 2010 to May 2011 at an island in southern Thailand. Five in 184 specimens were positive for NiV, using duplex nested RT-PCR assay on partial nucleocapsid fragment (357 bp). Whole sequences of nucleocapsid gene from four bats were characterized. All 5 partial fragments and 4 whole nucleocapsid genes formed a monophyletic with NiV-MY.
Our study showed that P. hypomelanus in southern Thailand and from Malaysia, a bordering country, harbored similar NiV. This finding indicates that NiV is not limited to central Thailand or P. lylei species, and it may be a source of inter-species transmission. This indicates a higher potential for a widespread NiV outbreak in Thailand. NiV surveillance in Pteropus bats, the major natural reservoirs, should be conducted continuously in countries or regions with high susceptibility to outbreaks.
Nipah virus (NiV) has been genetically characterized as Malaysia (NiV-MY) or Bangladesh (NiV-BD) strains, termed based on the country of their first outbreak, with fruit bats as the main reservoir . Based on currently available complete N gene sequences, NiV-MY was obtained from Malaysia (Pteropus hypomelanus, P. vampyrus, sick pigs and patients) and Cambodia (P. lylei), while NiV-BD was obtained from patients from Bangladesh and India. Outbreak in Malaysia of NiV-MY in 1998 contributed to encephalitis with 39 % fatality rate . Recent encephalitis outbreak (most likely caused by NiV) in Southern Philippines (2014) originated from fruit bats, transmitted to horses and then to humans, which was attributable to either horse slaughter or horse meat consumption . Secondary man-to-man spread was evident, with acute encephalitis, severe influenza-like illness, or meningitis (82 % fatality). NiV-BD has been associated with the Indian (2001) and Bangladeshi (since 2001) outbreaks of encephalitis and respiratory distress syndrome, with a fatality rate of over 70 % . Infectivity and pathogenicity of the two strains may be different, as NiV-BD seems to be associated with higher incidence of respiratory disease. A study found higher level of oral shedding in ferrets infected with NiV-BD than with NiV-MY . However, both strains caused similar respiratory tract lesions in Syrian hamsters . On the other hand, disease does not develop in Pteropid bats, whether infected naturally or experimentally [7–9].
Serological studies have demonstrated evidences of NiV infection in multiple bat species, including frugivorous and insectivorous bats. However, viral isolation and molecular characterization was mostly only successful in Pteropus species . NiV-MY was isolated from urine samples of P. hypomelanus  and P. vampyrus  from Malaysia. NiV-MY has also been found in other Pteropus bats from various countries; P. lylei in Cambodia, nucleocapsid (N) sequence of 1599 bp, shared 98 % identity with the Malaysian patient (AF212302)  and P. vampyrus in Indonesia, matrix gene of 251 bp had 100 % identity to NiV-MY in the Malaysian patient (AF212302) and 99.6 % homology with the Malaysian P. vampyrus bat (FN869553) . P. giganteus is believed to be the source of NiV infection in India and Bangladesh via contaminated date palm sap. NiV-BD has been detected from a liver homogenate of P. giganteus captured in India with its partial N gene (205 bp), showing 100 % homology with NiV infected patients in India (FJ513078.1) and Bangladesh (AY988601) . On the other hand NiV was detected in P. lylei from Thailand, with sequences (357 bp) sharing 98–99 % nucleotide homology with Bangladeshi patients (AY988601) , but there has been no report of an outbreak in humans or other animals.
The intensive surveillance of NiV in P. lylei in central Thailand has been conducted yearly since 2002. To date, NiV RNA was only detected in P. lylei among two other Pteropus species in Thailand. Fifty two (88.1 %) NiV-BD and 7 (11.9 %) NiV-MY were reported between 2002 and 2008 [15–17]. NiV IgG antibody has been found in P. hypomelanus and P. vampyrus from southern Thailand, but NiV RNA could not be successfully detected . The objective of this study was to identify and characterize the strains of NiV found in P. hypomelanus, in southern Thailand.
Survey location was an island in Southern Thailand, a National Park, where only P. hypomelanus bats roosted in trees along the beach (Fig. 1). Numbers of bat population varied each month due to rain storms (Table 1). Bat urine samples were collected under roost trees with permission from the Department of National Parks, Wildlife and Plant Conservation (No. 0907.1/20713). The samples were collected once in April 2010, and then monthly from December 2010 to May 2011, using a plastic sheet protocol as previously described . Two swabs of pooled bat urine collected from plastic sheet were stored in Lysis buffer (bioMérieux), transported on ice to the laboratory within 48 h and stored at −80 °C until analysis. The nucleic acid was extracted from bat urine using the NucliSENS easyMAG® extraction kit (bioMérieux). NiV RNA was first screened by duplex nested RT-PCR as previously described . Only PCR positive specimens were further characterized for whole N gene by 3 additional PCR assays. All PCR positive specimens were directly sequenced on hemi-nested PCR products using an automated ABI PRISM 377 model sequencer. The sequence segments were assembled by BioEdit program . Phylogenetic trees were generated by using maximum-likelihood method based on 357 bp (Fig. 2a) and 1599 bp (Fig. 2b) of N gene.
In addition, whole N gene sequences from six urine specimens of P. lylei, collected during 2010–2011 from central Thailand (see the map in Fig. 1), were included in the study. The sample collection and amplification methods were the same as for P. hypomelanus.”
Between 24 and 40 pooled bat urine samples were collected each trip (Table 1). NiV RNA was found in 5/184 (2.7 %) samples; in April 2010 (2 of 28, 7.1 %), February 2011 (1 of 24, 4.1 %) and May 2011 (2 of 29, 6.9 %) (Table 1). All NiV sequences (357 bp) shared 98.6–100 % identity with NiV-MY from the Malaysian P. hypomelanus (AF376747) . Four of 5 PCR positive specimens were successfully sequenced further for the complete coding domains of the N gene (1599 bp). They showed 99.1–99.4 % identity to NiV-MY from P. hypomelanus (Table 1). Whole N gene sequencing of sample no. 18254 (KT163257) with 100 % identity to Malaysian P. hypomelanus NiV (357 bp) was not achieved due to low amount of viral RNA.
Phylogenetic analyses using maximum-likelihood methods from 357 and 1599 bp of N gene are shown in Fig. 1a and b respectively. Sequences were analyzed with available NiV sequences from bats, humans and pigs in GenBank, and six additional whole N gene sequences of NiV-BD obtained in 2010–2011 from urine of P. lylei bats in central Thailand. The NiVs from P. hypomelanus in this study form a monophyletic clade with other NiV-MY (Fig. 1a and b). All of NiV-MY (partial N gene, Fig. 2a) from P. lylei clustered with P. hypomelanus NiVs from southern Thailand (this study) and Malaysia.
Further, the similarity between the whole NiV-BD N gene from Thai P. lylei (KT163251- KT163256) and Bangladeshi patient (AY988601) was 99.1–99.2 % (Fig. 2b). The nucleotide identity between NiV-N gene (205 bp) from P. giganteus in India  (JF899339, not included in the phylogenetic tree) and NiVs found in Thailand from P. hypomelanus (no.s 1353, 1753, 8554, 16454) and P. lylei (no.s 1054, 3054, 15053, 15553, 14653, 33654) were 97.1 and 99.5 % respectively.
The visual alignment of the amino acid (aa) sequences was constructed from 533 aa of 1599 nucleotides of N gene open reading frame sequences (Fig. 3). The aa difference of N gene between strains NiV-MY and -BD sequences were found at 5 domains (positions 387, 505, 506, 508 and 521). The N-aa sequences between sick pigs and patients from Malaysia were identical, while one aa sequence difference was observed in NiV isolate from P. hypomelanus in Malaysia (Fig. 3). There were 5–6 aa sequence differences between full-length N in NiV isolates from P. hypomelanus from southern Thailand and P. hypomelanus in Malaysia. However, they had only 1–2 aa sequence differences when compared with NiV isolate from P. vampyrus in Malaysia, and 4–5 aa sequence differences with NiV from P. lylei in Cambodia. Interestingly, N-aa sequences among four NiV-BD isolates from P. lylei in central Thailand (GenBank accession nos. KT163251- KT163254) and from Bangladeshi patient (GenBank accession no. JN808861) were identical (Fig. 3).
Discussion and conclusion
Whether P. hypomelanus or P. vampyrus bats were the spill-over hosts in the Malaysian outbreak has been unclear. Evidence in support of P. hypomelanus as the spill-over host is based on the near identical NiV sequence from this species to the other three NiVs isolated from human patients, and it was isolated shortly after the Malaysian outbreak . In addition, the amino acid sequences obtained from NiV’s whole genome differed from NiV from two infected pigs at merely 3 positions . Nonetheless, NiV sequences from P. vampyrus bats have more differing nucleotides to the isolates from the patients during the outbreak than NiV from P. hypomelanus, 98 and 6 nucleotides, respectively.
Our study at the southern island confirmed the similarity between NiV-MY from P. hypomelanus from Thailand and the Malaysian outbreak (in humans and pigs), and was in accord with the finding from one P. hypomelanus from Tioman island, Malaysia, about 960 km from our surveyed site . From our finding, it remains possible that P. hypomelanus is restricted to NiV-MY, as only this strain was found in the 5 pooled urine specimens positive for NiV. Unfortunately, there is limited data to support host specificity as only one NiV-MY from P. hypomelanus has been reported from Malaysia . A comprehensive analysis of bats across their migratory routes in Indo-Australian region should be conducted to clearly understand the ecology and transmission of NiV in P. hypomelanus.
The co-circulation of NiV-BD and NiV-MY strains within population of P. lylei in central Thailand  remains a mystery. High percent similarity of N nucleotide sequences from P. lylei and P. giganteus, and NiV patients in Bangladesh and the identity of whole N-aa among NiV isolates from P. lylei in Thailand and NiV patients in Bangladesh (Fig. 3), extrapolates the possibility of P. lylei bat in Thailand as a source of NiV infection to humans or other animals. The phylogenetic analysis on 357 bp (Fig. 2a) revealed that NiV-MY in P. lylei from Thailand, and P. hypomelanus from Malaysia and Thailand clustered together. This finding introduces the possibility of NiV transmission between P. lylei and P. hypomelanus, in addition to the previously suggested transmission between P. vampyrus and P. lylei . Furthermore, the higher percent similarity of NiV N-aa between P. hypomelanus in Thailand and P. vampyrus in Malaysia than P. hypomelanus in Malaysia, raise the possibility of transmission between P. hypomelanus and P. vampyrus. Molecular characterization of NiV N and other genes, particularly those encoding the envelop glycoprotein from P. vampyrus in Thailand is required to investigate this possibility.
The alarming issue is that many bat species can harbor NiV. NiV-BD has been found in P. giganteus (India) , and P. lylei (Thailand) . NiV-MY has been found in several bats species in many countries; P. hypomelanus (Malaysia  and Thailand [this study]), P. vampyrus (Malaysia  and Indonesia ), P. lylei (Cambodia  and Thailand ), Hipposideros larvatus (Thailand) , Taphozous species (Thailand)  and Rousettus amplexicaudatus (East Timor) . The wide susceptibility of NiV to multiple host species allows for expedition of a larger-scale outbreak.
In conclusion, detecting NiV-MY strain in P. hypomelanus in southern Thailand, whose sequences share similarities to that identified in the Malaysian outbreak, and to P. lylei bats in central Thailand, suggests a potential role of P. hypomelanus as a possible source for inter-species transmission. This may result in the emergence of NiV infection to other animals and humans in both southern and central Thailand. Continuing surveillance of NiV infection in bats, humans, and pigs is vital for early detection and can potentially reduce the scale of an outbreak.
Ethics approval and consent to participate
Bat urine samples were collected with permission from the Department of National Parks, Wildlife and Plant Conservation (No. 0907.1/20713). The urine specimens were collected under roost trees, no bat was caught for sample collection.
Consent for publication
Availability of data and materials
The GenBank accession numbers for Nipah virus sequences reported in this paper are: KT163247- KT163250 for whole N gene and KT163257 for 357 bp partial N gene from Pteropus hypomelanus and KT163251- KT163256 for whole N gene from P. lylei.
- N fragment/gene:
NiV Bangladesh strain [NiV strain from the Bangladeshi and Indian outbreaks]
NiV Malaysia strain [NiV strain from the Malaysian outbreak]
Polymerase Chain Reaction
Reverse Transcription - Polymerase Chain Reaction
Rota PA, Lo MK. Molecular virology of the henipaviruses. Curr Top Microbiol Immunol. 2012;359:41–58.
Chua KB. Nipah virus outbreak in Malaysia. J Clin Virol. 2003;26:265–75.
Ching PKG, de los Reyes VC, Sucaldito MN, Tayag E, Columna-Vingno AB, Malbas Jr FF, et al. Outbreak of Henipavirus Infection, Philippines, 2014. Emerg Infect Dis. 2015;21:328–31.
Luby SP, Hossain MJ, Gurley ES, Ahmed BN, Banu S, Khan SU, et al. Recurrent zoonotic transmission of Nipah virus into humans, Bangladesh, 2001–2007. Emerg Infect Dis. 2009;15:1229–35.
Clayton BA, Middleton D, Bergfeld J, Haining J, Arkinstall R, Wang L, et al. Transmission routes for Nipah virus from Malaysia and Bangladesh. Emerg Infect Dis. 2012;18:1983–93.
Baseler L, de Wit E, Scott DP, Munster VJ, Feldmann H. Syrian hamsters (Mesocricetus auratus) oronasally inoculated with a Nipah virus isolate from Bangladesh or Malaysia develop similar respiratory tract lesions. Vet Pathol. 2015;52:38–45.
Chua KB, Koh CL, Hooi PS, Wee KF, Khong JH, Chua BH, et al. Isolation of Nipah virus from Malaysian Island flying-foxes. Microb Infect. 2002;4:145–51.
Halpin K, Hyatt AD, Fogarty R, Middleton D, Bingham J, Epstein JH, et al. Pteropid bats are confirmed as the reservoir hosts of henipaviruses: a comprehensive experimental study of virus transmission. Am J Trop Med Hyg. 2011;85:946–51.
Rahman SA, Hassan L, Epstein JH, Mamat ZC, Yatim AM, Hassan SS, et al. Risk factors for Nipah virus infection among pteropid bats, Peninsular Malaysia. Emerg Infect Dis. 2013;19:51–60.
Simons RR, Gale P, Horigan V, Snary EL, Breed AC. Potential for Introduction of Bat-Borne Zoonotic Viruses into the EU: A Review. Viruses. 2014;6:2084–121.
Rahman SA, Hassan SS, Olival KJ, Mohamed M, Chang LY, Hassan L, et al. Characterization of Nipah virus from naturally infected Pteropus vampyrus bats, Malaysia. Emerg Infect Dis. 2010;16:1990–3.
Reynes JM, Counor D, Ong S, Faure C, Seng V, Molia S, et al. Nipah virus in Lyle’s flying foxes, Cambodia. Emerg Infect Dis. 2005;11:1042–7.
Sendow I, Ratnawati A, Taylor T, Adjid RM, Saepulloh M, Barr J, et al. Nipah virus in the fruit bat Pteropus vampyrus in Sumatera, Indonesia. PLoS One. 2013;8:e69544. doi:10.1371/journal.pone.0069544.
Yadav PD, Raut CG, Shete AM, Mishra AC, Towner JS, Nichol ST, et al. Detection of Nipah virus RNA in fruit bat (Pteropus giganteus) from India. Am J Trop Med Hyg. 2012;87:576–8.
Wacharapluesadee S, Ngamprasertwong T, Kaewpom T, Kattong P, Rodpan A, Wanghongsa S, et al. Genetic characterization of Nipah virus from Thai fruit bats (Pteropus lylei). Asian Biomedicine. 2013;7:813–9.
Wacharapluesadee S, Lumlertdacha B, Boongird K, Wanghongsa S, Chanhome L, Rollin P, et al. Bat Nipah virus, Thailand. Emerg Infect Dis. 2005;11:1949–51.
Wacharapluesadee S, Hemachudha T. Duplex nested RT-PCR for detection of Nipah virus RNA from urine specimens of bats. J Virol Methods. 2007;141:97–101.
Chua KB. A novel approach for collecting samples from fruit bats for isolation of infectious agents. Microb Infect. 2003;5:487–90.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.
AbuBakar S, Chang LY, Ali AR, Sharifah SH, Yusoff K, Zamrod Z. Isolation and molecular identification of Nipah virus from pigs. Emerg Infect Dis. 2004;12:2228–30.
Breed AC, Meers J, Sendow I, Bossart KN, Barr JA, Smith I, et al. The distribution of henipaviruses in Southeast Asia and Australasia: Is Wallace’s line a barrier to Nipah virus. PLoS One. 2013;8:e61316. doi:10.1371/journal.pone.0061316.
Wacharapluesadee S, Boongird K, Wanghongsa S, Ratanasetyuth N, Supavonwong P, Saengsen D, et al. A longitudinal study of the prevalence of Nipah virus in Pteropus lylei bats in Thailand: evidence for seasonal preference in disease transmission. Vector Borne Zoonotic Dis. 2010;10:183–90.
We gratefully acknowledge the local support from the Thai Red Cross Society, Chulalongkorn University, Kasetsart University and the Department of National Parks Wildlife and Plant Conservation. We would also like to thank Ms. Siriporn Ghai for her intensive editing of the manuscript.
This study was supported by research grants from the Higher Education Research Promotion and National Research University Project of Thailand, Office of the Higher Education Commission (HR1160A), the Ratchadaphiseksomphot Endowment Fund of Chulalongkorn University (RES560530148-HR), Health and Biomedical Science Research Program by the National Research Council of Thailand (NRCT) and Health System Research Institute (HSRI), the Research Chair Grant and the Cluster and Program Management Office (CPMO), from the National Science and Technology Development Agency (NSTDA), Thailand.
The authors declare that they have no competing interests.
SW and TH participated in the design of the study and drafted the manuscript. SW, PS, MP, PM, and BK conducted sampling of bat urine samples. TK carried out the molecular genetic studies. AR participated in the sequence alignment. All authors read and approved the final manuscript.