Positive selection operates continuously on hemagglutinin during evolution of H3N2 human influenza A virus

It has been proposed that antigenic evolution of hemagglutinin 1 (HA1) for H3N2 human influenza A virus was punctuated. In the population genetic analysis, however, it was controversial whether positive selection operated on HA1 in a punctuated manner for the branches of the phylogenetic tree where transitions to new antigenic clusters occurred (C branches), or continuously. In the molecular evolutionary analysis, positive selection was detected for the trunk (T) branches but the relationship between antigenic evolution and positive selection was unclear. Here molecular evolutionary analysis was conducted to examine natural selection operating on HA1 of H3N2 human influenza A virus by dividing HA1 into epitopes A-E and other sites, as well as dividing the phylogenetic tree into the C branches overlapping with the T branches (C-T branches), those not overlapping with the T branches (C-NT branches), the T branches not overlapping with the C branches (NC-T branches), and other branches (NC-NT branches). Positive selection was detected for C, T, and NC-T branches, whereas evolution for the NC-NT branches appeared to be mainly neutral. Positive selection appeared to have operated throughout the trunk, which covered the entire time period of the phylogenetic tree, suggesting that positive selection operated continuously on HA1 during evolution of H3N2 human influenza A virus.     

Genetic Requirement for Hemagglutinin Glycosylation and Its Implications for Influenza A H1N1 Virus Evolution

Influenza A virus has evolved and thrived in human populations. Since the 1918 influenza A pandemic, human H1N1 viruses had acquired additional N-linked glycosylation (NLG) sites within the globular head region of hemagglutinin (HA) until the NLG-free HA head pattern of the 1918 H1N1 virus was renewed with the swine-derived 2009 pandemic H1N1 virus. Moreover, the HA of the 2009 H1N1 virus appeared to be antigenically related to that of the 1918 H1N1 virus. Hence, it is possible that descendants of the 2009 H1N1 virus might recapitulate the acquisition of HA head glycosylation sites through their evolutionary drift as a means to evade preexisting immunity. We evaluate here the evolution signature of glycosylations found in the globular head region of H1 HA in order to determine their impact in the virulence and transmission of H1N1 viruses. We identified a polymorphism at HA residue 147 associated with the acquisition of glycosylation at residues 144 and 172. By in vitro and in vivo analyses using mutant viruses, we also found that the polymorphism at HA residue 147 compensated for the loss of replication, virulence, and transmissibility associated with the presence of the N-linked glycans. Our findings suggest that the polymorphism in H1 HA at position 147 modulates viral fitness by buffering the constraints caused by N-linked glycans and provide insights into the evolution dynamics of influenza viruses with implications in vaccine immunogenicity.     

Antigenic structure of the influenza B virus hemagglutinin: nucleotide sequence analysis of antigenic variants selected with monoclonal antibodies.

We report here the complete nucleotide sequence of the hemagglutinin (HA) gene of influenza B virus B/Oregon/5/80 and, through comparative sequence analysis, identify amino acid substitutions in the HA1 polypeptide responsible for the antigenic alterations in laboratory-selected antigenic variants of this virus. The complete nucleotide sequence of the B/Oregon/5/80 HA gene was established by a combination of chemical sequencing of a full-length cDNA clone and dideoxy sequencing of the virion RNA. The nucleotide sequence is very similar to previously reported influenza B virus HA gene sequences and differs at only nine nucleotide positions from the B/Singapore/222/79 HA gene (Verhoeyen et al., Nucleic Acids Res. 11:4703-4712, 1983). The nucleotide sequences of the HA1 portions of the HA genes of 18 laboratory-selected antigenic variants were determined by the dideoxy method. Comparison of the deduced amino acid sequences of the parental and variant HA1 polypeptides revealed 16 different amino acid substitutions at nine positions. All amino acid substitutions resulted from single-point mutations, and no double mutants were detected, demonstrating that as in the influenza A viruses, single amino acid substitutions are sufficient to alter the antigenicity of the HA molecule. Many of the amino acid substitutions in the variants occurred at positions also observed to change in natural drift strains. The substitutions appear to identify at least two immunodominant regions which correspond to proposed antigenic sites A and B on the influenza A virus H3 HA.     

Epitope mapping of the hemagglutinin molecule of a highly pathogenic H5N1 influenza virus by using monoclonal antibodies

