In July 2004, an outbreak of influenza A (H3N2) was detected at 3 Bhutanese refugee camps in southeastern Nepal. Hemagglutination inhibition showed that approximately 40% of the viruses from this outbreak were antigenically distinct from the A/Wyoming/3/03 vaccine strain. Four amino acid differences were observed in most of the 26 isolates compared with the A/Wyoming/3/2003 vaccine strain. All 4 substitutions are located within or adjacent to known antibody-binding sites. Several isolates showed a lysine-to-asparagine substitution at position 145 (K145N) in the hemagglutinin molecule, which may be noteworthy since position 145 is located within a glycosylation site and adjacent to an antibody-binding site. H3N2 viruses continue to drift from the vaccine strain and may remain as the dominant strains during the 2005-2006 influenza season. Thus, the 2005-2006 Northern Hemisphere vaccine strain was changed to A/California/7/2004, a virus with all 4 amino acid substitutions observed in these Nepalese isolates.

BACKGROUND: Influenza A viruses are divided into subtypes based on their hemagglutinin (H1 to H15) and neuraminidase (N1 to N9) glycoproteins. Of these, three A subtypes H1N1, H3N2 and H1N2 circulate in the human population. Influenza A viruses display a high antigenic variability called "antigenic drift" which allows the virus to escape antibody neutralization. OBJECTIVES: Evaluate the mutations apparition that might predict a divergent antigenic evolution of hemagglutinin in influenza A H1N1 and A H1N2 viruses. STUDY DESIGN: During the three winters of 2001-2002 to 2003-2004, 58 A H1N1 and 23 A H1N2 subtypes have been isolated from patients with influenza-like illness in the south of France. The HA1 region was analyzed by RT-PCR and subsequently sequenced to compare the HA1 genetic evolution of influenza A H1N1 and A H1N2 subtypes. RESULTS: Our results showed that 28 amino acid substitutions have accumulated in the HA1 region since the circulation of A/New Caledonia/20/99-like viruses in France. Of these, fifteen were located in four antigenic sites (B, C, D and E). Six of them were observed only in the A H1N2 isolates, six only in the A H1N1 isolates and three in both subtypes. Furthermore, nine of twenty two A H1N2 isolates from the winter of 2002-2003 shared a T90A amino acid change which has not been observed in any A H1N1 isolate; resulting in the introduction of a new glycosylation site close to the antigenic site E. This might mask some antigenic E determinants and therefore, modify the A H1N2 antigenicity. CONCLUSIONS: The divergent genetic evolution of hemagglutinin may ultimately lead to a significant different antigenicity between A H1N1 and A H1N2 subtypes that would require the introduction of a new subtype in the vaccine batches.

Although vaccine-induced immunity to influenza A virus is continually challenged by progressively selected mutations in the virus?s major antigens (antigenic drift), virus strains within a subtype (e.g., H1N1) are antigenically cross-reactive. Although cross-immunity diminishes as further mutations accumulate, necessitating frequent changes in vaccine strains, older vaccines are usually partially protective. The post-World War II epidemic of 1947 is notable for the total failure of a vaccine previously effective in the 1943-44 and 1944-45 seasons. We have combined extensive antigenic characterization of the hemagglutinin and neuraminidase antigens of the 1943 and 1947 viruses with analysis of their nucleotide and amino acid sequences and have found marked antigenic and amino acid differences in viruses of the two years. Furthermore, in a mouse model, vaccination with the 1943 vaccine had no effect on infection with the 1947 strain. These findings are important, because complete lack of cross-immunogenicity has been found previously only with antigenic shift, in which antigenically novel antigens have been captured by reassortment of human and animal strains, sometimes leading to pandemics. Although the 1947 epidemic lacked the usual hallmarks of pandemic disease, including an extensive increase in mortality, it warns of the possibility that extreme intrasubtypic antigenic variation (if coupled with an increase in disease severity) could produce pandemic disease without the introduction of animal virus antigens.

We have determined the structure of the HA of an avian influenza virus, A/duck/Ukraine/63, a member of the same antigenic subtype, H3, as the virus that caused the 1968 Hong Kong influenza pandemic, and a possible progenitor of the pandemic virus. We find that structurally significant differences between the avian and the human HAs are restricted to the receptor-binding site particularly the substitutions Q226L and G228S that cause the site to open and residues within it to rearrange, including the conserved residues Y98, W153, and H183. We have also analyzed complexes formed by the HA with sialopentasaccharides in which the terminal sialic acid is in either alpha2,3- or alpha2,6-linkage to galactose. Comparing the structures of complexes in which an alpha2,3-linked receptor analog is bound to the H3 avian HA or to an H5 avian HA leads to the suggestion that all avian influenza HAs bind to their preferred alpha2,3-linked receptors similarly, with the analog in a trans conformation about the glycosidic linkage. We find that alpha2,6-linked analogs are bound by both human and avian HAs in a cis conformation, and that the incompatibility of an alpha2,6-linked receptor with the alpha2,3-linkage-specific H3 avian HA-binding site is partially resolved by a small change in the position and orientation of the sialic acid. We discuss our results in relation to the mechanism of transfer of influenza viruses between species.

Phylogenetic patterns of the three polymerase (PB2, PB1 and PA) genes of a total of 20 influenza B viruses isolated during a 58 year period, 1940-1998, were analysed in detail in a parallel manner. All three polymerase genes consistently showed evolutionary divergence into two major distinct lineages and their amino acid profiles demonstrated conserved lineage-specific substitutions. Dendrogram topologies of the PB2 and PB1 genes were very similar and contrasted with that of the PA gene. It was of particular interest to reveal that even though the PA gene evolved into two major lineages, that of three recent Asian Victoria/1/87-like strains formed a branch cluster located in the same lineage as that of recent Yamagata/16/88-like isolates. Differences in the phylogenetic pathways of three polymerase genes were not only a reflection of genetic reassortment between co-circulating influenza B viruses, but also an indication that the polymerase genes were not co-evolving as a unit. As a result, comparison of the phylogenetic patterns of the three polymerase genes with previously determined patterns of the HA, NP, M and NS genes of 18 viruses defined the existence of eight distinct genome constellations. Also, similar phylogenetic profiles among the PA, NP and M genes, as well as between the PB2 and PB1 genes, were observed, suggesting possible functional interactions among these proteins. Completion of evolutionary analysis of the six internal genes and the HA gene of influenza B viruses revealed frequent genetic reassortment among co-circulating variable strains and suggested co-dependent evolution of genes.

BACKGROUND: Many RNA viruses do not have a single, representative genome but instead form a set of related variants that has been called a quasispecies. The sequence variability of such viruses presents a significant bioinformatics challenge. In order for the sequence information to be understood, the complete mutational spectrum needs to be distilled to a biologically relevant and analyzable representation. RESULTS: Here, we develop a "selection mapping" algorithm?QUASI?that identifies the positively selected variants of viral proteins. The key to the selection mapping algorithm is the identification of particular replacement mutations that are overabundant relative to silent mutations at each codon (e.g., threonine at hemagglutinin position 262). Selection mapping identifies such replacement mutations as positively selected. Conversely, selection mapping recognizes negatively selected variants as mutational "noise" (e.g., serine at hemagglutinin position 262). CONCLUSION: Selection mapping is a fundamental improvement over earlier methods (e.g., dN/dS) that identify positive selection at codons but do not identify which amino acids at these codons confer selective advantage. Using QUASI?s selection maps, we characterize the selected mutational landscapes of influenza A H3 hemagglutinin, HIV-1 reverse transcriptase, and HIV-1 gp120.

Algorithmic details to obtain maximum likelihood estimates of parameters on a large phylogeny are discussed. On a large tree, an efficient approach is to optimize branch lengths one at a time while updating parameters in the substitution model simultaneously. Codon substitution models that allow for variable nonsynonymous/synonymous rate ratios (omega = d(N)/d(S)) among sites are used to analyze a data set of human influenza virus type A hemagglutinin (HA) genes. The data set has 349 sequences. Methods for obtaining approximate estimates of branch lengths for codon models are explored, and the estimates are used to test for positive selection and to identify sites under selection. Compared with results obtained from the exact method estimating all parameters by maximum likelihood, the approximate methods produced reliable results. The analysis identified a number of sites in the viral gene under diversifying Darwinian selection and demonstrated the importance of including many sequences in the data in detecting positive selection at individual sites.