AbstractHuman rhinovirus species C (HRV-C) was recently discovered and this virus has been associated with various acute respiratory illnesses (ARI). However, the molecular evolution of the major antigens of this virus, including VP1, VP2 and VP3, is unknown. Thus, we performed complete VP1, VP2, and VP3 gene analyses of 139 clinical HRV-C strains using RT-PCR with newly designed primer sets and next-generation sequencing. We assessed the time-scale evolution and evolutionary rate of these genes using the Bayesian Markov chain Monte Carlo method. In addition, we calculated the pairwise distance and confirmed the positive/negative selection sites in these genes. The phylogenetic trees showed that the HRV-C strains analyzed using these genes could be dated back approximately 400 to 900 years and these strains exhibited high evolutionary rates (1.35 to 3.74 × 10−3 substitutions/site/year). Many genotypes (>40) were confirmed in the phylogenetic trees. Furthermore, no positively selected site was found in the VP1, VP2 and VP3 protein. Molecular modeling analysis combined with variation analysis suggested that the exterior surfaces of the VP1, VP2 and VP3 proteins are rich in loops and are highly variable. These results suggested that HRV-C may have an old history and unique antigenicity as an agent of various ARI.