Minor changes in the hemagglutinin of influenza A(H1N1)2009 virus alter its antigenic properties

Background

The influenza A(H1N1)2009 virus has been the dominant type of influenza A virus in Finland during the 2009–2010 and 2010–2011 epidemic seasons. We analyzed the antigenic characteristics of several influenza A(H1N1)2009 viruses isolated during the two influenza seasons by analyzing the amino acid sequences of the hemagglutinin (HA), modeling the amino acid changes in the HA structure and measuring antibody responses induced by natural infection or influenza vaccination.

Methods/Results

Based on the HA sequences of influenza A(H1N1)2009 viruses we selected 13 different strains for antigenic characterization. The analysis included the vaccine virus, A/California/07/2009 and multiple California-like isolates from 2009–2010 and 2010–2011 epidemic seasons. These viruses had two to five amino acid changes in their HA1 molecule. The mutation(s) were located in antigenic sites Sa, Ca1, Ca2 and Cb region. Analysis of the antibody levels by hemagglutination inhibition test (HI) indicated that vaccinated individuals and people who had experienced a natural influenza A(H1N1)2009 virus infection showed good immune responses against the vaccine virus and most of the wild-type viruses. However, one to two amino acid changes in the antigenic site Sa dramatically affected the ability of antibodies to recognize these viruses. In contrast, the tested viruses were indistinguishable in regard to antibody recognition by the sera from elderly individuals who had been exposed to the Spanish influenza or its descendant viruses during the early 20^th century.

Conclusions

According to our results, one to two amino acid changes (N125D and/or N156K) in the major antigenic sites of the hemagglutinin of influenza A(H1N1)2009 virus may lead to significant reduction in the ability of patient and vaccine sera to recognize A(H1N1)2009 viruses.     

Incidence of amantadine-resistant influenza A viruses in sentinel surveillance sites and nursing homes in Niigata, Japan

We surveyed the incidence of amantadine-resistant influenza A viruses both at sentinel surveillance sites and at nursing homes, and verified their types of change by partial nucleotide sequence analysis of the M2 protein. Fifty-five influenza A viruses from 27 sentinel surveillance sites during six influenza seasons from 1993 to 1999, and 26 influenza A viruses from 5 nursing homes from 1996 to 1999 were examined for susceptibility to the drug by virus titration in the presence or absence of amantadine. While amantadine-resistant viruses were not found in sentinel surveillance sites, a high frequency of resistance (8/26, 30.8%) in nursing homes was observed. Resistant viruses can occur quickly and be transmitted when used in an outbreak situation at nursing homes, where amantadine is used either for neurologic indications or for influenza treatment. Eight resistant viruses had a single amino acid change of the M2 protein at residue 30 or 31. In vitro, all 11 sensitive viruses turned resistant after 3 or 5 passages in the presence of 2 microg/ml amantadine, and they showed an amino acid change at residue 27, 30, or 31. The predominant amino acid substitution in the M2 protein of resistant viruses is Ser-31-Asp (a change at 31, serine to asparagine). The results indicate that a monitoring system for amantadine-resistant influenza viruses should be established without delay in Japan.     

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.     

Evolutionarily conserved protein sequences of influenza a viruses, avian and human, as vaccine targets

Background

Influenza A viruses generate an extreme genetic diversity through point mutation and gene segment exchange, resulting in many new strains that emerge from the animal reservoirs, among which was the recent highly pathogenic H5N1 virus. This genetic diversity also endows these viruses with a dynamic adaptability to their habitats, one result being the rapid selection of genomic variants that resist the immune responses of infected hosts. With the possibility of an influenza A pandemic, a critical need is a vaccine that will recognize and protect against any influenza A pathogen. One feasible approach is a vaccine containing conserved immunogenic protein sequences that represent the genotypic diversity of all current and future avian and human influenza viruses as an alternative to current vaccines that address only the known circulating virus strains.

Methodology/Principal Findings

Methodologies for large-scale analysis of the evolutionary variability of the influenza A virus proteins recorded in public databases were developed and used to elucidate the amino acid sequence diversity and conservation of 36,343 sequences of the 11 viral proteins of the recorded virus isolates of the past 30 years. Technologies were also applied to identify the conserved amino acid sequences from isolates of the past decade, and to evaluate the predicted human lymphocyte antigen (HLA) supertype-restricted class I and II T-cell epitopes of the conserved sequences. Fifty-five (55) sequences of 9 or more amino acids of the polymerases (PB2, PB1, and PA), nucleoprotein (NP), and matrix 1 (M1) proteins were completely conserved in at least 80%, many in 95 to 100%, of the avian and human influenza A virus isolates despite the marked evolutionary variability of the viruses. Almost all (50) of these conserved sequences contained putative supertype HLA class I or class II epitopes as predicted by 4 peptide-HLA binding algorithms. Additionally, data of the Immune Epitope Database (IEDB) include 29 experimentally identified HLA class I and II T-cell epitopes present in 14 of the conserved sequences.

Conclusions/Significance

This study of all reported influenza A virus protein sequences, avian and human, has identified 55 highly conserved sequences, most of which are predicted to have immune relevance as T-cell epitopes. This is a necessary first step in the design and analysis of a polyepitope, pan-influenza A vaccine. In addition to the application described herein, these technologies can be applied to other pathogens and to other therapeutic modalities designed to attack DNA, RNA, or protein sequences critical to pathogen function.     

Loss of a single N-linked glycan from the hemagglutinin of influenza virus is associated with resistance to collectins and increased virulence in mice

Background

Glycosylation on the globular head of the hemagglutinin (HA) protein of influenza virus acts as an important target for recognition and destruction of virus by innate immune proteins of the collectin family. This, in turn, modulates the virulence of different viruses for mice. The role of particular oligosaccharide attachments on the HA in determining sensitivity to collectins has yet to be fully elucidated.

Methods

When comparing the virulence of H3N2 subtype viruses for mice we found that viruses isolated after 1980 were highly glycosylated and induced mild disease in mice. During these studies, we were surprised to find a small plaque variant of strain A/Beijing/353/89 (Beij/89) emerged following infection of mice and grew to high titres in mouse lung. In the current study we have characterized the properties of this small plaque mutant both in vitro and in vivo.

Results

Small plaque mutants were recovered following plaquing of lung homogenates from mice infected with influenza virus seed Beij/89. Compared to wild-type virus, small plaque mutants showed increased virulence in mice yet did not differ in their ability to infect or replicate in airway epithelial cells in vitro. Instead, small plaque variants were markedly resistant to neutralization by murine collectins, a property that correlated with the acquisition of an amino acid substitution at residue 246 on the viral HA. We present evidence that this substitution was associated with the loss of an oligosaccharide glycan from the globular head of HA.

Conclusion

A point mutation in the gene encoding the HA of Beij/89 was shown to ablate a glycan attachment site. This was associated with resistance to collectins and increased virulence in mice.     

Incorporating background frequency improves entropy-based residue conservation measures

Background  

Several entropy-based methods have been developed for scoring sequence conservation in protein multiple sequence alignments. High scoring amino acid positions may correlate with structurally or functionally important residues. However, amino acid background frequencies are usually not taken into account in these entropy-based scoring schemes.     

Molecular Epidemiology of Influenza A(H1N1)pdm09 Viruses from Pakistan in 2009-2010

Background

In early 2009, a novel influenza A(H1N1) virus that emerged in Mexico and United States rapidly disseminated worldwide. The spread of this virus caused considerable morbidity with over 18000 recorded deaths. The new virus was found to be a reassortant containing gene segments from human, avian and swine influenza viruses.

Methods/Results

The first case of human infection with A(H1N1)pdm09 in Pakistan was detected on 18^th June 2009. Since then, 262 laboratory-confirmed cases have been detected during various outbreaks with 29 deaths (as of 31^st August 2010). The peak of the epidemic was observed in December with over 51% of total respiratory cases positive for influenza. Representative isolates from Pakistan viruses were sequenced and analyzed antigenically. Sequence analysis of genes coding for surface glycoproteins HA and NA showed high degree of high levels of sequence identity with corresponding genes of regional viruses circulating South East Asia. All tested viruses were sensitive to Oseltamivir in the Neuraminidase Inhibition assays.

Conclusions

Influenza A(H1N1)pdm09 viruses from Pakistan form a homogenous group of viruses. Their HA genes belong to clade 7 and show antigenic profile similar to the vaccine strain A/California/07/2009. These isolates do not show any amino acid changes indicative of high pathogenicity and virulence. It is imperative to continue monitoring of these viruses for identification of potential variants of high virulence or drug resistance.     

Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified

Background  

In recent years, model based approaches such as maximum likelihood have become the methods of choice for constructing phylogenies. A number of authors have shown the importance of using adequate substitution models in order to produce accurate phylogenies. In the past, many empirical models of amino acid substitution have been derived using a variety of different methods and protein datasets. These matrices are normally used as surrogates, rather than deriving the maximum likelihood model from the dataset being examined. With few exceptions, selection between alternative matrices has been carried out in an ad hoc manner.     

Conserved amino acid markers from past influenza pandemic strains

Background  

Finding the amino acid mutations that affect the severity of influenza infections remains an open and challenging problem. Of special interest is better understanding how current circulating influenza strains could evolve into a new pandemic strain. Influenza proteomes from distinct viral phenotype classes were searched for class specific amino acid mutations conserved in past pandemics, using reverse engineered linear classifiers.     

Long intervals of stasis punctuated by bursts of positive selection in the seasonal evolution of influenza A virus

Background  

The interpandemic evolution of the influenza A virus hemagglutinin (HA) protein is commonly considered a paragon of rapid evolutionary change under positive selection in which amino acid replacements are fixed by virtue of their effect on antigenicity, enabling the virus to evade immune surveillance.     

A context-free encoding scheme of protein sequences for predicting antigenicity of diverse influenza A viruses

Background

The evolution of influenza A viruses leads to the antigenic changes. Serological diagnosis of the antigenicity is usually labor-intensive, time-consuming and not suitable for early-stage detection. Computational prediction of the antigenic relationship between emerging and old strains of influenza viruses using viral sequences can facilitate large-scale antigenic characterization, especially for those viruses requiring high biosafety facilities, such as H5 and H7 influenza A viruses. However, most computational models require carefully designed subtype-specific features, thereby being restricted to only one subtype.

Methods

In this paper, we propose a Context-FreeEncoding Scheme (CFreeEnS) for pairs of protein sequences, which encodes a protein sequence dataset into a numeric matrix and then feeds the matrix into a downstream machine learning model. CFreeEnS is not only free from subtype-specific selected features but also able to improve the accuracy of predicting the antigenicity of influenza. Since CFreeEnS is subtype-free, it is applicable to predicting the antigenicity of diverse influenza subtypes, hopefully saving the biologists from conducting serological assays for highly pathogenic strains.

Results

The accuracy of prediction on each subtype tested (A/H1N1, A/H3N2, A/H5N1, A/H9N2) is over 85%, and can be as high as 91.5%. This outperforms existing methods that use carefully designed subtype-specific features. Furthermore, we tested the CFreeEnS on the combined dataset of the four subtypes. The accuracy reaches 84.6%, much higher than the best performance 75.1% reported by other subtype-free models, i.e. regional band-based model and residue-based model, for predicting the antigenicity of influenza. Also, we investigate the performance of CFreeEnS when the model is trained and tested on different subtypes (i.e. transfer learning). The prediction accuracy using CFreeEnS is 84.3% when the model is trained on the A/H1N1 dataset and tested on the A/H5N1, better than the 75.2% using a regional band-based model.

Conclusions

The CFreeEnS not only improves the prediction of antigenicity on datasets with only one subtype but also outperforms existing methods when tested on a combined dataset with four subtypes of influenza viruses.

     

Selective constraints on amino acids estimated by a mechanistic codon substitution model with multiple nucleotide changes

Background

Empirical substitution matrices represent the average tendencies of substitutions over various protein families by sacrificing gene-level resolution. We develop a codon-based model, in which mutational tendencies of codon, a genetic code, and the strength of selective constraints against amino acid replacements can be tailored to a given gene. First, selective constraints averaged over proteins are estimated by maximizing the likelihood of each 1-PAM matrix of empirical amino acid (JTT, WAG, and LG) and codon (KHG) substitution matrices. Then, selective constraints specific to given proteins are approximated as a linear function of those estimated from the empirical substitution matrices.

Results

Akaike information criterion (AIC) values indicate that a model allowing multiple nucleotide changes fits the empirical substitution matrices significantly better. Also, the ML estimates of transition-transversion bias obtained from these empirical matrices are not so large as previously estimated. The selective constraints are characteristic of proteins rather than species. However, their relative strengths among amino acid pairs can be approximated not to depend very much on protein families but amino acid pairs, because the present model, in which selective constraints are approximated to be a linear function of those estimated from the JTT/WAG/LG/KHG matrices, can provide a good fit to other empirical substitution matrices including cpREV for chloroplast proteins and mtREV for vertebrate mitochondrial proteins.

Conclusions/Significance

The present codon-based model with the ML estimates of selective constraints and with adjustable mutation rates of nucleotide would be useful as a simple substitution model in ML and Bayesian inferences of molecular phylogenetic trees, and enables us to obtain biologically meaningful information at both nucleotide and amino acid levels from codon and protein sequences.     

Generation and Characterization of Live Attenuated Influenza A(H7N9) Candidate Vaccine Virus Based on Russian Donor of Attenuation

Background

Avian influenza A (H7N9) virus has emerged recently and continues to cause severe disease with a high mortality rate in humans prompting the development of candidate vaccine viruses. Live attenuated influenza vaccines (LAIV) are 6:2 reassortant viruses containing the HA and NA gene segments from wild type influenza viruses to induce protective immune responses and the six internal genes from Master Donor Viruses (MDV) to provide temperature sensitive, cold-adapted and attenuated phenotypes.

Methodology/Principal Findings

LAIV candidate A/Anhui/1/2013(H7N9)-CDC-LV7A (abbreviated as CDC-LV7A), based on the Russian MDV, A/Leningrad/134/17/57 (H2N2), was generated by classical reassortment in eggs and retained MDV temperature-sensitive and cold-adapted phenotypes. CDC-LV7A had two amino acid substitutions N123D and N149D (H7 numbering) in HA and one substitution T10I in NA. To evaluate the role of these mutations on the replication capacity of the reassortants in eggs, the recombinant viruses A(H7N9)RG-LV1 and A(H7N9)RG-LV2 were generated by reverse genetics. These changes did not alter virus antigenicity as ferret antiserum to CDC-LV7A vaccine candidate inhibited hemagglutination by homologous A(H7N9) virus efficiently. Safety studies in ferrets confirmed that CDC-LV7A was attenuated compared to wild-type A/Anhui/1/2013. In addition, the genetic stability of this vaccine candidate was examined in eggs and ferrets by monitoring sequence changes acquired during virus replication in the two host models. No changes in the viral genome were detected after five passages in eggs. However, after ten passages additional mutations were detected in the HA gene. The vaccine candidate was shown to be stable in the ferret model; post-vaccination sequence data analysis showed no changes in viruses collected in nasal washes present at day 5 or day 7.

Conclusions/Significance

Our data indicate that the A/Anhui/1/2013(H7N9)-CDC-LV7A reassortant virus is a safe and genetically stable candidate vaccine virus that is now available for distribution by WHO to vaccine manufacturers.     

Frequency of D222G and Q223R hemagglutinin mutants of pandemic (H1N1) 2009 influenza virus in Japan between 2009 and 2010

Background

In April 2009, a novel swine-derived influenza A virus (H1N1pdm) emerged and rapidly spread around the world, including Japan. It has been suggested that the virus can bind to both 2,3- and 2,6-linked sialic acid receptors in infected mammals, in contrast to contemporary seasonal H1N1 viruses, which have a predilection for 2,6-linked sialic acid.

Methods/Results

To elucidate the existence and transmissibility of α2,3 sialic acid-specific viruses in H1N1pdm, amino acid substitutions within viral hemagglutinin molecules were investigated, especially D187E, D222G, and Q223R, which are related to a shift from human to avian receptor specificity. Samples from individuals infected during the first and second waves of the outbreak in Japan were examined using a high-throughput sequencing approach. In May 2009, three specimens from mild cases showed D222G and/or Q223R substitutions in a minor subpopulation of viruses infecting these individuals. However, the substitutions almost disappeared in the samples from five mild cases in December 2010. The D187E substitution was not widespread in specimens, even in May 2009.

Conclusions

These results suggest that α2,3 sialic acid-specific viruses, including G222 and R223, existed in humans as a minor population in the early phase of the pandemic, and that D222 and Q223 became more dominant through human-to-human transmission during the first and second waves of the epidemic. These results are consistent with the low substitution rates identified in seasonal H1N1 viruses in 2008.     

A class frequency mixture model that adjusts for site-specific amino acid frequencies and improves inference of protein phylogeny

Background  

Widely used substitution models for proteins, such as the Jones-Taylor-Thornton (JTT) or Whelan and Goldman (WAG) models, are based on empirical amino acid interchange matrices estimated from databases of protein alignments that incorporate the average amino acid frequencies of the data set under examination (e.g JTT + F). Variation in the evolutionary process between sites is typically modelled by a rates-across-sites distribution such as the gamma (Γ) distribution. However, sites in proteins also vary in the kinds of amino acid interchanges that are favoured, a feature that is ignored by standard empirical substitution matrices. Here we examine the degree to which the pattern of evolution at sites differs from that expected based on empirical amino acid substitution models and evaluate the impact of these deviations on phylogenetic estimation.     

Potentials 'R'Us web-server for protein energy estimations with coarse-grained knowledge-based potentials

Background  

Knowledge-based potentials have been widely used in the last 20 years for fold recognition, protein structure prediction from amino acid sequence, ligand binding, protein design, and many other purposes. However generally these are not readily accessible online.     

Phenotypic and Genetic Characterization of Avian Influenza H5N2 Viruses with Intra- and Inter-Duck Variations in Taiwan

Background

Human infections with avian influenza viruses (AIVs) have frequently raised global concerns of emerging, interspecies-transmissible viruses with pandemic potential. Waterfowl, the predominant reservoir of influenza viruses in nature, harbor precursors of different genetic lineages that have contributed to novel pandemic influenza viruses in the past.

Methods

Two duck influenza H5N2 viruses, DV518 and DV413, isolated through virological surveillance at a live-poultry market in Taiwan, showed phylogenetic relatedness but exhibited different replication capabilities in mammalian Madin-Darby Canine Kidney (MDCK) cells. This study characterizes the replication properties of the two duck H5N2 viruses and the determinants involved.

Results

The DV518 virus replicated more efficiently than DV413 in both MDCK and chicken DF1 cells. Interestingly, the infection of MDCK cells by DV518 formed heterogeneous plaques with great differences in size [large (L) and small (S)], and the two viral strains (p518-L and p518-S) obtained from plaque purification exhibited distinguishable replication kinetics in MDCK cells. Nonetheless, both plaque-purified DV518 strains still maintained their growth advantages over the plaque-purified p413 strain. Moreover, three amino acid substitutions in PA (P224S), PB2 (E72D), and M1 (A128T) were identified in intra-duck variations (p518-L vs p518-S), whereas other changes in HA (N170D), NA (I56T), and NP (Y289H) were present in inter-duck variations (DV518 vs DV413). Both p518-L and p518-S strains had the N170D substitution in HA, which might be related to their greater binding to MDCK cells. Additionally, polymerase activity assays on 293T cells demonstrated the role of vRNP in modulating the replication capability of the duck p518-L viruses in mammalian cells.

Conclusion

These results demonstrate that intra-host phenotypic variation occurs even within an individual duck. In view of recent human infections by low pathogenic AIVs, this study suggests possible determinants involved in the stepwise selection of virus variants from the duck influenza virus population which may facilitate inter-species transmission.     

Classification of viral zoonosis through receptor pattern analysis

Background  

Viral zoonosis, the transmission of a virus from its primary vertebrate reservoir species to humans, requires ubiquitous cellular proteins known as receptor proteins. Zoonosis can occur not only through direct transmission from vertebrates to humans, but also through intermediate reservoirs or other environmental factors. Viruses can be categorized according to genotype (ssDNA, dsDNA, ssRNA and dsRNA viruses). Among them, the RNA viruses exhibit particularly high mutation rates and are especially problematic for this reason. Most zoonotic viruses are RNA viruses that change their envelope proteins to facilitate binding to various receptors of host species. In this study, we sought to predict zoonotic propensity through the analysis of receptor characteristics. We hypothesized that the major barrier to interspecies virus transmission is that receptor sequences vary among species--in other words, that the specific amino acid sequence of the receptor determines the ability of the viral envelope protein to attach to the cell.     

A novel peptide ELISA for universal detection of antibodies to human H5N1 influenza viruses

Background

Active serologic surveillance of H5N1 highly pathogenic avian influenza (HPAI) virus in humans and poultry is critical to control this disease. However, the need for a robust, sensitive and specific serologic test for the rapid detection of antibodies to H5N1 viruses has not been met.

Methodology/Principal Findings

Previously, we reported a universal epitope (CNTKCQTP) in H5 hemagglutinin (HA) that is 100% conserved in H5N1 human isolates and 96.9% in avian isolates. Here, we describe a peptide ELISA to detect antibodies to H5N1 virus by using synthetic peptide that comprises the amino acid sequence of this highly conserved and antigenic epitope as the capture antigen. The sensitivity and specificity of the peptide ELISA were evaluated using experimental chicken antisera to H5N1 viruses from divergent clades and other subtype influenza viruses, as well as human serum samples from patients infected with H5N1 or seasonal influenza viruses. The peptide ELISA results were compared with hemagglutinin inhibition (HI), and immunofluorescence assay and immunodot blot that utilize recombinant HA1 as the capture antigen. The peptide ELISA detected antibodies to H5N1 in immunized animals or convalescent human sera whereas some degree of cross-reactivity was observed in HI, immunofluorescence assay and immunodot blot. Antibodies to other influenza subtypes tested negative in the peptide-ELISA.

Conclusion/Significance

The peptide-ELISA based on the highly conserved and antigenic H5 epitope (CNTKCQTP) provides sensitive and highly specific detection of antibodies to H5N1 influenza viruses. This study highlighted the use of synthetic peptide as a capture antigen in rapid detection of antibodies to H5N1 in human and animal sera that is robust, simple and cost effective and is particularly beneficial for developing countries and rural areas.