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

During protein evolution the amino acid substitutions accumulate with time. However, the effect of accumulation of the amino acid substitutions to structural changes has not been estimated well. We will propose that the discordance of amino acid substitution on the HA protein of influenza A virus is useful for the assessment of structural changes during evolution. Discordance value can be obtained from the experimental data of tolerance or intolerance by introducing site directed mutagenesis at the homologous positions of two HA proteins holding the same amino acid residues. The value of discordance correlated to the number of amino acid differences among proteins. In the H3HA discordance rate was calculated to be 0.45% per one amino acid change. Furthermore, discordance of amino acid substitutions suggests that tolerable amino acid substitutions in different order have a probability of promoting irreversible divergence of the HA protein to different subtypes.     

Frequent and widespread parallel evolution of protein sequences

Understanding the patterns and causes of protein sequence evolution is a major challenge in evolutionary biology. One of the critical unresolved issues is the relative contribution of selection and genetic drift to the fixation of amino acid sequence differences between species. Molecular homoplasy, the independent evolution of the same amino acids at orthologous sites in different taxa, is one potential signature of selection; however, relatively little is known about its prevalence in eukaryotic proteomes. To quantify the extent and type of homoplasy among evolving proteins, we used phylogenetic methodology to analyze 8 genome-scale data matrices from clades of different evolutionary depths that span the eukaryotic tree of life. We found that the frequency of homoplastic amino acid substitutions in eukaryotic proteins was more than 2-fold higher than expected under neutral models of protein evolution. The overwhelming majority of homoplastic substitutions were parallelisms that involved the most frequently exchanged amino acids with similar physicochemical properties and that could be reached by a single-mutational step. We conclude that the role of homoplasy in shaping the protein record is much larger than generally assumed, and we suggest that its high frequency can be explained by both weak positive selection for certain substitutions and purifying selection that constrains substitutions to a small number of functionally equivalent amino acids.     

Gnarled-trunk evolutionary model of influenza A virus hemagglutinin

Human influenza A viruses undergo antigenic changes with gradual accumulation of amino acid substitutions on the hemagglutinin (HA) molecule. A strong antigenic mismatch between vaccine and epidemic strains often requires the replacement of influenza vaccines worldwide. To establish a practical model enabling us to predict the future direction of the influenza virus evolution, relative distances of amino acid sequences among past epidemic strains were analyzed by multidimensional scaling (MDS). We found that human influenza viruses have evolved along a gnarled evolutionary pathway with an approximately constant curvature in the MDS-constructed 3D space. The gnarled pathway indicated that evolution on the trunk favored multiple substitutions at the same amino acid positions on HA. The constant curvature was reasonably explained by assuming that the rate of amino acid substitutions varied from one position to another according to a gamma distribution. Furthermore, we utilized the estimated parameters of the gamma distribution to predict the amino acid substitutions on HA in subsequent years. Retrospective prediction tests for 12 years from 1997 to 2009 showed that 70% of actual amino acid substitutions were correctly predicted, and that 45% of predicted amino acid substitutions have been actually observed. Although it remains unsolved how to predict the exact timing of antigenic changes, the present results suggest that our model may have the potential to recognize emerging epidemic strains.     

Reactivity of human convalescent sera with influenza virus hemagglutinin protein mutants at antigenic site A

How the antibodies of individual convalescent human sera bind to each amino acid residue at the antigenic sites of hemagglutinin (HA) of influenza viruses, and how the antigenic drift strains of influenza viruses are selected by human sera, is not well understood. In our previous study, it was found by a binding assay with a chimeric HA between A/Kamata/14/91 (Ka/91) and A/Aichi/2/68 that convalescent human sera, following Ka/91 like (H3N2) virus infection, bind to antigenic site A of Ka/91 HA. Here using chimeric HAs possessing single amino acid substitutions at site A, it was determined how those human sera recognize each amino acid residue at antigenic site A. It was found that the capacity of human sera to recognize amino acid substitutions at site A differs from one person to another and that some amino acid substitutions result in all convalescent human sera losing their binding capacity. Among these amino acid substitutions, certain ones might be selected by chance, thus creating successive antigenic drift. Phylogenetic analysis of the drift strains of Ka/91 showed amino acid substitutions at positions 133, 135 and 145 were on the main stream of the phylogenetic tree. Indeed, all of the investigated convalescent sera failed to recognize one of them.     

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.     

新发甲型H1N1流感病毒HA分子的变异分析

