BackgroundStudies regarding coxsackievirus (CV) tend to focus on epidemic outbreaks, an imbalanced topology is considered to be an indication of acute infection with partial cross-immunity. In enteroviruses, a clear understanding of the characteristics of tree topology, transmission, and its demographic dynamics in viral succession and circulation are essential for identifying prevalence trends in endemic pathogens such as coxsackievirus B2 (CV-B2). This study applied a novel Bayesian evolutionary approach to elucidate the phylodynamic characteristics of CV-B2. A dataset containing 51 VP1 sequences and a dataset containing 34 partial 3Dpol sequencing were analyzed, where each dataset included Taiwan sequences isolated during 1988–2013.ResultsFour and five genotypes were determined based on the 846-nucleotide VP1 and 441-nucleotide 3Dpol (6641–7087) regions, respectively, with spatiotemporally structured topologies in both trees. Some strains with tree discordance indicated the occurrence of recombination in the region between the VP1 and 3Dpol genes. The similarities of VP1 and 3Dpol gene were 80.0 %–96.8 % and 74.7 %–91.9 %, respectively. Analyses of population dynamics using VP1 dataset indicated that the endemic CV-B2 has a small effective population size. The balance indices, high similarity, and low evolutionary rate in the VP1 region indicated mild herd immunity selection in the major capsid region.ConclusionsPhylodynamic analysis can reveal demographic trends and herd immunity in endemic pathogens.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-015-0738-2) contains supplementary material, which is available to authorized users.