Abstract
Hand, foot and mouth disease (HFMD) is a public health threat worldwide, particularly in the Asia-Pacific region. Enterovirus A71 (EV-A71), coxsackievirus A16 (CVA16), and CVA6 are the major pathogens causing HFMD outbreaks in several countries, including Thailand. We retrieved 385 VP1 nucleotide sequences, comprising 228 EV-A71, 33 CVA16, and 124 CVA6, deposited in the databases between 2000 and 2022 for molecular evolutionary characterization using Bayesian phylogeny. All EV-A71 identified belonged to genotype B, subgenotypes B4, and B5, and to genotype C, subgenotypes C1, C2, C4a, C4b, and C5. The analyzes demonstrated these viruses' co-circulation and subgenotypic changes throughout the past two decades. The CVA16 was grouped in genotype B1, predominantly subgenotype B1a, and the CVA6 was grouped in subgenotype D3, clades 1-4. The tMRCA of EV-A71 genotypes B and C, CVA16 B1, and CVA6 D3 dated 1993.79, 1982.62, 1995.86, and 2007.31, respectively, suggesting that the viruses were likely introduced and cryptically circulated in Thailand before the HFMD cases were recognized. We demonstrated these viruses' fluctuation and cyclical pattern throughout the two decades of observation. This study provided insight into evolutionary dynamics concerning molecular epidemiology and supported the selection of current genotype-matched vaccines, vaccine development, and implementation.