Modeling Dengue Vector Dynamics under Imperfect Detection: Three Years of Site-Occupancy by Aedes aegypti and Aedes albopictus in Urban Amazonia

Abstract

Aedes aegypti and Ae. albopictus are the vectors of dengue, the most important arboviral disease of humans. To date, Aedes ecology studies have assumed that the vectors are truly absent from sites where they are not detected; since no perfect detection method exists, this assumption is questionable. Imperfect detection may bias estimates of key vector surveillance/control parameters, including site-occupancy (infestation) rates and control intervention effects. We used a modeling approach that explicitly accounts for imperfect detection and a 38-month, 55-site detection/non-detection dataset to quantify the effects of municipality/state control interventions on Aedes site-occupancy dynamics, considering meteorological and dwelling-level covariates. Ae. aegypti site-occupancy estimates (mean 0.91; range 0.79–0.97) were much higher than reported by routine surveillance based on ‘rapid larval surveys’ (0.03; 0.02–0.11) and moderately higher than directly ascertained with oviposition traps (0.68; 0.50–0.91). Regular control campaigns based on breeding-site elimination had no measurable effects on the probabilities of dwelling infestation by dengue vectors. Site-occupancy fluctuated seasonally, mainly due to the negative effects of high maximum (Ae. aegypti) and minimum (Ae. albopictus) summer temperatures (June-September). Rainfall and dwelling-level covariates were poor predictors of occupancy. The marked contrast between our estimates of adult vector presence and the results from ‘rapid larval surveys’ suggests, together with the lack of effect of local control campaigns on infestation, that many Aedes breeding sites were overlooked by vector control agents in our study setting. Better sampling strategies are urgently needed, particularly for the reliable assessment of infestation rates in the context of control program management. The approach we present here, combining oviposition traps and site-occupancy models, could greatly contribute to that crucial aim.

About the Authors:

Samael D. Padilla-Torres

Affiliation: Graduate Program in Ecology, Instituto Nacional de Pesquisas da Amazônia, Manaus, Amazonas, Brazil

Gonçalo Ferraz

Affiliations Graduate Program in Ecology, Instituto Nacional de Pesquisas da Amazônia, Manaus, Amazonas, Brazil, Biological Dynamics of Forest Fragments Project, Smithsonian Tropical Research Institute/Instituto Nacional de Pesquisas da Amazônia, Manaus, Amazonas, Brazil

Sergio L. B. Luz

Affiliation: Instituto Leônidas e Maria Deane – Fiocruz Amazônia, Manaus, Amazonas, Brazil

Elvira Zamora-Perea

Affiliation: Instituto Leônidas e Maria Deane – Fiocruz Amazônia, Manaus, Amazonas, Brazil

Fernando Abad-Franch

* E-mail: [email protected]

Affiliation: Instituto Leônidas e Maria Deane – Fiocruz Amazônia, Manaus, Amazonas, Brazil

Introduction

Dengue is the most common arboviral disease of humans [1]–[3]. About 50 million people contract dengue annually, and an estimated 22,000 die from severe forms of the disease [3], [4]. Dengue virus is transmitted by mosquitoes of the genus Aedes, particularly Aedes aegypti and Ae. albopictus [5]. In the absence of effective drugs or vaccines, prevention of dengue infections and severe dengue forms heavily relies upon vector control. However, despite massive spending and some encouraging results (e.g., [6]–[9]), neither vector populations nor, consequently, dengue transmission are currently under control; on the contrary, they are both clearly expanding [2], [10]. In South America, dengue incidence increased from ∼16/100,000 population in the 1980s to ∼72/100,000 in 2000–2007 [11].

Aedes aegypti, a species native to Africa, has successfully adapted to urban environments around the world; it preferentially breeds in artificial containers (where desiccated eggs can remain viable for months), rests within houses, and feeds on human blood [12], [13]. These traits have favored its man-mediated dispersal throughout the tropics [14], [15], and, together with its capacity to transmit dengue virus, have transformed Ae. aegypti in a major public health concern [16]. Ae. albopictus is more eclectic: it exploits both urban and rural tropical-subtropical habitats, makes use of natural and artificial breeding sites, and feeds on either humans or non-human vertebrates [17], [18]; this species, however, is less efficient than Ae. aegypti at transmitting dengue virus [18].

Dengue vector control is largely based on a combination of strategies aimed at eliminating Aedes breeding sites (either physically or by means of larvicides) and reducing adult mosquito populations (through environmental insecticide application) [6]–[9]. The design, implementation, and assessment of such strategies require detailed knowledge of vector population ecology, including the estimation of dwelling infestation rates [19], [20]. In general, vector control interventions are expected to have a negative effect on infestation by Ae. aegypti and Ae. albopictus at the local scale. Measuring such an effect requires reliable methods for ascertaining infestation; yet, detection of most animal species, including disease vectors, is rarely, if ever, perfect [21], [22]. Here, we treat infestation as the probability that a dwelling is occupied by vectors (i.e. site-occupancy) and use a hierarchical modeling approach to analyze the dynamics of site-occupancy by Ae. aegypti and Ae. albopictus. Our analysis is based on three years of oviposition trap (ovitrap) data from a central-Amazon urban setting. Taking imperfect detection into account, we quantify the effects of routine control interventions and selected environmental variables on the main indicator used in vector control program management – dwelling infestation rates – and on its temporal change.

Materials and Methods

Study Setting

