Skip to main content

Spatial and temporal patterns of dengue incidence in northeastern Thailand 2006–2016

Abstract

Background

Dengue, a viral disease transmitted by Aedes mosquitoes, is an important public health concern throughout Thailand. Climate variables are potential predictors of dengue transmission. Associations between climate variables and dengue have usually been performed on large-scale first-level national administrative divisions, i.e. provinces. Here we analyze data on a finer spatial resolution in one province, which is often more relevant for effective disease control design. The objective of this study was to investigate the effect of seasonal variations, monthly climate variability, and to identify local clusters of symptomatic disease at the sub-district level based on reported dengue cases.

Methods

Data on dengue cases were retrieved from the national communicable disease surveillance system in Thailand. Between 2006 and 2016, 15,167 cases were recorded in 199 sub-districts of Khon Kaen Province, northeastern Thailand. Descriptive analyses included demographic characteristics and temporal patterns of disease and climate variables. The association between monthly disease incidence and climate variations was analyzed at the sub-district level using Bayesian Poisson spatial regression. A hotspot analysis was used to assess the spatial patterns (clustered/dispersed/random) of dengue incidence.

Results

Dengue was predominant in the 5–14 year-old age group (51.1%). However, over time, dengue incidence in the older age groups (> 15 years) gradually increased and was the most affected group in 2013. Dengue outbreaks coincide with the rainy season. In the spatial regression model, maximum temperature was associated with higher incidence. The hotspot analysis showed clustering of cases around the urbanized area of Khon Kaen city and in rural areas in the southwestern portion of the province.

Conclusions

There was an increase in the number of reported dengue cases in older age groups over the study period. Dengue incidence was highly seasonal and positively associated with maximum ambient temperature. However, climatic variables did not explain all the spatial variation of dengue in the province. Further analyses are needed to clarify the detailed effects of urbanization and other potential environmental risk factors. These results provide useful information for ongoing prediction modeling and developing of dengue early warning systems to guide vector control operations.

Peer Review reports

Background

The annual global burden of dengue is estimated at 390 million infections, of which 96 million present clinically [1]. Four closely related RNA viruses in the family Flaviviridae (DENV1 to DENV4) are responsible for dengue disease. They are transmitted by Aedes (primarily subgenus Stegomyia) mosquitoes, particularly Aedes aegypti (L.) and Aedes albopictus (Skuse) [2]. Dengue has developed from a sporadically occurring disease to a major and re-emerging global public health problem over recent decades causing substantial economic disruption and social burden in endemic areas in Asia, Africa, and the Americas. There is no effective treatment for dengue and vaccination, so far, offers only incomplete protection [3, 4]. Therefore, vector control remains the most important means of prevention [5]. Effective vaccine or not, vector control will remain the cornerstone of dengue control for years to come [3].

Due to increasing incidence and rapid geographical expansion, dengue is the most common vector-borne disease in Thailand [6]. From 2000 to 2011, the number of reported cases varied from 20,000 to 140,000 cases each year [7]. Both Ae. aegypti and Ae. albopictus are common species and widely distributed in Thailand [8]. All four serotypes co-circulated in each of the major outbreaks that occurred in 1958, 1987, 1998, 2001, 2013, and 2015 [9,10,11,12,13,14]. The highest incidence typically occurs in 13–24 year-old age group with case clustering seen predominately in urban areas [15]. Males represent the majority of reported dengue cases in several Asian countries [16]. A study in Singapore showed that men were more exposed to infected mosquitoes than women, during daytime hours, at the workplace or while travelling to and from work. A forceful public health policy in Singapore [17] has greatly reduced the number of mosquitoes in and around homes, potentially rendering the larger male labor force more exposed to mosquito bites during working hours [16, 18]. Other causes for these apparent gender differences could be different health seeking behaviors or male-female differences in disease severity [19]. In the Lao People’s Democratic Republic male-female ratios in dengue cases varied between years and provinces [16]. We are not aware of similar spatio-temporal or socioeconomic differences in Thailand.

Thailand has adapted the dengue control strategy of the World Health Organization (WHO) [2], which consists of three main pillars: 1) patients diagnosed with dengue are required to avoid mosquito bites to prevent dengue transmission; 2) active community case detection of cases which do not result in clinical consultation; and 3) vector control, consisting of environmental management, source reduction, and chemical interventions using insecticide fogging against adult vectors and larvicides to control immature stages in containers [20]. Follow-up interventions are conducted by health officers or village health volunteers [20]. To determine the most appropriate and feasible intervention or combination of interventions, health officers need to consider local environmental, resource, and contextual factors that may influence effectiveness [21].

Climate variables are predictors of dengue infection [4, 22, 23]. Seasonal variation in climate shows a strong relationship with Ae. aegypti abundance and historical dengue incidence [24]. Temperature affects population biology of Aedes mosquitoes [25]. Higher temperatures increase larval development [26] and rates of multiple feeding, but reduce mosquito size [27]. The extrinsic incubation period declines as temperature rises, thus increasing the proportion of infected vectors, and enhancing the transmission potential of the vector [27,28,29].

As ambient temperature increases, so does dengue epidemic potential, peaking at around 29 °C and then decreases [29]. In subtropical and tropical regions such as Thailand, with mean diel temperatures of 26 °C (20 °C ≤ T ≤ 32 °C), an increase in diurnal temperature range can enhance transmission [29]. An analysis of data from Thailand (1978–1997) showed the incidence of dengue hemorrhagic fever (DHF) was negatively associated with higher rainfall in the southern region of the country, but positively associated with elevated ambient temperatures in the central and northern regions [30]. Another study using provincial monthly dengue data from 1983 to 2001 concluded that the relationships between weather variables and dengue transmission are very complex in Thailand [31]. The study found that transmission occurs within a specific temperature range, but that changes in humidity within this range can amplify the transmission potential with 80% of dengue cases occurring at a mean temperature of between 27.0 and 29.5 °C and a mean relative humidity of > 75%. They further found that large epidemics begin earlier, develop faster and can be predicted at a defined onset time. Non-linear modeling of more than 30 years (1982–2013) of monthly data by province in Thailand showed that inter-annual variations in rainfall and temperature with a lag time of one month can improve the explanation of dengue relative risk compared to a seasonal-spatial model [32]. The relationship between rainfall and dengue is complex, as it may create abundant breeding sites for the vector [33], but can also flush out sites if rain is too intense [33, 34]. Because household water storage may increase in the dry season, the resulting breeding habitats may weaken, or even reverse, the positive association between dengue and rainfall [35,36,37,38,39].