目的:从分子进化水平上分析流感的起源及发展问题,研究目前爆发的H1N1病毒的HA分子的变异行为.方法:以GenBank公布的甲型流感H1N1病毒血凝素(hemagglutinin,HA)核酸序列和我国及世界范围内近几年来报告的H1N1流感病毒HA的核酸及氨基酸序列为研究对象,利用CLUSTAL 1.83和NetNGlyc 1.0等生物信息学软件对HA核酸和氨基酸序列进行了比对分析;将其糖基化位点、氨基酸序列和抗原决定簇与以往流感病毒进行了比较.同时,还将人源和猪源甲型H1N1流感病毒的HA氨基酸序列进行了序列比对和系统发育分析.结果:最新爆发的甲型H1N1流感病毒的HA除了在60,259,453,512位点高度保守区域与之前爆发的流感病毒一致外,在249位点新出现1个"-NTT-"的糖基化位点.发现所有的甲型病毒的氨基酸序列在8个氨基酸位点均发生改变,而8个氨基酸位点位于6个抗原抗原决定簇上.结论:糖基化位点的增加,氨基酸位点的改变导致抗原决定簇的改变,即抗原性漂移现象,都成为引起其传染性改变的重要原因.     

Comparison of epitope structures of H3HAs through protein modeling of influenza A virus hemagglutinin: mechanism for selection of antigenic variants in the presence of a monoclonal antibody

Starting with nine plaques of influenza A/Kamata/14/91(H3N2) virus, we selected mutants in the presence of monoclonal antibody 203 (mAb203). In total, amino acid substitutions were found at nine positions (77, 80, 131, 135, 141, 142, 143, 144 and 146), which localized in the antigenic site A of the hemagglutinin (HA). The escape mutants differed in the extent to which they had lost binding to mAb203. HA protein with substitutions of some amino acid residues created by site-directed mutagenesis in the escape mutants retained the ability to bind to mAb203. Changes in the amino acid character affecting charge or hydrophobicity accounted for the binding capacity to the antibody of the HA with most of the substitutions in the escape mutants and binding-positive mutants. However, the effect of some amino acid substitutions remained unexplained. A three-dimensional model of the 1991 HA was constructed and used to analyze substituted amino acids in these mutants for the accessible surface hydrophobic and hydrophilic characters. One amino acid substitution in an escape mutant and another amino acid substitution in a binding-positive mutant seemed to be explained by the changes noted on this model.     

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.     

Role of Substitutions in the Hemagglutinin in the Emergence of the 1968 Pandemic Influenza Virus

Hemagglutinin (HA) of H3N2/1968 pandemic influenza viruses differs from the putative avian precursor by seven amino acid substitutions. Substitutions Q226L and G228S are known to be essential for adaptation of avian HA to mammals. We found that introduction of avian-virus-like amino acids at five other HA positions (positions 62, 81, 92, 144, and 193) of A/Hong Kong/1/1968 virus decreased viral replication in human cells and transmission in pigs. Thus, substitutions at some of these positions facilitated emergence of the pandemic virus.     

2009−2015年中国甲型H1N1流感病毒血凝素基因进化分析

【目的】为了解中国地区2009?2015年甲型H1N1流感病毒流行态势,分析血凝素(Hemagglutinin,HA)基因的变异情况及其遗传进化特征。【方法】汇集国家流感中心2009?2015年流感周报的流感流行数据,分析甲型H1N1流感的流行病学特征;从全球共享禽流感数据倡议组织数据库及美国国家生物技术中心数据库下载甲型H1N1流感病毒HA基因序列,采用生物学软件进行系统进化和遗传特性的分析。【结果】2009?2015年全国共发生4次甲型H1N1流感的流行高峰。2009?2015年毒株与参考毒株A/California/07/2009(H1N1)的HA基因同源性逐年降低。遗传进化分析显示同一年份的毒株在系统进化树上基本呈现集中分布,2011年的毒株独立形成2个分支。分子特征表现为HA基因的4个抗原决定簇氨基酸位点均有变异,其中Ca区的203位、Sa区的163位和Sb区的185位氨基酸位点逐渐替换为新的氨基酸。除2010年与2012年,其他年份的毒株通过不同模型均得到正向压力选择HA氨基酸位点240。【结论】甲型H1N1流感在中国地区成为主要流行的亚型之一。HA基因与其编码的氨基酸逐年变异,未来进一步的流感监测能力还需加强。     

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.     

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.     

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.     

Distinct glycoprotein inhibitors of influenza A virus in different animal sera.

