Impact of the COVID-19 outbreak and interventions on hand, foot and mouth disease in Zhengzhou, China, 2014-2022: a retrospective study.

BACKGROUND: Since December 2019, COVID-19 has spread rapidly around the world, and studies have shown that measures to prevent COVID-19 can largely reduce the spread of other infectious diseases. This study explored the impact of the COVID-19 outbreak and interventions on the incidence of HFMD. METHODS: We gathered data on the prevalence of HFMD from the Children's Hospital Affiliated to Zhengzhou University. An autoregressive integrated moving average model was constructed using HFMD incidence data from 2014 to 2019, the number of cases predicted from 2020 to 2022 was predicted, and the predicted values were compared with the actual measurements. RESULTS: From January 2014 to October 2022, the Children's Hospital of Zhengzhou University admitted 103,995 children with HFMD. The average number of cases of HFMD from 2020 to 2022 was 4,946, a significant decrease from 14,859 cases from 2014 to 2019. We confirmed the best ARIMA (2,0,0) (1,1,0)12 model. From 2020 to 2022, the yearly number of cases decreased by 46.58%, 75.54%, and 66.16%, respectively, compared with the forecasted incidence. Trends in incidence across sexes and ages displayed patterns similar to those overall. CONCLUSIONS: The COVID-19 outbreak and interventions reduced the incidence of HFMD compared to that before the outbreak. Strengthening public health interventions remains a priority in the prevention of HFMD.

Epidemiological and Genetic Characterization of Coxsackievirus A6-Associated Hand, Foot, and Mouth Disease in Gwangju, South Korea, in 2022.

Coxsackievirus A6 (CV-A6) has emerged as the predominant causative agent of hand, foot, and mouth disease (HFMD) in young children. Since the declaration of coronavirus disease 2019 (COVID-19) as a global pandemic, the incidence of infectious diseases, including HFMD, has decreased markedly. When social mitigation was relaxed during the COVID-19 pandemic in 2022, the re-emergence of HFMD was observed in Gwangju, South Korea, and seasonal characteristics of the disease appeared to have changed. To investigate the molecular characteristics of enterovirus (EV) associated with HFMD during 2022, 277 specimens were collected. Children aged younger than 5 years accounted for the majority of affected individuals. EV detection and genotyping were performed using real-time RT-PCR and nested RT-PCR followed by sequence analysis. The EV detection rate was found to be 82.3%, and the main genotype identified was CV-A6. Sixteen CV-A6 samples were selected for whole genome sequencing. According to phylogenetic analysis, all CV-A6 strains from this study belonged to the sub-genotype D3 clade based on VP1 sequences. Analysis of 3D polymerase phylogeny showed that only the recombinant RF-A group was identified. In conclusion, circulating EV types should be continuously monitored to understand pathogen emergence and evolution during the post-pandemic era.

Age-time-specific transmission of hand-foot-and-mouth disease enterovirus serotypes in Vietnam: A catalytic model with maternal immunity.

Hand, foot and mouth disease (HFMD) is highly prevalent in the Asia Pacific region, particularly in Vietnam. To develop effective interventions and efficient vaccination programs, we inferred the age-time-specific transmission patterns of HFMD serotypes enterovirus A71 (EV-A71), coxsackievirus A6 (CV-A6), coxsackievirus A10 (CV-A10), coxsackievirus A16 (CV-A16) in Ho Chi Minh City, Vietnam from a case data collected during 2013-2018 and a serological survey data collected in 2015 and 2017. We proposed a catalytic model framework with good adaptability to incorporate maternal immunity using various mathematical functions. Our results indicate the high-level transmission of CV-A6 and CV-A10 which is not obvious in the case data, due to the variation of disease severity across serotypes. Our results provide statistical evidence supporting the strong association between severe illness and CV-A6 and EV-A71 infections. The HFMD dynamic pattern presents a cyclical pattern with large outbreaks followed by a decline in subsequent years. Additionally, we identify the age group with highest risk of infection as 1-2 years and emphasise the risk of future outbreaks as over 50% of children aged 6-7 years were estimated to be susceptible to CV-A16 and EV-A71. Our study highlights the importance of multivalent vaccines and active surveillance for different serotypes, supports early vaccination prior to 1 year old, and points out the potential utility for vaccinating children older than 5 years old in Vietnam.

Modeling and analysis of the effect of optimal virus control on the spread of HFMD.

