Climate variation drives dengue dynamics

Lei Xu, Leif C. Stige, Kung-Sik Chan, Jie Zhou, Jun Yang, Shaowei Sang, Ming Wang, Zhicong Yang, Ziqiang Yan, Tong Jiang, Liang Lu, Yujuan Yue, Xiaobo Liu, Hualiang Lin, Jianguo Xu, Qiyong Liu, and Nils Chr. Stenseth
PNAS January 3, 2017 114 (1) 113-118; published ahead of print December 9, 2016 https://doi.org/10.1073/pnas.1618558114

  1. Contributed by Nils C. Stenseth, November 10, 2016 (sent for review September 10, 2016; reviewed by Robert A. Cheke and Kenneth L. Gage)

Significance

Dengue is a vector-borne infectious disease that threatens human health on a global scale. Because of climate change, globalization, and other factors, dengue has increasingly spread to new countries and over larger areas, from tropical to temperate zones. In this study, we found that climate has both direct effects on dengue incidence and indirect effects mediated by mosquito density, as mosquitoes are the vectors of dengue. The quantitative results derived from this study may be helpful toward advancing our understanding of how climate influences vector-borne diseases and prove useful for the control and prevention of dengue fever.

Abstract

Dengue, a viral infection transmitted between people by mosquitoes, is one of the most rapidly spreading diseases in the world. Here, we report the analyses covering 11 y (2005–2015) from the city of Guangzhou in southern China. Using the first 8 y of data to develop an ecologically based model for the dengue system, we reliably predict the following 3 y of dengue dynamics—years with exceptionally extensive dengue outbreaks. We demonstrate that climate conditions, through the effects of rainfall and temperature on mosquito abundance and dengue transmission rate, play key roles in explaining the temporal dynamics of dengue incidence in the human population. Our study thus contributes to a better understanding of dengue dynamics and provides a predictive tool for preventive dengue reduction strategies.

Dengue is one of the most rapidly spreading diseases in the world (1), including within the Guangdong province of southern China (2). During the last 50 y, the incidence of dengue has increased 30-fold with increasing geographic expansion to new countries (1). In 2010, an estimated 390 million dengue infections occurred, of which around 96 million showed symptoms (3). Dengue outbreaks in China were previously thought to be imported and initiated by people traveling to China from dengue-endemic areas elsewhere (4); however, recent studies suggest that dengue may now be endemic to China as well (2). The epidemiological triangle of both dengue fever and dengue hemorrhagic fever, which is the more serious form of dengue, includes hosts (humans), pathogens (one or more of five dengue virus serotypes) (5), and mosquito vectors (Aedes albopictus and Aedes aegypti) with their ecological interactions (6). The dengue outbreaks are qualitatively known to be strongly influenced by temperature (7), humidity, rainfall, and socioeconomic factors like urbanization (8). However, a full understanding of the quantitative nature of such effects is largely lacking. With this paper, we provide such a quantitative understanding of dengue dynamics.

In 2014, an extensive dengue outbreak hit China, with 47,127 dengue cases diagnosed, a new record since 1986 (9). Since the 1990s, dengue epidemics have gradually spread from Guangdong, Hainan, and Guangxi provinces (9). We present here a time series analysis of dengue dynamics, using dengue surveillance data for the years 2005−2015 from Guangzhou, the largest city in Guangdong and the city with the most dengue cases in China. We split the main analysis by using the first 8 y to develop a model, and the three remaining years to test that model, as these latter years encompass exceptionally extensive dengue outbreaks.

Monthly human dengue incidence data (number of diagnosed cases) were obtained from the China National Notifiable Disease Surveillance System (10) (Fig. 1). Monthly surveillance data of A. albopictus density, the only dengue vector species in Guangzhou, were obtained from local Centers for Disease Control and Prevention (11) (Methods). We considered indices of both larval and adult densities. Temperature and precipitation data were obtained from the China Meteorological Data Sharing Service System (www.cma.gov.cn/). We considered several alternative monthly temperature (degrees Celsius) and precipitation indices in the subsequent analyses, including mean temperature, monthly average of daily minimum temperature, monthly average of daily maximum temperature, cumulative precipitation (millimeters), and number of days with rainfall.

Fig. 1.
Fig. 1.

