Investigating Rare Risk Factors for Nipah Virus in Bangladesh: 2001-2012 HHS Public Access
Abstract
Human Nipah encephalitis outbreaks have been identified almost yearly in Bangladesh since 2001. Though raw date palm sap consumption and person-to-person contact are recognized as major transmission pathways, alternative pathways of transmission are plausible and may not have been identified due to limited statistical power in each outbreak. We conducted a risk factor analysis using all 157 cases and 632 controls surveyed in previous investigations during 2004-2012 to identify exposures independently associated with Nipah, since date palm sap was first asked about as an exposure in 2004. To further explore possible rare exposures, we also conducted in-depth interviews with all cases, or proxies, since 2001 that reported no exposure to date palm sap or contact with another case. Cases were 4.9 (95% 3.2-7.7) times more likely to consume raw date palm sap and 7.3 (95% 4.0-13.4) times more likely to have contact with a Nipah case than controls. In-depth interviews revealed that 39/182 (21%) of Nipah cases reporting neither date palm sap consumption nor contact with another case were misclassified. Prevention efforts should be focused on interventions to interrupt transmission through date palm sap consumption and person-to-person contact. Furthermore, pooling outbreak investigation data is a good method for assessing rare exposures.
Nipah virus is a paramyxovirus, first discovered during an outbreak involving pigs and people in Malaysia and Singapore in 1998 . Pig consumption of dropped, bat-bitten fruit was the most likely transmission pathway from bats to pigs, and contact with ill pigs was the source of infection for most human cases Chua 2010a; Chua 2010b; Luby and Gurley 2012) . The pigs in the outbreak were culled to prevent further transmission and no human cases have been reported from Malaysia or Singapore since (Chua 2010b; Shapshak et al. 2015) .
In Bangladesh, however, recurrent outbreaks of Nipah virus have been identified almost yearly since 2001 (ICDDRB 2003; ICDDRB 2005; Luby et al. 2009b; Luby and Gurley 2012; Sazzad et al. 2013; Luby 2013) . Case fatality among patients in Bangladesh has been >70% (Luby et al. 2009b; Luby and Gurley 2012) with one-third of survivors suffering permanent neurological deficits (Sejvar et al. 2007 ). The most common risk factors for human disease in Bangladesh have been consumption of raw date palm sap, a national delicacy (Luby et al. 2006; Montgomery et al. 2008; Luby et al. 2009b; Rahman et al. 2012) , and contact with another human Nipah case, primarily through caregiving (Hsu et al. 2004; Gurley et al. 2007b; Gurley et al. 2007a; Luby et al. 2009a; Homaira et al. 2010) . Some cases have neither reported exposure to date palm sap nor contact with another Nipah case during epidemiologic investigations, however, suggesting other possible unidentified pathways of transmission (Luby et al. 2006; Luby et al. 2009b; Luby et al. 2009a ).
Identifying alternative routes of Nipah transmission would provide public health officials additional opportunities for prevention and further our understanding of drivers of crossspecies transmission of this highly fatal emerging pathogen. During each individual outbreak the number of cases is small, typically less than 20, which limits the statistical power to identify associations with less common exposures (Luby et al. 2009b ). Therefore, this study aimed to identify rare risk factors for Nipah virus disease by combining all outbreak data from 2004 to 2012. We also aimed to identify other possible exposures to Pteropus fruit bats, the reservoir hosts for Nipah virus, or sick animals that were not adequately addressed in epidemiologic investigations but may have contributed to transmission ( Fig. 1) (Halpin et al. 2011 ).
Two kinds of Nipah cases were included in our analysis: laboratory-confirmed and probable cases. Laboratory-confirmed cases tested IgM positive against Nipah virus in serum via an enzyme-linked immunosorbent assay (ELISA) (Kashiwazaki et al. 2004; Gurley et al. 2007a) . Probable cases presented with febrile illness and altered mental status and had an epidemiologic link to a laboratory-confirmed case, but died before a blood specimen could be collected or before antibodies may have developed (Ramasundram et al. 2000) . Trained interviewers used structured questionnaires to collect basic demographics, household characteristics, travel histories, exposures related to animal, fruit, and contact with a Nipah case. Questionnaires were administered to cases who survived without neurological impairment or, for most cases, two to three proxy respondents. Proxy respondents were selected from within the case household or from others who had daily interactions with the case. Three, sometimes four, geographically matched controls per case were also interviewed. Consumption of raw date palm sap was added to the questionnaire in 2004 and fermented date palm sap, locally known as tari, was added in 2012 (Islam 2011) . Prior to 2004, questionnaires varied and controls were matched for age, sex, and geography. As a result, we restricted our analysis to data collected since 2004.
From 2004 to 2012, epidemiologists identified 157 human Nipah cases, including 80 confirmed and 77 probable cases, and enrolled 632 controls. We combined these data and calculated summary statistics of frequencies and means for all possible exposures. We recoded "don't know" responses as "no" to not reduce the sample size and delete observations with any "don't know" entry; the proportion of such responses did not exceed 10% for any exposure. This recoding would result in an underestimate in the effect of any one exposure. We used conditional logistic regression to build univariable models to investigate associations between each exposure and being a Nipah case, after adjusting for age and sex. We adjusted P values for multiple comparisons using the Bonferroni correction method (Holm 1979) . We capped adjusted P values at 1, so that if an adjusted P value was greater than 1, we only reported 1. We included contact with another Nipah case and history of consuming raw date palm sap in all possible models as we were interested in knowing if other exposures were independently associated with Nipah infection. A multivariable model of best fit was then built using the Akaike Information Criterion and biological plausibility, while still adjusting for age and sex. We calculated odds ratios (ORs) and 95% confidence intervals (CIs), qualifying an association as significant based on a P value < 0.05. Since our probable case definition did not include laboratory confirmation, we conducted a sensitivity analysis using only lab-confirmed Nipah cases.
For significant factors in the model of best fit, we also calculated the attributable risk using the multivariable odds ratio from the model and the percentage of exposed among cases, using the Bruzzi et al. method (Bruzzi et al. 1985; Rockhill et al. 1998) . Though calculating the attributable risk from an odds ratio is problematic because odds ratios are not collapsible (Greenland et al. 1999b) , because the disease of interest is rare, the odds ratio is a good approximation of the relative risk of Nipah virus. Though the risk of Nipah virus is geographically variable and stronger in the areas where cases and controls are selected, we assume the disease of interest is rare as few Nipah cases occur in the country of Bangladesh (Cummings 2009 ).
To further explore possible rare exposures to Nipah, in 2013 we conducted in-depth interviews with 34 Nipah cases identified since 2004 who themselves or their proxies reported neither exposure to date palm sap nor a Nipah case on the case-control study structured questionnaire. Date palm sap consumption was not investigated as an exposure during 2001-2003, so for the 11 cases during that time period who did not report person-toperson exposure, we sought to retrospectively investigate both date palm sap and other possible bat exposures. Of the total 45 cases investigated, trained anthropologists interviewed 11 cases directly and at least three proxy respondents per case for the 34 cases who did not survive. The interviewers collected detailed field notes on general practices of date palm sap and tari consumption, including whether cases engaged in these exposures prior to illness, as these exposures were not investigated in the first questionnaire. Dropped fruit consumption, typical animal contact, contact with a Nipah case, occupation, and travel were additionally asked about, both as general practices and exposures prior to illness. We reviewed the field notes and calculated the proportion of cases with plausible exposures to bats, bat secretions or excretions based on interview responses.
The Bangladesh Ministry of Health and Family Welfare authorized and actively participated in all reviewed and approved all protocols. Participation in these studies was strictly voluntary. Informed consent was obtained from all individual participants included in the study. For those <18 years of age, individual and parental consent were obtained. icddr,b's Institutional Review Board reviewed and approved the protocol for the in-depth interviews.
From 2004 to 2012, the median age of Nipah cases was 25 years with a range from 6 months to 75 years; 39% of cases were female. The primary occupation was student or homemaker; the mortality rate was 77%. In the univariable case-control analysis, Nipah cases were significantly more likely than controls to live in a tin-shed house than in a thatched house, to visit an area outside their district or a hospital in the one month prior to becoming ill, and to see bats in trees around the house at night (Table 1) . Contact with domestic animals, tree climbing, or other hypothesized risk factors were not significantly associated with being a Nipah case. However, cases were significantly more likely to consume raw date palm sap and have contact with another Nipah case. Controls were significantly more likely to eat washed fruit, either picked from the ground or not, before consuming it.
In the multivariable model that best explained the occurrence of Nipah infection, after controlling for age, sex, person-to-person contact, and date palm sap consumption, those infected with Nipah were significantly more likely to report seeing bats in trees around the house at night (OR 3.3, 95% CI 1.2-8.8) and to have visited an area outside one's own subdistrict (OR 3.5, 95% CI 1.5-8.2) in the month prior to becoming ill than healthy controls ( Table 2 ). The population attributable fraction for raw date palm sap consumption was 50%; that is, half of the Nipah risk in the population could be eliminated if exposure to raw date palm sap were prevented. Similarly, 30% of Nipah cases could be prevented if exposure to another Nipah case were prevented. The population attributable fraction of Nipah infection was 40% for reporting bats feeding in trees around the house at night and 22% for traveling to another sub-district.
During the in-depth interviews, 29 (85%) of the 34 cases or their proxies since 2004 with no reported exposure to date palm sap or a Nipah case on the structured questionnaire reported being exposed to date palm sap or a Nipah case. Of the 34 cases, 12 cases or proxies reported having consumed raw date palm sap while four consumed tari, which was not asked about in questionnaires as an exposure until 2012. No cases or proxies reported consuming both tari and raw data palm sap. Thirteen cases reported exposure to a Nipah case within 5-11 days before their own illness onset. One example of a misclassified case exposure was a patient from the 2004 outbreak in Rajbari, Bangladesh. This patient was breastfed by his mother while the mother was ill with Nipah virus. Although the epidemiologic questionnaire stated that the child had no contact with a febrile patient, upon in-depth interview family members reported that the mother had contact with her 24-month-old child at least during the early stage of her illness. In the 2011 Rangpur outbreak, for a case who was a professional driver, there was no reported contact with a Nipah patient on the questionnaire. The case-patient's colleagues and wife did not report any exposure to date palm sap or tari either; however, did report that the case transported a sick person in his taxi, possibly a Nipah case, during the 15 days before his death.
Five cases from 2004 to 2012 reported not having exposures to date palm sap or a Nipah case during the in-depth interview. One case, a 13-year-old girl, was reported by proxies to regularly eat bat-dropped fruit from the ground as well as climb trees to fetch the fruit. Relatives of another case reported that he played with bat-dropped fruit, such as boroi, olive, or banana. Neighbors of a 65-year-old woman reported that she may have consumed raw date palm sap available at the local market and that she regularly ate dropped star fruit from the grounds of her neighbor's house. The fourth case, an electrician, would climb electric pillars and trees to repair cable lines and may have had direct contact with bats who rested on cable lines. However, family members also reported he consumed bat-dropped fruit 6-7 days prior to his onset of illness. The last case, a 5-year-old boy, was fed bat meat as a traditional remedy for his breathing problem.
The primary objective of our analysis was to identify rare exposures associated with Nipah virus infection by combining data from multiple outbreaks to increase our power to find associations. What we found, however, was that the exposures repeatedly identified in individual outbreaks-consumption of raw sap and contact with another case-remained the best explanation for infection, even with greater power to find new associations. This is an important finding for the prevention of Nipah virus in Bangladesh because by ruling out other major contributors of infection, it allows the public health community to focus its efforts on interrupting these two transmission pathways.
While raw date palm sap consumption and contact with an infected person remained the predominant transmission pathways of human Nipah virus infection in Bangladesh, casepatients were more likely than controls to report that Pteropus bats visit trees near their homes at night. This association could represent a proxy for bats contaminating raw date palm sap at night and a bias of case-patients having prior knowledge that the disease is spread by bats. As the cumulative effects of multiple factors are not necessarily additive, the attributable risk of Nipah infection for bats being around the house at night (40%) is likely largely accounted for in the attributable risk for drinking date palm sap (50%) (Benichou 2001) . Additionally, while case-patients were more likely to travel outside their sub-district, this association could represent our inability to capture exposures to date palm sap or a Nipah case that occurred during travels, as village proxies likely did not travel with casepatients and, therefore, could not report on exposures during travel. In support of this, other studies suggest that date palm sap is not available in every village and that availability of the sap is a primary determinant of whether individuals consume it or not (Nahar et al. 2015) . Although we do not know of outbreaks occurring in areas where these people traveled, they may have encountered cases that went unobserved by our surveillance system. Therefore, the attributable risk of Nipah infection for traveling to another sub-district (22%) could largely be accounted for in the attributable risk of having contact with a Nipah case (30%) or drinking date palm sap (50%).
Zoonotic transmission of Nipah from bat-dropped fruit to pigs, and subsequently from pigs to humans, was the most likely transmission route in the 1999 Malaysian outbreak. Although there is anecdotal evidence of one goat to human transmission event in Bangladesh (Luby et al. 2009a ), we found no evidence of an association between exposure to sick animals and being a human case (Chua 2010a; Luby and Gurley 2012) . Although the in-depth interviews identified some potential exposures to bat secretions, like exposures to dropped fruit, it did not identify any common exposures that we have not evaluated in previous studies. Though we cannot rule out contact with sick livestock or bat-dropped fruit as occasional modes of Nipah transmission, this study suggests these exposures have not been important contributors to spillover infections in Bangladesh.
In the univariable analysis, eating washed fruit was significantly protective; however, less so after adjusting for multiple comparisons. This could be due to chance variation, as based on a 95% confidence interval, 1 in 20 significant findings will be spurious (Rothman et al. 2012) . The more likely explanation though may be that due to the high case fatality of Nipah, controls were able to report exposures more completely than cases. While it is possible that sour fruits with high vitamin C content may have a protective effect against viral infection through ascorbic acid, there was no observed association between eating fruit and Nipah infection in the multivariable model of best fit (Atherton et al. 1978 ).
In the case-control studies conducted during outbreaks, there was evidence of misclassification of exposures. Since most Nipah case-patients die, proxy respondents were asked about case-patient exposures. As proxies may have had limited knowledge of casepatients' activities prior to illness, this likely leads to inaccurate reporting and an underestimate of relevant exposures. This differential misclassification of exposures would lead to a bias towards the null, however, because the reported exposures are likely an underestimate. Therefore, the reported associations are likely an underestimate of true associations. Furthermore, we could not control for having a proxy respondent in the analysis since, as is inherent in case-control studies, selection of study participants is related to the outcome. Because the exposures of consuming date palm sap and having contact with a Nipah case are strongly associated with the outcome, which has a high case fatality, these exposures are spuriously linked to having a proxy respondent (Greenland et al. 1999a; Cole et al. 2010; Rothman et al. 2012 ). This limitation is inherent to investigations of diseases with high case fatality and should be considered in data analysis and interpretation.
While reporting bias likely contributed to the misclassification of exposures, there are numerous examples of discrepancies in quantitative and qualitative findings from outbreak investigations that used both methods. In the 2005 Tangail outbreak of Nipah, qualitative indepth investigations reported much higher date palm sap consumption than the quantitative survey (Luby et al. 2006 ). In the Bangladesh 2007 outbreak of illness from consumption of Xanthium stumarium, a plant also known as ghagra shak, qualitative investigators also found higher levels of reported exposure during in-depth interviews than in quantitative surveys .
Qualitative in-depth interviews are open-ended discussions between study participants and the investigators that aim to build rapport. By spending more time in case villages and building rapport, participants may be more willing to mention certain exposures that are socially stigmatizing or that participants suspect may be a cause of the outbreak. Interviews that use structured questionnaires are less likely to elicit honest responses of sensitive information, as they are relatively shorter by comparison and less likely to build rapport (Steckler et al. 1992) . Qualitative interviews are also imperfect, however; they can result in confirmation bias as a result of prodding by the interviewer, as discussed in anthropologic research literature (Bernard 1984) . Furthermore, when multiple proxy respondents are interviewed, it is important to note how discrepancies between answers will be interpreted. Thus, outbreak investigations could benefit from the use of both approaches, particularly when the hypotheses about exposures causing disease are still being formed.
A limitation of studying Nipah virus and other rare, emerging infections is that our knowledge of relevant exposures changes over time and subsequent changes in questionnaires compromise our ability to compare data across years. This is true for our analysis where date palm sap was not included in the questionnaire until 2004 and tari consumption not until 2012.
Furthermore, as some case-patients die before a blood sample can be collected, cases are classified as either probable, having an epidemiologic link with another case, fever, but no lab-confirmation, as well as lab-confirmed. An epidemiologic link can signify contact, or living in the same village, or sharing date palm sap. A concern with the use of probable cases in our analysis is that probable cases are not laboratory confirmed. In a sensitivity analysis conducted with only lab-confirmed Nipah cases, we find similar results and that our model is still significant with slightly reduced effect estimates for drinking raw date palm sap and larger effect estimates for having contact with another case (Supplement Table 1 ).
Additionally, cases were matched to controls initially for geographic location, age, and sex.
Only geographic location was used as a matching variable in later years to avoid complications of overmatching. To control for this change in methods, we adjusted for age and sex in all models. This likely resulted in a loss of statistical efficiency as some cases were matched for age and sex in addition to controlling for age and sex in the models. Lastly, imperfect recall is also a concern because in-depth interviews in this study were conducted 1-10 years after the exposures occurred, which may have further influenced the misclassification of exposures.
While this analysis did not identify new risk factors for Nipah infection, it confirmed the centrality of transmission through contaminated sap and contact with infected patients, allowing public health officials to further focus efforts to interrupt these transmission pathways. Our study also highlights the benefits of using both qualitative and quantitative methods in outbreak investigations; using the strengths of each approach is advantageous particularly when knowledge of the transmission of a disease is still evolving.
Efforts to prevent this highly fatal emerging infection should focus on interrupting transmission from bats to humans through date palm sap and reducing the exposure to Nipah-infected patients in Bangladesh. The use of bamboo skirts around the date palm sap collection area on trees prevents bat contamination of sap and can be practiced more widely, thereby reducing population risk Nahar et al. 2014) . For Nipah infection spread, handwashing with soap has had a significantly protective effect (Gurley et al. 2007a) and interventions to promote handwashing and increase access to handwashing stations in hospitals are urgently needed. Though some infection control has been ongoing, limited infection control within hospitals and by family caregivers continue to be major contributors to disease transmission (Ram et al. 2011; Islam et al. 2014) . The recent West Africa Ebola outbreak made evident the difficulty of infection control in limited resource settings and the need for interventions to prevent person-to-person transmission, especially for similarly fatal, emerging infections (Chertow et al. 2014) . Importantly, improving such interventions could cause additional benefits by reducing other bat-borne pathogens, including emerging and endemic pathogens, which could transmit via contaminated fluids or from person-toperson.
Refer to Web version on PubMed Central for supplementary material.
|