We mapped the hemagglutinin (HA) antigenic epitopes of a highly pathogenic H5N1 influenza virus on the three-dimensional HA structure by characterizing escape mutants of a recombinant virus containing A/Vietnam/1203/04 (H5N1) ΔHA and neuraminidase genes in the genetic background of A/Puerto Rico/8/34 (H1N1) virus. The mutants were selected with a panel of eight anti-HA monoclonal antibodies (MAbs), seven to A/Vietnam/1203/04 (H5N1) virus and one to A/Chicken/Pennsylvania/8125/83 (H5N2) virus, and the mutants’ HA genes were sequenced. The amino acid changes suggested three MAb groups: four MAbs reacted with the complex epitope comprising parts of the antigenic site B of H3 HA and site Sa of H1 HA, two MAbs reacted with the epitope corresponding to the antigenic site A in H3 HA, and two MAbs displayed unusual behavior: each recognized amino acid changes at two widely separate antigenic sites. Five changes were detected in amino acid residues not previously reported as changed in H5 escape mutants, and four others had substitutions not previously described. The HA antigenic structure differs substantially between A/Vietnam/1203/04 (H5N1) virus and the low-pathogenic A/Mallard/Pennsylvania/10218/84 (H5N2) virus we previously characterized (N. V. Kaverin et al., J. Gen. Virol. 83:2497-2505, 2002). The hemagglutination inhibition reactions of the MAbs with recent highly pathogenic H5N1 viruses were consistent with the antigenic-site amino acid changes but not with clades and subclades based on H5 phylogenetic analysis. These results provide information on the recognition sites of the MAbs widely used to study H5N1 viruses and demonstrate the involvement of the HA antigenic sites in the evolution of highly pathogenic H5N1 viruses, findings that can be critical for characterizing pathogenesis and vaccine design.     

Predicting the Antigenic Structure of the Pandemic (H1N1) 2009 Influenza Virus Hemagglutinin

The pandemic influenza virus (2009 H1N1) was recently introduced into the human population. The hemagglutinin (HA) gene of 2009 H1N1 is derived from “classical swine H1N1” virus, which likely shares a common ancestor with the human H1N1 virus that caused the pandemic in 1918, whose descendant viruses are still circulating in the human population with highly altered antigenicity of HA. However, information on the structural basis to compare the HA antigenicity among 2009 H1N1, the 1918 pandemic, and seasonal human H1N1 viruses has been lacking. By homology modeling of the HA structure, here we show that HAs of 2009 H1N1 and the 1918 pandemic virus share a significant number of amino acid residues in known antigenic sites, suggesting the existence of common epitopes for neutralizing antibodies cross-reactive to both HAs. It was noted that the early human H1N1 viruses isolated in the 1930s–1940s still harbored some of the original epitopes that are also found in 2009 H1N1. Interestingly, while 2009 H1N1 HA lacks the multiple N-glycosylations that have been found to be associated with an antigenic change of the human H1N1 virus during the early epidemic of this virus, 2009 H1N1 HA still retains unique three-codon motifs, some of which became N-glycosylation sites via a single nucleotide mutation in the human H1N1 virus. We thus hypothesize that the 2009 H1N1 HA antigenic sites involving the conserved amino acids will soon be targeted by antibody-mediated selection pressure in humans. Indeed, amino acid substitutions predicted here are occurring in the recent 2009 H1N1 variants. The present study suggests that antibodies elicited by natural infection with the 1918 pandemic or its early descendant viruses play a role in specific immunity against 2009 H1N1, and provides an insight into future likely antigenic changes in the evolutionary process of 2009 H1N1 in the human population.     

Molecular evolution of hemagglutinin genes of H1N1 swine and human influenza A viruses

Summary The hemagglutinin (HA) genes of influenza type A (H1N1) viruses isolated from swine were cloned into plasmid vectors and their nucleotide sequences were determined. A phylogenetic tree for the HA genes of swine and human influenza viruses was constructed by the neighbor-joining method. It showed that the divergence between swine and human HA genes might have occurred around 1905. The estimated rates of synonymous (silent) substitutions for swine and human influenza viruses were almost the same. For both viruses, the rate of synonymous substitution was much higher than that of nonsynonymous (amino acid altering) substitution. It is the case even for only the antigenic sites of the HA. This feature is consistent with the neutral theory of molecular evolution. The rate of nonsynonymous substitution for human influenza viruses was three times the rate for swine influenza viruses. In particular, nonsynonymous substitutions at antigenic sites occurred less frequently in swine than in humans. The difference in the rate of nonsynonymous substitution between swine and human influenza viruses can be explained by the different degrees of functional constraint operating on the amino acid sequence of the HA in both hosts.     