A within-host and between-host hand, foot and mouth disease (HFMD) mathematical model is established and the affect of optimal control in its within-host part on HFMD transmission is studied. Through define two basic reproduction numbers, by using the fast-slow system analysis method of time scale, the global stabilities of the between-host (slow) system and within-host (fast) system are researched, respectively. An optimal control problem with drug-treatment control on coupled within-host and between-host HFMD model is formulated and analysed theoretically. Finally, the purposed optimal control measures are applied to the actual HFMD epidemic analysis in Zhejiang Province, China from April 1, 2021 to June 30, 2021. The numerical results show that the drug control strategies can reduce the virus load per capita and can effectively prevent large-scale outbreaks of HFMD.

Modification effects of long-term air pollution levels on the relationship between short-term exposure to meteorological factors and hand, foot, and mouth disease: A distributed lag non-linear model-based study in Shandong Province, China.

The occurrence of hand, foot, and mouth disease (HFMD) is closely related to meteorological factors. However, location-specific characteristics, such as persistent air pollution, may increase the complexity of the impact of meteorological factors on HFMD, and studies across different areas and populations are largely lacking. In this study, a two-stage multisite time-series analysis was conducted using data from 16 cities in Shandong Province from 2015 to 2019. In the first stage, we obtained the cumulative exposure-response curves of meteorological factors and the number of HFMD cases for each city. In the second stage, we merged the estimations from the first stage and included city-specific air pollution variables to identify significant effect modifiers and how they modified the short-term relationship between HFMD and meteorological factors. High concentrations of air pollutants may reduce the risk effects of high average temperature on HFMD and lead to a distinct peak in the cumulative exposure-response curve, while lower concentrations may increase the risk effects of high relative humidity. Furthermore, the effects of average wind speed on HFMD were different at different levels of air pollution. The differences in modification effects between subgroups were mainly manifested in the diversity and quantity of significant modifiers. The modification effects of long-term air pollution levels on the relationship between sunshine hours and HFMD may vary significantly depending on geographical location. The people in age＜3 and male groups were more susceptible to long-term air pollution. These findings contribute to a deepening understanding of the relationship between meteorological factors and HFMD and provide evidence for relevant public health decision-making.

Spatiotemporal patterns and risk mapping of provincial hand, foot, and mouth disease in mainland China, 2014-2017.

Background: Hand, foot, and mouth disease (HFMD) has remained a serious public health threat since its first outbreak in China. Analyzing the province-level spatiotemporal distribution of HFMD and mapping the relative risk in mainland China will help determine high-risk provinces and periods of infection outbreaks for use in formulating new priority areas for prevention and control of this disease. Furthermore, our study examined the effect of air pollution on HFMD nationwide, which few studies have done thus far. Methods: Data were collected on the number of provincial monthly HFMD infections, air pollution, meteorological variables, and socioeconomic variables from 2014 to 2017 in mainland China. We used spatial autocorrelation to determine the aggregate distribution of HFMD incidence. Spatiotemporal patterns of HFMD were analyzed, risk maps were developed using the Bayesian spatiotemporal model, and the impact of potential influencing factors on HFMD was assessed. Results: In our study, from 2014 to 2017, the HFMD annual incidence rate in all provinces of mainland China ranged from 138.80 to 203.15 per 100,000 people, with an average annual incidence rate of 165.86. The temporal risk of HFMD for 31 Chinese provinces exhibited cyclical and seasonal characteristics. The southern and eastern provinces had the highest spatial relative risk (RR > 3) from 2014 to 2017. The HFMD incidence risk in provinces (Hunan, Hubei, and Chongqing) located in central China increased over time. Among the meteorological variables, except for the mean two-minute wind speed (RR 0.6878; 95% CI 0.5841, 0.8042), all other variables were risk factors for HFMD. High GDP per capita (RR 0.9922; 95% CI 0.9841, 0.9999) was a protective factor against HFMD. The higher the birth rate was (RR 1.0657; 95% CI 1.0185, 1.1150), the higher the risk of HFMD. Health workers per 1,000 people (RR 1.2010; 95% CI 1.0443, 1.3771) was positively correlated with HFMD. Conclusions: From 2014 to 2017, the central provinces (Hunan, Hubei, and Chongqing) gradually became high-risk regions for HFMD. The spatiotemporal pattern of HFMD risk may be partially attributed to meteorological and socioeconomic factors. The prevalence of HFMD in the central provinces requires attention, as prevention control efforts should be strengthened there.

