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.

Molecular evidence of persistent epidemic and evolution of subgenotype B1 coxsackievirus A16-associated hand, foot, and mouth disease in China.

Analysis of the Complete Genomes of Enterovirus 71 Subtypes in China.

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

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

VP1 Sequence Determination and Phylogenetic Analysis of Coxsackievirus A16 Strains

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

Genotyping and phylogeographic dynamics of coxsackievirus A16

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

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

Synergetic association between coxsackievirus A16 genotype evolution and recombinant form shifts

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

Genetic Divergence of Enterovirus D68 in China and the United States

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

Epidemiology of Hand, Foot, and Mouth Disease and Genetic Characterization of Coxsackievirus A16 in Shenyang, Liaoning Province, China, 2013–2023

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

Molecular Phylogeny of Coxsackievirus A16

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

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

Evolution, recombination and geographic spreading of global Coxsackievirus A6

Genomic surveillance reveals low-level circulation of two subtypes of genogroup C coxsackievirus A10 in Nanchang, Jiangxi Province, China, 2015-2023

Emerging Enteroviruses Causing Hand, Foot and Mouth Disease, China, 2010-2016.