# High Pathogenicity of Nipah Virus from Pteropus lylei Fruit Bats, Cambodia

Maria Gaudino, Noémie Aurine, Claire Dumont, Julien Fouret, Marion Ferren, Cyrille Mathieu, Olivier Reynard, Viktor E. Volchkov, Catherine Legras-Lachuer, Marie Claude Georges-Courbot, et al.

<figure><img src="image_1.png" /></figure>

► To cite this version:

Maria Gaudino, Noémie Aurine, Claire Dumont, Julien Fouret, Marion Ferren, et al. High Pathogenicity of Nipah Virus from Pteropus lylei Fruit Bats, Cambodia. Emerging Infectious Diseases, 2020, 26 (1), pp.104 - 113. 10.3201/eid2601.191284. hal-03091393

HAL Id: hal-03091393

https://cnrs.hal.science/hal-03091393

Submitted on 31 Dec 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# High Pathogenicity of Nipah Virus from *Pteropus lylei* Fruit Bats, Cambodia

Maria Gaudino, Noémie Aurine, Claire Dumont, Julien Fouret, Marion Ferren

Cyrille Mathieu, Olivier Reynard, Viktor E. Volchkov, Catherine Legras-Lachuer, Marie-Claude Georges-Courbot, Branka Horvat

We conducted an in-depth characterization of the Nipah virus (NiV) isolate previously obtained from a *Pteropus lylei* bat in Cambodia in 2003 (CSUR381). We performed full-genome sequencing and phylogenetic analyses and confirmed CSUR381 is part of the NiV-Malaysia genotype. In vitro studies revealed similar cell permissiveness and replication of CSUR381 (compared with 2 other NiV isolates) in both bat and human cell lines. Sequence alignments indicated conservation of the ephrin-B2 and ephrin-B3 receptor binding sites, the glycosylation site on the G attachment protein, as well as the editing site in phosphoprotein, suggesting production of nonstructural proteins V and W, known to counteract the host innate immunity. In the hamster animal model, CSUR381 induced lethal infections. Altogether, these data suggest that the Cambodia bat-derived NiV isolate has high pathogenic potential and, thus, provide insight for further studies and better risk assessment for future NiV outbreaks in Southeast Asia.

Nipah virus (NiV) is a zoonotic paramyxovirus that was first identified as the cause of an outbreak of encephalitis in humans in Malaysia and Singapore during 1998–1999 (1). Although NiV infection remains rare in humans, this virus has captured the attention of the public health community and scientists because of its high case-fatality rate, ranging from 40% in Malaysia to >90% in Bangladesh and

India (2,3). Its high pathogenicity, potential for interspecies transmission, and lack of validated medical countermeasures led to the classification of NiV as a Biosafety Level 4 (BSL-4) pathogen. In 2015, the World Health Organization (WHO) listed NiV as a priority pathogen because of its probability of causing severe outbreaks and subsequently placed NiV on the WHO Blueprint list of priority diseases (4). This designation was strengthened by the NiV outbreak in Kerala, India, where the virus had not previously been reported (3).

NiV is a member of the *Henipavirus* genus, together with Hendra virus, which first emerged in Brisbane, Queensland, Australia, in 1994 (5), and the nonpathogenic Cedar virus, which was discovered in Australia in 2009 (6). In addition, full-length henipalike viral sequences were found in fruit bats in Africa (7) and in rats in China (Moijang virus) (8). As part of the Mononegavirales, NiV has a nonsegmented, negative-sense, single-stranded RNA genome. *Pteropus* fruit bats, commonly known as flying foxes, are considered the henipavirus natural reservoir; these bats, when infected with NiV, do not seem to display any apparent clinical signs of disease (9).

Only 2 NiV lineages are known to circulate in Asia and cause disease in humans: NiV-Malaysia and NiV-Bangladesh. In India and Bangladesh, NiV transmission from bats to humans was shown to occur through the consumption of raw date palm juice or fruits contaminated with the saliva or urine from fruit bats (10). In addition, as observed during the outbreaks in Malaysia in 1998 and 1999, transmission can occur via contact with infected domestic animals, such as pigs, that act as amplifying hosts of the virus (11). Interhuman transmission has been reported in Bangladesh and India in >50% of NiV outbreaks (12).

Clinical manifestations of NiV infection in humans can range from asymptomatic to acute

Author affiliations: Centre International de Recherche en Infectiologie, CIRI, INSERM U1111, CNRS, UMR5308, Univ Lyon, University Claude Bernard Lyon 1, École Normale Supérieure de Lyon, Lyon, France (M. Gaudino, N. Aurine, C. Dumont, J. Fouret, M. Ferren, C. Mathieu, O. Reynard, V.E. Volchkov, M.-C. Georges-Courbot, B. Horvat); ViroScan 3D, Trévoux, France (J. Fouret, C. Legras-Lachuer); University Claude Bernard Lyon 1, LEM, UMR5557, CNRS, INRA, VetAgro Sup, Lyon (C. Legras-Lachuer); Unité de Biologie des Infections Virales Emergentes, Institut Pasteur, INSERM P4, Jean Mérieux, Lyon (M.-C. Georges-Courbot)

DOI: https://doi.org/10.3201/eid2601.191284
respiratory syndrome, generalized vasculitis, and fatal encephalitis. Among the few survivors of NiV outbreaks, long-term neurologic problems have been reported; 20% of patients have residual neurologic sequelae (13), and NiV-Malaysia-infected patients experienced relapse and late-onset encephalitis (14).

NiV and henipa-like viruses have been detected molecularly or serologically in *Pteropus* bats in different countries of Asia (15) and Africa (7), Australia (16), and Brazil (17), and the worldwide distribution of fruit bats poses a continuous threat to another spillover with possible pandemic potential (18). However, since 1998, all NiV cases in humans have been identified in Malaysia, India, Bangladesh, and the Philippines (19). Human cases of NiV have not been reported in Cambodia or neighboring countries since the first serologic detection of NiV in Cambodia and isolation of CSUR381 in *Pteropus lylei* bats in Cambodia in 2003 (20,21). Initial phylogenetic analyses of the nucleoprotein and attachment glycoprotein of CSUR381 suggested the virus was part of the NiV-Malaysia genotype (21). However, a full-genome characterization and phylogenetic analysis have not been performed. In addition, the growth dynamics and virulence of this virus have not been analyzed, thus limiting more comprehensive evaluation of this virus's pathogenic potential. In this study, we performed an in-depth characterization of CSUR381, including its pathogenicity both in vitro and in vivo, ultimately to assess the outbreak risk that isolates circulating in Cambodia pose in Southeast Asia.

