Surveillance, epidemiology, and pathogen spectrum of hand, foot, and mouth disease in mainland of China from 2008 to 2017

Exploring the spatiotemporal effects of meteorological factors on hand, foot and mouth disease: a multiscale geographically and temporally weighted regression study

The influence of meteorological factors on hand, foot, and mouth disease (HFMD) is not on the same scale, it's rare for previous studies to measure and recognize the independent regression relationship between each variable in space and time scale. This study used a multiscale geographically and temporally weighted regression (MGTWR) model to explore the relationship between the incidence of HFMD and related meteorological factors in Shandong Province, China, during 2015-2019 and attempted to quantify the influence of meteorological factors on HFMD under different spatiotemporal effects. Meanwhile, we used the Global Moran's I statistic and Local Moran's I statistic to test the spatial autocorrelation of the incidence of HFMD. HFMD had spatial autocorrelation at the county level in Shandong Province. The MGTWR model outperformed the OLS and GTWR models in determining the relationship between meteorological factors and HFMD. The study highlights significant spatiotemporal non-stationarity in the relationship between meteorological factors and HFMD. Temperature was predominantly positively correlated with HFMD, especially in the peninsula region during spring and summer. Humidity exhibited a predominantly positive correlation, especially in the Shandong Peninsula. Precipitation also showed a positive correlation with HFMD, particularly in western regions and during the winter months. Wind speed had a predominantly negative correlation with HFMD in the central and southwestern regions. The results might help public health authorities set priorities for targeted prevention and control measures in different regions and weather conditions, and provide guidance for the government to rationally allocate public health resources.

Effect of EV71 Vaccination on Transmission Dynamics of Hand, Foot, and Mouth Disease and Its Epidemic Prevention Threshold

OBJECTIVE

To investigate the effect of Enterovirus A71 (EV71) vaccination on the transmissibility of different enterovirus serotypes of hand, foot, and mouth disease (HFMD) in Zhejiang, China.

METHODS

Daily surveillance data of HFMD and EV71 vaccination from August 2016 to December 2019 were collected. Epidemic periods for each HFMD type were defined, and the time-varying effective reproduction number (Rt) was estimated, which could provide more direct evidence of disease epidemics than case number. General additive models (GAMs) were employed to analyze associations between EV71 vaccination quantity and rate and HFMD transmissibility. The epidemic prevention threshold, represented by required vaccination numbers and rates, was also estimated.

RESULTS

Vaccinating every 100,000 children ≤ 5 years could lead to a decrease in the Rt of EV71-associated HFMD by 14.44% (95%CI: 6.76%, 21.42%). Additionally, a positive correlation was observed between vaccinations among children ≤ 5 years old (per 100,000) and the increased transmissibility of other HFMD types (caused by enteroviruses other than EV71 and CA16) at 1.82% (95%CI: 0.80%, 2.84%). It was estimated that an additional 362,381 vaccinations, corresponding to increased vaccine coverage to 54.51% among children ≤ 5 years could effectively prevent EV71 epidemics in Zhejiang.

CONCLUSIONS

Our findings highlight the importance of enhancing EV71 vaccine coverage for controlling the epidemic of EV71-HFMD and assisting government officials in developing strategies to prevent HFMD.

Epidemiology of Enterovirus Genotypes in Association with Human Diseases

Enteroviruses (EVs) are well-known causes of a wide range of infectious diseases in infants and young children, ranging from mild illnesses to severe conditions, depending on the virus genotypes and the host's immunity. Recent advances in molecular surveillance and genotyping tools have identified over 116 different human EV genotypes from various types of clinical samples. However, the current knowledge about most of these genotypes, except for those of well-known genotypes like EV-A71 and EV-D68, is still limited due to a lack of comprehensive EV surveillance systems. This limited information makes it difficult to understand the true burden of EV-related diseases globally. Furthermore, the specific EV genotype associated with diseases varies according to country, population group, and study period. The same genotype can exhibit different epidemiological features in different areas. By integrating the data from established EV surveillance systems in the USA, Europe, Japan, and China, in combination with other EV infection studies, we can elaborate a better understanding of the distribution of prevalent EV genotypes and the diseases associated with EV. This review analyzed the data from various EV surveillance databases and explored the EV seroprevalence and the association of specific EV genotypes with human diseases.

The Milk of Cows Immunized with Trivalent Inactivated Vaccines Provides Broad-Spectrum Passive Protection against Hand, Foot, and Mouth Disease in Neonatal Mice

Hand, foot, and mouth disease (HFMD) is a contagious viral infection predominantly affecting infants and young children, caused by multiple enteroviruses, including Enterovirus 71 (EV71), Coxsackievirus A16 (CA16), Coxsackievirus A10 (CA10), and Coxsackievirus A6 (CA6). The high pathogenicity of HFMD has garnered significant attention. Currently, there is no specific treatment or broad-spectrum preventive measure available for HFMD, and existing monovalent vaccines have limited impact on the overall incidence or prevalence of the disease. Consequently, with the emergence of new viral strains driven by vaccine pressure, there is an urgent need to develop strategies for the rapid response and control of new outbreaks. In this study, we demonstrated the broad protective effect of maternal antibodies against three types of HFMD by immunizing mother mice with a trivalent inactivated vaccine targeting EV71, CA16, and CA10, using a neonatal mouse challenge model. Based on the feasibility of maternal antibodies as a form of passive immunization to prevent HFMD, we prepared a multivalent antiviral milk by immunizing dairy cows with the trivalent inactivated vaccine to target multiple HFMD viruses. In the neonatal mouse challenge model, this immunized milk exhibited extensive passive protection against oral infections caused by the three HFMD viruses. Compared to vaccines, this strategy may offer a rapid and broadly applicable approach to providing passive immunity for the prevention of HFMD, particularly in response to the swift emergence and spread of new variants.

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.

Synergetic association between coxsackievirus A16 genotype evolution and recombinant form shifts

Coxsackievirus A16 (CVA16) is a major pathogen that causes hand, foot, and mouth disease (HFMD). The recombination form (RF) shifts and global transmission dynamics of CVA16 remain unknown. In this retrospective study, global sequences of CVA16 were retrieved from the GenBank database and analyzed using comprehensive phylogenetic inference, RF surveys, and population structure. A total of 1,663 sequences were collected, forming a 442-sequences dataset for VP1 coding region analysis and a 345-sequences dataset for RF identification. Based on the VP1 coding region used for serotyping, three genotypes (A, B, and D), two subgenotypes of genotype B (B1 and B2), and three clusters of subgenotype B1 (B1a, B1b, and B1c) were identified. Cluster B1b has dominated the global epidemics, B2 disappeared in 2000, and D is an emerging genotype dating back to August 2002. Globally, four oscillation phases of CVA16 evolution, with a peak in 2013, and three migration pathways were identified. Europe, China, and Japan have served as the seeds for the global transmission of CVA16. Based on the 3D coding region of the RFs, five clusters of RFs (RF-A to -E) were identified. The shift in RFs from RF-B and RF-C to RF-D was accompanied by a change in genotype from B2 to B1a and B1c and then to B1b. In conclusion, the evolution and population dynamics of CVA16, especially the coevolution of 3D and VP1 genes, revealed that genotype evolution and RF replacement were synergistic rather than stochastic.

Effects of behaviour change on HFMD transmission

We propose a hand, foot and mouth disease (HFMD) transmission model for children with behaviour change and imperfect quarantine. The symptomatic and quarantined states obey constant behaviour change while others follow variable behaviour change depending on the numbers of new and recent infections. The basic reproduction number R 0 of the model is defined and shown to be a threshold for disease persistence and eradication. Namely, the disease-free equilibrium is globally asymptotically stable if R 0 ≤ 1 whereas the disease persists and there is a unique endemic equilibrium otherwise. By fitting the model to weekly HFMD data of Shanghai in 2019, the reproduction number is estimated at 2.41. Sensitivity analysis for R 0 shows that avoiding contagious contacts and implementing strict quarantine are essential to lower HFMD persistence. Numerical simulations suggest that strong behaviour change not only reduces the peak size and endemic level dramatically but also impairs the role of asymptomatic transmission.

The impact of EV71 vaccination program on hand, foot and mouth disease in Zhejiang Province, China: A negative control study

Objective

To estimate the potential causal impact of Enterovirus A71 (EV71) vaccination program on the reduction of EV71-infected hand, foot, and mouth disease (HFMD) in Zhejiang Province.

Methods

We utilized the longitudinal surveillance dataset of HFMD and EV71 vaccination in Zhejiang Province during 2010-2019. We estimated vaccine efficacy using a Bayesian structured time series (BSTS) model, and employed a negative control outcome (NCO) model to detect unmeasured confounding and reveal potential causal association.

Results

We estimated that 20,132 EV71 cases (95% CI: 16,733, 23,532) were prevented by vaccination program during 2017-2019, corresponding to a reduction of 29% (95% CI: 24%, 34%). The effectiveness of vaccination increased annually, with reductions of 11% (95% CI: 6%, 16%) in 2017 and 66% (95% CI: 61%, 71%) in 2019. Children under 5 years old obtained greater benefits compared to those over 5 years. Cities with higher vaccination coverage experienced a sharper EV71 reduction compared to those with lower coverage. The NCO model detected no confounding factors in the association between vaccination and EV71 cases reduction.