Genotype F of Echovirus 25 with multiple recombination pattern have been persistently and extensively circulating in Chinese mainland.

Echovirus 25 (E25), a member of the Enterovirus B (EV-B) species, can cause aseptic meningitis (AM), viral meningitis (VM), and acute flaccid paralysis (AFP). However, systematic studies on the molecular epidemiology of E25, especially those concerning its evolution and recombination, are lacking. In this study, 18 strains of E25, isolated from seven provinces of China between 2009 and 2018, were collected based on the Chinese hand, foot, and mouth disease (HFMD) surveillance network, and 95 sequences downloaded from GenBank were also screened. Based on the phylogenetic analysis of 113 full-length VP1 sequences worldwide, globally occurring E25 strains were classified into 9 genotypes (A-I), and genotype F was the dominant genotype in the Chinese mainland. The average nucleotide substitution rate of E25 was 6.08 × 10-3 substitutions/site/year, and six important transmission routes were identified worldwide. Seventeen recombination patterns were determined, of which genotype F can be divided into 9 recombination patterns. A positive selector site was found in the capsid protein region of genotype F. Recombination analysis and pressure selection analysis for genotype F showed multiple recombination patterns and evolution characteristics, which may be responsible for it being the dominant genotype in the Chinese mainland. This study provides a theoretical basis for the subsequent prevention and control of E25.

Emerging Enterovirus A71 Subgenogroup B5 Causing Severe Hand, Foot, and Mouth Disease, Vietnam, 2023.

We report on a 2023 outbreak of severe hand, foot, and mouth disease in southern Vietnam caused by an emerging lineage of enterovirus A71 subgenogroup B5. Affected children were significantly older than those reported during previous outbreaks. The virus should be closely monitored to assess its potential for global dispersal.

Atypical cutaneous findings of hand-foot-mouth disease in children: A systematic review.

INTRODUCTION: Hand-foot-mouth disease (HFMD) is a common childhood infectious disease. Atypical skin findings of HFMD, often associated with coxsackievirus A6 (CVA6), were first reported in 2008, with increasing reports worldwide since. Atypical lesions of HFMD often involve sites beyond the palms and soles and tend to have unusual, polymorphic morphology. METHODS: A systematic review was conducted on clinical features and outcomes of pediatric HFMD with atypical cutaneous manifestations. RESULTS: Eighty-five studies were included, representing 1359 cases with mean age 2.4 years and a male predominance of 61%. The most reported morphologies were vesicles (53%), papules (49%), and bullae (36%). Other morphologies included eczema herpeticum-like (19%), purpuric/petechial (7%), and Gianotti Crosti-like (4%). Common atypical sites included the arms and/or legs (47%), face (45%), and trunk (27%). CVA6 was identified in 63% of cases. Symptoms resolved in a mean of 10 days. Overall, 16% of cases received treatment, most commonly with acyclovir, intravenous antibiotics, or topical steroids. The most common complications were nail changes (21%) and desquamation (4%) which occurred a mean of 3 and 2 weeks after symptoms, respectively. CONCLUSION: Due to unusual morphologies resembling other conditions, HFMD with atypical cutaneous findings may be misdiagnosed, leading to inappropriate and unnecessary investigations, hospitalization, and treatment. Greater awareness of atypical presentations of HFMD is warranted to improve patient care and counseling on infection control precautions.

Detection of enteroviruses related to hand foot and mouth disease in wastewater of Asian communities.

Hand, foot, and mouth disease (HFMD) is contagious and predominantly affects children below the age of five. HFMD-associated serotypes of Enterovirus A (EVA) family include EVA71, Coxsackievirus A type 6 (CVA6), 10 (CVA10), and 16 (CVA16). Although prevalent in numerous Asian countries, studies on HFMD-causing agents in wastewater are scarce. This study aimed to conduct wastewater surveillance in various Asian communities to detect and quantify serotypes of EVA associated with HFMD. In total, 77 wastewater samples were collected from Indonesia, the Philippines, Thailand, and Vietnam from March 2022 to February 2023. The detection ratio for CVA6 RNA in samples from Vietnam was 40 % (8/20). The detection ratio for CVA6 and EVA71 RNA each was 25 % (5/20) for the Indonesian samples, indicating the need for clinical surveillance of CVA6, as clinical reports have been limited. For the Philippines, 12 % (2/17) of the samples were positive for CVA6 and EVA71 RNA each, with only one quantifiable sample each. Samples from Thailand had a lower detection ratio (1/20) for CVA6 RNA, and the concentration was unquantifiable. Conversely, CVA10 and CVA16 RNAs were not detected in any of the samples. The minimum and maximum concentrations of CVA6 RNA were 2.7 and 3.9 log10 copies/L and those for EVA71 RNA were 2.5 and 4.9 log10 copies/L, respectively. This study underscores the importance of wastewater surveillance in understanding the epidemiology of HFMD-associated EVA serotypes in Asian communities. Long-term wastewater surveillance is recommended to monitor changes in dominant serotypes, understand seasonality, and develop effective prevention and control strategies for HFMD.

