Abstract
Human parainfluenza virus 2 (HPIV-2) and human parainfluenza virus 4 (HPIV-4) are significant but underappreciated respiratory pathogens, particularly among high-risk populations including children, the elderly, and immunocompromised individuals. In this study, we sequenced 101 HPIV-2 and HPIV-4 genomes from respiratory samples collected in western Washington State and performed comprehensive evolutionary analyses using both new and publicly available sequences. Phylogenetic and phylodynamic analyses revealed that both HPIV-2 and HPIV-4 evolve at significantly faster rates compared to the mumps virus, a reference human orthorubulavirus. Notably, while HPIV-2 demonstrated the highest evolutionary rates in the surface glycoprotein HN, consistent with humoral immune-driven selection, the innate immune antagonist V/P gene evolved fastest in HPIV-4. We identified a hypervariable region within the HPIV-4V/P protein (residues 35 to 75), which structural modeling placed in a loop overlapping a known interferon antagonism domain in other paramyxovirus V proteins, though HPIV-4 is functionally incompetent in this activity. Expanded phylogenetic analysis across the Paramyxoviridae family uncovered a striking evolutionary discordance: while the HN glycoprotein and L polymerase of HPIV-4 and its 2 closest bat-derived viruses clustered within the Orthorubulavirus genus, their nucleoprotein (N), phosphoprotein (P), matrix (M), and fusion (F) proteins formed a distinct lineage outside the Rubulavirinae subfamily. Together, these findings highlight the distinct evolutionary trajectories of HPIV-2 and HPIV-4, raise hypotheses around complex Paramyxoviridae zoonotic events including recombination-like patterns, and demonstrate limitations of current L protein-based taxonomic classification schemes.