Molecular Evolution of the <i>Fusion</i> (<i>F</i>) Genes in Human Metapneumovirus Genotype B

Human metapneumovirus genotype B (HMPV-B) is an important respiratory pathogen, requiring detailed elucidation of the evolutionary and antigenic features of its fusion (F) gene. Using 500 sequences collected between 1982 and 2024, we investigated the molecular evolution, phylodynamics, and structural epitope landscape of the HMPV-B F gene. Time-scaled phylogeny dated the divergence of sublineages B1 and B2 to around 1938, and Bayesian Skyline Plot analysis showed that these sublineages exhibited distinct demographic trajectories over time. The F gene evolved at a rate of 1.01 × 10⁻³ substitutions/site/year; however, amino acid variation remained limited, consistent with pervasive purifying selection, with 39% of codons under strong negative selection and little consensus evidence for positive selection. Conformational B-cell epitope prediction demonstrated a high degree of conservation across neutralizing antibody binding regions (sites Ø and I–V), and amino acid substitutions occurring within these sites were not predicted to substantially alter epitope architecture. Together, these findings indicate that the HMPV-B F gene evolves under strong evolutionary constraint while maintaining stable antigenic features, supporting the potential for antibody-based strategies that target neutralizing antibody binding regions of the F protein.