Natural Selection on the Influenza Virus Genome

Detecting patches of protein sites of influenza A viruses under positive selection

Influenza A viruses are single-stranded RNA viruses capable of evolving rapidly to adapt to environmental conditions. Examples include the establishment of a virus in a novel host or an adaptation to increasing immunity within the host population due to prior infection or vaccination against a circulating strain. Knowledge of the viral protein regions under positive selection is therefore crucial for surveillance. We have developed a method for detecting positively selected patches of sites on the surface of viral proteins, which we assume to be relevant for adaptive evolution. We measure positive selection based on dN/dS ratios of genetic changes inferred by considering the phylogenetic structure of the data and suggest a graph-cut algorithm to identify such regions. Our algorithm searches for dense and spatially distinct clusters of sites under positive selection on the protein surface. For the hemagglutinin protein of human influenza A viruses of the subtypes H3N2 and H1N1, our predicted sites significantly overlap with known antigenic and receptor-binding sites. From the structure and sequence data of the 2009 swine-origin influenza A/H1N1 hemagglutinin and PB2 protein, we identified regions that provide evidence of evolution under positive selection since introduction of the virus into the human population. The changes in PB2 overlap with sites reported to be associated with mammalian adaptation of the influenza A virus. Application of our technique to the protein structures of viruses of yet unknown adaptive behavior could identify further candidate regions that are important for host–virus interaction.

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.

Is there a twelfth protein-coding gene in the genome of influenza A? A selection-based approach to the detection of overlapping genes in closely related sequences

Protein-coding genes often contain long overlapping open-reading frames (ORFs), which may or may not be functional. Current methods that utilize the signature of purifying selection to detect functional overlapping genes are limited to the analysis of sequences from divergent species, thus rendering them inapplicable to genes found only in closely related sequences. Here, we present a method for the detection of selection signatures on overlapping reading frames by using closely related sequences, and apply the method to several known overlapping genes, and to an overlapping ORF on the negative strand of segment 8 of influenza A virus (NEG8), for which the suggestion has been made that it is functional. We find no evidence that NEG8 is under selection, suggesting that the intact reading frame might be non-functional, although we cannot fully exclude the possibility that the method is not sensitive enough to detect the signature of selection acting on this gene. We present the limitations of the method using known overlapping genes and suggest several approaches to improve it in future studies. Finally, we examine alternative explanations for the sequence conservation of NEG8 in the absence of selection. We show that overlap type and genomic context affect the conservation of intact overlapping ORFs and should therefore be considered in any attempt of estimating the signature of selection in overlapping genes.

Overestimation of nonsynonymous/synonymous rate ratio by reverse-translation of aligned amino acid sequences

In the analysis of protein-coding nucleotide sequences, the ratio of the number of nonsynonymous substitutions to that of synonymous substitutions (d(N)/d(S)) is used as an indicator for the direction and magnitude of natural selection operating at the amino acid sequence level. The d(S) and d(N) values are estimated based on the comparison of homologous codons, which are often identified by converting (reverse-translating) aligned amino acid sequences into codon sequences. In this method, however, homologous codons may be mis-identified when frame-shifts occurred or amino acid sequences were mis-aligned, which may lead to overestimation of the d(N)/d(S) ratio. Here the effect of reverse-translating aligned amino acid sequences on the estimation of d(N)/d(S) ratio was examined through a large-scale analysis of protein-coding nucleotide sequences from vertebrate species. Apparently, 1-9% of codon sites that were identified as homologous with reverse-translation contained non-homologous codons, where the d(N)/d(S) ratio was unduly high. By correcting the d(N)/d(S) ratio for these codon sites, it was inferred that the ratio was 5-43% overestimated with reverse-translation. These results suggest that caution should be exerted in the study of natural selection using the d(N)/d(S) ratio by reverse-translating aligned amino acid sequences.

Universal vaccines: shifting to one for many or shooting too high too soon!

Switching from conventional strain-specific vaccines to multi-strain or multi-species universal vaccines is both justified and scientifically merited. Long-term cross-protective universal vaccines eliminate the need for repetitive short-term vaccination campaigns and short-notice vaccine redesign during impending epidemics. They also have the potential to be cost-effective, convenient, and amenable to stockpiling. Ongoing advances in genomics and reverse vaccinology along with the perceived ability of vaccines, if properly formulated, to induce cross-protective adaptive immunity and long-term T cell memory are at the heart of this trend. Consequently, the search for universal vaccines against influenza, HIV, and many other viral, bacterial, and fungal pathogens has intensified in recent years. Currently, several universal influenza vaccines are at different phases of clinical evaluation. That said, vaccine-related differential effectiveness, escape mutants, pathogen strain replacement, limited scope of cross-protective immunity, and diminished potential to reach optimal herd immunity thresholds present serious challenges to the concept and applicability of universal vaccines. Herein, the case for and the case against universal vaccines are investigated to realistically appreciate their prospects of success.

Phylodynamics and molecular evolution of influenza A virus nucleoprotein genes in Taiwan between 1979 and 2009

Background

Many studies concentrate on variation in the hemagglutinin glycoprotein (HA) because of its significance in host immune response, the evolution of this virus is even more complex when other genome segments are considered. Recently, it was found that cytotoxic T lymphocytes (CTL) play an important role in immunity against influenza and most CTL epitopes of human influenza viruses were remarkably conserved. The NP gene has evolved independently in human and avian hosts after 1918 flu pandemic and it has been assigned a putative role as a determinant of host range.

Methods and Findings

Phylodynamic patterns of the genes encoding nucleoprotein (NP) of influenza A viruses isolated from 1979–2009 were analyzed by applying the Bayesian Markov Chain Monte Carlo framework to better understand the evolutionary mechanisms of these Taiwanese isolates. Phylogenetic analysis of the NP gene showed that all available H3 worldwide isolates collected so far were genetically similar and divided into two major clades after the year 2004. We compared the deduced amino acid sequences of the NP sequences from human, avian and swine hosts to investigate the emergence of potential adaptive mutations. Overall, selective pressure on the NP gene of human influenza A viruses appeared to be dominated by purifying selection with a mean d_N/d_S ratio of 0.105. Site-selection analysis of 488 codons, however, also revealed 3 positively selected sites in addition to 139 negatively selected ones.

