 Research
 Open Access
 Published:
The burden of hand, foot, and mouth disease among children under different vaccination scenarios in China: a dynamic modelling study
BMC Infectious Diseases volume 21, Article number: 650 (2021)
Abstract
Background
Hand, foot, and mouth disease (HFMD) is a common illness in young children. A monovalent vaccine has been developed in China protecting against enterovirus71, bivalent vaccines preventing HFMD caused by two viruses are under development.
Objective
To predict and compare the incidence of HFMD under different vaccination scenarios in China.
Methods
We developed a compartmental model to capture enterovirus transmission and the natural history of HFMD in children aged 0–5, and calibrated to reported cases in the same agegroup from 2015 to 2018. We compared the following vaccination scenarios: different combinations of monovalent and bivalent vaccine; a program of constant vaccination to that of pulse vaccination prior to seasonal outbreaks.
Results
We estimate 1,982,819, 2,258,846, 1,948,522 and 2,398,566 cases from 2015 to 2018. Increased coverage of monovalent vaccine from 0 to 80% is predicted to decrease the cases by 797,262 (49.1%). Use of bivalent vaccine at an 80% coverage level would decrease the cases by 828,560. Use of a 2.0× pulse vaccination for the bivalent vaccine in addition to 80% coverage would reduce cases by over one million. The estimated R_{0} for HFMD in 2015–2018 was 1.08, 1.10, 1.35 and 1.17.
Conclusions
Our results point to the benefit of bivalent vaccine and using a pulse vaccination in specific months over routine vaccination. Other ways to control HFMD include isolation of patients in the early stage of dissemination, more frequent handwashing and ventilation, and better treatment options for patients.
Introduction
Hand, foot, and mouth disease (HFMD) is a common infectious disease mainly caused by various enteroviruses. HFMD usually affects children under age of five, with a incidence rate of approximately 2400 cases per 100,000 in 2018 in China [1, 2], and occurs more often in children under three [3]. After 2007, the HFMD epidemic was in an uptrend in China accompanied by serious outbreaks [4]. Although HFMD is usually selflimiting, it can result in complications associated with the central nervous system or death once progressing to severe cases [5, 6]. There are more than 20 types of enterovirus leading to HFMD and Enterovirus 71 (EV71) and Coxsackie virus A16 (CVA16) are the most commonly reported [3]. EV71 accounts for 70% severe HFMD cases and 90% HFMDrelated deaths in mainland China [7].
Mathematical models of enterovirus transmission are increasingly being used to design and inform public health interventions to reduce the spread of HFMD. For example, Chuo et al. developed a deterministic susceptibleinfectiousrecovered (SIR) model to predict the number of infected people during an HFMD occurred in Sarawak, Malaysia and found that the disease spread quite rapidly and the parameter that may be able to control the number of susceptible persons [8]. Yang et al. used a SEIR model, adding a compartment to represent the exposed people (E) to capture the incubation period of HFMD and showed the spread trend of HFMD in China, and the basic reproduction number indicated that it was an outbreak [9]. By adding in a quarantined compartment (Q) in a SEIQR model, Liu simulated HFMD transmission and found that quarantine is the optimal strategy for HFMD protection [10].
EV71 vaccine is a new vaccine recently developed in China. The Vero cellbased EV71 inactivated vaccine has been shown to induce immune responses to EV71 among infants and young children between 6 and 59 months of age [11]. By March 2013, this inactivated vaccines had completed phase IIII clinical trials in Mainland China, showing good safety and greater than 95% protection against HFMD caused by EV71 [11,12,13]. Three manufactures in mainland China have recently completed phase III trials for EV71 vaccines and have been approved by China Food and Drug Administration in December 2015 [14]. In total, three monovalent EV71 vaccines have been licensed while bivalent EV71/CA16 vaccines are also under development and have been shown to induce protective immunity against both EV71 and CVA16 in vivo and in mice [15]. Further empirical studies are needed to assess the safety and effectiveness of these different vaccines in preventing HFMD.
As empiric data emerge on the safety and individuallevel effectiveness of these vaccines, what remains unknown is their potential populationlevel impact on outbreaks and incidence of HFMD. Existing and future EV71 and CVA16 vaccines could differentially affect the incidence of HFMD. Thus, the aim of our study is to compare the predicted incidence of HFMD under different vaccination scenarios in China.
Methods
We developed and analyzed a deterministic, compartmental, mathematical model of enterovirus transmission to reproduce the observed HFMD epidemic in China from 2015 to 2018. The modelled population includes one agegroup to reflect all children from age zero to 5 years of age.
Data sources
Literature review
We conducted a comprehensive review of peerreviewed and grey literature including government reports to obtain the following model parameters: transmission rate, duration of infectiousness, serotypes distribution, and vaccine efficacy.
National HFMD surveillance dataset
HFMD refer to acute infectious diseases caused by enteroviruses such as human enterovirus 71 (EVA71) and Coxsackie virus A group 16 (CVA16), which are more common in preschool children. Severe cases are mostly caused by EVA71 infection. The diagnostic criteria refer to the People’s Republic of China Hygiene Industry Standard for Hand, Foot and Mouth Disease Diagnosis (WS 588–2018). We obtained the number of reported cases of HFMD per month from 2015 to 2018 in China from the official infectious disease report from Chinese Centre for Disease Control and Prevention (Table S1) [2]. Since January 1, 2008, all probable and laboratoryconfirmed HFMD cases have been reported to Chinese Centre for Disease Control and Prevention and by May 2, 2008, HFMD was a mandatory notifiable disease (III).
Demographic data
We obtained nationallevel estimates of the total number of children under 5, the birth rate, and the under5 mortality rate from 2015 to 2018 from the National Bureau of Statistics and Provincial Bureau of Statistics [1].
Model description
We developed an extended SEIR 6/8compartment model to evaluate the spread and incidence of HFMD. The population is divided into six or eight compartments as shown in Fig. 1, based on the natural history of HFMD and intervention with vaccination. Specifically, the total population (N) was divided into four categories: susceptible to HFMD (S), exposed to HFMDrelated pathogens (E), infectious (I) and recovered (R). With the exception of the recovered compartment, the other three compartments (S, E, and I) were further divided into 2 or 3 compartments based on vaccination, type of pathogens and hospitalization: vaccinated (V); exposed (E_{s}); vaccinated and exposed (E_{v}); infectious but not hospitalized including mildly ill patients and asymptomatic (I_{N}); infectious hospitalized with pathogens other than EV71 or CVA16 (I_{2}); infectious hospitalized with EV71HFMD or CVA16HFMD (I_{3}) to capture the bivalent vaccine. Considering that the EV71 monovalent vaccine was actually introduced in 2017 in China, we should modify the meaning of I_{2} and I_{3} compartments in 2017 and 2018, where \( {\mathrm{I}}_{\mathrm{h}}^{\mathrm{e}} \) indicates people who are infectious with pathogens other than EV71 and \( {\mathrm{I}}_{\mathrm{h}}^{\mathrm{e}} \) indicates people who are infectious with EV71HFMD. The total population (N) is therefore given by N = S + V + E_{s} + E_{v} + I_{N} + I_{2} + I_{3} + R. Fig. 1 and Table 1 describe the model parameters.
The dynamics of HFMD is described in the following set of coupled ordinary differential equations to represent the dynamics of HFMD without vaccination.
The following set of eqs. (2) represent the dynamics of HFMD with the bivalent vaccine (protective against EV71 and CVA16 simultaneously) and the monovalent vaccine (protective against EV71).
Formula (1) indicates the dynamic of HFMD under no vaccination. In formula (2), the vaccines both indicate the bivalent vaccine which can prevent children from being infected by EV71 and CVA16 simultaneously and the monovalent vaccine which can only prevent children from being infected by EV71. We explicitly simulate the introduction of the EV71 vaccine in 2017 because the monovalent vaccine was introduced at that time in China (in December 2016) [14] .
Model assumptions
We assumed that:

a)
All the newborns are susceptible.

b)
The epidemic was established by 2015 with the initial number of infectious [I_{2}(0) + I_{3}(0)] on January 1, 2015 equal to the number of reported cases in January 2015 divided by 31.

c)
We assumed that a vaccinated person has received the entire vaccination procedure (two doses of vaccination) and will be protected against EV71 or CVA16 immediately after completing the full course of inoculation with EV71 or CVA16 vaccine, and there is no crossprotection effect with other virus subtypes. The possibility of dominant virus subtypes of HFMD changing after vaccination is not considered.

d)
There is no HFMDattributable mortality. According to the HFMD report of Chinese Centers for Disease Control, the number of deaths due to HFMD is extremely small, about one in ten thousand [2].

e)
Since the EV71 vaccines were launched in China in December 2016, the coverage of the vaccines is still very low. We assume that there is no impact of vaccination on the transmission rate β, and everyone was susceptible before and including 2016.

f)
We assumed contact patterns reflect mass action without heterogeneity by age nor mixing by agegroups. Thus, the contact rate is subsumed within the biological transmission probability β.

g)
We assumed I_{N} include asymptomatic
Parameter estimation and model fitting
We estimate parameters by minimizing the sum of squares (MSS) [17] using the MATLAB R2018a (version 9.4) tool fminsearch [21] which is used to perform unconstrained nonlinear minimization. Posterior parameter values were obtained when the results of fminsearch converged.
Where I_{2} indicates the number of no EV71 or CVA16 HFMD \( \Big({\mathrm{I}}_{\mathrm{h}}^{\mathrm{ec}} \)) / no EV71 HFMD \( \left({\mathrm{I}}_{\mathrm{h}}^{\mathrm{e}}\right) \), I_{3} indicates the number of EV71 or CVA16 HFMD \( \left({\mathrm{I}}_{\mathrm{h}}^{\mathrm{ec}}\right) \) / only EV71 HFMD \( \left({\mathrm{I}}_{\mathrm{h}}^{\mathrm{e}}\right) \).
In order to test the goodness of fit between our model and the observed data, we used a ChiSquare Test for the following hypotheses: Null hypothesis H_{0}: The modeled results (I_{2} + I_{3}) are equal to the observed number of HFMD cases shown in Table S1. Alternative hypothesis H_{1}: The modeled results are not equal to the observed values shown in Table S1.
Derivation of basic reproductive number (R_{0})
The basic reproductive number (R_{0}) is of a great significance in Epidemiological studies. It is the expected number of secondary cases produced, in a completely susceptible population, after introducing a typical infective individual [22]. In this study, we can compute the basic reproduction number (R_{0}) in 2015–2018 using the method of van den Driessche and assuming that each year represents an independent epidemic with a fully susceptible population on January 1 of each year [23]. See the Additional file 2 for the derivation process and results. In this model, R_{0} is related to transmission rate (β), progression rate (α), the proportion of I_{N}, I_{2}, I_{3} (P_{1}, P_{2}, P_{3}), the recovery rate (γ_{1}, γ_{2}, γ_{3}), the mortality rate (μ) and the vaccine efficacy (VE).
Sensitivity analysis
Considering the formula of R_{0}(Formula (2) in the Additional file 2, the expression of R_{0} in 2015 and 2016), the basic reproductive number is affected by β, α, P_{1}, γ_{1}, γ_{2} and γ_{3} while α is fixed to 1. Therefore, we take a derivation with respect to β, P_{1}, γ_{1}, γ_{2} and γ_{3}, respectively. Derivation results are as follows:
We further conducted a sensitivity analysis to measure the impact of these five parameters (β, P_{1}, γ_{1}, γ_{2} and γ_{3}) on the number of EV71 or CVA16 HFMD (I_{3}) cases using Vensim PLP software (version 7.3.5, Ventana Systems, Inc.)
Analyses of different vaccination type and coverage scenarios
We used the model to simulate the HFMD incidence of I_{2} and I_{3} under different vaccine coverage scenarios. To compare HFMD incidence under the monovalent EV71 vaccine and bivalent vaccine coverage scenarios, we compared scenarios via a 5 × 6 vaccination matrix as follows: a) No EV71 vaccination throughout the year; b) 10% EV71 vaccination coverage (vaccination rate: 26,300 person/day); c) 30% EV71 vaccination coverage (vaccination rate: 78,900 person/day); d) 50% EV71 vaccination coverage (vaccination rate: 131,500 person/day); e) 80% EV71 vaccination coverage (vaccination rate: 210,400 person/day); In each of the above cases, the coverage rates of bivalent vaccine coverage were 0, 10, 30, 50, 80 and 100%. We estimated the rate of vaccination for a certain vaccination coverage level to be the number of children < 5 divided by 365 and multiplied by the coverage level. We assumed vaccine protection rate was 95% and there was no loss of protective immunity following vaccination.
Analyses of pulse vaccination
We then compared pulse vaccination strategies via increasing the vaccination rate in March, April and September (before the two peaks of HFMD incidence of a year) from 26,300 person/day (10% coverage rate), 78,900 person/day (30% coverage rate), 131,500 person/day (50% coverage rate) and 210,400 person/day (80% coverage rate) to 1.2 times, 1.4 times, 1.6 times, 1.8 times and 2.0 times, respectively. The annual vaccine coverage rates under different pulse vaccination strategies mentioned above should not be more than 100%. Therefore, we just introduced this strategy into 10, 30, 50 and 80% coverage rate baselines.
Results
Model fits
The model estimates diagnosed HFMD cases in China for 2015 (n = 1,982,819), 2016 (n = 2,258,846), 2017 (n = 1,948,522) and 2018 (n = 2,398,566), which closely approximates the observed data (Fig. 2). The monthly fitting and observed results are shown in Supplementary Table 1. The model provides an adequate fit to the observed data as shown in Table 2. All the p values are greater than 0.05, therefore, we cannot reject the null hypothesis at the 5% significant level. The modeled result is a good estimate of the observed value. The model also captures the seasonality of disease incidence with two peaks of HFMD incidences in a given year: March to May and September to October (Fig. 2).
Sensitivity analysis
As we can see from the formula (4), the value of \( \frac{d{R}_0}{d\beta} \) is greater than 0 and the values of \( \frac{d{R}_0}{d\gamma 1} \), \( \frac{d{R}_0}{d\gamma 2} \) and \( \frac{d{R}_0}{d\gamma 3} \) are less than 0. But whether the value of \( \frac{d{R}_0}{d{P}_1} \) is greater than 0 remains uncertain. These five parameters all range from 5% below to above of their fitting values. The influence of of β, P_{1} and γ_{1} on the number of cases are large, while the influence of γ_{2} and γ_{3} are minimal (Fig. 3).
Impact of different vaccination type and coverage
Figures 4 and 5 and Table 3 show the impact of bivalent and monovalent vaccines on the incidence of HFMD in one year. Higher vaccine coverage led to reduced incidence but the efficiency also declined as the NNV (number needed to be vaccinated per case reduction) increased. Increased coverage of the monovalent vaccine from 0 to 80% decreased the number of cases by 797,262(49.1%) in one year. Use of a bivalent vaccine at an 80% coverage level decreased the number of HFMD cases by 828,560(51.0%) in 1 y. Figure 4 shows that when both the EV71 vaccine and the bivalent vaccine reach high coverage, incidence may be fall to zero. As shown in Figs. 4 and 5, a change in bivalent vaccine coverage has a greater impact on I_{3} than on I_{2} as expected because I_{3} reflects the number of cases of HFMD due to both serotypes. Besides, vaccine coverage also has a slight influence on the number of I_{2} (Fig. 5). This may indicate that bivalent vaccines are more effective than the sum of their individual effects.
Impact of pulse vaccination strategies
Figure 6 and Table 4 show the impact of different strategies of pulse vaccination on the incidence of HFMD. With baseline 10% coverage rate (vaccination rate is 26,300 person/day), increasing pulse vaccination intensity from 1.0 time to 2.0 times would increase the impact from 130,639 to 167,660 cases prevented, and the NNV would decrease from 73.5 to 71.9. Incidence consistently declined with increasing pulse vaccination intensity across different baseline coverage of the bivalent vaccine. However, the value of NNV showed a slight increase from 92.7 to 93.9 at high baseline coverage of 80%, which is the opposite of other coverage scenarios.
Discussion
There have been many papers on how to control and prevent HFMD from public health or statistical model perspectives. However, there is few work based on mathematical models to simulate dynamics of HFMD and to predict various burdens of HFMD under different scenarios of bivalent vaccine coverage. Our model extends prior work by including vaccination and reproduced the seasonal pattern of HFMD observed among children under 5 years old in China. The findings have important implications for our understanding of the spread of HFMD in China, and consideration of vaccination as a control measure to reduce the spread of HFMD.
R_{0} is an important indicator of disease outbreak. Its value indicates the number of newly infected persons during one infectious period after one infected person is introduced to an totally susceptible population. If R_{0} is larger than 1, the outbreak will take place. If R_{0} is under 1, infection will die out. In our study, the formula of R_{0} consists of four parts. This is consistent with CC Lai’s research [24]. The values of R_{0} in China from 2015 to 2018 are all greater than 1. Theoretically it is already over the threshold value for HFMD outbreaks in this period if the virus is introduced and transmitted among the children. There are other studies about the basic reproduction numbers of HFMD, and its numerical range is 1.014–1.742. Wang et al. [25] added the density of pathogen in the contaminated environments as a compartment in their model and they calculated the basic reproduction number R_{0} = 1.7424. Li et al. [26] constructed a twostage structured model to fit the HFMD data from 2009 to 2014 in China and used the same way with us to estimate the basic reproduction number R0. In their study, R_{0} fluctuated in the range of 1.06–1.57 in 2019–2014, which is similar to our estimated range.
We also performed a sensitivity analysis to study the effects of different parameters on the incidence of HFMD. Derivative results show that R_{0} increases as β increases and γ_{1}, γ_{2}, γ_{3} decrease. Moreover, it’s difficult to figure out whether the recovery time courses for people who are hospitalized (1/γ_{2} and 1/γ_{3}) are longer than that of people who are not hospitalized (1/γ_{1}). For one, people who are hospitalized may receive effective medical treatment which would shorten the recovery time. For another, people usually go to hospital because of the severity of the condition, which means they may need more time to recover. Therefore, the value of \( \frac{d{R}_0}{d{P}_1} \) depends on medical interventions and disease conditions. Considering the biological meanings of these parameters, the trend of R_{0} change is consistent with reality. More specifically, the sensitivity analysis curve shows that these five parameters (β, P_{1}, γ_{1}, γ_{2}, γ_{3}) influence the incidence of EV71 and CVA16 HFMD in different directions, to varying degrees. Considering the biological meanings of β and γ_{1}, β, also known as the transmission rate, is composed of two parts: Contact frequency × Infection efficiency. Infection efficiency is determined by the characteristics of different pathogens themselves, and it is difficult for us to intervene. However, we can change the value of contact frequency in a variety of ways. γ_{1} indicates the recovery rate of people who are infectious but not hospitalized, which is equal to the reciprocal of recovery time. Most HFMD is a selflimited disease. Therefore, this part of patients who are infectious but not hospitalized accounts for the vast majority of all HFMD patients. Considering that they don’t go to hospital, home care and selfpurchasing antiviral drugs are key factors influencing the value of γ_{1} that we can intervene. Inspiringly, this provides us some new ideas to control HFMD, such as isolation of patients in the early stage of dissemination, more frequent handwashing and ventilation [27] to lower the β value, or better treatment conditions to increase the γ1, γ2 and γ3 value.
Apart from constructing the SEIR  VEIR model to fit the reported HFMD data, we introduced vaccination to predict how HFMD incidence changes along with different vaccine coverage rates. It is generally believed that EV71 vaccine can protect against EV71 but not against CVA16 infection, which is also a major cause of HFMD [28]. Therefore, the monovalent vaccine may be less acceptable as a longterm public health strategy to control HFMD. Our findings suggest that bivalent vaccines may be the best strategy to reduce HFMD incidence. Currently, the bivalent vaccine for EV71 and CVA16 HFMD is under development, the bivalent vaccine will be available in the nearly future. However, the number needed to be vaccinated (NNV) per case reduction suggests that the efficiency of vaccination is low. There are some possible reasons including the low probability of suffering from HFMD among whole population, which is lower than 0.2%. Moreover, the incidence rate in susceptible population (children under 5 years old) is only 2.5% [2]. Secondly, we underestimated the effectiveness of vaccines. Diseases caused by EV71 and CVA16 are not only HFMD. We did not consider the vaccine protection effect against other diseases.
The pulse vaccination strategy consists of repeated application of vaccine at discrete time with equal interval in a population in contrast to the traditional constant vaccination [20, 29]. Pulse vaccination is considered to be a good complement to routine vaccination. Liu [16] integrated the pulse vaccination into a tuberculosis model to vaccinate the susceptible individuals periodically to overcome the shortage of the constant vaccination strategy and the effect of mixed vaccination is obvious, where the proportion of infectious individuals quickly falls down to a very low level. In our study, we introduced pulse vaccination strategy with 6 different intensities, in 4 different baselines in three specific months in view of the apparent seasonal variation of HFMD. There are two peaks of HFMD incidences in one year: May to July and September to October. From the results, we can see the decrease in NNV (number needed to be vaccinated per case reduction) value in 10, 30, and 50% coverage rates along with the pulse intensity increasing, which means as the pulse intensity increases, the efficiency of vaccine also rises up. But in 80% coverage baseline, we find the opposite result, possibly due to the inefficiency of excessive coverage itself.
Although it is somewhat unsatisfying that our analysis could not pinpoint. First, national HFMD epidemic information is calculated monthly. We cannot obtain the daily data which may help us improve the accuracy of the model. Second, serotype replacement indeed exists as the proportions of different HFMD serotypes having been changing in recent years [18], but our model in current study does not account for this transition. We assume that the proportions of different serotypes of HFMD are constant from 2015 to 2018, with reference to baseline levels [19]. Therefore, this may lead to differences from the real HFMD epidemic. Third, EV71 vaccine was launched in China at the end of 2016. But during our investigation, we find people’s acceptance of this vaccine is generally low due to its limited effectiveness and its high price. Therefore, some results are just based on our hypothesis but may not be realized in the future considering the costeffectiveness. Apart from these, we conducted a simple pulse vaccination analysis and just introduced 6 pulse intensities in 3 specific months. However, we have no idea about what intensity of the pulse or pulse vaccination in which months can achieve the maximum of costeffectiveness.
Overall, as HFMD can be caused by numerous serotypes of enteroviruses, the existing EV71 vaccine has no crossprotection against other subtypes. The epidemic of HFMD is still developing in China. Although EV71 vaccine provides us a more specific way to control the HFMD caused by EV71, which is the leading cause of severe and fatal HFMD, its preventive effect and costeffectiveness on the overall HFMD epidemic is very limited. Other researchers have found that the proportion of EV71caused HFMD decreased with the popularization of EV71 vaccination, while the proportion of HFMD caused by other serotypes such as CVA16 and CVA10 increased [18]. This suggests multivalent vaccine and more effective vaccination strategy different from routine inoculation, such as pulse vaccination, would play an important role in preventing the transmission of alltypeHFMD in the future.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
References
China Statistical Yearbook Database. http://www.stats.gov.cn/tjsj/ndsj/. Accessed 4 Apr 2019.
Online Direct Reporting of Official Infectious Diseases. http://www.nhc.gov.cn/jkj/new_index.shtml. Accessed 4 Apr 2019.
Xing W, Liao Q, Viboud C, Zhang J, Sun J, Wu JT, et al. Hand, foot, and mouth disease in China, 200812: an epidemiological study. Lancet Infect Dis. 2014;14(4):308–18. https://doi.org/10.1016/S14733099(13)703426.
Liu SL, Pan H, Liu P, Amer S, Chan TC, Zhan J, et al. Comparative epidemiology and virology of fatal and nonfatal cases of hand, foot and mouth disease in mainland China from 2008 to 2014. Rev Med Virol. 2015;25(2):115–28. https://doi.org/10.1002/rmv.1827.
Chang LY, Huang LM, Gau SS, Wu YY, Hsia SH, Fan TY, et al. Neurodevelopment and cognition in children after enterovirus 71 infection. N Engl J Med. 2007;356(12):1226–34. https://doi.org/10.1056/NEJMoa065954.
Huang MC, Wang SM, Hsu YW, Lin HC, Chi CY, Liu CC. Longterm cognitive and motor deficits after enterovirus 71 brainstem encephalitis in children. Pediatrics. 2006;118(6):e1785–8. https://doi.org/10.1542/peds.20061547.
Mao QY, Wang Y, Bian L, Xu M, Liang Z. EV71 vaccine, a new tool to control outbreaks of hand, foot and mouth disease (HFMD). Expert Rev Vaccines. 2016;15(5):599–606. https://doi.org/10.1586/14760584.2016.1138862.
Chuo F, Tiing S, Labadin J: A simple deterministic model for the spread of hand, foot and mouth disease (HFMD) in Sarawak. In: 2008 Second Asia International Conference on Modelling & Simulation (AMS): 2008: IEEE; 2008: 947–952.
Yang JY, Chen Y, Zhang FQ. Stability analysis and optimal control of a handfootmouth disease (HFMD) model. J Appl Math Comput. 2013;41(1–2):99–117. https://doi.org/10.1007/s1219001205971.
Liu JL. Threshold dynamics for a HFMD epidemic model with periodic transmission rate. Nonlinear Dynam. 2011;64(1–2):89–95. https://doi.org/10.1007/s1107101098486.
Hu YM, Wang X, Wang JZ, Wang L, Zhang YJ, Chang L, et al. Immunogenicity, safety, and lot consistency of a novel inactivated enterovirus 71 vaccine in Chinese children aged 6 to 59 months. Clin Vaccine Immunol. 2013;20(12):1805–11. https://doi.org/10.1128/CVI.0049113.
Li R, Liu L, Mo Z, Wang X, Xia J, Liang Z, et al. An inactivated enterovirus 71 vaccine in healthy children. N Engl J Med. 2014;370(9):829–37. https://doi.org/10.1056/NEJMoa1303224.
Zhu F, Xu W, Xia J, Liang Z, Liu Y, Zhang X, et al. Efficacy, safety, and immunogenicity of an enterovirus 71 vaccine in China. N Engl J Med. 2014;370(9):818–28. https://doi.org/10.1056/NEJMoa1304923.
Li L, Yin H, An Z, Feng Z. Considerations for developing an immunization strategy with enterovirus 71 vaccine. Vaccine. 2015;33(9):1107–12. https://doi.org/10.1016/j.vaccine.2014.10.081.
Yicun C, Zhiqiang K, Qingwei L, Qibin L, Zhong HJV. A combination vaccine comprising of inactivated enterovirus 71 and coxsackievirus A16 elicits balanced protective immunity against both viruses. Vaccine. 2014;32(21):2406–12.
Liu S, Li Y, Bi Y, Huang Q. Mixed vaccination strategy for the control of tuberculosis: a case study in China. Math Biosci Eng. 2017;14(3):695–708. https://doi.org/10.3934/mbe.2017039.
Zhang X, Zhao Y, Neumann AU. Partial immunity and vaccination for influenza. J Comput Biol. 2010;17(12):1689–96. https://doi.org/10.1089/cmb.2009.0003.
Li J, Pan H, Wang X, Zhu Q, Ge Y, Cai J, et al. Epidemiological surveillance of hand, foot and mouth disease in Shanghai in 2014–2016, prior to the introduction of the enterovirus 71 vaccine. Emerg Microbes Infect. 2018;7(1):1–7.
An ZJ, Liu Y, Liao QH, Zhang Y, Li KL, Yue CY, et al. Technical guide for the use of enterovirus type 71 inactivated vaccine. Chin J Vaccines Immun. 2016;22(4):458l464.
Zhou YC, Liu HW. Stability of periodic solutions for an SIS model with pulse vaccination. Math Comput Model. 2003;38(3–4):299–308. https://doi.org/10.1016/S08957177(03)900884.
Li Y, Zhang J, Zhang X. Modeling and preventive measures of hand, foot and mouth disease (HFMD) in China. Int J Environ Res Public Health. 2014;11(3):3108–17. https://doi.org/10.3390/ijerph110303108.
Wang WD, Zhao XQ. Threshold dynamics for compartmental epidemic models in periodic environments. J Dyn Differ Equ. 2008;20(3):699–717. https://doi.org/10.1007/s1088400891118.
van den Driessche P, Watmough J. Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2002;180(1–2):29–48. https://doi.org/10.1016/S00255564(02)001086.
Lai CC, Jiang DS, Wu HM, Chen HH. A dynamic model for the outbreaks of hand, foot, and mouth disease in Taiwan. Epidemiol Infect. 2016;144(7):1500–11. https://doi.org/10.1017/S0950268815002630.
Wang J, Xiao Y, Peng Z. Modelling seasonal HFMD infections with the effects of contaminated environments in mainland China. Appl Math Comput. 2016;274:615–27. https://doi.org/10.1016/j.amc.2015.11.035.
Li Y, Wang L, Pang L, Liu S. The data fitting and optimal control of a hand, foot and mouth disease (HFMD) model with stage structure. Appl Math Comput. 2016;276:61–74. https://doi.org/10.1016/j.amc.2015.11.090.
Wang J, Cao Z, Zeng DD, Wang Q, Wang X: Assessment for spatial driving forces of HFMD prevalence in Beijing, China. In: Proceedings of the Second ACM SIGSPATIALInternational Workshop on the Use of GIS in Emergency Management: 2016: ACM; 2016: 6.
Shiyang S, Liping J, Zhenglun L, Qunying M, Weiheng S, Huafei Z, et al. Evaluation of monovalent and bivalent vaccines against lethal Enterovirus 71 and Coxsackievirus A16 infection in newborn mice. Vaccine. 2014;10(10):2885–95.
Gakkhar S, Negi K. Pulse vaccination in SIRS epidemic model with nonmonotonic incidence rate. Chaos Soliton Fract. 2008;35(3):626–38. https://doi.org/10.1016/j.chaos.2006.05.054.
Acknowledgements
Not applicable.
Funding
This work was supported by the World Health Organization (WHO Reference 2018/805784–0) and National Natural Science Foundation of China (Grant No. 81673233).
Author information
Authors and Affiliations
Contributions
Z.L., J.T., Y.L. and W.W. designed the study. Y.L. and Y.L. collected the HFMD incidence data and gained eyesight on the biology and natural history of HFMD. Z.L., J.T. and J.L.H. built the model and obtained the related parameters. Y.W., Y.L. and S.M. revised and improved the model, conducted the derivation of basic reproductive number. A.L.W. revised the structure of this paper and polished the language. J.T., and Z.L. drafted the manuscript. All authors critically reviewed and approved the final version of the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Table S1.
The monthly fitting and observed cases of hand, foot and mouth disease in China from 2015 to 2018.
12879_2021_6157_MOESM2_ESM.docx
Additional file 2
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Liu, Z., Tian, J., Wang, Y. et al. The burden of hand, foot, and mouth disease among children under different vaccination scenarios in China: a dynamic modelling study. BMC Infect Dis 21, 650 (2021). https://doi.org/10.1186/s1287902106157w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1287902106157w
Keywords
 Hand, foot and mouth disease
 SEIR model
 Vaccine
 Basic reproductive number
 Pulse vaccination