Next Article in Journal
Aggregation Site Choice by Gregarious Nymphs of the Desert Locust, Schistocerca gregaria, in the Sahara Desert of Mauritania
Next Article in Special Issue
Oviplate: A Convenient and Space-Saving Method to Perform Individual Oviposition Assays in Aedes aegypti
Previous Article in Journal
Nest Population Structure and Wood Litter Consumption by Microcerotermes indistinctus (Isoptera) in a Seasonally Dry Tropical Forest, Northeastern Brazil
Previous Article in Special Issue
Mosquito Innate Immunity
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 9
Issue 3
10.3390/insects9030098
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Increased Adult Aedes aegypti and Culex quinquefasciatus (Diptera: Culicidae) Abundance in a Dengue Transmission Hotspot, Compared to a Coldspot, within Kaohsiung City, Taiwan

by
Ka-Chon Ng
1,†,
Luis Fernando Chaves
2,3,†,
Kun-Hsien Tsai
1 and
Ting-Wu Chuang
4,*
1
College of Public Health, National Taiwan University, Taipei 10055, Taiwan
2
Instituto Costarricense de Investigación y Enseñanza en Nutrición y Salud (INCIENSA), Apartado Postal 4-2250, Tres Ríos, Cartago, Costa Rica
3
Programa de Investigación en Enfermedades Tropicales (PIET), Escuela de Medicina Veterinaria, Universidad Nacional, Apartado Postal 304-3000, Heredia, Costa Rica
4
Department of Molecular Parasitology and Tropical Diseases, School of Medicine, College of Medicine, Taipei Medical University, No. 250, Wuxing Street, Xinyi District, Taipei 11031, Taiwan
*
Author to whom correspondence should be addressed.
†
Contributed equally.
Insects 2018, 9(3), 98; https://doi.org/10.3390/insects9030098
Submission received: 20 June 2018 / Revised: 30 July 2018 / Accepted: 10 August 2018 / Published: 13 August 2018
(This article belongs to the Special Issue Integrative Mosquito Biology: From Molecules to Ecosystems)
Browse Figures
Versions Notes

Abstract

:
The assumption that vector abundance differences might drive spatial and temporal heterogeneities in vector-borne disease transmission is common, though data supporting it is scarce. Here, we present data from two common mosquito species Aedes aegypti (Linnaeus) and Culex quinquefasciatus Say, biweekly sampled as adults, from March 2016 through December 2017, with BG-sentinel traps in two neighboring districts of Kaohsiung City (KC), Taiwan. One district has historically been a dengue transmission hotspot (Sanmin), and the other a coldspot (Nanzih). We collected a total 41,027 mosquitoes, and we found that average mosquito abundance (mean ± S.D.) was higher in Sanmin (Ae. aegypti: 9.03 ± 1.46; Cx. quinquefasciatus: 142.57 ± 14.38) than Nanzih (Ae. aegypti: 6.21 ± 0.47; Cx. quinquefasciatus: 63.37 ± 8.71) during the study period. In both districts, Ae. aegypti and Cx. quinquefasciatus population dynamics were sensitive to changes in temperature, the most platykurtic environmental variable at KC during the study period, a pattern predicted by Schmalhausen’s law, which states that organisms are more sensitive to small changes in environmental variables whose average value is more uncertain than its extremes. Our results also suggest that differences in Ae. aegypti abundance might be responsible for spatial differences in dengue transmission at KC. Our comparative approach, where we also observed a significant increment in the abundance of Cx. quinquefasciatus in the dengue transmission hotspot, suggests this area might be more likely to experience outbreaks of other vector borne diseases and should become a primary focus for vector surveillance and control.

1. Introduction

The yellow fever mosquito, Aedes aegypti (Linnaeus), and the southern house mosquito, Culex quinquefasciatus Say, are two of the most common urban mosquitoes in tropical and subtropical environments around the world [1,2,3]. Ae. aegypti is the main vector of dengue, chikungunya, and Zika in subtropical and tropical regions [4,5,6], while Cx. quinquefasciatus is a vector of worms causing lymphatic filariasis [7], West Nile Virus, WNV [8], and St Louis Encephalitis Virus [9], and both species are major causes of nuisance to humans [10]. Cx. quinquesfaciatus is one of the few species where impacts of global warming on its overwintering patterns have been observed in subtropical environments [11], while Ae. aegypti is a species whose distribution seems to be expanding as a result of global warming [12,13,14]. With changing climatic patterns, increased movement of people and goods, and unplanned urbanization, the transmission risk of vector-borne diseases has seen an unprecedented growth across the globe [15,16,17]. Dengue fever is one of the vector-borne diseases that has seen a significant worldwide geographical expansion [18,19] and the most important vector-borne disease in Taiwan [20].
Mosquito abundance is a key parameter for the infection risk of vector-borne diseases [21,22,23]. Various climatic factors impact mosquito ecology and pathogen development and transmission [24,25]. An ubiquitous pattern is that high temperatures reduce mosquito developmental time [26,27,28,29] and also reduce the extrinsic incubation period of pathogens—for example, dengue virus in Ae. aegypti [30] and WNV in Culex spp. mosquitoes [8,31]. Depending in the environmental context, high temperatures can also drive mosquito abundance outbreaks, i.e., sudden increases in the abundance of an insect [32,33]. Relative humidity might influence the survival and flight ability of vectors [34,35,36]. The effect of precipitation, by contrast, depends on the environmental context where mosquitoes live [24]. For example, in dry environments moderate rainfall is required to create more habitats for vectors [37], rainfall can drive vector abundance outbreaks in wet environments [38,39], and rainfall can drive the recruitment of larvae from diapausing eggs [40]. However, rainfall impacts are highly nonlinear [39], given that extreme rainfall might wash-out larvae from aquatic habitats depending on the mosquito species [41]. In addition, droughts might trigger water accumulation in artificial containers by humans, which can increase the recruitment of adult mosquito vectors [42,43,44].
Although no single mosquito sampling method can capture all mosquito species in a given area [45,46,47] nor satisfactorily allow the estimation of mosquito population size [48], systematic sampling with a single method could help to understand the impact of environmental fluctuations on mosquito species abundance [24], and vector infection at a given place [49,50]. Moreover, insights about mosquito population dynamics [51,52,53] and pathogen transmission [54,55] become increasingly robust as observations are not restricted to a single season. In addition, systematic sampling could also help to understand differences in the entomological risk for disease transmission over space and time [50,55,56].
Kaohsiung City (KC) has been the epidemic center of dengue outbreaks in 2014 (14,971 confirmed cases) and 2015 (19,769 confirmed cases) in southern Taiwan [57,58,59]. Our previous work in southern Taiwan demonstrated important associations between both regional and local climatic fluctuations on dengue transmission [20] and spatio-temporal analyses of dengue incidence has shown the existence of persistent dengue transmission hotspots and coldspots in KC [60,61,62]. One hypothesis for spatial heterogeneity in dengue transmission is that hot and cold spots are associated with different patterns of mosquito abundance and response to environmental fluctuations, which can be specific to dengue vectors or shared with other mosquito species commonly co-occurring with dengue vectors. To test this hypothesis, we biweekly collected mosquitoes during 2016 and 2017 with BG-sentinel traps [63], placed in two highly urbanized districts of KC, one part of a dengue transmission hotspot and one of a dengue transmission coldspot. We used time series analysis to evaluate differences in the association between vector abundance and weather fluctuations in Ae. aegypti, the dominant dengue vector in southern Taiwan, and Cx. quinquefasciatus, a common urban nuisance mosquito, and a species also competent to transmit several pathogens.

2. Materials and Methods

2.1. Study Area and Weather Data