Analysis of a SARIMA-XGBoost model for hand, foot, and mouth disease in Xinjiang, China.

Hand, foot, and mouth disease (HFMD) is a common childhood infectious disease. The incidence of HFMD has a pronounced seasonal tendency and is closely related to meteorological factors such as temperature, rainfall, and wind speed. In this paper, we propose a combined SARIMA-XGBoost model to improve the prediction accuracy of HFMD in 15 regions of Xinjiang, China. The SARIMA model is used for seasonal trends, and the XGBoost algorithm is applied for the nonlinear effects of meteorological factors. The geographical and temporal weighted regression model is designed to analyze the influence of meteorological factors from temporal and spatial perspectives. The analysis results show that the HFMD exhibits seasonal characteristics, peaking from May to August each year, and the HFMD incidence has significant spatial heterogeneity. The meteorological factors affecting the spread of HFMD vary among regions. Temperature and daylight significantly impact the transmission of the disease in most areas. Based on the verification experiment of forecasting, the proposed SARIMA-XGBoost model is superior to other models in accuracy, especially in regions with a high incidence of HFMD.

A hybrid model for hand-foot-mouth disease prediction based on ARIMA-EEMD-LSTM.

BACKGROUND: Hand, foot, and mouth disease (HFMD) is a common infectious disease that poses a serious threat to children all over the world. However, the current prediction models for HFMD still require improvement in accuracy. In this study, we proposed a hybrid model based on autoregressive integrated moving average (ARIMA), ensemble empirical mode decomposition (EEMD) and long short-term memory (LSTM) to predict the trend of HFMD. METHODS: The data used in this study was sourced from the National Clinical Research Center for Child Health and Disorders, Chongqing, China. The daily reported incidence of HFMD from 1 January 2015 to 27 July 2023 was collected to develop an ARIMA-EEMD-LSTM hybrid model. ARIMA, LSTM, ARIMA-LSTM and EEMD-LSTM models were developed to compare with the proposed hybrid model. Root mean square error (RMSE), mean absolute error (MAE) and coefficient of determination (R2) were adopted to evaluate the performances of the prediction models. RESULTS: Overall, ARIMA-EEMD-LSTM model achieved the most accurate prediction for HFMD, with RMSE, MAPE and R2 of 4.37, 2.94 and 0.996, respectively. Performing EEMD on the residual sequence yields 11 intrinsic mode functions. EEMD-LSTM model is the second best, with RMSE, MAPE and R2 of 6.20, 3.98 and 0.996. CONCLUSION: Results showed the advantage of ARIMA-EEMD-LSTM model over the ARIMA model, the LSTM model, the ARIMA-LSTM model and the EEMD-LSTM model. For the prevention and control of epidemics, the proposed hybrid model may provide a more powerful help. Compared with other three models, the two integrated with EEMD method showed significant improvement in predictive capability, offering novel insights for modeling of disease time series.

Molecular Epidemiology and Evolution of Coxsackievirus A14.