Spatiotemporal dynamics of dengue in mainland China from the beginning of 2005 to the end of 2015. (A) Geographical distribution of dengue incidence in China. Radiuses of circles indicate incidence of dengue. Red circles represent dengue in Guangzhou, black circles represent dengue in other cities in Guangdong province, and gray circles represent dengue in other provinces. (B) Monthly time series of dengue incidence in Guangzhou. (C) Monthly time series of adult population density of the dengue vector A. albopictus. (D) Monthly time series of average maximum temperature in Guangzhou. (E) Monthly time series of number of days with rainfall in Guangzhou.

Results and Discussion

We first explored the direct and indirect effects of temperature and precipitation on dengue outbreaks for all years from 2005−2015 using a structural equation model (SEM) (12). The SEM results suggest significant direct, positive effects of temperature and precipitation on dengue incidence as well as indirect, positive effects of both variables through vector density (Fig. 2). Model diagnostics suggest that the SEM was robust: The normalized χ2 test (i.e., model χ2 divided by the degrees of freedom) = 0.149, root-mean-square error of approximation < 0.001, and comparative fit index = 1 (12). However, due to the presence of many zeros in the data (zero inflation), and the inability of the SEM to capture nonlinear effects, some caution in the interpretation of these results is needed. Nevertheless, this analysis serves as a good stepping stone for the full time series model development using generalized additive models (GAMs) (13) and zero-inflated GAMs (ZIGAMs) (14).

Fig. 2.
Fig. 2.

SEM analysis revealed direct and indirect climate effects on dengue incidence in Guangzhou from 2005 to 2015. Arrows with numbers indicate ecological effects and standardized coefficients. Temperature and precipitation are correlated (correlation coefficient = 0.39). Asterisks indicate statistically significant pathways (P < 0.05).

The generic model formulas for vector and dengue, respectively, are given by Eqs. 1 and 2,ln(Vt+1)=a0+f1[ln(Vt−1+1)]+f2(Tt−1)+f3(Pt−1)+εt.[1]Dt=I{b0+g1[ln(Dt−1+1)]+g2(Vt−1)+g3(Tt−1)+g4(Pt−1)}× exp{b1+g5[ln(Dt−1+1)]+g6(Vt−1)+g7(Tt−1)+g8(Pt−1)+δt}.[2]Here, Vt and Dt indicate, respectively, monthly vector density and dengue incidence in month t. Tt−1 and Pt−1 indicate temperature and precipitation for month t − 1. The parameters a0, b0, and b1 are intercepts. The functions f1, f2, f3, g1, ..., and g8 are either smooth (natural cubic splines) or linear functions. By optimizing model selection criteria that quantified the trade-off between model fit and parsimony (Tables S1 and S2), the temperature, precipitation, and vector indices to include in the models and the functional forms (linear or smooth) were chosen. Vector dynamics were analyzed with a GAM fitted on ln(Vt+1) transformed data. The dengue data (but not the vector data) were zero-inflated (Fig. 1 and Fig. S1). Hence, dengue incidence dynamics were analyzed with a ZIGAM that consisted of a binomial and a lognormal part. The function I(·) in Eq. 2 represents a Bernoulli variable that is either 0 (absent) or 1 (present), whose probability of success is given by the predictor expression enclosed within the parentheses, on the logit scale. The same predictor variables enter into the two submodels of the ZIGAM, although with potentially different functional forms.

View this table:
Table S1.

Model selection results for adult mosquito density

View this table:
Table S2.

Model selection results (log of marginal likelihood, logE) for dengue incidence

Fig. S1.
Fig. S1.

Test of zero inflation. Histogram plots of monthly mosquito density and dengue incidence data.

Results of ZIGAM analyses for the first 8 y (2005−2012) suggested that dengue incidence was best predicted as positive smooth functions of dengue incidence in the previous month, adult mosquito density, monthly average of daily maximum temperature, and number of days with rainfall (Fig. 3 and Table S2). We found that adult mosquito (i.e., vector) density was best predicted as linear effects of adult mosquito density in the previous month and number of days with rainfall (Table 1 and Table S1).

Fig. 3.
Fig. 3.