With a population of about 1.8 million, Manaus (3°6′S, 60°1′W) is the largest urban center of the Amazon basin (Figure 1). The city lies on the north bank of the Negro river and is surrounded by rainforest. The climate is tropical, warm and humid, with a relatively strong seasonality of rainfall and, to a lesser extent, temperature (Figure 2). After being declared eradicated from Brazil in the 1950s [15], Ae. aegypti reinfested Manaus in the late 1990s [23] and is currently widespread across all its neighborhoods [24]. Ae. albopictus was first recorded in 2002 [25], and is now also widespread [24]. Dengue transmission is endemic (i.e., occurs continuously) in the city, with recurrent epidemics and records of all known dengue virus serotypes [26]. As in other settings, dengue control in Manaus relies on dwelling visits by municipal or state agents, who physically eliminate breeding sites or treat them with larvicides; in “emergency” situations (in practice, when dengue cases begin to soar), environmental insecticide spraying aimed at reducing adult mosquito density is also used [27]. Vector control agencies also conduct regular infestation surveys on a random sample of dwellings in each neighborhood (see details in ref. [28]). The results of these ‘rapid larval surveys’ are used to set priorities and make decisions about control interventions, with control teams usually deployed to a neighborhood when dwelling infestation rises above 2%; officially, the Brazilian control program aims to keep dwelling infestation below 1% [20].

[Figure omitted. See PDF.]

Figure 1. Study area.

Manaus, state of Amazonas, Brazil (A) and Tancredo Neves neighborhood (B).

https://doi.org/10.1371/journal.pone.0058420.g001

[Figure omitted. See PDF.]

Figure 2. Observed dwelling infestation and meteorological variables during the study period.

Dwelling infestation (%; left y axis) by Aedes aegypti and Ae. albopictus; total monthly rainfall (mm; right y axis); and monthly averages of daily mean, minimum, and maximum temperatures (°C; left y axis).

https://doi.org/10.1371/journal.pone.0058420.g002

Sampling Strategy

We selected an area of ∼250,000 m² within the Manaus neighborhood of Tancredo Neves (Figure 1) for vector monitoring. This neighborhood is frequently infested by both target mosquito species, and about 1500 dengue cases were notified during the study period (refs. [24], [29], [30] and unpublished Municipal Health Department data). The typical Tancredo Neves dwelling – our unit of occupancy analysis – consists of a brick-walled house with a courtyard in a ∼10×20 m plot. In 2008, we randomly selected 50 dwellings for monthly sampling, and in 2010 we added five more dwellings, which were also sampled on a monthly basis; all residents gave informed consent. During the first 25 months, we used a combination of 2–4 ovitraps and 0–2 Adultraps [31]; afterwards, only the more sensitive ovitraps (3 per dwelling and month) were used. Traps contained hay infusion [32] and were operated for six days/month. In total, our analyses use data from nearly 5800 trap-weeks. Each month, mosquito larvae were identified to species, with the result of each individual trap recorded separately. Thus, for each dwelling and month between September 2008 and October 2011, we have a ‘detection history’ consisting of a series of binary outcomes (present = 1 and absent = 0), with one outcome per trap and per mosquito species. These ‘detection histories’ are central to our approach: by repeatedly sampling each dwelling in each month and recording each trap’s result separately, we can derive an estimate of the rate of false-negative trap results and use it to correct infestation estimates (see below). This summarizes the hierarchical nature of our models – detection histories result from the interplay of two hierarchically juxtaposed processes, where the sampling process of detecting a species in one trap is conditioned on the biological process of dwelling occupancy by that species. The use of multiple traps per dwelling enables formal separation of the two processes.

Covariates

In order to model the relation between environmental variables and infestation we obtained daily data on total rainfall, as well as on maximum, mean, and minimum temperature, from the Brazilian National Meteorological Institute (INMET). We chose these environmental metrics, or covariates, because we considered them potentially relevant for the spatial-temporal distribution of our two target species [33]–[36]. Because we had no prior information on possible time lags between meteorological changes and their effects on local mosquito populations, we decided to relate meteorological information and each month’s occupancy in three different ways: (i) looking at meteorological covariates measured, for each month, during sampling and the previous week (denoted 1-week-lag below); (ii) looking at covariates measured during sampling and the previous two weeks (2-week-lag); and (iii) looking at covariates measured during the four weeks before sampling (4-week-lag). All meteorological measurements were standardized to mean zero and standard deviation one before analysis.

Apart from rainfall and temperature, we also registered dwelling-level traits throughout the last 13 months of monitoring. Following criteria from Tun-Lin et al. [36] adapted to our setting, we separately assessed houses and courtyards; for each of these, we defined a covariate with values of 1 (poor overall maintenance, garbage accumulation, and, for courtyards, overgrown vegetation) or 0 (well-maintained houses or courtyards). Finally, we noted whether routine control interventions were or were not performed in our study area in each of the last 13 months of monitoring. These interventions were carried out by municipal/state agents and military staff, and involved elimination of breeding sites, physically or with larvicides [20], [27]; while campaigns are designed to target all dwellings, interventions are effectively limited to houses whose owners are present at the time of the visits and allow control agents to inspect their property.

Data Analyses

