Collections by aspirators
A total of 2,874 mosquitoes were collected by aspirators (Table 1). Ae. aegypti and Culex spp. accounted for the majority of mosquitoes, 44.8% and 53.2% respectively, with Ae. albopictus contributing only 2.1% (59 mosquitoes, mainly from the rural sites). There were significantly more Ae. aegypti collected in the rural villages than the urban villages (959 vs. 328, χ2 = 61.8, df = 1, P < 0.001). In addition, the second urban village (Non Tan) had even lower abundance than the first urban village (Nong Hai) (105 mosquitoes vs. 223). More female than male Ae. aegypti were collected overall (708 vs. 579, χ2 = 7.99, df = 1, P = 0.005). The abundance of Culex spp. was also higher in rural villages (992 vs. 615, χ2 = 17.5, df = 1, P < 0.001), but there was significant variation between villages in the rural and urban sites (χ2 = 9.96, df = 2, P = 0.008).
Table 1
Number of mosquitoes (%) collected by mechanical battery-driven aspirator and sticky traps in rural and urban areas in northeastern Thailand in 2019.
Species | Methods | Areas | No. of mosquitoes | |
Female | Male | Sum | Total |
Ae. aegypti | Aspirator | Rural | 524 (74.0) | 435 (75.1) | 959 (74.5) | |
Urban | 184 (26.0) | 144 (24.9) | 328 (25.5) | |
Sum | 708 (100.0) | 579 (100.0) | 1,287 (100.0) | |
Sticky trap | Rural | 334 (64.1) | 280 (64.7) | 614 (64.4) | |
Urban | 187 (35.9) | 153 (35.3) | 340 (35.6) | |
Sum | 521 (100.0) | 433 (100.0) | 954 (100.0) | 2,241 |
Ae. albopictus | Aspirator | Rural | 36 (83.7) | 12 (92.3) | 51 (91.5) | |
Urban | 7 (16.3) | 1 (7.7) | 8 (8.5) | |
Sum | 43 (100.0) | 13 (100.0) | 59 (100.0) | |
Sticky trap | Rural | 0 | 4 (100.0) | 4 (100.0) | |
Urban | 0 | 0 | 0 | |
Sum | 0 | 4 (100.0) | 4 (100.0) | 63 |
Culex spp. | Aspirator | Rural | 515 (71.8) | 414 (52.3) | 929 (60.8) | |
Urban | 202 (28.2) | 378 (47.7) | 599 (39.2) | |
Sum | 717 (100.0) | 792 (100.0) | 1,528 (100.0) | |
Sticky trap | Rural | 293 (62.1) | 262 (57.7) | 555 (63.6) | |
Urban | 179 (37.9) | 192 (42.3) | 317 (36.4) | |
Sum | 472 (100.0) | 454 (100.0) | 872 (100.0) | 2,400 |
Overall | Aspirator | Sum | 1,468 (59.7) | 1,384 (60.8) | 2,874 (61.1) | |
| Sticky trap | Sum | 993 (40.3) | 891 (39.2) | 1,830 (38.9) | |
| | Total | 2,461 (100.0) | 2,275 (100.0) | 4,704 (100.0) | 4,704 |
Male and female Ae. aegypti were most abundant at intermediate heights in rooms where they were most prevalent (i.e., all rooms except the kitchen) (Fig. 2A). Overall, in decreasing order for abundance by height: 0.75–1.5 m > 0-0.75 m > 1.5 + m (all t statistics for the three among-height comparisons, t > 8.5, P < 0.001). In decreasing order of abundance were bedroom > bathroom = living room > kitchen (t statistics for the five significant comparisons t > 3.1, P < 0.01). Both males and females showed the same distribution in heights and rooms.
Culex mosquitoes were more abundant at the two lower heights than at the higher height (t statistics for the two significant comparisons t > 5.8, P < 0.001). There was no difference in the number of male and female Culex mosquitoes caught (χ2 = 0.01, df = 1, P = 0.903). Although female Culex were found predominantly in the bedroom and bathroom, males showed very different room and height distributions (Fig. 2B). The full data set summaries of numbers of mosquitoes caught according to place, time and height for Ae. aegypti and Culex spp. are shown in Additional file 1: Table S1 and Additional file 2: Table S2.
There were very few Ae. albopictus mosquitoes overall, but there was a significant association with more Ae. albopictus present in the living room and kitchen than in the bedroom and bathroom (living room vs. bathroom t = 7.39, P < 0.001; living room vs. bedroom t = 5.99, P < 0.001; kitchen vs. bathroom t = 4.24, P < 0.001). There was no significant association with height (χ22 = 0.59, P = 0.745).
The abundance of Ae. aegypti decreased from morning to afternoon in the bathroom, bedroom and kitchen but increased in the living room (Fig. 3A). These changes were significant for the female mosquitoes in the bedroom (t = 2.25, P < 0.05), the kitchen (t = 2.71, P < 0.01) and the living room (t = 2.42, P < 0.05). By contrast, the abundance of Culex spp. decreased from morning to afternoon in all rooms (χ2 = 15.24, df = 1, P < 0.001) (Fig. 3B). For Ae. albopictus, numbers decreased from morning to afternoon, but this was not significant (χ2 = 3.70, df = 1, P = 0.053).
There were more bloodfed Ae. aegypti collected than unfed (χ2 = 9.35, df = 1, P = 0.002), but they did not differentially distribute by height (χ2 = 0.47, df = 2, P = 0.791) or room (χ2 = 2.75, df = 3, P = 0.432) (Fig. 4A). By contrast, fewer bloodfed Culex spp. female mosquitoes were collected than unfed (χ2 = 86.6, df = 1, P < 0.001) (Fig. 4B), there was a distinctly different height distribution for bloodfed mosquitoes in the bedroom as compared to the other rooms. There were too few bloodfed Ae. albopictus for analysis.
Sticky traps
A total of 1,830 mosquitoes were collected by sticky traps (Table 1). Aedes aegypti and Culex spp. accounted for 52.1% and 47.7% respectively, and Ae. albopictus contributing only 0.2% (only 4 mosquitoes, all from the urban sites). The full data set summaries of numbers of mosquitoes caught according to place, time and height for Ae. aegypti and Culex spp. are shown in Additional file 3: Table S3 and Additional file 4: Table S4.
Aedes aegypti abundance, as measured by sticky traps, was higher in the rural areas than the urban areas (χ2 = 60.46, df = 1, P < 0.001), but did not differ further at the village level (χ2 = 0.92, df = 2, P = 0.63). There were more female than male mosquitoes (χ2 = 10.2, df = 1, P = 0.002), but no sex differences in height or room resting preference were found. There were significant differences in overall resting behaviour, irrespective of sex, for both rooms (χ2 = 195.4, df = 3, P < 0.001) and heights (χ2 = 212.02, df = 2, P < 0.001) (Fig. 5A). The mosquito abundance in rooms and by height was in decreasing order, bedroom > bathroom > living room > kitchen (all t statistics t > 3.5, P < 0.001) and 0.75-1.5m > 0-0.75m > 1.5 + m (all t statistics t > 7.0, P < 0.001), respectively. There was no interaction between room and height for mosquito abundance.
Culex spp. sticky trap abundances were higher in the rural areas (χ2 = 18.06, df = 1, P < 0.001), but did not differ further at the village level (χ2 = 0.92, df = 2, P = 0.63). The numbers of female and male Culex mosquitoes did not significantly differ (χ2 = 0.17, df = 1, P = 0.68). There was a small interaction effect between room and height (χ2 = 14.16, df = 6, P = 0.03), reflecting the relatively lower mid-height abundance in the bathroom as compared to the other rooms (Fig. 5B). Otherwise, there were no significant differences in abundance by height (χ2 = 2.01, df = 2, P = 0.38) or room (χ2 = 5.27, df = 3, P = 0.154).
Only four Ae. albopictus were collected and thus were not analysed.
Comparison of mosquito collections by an aspirator and sticky traps
Overall, in houses where mosquitoes were collected both by aspiration and sticky traps, the sticky traps collected more mosquitoes (Ae. aegypti: aspirator mean: 0.58, standard error of the mean (SEM) 0.13 vs. sticky trap mean: 3.98, SEM: 0.22. Culex spp. aspirator mean: 0.91, SEM: 0.19 vs. sticky trap mean: 3.86, SEM: 0.19). However, per sampling sticky traps were deployed for a lot longer (i.e. 3 times 7 days vs. 30 minutes). For the comparison, only aspirator mosquito collections from houses that had sticky traps deployed were used. Taking into account the differences between area, village, room and height, the aspirator method was more efficient per sampling time effort in collecting Ae. aegypti (χ2 = 1165.0, df = 1, P < 0.001) and Culex spp. (χ2 = 1372.6, df = 1, P < 0.001) mosquitoes. Aspirators collected a mean of 0.019 Ae. aegypti/minute and 0.029 Culex spp./minute, whereas sticky traps collected a mean of 1.3 x 10− 4 Ae. aegypti/minute and 1.28 x 10− 4 Culex spp./minute. Although the efficiency of capture was higher at all heights, for Ae. aegypti there was a significant interaction effect where aspiration was even better at higher heights (vs. lowest height, t = 1.98, P = 0.048; vs. intermediate height, t = 3.38, P < 0.001).
DENV positive mosquitoes
Of the 422 Ae. aegypti females assessed for DENV infection only five specimens (1.7%) were positive for the virus and too few to analyze with respect to their distribution in rooms and across heights. All DENV positives were, however, collected in the rural area. Two specimens were positive for DENV-1 (0.7%), one for DENV-3 (0.3%), and one had a mixed DENV-1 and DENV-3 infection (0.3%). One specimen was positive for DENV-1, DENV-2 and DENV-3 combined (0.3%).
Association of socio-environmental variables
Summary household information collected using the questionnaire is shown in Additional file 5: Table S5. There were several notable differences between rural and urban areas. Overall, there were more people per household in urban settings (urban mean: 4.08 standard deviation (SD): 1.71 vs. rural mean: 3.37, SD 1.58). Urban household heads were predominantly people with permanent jobs, shopkeepers and casual labourers (87%), whereas rural household heads were predominantly farmers and casual labourers. Almost no urban houses had any livestock and relatively few rural houses had livestock (15%). Approximately half of the rural houses had wooden walls (53%), whereas 82% of urban houses had plaster walls. 87% of rural houses had squat toilets vs. 67% of urban houses having bowl toilets. 87% of urban houses had window screens, whereas only 12% of rural houses had screens. Larval control either with temephos or via cleaning of containers was carried out more frequently (> 90% on a weekly basis) in urban areas than in rural areas (50% of houses at a monthly rate).
The P-values for the univariable analyses of the association of socio-environmental variables with Ae. aegypti or Culex spp. mosquito numbers are shown in Additional file: Table S6 and Additional file: Table S7. In urban settings, in the final adequate multivariable model, only use of repellent (Yes/No) was associated with total Ae. aegypti number and the Yes category was found to associate with increased numbers of mosquitoes (χ2 = 5.29, df = 1, P = 0.022). In rural settings, the following variables were found to be associated with higher Ae. aegypti number: use of temephos for larval control with low use (only every 3 months) was associated with higher numbers of mosquitoes (χ2 = 9.55, df = 3, P = 0.024); cement walls (χ2 = 8.62, df = 2, P = 0.014), outdoor toilets (χ2 = 5.75, df = 1, P = 0.017) and increasing number of rooms (χ2 = 8.91, df = 1, P = 0.003). Clothing location was also associated with mosquito numbers (χ2 = 12.78, df = 2, P = 0.002), with intermediate height of hung clothing being associated with more mosquitoes than lower or higher levels (intermediate vs. low levels: t = 2.19, P < 0.05; intermediate vs. high levels: t = 3.25, P < 0.01).
For Culex spp. mosquitoes, in urban settings the following variables were found to be associated with higher mosquito numbers: absence of open eaves (χ2 = 6.34, df = 1, P = 0.012), infrequent use of fogging (every 3 months) (χ2 = 25.97, df = 2, P < 0.001) and low levels of wind flow (χ2 = 9.83, df = 2, P = 0.008). In rural settings, use of fogging (χ2 = 16.9, df = 1, P < 0.001) was associated with lower numbers of mosquitoes and quantity of hung clothing (χ2 = 6.23, df = 2, P = 0.045; none/little vs. a lot, t = 2.46, P < 0.05) was associated with higher numbers of mosquitoes. Shopkeeping as the household occupation was also associated with higher mosquito numbers (χ2 = 11.29, df = 4, P = 0.024).
Geckos and mosquito abundance
A total of 297 geckos were caught, of which 133 (45%) were in the urban area and 164 (55%) in the rural area. Close to 28% of the geckos were collected in bathrooms (n = 82), followed by bedrooms (n = 75; 25%), living rooms (n = 74; 25%), and kitchens (n = 66; 22%). Geckos were caught on the lower (n = 102; 34%), intermediate (n = 94; 32%) and upper levels (n = 101; 34%).
Univariable fitting of the gecko number caught revealed no significant association of gecko number with either the Ae. aegypti or Culex spp. number caught on sticky traps. Using the total number of mosquitoes caught by aspiration again revealed no association of gecko number with Ae. aegypti or Culex spp. mosquito number in rural settings. However, there were significant associations of mosquito numbers with gecko numbers in the urban sites. There was a negative association between gecko number and total number of Ae. aegypti caught by aspiration (χ2 = 12.61, df = 1, P < 0.001). By contrast, there was a positive association between gecko number and the total number of Culex spp. caught by aspiration (χ2 = 7.32, df = 1, P = 0.008). These associations remained in a multivariable analysis that included above-identified significant socio-environmental variables.