Spatio-temporal analysis can detect clusters of dengue disease and is useful for a better understanding of the dynamics of disease dispersion. Analysis of spatial and temporal variations is also useful in identifying high-risk locations and times of higher transmission risk, which are important for disease surveillance and control [15, 40].

The above-mentioned research on climate and dengue focused on larger spatiotemporal scales, such as monthly dengue surveillance and climate records at the provincial level [31, 41, 42]. The current study is novel because it uses data on the lowest administrative level, the sub-district, in one province to understand fine-scale spatial dengue-climate relationships. This is useful for developing more reliable prediction models for future projections applied in early warning and response systems, thus ultimately improving timely control interventions.

We analyzed data on reported dengue cases in Khon Kaen Province, northeastern Thailand collected between 2006 and 2016 to 1) describe demographic characteristics and seasonal variations of dengue cases; 2) determine the potential impact of climate variability on dengue incidence; and 3) identify clusters of dengue cases at the sub-district level.

Methods

Study area

The study was conducted in Khon Kaen Province, an area of approximately 10,900 km2 (16°25′12″N to 16°42′12″N and 102°49′48″E to 102°83′48″E). The province has 26 districts, 199 sub-districts, and 2139 villages. In 2010, the population was 1,767,601, of which 387,279 people lived in Mueang District that includes the provincial capital Khon Kaen (see Additional file 1). This province was selected as the study area because dengue is endemic with typical seasonal increases and occasional outbreaks. The province is primarily rural with a few large urban centers. Mueang District, the most densely populated area in the province, is a regional center for education, health, finance and commerce. The northern and southern parts of the district, along the major highway linking Bangkok with Lao People’s Democratic Republic, are rapidly developing. The districts in the northwestern and southeastern parts of the province are rural and agricultural. Classification of urban and rural areas depends on population density. An urban area is defined as a municipality or town with a population over 100,000 and a population density above 300 persons per square kilometer [43]. The average minimum and maximum seasonal temperatures are 16.7 °C (December–January) and 36.4 °C (April–May). The monthly minimum and maximum rainfall vary from 0 mm (dry season: November–April) to 240 mm (wet season: May–October).

Data collection

The Office of Disease Prevention and Control, Region 7 Khon Kaen (ODPC7), Department of Disease Control, Ministry of Public Health, Thailand provided data on the weekly number of reported dengue cases in Khon Kaen Province from 1 January 2006 to 31 December 2016. Dengue is a notifiable disease based on the National Communicable Disease Control Law, i.e., all government and private hospitals, clinics and other healthcare facilities must report all cases (confirmed and suspected) to the local health authority within 24 h of diagnosis [12]. Cases are recorded by degree of disease severity into one of three categories (at peak of illness): 1) dengue fever (DF), 2) dengue hemorrhagic fever (DHF), and 3) dengue shock syndrome (DSS), but the serotype is not recorded (or typically known except retrospectively). A patient is diagnosed with suspected DF when the following criteria are met and signs and symptoms are present: residence or recent travel to a dengue endemic area, acute fever accompanied by any two of the following: headache, myalgia, arthralgia, rash, positive tourniquet test and leucopenia, with no evidence of plasma leakage. DHF is recorded in patients with a temperature ≥ 38 °C, petechiae, ecchymosis, or a positive tourniquet test, thrombocytopenia (platelets < 100,000 cells/mm3), and evidence of plasma leakage. DSS, the most severe disease manifestation, is defined as having the same signs and symptoms as DHF, but progressing to circulatory failure. The Provincial Health Offices enter patient data into the standardized Disease Surveillance Report (Report 506) for recording communicable diseases in Thailand. The form provides the patient’s age, gender, house address, signs and symptoms, and date of medical consultation. DHF and DSS are based on both clinical symptoms and laboratory tests (usually complete blood count), and sometimes accompanied with a rapid diagnostic test (RDT); whereas, DF is seldom based on additional laboratory tests or by RDT.