Analysis of potentially nonlinear influences on dengue incidence in Guangzhou based on data from years 2005−2012, i.e., excluding the three last years of data. To account for zero inflation, (A−D) a separate binomial submodel quantifies predictor effects on outbreak risk (logit scale probability of incidence > 0) and (E−H) a lognormal submodel quantifies predictor effects on outbreak intensity when an outbreak occurs [ln(incidence)].

View this table:
Table 1.

Final model for adult mosquito density

We then used the selected GAM for climate effects on adult mosquito density (Table 1) and the selected ZIGAM for adult mosquito and climate effects on dengue (Fig. 3) to make 1-mo-ahead predictions for 2013−2015. We found that the selected models were indeed able to predict the high peaks in dengue cases for these years (Fig. 4), although with broad prediction bands. The two highest mosquito peaks, in years 2007 and 2014, were not predictable from the model. The 2014 mosquito peak in Guangzhou has recently been associated with the seasonal distribution of precipitation, the daily temperature range, and nonadditive effects of temperature and precipitation (15), pointing to more complex species−environment relationships than accommodated by our model. As our focus was on prediction, we intentionally restricted model complexity to reduce the risk of overfitting.

Fig. 4.
Fig. 4.

Observations and 1-mo-ahead predictions. (A) Adult mosquito density. (B) Dengue incidence. The vertical lines separate the years 2005−2012 from the years 2013−2015 over which out-of-sample predictions were made.

Finally, we repeated model selection, using the whole dataset from 2005 to 2015 (Tables S1 and S2). The same variables were selected in these models, and their estimated effects were similar, indicating that the same biological processes were driving the dynamics in the peak years from 2013 to 2015 as in the earlier period (Fig. 5 and Table 1). No additional factors, for example, related to dengue possibly becoming endemic and potentially occurring earlier in the season (2, 4), were thus needed to explain these peaks in dengue incidence.

Fig. 5.
Fig. 5.

Analysis of potentially nonlinear influences on dengue incidence in Guangzhou based on data from all years 2005−2015. To account for zero inflation, (A−D) a separate binomial submodel quantifies predictor effects on outbreak risk (logit scale probability of incidence > 0) and (E−H) a lognormal submodel quantifies predictor effects on outbreak intensity when an outbreak occurs [ln(incidence)].

This final ZIGAM of dengue dynamics shows that previous-month dengue incidence has positive effects on outbreak risk (i.e., on the probability of the incidence being larger than zero, Fig. 5A) and outbreak intensity [i.e., on ln(incidence), given that incidence is larger than zero, Fig. 5E]. Adult mosquito density mainly seems to affect outbreak intensity, but only after having reached a threshold of around 1 to 1.5 units on the mosquito index scale (monthly averaged A. albopictus captured per trap per night) (Fig. 5F). Temperature has significantly positive effects on both outbreak risk and intensity (Fig. 5 C and G), with some indications of leveling off at high temperatures (Fig. 5G). Precipitation has a nonsignificant positive association with outbreak risk (Fig. 5D) and a significant, approximately linear, and positive association with outbreak intensity (Fig. 5H).

The positive association between temperature and dengue incidence in our study is similar to, for example, Taiwan (16, 17), which is at the same latitude as Guangzhou. Our findings show that this association is mainly related to direct effects of temperature on dengue transmission rate. Mechanisms for such direct effects include, for example, temperature effects on the biting rate of mosquitoes (18), incubation period of pathogens (19), and human exposure to mosquitoes (e.g., by influencing the time spent outdoors or with open windows). In contrast, GAM results show no significant association of the temperature variables considered with mosquito population dynamics. The temperature−mosquito association suggested by the SEM likely reflects the mosquitoes’ seasonal cycle, with higher abundances in summer, as this model analyzes the total abundance rather than the change in abundance from one month to the next.

Our results show that precipitation influences dengue incidence through both direct effects on transmission rate [e.g., through effects on mosquito biting rate (18)], and indirect effects through mosquito population dynamics. Because A. albopictus breed in small pools of water, rainfall could result in an increase of breeding sites (20). We found the number of days with rainfall to be a better predictor than the monthly cumulative precipitation, which suggests that the underlying mechanism may be related to the maintenance of humid conditions over time. These findings complement a recent study showing a positive association between interannual variability in dengue incidence and surface water area (21). Together, our studies show that changes in breeding area of A. albopictus are important drivers of dengue dynamics at both monthly and yearly time scales. However, although precipitation appears to be a key driver of the shorter-term changes (this study), urban landscape features, such as recently constructed artificial lakes in Guangzhou, appear to be key drivers of longer-term trends (21).

