Skip to main content

Species-specific climate Suitable Conditions Index and dengue transmission in Guangdong, China

Abstract

Background

Optimal climatic conditions for dengue vector mosquito species may play a significant role in dengue transmission. We previously developed a species-specific Suitable Conditions Index (SCI) for Aedes aegypti and Aedes albopictus, respectively. These SCIs rank geographic locations based on their climatic suitability for each of these two dengue vector species and theoretically define parameters for transmission probability. The aim of the study presented here was to use these SCIs together with socio-environmental factors to predict dengue outbreaks in the real world.

Methods

A negative binomial regression model was used to assess the relationship between vector species-specific SCI and autochthonous dengue cases after accounting for potential confounders in Guangdong, China. The potential interactive effect between the SCI for Ae. albopictus and the SCI for Ae. aegypti on dengue transmission was assessed.

Results

The SCI for Ae. aegypti was found to be positively associated with autochthonous dengue transmission (incidence rate ratio: 1.06, 95% confidence interval: 1.03, 1.09). A significant interaction effect between the SCI of Ae. albopictus and the SCI of Ae. aegypti was found, with the SCI of Ae. albopictus significantly reducing the effect of the SCI of Ae. aegypti on autochthonous dengue cases. The difference in SCIs had a positive effect on autochthonous dengue cases.

Conclusions

Our results suggest that dengue fever is more transmittable in regions with warmer weather conditions (high SCI for Ae. aegypti). The SCI of Ae. aegypti would be a useful index to predict dengue transmission in Guangdong, China, even in dengue epidemic regions with Ae. albopictus present. The results also support the benefit of the SCI for evaluating dengue outbreak risk in terms of vector sympatry and interactions in the absence of entomology data in future research.

Graphical Abstract

Background

The worldwide incidence of dengue has risen rapidly in recent decades. According to model estimates, 390 million dengue virus infections occur yearly [1] and 3.9 billion people are at risk of infection [2]. The number of dengue cases reported to the WHO have climbed from 505,430 in 2000 to over 2.4 million in 2010 and 5.2 million in 2019 [3]. To date, dengue is the most prevalent human arboviral infection worldwide. Dengue is an arthropod-borne virus (arbovirus) infection caused by four dengue virus (DENV) serotypes (DENV-1, -2, -3 and -4), with short-lived cross-protection between serotypes [4]. DENV is carried and spread by Aedes mosquito species, among which Aedes aegypti is the principal dengue vector that transmits the virus and causes epidemics. Other Aedes mosquito species known to be vectors of dengue include Ae. albopictus, Ae. polynesiensis and Ae. scutellaris; however, the latter two species are considered to have a limited capacity to serve as dengue vectors [5]. Aedes aegypti and Ae. albopictus, once zoophilic forest species from sub-Saharan Africa and Asia, respectively, are now widespread globally due to the rapid growth of the human population, intensified international travel in the last century and, more recently, climate change [6, 7].

Aedes aegypti is well-adapted to urban areas in tropical regions, while Ae. albopictus can thrive in peridomestic and rural habitats in tropical, subtropical and temperate climates [8]. In China, Ae. aegypti was once only present in the areas of the tropical zone below 22°N, such as Hainan Island, the Leizhou peninsula in Guangdong province and the coastal area along Beibu Gulf, Guangxi province. However, since the beginning of the present century, Ae. aegypti has also been found in some bordering counties in Yunnan province, located between 22°N and 25°N [9]. In contrast, Ae. albopictus has a wide distribution in both tropical and temperate regions of China, ranging from 41°N to the country’s southern reaches [10]. Based on historical climatic data from the last decade, most areas where Ae. albopictus has become established had an annual mean temperature > 11 °C and annual mean precipitation > 500 mm [11].

Aedes aegypti is considered to be the principal dengue vector and has contributed to a global resurgence of dengue epidemics in the past three decades, supported by evidence that all major epidemics of dengue have occurred only in areas where Ae. aegypti is found [12]. Conversely, Ae. albopictus is considered to play a relatively minor role in the transmission of DENV compared to Ae. Aegypti and is not considered an ‘efficient’ dengue epidemic vector [12, 13]. In China, dengue outbreaks in Hainan province in 1980 and 1985–1986 and in Xishuangbanna of Yunnan province in 2013 were typically attributed to Ae. aegypti [14, 15]. However, outbreaks, largely of classical dengue (i.e. few cases of severe disease), in Guangdong, Fujian and Zhejiang provinces between 2004 and 2010 were caused solely by Ae. albopictus in the absence of Ae. aegypti [16, 17]. In 2014, China’s largest dengue outbreak recorded to date (or at least since dengue became a notifiable disease in 1989) occurred in Guangdong, where Ae. albopictus appears to be the sole dengue vector [18].

