Skip to main content

The spatial and temporal scales of local dengue virus transmission in natural settings: a retrospective analysis

Abstract

Background

Dengue is a vector-borne disease caused by the dengue virus (DENV). Despite the crucial role of Aedes mosquitoes in DENV transmission, pure vector indices poorly correlate with human infections. Therefore there is great need for a better understanding of the spatial and temporal scales of DENV transmission between mosquitoes and humans. Here, we have systematically monitored the circulation of DENV in individual Aedes spp. mosquitoes and human patients from Caratinga, a dengue endemic city in the state of Minas Gerais, in Southeast Brazil. From these data, we have developed a novel stochastic point process pattern algorithm to identify the spatial and temporal association between DENV infected mosquitoes and human patients.

Methods

The algorithm comprises of: (i) parameterization of the variogram for the incidence of each DENV serotype in mosquitoes; (ii) identification of the spatial and temporal ranges and variances of DENV incidence in mosquitoes in the proximity of humans infected with dengue; and (iii) analysis of the association between a set of environmental variables and DENV incidence in mosquitoes in the proximity of humans infected with dengue using a spatio-temporal additive, geostatistical linear model.

Results

DENV serotypes 1 and 3 were the most common virus serotypes detected in both mosquitoes and humans. Using the data on each virus serotype separately, our spatio-temporal analyses indicated that infected humans were located in areas with the highest DENV incidence in mosquitoes, when incidence is calculated within 2.5–3 km and 50 days (credible interval 30–70 days) before onset of symptoms in humans. These measurements are in agreement with expected distances covered by mosquitoes and humans and the time for virus incubation. Finally, DENV incidence in mosquitoes found in the vicinity of infected humans correlated well with the low wind speed, higher air temperature and northerly winds that were more likely to favor vector survival and dispersal in Caratinga.

Conclusions

We have proposed a new way of modeling bivariate point pattern on the transmission of arthropod-borne pathogens between vector and host when the location of infection in the latter is known. This strategy avoids some of the strong and unrealistic assumptions made by other point-process models. Regarding virus transmission in Caratinga, our model showed a strong and significant association between high DENV incidence in mosquitoes and the onset of symptoms in humans at specific spatial and temporal windows. Together, our results indicate that vector surveillance must be a priority for dengue control. Nevertheless, localized vector control at distances lower than 2.5 km around premises with infected vectors in densely populated areas are not likely to be effective.

Background

Dengue is a mosquito-borne viral infection that affects 500 million people every year [1]. Considering the large number of infections, dengue causes significant mortality and morbidity worldwide [2, 3]. Three quarters of human infections are estimated to be asymptomatic but a small percentage of patients can develop a severe and deadly form of dengue [1]. According to the World Health Organization, Brazil currently occupies the first place in the ranking of reported dengue cases in the world with incidence rates increasing since 2004 [4].

Dengue virus (DENV) is a positive single-stranded RNA virus belonging to the genus Flavivirus of the family Flaviviridae [5]. There are at least four distinct types of DENV that are genetically related by sharing at least 60% sequence similarity (here simply denoted as DENV1, DENV2, DENV3 and DENV4) [6], and co-circulating in several parts of the world including Brazil [7,8,9].

The DENV transmission cycle involves mosquitoes and humans, although virus circulation in other vertebrate hosts has also been reported [10, 11]. Aedes aegypti mosquitoes are mostly responsible for the urban transmission cycle [3] while Ae. albopictus have been suggested to have a role communicating sylvatic and urban cycles of DENV [12, 13]. Despite the essential role of mosquitoes in DENV transmission, there is little correlation between vector indices and human outbreaks [14,15,16,17] and hence the importance of this study. DENV transmission is highly dependent on environmental variables that influence vector population dynamics and virus-vector interactions. Macro-scale studies indicated that temperature, wind speed and precipitation are related to changes in DENV incidence (reviewed in [18]). Vector survival depends on water availability and therefore on precipitation, since larval stages of the mosquito are aquatic [19]. Wind influences flight-related activities of vector mosquitoes such as host-seeking. In fact, it has been shown that wind speed above 0.9 m/s seems to discourage nectar-feeding and biting activity [18]. Biting activity, survival rate and extrinsic incubation period (EIP) of the virus in the insect vector are all temperature-dependent. Higher temperatures reduce the duration of the EIP, increasing the likelihood of DENV transmission [20]. However, DENV transmission can only happen if there is contact between humans and mosquitoes, which depends on host-seeking activities of the vector as well as human mobility [21]. In particular, Ae. aegypti is considered a localized vector, with low dispersal range, i.e. around 60 m or less, and rarely beyond 500 m [22,23,24,25]. On this basis, Ae. aegypti females are expected to visit no more than two or three houses in their lifetime thus creating local clusters of DENV infections [26]. Finally, Stoddard et al. [27] found that house-to-house human mobility (social connections) played a key role in defining individual infection risk. Therefore, the contact rate between humans and infected mosquitoes, and the ecological, socio-economic and cultural factors must be considered when inferring the spatial and temporal scales of DENV transmission [28]. However, most of the studies on DENV transmission have focused only on one dimension, either the spatial scale of mosquito distribution or the temporal range of the transmission. Thus, the spatio-temporal dynamics of DENV transmission between mosquitoes and humans in natural and urban settings remains largely unknown. Nevertheless more effective public health policies require that we understand the scales of DENV circulation, which are essential to help predict and prevent outbreaks by vector control strategies [29]. These are at the moment the only effective intervention against outbreaks due to the absence of efficient treatments and vaccines for DENV.

The aim of this study was to estimate (i) the spatial and temporal scales of local dengue transmission in an urban area, and (ii) the effect of the environmental factors. For the first objective, we first analyzed DENV presence in individual mosquitoes and human patients. Then, in order to estimate the values and significance of the spatial and temporal scales of DENV infections in mosquitoes and humans, we developed a bivariate spatio-temporal model. This combines local ordinary kriging (a common interpolator, often used to predict the distribution, dispersal, and abundance of mosquito vectors [22]) of DENV incidence in mosquitoes and permutation of the locations of human DENV cases. For the second aim, we applied a geostatistical additive linear model [30] to identify the important environmental factors that correlated with DENV incidence in mosquitoes found in the proximity of human patients. Importantly, our model does not require the restrictive assumptions made by other spatio-temporal point processes such as marked-point process or inhomogeneous Poisson point process [31]. Hence, we have proposed and implemented a novel methodological framework that will help our understanding on the ecology and scales of local DENV transmission in natural settings.

Results

Mosquito collection and analysis

Between August 2010 and July 2011, we performed a comprehensive spatiotemporal survey of mosquitoes in the city of Caratinga utilizing 158 mosquito traps evenly distributed throughout the urban area (Fig. 1a). Aedes spp. females represented 95.9% of 763 mosquitoes captured, which confirms the preference described for the traps we utilized [32]. We observed higher numbers of Ae. aegypti (655) females compared to Ae. albopictus (73) as previously observed for urban areas [33,34,35,36]. We also captured a small number of male Aedes mosquitoes (30 Ae. aegypti and one Ae. albopictus) and four Culex mosquitoes. These results are in accordance with the types of traps we utilized that were designed for highly specific and detailed surveillance of Aedes mosquitoes and not population control [29, 37]. Spatially, Aedes mosquitoes were concentrated in the city center close to the city hall (Fig. 1b). Temporally, there was a sharp increase in mosquito numbers in the transition from spring to summer and a decrease in the fall (Fig. 1c).

Fig. 1
figure 1

