Abstract:Background: Hand, foot, and mouth disease (HFMD) is a global public health concern, notably within the Asia-Pacific region. Recently, the primary pathogen causing HFMD outbreaks across numerous countries, including China, is coxsackievirus (CV) A6, one of the most prevalent enteroviruses in the world. It is a new variant that has undergone genetic recombination and evolution, which might not only induce modifications in the clinical manifestations of HFMD but also heighten its pathogenicity because of nucleotide mutation accumulation. Objective: The study assessed the epidemiological characteristics of HFMD in China and characterized the molecular epidemiology of the major pathogen (CV-A6) causing HFMD. We attempted to establish the association between disease progression and viral genetic evolution through a molecular epidemiological study. Methods: Surveillance data from the Chinese Center for Disease Control and Prevention from 2021 to 2023 were used to analyze the epidemiological seasons and peaks of HFMD in Henan, China, and capture the results of HFMD pathogen typing. We analyzed the evolutionary characteristics of all full-length CV-A6 sequences in the NCBI database and the isolated sequences in Henan. To characterize the molecular evolution of CV-A6, time-scaled tree and historical population dynamics regarding CV-A6 sequences were estimated. Additionally, we analyzed the isolated strains for mutated or missing amino acid sites compared to the prototype CV-A6 strain. Results: The 2021-2023 epidemic seasons for HFMD in Henan usually lasted from June to August, with peaks around June and July. The monthly case reporting rate during the peak period ranged from 20.7% (4854/23,440) to 35% (12,135/34,706) of the total annual number of cases. Analysis of the pathogen composition of 2850 laboratory-confirmed cases identified 8 enterovirus serotypes, among which CV-A6 accounted for the highest proportion (652/2850, 22.88%). CV-A6 emerged as the major pathogen for HFMD in 2022 (203/732, 27.73%) and 2023 (262/708, 37.01%). We analyzed all CV-A6 full-length sequences in the NCBI database and the evolutionary features of viruses isolated in Henan. In China, the D3 subtype gradually appeared from 2011, and by 2019, all CV-A6 virus strains belonged to the D3 subtype. The VP1 sequences analyzed in Henan showed that its subtypes were consistent with the national subtypes. Furthermore, we analyzed the molecular evolutionary features of CV-A6 using Bayesian phylogeny and found that the most recent common ancestor of CV-A6 D3 dates back to 2006 in China, earlier than the 2011 HFMD outbreak. Moreover, the strains isolated in 2023 had mutations at several amino acid sites compared to the original strain. Conclusions: The CV-A6 virus may have been introduced and circulating covertly within China prior to the large-scale HFMD outbreak. Our laboratory testing data confirmed the fluctuation and periodic patterns of CV-A6 prevalence. Our study provides valuable insights into understanding the evolutionary dynamics of CV-A6.

Genetic characteristics of Coxsackievirus A6 from children with hand, foot and mouth disease in Beijing, China, 2017–2019

Epidemic Characteristics of Pathogen Spectrum and Cerebrospinal Fluid Analysis of Severe Hand , Foot , and Mouth Disease in Hangzhou ,2016

Molecular Epidemiology of Enterovirus from Children with Herpangina or Hand, Foot, and Mouth Disease in Hangzhou, 2016

Detection of human enterovirus 71 and coxsackievirus A16 in children with hand, foot and mouth disease in China.

Genetic Characteristics of the P1 Coding Region of Coxsackievirus A16 Associated with Hand, Foot, and Mouth Disease in China.

Genetic characteristics and phylogenetic analysis of coxsackievirus A6 isolated in Linyi, China, 2022-2023

Evolutionary Diversity of Coxsackievirus A6 Causing Severe Hand, Foot, and Mouth Disease - China, 2012-2023

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

The Epidemiological Study of Coxsackievirus A6 revealing Hand, Foot and Mouth Disease Epidemic patterns in Guangdong, China.

Circulation of Coxsackievirus A6 in hand-foot-mouth disease in Guangzhou, 2010-2012

Hand, Foot and Mouth Disease Caused by Coxsackievirus A6, Beijing, 2013.

Prevalence of Coxsackievirus A6 and Enterovirus 71 in Hand, Foot and Mouth Disease in Nanjing, China in 2013

An emerging and expanding clade accounts for the persistent outbreak of Coxsackievirus A6-associated hand, foot, and mouth disease in China since 2013.

Evolution, recombination and geographic spreading of global Coxsackievirus A6

VP1 Sequence Determination and Phylogenetic Analysis of Coxsackievirus A16 Strains

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

Molecular Epidemiology of Coxsackievirus A16: Intratype and Prevalent Intertype Recombination Identified

Coxsackievirus A6 Recombinant Subclades D3/A and D3/H Were Predominant in Hand-Foot-And-Mouth Disease Outbreaks in the Paediatric Population, France, 2010–2018

Co-Circulation And Genomic Recombination Of Coxsackievirus A16 And Enterovirus 71 During A Large Outbreak Of Hand, Foot, And Mouth Disease In Central China

Emergence, Circulation, and Spatiotemporal Phylogenetic Analysis of Coxsackievirus A6- and Coxsackievirus A10-Associated Hand, Foot, and Mouth Disease Infections from 2008 to 2012 in Shenzhen, China

Molecular Evolutionary Dynamics of Coxsackievirus A6 Causing Hand, Foot, and Mouth Disease From 2021 to 2023 in China: Genomic Epidemiology Study