 Research
 Open Access
 Published:
Estimating absolute indoor density of Aedes aegypti using removal sampling
Parasites & Vectors volume 12, Article number: 250 (2019)
Abstract
Background
Quantification of adult Aedes aegypti abundance indoors has relied on estimates of relative density (e.g. number of adults per unit of sampling or time), most commonly using traps or timed collections using aspirators. The lack of estimates of the sensitivity of collections and lack of a numerical association between relative and the absolute density of adult Ae. aegypti represent a significant gap in vector surveillance. Here, we describe the use of sequential removal sampling to estimate absolute numbers of indoor resting Ae. aegypti and to calculate calibration coefficients for timed Prokopack aspirator collections in the city of Merida, Yucatan State, Mexico. The study was performed in 200 houses that were selected based on recent occurrence of Aedesborne viral illness in residents. Removal sampling occurred in 10minute sampling rounds performed sequentially until no Ae. aegypti adult was collected for 3 hours or over 2 consecutive 10minute periods.
Results
A total of 3439 Ae. aegypti were collected. The sensitivity of detection of positive houses in the first sampling round was 82.5% for any adult Ae. aegypti, 78.5% for females, 75.5% for males and 73.3% for bloodfed females. The total number of Ae. aegypti per house was on average ~5 times higher than numbers collected for the first sampling round. There was a positive linear relationship between the relative density of Ae. aegypti collected during the first 10min round and the absolute density for all adult metrics. Coefficients from the linear regression were used to calibrate numbers from 10min collections into estimates of absolute indoor Ae. aegypti density for all adults, females and males.
Conclusions
Exhaustive removal sampling represents a promising method for quantification of absolute indoor Ae. aegypti density, leading to improved entomological estimates of mosquito distribution, a key measure in the assessments of the risk pathogen transmission, disease modeling and the evaluation of vector control interventions.
Background
If all individuals in a population cannot be counted, they must be sampled [1]. This basic principle represents the cornerstone of ecological field studies and provides the basis for estimating the relative or absolute numbers of individuals within a population or community [1, 2]. Many factors, including ecological, economic and statistical, influence the methodologies employed to estimate population size, which broadly include markrecapture, catchperuniteffort (CPUE, e.g. trapping, timed collections), removal sampling, distance sampling and quadrat methods [1, 2]. When populations have to be monitored regularly or at multiple locations, simpler (and often less expensive) and potentially less accurate methods such as CPUE are generally preferred to the more costly markrecapture approaches [1]. The use of CPUE is common practice in medical entomology [3]. Some examples include the use of passive or active traps (e.g. CDC light traps, ovitraps, adult resting boxes), timed collections (e.g. adult mosquito aspiration, kissing bug timed manual collections) and immature habitat sampling (e.g. Aedes pupal surveys, tick dragging). Irrespective of the method, all such approaches provide a measure of relative abundance (also called density), in which the number of collected individuals is a function of the time or effort employed to collect them. Such estimates are prone to bias for multiple reasons, including differences between the collectors themselves, heterogeneity of captures across space, and sensitivity of the sampling methods in situations of low population abundance [1, 2]. More often than not, such biases are not appropriately quantified, which limits the validity of estimates of relative density, particularly in situations in which vector abundance is low or spatially heterogeneous. In order for a sampling methodology to be robust and valid, measures of the association between the relative sample and the absolute density are desired [1, 2, 4, 5].
Aedes aegypti rest primarily indoors [6], where they frequently and preferentially bite humans [7]. Compared to other vector mosquitoes (e.g. Culex quinquefasciatus), Ae. aegypti is considered a low abundance species [8]. Multiple approaches have been implemented to estimate Ae. aegypti density and population size, with markrecapture experiments [9,10,11], the innovative use of the ratio of wildtype to Wolbachiainfected released mosquitoes [12, 13] or statistical/mathematical models fitted to field data [14, 15] as the most commonly used. Markrecapture (whether with dust or Wolbachia) studies point to an average of 5–10 females per premise as the density of Ae. aegypti during peaktransmission periods [11,12,13]. Extensive indoor aspiration sampling performed in Iquitos, Peru, collected <10 adults per house, on average, during periods of high virus transmission [16]. In Cairns, Australia, Williams et al. [14] used pupal productivity data fitted to a mathematical model of Ae. aegypti productivity to estimate an average of 4–23 female Ae. aegypti per premise during the period of arbovirus transmission (wet season). All such approaches relied on sampling methods such as BG sentinel traps, which may also be prone to error. For instance, using a mini markrelease study design, Johnson et al. [17] estimated that the BG sentinel trap captures ~20–30% of adults outdoors, providing a measure of the sensitivity or sampling rate for the traps. This is, to our knowledge, the only published attempt of calibration of relative to absolute abundance for any adult Ae. aegypti sampling method.
Timed adult aspiration (the use of motorized vacuums by an operator to capture resting and flying mosquitoes indoors during a defined period) is considered a gold standard for indoor adult Ae. aegypti sampling [8, 18]. Used primarily in research, adult aspiration provides a relatively unbiased measure of mosquito relative abundance (i.e. collecting both males and females as well as fed and unfed females) when conducted for ~10minutes per house [19,20,21,22]. This rapid measure of relative density may be prone to collector variability, or to differential capture rate as a function of vector density or the size and complexity of premises. Given no study has yet quantified the sensitivity of sampled and absolute estimates of Ae. aegypti density using aspiration devices, the magnitude of such potential sources of bias is unknown. To address this gap, we performed a field study to calibrate Prokopack [22] aspirator collections using correction coefficients derived from comparing sampled to absolute density estimated from sequential removal sampling. Removal sampling, where the sequential removal of individuals from the population using constant effort leads to a reduction in the catch per unit effort, can be used to estimate total population size [1]. The method assumes that: (i) the population is closed; (ii) the probability of each individual being caught is constant; and (iii) all individuals have the same probability of being collected on any given sample [1]. The original method involved fitting simple maximum likelihood regression models to data on the number caught on the t occasion versus the total catch up to occasion t − 1, allowing estimating initial population size when t = 0 [1]. However, if the population is depleted through the sampling procedure, total catch (rather than linear regression estimates) can be used to estimate absolute abundance. The removal method has been useful for estimating population size for fisheries (e.g. [23,24,25]) and was applied to estimate the population size of Ae. aegypti in Kenya [9], Anopheles spp. in Pakistan [26] and the efficiency of Ixodes scapularis population size and dragging efficiency [27], but not yet to associate indices of relative and absolute mosquito density.
Methods
The study was performed in the city of Merida (population ~1 million), Yucatan State, Mexico. Merida is highly endemic for dengue [28] and other Aedesborne viruses [29], which in the city are transmitted solely by Ae. aegypti [29, 30]. A total of 200 houses located in an Aedesborne disease transmission hotspot area of Merida [29] were included in this study. Only houses with a recent (within 1 month) occurrence of a symptomatic case of dengue, chikungunya or Zika (based on information from the local Ministry of Health, MOH) were included. The study ran for two transmission seasons to achieve the desired target number of 200 houses. Indoor mosquito collections were performed during the period of most intense virus transmission (June–December). After obtaining informed consent from householders, exhaustive adult mosquito collections using Prokopack aspirators [22] were conducted using removal sampling. Mosquitoes were sequentially collected from each house using constant effort and at predefined intervals. Each house was visited by a team of three entomologists trained in the collection of adult mosquitoes using the Prokopack aspirator. The removal sampling collection sequence was as follows: one of the entomologists entered the house and collected mosquitoes resting in all rooms, including the kitchen and bathroom, for a period of 10 min. A timer was used to make sure collections ended at the 10min mark, at which point the collector exited the house and gave the aspirator collection cup to a supervisor, who immediately sent a second collector inside the house to perform another 10min collection round. This collection sequence continued for 3 h or until no Ae. aegypti were collected for two consecutive rounds, whichever occurred first. The supervisor was the only team member who knew whether Ae. aegypti was collected during each sampling round. This procedure blinded collectors and provided less opportunity for bias. We chose 10 min as the sampling time because it is the average sampling duration of a standard urban home in Merida [30] and elsewhere [16]. As collection proceeded, each aspirator collection cup was labeled with the house code and the collection sequence number. The collectors alternated who performed the first 10min sampling at each house. After all entomological collections were complete, a survey was administered to the household head asking for basic information about the home (size, number of residents, presence of screens, etc). All collected mosquitoes were transferred to the laboratory, where they were knocked down in a −20 °C freezer for 10 min, and then sexed and identified to species following standard keys. Female Ae. aegypti were classified by their engorgement status following a categorical score [31] and then dissected to separate their head from the rest of the body for future virus testing.
Data analysis
We focused analyses on measures of total Ae. aegypti adults, Ae. aegypti females, bloodfed Ae. aegypti females and Ae. aegypti males. A binomial generalized additive mixed model (GAMM) was fitted to the presence of Ae. aegypti in each house for each 10min collection round. The model included infestation on each sampling round (presence = 1, absence = 0), with sampling round as a nonlinear term (continuous variable, set in 10min increments) and house ID as random intercept.
Absolute density per house was calculated in two ways: as the sum of Ae. aegypti collected across all sampling rounds (named total catch) or following the removal sampling equation method developed by Carle and Strub [16]. The regression method, also called maximum weighted likelihood method, fits a regression line to the number of mosquitoes caught on the ith sample, as a function of the total catch up to ith − 1. The equations of the method are: \(M = \mathop \sum \nolimits_{i = 1}^{k} \left( {k  i} \right)_{{u_{i} }}\) where k is the number of samples taken and u_{i} the number of animals caught in the ith sample. The population size, N, is calculated as the smallest integer greater than the total catch, T, that satisfies the following inequality: \(\left( {\frac{N + 1}{N  T + 1}} \right)\left( {\frac{kN  M  T + 0.5k}{kn  M + 1 + 0.5k}} \right)^{k} \le 1\). The standard error in this inequality is calculated using maximum likelihood [32]. We used the removal function from the R package FSA [33] to calculate N using the Carle and Strub method. We performed simple linear regression analysis to compare N (regression method) to the total catch, and assess the relevance of both approaches to estimate absolute density.
MannWhitney test was used to evaluate the difference in capture rate or sensitivity (percentage of Ae. aegypti collected in first 10min collection round divided by the total catch across two levels of vector abundance; low was defined as ≤10 Ae. aegypti/house, high was defined as >10 Ae. Aegypti/house). Simple linear regressions calibrated the relative abundance to absolute abundance values. Maximum likelihood was used to fit the regression equation to the data. To assess the effect of any household characteristics influencing model fit, we performed a multiple linear regression including variables such as household size, presence of mosquito screens and recent use of insecticides. All analyses were performed within the R programing environment (https://www.rproject.org/) and GAMMs were run using the lme4 package [34].
Results
A total of 3439 Ae. aegypti (female to male ratio, 1.7:1) were collected in 200 houses employing a total catch effort of 269 h (a total of 1615 10min collection rounds). All mosquitoes were collected within 16 10min sampling rounds (up to 160 min). Most houses were infested with Ae. aegypti adults in at least one 10min sampling round (n = 179, 89.5%), whereas 84.5% were infested with Ae. aegypti females, 71.5% with Ae. aegypti males and 79.5% with bloodfed Ae. aegypti females in at least one round. When analyzed by house, the probability of detecting a positive house (derived from a binomial GAMM) decreased as the population was sampled in sequential rounds for all metrics (Fig. 1), confirming the utility of sequential sampling in detecting finite populations inside houses. Sensitivity of detection of positive houses in the first sampling round was 82.5% for adults, 78.5% for females, 75.5% for males and 73.3% for bloodfed females (Fig. 2). Cumulative sensitivity increased slightly with increases in capture effort, reaching an asymptote at 40 min for any adults, 50 min for males and 60 min for females (Fig. 2). Aggregating data from the first two sampling rounds (i.e. equivalent to performing a 20min collection) was associated with a modest increase in sensitivity (+10.5% for any adults, +9.5% for females, +16.5% for males and 13% for bloodfed females) (Fig. 2).
The association between the total catch (sum of Ae. aegypti collected across all sampling rounds) and the number of Ae. aegypti estimated by the regression method is shown in Fig. 3. Across all levels of Ae. aegypti density, both methods provided nearly equivalent results (Fig. 3). A linear regression showed that absolute adult Ae. aegypti density estimates from total catch and the regression method are significantly and positively associated (beta = 0.826; standard error of beta = 0.03; tvalue = 31.6, P < 0.001, R^{2} = 0.85). As such, and for simplicity, the remainder of the manuscript will utilize total catch estimates for calibrating relative and absolute densities. Total catch was skewed (Fig. 4), averaging a total number per positive house of 19.3 (range, 1–244) Ae. aegypti adults, 12.9 (1–169) females, 14.3 (1–175) males and 11 (1–159) bloodfed females. These total catch averages were ~5 times higher than what was captured during the first sampling round, which was 4.4 (1–112) adults, 3.2 (1–86) females, 1.8 (1–49) males and 2.3 (1–65) bloodfed females. The absolute density of Ae. aegypti females and bloodfed females per house were strongly and significantly associated (linear regression, F_{(1,176)} = 10188, P < 0.0001, R^{2} = 0.98). Of the total Ae. aegypti collected per house, a median of 26% were captured during the first 10min round in low density houses and 24% in high density houses (Fig. 5). While the distribution of catchability values was wider for low than for high density houses (Fig. 5), the median did not vary statistically between the two density strata for Ae. aegypti adults (twosample Wilcoxon, W = 3568, P = 0.2379), females (W = 2932, P = 0.2532) or males (W = 1870, P = 0.7341) (Fig. 4). Differences among houses in percent catchability were not explained by neither the percentage of doors nor windows with mosquito screens (Gaussian generalized linear model, beta = −0.15, z = −0.368, P = 0.71), the area of the house (beta = 0.003, z = −0.73, P = 0.47) the number of residents in the house (beta = −0.07, z = −0.90, P = 0.38) and whether the house had low (≤10/house) or high (>10/house) absolute Ae. aegypti density (beta = −0.49, z = −1.42, P = 0.16). Given the lack of significant association between abundance and capture rate, further analyses considered abundance as a continuous variable.
There was a significant positive linear relationship between the relative abundance of Ae. aegypti collected in the first 10min round and the total catch for all adult metrics (Table 1, Fig. 6). The strong association between relative and absolute abundance allowed for the estimation of calibration coefficients from simple linear regressions (Table 1). Given the presence of a high capture location where 243 adults were collected, we implemented regression models including and excluding this high infestation location. Excluding the highest collection point from the regressions (Fig. 6, right panel) reduced model fit from an R^{2} of 0.73 to an R^{2} of 0.62 for adults and 0.76 to 0.59 for females while R^{2} for males remained unchanged. The intercept in the models including all data points ranged from 7.1 ± 1.2 for adults, 5.5 ± 0.8 for females and 3.3 ± 0.8 for males (Table 1). These metrics are indicative of the sampling error associated with a single 10min round of adult collection and allow for adjusting, through the regression equation, the capture rate of indoor Ae. aegypti when at least one mosquito is collected within the initial 10min collection.
Discussion
While entomological sampling is a cornerstone of vectorborne disease research and surveillance, there is a paucity of information about the sensitivity of existing sampling methods. Our research provides detailed estimates of the sensitivity of using Prokopack adult mosquito aspirators to collect indoor Ae. aegypti. Our estimates show that, regardless of the absolute abundance and housing characteristics, a 10minute collection detects ~80% of houses that are infested with Ae. aegypti and captures 24–26% of all Ae. aegypti. Our data was applied to a simple mathematical equation to provide the basis for future calibrations of relative density obtained from 10minute collections into absolute abundance estimates of indoor Ae. aegypti.
While relative measures of vector density can be informative and allow for comparisons between different time points and locations, they can also suffer from collector or sampling errors and are prone to spatial and temporal bias [1,2,3]. Identifying a functional relationship between estimates of relative and absolute density has been common practice in agriculture [4] and wildlife research [5]. To date, the advantage and power of quantifying absolute population abundance has not been widely employed in the field of medical entomology. For Ae. aegypti, the only example of calibration is the one performed by Johnson et al. [17] using a mini markrelease study to calibrate BG sentinel traps. Our study introduces a simpler methodology, sequential removal sampling, to reliably estimate calibration coefficients to quantify absolute Ae. aegypti density indoors. This methodology has the potential to be expanded to other disease vectors readily collected using Prokopack aspirators such as endophilic Anopheles spp. [19,20,21] and Culex spp. [22]. While our study shows the benefit of sequential removal sampling to quantify indoor absolute density, the applicability of this methodology for outdoor Ae. aegypti collections will need to be thoroughly assessed. A study in Cairns, Australia, found that BGSentinel traps collected more Ae. aegypti than CDC Backpack aspirators outdoors [35]. As such, Prokopack removal sampling outdoors will have to be implemented at a much higher effort than indoors, and probably with lower collection efficiency compared with traps. Given the expected low yield per sample, implementation of the regression method for estimating absolute density would be a necessary step, provided sequential collections lead to a reduction in the catch per unit effort over time [1].
The finding of a decay in catchability with each sampling round (Fig. 1) supports a key removal sampling assumption: populations are closed at the time of sampling. Using simple linear regressions, we estimated correction coefficients, which allowed calculation of absolute Ae. aegypti numbers indoors at the time of collection (together with 95% confidence intervals). We generated two fits to the data, one including and another excluding a house where more than 240 Ae. aegypti were collected. The model fit the data better when using all data points, but we suggest that our equations could be applied as follows: for houses where fewer than 30 total Ae. aegypti are collected during the tenminute collection round, the equations excluding the extreme data point should be used. For houses where more than a total of 30 Ae. aegypti are collected, the equation including all data points should be used. In either occasion, houses with no Ae. aegypti should be excluded from the calculation, as including them would lead to estimates of density equal to the intercept of the model. We believe that applying the regression equation estimated from our study would allow a more precise estimation of absolute abundance from a single 10minute collection round. As similar sequential removal studies are performed in other locations (ideally with different housing characteristics and vector densities), the information generated will help refine model equations and increase the value of calibrated adult Ae. aegypti indices.
Using absolute (rather than relative) density estimates has the potential to change the way entomological data are interpreted when evaluating public health interventions. As an example, we used the calibration coefficient generated by this study to calculate total indoor adult Ae. aegypti absolute abundance at baseline and at four sampling time points after applying indoor residual spraying (Fig. 7, data from [36]). Using the estimates of absolute abundance provided a much more dramatic measure of the absolute impact of indoor residual spraying, with an estimated reduction from ~10.0 to ~2.0 Ae. aegypti per house, compared to 2.0 to 0.25 Ae. aegypti per house when using the relative abundances before calibration (Fig. 7). While the relative change is still proportional among both sets of data, the absolute impact may be more relevant at the time of identifying entomological thresholds of protection or intervention impact. Rearandrelease methods of vector population suppression or modification (e.g. sterile insect technique, Wolbachiarelated interventions, transgenic mosquitoes) would also benefit from more accurate estimates of absolute mosquito abundance. As mosquito releases are informed by the relative density of adults sampled (e.g. 10 times more adults than the ones collected in weekly BGsentinel trapping [37]), estimates of absolute population density are key for successful rearandrelease implementation. Our estimates indicate that houses can have up to five times more adult Ae. aegypti than estimated with a 10minute collection round. As such, estimates of released mosquitoes based on relative density would significantly underestimate the effort needed for successful population control or modification.
Compared to other species, Ae. aegypti is considered a low abundance vector [8]. This assessment has been inferred from the low numbers of adult Ae. aegypti generally collected indoors using aspirators, with collections during high density periods averaging four to five females per house [11, 16]. Our estimates of absolute indoor Ae. aegypti density show a different picture. On average, ~13 female Ae. aegypti were collected per house, with a highly overdispersed distribution that, in our sample of 200 houses, reached a maximum of 169 females. In fact, 21.4% of Ae. aegypti positive houses had 20 or more females indoors. In Merida, the average household size is four people. Assuming equal biting probability, we can hypothesize that in 21% of houses with Ae. aegypti, there are a minimum of five Ae. aegypti females per person. Mathematical models have been a useful tool for linking such entomological measures to virus transmission estimates. However, most models have relied on relative density data to parameterize mosquito abundance [38]. As more reliable estimates of Ae. aegypti absolute density are generated, there will be a need to recalculate arbovirus transmission risk under scenarios that, based on our estimates, could include up to five times more Ae. aegypti females than previously assumed. Whether these updated estimates of vector abundance will result in profound epidemiological changes will have to be further explored.
For mosquitoborne pathogens, there is a renewed interest in comprehensively evaluating the epidemiological impact of vector control [39, 40]. An emphasis on epidemiological endpoints for evaluating interventions has emerged from decades of research pointing to a poor association between entomological measures and vectorborne disease transmission risk [39, 41]. Particularly for Aedesborne viruses, poor estimates of Ae. aegypti adult abundance, combined with the focal and local nature of virus transmission [42, 43] have challenged the finding of an association between entomologic data and epidemiological risk [8]. Our findings suggest that these poor correlations may also be the result of a lack of accuracy in estimating the absolute abundance of Ae. aegypti. As such, we recommend a greater emphasis on applying methods to estimate total adult Ae. aegypti density in future entomological studies and vector control evaluations.
Conclusions
While myriad sampling methods are used to quantify the density of vectors (e.g. traps, active collections, landing counts), estimates of their sensitivity are rare. This knowledge gap impacts the reliability of entomological estimates, and limits the quantification of the entomological impact of vector control. In this article, we introduce a methodology for rapidly estimating the absolute density of vectors indoors, and use it to estimate the sensitivity of the Prokopack mosquito collector on Aedes aegypti populations from Merida, Mexico.
Availability of data and materials
Data supporting the conclusions of this article are included within the article. The datasets used and/or analyzed during the present study are available from the corresponding author upon reasonable request.
Abbreviations
 CPUE:

catch per unit effort
 MOH:

Ministry of Health
References
 1.
Krebs CJ. Ecological methodology. Menlo Park: Benjamin/Cummings; 1999.
 2.
Henderson PA, Southwood TRE. Ecological methods. Hoboken: Wiley; 2016.
 3.
Silver JB. Mosquito ecology: field sampling methods. Dordrecht: Springer Netherlands; 2007.
 4.
Pedigo LP, Buntin GD. Handbook of sampling methods for arthropods in agriculture. Boca Raton: Taylor & Francis; 1993.
 5.
Falcy MR, McCormick JL, Miller SA. Proxies in practice: calibration and validation of multiple indices of animal abundance. J Fish Wild Manag. 2016;7:117–28.
 6.
DzulManzanilla F, IbarraLopez J, Bibiano Marin W, MartiniJaimes A, Leyva JT, CorreaMorales F, et al. Indoor resting behavior of Aedes aegypti (Diptera: Culicidae) in Acapulco, Mexico. J Med Entomol. 2017;54:501–4.
 7.
Liebman KA, Stoddard ST, Reiner RC Jr, Perkins TA, Astete H, Sihuincha M, et al. Determinants of heterogeneous blood feeding patterns by Aedes aegypti in Iquitos, Peru. PLoS Negl Trop Dis. 2014;8:e2702.
 8.
Scott TW, Morrison AC, Takken W. Aedes aegypti density and the risk of denguevirus transmission. In: Takken W, Scott TW, editors. Ecological aspects for application of genetically modified mosquitoes. Frontis: Wageningen; 2003.
 9.
Trpis M, Hausermann W, Craig GB Jr. Estimates of population size, dispersal, and longevity of domestic Aedes aegypti aegypti (Diptera: Culicidae) by markreleaserecapture in the village of Shauri Moyo in eastern Kenya. J Med Entomol. 1995;32:27–33.
 10.
Neira M, Lacroix R, Caceres L, Kaiser PE, Young J, Pineda L, et al. Estimation of Aedes aegypti (Diptera: Culicidae) population size and adult male survival in an urban area in Panama. Mem Inst Oswaldo Cruz. 2014;109:879–86.
 11.
Villela DA, Codeco CT, Figueiredo F, Garcia GA, MacieldeFreitas R, Struchiner CJ. A Bayesian hierarchical model for estimation of abundance and spatial density of Aedes aegypti. PLoS ONE. 2015;10:e0123794.
 12.
Garcia Gde A, Dos Santos LM, Villela DA, MacieldeFreitas R. Using Wolbachia releases to estimate Aedes aegypti (Diptera: Culicidae) population size and survival. PLoS ONE. 2016;11:e0160196.
 13.
Ritchie SA, Montgomery BL, Hoffmann AA. Novel estimates of Aedes aegypti (Diptera: Culicidae) population size and adult survival based on Wolbachia releases. J Med Entomol. 2013;50:624–31.
 14.
Williams CR, Johnson PH, Ball TS, Ritchie SA. Productivity and population density estimates of the dengue vector mosquito Aedes aegypti (Stegomyia aegypti) in Australia. Med Vet Entomol. 2013;27:313–22.
 15.
Rodrigues Mde M, Marques GR, Serpa LL, Arduino Mde B, Voltolini JC, Barbosa GL, et al. Density of Aedes aegypti and Aedes albopictus and its association with number of residents and meteorological variables in the home environment of dengue endemic area, Sao Paulo, Brazil. Parasites Vectors. 2015;8:115.
 16.
LaCon G, Morrison AC, Astete H, Stoddard ST, PazSoldan VA, Elder JP, et al. Shifting patterns of Aedes aegypti fine scale spatial clustering in Iquitos, Peru. PLoS Negl Trop Dis. 2014;8:e3038.
 17.
Johnson PH, Spitzauer V, Ritchie SA. Field sampling rate of BGsentinel traps for Aedes aegypti (Diptera: Culicidae) in suburban Cairns, Australia. J Med Entomol. 2012;49:29–34.
 18.
Gubler DJ, Ooi EE, Vasudevan S, Farrar J. Dengue and dengue hemorrhagic fever. 2nd ed. Wallingford: CABI; 2014.
 19.
Charlwood JD, Kessy E, Yohannes K, Protopopoff N, Rowland M, LeClair C. Studies on the resting behaviour and host choice of Anopheles gambiae and An. arabiensis from Muleba, Tanzania. Med Vet Entomol. 2018;32:263–70.
 20.
Maia MF, Robinson A, John A, Mgando J, Simfukwe E, Moore SJ. Comparison of the CDC Backpack aspirator and the Prokopack aspirator for sampling indoor and outdoorresting mosquitoes in southern Tanzania. Parasites Vectors. 2011;4:124.
 21.
Onyango SA, Kitron U, Mungai P, Muchiri EM, Kokwaro E, King CH, et al. Monitoring malaria vector control interventions: effectiveness of five different adult mosquito sampling methods. J Med Entomol. 2013;50:1140–51.
 22.
VazquezProkopec GM, Galvin WA, Kelly R, Kitron U. A new, costeffective, batterypowered aspirator for adult mosquito collections. J Med Entomol. 2009;46:1256–9.
 23.
Bord S, Bioche C, Druilhet P. A cautionary note on Bayesian estimation of population size by removal sampling with diffuse priors. Biom J. 2018;60:450–62.
 24.
Stewart DR, Butler MJ, Johnson LA, Cajero A, Young AN, Harris GM. Efficacy of depletion models for estimating abundance of endangered fishes in streams. Fish Res. 2019;209:208–17.
 25.
Mitro MG, Zale AV. Predicting fish abundance using singlepass removal sampling. Can J Fish Aquat Sci. 2000;57:951–61.
 26.
Reisen WK, Mahmood F. Relative abundance, removal sampling, and markreleaserecapture estimates of population size of Anopheles culicifacies and Anopheles stephensi at diurnal resting sites in rural Punjab Province, Pakistan. Mosq News. 1981;41:22–30.
 27.
Daniels TJ, Falco RC, Fish D. Estimating population size and drag sampling efficiency for the blacklegged tick (Acari: Ixodidae). J Med Entomol. 2000;37:357–63.
 28.
PaviaRuz N, Rojas DP, Salha V, Granja P, BalamMay A, Longini IM, et al. Seroprevalence of dengue antibodies in three urban settings in Yucatan, Mexico. Am J Trop Med Hyg. 2018;98:1202–8.
 29.
Bisanzio D, DzulManzanilla F, GomezDantes H, PaviaRuz N, Hladish TJ, Lenhart A, et al. Spatiotemporal coherence of dengue, chikungunya and Zika outbreaks in Merida, Mexico. PLoS Negl Trop Dis. 2018;12:e0006298.
 30.
ManriqueSaide P, Coleman P, McCall PJ, Lenhart A, VazquezProkopec G, Davies CR. Multiscale analysis of the associations among egg, larval and pupal surveys and the presence and abundance of adult female Aedes aegypti (Stegomyia aegypti) in the city of Merida, Mexico. Med Vet Entomol. 2014;28:264–72.
 31.
Detinova TS, Beklemishev WN, Bertram DS. World Health Organization Agegrouping methods in Diptera of medical importance, vol. 47., World Health Organization. Monograph SeriesGeneva: World Health Organization; 1962. p. 1–213.
 32.
Carle F, Strub M. A new method for estimating population size from removal data. Biometrics. 1978;34:621–30.
 33.
Ogle D, Wheeler P, Dinno A. FSA: Fisheries stock analysis. R package version 0822; 2018. https://github.com/droglenc/FSA.
 34.
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixedeffects models using lme4. J Stat Softw. 2015;67:48.
 35.
Williams CR, Long SA, Russell RC, Ritchie SA. Field efficacy of the BGSentinel compared with CDC Backpack Aspirators and CO2baited EVS traps for collection of adult Aedes aegypti in Cairns, Queensland, Australia. J Am Mosq Control Assoc. 2006;22:296–300.
 36.
VazquezProkopec GM, MedinaBarreiro A, CheMendoza A, DzulManzanilla F, CorreaMorales F, GuillermoMay G, et al. Deltamethrin resistance in Aedes aegypti results in treatment failure in Merida, Mexico. PLoS Negl Trop Dis. 2017;11:e0005656.
 37.
Ritchie SA, Johnson BJ. Advances in vector control science: rearandrelease strategies show promise… but donʼt forget the basics. J Infect Dis. 2017;215(Suppl. 2):S103–8.
 38.
Andraud M, Hens N, Marais C, Beutels P. Dynamic epidemiological models for dengue transmission: a systematic review of structural approaches. PLoS ONE. 2012;7:e49085.
 39.
Wilson AL, Boelaert M, Kleinschmidt I, Pinder M, Scott TW, Tusting LS, et al. Evidencebased vector control? Improving the quality of vector control trials. Trends Parasitol. 2015;31:380–90.
 40.
Reiner RC Jr, Achee N, Barrera R, Burkot TR, Chadee DD, Devine GJ, et al. Quantifying the epidemiological impact of vector control on dengue. PLoS Negl Trop Dis. 2016;10:e0004588.
 41.
Bowman LR, RungeRanzinger S, McCall PJ. Assessing the relationship between vector indices and dengue transmission: a systematic review of the evidence. PLoS Negl Trop Dis. 2014;8:e2848.
 42.
Anders KL, le Nga H, Thuy NT, Ngoc TV, Tam CT, Tai LT, et al. Households as foci for dengue transmission in highly urban Vietnam. PLoS Negl Trop Dis. 2015;9:e0003528.
 43.
Stoddard ST, Forshey BM, Morrison AC, PazSoldan VA, VazquezProkopec GM, Astete H, et al. Housetohouse human movement drives dengue virus transmission. Proc Natl Acad Sci USA. 2013;110:994–9.
Acknowledgements
We would like to thank the residents of Merida who graciously agreed to be part of this study.
Funding
Research funding was provided by the Bureau for Global Health, U.S. Agency for International Development (USAID:AID‐OAA‐F‐16‐00094; PI, Devine G.), and by an Interagency Agreement between USAID and the US Centers for Disease Control and Prevention (CDC: OADS BAA 2016N17844; PI, VazquezProkopec G.M.), by the Canadian Institutes of Health Research (CIHR) and IDRC (Preventing Zika disease with novel vector control approaches, Project 108412) and by Fondo Mixto CONACyT (Mexico)Gobierno del Estado de Yucatan (Project YUC20170301556). The opinions expressed herein are those of the authors and do not necessarily reflect the views of the U.S. Agency for International Development. The findings and conclusions in this paper are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention.
Author information
Affiliations
Contributions
Conceptualized the research: EKC, GVP, AL and PMS. Conducted the research: EKC, AMB, ACR, NPR, GAT and MD. Wrote the manuscript: EKC, AL, NPR, MD, PMS and GVP. Contributed with reagents/data: NPR. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Protocols for this study were approved by Emory University’s ethics committee under protocol ID: IRB00082848. The protocol was also approved by the Ethics and Research Committee from the OʼHoran General Hospital from the state Ministry of Health, Register no. CEI034114. Involvement of CDC investigators was approved by the Office of the Associate Director for Science, Center for Global Health at CDC. Written informed consent was obtained from the head of household prior to mosquito collection.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
About this article
Cite this article
KoyocCardeña, E., MedinaBarreiro, A., CohuoRodríguez, A. et al. Estimating absolute indoor density of Aedes aegypti using removal sampling. Parasites Vectors 12, 250 (2019). https://doi.org/10.1186/s130710193503y
Received:
Accepted:
Published:
Keywords
 Sampling
 Entomology
 Population abundance