- Research article
- Open access
- Published:
Associations between household environmental factors and immature mosquito abundance in Quetzaltenango, Guatemala
BMC Public Health volume 19, Article number: 1729 (2019)
Abstract
Background
Aedes aegypti-borne diseases are becoming major public health problems in tropical and sub-tropical regions. While socioeconomic status has been associated with larval mosquito abundance, the drivers or possible factors mediating this association, such as environmental factors, are yet to be identified. We examined possible associations between proximity to houses and roads and immature mosquito abundance, and assessed whether these factors and mosquito prevention measures mediated any association between household environmental factors and immature mosquito abundance.
Methods
We conducted two cross-sectional household container surveys in February–March and November–December, 2017, in urban and rural areas of Quetzaltenango, Guatemala. We used principal components analysis to identify factors from 12 variables to represent the household environment. One factor which included number of rooms in house, electricity, running water, garbage service, cable, television, telephone, latrine, well, and sewer system, was termed “environmental capital.” Environmental capital scores ranged from 0 to 5.5. Risk factors analyzed included environmental capital, and distance from nearest house/structure, paved road, and highway. We used Poisson regression to determine associations between distance to nearest house/structure, roads, and highways, and measures of immature mosquito abundance (total larvae, total pupae, and positive containers). Using cubic spline generalized additive models, we assessed non-linear associations between environmental capital and immature mosquito abundance. We then examined whether fumigation, cleaning containers, and distance from the nearest house, road, and highway mediated the relationship between environmental capital and larvae and pupae abundance.
Results
We completed 508 household surveys in February–March, and we revisited 469 households in November–December. Proximity to paved roads and other houses/structures was positively associated with larvae and pupae abundance and mediated the associations between environmental capital and total numbers of larvae/pupae (p ≤ 0.01). Distance to highways was not associated with larval/pupal abundance (p ≥ 0.48). Households with the lowest and highest environmental capital had fewer larvae/pupae than households in the middle range (p < 0.01).
Conclusions
We found evidence that proximity to other houses and paved roads was associated with greater abundance of larvae and pupae. Understanding risk factors such as these can allow for improved targeting of surveillance and vector control measures in areas considered at higher risk for arbovirus transmission.
Background
Approximately 6.01 billion people currently live in areas suitable for Aedes aegypti disease transmission [1]. Ae. aegypti-borne diseases, such as dengue (DENV), chikungunya (CHIKV), and Zika (ZIKV) viruses, are found in tropical and subtropical zones with an abundance of these species, including Central America [2,3,4]. Other than for yellow fever vaccine [5], no broadly licensed commercial vaccines are available for the principal Ae. aegypti-borne arboviruses, so vector control remains the primary strategy to limit their spread [6]. Climate change, urbanization, migration, human behavior, and ecosystem modification are among the myriad factors influencing the geographic spread of Ae. aegypti and their associated viruses [1, 7, 8].
Ae. aegypti are highly productive in urban environments and have a strong preference for human blood [9]. Ae. aegypti spend the majority of their lives in the houses where they emerged, flying an average of 40–80 m during the course of their lifetimes [10]. Oviposition sites are selected based on their physical, chemical, and biological characteristics, such as container type, depth, water quality, and sun exposure [11, 12]. Ideal larval habitats for Ae. aegypti are dark-colored containers filled with stagnant water and organic material in shaded areas around houses [11, 13, 14]. Productive container types include flower pots, tires, vases, buckets, cans, rain gutters, fountains, bottles, and birdbaths [11, 13, 14]. Greater human population densities provide more feeding opportunities for Ae. aegypti [15].
Studies of socioeconomic status (SES) impacts on Ae. aegypti abundance mostly report greater Ae. aegypti population densities in low SES areas [16,17,18,19,20,21,22]. Most studies have only considered income, occupation, and education as the SES factors. Few studies have evaluated associations between household environmental measures as attributes of SES and mosquito abundance. The household environmental factors that can influence mosquito infestation are quite heterogeneous. These include piles of garbage [21], open wells [23, 24], storm sewers [25], and septic tanks [26]. Less information is available on spatial risk factors, but proximity to vacant lots [27, 28], vegetation or green spaces [29], other houses/structures [30], and roads [31, 32], have been shown to be predictive of mosquito abundance. Household infrastructure may also influence the mosquito microenvironment [33,34,35]. For example, the premise condition index has been shown to be an effective tool at classifying houses according to risk of having mosquito breeding sites [33,34,35]. This index can be used to prioritize neighborhoods for vector control interventions.
For this study, we evaluated whether proximity to other houses/structures and roads, and household environmental factors were associated with immature mosquito abundance. A secondary objective was to determine how mosquito abatement interventions, including fumigation and cleaning possible larval habitat containers, influence immature mosquito abundance. It is particularly important to examine these relationships in Central America, which has been host to large outbreaks of arbovirus infection and where vector control resources are limited [36].
Methods
Study site
We selected two municipalities in the Guatemalan department of Quetzaltenango, Coatepeque and Génova (Fig. 1), as study sites based on their high risk for arboviral disease transmission and high mosquito pupal index (> 25% of houses with pupal infestations) [37]. Coatepeque (14°42′00″N 91°52′00″O) and Génova (14°37′00″N 91°50′00″O) are located in the south-western region of the Republic of Guatemala and have a tropical climate. The mean annual temperatures for Coatepeque and Génova are 25.7 °C and 26.2 °C, respectively, the mean annual precipitations are 308 mm and 285 mm, and the mean elevations are 498 m and 350 m [38]. This study included two communities in Coatepeque (La Unión, El Jardín) and six communities in Génova (30 de Junio, Robles, Nueva Italia, Génova, San Jose, Guadalupe) (Additional file 1: Figure S1). The communities in Coatepeque were selected based on the presence of Vigilancia Integrada Comunitaria (Integrated Community Surveillance), a prospective public health syndromic surveillance system for diarrheal, respiratory, and febrile illnesses of the Centro de Estudios en Salud/Universidad del Valle de Guatemala in collaboration with the Guatemalan Ministry of Health and the United States Centers for Disease Control and Prevention (CDC). In Génova, all of the communities reporting a high pupal index were included, with the exception of one community that posed a security risk for field personnel. Six sites in Génova were selected to achieve comparable population size to the two sites in Coatepeque. We remotely identified each probable house structure within each community using Google satellite imagery for 2016 in QGIS 2.2 (QGIS Development Team, 2019). The Ministerio de Salud Pública y Asistencia Social (MSPAS) provided detailed maps of each community in order to demonstrate community boundaries. All probable houses were identified and verified on-site to confirm classification of structures [39]. Houses were then randomly selected in each village using a two-stage sampling procedure based on a geographic 100 × 100 m grid. We first randomly selected grids, enumerated households, and then used a random number generator to select one house within each grid. In both Coatepeque and Génova, selected houses accounted for 10% of the total community population (n = 250 and n = 258, respectively). If no one was at home during recruitment, if householders chose not to participate, or if the selected structure was not a house, we selected the nearest house to the right of the front door as a replacement.
Container inspection and questionnaire
After obtaining informed consent from homeowners, we conducted cross-sectional surveys for container-inhabiting mosquitoes in February–March, 2017 (the local dry season) and November–December, 2017 (the local rainy season) in both Génova and Coatepeque. We conducted two surveys to capture immature mosquito abundance in Guatemala’s two seasons. All containers ≥3 L inside and outside the houses were inspected for any genera of mosquito larvae and pupae, and total numbers of mosquito larvae and pupae from all containers in each house and the containers with any mosquito larvae or pupae were recorded. Larvae and pupae were analyzed separately, as pupal counts are considered more representative of local adult mosquito populations [40, 41]. We did not identify larval and pupal genus or species. We interviewed the heads of household or another adult residing in the house, and responses were transcribed onto Excel spreadsheets. Questions covered mosquito control measures, waste disposal, and socioeconomic indicators.
Variables
We assessed household environment factors and distance from nearest house/structure, paved road, and main transportation corridor running through the city/village as risk factors for vector concentrations. We assumed the main transportation corridor was the nearest highway or the only paved road in villages that did not have highway access.
We used principal components factor analysis to identify factors based on 12 variables from the first household survey to represent household attributes of SES. These included: number of rooms in the house (1–4, > 5), electricity (yes, no), running water (yes, no), a television (yes, no), a landline telephone (yes, no), a latrine (yes, no), cable television service (yes, no), a mobile phone (yes, no), trash disposal service (yes, no), a water well (yes, no), sewer system (yes, no), and a rainwater collection system (yes, no). The resultant compound factor, which we termed “environmental capital,” included all of the variables except a mobile phone and rainwater collection system (Additional file 1: Table S1). Variables highly correlated with the factor were weighted against their eigenvector. This factor reflects some of the attributes of the Encuesta Nacional de Salud Materno Infantil (National Survey of Maternal and Child Health), which focuses on the health of children and adults in Guatemala [42]. This household factor from the first survey explained 32% of the variability in the data and was used to represent environmental capital in the second survey as well. Higher environmental capital scores indicated higher SES and ranged from 0 to 5.5.
The measures of immature mosquito abundance were the total number larvae (continuous), total number of pupae (continuous), and positive containers (continuous). Categorical covariates included survey period (February–March vs. November–December), residence (urban vs. rural), self-reported cleaned (scrubbed, treated, or emptied standing water) containers (barrels, pots, tires, etc.) at least once in the last 6 months (yes, no), and self-reported homeowner or vector control authority fumigation inside/outside house at least once in the last 6 months (yes, no). Continuous covariates included the number of people in a household and the total number of containers ≥3 L with water at the time of the visit per household (e.g., buckets, barrels, flower pots, etc.). ‘Urban’ residences were those in El Jardín, Coatepeque, whereas ‘rural’ residences were all other communities, as defined by the census [43].
Spatial analysis
Coordinates of each house were entered into geographical information system software (ArcGIS Pro 2.2.4 software; ESRI, Redlands, CA) and overlaid on basemaps and satellite images from December 8, 2018, of Coatepeque and Génova [44]. These maps were used to locate and visualize households and roads. We collected ground truth data through site visits during both survey periods. The distance between a house and its closest neighboring house or other structure (e.g., store, church) or road was ascertained by measuring the Euclidean distance between points taken from the front door of the house to the closest edge of lines representing roads [45,46,47]. Within the sub-set of sampled houses in each community, we also attempted to detect spatial clusters of houses with larval infestations.
Statistical analysis
Medians and interquartile ranges were reported for continuous variables (total number of larvae; total number of pupae; number of positive containers; number of containers ≥3 l; number of people in household; distance to nearest paved road, highway, and house/structure; environmental capital). Frequency distributions were reported for categorical variables (cleaned containers, fumigation, urban/rural residence).
We used Poisson regression, which is used to model count data, to analyze unadjusted (Model 1) and adjusted (Model 2) associations between hypothesized risk factors (distance to nearest house/structure, paved road, highway), and immature mosquito abundance (number of larvae, pupae, and positive containers), with household as a repeated measure (two time points). We used generalized estimating equations to estimate the population-averaged effect and used compound symmetry as the covariance structure to account for correlations resulting from two measurements (February–March, November–December) of immature mosquito abundance on the same houses within each site. In Model 2, we used directed acyclic graphs [48, 49] to select each covariate for model inclusion based on a priori importance and evidence from the scientific literature of being potential confounders of associations between our exposures of interest and mosquito larvae and pupae abundance (Additional file 1: Figure S2). The adjusted models included environmental capital (categorized by tertiles) [50, 51], survey period [52], urban/rural residence [53], the number of people per household [54], cleaned containers [55], fumigated inside/outside house [56], and the total number of containers ≥3 l per household [50, 57]. Tolerance values were used to assess potential collinearity between all independent variables [58]. Due to the potential over-dispersion of larvae and pupae abundance, negative binomial regression models were fitted to evaluate the same associations as a sensitivity analysis [59].
We then used cubic spline generalized additive models to explore potential non-linear relationships between environmental capital and immature mosquito abundance (number of larvae, pupae, and positive containers) separately for both survey periods.
Finally, we assessed whether factors including fumigation, cleaned containers, and distance to nearest paved road, highway, and household/structure mediated the relationship between environmental capital and the total number of larvae, pupae, and positive containers. This analysis followed causal mediation analysis methods as previously described by VanderWeele [60]. The mediation models were Poisson models to estimate the association between environmental capital and the distance to the nearest house/structure, paved road, and highway, and binomial models to estimate the association between environmental capital and cleaned containers and fumigation history, which are dichotomous variables. The outcome models were Poisson models that estimated the association between environmental capital and immature mosquito indicators (number of larvae, pupae, and positive containers), adjusting for the mediators. All hypothesized mediators were included in outcome models. The “mediation” package in R 3.5.2 statistical software (R Development Core Team, Vienna, Austria) was used for multilevel causal mediation analyses [61]. We ran one thousand Monte Carlo simulations in this analysis for variance estimation. Estimates, standard errors, and the proportion mediated were reported. All analyses, other than mediation, were calculated using SAS V.9.4 (SAS Institute, Inc., Cary, North Carolina).
Results
Household characteristics
In February–March, 508 household inspections were completed. In November–December, 469 of those households (92.3%) were revisited for a second survey (some houses were not revisited because the homeowner was unavailable). An additional 18 households that were eligible but unavailable during the first survey were included in the second survey. Of all houses, 72.7% were in rural areas (Table 1). There was a median of five people per household. The median distances to the nearest house/structure, paved road, and highway were 3.1 m, 13.9 m and 244.1 m for rural residences and 1 m, 4.9 m, and 144.3 m for urban residences, respectively. The median numbers of larvae, pupae, and positive containers were 8, 1, and 1 in rural residences and 20, 2, and 1 in urban residences, respectively.
Geographical distances
Distance to the nearest paved road was inversely associated with the total number of larvae, pupae, and positive containers per house in Models 1 and 2 (p ≤ 0.01) (Table 2). For every 10-m increase in distance from the nearest paved road, the total number of larvae and positive containers decreased by a factor of 0.96 and the number of pupae decreased by a factor of 0.93, adjusting for environmental capital, urban/rural residence, the number of people per household, cleaned containers, fumigation history, and the total number of containers. Tolerance values were above 0.50, so there was no evidence of collinearity among any of the independent variables.
Distance to the nearest highway was not associated with the number of larvae, pupae, or positive containers per household in Models 1 and 2 (p ≥ 0.28) (Table 2).
Distance from the nearest household/structure was inversely associated with the total number of larvae and pupae and number of positive containers per house in Models 1 and 2 (p < 0.01) (Table 2). For every 1-m increase in distance from the nearest house/structure, the total number of larvae and positive containers decreased by a factor of 0.97 and the number of pupae decreased by a factor of 0.95, adjusting for relevant covariates. Full model outputs are reported in Additional file 1: Tables S2–S4, but these estimates should be interpreted with caution, because the relationships between the covariates and outcomes are not adjusted for confounders [49]. Results from negative binomial models were similar for distance to the nearest paved road, highway, and house/structure (Additional file 1: Table S5).
We did not verify measurements obtained using ArcGIS between houses and roads on the ground, but the ground resolution of the ArcGIS world imagery for our study sites is 0.46 m and objects in the map are within 5 m of their true location [62].
Spatial clusters of larvae and pupae
High/Low clustering (Getis-Ord General G) analyses did not reveal spatially dependent clusters for immature mosquito abundance indicators (number of larvae, pupae, and positive containers) for either time point (p ≥ 0.40).
Environmental capital
Cubic splines demonstrated significant non-linear relationships between environmental capital and the number of larvae and pupae per house that were similar for both survey periods (p < 0.01) (Fig. 2). For both surveys, households with the lowest and highest environmental capital had significantly fewer larvae and pupae compared to those in the middle (p < 0.01). Results for the number of positive containers were similar (Additional file 1: Fig. S3).
Distance to the nearest paved road and house/structure were significant mediators of the relationship between environmental capital and the number of larvae and pupae (p < 0.01) (Table 3). A one-unit increase in environmental capital was associated with a significant decrease in distance from the nearest paved road or house/structure, which in turn was associated with more larvae and pupae when environmental capital was held constant (p < 0.01). Fumigated houses, cleaned containers, and distance to the nearest highway were not significant mediators of the association between environmental capital and the number of larvae and pupae. Results for the number of positive containers were similar (Additional file 1: Table S6).
Discussion
This study identified environmental factors and SES attributes that were associated with mosquito larvae and pupae abundance. Distance to the nearest paved road and house/structure were inversely associated with larvae and pupae abundance and were significant mediators of the relationship between environmental capital and the number of larvae and pupae per house. Cubic splines revealed that households of middle environmental capital had significantly more larvae and pupae than those with the lowest and highest environmental capital.
Our finding that households closer to paved roads had more larvae and pupae is consistent with previous studies from Kansas and Bermuda, which found greater numbers of adult mosquitoes and eggs closer to roads [31, 32]. Proximity to paved roads may indicate greater population density, which would include more containers and greater availability of blood meals. The association remained significant after adjusting for the total number of containers ≥3 L per household, which may suggest a greater presence of smaller containers like cups, cans, and bottles, in areas closer to roads [31]. These containers are also conceivably productive larval habitats. This association was further supported by mediation analyses, which showed that distance to the nearest paved road was a significant mediator of the relationship between environmental capital and number of larvae and pupae. As environmental capital increased, distance to the nearest paved road decreased. Households closer to paved roads had significantly more larvae and pupae, holding environmental capital constant. It is conceivable that households with greater environmental capital, which are closer to roads, are more likely to own barrels and other large water storage containers, which may support larger mosquito populations if they are not properly managed. More mosquitoes in areas closer to paved roads may also increase the risk of the spread of arboviral infections, which was reported in a CHIKV study in Pakistan [63].
Distance to the nearest highway was not a significant predictor of larvae and pupae abundance. One study in Taiwan reported that the number of dengue fever cases corresponded inversely with distance from highways, further indicating that Ae. aegypti abundance may be associated with population density [64]. Proximity to highways in our study was not necessarily suggestive of greater human population density, which may have greater influence on mosquito abundance [65, 66]. These results may suggest that the immediate household environment contributes more to larvae and pupae abundance than more distant neighborhood factors [67,68,69]. This is particularly important for Ae. aegypti, as immatures tend to be highly aggregated in space and time, rarely dispersing beyond 30–40 m of the household where they developed as larvae [67, 69].
Distance to the nearest house/structure was inversely associated with larvae and pupae abundance. Furthermore, mediation analyses revealed that households with higher environmental capital were closer to other houses/structures and had significantly more larvae and pupae. We are unaware of other studies assessing distance to the nearest structure as a mediator between SES and mosquito abundance. Previous studies of associations between distance to the nearest building and mosquito abundance are inconsistent. Some report greater Anopheles and Aedes abundance in houses/structures closer together [30, 70, 71], whereas others do not [31, 72]. Urbanization and greater human population density lead to a greater number of artificial containers, which creates an abundance of potential habitats for mosquitoes, including tires, flowerpots, and cans [15]. Urban environments may also be more favorable for Ae. aegypti due to the absence of natural vegetation, competition, and predation [12, 15, 73, 74]. These results reinforce the premise that mosquito control requires community-wide efforts, as individual houses with disproportionately high numbers of mosquitoes may pose risks to their closest neighbors, and indeed the entire community [68].
Recent history of fumigation inside/outside of the house and containers that had been cleaned but could still serve as immature habitats for mosquitoes were not significant mediators between environmental capital and the number of larvae and pupae. Fumigating and cleaning containers with standing water are established mosquito control measures [55, 56, 75]. Fumigation is only provided by MSPAS in Guatemala. It could be that our measure of environmental capital was not predictive of these preventive measures in these communities or that fumigation may not have been effective in these areas. Alternatively, our cross-sectional survey that asked whether participants performed these prevention measures in the last 6 months may have been insufficient to assess the efficacy of these interventions, which require repeated application. Fumigation frequency and insecticide resistance should be also considered.
Households of middle environmental capital had significantly more larvae and pupae than households with the lowest and highest environmental capital for both surveys. In this study, environmental capital included access to running water, improved sanitation, a sewer system, and trash disposal service, which are typically associated with reduced mosquito populations [23, 24, 52, 76,77,78]. Greater environmental capital may also indicate higher values of other SES indicators, including income, occupation, and education, which are associated with greater mosquito prevention measures, such as removing containers with standing water [17,18,19, 21]. Conversely, low environmental capital was associated with greater distance to the nearest paved road, which was associated with fewer mosquitoes. It is conceivable that these distances exceeded the typical flight range for mosquitoes [79]. Moreover, houses with low environmental capital in this study had fewer barrels and other large water storage containers that were most productive for mosquitoes.
Our study did not characterize larval genus or species, but multiple species of Aedes, Anopheles, and Culex mosquitoes have been reported in Quetzaltenango department, where our study was conducted [80,81,82,83]. Specific species in Quetzaltenango include Ae. aegypti and Ae. albopictus [80, 81], which preferentially lay eggs in household containers [84]; An. hectoris, An. parapunctipennis, and An. xelajuensis, which prefer marshes, trees, swamps, fields, streams, and rivers [85]; and Cx. corniger, Cx. peus, and Cx. quinquefasciatus, whose breeding sites include storm sewers, cesspits, and polluted water [26, 86]. Given our container surveys occurred exclusively in households, we suspect that the majority of the immatures that we collected were either Ae. aegypti or Ae. albopictus.
Our study had several limitations. First, we sampled communities based on high entomological indices and are thus these are not representative of all communities in Guatemala. However, the households are representative of the local communities. Second, cross-sectional surveys of mosquitoes are time sensitive [41] and our two surveys points were insufficient to fully capture the temporal variability of mosquito larvae and pupae, despite including both dry and rainy seasons. Third, our survey assessments of whether participants fumigated inside/outside the house or cleaned their containers in the last 6 months were likely inadequate to assess the efficacy of these prevention strategies. Fourth, we did not include containers < 3 L on household premises such as discarded cups and cans, which could also serve as immature mosquito habitats.
Conclusions
The global human population is expected to peak around 9.6 billion by 2050, favoring the spread of vector-borne diseases [87, 88]. With climate change, increasing temperatures, and more frequent flooding, the geographic range of Ae. aegypti and Ae. albopictus is increasing [1, 89]. The findings reported here provide evidence that proximity to other houses/structures and paved roads was associated with more mosquito larvae and pupae in containers around households. Furthermore, households with higher environmental capital were closer to other houses/structures and paved roads, and had significantly greater larvae and pupae abundance. Finally, households with middle environmental capital had significantly more larvae and pupae than the lowest and highest tiers. In resource-limited vector control programs, findings such as these can be used to focus efforts on areas of greater population density closer to roads. The findings also highlight the importance of programs that take into account neighborhood-level risks and mitigation strategies when promoting the prevention of vector-borne diseases.
Availability of data and materials
The surveys and data that support the findings of this study are not in English, but are available from the Centro de Estudios en Salud, Universidad del Valle de Guatemala, with permissions from the Centers for Disease Control and Prevention.
Abbreviations
- Ae. :
-
Aedes
- An. :
-
Anopheles
- CDC:
-
Centers for Disease Control and Prevention
- CHIKV:
-
chikungunya virus
- Cx. :
-
Culex
- DENV:
-
dengue virus
- IQR:
-
interquartile range
- MSPAS:
-
Ministerio de Salud Pública y Asistencia Social
- SE:
-
standard error
- SES:
-
socioeconomic status
- ZIKV:
-
Zika virus
References
Ryan SJ, Carlson CJ, Mordecai EA, Johnson LR. Global expansion and redistribution of Aedes-borne virus transmission risk with climate change. PLoS Negl Trop Dis. 2019;13(3):e0007213.
Pan American Health Organization. Reported cases of dengue fever in the Americas: PLISA: health information platform for the Americas; 2018 [cited 2018 Oct 3]. Available from: http://www.paho.org/data/index.php/en/mnu-topics/indicadores-dengue-en/dengue-nacional-en/252-dengue-pais-ano-en.html.
Pan American Health Organization. Zika suspected and confirmed cases reported by countries and territories in the Americas: cumulative cases, 2015–2018. 2018 [cited 2018 Sept 16]. Available from: https://www.paho.org/hq/index.php?option=com_content&view=article&id=12390:zika-cumulative-cases&Itemid=42090&lang=es.
Pan American Health Organization. Chikungunya: data, maps and statistics 2018 [cited 2018 Sept 16]. Available from: https://www.paho.org/hq/index.php?option=com_topics&view=rdmore&cid=5927&item=chikungunya&type=statistics&Itemid=40931&lang=en.
World Health Organization. Yellow fever: key facts 2018 [cited 2018 Sept 16]. Available from: https://www.who.int/news-room/fact-sheets/detail/yellow-fever.
Ritchie SA, van den Hurk AF, Smout MJ, Staunton KM, Hoffmann AA. Mission accomplished? We need a guide to the ‘Post Release’World of Wolbachia for Aedes-borne disease control. Trends Parasitol. 2018;34:217–26.
Kilpatrick AM, Randolph SE. Drivers, dynamics, and control of emerging vector-borne zoonotic diseases. Lancet. 2012;380(9857):1946–55.
Paixão ES, Teixeira MG, Rodrigues LC. Zika, chikungunya and dengue: the causes and threats of new and re-emerging arboviral diseases. BMJ Glob Health. 2018;3(Suppl 1):e000530.
World Health Organization. Comprehensive guideline for prevention and control of dengue and dengue haemorrhagic fever 2011 [cited 2019 July 29]. Available from: https://apps.who.int/iris/handle/10665/204894.
Teurlai M, Menkes CE, Cavarero V, Degallier N, Descloux E, Grangeon J-P, et al. Socio-economic and climate factors associated with dengue fever spatial heterogeneity: a worked example in New Caledonia. PLoS Negl Trop Dis. 2015;9(12):e0004211.
Getachew D, Tekie H, Gebre-Michael T, Balkew M, Mesfin A. Breeding sites of Aedes aegypti: potential dengue vectors in Dire Dawa, East Ethiopia. Interdiscip Perspect Infect Dis. 2015;2015:706276.
Zahouli JB, Koudou BG, Müller P, Malone D, Tano Y, Utzinger J. Urbanization is a main driver for the larval ecology of Aedes mosquitoes in arbovirus-endemic settings in South-Eastern Côte d'Ivoire. PLoS Negl Trop Dis. 2017;11(7):e0005751.
Centers for Disease Control and Prevention. Dengue and the Aedes aegypti mosquito 2012 [cited 2018 Oct 4]. Available from: https://www.cdc.gov/dengue/resources/30Jan2012/aegyptifactsheet.pdf.
Ferdousi F, Yoshimatsu S, Ma E, Sohel N, Wagatsuma Y. Identification of essential containers for Aedes larval breeding to control dengue in Dhaka, Bangladesh. Trop Med Health. 2015;43:253–64.
Li Y, Kamara F, Zhou G, Puthiyakunnon S, Li C, Liu Y, et al. Urbanization increases Aedes albopictus larval habitats and accelerates mosquito development and survivorship. PLoS Negl Trop Dis. 2014;8(11):e3301.
Chambers D, Young L, Hill JH. Backyard mosquito larval habitat availability and use as influenced by census tract determined resident income levels. J Am Mosq Control Assoc. 1986;2(4):539–44.
Quintero J, Carrasquilla G, Suárez R, González C, Olano VA. An ecosystemic approach to evaluating ecological, socioeconomic and group dynamics affecting the prevalence of Aedes aegypti in two Colombian towns. Cad Saude Publica. 2009;25:s93–s103.
Unlu I, Farajollahi A, Healy SP, Crepeau T, Bartlett-Healy K, Williges E, et al. Area-wide management of Aedes albopictus: choice of study sites based on geospatial characteristics, socioeconomic factors and mosquito populations. Pest Manag Sci. 2011;67(8):965–74.
Spiegel JM, Bonet M, Ibarra AM, Pagliccia N, Ouellette V, Yassi A. Social and environmental determinants of Aedes aegypti infestation in Central Havana: results of a case–control study nested in an integrated dengue surveillance programme in Cuba. Tropical Med Int Health. 2007;12(4):503–10.
Dhar-Chowdhury P, Haque CE, Lindsay R, Hossain S. Socioeconomic and ecological factors influencing Aedes aegypti prevalence, abundance, and distribution in Dhaka, Bangladesh. Am J Trop Med Hyg. 2016;94(6):1223–33.
Syed M, Saleem T, Syeda U-R, Habib M, Zahid R, Bashir A, et al. Knowledge, attitudes and practices regarding dengue fever among adults of high and low socioeconomic groups. J Pak Med Assoc. 2010;60(3):243.
Caprara A. Lima JWdO, Marinho ACP, Calvasina PG, Landim LP, Sommerfeld J. irregular water supply, household usage and dengue: a bio-social study in the Brazilian northeast. Cad Saude Publica. 2009;25:S125–S36.
Yamamoto S, Louis V, Sié A, Sauerborn R. Household risk factors for clinical malaria in a semi-urban area of Burkina Faso: a case–control study. Trans R Soc Trop Med Hyg. 2010;104(1):61–5.
Pålsson K, Jaenson TG, Dias F, Laugen AT, Björkman A. Endophilic Anopheles mosquitoes in Guinea Bissau, West Africa, in relation to human housing conditions. J Med Entomol. 2004;41(4):746–52.
Manrique-Saide P, Uc V, Prado C, Carmona C, Vadillo J, Chan R, et al. Storm sewers as larval habitats for Aedes aegypti and Culex spp. in a neighborhood of Merida, Mexico. J Am Mosq Control Assoc. 2012;28(3):255–7.
Burke R, Barrera R, Lewis M, Kluchinsky T, Claborn D. Septic tanks as larval habitats for the mosquitoes Aedes aegypti and Culex quinquefasciatus in playa-Playita, Puerto Rico. Med Vet Entomol. 2010;24(2):117–23.
Baak-Baak CM, Arana-Guardia R, Cigarroa-Toledo N, Loroño-Pino MA, Reyes-Solis G, Machain-Williams C, et al. Vacant lots: productive sites for Aedes (Stegomyia) aegypti (Diptera: Culicidae) in Merida City, Mexico. J Med Entomol. 2014;51(2):475–83.
Costa F, Fattore G, Abril M. Diversity of containers and buildings infested with Aedes aegypti in Puerto Iguazú, Argentina. Cad Saude Publica. 2012;28:1802–6.
Estallo E, Sangermano F, Grech M, Ludueña-Almeida F, Frías-Cespedes M, Ainete M, et al. Modelling the distribution of the vector Aedes aegypti in a central Argentine city. Med Vet Entomol. 2018;32:451–61.
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.
Kaplan L, Kendell D, Robertson D, Livdahl T, Khatchikian C. Aedes aegypti and Aedes albopictus in Bermuda: extinction, invasion, invasion and extinction. Biol Invasions. 2010;12(9):3277–88.
Ganser C, Wisely SM. Patterns of spatio-temporal distribution, abundance, and diversity in a mosquito community from the eastern Smoky Hills of Kansas. J Vector Ecol. 2013;38(2):229–36.
Andrighetti MTM, Galvani KC, da Graça Macoris ML. Evaluation of premise condition index in the context of Aedes aegypti control in Marília, São Paulo, Brazil. Dengue Bull. 2009;33:167–75.
Peres RC, Rego R, Maciel-de-Freitas R. The use of the premise condition index (PCI) to provide guidelines for Aedes aegypti surveys. J Vector Ecol. 2013;38(1):190–3.
Tun-Lin W, Kay B, Barnes A. The premise condition index: a tool for streamlining surveys of Aedes aegypti. Am J Trop Med Hyg. 1995;53(6):591–4.
Benelli G, Mehlhorn H. Declining malaria, rising of dengue and Zika virus: insights for mosquito vector control. Parasitol Res. 2016;115(5):1747–54.
Ministerio de Salud Pública y Asistencia Social. Sistema de Información Gerencial de Salud 2019 [cited 2019 March 8]. Available from: https://www.mspas.gob.gt/.
Anacafe. Estaciones meteorológicas - clima 2018 [cited 2019 March 8]. Available from: http://meteorologia.anacafe.org/Clima/.
Paquet C, Daniel M, Kestens Y, Léger K, Gauvin L. Field validation of listings of food stores and commercial physical activity establishments from secondary data. Int J Behav Nutr Phys Act. 2008;5(1):58.
Focks DA, Chadee DD. Pupal survey: an epidemiologically significant surveillance method for Aedes aegypti: an example using data from Trinidad. Am J Trop Med Hyg. 1997;56(2):159–67.
Morrison AC, Gray K, Getis A, Astete H, Sihuincha M, Focks D, et al. Temporal and geographic patterns of Aedes aegypti (Diptera: Culicidae) production in Iquitos, Peru. J Med Entomol. 2004;41(6):1123–42.
Ministerio de Salud Pública y Asistencia Social INdE, ICF International. Encuesta Nacional de Salud Materno Infantil 2014–2015: informe Final [Internet]. 2017.
Instituto Nacional de Estadística - INE/Guatemala. Catálogo de departamentos, municipios y lugares poblados de Guatemala (Censo 2002) 2002 [cited 2019 Feb 20]. Available from: https://www.ine.gob.gt/sistema/uploads/2014/02/20/jZqeGe1H9WdUDngYXkWt3GIhUUQCukcg.pdf.
ESRI. World topographic map 2018 [cited 2019 Feb 20]. Available from: https://www.arcgis.com/home/webmap/viewer.html.
Sander HA, Ghosh D, van Riper D, Manson SM. How do you measure distance in spatial models? An example using open-space valuation. Environ Plann B: Plann Des. 2010;37(5):874–94.
Mahabir R, Severson D, Chadee D. Impact of road networks on the distribution of dengue fever cases in Trinidad, West Indies. Acta Trop. 2012;123(3):178–83.
ESRI (Earth Science Resource Institute). How average nearest neighbor works 2019 [cited 2019 Nov 11]. Available from: https://pro.arcgis.com/en/pro-app/tool-reference/spatial-statistics/h-how-average-nearest-neighbor-distance-spatial-st.htm.
Weng H-Y, Hsueh Y-H, Messam LLM, Hertz-Picciotto I. Methods of covariate selection: directed acyclic graphs and the change-in-estimate procedure. Am J Epidemiol. 2009;169(10):1182–90.
Westreich D, Greenland S. The table 2 fallacy: presenting and interpreting confounder and modifier coefficients. Am J Epidemiol. 2013;177(4):292–8.
Velasco-Salas ZI, Sierra GM, Guzmán DM, Zambrano J, Vivas D, Comach G, et al. Dengue seroprevalence and risk factors for past and recent viral transmission in Venezuela: a comprehensive community-based study. Am J Trop Med Hyg. 2014;91(5):1039–48.
Sissoko D, Moendandze A, Malvy D, Giry C, Ezzedine K, Solet JL, et al. Seroprevalence and risk factors of chikungunya virus infection in Mayotte, Indian Ocean, 2005-2006: a population-based survey. PLoS One. 2008;3(8):e3066.
Fuller TL, Calvet G, Estevam CG, Angelo JR, Abiodun GJ, Halai U-A, et al. Behavioral, climatic, and environmental risk factors for Zika and Chikungunya virus infections in Rio de Janeiro, Brazil, 2015-16. PLoS One. 2017;12(11):e0188002.
Ali A, ur Rehman H, Nisar M, Rafique S, Ali S, Hussain A, et al. Seroepidemiology of dengue fever in Khyber Pakhtunkhawa, Pakistan. Int J Infect Dis. 2013;17(7):e518–e23.
Vincenti-Gonzalez MF, Grillet M-E, Velasco-Salas ZI, Lizarazo EF, Amarista MA, Sierra GM, et al. Spatial analysis of dengue seroprevalence and modeling of transmission risk factors in a dengue hyperendemic city of Venezuela. PLoS Negl Trop Dis. 2017;11(1):e0005317.
Phuanukoonnon S, Mueller I, Bryan JH. Effectiveness of dengue control practices in household water containers in Northeast Thailand. Tropical Med Int Health. 2005;10(8):755–63.
Kenneson A, Beltrán-Ayala E, Borbor-Cordova MJ, Polhemus ME, Ryan SJ, Endy TP, et al. Social-ecological factors and preventive actions decrease the risk of dengue infection at the household-level: results from a prospective dengue surveillance study in Machala, Ecuador. PLoS Negl Trop Dis. 2017;11(12):e0006150.
Brunkard JM, López JLR, Ramirez J, Cifuefntes E, Rothenberg SJ, Hunsperger EA, et al. Dengue fever seroprevalence and risk factors, Texas–Mexico border, 2004. Emerg Infect Dis. 2007;13(10):1477.
Midi H, Sarkar SK, Rana S. Collinearity diagnostics of binary logistic regression model. J Interdiscip Math. 2010;13(3):253–67.
Lawless JF. Negative binomial and mixed poisson regression. Can J Stat. 1987;15(3):209–25.
VanderWeele T. Explanation in causal inference: methods for mediation and interaction. Oxford: Oxford University Press; 2015.
Tingley D, Yamamoto T, Hirose K, Keele L, Imai K. Mediation: R package for causal mediation analysis. 2014.
ESRI. ArcGIS Online Standard Service: World Imagery Collection, Map Server. Maps throughout this book were created using ArcGIS® software by ESRI. ArcGIS® and ArcMapTM Are the Intellectual Property of ESRI and Are Used Herein under License, Copyright ©ESRI. 2018. 2018 [cited 2019 April 28]. Available from: https://www.arcgis.com/index.html.
Hira FS, Asad A, Farrah Z, Basit RS, Mehreen F, Muhammad K. Patterns of occurrence of dengue and chikungunya, and spatial distribution of mosquito vector Aedes albopictus in Swabi district, Pakistan. Tropical Med Int Health. 2018;23(9):1002–13.
Hsueh Y-H, Lee J, Beltz L. Spatio-temporal patterns of dengue fever cases in Kaoshiung City, Taiwan, 2003–2008. Appl Geogr. 2012;34:587–94.
Obenauer JF, Joyner TA, Harris JB. The importance of human population characteristics in modeling Aedes aegypti distributions and assessing risk of mosquito-borne infectious diseases. Trop Med Health. 2017;45(1):38.
Padmanabha H, Durham D, Correa F, Diuk-Wasser M, Galvani A. The interactive roles of Aedes aegypti super-production and human density in dengue transmission. PLoS Negl Trop Dis. 2012;6(8):e1799.
LaCon G, Morrison AC, Astete H, Stoddard ST, Paz-Soldan VA, Elder JP, et al. Shifting patterns of Aedes aegypti fine scale spatial clustering in Iquitos, Peru. PLoS Negl Trop Dis. 2014;8(8):e3038.
Yoon I-K, Getis A, Aldstadt J, Rothman AL, Tannitisupawong D, Koenraadt CJ, et al. Fine scale spatiotemporal clustering of dengue virus transmission in children and Aedes aegypti in rural Thai villages. PLoS Negl Trop Dis. 2012;6(7):e1730.
Schafrick NH, Milbrath MO, Berrocal VJ, Wilson ML, Eisenberg JN. Spatial clustering of Aedes aegypti related to breeding container characteristics in coastal Ecuador: implications for dengue control. Am J Trop Med Hyg. 2013;89(4):758–65.
Minakawa N, Mutero CM, Githure JI, Beier JC, Yan G. Spatial distribution and habitat characterization of anopheline mosquito larvae in western Kenya. Am J Trop Med Hyg. 1999;61(6):1010–6.
Mwangangi JM, Shililu J, Muturi EJ, Muriu S, Jacob B, Kabiru EW, et al. Anopheles larval abundance and diversity in three rice agro-village complexes Mwea irrigation scheme, Central Kenya. Malar J. 2010;9(1):228.
Minakawa N, Sonye G, Yan G. Relationships between occurrence of Anopheles gambiae sl (Diptera: Culicidae) and size and stability of larval habitats. J Med Entomol. 2005;42(3):295–300.
Matysiak A, Roess A. Interrelationship between climatic, ecologic, social, and cultural determinants affecting dengue emergence and transmission in Puerto Rico and their implications for Zika response. J Trop Med. 2017;2017:1–14.
Beaty BJ, Black W, Eisen L, Flores AE, García-Rejón JE, Loroño-Pino M, et al. The intensifying storm: domestication of Aedes aegypti, urbanization of arboviruses, and emerging insecticide resistance. Washington, D.C.: National Academies Press; 2016. Global Health Im-pacts of Vector-Borne Diseases: Workshop Summary
Centers for Disease Control & Prevention. Controlling mosquitoes at home 2018 [cited 2018 Oct 4]. Available from: https://www.cdc.gov/zika/prevention/controlling-mosquitoes-at-home.html.
Nakkhara P, Chongsuvivatwong V, Thammapalo S. Risk factors for symptomatic and asymptomatic chikungunya infection. Trans R Soc Trop Med Hyg. 2013;107(12):789–96.
Aguiar BS, Lorenz C, Virginio F, Suesdek L, Chiaravalloti-Neto F. Potential risks of Zika and chikungunya outbreaks in Brazil: a modeling study. Int J Infect Dis. 2018;70:20–9.
Thai KT, Binh TQ, Giao PT, Phuong HL, Hung LQ, Nam NV, et al. Seroprevalence of dengue antibodies, annual incidence and risk factors among children in southern Vietnam. Tropical Med Int Health. 2005;10(4):379–86.
World Health Organization. Dengue control: the mosquito 2019 [cited 2019 April 9]. Available from: https://www.who.int/denguecontrol/mosquito/en/.
López MAL, Dávila M, Canet M, López Y, Flores E, Dávila A, et al. Distribución de Aedes aegypti y Aedes albopictus en Guatemala 2016. Cienc Tecnol Salud. 2017;4(1):21–31.
Clark-Gil S, Darsie R. The mosquitoes of Guatemala, their identification, distribution and bionomics, with keys to adult females and larvae. Mosq Syst. 1983;15(3):151–284.
Pérez Esquivel WA. Estudio sobre la prevalencia de casos positivos de enfermedades metaxénicas, en la población de Coatepeque, Quetzaltenango, de Enero a Octubre del 2016 según registros del Ministerio de Salud Pública y Asistencia Social Galileo Universidad 2019 [cited 2019 April 11]. Available from: http://biblioteca.galileo.edu/tesario/handle/123456789/761.
Arevalo-Herrera M, Quiñones ML, Guerra C, Céspedes N, Giron S, Ahumada M, et al. Malaria in selected non-Amazonian countries of Latin America. Acta Trop. 2012;121(3):303–14.
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.
Centers for Disease Control and Prevention. About Malaria: Anopheles Mosquitoes 2018 [cited 2018 Oct 4]. Available from: https://www.cdc.gov/malaria/about/biology/index.html.
Nchoutpouen E, Talipouo A, Djiappi-Tchamen B, Djamouko-Djonkam L, Kopya E, Ngadjeu CS, et al. Culex species diversity, susceptibility to insecticides and role as potential vector of lymphatic filariasis in the city of Yaoundé, Cameroon. PLoS Negl Trop Dis. 2019;13(4):e0007229.
United Nations. World population prospects: the 2015 revision. United Nations Dep Econ Soc Aff. 2015;33(2):1–66.
Kraemer MU, Reiner RC, Brady OJ, Messina JP, Gilbert M, Pigott DM, et al. Past and future spread of the arbovirus vectors Aedes aegypti and Aedes albopictus. Nat Microbiol. 2019;4(5):854.
Franklinos LH, Jones KE, Redding DW, Abubakar I. The effect of global change on mosquito-borne disease. Lancet Infect Dis. 2019;19:e302–12.
Acknowledgements
The authors are grateful to Dra. Beatriz Santamarina, Dra. Analuisa Castro de Quan, Juan López, Juventino Cobon, Negly Castillo, Jorge Sincal, Nazario López, Jorge Luis López Díaz, Carlos Enrique Vásquez Ramos, and Jonathan Josue López Barrios, as well as local authorities, vector control personnel, and all field staff from both Coatepeque and Génova.
Funding
This study was supported by the Centers for Disease Control and Prevention (CDC) through cooperative agreement number 5U01GH001003–04. ZM was supported by the Fogarty International Center of the National Institutes of Health (NIH) under Award Number D43TW009343 and the University of California Global Health Institute (UCGHI). The content is solely the responsibility of the authors and does not necessarily represent the official views of UVG, CDC, NIH or UCGHI.
Author information
Authors and Affiliations
Contributions
ZJM contributed to study design, performed analyses, interpreted data, drafted the initial manuscript, and coordinated coauthor contributions. SS, JGJ, and CR collected data, interpreted results, and reviewed the manuscript. KCB and AL contributed to study design, data interpretation, and manuscript review. CCR contributed to study design, fieldwork coordination, data interpretation, and reviewing the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study was reviewed and approved by the Universidad del Valle de Guatemala Comité de Ética en Investigación. Participants provided written informed consent prior to data collection, according to the norms for Human Experimentation, Code of Ethics of the World Medical Association (Declaration of Helsinki, 1995).
Consent for publication
Not applicable.
Competing interests
The authors declare no conflict of interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1 Table S1.
Principal components factor analysis of household environment variables, Coatepeque and Génova, Guatemala, 2017 (n = 508). Table S2. Full model output of adjusted associations between geographical distances to paved roads and immature mosquito abundance, Poisson regression, Coatepeque and Génova, Guatemala, 2017. Table S3. Full model output of adjusted associations between geographical distances to highways and immature mosquito abundance, Poisson regression, Coatepeque and Génova, Guatemala, 2017. Table S4. Full model output of adjusted associations between geographical distances to houses or structures and immature mosquito abundance, Poisson regression, Coatepeque and Génova, Guatemala, 2017. Table S5. Associations between geographical distances to roads/structures and immature mosquito abundance, negative binomial regression, Coatepeque and Génova, Guatemala, 2017. Table S6. Mediation of distances to roads/structures and mosquito prevention measures on the association between environmental capital and the number of containers with any mosquito larvae or pupae per household, Coatepeque and Génova, Guatemala, 2017. Figure S1. Aerial view of communities in Coatepeque and Génova, Guatemala, 2017. Figure S2. Directed acyclic graphs of associations between geographical distances to roads and houses/structures and immature mosquito abundance. Figure S3. Cubic splines of associations between environmental capital and the number of containers with any mosquito larvae or pupae per household, Coatepeque and Génova, Guatemala, 2017.
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
Madewell, Z.J., Sosa, S., Brouwer, K.C. et al. Associations between household environmental factors and immature mosquito abundance in Quetzaltenango, Guatemala. BMC Public Health 19, 1729 (2019). https://doi.org/10.1186/s12889-019-8102-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12889-019-8102-5