Kaohsiung City (KC) is a city with 2.8 million habitants and a density of 942 people/km2 (Department of Household Registration, Ministry of Interior). KC is the main dengue epidemic center in Taiwan and accounts for more than 70% of dengue cases in the whole country (Taiwan Centers for Disease Control). Two highly urbanized (i.e., with over 40% of land cover used for housing and roads) districts of KC—Sanmin and Nanzih—were selected for the study. Sanmin district has historically been a dengue transmission hotspot with a high annual incidence rate (209.5 per 100,000). In contrast, Nanzih has been a dengue transmission coldspot, where dengue incidence rate is about 47.9 per 100,000. The location of the districts and sampling locations within each district are shown in Figure 1, and coordinates and pictures of the places where mosquito traps were located are presented in Supplementary Table S1 [online only]. For each sampling location, we evaluated the percent land cover in a circular buffer of 500 m radius around each trap, the maximum recorded dispersal distance of Ae. aegypti in urban landscapes [64], and a radius were most host seeking Cx. quinquefasciatus occur in relation to their larval habitats of origin [65]. The percent of urbanized land surrounding the traps was (mean ± S.D.): 55.75 ± 5.32 in Sanmin and 49.61 ± 8.49 in Nanzih, values not significantly different according to a Welch’s t test (t = 1.501, df = 8.408, p-value = 0.17). Moreover, landscape configuration, assessed through a Principal Components Analysis (PCA) [65,66] of the percent of land used for housing, roads, parks, and other purposes, described in Supplementary Figure S1 [online only], also showed no significant statistical differences between the two sites (t = 2.314, df = 7.195, p-value = 0.05289) when compared with a Welch’s t test, a t test that corrects df considering variance differences between the groups compared [67]. Despite the similarities in land cover, the studied districts have different population density per square KM, with Sanmin (17,321 habitants/KM2) more than doubling the population density of Nanzih (7130 habitants/KM2).
Only one weather station (Kaohsiung Weather Service Office; ID:46744 (120.30° E, 22.57° N)) (Figure 1) exists within the limits of KC, so all the climate variables used in this analysis, including daily average temperature, daily average precipitation, and daily average relative humidity, were summarized from the hourly measurements of the weather station. The study period coincided with one of the strongest El Niño events on record [66,67], and as consequence of this phenomenon, weather patterns were anomalous in KC. During the last quarter of 2016, precipitation was higher than usual in KC (Supplementary Figure S2A [Online Only]). Meanwhile, mean temperature (Supplementary Figure S2B [Online Only]) was hotter than usual during the whole study period, a pattern also observed for maximum (Supplementary Figure S2C [Online Only]) and minimum temperature (Supplementary Figure S2D [Online Only]). For the analysis, we employed daily data from 1 January 2016 to 31 December 2017, which was processed at biweekly intervals coinciding with mosquito sampling. For each time interval, we computed average values and the standard deviation and kurtosis of the above mentioned weather variables to understand the impacts of environmental variability on mosquito population dynamics [40]. The average distance between the weather station and mosquito traps was 12.13 KM (SD: 3.96 KM, Range: 8.1–16.9 KM). The resulting biweekly weather data can be observed in Figure 2. In Kaohsiung City, winters, from December to February, have low rainfall (Figure 2A), low relative humidity (Figure 2B), and temperature also drops to a minimum (Figure 2C). By contrast, during the summer, i.e., June to August, rainfall (Figure 2A), relative humidity (Figure 2B), and temperature (Figure 2C) reach their maximum values. With the exception of rainfall, which had a strong variation and which was followed in magnitude by relative humidity, the SD of all temperature variables was similar and low during the study period (Figure 2D). Meanwhile, the most leptokurtic variable, i.e., with high kurtosis values which represent low variation around the mean [68], where rainfall and maximum temperature, which reached maximum values from November to March, while the most platykurtic variable, more uncertain around the mean value than the extremes of its distribution [69], was minimum temperature (Figure 2E).

2.2. Mosquito Sampling

In this study, we collected mosquitoes with a total of 12 BG-Sentinel traps (Biogents, A.G. Regensburg, Germany). Traps were set in the two small area districts of Kaohsiung City: Sanmin (19.79 KM2) and Nanzih (25.83 KM2), as shown in Figure 1. Six traps were placed in each district, from 3 March 2016 to 26 December 2017, a number constrained by our limited budget and manpower but appropriate for the small area districts that we studied. Traps were baited with BG lure, a substance that mimics compounds released by the human skin [70]. Baits were replaced every two months following manufacturer recommendations [63,71]. Traps were operated continuously and mosquito samples were collected by laboratory staff on a bi-weekly basis. All traps were set in sheltered open household spaces, e.g., yards and/or garages, and household owners kindly provided electricity for trap operation. BG traps baited with BG lure were primarily chosen given their ability to capture Aedes (Stegomyia) spp., for example, Ae. aegypti and Aedes albopictus (Skuse) [62,72], but BG traps are known to capture a diverse set of mosquito species, including Culex (Culex) spp. [62,73,74] with medical importance such as Cx. quinquefasciatus in subtropical environments [63,73]. At the landscape level, all traps were located in households with similar environmental characteristics, in urban areas dominated by residential use. Mosquito samples were collected by the laboratory staff and mosquitoes were identified according to morphological characteristics presented in the taxonomic key by JC Lien [75]. Both Aedes spp. and Culex spp. were included in the analysis. Only female mosquitoes were considered in the study due to their importance for disease transmission.

2.3. Statistical Analysis

We analyzed average time series for Sanmin (dengue transmission hotspot) and Nanzih (dengue transmission coldspot) districts of KC. We analyzed average time series because a couple of the traps broke during the first year of the study, which we could not replace because of limited funding. When estimating the biweekly mean abundance, we considered all traps available each week, since there was no significant difference on the estimated mean abundance for each mosquito species and study location when these traps were excluded (Appendix A), thus ensuring mean estimates were not biased by including data from the traps that broke at the start of the study. For the time series analysis, we employed autoregressive (AR) models [76]. On seven occasions, it was not possible to collect samples on the exact date scheduled for sampling. Therefore, before the time series analysis, we used the lowess method to interpolate observations for the missing dates [77] in order to fulfill the assumption of data sampled at equal time intervals of AR models [76]. To fit the AR time series models, we followed a protocol described in detail by T Hoshi, Y Higa and LF Chaves [78], where null AR models are fit to mosquito abundance time series after selecting lags by inspecting the autocorrelation function (ACF) and then confirmed by inspecting the partial autocorrelation function (PACF) of the time series. Then, the time series of environmental covariates are pre-whitened using the coefficients of the null AR models. The resulting time series are used to estimate cross correlation functions (CCF), which are employed to select environmental covariates and their lags. The selected variables are then included in a full model for each mosquito time series, which is then simplified by a process of backward elimination based on the minimization of the Akaike Information Criterion (AIC). Prior to fitting each time series model, the mean of each environmental variable is removed in order to ease the interpretation of the time series model intercept as the mean value of the time series [23]. For further details about the methods, please refer to T Hoshi, Y Higa and LF Chaves [78] and LF Chaves, N Imanishi, and T Hoshi [23].

3. Results

