Hostname: page-component-788cddb947-kc5xb Total loading time: 0 Render date: 2024-10-13T17:07:48.117Z Has data issue: false hasContentIssue false

Seasonality of the transmissibility of hand, foot and mouth disease: a modelling study in Xiamen City, China

Published online by Cambridge University Press:  30 December 2019

Zehong Huang
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Xiamen City, Fujian Province, People's Republic of China
Mingzhai Wang
Affiliation:
Xiamen Centre for Disease Control and Prevention, Xiamen City, Fujian Province, People's Republic of China
Luxia Qiu
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Xiamen City, Fujian Province, People's Republic of China
Ning Wang
Affiliation:
Respiratory Department, Shanghai General Hospital, Shanghai, People's Republic of China
Zeyu Zhao
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Xiamen City, Fujian Province, People's Republic of China
Jia Rui
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Xiamen City, Fujian Province, People's Republic of China
Yao Wang
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Xiamen City, Fujian Province, People's Republic of China
Xingchun Liu
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Xiamen City, Fujian Province, People's Republic of China
Mikah Ngwanguong Hannah
Affiliation:
Medical College, Xiamen University, Xiamen City, Fujian Province, People's Republic of China
Benhua Zhao
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Xiamen City, Fujian Province, People's Republic of China
Yanhua Su*
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Xiamen City, Fujian Province, People's Republic of China
Bin Zhao*
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, Laboratory Department, Xiang'an Hospital of Xiamen University, Xiamen City, Fujian Province, People's Republic of China
Tianmu Chen*
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Xiamen City, Fujian Province, People's Republic of China
*
Authors for correspondence: Tianmu Chen, E-mail: 13698665@qq.com; Yanhua Su, E-mail: suyanhua813@xmu.edu.cn; Bin Zhao, E-mail: 393603468@qq.com
Authors for correspondence: Tianmu Chen, E-mail: 13698665@qq.com; Yanhua Su, E-mail: suyanhua813@xmu.edu.cn; Bin Zhao, E-mail: 393603468@qq.com
Authors for correspondence: Tianmu Chen, E-mail: 13698665@qq.com; Yanhua Su, E-mail: suyanhua813@xmu.edu.cn; Bin Zhao, E-mail: 393603468@qq.com
Rights & Permissions [Opens in a new window]

Abstract

This study attempts to figure out the seasonality of the transmissibility of hand, foot and mouth disease (HFMD). A mathematical model was established to calculate the transmissibility based on the reported data for HFMD in Xiamen City, China from 2014 to 2018. The transmissibility was measured by effective reproduction number (Reff) in order to evaluate the seasonal characteristics of HFMD. A total of 43 659 HFMD cases were reported in Xiamen, for the period 2014 to 2018. The median of annual incidence was 221.87 per 100 000 persons (range: 167.98/100,000–283.34/100 000). The reported data had a great fitting effect with the model (R2 = 0.9212, P < 0.0001), it has been shown that there are two epidemic peaks of HFMD in Xiamen every year. Both incidence and effective reproduction number had seasonal characteristics. The peak of incidence, 1–2 months later than the effective reproduction number, occurred in Summer and Autumn, that is, June and October each year. Both the incidence and transmissibility of HFMD have obvious seasonal characteristics, and two annual epidemic peaks as well. The peak of incidence is 1–2 months later than Reff.

Type
Original Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s) 2019

Introduction

Hand, foot and mouth disease (HFMD), resulting from a group of enterovirus (mainly include enterovirus 71 (EV-A71) and coxsackievirus A16 (CV-A16) [Reference Klein and Chong1]), is an infectious disease common in children under 5 years of age. Clinical manifestations of HFMD are characterised by self-limited fever, typical skin rash, with or without mouth ulcers, complications of the nervous or cardiopulmonary systems which can also be found in several cases, and even results in death sometimes [Reference Koh2].

HFMD has spread widely in the Western Pacific Region [3]. Millions of infections were reported to have caused 96 900 Disability Adjusted of Life Years in some Asian countries [Reference Koh4], and had led to an extremely heavy disease burden in Asia and to the world at large.

