Highly pathogenic avian influenza virus subtype H5N1 escaping neutralization: more than HA variation

Influenza A viruses are one of the major threats in modern health care. Novel viruses arise due to antigenic drift and antigenic shift, leading to escape from the immune system and resulting in a serious problem for disease control. In order to investigate the escape process and to enable predictions of escape, we serially passaged influenza A H5N1 virus in vitro 100 times under immune pressure. The generated escape viruses were characterized phenotypically and in detail by full-genome deep sequencing. Mutations already found in natural isolates were detected, evidencing the in vivo relevance of the in vitro-induced amino acid substitutions. Additionally, several novel alterations were triggered. Altogether, the results imply that our in vitro system is suitable to study influenza A virus evolution and that it might even be possible to predict antigenic changes of influenza A viruses circulating in vaccinated populations.     

Antigenic drift in influenza virus H3 hemagglutinin from 1968 to 1980: multiple evolutionary pathways and sequential amino acid changes at key antigenic sites

Surveys of the antigenic properties of a wide range of variants of the H3N2 (Hong Kong) influenza virus subtype have revealed complex patterns of variants cocirculating during each of the main epidemic eras of the subtype. We determined hemagglutinin (HA) gene sequences for 14 isolates chosen to give the wildest possible spread of variant types. The addition of these data to existing HA gene sequence information for other variants provides a comprehensive picture of HA gene evolution during antigenic drift among H3N2 subtype viruses. The data reveal the existence of multiple evolutionary pathways during at least one period of development of the subtype and strikingly demonstrate that amino acid changes are limited to a small number of locations on the HA molecule during antigenic drift. The occurrence of sequential amino acid changes at key positions within these variable regions suggests that the HA structure has remained constant during subtype evolution so that only limited possibilities remain for further antigenic drift among H3N2 viruses.     

Antigenic drift in H5N1 avian influenza virus in poultry is driven by mutations in major antigenic sites of the hemagglutinin molecule analogous to those for human influenza virus

H5N1 highly pathogenic avian influenza virus has been endemic in poultry in Egypt since 2008, notwithstanding the implementation of mass vaccination and culling of infected birds. Extensive circulation of the virus has resulted in a progressive genetic evolution and an antigenic drift. In poultry, the occurrence of antigenic drift in avian influenza viruses is less well documented and the mechanisms remain to be clarified. To test the hypothesis that H5N1 antigenic drift is driven by mechanisms similar to type A influenza viruses in humans, we generated reassortant viruses, by reverse genetics, that harbored molecular changes identified in genetically divergent viruses circulating in the vaccinated population. Parental and reassortant phenotype viruses were antigenically analyzed by hemagglutination inhibition (HI) test and microneutralization (MN) assay. The results of the study indicate that the antigenic drift of H5N1 in poultry is driven by multiple mutations primarily occurring in major antigenic sites at the receptor binding subdomain, similarly to what has been described for human influenza H1 and H3 subtype viruses.     

Glycosylation site alteration in the evolution of influenza A (H1N1) viruses

Influenza virus typically alters protein glycosylation in order to escape immune pressure from hosts and hence to facilitate survival in different host environments. In this study, the patterns and conservation of glycosylation sites on HA and NA of influenza A/H1N1 viruses isolated from various hosts at different time periods were systematically analyzed, by employing a new strategy combining genome-based glycosylation site prediction and 3D modeling of glycoprotein structures, for elucidation of the modes and laws of glycosylation site alteration in the evolution of influenza A/H1N1 viruses. The results showed that influenza H1N1 viruses underwent different alterations of protein glycosylation in different hosts. Two alternative modes of glycosylation site alteration were involved in the evolution of human influenza virus: One was an increase in glycosylation site numbers, which mainly occurred with high frequency in the early stages of evolution. The other was a change in the positional conversion of the glycosylation sites, which was the dominating mode with relatively low frequency in the later evolutionary stages. The mechanisms and possibly biological functions of glycosylation site alteration for the evolution of influenza A/H1N1 viruses were also discussed. Importantly, the significant role of positional alteration of glycosylation sites in the host adaptation of influenza virus was elucidated. Although the results still need to be supported by experimental data, the information here may provide some constructive suggestions for research into the glycosylation of influenza viruses as well as even the design of surveillance and the production of viral vaccines.     

Computational Identification of Antigenicity-Associated Sites in the Hemagglutinin Protein of A/H1N1 Seasonal Influenza Virus

