Exploring the utility of social-ecological and entomological risk factors for dengue infection as surveillance indicators in the dengue hyper-endemic city of Machala, Ecuador.

The management of mosquito-borne diseases is a challenge in southern coastal Ecuador, where dengue is hyper-endemic and co-circulates with other arboviral diseases. Prior work in the region has explored social-ecological factors, dengue case data, and entomological indices. In this study, we bring together entomological and epidemiological data to describe links between social-ecological factors associated with risk of dengue transmission at the household level in Machala, Ecuador. Households surveys were conducted from 2014-2017 to assess the presence of adult Aedes aegypti (collected via aspiration) and to enumerate housing conditions, demographics, and mosquito prevention behaviors. Household-level dengue infection status was determined by laboratory diagnostics in 2014-2015. Bivariate analyses and multivariate logistic regression models were used to identify social-ecological variables associated with household presence of female Ae. aegypti and household dengue infection status, respectively. Aedes aegypti presence was associated with interruptions in water service and weekly trash collection, and household air conditioning was protective against mosquito presence. Presence of female Ae. aegypti was not associated with household dengue infections. We identified shaded patios and head of household employment status as risk factors for household-level dengue infection, while window screening in good condition was identified as protective against dengue infection. These findings add to our understanding of the systems of mosquito-borne disease transmission in Machala, and in the larger region of southern Ecuador, aiding in the development of improved vector surveillance efforts, and targeted interventions.

Introduction

The management of arthropod-borne viruses (arboviruses) is a fundamental task of public health agencies throughout much of Latin America and the tropics. Pathogens transmitted by the mosquito Aedes aegypti, such as dengue virus, impose some of the greatest growing disease burdens in this region [1]. Dengue has increased 30-fold over the past 50-years, globally and in South America. Over the past decade more than 19.6 million cases of dengue infection were reported in the Americas [2-4]. Dengue, like many arboviral infections, has limited options for clinical treatment and marketable vaccines, leaving vector control as a primary means of preventing and mitigating large outbreaks [5-7]. While there are many established protocols for designing and conducting public health vector control activities, these are subject to logistical constraints, resource availability, and fluctuating funding streams [6-11]. Development of targeted vector control strategies that focus on surveillance indicators, or known risk factors for transmission, offer a possible solution for implementing effective control measures in resource-limited communities (e.g., Lippi et al. 2020) [12]. Yet, social and ecological risk of exposure to arboviral infections is often understudied at fine spatial scales, such as individual households, at which many vector control interventions operate. Household-level features and behaviors can vary greatly within the built environment leading to varying risk of exposure to vectors, and thus, transmission [13-15].

Entomological and epidemiological surveillance systems can be integral components of public health vector control programs, where incoming data are key for informing management strategies and developing policy [11,16,17]. Vector abundance data are often used by public health agencies to trigger vector control responses, though the utility of directly inferring transmission risk from Ae. aegypti densities is difficult to establish, and varies widely across surveillance programs [17-19]. Beyond informing public agency responses, surveillance indicators can be used in the broader context of modeling risk by establishing correlative links between health outcomes and the social-ecological conditions in which they occur [20,21]. Ecological studies that establish such connections between dengue and factors such as population density, socioeconomic status, and housing conditions are routinely seen in the public health literature [14,20,22-24]. While observational studies such as these are inappropriate for inferring causality, they remain useful for informing current decision-making and guiding further epidemiological investigations.

In Ecuador, the national public health agency, the Ministerio de Salud Pública (MSP), oversees vector control, entomological surveillance, and epidemiological surveillance initiatives [25,26]. Ecuador’s southern coastal region is an area of historically high mosquito-borne disease activity [27,28]. El Oro province, located on the border with Peru, is consistently burdened with dengue fever outbreaks, where four serotypes of dengue virus are in circulation (DENV1-4), and more recently, chikungunya and Zika viruses [27]. The port city of Machala, with an estimated population of over 280,000, is Ecuador’s fourth largest city, the capital of El Oro province, and a major epicenter of commercial agricultural trade (Fig 1) [29]. Machala has notably high dengue fever incidence compared to other municipalities in Ecuador, where annual incidence rates reported to the MSP were 42.6 cases per 10,000 people in 2014 [27].

Fig 1