The significance of prevention and control of HFMD is particularly important in China. More than 2.4 million infections are reported in China per year [3], with the enormous population base resulting in severe complications and even deaths. The outbreak of HFMD occurs almost every year [Reference Shu-Chun, En-Tzu and Ding-Ping5]. All of the above make HFMD into a non-negligible emergency.

Since the first report in Shanghai, 1981, varying degrees of outbreaks have taken place in Dandong, Wuhan, Hong Kong, Taiwan etc. [Reference Shu-Chun, En-Tzu and Ding-Ping5], which has brought bitter suffering to the state. Until May 2008, HFMD was defined as a C-class disease in Mainland China, which opened its surveillance and brought a new stage for the study of anti-HFMD measures.

Finding the potential factors affecting transmissibility of HFMD and establishing a suitable HFMD transmission model to predict the epidemic trend will undoubtedly bring important clues to the prevention and control of HFMD [Reference Zhu6]. In recent years, several mathematical models have been developed to calculate the transmissibility of HFMD including spatial factors [Reference Ventarola, Bordone and Silverberg7], climatic factors [Reference Bo8] and even environmental and socioeconomic factors [Reference Song9], which have a great effect in predicting its epidemic trend. But there are still some deficiencies among the studies of seasonal transmissibility of HFMD. The HFMD has obvious seasonality, which is in the cross of spring and summer, summer and autumn and the seasonal fluctuation increases as the location approaches the equator [Reference Li10]. The seasonal study of HFMD can provide an effective basis for preventing and controlling HFMD.

Therefore, this study was based on the susceptible-exposed-infectious-asymptomatic-removed (SEIAR) model, combined with the epidemiological characteristics of HFMD transmission, and established a SEIAR dynamic model of HFMD epidemic with seasonality. The daily incidence data from 2014 to 2018 in Xiamen were collected to evaluate the seasonal transmissibility of HFMD, so as to provide clues for the development of disease prevention and control measures.

Materials and methods

Study design

In this study, the mathematical epidemiological method was used to mine the data by the SEIAR model in order to master the characteristics and transmissibility of HFMD.

Date collection

Xiamen, a large city in Fujian Province, is located on the south-east coast of China and belongs to the subtropical monsoon climate area. According to the 2018 census, the resident population of Xiamen has exceeded 4 million. The data of reported HFMD cases from Xiamen were extracted from the China Information System for Disease Control and Prevention from January 2014 to December 2018. All cases included were diagnosed in accordance with the Hand, Foot and Mouth Disease Diagnosis and Treatment Guide (2010 Edition) and Diagnosis for Hand, Foot and Mouth Disease (WS-588 2018), issued by the National Health Commission of the People's Republic of China.

Transmission model

Our previous study shows that the susceptible-infectious-removed model could be used to calculate the transmissibility of HFMD based on weekly incidence data [Reference Chen11]. In this study, daily incidence data were employed. Therefore, the incubation period and asymptomatic infection were considered. Thus, a SEIAR model [Reference Chen11, Reference Wang, Xiao and Cheke12] was established in this study. According to the natural history of HFMD [3, Reference Merzel Sabovic13], the population was divided into susceptible individuals (S), exposed but not yet infected individuals (E), infectious and symptomatic individuals (I), asymptomatic individuals (A) which means infectious but not symptomatic and recovery individuals (R). The model began with the state of susceptible individuals (S), and then entered the incubation state (E) when infected. The subsequent transitions included state of infection with symptomatic (I) and asymptomatic state (A), both of them would end up with recovery state (R). Besides, all the states above had specific death rates, and newborns joined the susceptible individuals (S) with a specific birth ratio. The model diagram is shown in Fig. 1.

Fig. 1. Flow chart of the SEIAR model for HFMD.