2017年河北省蛋鸡养殖场H9N2亚型禽流感病毒的遗传演化和抗原性分析

【背景】H9N2亚型禽流感病毒在鸡群中广泛流行,引起巨大损失。【目的】了解河北省蛋鸡养殖场H9N2亚型禽流感病毒(avian influenza virus,AIV)的基因序列和抗原性的变异情况,为该病原的科学防控提供理论依据。【方法】于2017年从河北省部分蛋鸡养殖场分离鉴定出7株H9N2亚型AIV,对其HA基因进行序列测定,并进行遗传演化、关键氨基酸位点及抗原性分析。【结果】7株分离毒株HA基因同源性在95.5%?97.2%之间;与2016年前的流行毒株相比,分离病毒HA裂解位点均为典型低致病性AIV特征,在受体结合区域出现变异,潜在糖基化位点无明显差异;抗原分析结果显示分离毒株与早期分离株相比抗原性发生了变异,形成了新的抗原群;抗原性相关位点分析显示,分离毒株在9个位点发生了较为明显的突变,可能是导致抗原性变异的分子基础。【结论】河北省蛋鸡养殖场H9N2亚型AIV中的流行毒株在关键功能区发生基因突变,并且抗原性发生变异,提示应持续监测H9N2亚型AIV的遗传变异情况,并及时更换疫苗株。     

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.     

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.     

Accumulation of amino acid substitutions promotes irreversible structural changes in the hemagglutinin of human influenza AH3 virus during evolution

In order to clarify the effect of an accumulation of amino acid substitutions on the hemadsorption character of the influenza AH3 virus hemagglutinin (HA) protein, we introduced single-point amino acid changes into the HA1 domain of the HA proteins of influenza viruses isolated in 1968 (A/Aichi/2/68) and 1997 (A/Sydney/5/97) by using PCR-based random mutation or site-directed mutagenesis. These substitutions were classified as positive or negative according to their effects on the hemadsorption activity. The rate of positive substitutions was about 50% for both strains. Of 44 amino acid changes that were identical in the two strains with regard to both the substituted amino acids and their positions in the HA1 domain, 22% of the changes that were positive in A/Aichi/2/68 were negative in A/Sydney/5/97 and 27% of the changes that were negative in A/Aichi/2/68 were positive in A/Sydney/5/97. A similar discordance rate was also seen for the antigenic sites. These results suggest that the accumulation of amino acid substitutions in the HA protein during evolution promoted irreversible structural changes and therefore that antigenic changes in the H3HA protein may not be limited.     

Inferring natural selection operating on conservative and radical substitution at single amino acid sites

Natural selection operating on amino acid substitution at single amino acid sites can be detected by comparing the rates of synonymous (r(S)) and nonsynonymous (r(N)) nucleotide substitution at single codon sites. Amino acid substitutions can be classified as conservative or radical according to whether they retain the properties of the substituted amino acid. Here methods for comparing the rates of conservative (r(C)) and radical (r(R)) nonsynonymous substitution with r(S) at single codon sites were developed to detect natural selection operating on these substitutions at single amino acid sites. A method for comparing r(C) and r(R) at single codon sites was also developed to detect biases toward these substitutions at single amino acid sites. Charge was used as the property of the amino acids. In a computer simulation, false-positive rates of these methods were always < 5%, unless termination sites were included in the computation of the numbers of sites and estimates of transition/transversion rate ratio were highly biased. The frequency of detection of natural selection operating on conservative substitution was almost independent of the presence of natural selection operating on radical substitution, and vice versa. Natural selection operating specifically on conservative and radical substitution was detected more efficiently by comparing r(S) with r(C) and r(S) with r(R) than by comparing r(S) with r(N). These methods also appeared to be robust against the occurrence of recombination during evolution. In an analysis of class I human leukocyte antigen, negative selection operating on conservative substitution, but not positive selection operating on radical substitution, was observed at some of the codon sites with r(R) > r(C), suggesting that r(R) > r(C) may not necessarily be an indicator of positive selection operating on radical substitution.     

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.     

