Skip to main content

Modeling gene sequences over time in 2009 H1N1 Influenza A Virus populations

Abstract

Background

A sudden emergence of Influenza A Virus (IAV) infections with a new pandemic H1N1 IAV is taking place since April of 2009. In order to gain insight into the mode of evolution of these new H1N1 strains, we performed a Bayesian coalescent Markov chain Monte Carlo (MCMC) analysis of full-length neuraminidase (NA) gene sequences of 62 H1N1 IAV strains (isolated from March 30th to by July 28th, 2009).

Results

The results of these studies revealed that the expansion population growth model was the best to fit the sequence data. A mean of evolutionary change of 7.84 × 10-3 nucleotide substitutions per site per year (s/s/y) was obtained for the NA gene. A significant contribution of first codon position to this mean rate was observed. Maximum clade credibility trees revealed a rapid diversification of NA genes in different genetic lineages, all of them containing Oseltamivir-resistant viruses of very recent emergence. Mapping of naturally occurring amino acid substitutions in the NA protein from 2009 H1N1 IAV circulating in 62 different patients revealed that substitutions are distributed all around the surface of the molecule, leaving the hydrophobic core and the catalytic site essentially untouched.

Conclusion

High evolutionary rates and fast population growth have contributed to the initial transmission dynamics of 2009 H1N1 IAV. Naturally occurring substitutions are preferentially located at the protein surface and do not interfere with the NA active site. Antigenic regions relevant for vaccine development can differ from previous vaccine strains and vary among patients.

Background

Influenza A virus (IAV) is a member of the family Orthomyxoviridae and contains eight segments of a single-stranded RNA genome with negative polarity [1]. IAV causes 300,000-500,000 deaths worldwide each year, and in pandemic years, this number can increase to 1 million (in 1957-1958) or as high as 50 million, as was seen in 1918-1919 [2]. Unlike most pathogens where exposure leads to lasting immunity in the host, IAV presents a moving antigenic target [3], evading specific immunity triggered by previous infections. This process, called antigenic drift, is the result of the selective fixation of mutations in the gene encoding the hemagglutinin (HA) protein, the major target for the host immune response [4]. Variants that best escape the host immune response are thought to have a significant reproductive advantage [5].

Another process, called antigenic shift, is also considered a major force in the evolution of IAV [4, 5]. Antigenic shift occurs when the virus acquires an HA of a different IAV subtype via reassortment of one or more gene segments and is thought to be the basis for the more devastating influenza pandemics that occurred several times in the last century [6]. New IAV pandemics may emerge through reassortation with strains from swine or avian reservoirs [7].

There have been three pandemics in the last hundred years: in 1918 (H1N1 subtype) [8], 1957 (H2N2 subtype) [9], and in 1968 (H3N2 subtype) [10]. During each of these pandemics, the new virus drove the previous pandemic subtype out of circulation [3]. In 1977, the H1N1 subtype reappeared, and has been co-circulating with H3N2 since then [11, 12].

IAV H3N2 viruses have been the predominant strains during the last 20 years, with the exception of the 1988-1989 and 2000-2001 seasons where H1N1 infections dominated [13].

A sudden emergence of IAV infections with new H1N1 strains of pandemic potential is taking place since April of 2009, starting in Mexico and spreading to several other countries around the world [14]. The World Health Organization (WHO) has raised the Influenza pandemic alert to the maximum level 6 [15].

Oseltamivir phosphate is a prodrug of oseltamivir carboxylate, a highly specific inhibitor of IAV neuraminidases. Oseltamivir carboxylate binds to highly conserved, essential amino acids in the catalytic site of neuraminidase (NA), preventing virus release from infected cells and subsequent virus spread [16]. An amino acid substitution at position 275 (H275Y) of the NA protein has been associated to resistance to Oseltamivir [17].

Initial testing of the 2009 pandemic H1N1 IAV strains found the viruses to be susceptible to neuraminidase inhibitors (oseltamivir and zanamivir).