The SEIAR model used in this study was based on the following assumptions:

  1. 1) The model assumed that HFMD cannot propagate vertically, so the new individuals born in all kinds of people are susceptible. Set the birth rate to br, and the natural mortality rate to dr. The mortality rate of infectious individual f was obtained by references, only 0.03%. Therefore, it is known that the mortality rate of all kinds of people in the disease spectrum is low, and the mortality rate of population attributable to HFMD is even lower. So the mortality rate of the whole population was approximately replaced by the mortality rate of several types of people except infectious individual, whose mortality rate was set as the sum of the mortality of the whole population and the mortality of HFMD.

  2. 2) Transmission of HFMD occurs via person–person, and the transmissibility between infectious individual and asymptomatic one may be different. So, the k was defined as the relative transmissibility rate of asymptomatic to symptomatic individuals. At the same time, we assumed the susceptible individuals will be potentially infectious (will be infected) as long as they are in contact with infectious individuals or asymptomatic individuals, and the coefficient of the infection rate was set as β.

  3. 3) We considered that there was a certain proportion of exposed individuals pE (0⩽p⩽1) transformed into infectious individuals I after incubation, another part of exposed individuals (1 − p) E were transformed into asymptomatic individuals after incubation as well. At time t, the development speed from the E to I pathway is same as the E to A pathway and both of them were proportional to the number of exposed individuals. The proportional coefficient of the former was , and the latter was (1 − p) ω.

  4. 4) Infectious and asymptomatic individuals may move to R, with speed of recovering being in direct proportion to the number of individuals. The proportional coefficients were γ and γ′ respectively.

  5. 5) Infectious and asymptomatic individuals can produce antibodies after recovery. This study assumed that the recovered individuals have permanent immunity which means they would no longer be infected, so R was set as the end of the model.

  6. 6) Perinatal transmission of enteroviruses is clearly documented [Reference Tassin14, Reference Bendig15]. The main modes of perinatal transmission are by intrapartum exposure to maternal blood and/or genital secretions, as well as the faecal–oral and respiratory routes after delivery [Reference Harik and DeBiasi16]. However, evidence of vertical transmission of enteroviruses associated with HFMD has not been confirmed by specific studies. Serological data on enterovirus infections during pregnancy were also lacking, because it's hard to diagnose enterovirus infection in adult [Reference Méreaux17]. So we assumed that HFMD cannot propagate vertically.

The differential equations of the model were used to describe the dynamic changes of each state. The corresponding model equations are as follows:

$$\displaystyle{{{\rm d}S} \over {{\rm d}t}} = brN-\beta S(I + kA)-drS$$
$$\displaystyle{{{\rm d}E} \over {{\rm d}t}} = \beta S(I + kA)-\omega E-drE$$
$$\displaystyle{{{\rm d}I} \over {{\rm d}t}} = p\omega E-\gamma I-drI-fI$$
$$\displaystyle{{{\rm d}A} \over {{\rm d}t}} = (1-p)\omega E-{\gamma} ^{\prime}A-drA$$
$$\displaystyle{{{\rm d}R} \over {{\rm d}t}} = \gamma I + \gamma ^{\prime}A-drR$$

The left side of the equations represents the instantaneous velocities of S, E, I, A and R at time t.

The seasonality of the transmission

In this study, the seasonality of the transmission was considered. According to the mechanism of the SEIAR model we built, the seasonality should be dynamic, focusing on β. Therefore, a trigonometric function was adopted and shown as follows:

$$\beta = \beta _0\left[ {1 + \sin \left( {\displaystyle{{2\pi (t + \alpha )} \over T}} \right)} \right]$$

In the equation, β 0, t, α and T refer to the baseline of the transmission relative rate, time, a constant which adjusts the position of time, and the time span of the season cycle respectively.

Estimation of parameters