Conclusions

This study suggested a potential causal effect of the EV71 vaccination, highlighting the importance of achieving higher vaccine coverage to control the HFMD.

Surveillance, epidemiology, and impact of the coronavirus disease 2019 interventions on the incidence of enterovirus infections in Nanchang, China, 2010-2022

Introduction

Pathogen spectrum of Hand, foot and mouth disease (HFMD) has substantially changed in the past decade in China. Growing evidence has indicated that anti-COVID-19 nonpharmaceutical interventions (NPIs) can support control of various infectious diseases, including intestinal diseases.

Methods

In this study, HFMD cases were enrolled from sentinel hospitals of Nanchang, Jiangxi province, and enteroviruses were genotyped using specific real time RT-PCR. We systematically characterized the epidemiology of HFMD based on the continuous molecular surveillance and estimated the impact of COVID-19 intervention on HFMD incidence using seasonal autoregressive integrated moving average (ARIMA) models.

Results

A total of 10247 HFMD cases were included during 2010-2022, of which 6121 enterovirus (EV)-positive cases (59.7%) were identified by real-time RT-PCR. Over 80% cases were associated with EV-A71 and coxsackievirus A16 (CVA16) during 2010-2012, while the type distribution significantly changed as CVA6 emerged to be dominant, accounting for 22.6%-59.6% during 2013-2022. It was observed that the prevalence patterns of EV-A71 and CVA16 were similar and both of them peaked in the second quarter and then leveled off. However, CVA6 was generally prevalent around the fourth quarter, demonstrating a staggered prevalence during 2010-2019. During the COVID-19 epidemic, the seasonal HFMD epidemic peak was restrained, and the ARIMA analysis indicated that the COVID-19 intervention had mitigated EV transmission during the first COVID-19 outbreak in early 2020. In addition, bivariate Spearman's cross-correlation coefficients were estimated for the major types CVA6, CVA16 and EV-A71. Our analyses indicated the possible existence of correlations among CVA6, CVA16 and EV-A71 prevalence in the epidemiological level.

Discussion

Taken together, the type distribution of HFMD has substantially changed over the last decade and CVA6 and CVA16 are currently the most predominant types co-circulating in Nanchang. The anti-COVID-19 NPIs significantly reduced the incidence of EV infections.

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.

Phylogenetics and phylogeographic characteristics of coxsackievirus A16 in hand foot and mouth disease and herpangina cases collected in Beijing, China from 2019 to 2021

Coxsackievirus A16 (CV-A16) is a significant pathogen responsible for causing hand foot and mouth disease (HFMD) and herpangina (HA). This study aimed to investigate the recent evolution and spread of CV-A16 by monitoring HFMD and HA cases in 29 hospitals across 16 districts in Beijing from 2019 to 2021. The first five cases of HFMD and the first five cases of HA each month in each hospital were included in the study. Real-time reverse transcription polymerase chain reaction was used to identify CV-A16, CV-A6, and EV-A71. From each district, two to four CV-A16 positive samples with a relatively long sampling time interval every month were selected for sequencing. A total of 3344 HFMD cases and 2704 HA cases were enrolled in this study, with 76.0% (2541/3344) of HFMD and 45.4% (1227/2704) of HA cases confirmed to be infected by enterovirus. Among the EV-positive samples, CV-A16 virus was detected in 33.61% (854/2541) of HFMD cases and 13.4% (165/1227) of HA cases, with the predominant cluster being B1a. Both B1a and B1b had a co-circulation of local and imported strains, with different origin time (1993 vs. 1995), different global distribution (14 countries vs. 10 countries), and different transmission centers but mainly distributed in the southern and eastern regions of Beijing. Strengthening surveillance of HFMD in southern and eastern regions will improve the prevention and control efficiency of enterovirus infections.

Prevalence of Non-Polio Enteroviruses in the Sewage of Guangzhou City, China, from 2013 to 2021

Continuous surveillance of enteroviruses (EVs) in urban domestic sewage can timely reflect the circulation of EVs in the environment and crowds, and play a predictive and early warning role in EV-related diseases. To better understand the long-term epidemiological trends of circulating EVs and EV-related diseases, we conducted a 9-year (2013 to 2021) surveillance study of non-polio EVs (NPEVs) in urban sewage in Guangzhou city, China. After concentrating and isolating the viruses from the sewage samples, NPEVs were detected and molecular typing was performed. Twenty-one different NPEV serotypes were identified. The most isolated EVs were echovirus 11 (E11), followed by coxsackievirus (CV) B5, E6, and CVB3. EV species B prevailed in sewage samples, but variations in the annual frequency of different serotypes were also observed in different seasons, due to spatial and temporal factors. E11 and E6 were detected continuously before 2017, and the number of isolates was relatively stable during the surveillance period. However, after their explosive growth in 2018 and 2019, their numbers suddenly decreased significantly. CVB3 and CVB5 had alternating trends; CVB5 was most frequently detected in 2013 to 2014 and 2017 to 2018, while CVB3 was most frequently detected in 2015 to 2016 and 2020 to 2021. Phylogenetic analysis showed that at least two different transmission chains of CVB3 and CVB5 were prevalent in Guangzhou City. Our results show that in the absence of a comprehensive and systematic EV-related disease surveillance system in China, environmental surveillance is a powerful and effective tool to strengthen and further investigate the invisible transmission of EVs in the population. IMPORTANCE This study surveilled urban sewage samples from north China for 9 years to monitor enteroviruses. Samples were collected, processed, and viral identification and molecular typing were performed. We detected 21 different non-polio enteroviruses (NPEVs) with yearly variations in prevalence and peak seasons. In addition, this study is very important for understanding the epidemiology of EVs during the COVID-19 pandemic, as the detection frequency and serotypes of EVs in sewage changed considerably around 2020. We believe that our study makes a significant contribution to the literature because our results strongly suggest that environmental surveillance is an exceptionally important tool, which can be employed to detect and monitor organisms of public health concern, which would otherwise be missed and under-reported by case-based surveillance systems alone.

Epidemiology of Hand, Foot, and Mouth Disease and Genetic Evolutionary Characteristics of Coxsackievirus A10 in Taiyuan City, Shanxi Province from 2016 to 2020

In recent years, the prevalence of hand, foot, and mouth disease (HFMD) caused by enteroviruses other than enterovirus A71 (EV-A71) and coxsackievirus A16 (CVA16) has gradually increased. The throat swab specimens of 2701 HFMD cases were tested, the VP1 regions of CVA10 RNA were amplified using RT-PCR, and phylogenetic analysis of CVA10 was performed. Children aged 1–5 years accounted for the majority (81.65%) and boys were more than girls. The positivity rates of EV-A71, CVA16, and other EVs were 15.22% (219/1439), 28.77% (414/1439), and 56.01% (806/1439), respectively. CVA10 is one of the important viruses of other EVs. A total of 52 CVA10 strains were used for phylogenetic analysis based on the VP1 region, 31 were from this study, and 21 were downloaded from GenBank. All CVA10 sequences could be assigned to seven genotypes (A, B, C, D, E, F, and G), and genotype C was further divided into C1 and C2 subtypes, only one belonged to subtype C1 and the remaining 30 belonged to C2 in this study. This study emphasized the importance of strengthening the surveillance of HFMD to understand the mechanisms of pathogen variation and evolution, and to provide a scientific basis for HFMD prevention, control, and vaccine development.

Genomic Epidemiology and Transmission Dynamics of Global Coxsackievirus B4

The aim of this study was to determine the global genetic diversity and transmission dynamics of coxsackievirus B4 (CVB4) and to propose future directions for disease surveillance. Next-generation sequencing was performed to obtain the complete genome sequence of CVB4, and the genetic diversity and transmission dynamics of CVB4 worldwide were analyzed using bioinformatics methods such as phylogenetic analysis, evolutionary dynamics, and phylogeographic analysis. Forty complete genomes of CVB4 were identified from asymptomatic infected individuals and hand, foot, and mouth disease (HFMD) patients. Frequent recombination between CVB4 and EV-B multiple serotypes in the 3D^pol region was found and formed 12 recombinant patterns (A-L). Among these, the CVB4 isolated from asymptomatic infected persons and HFMD patients belonged to lineages H and I, respectively. Transmission dynamics analysis based on the VP1 region revealed that CVB4 epidemics in countries outside China were dominated by the D genotype, whereas the E genotype was dominant in China, and both genotypes evolved at a rate of > 6.50 × 10^-3 substitutions/site/year. CVB4 spreads through the population unseen, with the risk of disease outbreaks persisting as susceptible individuals accumulate. Our findings add to publicly available CVB4 genomic sequence data and deepen our understanding of CVB4 molecular epidemiology.

Coxsackievirus A6 Infection Causes Neurogenic Pathogenesis in a Neonatal Murine Model