Normal horse and guinea pig sera contain the glycoprotein inhibitor alpha 2-macroglobulin, which inhibits the infectivity and hemagglutinating activity of influenza A viruses of the H2 and H3 subtypes. In the current study, the presence of inhibitors of influenza A virus in pig and rabbit sera was investigated. Variants of influenza virus type A/Los Angeles/2/87(H3N2) that were resistant to horse, pig, or rabbit serum were isolated. Analysis of the variant viruses with anti-hemagglutinin (HA) monoclonal antibodies revealed that antigenic changes occurred with the development of serum inhibitor resistance. Characterization of the inhibitors in pig and rabbit sera by using periodate and receptor-destroying enzyme demonstrated that carbohydrate is an important constituent of the active portion of both inhibitor molecules and that sialic acid is involved in the interaction of the inhibitors with influenza virus HA. Nucleotide sequence analysis of the HA molecule revealed that the serum-resistant variants each acquired a different set of amino acid alterations. The multiply resistant variants maintained the original amino acid changes and acquired additional changes. Sequence modifications in the HA involved the conserved amino acids within the receptor binding site (RBS) at position 137 and the second-shell RBS residues at positions 155 and 186. Amino acid changes also occurred within antigenic site A (position 145) and directly behind the receptor binding pocket (position 220). Amino acid alterations resulted in the acquisition of a potential glycosylation site at position 128 and the loss of potential glycosylation sites at positions 246 and 248. The localization of the amino acid changes in HA1 to the region of the RBS supports the concept of serum inhibitors as receptor analogs. The unique set of mutations acquired by the serum inhibitor-resistant variants strongly suggests that horse, pig, and rabbit sera each contain distinct glycoprotein inhibitors of influenza A virus.     

Variant influenza virus hemagglutinin that induces fusion at elevated pH.

The hemagglutinin (HA) glycoprotein of influenza virus performs two critical roles during infection: it binds virus to cell surface sialic acids, and under mildly acidic conditions it induces fusion of the virion with intracellular membranes, liberating the genome into the cytoplasm. The pH dependence of fusion varies for different influenza virus strains. Here we report the isolation and characterization of a naturally occurring variant of the X31 strain that fuses at a pH 0.2 units higher than the parent strain does and that is less sensitive to the effects of ammonium chloride, a compound known to elevate endosomal pH. The bromelain-solubilized ectodomain of the variant HA displayed a corresponding shift in the pH at which it changed conformation and bound to liposomes. Cloning and sequencing of the variant HA gene revealed amino acid substitutions at three positions in the polypeptide. Two substitutions were in antigenic determinants in the globular region of HA1, and the third occurred in HA2 near the base of the molecule. By using chimeric HA molecules expressed in CV-1 cells from simian virus 40-based vectors, we demonstrated that the change in HA2 was solely responsible for the altered fusion phenotype. This substitution, asparagine for aspartic acid at position 132, disrupted a highly conserved interchain salt bridge between adjacent HA2 subunits. The apparent role of this residue in stabilizing the HA trimer is consistent with the idea that the trimer dissociates at low pH. Furthermore, the results demonstrate that influenza virus populations contain fusion variants, raising the possibility that such variants may play a role in the evolution of the virus.     

H5N2亚型禽流感病毒血凝素裂解位点碱性氨基酸对应mRNA核苷酸的二级结构分析

以H5N2亚型禽流感病毒毒株血凝素蛋白裂解位点碱性氨基酸为研究对象,对其密码子偏好性和对应mRNA序列的折叠二级结构特点进行研究和分析。旨在探讨裂解位点氨基酸对应mRNA核苷酸片段的二级结构与病毒致病力的关系,希望能对禽流感病毒的研究提供一些基础性信息。将mRNA样本按照序列等步长递增的方法,用RNAstructure 4．1程序预测这些样本的动态延伸折叠二级结构。序列和结构的分析结果：裂解位点的碱性氨基酸对富含腺嘌呤的密码子有强烈偏好;与碱性氨基酸对应的mRNA片段上的核苷酸主要位于折叠二级结构的单链环区,少数位于配对螺旋区。结果表明：裂解位点氨基酸对应的mRNA核苷酸形成发夹端环的大小与其碱性氨基酸的多少具有正相关性。     

Structural features influencing hemagglutinin cleavability in a human influenza A virus.