Conclusions

The demographic history inferred by Bayesian skyline plot showed that the effective number of infections underwent a period of smooth and steady growth from 1998 to 2001, followed by a more recent rise in the rate of spread. Further understanding the correlates of interspecies transmission of influenza A virus genes from other host reservoirs to the human population may help to elucidate the mechanisms of variability among influenza A virus.

Identification of amino acids in highly pathogenic avian influenza H5N1 virus hemagglutinin that determine avian influenza species specificity

To test the role of neutralizing antibodies (nAbs) and receptor adaptation in interspecies transmission of influenza virus, two H5N1 strains, isolated from human and avian hosts, with four amino acid differences in hemagglutinin (HA) and seven HA mutations were studied. We found that a mutation at amino acid position 90 in the H5N1 HA, outside the receptor-binding domain (RBD), could simultaneously induce changes in the RBD conformation to escape from nAb binding and alter the receptor preference through long-range regulation. This mutation was deemed a “key event” for interspecies transmission. It is likely a result of positive selection caused by antibodies, allowing the original invasion by new species-specific variants. A mutation at amino acid position 160 in the RBD only induced a change in receptor preference. This mutation was deemed a “maintaining adaptation”, which ensured that influenza virus variants would be able to infect new organisms of a different species successfully. The mutation is the result of adaptation caused by the receptor. Our results suggest that continuing occurrence of these two types of mutations made the variants persist in the new host species.

Evolutionary dynamics of influenza A nucleoprotein (NP) lineages revealed by large-scale sequence analyses

Influenza A viral nucleoprotein (NP) plays a critical role in virus replication and host adaptation, however, the underlying molecular evolutionary dynamics of NP lineages are less well-understood. In this study, large-scale analyses of 5094 NP nucleotide sequences revealed eight distinct evolutionary lineages, including three host-specific lineages (human, classical swine and equine), two cross-host lineages (Eurasian avian-like swine and swine-origin human pandemic H1N1 2009) and three geographically isolated avian lineages (Eurasian, North American and Oceanian). The average nucleotide substitution rate of the NP lineages was estimated to be 2.4 × 10(-3) substitutions per site per year, with the highest value observed in pandemic H1N1 2009 (3.4 × 10(-3)) and the lowest in equine (0.9 × 10(-3)). The estimated time of most recent common ancestor (TMRCA) for each lineage demonstrated that the earliest human lineage was derived around 1906, and the latest pandemic H1N1 2009 lineage dated back to December 17, 2008. A marked time gap was found between the times when the viruses emerged and were first sampled, suggesting the crucial role for long-term surveillance of newly emerging viruses. The selection analyses showed that human lineage had six positive selection sites, whereas pandemic H1N1 2009, classical swine, Eurasian avian and Eurasian swine had only one or two sites. Protein structure analyses revealed several positive selection sites located in epitope regions or host adaptation regions, indicating strong adaptation to host immune system pressures in influenza viruses. Along with previous studies, this study provides new insights into the evolutionary dynamics of influenza A NP lineages. Further lineage analyses of other gene segments will allow better understanding of influenza A virus evolution and assist in the improvement of global influenza surveillance.

Statistical methods for detecting natural selection from genomic data

In the study of molecular and phenotypic evolution, understanding the relative importance of random genetic drift and positive selection as the mechanisms for driving divergences between populations and maintaining polymorphisms within populations has been a central issue. A variety of statistical methods has been developed for detecting natural selection operating at the amino acid and nucleotide sequence levels. These methods may be largely classified into those aimed at detecting recurrent and/or recent/ongoing natural selection by utilizing the divergence and/or polymorphism data. Using these methods, pervasive positive selection has been identified for protein-coding and non-coding sequences in the genomic analysis of some organisms. However, many of these methods have been criticized by using computer simulation and real data analysis to produce excessive false-positives and to be sensitive to various disturbing factors. Importantly, some of these methods have been invalidated experimentally. These facts indicate that many of the statistical methods for detecting natural selection are unreliable. In addition, the signals that have been believed as the evidence for fixations of advantageous mutations due to positive selection may also be interpreted as the evidence for fixations of deleterious mutations due to random genetic drift. The genomic diversity data are rapidly accumulating in various organisms, and detection of natural selection may play a critical role for clarifying the relative role of random genetic drift and positive selection in molecular and phenotypic evolution. It is therefore important to develop reliable statistical methods that are unbiased as well as robust against various disturbing factors, for inferring natural selection.

Variation in the analysis of positively selected sites using nonsynonymous/synonymous rate ratios: an example using influenza virus

Sites in a gene showing the nonsynonymous/synonymous rate ratio (ω) >1 have been frequently identified to be under positive selection. To examine the performance of such analysis, sites of the ω ratio >1 in the HA1 gene of H3N2 subtype human influenza viruses were identified from seven overlapping sequence data sets in this study. Our results showed that the sites of the ω ratio >1 were of significant variation among the data sets even though they targeted similar clusters, indicating that the analysis is likely to be either of low sensitivity or of low specificity in identifying sites under positive selection. Most (43/45) of the sites showing ω >1 calculated from at least one data set are involved in B-cell epitopes which cover less than a half sites in the protein, suggesting that the analysis is likely to be of low sensitivity rather than of low specificity. It was further found that the analysis sensitivity could not be enhanced by including more sequences or covering longer time intervals. Previously some reports also likely identified only a portion of the sites under positive selection in the viral gene using the ω ratio. Low sensitivity of the analysis may result from that some sites under positive selection in the gene are also under negative (purifying) selection simultaneously for functional constrains, and so their ω ratios could be <1. Theoretically, the sites under the two opposite selection forces at the same time favor only certain nonsynonymous changes, e.g. those changing the antigenicity of the gene and maintaining the gene function. This study also suggested that sometimes we can identify more sites under positive selection using the ω ratio by integrating the positively selected sites estimated from multiple data sets.

