Abstract

Dengue is an important vector-borne disease transmitted by the mosquitoes Aedes aegypti and Ae. albopictus. In the absence of an effective vaccine, vector control has become the key intervention tool in controlling the disease. Vector densities are significantly affected by the changing weather patterns of a region. The present study was conducted in three selected localities, i.e., urban Bandaranayakapura, semiurban Galgamuwa, and rural Buluwala in the Kurunegala district of Sri Lanka to assess spatial and temporal distribution of dengue vector mosquitoes and to predict vector prevalence with respect to changing weather parameters. Monthly ovitrap surveys and larval surveys were conducted from January to December 2019 and continued further in the urban area up to December 2021. Aedes aegypti was found moderately in the urban area and to a lesser extent in semiurban but not in the rural area. Aedes albopictus had the preference for rural over urban areas. Aedes aegypti preferred indoor breeding, while Ae. albopictus preferred both indoor and outdoor. For Ae. albopictus, ovitrap index (OVI), premise index (PI), container index (CI), and Breteau index (BI) correlated with both the rainfall (RF) and relative humidity (RH) of the urban site. Correlations were stronger between OVI and RH and also between BI and RF. Linear regression analysis was fitted, and a prediction model was developed using BI and RF with no lag period (R2 (sq) = 86.3%; F = 53.12; R2 (pred) = 63.12%; model: Log10 (BI) = 0.153 + 0.286 Log10 (RF); RMSE = 1.49). Another prediction model was developed using OVI and RH with one month lag period (R2 (sq) = 70.21%; F = 57.23; model: OVI predicted = 15.1 + 0.528 Lag 1 month RH; RMSE = 2.01). These two models can be used to monitor the population dynamics of Ae. albopictus in urban settings to predict possible dengue outbreaks.

1. Introduction

Dengue is a fast-growing mosquito borne disease, found especially in the countries located within the equatorial zone, and has been ranked as one of the top ten threats to global health in 2019 by the World Health Organization [1, 2]. It is an arboviral disease caused by infection of one or more of the four dengue virus serotypes. Dengue virus is transmitted from human to human through the bites of infective females of Aedes aegypti Linnaeus and Ae. albopictus Skuse mosquitoes (Diptera: Culicidae) [3, 4]. Aedes aegypti is mainly endophilic and considered the primary vector of the disease. It has adapted to human habitats and breeds primarily in artificial water containers such as water storage barrels, old tires, and flower pots. Aedes albopictus is mainly exophilic, less adapted to human habitats, and mainly breeds in natural containers such as tree stumps and coconut shells and to a lesser extent in artificial containers [57]. Although Ae. aegypti has a higher distribution range and a higher disease transmission capacity, Ae. albopictus is considered one of the most invasive mosquito species emerging as a global public-health threat [8].

In Sri Lanka, the prevalence of high-risk level dengue is mainly due to favorable meteorological conditions which promote high vector abundance. Previous workers have shown unique distribution patterns of the two vector species in different geographical regions of the country. Aedes aegypti is the predominant Aedes species in northern [9] and western [10] provinces, whereas Ae. albopictus is dominant in western, north-western, central, and southern provinces [1113]. Since an effective vaccine against the virus is not available, dengue prevention and control mainly depend on effective vector control measures against immature aquatic stages and adults. Many studies have shown that the creation of new larval habitats, as a result of rapid urbanization programs, has resulted in the expansion of dengue vector populations into new areas playing a significant role in vector distribution patterns [6, 14, 15]. Therefore, urbanization substantially increases adult mosquito density, larval development rate, and adult survival time of the vector species which, in turn, potentially increase disease transmission.

In Sri Lanka, dengue has been reported as a public health problem since the 1960s. A major upward shift in the disease incidence has been reported since 2009. Sri Lanka faced a massive dengue epidemic in 2017, with 186,101 reported cases and 440 deaths [16]. Kurunegala district, situated in the north-western province of Sri Lanka, annually reports >2000 dengue cases, and dengue fever outbreaks have frequently occurred in the district since 2014 (Epidemiological data, 2020). As in several other districts of the country, rapid urbanization processes have been taken place in the district in the recent past. However, the district reports high case incidence rates of dengue not only from urban areas but also from semiurban and rural areas. Out of 29 Medical Officer of Health (MOH) areas in Kurunegala district, 2 have been classified as priority high risk MOH areas and 12 as high risk for dengue. Vector control in the district is mainly by source reduction, application of insecticides, public health education, and legislations. The present study aimed to assess spatial and temporal distribution of dengue vector mosquitoes in some selected localities of the Kurunegala district of Sri Lanka over a three-year period and to predict vector prevalence with respect to changing weather parameters.