Coxsackievirus A6 (CVA6), a member of species A enterovirus, is associated with outbreaks of hand-foot-and-mouth disease and causes a large nationwide burden of disease. However, the molecular pathogenesis of CVA6 remains unclear. In the present study, we established a suckling Institute of Cancer Research (ICR) mouse infection model to explore the neural pathogenicity of CVA6. Five-day-old mice infected with CVA6 strain F219 showed lethargy and paralysis, and died 5 or 6 days after infection via IM injection. Cerebral edema and neuronal cell swelling were observed in the infected brain tissue, and we found that the CVA6 VP1 antigen could co-localize with GFAP-positive astrocytes in infected mouse brain using an immunofluorescence assay. CVA6 strain F219 can also infect human glioma (U251) cells. Transcriptome analysis of brain tissues from infected mice and infected U251 cells showed that significantly differentially expressed genes were enriched in antiviral and immune response and neurological system processes. These results indicate that CVA6 could cause neural pathogenesis and provide basic data for exploring the mechanism of how host-cell interactions affect viral replication and pathogenesis. Importance: Coxsackievirus A6 (CVA6) surpasses the two main pathogens, enterovirus 71 (EV-A71) and coxsackievirus A16 (CVA16), which are the leading pathogens causing HFMD in many provinces of China. In our study, CVA6 infection caused neurogenic pathogenesis in a neonatal murine model, manifesting as cerebral edema and neuronal cell swelling, CVA6 VP1 antigen could co-localize with GFAP-positive astrocytes in the infected mouse brain. Based on CVA6-infected brain tissue and U251 cell transcriptome analysis, we found upregulated antiviral and immune response-related genes such as Zbp1, Usp18, Oas2, Irf7, Ddx60, Ifit3, Ddx58, and Isg15, while the neurological system process-related genes were downregulated, including Fcrls, Ebnrb, Cdk1, and Anxa5.

Antiviral effects of the petroleum ether extract of Tournefortia sibirica L. against enterovirus 71 infection in vitro and in vivo

Enterovirus 71 (EV71) is the major cause of severe hand, foot, and mouth disease (HFMD). Compared to other HFMD pathogens, like coxsackievirus A16 (CVA16), EV71 can invade the central nervous system and cause permanent damage. At present, there are no available antivirals against EV71 for clinical treatment. Herein, multiple Chinese botanical drugs were collected, and 47 types of botanical extracts were extracted using aqueous solutions and organic solvents. Based on the cytopathic effect inhibition assay, petroleum ether extract of Tournefortia sibirica L. (PE-TS) demonstrated 97.25% and 94.75% inhibition rates for EV71 infection (at 250 μg/ml) and CVA16 infection (at 125 μg/ml), respectively, with low cytotoxicity. Preliminary mechanistic studies showed that PE-TS inhibits replication of EV71 genomic RNA and synthesis of the EV71 protein. The released extracellular EV71 progeny virus titer decreased by 3.75 lg under PE-TS treatment. Furthermore, using a newborn mouse model, PE-TS treatment protected 70% and 66.7% of mice from lethal dose EV71 intracranial challenge via administration of intraperitoneal injection at 0.4 mg/g and direct lavage at 0.8 mg/g, respectively. The chemical constituents of the PE-TS were analyzed by Gas Chromatography-Mass Spectrometer (GC-MS), and a total of 60 compounds were identified. Compound-target network analysis and molecular docking implied potential bioactive compounds and their protein targets against EV71 associated pathology. The present study identified antiviral effects of PE-TS against EV71/CVA16 infection in vitro and EV71 infection in vivo, providing a potential antiviral botanical drug extract candidate for HFMD drug development.

Epidemiology of hand, foot, and mouth disease and the genetic characteristics of Coxsackievirus A16 in Taiyuan, Shanxi, China from 2010 to 2021

Hand, foot, and mouth disease (HFMD) is a common childhood infectious disease caused by human enteroviruses (EV). This study aimed to describe the epidemiological features of HFMD and the genetic characteristics of Coxsackievirus A16 (CVA16) in Taiyuan, Shanxi, China, from 2010 to 2021. Descriptive epidemiological methods were used to analyze the time and population distribution of HFMD and the genetic characteristics of CVA16. Except being affected by the COVID-19 epidemic in 2020, HFMD epidemics were sporadic from January to March each year, and began to increase in April, with a major epidemic peak from May to August, which declined in September, followed by a secondary peak from October to December. The prevalence of EV infection was the highest in children aged one to five years (84.42%), whereas its incidence was very low in children under one year of age (5.48%). Enterovirus nucleic acid was detected by real-time reverse transcription polymerase chain reaction in 6641 clinical specimens collected from patients with HFMD from 2010 to 2021, and 4236 EV-positive specimens were detected, including 988 enterovirus A71 (EV-A71), 1488 CVA16, and 1760 other enteroviruses. CVA16 remains prevalent and has co-circulated with other EVs in Taiyuan from 2010 to 2021. A phylogenetic tree constructed based on the VP1 region showed that all CVA16 strains belonged to two different clades of the B1 genotype, B1a and B1b. They showed a nucleotide similarity of 86.5-100%, and an amino acid similarity of 96.9-100%. Overall, these findings add to the global genetic resources of CVA16, demonstrate the epidemiological characteristics of HFMD as well as the genetic features of CVA16 in Taiyuan City during 2010-2021, and provide supporting evidence for the prevention and control of HFMD.

Molecular epidemiology of enteroviruses associated with hand, foot, and mouth disease in South India from 2015 to 2017

Hand, foot, and mouth disease (HFMD) is a common childhood infection caused by human enteroviruses and is clinically characterised by fever with vesicular rash on the hands, feet, and mouth. While enterovirus A71 (EV-A71) and coxsackievirus A16 (CVA16) were the major etiological agents of HFMD in India earlier, the data on recently circulating enteroviruses associated with HFMD are sparse. Here, we describe the molecular epidemiology of enteroviruses associated with HFMD in South India from 2015 to 2017. We used archived enterovirus real-time reverse transcription (RT) PCR-positive vesicle swab and/or throat swab specimens from clinically suspected HFMD cases collected from four secondary-care hospitals in South India between July 2015 and December 2017. PCR amplification and sequencing were done based on the 5'VP1, 3'VP1, VP2, or 5´NCR regions to identify enterovirus types. Genetic diversity among enteroviruses was inferred by phylogenetic analysis. Of the 107 enterovirus RNA real-time RT-PCR-positive HFMD cases, 69 (64%) were typed as CVA6, 16 (15%) were CVA16, and one (1%) was CVA10, whereas in 21 (20%) cases, the virus was not typeable by any of the methods used in the study. The majority of HFMD cases (89, 83%) were in children less than five years old, while 11 (10.3%) were in adults. 5'VP1 yielded the maximum number of enteroviruses genotyped, and phylogenetic analysis showed that the CVA6 strains belonged to subclade D3, while the subclades of CVA16 and CVA10 were B1c and D, respectively. The predominant etiological agent of HFMD in South India during 2015-2017 was CVA6, followed by CVA16 and CVA10.

Genomic Epidemiology and Phylodynamic Analysis of Enterovirus A71 Reveal Its Transmission Dynamics in Asia

Enterovirus A71 (EV-A71) is one of the main pathogens causing hand, foot, and mouth disease (HFMD) outbreaks in Asian children under 5 years of age. In severe cases, it can cause neurological complications and be life-threatening. In this study, 200 newly sequenced EV-A71 whole-genome sequences were combined with 772 EV-A71 sequences from GenBank for large-scale analysis to investigate global EV-A71 epidemiology, phylogeny, and Bayesian phylodynamic characteristics. Based on the phylogenetic analysis of the EV-A71 3D^pol region, six new evolutionary lineages (lineages B, J, K, O, P, and Q) were found in this study, and the number of evolutionary lineages was expanded from 11 to 17. Temporal dynamics and recombination breakpoint analyses based on genotype C revealed that recombination of nonstructural protein-coding regions, including 3D^pol, is an important reason for the emergence of new lineages. The EV-A71 epidemic in the Asia-Pacific region is complex, and phylogeographic analysis found that Vietnam played a key role in the spread of subgenotypes B5 and C4. The origin of EV-A71 subgenotype C4 in China is East China, which is closely related to the prevalence of subgenotype C4 in the south and throughout China. Selection pressure analysis revealed that, in addition to VP1 amino acid residues VP1-98 and VP1-145, which are associated with EV-A71 pathogenicity, amino acid residues VP1-184 and VP1-249 were also positively selected, and their functions still need to be determined by biology and immunology. This study aimed to provide a solid theoretical basis for EV-A71-related disease surveillance and prevention, antiviral research, and vaccine development through a comprehensive analysis. IMPORTANCE EV-A71 is one of the most important pathogens causing HFMD outbreaks; however, large-scale studies of EV-A71 genomic epidemiology are currently lacking. In this study, 200 new EV-A71 whole-genome sequences were determined. Combining these with 772 EV-A71 whole-genome sequences in the GenBank database, the evolutionary and transmission characteristics of global and Asian EV-A71 were analyzed. Six new evolutionary lineages were identified in this study. We also found that recombination in nonstructural protein-coding regions, including 3D^pol, is an important cause for the emergence of new lineages. The results provided a solid theoretical basis for EV-A71-related disease surveillance and prevention, antiviral research, and vaccine development.