Our analytical approach involved two main steps. First, we used descriptive statistics, tables, and graphs to explore the data [37], and calculated naïve infestation rates (i.e., rates that assume perfect detection of vectors) of both target species for later comparison with model-based estimates (see below). Second, we implemented a set of hierarchical models of occupancy dynamics. These models explicitly account for imperfect detection, providing estimates of detection probability (denoted p) conditioned on occupancy (denoted ψ), and treat temporal changes in occupancy as a first-order Markov process [38]–[40]. Thus, the probability of a site being occupied in month t depends on the occupancy state of that site in the previous month, t−1. This also accounts for a form of temporal autocorrelation: when observations on the same sampling unit are positively correlated, values close in time are more similar than those separated by longer periods [40]. Apart from p and ψ, our occupancy dynamics models also provide information about the probability that a dwelling that is infested at time t−1 will become uninfested at time t (named local extinction probability, or ε) or that an uninfested dwelling at time t−1 becomes infested at time t (colonization probability, γ). Models can use a variety of parameterizations to represent the same occupancy-dynamic process: we chose to focus on the estimation of ψ and γ, which, combined, can yield information on ε. Parameters p, ψ, ε, and γ can be modeled as functions of sampling or environmental covariates. For simplicity, we model covariates on p and ψ.

The main assumption of the models is population closure – i.e., site-specific occupancy remains unchanged within each sampling month. In our case, closure was guaranteed by the simultaneous assessment of all traps set within each dwelling. This design could however result in nonindependence of traps set within a single dwelling: if detection in one trap increases detection probabilities in the others, this would result in some overestimation of p and, therefore, a negatively biased ψ estimate. Finally, the usual assumption of independence of dwellings with regard to infestation also applies; if violated, this would inflate the precision of occupancy estimates.

As summarized above, estimation of p requires detection histories from repeated samples or multiple traps. A ‘011’ detection history for a dwelling and month, for example, indicates that one trap was false-negative; on the other hand, all-zero histories may result from two scenarios: either the dwelling was truly uninfested or it was infested yet all traps failed to detect the vectors [38]–[40]. Since detection probability (p) is conditional on site-occupancy, it can be interpreted as the sensitivity of the vector-detection method – its ability to detect the presence of the vectors in dwellings that are actually infested.

Models were fit by likelihood maximization and ranked according to the Akaike information criterion corrected for small sample size (AICc) [41]. Model fitting and ranking were carried out with the freely-available software PRESENCE 4.0 [42]. To avoid repetition, further details on model specification, comparison, and selection are presented in the Results section and in Tables 1, 2, and S1.

[Figure omitted. See PDF.]

Table 1. Effects of control interventions, meteorological variables, and dwelling traits on infestation rates by dengue vectors: dynamic site-occupancy models fitted to a 13-month dataset.

https://doi.org/10.1371/journal.pone.0058420.t001

[Figure omitted. See PDF.]

Table 2. Meteorological covariate effects on dwelling infestation rates by Aedes aegypti and Ae. albopictus: dynamic site-occupancy models fitted to a 38-month dataset.

https://doi.org/10.1371/journal.pone.0058420.t002

We fit occupancy dynamic models separately to (i) the 13-month subset of data for which we recorded vector control interventions and house/courtyard covariates and (ii) the full 38-month dataset. This resulted in a two-stage analysis. On the first stage, we focused on modeling the effects of control interventions on ψ, looking both at interventions that took place in the same month as sampling (denoted control_no-lag) and during the preceding month (lagged effect, control_4-week-lag). These models also consider meteorological and dwelling conditions. Since two teams were involved in vector monitoring, we also modeled detection probability (p) as a function of the observer team to account for possible differences in team performance [39], [40].

On the second stage, we set aside control interventions and focused on estimating time-dependent occupancy for the whole 38-month dataset, including the final 13 months of control-intervention monitoring. This second set of models also considered the effects of meteorological covariates on occupancy, albeit with a larger amount of data. Since we used two trapping devices during the first phase of monitoring, detection probabilities were modeled as a function of trap type, and, once again, as a function of the observer team. We also assessed the amount of bias present in naïve vs. model-derived infestation rate estimates (bias = 1 − [naïve/model-derived values]).

Ethics Statement

Sampling was carried out with permission from dwelling owners, and did not involve endangered or protected species. SLBL holds a permanent license (27733-1) from the Brazilian Institute for the Environment and Natural Resources (IBAMA) for sampling insect vectors such as the Aedes species we studied.

Results

Descriptive Results: Observed Infestation

Both vector species were detected in a high proportion of dwellings throughout the study period (Figure 2), with harmonic mean values of 0.68 for Ae. aegypti (range, 0.50–0.91) and 0.61 for Ae. albopictus (range, 0.28–0.86). There was an apparent relationship between site-occupancy and weather (Figures 2 and 3). The particularly hot and dry period of June-September 2009 coincided with a sharp decrease of Ae. albopictus infestation: observed values fell from ∼0.70–0.80 to ∼0.30–0.50. A less marked decline was also apparent for Ae. aegypti. Both species, however, quickly recovered with the onset of the rainy season. Dwelling infestation indices (the World Health Organization ‘house index’) reported by routine municipal surveillance for our study neighborhood, based on 13 ‘rapid larval surveys’ [28] carried out between October 2008 and October 2011 (Figures 3 and 4), yielded a harmonic mean of just 0.033 (range, 0.015–0.089). These descriptive results rely on the assumption that vectors were absent from sites where they were not observed; however, no perfect vector-detection method is available. The modeling results summarized in the next section address this key limitation.

[Figure omitted. See PDF.]

Figure 3. Observed and model-estimated dwelling infestation by Aedes aegypti (A) and Ae. albopictus (B).