2. Methods

2.1. Study Sites

The study was conducted in three selected localities in the Kurunegala district (7°45′ N, 80° 15′ E) which is located in the North-Western Province of Sri Lanka covering 4812.7 km2 with a population of approximately 1,676,000 living in nearly 439,065 households. About 32.6% of the population depends on agriculture-related employment. Average maximum and minimum temperatures were approximately 32.83°C and 23.43°C, respectively, during the study period. Annual accumulated rainfall and average humidity were approximately 164.65 mm and 67%. The three selected localities i.e., urban Bandaranayakapura (BAN) in Kurunegala MOH area, semiurban Galgamuwa (GAL) in Galgamuwa MOH area, and rural Buluwala (BUL) in Rideegama MOH area (Figure 1), were from dengue high risk areas based on the dengue and dengue epidemic prevalence data for the previous five-year period (2013–2018). The study was carried out only in residential premises (premises with care of house owners) of all study sites. Bandaranayakapura (96 km2) is an urbanized locality with about 600 households with water supply, sanitation facilities, electricity services, and rubbish collection within an area of. Galgamuwa (278.4 km2) is a suburban locality with about 500 households, poor sanitary facilities, and no municipal water supply. Buluwala (220 km2) is a rural area with 400 householders, dense vegetation, and mountains.

2.2. Egg Surveys

Monthly ovitrap surveys were conducted from January to December 2019 in 50 randomly selected house premises at each study site except for urban BAN area where the time period of the survey was extended up to December 2021. Following WHO guidelines, mosquito eggs were collected using standard ovitraps which were made with black plastic cups (8 cm height × 7 cm diameter) filled with dechlorinated tap water.

A rectangular filter paper strip (3 × 22 cm) was placed inside the cup as the oviposition substrate. Two traps were kept at each house premise and the filter paper strips were labeled with the date, house number, and the location. One ovitrap was kept outdoor 5–10 m away from the house (depending on the size of the garden) and was labeled as “A.” The other trap “B” was kept inside the house. Filter paper strip of each trap was replaced with a new one once in every 5 days, and the removed filter papers were brought to the laboratory at the entomological surveillance unit, RDHS office, Kurunegala. In the laboratory, filter paper strips with the eggs were air-dried and stored at room temperature. They were examined under a stereomicroscope and the eggs of the two dengue vector species were identified using the shape and the size. To confirm the identity, randomly selected eggs of the two types were reared in separate containers and the larvae were identified using the key introduced by Tanaka and Mizusawa [17]. Number of eggs of the two dengue vector species were counted and recorded separately. Ovitrap index (OVI) was calculated as the percentage of positive ovitraps, and the egg density index (EDI) was calculated as the average number of eggs per trap for each survey.

2.3. Larval Surveys

Monthly larval surveys were conducted from January to December 2019 at each study sites except for BAN area where the survey was continued up to December 2021. Following the systematic sampling approach, every 6th premise in the inspection route of the locality was chosen for sampling and a minimum of 100 premises was selected from each study site. At each premise, all the potential containers in both indoor and outdoor areas were examined, and sampling was performed from all positive containers. Larvae were collected using dippers and pipets (50 ml) depending on the nature of the container. Collected larvae were transferred to labeled plastic vials and transported to the laboratory for further analysis. All the larvae were reared to 3rd and 4th instar levels for identification using standard larval identification keys [17]. Larval indices, i.e., Premise/house Index (PI/HI), Container Index (CI), and Breteau Index (BI) were calculated for both vector species to determine the larval abundance at each site [2, 18].

2.4. Meteorological Data

Maximum and minimum daily temperature, daily rainfall, and daily relative humidity for all study sites were obtained from the Department of Meteorology, Colombo, Sri Lanka.

2.5. Statistical Analyses

Data obtained were formatted for layouts using MS-Excel. As a preliminary analysis, meteorological data were subjected to the Anderson–Darling Normality Test at a significance level of 5% using Minitab 16. When distributions were not in accordance with the criteria of normality, data were transformed into Log10 to obtain a normal distribution. Study results were descriptively summarized using frequencies, proportions, means, and standard deviations. The larval indices (PI, CI, and BI) and egg indices (OVI and EDI) of urban, semiurban, and rural settings were comparatively analyzed by mean comparisons using the independent sample t-test and one-way ANOVA.

