Effects of passage history and sampling bias on phylogenetic reconstruction of human influenza A evolution

In this paper we determine the extent to which host-mediated mutations and a known sampling bias affect evolutionary studies of human influenza A. Previous phylogenetic reconstruction of influenza A (H3N2) evolution using the hemagglutinin gene revealed an excess of nonsilent substitutions assigned to the terminal branches of the tree. We investigate two hypotheses to explain this observation. The first hypothesis is that the excess reflects mutations that were either not present or were at low frequency in the viral sample isolated from its human host, and that these mutations increased in frequency during passage of the virus in embryonated eggs. A set of 22 codons known to undergo such “host-mediated” mutations showed a significant excess of mutations assigned to branches attaching sequences from egg-cultured (as opposed to cell-cultured) isolates to the tree. Our second hypothesis is that the remaining excess results from sampling bias. Influenza surveillance is purposefully biased toward sequencing antigenically dissimilar strains in an effort to identify new variants that may signal the need to update the vaccine. This bias produces an excess of mutations assigned to terminal branches simply because an isolate with no close relatives is by definition attached to the tree by a relatively long branch. Simulations show that the magnitude of excess mutations we observed in the hemagglutinin tree is consistent with expectations based on our sampling protocol. Sampling bias does not affect inferences about evolution drawn from phylogenetic analyses. However, if possible, the excess caused by host-mediated mutations should be removed from studies of the evolution of influenza viruses as they replicate in their human hosts.

Pattern of mutation in the genome of influenza A virus on adaptation to increased virulence in the mouse lung: Identification of functional themes

The genetic basis for virulence in influenza virus is largely unknown. To explore the mutational basis for increased virulence in the lung, the H3N2 prototype clinical isolate, A/HK/1/68, was adapted to the mouse. Genomic sequencing provided the first demonstration, to our knowledge, that a group of 11 mutations can convert an avirulent virus to a virulent variant that can kill at a minimal dose. Thirteen of the 14 amino acid substitutions (93%) detected among clonal isolates were likely instrumental in adaptation because of their positive selection, location in functional regions, and/or independent occurrence in other virulent influenza viruses. Mutations in virulent variants repeatedly involved nuclear localization signals and sites of protein and RNA interaction, implicating them as novel modulators of virulence. Mouse-adapted variants with the same hemagglutinin mutations possessed different pH optima of fusion, indicating that fusion activity of hemagglutinin can be modulated by other viral genes. Experimental adaptation resulted in the selection of three mutations that were in common with the virulent human H5N1 isolate A/HK/156/97 and that may be instrumental in its extreme virulence. Analysis of viral adaptation by serial passage appears to provide the identification of biologically relevant mutations.

Sequence of the 1918 pandemic influenza virus nonstructural gene (NS) segment and characterization of recombinant viruses bearing the 1918 NS genes

The influenza A virus pandemic of 1918–1919 resulted in an estimated 20–40 million deaths worldwide. The hemagglutinin and neuraminidase sequences of the 1918 virus were previously determined. We here report the sequence of the A/Brevig Mission/1/18 (H1N1) virus nonstructural (NS) segment encoding two proteins, NS1 and nuclear export protein. Phylogenetically, these genes appear to be close to the common ancestor of subsequent human and classical swine strain NS genes. Recently, the influenza A virus NS1 protein was shown to be a type I IFN antagonist that plays an important role in viral pathogenesis. By using the recently developed technique of generating influenza A viruses entirely from cloned cDNAs, the hypothesis that the 1918 virus NS1 gene played a role in virulence was tested in a mouse model. In a BSL3+ laboratory, viruses were generated that possessed either the 1918 NS1 gene alone or the entire 1918 NS segment in a background of influenza A/WSN/33 (H1N1), a mouse-adapted virus derived from a human influenza strain first isolated in 1933. These 1918 NS viruses replicated well in tissue culture but were attenuated in mice as compared with the isogenic control viruses. This attenuation in mice may be related to the human origin of the 1918 NS1 gene. These results suggest that interaction of the NS1 protein with host-cell factors plays a significant role in viral pathogenesis.

Influenza hemagglutinin is spring-loaded by a metastable native conformation

Enveloped viruses enter cells by protein-mediated membrane fusion. For influenza virus, membrane fusion is regulated by the conformational state of the hemagglutinin (HA) protein, which switches from a native (nonfusogenic) structure to a fusion-active (fusogenic) conformation when exposed to the acidic environment of the cellular endosome. Here we demonstrate that destabilization of HA at neutral pH, with either heat or the denaturant urea, triggers a conformational change that is biochemically indistinguishable from the change triggered by low pH. In each case, the conformational change is coincident with induction of membrane-fusion activity, providing strong evidence that the fusogenic structure is formed. These results indicate that the native structure of HA is trapped in a metastable state and that the fusogenic conformation is released by destabilization of native structure. This strategy may be shared by other enveloped viruses, including those that enter the cell at neutral pH, and could have implications for understanding the membrane-fusion step of HIV infection.

Long term trends in the evolution of H(3) HA1 human influenza type A

We have studied the HA1 domain of 254 human influenza A(H3N2) virus genes for clues that might help identify characteristics of hemagglutinins (HAs) of circulating strains that are predictive of that strain’s epidemic potential. Our preliminary findings include the following. (i) The most parsimonious tree found requires 1,260 substitutions of which 712 are silent and 548 are replacement substitutions. (ii) The HA1 portion of the HA gene is evolving at a rate of 5.7 nucleotide substitutions/year or 5.7 × 10−3 substitutions/site per year. (iii) The replacement substitutions are distributed randomly across the three positions of the codon when allowance is made for the number of ways each codon can change the encoded amino acid. (iv) The replacement substitutions are not distributed randomly over the branches of the tree, there being 2.2 times more changes per tip branch than for non-tip branches. This result is independent of how the virus was amplified (egg grown or kidney cell grown) prior to sequencing or if sequencing was carried out directly on the original clinical specimen by PCR. (v) These excess changes on the tip branches are probably the result of a bias in the choice of strains to sequence and the detection of deleterious mutations that had not yet been removed by negative selection. (vi) There are six hypervariable codons accumulating replacement substitutions at an average rate that is 7.2 times that of the other varied codons. (vii) The number of variable codons in the trunk branches (the winners of the competitive race against the immune system) is 47 ± 5, significantly fewer than in the twigs (90 ± 7), which in turn is significantly fewer variable codons than in tip branches (175 ± 8). (viii) A minimum of one of every 12 branches has nodes at opposite ends representing viruses that reside on different continents. This is, however, no more than would be expected if one were to randomly reassign the continent of origin of the isolates. (ix) Of 99 codons with at least four mutations, 31 have ratios of non-silent to silent changes with probabilities less than 0.05 of occurring by chance, and 14 of those have probabilities <0.005. These observations strongly support positive Darwinian selection. We suggest that the small number of variable positions along the successful trunk lineage, together with knowledge of the codons that have shown positive selection, may provide clues that permit an improved prediction of which strains will cause epidemics and therefore should be used for vaccine production.

Negative-strand RNA viruses: genetic engineering and applications.

The negative-strand RNA viruses are a broad group of animal viruses that comprise several important human pathogens, including influenza, measles, mumps, rabies, respiratory syncytial, Ebola, and hantaviruses. The development of new strategies to genetically manipulate the genomes of negative-strand RNA viruses has provided us with new tools to study the structure-function relationships of the viral components and their contributions to the pathogenicity of these viruses. It is also now possible to envision rational approaches--based on genetic engineering techniques--to design live attenuated vaccines against some of these viral agents. In addition, the use of different negative-strand RNA viruses as vectors to efficiently express foreign polypeptides has also become feasible, and these novel vectors have potential applications in disease prevention as well as in gene therapy.

Effective inhibition of influenza virus production in cultured cells by external guide sequences and ribonuclease P

The polymerase (PB2) and nucleocapsid (NP) genes encoded by the genome of influenza virus are essential for replication of the virus. When synthetic genes that express RNAs for external guide sequences targeted to the mRNAs of the PB2 and NP genes are stably incorporated into mouse cells in tissue culture, infection of these cells with influenza virus is nonproductive. Endogenous RNase P cleaves the targeted influenza virus mRNAs when they are in a complex with the external guide sequences. Targeting two different mRNAs simultaneously inhibits viral particle production more efficiently than does targeting only one mRNA.

Hepatitis C virus and other Flaviviridae viruses enter cells via low density lipoprotein receptor

Endocytosis of the Flaviviridae viruses, hepatitis C virus, GB virus C/hepatitis G virus, and bovine viral diarrheal virus (BVDV) was shown to be mediated by low density lipoprotein (LDL) receptors on cultured cells by several lines of evidence: by the demonstration that endocytosis of these virus correlated with LDL receptor activity, by complete inhibition of detectable endocytosis by anti-LDL receptor antibody, by inhibition with anti-apolipoprotein E and -apolipoprotein B antibodies, by chemical methods abrogating lipoprotein/LDL receptor interactions, and by inhibition with the endocytosis inhibitor phenylarsine oxide. Confirmatory evidence was provided by the lack of detectable LDL receptor on cells known to be resistant to BVDV infection. Endocytosis via the LDL receptor was shown to be mediated by complexing of the virus to very low density lipoprotein or LDL but not high density lipoprotein. Studies using LDL receptor-deficient cells or a cytolytic BVDV system indicated that the LDL receptor may be the main but not exclusive means of cell entry of these viruses. Studies on other types of viruses indicated that this mechanism may not be exclusive to Flaviviridae but may be used by viruses that associate with lipoprotein in the blood. These findings provide evidence that the family of LDL receptors may serve as viral receptors.

A functionally defined model for the M2 proton channel of influenza A virus suggests a mechanism for its ion selectivity

The M2 protein from influenza A virus forms proton-selective channels that are essential to viral function and are the target of the drug amantadine. Cys scanning was used to generate a series of mutants with successive substitutions in the transmembrane segment of the protein, and the mutants were expressed in Xenopus laevis oocytes. The effect of the mutations on reversal potential, ion currents, and amantadine resistance were measured. Fourier analysis revealed a periodicity consistent with a four-stranded coiled coil or helical bundle. A three-dimensional model of this structure suggests a possible mechanism for the proton selectivity of the M2 channel of influenza virus.

Structural evidence for a second sialic acid binding site in avian influenza virus neuraminidases

The x-ray structure of a complex of sialic acid (Neu5Ac) with neuraminidase N9 subtype from A/tern/Australia/G70C/75 influenza virus at 4°C has revealed the location of a second Neu5Ac binding site on the surface of the enzyme. At 18°C, only the enzyme active site contains bound Neu5Ac. Neu5Ac binds in the second site in the chair conformation in a similar way to which it binds to hemagglutinin. The residues that interact with Neu5Ac at this second site are mostly conserved in avian strains, but not in human and swine strains, indicating that it has some as-yet-unknown biological function in birds.

Influenza C virus CM2 integral membrane glycoprotein is produced from a polypeptide precursor by cleavage of an internal signal sequence

The influenza C virus CM2 protein is a small glycosylated integral membrane protein (115 residues) that spans the membrane once and contains a cleavable signal sequence at its N terminus. The coding region for CM2 (CM2 ORF) is located at the C terminus of the 342-amino acid (aa) ORF of a colinear mRNA transcript derived from influenza C virus RNA segment 6. Splicing of the colinear transcript introduces a translational stop codon into the ORF and the spliced mRNA encodes the viral matrix protein (CM1) (242 aa). The mechanism of CM2 translation was investigated by using in vitro and in vivo translation of RNA transcripts. It was found that the colinear mRNA derived from influenza C virus RNA segment 6 serves as the mRNA for CM2. Furthermore, CM2 translation does not depend on any of the three in-frame methionine residues located at the beginning of CM2 ORF. Rather, CM2 is a proteolytic cleavage product of the p42 protein product encoded by the colinear mRNA: a cleavage event that involves the recognition and cleavage of an internal signal peptide presumably by signal peptidase resident in the endoplasmic reticulum. Alteration of the predicted signal peptidase cleavage site by mutagenesis blocked generation of CM2. The other polypeptide species resulting from the cleavage of p42, designated p31, contains the CM1 coding region and an additional C-terminal 17 aa (formerly the CM2 signal peptide). Protein p31, in comparison to CM1, displays characteristics of an integral membrane protein.

Studies using double mutants of the conformational transitions in influenza hemagglutinin required for its membrane fusion activity