Spatial and temporal distribution of Aedes mosquitoes in the city of Caratinga. a Localization of mosquito traps in the urban area of Caratinga. Main roads are in yellow, orange and red whereas local streets are in white. Areas shaded in grey are residential and green marks are major forests (source: OpenStreetMap). Filled circles represent traps with at least one mosquito capture. b Same map as in a showing the density of Ae. aegypti and Ae. albopictus mosquitoes based on the number of captures in Caratinga. Mosquito density follows the color scale shown in the figure. c Total number of female Aedes mosquitoes (Ae. aegypti and Ae. albopictus) captured during the period of the study. Mosquito numbers show an increase in the summer and remain high until the fall. d Number of DENV infected Aedes females captured during this study where virus types are indicated by color. e Number of human dengue cases that were positive for viral RNA in the blood with the DENV types indicated by color

Out of the total, we detected 43 Aedes mosquitoes (36 Ae. aegypti females, two Ae. aegypti males and five Ae. albopictus females) that were positive for DENV. Infection of male mosquitoes is indicative of DENV circulation between mosquitoes but does not directly contribute to virus transmission to humans. Considering only Aedes females, we observed a total DENV prevalence of 5.6%, or 41 individuals out of 728 (Fig. 1d). This number is higher than observed in other studies where DENV prevalence in mosquitoes ranged between 0.58 and 1.8% [38,39,40]. However, our strategy involved sample preparation and testing of individual mosquitoes, which increases the sensitivity compared to the analysis of mosquito pools in previous studies [41,42,43,44].

Regarding the diversity of DENV serotypes, we observed that the majority of wild Aedes female mosquitoes carried DENV1 or DENV3 (Fig. 1d). One A. albopictus female was positive for DENV4 but was not included in further analysis due to the small numbers. For the subsequent spatio-temporal analyses of DENV transmission, we only considered female Aedes mosquitoes infected with DENV1 and DENV3.

Identification of the parameters of the variograms for the spatial incidence of DENV1 and DENV3 in mosquitoes

We next analyzed the spatial DENV incidence in mosquitoes by fitting a variogram. A variogram is a function of the variance of an attribute between point locations, in this case DENV incidence between trap locations. This variance is usually lower at shorter distances so that closer locations are more likely to have similar values than locations far apart. The variograms of DENV1 and DENV3 incidence in mosquitoes were fitted with a restricted maximum likelihood approach using an exponential function that returned lowest errors in cross validation (via kriging) (mean error 0.0002; mean squared error 0.006). For both DENV1 and DENV3, we also calculated the variogram confidence envelopes, which indicate the amount of variation that would be expected in the case of independence between incidence and location. Figure 2a shows an almost identical autocorrelation pattern and envelopes for the two DENV serotypes (from 1000 permutations).

Fig. 2
figure 2

Similarities in the spatial dependence of DENV1 and DENV3 incidence in mosquitoes. a Variogram envelopes where dashed lines are the modelled semivariance (variance between points at a varying distance) for DENV1 (green) and DENV3 (red) incidence in mosquitoes. The continuous lines are the envelopes for DENV1 (green) and DENV3 (red). If the modelled variogram crosses the simulation envelopes it is considered significantly different which was not the case. b Variogram isotropy test where ΔM, DENV incidence, variograms at 0 (green) and 90 (red) degrees. The two curves shows that the variogram changes along distance are similar in the two orthogonal directions in space

Since the two variograms are statistically very similar, we merged the data for DENV1 and DENV3 in a single variable, ΔM. This new variable contains the values of DENV1 incidence in mosquitoes for humans infected with DENV1, and DENV3 incidence in mosquitoes for humans carrying DENV3; in this manner, we are not averaging the incidence, but unifying them in a single variable. We then tested if this variable is isotropic in space or, in other words, if the variogram changes with direction. Figure 2b shows that the variograms at 0 and 90 degrees are very similar suggesting that an isotropic model (that does not account for direction) can be used for the next analyses.

In summary, our results indicate that the incidence of DENV1 and DENV3 in mosquitoes is characterized by a similar spatial process, which is invariant with direction (North to South and East to West). The common estimated variogram parameters are a spatial range of 1.469 km (1.731 km and 1.261 km at 95% permutation interval), a noise of 0.010 (0.002 and 0.017 at 95% permutation interval) and spatially dependent variance of 0.006 (0.002 and 0.006 at 95% permutation interval). The observed spatial range of 1.469 km is larger than the average dispersal of Aedes mosquitoes [25], but lower than the known spatial scales of Aedes urban geographical distribution [45].

Human data collection and analysis

In total, we analyzed 44 human blood samples representing 47% of the dengue patients notified based on symptomatology by the city of Caratinga for the 2010/2011 outbreak. We first analyzed human samples for presence of DENV-specific antibodies (IgM and IgG) and observed that 24 out of the 44 (54.5%) were positive for anti-DENV IgG (Table 1). Detection of DENV specific IgG suggests that the population in Caratinga had previous contact with DENV. Fifteen of these 24 samples were also positive for anti-DENV IgM suggesting a primary response to an ongoing infection. We were able to detect DENV RNA by RT-PCR in 15 out of the 44 blood samples suggesting that the majority of patients did not present detectable viraemia at the time of the sample collection (Fig. 1e). Indeed, viraemia is detectable for only a short period during the early stages of infection, which usually coincides with the appearance of symptoms [46].

Table 1 Data summary of human patients

Regarding the diversity of DENV types, DENV1 and DENV3 were detected in the majority of human samples from Caratinga (Fig. 1e). The exception was DENV2 that was detected in one human patient but not in wild mosquitoes. We note that DENV1 and DENV3 were similarly observed in human and mosquito samples (as reported elsewhere [47, 48]). Thus, we then restricted our spatio-temporal analysis of human patient data to the 15 individuals (see the following sections) with confirmed presence of DENV1 and DENV3 by RT-PCR (rows indicated in Table 1).

Combining human and mosquito data: identification of the spatial and temporal scales of DENV transmission

To identify the spatial and temporal scales of virus transmission we tested if DENV incidence in mosquitoes in the proximity of patient locations was significantly higher at patient locations (as measured by the probability 1 – mcp, with mcp being the number of permutations with higher dengue incidence in mosquitoes than any other configuration of DENV incidence in mosquitoes). Therefore, we identified the most significant scenario as the one with lowest mcp value. This analysis showed that the probability of DENV incidence in mosquitoes was significantly higher within 2.5 km from patient locations and 50 days before the appearance of dengue symptoms in humans (“before human infection”) compared to a random permutation of human cases (mcp = 0.2, which means that only 0.2% randomizations returned a higher value of DENV incidence) (Fig. 3a). Notably, other combinations of 2.5 and 3 km with 30 to 70 days “before human infection” returned values of mcp lower than 5%, and therefore epidemiologically important. In terms of a possible directionality of transmission, the analysis showed that spatial pattern in human infections was more likely the result of increased incidence of DENV in mosquitoes “before human infection” (which may suggest a flow of DENV from mosquitoes to humans), since a 0.2% mcp was lower than 19% mcp for “after human infection” (which may represent DENV flow from humans to mosquitoes) and 1.4% “before and after human infection” (representing an increase in DENV incidence in mosquitoes from before to after the human infections) (Fig. 3a-c). In other words, the chance of having higher incidence of DENV in mosquitoes in proximity of patient locations was almost 100 times lower for “after human infection” than “before human infection”. For “after human infection” the optimal spatial and temporal ranges were identical to “before after infection”, but with none of the combinations resulting in an mcp < 5% (Fig. 3b). Finally, for “before and after human infection” test, the optimal ranges were 70 days and 2.5 km, with all the combinations between 2.5 km and 50 to 90 days lower than 5% mcp (Fig. 3c). In our analysis, we did not differentiate between infected A. aegypti and A. albopictus females, because the latter represented just a small proportion of the dataset, and the main interest in this work is on DENV transmission independently from the vector species. However, for completeness running the same framework only for A. aegypti mosquitoes, we obtained the same optimal spatial and temporal ranges, and spatial credible interval (2.5–3 km), but with larger credible interval for the temporal range (20–70 days) than “before human infection”.