## Materials and Methods

### Viruses

In this study, we used 3 different NiV isolates: the NiV isolate CSUR381 from Cambodia (GenBank accession no. MK801755), NiV-Malaysia isolate UMMC1 (GenBank accession no. AY029767), and NiV-Bangladesh isolate SPB200401066 (GenBank accession no. AY988601). CSUR381 was isolated from *P. lylei* bat urine at the Institut Pasteur in Battambang, Cambodia, in 2003 (21), and the other 2 isolates were obtained from infected patients. We produced and titrated all viruses on Vero E6 cells.

### Full-Genome Sequencing

We amplified and titrated the Cambodia NiV isolate on Vero cells and, after the second cell passage, extracted viral RNA from supernatant using the QIAamp Viral RNA Mini Kit (QIAGEN, https://www.qiagen.com) according to the manufacturer's instructions. We treated samples with DNase, purified and quantified

RNA using the QuantiFluor RNA System (Promega, https://www.promega.com), and analyzed using the AATI High Sensitivity Genomic DNA Analysis Fragment Analyzer (Advanced Analytical Technologies Inc., https://www.agilent.com). Then, we amplified viral RNA using the Single Primer Isothermal Amplification Kit (NuGEN, https://www.nugen.com). We prepared a library using Ovation Ultra Low (NuGen), which gave us average DNA fragment sizes of 382-426 bp. We then sequenced the whole genome of CSUR381 using MiSeq Nano v2 (Illumina, https://www.illumina.com), which produced read lengths of 2 × 150 nt.

We carried out genome assembly de novo using SPAdes (http://cab.spbu.ru/software/spades) and predicted open reading frames with Prodigal (https://omictools.com/prodigal-tool). We reconstructed the missing 5' and 3' extremities using the 5' RACE System for Rapid Amplification of cDNA Ends (Invitrogen, https://www.thermofisher.com) according to the manufacturer's recommendations. We performed reverse transcription with SuperScript II Reverse Transcriptase (Invitrogen) using the primer GSP1-leader (5'-GACCATTGATCCAACATC-3') to recover the viral leader sequence and GSP-trailer (5'-AAAGTGATTGTCTACTCACT-3') to recover the trailer sequence. After column purification, we tailed cDNA sequences with cytidine triphosphate and terminal deoxynucleotidyl transferase. Last, we amplified dC-tailed cDNA using the Abridged Anchor Primer provided in the 5' RACE System for Rapid Amplification of cDNA Ends Kit and primers nested-GSP2-leader (5'-TACAGCTTCAATGTCTGGGTCATT-3') to amplify the viral leader sequence, and nested-GSP2-trailer (5'-CAAGTTCAAGGACACCAAAAGT-3') to amplify the viral trailer sequence. We sequenced PCR products using Sanger technology and submitted the complete genome sequence of CSUR381 to GenBank (accession no. MK801755).

### Phylogenetic Analyses

Using the ClustalW algorithm (http://www.clustal.org), we performed multiple alignments for complete genomes and individual gene sequences. We implemented and manually checked the quality of alignments using BioEdit version 7.2.6 (22) and conducted genomic characterization and evolutionary analyses in MEGA version 7.0.26 (23). After determining the best DNA model to use for each alignment, we constructed maximum-likelihood phylogenetic trees for complete NiV genomes and all virus coding sequences. For statistical support, we used 500 bootstrap replicates for the analysis of the complete genome and 1,000 replicates for analyses of each gene.
## Cell Lines and Infection

We cultured NCI-H358 (human bronchioalveolar carcinoma) and Vero E6 (African green monkey kidney) cells in Dulbecco's modified Eagle medium (DMEM) with GlutaMAX (Thermo Fisher Scientific, https://www.thermofisher.com) supplemented with 1% penicillin-streptomycin (10,000 U/mL), 1% L-glutamine, and 10% heat-inactivated (56°C for 30 minutes) fetal calf serum (FCS). We cultured human pulmonary microvascular endothelial cells (HPMECs) (24) in endothelial cell growth medium (Growth Medium MV 2 Kit; PromoCell, https://www.promocell.com). We incubated all these cell lines at 37°C with 5% carbon dioxide; all cell lines, including the *Pteropus* cell line described in the next paragraph, tested negative for *Mycoplasma* spp. by the MycoAlert kit (Lonza, https://www.lonza.com).

We generated a *Pteropus* flying fox cell line using a skin biopsy from the wing membrane of a female *P. giganteus* (also known as *P. medius* and flying fox) bat (25) of the order Yinpterochiroptera. Biopsies were collected from bats by Tiergarten Schönbrunn (Vienna, Austria) staff during regular veterinary checkups following appropriate guidelines to minimize animal stress. The biopsies were washed with sterile phosphate-buffered saline and transferred into Freezing Medium Cryo-SFM (PromoCell), and sample vials were put on dry ice for shipment to Centre International de Recherche en Infectiologie in Lyon, France. To obtain primary cell cultures, we fractionated the biopsies into petri dishes, harvested the homogenates, and incubated them at 37°C with 5% carbon dioxide in DMEM/F-12 (Gibco, https://www.thermofisher.com) supplemented with 10% fetal calf serum FCS, 1% L-glutamine (200 mM), 1,000 U/mL of penicillin, 1,000 U/mL of streptomycin, and 2.50 µg/mL amphotericin B (Gibco). We subsequently immortalized primary cells using the lentiviral vector SV40 large T-antigen produced at Genetic Analysis and Vectorology Platform (AniRA, École Normale Supérieure de Lyon, Lyon). We evaluated different clones on the basis of their morphologic stability and transfectability using jetPRIME kit (Polyplus, https://www.polyplus-transfection.com). We confirmed immortalization of clones by detecting large T-antigen inserts by reverse transcription PCR (RT-PCR). We cultured the final *Pteropus* cell line, which we designated PATGV1.12, in DMEM GlutaMAX supplemented with 10% heat-inactivated FCS. We additionally confirmed that this cell line was derived from *P. giganteus* bats by sequencing the mitochondrial region D-loop (26) and nuclear introns ACOX2, COPS7A, BGN, ROGD1, and STAT5A, which has

been suggested to be pertinent for distinguishing among closely related bat species (27).

We infected cells in 12-well plates at 80% confluence with a multiplicity of infection (MOI) of 0.3. For virus replication kinetics studies, we took 4 time points postinfection into consideration: 0 h, 24 h, 48 h, and 72 h. We performed infections in BLS-4 facility Jean Mérieux (Lyon). For each time point, we collected cell lysates according to validated BSL-4 procedures. We collected supernatants and kept them at -80°C until titration by plaque assay on Vero E6 cells.

## Pseudotyping of Vesicular Stomatitis Virus and Evaluation of Cell Permissiveness

We used rVSVΔG-RFP (a recombinant vesicular stomatitis virus [VSV] in which the envelope glycoprotein G gene is replaced with the red fluorescent protein gene) (28,29) to generate pseudotyped VSVs harboring different combinations of NiV envelope glycoprotein G (attachment protein) and F (fusion protein) on their surfaces. Complementing rVSVΔG-RFP-infected cells with NiV glycoproteins expressed in trans, we were able to produce stocks of pseudotyped VSVs identical in their genetic background and differing only in the nature of their surface glycoproteins. Because the infectivity of rVSVΔG-RFP pseudotypes is restricted to a single round of replication, this tool is largely used for studying viral entry for a broad range of highly pathogenic viruses (30).

To create the pseudotypes, we cloned the NiV glycoprotein G and F genes from RNA isolated from CSUR381, UMMC1, and SPB200401066 into 6 separate pCAGGS plasmid vectors. We transfected these 3 plasmid pairs separately into BSR-T7 cells using TransIT-LT1 Transfection Reagent (Mirus Bio, https://www.mirusbio.com). We infected cells with rVSVΔG-RFP 16 h after transfection to produce a pseudotyped VSV for each NiV isolate. We collected supernatants at 24 h postinfection and concentrated pseudotyped VSVs by ultracentrifugation (28,000 rpm for 2 h at 4°C). We titrated these viruses on Vero cells. To evaluate viral entry into different cell lines, we performed infections in 24-well plates using 80% confluent, adherent cells and a 1-h contact between virus and cells. We determined the percentage of cells infected 6 h postinfection by quantifying cells expressing RFP via flow cytometry on a BD LSRFortessa (https://www.bd.com).

## RNA Extraction and Real-Time RT-PCR

At the indicated time points, we collected cells and extracted RNA using the NucleoSpin RNA Kit (Macherey-Nagel, https://www.mn-net.com) according
to the manufacturer's instructions. We assessed the yield and purity of extracted RNA using the DS-11-FX spectrophotometer (DeNovix, https://www.denovix.com). We reverse transcribed extracted RNA using the iScript Select cDNA Synthesis Kit (Bio-Rad, https://www.bio-rad.com) and performed real-time PCR using Platinum SYBR Green qPCR SuperMix-UDG (Invitrogen) on a StepOnePlus Real-Time PCR System (Applied Biosystems, https://www.thermofisher.com). As previously described (31), we amplified the NiV nucleoprotein gene and the *Pteropus* glyceraldehyde 3-phosphate dehydrogenase housekeeping gene using forward primer 5'-ATCATCCCTGCTTCTACT-3' and reverse primer 3'AGGTCAGATCCACAACT-5'. We analyzed quantitative RT-PCR results using StepOne version 2.3 (Applied Biosystems).

## Experimental Infection of Hamsters

We obtained 2-month-old male golden hamsters (*Mesocricetus auratus*) from Janvier Labs (https://www.janvier-labs.com). We housed hamsters in a BSL-4 containment facility (INSERM P4, Jean Mérieux, Lyon) and handled them according to the regulations for animal maintenance of France. We treated hamsters with isoflurane anesthesia before manipulations. We subcutaneously infected 2 groups of 6 hamsters with a high dose (13,500 PFU/animal) of either the NiV-Malaysia or Cambodia NiV isolate and followed hamsters daily to record their body temperature and weight. The regional ethics committee for animal experimentation (Lyon) approved these animal experiments.

## Results

### Full-Genome Characterization and Phylogenetic Analyses

Analysis of the assembled viral sequence of CSUR381 showed a total genome length of 18,246 nt, similar to the lengths of NiV isolates reported in Malaysia. The nucleotide composition was 27.8% T or U, 18.4% C, 33.6% A, and 20.2% G; total GC content was 38.6%. To

<table><caption>**Table 1.** Whole-genome pairwise nucleotide identity comparisons between Nipah virus CSUR381, Cambodia, 2003, and other available henipaviruses</caption><thead><tr><th>Henipavirus (GenBank accession no.)</th><th>Nucleotide identity, %*</th></tr></thead><tbody><tr><td>Nipah/Malaysia/2000/human (NC002728.1)</td><td>97.7</td></tr><tr><td>Nipah/Malaysia/2001/human (AY029768.1)</td><td>97.7</td></tr><tr><td>Nipah/Malaysia/1999/swine (AJ627196.1)</td><td>97.6</td></tr><tr><td>Nipah/Bangladesh/2004/human (AY988601.1)</td><td>91.6</td></tr><tr><td>Nipah/India/2007/human (FJ513078.1)</td><td>91.4</td></tr><tr><td>Nipah/Bangladesh/2008/human (JN808863.1)</td><td>91.4</td></tr><tr><td>Nipah/India/2018/human (MH396625.1)</td><td>91.3</td></tr><tr><td>Hendra/Australia/2008/human (JN255805.1)</td><td>70.0</td></tr><tr><td>Hendra/Australia/2009/bat (JN255803.1)</td><td>69.9</td></tr><tr><td>Hendra/Australia/2007/horse (HM044321.1)</td><td>69.9</td></tr><tr><td>Cedar/Australia/2009/Ptalecto (JQ001776.1)</td><td>56.7</td></tr><tr><td>Paramyxo/Ghana/2009/Eidolonhelvum (HQ660129.1)</td><td>52.9</td></tr><tr><td>Mojiang/China/2012/Rattusflavipectus (NC025352.1)</td><td>48.9</td></tr></tbody></table>

*Calculated by using the p-distance method.

investigate genetic relationships between CSUR381 and other henipaviruses, we constructed distance matrices for the complete genome and for each gene using the p-distance method. When we compared the sequence of CSUR381 with those of other NiVs available in GenBank, the most similar sequences (with 97.7% nucleotide identity) were from NiV-Malaysia human isolates (GenBank accession nos. NC002728.1 and AY029768.1; Table 1). We also calculated nucleotide identity and amino acid homology for each of the 6 structural genes (Table 2). Genetic pairwise comparisons with other NiV isolates showed the lowest nucleotide identity and amino acid homology for phosphoprotein (087.1/82.7%) and the highest for matrix protein (98.9/99.4%).

Using the maximum-likelihood method, we constructed phylogenetic trees on the basis of the complete genome (Figure 1, panel A) and the nucleocapsid gene (Figure 1, panel B). The general time-reversible model for the complete genome and the Kimura 2-parameter model for the nucleocapsid gene were predicted to be the best for performing those particular phylogenetic analyses. CSUR381 clustered with the monophyletic group of the NiV-Malaysia genotype

<table><caption>**Table 2.** Pairwise comparison of NiV CSUR381, Cambodia, 2003, and other available NiVs, by NiV gene*</caption><thead><tr><th rowspan="2">NiV (GenBank accession no.)</th><th colspan="6">NiV gene, % nucleotide identity/% homology of deduced amino acid†</th></tr><tr><th>N</th><th>P</th><th>M</th><th>F</th><th>G</th><th>L</th></tr></thead><tbody><tr><td>Nipah/Malaysia/2010/Pvampyrus (FN869553.1)</td><td>98.3/99.2</td><td>96.3/94.9</td><td>98.9/99.4</td><td>98.4/98.7</td><td>97.1/98.5</td><td>98.0/99.3</td></tr><tr><td>Nipah/Malaysia/2000/human (NC002728.1)</td><td>97.9/98.7</td><td>95.9/94.1</td><td>98.5/99.1</td><td>98.1/98.9</td><td>97.1/98.5</td><td>98.1/99.4</td></tr><tr><td>Nipah/Bangladesh/2004/human (AY988601.1)</td><td>93.8/98.5</td><td>87.8/84.8</td><td>93.0/99.1</td><td>93.0/98.2</td><td>88.3/95.5</td><td>91.7/98.2</td></tr><tr><td>Nipah/Bangladesh/2008/human (JN808863.1)</td><td>93.9/98.7</td><td>87.6/84.4</td><td>93.0/99.1</td><td>93.3/98.5</td><td>88.1/95.5</td><td>91.8/98.4</td></tr><tr><td>Nipah/India/2007/human (FJ513078.1)</td><td>93.5/98.5</td><td>87.4/84.3</td><td>92.9/98.9</td><td>93.1/98.4</td><td>88.0/95.3</td><td>91.9/98.4</td></tr><tr><td>Nipah/India/2018/human (MH396625.1)</td><td>93.3/98.5</td><td>87.1/82.7</td><td>92.6/99.1</td><td>93.0/98.5</td><td>87.2/95.3</td><td>91.7/98.5</td></tr><tr><td>Nipah/Thailand/2010/Plylei (KT163252.1)</td><td>93.7/98.7</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Nipah/Thailand/2010/Phypomelaneus (KT163247.1)</td><td>97.8/99.0</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

*F, fusion protein; G, attachment protein; L, polymerase; M, matrix protein; N, nucleoprotein; NiV, Nipah virus; P, phosphoprotein.

†Both calculated by using the p-distance method.
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Maximum-likelihood phylogenetic analysis of NiV CSUR381, Cambodia, 2003 (red triangle), compared with other henipaviruses and NiVs. A) Phylogenetic tree constructed with complete genome sequences. A general time-reversible model was calculated as the best DNA model to conduct this analysis. B) Phylogenetic tree constructed by using the nucleocapsid gene. The Kimura 2-parameter model was calculated as the best DNA model to conduct this analysis. Bootstrap statistical support is marked on branch nodes. GenBank accession numbers of isolates are provided in branches, and NiV lineages of isolates are indicated. Phylogenetic trees are drawn to scale; scale bars represent branch lengths measured in the number of substitutions per site. NiV, Nipah virus.</figcaption></figure>

for both the whole genome and nucleocapsid gene; bootstrap support was >98% in all cases, confirming the previous partial genomic characterization of CSUR381 (24). We then generated phylogenetic trees for each of the coding sequences of the 6 NiV structural proteins, which gave equivalent results (Appendix Figure 1, https://wwwnc.cdc.gov/EID/article/26/1/19-1284-App1.pdf).

Multiple alignment of the henipavirus phosphoprotein gene (Figure 2) revealed high conservation of

the editing site (5'-AAAAAGGG-3') in CSUR381, similar to other NiV and Hendra virus isolates and different from Cedar virus, a nonpathogenic virus isolated from a *P. alecto* bat in Australia (6). This finding suggests that CSUR381 might produce the nonstructural proteins V and W, capable of interacting with the host innate cellular immune response (32). Comparisons of the deduced V, W, and C amino acid homologies between CSUR381 and other known NiVs showed a variation of 88%-100% (Appendix Table).
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Multiple alignment of the phosphoprotein gene of Nipah virus CSUR381, Cambodia, 2003, and other henipavirus isolates. The highly conserved editing site (5'-AAAAAGGG-3', red outline) is present in all Nipah and Hendra virus sequences but absent in the nonpathogenic Cedar virus sequence. GenBank accession numbers are provided for all isolates.</figcaption></figure>

## Evaluation of Virus Entry

We next determined the cellular permissiveness of a human endothelial cell line (HPMEC), a human respiratory epithelial cell line (NCI-H358), the newly generated *Pteropus* bat cell line (PAT-GV1.12), and Vero cells to CSUR381 compared with the NiV-Malaysia (UMMC1) and NiV-Bangladesh (SPB200401066) isolate using pseudotyped rVSVΔG-RFP viruses. Cell lines were infected for 1 h at an MOI of 0.3. The percentages of cells infected were analyzed by flow cytometry 6 h after infection (Figure 3), and results from HPMEC, NCI-H358, and PAT-GV1.12 were normalized to the findings from Vero. All tested cell lines were permissive to infection with all 3 viruses tested. Entry of NiV pseudotypes into the bat cell line PATGV1.12 and human respiratory epithelial cell line was similar. Compared with the NiV-Malaysia and NiV-Bangladesh pseudotypes, the CSUR381 pseudotyped virus showed higher but not significantly increased entry into the 3 tested cell lines (1-way analysis of variance).

We further analyzed the amino acid sequences of the F and G proteins of the 3 viruses by multiple alignment. The glycosylation site (N529/Q530/T531) (33) and ephrin-B2 and ephrin-B3 binding sites (34) in the G attachment protein were preserved (Appendix Figure 2). In addition, multiple alignments showed that the F cleavage site was preserved among all analyzed NiV isolates (Appendix Figure 3). Last, an analysis of the predicted N-terminal and C-terminal heptad-repeat regions within the F protein, which are needed for NiV fusion (35), showed high conservation, and compared with NiV-Malaysia and NiV-Bangladesh, only 1 aa difference (V159→I) was detected in CSUR381 (Appendix Figure 3). Altogether,

the high conservation of the NiV glycoproteins and results from pseudotype virus studies suggest that CSUR381 can enter target cells at least as well as NiV-Malaysia and NiV-Bangladesh.

## Replication of NiV Isolates in Different Cell Types

To further evaluate the virulence of CSUR381, we compared the replication kinetics of this virus with those of the NiV-Malaysia and NiV-Bangladesh isolates.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Evaluation of entry of VSVΔG-RFPs (vesicular stomatitis virus in which the envelope glycoprotein G gene is replaced with the red fluorescent protein gene) pseudotyped with the surface glycoproteins of NiVs CSUR381 (Cambodia 2003 isolate), UMMC1 (NiV-Malaysia isolate), and SPB200401066 (NiV-Bangladesh isolate) in different cell types. Infections of HPMEC, NCI-H358 (human bronchioalveolar cells), PATGV1.12 (bat cells), and Vero cells were performed at a multiplicity of infection of 0.3, and the percentages of infected cells were evaluated 6 hours postinfection by measuring RFP by flow cytometry and normalizing values to those from Vero cells. Histograms indicate the mean of 3 independent experiments, and error bars indicate upper half of SD. HPMEC, human pulmonary microvascular endothelial cell; NiV, Nipah virus.</figcaption></figure>
We infected cell types known to be primary targets of NiV in humans, pulmonary endothelial (HPMEC) and bronchioalveolar epithelial (NCI-H358) cells, and the bat cell line PATGV1.12 at an MOI of 0.3 (Figure 4). NiV RNA synthesis was highest in HPMEC, where NiV-Bangladesh replicated the best, although a similar level of RNA and infectious virus particle production was observed for all 3 viruses (Figure 4, panel A). In accordance with virus entry studies (Figure 3), virus replication was also observed in PATGV1.12 (Figure 4, panels A and B). Differences among the 3

tested NiV isolates were observed only in NCI-H358, where NiV-Malaysia RNA synthesis was significantly increased (p<0.001 by 2-way analysis of variance) compared with NiV-Bangladesh, provoking remarkable cytopathic effects (Figure 4, panel C). The formation of giant multinucleated cells, a hallmark of NiV infection, were already visible at 24 hours postinfection in all cell types and further developed during the course of the infection (Figure 4, panels C-E). Vero cells showed the most visible cytopathic effects, probably because of their interferon incompetence (36).

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Replication of NiVs CSUR381 (Cambodia 2003 isolate), UMMC1 (NiV-Malaysia isolate), and SPB200401066 (NiV-Bangladesh isolate) in Vero, HPMEC, NCI-H358 (human bronchioalveolar cells), and PATGV1.12 (bat cells). A) Real-time reverse transcription PCR titer. Cells were infected at a multiplicity of infection of 0.3, and the production of the nucleocapsid gene was measured. Significance was measured by 2-way analysis of variance. B) Kinetics of infectious virus particle production in supernatant measured by Vero plaque assay. The average of 2 independent experiments is presented. C-E) Cytopathic effect of observed by light microscopy 48 h after infection with the NiV-Malaysia isolate (C), the NiV-Bangladesh isolate (D), and CSUR381 (E). Giant multinucleated cells are indicated with arrowheads. Scale bar indicates 25 µm. HPMEC, human pulmonary microvascular endothelial cell; NiV, Nipah virus. **p<0.01; ***p<0.001.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Pathogenicity of NiV CSUR381 (Cambodia 2003 isolate) and UMMC1 (NiV-Malaysia isolate) in golden hamsters (6 hamsters/group). A) Body weight. B) Body temperature. C) Survival. Survival between groups was not significantly different (Mantel-Cox test). NiV, Nipah virus.</figcaption></figure>

## Experimental Infection of Hamsters

We compared the pathogenicities of CSUR381 and the NiV-Malaysia isolate using the golden hamster animal model (37). We infected 6 hamsters with either CSUR381 or the NiV-Malaysia isolate and followed them for clinical signs of infection. At 6 days postinfection, the first neurologic signs (which included paralysis and trembling limbs) were observed in both groups; their presentation rapidly evolved toward breathing difficulties and prostration. Weight reductions were evident in several animals in the late stages of infection (Figure 5, panel A), and decreases in body temperature were found in a few hamsters (Figure 5, panel B). At 7 days postinfection, 100% lethality was observed in the CSUR381 group. In the Malaysia group, 1 animal survived until 10 days postinfection (Figure 5, panel C); however, the difference between the 2 groups was not significant. These results demonstrate similar lethality of the 2 analyzed NiV isolates, supporting our other data and suggesting CSUR381 has a high pathogenic potential.

## Discussion

In this study, we performed a molecular and genetic characterization of CSUR381, a NiV isolated from *P. lylei* bats in Cambodia. Furthermore, we analyzed its pathogenicity compared with those of 2 other NiV isolates derived from human patients from the Malaysia and Bangladesh outbreaks. Our results highly suggest that CSUR381 is part of the NiV-Malaysia genotype. Further phylogenetic comparisons with other NiV isolates demonstrated 83%-99% amino acid homology for each of the 6 structural proteins. In addition, the editing site of the phosphoprotein gene was preserved, suggesting possible production of the nonstructural V and W proteins known to be involved in counteracting the host innate immune system and thus contributing to pathogenicity of CSUR381 (32).

Our virus entry studies showed highly similar results among the NiVs tested. All isolates

entered *Pteropus* bat and human cell lines at similar levels; high conservation of the NiV entry receptors (ephrin-B2 and ephrin-B3) (38) might be responsible for the observed results. Our data also indicate that CSUR381 enters all tested cell types as well as the other 2 NiV isolates tested, suggesting that virus entry is not a limiting factor preventing CSUR381 spillover from bats to humans. In addition, all 3 tested NiV isolates infected cells and replicated in bat and human cell lines at similar levels. Results of infections with CSUR381 in hamsters additionally strengthened the notion that CSUR381 is possibly similar pathogenically to the tested NiV-Malaysia strain, which caused fatal outbreaks in Malaysia (1).

Although NiV has been shown to circulate in Cambodia (20,21), Thailand (39), and Vietnam (40), transmission to humans or domestic animals has not been reported in these countries. According to our results, the absence of detected outbreaks in this region cannot be attributed to lower pathogenicity of the circulating NiVs; our results suggest that other factors probably contribute. However, the NiV isolate presented in this report has been the only live NiV isolated in this region, and the existence of other NiVs with different pathogenic potentials cannot be excluded.

In Cambodia, *P. lylei* bats were found to often forage in residential areas and visit palm trees used in the region as a source of date palm sap; thus, opportunities abound for bats to interact with humans and livestock in this country (41). Bat colony migration toward urban sites is further enhanced by the presence of hunters in rural areas (42) and deforestation (causing consequent damage to roosting trees and food sources) (43). Contamination of palm sap, which is consumed raw by persons in the region, with bat urine, saliva, or feces was found to be a major route of NiV transmission to humans during annual outbreaks in Bangladesh (10).

Diverse agricultural practices in Southeast Asia could also play a role in NiV regional ecodynamics, potentially favoring easier NiV spillover in some countries over others. High intensity pig farming
was recognized as a major risk factor for outbreaks in Malaysia during 1998–1999; because of the low-scale pig production ongoing in Cambodia (44), the risk for NiV transmission from *Pteropus* spp. to domestic animals and humans in this country might be reduced.

Unrecognized NiV outbreaks might have occurred in Cambodia and neighboring countries; hospital-based surveillance in Bangladesh was shown to have missed nearly half of the NiV outbreaks in that country since the first reported virus emergence (45). Interdisciplinary approaches are certainly required to identify these outbreaks and the drivers of NiV emergence (46), and regular testing of patients with encephalitis in Cambodia and neighboring countries could provide additional insight. Our study contributes to the assessment of the risk for NiV outbreaks in Asia. Our findings can be used to help target adequate preventive measures, which could ultimately help reduce the risk for NiV emergence.

## Acknowledgments

We thank Jean-Marc Reynes and Pasteur Institute staff for providing us with the NiV Cambodia isolate, Pierre E. Rollin for the NiV Bangladesh isolate, and Doris Preininger, Anton Weissenbacher, and Tiergarten Schönbrunn for *P. giganteus* bat sampling. We thank Amelia Charlotte Coggon for English proofreading of the manuscript, and we also thank François Enchéry, Kévin Dhondt, Mathieu Iampietro, Sylvain Baize, and Géraldine Gourru-Lesimple for help initiating and finalizing this work.

This work was supported by LABEX ECOFECT (ANR-11-LABX-0048) of Lyon University, within the program Investissements d'Avenir (ANR-11-IDEX-0007) operated by the French National Research Agency, ANR-18-CE11-0014-02, Aviesan Sino-French Agreement on Nipah Virus Study, and the International Division of the Institut Pasteur in Paris (Actions Concertées Inter-Pasteurienne). J.F. was supported by the doctoral fellowship CIFRE-Défense operated by the Direction Générale de l'Armement.

## About the Author

Ms. Gaudino is a graduate student at the International Centre for Infectiology Research in Lyon, France. Her main research interests include the study of mechanisms of pathogenicity and epidemiology of viral zoonoses.

## References

1. Chua KB, Bellini WJ, Rota PA, Harcourt BH, Tamin A, Lam SK, et al. Nipah virus: a recently emergent deadly paramyxovirus. Science. 2000;288:1432–5. https://doi.org/10.1126/science.288.5470.1432

2. Arunkumar G, Chandni R, Mourya DT, Singh SK, Sadanandan R, Sudan P, et al.; Nipah Investigators People and Health Study Group. Outbreak investigation of Nipah virus disease in Kerala, India, 2018. J Infect Dis. 2019;219:1867–78. https://doi.org/10.1093/infdis/jiy612

3. Luby SP, Gurley ES, Hossain MJ. Transmission of human infection with Nipah virus. Clin Infect Dis. 2009;49:1743–8. https://doi.org/10.1086/647951

4. Sweileh WM. Global research trends of World Health Organization's top eight emerging pathogens. Global Health. 2017;13:9. https://doi.org/10.1186/s12992-017-0233-9

5. Murray K, Selleck P, Hooper P, Hyatt A, Gould A, Gleeson L, et al. A morbillivirus that caused fatal disease in horses and humans. Science. 1995;268:94–7. https://doi.org/10.1126/science.7701348

6. Marsh GA, de Jong C, Barr JA, Tachedjian M, Smith C, Middleton D, et al. Cedar virus: a novel henipavirus isolated from Australian bats. PLoS Pathog. 2012;8:e1002836. http://dx.doi.org/10.1371/journal.ppat.1002836

7. Drexler JF, Corman VM, Müller MA, Maganga GD, Vallo P, Binger T, et al. Bats host major mammalian paramyxoviruses. Nat Commun. 2012;3:796. https://doi.org/10.1038/ncomms1796

8. Wu Z, Yang L, Yang F, Ren X, Jiang J, Dong J, et al. Novel henipa-like virus, Mojiang paramyxovirus, in rats, China, 2012. Emerg Infect Dis. 2014;20:1064–6. https://doi.org/10.3201/eid2006.131022

9. Middleton DJ, Morrissy CJ, van der Heide BM, Russell GM, Braun MA, Westbury HA, et al. Experimental Nipah virus infection in pteropid bats (*Pteropus poliocephalus*). J Comp Pathol. 2007;136:266–72. https://doi.org/10.1016/j.jcpa.2007.03.002

10. Luby SP, Rahman M, Hossain MJ, Blum LS, Husain MM, Gurley E, et al. Foodborne transmission of Nipah virus, Bangladesh. Emerg Infect Dis. 2006;12:1888–94. https://doi.org/10.3201/eid1212.060732

11. Middleton DJ, Westbury HA, Morrissy CJ, van der Heide BM, Russell GM, Braun MA, et al. Experimental Nipah virus infection in pigs and cats. J Comp Pathol. 2002;126:124–36. https://doi.org/10.1053/jcpa.2001.0532

12. Luby SP, Gurley ES, Hossain MJ. Transmission of human infection with Nipah virus. Clin Infect Dis. 2009;49:1743–8. https://doi.org/10.1086/647951

13. Sejvar JJ, Hossain J, Saha SK, Gurley ES, Banu S, Hamadani JD, et al. Long-term neurological and functional outcome in Nipah virus infection. Ann Neurol. 2007; 62:235–42. https://doi.org/10.1002/ana.21178

14. Tan CT, Goh KJ, Wong KT, Sarji SA, Chua KB, Chew NK, et al. Relapsed and late-onset Nipah encephalitis. Ann Neurol. 2002;51:703–8. https://doi.org/10.1002/ana.10212

15. Enchéry F, Horvat B. Understanding the interaction between henipaviruses and their natural host, fruit bats: paving the way toward control of highly lethal infection in humans. Int Rev Immunol. 2017;36:108–21. https://doi.org/10.1080/08830185.2016.1255883

16. Halpin K, Young PL, Field HE, Mackenzie JS. Isolation of Hendra virus from pteropid bats: a natural reservoir of Hendra virus. J Gen Virol. 2000;81:1927–32. https://doi.org/10.1099/0022-1317-81-8-1927

17. de Araujo J, Lo MK, Tamin A, Ometto TL, Thomazelli LM, Nardi MS, et al. Antibodies against henipa-like viruses in Brazilian bats. Vector Borne Zoonotic Dis. 2017;17:271–4. https://doi.org/10.1089/vbz.2016.2051

18. Luby SP. The pandemic potential of Nipah virus. Antiviral Res. 2013;100:38–43. https://doi.org/10.1016/j.antiviral.2013.07.011
19. Mathieu C, Horvat B. Henipavirus pathogenesis and antiviral approaches. Expert Rev Anti Infect Ther. 2015;13:343-54. https://doi.org/10.1586/14787210. 2015.1001838

20. Olson JG, Rupprecht C, Rollin PE, An US, Niezgoda M, Clemins T, et al. Antibodies to Nipah-like virus in bats (*Pteropus lylei*), Cambodia. Emerg Infect Dis. 2002;8:987-8. https://doi.org/10.3201/eid0809.010515

21. Reynes J-M, 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. https://doi.org/10.3201/ eid1107.041350

22. Thomas HA. BioEdit, a user friendly biological sequence alignment editor and analysis program for Windows 95/98/ NT. Nucleic Acids Symp Ser. 1999;41:95-8.

23. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870-4. https://doi.org/ 10.1093/molbev/msw054

24. Krump-Konvalinkova V, Bittinger F, Unger RE, Peters K, Lehr HA, Kirkpatrick CJ. Generation of human pulmonary microvascular endothelial cell lines. Lab Invest. 2001;81:1717-27. https://doi.org/10.1038/labinvest.3780385

25. Mlíkovský J. Correct name for the Indian flying fox (*Pteropodidae*). Vespertilio. 2012;16:203-4.

26. Saccone C, Attimonelli M, Sbisà E. Structural elements highly preserved during the evolution of the D-loop-containing region in vertebrate mitochondrial DNA. J Mol Evol. 1987;26:205-11. https://doi.org/10.1007/BF02099853

27. Dool SE, Puechmaille SJ, Foley NM, Allegrini B, Bastian A, Mutumi GL, et al. Nuclear introns outperform mitochondrial DNA in inter-specific phylogenetic reconstruction: lessons from horseshoe bats (*Rhinolophidae: Chiroptera*). Mol Phylogenet Evol. 2016;97:196-212. https://doi.org/10.1016/ j.ympev.2016.01.003

28. Reynard O, Volchkov VE. Characterization of a novel neutralizing monoclonal antibody against Ebola virus GP. J Infect Dis. 2015;212(Suppl 2):S372-8. https://doi.org/ 10.1093/infdis/jiv303

29. Talekar A, Pessi A, Glickman F, Sengupta U, Briese T, Whitt MA, et al. Rapid screening for entry inhibitors of highly pathogenic viruses under low-level biocontainment. PLoS One. 2012;7:e30538. https://doi.org/10.1371/ journal.pone.0030538

30. Whitt MA. Generation of VSV pseudotypes using recombinant ΔG-VSV for studies on virus entry, identification of entry inhibitors, and immune responses to vaccines. J Virol Methods. 2010;169:365-74. http://dx.doi.org/10.1016/j.jviromet.2010.08.006

31. Mathieu C, Guillaume V, Sabine A, Ong KC, Wong KT, Legras-Lachuer C, et al. Lethal Nipah virus infection induces rapid overexpression of CXCL10. PLoS One. 2012;7:e32157. https://doi.org/10.1371/journal.pone.0032157

32. Satterfield BA, Cross RW, Fenton KA, Agans KN, Basler CF, Geisbert TW, et al. The immunomodulating V and W proteins of Nipah virus determine disease course. Nat Commun. 2015;6:7483. https://doi.org/10.1038/ ncomms8483

33. Guillaume V, Aslan H, Ainouze M, Guerbois M, Wild TF, Buckland R, et al. Evidence of a potential receptor-binding site on the Nipah virus G protein (NiV-G): identification of globular head residues with a role in fusion promotion and their localization on an NiV-G structural model. J Virol. 2006;80:7546-54. https://doi.org/10.1128/JVI.00190-06

34. Xu K, Broder CC, Nikolov DB. Ephrin-B2 and ephrin-B3 as functional henipavirus receptors. Semin Cell Dev Biol. 2012;23:116-23. https://doi.org/10.1016/ j.semcdb.2011.12.005

35. Porotto M, Rockx B, Yokoyama CC, Talekar A, Devito I, Palermo LM, et al. Inhibition of Nipah virus infection in vivo: targeting an early stage of paramyxovirus fusion activation during viral entry. PLoS Pathog. 2010;6:e1001168. https://doi.org/10.1371/journal.ppat.1001168

36. Desmyter J, Melnick JL, Rawls WE. Defectiveness of interferon production and of rubella virus interference in a line of African green monkey kidney cells (Vero). J Virol. 1968;2:955-61.

37. Wong KT, Grosjean I, Brisson C, Blanquier B, Fevre-Montange M, Bernard A, et al. A golden hamster model for human acute Nipah virus infection. Am J Pathol. 2003;163:2127-37. https://doi.org/10.1016/S0002-9440 (10)63569-9

38. Bossart KN, Tachedjian M, McEachern JA, Crameri G, Zhu Z, Dimitrov DS, et al. Functional studies of host-s pecific ephrin-B ligands as henipavirus receptors. Virology. 2008;372:357-71. https://doi.org/10.1016/j.virol.2007.11.011

39. 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. https://doi.org/10.3201/eid1112.050613

40. Hasebe F, Thuy NTT, Inoue S, Yu F, Kaku Y, Watanabe S, et al. Serologic evidence of Nipah virus infection in bats, Vietnam. Emerg Infect Dis. 2012;18:536-7. https://doi.org/ 10.3201/eid1803.111121

41. Choden K, Ravon S, Epstein JH, Hoem T, Furey N, Gely M, et al. *Pteropus lylei* primarily forages in residential areas in Kandal, Cambodia. Ecol Evol. 2019;9:4181-91. https://doi.org/10.1002/ece3.5046

42. Thanapongtharm W, Linard C, Wiriyarat W, Chinsorn P, Kanchanasaka B, Xiao X, et al. Spatial characterization of colonies of the flying fox bat, a carrier of Nipah virus in Thailand. BMC Vet Res. 2015;11:81. http://dx.doi.org/ 10.1186/s12917-015-0390-0

43. Thanapongtharm W, Paul MC, Wiratsudakul A, Wongphruksasoong V, Kalpravidh W, Wongsathapornchai K, et al. A spatial assessment of Nipah virus transmission in Thailand pig farms using multi-criteria decision analysis. BMC Vet Res. 2019;15:73. https://doi.org/10.1186/ s12917-019-1815-y

44. Huynh TTT, Aarnink AJA, Drucker A, Verstegen MWA. Pig production in Cambodia, Laos, Philippines, and Vietnam: a review. Asian J Agric Dev. 2006;3:69-90.

45. Hegde ST, Salje H, Sazzad HMS, Hossain MJ, Rahman M, Daszak P, et al. Using healthcare-seeking behaviour to estimate the number of Nipah outbreaks missed by hospital-based surveillance in Bangladesh. Int J Epidemiol. 2019;48:1219-27. https://doi.org/10.1093/ije/dyz057

46. Daszak P, Zambrana-Torrelio C, Bogich TL, Fernandez M, Epstein JH, Murray KA, et al. Interdisciplinary approaches to understanding disease emergence: the past, present, and future drivers of Nipah virus emergence. Proc Natl Acad Sci U S A. 2013;110(Suppl 1):3681-8. https://doi.org/10.1073/ pnas.1201243109

Address for correspondence: Branka Horvat, CIRI, INSERM U1111, 21 Ave Tony Garnier, 69007, Lyon, France; email: branka.horvat@inserm.fr
# High Pathogenicity of Nipah Virus from *Pteropus lylei* Fruit Bats, Cambodia

## Appendix

<table><caption>**Appendix Table.** Genetic pairwise comparisons of NiV nonstructural proteins V, W, and C between NiV/Cambodia and other NiV isolates*</caption><thead><tr><th rowspan="2">NiV, GenBank accession no.</th><th colspan="3">Nonstructural NiV proteins, % nucleotide identity/% homology of deduced amino acid</th></tr><tr><th>V</th><th>W</th><th>C</th></tr></thead><tbody><tr><td>Nipah/Malaysia/2000/human, NC_002728.1</td><td>97.2/95.7</td><td>97.3/95.8</td><td>98.6/98.2</td></tr><tr><td>Nipah/Malaysia/2010/Pvampyrus, FN869553.1</td><td>97.1/95.5</td><td>96.6/94.0</td><td>99.4/100</td></tr><tr><td>Nipah/Bangladesh/2004/human, AY988601.1</td><td>92.4/89.4</td><td>92.5/89.6</td><td>97.6/97.0</td></tr><tr><td>Nipah/Bangladesh/2008/human, JN808863.1</td><td>91.8/88.5</td><td>91.9/88.7</td><td>97.4/97.0</td></tr><tr><td>Nipah/India/2007/human, FJ513078.1</td><td>92.1/88.9</td><td>92.2/89.1</td><td>97.6/97.0</td></tr><tr><td>Nipah/India/2018/human, MH396625.1</td><td>91.7/88.0</td><td>91.8/88.0</td><td>97.2/96.4</td></tr></tbody></table>

A

<figure><img src="image_9.png" /></figure>

0.1

B

<figure><img src="image_10.png" /></figure>

0.02

C

<figure><img src="image_11.png" /></figure>

0.2

D

<figure><img src="image_12.png" /></figure>

0.05

E

<figure><img src="image_13.png" /></figure>

0.05

**Appendix Figure 1.** Maximum-likelihood phylogenetic trees based on genes coding for *Henipavirus*. A) Phosphoprotein, B) Matrix protein, C) Fusion protein, D) Attachment glycoprotein, and E) Polymerase. The Cambodia 2003 Nipah virus isolate CSUR381 (GenBank MK801755, red triangle) clusters with NiV-Malaysia isolates. Bootstrap statistical support is shown on branch nodes. Hasegawa-Kishino-Yano model and Tamura-3-parameter were calculated as the best DNA model to conduct the analysis for these genes. The phylogenetic trees are drawn to scale and scale bar represents branch lengths measured in the number of substitutions per site.