The study was conducted Machala, a city in the South American (A) country of Ecuador (B), located in El Oro province (B, shown in yellow). Households sampled during this study were located throughout the city of Machala, consisting of 460 households in 94 sampling clusters surveyed for Aedes aegypti mosquitoes (C), and 141 households in 33 sampling clusters surveyed for both Ae. aegypti mosquitoes and dengue infection status (D). Household locations were aggregated to census block for de-identification purposes in these figures. This figure was produced in ArcMap 10.6.1 (ESRI, Redlands, CA) using shapefiles freely available from the Natural Earth dataset ver. 4.1.0 (naturalearthdata.com) and georeferenced surveillance data provided by the MSP and edited by CAL.

Dengue is hyperendemic in Machala, where Ae. aegypti is prevalent and transmission can occur throughout the year [30,31]. The high incidence of dengue, large population, and longstanding history of vector control and surveillance activities make Machala an ideal setting for studying arboviral transmission systems. Several previous studies conducted in Machala have assessed risk factors for dengue fever outbreaks at the mesoscale (e.g. neighborhoods or census blocks), including larval vector densities, environmental drivers, sociodemographic characteristics, community perceptions, and service accessibility [12,14,25,30,32-34]. The findings of these studies are useful in a public health advisory capacity, providing information that guides decision-making and research. Nevertheless, studies conducted on data aggregated to administrative units face limitations in terms of guiding spatially explicit interventions on fine scales, and we cannot assume that neighborhood-level associations will be observable at the level of households, where people are exposed to mosquitoes [35-37], and at which scale interventions such as fogging (household insecticidal spraying), or individual prevention behaviors occur.

The objective of this study is to identify relationships between household social-ecological factors and adult Ae. aegypti presence and dengue cases, respectively, in Machala, Ecuador. We also test if Ae. aegypti presence is an indicator of household dengue fever cases. Household characteristics associated with the dengue fever transmission cycle are potentially identifiable targets for public health vector control intervention to reduce mosquitoes.

Methods

Cluster study design

The data used in this study were collected as part of a SUNY Upstate Medical University epidemiological study of dengue in Machala, Ecuador conducted from 2014–2017. The study protocol for collection of human case data was reviewed and approved by the Institutional Review Boards at the State University of New York (SUNY) Upstate Medical University, Cornell University, the Human Research Protection Office of the U.S. Department of Defense, the Luis Vernaza Hospital in Guayaquil, Ecuador, and the Ecuadorian MSP. Four sentinel clinics operated by the MSP were selected to enroll households into the study. These clinics were chosen based on historical patterns of dengue detection and available medical resources. A maximum of four index cases were randomly selected each week from clinically diagnosed dengue fever cases passively detected in sentinel healthcare facilities. A cluster study design was used to enroll additional participants. After identifying the households of index cases, up to four additional households were enrolled within a 200m radius of each index house, representing the approximate flight range of Ae. aegypti mosquitoes; this study protocol has been described in detail previously [38]. Each household was visited once during the study, within two days of enrollment, at which time household survey data, entomological data, and blood samples from individuals were collected. Locations of enrolled households (i.e. latitude and longitude coordinates) were recorded on site using Garmin handheld global positioning system (GPS) units.

Household survey–Field teams comprising local Ecuadorians visited each enrolled household, administering a survey to the head of household. Survey questions were developed based on previous studies from the region, which identified potential local factors of interest as they relate to risk of dengue transmission [30]. The survey tool used in this study can be found in the supporting information from Kenneson et. al [39,40]. The survey was used to collect information on household demographics, socioeconomic status, access to municipal services, and human behaviors that may influence the presence of mosquitoes (e.g. water storage practices and use of pesticides). Field teams also visually assessed housing conditions and structures, collecting data on patio conditions, window screening, ventilation method (i.e. fan or air conditioning unit), and construction materials. In addition to presence, the quality of certain structures was categorized by technicians on site, where features were considered to be in “poor” condition if they could not exclude mosquitoes from the household (e.g. window screens with visible holes) or provided ovipositional sites for Ae. aegypti mosquitoes (e.g. patios with collected garbage). Patios were additionally categorized as “shaded” if 50% of more of the area was covered by large tree canopies or greenhouse cloth. General housing condition was categorically ranked by technicians on site, ranging from “good” housing, which consisted of new construction in good repair (e.g. no holes in walls, well-maintained roof, windows and doors that shut, etc.), to “poor” conditions where structures were not maintained (e.g. older buildings in disrepair, observable damage to structure, holes in walls, unpainted exteriors, etc.). The complete list of household variables used in this study, with definitions, is provided in S1 Table.