Molecular epidemiology and antigenic analyses of influenza A viruses H3N2 in Taiwan

The severity of an influenza epidemic season may be influenced not only by variability in the surface glycoproteins, but also by differences in the internal proteins of circulating influenza viruses. To better understand viral antigenic evolution, all eight gene segments from 44 human H3N2 epidemic strains isolated during 2004-2008 in Taiwan were analyzed to provide a profile of protein variability. Comparison of the evolutionary profiles of the HA, NA and PB2 genes of influenza A (H3N2) viruses indicated that they were derived from a group of H3N2 isolates first seen in 2004. However, the PA, M and PB1 genes were derived from a different group of H3N2 isolates from 2004. Tree topology revealed the NP and NS genes could each be segregated into two groups similar to those for the polymerase genes. In addition, new genetic variants occurred during the non-epidemic period and become the dominant strain after one or two seasons. Comparison of evolutionary patterns in consecutive years is necessary to correlate viral genetic changes with antigenic changes as multiple lineages co-circulate.

Evolution of the influenza A virus untranslated regions

In case of the influenza virus, untranslated regions are considered to possess plural functions. The specific packaging model suggests that each of eight segments contains a unique "packaging signal", situated within the untranslated region and coding regions at the 5' and 3' ends of the genomic RNA. In addition to packaging signal, UTRs play key roles in polyadenylation and signals for genomic replication. In the present study, we report the evolutionary characteristics of untranslated regions of influenza virus. We analyzed 574 sequences of the human influenza A (H3N2) virus. The sequence of the untranslated region is highly conserved. Our analysis produced several observations regarding the untranslated region: (1) it has an extremely low divergence; (2) its evolutionary speed is not associated with the protein(s) it encodes; (3) its evolutionary rate is smaller than that of the open reading frame. We revealed the evolutionary characteristics of the untranslated region; this is the first study on the non-coding region of the influenza virus from the perspective of evolutionary molecular biology. In the future, clarifying the detailed mechanisms of packaging as well as the function of the RNA sequence in the untranslated region will increase the understanding of viral life cycles and evolution.

The genomic rate of molecular adaptation of the human influenza A virus

Quantifying adaptive evolution at the genomic scale is an essential yet challenging aspect of evolutionary biology. Here, we develop a method that extends and generalizes previous approaches to estimate the rate of genomic adaptation in rapidly evolving populations and apply it to a large data set of complete human influenza A virus genome sequences. In accord with previous studies, we observe particularly high rates of adaptive evolution in domain 1 of the viral hemagglutinin (HA1). However, our novel approach also reveals previously unseen adaptation in other viral genes. Notably, we find that the rate of adaptation (per codon per year) is higher in surface residues of the viral neuraminidase than in HA1, indicating strong antibody-mediated selection on the former. We also observed high rates of adaptive evolution in several nonstructural proteins, which may relate to viral evasion of T-cell and innate immune responses. Furthermore, our analysis provides strong quantitative support for the hypothesis that human H1N1 influenza experiences weaker antigenic selection than H3N2. As well as shedding new light on the dynamics and determinants of positive Darwinian selection in influenza viruses, the approach introduced here is applicable to other pathogens for which densely sampled genome sequences are available, and hence is ideally suited to the interpretation of next-generation genome sequencing data.

Prevalence of epistasis in the evolution of influenza A surface proteins

The surface proteins of human influenza A viruses experience positive selection to escape both human immunity and, more recently, antiviral drug treatments. In bacteria and viruses, immune-escape and drug-resistant phenotypes often appear through a combination of several mutations that have epistatic effects on pathogen fitness. However, the extent and structure of epistasis in influenza viral proteins have not been systematically investigated. Here, we develop a novel statistical method to detect positive epistasis between pairs of sites in a protein, based on the observed temporal patterns of sequence evolution. The method rests on the simple idea that a substitution at one site should rapidly follow a substitution at another site if the sites are positively epistatic. We apply this method to the surface proteins hemagglutinin and neuraminidase of influenza A virus subtypes H3N2 and H1N1. Compared to a non-epistatic null distribution, we detect substantial amounts of epistasis and determine the identities of putatively epistatic pairs of sites. In particular, using sequence data alone, our method identifies epistatic interactions between specific sites in neuraminidase that have recently been demonstrated, in vitro, to confer resistance to the drug oseltamivir; these epistatic interactions are responsible for widespread drug resistance among H1N1 viruses circulating today. This experimental validation demonstrates the predictive power of our method to identify epistatic sites of importance for viral adaptation and public health. We conclude that epistasis plays a large role in shaping the molecular evolution of influenza viruses. In particular, sites with , which would normally not be identified as positively selected, can facilitate viral adaptation through epistatic interactions with their partner sites. The knowledge of specific interactions among sites in influenza proteins may help us to predict the course of antigenic evolution and, consequently, to select more appropriate vaccines and drugs.

Epistasis describes non-additive interactions among genetic sites: the consequence of a mutation at one site may depend on the status of the genome at other sites. In an extreme case, a mutation may have no effect if it arises on one genetic background, but a strong effect on another background. Epistatic mutations in viruses and bacteria that live under severe conditions, such as antibiotic treatments or immune pressure, often allow pathogens to develop drug resistance or escape the immune system. In this paper we develop a new phylogenetic method for detecting epistasis, and we apply this method to the surface proteins of the influenza A virus, which are important targets of the immune system and drug treatments. The authors identify and characterize hundreds of epistatic mutations in these proteins. Among those identified, we find the specific epistatic mutations that were recently shown, experimentally, to confer resistance to the drug Tamiflu. The results of this study may help to predict the course of influenza's antigenic evolution and to select more appropriate vaccines and drugs.