Fig. 3
figure 3

Probability of DENV incidence in mosquitoes around patient locations. We calculated whether the probability that DENV incidence in mosquitoes around patient locations was larger than in other permutation of the locations “before human infection” (a), “after human infection” (b) or “before and after human infection” (c). The highest observed probability (1-mcp) was concentrated between 2.5 and 3 km spatial range and 30 to 70 days temporal range, with a maximum at 2.5 km and 50 days (red area) in the model “before human infection”

Modeling the association between environmental variables and DENV incidence in mosquitoes around patient locations: a spatio-temporal additive, geostatistical, linear model

The spatio-temporal modeling of the association between environmental variables and DENV incidence in mosquitoes 50 days before human infections and 2.5 km around patient locations was then analyzed. We note that this analysis is limited by the small sample size since we had to use only data for DENV1 and DENV3 infections in mosquitoes and humans. However, the combinations of mosquito survey sites and human locations at 2.5 km and 50 days, generated a dataset of slightly more than 1000 records.

Our results showed that there were three environmental variables, wind speed (P-value of 0.03), wind direction (P-value of 10−7) and air temperature (P-value of 0.002) that significantly correlated with spatiotemporal DENV incidence in mosquitoes in proximity of patient locations (Table 2). The cross validation test for spatio-temporal additive, geostatistical and linear model considering the environmental variables returned a mean error of 0.0008, a mean squared error of 0.003 and a mean squared deviation ratio of 0.999 (ideal value 1). The mean error and mean squared error are 2 and 8% of the incidence mean (0.036), respectively.

Table 2 Estimated coefficients for important variables associated with DENV incidence in mosquitoes around patient locations

Discussion

To our knowledge, this is the first report on the analysis on the spatial and temporal scales of DENV transmission between mosquitoes and humans in a small urban setting. We surveyed individual mosquitoes over a one-year period throughout most of a small urban area and analyzed a significant number of human samples collected during the concurrent dengue outbreak providing important information about the complexity of DENV circulation in mosquitoes and humans. Previous attempts focused on analyzing the clustering pattern of human dengue cases without taking into account vector populations or DENV incidence in mosquitoes (e.g. [49,50,51]). An exception to this can be found in [52], where the authors investigated the number of human dengue cases per district and the relative mosquito density (but not DENV incidence in mosquitoes), using a generalized linear mixed model. In addition, only a few studies have characterized the distribution of DENV types detected in mosquitoes and humans during the same outbreak [53, 54].

Despite the lack of research in the spatio-temporal scales of DENV transmission, it is well accepted that early warning systems and disease surveillance may benefit from this information [55, 56]. Dengue epidemiological surveillance in Brazil has often taken into account vector density and the number of people infected with DENV [57]. The former is estimated by the Breteau Index (BI), which is the number of positive containers per 100 houses inspected. For the human dengue cases, the WHO suggested the “endemic channel” rule, which identifies an outbreak if the number of cases reported is larger than two standard deviations above the “endemic channel” (weekly or monthly average incidence in the previous 5 to 7 years) in weekly or monthly reports [58]. During an outbreak, a perifocal spraying strategy is often utilized where insecticide is applied around 400 m radius of the residence where a human was detected with dengue. As pointed out by the WHO guidelines, this approach does not account for the time delay between when a person is infective to mosquitoes and the appearance of symptoms. In addition, the large number of asymptomatic cases and the capacity of infected people to transport DENV over longer distances more rapidly than infected mosquitoes pose a serious limitation to the efficacy of perifocal spraying [26]. It is noteworthy that our sampling design and statistical method do not allow us to estimate what proportion of the variance in incidence falls within the 400 m range. However, the spatial range of 1.5 km suggests that transmission occurs more at neighborhood than household level. Therefore, the radius of perifocal spraying should be extended, unless areas at higher risk of infections are known prior to an intervention.

We tested if DENV incidence in mosquitoes around locations where humans were detected with dengue was higher “before human infection”, “after human infection” or “before and after human infection”. We found that the strongest association of an increase in DENV incidence in mosquitoes occurred “before human infection”. The weaker association between DENV incidence in mosquitoes “after human infection” may suggest saturation of the infection in the mosquito population. In addition, vector control interventions after the start of the human outbreak can contribute to suppress DENV incidence in mosquitoes resulting in low values “after human infection”.

Our results indicated a strong spatiotemporal correlation between incidence of DENV in mosquitoes and the appearance of human dengue cases in Caratinga. The approach described in this paper was accurate in identifying the spatial range of transmission but proved less robust in identifying the temporal range. Human infections were strongly associated with higher incidence of DENV in mosquitoes within 2.5–3 km around patient location (optima at 2.5 km) and 30–70 days (optima at 50 days) “before human infection”. Difficulties in sampling infected mosquitoes and detecting infected individuals could have affected our ability to determine a temporal range closer to the extrinsic incubation period, which would be more intuitive. Nevertheless, temporal heterogeneity has been shown in other studies [22]. For example, Pepin et al. [52] found a temporal range of two to 13 weeks between the peak in mosquito density and human infections. Other analyses, focused only on the clusters of human cases, found a temporal lag from 2 to 10 weeks [59, 60]. Certainly, the minimum amount of time necessary from when the a mosquito acquires DENV to the point when it transmits it to a human and the human becoming infectious is generally 2 weeks under favorable climate conditions. However, host-seeking activities and the relative low rate of successful transmission can extend the time lag to well above the two to 3 weeks window [52]. Finally, the temporal accuracy of our results is probably affected by the technique employed for DENV detection in whole mosquitoes, which does not distinguish whether the vector is infectious to humans at the time of collection.

Regarding spatial ranges, we observed 1.5 km for infected mosquitoes and 2.5 km for DENV transmission between mosquitoes and humans. These ranges are larger than some vector distribution and local transmission studies, but are certainly smaller than other medium and large-scale analyses [19, 60]. The distribution of the traps cannot have affected the optimal spatial range since they are homogeneously distributed every ~ 0.2 km2 (this accounts for Aedes dispersal, estimated to be 800 m every 6 days [61]) and abundant around each patient location (32 to 82 traps) [62]. The size of the urban area, human mobility, asymptomatic cases and vector control interventions (elimination of potential breeding containers, spraying insecticide which tend to disperse the clustered distribution of mosquitoes) could explain why DENV infections in humans are not clustered. The sparsity of human dengue cases is confirmed by the negatively skewed distribution of the spatial ranges towards large values (Fig. 3). All these factors may have contributed to identify a medium-scale of 2.5 km range as most significant [21, 27, 63]. A priori knowledge about household commuting behavior could improve the identification of the spatial range for transmission and help defining zones for targeted control, especially if the scales of human movement exceed that of the vectors.

The ecological analysis (Table 2) confirmed the importance of warm air temperature (in agreement with [64]) and low wind speed in DENV transmission. Indeed, winds are generally an important component in mosquito dispersal and host-seeking strategy [65]. The importance of low wind speed in increasing DENV incidence has been found in Brazil [66, 67] and many other countries including Australia [68], Sri Lanka [69], Barbados [70], Pakistan [71], China [72], Malaysia [73] and Vietnam [74]. In our dataset, wind direction was positively associated to DENV incidence in mosquitoes. This means that winds from the North and Northeast (which are the predominant winds in Caratinga, according to the climatological normal from INMET) were more likely associated with high DENV incidence in mosquitoes. Remarkably, wind direction is the only variable present in all the five best models (in terms of mean error in cross validation), followed by temperature (four out of five best models) and wind speed (two out of five models) (Additional file 1: Table S1). Relative humidity was not selected in the first two models, but has higher frequency (three out of five best models) than wind speed. The total number of mosquitoes caught by the traps was not significant in the best five models.