Entomological survey

Field teams of entomology technicians visited households enrolled in the cluster study from 2014–2017. Households were surveyed for mosquitoes within two days of enrollment to assess abundance and presence of Ae. aegypti. Sampling consisted of backpack aspirator (Prokopack) collection of adult mosquitoes in and around households [41]. The entomological sampling protocol consisted of one technician operating the backpack aspirator and sampling the intradomicile (i.e. within the home) for 20 minutes and the peridomicile (i.e. courtyard or patio associated with home) for 10 minutes. Each room within a given household was sampled, starting at the floor, sampling under furniture, and working up to the ceiling. Adult mosquitoes were stored in coolers in the field, then identified, counted, and sorted by sex and species (i.e. Ae. aegypti versus other) at the entomological lab at the Technical University of Machala. Abundance counts of female Ae. aegypti for each year of the study were aggregated by households into a binary (i.e. presence/absence) dataset for statistical modeling (S1 Data).

Human case data

Data on dengue infection status were collected for enrolled households from 2014–2015; diagnostic results were not available for 2016–2017. Prior to initiating the study, all participants, or adult legal guardians, engaged in a written informed consent process. Cases detected via passive clinical surveillance were informed that they may be randomly selected to participate in the cluster study prior to signing the informed consent. The study included children (> 6 months) to adults who were evaluated for DENV infection.

Laboratory diagnostic tests for DENV were run on blood samples collected from study participants. These included NS1 rapid strip tests (Panbio Dengue Early Rapid), ELISA assays for NS1 (Panbio Dengue Early ELISA), IgM immunoassays (Panbio Dengue Capture IgM), and RT-PCR. Full details of the diagnostic testing used in this study are outlined in Stewart-Ibarra et al. (2018) [27]. Household members were considered positive for DENV infection if at least two diagnostic tests were positive. Case counts for each given year of the study were aggregated by households into a binary dataset for statistical modeling.

Statistical analysis

Statistical analyses and modeling were performed in R (ver. 3.6.2). Intraclass correlation coefficient values were calculated for the outcomes of Ae. aegypti presence and DENV infections, and indicated that household clusters did not impact observed variance in the data [42]. Bivariate analysis of social-ecological survey variables was conducted via Fisher’s exact test, testing for significant factors influencing the household presence of female Ae. aegypti mosquitoes or dengue cases, respectively. Statistically significant variables were used as predictors in a generalized linear model (GLM) framework, specifying a logistic modeling distribution (GLM, family = binomial, link = logit), to assess combined household-level influences on outcomes of interest (i.e. mosquito presence or presence of dengue cases).

Variance inflation factors (VIF) were calculated to evaluate multi-collinearity in household survey variables, with values below 10 indicating low collinearity [43]. Stability of logistic regression models was assessed via condition numbers (k), where values below 30 indicate stable models [44]. Global Moran’s I with inverse distance weighting was performed in ArcMap (ver. 10.6.1, ESRI, Redlands, CA, USA) on model residuals to test for presence of spatial autocorrelation, a violation of GLM assumptions [45].

Results

Ninety-four cases with lab-confirmed acute symptomatic DENV infections were randomly selected for this study via passive surveillance at MSP sentinel facilities during the 2014–2017 study period. These index cases were used to subsequently enroll an additional 366 households into the cluster study, for a total of 460 households in Machala, spanning 94 sampling clusters surveyed for household risk factors and mosquito presence during the study (Table 1, Figs 1C and 2). Aedes aegypti mosquitoes were collected in 54.1% of households (n = 249), and 43.9% of households (n = 202) had female Ae. aegypti present.

Table 1

Summary of social-ecological variables collected from 2014–2017 household surveys in Machala, Ecuador.

Dengue survey households represent a subset of households with valid diagnostic test results and complete household and entomological survey data.

Parameter	Households (n = 460)	%	Households (n = 141)	%
	Entomological Survey		Dengue Survey
