Introduction
Dengue fever (DF) is a mosquito-borne viral infection caused by any one of four serotypes of dengue virus (DENV 1–4) [1]. It is a major public-health concern throughout tropical and sub-tropical regions. In recent decades, the incidence of DF has grown dramatically worldwide, making it the most important mosquito-borne viral disease [Reference Shepard2, Reference Bhatt3]. It has been estimated that 390 million dengue infections occur per year, of which 96 million manifest clinically (with any severity of disease) [Reference Bhatt3]. Brady et al. have estimated that more than 3.9 billion people in 128 countries are at risk of infection with dengue viruses [Reference Brady4]. The most seriously affected regions include Africa, the Americas, the Eastern Mediterranean, South-East Asia and the Western Pacific [Reference Bhatt3–5]. The drivers of the rapid dengue expansion include urbanisation, globalisation (travel and trade), the lack of effective mosquito control and climate change [Reference Murray, Quam and Wilder-Smith6, Reference Gubler7].
The first outbreak of DF in mainland China occurred in 1978 in the city of Foshan in Guangdong Province. Since then, DF outbreaks have been recorded sequentially in Hainan, Guangdong, Guangxi, Fujian, Zhejiang and Yunnan provinces [Reference Wu8]; the highest number of DF cases reported in China occurred in Guangdong Province [Reference Lai9]. In 2014, the most severe DF outbreak in history occurred in Guangdong Province, with more than 45 000 cases of infection [Reference Wu8–Reference Sun10]. In this outbreak, 37 359 cases of infection were reported in Guangzhou, the provincial capital of Guangdong Province, seven times the historical record, confirming Guangzhou as the most affected city [Reference Yang11]. DF outbreaks in China were previously thought to be imported [Reference Lai9], but recent studies suggest that DF may be endemic to China [Reference Lin12, Reference Luo13].
DENV are transmitted by Aedes mosquitoes, which are highly sensitive to climate, as temperature, precipitation and humidity, for example, can influence dengue transmission both directly and indirectly mediated by mosquito density [Reference Xu14, Reference Morin, Comrie and Ernst15]. The effects of climatic factors on DF incidence at different times and in different regions have been studied [Reference Earnest, Tan and Wilder-Smith16–Reference Chien and Yu20] and can vary dramatically according to region [Reference Morin, Comrie and Ernst15, Reference Vu21]. Many studies have investigated the relationship between DF incidence and climatic factors in Guangzhou, most have employed regression models, such as the Poisson model [Reference Sang22], the negative binomial model [Reference Shen23], the zero-inflated regression model [Reference Wang24] and the distributed lag non-linear models [Reference Xiang25]. Although there are some inconsistencies, the results of these studies have provided valuable implications for dengue risk assessment, prediction and prevention. However, the relationship between climate factors and dengue transmission is complex [Reference Morin, Comrie and Ernst15]; therefore, further research is needed.
Socioeconomic factors, such as urbanisation, also influence the incidence of DF [Reference Morin, Comrie and Ernst15, Reference Lowe18, Reference Lowe26]. One study investigated the effects of socioeconomic and environmental factors on the incidence DF at the township level in the Pearl River Delta economic zone in 2013 [Reference Qi27]. Another study examined the effects of socioecological factors on the spatial distribution of DF in the unprecedented 2014 outbreak in Guangzhou [Reference Cao28]. Both these studies employed annual data. However, some socioeconomic factors, such as the number of tourists, may vary dramatically in a year. To our knowledge, no study on the effects of such factors on the incidence of DF has been published.
In this study, we use a novel spatial-temporal method, the GeoDetector, to identify the potential climatic and socioeconomic factors associated with the DF incidence in Guangzhou. The results may be useful in understanding the relationship between DF incidence and risk factors in Guangzhou, in providing information to predict DF outbreaks, and in developing preventive measures.
Methods
Study area
Guangzhou, the city most affected by dengue in mainland China, is located on the south-east coast of China (22°26′–23°56′N, 112°57′–114°3′E, Fig. 1) and had a population of 13.08 million in 2014 [29]. As the capital of Guangdong Province, Guangzhou serves as the political, economic, scientific, educational, tourism and cultural centre in southern China. The city has a humid subtropical climate; summers are wet with high temperatures and high humidity, and winters are mild and comparatively dry.
Fig. 1. Geographic location of Guangdong Province and Guangzhou city in China, and meteorological stations used to interpolate the meteorological data in Guangzhou.
Materials
Guangzhou DF cases
In this study, monthly DF cases from January 2006 to December 2014 (Fig. 2) were obtained from the Chinese Center for Disease Control and Prevention (CDC, http://www.chinacdc.cn/). DF incidence was calculated by determining the ratio between the number of DF cases and the population size of the year. All four serotypes had been detected in Guangzhou, which was dominated by DENV 1–2; DENV 3 was first detected in 2009 and DENV 4 re-emerged in Guangzhou in 2010 [Reference Luo13, Reference Sang30].
Fig. 2. Guangzhou DF monthly cases from January 2006 to December 2014: (A) January 2007 to January 2013; (B) June 2006 to December 2006; (C) June 2013 to December 2013 and (D) June 2014 to December 2014.
Potential risk factors
Transmission of DENV requires multiple factors: virus must be present (for example, imported cases), sufficient susceptible humans to the virus and contact between humans and mosquito vectors. While the first and second factors are usually related to socio-economic factors (for example, population density, travellers go or from abroad), the third factor is influenced by environmental factors, including climatic factors, living conditions, land cover type, etc. In this study, we examined the effects of both climatic and socioeconomic risk factors (Fig. 3) on the DF incidence in Guangzhou.
Fig. 3. Potential risk factors to DF considered.
In Guangzhou, Aedes albopictus was the only vector species that transmitted DENV [Reference Lai9, Reference Li31]. The Breteau index (BI), calculated as the number of containers positive for Aedes mosquito larvae per 100 houses [1], was used to measure the mosquito population. The monthly BI of Guangzhou is shown in Figure 4.
Fig. 4. Guangzhou monthly BI between January 2006 and December 2014.
The climatic variables influencing the mosquito population considered include the monthly mean temperature (T), mean maximum temperature (T max), mean minimum temperature (T min), precipitation (P), the number of days with rainfall (R) and maximum daily precipitation (P max), as well as their lags of 1–3 months (V −j, j = 1, 2, 3, V denotes one of the former symbols). These data from meteorological stations located in Guangzhou and within 100 km of the borders of Guangzhou (Fig. 1) between 2006 and 2014 were obtained from the China Meteorological Data Sharing Service System (http://data.cma.cn/). We then employed the inverse distance weighted interpolation method to obtain the corresponding data for Guangzhou. The monthly mean, mean maximum and mean minimum temperatures are shown in Figure 5, the other climatic factors are shown in Figure 6.
Fig. 5. Guangzhou monthly mean, mean maximum and mean minimum temperatures between January 2006 and December 2014.
Fig. 6. Monthly precipitation (A), number of days with rainfall (B) and maximum daily precipitation (C) in Guangzhou between January 2006 and December 2014.
In addition to the climatic factors, socioeconomic factors were also considered to be potentially important drivers in DENV transmission [Reference Lowe18, Reference Lowe26]. In this study, we collected the data pertaining to monthly domestic (D) and overseas tourists (O) coming to Guangzhou from the Guangzhou Statistics Bureau (http://www.gzstats.gov.cn/). Figure 7 shows the variations in the number of tourists over time. Population density and the proportion of urban and rural dwellers were also considered potential risk factors to DENV transmission.
Fig. 7. Tourists coming to Guangzhou between January 2006 and December 2014.
GeoDetector
Several studies have reported the associations between DF incidence in Guangzhou and risk factors [Reference Sang22–Reference Xiang25]. However, these methods usually limited to some shortcomings, for example, the relationship between DF incidence and risk factors are usually nonlinear; the risk factors may have an interactive effect on DF incidence but cannot be fully reflected by regression models. The above limitations do not exist in the case of the GeoDetector model [Reference Wang32], of which, the basic idea is to measure the correspondence of the spatial temporal distribution of response variables (e.g. DF incidence) to that of suspected determinants (e.g. climatic factors). It assumes that if the suspected determinant is a disease risk, the spatial temporal distribution of the disease should be similar to that of the factor. That is to say, if a potential factor (X) causes a disease (Y), their temporal-spatial distributions tend to be consistent [Reference Wang32], as measured by the power of determinant:
$$q = 1-\displaystyle{1 \over {N\sigma ^2}}\mathop \sum \limits_{h = 1}^L N_h\sigma _h^2 $$where σ 2 denotes the variance of Y in the study area; N is the size of the population Y, which is composed of L strata (h = 1, 2, …, L) and 
$\sigma _h^2 $ stands for the variance within stratum h. The value of the q-statistic, ranging from 0 to 1, denotes the determinant power of a potential risk factor X. This means that the factor explains 
$q \times 100{\rm \%} $ of the DF incidence. The bigger it is, the more determinant power of the factor X. If q = 1, Y is determined by X completely. On the contrary, if q = 0, the factor X is unrelated to Y completely. In addition to extracting the determinant power of a single factor, the q-statistic can also be used to probe the determinant power of the interactive effect of two risk factors, X 1 and X 2. By comparing the interactive effect between two factors with the effect of their individual contributions to Y, the interactive effects are determined (shown in Table 1). Detailed information about the GeoDetector model can be found in [Reference Wang32–Reference Wang and Xu34].
Table 1. Types of interaction between two factors
The symbol of “
$\cap $” means the interaction of two factors, and 
$q(X_1\cap X_2)$ represents the q-value for the interaction between two factors X 1 and X 2.
The GeoDetector model can handle both quantitative and nominal data, and thus offer a novel approach for detecting interactions between the risk factors for infectious diseases, such as bacillary dysentery [Reference Xu35] and hand-foot-and-mouth disease [Reference Xu36–Reference Wang38]. In this study, we used three geographical detectors, i.e. factor detector, ecological detector and interaction detector to identify, respectively, the responsible factors, the relative importance factors and the interaction relationship for each two factors. The method was implemented using software downloaded from the GeoDetector website at http://www.geodetector.org/. All variables are temporal variables but without spatial variation. The population density and the proportion of urban and rural dwellers were stable in time within a year, while the climatic variables and tourists coming to Guangzhou varied monthly. To use the model, the risk factors measured by continuous data need to be discretised. The most commonly used discretisation methods include the equal interval method, the quantile method, the natural breaks method, the geometrical interval method, the standard deviation method, etc. [Reference Cao, Ge and Wang39]. Here, we used the q-statistic as the indicator for the assessment of the discretisation method and obtained the stratification information [Reference Cao, Ge and Wang39].
Results
Descriptive statistics
There were a total of 39 695 DF cases between January 2006 and December 2014. An obvious seasonal pattern in the number of DF cases can be seen. Most cases occurred in the July–November period, peaking in September and October (Fig. 2). The DF incidence also shows striking annual variations, three large outbreaks occurred in 2006, 2013 and 2014, with the highest number of cases recorded in 2014 (37 359 cases). The DF cases and the climatic and socioeconomic factors have been summarised (Table 2).
Table 2. Statistical description of incidence and potential risk factors of DF
‘s.d.’ denotes the standard deviation of the corresponding variable.
Results from the GeoDetector
We used the GeoDetector method to detect the determinant power of each risk factor and their interactive effects to the DF incidence in Guangzhou. The potential risk factors include six climatic factors (the monthly mean temperature, mean maximum temperature, mean minimum temperature, precipitation, number of days with rainfall and maximum daily precipitation) and four socioeconomic factors (the proportion of urban and rural dwellers, population density and overseas and domestic tourists). The determinant power of each factor and their interactive effects were qualified by the q-value calculated using the GeoDetector model. The association between the DF incidence and each factor – for example, a positive or negative relationship – was assessed using Pearson correlation coefficients.
The DF incidence demonstrates an apparent seasonal variation. This variation is correlated with the mosquito density and climatic factors and/or their lags of 1–3 months with different determinant powers. The q-values of BI and the climatic factors are shown in Table 3. The results show that the determinant power of BI is almost 60%, indicating a close relationship between BI and DF incidence. Among the climatic factors, the number of days with rainfall of the month before last has the highest determinant power, with a q-value of 0.898 (P < 0.01). However, the Pearson correlation was not significant (different from zero). We suppose that a nonlinear relationship exists between the DF incidence and the number of rainy days of the month before last. Regarding precipitation, the determinant power with the 2-month lag is the highest among the other lags, with a q-value of 0.486 (P < 0.01). The maximum daily precipitation in the last 2 months had the same q-values, 0.486. However, the Pearson correlation coefficients between DF incidence and them were not significant, meaning nonlinear relationships may exist between them. The interpretation for the nonlinear relationships can be found in the next section.
Table 3. q-values for mosquito density index and climatic factors and their lags of month
Note: all the q-values are significant from zero with P < 0.05.
The determinant power of the mean temperature, mean maximum temperature, mean minimum temperature and their lags of 1–3 months are similar, with q-values within the range of 0.38–0.50 (P < 0.01). The Pearson correlations between these variables and the DF incidence were not significant (P > 0.05). This may be because the impact of these climatic variables and these lags on the DF incidence were nonlinear. The Pearson correlation coefficients between DF incidence and each risk variable considered are shown in Table 4.
Table 4. Pearson correlation coefficients between DF incidence and variables related to climate
#P > 0.1; *0.05 < P < 0.1.
The determinant powers of socioeconomic factors were also calculated using the GeoDetector method. Both domestic and overseas tourists coming to Guangzhou have significant effects on DF incidence, with q values of 0.489 and 0.487, respectively (both P < 0.01). More tourists were associated with higher DF incidence. The Pearson correlation coefficients between domestic and overseas tourists and the DF incidence are 0.207 (P < 0.05) and 0.174 (0.05 <P < 0.1), respectively. For the population density and the proportion of urban and rural dwellers, neither has significant explanatory power, with q-values of 0.0723 (P = 0.1234) and 0.0671 (P = 0.1503).
From the former results, we found that, for each risk factor related to temperature, the factor with 3-months lag has the largest q-value; for the factors related to precipitation, is the one with 2-months lag. We analysed their interactive effects with BI, the domestic and overseas tourists. The results show that the coupled impact of every two factors considered was enhanced compared with the effect of their individual influences (Table 5). As we know, both higher temperature and more precipitation would provide more suitable conditions for mosquitoes, the results show that their interactive effect is enhanced. DF in Guangzhou was inspired by imported cases [Reference Sang30], travellers go to or from countries where DF is endemic have chance to import DENV. Coupled with greater mosquito population, or suitable environmental and climatic factors, more travellers would increase the risk of DF outbreak, caused enhanced interactive effects.
Table 5. The q-values for the interactive effect of different factors
EB, enhance (bivariate); EN, enhance (nonlinear).
Discussion
DF has been a growing public-health concern in mainland China since its emergence in Guangdong Province in 1978 [Reference Jin, Lee and Shu40]. The city of Guangzhou, the provincial capital of Guangdong Province and one of the largest modern cities in China, has had the highest DF incidence in mainland China [Reference Yang11]. In this study, the determinant powers of potential climatic and socioeconomic factors were examined. The results indicate that the number of rainy days of the month before last had the greatest determinant power; the BI, the temperature and its lags of 1–3 months, precipitation, the maximum daily precipitation and the domestic and overseas tourists had significant impact on DF incidence in Guangzhou. The interactive effects of the mean minimum temperature and the other climatic variables were enhanced compared with the effect of their individual effects.
There is an obvious seasonal pattern in the Guangzhou DF incidence. This seasonal variation is associated with climatic factors through direct and indirect effects [Reference Xu14, Reference Morin, Comrie and Ernst15]. There is an opposite effect of precipitation on mosquito populations. While precipitation provides habitats for the aquatic stages of the mosquito life cycle, and higher precipitation is associated with increased mosquito populations; intense rainfall may wash out breeding sites and therefore have a negative effect on mosquito populations [Reference Morin, Comrie and Ernst15]. This study shows that the number of days with rainfall of the month before last has great impact on the DF incidence, with a q-value of 0.898. This is consistent with the findings in [Reference Xu14], which show the number of days with rainfall to be a better predictor. The determinant powers of the other variables related to rainfall, the precipitation, the maximum daily precipitation and their lags of 1–3 months, were around 0.38 to 0.50. A non-significant positive association between DF outbreak risk and precipitation was detected in [Reference Xu14]. In this study, integrating the determinant powers calculated using the GeoDetector model with the Pearson correlation analysis, we found that nonlinear relationships exist between the DF incidence and the variables related to rainfall.
While precipitation plays an essential role in providing habitats for the aquatic stages of the mosquito, temperature has a direct effect on DENV replication within vectors and indirect effects on DENV transmission by influencing the mosquito development and survival [Reference Morin, Comrie and Ernst15]. This study also shows that temperature and its lags of 1–3 months are important factors that influence the DF incidence in Guangzhou. Higher temperature is associated with higher DF incidence, which is a positive association. Like precipitation, the relationships between DF incidence and the three variables related to temperature, as well as their lags of 1–3 months, may be nonlinear. It is known that mosquito population dynamics, the oviposition rates and transition and mortality in different stages of the life cycle are considerably influenced by temperature [Reference Morin, Comrie and Ernst15, Reference Yang41–Reference Delatte43]. In short, an ideal range for survival and development through all phases of the mosquito life cycle occurs, for example, at 20–30 °C, as stated in [Reference Tun-Lin, Burkot and Kay44]. The development of mosquitoes is often slower in both cooler and higher temperatures [Reference Yang41, Reference Rossi, Ólivêr, Massad, Ferreira and Godoy42], although the critical limiting values may be different [Reference Morin, Comrie and Ernst15]. Temperature is also a key component in the ecology of DENV in that increasing temperature can accelerate viral replication within the vector and shorten the extrinsic incubation period (EIP) [Reference Chan and Johansson45, Reference Tjaden46]. Because of the variation in temperature throughout the day in nature, the effect of diurnal temperature ranges (DTRs) on the susceptibility of Ae. aegypti to DENV to DENV has also been explored in [Reference Lambrechts47]. The authors found that compared with moderate DTRs or constant temperature, larger DTRs decreased the probability of infection of Ae. aegypti, but the EIP was unchanged [Reference Lambrechts47]. Synthesising the above results, the overall effect of temperature on DENV ecology is complex; increasing temperature may accelerate parts of the viral transmission cycle, while other variables may become limited by higher temperatures [Reference Morin, Comrie and Ernst15]. Therefore, the nonlinear relationship between DF incidence and temperature is readily comprehensible.
DF incidence is also influenced by socioeconomic factors, such as population density, road density and land cover [Reference Qi27, Reference Cao28, Reference Ren48]. In this study, we examined the determinant power of domestic and overseas tourists coming to Guangzhou, the population density and the proportion of urban and rural dwellers. The results show that both domestic and overseas tourists coming to Guangzhou have a high determinant power of 0.489 and 0.487, respectively. There was a significant positive correlation between the tourists and DF incidence. The situation is similar to that in Australia, where there is an increasing trend in DF incidence as the number of overseas visitors increases [Reference Akter49]. Residents in Guangdong Province keep close connections with South East Asian countries that are in the DF-endemic regions [Reference Fan50]. Therefore, more travellers from these countries increase the risk of DENV transmission. We also explored the influence of population density and the proportion of urban and rural dwellers. The results show that these factors had non-significant determinant power in relation to DF incidence. This is inconsistent with the findings in [Reference Qi27], which analysed the relationships on a finer spatial scale. This inconsistency may be because the data we used were aggregated over large spatial scales.
The findings in this study have implication for the risk assessment, early warning of DF in Guangzhou, and to implement preventive strategies. Special attention should be paid to the rainy days, which has the largest determinant power on DF incidence. Besides mosquito density and climatic factors, which have been emphasised in the literatures [Reference Sang22–Reference Xiang25], travellers to the city, especially to or from countries where DF is endemic, should also be noted.
Our study aims to detect the impact of climatic and socioeconomic factors on DF incidence in Guangzhou. However, there were some limitations. First, imported DF cases were not excluded due to a lack of information. Nevertheless, this problem has little influence on the results because the imported cases only accounted for <1% of all DF cases in Guangzhou in 2006–2014 [Reference Sang22, Reference Xiang25]. Second, we used aggregated data for the whole city, while the associations between DF and risk factors are site-specific [Reference Morin, Comrie and Ernst15, Reference Vu21]. This may result in some detailed information been submerged; therefore, finer-scale research is needed in the future. Third, unreported and inapparent cases are an inherent limitation in the DF surveillance data. However, our results still have significance in the risk assessment of DENV transmission, in the prediction of DF outbreaks, and in implementing preventive DF reduction strategies.
Acknowledgements
This study was supported by the Fundamental Research Funds for the Central Universities (no. 9161017003 and 9161018003) and the State Key Laboratory of Resources and Environmental Information System, the Undergraduate Innovative Test Program (no. 201710079087).
Conflict of interest
The authors declare no conflict of interest.
Open access