We also applied a non-spatial model (simple general linear regression) to analyze DENV transmission between mosquitoes and humans. This model showed similar fitting results for wind speed and direction but did not perform better than the spatial model (Additional file 1: Table S1). However, the effect of temperature is inverted (lower temperature associated with larger incidence) due to the fact that the variance is not adjusted for the spatial autocorrelation.

Finally, we found a weak trend (P = 0.08) on the question whether infections were simply due to a larger mosquito density around patient locations “before human infection” as suggested by several studies [57, 75]. This result confirms the unclear association between pure vector indices and DENV incidence in humans and reinforces the importance of directly monitoring mosquito infection [76]. However, considering the difficulties in sampling and identifying infected mosquitoes, surveillance based solely on mosquito abundance may still be useful in assessing overall risk, especially when used with other epidemiological parameters.

It is important to point out that our work has some limitations including: (i) small number of samples for human patients and DENV infected mosquitoes; (ii) the fact we did not analyze asymptomatic humans that contribute to DENV circulation and can represent 50–90% of all infections [1]; (iii) the possibility that DENV circulates in other vertebrate hosts [10]; and (iv) the effect of human mobility [27, 63]. Due to the spatial and temporal scale of the DENV epidemic in Caratinga, it is possible that traps with recorded infected mosquitoes, are considered at the same time in different patient locations. Although this is true, we have not found particular traps recurring more often than others. In fact, 60% of the traps are recurring in more than 50% of the patient locations. In addition, the fact that we tested for randomness of spatial and temporal scales within different scenarios (“before human infection”, “after human infection” and “before and after human infection”) has certainly reduced the risk of bias. However, since our analysis only focused on the comparison between scenarios and the accuracy of the obtained parameters, it is possible that additional transmission components, such as human mobility, could significantly change our parameterization.

We observed no evidence that DENV serotypes influence the spatial process (autocorrelation) of infected mosquitoes (Fig. 2a). This result must be confirmed (especially due to the small sample size), since some serotypes may influence mosquito behavior [77]. In contrast, as found elsewhere (e.g. [78]), the distribution of mosquitoes infected with different DENV serotypes is patchy and localized to a small proportion of the territory when the spatial autocorrelation is accounted for.

Conclusions

Dengue is currently the most prevalent viral disease in the world. Mosquito control strategies can have a significant impact on transmission although vector indices show little correlation with human outbreaks [14]. Thus, better understanding of DENV circulation between mosquitoes and humans in natural settings is required. More accurate estimations of the spatial scales of DENV transmission can help design more efficient surveillance and vector control systems. The method proposed in this work was able to infer parameters of DENV transmission that are coherent with experimental analysis in the laboratory. Our spatio-temporal analysis indicates that local incidence of DENV in mosquitoes is a great indicator of the probability of human dengue cases and should be a major target for surveillance strategies meant to raise preparedness since localized vector control interventions are not likely to work. Altogether, this work contributes to the understanding of the spatial and temporal transmission of DENV in endemic urban areas. The developed framework can be used for any spatiotemporal process with reference points (in statistics, those points containing only the information of presence, like the patient locations) when no initial assumptions about reference (patient) and target (mosquito traps) points can be made.

Methods

Study area

The study was conducted in the municipality of Caratinga located in the east of Minas Gerais state, Southeastern Brazil (19°47′24"S, 42°08′20"W). Caratinga has a total territorial extension of 1,258,778 km2 with approximately 90,000 inhabitants. Atlantic forest is the biome of the municipality with 21 °C median annual temperature and 80% humidity (data from the Brazilian Ministry of Agriculture, INMET) [79]. Caratinga is a dengue endemic region with 1274, 91, 409 and 51 confirmed human cases each year from 2007 to 2010 prior to this study (data from the Brazilian Ministry of Health, SINAN) [80].

Mosquito collection

We took advantage of a mosquito monitoring system that has been utilized by Caratinga and other municipalities in the state of Minas Gerais [37]. This system uses mosquito traps known as MosquiTRAP [29, 81] that are placed 300 m apart, following a 300 × 300 m grid design, in a study area. Traps are placed outdoors, usually in the front or backyard of the property where access is restricted. Traps are never placed in open areas accessible to the general public. Whenever the 300 m distance between traps could not be achieved, traps were placed in a distance of less than 300 m, never more than this limit. In Caratinga, 158 MosquiTRAPs were installed at an average density of 5.15 traps/km2 (3.2–12.5 traps/km2, 95% confidence interval) [29]. These traps covered 14.2 km2 of the urban and peri-urban areas of the city of Caratinga (Fig. 1). Between August 2010 and July 2011, traps were inspected by local health authorities once a week and each captured mosquito was collected and visually classified by species and gender [82]. Mosquitoes were individually stored in liquid N2 before being processed for RNA extraction using Trizol reagent (Invitrogen, Carlsbad, USA).

Human blood samples

In order to compare DENV circulation in mosquitoes to human dengue, we obtained blood samples from patients that sought medical attention in the city of Caratinga between February and July of 2011, concomitantly with mosquito captures. Human samples were collected by professional nurses employed by the health system of Caratinga as part of a city surveillance plan. Blood samples were collected at home from patients that had any dengue related symptoms such as fever, anorexia, myalgia, arthralgia and headaches and requested assistance. We obtained 44 samples representing 47% of the 94 patients that sought medical attention by health professionals from the city of Caratinga. Samples were collected at different times after the appearance of symptoms ranging from 3 to 87 days. Blood samples were mixed with EDTA as an anticoagulant and stored at 4 °C. Serum was obtained from blood samples and inactivated at 56 °C for 30 min. Inactivated serum was tested for the presence of DENV specific IgM and IgG utilizing the Panbio Dengue IgM Capture ELISA (Alere, Waltham, USA) and Dengue IgG ELISA Test (HUMAN Diagnostics, Wiesbaden, Germany). For each patient, we utilized information of time and location based the date of onset of symptoms and home address reported to health authorities. Total RNA extraction from human blood samples was performed using Trizol LS Reagent (Invitrogen). Unlinked anonymous testing of human blood samples was approved by the ethics committee in research (COEP) of Universidade Federal de Minas Gerais (number 415/04 to EGK).

DENV detection by RT-qPCR and RT-PCR

RNA extracted from whole mosquitoes and human blood samples were screened for the presence of DENV RNA by RT-qPCR. For this purpose, we designed specific primers targeting a conserved region in the 5’ UTR of the DENV genome of different virus types (Additional file 2: Table S2). Samples that showed DENV amplification by RT-qPCR were subsequently tested using a previously described nested RT-PCR strategy [83]. This strategy utilizes DENV generic primers on the first PCR followed by a second reaction using primers targeting a region of the NS5 gene that are specific for each virus type (Additional file 2: Table S2). Reverse transcription was performed using total RNA (200 ng to 1 μg) and random primers (300 ng/μl) using MMLV reverse transcriptase. qPCR was performed using Power SYBR Green I (Applied Biosystems, Foster City, USA) and 2 pmol of specific primers. Conventional PCR was performed using 10 pmol of gene specific primers. Importantly, in the case of mosquitoes, our strategy is not able to infer whether individuals are infective at the time of collection.

Statistical analyses

Statistical analyses comprised: (A) parameterization of the variogram for incidence of each DENV serotype in mosquitoes and analysis of their differences; (B) estimation of the spatial and temporal scales in local DENV transmission and calculation of their significance by 10,000 permutations of patient locations (home addresses of each human case was used as patient location); and (C) spatio-temporal additive, geostatistical, linear model on the incidence of DENV in mosquitoes found around patient locations and using environmental data.