We collected a total of 41,027 mosquitoes, over a total 746 biweekly trap collections. We collected a total of eight Aedes albopictus (Skuse), 2729 Ae. Aegypti, and 38,290 Cx. quinquefasciatus. In each one of the districts, we collected 4 Ae. albopictus during the full study period. In contrast, the average (±S.D.) abundance per trap and biweek of Ae. aegypti (9.03 ± 1.46) and Cx. quinquefasciatus (142.57 ± 14.38) were higher in Sanmin (Figure 3A,C) than Nanzih (Ae. aegypti: 6.21 ± 0.47 see Figure 3B; Cx. quinquefasciatus: 63.37 ± 8.71 see Figure 3D). Figure 3 also shows that starting in 26 August 2016 we lost two of the traps, reducing the number of traps to five at each site. Overall, the population dynamics of Ae. aegypti at Sanmin (Figure 3A) and Nanzih (Figure 3B) did not show major fluctuations, on average, through the study period. Nevertheless, some locations showed “outbreaks” in the sense, that often, sudden extraordinary changes in the abundance of Ae. aegypti were observed, which were larger in Sanmin (Figure 3A) when compared to Nanzih (Figure 3B).
In general, all time series were autocorrelated (Figure 4), and patterns of autocorrelation were mainly suggestive of first order AR processes, i.e., significant autocorrelation at one time lag, for Ae. aegypti in Sanmin (Figure 4A) and Nanzih (Figure 4B) and Cx. quinquefasciatus in Sanmin (Figure 4C) but not in Nanzih (Figure 4D). Using this information, first order autoregressive models were fitted as null models to all time series but Cx. quinquefasciatus at Nanzih (Figure 4D).
The PACF of the Ae. aegypti time series from Sanmin furtherly suggested a first order AR process (Figure 5A). This time series was not associated with the mean (Figure 5B) and SD (Figure 5C) of any weather variable, but it was associated positively at lag 1 and negatively at lag 2 with minimum temperature kurtosis (Figure 5D). By contrast, the PACF of Ae. aegypti in Nanzih did not show any significant autoregression (Figure 5E) and like the pattern observed in Sanmin the Nanzih time series was not correlated the mean (Figure 5F) and SD (Figure 5G) of any weather variable. However, it was associated positively at lag 1 with the minimum, mean and maximum temperature kurtosis (Figure 5H). Meanwhile, the PACF of the Cx. quinquefasciatus time series from Sanmin confirmed a first order AR process, with a significant correlation at lag 7 that was ignored when fitting the time series models (Figure 5I). This time series was significantly associated at lag 1 with the mean (Figure 5J) and SD (Figure 5K) of temperature (mean, maximum and minimum) and relative humidity and rainfall kurtosis (Figure 5L). The PACF of Cx. quinquefasciatus at Nanzih, by contrast, showed no significant correlation (Figure 5M) and was associated with the mean (Figure 5N) and SD (Figure 5O) of rainfall with five biweeks of lag, with temperature (minimum, mean and maximum) at lag 0 (Figure 5N), and one biweek of lag with rainfall kurtosis (Figure 5P).
Using the information from Figure 5, we then fitted full models, including all significantly correlated lags of environmental time series, for each mosquito time series. These full models were simplified by backward elimination based on the minimization of the AIC, and each process of model selection is presented as Supplementary online materials: Ae. aegypti at Sanmin (Supplementary Table S2 [online only]) and Nanzih (Supplementary Table S3 [online only]); Cx. quinquefasciatus at Sanmin (Supplementary Table S4 [online only]) and Nanzih (Supplementary Table S5 [online only]). Parameter estimates for the best time series models for Ae. aegypti are presented in Table 1.
The biweekly average number (Intercept in Table 1) of Ae. aegypti mosquitoes was 50% higher in Sanmin (around 9) when compared with Nanzih (around 6) and in both locations the dynamics were sensitive to changes in minimum and mean temperature kurtosis, which were the more platykurtic environmental variables at KC during the study period (Figure 2E). Parameter estimates for the best time series models for Cx. quinquefasciatus are presented in Table 2. Cx. quinquefasciatus average abundance at Sanmin (around 143) more than doubled what we found at Nanzih (around 63). In both districts increasing temperature had a negative impact on Cx. quinquefasciatus abundance, with a larger coefficient in Sanmin (Table 2), but only at Nanzih did the SD of Rainfall had a positive impact (Table 2) on Cx. quinquefasciatus abundance.

4. Discussion

Our data shows that, during the study period (March 2016–December 2017), a dengue transmission hotspot consistently had more mosquitoes than a neighboring transmission coldspot at KC, Taiwan. This fact supports the hypothesis that differences in transmission between these two districts might be related to differences in vector abundance. Moreover, the fact that both Ae. aegypti and Cx. quinquefasciatus were more abundant in the hotspot than the coldspot highlights the need to focus vector surveillance and control in hotspot areas. Specific to dengue transmission ecology in KC, our data showed the presence of Ae. aegypti near households, an observation furtherly supporting that Ae. aegypti might play the most important role in dengue virus transmission in southern Taiwan [79,80].
Our longitudinal mosquito data also highlights the importance of high order environmental variability, measured by kurtosis, on forcing mosquito abundance fluctuations, a phenomenon observed elsewhere for Ae. aegypti [32,39], Cx. quinquefasciatus [81,82,83] and other mosquito species [23,38,71,84], with patterns following the prediction of Schmalhausen’s law, the biological principle stating that organisms are sensitive not only to average environmental conditions but also to the environmental variability per se [25]. In our study, both Ae. aegypti and Cx. quinquefasciatus were sensitive to temperature, the most platykurtic environmental variable during the study period, yet Ae. aegypti fluctuations were associated with changes in temperature kurtosis, while Cx. quinquefasciatus tracked changes in maximum temperature. These differences might be related to different impacts of temperature on the life history of these two mosquito species. For example, Ae. aegypti abundance changes might reflect an accelerated larval development to avoid death by extreme hot temperatures, and/or by larval habitat desiccation [32,39,85]. Meanwhile, abundance changes in Cx. quinquefasciatus could be related to decreased oviposition [86,87] and reduced aquatic survival at extreme high temperatures [88]. The different feeding behaviors of Ae. aegypti (day time biter) and Cx. quinquefasciatus (night time biter) might also be associated with the environmental conditions at their active time [10,89].
Although our longitudinal study was long enough to show the importance of platykurtic environmental factors on mosquito abundance fluctuations, it is an open question if impacts of the El Niño Southern Oscillation (ENSO) on dengue transmission [20,90,91,92] could be related to changes in vector abundance. For example, in leishmaniasis we have shown that vector abundance and transmission increase during the cold ENSO phase in Panamá [93,94]. One limitation of this study is that given the lack of baseline longitudinal studies it is impossible to assess if ENSO impacted the population dynamics of both Ae. aegypti and Cx. quinquefasciatus, especially considering this study was done during one of the strongest ENSO events on record [66,67], which altered weather patterns at KC (Supplementary Figure S2 [Online Only]), and also altered weather and phenological patterns of mosquitoes elsewhere in east Asia, as documented for Tripteroides bambusa (Yamada) in Nagasaki, Japan, where the comparison of records from the 1970’s and 1980’s with data from 2016 and 2017 showed an increase in winter temperatures and that pupal recruitment season has lengthened its duration [95].
One potential limitation of this study was the use of weather data from KC’s weather station, which was around 12 km apart from the study sites. Our previous experience has been that over spatial scales similar to the one of this study [68,95], and even larger [71] temperature fluctuations are more synchronous than mosquito abundance, and very similar to observations from weather stations located at a distance similar to the one between KC’s weather station and our study locations [95]. Based on that, we believe data from KC’s weather station was appropriate to evaluate the association between weather fluctuations and mosquito abundance at the two study sites. Nevertheless, it is an open question whether more finely grained weather data, from weather sensors attached to the traps, could help to improve the understanding of weather impacts on vector abundance at finer spatial scales.
Finally, our findings, in conjunction with observations about the high correlation between entomological indicators and dengue transmission [54], support the need to strengthen systematic vector surveillance at sites that consistently behave as dengue transmission hotspots and the importance of mosquito abundance for transmission, since we significantly found more mosquitoes in Sanmin (the dengue transmission hotspot) than in Nanzih (the dengue transmission coldspot), though these patterns might also be affected by other environmental, demographic, and social aspects important for vector-borne disease transmission [24,25]. In that sense, the development of intelligent mosquito traps that automatically count and identify mosquitoes is a very promising tool that could increase the reliability of early warning systems for vector-borne disease transmission, which primarily rely on weather data [96,97] given delays associated with the processing of entomological information [49,98].

5. Conclusions

Our data suggest that spatially heterogeneous Ae. aegypti abundance might be a driver of spatial differences in dengue transmission at KC. Our comparative approach, where we also observed a significant increment in the abundance of Cx. quinquefasciatus in the dengue transmission hotspot, suggests this area might be more likely to experience outbreaks of other vector borne diseases and should become a primary focus for vector ecology research, surveillance and control.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4450/9/3/98/s1, Figure S1: Principal Component Analysis (PCA) to assess landscape configuration, Figure S2: Seasonal boxplots of climatic variables at Kaohsiung City, Taiwan. Table S1: The characteristics of BG-Traps in the sampling area. Table S2: Model selection results for Aedes aegypti in Sanmin district. Table S3: Model selection results for Aedes aegypti in Nanzih district. Table S4: Model selection results for Culex quinquefasciatus in Sanmin district, Table S4: Model selection results for Culex quinquefasciatus in Nanzih district.

Author Contributions

Data curation, K.-C.N.; Formal analysis, L.F.C. and T.-W.C.; Investigation, K.-C.N. and K.-H.T.; Methodology, L.F.C.; Resources, K.-H.T.; Supervision, K.-H.T. and T.-W.C.; Visualization, T.-W.C.; Writing—original draft, K.-C.N., L.F.C. and T.-W.C.

Funding

This study was funded by the Taiwan Ministry of Science and Technology (MOST 104-2119-M-038-002). LFC was funded by INCIENSA.

