Recent phylodynamic studies have focused on using tree topology patterns to elucidate interactions among the epidemiological, evolutionary, and demographic characteristics of infectious agents. However, because studies of viral phylodynamics tend to focus on epidemic outbreaks, tree topology signatures of tissue-tropism pathogens might not be clearly identified. Therefore, this study used a novel Bayesian evolutionary approach to analyze the A24 variant of coxsackievirus (CV-A24v), an ocular-tropism agent of acute hemorrhagic conjunctivitis. Analyses of the 915-nucleotide VP1 and 690-nt 3Dpol regions of 21 strains isolated in Taiwan and worldwide during 1985–2010 revealed a clear chronological trend in both the VP1 and 3Dpol phylogenetic trees: the emergence of a single dominant cluster in each outbreak. The VP1 sequences included three genotypes: GI (prototype), GIII (isolated 1985–1999), and GIV (isolated after 2000); no VP1 sequences from GII strains have been deposited in GenBank. Another five genotypes identified in the 3Dpol region had support values >0.9. Geographic and demographic transitions among CV-A24v clusters were clearly identified by Bayes algorithm. The transmission route was mapped from India to China and then to Taiwan, and each prevalent viral population declined before new clusters emerged. Notably, the VP1 and 3Dpol genes had high nucleotide sequence similarities (94.1% and 95.2%, respectively). The lack of co-circulating lineages and narrow tissue tropism affected the CV-A24v gene pool.

Coxsackievirus (CV)-B5 is a common human enterovirus reported worldwide; swine vesicular disease virus (SVDV) is a porcine variant of CV-B5. To clarify the transmission dynamics and molecular basis of host switching between CV-B5 and SVDV, we analysed and compared the VP1 and partial 3Dpol gene regions of these two viruses. Spatiotemporal dynamics of viral transmission were estimated using a Bayesian statistical inference framework. The detected selection events were used to analyse the key molecules associated with host switching. Analyses of VP1 sequences revealed six CV-B5 genotypes (A1–A4 and B1–B2) and three SVDV genotypes (I–III). Analyses of partial 3Dpol revealed five clusters (A–E). The genotypes evolved sequentially over different periods, albeit with some overlap. The major hub of CV-B5 transmission was in China whereas the major hubs of SVDV transmission were in Italy. Network analysis based on deduced amino acid sequences showed a diverse extension of the VP1 structural protein, whereas most sequences were clustered into two haplotypes in the partial 3Dpol region. Residue 178 of VP1 showed four epistatic interactions with residues known to play essential roles in viral host tropism, cell entry, and viral decoating.