Meteorological data from 1 January 2006 to 31 December 2016 were downloaded from the data library of the International Research Institute for Climate and Society [44], which contains specific climate data from different sources, such as The National Centers for Environmental Prediction (NCEP), Climate Forecast System Reanalysis (CFSR) [45], and Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS) global rainfall datasets [46]. For each sub-district, daily temperatures (°C) were retrieved from NCEP and daily rainfall (mm) from CHIRPS. These data generated the monthly means used in the analysis (see Additional files 2 and 3). The spatial resolution of rainfall is 0.05 × 0.05 degrees (CHIRPS) and for temperature 0.2 × 0.2 degrees (NCEP CFSR v2, https://rda.ucar.edu/datasets/ds094.1/). A centroid was created for each sub-district. Rainfall and temperature data for each sub-district was determined based on the grid cell in which the centroid was located.

Analysis

Monthly data on dengue cases and climate (rainfall and temperature) from the study period were combined to visualize seasonal patterns and temporal trends. Dengue incidence was calculated using the monthly number of reported cases and sub-district population size in 2010 reflecting the mid-study denominator [47].

Bayesian Poisson regression models were used to assess associations with the number of monthly cases in 199 sub-districts. Population was used as the denominator in the model (i.e. log-population as an offset). The neighborhood relationship between the sub-districts were defined using their adjacency matrix; ‘1’ for a pair of sub-districts sharing a border, otherwise ‘0’. Hence the following model was used:

$$ {Y}_{ij}\sim Poisson\ \left({\mu}_{ij}\right) $$
$$ \log\ \left({\mu}_{ij}\right)=\log\ \left({P}_i\right)+{\theta}_{ij} $$
$$ {\theta}_{ij}=\alpha +{\beta}_k{x}_{ij k}+{u}_{ij}, $$

where Yij is the observed mean number of cases for the ith sub-district in jth month (i = 1, …,199; j = 1, …,12), Pi is the sub-district population size, α is the intercept, and βk is the regression coefficient for covariate k. For the main model, the covariates (xk) were: population density per square kilometer; gender (proportion of males among the cases), mean age in years of the cases; mean rainfall; and minimum and maximum temperature. As a non-mechanistic way of measuring the seasonality of incidence, a second set of covariates was obtained by replacing three meteorological variables by sine and cosine terms with period 12 months. Finally, uij is the random effect that captures the spatio-temporal autocorrelation in response data Yij, whose variance depends on the adjacency matrix.

Conditional autoregressive (CAR) priors [48] structure were used on uij and for (α, βk), non-informative normal prior distributions was used. Flat and conjugate priors were specified for uij using inverse gamma distributions with shape and scale parameters equal to 0.001. Markov chain Monte Carlo simulation was used to estimate the model parameters, sampling 300,000 times, with the first 150,000 as the burn-in, and keeping the results from every tenth iteration. The “ST.CARar” function of the R statistical software package CARBayesST (www.r-project.org) was used to fit the model. Convergence was assessed by trace plots and checked by the convergence Z-score diagnostic function [49]. The Watanabe-Akaike Information Criterion (WAIC) was used as a measure of goodness of fit [50].

Local Indicators of Spatial Association (LISA) were used to identify significant hotspots, coldspots, and outliers of dengue incidence at the sub-district level [51]. A hotspot is defined as an area that is surrounded by other high incidence areas, i.e. incidence is higher than the expected number given a random distribution of cases (so called high-high cluster). A coldspot is defined as an area surrounded by other low incidence areas (low-low cluster). Hotspot detection can be useful, even if the global pattern is not clustered. Moreover, case clusters that occur randomly can also have an influence on the spread of an infectious disease [52].

Results

General results

Dengue cases numbering 15,167 were reported over the 11-year period by all hospitals and clinics in Khon Kaen Province. Of these, there were 7461 dengue fever cases (49.2%) and 7706 severe dengue cases (50.8%), comprising both DHF and DSS. The demographic characteristics of patients are summarized in Table 1. Males represented the majority of patients (8057; 53.1%). Ages ranged from 4 months to 92 years old (median 13 years). The highest number of patients was in the 5–14 year-old age group (7758; 51.1%), followed by 15–29 years (5026; 33.1%) and 30–44 years (937; 6.2%). The proportion of older age groups (> 15 years), increased from nearly 20% of all cases in 2006 to more than 50% in 2016 (Fig. 1). The highest recorded disease incidence was in 2013, approximately 80 per 100,000 population (Fig. 2). Incidence was high during the rainy season (May–September), with July having the highest incidence (Fig. 3).

Table 1 Demographic characteristics of dengue reported cases in Khon Kaen Province, Thailand, 2006–2016
Fig. 1
figure1

Age distribution of reported dengue cases (DF, DHF and DSS) in Khon Kaen Province, Thailand, 2006–2016

Fig. 2
figure2

Monthly dengue incidence (a), dengue anomaly (b), rainfall (c), rainfall anomaly (d), temperature (e) and temperature anomaly (f) in Khon Kaen Province, Thailand, 2006–2016. DF = dengue fever, DHF = dengue hemorrhagic fever, DSS = dengue shock syndrome

Fig. 3
figure3

Mean monthly dengue incidence per 100,000 persons (a) and monthly average of rainfall (bar) and temperature (line) (b) in Khon Kaen Province, Thailand, 2006–2016

Association between dengue cases and climatic factors

Mean rainfall and maximum temperature were positively associated with dengue incidence, and minimum temperature was negatively associated, in terms of their point estimates (Table 2). However, among the three 95% credible intervals (CIs), only the one for maximum temperature excluded 1 (null effect). The rate ratio for maximum temperature was 1.055, implying 5.5% (95% CI 0.9–11.5%) increase in cases with an increase of 1 °C per month. The range of this variable was from 30.7 °C to 44.9 °C. The rate ratio for mean rainfall was 1.004, indicating that increasing rainfall by one unit (1 cm) per month would increase dengue incidence by about 0.4%. The Watanabe-Akaike Information Criterion (WAIC) for this model was 10,028.75. For the model with two sinusoid terms replacing the three meteorological variables, the WAIC was very similar, at 10028.23. This sinusoid terms had a peak to trough rate ratio of 5.8, and a peak in mid-July, i.e. a roughly six-fold difference in fitted incidence from mid-July to mid-January.

Table 2 Point estimates and 95% credible interval of the Bayesian Poisson regression model on number of all monthly dengue cases (DF, DHF and DSS) and covariates in Khon Kaen Province, Thailand, 2006–2016

The mean dengue incidence was high in the central northeastern sub-districts, around Khon Kaen city, and in the southwestern sub-districts of the province (red and orange in Fig. 4a). The distribution of the posterior means of the random effects (from the CAR model with meteorological variables) show some clustering, indicating that the variables in the model did not account fully for the spatial variation in the data (Fig. 4b). Posterior distribution plots are shown in Additional file 4. High clusters were present around Khon Kaen city and the southwestern portion of the province and low clusters were present in the northwestern area (Fig. 5), from the LISA analysis. When broken down by month, the incidences show the same clustering patterns, especially during July–August (Additional file 5).

Fig. 4
figure4

Mean dengue prevalence by sub-district (a) and spatial distribution of the posterior means of random effects for dengue (b) in Khon Kaen Province, Thailand, 2006–2016

Fig. 5
figure5

High and low clustering of dengue incidence in Khon Kaen Province, Thailand, 2006–2016

Discussion

The majority (~ 90%) of patients were below the age of 30 years. The trend during the study period showed that the proportion of dengue cases younger than 15 years declined from almost 80% in 2006 to below 50% in 2016. Dengue fever is generally more common in younger age groups [53], although there is evidence showing increasing incidence of more severe disease and outcomes among older age groups [54]. Our observations are also consistent with a population age shift, potentially influenced by demographic changes, such as the birth and death rates that show decreasing trends during 2011 and 2015 [55]. Thailand, in general, is undergoing a demographic transition where the proportion older adults are gradually increasing with an increase in median age of the general population. A higher proportion of adults will also increase the number of immune individuals (those with previous exposure to dengue virus) in the population, which might theoretically decrease the risk of dengue infection in younger people by providing alternative blood sources for infectious mosquitoes [56]. This age shift has also been observed in other Asian countries with a higher frequency of dengue cases among people 15 years of age and older [16]. Increases in disease incidence in older age groups may be explained by an increase in secondary infections and changes in circulating dengue virus serotypes [57], which have been shown to be important risk factors for severe clinical presentations [58,59,60,61,62].

There were clear seasonal patterns of dengue incidence in Khon Kaen Province during the study period. Dengue occurs throughout the rainy season, with 73% of cases reported between May and September. Although maximum temperature was associated with higher incidence (Table 2), the model with meteorological covariates had similar performance (in terms of the WAIC) to a non-mechanistic model, which simply fitted a sinusoidal pattern with a period of 12 months. In our study, a 1 cm increase in monthly rainfall was associated with a 0.4% increase in dengue incidence. In Timor Leste, results from similar modeling analyses showed a far larger effect: a 47% increase in incidence per 1 mm increase in annual rainfall [63]. Different climate patterns between Timor Leste and Thailand might explain these differences. Rainfall can affect the availability of mosquito larval habitats [34]. During rainy and dry periods of the year, permanent water containers are common in and around households; some located in toilet or bathroom spaces providing continuous year round mosquito production [35,36,37,38,39, 64]. Large water storage jars and tanks are the most commonly used containers in Thailand [64]. A study correlating rainfall and clinical dengue cases in Thailand from 2002 to 2003 also found that the dengue incidence was closely related with rainfall [65].

Temperature is another primary environmental risk factor for dengue transmission. Sea surface temperature (SST) changes, generally related to periodic El Niño Southern Oscillation effects, and air temperature, having more direct short-term effects, have both been shown to influence dengue incidence [63, 66]. Dengue incidence increased by 19.4% with a 1 °C increase in SST and 2.6% with a 1 °C increase in weekly maximum temperature in the Texas-Mexico border region [66]. Another study found that a 1 °C monthly increase in mean ambient temperature, dengue incidence increased by 0.7% [63]. In our study, the rate ratio for maximum temperature was 1.055 per °C, within the range from 30.7 °C to 44.9 °C. Higher temperatures enhance viral replication in the vector mosquito in a shorter amount of time and thus increase transmission potential of dengue viruses. A study of the extrinsic incubation period (EIP) of dengue serotype 2 in Aedes albopictus found that the virus remained in the midgut at 18 °C but could disseminate and invade the salivary glands at temperatures between 23 °C and 32 °C [67], thereby showing higher temperatures produce a shorter EIP and greater transmission potential. The strong and consistent relationships between climate, particularly rainfall and temperature, and the number of dengue cases have been used to develop prediction models to implement more timely dengue control measures [68, 69]. Relationships between dengue transmission and climatic variables have been examined in numerous studies, as shown above, but the question remains how to use such relationships in predicting impending outbreaks and applying effective interventions in time to avert them. User-friendly tools, such as the operational guide on Early Warning and Response System developed with support from the WHO/TDR and the European Union [70], are needed and will be tested in forthcoming work in Khon Kaen Province.

The highest dengue incidence seen in this study occurred in two areas of the province: around Khon Kaen Mueang District in the northeast, and in Manchakhiri and Khokphochai districts in the southwest. Mueang District includes the provincial capital and has the highest human population density, and in general, more conducive to dengue transmission. Manchakhiri and Khokphochai districts have lower population densities, but are, from our observations, seemingly similar to other districts in the province, i.e. vector species are present, larval habitats are plentiful, with a susceptible human population; therefore there must be other yet unexplored factors that support high dengue transmission in these two districts.

Although dengue incidence is influenced by rainfall and temperature, in our data there is no apparent spatial clustering of cases associated with the spatial variability in these environmental parameters. Rather, other factors such as urbanization are likely causes of the observed clustering effect [71]. However, population density, which was included in the regression model as a measure of urbanization, was not independently associated with dengue incidence. The residual spatial variation visible in Fig. 4b suggests that variables beyond those included in the spatial regression model are needed to explain differences in incidence between urban and rural subdistricts. Moreover, hotspots in more rural areas of southwestern Khon Kaen Province, further corroborate the influence of factors other than urbanization driving transmission. We do not know of any specific reasons for why these rural areas should have elevated dengue prevalence. One speculation could be that the lakes and swamps that are common in this area may provide suitable humidity for mosquitoes to thrive, but this was not studied here. Large changes in population size over time will affect outcomes. However, during 2000 and 2015, the average annual population growth rate in Thailand was less than 0.5% [72], which might not have affected the results substantially. Rural-urban migration is common in Thailand, with people drawn by, for example, better education, job opportunities, health facilities, standard of living, and wages [73]. Human movement is also an important factor in the dynamics of dengue transmission [74]. Adults are more likely to have greater mobility than younger age groups; therefore, to understand the circulation of the virus information on recent travel history and working conditions (location, time of work, etc.) is required. Elsewhere in Thailand, greater vulnerability to dengue infection has been observed in villages situated closer to urban centers [75]. Such neighboring effects are related to similarities in human behavior, development infrastructure, and ecological surroundings. Moreover, similar lifestyles and social interactions between neighboring areas are evident between villages that share social and religious centers such as schools, temples, mosques and community halls [75]. Hence, the results presented here are generalizable to most of northern Thailand, Laos, and Cambodia, and potentially Vietnam and Myanmar as well, under similar epidemiological settings.

Data collected from national surveillance systems come with inherent limitations, including underreporting and misreporting of symptomatic cases as well as the absence of subclinical and asymptomatic infections [76]. Moreover, dengue cases are seldom laboratory confirmed or identified to serotype. Another limitation of this study is inaccuracy, albeit minor, of the population denominators within sub-districts, as these were taken as fixed values from a single census (2010). Lastly, the possibility of travel-related infections was not determined in this study, which would provide potential misclassification bias. Nationally, the importance of travel-related dengue would vary by locality based on mobility. Obviously, we cannot exclude the possibility that some dengue infections were acquired outside the study area, thus potentially affecting the analysis and conclusions. However, if the general travel patterns had not changed significantly over the 11-year observation period, the dengue disease trends reported in this study would remain valid.

Conclusion

We examined the epidemiology of dengue in Khon Kaen Province, Thailand between 2006 and 2016. There was an increase in older age groups reporting symptomatic dengue. Symptomatic dengue disease in people > 15 years of age is now more common than in children in this province, an observation that has been seen in other Asian countries. This study used monthly sub-district level data to show that rainfall and temperature have significant effects on dengue transmission in the province. Spatial clustering of cases is partly associated with urban areas closer to Khon Kaen city and rural areas in the southwest of the province. However, the current analysis was not able to detect a close proxy factor to quantify a relationship between urbanization and dengue incidence. The data set awaits further analysis for temporal patterns of infection for use in disease prediction modeling and developing dengue early warning systems to guide vector control operations.

Availability of data and materials

The datasets used during the current study are available from the corresponding author on reasonable request.

Abbreviations

CAR:

Conditional autoregressive

DF:

Dengue fever

DHF:

Dengue hemorrhagic fever

DSS:

Dengue shock syndrome

EIP:

Extrinsic incubation period

ODPC7:

Office of Disease Prevention and Control 7

SST:

Sea surface temperature

References

  1. 1.

    Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496(7446):504–7.

    CAS  Article  Google Scholar 

  2. 2.

    WHO. Dengue: Guidelines for Diagnosis, Treatment, Prevention and Control.: New Edition. edn. World Health Organization, Geneva; 2009.

  3. 3.

    Christofferson RC, Mores CN. A role for vector control in dengue vaccine programs. Vaccine. 2015;33(50):7069–74.

    Article  Google Scholar 

  4. 4.

    Descloux E, Mangeas M, Menkes CE, Lengaigne M, Leroy A, Tehei T, et al. Climate-based models for understanding and forecasting dengue epidemics. PLoS Negl Trop Dis. 2012;6(2):e1470.

    Article  Google Scholar 

  5. 5.

    Wen TH, Lin NH, Chao DY, Hwang KP, Kan CC, Lin KC, et al. Spatial-temporal patterns of dengue in areas at risk of dengue hemorrhagic fever in Kaohsiung, Taiwan, 2002. Int J Infect Dis. 2010;14(4):e334–43.

    Article  Google Scholar 

  6. 6.

    WHO. Vector-borne diseases (Dengue). World Health Organization. Available from: http://www.who.int/mediacentre/factsheets/fs387/en/index2.html.

  7. 7.

    Limkittikul K, Brett J, L'Azou M. Epidemiological trends of dengue disease in Thailand (2000-2011): a systematic literature review. PLoS Negl Trop Dis. 2014;8(11):e3241.

    Article  Google Scholar 

  8. 8.

    Ponlawat A, Scott JG, Harrington LC. Insecticide susceptibility of Aedes aegypti and Aedes albopictus across Thailand. J Med Entomol. 2005;42(5):821–5.

    CAS  Article  Google Scholar 

  9. 9.

    Halstead SB, Yamarat C. Recent epidemics of hemorrhagic fever in Thailand. Observations related to pathogenesis of a “new” dengue disease. Am J Public Health Nations Health. 1965;55:1386–95.

    CAS  Article  Google Scholar 

  10. 10.

    Ministry of Public Health, Thailand. Epidemiology. Annual epidemiological surveillance report. Ministry of Public Health, Nonthaburi, Thailand. 1987.

  11. 11.

    Ministry of Public Health, Thailand. Epidemiology. Annual epidemiological surveillance report. Ministry of Public Health, Nonthaburi, Thailand. 1998.

  12. 12.

    Ministry of Public Health, Thailand.Epidemiology. Annual epidemiological surveillance report. Ministry of Public Health, Nonthaburi, Thailand. 2001.

  13. 13.

    Ministry of Public Health, Thailand. Epidemiology. Annual epidemiological surveillance report. Ministry of Public Health, Nonthaburi, Thailand. 2013.

  14. 14.

    Ministry of Public Health, Thailand. Annual epidemiological surveillance report. 2015. Ministry of Public Health, Nonthaburi, Thailand. Available from: http://www.boe.moph.go.th/Annual/AESR2015/.

  15. 15.

    Jeefoo P, Tripathi NK, Souris M. Spatio-temporal diffusion pattern and hotspot detection of dengue in Chachoengsao province, Thailand. Int J Environ Res Public Health. 2011;8(1):51–74.

    Article  Google Scholar 

  16. 16.

    Anker M, Arima Y. Male-female differences in the number of reported incident dengue fever cases in six Asian countries. Western Pac Surveill Response J. 2011;2(2):17–23.

    PubMed  PubMed Central  Google Scholar 

  17. 17.

    Ooi EE. Changing pattern of dengue transmission in Singapore. Dengue Bulletin. 2001;25.

  18. 18.

    Ooi EE, Goh KT, Gubler DJ. Dengue prevention and 35 years of vector control in Singapore. Emerg Infect Dis. 2006;12(6):887–93.

    Article  Google Scholar 

  19. 19.

    Yew YW, Ye T, Ang LW, Ng LC, Yap G, James L, et al. Seroepidemiology of dengue virus infection among adults in Singapore. Ann Acad Med Singap. 2009;38(8):667–75.

    PubMed  Google Scholar 

  20. 20.

    Ministry of Public Health, Thailand. Reporting of Priority Diseases Guideline. Ministry of Public Health, Nonthaburi, Thailand.Thailand. Available from: http://www.boe.moph.go.th/files/report/20121008_18818829.pdf.

  21. 21.

    Thawillarp S. Evaluation of possible dengue outbreak detection methodologies for Thailand, which one should be implemented? Johns Hopkins University; 2017.

    Google Scholar 

  22. 22.

    Gould EA, Higgs S. Impact of climate change and other factors on emerging arbovirus diseases. Trans R Soc Trop Med Hyg. 2009;103(2):109–21.

    CAS  Article  Google Scholar 

  23. 23.

    Sang S, Gu S, Bi P, Yang W, Yang Z, Xu L, et al. Predicting unprecedented dengue outbreak using imported cases and climatic factors in Guangzhou, 2014. PLoS Negl Trop Dis. 2015;9(5):e0003808.

    Article  Google Scholar 

  24. 24.

    Hopp MJ, Foley JA. Worldwide fluctuations in dengue fever cases related to climate variability. Clim Res. 2003;25:85–94.

    Article  Google Scholar 

  25. 25.

    Alto BW, Juliano SA. Temperature effects on the dynamics of Aedes albopictus (Diptera: Culicidae) populations in the laboratory. J Med Entomol. 2001;38(4):548–56.

    CAS  Article  Google Scholar 

  26. 26.

    Rueda LM, Patel KJ, Axtell RC, Stinner RE. Temperature-dependent development and survival rates of Culex quinquefasciatus and Aedes aegypti (Diptera: Culicidae). J Med Entomol. 1990;27(5):892–8.

    CAS  Article  Google Scholar 

  27. 27.

    Scott TW, Amerasinghe PH, Morrison AC, Lorenz LH, Clark GG, Strickman D, et al. Longitudinal studies of Aedes aegypti (Diptera: Culicidae) in Thailand and Puerto Rico: blood feeding frequency. J Med Entomol. 2000;37(1):89–101.

    CAS  Article  Google Scholar 

  28. 28.

    Rohani A, Wong YC, Zamre I, Lee HL, Zurainee MN. The effect of extrinsic incubation temperature on development of dengue serotype 2 and 4 viruses in Aedes aegypti (L.). Southeast Asian J Trop Med Public Health. 2009;40(5):942–50.

    CAS  PubMed  Google Scholar 

  29. 29.

    Liu-Helmersson J, Stenlund H, Wilder-Smith A, Rocklov J. Vectorial capacity of Aedes aegypti: effects of temperature and implications for global dengue epidemic potential. PLoS One. 2014;9(3):e89783.

    Article  Google Scholar 

  30. 30.

    Thammapalo S, Chongsuwiwatwong V, McNeil D, Geater A. The climatic factors influencing the occurrence of dengue hemorrhagic fever in Thailand. Southeast Asian J Trop Med Public Health. 2005;36(1):191–6.

    PubMed  Google Scholar 

  31. 31.

    Campbell KM, Lin CD, Iamsirithaworn S, Scott TW. The complex relationship between weather and dengue virus transmission in Thailand. Am J Trop Med Hyg. 2013;89(6):1066–80.

    Article  Google Scholar 

  32. 32.

    Lowe R, Cazelles B, Paul R, Rodo X. Quantifying the added value of climate information in a spatio-temporal dengue model. Stoch Env Res Risk A. 2016;30(8):2067–78.

    Article  Google Scholar 

  33. 33.

    Rozilawati H, Zairi J, Adanan CR. Seasonal abundance of Aedes albopictus in selected urban and suburban areas in Penang. Malaysia Trop Biomed. 2007;24(1):83–94.

    CAS  PubMed  Google Scholar 

  34. 34.

    Li CF, Lim TW, Han LL, Fang R. Rainfall, abundance of Aedes aegypti and dengue infection in Selangor, Malaysia. Southeast Asian J Trop Med Public Health. 1985;16(4):560–8.

    CAS  PubMed  Google Scholar 

  35. 35.

    Lambdin BH, Schmaedick MA, McClintock S, Roberts J, Gurr NE, Marcos K, et al. Dry season production of filariasis and dengue vectors in American Samoa and comparison with wet season production. Am J Trop Med Hyg. 2009;81(6):1013–9.

    Article  Google Scholar 

  36. 36.

    Bermudi PMM, Kowalski F, Menzato MM, Ferreira MDC, Passos W, Oku VJA, et al. Aedes aegypti breeding site in an underground rainwater reservoir: a warning. Rev Saude Publica. 2017;51:122.

    Article  Google Scholar 

  37. 37.

    Serpa LL, Costa KV, Voltolini JC, Kakitani I. Seasonal variation of Aedes aegypti and Aedes albopictus in a city of southeastern Brazil. Rev Saude Publica. 2006;40(6):1101–5.

    Article  Google Scholar 

  38. 38.

    Kittayapong P, Strickman D. Distribution of container-inhabiting Aedes larvae (Diptera: Culicidae) at a dengue focus in Thailand. J Med Entomol. 1993;30(3):601–6.

    CAS  Article  Google Scholar 

  39. 39.

    Tonn RJ, Sheppard PM, Macdonald WW, Bang YH. Replicate surveys of larval habitats of Aedes aegypti in relation to dengue haemorrhagic fever in Bangkok. Thailand Bull World Health Organ. 1969;40(6):819–29.

    CAS  PubMed  Google Scholar 

  40. 40.

    Sharma KD, Mahabir RS, Curtin KM, Sutherland JM, Agard JB, Chadee DD. Exploratory space-time analysis of dengue incidence in Trinidad: a retrospective study using travel hubs as dispersal points, 1998-2004. Parasit Vectors. 2014;7:341.

    Article  Google Scholar 

  41. 41.

    Cummings DA, Irizarry RA, Huang NE, Endy TP, Nisalak A, Ungchusak K, Burke DS. Travelling waves in the occurrence of dengue haemorrhagic fever in Thailand. Nature. 2004;427(6972):344–7.

    CAS  Article  Google Scholar 

  42. 42.

    van Panhuis WG, Choisy M, Xiong X, Chok NS, Akarasewi P, Iamsirithaworn S, et al. Region-wide synchrony and traveling waves of dengue across eight countries in Southeast Asia. Proc Natl Acad Sci U S A. 2015;112(42):13069–74.

    Article  Google Scholar 

  43. 43.

    Lolekha S, Tanthiphabha W, Sornchai P, Kosuwan P, Sutra S, Warachit B, Chup-Upprakarn S, Hutagalung Y, Weil J, Bock HL. Effect of climatic factors and population density on varicella zoster virus epidemiology within a tropical country. Am J Trop Med Hyg. 2001;64(3–4):131–6.

    CAS  Article  Google Scholar 

  44. 44.

    International Research Institute for Climate and Society: Climate Data Library. International Research Institute for Climate and Society. 2017.

  45. 45.

    NCEP Climate Forecast System Reanalysis (CFSR) Selected Hourly Time-Series Products, January 1979 To December 2010. Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory.

  46. 46.

    Funk C, Peterson P, Landsfeld M, Pedreros D, Verdin J, Shukla S, et al. The climate hazards infrared precipitation with stations-a new environmental record for monitoring extremes. Sci Data. 2015;2:150066.

    Article  Google Scholar 

  47. 47.

    National Statistical Office. The 2010 Population and housing census. Bangkok, Thailand 2010.

  48. 48.

    Rushworth A, Lee D, Sarran C. An adaptive spatiotemporal smoothing model for estimating trends and step changes in disease risk. J Royal Stat Soc Series C-Applied Stat. 2017;66(1):141–57.

    Article  Google Scholar 

  49. 49.

    Geweke J. In: Bernardo JM, Berger JO, Dawid AP, Smith AFM, editors. Evaluating the Accuracy of Sampling-Based Approaches to Calculating Posterior Moments. In Bayesian Statistics 4. Oxford: Clarendon Press; 1992.

    Google Scholar 

  50. 50.

    Burnham KP, Anderson DR. Multimodel inference: understanding AIC and BIC in model selection. Sociol Methods Res. 2004;33(2):261–304.

    Article  Google Scholar 

  51. 51.

    Anselin L. Local indicators of spatial association-LISA. Geogr Anal. 1995;27:93–115.

    Article  Google Scholar 

  52. 52.

    Pessanha JE, Caiaffa WT, Almeida MC, Brandao ST, Proietti FA. Diffusion pattern and hotspot detection of dengue in Belo Horizonte, Minas Gerais. Brazil J Trop Med. 2012;2012:760951.

    PubMed  Google Scholar 

  53. 53.

    Rigau-Perez JG, Vorndam AV, Clark GG. The dengue and dengue hemorrhagic fever epidemic in Puerto Rico, 1994-1995. Am J Trop Med Hyg. 2001;64(1–2):67–74.

    CAS  Article  Google Scholar 

  54. 54.

    Guha-Sapir D, Schimmer B. Dengue fever: new paradigms for a changing epidemiology. Emerg Themes Epidemiol. 2005;2(1):1.

    Article  Google Scholar 

  55. 55.

    National Statistical Office. Population and Housing census (Births and Deaths by Sex: 2011–2015, Khon Kaen province). Khon Kaen, Thailand 2018.

  56. 56.

    Cummings DA, Iamsirithaworn S, Lessler JT, McDermott A, Prasanthong R, Nisalak A, et al. The impact of the demographic transition on dengue in Thailand: insights from a statistical analysis and mathematical modeling. PLoS Med. 2009;6(9):e1000139.

    Article  Google Scholar 

  57. 57.

    Kanakaratne N, Wahala WM, Messer WB, Tissera HA, Shahani A, Abeysinghe N, et al. Severe dengue epidemics in Sri Lanka, 2003-2006. Emerg Infect Dis. 2009;15(2):192–9.

    Article  Google Scholar 

  58. 58.

    Guzman MG, Kouri G. Dengue and dengue hemorrhagic fever in the Americas: lessons and challenges. J Clin Virol. 2003;27(1):1–13.

    Article  Google Scholar 

  59. 59.

    Deen JL, Harris E, Wills B, Balmaseda A, Hammond SN, Rocha C, et al. The WHO dengue classification and case definitions: time for a reassessment. Lancet. 2006;368(9530):170–3.

    Article  Google Scholar 

  60. 60.

    Karyanti MR, Uiterwaal CS, Kusriastuti R, Hadinegoro SR, Rovers MM, Heesterbeek H, et al. The changing incidence of dengue haemorrhagic fever in Indonesia: a 45-year registry-based analysis. BMC Infect Dis. 2014;14:412.

    Article  Google Scholar 

  61. 61.

    Guzman MG, Alvarez M, Halstead SB. Secondary infection as a risk factor for dengue hemorrhagic fever/dengue shock syndrome: an historical perspective and role of antibody-dependent enhancement of infection. Arch Virol. 2013;158(7):1445–59.

    CAS  Article  Google Scholar 

  62. 62.

    Mohd-Zaki AH, Brett J, Ismail E, L'Azou M. Epidemiology of dengue disease in Malaysia (2000-2012): a systematic literature review. PLoS Negl Trop Dis. 2014;8(11):e3159.

    Article  Google Scholar 

  63. 63.

    Wangdi K, Clements ACA, Du T, Nery SV. Spatial and temporal patterns of dengue infections in Timor-Leste, 2005-2013. Parasit Vectors. 2018;11(1):9.

    Article  Google Scholar 

  64. 64.

    Vannavong N, Seidu R, Stenstrom TA, Dada N, Overgaard HJ. Effects of socio-demographic characteristics and household water management on Aedes aegypti production in suburban and rural villages in Laos and Thailand. Parasit Vectors. 2017;10(1):170.

    Article  Google Scholar 

  65. 65.

    Wiwanitkit V. An observation on correlation between rainfall and the prevalence of clinical cases of dengue in Thailand. J Vector Borne Dis. 2006;43(2):73–6.

    PubMed  Google Scholar 

  66. 66.

    Brunkard JM, Cifuentes E, Rothenberg SJ. Assessing the roles of temperature, precipitation, and ENSO in dengue re-emergence on the Texas-Mexico border region. Salud Publica Mex. 2008;50(3):227–34.

    Article  Google Scholar 

  67. 67.

    Liu Z, Zhang Z, Lai Z, Zhou T, Jia Z, Gu J, et al. Temperature increase enhances Aedes albopictus competence to transmit dengue virus. Front Microbiol. 2017;8:2337.

    Article  Google Scholar 

  68. 68.

    Phung D, Talukder MR, Rutherford S, Chu C. A climate-based prediction model in the high-risk clusters of the Mekong Delta region, Vietnam: towards improving dengue prevention and control. Tropical Med Int Health. 2016;21(10):1324–33.

    Article  Google Scholar 

  69. 69.

    Karim MN, Munshi SU, Anwar N, Alam MS. Climatic factors influencing dengue cases in Dhaka city: a model for dengue prediction. Indian J Med Res. 2012;136(1):32–9.

    PubMed  PubMed Central  Google Scholar 

  70. 70.

    WHO/TDR. Operational Guide: Early Warning and Response System (EWARS) for dengue outbreaks. Training in Tropical Diseases (TDR). World Health Organization, Special Programme for Research and Training in Tropical Diseases 2017;52.

  71. 71.

    Gubler DJ. Dengue, urbanization and globalization: the unholy trinity of the 21(st) century. Trop Med Health. 2011;39(4 Suppl):3–11.

    Article  Google Scholar 

  72. 72.

    United Nations, Department of Economic and Social Affairs, Population Division 2017. World Population Prospects: The 2017 Revision, custom data acquired via website: https://population.un.org/wpp/DataQuery/.

  73. 73.

    Katewongsa P. Benefits of rural-urban migration for migrants’ better life: a case study in Nang Rong, Buriram Province. Thailand Thammasat Review. 2015;18(1):63–81.

    Google Scholar 

  74. 74.

    Stoddard ST, Forshey BM, Morrison AC, Paz-Soldan VA, Vazquez-Prokopec GM, Astete H, Vilcarromero S, Elder JP, Halsey ES, et al. House-to-house human movement drives dengue virus transmission. Proc Natl Acad Sci U S A. 2013;110(3):994–9.

    CAS  Article  Google Scholar 

  75. 75.

    Tipayamongkholgul M, Lisakulruk S. Socio-geographical factors in vulnerability to dengue in Thai villages: a spatial regression analysis. Geospat Health. 2011;5(2):191–8.

    Article  Google Scholar 

  76. 76.

    da Fonseca BA, Fonseca SN. Dengue virus infections. Curr Opin Pediatr. 2002;14(1):67–71.

    Article  Google Scholar 

Download references

Acknowledgements

We thank the Office of Disease Prevention and Control, Region 7, Khon Kaen, Department of Disease Control, Ministry of Public Health, Thailand, for their kind cooperation and providing the dengue case data used in this study.

Consent for publications

Not applicable.

Funding

This project is funded by the Research Council of Norway (project no. 250443). The funding body had no role in study design, data collection, analysis, interpretation or preparation of the manuscript.

Author information

Affiliations

Authors

Contributions

TP designed the study, performed data analysis, interpretation and drafted the manuscript. BZ and UH contributed to study design, analysis and writing. CP, TE, SA, BF, MJB and NA guided the research and reviewed the manuscript. KT provided dengue data and guided the research. HJO proposed the research, made substantial contributions to conception, design and coordination of the study, participated in interpretation of data and drafting the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Hans J. Overgaard.

Ethics declarations

Ethics approval and consent to participate

This study is part of the DENGUE-INDEX Project funded by the Research Council of Norway (Project no. 250443), which received ethical approval from Khon Kaen University Ethics Committee for Human Research (HE591099–20/04/2016 and HE601221–01/09/2017), the London School of Hygiene and Tropical Medicine Ethical Committee, UK (LSHTM Ethics Ref: 10534–19/04/2016 and LSHTM Ethics Ref: 14275–16/08/2017), and the Regional Committee for Medical and Health Research Ethics, South East Norway (2016/357–15/04/2016 and 2017/1826–13/03/2018). Permission to use these data was approved by the Ministry of Public Health, Thailand.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Population density per sub-district, Khon Kaen province, Thailand, 2006 to 2016. Location of province in northeastern Thailand (inset). (PDF 128 kb)

Additional file 2:

Average monthly temperature (°C) per sub-district, Khon Kaen province, Thailand, January to December 2006–2016. (PDF 51 kb)

Additional file 3:

Average monthly rainfall (mm) per sub-district, Khon Kaen province, Thailand, January to December 2006–2016. (PDF 49 kb)

Additional file 4:

Posterior distribution plots of A) Mean rainfall, B) Minimum temperature, C) Maximum temperature, D) Age, E) Gender, and F) Population density. (PDF 265 kb)

Additional file 5:

Average monthly incidence of dengue (DF, DHF, and DSS) per 10,000 persons in Khon Kaen province, Thailand, January to December 2006–2016. (PDF 58 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Phanitchat, T., Zhao, B., Haque, U. et al. Spatial and temporal patterns of dengue incidence in northeastern Thailand 2006–2016. BMC Infect Dis 19, 743 (2019). https://doi.org/10.1186/s12879-019-4379-3

Download citation

Keywords

  • Dengue
  • Climate
  • Seasonal
  • Temperature
  • Rainfall
  • Thailand