Detailed studies on the mode of evolution of these new H1N1 IAV strains are extremely important for our understanding of the molecular mechanisms involved in the emergence, spread and resistance of new H1N1 IAV strains of pandemic potential. In order gain insight into these matters, we have performed a Bayesian coalescent Markov chain Monte Carlo analysis of full-length NA gene sequences of 62 emerging 2009 H1N1 IAV strains (isolated from March 30th to July 28th, 2009). The results of these studies revealed high rate of evolutionary change of NA genes, fast expansion of the H1N1 IAV populations and emergence of anti-viral resistant viruses. Naturally occurring amino acid substitutions in the NA of H1N1 IAV strains circulating in 62 different patients preferentially located at the protein surface and do not interfere with the NA active site.

Methods

Neuraminidase sequences

Full-length NA sequences from the 2009 emerging H1N1 IAV strains, were obtained from The Influenza Virus Resource at the National Center for Biotechnological Information [18]. For strain names, dates of isolation and accession numbers see Table S1, Additional file 1.

Sequence alignment

NA sequences were aligned using the MUSCLE program [19].

Evolutionary Model analysis

Once aligned, the FindModel program [20] was used to identify the optimal evolutionary model that best fitted our sequence dataset. Akaike Information Criteria revealed that the General Time Reversible (GTR) model was the best fit to the data (Table S2, Additional file 2).

Recombination Detection Tests

To test whether a recombination event occurred on any of the sequences included in these studies, two different approaches implemented in the SimPlot program [21] were used: (1) a sliding window analysis of distances and (2) the bootscanning [22]. No recombinant strains were found in the datasets (not shown).

Bayesian Coalescent Inference Studies

