Skip to main content

Evolutionary genomics of the pandemic 2009 H1N1 influenza viruses (pH1N 1v)

Abstract

Background

A new strain of human H1N1 influenza A viruses was broken out in the April 2009 and caused worldwide pandemic emergency. The present study is trying to estimate a temporal reassortment history of 2009 H1N1 viruses by phylogenetic analysis based on a total 394 sequences of H1N1viruses isolated from swine, human and avian.

Results

Phylogenetic trees of eight gene segments showed that viruses sampled from human formed a well-supported clade, whereas swine and avian lineages were intermixed together. A new divergence swine sublineage containing gene segments of 2009 H1N1 viruses was characterized, which were closely related with swine viruses collected from USA and South Korea during 2004 to 2007 in six segments (PB2, PB1, PA, HA, NP and NS), and to swine viruses isolated from Thailand during 2004 to 2005 in NA and M. Substitution rates were varied drastically among eight segments and the average substitution rate was generally higher in 2009 H1N1 than in swine and human viruses (F 2 ,23 = 5.972, P < 0.01). Similarly, higher dN/dS substitution ratios were identified in 2009 H1N1 than in swine and human viruses except M2 gene (F2, 25 = 3.779, P < 0.05). The ages of 2009 H1N1 viruses were estimated around 0.1 to 0.5 year, while their common ancestors with closest related swine viruses existed between 9.3 and 17.37 years ago.

Conclusion

Our results implied that at least four reassortments or transmissions probably occurred before 2009 H1N1 viruses. Initial reassortment arose in 1976 and avian-like Eurasian swine viruses emerged. The second transmission happened in Asia and North America between 1988 and 1992, and mostly influenced six segments (PB2, PB1, PA, HA, NP and NS). The third reassortment occurred between North American swine and avian viruses during 1998 to 2000, which involved PB2 and PA segments. Recent reassortments occurred among avian-to-swine reassortant, Eurasian and classical swine viruses during 2004 to 2005. South Korea, Thailand and USA, were identified as locations where reassortments most likely happened. The co-circulation of multiple swine sublineages and special lifestyle in Asia might have facilitated mixing of diverse influenza viruses, leading to generate a novel virus strain.

Background

In April 2009, a new strain of human H1N1 influenza A viruses was identified in Mexico and USA, which caused a rapid global spread and moved the global pandemic alert level to phase 6 on 11 June 2009 [1]. This virus has been widely spread over 170 countries from April to June 2009 [1], and has the ability to persist in the human population, potentially with more severe clinical consequences [2]. Initial genetic characterization of the 2009 H1N1 viruses suggested swine as their probable source, on the basis of sequence similarity to previously reported swine influenza viruses [3–5]. The NA and M gene segments are in the Eurasian swine genetic lineage. The other six gene segments are well clustered into North American triple-reassortant swine lineage. More specially, the HA, NP and NS gene segments belong to classical swine lineage; while the PB2, PB1 and PA gene segments are derived from the North American swine, avian and human H3N2 reassortant lineage [4–6].

Swine has been believed to play a vital role in interspecies transmission of influenza viruses, since it harbors receptors to both avian and human influenza virus strains. Pigs have been considered as a possible "mixing vessel" in which genetic material can be exchanged, with the potential to result in novel progeny viruses to which humans are immunological naive and highly susceptible [7]. Sporadic cases of swine influenza in human, combined with seroepidemiological studies demonstrate increased risk of swine influenza in occupationally exposed workers, highlighting the crucial role that swine may play in the development of new strains of influenza viruses [8–10]. Notably, human infection with H1N1 swine influenza has been a nationally notifiable disease in United State since 2007 [11].