Molecular epidemiology and clinical features of hand, foot and mouth disease requiring hospitalization after the use of enterovirus A71 inactivated vaccine in Chengdu, China, 2017-2022: a descriptive study

Three inactivated enterovirus A71 (EV-A71) vaccines have been widely vaccinated among children in the targeted age group in mainland China since mid-2016. However, comprehensive virological surveillance of hand, foot and mouth disease (HFMD) over multiple years after the use of EV-A71 vaccines has rarely been conducted. Using long-term data extracted from the Public Health and Clinical Center of Chengdu, we described the clinical, aetiological, and epidemiological characteristics of HFMD inpatients after the use of EV-A71 vaccines from 2017 through 2022. A total of 5115 patients were selected for analysis with a male-to-female ratio of 1.63:1 and were mostly under 5 years of age (97.6%). Among these cases, 4.3% presented with severe symptoms, and 4.1% of severe cases experienced significant complications. EV-A71 was no longer the major serotype for laboratory-confirmed HFMD, responsible for 15.6% of severe cases and 1.2% of mild cases. A significant downwards trend of EV-A71 infections was observed after the use of EV-A71 vaccines (P for trend < 0.001). Coxsackievirus A6 was the predominant pathogen, accounting for 63.5% of mild cases and 36.2% of severe cases. Coxsackievirus A10 (CV-A10) and A16 were sporadically detected, and an upwards trend was observed in the proportion of CV-A10 infections. This study provides baseline molecular epidemiology for the evaluation of EV-A71 vaccination impact and potential serotype replacement based on HFMD inpatients. Additional nationwide and population-based epidemiologic and serologic studies are essential to elucidate HFMD dynamics after the use of EV-A71 vaccines, and to inform public health authorities to introduce optimized intervention strategies.

The impact of enterovirus A71 vaccination program on hand, foot, and mouth disease in Guangdong, China:a longitudinal surveillance study

Enterovirus A71 (EV71) vaccination program was introduced in 2016 in China. Based on a longitudinal surveillance dataset from 2012 to 2019 in Guangdong, China, we estimated the impact of the EV71 vaccination program on hand, foot, and mouth disease (HFMD) incidence, by using a counterfactual prediction made from synthetic control approach integrated with a Bayesian time-series model. We observed a relative reduction of 41.4% for EV71-associated HFMD cases during the post-vaccination period of 2017-2019, corresponding to 26,226 cases averted. The reduction of EV71-associated HFMD cases raised with the elevation of EV71 vaccine coverage by year. We found an indirect effect for the children aged 6-14 years who were less likely to be vaccinated. Whereas, the EV71 vaccine may not protect against non-EV71-associated HFMD. This study provides a template for ongoing public health surveillance of EV71 vaccine effectiveness with a counterfactual study design. Our results show strong evidence of the EV71 vaccination program working on reducing EV71-associated HFMD in real-world settings. The finding will benefit policy-making of EV71 vaccination and the prevention of HFMD.

Multiple molecular characteristics of circulating enterovirus types among pediatric hand, foot and mouth disease patients after EV71 vaccination campaign in Wuxi, China

The molecular properties of the circulating causative agents of hand, foot and mouth disease (HFMD) in Wuxi remain unclear, posing diagnostic and prevention challenges. Additionally, in several regions of mainland China, the EV71 immunisation drastically reduced related cases and altered the HFMD pathogen spectrum, while the precise situation in Wuxi remained unknown. To address these issues, paediatric HFMD cases diagnosed in the clinic were enrolled and anal swabs were acquired in the spring of 2019. The 5′-UTR and VP1 genes were interpreted using RT-nPCR with degenerate primers to confirm their genotypes. Following that, the entire genome sequences of each viral type were recovered, allowing for the interpretation of several molecular properties. A total of 249 clinically confirmed HFMD cases had their anal swabs taken for viral identification, from which the genome sequences of seven genotypes were recovered. Coxsackievirus A16 is the most prevalent type, followed by Coxsackievirus A6, A10, A2, A4, A5 and Echovirus 11, all of which were genetically determined for the first time in Wuxi. Phylogenetic and recombination analyses were used to evaluate their evolutionary relationships with other strains found in other regions. Noticeably, a CVA16 subtype, responsible for a large proportion of the observed cases, phylogenetically clustered within clade B1a along with some strains from other countries, was the first one to be reported in China. Furthermore, some recombination events were inferred from strains detected in sporadic cases, particularly the recombination between CVA2 and CVA5 strains. Our investigation elucidated the multiple molecular characteristics of the HFMD causal enterovirus strains in Wuxi, underlining the potential hazards associated with these circulating viral types in the population and aiding in future surveillance and prevention of this disease.

Computational Screening of Newly Designed Compounds against Coxsackievirus A16 and Enterovirus A71

Outbreaks of hand, foot, and mouth disease (HFMD) that occur worldwide are mainly caused by the Coxsackievirus-A16 (CV-A16) and Enterovirus-A71 (EV-A71). Unfortunately, neither an anti-HFMD drug nor a vaccine is currently available. Rupintrivir in phase II clinical trial candidate for rhinovirus showed highly potent antiviral activities against enteroviruses as an inhibitor for 3C protease (3Cpro). In the present study, we focused on designing 50 novel rupintrivir analogs against CV-A16 and EV-A71 3Cpro using computational tools. From their predicted binding affinities, the five compounds with functional group modifications at P1′, P2, P3, and P4 sites, namely P1′-1, P2-m3, P3-4, P4-5, and P4-19, could bind with both CV-A16 and EV-A71 3Cpro better than rupintrivir. Subsequently, these five analogs were studied by 500 ns molecular dynamics simulations. Among them, P2-m3, the derivative with meta-aminomethyl-benzyl group at the P2 site, showed the greatest potential to interact with the 3Cpro target by delivering the highest number of intermolecular hydrogen bonds and contact atoms. It formed the hydrogen bonds with L127 and K130 residues at the P2 site stronger than rupintrivir, supported by significantly lower MM/PB(GB)SA binding free energies. Elucidation of designed rupintrivir analogs in our study provides the basis for developing compounds that can be candidate compounds for further HFMD treatment.

The impact of anti‐COVID‐19 non‐pharmaceutical interventions on hand, foot, and mouth disease—a spatiotemporal perspective in Xi'an, northwestern China

Growing evidence has shown that anti‐COVID‐19 nonpharmaceutical interventions (NPIs) can support prevention and control of various infectious diseases, including intestinal diseases. However, most studies focused on the short‐term mitigating impact and neglected the dynamic impact over time. This study is aimed to investigate the dynamic impact of anti‐COVID‐19 NPIs on hand, foot, and mouth disease (HFMD) over time in Xi'an City, northwestern China. Based on the surveillance data of HFMD, meteorological and web search data, Bayesian Structural Time Series model and interrupted time series analysis were performed to quantitatively measure the impact of NPIs in sequent phases with different intensities and to predict the counterfactual number of HFMD cases. From 2013 to 2021, a total number of 172,898 HFMD cases were reported in Xi'an. In 2020, there appeared a significant decrease in HFMD incidence (−94.52%, 95% CI: −97.54% to −81.95%) in the first half of the year and the peak period shifted from June to October by a small margin of 6.74% compared to the previous years of 2013 to 2019. In 2021, the seasonality of HFMD incidence gradually returned to the bimodal temporal variation pattern with a significant average decline of 61.09%. In particular, the impact of NPIs on HFMD was more evident among young children (0–3 years), and the HFMD incidence reported in industrial areas had an unexpected increase of 51.71% in 2020 autumn and winter. Results suggested that both direct and indirect NPIs should be implemented as effective public health measures to reduce infectious disease and improve surveillance strategies, and HFMD incidence in Xi'an experienced a significant rebound to the previous seasonality after a prominent decline influenced by the anti‐COVID‐19 NPIs.

HFMD transmission changed during the COVID‐19 pandemic; The impact of anti‐COVID‐19 NPIs on HFMD varied between different populations; The responsive NPIs indeed affected the HFMD incidence at different stages with potential long‐term impact.

Surveillance, Epidemiology and Impact of EV-A71 Vaccination on Hand, Foot, and Mouth Disease in Nanchang, China, 2010-2019