Monthly model-derived site-occupancy estimates (solid circles, with 95% confidence intervals); monthly observed infestation (empty circles); and Ae. aegypti infestation indices derived from 13 ‘rapid larval surveys’ [28] (red circles in panel A). On the x axis, grey boxes highlight the periods in which city-wide, massive Aedes control campaigns, called Operação Impacto [29], [30], took place. Arrows indicate months in which control activities were performed in our study neighborhood (red arrows, interventions included as model covariates).

https://doi.org/10.1371/journal.pone.0058420.g003

[Figure omitted. See PDF.]

Figure 4. Bias in Aedes aegypti (left) and Ae. albopictus (right) observed infestation.

Model-derived point estimates (“Model”) correspond to the top-ranking, 38-month dynamic model for each species; observed dwelling infestation recorded during our surveys (“Observed”) and indices of dwelling infestation by Ae. aegypti reported by the regular vector surveillance system, derived from ‘rapid larval surveys’ [28] (“RLS”). Monthly values (empty circles) and quartiles 50% (horizontal line within box), 25%–75% (box lower-upper limits), 10%–90% (short lines), and 0%–100% (bottom-top lines) are shown. All values are presented as proportions.

https://doi.org/10.1371/journal.pone.0058420.g004

Modeling Results I: Effects of Control Interventions

On the first stage of our analysis we modeled the effects of vector control interventions carried out by local health authorities on site-occupancy by Ae. aegypti and Ae. albopictus. These models used data from 55 dwellings monitored from October 2010 to October 2011 with up to three ovitraps per dwelling and month. Overall, the data encompass results from 1907 ovitraps, of which 849 detected Ae. aegypti and 828 detected Ae. albopictus.

Aedes aegypti detection/non-detection data are best explained by a model with just one covariate on ψ, the average of maximum daily temperatures measured with a 2-week-lag (t_{max-2-week-lag}), which had a negative effect on site-occupancy (Table 1). The second-ranking model is also substantially supported by the data (ΔAICc = 0.73); it includes the additive effects of t_{max-2-week-lag} and control interventions carried out during the same month (control_no-lag) on ψ. The effect of temperature was again negative; this model also yielded a negative point estimate of the control coefficient, but uncertainty about this estimate is large and the 95% confidence interval overlaps zero (Table 1). Among candidate models including dwelling covariates, the one with the lowest AICc estimates a weak, positive effect of poor house condition on infestation, but, again, the estimate of this effect is too uncertain to draw any strong conclusions (Table 1).

The top-ranking model for Ae. albopictus estimates a negative effect of 1-week-lagged minimum temperatures (t_{min-1-week-lag}) on infestation; in addition, the model suggests that houses in poor condition might have been at a slightly higher risk of infestation, albeit the estimated coefficient’s 95% confidence interval includes zero (Table 1). Adding control interventions carried out the month before (control_4-week-lag) resulted in a model that fits reasonably well (ΔAICc<1). For this second model, the negative coefficient of control_4-week-lag on ψ is nevertheless small and imprecise, with confidence intervals including zero (Table 1). Finally, a model with t_{min-1-week-lag}, house condition, and control_no-lag as covariates also performed reasonably well; it estimated a small positive coefficient, again not different from zero, of control interventions (Table 1). The null model (without any covariates) and those models exploring the effects of courtyard covariates all had ΔAICc≥3 (see Table S1).

Modeling Results II: Long-term Site-occupancy Dynamics

The results in the previous section show that modeling time-specific occupancy as a function of control interventions or dwelling-level covariates did not improve the ranking of the models. Therefore, we felt justified to extend modeling to the full dataset – focusing on the potential effects of meteorological variables on site-occupancy by Aedes species. The analyses make use of a 38-month dataset including individual results of 5799 trap-weeks, which detected Ae. aegypti on 2641 occasions and Ae. albopictus on 2538 occasions. In these models, colonization probability (denoted γ) was constrained to be constant across months, while monthly local (dwelling-level) extinction probabilities (ε), of primary interest in the context of vector control, were derived from ψ and γ estimates as described in MacKenzie et al. [39], [40].

The full Aedes aegypti data were best explained by a model including the 2-week-lagged average of daily maximum temperatures (t_{max-2-week-lag}), which had a negative effect on site-occupancy probabilities (Table 2). The model with t_{max-1-week-lag} as a covariate on ψ also fitted the data well, and estimated a similar effect to that of t_{max-2-week-lag} (Table 2). The remaining models that we examined, including a null model without any covariates, performed substantially worse than these two top-ranking models (see Table S1). Among models that included rainfall covariates, the best-performing one had a ΔAICc = 6.29 and estimated a positive effect of total rainfall (r_1-week-lag) on ψ (Table 2).

Figure 3A shows monthly site-occupancy estimates for Ae. aegypti derived from the lowest-AICc model. With few exceptions, point estimates were consistently >90% (harmonic mean 0.91; range, 0.79–0.97) and showed a weak seasonal pattern (Figure 3A). Model-based infestation estimates are about 30% higher than observed values based on ovitrap results (median bias, 0.29) (Figure 4). The estimated average sensitivity of ovitraps at detecting infestation by Ae. aegypti varied from p = 0.48 (95% confidence interval 0.45–0.51) to p = 0.65 (0.63–0.67), depending on which field team performed monitoring (details not shown); 1355 out of 4553 ovitrap-weeks and 450 out of 1246 Adultrap-weeks did not detect Ae. aegypti in dwellings where other traps yielded evidence of infestation. Local extinction probability estimates were overall very low (harmonic mean ε = 0.04; range, 0.01–0.18), reaching higher values in hotter months (Figure 5A); mean site-colonization probabilities were estimated as γ = 0.66 (95% confidence interval 0.54–0.76) over the study period.