Abundance of adult female mosquitoes is considered the most appropriate measure of entomological risk [26]. Therefore, we did not use two DENV infected and 31 uninfected Aedes males in our spatiotemporal analysis. We considered females of Ae. aegypti combined with Ae. albopictus, since we were more interested in the regional dynamics of DENV transmission and not specific contributions of each vector. However, we recognize that the two vectors have a different ecology and behavior that can affect the spatial and temporal scales of DENV transmission [84, 85]. We compared our results using both mosquito species to Ae. aegypti alone and observed no major differences. A separate analysis of Ae. albopictus was not possible due to the limited amount of traps with infected (only five found in two traps) and uninfected Ae. albopictus, which are predominantly located on the edges of the city.

We used the variogram parameterization from (A) to inform kriging (value of a function at a given point by computing a weighted average of the known values of the function in the neighborhood of the point) in (B). Then we used the spatial and temporal ranges obtained in (ii) to produce the response variable for (C).

Parameterization of the variogram

For each mosquito trap we calculated the incidence of DENV in mosquitoes (which is the ratio between number of infected mosquitoes for each DENV serotype and the total number of captured individuals in a trap in a week). DENV incidence in mosquito variogram parameters (temporal range, t, and spatial range s) are fitted with a restricted maximum likelihood approach [86] using a spatio-temporal exponential function, where d and h are the spatial and temporal Euclidian distances, respectively, between trap locations.

$$ \widehat{\gamma}\left(s,t\right)=\kern0.5em \exp \left(-\frac{d}{s}\kern0.5em +\frac{d}{s}\kern0.5em \frac{h}{t}-\frac{h}{t}\kern0.5em \right) $$

The exponential function form was chosen since it returned the lowest errors in cross validation (via kriging) compared to spherical and Gneiting functions [87]. We tested if there were any differences in the spatio-temporal process (in terms of variogram parameterization) of the incidence for each DENV serotype in mosquitoes. To evaluate this, we permuted data values of DENV incidence at the spatial locations to produce upper and lower envelopes of the variogram for each virus type [88, 89]. We then assumed that the variograms of DENV serotypes are identical if the area within the lower and upper envelopes overlap by 95%. Finally we explored the presence of anisotropy (i.e. when the spatial variation around each point have a component of direction) by plotting the variograms at different directions [90].

Spatial and temporal scales in local DENV transmission and their statistical significance

The spatial and temporal ranges of DENV transmission were estimated by finding those scales that maximize the predicted mean incidence of DENV in mosquitoes at the patient location, i.e. the geographical coordinates at the exact home address, compared to the predicted mean from a set of permutations of patient locations. In practice:

  1. 1.

    Define a grid of spatial and temporal ranges.

  2. 2.

    For each spatial and temporal range combination:

    1. a.

      For each patient location, select those mosquito traps (and relative DENV incidence) within the selected spatial and temporal ranges around the patient, and considering only the mosquitoes infected with the same DENV serotype detected in the patient.

    2. b.

      For each patient location, v, apply local ordinary kriging [using the parameterization obtained from section (A)] to the DENV incidences selected in the previous step to estimate the value of incidence, \( {\widehat{I}}_v \) (kriging applied individually). The mean of the estimated values, F1, is:

$$ F1=\frac{1}{N}{\sum}_{v=1}^N{\widehat{I}}_v $$

where N is the total number of patients.

  1. c.

    Create a variable “count” equal to 0.

  2. d.

    Repeat for 9999 times:

  1. i.

    Permute the location of human patients in time (e.g. swapping the time of human infection between locations)

  2. ii.

    Recalculate (a) and (b)

  3. iii.

    If the new mean of the kriging fitted values is larger than F1, then add 1 to count.

  1. e.

    Calculate mcp, the probability that DENV incidence in mosquitoes in permuted locations is higher than the one at the original patient locations:

$$ mcp=\mathrm{count}/10,000 $$
  1. f.

    Go to (a) with a new spatial and temporal range combination. If all the spatial and temporal range combinations are explored, then stop.

A lower mcp indicates lower probability that there is an alternative spatio-temporal configuration of the human patient locations with larger DENV incidence in mosquitoes. This algorithm was repeated for temporal ranges before the human infection (thereafter “before human infection”), after human infection (thereafter “after human infection”) or both (thereafter “before and after human infection”).

Preliminary analyses on mosquito spatial and temporal variogram restricted the exploration of spatial ranges between 0.5 and 5 km and temporal ranges between 0 and 100 days [30]. These spatial and temporal ranges are at the low and high limits expected for mosquito-human interactions in DENV transmission [52]. Finally mosquito DENV incidence within the optimized spatial and temporal range from human patient location, ΔM, were selected as response variable for the next step of the analysis.

Spatio-temporal additive, geostatistical, linear model on the incidence of DENV incidence in mosquitoes around human patient locations

An additive, geostatistical, linear model was applied to the incidence of DENV in mosquitoes at survey locations in proximity (in space and time) of infected patient locations, ΔM, as described in [30] but without autoregressive term and random effects. The additive, geostatistical, linear model contains a matrix of covariates X (fixed effects); a spatiotemporal correlation effect, Z; and an error component ε:

$$ \Delta {M}_{q,u}=\beta X+Z+\varepsilon $$

\( Z\sim N\left(0,{\sigma}_z^2\widehat{\gamma}\left(s,t\right)\right) \)

$$ \varepsilon \sim N\left(0,{\sigma}_{\varepsilon}^2I\right) $$

therefore:

$$ \Delta {M}_{q,u}\sim N\left(\beta X,{\sigma}_z^2\widehat{\gamma}\left(s,t\right)+{\sigma}_{\varepsilon}^2I\right) $$

where the subscripts q and u indicate the location and the time of the patient event, respectively; β is a vector of coefficients for X; Z is a one column vector with spatio-temporal normally distributed random effects with mean zero and a covariance matrix given by the product of the spatial variance, \( {\sigma}_z^2 \), and the correlation matrix, \( \widehat{\gamma} \). As shown above, \( \widehat{\gamma} \) is expressed as a function of the spatial correlation parameter, s, defining the spatial range of ΔM; and the temporal correlation, t, defining the temporal range of ΔM. Finally, ε is the independent and identically normally distributed error, with error variance \( {\sigma}_{\varepsilon}^2 \)I. The candidate covariates in X are total number of mosquitoes, temperature, atmospheric pressure, wind speed, wind direction, and relative humidity. Important variables were selected based on the improvement of the mean standard error in cross-validation from the total 63 possible combinations between variables.

Model validation was performed by leave-one-out cross validation, for which we produced three statistics: the mean error, the mean squared error and the mean squared deviation ratio (full description in [30]).

Abbreviations

DENV:

Dengue virus

EDTA:

Ethylenediaminetetraacetic acid

ELISA:

Enzyme-linked immunosorbent assay

IgG:

Immunoglobulin G

IgM:

Immunoglobulin M

RNA:

Ribonucleic acid

RT-PCR:

Reverse transcription polymerase chain reaction