Acknowledgments

We thank the Department of Health, Kaohsiung City Government staff and the citizens who assisted with data collection.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Results of the Welch’s t test comparing the mean number of Culex quinquefasciatus and Aedes aegypti from biweeks 1 to 11 with and without the broken traps.
Cx. quinquefasciatus from Sanmin
Welch’s t = −0.7647, d.f. = 19.40, p-value = 0.4536
Cx. quinquefasciatus from Nanzih
Welch’s t = 0.1992, d.f. = 19.64, p-value = 0.8441
Ae. aegypti from Sanmin
Welch’s t = −0.1060, d.f. = 19.47, p-value = 0.9167
Ae. aegypti from Nanzih
Welch’s t = 0.4463, d.f. = 19.58, p-value = 0.6603

References

  1. Silver, J.B. Mosquito Ecology: Field Sampling Methods, 3rd ed.; Springer: New York, NY, USA, 2008. [Google Scholar]
  2. Burke, R.; Barrera, R.; Lewis, M.; Kluchinsky, T.; Claborn, D. Septic tanks as larval habitats for the mosquitoes Aedes albopictus and Culex quinquefasciatus in Playa-Playita, Puerto Rico. Med. Vet. Èntomol. 2010, 24, 117–123. [Google Scholar] [CrossRef] [PubMed]
  3. Reisen, W.K. The contrasting bionomics of Culex mosquitoes in western North America. J. Am. Mosq. Control. Assoc. 2012, 28, 82–91. [Google Scholar] [CrossRef] [PubMed]
  4. Russo, F.B.; Jungmann, P.; Beltrao-Braga, P.C.B. Zika infection and the development of neurological defects. Cell. Microbiol. 2017, 19, e12744. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Weaver, S.C.; Lecuit, M. Chikungunya virus and the global spread of a mosquito-borne disease. N. Engl. J. Med. 2015, 372, 1231–1239. [Google Scholar] [CrossRef] [PubMed]
  6. Simmons, C.P.; Farrar, J.J.; van Nguyen, V.; Wills, B. Dengue. N. Engl. J. Med. 2012, 366, 1423–1432. [Google Scholar] [CrossRef] [PubMed]
  7. Sasa, M. Mosquito as the vector of filariasis. Jpn. J. Sanit. Zoöl. 1965, 16, 171. [Google Scholar]
  8. Richards, S.L.; Lord, C.C.; Pesko, K.N.; Tabachnick, W.J. Environmental and Biological Factors Influencing Culex pipiens quinquefasciatus (Diptera: Culicidae) Vector Competence for West Nile Virus. Am. J. Trop. Med. Hyg. 2010, 83, 126–134. [Google Scholar] [CrossRef] [PubMed]
  9. Richards, S.L.; Lord, C.C.; Pesko, K.; Tabachnick, W.J. Environmental and Biological Factors Influencing Culex pipiens quinquefasciatus Say (Diptera: Culicidae) Vector Competence for Saint Louis Encephalitis Virus. Am. J. Trop. Med. Hyg. 2009, 81, 264–272. [Google Scholar] [PubMed]
  10. Edman, J.D. Disease control through manipulation of vector-host interaction: Some historical and evolutionary perspectives. In Proceedings of a Symposium: The Role of Vector-Host Interactions in Disease Tranmission; Scott, T.W., Grumstrup-Scott, J., Eds.; Entomological Society of America: Washington, DC, USA, 1988; pp. 43–50. [Google Scholar]
  11. Reisen, W.K.; Thiemann, T.; Barker, C.M.; Lu, H.L.; Carroll, B.; Fang, Y.; Lothrop, H.D. Effects of Warm Winter Temperature on the Abundance and Gonotrophic Activity of Culex (Diptera: Culicidae) in California. J. Med. Èntomol. 2010, 47, 230–237. [Google Scholar] [CrossRef] [PubMed]
  12. Campbell, L.P.; Luther, C.; Moo-Llanes, D.; Ramsey, J.M.; Danis-Lozano, R.; Peterson, A.T. Climate change influences on global distributions of dengue and chikungunya virus vectors. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2015, 370. [Google Scholar] [CrossRef] [PubMed]
  13. Metzger, M.E.; Hardstone Yoshimizu, M.; Padgett, K.A.; Hu, R.; Kramer, V.L. Detection and Establishment of Aedes albopictus and Aedes albopictus (Diptera: Culicidae) Mosquitoes in California, 2011–2015. J. Med. Èntomol. 2017, 54, 533–543. [Google Scholar] [CrossRef] [PubMed]
  14. Johnson, T.L.; Haque, U.; Monaghan, A.J.; Eisen, L.; Hahn, M.B.; Hayden, M.H.; Savage, H.M.; McAllister, J.; Mutebi, J.-P.; Eisen, R.J. Modeling the Environmental Suitability for Aedes (Stegomyia) aegypti and Aedes (Stegomyia) albopictus (Diptera: Culicidae) in the Contiguous United States. J. Med. Èntomol. 2017, 54, 1605–1614. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Ebi, K.L.; Nealon, J. Dengue in a changing climate. Environ. Res. 2016, 151, 115–123. [Google Scholar] [CrossRef] [PubMed]
  16. Messina, J.P.; Brady, O.J.; Scott, T.W.; Zou, C.; Pigott, D.M.; Duda, K.A.; Bhatt, S.; Katzelnick, L.; Howes, R.E.; Battle, K.E.; et al. Global spread of dengue virus types: Mapping the 70 year history. Trends Microbiol. 2014, 22, 138–146. [Google Scholar] [CrossRef] [PubMed]
  17. Smith, D.L.; Perkins, T.A.; Reiner, J.R.C.; Barker, C.M.; Niu, T.; Chaves, L.F.; Ellis, A.M.; George, D.B.; Le Menach, A.; Pulliam, J.R.C.; et al. Recasting the theory of mosquito-borne pathogen transmission dynamics and control. Trans. R. Soc. Trop. Med. Hyg. 2014, 108, 185–197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Wallace, R.; Chaves, L.F.; Bergmann, L.; Ayres Lopes, C.F.J.; Hogerwerf, L.; Kock, R.; Wallace, R.G. Clear-Cutting Disease Control: Capital-Led Deforestation, Public Health Austerity, and Vector-Borne Infection; Springer: New York, NY, USA, 2018. [Google Scholar]
  19. Weaver, S.C.; Reisen, W.K. Present and future arboviral threats. Antivir. Res. 2010, 85, 328–345. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Chuang, T.W.; Chaves, L.F.; Chen, P.J. Effects of local and regional climatic fluctuations on dengue outbreaks in southern Taiwan. PLoS ONE 2017, 12, e0178698. [Google Scholar] [CrossRef] [PubMed]
  21. Bolling, B.G.; Barker, C.M.; Moore, C.G.; Pape, W.J.; Eisen, L. Seasonal patterns for entomological measures of risk for exposure to Culex vectors and West Nile virus in relation to human disease cases in northeastern Colorado. J. Med. Entomol. 2009, 46, 1519–1531. [Google Scholar] [CrossRef] [PubMed]
  22. Garrett-Jones, C. Prognosis for Interruption of Malaria Transmission through Assessment of Mosquitos Vectorial Capacity. Nature 1964, 204, 1173–1175. [Google Scholar] [CrossRef] [PubMed]
  23. Chaves, L.F.; Imanishi, N.; Hoshi, T. Population dynamics of Armigeres subalbatus (Diptera: Culicidae) across a temperate altitudinal gradient. Bull. Èntomol. Res. 2015, 105, 589–597. [Google Scholar] [CrossRef] [PubMed]
  24. Chaves, L.F. Climate change and the biology of insect vectors of human pathogens. In Invertebrates and Global Climate Change; Johnson, S., Jones, H., Eds.; Wiley: Chichester, UK, 2017; pp. 126–147. [Google Scholar]
  25. Chaves, L.F.; Koenraadt, C.J.M. Climate Change and Highland Malaria: Fresh Air for a Hot Debate. Q. Rev. Boil. 2010, 85, 27–55. [Google Scholar] [CrossRef] [Green Version]
  26. Bar-Zeev, M. The effect of temperature on the growth rate and survival of the immature stages of Aedes albopictus (L.). Bull. Èntomol. Res. 1958, 49, 157–163. [Google Scholar] [CrossRef]
  27. Couret, J.; Benedict, M.Q. A meta-analysis of the factors influencing development rate variation in Aedes albopictus (Diptera: Culicidae). BMC Ecol. 2014, 14, 3. [Google Scholar] [CrossRef] [PubMed]
  28. Couret, J. Meta-Analysis of Factors Affecting Ontogenetic Development Rate in the Culex pipiens (Diptera: Culicidae) Complex. Environ. Èntomol. 2013, 42, 614–626. [Google Scholar] [CrossRef] [PubMed]
  29. Rueda, L.M.; Patel, K.J.; Axtell, R.C.; Stinner, R.E. Temperature-Dependent Development and Survival Rates of Culex quinquefasciatus and Aedes albopictus (Diptera: Culicidae). J. Med. Èntomol. 1990, 27, 892–898. [Google Scholar] [CrossRef] [PubMed]
  30. Watts, D.M.; Burke, D.S.; Harrison, B.A.; Whitmire, R.E.; Nisalak, A. Effect of temperature on the vector efficiency of Aedes albopictus for dengue 2 virus. Am. J. Trop. Med. Hyg. 1986, 36, 143–152. [Google Scholar] [CrossRef]
  31. Reisen, W.K.; Fang, Y.; Martinez, V.M. Effects of Temperature on the Transmission of West Nile Virus by Culex tarsalis (Diptera: Culicidae). J. Med. Èntomol. 2006, 43, 309–317. [Google Scholar] [CrossRef] [PubMed]
  32. Chaves, L.F.; Scott, T.W.; Morrison, A.C.; Takada, T. Hot temperatures can force delayed mosquito outbreaks via sequential changes in Aedes albopictus demographic parameters in autocorrelated environments. Acta Trop. 2014, 129, 15–24. [Google Scholar] [CrossRef] [PubMed]
  33. Berryman, A.; Barbosa, P.; Schultz, J. The theory and classification of outbreaks. In Insect Outbreaks; Barbosa, P., Schultz, J.C., Eds.; Academic Press: San Diego, CA, USA, 1987; pp. 3–30. [Google Scholar]
  34. Ye, Y.; Louis, V.R.; Simboro, S.; Sauerborn, R. Effect of meteorological factors on clinical malaria risk among children: An assessment using village-based meteorological stations and community-based parasitological survey. BMC Public Health 2007, 7, 101. [Google Scholar] [CrossRef] [PubMed]
  35. Day, J.F.; Curtis, G.A. When it rains they soar—And that makes Culex nigripalpus a dangerous mosquito. Am. Èntomol. 1994, 40, 162–167. [Google Scholar] [CrossRef]
  36. Dow, R.P.; Gerrish, G.M. Day-to-day change in relative humidity and activity of Culex nigripalpus (Diptera: Culicidae). Ann. Èntomol. Soc. Am. 1970, 63, 995–999. [Google Scholar] [CrossRef] [PubMed]
  37. Edman, J.D.; Scott, T.W.; Costero, A.; Morrison, A.C.; Harrington, L.C.; Clark, G.G. Aedes albopictus (Diptera: Culicidae) movement influenced by availability of oviposition sites. J. Med. Èntomol. 1998, 35, 578–583. [Google Scholar] [CrossRef] [PubMed]
  38. Chaves, L.F. Globally invasive, withdrawing at home: Aedes albopictus and Aedes japonicus facing the rise of Aedes Flavopictus. Int. J. Biometeorol. 2016, 60, 1727–1738. [Google Scholar] [CrossRef] [PubMed]
  39. Chaves, L.F.; Morrison, A.C.; Kitron, U.D.; Scott, T.W. Nonlinear impacts of climatic variability on the density-dependent regulation of an insect vector of disease. Glob. Chang. Boil. 2012, 18, 457–468. [Google Scholar] [CrossRef]
  40. Chaves, L.F.; Moji, K. Density Dependence, Landscape, and Weather Impacts on Aquatic Aedes japonicus japonicus (Diptera: Culicidae) Abundance along an Urban Altitudinal Gradient. J. Med. Èntomol. 2018, 55, 329–341. [Google Scholar] [CrossRef] [PubMed]
  41. Koenraadt, C.J.M.; Harrington, L.C. Flushing Effect of Rain on Container-Inhabiting Mosquitoes Aedes albopictus and Culex pipiens (Diptera: Culicidae). J. Med. Èntomol. 2008, 45, 28–35. [Google Scholar] [CrossRef] [PubMed]
  42. Baak-Baak, C.M.; Moo-Llanes, D.A.; Cigarroa-Toledo, N.; Puerto, F.I.; Machain-Williams, C.; Reyes-Solis, G.; Nakazawa, Y.J.; Ulloa-Garcia, A.; Garcia-Rejon, J.E. Ecological Niche Model for Predicting Distribution of Disease-Vector Mosquitoes in Yucatan State, Mexico. J. Med. Entomol. 2017, 54, 854–861. [Google Scholar] [CrossRef] [PubMed]
  43. Trewin, B.J.; Kay, B.H.; Darbro, J.M.; Hurst, T.P. Increased container-breeding mosquito risk owing to drought-induced changes in water harvesting and storage in Brisbane, Australia. Int. Health 2013, 5, 251–258. [Google Scholar] [CrossRef] [PubMed]
  44. Barrera, R.; Avila, J.; Gonzalez-Tellez, S. Unreliable supply of potable water and elevated Aedes albopictus larval indices—A casual relationship. J. Am. Mosq. Control. Assoc. 1993, 9, 189–195. [Google Scholar] [PubMed]
  45. Hoshi, T.; Imanishi, N.; Higa, Y.; Chaves, L.F. Mosquito Biodiversity Patterns around Urban Environments in South-Central Okinawa Island, Japan. J. Am. Mosq. Control. Assoc. 2014, 30, 260–267. [Google Scholar] [CrossRef] [PubMed]
  46. Torres, R.; Samudio, R.; Carrera, J.-P.; Young, J.; Márquez, R.; Hurtado, L.; Weaver, S.; Chaves, L.F.; Tesh, R.; Cáceres, L. Enzootic mosquito vector species at equine encephalitis transmission foci in the República de Panamá. PLoS ONE 2017, 12, e0185491. [Google Scholar] [CrossRef] [PubMed]
  47. Chaverri, L.G.; Dillenbeck, C.; Lewis, D.; Rivera, C.; Romero, L.M.; Chaves, L.F. Mosquito Species (Diptera: Culicidae) Diversity from Ovitraps in a Mesoamerican Tropical Rainforest. J. Med. Èntomol. 2018, 55, 646–653. [Google Scholar] [CrossRef] [PubMed]
  48. Reisen, W.K.; Lothrop, H.D. Effects of sampling design on the estimation of adult mosquito abundance. J. Am. Mosq. Control. Assoc. 1999, 15, 105–114. [Google Scholar] [PubMed]
  49. Shand, L.; Brown, W.M.; Chaves, L.F.; Goldberg, T.L.; Hamer, G.L.; Haramis, L.; Kitron, U.; Walker, E.D.; Ruiz, M.O. Predicting West Nile Virus Infection Risk from the Synergistic Effects of Rainfall and Temperature. J. Med. Èntomol. 2016, 53, 935–944. [Google Scholar] [CrossRef] [PubMed]
  50. Ruiz, M.O.; Chaves, L.F.; Hamer, G.L.; Sun, T.; Brown, W.M.; Walker, E.D.; Haramis, L.; Goldberg, T.L.; Kitron, U.D. Local impact of temperature and precipitation on West Nile virus infection in Culex species mosquitoes in northeast Illinois, USA. Parasites Vectors 2010, 3, 19. [Google Scholar] [CrossRef] [PubMed]
  51. Scott, T.W.; Morrison, A.C.; Lorenz, L.H.; Clark, G.G.; Strickman, D.; Kittayapong, P.; Zhou, H.; Edman, J.D. Longitudinal studies of Aedes albopictus (Diptera: Culicidae) in Thailand and Puerto Rico: Population dynamics. J. Med. Èntomol. 2000, 37, 77–88. [Google Scholar] [CrossRef] [PubMed]
  52. Reisen, W.K.; Meyer, R.P.; Shields, J.; Arbolante, C. Population ecology of preimaginal Culex tarsalis (Diptera: Culicidae) in Kern County, California. J. Med. Entomol. 1989, 26, 10–22. [Google Scholar] [CrossRef] [PubMed]
  53. Reisen, W.K.; Milby, M.M.; Meyer, R.P. Population dynamics of adult Culex mosquitoes (Diptera: Culicidae) along the Kern River, Kern County, California, in 1990. J. Med. Entomol. 1992, 29, 531–543. [Google Scholar] [CrossRef] [PubMed]
  54. Cromwell, E.A.; Stoddard, S.T.; Barker, C.M.; Van Rie, A.; Messer, W.B.; Meshnick, S.R.; Morrison, A.C.; Scott, T.W. The relationship between entomological indicators of Aedes albopictus abundance and dengue virus infection. PLoS Negl. Trop. Dis. 2017, 11, e0005429. [Google Scholar] [CrossRef] [PubMed]
  55. Chaves, L.F.; Hamer, G.L.; Walker, E.D.; Brown, W.M.; Ruiz, M.O.; Kitron, U.D. Climatic variability and landscape heterogeneity impact urban mosquito diversity and vector abundance and infection. Ecosphere 2011, 2, 1–21. [Google Scholar] [CrossRef]
  56. Reisen, W.K. Landscape epidemiology of vector-borne diseases. Annu. Rev. Entomol. 2010, 55, 461–483. [Google Scholar] [CrossRef] [PubMed]
  57. Wang, S.F.; Chang, K.; Loh, E.W.; Wang, W.H.; Tseng, S.P.; Lu, P.L.; Chen, Y.H.; Chen, Y.A. Consecutive large dengue outbreaks in Taiwan in 2014–2015. Emerg. Microbes Infect. 2016, 5, e123. [Google Scholar] [CrossRef] [PubMed]
  58. Wang, S.F.; Wang, W.H.; Chang, K.; Chen, Y.H.; Tseng, S.P.; Yen, C.H.; Wu, D.C.; Chen, Y.M. Severe Dengue Fever Outbreak in Taiwan. Am. J. Trop. Med. Hyg. 2016, 94, 193–197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Chang, K.; Chen, C.D.; Shih, C.M.; Lee, T.C.; Wu, M.T.; Wu, D.C.; Chen, Y.H.; Hung, C.H.; Wu, M.C.; Huang, C.C.; et al. Time-Lagging Interplay Effect and Excess Risk of Meteorological/Mosquito Parameters and Petrochemical Gas Explosion on Dengue Incidence. Sci. Rep. 2016, 6, 35028. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Li, Y.; Su, X.; Zhou, G.; Zhang, H.; Puthiyakunnon, S.; Shuai, S.; Cai, S.; Gu, J.; Zhou, X.; Yan, G.; et al. Comparative evaluation of the efficiency of the BG-Sentinel trap, CDC light trap and Mosquito-oviposition trap for the surveillance of vector mosquitoes. Parasit Vectors 2016, 9, 446. [Google Scholar] [CrossRef] [PubMed]
  61. Wright, J.A.; Larson, R.T.; Richardson, A.G.; Cote, N.M.; Stoops, C.A.; Clark, M.; Obenauer, P.J. Comparison of BG-Sentinel(R) Trap and Oviposition Cups for Aedes albopictus and Aedes albopictus Surveillance in Jacksonville, Florida, USA. J. Am. Mosq. Control. Assoc. 2015, 31, 26–31. [Google Scholar] [CrossRef] [PubMed]
  62. Degener, C.M.; Azara, T.M.; Roque, R.A.; Codeco, C.T.; Nobre, A.A.; Ohly, J.J.; Geier, M.; Eiras, A.E. Temporal abundance of Aedes albopictus in Manaus, Brazil, measured by two trap types for adult mosquitoes. Mem. Inst. Oswaldo Cruz 2014, 109, 1030–1040. [Google Scholar] [CrossRef] [PubMed]
  63. Maciel-de-Freitas, R.; Eiras, Á.E.; Lourenço-de-Oliveira, R. Field evaluation of effectiveness of the BG-Sentinel, a new trap for capturing adult Aedes albopictus (Diptera: Culicidae). Mem. Inst. Oswaldo Cruz 2006, 101, 321–325. [Google Scholar] [CrossRef] [PubMed]
  64. Harrington, L.C.; Scott, T.W.; Lerdthusnee, K.; Coleman, R.C.; Costero, A.; Clark, G.G.; Jones, J.J.; Kitthawee, S.; Kittayapong, P.; Sithiprasasna, R. Dispersal of the dengue vector Aedes albopictus within and between rural communities. Am. J. Trop. Med. Hyg. 2005, 72, 209–220. [Google Scholar] [PubMed]
  65. Medeiros, M.C.I.; Boothe, E.C.; Roark, E.B.; Hamer, G.L. Dispersal of male and female Culex quinquefasciatus and Aedes albopictus mosquitoes using stable isotope enrichment. PLoS Negl. Trop. Dis. 2017, 11, e0005347. [Google Scholar] [CrossRef] [PubMed]
  66. Levine, A.F.; McPhaden, M.J. How the July 2014 easterly wind burst gave the 2015–2016 El Niño a head start. Geophys. Res. Lett. 2016, 43, 6503–6510. [Google Scholar] [CrossRef]
  67. McPhaden, M.J. Playing hide and seek with El Nino. Nat. Clim Chang. 2015, 5, 791–795. [Google Scholar] [CrossRef]
  68. Chaves, L.F. Mosquito Species (Diptera: Culicidae) Persistence and Synchrony across an Urban Altitudinal Gradient. J. Med. Èntomol. 2017, 54, 329–339. [Google Scholar] [CrossRef] [PubMed]
  69. Ross, S.M. Introduction to Probability Models; Academic Press: Cambridge, MA, USA, 2014. [Google Scholar]
  70. Irish, S.R.; Chandre, F.; N’Guessan, R. Comparison of Octenol- And BG Lure®-Baited Biogents Sentinel Traps and an Encephalitis Virus Surveillance Trap in Portland, OR. J. Am. Mosq. Control. Assoc. 2008, 24, 393–397. [Google Scholar] [CrossRef] [PubMed]
  71. Chaves, L.F.; Higa, Y.; Lee, S.H.; Jeong, J.Y.; Heo, S.T.; Kim, M.; Minakawa, N.; Lee, K.H. Environmental Forcing Shapes Regional House Mosquito Synchrony in a Warming Temperate Island. Environ. Èntomol. 2013, 42, 605–613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Ritchie, S.A.; Moore, P.; Carruthers, M.; Williams, C.; Montgomery, B.; Foley, P.; Ahboo, S.; Van Den Hurk, A.F.; Lindsay, M.D.; Cooper, B.; et al. Discovery of a Widespread Infestation of Aedes albopictus in the Torres Strait, Australia. J. Am. Mosq. Control. Assoc. 2006, 22, 358–365. [Google Scholar] [CrossRef]
  73. Kröckel, U.; Rose, A.; Eiras, Á.E.; Geier, M. New tools for surveillance of adult yellow fever mosquitoes: Comparison of trap catches with human landing rates in an urban environment. J. Am. Mosq. Control. Assoc. 2006, 22, 229–238. [Google Scholar] [CrossRef]
  74. Farajollahi, A.; Kesavaraju, B.; Price, D.C.; Williams, G.M.; Healy, S.P.; Gaugler, R.; Nelder, M.P. Field efficacy of BG-Sentinel and industry-standard traps for Aedes albopictus (Diptera: Culicidae) and West Nile virus surveillance. J. Med. Èntomol. 2009, 46, 919–925. [Google Scholar] [CrossRef] [PubMed]
  75. Lien, J.C. Pictorial Keys to Mosquitoes of Taiwan; Yi Hsien Publishing Co., Ltd.: New Taipei, Taiwan, 2004. [Google Scholar]
  76. Shumway, R.H.; Stoffer, D.S. Time Series Analysis and Its Applications, 3rd ed.; Springer: New York, NY, USA, 2011. [Google Scholar]
  77. Cleveland, W.S.; Devlin, S.J. Locally Weighted Regression: An Approach to Regression Analysis by Local Fitting. J. Am. Stat. Assoc. 1988, 83, 596–610. [Google Scholar] [CrossRef]
  78. Hoshi, T.; Higa, Y.; Chaves, L.F. Uranotaenia novobscura ryukyuana (Diptera: Culicidae) Population Dynamics Are Denso-Dependent and Autonomous from Weather Fluctuations. Ann. Èntomol. Soc. Am. 2014, 107, 136–142. [Google Scholar] [CrossRef]
  79. Tsai, P.J.; Teng, H.J. Role of Aedes albopictus (Linnaeus) and Aedes albopictus (Skuse) in local dengue epidemics in Taiwan. BMC Infect. Dis. 2016, 16, 662. [Google Scholar] [CrossRef] [PubMed]
  80. Yang, C.F.; Hou, J.N.; Chen, T.H.; Chen, W.J. Discriminable roles of Aedes albopictus and Aedes albopictus in establishment of dengue outbreaks in Taiwan. Acta Trop. 2014, 130, 17–23. [Google Scholar] [CrossRef] [PubMed]
  81. Hayes, J. Seasonal changes in population structure of Culex pipiens quinquefasciatus Say (Diptera: Culicidae): Study of an isolated population. J. Med. Èntomol. 1975, 12, 167–178. [Google Scholar] [CrossRef] [PubMed]
  82. Hayes, J.; Hsi, B.P. Interrelationships between selected meteorologic phenomena and immature stages of Culex pipiens quinquefasciatus Say: Study of an isolated population. J. Med. Èntomol. 1975, 12, 299–308. [Google Scholar] [CrossRef] [PubMed]
  83. Hayes, J.; Downs, T.D. Seasonal changes in an isolated population of Culex pipiens quinquefasciatus (Diptera: Culicidae): A time series analysis. J. Med. Èntomol. 1980, 17, 63–69. [Google Scholar] [CrossRef]
  84. Shaman, J.; Day, J.F. Reproductive Phase Locking of Mosquito Populations in Response to Rainfall Frequency. PLoS ONE 2007, 2, e331. [Google Scholar] [CrossRef] [PubMed]
  85. Focks, D.A.; Haile, D.G.; Daniels, E.; Mount, G.A. Dynamic Life Table Model for Aedes albopictus (Diptera, Culcidae)—Analysis of the Literature and Model Development. J. Med. Èntomol. 1993, 30, 1003–1017. [Google Scholar] [CrossRef] [PubMed]
  86. Nguyen, A.T.; Williams-Newkirk, A.J.; Kitron, U.D.; Chaves, L.F. Seasonal weather, nutrient dynamics and conspecific presence impacts on the southern house mosquito oviposition dynamics in combined sewage overflows. J. Med. Èntomol. 2012, 49, 1328–1338. [Google Scholar] [CrossRef] [PubMed]
  87. Chaves, L.F.; Kitron, U.D. Weather variability impacts on oviposition dynamics of the southern house mosquito at intermediate time scales. Bull. Èntomol. Res. 2011, 101, 633–641. [Google Scholar] [CrossRef] [PubMed]
  88. Chaves, L.F.; Keogh, C.L.; Nguyen, A.M.; Decker, G.M.; Vazquez-Prokopec, G.M.; Kitron, U.D. Combined sewage overflow accelerates immature development and increases body size in the urban mosquito Culex Quinquefasciatus. J. Appl. Èntomol. 2011, 135, 611–620. [Google Scholar] [CrossRef]
  89. Scott, T.W.; Amerasinghe, P.H.; Morrison, A.C.; Lorenz, L.H.; Clark, G.G.; Strickman, D.; Kittayapong, P.; Edman, J.D. Longitudinal studies of Aedes albopictus (Diptera: Culicidae) in Thailand and Puerto Rico: Blood feeding frequency. J. Med. Èntomol. 2000, 37, 89–101. [Google Scholar] [CrossRef] [PubMed]
  90. Cazelles, B.; Chavez, M.; McMichael, A.J.; Hales, S. Nonstationary influence of El Nino on the synchronous dengue epidemics in Thailand. PLoS Med. 2005, 2, e106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Johansson, M.A.; Cummings, D.A.; Glass, G.E. Multiyear climate variability and dengue—El Nino southern oscillation, weather, and dengue incidence in Puerto Rico, Mexico, and Thailand: A longitudinal data analysis. PLoS Med. 2009, 6, e1000168. [Google Scholar] [CrossRef] [PubMed]
  92. Thai, K.T.D.; Cazelles, B.; Nguyen, N.V.; Vo, L.T.; Boni, M.F.; Farrar, J.; Simmons, C.P.; van Doorn, H.R.; de Vries, P.J. Dengue Dynamics in Binh Thuan Province, Southern Vietnam: Periodicity, Synchronicity and Climate Variability. PLoS Negl. Trop. Dis. 2010, 4, e747. [Google Scholar] [CrossRef] [PubMed]
  93. Yamada, K.; Valderrama, A.; Gottdenker, N.; Cerezo, L.; Minakawa, N.; Saldaña, A.; Calzada, J.E.; Chaves, L.F. Macroecological patterns of American Cutaneous Leishmaniasis transmission across the health areas of Panamá (1980–2012). Parasite Epidemiol. Control. 2016, 1, 42–55. [Google Scholar] [CrossRef] [PubMed]
  94. Chaves, L.F.; Calzada, J.E.; Valderama, A.; Saldaña, A. Cutaneous Leishmaniasis and Sand Fly fluctuations are associated with El Niño in Panamá. PLoS Negl. Trop. Dis. 2014, 8, e3210. [Google Scholar] [CrossRef] [PubMed]
  95. Chaves, L.F.; Jian, J.-Y.; Moji, K. Overwintering in the Bamboo Mosquito Tripteroides bambusa (Diptera: Culicidae) During a Warm, But Unpredictably Changing, Winter. Environ. Èntomol. 2018, 47, 148–158. [Google Scholar] [CrossRef] [PubMed]
  96. Chaves, L.F.; Pascual, M. Comparing Models for Early Warning Systems of Neglected Tropical Diseases. PLoS Negl. Trop. Dis. 2007, 1, e33. [Google Scholar] [CrossRef] [PubMed]
  97. Shi, Y.; Liu, X.; Kok, S.-Y.; Rajarethinam, J.; Liang, S.; Yap, G.; Chong, C.-S.; Lee, K.-S.; Tan, S.S.; Chin, C.K.Y. Three-month real-time dengue forecast models: An early warning system for outbreak alerts and policy decision support in Singapore. Environ. Health Perspect. 2016, 124, 1369–1375. [Google Scholar] [CrossRef] [PubMed]
  98. Hii, Y.L.; Rocklöv, J.; Wall, S.; Ng, L.C.; Tang, C.S.; Ng, N. Optimal Lead Time for Dengue Forecast. PLoS Negl. Trop. Dis. 2012, 6, e1848. [Google Scholar] [CrossRef] [PubMed]