As in other dengue epidemic areas (11, 22, 23), we find a strong link between dengue incidence and density of Aedes mosquitoes, the vectors of dengue. These findings support the utility of vector population control through, for example, pesticides, Wolbachia bacteria infection of mosquitoes, or habitat measures, for reducing the intensity of dengue outbreaks (24, 25). Interestingly, such control efforts may be the reason why mosquito density was lower in 2015 than we had predicted based on climate conditions, and why dengue incidence was lower than predicted for 3 mo during the buildup phase of the outbreak that year.

Our study builds on previous research that links dengue and vector dynamics to climate variables (26, 27) and climate change (28), and extends these by quantifying relationships based on state-of-the-art statistical analysis of multiyear time series. Our findings strongly suggest that climate conditions drive the outbreaks of dengue in Guangzhou. However, our modeling approaches have some limitations, such as the low number of years with very high dengue incidence during the study period (2013 and 2014). We did not consider potentially relevant local factors, such as population movements, the occurrence of household water tanks, vegetation cover, and distance to rivers or other water bodies. Further, although our results point to climate factors associated with increased risk of dengue transmission and outbreaks, human activities and their impact on local ecology may be more important than climate factors in driving the long-term trends (29⇓–31). Moreover, the role of climate is likely to vary among locations and even over time in the same location (32). However, the quantitative results derived from this study may be helpful toward advancing our understanding of how climate influences vector-borne diseases and prove useful for the control and prevention of dengue fever.

Methods

Dengue Incidence Data.

Daily records of human dengue cases from 2005 to 2015 were obtained from the China National Notifiable Disease Surveillance System. All human dengue cases were diagnosed according to the diagnostic criteria for Dengue Fever (WS216–2008) enacted by Chinese Ministry of Health (11, 33). The information on individual dengue cases included whether they were local or imported (local cases were defined as dengue cases where infection most likely occurred due to local transmission as there was no documentation of travel to dengue-endemic foreign regions; imported cases were those where infection most likely occurred outside of China). Human dengue incidence data were aggregated per month (denoted as Dt) and matched to monthly surveillance data of adult mosquito population density.

Vector Density Data.

We considered both the adult mosquito density (denoted as Vta) and the larval density (denoted as Vtl). Adult mosquitoes (A. albopictus) were sampled by light traps. Larval density was measured by the Breteau index, which is defined as the number of positive containers per 100 households inspected (1). Both adult density and larval density comprised the information on mosquito density sampled in residential area households (>50 households sampled per month), parks, construction sites, and other urban areas.

Climate Variables in the Models.

We considered several different monthly temperature and precipitation variables for model selection. The temperature variables (degrees Celsius) included average temperature (denoted as Tt−1a), average of daily lowest temperature (denoted as Tt−1l), and average of daily highest temperature (denoted as Tt−1h). Precipitation variables included accumulated precipitation (millimeters, denoted as Pt−1), and number of days with rainfall (days, denoted as Pt−1d). We modeled the effects of climate variables on dengue incidence with t = 1-mo time lag to be consistent with the survival time of A. albopictus under field conditions: 50% of the population has been found to survive ≥10 d, and 5% have been found to survive ≥40 d; the time for A. albopictus transmission of dengue may be less than 40 d (34, 35).

Statistical Modeling.

SEM: Rationale, model construction, and selection.

SEM estimates the structural correlation among the variables (36, 37). SEM was used to analyze the structural linear correlation (as shown in Fig. 2) between temperature (average maximum temperature), precipitation (number of days with rainfall), vector density (monthly averaged mosquitoes captured per trap per night), and dengue incidence (natural log transformed). The effects of temperature and precipitation on dengue were divided into direct effects and effects mediated by vector density. At the same time, the correlation between the variables of temperature and precipitation was assessed in the model. A dummy variable YEAR2014 (year 2014 = 1, other years = 0) was added into the SEM model as a control variable because the large outbreak of dengue in year 2014 (38) otherwise hindered convergence.

GAM of vector dynamics.

