- Research article
- Open Access
Phylogeography of Avian influenza A H9N2 in China
© Jin et al.; licensee BioMed Central. 2014
- Received: 1 August 2014
- Accepted: 8 December 2014
- Published: 15 December 2014
During the past two decades, avian influenza A H9N2 viruses have spread geographically and ecologically in China. Other than its current role in causing outbreaks in poultry and sporadic human infections by direct transmission, H9N2 virus could also serve as an progenitor for novel human avian influenza viruses including H5N1, H7N9 and H10N8. Hence, H9N2 virus is becoming a notable threat to public health. However, despite multiple lineages and genotypes that were detected by previous studies, the migration dynamics of the H9N2 virus in China is unclear. Increasing such knowledge would help us better prevent and control H9N2 as well as other future potentially threatening viruses from spreading across China. The objectives of this study were to determine the source, migration patterns, and the demography history of avian influenza A H9N2 virus that circulated in China.
Using Bayesian phylogeography framework, we showed that the H9N2 virus in mainland China may have originated from the Hong Kong Special Administrative Region (SAR). Southern China, most likely the Guangdong province acts as the primary epicentre for multiple H9N2 strains spreading across the whole country, and eastern China, most likely the Jiangsu province, acts as an important secondary source to seed outbreaks. Our demography inference suggests that during the long-term migration process, H9N2 evolved into multiple diverse lineages and then experienced a selective sweep, which reduced its genetic diversity. Importantly, such a selective sweep may pose a greater threat to public health because novel strains confer higher fitness advantages than strains being replaced and could generate new viruses through reassortment.
Our analyses indicate that migratory birds, poultry trade and transportation have all contributed to the spreading of the H9N2 virus in China. The ongoing migration and evolution of H9N2, which poses a constant threat to the human population, highlights the need for a more comprehensive surveillance of wild birds and for the enhancement of biosafety for China’s poultry industry.
- Influenza A virus
In Asia, the first isolation of H9N2 avian influenza virus was in Hong Kong in the mid 1970s. Since then, H9N2 avian influenza viruses have been primarily detected in ducks, at live poultry markets, in Hong Kong . In the mainland of China, the first outbreak of H9N2 chicken influenza occurred in the Guangdong province in 1994 . Subsequently, the virus rapidly spread across many different regions of mainland China and became the most prevalent influenza virus circulating in poultry such as chicken, duck, and quail [3, 4]. While causing sporadic and endemic outbreaks in China, H9N2 viruses have evolved into distinct lineages in land-based poultry. These lineages include the G1-like (represented by Quail/Hong Kong/G1/97) and Ck/Bei-like or Y280-like (represented by Chicken/Beijing/1/94 or Duck/Hong Kong/Y280/97) lineages, which have been predominant in China since the mid-1990s [5, 6]. In addition, the H9N2 virus has expanded its host range to mammalian species. Ck/Bei/94-like viruses were isolated from domestic pigs in Hong Kong in 1998 . H9N2 virus infection in pig farms was also confirmed in Shandong and several other provinces in mainland China . Most importantly, several infectious cases in humans exhibiting mild respiratory disease have been reported since 1997 from Hong Kong and China [9, 10].
Prior phylogenetic analysis has revealed that H9N2 has donated its six internal genes to highly pathogenic avian influenza (HPIV) H5N1 viruses, which caused the Hong Kong H5N1 outbreak in 1997 . Recently, there is increasing evidence that shows that avian H9N2 virus may act as a source for novel human avian influenza viruses. In March, 2013, a novel avian-origin H7N9 influenza virus was identified in eastern China and quickly spread to many provinces and cities, resulting in many human infections with a high case fatality rate . Studies have demonstrated that novel H7N9 viruses derive six internal genes from various avian H9N2 lineages through continuing reassorment [13–15]. On November, 2013, another newly reassorted avian H10N8 virus, which possesses internal gene cassettes recruited from the poultry H9N2 virus, was isolated from a 73 years old woman . Together, these data indicate that avian H9N2 viruses maintain a genetic resource to cause human infections, either by direct transmission or by the generation of novel strains.
The H9N2 avian influenza virus was widely distributed in different regions of China and occasionally jumps hosts and reassorts with other subtypes of influenza virus, posing a severe public health threat . In the past decade, the genetic and antigenic evolution of the H9N2 virus in China is well documented [18, 19]. However, the geographic diffusion of H9N2 in China is not fully understood. Understanding its geographic spread and migration patterns could help us better prepare for epidemics as well as endemics and provide useful guidelines for epidemiological control. In this study, we track the migration of the H9N2 virus over different regions in mainland China and the Hong Kong SAR using a Bayesian phylogeography approach . With a Bayesian phylogeography analysis, we determined the source, migration patterns, and corresponding demography history of the avian influenza A H9N2 virus that circulated in China.
Phylogeography reconstruction of the H9N2 virus in China
The genealogical tree of HA shows that the early Hong Kong clade and their recent descendants are positioned basal to all G1 and Ck/Bei lineages (Figure 1). The NA gene shows a similar pattern with the exception of some topology differences (Figure 1), which was most likely caused by reassortment . Molecular clock analyses of the HA gene shows that the estimated TMRCA (time of most recent common ancestor) for the root of the H9N2 virus in Hong Kong and mainland China was around Nov 1967 (95% HPD [Apr 1960, Jan 1974]), whereas the TMRCA of NA was around Dec 1963 (95% HPD [Mar 1952, Nov 1973]). The time estimation encompasses the earliest records of H9N2 isolation from a previous report .
Spatial reconstruction for HA shows that Hong Kong has the highest root state posterior probability (0.9448). However, as for NA, most probable root state agrees with the Guangdong province (posterior probability 0.3788) instead of Hong Kong (Figure 2). We calculated the Kullback–Leibler (KL) divergence, which estimates the divergence of prior and posterior probabilities of the root state for HA and NA. The NA tree yields lower KL divergence than that of the HA tree (1.3925 vs 2.7747), suggesting a larger uncertainty in the root state estimation. This difference is most likely the result of the deeper root estimation of NA or reassortment in the influenza genome, as these factors increase uncertainty on the root state [20, 24, 25]. Despite the difference, the significantly high posterior root state probability of Hong Kong and larger KL divergence in HA suggests that Hong Kong is most likely the source of the influenza H9N2 migration. Furthermore, we observed that most inferred locations of the trunk nodes in HA and NA trees were the Hong Kong SAR or the Guangdong province, implying that these two places might be responsible for the wide spread of the H9N2 viruses. In addition, either in HA or NA, Hong Kong SAR and Guangdong province share the majority of posterior mass (Figure 2), which highlights both Hong Kong and Guangdong having important roles in H9N2 migration to other locations on mainland China. Thus, our Bayesian phylogeograph analysis indicates that the H9N2 virus in Hong Kong SAR is the source of the influenza H9N2 migration and that neighboring province Guangdong also plays an important role in the virus’s wide geographic spread. Moreover, our time estimation also demonstrates that from approximately 1995 to 2000 (Figure 1), wide migration and vast genetic diversification of the H9N2 virus occurred, during which the virus migrated from Hong Kong and Guangdong to many other interior provinces (Jiangsu, Anhui, Shandong etc.) of China and formed G1 and Ck/Bei lineages (Figure 1).
Spatial and temporal dynamics of H9N2 virus geographic dispersal
Population dynamics during the H9N2 virus geographic diffusion process
The above Bayesian phylogeographic analysis suggests that the Hong Kong SAR is the original source of the H9N2 virus circulating in China. In addition, the results indicate that the Guangdong province in southern China is linked to multiple introductions of the H9N2 virus in other geographical areas of mainland China (Figures 1, 2 and 3). While migrating through the entire country, major lineages, such as G1 and Ck/Bei of the H9N2 virus were generated from 1995 to 2000, as demonstrated by TMRCA estimation. This finding is further supported by the BSP analysis in which the effective population size increased rapidly during this time period (Figure 5). The inference is consistent with previous reports and is in accordance with the time in which the H9N2 virus became widespread in poultry throughout China [5, 28].
Our spatial and temporal dynamics for the geographic spread of the H9N2 virus in China revealed more details about its migration patterns. In the map, extensive migration links were established repeatedly from Guangdong, further proving that Guangdong acted as the primary epicentre. This observation supports the idea that southern China is an influenza epicenter because half of the four human influenza pandemics first appeared here; furthermore, several epidemics, including H5N1 and H9N2, occurred here as well [29–31]. Additionally, the diffusion process involves much more than spread from a primary epicenter, and we found that Jiangsu, in eastern China, may have served as a secondary epicentre, seeding tertiary outbreaks in near and remote locations. Hence, southern China, especially Guangdong, perhaps served as the critical source for the spread of the H9N2 virus in China, whereas eastern China, for example, the Jiangsu province, was also an important secondary source to seed outbreaks.
Guangdong and Jiangsu provinces are near the Pearl River Delta Economic Zone and the Yangtze River Delta Economic Zone, respectively, and both are booming markets with advanced transport infrastructures, such as highways. China is the second-largest producer of poultry meat in the world, after the United States, and these two provinces produce the largest amount of poultry in China [32, 33]. Intensive trade and transportation of poultry from the two provinces may largely facilitate the expansion of H9N2’s geographic distribution. In addition, at present, poultry farming in China mainly comes from small individual backyard farms that are free-range and high-density. Additionally, wet markets are important retail channels in China due to the Chinese preference for fresh meat . Such enormous old farms and live poultry markets in China often lack biosafety measures, creating a high risk environment for the spreading of diseases [35–37].
Wild bird migration seemed to impact H9N2 dispersion through the East Asian-Australian flyway. A large majority of the diffusion links, which were established between areas of southern, eastern, northern and central China (Figure 3), are covered by the East Asian-Australian flyway . Every migrating season, in this large region, hundreds of thousands of waterfowl move to southern China in the winter and return to northern and north-eastern China in spring [39, 40]. We found that the initial spread of the H9N2 virus from south to north coincides with birds’ migration pattern (Figure 3). Moreover, we observed several long distance epidemiological links directly connecting Guangdong and Northeast China (Jilin and Heilongjiang) (Figure 3). All of these observations increase the likelihood that migratory wild waterfowl carrying H9N2 may play a prominent role in the dispersal of the virus.
Although multiple lineages of H9N2 occurred as it spread through the whole country via poultry and wild birds, demographic inference suggests that its genetic diversity dropped off after 2000. In addition, from genealogical trees, we observed that the Ck/Bei lineage is the dominant lineages in China, and it continually evolves into antigenically novel strains (Figure 3). Therefore, we believe that the H9N2 virus in China underwent a selective sweep after H9N2’s wide spread and genetic diversification. The ongoing sweep replaced old virus strains with antigenically novel strains, resulting in the reduction of genetic diversity. Importantly, such a selective sweep may pose a great threat to public health because novel strains confer higher fitness advantages than the strains they replace and could generate new virus through reassortment; for example, the A/brambling/Beijing/16/2012 like H9N2 virus, which evolved recently, generated the novel avian H7N9 virus in China .
Complex migration dynamics, selective sweep, and co-circulation with other types of avian influenza viruses have made H9N2 in China a notable threat to public health [17, 42]. Understanding its diffusion process and migration patterns is a crucial step for developing effective prevention and control strategies for not only epidemic but also pandemic spreads. Current data shows H9N2 mainly circulates in poultry; therefore, enhancing the biosafety of China’s poultry industry would largely halt the spread of H9N2 as well as future potentially threatening viruses. As large areas of China are covered by migratory flyways, the contribution of wild birds in spreading viruses cannot be ignored; more comprehensive surveillances of H9N2 in wild birds needs to be conducted. In the future, we will continue to gain new insights into the viral migration dynamics and evolutionary ecology in China.
In this work, Bayesian phylogeography framework was applied to determine the source, migration patterns, and corresponding demography history of the avian influenza A H9N2 virus that circulated in China. Our results suggest that H9N2 in China originated from Hong Kong and subsequently spread to the Guangdong province and other regions. Southern China acts as the primary epicenter for multiple H9N2 strains migrating throughout the country, while eastern China acts as an important secondary source to seed outbreaks. Our demography inference suggests that during the long-term migration process, H9N2 evolved into multiple diverse lineages and then underwent a selective sweep, which reduced its genetic diversity. Our analyses indicate that migratory birds, poultry trade, and transportation may all have contributed to spreading the H9N2 virus in China.
Sequence data preparation and alignment
All hemagglutinin and neuraminidase gene sequences of H9N2 viruses were downloaded from the NCBI Influenza Virus Resource . To reduce the number of sequences, first, HA sequences less than 95% full length were removed, and the resulting HA sequences were clustered using CD-HIT v4.6  with a threshold level of 0.95. Next, we sampled sequences from each cluster and retained strains whose corresponding neuraminidase sequences were also available and took one strain per year, per location, and per host. Thus, our final dataset consisted of 209 H9N2 isolates for which hemagglutinin and neuraminidase sequences were both available. The coding region of the final HA and NA sequences were aligned using MAFFT v7.058  and then were inspected manually according to the amino acid sequences using Mega v5.05 .
Bayesian phylogeography reconstruction
We inferred time-scaled phylogenies of HA and NA by Bayesian Markov Chain Monte Carlo (MCMC) sampling using BEAST v1.8.0 . We used the SRD06 codon position model and the uncorrelated log-normal relaxed clock model under a Bayesian skyline coalescent tree prior in the MCMC simulations [27, 48, 49]. We used the Bayesian sky line plot with a Piecewise constant model to elucidate the population dynamics of H9N2 viruses.
To infer ancestor location and migration events, we grouped the 209 H9N2 isolates into 21 localities including the Hong Kong SAR of China and the Chinese provinces of Anhui, Fujian, Gansu, Guangdong, Guangxi, Hebei, Henan, Heilongjiang, Hubei, Hunan, Jiangsu, Jiangxi, Jilin, Ningxia, Shaanxi, Sichuan, Shandong, Tibet, Yunnan, and Zhejiang. The Chinese municipalities of Beijing and Shanghai were grouped to Hebei and Jiangsu provinces, respectively as they lie in close proximity to these two provinces. Then, the spatial location reconstruction and viral migration were estimated using the discrete Bayesian phylogeographic method implemented in BEAST v1.8.0, which utilised a continuous time Markov Chain (CTMC) over discrete sampling locations, and applied a Bayesian stochastic search variable selection (BSSVS) model .
For each data set (HA and NA), we performed four independent runs for 50 million generations with sampling every 5000 steps. Convergence and effective sampling size (ESS) of estimates were assessed by visual inspection using Tracer v1.6 . Multiple chains were then combined after a 10% burn-in using LogCombiner v1.8.0 included in the BEAST package. The maximum clade credibility (MCC) trees with temporal and spatial annotation were summarised with a 10% burn-in removed using TreeAnnotator v1.8.0 in the BEAST package and presentation figures were generated with FigTree v1.4.2 .
To provide statistical power for the phylogeographic model, we used the Kullback–Leibler (KL) divergence, which measures the deviation between the posterior and prior probabilities of the root state. To calculate KL divergence, we employed a fixed prior of 1/K, where K is the number of unique states, and with the posterior estimates in Figures 2. Small KL values indicated that the phylogeographic models were not able to produce root state posteriors with strong statistical power .
We also conducted Bayes factor (BF) tests to provide statistical support for transmission routes between different geographic locations using SPREAD v1.0.6 with a BF cutoff of three . The Bayes factor (BF) values represent the difference between the posterior and prior probabilities that the rates between two locations are non-zero. Thus, routes with high BF have large odds that a migration exists between two locations.
Visualizing phylogeographic diffusion
To animate viral dispersal over the time, we converted annotated MCC trees into a keyhole markup language (KML) file using SPREAD v1.0.6 , which can be visualised via Google Earth . Example KML files showing avian influenza H9N2 migration processes inferred from HA and NA genes can be found in the Additional file 1 and Additional file 2.
Availability of supporting data
All gene sequences of H9N2 viruses used in this study were downloaded from NCBI Influenza Virus Resource. Accession numbers of hemagglutinin and neuraminidase gene sequences can be found in additional file 3.
The study was supported by grants from China Ministry of Science and Technology Project 973 (number 2013CB910800) and Military Medical Innovation Project of Academy of Military Medicine Science (number 2012CXJJ023).
- Shortridge KF: Pandemic influenza: a zoonosis?. Semin Respir Infect. 1992, 7 (1): 11-25.PubMedGoogle Scholar
- Chen B, Zhang Z, Chen W: Isolation and preliminary serological characterization of type A influenza viruses from chickens. China J Vet Med. 1994, 22: 3-5.Google Scholar
- Li CJ, Yu KZ, Tian GB, Yu DD, Liu LL, Jing B, Ping JH, Chen HL: Evolution of H9N2 influenza viruses from domestic poultry in Mainland China. Virology. 2005, 340 (1): 70-83. 10.1016/j.virol.2005.06.025.PubMedView ArticleGoogle Scholar
- Li KS, Xu KM, Peiris JSM, Poon LLM, Yu KZ, Yuen KY, Shortridge KF, Webster RG, Guan Y: Characterization of H9 subtype influenza viruses from the ducks of southern China: A candidate for the next influenza pandemic in humans?. J Virol. 2003, 77 (12): 6988-6994. 10.1128/JVI.77.12.6988-6994.2003.PubMed CentralPubMedView ArticleGoogle Scholar
- Xu KM, Smith GJD, Bahl J, Duan L, Tai H, Vijaykrishna D, Wang J, Zhang JX, Li KS, Fan XH, Webster RG, Chen H, Peiris JSM, Guan Y: The genesis and evolution of H9N2 influenza viruses in poultry from southern china, 2000 to 2005. J Virol. 2007, 81 (19): 10389-10401. 10.1128/JVI.00979-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Guan Y, Shortridge KF, Krauss S, Chin PS, Dyrting KC, Ellis TM, Webster RG, Peiris M: Two lineages of H9N2 influenza viruses continue to circulate in land-based poultry in southeastern China. International Congress Series. 2001, 1219: 187-193.View ArticleGoogle Scholar
- 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?. J Virol. 2001, 75 (20): 9679-9686. 10.1128/JVI.75.20.9679-9686.2001.PubMed CentralPubMedView ArticleGoogle Scholar
- Xu CT, Fan WX, Wei R, Zhao HK: Isolation and identification of swine influenza recombinant A/Swine/Shandong/1/2003(H9N2) virus. Microbes Infect. 2004, 6 (10): 919-925. 10.1016/j.micinf.2004.04.015.PubMedView ArticleGoogle Scholar
- Peiris M, Yuen KY, Leung CW, Chan KH, Ip PLS, Lai RWM, Orr WK, Shortridge KF: Human infection with influenza H9N2. Lancet. 1999, 354 (9182): 916-917. 10.1016/S0140-6736(99)03311-5.PubMedView ArticleGoogle Scholar
- Butt KM, Smith GJD, Chen HL, Zhang LJ, Leung YHC, Xu KM, Lim W, Webster RG, Yuen KY, Peiris JSM, Guan Y: Human infection with an avian H9N2 influenza A virus in Hong Kong in 2003. J Clin Microbiol. 2005, 43 (11): 5760-5767. 10.1128/JCM.43.11.5760-5767.2005.PubMed CentralPubMedView ArticleGoogle Scholar
- Guan Y, Shortridge KF, Krauss S, Webster RG: Molecular characterization of H9N2 influenza viruses: Were they the donors of the "internal" genes of H5N1 viruses in Hong Kong?. Proc Natl Acad Sci U S A. 1999, 96 (16): 9363-9367. 10.1073/pnas.96.16.9363.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu J, Xiao H, Wu Y, Liu D, Qi X, Shi Y, Gao GF: H7N9: a low pathogenic avian influenza A virus infecting humans. Curr Opin Virol. 2014, 5: 91-97.PubMedView ArticleGoogle Scholar
- Liu D, Shi W, Shi Y, Wang D, Xiao H, Li W, Bi Y, Wu Y, Li X, Yan J, Liu W, Zhao G, Yang W, Wang Y, Ma J, Shu Y, Lei F, Gao GF: Origin and diversity of novel avian influenza A H7N9 viruses causing human infection: phylogenetic, structural, and coalescent analyses. Lancet. 2013, 381 (9881): 1926-1932. 10.1016/S0140-6736(13)60938-1.PubMedView ArticleGoogle Scholar
- Zhang L, Zhang Z, Weng Z: Rapid Reassortment of Internal Genes in Avian Influenza A(H7N9) Virus. Clin Infect Dis. 2013, 57 (7): 1059-1061. 10.1093/cid/cit414.PubMedView ArticleGoogle Scholar
- Lu J, Wu J, Zeng X, Guan D, Zou L, Yi L, Liang L, Ni H, Zhang X, Zhong H: Continuing reassortment leads to the genetic diversity of Influenza H7N9 in Guangdong China. J Virol. 2014, 88 (15): 8297-306. 10.1128/JVI.00630-14.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen H, Yuan H, Gao R, Zhang J, Wang D, Xiong Y, Fan G, Yang F, Li X, Zhou J, Zou S, Yang L, Chen T, Dong L, Bo H, Zhao X, Zhang Y, Lan Y, Bai T, Dong J, Li Q, Wang S, Zhang Y, Li H, Gong T, Shi Y, Ni X, Li J, Zhou J, Fan J: Clinical and epidemiological characteristics of a fatal case of avian influenza A H10N8 virus infection: a descriptive study. Lancet. 2014, 383 (9918): 714-721. 10.1016/S0140-6736(14)60111-2.PubMedView ArticleGoogle Scholar
- Liu D, Shi W, Gao GF: Poultry carrying H9N2 act as incubators for novel human avian influenza viruses. Lancet. 2014, 383 (9920): 869-869. 10.1016/S0140-6736(14)60386-X.PubMedView ArticleGoogle Scholar
- Xu KM, Li KS, Smith GJD, Li JW, Tai H, Zhang JX, Webster RG, Peiris JSM, Chen H, Guan Y: Evolution and molecular epidemiology of H9N2 influenza A viruses from Quail in southern China, 2000 to 2005. J Virol. 2007, 81 (6): 2635-2645. 10.1128/JVI.02316-06.PubMed CentralPubMedView ArticleGoogle Scholar
- Dong G, Luo J, Zhang H, Wang C, Duan M, Deliberto TJ, Nolte DL, Ji G, He H: Phylogenetic Diversity and Genotypical Complexity of H9N2 Influenza A Viruses Revealed by Genomic Sequence Analysis. Plos One. 2011, 6 (2): e17212-10.1371/journal.pone.0017212.PubMed CentralPubMedView ArticleGoogle Scholar
- Lemey P, Rambaut A, Drummond AJ, Suchard MA: Bayesian Phylogeography finds its roots. Plos Comput Biol. 2009, 5 (9): e1000520-10.1371/journal.pcbi.1000520.PubMed CentralPubMedView ArticleGoogle Scholar
- Wallace RG, Hodac H, Lathrop RH, Fitch WM: A statistical phylogeography of influenza A H5N1. Proc Natl Acad Sci U S A. 2007, 104 (11): 4473-4478. 10.1073/pnas.0700435104.PubMed CentralPubMedView ArticleGoogle Scholar
- Wallace RG, Fitch WM: Influenza A H5N1 Immigration Is Filtered Out at Some International Borders. Plos One. 2008, 3 (2): e1697-10.1371/journal.pone.0001697.PubMed CentralPubMedView ArticleGoogle Scholar
- Homme PJ, Easterday BC: Avian influenza virus infections. I. Characteristics of influenza A-turkey-Wisconsin-1966 virus. Avian Dis. 1970, 14 (1): 66-74. 10.2307/1588557.PubMedView ArticleGoogle Scholar
- Wei K, Chen Y, Xie D: Genome-scale evolution and phylodynamics of H5N1 influenza virus in China during 1996–2012. Vet Microbiol. 2013, 167 (3–4): 383-393.PubMedView ArticleGoogle Scholar
- Scotch M, Mei C, Makonnen YJ, Pinto J, Ali A, Vegso S, Kane M, Sarkar IN, Rabinowitz P: Phylogeography of influenza A H5N1 clade 220.127.116.11 in Egypt. BMC Genomics. 2013, 14: 871-10.1186/1471-2164-14-871.PubMed CentralPubMedView ArticleGoogle Scholar
- Google Earth. [http://earth.google.com]
- Drummond AJ, Rambaut A, Shapiro B, Pybus OG: Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005, 22 (5): 1185-1192. 10.1093/molbev/msi103.PubMedView ArticleGoogle Scholar
- Sun Y, Pu J, Jiang Z, Guan T, Xia Y, Xu Q, Liu L, Ma B, Tian F, Brown EG, Liu J: Genotypic evolution and antigenic drift of H9N2 influenza viruses in China from 1994 to 2008. Vet Microbiol. 2010, 146 (3–4): 215-225.PubMedView ArticleGoogle Scholar
- Shortridge KF: Is China an influenza epicentre?. Chin Med J (Engl). 1997, 110 (8): 637-641.Google Scholar
- Shortridge KF, Stuart-Harris CH: An influenza epicentre?. Lancet. 1982, 2 (8302): 812-813.PubMedView ArticleGoogle Scholar
- Webster RG, Bean WJ, Gorman OT, Chambers TM, Kawaoka Y: Evolution and ecology of influenza A viruses. Microbiol Rev. 1992, 56 (1): 152-179.PubMed CentralPubMedGoogle Scholar
- USDA Foreign Agricultural Services: Livestock and Poultry: World Markets and Trade. 2006, [http://usda.mannlib.cornell.edu/usda/fas/livestock-poultry-ma/2000s/2006/livestock-poultry-ma-10-02-2006.pdf]Google Scholar
- Ministry of Agriculture of the People's Republic of China: China Agricultural Development Report. 2005, Beijing: China Agricultural PressGoogle Scholar
- Poon Coen: An Overview of China’s Poultry industry. [http://edepot.wur.nl/118587]
- Liu M, He SQ, Walker D, Zhou NN, Perez DR, Mo B, Li F, Huang XT, Webster RG, Webby RJ: The influenza virus gene pool in a poultry market in South Central China. Virology. 2003, 305 (2): 267-275. 10.1006/viro.2002.1762.PubMedView ArticleGoogle Scholar
- Wan X-F, Dong L, Lan Y, Long L-P, Xu C, Zou S, Li Z, Wen L, Cai Z, Wang W, Li X, Yuan F, Sui H, Zhang Y, Dong J, Sun S, Gao Y, Wang M, Bai T, Yang L, Li D, Yang W, Yu H, Wang S, Feng Z, Wang Y, Guo Y, Webby RJ, Shu Y: Indications that Live Poultry Markets Are a Major Source of Human H5N1 Influenza Virus Infection in China. J Virol. 2011, 85 (24): 13432-13438. 10.1128/JVI.05266-11.PubMed CentralPubMedView ArticleGoogle Scholar
- Huang K, Zhu H, Fan X, Wang J, Cheung C-L, Duan L, Hong W, Liu Y, Li L, Smith DK, Chen H, Webster RG, Webby RJ, Peiris M, Guan Y: Establishment and Lineage Replacement of H6 Influenza Viruses in Domestic Ducks in Southern China. J Virol. 2012, 86 (11): 6075-6083. 10.1128/JVI.06389-11.PubMed CentralPubMedView ArticleGoogle Scholar
- Boere GC, Stroud DA: The flyway concept: what it is and what it isn’t. Waterbirds Around the World. Edited by: Boere GC, Galbraith CA, Stroud DA. 2006, Edinburgh: The Stationery Office, 40-47.Google Scholar
- Takekawa JY, Newman SH, Xiao X, Prosser DJ, Spragens KA, Palm EC, Yan B, Li T, Lei F, Zhao D, Douglas DC, Muzaffar SB, Ji W: Migration of Waterfowl in the East Asian Flyway and Spatial Relationship to HPAI H5N1 Outbreaks. Avian Dis. 2010, 54 (1): 466-476.PubMedView ArticleGoogle Scholar
- Shi B, Xia S, Yang G-J, Zhou X-N, Liu J: Inferring the potential risks of H7N9 infection by spatiotemporally characterizing bird migration and poultry distribution in eastern China. Infect Dis Poverty. 2013, 2 (1): 8-8. 10.1186/2049-9957-2-8.PubMed CentralPubMedView ArticleGoogle Scholar
- Gao R, Cao B, Hu Y, Feng Z, Wang D, Hu W, Chen J, Jie Z, Qiu H, Xu K, Xu X, Lu H, Zhu W, Gao Z, Xiang N, Shen Y, He Z, Gu Y, Zhang Z, Yang Y, Zhao X, Zhou L, Li X, Zou S, Zhang Y, Li X, Yang L, Guo J, Dong J, Li Q: Human Infection with a Novel Avian-Origin Influenza A (H7N9) Virus. New Engl J Med. 2013, 368 (20): 1888-1897. 10.1056/NEJMoa1304459.PubMedView ArticleGoogle Scholar
- Shi W, Li W, Li X, Haywood J, Ma J, Gao GF, Liu D: Phylogenetics of varied subtypes of avian influenza viruses in China: potential threat to humans. Protein Cell. 2014, 5 (4): 253-257. 10.1007/s13238-014-0036-1.PubMed CentralPubMedView ArticleGoogle Scholar
- Bao Y, Bolotov P, Dernovoy D, Kiryutin B, Zaslavsky L, Tatusova T, Ostell J, Lipman D: The influenza virus resource at the national center for biotechnology information. J Virol. 2008, 82 (2): 596-601. 10.1128/JVI.02005-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Li WZ, Jaroszewski L, Godzik A: Clustering of highly homologous sequences to reduce the size of large protein databases. Bioinformatics. 2001, 17 (3): 282-283. 10.1093/bioinformatics/17.3.282.PubMedView ArticleGoogle Scholar
- Katoh K, Misawa K, Kuma K, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30 (14): 3059-3066. 10.1093/nar/gkf436.PubMed CentralPubMedView ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.PubMed CentralPubMedView ArticleGoogle Scholar
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.PubMed CentralPubMedView ArticleGoogle Scholar
- Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. Plos Biol. 2006, 4 (5): e88-10.1371/journal.pbio.0040088.PubMed CentralPubMedView ArticleGoogle Scholar
- Shapiro B, Rambaut A, Drummond AJ: Choosing appropriate substitution models for the phylogenetic analysis of protein-coding sequences. Mol Biol Evol. 2006, 23 (1): 7-9.PubMedView ArticleGoogle Scholar
- Rambaut A, Suchard MA, Xie D, Drummond AJ: Tracer v1.6. 2014, [http://beast.bio.ed.ac.uk/Tracer]Google Scholar
- Rambaut A: FigTree v1.4.2. 2014, [http://beast.bio.ed.ac.uk/figtree]Google Scholar
- Bielejec F, Rambaut A, Suchard MA, Lemey P: SPREAD: spatial phylogenetic reconstruction of evolutionary dynamics. Bioinformatics. 2011, 27 (20): 2910-2912. 10.1093/bioinformatics/btr481.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.