Dengue fever is a major public health problem worldwide, caused by any of four virus (DENV-1, DENV-2, DENV-3 and DENV-4; Flaviviridae: Flavivirus), transmitted by Aedes aegypti mosquito. Reducing the levels of infestation by A. aegypti is one of the few current strategies to control dengue fever. Entomological indicators are used by dengue national control program to measure the infestation of A. aegypti, but little is known about predictive power of these indicators to measure dengue risk. In this spatial case-control study, we analyzed the spatial distribution of the risk of dengue and the influence of entomological indicators of A. aegypti in its egg, larva-pupa and adult stages occurring in a mid-size city in the state of São Paulo. The dengue cases were those confirmed by the city' s epidemiological surveillance system and the controls were obtained through random selection of points within the perimeter of the inhabited area. The values of the entomological indicators were extrapolated for the entire study area through the geostatistical ordinary kriging technique. For each case and control, the respective indicator values were obtained, according with its geographical coordinates and analyzed by using a generalized additive model. Dengue incidence demonstrated a seasonal behavior, as well as the entomological indicators of all mosquito' s evolutionary stages. The infestation did not present a significant variation in intensity and was not a limiting or determining factor of the occurrence of cases in the municipality. The risk maps of the disease from crude and adjusted generalized additive models did not present differences, suggesting that areas with the highest values of entomological indicators were not associated with the incidence of dengue. The inclusion of other variables in the generalized additive models may reveal the modulatory effect for the risk of the disease, which is not found in this study.
Dengue is a disease caused by any of four virus (DENV-1, DENV-2, DENV-3 and DENV-4; Flaviviridae: Flavivirus) and has become a major global health issue. An epidemic similar to dengue was registered in 1699 in Central America. In Philadelphia, USA, a major epidemic occurred in 1780 and epidemics became common in the early 20th century [1]. There have been references of dengue epidemics, in Brazil, since 1923 but with no laboratorial confirmation until 1986. Dengue fever outbreaks have occurred in several states since 1986, with the identification of the DENV-1 virus in 1986 [2], DENV-2 in 1990 [3], DENV-3 in 2001 [4] and DENV-4, isolated in Manaus from 2005 to 2007 [5] and in São Paulo in 2011 [6]. All four serotypes are currently circulating in Brazil [6]. In 2010, with over one million cases of the disease, the highest reported incidence occurred in the Northern and Central-Western regions, with 621. 7 and 1,536. 8 cases per 100,000 inhabitants, respectively. In these regions, the majority of municipalities presented rates higher than 300 per 100,000 inhabitants [7], [8]. In the state of São Paulo, the incidence in 2010 was of 503 cases per 100,000 inhabitants [8]. A. aegypti is still the only vector of epidemiological importance as the transmitter of the dengue virus (DENV) in the Americas [9], [10]. Reducing levels of A. aegypti infestation is one of just a few strategies for disease control, since there is no vaccine available yet. Other undirected strategies refer to basic sanitation, garbage collection and proper water supply, which would eliminate the need for water storage. Besides, we have education activities to improve population commitment in order to eliminate breeding sites. The dengue national control program uses sampling methods to collect data in the field and to build indicators of A. aegypti presence, in the various stages of the vector life cycle, mainly traditional Breteau index and House index [11]. Although adult forms have a direct impact on virus transmission, the most used indicators to measure vector infestation are based on larvae, pupae and eggs [12]. Some studies demonstrate the direct relationship among infestation levels with the risk of epidemics in various regions of the world [13]–[16], though epidemic transmission is also reported in the presence of very low infestation levels [17]–[20]. In addition, factors associated to human population organization have a decisive role in the circulation of the virus and in the establishment of breeding sites of A. aegypti [21], which has demonstrated a great ability to adapt to different environmental contexts [9], [22], [23]. The multiple factors involved in the transmission of the disease therefore, require different approaches for understanding the forms of transmission. Over the last few decades, geoprocessing and digital mapping techniques were incorporated to the analysis of public health issues, as well as the use of spatial analysis programs to visualize the spatial distribution and spatio-temporal patterns of epidemiological data [24], [25]. These techniques allow the development of models to predict the risk of disease and territorial infestation, mapping environmental and social conditions associated to such patterns [26]. A number of studies using these techniques have analyzed the spatial and temporal distribution of A. aegypti [27], [28] and dengue transmission [29]–[32], as well as their relationships between each other [33], [34]. One of the main instruments for the operationalization of the vector control program is by monitoring the dispersion and abundance of mosquitoes through entomological indicators. It would be useful to understand the spatial and temporal patterns in small geographical scales [35]. These indicators are currently part of service routines in Brazil, but suffer from a number of limitations and particularities identified from control program perspective, which have not been studied yet. This study has the objective of evaluating the association between the spatial distribution of incidence of dengue and the entomological indicators in a middle-size city in the state of São Paulo, Brazil.
The fieldwork was carried out with the consent of the residents, and no references were made to the names or addresses of the residents, dengue patients, and control individuals. This project was approved by the Research Ethics Committee - CEP/UNITAU, declaration No. 302/12, protocol 459/09.
From February to December of 2011,4, 688 buildings were inspected and A. aegypti larvae and/or pupae were found in 186 (4%) buildings. We surveyed 1,711 households for the presence of A. aegypti adult forms. We captured 582 specimens of A. aegypti in 17. 7% of the domiciles. The percentage of adult forms was greater than that of immature forms in the buildings. We collected 11,395 eggs; 65. 4% were in the peridomicile area and 34. 6% in the intra-domicile. Whilst a larger number of eggs were found outdoors than indoors during the entire period, the number of eggs was larger indoors than outdoors during the driest months (Table 1). A similar seasonal pattern was observed for the three indicators during the hottest and wettest period of the year (February to April) (Figure 2). During this period, 195 dengue cases were reported, all of which practically occurred during the first semester, the hottest and most humid period of the year. The dengue epidemic curve follows the larva-pupa indicator curve, where there is an increase up to April and a decrease up to the last case reported in September. We also observed a similarity between the curves for the adult forms and egg indicators, i. e. , a decreasing curve up to the end of the dengue epidemic. Figure 3 presents the monthly dengue cases (points) and the quarterly entomological indicators estimated by the ordinary kriging, showing areas of different infestation intensity, with a gradient of colors varying from green (lower values) to red (higher values). We noticed a similarity among indices for all quarters, although some differences among the three indicators in the quarters can be pointed out. The Q2 quarter adult index did not agree with the larva-pupa and egg indices; in Q3 quarter, the larva and pupae indices did not agree with the adult and egg indices; in the last quarter, Q4, all the three indices were discordant (Figure 3). We observed that during Q1, the cases in March were highly concentrated in the southern region of the map, where the three indicators presented low values. During Q2, the cases in April were highly concentrated in the eastern region, with low values for eggs and adult indicators. The intensity of the epidemic was reduced during Q3 and Q4, when dengue cases were reported all over the city, albeit with different patterns of infestation. We observed the incidence of dengue as well as the presence of the vector in different stages of its biological cycle throughout the area during the entire period of the study. However, there was little spatial coincidence amongst the incidence of dengue and the intensity of the entomological indicators. The exceptions were the presence of case clusters and high adult infestation in Q2 and Q3, and also a degree of agreement between high levels of larva-pupa infestation and incidence of dengue in Q3 (Figure 3). Table 2 shows the estimated risk of dengue and its respective confidence interval obtained from the adjusted model (GAM). We noted that in most of the analyzed period the highest values of the entomological indicators did not match with the dengue cases areas. An exception was viewed from February to April with the presence of larvae-pupae related with an elevated risk of disease. In the period from May to July the indicators related to eggs and larvae-pupae were associated with a mild risk of dengue. Egg indicator did not present statistical association with cases in the multiple model from February to April, so like the adult form indicator from April to June. All the entomological indicators presented negative association with the occurrence of dengue cases from March to May, showing that areas with elevated indicators presented minor occurrence of dengue cases in the period. Model results corroborating with those observed in Figure 3 in which higher occurrence of dengue cases did not match with higher entomological indicators. The crude relative spatial risk model for dengue (space only) and adjusted model (space + entomological indicators) for the four quarters, obtained by GAM model are presented in Figure 4. Maps show a major variation in the transmission risk at different locations throughout the quarters, high in the south during Q1, in northeast during Q2, expanding further east during Q3 and to the north and south-central regions during Q4. Both crude and adjusted models present similar spatial feature, indicating lack of effects of the entomological indicators in the risk of disease incidence. In Figure 3 we observed dengue cases in different areas during the study period, like “moving” through the region. Also the higher or lower intensity of the vector are identified in several areas of the study region, however it did not follow the disease displacement, suggesting that there is not a direct spatial relationship.
In this study, we examined the spatial and temporal distribution of three entomological indicators A. aegypti mosquitos – eggs, larvae-pupae, and adults – in Sumaré. Simultaneously, we analyzed the spatial distribution of dengue incidence and its association with the indicators. The indicators presented similar seasonal behavior throughout the period and vectors were present in the entire area of the study during the dengue transmission period. Nevertheless we observed that the clusters of dengue incidence did not correspond with the distribution of entomological indicators Incidence maps of the disease obtained from crude and adjusted models did not show significant differences in dengue incidence among the areas, suggesting that the vector indicators did not influence dengue incidence in a given space in the municipality. Although we observed a few coincidences between indicators and dengue cases in some areas, in general, the case clusters did not correspond to the highest values of vectors indices. However, caution is required when analyzing these results, especially considering the environmental interventions carried out by the health services team during epidemic periods, such as the chemical treatment of the locations with confirmed cases, the intensification of measures to control breeding areas, and health education. In addition it should be noted that the response of the population to campaign pleas and dengue control measures, especially during an epidemic, may led to the elimination of breeding grounds of immature forms. It is important to consider that the study area experienced infestation for over 15 years and deficient public services, besides the intense traffic between neighborhoods. The circulation of serotype DENV-1, DENV-2 and DENV-3 in the previous years may have interfered in the transmission of the virus, once a population portion may have acquired immunity to one or more serotypes. These factors may be relevant in explaining the circulation of the virus during the entire period of the study. The longstanding A. aegypti infestation rates in the city may have had a limited role in the circulation of the virus in Sumaré in 2011. Thus it is plausible that the spatial distribution of the disease depends on the vector distribution, but not necessarily on its highest concentration. Despite the evident seasonality of vector infestation and viral transmission (temporal association), the spatial distribution of disease risk and high rates of entomological indicators did not coincide in multiple models. The spatial analysis could detect risk (incidence) on a scale of blocks and neighborhoods. In this way we obtained local estimators of risk instead of mean values of indicators of large areas. The apparent temporal correlation observed in the indicators curve over the period of transmission, was not spatially evidenced, as seen in the results of GAM models. They showed in some periods, negative association, in others, lack of statistical significance, and others risk of dengue. Furthermore, the vector adaptability to environment inhabited by human [47] makes it difficult to control the dengue transmission, because the vectors are abundant enough for triggering and maintaining the circulation of the virus. Other variables like population size, epidemic duration and climate variables seem to determine the spread of the epidemic in longstanding infestation areas, as observed by Siqueira Junior et al. [48] and Chowell et al. [49]. The spatial association among entomological indicators has been the object of a few studies conducted at various locations. In Rio de Janeiro, Honório et al. (2009) recently reported dengue infection among residents living in areas with a low mosquito density, suggesting that the infection did not occur inside the residence [31]. In Tupã, São Paulo, Barbosa and Lourenço (2010) did not find a spatial relationship between dengue incidence and larvae infestation [28]. In Bangladesh, Ali et al. (2003) found a positive association between the dengue incidence and vector infestation in areas close to hospitals [50]. In Campinas, São Paulo, Cordeiro et al. (2011) showed a positive association between the increase in the larvae density and the incidence of dengue [34]. Chowell et al. (2008) also recorded a major variation in the weekly dengue incidence among the provinces of Peru, probably because of the level of infection spread by the mosquitoes, climate variation, circulation of different serotypes, and the population' s immunological history [49]. Studies on the association between dengue transmission and entomological, environmental, socio-economical, and other indicators have presented conflicting results but reach the consensus that the dynamics of dengue transmission are complex and difficult to understand. Also authors agree that transmission depends on the environmental context and on variables that were not taken into consideration in this study, e. g. , populational immunity, circulating serotypes, and control measures [14]–[16], [27], [47], [51]. Furthermore, underreporting of cases is also assumed because some patients were either asymptomatic or had mild symptoms, which were not reported to the health services, as mentioned in various studies [20], [31], [52], [53]. According to Focks et al. , in order to prevent the transmission of the dengue, it is necessary to maintain vector infestation at critically low levels [54]. The Pan American Health Organization describes that a building index up to 0. 1% implies low risk of dengue transmission, between 0. 1 and 5%, medium risk, and above 5%, high risk [15]. However, some authors have reported dengue transmission even when the indices of entomological indicators were relatively low [14], [18]–[20]. Currently, there is no threshold for the suspicion of the risk of dengue. This study showed that the household infestation was above 5% until May, period that occurred more than 90% of the cases. Honório et al. suggested that information on the patterns of populational movement help improve the understanding about the transmission dynamics of the disease and possible locations of its incidence [31]. Getis states that only one or few infected A. aegypti mosquitoes transmit the virus to several susceptible humans within a period of a few days [55]. According to Kan, the traffic of people and mosquitoes from neighborhoods where dengue incidences have been reported can explain the shift in the pattern of the epidemics [30]. This study was the first to evaluate the spatial relationship between entomological indicators of all stages of the A. aegypti mosquito: egg, larvae and pupae, and adult, and monthly and concurrent measurements as well as dengue incidence. One of the limitations of the study is the lack of data about populational immunity and the interventions implemented in response to the autochthonous transmission that could have influenced the results. Besides, the population' s movement patterns and the elements indicating the main transmission locations are unknown [31]. By only analyzing the reported cases during the period without considering the asymptomatic patients and those who did not seek health services, the findings were certainly underestimated. The other limitation was the use of data from a surface smoothed by ordinary kriging, which despite being a linear unbiased estimator, promotes the smoothing of results, thereby overestimating the lowest and underestimating the highest values [44]. Many factors are involved in the spatial spread of an epidemic. Certain factors, e. g. , vector control programs and populational immunity to the circulating virus, may have a modulating effect on the dengue incidence. Besides, the introduction of a virus, its establishment and propagation, and the concomitant movement of various serotypes, owing to the population traffic, may also be part of the spatial dimension of the epidemic spread. In this case, infestation was not a limiting variable for transmission [32]. In this study, we were able to simultaneously analyze the incidence of dengue and conduct a survey on the entomological indicators of A. aegypti; however, we did not find a spatial correlation between the indicators and disease incidence. The infestation did not present a major variation in intensity and was not a limiting or determining factor of dengue incidences in a given space in the municipality. None of the different stages entomological indicators in the vector' s lifecycle was a predictor of disease occurrence in areas at risk of dengue transmission. The inclusion of other variables in generalized additive models could eventually reveal other modulating factors that have an influence on spatial pattern of the disease.
Dengue is a disease caused by a virus which has four serotypes DENV-1, DENV-2, DENV-3 and DENV-4. In Americas, A. aegypti is acknowledged as the only dengue vector in America. Currently the only strategy to prevent dengue is controlling A. aegypti mosquitoes. The generalized additive model was used to understand the relationship of the indicators of the presence of eggs, larvae-pupae and adult stages of A. aegypti with the occurrence of dengue cases in a medium sized city of São Paulo state, Brazil. Dengue incidence as well as the entomological indicators in all stages of the mosquito showed a seasonal behavior. The infestation level was not a limiting or a determinant factor of the occurrence of cases in the municipality. Risk maps of the disease, from the crude and adjusted by generalized additive models, showed no differences, suggesting that the entomological indicators did not influence the incidence of dengue in the city. The inclusion of other variables in the generalized additive model could reveal the modulating effect on the disease risk, not found in this study.
medicine and health sciences
environmental epidemiology
infectious disease epidemiology
epidemiology
spatial epidemiology
Spatial Distribution of the Risk of Dengue and the Entomological Indicators in Sumaré, State of São Paulo, Brazil