The antigenic variability of influenza viruses has always made influenza vaccine development challenging. The punctuated nature of antigenic drift of influenza virus suggests that a relatively small number of genetic changes or combinations of genetic changes may drive changes in antigenic phenotype. The present study aimed to identify antigenicity-associated sites in the hemagglutinin protein of A/H1N1 seasonal influenza virus using computational approaches. Random Forest Regression (RFR) and Support Vector Regression based on Recursive Feature Elimination (SVR-RFE) were applied to H1N1 seasonal influenza viruses and used to analyze the associations between amino acid changes in the HA1 polypeptide and antigenic variation based on hemagglutination-inhibition (HI) assay data. Twenty-three and twenty antigenicity-associated sites were identified by RFR and SVR-RFE, respectively, by considering the joint effects of amino acid residues on antigenic drift. Our proposed approaches were further validated with the H3N2 dataset. The prediction models developed in this study can quantitatively predict antigenic differences with high prediction accuracy based only on HA1 sequences. Application of the study results can increase understanding of H1N1 seasonal influenza virus antigenic evolution and accelerate the selection of vaccine strains.     

Positive selection for gains of N-linked glycosylation sites in hemagglutinin during evolution of H3N2 human influenza A virus

The number of N-linked glycosylation sites in the globular head of hemagglutinin (HA) has increased during evolution of H3N2 human influenza A virus. Here natural selection operating on the gains of N-linked glycosylation sites was examined by using the single-site analysis and the single-substitution analysis. In the single-site analysis, positive selection was not inferred at the amino acid sites where the substitutions generating N-linked glycosylation sites were observed, but was detected at antigenic sites. In contrast, in the single-substitution analysis, positive selection was detected for the amino acid substitutions generating N-linked glycosylation sites. The single-site analysis and the single-substitution analysis appeared to be suitable for detecting recurrent and episodic natural selection, respectively. The gains of N-linked glycosylation sites were likely to be positively selected for the function of shielding antigenic sites from immune responses. At the antigenic sites, positive selection appeared to have operated not only on the radical substitution but also on the conservative substitution in terms of the charge of amino acids, suggesting that the antigenic drift is not a by-product of the evolution of receptor binding avidity in HA of human H3N2 virus.     

Subtype- and antigenic site-specific differences in biophysical influences on evolution of influenza virus hemagglutinin

ABSTRACT: BACKGROUND: Influenza virus undergoes rapid evolution by both antigenic shift and antigenic drift. Antibodies, particularly those binding near the receptor-binding site of hemagglutinin (HA) or the neuraminidase (NA) active site, are thought to be the primary defense against influenza infection, and mutations in antibody binding sites can reduce or eliminate antibody binding. The binding of antibodies to their cognate antigens is governed by such biophysical properties of the interacting surfaces as shape, non-polar and polar surface area, and charge. Methods: To understand forces shaping evolution of influenza virus, we have examined HA sequences of human influenza A and B viruses, assigning each amino acid values reflecting total accessible surface area, non-polar and polar surface area, and net charge due to the side chain. Changes in each of these values between neighboring sequences were calculated for each residue and mapped onto the crystal structures. Results: Areas of HA showing the highest frequency of changes agreed well with previously identified antigenic sites in H3 and H1 HAs, and allowed us to propose more detailed antigenic maps and novel antigenic sites for H1 and influenza B HA. Changes in biophysical properties differed between HAs of different subtypes, and between different antigenic sites of the same HA. For H1, statistically significant differences in several biophysical quantities compared to residues lying outside antigenic sites were seen for some antigenic sites but not others. Influenza B antigenic sites all show statistically significant differences in biophysical quantities for all antigenic sites, whereas no statistically significant differences in biophysical quantities were seen for any antigenic site is seen for H3. In many cases, residues previously shown to be under positive selection at the genetic level also undergo rapid change in biophysical properties. Conclusions: The biophysical consequences of amino acid changes introduced by antigenic drift vary from subtype to subtype, and between different antigenic sites. This suggests that the significance of antibody binding in selecting new variants may also be variable for different antigenic sites and influenza subtypes.     

Patterns of predicted T-cell epitopes associated with antigenic drift in influenza H3N2 hemagglutinin