Of p, ω, γ and f, the four parameters in the model (Table 1) came from the literatures. The range of the symptomatic infection ratio was from 19% to 47% [Reference Koh2, Reference Xing18, Reference Chang19], which was chosen as the median 44.23%, that is, P = 0.4423. The range of incubation 1/ω was from 3 to 7 days [Reference Koh2, Reference Chang20], with 5 days chosen as median (with a median of 5 days), that is, ω = 0.2. The course of symptomatic infection was 2 weeks [Reference Du21, Reference Takahashi22], that is, recovery rate of the infection γ = 0.0714. The range of course of asymptomatic infection 1/γ′ was 2 to 4 weeks [3, Reference Xing18, Reference Chang19], with 3 weeks chosen as median (with a median of 3 weeks), that is, the recovery rate of the asymptomatic γ′ = 0.0476. The range of the fatality rate was from 0.0001 to 0.0005 [Reference Yang23-Reference Zhang25], with 0.0003 chosen as median (with a median of 5 days). The parameters β were estimated by curve fitting. There was no clearly relevant data or reference to support the parameter k which still remained uncertain. Thus, k = 1 was selected for calculation, and sensitivity analysis was carried out to figure its influence on the model. The birth rate (br) and the mortality rate (dr) were collected from Xiamen Centre for Disease Control and Prevention.

Table 1. Parameter definitions and values

CDC: Center for Disease Control and Prevention.

Transmissibility of HFMD

In this study, the effective reproduction number R eff was used to measure the transmissibility of HFMD, thought of as the number of new cases which was produced by a typical case during the period of infection. If it is more than 1, it indicates that the transmissibility of the disease is strong and the disease will be able to spread in the population and the opposite is true when the effective reproduction number is less than 1.

$$R_{{\rm eff}} = \beta S\left( {\displaystyle{\,p \over \gamma} + \displaystyle{{({\rm 1}-p)k} \over {\gamma^{\prime}}}} \right)$$

Sensitivity analysis

There was no clearly relevant data or literature to support the parameter k which is still uncertain. Therefore, the sensitivity was tested by dividing k into six values within the range of 0 to 1, that is, k = 0, k = 0. 2, k = 0. 4, k = 0. 6, k = 0. 8 and k = 1 (indicates transmissibility rate of asymptomatic individuals is the same as symptomatic ones).

Simulation methods and statistical analysis

In this study, Berkeley Madonna 8.3.18 (developed by Robert Macey and George Oster of the University of California at Berkeley. Copyright ©1993–2001 Robert I. Macey & George F. Oster) was employed for fitting the model with the data to calculate the parameters, while the transmissibility was calculated based on the parameters. The determination coefficient (R 2) was used to evaluate the goodness of fit of the model simulation. The regression analysis, which was used to calculate R 2, was performed by the software SPSS 21.0.0.0 (IBM Corp., Armonk, NY, USA).

Results

Incidence of HFMD in Xiamen

A total of 43 659 cases of HFMD were reported in Xiamen for the period 2014 to 2018 (Fig. 2). The median annual reported incidence was 221.87 per 100 000 persons (range: 167.98 per 100 000 persons–283.34 per 100 000 persons). The number of reported cases and the incidence showed an upward trend over the past 5 years except in 2015, which was shown to have a slight decrease. The upward trend was the most obvious between 2015 and 2017, and slowed down after 2017.

Fig. 2. Number of reported cases and incidence of HFMD in Xiamen City from 2014 to 2018.

Simulation results of HFMD transmission by the SEIAR model

The results of curve fitting between the model and the reported cases are illustrated in Fig. 3. It was found that there are two epidemic cycles in summer and autumn each year, whose peaks were in June and October, through the analysis of the reported cases and curve of fitting.

Fig. 3. Curve fitting results based on the data in Xiamen, 2014–2018.

As was depicted in the result of curve fitting, it was shown that the number of cases had a periodic increase trend except in 2015, which had a slight decrease. After 2015, the increasing speed of reported cases was apparent in 2017 and 2018, whose number of reported cases increased obviously than the previous 3 years.

Table 2 lists the fitting values of β 0, R 2 and P based on the data from 2014 to 2015. The mean of β 0 is 2.9270 × 10−8, which means that the model has a great fitting effect with the data of HFMD (R 2 = 0.9212, P < 0.0001).

