This paper is in the following e-collection/theme issue:

Infectious Diseases (non-STD/STI) (1592) Public (e)Health, Digital Epidemiology and Public Health Informatics (750) Infodemiology and Infoveillance (1369) Prevention and Health Promotion (1015) Surveillance Systems (466) COVID-19 Vaccination (155) Hygiene and Infection Prevention (325)

Published on in Vol 10 (2024)

Preprints (earlier versions) of this paper are available at https://preprints.jmir.org/preprint/63146, first published .
Influence of the Enterovirus 71 Vaccine and the COVID-19 Pandemic on Hand, Foot, and Mouth Disease in China Based on Counterfactual Models: Observational Study

Influence of the Enterovirus 71 Vaccine and the COVID-19 Pandemic on Hand, Foot, and Mouth Disease in China Based on Counterfactual Models: Observational Study

Influence of the Enterovirus 71 Vaccine and the COVID-19 Pandemic on Hand, Foot, and Mouth Disease in China Based on Counterfactual Models: Observational Study

Authors of this article:

Jia Nie1, 2 Author Orcid Image ;   Tian Huang3 Author Orcid Image ;   Yuhong Sun4 Author Orcid Image ;   Zutong Peng5 Author Orcid Image ;   Wenlong Dong1 Author Orcid Image ;   Jiancheng Chen6 Author Orcid Image ;   Di Zheng7 Author Orcid Image ;   Fuyin Guo8 Author Orcid Image ;   Wenhui Shi9 Author Orcid Image ;   Yuewei Ling10 Author Orcid Image ;   Weijia Zhao5 Author Orcid Image ;   Haijun Yang11 Author Orcid Image ;   Tiejun Shui3 Author Orcid Image ;   Xiangyu Yan1 Author Orcid Image
Article Authors Cited by Tweetations Metrics

    1School of Disaster and Emergency Medicine, Tianjin University, No. 92, Weijin Road, Nankai District, Tianjin, China

    2Department of Management Science and Engineering, Stanford University, Stanford, CA, United States

    3Yan’an Hospital (affiliated with Kunming Medical University), Yunnan, China

    4Institute for Risk Assessment Sciences, Utrecht University, Utrecht, Netherlands

    5Yunnan Center for Disease Control and Prevention, Yunnan, China

    6College of Animal Science and Technology, China Agricultural University, Beijing, China

    7Department of Dermatology and Venereology, the First Affiliated Hospital of Kunming Medical University, Yunnan, China

    8Xiamen Peiyang BCI & Smart Health Innovation Research Institution, Xiamen, China

    9People's Hospital of Pu'er, Yunnan, China

    10People's Hospital of Lincang, Yunnan, China

    11Lanke Medical Technology Nanjing Research Institution, Nanjing, China

    *these authors contributed equally

    Corresponding Author:

    Xiangyu Yan, PhD


    Background: Hand, foot, and mouth disease (HFMD) is a highly contagious viral illness. Understanding the long-term trends of HFMD incidence and its epidemic characteristics under the circumstances of the enterovirus 71 (EV71) vaccination program and the outbreak of COVID-19 is crucial for effective disease surveillance and control.

    Objective: We aim to give an overview of the trends of HFMD over the past decades and evaluate the impact of the EV71 vaccination program and the COVID-19 pandemic on the epidemic trends of HFMD.

    Methods: Using official surveillance data from the Yunnan Province, China, we described long-term incidence trends and severity rates of HFMD as well as the variation of enterovirus proportions among cases. We conducted the autoregressive integrated moving average (ARIMA) of time series analyses to predict monthly incidences based on given subsets. The difference between the actual incidences and their counterfactual predictions was compared using absolute percentage errors (APEs) for periods after the EV71 vaccination program and the COVID-19 pandemic, respectively.

    Results: The annual incidence of HFMD fluctuated between 25.62 cases per 100,000 people in 2008 and 221.52 cases per 100,000 people in 2018. The incidence for men ranged from 30 to 250 cases per 100,000 people from 2008 to 2021, which was constantly higher than that for women. The annual incidence for children aged 1 to 2 years old ranged from 54.54 to 630.06 cases per 100,000 people, which was persistently higher than that for other age groups. For monthly incidences, semiannual peaks were observed for each year. All actual monthly incidences of 2014 to 2015 fell within the predicted 95% CI by the ARIMA(1,0,1)(1,1,0)[12] model. The average APE was 19% for a 2-year prediction. After the EV71 vaccination program, the actual monthly incidence of HFMD was consistently lower than the counterfactual predictions by ARIMA(1,0,1)(1,1,0)[12], with negative APEs ranging from −11% to −229% from January 2017 to April 2018. In the meantime, the proportion of EV71 among the enteroviruses causing HFMD decreased significantly, and the proportion was highly correlated (r=0.73, P=.004) with the severity rate. After the onset of the COVID-19 pandemic in 2020, the actual monthly incidence of HFMD consistently maintained a lower magnitude compared to the counterfactual predictions—ARIMA(1,0,1)(0,1,0)[12]—from February to September 2020, with considerable negative APEs (ranging from −31% to −2248%).

    Conclusions: EV71 vaccination alleviated severe HFMD cases and altered epidemiological trends. The HFMD may also benefit from nonpharmaceutical interventions during outbreaks such as the COVID-19 pandemic. Further development of a multivalent virus vaccine is crucial for effectively controlling HFMD outbreaks. Policymakers should implement nonpharmaceutical interventions and emphasize personal hygiene for routine prevention when appropriate.

    JMIR Public Health Surveill 2024;10:e63146

    doi:10.2196/63146

    Keywords

    hand, foot, and mouth disease; vaccination; enterovirus 71; COVID-19; epidemical trend; HFMD; EV71 


    Hand, foot, and mouth disease (HFMD) is a viral illness commonly seen in children under 5 years old, characterized by typical manifestations such as oral herpes and rashes on the hands and feet [1]. The incidence of HFMD has been noted in China since 2007, with multiple severe outbreaks. This contagious disease was officially integrated into the National Notifiable Disease Reporting System in 2008 [2]. Laboratory surveillance results showed that the causes of HFMD in China were enterovirus 71 (EV71) and coxsackievirus A16 (CV-A16) in most cases, and EV71 was the most frequently identified serotype among both severe and fatal cases [3]. The first inactivated EV71 whole virus vaccine for preventing severe HFMD was approved by the China Food and Drug Administration in December 2015 [2]. Most regions in China implemented the vaccine in the second half of 2016, and changes in the epidemiology of HFMD have been noted ever since [4].

    The changes in epidemiological characteristics were observed in many infectious diseases following the COVID-19 outbreak in China [5]. Under the policy of nonpharmaceutical interventions (NPIs) [6], which includes travel restrictions, contact reductions, and social distancing, diseases caused by respiratory pathogens were found to be particularly impacted by NPIs [7]. Other infectious diseases with alternative transmission modes, such as HFMD, were also influenced by the interventions [8]. Shen et al [9] explored the influence of the COVID-19 pandemic in Xi’an, using the number of HFMD cases per month from 2013 to 2019 as training datasets, and compared the actual number of cases per month with the predictions in 2020 to 2021 generated by a Bayesian structural time series model. A significant decrease (−94.52%) in the number of HFMD cases in the first half-year of 2020 was observed. Another study focusing on the effects of different NPI levels on HFMD in Guangzhou, using cases from 2015 to 2021, elaborated a significant benefit of strong NPIs on HFMD case reduction [10]. Although the implementation of the EV71 vaccination and the COVID-19 pandemic are both potential factors that could alter the epidemiological trend of HFMD, previous investigations have predominantly focused on evaluating the impact of either one of the two factors [9-12]. Therefore, it is essential to conduct quantitative evaluations of both the vaccine and the COVID-19 pandemic on the prevailing trends of HFMD epidemiology.

    Yunnan, situated in the southwestern periphery of China, covers an approximate area of 394,000 km2, with approximately 47 million inhabitants in 2021 [13]. It plays a pivotal role because of its unique location within the Greater Mekong Subregion, which is an assemblage encompassing sovereign entities including the People’s Republic of China, Cambodia, the Lao People’s Democratic Republic, Myanmar, Thailand, and Vietnam [14]. Collaborations among these regions are closely tied, including public health decisions [15]. Therefore, exploring the epidemiological characteristics of major infectious diseases represented by HFMD in the Yunnan Province is crucial for global health.

    This study aims to depict the long-term epidemiological trend of HFMD after its outbreak in the Yunnan Province, China, with quantitative evaluations of the impact of both the EV71 vaccination and the COVID-19 pandemic, using a counterfactual time series modeling approach.


    Data Sources

    Monthly HFMD surveillance data and demographic statistics from January 1, 2008, to December 31, 2021, in Yunnan were obtained from the National Notifiable Disease Reporting System established by the Chinese Center for Disease Control and Prevention [16]. Most HFMD cases identified in the system were clinically diagnosed [17] and were reported to the system within 24 hours after diagnosis. The criteria for severe cases include those with symptoms such as nervous system involvement; respiratory and circulatory dysfunction; abnormal laboratory test results, implicating increased peripheral blood leukocytes; abnormal cerebrospinal fluid; elevated blood sugar levels; and abnormalities in electroencephalogram, brain and spinal cord magnetic resonance imaging, chest x-ray, and echocardiography findings [18]. Sampling is required for all severe and fatal cases. If no severe case occurs in the county, a random pathogen test of at least 5 cases per month is required. If there are fewer than 5 cases during that month, all cases should be sampled for a pathogen test. Age, gender, and date of illness onset for each case were obtained from the system.

    Ethical Considerations

    This study was approved by the Research Ethics Committee of the Yunnan Center for Disease Control and Prevention (2023‐19). The research ethics committee waived informed consent because the study did not involve identifiable personal information.

    Descriptive Analysis

    We calculated annual HFMD incidences from 2008 to 2021 and incidences by different groups of gender (male and female) and age (<1, 1 to <2, 2 to <5, 5 to <10, 10 to <15, and ≥15 years). The severity rate was calculated as the number of severe and fatal cases divided by the total number of cases per year. The proportions of HFMD enterovirus subtypes (ie, CV-A16, EV71, and other enteroviruses) per year for the sampled cases were also calculated. The association between the severity rate and the proportion of EV71 among the enterovirus subtypes was investigated using the Spearman correlation analysis.

    Autoregressive Integrated Moving Average Modeling and Counterfactual Analysis

    Autoregressive integrated moving average (ARIMA) models (p,d,q)(P,D,Q)s were constructed to predict the subsequent monthly HFMD incidences based on given incidences [19]. The parameters p and P were the autoregressive term and seasonal autoregressive term, respectively; d and D were the respective nonseasonal difference and seasonal difference terms; q and Q were the moving average term and seasonal moving average term, respectively; and s was the periodic term. The model predicted the HFMD incidence despite the EV71 vaccination program and the occurrence of the COVID-19 pandemic, which was also called the counterfactual model. The counterfactual model predicted the incidence and analyzed how the actual incidence diverged from the predictions [20]. Specifically, using the auto.arima () function from the forecast package in R [21], the ARIMA model was built based on the training set, and an optimal set of parameters was automatically chosen by comparing the combinatorial spectrum of parameters according to the rule of the minimum Akaike information criterion or Bayesian information criterion. A Ljung–Box test was performed to check whether the residual sequence was a white noise sequence (P>.05 if the residual sequence was a white noise sequence) [22]. The absolute percent error (APE) was calculated based on the actual and predicted monthly incidence ([actual incidence−predicted incidence]/actual incidence).

    Several subsets of time series data were used for different purposes: (1) to test the predicting ability of the counterfactual models, monthly incidences of HFMD from January 2008 to December 2013 were used as the training set. The predicted incidence and the actual one from January 2014 to December 2015 were compared using APE. (2) To evaluate the influence of the EV71 vaccination program, monthly incidences of HFMD from January 2008 to December 2016 were used as the training set, and actual incidences from January 2017 to December 2018 were used for comparison. The training set was selected because the vaccination program did not start until late 2016 due to practical issues. (3) To evaluate the influence of the COVID-19 pandemic, monthly incidences from January 2017 to December 2019 were used as the training set, and incidences from January 2020 to December 2021 were compared with its counterpart. The selected ARIMA models for each training set are shown in Table S1 in Multimedia Appendix 1.

    The data were analyzed using R, version 4.2.3.


    Overview

    A total of 891,961 HFMD cases were reported in the Yunnan Province, China, from January 1, 2008, to December 31, 2021. Men accounted for 58.13% (n=518,513) of the reported cases. The proportions of cases of children aged <1, 1 to <2, 2 to <5, 5 to <10, 10 to <15, and ≥15 years old were 8.99% (n=80,188), 29.94% (n=267,702), 49.31% (n=439,798), 10.27% (n=91,583), 1.09% (n=9761), and 0.40% (n=3557), respectively. Out of all cases, 10,574 severe and 150 fatal cases were identified. The pathogen test from the samples (n=77,660) showed the proportions of CV-A16, EV71, and other kinds of enteroviruses were 31.96% (n=24,822), 27.72% (n=21,525), and 40.32% (n=31,313), respectively.

    Epidemiological Trends of HFMD

    The annual incidence fluctuated between 25.62 cases per 100,000 people in 2008 and 221.52 cases per 100,000 people in 2018. The incidence for men ranged from 30.06 to 249.89 cases per 100,000 people, which were constantly higher than that of women (20.80-190.92 cases per 100,000 people) from 2008 to 2021 (Figure 1). Although children between 2 and 5 years old accounted for more than half of the cases identified, the annual incidence for children aged 1 to 2 years old, ranging from 54.54 to 630.06 cases per 100,000 people, was persistently higher than that for other age groups (Figure 2).

    ‎
    Figure 1. Annual incidences of hand, foot, and mouth disease (per 100,000 people) by sex in the Yunnan Province, China, 2008 to 2021.
    ‎
    Figure 2. Annual incidences of hand, foot, and mouth disease (per 10,000 people) by age group (in years) in the Yunnan Province, China, 2008 to 2021.

    The monthly incidence of HFMD was characterized by 2 peaks each year, showing semiannual peaks (Figure S1 in Multimedia Appendix 1). The peaks included a major one around May, followed by a smaller one in autumn, except for 2020, during which only 1 peak was seen. Notable decreases in incidence were seen in 2017 and 2020. However, an increase was seen in 2018 to 2019.

    The severity rate and the proportion of enterovirus subtypes varied over the years (Figure 3). The highest severity rate appeared in 2011 at 405.11 per 10,000 cases, and a notable decline was observed in 2016, with a continuous decrease till the lowest point was reached in 2020 at 22.76 per 10,000 cases. The proportion of EV71 and CV-A16 among all tested enterovirus increased substantially from 2009 to 2012, reaching a peak of 98.18% (5985/6096) in 2012. During this period, the highest proportion of EV71 was found in 2009, accounting for 84% (244/289) of tested cases. On the other hand, the proportion of other enteroviruses generally decreased and remained relatively low during 2009‐2012. The proportion of EV71 and CV-A16 accounted for the major enteroviruses of the tested samples from 2013 to 2016. A turning point was seen in 2017, with a significant increase in the proportion of other types of enteroviruses, accounting for more than half of the tested cases. A continuing decline in the proportion of EV71 was observed till 2020, with it only accounting for 4.49% (279/6210) of cases in that year. In 2020 and 2021, the majority of enteroviruses were other enterovirus subtypes (excluding EV71 and CV-A16), representing 81.82% (5081/6,210) and 66.06% (3972/6013), respectively. Annual severity rates and the proportions of EV71 were highly correlated with a Spearman correlation coefficient of 0.73 (P=.004).

    ‎
    Figure 3. Annual severity rates (per 10,000 cases) of reported hand, foot, and mouth disease patients and proportions of enterovirus subtypes in the Yunnan Province, China, 2008 to 2021.

    ARIMA Modeling and Counterfactual Analysis

    Ability of Prediction

    Using actual incidences from 2008 to 2013 as the training set, the ARIMA(1,0,1)(1,1,0)[12] model predicted the incidences and 95% CIs from January 2014 to December 2015, which are shown in Figure 4A. All actual monthly incidences of 2014 to 2015 fell within the predicted 95% CI (Table S2 in Multimedia Appendix 1). The average APE was 19% for a 2-year prediction.

    ‎
    Figure 4. Actual monthly incidences (per 100,000 people) of hand, foot, and mouth disease and the predictions by autoregressive integrated moving average analysis in the Yunnan Province, China, 2008 to 2021.
    Influence of the EV71 Vaccination Program

    Using actual incidences from 2008 to 2016 as the training set, the ARIMA(1,0,1)(1,1,0)[12] model predicted incidences of HFMD for 2017 to 2018, which are are shown in Figure 4B and Table 1. The actual HFMD incidences were continuously lower than the predicted counterparts from January 2017 to April 2018. APEs exceeding –100% were identified during this period, indicating that the actual incidence was less than half of the predicted incidence. The phenomenon was more obvious in April 2017 (APE −107%), December 2017 (APE −229%), February 2018 (APE −127%), and March 2018 (APE −118%). However, an upward trend of the actual incidences was observed from May 2018 onward, surpassing the levels predicted by the model, and fluctuations in monthly incidence resurfaced subsequently. In the end, the actual incidence of HFMD reached nearly 1.5 times higher that of the predicted ones from October to December 2018.

    Table 1. Actual monthly incidences (per 100,000 people) of hand, foot, and mouth disease and the predictions by autoregressive integrated moving average analysis in the Yunnan Province, China, after the EV71 vaccination program in 2017‐2019.
    MonthActual incidencePredicted incidence95% CIAbsolute percentage error, %
    Jan 201714.7116.38(11.61 to 21.14)−11
    Feb 20177.7311.06(4.70 to 17.42)−43
    Mar 20177.4014.20(7.64 to 20.76)−92
    Apr 201710.3921.54(14.95 to 28.13)−107
    May 201715.1828.64(22.05 to 35.24)−89
    Jun 201718.6231.93(25.33 to 38.52)−71
    Jul 201719.5522.24(15.65 to 28.84)−14
    Aug 20179.3011.21(4.61 to 17.81)−21
    Sep 20178.3911.21(4.62 to 17.81)−34
    Oct 201711.8013.31(6.71 to 19.90)−13
    Nov 201710.8613.27(6.67 to 19.87)−22
    Dec 20174.4014.46(7.87 to 21.06)−229
    Jan 20188.3914.95(8.22 to 21.67)−78
    Feb 20184.4210.02(3.20 to 16.84)−127
    Mar 20185.8412.70(5.87 to 19.53)−117
    Apr 201814.6618.18(11.34 to 25.01)−24
    May 201826.4726.13(19.29 to 32.96)1
    Jun 201831.9036.24(29.40 to 43.07)−14
    Jul 201829.9327.93(21.09 to 34.77)7
    Aug 201812.4312.65(5.82 to 19.49)−2
    Sep 201813.3512.13(5.29 to 18.96)9
    Oct 201822.8416.31(9.48 to 23.15)29
    Nov 201826.4918.32(11.48 to 25.15)31
    Dec 201824.8018.65(11.82 to 25.49)25
    Influence of the COVID-19 Pandemic

    Concerning the influence of the COVID-19 pandemic, monthly incidences of HFMD from 2017 to 2019 were used as the training set. Comparisons were made between the actual monthly incidences of 2020‐2021 and the predicted values by the selected ARIMA(1,0,1)(0,1,0)[12] model (Figure 4C and Table 2). After the outbreak of COVID-19, the expected semiannual peaks of monthly incidences disappeared, evidenced by a singular peak occurring later in 2020. In addition, the actual incidence maintained a substantially diminished level compared with the prediction from February to September 2020. The actual HFMD incidences were 0.74, 0.52, and 1.33 per 100,000 people for April, May, and June 2020, respectively, whereas the corresponding predictions were 9.71, 12.21, and 23.57 per 100,000 people, respectively. More than tenfold discrepancies were observed between the actual incidences and the predictions for April, May, and June 2020 (APEs −1212%, −2248%, −1672%, respectively). The actual incidences increased following September 2020, and the semiannual peaks of monthly incidences re-emerged in 2021.

    Table 2. Actual monthly incidences (per 100,000 people) of hand, foot, and mouth disease and the predictions by autoregressive integrated moving average analysis in the Yunnan Province, China, after the COVID-19 outbreak in 2020‐2021.
    MonthActual incidencePredicted incidence95% CIAbsolute percentage error, %
    Jan 20208.529.42(0.23 to18.61)11
    Feb 20201.362.44(−13.73 to 18.61)−79
    Mar 20200.854.87(−12.98 to 22.73)−473
    Apr 20200.749.71(−8.66 to 28.07)−1212
    May 20200.5212.21(−6.31 to 30.74)−2248
    Jun 20201.3323.57(4.99 to 42.15)−1672
    Jul 20204.9932.77(14.17 to 51.37)−557
    Aug 20208.2117.58(−1.02 to 36.19)−114
    Sep 202017.0522.40(3.80 to 41.01)−31
    Oct 202036.9424.37(5.77 to 42.97)34
    Nov 202038.8319.38(0.78 to 37.98)50
    Dec 202023.5214.50(−4.10 to 33.11)38
    Jan 202114.119.41(−11.35 to 30.17)33
    Feb 20214.402.44(−22.23 to 27.10)45
    Mar 20214.824.87(−20.93 to 30.67)−1
    Apr 20216.669.70(−16.46 to 35.87)−46
    May 20219.3212.21(−14.06 to 38.49)−31
    Jun 202114.4323.57(−2.74 to 49.89)−63
    Jul 202114.6032.77(6.44 to 59.09)−124
    Aug 20215.5617.58(−8.75 to 43.91)−216
    Sep 20215.7722.40(−3.93 to 48.73)−288
    Oct 202110.2624.37(−1.96 to 50.70)−138
    Nov 202116.5019.38(−6.95 to 45.71)−17
    Dec 202115.5914.50(−11.83 to 40.84)7

    Principle Findings

    With in-depth epidemiological and time series analyses of HFMD in the Yunnan Province, we depicted long-term epidemiological trends from 2008 to 2021, with semiannual peaks in its incidence and a notable decrease after the implementation of the EV71 vaccination program in 2017 and the outbreak of COVID-19 in 2020. Moreover, the vaccination campaign reduced the number of severe HFMD cases and shifted the pathogen composition of HFMD. The counterfactual model predicted the incidence without considering the vaccination program and the COVID-19 pandemic, indicating that the vaccination program had alleviated HFMD incidence and that the NPIs during the COVID-19 pandemic may have had a significant beneficial effect on mitigating HFMD.

    The annual incidence rate ranged from 25.79 to 653.26 cases per 100,000 people in most Chinese provinces from 2015 to 2019 [23]. The incidence in the Yunnan Province fluctuated between 25.62 and 221.52 cases per 100,000 people, representing a moderate HFMD incidence level across China. We observed a characteristic semiannual pattern in HFMD incidence, with major peaks occurring around May and smaller peaks in autumn. The observed semiannual pattern was also seen in other provinces, like Jiangsu [24]. However, the incidence of HFMD in 2017 and 2020 deviated from this pattern. A significant decrease was seen in early 2017 and 2020, indicating the immediate benefits of vaccination and NPIs in mitigating HFMD transmission.

    Similar epidemiological trends regarding the pathogen composition of HFMD were also found in other regions [4]. For example, the number of cases associated with EV71 in the Guangxi Province also decreased significantly following 2017. Cases that tested positive for EV71 comprised 28% of all cases in Guangxi from 2013 to 2015, while this percentage decreased to 13.3% in 2017 to 2019. However, EV71 is still the most prevalent serotype in unvaccinated areas, such as Vietnam and Korea [25,26]. The initial application of the EV71 vaccination program resulted in a reduction in incidence; however, a subsequent resurgence of HFMD incidence was observed in this study. This phenomenon can be attributed to the vaccine’s specific targeting of EV71, which was predominant in the earlier years, leading to a decline in overall incidence. However, as the vaccine induced changes in the pathogen composition of HFMD, its effectiveness in controlling the incidence of other nontargeted enterovirus subtypes was limited, resulting in the observed resurgence in incidence. Nonetheless, the vaccination program did demonstrate significant efficacy in reducing the severity rate among HFMD cases. Although the EV71 vaccination effectively reduced the severity rate and incidence in the early stage, its capacity to control the incidence diminished over time. Thus, a more comprehensive vaccine targeting a broader spectrum of pathogens warrants consideration.

    In the early stages of the COVID-19 outbreak, the Chinese government implemented stringent NPIs, including contact reductions, social distancing, and the closure of public spaces like schools [27]. In response to these interventions, the typical semiannual peaks of HFMD incidences shifted into 1 peak per year in 2020, maintaining an extremely low level in Yunnan. The biggest difference between the actual and the predicted incidence of the ARIMA model occurred in early 2020 when the strictest lockdown was mandated by the Chinese government. The most extreme APE of –2248% was observed in May 2020. Previous reports indicated that the incidences of respiratory-transmitted diseases, vector-borne diseases, and other gastrointestinal infectious diseases decreased rapidly during the COVID-19 pandemic in China, especially during the “first-level response” phase in 2020 [28,29]. However, with the implementation of the “normalized control” and “dynamic COVID-zero” measures of less strict contact control in populations [30,31], a resurgence in the incidence of these diseases was observed. This parallels the long-term trends of HFMD observed in our study. Following the reopening of primary schools in September 2020, the incidence rate of HFMD experienced a swift rebound to the single peak in the year, before returning to the pattern similar to previous years.

    This study provides valuable insights into the long-term patterns of HFMD incidence and severity rate in the Yunnan Province. By utilizing official surveillance data and applying an ARIMA analysis, this study quantitively assessed the impact of the EV71 vaccination program and the COVID-19 pandemic on the HFMD epidemic. Our study is enhanced by the counterfactual model, which allows for a comprehensive evaluation of the long-term incidence trend of HFMD and provides us with reliable predictions. Moreover, a broad time span and the consideration of 2 possible impactors sets this study apart from previous investigations.

    Limitations

    Some limitations of this study should be acknowledged. First, this study assessed the impact of the COVID-19 pandemic on HFMD prevalence regardless of the different levels of NPIs. The influence of other public health interventions beyond vaccination and NPIs during the COVID-19 pandemic was not explicitly examined, potentially confounding the observed associations. Moreover, the ARIMA model, though useful for forecasting HFMD incidence, might not account for all potential confounding factors affecting disease dynamics.

    Conclusion

    In conclusion, this is the first study to describe the long-term epidemic trend of HFMD in the Yunnan Province. We quantitatively explored the impact of the EV71 vaccination program and the outbreak of COVID-19 on the epidemic trends of HFMD. The vaccination program effectively reduced severe HFMD cases and altered the composition of enterovirus subtypes. NPIs during the COVID-19 pandemic showed potential benefits in mitigating HFMD transmission. Public health authorities may consider the combined implementation of multivirus vaccination programs and appropriate NPIs to control HFMD outbreaks.

    Acknowledgments

    We would like to express our gratitude to the staff of the Yunnan Center for Disease Control and Prevention and health care professionals from related health care institutes for their contributions to the data collection process. This work was supported by the Tianjin University Science and Technology Innovation Leading Talent Cultivation Project (2024XQM-0024) and the Chunshui discipline construction fund (302-0704000001).

    Data Availability

    The datasets supporting the findings of this study are available from the corresponding author upon reasonable request.

    We did not use generative artificial intelligence and artificial intelligence–assisted technologies in the writing process.

    Authors' Contributions

    JN, TH, HY, TS, and XY designed the study. ZP, WD, JC, DZ, FG, WS, YL, and WZ collected the data. JN, TH, and YS performed the analysis and interpreted the results. JN, TH, HY, TS, and XY coordinated and drafted the manuscript. All authors read and approved the final manuscript.

    Conflicts of Interest

    None declared.

    Multimedia Appendix 1

    Supplementary files.

    DOCX File, 38 KB
    1. Zhu P, Ji W, Li D, et al. Current status of hand-foot-and-mouth disease. J Biomed Sci. Feb 24, 2023;30(1):15. [CrossRef] [Medline]
    2. 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. May 2016;15(5):599-606. [CrossRef] [Medline]
    3. Yang B, Liu F, Liao Q, et al. Epidemiology of hand, foot and mouth disease in China, 2008 to 2015 prior to the introduction of EV-A71 vaccine. Euro Surveill. Dec 2017;22(50):16-00824. [CrossRef] [Medline]
    4. Wang J, Jiang L, Zhang C, He W, Tan Y, Ning C. The changes in the epidemiology of hand, foot, and mouth disease after the introduction of the EV-A71 vaccine. Vaccine. Jun 8, 2021;39(25):3319-3323. [CrossRef] [Medline]
    5. Chen L, Wang L, Xing Y, et al. Persistence and variation of the indirect effects of COVID-19 restrictions on the spectrum of notifiable infectious diseases in China: analysis of national surveillance among children and adolescents from 2018 to 2021. JMIR Public Health Surveill. May 15, 2024;10:e47626. [CrossRef] [Medline]
    6. Majeed A, Quint JK, Bhatt S, Davies F, Islam N. Non-pharmaceutical interventions: evaluating challenges and priorities for future health shocks. BMJ. Oct 7, 2024;387:e080528. [CrossRef] [Medline]
    7. Fricke LM, Glöckner S, Dreier M, Lange B. Impact of non-pharmaceutical interventions targeted at COVID-19 pandemic on influenza burden - a systematic review. J Infect. Jan 2021;82(1):1-35. [CrossRef] [Medline]
    8. Geng Y, Zhang L. Impact of non-pharmaceutical interventions during COVID-19 pandemic on pertussis, scarlet fever and hand-foot-mouth disease in China. J Infect. Feb 2022;84(2):e13-e15. [CrossRef] [Medline]
    9. Shen L, Sun M, Song S, et al. The impact of anti-COVID-19 nonpharmaceutical interventions on hand, foot, and mouth disease-a spatiotemporal perspective in Xi’an, northwestern China. J Med Virol. Jul 2022;94(7):3121-3132. [CrossRef] [Medline]
    10. Wu K, Ma X, Liu H, et al. Effects of different levels of non-pharmaceutical interventions on hand, foot and mouth disease in Guangzhou, China. BMC Public Health. Dec 20, 2022;22(1):2398. [CrossRef] [Medline]
    11. He F, Rui J, Deng Z, et al. Surveillance, epidemiology and impact of EV-A71 vaccination on hand, foot, and mouth disease in Nanchang, China, 2010-2019. Front Microbiol. 2022;12:811553. [CrossRef] [Medline]
    12. Han Y, Chen Z, Zheng K, et al. Epidemiology of hand, foot, and mouth disease before and after the introduction of enterovirus 71 vaccines in Chengdu, China, 2009-2018. Pediatr Infect Dis J. Oct 2020;39(10):969-978. [CrossRef] [Medline]
    13. Government information disclosure [Webpage in Chinese]. Yunnan Provincial Bureau of Statistics. URL: http://stats.yn.gov.cn/List22.aspx [Accessed 2024-12-04]
    14. Thirty years of cooperation. Greater Mekong Subregion. URL: https://greatermekong.org [Accessed 2024-12-01]
    15. Songserm N, Woradet S, Kankarn W, et al. Cholangiocarcinoma protective factors in Greater Mekong Subregion: critical issues for joint planning to sustainably solve regional public health problems. PLoS One. 2022;17(1):e0262589. [CrossRef] [Medline]
    16. Ren X, Wang L, Cowling BJ, et al. Systematic review: National Notifiable Infectious Disease Surveillance System in China. OJPHI. 2019;11(1). [CrossRef]
    17. Li XW, Ni X, Qian SY, et al. Chinese guidelines for the diagnosis and treatment of hand, foot and mouth disease (2018 edition). World J Pediatr. Oct 2018;14(5):437-447. [CrossRef] [Medline]
    18. The Ministry of Health Office issued the hand, foot, and mouth disease diagnosis and treatment guidelines (2008 edition) [Article in Chinese]. The State Council of the People’s Republic of China. 2008. URL: https://www.gov.cn/gzdt/2008-12/12/content_1176057.htm [Accessed 2024-12-02]
    19. Sato RC. Disease management with ARIMA model in time series. Einstein (Sao Paulo). 2013;11(1):128-131. [CrossRef] [Medline]
    20. Schaffer AL, Dobbins TA, Pearson SA. Interrupted time series analysis using autoregressive integrated moving average (ARIMA) models: a guide for evaluating large-scale health interventions. BMC Med Res Methodol. Mar 22, 2021;21(1):58. [CrossRef] [Medline]
    21. Hyndman RJ, Khandakar Y. Automatic time series forecasting: the forecast package for R. J Stat Soft. 2008;27(3). [CrossRef]
    22. Yan X, Wang X, Zhang X, Wang L, Zhang B, Jia Z. The epidemic of sexually transmitted diseases under the influence of COVID-19 in China. Front Public Health. 2021;9:737817. [CrossRef] [Medline]
    23. Chen B, Yang Y, Xu X, et al. Epidemiological characteristics of hand, foot, and mouth disease in China: a meta-analysis. Medicine (Baltimore). May 21, 2021;100(20):e25930. [CrossRef] [Medline]
    24. Han Y, Ji H, Shen W, et al. Disease burden in patients with severe hand, foot, and mouth disease in Jiangsu Province: a cross-sectional study. Hum Vaccin Immunother. Nov 30, 2022;18(5):2049168. [CrossRef] [Medline]
    25. Hoa-Tran TN, Dao ATH, Nguyen AT, et al. Coxsackieviruses A6 and A16 associated with hand, foot, and mouth disease in Vietnam, 2008-2017: essential information for rational vaccine design. Vaccine. Dec 14, 2020;38(52):8273-8285. [CrossRef] [Medline]
    26. Kang HJ, Yoon Y, Lee YP, et al. A different epidemiology of enterovirus A and enterovirus B co-circulating in Korea, 2012-2019. J Pediatr Infect Dis Soc. Apr 30, 2021;10(4):398-407. [CrossRef] [Medline]
    27. Lai S, Ruktanonchai NW, Zhou L, et al. Effect of non-pharmaceutical interventions to contain COVID-19 in China. Nature. Sep 2020;585(7825):410-413. [CrossRef] [Medline]
    28. Zhao T, Liu H, Bulloch G, Jiang Z, Cao Z, Wu Z. The influence of the COVID-19 pandemic on identifying HIV/AIDS cases in China: an interrupted time series study. Lancet Reg Health West Pac. Apr 5, 2023;36:100755. [CrossRef] [Medline]
    29. Zuo Z, Yang C, Ye F, et al. Trends in respiratory diseases before and after the COVID-19 pandemic in China from 2010 to 2021. BMC Public Health. Feb 1, 2023;23(1):217. [CrossRef] [Medline]
    30. Liang WN, Liu M, Liu J, Wang YD, Wu J, Liu X. The dynamic COVID-zero strategy on prevention and control of COVID-19 in China [Article in Chinese]. Zhonghua Yi Xue Za Zhi. Jan 25, 2022;102(4):239-242. [CrossRef] [Medline]
    31. Liu J, Liu M, Liang W. The dynamic COVID-zero strategy in China. China CDC Wkly. Jan 28, 2022;4(4):74-75. [CrossRef] [Medline]

    ‎
    APE: absolute percentage error
    ARIMA: autoregressive integrated moving average
    CV-A16: coxsackievirus A16
    EV71: enterovirus 71
    HFMD: hand, foot, and mouth disease
    NPI: nonpharmaceutical intervention

    Edited by Amaryllis Mavragani; submitted 11.06.24; peer-reviewed by Bo Zhang, Peng Wang; final revised version received 20.10.24; accepted 29.10.24; published 17.12.24.

    Copyright

    © Jia Nie, Tian Huang, Yuhong Sun, Zutong Peng, Wenlong Dong, Jiancheng Chen, Di Zheng, Fuyin Guo, Wenhui Shi, Yuewei Ling, Weijia Zhao, Haijun Yang, Tiejun Shui, Xiangyu Yan. Originally published in JMIR Public Health and Surveillance (https://publichealth.jmir.org), 17.12.2024.

    This is an open-access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIR Public Health and Surveillance, is properly cited. The complete bibliographic information, a link to the original publication on https://publichealth.jmir.org, as well as this copyright and license information must be included.