The mgcv package (version 1.7-29) of R was used for GAM analyses (39). Because the different temperature and precipitation variables were strongly correlated, we considered all combinations of the alternative temperature, precipitation, and vector density variables (but no more than one variable for each factor in the model). To account for density dependence, we included ln(1 + density) in the previous month as a predictor variable in the model. We also compared linear and nonlinear model formulations, by letting the covariate effects in Eq. 1 to be either linear or smooth (natural cubic splines with a maximum of 3 degrees of freedom). As selection criterion, we used the generalized cross-validation criterion (GCV). A model with lower GCV has more predictive power and was hence selected (Table S1). To evaluate an additional conservative criterion, we used genuine cross-validation. Specifically, we removed one “test” year at a time, fitted the model on the remaining dataset, made predictions for the test year, repeated this for all years, and calculated log-scale root-mean-square prediction error [i.e., the square root of the sum of squared differences ln(observed density + 1) – ln(predicted density + 1)]. Deviance residuals from the final modes were checked for temporal autocorrelation using autocorrelation and partial autocorrelation functions.

ZIGAM of Dengue dynamics.

Histograms of dengue incidence clearly showed zero inflation problems (Fig. S1). The COZIGAM package (version 2.0-4) of R was used for Constrained and Unconstrained Zero-Inflated GAMs (14). The constrained (CoZIGAM) and the unconstrained (UnZIGAM) models differ in that CoZIGAM constrains the fitted values of the binomial part of the model to be proportional to the fitted values of the lognormal part (on the linear scales of the predictors). As model selection criterion, we used an approximation of the log of the marginal likelihood (logE), which is related to the Bayesian Information Criterion. Models with higher logE are preferred. As response variable, only the incidence of locally infected dengue cases was considered. For the autoregressive term, ln(1 + incidence) in the previous month, both local and imported dengue cases were considered.

The maximum degree of freedom for the smooth terms was set to 3 in the model selection. By inspecting the model providing the highest logE for the first 8 y (Table S2), we found that the effect of average maximum temperature was not biologically interpretable (Fig. S2). Most likely, this result reflects overfitting, and hence we modified the model by reducing the maximum degree of freedom for this term to 2. This rendered an interpretable temperature effect (Fig. 3). The logE of this model (−60.515) remained a better fit than those of other candidate models from Table S2. Hence, this model was used for predictions.

Fig. S2.
Fig. S2.

The initially selected ZIGAM of dengue incidence based on data from 2005 to 2012. To account for zero inflation, (A–D) a separate binomial submodel quantifies predictor effects on outbreak risk (logit scale probability of incidence > 0) and (E–H) a lognormal submodel quantifies predictor effects on outbreak intensity when an outbreak occurs [ln(incidence)]. The effect of average maximum temperature in G is likely overfitted.

Predictions of Dengue incidence for 2013–2015.

The selected vector and dengue models based on data from 2005 to 2012 were used to calculate 1-mo-ahead predictions for 2013–2015. Predictions of mosquito density and dengue incidence for month t were calculated based on observed temperature, precipitation, dengue incidence, and mosquito density in month t − 1. The predicted mosquito density for month t was used as input in the prediction of dengue for month t.

Uncertainty estimation.

Prediction bands were calculated using bootstrap, as follows: For vector dynamics, we used a wild bootstrap approach (40). The residuals from the model for the reduced dataset (2005−2012) were first rescaled to have the same variance as estimated by the model. We did this by multiplying the residuals with the square root of the scale factor estimated from the model, divided by the standard deviation of the residuals. Ten thousand bootstrap datasets of the same size as the original reduced dataset were generated by adding to the fitted values their corresponding (rescaled) residuals, with the sign of the residuals randomly flipped. Specifically, each residual was multiplied by a random variable that equaled 1 or −1 with equal probability, before being added back to the fitted values. A bootstrap distribution for vector density predictions was generated by refitting the vector model on each bootstrap dataset and using these refitted models to make predictions for the whole period (2005−2015).

For dengue dynamics, we generated 10,000 bootstrap datasets in a two-step approach. In the first step, we sampled with replacement the rescaled residuals from the lognormal part of the ZIGAM and added them back to the fitted values. In the second step, each data point was multiplied by 0 or 1, with probability p being the fitted value from the binomial part of the ZIGAM. We then refitted the ZIGAM with the bootstrap data. Finally, we used these refitted ZIGAMs and the bootstrap time series of mosquito density to generate a bootstrap distribution of dengue incidence for the whole period.

