Coxsackievirus A16 (CVA16) is well known for causing hand-foot-and-mouth disease (HFMD) and outbreaks were frequently reported in Taiwan in the past twenty years. The epidemiology and genetic variations of CVA16 in Taiwan from 1998 to 2021 were analyzed in this study. CVA16 infections usually occurred in early summer and early winter, and showed increased incidence in 1998, 2000-2003, 2005, 2007-2008, and 2010 in Taiwan. Little or no CVA16 was detected from 2017 to 2021. CVA16 infection was prevalent in patients between 1 to 3 years old. A total of 69 isolates were sequenced. Phylogenetic analysis based on the VP1 region showed that CVA16 subgenotype B1 was dominantly isolated in Taiwan from 1998 to 2019, and B2 was identified only from isolates collected in 1999 and 2000. There was a high frequency of synonymous mutations in the amino acid sequences of the VP1 region among CVA16 isolates, with the exception of position 145 which showed positive selection. The recombination analysis of the whole genome of CVA16 isolates indicated that the 5'-untranslated region and the non-structural protein region of CVA16 subgenotype B1 were recombined with Coxsackievirus A4 (CVA4) and enterovirus A71 (EVA71) genotype A, respectively. The recombination pattern of subgenotype B2 was similar to B1, however, the 3D region was similar to EVA71 genotype B. Cross-neutralization among CVA16 showed that mouse antisera from various subgenotypes viruses can cross-neutralize different genotype with high neutralizing antibody titers. These results suggest that the dominant CVA16 genotype B1 can serve as a vaccine candidate for CVA16.
