Porcine circovirus type 2 (PCV2) is an important pathogen that causes significant economic losses in the swine industry worldwide. Five major PCV2 genotypes have been identified, including PCV2a, PCV2b, PCV2c, PCV2d, and PCV2e. To investigate the prevalence and phylodynamics of the different PCV2 genotypes in Taiwan, 214 PCV2 ORF2 sequences from Taiwan and other countries were analyzed. Genotypic differences were observed among PCV2a, 2b, and 2d at amino acid position 89 in ORF2, with isoleucine (I), arginine (R), and leucine (L), respectively. Similar to other countries, a genotypic shift was also observed in Taiwan, where the predominant genotype shifted from PCV2b to 2d after 2010. The estimated nucleotide substitution rate of Taiwanese strains in the ORF2 region was 8.467 × 10⁻⁴ substitutions per site per year. This rapid evolution rate of PCV2 may lead to the genotypic shift observed in Taiwan. The times to the most recent common ancestor (TMRCA) for PCV2a, -2b, and -2d-2 was dated to 1970, 1992 and 2004, respectively. Thus, the PCV2a, -2b, and -2d genotypes were already present in Taiwan before the introduction of the PCV2 vaccine.

Canine parvovirus type 2c (CPV-2c) emerged in 2000 and is known for causing a more severe disease than other CPV-2 variants in puppies. In 2015, the emerging CPV-2c variant was isolated in Taiwan and it subsequently became the predominant variant. To trace the evolution of Taiwanese CPV-2c, we compared complete VP2 genes of CPV-2c from Taiwan and sequences obtained from GenBank. The evolutionary rate of CPV-2c was estimated to be 4.586 × 10⁻⁴ substitutions per site per year (95% highest posterior density (HPD) was 3.284–6.076 × 10⁻⁴). The time to the most recent common ancestor (TMRCA) dated to 1990 (95% HPD: 1984–1996) and 2011 (95% HPD: 2010–2013) for the CPV-2c variant and Taiwanese isolates, respectively. The CPV-2c variant isolated from Taiwan was clustered with CPV-2c from China. This phylogenetic clade began to branch off in approximately 2010 (95% HPD was 3.823–6.497). Notably, two unique mutations of Taiwanese CPV-2c were found, Q383R and P410L. In summary, this is the first report on the genome evolution of CPV-2c in Taiwan, revealing that this CPV-2c variant shares a common evolutionary origin with strains from China. The demographic history inferred by the Bayesian skyline plot showed that the effective population of CPV-2c increased until 2006 and then slowly declined until 2011.

Taiwan has been recognized by the World Organization for Animal Health as rabies-free since 1961. Surprisingly, rabies virus (RABV) was identified in a dead Formosan ferret badger in July 2013. Later, more infected ferret badgers were reported from different geographic regions of Taiwan. In order to know its evolutionary history and spatial temporal dynamics of this virus, phylogeny was reconstructed by maximum likelihood and Bayesian methods based on the full-length of glycoprotein (G), matrix protein (M), and nucleoprotein (N) genes. The evolutionary rates and phylogeographic were determined using Beast and SPREAD software. Phylogenetic trees showed a monophyletic group containing all of RABV isolates from Taiwan and it further separated into three sub-groups. The estimated nucleotide substitution rates of G, M, and N genes were between 2.49 × 10⁻⁴–4.75 × 10⁻⁴ substitutions/site/year, and the mean ratio of dN/dS was significantly low. The time of the most recent common ancestor was estimated around 75, 89, and 170 years, respectively. Phylogeographic analysis suggested the origin of the epidemic could be in Eastern Taiwan, then the Formosan ferret badger moved across the Central Range of Taiwan to western regions and separated into two branches. In this study, we illustrated the evolution history and phylogeographic of RABV in Formosan ferret badgers.

Studies 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. Four 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. Phylodynamic analysis can reveal demographic trends and herd immunity in endemic pathogens.

Human bocavirus (HBoV) is a causative agent of respiratory and gastrointestinal diseases worldwide. Four HBoV species (HBoV1-4) have been identified so far. Although a previous report has documented the HBoV association with acute gastroenteritis (AGE) in Taiwan, their epidemiology, genetic diversity, and phylogenetic relationships remain unclear. In this study, we focused on an investigation of these unsolved issues, which will help to reveal molecular epidemiology and phylogeny of the circulating HBoV2 in Taiwan. A total of 176 stool samples were collected from children with AGE for this study. PCR amplification and sequencing on the VP1 gene region were used to identify species. Phylogenetic analysis was conducted by maximum-likelihood and neighbor-joining methods. Selection pressure was also estimated to obtain HBoV evolutionary information. Our results showed the prevalence of HBoV in AGE children was 8.5%, of which HBoV1 was the predominant species (6.3%), followed by HBoV2 (2.3%). Phylogenetic analysis showed those Taiwanese HBoV2 strains have significant genetic variability and can be divided into two clusters. One belongs to HBoV2 genotype A and the other forms an independent unclassified cluster. The nucleotide distance between that independent cluster and the known HBoV2 genotypes was more than 5%, suggesting a new HBoV2 genotype. No positive selection site was found and the virus was under purifying selection. This is the first report to reveal HBoV2 genetic diversity and phylogenetic relationships among AGE children in Taiwan. We find that HBoV2 may have been introduced into the country by multiple origins, and a potential new HBoV2 genotype is proposed.