Amino acid substitutions widely distributed throughout the influenza hemagglutinin (HA) influence the pH of its membrane fusion activity. We have combined a number of these substitutions in double mutants and determined the effects on the pH of fusion and on the pH at which the refolding of HA required for fusion occurs. By analyzing combinations of mutations in three regions of the metastable neutral-pH HA that are rearranged at fusion pH we obtain evidence for both additive and nonadditive effects and for an apparent order of dominance in the effects of amino acid substitutions in particular regions on the pH of fusion. We conclude that there are at least three components in the structural transition required for membrane fusion activity and consider possible pathways for the transition in relation to the known differences between neutral and fusion pH HA structures.

Mutation rates among RNA viruses

The rate of spontaneous mutation is a key parameter in modeling the genetic structure and evolution of populations. The impact of the accumulated load of mutations and the consequences of increasing the mutation rate are important in assessing the genetic health of populations. Mutation frequencies are among the more directly measurable population parameters, although the information needed to convert them into mutation rates is often lacking. A previous analysis of mutation rates in RNA viruses (specifically in riboviruses rather than retroviruses) was constrained by the quality and quantity of available measurements and by the lack of a specific theoretical framework for converting mutation frequencies into mutation rates in this group of organisms. Here, we describe a simple relation between ribovirus mutation frequencies and mutation rates, apply it to the best (albeit far from satisfactory) available data, and observe a central value for the mutation rate per genome per replication of μg ≈ 0.76. (The rate per round of cell infection is twice this value or about 1.5.) This value is so large, and ribovirus genomes are so informationally dense, that even a modest increase extinguishes the population.

Variable efficacy of repeated annual influenza vaccination

Conclusions have differed in studies that have compared vaccine efficacy in groups receiving influenza vaccine for the first time to efficacy in groups vaccinated more than once. For example, the Hoskins study [Hoskins, T. W., Davis, J. R., Smith, A. J., Miller, C. L. & Allchin, A. (1979) Lancet i, 33–35] concluded that repeat vaccination was not protective in the long term, whereas the Keitel study [Keitel, W. A., Cate, T. R., Couch, R. B., Huggins, L. L. & Hess, K. R. (1997) Vaccine 15, 1114–1122] concluded that repeat vaccination provided continual protection. We propose an explanation, the antigenic distance hypothesis, and test it by analyzing seven influenza outbreaks that occurred during the Hoskins and Keitel studies. The hypothesis is that variation in repeat vaccine efficacy is due to differences in antigenic distances among vaccine strains and between the vaccine strains and the epidemic strain in each outbreak. To test the hypothesis, antigenic distances were calculated from historical hemagglutination inhibition assay tables, and a computer model of the immune response was used to predict the vaccine efficacy of individuals given different vaccinations. The model accurately predicted the observed vaccine efficacies in repeat vaccinees relative to the efficacy in first-time vaccinees (correlation 0.87). Thus, the antigenic distance hypothesis offers a parsimonious explanation of the differences between and within the Hoskins and Keitel studies. These results have implications for the selection of influenza vaccine strains, and also for vaccination strategies for other antigenically variable pathogens that might require repeated vaccination.

Suppression of gene silencing: A general strategy used by diverse DNA and RNA viruses of plants

In transgenic and nontransgenic plants, viruses are both initiators and targets of a defense mechanism that is similar to posttranscriptional gene silencing (PTGS). Recently, it was found that potyviruses and cucumoviruses encode pathogenicity determinants that suppress this defense mechanism. Here, we test diverse virus types for the ability to suppress PTGS. Nicotiana benthamiana exhibiting PTGS of a green fluorescent protein transgene were infected with a range of unrelated viruses and various potato virus X vectors producing viral pathogenicity factors. Upon infection, suppression of PTGS was assessed in planta through reactivation of green fluorescence and confirmed by molecular analysis. These experiments led to the identification of three suppressors of PTGS and showed that suppression of PTGS is widely used as a counter-defense strategy by DNA and RNA viruses. However, the spatial pattern and degree of suppression varied extensively between viruses. At one extreme, there are viruses that suppress in all tissues of all infected leaves, whereas others are able to suppress only in the veins of new emerging leaves. This variation existed even between closely related members of the potexvirus group. Collectively, these results suggest that virus-encoded suppressors of gene silencing have distinct modes of action, are targeted against distinct components of the host gene-silencing machinery, and that there is dynamic evolution of the host and viral components associated with the gene-silencing mechanism.