Table 2. Estimated values of β 0 and goodness of fit of the model to reported HFMD data from 2014 to 2018 in Xiamen

The transmissibility of HFMD

The median of R eff of HFMD in Xiamen from 2014 to 2018 was 1.1048 (range: 1.4045 × 10−5–3.7459). The difference was significant in each year, showing obvious periodic characteristic. There were two transmission peaks each year (from April to May and August to September), and the transmissibility decreased rapidly then. During the peak period, the transmissibility of HFMD can reach 200 000 times more than the lowest. Table 3 shows the details.

Table 3. Effective reproduction number (R eff)

As illustrated in Fig. 4, R eff was consistent with the incidence data, whose peaks were mostly 1–2 months before the peaks of the incidence. Both R eff and incidence revealed the obvious periodic characteristic, with two transmission peaks each year.

Fig. 4. The daily incidence rate and R eff in Xiamen, 2014–2018.

Sensitivity analysis

The fitting curve with six values of parameter k in the range of 0–1 had high coincidence degree, which indicates that the SEIAR model used in this study was not sensitive to it. The specific results are elucidated in Fig. 5.

Fig. 5. Sensitivity analysis of k (range: 0–1).

Discussion

In this study, we applied the SEIAR model with seasonal correction for the first time to clarify the transmissibility of HFMD, which would offer a guideline for the study of transmissibility, prevention and control measures of HFMD.

Validity of model

The following three points guaranteed the validity of the model: the fitting effect of the model was tested by regression analysis (R 2 = 0.9212). At the same time, sensitivity was examined by parameter k, which indicates that the SEIAR model used in this study is not sensitive to it. Besides, the model was based on daily data, which increase the sample and accuracy of the study.

Not only was the above enough to guarantee the validity of the model, but also ensure the reliability of the research, which paves the way for further research.

The epidemiological characteristics of HFMD

Our finding showed that the number of HFMD cases increased in Xiamen from 2014 to 2018, except 2015. The decrease in 2015 might be due to the launch of the vaccine or the decreasing of annual rainfall (1409 mm, 289 mm less than the usual, data from Xiamen Meteorological Bureau). The incidence increased remarkably from 2015 to 2017 and slightly from 2017 to 2018. The reason of the incidence changing remains unclear and more research studies are needed to explore the mechanism of the transmission.

Obvious seasonality and the two peaks twice a year of HFMD had been shown in the study, which were consistent with the previous studies [Reference Zou26Reference Wang29], and the conclusion on time distribution of cases in summer and autumn also fits basically the previous studies [Reference Zhang28, Reference Liu30]. However, it can be found that there was no significant difference between the two incidence peaks, which was different from the point of view of Zhang et al. [Reference Zhang28].

Assessment of transmissibility of HFMD

The transmissibility of HFMD varies greatly within a year, the highest peak value was 200 000 times more than the lowest, and the annual median of R eff was 2.2–3.7, which indicated that the transmission rate of HFMD in Xiamen is high that will result in the increase in the number of patients with HFMD. The range of R eff in previous studies was obviously disparate because of different experimental designs. Such as 1.14–1.24 by Yang et al. [Reference Yang31], which was slightly lower than this study, it might be influenced by its focus that placed more emphasis on community. The range of R eff was between 2.5 to 5.5 from two studies in Hong Kong and Singapore [Reference Lim32, Reference Ma33], which was approximately consistent with this study.

The seasonal characteristics of the transmissibility of HFMD are notable. Its first and second peaks are from March to April and August to September every year. The transmissibility during the other months is relatively low. The peak of transmission in autumn was stronger than in spring. Some studies showed that CA-V16 might lead to the peak of HFMD in autumn [Reference Zhuang34]. The intensity difference between the peaks of HFMD in Xiamen may be related to the etiological characteristics of HFMD.

The periodicity of R eff is earlier than the incidence of HFMD, and its peak is 1–2 months earlier than the latter, which indicates that the transmissibility of HFMD is beforehand the outbreak of HFMD and offers some indication for the prevention and control of HFMD.