Rapid induction of IgE responses to a worm cysteine protease during murine pre-patent schistosome infection

BACKGROUND

During the pre-patent stage of infection, juvenile Schistosoma blood flukes co-opt signals from the adaptive immune system to facilitate parasite development, but the types of responses that are induced at this early stage of infection, and the parasite antigens they target, have not been characterized.

RESULTS

Through analysis of experimental pre-patent infections, we show that the S. mansoni cysteine protease SmCB1 is rapidly targeted by an antigen-specific IgE response. The induction of this response is independent of schistosome eggs as infection with male or female worms alone also induced SmCB1-specific IgE. We also show that the SmCB1-specific IgE response is dependent on cognate CD4+ T cell help and IL-4, suggesting that pre-patent Th2 responses provide T cell help for the SmCB1-specific IgE response. Finally, exposed human subjects also produced IgE against SmCB1.

CONCLUSIONS

Our data demonstrate that, like eggs, schistosome worms also induce functional type 2 responses and that a parasite cysteine protease is an inducer of type 2 responses during the early stages of schistosome infection.

Influenza A virus protein PB1-F2 exacerbates IFN-beta expression of human respiratory epithelial cells

The PB1-F2 protein of the influenza A virus (IAV) contributes to viral pathogenesis by a mechanism that is not well understood. PB1-F2 was shown to modulate apoptosis and to be targeted by the CD8(+) T cell response. In this study, we examined the downstream effects of PB1-F2 protein during IAV infection by measuring expression of the cellular genes in response to infection with wild-type WSN/33 and PB1-F2 knockout viruses in human lung epithelial cells. Wild-type virus infection resulted in a significant induction of genes involved in innate immunity. Knocking out the PB1-F2 gene strongly decreased the magnitude of expression of cellular genes implicated in antiviral response and MHC class I Ag presentation, suggesting that PB1-F2 exacerbates innate immune response. Biological network analysis revealed the IFN pathway as a link between PB1-F2 and deregulated genes. Using quantitative RT-PCR and IFN-β gene reporter assay, we determined that PB1-F2 mediates an upregulation of IFN-β expression that is dependent on NF-κB but not on AP-1 and IFN regulatory factor-3 transcription factors. Recombinant viruses knocked out for the PB1-F2 and/or the nonstructural viral protein 1 (the viral antagonist of the IFN response) genes provide further evidence that PB1-F2 increases IFN-β expression and that nonstructural viral protein 1 strongly antagonizes the effect of PB1-F2 on the innate response. Finally, we compared the effect of PB1-F2 variants taken from several IAV strains on IFN-β expression and found that PB1-F2-mediated IFN-β induction is significantly influenced by its amino acid sequence, demonstrating its importance in the host cell response triggered by IAV infection.

Genetic characterization of circulating seasonal Influenza A viruses (2005-2009) revealed introduction of oseltamivir resistant H1N1 strains during 2009 in eastern India

Influenza surveillance was implemented in Kolkata, eastern India in 2005 to identify the circulating subtypes and characterize their genetic diversity. Throat and nasal swabs were collected from outpatients with influenza-like illness (ILI). Of 2844 ILI cases identified at two referral hospitals during October 2005-September 2009, 309 (10.86%) were positive for Influenza A by real time RT-PCR, of which 110 (35.60%) were subtyped as H1N1 and 199 (64.40%) as H3N2. Comparison of the nucleotide (nt) and amino acid (aa) sequences of the HA1 gene for H1N1 and H3N2 strains showed that a subset of strains precede WHO recommended contemporary strains by 1-2 years. The Kolkata H1N1 strains clustered in Clade II, subgroup 2B with A/Brisbane/59/2007 but were distant from the corresponding vaccine strains (New Caledonia/20/99 and A/Solomon Island/3/06). The 2005-06 and 2007 H3N2 strains (15/17) clustered either A/Brisbane/10/2007-like (n=8) or A/Nepal/921/2006 like (n=7) strains, whereas 2008 strains (8/12) and 2009 strains (4/4) were similar to the 2010-11 vaccine strain A/Perth/16/2009. More aa substitutions were found in HA or NA genes of H3N2 than in H1N1 strains. No mutation conferring neuraminidase resistance was observed in any of the strain during 2005-08, however in 2009, drug resistant marker (H275Y) was present in seasonal H1N1, but not in co-circulating H3N2 strains. This is the first report of genetic characterization of circulating Influenza A strains from India. The results also highlight the importance of continuing Influenza surveillance in developing countries of Asia for monitoring unusual strains with pandemic potential and mutations conferring antiviral resistance.

Differential localization and function of PB1-F2 derived from different strains of influenza A virus

PB1-F2 is a viral protein that is encoded by the PB1 gene of influenza A virus by alternative translation. It varies in length and sequence context among different strains. The present study examines the functions of PB1-F2 proteins derived from various human and avian viruses. While H1N1 PB1-F2 was found to target mitochondria and enhance apoptosis, H5N1 PB1-F2, surprisingly, did not localize specifically to mitochondria and displayed no ability to enhance apoptosis. Introducing Leu into positions 69 (Q69L) and 75 (H75L) in the C terminus of H5N1 PB1-F2 drove 40.7% of the protein to localize to mitochondria compared with the level of mitochondrial localization of wild-type H5N1 PB1-F2, suggesting that a Leu-rich sequence in the C terminus is important for targeting of mitochondria. However, H5N1 PB1-F2 contributes to viral RNP activity, which is responsible for viral RNA replication. Lastly, although the swine-origin influenza virus (S-OIV) contained a truncated form of PB1-F2 (12 amino acids [aa]), potential mutation in the future may enable it to contain a full-length product. Therefore, the functions of this putative S-OIV PB1-F2 (87 aa) were also investigated. Although this PB1-F2 from the mutated S-OIV shares only 54% amino acid sequence identity with that of seasonal H1N1 virus, it also increased viral RNP activity. The plaque size and growth curve of the viruses with and without S-OIV PB1-F2 differed greatly. The PB1-F2 protein has various lengths, amino acid sequences, cellular localizations, and functions in different strains, which result in strain-specific pathogenicity. Such genetic and functional diversities make it flexible and adaptable in maintaining the optimal replication efficiency and virulence for various strains of influenza A virus.