Housing conditions
Good Condition Housing	144	32.0%	55	39.0%
Poor Condition Housing	55	12.0%	16	11.3%
Cane House Construction	32	7.2%	9	6.4%
Wood House Construction	9	2.0%	2	1.4%
House is Rented	40	9.5%	15	10.6%
Piped Water in Household	320	76.0%	100	70.9%
Interruptions in Piped Water	245	53.3%	78	55.3%
Municipal Garbage Collection	400	95%	130	92.2%
Municipal Sewage	356	84.8%	119	84.4%
Septic Tank	54	12.8%	20	14.2%
Air Conditioning	30	7.1%	13	9.2%
Uses a Fan	164	39.0%	70	49.6%
Screens on Windows	102	22.8%	44	31.2%
Window Screen in Good Condition	88	19.7%	35	24.8%
Access to Paved Roads	278	62.2%	68	48.2%
Standing Water Present	213	49.5%	72	51.1%
Adjacent to Abandoned Housing	119	26.6%	49	34.8%
Patio Present	367	82.1%	111	78.7%
Patio in Bad Condition	81	18.1%	21	14.9%
Patio Shaded (> 50%)	56	12.5%	19	13.5%
Household Demographics and Practices
Head of Household (HOH) Employed	359	85.7%	120	85.1%
HOH male	323	76.7%	103	73.0%
HOH Earns Less than Minimum Wage	74	17.6%	11	7.8%
HOH Has Secondary Education	161	38.2%	57	40.4%
Stores Water	283	67.2%	96	68.1%
Know mosquitoes transmit dengue	372	88.8%	124	87.9%
Know standing water produces mosquitoes	402	95.9%	136	95.7%
Uses Larvicide	98	23.3%	47	33.3%
Uses Abate	151	35.9%	2	1.4%

Fig 2

Diagram of household enrollment and data collection for cluster study design in Machala, Ecuador.

Housing conditions vary considerably in Machala, with approximately one third (32%) of surveyed housing classified in good condition, and 12% of housing structures classified in poor condition. Access to municipal services among sampled households was high, with 62.2% accessible by paved roads, 76.0% having access to piped water, 84.8% with access to sewage, and 95.0% receiving garbage collection services. The head of household was employed in 85.3% of surveyed households, and 17.6% of heads of household earned less than the minimum wage.

Four hundred study participants from enrolled households were tested for dengue infection status during 2014–2015, and 311 of these had valid diagnostic test results and complete household and entomological survey data. Infection status during 2014–2015 was aggregated to the household level for 141 households in 33 sampling clusters (Table 1, Fig 1D). Of these households, 29.08% had one or more household members positive for dengue. All four DENV serotypes were detected in Machala during the study period, and of the thirty-two households tested for dengue serotype, DENV-2 (56.3%) and DENV-1 (40.6%) were the most prevalent infections. Aedes aegypti mosquitoes were collected in 57.4% of the subset of households (n = 81) that were tested for dengue infection, and 49.6% of these households had female Ae. aegypti present (n = 70).

Housing conditions, demographics, and practices were statistically compared for households where female Ae. aegypti were found, or were not found, in the entomological survey (Table 2). Significant (p<0.05) factors positively associated with mosquito presence were interruptions in piped water availability, municipal garbage collection, access to paved roads, presence of standing water, and having a patio in poor condition. Home air conditioning, window screening, window screening in good condition, and adjacency to abandoned housing were negatively associated with female Ae. aegypti presence.

Table 2

Social-ecological factors in households with versus without female Aedes aegypti present.

Parameter	Households (n = 202)	%	Households (n = 258)	%	p-value
	Ae. aegypti Present		Ae. aegypti Absent
