Phylogenetic Analyses of RdRp Region and VP1 Gene in Human Norovirus Genotype GII.17[P17] Variants

In this study, we investigated the long-term evolutionary dynamics of human norovirus GII.17[P17] using the RNA-dependent RNA polymerase (RdRp) region and the VP1 capsid gene, integrating phylogenetics, time-scaled inference, phylodynamics, and structure-based analyses. Maximum-likelihood phylogenies of both genomic regions consistently resolved four major clades (Clades 1-4). VP1 patristic-distance distributions indicated higher within-clade diversity in the phylogenetically basal Clades 1 and 3, whereas Clades 2 and 4 showed lower diversity, consistent with recent demographic expansion. Similarity-plot analysis identified pronounced variability in the VP1 P2 domain, while the S and P1 domains remained comparatively conserved, supporting P2 as the primary hotspot of diversification. Bayesian time-scaled analyses estimated the most recent common ancestor around 1993 (VP1) and 2000 (RdRp) and revealed two major lineages (Clade 1/2 and Clade 3/4), with the split between Clades 3 and 4 occurring around 2016-2017. Bayesian skyline plots showed a marked increase in effective population size after 2013, and substitution-rate estimates indicated faster evolution in VP1 than in RdRp, with higher VP1 rates in the Clade 3/4 lineage than in Clade 1/2. Capsid dimer modeling further mapped high-confidence conformational B-cell epitopes and positively selected residues predominantly to the distal surface of P2, with broadly conserved spatial patterns across clades. Compared with the Clade 1 reference (Kawasaki323), Clade 2 accumulated numerous P2 substitutions, whereas Clades 3 and 4 retained fewer changes and remained closer to Clade 1 at the amino-acid level. Together, these results suggest lineage turnover within GII.17[P17] driven by constrained diversification at the P2 surface, potentially contributing to the recent predominance of the Clade 3/4 lineage.