Utility of mosquito surveillance data for spatial prioritization of vector control against dengue viruses in three Brazilian cities
- Kim M Pepin^{1, 2, 3}Email author,
- Clint B Leach^{3},
- Cecilia Marques-Toledo^{4},
- Karla H Laass^{5},
- Kelly S Paixao^{5},
- Angela D Luis^{1, 3, 8},
- David TS Hayman^{3, 6, 10},
- Nels G Johnson^{3},
- Michael G Buhnerkempe^{3, 9},
- Scott Carver^{3, 7},
- Daniel A Grear^{3},
- Kimberly Tsao^{3},
- Alvaro E Eiras^{5}Email author and
- Colleen T Webb^{1, 3}
https://doi.org/10.1186/s13071-015-0659-y
© Pepin et al.; licensee BioMed Central. 2015
Received: 17 December 2014
Accepted: 12 January 2015
Published: 15 February 2015
Abstract
Background
Vector control remains the primary defense against dengue fever. Its success relies on the assumption that vector density is related to disease transmission. Two operational issues include the amount by which mosquito density should be reduced to minimize transmission and the spatio-temporal allotment of resources needed to reduce mosquito density in a cost-effective manner. Recently, a novel technology, MI-Dengue, was implemented city-wide in several Brazilian cities to provide real-time mosquito surveillance data for spatial prioritization of vector control resources. We sought to understand the role of city-wide mosquito density data in predicting disease incidence in order to provide guidance for prioritization of vector control work.
Methods
We used hierarchical Bayesian regression modeling to examine the role of city-wide vector surveillance data in predicting human cases of dengue fever in space and time. We used four years of weekly surveillance data from Vitoria city, Brazil, to identify the best model structure. We tested effects of vector density, lagged case data and spatial connectivity. We investigated the generality of the best model using an additional year of data from Vitoria and two years of data from other Brazilian cities: Governador Valadares and Sete Lagoas.
Results
We found that city-wide, neighborhood-level averages of household vector density were a poor predictor of dengue-fever cases in the absence of accounting for interactions with human cases. Effects of city-wide spatial patterns were stronger than within-neighborhood or nearest-neighborhood effects. Readily available proxies of spatial relationships between human cases, such as economic status, population density or between-neighborhood roadway distance, did not explain spatial patterns in cases better than unweighted global effects.
Conclusions
For spatial prioritization of vector controls, city-wide spatial effects should be given more weight than within-neighborhood or nearest-neighborhood connections, in order to minimize city-wide cases of dengue fever. More research is needed to determine which data could best inform city-wide connectivity. Once these data become available, MI-dengue may be even more effective if vector control is spatially prioritized by considering city-wide connectivity between cases together with information on the location of mosquito density and infected mosquitos.
Keywords
Vector control Dengue Surveillance Vector density Mosquito-human interactions Gravity model INLABackground
Understanding the relationship between Aedes aegypti vectors and the patterns of dengue fever they cause is important in the design of vector-based disease control strategies. Because it is often not feasible or possible to eradicate the mosquito vectors [1], quantitative knowledge of how vector density relates to disease incidence is essential for deciding how much vector populations need to be reduced in order to decrease disease incidence adequately. Mechanistic knowledge of transmission is also important because methods of vector control that are designed based on perceived spatial patterns of cases are often not effective [1,2]. Identifying how vectors are connected to disease incidence in space and time would allow for more cost-effective strategies of implementing vector controls.
The strength and direction of the relationship between mosquito density and dengue infection varies depending on the spatial scale at which data are collected and community characteristics [3-8]. For example, when comparing adult vector densities with prevalence of human infections across three sets of community conditions (urban, suburban, slum) within Rio de Janeiro, Brazil, Honorio et al. [3] found higher infection prevalence in the slum where vector density was lowest. This negative relationship was hypothesized to be because living conditions in the slum facilitated greater rates of vector-human contact relative to the highly developed urban area. At the household scale, no relationship between vector density and disease prevalence was found [3], although it was acknowledged that larger numbers of infections are required at this scale before appropriate conclusions can be drawn. On the other hand, in rural villages in Thailand, a non-significant but positive trend in the relationship between adult vector density and child infection prevalence was found at the household and between-house levels [7]. Considering that within- and between-house transmission have been shown to be important [9], the weak relationship between adult vector density and human infections at the household level is surprising.
One potential explanation for the weak relationship is sampling – the number of replicate samples in space and time, or techniques used for vector collection, may not be adequate for estimating household mosquito density at a level of precision that is smaller than the ecologically-determined variation in vector density. A second explanation could be human movement [10] – human contact patterns at different spatial scales (local and long-distance) can explain spatial dengue transmission [9,11], highlighting that movement at multiple spatial scales is important to consider when linking vector densities to human cases. Theoretical work has demonstrated that the rate of within-city transmission of dengue virus depends on the type of human movements: regular movement patterns due to commuting patterns, for example, can slow the rate of disease spread by up to 25% in comparison to random movement patterns [12]. In contrast, temporally unstructured movements, such as those found in resource-poor settings, can increase the size of an epidemic by up to 20% [13]. Thus, consideration of human movement at different spatial scales is important for understanding how mosquito density data can be used for targeting vector controls.
Several cities in Brazil have implemented a city-wide mosquito trapping system, MI-Dengue, which monitors weekly prevalence of gravid Ae. aegypti and Ae. Albopictus city-wide in real time [14-16]. Traps are associated with households and spaced in a grid-like manner at ~200-300 m intervals, depending on the city. Vector density data are automatically available for control personnel who respond by focusing source reduction, larvicide and, more rarely, adulticide activities to neighborhood blocks with high mosquito density. The MI-dengue system – based mainly on the idea that spatially targeting areas with higher densities of gravid female mosquitos will decrease case loads using fewer resources – has been shown to be effective and cost-effective for reducing human infections [16]. It has been demonstrated that confirmed cases in humans cluster with high mosquito density in space and time [14], but rigorous quantitative analyses that identify how to best use the surveillance data have not been conducted. Although information on infected mosquitos and confirmed cases in humans are given the highest weight in spatial prioritization of vector control, these data are rarer and often not available until well after transmission has occurred, emphasizing the importance of identifying the best method of using mosquito density data in spatial prioritization of vector control.
While experiments to determine appropriate spatial scales for estimating vector density are still ongoing, the available data are numerous (~5,726 - 43,467 mosquitoes surveyed annually per city) and could reveal useful insight on the spatio-temporal relationship between vector densities and human cases within entire cities. Here, we sought to better understand the city-wide relationship between vector densities and human cases to provide further guidance for spatially targeting vector control work. Our analysis has the following four aims, to: 1) quantify the relative role of city-wide mosquito surveillance data in predicting city-wide cases of dengue, 2) identify the spatial scale at which case data from other neighborhoods are important, 3) identify whether readily available data related to urban characteristics can be used to approximate spatial patterns of human cases, and 4) understand how city-wide mosquito surveillance data can be used to spatially prioritize vector control activities in order to have the maximum effect on preventing cases of dengue fever. We base our analyses on data from Vitoria city, Brazil, because it had the longest time series of surveillance data (~5 years), but we use data from two other cities for validation of model structure and a deeper understanding of model parameters.
Methods
Study site
Models were developed using data from Vitoria city, Brazil, an economically prosperous coastal city that is the largest city (348,265 inhabitants) in the state of Espirito Santo in southeastern Brazil. Among the 27 major cities in Brazil, Vitoria has the 4^{th} highest human development index (HDI; 0.85), the highest gross domestic product per capita and an unemployment rate of 7.25% (Brazilian Institute of Statistics and Geography, 2010 Census). The climate is tropical with an annual mean temperature of 23°C and a rainy season between October and January (National Institute of Meteorology, Brazil). Due to its prosperity, size and port capabilities, there is frequent movement of people and merchandise to and from nearby and more distant cities that are less developed.
Data from two other cities, Governador Valadares (GV; population 263,594) and Sete Lagoas (SL; population 208,847), both in the state of Minas Gerais, Brazil, were used for model testing and validation. Both cities have a history of dengue fever outbreaks and are similarly economically prosperous with HDI and unemployment rates of 0.77 and 6.8% (GV), and 0.76 and 6.8% (SL) (Brazilian Institute of Statistics and Geography, 2010 Census). The river Doce bisects GV acting as a gateway between major marine ports. Annual mean temperatures are 24.6°C (GV) and 20.9°C (SL), with a rainy season between October and March (National Institute of Meteorology, Brazil).
Neighborhood-level population sizes, areas (Additional file 1, spreadsheet "Neighborhoods") and economic data were obtained from the 2010 census (mentioned above), from the local vector control managers and the Ministry of Health Secretaries. For Vitoria, economic values were the sum of the registered commercial (including industry and service) units for each neighborhood. For GV and SL, neighborhood economic data were either the number of registered residences or commercial units per neighborhood.
Case data
Notified cases of dengue fever were obtained from each city’s Ministry of Health Secretary’s official database, which lists dengue-fever cases by their residential address and date of first symptoms. In Brazil, dengue is a mandatory notifiable disease and thus the database represents all cases where any kind of medical care was sought. However, only samples at the start of an epidemic are validated for the presence of dengue virus. Once an epidemic is deemed started, most other cases are diagnosed symptomatically, such that consistent serotype information is unavailable. Although neighborhood assignments were complete, street address information was often lacking, thus we aggregated the case data to the neighborhood level - a political boundary defined by the city. The numbers of neighborhoods in each city were: Vitoria – 75, GV– 65, SL – 98. Neighborhood population sizes and areas were variable both within and between cities (mean ± 2 standard errors for population sizes and areas in km^{2} were: Vitoria – 4,080 ± 1,614, 0.47 ± 0.11; GV – 3,435 ± 862, 3.14 ± 4.69; SL – 2,000 ± 327, 0.37 ± 0.06; Additional file 1: Table S1). We summed the cases in weekly intervals to match the temporal scale of the mosquito data.
Mosquito surveillance data
Mosquito data were obtained from a city-wide surveillance system (MI-Dengue) [15] managed by the company, Ecovec, which originated from an academic setting and is located in Belo Horizonte, MG, Brazil. The system is comprised of a network of sticky traps, called MosquiTRAP, which have been extensively tested and described elsewhere [17-21]. Briefly, traps are placed in a lattice throughout the entire city. Each trap is checked weekly for mosquitos, which are identified to species level. The data are entered by cell phones to a database that automatically generates maps of mosquito density for control personnel, who target control to highly infested areas. We obtained weekly counts of the gravid female Ae. aegypti (93.2– 98.4% of all mosquitos depending on city) and Ae. albopictus species, the primary and secondary vectors of dengue fever. Because each trap was located on the inside or outside of a residence, we expressed the mosquito data as average household mosquito density per neighborhood (mosquitos/traps per neighborhood per week; 18.6 traps/neighborhood on average) to match the spatial scale of the available case data. Using an average household abundance estimate also has the advantage of reducing the uncertainty in household mosquito density compared with using single-replicate trap-level counts for each time point. Mean number of mosquitos and traps counted per week across the three cities were: Vitoria – 716.8 ± 342.5 standard deviation (SD) and 1391.6 ± 32.0 SD, respectively; GV – 212.5 ± 81.8 SD and 373.0 ± 50.4 SD, respectively; SL – 95.4 ± 72.3 SD and 411.2 ± 123.8 SD, respectively (Table S1).The area monitored per city was: 33 km^{2} (Vitoria), 27 km^{2} (GV) and 31 km^{2} (SL), which yields mean weekly trap monitoring densities of 42.2, 13.8 and 13.3 traps per km^{2} in the three cities respectively (Additional file 1, Spreadsheet "Traps").
In all three cities, routine vector control occurs following guidelines of the Brazilian Dengue Control Program. This includes mainly larvacide and source reduction activities that occur systematically (moving from block to block) throughout each city year-round. In addition to these activities, adulticide is conducted in blocks where high numbers of mosquitoes are identified, following the recommendations by Ecovec (www.ecovec.com). The effects of these controls, and other factors that affect mosquito populations such as weather, are implicit in the mosquito density data. Thus, although mosquito populations are altered by several biotic and abiotic factors, the mosquito surveillance data are a means of directly examining effects of mosquito density on disease incidence.
Statistical model structure and parameter estimation
where Y(t) represents mean mosquito density in the entire city during week t. Approximate Bayesian inference by integrated nested Laplace approximations was used for parameter estimation. R software Version 3.0.1 and the package R-INLA (www.r-inla.org) were used to perform the analyses [22].
Description of covariates
The importance and structure of spatial coupling between neighborhoods (a proxy for human movement) was examined as a main effect using a modified gravity model (described below). All covariate data were normalized in order to compare the strength of parameter estimates. A term for spatial autocorrelation was not included in the final models because it was not significant (according to a Moran’s I test on shifted residuals) in preliminary fitting of gravity model terms. We also compared our models, which included a covariate-based exploration of the case data, with autoregressive lag 1 models (AR1) and found similar levels of predictive power (data not shown).
Candidate structures for the components of the f(x _{ j } )
Components of f(x _{ j } ) | Description of hypothesis tested |
---|---|
1) 1 | H_{1}: Dengue cases occur in a random spatial pattern based on the number of cases in all other neighborhoods. |
2) d _{ ij } ^{ α } | H_{2}: The spatial pattern of dengue-case occurrence correlates with distance between other neighborhoods; cases are more likely to occur in neighborhoods that are closer to neighborhoods experiencing cases. |
3) (1/E ^{ * } _{ j })^{ α } | H_{3}: The spatial pattern of dengue-case occurrence correlates with neighborhood economy values; cases are more likely to occur if neighborhoods with high-economy values are experiencing cases. |
4) (1/D _{ j })^{ α } | H_{4}: The spatial pattern of dengue-case occurrence correlates with neighborhood population density values; cases are more likely to occur if neighborhoods with high-density values are experiencing cases. |
Model selection
- 1)
Selection of lags. For each possible covariate (as shown in Equation 1) individually, we identified the best lag time between it and the response variable (x1, x2, y1, y2, z1 and z2 in Equation 1). Lags were calculated as a 3-week average because we hypothesized that a window of time in the past may best explain the relationship (preliminary analyses confirmed this hypothesis). The 3-week window was chosen because 2–3 weeks is the combined amount of time from an infectious mosquito bite to a case report, on average [27]. This is simply the combination of average incubation periods in vectors and humans and assumes that an infectious vector would transmit immediately upon becoming infectious. Thus, lag 1 was the average of weeks 1 to 3 in the past. The longest lag we investigated was 18–20 weeks.
- 2)
Selection of scaling factors. Similar to previous work [23], we hypothesized a scaling factor on the gravity terms would be important because these covariates described interactions that could be non-linear. Because initial attempts to fit this parameter were unsuccessful due to the effects of its non-linearity on convergence, we identified the best scaling factor (α1-α3 in Eqn. 3) for each possible between-neighborhood covariate (Eqn. 3, last 3 covariates) by fitting models using a range of fixed scaling factors (α = 0.001, 0.01, 0.1, 0.5, 1, 2). These values were chosen because they represent a range of biologically realistic functions for the relationship between gravity components (concave-up, concave-down or linear). The lowest value (i.e., 0.0001) was chosen based on convergence to the lowest DIC (representing asymptotic behavior of the best value) and for values above the highest (i.e., 2) the DIC continued to increase in the exponential part of the curve (i.e., values > 2 did not produce good fits).
- 3)
Mosquito and human case terms. We compared models with only mosquito density data (M _{ i,t-x1 } and ∑_{j} M _{ j,t-x2 } ^{ α1 }) to those with only human-case notifications (Y _{ i,t-y1 } and ∑_{j}(Y _{ j,t-y2}/f(x _{j}))^{ α2 }), and those with both types of covariate data (i.e., Eqn. 3), to investigate the role of mosquito density data.
- 4)
Spatial scale of between neighborhood interactions. For the between-neighborhood effects, we compared two scales: 1) nearest-neighbor effects (i.e., local) – where only covariate data from immediately adjacent neighborhoods were used to predict cases and 2) global effects – where data from all other neighborhoods city-wide were used to predict cases.
- 5)
Proxies describing between-neighborhood weights. For the global between-neighborhood covariates, we compared different functions for weighting between-neighborhood effects (f(x _{ j }) in Eqn. 1), including economic value (1/E _{ j }), population density (1/D _{ j }) and travel distance between neighborhoods (1/d _{ ij }; Table 1). We hypothesized that high-economy or high-density neighborhoods would attract more people on a regular basis, creating hubs for disease transmission and spatial spread. Similarly, we hypothesized that disease transmission from other neighborhoods would be more likely between neighborhoods with faster road travel. These ideas are similar to a recent study showing that dengue hotspots occur along major roads and transportation hubs [28]. Because mosquitos rarely travel beyond 200 m [29], which is mainly within a neighborhood, the weightings were only applied to the terms with case notification data, β_{5} and β_{6} (Eqn. 3), and not the global mosquito term, β_{4} (Eqn. 3).
Because Steps 1 and 2 were not the focus of our analysis, results from these analyses are presented in the Supplementary Material (Additional file 2: Figure S1, Additional file 3: Figure S2, Additional file 4: Figure S3, Additional file 5: Figure S4, Additional file 6: Figure S5, Additional file 7: Figure S6). Results from Steps 3–5 are reported in the main text.
Model evaluation
All steps were conducted using data from Vitoria from week 45 of 2007 through 2011, thus withholding data from 2012 for evaluation of the final model by out-of-sample prediction (i.e., forecasting). As a second means of model validation, we used the best model selected from Vitoria data on data from two other cities: GV and SL. For this, we re-estimated parameters using the best fit Vitoria-derived model structure from our model selection procedure and covariate data from each other city. Again, we only used a portion of the data for parameter estimation and predicted both this in-sample data as well as the remaining (out-of-sample) data. Because the magnitude and direction of parameter values in the three cities were so different, we did not attempt to predict data in the other two cities using parameters estimated from Vitoria covariate data. Instead, we compared the city-specific parameters.
We also conducted Steps 1, 2, 3 and 5 (above) on data from GV and SL in order to evaluate the generality of conclusions drawn based on the Vitoria time series and to gain a better understanding of how the best model may differ due to city-specific circumstances. The latter two cities did not have as much data: GV in-sample – 90 weeks, GV out-of-sample – 30 weeks, SL in-sample – 86 weeks, and SL out-of-sample – 13 weeks. In-sample data were from 2009 and 2010 while out-of-sample data were from 2011.
Results
Role of mosquito data
Goodness-of-fit for the best model selected using data from Vitoria
City | In-sample | Out-of-sample | ||||
---|---|---|---|---|---|---|
Neigh. | City | Weeks | Neigh. | City | Weeks | |
Vitoria | 0.63 | 0.91 | 217 | 0.49 | 0.30 | 34 |
GV | 0.58 | 0.91 | 90 | 0.49 | 0.63 | 30 |
SL | 0.51 | 0.86 | 86 | 0.33 | 0.70 | 16 |
Utility of proxies for weighting between-neighborhood case data
Although the DIC score was lowest for this full, “best” model, the mlCPO’s (data not shown) and r values were very similar for all proxies of neighborhood connectivity (Figure 2C). Thus, although the mlCPO’s and r values followed the same rank order as the DIC values, the high similarity of r values from models with different proxies for neighborhood connectivity did not indicate biologically important differences between the models in any of the cities (Figure 4). In summary, we found that models including global between-neighborhood effects in addition to within-neighborhood effects performed best and that all 3 types of covariates (mosquito density, case notifications and the interaction of these two covariates), but specific proxies for weighting global connectivity were similar to one another.
Generality of the vitoria model
The best model for each city included quite different lag times between mosquito density and cases: 13–15 weeks for Vitoria, 1–3 or 2–4 weeks for GV and 6–8 weeks for SL (Additional file 3: Figure S2, Additional file 5: Figure S4 and Additional file 7: Figure S6). However, a lag of 1–3 weeks was always best for the case notification data.
Discussion
Role of mosquito data
We found that even with city-wide household-level mosquito surveillance data, the relationship between mosquito density and cases is weak. Although MI-dengue has been effective at decreasing cases city-wide by basing spatial prioritization on within-neighborhood data on mosquito density and recent infections in humans [16], our results highlight that additional data may be useful for further improvements in preventing cases of dengue city-wide. Previous work has similarly found a weak [5-8] or even negative [3,4] relationship between household mosquito density and cases. Part of the reason for the obscured understanding of the role of mosquito densities can likely be attributed to high variation in vector competence across relatively fine spatial and temporal scales [30], emphasizing that surveillance for infected mosquitos should be prioritized. In fact, a strong relationship between the density of infected mosquitos and cases has been observed [8]. In our system, a new technology which monitors the density of infected mosquitoes by serotype, MI-Virus, was recently developed but has not been implemented long enough or city-wide in order for us to have evaluated the data in this study (although the information provided by MI-virus is already being used for spatial targeting of vector control where it is available). As city-wide MI-Virus data become available at the same spatial scale as the density estimates, analyses should be extended to include these data, which may lead to more accurate guidance for spatial prioritization of vector controls. Similarly, to extract the most information from the MI-virus data, it will be important to obtain data on human diagnoses at the level of serotype because the relationship between mosquito density and human cases depends on the interaction of serotype-specific pre-existing immunity and the prevalence of different serotypes [31,32].
Understanding the role of mosquito density in predicting human cases of dengue fever under any experimental design is complicated by sampling scale and variability. It is thought that the mosquito density required to sustain transmission is in fact very low [33]. If the sampling techniques used to enumerate mosquito density are too coarse to distinguish prevalence values around the transmission threshold, then it is possible that a sampling protocol with more replication, or a trapping technology that captures more mosquitos, is required. Studies that aim to determine the precision needed to distinguish low mosquito densities (i.e., near the transmission threshold) with adequate precision, such as mass trapping in enclosed mosquito populations of known sizes using various levels of replication and spatial arrangements, are needed to assess accuracy and precision of mosquito surveillance data. Likewise, better quantification of mosquito thresholds that permit transmission among humans is important for choosing appropriate trapping parameters.
In Brazil, routine vector control occurs city-wide throughout the year following national vector control guidelines [34]. Very broadly, personnel move through entire cities, block by block, neighborhood by neighborhood, in a systematic manner over the course of several months, mainly applying larvicide and source reduction. Documentation of these efforts was too sparse to be included in our models, but we do not expect that they would have obscured our ability to quantify the relationship between mosquito density and cases because they target immature stages and our system quantifies gravid adult females at a weekly scale. Additional controls are spatially targeted based on mosquito surveillance data, dengue cases data and data on infected mosquitos when they are available. However, controls based on human cases are often too late to prevent transmission because suspected cases are not confirmed until 6–8 weeks after notification. In cities where MI-dengue surveillance is conducted, the additional control activities can be targeted to blocks with the highest mosquito densities (or blocks with infected mosquitos where MI-Virus data are available) very rapidly after the mosquito populations achieve high numbers because the longest time lag between trap checking is one week and data can be visualized on the on-line MI-dengue mapping system immediately after a trap is examined [15]. If MI-dengue-based vector-control work varies in intensity non-randomly, as is likely the case due to variability in efficacy that depends on urban structures, and some blocks are responsible for more transmission than others, the combined effect could be a weak relationship between mosquito density and human cases. Moreover, the MI-dengue-based vector-control activities could explain the different best fitting lags for mosquito data among cities, if for example, in some cities the lag between transmission and available data/response is consistently longer than in other cities. Better documentation of the timing, intensity and effectiveness of vector control work in response to MI-dengue surveillance data is needed to investigate how these activities may affect interpretation of how to use mosquito density data for strategic planning of vector control work.
Although models including only mosquito data performed more poorly than those containing only case-notification data, the interaction between mosquito density and case notifications was strongly significant in all three cities. Thus, consideration of the mosquito-human interaction is important in order to more accurately predict cases in space and time. Theoretical work has similarly found a low correlation between R _{ 0 } (the average number of secondary cases in a naïve population) and mosquito density within an area due to human movement [10]. Also, when the mosquito population is highly heterogeneous, frequent travel to areas with high mosquito density can cause an epidemic or sustain low levels of transmission (depending on connectivity levels) [35], providing mechanistic insight into why mosquito density alone may not be a good predictor of human cases. The importance of between-neighborhood effects in our models suggests that movement among neighborhoods is an important driver of dengue dynamics and that the neighborhood scale, given appropriate movement data, may be effective at capturing mosquito-human interactions.
Utility of proxies for weighting between-neighborhood case data
We found that global between-neighborhood effects were stronger than either nearest-neighbor effects or within-neighborhood effects, suggesting that many infections occurred distant from the home neighborhood. Our finding that non-local effects within a city impact spatial dynamics is similar to previous work where significant spatio-temporal clustering occurs at distances up to 2.8 km [36] and where 34.7% of cases did not show any spatio-temporal clustering [11]. However, the stronger role of non-local relative to local spatial coupling in our study contrasts the finding that house-to-house human movement may predominantly drive spatial spread [9]. This discrepancy in the relative role case data from further distances may at least partly be due to differences in urban characteristics and human behavior.
While our results show that city-wide cases impact how mosquito density translates to human cases, we were not able to understand its mechanistic nature more fully given the available data. We hypothesized that economic values, population densities or travel time on roads may be good approximations to commuting patterns, but weighting between-neighborhood effects by these metrics did not explain significantly more variation than in unweighted mixing between neighborhoods. This may be because when the force of infection is high in several neighborhoods simultaneously, the probability of contact (and hence transmission) is increased to most neighborhoods, thus diluting the role of more specific patterns of connectivity (similar idea to theoretical work showing that high rates of movement increase overall transmission [35]). However, because of the importance of the term for the interaction between mosquito density and cases, it is possible that a more direct measure of neighborhood connectivity (e.g. measurements of between-neighborhood human movement) would improve the predictive power of our model by making more accurate spatial predictions when transmission rates are lower.
Generality of the vitoria model
Our analysis showed that the best model (as determined using Vitoria data) performed quite well at multiple tests of predictive power: forecasting future data in Vitoria as well as prediction of in- and out-of-sample data in two additional cities. This emphasizes that the general structure of the Vitoria model is a useful framework for quantifying different scales of spatial coupling in different cities. However, because the operational scale of vector control is the city block, using our model structure with block-level case and mosquito surveillance data will be most useful for directing operational work spatially.
We expected that the lag-time between mosquito density and human cases would approximate the virus life cycle (i.e., extrinsic incubation period + search time + intrinsic incubation period). While this was true for GV and SL (2 and 6 week lags), for Vitoria, the strongest signal was at a 13-week lag (although a strong signal was also observed at 4 weeks). The difference between cities in the most significant lag time between mosquito density and cases could be due to differences in the temporal patterns of vector control work (i.e., variable resources over time), the relative emphasis of different types of control (i.e., response-based versus prevention-based), or the total amount of resources available to conduct vector control (i.e., ability to respond to some versus many high-risk sites). We attempted to investigate these factors using vector control data from the 3 cities but we discovered that much of the control work was unrecorded. A study that includes a standardized method for recording the dates, times, location, type and amount of vector control - alongside MI-dengue surveillance - will be instrumental in interpreting the effects of control on the lag between mosquito density and human cases, and ultimately on reducing uncertainty on how to spatially prioritize vector control work.
The strong predictive ability of the case data alone shows that reasonable quantitative neighborhood-level predictions, especially with regards to the timing and magnitude of outbreaks, can be made from case notification data in the absence of mosquito surveillance data. Additionally, although results from Vitoria are based on almost 5 years of weekly data, similarly good fits and forecasts were possible in GV where less than 2 years of weekly data were available. However, the best forecasts were from models that included only space-time autocorrelation, instead of biological covariates (models not presented here). Thus, if the interest is in prediction for response planning, a non-mechanistic saturated model based on autocorrelation is likely to be the best approach. We did not present these models because our interest was in gaining an understanding of the relative role of biological factors and their spatial scales. Furthermore, the case data are often not available until about 6–8 weeks after diagnosis, which is why it is important to explore the utility of other data sources that may be available sooner.
Conclusion
A mechanistic understanding of how mosquito density maps to disease transmission among humans is crucial for the development of quantitative tools that could guide spatial prioritization of vector control [10,35,37]. Despite the demonstrated efficacy of MI-dengue at preventing cases of dengue fever in several cities [16], our current work emphasizes that even further case reductions may be achieved if spatial prioritization occurred by additionally considering city-wide neighborhood connectivity – i.e., prioritizing highly connected areas with high mosquito density. As we did not find that readily available proxies of neighborhood connectivity explained spatial coupling, direct measures of city-wide connectivity (e.g., space use by humans [13]) seem important for maximizing the preventative utility of mosquito surveillance data. Once these data are available, they can be used to identify which areas with high densities of mosquitos are most critical for targeting vector control in order to minimize transmission of dengue viruses among humans. A complimentary approach is to develop a spatially-explicit disease dynamic model that could be used to estimate city-wide connectivity, identify transmission hotspots and identify strategies of vector control that minimize city-wide cases. These are the goals of our ongoing research. Future research should also include city-wide MI-virus data as they become available. Ideally, case data should be collected at the block level, the operational unit, and serotype-specific case data are important for a better understanding of how to employ mosquito density data for spatial prioritization of vector control.
Declarations
Acknowledgements
Thanks to the dengue control program managers in the cities of Vitoria, GV and SL, especially to André Capezzuto, José Batista and Maria Jose Lanza, respectively, for help in obtaining the data for analyses. KMP, ADL, DTSH and CTW were funded by the RAPIDD program of the Science and Technology Directorate, U.S. Department of Homeland Security, and the Fogarty International Center, NIH. KMP was also funded by USDA-APHIS-WS during the latter part of the study. AEE, CMT, KHL and KSP were funded by CNPq (Pronex-Dengue – grants # 550131/2010-8 and Doenças Neglenciadas # 404211/2012-7) and ICNT-Dengue. DTSH acknowledges the David H Smith postdoctoral fellowship for funding. SC acknowledges funding from the University of Tasmania Research Enhancement Grant Scheme (grant # C 20897). CBL was supported by the NSF Graduate Research Fellowship Program under grant DGE-1321845. The views and conclusions contained in this document are those of the authors and should not be interpreted as necessarily representing the official policies, either expressed or implied, of USDA-APHIS-VS, USDA-APHIS-WS, DHS or NIH.
Authors’ Affiliations
References
- Morrison AC, Zielinski-Gutierrez E, Scott TW, Rosenberg R. Defining challenges and proposing solutions for control of the virus vector Aedes aegypti. PLoS Med. 2008;5(3):e68.View ArticlePubMed CentralPubMedGoogle Scholar
- Huy R, Buchy P, Conan A, Ngan C, Ong S, Ali R, et al. National dengue surveillance in Cambodia 1980–2008: epidemiological and virological trends and the impact of vector control. B World Health Organ. 2010;88(9):650–7.View ArticleGoogle Scholar
- Honorio NA, Nogueira RMR, Codeco CT, Carvalho MS, Cruz OG, Magalhaes M, et al. Spatial evaluation and modeling of dengue seroprevalence and vector density in Rio de Janeiro, Brazil. Plos Neglect Trop Dis. 2009;3(11):e545.View ArticleGoogle Scholar
- Lin CH, Wen TH. Using Geographically Weighted Regression (GWR) to explore spatial varying relationships of immature mosquitoes and human densities with the incidence of dengue. Int J Environ Res Public Health. 2011;8(7):2798–815.View ArticlePubMed CentralPubMedGoogle Scholar
- Scott TW, Morrison AC. Vector dynamics and transmission of dengue virus: implications for dengue surveillance and prevention strategies vector dynamics and dengue prevention. In: Rothman AL, editor. Dengue virus. 3382010th ed. 2010. p. 115–28.View ArticleGoogle Scholar
- Thammapalo S, Nagao Y, Sakamoto W, Saengtharatip S, Tsujitani M, Nakamura Y, et al. Relationship between transmission intensity and incidence of dengue hemorrhagic fever in Thailand. Plos Neglect Trop Dis. 2008;2(7):e263.View ArticleGoogle Scholar
- Mammen MP, Pimgate C, Koenraadt CJ, Rothman AL, Aldstadt J, Nisalak A, et al. Spatial and temporal clustering of dengue virus transmission in Thai villages. PLoS Med. 2008;5(11):e205.View ArticlePubMed CentralPubMedGoogle Scholar
- Yoon IK, Getis A, Aldstadt J, Rothman AL, Tannitisupawong D, Koenraadt CJ, et al. Fine scale spatiotemporal clustering of dengue virus transmission in children and Aedes aegypti in rural Thai villages. PLoS Negl Trop Dis. 2012;6(7):e1730.View ArticlePubMed CentralPubMedGoogle Scholar
- Stoddard ST, Forshey BM, Morrison AC, Paz-Soldan VA, Vazquez-Prokopec GM, Astete H, et al. House-to-house human movement drives dengue virus transmission. Proc Natl Acad Sci U S A. 2013;110(3):994–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Stoddard ST, Morrison AC, Vazquez-Prokopec GM, Paz Soldan V, Kochel TJ, Kitron U, et al. The role of human movement in the transmission of vector-borne pathogens. PLoS Negl Trop Dis. 2009;3(7):e481.View ArticlePubMed CentralPubMedGoogle Scholar
- Vazquez-Prokopec GM, Kitron U, Montgomery B, Horne P, Ritchie SA. Quantifying the spatial dimension of dengue virus epidemic spread within a tropical urban environment. PLoS Negl Trop Dis. 2010;4(12):e920.View ArticlePubMed CentralPubMedGoogle Scholar
- Danon L, House T, Keeling MJ. The role of routine versus random movements on the spread of disease in Great Britain. Epidemics. 2009;1(4):250–8.View ArticlePubMedGoogle Scholar
- Vazquez-Prokopec GM, Bisanzio D, Stoddard ST, Paz-Soldan V, Morrison AC, Elder JP, et al. Using GPS technology to quantify human mobility, dynamic contacts and infectious disease dynamics in a resource-poor urban environment. PLoS One. 2013;8(4):e58802.View ArticlePubMed CentralPubMedGoogle Scholar
- De Melo DPO, Scherrer LR, Eiras AE. Dengue fever occurrence and vector detection by larval survey, ovitrap and MosquiTRAP: a space-time clusters analysis. PLoS One. 2012;7(7):e42125.View ArticlePubMed CentralPubMedGoogle Scholar
- Eiras AE, Resende MC. Preliminary evaluation of the “Dengue-MI” technology for Aedes aegypti monitoring and control. Cad Saude Publica. 2009;25:S45–58.View ArticlePubMedGoogle Scholar
- Pepin KM, Marques-Toledo C, Scherer L, Morais MM, Ellis B, Eiras AE. Cost-effectiveness of novel system of mosquito surveillance and control, Brazil. Emerg Infect Dis. 2013;19(4):542–50.View ArticlePubMed CentralPubMedGoogle Scholar
- de Resende MC, de Azara TMF, Costa IO, Heringer LC, de Andrade MR, Acebal JL, et al. Field optimisation of MosquiTRAP sampling for monitoring Aedes aegypti Linnaeus (Diptera: Culicidae). Mem I Oswaldo Cruz. 2012;107(3):294–302.View ArticleGoogle Scholar
- Honorio NA, Codeco CT, Alvis FC, Magalhaes M, Lourenco-De-Oliveira R. Temporal distribution of Aedes aegypti in different districts of Rio De Janeiro, Brazil, measured by two types of traps. J Med Entomol. 2009;46(5):1001–14.View ArticlePubMedGoogle Scholar
- Lourenco-de-Oliveira R, Lima JBP, Peres R, Alves FD, Eiras AE, Codeco CT. Comparison of different uses of adult traps and ovitraps for assessing dengue vector infestation in endemic areas. J Am Mosquito Contr. 2008;24(3):387–92.View ArticleGoogle Scholar
- Maciel-de-Freitas R, Peres RC, Alves F, Brandolini MB. Mosquito traps designed to capture Aedes aegypti (Diptera: Culicidae) females: preliminary comparison of Adultrap, MosquiTRAP and backpack aspirator efficiency in a dengue-endemic area of Brazil. Mem I Oswaldo Cruz. 2008;103(6):602–5.View ArticleGoogle Scholar
- Favaro EA, Mondini A, Dibo MR, Barbosa AAC, Eiras AE, Neto FC. Assessment of entomological indicators of Aedes aegypti (L.) from adult and egg collections in Sao Paulo, Brazil. J Vector Ecol. 2008;33(1):8–16.View ArticlePubMedGoogle Scholar
- Rue H, Martino S, Chopkin N. Approximate bayesian inference for latent gaussian models using integrated nested laplace approximations. J Roy Stat Soc B. 2009;71(2):319–92.View ArticleGoogle Scholar
- Xia Y, Bjornstad ON, Grenfell BT. Measles metapopulation dynamics: a gravity model for epidemiological coupling and dynamics. Am Nat. 2004;164(2):267–81.View ArticlePubMedGoogle Scholar
- Ferrari MJ, Grais RF, Bharti N, Conlan AJ, Bjornstad ON, Wolfson LJ, et al. The dynamics of measles in sub-Saharan Africa. Nature. 2008;451(7179):679–84.View ArticlePubMedGoogle Scholar
- Spiegelhalter DJ, Best NG, Carlin BR, van der Linde A. Bayesian measures of model complexity and fit. J Roy Stat Soc B. 2002;64:583–616.View ArticleGoogle Scholar
- Held L, Schrodle B, Rue H. Posterior and cross-validatory predictive checks: a comparison of MCMC and INLA. In: Kneib T, Tutz G, editors. Statistical modelling and regression structures. Berlin, Germany: Springer,verlag; 2010.Google Scholar
- Chan M, Johansson MA. The incubation periods of dengue viruses. PLoS One. 2012;7(11):e50972.View ArticlePubMed CentralPubMedGoogle Scholar
- Sharma KD, Mahabir RS, Curtin KM, Sutherland JM, Agard JB, Chadee DD. Exploratory space-time analysis of dengue incidence in Trinidad: a retrospective study using travel hubs as dispersal points, 1998–2004. Parasite Vector. 2014;7:341.View ArticleGoogle Scholar
- David MR, Lourenco-de-Oliveira R, Maciel de Freitas R. Container productivity, daily survival rates and dispersal of Aedes aegypti mosquitoes in a high income dengue epidemic neighborhood of Rio de Janeiro: presumed influence of differential urban structure on mosquito biology. Mem I Oswaldo Cruz. 2009;104(6):927–32.View ArticleGoogle Scholar
- Gonçalves CM, Melo FF, Bezerra JMT, Chaves BA, Silva BM, Silva LD, et al. Distinct variation in vector competence among nine field populations of Aedes aegypti from a Brazilian dengue-endemic risk city. Parasite Vector. 2014;7:320.View ArticleGoogle Scholar
- Luo L, Liang HY, Hu YS, Liu WJ, Wang YL, Jing QL, et al. Epidemiological, virological, and entomological characteristics of dengue from 1978 to 2009 in Guangzhou, China. J Vector Ecol. 2012;37(1):230–40.View ArticlePubMedGoogle Scholar
- Liebman KA, Stoddard ST, Morrison AC, Rocha C, Minnick S, Sihuincha M, et al. Spatial dimensions of dengue virus transmission across interepidemic and epidemic periods in Iquitos, Peru (1999–2003). Plos Neglect Trop Dis. 2012;6(2):e1472.View ArticleGoogle Scholar
- Kuno G. Review of the factors modulating dengue transmission. Epidemiol Rev. 1995;17(2):321–35.PubMedGoogle Scholar
- Ministerio da saude SdVeS. Diretrizes Nacionais para a Prevenção e Controle de Epidemias de Dengue. Edited by épidemiológica DdV. http://www.dengue.pr.gov.br.
- Adams B, Kapan DD. Man bites mosquito: understanding the contribution of human movement to vector-borne disease dynamics. PLoS One. 2009;4(8):e6763.View ArticlePubMed CentralPubMedGoogle Scholar
- Rotela C, Fouque F, Lamfri M, Sabatier P, Introini V, Zaidenberg M, et al. Space-time analysis of the dengue spreading dynamics in the 2004 Tartagal outbreak, Northern Argentina. Acta Trop. 2007;103(1):1–13.View ArticlePubMedGoogle Scholar
- Cosner C, Beier JC, Cantrell RS, Impoinvil D, Kapitanski L, Potts MD, et al. The effects of human movement on the persistence of vector-borne diseases. J Theor Biol. 2009;258(4):550–60.View ArticlePubMed CentralPubMedGoogle Scholar
Copyright
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.