The significance for the disease control

At present, the main measures to prevent and control HFMD are by enterovirus vaccine, as well as the early identification and intervention of disease outbreaks [3]. According to the conclusion of this study, the identification time of epidemic situation should be adjusted to the time before the peak of incidence, that is, the peak of R eff, in order to accord with the epidemiological characteristics of the disease, and to prevent and control it earlier, whose specific time should be 1–2 months ahead of the peak of incidence.

Greater effect will be received if the correct time of intervention is chosen ahead of the peak in transmissibility of HFMD. The measures should be applied to the susceptible population, that is, from February to March, July to August of each year, so as to control the spread of epidemic situation to the maximum extent and get the greatest prevention and control effect. For example, the school term for the children under the age of 10 years should be changed to later dates in the areas with a high incidence of HFMD, to cause a change in the peak and interfere with the spread of infection.

The results of this study also have guiding significance for vaccination. The previous study showed that the durability of neutralising antibody induced by EV-A71 vaccine was limited, and it was necessary to strengthen the injection for long-term protection [Reference Aswathyraj35]. The seasonal transmissibility of HFMD should be considered into vaccine injection time and strengthening injection time, which can be placed in February and July of each year, to increase the income of vaccine and achieve more effective prevention and control.

Limitations

Limited by the availability of data, there are some limitations in this study. Some factors (gender and age) which might affect the transmissibility of the disease were not considered in the model. In addition, there may be some difference in the fitting of the epidemic situation and the estimation of the transmissibility in other high-latitude areas, as this study is confined to Xiamen. Besides, for the reinfection of HFMD, there could be a 2.02% probability according to the previous literature [Reference Chen36]. Considering that the scale of the proportion, the maintenance time of the antibody is greatly affected by various factors as well. Therefore, this study abandoned pathway of recovered/removed individuals transforming into susceptible individuals when establishing the model.

The vertical transmission of HFMD was not considered, which might slightly impact the result of our study.

Conclusions

Based on the incidence data of Xiamen from 2014 to 2018, the SEIAR model has a great fitting effect in the incidence situation of HFMD. It's found that the incidence and transmissibility of HFMD are seasonal, which have two seasonal peaks. The peak of incidence is in summer and autumn every year, and the peak of transmissibility is 1–2 months ahead of the peak of incidence.

Data

Extra data is available by emailing to Dr Mingzhai Wang () on reasonable request.

Acknowledgements

The authors thank all study participants for providing the data and field investigators for collecting the data.

Author contributions

TC, ZH and LQ designed the study. MW collected data. TC, ZH, LQ, NW, ZZ, JR, YW, XL, BZ and YS performed the analysis. TC, ZH, LQ, NW and MNH wrote the first draft of this paper. All authors contributed to the writing of the manuscript.

Financial support

This study was partly supported by the Open Research Fund of State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics (SKLVD2018KF001 and SKLVD2018KF002).

Conflict of interest

None declared.

Ethical standards

This effort of disease control was part of CDC's routine responsibility in Xiamen City, China. Therefore, institutional review and informed consent were not required for this study. All data analysed were anonymised.

Footnotes

*

These authors contributed equally to this study

References