As the proportion of non-enterovirus 71 and non-coxsackievirus A16 which proportion of composition in the hand, foot, and mouth pathogenic spectrum gradually increases worldwide, the attention paid to other enteroviruses has increased. As a member of the species enterovirus A, coxsackievirus A14 (CVA14) has been epidemic around the world until now since it has been isolated. However, studies on CVA14 are poor and the effective population size, evolutionary dynamics, and recombination patterns of CVA14 are not well understood. In this study, 15 CVA14 strains were isolated from HFMD patients in mainland China from 2009 to 2019, and the complete sequences of CVA14 in GenBank as research objects were analyzed. CVA14 was divided into seven genotypes A-G based on an average nucleotide difference of the full-length VP1 coding region of more than 15%. Compared with the CVA14 prototype strain, the 15 CVA14 strains showed 84.0-84.7% nucleotide identity in the complete genome and 96.9-97.6% amino acid identity in the encoding region. Phylodynamic analysis based on 15 CVA14 strains and 22 full-length VP1 sequences in GenBank showed a mean substitution rate of 5.35 × 10-3 substitutions/site/year (95% HPD: 4.03-6.89 × 10-3) and the most recent common ancestor (tMRCA) of CVA14 dates back to 1942 (95% HPD: 1930-1950). The Bayesian skyline showed that the effective population size had experienced a decrease-increase-decrease fluctuation since 2004. The phylogeographic analysis indicated two and three possible migration paths in the world and mainland China, respectively. Four recombination patterns with others of species enterovirus A were observed in 15 CVA14 strains, among which coxsackievirus A2 (CVA2), coxsackievirus A4 (CVA4), coxsackievirus A6 (CVA6), coxsackievirus A8 (CVA8), and coxsackievirus A12 (CVA12) may act as recombinant donors in multiple regions. This study has filled the gap in the molecular epidemiological characteristics of CVA14, enriched the global CVA14 sequence database, and laid the epidemiological foundation for the future study of CVA14 worldwide.

Molecular epidemiology of hand, foot, and mouth disease in Karnataka, India in 2022.

BACKGROUND: Hand, foot, and mouth disease (HFMD) is an enteroviral disease that occurs as outbreaks and sporadic cases in India. In this study, we investigated and characterized the aetiology of HFMD cases that occurred in Karnataka, South India from April to October 2022. METHODS: Throat swabs, vesicular swabs, urine, and blood samples from suspected cases were analysed by reverse transcription polymerase chain reaction (RT-PCR) for the detection of enteroviruses. Molecular typing of the enterovirus-positive samples was carried out by amplifying the partial virion protein 1(VP1) gene sequence, followed by sequencing and phylogenetic analysis. RESULTS: Out of the 187 samples received from 82 cases, 93 (50%) tested positive (55/82 cases, 67%) for enteroviruses, with the majority of the HFMD cases reported in paediatric population of less than 5 years (36/55, 65.4%), while 3 cases (3/55, 5.4%) were adults. Out of the 55 enterovirus-positive cases, 31 showed partial VP1 region amplification and 19 of these cases were typed as coxsackievirus A16 (CV-A16) (13/19, 68.4%) and CV-A6 (6/19, 31.6%). The CV-A16 strains identified belonged to subclade B1c while two CV-A6 strains belonged to subclade E2. On molecular testing for other viruses causing fever-rash symptoms, 4/27 (15%) enterovirus-negative cases were detected as herpes simplex virus (1 case) and varicella zoster virus (3 cases) positive. CONCLUSION: The main causative agent of HFMD in Karnataka in 2022 was CV-A16, followed by CV-A6. Apart from the common paediatric HFMD cases, adult cases were also reported during this period. Further studies involving laboratory and clinical investigations are essential for monitoring and managing HFMD in the community.

Non-linear association between daily mean temperature and children's hand foot and mouth disease in Chongqing, China.

Chongqing was seriously affected by hand, foot, and mouth disease (HFMD), but the relationships between daily mean temperature and the incidence of HFMD remain unclear. This study used distributed lag nonlinear model to evaluate the effect of daily mean temperature on the incidence of HFMD in children aged < 5 years in Chongqing. Daily HFMD data from 2012 to 2019 in Chongqing were retrieved from the notifiable infectious disease surveillance system. A total of 413,476 HFMD cases aged < 5 years were reported in Chongqing from 2012 to 2019. The exposure-response curve of daily mean temperature and daily HFMD cases was wavy-shaped. The relative risks (RRs) increased as daily mean temperature below 5.66 °C or above 9.43 °C, with two peaks at 16.10 °C and 26.68 °C. The RRs reached the highest when the daily mean temperature at 26.68 °C on the current day (RR = 1.20, 95% CI 1.09-1.32), followed by the daily mean temperature at 16.10 °C at lag 5 days (RR = 1.07, 95% CI 1.05-1.08). The RRs for girls and daycare children were much higher than those for boys and scattered children, respectively. Taken together, daily mean temperature has strong effect on HFMD in children aged < 5 years old in Chongqing, particularly for girls and daycare children.