[Figure omitted. See PDF.]

Figure 5. Derived estimates of local extinction probabilities (ε) for Aedes aegypti (A) and Ae. albopictus (B).

For each species, ε estimates (bold black lines) and 95% confidence intervals (thin grey lines) were derived from the best-performing (lowest AICc) 38-month model. We also plot variation (z-scores) of average maximum temperatures during sampling days and the previous two weeks (t_{max-2-week-lag}, right y axes in each panel; colored areas); this was the meteorological covariate in the best Aedes aegypti model. On the x axis, grey boxes highlight the periods in which city-wide, massive Aedes control campaigns, called Operação Impacto [29], [30], took place; note that they coincide with months of very low ε values.

https://doi.org/10.1371/journal.pone.0058420.g005

The best-ranking Ae. albopictus model included only one site-occupancy covariate, t_{min-1-week-lag}, which had a negative effect on ψ (Table 2). The remaining models performed substantially worse (ΔAICc>20), but several of the candidate specifications we tested had convergence problems. The only model with a rain covariate estimates a positive effect of 4-week-lagged rainfall on site-occupancy by Ae. albopictus (Table 2). Site-occupancy estimates derived from the best-ranking model are presented in Figure 3B. As with Ae. aegypti, monthly ψ values were always high (harmonic mean 0.83; range, 0.66–0.94), with minimum ψ = 0.66 (95% confidence interval 0.59–0.72) in October 2011. Monthly Ae. albopictus ψ estimates were more unstable than those of Ae. aegypti, with relatively strong fluctuations after the dry-hot summer of 2009 (Figure 3B). Observed infestation (based on ovitraps) was also biased downwards (by ∼26%) in our Ae. albopictus data (Figure 4), yet ovitraps were fairly sensitive at detecting Ae. albopictus (p = 0.63, 95% confidence interval 0.62–0.65). Monthly local extinction probabilities were low: harmonic mean ε = 0.07, range 0.02–0.32, with the maximum value estimated for October 2011 (Figure 5B). Mean dwelling colonization probability was estimated as γ = 0.59 (95% confidence interval 0.51–0.66).

Discussion

Reliable dwelling infestation estimates are critical for decision-making in the context of dengue vector surveillance and control. The definition of programmatic goals, the management of resources, and the assessment of intervention effects all rely heavily upon such estimates. Using a large dataset and a sound analytical approach we have shown that routine vector surveillance and control are both performing poorly: at least in our study setting, (i) ‘rapid larval surveys’ yielded dwelling infestation indices that were markedly lower than the site-occupancy rates based on ovitrap data, and (ii) control campaigns had negligible effects on site-occupancy. Our results suggest that combining ovitrap-based surveillance (e.g., [43]–[45]) with analytical methods that account for imperfect detection (e.g., [21], [22], [38]–[40]) would help quantitatively assess, and likely enhance, dengue control programs. Moreover, from a disease transmission perspective, the presence of foraging gravid females in a dwelling, which ovitraps detect with reasonable sensitivity, is arguably more important than the presence of larvae, which is what ‘rapid larval surveys’ aim to detect.

Before discussing our findings any further, we identify several study limitations to keep in mind when interpreting the present results. First, we used detection/non-detection data, ignoring variations in vector abundance. However, presence-absence and abundance data seem to correlate well (e.g., refs. [46], [47]), and, importantly, both empirical and modeling results suggest that Ae. aegypti abundance thresholds (above which dengue transmission is maintained) are typically very low, i.e., ∼0.5–1 female per person or ∼0.5–1.5 pupae per person [48], [49]; therefore, the probability that at least one gravid female is present in a dwelling, which we modeled here, is clearly an epidemiologically relevant parameter. Second, some of our data may violate the assumption of independence of traps with regard to detection and of dwellings with regard to infestation; this may result in negatively biased ψ estimates with overly narrow confidence intervals. The high ψ estimates we found suggest that this problem was, in practice, negligible – our conclusions would not change because of somewhat broader confidence intervals. In addition, we measured only, and coarsely, a small number of covariates, but these were selected because of their known importance for our target species (e.g., [33]–[36], [50]). Our ‘control’ covariate included control interventions in just three out of 13 months of assessment, and this clearly lowered the precision of effect-size estimates: it seems possible that with more intervention data we might be able to detect a small effect whose 95% confidence interval could exclude zero. Yet, since ∼70–90% of dwellings remained infested despite control interventions, ‘statistical significance’ would in this case be of no practical importance [49]. Acknowledging these caveats, we feel confident that our models adequately estimate infestation rates as well as some of the major determinants of those rates in our study area. The main difference between our approach and previous attempts to assess infestation by dengue vectors is that we go beyond measuring indirect indices of infestation (i.e., adult or larvae presence/absence or counts) to produce statistical estimates of the probability that our study units (dwellings) are occupied by the target vector species.

We found little evidence that dwelling infestation rates decreased measurably as a result of the vector control campaigns carried out by local health authorities in our study neighborhood. These campaigns involved the elimination/treatment of thousands of artificial breeding containers [29], [30], and were expected to have larger effects on Ae. aegypti, which unlike Ae. albopictus rarely breeds in natural water collections [5], [17]. Our results show, however, no measurable effect of control interventions on any of the two vectors (Table 1); indeed, females of both species consistently continued to lay eggs, and probably forage, in most of the dwellings we surveyed, irrespective of whether control interventions had or had not taken place in the neighborhood. Our models suggest that this lack of effect could be related to the fact that interventions are usually planned to coincide with the wet-cool season, which is when local extinction probabilities drop to their lowest values (Figure 5). Summer interventions might perhaps be more effective [51], since they could synergize the negative effects of high temperatures on Ae. albopictus and Ae. aegypti detected by our models and in previous studies (e.g., [50]–[53]). These negative effects of high temperature, however, have to be considered in the particular context of our study. First, the extrinsic incubation period of dengue virus and the vector’s gonotrophic cycle can both be shortened by warmer weather, increasing transmission risk; second, relatively high temperatures probably favor vector development in overall cooler climates [10].