References

  1. Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496(7446):504–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Gubler DJ. Dengue, urbanization and globalization: the unholy trinity of the 21st century. Trop Med Health. 2011;39(Suppl. 4):3–11.

  3. Kraemer MU, Sinka ME, Duda KA, Mylne A, Shearer FM, Brady OJ, et al. The global compendium of Aedes aegypti and Ae. albopictus occurrence. Sci Data. 2015;2:150035.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Fares RC, Souza KP, Anez G, Rios M. Epidemiological scenario of dengue in Brazil. Biomed Res Int. 2015;2015:321873.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Weaver SC, Vasilakis N. Molecular evolution of dengue viruses: contributions of phylogenetics to understanding the history and epidemiology of the preeminent arboviral disease. Infec Genet Evol. 2009;9(4):523–40.

    Article  CAS  Google Scholar 

  6. Chambers TJ, Hahn CS, Galler R, Rice CM. Flavivirus genome organization, expression, and replication. Annu Rev Microbiol. 1990;44(1):649–88.

    Article  CAS  PubMed  Google Scholar 

  7. Villabona-Arenas CJ, Zanotto PM. Worldwide spread of dengue virus type 1. PLoS One. 2013;8(5):e62649.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Heringer M, Nogueira RM, de Filippis AM, Lima MR, Faria NR, Nunes PC, et al. Impact of the emergence and re-emergence of different dengue viruses’ serotypes in Rio de Janeiro, Brazil, 2010 to 2012. Trans R Soc Trop Med Hyg. 2015;109(4):268–74.

    Article  PubMed  Google Scholar 

  9. Drumond BP, Mondini A, Schmidt DJ, de Morais Bronzoni RV, Bosch I, Nogueira ML. Circulation of different lineages of dengue virus 2, genotype American/Asian in Brazil: dynamics and molecular and phylogenetic characterization. PLoS One. 2013;8(3):e59422.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. de Thoisy B, Lacoste V, Germain A, Munoz-Jordan J, Colon C, Mauffrey JF, et al. Dengue infection in neotropical forest mammals. Vector Borne Zoonotic Dis. 2009;9(2):157–70.

    Article  PubMed  Google Scholar 

  11. Diallo M, Ba Y, Sall AA, Diop OM, Ndione JA, Mondo M, et al. Amplification of the sylvatic cycle of dengue virus type 2, Senegal, 1999-2000: entomologic findings and epidemiologic considerations. Emerg Infect Dis. 2003;9(3):362–7.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Whitehorn J, Kien DT, Nguyen NM, Nguyen HL, Kyrylos PP, Carrington LB, et al. Comparative susceptibility of Aedes albopictus and Aedes aegypti to dengue virus infection after feeding on blood of viremic humans: implications for public health. J Infec Dis. 2015;212(8):1182–90.

    Article  CAS  Google Scholar 

  13. Richards SL, Anderson SL, Alto BW. Vector competence of Aedes aegypti and Aedes albopictus (Diptera: Culicidae) for dengue virus in the Florida keys. J Med Entomol. 2012;49(4):942–6.

    Article  PubMed  Google Scholar 

  14. Bowman LR, Runge-Ranzinger S, McCall PJ. Assessing the relationship between vector indices and dengue transmission: a systematic review of the evidence. PLoS Negl Trop Dis. 2014;8(5):e2848.

    Article  PubMed  PubMed Central  Google Scholar 

  15. de Melo DP, 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.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Sanchez L, Vanlerberghe V, Alfonso L, Marquetti Mdel C, Guzman MG, Bisset J, et al. Aedes aegypti larval indices and risk for dengue epidemics. Emerg Infect Dis. 2006;12(5):800–6.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Chadee DD. Dengue cases and Aedes aegypti indices in Trinidad, West Indies. Acta Trop. 2009;112(2):174–80.

    Article  CAS  PubMed  Google Scholar 

  18. Wong GK, Jim C. Urban-microclimate effect on vector mosquito abundance of tropical green roofs. Build Environ. 2017;112:63–76.

    Article  Google Scholar 

  19. Costa ACC, Codeco CT, Honorio NA, Pereira GR, Pinheiro CFN, Nobre AA. Surveillance of dengue vectors using spatio-temporal Bayesian modeling. BMC Med Inform Decis. 2015;15(93):1–12.

  20. Wong GK, Jim C. Do vegetated rooftops attract more mosquitoes? Monitoring disease vector abundance on urban green roofs. Sci Total Environ. 2016;573:222–32.

    Article  CAS  PubMed  Google Scholar 

  21. Karl S, Halder N, Kelso JK, Ritchie SA, Milne GJ. A spatial simulation model for dengue virus infection in urban areas. BMC Infect Dis. 2014;14:447.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Azil AH, Bruce D, Williams CR. Determining the spatial autocorrelation of dengue vector populations: influences of mosquito sampling method, covariables, and vector control. J Vector Ecol. 2014;39(1):153–63.

    Article  PubMed  Google Scholar 

  23. Nino L. Use of the function semivariogram and kriging estimation in the spacial analysis of Aedes aegypti (Diptera: Culicidae) distributions. Biomedica. 2008;28(4):578–86.

    Article  PubMed  Google Scholar 

  24. Mercado-Hernandez R, Fernandez-Salas I, Villarreal-Martinez H. Spatial distribution of the larval indices of Aedes aegypti in Guadalupe, Nuevo Leon, Mexico, with circular distribution analysis. J Am Mosq Control Assoc. 2003;19(1):15–8.

    PubMed  Google Scholar 

  25. Magori K, Legros M, Puente ME, Focks DA, Scott TW, Lloyd AL, et al. Skeeter buster: a stochastic, spatially explicit modeling tool for studying Aedes aegypti population replacement and population suppression strategies. PLoS Negl Trop Dis. 2009;3(9):e508.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Getis A, Morrison AC, Gray K, Scott TW. Characteristics of the spatial pattern of the dengue vector, Aedes aegypti, in Iquitos, Peru. Am J Trop Med Hyg. 2003;69(5):494–505.

    PubMed  Google Scholar 

  27. 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 USA. 2013;110(3):994–9.

  28. Viennet E, Ritchie SA, Faddy HM, Williams CR, Harley D. Epidemiology of dengue in a high-income country: a case study in Queensland, Australia. Parasit Vectors. 2014;7:379.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Sedda L, Qi Q, Tatem AJ. A geostatistical analysis of the association between armed conflicts and Plasmodium falciparum malaria in Africa, 1997–2010. Malar J. 2015;14:500.

  31. Diggle PJ. Statistical analysis of spatial and spatio-temporal point patterns. England: CRC Press; 2013.

    Google Scholar 

  32. Favaro EA, Dibo MR, Mondini A, Ferreira AC, Barbosa AA, Eiras AE, et al. Physiological state of Aedes (Stegomyia) aegypti mosquitoes captured with MosquiTRAPs in Mirassol, Sao Paulo, Brazil. J Vector Ecol. 2006;31(2):285–91.

    Article  PubMed  Google Scholar 

  33. Cecilio SG, Junior WF, Totola AH, de Brito Magalhaes CL, Ferreira JM, de Magalhaes JC. Dengue virus detection in Aedes aegypti larvae from southeastern Brazil. J Vector Ecol. 2015;40(1):71–4.

    Article  PubMed  Google Scholar 

  34. da Rocha Taranto MF, Pessanha JE, dos Santos M, Ddos Santos Pereira Andrade AC, Camargos VN, Alves SN, et al. Dengue outbreaks in Divinopolis, south-eastern Brazil and the geographic and climatic distribution of Aedes albopictus and Aedes aegypti in 2011-2012. Trop Med Int Health: TM & IH. 2015;20(1):77–88.

    Article  Google Scholar 

  35. Piovezan R, Rosa SL, Rocha ML, de Azevedo TS, Von Zuben CJ. Entomological surveillance, spatial distribution, and diversity of Culicidae (Diptera) immatures in a rural area of the Atlantic Forest Biome, State of Sao Paulo, Brazil. J Vector Ecol. 2013;38(2):317–25.

  36. Mwangangi JM, Midega J, Kahindi S, Njoroge L, Nzovu J, Githure J, et al. Mosquito species abundance and diversity in Malindi, Kenya and their potential implication in pathogen transmission. Parasitol Res. 2012;110(1):61–71.

    Article  PubMed  Google Scholar 

  37. Eiras AE, Resende MC. Preliminary evaluation of the “dengue-MI” technology for Aedes aegypti monitoring and control. Cad Saude Publica. 2009;25(1):14.

    Google Scholar 

  38. Vilela AP, Figueiredo LB, dos Santos JR, Eiras AE, Bonjardim CA, Ferreira PC, et al. Dengue virus 3 genotype I in Aedes aegypti mosquitoes and eggs, Brazil, 2005–2006. Emerg Infect Dis. 2010;16(6):989–92.

  39. Vikram K, Nagpal BN, Pande V, Srivastava A, Saxena R, Singh H, et al. Detection of dengue virus in individual Aedes aegypti mosquitoes in Delhi, India. J Vector Borne Dis. 2015;52(2):129–33.

    PubMed  Google Scholar 

  40. Urdaneta L, Herrera F, Pernalete M, Zoghbi N, Rubio-Palis Y, Barrios R, et al. Detection of dengue viruses in field-caught Aedes aegypti (Diptera: Culicidae) in Maracay, Aragua state, Venezuela by type-specific polymerase chain reaction. Infec Genet Evol. 2005;5(2):177–84.

    Article  CAS  Google Scholar 

  41. Chow VT, Chan YC, Yong R, Lee KM, Lim LK, Chung YK, et al. Monitoring of dengue viruses in field-caught Aedes aegypti and Aedes albopictus mosquitoes by a type-specific polymerase chain reaction and cycle sequencing. Am J Trop Med Hyg. 1998;58(5):578–86.

    Article  CAS  PubMed  Google Scholar 

  42. Le Goff G, Revollo J, Guerra M, Cruz M, Barja Simon Z, Roca Y, et al. Natural vertical transmission of dengue viruses by Aedes aegypti in Bolivia. Parasite. 2011;18(3):277–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Figueiredo ML, Alfonso HL, Amarilla AA, Figueiredo LT, Aquino VH, Costa CA, et al. Detection of DENV-4 genotype I from mosquitoes collected in the city of Manaus, Brazil. Virol J. 2013;10:60.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Chen CF, Shu PY, Teng HJ, Su CL, Wu JW, Wang JH, et al. Screening of dengue virus in field-caught Aedes aegypti and Aedes albopictus (Diptera: Culicidae) by one-step SYBR green-based reverse transcriptase-polymerase chain reaction assay during 2004–2007 in southern Taiwan. Vector Borne Zoonotic Dis. 2010;10(10):1017–25.

  45. Higa Y. Dengue vectors and their spatial distribution. Trop Med Health. 2011;39(4 Suppl):17–27.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Guzman MG, Halstead SB, Artsob H, Buchy P, Farrar J, Gubler DJ, et al. Dengue: a continuing global threat. Nat Rev Microbiol. 2010;8(Suppl. 12):S7–16.

  47. Gaye A, Faye O, Diagne CT, Faye O, Diallo D, Weaver SC, et al. Oral susceptibility of Aedes aegypti (Diptera: Culicidae) from Senegal for dengue serotypes 1 and 3 viruses. Trop Med Int Health: TM & IH. 2014;19(11):1355–9.

    Article  CAS  Google Scholar 

  48. Pepin KM, Lambeth K, Hanley KA. Asymmetric competitive suppression between strains of dengue virus. BMC Microbiol. 2008;8:28.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Wijayanti SPM, Porphyre T, Chase-Topping M, Rainey SM, McFarlane M, Schnettler E, et al. The importance of socio-economic versus environmental risk factors for reported dengue cases in Java, Indonesia. PLoS Negl Trop Dis. 2016;10(9):e0004964.

  50. Telle O, Vaguet A, Yadav NK, Lefebvre B, Daude E, Paul RE, et al. The spread of dengue in an endemic urban milieu-the case of Delhi, India. PLoS One. 2016;11(1):e0146539.

  51. Naish S, Dale P, Mackenzie JS, McBride J, Mengersen K, Tong SL. Spatial and temporal patterns of locally-acquired dengue transmission in northern Queensland, Australia, 1993–2012. PLoS One. 2014;9(4):e92524.

  52. Pepin KM, Leach CB, Marques-Toledo C, Laass KH, Paixao KS, Luis AD, et al. Utility of mosquito surveillance data for spatial prioritization of vector control against dengue viruses in three Brazilian cities. Parasit Vectors. 2015;8:98.

    Article  PubMed  PubMed Central  Google Scholar 

  53. de Castro MG, de Nogueira FB, Nogueira RM, Lourenco-de-Oliveira R, dos Santos FB. Genetic variation in the 3′ untranslated region of dengue virus serotype 3 strains isolated from mosquitoes and humans in Brazil. Virol J. 2013;10:3.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Mondini A, de Moraes Bronzoni RV, Nunes SH, Chiaravalloti Neto F, Massad E, Alonso WJ, et al. Spatio-temporal tracking and phylodynamics of an urban dengue 3 outbreak in Sao Paulo, Brazil. PLoS Negl Trop Dis. 2009;3(5):e448.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Racloz V, Ramsey R, Tong SL, Hu WB. Surveillance of dengue fever virus: a review of epidemiological models and early warning systems. PLoS Negl Trop Dis. 2012;6(5):e1648

  56. Wen T-H, Lin M-H, Teng H-J, Chang N-T. Incorporating the human-Aedes mosquito interactions into measuring the spatial risk of urban dengue fever. Appl Geog. 2015;62:256–66.

  57. Vargas WP, Kawa H, Sabroza PC, Soares VB, Honorio NA, de Almeida AS. Association among house infestation index, dengue incidence, and sociodemographic indicators: surveillance using geographic information system. BMC Public Health. 2015;15:746.

    Article  PubMed  PubMed Central  Google Scholar 

  58. WHO. Dengue: guidelines for diagnosis, treatment, prevention and control. Geneva: WHO; 2009.

    Google Scholar 

  59. Kan CC, Lee PF, Wen TH, Chao DY, Wu MH, Lin NH, et al. Two clustering diffusion patterns identified from the 2001–2003 dengue epidemic, Kaohsiung, Taiwan. Am J Trop Med Hyg. 2008;79(3):344–52.

  60. Estallo EL, Carbajo AE, Grech MG, Frias-Cespedes M, Lopez L, Lanfri MA, et al. Spatio-temporal dynamics of dengue 2009 outbreak in Cordoba City, Argentina. Acta Trop 2014; 136:129–36.

  61. Honorio NA, Da Costa SW, Leite PJ, Goncalves JM, Lounibos LP, Lourenco-de-Oliveira R. Dispearsal of Aedes aegypti and Aedes albopictus (Diptera: Culicidae) in an urban endemic dengue area in the state of Rio de Janeiro, Brazil. Mem Inst Oswaldo Cruz. 2003;98(2):191–8.

  62. Sun CC, Fuller AK, Royle JA. Trap configuration and spacing influences parameter estimates in spatial capture-recapture models. PLoS One. 2014;9(2):e88025.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Barmak DH, Dorso CO, Otero M. Modelling dengue epidemic spreading with human mobility. Physica A. 2016;447:129–40.

    Article  Google Scholar 

  64. Halstead SB. Dengue virus-mosquito interactions. Annu Rev Entomol. 2008;53:273–91.

    Article  CAS  PubMed  Google Scholar 

  65. Cummins B, Cortez R, Foppa IM, Walbeck J, Hyman JM. A spatial model of mosquito host-seeking behavior. PLoS Comput Biol. 2012;8(5):e1002500.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Simoes TC, Codeco CT, Nobre AA, Eiras AE. Modeling the non-stationary climate dependent temporal dynamics of Aedes aegypti. PLoS One. 2013;8(8):e64773.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Rosa-Freitas MG, Schreiber KV, Tsouris P, de Souza Weimann ET, Luitgards-Moura JF. Associations between dengue and combinations of weather factors in a city in the Brazilian Amazon. Rev Panam Salud Publica. 2006;20(4):256–67.

    Article  PubMed  Google Scholar 

  68. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Talagala T. Distributed lag nonlinear modelling appraoch to identify relationship between climatic factors and dengue incidence in Colombo disrict, Sri Lanka. Epidemiol Biostat Public Health. 2015;12(4):e11522.

    Google Scholar 

  70. Depradine C, Lovell E. Climatological variables and the incidence of dengue fever in Barbados. Int J Environ Health Res. 2004;14(6):429–41.

    Article  PubMed  Google Scholar 

  71. Khalid B, Ghaffar A. Dengue transmission based on urban environmental gradients in different cities of Pakistan. Int J Biometeorol. 2015;59(3):267–83.

    Article  PubMed  Google Scholar 

  72. Wang C, Jiang B, Fan J, Wang F, Liu Q. A study of the dengue epidemic and meteorological factors in Guangzhou, China, by using a zero-inflated Poisson regression model. Asia-Pac J Public Health. 2014;26(1):48–57.

    Article  PubMed  Google Scholar 

  73. Cheong YL, Burkart K, Leitao PJ, Lakes T. Assessing weather effects on dengue disease in Malaysia. Int J Environ Res Public Health. 2013;10(12):6319–34.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Cuong HQ, Hien NT, Duong TN, Phong TV, Cam NN, Farrar J, et al. Quantifying the emergence of dengue in Hanoi, Vietnam: 1998–2009. PLoS Negl Trop Dis. 2011;5(9):e1322.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Taber ED, Hutchinson ML, Smithwick EAH, Blanford JI. A decade of colonization: the spread of the Asian tiger mosquito in Pennsylvania and implications for disease risk. J Vector Ecol. 2017;42(1):3–12.

    Article  PubMed  Google Scholar 

  76. Morin CW, Comrie AC, Ernst K. Climate and dengue transmission: evidence and implications. Environ Health Perspect. 2013;121(11–12):1264–72.

    PubMed  PubMed Central  Google Scholar 

  77. Maciel-de-Freitas R, Sylvestre G, Gandini M, Koella JC. The influence of dengue virus serotype-2 infection on Aedes aegypti (Diptera: Culicidae) motivation and avidity to blood feed. PLoS One. 2013;8(6):e65252.

  78. Méndez-Galván J, Sánchez-Casas RM, Gaitan-Burns A, Díaz-González EE, Ibarra-Juarez LA, Garza CEM, et al. Detection of Aedes aegypti mosquitoes infected with dengue virus as a complementary method for increasing the sensitivity of surveillance: identification of serotypes 1, 2, and 4 by RT-PCR in Quintana Roo, Mexico. Southwest Entomol. 2014;39(2):307–16.

    Article  Google Scholar 

  79. Instituto Nacional de Metereologia. Ministerio da Agricultura, Pecuaria e Abastecimento do Brasil. 2015. http://www.inmet.gov.br/portal. Accessed 1 Oct 2015.

  80. DATASUS–Sinan. Ministerio da Saude. 2016. http://tabnet.datasus.gov.br/cgi/deftohtm.exe?sinannet/dengue/bases/denguebrnet.def. Accessed 15 Jan 2016.

  81. Degener CM, de Azara TMF, Roque RA, Rosner S, Rocha ESO, Kroon EG, et al. Mass trapping with MosquiTRAPs does not reduce Aedes aegypti abundance. Mem Inst Oswaldo Cruz. 2015;110(4):517–27.

    Article  PubMed  PubMed Central  Google Scholar 

  82. Consoli RAGB, Oliveira RL. Principais mosquitos de importância sanitária no Brasil. Rio de Janeiro: Editora FIOCRUZ; 1994.

    Book  Google Scholar 

  83. de Morais Bronzoni RV, Baleotti FG, Ribeiro Nogueira RM, Nunes M, Moraes Figueiredo LT. Duplex reverse transcription-PCR followed by nested PCR assay for detection and identification of Brazilian alphaviruses and flaviviruses. J Clin Microbiol. 2005;43(2):7.

    Google Scholar 

  84. Kraemer MU, Sinka ME, Duda KA, Mylne AQ, Shearer FM, Barker CM, et al. The global distribution of the arbovirus vectors Aedes aegypti and Ae. albopictus. elife. 2015;4:e08347.

    Article  PubMed  PubMed Central  Google Scholar 

  85. Honorio NA, Castro MG, Barros FS, Magalhaes Mde A, Sabroza PC. The spatial distribution of Aedes aegypti and Aedes albopictus in a transition zone, Rio de Janeiro, Brazil. Cad Saude Publica. 2009;25(6):1203–14.

    Article  PubMed  Google Scholar 

  86. Pinheiro JC, Bates DJ, D DebRoy S, Sakar D. The Nlme package: linear and nonlinear mixed effects models, R version 3, vol. 6; 2012.

    Google Scholar 

  87. Cressie N. Statistics for spatial data. Hoboken: Wiley; 2015.

    Google Scholar 

  88. Schlather M. Construction of covariance functions and unconditional simulation of random fields. In: Porcu E, Montero JM, Schlather M, editors. Advances and challenges in space-time modelling of natural events. Berlin: Springer Berlin Heidelberg; 2012. p. 25–54.

  89. Diggle P, Ribeiro PJ. Model-based Geostatistics. New York: Springer; 2007.

    Google Scholar 

  90. Zimmerman DL. Another look at anisotropy in geostatistics. Math Geol. 1993;25(4):453–70.

    Article  Google Scholar 

