key: cord-276254-q04hqra2
authors: Paul, Kishor Kumar; Salje, Henrik; Rahman, Muhammad W.; Rahman, Mahmudur; Gurley, Emily S.
title: Comparing insights from clinic-based versus community-based outbreak investigations: a case study of chikungunya in Bangladesh
date: 2020-06-02
journal: Int J Infect Dis
DOI: 10.1016/j.ijid.2020.05.111
sha: 
doc_id: 276254
cord_uid: q04hqra2

Abstract Background Outbreak investigations typically focus their efforts on identifying cases that present at healthcare facilities. However, these cases rarely represent all cases in the wider community. In this context, community-based investigations may provide additional insight into key risk factors for infection, however, the benefits of these more laborious data collection strategies remains unclear. Methods We used different subsets of the data from a comprehensive outbreak investigation to compare the inferences we make in alternative investigation strategies. Results The outbreak investigation team interviewed 1,933 individuals from 460 homes. 364 (18%) of individuals had symptoms consistent with chikungunya. A theoretical clinic-based study would have identified 26% of the cases. Adding in community-based cases provided an overall estimate of the attack rate in the community. Comparison with controls from the same household revealed that those with at least secondary education had a reduced risk. Finally, enrolling residents from households across the community allowed us to characterize spatial heterogeneity of risk and identify the type of clothing usually worn and travel history as risk factors. This also revealed that household-level use of mosquito control was not associated with infection. Conclusions These findings highlight that while clinic-based studies may be easier to conduct, they only provide limited insight into the burden and risk factors for disease. Enrolling people who escaped from infection, both in the household and in the community allows a step change in our understanding of the spread of a pathogen and maximizes opportunities for control.