Housing conditions
Good Condition Housing	53	26.24%	91	35.27%	0.324
Poor Condition Housing	28	13.86%	27	10.47%	0.313
Cane House Construction	9	4.46%	23	8.91%	0.065
Wood House Construction	3	1.49%	6	2.33%	0.737
House is Rented	20	9.90%	20	7.75%	0.617
Piped Water in Household	147	72.77%	173	67.05%	0.731
Interruptions in Piped Water	125	61.88%	120	46.51%	0.001
Municipal Garbage Collection	187	92.57%	213	82.56%	0.013
Municipal Sewage	166	82.18%	190	73.64%	0.278
Septic Tank	20	9.90%	34	13.18%	0.241
Air Conditioning	5	2.48%	25	9.69%	0.001
Uses a Fan	84	41.58%	80	31.01%	0.057
Screens on Windows	34	16.83%	68	26.36%	0.013
Window Screen in Good Condition	28	13.86%	60	23.26%	0.009
Access to Paved Roads	137	67.82%	141	54.65%	0.008
Standing Water Present	108	53.47%	105	40.70%	0.010
Adjacent to Abandoned Housing	42	20.79%	77	29.84%	0.024
Patio Present	167	82.67%	200	77.52%	0.321
Patio in Bad Condition	44	21.78%	37	14.34%	0.049
Patio Shaded (> 50%)	26	12.87%	30	11.63%	0.775
Household Demographics and Practices
Head of Household (HOH) Employed	161	79.70%	198	76.74%	0.679
HOH male	145	71.78%	178	68.99%	0.730
HOH Earns Less than Minimum Wage	35	17.33%	39	15.12%	0.797
HOH Has Secondary Education	75	37.13%	86	33.33%	0.763
Stores Water	125	61.88%	158	61.24%	0.532
Know mosquitoes transmit dengue	168	83.17%	204	79.07%	0.645
Know standing water produces mosquitoes	180	89.11%	222	86.05%	0.346
Uses Larvicide	42	20.79%	56	21.71%	0.643
Uses Abate	72	35.64%	79	30.62%	0.477

P-values in bold denote statistical significance found via Fisher’s exact test

Statistically significant variables from bivariate analyses were put into a multivariate logistic regression model for female Ae. aegypti presence (Table 3). Water interruptions (OR = 1.67) and garbage collection (OR = 3.29) remained significant risk factors for household Ae. aegypti, while air conditioning (OR = 0.28) was associated with reduced presence. Other variables were no longer significant in the multivariate model. Low condition number indicated model stability (k = 7.15).

Table 3

Logistic regression model of household female Aedes aegypti presence in Machala, Ecuador.

Model	Estimate	SE	P-Value	OR	95% CI
Intercept	-1.81	0.63	0.004	-	-
Water Interruptions	0.51	0.21	0.017	1.67	1.10–2.55
Trash Collection	1.19	0.59	0.043	3.29	1.13–11.98
Air Conditioning	-1.27	0.52	0.014	0.28	0.09–0.72
Screens Present	-0.13	0.34	0.697	0.88	0.45–1.71
Screens in Good Condition	-0.40	0.36	0.263	0.67	0.33–1.35
Paved Roads	0.45	0.23	0.053	1.56	0.10–2.46
Standing Water Present	0.25	0.23	0.264	1.29	0.83–2.01
Adjacent to Abandoned Property	-0.37	0.25	0.141	0.69	0.42–1.13
Patio in Poor Condition	0.23	0.29	0.426	1.26	0.72–2.22

k = 7.15, pseudo R2 = 0.09

Housing conditions, demographics, and practices were statistically compared for households with DENV infections versus households without DENV infections (Table 4). Significant (p<0.05) factors positively associated with household DENV infections were having a shaded patio and HOH employment. Households with window screening in good condition were negatively associated with presence of DENV infections. Female Ae. aegypti presence in households was not significantly associated with presence of DENV infections.

Table 4

Social-ecological factors in households with versus without dengue infections.

Parameter	Households (n = 41)	%	Households (n = 100)	%	P-value
	Dengue Positive		Dengue Negative
Housing conditions
Good Condition Housing	11	26.83%	44	44.0%	0.061
Poor Condition Housing	7	17.07%	9	9.0%	0.240
Cane House Construction	5	12.20%	4	4.0%	0.122
Wood House Construction	1	2.44%	1	1.0%	0.499
House is Rented	6	14.63%	9	9.0%	0.371
Piped Water in Household	26	63.41%	74	74.0%	0.225
Interruptions in Piped Water	19	46.34%	59	59.0%	0.194
Municipal Garbage Collection	37	90.24%	93	93.0%	0.730
Municipal Sewage	32	78.05%	87	87.0%	0.206
Septic Tank	8	19.51%	12	12.0%	0.290
Air Conditioning	6	14.63%	7	7.0%	0.120
Uses a Fan	18	43.90%	52	52.0%	0.459
Screens on Windows	8	19.51%	36	36.0%	0.072
Window Screen in Good Condition	5	23.81%	30	30.0%	0.032
Access to Paved Roads	17	41.46%	51	51.0%	0.355
Standing Water Present	23	56.10%	49	49.0%	0.464
Adjacent to Abandoned Housing	18	43.90%	31	31.0%	0.174
Patio Present	33	80.49%	78	78.0%	0.824
Patio in Bad Condition	8	19.51%	13	13.0%	0.434
Patio Shaded (> 50%)	10	24.39%	9	9.0%	0.027
Household Demographics and Practices
Head of Household (HOH) Employed	39	95.12%	81	81.0%	0.037
HOH male	31	75.61%	72	72.0%	0.835
HOH Earns Less than Minimum Wage	2	4.88%	9	9.0%	0.510
HOH Has Secondary Education	17	41.46%	40	40.0%	1.000
Stores Water	27	65.85%	69	69.0%	0.843
Know mosquitoes transmit dengue	38	92.68%	86	86.0%	0.395
Know standing water produces mosquitoes	40	97.56%	96	96.0%	1.000
Uses Larvicide	14	34.15%	33	33.0%	1.000
Uses Abate	1	2.44%	1	1.0%	0.499
Entomological Survey
Female Aedes aegypti present	22	53.66%	48	48.0%	0.582