Download references

Acknowledgments

We thank Luciano Moreira for providing reagents, Karla Pollyanna Vieira de Oliveira for technical help, Ecovec LTDA for installing the mosquito trap system in Caratinga, and Viv Dillon for reviewing the manuscript. We are also indebted to Eduardo Pessanha (in memoriam) and Health Surveillance authorities in the city of Caratinga for mosquito captures and blood sample collections.

Funding

This work was supported with grants from CAPES, CNPq and FAPEMIG to EGK and JTM. EGK and JTM are researchers from CNPq.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author upon request.

Author information

Authors and Affiliations

Authors

Contributions

Conceived and designed the experiments: APPV, EGK and JTM. Performed the experiments: APPV, CHPG, RPO, ATSS and LCS. Conceived and designed the statistical analysis: LS. Analyzed the data: LS, APPV, ERGRA, ANAG, BPD, EGK and JTM. Contributed reagents/materials/analysis tools: AEE, EGK and JTM. Wrote the paper: LS, APPV, ERGRA, EGK and JTM. All authors read and approved the final manuscript.

Corresponding author

Correspondence to João Trindade Marques.

Ethics declarations

Ethics approval and consent to participate

Unlinked anonymous testing of human blood samples was approved by the ethics committee in research (COEP) of Universidade Federal de Minas Gerais (number 415/04 to EGK). Human blood samples were utilized in accordance with the approved guidelines from the ethics committee.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1: Table S1.

Selection results for linear and non-linear models. (PDF 76 kb)

Additional file 2: Table S2.

Oligonucleotides primers utilized for PCR amplification of specific targets. (PDF 91 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sedda, L., Vilela, A.P.P., Aguiar, E.R.G.R. et al. The spatial and temporal scales of local dengue virus transmission in natural settings: a retrospective analysis. Parasites Vectors 11, 79 (2018). https://doi.org/10.1186/s13071-018-2662-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-018-2662-6

Keywords