Abstract
Clade 2.2 H5N1 influenza viruses have caused an unusually high number of human infections, providing a unique opportunity to investigate early molecular steps associated with host adaptation. Although most work has focused on hemagglutinin (HA), the contribution of neuraminidase (NA) to these early adaptive events has remained unclear. By analyzing publicly available sequences from clade 2.2-infected patients, we identified 20 NA mutations and compared their phenotypes to 20 mutations acquired during diversification in primary human airway cells under drug-free conditions. Most patient-derived NA mutations resulted in modest reductions in sialidase activity, keeping activity within a functional range that supported improved replication in α2,6 sialylglycan (α2,6 Sia)-dominant environments, whereas excessive reduction impaired fitness. Notably, the phenotypes of culture-selected and patient-derived mutations were highly concordant, suggesting that these NA changes arose through natural selection rather than antiviral pressure. Re-analysis of patient sequences further revealed that many adaptive NA mutations co-occur with HA mutations that confer only weak, partial α2,6 Sia binding. Using reverse genetics, we found that such naturally occurring HA/NA mutation pairs acted cooperatively in a receptor-context-dependent manner to support α2,6-associated replication relative to HA-only mutants, placing these variants within a constrained "early-adaptation space" characterized by limited α2,6 engagement and moderately reduced NA activity. Together, these findings indicate that early human adaptation of clade 2.2 H5N1 involves not only HA and PB2, but also incremental, cooperative tuning of NA function. Monitoring coordinated HA-NA evolution may therefore improve risk assessment frameworks for zoonotic influenza viruses poised at early stages of human host adaptation.