After the first national-scale outbreak of Hand, foot, and mouth disease (HFMD) in China, a national surveillance network was established. Here we described the epidemiology and pathogenic profile of HFMD and the impact of EV-A71 vaccination on pathogen spectrum of enteroviruses in the southeastern Chinese city of Nanchang during 2010–2019. A total of 7,951 HFMD cases from sentinel hospitals were included, of which 4,800 EV-positive cases (60.4%) were identified by real-time RT-PCR. During 2010–2012, enterovirus 71 (EV-A71) was the main causative agent of HFMD, causing 63.1% of cases, followed by 19.3% cases associated with coxsackievirus A16 (CV-A16). Since 2013, the proportion of other enteroviruses has increased dramatically, with the sub genotype D3 strain of Coxsackievirus A6 (CV-A6) replacing the dominance of EV-A71. These genetically diverse native strains of CV-A6 have co-transmitted and co-evolved in Nanchang. Unlike EV-A71 and CV-A16, most CV-A6 infections were concentrated in autumn and winter. The incidence of EV-A71 infection negatively correlated with EV-A71 vaccination (r = −0.990, p = 0.01). And severe cases sharply declined as the promotion of EV-A71 vaccines. After 2-year implementation of EV-A71 vaccination, EV-A71 is no longer detected from the reported HFMD cases in Nanchang. In conclusion, EV-A71 vaccination changed the pattern of HFMD epidemic, and CV-A6 replaced the dominance of EV-A71 over time.

Long-term neurodevelopment outcomes of hand, foot and mouth disease inpatients infected with EV-A71 or CV-A16, a retrospective cohort study

Hand, foot and mouth disease (HFMD) is a common infectious disease in western Asia area and the full range of the long-term sequelae of HFMD remains poorly described. We conducted a retrospective hospital-based cohort study of HFMD patients with central nervous system (CNS) complications caused by EV-A71 or CV-A16 between 2010 and 2016. Patients were classified into three groups, including CNS only, autonomic nervous system (ANS) dysregulation, and cardiorespiratory failure. Neurologic examination, neurodevelopmental assessments, Magnetic Resonance Imaging (MRI) and lung function, were performed at follow up. Of the 176 patients followed up, 24 suffered CNS only, 133 ANS dysregulation, and 19 cardiorespiratory failure. Median follow-up period was 4.3 years (range [1.4-8.3]). The rate of neurological abnormalities was 25% (43 of 171) at discharge and 10% (17 of 171) at follow-up. The rates of poor outcome were significantly different between the three groups of complications in motor (28%, 38%, 71%) domain (p=0.020), but not for cognitive (20%, 24%, 35%), language (25%, 36%, 41%) and adaptive (24%, 16%, 26%) domains (p = 0.537, p = 0.551, p = 0.403). For children with ventilated during hospitalization, 41% patients (14 of 34) had an obstructive ventilatory defect, and one patient with scoliosis had mixed ventilatory dysfunction. Persistent abnormalities on brain MRI were 0% (0 of 7), 9% (2 of 23) and 57% (4 of 7) in CNS, ANS and cardiorespiratory failure group separately. Patients with HFMD may have abnormalities in neurological, motor, language, cognition, adaptive behaviour and respiratory function. Long-term follow-up programmes for children's neurodevelopmental and respiratory function may be warranted.

Using Geographically Weighted Regression to Study the Seasonal Influence of Potential Risk Factors on the Incidence of HFMD on the Chinese Mainland

Hand, foot, and mouth disease (HFMD) is an epidemic infectious disease in China. Its incidence is affected by a variety of natural environmental and socioeconomic factors, and its transmission has strong seasonal and spatial heterogeneity. To quantify the spatial relationship between the incidence of HFMD (I-HFMD) and eight potential risk factors (temperature, humidity, precipitation, wind speed, air pressure, altitude, child population density, and per capita GDP) on the Chinese mainland, we established a geographically weighted regression (GWR) model to analyze their impacts in different seasons and provinces. The GWR model successfully describes the spatial changes of the influence of potential risks, and shows greatly improved estimation performance compared with the ordinary linear regression (OLR) method. Our findings help to understand the seasonally and spatially relevant effects of natural environmental and socioeconomic factors on the I-HFMD, and can provide information to be used to develop effective prevention strategies against HFMD at different locations and in different seasons.

Coxsackievirus B4: an underestimated pathogen associated with a hand, foot, and mouth disease outbreak

In order to discover the causes of a coxsackievirus B4 (CV-B4)-associated hand, foot, and mouth disease (HFMD) outbreak and to study the evolutionary characteristics of the virus, we sequenced isolates obtained during an outbreak for comparative analysis with previously sequenced strains. Phylogenetic and evolutionary dynamics analysis was performed to examine the genetic characteristics of CV-B4 in China and worldwide. Phylogenetic analysis showed that CV-B4 originated from a common ancestor in Shandong. CV-B4 strains isolated worldwide could be classified into genotypes A-E based on the sequence of the VP1 region. All CV-B4 strains in China belonged to genotype E. The global population diversity of CV-B4 fluctuated substantially over time, and CV-B4 isolated in China accounted for a significant increase in the diversity of CV-B4. The average nucleotide substitution rate in VP1 of Chinese CV-B4 (5.20 × 10^-3 substitutions/site/year) was slightly higher than that of global CV-B4 (4.82 × 10^-3 substitutions/site/year). This study is the first to investigate the evolutionary dynamics of CV-B4 and its association with an HFMD outbreak. These findings explain both the 2011 outbreak and the global increase in CV-B4 diversity. In addition to improving our understanding of a major outbreak, these findings provide a basis for the development of surveillance strategies.

Epidemiological characteristics of hand, foot, and mouth disease in China: A meta-analysis

BACKGROUND

To study the epidemic features of hand-foot-mouth disease (HFMD) in mainland China through systematic review and meta-analysis so as to provide evidence for the future prevention and control of HFMD.

METHODS

Articles on the epidemic features of HFMD in mainland China, written in English or Chinese and released between January 1, 2015 and January 1, 2020, were searched from English literature databases including Embase, Web of Science, PubMed, Cochrane library, Google academic, and Chinese literature databases including China national knowledge infrastructure (CNKI), Wanfang, and China Biology Medicine (CBM). Papers were selected according to the inclusion and exclusion criteria, and quality scoring was performed. Meta-analysis, sensitivity analysis, and identification of publication bias were finished through STATA version 12.0 software.

RESULTS

A total of 23 articles were included in this study, the total number of cases was 377,083, of which the total number of male cases was 231,798 and the total number of female cases was 145,285, the sex ratio was about 1.6:1, and the incidence of HFMD in China was 1.61‰ (95% confidence interval [CI]: 1.21‰-1.94‰). The results of the subgroup analysis showed that the incidence of HFMD in mainland China was the highest in South China, in 2014, in 1-year-old group and in other types of enteroviruses, respectively, with the rate of 3.48‰ (95% CI: 1.22‰-5.73‰), 1.81‰ (95% CI: 1.06‰-2.57‰), 15.20‰ (95% CI: 5.00‰-25.30‰), and 1.83‰ (95% CI: 1.32‰-2.33‰), respectively. The differences among the above 4 subgroups were statistically significant (P < .05). There were no publication bias in this study, and the sensitivity analysis results suggested that the meta-analysis results were robust.

CONCLUSION

There were differences in the distribution of region, time, population, and etiology of HFMD in mainland China. Health departments should adopt key strategies and measures for key populations in key areas to prevent and control the development of HFMD, and improve the ability of pathogen detection and typing in laboratories.

A novel interspecies recombinant enterovirus (Enterovirus A120) isolated from a case of acute flaccid paralysis in China

EV-A120 is a recently identified serotype of the enterovirus A species. Only one full-length genomic sequence is currently available in GenBank, and very few studies have been conducted on EV-A120 globally. Thus, additional information and research on EV-A120 are needed to explore its genetic characteristics, phylogeny, and relationship with enteroviral disease. In this study, we report the phylogenetic characteristics of a EV-A120 strain (Q0082/XZ/CHN/2000) from Tibet, China. The amino acid sequence similarity and nucleotide sequence similarity of the full-length genomic sequence of this EV-A120 strain and the EV-A120 prototype strain were 96.3% and 79.9%, respectively, showing an evolutionary trend. Recombination analysis found intraspecies recombination in the 5′ -UTR, 2B, 2C, and 3D regions. Serum neutralization testing of the EV-A120 (Q0082) strain was also carried out. Low serum-positive rates and geometric mean titres (GMTs) indicated that the extent of EV-A120 transmission and exposure in the population was very limited compared with that in the outbreaks of EV-A71 and CV-A16 in China since 2008. The EV-A120 strain (Q0082) is non-temperature sensitive, indicating its potential to spread in the population. In summary, this study reports the full-length genomic sequence of EV-A120 and provides important information for its global molecular epidemiology.

The epidemiological characteristics of enterovirus infection before and after the use of enterovirus 71 inactivated vaccine in Kunming, China

