- Open Access
Evolutionary genomics of the pandemic 2009 H1N1 influenza viruses (pH1N 1v)
Virology Journalvolume 8, Article number: 250 (2011)
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.
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.
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.
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 . This virus has been widely spread over 170 countries from April to June 2009 , and has the ability to persist in the human population, potentially with more severe clinical consequences . 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 . 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 .
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 .
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 . Eight genome segment alignment datasets (PB2, PB1, PA, HA, NP, NA, M and NS) were generated.
Phylogenetic trees of eight alignment datasets were reconstructed using the ML approach implemented in PyhML 2.4 , 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 . 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 , utilizing the number and temporal distribution of genetic differences among viruses sampled at different times . 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 , 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.
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).
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).
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).
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).
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).
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 .
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 . 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 , 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 . 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 .
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' . 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 . 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.
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.
WHO: Pandemic (H1N1). 2009, [http://www.wpro.who.int/media_centre/press_releases/pr_20091206.htm]
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.
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.
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.
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.
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.
Webster RG, Bean WJ, Gorman OT, Chambers TM, Kawaoka Y: Evolution and ecology of influenza-A viruses. Microbiological Reviews. 1992, 56: 152-179.
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.
Trifonov V, Khiabanian H, Greenbaum B, Rabadan R: The origin of the recent swine influenza A(H1N1) virus infecting humans. Eurosurveillance. 2009, 14:
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.
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.
Shope R, Lewis P: Swine influenza: experimental transmission and pathology. Journal of Experimental Medcine. 1931, 54: 349-359. 10.1084/jem.54.3.349.
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.
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.
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.
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.
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.
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.
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.
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.
Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Nylander J: MrModeltest v.1.0. Program distributed by the author. 2002, Uppsala: Evolutionary Biology Centre, Uppsala University
Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolutionary Biology. 2007, 7:
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.
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.
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.
Shortridge KF, Stuartharris CH: An influenza epicenter. Lancet. 1982, 2: 812-813.
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.
Horimoto T, Kawaoka Y: Influenza: Lessons from past pandemics, warnings from current incidents. Nature Reviews Microbiology. 2005, 3: 591-600. 10.1038/nrmicro1208.
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.
Suarez DL: Evolution of avian influenza viruses. Veterinary Microbiology. 2000, 74: 15-27. 10.1016/S0378-1135(00)00161-9.
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.
The authors declare that they have no competing interests.
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.