Abstract
Background We report here a Nipah virus (NiV) outbreak in Kozhikode district of Kerala state, India which had caused fatal encephalitis in an adolescent male and the outbreak response which led to the successful containment of the disease and the related investigations.
Methods Quantitative real-time RT-PCR, ELISA based antibody detection and whole genome sequencing were performed to confirm the Nipah virus infection. Contacts of the index case were traced and isolated based on risk categorization. Bats from the areas near the epicenter of the outbreak were sampled for throat swabs, rectal swabs and blood samples for Nipah virus screening by real time RT-PCR and anti-Nipah virus bat IgG ELISA. Plaque reduction neutralization test was performed for the detection of neutralizing antibodies.
Results Nipah viral RNA and anti-NiV IgG antibodies were detected in the serum of the index case. Rapid establishment of an onsite NiV diagnostic facility and contact tracing helped in quick containment of the outbreak. NiV sequences retrieved from the clinical specimen of the index case formed a sub-cluster with the earlier reported Nipah I genotype sequences from India with more than 95% similarity. Anti-NiV IgG positivity could be detected in 21% of Pteropus medius and 37.73% of Rousettus leschenaultia. Neutralizing antibodies against NiV could be detected in P.medius.
Conclusions Stringent surveillance and awareness campaigns needs to be implemented in the area to reduce human-bat interactions and minimize spill over events which can lead to sporadic outbreaks of NiV.
Introduction
Nipah virus (NiV) causes a highly lethal disease with acute severe encephalitis and acute respiratory distress syndrome in humans. NiV is a Paramyxovirus that was identified for the first time during an outbreak of severe encephalitis among the pig-farmers in Malaysia during 19981. The virus is transmitted to humans by direct contact with the respiratory secretions or body fluids of infected animals such as bats and pigs or by consumption of contaminated fruits/palm sap. Both animal-to-human and human-to-human transmission have been documented 2.
The disease has been enlisted as a priority disease in the Research and Development blueprint of the World Health Organization (WHO) from the year 2015. South-East Asia is considered as an endemic region for the virus2. In Malaysia and Singapore outbreak, humans acquired the infection through the contact with infected pigs1,3. Subsequently, an outbreak with encephalitis among human population was observed in Meherpur, Bangladesh during 2001 where the source of infection was traced down to drinking the contaminated raw palm sap or climbing the trees coated with bat excrement4. India has witnessed two outbreaks of Nipah virus encephalitis in the eastern state of West Bengal, bordering Bangladesh5,6. The first outbreak of NiV was observed among hospital visitors and health workers after coming in contact with a hospitalized patient in Siliguri during January and February 20015. In the year 2007, a second NiV outbreak was reported from Nadia district of West Bengal6. A case fatality of 70-100% was observed during these two outbreaks5,6. Since 2010, Indian Council of Medical Research-National Institute of Virology (ICMR-NIV), Pune has taken up the surveillance of NiV in bat populations across the country. During this, the presence of NiV was detected among Pteropus medius from Myanaguri, West Bengal in 2010 and Cooch Bihar district, West Bengal and Dhubri district, Assam in 20157,8.
After a decade of the last outbreak, a dreadful emergence of NiV was observed in Kozhikode district, Kerala State during May 2018 with case fatality rate of 89%9. The outbreak was contained with the quick actions of national and state health system10. A year later, a young male from Ernakulam district, Kerala got infected with NiV. Due to the early detection; the further spread of the virus was quickly curtailed. A detailed outbreak investigation to find the source of NiV infection was carried out by ICMR-NIV, Pune during these outbreaks which showed the presence of NiV and anti-NiV antibodies in Pteropus medius11. The unprecedented emergence of NiV in the Kerala state in 2018 and 2019 raised a serious concern. Recently a NiV outbreak was reported in Kozhikode Kerala in September 2021, where an adolescent male who presented with acute encephalitis and tested positive for NiV, succumbed to the infection.
Here, we describe the 2021 NiV outbreak management in Kozhikode district, Kerala, India with emphasis on the field laboratory set up and quick diagnosis along with the bat survey to trace the source of infection. Immediately after confirming the NiV infection in September 2021, outbreak containment response was initiated in the state of Kerala.
Materials and Methods
Case history
In August 2021, an adolescent male (index case) resident of Kozhikode district, Kerala state, India developed low grade fever (Supplementary Figure.1). His family took him to a nearby private clinic (Hospital-1) and sought treatment for fever. In August 2021 his condition deteriorated and he was transferred to another hospital (Hospital-2). In September 2021, the patient’s condition deteriorated further and he developed symptoms of acute encephalitis and myocarditis. On the request of his family, he was transferred to a tertiary care hospital in Kozhikode (Hospital-3). Magnetic resonance imaging of the brain showed multiple small infarcts in the cerebellum, cerebrum, medulla oblongata and pons. With a high suspicion of NiV infection, his clinical samples including plasma, EDTA blood, serum, and cerebrospinal fluid (CSF) samples, endotracheal (ET) secretion and bronchial wash were sent to ICMR-NIV in September 2021.
Bat trapping and sample collection
For understanding the source of NiV infection and considering a brief history given by parent of the index case regarding the consumption of fruit from the orchard near the house, bat sampling was done during the month of September 2021. The sampling was performed with the prior approval from the Institutional Animal Ethics Committee, Institutional Biosafety Committee of ICMR-NIV, Pune and the Principal Chief Conservator of Forests, Government of Kerala. Four roosting sites were chosen from nearby vicinity of the index case house for sample collection (Supplementary Figure.1). These sites were in Kodiyathur, Cheruvadi, Omassery and Thamarassery) in the state of Kerala, India. Bats were trapped using the mist nets as described earlier11. Body weight, sex, secondary sexual characters and forearm length were noted. Blood (n=91), throat (n=102) and rectal swab (n=102) samples were collected from the trapped bats following isoflurane anesthesia. Species identification was performed by mitochondrial Cytochrome b gene PCR as described earlier.12
Real time quantitative RT-PCR
Real-time reverse transcriptase polymerase chain reaction (qRT-PCR) was performed on the samples for the NiV diagnosis as described earlier13. One ml of serum/swab samples were used for RNA extraction in an extraction machine using Magmax Viral RNA isolation kit as per manufacturer’s instructions. Five microlitre extracted RNA was used for qRT-PCR using primers designed for Nipah NP gene.13 For SARS-CoV-2 detection, qRT-PCR was performed for throat/nasal swab samples using primers for E gene as described earlier14.
Anti-Nipah human IgM ELISA
The assay was performed as described earlier11. Briefly, the microtiter plates were coated with 1: 100 diluted anti-human IgM antibodies (Sigma, USA) overnight at 40C. Following blocking of the wells with 2% bovine serum albumin for 2 hours at room temperature, 100 µl of 1:100 diluted heat inactivated serum samples were added to duplicate wells and was incubated at 37°C for one hour. After washing with wash buffer for four times, inactivated NiV antigen (positive antigen) and Vero-CCL81 cell lysate (negative antigen) were added and incubated for one hr at 37°C. The plates were washed and anti-NiV biotinylated antibodies were added and incubated for one hour. This was followed by washing and addition of Streptavidine Horseradish peroxidase (HRP) and incubation for 30 min at 37°C. As substrate, one hundred µl of 3, 3’,5,5’-Tetramethylbenzidine (TMB) was added and incubated for 10 minutes. The reaction was stopped by 1 N sulphuric acid, and the absorbance was measured at 450 nm. The assay cut off were an Optical Density (OD) value more than or equal to 0.20 and the positive/negative (P/N) ratio more than or equal to 1.5.
Anti-Nipah human IgG ELISA
Briefly, microtiter plates were coated with 0.1 ml of inactivated NiV antigen and normal Vero CCL81 cell lysate as negative antigen in carbonate buffer (pH 9.2, 0.025 M) overnight at 4°C. The wells were blocked with liquid plate sealer (Candor, USA) for 2 hours at room temperature. The plates were washed five times using 10 mM PBS (pH 7.4) with 0.1% Tween-20 (Sigma, USA) at the end of each step. Hundred micro litre of 1:100 diluted heat inactivated human serum samples were added to both positive and negative antigen wells and was incubated at 37°C for one hour. After washing, 100 micro litre of 1:15000 dilution of anti-human IgG HRP antibodies were added and incubated for one hr at 37°C. After washing, 100 µl of TMB substrate was added and incubated for 10 minutes. The reaction was stopped by 1 N sulphuric acid and plates were read at 450 nm. OD value more than or equal to 0.35 and the P/N ratio more than or equal to 1.5 was considered as positive.
Whole Genome Sequencing (WGS)
To identify the NiV genotype, WGS was carried out on different clinical samples obtained from the patient. RNA was extracted from the clinical samples and WGS was performed using the methods descried earlier15. The viral reads generated were analyzed using reference-based mapping, performed in CLC Genomics Workbench version 21.0.4. In order to retrieve the complete genome sequence of the virus all the generated reads were mapped to the reference genome. Phylogenetic analyses and the amino acid variations of the retrieved NiV sequence with the other representative NiV sequences of earlier outbreaks were performed. SARS-CoV-2 positive samples were sequenced as per the method described earlier to identify the lineage in circulation16.
Anti-NiV bat IgG ELISA
The assay was performed as described earlier15. NiV antigen/uninfected Vero CCL81 cell lysate antigen coated microtiter plates were used. Diluted bat sera samples were added and incubated for one hour in the microtiter plates at 37°C. The plates were washed four times following the incubation period using wash buffer. Hundred microliter of anti-bat IgG HRP conjugate (1:2000 dilution) was added and incubated. TMB substrate (Clinical Science Product incorporation NeA blue, USA) was added after washing and was again incubated at 37°C for 15 minutes. The reaction was stopped and plates were read at 450 nm.
Plaque reduction neutralization test
Heat inactivated bat serum samples were mixed with NiV (GenBank accession number: MH523642) containing 50 plaque forming units in 1:1 ratio so as to make a final 10-fold dilution of the serum virus mixture. Anti-Nipah IgG positive mice serum was used as positive control and Anti-Nipah IgG negative mice serum was used as negative control. The mixture was incubated for 1 hour and was added to 24-well tissue culture plate containing a confluent monolayer of Vero CCL-81 cells. The plate was incubated in a CO2 incubator at 37°C for one hour and an overlay medium containing 2% Carboxy methyl cellulose in 2X Minimal Essential Media with 2% fetal bovine serum was added after removing the inoculum. The plate was further incubated at 37°C in CO2 incubator for 3 days. After removing the overlay medium, the plate was washed and stained with amido black. The plaques were counted. The titre was defined as the highest serum dilution that resulted in 50 per cent (PRNT50) reduction in the number of plaques.
Results
Detection and confirmation of NiV infection
NiV infection was confirmed by qRT-PCR and by the detection of anti-Nipah IgM antibodies in serum sample of the patient. Virus isolation attempts from the samples in Vero CCL-81 cells were not successful. In September 2021, NiV outbreak in Kerala was declared by the Ministry of Health and Family Welfare, Government of India. The patient succumbed to the infection in September 2021. Contact tracing and isolation of high-risk contacts, augmentation of laboratory testing capacity, bat sampling and laboratory investigations were undertaken following the event as described below.
Establishment of an on-site field NiV diagnostic facility
A team from ICMR-NIV Pune had set up a field diagnostic laboratory in the Department of Microbiology, Government Medical College, Kozhikode, Kerala by September 2021 following all the essential biosafety guidelines and standard operating procedures. The field laboratory screened a total of 125 NiV suspected asymptomatic contacts as well as non-epidemiologically linked suspected NiV cases by qRT-PCR.10
Risk categorization and contact tracing
After declaration of NiV outbreak, systematic field investigations were undertaken to identify the epidemiologically linked close contacts, including health care workers, family members, neighbours, bystanders, etc. The close contacts were classified into primary contacts and secondary contacts and were further grouped into high-risk and low-risk contacts. The high-risk category included individuals with either a history of direct contact with body fluids (blood, urine, saliva, vomitus, etc.) of the confirmed NiV case/a probable case that died without laboratory confirmation or having spent about 12 hour nearby or in closed space with confirmed NiV case. The low-risk contacts were categorized as those having contact with the confirmed NiV case through touching or contact with clothes, linen, or any other items. Samples of symptomatic contacts were shipped to the Biosafety Level 4 (BSL-4) laboratory of ICMR-NIV, Pune for diagnosis. All asymptomatic contacts were tested on site at the field laboratory. Considering the ongoing COVID-19 pandemic, all the close contacts were also screened for SARS-CoV-2.
A total of 240 contacts were listed and among them 64 close contacts [33 female/31 male] were identified and grouped into primary high-risk (n=50) and low-risk (n=9) contacts; secondary high-risk (n=3) and low-risk (n=2) contacts (Figure. 1). Out of 59 primary contacts, 40 were asymptomatic while all the 5 secondary contacts were asymptomatic (Supplementary Table.1). All the 64 close contacts tested negative for NiV by qRTPCR and anti-NiV IgM and IgG ELISA. The throat/nasal swab of the 12 close contacts (symptomatic-8 and asymptomatic-4) were found positive for SARS-CoV-2 by qRT-PCR (Supplementary Table.1). On sequencing of SARS-CoV-2 positive samples (n=12), Delta variant (B.1.167.2) and its derivatives (AY.26) were detected in ten and two cases respectively.
Genomic characterization of NiV from clinical specimens of the index case
A phylogenetic analysis was performed for the retrieved NiV sequence (17066 nucleotide) with the other representative NiV sequences of the earlier outbreaks (Figure.2, Supplementary table 2). The retrieved sequence clearly segregated from the Bangladesh NiV sequences and clustered into earlier described Indian (‘I’) genotype. The retrieved sequence showed 99.62 and 99.51 percent nucleotide similarity (PNS) with sequences obtained from human samples during 2018 NiV outbreak and P. medius samples during 2019 outbreak respectively. A detailed PNS of the representative NiV sequences with that of the NiV sequence from Malaysia (Accession Number: NC_002728.1) is provided in the Table.1.
The different genes of NiV showed amino acid variation between the Bangladesh and I genotype cluster. The changes observed were in the N gene (S503N, P520S, E752G, Q758E, R818H, I820L, T919N, Q982R, A1162T, and G1216D), F gene (I15L) and G gene (R344M, I384V, V427I) of I genotype of the NiV compared to the NC_002728.1. Similarly, the changes were observed in the N gene (R505K, S900G, D921N), G gene (R344K, K386E, T498K) and L gene (R1262K, N 237D) when the Bangladesh (BD) NiV sequences were compared to the NC_002728.1 (Table.2).
Detection of anti-NiV IgG antibodies in bat samples
The bat species sampled in the study included Pteropus medius (n=38), Rousettus leschenaultia (n=63) and Pipistrellus sp. (n=1) (Supplementary Table.3). All the bat samples were found to be negative for Nipah viral RNA (Table.3). The serum samples of P. medius (n=8) and R. leschenaultia (n=20) were tested positive for anti-NiV IgG antibodies. The P. medius samples from two sites ie., Kodiyathur and Thamarassery showed positivity of 20% and 56% by ELISA respectively and was further confirmed by PRNT. Two samples were excluded from the assay due to insufficient sample quantity. R leschenaulti samples which showed seropositivity by ELISA were found negative for neutralizing antibodies (Table.3).
Discussion
Nipah virus outbreaks have been reported from Malaysia, Singapore, Bangladesh and India with a range of clinical presentations and case fatality rate of 40 to 100%.1-7 Malaysia and Singapore had reported single NiV outbreak episode whereas Bangladesh reports annual outbreaks of NiV4. In India NiV outbreaks have been localized to two regions ie., in the North eastern region of West Bengal and in the southernmost region of Kerala which are separated by a distance of more than 2000 km.5,6,9-11 NiV outbreaks are mostly sporadic and the magnitude of the outbreak can be restricted by prompt public health response and containment measures. In 2018, the first NiV outbreak was reported in Kerala state with significant number of losses of lives (91% fatality) and in the year 2019 and 2021 outbreaks in Kerala reported only single case without any further human to human transmission.9,11 The diagnosis of the NiV infection is difficult considering the overlapping presentations of the acute respiratory infections and encephalitis syndromes. In the current SARS-CoV-2 pandemic scenario, the diagnosis became even more challenging considering the overlapping clinical features. The quick outbreak containment response and the biosafety practices followed during the COVID-19 pandemic period in Kerala state might have helped in preventing further transmission and restricting the outbreak to a single case.
The role of intermediate hosts like pigs and bats has been demonstrated in the NiV transmission cycle in the previous reported outbreaks.1,3,4 Pteropus genus of bats appears to be the major reservoir of the NiV.8,18,20 P. medius is the only Pteropus genus bat present in the Indian subcontinent.19 Pteropus bats has shown NiV RNA positivity or seropositivity from Indian states like West Bengal, Assam, Haryana and Kerala indicating the risk of spill-over.7,8,11,20 Viral genome recovered from the current outbreak also clustered with the previously reported NiV genotype from Kerala in bats and humans samples. This suggests a stable genotype circulating locally in bat population in Kerala. Research has shown NiV evolves at a slower rate compared to other similar RNA viruses.21,22 In 2018, the outbreak occurred in Perambra, Kozhikode which is about 40 km far from the location of the present outbreak.9 A total of 23 cases and 21 deaths were reported. 21% of the Pteropus medius bats surveyed in the vicinity of the index case residence. In 2019, NiV RNA positivity was documented in P medius bats from Thodupuzha which is more than 200 km from Kozhikode district.11 Pteropus species movement could also play an important role in virus spread, but in India such records are not available. Home ranges of P. medius appear to be smaller than P. vampyrus, another frugivorous bat in Malaysia and this depend on food availability18. Studies on the Pteropus species movement and connectivity among the bat populations could help us in understanding the potential of virus spread to bat colonies of adjacent areas.
NiV RNA could not be detected in any of the bats samples in the present study. Previous studies have also reported a low PCR positivity in P medius. 18,23,24. Virus shedding in bats is driven by multiple factors like individual immune status, pregnancy, virus recrudescence, stress etc.18,24 Re-infection in adult bats are possible in approximately 7 years24. We detected 25% and 3.2% NiV positivity in bats from outbreak regions of Kerala in 2018 and 2019.9,11 High seroprevalence in the bat colonies can dampen the virus transmission.18 Seasonality of NiV outbreaks have been reported is linked to the breeding season and fruit harvesting season.25 Unlike the previous reported outbreaks in the month of May and June in Kerala, the present outbreak was in August end. A cross sectional spatial study conducted between 2006 and 2012 in Bangladesh reports the ability of NiV shedding by P. medius throughout the year.18
During the current outbreak investigations, we could detect antibodies in a total of 21% P. medius bats. If we further classify this positivity, site wise, a 20% positivity was observed in Kodiyathur, a site within 1 km distance from the index case house and 55.6% in Thamarassery, a site 18 km far from this outbreak. During the investigation, 2 juvenile P. medius bats were found positive for anti-NiV IgG antibodies. The seroprevalence in the juveniles is an indicator of enzootic cycle of NiV in recent past. Waning of immunity with time and loss of maternal antibodies could affect the transmission dynamics within the bat colonies.18,24 Local NiV epizootics in bats contribute to the outbreaks in Bangladesh.18 Sporadic natures of outbreaks in case of NiV are explained by the NiV dynamics in the resident bat colonies and the bat-human interface resulting in spill over.
Non-pteropid bats have shown NiV seropositivity in many countries.17,23,26,27 In India, other than P medius, bats like R. leschenaultia and P. pipistrellus have shown seropositivity by ELISA.17 Interspecies transmission could be possible with other frugivorous bats as they share same habitat for food resources. Similarly, in the present study, we could detect seropositivity in R. leschenaultia bats by ELISA. But absence of NAb against NiV in these bats points to the possibility of detection of related Henipaviruses in these bats. Similar observations are reported in Vietnam where the ELISA positive samples of R. leschenaultia failed to show any NAb, indicating the probability of cross reacting antibodies against non-neutralizing epitopes26. Antibody positivity by ELISA and Western Blot in R leschenaultia has been reported from China, but neutralization studies were not performed in the samples27.
The findings from the current outbreak, suggest that spill-over of NiV infection in humans is sporadic and the seroprevalence in bats indicates the prevalence of the NiV infections in the P. medius population. Stringent surveillance measures and awareness campaigns to limit interactions with bats needs to be implemented in the area as the NiV transmission dynamics depends on multiple host factors including human behavior and human-bat interface. For early detection and containment of NiV outbreaks, it is critical to strengthen human surveillance for Acute Encephalitis Syndrome and Severe Acute Respiratory Infection, including testing for NiV in susceptible areas.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Funding
ICMR supported the funding for this study under the mission-oriented project to the ICMR-National Institute of Virology, Pune. The funders had no role in study design, data collection or interpretation, or the decision to submit the work for publication. The findings and conclusions in this study are of the authors.
Acknowledgments
Authors extend gratitude to Veena George [Hon’ble Minister for Health and Social Justice, Kerala], for her exemplary leadership and efficient coordination of the Nipah virus disease control activities, and “The Team Kerala Health,” the district administration, press, media, and people for their efforts in streamlining the public health responses. The authors are thankful to Prof. Balram Bhargava, Director General, ICMR and Dr Samiran Panda, Scientist G & Head, Epidemiology and Communicable Diseases Division, ICMR for their support. Authors are also thankful for the staff of the Administration section, ICMR-NIV Pune for their support during the outbreak response. We also acknowledge the help and support received from Dr. Arun Sacharia, Chief Veterinarian (Wild life) and Dr. Arun Sathyan, Assistant Veterinarian, Thamaraserry during bat survey.