Insects 09 00098 g001 550
Figure 1. Study Locations. (A) Sanmin—dengue transmission hotspot (B) Nanzih—dengue transmission coldspot. The small map shows the location of Kaohsiung City, southern Taiwan, and also Sanmin and Nanzih districts.
Figure 1. Study Locations. (A) Sanmin—dengue transmission hotspot (B) Nanzih—dengue transmission coldspot. The small map shows the location of Kaohsiung City, southern Taiwan, and also Sanmin and Nanzih districts.
Insects 09 00098 g001
Insects 09 00098 g002 550
Figure 2. Weather Time series data at Kaohsiung, Taiwan. (A) Rainfall in mm; (B) Relative humidity, %; (C) Temperature including Maximum (Max), Minimum (Min), and Average (Mean), see inset legend for details about the lines. (D) Standard deviation, SD, of the environmental variables, lines are color coded according to panels A, B, and C. (E) Kurtosis of the environmental variables, lines are color coded according to panels A, B, and C, for rainfall kurtosis black dots indicate infinite values, a perfect leptokurtic distribution estimated during periods without rainfall.
Figure 2. Weather Time series data at Kaohsiung, Taiwan. (A) Rainfall in mm; (B) Relative humidity, %; (C) Temperature including Maximum (Max), Minimum (Min), and Average (Mean), see inset legend for details about the lines. (D) Standard deviation, SD, of the environmental variables, lines are color coded according to panels A, B, and C. (E) Kurtosis of the environmental variables, lines are color coded according to panels A, B, and C, for rainfall kurtosis black dots indicate infinite values, a perfect leptokurtic distribution estimated during periods without rainfall.
Insects 09 00098 g002
Insects 09 00098 g003 550
Figure 3. Mosquito Abundance Time Series. Aedes aegypti per BG Sentinel trap collection in (A) Sanmin and (B) Nanzih. Culex quinquefasciatus per BG Sentinel trap collection in (C) Sanmin and (D) Nanzih. At each site line colors correspond to the trap locations in Figure 1, with the exception of black lines which represent average values, see inset legend in (A). Grey circles indicate abundance estimates interpolated with lowess, see inset legend in (B). To ease visualization the y axis in (C,D) is shown in logarithmic scale.
Figure 3. Mosquito Abundance Time Series. Aedes aegypti per BG Sentinel trap collection in (A) Sanmin and (B) Nanzih. Culex quinquefasciatus per BG Sentinel trap collection in (C) Sanmin and (D) Nanzih. At each site line colors correspond to the trap locations in Figure 1, with the exception of black lines which represent average values, see inset legend in (A). Grey circles indicate abundance estimates interpolated with lowess, see inset legend in (B). To ease visualization the y axis in (C,D) is shown in logarithmic scale.
Insects 09 00098 g003
Insects 09 00098 g004 550
Figure 4. Autocorrelation functions (ACF) Aedes aegypti per BG Sentinel trap collection in (A) Sanmin and (B) Nanzih. Culex quinquefasciatus per BG Sentinel trap collection in (C) Sanmin and (D) Nanzih. In the panels, dashed blue lines are 95% confidence limits within which correlations are expected to occur by chance.
Figure 4. Autocorrelation functions (ACF) Aedes aegypti per BG Sentinel trap collection in (A) Sanmin and (B) Nanzih. Culex quinquefasciatus per BG Sentinel trap collection in (C) Sanmin and (D) Nanzih. In the panels, dashed blue lines are 95% confidence limits within which correlations are expected to occur by chance.
Insects 09 00098 g004
Insects 09 00098 g005 550
Figure 5. Partial autocorrelation (PACF) and cross correlation functions (CCF). Aedes aegypti in Sanmin: (A) PACF and (B) CCF with mean values of environmental variables (C) CCF with SD values of environmental variables (D) CCF with kurtosis values of environmental variables. Ae. aegypti in Nanzih: (E) PACF and (F) CCF with mean values of environmental variables (G) CCF with SD values of environmental variables (H) CCF with kurtosis values of environmental variables. Culex quinquefasciatus in Sanmin: (I) PACF and (J) CCF with mean values of environmental variables (K) CCF with SD values of environmental variables (L) CCF with kurtosis values of environmental variables. Cx. quinquefascuatus in Nanzih: (M) PACF and (N) CCF with mean values of environmental variables (O) CCF with SD values of environmental variables (P) CCF with kurtosis values of environmental variables. In all panels dashed blue lines are 95% confidence limits within which correlations are expected to occur by chance. In the CCF plots lines, environmental variables are color coded following the pattern described in Figure 3.
Figure 5. Partial autocorrelation (PACF) and cross correlation functions (CCF). Aedes aegypti in Sanmin: (A) PACF and (B) CCF with mean values of environmental variables (C) CCF with SD values of environmental variables (D) CCF with kurtosis values of environmental variables. Ae. aegypti in Nanzih: (E) PACF and (F) CCF with mean values of environmental variables (G) CCF with SD values of environmental variables (H) CCF with kurtosis values of environmental variables. Culex quinquefasciatus in Sanmin: (I) PACF and (J) CCF with mean values of environmental variables (K) CCF with SD values of environmental variables (L) CCF with kurtosis values of environmental variables. Cx. quinquefascuatus in Nanzih: (M) PACF and (N) CCF with mean values of environmental variables (O) CCF with SD values of environmental variables (P) CCF with kurtosis values of environmental variables. In all panels dashed blue lines are 95% confidence limits within which correlations are expected to occur by chance. In the CCF plots lines, environmental variables are color coded following the pattern described in Figure 3.
Insects 09 00098 g005
Table
Table 1. Parameter estimates for the best time series explaining abundance changes in Aedes aegypti sampled with BG-sentinel traps at Sanmin (dengue transmission hotsop) and Nanzih (dengue transmission coldspot), Kaohsiung City, Taiwan. In parameters, Intercept is the mean value of the time series and AR = autoregressive term, K = kurtosis, T = temperature, mean = average and min = minimum. Lags are in biweeks.
Table 1. Parameter estimates for the best time series explaining abundance changes in Aedes aegypti sampled with BG-sentinel traps at Sanmin (dengue transmission hotsop) and Nanzih (dengue transmission coldspot), Kaohsiung City, Taiwan. In parameters, Intercept is the mean value of the time series and AR = autoregressive term, K = kurtosis, T = temperature, mean = average and min = minimum. Lags are in biweeks.
Parameters (lag)SanminNanzih
EstimatesS.E.zp-ValueEstimatesS.E.zp-Value
Intercept9.03371.45716.196.65 × 10−10 *6.21490.471613.181.18 × 10−39 *
AR (1)0.55370.12874.301.69 × 10−5 *----
KTmin (1)2.00870.89282.250.024 *----
KTmin (2)−2.08860.8976−2.330.019 *----
KTmean (1) ----2.65120.5914.497.26 × 10−6 *
Error Variance18.77---9.33---
* Statistically significant (p < 0.05).
Table
Table 2. Parameter estimates for the best time series explaining abundance changes in Culex quinquefasciatus sampled with BG-sentinel traps at Sanmin (dengue transmission hotsop) and Nanzih (dengue transmission coldspot), Kaohsiung City, Taiwan. In parameters Intercept is the mean value of the time series and SD = Standard deviation, T = temperature, max = maximum, R = Rainfall. Lags are in biweeks.
Table 2. Parameter estimates for the best time series explaining abundance changes in Culex quinquefasciatus sampled with BG-sentinel traps at Sanmin (dengue transmission hotsop) and Nanzih (dengue transmission coldspot), Kaohsiung City, Taiwan. In parameters Intercept is the mean value of the time series and SD = Standard deviation, T = temperature, max = maximum, R = Rainfall. Lags are in biweeks.
Parameter (lag)SanminNanzih
EstimatesS.E.zp-ValueEstimatesS.E.zp-Value
Intercept142.573814.37829.923.55 × 10−23 *63.37028.717.283.45 × 10−13 *
Tmax (1)−22.09495.1298−4.311.65 × 10−5 *----
SDR (5)----1.73110.48843.540.00039 *
Tmax (0)----−11.85352.9783−3.986.89 × 10−5 *
Error Variance7743---2882---
* Statistically significant (p < 0.05).