江苏首例甲型H1N1(2009)流感病毒分离及其血凝素分子遗传特征分析

2009年6月12日,江苏确诊首例甲型H1N1(2009)病例。通过细胞和鸡胚分离系统,我们分离到一株具有较高血凝活性的病毒,命名为A/Jiangsu/1/2009。为了跟踪病毒的变异情况,我们开展了病毒的全基因组测序工作,在此基础上对其血凝素基因(Haemagglutinin,HA)的遗传特性进行了详细研究。分离株HA蛋白不具有多碱基HA裂解位点,具有低致病性流感病毒特点。与参考株A/California/04/2009相比,分离株A/Jiangsu/1/2009HA蛋白的有4个氨基酸发生了突变,但都不在已知的抗原位点上。分离株有5个潜在糖基化位点,这与近年来古典猪H1N1和北美三源重配猪H1病毒完全一致,保留了古典猪H1的特点。与禽流感H1病毒相比,分离株HA蛋白受体结合位点上的E190D和G225D发生突变,这可能成为新甲型H1N1(2009)在人际间传播的一个重要分子基础。此外,其它受体结合位点上相关氨基酸同时具有人和猪流感病毒的特点。本研究首次对早期流行的甲型H1N1(2009)流感病毒的HA蛋白的分子遗传特征进行了详细研究,对进一步监测病原变异具有重要指导意义。     

Directionality in the evolution of influenza A haemagglutinin

The evolution of haemagglutinin (HA), an important influenza virus antigen, has been the subject of intensive research for more than two decades. Many characteristics of HA's sequence evolution are captured by standard Markov chain substitution models. Such models assign equal fitness to all accessible amino acids at a site. We show, however, that such models strongly underestimate the number of homoplastic amino acid substitutions during the course of HA's evolution, i.e. substitutions that repeatedly give rise to the same amino acid at a site. We develop statistics to detect individual homoplastic events and find that they preferentially occur at positively selected epitopic sites. Our results suggest that the evolution of the influenza A HA, including evolution by positive selection, is strongly affected by the long-term site-specific preferences for individual amino acids.     

Compensatory evolution of net-charge in influenza A virus hemagglutinin

The propagation of influenza A virus depends on the balance between the activities of hemagglutinin (HA) for binding to host cells and neuraminidase (NA) for releasing from infected cells (HA-NA balance). Since the host cell membrane and the sialic acid receptor are negatively charged, the amino acid substitutions increasing (charge+) and decreasing (charge-) the positive charge of HA subunit 1 (HA1) enhance and reduce, respectively, the binding avidity and affinity. The positive charge of HA1 in human influenza A virus bearing subtype H3N2 (A/H3N2 virus) was observed to have increased during evolution, but the evolutionary mechanism for this observation was unclear because this may disrupt the HA-NA balance. Here we show, from the phylogenetic analysis of HA for human A/H3N2 and A/H1N1 viruses, that the relative frequencies of charge+ and charge- substitutions were elevated on the branches where the number of N-glycosylation sites (NGS) increased and decreased, respectively, compared to those where the number of NGS did not change. On the latter branches, the net-charge of HA1 appeared to have been largely maintained to preserve its structure and function. Since the charge+ and charge- substitutions in HA1 have opposite effects to the gain and loss of NGS on the binding and release of the virus, the net-charge of HA1 may have evolved to compensate for the effect of the gain and loss of NGS, probably through changing the avidity. Apparently, the relative frequency of charge- substitutions in HA1 of A/H3N2 virus was elevated after the introduction of oseltamivir, and that of charge+ substitutions in HA1 of A/H1N1 virus was elevated after the spread of oseltamivir resistance. These observations may also be explained by the compensatory effect of the net-charge in HA1 on the NA activity for keeping the HA-NA balance.     

Evidence for N-glycan shielding of antigenic sites during evolution of human influenza A virus hemagglutinin