One practical implication of our findings is that Aedes breeding sites appear to be often overlooked by vector control agents during active surveillance and, principally, in control campaigns. This suggests a key drawback to be addressed in the development of novel Aedes control strategies, which should not heavily depend on the ability of control agents to detect breeding sites while inspecting premises. Two major candidate strategies address this problem from very different, but complementary, perspectives: (i) the use of adult mosquitoes to transfer potent larvicidal particles from contaminated ‘dissemination stations’ to clean breeding sites [54], and (ii) the release of mosquitoes carrying transgenes [55], [56] or specific Wolbachia strains [57] that impair reproduction and/or reduce competence to transmit dengue virus.

Conclusions

The reported bias of infestation indices that do not account for imperfect detection suggests that the findings of most dengue vector ecology studies must be interpreted with caution. Even ovitraps, which performed relatively well, yielded naïve infestation rates that were consistently biased downwards by about 30%. The analytical strategy we used here incorporates this sampling-process uncertainty, and could therefore substantially contribute to this field of inquiry.

Finally, our results suggest two promising avenues for the much-needed improvement of dengue vector surveillance [58]. First, simple hay infusion-baited ovitraps [32] should be preferred to ‘rapid larval surveys’: they are more sensitive and provide a measure of dwelling infestation by foraging gravid females (see also, e.g., refs. [43]–[45], [59]). Second, the repeated-sampling approach we used considerably improves infestation rate estimates by explicitly taking imperfect detection into account. Enhanced entomological surveillance systems and data analyses that explicitly account for the detection process would, in turn, allow for reliably assessing the effects of control interventions, irrespective of the specific tactics employed. Without such an assessment, the grounds on which massive public spending is directed towards dengue vector control (e.g., [60]) remain questionable.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

The complete sets of site-occupancy dynamic models.

https://doi.org/10.1371/journal.pone.0058420.s001

(XLSX)

Acknowledgments

We thank RM Mota, DLN Leite, and AC Barbosa for field assistance. RS Rocha and the Manaus municipal and state Health Departments provided data on routine control-surveillance activities. We particularly thank the residents who participated in the study – and patiently let us keep running our vector monitoring system. SA Juliano, LL Bailey, RE Gürtler, and BW Nelson provided insightful comments on earlier versions of this manuscript. This work is contribution number 17 of the Research Program on Infectious Disease Ecology in the Amazon (RP-IDEA) of the Instituto Leônidas e Maria Deane – Fiocruz Amazônia.

Author Contributions

Contributed to the last version of the manuscript: FA-F SDP-T GF SLBL EZ-P. Conceived and designed the experiments: FA-F GF SDP-T SLBL. Performed the experiments: SLBL SDP-T EZ-P. Analyzed the data: SDP-T GF FA-F. Contributed reagents/materials/analysis tools: FA-F GF SDP-T EZ-P SLBL. Wrote the paper: SDP-T FA-F GF.

Citation: Padilla-Torres SD, Ferraz G, Luz SLB, Zamora-Perea E, Abad-Franch F (2013) Modeling Dengue Vector Dynamics under Imperfect Detection: Three Years of Site-Occupancy by Aedes aegypti and Aedes albopictus in Urban Amazonia. PLoS ONE 8(3): e58420. https://doi.org/10.1371/journal.pone.0058420

References

1. WHO (2006) Scientific Working Group Report on Dengue. Available: http://www.who.int/tdr/publications/documents/swg_dengue_2.pdf. Accessed 2012 February 14.

2. WHO (2012) Global Alert and Response: Impact of Dengue. Available: http://www.who.int/csr/disease/dengue/impact/en/. Accessed 2012 February 14.

3. Guzmán MG, Halstead SB, Artsob H, Buchy P, Farrar J, et al. (2010) Dengue: a continuing global threat. Nat Rev Microbiol 8(12 Suppl.): S7–S16.

4. Guzmán A, Istúriz RE (2010) Update on the global spread of dengue. Int J Antimicrob Agents (Suppl. 1): S40–S42.

5. Gubler DJ (1998) Dengue and dengue hemorrhagic fever. Clin Microbiol Rev 11(3): 480–496.

6. Heintze C, Garrido MV, Kroeger A (2007) What do community-based dengue control programmes achieve? A systematic review of published evaluations. Trans R Soc Trop Med Hyg 101(4): 317–325.

7. Ballenger-Browning KK, Elder JP (2009) Multi-modal Aedes aegypti mosquito reduction interventions and dengue fever prevention. Trop Med Int Health 14(12): 1542–1551.

8. Gürtler RE, Garelli FM, Coto HD (2009) Effects of a five-year citywide intervention program to control Aedes aegypti and prevent dengue outbreaks in northern Argentina. PLoS Negl Trop Dis 3(4): e427.

9. Esu E, Lenhart A, Smith L, Horstick O (2010) Effectiveness of peridomestic space spraying with insecticide on dengue transmission; systematic review. Trop Med Int Health 15(5): 619–631.