<footnote>* NiV, Nipah virus.</footnote>

<footnote>†The nucleotide identity and deduced amino acid homology were calculated using p-distance method.</footnote>
<figure><img src="image_14.png" /><figcaption>**Appendix Figure 2.** Amino acid multiple alignment of the attachment glycoprotein (G) of three Nipah virus (NiV) isolates used in the study. Specific amino acid changes were found in NiV/Cambodia G sequence (V24→I, R248→K, G327→D, I408→V, N478→S). The glycosylation site in G attachment protein (N529/Q530/T531) is conserved among the three tested isolates (marked in red). Polymorphism at Ephrin B2 and B3 potential binding sites was not observed (W504, Q530, T531, A532, and N557, E505, E533, respectively).</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Appendix Figure 3.** Amino acid multiple alignment of the fusion protein (F) of three Nipah virus (NiV) isolates used in the study. Specific NiV/C amino acid changes were found in the beginning of F2 subunit (I4→V, I15→M). Fusion protein cleavage region (N104-R109, marked in red) and fusion peptide (L110-V133, marked in blue) are preserved among all analyzed NiV isolates. Predicted N-terminal and C-</figcaption></figure>
terminal heptad repeat regions (HRN and HRC) are marked in green and orange respectively. V159→I mutation was observed in HRN region, while no variability was found within the HRC.