Comparison of selection pressures on the HA gene of pandemic (2009) and seasonal human and swine influenza A H1 subtype viruses

The 2009 human pandemic influenza (H1N1) virus possesses the HA gene of the H1 subtype. The evolutionary process of the 2009 H1N1 virus remains to be defined. We performed genetic analyses of the HA gene by comparing the 2009 H1N1 virus with seasonal human and swine viruses. We analyzed sequences of 116 2009 H1N1 viruses, and obtained 1457 seasonal H1N1, 365 swine H1, and 1332 2009 H1N1 viruses from the database. Selection pressure for the 2009 H1N1 virus was higher than that for the swine virus and equivalent to that for the seasonal virus. Positions 206 and 264 were found to be positively selected sites. We also identified sites under different selection pressures from the seasonal or swine virus that may be involved in imparting significant biological characteristics. The evolutionary characteristics of the H1 gene of the 2009 H1N1 virus differed from those of seasonal and swine viruses.

Use of phylogenetics in the molecular epidemiology and evolutionary studies of viral infections

Since DNA sequencing techniques first became available almost 30 years ago, the amount of nucleic acid sequence data has increased enormously. Phylogenetics, which is widely applied to compare and analyze such data, is particularly useful for the analysis of genes from rapidly evolving viruses. It has been used extensively to describe the molecular epidemiology and transmission of the human immunodeficiency virus (HIV), the origins and subsequent evolution of the severe acute respiratory syndrome (SARS)-associated coronavirus (SCoV), and, more recently, the evolving epidemiology of avian influenza as well as seasonal and pandemic human influenza viruses. Recent advances in phylogenetic methods can infer more in-depth information about the patterns of virus emergence, adding to the conventional approaches in viral epidemiology. Examples of this information include estimations (with confidence limits) of the actual time of the origin of a new viral strain or its emergence in a new species, viral recombination and reassortment events, the rate of population size change in a viral epidemic, and how the virus spreads and evolves within a specific population and geographical region. Such sequence-derived information obtained from the phylogenetic tree can assist in the design and implementation of public health and therapeutic interventions. However, application of many of these advanced phylogenetic methods are currently limited to specialized phylogeneticists and statisticians, mainly because of their mathematical basis and their dependence on the use of a large number of computer programs. This review attempts to bridge this gap by presenting conceptual, technical, and practical aspects of applying phylogenetic methods in studies of influenza, HIV, and SCoV. It aims to provide, with minimal mathematics and statistics, a practical overview of how phylogenetic methods can be incorporated into virological studies by clinical and laboratory specialists.

Differentiation of human influenza A viruses including the pandemic subtype H1N1/2009 by conventional multiplex PCR

April 2009 witnessed the emergence of a novel H1N1 influenza A virus infecting the human population. Currently, pandemic and seasonal influenza viruses are co-circulating in human populations. Understanding the course of the emerging pandemic virus is important. It is still unknown how the novel virus co-circulates with or outcompetes seasonal viruses. Sustainable and detailed influenza surveillance is required throughout the world including developing countries. In the present study, a multiplex PCR using four primers was developed, which was designed to differentiate the pandemic H1N1 virus from the seasonal H1N1 and H3N2 viruses, to obtain amplicons of different sizes. Multiplex PCR analysis could clearly differentiate the three subtypes of human influenza A virus. This assay was performed using 206 clinical samples collected in 2009 in Japan. Between February and April, four samples were subtyped as seasonal H1N1 and four as seasonal H3N2. All samples collected after July were subtyped as pandemic H1N1. Currently, pandemic viruses seem to have replaced seasonal viruses almost completely in Japan. This is a highly sensitive method and its cost is low. Influenza surveillance using this assay would provide significant information on the epidemiology of both pandemic and seasonal influenza.

Reassortment between swine influenza A viruses increased their adaptation to humans in pandemic H1N1/09

In April 2009, pandemic H1N1/09 influenza, which originated from swine influenza, appeared in North America, and it has since spread globally among humans. It is important to know how swine influenza A virus broke the host barrier to cause a pandemic. We analyzed 673 strains of human, avian, and swine influenza viruses and assessed the internal genes PB2, PB1, PA, NP, M, and NS. Here we found accumulation of mutations in segments that were retained as well as introduced due to genetic reassortment of viruses. The retained segments may have to mutate to accommodate new segments. The mutations caused by interaction among segments retained and introduced due to reassortment between swine influenza viruses may have increased the adaptation of the virus to humans, leading to pandemic H1N1/09. We indicate the sites that probably contributed to the acquisition of efficient human-to-human transmission.

Hitchhiking and the population genetic structure of avian influenza virus

Previous studies have revealed a major difference in the phylogenetic structure, extent of genetic diversity, and selection pressure between the surface glycoproteins and internal gene segments of avian influenza viruses (AIV) sampled from wild birds. However, what evolutionary processes are responsible for these strikingly different evolutionary patterns is unclear. To address this issue, we estimated the rate of evolutionary change and time of origin of each segment of AIV sampled globally. Strikingly, the internal segments of the sampled AIV strains possess common ancestors that existed less than 200 years ago. Similarly recent times of origin were observed for each of the individual subtypes within the HA, NA, and NS gene segments. Such a shallow history of genetic diversity suggests an evolutionary model in which the genetic structure of AIV is shaped by a combination of occasional selective sweeps in the HA and NA (and possibly NS) segments, coupled with transient genetic linkage to the internal gene segments.