2.5.1. Correlation and Regression Analyses

Correlation analysis was performed using Spearman’s rank correlation analysis. Relationship between the indices and the weather parameters (monthly rainfall (mm), minimum and maximum temperature (°C), and relative humidity (%)) in different sites was analyzed. Spatial autocorrelations were tested using Moran’s I statistic, and the variables selected for the final model satisfied the independency at 5% significant level. Correlation analysis was done between the selected meteorological variable with different time lags (no lag, 1-week lag, 2-week lag etc.). The lag period that had the highest correlation with the index was used to create the models.

Stepwise multiple regression analysis (SMR) was done to develop a prediction model to evaluate the association between weather variables with the larval indices (PI, CI and BI) and egg indices (OVI and EDI). A forecasting regional model (FRM) was then developed and evaluated to predict the larval indices. The predicted larval indices could then be used to forecast possible dengue outbreaks by predicting the patterns of larval abundance of dengue vectors. The general model was based on the following equation:where y is the estimated index value, is a constant, are coefficients of environmental variables .

R-squared and F-test were used to statistically measure how close the data are to the fitted regression line. It is also known as the coefficient of determination, or the coefficient of multiple determinations for multiple regressions. To validate models, actual and predicted values were used.

3. Results

3.1. Abundance of Dengue Vector Species in Different Urban Settings
3.1.1. Egg Survey

A total of 24,403 eggs were collected from the ovitrap survey. Out of 3550 ovitraps placed, 1535 were positive for Aedes eggs in the three study sites during January to December 2019 (urban = 421 (27.46%); semiurban = 513 (33.42%); rural = 601 (39.15%)). Aedes aegypti eggs were less abundant (n = 497 (2.03%)) while Ae. albopictus eggs were dominant (n = 23,872 (97.82%)) in positive ovitraps. Mixed infestation of both Ae. aegypti and Ae. albopictus was not common (n = 34, 0.14%).

Aedes aegypti eggs were present only in urban and semiurban areas, and OVI at urban site was significantly higher than that at semiurban site () (Table 1). Aedes albopictus eggs were predominant in all study sites. Although the highest number of eggs were received from rural followed by semiurban sites, the OVI of Ae. albopictus was not significantly different between study sites (). Aedes aegypti eggs showed a significantly higher abundance in indoor than outdoor collections (), and Ae. albopictus eggs were more abundant in outdoor collections () (Table 1).

3.1.2. Larval Survey

During the period January to December 2019, a total of 3693 house premises were inspected. Out of them, 240 (6.4%) premises were positive for dengue vector larvae. A total of 2172 Aedes larvae were collected from these premises [Ae. aegypti larvae n = 94 (4.33%); Ae. albopictus larvae n = 2078 (95.67%)]. Aedes aegypti larvae were obtained only from the urban locality (n = 94, 4.33%) while Ae. albopictus larvae were reported from all three study sites [Urban: 531 (24.4%); semi-urban: 567 (26.1%); rural: 981 (45.16%)]. However, the larval indices PI, CI and BI for Ae. albopictus were not significantly different among the study sites () (Table 2).

3.2. Seasonal Fluctuation of OVI and BI

Larval numbers of BAN Ae. aegypti were too low to show the seasonal fluctuation of BI. Monthly variations of OVI for BAN Ae. aegypti and, OVI and BI for Ae. albopictus from all three study sites were plotted against the rain fall and relative humidity for the period January to December 2019 (Figure 2). An increasing trend could be observed both in OVI and BI with the rainfall in all three areas. However, a decline of BI was observed with heavy rain falls in September-November, probably due to flush over of the breeding sites (Figure 2).

3.3. Correlation of OVI and BI with Weather Parameters

Sessional variation of egg and larval indices of both species were analyzed against rain fall (RF), relative humidity (RH), and temperature data for the period January to December 2019. For Ae. aegypti, no correlation was obtained between the indices and weather parameters since the egg and larval numbers were very low. All the variables indicated Moran’s I index was close to zero. Container Index and Premise Index of Ae. albopictus showed a significant spatial autocorrelation across the study period, while BI indicated a nonsignificant () and weak spatial autocorrelation (Moran’s I = 0.125). We found BI to be suitable for developing the prediction model. Ovitrap Index of Ae. albopictus showed lower spatial autocorrelation (Moran’s I = 0.002) across the study period. Therefore, development of a prediction model for mosquito abundance using ovitrap collections was suggested. For Ae. albopictus, OVI, PI, CI, and BI correlated with both the RF and RH only for the urban site. However, correlations were strongly positive between OVI and RH (r = 0.680, ), and between BI and RF (r = .970, ). Therefore, it was decided to continue the study in urban (BAN) area for further two years to develop prediction models for these indices.