P-values in bold denote statistical significance found via Fisher’s exact test

Statistically significant variables from bivariate analyses were put into a multivariate logistic regression model for presence of household DENV infection (Table 5). Head of household employment (OR = 4.96) and households with shaded patios (OR = 3.81) were significant risk factors for household DENV infection. Window screening in good condition (OR = 0.30) was protective against presence of DENV infections. Low condition number indicated model stability (k = 4.75).

Table 5

Logistic regression model of household presence of dengue cases in Machala, Ecuador.

Model	Estimate	SE	P-Value	OR	95% CI
Intercept	-2.30	0.78	0.003	-	-
Screens in Good Condition	-1.20	0.55	0.028	0.30	0.09–0.82
Patio Shaded	1.34	0.54	0.014	3.81	1.33–11.47
HOH Employed	1.60	0.79	0.043	4.96	1.28–33.30

k = 4.75, pseudo R2 = 0.11

Intraclass correlation (ICC) values for both entomological and human case datasets were extremely low (i.e. near or equal to zero), indicating that the clustered study design did not meaningfully contribute to variance in the data. Despite the clustered arrangement of households in the study design, global clustering was not detected in model residuals for mosquito presence (Moran’s I = 0.04, p-value = 0.07) or dengue presence (Moran’s I = -0.14, p-value = 0.15), indicating no impact of spatial autocorrelation on model results.

Discussion

This study provides a household-level assessment of dengue fever risk in Machala, Ecuador, a city with historically high burden of the disease. Indicators of mosquito presence are at times used by public health agencies as a proxy for risk of disease exposure and basis for vector control operations, particularly in absence of active human surveillance data [2,6]. This is an intuitive approach, as widespread surveillance of human cases is resource intensive, and vector presence is a necessary component of transmission cycles for pathogens like dengue virus [46,47]. Yet, despite the close association of Ae. aegypti with built urban and human environments, we could not find a distinct link between adult mosquito presence and dengue risk in Machala. Female Ae. aegypti presence was not identified as a significant household risk factor for dengue infection status. Interestingly, this finding parallels a previous study conducted by Stewart-Ibarra et. al during the 2010 dengue outbreak in Machala, where neighborhood larval mosquito indices were not predictive of dengue cases in the study’s best-fit model [30].

The lack of association between mosquito presence and dengue infections indicated in this study may reflect a disconnect between household conditions and locations where transmission is occurring, particularly in the context of hyperendemic dengue transmission in Machala. Studies conducted in other Latin American cities with hyperendemic dengue transmission have shown that within cities, transmission risk is influenced by neighboring census tracts, and no specific locations within a city were driving dengue cases [48,49]. Daily human movements and social connections within hyperendemic spaces can also drive dengue transmission, particularly on fine spatial scales [50,51]. Given the high prevalence of the disease in the community, risk of transmission in Machala may be more dependent on social and behavioral factors that increase exposure to infective mosquito bites, rather than mere presence or absence of the vector in the immediate vicinity of the household. Although not informative for dengue risk at the household level in this study, municipal and sub-regional entomological surveillance is still a cornerstone of successful public health vectors control programs, where mosquito counts are necessary to set goals and proactively trigger the deployment of vector control services ahead of epidemic peaks [6,52,53]. However, in locations with hyperendemic arboviral transmission, like Machala, surveillance of human cases may provide better primary information for guiding policy and developing direct community interventions on fine spatial scales.