A variety of factors affect the spatial and temporal dynamics of mosquito populations and therefore of dengue fever transmission patterns [19]. Temperature, rainfall and humidity affect the emergence and viability of eggs, and the size and longevity of adult mosquitoes and their dispersal [20]. Unplanned urbanisation, high population density and overcrowding favour the spread of mosquito breeding sites and infection [19]. The WHO recommends that many dengue-endemic nations use vector surveillance to provide a quantitative evaluation of fluctuations in the number and geographical distribution of dengue vector populations, with the purpose of forecasting outbreaks and assessing management options [21]. However, dengue transmission usually depends on a range of factors, including susceptible human host population(s), mosquito density and climate, rather than a fixed entomologic threshold [22]. Additionally, to date there has been minimal evidence of meaningful correlations between vector biomass indices and dengue transmission, which might be used to forecast epidemics consistently [23]. Furthermore, the geographical variability in the association between dengue vector abundance and transmission or outbreak occurrence is poorly understood, which suggests that vector population indices may have limited capacity in forecasting dengue transmission increases [24]. Therefore, other locality-specific data may have better potential as predictors.

Climate is a crucial environmental determinant of vector geographical distribution and vectorial capacity. Climate factors such as temperature and precipitation and their seasonal patterns can fundamentally affect mosquito population dynamics and individual features relevant to vector biology, such as development, reproduction and survival, and consequently DENV transmission patterns [19, 20, 25]. A sufficiently warm ambient temperature is also critical for virus replication and dissemination to the salivary glands in female mosquitoes. Conversely, cooler temperatures slow this viral amplification, and if the development time of the pathogen exceeds the lifespan of the infected mosquito, transmission cannot occur [8, 26]. Previous studies have estimated suitable temperatures that support dengue transmission [2730]. Mordecai et al. synthesised existing empirical evidence on temperature suitability for vector and pathogen survival and used a mechanistic approach to ascertain suitable temperatures for vector capacity in dengue transmission. According to empirical data and probability models, suitable temperature ranges are between 17.05 °C and 34.61 °C for Ae. aegypti and between 15.84 °C and 31.51 °C for Ae. albopictus [31]. Based on these identified suitable temperature ranges for dengue vectorial capacity, we previously developed the Suitable Conditions Index (SCI) for Ae. aegypti and Ae. albopictus, respectively, to define their potential vector capacity. The SCI ranks geographical locations based on their climatic suitability for each vector species that transmits dengue and theoretically defines the spatial parameters for transmission [32]. By incorporating the SCI into an outbreak epidemiologic study, it is possible to investigate the potential value of SCI data for dengue transmission or outbreak prediction. In the study reported here, we used the autochthonous dengue case of the Guangdong dengue outbreak in 2014 to assess the potential of a climate-based SCI as a local dengue outbreak early warning index.

Methods

Study area

Guangdong (109°39′53″E to 117°18′51″E and 20°13′24″N to 25°31′10″N) Province, southern China, covers approximately 179,800 km2 in area. It is the most populous province in China, with 126.01 million inhabitants (7th National Population Census), accounting for 8.93% of Mainland China’s population [33] (Fig. 1). Guangdong has a humid subtropical climate, with short, mild and relatively dry winters and long, hot and very wet summers, influenced by the Asian monsoon season [34].

Fig. 1
figure 1

Study area and population density in Guangdong province, 2010

Data collection

Data on autochthonous and imported dengue cases reported from Guangdong in 2014 were collected from the National Notifiable Infectious Diseases Reporting Information System. Dengue has been a statutory notifiable communicable disease in China since 1989 [35]. Dengue is diagnosed according to the national surveillance protocol with standardised case definitions, including clinically diagnosed and laboratory-confirmed dengue cases. An imported dengue case was defined when the patient had travelled abroad to a dengue-endemic country within 15 days of the onset of illness. In some cases, importation was defined based on laboratory results showing that the infecting DENV had genomic sequence similarity with viruses isolated from the putative source region where the patient had travelled [36]. In the absence of meeting the criteria for an imported case, a dengue case was considered to be autochthonous dengue case.

The SCI was developed in our previous study [32] to define climate suitability for Ae. aegypti and Ae. albopictus, respectively. Based on the optimal temperature ranges for vectorial capacity to transmit DENV identified by Mordecai et al. using robust mechanistic models [31], the SCI was defined as the product of the number of suitable months (where the monthly minimum and maximum temperatures fall within the determined temperature range) and the average monthly precipitation (as a weight). The input data for the SCI model were available on an annual interval and at the county scale. Therefore, we produced SCI values for each of Ae. aegypti and Ae. albopictus in 2014 for 124 counties/districts in Guangdong Province. In this study we used 2014 dengue outbreak data in Guangdong, which is the worst dengue outbreak in the last few decades in China. The dengue outbreak data across 124 districts/county levels provides a unique opportunity to identity the relationship between SCI and dengue cases at spatial distribution.