3.4. Prediction Model for OVI

When considering three-year data from Jan 2019 to Dec 2021, the correlation between OVI and RH was still significant and was strongly positive (r = 0.786, ). The association between OVI and RH was then evaluated in 1 week, 2 week, 3 week, 4 week, and 8 week lag periods. Significant positive correlations were obtained for 4 week (r = 0.610, ) and 8 week (r = 0.624, ) lag periods.

Linear regression analysis was fitted and a prediction model was developed using OVI and RH with no lag period (R2 (sq) = 77.24%; F = 23.7; model: OVI predicted = −28.9 + 0.931 RH; RMSE = 3.47), with 1 month lag period (R2 (sq) = 70.21%; F = 57.23; model: OVI predicted = 15.1 + 0.528 Lag 1 month RH; RMSE = 2.01), with 2-month lag period (OVI predicted = 14.7 + 0.536 Lag 2 month RH; RMSE = 2.28).

Since the lowest root mean square error (RSME) value (2.01) was with actual and lag 1-month predicted values, a model was generated for RH/actual OVI and RH/1-month lag OVI for urban (BAN) Ae. albopictus (Figure 3). It is interesting to see that the linear predicted 1-month lag OVI line directly overlaps with the linear actual OVI line.

3.5. Prediction Model for BI

Three-year data (2019–2021) showed a positive correlation between BI and RF (r = 0.789, ) for BAN Ae. albopictus. The association between BI and RF was evaluated with actual time RF and lag periods after RF. Significant correlations were obtained up to the 3-week lag period (r = 0.970 for no lag, 0.495 for 1-week lag, 0.381 for 2-week lag, and 0.306 for 3-week lag).

Linear regression analysis was fitted, and a prediction model was developed using BI and RF with no lag period (R2 (sq) = 86.3%; F = 53.12; R2 (pred) = 63.12%; model: Log10 (BI) = 0.153 + 0.286∗ Log10 (RF); RMSE = 1.49) (Figure 4).

4. Discussion

Dengue is considered the most prevalent vector borne disease in tropics. Also, it is reportedly sensitive to environmental changes such as climate change and urbanization. Spatial distribution of vectors along urban-rural gradients, with Ae. aegypti preferring urban and Ae. albopictus preferring rural environments has been observed previously [19, 20], and the present results further confirm this observation. During the present study Ae. aegypti was found in the urban area, and to a lesser extent, in the semiurban area but not in the rural area. Aedes albopictus was dominating from rural to the urban area indicating its invasion of urban areas. Higher adaptability and invasive behavior of this species have been reported elsewhere [15, 21]. Our results confirm previous findings that Ae. aegypti prefers indoor breeding to outdoor while Ae. albopictus prefers both indoor and outdoor breeding [12, 22, 23]. Both the feeding and breeding behavior of Ae. aegypti are more or less limited to domestic environments whereas the activity pattern of Ae. albopictus has a wider range as it equally prefers to feed on both wild animals and humans and is more adaptive to breed in variety of environments [3, 19, 20, 24]. As far as the transmission of the disease is concerned, both these vectors are equally important, as DENV has been isolated from both these vectors [25, 26]. Some studies have suggested that the dengue virus transmission risk is higher with Ae. albopictus due to its wider distribution [27]. High prevalence of both these vectors has been reported during the 2017-2018 dengue outbreak in the Kurunegala district (Entomological report data, RDHS, Kurunegala).

Sri Lanka receives southwest monsoon rains during May to September and northeast monsoon rains during December to February. Surendran et al. [28] has reported egg density peaks corresponding to increased precipitation in monsoonal periods from the northern Jaffna district of Sri Lanka. Kurunegala district is situated in northwest region and hardly receives monsoon rains. The district receives rains mainly from intermonsoons, which occur throughout the country during March to April and October to November periods and a mild rainfall from conventional rains. All three study sites reported increased precipitation towards the end of the year. Studies have emphasized that an OVI above 10% for Aedes vector species is an indication of a possible risk of dengue outbreak, provided that dengue positive cases are present in the area [29]. If so, all three study sites used for the present study were at dengue outbreak-risk state if the virus is available. Entomological surveillances are mainly based on different larval indices and the House Index (HI-percentage of houses positive for larvae) and the Breteau index (BI-number of positive containers per 100 houses) are the most widely used larval indices [4]. However, discriminating thresholds have never been determined for larval indices for dengue fever transmission. Values of HI > 1% and BI > 5 have been proposed to indicate high risk levels for yellow fever transmission [4]. If a similar level of discrimination is applied for dengue, our study areas are at high risk of dengue if the virus is available.