1.Klein, M and Chong, P (2015) Is a multivalent hand, foot, and mouth disease vaccine feasible? Human Vaccines & Immunotherapeutics 11, 26882704.CrossRefGoogle ScholarPubMed
2.Koh, WM et al. (2016) The epidemiology of hand, foot and mouth disease in Asia: a systematic review and analysis. The Pediatric Infectious Disease Journal 35, e285e300.CrossRefGoogle ScholarPubMed
3.World Health Organization (2011) A Guide to Clinical Management and Public Health Response for Hand, Foot and Mouth Disease (HFMD). Manila: WHO Regional Office for the Western Pacific.Google Scholar
4.Koh, WM et al. (2018) Severity and burden of hand, foot and mouth disease in Asia: a modelling study. BMJ Global Health 3, e000442e000442.CrossRefGoogle ScholarPubMed
5.Shu-Chun, C, En-Tzu, W and Ding-Ping, L (2013) A review of prevention and control for enterovirus infections in Asia. Epidemiology Bulletin 29, 101117.Google Scholar
6.Zhu, Q et al. (2011) Surveillance of hand, foot, and mouth disease in mainland China (2008–2009). Biomedical and Environmental Sciences 24, 349356.Google Scholar
7.Ventarola, D, Bordone, L and Silverberg, N (2015) Update on hand-foot-and-mouth disease. Clinics in Dermatology 33, 340346.CrossRefGoogle ScholarPubMed
8.Bo, Y-C et al. (2014) Using an autologistic regression model to identify spatial risk factors and spatial risk patterns of hand, foot and mouth disease (HFMD) in Mainland China. BMC Public Health 14, 358.CrossRefGoogle Scholar
9.Song, C et al. (2019) Exploring spatiotemporal nonstationary effects of climate factors on hand, foot, and mouth disease using Bayesian spatiotemporally varying coefficients (STVC) model in Sichuan, China. Science of The Total Environment 648, 550560.CrossRefGoogle ScholarPubMed
10.Li, L et al. (2018) A spatiotemporal mixed model to assess the influence of environmental and socioeconomic factors on the incidence of hand, foot and mouth disease. BMC Public Health 18, 274.CrossRefGoogle ScholarPubMed
11.Chen, S et al. (2019) Estimating the transmissibility of hand, foot, and mouth disease by a dynamic model. Public Health 174, 4248.CrossRefGoogle ScholarPubMed
12.Wang, J, Xiao, Y and Cheke, RA (2016) Modelling the effects of contaminated environments on HFMD infections in mainland China. Biosystems 140, 17.CrossRefGoogle ScholarPubMed
13.Merzel Sabovic, EK et al. (2019) Atypical hand, foot, and mouth disease in an adult patient: a case report and literature review. Acta dermatovenerologica Alpina, Pannonica, et Adriatica 28, 8588.CrossRefGoogle Scholar
14.Tassin, M et al. (2014) A case of congenital echovirus 11 infection acquired early in pregnancy. Journal of Clinical Virology 59, 7173.CrossRefGoogle ScholarPubMed
15.Bendig, JW et al. (2003) Coxsackievirus B3 sequences in the blood of a neonate with congenital myocarditis, plus serological evidence of maternal infection. Journal of Medical Virology 70, 606609.CrossRefGoogle ScholarPubMed
16.Harik, N and DeBiasi, RL (2018) Neonatal nonpolio enterovirus and parechovirus infections. Seminars in Perinatology 42, 191197.CrossRefGoogle ScholarPubMed
17.Méreaux, J et al. (2017) Enterovirus infection during pregnancy: underestimated cause of fetal and neonatal complications? Gynecologie, Obstetrique, Fertilite & Senologie 45, 231237.CrossRefGoogle ScholarPubMed
18.Xing, W et al. (2014) Hand, foot, and mouth disease in China, 2008–12: an epidemiological study. The Lancet Infectious Diseases 14, 308318.CrossRefGoogle Scholar
19.Chang, L-Y et al. (2002) Risk factors of enterovirus 71 infection and associated hand, foot, and mouth disease/herpangina in children during an epidemic in Taiwan. Pediatrics 109, e88.CrossRefGoogle ScholarPubMed
20.Chang, L-Y et al. (2004) Transmission and clinical features of enterovirus 71 infections in household contacts in Taiwan. Jama 291, 222227.CrossRefGoogle ScholarPubMed
21.Du, Z et al. (2017) Estimating the basic reproduction rate of HFMD using the time series SIR model in Guangdong, China. PLoS ONE 12, e0179623.CrossRefGoogle ScholarPubMed
22.Takahashi, S et al. (2016) Hand, foot, and mouth disease in China: modeling epidemic dynamics of enterovirus serotypes and implications for vaccination. PLoS Medicine 13, e1001958.CrossRefGoogle ScholarPubMed
23.Yang, Z et al. (2017) Estimating the incubation period of hand, foot and mouth disease for children in different age groups. Scientific Reports 7, 16464.CrossRefGoogle ScholarPubMed
24.Liu, S-L et al. (2015) Comparative epidemiology and virology of fatal and nonfatal cases of hand, foot and mouth disease in mainland China from 2008 to 2014. Reviews in Medical Virology 25, 115128.CrossRefGoogle ScholarPubMed
25.Zhang, J et al. (2015) Statistical monitoring of the hand, foot and mouth disease in China. Biometrics 71, 841850.CrossRefGoogle Scholar
26.Zou, X-N et al. (2012) Etiologic and epidemiologic analysis of hand, foot, and mouth disease in Guangzhou city: a review of 4,753 cases. The Brazilian Journal of Infectious Diseases 16, 457465.CrossRefGoogle Scholar
27.Deng, T et al. (2013) Spatial-temporal clusters and risk factors of hand, foot, and mouth disease at the district level in Guangdong Province, China. PLoS ONE 8, e56943.CrossRefGoogle ScholarPubMed
28.Zhang, W et al. (2015) An epidemic analysis of hand, foot, and mouth disease in Zunyi, China between 2012 and 2014. Saudi Medical Journal 36, 593598.CrossRefGoogle ScholarPubMed
29.Wang, Y et al. (2011) Hand, foot, and mouth disease in China: patterns of spread and transmissibility. Epidemiology 22, 781792.CrossRefGoogle ScholarPubMed
30.Liu, W et al. (2014) Co-circulation and genomic recombination of coxsackievirus A16 and enterovirus 71 during a large outbreak of hand, foot, and mouth disease in Central China. PLoS ONE 9, e96051.CrossRefGoogle ScholarPubMed
31.Yang, B et al. (2016) Transmission of hand, foot and mouth disease and its potential driving factors in Hong Kong. Scientific Reports 6, 27500.CrossRefGoogle ScholarPubMed
32.Lim, CTK et al. (2016) Basic reproduction number of coxsackievirus type A6 and A16 and enterovirus 71: estimates from outbreaks of hand, foot and mouth disease in Singapore, a tropical city-state. Epidemiology and Infection 144, 10281034.CrossRefGoogle ScholarPubMed
33.Ma, E et al. (2011) Estimation of the basic reproduction number of enterovirus 71 and coxsackievirus A16 in hand, foot, and mouth disease outbreaks. The Pediatric Infectious Disease Journal 30, 675679.CrossRefGoogle ScholarPubMed
34.Zhuang, ZC et al. (2015) Epidemiological research on hand, foot, and mouth disease in mainland China. Viruses 7, 64006411.CrossRefGoogle ScholarPubMed
35.Aswathyraj, S et al. (2016) Hand, foot and mouth disease (HFMD): emerging epidemiology and the need for a vaccine strategy. Medical Microbiology and Immunology 205, 397407.CrossRefGoogle ScholarPubMed
36.Chen, GP et al. (2016) Epidemiological characteristics and influential factors of hand, foot and mouth disease (HFMD) reinfection in children in Anhui province. Epidemiology and Infection 144, 153160.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1. Flow chart of the SEIAR model for HFMD.

Figure 1

Table 1. Parameter definitions and values

Figure 2

Fig. 2. Number of reported cases and incidence of HFMD in Xiamen City from 2014 to 2018.

Figure 3

Fig. 3. Curve fitting results based on the data in Xiamen, 2014–2018.

Figure 4

Table 2. Estimated values of β0 and goodness of fit of the model to reported HFMD data from 2014 to 2018 in Xiamen

Figure 5

Table 3. Effective reproduction number (Reff)

Figure 6

Fig. 4. The daily incidence rate and Reff in Xiamen, 2014–2018.

Figure 7

Fig. 5. Sensitivity analysis of k (range: 0–1).