Abstract
Avian influenza virus (AIV) H7N9 that emerged in 2013 in eastern China is a novel zoonotic agent mainly circulating in poultry without clinical signs but causing severe disease with high fatality in humans in more than 5 waves. Since the emergence of highly pathogenic (HP) H7N9 variants in 2016, it has induced heavy losses in the poultry industry leading to the implementation of an intensive nationwide vaccination program at the end of wave 5 (September 2017). To characterize the ongoing evolution of H7N9 AIV, we conducted analyses of H7N9 glycoprotein genes obtained from 2013 to 2019. Bayesian analyses revealed a decreasing population size of HP H7N9 variants post wave 5. Phylogenetic topologies revealed that two novel small subclades were formed and carried several fixed amino acid mutations that were along HA and NA phylogenetic trees since wave 5. Some of the mutations were located at antigenic sites or receptor binding sites. The antigenic analysis may reveal a significant antigenic drift evaluated by hemagglutinin inhibition (HI) assay and the antigenicity of H7N9 AIV might evolute in large leaps in wave 7. Molecular simulations found that the mutations (V135T, S145P, and L226Q) around the HA receptor pocket increased the affinity to α2,3-linked sialic acid (SIA) while decreased to α2,6-linked SIA. Altered affinity may suggest that HP H7N9 variations aggravate the pathogenicity to poultry but lessen the threat to public health. Selection analyses showed that the HP H7N9 AIV experienced an increasing selection pressure since wave 5, and the national implementation of vaccination might intensify the role of natural selection during the evolution waves 6 and 7. In summary, our data provide important insights about the genetic and antigenic diversity of circulating HP H7N9 viruses from 2017 to 2019. Enhanced surveillance is urgently warranted to understand the current situation of HP H7N9 AIV.