Weather factors such as rainfall (RF), temperature, and relative humidity (RH) have been widely studied for their usability in predicting dengue incidence trends [30, 31]. Previous workers have reported that temperature and humidity have significant associations with the dengue incidence [3236]. It has been suggested that RH can affect larval density by extending adult mosquito survival [37]. Wu et al. [32] and Campbell et al. [34] identified RF as having a strong positive correlation with dengue incidence. There can be several factors, other than weather, which may act critically in determining spatial distribution of dengue mosquitoes. Urbanization, waste management, sanitary facilities, community awareness, and vector control strategies are some of these. However, some researchers have shown that the value of entomological indices as predictors of Aedes-borne disease outbreaks is inconclusive and limited due to the complex epidemiology of dengue, which involves dynamic interplay of multiple factors such as herd immunity within a population, distinct serotypes of the virus, intervention programs, mosquito-human interaction, and mosquito-virus interaction [38, 39]. Nevertheless, significant influence of meteorological factors on the development, survival, density, and oviposition rate of the vector mosquitoes, expressed as larval [40], pupal [41], and adult indices [38] can certainly make an indication on dengue disease transmission.

Forecasting vector abundance with weather variables have been attempted by several workers. Moore et al. [42] in Puerto Rico and Pontes et al. [40] in Fortaleza, Brazil, used temporal graphics to compare seasonal fluctuation of RF, Aedes larval indices, and dengue incidence. They observed a strong relationship between larval indices and RF patterns. Present study demonstrated that OVI, PI, CI, and BI of Ae. albopictus correlates with both the RF and RH in urban areas. Strong positive association between OVI and RH and between BI and RF for Ae. albopictus in urban areas compelled us to continue sampling for two more years and to develop prediction models to forecast OVI and BI by using respective RH and RF data. However, as it has been reported previously, heavy precipitation flushed the water containers removing mosquito larvae away from breeding sites resulting a decline of BI [43, 44].

A significant relationship, giving the lowest RSME value (2.01), was observed between RH and Ae. albopictus OVI with a lag of 1 month in urban area. In the generated model, the linear predicted 1-month lag OVI line directly overlaps with the linear actual OVI line showing the accuracy in predicting. Several workers have shown the influence of RH on adult mortality [45, 46], egg maturity process, and egg laying inducement in Aedes mosquitoes [47]. Prediction model for BI against RF gave a low RMSE (1.49) indicating that the model can be successfully used to forecast mosquito populations with respect to RF changes in urban areas. It was observed that BI significantly correlates with RF even after a period of up to three weeks, for Ae. Albopictus. Assessment of RMSE and R square values suggested that our model can provide accurate predictions on the abundance of Ae. Albopictus. Therefore, these models can be successfully utilized to predict vector prevalence in urban sites, and successful vector control interventions can be adopted during the lag periods to minimize disease incidence.

5. Conclusions

Aedes aegypti population was mainly found in urban area, and to a lesser extent, in semiurban area. It was not found in rural areas. Aedes albopictus was present in all three areas with a preference from rural to urban but dominating even the urban settings. Aedes aegypti preferred indoor breeding while Ae. albopictus preferred both indoor and outdoor. Strong correlations were observed between the ovitrap index (OVI) and the relative humidity (RH) and between Breteau index (BI) and the rainfall (RF). Prediction models were developed using BI and RF with no lag period and using OVI and RH with 1 month lag period. These two models will be able to detect population dynamics of Ae. albopictus vector mosquitoes in urban settings allowing to predict possible dengue outbreaks.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

Assistance given by the staff of the Entomological Surveillance Unit, Office of Regional Director of Health Services, Kurunegala, and the Department of Zoology, University of Peradeniya, Sri Lanka, is acknowledged. This work was supported by the Office of Regional Director of Health Services, Kurunegala, and the University of Peradeniya, Sri Lanka.