Share and Cite

MDPI and ACS Style

Ng, K.-C.; Chaves, L.F.; Tsai, K.-H.; Chuang, T.-W. Increased Adult Aedes aegypti and Culex quinquefasciatus (Diptera: Culicidae) Abundance in a Dengue Transmission Hotspot, Compared to a Coldspot, within Kaohsiung City, Taiwan. Insects 2018, 9, 98. https://doi.org/10.3390/insects9030098

AMA Style

Ng K-C, Chaves LF, Tsai K-H, Chuang T-W. Increased Adult Aedes aegypti and Culex quinquefasciatus (Diptera: Culicidae) Abundance in a Dengue Transmission Hotspot, Compared to a Coldspot, within Kaohsiung City, Taiwan. Insects. 2018; 9(3):98. https://doi.org/10.3390/insects9030098

Chicago/Turabian Style

Ng, Ka-Chon, Luis Fernando Chaves, Kun-Hsien Tsai, and Ting-Wu Chuang. 2018. "Increased Adult Aedes aegypti and Culex quinquefasciatus (Diptera: Culicidae) Abundance in a Dengue Transmission Hotspot, Compared to a Coldspot, within Kaohsiung City, Taiwan" Insects 9, no. 3: 98. https://doi.org/10.3390/insects9030098

APA Style

Ng, K. -C., Chaves, L. F., Tsai, K. -H., & Chuang, T. -W. (2018). Increased Adult Aedes aegypti and Culex quinquefasciatus (Diptera: Culicidae) Abundance in a Dengue Transmission Hotspot, Compared to a Coldspot, within Kaohsiung City, Taiwan. Insects, 9(3), 98. https://doi.org/10.3390/insects9030098

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop