Skip to main content
  • Research article
  • Open access
  • Published:

Short-term effects of meteorological factors on pediatric hand, foot, and mouth disease in Guangdong, China: a multi-city time-series analysis

Abstract

Background

Literature shows inconsistency in meteorological effects on Hand, foot, and mouth disease (HFMD) in different cities. This multi-city study aims to investigate the meteorological effects on pediatric HFMD occurrences and the potential effect modification by geographic factors.

Methods

Based on daily time-series data in eight major cities in Guangdong, China during 2009–2013, mixed generalized additive models were employed to estimate city-specific meteorological effects on pediatric HFMD. Then, a random-effect multivariate meta-analysis was conducted to obtain the pooled risks and to explore heterogeneity explained by city-level factors.

Results

There were a total of 400,408 pediatric HFMD cases (children aged 0–14 years old) with an annual incidence rate of 16.6 cases per 1,000 children, clustered in males and children under 3 years old. Daily average temperature was positively associated with pediatric HFMD cases with the highest pooled relative risk (RR) of 1.52 (95 % CI: 1.30–1.77) at the 95th percentile of temperature (30.5 °C) as compared to the median temperature (23.5 °C). Significant non-linear positive effects of high relative humidity were also observed with a 13 % increase (RR = 1.13, 95 % CI: 1.00–1.28) in the risk of HFMD at the 99th percentile of relative humidity (86.9 %) as compared to the median value (78 %). The effect estimates showed geographic variations among the cities which was significantly associated with city’s latitude and longitude with an explained heterogeneity of 32 %.

Conclusions

Daily average temperature and relative humidity had non-linear and delayed effects on pediatric HFMD and the effects varied across different cities. These findings provide important evidence for comprehensive understanding of the climatic effects on pediatric HFMD and for the authority to take targeted interventions and measures to control the occurrence and transmission of HFMD.

Peer Review reports

Background

Hand, foot, and mouth disease (HFMD) is a common viral illness that usually affects infants and children under five, characterized with typical symptoms of fever, skin eruptions on hands and feet, and vesicles in oropharynx [1]. Because of non-effective therapy or vaccine, there is wide transition and periodical outbreak. Most cases are mild and self-limited, however, rare but severe complications or even deaths may occur. HFMD burdens the entire society and constitutes a worldwide public health threat. Particularly in recent years, the morbidity of HFMD keeps growing. In mainland China, HFMD is monitored as a class “C” notifiable disease with 2,712,925 reported cases and 384 deaths in 2014 [2]. Guangdong, the biggest province in South China, has suffered frequently from HFMD epidemics, because of the large population density and mobility and the typical tropical or sub-tropical climate.

Previous literature has shown that HFMD morbidity has a 1-year periodicity, with a major peak between spring and early summer and a smaller peak during autumn [35]. Some single-city studies have documented that meteorological factors play an important role in the seasonal pattern of HFMD, but there are inconsistencies in the identified factors, the exposure-response relationships, the lag patterns and the time scales. For example, Chang et al [6] found HFMD cases increased with daily mean temperature in Taiwan, while such effects were not found in island-type territory in East Asia [7]; Wu et al [8] pointed out a non-linear relationship between humidity and HFMD cases in Rizhao, China, but a positive linear association was identified in Taiwan [6]; The maximum lag of delayed meteorological effects was identified from several days to 3 months [8, 9].

Previous assessments of meteorological effects on HFMD occurrences were mainly based on an Autoregressive Integrated Moving Average (ARIMA) model and Seasonal Autoregressive Integrated Moving Average (SARIMA) model, which clarify the effects of meteorological factors at temporal scale with the hypothesis of linear association [9, 10]. It is required to further explicate potential lag patterns and nonlinear features of meteorological effects. Recently, a few studies have applied generalized additive model (GAM) or distributed lag non-linear models (DLNM) to quantify the temporal effects of meteorological factors on HFMD [10, 11]. In addition to temporal variations, there were substantial geographic differences in HFMD occurrences, even in the same climate [12], however, the influence of geographic factors on the climate-HFMD association remains unknown.

This study aims to investigate the meteorological effects on pediatric HFMD occurrences, and to examine potential effect modification by geographic factors in eight major cities in Guangdong, China.

Methods

Data sources

In mainland China, the web-based infection diseases monitor information system has been established in 2004, and HFMD has been monitored and reported as a class C notifiable disease since 2008. This system was managed and maintained by the Chinese Centre for Disease Control and Prevention (China CDC). China CDC provided all pediatric HFMD cases (patients aged 0–14 years old) in eight major cities of Guangdong Province from 1 January 2009 to 31 December, 2013. The eight cities included Guangzhou, Shaoguan, Shantou, Heyuan, Yangjiang, Guangning, Luoding, and Xuwen (Fig. 1).

Fig. 1
figure 1

Map of Guangdong and annual morbidity of pediatric HFMD (per 1,000 children) in eight cities. The stars represent the locations of the cities

Daily meteorological data during the study period for each city were obtained from the China Meteorological Data Sharing Service System, including daily average temperature, precipitation, atmospheric pressure, relative humidity, wind speed and sunshine hours. Demographical and socio-economic data in each city, including population density (persons per square km), the ratio of male to female population (%), the ratio of children under 5 years old (%), Gross Domestic Product (GDP) per capita (RMB), living space per capita (square meter), and the ratio of population flow (%) were collected from the Sixth National Population Census of China in 2010.

Statistical analyses

A two-stage analysis was performed to quantify the meteorological effects on HFMD occurrences. In the first stage, Spearman rank correlation analysis was conducted to check the collinearity between meteorological factors. There was a strong negative correlation between atmospheric pressure and average temperature with a Spearman correlation coefficient ranging from -0.885 to -0.821 in eight cities, and precipitation was highly correlated with relative humidity (Additional file 1). Considering no evidence of the effects of atmospheric pressure and precipitation on HFMD, and avoiding a collinearity problem, we did not include atmospheric pressure and precipitation in the final model.

A mixed generalized additive model (MGAM) was used for each of eight cities to derive parameter estimates of the exposure–response associations between meteorological factors and HFDM. In the second stage, these city-specific estimates were pooled as outcomes using a random-effect multivariate meta-analytical model. We examined the potential heterogeneity between cities and identified city-level factors associated with the heterogeneity.

First-stage model

As an extension of GAM, MGAM incorporates additive parametric functions of covariates and autoregressive terms into models [13]. MGAM has very good performance to control for the autocorrelation of residuals and lead to a robust estimate of the estimates and standard errors. We firstly examined the effects of four meteorological measures and did not find statistically significant effects of average wind speed and sunshine hours on HFMD (Additional file 2). To avoid the redundancy or the inappropriate parameters estimates, these two factors were excluded from the final model. That is, we only considered the effects of mean temperature and relative humidity in the final model. The first-stage model can be specified as follows:

$$ \ln \left(\mathrm{E}\left({y}_{it}\right)\right) = {l}_{it} + {\tau}_{it} $$
(1)

fixed effects:

$$ \begin{array}{l}\begin{array}{l}{l}_{it}=f\left(tem{p}_{it},4\right)+f\left( hum{i}_{it},4\right)\hfill \\ {}\kern3em + ns\left(tim{e}_i,5*6\right)+I\left( Holida{y}_{it}\right)+I\left( DO{W}_{it}\right)\hfill \end{array}\\ {}\ \end{array} $$
(2)

autoregressive terms:

$$ {\tau}_{it}={\displaystyle \sum_{\mathrm{k}=1}^3}{\beta}_{ik}\left( \ln \left({\mathrm{y}}_{i,t-k}^{*}\right)-{l}_{i,t-k}\right) $$
(3)

where y* i,t-k  = max (y i,t-k , 0.5) and β ik is the coefficient of autoregressive random effect. y it is the daily number of HFMD cases in city i (1,2, …, 8) on calendar day t (1, 2, …,1826), supposed to follow the quasi-Poisson distribution allowing for over-dispersion. Function f Sj represents the two-dimensional space of cross-basis functions of natural cubic splines with unified lag to determine their distributed lag non-linear effects [14, 15]. Temp it and humi it indicate daily average temperature and relative humidity with 4° of freedom (df). Consist with previous study, we specified a maximum lag of 14 days [16] to adequately examine the lagged effects. A natural cubic spline of time with 6 df per year was applied to control for he seasonality and long-term trend of HFMD occurrences. Holiday and DOW are the categorical variables indicating the public holiday and day of the week, respectively. The relative risk (RR) of HFMD occurrences was calculated as compared to the median values of meteorological factors.

To further examine whether the meteorological effects on HFMD differed in two peak seasons (i.e. spring and autumn), we conducted stratified analyses for spring (March-May) and autumn (September-November). Sensitivity analyses were performed by varying df (3–5) for lags, df (3–5) for meteorological factors, df (4–7) for time as well as the maximum lag (14–16).

Second-stage model

The second-stage model can be specified as follows [17]:

$$ {\theta}_i^{\hbox{'}}\sim {\mathrm{N}}_k\left(\theta, {S}_i+\varPsi \right) $$
(4)

where θ is the coefficients of the function f. θ i is the valid parameter in city i. θ i is the estimate of θ i as the outcome parameters acquired in the first stage and follows the multivariate normal distribution Nk(θ i , S i ) with k dimension. Ψ represents the between-city variance-covariance matrix. The restricted maximum likelihood (REML) was utilized to estimate the parameters in the random-effect multivariate meta-analytical model. The heterogeneity was identified by the Cochran Q test and quantified by I2. Best linear unbiased prediction (BLUP) was applied to project the exposure–response associations [18]. Wald test was used to identify the statistical significance of the associations between the city-level variables and the outcome parameters. All data analyses were performed in R 3.1.2.

Results

During the 5-year period from 1 January, 2009 to 31 December, 2013, there were a total of 400,408 HFMD cases, with an annual morbidity of 16.6 cases per 1,000 children. The majority of cases (95.3 %) were diagnosed by clinical symptoms and signs. The most common causative agent are coxsackievirus A16 (CV-A16) (4855 cases, 25.8 %) and enterovirus71 (EV71) (4407 cases, 23.4 %). 397 cases (0.1 %) experienced severe situation among which a considerable proportion of patients (14.4 %) died. The annual morbidity among male was approximately 1.5 times as high as that among female. Children under 3 years old were at the highest risk of HFMD with an annual morbidity of 58.0 cases per 1,000 children. Most cases (80.11 %) were scattered children who did not attend school or kindergarten, and among them 93.01 % were under 3 years old. The summarized statistics of HFMD cases and meteorological factors in each city were shown in Additional file 3.

Figure 2 shows the exposure-response associations between meteorological factors and HFMD cases in each of eight cities. Although there were some differences in the magnitude of effect estimates (i.e. RR) among eight cities (Fig. 2a), in general, daily average temperature was positively associated with HFMD. The highest pooled cumulative RR was 1.52 (95 % confidence interval (CI): 1.30–1.77) over lag 0–14 days at the 95th percentile (30.5 °C) as compared to the median temperature (23.5 °C). The effects of high temperatures peaked at lag 6 days and nearly disappeared at lag 14 days (Fig. 3-A3 and 3-A4). There was 37 % increase (RR = 1.37, 95 % CI: 1.18–1.59) in HFMD cases over lag 0–14 days by comparing the 75th percentile of temperature to the median temperature, and significant cumulative effect estimates were also observed for lag 0–7 days. The risk of HFMD increased 19 % (RR = 1.19, 95 % CI: 1.07–1.31) and 51 % (RR = 1.51, 95 % CI: 1.29–1.77) over lag 0–7 days and 0–14 days, respectively, by comparing the 95th percentile to the median temperature (Table 1).

Fig. 2
figure 2

The cumulative effects of meteorological factors on HFMD over lag 0–14 days in Guangdong, 2009–2013. a and b show the cumulative relative risks of HFMD associated with temperature and relative humidity, respectively. The bold red line represents the pooled effects, and the dashed lines represent the city-specific estimates. Reference values were the medians, that is, 23.5 °C for temperature and 78.0 % for relative humidity

Fig. 3
figure 3

The lag patterns of pooled effects of meteorological factors on HFMD across lag 0–14 days. A1-A4 show the pooled relative risks (RRs) of HFMD associated with the 5th, 25th, 75th, and 95th percentiles of average temperature compared to the median temperature, respectively. B1-B4 show the corresponding effects of relative humidity

Table 1 The pooled effects of meteorological factors on HFMD occurrences in eight cities

The cumulative RR of HFMD increased with the relative humidity and peaked at a relative humidity of the 99th percentile of humidity (86.9 %) with a RR of 1.13 (95 % CI: 1.00–1.28) compared to the median humidity (Fig. 2b). The lag-response curve shows that low humidity was associated with a decrease in HFMD occurrences (Fig. 3-B1 and 3-B2). High relative humidity was associated with an increase in HFMD and the effects grew during the lag period of 0–6 days (Fig. 3-B3 and 3-B4). An increase of 3 % (RR = 1.03, 95 % CI: 1.00–1.06) and 12 % (RR = 1.12, 95 % CI:1.02–1.24) in HFMD cases was associated with the 75th percentile compared to the median humidity over lag 0–7 days and lag 0–14 days, respectively (Table 1).

The HFMD incidences showed slightly differences with day of the week. An increase of 3.2 % (RR = 1.032, 95 % CI: 1.00, 1.065) was found in HFMD cases on Monday as compared to Sunday. The RRs (95 % CI) were 0.996 (0.96–1.025), 0.977 (0.957–0.998), 0.961 (0.923–0.999), 0.975 (0.926–1.027), and 0.96 (0.936–0.984) on Tuesday-Saturday, respectively, as compared to Sunday. HFMD incidence tended to be higher on holidays (RR = 1.02, 95 % CI: 0.960–1.084) than non-holidays.

A significant heterogeneity (Q = 75.93, p < 0.001) was identified between the eight cities in the intercept-only random-effect multivariate model, which contributed a moderate part (63.12 %) to the overall variations (Table 2). Multivariate meta-regression analyses revealed that the heterogeneity was associated with city’s latitude and longitude, while no significant association was found with other factors (P > 0.05) (Table 2). The inclusion of these two factors in the model showed a substantial decrease in heterogeneity (31.38 %), compared to the intercept-only model.

Table 2 Heterogeneity of meteorological effects and its association with city-level factors

The stratified analyses revealed that HFMD risk increased with temperature and relative humidity both in spring and autumn, and there were generally not substantial differences in the effect estimates (Additional file 4). Sensitivity analyses on df for meteorological factors, lag, and time, and maximum lag demonstrated that the final model was appropriate and stable without substantial variations in Likelihood Akaike information criteria for quasi-Poison (Q-AIC) value and the effect estimates.

Discussion

We observed the annual morbidity of pediatric HFMD was 16.6 cases per 1,000 children in eight cities in Guangdong, China. The HFMD morbidity in the whole population in Guangdong (2.2 cases per 1,000 persons) was approximately 1.6 times higher than the national average in China (1.4 cases per 1,000 persons) in 2010 and 2.8 times than that reported in Japan in 2010 (0.8 cases per 1,000 persons) [3, 19, 20]. The morbidity of HFMD in males was approximately 1.5 times as high as that among females, consistent with the reports in Korea and Japan [21, 22]. A possible explanation could be gender difference in susceptibility at host genetic level [23]. Also, children under three years old were more susceptible to HFMD, which was observed in other regions [24, 25]. Their poor immune systems that are lack of neutralizing antibodies can also contribute to high susceptibility to HFMD infection [26]. By contrast, a study in Taiwan found that infant had the least HFMD morbidity probably because of the preexisting neutralizing antibody to EV71, acquired through transplacental transfer [27].

In this study, we observed non-linear and delayed effects of daily average temperature, and relative humidity on HFMD cases, and such effects varied geographically. There is no doubt that temperature, especially the average temperature, is the generally acknowledged meteorological factor affecting HFMD occurrences [10, 11, 19, 28]. Several researches considered maximum or minimum temperature as an indicator [10, 19, 29], but average temperature could be the most familiar and powerful predictor to the public. The positive and significant non-liner effects of average temperature on HFMD cases were found in this study, supported by preceding findings [10, 19, 22, 29]. The effects peaked at 30.5 °C (the 98th percentile of temperature), which was similar to the finding in Japan (about 29 °C, the 95th percentile of temperature) [22]. This is reasonable according to an experimental finding in cynomolgus monkeys that the infectivity and activity of EV71 could be inhibited when the temperature is higher than 25 °C due to the potential properties of circulating strains [30, 31]. Besides, a vitro experiment noted that enterovirus replication was restricted approximately 90 % at 39 °C, as compared to the replication at 37 °C [32]. Additionally, the serological antibody in human body may be adjusted by temperature. An experimental study in mice has shown that the antibody to EV71 increases with temperature by accelerating the DNA replication [33], suggesting that high ambient temperature may limit HFMD occurrences. However, the mechanism still needs to be explored in future researches.

Relative humidity was another significant factor positively affecting the HFMD occurrences and the associated risk did not decrease until it was 86.9 % (the 99th percentile of humidity). The potential reason for the effects of relative humidity is possibly due to its profound effects on immunity-oriented problems. The metabolism rate of children decrease when the relatively humidity is relatively high, facilitating the HFMD infection. Besides, previous experiment found that under high relative humidity, enterovirus 70 can be recovered from non-porous surfaces even after 24 h and high humid condition was important for the transmission of enterovirus 70 [34]. We observed high relative humidity effect delayed 2 days and peaked at about 6 days later, which was supported by the possible period of HFMD incubation (2–7 days). A longer lag up to 3 weeks was recognized in Japan despite its insignificant effect [22]. Limited literature has demonstrated the potential explanation of this discrepancy in lag patterns. The longer lag days seemed to be irrational considering the incubation period.

In this study, we found that Sunday and Monday had relatively higher HFMD incidence than the rest days of the week. The HFMD incidence tended to be higher on holidays than on non-holidays. Ma et al [13] pointed out a similar trend in bacillary dysentery incidence. This can be explained by patients’ behaviors. On public holidays, children have more outdoor activities and cluster in public places, which could increase the infection rate of HFMD. Moreover, infected children are more likely to be taken to see doctors on weekends rather than on weekdays.

There are substantial geographic variations in HFMD. Having an understanding of the geographical heterogeneities of meteorological effects between different areas may provide an appropriate way to estimate HFMD risks and powerful evidence for local authorities to establish precise warming system of HFMD occurrences. A study in East Asia observed variations in peak timing of HFMD occurrences among cities at different latitudes [35], which was caused mainly by meteorological factors. Hu et al [36] found that the morbidity of HFMD was significantly associated with population density and meteorological factors based on a Geographically Weighted Regression Model. However, to our knowledge, little literature has investigated potential between-cities discrepancies in meteorological effects on HFMD occurrences. In this study, both the exposure-response curves and the heterogeneity analysis indicated significant geographical differences in meteorological effects, which were significantly associated with the latitude and longitude with an explained heterogeneity of 31.38 %. We found that HFMD occurrences were generally more sensitive to temperature effects at lower latitude, but more sensitive to relative humidity effects at higher latitude. The effect discrepancies are mainly because the cities at lower attitudes generally have tropical climate and are also located nearby the Pacific Ocean, which may provide higher humid climate so that people living in these areas are less sensitive to the effects of increased relative humidity. Meanwhile, sea water helps to cool down for people living around the Pacific Ocean at lower latitude and increase their sensitivity to temperature to some extent. Additionally, HFMD occurrences were more sensitive to temperature and relative humidity effects at higher longitude, while the effects variations at lower longitude like Xuwen were narrow. The effect of relative humidity on HFMD in Guangning seems to be relatively higher than the effects in other cities when the humidity was below 50 % (Fig. 2b). The common parameter specifications were used to make the estimates comparable among cities but may not be optimal for Guangning. The daily average cases were 4.2 in Guangning, much less than that (30.7) in other seven cities. The small number of daily cases may lead to biased effect estimates in Guangning.

There are several limitations in the present study. Firstly, although the 5-year data in eight cities in Guangdong were used as possible as we can, the number of daily HFMD cases for subgroups in each city was too small to support subgroup analyses for meteorological effects. Secondly, we considered many geographic, demographic, and socio-economic factors at the city level but the situations are relatively homogeneous in the same province and there may be other HFMD-related information, such as the level of health care utilization, which would help to explore the potential explanation for heterogeneity. Further analyses for a diverse range of cities in whole China with varying socio-economic conditions would provide an overall understanding of meteorological effects on HFMD occurrences in China, which may be very valuable for the prevention and control of HFMD.

Conclusions

HFMD remains a crucial public health problem in Guangdong, China. HFMD cases clustered in male and children under 3 years old. Meteorological factors including temperature and relative humidity had delayed and non-linear effects on HFMD occurrences. The meteorological effects were closely linked to geographic longitude and latitude. These findings provide important evidence for further understanding of geographic variations of HFMD. The results can also help local authorities take corresponding interventions and measures to control HFMD development before reaching the peak risks. Further combination of HFMD data from other regions with more city-level information is warranted to construct the HFMD warning system.

Abbreviations

ARIMA:

Autoregressive integrated moving average model

BLUP:

Best linear unbiased prediction

China CDC:

Chinese Centre for Disease Control and Prevention

CI:

Confidence interval

CV-A16 :

Coxsackievirus A16

DLNM:

Distributed lag non-linear model

EV71 :

Enterovirus71

GAM:

Generalized additive model

HFMD:

Hand, foot, and mouth disease

MGAM:

Mixed generalized additive model

Q-AIC:

Likelihood Akaike information criteria for quasi-Poison

REML:

Restricted maximum likelihood

RR:

Relative risk

SARIMA:

Seasonal autoregressive integrated moving average

References

  1. Zou XN, Zhang XZ, Wang B, Qiu YT. Etiologic and epidemiologic analysis of hand, foot, and mouth disease in Guangzhou city: a review of 4,753 cases. Braz J Infect Dis. 2012;16:457–65.

    Article  PubMed  Google Scholar 

  2. World Health Organization. Hand, foot, and mouth disease surveillance summary. Available at: http://www.wpro.who.int/emerging_diseases/HFMD/en/. Accessed 1st Aug 2015.

  3. Xing W, Liao Q, Viboud C, Zhang J, Sun J, Wu JT, et al. Hand, foot, and mouth disease in China, 2008–12: an epidemiological study. Lancet Infect Dis. 2014;14:308–18.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Ang LW, Koh BK, Chan KP, Chua LT, James L, Goh KT. Epidemiology and control of hand, foot and mouth disease in Singapore, 2001-2007. Ann Acad Med Singapore. 2009;38:106–12.

    PubMed  Google Scholar 

  5. Podin Y, Gias ELM, Ong F, Leong Y-W, Yee S-F, Yusof MA, et al. Sentinel surveillance for human enterovirus 71 in Sarawak, Malaysia: lessons from the first 7 years. BMC Public Health. 2006;6:1471–2458.

    Article  Google Scholar 

  6. Chang HL, Chio CP, Su HJ, Liao CM, Lin CY, Shau WY, et al. The association between enterovirus 71 infections and meteorological parameters in Taiwan. PLoS One. 2012;7:e46845.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Lee C-CD, Tang J-H, Hwang J-S, Shigematsu M, Chan T-C. Effect of meteorological and geographical factors on the epidemics of hand, foot, and mouth diseasein island-type territory, East Asia. J Biomed Biotechnol. 2014;23.

  8. Wu H, Wang H, Wang Q, Xin Q, Lin H. The effect of meteorological factors on adolescent hand, foot, and mouth disease and associated effect modifiers. Glob Health Action. 2014;7:24664.

    Article  PubMed  Google Scholar 

  9. Yu L, Zhou L, Tan L, Jiang H, Wang Y, Wei S, et al. Application of a new hybrid model with seasonal auto-regressive integrated moving average (ARIMA) and nonlinear auto-regressive neural network (NARNN) in forecasting incidence cases of HFMD in Shenzhen, China. PLoS One. 2014;9:e98241.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Wei J, Hansen A, Liu Q, et al. The effect of meteorological variables on the transmission of hand, foot and mouth disease in four major cities of shanxi province, china: a time series data analysis (2009-2013). PLoS Negl Trop Dis. 2015;9:1–19.

    Article  Google Scholar 

  11. Chen C, Lin H, Li X, Lang L, Xiao X, Ding P, et al. Short-term effects of meteorological factors on children hand, foot and mouth disease in Guangzhou, China. Int J Biometeorol. 2014;58:1605–14.

    Article  PubMed  Google Scholar 

  12. Xie Y-h, Chongsuvivatwong V, Tang Z, McNeil EB, Tan Y. Spatio-temporal clustering of hand, foot, and mouth disease at the county level in Guangxi, China. PLoS One. 2014;9:e88065.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Ma W, Sun X, Song Y, Tao F, Feng W, He Y, et al. Applied mixed generalized additive model to assess the effect of temperature on the incidence of bacillary dysentery and its forecast. PLoS One. 2013;8:e62122.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Gasparrini A, Armstrong B, Kenward MG. Distributed lag non-linear models. Stat Med. 2010;29:2224–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Guo C, Yang L, Ou CQ, Li L, Zhuang Y, Yang J, et al. Malaria incidence from 2005-2013 and its associations with meteorological factors in Guangdong, China. Malar J. 2015;14:116.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Hii YL, Rocklöv J, Ng N. Short term effects of weather on hand, foot and mouth disease. PLoS One. 2011;6:e16796.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Jackson D, Riley R, White IR. Multivariate meta-analysis: potential and promise. Stat Med. 2011;30:2481–98.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Gasparrini A, Armstrong B, Kenward MG. Multivariate meta-analysis for non-linear and other multi-parameter associations. Stat Med. 2012;31:3821–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Deng T, Huang Y, Yu S, Gu J, Huang C, Xiao G, et al. Spatial-temporal clusters and risk factors of hand, foot, and mouth disease at the district level in Guangdong Province, China. PLoS One. 2013;8:e56943.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. World Health Organization. Western Pacific Regional Office of the World Health. Available at: http://www.wpro.who.int/emerging_diseases/HFMD/en/index.html. Accessed 1st Aug 2015.

  21. Kim SJ, Kim JH, Kang JH, Kim DS, Kim KH, Kim KH, et al. Risk factors for neurologic complications of hand, foot and mouth disease in the Republic of Korea, 2009. J Korean Med Sci. 2013;28:120–7.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Onozuka D, Hashizume M. The influence of temperature and humidity on the incidence of hand, foot, and mouth disease in Japan. Sci Total Environ. 2011;410–411:119–25.

    Article  PubMed  Google Scholar 

  23. Chen K-T, Chang H-L, Wang S-T, Cheng Y-T, Yang J-Y. Epidemiologic features of hand-foot-mouth disease and herpangina caused by enterovirus 71 in Taiwan, 1998–2005. Pediatrics. 2007;120:e244–52.

    Article  PubMed  Google Scholar 

  24. Chan KP, Goh KT, Chong CY, Teo ES, Lau G, Ling AE. Epidemic hand, foot and mouth disease caused by human enterovirus 71, Singapore. Emerg Infect Dis. 2003;9:78–85.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Ishimaru Y, Nakano S, Yamaoka K, Takami S. Outbreaks of hand, foot, and mouth disease by enterovirus 71. High incidence of complication disorders of central nervous system. Arch Dis Child. 1980;55:583–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Zhu Z, Zhu S, Guo X, Wang J, Wang D, Yan D, et al. Retrospective seroepidemiology indicated that human enterovirus 71 and coxsackievirus A16 circulated wildly in central and southern China before large-scale outbreaks from 2008. Virol J. 2010;7:300–5.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Chang LY, King CC, Hsu KH, Ning HC, Tsao KC, Li CC, et al. Risk factors of enterovirus 71 infection and associated hand, foot, and mouth disease/herpangina in children during an epidemic in Taiwan. Pediatrics. 2002;109:e88.

    Article  PubMed  Google Scholar 

  28. Li T, Yang Z, Di B, Wang M. Hand-foot-and-mouth disease and weather factors in Guangzhou, southern China. Epidemiol Infect. 2014;142:1741–50.

    Article  CAS  PubMed  Google Scholar 

  29. Huang Y, Deng T, Yu S, Gu J, Huang C, Xiao G, et al. Effect of meteorological variables on the incidence of hand, foot, and mouth disease in children: a time-series analysis in Guangzhou, China. BMC Infect Dis. 2013;13:134–42.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Lu P, Zhou B. Hand-foot-mouth disease 17. Radiology of Infectious Diseases. 2015;1:127.

    Article  Google Scholar 

  31. Arita M, Shimizu H, Nagata N, et al. Temperature-sensitive mutants of enterovirus 71 show attenuation in cynomolgus monkeys. J Gen Virol. 2005;86:1391–401.

    Article  CAS  PubMed  Google Scholar 

  32. Stanton GJ, Langford M, Baron S. Effect of interferon, elevated temperature, and cell type on replication of acute hemorrhagic conjunctivitis viruses. Infect Immun. 1977;18:370–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. Yu C-K, Chen C-C, Chen C-L, et al. Neutralizing antibody provided protection against enterovirus type 71 lethal challenge in neonatal mice. J Biomed Sci. 2000;7:523–8.

    Article  CAS  PubMed  Google Scholar 

  34. Sheikh S. Cloning and nucleotide sequence analysis of the P3 region of enterovirus 70. University of Ottawa (Canada); 1990

  35. Lee C-CD, Tang J-H, Hwang J-S, Shigematsu M, Chan T-C. Effect of meteorological and geographical factors on the epidemics of hand, foot, and mouth disease in island-type territory, East Asia. Biomed Res Int. 2015;2015:1–23.

    Google Scholar 

  36. Hu M, Li Z, Wang J, Jia L, Liao Y, Lai S, et al. Determinants of the incidence of hand, foot and mouth disease in China using geographically weighted regression models. PLoS One. 2012;7:e38978.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

None.

Funding

This work was supported by National Basic Research Program of China (973 Program) (Grant No. 2012CB955504), National Natural Science Foundation of China [81573249], and Science and Technology Planning Project of Guangdong Province [2013B021800041]. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and Additional file 5. Individual raw data are available from the authors upon reasonable request and with permission of Chinese Center for Disease Control and Prevention.

Authors’ contributions

QY and CQ initiated the study. QY and JY contributed to data collection. CG, YG, QQ and SQ cleaned and analyzed the data. CG, CQ, and QQ drafted the manuscript. All authors revised the manuscript critically and approved the final manuscript.

Authors’ information

None to provide.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

The study proposal was approved by Ethics Committee of Chinese Center for Disease Control and Prevention and Ethics Committee of Southern Medical University. Informed consent was not required because the data used in this study were routine surveillance information derived from Chinese Infection Diseases Monitor Information System and all data were recorded and analyzed anonymously.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Chun-Quan Ou or Qi-Yong Liu.

Additional files

Additional file 1:

Spearman rank analysis between HFMD cases and meteorological factors in eight cities in Guangdong, 2009-2013. ** denotes significant correlations at the 0.01 level, and * represent significant correlations at the 0.05 level. (DOCX 20 kb)

Additional file 2:

The cumulative effects of sunshine and wind speed on HFMD over lag 0–14 days in Guangdong, 2009–2013. The bold red line represents the pooled effects, and the dashed lines represent the city-specific estimates. Reference values were the medians, that is, 4.78 h for sunshine and 2.70 m/s for wind speed. (DOC 52 kb)

Additional file 3:

Descriptive analysis for study-level variables in eight cities in Guangdong, 2009–2013. (DOCX 18 kb)

Additional file 4:

The stratified analyses of meteorological effects by seasons. The bold red line represents the pooled effects, and the dashed lines represent the city-specific estimates. Reference values were the medians, that is, 4.78 h for sunshine and 2.70 m/s for wind speed. (DOCX 1400 kb)

Additional file 5:

Daily aggregated data of HFMD cases and the demographic, socio-economic, and meteorological factors in each of eight cities in Guangdong, 2009–2013. The variables represent the date under study, city code, the number of daily HFMD cases, daily mean temperature, mean wind, mean air atmospheric pressure, mean humidity, rain, sunshine hour, and city-level information including longitude, latitude, GDP per capital, population, population density, the ratio of male to female population, day of the week, holiday, the ratio of children under 5 years old, living space per capita, the ratio of population flow. (CSV 1389 kb)

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Guo, C., Yang, J., Guo, Y. et al. Short-term effects of meteorological factors on pediatric hand, foot, and mouth disease in Guangdong, China: a multi-city time-series analysis. BMC Infect Dis 16, 524 (2016). https://0-doi-org.brum.beds.ac.uk/10.1186/s12879-016-1846-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12879-016-1846-y

Keywords