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

Introduction In recent years, coxsackievirus (CV) A10 has been associated with increasing sporadic hand, foot, and mouth disease (HFMD) cases and outbreaks globally. In addition to mild symptoms such as pharyngitis and herpangina, CVA10-related complications or even fatality can occur. Currently, systematic phylogenetic studies of CVA10 are limited. Methods In this study, we first explored the epidemiological and genetic characteristics of CVA10 in Nanchang, an inland southeastern city of China, based on the HFMD surveillance network from 2015-2023. Results Among 3429 enterovirus-positive cases, 110 (3.04%) were associated with CVA10, with a male-to-female ratio of 1.62. The median age of the CVA10 patients was 2.3 years (interquartile range, IQR 1.0-4.0), with 94.55% (104/110) of the patients aged less than 5 years. Phylogenetic analyses using the full-length VP1, 5’UTR, P1, P2, P3 sequences and near full-length genomes indicated that CVA10 strains ( n = 57) isolated in Nanchang belonged to genogroup C; two strains identified in 2017 belonged to C1 subtypes clustered with strains from Vietnam, Madagascar, France and Spain; and the others belonged to C2 subtypes interdigitating with CVA10 isolates from mainland China, the United States and Australia. Through extensive analysis, we identified a rare F168Y mutation in epitope 4 of VP1 in a Madagascar strain of genogroup F and a Chinese strain of genogroup C. Based on Bayesian evolutionary analyses, the average nucleotide substitution rate for the VP1 gene of CV10 strains was 3.07×10 –3 substitutions/site/year. The most recent common ancestor (tMRCA) of genogroup C was dated 1990.84, and the tMRCA of CVA10 strains from Nanchang was dated approximately 2003.16, similar to strains circulating in other regions of China, suggesting that the viruses were likely introduced and cryptically circulated in China before the establishment of the HFMD surveillance network. Recombination analysis indicated intertypic recombination of the Nanchang strain with the genogroup G strain in the 3D region. Discussion Given the shifting dominance of viral genotypes and frequent recombination events, the existing surveillance system needs to be regulated to enhance genomic surveillance efforts on a more diverse spectrum of genotypes in the future.