Infectious disease outbreaks have the potential to place a significant burden on public health resources. Understanding who is at risk of becoming infected is critical for the focused targeting of interventions. Due to relative ease of access and limited cost requirements, outbreak investigations typically focus on cases that present at formal healthcare centers such as hospitals or community clinics. For example, data collection performed as part of epidemiological investigations during the recent epidemics of Ebola, Zika and MERS focused on quantifying the number of cases and their characteristics. (Al-Abdallat et al., 2014 , Lu et al., 2015 , Teixeira et al., 2016 These case-counting exercises provided key insights into fundamental epidemiological parameters such as the basic reproductive number and case fatality rates, and allowed the projection of the future course of the epidemic. (Aylward et al., 2014 , Lessler et al., 2014 , Lewnard et al., 2014 , Yamin et al., 2015 However, without information on the underlying population, and especially characteristics of individuals who avoid infection, these approaches limit our ability to make mechanistic insights, quantify burden of disease, and identify risk factors for infection, hampering efforts to develop targeted control strategies.

Cases that present at healthcare centers may only represent a small minority of all cases. In addition, some individuals are more likely to visit formal healthcare providers than others, including those with more severe illness, and differences in healthcare seeking can vary by age, gender and socioeconomic status. (Chowdhury et al., 2007 , Nikolay et al., 2017 , Pandey et al., 2002 Household-based outbreak investigations, where investigation teams visit affected J o u r n a l P r e -p r o o f communities, permit a more comprehensive understanding of pathogen spread that limits the impact of healthcare seeking patterns. (France et al., 2010) However, these investigations are usually still focused on identifying individuals that got sick. (Boore et al., 2013 , France et al., 2010 Without also understanding who is avoiding infection in a community, it is difficult to identify the key risk factors for infection, limiting potential inferences. The possible insights from alternative investigation strategies have not previously been systematically compared. Here, we use the results of a detailed chikungunya outbreak investigation from Bangladesh as an example to consider the inferences made under different investigation scenarios.

Chikungunya virus is a mosquito-borne alphavirus transmitted to humans by Aedes mosquitoes causing acute fever, joint pain, and skin rash. (Aubry et al., 2015) Chikungunya fever was first recognized in 1952 in Tanzania. (Lumsden, 1955) Since then, outbreaks of chikungunya have been regularly identified across the tropics and sub-tropics. The first chikungunya outbreak in Bangladesh was identified in 2008 in two northwestern districts bordering India.(ICDDR) Since then regular outbreaks have been detected. (Khatun et al., 2015 , Salje et al., 2016b Here we use the results from a detailed investigation of an outbreak of chikungunya virus in a village in Tangail, Bangladesh where the outbreak team visited every household in the community and interviewed all members in each household. The comprehensive household investigation captured both those who did get infected and those that escaped from infection. The objective of this study was to compare our approach, in terms of the inferences about the outbreak, to more limited investigation strategies.

In late November 2012, a local health official of Gopalpur sub-district in Tangail households in the village consented to being enrolled in the study. Questionnaires were administered in all households to identify suspected cases, identify demographic characteristics, and travel histories of individuals within households. Suspect cases were defined as residents with acute onset of fever with rash or joint pain within 6 months prior to beginning the investigation.

Study staff administered questionnaires to household heads about household members' demographic data and history of illness, water source, construction materials, and mosquito control measures in the household. Potential mosquito breeding containers in and around the participating households with stored water were inspected for presence of larvae. Suspected cases were asked about their symptoms with onset date and specifics about their treatment seeking behavior. The GPS location of all homes was also recorded.

Determining the etiology of the outbreak All household members, irrespective of their suspected case status, were asked to provide a single 5 ml blood specimen for laboratory testing. Blood specimens were spun in the field to separate serum, which were then stored on ice and transported to the virology laboratory of IEDCR. The serum samples were tested for IgM antibodies against chikungunya by enzyme linked immunosorbant assay (ELISA) (Standard Diagnostics, Inc., South Korea). Suspected casepatients who had IgM antibodies against chikungunya in their serum were termed laboratory confirmed cases.

We created four different datasets that allowed us to consider different outbreak investigation strategies:

This dataset consisted of all suspect cases that reported that they visited a formal healthcare setting (defined as government or non-government primary healthcare center/clinic/hospital) following the onset of symptoms.

This dataset consisted of all suspect cases, irrespective of their healthcare seeking behaviors.

This dataset consisted of all suspect cases plus controls consisting of household members of these cases.

This dataset consisted of all members of all households in the community, regardless of symptoms.

The epidemic curve was constructed using symptom onset date of chikungunya cases. GPS locations of households with and without chikungunya cases were used to prepare spatial distribution maps. For the case-only datasets (datasets A and B), we compared the age and sex distribution of the cases with that for the district from the 2011 census (Bangladesh Bureau of Statistics, 2011) . For the datasets with information on individuals who escaped from infection (datasets C and D), we initially used simple logistic regression to compare the demographics, typical apparel worn, travel history within the last six months, and household characteristics of cases with non-cases. We then built multivariable logistic regression models to identify adjusted risk factors for chikungunya fever. We initially placed all variables with a p-value of <0.05 in the unadjusted analysis into a multivariable model. We then used backward stepwise selection using the Akaike information criterion (AIC) (Sauerbrei et al., 2007) to identify the best model.

Not all individuals who get infected will present with symptoms. We attempted to capture these individuals by asking for blood samples from all community members. To assess the impact of J o u r n a l P r e -p r o o f misclassifying asymptomatically infected individuals as controls in datasets C and D, we conducted sensitivity analyses where these individuals were reclassified as cases.

All participants provided written informed consent prior to interviews and blood specimen collection and the Ministry of Health and Family Welfare, Government of Bangladesh reviewed and approved the outbreak investigation plan.

The 

Ninety-five suspect cases reported visiting a formal healthcare facility for symptoms consistent with chikungunya between July and November, with the peak number of cases occurring in J o u r n a l P r e -p r o o f October (Figure 1 : Panel A). Cases sought care in three different centers: 25 sought care in a government run community clinic, 8 in a government run sub-district health complex and 62 in a private clinic. The median age was 28 years (interquartile range (IQR) = 14-45 years) and the majority (60%) were female (Table 1 ). If we used the age and sex distribution of the district from the 2011 national census, we find that there is an increased risk of disease in those between the ages of 30-49 compared to those aged below 10 years (OR 2.51, 95% 1.24-5.51) and that females were at increased risk of infection compared to males (OR of 1.43, 95% CI: 0.94-2.23) ( Table 2) .

An additional 269 suspect cases were identified in the community who did not seek care in formal healthcare facilities. Of these, 246 individuals visited a local pharmacy and 21 individuals visited the informal sector (unlicensed medical practitioner, traditional healer, and homeopath).

The distribution of dates of symptom onset for all cases was nearly identical to the distribution for those that visited clinics (Spearman correlation of 0.95) ( Figure 1A ). The proportion of suspect cases visiting a clinic varied between 18% in 51-55 years age group and 44% in ≤5 years age group ( Figure 1B) . The conclusions about age, sex, educational levels, use of mosquito controls and clinical presentation of suspect cases were similar when using datasets of all cases or only those that sought care in clinics (Table 1) , however, those who presented to clinics were more likely to travel outside the district (41% vs 28%, p-value 0.001). Cases who attended formal healthcare settings also appeared to come from similar parts of the community as cases who did not (Figure 2A-B) . Similar to the analysis using clinical cases only, using data from the national census identified increasing risk among females for being a case (Table 2) .

Inference from community cases plus controls from same household J o u r n a l P r e -p r o o f Incorporating controls from the households where cases reside allowed us to assess additional potential risk factors for being a case. Consistent with inferences using census data, logistic regression models that used household controls also identified increased risk among females (aOR 1.75, 95% CI 1.31-2.36) ( Table 2 ). In addition, this analysis showed that cases were significantly less likely to have secondary (aOR 0.64, 95% CI 0.42-0.98) or more formal (higher secondary) education (aOR 0.28, 95% CI 0.17-0.47) compared to the household controls.

Incorporating data from the entire community showed that the chikungunya outbreak was largely constrained to the center of the village, with few households affected on the east and west borders but virtually all households affected in the center ( Figure 2C ). This is despite the entire community only being a few hundred meters wide. The expanded dataset also allowed us to understand the risk factors for infection in the wider community. As with the previous analyses, females had an increased risk of being a case (OR: 1.62, 95% CI 1.29-2.03) ( Table 2 ), although the difference by sex was concentrated in adults with no difference among children ( Figure 1C ).

Further individuals who reported usually wearing clothing that exposed both limbs had 1.80 the odds of being a case compared to individuals wearing clothing that exposed upper limbs only (95% CI 1.30-2.50). Those who had travelled outside Tangail district within the last six months also had increased odds of being a case (aOR 1.47, 95% CI 1.06-2.03). Individuals who had higher secondary or more formal education (aOR 0.38, 95% CI 0.24-0.62) were less likely to be a case than individuals without formal education. We did not identify any household characteristics that were associated with being a case, including presence of mosquito larvae (aOR 0.87, 95% CI: 0.66-1.16), daily use of anti-mosquito coil (aOR 1.03, 95% CI: 0.82-1.31), number of J o u r n a l P r e -p r o o f household members (aOR 1.00, 95% CI: 0.79-1.26), and number of rooms in the household (aOR 1.11, 95% CI: 0.82-1.51).

Fifty-two individuals without symptoms tested positive for CHIKV. We found no significant demographic differences between symptomatic suspected cases and IgM-confirmed asymptomatic cases in those who gave blood (Table S1 ). In sensitivity analysis, we removed these individuals from the 'control' population and included them in the 'case' populations. Risk factors for being a case identified in the previous analysis remained similar in both scenarios where we considered household contacts as controls and individuals from all community households as controls (supplementary information, Table S4 ). However, we found important differences in the probability of providing blood. Those with symptoms were 6.1 times more likely to provide blood than those without symptoms. Further, among asymptomatic individuals, only 2% of children 0-9 years provided a sample compared to 15% among those 30-49 (Table   S2 ). There were also significant differences by sex (14% of asymptomatic males gave blood compared to 8% of females, p-value <0.001) and educational level with more educated people less likely to provide samples (Table S2 ).

Outbreak investigations are central to informed responses to public health emergencies caused by the emergence of an infectious pathogen. However, outbreak investigations currently largely revolve around case-counting exercises that limit our ability to identify who is at risk for J o u r n a l P r e -p r o o f infection and who is not. Here, by using the results of a comprehensive outbreak investigation, we have been able to explicitly explore the impact of different investigation strategies in the same outbreak. We found that a clinic-based study that used data from all the formal healthcare settings would have identified a quarter of all cases and, using census data, have correctly identified female sex as an important risk factor for disease. However, it is only through the recruitment of people who did not get sick that we could identify the importance of travel history, educational level and apparel usage in determining who gets sick. Controls from the wider community were also required to demonstrate which household-level characteristics were important for risk,

showing that the use of mosquito coils was not protective, and to map spatial heterogeneity in risk, key to intervention development and deployment.

This study highlights the significant heterogeneity in healthcare seeking. Even in a small community such as this, cases visited nine different sources of healthcare, three of which could be considered formal healthcare settings. Infectious disease surveillance activities are unlikely to be able to collate datasets from this diverse range of healthcare sources, even among only those within formal sector, suggesting that outbreak investigations that rely on cases that seek healthcare likely substantially underestimate the magnitude of outbreaks.

Using the results of our study, we provide our assessment of the ability of different investigation strategies to capture key characteristics of an outbreak (Table 3 ). In practice, the decision to expand outbreak investigations beyond information available from healthcare systems will depend on the resources available. Where outbreak teams are already performing communitybased case-investigations, the additional time and effort to also collect data on those without symptomsboth from case-households as well as neighboring householdsmay be marginal. This comprehensive outbreak investigation employed ten field-based investigators and took seven days to complete. An investigation strategy only focused on cases in the community would have taken only marginally less person-time as finding cases in the community anyway typically requires comprehensive door-to-door surveys. Our findings highlight how this additional data collection effort can help reveal the drivers of transmission, allowing mechanistic insight into pathogen spread and maximizing opportunities to control, many of which would not be possible from case-based investigations (Table 3 )..Where it is collected, an additional major benefit of the comprehensive dataset is that it can inform mathematical models that reconstruct entire outbreaks, allowing us to estimate the mean transmission distance (previously estimated here at 95 meters) (Salje et al., 2016b) .

Travelling outside Tangail district within the six months before the outbreak was associated with increased chikungunya fever risk. Human movement can introduce CHIKV into new areas, causing epidemics (Chretien and Linthicum, 2007) . No other areas of Bangladesh were reporting outbreaks of CHIKV at this time, though outbreaks may have been missed due to poor surveillance. Although individuals of all ages were affected by chikungunya in this outbreak, incidence increased with age among females, potentially linked to increased time women spend at home compared to males, increasing their risk of being bitten by the largely home-dwelling Aedes mosquito (Salje et al., 2016b) . In this outbreak, household use of mosquito coils was not protective against chikungunya, which is consistent with the findings from a recent meta-analysis on household level risk factors for dengue, which is also spread through Aedes mosquitoes (Bowman et al., 2016) .

Serum samples have the potential to provide important information about the level of asymptomatic infection during an outbreak, as has previously been shown during previous CHIKV outbreaks (Salje et al., 2016a , Sissoko et al., 2008 . In addition, this outbreak investigation was carried out six months after the outbreak began and community members may have been unable to reliably recall their symptoms or the date their symptoms started, particularly for milder illnesses, which may have led to an underestimation of suspected cases. Serological confirmation could help detect any missing infections. However, our study highlights how some caution needs to be taken when interpreting serological data. Firstly, while we sought to obtain blood samples from all participants, only one in five individuals agreed. We found that the probability of agreeing to provide blood depended strongly on having had chikungunya symptoms (individuals who had symptoms were more likely to provide samples). Children, women and those with a high educational level were less likely to give blood. Secondly, the sensitivity of the commercial assay we used has been estimated to be <40% in individuals where IgM is still circulating (Johnson et al., 2016) and is likely to be even lower here, as the blood draw occurred after IgM antibodies would have waned to undetectable levels for many infected individuals (Kam et al., 2012) . Future studies should consider underlying biases in who is providing blood as well as considering the use of complementary IgG assays to help improve the interpretability of serological findings.

This investigation suggests that chikungunya virus has become an emerging public health problem in Bangladesh, and outbreak investigations of emerging infections often have the objective of estimating attack rates of diseases and identifying the risk factors that lead to infection. Our analysis suggests that the optimal strategy for attaining these objectives during an outbreak is to conduct case finding, testing, and data collection in communities. Many recent J o u r n a l P r e -p r o o f outbreaks of emerging infections have suffered due to a lack of detailed information about attack rates and risk for infection, due to their limited investigation strategies (Ahmed et al., 2015 , Ballera et al., 2015 , Khatun et al., 2015 . Future investigations of emerging infection outbreaks should consider using these more intensive strategies, at least in a subset of investigations, to improve our understanding of these infections and our public health response.

According to institutional data policy of the International Centre for Diarrhoeal Disease Research, Bangladesh (icddr,b), summary of data can be publicly displayed or can be made publicly accessible.

To protect intellectual property rights of primary data, icddr,b cannot make primary data publicly available. However, upon request, Institutional Data Access Committee of icddr,b can provide access to primary data to any individual, upon reviewing the nature and potential use of the data. Requests for data can be forwarded to:

This work was supported by Centers for Disease Control and Prevention (CDC), Atlanta, USA [cooperative agreement no: 5U01CI000628]. In addition, the Government of Bangladesh, Canada, Sweden and the UK provided core/unrestricted funding support for this work. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. 

Chikungunya Virus Outbreak

Hospital-associated outbreak of Middle East respiratory syndrome coronavirus: a serologic, epidemiologic, and clinical description

Chikungunya outbreak

Ebola virus disease in West Africa--the first 9 months of the epidemic and forward projections

Investigation of chikungunya fever outbreak in Laguna

Bangladesh population and housing census 2011

Added value of a household-level study during an outbreak investigation of Salmonella serotype Saintpaul infections

Is dengue vector control deficient in effectiveness or evidence?: Systematic review and meta-analysis

Delivery complications and healthcare-seeking behaviour: the Bangladesh Demographic Health Survey

Chikungunya in Europe: what's next?

Household transmission of 2009 influenza A (H1N1) virus after a school-based outbreak

First identified outbreak of chikungunya in Bangladesh

Laboratory Diagnosis of Chikungunya Virus Infections and Commercial Sources for Diagnostic Assays

Early appearance of neutralizing immunoglobulin G3 antibodies is associated with chikungunya virus clearance and long-term clinical protection

An Outbreak of Chikungunya in Rural Bangladesh

Estimating Potential Incidence of MERS-CoV Associated with Hajj Pilgrims to Saudi Arabia

Dynamics and control of Ebola virus transmission in Montserrado, Liberia: a mathematical modelling analysis

Ebola Virus Outbreak Investigation

An epidemic of virus disease in Southern Province, Tanganyika Territory

General description and epidemiology

Evaluating Hospital-Based Surveillance for Outbreak Detection in Bangladesh: Analysis of Healthcare Utilization Data

Gender differences in healthcare-seeking during common illnesses in a rural community of West Bengal

Reconstruction of 60 Years of Chikungunya Epidemiology in the Philippines Demonstrates Episodic and Focal Transmission

How social structures, space, and behaviors shape the spread of infectious diseases using chikungunya as a case study

Selection of important variables and determination of functional form for continuous predictors in multivariable model building

Seroprevalence and risk factors of chikungunya virus infection in Mayotte

The Epidemic of Zika Virus-Related Microcephaly in Brazil: Detection, Control, Etiology, and Future Scenarios

EFfect of ebola progression on transmission and control in liberia

The authors have no competing interests to declare. J o u r n a l P r e -p r o o f