Normalised Difference Vegetation Index (NDVI) data were extracted from the Moderate Resolution Imaging Spectroradiometer (MODIS) website of the US National Aeronautics and Space Administration (NASA; https://www.nasa.gov/). NDVI data are used to characterise and monitor the global distribution of vegetation conditions and are often used to model global and regional climate [37]. We utilised MODIS Vegetation Indices (MOD13A3 product) version 6 data provided monthly at the 1-km spatial resolution, then aggregated annually and at the county scale [38]. In addition, human population data were obtained from the Sixth National Population Census data (2010), which is available from the National Bureau of Statistics of China (http://www.stats.gov.cn/enGliSH/).

Negative binomial regression models

The negative binomial (NB) regression model was used to assess the relationship between autochthonous dengue cases and potential predictors, including SCI for each vector species, NDVI, human population density and the number of imported dengue cases. The NB regression model was also used to account for the over-dispersed count of dengue case (Table 1). A log-scale population was used as an offset to control for variation in human population size. Univariable and multivariable analyses were implemented successively. In the univariable analyses, we considered each potential predictor in an NB regression model to identify significant predictors. In a multivariable analysis, we examined the potential effect of significant predictors on dengue case occurrence (presented as incidence rate ratio [IRR]). Model goodness-of-fit was assessed by the Bayesian information criterion (BIC), with smaller values of BIC indicating a better fit.

We also examined the potential interaction between the SCI for Ae. albopictus and the SCI for Ae. aegypti (SCI-Ae. aegypti × SCI-Ae. albopictus). We assumed that the SCI was a proxy for vector density, with a high SCI value representing high vector density and vice versa. The interaction between the two vectors was examined by assessing how the effect of the SCI of a specific vector on autochthonous dengue case occurrence changed when the SCI of another vector changed. The difference in species-specific SCI (defined as ‘SCI-Ae. Aegypti—SCI-Ae. albopictus’) was also fitted to the NB regression model. As Ae. albopictus has a different niche to Ae. aegypti in terms of feeding preferences and, critically for the SCI, a different thermal range, it was reasonable to assume that a more considerable difference may indicate a reduced probability of ecological constraints between the two mosquito species. We therefore used the SCI difference to explore the correlation between the probability of ecological constraints and dengue transmission risk. The regression analysis was implemented in Stata/SE 16.0 (StataCorp, College Station, TX, USA) using “xtnbreg” (fixed-effects, random-effects and population-averaged negative binomial models, respectively) [39].

Results

Descriptive analyses

Summary statistics for autochthonous dengue cases, imported cases, NDVI, human population density, the SCI for Ae. aegypti and Ae. albopictus, respectively, are presented in Table 1. In 2014, there were 44,939 (44,862 autochthonous and 77 imported) dengue cases recorded from 107 of the 124 counties in Guangdong Province, with the number of all dengue cases per county/district ranging from 0 to 11,846 (mean number of cases per county/district: 362.41, standard deviation [SD]: 1388.57). The mean annual dengue incidence averaged 28.9 cases per 100,000 population, with the highest incidence (532.6/100,000) being reported from Baiyun District, followed by Liwan (497.3/100,000), Yuexiu (414.1/100,000) and Haizhu (385.2/100,000), all of which are in Guangzhou, the capital city of Guangdong (Fig. 2a).

Table 1 Summary statistics for autochthonous dengue cases, imported dengue cases, Normalised Difference Vegetation Index, human population density, Suitable Conditions Index for Ae. aegypti and Ae. albopictus, respectively, in 124 counties/districts of Guangdong Province, 2014
Fig. 2
figure 2

Geographical distribution of autochthonous dengue incidence (a), Suitable Conditions Index (SCI) for Aedes aegypti (b) and SCI for Aedes albopictus (c) at a county level in Guangdong Province, China, 2014

Guangdong was found to have the most suitable climatic conditions for dengue transmission for each vector (Fig. 2b, c). The SCI for Ae. aegypti was highest in the Leizhou peninsula below 22°N and in western coastal areas of the province between 22°N and 23°N, with a gradual reduction in suitability moving further north, and was lowest in the most northern part of the study area (Fig. 2b). The SCI for Ae. albopictus was less likely to be distributed evenly across low- to high-latitude areas (Fig. 2c). Regions with the highest SCI for both vectors overlapped in the western coastal areas below 23°N (Fig. 2b, c).

The relationship between autochthonous dengue incidence, latitude, imported cases and SCI for the two vector species is presented in Fig. 3a–d. The correlation between autochthonous dengue incidence and latitude was assessed against the SCI for Ae. aegypti and SCI for Ae. albopictus, respectively; a consistently higher value of autochthonous dengue incidence and SCI was found in areas at 23°N (Fig. 3a, b). The relationship between imported dengue cases, SCI for two mosquito vectors and autochthonous dengue cases indicated a strong relationship between autochthonous and imported dengue cases, which may remain independent over various levels of SCI for each vector species (Fig. 3c, d).

Fig. 3
figure 3

The relationship between autochthonous dengue incidence, latitude, imported dengue cases and Suitable Conditions Index (SCI) in Guangdong Province, China, 2014. Relationship between autochthonous dengue incidence, latitude and SCI for Aedes aegypti (a), relationship between autochthonous dengue incidence, latitude and SCI for Aedes albopictus (b), relationship between autochthonous dengue incidence and imported dengue cases and SCI for Aedes aegypti (c), relationship between autochthonous dengue incidence and imported dengue cases and SCI for Aedes albopictus (d), and relationship between autochthonous dengue incidence and SCI for Aedes aegypti over various levels of SCI for Aedes albopictus (e).

The relationship between autochthonous dengue incidence and SCI for Ae. aegypti was assessed when the SCI for Ae. albopictus was held constant at four levels ranging from low to high (34.9–58.1, 58.2–63.6, 63.8–84.5 and 85.5–141.8), and the results are shown in Fig. 3e. There was a moderate relationship between dengue incidence and SCI for Ae. aegypti when the SCI of Ae. albopictus was held constant at 34.9–58.1 (e.g. minimum to mean, with 1 SD) (Fig. 3e).

Spearman’s correlation and bivariate Moran’s I statistic

Spearman’s correlation coefficients between any two variables considered in this study are presented in Table 2. There was a moderate relationship between autochthonous dengue incidence and the SCI for Ae. aegypti (rs = 0.41, P = 0.000), while there was a weak correlation between incidence and the SCI for Ae. albopictus (rs = -0.18, P = 0.000). There was a moderate correlation between autochthonous dengue incidence and difference in SCI (rs = 0.53, P = 0.000). There was a moderate correlation between autochthonous dengue incidence and the other variables, including imported cases, NDVI and human population density. There was a moderate relationship between the two SCIs (rs = 0.43, P = 0.000). We used bivariate Moran’s I [40] to further test the spatial relationship between the two SCIs. The global and local bivariate Moran’s I between the SCI for Ae. albopictus was significantly associated with the SCI for Ae. aegypti (Moran’s I = 0.426, P = 0.001) (Additional file 1: Figure S1).

Table 2 Correlation between autochthonous dengue annual incidence and exploratory variables

NB Regression models

The results of univariable analyses (Model I in Table 3) indicated that the SCI for Ae. aegypti, the difference between the SCI for Ae. aegypti and the SCI for Ae. albopictus, NDVI and imported cases were significant predictors of dengue transmission. The results of the multivariable models without the interaction terms (Models II and III) indicated that autochthonous dengue cases increased by 4% (IRR: 1.04, 95% confidence interval [CI] 1.02, 1.06) for a 1-unit increase in SCI for Ae. aegypti, but decreased by 3% (IRR: 0.97, 95% CI: 0.96, 0.99) for every 1-unit increase in SCI for Ae. albopictus. Autochthonous dengue cases increased by 28% (IRR: 1.28, 95% CI: 1.02, 1.60) and by 4% (IRR: 1.04, 95% CI: 1.02, 1.06) for every 1-unit increase in imported cases and difference in SCI, respectively. Dengue cases decreased by 3% (IRR: 0.97, 95% CI: 0.96, 0.99) for a 1-unit increase in NDVI. A significant interaction effect between the SCI of Ae. albopictus and the SCI of Ae. aegypti was found, with the SCI of Ae. albopictus significantly reducing the effect of the SCI of Ae. aegypti on autochthonous dengue cases (Table 3: Model IV). To avoid possible collinearity between NDVI and population density (rs = −0.73, P = 0.000; Table 2; more importantly, the population density was not a significant predictor in multivariable models), the population density was excluded from the NB regression analyses.

Table 3 Incidence rate ratio and 95% confidence interval for potential predictors of autochthonous dengue case in Guangdong, 2014

Discussion

In the present study we investigated the association between vector species-specific SCI and autochthonous dengue cases in 2014, Guangdong province, China, and the value of SCI as a predictor of dengue outbreak risk in this geographical region. We found that the SCI was statistically associated with autochthonous dengue cases, implying that the SCI might be beneficial in forecasting the probability of local dengue transmission or outbreak occurrence. We also identified a statistically significant interaction between the SCI for Ae. albopictus and the SCI for Ae. aegypti, which invites the hypothesis that vector species interactions may also influence DENV transmission.

The SCI for Ae. aegypti was positively associated with autochthonous dengue transmission. However, dengue incidence was also negatively associated with the interactive effect between the SCI for Ae. albopictus and that for Ae. aegypti, which indicates that the effect of SCI for Ae. aegypti may change over the various levels of the SCI for Ae. albopictus. As shown in Fig. 3e, the region has a higher risk of dengue transmission if the region has a low SCI for Ae. albopictus but a high SCI for Ae. aegypti.

Aedes aegypti lives within or near human dwellings and tends to breed in urbanised ecological niches, whereas Ae. albopictus is more frequently found in rural areas and breeds outdoors. In addition, Ae. aegypti almost always feeds on humans, while Ae. albopictus usually bites humans and animals opportunistically [3]. These behavioural differences may indicate that Ae. aegypti has more opportunities to increase human-mosquito contact and thus increase the probability of DENV transmission. Although these two mosquito species are susceptible to infection by DENV, it the authors of a study conducted in Singapore reported that 6.9% of field-caught Ae. aegypti but only 2.9% of Ae. albopictus were DENV virus-positive [41]. Aedes aegypti is the primary dengue vector that can efficiently establish an outbreak. By contrast, Ae. albopictus plays a modest role in DENV transmission and is not an ‘efficient’ epidemic vector [12, 13]. The authors of another study conducted in Singapore reported that dengue-positive Ae. aegypti were detected 6 weeks before a dengue outbreak, but that infected Ae. albopictus did not appear until the number of cases was rising [42]. These differences in crucial characteristics, including behaviour and virus transmission capacity, may cause discrepancies in the roles of these two mosquito species in dengue outbreaks. There are still several significant gaps in our understanding of the correlations between climate suitability and virus transmission probability and intensity. Although temperature-dependent transmission models have been proven to be valuable for predicting vector spread and dengue transmission, most of these modelling studies found a very broad correlation of transmission occurring within a temperature range [25, 43]. Few models can sufficiently capture transmission dynamics to provide reliable predictions. Human mobility and behaviour, urban development and microclimates, aquatic habitat supply and vector control techniques may all impact dengue transmission [19, 4449]. These characteristics may be as important as climate suitability in determining dengue transmission risk. Moreover, transmission dynamics and bionomic responses differ across vector species, mosquito populations and DENV serotypes [23, 50, 51]. While some potential confounders have been adjusted in our study, other elements may contribute to dengue transmission that we have not considered. Finally, whether or not primary mosquito vectors are present, optimal climatic conditions for dengue vector species may play a substantial role in an outbreak. This suggests the possibility that vector species are a prerequisite for initiating and establishing a dengue outbreak while suitable conditions may expand and maintain the outbreak through complicated effects on vector sympatry or interactions.

Although the SCI was not equivalent to mosquito species absence/presence or densities revealed only by entomological surveillance data, specific vector mosquitoes should only occur at a suitable habitat with favourable meteorological conditions. Our findings indicate the possibility that the SCI might be indirectly associated with autochthonous dengue cases through the variation in dengue vector population, which means that the SCI might be beneficial in the prediction of dengue transmission or outbreaks in the absence of entomological data. The SCI may therefore be meaningful in some dengue-endemic countries where entomological surveillance data are usually unavailable and, more importantly, where there has been no reliable indication of any consistent association between entomological indices and dengue cases [23].

Our results suggest a significant interaction between the SCI for Ae. aegypti (main effect) and the SCI for Ae. albopictus, which means the expected effect of the SCI for Ae. aegypti on autochthonic dengue cases may change over the various levels of the SCI for Ae. albopictus. Specifically, a 1-unit increase in the SCI for Ae. aegypti may indicate a higher risk of dengue in a region with a low SCI for Ae. albopictus than in a region with a high SCI for Ae. albopictus. Given that regions with a high SCI for each of these two Aedes vectors have some overlaps in Guangdong, the coexistence of these two vectors is theoretically possible. In areas sympatric for the two species, declining population and displacement of Ae. aegypti have been reported because of the superiority of Ae. albopictus in the competition for resources at the larval stage and asymmetric sterilisation at the adult stage following interspecific mating [52]. It has also been supposed that Ae. aegypti prefers to reproduce in areas devoid of Ae. albopictus [53]. As far as our knowledge, there is no sufficient evidence of vector competitiveness from field study in Guangdong. However, recent evidence from a laboratory experiment suggests competitive displacement is theoretically possible. In that laboratory experiment, larvae of the two species were mixed, allowed to emerge and then cycled through six generations, with the results showing that Ae. aegypti from the Leizhou peninsula were suppressed (a combination of generation time and abundance was used to judge competitiveness) by Ae. albopictus from Guangzhou but not by Ae. albopictus from other cities [54]. This finding warrants further exploration because the sympatric habit for the two species is expected to expand to additional areas in China [55].

We also observed that the difference in SCIs was positively associated with dengue transmission. The two Aedes species are usually distinct in terms of ecological environment, feeding preferences and, critically for the SCI, thermal range. Therefore, the high SCI for Ae. aegypti correlating with the low SCI for Ae. albopictus is meaningful, indicating that a larger difference in SCI may indicate less competition between the two mosquito species. In other words, areas that are optimally suitable for both species (i.e. low difference in SCI) may, somewhat counterintuitively, be associated with a decreased dengue transmission risk, potentially through competitive displacement of one species by the other.

There are several limitations to our study. The SCI is estimated on an annual basis, so it may not accurately capture seasonal variation and time-lag effects of climate factors in dengue case occurrence. A finer scale of SCI is thus needed in any subsequent spatiotemporal analysis in this geographical region. Furthermore, climate-based SCI may not accurately reflect dengue outbreak risk in tropical regions where temperature fluctuations have a weaker correlation with DENV transmission. In these tropical endemic countries, dengue outbreaks are most common in areas where multiple DENV serotypes are endemic and heterologous DENV infections are common [56, 57]. In tropical epidemic regions, the SCI should be validated against the locality-specific factor data to increase the predictive efficacy of the SCI. Although mosquito proliferation associated with high temperature and humidity, and the rainy season are major driving factors for dengue outbreaks, many socio-economic and environmental variables that are largely independent of climate, such as urbanisation and overcrowding, also facilitate the spread of dengue. These factors are not considered and accounted for in this study. However, as a first step, our SCI provides valuable insight into the possible relationship between optimal transmission environments and predictions of broad disease occurrence. To improve the prediction efficiency of the SCI, additional parameters, such as environmental variation (e.g. oviposition habitat availability, seasonal and daily temperature fluctuation) and socio-economic variables should be included. An enhanced SCI may be beneficial for extrapolating the possible geographical range of transmission beyond the existing environmental context (e.g. under climate change and for newly invading pathogens).

Conclusions

The results of this study show that the climate-based SCI is beneficial for evaluating dengue outbreak risk and the potential effect of possible vector competition on dengue transmission. We suggest that the DENV may be more transmissible in a region with warmer weather conditions (high SCI for Ae. aegypti). The SCI for Ae. aegypti could be used to predict dengue transmission even in locations lacking Ae. aegypti but with Ae. albopictus present. Moreover, the SCI can be used to investigate the relative contributions of Ae. albopictus and Ae. aegypti to DENV outbreaks in the absence of entomology data in future research.

Availablility data and materials

Data supporting the conclusions of this article are included within the article. In addition, the datasets used and analysed during the study are available from the corresponding author upon reasonable request.

Abbreviations

BIC:

Bayesian information criterion

DENV:

Dengue virus

IRR:

Incidence rate ratio

NB:

Negative binomial

NDVI:

Normalised Difference Vegetation Index

SCI:

Suitable Conditions Index

References

  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:504–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Brady OJ, Gething PW, Bhatt S, Messina JP, Brownstein JS, Hoen AG, et al. Refining the global spatial limits of dengue virus transmission by evidence-based consensus. PLoS Negl Trop Dis. 2012;6:e1760.

    Article  PubMed  PubMed Central  Google Scholar 

  3. WHO. Dengue and severe dengue. 2022. https://www.who.int/news-room/fact-sheets/detail/dengue-and-severe-dengue. Accessed 15 October 2021

  4. Reich NG, Shrestha S, King AA, Rohani P, Lessler J, Kalayanarooj S, et al. Interactions between serotypes of dengue highlight epidemiological impact of cross-immunity. J R Soc Interface. 2013;10:20130414.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Guzman MG, Harris E. Dengue. Lancet. 2015;385:453–65.

    Article  PubMed  Google Scholar 

  6. Guzman MG, Halstead SB, Artsob H, Buchy P, Farrar J, Gubler DJ, et al. Dengue: a continuing global threat. Nat Rev Microbiol. 2010;8:S7–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Simmons CP, Farrar JJ, van Vinh CN, Wills B. Dengue. N Engl J Med. 2012;366:1423–32.

    Article  CAS  PubMed  Google Scholar 

  8. Guzman MG, Gubler DJ, Izquierdo A, Martinez E, Halstead SB. Dengue infection. Nat Rev Dis Primers. 2016;2:16055.

    Article  PubMed  Google Scholar 

  9. Xie H, Zhou H, Yang Y. Advances in the research on the primary dengue vector Aedes aegypti in China. Chin J Vector Biol Control. 2011;22:194–7 (in Chinese).

    Google Scholar 

  10. Dao-fang G, Hong-ning Z. Progress in Dengue fever important vector Aedes albopctus in China. Chin J Vector Biol Control. 2009;20:607–10 (in Chinese).

    Google Scholar 

  11. Wu F, Liu QY, Lu L, Wang JF, Song XP, Ren DS. Distribution of Aedes albopictus (Diptera: Culicidae) in Northwestern China. Vector-Borne Zoonotic Dis. 2011;11:1181–6.

    Article  PubMed  Google Scholar 

  12. Brady OJ, Hay SI. The global expansion of dengue: how Aedes aegypti mosquitoes enabled the first pandemic arbovirus. Annu Rev Entomol. 2020;65:191–208.

    Article  CAS  PubMed  Google Scholar 

  13. Lambrechts L, Scott TW, Gubler DJ. Consequences of the expanding global distribution of Aedes albopictus for dengue virus transmission. PLoS Negl Trop Dis. 2010;4:e646.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Qiu F-X, Chen Q-Q, Ho Q-Y, Chen W-Z, Zhao Z-G, Zhao B-W. The first epidemic of dengue hemorrhagic fever in the People’s Republic of China. Am J Trop Med Hyg. 1991;44:364–70.

    Article  CAS  PubMed  Google Scholar 

  15. Zhang F-C, Zhao H, Li L-H, Jiang T, Hong W-X, Wang J, et al. Severe dengue outbreak in Yunnan, China, 2013. Int J Infect Dis. 2014;27:4–6.

    Article  PubMed  Google Scholar 

  16. Wu J-Y, Lun Z-R, James AA, Chen X-G. Dengue fever in mainland China. Am J Trop Med Hyg. 2010;83:664.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Peng H-J, Lai H-B, Zhang Q-L, Xu B-Y, Zhang H, Liu W-H, et al. A local outbreak of dengue caused by an imported case in Dongguan China. BMC Public Health. 2012;12:1–8.

    Article  Google Scholar 

  18. Cheng Q, Jing Q, Spear RC, Marshall JM, Yang Z, Gong P. Climate and the timing of imported cases as determinants of the dengue outbreak in Guangzhou, 2014: evidence from a mathematical model. PLoS Negl Trop Dis. 2016;10:e0004417.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Barrera R, Amador M, MacKay AJ. Population dynamics of Aedes aegypti and dengue as influenced by weather and human behavior in San Juan. Puerto Rico PLoS Negl Trop Dis. 2011;5:e1378.

    Article  PubMed  Google Scholar 

  20. Campbell KM, Lin C, Iamsirithaworn S, Scott TW. The complex relationship between weather and dengue virus transmission in Thailand. Am J Trop Med Hyg. 2013;89:1066.

    Article  PubMed  PubMed Central  Google Scholar 

  21. World Health Organization & Special Programme for Research and Training in Tropical Diseases & World Health Organization. Dept. of Control of Neglected Tropical Diseases & World Health Organization. Dept. of Epidemic and Pandemic Alert and Response Dengue : guidelines for diagnosis, treatment, prevention and control (New ed). World Health Organization, Geneva, 2009.

  22. Chadee D. Dengue cases and Aedes aegypti indices in Trinidad West Indies. Acta Trop. 2009;112:174–80.

    Article  CAS  PubMed  Google Scholar 

  23. Bowman LR, Runge-Ranzinger S, McCall P. Assessing the relationship between vector indices and dengue transmission: a systematic review of the evidence. PLoS Negl Trop Dis. 2014;8:e2848.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Lin C-H, Wen T-H. Using geographically weighted regression (GWR) to explore spatial varying relationships of immature mosquitoes and human densities with the incidence of dengue. Int J Environ Res Public Health. 2011;8:2798–815.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Mordecai EA, Cohen JM, Evans MV, Gudapati P, Johnson LR, Lippi CA, et al. Detecting the impact of temperature on transmission of Zika, dengue, and chikungunya using mechanistic models. PLoS Negl Trop Dis. 2017;11:e0005568.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Reiter P. Climate change and mosquito-borne disease. Environ Health Perspect. 2001;109:141–61.

    PubMed  PubMed Central  Google Scholar 

  27. Yang H, Macoris M, Galvani K, Andrighetti M, Wanderley D. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiol Infect. 2009;137:1188–202.

    Article  CAS  PubMed  Google Scholar 

  28. Tun-Lin W, Burkot T, Kay B. Effects of temperature and larval diet on development rates and survival of the dengue vector Aedes aegypti in north Queensland Australia. Med Vet Entomol. 2000;14:31–7.

    Article  CAS  PubMed  Google Scholar 

  29. Kamimura K, Matsuse IT, Takahashi H, Komukai J, Fukuda T, Suzuki K, et al. Effect of temperature on the development of Aedes aegypti and Aedes albopictus. Med Entomol Zool. 2002;53:53–8.

    Article  Google Scholar 

  30. Couret J, Dotson E, Benedict MQ. Temperature, larval diet, and density effects on development rate and survival of Aedes aegypti (Diptera: Culicidae). PLoS ONE. 2014;9:e87468.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Mordecai EA, Cohen JM, Evans MV, Gudapati P, Johnson LR, Lippi CA, et al. Detecting the impact of temperature on transmission of Zika, dengue, and chikungunya using mechanistic models. PLoS Negl Trop Dis. 2017;11:5568.

    Article  Google Scholar 

  32. Davis C, Murphy AK, Bambrick H, Devine GJ, Frentiu FD, Yakob L, et al. A regional suitable conditions index to forecast the impact of climate change on dengue vectorial capacity. Environ Res. 2021;195:110849.

    Article  CAS  PubMed  Google Scholar 

  33. National Bureau of Statistics of China. Communiqué of the Seventh National Population Census [1] (No. 3)—population ry Region. 2021. http://www.stats.gov.cn/english/PressRelease/202105/t20210510_1817188.html. Accessed 15 October 2021

  34. Liu K, Hou X, Ren Z, Lowe R, Wang Y, Li R, et al. Climate factors and the East Asian summer monsoon may drive large outbreaks of dengue in China. Environ Res. 2020;183:109190.

    Article  CAS  PubMed  Google Scholar 

  35. Lai S, Huang Z, Zhou H, Anders KL, Perkins TA, Yin W, et al. The changing epidemiology of dengue in China, 1990–2014: a descriptive analysis of 25 years of nationwide surveillance data. BMC Med. 2015;13:1–12.

    Article  CAS  Google Scholar 

  36. Ministry of Health of the People’s Republic of China. Diagnostic criteria for dengue fever (WS 216–2008). Beijing: People’s Medical Publishing House; 2008. p. 1–17 (in Chinese).

    Google Scholar 

  37. Zhengxing W, Chuang L, Alfredo H. From AVHRR-NDVI to MODIS-EVI: advances in vegetation index research. Acta Ecol Sin. 2003;23:979–87.

    Google Scholar 

  38. Didan K, Munoz AB, Solano R, Huete A. MODIS vegetation index user’s guide (MOD13 series). Vegetation Index and Phenology Lab. Tucson: University of Arizona; 2015.

    Google Scholar 

  39. Stata.com. xtnbreg—fixed-effects, random-effects, & population-averaged negative binomial models. https://www.stata.com/manuals/xtxtnbreg.pdf. Accessed 12 May 2021

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

    Article  Google Scholar 

  41. Chung YK, Pang FY. Dengue virus infection rate in field populations of female Aedes aegypti and Aedes albopictus in Singapore. Trop Med Int Health. 2002;7:322–30.

    Article  CAS  PubMed  Google Scholar 

  42. Chow V, Chan Y, Yong R, Lee K, Lim L, Chung Y, et al. Monitoring of dengue viruses in field-caught Aedes aegypti and Aedes albopictus mosquitoes by a type-specific polymerase chain reaction and cycle sequencing. Am J Trop Med Hyg. 1998;58:578–86.

    Article  CAS  PubMed  Google Scholar 

  43. Xu L, Stige LC, Chan K-S, Zhou J, Yang J, Sang S, et al. Climate variation drives dengue dynamics. Proc Natl Acad Sci USA. 2017;114:113–8.

    Article  CAS  PubMed  Google Scholar 

  44. Wesolowski A, Qureshi T, Boni MF, Sundsøy PR, Johansson MA, Rasheed SB, et al. Impact of human mobility on the emergence of dengue epidemics in Pakistan. Proc Natl Acad Sci USA. 2015;112:11887–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Zhu G, Liu T, Xiao J, Zhang B, Song T, Zhang Y, et al. Effects of human mobility, temperature and mosquito control on the spatiotemporal transmission of dengue. Sci Total Environ. 2019;651:969–78.

    Article  CAS  PubMed  Google Scholar 

  46. Wu P-C, Lay J-G, Guo H-R, Lin C-Y, Lung S-C, Su H-J. Higher temperature and urbanization affect the spatial patterns of dengue fever transmission in subtropical Taiwan. Sci Total Environ. 2009;407:2224–33.

    Article  CAS  PubMed  Google Scholar 

  47. Li Y, Kamara F, Zhou G, Puthiyakunnon S, Li C, Liu Y, et al. Urbanization increases Aedes albopictus larval habitats and accelerates mosquito development and survivorship. PLoS Negl Trop Dis. 2014;8:e3301.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Wimberly MC, Davis JK, Evans MV, Hess A, Newberry PM, Solano-Asamoah N, et al. Land cover affects microclimate and temperature suitability for arbovirus transmission in an urban landscape. PLoS Negl Trop Dis. 2020;14:e0008614.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Barrera R, Amador M, Diaz A, Smith J, Munoz-Jordan J, Rosario Y. Unusual productivity of Aedes aegypti in septic tanks and its implications for dengue control. Med Vet Entomol. 2008;22:62–9.

    Article  CAS  PubMed  Google Scholar 

  50. Thongsripong P, Green A, Kittayapong P, Kapan D, Wilcox B, Bennett S. Mosquito vector diversity across habitats in central Thailand endemic for dengue and other arthropod-borne diseases. PLoS Negl Trop Dis. 2013;7:e2507.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Endy TP, Nisalak A, Chunsuttiwat S, Libraty DH, Green S, Rothman AL, et al. Spatial and temporal circulation of dengue virus serotypes: a prospective study of primary school children in Kamphaeng Phet Thailand. Am J Epidemiol. 2002;156:52–9.

    Article  PubMed  Google Scholar 

  52. Yang BY, Borgert BA, Alto BW, Boohene CK, Brew J, Deutsch K, et al. Modelling distributions of Aedes aegypti and Aedes albopictus using climate, host density and interspecies competition. Plos Negl Trop Dis. 2021;15:e0009063.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Tsai P-J, Teng H-J. Role of Aedes aegypti (Linnaeus) and Aedes albopictus (Skuse) in local dengue epidemics in Taiwan. BMC Infect Dis. 2016;16:1–20.

    Article  CAS  Google Scholar 

  54. Chen B-J, Qin B, Bai A-y, Wu J, Deng H, Duan J-h, et al. An experimental study of interspecific competition between Aedes aegypti from Wushi town of Leizhou and Ae. albopictus from different places in Guangdong province, China. Chin J Vector Biol Control. 2020;31:486–9.

  55. Liu B, Gao X, Ma J, Jiao Z, Xiao J, Hayat MA, et al. Modeling the present and future distribution of arbovirus vectors Aedes aegypti and Aedes albopictus under climate change scenarios in Mainland China. Sci Total Environ. 2019;664:203–14.

    Article  CAS  PubMed  Google Scholar 

  56. Nisalak A, Endy TP, Nimmannitya S, Kalayanarooj S, Thisayakorn U, Scott RM, et al. Serotype-specific dengue virus circulation and dengue disease in Bangkok, Thailand from 1973 to 1999. Am J Trop Med Hyg. 2003;68:191–202.

    Article  PubMed  Google Scholar 

  57. Phan DQ, Nguyen LDN, Pham ST, Nguyen T, Pham PTT, Nguyen STH, et al. The distribution of dengue virus serotype in Quang Nam Province (Vietnam) during the outbreak in 2018. Int J Environ Res Public Health. 2022;19:1285.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful to the Chinese Centre for Disease Control and Prevention (CDC) for providing dengue case data in Guangdong. We also thank NASA’s Moderate Resolution Imaging Spectroradiometer (MODIS) for providing the Normalised Difference Vegetation Index (NDVI) data, and the National Bureau of Statistics of China for providing population census data. We would like to thank Jingli Yang for her assistance in preparing the figures.

Funding

This study was supported by the National Health and Medical Research Council, Australia (APP 1138622). The content of this publication is solely the authors' responsibility and does not reflect the views of the Australian National Health and Medical Research Council. The funder had no role in study design, data collection and analysis, decision to publish, or manuscript preparation.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualisation: WH. Methodology: WH, HB, FF, XH, GD. Formal analysis: WH, XL, XH. Writing–original draft preparation: XL. Writing–review & editing: XL, WH, HB, FF, XH, CD, GD. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Wenbiao Hu.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the University Human Research Ethics Committee of Queensland University of Technology (Ref: 1800000058).

Consent for publication

Not applicable.

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.

Supplementary Information

Additional file 1:

Figure S1. Bivariate Moran’s I scatter plot between the SCI for Ae. aegypti (original variable) and the SCI for Ae. albopictus (spatial lag as the second variable) in Guangdong, 2014 (a) and significant determinator (b).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lu, X., Bambrick, H., Frentiu, F.D. et al. Species-specific climate Suitable Conditions Index and dengue transmission in Guangdong, China. Parasites Vectors 15, 342 (2022). https://doi.org/10.1186/s13071-022-05453-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-022-05453-x

Keywords