After the emergence of influenza A viruses in the human population, the number of N-glycosylation sites (NGS) in the globular head region of hemagglutinin (HA) has increased continuously for several decades. It has been speculated that the addition of NGS to the globular head region of HA has conferred selective advantages to the virus by preventing the binding of antibodies (Ab) to antigenic sites (AS). Here, the effect of N-glycosylation on the binding of Ab to AS in human influenza A virus subtype H3N2 (A/H3N2) was examined by inferring natural selection at AS and other sites (NAS) that are located close to and distantly from the NGS in the three-dimensional structure of HA through a comparison of the rates of synonymous (d(S)) and nonsynonymous (d(N)) substitutions. When positions 63, 122, 126, 133, 144, and 246 in the globular head region of HA were non-NGS, the d(N)/d(S) was >1 and positive selection was detected at the AS located near these positions. However, the d(N)/d(S) value decreased and the evidence of positive selection disappeared when these positions became NGS. In contrast, d(N)/d(S) at the AS distantly located from the positions mentioned above and at the NAS of any location were generally <1 and did not decrease when these positions changed from non-NGS to NGS. These results suggest that the attachment of N-glycans to the NGS in the globular head region of HA prevented the binding of Ab to AS in the evolutionary history of human A/H3N2 virus.     

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.     

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.     

Comparative characterization of the glycosylation profiles of an influenza hemagglutinin produced in plant and insect hosts

The main objective of this study was to characterize the N-linked glycosylation profiles of recombinant hemagglutinin (HA) proteins expressed in either insect or plant hosts, and to develop a mass spectrometry based workflow that can be used in quality control to assess batch-to-batch reproducibility for recombinant HA glycosylation. HA is a surface glycoprotein of the influenza virus that plays a key role in viral infectivity and pathogenesis. Characterization of the glycans for plant recombinant HA from the viral strain A/California/04/09 (H1N1) has not yet been reported. In this study, N-linked glycosylation patterns of the recombinant HAs from both insect and plant hosts were characterized by precursor ion scan-driven data-dependent analysis followed by high-resolution MS/MS analysis of the deglycosylated tryptic peptides. Five glycosylation sites (N11, N23, N276, N287, and N481) were identified containing high mannose type glycans in plant-expressed HAs, and complex type glycoforms for the insect-expressed HA. More than 95% site occupancy was observed for all glycosylation sites except N11, which was 60% occupied. Multiple-reaction monitoring based quantitation analysis was developed for each glycopeptide isoform and the quantitative results indicate that the Man(8) GlcNAc(2) is the dominant glycan for all sites in plant-expressed HAs. The relative abundance of the glycoforms at each specific glycosylation site and the relative quantitation for each glycoform among three HAs were determined. Few differences in the glycosylation profiles were detected between the two batches of plant HAs studied, but there were significant differences between the glycosylation patterns in the HAs generated in plant and insect expression hosts.     

Changing selective pressure during antigenic changes in human influenza H3

The rapid evolution of influenza viruses presents difficulties in maintaining the optimal efficiency of vaccines. Amino acid substitutions result in antigenic drift, a process whereby antisera raised in response to one virus have reduced effectiveness against future viruses. Interestingly, while amino acid substitutions occur at a relatively constant rate, the antigenic properties of H3 move in a discontinuous, step-wise manner. It is not clear why this punctuated evolution occurs, whether this represents simply the fact that some substitutions affect these properties more than others, or if this is indicative of a changing relationship between the virus and the host. In addition, the role of changing glycosylation of the haemagglutinin in these shifts in antigenic properties is unknown. We analysed the antigenic drift of HA1 from human influenza H3 using a model of sequence change that allows for variation in selective pressure at different locations in the sequence, as well as at different parts of the phylogenetic tree. We detect significant changes in selective pressure that occur preferentially during major changes in antigenic properties. Despite the large increase in glycosylation during the past 40 years, changes in glycosylation did not correlate either with changes in antigenic properties or with significantly more rapid changes in selective pressure. The locations that undergo changes in selective pressure are largely in places undergoing adaptive evolution, in antigenic locations, and in locations or near locations undergoing substitutions that characterise the change in antigenicity of the virus. Our results suggest that the relationship of the virus to the host changes with time, with the shifts in antigenic properties representing changes in this relationship. This suggests that the virus and host immune system are evolving different methods to counter each other. While we are able to characterise the rapid increase in glycosylation of the haemagglutinin during time in human influenza H3, an increase not present in influenza in birds, this increase seems unrelated to the observed changes in antigenic properties.