Enterovirus A71 (EV-A71) inactivated vaccines have been widely inoculated among children in Kunming City after it was approved. However, there was a large-scale outbreak of Enteroviruses (EVs) infection in Kunming, 2018. The epidemiological characteristics of HFMD and EVs were analysed during 2008–2018, which are before and three years after EV-A71 vaccine starting to use. The changes in infection spectrum were also investigated, especially for severe HFMD in 2018. The incidence of EV-A71 decreased dramatically after the EV-A71 vaccine starting use. The proportion of non-CV-A16/EV-A71 EVs positive patients raised to 77.17–85.82%, while, EV-A71 and CV-A16 only accounted for 3.41–7.24% and 6.94–19.42% in 2017 and 2018, respectively. CV-A6 was the most important causative agent in all clinical symptoms (severe HFMD, HFMD, Herpangina and fever), accounting from 42.13% to 62.33%. EV-A71 only account for 0.36–2.05%. In severe HFMD, CV-A6 (62.33%), CV-A10 (11.64%), and CV-A16 (10.96%) were the major causative agent in 2018. EV-A71 inactivated vaccine has a good protective effect against EV-A71 and induced EVs infection spectrum changefully. EV-A71 vaccine has no or insignificant cross-protection effect on CV-A6, CV-A10, and CV-A16. Herein, developing 4-valent combined vaccines is urgently needed.

Bioinformatics-based prediction of conformational epitopes for Enterovirus A71 and Coxsackievirus A16

Enterovirus A71 (EV-A71), Coxsackievirus A16 (CV-A16) and CV-A10 are the major causative agents of hand, foot and mouth disease (HFMD). The conformational epitopes play a vital role in monitoring the antigenic evolution, predicting dominant strains and preparing vaccines. In this study, we employed a Bioinformatics-based algorithm to predict the conformational epitopes of EV-A71 and CV-A16 and compared with that of CV-A10. Prediction results revealed that the distribution patterns of conformational epitopes of EV-A71 and CV-A16 were similar to that of CV-A10 and their epitopes likewise consisted of three sites: site 1 (on the "north rim" of the canyon around the fivefold vertex), site 2 (on the "puff") and site 3 (one part was in the "knob" and the other was near the threefold vertex). The reported epitopes highly overlapped with our predicted epitopes indicating the predicted results were reliable. These data suggested that three-site distribution pattern may be the basic distribution role of epitopes on the enteroviruses capsids. Our prediction results of EV-A71 and CV-A16 can provide essential information for monitoring the antigenic evolution of enterovirus.

Antiviral peptides against Enterovirus A71 causing hand, foot and mouth disease

The emergence of new and resistant viruses is a serious global burden. Conventional antiviral therapy with small molecules has led to the development of resistant mutants. In the case of hand, foot and mouth disease (HFMD), the absence of a US-FDA approved vaccine calls for urgent need to develop an antiviral that could serve as a safe, potent and robust therapy against the neurovirulent Enterovirus A71 (EV-A71). Natural peptides such as lactoferrin, melittin and synthetic peptides such as SP40, RGDS and LVLQTM have been studied against EV-A71 and have shown promising results as potent antivirals in pre-clinical studies. Peptides are considered safe, efficacious and pose fewer chances of resistance. Poor pharmacokinetic features of peptides can be overcome by the use of chemical modifications to improve in vivo delivery particularly by oral route. The use of nanotechnology can remarkably assist in the oral delivery of peptides and enhance stability in vivo. This can greatly increase patient compliance and make it more attractive as antiviral therapy.

Early Evidence of Inactivated Enterovirus 71 Vaccine Impact Against Hand, Foot, and Mouth Disease in a Major Center of Ongoing Transmission in China, 2011-2018: A Longitudinal Surveillance Study

BACKGROUND

Enterovirus 71 (EV71) is a major causative agent of hand, foot, and mouth disease (HFMD), associated with severe manifestations of the disease. Pediatric immunization with inactivated EV71 vaccine was initiated in 2016 in the Asia-Pacific region, including China. We analyzed a time series of HFMD cases attributable to EV71, coxsackievirus A16 (CA16), and other enteroviruses in Chengdu, a major transmission center in China, to assess early impacts of immunization.

METHODS

Reported HFMD cases were obtained from China's notifiable disease surveillance system. We compared observed postvaccination incidence rates during 2017-2018 with counterfactual predictions made from a negative binomial regression and a random forest model fitted to prevaccine years (2011-2015). We fit a change point model to the full time series to evaluate whether the trend of EV71 HFMD changed following vaccination.

RESULTS

Between 2011 and 2018, 279 352 HFMD cases were reported in the study region. The average incidence rate of EV71 HFMD in 2017-2018 was 60% (95% prediction interval [PI], 41%-72%) lower than predicted in the absence of immunization, corresponding to an estimated 6911 (95% PI, 3246-11 542) EV71 cases averted over 2 years. There were 52% (95% PI, 42%-60%) fewer severe HFMD cases than predicted. However, the incidence rate of non-CA16 and non-EV71 HFMD was elevated in 2018. We identified a significant decline in the trend of EV71 HFMD 4 months into the postvaccine period.

CONCLUSIONS

We provide the first real-world evidence that programmatic vaccination against EV71 is effective against childhood HFMD and present an approach to detect early vaccine impact or intended consequences from surveillance data.

Distribution of enterovirus genotypes detected in clinical samples in Hungary, 2010-2018

This report provides the findings of a retrospective surveillance study on the emergence and circulation of enteroviruses with their associated clinical symptoms over a nine-year period detected at the National Enterovirus Reference Laboratory in Hungary between 2010-2018.Enterovirus (EV) detection and genotyping were performed directly from clinical samples. From 4,080 clinical specimens 25 EV types were identified with a median age of patients of 5 years and 68% of all cases affected children aged 10 years or younger, although infections occurred in all age-groups. In 130 cases neurological symptoms were recorded, in 123 cases the infection presented in skin related signs including hand, foot, and mouth disease (HFMD), herpangina and rash. In 2010 EV-A71 was found to cause the majority of diagnosed EV infections while in 2011 and from 2014-2018, Coxsackievirus (CV)-A6 was identified most often. Echovirus E6 accounted for the most cases in 2012 and Echovirus 30 dominated in 2013. EV-D68 was identified only in 2010 and 2013.Widespread circulation of several EV-A and EV-B viruses with occasional occurrence of EV-C and EV-D was detected. The ability of EVs to cause severe infections in sporadic cases and regular outbreaks highlight the importance of continued monitoring of circulating EV types.

Amino acid substitutions in VP2, VP1, and 2C attenuate a Coxsackievirus A16 in mice

Coxsackievirus A16 (CVA16) is one of the major etiological agents of hand, foot and mouth disease (HFMD), a common acute infectious disease affecting infants and young children. Severe symptoms of the central nervous system may develop and even lead to death. Here, a plaque-purified CVA16 strain, L731-P1 (P1), was serially passaged in Vero cells for six times and passage 6 (P6) stock became highly attenuated in newborn mice. Genomic sequencing of the P1 and P6 revealed seven nucleotide substitutions at positions 1434 (C to U), 2744 (A to G), 2747 (A to G), 3161 (G to A), 3182 (A to G), 4968 (C to U), and 6064 (C to U). Six of these substitutions resulted in amino acid changes at VP2-T161 M, VP1-N102D, VP1-T103A, VP1-E241K, VP1-T248A, and 2C-S297F, respectively. P1-based infectious cDNA was generated to further investigate these virulent determinants. Independent reverse transcription-polymerase chain reaction (RT-PCR) amplifications for mutant constructions and plaque-purification of the P6 for isolation of variants were performed to determine dominant mutations and strains more related to attenuation. The virulent P1, attenuated P6, as well as a plaque purified strain (PP) and other four recombinant mutants, were inoculated into one-day-old BALB/c mice and the 50% lethal dose of each strain was determined. Comparison of virulence among these strains indicated that amino acid changes of VP1-N102D, VP1-E241K and 2C-S297F might be associated more closely with a high level attenuation of CVA16-L731-P6 than other mutations. Identification of novel residues associated with virulence may contribute to understanding of molecular basis of virulence of CVA16 and other enteroviruses.

Neutralizing antibody against Enterovirus-A71 in Thai children: A longitudinal study from birth to age 4 years

Thailand is one of the countries in the Asia-pacific region that has been most affected by the Enterovirus-A71 (EV-A71) epidemic. An individual who is susceptible to EV-A71 may also be infected asymptomatically, thus, a serological assay is a useful tool to estimate the cumulative incidence of infection in the community and to provide guidance for vaccination scheduling. There have been several candidate EV-A71 vaccines, of which three have been approved and licensed in China. The population target for EV-A71 vaccine is children younger than three years of age. In Thailand, there are limited data available on the seroprevalence of EV-A71 neutralizing (NT) antibodies and the timing of seroconversion in children. This study aims to investigate the seroprevalence and seroconversion rate of EV-A71 NT antibody in a cohort of Thai children. Sera were collected at the King Chulalongkorn Memorial Hospital in Bangkok, Thailand from 100 children between 2015 and 2020. Maternal sera were collected on the day of delivery. Serum samples from children were collected at birth (month 0) and at 2, 7, 18, 24, 36, and 48 months of age to test for EV-A71 NT antibody titers using an enzyme-linked immunosorbent assay (ELISA)-based microneutralization test. The seroprotection rate (NT antibody ≥1:16) in children at months 0, 2, 7, 18, 24, 36, and 48 was 81.0%, 60.0%, 9.0%, 10.0%, 13.0%, 17.0%, and 37.1%, respectively. The seroprotection rate was lowest at month 7 due to waning of the maternal antibody and the immunity of children increased with increasing age. At 48 months of age, less than 40% of children were seroprotected. Children at the age of 6 months should be considered a primary target for vaccination.