Molecular epidemiology and phylodynamic analysis of enterovirus 71 in Beijing, China, 2009-2019.

BACKGROUND: Enterovirus 71(EV71)-associated hand, foot and mouth disease (HFMD) decreased dramatically in Beijing from 2009 to 2019. This study was to investigate the epidemiological characteristics, evolutionary dynamics, geographic diffusion pathway, and other features of EV71 in Beijing, China. METHODS: We conducted a retrospective study of EV71-associated HFMD and its causative agent in Beijing, China, from 2009 to 2019. Phylogenetic and phylogeographic methods based on the EV71 genome were used to determine the evolution features, origin, and spatiotemporal dynamics. Positive selection sites in the VP1 gene were identified and exhibited in the tertiary structure. Bayesian birth-death skyline model was used to estimate the effective reproductive number (Re). RESULTS: EV71-associated HFMD decreased greatly in Beijing. From 2009 to 2019, EV71 strains prevalent in Beijing shared high homology in each gene segment and evolved with a rate of 4.99*10- 3 substitutions per site per year. The genetic diversity of EV71 first increased and peaked in 2012 and then decreased with fluctuations. The time to the most recent common ancestor (TMRCA) of EV71 in Beijing was estimated around 2003 when the EV71 strains were transmitted to Beijing from east China. Beijing played a crucial role in seeding EV71 to central China as well. Two residues (E145Q/G, A293S) under positive selection were detected from both the VP1 dataset and the P1 dataset. They were embedded within the loop of the VP1 capsid and were exposed externally. Mean Re estimate of EV71 in Beijing was about 1.007. CONCLUSION: In recent years, EV71 was not the primary causative agent of HFMD in Beijing. The low Re estimate of EV71 in Beijing implied that strategies for preventing and controlling HFMD were performed effectively. Beijing and east China played a crucial role in disseminating EV71 to other regions in China.

Molecular evolutionary dynamics of enterovirus A71, coxsackievirus A16 and coxsackievirus A6 causing hand, foot and mouth disease in Thailand, 2000-2022.

Hand, foot and mouth disease (HFMD) is a public health threat worldwide, particularly in the Asia-Pacific region. Enterovirus A71 (EV-A71), coxsackievirus A16 (CVA16), and CVA6 are the major pathogens causing HFMD outbreaks in several countries, including Thailand. We retrieved 385 VP1 nucleotide sequences, comprising 228 EV-A71, 33 CVA16, and 124 CVA6, deposited in the databases between 2000 and 2022 for molecular evolutionary characterization using Bayesian phylogeny. All EV-A71 identified belonged to genotype B, subgenotypes B4, and B5, and to genotype C, subgenotypes C1, C2, C4a, C4b, and C5. The analyzes demonstrated these viruses' co-circulation and subgenotypic changes throughout the past two decades. The CVA16 was grouped in genotype B1, predominantly subgenotype B1a, and the CVA6 was grouped in subgenotype D3, clades 1-4. The tMRCA of EV-A71 genotypes B and C, CVA16 B1, and CVA6 D3 dated 1993.79, 1982.62, 1995.86, and 2007.31, respectively, suggesting that the viruses were likely introduced and cryptically circulated in Thailand before the HFMD cases were recognized. We demonstrated these viruses' fluctuation and cyclical pattern throughout the two decades of observation. This study provided insight into evolutionary dynamics concerning molecular epidemiology and supported the selection of current genotype-matched vaccines, vaccine development, and implementation.

Dual-attention-based recurrent neural network for hand-foot-mouth disease prediction in Korea.

Hand-foot-mouth disease (HFMD) is a viral disease that occurs primarily in children. Meteorological factors have a significant impact on its popularity annually in Korea. This study proposes a new HFMD prediction model using a dual-attention-based recurrent neural network (DA-RNN) and important weather factors for HFMD in Korea. First, suspected cases of HFMD in each state were predicted using meteorological factors from the DA-RNN. Second, the weather factors were divided into six categories: temperature, wind, rainfall, day length, humidity, and air pollution to conduct sensitivity analysis. Because of this prediction, the proposed model showed the best performance in predicting the number of suspected HFMD cases in a week compared with other RNN methods. Sensitivity analysis showed that air pollution and rainfall play an important role in HFMD in Korea. This model provides information for HFMD prevention and control and can be extended to predict other infectious diseases.