Evidence for a novel gene associated with human influenza A viruses

BACKGROUND

Influenza A virus genomes are comprised of 8 negative strand single-stranded RNA segments and are thought to encode 11 proteins, which are all translated from mRNAs complementary to the genomic strands. Although human, swine and avian influenza A viruses are very similar, cross-species infections are usually limited. However, antigenic differences are considerable and when viruses become established in a different host or if novel viruses are created by re-assortment devastating pandemics may arise.

RESULTS

Examination of influenza A virus genomes from the early 20th Century revealed the association of a 167 codon ORF encoded by the genomic strand of segment 8 with human isolates. Close to the timing of the 1948 pseudopandemic, a mutation occurred that resulted in the extension of this ORF to 216 codons. Since 1948, this ORF has been almost totally maintained in human influenza A viruses suggesting a selectable biological function. The discovery of cytotoxic T cells responding to an epitope encoded by this ORF suggests that it is translated into protein. Evidence of several other non-traditionally translated polypeptides in influenza A virus support the translation of this genomic strand ORF. The gene product is predicted to have a signal sequence and two transmembrane domains.

CONCLUSION

We hypothesize that the genomic strand of segment 8 of encodes a novel influenza A virus protein. The persistence and conservation of this genomic strand ORF for almost a century in human influenza A viruses provides strong evidence that it is translated into a polypeptide that enhances viral fitness in the human host. This has important consequences for the interpretation of experiments that utilize mutations in the NS1 and NEP genes of segment 8 and also for the consideration of events that may alter the spread and/or pathogenesis of swine and avian influenza A viruses in the human population.

Large-scale sequence analysis of M gene of influenza A viruses from different species: mechanisms for emergence and spread of amantadine resistance

Influenza A virus infects many species, and amantadine is used as an antiviral agent. Recently, a substantial increase in amantadine-resistant strains has been reported, most of which have a substitution at amino acid position 31 in the M2 gene. Understanding the mechanism responsible for the emergence and spread of antiviral resistance is important for developing a treatment protocol for seasonal influenza and for deciding on a policy for antiviral stockpiling for pandemic influenza. The present study was conducted to identify the existence of drug pressure on the emergence and spread of amantadine-resistant influenza A viruses. We analyzed data on more than 5,000 virus sequences and constructed a phylogenetic tree to calculate selective pressures on sites in the M2 gene associated with amantadine resistance (positions 26, 27, 30, and 31) among different hosts. The phylogenetic tree revealed that the emergence and spread of the drug-resistant M gene in different hosts and subtypes were independent and not through reassortment. For human influenza virus, positive selection was detected only at position 27. Selective pressures on the sites were not always higher for human influenza virus than for viruses of other hosts. Additionally, selective pressure on position 31 did not increase after the introduction of amantadine. Although there is a possibility of drug pressure on human influenza virus, we could not find positive pressure on position 31. Because the recent rapid increase in drug-resistant virus is associated with the substitution at position 31, the resistance may not be related to drug use.

Evolution of the M gene of the influenza A virus in different host species: large-scale sequence analysis

Background

Influenza A virus infects not only humans, but also other species including avian and swine. If a novel influenza A subtype acquires the ability to spread between humans efficiently, it could cause the next pandemic. Therefore it is necessary to understand the evolutionary processes of influenza A viruses in various hosts in order to gain better knowledge about the emergence of pandemic virus. The virus has segmented RNA genome and 7th segment, M gene, encodes 2 proteins. M1 is a matrix protein and M2 is a membrane protein. The M gene may be involved in determining host tropism. Besides, novel vaccines targeting M1 or M2 protein to confer cross subtype protection have been under development. We conducted the present study to investigate the evolution of the M gene by analyzing its sequence in different species.

Results

Phylogenetic tree revealed host-specific lineages and evolution rates were different among species. Selective pressure on M2 was stronger than that on M1. Selective pressure on M1 for human influenza was stronger than that for avian influenza, as well as M2. Site-by-site analyses identified one site (amino acid position 219) in M1 as positively selected in human. Positions 115 and 121 in M1, at which consensus amino acids were different between human and avian, were under negative selection in both hosts. As to M2, 10 sites were under positive selection in human. Seven sites locate in extracellular domain. That might be due to host's immune pressure. One site (position 27) positively selected in transmembrane domain is known to be associated with drug resistance. And, two sites (positions 57 and 89) locate in cytoplasmic domain. The sites are involved in several functions.

Conclusion

The M gene of influenza A virus has evolved independently, under different selective pressure on M1 and M2 among different hosts. We found potentially important sites that may be related to host tropism and immune responses. These sites may be important for evolutional process in different hosts and host adaptation.

Interspecies transmission and host restriction of avian H5N1 influenza virus

Long-term endemicity of avian H5N1 influenza virus in poultry and continuous sporadic human infections in several countries has raised the concern of another potential pandemic influenza. Suspicion of the avian origin of the previous pandemics results in the close investigation of the mechanism of interspecies transmission. Entry and fusion is the first step for the H5N1 influenza virus to get into the host cells affecting the host ranges. Therefore receptor usage study has been a major focus for the last few years. We now know the difference of the sialic acid structures and distributions in different species, even in the different parts of the same host. Many host factors interacting with the influenza virus component proteins have been identified and their role in the host range expansion and interspecies transmission is under detailed scrutiny. Here we review current progress in the receptor usage and host factors.

A method for the simultaneous estimation of selection intensities in overlapping genes