10. Kyle JL, Harris E (2008) Global spread and persistence of dengue. Annu Rev Microbiol 62: 71–92.

11. San Martín JL, Brathwaite O, Zambrano B, Solórzano JO, Bouckenooghe A, et al. (2010) The epidemiology of dengue in the Americas over the last three decades: a worrisome reality. Am J Trop Med Hyg 82(1): 128–135.

12. Gubler DJ (1988) Dengue. In: Monath TP, editor. Epidemiology of Arthropod-Borne Viral Diseases. Boca Raton: CRC Press, Inc. 223–260.

13. Reiter P (2007) Oviposition, dispersal, and survival in Aedes aegypti: implications for the efficacy of control strategies. Vector Borne Zoonotic Dis 7(2): 261–273.

14. Gubler DJ (2002) Epidemic dengue/dengue hemorrhagic fever as a public health, social and economic problem in the 21st century. Trends Microbiol 10(2): 100–103.

15. Gonçalves da Silva A, Cunha ICL, Santos WS, Luz SLB, Ribolla PEM, et al. (2012) Gene flow networks among American Aedes aegypti populations. Evol Appl 5(7): 664–676.

16. Gubler DJ (2002) The global emergence/resurgence of arboviral diseases as public health problems. Arch Med Res 33(4): 330–342.

17. Gratz NG (2004) Critical review of the vector status of Aedes albopictus. Med Vet Entomol 18(3): 215–227.

18. Lambrechts L, Scott TW, Gubler DJ (2010) Consequences of the expanding global distribution of Aedes albopictus for dengue virus transmission. PLoS Negl Trop Dis 4(5): e646.

19. Reiter P, Gubler DJ (1997) Surveillance and control of urban dengue vectors. In: Gubler DJ, Kuno G, editors. Dengue and Dengue Hemorrhagic Fever. London: CAB International. 425–462.

20. FUNASA (2002) Programa Nacional de Controle da Dengue. Brasília: Ministério da Saúde/Fundação Nacional de Saúde.

21. Mackenzie DI (2005) Was it there? Dealing with imperfect detection for species presence/absence data. Aust NZ J Stat 47(1): 65–74.

22. Abad-Franch F, Ferraz G, Campos C, Palomeque FS, Grijalva MJ, et al. (2010) Modeling disease vector occurrence when detection is imperfect: infestation of Amazonian palm trees by triatomine bugs at three spatial scales. PLoS Negl Trop Dis 4(3): e620.

23. Figueiredo RMP, Thatcher BD, Lima ML, Almeida TC, Alecrim WD, et al. (2004) Doenças exantemáticas e primeira epidemia de dengue ocorrida em Manaus, Amazonas, no período de 1998–1999. Rev Soc Bras Med Trop 37(6): 476–479.

24. Ríos-Velásquez CM, Codeço CT, Honório NA, Sabroza PS, Moresco M, et al. (2007) Distribution of dengue vectors in neighborhoods with different urbanization types of Manaus, state of Amazonas, Brazil. Mem Inst Oswaldo Cruz 102(5): 617–623.

25. Fé NF, Barbosa MGV, Alecrim WD, Guerra MVF (2003) Registro da ocorrência de Aedes albopictus em área urbana do município de Manaus, Amazonas. Rev Saude Publica 37(5): 674–675.

26. Figueiredo RMP, Naveca FG, Souza MB, Melo M, Viana SS, et al. (2008) Dengue virus type 4, Manaus, Brazil. Emerg Infect Dis 14(4): 667–669.

27. FUNASA (2001) Dengue: Instruções para Pessoal de Combate ao Vetor. Manual de Normas Técnicas, 3a ed. Brasília: Ministério da Saúde/Fundação Nacional de Saúde.

28. Ministério da Saúde (2005) Diagnóstico Rápido nos Municípios para Vigilância Entomológica do Aedes aegypti no Brasil – LIRAa. Metodologia para Avaliação dos Índices Breteau e Predial. Brasília: Ministério da Saúde.

29. Fundação de Vigilância em Saúde do Estado do Amazonas (2008) Relatório Final: Operação Impacto – Controle da Dengue. Manaus: Fundação de Vigilância em Saúde do Estado do Amazonas.

30. Fundação de Vigilância em Saúde do Estado do Amazonas (2009) Relatório Final: Operação Impacto II – Controle da Dengue. Manaus: Fundação de Vigilância em Saúde do Estado do Amazonas.

31. Donatti JE, Gomes AC (2007) Adultrap: descrição de armadilha para adultos de Aedes aegypti (Diptera, Culicidae). Rev Bras Entomol 51(2): 255–256.

32. Reiter P, Amador MA, Colon N (1991) Enhancement of the CDC ovitrap with hay infusions for daily monitoring of Aedes aegypti populations. J Am Mosquito Contr Assoc 7(1): 52–55.

33. Scott TW, Morrison AC, Lorenz LH, Clark GG, Strickman D, et al. (2000) Longitudinal studies of Aedes aegypti (Diptera: Culicidae) in Thailand and Puerto Rico: population dynamics. J Med Entomol 37(1): 77–88.

34. Alto BW, Juliano SA (2001) Temperature effects on the dynamics of Aedes albopictus (Diptera: Culicidae) populations in the laboratory. J Med Entomol 38(4): 548–556.

35. Favier C, Degallier N, Vilarinhos PDTR, de Carvalho MDSL, Yoshizawa MAC, et al. (2006) Effects of climate and different management strategies on Aedes aegypti breeding sites: a longitudinal survey in Brasília (DF, Brazil). Trop Med Int Health 11(7): 1104–1118.

