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.
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.