Inferring the intensity of positive selection in protein-coding genes is important since it is used to shed light on the process of adaptation. Recently, it has been reported that overlapping genes, which are ubiquitous in all domains of life, seem to exhibit inordinate degrees of positive selection. Here, we present a new method for the simultaneous estimation of selection intensities in overlapping genes. We show that the appearance of positive selection is caused by assuming that selection operates independently on each gene in an overlapping pair, thereby ignoring the unique evolutionary constraints on overlapping coding regions. Our method uses an exact evolutionary model, thereby voiding the need for approximation or intensive computation. We test the method by simulating the evolution of overlapping genes of different types as well as under diverse evolutionary scenarios. Our results indicate that the independent estimation approach leads to the false appearance of positive selection even though the gene is in reality subject to negative selection. Finally, we use our method to estimate selection in two influenza A genes for which positive selection was previously inferred. We find no evidence for positive selection in both cases.

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.

Evolution of drug resistance in multiple distinct lineages of H5N1 avian influenza

Some predict that influenza A H5N1 will be the cause of a pandemic among humans. In preparation for such an event, many governments and organizations have stockpiled antiviral drugs such as oseltamivir (Tamiflu). However, it is known that multiple lineages of H5N1 are already resistant to another class of drugs, adamantane derivatives, and a few lineages are resistant to oseltamivir. What is less well understood is the evolutionary history of the mutations that confer drug resistance in the H5N1 population. In order to address this gap, we conducted phylogenetic analyses of 676 genomic sequences of H5N1 and used the resulting hypotheses as a basis for asking 3 molecular evolutionary questions: (1) Have drug-resistant genotypes arisen in distinct lineages of H5N1 through point mutation or through reassortment? (2) Is there evidence for positive selection on the codons that lead to drug resistance? (3) Is there evidence for covariation between positions in the genome that confer resistance to drugs and other positions, unrelated to drug resistance, that may be under selection for other phenotypes? We also examine how drug-resistant lineages proliferate across the landscape by projecting or phylogenetic analysis onto a virtual globe. Our results for H5N1 show that in most cases drug resistance has arisen by independent point mutations rather than reassortment or covariation. Furthermore, we found that some codons that mediate resistance to adamantane derivatives are under positive selection, but did not find positive selection on codons that mediate resistance to oseltamivir. Together, our phylogenetic methods, molecular evolutionary analyses, and geographic visualization provide a framework for analysis of globally distributed genomic data that can be used to monitor the evolution of drug resistance.

A maximum likelihood method for detecting directional evolution in protein sequences and its application to influenza A virus

We develop a model-based phylogenetic maximum likelihood test for evidence of preferential substitution toward a given residue at individual positions of a protein alignment--directional evolution of protein sequences (DEPS). DEPS can identify both the target residue and sites evolving toward it, help detect selective sweeps and frequency-dependent selection--scenarios that confound most existing tests for selection, and achieve good power and accuracy on simulated data. We applied DEPS to alignments representing different genomic regions of influenza A virus (IAV), sampled from avian hosts (H5N1 serotype) and human hosts (H3N2 serotype), and identified multiple directionally evolving sites in 5/8 genomic segments of H5N1 and H3N2 IAV. We propose a simple descriptive classification of directionally evolving sites into 5 groups based on the temporal distribution of residue frequencies and document known functional correlates, such as immune escape or host adaptation.

The evolution of human influenza A viruses from 1999 to 2006: a complete genome study

Background

Knowledge about the complete genome constellation of seasonal influenza A viruses from different countries is valuable for monitoring and understanding of the evolution and migration of strains. Few complete genome sequences of influenza A viruses from Europe are publicly available at the present time and there have been few longitudinal genome studies of human influenza A viruses. We have studied the evolution of circulating human H3N2, H1N1 and H1N2 influenza A viruses from 1999 to 2006, we analysed 234 Danish human influenza A viruses and characterised 24 complete genomes.

Results

H3N2 was the prevalent strain in Denmark during the study period, but H1N1 dominated the 2000–2001 season. H1N2 viruses were first observed in Denmark in 2002–2003. After years of little genetic change in the H1N1 viruses the 2005–2006 season presented H1N1 of greater variability than before. This indicates that H1N1 viruses are evolving and that H1N1 soon is likely to be the prevalent strain again. Generally, the influenza A haemagglutinin (HA) of H3N2 viruses formed seasonal phylogenetic clusters. Different lineages co-circulating within the same season were also observed. The evolution has been stochastic, influenced by small "jumps" in genetic distance rather than constant drift, especially with the introduction of the Fujian-like viruses in 2002–2003. Also evolutionary stasis-periods were observed which might indicate well fit viruses. The evolution of H3N2 viruses have also been influenced by gene reassortments between lineages from different seasons. None of the influenza genes were influenced by strong positive selection pressure. The antigenic site B in H3N2 HA was the preferred site for genetic change during the study period probably because the site A has been masked by glycosylations. Substitutions at CTL-epitopes in the genes coding for the neuraminidase (NA), polymerase acidic protein (PA), matrix protein 1 (M1), non-structural protein 1 (NS1) and especially the nucleoprotein (NP) were observed. The N-linked glycosylation pattern varied during the study period and the H3N2 isolates from 2004 to 2006 were highly glycosylated with ten predicted sequons in HA, the highest amount of glycosylations observed in this study period.

Conclusion

The present study is the first to our knowledge to characterise the evolution of complete genomes of influenza A H3N2, H1N1 and H1N2 isolates from Europe over a time period of seven years from 1999 to 2006. More precise knowledge about the circulating strains may have implications for predicting the following season strains and thereby better matching the vaccine composition.

Natural selection for nucleotide usage at synonymous and nonsynonymous sites in influenza A virus genes