Antigenic drift allowing escape from neutralizing antibodies is an important feature of transmission and survival of influenza viruses in host populations. Antigenic drift has been studied in particular detail for influenza A H3N2 and well defined antigenic clusters of this virus documented. We examine how host immunogenetics contributes to determination of the antibody spectrum, and hence the immune pressure bringing about antigenic drift. Using uTOPE™ bioinformatics analysis of predicted MHC binding, based on amino acid physical property principal components, we examined the binding affinity of all 9-mer and 15-mer peptides within the hemagglutinin 1 (HA1) of 447 H3N2 virus isolates to 35 MHC-I and 14 MHC-II alleles. We provide a comprehensive map of predicted MHC-I and MHC-II binding affinity for a broad array of HLA alleles for the H3N2 influenza HA1 protein. Each HLA allele exhibited a characteristic predicted binding pattern. Cluster analysis for each HLA allele shows that patterns based on predicted MHC binding mirror those described based on antibody binding. A single amino acid mutation or position displacement can result in a marked difference in MHC binding and hence potential T-helper function. We assessed the impact of individual amino acid changes in HA1 sequences between 10 virus isolates from 1968–2002, representative of antigenic clusters, to understand the changes in MHC binding over time. Gain and loss of predicted high affinity MHC-II binding sites with cluster transitions were documented. Predicted high affinity MHC-II binding sites were adjacent to antibody binding sites. We conclude that host MHC diversity may have a major determinant role in the antigenic drift of influenza A H3N2.     

The evolutionary dynamics of receptor binding avidity in influenza A: a mathematical model for a new antigenic drift hypothesis

The most salient feature of influenza evolution in humans is its antigenic drift. This process is characterized by structural changes in the virus''s B-cell epitopes and ultimately results in the ability of the virus to evade immune recognition and thereby reinfect previously infected hosts. Until recently, amino acid substitutions in epitope regions of the viral haemagglutinin were thought to be positively selected for their ability to reduce antibody binding and therefore were thought to be responsible for driving antigenic drift. However, a recent hypothesis put forward by Hensley and co-workers posits that cellular receptor binding avidity is the dominant phenotype under selection, with antigenic drift being a side effect of these binding avidity changes. Here, we present a mathematical formulation of this new antigenic drift model and use it to show how rates of antigenic drift depend on epidemiological parameters. We further use the model to evaluate how two different vaccination strategies can impact antigenic drift rates and ultimately disease incidence levels. Finally, we discuss the assumptions present in the model formulation, predictions of the model, and future work that needs to be done to determine the consistency of this hypothesis with known patterns of influenza''s genetic and antigenic evolution.     

The evolution of human influenza viruses.

The evolution of influenza viruses results in (i) recurrent annual epidemics of disease that are caused by progressive antigenic drift of influenza A and B viruses due to the mutability of the RNA genome and (ii) infrequent but severe pandemics caused by the emergence of novel influenza A subtypes to which the population has little immunity. The latter characteristic is a consequence of the wide antigenic diversity and peculiar host range of influenza A viruses and the ability of their segmented RNA genomes to undergo frequent genetic reassortment (recombination) during mixed infections. Contrasting features of the evolution of recently circulating influenza AH1N1, AH3N2 and B viruses include the rapid drift of AH3N2 viruses as a single lineage, the slow replacement of successive antigenic variants of AH1N1 viruses and the co-circulation over some 25 years of antigenically and genetically distinct lineages of influenza B viruses. Constant monitoring of changes in the circulating viruses is important for maintaining the efficacy of influenza vaccines in combating disease.     

Studies on the adaptation of influenza viruses to MDCK cells

The amino acid sequences and biological properties of the haemagglutinin of three variants of the influenza virus X-31 (H3N2) selected for their capacity to grow in MDCK cells are reported. In two variants, amino acid substitutions at HA1 residues 8 and 144 correlated with the loss of a site for glycosylation and specific changes in antigenicity, respectively. In all three variants substitution of an arginine residue for histidine at HA1 position 17 was correlated with increased pH optima of haemolysis. The importance of this substitution for cleavage of the haemagglutinin precursor required to produce infectious virus is discussed in relation to the three-dimensional structure of X-31 haemagglutinin.     

Antigenic and genetic evolution of swine influenza A (H3N2) viruses in Europe

In the early 1970s, a human influenza A/Port Chalmers/1/73 (H3N2)-like virus colonized the European swine population. Analyses of swine influenza A (H3N2) viruses isolated in The Netherlands and Belgium revealed that in the early 1990s, antigenic drift had occurred, away from A/Port Chalmers/1/73, the strain commonly used in influenza vaccines for pigs. Here we show that Italian swine influenza A (H3N2) viruses displayed antigenic and genetic changes similar to those observed in Northern European viruses in the same period. We used antigenic cartography methods for quantitative analyses of the antigenic evolution of European swine H3N2 viruses and observed a clustered virus evolution as seen for human viruses. Although the antigenic drift of swine and human H3N2 viruses has followed distinct evolutionary paths, potential cluster-differentiating amino acid substitutions in the influenza virus surface protein hemagglutinin (HA) were in part the same. The antigenic evolution of swine viruses occurred at a rate approximately six times slower than the rate in human viruses, even though the rates of genetic evolution of the HA at the nucleotide and amino acid level were similar for human and swine H3N2 viruses. Continuous monitoring of antigenic changes is recommended to give a first indication as to whether vaccine strains may need updating. Our data suggest that humoral immunity in the population plays a smaller role in the evolutionary selection processes of swine H3N2 viruses than in human H3N2 viruses.     