The evolutionary rate and mode of evolution of the newly emerging 2009 H1N1 IAV strains were determined using a coalescent Bayesian Markov chain Monte Carlo (MCMC) approach as implemented in the BEAST package [23]. Sixty-two full-length NA gene sequences were included in these analyses. For names, accession numbers and date of isolation of strains included in these studies, see Table S1, Additional file 1,. Using the GTR model, 60 million steps of MCMC and dates introduced by day of isolation, different population dynamic models were tested (constant population size, exponential population growth, expansion population growth, logistic population growth and Bayesian Skyline). Statistical uncertainty in the data was reflected by the 95% highest probability density (HPD) values. Results were examined using the TRACER program from the BEAST package [24]. Convergence was assessed with ESS (Effective Sample Size) values, after a burning of 6 million steps. Maximum clade credibility trees were generated using Tree Annotator from the BEAST package and the FigTree v1.2.2 (available at: http://tree.bio.ed.ac.uk/) was used for the visualization of the annotated trees.

Results

Modelling gene sequences changes over time in NA gene of 2009 H1N1 emerging strains

In order to gain insight into the evolutionary rate and mode of evolution of 2009 H1N1 IAV strains, we used a Bayesian Markov Chain Montecarlo (MCMC) approach to analyze 62 full-length NA gene sequences from 2009 H1N1 IAV strains isolated from March 30th to July 28th, 2009 (for strains names, accession numbers and dates of isolation, see Table S1, Additional file 1,).

Using the GTR model and 60 million steps of MCMC, different population dynamics models were tested (constant population size, exponential population growth, expansion population growth, logistic population growth and Bayesian skyline). Statistical uncertainty in the data was reflected by the 95% highest probability density (HDP) values. Convergence was assessed with Effective Sample Size (ESS) values, after a burning of 6 million steps. Comparison of the values obtained for marginal likelihoods as well as ESS of these models revealed that the Expansion Population Growth model was the best to fit the data.

The results shown in Table 1 are the outcome of the analysis for 60 million steps of the MCMC, using the GTR model, a relaxed clock [24] and the Expansion Population Growth model [25].

As can be seen in Table 1, our results suggest that the NA gene of the 2009 H1N1 emerging IAV strains evolved from ancestors that existed around August 17th, 2008. This is in agreement with previous results situating the most recent common ancestor (MRCA) for the NA gene of 2009 H1N1 IAV around August 8th, 2008 [26].

When the GTR model is used, a mean of 7.84 × 10-3 nucleotide substitutions per site per year (s/s/y) was obtained for the NA gene (Table 1). This rate is roughly comparable to previous estimations of IAV NA evolutionary rates (3.6 × 10-3 s/s/y) [26]. Interestingly, a significant contribution of the first codon position to the evolutionary rate was also found (Table 1). Moreover, an important expansion growth rate was observed (see Table 1).

Table 1 Bayesian coalescent inference of full-length NA sequences from 2009 H1N1 Influenza A virus strains.

Phylogenetic tree analysis of NA genes from 2009 H1N1 IAV strains

To study the phylogenetic relations among the NA genes from the 62 H1N1 IAV strains enrolled in these studies, maximum clade credibility trees were generated using software from the BEAST package [23]. The results of these studies are shown in Figure 1.

Figure 1
figure 1

Bayesian MCMC phylogenetic tree analysis of 62 NA genes from 2009 H1N1 IAV strains. A maximum credibility clade obtained using the GTR model, the expansion population growth model and a relaxed clock (uncorrelated exponential) is shown. Strains in the tree are shown by name. Main genetic sub-branches are indicated by capital letters (A through D). Node ages are shown in days at the nodes of the tree. The tree is rooted to theirMRCA. Bar at the bottom of the tree show time in days. Strains carrying the H275Y, that confers resistance to Oseltamivir, are shown in red.

As it can be seen in the figure, different genetic sub-branches can be observed. Interestingly, Oseltamivir-resistant viruses can be observed in all main genetic sub-branches. These viruses are situated on the tip of the trees suggesting a recent emergence from the 2009 H1N1 IAV populations (see Figure 1). This is in agreement with the initial studies revealing that 2009 H1N1 IAV strains were susceptible to Oseltamivir and the recent selection of resistant viruses from these viral populations [17].

Mapping of positive-selected and co-evolving sites in a 3D NA protein model

An homology-based 3D structure model of the NA protein of 2009 H1N1 IAV strains have been very recently obtained [27] (available at http://mendel.bii.a-star.edu.sg/SEQUENCES/H1N1/). In order to observe if the amino acids substitutions naturally occurring in the NA genes of the 62 H1N1 IAV studied were associated to previously identified antigenic regions or the active site of the NA protein (being the latter the binding cavity of Oseltamivir and other NA inhibitors drugs), we mapped all substitutions found in NA proteins of all IAV enrolled in these studies in a temporal order, according to the date of isolation of each strain. The results of these studies are shown in Figure 2.

Figure 2
figure 2

Mapping of naturally occurring amino acid substitutions in a NA protein 3D structure. The 3D structure model of the NA protein from 2009 H1N1 IAV shown in the figure was obtained by Mauer-Stroh et al. [27] (Bioinformatic Institute, A*STAR's Biomedical Sciences Institutes, Singapore). Oseltamivir atoms are shown in red. Antibodies binding sites are shown in green. Position 275, where substitution H275Y confers resistance to Oseltamivir, is shown in pink. The substitutions found among strains isolated during 30, 60, 90 (where viruses with H275Y substitution also arise) and 119 days (from March 30th, 2009) are shown in yellow, blue, orange and white in A through D, respectively. Dotted white lines show distances in Å.

As it can be seen in the figure, no substitution was found to be related to the active site of the NA protein (see Figure 2). Importantly, none of the substitutions found in our dataset appears sufficiently close to affect the drug binding pocket (see Figure 2).

Discussion

The antigenic variability of IAV is the basis for the recurring epidemics each year [28]. IAV presents a moving antigenic target, evading immunity triggered by previous infections. For these reasons, efforts to characterize epidemic variants [29] are deemed important for improving influenza vaccine formulation, since the closer the vaccine strain is to the dominant variant, the more effective the vaccine [30].

A sudden emergence of new H1N1 IAV of swine origin is taking place since April of 2009 [15]. This pandemic started in Mexico and it is currently spreading to all regions of the world [14]. On June 11th, the WHO officially raised the phase of pandemic alert to level 6. As of July 19th, 137,232 cases of the 2009 H1N1 IAV emerging strains have been officially confirmed in 142 countries [31].

Different approaches have been extremely useful in increasing our understanding of the spatial-temporal transmission dynamics of influenza. They have also provided assistance in evaluating the potential severity of IAV pandemics, where severity was defined by the value of the Basic Reproduction Number (R0) [32].

The R0 for novel influenza A (H1N1) has recently been estimated to be between 1.4 and 1.6 [33], revealing an important expansion of this IAV population. Fortunately, this value is below values of R0 estimated for the 1918-1919 pandemic strain (mean R0~2, range 1.4 to 2.8) [32].

These results are in agreement with the results found in this work using a Bayesian coalescent MCMC approach (see Table 1). A high expansion growth rate (66.43 new infections/individual/year) was achieved, particularly considering the short period of time studied (March 30th to July 28th, 2009). These results suggest that the pandemic caused by the 2009 H1N1 IAV will continue its expansion phase at a significant rate.

We estimated that the NA of the 2009 H1N1 IAV evolved from ancestors that existed around August 17th, 2008 (Table 1). Interestingly, this date is in agreement with first estimations of the MRCA for that gene of these pandemic strains (August 8th, 2008) [26]. This result suggests that the NA gene segment of these viruses were presumably circulating in the swine reservoir before emerging into the human population, in agreement with recent results [34].

The first estimation of evolutionary rate for the NA gene of the 2009 H1N1 IAV strains established a rate of 3.65 × 10-3[26]. In this study, a mean evolutionary rate of 7.84 × 10-3 s/s/y was obtained (see Table 1). Although not entirely dissimilar rates are found, the possible differences among the two estimations may be due to the fact that the first estimations were carried out at the beginning of the pandemic outbreak, where only 30 NA sequences isolated over a shorter time span (from March to May) were available [26]. In this work, 62 full-length NA sequences, isolated from March to the end of July, were employed (see Table S1, Additional file 1). Importantly, a contribution of first codon position of 0.97 (from a total of 3.0) to the mean evolutionary rate was found (Table 1). This speaks of a comparatively higher contribution of non-synonymous substitutions to the mean substitution rate. This result is in agreement with previous reports showing a comparatively higher non-synonymous to synonymous (dn/ds) substitution rate ratio in the 2009 H1N1 IAV strains [26]. Moreover, a contribution of first codon position to main evolutionary rate like the one found in this study is significantly higher than the ones previously found in other RNA viruses, like Hepatitis A virus (0.33, VP1 gene) [35] and Noroviruses (0.55, VP1 gene) [36].

Maximum clade credibility trees revealed a rapid diversification of NA genes in at least four main phylogenetic lineages (Figure 1). Nevertheless, due to the fact that the degree of genetic variation among all strains included in these analysis is roughly low (with a maximum degree of variation of 0.64%), more studies will be needed to confirm these findings. Interestingly, Oseltamivir resistance was found only in the more recent samples and two of them, the A/Washington/28/2009 and the A/Washington/29/2009, appear to be phylogenetically and geographically linked (Figure 1). This finding also suggests that anti-viral resistant viruses can emerge in any genetic lineage, as a result of selection of mutant viruses from the viral population. Oseltamivir-resistant viruses are situated on the tips of the tree. This reveals a recent emergence from previously susceptible viruses (Figure 1). This result is in agreement with initial studies showing the susceptibility of the H1N1 IAV emerging strains to this drug and its widespread use to combat the spread of these viruses all around the world [17].

Mapping of substitutions found in the 62 NA proteins during the four month period covered by this study revealed that substitutions are distributed all around the surface of the molecule, leaving the hydrophobic core and the catalytic site essentially untouched (see Fig. 2). Nevertheless, strains carrying the H275Y substitution were also observed in these studies (see Figs 1 and 2). Very recent studies on the structure of the NA of mutant IAV strains carrying this substitution, revealed that the bulkier Tyr residue alters the orientation of the key Glu 277 residue [37]. On binding Oseltamivir, the conformation of the Glu 277 side chain of the wild type enzyme is altered such that it exposes a hydrophobic site with which the pentyloxy group of Oseltamivir interacts [37]. In the mutant enzyme, the bulkier Tyr residue at position 275 displaces the carboxyl group of Glu 277 into the binding site, such that it disrupts the hydrophobic pocket and causes a change in conformation of the pentyloxy substituent of Oseltamivir, with consequent reduction in affinity of binding of some 300-fold or greater [38].

Interestingly, this is not the case of Zanamivir, since the H275Y substitution causes only a small shift in the position of Glu 277, without disrupting the H-bonds between Glu 277 and the glycerol moiety of the drug [38]. This suggests that other NA inhibitors, like Zanamivir or Peramivir should still be effective against this H1N1 IAV strains. Importantly, substitutions observed in the 62 patients enrolled in this study suggest that changes at possible antigenic sites at NA protein surface may indeed occur (see Figure 2). For that reason, vaccines to previous strains or acquired immunity from previous IAV infections are expected to be less effective. More detailed studies on the 2009 H1N1 IAV evolution are extremely needed in order to select appropriate IAV vaccine strains.

Conclusion

A coalescent Bayesian Markov Chain Montecarlo (MCMC) approach was used to analyze 62 full-length NA gene sequences from 2009 H1N1 IAV strains, isolated from March 30th to July 28th, 2009. When the Expansion Population Growth model was employed a high rate of evolutionary change of 7.84 × 10-3 s/s/y was obtained for the NA gene. Importantly, a significant contribution of the first codon position to the mean evolutionary rate was also found. Moreover, an important expansion growth rate of 66.43 new infections/individual/year was also observed. Taking these results together, high evolutionary rates and fast population growth have contributed to the initial transmission dynamics of 2009 H1N1 IAV. Naturally occurring substitutions are preferentially located at the protein surface and do not interfere with the NA active site. Antigenic regions relevant for vaccine development can differ from previous vaccine strains and vary among patients.

References

  1. Neumann G, Brownlee GG, Fodor E, Kawaoka Y: Orthomyxovirus replication, transcription, and polyadelynation. Curr Top Microbiol Immunol 2004, 283: 121-143.

    CAS  PubMed  Google Scholar 

  2. Nguyen-Van-Tam JS, Hampson AW: The epidemiology and clinical impact of pandemic influenza. Vaccine 2003, 21: 1762-1768. 10.1016/S0264-410X(03)00069-0

    Article  PubMed  Google Scholar 

  3. Wolf YI, Viboud C, Holmes EC, Koonin EV, Lipman DJ: Long intervals of stasis punctuated by bursts of positive selection in the seasonal evolution of influenza A virus. Biology Direct 2006, 1: 34. 10.1186/1745-6150-1-34

    Article  PubMed Central  PubMed  Google Scholar 

  4. Hillerman MR: Realities and enigmas of human viral influenza: pathogenesis, epidemiology and control. Vaccine 2002, 20: 3068-3087. 10.1016/S0264-410X(02)00254-2

    Article  Google Scholar 

  5. De Jong JC, Rimmelzwaan GF, Fouchier RA, Osterhaus AD: Influenza virus: a master of metamorphosis. J Infect 2000, 40: 218-228. 10.1053/jinf.2000.0652

    Article  CAS  PubMed  Google Scholar 

  6. Ferguson NM, Galvani AP, Bush RM: Ecological and immunological determinants of influenza evolution. Nature 2003, 422: 428-433. 10.1038/nature01509

    Article  CAS  PubMed  Google Scholar 

  7. Kilboune ED, Smith C, Brett I, Pokorny BA, Johansson B, Cox N: The total influenza vaccine failure of 1947 revisited: major intrasubtypic antigenic change can explain failure of vaccine in a post-World War II epidemic. Proc Natl Acad Sci USA 2002, 99: 10748-10752. 10.1073/pnas.162366899

    Article  Google Scholar 

  8. Smith GJD, Bahl J, Vijaykrishna D, Zhang J, Poon LLM, Chen H, Webster RG, Malik JS, Guan Y: Dating the emergence of pandemic influenza viruses. Proc Natl Acad Sci USA 2009, 106: 11709-11712. 10.1073/pnas.0904991106

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Scholtissek C, Rohde W, Von Hoyningen V, Rott R: On the origin of the human influenza virus subtype H2N2 and H3N2. Virology 1978, 87: 13-20. 10.1016/0042-6822(78)90153-8

    Article  CAS  PubMed  Google Scholar 

  10. Kawaoka Y, Krauss S, Webster RG: Avian-to-human transmission of the PB1 gene of influenza A viruses in the 1957 and 1968 pandemics. J Virol 1989, 63: 4603-4608.

    PubMed Central  CAS  PubMed  Google Scholar 

  11. Scholtissek C, von Hoyningen V, Rott R: Genetic relatedness between the new 1977 epidemic strains (H1N1) of influenza and human influenza strains isolated between 1947 to 1986 as determined by oligonucleotide mapping and sequencing studies. J Gen Virol 1978, 70: 299-313.

    Google Scholar 

  12. Cox NJ, Black RA, Kendal AP: Pathways of evolution of influenza A (H1N1) viruses from 1977 to 1986 as determined by oligonucleotide mapping and sequencing studies. J Gen Virol 1978, 70: 299-313. 10.1099/0022-1317-70-2-299

    Article  Google Scholar 

  13. Lin YP, Gregory V, Bennett M, Hay A: Recent changes among human influenza viruses. Virus Res 2004, 103: 47-52. 10.1016/j.virusres.2004.02.011

    Article  CAS  PubMed  Google Scholar 

  14. Centers for Disease Control and Prevention: Update: infections with a swine-origin influenza A (H1N1) virus - United States and other countries, April 28, 2009. MMWR Morb Mortal Wkly Rep 2009, 58: 431-433.

    Google Scholar 

  15. Zaracostas J: World Health Organization declares A (H1N1) influenza pandemic. BMJ 2009, 338: b2425. 10.1136/bmj.b2425

    Article  Google Scholar 

  16. Aoki FY, Boivin G, Roberts N: Influenza virus susceptibility and resistance to oseltamivir. Antivir Ther 2007, 12: 603-616.

    CAS  PubMed  Google Scholar 

  17. Centers for Disease Control and Prevention: Oseltamivir-resistant 2009 pandemic influenza A (H1N1) virusinfection in two summer campers receiving prophylaxis--North Carolina, 2009. MMWR Morb Mortal Wkly Rep 2009, 58: 969-972.

    Google Scholar 

  18. Bao Y, Bolotov D, Dernovoy B, Kiryutin L, Zaslavsky L, Tatusova T, Ostell J, Lipman D: The Influenza Virus Resource at the National Center for Biotechnology Information. J Virol 2008, 82: 596-601. 10.1128/JVI.02005-07

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics 2004, 5: 113. 10.1186/1471-2105-5-113

    Article  PubMed Central  PubMed  Google Scholar 

  20. Posada D, Crandall KA: Selecting the best-fit model of nucleotide substitution. Syst Biol 2001, 50: 580-601. 10.1080/106351501750435121

    Article  CAS  PubMed  Google Scholar 

  21. Lole KS, Bollinger RC, Parnjape RS, Gadkari D, Kulkarni SS: Full length human immunodeficiency virus type I genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol 1999, 73: 152-160.

    PubMed Central  CAS  PubMed  Google Scholar 

  22. Salminen MO, Carr JK, Burke DS, McCutchan FE: Identification of break-points in intergenotypic recombinants of HIV type I by bootscanning. AIDS Res Hum Retroviruses 1995, 11: 1423-1425. 10.1089/aid.1995.11.1423

    Article  CAS  PubMed  Google Scholar 

  23. Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 2007, 7: 214. 10.1186/1471-2148-7-214

    Article  PubMed Central  PubMed  Google Scholar 

  24. Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol 2006, 4: 88. 10.1371/journal.pbio.0040088

    Article  Google Scholar 

  25. Drummond AJ, Rambaut A, Shapiro B, Pybus OG: Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol 2005, 22: 1185-1192. 10.1093/molbev/msi103

    Article  CAS  PubMed  Google Scholar 

  26. Smith GJ, Vijaykrishna D, Bahl J, Lycett SJ, Worobey M, Pybus OG, Ma SK, Cheung CL, Raghwani J, Bhatt S, Peiris JS, Guan Y, Rambaut A: Origins and evolutionary genomics of the 2009 swine-origin H1N1 influenza A epidemic. Nature 2009, 459: 1122-1125. 10.1038/nature08182

    Article  CAS  PubMed  Google Scholar 

  27. Maurer-Stroh S, Ma J, Tze Chuen Lee R, Sirota FL, Eisenhaver F: Mapping the sequence mutations of the 2009 H1N1 influenza A virus neuraminidase relative to drug and antibody binding sites. Biology Direct 2009, 4: 18. 10.1186/1745-6150-4-18

    Article  PubMed Central  PubMed  Google Scholar 

  28. De Jong JC, Rimmelzwaan GF, Fouchier RA, Osterhaus AD: Influenza virus: a master of metamorphosis. J Infect 2000, 40: 218-228. 10.1053/jinf.2000.0652

    Article  CAS  PubMed  Google Scholar 

  29. Wolf YI, Viboud C, Holmes EC, Koonin EV, Lipman DJ: Long intervals of stasis punctuated by bursts of positive selection in the seasonal evolution of influenza A virus. Biology Direct 2006, 1: 34. 10.1186/1745-6150-1-34

    Article  PubMed Central  PubMed  Google Scholar 

  30. Nelson MI, Viboud C, Simonsen L, Bennett RT, Grieserner SB, St George K, Taylor J, Spiro DJ, Sengamalay NA, Ghedin E, Taubenberger JK, Holmes EC: Multiple reassortment events in the evolutionary history of H1N1 Influenza A Virus since 1918. PLoS Pathog 2008, 4: e1000012. 10.1371/journal.ppat.1000012

    Article  PubMed Central  PubMed  Google Scholar 

  31. Balcan D, Hu H, Goncalves B, Bajardi P, Poletto C, Ramasco JJ, Paolotti D, Perra N, Tizzoni M, Broeck W, Colizza V, Vestignani A: Seasonal transmission potential and activity peaks of the new influenza A (H1N1): a Monte Carlo likelihood analysis based on human morbility. BMC Med 2009, 7: e45. 10.1186/1741-7015-7-45

    Article  Google Scholar 

  32. Coburn BJ, Wangner BG, Blower S: Modeling influenza epidemics and pandemics: insights into the future of swine flu (H1N1). BMC Med 2009, 7: e30. 10.1186/1741-7015-7-30

    Article  Google Scholar 

  33. Fraser C, Donnelly CA, Cauchmez S, Hanage WP, Van Kerkhove MD, Hollingsworth TD, Griffin J, Baggaley RF, Jenkins HE, Lyons EJ, Jombart T, Hinsley WR, Grassly NC, Balloux F, Ghani AC, Ferguson NM, Rambaut A, Pybus OG, Lopez-Gatell H, Alpuche-Aranda CM, Chapela IB, Zavala EP, Guevara DM, Checchi F, Garcia E, Hugonnet S, Roth C, WHO Rapid Pandemic Assessment Collaboration: Pandemic potential of a strain of Influenza A (H1N1): early findings. Science 2009, 324: 1557-1561. 10.1126/science.1176062

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Garten RJ, Davis CT, Russell CA, Shu B, Lindstrom S, Balish A, Sessions WM, Xu X, Skepner E, Deyde V, Okomo-Adhiambo M, Gubareva L, Barnes J, Smith CB, Emerly SL, Hillman MJ, Rivailler P, Smagala J, de Graaf M, Burke DF, Fouchier RA, Pappas C, Alpuche-Aranda CM, Lopez-Gattel H, Olivera H, Lopez I, Myers CA, Faix D, Blair PJ, Yu C, Keene KM, Dotson PD Jr, Boxrud D, Sambol AR, Abid SH, St George K, Bannerman T, Moore AL, Stringer DJ, Blevins P, Demmler-Harrison GJ, Ginsberg M, Kriner P, Waterman S, Smole S, Guevara HF, Belongia EA, Clark PA, Beatrice ST, Donis R, Katz J, Finelli L, Bridges CB, Shaw M, Jernigan DB, Uyeki TM, Smith DJ, Klimov AI, Cox NJ: Antigenic and genetic characteristics of swine-origin 2009 A (H1N1) influenza viruses circulating in humans. Science 2009, 325: 197-201. 10.1126/science.1176225

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Moratorio G, Costa-Mattioli M, Piovani R, Romero H, Musto H, Cristina J: Bayesian coalescent inference of hepatitis A virus populations: evolutionary rates and patterns. J Gen Virol 2007, 88: 3039-3042. 10.1099/vir.0.83038-0

    Article  CAS  PubMed  Google Scholar 

  36. Victoria M, Miagostovich MP, Ferreira MS, Vieira CB, Fioretti JM, Leite JP, Colina R, Cristina J: Bayesian coalescent inference reveals high evolutionary rates and expansion of Norovirus populations. Infect Genet Evol 9: 927-932. 10.1016/j.meegid.2009.06.014

  37. Collins PJ, Haire LF, Lin YP, Liu J, Russell RJ, Walker PA, Skehel JJ, Martin SR, Hay AJ, Gamblin SJ: Crystal structures of oseltamivir-resistant influenza virus neuraminidase mutants. Nature 2008, 453: 1258-1261. 10.1038/nature06956

    Article  CAS  PubMed  Google Scholar 

  38. Collins PJ, Haire LF, Lin YP, Liu J, Russell RJ, Walker PA, Martin SR, Daniels RS, Gregory V, Skehel JJ, Gamblin SJ, Hay AJ: Structural basis for oseltamivir resistance of influenza viruses. Vaccine 2009, 27: 6317-6323. 10.1016/j.vaccine.2009.07.017

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We acknowledge support by International Atomic Energy Agency, through Research Contract No. 15792 and Comisión Sectorial de Investigación Científica (CSIC), Universidad de la República, Uruguay, through I+D Project "Variabilidad Genética y Evolución Viral de Virus Influenza A en Uruguay". NG and JC acknowledge support by Agencia Nacional de Investigación e Innovación (ANII) and PEDECIBA, Uruguay.

We thank anonymous reviewers for important advice in improvement of this work.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Juan Cristina.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

NG and JC conceived the study. AF, GM and JC designed and performed the Bayesian coalescent studies and phylogenetic analysis. NG, AF, GM and RC contributed to the discussion and interpretation of the results. JC wrote the paper. All authors read and approved the final manuscript.

Electronic supplementary material

12985_2009_729_MOESM1_ESM.DOC

Additional file 1:Origins of the NA sequences from 2009 H1N1 IAV strains. A table describing the names, date of isolation and accession numbers of all IAV strains included in this study. (DOC 32 KB)

12985_2009_729_MOESM2_ESM.DOC

Additional file 2:FindModel results for NA genes of 2009 H1N1 IAV strains. A table describing the results found for different evolutionary models tested in this study. (DOC 33 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Goñi, N., Fajardo, A., Moratorio, G. et al. Modeling gene sequences over time in 2009 H1N1 Influenza A Virus populations. Virol J 6, 215 (2009). https://doi.org/10.1186/1743-422X-6-215

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1743-422X-6-215

Keywords