Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Reassortment patterns of avian influenza virus internal segments among different subtypes</title>
<author><name sortKey="Lu, Lu" sort="Lu, Lu" uniqKey="Lu L" first="Lu" last="Lu">Lu Lu</name>
<affiliation><nlm:aff id="I1">University of Edinburgh, Institute of Evolutionary Biology, Ashworth Laboratories, West Mains Road, Edinburgh EH9 3JT, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Lycett, Samantha J" sort="Lycett, Samantha J" uniqKey="Lycett S" first="Samantha J" last="Lycett">Samantha J. Lycett</name>
<affiliation><nlm:aff id="I1">University of Edinburgh, Institute of Evolutionary Biology, Ashworth Laboratories, West Mains Road, Edinburgh EH9 3JT, UK</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="I2">Current address: University of Glasgow, Institute of Biodiversity, Animal Health and Comparative Medicine, Glasgow G12 8QQ, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Leigh Brown, Andrew J" sort="Leigh Brown, Andrew J" uniqKey="Leigh Brown A" first="Andrew J" last="Leigh Brown">Andrew J. Leigh Brown</name>
<affiliation><nlm:aff id="I1">University of Edinburgh, Institute of Evolutionary Biology, Ashworth Laboratories, West Mains Road, Edinburgh EH9 3JT, UK</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">24456010</idno>
<idno type="pmc">3905155</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3905155</idno>
<idno type="RBID">PMC:3905155</idno>
<idno type="doi">10.1186/1471-2148-14-16</idno>
<date when="2014">2014</date>
<idno type="wicri:Area/Pmc/Corpus">000587</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000587</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Reassortment patterns of avian influenza virus internal segments among different subtypes</title>
<author><name sortKey="Lu, Lu" sort="Lu, Lu" uniqKey="Lu L" first="Lu" last="Lu">Lu Lu</name>
<affiliation><nlm:aff id="I1">University of Edinburgh, Institute of Evolutionary Biology, Ashworth Laboratories, West Mains Road, Edinburgh EH9 3JT, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Lycett, Samantha J" sort="Lycett, Samantha J" uniqKey="Lycett S" first="Samantha J" last="Lycett">Samantha J. Lycett</name>
<affiliation><nlm:aff id="I1">University of Edinburgh, Institute of Evolutionary Biology, Ashworth Laboratories, West Mains Road, Edinburgh EH9 3JT, UK</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="I2">Current address: University of Glasgow, Institute of Biodiversity, Animal Health and Comparative Medicine, Glasgow G12 8QQ, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Leigh Brown, Andrew J" sort="Leigh Brown, Andrew J" uniqKey="Leigh Brown A" first="Andrew J" last="Leigh Brown">Andrew J. Leigh Brown</name>
<affiliation><nlm:aff id="I1">University of Edinburgh, Institute of Evolutionary Biology, Ashworth Laboratories, West Mains Road, Edinburgh EH9 3JT, UK</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">BMC Evolutionary Biology</title>
<idno type="eISSN">1471-2148</idno>
<imprint><date when="2014">2014</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><sec><title>Background</title>
<p>The segmented RNA genome of avian Influenza viruses (AIV) allows genetic reassortment between co-infecting viruses, providing an evolutionary pathway to generate genetic innovation. The genetic diversity (16 haemagglutinin and 9 neuraminidase subtypes) of AIV indicates an extensive reservoir of influenza viruses exists in bird populations, but how frequently subtypes reassort with each other is still unknown. Here we quantify the reassortment patterns among subtypes in the Eurasian avian viral pool by reconstructing the ancestral states of the subtypes as discrete states on time-scaled phylogenies with respect to the internal protein coding segments. We further analyzed how host species, the inferred evolutionary rates and the d<sub>N</sub>
/d<sub>S</sub>
ratio varied among segments and between discrete subtypes, and whether these factors may be associated with inter-subtype reassortment rate.</p>
</sec>
<sec><title>Results</title>
<p>The general patterns of reassortment are similar among five internal segments with the exception of segment 8, encoding the Non-Structural genes, which has a more divergent phylogeny. However, significant variation in rates between subtypes was observed. In particular, hemagglutinin-encoding segments of subtypes H5 to H9 reassort at a lower rate compared to those of H1 to H4, and Neuraminidase-encoding segments of subtypes N1 and N2 reassort less frequently than N3 to N9. Both host species and d<sub>N</sub>
/d<sub>S</sub>
ratio were significantly associated with reassortment rate, while evolutionary rate was not associated. The d<sub>N</sub>
/d<sub>S</sub>
ratio was negatively correlated with reassortment rate, as was the number of negatively selected sites for all segments.</p>
</sec>
<sec><title>Conclusions</title>
<p>These results indicate that overall selective constraint and host species are both associated with reassortment rate. These results together identify the wild bird population as the major source of new reassortants, rather than domestic poultry. The lower reassortment rates observed for H5N1 and H9N2 may be explained by the large proportion of strains derived from domestic poultry populations. In contrast, the higher rates observed in the H1N1, H3N8 and H4N6 subtypes could be due to their primary origin as infections of wild birds with multiple low pathogenicity strains in the large avian reservoir.</p>
</sec>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Webster, Rg" uniqKey="Webster R">RG Webster</name>
</author>
<author><name sortKey="Bean, Wj" uniqKey="Bean W">WJ Bean</name>
</author>
<author><name sortKey="Gorman, Ot" uniqKey="Gorman O">OT Gorman</name>
</author>
<author><name sortKey="Chambers, Tm" uniqKey="Chambers T">TM Chambers</name>
</author>
<author><name sortKey="Kawaoka, Y" uniqKey="Kawaoka Y">Y Kawaoka</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mchardy, Ac" uniqKey="Mchardy A">AC McHardy</name>
</author>
<author><name sortKey="Adams, B" uniqKey="Adams B">B Adams</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Webster, Rg" uniqKey="Webster R">RG Webster</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Parrish, Cr" uniqKey="Parrish C">CR Parrish</name>
</author>
<author><name sortKey="Kawaoka, Y" uniqKey="Kawaoka Y">Y Kawaoka</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kalthoff, D" uniqKey="Kalthoff D">D Kalthoff</name>
</author>
<author><name sortKey="Globig, A" uniqKey="Globig A">A Globig</name>
</author>
<author><name sortKey="Beer, M" uniqKey="Beer M">M Beer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lindstrom, Se" uniqKey="Lindstrom S">SE Lindstrom</name>
</author>
<author><name sortKey="Cox, Nj" uniqKey="Cox N">NJ Cox</name>
</author>
<author><name sortKey="Klimov, A" uniqKey="Klimov A">A Klimov</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Belshe, Rb" uniqKey="Belshe R">RB Belshe</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Smith, Gjd" uniqKey="Smith G">GJD Smith</name>
</author>
<author><name sortKey="Vijaykrishna, D" uniqKey="Vijaykrishna D">D Vijaykrishna</name>
</author>
<author><name sortKey="Bahl, J" uniqKey="Bahl J">J Bahl</name>
</author>
<author><name sortKey="Lycett, Sj" uniqKey="Lycett S">SJ Lycett</name>
</author>
<author><name sortKey="Worobey, M" uniqKey="Worobey M">M Worobey</name>
</author>
<author><name sortKey="Pybus, Og" uniqKey="Pybus O">OG Pybus</name>
</author>
<author><name sortKey="Ma, Sk" uniqKey="Ma S">SK Ma</name>
</author>
<author><name sortKey="Cheung, Cl" uniqKey="Cheung C">CL Cheung</name>
</author>
<author><name sortKey="Raghwani, J" uniqKey="Raghwani J">J Raghwani</name>
</author>
<author><name sortKey="Bhatt, S" uniqKey="Bhatt S">S Bhatt</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Taubenberger, Jk" uniqKey="Taubenberger J">JK Taubenberger</name>
</author>
<author><name sortKey="Kash, Jc" uniqKey="Kash J">JC Kash</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Webster, Rg" uniqKey="Webster R">RG Webster</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cheung, Cl" uniqKey="Cheung C">CL Cheung</name>
</author>
<author><name sortKey="Vijaykrishna, D" uniqKey="Vijaykrishna D">D Vijaykrishna</name>
</author>
<author><name sortKey="Smith, Gj" uniqKey="Smith G">GJ Smith</name>
</author>
<author><name sortKey="Fan, Xh" uniqKey="Fan X">XH Fan</name>
</author>
<author><name sortKey="Zhang, Jx" uniqKey="Zhang J">JX Zhang</name>
</author>
<author><name sortKey="Bahl, J" uniqKey="Bahl J">J Bahl</name>
</author>
<author><name sortKey="Duan, L" uniqKey="Duan L">L Duan</name>
</author>
<author><name sortKey="Huang, K" uniqKey="Huang K">K Huang</name>
</author>
<author><name sortKey="Tai, H" uniqKey="Tai H">H Tai</name>
</author>
<author><name sortKey="Wang, J" uniqKey="Wang J">J Wang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huang, K" uniqKey="Huang K">K Huang</name>
</author>
<author><name sortKey="Bahl, J" uniqKey="Bahl J">J Bahl</name>
</author>
<author><name sortKey="Fan, Xh" uniqKey="Fan X">XH Fan</name>
</author>
<author><name sortKey="Vijaykrishna, D" uniqKey="Vijaykrishna D">D Vijaykrishna</name>
</author>
<author><name sortKey="Cheung, Cl" uniqKey="Cheung C">CL Cheung</name>
</author>
<author><name sortKey="Webby, Rj" uniqKey="Webby R">RJ Webby</name>
</author>
<author><name sortKey="Webster, Rg" uniqKey="Webster R">RG Webster</name>
</author>
<author><name sortKey="Chen, H" uniqKey="Chen H">H Chen</name>
</author>
<author><name sortKey="Smith, Gj" uniqKey="Smith G">GJ Smith</name>
</author>
<author><name sortKey="Peiris, Js" uniqKey="Peiris J">JS Peiris</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huang, K" uniqKey="Huang K">K Huang</name>
</author>
<author><name sortKey="Zhu, H" uniqKey="Zhu H">H Zhu</name>
</author>
<author><name sortKey="Fan, X" uniqKey="Fan X">X Fan</name>
</author>
<author><name sortKey="Wang, J" uniqKey="Wang J">J Wang</name>
</author>
<author><name sortKey="Cheung, Cl" uniqKey="Cheung C">CL Cheung</name>
</author>
<author><name sortKey="Duan, L" uniqKey="Duan L">L Duan</name>
</author>
<author><name sortKey="Hong, W" uniqKey="Hong W">W Hong</name>
</author>
<author><name sortKey="Liu, Y" uniqKey="Liu Y">Y Liu</name>
</author>
<author><name sortKey="Li, L" uniqKey="Li L">L Li</name>
</author>
<author><name sortKey="Smith, Dk" uniqKey="Smith D">DK Smith</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wood, Gw" uniqKey="Wood G">GW Wood</name>
</author>
<author><name sortKey="Mccauley, Jw" uniqKey="Mccauley J">JW McCauley</name>
</author>
<author><name sortKey="Bashiruddin, Jb" uniqKey="Bashiruddin J">JB Bashiruddin</name>
</author>
<author><name sortKey="Alexander, Dj" uniqKey="Alexander D">DJ Alexander</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Klenk, Hd" uniqKey="Klenk H">HD Klenk</name>
</author>
<author><name sortKey="Garten, W" uniqKey="Garten W">W Garten</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Stienekegrober, A" uniqKey="Stienekegrober A">A Stienekegrober</name>
</author>
<author><name sortKey="Vey, M" uniqKey="Vey M">M Vey</name>
</author>
<author><name sortKey="Angliker, H" uniqKey="Angliker H">H Angliker</name>
</author>
<author><name sortKey="Shaw, E" uniqKey="Shaw E">E Shaw</name>
</author>
<author><name sortKey="Thomas, G" uniqKey="Thomas G">G Thomas</name>
</author>
<author><name sortKey="Roberts, C" uniqKey="Roberts C">C Roberts</name>
</author>
<author><name sortKey="Klenk, Hd" uniqKey="Klenk H">HD Klenk</name>
</author>
<author><name sortKey="Garten, W" uniqKey="Garten W">W Garten</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Duan, L" uniqKey="Duan L">L Duan</name>
</author>
<author><name sortKey="Campitelli, L" uniqKey="Campitelli L">L Campitelli</name>
</author>
<author><name sortKey="Fan, Xh" uniqKey="Fan X">XH Fan</name>
</author>
<author><name sortKey="Leung, Yh" uniqKey="Leung Y">YH Leung</name>
</author>
<author><name sortKey="Vijaykrishna, D" uniqKey="Vijaykrishna D">D Vijaykrishna</name>
</author>
<author><name sortKey="Zhang, Jx" uniqKey="Zhang J">JX Zhang</name>
</author>
<author><name sortKey="Donatelli, I" uniqKey="Donatelli I">I Donatelli</name>
</author>
<author><name sortKey="Delogu, M" uniqKey="Delogu M">M Delogu</name>
</author>
<author><name sortKey="Li, Ks" uniqKey="Li K">KS Li</name>
</author>
<author><name sortKey="Foni, E" uniqKey="Foni E">E Foni</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Campitelli, L" uniqKey="Campitelli L">L Campitelli</name>
</author>
<author><name sortKey="Mogavero, E" uniqKey="Mogavero E">E Mogavero</name>
</author>
<author><name sortKey="De Marco, Ma" uniqKey="De Marco M">MA De Marco</name>
</author>
<author><name sortKey="Delogu, M" uniqKey="Delogu M">M Delogu</name>
</author>
<author><name sortKey="Puzelli, S" uniqKey="Puzelli S">S Puzelli</name>
</author>
<author><name sortKey="Frezza, F" uniqKey="Frezza F">F Frezza</name>
</author>
<author><name sortKey="Facchini, M" uniqKey="Facchini M">M Facchini</name>
</author>
<author><name sortKey="Chiapponi, C" uniqKey="Chiapponi C">C Chiapponi</name>
</author>
<author><name sortKey="Foni, E" uniqKey="Foni E">E Foni</name>
</author>
<author><name sortKey="Cordioli, P" uniqKey="Cordioli P">P Cordioli</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lei, F" uniqKey="Lei F">F Lei</name>
</author>
<author><name sortKey="Shi, W" uniqKey="Shi W">W Shi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Claas, Ec" uniqKey="Claas E">EC Claas</name>
</author>
<author><name sortKey="Osterhaus, Ad" uniqKey="Osterhaus A">AD Osterhaus</name>
</author>
<author><name sortKey="Van Beek, R" uniqKey="Van Beek R">R van Beek</name>
</author>
<author><name sortKey="De Jong, Jc" uniqKey="De Jong J">JC De Jong</name>
</author>
<author><name sortKey="Rimmelzwaan, Gf" uniqKey="Rimmelzwaan G">GF Rimmelzwaan</name>
</author>
<author><name sortKey="Senne, Da" uniqKey="Senne D">DA Senne</name>
</author>
<author><name sortKey="Krauss, S" uniqKey="Krauss S">S Krauss</name>
</author>
<author><name sortKey="Shortridge, Kf" uniqKey="Shortridge K">KF Shortridge</name>
</author>
<author><name sortKey="Webster, Rg" uniqKey="Webster R">RG Webster</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chen, H" uniqKey="Chen H">H Chen</name>
</author>
<author><name sortKey="Smith, Gj" uniqKey="Smith G">GJ Smith</name>
</author>
<author><name sortKey="Zhang, Sy" uniqKey="Zhang S">SY Zhang</name>
</author>
<author><name sortKey="Qin, K" uniqKey="Qin K">K Qin</name>
</author>
<author><name sortKey="Wang, J" uniqKey="Wang J">J Wang</name>
</author>
<author><name sortKey="Li, Ks" uniqKey="Li K">KS Li</name>
</author>
<author><name sortKey="Webster, Rg" uniqKey="Webster R">RG Webster</name>
</author>
<author><name sortKey="Peiris, Js" uniqKey="Peiris J">JS Peiris</name>
</author>
<author><name sortKey="Guan, Y" uniqKey="Guan Y">Y Guan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Feare, Cj" uniqKey="Feare C">CJ Feare</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Capua, I" uniqKey="Capua I">I Capua</name>
</author>
<author><name sortKey="Alexander, Dj" uniqKey="Alexander D">DJ Alexander</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Horimoto, T" uniqKey="Horimoto T">T Horimoto</name>
</author>
<author><name sortKey="Kawaoka, Y" uniqKey="Kawaoka Y">Y Kawaoka</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Drummond, Aj" uniqKey="Drummond A">AJ Drummond</name>
</author>
<author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lemey, P" uniqKey="Lemey P">P Lemey</name>
</author>
<author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name>
</author>
<author><name sortKey="Drummond, Aj" uniqKey="Drummond A">AJ Drummond</name>
</author>
<author><name sortKey="Suchard, Ma" uniqKey="Suchard M">MA Suchard</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lemey, P" uniqKey="Lemey P">P Lemey</name>
</author>
<author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name>
</author>
<author><name sortKey="Welch, Jj" uniqKey="Welch J">JJ Welch</name>
</author>
<author><name sortKey="Suchard, Ma" uniqKey="Suchard M">MA Suchard</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Faria, Nr" uniqKey="Faria N">NR Faria</name>
</author>
<author><name sortKey="Suchard, Ma" uniqKey="Suchard M">MA Suchard</name>
</author>
<author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name>
</author>
<author><name sortKey="Lemey, P" uniqKey="Lemey P">P Lemey</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lycett, Sj" uniqKey="Lycett S">SJ Lycett</name>
</author>
<author><name sortKey="Baillie, G" uniqKey="Baillie G">G Baillie</name>
</author>
<author><name sortKey="Coulter, E" uniqKey="Coulter E">E Coulter</name>
</author>
<author><name sortKey="Bhatt, S" uniqKey="Bhatt S">S Bhatt</name>
</author>
<author><name sortKey="Kellam, P" uniqKey="Kellam P">P Kellam</name>
</author>
<author><name sortKey="Mccauley, Jw" uniqKey="Mccauley J">JW McCauley</name>
</author>
<author><name sortKey="Wood, Jl" uniqKey="Wood J">JL Wood</name>
</author>
<author><name sortKey="Brown, Ih" uniqKey="Brown I">IH Brown</name>
</author>
<author><name sortKey="Pybus, Og" uniqKey="Pybus O">OG Pybus</name>
</author>
<author><name sortKey="Leigh Brown, Aj" uniqKey="Leigh Brown A">AJ Leigh Brown</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, Ks" uniqKey="Li K">KS Li</name>
</author>
<author><name sortKey="Guan, Y" uniqKey="Guan Y">Y Guan</name>
</author>
<author><name sortKey="Wang, J" uniqKey="Wang J">J Wang</name>
</author>
<author><name sortKey="Smith, Gj" uniqKey="Smith G">GJ Smith</name>
</author>
<author><name sortKey="Xu, Km" uniqKey="Xu K">KM Xu</name>
</author>
<author><name sortKey="Duan, L" uniqKey="Duan L">L Duan</name>
</author>
<author><name sortKey="Rahardjo, Ap" uniqKey="Rahardjo A">AP Rahardjo</name>
</author>
<author><name sortKey="Puthavathana, P" uniqKey="Puthavathana P">P Puthavathana</name>
</author>
<author><name sortKey="Buranathai, C" uniqKey="Buranathai C">C Buranathai</name>
</author>
<author><name sortKey="Nguyen, Td" uniqKey="Nguyen T">TD Nguyen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dong, G" uniqKey="Dong G">G Dong</name>
</author>
<author><name sortKey="Luo, J" uniqKey="Luo J">J Luo</name>
</author>
<author><name sortKey="Zhang, H" uniqKey="Zhang H">H Zhang</name>
</author>
<author><name sortKey="Wang, C" uniqKey="Wang C">C Wang</name>
</author>
<author><name sortKey="Duan, M" uniqKey="Duan M">M Duan</name>
</author>
<author><name sortKey="Deliberto, Tj" uniqKey="Deliberto T">TJ Deliberto</name>
</author>
<author><name sortKey="Nolte, Dl" uniqKey="Nolte D">DL Nolte</name>
</author>
<author><name sortKey="Ji, G" uniqKey="Ji G">G Ji</name>
</author>
<author><name sortKey="He, H" uniqKey="He H">H He</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dugan, Vg" uniqKey="Dugan V">VG Dugan</name>
</author>
<author><name sortKey="Chen, R" uniqKey="Chen R">R Chen</name>
</author>
<author><name sortKey="Spiro, Dj" uniqKey="Spiro D">DJ Spiro</name>
</author>
<author><name sortKey="Sengamalay, N" uniqKey="Sengamalay N">N Sengamalay</name>
</author>
<author><name sortKey="Zaborsky, J" uniqKey="Zaborsky J">J Zaborsky</name>
</author>
<author><name sortKey="Ghedin, E" uniqKey="Ghedin E">E Ghedin</name>
</author>
<author><name sortKey="Nolting, J" uniqKey="Nolting J">J Nolting</name>
</author>
<author><name sortKey="Swayne, De" uniqKey="Swayne D">DE Swayne</name>
</author>
<author><name sortKey="Runstadler, Ja" uniqKey="Runstadler J">JA Runstadler</name>
</author>
<author><name sortKey="Happ, Gm" uniqKey="Happ G">GM Happ</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Baba, K" uniqKey="Baba K">K Baba</name>
</author>
<author><name sortKey="Shibata, R" uniqKey="Shibata R">R Shibata</name>
</author>
<author><name sortKey="Sibuya, M" uniqKey="Sibuya M">M Sibuya</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name>
</author>
<author><name sortKey="Pybus, Og" uniqKey="Pybus O">OG Pybus</name>
</author>
<author><name sortKey="Nelson, Mi" uniqKey="Nelson M">MI Nelson</name>
</author>
<author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name>
</author>
<author><name sortKey="Taubenberger, Jk" uniqKey="Taubenberger J">JK Taubenberger</name>
</author>
<author><name sortKey="Holmes, Ec" uniqKey="Holmes E">EC Holmes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ludwig, S" uniqKey="Ludwig S">S Ludwig</name>
</author>
<author><name sortKey="Schultz, U" uniqKey="Schultz U">U Schultz</name>
</author>
<author><name sortKey="Mandler, J" uniqKey="Mandler J">J Mandler</name>
</author>
<author><name sortKey="Fitch, Wm" uniqKey="Fitch W">WM Fitch</name>
</author>
<author><name sortKey="Scholtissek, C" uniqKey="Scholtissek C">C Scholtissek</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Treanor, Jj" uniqKey="Treanor J">JJ Treanor</name>
</author>
<author><name sortKey="Snyder, Mh" uniqKey="Snyder M">MH Snyder</name>
</author>
<author><name sortKey="London, Wt" uniqKey="London W">WT London</name>
</author>
<author><name sortKey="Murphy, Br" uniqKey="Murphy B">BR Murphy</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, S" uniqKey="Wang S">S Wang</name>
</author>
<author><name sortKey="Shi, Wm" uniqKey="Shi W">WM Shi</name>
</author>
<author><name sortKey="Mweene, A" uniqKey="Mweene A">A Mweene</name>
</author>
<author><name sortKey="Wei, Hl" uniqKey="Wei H">HL Wei</name>
</author>
<author><name sortKey="Bai, Gr" uniqKey="Bai G">GR Bai</name>
</author>
<author><name sortKey="Liu, Jh" uniqKey="Liu J">JH Liu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Suarez, Dl" uniqKey="Suarez D">DL Suarez</name>
</author>
<author><name sortKey="Perdue, Ml" uniqKey="Perdue M">ML Perdue</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lemey, P" uniqKey="Lemey P">P Lemey</name>
</author>
<author><name sortKey="Minin, Vn" uniqKey="Minin V">VN Minin</name>
</author>
<author><name sortKey="Bielejec, F" uniqKey="Bielejec F">F Bielejec</name>
</author>
<author><name sortKey="Kosakovsky Pond, Sl" uniqKey="Kosakovsky Pond S">SL Kosakovsky Pond</name>
</author>
<author><name sortKey="Suchard, Ma" uniqKey="Suchard M">MA Suchard</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kosakovsky Pond, Sl" uniqKey="Kosakovsky Pond S">SL Kosakovsky Pond</name>
</author>
<author><name sortKey="Frost, Sd" uniqKey="Frost S">SD Frost</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bahl, J" uniqKey="Bahl J">J Bahl</name>
</author>
<author><name sortKey="Vijaykrishna, D" uniqKey="Vijaykrishna D">D Vijaykrishna</name>
</author>
<author><name sortKey="Holmes, Ec" uniqKey="Holmes E">EC Holmes</name>
</author>
<author><name sortKey="Smith, Gj" uniqKey="Smith G">GJ Smith</name>
</author>
<author><name sortKey="Guan, Y" uniqKey="Guan Y">Y Guan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chen, R" uniqKey="Chen R">R Chen</name>
</author>
<author><name sortKey="Holmes, Ec" uniqKey="Holmes E">EC Holmes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cauthen, An" uniqKey="Cauthen A">AN Cauthen</name>
</author>
<author><name sortKey="Swayne, De" uniqKey="Swayne D">DE Swayne</name>
</author>
<author><name sortKey="Schultz Cherry, S" uniqKey="Schultz Cherry S">S Schultz-Cherry</name>
</author>
<author><name sortKey="Perdue, Ml" uniqKey="Perdue M">ML Perdue</name>
</author>
<author><name sortKey="Suarez, Dl" uniqKey="Suarez D">DL Suarez</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, C" uniqKey="Li C">C Li</name>
</author>
<author><name sortKey="Yu, K" uniqKey="Yu K">K Yu</name>
</author>
<author><name sortKey="Tian, G" uniqKey="Tian G">G Tian</name>
</author>
<author><name sortKey="Yu, D" uniqKey="Yu D">D Yu</name>
</author>
<author><name sortKey="Liu, L" uniqKey="Liu L">L Liu</name>
</author>
<author><name sortKey="Jing, B" uniqKey="Jing B">B Jing</name>
</author>
<author><name sortKey="Ping, J" uniqKey="Ping J">J Ping</name>
</author>
<author><name sortKey="Chen, H" uniqKey="Chen H">H Chen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lebarbenchon, C" uniqKey="Lebarbenchon C">C Lebarbenchon</name>
</author>
<author><name sortKey="Feare, Cj" uniqKey="Feare C">CJ Feare</name>
</author>
<author><name sortKey="Renaud, F" uniqKey="Renaud F">F Renaud</name>
</author>
<author><name sortKey="Thomas, F" uniqKey="Thomas F">F Thomas</name>
</author>
<author><name sortKey="Gauthier Clerc, M" uniqKey="Gauthier Clerc M">M Gauthier-Clerc</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gu, M" uniqKey="Gu M">M Gu</name>
</author>
<author><name sortKey="Liu, W" uniqKey="Liu W">W Liu</name>
</author>
<author><name sortKey="Cao, Y" uniqKey="Cao Y">Y Cao</name>
</author>
<author><name sortKey="Peng, D" uniqKey="Peng D">D Peng</name>
</author>
<author><name sortKey="Wang, X" uniqKey="Wang X">X Wang</name>
</author>
<author><name sortKey="Wan, H" uniqKey="Wan H">H Wan</name>
</author>
<author><name sortKey="Zhao, G" uniqKey="Zhao G">G Zhao</name>
</author>
<author><name sortKey="Xu, Q" uniqKey="Xu Q">Q Xu</name>
</author>
<author><name sortKey="Zhang, W" uniqKey="Zhang W">W Zhang</name>
</author>
<author><name sortKey="Song, Q" uniqKey="Song Q">Q Song</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Deng, G" uniqKey="Deng G">G Deng</name>
</author>
<author><name sortKey="Tan, D" uniqKey="Tan D">D Tan</name>
</author>
<author><name sortKey="Shi, J" uniqKey="Shi J">J Shi</name>
</author>
<author><name sortKey="Cui, P" uniqKey="Cui P">P Cui</name>
</author>
<author><name sortKey="Jiang, Y" uniqKey="Jiang Y">Y Jiang</name>
</author>
<author><name sortKey="Liu, L" uniqKey="Liu L">L Liu</name>
</author>
<author><name sortKey="Tian, G" uniqKey="Tian G">G Tian</name>
</author>
<author><name sortKey="Kawaoka, Y" uniqKey="Kawaoka Y">Y Kawaoka</name>
</author>
<author><name sortKey="Li, C" uniqKey="Li C">C Li</name>
</author>
<author><name sortKey="Chen, H" uniqKey="Chen H">H Chen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kida, H" uniqKey="Kida H">H Kida</name>
</author>
<author><name sortKey="Yanagawa, R" uniqKey="Yanagawa R">R Yanagawa</name>
</author>
<author><name sortKey="Matsuoka, Y" uniqKey="Matsuoka Y">Y Matsuoka</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sharp, Gb" uniqKey="Sharp G">GB Sharp</name>
</author>
<author><name sortKey="Kawaoka, Y" uniqKey="Kawaoka Y">Y Kawaoka</name>
</author>
<author><name sortKey="Jones, Dj" uniqKey="Jones D">DJ Jones</name>
</author>
<author><name sortKey="Bean, Wj" uniqKey="Bean W">WJ Bean</name>
</author>
<author><name sortKey="Pryor, Sp" uniqKey="Pryor S">SP Pryor</name>
</author>
<author><name sortKey="Hinshaw, V" uniqKey="Hinshaw V">V Hinshaw</name>
</author>
<author><name sortKey="Webster, Rg" uniqKey="Webster R">RG Webster</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Olsen, B" uniqKey="Olsen B">B Olsen</name>
</author>
<author><name sortKey="Munster, Vj" uniqKey="Munster V">VJ Munster</name>
</author>
<author><name sortKey="Wallensten, A" uniqKey="Wallensten A">A Wallensten</name>
</author>
<author><name sortKey="Waldenstrom, J" uniqKey="Waldenstrom J">J Waldenstrom</name>
</author>
<author><name sortKey="Osterhaus, Ad" uniqKey="Osterhaus A">AD Osterhaus</name>
</author>
<author><name sortKey="Fouchier, Ra" uniqKey="Fouchier R">RA Fouchier</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Edgar, Rc" uniqKey="Edgar R">RC Edgar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Stamatakis, A" uniqKey="Stamatakis A">A Stamatakis</name>
</author>
<author><name sortKey="Ludwig, T" uniqKey="Ludwig T">T Ludwig</name>
</author>
<author><name sortKey="Meier, H" uniqKey="Meier H">H Meier</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Raftery, Manae" uniqKey="Raftery M">MANAE Raftery</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Drummond, Aj" uniqKey="Drummond A">AJ Drummond</name>
</author>
<author><name sortKey="Ho, Sy" uniqKey="Ho S">SY Ho</name>
</author>
<author><name sortKey="Phillips, Mj" uniqKey="Phillips M">MJ Phillips</name>
</author>
<author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tavare, S" uniqKey="Tavare S">S Tavaré</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shapiro, B" uniqKey="Shapiro B">B Shapiro</name>
</author>
<author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name>
</author>
<author><name sortKey="Drummond, Aj" uniqKey="Drummond A">AJ Drummond</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hasegawa, M" uniqKey="Hasegawa M">M Hasegawa</name>
</author>
<author><name sortKey="Kishino, H" uniqKey="Kishino H">H Kishino</name>
</author>
<author><name sortKey="Yano, T" uniqKey="Yano T">T Yano</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pagel, M" uniqKey="Pagel M">M Pagel</name>
</author>
<author><name sortKey="Meade, A" uniqKey="Meade A">A Meade</name>
</author>
<author><name sortKey="Barker, D" uniqKey="Barker D">D Barker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ragonnet Cronin, M" uniqKey="Ragonnet Cronin M">M Ragonnet-Cronin</name>
</author>
<author><name sortKey="Hodcroft, E" uniqKey="Hodcroft E">E Hodcroft</name>
</author>
<author><name sortKey="Hue, S" uniqKey="Hue S">S Hue</name>
</author>
<author><name sortKey="Fearnhill, E" uniqKey="Fearnhill E">E Fearnhill</name>
</author>
<author><name sortKey="Delpech, V" uniqKey="Delpech V">V Delpech</name>
</author>
<author><name sortKey="Leigh Brown, Aj" uniqKey="Leigh Brown A">AJ Leigh Brown</name>
</author>
<author><name sortKey="Lycett, S" uniqKey="Lycett S">S Lycett</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="O Brien, Jd" uniqKey="O Brien J">JD O'Brien</name>
</author>
<author><name sortKey="Minin, Vn" uniqKey="Minin V">VN Minin</name>
</author>
<author><name sortKey="Suchard, Ma" uniqKey="Suchard M">MA Suchard</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ghedin, E" uniqKey="Ghedin E">E Ghedin</name>
</author>
<author><name sortKey="Sengamalay, Na" uniqKey="Sengamalay N">NA Sengamalay</name>
</author>
<author><name sortKey="Shumway, M" uniqKey="Shumway M">M Shumway</name>
</author>
<author><name sortKey="Zaborsky, J" uniqKey="Zaborsky J">J Zaborsky</name>
</author>
<author><name sortKey="Feldblyum, T" uniqKey="Feldblyum T">T Feldblyum</name>
</author>
<author><name sortKey="Subbu, V" uniqKey="Subbu V">V Subbu</name>
</author>
<author><name sortKey="Spiro, Dj" uniqKey="Spiro D">DJ Spiro</name>
</author>
<author><name sortKey="Sitz, J" uniqKey="Sitz J">J Sitz</name>
</author>
<author><name sortKey="Koo, H" uniqKey="Koo H">H Koo</name>
</author>
<author><name sortKey="Bolotov, P" uniqKey="Bolotov P">P Bolotov</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article" xml:lang="en"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">BMC Evol Biol</journal-id>
<journal-id journal-id-type="iso-abbrev">BMC Evol. Biol</journal-id>
<journal-title-group><journal-title>BMC Evolutionary Biology</journal-title>
</journal-title-group>
<issn pub-type="epub">1471-2148</issn>
<publisher><publisher-name>BioMed Central</publisher-name>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">24456010</article-id>
<article-id pub-id-type="pmc">3905155</article-id>
<article-id pub-id-type="publisher-id">1471-2148-14-16</article-id>
<article-id pub-id-type="doi">10.1186/1471-2148-14-16</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Reassortment patterns of avian influenza virus internal segments among different subtypes</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" id="A1"><name><surname>Lu</surname>
<given-names>Lu</given-names>
</name>
<xref ref-type="aff" rid="I1">1</xref>
<email>L.Lu-3@sms.ed.ac.uk</email>
</contrib>
<contrib contrib-type="author" id="A2"><name><surname>Lycett</surname>
<given-names>Samantha J</given-names>
</name>
<xref ref-type="aff" rid="I1">1</xref>
<xref ref-type="aff" rid="I2">2</xref>
<email>samantha.lycett@glasgow.ac.uk</email>
</contrib>
<contrib contrib-type="author" corresp="yes" id="A3"><name><surname>Leigh Brown</surname>
<given-names>Andrew J</given-names>
</name>
<xref ref-type="aff" rid="I1">1</xref>
<email>A.Leigh-Brown@ed.ac.uk</email>
</contrib>
</contrib-group>
<aff id="I1"><label>1</label>
University of Edinburgh, Institute of Evolutionary Biology, Ashworth Laboratories, West Mains Road, Edinburgh EH9 3JT, UK</aff>
<aff id="I2"><label>2</label>
Current address: University of Glasgow, Institute of Biodiversity, Animal Health and Comparative Medicine, Glasgow G12 8QQ, UK</aff>
<pub-date pub-type="collection"><year>2014</year>
</pub-date>
<pub-date pub-type="epub"><day>24</day>
<month>1</month>
<year>2014</year>
</pub-date>
<volume>14</volume>
<fpage>16</fpage>
<lpage>16</lpage>
<history><date date-type="received"><day>9</day>
<month>10</month>
<year>2013</year>
</date>
<date date-type="accepted"><day>14</day>
<month>1</month>
<year>2014</year>
</date>
</history>
<permissions><copyright-statement>Copyright © 2014 Lu et al.; licensee BioMed Central Ltd.</copyright-statement>
<copyright-year>2014</copyright-year>
<copyright-holder>Lu et al.; licensee BioMed Central Ltd.</copyright-holder>
<license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by/2.0"><license-p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0">http://creativecommons.org/licenses/by/2.0</ext-link>
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/publicdomain/zero/1.0/">http://creativecommons.org/publicdomain/zero/1.0/</ext-link>
) applies to the data made available in this article, unless otherwise stated.</license-p>
</license>
</permissions>
<self-uri xlink:href="http://www.biomedcentral.com/1471-2148/14/16"></self-uri>
<abstract><sec><title>Background</title>
<p>The segmented RNA genome of avian Influenza viruses (AIV) allows genetic reassortment between co-infecting viruses, providing an evolutionary pathway to generate genetic innovation. The genetic diversity (16 haemagglutinin and 9 neuraminidase subtypes) of AIV indicates an extensive reservoir of influenza viruses exists in bird populations, but how frequently subtypes reassort with each other is still unknown. Here we quantify the reassortment patterns among subtypes in the Eurasian avian viral pool by reconstructing the ancestral states of the subtypes as discrete states on time-scaled phylogenies with respect to the internal protein coding segments. We further analyzed how host species, the inferred evolutionary rates and the d<sub>N</sub>
/d<sub>S</sub>
ratio varied among segments and between discrete subtypes, and whether these factors may be associated with inter-subtype reassortment rate.</p>
</sec>
<sec><title>Results</title>
<p>The general patterns of reassortment are similar among five internal segments with the exception of segment 8, encoding the Non-Structural genes, which has a more divergent phylogeny. However, significant variation in rates between subtypes was observed. In particular, hemagglutinin-encoding segments of subtypes H5 to H9 reassort at a lower rate compared to those of H1 to H4, and Neuraminidase-encoding segments of subtypes N1 and N2 reassort less frequently than N3 to N9. Both host species and d<sub>N</sub>
/d<sub>S</sub>
ratio were significantly associated with reassortment rate, while evolutionary rate was not associated. The d<sub>N</sub>
/d<sub>S</sub>
ratio was negatively correlated with reassortment rate, as was the number of negatively selected sites for all segments.</p>
</sec>
<sec><title>Conclusions</title>
<p>These results indicate that overall selective constraint and host species are both associated with reassortment rate. These results together identify the wild bird population as the major source of new reassortants, rather than domestic poultry. The lower reassortment rates observed for H5N1 and H9N2 may be explained by the large proportion of strains derived from domestic poultry populations. In contrast, the higher rates observed in the H1N1, H3N8 and H4N6 subtypes could be due to their primary origin as infections of wild birds with multiple low pathogenicity strains in the large avian reservoir.</p>
</sec>
</abstract>
</article-meta>
</front>
<body><sec><title>Background</title>
<p>Influenza A viruses are highly adaptable, able to evade host immune responses and to transmit among different host species. The segmented genome is composed of eight negative sense RNA strands. Segments 4 and 6 encode haemagglutinin (HA) and neuraminidase (NA), respectively. The external domains of HA and NA are the major targets for host neutralizing antibodies and are used in the subtyping of influenza viruses. The other six are internal gene segments: the first three segments encode the viral polymerase complex: PB2, PB1 and PA, and segment 5 encodes the nucleoprotein (NP). These four proteins comprise the ribonucleoproteins (vRNPs), which are the core of the virion, with the viral negative stranded RNAs packaged inside. Segment 7 encodes both matrix protein (M1) and ion channel protein (M2). Two separate Non-Structural proteins NS1 and NS2 are encoded by segment 8 [<xref ref-type="bibr" rid="B1">1</xref>
,<xref ref-type="bibr" rid="B2">2</xref>
].</p>
<p>The structure of the virus genome allows for exchange among the eight RNA segments between viruses co-infecting a cell, a process termed reassortment [<xref ref-type="bibr" rid="B3">3</xref>
]. The viral reassortants derived from different virus subtypes can acquire completely new antigens (antigenic shift) thus avoiding recognition by previously infected hosts and allowing for efficient transmission, leading to a pandemic [<xref ref-type="bibr" rid="B4">4</xref>
]. Three human influenza pandemics in the 20th century had either complete or partial avian origins: the H1N1 Spanish flu in 1918 appears to have been derived from an earlier avian source and both the H2N2 Asian flu in 1957 and the H3N2 Hong Kong Flu in 1968 were reassortants between human and avian strains [<xref ref-type="bibr" rid="B5">5</xref>
-<xref ref-type="bibr" rid="B7">7</xref>
]. Genetic reassortment also played a prominent role in the origin of the swine-origin influenza A (H1N1) 2009 virus [<xref ref-type="bibr" rid="B8">8</xref>
].</p>
<p>Wild avian species of wetlands and aquatic environments such as the Anseriformes (particularly ducks, geese, and swans) and Charadriiformes (particularly gulls, terns, and waders) have been recognized as the major natural reservoirs of influenza A viruses where the viruses are maintained predominantly by asymptomatic birds, indicating a long-standing history of endemic infection [<xref ref-type="bibr" rid="B1">1</xref>
]. Terrestrial poultry, such as chickens and turkeys are not traditionally seen as reservoir hosts of avian influenza virus (AIV), but are susceptible to infection with wild-bird-derived AIV [<xref ref-type="bibr" rid="B9">9</xref>
] through exposure to fecal material from wild birds when they migrate through an area [<xref ref-type="bibr" rid="B10">10</xref>
]. Previous studies indicated that domestic ducks may act as intermediaries between migratory ducks and terrestrial poultry in southern China [<xref ref-type="bibr" rid="B11">11</xref>
,<xref ref-type="bibr" rid="B12">12</xref>
], thus they could be thought as part of the larger reservoir within which different subtypes of influenza interact with each other [<xref ref-type="bibr" rid="B13">13</xref>
].</p>
<p>Avian influenza A virus strains are classified as low pathogenicity (LP) or high pathogenicity (HP) [<xref ref-type="bibr" rid="B5">5</xref>
,<xref ref-type="bibr" rid="B14">14</xref>
]. Highly pathogenic strains primarily arise in subtypes H5 and H7 and are distinguished by the appearance of a highly basic cleavage site in HA [<xref ref-type="bibr" rid="B15">15</xref>
,<xref ref-type="bibr" rid="B16">16</xref>
]. Genetic exchanges between viruses circulating in wild and domestic birds have been documented, usually involving LPAI viruses in domestic poultry [<xref ref-type="bibr" rid="B17">17</xref>
,<xref ref-type="bibr" rid="B18">18</xref>
]. More rarely, HPAI virus can be found to infect both wild and domestic birds, as seen in the case of HPH5N1 avian influenza. Since 2003, HPAI H5N1 AIV spread from south China into Southeast Asia, causing sporadic transmissions to humans and tremendous damage to the poultry industry [<xref ref-type="bibr" rid="B19">19</xref>
,<xref ref-type="bibr" rid="B20">20</xref>
]. Repeated spill-over of Asian HPAI H5N1 viruses from domestic to wild birds has also been observed, which demonstrates the potential for reverse flow [<xref ref-type="bibr" rid="B21">21</xref>
,<xref ref-type="bibr" rid="B22">22</xref>
], and supports a close interaction between influenza strains in wild and domestic birds [<xref ref-type="bibr" rid="B23">23</xref>
].</p>
<p>The 16 HA subtypes and 9 NA subtypes have allowed at least 103 of the possible 144 type A influenza A virus HA-NA combinations to arise [<xref ref-type="bibr" rid="B24">24</xref>
], which indicates a high frequency of reassortment among different HA and NA subtypes. However some subtypes are rarely detected, indicating that some restrictions on possible combinations may exist. Thus, two important questions, regarding the general pattern of genetic interaction between different influenza A virus subtypes remain: whether certain subtypes are more likely to reassort than the other subtypes, and why this might be so.</p>
<p>Recent advances in Bayesian phylogenetics have provided methods which directly link patterns of genetic diversity to ecological processes, including changes in population size and substructure through time, allowing simultaneous insights into the spatial, temporal, and demographic dynamics of rapidly evolving pathogen populations, collectively described as “phylodynamics” [<xref ref-type="bibr" rid="B25">25</xref>
-<xref ref-type="bibr" rid="B27">27</xref>
]. The spatial phylodynamic process can be recovered from genomic data using phylogeographic analyses, allowing us to quantitatively study how and to what extent ecological features shape pathogen genetic diversity [<xref ref-type="bibr" rid="B28">28</xref>
].</p>
<p>Unlike the dynamics of seasonal human influenza strains, characterized by a few co-circulating subtypes, little inter-subtype reassortment and broadly correlated evolution of lineages of the individual segments [<xref ref-type="bibr" rid="B29">29</xref>
], avian influenza has many co-circulating subtypes which caused a complex distribution of internal segment lineages [<xref ref-type="bibr" rid="B30">30</xref>
,<xref ref-type="bibr" rid="B31">31</xref>
]. Conversely internal segment lineages have their own continuous histories, but with many different external segment partners [<xref ref-type="bibr" rid="B32">32</xref>
]. Here we will consider the patterns of acquisition and loss of the external protein coding segments on the backbone of each internal protein coding segment, since this represents the mechanism for generating antigenic novelty in a circulating population of viruses from the perspective of the internal protein coding segments.</p>
<p>With the aim of exploring the rate of genetic exchange among viral subtypes during evolution, we carried out a phylodynamic analysis of the influenza A viruses circulating in avian populations (in both wild birds and domestic birds) from 1956 – 2011. In particular, we focused on investigating how avian influenza segments specific to the HA, NA and combined HA-NA subtypes evolve on the backbones of the six internal gene segments, by quantifying the inter-subtype reassortment rates and testing possible ecological and evolutionary factors that are associated with such pattern of reassortment: the relative proportion of each host species, the inferred evolutionary relaxed clock rates and the selective constraint acting on the genes.</p>
</sec>
<sec sec-type="results"><title>Results</title>
<sec><title>Time-scaled phylogenies of internal genes of Eurasian AIV</title>
<p>To explore how the HA and NA subtypes reassorted with the internal segments of Eurasian AIV, we used discrete trait models (in BEAST) upon empirical phylogenetic trees of the six internal gene segments. For reasons of computational tractability, we analyzed subsamples of 344 sequences generated using a stratified approach which maintained the range of genetic and subtype diversity. These were composed of 7 HA and 7 NA subtypes and 13 combined HA-NA subtypes respectively (see Methods and Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S2).</p>
<p>The reassortment history of Eurasian AIV was inferred by discrete trait models (HA, NA and combined HA-NA subtypes) on phylogenetic trees. Figure <xref ref-type="fig" rid="F1">1</xref>
shows the MCC tree of PB2 with branches colored by the HA subtype at the child nodes (or tips), and the other trees can be found in Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S2. The two most over-represented subtypes can be easily seen to form single subtype clades; colored green (H5N1) and orange (H9N2). However, it is apparent that some different subtypes were interspersed with each other, and the internal segments did not form monophyletic lineages with respect to antigenic subtype, indicating extensive reassortment between external and internal gene segments. A separation of two main clades can be seen in the PB2 phylogeny (Figure <xref ref-type="fig" rid="F1">1</xref>
A-C). The sequences in each clade cover a range of subtypes, host, sample location and time and do not correspond to the strains belonging to allele A and B in NS (fewer than 40% of the sequences in either lineage are the same). In contrast, lower genetic diversity was observed in the remaining 4 gene segments (PB1, PA, NP and M) (Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S2).</p>
<fig id="F1" position="float"><label>Figure 1</label>
<caption><p><bold>Bayesian maximum clade credibility phylogenies for the internal segments encoding PB2 and NS. A</bold>
: PB2 phylogeny, colored by HA subtype. <bold>B</bold>
: PB2 phylogeny colored by NA subtype. <bold>C</bold>
: PB2 phylogeny colored by combined HA-NA subtype. <bold>D</bold>
: NS phylogeny colored by combined HA-NA subtype. Branches are colored according to the different subtypes of their descendent nodes (C and D share the same key). The full tree with all taxon names is given in Additional file <xref ref-type="supplementary-material" rid="S1">1</xref>
: Figure S2.</p>
</caption>
<graphic xlink:href="1471-2148-14-16-1"></graphic>
</fig>
</sec>
<sec><title>Reassortment rates of discrete subtypes of Eurasian AIV</title>
<p>A range of inter-subtype reassortment rates was found with respect to the internal segments of Eurasian AIV. For both HA and NA subtypes, there were 7 distinct states and thus each discrete trait model contained 42 (asymmetric) transition pairs (e.g. H1 to H3, H3 to H1, N1 to N2, N2 to N1). For the combined HA-NA subtypes, there were 156 different possible transition pairs based on 13 distinct states (e.g. H1N1 to H3N8, H3N8 to H1N1). We analyzed the inter-subtype reassortment rate at three levels: rate per segment (Table <xref ref-type="table" rid="T1">1</xref>
), mean rate per subtype (Table <xref ref-type="table" rid="T2">2</xref>
) and absolute transition rate between pairs of subtypes (Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S2).</p>
<table-wrap position="float" id="T1"><label>Table 1</label>
<caption><p>Overall reassortment rate of six internal segments for H, N and HN combined subtype</p>
</caption>
<table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
</colgroup>
<thead valign="top"><tr><th align="left"> </th>
<th align="center"><bold>PB2</bold>
</th>
<th align="center"><bold>PB1</bold>
</th>
<th align="center"><bold>PA</bold>
</th>
<th align="center"><bold>NP</bold>
</th>
<th align="center"><bold>M</bold>
</th>
<th align="center"><bold>NS</bold>
</th>
</tr>
</thead>
<tbody valign="top"><tr><td align="left" valign="bottom">H<hr></hr>
</td>
<td align="center" valign="bottom">0.30 (0.13,0.39)<hr></hr>
</td>
<td align="center" valign="bottom">0.30 (0.13,0.45)<hr></hr>
</td>
<td align="center" valign="bottom">0.31 (0.11,0.44)<hr></hr>
</td>
<td align="center" valign="bottom">0.28 (0.13,0.38)<hr></hr>
</td>
<td align="center" valign="bottom">0.37 (0.17,0.46)<hr></hr>
</td>
<td align="center" valign="bottom">0.63 (0.28,0.70)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">N<hr></hr>
</td>
<td align="center" valign="bottom">0.46 (0.22,0.57)<hr></hr>
</td>
<td align="center" valign="bottom">0.47 (0.18,0.63)<hr></hr>
</td>
<td align="center" valign="bottom">0.49 (0.19,0.64)<hr></hr>
</td>
<td align="center" valign="bottom">0.29 (0.14,0.35)<hr></hr>
</td>
<td align="center" valign="bottom">0.40 (0.29,0.49)<hr></hr>
</td>
<td align="center" valign="bottom">0.72 (0.43,0.97)<hr></hr>
</td>
</tr>
<tr><td align="left">H-N</td>
<td align="center">0.37 (0.21,0.45)</td>
<td align="center">0.46 (0.26,0.52)</td>
<td align="center">0.37 (0.21,0.45)</td>
<td align="center">0.34 (0.21,0.40)</td>
<td align="center">0.35 (0.20,0.43)</td>
<td align="center">0.68 (0.41,0.76)</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap position="float" id="T2"><label>Table 2</label>
<caption><p>Inter-subtype reassortment rate per subtype for H, N and H-N combined subtypes of 6 internal segments</p>
</caption>
<table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
</colgroup>
<thead valign="top"><tr><th align="left"> </th>
<th align="center"><bold>PB2</bold>
</th>
<th align="center"><bold>PB1</bold>
</th>
<th align="center"><bold>PA</bold>
</th>
<th align="center"><bold>NP</bold>
</th>
<th align="center"><bold>M</bold>
</th>
<th align="center"><bold>NS</bold>
</th>
</tr>
</thead>
<tbody valign="top"><tr><td align="left" valign="bottom"><bold>H</bold>
<hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H1<hr></hr>
</td>
<td align="center" valign="bottom">0.47<hr></hr>
</td>
<td align="center" valign="bottom">0.54<hr></hr>
</td>
<td align="center" valign="bottom">0.39<hr></hr>
</td>
<td align="center" valign="bottom">0.45<hr></hr>
</td>
<td align="center" valign="bottom">0.54<hr></hr>
</td>
<td align="center" valign="bottom">1.15<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H3<hr></hr>
</td>
<td align="center" valign="bottom">0.35<hr></hr>
</td>
<td align="center" valign="bottom">0.35<hr></hr>
</td>
<td align="center" valign="bottom">0.35<hr></hr>
</td>
<td align="center" valign="bottom">0.29<hr></hr>
</td>
<td align="center" valign="bottom">0.41<hr></hr>
</td>
<td align="center" valign="bottom">0.65<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H4<hr></hr>
</td>
<td align="center" valign="bottom">0.44<hr></hr>
</td>
<td align="center" valign="bottom">0.46<hr></hr>
</td>
<td align="center" valign="bottom">0.53<hr></hr>
</td>
<td align="center" valign="bottom">0.41<hr></hr>
</td>
<td align="center" valign="bottom">0.53<hr></hr>
</td>
<td align="center" valign="bottom">0.78<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H5<hr></hr>
</td>
<td align="center" valign="bottom">0.06<hr></hr>
</td>
<td align="center" valign="bottom">0.04<hr></hr>
</td>
<td align="center" valign="bottom">0.05<hr></hr>
</td>
<td align="center" valign="bottom">0.02<hr></hr>
</td>
<td align="center" valign="bottom">0.03<hr></hr>
</td>
<td align="center" valign="bottom">0.06<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H6<hr></hr>
</td>
<td align="center" valign="bottom">0.09<hr></hr>
</td>
<td align="center" valign="bottom">0.10<hr></hr>
</td>
<td align="center" valign="bottom">0.11<hr></hr>
</td>
<td align="center" valign="bottom">0.09<hr></hr>
</td>
<td align="center" valign="bottom">0.12<hr></hr>
</td>
<td align="center" valign="bottom">0.17<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H7<hr></hr>
</td>
<td align="center" valign="bottom">0.16<hr></hr>
</td>
<td align="center" valign="bottom">0.23<hr></hr>
</td>
<td align="center" valign="bottom">0.14<hr></hr>
</td>
<td align="center" valign="bottom">0.18<hr></hr>
</td>
<td align="center" valign="bottom">0.14<hr></hr>
</td>
<td align="center" valign="bottom">0.26<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H9<hr></hr>
</td>
<td align="center" valign="bottom">0.02<hr></hr>
</td>
<td align="center" valign="bottom">0.05<hr></hr>
</td>
<td align="center" valign="bottom">0.04<hr></hr>
</td>
<td align="center" valign="bottom">0.04<hr></hr>
</td>
<td align="center" valign="bottom">0.03<hr></hr>
</td>
<td align="center" valign="bottom">0.04<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom"><bold>N</bold>
<hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">N1<hr></hr>
</td>
<td align="center" valign="bottom">0.06<hr></hr>
</td>
<td align="center" valign="bottom">0.07<hr></hr>
</td>
<td align="center" valign="bottom">0.07<hr></hr>
</td>
<td align="center" valign="bottom">0.04<hr></hr>
</td>
<td align="center" valign="bottom">0.06<hr></hr>
</td>
<td align="center" valign="bottom">0.05<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">N2<hr></hr>
</td>
<td align="center" valign="bottom">0.09<hr></hr>
</td>
<td align="center" valign="bottom">0.07<hr></hr>
</td>
<td align="center" valign="bottom">0.10<hr></hr>
</td>
<td align="center" valign="bottom">0.07<hr></hr>
</td>
<td align="center" valign="bottom">0.12<hr></hr>
</td>
<td align="center" valign="bottom">0.16<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">N3<hr></hr>
</td>
<td align="center" valign="bottom">0.29<hr></hr>
</td>
<td align="center" valign="bottom">0.28<hr></hr>
</td>
<td align="center" valign="bottom">0.29<hr></hr>
</td>
<td align="center" valign="bottom">0.20<hr></hr>
</td>
<td align="center" valign="bottom">0.35<hr></hr>
</td>
<td align="center" valign="bottom">0.59<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">N5<hr></hr>
</td>
<td align="center" valign="bottom">0.58<hr></hr>
</td>
<td align="center" valign="bottom">0.48<hr></hr>
</td>
<td align="center" valign="bottom">0.64<hr></hr>
</td>
<td align="center" valign="bottom">0.36<hr></hr>
</td>
<td align="center" valign="bottom">0.75<hr></hr>
</td>
<td align="center" valign="bottom">1.01<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">N6<hr></hr>
</td>
<td align="center" valign="bottom">0.72<hr></hr>
</td>
<td align="center" valign="bottom">0.58<hr></hr>
</td>
<td align="center" valign="bottom">0.69<hr></hr>
</td>
<td align="center" valign="bottom">0.37<hr></hr>
</td>
<td align="center" valign="bottom">0.76<hr></hr>
</td>
<td align="center" valign="bottom">1.06<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">N7<hr></hr>
</td>
<td align="center" valign="bottom">0.57<hr></hr>
</td>
<td align="center" valign="bottom">0.84<hr></hr>
</td>
<td align="center" valign="bottom">0.59<hr></hr>
</td>
<td align="center" valign="bottom">0.39<hr></hr>
</td>
<td align="center" valign="bottom">0.56<hr></hr>
</td>
<td align="center" valign="bottom">0.99<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">N8<hr></hr>
</td>
<td align="center" valign="bottom">0.51<hr></hr>
</td>
<td align="center" valign="bottom">0.42<hr></hr>
</td>
<td align="center" valign="bottom">0.49<hr></hr>
</td>
<td align="center" valign="bottom">0.28<hr></hr>
</td>
<td align="center" valign="bottom">0.60<hr></hr>
</td>
<td align="center" valign="bottom">0.86<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom"><bold>H-N</bold>
<hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H1N1<hr></hr>
</td>
<td align="center" valign="bottom">0.48<hr></hr>
</td>
<td align="center" valign="bottom">0.60<hr></hr>
</td>
<td align="center" valign="bottom">0.41<hr></hr>
</td>
<td align="center" valign="bottom">0.45<hr></hr>
</td>
<td align="center" valign="bottom">0.46<hr></hr>
</td>
<td align="center" valign="bottom">0.90<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H3N8<hr></hr>
</td>
<td align="center" valign="bottom">0.34<hr></hr>
</td>
<td align="center" valign="bottom">0.39<hr></hr>
</td>
<td align="center" valign="bottom">0.34<hr></hr>
</td>
<td align="center" valign="bottom">0.29<hr></hr>
</td>
<td align="center" valign="bottom">0.33<hr></hr>
</td>
<td align="center" valign="bottom">0.64<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H4N6<hr></hr>
</td>
<td align="center" valign="bottom">0.45<hr></hr>
</td>
<td align="center" valign="bottom">0.54<hr></hr>
</td>
<td align="center" valign="bottom">0.47<hr></hr>
</td>
<td align="center" valign="bottom">0.38<hr></hr>
</td>
<td align="center" valign="bottom">0.43<hr></hr>
</td>
<td align="center" valign="bottom">0.77<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H5N1<hr></hr>
</td>
<td align="center" valign="bottom">0.03<hr></hr>
</td>
<td align="center" valign="bottom">0.05<hr></hr>
</td>
<td align="center" valign="bottom">0.04<hr></hr>
</td>
<td align="center" valign="bottom">0.03<hr></hr>
</td>
<td align="center" valign="bottom">0.03<hr></hr>
</td>
<td align="center" valign="bottom">0.06<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H5N2<hr></hr>
</td>
<td align="center" valign="bottom">0.35<hr></hr>
</td>
<td align="center" valign="bottom">0.46<hr></hr>
</td>
<td align="center" valign="bottom">0.38<hr></hr>
</td>
<td align="center" valign="bottom">0.31<hr></hr>
</td>
<td align="center" valign="bottom">0.32<hr></hr>
</td>
<td align="center" valign="bottom">0.66<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H5N3<hr></hr>
</td>
<td align="center" valign="bottom">0.35<hr></hr>
</td>
<td align="center" valign="bottom">0.54<hr></hr>
</td>
<td align="center" valign="bottom">0.45<hr></hr>
</td>
<td align="center" valign="bottom">0.36<hr></hr>
</td>
<td align="center" valign="bottom">0.42<hr></hr>
</td>
<td align="center" valign="bottom">0.75<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H6N1<hr></hr>
</td>
<td align="center" valign="bottom">0.10<hr></hr>
</td>
<td align="center" valign="bottom">0.16<hr></hr>
</td>
<td align="center" valign="bottom">0.15<hr></hr>
</td>
<td align="center" valign="bottom">0.09<hr></hr>
</td>
<td align="center" valign="bottom">0.12<hr></hr>
</td>
<td align="center" valign="bottom">0.16<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H6N2<hr></hr>
</td>
<td align="center" valign="bottom">0.21<hr></hr>
</td>
<td align="center" valign="bottom">0.29<hr></hr>
</td>
<td align="center" valign="bottom">0.18<hr></hr>
</td>
<td align="center" valign="bottom">0.18<hr></hr>
</td>
<td align="center" valign="bottom">0.18<hr></hr>
</td>
<td align="center" valign="bottom">0.51<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H6N5<hr></hr>
</td>
<td align="center" valign="bottom">0.46<hr></hr>
</td>
<td align="center" valign="bottom">0.44<hr></hr>
</td>
<td align="center" valign="bottom">0.40<hr></hr>
</td>
<td align="center" valign="bottom">0.29<hr></hr>
</td>
<td align="center" valign="bottom">0.42<hr></hr>
</td>
<td align="center" valign="bottom">0.84<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H7N1<hr></hr>
</td>
<td align="center" valign="bottom">0.36<hr></hr>
</td>
<td align="center" valign="bottom">0.41<hr></hr>
</td>
<td align="center" valign="bottom">0.29<hr></hr>
</td>
<td align="center" valign="bottom">0.34<hr></hr>
</td>
<td align="center" valign="bottom">0.29<hr></hr>
</td>
<td align="center" valign="bottom">0.59<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H7N3<hr></hr>
</td>
<td align="center" valign="bottom">0.22<hr></hr>
</td>
<td align="center" valign="bottom">0.16<hr></hr>
</td>
<td align="center" valign="bottom">0.17<hr></hr>
</td>
<td align="center" valign="bottom">0.21<hr></hr>
</td>
<td align="center" valign="bottom">0.16<hr></hr>
</td>
<td align="center" valign="bottom">0.25<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H7N7<hr></hr>
</td>
<td align="center" valign="bottom">0.40<hr></hr>
</td>
<td align="center" valign="bottom">0.60<hr></hr>
</td>
<td align="center" valign="bottom">0.41<hr></hr>
</td>
<td align="center" valign="bottom">0.41<hr></hr>
</td>
<td align="center" valign="bottom">0.33<hr></hr>
</td>
<td align="center" valign="bottom">0.81<hr></hr>
</td>
</tr>
<tr><td align="left">H9N2</td>
<td align="center">0.02</td>
<td align="center">0.05</td>
<td align="center">0.04</td>
<td align="center">0.03</td>
<td align="center">0.03</td>
<td align="center">0.05</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In general, the overall reassortment rate at the segment level was slightly higher for the NA subtypes than for the HA and HA-NA combined subtypes (Table <xref ref-type="table" rid="T1">1</xref>
). In addition, the reassortment rates of the combined HA-NA subtypes (e.g. H5N1) were related to the rate of the corresponding independent HA (e.g. H5) and NA (e.g.N1) subtypes with P < 0.01 in all six segments.</p>
<p>When comparing the reassortment rates among the six internal segments, we found the rates for five (PB2, PB1, PA, NP and M) were similar. Among these NP had a slightly lower and PB1 slightly higher overall rate for HA, NA and the HA-NA combined subtypes. In contrast to these five segments, the NS segment possessed significantly higher overall reassortment rates than those of the other segments (Table <xref ref-type="table" rid="T1">1</xref>
).</p>
<p>We compared the mean reassortment rates of different subtypes, and found a substantial range among each subtype in all six internal segments and all three subtype traits (Table <xref ref-type="table" rid="T2">2</xref>
). Specifically, of the HA subtypes, H5, H6, H7 and H9 showed lower reassortment rates relative to the others, and H5 and H9 had the lowest rates in the six internal segments. By contrast, the reassortment rates of the H1, H4, and H3 subtypes were high, especially for H1 with the NS segment, which had the highest rate (1.148 exchanges per year) compared to all the others. For NA subtypes, N1 and N2 presented lower reassortment rates compared to those of N3, N5, N6, N7 and N9. For the HA-NA joint subtypes, reassortment rates of H5N1 and H9N2 were considerably lower than with the other subtypes in all segments, and H9N2 had the lowest rates in the PB2 segment (0.0220 exchanges per year) and NS segment (0.0487 exchanges per year), while H5N1 had the lowest rates in the other segments (PB1: 0.0449; PA: 0.0352; NP: 0.0249; and M: 0.0266 exchanges per year); other subtypes with low rates were H6N1, H6N2 and H7N3 subtypes. In contrast, H1N1, H3N8 and H4N6 and others exhibited high reassortment rates.</p>
<p>For clarity, the patterns of pairwise transition between individual subtypes are represented as a series of heat maps for the six internal segments in Figure <xref ref-type="fig" rid="F2">2</xref>
(and Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S5). The inferred asymmetric reassortment transition rate matrixes are shown in Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S2. The color intensity in each column reflects the rate with which each subtype (labeled on the y-axis as donors) was found to reassort into a clade of other subtypes. The color intensity in each row reflects the rate with which the reciprocal events were inferred (labeled on x-axis as recipients).</p>
<fig id="F2" position="float"><label>Figure 2</label>
<caption><p><bold>Reassortment rates between internal segments of different subtypes. A</bold>
: PB2 segments associated with HA subtypes shown, which means PB2 that are associated with H1, H3 and H4 subtypes rapidly reassort to become associated with any HA, while PB2 associated with H5, H6, H7 and H9 subtypes reassort more slowly to become associated with any HA. <bold>B</bold>
: PB2 segments associated with NA subtypes shown; <bold>C</bold>
: PB2 segments associated with combined HA- NA subtypes shown; <bold>D</bold>
: NS segments associated with combined HA- NA subtypes shown. Reassortment rates are shown as a heat map. Both forward and reverse rates are shown, with subtypes labeled on the left side acting as donors and subtypes labeled at the bottom acting as recipients. The values of the mean reassortment rates, in exchanges per year, are represented as colors from yellow (low) to red (high). The scale of rates is on the left hand side of each plot.</p>
</caption>
<graphic xlink:href="1471-2148-14-16-2"></graphic>
</fig>
<p>A range of reassortment rates were observed among subtypes as recipients, which indicates that each subtype possesses a different probability of giving the viral gene to viruses of other subtypes during co-infection. In comparison, the rates of subtypes as donors from the other subtypes were consistent which gives rise to the horizontal striping seen in the heat maps. The distinction between source and sink rates were further quantified by estimating the variance of rates in columns to the variance of rates in rows (Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S2, Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S6), which showed the variances were significantly smaller in the source to sink direction than the reverse. This pattern can be seen in the HA, NA subtypes independently and in the HA-NA combined subtypes, and is also similar in all six internal segments.</p>
</sec>
<sec><title>Correlation of reassortment with genetic diversity</title>
<p>To further explore what factors are required to explain the pattern of reassortment rate in the segment level, we first tested the genetic diversity at the segment level which can be represented by the mean time to the most recent common ancestor (tMRCA) based on each phylogeny. The tMRCA of three internal gene segments were quite similar and recent: PB2 (75.71 years, with 95% HPD (Highest Posterior Density): 64.70-87.21), PA (77.95 years with 95% HPD: 64.42-93.01) and the M segment (73.17 years, with 95% HPD of 61.56-87.87). The PB1 segment had an older mean tMRCA, of 105.12 years (95% HPD: 80.21-135.71) while NP was less divergent than the other segments (tMRCA 63.22 years; 95% HPD: 58.32-69.33). Compared to the other segments, the tree topology of the NS segment was characterized by a very deep divergence between the two sub-lineages formed by the A and B alleles which represent two distinct gene pools, and consequently has a much older tMRCA of 226 years (95% HPD: 112–273) (Figure <xref ref-type="fig" rid="F1">1</xref>
D).</p>
<p>The association between genetic diversity and the overall reassortment rate of six internal segments of Eurasian AIV was tested for HA, NA and HA-NA combined subtypes, respectively. The results (Figure <xref ref-type="fig" rid="F3">3</xref>
) exhibited a strong positive correlation for the H subtype (p = 4.7 × 10<sup>-4</sup>
, R = 0.982), the N subtype (p = 0.01, R = 0.879) and the H-N combined subtype (p = 7.569 × 10<sup>-5</sup>
, R = 0.99).</p>
<fig id="F3" position="float"><label>Figure 3</label>
<caption><p><bold>Correlation between TMRCA (years) and reassortment rate (exchanges per year).</bold>
For HA subtype, NA subtype and HA-NA combined subtype, respectively. Scatterplots of mean TMRCA (95% HPD intervals were shown in solid vertical lines) against the mean reassortment rate (95% HPD intervals shown in dotted horizontal lines) of each internal segment. The related segment names were labelled alongside.</p>
</caption>
<graphic xlink:href="1471-2148-14-16-3"></graphic>
</fig>
</sec>
<sec><title>Correlation of reassortment rates with host</title>
<p>The prevalence of AIV in different avian hosts could also influence the reassortment pattern. The host prevalence of different subtypes was calculated from the subsampled data sets used to make the phylogenetic trees, which was proportional to the original data. The primary host difference lies between wild birds and domestic birds which live in completely different ecosystems. Subtypes H1N1, H3N8 and H4N6 had the highest prevalence in wild birds (Table <xref ref-type="table" rid="T3">3</xref>
). We tested the association between individual reassortment rates and the proportion of domestic birds for the HA-NA combined subtypes. We found that there was a significant positive correlation between the proportion of wild bird strains per subtype and the reassortment rates for those subtypes, and that the correlation held in each of the six segments (PB2: p = 0.009, R = 0.69; PB1: p = 0.004, R = 0.73; PA: P = 0.006, R = 0.71; NP: P = 0.02, R = 0.62; M: P = 0.06, R = 0.71; NS: p = 0.002, R = 0.76; Figure <xref ref-type="fig" rid="F4">4</xref>
A and Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S8).</p>
<table-wrap position="float" id="T3"><label>Table 3</label>
<caption><p>Proportion of domestic birds of each subtype analyzed</p>
</caption>
<table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
</colgroup>
<thead valign="top"><tr><th align="left"><bold>Subtype</bold>
</th>
<th align="center"><bold>Wild</bold>
<sup><bold>a</bold>
</sup>
</th>
<th align="center"><bold>Ans</bold>
<sup><bold>b</bold>
</sup>
</th>
<th align="center"><bold>Ans-domestic</bold>
<sup><bold>c</bold>
</sup>
</th>
<th align="center"><bold>Gal-domestic</bold>
<sup><bold>d</bold>
</sup>
</th>
</tr>
</thead>
<tbody valign="top"><tr><td align="left" valign="bottom">H1N1<hr></hr>
</td>
<td align="center" valign="bottom">0.50<hr></hr>
</td>
<td align="center" valign="bottom">1.00<hr></hr>
</td>
<td align="center" valign="bottom">0.50<hr></hr>
</td>
<td align="center" valign="bottom">0.00<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H3N8<hr></hr>
</td>
<td align="center" valign="bottom">0.44<hr></hr>
</td>
<td align="center" valign="bottom">0.89<hr></hr>
</td>
<td align="center" valign="bottom">0.50<hr></hr>
</td>
<td align="center" valign="bottom">0.06<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H4N6<hr></hr>
</td>
<td align="center" valign="bottom">0.58<hr></hr>
</td>
<td align="center" valign="bottom">0.92<hr></hr>
</td>
<td align="center" valign="bottom">0.42<hr></hr>
</td>
<td align="center" valign="bottom">0.00<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H5N1<hr></hr>
</td>
<td align="center" valign="bottom">0.35<hr></hr>
</td>
<td align="center" valign="bottom">0.42<hr></hr>
</td>
<td align="center" valign="bottom">0.25<hr></hr>
</td>
<td align="center" valign="bottom">0.40<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H5N2<hr></hr>
</td>
<td align="center" valign="bottom">0.38<hr></hr>
</td>
<td align="center" valign="bottom">0.81<hr></hr>
</td>
<td align="center" valign="bottom">0.38<hr></hr>
</td>
<td align="center" valign="bottom">0.19<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H5N3<hr></hr>
</td>
<td align="center" valign="bottom">0.42<hr></hr>
</td>
<td align="center" valign="bottom">0.92<hr></hr>
</td>
<td align="center" valign="bottom">0.58<hr></hr>
</td>
<td align="center" valign="bottom">0.00<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H6N1<hr></hr>
</td>
<td align="center" valign="bottom">0.34<hr></hr>
</td>
<td align="center" valign="bottom">0.45<hr></hr>
</td>
<td align="center" valign="bottom">0.24<hr></hr>
</td>
<td align="center" valign="bottom">0.41<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H6N2<hr></hr>
</td>
<td align="center" valign="bottom">0.50<hr></hr>
</td>
<td align="center" valign="bottom">0.88<hr></hr>
</td>
<td align="center" valign="bottom">0.42<hr></hr>
</td>
<td align="center" valign="bottom">0.08<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H6N5<hr></hr>
</td>
<td align="center" valign="bottom">0.50<hr></hr>
</td>
<td align="center" valign="bottom">1.00<hr></hr>
</td>
<td align="center" valign="bottom">0.50<hr></hr>
</td>
<td align="center" valign="bottom">0.00<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H7N1<hr></hr>
</td>
<td align="center" valign="bottom">0.30<hr></hr>
</td>
<td align="center" valign="bottom">0.60<hr></hr>
</td>
<td align="center" valign="bottom">0.40<hr></hr>
</td>
<td align="center" valign="bottom">0.30<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H7N3<hr></hr>
</td>
<td align="center" valign="bottom">0.20<hr></hr>
</td>
<td align="center" valign="bottom">0.40<hr></hr>
</td>
<td align="center" valign="bottom">0.20<hr></hr>
</td>
<td align="center" valign="bottom">0.60<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H7N7<hr></hr>
</td>
<td align="center" valign="bottom">0.50<hr></hr>
</td>
<td align="center" valign="bottom">0.75<hr></hr>
</td>
<td align="center" valign="bottom">0.25<hr></hr>
</td>
<td align="center" valign="bottom">0.25<hr></hr>
</td>
</tr>
<tr><td align="left">H9N2</td>
<td align="center">0.15</td>
<td align="center">0.26</td>
<td align="center">0.18</td>
<td align="center">0.67</td>
</tr>
</tbody>
</table>
<table-wrap-foot><p><sup>a</sup>
The overall proportion of AIV with a wild bird host.</p>
<p><sup>b</sup>
The overall proportion of AIV with a host species belonging to Anseriformes.</p>
<p><sup>c</sup>
The proportion of AIV with a domestic Anseriformes host.</p>
<p><sup>d</sup>
The proportion of AIV with a domestic Galliformes host.</p>
</table-wrap-foot>
</table-wrap>
<fig id="F4" position="float"><label>Figure 4</label>
<caption><p><bold>Correlation between prevalence in host and reassortment rate per subtype in PB2 segment.</bold>
The scatterplot of the mean reassortment rates of 13H-N combined subtypes of each internal segment against the proportion of AIV 1) in wild birds, 2) in Anseriformes 3) in domestic birds of Anseriformes and 4) in domestic birds of Galliformes. The plots of the remaining five segments are shown in supplementary figures (Additional file <xref ref-type="supplementary-material" rid="S1">1</xref>
: Figure S8-S11).</p>
</caption>
<graphic xlink:href="1471-2148-14-16-4"></graphic>
</fig>
<p>The prevalence of AIV in different host species (by Order) was also analyzed. In this dataset, subtypes with a high reassortment rate also have the highest prevalence in Anseriformes, which is the Order representing birds living in predominantly freshwater habitats (Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S4). We found a very strong positive correlation between reassortment rates and the proportion of Anseriformes (PB2: p = 3 × 10<sup>-4</sup>
, R = 0.85; PB1: p = 2 × 10<sup>-4</sup>
, R = 0.84; PA: p = 2 × 10<sup>-4</sup>
, R = 0.84; NP: p = 0.002, R = 0.77; M: p = 4 × 10<sup>-5</sup>
, R = 0.89; NS: p = 9 × 10<sup>-6</sup>
, R = 0.90) (Figure <xref ref-type="fig" rid="F4">4</xref>
B and Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S9), which may indicate the subtype with high reassortment rate could be better explained by their high proportion in waterfowl than wild birds. We further tested wild (mainly mallard in our study) and domestic Anseriformes (mainly domestic duck in our study) separately (Table <xref ref-type="table" rid="T3">3</xref>
). Both of them show a strong positive correlation with reassortment rate per subtype (Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S10 and Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S11).</p>
<p>On the other hand, subtypes H9N2, H7N3, H7N1 and H5N1 had a high prevalence in domestic birds of the order Galliformes (chickens, turkeys, etc.). There was a strong negative correlation between proportion of AIV in domestic birds of the order Galliformes and rate of reassortment per subtype (PB2: p = 0.002, R = −0.77; PB1: p = 0.001, R = −0.79; PA: p = 0.001, R = −0.70; NP: p = 0.008, R = −0.70; M: p = 2 × 10<sup>-4</sup>
, R = −0.84; NS: p = 1 × 10<sup>-4</sup>
, R = − 0.85) (Figure <xref ref-type="fig" rid="F4">4</xref>
C and Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S12).</p>
</sec>
<sec><title>Correlation of reassortment rate with evolutionary rate</title>
<p>We further tested whether the reassortment rates were correlated with evolutionary rates of AIV at both segment and subtype level. The overall evolutionary rates of Eurasian AIV in different segments were the relaxed clock rates per segment from the MCMC tree (Table <xref ref-type="table" rid="T4">4</xref>
) [<xref ref-type="bibr" rid="B25">25</xref>
]. No significant difference in evolutionary relaxed clock rate was found among six internal genes of the analyzed AIV sequences, with estimated mean rates from 2.27 to 2.85 × 10<sup>-3</sup>
substitutions per site per year (subs/site/year). In addition, no significant correlation was identified between overall reassortment per segment and evolutionary rate per segment (P = 0.7).</p>
<table-wrap position="float" id="T4"><label>Table 4</label>
<caption><p>Evolutionary Rate of internal segment of Eurasian AIV</p>
</caption>
<table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
</colgroup>
<thead valign="top"><tr><th align="left"> </th>
<th align="center"><bold>PB2</bold>
</th>
<th align="center"><bold>PB1</bold>
</th>
<th align="center"><bold>PA</bold>
</th>
<th align="center"><bold>NP</bold>
</th>
<th align="center"><bold>M</bold>
</th>
<th align="center"><bold>NS</bold>
</th>
</tr>
</thead>
<tbody valign="top"><tr><td align="left" valign="bottom">segment<sup>a</sup>
<hr></hr>
</td>
<td align="center" valign="bottom">2.85 (2.54,3.19)<hr></hr>
</td>
<td align="center" valign="bottom">2.27 (1.96,2.61)<hr></hr>
</td>
<td align="center" valign="bottom">2.62 (2.32,2.91)<hr></hr>
</td>
<td align="center" valign="bottom">2.65 (2.36,2.95)<hr></hr>
</td>
<td align="center" valign="bottom">2.41 (2.07,2.76)<hr></hr>
</td>
<td align="center" valign="bottom">2.56 (2.18,2.95)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">subtype<sup>b</sup>
<hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H1N1<hr></hr>
</td>
<td align="center" valign="bottom">2.63 (1.88,3.38)<hr></hr>
</td>
<td align="center" valign="bottom">2.36 (1.62,3.10)<hr></hr>
</td>
<td align="center" valign="bottom">2.52 (1.41,3.64)<hr></hr>
</td>
<td align="center" valign="bottom">2.70 (1.69,3.72)<hr></hr>
</td>
<td align="center" valign="bottom">2.62 (1.32,3.93)<hr></hr>
</td>
<td align="center" valign="bottom">2.28 (1.46,3.11)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H3N8<hr></hr>
</td>
<td align="center" valign="bottom">2.40 (1.87,2.93)<hr></hr>
</td>
<td align="center" valign="bottom">2.86 (2.12,3.60)<hr></hr>
</td>
<td align="center" valign="bottom">2.38 (1.62,3.13)<hr></hr>
</td>
<td align="center" valign="bottom">2.30 (1.61,2.99)<hr></hr>
</td>
<td align="center" valign="bottom">2.13 (1.33,2.92)<hr></hr>
</td>
<td align="center" valign="bottom">2.20 (1.60,2.81)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H4N6<hr></hr>
</td>
<td align="center" valign="bottom">2.80 (2.42,3.18)<hr></hr>
</td>
<td align="center" valign="bottom">2.59 (1.85,3.33)<hr></hr>
</td>
<td align="center" valign="bottom">2.33 (1.50,3.15)<hr></hr>
</td>
<td align="center" valign="bottom">2.62 (1.73,3.51)<hr></hr>
</td>
<td align="center" valign="bottom">2.28 (1.39,3.17)<hr></hr>
</td>
<td align="center" valign="bottom">2.56 (1.80,3.31)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H5N1<hr></hr>
</td>
<td align="center" valign="bottom">2.93 (2.61,3.25)<hr></hr>
</td>
<td align="center" valign="bottom">2.84 (2.10,3.58)<hr></hr>
</td>
<td align="center" valign="bottom">2.62 (1.85,3.38)<hr></hr>
</td>
<td align="center" valign="bottom">2.60 (1.94,3.26)<hr></hr>
</td>
<td align="center" valign="bottom">2.32 (1.60,3.04)<hr></hr>
</td>
<td align="center" valign="bottom">2.73 (2.30,3.16)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H5N2<hr></hr>
</td>
<td align="center" valign="bottom">2.91 (2.28,3.54)<hr></hr>
</td>
<td align="center" valign="bottom">2.78 (2.04,3.52)<hr></hr>
</td>
<td align="center" valign="bottom">2.39 (1.57,3.21)<hr></hr>
</td>
<td align="center" valign="bottom">2.74 (1.86,3.62)<hr></hr>
</td>
<td align="center" valign="bottom">2.50 (1.57,3.43)<hr></hr>
</td>
<td align="center" valign="bottom">2.45 (1.83,3.07)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H5N3<hr></hr>
</td>
<td align="center" valign="bottom">3.01 (2.29,3.72)<hr></hr>
</td>
<td align="center" valign="bottom">2.44 (1.70,3.18)<hr></hr>
</td>
<td align="center" valign="bottom">2.42 (1.54,3.30)<hr></hr>
</td>
<td align="center" valign="bottom">2.54 (1.66,3.43)<hr></hr>
</td>
<td align="center" valign="bottom">2.68 (1.51,3.85)<hr></hr>
</td>
<td align="center" valign="bottom">2.25 (1.61,2.88)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H6N1<hr></hr>
</td>
<td align="center" valign="bottom">3.20 (2.67,3.74)<hr></hr>
</td>
<td align="center" valign="bottom">2.78 (2.04,3.52)<hr></hr>
</td>
<td align="center" valign="bottom">2.54 (1.71,3.36)<hr></hr>
</td>
<td align="center" valign="bottom">2.53 (1.81,3.24)<hr></hr>
</td>
<td align="center" valign="bottom">2.16 (1.40,2.93)<hr></hr>
</td>
<td align="center" valign="bottom">2.61 (2.03,3.19)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H6N2<hr></hr>
</td>
<td align="center" valign="bottom">2.85 (2.32,3.38)<hr></hr>
</td>
<td align="center" valign="bottom">2.90 (2.16,3.64)<hr></hr>
</td>
<td align="center" valign="bottom">2.61 (1.78,3.45)<hr></hr>
</td>
<td align="center" valign="bottom">2.52 (1.79,3.25)<hr></hr>
</td>
<td align="center" valign="bottom">2.51 (1.59,3.43)<hr></hr>
</td>
<td align="center" valign="bottom">2.51 (1.88,3.14)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H6N5<hr></hr>
</td>
<td align="center" valign="bottom">2.72 (1.97,3.47)<hr></hr>
</td>
<td align="center" valign="bottom">2.99 (2.25,3.73)<hr></hr>
</td>
<td align="center" valign="bottom">2.45 (1.55,3.34)<hr></hr>
</td>
<td align="center" valign="bottom">2.54 (1.68,3.41)<hr></hr>
</td>
<td align="center" valign="bottom">2.18 (1.16,3.20)<hr></hr>
</td>
<td align="center" valign="bottom">2.30 (1.47,3.13)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H7N1<hr></hr>
</td>
<td align="center" valign="bottom">2.77 (2.08,3.45)<hr></hr>
</td>
<td align="center" valign="bottom">2.68 (1.94,3.42)<hr></hr>
</td>
<td align="center" valign="bottom">2.35 (1.47,3.24)<hr></hr>
</td>
<td align="center" valign="bottom">2.57 (1.71,3.44)<hr></hr>
</td>
<td align="center" valign="bottom">2.37 (1.32,3.41)<hr></hr>
</td>
<td align="center" valign="bottom">2.33 (1.57,3.09)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H7N3<hr></hr>
</td>
<td align="center" valign="bottom">2.54 (1.89,3.18)<hr></hr>
</td>
<td align="center" valign="bottom">2.42 (1.68,3.16)<hr></hr>
</td>
<td align="center" valign="bottom">2.08 (1.31,2.85)<hr></hr>
</td>
<td align="center" valign="bottom">2.47 (1.66,3.28)<hr></hr>
</td>
<td align="center" valign="bottom">2.17 (1.19,3.16)<hr></hr>
</td>
<td align="center" valign="bottom">1.92 (1.25,2.59)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H7N7<hr></hr>
</td>
<td align="center" valign="bottom">2.98 (2.15,3.80)<hr></hr>
</td>
<td align="center" valign="bottom">2.55 (1.81,3.29)<hr></hr>
</td>
<td align="center" valign="bottom">2.31 (1.39,3.22)<hr></hr>
</td>
<td align="center" valign="bottom">2.46 (1.59,3.32)<hr></hr>
</td>
<td align="center" valign="bottom">2.06 (1.07,3.05)<hr></hr>
</td>
<td align="center" valign="bottom">2.32 (1.35,3.28)<hr></hr>
</td>
</tr>
<tr><td align="left">H9N2</td>
<td align="center">3.11 (2.71,3.52)</td>
<td align="center">2.87 (2.13,3.61)</td>
<td align="center">2.74 (1.92,3.56)</td>
<td align="center">2.72 (2.00,3.45)</td>
<td align="center">2.36 (1.61,3.11)</td>
<td align="center">2.62 (2.18,3.05)</td>
</tr>
</tbody>
</table>
<table-wrap-foot><p><sup>a</sup>
Estimated overall evolutionary rate (×10<sup>3</sup>
subst/site/year) per internal segment. Mean rate under the uncorrelated log-normal relaxed molecular clock with 95% HPD of all MCMC samples are shown in brackets.</p>
<p><sup>b</sup>
Estimated evolutionary rate (×10<sup>3</sup>
subst/site/year) per H-N combined subtype in each internal segment. Mean clock rate of certain subtype on MCMC tree nodes with 95% confidential interval (CI) of all MCMC samples are shown in brackets.</p>
</table-wrap-foot>
</table-wrap>
<p>The evolutionary rates in different subtypes were estimated by considering the relaxed clock rates corresponding to branches of particular subtypes. The rates per subtype were also similar between different subtypes in the same segment and not significantly different across different segments (Table <xref ref-type="table" rid="T4">4</xref>
). There was no significant correlation between evolutionary rate and reassortment rate per subtype in each internal segment (P > 0.1) (Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S13), which is also confirmed by the distribution of correlation coefficients between the evolutionary rates (relaxed clock rates from 1000 random MCMC samples post-burnin) and reassortment rates (rates were taken from 1000 random MCMC samples post-burnin). In addition, there was no indication of any interaction between evolutionary rate and the proportion of domestic birds (P > 0.5).</p>
</sec>
<sec><title>Correlation of reassortment with selection</title>
<p>Finally, we estimated the total d<sub>N</sub>
/d<sub>S</sub>
ratio for each internal segment as well as the individual ratio per subtype in each segment (Table <xref ref-type="table" rid="T5">5</xref>
). Initial analysis of the coding regions in the five segments PB2, PB1, PA, NP and M1 showed the d<sub>N</sub>
/d<sub>S</sub>
ratios were similar to each other with mean values from 0.11 to 0.17, indicating, as expected, that at the amino acid level purifying selection predominates in these internal segments; however, the ratio for the NS segment (mean 0.439) was significantly higher than those of the other segments. This suggested a strong correlation between reassortment rate and d<sub>N</sub>
/d<sub>S</sub>
ratio at the segment level (P = 0.006, R = 0.93). However, this result is entirely dependent on the high value for NS1 and the correlation became negative and non-significant when NS1 was left out (R = −0.23, p = 0.83). Compared to the evolutionary rate, the d<sub>N</sub>
/d<sub>S</sub>
ratio for each subtype differs noticeably between different subtypes in the same segment and is very different across different segments (Table <xref ref-type="table" rid="T4">4</xref>
). For instance, in PB2, H4N6 had the lowest d<sub>N</sub>
/d<sub>S</sub>
ratio (0.22, 95% HPD: 0.13-0.30) and H6N1 had the highest (0.76, 95% HPD: 0.55-0.96). In NS1, H7N7 had the lowest ratio (0.351, with 95% HPD: 0–0.842) and H7N3 had the highest ratio which was equal to 1, but this difference was not significant given a wide range of statistical uncertainty for these subtypes in the NS1 segment. We found strong negative correlations between mean d<sub>N</sub>
/d<sub>S</sub>
and the reassortment rate per subtype in each of the six segments (PB2: p = 0.01, R = −0.67; PB1: p = 0.009, R = −0.69; PA: P = 2 × 10<sup>-4</sup>
, R = −0.86; NP: P = 0.01, R = −0.67; M: P = 0.002, R = −0.78; NS: p = 0.001, R = −0.80) (Figure <xref ref-type="fig" rid="F5">5</xref>
).</p>
<table-wrap position="float" id="T5"><label>Table 5</label>
<caption><p><bold>d</bold>
<sub><bold>N</bold>
</sub>
<bold>/d</bold>
<sub><bold>S </bold>
</sub>
<bold>ratio of internal segments of Eurasian AIV</bold>
</p>
</caption>
<table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
</colgroup>
<thead valign="top"><tr><th align="left"> </th>
<th align="center"><bold>PB2</bold>
</th>
<th align="center"><bold>PB1</bold>
</th>
<th align="center"><bold>PA</bold>
</th>
<th align="center"><bold>NP</bold>
</th>
<th align="center"><bold>M1</bold>
</th>
<th align="center"><bold>NS1</bold>
</th>
</tr>
</thead>
<tbody valign="top"><tr><td align="left" valign="bottom">segment<sup>a</sup>
<hr></hr>
</td>
<td align="center" valign="bottom">0.14 (0.14,0.15)<hr></hr>
</td>
<td align="center" valign="bottom">0.12 (0.12,0.13)<hr></hr>
</td>
<td align="center" valign="bottom">0.17 (0.17,0.17)<hr></hr>
</td>
<td align="center" valign="bottom">0.13 (0.13,0.14)<hr></hr>
</td>
<td align="center" valign="bottom">0.11 (0.11,0.12)<hr></hr>
</td>
<td align="center" valign="bottom">0.59 (0.15,1.03)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">subtype<sup>b</sup>
<hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H1N1<hr></hr>
</td>
<td align="center" valign="bottom">0.31 (0.15,0.46)<hr></hr>
</td>
<td align="center" valign="bottom">0.23 (0.09,0.37)<hr></hr>
</td>
<td align="center" valign="bottom">0.30 (0.16,0.43)<hr></hr>
</td>
<td align="center" valign="bottom">0.64 (0.11,0.38)<hr></hr>
</td>
<td align="center" valign="bottom">0.02 (0.02,0.11)<hr></hr>
</td>
<td align="center" valign="bottom">0.49 (0.06,0.91)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H3N8<hr></hr>
</td>
<td align="center" valign="bottom">0.31 (0.19,0.43)<hr></hr>
</td>
<td align="center" valign="bottom">0.25 (0.16,0.34)<hr></hr>
</td>
<td align="center" valign="bottom">0.31 (0.21,0.42)<hr></hr>
</td>
<td align="center" valign="bottom">0.35 (0.14,0.24)<hr></hr>
</td>
<td align="center" valign="bottom">0.15 (0.03,0.28)<hr></hr>
</td>
<td align="center" valign="bottom">0.77 (0.39,1.15)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H4N6<hr></hr>
</td>
<td align="center" valign="bottom">0.22 (0.13,0.30)<hr></hr>
</td>
<td align="center" valign="bottom">0.41 (0.26,0.55)<hr></hr>
</td>
<td align="center" valign="bottom">0.31 (0.19,0.43)<hr></hr>
</td>
<td align="center" valign="bottom">0.47 (0.16,0.31)<hr></hr>
</td>
<td align="center" valign="bottom">0.34 (0.06,0.62)<hr></hr>
</td>
<td align="center" valign="bottom">0.46 (0.15,0.78)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H5N1<hr></hr>
</td>
<td align="center" valign="bottom">0.56 (0.48,0.64)<hr></hr>
</td>
<td align="center" valign="bottom">0.56 (0.48,0.65)<hr></hr>
</td>
<td align="center" valign="bottom">0.92 (0.80,1.04)<hr></hr>
</td>
<td align="center" valign="bottom">0.80 (0.57,0.69)<hr></hr>
</td>
<td align="center" valign="bottom">0.65 (0.48,0.82)<hr></hr>
</td>
<td align="center" valign="bottom">0.97 (0.78,1.15)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H5N2<hr></hr>
</td>
<td align="center" valign="bottom">0.42 (0.29,0.54)<hr></hr>
</td>
<td align="center" valign="bottom">0.29 (0.20,0.37)<hr></hr>
</td>
<td align="center" valign="bottom">0.35 (0.24,0.46)<hr></hr>
</td>
<td align="center" valign="bottom">0.49 (0.24,0.36)<hr></hr>
</td>
<td align="center" valign="bottom">0.23 (0.09,0.37)<hr></hr>
</td>
<td align="center" valign="bottom">0.41 (0.16,0.65)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H5N3<hr></hr>
</td>
<td align="center" valign="bottom">0.36 (0.22,0.49)<hr></hr>
</td>
<td align="center" valign="bottom">0.34 (0.20,0.48)<hr></hr>
</td>
<td align="center" valign="bottom">0.22 (0.13,0.31)<hr></hr>
</td>
<td align="center" valign="bottom">0.51 (0.19,0.35)<hr></hr>
</td>
<td align="center" valign="bottom">0.15 (0.02,0.29)<hr></hr>
</td>
<td align="center" valign="bottom">0.50 (0.15,0.84)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H6N1<hr></hr>
</td>
<td align="center" valign="bottom">0.76 (0.55,0.96)<hr></hr>
</td>
<td align="center" valign="bottom">0.60 (0.45,0.75)<hr></hr>
</td>
<td align="center" valign="bottom">0.57 (0.39,0.74)<hr></hr>
</td>
<td align="center" valign="bottom">0.64 (0.36,0.50)<hr></hr>
</td>
<td align="center" valign="bottom">0.73 (0.38,1.08)<hr></hr>
</td>
<td align="center" valign="bottom">0.83 (0.56,1.11)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H6N2<hr></hr>
</td>
<td align="center" valign="bottom">0.42 (0.29,0.54)<hr></hr>
</td>
<td align="center" valign="bottom">0.27 (0.20,0.35)<hr></hr>
</td>
<td align="center" valign="bottom">0.33 (0.24,0.41)<hr></hr>
</td>
<td align="center" valign="bottom">0.43 (0.21,0.32)<hr></hr>
</td>
<td align="center" valign="bottom">0.50 (0.25,0.75)<hr></hr>
</td>
<td align="center" valign="bottom">0.56 (0.32,0.80)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H6N5<hr></hr>
</td>
<td align="center" valign="bottom">0.53 (0.14,0.92)<hr></hr>
</td>
<td align="center" valign="bottom">0.46 (0.22,0.69)<hr></hr>
</td>
<td align="center" valign="bottom">0.41 (0.20,0.62)<hr></hr>
</td>
<td align="center" valign="bottom">0.77 (0.15,0.46)<hr></hr>
</td>
<td align="center" valign="bottom">0.04 (0.04,0.25)<hr></hr>
</td>
<td align="center" valign="bottom">0.72 (0.12,1.32)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H7N1<hr></hr>
</td>
<td align="center" valign="bottom">0.40 (0.21,0.59)<hr></hr>
</td>
<td align="center" valign="bottom">0.43 (0.19,0.67)<hr></hr>
</td>
<td align="center" valign="bottom">0.37 (0.18,0.56)<hr></hr>
</td>
<td align="center" valign="bottom">0.98 (0.24,0.61)<hr></hr>
</td>
<td align="center" valign="bottom">0.53 (0.08,1.14)<hr></hr>
</td>
<td align="center" valign="bottom">0.76 (0.12,1.40)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H7N3<hr></hr>
</td>
<td align="center" valign="bottom">0.73 (0.08,1.39)<hr></hr>
</td>
<td align="center" valign="bottom">0.57 (0.33,0.80)<hr></hr>
</td>
<td align="center" valign="bottom">0.57 (0.28,0.86)<hr></hr>
</td>
<td align="center" valign="bottom">0.95 (0.17,0.56)<hr></hr>
</td>
<td align="center" valign="bottom">1.05 (0.21,1.89)<hr></hr>
</td>
<td align="center" valign="bottom">1.01 (0.07,2.43)<hr></hr>
</td>
</tr>
<tr><td align="left" valign="bottom">H7N7<hr></hr>
</td>
<td align="center" valign="bottom">0.40 (0.19,0.61)<hr></hr>
</td>
<td align="center" valign="bottom">0.45 (0.20,0.71)<hr></hr>
</td>
<td align="center" valign="bottom">0.29 (0.16,0.42)<hr></hr>
</td>
<td align="center" valign="bottom">0.52 (0.08,0.30)<hr></hr>
</td>
<td align="center" valign="bottom">0.03 (0.05,0.21)<hr></hr>
</td>
<td align="center" valign="bottom">0.35 (0.03,0.84)<hr></hr>
</td>
</tr>
<tr><td align="left">H9N2</td>
<td align="center">0.55 (0.46,0.64)</td>
<td align="center">0.59 (0.52,0.67)</td>
<td align="center">0.64 (0.56,0.72)</td>
<td align="center">0.80 (0.56,0.68)</td>
<td align="center">0.57 (0.40,0.74)</td>
<td align="center">0.92 (0.72,1.12)</td>
</tr>
</tbody>
</table>
<table-wrap-foot><p><sup>a</sup>
Estimated overall d<sub>N</sub>
/d<sub>S</sub>
ratio per internal segment. Mean ratio of each internal segments with 95% (CI) of all MCMC samples shown in brackets.</p>
<p><sup>b</sup>
Estimated d<sub>N</sub>
/d<sub>S</sub>
ratio per H-N combined subtype in each internal segment. Mean d<sub>N</sub>
/d<sub>S</sub>
of each subtype on MCMC tree nodes with 95% CI of all MCMC samples are shown in brackets.</p>
</table-wrap-foot>
</table-wrap>
<fig id="F5" position="float"><label>Figure 5</label>
<caption><p><bold>Correlation between selection and reassortment rate.</bold>
The scatterplot of reassortment rate of HA-NA combined subtype against dN/dS ratio of 6 internal segments were represented by different colors: PB2 (red); PB1 (green); PA (blue); NP (orange); M1 (purple); NS1 (black).</p>
</caption>
<graphic xlink:href="1471-2148-14-16-5"></graphic>
</fig>
<p>In order to confirm the negative correlation between the d<sub>N</sub>
/d<sub>S</sub>
ratio and inter-subtype reassortment rate seen in BEAST, we used the independent approach of Single Likelihood Ancestor Counting (SLAC) in the HyPhy package (<ext-link ext-link-type="uri" xlink:href="http://www.datamonkey.org/">http://www.datamonkey.org/</ext-link>
) to estimate d<sub>N</sub>
/d<sub>S</sub>
per subtype. We found a strong negative correlation between inter-subtype reassortment and d<sub>N</sub>
/d<sub>S</sub>
ratio in all six internal segments using SLAC as well (Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S14). SLAC estimates the numbers of codons subject to positive and negative selection and we were thus able to identify the nature of the differences in selection. There was no evidence for positive selection contributing to the variation among subtypes, however there was substantial variation in the number of negatively selected sites, which was strongly negatively correlated with reassortment rate with R < −0.65, p ≤ 0.01, for all segments except NS1 for which R = −0.55, p = 0.05 (Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S4). Therefore the intensity of purifying selection, or selective constraint, on the internal segments is strongly associated with the rate of reassortment.</p>
<p>In addition, we also found a strong negative correlation between the proportion of wild birds and the d<sub>N</sub>
/d<sub>S</sub>
ratio of different subtypes per segment (PB2: p = 0.04, R = −0.58; PB1: p = 0.01, R = −0.46; PA: p = 0.04, R = −0.57; NP: p = 0.02, R = −0.64; M: p = 0.007, R = −0.71; NS: p = 0.002, R = −0.77). We then jointly estimated the partial correlation coefficients among d<sub>N</sub>
/d<sub>S</sub>
, host species and reassortment rates [<xref ref-type="bibr" rid="B33">33</xref>
]. This showed that reassortment rate was still significantly negatively correlated with d<sub>N</sub>
/d<sub>S</sub>
independent of the proportion of wild birds as host in all six internal segments (R: -0.79 to −0.44, P < 0.01; Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S3 column 5 and 7); and independently positively correlated with proportion of wild birds as host (R:-0.66 to −0.29, P < 0.07; Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S3 column 8 and 10) with borderline significance, thus both host and selective constraint are associated with reassortment rate.</p>
</sec>
</sec>
<sec sec-type="discussion"><title>Discussion</title>
<p>We have systematically analyzed the evolution of internal segments of avian influenza viruses use a comprehensive avian influenza A virus dataset. Using Bayesian analysis of phylogenetic discrete traits within the integrated statistical framework in the BEAST software package, we quantified evolutionary parameters including reassortment rate, evolutionary relaxed clock rate and selective constraint for discrete. Notably, the quantified pattern of reassortment illustrates how and to what extent the subtypes of AIV interact in the internal segments in bird populations. Firstly, the overall inter-subtype reassortment rates are similar in five segments, while NS has higher rates overall, and secondly, in each segment, reassortment rates of different subtypes showed significant diversity and were asymmetric.</p>
<p>The reassortment rate of the six internal segments of Eurasian AIV is significantly associated with gene diversity, which is probably due to the phylogenetic history and the biological function of the segments during virus infection [<xref ref-type="bibr" rid="B34">34</xref>
]. Here we identified that NS has a much more diverse phylogeny and higher inter-subtype reassortment rate than the other five segments of the avian influenza viruses of Eurasian lineage. This result is in agreement with those of a previous study on AIV from the wild bird population of North America, which found extensive phylogenetic incongruence in the NS gene [<xref ref-type="bibr" rid="B32">32</xref>
]. The NS phylogeny of avian influenza viruses is characterized by a deep divergence between the A and B alleles suggesting that the two alleles are subject to some form of balancing selection [<xref ref-type="bibr" rid="B32">32</xref>
,<xref ref-type="bibr" rid="B35">35</xref>
-<xref ref-type="bibr" rid="B38">38</xref>
].</p>
<p>While we found that reassortment rate was not associated with evolutionary rate, the segment-wide d<sub>N</sub>
/d<sub>S</sub>
ratio per subtype was strongly negatively correlated with the inter-subtype reassortment. The variation in d<sub>N</sub>
/d<sub>S</sub>
ratio was revealed by two independent methods: robust counting in BEAST [<xref ref-type="bibr" rid="B39">39</xref>
], and SLAC in HyPhy [<xref ref-type="bibr" rid="B40">40</xref>
]. The latter gives direct estimates of the numbers of selected sites from which it was seen that the variation between subtypes in d<sub>N</sub>
/d<sub>S</sub>
appears to arise from differences in the extent of purifying selection, as indicated by differences in the number of negatively selected sites for each subtype, which itself yielded a strongly negative association with reassortment rate (SLAC identified only 3 positively selected sites; Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S4). This suggests that the overall level of purifying selection, or selective constraint, differs between subtypes: more negatively selected sites indicating higher overall constraints leading to lower d<sub>N</sub>
/d<sub>S</sub>
.</p>
<p>As mentioned above, the evolutionary relaxed clock rates (of the internal gene segments) in Eurasian AIV show little difference among viral subtypes and also between the internal segments. In our results therefore, evolutionary rates were not correlated with inter-subtype reassortment rates and also not affected much by host species. This agrees with the conclusions of a previous study which indicated that the evolutionary rates in wild aquatic birds are not radically lower than those seen in poultry species [<xref ref-type="bibr" rid="B41">41</xref>
,<xref ref-type="bibr" rid="B42">42</xref>
]. Pathogenicity status might also be an explanation of the lower reassortment rates of H5N1: HPAI H5N1 is lethal to poultry [<xref ref-type="bibr" rid="B43">43</xref>
] thus could have had less chance to undergo reassortment in the host due to the short duration of infection of the birds as compared to the less symptomatic infection, typically of longer duration, in low pathogenicity strains. However, our results indicate that pathogenicity status (which is defined in chickens) is not the major factor affecting the overall reassortment rate in all species; of those subtypes with high pathogenicity, such as H5N1, H7N3, H7N1 and H7N7, only H5N1 has low reassortment rate; and there are other AIV subtypes with low pathogenicity that have low reassortment rate, most notably H9N2 [<xref ref-type="bibr" rid="B44">44</xref>
].</p>
<p>We show that the reassortment rates among viral subtypes are correlated with their prevalence in different types of host. High reassortment rates could be associated with high prevalence of certain subtypes in wild birds and particularly among Anseriformes host species. In contrast, the lower reassortment rate of certain subtypes could be explained by their high prevalence in domestic birds (mainly in Galliformes). Wild bird and domestic bird populations represent two different ecosystems: natural and artificial respectively, in which influenza viruses face different ecologic constraints such as host population structure density, and opportunities for virus transmission [<xref ref-type="bibr" rid="B45">45</xref>
]. For instance, transmission by feces provides a way for wild birds to spread viruses to other wild birds as they migrate through an area [<xref ref-type="bibr" rid="B10">10</xref>
], while AIV in domestic birds, especially domestic birds of the order Galliformes may exist in a more isolated population (within a physically confined flock) than might occur in wild birds, thus with less possibility of reassortment.</p>
<p>In comparison to the AIV in domestic Galliformes, AIV at high prevalence in the domestic duck (the H6 subtype in our study) has a higher reassortment rate. Domestic ducks or geese farming occurs in high-density but free-range settings in Asia [<xref ref-type="bibr" rid="B46">46</xref>
]. This farming style creates an environment in which migratory birds and domestic ducks are in close contact, and provide more opportunities for AIV to reassort [<xref ref-type="bibr" rid="B47">47</xref>
]. In addition, different responses to virus infection may also influence the opportunity of AIV reassortment in certain host species. Avian influenza virus in chickens (Galliformes) often leads to severe respiratory symptoms and high mortality rates of up to 100%. In contrast, AIV infection generally causes no major clinical signs and immune response in ducks (Anseriformes) [<xref ref-type="bibr" rid="B48">48</xref>
]. Different influenza virus subtypes can also infect wild ducks concomitantly, creating the opportunity for genetic mixing [<xref ref-type="bibr" rid="B49">49</xref>
]. Those subtypes showing the high reassortment rates in our study, especially H1N1, H3N8 and H4N6 represent the most common subtype combinations that are found in wild ducks during periods of peak prevalence and have a maximum detectable subtype diversity in the late summer and fall [<xref ref-type="bibr" rid="B32">32</xref>
,<xref ref-type="bibr" rid="B50">50</xref>
].</p>
</sec>
<sec sec-type="conclusions"><title>Conclusions</title>
<p>In summary, we quantified the inter-subtype reassortment rates using discrete trait modeling upon time resolved phylogenies with respect to all six internal segments of Eurasian AIV. The overall reassortment rates for the external protein coding segments on internal protein coding segments were highly correlated with the internal segment genetic diversity (divergence time) and negatively correlated with d<sub>N</sub>
/d<sub>S</sub>
ratio and selective constraint. In addition, we found that different subtypes are characterized by reassortment rates that are strongly dependent on host type (wild <italic>v.</italic>
domestic). These results are consistent with the hypothesis that lack of immunological cross-protection at the subtype level in avian species allows frequent mixed infections in wild bird populations and provide novel insights into the evolutionary dynamics of avian influenza viruses.</p>
</sec>
<sec sec-type="methods"><title>Methods</title>
<sec><title>Data preparation and Preliminary phylogenetic analysis</title>
<p>Initial analyses included 3226 full length complete genome avian influenza sequences of any subtype and sampling year from the NCBI Influenza Virus Resource (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/genomes/FLU">http://www.ncbi.nlm.nih.gov/genomes/FLU</ext-link>
). Laboratory recombinants or highly cultured sequences were excluded.</p>
<p>The locations of samples were classified by countries and states/provinces separately. The host species were categorized by taxonomic bird Orders, and were further classified into wild birds and domestic birds: for instance, aquatic birds of the orders of Anseriformes (ducks, swans, etc.) and Charadriiformes (gulls, terns, etc.). Domestic birds can be further separated into domestic Anseriformes (ducks, geese etc.) and domestic Galliformes (chickens, turkeys etc.). The date of collection of each sequence was converted into fractional years (decimal date) for the estimation of divergence time, however for isolates where only year information was available the date was taken to be half way through the year, and if only month and year information was available the date was taken to be the 15<sup>th</sup>
of the month (shown in Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S5)”. For each internal segment, sequence alignments of coding regions were created using MUSCLE [<xref ref-type="bibr" rid="B51">51</xref>
] and then adjusted manually using BioEdit version 7.2.3 (<ext-link ext-link-type="uri" xlink:href="http://www.mbio.ncsu.edu/bioedit/page2.html">http://www.mbio.ncsu.edu/bioedit/page2.html</ext-link>
). The gene segment alignments were manually concatenated to generate a single alignment of complete genome sequences, and duplicate genomes (having 100% nucleotide identity) were removed.</p>
<p>A preliminary phylogenetic tree of each segment was generated using the Neighbor Joining method in MEGA 5. These phylogenies were analyzed using Path-O-Gen version 1.2 (<ext-link ext-link-type="uri" xlink:href="http://tree.bio.ed.ac.uk/software/pathogen/">http://tree.bio.ed.ac.uk/software/pathogen/</ext-link>
), which is a tool to investigate the temporal signal and ‘clock-likeness’ of molecular phylogenies. The plots of root-to-tip showed a good linear relationship between genetic divergence and sampling time of the AIV sequences, except two strains (A/pintail/ALB/179/1993 in PB2 and A/shorebird/DE/12/2004 and A/duck/Yokohama/aq10/2003 in HA, respectively), indicating that they might contain incorrect date information and were removed from the dataset. Thus, our final dataset comprised 2996 full genome sequences sampled over 55 years from 1956 to 2011. Phylogenetic trees for each segment were generated using RAxML [<xref ref-type="bibr" rid="B52">52</xref>
], employing maximum likelihood (ML) analyses with 500 bootstraps.</p>
<p>Phylogenetic analysis of the 2996 AIV sequences datasets confirms a clear separation of AIV sequences from east and west hemispheres, although some sequences are found to cross the boundaries, which confirmed that geographical separation of host species has shaped the avian influenza A virus gene pool into almost completely independently evolving Eurasian and American lineages (Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S1). We focused our study on the avian influenza sequences of the Eurasian lineage, which included 1140 strains and composed of 15 HA subtype (H1-16, without H14), 9 NA subtype and 60 combined HA-NA subtypes in total.</p>
</sec>
<sec><title>Time-scaled phylogenetic analysis</title>
<p>Time-scaled trees were generated for each internal protein-coding segment with a subsample of the 1140 sequences dataset with known isolation dates using BEAST (version 1.7.3). Different substitution models, clock models and tree models were evaluated for each segment by the calculation of the Bayes factor (BF), which is the ratio of the marginal likelihoods of the two models, and log BF > 3 was taken as indicating support for one model over another [<xref ref-type="bibr" rid="B53">53</xref>
] in Tracer v1.4 (<ext-link ext-link-type="uri" xlink:href="http://tree.bio.ed.ac.uk/software/tracer/">http://tree.bio.ed.ac.uk/software/tracer/</ext-link>
) [<xref ref-type="bibr" rid="B54">54</xref>
]. A general-time-reversal (GTR) model [<xref ref-type="bibr" rid="B55">55</xref>
] with gamma distributed rate heterogeneity of 4 rate categories (C4) was chosen for the M and NS segments for both of which have overlapping reading frames. For the PB2, PB1, PA and NP segments, the SRD06 nucleotide substitution model [<xref ref-type="bibr" rid="B56">56</xref>
] was chosen as the default nucleotide substitution model (preferred over GTR model with log BF from 3242 to 3289), which allows one Hasegawa, Kishino and Yano (HKY) model [<xref ref-type="bibr" rid="B57">57</xref>
] for codon positions 1 and 2 and a different HKY model for position 3, together with gamma-distributed rates across sites. A constant population size model (preferred to either an exponential model or a Bayesian skyride model with log BF from 7.23 to 7.52) and a relaxed uncorrelated log-normal molecular clock model (preferred to a strict clock model with log BF from 5.66 to 15.32) were chosen for all eight segments. The MCMC was run for 10<sup>8</sup>
steps and sampled every 10<sup>4</sup>
steps. Two independent runs were used in each segment to confirm the convergence and then combined. The Maximum Clade Credibility (MCC) trees were obtained and the mean ages of each node and the corresponding HPD were calculated by using Tree Annotator v1.7.3 in BEAST.</p>
</sec>
<sec><title>Discrete trait mapping and estimation of reassortment rates</title>
<p>To quantify the genetic association among subtypes (HA, NA and HA-NA) in AIV from Eurasia, We implemented the asymmetric discrete traits model with irreversible transitions (with the uniform prior with mean = 1 for the trait.rate) to reconstruct the ancestral subtype states of the internal nodes from the posterior time-scaled tree distribution [<xref ref-type="bibr" rid="B29">29</xref>
], using a relaxed uncorrelated log-normal molecular clock model, SRD60 model of nucleotide substitution and a constant-population coalescent process prior across the phylogeny of each segment. In addition, a BF test was constructed to identify the significant transition rates between discrete traits with the Bayesian stochastic search variable selection (BSSVS) extension of the discrete model. For the rates calculated from BSSVS, a BF test was adopted to identify significant non-zero transition rates (BF = 3) [<xref ref-type="bibr" rid="B28">28</xref>
]. For each of the six internal segments, we did 3 separate analyses on sequences labelled by HA subtype, NA subtype, and HA-NA combined subtype. The reassortment rate being analysed in this study was classified into three levels: a) The overall reassortment rate for an internal segment (e.g. PB2) reported for HA subtype and for NA and for combined HA-NA was the overall mutation rate of HA,NA or combined HA-NA traits in substitutions per site per year; b) The rate of reassortment for an internal segment (e.g. PB2) starting from each different subtype (e.g. H1, N1, or H1N1) is the mean rate that averaged over the individual transition rates from certain subtype to all the other subtypes; c) The asymmetric rate of reassortment for an internal segment (e.g.PB2) starting from one subtype to another (e.g. H1 to H2, H2 to H1) was estimated from each relative transition rate multiplied by the corresponding clock rate per Monte Carlo Markov Chain (MCMC) sample.</p>
<p>For comparison, the reassortment rates were also estimated by MultiState in BayesTraits (version 1.0) [<xref ref-type="bibr" rid="B58">58</xref>
], using a maximum likelihood method with 1000 bootstrap replications. Only the rates of HA and NA subtypes (not HA-NA joint subtype) could be estimated using Multistate, as only a maximum of 10 distinct states were allowed. The mean rates estimated by two different reconstruction methods (MCMC and ML) are similar (data not shown), although there was a significant shrinkage in uncertainty using the formal method (see Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S3). This result suggested that reassortment rates estimated with discrete trait models on the MCMC-based Bayesian framework were preferred.</p>
</sec>
<sec><title>Robustness of random sampling</title>
<p>In order to generate datasets suitable for obtaining time-resolved trees, stratified subsampling was conducted using custom R scripts. Subtypes with more than 7 sequences were subjected to further analysis. In addition, to keep the pattern of the original phylogeny, sequences in isolated branches with extremely long branches were ignored. To retain a similar prevalence of isolates for each subtype, random subsampling was conducted to keep one sequence per subtype per host species per location per year, and was repeated 3 times (Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S1, rows A1-A3). Among the three subgroups (A1, A2, and A3), the mean individual reassortment rate of each single transmission were similar to each other and with overlapping credible intervals (Cls) , which indicated that the estimated reassortment rates are robust by random sampling (Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S4). Thus, the reassortment rates shown in the results were the mean rate of the three subgroups.</p>
</sec>
<sec><title>Test of sample size effect</title>
<p>The dataset in this study was obtained from the NCBI database, where the number of sequences is unevenly distributed among different subtypes. HPAI H5N1 AIV in particular accounts for a very high proportion of sequences due to large-scale surveillance efforts. To investigate a possible sample size effect on the reassortment rate, we estimated reassortment rates using two other subsampled datasets which comprised different numbers of representative sequences from the original datasets of 344. To identify whether the sample size of each subtypes would affect the reassortment rate, we further select representative sequences of different subtypes in two other different ways: 1) by reducing the number of the most over-represented subtype to equal the number of the rare subtype in the original dataset (Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S1, row B); and 2) by reducing the number of sequences of all the subtypes to be the same as the rare subtype (Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S1, row C). Representative isolates of the down-sampled datasets were selected according to the tree phylogeny of the original maximum likelihood tree generated with all the Eurasian AIV, using a software tool which searches for sequence clusters [<xref ref-type="bibr" rid="B59">59</xref>
] (<ext-link ext-link-type="uri" xlink:href="http://hiv.bio.ed.ac.uk/software.html">http://hiv.bio.ed.ac.uk/software.html</ext-link>
), with 70% bootstrap support and a genetic distance lower than 0.045 specified. In this comparison, the reassortment rates among subtypes were similar in all the three datasets (Additional file <xref ref-type="supplementary-material" rid="S2">1</xref>
: Figure S7). Thus, we can conclude that the sample size does not affect the reassortment pattern of the Eurasian AIV in our study.</p>
</sec>
<sec><title>Estimation of evolutionary rates</title>
<p>The time scaled Bayesian phylogenies were generated with an uncorrelated log-normal relaxed molecular clock, and each branch was annotated with its mean substitution rate (ucld.mean). Additionally the nodes of the trees were annotated with the inferred HA, NA or HA-NA subtype states from the respective discrete trait models. The mean evolutionary rates per subtype were extracted from the trees, averaging over the rates on only those branches for which both the parent and child nodes were annotated with that subtype. The overall evolutionary rate of each segment was the ucld.mean rate of the MCMC tree phylogeny of each segment with 95% HPD intervals.</p>
</sec>
<sec><title>Estimation of the dN/dS ratio</title>
<p>Selection of certain subtype on the tree phylogeny was analyzed with implementation of robust counting on the trees generated with Bayesian discrete models in Beast, which reconstructed the synonymous and non-synonymous change counts using a 3-partition codon model [<xref ref-type="bibr" rid="B60">60</xref>
]. For segment 6 and 7, M1 and NS1 were analyzed instead of total length due to overlapping reading frames [<xref ref-type="bibr" rid="B61">61</xref>
]. For either conditional (C) and unconditional (U) stochastic mapping realizations, the numbers of synonymous (S) substitutions <italic>C</italic>
(S-C), <italic>C</italic>
(S-U) and nonsynonymous (N) substitutions <italic>C</italic>
(N-C), <italic>C</italic>
(N-U) for each state on each branch were recorded, and then converted into d<sub>N</sub>
/d<sub>S</sub>
values (indicated by certain state/subtype) by computing the ratio <italic>C</italic>
(N-C)<italic>/C</italic>
(S-C)<italic>/C</italic>
(N-U)<italic>/C</italic>
(S-U) [<xref ref-type="bibr" rid="B39">39</xref>
]. The ratio per subtype was estimated from the ratio per branch across all the MCMC trees. The overall selective pressure of each segment was the mean difference between the mean number of total N and total S of across all the MCMC trees of certain segment with 95%HPD intervals.</p>
<p>We also used Single Likelihood Ancestor Counting (SLAC) in HyPhy (<ext-link ext-link-type="uri" xlink:href="http://www.datamonkey.org/">http://www.datamonkey.org/</ext-link>
) to estimate dn/ds per subtype. Here we used the original datasets (1185 AIV in total) which have more sequences of each subtype than the subsampling dataset. The HKY85 nucleotide substitution model was used in SLAC-based analyses.</p>
</sec>
</sec>
<sec><title>Additional files</title>
<p>Models, priors and settings used in this analysis are provided in the XML used as input for BEAST 1.7.3, which are available from the Dryad Digital Repository: <ext-link ext-link-type="uri" xlink:href="http://doi.org/10.5061/dryad.246pr">http://doi.org/10.5061/dryad.246pr</ext-link>
.</p>
</sec>
<sec><title>Competing interests</title>
<p>The authors declare that no competing interests exist.</p>
</sec>
<sec><title>Authors’ contributions</title>
<p>LL performed the phylogenetic and discrete trait mapping analyses, was involved in the study design and drafted the manuscript. SJL provided the major scripts for R programs used in this study, helped with the study design and provided guidance on the analysis. AJLB conceived the study, provided guidance on its design and drafted the manuscript. All authors read and approved the final manuscript.</p>
</sec>
<sec sec-type="supplementary-material"><title>Supplementary Material</title>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title>
<p><bold>Figure S1.</bold>
Maximum likelihood tree of 2996 AIV of PB2 segment. <bold>Figure S2.</bold>
Bayesian maximum clade credibility phylogenies for the five internal segments (A: HA subtype, B: NA subtype, C: Combined HA-NA subtype) of Eurasian Avian influenza viruses. <bold>Figure S3.</bold>
Distribution of variance of reassortment rate using BSSVS and Bayes Traits. <bold>Figure S4.</bold>
Robustness of reassortment rate (HA, NA and HA-NA combined subtype) estimates using BSSVS among 3 random sampling subgroups. <bold>Figure S5.</bold>
Heat map of reassortment rates of remain 5 internal segments (PB1, PA, NP, M, NS) of different subtypes (HA, NA and both HA and NA). <bold>Figure S6.</bold>
Distribution of reassortment rates in the PB2 heat maps for H, N, and H-N subtypes in rows (blue) and columns (red) (Figure <xref ref-type="fig" rid="F2">2</xref>
and Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S2). <bold>Figure S7.</bold>
Heat map of reassortment rates of HA, NA and combined HA-NA subtypes of PB2 segment on three AIV datasets of different sample sizes (Group A, B and C are correspondent to Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S1). <bold>Figure S8.</bold>
Correlation between proportion of AIV in wild birds and reassortment rate of different subtypes. <bold>Figure S9.</bold>
Correlation between proportion of AIV in Anseriformes and reassortment rate of different subtypes. <bold>Figure S10.</bold>
Correlation between the proportion of AIV in wild anseriformes and the reassortment rate of different subtypes. <bold>Figure S11.</bold>
Correlation between proportion of AIV in domestic anseriformes and reassortment rate of different subtypes. <bold>Figure S12.</bold>
Correlation between proportion of AIV in domestic galliformes and reassortment rate per subtype. <bold>Figure S13.</bold>
Correlation between evolutionary rate and reassortment rate per subtype. <bold>Figure S14.</bold>
Correlation between inter-subtype reassortment rate and dN/dS (using SLAC).</p>
</caption>
<media xlink:href="1471-2148-14-16-S1.pdf"><caption><p>Click here for file</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title>
<p><bold>Table S1.</bold>
Number of sequences by subtype in each subgroup. <bold>Table S2.</bold>
Matrix of reassortment rate between pairs of subtypes (HA, NA and HA-NA combined) of 6 internal segments. <bold>Table S3.</bold>
Partial correlation between reassortment and domestic proportion or dN/dS ratio. <bold>Table S4.</bold>
Number of selected sites for 13 subtypes and the correlation with inter-subtype reassortment rate in each internal segment. <bold>Table S5.</bold>
Sequences for the discrete traits analysis. Information for subsampling data set A2 (n = 344, Additional file <xref ref-type="supplementary-material" rid="S1">2</xref>
: Table S1) of Eurasian AIV sequences used in this discrete traits analysis. Isolate name, subtype, host species, date of collection and the estimated decimal date are listed.</p>
</caption>
<media xlink:href="1471-2148-14-16-S2.pdf"><caption><p>Click here for file</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back><sec><title>Acknowledgments</title>
<p>This work was funded by a China Scholarships Council/University of Edinburgh Scholarship to L.L. SJL was supported by the Wellcome Trust (grant number 092807).</p>
</sec>
<ref-list><ref id="B1"><mixed-citation publication-type="journal"><name><surname>Webster</surname>
<given-names>RG</given-names>
</name>
<name><surname>Bean</surname>
<given-names>WJ</given-names>
</name>
<name><surname>Gorman</surname>
<given-names>OT</given-names>
</name>
<name><surname>Chambers</surname>
<given-names>TM</given-names>
</name>
<name><surname>Kawaoka</surname>
<given-names>Y</given-names>
</name>
<article-title>Evolution and ecology of influenza A viruses</article-title>
<source>Microbiol Rev</source>
<year>1992</year>
<volume>56</volume>
<issue>1</issue>
<fpage>152</fpage>
<lpage>179</lpage>
<pub-id pub-id-type="pmid">1579108</pub-id>
</mixed-citation>
</ref>
<ref id="B2"><mixed-citation publication-type="journal"><name><surname>McHardy</surname>
<given-names>AC</given-names>
</name>
<name><surname>Adams</surname>
<given-names>B</given-names>
</name>
<article-title>The role of genomics in tracking the evolution of influenza A virus</article-title>
<source>PLoS Pathog</source>
<year>2009</year>
<volume>5</volume>
<issue>10</issue>
<fpage>e1000566</fpage>
<pub-id pub-id-type="doi">10.1371/journal.ppat.1000566</pub-id>
<pub-id pub-id-type="pmid">19855818</pub-id>
</mixed-citation>
</ref>
<ref id="B3"><mixed-citation publication-type="journal"><name><surname>Webster</surname>
<given-names>RG</given-names>
</name>
<article-title>Influenza: an emerging disease</article-title>
<source>Emerg Infect Dis</source>
<year>1998</year>
<volume>4</volume>
<issue>3</issue>
<fpage>436</fpage>
<lpage>441</lpage>
<pub-id pub-id-type="doi">10.3201/eid0403.980325</pub-id>
<pub-id pub-id-type="pmid">9716966</pub-id>
</mixed-citation>
</ref>
<ref id="B4"><mixed-citation publication-type="journal"><name><surname>Parrish</surname>
<given-names>CR</given-names>
</name>
<name><surname>Kawaoka</surname>
<given-names>Y</given-names>
</name>
<article-title>The origins of new pandemic viruses: the acquisition of new host ranges by canine parvovirus and influenza A viruses</article-title>
<source>Annu Rev Microbiol</source>
<year>2005</year>
<volume>59</volume>
<fpage>553</fpage>
<lpage>586</lpage>
<pub-id pub-id-type="doi">10.1146/annurev.micro.59.030804.121059</pub-id>
<pub-id pub-id-type="pmid">16153179</pub-id>
</mixed-citation>
</ref>
<ref id="B5"><mixed-citation publication-type="journal"><name><surname>Kalthoff</surname>
<given-names>D</given-names>
</name>
<name><surname>Globig</surname>
<given-names>A</given-names>
</name>
<name><surname>Beer</surname>
<given-names>M</given-names>
</name>
<article-title>(Highly pathogenic) avian influenza as a zoonotic agent</article-title>
<source>Vet Microbiol</source>
<year>2010</year>
<volume>140</volume>
<issue>3–4</issue>
<fpage>237</fpage>
<lpage>245</lpage>
<pub-id pub-id-type="pmid">19782482</pub-id>
</mixed-citation>
</ref>
<ref id="B6"><mixed-citation publication-type="journal"><name><surname>Lindstrom</surname>
<given-names>SE</given-names>
</name>
<name><surname>Cox</surname>
<given-names>NJ</given-names>
</name>
<name><surname>Klimov</surname>
<given-names>A</given-names>
</name>
<article-title>Genetic analysis of human H2N2 and early H3N2 influenza viruses, 1957–1972: evidence for genetic divergence and multiple reassortment events</article-title>
<source>Virology</source>
<year>2004</year>
<volume>328</volume>
<issue>1</issue>
<fpage>101</fpage>
<lpage>119</lpage>
<pub-id pub-id-type="doi">10.1016/j.virol.2004.06.009</pub-id>
<pub-id pub-id-type="pmid">15380362</pub-id>
</mixed-citation>
</ref>
<ref id="B7"><mixed-citation publication-type="journal"><name><surname>Belshe</surname>
<given-names>RB</given-names>
</name>
<article-title>The origins of pandemic influenza–lessons from the 1918 virus</article-title>
<source>N Engl J Med</source>
<year>2005</year>
<volume>353</volume>
<issue>21</issue>
<fpage>2209</fpage>
<lpage>2211</lpage>
<pub-id pub-id-type="doi">10.1056/NEJMp058281</pub-id>
<pub-id pub-id-type="pmid">16306515</pub-id>
</mixed-citation>
</ref>
<ref id="B8"><mixed-citation publication-type="journal"><name><surname>Smith</surname>
<given-names>GJD</given-names>
</name>
<name><surname>Vijaykrishna</surname>
<given-names>D</given-names>
</name>
<name><surname>Bahl</surname>
<given-names>J</given-names>
</name>
<name><surname>Lycett</surname>
<given-names>SJ</given-names>
</name>
<name><surname>Worobey</surname>
<given-names>M</given-names>
</name>
<name><surname>Pybus</surname>
<given-names>OG</given-names>
</name>
<name><surname>Ma</surname>
<given-names>SK</given-names>
</name>
<name><surname>Cheung</surname>
<given-names>CL</given-names>
</name>
<name><surname>Raghwani</surname>
<given-names>J</given-names>
</name>
<name><surname>Bhatt</surname>
<given-names>S</given-names>
</name>
<etal></etal>
<article-title>Origins and evolutionary genomics of the 2009 swine-origin H1N1 influenza A epidemic</article-title>
<source>Nature</source>
<year>2009</year>
<volume>459</volume>
<issue>7250</issue>
<fpage>1122</fpage>
<lpage>1125</lpage>
<pub-id pub-id-type="doi">10.1038/nature08182</pub-id>
<pub-id pub-id-type="pmid">19516283</pub-id>
</mixed-citation>
</ref>
<ref id="B9"><mixed-citation publication-type="journal"><name><surname>Taubenberger</surname>
<given-names>JK</given-names>
</name>
<name><surname>Kash</surname>
<given-names>JC</given-names>
</name>
<article-title>Influenza virus evolution, host adaptation, and pandemic formation</article-title>
<source>Cell Host Microbe</source>
<year>2010</year>
<volume>7</volume>
<issue>6</issue>
<fpage>440</fpage>
<lpage>451</lpage>
<pub-id pub-id-type="doi">10.1016/j.chom.2010.05.009</pub-id>
<pub-id pub-id-type="pmid">20542248</pub-id>
</mixed-citation>
</ref>
<ref id="B10"><mixed-citation publication-type="journal"><name><surname>Webster</surname>
<given-names>RG</given-names>
</name>
<article-title>The importance of animal influenza for human disease</article-title>
<source>Vaccine</source>
<year>2002</year>
<volume>20</volume>
<issue>Suppl 2</issue>
<fpage>S16</fpage>
<lpage>20</lpage>
<pub-id pub-id-type="pmid">12110250</pub-id>
</mixed-citation>
</ref>
<ref id="B11"><mixed-citation publication-type="journal"><name><surname>Cheung</surname>
<given-names>CL</given-names>
</name>
<name><surname>Vijaykrishna</surname>
<given-names>D</given-names>
</name>
<name><surname>Smith</surname>
<given-names>GJ</given-names>
</name>
<name><surname>Fan</surname>
<given-names>XH</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>JX</given-names>
</name>
<name><surname>Bahl</surname>
<given-names>J</given-names>
</name>
<name><surname>Duan</surname>
<given-names>L</given-names>
</name>
<name><surname>Huang</surname>
<given-names>K</given-names>
</name>
<name><surname>Tai</surname>
<given-names>H</given-names>
</name>
<name><surname>Wang</surname>
<given-names>J</given-names>
</name>
<etal></etal>
<article-title>Establishment of influenza a virus (H6N1) in minor poultry species in southern China</article-title>
<source>J Virol</source>
<year>2007</year>
<volume>81</volume>
<issue>19</issue>
<fpage>10402</fpage>
<lpage>10412</lpage>
<pub-id pub-id-type="doi">10.1128/JVI.01157-07</pub-id>
<pub-id pub-id-type="pmid">17652385</pub-id>
</mixed-citation>
</ref>
<ref id="B12"><mixed-citation publication-type="journal"><name><surname>Huang</surname>
<given-names>K</given-names>
</name>
<name><surname>Bahl</surname>
<given-names>J</given-names>
</name>
<name><surname>Fan</surname>
<given-names>XH</given-names>
</name>
<name><surname>Vijaykrishna</surname>
<given-names>D</given-names>
</name>
<name><surname>Cheung</surname>
<given-names>CL</given-names>
</name>
<name><surname>Webby</surname>
<given-names>RJ</given-names>
</name>
<name><surname>Webster</surname>
<given-names>RG</given-names>
</name>
<name><surname>Chen</surname>
<given-names>H</given-names>
</name>
<name><surname>Smith</surname>
<given-names>GJ</given-names>
</name>
<name><surname>Peiris</surname>
<given-names>JS</given-names>
</name>
<etal></etal>
<article-title>Establishment of an H6N2 influenza virus lineage in domestic ducks in southern China</article-title>
<source>J Virol</source>
<year>2010</year>
<volume>84</volume>
<issue>14</issue>
<fpage>6978</fpage>
<lpage>6986</lpage>
<pub-id pub-id-type="doi">10.1128/JVI.00256-10</pub-id>
<pub-id pub-id-type="pmid">20463062</pub-id>
</mixed-citation>
</ref>
<ref id="B13"><mixed-citation publication-type="journal"><name><surname>Huang</surname>
<given-names>K</given-names>
</name>
<name><surname>Zhu</surname>
<given-names>H</given-names>
</name>
<name><surname>Fan</surname>
<given-names>X</given-names>
</name>
<name><surname>Wang</surname>
<given-names>J</given-names>
</name>
<name><surname>Cheung</surname>
<given-names>CL</given-names>
</name>
<name><surname>Duan</surname>
<given-names>L</given-names>
</name>
<name><surname>Hong</surname>
<given-names>W</given-names>
</name>
<name><surname>Liu</surname>
<given-names>Y</given-names>
</name>
<name><surname>Li</surname>
<given-names>L</given-names>
</name>
<name><surname>Smith</surname>
<given-names>DK</given-names>
</name>
<etal></etal>
<article-title>Establishment and lineage replacement of H6 influenza viruses in domestic ducks in southern China</article-title>
<source>J Virol</source>
<year>2012</year>
<volume>86</volume>
<issue>11</issue>
<fpage>6075</fpage>
<lpage>6083</lpage>
<pub-id pub-id-type="doi">10.1128/JVI.06389-11</pub-id>
<pub-id pub-id-type="pmid">22438558</pub-id>
</mixed-citation>
</ref>
<ref id="B14"><mixed-citation publication-type="journal"><name><surname>Wood</surname>
<given-names>GW</given-names>
</name>
<name><surname>McCauley</surname>
<given-names>JW</given-names>
</name>
<name><surname>Bashiruddin</surname>
<given-names>JB</given-names>
</name>
<name><surname>Alexander</surname>
<given-names>DJ</given-names>
</name>
<article-title>Deduced amino acid sequences at the haemagglutinin cleavage site of avian influenza a viruses of H5 and H7 subtypes</article-title>
<source>Arch Virol</source>
<year>1993</year>
<volume>130</volume>
<issue>1–2</issue>
<fpage>209</fpage>
<lpage>217</lpage>
<pub-id pub-id-type="pmid">8503786</pub-id>
</mixed-citation>
</ref>
<ref id="B15"><mixed-citation publication-type="journal"><name><surname>Klenk</surname>
<given-names>HD</given-names>
</name>
<name><surname>Garten</surname>
<given-names>W</given-names>
</name>
<article-title>Host cell proteases controlling virus pathogenicity</article-title>
<source>Trends Microbiol</source>
<year>1994</year>
<volume>2</volume>
<issue>2</issue>
<fpage>39</fpage>
<lpage>43</lpage>
<pub-id pub-id-type="doi">10.1016/0966-842X(94)90123-6</pub-id>
<pub-id pub-id-type="pmid">8162439</pub-id>
</mixed-citation>
</ref>
<ref id="B16"><mixed-citation publication-type="journal"><name><surname>Stienekegrober</surname>
<given-names>A</given-names>
</name>
<name><surname>Vey</surname>
<given-names>M</given-names>
</name>
<name><surname>Angliker</surname>
<given-names>H</given-names>
</name>
<name><surname>Shaw</surname>
<given-names>E</given-names>
</name>
<name><surname>Thomas</surname>
<given-names>G</given-names>
</name>
<name><surname>Roberts</surname>
<given-names>C</given-names>
</name>
<name><surname>Klenk</surname>
<given-names>HD</given-names>
</name>
<name><surname>Garten</surname>
<given-names>W</given-names>
</name>
<article-title>Influenza-Virus hemagglutinin with multibasic cleavage site is activated by Furin, a subtilisin-like endoprotease</article-title>
<source>Embo J</source>
<year>1992</year>
<volume>11</volume>
<issue>7</issue>
<fpage>2407</fpage>
<lpage>2414</lpage>
<pub-id pub-id-type="pmid">1628614</pub-id>
</mixed-citation>
</ref>
<ref id="B17"><mixed-citation publication-type="journal"><name><surname>Duan</surname>
<given-names>L</given-names>
</name>
<name><surname>Campitelli</surname>
<given-names>L</given-names>
</name>
<name><surname>Fan</surname>
<given-names>XH</given-names>
</name>
<name><surname>Leung</surname>
<given-names>YH</given-names>
</name>
<name><surname>Vijaykrishna</surname>
<given-names>D</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>JX</given-names>
</name>
<name><surname>Donatelli</surname>
<given-names>I</given-names>
</name>
<name><surname>Delogu</surname>
<given-names>M</given-names>
</name>
<name><surname>Li</surname>
<given-names>KS</given-names>
</name>
<name><surname>Foni</surname>
<given-names>E</given-names>
</name>
<etal></etal>
<article-title>Characterization of low-pathogenic H5 subtype influenza viruses from Eurasia: implications for the origin of highly pathogenic H5N1 viruses</article-title>
<source>J Virol</source>
<year>2007</year>
<volume>81</volume>
<issue>14</issue>
<fpage>7529</fpage>
<lpage>7539</lpage>
<pub-id pub-id-type="doi">10.1128/JVI.00327-07</pub-id>
<pub-id pub-id-type="pmid">17507485</pub-id>
</mixed-citation>
</ref>
<ref id="B18"><mixed-citation publication-type="journal"><name><surname>Campitelli</surname>
<given-names>L</given-names>
</name>
<name><surname>Mogavero</surname>
<given-names>E</given-names>
</name>
<name><surname>De Marco</surname>
<given-names>MA</given-names>
</name>
<name><surname>Delogu</surname>
<given-names>M</given-names>
</name>
<name><surname>Puzelli</surname>
<given-names>S</given-names>
</name>
<name><surname>Frezza</surname>
<given-names>F</given-names>
</name>
<name><surname>Facchini</surname>
<given-names>M</given-names>
</name>
<name><surname>Chiapponi</surname>
<given-names>C</given-names>
</name>
<name><surname>Foni</surname>
<given-names>E</given-names>
</name>
<name><surname>Cordioli</surname>
<given-names>P</given-names>
</name>
<etal></etal>
<article-title>Interspecies transmission of an H7N3 influenza virus from wild birds to intensively reared domestic poultry in Italy</article-title>
<source>Virology</source>
<year>2004</year>
<volume>323</volume>
<issue>1</issue>
<fpage>24</fpage>
<lpage>36</lpage>
<pub-id pub-id-type="doi">10.1016/j.virol.2004.02.015</pub-id>
<pub-id pub-id-type="pmid">15165816</pub-id>
</mixed-citation>
</ref>
<ref id="B19"><mixed-citation publication-type="journal"><name><surname>Lei</surname>
<given-names>F</given-names>
</name>
<name><surname>Shi</surname>
<given-names>W</given-names>
</name>
<article-title>Prospective of genomics in revealing transmission, reassortment and evolution of wildlife-borne avian influenza A (H5N1) Viruses</article-title>
<source>Curr Genomics</source>
<year>2011</year>
<volume>12</volume>
<issue>7</issue>
<fpage>466</fpage>
<lpage>474</lpage>
<pub-id pub-id-type="doi">10.2174/138920211797904052</pub-id>
<pub-id pub-id-type="pmid">22547954</pub-id>
</mixed-citation>
</ref>
<ref id="B20"><mixed-citation publication-type="journal"><name><surname>Claas</surname>
<given-names>EC</given-names>
</name>
<name><surname>Osterhaus</surname>
<given-names>AD</given-names>
</name>
<name><surname>van Beek</surname>
<given-names>R</given-names>
</name>
<name><surname>De Jong</surname>
<given-names>JC</given-names>
</name>
<name><surname>Rimmelzwaan</surname>
<given-names>GF</given-names>
</name>
<name><surname>Senne</surname>
<given-names>DA</given-names>
</name>
<name><surname>Krauss</surname>
<given-names>S</given-names>
</name>
<name><surname>Shortridge</surname>
<given-names>KF</given-names>
</name>
<name><surname>Webster</surname>
<given-names>RG</given-names>
</name>
<article-title>Human influenza A H5N1 virus related to a highly pathogenic avian influenza virus</article-title>
<source>Lancet</source>
<year>1998</year>
<volume>351</volume>
<issue>9101</issue>
<fpage>472</fpage>
<lpage>477</lpage>
<pub-id pub-id-type="doi">10.1016/S0140-6736(97)11212-0</pub-id>
<pub-id pub-id-type="pmid">9482438</pub-id>
</mixed-citation>
</ref>
<ref id="B21"><mixed-citation publication-type="journal"><name><surname>Chen</surname>
<given-names>H</given-names>
</name>
<name><surname>Smith</surname>
<given-names>GJ</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>SY</given-names>
</name>
<name><surname>Qin</surname>
<given-names>K</given-names>
</name>
<name><surname>Wang</surname>
<given-names>J</given-names>
</name>
<name><surname>Li</surname>
<given-names>KS</given-names>
</name>
<name><surname>Webster</surname>
<given-names>RG</given-names>
</name>
<name><surname>Peiris</surname>
<given-names>JS</given-names>
</name>
<name><surname>Guan</surname>
<given-names>Y</given-names>
</name>
<article-title>Avian flu: H5N1 virus outbreak in migratory waterfowl</article-title>
<source>Nature</source>
<year>2005</year>
<volume>436</volume>
<issue>7048</issue>
<fpage>191</fpage>
<lpage>192</lpage>
<pub-id pub-id-type="doi">10.1038/nature03974</pub-id>
<pub-id pub-id-type="pmid">16007072</pub-id>
</mixed-citation>
</ref>
<ref id="B22"><mixed-citation publication-type="journal"><name><surname>Feare</surname>
<given-names>CJ</given-names>
</name>
<article-title>Role of wild birds in the spread of highly pathogenic avian influenza virus H5N1 and implications for global surveillance</article-title>
<source>Avian Dis</source>
<year>2010</year>
<volume>54</volume>
<issue>1 Suppl</issue>
<fpage>201</fpage>
<lpage>212</lpage>
<pub-id pub-id-type="pmid">20521633</pub-id>
</mixed-citation>
</ref>
<ref id="B23"><mixed-citation publication-type="journal"><name><surname>Capua</surname>
<given-names>I</given-names>
</name>
<name><surname>Alexander</surname>
<given-names>DJ</given-names>
</name>
<article-title>Ecology, epidemiology and human health implications of avian influenza viruses: why do we need to share genetic data?</article-title>
<source>Zoonoses Public Health</source>
<year>2008</year>
<volume>55</volume>
<issue>1</issue>
<fpage>2</fpage>
<lpage>15</lpage>
<pub-id pub-id-type="doi">10.1111/j.1863-2378.2007.01081.x</pub-id>
<pub-id pub-id-type="pmid">18201321</pub-id>
</mixed-citation>
</ref>
<ref id="B24"><mixed-citation publication-type="journal"><name><surname>Horimoto</surname>
<given-names>T</given-names>
</name>
<name><surname>Kawaoka</surname>
<given-names>Y</given-names>
</name>
<article-title>Influenza: lessons from past pandemics, warnings from current incidents</article-title>
<source>Nat Rev Microbiol</source>
<year>2005</year>
<volume>3</volume>
<issue>8</issue>
<fpage>591</fpage>
<lpage>600</lpage>
<pub-id pub-id-type="doi">10.1038/nrmicro1208</pub-id>
<pub-id pub-id-type="pmid">16064053</pub-id>
</mixed-citation>
</ref>
<ref id="B25"><mixed-citation publication-type="journal"><name><surname>Drummond</surname>
<given-names>AJ</given-names>
</name>
<name><surname>Rambaut</surname>
<given-names>A</given-names>
</name>
<article-title>BEAST: Bayesian evolutionary analysis by sampling trees</article-title>
<source>BMC Evol Biol</source>
<year>2007</year>
<volume>7</volume>
<issue>1</issue>
<fpage>214</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2148-7-214</pub-id>
<pub-id pub-id-type="pmid">17996036</pub-id>
</mixed-citation>
</ref>
<ref id="B26"><mixed-citation publication-type="journal"><name><surname>Lemey</surname>
<given-names>P</given-names>
</name>
<name><surname>Rambaut</surname>
<given-names>A</given-names>
</name>
<name><surname>Drummond</surname>
<given-names>AJ</given-names>
</name>
<name><surname>Suchard</surname>
<given-names>MA</given-names>
</name>
<article-title>Bayesian phylogeography finds its roots</article-title>
<source>PLoS Comput Biol</source>
<year>2009</year>
<volume>5</volume>
<issue>9</issue>
<fpage>e1000520</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1000520</pub-id>
<pub-id pub-id-type="pmid">19779555</pub-id>
</mixed-citation>
</ref>
<ref id="B27"><mixed-citation publication-type="journal"><name><surname>Lemey</surname>
<given-names>P</given-names>
</name>
<name><surname>Rambaut</surname>
<given-names>A</given-names>
</name>
<name><surname>Welch</surname>
<given-names>JJ</given-names>
</name>
<name><surname>Suchard</surname>
<given-names>MA</given-names>
</name>
<article-title>Phylogeography takes a relaxed random walk in continuous space and time</article-title>
<source>Mol Biol Evol</source>
<year>2010</year>
<volume>27</volume>
<issue>8</issue>
<fpage>1877</fpage>
<lpage>1885</lpage>
<pub-id pub-id-type="doi">10.1093/molbev/msq067</pub-id>
<pub-id pub-id-type="pmid">20203288</pub-id>
</mixed-citation>
</ref>
<ref id="B28"><mixed-citation publication-type="journal"><name><surname>Faria</surname>
<given-names>NR</given-names>
</name>
<name><surname>Suchard</surname>
<given-names>MA</given-names>
</name>
<name><surname>Rambaut</surname>
<given-names>A</given-names>
</name>
<name><surname>Lemey</surname>
<given-names>P</given-names>
</name>
<article-title>Toward a quantitative understanding of viral phylogeography</article-title>
<source>Curr Opin Virol</source>
<year>2011</year>
<volume>1</volume>
<issue>5</issue>
<fpage>423</fpage>
<lpage>429</lpage>
<pub-id pub-id-type="doi">10.1016/j.coviro.2011.10.003</pub-id>
<pub-id pub-id-type="pmid">22440846</pub-id>
</mixed-citation>
</ref>
<ref id="B29"><mixed-citation publication-type="journal"><name><surname>Lycett</surname>
<given-names>SJ</given-names>
</name>
<name><surname>Baillie</surname>
<given-names>G</given-names>
</name>
<name><surname>Coulter</surname>
<given-names>E</given-names>
</name>
<name><surname>Bhatt</surname>
<given-names>S</given-names>
</name>
<name><surname>Kellam</surname>
<given-names>P</given-names>
</name>
<name><surname>McCauley</surname>
<given-names>JW</given-names>
</name>
<name><surname>Wood</surname>
<given-names>JL</given-names>
</name>
<name><surname>Brown</surname>
<given-names>IH</given-names>
</name>
<name><surname>Pybus</surname>
<given-names>OG</given-names>
</name>
<name><surname>Leigh Brown</surname>
<given-names>AJ</given-names>
</name>
<article-title>Estimating reassortment rates in co-circulating Eurasian swine influenza viruses</article-title>
<source>J Gen Virol</source>
<year>2012</year>
<volume>93</volume>
<issue>Pt 11</issue>
<fpage>2326</fpage>
<lpage>2336</lpage>
<pub-id pub-id-type="pmid">22971819</pub-id>
</mixed-citation>
</ref>
<ref id="B30"><mixed-citation publication-type="journal"><name><surname>Li</surname>
<given-names>KS</given-names>
</name>
<name><surname>Guan</surname>
<given-names>Y</given-names>
</name>
<name><surname>Wang</surname>
<given-names>J</given-names>
</name>
<name><surname>Smith</surname>
<given-names>GJ</given-names>
</name>
<name><surname>Xu</surname>
<given-names>KM</given-names>
</name>
<name><surname>Duan</surname>
<given-names>L</given-names>
</name>
<name><surname>Rahardjo</surname>
<given-names>AP</given-names>
</name>
<name><surname>Puthavathana</surname>
<given-names>P</given-names>
</name>
<name><surname>Buranathai</surname>
<given-names>C</given-names>
</name>
<name><surname>Nguyen</surname>
<given-names>TD</given-names>
</name>
<etal></etal>
<article-title>Genesis of a highly pathogenic and potentially pandemic H5N1 influenza virus in eastern Asia</article-title>
<source>Nature</source>
<year>2004</year>
<volume>430</volume>
<issue>6996</issue>
<fpage>209</fpage>
<lpage>213</lpage>
<pub-id pub-id-type="doi">10.1038/nature02746</pub-id>
<pub-id pub-id-type="pmid">15241415</pub-id>
</mixed-citation>
</ref>
<ref id="B31"><mixed-citation publication-type="journal"><name><surname>Dong</surname>
<given-names>G</given-names>
</name>
<name><surname>Luo</surname>
<given-names>J</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>H</given-names>
</name>
<name><surname>Wang</surname>
<given-names>C</given-names>
</name>
<name><surname>Duan</surname>
<given-names>M</given-names>
</name>
<name><surname>Deliberto</surname>
<given-names>TJ</given-names>
</name>
<name><surname>Nolte</surname>
<given-names>DL</given-names>
</name>
<name><surname>Ji</surname>
<given-names>G</given-names>
</name>
<name><surname>He</surname>
<given-names>H</given-names>
</name>
<article-title>Phylogenetic diversity and genotypical complexity of H9N2 influenza A viruses revealed by genomic sequence analysis</article-title>
<source>PLoS One</source>
<year>2011</year>
<volume>6</volume>
<issue>2</issue>
<fpage>e17212</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0017212</pub-id>
<pub-id pub-id-type="pmid">21386964</pub-id>
</mixed-citation>
</ref>
<ref id="B32"><mixed-citation publication-type="journal"><name><surname>Dugan</surname>
<given-names>VG</given-names>
</name>
<name><surname>Chen</surname>
<given-names>R</given-names>
</name>
<name><surname>Spiro</surname>
<given-names>DJ</given-names>
</name>
<name><surname>Sengamalay</surname>
<given-names>N</given-names>
</name>
<name><surname>Zaborsky</surname>
<given-names>J</given-names>
</name>
<name><surname>Ghedin</surname>
<given-names>E</given-names>
</name>
<name><surname>Nolting</surname>
<given-names>J</given-names>
</name>
<name><surname>Swayne</surname>
<given-names>DE</given-names>
</name>
<name><surname>Runstadler</surname>
<given-names>JA</given-names>
</name>
<name><surname>Happ</surname>
<given-names>GM</given-names>
</name>
<etal></etal>
<article-title>The evolutionary genetics and emergence of avian influenza viruses in wild birds</article-title>
<source>PLoS Pathog</source>
<year>2008</year>
<volume>4</volume>
<issue>5</issue>
<fpage>e1000076</fpage>
<pub-id pub-id-type="doi">10.1371/journal.ppat.1000076</pub-id>
<pub-id pub-id-type="pmid">18516303</pub-id>
</mixed-citation>
</ref>
<ref id="B33"><mixed-citation publication-type="journal"><name><surname>Baba</surname>
<given-names>K</given-names>
</name>
<name><surname>Shibata</surname>
<given-names>R</given-names>
</name>
<name><surname>Sibuya</surname>
<given-names>M</given-names>
</name>
<article-title>Partial correlation and conditional correlation as measures of conditional independence</article-title>
<source>Aust NZ J Stat</source>
<year>2004</year>
<volume>46</volume>
<issue>4</issue>
<fpage>657</fpage>
<lpage>664</lpage>
<pub-id pub-id-type="doi">10.1111/j.1467-842X.2004.00360.x</pub-id>
</mixed-citation>
</ref>
<ref id="B34"><mixed-citation publication-type="journal"><name><surname>Rambaut</surname>
<given-names>A</given-names>
</name>
<name><surname>Pybus</surname>
<given-names>OG</given-names>
</name>
<name><surname>Nelson</surname>
<given-names>MI</given-names>
</name>
<name><surname>Viboud</surname>
<given-names>C</given-names>
</name>
<name><surname>Taubenberger</surname>
<given-names>JK</given-names>
</name>
<name><surname>Holmes</surname>
<given-names>EC</given-names>
</name>
<article-title>The genomic and epidemiological dynamics of human influenza A virus</article-title>
<source>Nature</source>
<year>2008</year>
<volume>453</volume>
<issue>7195</issue>
<fpage>615</fpage>
<lpage>619</lpage>
<pub-id pub-id-type="doi">10.1038/nature06945</pub-id>
<pub-id pub-id-type="pmid">18418375</pub-id>
</mixed-citation>
</ref>
<ref id="B35"><mixed-citation publication-type="journal"><name><surname>Ludwig</surname>
<given-names>S</given-names>
</name>
<name><surname>Schultz</surname>
<given-names>U</given-names>
</name>
<name><surname>Mandler</surname>
<given-names>J</given-names>
</name>
<name><surname>Fitch</surname>
<given-names>WM</given-names>
</name>
<name><surname>Scholtissek</surname>
<given-names>C</given-names>
</name>
<article-title>Phylogenetic relationship of the nonstructural (NS) genes of influenza A viruses</article-title>
<source>Virology</source>
<year>1991</year>
<volume>183</volume>
<issue>2</issue>
<fpage>566</fpage>
<lpage>577</lpage>
<pub-id pub-id-type="doi">10.1016/0042-6822(91)90985-K</pub-id>
<pub-id pub-id-type="pmid">1830182</pub-id>
</mixed-citation>
</ref>
<ref id="B36"><mixed-citation publication-type="journal"><name><surname>Treanor</surname>
<given-names>JJ</given-names>
</name>
<name><surname>Snyder</surname>
<given-names>MH</given-names>
</name>
<name><surname>London</surname>
<given-names>WT</given-names>
</name>
<name><surname>Murphy</surname>
<given-names>BR</given-names>
</name>
<article-title>The B allele of the NS gene of avian influenza viruses, but not the A allele, attenuates a human influenza A virus for squirrel monkeys</article-title>
<source>Virology</source>
<year>1989</year>
<volume>171</volume>
<issue>1</issue>
<fpage>1</fpage>
<lpage>9</lpage>
<pub-id pub-id-type="doi">10.1016/0042-6822(89)90504-7</pub-id>
<pub-id pub-id-type="pmid">2525836</pub-id>
</mixed-citation>
</ref>
<ref id="B37"><mixed-citation publication-type="journal"><name><surname>Wang</surname>
<given-names>S</given-names>
</name>
<name><surname>Shi</surname>
<given-names>WM</given-names>
</name>
<name><surname>Mweene</surname>
<given-names>A</given-names>
</name>
<name><surname>Wei</surname>
<given-names>HL</given-names>
</name>
<name><surname>Bai</surname>
<given-names>GR</given-names>
</name>
<name><surname>Liu</surname>
<given-names>JH</given-names>
</name>
<article-title>Genetic analysis of the nonstructural (NS) genes of H9N2 chicken influenza viruses isolated in China during 1998–2002</article-title>
<source>Virus Genes</source>
<year>2005</year>
<volume>31</volume>
<issue>3</issue>
<fpage>329</fpage>
<lpage>335</lpage>
<pub-id pub-id-type="doi">10.1007/s11262-005-3251-2</pub-id>
<pub-id pub-id-type="pmid">16175338</pub-id>
</mixed-citation>
</ref>
<ref id="B38"><mixed-citation publication-type="journal"><name><surname>Suarez</surname>
<given-names>DL</given-names>
</name>
<name><surname>Perdue</surname>
<given-names>ML</given-names>
</name>
<article-title>Multiple alignment comparison of the non-structural genes of influenza A viruses</article-title>
<source>Virus Res</source>
<year>1998</year>
<volume>54</volume>
<issue>1</issue>
<fpage>59</fpage>
<lpage>69</lpage>
<pub-id pub-id-type="doi">10.1016/S0168-1702(98)00011-2</pub-id>
<pub-id pub-id-type="pmid">9660072</pub-id>
</mixed-citation>
</ref>
<ref id="B39"><mixed-citation publication-type="journal"><name><surname>Lemey</surname>
<given-names>P</given-names>
</name>
<name><surname>Minin</surname>
<given-names>VN</given-names>
</name>
<name><surname>Bielejec</surname>
<given-names>F</given-names>
</name>
<name><surname>Kosakovsky Pond</surname>
<given-names>SL</given-names>
</name>
<name><surname>Suchard</surname>
<given-names>MA</given-names>
</name>
<article-title>A counting renaissance: combining stochastic mapping and empirical Bayes to quickly detect amino acid sites under positive selection</article-title>
<source>Bioinformatics</source>
<year>2012</year>
<volume>28</volume>
<issue>24</issue>
<fpage>3248</fpage>
<lpage>3256</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/bts580</pub-id>
<pub-id pub-id-type="pmid">23064000</pub-id>
</mixed-citation>
</ref>
<ref id="B40"><mixed-citation publication-type="journal"><name><surname>Kosakovsky Pond</surname>
<given-names>SL</given-names>
</name>
<name><surname>Frost</surname>
<given-names>SD</given-names>
</name>
<article-title>Not so different after all: a comparison of methods for detecting amino acid sites under selection</article-title>
<source>Mol Biol Evol</source>
<year>2005</year>
<volume>22</volume>
<issue>5</issue>
<fpage>1208</fpage>
<lpage>1222</lpage>
<pub-id pub-id-type="doi">10.1093/molbev/msi105</pub-id>
<pub-id pub-id-type="pmid">15703242</pub-id>
</mixed-citation>
</ref>
<ref id="B41"><mixed-citation publication-type="journal"><name><surname>Bahl</surname>
<given-names>J</given-names>
</name>
<name><surname>Vijaykrishna</surname>
<given-names>D</given-names>
</name>
<name><surname>Holmes</surname>
<given-names>EC</given-names>
</name>
<name><surname>Smith</surname>
<given-names>GJ</given-names>
</name>
<name><surname>Guan</surname>
<given-names>Y</given-names>
</name>
<article-title>Gene flow and competitive exclusion of avian influenza A virus in natural reservoir hosts</article-title>
<source>Virology</source>
<year>2009</year>
<volume>390</volume>
<issue>2</issue>
<fpage>289</fpage>
<lpage>297</lpage>
<pub-id pub-id-type="doi">10.1016/j.virol.2009.05.002</pub-id>
<pub-id pub-id-type="pmid">19501380</pub-id>
</mixed-citation>
</ref>
<ref id="B42"><mixed-citation publication-type="journal"><name><surname>Chen</surname>
<given-names>R</given-names>
</name>
<name><surname>Holmes</surname>
<given-names>EC</given-names>
</name>
<article-title>Avian influenza virus exhibits rapid evolutionary dynamics</article-title>
<source>Mol Biol Evol</source>
<year>2006</year>
<volume>23</volume>
<issue>12</issue>
<fpage>2336</fpage>
<lpage>2341</lpage>
<pub-id pub-id-type="doi">10.1093/molbev/msl102</pub-id>
<pub-id pub-id-type="pmid">16945980</pub-id>
</mixed-citation>
</ref>
<ref id="B43"><mixed-citation publication-type="journal"><name><surname>Cauthen</surname>
<given-names>AN</given-names>
</name>
<name><surname>Swayne</surname>
<given-names>DE</given-names>
</name>
<name><surname>Schultz-Cherry</surname>
<given-names>S</given-names>
</name>
<name><surname>Perdue</surname>
<given-names>ML</given-names>
</name>
<name><surname>Suarez</surname>
<given-names>DL</given-names>
</name>
<article-title>Continued circulation in China of highly pathogenic avian influenza viruses encoding the hemagglutinin gene associated with the 1997 H5N1 outbreak in poultry and humans</article-title>
<source>J Virol</source>
<year>2000</year>
<volume>74</volume>
<issue>14</issue>
<fpage>6592</fpage>
<lpage>6599</lpage>
<pub-id pub-id-type="doi">10.1128/JVI.74.14.6592-6599.2000</pub-id>
<pub-id pub-id-type="pmid">10864673</pub-id>
</mixed-citation>
</ref>
<ref id="B44"><mixed-citation publication-type="journal"><name><surname>Li</surname>
<given-names>C</given-names>
</name>
<name><surname>Yu</surname>
<given-names>K</given-names>
</name>
<name><surname>Tian</surname>
<given-names>G</given-names>
</name>
<name><surname>Yu</surname>
<given-names>D</given-names>
</name>
<name><surname>Liu</surname>
<given-names>L</given-names>
</name>
<name><surname>Jing</surname>
<given-names>B</given-names>
</name>
<name><surname>Ping</surname>
<given-names>J</given-names>
</name>
<name><surname>Chen</surname>
<given-names>H</given-names>
</name>
<article-title>Evolution of H9N2 influenza viruses from domestic poultry in Mainland China</article-title>
<source>Virology</source>
<year>2005</year>
<volume>340</volume>
<issue>1</issue>
<fpage>70</fpage>
<lpage>83</lpage>
<pub-id pub-id-type="doi">10.1016/j.virol.2005.06.025</pub-id>
<pub-id pub-id-type="pmid">16026813</pub-id>
</mixed-citation>
</ref>
<ref id="B45"><mixed-citation publication-type="journal"><name><surname>Lebarbenchon</surname>
<given-names>C</given-names>
</name>
<name><surname>Feare</surname>
<given-names>CJ</given-names>
</name>
<name><surname>Renaud</surname>
<given-names>F</given-names>
</name>
<name><surname>Thomas</surname>
<given-names>F</given-names>
</name>
<name><surname>Gauthier-Clerc</surname>
<given-names>M</given-names>
</name>
<article-title>Persistence of highly pathogenic avian influenza viruses in natural ecosystems</article-title>
<source>Emerg Infect Dis</source>
<year>2010</year>
<volume>16</volume>
<issue>7</issue>
<fpage>1057</fpage>
<lpage>1062</lpage>
<pub-id pub-id-type="doi">10.3201/eid1607.090389</pub-id>
<pub-id pub-id-type="pmid">20587174</pub-id>
</mixed-citation>
</ref>
<ref id="B46"><mixed-citation publication-type="journal"><name><surname>Gu</surname>
<given-names>M</given-names>
</name>
<name><surname>Liu</surname>
<given-names>W</given-names>
</name>
<name><surname>Cao</surname>
<given-names>Y</given-names>
</name>
<name><surname>Peng</surname>
<given-names>D</given-names>
</name>
<name><surname>Wang</surname>
<given-names>X</given-names>
</name>
<name><surname>Wan</surname>
<given-names>H</given-names>
</name>
<name><surname>Zhao</surname>
<given-names>G</given-names>
</name>
<name><surname>Xu</surname>
<given-names>Q</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>W</given-names>
</name>
<name><surname>Song</surname>
<given-names>Q</given-names>
</name>
<etal></etal>
<article-title>Novel reassortant highly pathogenic avian influenza (H5N5) viruses in domestic ducks</article-title>
<source>China. Emerg Infect Dis</source>
<year>2011</year>
<volume>17</volume>
<issue>6</issue>
<fpage>1060</fpage>
<lpage>1063</lpage>
</mixed-citation>
</ref>
<ref id="B47"><mixed-citation publication-type="journal"><name><surname>Deng</surname>
<given-names>G</given-names>
</name>
<name><surname>Tan</surname>
<given-names>D</given-names>
</name>
<name><surname>Shi</surname>
<given-names>J</given-names>
</name>
<name><surname>Cui</surname>
<given-names>P</given-names>
</name>
<name><surname>Jiang</surname>
<given-names>Y</given-names>
</name>
<name><surname>Liu</surname>
<given-names>L</given-names>
</name>
<name><surname>Tian</surname>
<given-names>G</given-names>
</name>
<name><surname>Kawaoka</surname>
<given-names>Y</given-names>
</name>
<name><surname>Li</surname>
<given-names>C</given-names>
</name>
<name><surname>Chen</surname>
<given-names>H</given-names>
</name>
<article-title>Complex reassortment of multiple subtypes of avian influenza viruses in domestic ducks at the Dongting Lake region of China</article-title>
<source>J Virol</source>
<year>2013</year>
<volume>87</volume>
<issue>17</issue>
<fpage>9452</fpage>
<lpage>9462</lpage>
<pub-id pub-id-type="doi">10.1128/JVI.00776-13</pub-id>
<pub-id pub-id-type="pmid">23804642</pub-id>
</mixed-citation>
</ref>
<ref id="B48"><mixed-citation publication-type="journal"><name><surname>Kida</surname>
<given-names>H</given-names>
</name>
<name><surname>Yanagawa</surname>
<given-names>R</given-names>
</name>
<name><surname>Matsuoka</surname>
<given-names>Y</given-names>
</name>
<article-title>Duck influenza lacking evidence of disease signs and immune response</article-title>
<source>Infect Immun</source>
<year>1980</year>
<volume>30</volume>
<issue>2</issue>
<fpage>547</fpage>
<lpage>553</lpage>
<pub-id pub-id-type="pmid">7439994</pub-id>
</mixed-citation>
</ref>
<ref id="B49"><mixed-citation publication-type="journal"><name><surname>Sharp</surname>
<given-names>GB</given-names>
</name>
<name><surname>Kawaoka</surname>
<given-names>Y</given-names>
</name>
<name><surname>Jones</surname>
<given-names>DJ</given-names>
</name>
<name><surname>Bean</surname>
<given-names>WJ</given-names>
</name>
<name><surname>Pryor</surname>
<given-names>SP</given-names>
</name>
<name><surname>Hinshaw</surname>
<given-names>V</given-names>
</name>
<name><surname>Webster</surname>
<given-names>RG</given-names>
</name>
<article-title>Coinfection of wild ducks by influenza A viruses: distribution patterns and biological significance</article-title>
<source>J Virol</source>
<year>1997</year>
<volume>71</volume>
<issue>8</issue>
<fpage>6128</fpage>
<lpage>6135</lpage>
<pub-id pub-id-type="pmid">9223507</pub-id>
</mixed-citation>
</ref>
<ref id="B50"><mixed-citation publication-type="journal"><name><surname>Olsen</surname>
<given-names>B</given-names>
</name>
<name><surname>Munster</surname>
<given-names>VJ</given-names>
</name>
<name><surname>Wallensten</surname>
<given-names>A</given-names>
</name>
<name><surname>Waldenstrom</surname>
<given-names>J</given-names>
</name>
<name><surname>Osterhaus</surname>
<given-names>AD</given-names>
</name>
<name><surname>Fouchier</surname>
<given-names>RA</given-names>
</name>
<article-title>Global patterns of influenza a virus in wild birds</article-title>
<source>Science</source>
<year>2006</year>
<volume>312</volume>
<issue>5772</issue>
<fpage>384</fpage>
<lpage>388</lpage>
<pub-id pub-id-type="doi">10.1126/science.1122438</pub-id>
<pub-id pub-id-type="pmid">16627734</pub-id>
</mixed-citation>
</ref>
<ref id="B51"><mixed-citation publication-type="journal"><name><surname>Edgar</surname>
<given-names>RC</given-names>
</name>
<article-title>MUSCLE: multiple sequence alignment with high accuracy and high throughput</article-title>
<source>Nucleic Acids Res</source>
<year>2004</year>
<volume>32</volume>
<issue>5</issue>
<fpage>1792</fpage>
<lpage>1797</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkh340</pub-id>
<pub-id pub-id-type="pmid">15034147</pub-id>
</mixed-citation>
</ref>
<ref id="B52"><mixed-citation publication-type="journal"><name><surname>Stamatakis</surname>
<given-names>A</given-names>
</name>
<name><surname>Ludwig</surname>
<given-names>T</given-names>
</name>
<name><surname>Meier</surname>
<given-names>H</given-names>
</name>
<article-title>RAxML-III: a fast program for maximum likelihood-based inference of large phylogenetic trees</article-title>
<source>Bioinformatics</source>
<year>2005</year>
<volume>21</volume>
<issue>4</issue>
<fpage>456</fpage>
<lpage>463</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/bti191</pub-id>
<pub-id pub-id-type="pmid">15608047</pub-id>
</mixed-citation>
</ref>
<ref id="B53"><mixed-citation publication-type="journal"><name><surname>Raftery</surname>
<given-names>MANAE</given-names>
</name>
<article-title>Approximating Bayesian inference with the weighted likelihood bootstrap</article-title>
<source>J R Stat Soc</source>
<year>1994</year>
<volume>56</volume>
<fpage>3</fpage>
<lpage>48</lpage>
</mixed-citation>
</ref>
<ref id="B54"><mixed-citation publication-type="journal"><name><surname>Drummond</surname>
<given-names>AJ</given-names>
</name>
<name><surname>Ho</surname>
<given-names>SY</given-names>
</name>
<name><surname>Phillips</surname>
<given-names>MJ</given-names>
</name>
<name><surname>Rambaut</surname>
<given-names>A</given-names>
</name>
<article-title>Relaxed phylogenetics and dating with confidence</article-title>
<source>PLoS Biol</source>
<year>2006</year>
<volume>4</volume>
<issue>5</issue>
<fpage>e88</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pbio.0040088</pub-id>
<pub-id pub-id-type="pmid">16683862</pub-id>
</mixed-citation>
</ref>
<ref id="B55"><mixed-citation publication-type="journal"><name><surname>Tavaré</surname>
<given-names>S</given-names>
</name>
<article-title>Some probabilistic and statistical problems on the analysis of DNA sequences</article-title>
<source>Lect Math Life Sci</source>
<year>1986</year>
<volume>17</volume>
<fpage>57</fpage>
<lpage>86</lpage>
</mixed-citation>
</ref>
<ref id="B56"><mixed-citation publication-type="journal"><name><surname>Shapiro</surname>
<given-names>B</given-names>
</name>
<name><surname>Rambaut</surname>
<given-names>A</given-names>
</name>
<name><surname>Drummond</surname>
<given-names>AJ</given-names>
</name>
<article-title>Choosing appropriate substitution models for the phylogenetic analysis of protein-coding sequences</article-title>
<source>Mol Biol Evol</source>
<year>2006</year>
<volume>23</volume>
<issue>1</issue>
<fpage>7</fpage>
<lpage>9</lpage>
<pub-id pub-id-type="pmid">16177232</pub-id>
</mixed-citation>
</ref>
<ref id="B57"><mixed-citation publication-type="journal"><name><surname>Hasegawa</surname>
<given-names>M</given-names>
</name>
<name><surname>Kishino</surname>
<given-names>H</given-names>
</name>
<name><surname>Yano</surname>
<given-names>T</given-names>
</name>
<article-title>Dating of the human-ape splitting by a molecular clock of mitochondrial DNA</article-title>
<source>J Mol Evol</source>
<year>1985</year>
<volume>22</volume>
<issue>2</issue>
<fpage>160</fpage>
<lpage>174</lpage>
<pub-id pub-id-type="doi">10.1007/BF02101694</pub-id>
<pub-id pub-id-type="pmid">3934395</pub-id>
</mixed-citation>
</ref>
<ref id="B58"><mixed-citation publication-type="journal"><name><surname>Pagel</surname>
<given-names>M</given-names>
</name>
<name><surname>Meade</surname>
<given-names>A</given-names>
</name>
<name><surname>Barker</surname>
<given-names>D</given-names>
</name>
<article-title>Bayesian estimation of ancestral character states on phylogenies</article-title>
<source>Syst Biol</source>
<year>2004</year>
<volume>53</volume>
<issue>5</issue>
<fpage>673</fpage>
<lpage>684</lpage>
<pub-id pub-id-type="doi">10.1080/10635150490522232</pub-id>
<pub-id pub-id-type="pmid">15545248</pub-id>
</mixed-citation>
</ref>
<ref id="B59"><mixed-citation publication-type="journal"><name><surname>Ragonnet-Cronin</surname>
<given-names>M</given-names>
</name>
<name><surname>Hodcroft</surname>
<given-names>E</given-names>
</name>
<name><surname>Hue</surname>
<given-names>S</given-names>
</name>
<name><surname>Fearnhill</surname>
<given-names>E</given-names>
</name>
<name><surname>Delpech</surname>
<given-names>V</given-names>
</name>
<name><surname>Leigh Brown</surname>
<given-names>AJ</given-names>
</name>
<name><surname>Lycett</surname>
<given-names>S</given-names>
</name>
<article-title>Automated analysis of phylogenetic clusters</article-title>
<source>BMC Bioinformatics</source>
<year>2013</year>
<volume>14</volume>
<issue>1</issue>
<fpage>317</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2105-14-317</pub-id>
<pub-id pub-id-type="pmid">24191891</pub-id>
</mixed-citation>
</ref>
<ref id="B60"><mixed-citation publication-type="journal"><name><surname>O'Brien</surname>
<given-names>JD</given-names>
</name>
<name><surname>Minin</surname>
<given-names>VN</given-names>
</name>
<name><surname>Suchard</surname>
<given-names>MA</given-names>
</name>
<article-title>Learning to count: robust estimates for labeled distances between molecular sequences</article-title>
<source>Mol Biol Evol</source>
<year>2009</year>
<volume>26</volume>
<issue>4</issue>
<fpage>801</fpage>
<lpage>814</lpage>
<pub-id pub-id-type="doi">10.1093/molbev/msp003</pub-id>
<pub-id pub-id-type="pmid">19131426</pub-id>
</mixed-citation>
</ref>
<ref id="B61"><mixed-citation publication-type="journal"><name><surname>Ghedin</surname>
<given-names>E</given-names>
</name>
<name><surname>Sengamalay</surname>
<given-names>NA</given-names>
</name>
<name><surname>Shumway</surname>
<given-names>M</given-names>
</name>
<name><surname>Zaborsky</surname>
<given-names>J</given-names>
</name>
<name><surname>Feldblyum</surname>
<given-names>T</given-names>
</name>
<name><surname>Subbu</surname>
<given-names>V</given-names>
</name>
<name><surname>Spiro</surname>
<given-names>DJ</given-names>
</name>
<name><surname>Sitz</surname>
<given-names>J</given-names>
</name>
<name><surname>Koo</surname>
<given-names>H</given-names>
</name>
<name><surname>Bolotov</surname>
<given-names>P</given-names>
</name>
<etal></etal>
<article-title>Large-scale sequencing of human influenza reveals the dynamic nature of viral genome evolution</article-title>
<source>Nature</source>
<year>2005</year>
<volume>437</volume>
<issue>7062</issue>
<fpage>1162</fpage>
<lpage>1166</lpage>
<pub-id pub-id-type="doi">10.1038/nature04239</pub-id>
<pub-id pub-id-type="pmid">16208317</pub-id>
</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/H2N2V1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 0005870 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 0005870 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= H2N2V1 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.33. | ![]() |