Prediction of biological functions on glycosylation site migrations in human influenza H1N1 viruses

Protein glycosylation alteration is typically employed by various viruses for escaping immune pressures from their hosts. Our previous work had shown that not only the increase of glycosylation sites (glycosites) numbers, but also glycosite migration might be involved in the evolution of human seasonal influenza H1N1 viruses. More importantly, glycosite migration was likely a more effectively alteration way for the host adaption of human influenza H1N1 viruses. In this study, we provided more bioinformatics and statistic evidences for further predicting the significant biological functions of glycosite migration in the host adaptation of human influenza H1N1 viruses, by employing homology modeling and in silico protein glycosylation of representative HA and NA proteins as well as amino acid variability analysis at antigenic sites of HA and NA. The results showed that glycosite migrations in human influenza viruses have at least five possible functions: to more effectively mask the antigenic sites, to more effectively protect the enzymatic cleavage sites of neuraminidase (NA), to stabilize the polymeric structures, to regulate the receptor binding and catalytic activities and to balance the binding activity of hemagglutinin (HA) with the release activity of NA. The information here can provide some constructive suggestions for the function research related to protein glycosylation of influenza viruses, although these predictions still need to be supported by experimental data.     

Evolutionary pathways of the pandemic influenza A (H1N1) 2009 in the UK

The emergence of the influenza (H1N1) 2009 virus provided a unique opportunity to study the evolution of a pandemic virus following its introduction into the human population. Virological and clinical surveillance in the UK were comprehensive during the first and second waves of the pandemic in 2009, with extensive laboratory confirmation of infection allowing a detailed sampling of representative circulating viruses. We sequenced the complete coding region of the haemagglutinin (HA) segment of 685 H1N1 pandemic viruses selected without bias during two waves of pandemic in the UK (April-December 2009). Phylogenetic analysis showed that although temporal accumulation of amino acid changes was observed in the HA sequences, the overall diversity was less than that typically seen for seasonal influenza A H1N1 or H3N2. There was co-circulation of multiple variants as characterised by signature amino acid changes in the HA. A specific substitution (S203T) became predominant both in UK and global isolates. No antigenic drift occurred during 2009 as viruses with greater than four-fold reduction in their haemagglutination inhibition (HI) titre ("low reactors") were detected in a low proportion (3%) and occurred sporadically. Although some limited antigenic divergence in viruses with four-fold reduction in HI titre might be related to the presence of 203T, additional studies are needed to test this hypothesis.     

Establishment and lineage replacement of H6 influenza viruses in domestic ducks in southern China

Domestic ducks in southern China act as an important reservoir for influenza viruses and have also facilitated the establishment of multiple H6 influenza virus lineages. To understand the continuing evolution of these established lineages, 297 H6 viruses isolated from domestic ducks during 2006 and 2007 were genetically and antigenically analyzed. Phylogenetic analyses showed that group II duck H6 viruses had replaced the previously predominant group I lineage and extended their geographic distribution from coastal to inland regions. Group II H6 virus showed that the genesis and development of multiple types of deletions in the neuraminidase (NA) stalk region could occur in the influenza viruses from domestic ducks. A gradual replacement of the N2 NA subtype with N6 was observed. Significant antigenic changes occurred within group II H6 viruses so that they became antigenically distinguishable from group I and gene pool viruses. Gene exchange between group II H6 viruses and the established H5N1, H9N2, or H6N1 virus lineages in poultry in the region was very limited. These findings suggest that domestic ducks can facilitate significant genetic and antigenic changes in viruses established in this host and highlight gaps in our knowledge of influenza virus ecology and even the evolutionary behavior of this virus family in its aquatic avian reservoirs.     

Analysis of Adaptation Mutants in the Hemagglutinin of the Influenza A(H1N1)pdm09 Virus