The 0.025 [0.05] and 0.975 [0.95] quantiles from this distribution delineate the 95% [90%] prediction limits.

Acknowledgments

We thank Sari C. Cunningham for improving the English of this manuscript. This research was supported by the National Basic Research Program of China (973 Program) (Grant 2012CB955504), the National Natural Science Foundation of China (Grant 81273139), State Key Laboratory for Infectious Disease Prevention and Control Independent Fund (Grant 2014SKLID106), and the Centre for Ecological and Evolutionary Synthesis of the University of Oslo.

Footnotes

  • ↵1L.X. and L.C.S. contributed equally to this work.

  • ↵2To whom correspondence may be addressed. Email: n.c.stenseth{at}ibv.uio.no, xujianguo{at}icdc.cn, or liuqiyong{at}icdc.cn.

  • Author contributions: L.X. and N.C.S. designed research; L.X., L.C.S., K.-S.C., J.Z., J.Y., S.S., M.W., Z. Yang, Z. Yan, T.J., L.L., Y.Y., X.L., H.L., J.X., Q.L., and N.C.S. performed research; L.X., L.C.S., K.-S.C., and J.Z. analyzed data; and L.X., L.C.S., and N.C.S. wrote the paper.

  • Reviewers: R.A.C., University of Greenwich at Medway; and K.L.G., Centers for Disease Control and Prevention.

  • The authors declare no conflict of interest.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1618558114/-/DCSupplemental.

Freely available online through the PNAS open access option.