We did not use mosquito abundance as a metric in modeling in this study due to low counts of adult mosquitoes (e.g. <10) recorded in household surveys via aspiration. Backpack aspiration is a logistically feasible surveillance method for many vector control agencies in Latin America. This method of sampling adult mosquitoes is relatively accessible, and there exists evidence of high sensitivity for detecting the presence of mosquitoes in households [18,41,54]. Regardless, indoor household Ae. aegypti densities can be notoriously low, and samples collected with backpack aspirators are often highly dependent on factors such as user skill and complexity of the indoor environment [55-57]. Some studies suggest that alternative trapping methods, such as the use of passive BG-Sentinel traps, may be better suited for targeting Ae. aegypti [9,56,57]. In contrast with binary presence data used here, adult mosquito abundance data have been shown in some instances to correlate with risk of dengue infection, particularly in longitudinal studies that can account for temporal variability and lagged effects [18,58].

Presence of Ae. aegypti was not predictive of dengue infections in Machala, yet mosquitoes are a necessary component of the dengue transmission cycle. Identification of factors that promote mosquito presence can provide valuable information for establishing mosquito control targets, validating or refining current MSP advisories and educational initiatives for the control of container-breeding mosquitoes [2,26]. We identified several factors related to water availability and housing conditions in our bivariate analyses that were associated with vector presence. These included interruptions in water service, presence of standing water around the household, patios in poor condition, municipal garbage collection, and paved road access. Air conditioning, proximity to abandoned housing, window screening, and window screening in good condition were negatively associated with mosquito presence. We did not find associations between mosquito presence and household demographics and practices, suggesting that measured social risk factors and personal protective behaviors for dengue did not reliably indicate local presence of adult female Ae. aegypti in Machala. Two risk factors, water interruptions and garbage collection, and one protective factor, household air conditioning, remained significant in the logistic regression model of female Ae. aegypti presence.

Associations between household characteristics and mosquito presence may be leveraged into actionable abatement strategies by the local public health agency, improving vector control and outreach efforts in Machala. Further, the overlap of risk factors with other studies may indicate candidate targets for regionally generalized mosquito control recommendations. For example, water availability and storage practices have previously been associated with immature mosquito presence in different cities, and at different spatial scales, in coastal Ecuador [14,59]. Protective measures may also be regionally transferrable. Factors that eliminate the need to open doors and windows for home ventilation, such as air conditioning, are typically found to be protective against exposure to mosquitoes across many systems [60-62]. Counterintuitive relationships observed at local scales, such as the positive association between weekly garbage collections and Ae. aegypti presence in Machala, are also informative to local mosquito control. Frequent municipal trash collection is generally associated with reduced debris, which translates to reduced ovipositional sites for container-breeding mosquitoes [63]. However, a census tract-level correlation between garbage collection and dengue fever was also observed in Guayaquil, Ecuador, suggesting that highly developed urban centers may provide ample mosquito production sites, and elevated disease risk, regardless of access to municipal services [20].

Household-level variables also influenced dengue infection status of homes in Machala. Window and house screening that was in good condition was identified as a significant protective factor against the presence of household dengue infections. Screened entryways can reduce exposure to potentially infective mosquitoes by establishing barriers to access of vectors to households [64]. The use of physical barriers to exclude insects is a cornerstone of integrated pest management, and has been effectively used to reduce dengue fever risk in many locations [64-67]. Yet, in Machala we found that presence of screening itself was not significantly protective against dengue. Screening in disrepair, where survey technicians reported visible damage or holes in screens, may not sufficiently exclude mosquitoes from households. We also identified the presence of shaded patios as a dengue risk factor in our study, a feature which has been linked to transmission risk in other studies in the region, promoting favorable habitat and oviposition sites for container-breeding mosquitoes [14,39].

The lack of association between household dengue cases and indicators of water availability found in this study is notable, and was unexpected. Frequent water interruptions promote water storage behaviors, such as collecting water in uncovered drums and basins, which increase suitable oviposition sites that can be exploited by Ae. aegypti. Mosquito reproduction and abundance are thus impacted by water storage practices, and interruption in water service was a significant predictor in the model for mosquito presence in this study. Access to piped water and water storage have been linked to both dengue infections and the proliferation of Ae. aegypti across multiple spatial scales, in Machala and other coastal Ecuadorian cities [14,20,30,31,68,69].