Hemagglutinin is the major surface glycoprotein of influenza viruses. It participates in the initial steps of viral infection through receptor binding and membrane fusion events. The influenza pandemic of 2009 provided a unique scenario to study virus evolution. We performed molecular dynamics simulations with four hemagglutinin variants that appeared throughout the 2009 influenza A (H1N1) pandemic. We found that variant 1 (S143G, S185T) likely arose to avoid immune recognition. Variant 2 (A134T), and variant 3 (D222E, P297S) had an increased binding affinity for the receptor. Finally, variant 4 (E374K) altered hemagglutinin stability in the vicinity of the fusion peptide. Variants 1 and 4 have become increasingly predominant, while variants 2 and 3 declined as the pandemic progressed. Our results show some of the different strategies that the influenza virus uses to adapt to the human host and provide an example of how selective pressure drives antigenic drift in viral proteins.     

Glycosylation and an amino acid insertion in the head of hemagglutinin independently affect the antigenic properties of H5N1 avian influenza viruses

Antigenic drift forces us to frequently update influenza vaccines; however, the genetic basis for antigenic variation remains largely unknown. In this study, we used clade 7.2 H5 viruses as models to explore the molecular determinants of influenza virus antigenic variation. We generated eight monoclonal antibodies(MAbs) targeted to the hemagglutinin(HA) protein of the index virus A/chicken/Shanxi/2/2006 and found that two representative antigenically drifted clade 7.2 viruses did not react with six of the eight MAbs. The E131 N mutation and insertion of leucine at position 134 in the HA protein of the antigenically drifted strains eliminated the reactivity of the virus with the MAbs. We also found that the amino acid N131 in the H5 HA protein is glycosylated. Our results provide experimental evidence that glycosylation and an amino acid insertion or deletion in HA influence antigenic variation.     

Glycosylation focuses sequence variation in the influenza A virus H1 hemagglutinin globular domain

Antigenic drift in the influenza A virus hemagglutinin (HA) is responsible for seasonal reformulation of influenza vaccines. Here, we address an important and largely overlooked issue in antigenic drift: how does the number and location of glycosylation sites affect HA evolution in man? We analyzed the glycosylation status of all full-length H1 subtype HA sequences available in the NCBI influenza database. We devised the “flow index” (FI), a simple algorithm that calculates the tendency for viruses to gain or lose consensus glycosylation sites. The FI predicts the predominance of glycosylation states among existing strains. Our analyses show that while the number of glycosylation sites in the HA globular domain does not influence the overall magnitude of variation in defined antigenic regions, variation focuses on those regions unshielded by glycosylation. This supports the conclusion that glycosylation generally shields HA from antibody-mediated neutralization, and implies that fitness costs in accommodating oligosaccharides limit virus escape via HA hyperglycosylation.     

Glycosylation at Asn91 of H1N1 haemagglutinin affects binding to glycan receptors

The glycoprotein HA (haemagglutinin) on the surface of influenza A virus plays a central role in recognition and binding to specific host cell-surface glycan receptors and in fusion of viral membrane to the host nuclear membrane during viral replication. Given the abundance of HA on the viral surface, this protein is also the primary target for host innate and adaptive immune responses. Although addition of glycosylation sites on HA are a part of viral evolution to evade the host immune responses, there are specific glycosylation sites that are conserved during most of the evolution of the virus. In the present study, it was demonstrated that one such conserved glycosylation site at Asn(91) in H1N1 HA critically governs the glycan receptor-binding specificity and hence would potentially impinge on the host adaptation of the virus.     

Identification of Low- and High-Impact Hemagglutinin Amino Acid Substitutions That Drive Antigenic Drift of Influenza A(H1N1) Viruses

Determining phenotype from genetic data is a fundamental challenge. Identification of emerging antigenic variants among circulating influenza viruses is critical to the vaccine virus selection process, with vaccine effectiveness maximized when constituents are antigenically similar to circulating viruses. Hemagglutination inhibition (HI) assay data are commonly used to assess influenza antigenicity. Here, sequence and 3-D structural information of hemagglutinin (HA) glycoproteins were analyzed together with corresponding HI assay data for former seasonal influenza A(H1N1) virus isolates (1997–2009) and reference viruses. The models developed identify and quantify the impact of eighteen amino acid substitutions on the antigenicity of HA, two of which were responsible for major transitions in antigenic phenotype. We used reverse genetics to demonstrate the causal effect on antigenicity for a subset of these substitutions. Information on the impact of substitutions allowed us to predict antigenic phenotypes of emerging viruses directly from HA gene sequence data and accuracy was doubled by including all substitutions causing antigenic changes over a model incorporating only the substitutions with the largest impact. The ability to quantify the phenotypic impact of specific amino acid substitutions should help refine emerging techniques that predict the evolution of virus populations from one year to the next, leading to stronger theoretical foundations for selection of candidate vaccine viruses. These techniques have great potential to be extended to other antigenically variable pathogens.