H1N1 influenza viruses were first isolated from swine in 1930 and have been shown to be antigenical highly similar to a recently reconstructed human 1918 A (H1N1) virus [12, 13], and likely share a common ancestor [14, 15]. From 1930 to late 1990s these classical swine influenza viruses circulated in swine and remained antigenically stable [4, 16, 17]. In 1998, a new reassortant swine virus, comprising genes from classical swine viruses and North American avian viruses, was reported as the cause of outbreaks in North American swine influenza, with subsequent establishment in pig populations. In Europe, avian H1N1 viruses were introduced to pigs and first detected in Belgium in 1979. This avian-like swine lineage emerged and subsequently prevailed in Asia after that time [18].

Given the evolutionary history of swine influenza viruses, it is likely that multiple reassortment events resulted in emergence of the novel 2009 human H1N1 viruses. However, the poor surveillance for swine influenza viruses fails to detect these potential reassortants, and where reassortment events most likely happened is also unclear. Considering the ongoing health risk posed by 2009 H1N1 viruses, it is very important to estimate evolutionary genomics of the novel pandemic viruses; to identify the potential reassortment events and where they probably occurred. We therefore undertook a detailed phylogenetic analysis in a total 394 sequences of H1N1 viruses taken from swine, human and avian from 1918 to 2009 to estimate a temporal reconstruction of reassortment history of 2009 H1N1 viruses. By employed a sophisticated Bayesian Markov Chain Monte Carlo (MCMC) approach, the time to the common ancestor and population dynamics of 2009 H1N1 viruses and their closest related swine viruses were also determined.

Materials and methods

Sequence data collection and alignment

H1N1 influenza viruses isolated from swine, human and avian hosts between 1918 and 2009 were studied. Their genomic sequences were accessed from NCBI before May 2009, and aligned using Mega 4.0 [19]. Eight genome segment alignment datasets (PB2, PB1, PA, HA, NP, NA, M and NS) were generated.

Phylogenetic analyses

Phylogenetic trees of eight alignment datasets were reconstructed using the ML approach implemented in PyhML 2.4 [20], and with bootstrapping analyses of 1,000 pseudo-replicate datasets. The robustness of the ML tree topology was assessed by comparing the ML topology with the topologies sampled in the MCMC analyses performed in MrBayes version 3.1.2 [21, 22]. The substitution model for each gene dataset was selected according to MrModeltest 2.2 [23]. Phylogenies were rooted with the H1N1 human virus in 1918.

Estimating the rate of evolution and genetic diversity dynamics

Overall rates of evolutionary change (nucleotide substitutions per site per year) for swine, human, and 2009 human H1N1 viruses were estimated using the BEAST program 1.4.6 [24], utilizing the number and temporal distribution of genetic differences among viruses sampled at different times [25]. Estimate of evolution rate was not performed for avian viruses because limited sample size (n = 47). Previous analyses of human influenza A viruses evolution found that the uncorrelated exponential relaxed clock model provided a better fit to the H1N1 influenza viruses data [24], thus it was employed in estimating the average rates of nucleotide substitution.

To estimate the date of origin for 2009 H1N1 viruses, we preformed a Bayesian skyline plot (BSP) for each gene segment to estimate the time of the most recent common ancestor (TMRCA). TMRCA of 2009 H1N1 viruses and their closet related swine viruses was also estimated for comparisons. All estimates incorporated the General Time Reversible (GTR) and Γ4 model of DNA substitution as this model was consistently with the best supported one in MrModeltest 2.2. The uncorrelated exponential relaxed clock model was selected in estimations. The uncertainty in the estimates reflected in the 95% highest probability density (HPD), and in each case chain lengths were run for sufficient time to achieve convergence.

Measurement of selection pressures