Genomic epidemiology of coxsackievirus A16 in mainland of China, 2000-18

Hand, foot, and mouth disease (HFMD), which is a frequently reported and concerning disease worldwide, is a severe burden on societies globally, especially in the countries of East and Southeast Asia. Coxsackievirus A16 (CV-A16) is one of the most important causes of HFMD and a severe threat to human health, especially in children under 5 years of age. To investigate the epidemiological characteristics, spread dynamics, recombinant forms (RFs), and other features of CV-A16, we leveraged the continuous surveillance data of CV-A16-related HFMD cases collected over an 18-year period. With the advent of the EV-A71 vaccine since 2016, which targeted the EV-A71-related HFMD cases, EV-A71-related HFMD cases decreased dramatically, whereas the CV-A16-related HFMD cases showed an upward trend from 2017 to October 2019. The CV-A16 strains observed in this study were genetically related and widely distributed in the mainland of China. Our results show that three clusters (B1a-B1c) existed in the mainland of China and that the cluster of B1b dominates the diffusion of CV-A16 in China. We found that eastern China played a decisive role in seeding the diffusion of CV-A16 in China, with a more complex and variant transmission trend. Although EV-A71 vaccine was launched in China in 2016, it did not affect the genetic diversity of CV-A16, and its genetic diversity did not decline, which confirmed the epidemiological surveillance trend of CV-A16. Two discontinuous clusters (2000-13 and 2014-18) were observed in the full-length genome and arranged along the time gradient, which revealed the reason why the relative genetic diversity of CV-A16 increased and experienced more complex fluctuation model after 2014. In addition, the switch from RFs B (RF-B) and RF-C co-circulation to RF-D contributes to the prevalence of B1b cluster in China after 2008. The correlation between genotype and RFs partially explained the current prevalence of B1b. This study provides unprecedented full-length genomic sequences of CV-A16 in China, with a wider geographic distribution and a long-term time scale. The study presents valuable information about CV-A16, aimed at developing effective control strategies, as well as a call for a more robust surveillance system, especially in the Asia-Pacific region.

Global Spread of the B5 Subgenotype EV-A71 and the Phylogeographical Analysis of Chinese Migration Events

The subgenotype B5 of EV-A71 is a widely circulating subgenotype that frequently spreads across the globe. Several outbreaks have occurred in nations, such as Malaysia, Thailand, Vietnam, and Japan. Appearing first in Taiwan, China, the subgenotype has been frequently reported in mainland of China even though no outbreaks have been reported so far. The current study reconstructed the migration of the B5 subgenotype of EV-A71 in China via phylogeographical analysis. Furthermore, we investigated its population dynamics in order to draw more credible inferences. Following a dataset cleanup of B5 subgenotype of EV-A71, we detected earlier B5 subgenotypes of EV-A71 sequences that had been circulating in Malaysia and Singapore since the year 2000, which was before the 2003 outbreak that occurred in Sarawak. The Bayesian inference indicated that the most recent common ancestor of B5 subgenotype EV-A71 appeared in September, 1994 (1994.75). With respect to the overall prevalence, geographical reconstruction revealed that the B5 subgenotype EV-A71 originated singly from single-source cluster and subsequently developed several active lineages. Based on a large amount of data that was accumulated, we conclude that the appearance of the B5 subgenotype of EV-A71 in mainland of China was mainly due to multiple migrations from different origins.

A fractional system of delay differential equation with nonsingular kernels in modeling hand-foot-mouth disease

In this article, we examine a computational model to explore the prevalence of a viral infectious disease, namely hand-foot-mouth disease, which is more common in infants and children. The structure of this model consists of six sub-populations along with two delay parameters. Besides, by taking advantage of the Atangana-Baleanu fractional derivative, the ability of the model to justify different situations for the system has been improved. Discussions about the existence of the solution and its uniqueness are also included in the article. Subsequently, an effective numerical scheme has been employed to obtain several meaningful approximate solutions in various scenarios imposed on the problem. The sensitivity analysis of some existing parameters in the model has also been investigated through several numerical simulations. One of the advantages of the fractional derivative used in the model is the use of the concept of memory in maintaining the substantial properties of the understudied phenomena from the origin of time to the desired time. It seems that the tools used in this model are very powerful and can effectively simulate the expected theoretical conditions in the problem, and can also be recommended in modeling other computational models in infectious diseases.

Genetic Diversity Analysis of Coxsackievirus A8 Circulating in China and Worldwide Reveals a Highly Divergent Genotype

Coxsackievirus A8 (CV-A8) is one of the pathogens associated with hand, foot and mouth disease (HFMD) and herpangina (HA), occasionally leading to severe neurological disorders such as acute flaccid paralysis (AFP). Only one study aimed at CV-A8 has been published to date, and only 12 whole-genome sequences are publicly available. In this study, complete genome sequences from 11 CV-A8 strains isolated from HFMD patients in extensive regions from China between 2013 and 2018 were determined, and all sequences from GenBank were retrieved. A phylogenetic analysis based on a total of 34 complete VP1 sequences of CV-A8 revealed five genotypes: A, B, C, D and E. The newly emerging genotype E presented a highly phylogenetic divergence compared with the other genotypes and was composed of the majority of the strains sequenced in this study. Markov chain Monte Carlo (MCMC) analysis revealed that genotype E has been evolving for nearly a century and somehow arose in approximately 2010. The Bayesian skyline plot showed that the population size of CV-A8 has experienced three dynamic fluctuations since 2001. Amino acid residues of VP1100N, 103Y, 240T and 241V, which were embedded in the potential capsid loops of genotype E, might enhance genotype E adaption to the human hosts. The CV-A8 whole genomes displayed significant intra-genotypic genetic diversity in the non-capsid region, and a total of six recombinant lineages were detected. The Chinese viruses from genotype E might have emerged recently from recombining with European CV-A6 strains. CV-A8 is a less important HFMD pathogen, and the capsid gene diversity and non-capsid recombination variety observed in CV-A8 strains indicated that the constant generation of deleterious genomes and a constant selection pressure against these deleterious mutations is still ongoing within CV-A8 quasispecies. It is possible that CV-A8 could become an important pathogen in the HFMD spectrum in the future. Further surveillance of CV-A8 is greatly needed.

Molecular typing and characterization of a novel genotype of EV-B93 isolated from Tibet, China

EV-B93 is a novel serotype within the Enterovirus B species and is uncommon worldwide. Currently, only one full-length genomic sequence (the prototype strain) has been deposited in the GenBank database. In this study, three EV-B93 were identified, including one from an acute flaccid paralysis (AFP) patient (named 99052/XZ/CHN/1999, hereafter XZ99052) and two from healthy children (named 99096/XZ/CHN/1999 and 99167/XZ/CHN/1999, hereafter XZ99096 and XZ99167, respectively) from Tibet in 1999 during the polio eradication program. The identity between the nucleotide and amino acid sequences of the Tibet EV-B93 strain and the EV-B93 prototype strain is 83.2%–83.4% and 96.8%–96.9%, respectively. The Tibet EV-B93 strain was found to have greater nucleotide sequence identity in the P3 region to another enterovirus EV-B107 as per a phylogenetic tree analysis, which revealed that recombination occurred. Seroepidemiology data showed that EV-B93 has not produced an epidemic in Tibet and there may be susceptible individuals. The three Tibet EV-B93 strains are temperature-resistant with prognosticative virulence, suggesting the possibility of a potential large-scale outbreak of EV-B93. The analyzed EV-B93 strains enrich our knowledge about this serotype and provide valuable information on global EV-B93 molecular epidemiology. What is more, they permit the appraisal of the serotype's potential public health impact and aid in understanding the role of recombination events in the evolution of enteroviruses.

Genetic recombination in fast-spreading coxsackievirus A6 variants: a potential role in evolution and pathogenicity