References

  1. ↵
    1. World Health Organization
    (2009) Dengue Guidelines for Diagnosis, Treatment, Prevention and Control: New Edition (World Health Org, Geneva), pp 3–13.
    .
  2. ↵
    1. Lin YP, et al.
    (2016) Clinical and epidemiological features of the 2014 large-scale dengue outbreak in Guangzhou city, China. BMC Infect Dis 16(1):102.
    .
    OpenUrl
  3. ↵
    1. Bhatt S, et al.
    (2013) The global distribution and burden of dengue. Nature 496(7446):504–507.
    .
    OpenUrlCrossRefPubMed
  4. ↵
    1. Lai S, et al.
    (2015) The changing epidemiology of dengue in China, 1990−2014: A descriptive analysis of 25 years of nationwide surveillance data. BMC Med 13(1):100.
    .
    OpenUrlCrossRefPubMed
  5. ↵
    1. Normile D
    (2013) Tropical medicine. Surprising new dengue virus throws a spanner in disease control efforts. Science 342(6157):415.
    .
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Rezza G
    (2012) Aedes albopictus and the reemergence of Dengue. BMC Public Health 12(1):72.
    .
    OpenUrlCrossRefPubMed
  7. ↵
    1. Yang HM,
    2. Macoris MLG,
    3. Galvani KC,
    4. Andrighetti MTM,
    5. Wanderley DMV
    (2009) Assessing the effects of temperature on dengue transmission. Epidemiol Infect 137(8):1179–1187.
    .
    OpenUrlCrossRefPubMed
  8. ↵
    1. Naish S, et al.
    (2014) Climate change and dengue: A critical and systematic review of quantitative modelling approaches. BMC Infect Dis 14(1):167.
    .
    OpenUrlCrossRefPubMed
  9. ↵
    1. Wu JY,
    2. Lun ZR,
    3. James AA,
    4. Chen XG
    (2010) Dengue fever in mainland China. Am J Trop Med Hyg 83(3):664–671.
    .
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Zhang H, et al.
    (2014) Spatiotemporal characterizations of dengue virus in mainland China: Insights into the whole genome from 1978 to 2011. PLoS One 9(2):e87630.
    .
    OpenUrl
  11. ↵
    1. Sang S, et al.
    (2014) Predicting local dengue transmission in Guangzhou, China, through the influence of imported cases, mosquito density and climate variability. PLoS One 9(7):e102755.
    .
    OpenUrl
  12. ↵
    1. Xu L, et al.
    (2015) The trophic responses of two different rodent-vector-plague systems to climate change. Proc R Soc B 282(1800):20141846.
    .
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Hastie T,
    2. Tibshirani R
    (1990) Generalized Additive Models (Chapman and Hall, London), 1st Ed.
    .
  14. ↵
    1. Liu H,
    2. Chan K-S
    (2010) Introducing COZIGAM: An R package for unconstrained and constrained zero-inflated generalized additive model analysis. J Stat Softw 35(11):26.
    .
    OpenUrl
  15. ↵
    1. Wang X,
    2. Tang S,
    3. Cheke RA
    (2016) A stage structured mosquito model incorporating effects of precipitation and daily temperature fluctuations. J Theor Biol 411:27–36.
    .
    OpenUrl
  16. ↵
    1. Wu PC,
    2. Guo HR,
    3. Lung SC,
    4. Lin CY,
    5. Su HJ
    (2007) Weather as an effective predictor for occurrence of dengue fever in Taiwan. Acta Trop 103(1):50–57.
    .
    OpenUrlCrossRefPubMed
  17. ↵
    1. Wu PC,
    2. Su HJ,
    3. Guo HR,
    4. Lung SC,
    5. Lin CY
    (2005) Temperature can be an effective predictor for dengue fever outbreak. Epidemiology 16(5):S72.
    .
    OpenUrl
  18. ↵
    1. Patz JA,
    2. Epstein PR,
    3. Burke TA,
    4. Balbus JM
    (1996) Global climate change and emerging infectious diseases. JAMA 275(3):217–223.
    .
    OpenUrlCrossRefPubMed
  19. ↵
    1. Xiao FZ, et al.
    (2014) The effect of temperature on the extrinsic incubation period and infection rate of dengue virus serotype 2 infection in Aedes albopictus. Arch Virol 159(11):3053–3057.
    .
    OpenUrl
  20. ↵
    1. Gage KL,
    2. Burkot TR,
    3. Eisen RJ,
    4. Hayes EB
    (2008) Climate and vectorborne diseases. Am J Prev Med 35(5):436–450.
    .
    OpenUrlCrossRefPubMed
  21. ↵
    1. Tian H, et al.
    (2016) Surface water areas significantly impacted 2014 dengue outbreaks in Guangzhou, China. Environ Res 150:299–305.
    .
    OpenUrl
  22. ↵
    1. Morin CW,
    2. Comrie AC,
    3. Ernst K
    (2013) Climate and dengue transmission: Evidence and implications. Environ Health Perspect 121(11-12):1264–1272.
    .
    OpenUrlCrossRefPubMed
  23. ↵
    1. Tapia-Conyer R,
    2. Betancourt-Cravioto M,
    3. Méndez-Galván J
    (2012) Dengue: An escalating public health problem in Latin America. Paediatr Int Child Health 32(Suppl 1):14–17.
    .
    OpenUrlPubMed
  24. ↵
    1. Hoffmann AA, et al.
    (2011) Successful establishment of Wolbachia in Aedes populations to suppress dengue transmission. Nature 476(7361):454–457.
    .
    OpenUrlCrossRefPubMed
  25. ↵
    1. Maciel-de-Freitas R,
    2. Valle D
    (2014) Challenges encountered using standard vector control measures for dengue in Boa Vista, Brazil. Bull World Health Organ 92(9):685–689.
    .
    OpenUrl
  26. ↵
    1. Hugo LE, et al.
    (2014) Adult survivorship of the dengue mosquito Aedes aegypti varies seasonally in central Vietnam. PloS Negl Trop Dis 8(2):e2669.
    .
    OpenUrlCrossRefPubMed
  27. ↵
    1. Otero M,
    2. Solari HG
    (2010) Stochastic eco-epidemiological model of dengue disease transmission by Aedes aegypti mosquito. Math Biosci 223(1):32–46.
    .
    OpenUrlCrossRefPubMed
  28. ↵
    1. Murray NEA,
    2. Quam MB,
    3. Wilder-Smith A
    (2013) Epidemiology of dengue: Past, present and future prospects. Clin Epidemiol 5(1):299–309.
    .
    OpenUrl
  29. ↵
    1. Gubler D
    (2012) Dengue, urbanization and globalization: The unholy trinity of the 21st century. Int J Infect Dis 16:E2.
    .
    OpenUrl
  30. ↵
    1. Ooi EE,
    2. Gubler DJ
    (2009) Dengue in Southeast Asia: Epidemiological characteristics and strategic challenges in disease prevention. Cad Saude Publica 25(Suppl 1):S115–S124.
    .
    OpenUrlCrossRefPubMed
  31. ↵
    1. Reiter P
    (2001) Climate change and mosquito-borne disease. Environ Health Perspect 109(Suppl 1):141–161.
    .
    OpenUrlCrossRefPubMed
  32. ↵
    1. Cazelles B,
    2. Chavez M,
    3. McMichael AJ,
    4. Hales S
    (2005) Nonstationary influence of El Niño on the synchronous dengue epidemics in Thailand. PLoS Med 2(4):e106.
    .
    OpenUrlCrossRefPubMed
  33. ↵
    1. Sun J, et al.
    (2011) Dengue virus serotype 3 subtype III, Zhejiang Province, China. Emerg Infect Dis 17(2):321–323.
    .
    OpenUrlPubMed
  34. ↵
    1. Brady OJ, et al.
    (2013) Modelling adult Aedes aegypti and Aedes albopictus survival at different temperatures in laboratory and field settings. Parasites Vectors 6(1):351.
    .
    OpenUrlCrossRefPubMed
  35. ↵
    1. Costanzo KS,
    2. Schelble S,
    3. Jerz K,
    4. Keenan M
    (2015) The effect of photoperiod on life history and blood-feeding activity in Aedes albopictus and Aedes aegypti (Diptera: Culicidae). J Vector Ecol 40(1):164–171.
    .
    OpenUrlCrossRefPubMed
  36. ↵
    1. Bentler PM
    (1990) Comparative fit indexes in structural models. Psychol Bull 107(2):238–246.
    .
    OpenUrlCrossRefPubMed
  37. ↵
    1. Steiger JH
    (1990) Structural model evaluation and modification: An interval estimation approach. Multivariate Behav Res 25(2):173–180.
    .
    OpenUrlCrossRefPubMed
  38. ↵
    1. Maassen GH,
    2. Bakker AB
    (2001) Suppressor variables in path models—Definitions and interpretations. Sociol Methods Res 30(2):241–270.
    .
    OpenUrlCrossRef
  39. ↵
    1. R Core Team
    (2014) R: A Language and Environment for Statistical Computing (R Found Stat Comput, Vienna).
    .
  40. ↵
    1. Mammen E
    (1993) Bootstrap and wild bootstrap for high-dimensional linear-models. Ann Stat 21(1):255–285.
    .
    OpenUrlCrossRef