Head of household employment was a risk factor for dengue presence in this study. This relationship is counterintuitive, as HoH employment is often used as a proxy for higher socioeconomic status, which is generally associated with improved housing quality and lower arboviral disease risk [6,70]. However, occupational exposures to mosquito-borne diseases, where infective bites occur at the workplace, are known to vary across employment fields [71,72]. Although the survey tool used in this study represents a major step forward in the collection of fine-scale data regarding risk of dengue transmission in Ecuador, the questions did not necessarily capture the nuance of some risk elements. Events occurring beyond the home, tied to employment, such as those associated with occupation type, daily commuting patterns, or other work-related travel, were not captured by the survey tool. Such an association was reported in the Galapagos, where travel patterns were associated with self-reported dengue illness [45].

A previous study, conducted by Kenneson et al. (2017), also examined household dengue fever risk in Machala during 2014–2015 [39]. In the current study, we restricted analyses to households with complementary dengue case and entomological survey data, and used a more conservative and proximal dengue case definition, requiring at least two positive diagnostic tests, excluding IgG immunoassays, indicative of longer term exposure [4]. The 2017 study, with broader case definitions and a larger group of households, also found that shaded patios were significant risk factors for dengue fever, but additionally identified adjacency to abandoned properties and frequent water interruptions as significant risk factors. We found that frequent water interruption was a risk factor in our multivariate model for Ae. aegypti presence, but was not significant in the model for dengue case presence. While Kenneson et al. did not examine overall housing conditions, making direct comparisons unfeasible, they instead demonstrated that preventive actions and perceptions were negatively associated with dengue fever outcomes. While these studies leverage overlapping sampling events, the differences in objectives, case definitions, and drivers explored, indicate a need for further examination of the implications for connections in the system between the transmission conditions, the scale of human behavior, and household impact, and which components of the lived environment interact with the vector-human transmission cycle.

Household-level data are useful in advancing our understanding of localized disease transmission and exposure risk, yet there are caveats that must be considered before findings are used in an applied capacity. The statistical models used in this study are correlative, and we therefore must avoid assigning causality to any significant relationships without further investigation. Furthermore, data collected via household surveys certainly provide information on a fine scale, but they do not capture the full social-ecological transmission environment. This is exemplified by the lack of information concerning potential workplace and school exposures. We also caution against reducing localized risk indicators into a dichotomy of “poor” versus “affluent” neighborhoods; this fails to account for nuances in risk due to personal behaviors, or exposure due to hyperendemic presence of mosquito vectors. Focusing solely on economic indicators of risk has the additional potential to generate stigma in communities, which can undermine public health initiatives [73].

The results of this study aid us in identifying areas for further study as potential targets for developing proactive mosquito control activities and intervention strategies in Machala. The indication of missing links between entomological exposure risk indicators at the household level, where vector control intervention is often carried out, and the indicators of dengue infection, point to a need for a more complete examination of the social-ecological transmission environment. Nonetheless, comparing our findings to previously published work, we have identified commonalities in risk factors across spatial scales. Thus, we indicate that some household-level risk factors, such as water interruptions, may be generalizable in a policy development capacity.

Conclusions

In this study, we sought to assess the utility of entomological surveillance, specifically backpack aspiration of adult Ae. aegypti mosquitoes, in the context of simultaneous epidemiological information about dengue, at the household level, in Machala, Ecuador. The use of backpack aspirators is an attractive option for entomological surveillance due to their relative accessibility. Yet, this method of surveillance faces a myriad of challenges including training technicians and gaining access to homes. Here, we found a mismatch between the entomological surveillance information for adult mosquitoes, the predictive capacity of those results in the context of risk, and identifiable social-ecological risk factors for dengue transmission. In situations where adult mosquito surveillance data are collected as an indicator of transmission risk, the utility of these data for identifying potential intervention targets should be evaluated. In Machala, where dengue is hyper-endemic, adult mosquito surveillance may not lead to major improvements in dengue risk assessment, particularly when compared to larval mosquito surveys and epidemiological surveillance conducted at broader scales.

Description of social-ecological variables collected in Machala via household surveys.

(DOCX)

Click here for additional data file.

Female Aedes aegypti presence detected in household entomological surveys conducted in Machala, Ecuador from 2014–2017.

(XLSX)

Click here for additional data file.