Influenza A virus is one of the best-studied viruses and a model organism for the study of molecular evolution; in particular, much research has focused on detecting natural selection on influenza virus proteins. Here, we study the dynamics of the synonymous and nonsynonymous nucleotide composition of influenza A virus genes. In several genes, the nucleotide frequencies at synonymous positions drift away from the equilibria predicted from the synonymous substitution matrices. We investigate possible reasons for this unexpected behavior by fitting several regression models. Relaxation toward a mutation-selection equilibrium following a host jump fails to explain the dynamics of the synonymous nucleotide composition, even if we allow for slow temporal changes in the substitution matrix. Instead, we find that deep internal branches of the phylogeny show distinct patterns of nucleotide substitution and that these branches strongly influence the dynamics of nucleotide composition, suggesting that the observed trends are at least in part a result of natural selection acting on synonymous sites. Moreover, we find that the dynamics of the nucleotide composition at synonymous and nonsynonymous sites are highly correlated, providing evidence that even nonsynonymous sites can be influenced by selection pressure for nucleotide composition.

Approaches for identifying targets of positive selection

Despite significant advancements in both empirical and theoretical population genetics throughout the past century, fundamental questions about the evolutionary forces that shape genomic diversity remain unresolved. Perhaps foremost among these are the strength and frequency of adaptive evolution. To quantify these parameters, statistical tools are needed that are capable of effectively identifying targets of positive selection throughout the genome in an unbiased manner, and functional approaches are needed that are capable of connecting these identified genotypes with the resulting adaptively significant phenotypes. Here we review recent advancements in both statistical and empirical methodology, and discuss important challenges and opportunities that remain as researchers continue to uncouple the relative importance of stochastic and deterministic factors in the evolution of natural populations.

More effective purifying selection on RNA viruses than in DNA viruses

Analysis of the pattern of nucleotide diversity in 222 independent viral sequence datasets showed the prevalence of purifying selection. In spite of the higher mutation rate of RNA viruses, our analyses revealed stronger evidence of the action of purifying selection in RNA viruses than in DNA viruses. The ratio of nonsynonymous to synonymous nucleotide diversity was significantly lower in RNA viruses than in DNA viruses, indicating that nonsynonymous mutations have been removed at a greater rate (relative to the mutation rate) in the former than in the latter. Moreover, statistics that measure the occurrence of rare polymorphisms revealed significantly a greater excess of rare nonsynonymous polymorphisms in RNA viruses than in DNA viruses but no difference with respect to synonymous polymorphisms. Since rare nonsynonymous polymorphisms are likely to be undergoing the effects of purifying selection acting to eliminate them, this result implies a stronger signature of ongoing purifying selection in RNA viruses than in DNA viruses. Across datasets from both DNA viruses and RNA viruses, we found a negatively allometric relationship between nonsynonymous and synonymous nucleotide diversity; in other words, nonsynonymous nucleotide diversity increased with synonymous nucleotide diversity at a less than linear rate. These findings are most easily explained by the occurrence of slightly deleterious mutations. The fact that the negative allometry was more pronounced in RNA viruses than in DNA viruses provided additional evidence that purifying selection is more effective in the former than in the latter.

Pattern of nucleotide substitution in the overlapping nonstructural genes of influenza A virus and implication for the genetic diversity of the H5N1 subtype

In viruses under strong pressure to minimize genome size, overlapping genes represent a fine strategy to condense a maximum amount of information into short nucleotide sequences. Here, we investigated the evolution of the genes encoding the nonstructural proteins NS1 and NS2 of influenza A virus (IAV), which are one of the best characterized cases of gene overlap. By a detailed analysis of about four hundred sequences grouped into 11 IAV subtypes, we found that the overlapping coding region of the NS1 gene shows a significant increase of the rate of nonsynonymous change, with respect to its nonoverlapping counterpart. The same feature was observed in the overlapping coding region of the NS2 gene. Such a variation pattern, which implies the occurrence of several amino acid substitutions in the protein regions encoded by overlapping frames, is different from the pattern of constrained evolution typical of other viral overlapping-gene systems. Amino acid sequence analysis of the NS1 and NS2 proteins revealed that some nonsynonymous substitutions, located in the region of gene overlap, play a critical role in shaping the genetic diversity of the highly pathogenic subtype H5N1. Since both proteins contribute to disease pathogenesis by affecting many virus and host-cell processes, information provided by this study should be useful to highlight the impact of nonstructural gene variation on the pathogenicity of H5N1 viruses.

Computational analysis and identification of amino acid sites in dengue E proteins relevant to development of diagnostics and vaccines

We have identified 72 completely conserved amino acid residues in the E protein of major groups of the Flavivirus genus by computational analyses. In the dengue species we have identified 12 highly conserved sequence regions, 186 negatively selected sites, and many dengue serotype-specific negatively selected sites. The flavivirus-conserved sites included residues involved in forming six disulfide bonds crucial for the structural integrity of the protein, the fusion motif involved in viral infectivity, and the interface residues of the oligomers. The structural analysis of the E protein showed 19 surface-exposed non-conserved residues, 128 dimer or trimer interface residues, and regions, which undergo major conformational change during trimerization. Eleven consensus T(h)-cell epitopes common to all four dengue serotypes were predicted. Most of these corresponded to dengue-conserved regions or negatively selected sites. Of special interest are six singular sites (N(37), Q(211), D(215), P(217), H(244), K(246)) in dengue E protein that are conserved, are part of the predicted consensus T(h)-cell epitopes and are exposed in the dimer or trimer. We propose these sites and corresponding epitopic regions as potential candidates for prioritization by experimental biologists for development of diagnostics and vaccines that may be difficult to circumvent by natural or man-made alteration of dengue virus.

The evolution of epidemic influenza

Recent developments in complete-genome sequencing, antigenic mapping and epidemiological modelling are greatly improving our knowledge of the evolution of human influenza virus at the epidemiological scale. In particular, recent studies have revealed a more complex relationship between antigenic evolution, natural selection and reassortment than previously realized. Despite these advances, there is much that remains to be understood about the epidemiology of influenza virus, particularly the processes that determine the virus's strong seasonality. We argue that a complete understanding of the evolutionary biology of this important human pathogen will require a genomic view of genetic diversity, including the acquisition of polymorphism data from within individual hosts and from geographical regions, particularly the tropics, which have been poorly surveyed to date.