To determine the overall selection pressure faced by each gene segment in 2009 H1N1 viruses, we compared the dN/dS rates of this lineage with those of swine and human virus lineages. Swine and human viruses after 1976 were selected in this analysis because either widely sampling was available or the earliest origin of the 2009 H1N1 viruses could go back to avian-to-swine reassortantment in 1976. We estimated the mean numbers of nonsynonymous substitutions (dN) and synonymous substitutions (dS) per site using the single-likelihood ancestor counting (SLAC) method accessed through the Datamonkey interface (http://www.datamonkey.org). In all cases, dN/dS estimates were based on neighbour-Joining trees under the substitution model selected by website.

Results

Phylogeny of H1N1 influenza viruses and genetic origin of 2009 H1N1 human viruses

Phylogenetic trees of all eight segments demonstrated viruses sampled from human formed a well-supported clade in all phylogenies with the maximum likelihood support rate of 80 to 100%, whereas swine and avian lineages were intermixed together (Figure 1, Additional file 1). Viruses sampled from swine in North America and Asia tended to cluster together (classical swine sublineage), and formed a paraphyly topology with human clade. Swine viruses isolated from Europe and Asia after 1977 clustered with European avian viruses and nested into avian lineage (Eurasian swine sublineage). In PB2 and PA, some swine viruses were nested into North American avian clade, of which were mostly sampled after 1998 and isolated from North America and Asia (avian-to-swine reassortant sublineage). The PB1 had undergone more complex reassortments, which occurred among North American avian, human and swine viruses (triple reassortant swine sublineage) (Additional file 2). Notably, many swine, avian and human H1N1 viruses isolated from 1988 to 1992 intermixed together in all eight gene segments (Figure 1 and Figure 2). These mixed clades were often linked by a long branch with their relatives (Figure 1 and Additional file 1).

Figure 1
figure 1

Phylogenetic trees of H1N1 influenza viruses. ML phylogenies reconstructed from (A) HA gene; (B) NA gene; (C) PA gene. Human lineage is showed in red; while swine and avian are characterized in green and blue, respectively. Only the topological supports for major lineage or sublineages and reassortant viruses are showed and summarized from 1,000 ML bootstrap replicates. Transmissions among swine, avian and human viruses between 1988 and 1992 are highlighted in blue. The closest related swine viruses of 2009 H1N1 viruses are highlighted in red. Phylogenetic trees for five other gene segments can be found in Additional file 1.

Figure 2
figure 2

Reconstructed reassortment events leading to the emergence of 2009 H1N1 influenza viruses.

Phylogenetic trees characterized a new divergence swine sublineage containing gene segments of 2009 H1N1 viruses, which were apparently mostly similar to classical swine sublineage in HA, NP and NS, to avian-like Eurasian swine sublineage in NA and M, to avian-to-swine reassortant swine sublineage in PB2 and PA, and to triple reassorant swine sublineage in PB1 (Figure 1 and Figure 2). Within new divergence swine sublineage, 2009 H1N1 viruses were most closely related with A/swine/USA/2005, A/swine/USA/2007, A/swine/South Korea/2004 and A/swine/South Korea/2005 in six gene segments (PB2, PB1, PA, HA, NP and NS). Comparisons of the nucleotide similarity showed 98 to 99% homology (Additional file 3). These ancestral relatives contained genes from both classical and avian-to-swine reassortant swine sublineages, indicating reassortments occurring in 2004 and 2005, particularly in North America and Asia (Figure 3). For NA and M, the closest related viruses were A/swine/Thailand/2005 and A/swine/Thailand/2004, and their similarity was around 97 to 98% (Additional file 4).

Figure 3
figure 3

Reassortment events involved with the genetic parents of the 2009 H1N1 influenza viruses. Asia, particularly Thailand, was identified as the place where reassortments between Eurasian and classical swine sublineages occurred; alternatively, Asia and North America, especially as South Korea and USA, was where reassortments between classical and avian-to-swine reassortant sublineages happened.

Substitution rate estimations and population dynamics

The substitution rates were varied among eight gene segments, ranging from 3.51* 10-3 for NS to 1.16* 10-3 for M in swine viruses; and from 2.88 * 10-3 for PB1 to 0.99* 10-3 for M in human viruses (Figure 4). The average substitution rate was generally higher in the 2009 H1N1 viruses than in swine and human viruses (F 2 ,23 = 5.972, P < 0.01). However, differences were more significant in HA, NP and NA and M (F 2 ,11 = 6.668, P < 0.05), but not in others (F 2 ,11 = 2.12, P = 0.176).

Figure 4
figure 4

Rates of nucleotide substitution in H1N1 influenza viruses. Mean substitution rates are showed for each gene segment, with 95% lower and upper HPD values shown as error bars. Blue shows the 2009 H1N1 viruses, while red and green indicate swine and human viruses, respectively.

The ages of the 2009 H1N1 viruses were estimated extremely recent (range 0.1 to 0.5 year). The common ancestors of 2009 H1N1 and closest related swine viruses existed between 9.3 and 17.37 years ago (Additional file 4, Table S1). BSP estimation showed an abrupt rise in population dynamic after a long-term constant population size in each gene segment (Figure 5 and Additional file 5).

Figure 5
figure 5

Estimations population dynamics of 2009 H1N1 viruses and their closest related swine viruses. Green, blue and yellow shows HA, NA and PA, respectively, which represented three different origins, especially as classical, Eurasian and avian-to-swine reassortant swine sublineages. Population dynamic estimates for classical, Eurasian and avian-to-swine reassortant swine sublineages can be found in Additional file 5. X-axis: time in years before present; Y-axis: estimated population size [units = Neτ, the product of effective population size and generation length in years (log-transformed)]. The median estimate and both 95% HPD limits are indicated.

Diversifying selection in 2009, swine and human H1N1 virus lineages

The higher dN/dS rates were observed in HA, NP and NS segments in 2009 H1N1 viruses (Figure 6). Similarly, HA and NS showed higher dN/dS values in swine and human virus lineages, in addition to NA and M (Figure 6). In general, 2009 H1N1 viruses showed a comparatively higher dN/dS substitution rate than those in swine and human viruses except in M2 gene (F2, 25 = 3.779, P < 0.05).

Figure 6
figure 6

Comparisons of d N / d S ratios among 2009 H1N1, swine and human viruses. Blue showed 2009 H1N1 viruses; red and green indicated swine and human viruses, respectively. The dN/dS ratio estimate could not be performed for NS1 gene of 2009 H1N1 viruses because only few variable amino acid sites are available.

Discussion

Multiple reassortment events resulting in new pandemic 2009 H1N1 human viruses

Previous studies have classified new emerged 2009 H1N1 viruses a swine origin [4, 5]. Sublineages of swine viruses in our phylogenetic trees belonged to different lineages of hosts, which suggest that multiple reassortment events occurred. It is thus important to recognize which reassortant groups are closely related to 2009 H1N1 viruses. These ancestral groups might give us clues of when and where these reassortant events have probably happened [4].

Phylogenetic analyses revealed that the initial reassortment event might occur in 1976, while gene segments of NA and M of 2009 H1N1 viruses derived from the Eurasian H1N1 swine sublineage. In 1976, an avian H1N1 virus strain was introduced to pigs and established a new avian-like H1N1 swine sublineage (H1N1 Eurasian swine) and circulated in European and Asian countries [4, 5]. Between 1988 and 1992, swine, avian and human viruses frequently intermixed together, which mostly occurred in Asia and North America. These transmissions might influence PB2, PB1, PA, HA, NP and NS genes of 2009 H1N1 viruses. Furthermore, 2009 H1N1 viruses have undergone another reassortment between North American swine and avian viruses. This event occurred between 1998 and 2000, and gene segments of PB2 and PA derived from this avian-to-swine reassortant sublineage. The PB1 has experienced more complex reassortments, first of which occurred between North American avian and human viruses during 1978 to 1979. The new avian-to-human reassortant viruses underwent subsequently reassortments with North American swine viruses during 2000s. Investigations for the closest related ancestors of 2009 H1N1 viruses revealed recent reassortments occurred both between avian-to-swine reassortant and classical swine sublineages, and between classical and Eurasian swine sublineages. These results thus imply that 2009 H1N1 viruses might have acquired gene segments from different groups of parental H1N1 viruses, which underwent at least four reassortments and transmissions among swine sublineages or swine, avian and human lineages.

Where reassortment events most likely happened

Detection of the closest ancestral genes for each of eight segments of 2009 H1N1 viruses may suggest where reassortment events most likely happened. Reassortments between classical and avian-to-swine reassortant swine sublineages occurred in South Korea and USA, while reassortments between Eurasian and classical swine sublineages were found in Thailand. These thus indicated North America and Asia might be the place where reassortments probably occurred. Classical and avian-to-swine reassortant swine sublineages have co-circulated in North America [26, 27] and frequent reassortments between two sublineages have been detected and resulted in the generation of triple reassortant swine H1N1 and H1N2 viruses [4]. In Asia, the classical, Eurasian avian-like and avian-to-swine reassortant swine influenza viruses co-circulate, and extensive reassortments among three sublineages have been observed [5], with reassortments between Eurasian avian-like and reassortant swine sublineages occurring as early as 2003 (for example, Hong Kong). The co-circulation among multiple swine sublineages might have promoted opportunities of reassortments and tended to create the new influenza virus lineage.

Asia, especially Southeast Asia, is thought to be the epicenter for the influenza pandemics throughout history [28]. The special environment and lifestyle in Southeast Asia provide more chances for wild aquatic birds, domestic poultry, pigs and human to contact closely, and create the opportunities for interspecies transmission and generation of new reassortment influenza viruses.

Population dynamics of 2009 H1N1 viruses and their closest related viruses

The common ancestors of 2009 H1N1 viruses and their closest related swine viruses might have been existed 9 to 17 years before; in contrast, the currently sampled 2009 H1N1 viruses shared a common ancestor around January 2009. It is possible that ancestors of 2009 H1N1 viruses have been cryptically circulating for several years before emergence in human. BSP analyses revealed a rapid growth in early 2001 after more continuous and relatively slow growth in a decade. In recent years, both classical and avian-to-swine reassortant swine influenza viruses have occasionally been isolated from humans [7, 8, 29], however, the lack of systematic swine virus surveillance allowed for undetected persistence and evolution of this potentially pandemic strain for many years. Thus systematic monitor and surveillance of H1N1 swine viruses for these regions, particularly Asian countries, should be required for future controlling of the outbreaks of new strain influenza viruses.

Different evolutionary backgrounds in 2009 H1N1, human and swine influenza viruses

2009 H1N1 viruses showed comparatively higher substitution rates than those of human and swine viruses. Similarly, relatively high ratios of non-synonymous to synonymous substitutions per site (dN/dS) were also observed for all segments of 2009 H1N1 viruses except M2 gene. The difference was more significant in HA, NS and NP genes, most likely reflecting different selection pressures at eight gene segments [30].

The difference among evolutionary rates of swine, human and 2009 viruses might reflect various evolutionary backgrounds. Human and swine H1N1 viruses were characterized by low rates of evolutionary change, particularly at amino acid-changing sites, thus proposed to have reached an 'evolutionary stasis' [31]. These observations demonstrate that selective pressure is less intense, so that there is little selective requirement to repeatedly fix amino acid changes that evade host's immune responses [32]. However, a high evolutionary rate in the newly emerged 2009 H1N1 viruses might be due to adaptation to the new host species, indicating strong selective advantages. Alternatively, it could also be a general feature of intensively sampled emerging epidemics, which should be considered in the other evolutionary analyses such as the surveillance of the 2009 H1N1 viruses breakout.

Conclusions

Considered ongoing health risk of 2009 H1N1 influenza viruses, it is necessary to detect potential reassortments involved in emergence of the novel viruses through the phylogenetic analyses. Results in the present study imply that at least four reassortments or transmissions probably occurred before H1N1 broke out. South Korea, Thailand and USA, were identified as locations where reassortments most likely happened. The co-circulation of multiple swine sublineages and special lifestyle in Asia might have facilitated mixing of diverse swine, avian and human viruses, leading to generate novel influenza viruses. 2009 H1N1 viruses showed comparatively higher selection than swine and human viruses in all gene segments except for M2. Our results also revealed that ancestors of 2009 H1N1 viruses existed between 9 and 17 years ago, suggesting that this virus strain has been circulating undetected for decade years and the importance of systematic surveillance of influenza viruses in swine, especially in Asian areas.

References

  1. WHO: Pandemic (H1N1). 2009, [http://www.wpro.who.int/media_centre/press_releases/pr_20091206.htm]

    Google Scholar 

  2. Munster VJ, de Wit E, van den Brand JMA, Herfst S, Schrauwen EJA, Bestebroer TM, van de Vijver D, Boucher CA, Koopmans M, Rimmelzwaan GF, Kuiken T, Osterhaus ADME, Fouchier RAM: Pathogenesis and transmission of swine-origin 2009 A(H1N1) influenza virus in ferrets. Science. 2009, 325: 481-483.

    CAS  PubMed  Google Scholar 

  3. Dawood FS, Jain S, Finelli L, Shaw MW, Lindstrom S, Garten RJ, Gubareva LV, Xu XY, Bridges CB, Uyeki TM: Emergence of a novel swine-origin influenza A (H1N1) virus in humans novel swine-origin influenza A (H1N1) virus investigation team. New England Journal of Medicine. 2009, 360: 2605-2615.

    Article  PubMed  Google Scholar 

  4. Garten RJ, Davis CT, Russell CA, Shu B, Lindstrom S, Balish A, Sessions WM, Xu XY, Skepner E, Deyde V, Okomo-Adhiambo M, Gubareva L, Barnes J, Smith CB, Emery SL, Hillman MJ, Rivailler P, Smagala J, Graaf M, Burke DF, Fouchier RAM, Pappas C, Alpuche-Aranda CM, Lopez-Gatell H, Olivera H, Lopez I, Myers CA, Faix D, Blair PJ, Yu C, et al: 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 

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

    Article  CAS  PubMed  Google Scholar 

  6. Bi YH, Fu GH, Peng JS, Sun YP, Wang JG, Pu J, Zhang Y, Gao HJ, Ma GP, Tian FL, Brown IH, Liu JH: Novel swine influenza virus reassortants in pigs, China. Emerging Infectious Disease. 2010, 16: 1162-1164.

    Article  Google Scholar 

  7. Webster RG, Bean WJ, Gorman OT, Chambers TM, Kawaoka Y: Evolution and ecology of influenza-A viruses. Microbiological Reviews. 1992, 56: 152-179.

    PubMed Central  CAS  PubMed  Google Scholar 

  8. Myers KP, Olsen CW, Gray GC: Cases of swine influenza in humans: A review of the literature. Clinical Infectious Diseases. 2007, 44: 1084-1088. 10.1086/512813.

    Article  PubMed Central  PubMed  Google Scholar 

  9. Trifonov V, Khiabanian H, Greenbaum B, Rabadan R: The origin of the recent swine influenza A(H1N1) virus infecting humans. Eurosurveillance. 2009, 14:

    Google Scholar 

  10. Van Reeth K: Avian and swine influenza viruses: our current understanding of the zoonotic risk. Veterinary Research. 2007, 38: 243-260. 10.1051/vetres:2006062.

    Article  CAS  PubMed  Google Scholar 

  11. Shinde V, Bridges CB, Uyeki TM, Shu B, Balish A, Xu XY, Lindstrom S, Gubareva LV, Deyde V, Garten RJ, Harris M, Gerber S, Vagoski S, Smith F, Pascoe N, Martin K, Dufficy D, Ritger K, Conover C, Quinlisk P, Klimov A, Bresee JS, Finelli L: Triple-reassortant swine influenza A (H1) in humans in the United States, 2005-2009. New England Journal of Medicine. 2009, 360: 2616-2625. 10.1056/NEJMoa0903812.

    Article  CAS  PubMed  Google Scholar 

  12. Shope R, Lewis P: Swine influenza: experimental transmission and pathology. Journal of Experimental Medcine. 1931, 54: 349-359. 10.1084/jem.54.3.349.

    Article  CAS  Google Scholar 

  13. Tumpey TM, Garcia-Sastre A, Taubenberger JK, Palese P, Swayne DE, Basler CF: Pathogenicity and immunogenicity of influenza viruses with genes from the 1918 pandemic virus. Proceedings of the National Academy of Sciences of the United States of America. 2004, 101: 3166-3171. 10.1073/pnas.0308391100.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Gorman OT, Bean WJ, Kawaoka Y, Donatelli I, Guo YJ, Webster RG: Evolution of influenza-A virus nucleoprotein genes - implications for the origins of H1N1 human and classical swine viruses. Journal of Virology. 1991, 65: 3704-3714.

    PubMed Central  CAS  PubMed  Google Scholar 

  15. Reid AH, Taubenberger JK: The origin of the 1918 pandemic influenza virus: a continuing enigma. Journal of General Virology. 2003, 84: 2285-2292. 10.1099/vir.0.19302-0.

    Article  CAS  PubMed  Google Scholar 

  16. Sheerar MG, Easterday BC, Hinshaw VS: Antigenic conservation of H1N1 swine influenza-viruses. Journal of General Virology. 1989, 70: 3297-3303. 10.1099/0022-1317-70-12-3297.

    Article  PubMed  Google Scholar 

  17. Vincent AL, Lager KM, Ma WJ, Lekcharoensuk P, Gramer MR, Loiacono C, Richt JA: Evaluation of hemagglutinin subtype 1 swine influenza viruses from the United States. Veterinary Microbiology. 2006, 118: 212-222. 10.1016/j.vetmic.2006.07.017.

    Article  CAS  PubMed  Google Scholar 

  18. Liu JH, Bi YH, Qin K, Fu GH, Yang J, Peng JS, Ma GP, Liu QF, Pu J, Tian FL: Emergence of European avian-like H1N1 swine influenza A viruses in China. Journal of Clinical Microbiology. 2009, 47: 2643-2646. 10.1128/JCM.00262-09.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Molecular Biology and Evolution. 2007, 24: 1596-1599. 10.1093/molbev/msm092.

    Article  CAS  PubMed  Google Scholar 

  20. Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Systematic Biology. 2003, 52: 696-704. 10.1080/10635150390235520.

    Article  PubMed  Google Scholar 

  21. Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.

    Article  CAS  PubMed  Google Scholar 

  22. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.

    Article  CAS  PubMed  Google Scholar 

  23. Nylander J: MrModeltest v.1.0. Program distributed by the author. 2002, Uppsala: Evolutionary Biology Centre, Uppsala University

    Google Scholar 

  24. Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolutionary Biology. 2007, 7:

    Google Scholar 

  25. Drummond AJ, Nicholls GK, Rodrigo AG, Solomon W: Estimating mutation parameters, population history and genealogy simultaneously from temporally spaced sequence data. Genetics. 2002, 161: 1307-1320.

    PubMed Central  CAS  PubMed  Google Scholar 

  26. Jung K, Song DS: Evidence of the co-circulation of influenza H1N1, H1N2 and H3N2 viruses in the pig population of Korea. Veterinary Record. 2007, 161: 104-105. 10.1136/vr.161.3.104.

    Article  CAS  PubMed  Google Scholar 

  27. Peiris JSM, Guan Y, Markwell D, Ghose P, Webster RG, Shortridge KF: Cocirculation of avian H9N2 and contemporary "human" H3N2 influenza A viruses in pigs in southeastern China: Potential for genetic reassortment?. Journal of Virology. 2001, 75: 9679-9686. 10.1128/JVI.75.20.9679-9686.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Shortridge KF, Stuartharris CH: An influenza epicenter. Lancet. 1982, 2: 812-813.

    Article  CAS  PubMed  Google Scholar 

  29. Hay AJ, Gregory V, Douglas AR, Lin YP: The evolution of human influenza viruses. Philosophical Transactions of The Royal Society of London Series B-Biological Sciences. 2001, 356: 1861-1870.

    Article  CAS  Google Scholar 

  30. Horimoto T, Kawaoka Y: Influenza: Lessons from past pandemics, warnings from current incidents. Nature Reviews Microbiology. 2005, 3: 591-600. 10.1038/nrmicro1208.

    Article  CAS  PubMed  Google Scholar 

  31. Plotkin JB, Dushoff J, Levin SA: Hemagglutinin sequence clusters and the antigenic evolution of influenza A virus. Proceedings of the National Academy of Sciences of the United States of America. 2002, 99: 6263-6268. 10.1073/pnas.082110799.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Suarez DL: Evolution of avian influenza viruses. Veterinary Microbiology. 2000, 74: 15-27. 10.1016/S0378-1135(00)00161-9.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The research was supported by CAS Innovation Program (KSCX2-YW-R-158), IOZ Innovation Program and NSFC (30925008) to Lei FM, and NSFC program J0930004 to Zhang RYand Cui P.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Ziyuan Duan or Fumin Lei.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

This study was elaborated and launched by FML and ZYD, who leading two research groups on ornithology and genetics. YHQ and RY Z accomplished the data analyses and YHQ narrated the paper. PC did much effort on data collection and editing. GS participated in the work of manuscript preparation and revision. All authors have read and approved the final manuscript.

Electronic supplementary material

12985_2011_1385_MOESM1_ESM.PDF

Additional file 1: Figure S1. Phylogenetic trees of H1N1 influenza viruses. ML phylogenies based on PB1, PB2, NP, MP and NS gene segments respectively. Human lineage is showed in red; while swine and avian are characterized in green and blue, respectively. Transmissions among swine, avian and human viruses between 1988 and 1992 are highlighted in blue. The closest related swine viruses of 2009 H1N1 viruses are highlighted in red. (PDF 9 MB)

12985_2011_1385_MOESM2_ESM.PDF

Additional file 2: Figure S2. Possible reassortments of PB1 in the emergence of 2009 H1N1 viruses. The PB1 has undergone more complex reassortments, which occurred among North American avian, human and swine viruses. (PDF 117 KB)

12985_2011_1385_MOESM3_ESM.PDF

Additional file 3: Table S1. Nucleotide identities (uncorrected Pdistance) between 2009 H1N1 viruses and their closest related swine viruses. (PDF 30 KB)

12985_2011_1385_MOESM4_ESM.DOC

Additional file 4: Table S2. Time of recent common ancestors (TMRCA) for 2009 H1N1 viruses and their closest related swine viruses. (DOC 28 KB)

12985_2011_1385_MOESM5_ESM.pdf

Additional file 5: Figure S3. Population dynamic estimates for classical, Eurasian and avian-to-swine reassortant swine sublineages. X-axis: time in years before present; Y-axis: estimated population size [units = Neτ, the product of effective population size and generation length in years (log-transformed)]. The median estimate and both 95% HPD limits are indicated. (PDF 247 KB)

Authors’ original submitted files for images

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

Qu, Y., Zhang, R., Cui, P. et al. Evolutionary genomics of the pandemic 2009 H1N1 influenza viruses (pH1N 1v). Virol J 8, 250 (2011). https://doi.org/10.1186/1743-422X-8-250

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1743-422X-8-250

Keywords