The cleavability of the hemagglutinin (HA) molecule is related to the virulence of avian influenza A viruses, but its influence on human influenza virus strains is unknown. Two structural features are involved in the cleavage of avian influenza A virus HAs: a series of basic amino acids at the cleavage site and an oligosaccharide side chain in the near vicinity. The importance of these properties in the cleavability of a human influenza A virus (A/Aichi/2/68) HA was investigated by using mutants that contained or lacked an oligosaccharide side chain and had either four or six basic amino acids. All mutants except the one that contains a single mutation at the glycosylation site were cleaved, although not completely, demonstrating that a series of basic amino acids confers susceptibility to cellular cleavage enzymes among human influenza virus HAs. The mutants containing six basic amino acids at the cleavage site showed limited polykaryon formation upon exposure to low pH, indicating that cleavage was adequate to impart fusion activity to the HA. Deletion of the potential glycosylation site had no effect on the cleavability of these mutants; hence, the oligosaccharide side chain appears to have no role in human influenza virus HA cleavage. The inability to induce high cleavability in a human influenza A virus HA by insertion of a series of basic amino acids at the cleavage site indicates that other, as yet unidentified structural features are needed to enhance the susceptibility of these HAs to cellular proteases.     

Restriction of amino acid change in influenza A virus H3HA: comparison of amino acid changes observed in nature and in vitro

We introduced 248 single-point amino acid changes into hemagglutinin (HA) protein of the A/Aichi/2/68 (H3N2) strain by a PCR random mutation method. These changes were classified as positive or negative according to their effect on hemadsorption activity. We observed following results. (i) The percentage of surviving amino acid changes on the HA1 domain that did not abrogate hemadsorption activity was calculated to be ca. 44%. In nature, it is estimated to be ca. 39.6%. This difference in surviving amino acid changes on the HA protein between natural isolates and in vitro mutants might be due to the immune pressure against the former. (ii) A total of 26 amino acid changes in the in vitro mutants matched those at which mainstream amino acid changes had occurred in the H3HA1 polypeptide from 1968 to 2000. Of these, 25 were positive. We suggest that the majority of amino acid changes on the HA protein during evolution might be restricted to those that were positive on the HA of A/Aichi/2/68. (iii) We constructed two-point amino acid changes on the HA protein by using positive mutants. These two-point amino acid changes with a random combination did not inhibit hemadsorption activity. It is possible that an accumulation of amino acid change might occur without order. (iv) From the analysis of amino acids participating in mainstream amino acid change, each antigenic site could be further divided into smaller sites. The amino acid substitutions in the gaps between these smaller sites resulted in mostly hemadsorption-negative changes. These gap positions may play an important role in maintaining the function of the HA protein, and therefore amino acid changes are restricted at these locations.     

Defects in antigen presentation of mutant influenza haemagglutinins are reversed by mutations in the MHC class II molecule

A functional analysis was undertaken of the effects of mutating single amino acid residues in the alpha chain of the I-Ak molecule (to alanine; residues 50-79) on the ability of I-Ak transfectants to process and present influenza haemagglutinin to CD4+ T cell clones specific for two major antigenic sites of the HA1 subunit. In each instance, T cells were insensitive to a majority of substitutions in Ak with the exception of a few critical residues that differed for individual T cell clones. But more significantly, the failure of T cell clones to respond to mutant influenza viruses, containing drift substitutions within a T cell recognition site, in association with wild type I-Ak, could be reversed by single substitutions in Ak alpha. A T cell clone specific for HA1 120-139 failed to respond to a laboratory mutant virus (HA1 135 Gly----Arg) whereas optimal responses were observed with a mutant Ak transfectant (Ak alpha 56 Arg----Ala). Similarly, mutant transfectant 62 (Ak alpha 62 Gly----Ala) was able to present a natural variant virus A/TEX/77 to a T cell clone specific for HA1 48-67. We propose that Ak alpha 56 and Ak alpha 62 increase the affinity of association of mutant HA1 peptides for class II and therefore confer T cell recognition of variant viruses.     

Identification of amino acids of influenza virus HA responsible for resistance to a fusion inhibitor, Stachyflin

We have recently described a novel hemagglutinin (HA) conformational change inhibitor of human influenza virus, Stachyflin (Yoshimoto et al, Arch. Virol., 144, 1-14, 1999). Stachyflin-resistant variants of human influenza A/WSN/33 (H1N1) virus were isolated in vitro and the nucleotide sequences of their HA genes were determined. The relation of amino acid substitutions and Stachyflin resistance was analyzed with in vitro membrane fusion between HA-expressing cells and octadecylrhodamine (R18)-labelled chick erythrocytes (RBC). The amino acid substitutions, lysine to arginine at position 51 or lysine to glutamic acid at position 121 of the HA2 subunit of the HA protein was enough to confer a Stachyflin-resistant phenotype of HA protein. The molecular mechanism of anti-HA conformational change activity of Stachyflin is discussed.