Hand, foot, and mouth disease (HFMD) is a common global epidemic. From 2008 onwards, many HFMD outbreaks caused by coxsackievirus A6 (CV-A6) have been reported worldwide. Since 2013, with a dramatically increasing number of CV-A6-related HFMD cases, CV-A6 has become the predominant HFMD pathogen in mainland China. Phylogenetic analysis based on the VP1 capsid gene revealed that subtype D3 dominated the CV-A6 outbreaks. Here, we performed a large-scale (near) full-length genetic analysis of global and Chinese CV-A6 variants, including 158 newly sequenced samples collected extensively in mainland China between 2010 and 2018. During the global transmission of subtype D3 of CV-A6, the noncapsid gene continued recombining, giving rise to a series of viable recombinant hybrids designated evolutionary lineages, and each lineage displayed internal consistency in both genetic and epidemiological features. The emergence of lineage –A since 2005 has triggered CV-A6 outbreaks worldwide, with a rate of evolution estimated at 4.17 × 10⁻³ substitutions site^-1 year⁻¹ based on a large number of monophyletic open reading frame sequences, and created a series of lineages chronologically through varied noncapsid recombination events. In mainland China, lineage –A has generated another two novel widespread lineages (–J and –L) through recombination within the enterovirus A gene pool, with robust estimates of occurrence time. Lineage –A, –J, and –L infections presented dissimilar clinical manifestations, indicating that the conservation of the CV-A6 capsid gene resulted in high transmissibility, but the lineage-specific noncapsid gene might influence pathogenicity. Potentially important amino acid substitutions were further predicted among CV-A6 variants. The evolutionary phenomenon of noncapsid polymorphism within the same subtype observed in CV-A6 was uncommon in other leading HFMD pathogens; such frequent recombination happened in fast-spreading CV-A6, indicating that the recovery of deleterious genomes may still be ongoing within CV-A6 quasispecies. CV-A6-related HFMD outbreaks have caused a significant public health burden and pose a great threat to children’s health; therefore, further surveillance is greatly needed to understand the full genetic diversity of CV-A6 in mainland China.

Phylogenetic characteristics and molecular epidemiological analysis of novel enterovirus EV-B83 isolated from Tibet, China

Enterovirus B83 (EV-B83) is a new member of the enterovirus B group. Currently, there are only two full-length genomic sequences of EV-B83 in the GenBank database and few VP1 region sequences. The aetiology and epidemiology of EV-B83 is unclear. 24 stool specimens were collected from twelve AFP patients and 298 stool specimens were collected from 298 healthy children in support of polio eradication activities in Tibet in 1999. Two polioviruses (isolated by L20B cell) and one non-polio enterovirus (isolated by RD cell) were isolated from AFP patients and nine polioviruses (isolated by L20B cell) and 90 non-polio enteroviruses (isolated by RD cell) were isolated from health children. Through molecular typing, we confirmed that the six of non-polio enteroviruses belong to EV-B83. The sequence similarity between the VP1 region of the Tibet isolates and that of the EV-B83 prototype strain was 80%. The maximum-likelihood phylogenetic tree of the partial VP1 region in EV-B83 demonstrated that EV-B83 formed four genotypes globally during the evolution process. The six Tibet EV-B83 strains formed the D genotype alone. Recombination analysis of Tibet EV-B83 showed that CV-B4, CV-A9, EV-B80, and EV-B106 may act as recombinant donors in multiple regions. The serum neutralization test showed that the antibody-positive rate was 58.8% and GMT was 1:19.70, which was higher than the previously reported results of EV-B106 and EV-B80. Temperature sensitivity test results showed that the six Tibet EV-B83 strains were temperature-insensitive with stronger virulence and potential infectivity, which was consistent with the results of the serum neutralization test. This study enriched the genome-wide sequence, epidemiological characteristics, and provided basic data for the follow-up study of EV-B83.

Integration of a Kalman filter in the geographically weighted regression for modeling the transmission of hand, foot and mouth disease

Background

Hand, foot and mouth disease (HFMD) is a common infectious disease whose mechanism of transmission continues to remain a puzzle for researchers. The measurement and prediction of the HFMD incidence can be combined to improve the estimation accuracy, and provide a novel perspective to explore the spatiotemporal patterns and determinant factors of an HFMD epidemic.

Methods

In this study, we collected weekly HFMD incidence reports for a total of 138 districts in Shandong province, China, from May 2008 to March 2009. A Kalman filter was integrated with geographically weighted regression (GWR) to estimate the HFMD incidence. Spatiotemporal variation characteristics were explored and potential risk regions were identified, along with quantitatively evaluating the influence of meteorological and socioeconomic factors on the HFMD incidence.

Results

The results showed that the average error covariance of the estimated HFMD incidence by district was reduced from 0.3841 to 0.1846 compared to the measured incidence, indicating an overall improvement of over 50% in error reduction. Furthermore, three specific categories of potential risk regions of HFMD epidemics in Shandong were identified by the filter processing, with manifest filtering oscillations in the initial, local and long-term periods, respectively. Amongst meteorological and socioeconomic factors, the temperature and number of hospital beds per capita, respectively, were recognized as the dominant determinants that influence HFMD incidence variation.

Conclusions

The estimation accuracy of the HFMD incidence can be significantly improved by integrating a Kalman filter with GWR and the integration is effective for exploring spatiotemporal patterns and determinants of an HFMD epidemic. Our findings could help establish more accurate HFMD prevention and control strategies in Shandong. The present study demonstrates a novel approach to exploring spatiotemporal patterns and determinant factors of HFMD epidemics, and it can be easily extended to other regions and other infectious diseases similar to HFMD.

Antiviral activity of silymarin in comparison with baicalein against EV-A71

BACKGROUND

The hand, foot and mouth disease (HFMD) is a febrile and exanthematous childhood disease mainly caused by Enterovirus 71 (EV-A71). In severe HFMD, virulent EV-A71 strains can cause acute flaccid paralysis and cardiopulmonary edema leading to death. Currently, no FDA approved antiviral treatment or vaccine is available for EV-A71. Flavonoids such as silymarin and baicalein are known to possess in vitro antiviral properties against viruses. In this study, the cytotoxicity and antiviral activity of silymarin, baicalein and baicalin were investigated.

METHODS

The cytotoxic effects of three flavonoids towards rhabdomyosarcoma (RD) cells were first examined using cell proliferation MTS [3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium] assay. Compounds found to be non-cytotoxic in RD cells were evaluated for their in vitro antiviral properties against the EV-A71 subgenotype B4 strain 41 (5865/SIN/000009) using antiviral assays. Viral infectivity was determined by reduction of the formation of plaques in RD cells. For the measurement of RNA copy number, the real time quantitative reverse transcription PCR (qRT-PCR) was used. The most potent compound was further evaluated to determine the mode of action of inhibition by time course, virus attachment and entry assays in Vero cells.

RESULTS

Silymarin was shown to exert direct extracellular virucidal effects against EV-A71 at 50% inhibitory concentration (IC₅₀) of 15.2 ± 3.53 μg/mL with SI of 10.53. Similarly, baicalein exhibited direct extracellular virucidal effects against EV-A71 at a higher IC₅₀ value of 30.88 ± 5.50 μg/mL with SI of 13.64. Besides virucidal activity, silymarin was shown to block both viral attachment and entry of EV-A71 to inhibit infection in Vero cells.

CONCLUSIONS

Silymarin has a stronger inhibition activity against EV-A71 in comparison to baicalein. It could serve as a promising antiviral drug to treat EV-A71 infections.

Pathogen Spectrum of Hand, Foot, and Mouth Disease Based on Laboratory Surveillance - China, 2018

WHAT IS ALREADY KNOWN ABOUT THIS TOPIC?

Enterovirus 71 (EV-A71) is the main causative pathogen for severe and fatal patients with Hand, Foot, and Mouth Disease (HFMD) in mainland China from 2008 to 2017. Non-EV-A71 and non-CV-A16 (other enterovirus) serotypes were the major causative-serotypes for mild HFMD in years of 2013, 2015, and 2017.

WHAT IS ADDED BY THIS REPORT?

In 2018, other enterovirus serotypes replaced EV-A71 for the first time as the major cause of severe HFMD with a proportion of 70.7%. However, at the national level, only a small proportion of the other enterovirus serotypes were further identified as CV-A6 and CV-A10.

WHAT ARE THE LIMITATIONS FOR PUBLIC HEALTH PRACTICE?

Further identification of other enterovirus serotypes is highly recommended for provincial CDCs, especially for severe HFMD. Studies contributing to a multivalent vaccine for HFMD should be prioritized.

Pathogenic characteristics of hand, foot and mouth disease in Shaanxi Province, China, 2010-2016

Hand, foot, and mouth disease (HFMD) is a common childhood illness caused by enteroviruses. We analyzed the pathogenic characteristics of HFMD in Shaanxi province, China, during 2010-2016. Clinical samples were collected from HFMD cases. Real-time PCR and RT-PCR were used to identify the enterovirus(EVs) serotypes. Viral RNA sequences were amplified using RT-PCR and compared by phylogenetic analysis. Descriptive epidemiological methods were used to analyze. A total of 16,832 HFMD positive cases were confirmed in the laboratory. EV-A71 and CV-A16 were the main pathogens in 2010. EV-A71 was the dominant pathogen in the periods of 2011 to 2012 and 2014, 2016. In 2013 and 2015, other EVs increased greatly, in which CV-A6 was the predominant pathogen. EV-A71 was more frequently detected in deaths and severe cases. Phylogenetic analysis revealed that EV-A71 belonged to the C4a evolution branch of C4 sub-genotype and CV-A16 belonged to the B1a or B1b evolution branch of B1 sub-genotype, whereas CV-A6 strains were assigned to D2 or D3 sub-genotype. The pathogen spectrum of HFMD has changed in 7 years, and the major serotypes EV-A71, CV- A16 and CV- A6 alternated or co-circulated. Long-term surveillance and research of EVs should be strengthened for the prevention and control of HFMD.