View Abstract
Back to top
Article Alerts
Email Article
Citation Tools
Request Permissions
Climate variation drives dengue dynamics
Lei Xu, Leif C. Stige, Kung-Sik Chan, Jie Zhou, Jun Yang, Shaowei Sang, Ming Wang, Zhicong Yang, Ziqiang Yan, Tong Jiang, Liang Lu, Yujuan Yue, Xiaobo Liu, Hualiang Lin, Jianguo Xu, Qiyong Liu, Nils Chr. Stenseth
Proceedings of the National Academy of Sciences Jan 2017, 114 (1) 113-118; DOI: 10.1073/pnas.1618558114
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo

You May Also be Interested in

A study uncovers genetic evidence of early colonization of Jamaica by primates and suggests that adaptation can shape species morphology in novel environments, even when species come from morphologically conservative families. Image courtesy of Lorraine Meeker (American Museum of Natural History, New York).
Primate colonization of the Caribbean
A study uncovers genetic evidence of early colonization of Jamaica by primates and suggests that adaptation can shape species morphology in novel environments, even when species come from morphologically conservative families.
Image courtesy of Lorraine Meeker (American Museum of Natural History, New York).
Researchers report that some moth species use sound absorbing-scales to avoid detection by bats and suggest that the findings could aid the design of biologically inspired sound absorbers for an array of applications, including noise mitigation and building acoustics. Image courtesy of Thomas R. Neil.
How moth scales help avoid detection by bats
Researchers report that some moth species use sound absorbing-scales to avoid detection by bats and suggest that the findings could aid the design of biologically inspired sound absorbers for an array of applications, including noise mitigation and building acoustics.
Image courtesy of Thomas R. Neil.
Proceedings of the National Academy of Sciences: 115 (50)
Current Issue

Submit

Jump to section