36. Tun-Lin W, Kay BH, Barnes A (1995) The premise condition index: a tool for streamlining surveys of Aedes aegypti. Am J Trop Med Hyg 53(6): 591–594.

37. Tukey JW (1977) Exploratory Data Analysis. Reading: Addison-Wesley.

38. MacKenzie DI, Nichols JD, Lachman GB, Droege S, Royle JA, et al. (2002) Estimating site occupancy rates when detection probabilities are less than one. Ecology 83(8): 2248–2255.

39. MacKenzie DI, Nichols JD, Hines JE, Knutson MG, Franklin AB (2003) Estimating site occupancy, colonization, and local extinction when a species is detected imperfectly. Ecology 84(8): 2200–2207.

40. MacKenzie DI, Nichols JD, Royle JA, Pollock KH, Bailey LL, et al. (2006) Occupancy Estimation and Modeling: Inferring Patterns and Dynamics of Species Occurrence. San Diego: Elsevier Academic Press.

41. Burnham KP, Anderson DR (2001) Kullback-Leibler information as a basis for strong inference in ecological studies. Wildlife Res 28(2): 111–119.

42. Hines JE (2006) PRESENCE 4.0 Software to estimate patch occupancy and related parameters. USGS-PWRC. http://www.mbr-pwrc.usgs.gov/software/presence.shtml.

43. Furlow BM, Young WW (1970) Larval surveys compared to ovitrap surveys for detecting Aedes aegypti and Aedes triseriatus. Mosq News 30(3): 468–470.

44. Chadee DD (1986) A comparison of three Aedes aegypti sampling methods in Trinidad, West-Indies. Cah ORSTOM Ser Ent Med Parasitol 24(3): 199–205.

45. Lourenço-de-Oliveira R, Lima JBP, Peres R, Alves FDC, Eiras AE, et al. (2008) Comparison of different uses of adult traps and ovitraps for assessing dengue vector infestation in endemic areas. J Am Mosquito Contr Assoc 24(3): 387–392.

46. Mogi M, Choochote W, Khamboonruang C, Suwanpanit P (1990) Applicability of presence-absence and sequential sampling for ovitrap surveillance of Aedes (Diptera: Culicidae) in Chiang Mai, Northern Thailand. J Med Entomol 27(4): 509–514.

47. Tun-Lin W, Kay BH, Barnes A, Forsyth S (1996) Critical examination of Aedes aegypti indices: correlations with abundance. Am J Trop Med Hyg 54(5): 543–547.

48. Focks DA, Brenner RJ, Hayes J, Daniels E (2000) Transmission thresholds for dengue in terms of Aedes aegypti pupae per person with discussion of their utility in source reduction efforts. Am J Trop Med Hyg 62(1): 11–18.

49. Medeiros LCdC, Castilho CAR, Braga C, de Souza WV, Regis L, et al. (2011) Modeling the dynamic transmission of dengue fever: investigating disease persistence. PLoS Negl Trop Dis 5(1): e942.

50. Barrera R, Amador M, Mackay A J (2011) Population dynamics of Aedes aegypti and dengue as influenced by weather and human behavior in San Juan, Puerto Rico. PLoS Negl Trop Dis 5(12): e1378.

51. Chadee DD (2009) Impact of pre-seasonal focal treatment on population densities of the mosquito Aedes aegypti in Trinidad, West Indies: A preliminary study. Acta Trop 109(3): 236–240.

52. Juliano SA, O’Meara GF, Morrill JR, Cutwa MM (2002) Desiccation and thermal tolerance of eggs and the coexistence of competing mosquitoes. Oecologia 130(3): 458–469.

53. Leisnham PT, Juliano SA (2009) Spatial and temporal patterns of coexistence between competing Aedes mosquitoes in urban Florida. Oecologia 160(2): 343–352.

54. Devine GJ, Zamora Perea E, Killeen GF, Stancil JD, Clark SJ, et al. (2009) Using adult mosquitoes to transfer insecticides to Aedes aegypti larval habitats. Proc Natl Acad Sci USA 106(28): 11530–11534.

55. Wise de Valdez MR, Nimmo D, Betz J, Gong HF, James AA, et al. (2011) Genetic elimination of dengue vector mosquitoes. Proc Natl Acad Sci USA 108(12): 4772–4775.

56. Harris AF, Nimmo D, McKemey AR, Kelly N, Scaife S, et al. (2011) Field performance of engineered male mosquitoes. Nat Biotechnol 29(11): 1034–1037.

57. Hoffmann AA, Montgomery BL, Popovici J, Iturbe-Ormaetxe I, Johnson PH, et al. (2011) Successful establishment of Wolbachia in Aedes populations to suppress dengue transmission. Nature 476(7361): 454–457.

58. Morrison AC, Zielinski-Gutierrez E, Scott TW, Rosenberg R (2008) Defining challenges and proposing solutions for control of the virus vector Aedes aegypti. PLoS Med 5(3): e68.

59. Pilger D, Lenhart A, Manrique-Saide P, Siqueira JB, Rocha WT, et al. (2011) Is routine dengue vector surveillance in central Brazil able to accurately monitor the Aedes aegypti population? Results from a pupal productivity survey. Trop Med Int Health 16(9): 1143–1150.

60. Taliberti H, Zucchi P (2010) Custos diretos do programa de prevenção e controle da dengue no Município de São Paulo em 2005. Rev Panam Salud Publica 27(3): 175–180.