Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Identifying Changes in Selective Constraints: Host Shifts in Influenza</title><author><name sortKey="Tamuri, Asif U" sort="Tamuri, Asif U" uniqKey="Tamuri A" first="Asif U." last="Tamuri">Asif U. Tamuri</name><affiliation><nlm:aff id="aff1"></nlm:aff></affiliation></author><author><name sortKey="Dos Reis, Mario" sort="Dos Reis, Mario" uniqKey="Dos Reis M" first="Mario" last="Dos Reis">Mario Dos Reis</name><affiliation><nlm:aff id="aff1"></nlm:aff></affiliation></author><author><name sortKey="Hay, Alan J" sort="Hay, Alan J" uniqKey="Hay A" first="Alan J." last="Hay">Alan J. Hay</name><affiliation><nlm:aff id="aff1"></nlm:aff></affiliation></author><author><name sortKey="Goldstein, Richard A" sort="Goldstein, Richard A" uniqKey="Goldstein R" first="Richard A." last="Goldstein">Richard A. Goldstein</name><affiliation><nlm:aff id="aff1"></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">19911053</idno><idno type="pmc">2770840</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2770840</idno><idno type="RBID">PMC:2770840</idno><idno type="doi">10.1371/journal.pcbi.1000564</idno><date when="2009">2009</date><idno type="wicri:Area/Pmc/Corpus">000A41</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000A41</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Identifying Changes in Selective Constraints: Host Shifts in Influenza</title><author><name sortKey="Tamuri, Asif U" sort="Tamuri, Asif U" uniqKey="Tamuri A" first="Asif U." last="Tamuri">Asif U. Tamuri</name><affiliation><nlm:aff id="aff1"></nlm:aff></affiliation></author><author><name sortKey="Dos Reis, Mario" sort="Dos Reis, Mario" uniqKey="Dos Reis M" first="Mario" last="Dos Reis">Mario Dos Reis</name><affiliation><nlm:aff id="aff1"></nlm:aff></affiliation></author><author><name sortKey="Hay, Alan J" sort="Hay, Alan J" uniqKey="Hay A" first="Alan J." last="Hay">Alan J. Hay</name><affiliation><nlm:aff id="aff1"></nlm:aff></affiliation></author><author><name sortKey="Goldstein, Richard A" sort="Goldstein, Richard A" uniqKey="Goldstein R" first="Richard A." last="Goldstein">Richard A. Goldstein</name><affiliation><nlm:aff id="aff1"></nlm:aff></affiliation></author></analytic><series><title level="j">PLoS Computational Biology</title><idno type="ISSN">1553-734X</idno><idno type="eISSN">1553-7358</idno><imprint><date when="2009">2009</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>The natural reservoir of Influenza A is waterfowl. Normally, waterfowl viruses are not adapted to infect and spread in the human population. Sometimes, through reassortment or through whole host shift events, genetic material from waterfowl viruses is introduced into the human population causing worldwide pandemics. Identifying which mutations allow viruses from avian origin to spread successfully in the human population is of great importance in predicting and controlling influenza pandemics. Here we describe a novel approach to identify such mutations. We use a sitewise non-homogeneous phylogenetic model that explicitly takes into account differences in the equilibrium frequencies of amino acids in different hosts and locations. We identify 172 amino acid sites with strong support and 518 sites with moderate support of different selection constraints in human and avian viruses. The sites that we identify provide an invaluable resource to experimental virologists studying adaptation of avian flu viruses to the human host. Identification of the sequence changes necessary for host shifts would help us predict the pandemic potential of various strains. The method is of broad applicability to investigating changes in selective constraints when the timing of the changes is known.</p></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="Schafer, Jr" uniqKey="Schafer J">JR Schafer</name></author><author><name sortKey="Kawaoka, Y" uniqKey="Kawaoka Y">Y Kawaoka</name></author><author><name sortKey="Bean, Wj" uniqKey="Bean W">WJ Bean</name></author><author><name sortKey="Suss, J" uniqKey="Suss J">J Suss</name></author><author><name sortKey="Senne, D" uniqKey="Senne D">D Senne</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Reid, Ah" uniqKey="Reid A">AH Reid</name></author><author><name sortKey="Taubenberger, Jk" uniqKey="Taubenberger J">JK Taubenberger</name></author><author><name sortKey="Fanning, Tg" uniqKey="Fanning T">TG Fanning</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Taubenberger, Jk" uniqKey="Taubenberger J">JK Taubenberger</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Taubenberger, Jk" uniqKey="Taubenberger J">JK Taubenberger</name></author><author><name sortKey="Reid, Ah" uniqKey="Reid A">AH Reid</name></author><author><name sortKey="Lourens, Rm" uniqKey="Lourens R">RM Lourens</name></author><author><name sortKey="Wang, R" uniqKey="Wang R">R Wang</name></author><author><name sortKey="Jin, G" uniqKey="Jin G">G Jin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Antonovics, J" uniqKey="Antonovics J">J Antonovics</name></author><author><name sortKey="Hood, Me" uniqKey="Hood M">ME Hood</name></author><author><name sortKey="Baker, Ch" uniqKey="Baker C">CH Baker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dos Reis, M" uniqKey="Dos Reis M">M dos Reis</name></author><author><name sortKey="Hay, Aj" uniqKey="Hay A">AJ Hay</name></author><author><name sortKey="Goldstein, Ra" uniqKey="Goldstein R">RA Goldstein</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gibbs, M" uniqKey="Gibbs M">M Gibbs</name></author><author><name sortKey="Gibbs, A" uniqKey="Gibbs A">A Gibbs</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Smith, Gjd" uniqKey="Smith G">GJD Smith</name></author><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="Zhang, Jz" uniqKey="Zhang J">JZ Zhang</name></author><author><name sortKey="Poon, Llm" uniqKey="Poon L">LLM Poon</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lin, Yp" uniqKey="Lin Y">YP Lin</name></author><author><name sortKey="Shaw, M" uniqKey="Shaw M">M Shaw</name></author><author><name sortKey="Gregory, V" uniqKey="Gregory V">V Gregory</name></author><author><name sortKey="Cameron, K" uniqKey="Cameron K">K Cameron</name></author><author><name sortKey="Lim, W" uniqKey="Lim W">W Lim</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gambaryan, As" uniqKey="Gambaryan A">AS Gambaryan</name></author><author><name sortKey="Matrosovich, Mn" uniqKey="Matrosovich M">MN Matrosovich</name></author><author><name sortKey="Yamnikova, Ss" uniqKey="Yamnikova S">SS Yamnikova</name></author><author><name sortKey="Lvov, Dk" uniqKey="Lvov D">DK Lvov</name></author><author><name sortKey="Robertson, Js" uniqKey="Robertson J">JS Robertson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Connor, Rj" uniqKey="Connor R">RJ Connor</name></author><author><name sortKey="Kawaoka, Y" uniqKey="Kawaoka Y">Y Kawaoka</name></author><author><name sortKey="Webster, Rg" uniqKey="Webster R">RG Webster</name></author><author><name sortKey="Paulson, Jc" uniqKey="Paulson J">JC Paulson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Matrosovich, M" uniqKey="Matrosovich M">M Matrosovich</name></author><author><name sortKey="Tuzikov, A" uniqKey="Tuzikov A">A Tuzikov</name></author><author><name sortKey="Bovin, N" uniqKey="Bovin N">N Bovin</name></author><author><name sortKey="Gambaryan, A" uniqKey="Gambaryan A">A Gambaryan</name></author><author><name sortKey="Klimov, A" uniqKey="Klimov A">A Klimov</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nobusawa, E" uniqKey="Nobusawa E">E Nobusawa</name></author><author><name sortKey="Aoyama, T" uniqKey="Aoyama T">T Aoyama</name></author><author><name sortKey="Kato, H" uniqKey="Kato H">H Kato</name></author><author><name sortKey="Suzuki, Y" uniqKey="Suzuki Y">Y Suzuki</name></author><author><name sortKey="Tateno, Y" uniqKey="Tateno Y">Y Tateno</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rogers, Gn" uniqKey="Rogers G">GN Rogers</name></author><author><name sortKey="Paulson, Jc" uniqKey="Paulson J">JC Paulson</name></author><author><name sortKey="Daniels, Rs" uniqKey="Daniels R">RS Daniels</name></author><author><name sortKey="Skehel, Jj" uniqKey="Skehel J">JJ Skehel</name></author><author><name sortKey="Wilson, Ia" uniqKey="Wilson I">IA Wilson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vines, A" uniqKey="Vines A">A Vines</name></author><author><name sortKey="Wells, K" uniqKey="Wells K">K Wells</name></author><author><name sortKey="Matrosovich, M" uniqKey="Matrosovich M">M Matrosovich</name></author><author><name sortKey="Castrucci, Mr" uniqKey="Castrucci M">MR Castrucci</name></author><author><name sortKey="Ito, T" uniqKey="Ito T">T Ito</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Baigent, Sj" uniqKey="Baigent S">SJ Baigent</name></author><author><name sortKey="Mccauley, Jw" uniqKey="Mccauley J">JW McCauley</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Naffakh, N" uniqKey="Naffakh N">N Naffakh</name></author><author><name sortKey="Tomoiu, A" uniqKey="Tomoiu A">A Tomoiu</name></author><author><name sortKey="Rameix Welti, Ma" uniqKey="Rameix Welti M">MA Rameix-Welti</name></author><author><name sortKey="Van Der Werf, S" uniqKey="Van Der Werf S">S van der Werf</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hatta, M" uniqKey="Hatta M">M Hatta</name></author><author><name sortKey="Gao, P" uniqKey="Gao P">P Gao</name></author><author><name sortKey="Halfmann, P" uniqKey="Halfmann P">P Halfmann</name></author><author><name sortKey="Kawaoka, Y" uniqKey="Kawaoka Y">Y Kawaoka</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Steel, J" uniqKey="Steel J">J Steel</name></author><author><name sortKey="Lowen, A" uniqKey="Lowen A">A Lowen</name></author><author><name sortKey="Mubareka, S" uniqKey="Mubareka S">S Mubareka</name></author><author><name sortKey="Palese, P" uniqKey="Palese P">P Palese</name></author><author><name sortKey="Baric, R" uniqKey="Baric R">R Baric</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Subbarao, Ek" uniqKey="Subbarao E">EK Subbarao</name></author><author><name sortKey="London, W" uniqKey="London W">W London</name></author><author><name sortKey="Murphy, Br" uniqKey="Murphy B">BR Murphy</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chen, Gw" uniqKey="Chen G">GW Chen</name></author><author><name sortKey="Chang, Sc" uniqKey="Chang S">SC Chang</name></author><author><name sortKey="Mok, Ck" uniqKey="Mok C">CK Mok</name></author><author><name sortKey="Lo, Yl" uniqKey="Lo Y">YL Lo</name></author><author><name sortKey="Kung, Yn" uniqKey="Kung Y">YN Kung</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Finkelstein, D" uniqKey="Finkelstein D">D Finkelstein</name></author><author><name sortKey="Mukatira, S" uniqKey="Mukatira S">S Mukatira</name></author><author><name sortKey="Mehta, P" uniqKey="Mehta P">P Mehta</name></author><author><name sortKey="Obenauer, J" uniqKey="Obenauer J">J Obenauer</name></author><author><name sortKey="Su, X" uniqKey="Su X">X Su</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Miotto, O" uniqKey="Miotto O">O Miotto</name></author><author><name sortKey="Heiny, A" uniqKey="Heiny A">A Heiny</name></author><author><name sortKey="Tan, T" uniqKey="Tan T">T Tan</name></author><author><name sortKey="August, J" uniqKey="August J">J August</name></author><author><name sortKey="Brusic, V" uniqKey="Brusic V">V Brusic</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mayr, E" uniqKey="Mayr E">E Mayr</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bruno, W" uniqKey="Bruno W">W Bruno</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Freire Maia, N" uniqKey="Freire Maia N">N Freire-Maia</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Blouin, C" uniqKey="Blouin C">C Blouin</name></author><author><name sortKey="Boucher, Y" uniqKey="Boucher Y">Y Boucher</name></author><author><name sortKey="Roger, Aj" uniqKey="Roger A">AJ Roger</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dorman, K" uniqKey="Dorman K">K Dorman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gu, X" uniqKey="Gu X">X Gu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gu, X" uniqKey="Gu X">X Gu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gu, X" uniqKey="Gu X">X Gu</name></author><author><name sortKey="Fu, Y X" uniqKey="Fu Y">Y-X Fu</name></author><author><name sortKey="Li, W H" uniqKey="Li W">W-H Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Knudsen, B" uniqKey="Knudsen B">B Knudsen</name></author><author><name sortKey="Miyamoto, Mm" uniqKey="Miyamoto M">MM Miyamoto</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Penn, O" uniqKey="Penn O">O Penn</name></author><author><name sortKey="Stern, A" uniqKey="Stern A">A Stern</name></author><author><name sortKey="Rubinstein, N" uniqKey="Rubinstein N">N Rubinstein</name></author><author><name sortKey="Dutheil, J" uniqKey="Dutheil J">J Dutheil</name></author><author><name sortKey="Bacharach, E" uniqKey="Bacharach E">E Bacharach</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pupko, T" uniqKey="Pupko T">T Pupko</name></author><author><name sortKey="Galtier, N" uniqKey="Galtier N">N Galtier</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Guindon, S" uniqKey="Guindon S">S Guindon</name></author><author><name sortKey="Rodrigo, Ag" uniqKey="Rodrigo A">AG Rodrigo</name></author><author><name sortKey="Dyer, Ka" uniqKey="Dyer K">KA Dyer</name></author><author><name sortKey="Huelsenbeck, Jp" uniqKey="Huelsenbeck J">JP Huelsenbeck</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, Zh" uniqKey="Yang Z">ZH Yang</name></author><author><name sortKey="Nielsen, R" uniqKey="Nielsen R">R Nielsen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, Jz" uniqKey="Zhang J">JZ Zhang</name></author><author><name sortKey="Nielsen, R" uniqKey="Nielsen R">R Nielsen</name></author><author><name sortKey="Yang, Zh" uniqKey="Yang Z">ZH Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Forsberg, R" uniqKey="Forsberg R">R Forsberg</name></author><author><name sortKey="Christiansen, Fb" uniqKey="Christiansen F">FB Christiansen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Blackburne, Bp" uniqKey="Blackburne B">BP Blackburne</name></author><author><name sortKey="Hay, Aj" uniqKey="Hay A">AJ Hay</name></author><author><name sortKey="Goldstein, Ra" uniqKey="Goldstein R">RA Goldstein</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Whelan, S" uniqKey="Whelan S">S Whelan</name></author><author><name sortKey="Goldman, N" uniqKey="Goldman N">N Goldman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Benjamini, Y" uniqKey="Benjamini Y">Y Benjamini</name></author><author><name sortKey="Hochberg, T" uniqKey="Hochberg T">T Hochberg</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Crow, Jf" uniqKey="Crow J">JF Crow</name></author><author><name sortKey="Kimura, M" uniqKey="Kimura M">M Kimura</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dawood, Fs" uniqKey="Dawood F">FS Dawood</name></author><author><name sortKey="Jain, S" uniqKey="Jain S">S Jain</name></author><author><name sortKey="Finelli, L" uniqKey="Finelli L">L Finelli</name></author><author><name sortKey="Shaw, Mw" uniqKey="Shaw M">MW Shaw</name></author><author><name sortKey="Lindstrom, Se" uniqKey="Lindstrom S">SE Lindstrom</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sheerar, M" uniqKey="Sheerar M">M Sheerar</name></author><author><name sortKey="Easterday, B" uniqKey="Easterday B">B Easterday</name></author><author><name sortKey="Hinshaw, V" uniqKey="Hinshaw V">V Hinshaw</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vincent, A" uniqKey="Vincent A">A Vincent</name></author><author><name sortKey="Lager, K" uniqKey="Lager K">K Lager</name></author><author><name sortKey="Ma, W" uniqKey="Ma W">W Ma</name></author><author><name sortKey="Lekcharoensuk, P" uniqKey="Lekcharoensuk P">P Lekcharoensuk</name></author><author><name sortKey="Gramer, M" uniqKey="Gramer M">M Gramer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Boni, Mf" uniqKey="Boni M">MF Boni</name></author><author><name sortKey="Zhou, Y" uniqKey="Zhou Y">Y Zhou</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="Chare, Er" uniqKey="Chare E">ER Chare</name></author><author><name sortKey="Gould, Ea" uniqKey="Gould E">EA Gould</name></author><author><name sortKey="Holmes, Ec" uniqKey="Holmes E">EC Holmes</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Krasnitz, M" uniqKey="Krasnitz M">M Krasnitz</name></author><author><name sortKey="Levine, Aj" uniqKey="Levine A">AJ Levine</name></author><author><name sortKey="Rabadan, R" uniqKey="Rabadan R">R Rabadan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, Z" uniqKey="Yang Z">Z Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Felsenstein, J" uniqKey="Felsenstein J">J Felsenstein</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Felsenstein, J" uniqKey="Felsenstein J">J Felsenstein</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bao, Ym" uniqKey="Bao Y">YM Bao</name></author><author><name sortKey="Bolotov, P" uniqKey="Bolotov P">P Bolotov</name></author><author><name sortKey="Dernovoy, D" uniqKey="Dernovoy D">D Dernovoy</name></author><author><name sortKey="Kiryutin, B" uniqKey="Kiryutin B">B Kiryutin</name></author><author><name sortKey="Zaslavsky, L" uniqKey="Zaslavsky L">L Zaslavsky</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="Suyama, M" uniqKey="Suyama M">M Suyama</name></author><author><name sortKey="Torrents, D" uniqKey="Torrents D">D Torrents</name></author><author><name sortKey="Bork, P" uniqKey="Bork P">P Bork</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Guindon, S" uniqKey="Guindon S">S Guindon</name></author><author><name sortKey="Gascuel, O" uniqKey="Gascuel O">O Gascuel</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, Ta" uniqKey="Yano T">TA Yano</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, Z" uniqKey="Yang Z">Z Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, Z" uniqKey="Yang Z">Z Yang</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS Comput Biol</journal-id><journal-id journal-id-type="iso-abbrev">PLoS Comput. Biol</journal-id><journal-id journal-id-type="publisher-id">plos</journal-id><journal-id journal-id-type="pmc">ploscomp</journal-id><journal-title-group><journal-title>PLoS Computational Biology</journal-title></journal-title-group><issn pub-type="ppub">1553-734X</issn><issn pub-type="epub">1553-7358</issn><publisher><publisher-name>Public Library of Science</publisher-name><publisher-loc>San Francisco, USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">19911053</article-id><article-id pub-id-type="pmc">2770840</article-id><article-id pub-id-type="publisher-id">09-PLCB-RA-0610R2</article-id><article-id pub-id-type="doi">10.1371/journal.pcbi.1000564</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline"><subject>Computational Biology/Comparative Sequence Analysis</subject><subject>Computational Biology/Evolutionary Modeling</subject><subject>Evolutionary Biology/Microbial Evolution and Genomics</subject><subject>Infectious Diseases/Viral Infections</subject><subject>Molecular Biology/Molecular Evolution</subject><subject>Virology/Virus Evolution and Symbiosis</subject></subj-group></article-categories><title-group><article-title>Identifying Changes in Selective Constraints: Host Shifts in Influenza</article-title><alt-title alt-title-type="running-head">Changes in Selective Constraints: Influenza</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Tamuri</surname><given-names>Asif U.</given-names></name><xref ref-type="aff" rid="aff1"></xref></contrib><contrib contrib-type="author"><name><surname>dos Reis</surname><given-names>Mario</given-names></name><xref ref-type="aff" rid="aff1"></xref></contrib><contrib contrib-type="author"><name><surname>Hay</surname><given-names>Alan J.</given-names></name><xref ref-type="aff" rid="aff1"></xref></contrib><contrib contrib-type="author"><name><surname>Goldstein</surname><given-names>Richard A.</given-names></name><xref ref-type="aff" rid="aff1"></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>National Institute for Medical Research, London, United Kingdom</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>Fraser</surname><given-names>Christophe</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1">Imperial College London, United Kingdom</aff><author-notes><corresp id="cor1">* E-mail: <email>Richard.Goldstein@nimr.mrc.ac.uk</email></corresp><fn fn-type="con"><p>Conceived and designed the experiments: MdR AJH RAG. Performed the experiments: AUT RAG. Analyzed the data: AUT MdR RAG. Contributed reagents/materials/analysis tools: MdR. Wrote the paper: MdR AJH RAG.</p></fn></author-notes><pub-date pub-type="collection"><month>11</month><year>2009</year></pub-date><pub-date pub-type="epub"><day>13</day><month>11</month><year>2009</year></pub-date><volume>5</volume><issue>11</issue><elocation-id>e1000564</elocation-id><history><date date-type="received"><day>27</day><month>5</month><year>2009</year></date><date date-type="accepted"><day>15</day><month>10</month><year>2009</year></date></history><permissions><copyright-statement>Tamuri et al.</copyright-statement><copyright-year>2009</copyright-year><license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are properly credited.</license-p></license></permissions><abstract><p>The natural reservoir of Influenza A is waterfowl. Normally, waterfowl viruses are not adapted to infect and spread in the human population. Sometimes, through reassortment or through whole host shift events, genetic material from waterfowl viruses is introduced into the human population causing worldwide pandemics. Identifying which mutations allow viruses from avian origin to spread successfully in the human population is of great importance in predicting and controlling influenza pandemics. Here we describe a novel approach to identify such mutations. We use a sitewise non-homogeneous phylogenetic model that explicitly takes into account differences in the equilibrium frequencies of amino acids in different hosts and locations. We identify 172 amino acid sites with strong support and 518 sites with moderate support of different selection constraints in human and avian viruses. The sites that we identify provide an invaluable resource to experimental virologists studying adaptation of avian flu viruses to the human host. Identification of the sequence changes necessary for host shifts would help us predict the pandemic potential of various strains. The method is of broad applicability to investigating changes in selective constraints when the timing of the changes is known.</p></abstract><abstract abstract-type="summary"><title>Author Summary</title><p>Influenza A's natural reservoir is waterfowl. Sometimes avian virus genomic segments are able to shift to a human host, either <italic>in toto</italic> or by combining with those that underwent a previous host shift event. Such host shift events can cause worldwide pandemics in their immunologically naive hosts. In order for these host shifts to establish a stable lineage, the virus has to adapt to the new host. Identifying the changes that have occurred in the past can provide important clues about how this process happens, and how surveillance for new influenza threats should be targeted. Unfortunately, it is difficult to determine whether an amino acid has changed due to adaptation to the new host or whether the change occurred through random drift. Here we describe a novel phylogenetic approach to identifying locations where the nature of the selective pressure exerted on the location has changed corresponding to the host shift event. We identify a set of locations on a number of the genomic segments. The approach we describe is of wide applicability when the timing of the change of selective constraints is known in advance.</p></abstract><counts><page-count count="14"></page-count></counts></article-meta></front><body><sec id="s1"><title>Introduction</title><p>Influenza A has the distinction of being an old disease, a recurring disease, and an ‘emerging’ disease. Influenza A viruses are found in humans as well as in other animals including swine, horses, sea mammals, and birds, of which waterfowl are considered the natural reservoir <xref rid="pcbi.1000564-Webster1" ref-type="bibr">[1]</xref>. Subtypes of influenza A are distinguished by two surface glycoproteins; haemagglutinin (HA), the primary target of the immune response, and neuraminidase (NA). There are sixteen known types of haemagglutinin (H1 to H16) and nine of neuraminidase (N1 to N9), all found in waterfowl. Only H1, H2, H3 and N1, N2, however, are known to have caused epidemic disease in humans. The predominant forms of influenza A currently circulating in humans are H1N1 and H3N2.</p><p>There are two distinct problems represented by influenza. Firstly, the various subtypes currently in circulation in humans cause significant morbidity and loss of life. Secondly, periodically a subtype of influenza can make the shift from aquatic birds to humans, possibly through an intermediate host, resulting in a widespread pandemic in an immunologically-naïve population. These ‘antigenic shifts’ can occur either through the transfer of an entire virus from one host to another, or through a re-assortment process where genomic segments of the avian virus mix with genomic segments currently circulating in humans. In 1957 three virus segments (HA, NA, and PB1) from an avian-like source were combined with the other five segments already circulating in humans to create the H2N2 ‘Asian flu’ pandemic, while in 1968 two segments (HA and PB1) from an avian-like source were combined with the other six from the already-present human H2N2 virus to form the H3N2 ‘Hong Kong flu’ pandemic <xref rid="pcbi.1000564-Schafer1" ref-type="bibr">[2]</xref>. It has been suggested that the 1918 H1N1 ‘Spanish flu’ virus was the result of a single host-shift event from birds to humans <xref rid="pcbi.1000564-Reid1" ref-type="bibr">[3]</xref>–<xref rid="pcbi.1000564-Taubenberger2" ref-type="bibr">[5]</xref> but this remains controversial <xref rid="pcbi.1000564-Antonovics1" ref-type="bibr">[6]</xref>–<xref rid="pcbi.1000564-Smith1" ref-type="bibr">[9]</xref>. In recent years, a number of different avian subtypes have caused sporadic human infections, including H5N1, H7N3, H7N7, and H9N2 <xref rid="pcbi.1000564-Lin1" ref-type="bibr">[10]</xref>. While there is evidence for sporadic transmissions of these avian viruses between humans, the genetic changes necessary for widespread human-human transmission have, so far, seemingly not occurred.</p><p>A number of different proteins have been implicated in determining host ranges. Influenza haemagglutinin binds to sialic acid linked to galactose on the surface of the targeted cell; the differing nature of the sialic acid-galactose linkages in birds and humans (α2,3 sialic acid linkages in the bird gut, α2,6 sialic acid linkages of the upper human respiratory tract <xref rid="pcbi.1000564-Gambaryan1" ref-type="bibr">[11]</xref>) provides an important barrier to host shift events. A number of amino acid substitutions have occurred in human influenza haemagglutinin (e.g. Q226L and G228S in H2 and H3, E190N/D and G225E/D in H1) to adjust to the different receptors <xref rid="pcbi.1000564-Connor1" ref-type="bibr">[12]</xref>–<xref rid="pcbi.1000564-Vines1" ref-type="bibr">[16]</xref>. Neuraminidase, the protein responsible for cleaving the haemagglutinin from the receptor surface, also seems adapted to the particular sialic acid linkages, as well as for the pH and temperature of the host tissues <xref rid="pcbi.1000564-Baigent1" ref-type="bibr">[17]</xref>. Proteins in the viral replication complex (PA, PB1, PB2, and NP) have also been implicated in limiting host range by restricting replication and intra-host spread in mammals (for a review, see <xref rid="pcbi.1000564-Naffakh1" ref-type="bibr">[18]</xref>.) Of particular note is the PB2 gene, where one specific substitution, E627K, was identified and characterised experimentally as crucial for replication and intra-host spread in mammals <xref rid="pcbi.1000564-Hatta1" ref-type="bibr">[19]</xref>–<xref rid="pcbi.1000564-Subbarao1" ref-type="bibr">[21]</xref>.</p><p>As part of the widespread surveillance effort, it is important to understand the process of host shifts, and to identify the important changes that are necessary for the shift to occur, or that make the shift more likely. We currently have many examples of both avian and human viruses, so there have been a number of efforts at identifying ‘genetic signatures’ that characterise the virus as adapted to one or the other host. The most common method is to identify sites where the distribution of amino acids found in the virus in one host are sufficiently different from the distribution of amino acids found in the same site in viruses that affect the other host <xref rid="pcbi.1000564-Chen1" ref-type="bibr">[22]</xref>–<xref rid="pcbi.1000564-Miotto1" ref-type="bibr">[24]</xref>. Unfortunately, there are two fundamental problems with this approach.</p><p>Firstly, the observed changes could represent the result of neutral drift rather than anything specific to the nature of the different hosts. As the human viruses are more closely related to each other than they are to the avian viruses, it would be expected that there would be characteristic amino acids found in the human lineages that are distinct from those found in the avian lineages because of the ‘founder effect’ <xref rid="pcbi.1000564-Mayr1" ref-type="bibr">[25]</xref>, that is, the maintenance of the idiosyncratic properties of the particular virus that first infected humans. Comparisons of amino acid frequencies in viruses from the two hosts cannot easily distinguish between those that accidentally accompanied the host shift event and those that were actually associated with different selective constraints acting on the viruses in the two hosts.</p><p>The second related problem is the use of inappropriate statistical tests to identify when these two distributions are sufficiently different. The statistical tests used generally assume that each of the observed sequences represent a set of independent measurements. But the underlying phylogenetic relationships will generate correlations in the amino acids at a site, confounding the signal due to the host shift event. This can be demonstrated by considering <xref ref-type="fig" rid="pcbi-1000564-g001">Figure 1</xref>, which shows two possible situations where the avian viruses all have a leucine in a given position where all of the human viruses have a valine in the same position. In example A the results are statistically significant, in that the positions are independent, and it is unlikely that the simultaneous parallel changes in sequence occurred at random in the human viruses but not in the avian viruses. In example B there is much less statistical signal, as only one change of amino acid on the branch connecting the human and avian viruses is needed to explain the multiple observations. By neglecting the underlying phylogenetic structure, a single change of amino acid can be interpreted as a large number of independent events, grossly exaggerating the statistical significance.</p><fig id="pcbi-1000564-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000564.g001</object-id><label>Figure 1</label><caption><title>Possible evolutionary scenarios.</title><p>Two possible phylogenetic trees representing the situation where eight different avian sequences have a L in a given position, while eight different human sequences have a V. (Branch lengths are not drawn to scale.) The situation shown on the right provides much weaker evidence for a shift in selective constraints compared with the situation shown on the left.</p></caption><graphic xlink:href="pcbi.1000564.g001"></graphic></fig><p>A number of the published approaches to this problem suffer from the above problems. For example, both Chen et al. <xref rid="pcbi.1000564-Chen1" ref-type="bibr">[22]</xref> and Miotto et al. <xref rid="pcbi.1000564-Miotto1" ref-type="bibr">[24]</xref> employed an information-based approach to identify sites where host-specific amino acids can be identified. Their computations of entropy (a measure of sequence diversity) and mutual information (the dependence of the observed residue distribution on host species) are based on considering every observed sequence as an independent data-point, ignoring correlations between the evolutionarily related sequences. Different distributions in the two hosts can be explained due to the founder effect described above, independent of any role these sites have in host adaptation. That is not to say that their results are incorrect, only that these problems make it impossible to determine their statistical significance. Finkelstein et al <xref rid="pcbi.1000564-Finkelstein1" ref-type="bibr">[23]</xref> looked at sites with a significantly higher degree of conservation in human lineages than avian lineages, and identified 32 markers within the M1, NP, NS, PA, and PB2 genes, 26 of them on the polymerase proteins NP, PA, and PB2. This analysis did not consider the phylogenetic relationships explicitly in their calculation of conservation, choosing instead to base their calculation on the frequency of the different amino acids observed in that site in the different hosts. While they employed strict tests for, for instance, multiple hypothesis testing, it is difficult to determine how much their results were affected by considering only frequencies of amino acids to represent the selective constraints, again ignoring the underlying phylogenetic relationships. It is known, for instance, that such counting methods produce very inaccurate amino acid frequencies compared with phylogenetically-based methods <xref rid="pcbi.1000564-Bruno1" ref-type="bibr">[26]</xref>, and can not generally identify the rate of substitutions in the tree, but only the range of acceptable amino acids.</p><p>As described above, the differences in the distribution of amino acids at a given site between avian and human viruses might represent neutral drift or, more interestingly, a change in the underlying selective pressure applied to the virus by the host. Rather than characterising only the difference in observed amino acid distributions, we can instead look directly for evidence of changes in the selective constraints by modelling the phylogenetics explicitly. These selective constraint changes will result in differences in the substitution process, as mutations that arise in one virus or another will have different probabilities of achieving fixation. Thus, changes in selection constraints will manifest themselves as changes in the observed substitution rates. This also allows rigorous statistical methods, such as the likelihood ratio test, to be used to establish statistical significance.</p><p>The selective pressure acting on a site can be positive, negative, or neutral. Positive selection, also called adaptive, or more misleadingly <xref rid="pcbi.1000564-FreireMaia1" ref-type="bibr">[27]</xref> ‘Darwinian’, refers to the acceptance of advantageous mutations; negative, or purifying selection involves the rejection of deleterious mutations. Neutral selection pressure involves the chance acceptance of mutations that do not have a significant effect on the fitness. Both positive and negative selection pressure represent strong constraints on the amino acids at a given site; the difference is that during purifying selection the current amino acids generally fulfil these constraints so change is restricted, while during adaptive evolution the current amino acids are not well suited, generally due to changes in the constraints or a selective advantage for diversification, enhancing the rate of evolution until more appropriate residues are found.</p><p>Changes in the selective constraints can result in changes in the rate of substitutions at that location. If the initial amino acids do not match the current requirements of that site, there may be an adaptive burst of faster substitutions until the constraints are satisfied. Modifications of the stringency of the constraints, causing a given site to be more or less restricted, may cause a longer-term change in the substitution rate without necessarily causing an adaptive burst. Previous phylogenetic methods have generally focused on identifying changes in the absolute substitution rate <xref rid="pcbi.1000564-Blouin1" ref-type="bibr">[28]</xref>–<xref rid="pcbi.1000564-Pupko1" ref-type="bibr">[35]</xref> or ratio of non-synonymous to synonymous changes <xref rid="pcbi.1000564-Guindon1" ref-type="bibr">[36]</xref>–<xref rid="pcbi.1000564-Zhang1" ref-type="bibr">[38]</xref>. The latter method was used, for instance, to identify twelve sites on the influenza A nucleoprotein that seem to have undergone a change in selective constraints corresponding to the switch from avian to human host <xref rid="pcbi.1000564-Forsberg1" ref-type="bibr">[39]</xref>. While these approaches are often useful, transient position-specific adaptive bursts are difficult to identify given the short duration of the effect. Sites can also undergo shifts in selective constraints without adaptive bursts or detectible changes in substitution rates, especially if the constraints in the two hosts overlap. Monitoring changes in the <italic>nature</italic> of the selective constraints has been much less common <xref rid="pcbi.1000564-Blackburne1" ref-type="bibr">[40]</xref> and has not been applied to host shift events.</p><p>In this paper we investigate the use of a phylogenetic method to detect changes in selective constraints that considers not only changes in the magnitude of selection constraints, but also changes in its nature, represented as the relative propensity for the different amino acids. We do this by considering two different models for each site, a homogeneous model where the selective constraints are independent of host, the other a non-homogeneous model where the selective constraints depend upon the host. The likelihood ratio test can then determine the level of statistical support for rejecting the null hypothesis of no such dependence.</p></sec><sec id="s2"><title>Results</title><p>We start our analysis with a set of human and avian influenza viral sequences and the associated phylogenetic trees for each influenza gene. We consider the different haemagglutinin and neuraminidase serotypes (e.g. H1, H2, H3, N1, N2) separately. For each non-conserved site, we apply increasingly complicated substitution models, using the Likelihood Ratio Test (LRT) to evaluate the statistical support for each further complication.</p><p>The substitution models are defined by a symmetric exchangeability matrix <bold>S</bold>, the equilibrium frequencies of the twenty amino acids <bold>π</bold>, and a rate scaling parameter ν representing the relative substitution rate at that site compared with other sites. The simplest model, Model 1, consists of the WAG exchangeability matrix combined with the associated equilibrium frequencies for the different amino acids <xref rid="pcbi.1000564-Whelan1" ref-type="bibr">[41]</xref>, with one adjustable parameter per site representing the scaling factor ν. We then consider Model 2 where the equilibrium frequencies of the amino acids are optimised individually for each site <xref rid="pcbi.1000564-Bruno1" ref-type="bibr">[26]</xref>. The likelihood ratio test demonstrated that the use of site-specific equilibrium frequencies was justified for all sites (<italic>P</italic> values ranging from 0.028 to 9.4×10<sup>−27</sup>).</p><p>We then created a non-homogeneous model, Model 3 where virus substitutions are modelled by one set of substitution rates in the avian host, and by a different set of substitution rates in the human host, as illustrated in <xref ref-type="fig" rid="pcbi-1000564-g002">Figure 2</xref>. The two different substitution models shared the WAG exchangeability matrix <bold>S</bold> and a site-specific rate-scaling factor ν, but now the equilibrium amino acid frequencies were both host- and site-specific. We identified sites with statistical support for different substitution rates in the two hosts, using a false discovery rate (FDR) method to account for multiple hypothesis testing <xref rid="pcbi.1000564-Benjamini1" ref-type="bibr">[42]</xref>. We identified 172 sites with an FDR<0.05 (<italic>i.e.</italic> we would expect 5% of these sites to be false positives), and 518 sites with an FDR<0.20. We will refer to the 172 higher-confidence locations as ‘A sites’ and the remaining 346 lower-confidence locations as ‘B sites’.</p><fig id="pcbi-1000564-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000564.g002</object-id><label>Figure 2</label><caption><title>Homogeneous and non-homogeneous substitution models.</title><p>Illustrative phylogenetic trees showing set of avian and human influenza sequences. A: In the homogeneous models (Models 1 and 2), the same substitution rates are used throughout the tree. B: In the non-homogeneous models (Models 3 and 4) different substitution rates are used for the avian (red) and human (blue) lineages. The root of the tree is assumed to be inside the avian lineage. (Because the model is reversible within the avian clade, the exact location of the root within this clade does not affect the calculation.) The host shift event is assumed to occur at the midpoint of the branch connecting the common ancestor of the human strains with its parent.</p></caption><graphic xlink:href="pcbi.1000564.g002"></graphic></fig><p>We then considered if modelling differences in the equilibrium amino acid frequencies was adequate, or whether we should include host-dependent rate scaling factors as well. We implemented a more complicated model (Model 4) where the substitution rates were still defined with the WAG exchangeability matrix, but now both the equilibrium frequencies and the scaling factor ν were host- and site-dependent. Of the 2143 sites considered, few (37) had <italic>P</italic> values less than 0.05; after correcting for multiple hypothesis testing using the false discovery rate method, no site yielded any statistically significant improvement. The results described below will be based on Model 3 above.</p><p>The list of 172 ‘A’ sites (FDR<0.05) is shown in <xref ref-type="table" rid="pcbi-1000564-t001">Table 1</xref>. Sites were found on all of the genes considered. Supporting <xref ref-type="supplementary-material" rid="pcbi.1000564.s001">Table S1</xref> shows the list of the 518 ‘A’ and ‘B’ sites with FDR<0.20. Sites that have been identified experimentally are detected using this method, notably PB2 627. HA sites H1 190 and 225 and H3 228 are also identified. Sites H2 226 and 228 are significant at the weaker FDR<0.20 level, while H3 226 is not statistically significant.</p><table-wrap id="pcbi-1000564-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000564.t001</object-id><label>Table 1</label><caption><title>Sites identified as undergoing changes in selective pressure during host shifts from birds to humans.</title></caption><alternatives><graphic id="pcbi-1000564-t001-1" xlink:href="pcbi.1000564.t001"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Location</td><td align="left" rowspan="1" colspan="1">P</td><td align="left" rowspan="1" colspan="1">d</td><td align="left" rowspan="1" colspan="1">Avian residues</td><td align="left" rowspan="1" colspan="1">Human residues</td><td align="left" rowspan="1" colspan="1">Avian selective constraint strength (d)</td><td align="left" rowspan="1" colspan="1">Human selective constraint strength (d)</td><td align="left" rowspan="1" colspan="1">Cal-4</td></tr></thead><tbody><tr><td colspan="8" align="left" rowspan="1"><bold>H1</bold></td></tr><tr><td align="left" rowspan="1" colspan="1">−5</td><td align="left" rowspan="1" colspan="1">3.05E-04</td><td align="left" rowspan="1" colspan="1">2.59E-03</td><td align="left" rowspan="1" colspan="1">E</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">2.74</td><td align="left" rowspan="1" colspan="1">2.85</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">5.09E-05</td><td align="left" rowspan="1" colspan="1">7.86E-04</td><td align="left" rowspan="1" colspan="1">F</td><td align="left" rowspan="1" colspan="1">L</td><td align="left" rowspan="1" colspan="1">3.30</td><td align="left" rowspan="1" colspan="1">2.63</td><td align="left" rowspan="1" colspan="1">L</td></tr><tr><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">1.58E-03</td><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">V((A))</td><td align="left" rowspan="1" colspan="1">A</td><td align="left" rowspan="1" colspan="1">2.59</td><td align="left" rowspan="1" colspan="1">2.91</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">4.47E-03</td><td align="left" rowspan="1" colspan="1">0.021</td><td align="left" rowspan="1" colspan="1">L</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">2.63</td><td align="left" rowspan="1" colspan="1">2.71</td><td align="left" rowspan="1" colspan="1">A</td></tr><tr><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">9.36E-03</td><td align="left" rowspan="1" colspan="1">0.038</td><td align="left" rowspan="1" colspan="1">V((I))</td><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">2.68</td><td align="left" rowspan="1" colspan="1">2.63</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">54</td><td align="left" rowspan="1" colspan="1">4.90E-05</td><td align="left" rowspan="1" colspan="1">7.86E-04</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">2.79</td><td align="left" rowspan="1" colspan="1">2.89</td><td align="left" rowspan="1" colspan="1">R</td></tr><tr><td align="left" rowspan="1" colspan="1">63</td><td align="left" rowspan="1" colspan="1">5.80E-03</td><td align="left" rowspan="1" colspan="1">0.025</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">2.89</td><td align="left" rowspan="1" colspan="1">2.79</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">70</td><td align="left" rowspan="1" colspan="1">5.19E-03</td><td align="left" rowspan="1" colspan="1">0.023</td><td align="left" rowspan="1" colspan="1">L</td><td align="left" rowspan="1" colspan="1">I((V))</td><td align="left" rowspan="1" colspan="1">2.63</td><td align="left" rowspan="1" colspan="1">2.57</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">77</td><td align="left" rowspan="1" colspan="1">4.58E-03</td><td align="left" rowspan="1" colspan="1">0.021</td><td align="left" rowspan="1" colspan="1">D</td><td align="left" rowspan="1" colspan="1">E(G)</td><td align="left" rowspan="1" colspan="1">3.07</td><td align="left" rowspan="1" colspan="1">2.25</td><td align="left" rowspan="1" colspan="1">E</td></tr><tr><td align="left" rowspan="1" colspan="1">80</td><td align="left" rowspan="1" colspan="1">1.86E-03</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">S((P))</td><td align="left" rowspan="1" colspan="1">2.73</td><td align="left" rowspan="1" colspan="1">2.26</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td align="left" rowspan="1" colspan="1">91</td><td align="left" rowspan="1" colspan="1">2.57E-03</td><td align="left" rowspan="1" colspan="1">0.014</td><td align="left" rowspan="1" colspan="1">S(T)</td><td align="left" rowspan="1" colspan="1">P</td><td align="left" rowspan="1" colspan="1">2.06</td><td align="left" rowspan="1" colspan="1">3.13</td><td align="left" rowspan="1" colspan="1">P</td></tr><tr><td align="left" rowspan="1" colspan="1">120</td><td align="left" rowspan="1" colspan="1">7.69E-04</td><td align="left" rowspan="1" colspan="1">5.68E-03</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">R</td><td align="left" rowspan="1" colspan="1">2.86</td><td align="left" rowspan="1" colspan="1">2.82</td><td align="left" rowspan="1" colspan="1">R</td></tr><tr><td align="left" rowspan="1" colspan="1">138</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">A</td><td align="left" rowspan="1" colspan="1">S((A))</td><td align="left" rowspan="1" colspan="1">2.91</td><td align="left" rowspan="1" colspan="1">2.24</td><td align="left" rowspan="1" colspan="1">A</td></tr><tr><td align="left" rowspan="1" colspan="1">141</td><td align="left" rowspan="1" colspan="1">1.55E-04</td><td align="left" rowspan="1" colspan="1">1.48E-03</td><td align="left" rowspan="1" colspan="1">Y</td><td align="left" rowspan="1" colspan="1">H</td><td align="left" rowspan="1" colspan="1">3.50</td><td align="left" rowspan="1" colspan="1">4.01</td><td align="left" rowspan="1" colspan="1">H</td></tr><tr><td align="left" rowspan="1" colspan="1">154</td><td align="left" rowspan="1" colspan="1">4.62E-05</td><td align="left" rowspan="1" colspan="1">7.86E-04</td><td align="left" rowspan="1" colspan="1">I((L))</td><td align="left" rowspan="1" colspan="1">L</td><td align="left" rowspan="1" colspan="1">2.51</td><td align="left" rowspan="1" colspan="1">2.63</td><td align="left" rowspan="1" colspan="1">L</td></tr><tr><td align="left" rowspan="1" colspan="1">155</td><td align="left" rowspan="1" colspan="1">1.01E-05</td><td align="left" rowspan="1" colspan="1">7.42E-04</td><td align="left" rowspan="1" colspan="1">T(I)</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">2.20</td><td align="left" rowspan="1" colspan="1">2.73</td><td align="left" rowspan="1" colspan="1">V</td></tr><tr><td align="left" rowspan="1" colspan="1">159</td><td align="left" rowspan="1" colspan="1">1.57E-04</td><td align="left" rowspan="1" colspan="1">1.48E-03</td><td align="left" rowspan="1" colspan="1">N(T)</td><td align="left" rowspan="1" colspan="1">G((S))</td><td align="left" rowspan="1" colspan="1">2.41</td><td align="left" rowspan="1" colspan="1">2.49</td><td align="left" rowspan="1" colspan="1">N</td></tr><tr><td align="left" rowspan="1" colspan="1">160</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">S</td><td align="left" rowspan="1" colspan="1">L((S))</td><td align="left" rowspan="1" colspan="1">2.50</td><td align="left" rowspan="1" colspan="1">2.41</td><td align="left" rowspan="1" colspan="1">S</td></tr><tr><td align="left" rowspan="1" colspan="1">163</td><td align="left" rowspan="1" colspan="1">1.89E-04</td><td align="left" rowspan="1" colspan="1">1.70E-03</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">N((T,S))</td><td align="left" rowspan="1" colspan="1">2.89</td><td align="left" rowspan="1" colspan="1">2.52</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">187</td><td align="left" rowspan="1" colspan="1">3.47E-03</td><td align="left" rowspan="1" colspan="1">0.017</td><td align="left" rowspan="1" colspan="1">T((N))</td><td align="left" rowspan="1" colspan="1">N((S))</td><td align="left" rowspan="1" colspan="1">2.61</td><td align="left" rowspan="1" colspan="1">2.71</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td align="left" rowspan="1" colspan="1">188</td><td align="left" rowspan="1" colspan="1">4.16E-05</td><td align="left" rowspan="1" colspan="1">7.86E-04</td><td align="left" rowspan="1" colspan="1">T((V,A))</td><td align="left" rowspan="1" colspan="1">I((S,T,M))</td><td align="left" rowspan="1" colspan="1">2.23</td><td align="left" rowspan="1" colspan="1">2.20</td><td align="left" rowspan="1" colspan="1">S</td></tr><tr><td align="left" rowspan="1" colspan="1">189</td><td align="left" rowspan="1" colspan="1">3.84E-04</td><td align="left" rowspan="1" colspan="1">2.96E-03</td><td align="left" rowspan="1" colspan="1">S(G)((D,N))</td><td align="left" rowspan="1" colspan="1">G(K,R)((T,E,D))</td><td align="left" rowspan="1" colspan="1">1.64</td><td align="left" rowspan="1" colspan="1">1.36</td><td align="left" rowspan="1" colspan="1">A</td></tr><tr><td align="left" rowspan="1" colspan="1">190</td><td align="left" rowspan="1" colspan="1">1.76E-09</td><td align="left" rowspan="1" colspan="1">2.99E-07</td><td align="left" rowspan="1" colspan="1">E</td><td align="left" rowspan="1" colspan="1">D(V)((N))</td><td align="left" rowspan="1" colspan="1">2.74</td><td align="left" rowspan="1" colspan="1">2.20</td><td align="left" rowspan="1" colspan="1">D</td></tr><tr><td align="left" rowspan="1" colspan="1">192</td><td align="left" rowspan="1" colspan="1">0.012</td><td align="left" rowspan="1" colspan="1">0.044</td><td align="left" rowspan="1" colspan="1">Q</td><td align="left" rowspan="1" colspan="1">(K,M,R)</td><td align="left" rowspan="1" colspan="1">3.34</td><td align="left" rowspan="1" colspan="1">1.98</td><td align="left" rowspan="1" colspan="1">Q</td></tr><tr><td align="left" rowspan="1" colspan="1">193</td><td align="left" rowspan="1" colspan="1">4.10E-03</td><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">N(E)((T,S))</td><td align="left" rowspan="1" colspan="1">(A,T,N)</td><td align="left" rowspan="1" colspan="1">1.80</td><td align="left" rowspan="1" colspan="1">1.84</td><td align="left" rowspan="1" colspan="1">S</td></tr><tr><td align="left" rowspan="1" colspan="1">197</td><td align="left" rowspan="1" colspan="1">5.62E-03</td><td align="left" rowspan="1" colspan="1">0.025</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">T(K)</td><td align="left" rowspan="1" colspan="1">2.79</td><td align="left" rowspan="1" colspan="1">2.15</td><td align="left" rowspan="1" colspan="1">N</td></tr><tr><td align="left" rowspan="1" colspan="1">198</td><td align="left" rowspan="1" colspan="1">3.23E-04</td><td align="left" rowspan="1" colspan="1">2.62E-03</td><td align="left" rowspan="1" colspan="1">T((V,A))</td><td align="left" rowspan="1" colspan="1">E((G,V))</td><td align="left" rowspan="1" colspan="1">2.28</td><td align="left" rowspan="1" colspan="1">2.47</td><td align="left" rowspan="1" colspan="1">A</td></tr><tr><td align="left" rowspan="1" colspan="1">214</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">T((N))</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">2.49</td><td align="left" rowspan="1" colspan="1">2.70</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">222</td><td align="left" rowspan="1" colspan="1">2.46E-03</td><td align="left" rowspan="1" colspan="1">0.013</td><td align="left" rowspan="1" colspan="1">K((R))</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">2.62</td><td align="left" rowspan="1" colspan="1">2.89</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">225</td><td align="left" rowspan="1" colspan="1">6.46E-05</td><td align="left" rowspan="1" colspan="1">9.15E-04</td><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">D((G,N))</td><td align="left" rowspan="1" colspan="1">2.60</td><td align="left" rowspan="1" colspan="1">2.83</td><td align="left" rowspan="1" colspan="1">D</td></tr><tr><td align="left" rowspan="1" colspan="1">238</td><td align="left" rowspan="1" colspan="1">1.82E-05</td><td align="left" rowspan="1" colspan="1">7.42E-04</td><td align="left" rowspan="1" colspan="1">D</td><td align="left" rowspan="1" colspan="1">E</td><td align="left" rowspan="1" colspan="1">3.07</td><td align="left" rowspan="1" colspan="1">2.70</td><td align="left" rowspan="1" colspan="1">E</td></tr><tr><td align="left" rowspan="1" colspan="1">239</td><td align="left" rowspan="1" colspan="1">4.90E-05</td><td align="left" rowspan="1" colspan="1">7.86E-04</td><td align="left" rowspan="1" colspan="1">Q</td><td align="left" rowspan="1" colspan="1">P</td><td align="left" rowspan="1" colspan="1">3.34</td><td align="left" rowspan="1" colspan="1">3.21</td><td align="left" rowspan="1" colspan="1">P</td></tr><tr><td align="left" rowspan="1" colspan="1">244</td><td align="left" rowspan="1" colspan="1">1.50E-03</td><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">I((M))</td><td align="left" rowspan="1" colspan="1">2.73</td><td align="left" rowspan="1" colspan="1">2.56</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td align="left" rowspan="1" colspan="1">248</td><td align="left" rowspan="1" colspan="1">2.08E-03</td><td align="left" rowspan="1" colspan="1">0.012</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">N((S))</td><td align="left" rowspan="1" colspan="1">2.73</td><td align="left" rowspan="1" colspan="1">2.67</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td align="left" rowspan="1" colspan="1">261</td><td align="left" rowspan="1" colspan="1">1.21E-03</td><td align="left" rowspan="1" colspan="1">8.54E-03</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">S((N))</td><td align="left" rowspan="1" colspan="1">2.79</td><td align="left" rowspan="1" colspan="1">2.45</td><td align="left" rowspan="1" colspan="1">E</td></tr><tr><td align="left" rowspan="1" colspan="1">262</td><td align="left" rowspan="1" colspan="1">1.46E-04</td><td align="left" rowspan="1" colspan="1">1.48E-03</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">R</td><td align="left" rowspan="1" colspan="1">2.89</td><td align="left" rowspan="1" colspan="1">2.82</td><td align="left" rowspan="1" colspan="1">R</td></tr><tr><td align="left" rowspan="1" colspan="1">271A</td><td align="left" rowspan="1" colspan="1">1.55E-04</td><td align="left" rowspan="1" colspan="1">1.48E-03</td><td align="left" rowspan="1" colspan="1">D</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">3.07</td><td align="left" rowspan="1" colspan="1">2.79</td><td align="left" rowspan="1" colspan="1">D</td></tr><tr><td align="left" rowspan="1" colspan="1">272</td><td align="left" rowspan="1" colspan="1">2.98E-03</td><td align="left" rowspan="1" colspan="1">0.015</td><td align="left" rowspan="1" colspan="1">A(T,V)</td><td align="left" rowspan="1" colspan="1">A</td><td align="left" rowspan="1" colspan="1">1.97</td><td align="left" rowspan="1" colspan="1">2.87</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td align="left" rowspan="1" colspan="1">274</td><td align="left" rowspan="1" colspan="1">2.18E-05</td><td align="left" rowspan="1" colspan="1">7.42E-04</td><td align="left" rowspan="1" colspan="1">V((I))</td><td align="left" rowspan="1" colspan="1">M</td><td align="left" rowspan="1" colspan="1">2.73</td><td align="left" rowspan="1" colspan="1">3.42</td><td align="left" rowspan="1" colspan="1">V</td></tr><tr><td align="left" rowspan="1" colspan="1">279</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">A((S))</td><td align="left" rowspan="1" colspan="1">2.73</td><td align="left" rowspan="1" colspan="1">2.77</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td align="left" rowspan="1" colspan="1">280</td><td align="left" rowspan="1" colspan="1">1.90E-03</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">R(K)</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">2.18</td><td align="left" rowspan="1" colspan="1">2.89</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td align="left" rowspan="1" colspan="1">285</td><td align="left" rowspan="1" colspan="1">1.48E-05</td><td align="left" rowspan="1" colspan="1">7.42E-04</td><td align="left" rowspan="1" colspan="1">H((Y,R))</td><td align="left" rowspan="1" colspan="1">Q</td><td align="left" rowspan="1" colspan="1">3.60</td><td align="left" rowspan="1" colspan="1">3.20</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">288</td><td align="left" rowspan="1" colspan="1">4.90E-05</td><td align="left" rowspan="1" colspan="1">7.86E-04</td><td align="left" rowspan="1" colspan="1">L</td><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">2.63</td><td align="left" rowspan="1" colspan="1">2.66</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">300</td><td align="left" rowspan="1" colspan="1">7.61E-05</td><td align="left" rowspan="1" colspan="1">9.95E-04</td><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">V</td><td align="left" rowspan="1" colspan="1">2.66</td><td align="left" rowspan="1" colspan="1">2.79</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">309</td><td align="left" rowspan="1" colspan="1">1.76E-03</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">V(I)</td><td align="left" rowspan="1" colspan="1">V</td><td align="left" rowspan="1" colspan="1">2.49</td><td align="left" rowspan="1" colspan="1">2.80</td><td align="left" rowspan="1" colspan="1">V</td></tr><tr><td align="left" rowspan="1" colspan="1">310</td><td align="left" rowspan="1" colspan="1">1.26E-04</td><td align="left" rowspan="1" colspan="1">1.48E-03</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">R</td><td align="left" rowspan="1" colspan="1">2.89</td><td align="left" rowspan="1" colspan="1">2.80</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">323</td><td align="left" rowspan="1" colspan="1">8.95E-03</td><td align="left" rowspan="1" colspan="1">0.037</td><td align="left" rowspan="1" colspan="1">V</td><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">2.83</td><td align="left" rowspan="1" colspan="1">2.61</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">HA2 72</td><td align="left" rowspan="1" colspan="1">3.32E-03</td><td align="left" rowspan="1" colspan="1">0.033</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">2.79</td><td align="left" rowspan="1" colspan="1">2.89</td><td align="left" rowspan="1" colspan="1">H</td></tr><tr><td align="left" rowspan="1" colspan="1">HA2 77</td><td align="left" rowspan="1" colspan="1">2.08E-04</td><td align="left" rowspan="1" colspan="1">4.17E-03</td><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">M</td><td align="left" rowspan="1" colspan="1">2.66</td><td align="left" rowspan="1" colspan="1">3.47</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">HA2 116</td><td align="left" rowspan="1" colspan="1">2.07E-04</td><td align="left" rowspan="1" colspan="1">4.17E-03</td><td align="left" rowspan="1" colspan="1">R</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">2.83</td><td align="left" rowspan="1" colspan="1">2.89</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">HA2 127</td><td align="left" rowspan="1" colspan="1">7.99E-04</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">R</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">2.83</td><td align="left" rowspan="1" colspan="1">2.89</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td colspan="8" align="left" rowspan="1"><bold>H2</bold></td></tr><tr><td align="left" rowspan="1" colspan="1">186</td><td align="left" rowspan="1" colspan="1">5.32E-04</td><td align="left" rowspan="1" colspan="1">0.018</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">I(N)((K))</td><td align="left" rowspan="1" colspan="1">2.78</td><td align="left" rowspan="1" colspan="1">1.95</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">197</td><td align="left" rowspan="1" colspan="1">3.62E-04</td><td align="left" rowspan="1" colspan="1">0.018</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">E(K,N)</td><td align="left" rowspan="1" colspan="1">2.78</td><td align="left" rowspan="1" colspan="1">2.01</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">205</td><td align="left" rowspan="1" colspan="1">1.11E-04</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">S(V)</td><td align="left" rowspan="1" colspan="1">2.60</td><td align="left" rowspan="1" colspan="1">1.96</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">HA2 45</td><td align="left" rowspan="1" colspan="1">5.36E-03</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">I(V)</td><td align="left" rowspan="1" colspan="1">F</td><td align="left" rowspan="1" colspan="1">2.28</td><td align="left" rowspan="1" colspan="1">3.30</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">HA2 130</td><td align="left" rowspan="1" colspan="1">5.38E-03</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">A</td><td align="left" rowspan="1" colspan="1">V((A))</td><td align="left" rowspan="1" colspan="1">2.91</td><td align="left" rowspan="1" colspan="1">2.57</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">HA2 169</td><td align="left" rowspan="1" colspan="1">5.47E-03</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">2.79</td><td align="left" rowspan="1" colspan="1">2.89</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">HA2 180</td><td align="left" rowspan="1" colspan="1">2.11E-03</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">N((S))</td><td align="left" rowspan="1" colspan="1">S</td><td align="left" rowspan="1" colspan="1">2.63</td><td align="left" rowspan="1" colspan="1">2.50</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td colspan="8" align="left" rowspan="1"><bold>H3</bold></td></tr><tr><td align="left" rowspan="1" colspan="1">−7</td><td align="left" rowspan="1" colspan="1">2.46E-03</td><td align="left" rowspan="1" colspan="1">0.037</td><td align="left" rowspan="1" colspan="1">C(Y)</td><td align="left" rowspan="1" colspan="1">Y</td><td align="left" rowspan="1" colspan="1">2.92</td><td align="left" rowspan="1" colspan="1">3.43</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">5.04E-05</td><td align="left" rowspan="1" colspan="1">5.73E-03</td><td align="left" rowspan="1" colspan="1">(G,S,C)</td><td align="left" rowspan="1" colspan="1">A(T)((S))</td><td align="left" rowspan="1" colspan="1">1.79</td><td align="left" rowspan="1" colspan="1">2.13</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">2.13E-03</td><td align="left" rowspan="1" colspan="1">0.035</td><td align="left" rowspan="1" colspan="1">S((P))</td><td align="left" rowspan="1" colspan="1">P(S)</td><td align="left" rowspan="1" colspan="1">2.30</td><td align="left" rowspan="1" colspan="1">2.70</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">57</td><td align="left" rowspan="1" colspan="1">1.15E-03</td><td align="left" rowspan="1" colspan="1">0.028</td><td align="left" rowspan="1" colspan="1">K((R))</td><td align="left" rowspan="1" colspan="1">Q(R)</td><td align="left" rowspan="1" colspan="1">2.60</td><td align="left" rowspan="1" colspan="1">2.72</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">63</td><td align="left" rowspan="1" colspan="1">7.32E-04</td><td align="left" rowspan="1" colspan="1">0.024</td><td align="left" rowspan="1" colspan="1">D</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">3.07</td><td align="left" rowspan="1" colspan="1">2.79</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">67</td><td align="left" rowspan="1" colspan="1">1.34E-03</td><td align="left" rowspan="1" colspan="1">0.028</td><td align="left" rowspan="1" colspan="1">(I,V,M)</td><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">1.86</td><td align="left" rowspan="1" colspan="1">2.60</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">92</td><td align="left" rowspan="1" colspan="1">6.95E-05</td><td align="left" rowspan="1" colspan="1">5.73E-03</td><td align="left" rowspan="1" colspan="1">N((S))</td><td align="left" rowspan="1" colspan="1">K((T,N))</td><td align="left" rowspan="1" colspan="1">2.49</td><td align="left" rowspan="1" colspan="1">2.48</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">145</td><td align="left" rowspan="1" colspan="1">1.37E-03</td><td align="left" rowspan="1" colspan="1">0.028</td><td align="left" rowspan="1" colspan="1">N(R,S)</td><td align="left" rowspan="1" colspan="1">K(N)((S))</td><td align="left" rowspan="1" colspan="1">1.98</td><td align="left" rowspan="1" colspan="1">2.00</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">213</td><td align="left" rowspan="1" colspan="1">1.17E-04</td><td align="left" rowspan="1" colspan="1">6.46E-03</td><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">V</td><td align="left" rowspan="1" colspan="1">2.66</td><td align="left" rowspan="1" colspan="1">2.83</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">228</td><td align="left" rowspan="1" colspan="1">5.59E-04</td><td align="left" rowspan="1" colspan="1">0.023</td><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">S</td><td align="left" rowspan="1" colspan="1">2.60</td><td align="left" rowspan="1" colspan="1">2.50</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">244</td><td align="left" rowspan="1" colspan="1">1.77E-03</td><td align="left" rowspan="1" colspan="1">0.033</td><td align="left" rowspan="1" colspan="1">V</td><td align="left" rowspan="1" colspan="1">L((I))</td><td align="left" rowspan="1" colspan="1">2.83</td><td align="left" rowspan="1" colspan="1">2.42</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td colspan="8" align="left" rowspan="1"><bold>N1</bold></td></tr><tr><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">5.96E-03</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">P</td><td align="left" rowspan="1" colspan="1">P((T,S))</td><td align="left" rowspan="1" colspan="1">3.21</td><td align="left" rowspan="1" colspan="1">3.04</td><td align="left" rowspan="1" colspan="1">P</td></tr><tr><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">7.25E-03</td><td align="left" rowspan="1" colspan="1">0.043</td><td align="left" rowspan="1" colspan="1">M(I)</td><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">2.81</td><td align="left" rowspan="1" colspan="1">2.64</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">34</td><td align="left" rowspan="1" colspan="1">4.03E-03</td><td align="left" rowspan="1" colspan="1">0.031</td><td align="left" rowspan="1" colspan="1">V((G,I,A))</td><td align="left" rowspan="1" colspan="1">(I,V,A)</td><td align="left" rowspan="1" colspan="1">2.36</td><td align="left" rowspan="1" colspan="1">1.71</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">42</td><td align="left" rowspan="1" colspan="1">3.43E-03</td><td align="left" rowspan="1" colspan="1">0.028</td><td align="left" rowspan="1" colspan="1">(G,N)((S,D))</td><td align="left" rowspan="1" colspan="1">S((N))</td><td align="left" rowspan="1" colspan="1">1.63</td><td align="left" rowspan="1" colspan="1">2.35</td><td align="left" rowspan="1" colspan="1">N</td></tr><tr><td align="left" rowspan="1" colspan="1">46</td><td align="left" rowspan="1" colspan="1">1.96E-03</td><td align="left" rowspan="1" colspan="1">0.017</td><td align="left" rowspan="1" colspan="1">(A,P,V,T,S)</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">1.35</td><td align="left" rowspan="1" colspan="1">2.72</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">47</td><td align="left" rowspan="1" colspan="1">4.34E-03</td><td align="left" rowspan="1" colspan="1">0.031</td><td align="left" rowspan="1" colspan="1">E(G)((D))</td><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">1.92</td><td align="left" rowspan="1" colspan="1">2.59</td><td align="left" rowspan="1" colspan="1">E</td></tr><tr><td align="left" rowspan="1" colspan="1">52</td><td align="left" rowspan="1" colspan="1">1.42E-03</td><td align="left" rowspan="1" colspan="1">0.013</td><td align="left" rowspan="1" colspan="1">S</td><td align="left" rowspan="1" colspan="1">R((G,N,K))</td><td align="left" rowspan="1" colspan="1">2.50</td><td align="left" rowspan="1" colspan="1">2.51</td><td align="left" rowspan="1" colspan="1">S</td></tr><tr><td align="left" rowspan="1" colspan="1">59</td><td align="left" rowspan="1" colspan="1">2.08E-03</td><td align="left" rowspan="1" colspan="1">0.018</td><td align="left" rowspan="1" colspan="1">N((K))</td><td align="left" rowspan="1" colspan="1">S((N,R))</td><td align="left" rowspan="1" colspan="1">2.65</td><td align="left" rowspan="1" colspan="1">2.19</td><td align="left" rowspan="1" colspan="1">N</td></tr><tr><td align="left" rowspan="1" colspan="1">67</td><td align="left" rowspan="1" colspan="1">6.62E-03</td><td align="left" rowspan="1" colspan="1">0.043</td><td align="left" rowspan="1" colspan="1">(L,I,V)</td><td align="left" rowspan="1" colspan="1">V</td><td align="left" rowspan="1" colspan="1">1.63</td><td align="left" rowspan="1" colspan="1">2.83</td><td align="left" rowspan="1" colspan="1">V</td></tr><tr><td align="left" rowspan="1" colspan="1">74</td><td align="left" rowspan="1" colspan="1">3.23E-04</td><td align="left" rowspan="1" colspan="1">6.33E-03</td><td align="left" rowspan="1" colspan="1">(L,F,S,V)((I))</td><td align="left" rowspan="1" colspan="1">V</td><td align="left" rowspan="1" colspan="1">1.28</td><td align="left" rowspan="1" colspan="1">2.81</td><td align="left" rowspan="1" colspan="1">F</td></tr><tr><td align="left" rowspan="1" colspan="1">80</td><td align="left" rowspan="1" colspan="1">7.93E-04</td><td align="left" rowspan="1" colspan="1">9.32E-03</td><td align="left" rowspan="1" colspan="1">V((R,A,M))</td><td align="left" rowspan="1" colspan="1">(I,V,K)((T,S))</td><td align="left" rowspan="1" colspan="1">2.21</td><td align="left" rowspan="1" colspan="1">1.32</td><td align="left" rowspan="1" colspan="1">V</td></tr><tr><td align="left" rowspan="1" colspan="1">157</td><td align="left" rowspan="1" colspan="1">2.11E-04</td><td align="left" rowspan="1" colspan="1">5.47E-03</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">A</td><td align="left" rowspan="1" colspan="1">2.71</td><td align="left" rowspan="1" colspan="1">2.87</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td align="left" rowspan="1" colspan="1">189</td><td align="left" rowspan="1" colspan="1">3.08E-08</td><td align="left" rowspan="1" colspan="1">6.04E-06</td><td align="left" rowspan="1" colspan="1">S((G))</td><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">2.40</td><td align="left" rowspan="1" colspan="1">2.60</td><td align="left" rowspan="1" colspan="1">N</td></tr><tr><td align="left" rowspan="1" colspan="1">214</td><td align="left" rowspan="1" colspan="1">4.73E-06</td><td align="left" rowspan="1" colspan="1">3.09E-04</td><td align="left" rowspan="1" colspan="1">D</td><td align="left" rowspan="1" colspan="1">E(G)</td><td align="left" rowspan="1" colspan="1">3.07</td><td align="left" rowspan="1" colspan="1">2.15</td><td align="left" rowspan="1" colspan="1">D</td></tr><tr><td align="left" rowspan="1" colspan="1">220</td><td align="left" rowspan="1" colspan="1">4.22E-04</td><td align="left" rowspan="1" colspan="1">6.82E-03</td><td align="left" rowspan="1" colspan="1">R((G))</td><td align="left" rowspan="1" colspan="1">K(R)</td><td align="left" rowspan="1" colspan="1">2.57</td><td align="left" rowspan="1" colspan="1">2.37</td><td align="left" rowspan="1" colspan="1">R</td></tr><tr><td align="left" rowspan="1" colspan="1">221</td><td align="left" rowspan="1" colspan="1">4.52E-04</td><td align="left" rowspan="1" colspan="1">6.82E-03</td><td align="left" rowspan="1" colspan="1">N(G)</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">2.41</td><td align="left" rowspan="1" colspan="1">2.87</td><td align="left" rowspan="1" colspan="1">N</td></tr><tr><td align="left" rowspan="1" colspan="1">264</td><td align="left" rowspan="1" colspan="1">8.09E-04</td><td align="left" rowspan="1" colspan="1">9.32E-03</td><td align="left" rowspan="1" colspan="1">(I,A,V)</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">1.71</td><td align="left" rowspan="1" colspan="1">2.68</td><td align="left" rowspan="1" colspan="1">V</td></tr><tr><td align="left" rowspan="1" colspan="1">274</td><td align="left" rowspan="1" colspan="1">2.80E-05</td><td align="left" rowspan="1" colspan="1">9.13E-04</td><td align="left" rowspan="1" colspan="1">Y</td><td align="left" rowspan="1" colspan="1">(S,F)((Y))</td><td align="left" rowspan="1" colspan="1">3.50</td><td align="left" rowspan="1" colspan="1">2.12</td><td align="left" rowspan="1" colspan="1">Y</td></tr><tr><td align="left" rowspan="1" colspan="1">288</td><td align="left" rowspan="1" colspan="1">5.84E-04</td><td align="left" rowspan="1" colspan="1">8.18E-03</td><td align="left" rowspan="1" colspan="1">I(V)</td><td align="left" rowspan="1" colspan="1">V</td><td align="left" rowspan="1" colspan="1">2.08</td><td align="left" rowspan="1" colspan="1">2.83</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">289</td><td align="left" rowspan="1" colspan="1">8.64E-04</td><td align="left" rowspan="1" colspan="1">9.41E-03</td><td align="left" rowspan="1" colspan="1">(I,T,M)</td><td align="left" rowspan="1" colspan="1">M</td><td align="left" rowspan="1" colspan="1">1.82</td><td align="left" rowspan="1" colspan="1">3.47</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td align="left" rowspan="1" colspan="1">309</td><td align="left" rowspan="1" colspan="1">6.95E-03</td><td align="left" rowspan="1" colspan="1">0.043</td><td align="left" rowspan="1" colspan="1">N(D)</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">2.37</td><td align="left" rowspan="1" colspan="1">2.79</td><td align="left" rowspan="1" colspan="1">N</td></tr><tr><td align="left" rowspan="1" colspan="1">311</td><td align="left" rowspan="1" colspan="1">2.12E-05</td><td align="left" rowspan="1" colspan="1">8.30E-04</td><td align="left" rowspan="1" colspan="1">E((D))</td><td align="left" rowspan="1" colspan="1">D</td><td align="left" rowspan="1" colspan="1">2.59</td><td align="left" rowspan="1" colspan="1">3.07</td><td align="left" rowspan="1" colspan="1">E</td></tr><tr><td align="left" rowspan="1" colspan="1">329</td><td align="left" rowspan="1" colspan="1">7.22E-03</td><td align="left" rowspan="1" colspan="1">0.043</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">K(E)((R))</td><td align="left" rowspan="1" colspan="1">2.79</td><td align="left" rowspan="1" colspan="1">2.43</td><td align="left" rowspan="1" colspan="1">N</td></tr><tr><td align="left" rowspan="1" colspan="1">339</td><td align="left" rowspan="1" colspan="1">1.27E-03</td><td align="left" rowspan="1" colspan="1">0.012</td><td align="left" rowspan="1" colspan="1">S((L))</td><td align="left" rowspan="1" colspan="1">T(Y)((N))</td><td align="left" rowspan="1" colspan="1">2.40</td><td align="left" rowspan="1" colspan="1">2.11</td><td align="left" rowspan="1" colspan="1">S</td></tr><tr><td align="left" rowspan="1" colspan="1">340</td><td align="left" rowspan="1" colspan="1">1.18E-05</td><td align="left" rowspan="1" colspan="1">5.80E-04</td><td align="left" rowspan="1" colspan="1">(L,S,P)((H))</td><td align="left" rowspan="1" colspan="1">V((A,H,P))</td><td align="left" rowspan="1" colspan="1">1.55</td><td align="left" rowspan="1" colspan="1">2.26</td><td align="left" rowspan="1" colspan="1">S</td></tr><tr><td align="left" rowspan="1" colspan="1">341</td><td align="left" rowspan="1" colspan="1">4.13E-04</td><td align="left" rowspan="1" colspan="1">6.82E-03</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">D</td><td align="left" rowspan="1" colspan="1">2.76</td><td align="left" rowspan="1" colspan="1">3.05</td><td align="left" rowspan="1" colspan="1">N</td></tr><tr><td align="left" rowspan="1" colspan="1">351</td><td align="left" rowspan="1" colspan="1">6.82E-04</td><td align="left" rowspan="1" colspan="1">8.91E-03</td><td align="left" rowspan="1" colspan="1">F((Y))</td><td align="left" rowspan="1" colspan="1">Y</td><td align="left" rowspan="1" colspan="1">3.01</td><td align="left" rowspan="1" colspan="1">3.49</td><td align="left" rowspan="1" colspan="1">F</td></tr><tr><td align="left" rowspan="1" colspan="1">365</td><td align="left" rowspan="1" colspan="1">1.20E-03</td><td align="left" rowspan="1" colspan="1">0.012</td><td align="left" rowspan="1" colspan="1">T(I,P)</td><td align="left" rowspan="1" colspan="1">N((S,T))</td><td align="left" rowspan="1" colspan="1">2.05</td><td align="left" rowspan="1" colspan="1">2.59</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">382</td><td align="left" rowspan="1" colspan="1">8.60E-08</td><td align="left" rowspan="1" colspan="1">8.42E-06</td><td align="left" rowspan="1" colspan="1">E((G,D))</td><td align="left" rowspan="1" colspan="1">D(N)</td><td align="left" rowspan="1" colspan="1">2.35</td><td align="left" rowspan="1" colspan="1">2.30</td><td align="left" rowspan="1" colspan="1">G</td></tr><tr><td align="left" rowspan="1" colspan="1">393</td><td align="left" rowspan="1" colspan="1">4.35E-03</td><td align="left" rowspan="1" colspan="1">0.031</td><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">V(I)</td><td align="left" rowspan="1" colspan="1">2.63</td><td align="left" rowspan="1" colspan="1">2.30</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">427</td><td align="left" rowspan="1" colspan="1">4.44E-03</td><td align="left" rowspan="1" colspan="1">0.031</td><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">V(I)</td><td align="left" rowspan="1" colspan="1">2.66</td><td align="left" rowspan="1" colspan="1">2.46</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">430</td><td align="left" rowspan="1" colspan="1">2.23E-04</td><td align="left" rowspan="1" colspan="1">5.47E-03</td><td align="left" rowspan="1" colspan="1">R((L))</td><td align="left" rowspan="1" colspan="1">L((Q,R))</td><td align="left" rowspan="1" colspan="1">2.59</td><td align="left" rowspan="1" colspan="1">2.30</td><td align="left" rowspan="1" colspan="1">R</td></tr><tr><td align="left" rowspan="1" colspan="1">455</td><td align="left" rowspan="1" colspan="1">2.60E-04</td><td align="left" rowspan="1" colspan="1">5.66E-03</td><td align="left" rowspan="1" colspan="1">G(S,D)</td><td align="left" rowspan="1" colspan="1">N(D)</td><td align="left" rowspan="1" colspan="1">1.65</td><td align="left" rowspan="1" colspan="1">2.47</td><td align="left" rowspan="1" colspan="1">W</td></tr><tr><td colspan="8" align="left" rowspan="1"><bold>N2</bold></td></tr><tr><td align="left" rowspan="1" colspan="1">41</td><td align="left" rowspan="1" colspan="1">1.58E-03</td><td align="left" rowspan="1" colspan="1">0.025</td><td align="left" rowspan="1" colspan="1">E((G))</td><td align="left" rowspan="1" colspan="1">E</td><td align="left" rowspan="1" colspan="1">2.43</td><td align="left" rowspan="1" colspan="1">2.74</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">4.12E-03</td><td align="left" rowspan="1" colspan="1">0.039</td><td align="left" rowspan="1" colspan="1">V(A)((T,I))</td><td align="left" rowspan="1" colspan="1">V((A))</td><td align="left" rowspan="1" colspan="1">1.98</td><td align="left" rowspan="1" colspan="1">2.69</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">51</td><td align="left" rowspan="1" colspan="1">5.36E-04</td><td align="left" rowspan="1" colspan="1">0.012</td><td align="left" rowspan="1" colspan="1">V((M,T))</td><td align="left" rowspan="1" colspan="1">M</td><td align="left" rowspan="1" colspan="1">2.39</td><td align="left" rowspan="1" colspan="1">3.43</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">60</td><td align="left" rowspan="1" colspan="1">5.79E-03</td><td align="left" rowspan="1" colspan="1">0.047</td><td align="left" rowspan="1" colspan="1">R(K)</td><td align="left" rowspan="1" colspan="1">R</td><td align="left" rowspan="1" colspan="1">2.30</td><td align="left" rowspan="1" colspan="1">2.83</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">62</td><td align="left" rowspan="1" colspan="1">4.28E-03</td><td align="left" rowspan="1" colspan="1">0.039</td><td align="left" rowspan="1" colspan="1">(I,T,M)((V))</td><td align="left" rowspan="1" colspan="1">I((T))</td><td align="left" rowspan="1" colspan="1">1.59</td><td align="left" rowspan="1" colspan="1">2.59</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">70</td><td align="left" rowspan="1" colspan="1">1.82E-04</td><td align="left" rowspan="1" colspan="1">6.14E-03</td><td align="left" rowspan="1" colspan="1">S(H,N)</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">1.96</td><td align="left" rowspan="1" colspan="1">2.78</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">72</td><td align="left" rowspan="1" colspan="1">2.10E-03</td><td align="left" rowspan="1" colspan="1">0.027</td><td align="left" rowspan="1" colspan="1">T(I)((V))</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">2.00</td><td align="left" rowspan="1" colspan="1">2.70</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">77</td><td align="left" rowspan="1" colspan="1">2.05E-03</td><td align="left" rowspan="1" colspan="1">0.027</td><td align="left" rowspan="1" colspan="1">(I,K,T)((V,L))</td><td align="left" rowspan="1" colspan="1">I((K))</td><td align="left" rowspan="1" colspan="1">1.36</td><td align="left" rowspan="1" colspan="1">2.32</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">81</td><td align="left" rowspan="1" colspan="1">2.83E-04</td><td align="left" rowspan="1" colspan="1">8.22E-03</td><td align="left" rowspan="1" colspan="1">A((V,M,L,I,T))</td><td align="left" rowspan="1" colspan="1">(P,L,A,T)</td><td align="left" rowspan="1" colspan="1">1.92</td><td align="left" rowspan="1" colspan="1">1.62</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">83</td><td align="left" rowspan="1" colspan="1">7.97E-05</td><td align="left" rowspan="1" colspan="1">3.24E-03</td><td align="left" rowspan="1" colspan="1">G(E,D)((K))</td><td align="left" rowspan="1" colspan="1">E</td><td align="left" rowspan="1" colspan="1">1.61</td><td align="left" rowspan="1" colspan="1">2.70</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">125</td><td align="left" rowspan="1" colspan="1">7.62E-04</td><td align="left" rowspan="1" colspan="1">0.014</td><td align="left" rowspan="1" colspan="1">(G,S,D)</td><td align="left" rowspan="1" colspan="1">D</td><td align="left" rowspan="1" colspan="1">1.63</td><td align="left" rowspan="1" colspan="1">3.07</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">126</td><td align="left" rowspan="1" colspan="1">3.57E-03</td><td align="left" rowspan="1" colspan="1">0.038</td><td align="left" rowspan="1" colspan="1">(L,P,T)((H,S))</td><td align="left" rowspan="1" colspan="1">P((S))</td><td align="left" rowspan="1" colspan="1">1.68</td><td align="left" rowspan="1" colspan="1">3.08</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">147</td><td align="left" rowspan="1" colspan="1">7.60E-04</td><td align="left" rowspan="1" colspan="1">0.014</td><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">D((N))</td><td align="left" rowspan="1" colspan="1">2.60</td><td align="left" rowspan="1" colspan="1">2.94</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">155</td><td align="left" rowspan="1" colspan="1">4.04E-03</td><td align="left" rowspan="1" colspan="1">0.039</td><td align="left" rowspan="1" colspan="1">H</td><td align="left" rowspan="1" colspan="1">Y(H)</td><td align="left" rowspan="1" colspan="1">3.96</td><td align="left" rowspan="1" colspan="1">3.19</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">192</td><td align="left" rowspan="1" colspan="1">3.44E-04</td><td align="left" rowspan="1" colspan="1">8.74E-03</td><td align="left" rowspan="1" colspan="1">V(I)</td><td align="left" rowspan="1" colspan="1">V</td><td align="left" rowspan="1" colspan="1">2.48</td><td align="left" rowspan="1" colspan="1">2.83</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">216</td><td align="left" rowspan="1" colspan="1">6.64E-03</td><td align="left" rowspan="1" colspan="1">0.048</td><td align="left" rowspan="1" colspan="1">(G,V,A)((S))</td><td align="left" rowspan="1" colspan="1">V(G,S)</td><td align="left" rowspan="1" colspan="1">1.45</td><td align="left" rowspan="1" colspan="1">1.92</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">283</td><td align="left" rowspan="1" colspan="1">8.85E-04</td><td align="left" rowspan="1" colspan="1">0.015</td><td align="left" rowspan="1" colspan="1">R(Q)</td><td align="left" rowspan="1" colspan="1">R</td><td align="left" rowspan="1" colspan="1">2.36</td><td align="left" rowspan="1" colspan="1">2.83</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">286</td><td align="left" rowspan="1" colspan="1">2.03E-03</td><td align="left" rowspan="1" colspan="1">0.027</td><td align="left" rowspan="1" colspan="1">(I,E,N)((D))</td><td align="left" rowspan="1" colspan="1">G((D))</td><td align="left" rowspan="1" colspan="1">1.51</td><td align="left" rowspan="1" colspan="1">2.45</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">315</td><td align="left" rowspan="1" colspan="1">6.05E-03</td><td align="left" rowspan="1" colspan="1">0.047</td><td align="left" rowspan="1" colspan="1">G(S,R)((N))</td><td align="left" rowspan="1" colspan="1">S(R)</td><td align="left" rowspan="1" colspan="1">1.53</td><td align="left" rowspan="1" colspan="1">2.14</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">328</td><td align="left" rowspan="1" colspan="1">4.65E-05</td><td align="left" rowspan="1" colspan="1">2.36E-03</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">K((R))</td><td align="left" rowspan="1" colspan="1">2.79</td><td align="left" rowspan="1" colspan="1">2.70</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">331</td><td align="left" rowspan="1" colspan="1">2.49E-03</td><td align="left" rowspan="1" colspan="1">0.028</td><td align="left" rowspan="1" colspan="1">(I,R,G,S)</td><td align="left" rowspan="1" colspan="1">S(R)((N))</td><td align="left" rowspan="1" colspan="1">1.38</td><td align="left" rowspan="1" colspan="1">1.86</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">338</td><td align="left" rowspan="1" colspan="1">4.47E-03</td><td align="left" rowspan="1" colspan="1">0.039</td><td align="left" rowspan="1" colspan="1">R(K)</td><td align="left" rowspan="1" colspan="1">L(Q,W)((K,R))</td><td align="left" rowspan="1" colspan="1">2.34</td><td align="left" rowspan="1" colspan="1">1.79</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">369</td><td align="left" rowspan="1" colspan="1">5.46E-03</td><td align="left" rowspan="1" colspan="1">0.046</td><td align="left" rowspan="1" colspan="1">D</td><td align="left" rowspan="1" colspan="1">K(E)</td><td align="left" rowspan="1" colspan="1">3.07</td><td align="left" rowspan="1" colspan="1">2.18</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">378</td><td align="left" rowspan="1" colspan="1">2.43E-03</td><td align="left" rowspan="1" colspan="1">0.028</td><td align="left" rowspan="1" colspan="1">R(K)</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">2.29</td><td align="left" rowspan="1" colspan="1">2.89</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">381</td><td align="left" rowspan="1" colspan="1">1.46E-05</td><td align="left" rowspan="1" colspan="1">9.88E-04</td><td align="left" rowspan="1" colspan="1">G((D,N))</td><td align="left" rowspan="1" colspan="1">E(D)</td><td align="left" rowspan="1" colspan="1">2.36</td><td align="left" rowspan="1" colspan="1">2.42</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">384</td><td align="left" rowspan="1" colspan="1">4.38E-06</td><td align="left" rowspan="1" colspan="1">4.45E-04</td><td align="left" rowspan="1" colspan="1">(A,T,I)((V,N,S))</td><td align="left" rowspan="1" colspan="1">V(I)</td><td align="left" rowspan="1" colspan="1">1.18</td><td align="left" rowspan="1" colspan="1">2.06</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">386</td><td align="left" rowspan="1" colspan="1">2.74E-06</td><td align="left" rowspan="1" colspan="1">4.45E-04</td><td align="left" rowspan="1" colspan="1">A((P))</td><td align="left" rowspan="1" colspan="1">P((S))</td><td align="left" rowspan="1" colspan="1">2.68</td><td align="left" rowspan="1" colspan="1">3.09</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">396</td><td align="left" rowspan="1" colspan="1">6.80E-03</td><td align="left" rowspan="1" colspan="1">0.048</td><td align="left" rowspan="1" colspan="1">V(I)</td><td align="left" rowspan="1" colspan="1">V</td><td align="left" rowspan="1" colspan="1">2.21</td><td align="left" rowspan="1" colspan="1">2.83</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">399</td><td align="left" rowspan="1" colspan="1">6.69E-03</td><td align="left" rowspan="1" colspan="1">0.048</td><td align="left" rowspan="1" colspan="1">D</td><td align="left" rowspan="1" colspan="1">E</td><td align="left" rowspan="1" colspan="1">3.06</td><td align="left" rowspan="1" colspan="1">2.74</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td colspan="8" align="left" rowspan="1"><bold>M1</bold></td></tr><tr><td align="left" rowspan="1" colspan="1">115</td><td align="left" rowspan="1" colspan="1">2.05E-04</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">V</td><td align="left" rowspan="1" colspan="1">I((V))</td><td align="left" rowspan="1" colspan="1">2.79</td><td align="left" rowspan="1" colspan="1">2.55</td><td align="left" rowspan="1" colspan="1">V</td></tr><tr><td align="left" rowspan="1" colspan="1">137</td><td align="left" rowspan="1" colspan="1">1.30E-04</td><td align="left" rowspan="1" colspan="1">0.011</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">A((T))</td><td align="left" rowspan="1" colspan="1">2.71</td><td align="left" rowspan="1" colspan="1">2.76</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td align="left" rowspan="1" colspan="1">174</td><td align="left" rowspan="1" colspan="1">1.06E-03</td><td align="left" rowspan="1" colspan="1">0.029</td><td align="left" rowspan="1" colspan="1">R</td><td align="left" rowspan="1" colspan="1">K(R)</td><td align="left" rowspan="1" colspan="1">2.80</td><td align="left" rowspan="1" colspan="1">2.20</td><td align="left" rowspan="1" colspan="1">R</td></tr><tr><td align="left" rowspan="1" colspan="1">231</td><td align="left" rowspan="1" colspan="1">4.64E-04</td><td align="left" rowspan="1" colspan="1">0.017</td><td align="left" rowspan="1" colspan="1">D</td><td align="left" rowspan="1" colspan="1">(N,D,S)</td><td align="left" rowspan="1" colspan="1">3.02</td><td align="left" rowspan="1" colspan="1">1.70</td><td align="left" rowspan="1" colspan="1">D</td></tr><tr><td colspan="8" align="left" rowspan="1"><bold>M2</bold></td></tr><tr><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">9.63E-04</td><td align="left" rowspan="1" colspan="1">0.039</td><td align="left" rowspan="1" colspan="1">(L,H,P)</td><td align="left" rowspan="1" colspan="1">P</td><td align="left" rowspan="1" colspan="1">2.20</td><td align="left" rowspan="1" colspan="1">3.21</td><td align="left" rowspan="1" colspan="1">P</td></tr><tr><td align="left" rowspan="1" colspan="1">93</td><td align="left" rowspan="1" colspan="1">4.26E-04</td><td align="left" rowspan="1" colspan="1">0.035</td><td align="left" rowspan="1" colspan="1">N((Y,I,S))</td><td align="left" rowspan="1" colspan="1">(S,I,Q)((N))</td><td align="left" rowspan="1" colspan="1">2.39</td><td align="left" rowspan="1" colspan="1">1.49</td><td align="left" rowspan="1" colspan="1">N</td></tr><tr><td colspan="8" align="left" rowspan="1"><bold>NS1</bold></td></tr><tr><td align="left" rowspan="1" colspan="1">81</td><td align="left" rowspan="1" colspan="1">5.08E-04</td><td align="left" rowspan="1" colspan="1">0.029</td><td align="left" rowspan="1" colspan="1">I(T)((V,M))</td><td align="left" rowspan="1" colspan="1">M((V))</td><td align="left" rowspan="1" colspan="1">1.93</td><td align="left" rowspan="1" colspan="1">3.39</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">84</td><td align="left" rowspan="1" colspan="1">4.84E-04</td><td align="left" rowspan="1" colspan="1">0.029</td><td align="left" rowspan="1" colspan="1">(V,M,G,S)((L,A,I,T))</td><td align="left" rowspan="1" colspan="1">T(A)((V))</td><td align="left" rowspan="1" colspan="1">1.17</td><td align="left" rowspan="1" colspan="1">1.97</td><td align="left" rowspan="1" colspan="1">V</td></tr><tr><td align="left" rowspan="1" colspan="1">215</td><td align="left" rowspan="1" colspan="1">8.84E-06</td><td align="left" rowspan="1" colspan="1">1.55E-03</td><td align="left" rowspan="1" colspan="1">(P,S,L)((T,A))</td><td align="left" rowspan="1" colspan="1">T((P))</td><td align="left" rowspan="1" colspan="1">1.65</td><td align="left" rowspan="1" colspan="1">2.60</td><td align="left" rowspan="1" colspan="1">P</td></tr><tr><td align="left" rowspan="1" colspan="1">227</td><td align="left" rowspan="1" colspan="1">6.65E-04</td><td align="left" rowspan="1" colspan="1">0.029</td><td align="left" rowspan="1" colspan="1">E((G,K))</td><td align="left" rowspan="1" colspan="1">R(G)((E))</td><td align="left" rowspan="1" colspan="1">2.43</td><td align="left" rowspan="1" colspan="1">2.19</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td colspan="8" align="left" rowspan="1"><bold>NP</bold></td></tr><tr><td align="left" rowspan="1" colspan="1">77</td><td align="left" rowspan="1" colspan="1">1.37E-04</td><td align="left" rowspan="1" colspan="1">8.33E-03</td><td align="left" rowspan="1" colspan="1">R(K)</td><td align="left" rowspan="1" colspan="1">K(R)</td><td align="left" rowspan="1" colspan="1">2.43</td><td align="left" rowspan="1" colspan="1">2.52</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">101</td><td align="left" rowspan="1" colspan="1">1.88E-07</td><td align="left" rowspan="1" colspan="1">3.45E-05</td><td align="left" rowspan="1" colspan="1">(E,D,N)</td><td align="left" rowspan="1" colspan="1">G(N)((D))</td><td align="left" rowspan="1" colspan="1">1.90</td><td align="left" rowspan="1" colspan="1">1.96</td><td align="left" rowspan="1" colspan="1">D</td></tr><tr><td align="left" rowspan="1" colspan="1">102</td><td align="left" rowspan="1" colspan="1">3.60E-04</td><td align="left" rowspan="1" colspan="1">0.013</td><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">G((R))</td><td align="left" rowspan="1" colspan="1">2.60</td><td align="left" rowspan="1" colspan="1">2.44</td><td align="left" rowspan="1" colspan="1">G</td></tr><tr><td align="left" rowspan="1" colspan="1">131</td><td align="left" rowspan="1" colspan="1">2.71E-04</td><td align="left" rowspan="1" colspan="1">0.012</td><td align="left" rowspan="1" colspan="1">A</td><td align="left" rowspan="1" colspan="1">A(R)</td><td align="left" rowspan="1" colspan="1">2.91</td><td align="left" rowspan="1" colspan="1">2.51</td><td align="left" rowspan="1" colspan="1">A</td></tr><tr><td align="left" rowspan="1" colspan="1">136</td><td align="left" rowspan="1" colspan="1">6.30E-04</td><td align="left" rowspan="1" colspan="1">0.019</td><td align="left" rowspan="1" colspan="1">(L,M,I)</td><td align="left" rowspan="1" colspan="1">I(M)</td><td align="left" rowspan="1" colspan="1">1.91</td><td align="left" rowspan="1" colspan="1">2.32</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">283</td><td align="left" rowspan="1" colspan="1">3.47E-03</td><td align="left" rowspan="1" colspan="1">0.049</td><td align="left" rowspan="1" colspan="1">L</td><td align="left" rowspan="1" colspan="1">P</td><td align="left" rowspan="1" colspan="1">2.63</td><td align="left" rowspan="1" colspan="1">3.21</td><td align="left" rowspan="1" colspan="1">L</td></tr><tr><td align="left" rowspan="1" colspan="1">305</td><td align="left" rowspan="1" colspan="1">1.42E-03</td><td align="left" rowspan="1" colspan="1">0.032</td><td align="left" rowspan="1" colspan="1">R(K)</td><td align="left" rowspan="1" colspan="1">K((R))</td><td align="left" rowspan="1" colspan="1">2.43</td><td align="left" rowspan="1" colspan="1">2.80</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">335</td><td align="left" rowspan="1" colspan="1">9.62E-04</td><td align="left" rowspan="1" colspan="1">0.025</td><td align="left" rowspan="1" colspan="1">S</td><td align="left" rowspan="1" colspan="1">S((F))</td><td align="left" rowspan="1" colspan="1">2.50</td><td align="left" rowspan="1" colspan="1">2.38</td><td align="left" rowspan="1" colspan="1">S</td></tr><tr><td align="left" rowspan="1" colspan="1">353</td><td align="left" rowspan="1" colspan="1">3.10E-03</td><td align="left" rowspan="1" colspan="1">0.047</td><td align="left" rowspan="1" colspan="1">(V,I,L)((A))</td><td align="left" rowspan="1" colspan="1">(C,S,F,L)((I,V))</td><td align="left" rowspan="1" colspan="1">1.45</td><td align="left" rowspan="1" colspan="1">1.48</td><td align="left" rowspan="1" colspan="1">I</td></tr><tr><td align="left" rowspan="1" colspan="1">357</td><td align="left" rowspan="1" colspan="1">1.96E-03</td><td align="left" rowspan="1" colspan="1">0.036</td><td align="left" rowspan="1" colspan="1">Q</td><td align="left" rowspan="1" colspan="1">K((R))</td><td align="left" rowspan="1" colspan="1">3.29</td><td align="left" rowspan="1" colspan="1">2.65</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">375</td><td align="left" rowspan="1" colspan="1">1.26E-04</td><td align="left" rowspan="1" colspan="1">8.33E-03</td><td align="left" rowspan="1" colspan="1">(V,D,E,G)((S,N))</td><td align="left" rowspan="1" colspan="1">G(V)((E))</td><td align="left" rowspan="1" colspan="1">1.30</td><td align="left" rowspan="1" colspan="1">1.97</td><td align="left" rowspan="1" colspan="1">D</td></tr><tr><td align="left" rowspan="1" colspan="1">425</td><td align="left" rowspan="1" colspan="1">1.88E-03</td><td align="left" rowspan="1" colspan="1">0.036</td><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">I(V)</td><td align="left" rowspan="1" colspan="1">2.66</td><td align="left" rowspan="1" colspan="1">2.11</td><td align="left" rowspan="1" colspan="1">V</td></tr><tr><td align="left" rowspan="1" colspan="1">472</td><td align="left" rowspan="1" colspan="1">2.54E-03</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">T(A)</td><td align="left" rowspan="1" colspan="1">2.73</td><td align="left" rowspan="1" colspan="1">2.31</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td colspan="8" align="left" rowspan="1"><bold>PA</bold></td></tr><tr><td align="left" rowspan="1" colspan="1">356</td><td align="left" rowspan="1" colspan="1">2.66E-04</td><td align="left" rowspan="1" colspan="1">0.035</td><td align="left" rowspan="1" colspan="1">K((R))</td><td align="left" rowspan="1" colspan="1">R((K))</td><td align="left" rowspan="1" colspan="1">2.70</td><td align="left" rowspan="1" colspan="1">2.68</td><td align="left" rowspan="1" colspan="1">R</td></tr><tr><td align="left" rowspan="1" colspan="1">552</td><td align="left" rowspan="1" colspan="1">1.94E-04</td><td align="left" rowspan="1" colspan="1">0.035</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">S</td><td align="left" rowspan="1" colspan="1">2.73</td><td align="left" rowspan="1" colspan="1">2.50</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td colspan="8" align="left" rowspan="1"><bold>PB1</bold></td></tr><tr><td align="left" rowspan="1" colspan="1">52</td><td align="left" rowspan="1" colspan="1">4.10E-04</td><td align="left" rowspan="1" colspan="1">0.032</td><td align="left" rowspan="1" colspan="1">K((R))</td><td align="left" rowspan="1" colspan="1">R(K)</td><td align="left" rowspan="1" colspan="1">2.74</td><td align="left" rowspan="1" colspan="1">2.17</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">517</td><td align="left" rowspan="1" colspan="1">3.63E-04</td><td align="left" rowspan="1" colspan="1">0.032</td><td align="left" rowspan="1" colspan="1">I((V))</td><td align="left" rowspan="1" colspan="1">V(I)</td><td align="left" rowspan="1" colspan="1">2.58</td><td align="left" rowspan="1" colspan="1">2.06</td><td align="left" rowspan="1" colspan="1">V</td></tr><tr><td align="left" rowspan="1" colspan="1">584</td><td align="left" rowspan="1" colspan="1">7.67E-07</td><td align="left" rowspan="1" colspan="1">1.81E-04</td><td align="left" rowspan="1" colspan="1">R((H))</td><td align="left" rowspan="1" colspan="1">Q(H)</td><td align="left" rowspan="1" colspan="1">2.73</td><td align="left" rowspan="1" colspan="1">2.93</td><td align="left" rowspan="1" colspan="1">Q</td></tr><tr><td colspan="8" align="left" rowspan="1"><bold>PB2</bold></td></tr><tr><td align="left" rowspan="1" colspan="1">44</td><td align="left" rowspan="1" colspan="1">6.18E-04</td><td align="left" rowspan="1" colspan="1">0.023</td><td align="left" rowspan="1" colspan="1">A(S)</td><td align="left" rowspan="1" colspan="1">L(S)</td><td align="left" rowspan="1" colspan="1">2.53</td><td align="left" rowspan="1" colspan="1">2.15</td><td align="left" rowspan="1" colspan="1">A</td></tr><tr><td align="left" rowspan="1" colspan="1">105</td><td align="left" rowspan="1" colspan="1">9.68E-05</td><td align="left" rowspan="1" colspan="1">0.013</td><td align="left" rowspan="1" colspan="1">T(A)((I,M))</td><td align="left" rowspan="1" colspan="1">V(M)((I))</td><td align="left" rowspan="1" colspan="1">2.01</td><td align="left" rowspan="1" colspan="1">2.42</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td align="left" rowspan="1" colspan="1">199</td><td align="left" rowspan="1" colspan="1">2.78E-04</td><td align="left" rowspan="1" colspan="1">0.023</td><td align="left" rowspan="1" colspan="1">A</td><td align="left" rowspan="1" colspan="1">S</td><td align="left" rowspan="1" colspan="1">2.88</td><td align="left" rowspan="1" colspan="1">2.50</td><td align="left" rowspan="1" colspan="1">A</td></tr><tr><td align="left" rowspan="1" colspan="1">475</td><td align="left" rowspan="1" colspan="1">5.46E-04</td><td align="left" rowspan="1" colspan="1">0.023</td><td align="left" rowspan="1" colspan="1">L((M))</td><td align="left" rowspan="1" colspan="1">M</td><td align="left" rowspan="1" colspan="1">2.51</td><td align="left" rowspan="1" colspan="1">3.50</td><td align="left" rowspan="1" colspan="1">L</td></tr><tr><td align="left" rowspan="1" colspan="1">493</td><td align="left" rowspan="1" colspan="1">1.78E-03</td><td align="left" rowspan="1" colspan="1">0.039</td><td align="left" rowspan="1" colspan="1">R((K))</td><td align="left" rowspan="1" colspan="1">K((R))</td><td align="left" rowspan="1" colspan="1">2.53</td><td align="left" rowspan="1" colspan="1">2.70</td><td align="left" rowspan="1" colspan="1">R</td></tr><tr><td align="left" rowspan="1" colspan="1">569</td><td align="left" rowspan="1" colspan="1">2.39E-03</td><td align="left" rowspan="1" colspan="1">0.048</td><td align="left" rowspan="1" colspan="1">T((A))</td><td align="left" rowspan="1" colspan="1">A((S))</td><td align="left" rowspan="1" colspan="1">2.51</td><td align="left" rowspan="1" colspan="1">2.69</td><td align="left" rowspan="1" colspan="1">T</td></tr><tr><td align="left" rowspan="1" colspan="1">613</td><td align="left" rowspan="1" colspan="1">1.11E-03</td><td align="left" rowspan="1" colspan="1">0.035</td><td align="left" rowspan="1" colspan="1">V(A)((I))</td><td align="left" rowspan="1" colspan="1">T(I,A)</td><td align="left" rowspan="1" colspan="1">2.33</td><td align="left" rowspan="1" colspan="1">1.82</td><td align="left" rowspan="1" colspan="1">V</td></tr><tr><td align="left" rowspan="1" colspan="1">627</td><td align="left" rowspan="1" colspan="1">1.20E-03</td><td align="left" rowspan="1" colspan="1">0.035</td><td align="left" rowspan="1" colspan="1">E(K)</td><td align="left" rowspan="1" colspan="1">K</td><td align="left" rowspan="1" colspan="1">2.20</td><td align="left" rowspan="1" colspan="1">2.89</td><td align="left" rowspan="1" colspan="1">E</td></tr><tr><td align="left" rowspan="1" colspan="1">661</td><td align="left" rowspan="1" colspan="1">5.91E-04</td><td align="left" rowspan="1" colspan="1">0.023</td><td align="left" rowspan="1" colspan="1">A(T)((V))</td><td align="left" rowspan="1" colspan="1">T((V))</td><td align="left" rowspan="1" colspan="1">2.28</td><td align="left" rowspan="1" colspan="1">2.51</td><td align="left" rowspan="1" colspan="1">A</td></tr><tr><td align="left" rowspan="1" colspan="1">682</td><td align="left" rowspan="1" colspan="1">1.48E-03</td><td align="left" rowspan="1" colspan="1">0.038</td><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">S(N)</td><td align="left" rowspan="1" colspan="1">2.60</td><td align="left" rowspan="1" colspan="1">1.94</td><td align="left" rowspan="1" colspan="1">G</td></tr><tr><td align="left" rowspan="1" colspan="1">684</td><td align="left" rowspan="1" colspan="1">7.63E-05</td><td align="left" rowspan="1" colspan="1">0.013</td><td align="left" rowspan="1" colspan="1">A((T))</td><td align="left" rowspan="1" colspan="1">S(T)</td><td align="left" rowspan="1" colspan="1">2.69</td><td align="left" rowspan="1" colspan="1">1.95</td><td align="left" rowspan="1" colspan="1">S</td></tr><tr><td align="left" rowspan="1" colspan="1">702</td><td align="left" rowspan="1" colspan="1">1.58E-03</td><td align="left" rowspan="1" colspan="1">0.038</td><td align="left" rowspan="1" colspan="1">K(R)</td><td align="left" rowspan="1" colspan="1">R</td><td align="left" rowspan="1" colspan="1">2.42</td><td align="left" rowspan="1" colspan="1">2.78</td><td align="left" rowspan="1" colspan="1">K</td></tr><tr><td align="left" rowspan="1" colspan="1">740</td><td align="left" rowspan="1" colspan="1">3.83E-04</td><td align="left" rowspan="1" colspan="1">0.023</td><td align="left" rowspan="1" colspan="1">D</td><td align="left" rowspan="1" colspan="1">D(N)</td><td align="left" rowspan="1" colspan="1">3.03</td><td align="left" rowspan="1" colspan="1">2.24</td><td align="left" rowspan="1" colspan="1">D</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt101"><label></label><p>Legend: Position: Sites in H1, H2, and H3 are identified with respect to their H3 positions. (Locations with negative positions represent an N-terminal insertion in H1 relative to H3. Similarly, location H1 271A represents an insertion between locations 271 and 272.) Numbering refer to HA1 unless specified HA2. Sites in N1 and N2 are identified with respect to their own indices. P: P value using LRT for non-homogeneous model (Model 3) compared with homogeneous model (Model 4). δ: Minimum FDR value required for this site to be identified. This should not be equated with the probability that this identification is a false positive. Residues Identities: Residue identity is shown unadorned if its equilibrium frequency is greater than 0.5, in single parentheses if its equilibrium frequency is between 0.1 and 0.5, and in double parentheses if its equilibrium frequency is between 0.01 and 0.1. Residues with equilibrium frequencies below 0.01 are not listed. Selective (Sel.) constraints represent the strength of the selective pressure in the two different hosts as defined in Equation 5. Cal-4: Identity of the amino acid at this location in a sample of the recent Swine flu outbreak (A/California/04/2009(H1N1)).</p></fn></table-wrap-foot></table-wrap><p>To assess the performance of the technique describe here, we simulated each one of the 264 variable sites in the PB2 gene ten times (2,640 simulations in total). All sites were simulated using the same fixed tree topology. The 22 ‘A’ and ‘B’ sites identified as undergoing selective constraint changes (FDR<0.20) were simulated under the non-homogeneous model, using the parameters obtained by optimizing model 3. Similarly, the 242 locations with no evidence for change in selective constraints were simulated under the homogeneous model (model 2). We then applied the analysis described above to identify locations in the synthetic datasets that had undergone changes in selective pressure. On average, we observed that 1.5% of the locations identified with FDR<0.05 were false positives (false positive rate of 0.08%); this increased to 3.6% (false positive rate of 0.2%) for FDR<0.20. This indicates that the FDR values are, at least for PB2, likely conservative. Of the 22 locations modelled with changing selective constraints, 12.9 were identified with FDR<0.05 (false negative rate of 41%), with 16.2 identified with FDR<0.20 (false negative rate of 26%). The 13 ‘A’ sites were identified more consistently, with 10.1 found with FDR<0.05 and 11.0 found with FDR<0.20. This suggests that there remain more locations undergoing changes in selective pressure than are being identified with the procedure described here.</p><p>Our approach relies on the prior construction of an appropriate phylogenetic tree. In order to estimate the effect of phylogenetic uncertainty, we repeated the analysis of the PB2 gene segment with ten different phylogenetic trees obtained through non-parametric bootstrapping. The 13 ‘A’ sites were identified on 79% of the bootstrap trees with FDR<0.05 and identified on 90% with FDR<0.20. 85% of the 22 ‘A’ and ‘B’ sites were similarly identified on the bootstrap trees with FDR<0.20. Conversely, the bootstrap trees identified on average 2% (with FDR<0.05) and 6% (with FDR<0.20) of alternative locations that were not identified on the original tree. These might be false positives for the alternative trees, suggesting a similar amount of false positives on the original tree. Some of these locations, however, may be locations with changes in selective constraints, and thus represent false negatives for the original tree; most of these locations would have been so identified with a higher FDR threshold of 0.50, although these points represent only about 12% of the otherwise unidentified locations.</p><p>We constructed a simple model to help explain the lack of statistically significant improvement with adding host-specific scaling factors. This was based on considering a protein site where two amino acids (<italic>A</italic> and <italic>B</italic>) are present, where an organism with residue <italic>B</italic> has a fitness equal to 1−<italic>s</italic> relative to an organism with residue <italic>A</italic>. We used Kimura's fixation rate theory <xref rid="pcbi.1000564-Crow1" ref-type="bibr">[43]</xref> to calculate the resulting substitution rates between <italic>A</italic> to <italic>B</italic>, and formulate these expressions in terms of a rate scaling factor ν and equilibrium frequencies π<italic><sub>A</sub></italic> and π<italic><sub>B</sub></italic> ( = 1−π<italic><sub>A</sub></italic>). We considered how ν, π<italic><sub>A</sub></italic> , and π<italic><sub>B</sub></italic> change as the relative fitness difference between <italic>A</italic> and <italic>B</italic> is altered. We also considered the overall rate at which substitutions occur in both directions, both for negative selection where the residues are at equilibrium (Γ<sub>−</sub>) as well as for positive selection (Γ<sub>+</sub>) where the location contains the unfavourable residue <italic>B</italic>. <xref ref-type="fig" rid="pcbi-1000564-g003">Figure 3</xref> shows the dependence of π<italic><sub>A</sub></italic>, π<italic><sub>B</sub></italic>, ν, Γ<sub>−</sub>, and Γ<sub>+</sub> (the latter three normalised by the mutation rate μ) on the relative fitness difference <italic>s</italic> (scaled by the effective population size <italic>N<sub>eff</sub></italic>). As shown, under conditions of negative selection, increasing fitness differences result in a decrease in the overall rate of substitutions, but an <italic>increase</italic> in the rate-scaling factor. There is a relatively weak dependence of ν on <italic>s</italic> as long as the latter is not large relative to 1/<italic>N</italic><sub>eff</sub>. Under conditions of positive selection, both quantities increase with larger fitness differences.</p><fig id="pcbi-1000564-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000564.g003</object-id><label>Figure 3</label><caption><title>Changing equilibrium frequencies and rates versus selective constraints.</title><p>Top: Dependence of rate scaling factor ν (solid line) and rate of substitutions for adaptive (positive) selection Γ<sub>+</sub> (dotted line) and purifying (negative) selection Γ<sub>−</sub> (dashed line), scaled by mutation rate μ, as a function of scaled selective disadvantage of residue <italic>B</italic> compared with residue <italic>A</italic> (2<italic>N</italic><sub>eff</sub><italic>s</italic>). Bottom: Equilibrium frequencies π<italic><sub>A</sub></italic> of <italic>A</italic> (solid line) and π<italic><sub>B</sub></italic> of <italic>B</italic> (dashed line) as a function of scaled selective disadvantage.</p></caption><graphic xlink:href="pcbi.1000564.g003"></graphic></fig><p>The theoretically predicted weak dependence of ν on selective pressure and the lack of statistical support for host-dependent values of this parameter indicate that ν is not a good measure of the degree of selective constraints. To generate a more appropriate measure, we calculated the relative entropy between the equilibrium frequencies and what would be expected under no selection, <bold>π<sup>0</sup></bold>, estimating the latter by averaging the amino acid frequencies over our entire database. This measure of selective constraint magnitudes for the various sites in avian and human hosts are presented in <xref ref-type="table" rid="pcbi-1000564-t001">Table 1</xref>, Supporting <xref ref-type="supplementary-material" rid="pcbi.1000564.s001">Table S1</xref>, and in <xref ref-type="fig" rid="pcbi-1000564-g004">Figure 4</xref>.</p><fig id="pcbi-1000564-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000564.g004</object-id><label>Figure 4</label><caption><title>Selective constraint strengths for viral sites in Avian and Human influenza.</title><p>Strength of selective constraints (measured as <italic>d</italic>, as described in Eqn. 5), for viral sites identified (FDR<0.05) as under different selective constraints in avian and human hosts. Colour coding refers to specific gene: HA (red), NA (blue), M1 (black), M2 (brown), NP (green), NS1 (orange), PA (cyan), PB1 (purple), PB2 (magenta). Selective sites are labelled.</p></caption><graphic xlink:href="pcbi.1000564.g004"></graphic></fig></sec><sec id="s3"><title>Discussion</title><p>As described in the <xref ref-type="sec" rid="s1">introduction</xref>, ignoring the underlying phylogenetic relationship often results in a gross over-estimation of statistical significance, as single evolutionary events are interpreted as a large number of independent measurements. Correspondingly, certain sites that have been identified by other methods that do not model the underlying phylogenetics lose their statistical significance when the phylogenetics is considered. For instance, site 271 in PB2 is identified as a significant site in three previous analyses <xref rid="pcbi.1000564-Chen1" ref-type="bibr">[22]</xref>–<xref rid="pcbi.1000564-Miotto1" ref-type="bibr">[24]</xref>; human viral sequences are most commonly alanine at this site, while avian viral sequences are predominantly threonine, although alanine also occurs. When each sequence is interpreted as an independent event, there is strong statistical support for host-specific amino acid distributions at this site. All of the alanines in the human lineage, however, can be explained by a single threonine to alanine substitution. In contrast, in the avian influenza there were at least three independent threonine to alanine substitutions. The single example of the substitution in human influenza is not significant given the relative frequency of this transition in avian influenza. Indeed, the more complex Model 3 incorporating host-dependent substitution rates a <italic>P</italic> value of 0.095 compared with Model 2 that assumes no such host-dependence, with an unimpressive false discovery rate of 0.48 after the correction for multiple hypothesis testing. More threonine to alanine substitutions in the human lineage, even if that meant more human sequences with a threonine at this site, would have provided more statistical support. The statistical support would also have been larger if the various avian strains with an alanine at this site represented the result of a single substitution.</p><p>The sites that are identified are those with a significant statistical signal given the available data; other sites might be undergoing shifts in selective constraints that are not detected for different reasons. As with all appropriate statistical methods applied to this problem, we require adequate evolutionary time and a suitable substitution rate for the substitution patterns to be detectable. In particular, there has to be sufficient evolutionary time in both the avian and human lineages for the parameters in the substitution models to be sufficiently well defined in each so that the differences in selective constraints are detectible. This will require longer evolutionary time when the selective constraint changes are smaller. As shown in the phylogenetic trees (Supporting <xref ref-type="supplementary-material" rid="pcbi.1000564.s002">Figures S1</xref>, <xref ref-type="supplementary-material" rid="pcbi.1000564.s003">S2</xref>, <xref ref-type="supplementary-material" rid="pcbi.1000564.s004">S3</xref>, <xref ref-type="supplementary-material" rid="pcbi.1000564.s005">S4</xref>, <xref ref-type="supplementary-material" rid="pcbi.1000564.s006">S5</xref>, <xref ref-type="supplementary-material" rid="pcbi.1000564.s007">S6</xref>, <xref ref-type="supplementary-material" rid="pcbi.1000564.s008">S7</xref>, <xref ref-type="supplementary-material" rid="pcbi.1000564.s009">S8</xref>, <xref ref-type="supplementary-material" rid="pcbi.1000564.s010">S9</xref>, <xref ref-type="supplementary-material" rid="pcbi.1000564.s011">S10</xref>, and <xref ref-type="supplementary-material" rid="pcbi.1000564.s012">S11</xref>), there is relatively little sequence evolution in the human H2 lineage; this is possibly the cause of the relatively few sites identified in this gene subtype. There are more H3 sequences, although most available avian H3 sequences are highly similar, reducing our ability to detect selective pressure changes in this gene subtype. In particular, we do not identify the H3 Q226L mutation whose importance has been determined experimentally, as the strict conservation of glutamine in the avian lineage is not highly informative given the lack of evolutionary divergence among the avian H3 sequences. Finally, the improvement in the log likelihood necessary for a given level of statistical significance is a function in the increase in the number of adjustable parameters between the two models, which is one minus the number of amino acids found in that location. Locations that are highly variable require more adjustable parameters, reducing the power of the likelihood ratio test. In particular, human H3 viruses contain glutamine, leucine, isoleucine, and valine at position 226, making identification of selective constraint changes at this location difficult.</p><p>The identified changes in selective constraints may not be the direct result of the host shift event. Selection constraint changes at one site might be a response to substitutions that occur at a different site, even if those changes were themselves the result of neutral drift. We have also assumed that the change in selection constraint occurs simultaneously with the host shift event. In reality this method has limited temporal resolution, and changes in the substitution rate occurring near the host shift event might also be identified.</p><p>We do not include ‘pre-selection’ in the model, that is, that the match between the avian sequence and the selective constraints in the human host does not influence the probability that that particular virus strain will undergo a shift to humans. This could be added to such a phylogenetic-based model by considering the probability that a host shift would occur on a given lineage as a function of the protein sequence. This would greatly increase the complexity of the model, increasing the number of adjustable parameters, reducing the statistical power of the method. It is important to note, however, that this would increase the number of false negatives, as these occurrences would look identical to founder effects. It is less likely that this process would produce false positives.</p><p>We have included information from the A/California/04/2009 (H1N1) sequences from the 2009 ‘Swine flu’ pandemic in <xref ref-type="table" rid="pcbi-1000564-t001">Table 1</xref>. This strain represents a reassortment of avian-like European swine genes (M1, M2, NA) with a triple-reassortment strain previously circulating in swine containing segments originally from the classical swine (NP, NS, HA), human (PB1), and avian lineages (PA, PB2) <xref rid="pcbi.1000564-Dawood1" ref-type="bibr">[44]</xref>. Considering the locations identified with a false discovery rate of 5%, most segments originally from classical or European swine (NA, M1, NP, NS1) mostly matched the human selective constraints, suggesting a similarity between the constraints in humans and swine. The exception is the HA gene, where many locations seemed to match the avian selective constraints despite its classical-swine origin, possibly reflecting the slow rate of antigenic change of the classical swine haemagglutinin <xref rid="pcbi.1000564-Sheerar1" ref-type="bibr">[45]</xref>,<xref rid="pcbi.1000564-Vincent1" ref-type="bibr">[46]</xref>. In the segments more recently from the avian lineages (PA, PB2), most locations more closely matched the avian constraints, while PB2 684 and PA 356 more closely matched the human. Interestingly, by comparing with avian sequences, it appears that PB2 A684S and PA K356R substitutions, both involving changes from an avian-like to a human-like amino acid, occurred in the interval between the host shift to swine and the subsequent transfer to humans, suggesting that these changes might be related to the ability of these viruses to infect humans.</p><sec id="s3a"><title>Changes in equilibrium frequencies versus changes in the rate-scaling factor</title><p>Most methods that look for changes in the substitution rates model this as changes in ν, the scaling parameter, or in the related ratio of synonymous to non-synonymous substitutions. In our analysis, we find that, when we allow the equilibrium frequencies <bold>π</bold> to vary, there is no statistically significant variation in ν. This seems initially counter-intuitive, as there are some sites where there seems to be substantial changes in the degree of conservation; in site 274 in N1, for instance, is almost universally tyrosine in avian viruses, while it varies between tyrosine, serine, and phenylalanine in human viruses. Yet the likelihood ratio test applied to this site rejects the inclusion of host-dependent scaling factors with a <italic>P</italic> value of 0.90, suggesting that the relationship between rate scaling factors and site variation are not simply related.</p><p>This observation motivated our simple model to try to gain insight into the relationship between equilibrium frequencies and rate scaling factors, by considering a protein site where two different amino acids, <italic>A</italic> and <italic>B</italic>, are found. We imagine that organisms possessing residue <italic>A</italic> at this location have a fitness advantage. Negative purifying selection would occur when the residues at this location are at their equilibrium value, while positive selection would occur when this location was filled by <italic>B</italic>, such as might occur when the selective pressure on the protein changes. By using Kimura's theory of fixation probability <xref rid="pcbi.1000564-Crow1" ref-type="bibr">[43]</xref>, we can calculate the values of the rate scaling factor ν, the overall rate of substitutions for purifying (Γ<sub>−</sub>) and positive selection (Γ<sub>+</sub>), and the equilibrium frequencies of <italic>A</italic> (π<italic><sub>A</sub></italic>) and <italic>B</italic> (π<italic><sub>B</sub></italic>), as a function of the different finesses provided to an organism with the two different possible amino acids at that location, as described in the <xref ref-type="sec" rid="s4">Methods</xref> section below. Normalised values of ν, Γ<sub>−</sub>, and Γ<sub>+</sub> are plotted as a function of 2<italic>N</italic><sub>eff</sub><italic>s</italic> in <xref ref-type="fig" rid="pcbi-1000564-g003">Figure 3</xref>. As shown, ν varies surprisingly little with <italic>s</italic> as long as <italic>s</italic> is not much more than 1/<italic>N</italic><sub>eff</sub>. This explains why including a host-specific ν never yielded statistically significant improvements with our data. When we consider adaptive substitutions, larger values of <italic>s</italic> correspond to higher selective constraints, larger values of ν, and faster evolution. The situation is quite different with purifying selection. As might be expected, larger values of <italic>s</italic> (corresponding to larger degree of purifying selection) result in a slower substitution rate, but this actually corresponds to <italic>larger</italic> values of ν. The reason why most phylogenetic programs use an inverted relationship, where larger values of ν correspond to faster substitution rates, is that they do not consider the value of <bold>π</bold> appropriate for each site. By assuming that the same values of <bold>π</bold> apply to all sites, a more extreme distribution of equilibrium frequencies, resulting in a decrease in the number of substitutions, is interpreted as a reduction in ν although this parameter is, in fact, increasing</p><p>The magnitude of the selective constraints for the various sites in avian and human hosts are presented in <xref ref-type="table" rid="pcbi-1000564-t001">Table 1</xref>, Supporting <xref ref-type="supplementary-material" rid="pcbi.1000564.s001">Table S1</xref>, and in <xref ref-type="fig" rid="pcbi-1000564-g004">Figure 4</xref>. It is interesting to note the number of positions under changing selection constraints where the magnitudes of the selection constraints are relatively constant. Such sites would be difficult to detect by looking for changes in the substitution rate, especially in cases where the distributions of amino acids found in the two hosts have significant overlap.</p><p>The methods described here are applicable for a wide range of problems involving changes in selective constraints. There are two particular factors, however, that make the technique especially well suited for influenza. Firstly, the branch along which the selective pressure changes can be identified <italic>a priori</italic>. Secondly, it is important to generate appropriate phylogenetic trees for the position under consideration. Generation of such trees can be complicated when there is incongruence between different locations. For influenza, incongruence between the various genomic segments results from the process of reassortment, where chimeric viruses containing genomic segments of different origin result from multiple infections. We are able to address this issue by considering each different genomic segment independently, constructing gene-specific phylogenetic trees. A more difficult problem is intra-gene homologous recombination, where different regions of a single genomic segment have different phylogenies. Such recombination is either extremely rare or non-existent in influenza (as well as other negative RNA viruses), and has never been observed experimentally <xref rid="pcbi.1000564-Boni1" ref-type="bibr">[47]</xref>–<xref rid="pcbi.1000564-Krasnitz1" ref-type="bibr">[49]</xref>.</p><p>We have assumed that the transitions from avian to human hosts did not go through an intermediate species, such as swine. There is no evidence of involvement of swine in the 1957 Asian flu and 1968 Hong Kong flu host shift events. Based on his analysis of the 1918 Spanish flu sequences and the relative timing of the 1918 influenza outbreaks in swine and humans, Taubenberger concluded that the Spanish flu transferred <italic>in toto</italic> from birds to humans and from humans to swine <xref rid="pcbi.1000564-Reid1" ref-type="bibr">[3]</xref>–<xref rid="pcbi.1000564-Taubenberger2" ref-type="bibr">[5]</xref>, although this conclusion has been challenged <xref rid="pcbi.1000564-Antonovics1" ref-type="bibr">[6]</xref>–<xref rid="pcbi.1000564-Smith1" ref-type="bibr">[9]</xref>. If an intermediate host species were involved, it would not be expected to affect the results if the selective constraints at any location in this intermediate host were to resemble either that of avian or human viruses, as this would only change the timing of the shift from one selective constraint to another. If there were an intermediate host and the selective constraints at some locations in this intermediate host were strong and substantially different from either avian or human viruses, the amount of evolutionary time in this intermediate host were sufficiently long, and the evolutionary time in humans sufficiently short so that the new equilibrium is not attained, the results of these calculations could be affected.</p><p>There are two other important assumptions made in this work. Firstly, we assume that the selective constraints in human and avian viruses are constant, and that each location can be considered independently. We do not consider, for instance, that there may be different selective constraints in low-pathogenic and high-pathogenic avian viruses, or that compensatory changes can occur elsewhere in the protein or even in other proteins. The observation (both here and experimentally <xref rid="pcbi.1000564-Connor1" ref-type="bibr">[12]</xref>–<xref rid="pcbi.1000564-Vines1" ref-type="bibr">[16]</xref>) that different hemagglutinin subtypes undergo different patterns of change of selective constraints indicates that this assumption is not strictly valid.</p></sec></sec><sec id="s4"><title>Methods</title><sec id="s4a"><title>Theory</title><p>For the following discussion we assume the evolution of a viral protein along a phylogenetic tree with two different host lineages, avian and human, where we consider the root of the tree to exist somewhere in the avian lineage. The evolution of amino acids in a site along a phylogenetic tree can be modelled as a continuous Markov process, described by a 20×20 substitution matrix <bold>Q</bold>. (Standard phylogenetic modelling techniques are described in <xref rid="pcbi.1000564-Yang2" ref-type="bibr">[50]</xref>.) In order to provide for time reversibility (that is, the expected number <italic>i</italic> to <italic>j</italic> transitions equalling the expected number of transitions from <italic>j</italic> to <italic>i</italic>), this is commonly represented as <inline-formula><inline-graphic xlink:href="pcbi.1000564.e001.jpg" mimetype="image"></inline-graphic></inline-formula> where <bold>S</bold> is a symmetric matrix representing the exchangeability of amino acids <italic>i</italic> and <italic>j</italic>, π<italic><sub>j</sub></italic> is the equilibrium frequency of amino acid <italic>j</italic> (<inline-formula><inline-graphic xlink:href="pcbi.1000564.e002.jpg" mimetype="image"></inline-graphic></inline-formula>) and ν is a scaling parameter that accounts for the overall rate of substitution at the site. <bold>S</bold> encodes the underlying codon structure as well as the relative similarities of the physicochemical properties of the amino acids, while the equilibrium frequencies represent the relative propensities for each of the amino acids at that site. We can calculate the likelihood of the data at this site given the model using Felsenstein's pruning algorithm <xref rid="pcbi.1000564-Felsenstein1" ref-type="bibr">[51]</xref>,<xref rid="pcbi.1000564-Felsenstein2" ref-type="bibr">[52]</xref>.</p><p>We first consider a standard substitution model where <bold>S</bold> and <bold>π</bold> are given by the WAG substitution matrix <xref rid="pcbi.1000564-Whelan1" ref-type="bibr">[41]</xref>, where each site in the set of proteins is characterised by a distinct substitution rate scaling factor ν whose value is determined by maximising the log likelihood given the sequence data at that site and the input phylogenetic tree. This we refer to as Model 1. We then considered the appropriateness of modelling each site in the set of proteins with a distinctive set of equilibrium amino acid frequencies <xref rid="pcbi.1000564-Bruno1" ref-type="bibr">[26]</xref>, what we refer to as single-site homogeneous Model 2. We adjust the values of <bold>π</bold> simultaneously with ν to maximise the likelihood. To avoid over parameterisation, we still use WAG <bold>S</bold> values for all sites. The tree topology is assumed fixed, and branch lengths are the same for all sites. In order to reduce the number of adjustable parameters, π<italic><sub>i</sub></italic> = 0 for any amino acids not found at that site. As the equilibrium frequencies of the amino acids not observed are set to zero, this results in an increase in the number of parameters equal to the number of amino acids present at that site minus one (due to the constraint that the equilibrium frequencies must sum to one). We then use the likelihood ratio test to see if site-dependent equilibrium frequencies can be justified with the data. As described in the <xref ref-type="sec" rid="s2">Results</xref> section, the site-dependence of the equilibrium frequencies could be justified for all sites.</p><p>Now let us imagine that upon inspection of the phylogenetic tree, we notice that amino acid preferences at a particular site seem different in the two host clades. We can incorporate this observation into our model by using two distinct <bold>Q</bold> matrices to describe the evolution of this site in the different hosts, as illustrated in <xref ref-type="fig" rid="pcbi-1000564-g002">Figure 2</xref>. For the reservoir avian host we write <inline-formula><inline-graphic xlink:href="pcbi.1000564.e003.jpg" mimetype="image"></inline-graphic></inline-formula> and for the new human host <inline-formula><inline-graphic xlink:href="pcbi.1000564.e004.jpg" mimetype="image"></inline-graphic></inline-formula> where <bold>π</bold> and <bold>π′</bold> represent the equilibrium amino acid frequencies at that site in avian and human viruses, respectively. (In principle we could also have <bold>S</bold> depend upon the host, but this would result in a large increase in the number of adjustable parameters. We will consider host-dependence of ν below.) The host shift event is defined as the midpoint of the branch connecting the common ancestor of the human viruses with its parent node. We can now calculate a new likelihood for this site using the same fixed topology, again adjusting <bold>π</bold>, <bold>π′</bold>, and ν to maximise the likelihood. We call this the single site non-homogeneous model, Model 3. Again, the increase in the number of adjustable parameters for Model 3 relative to Model 2 equals the number of amino acid types observed at that site minus one. Because the Model 2 is nested inside Model 3, we can again use the likelihood ratio test to test the hypothesis of different selective constraints in different hosts at that site.</p><p>In general, for a protein with <italic>N</italic> variable sites, we could repeat the procedure above for each site in the alignment, and perform <italic>N</italic> likelihood ratio tests. This would generate a list of those sites that show statistically different amino acid compositions, and hence distinctive selective constraints, in the different hosts. Following the calculation of the statistical significance for each site we can then use standard false discovery rate (FDR) methods to account for multiple hypothesis testing <xref rid="pcbi.1000564-Benjamini1" ref-type="bibr">[42]</xref>.</p><p>Finally, we consider if, in addition to host-dependent equilibrium frequencies, we also have statistical evidence for host-dependent rate scaling factors. We again use <inline-formula><inline-graphic xlink:href="pcbi.1000564.e005.jpg" mimetype="image"></inline-graphic></inline-formula> for the reservoir avian host but now use <inline-formula><inline-graphic xlink:href="pcbi.1000564.e006.jpg" mimetype="image"></inline-graphic></inline-formula> for the new human host where ν and ν′ represent the rate scaling factors at that site in avian and human viruses, respectfully. Again, Model 3 is nested inside Model 4 with an increase of one adjustable parameter, meaning that the statistical support for this extra factor can be evaluated with the likelihood ratio test. We do not observe support for this extra parameter in any of the sites after adjusting for multiple hypothesis testing.</p></sec><sec id="s4b"><title>Data and data analysis</title><p>Human and avian viral sequences were collected from the NCBI Influenza Virus Resource <xref rid="pcbi.1000564-Bao1" ref-type="bibr">[53]</xref>. Due to the frequency of reassortment, we cannot assume that the phylogenetic relationships for the various genomic segments are similar; they must be treated independently, including creating genetic-segment specific phylogenetic trees. The sequences for the various segments were treated as independent data sets, with separate datasets for the H1, H2, H3, N1, and N2 genes. Clusters of highly similar sequences (approximately >99.5%) were culled as to reduce the overall number of sequences to around 400 per dataset. It is common to find sporadic transmissions between avian, human, and other (e.g. swine) hosts; we eliminated all sequences resulting from such transmissions (e.g. human H5N1 sequences), leaving us with a single connected set of avian sequences and separate monophyletic human clades corresponding to the host shift events of 1918 (H1, N1, internal genes), 1957 (H2, N2, PB1), and 1968 (H3, PB1).</p><p>In order to generate more accurate phylogenetic trees, the culled sequences were aligned at the amino acid level (MUSCLE, <xref rid="pcbi.1000564-Edgar1" ref-type="bibr">[54]</xref>), with these alignments then used to create nucleotide codon alignments (PAL2NAL, <xref rid="pcbi.1000564-Suyama1" ref-type="bibr">[55]</xref>). The phylogenetic tree topologies were then created for the nucleotide data using PhyML (<xref rid="pcbi.1000564-Guindon2" ref-type="bibr">[56]</xref>; HKY85 model <xref rid="pcbi.1000564-Hasegawa1" ref-type="bibr">[57]</xref>, Gamma-distributed rates). The resulting trees are included as Supporting <xref ref-type="supplementary-material" rid="pcbi.1000564.s002">Figure S1</xref>–<xref ref-type="supplementary-material" rid="pcbi.1000564.s012">S11</xref>. Because amino acid distances are needed for the models developed here, branch lengths were then re-optimised for this fixed tree topology using the corresponding amino acid data (PAML <xref rid="pcbi.1000564-Yang3" ref-type="bibr">[58]</xref>,<xref rid="pcbi.1000564-Yang4" ref-type="bibr">[59]</xref>, WAG substitution matrix <xref rid="pcbi.1000564-Whelan1" ref-type="bibr">[41]</xref>, Gamma-distributed rates). The analysis was then performed with each gene set, based on the phylogenetic tree for the genomic segment in which the gene is located. A computer program written in Java that implements and optimises the various models described above is available from the authors.</p><p>The determinations of changes in selective constraints at each site is a separate hypothesis to be evaluated, so we must address the multiple-hypothesis question, that is, if we ask a suitably large number of statistical questions we are likely, at random, to obtain some statistically-significant results. We use the false discovery rate method, that is, specifying for each site the false positive rate that would have to be tolerated in order for that result to be statistically significant, following the Benjamini and Hochberg estimator <xref rid="pcbi.1000564-Benjamini1" ref-type="bibr">[42]</xref>. We first choose an acceptable false discovery rate δ. If <italic>P</italic>(<italic>k</italic>) is the <italic>k</italic>-th smallest <italic>P</italic> value for a set of <italic>n</italic> sites, we choose the largest value of <italic>k</italic> so that <inline-formula><inline-graphic xlink:href="pcbi.1000564.e007.jpg" mimetype="image"></inline-graphic></inline-formula>. As different genes are evolving in different circumstances, we would not expect the fraction of sites in each gene undergoing changes in selective constraints to be the same. Combining all of the genes together in one dataset would result in an increase in false positives for the genes with fewer changes in selective constraints, and an increase in false negatives for the genes with more changes in selective constraints. For this reason we analyse the false discovery rate for each gene individually. <xref ref-type="table" rid="pcbi-1000564-t001">Table 1</xref> and Supporting <xref ref-type="supplementary-material" rid="pcbi.1000564.s001">Table S1</xref> list, for each site, the smallest possible acceptable false discovery rate that would result in that site being labelled as statistically significant. These should not be interpreted as the probability that that given site is a false positive.</p></sec><sec id="s4c"><title>Parametric bootstrapping</title><p>Each site was simulated under the homogeneous (Model 2) and non-homogeneous (Model 3) models 10 times using the program <italic>Evolver</italic><xref rid="pcbi.1000564-Yang4" ref-type="bibr">[59]</xref> using the estimated tree topology and the WAG+F substitution matrix <xref rid="pcbi.1000564-Yang2" ref-type="bibr">[50]</xref>. For each site, the tree was scaled according to the site-specific estimated rate-scaling parameter ν. Simulation under the non-homogeneous model was performed in two steps: the avian part of the tree was simulated using a randomly generated root sequence following the avian equilibrium frequencies for that location. The avian subtree contained a host shift tip that served as the root of the human clade. The human subtree was then simulated according the human equilibrium frequencies using the simulated avian sequence at the host shift.</p></sec><sec id="s4d"><title>Alternative tree topologies</title><p>The PB2 sequence was bootstrapped 10 times and tree topology re-estimated for each boot sample. The homogeneous and non-homogeneous models were optimised for the observed data at each location, and the LRT was performed again for each one of the 10 new tree topologies so as to assess the effect of tree topology uncertainty on the identification of adaptive sites.</p></sec><sec id="s4e"><title>Simple model for relationship between equilibrium frequencies and scaling factors</title><p>Consider a protein site where two amino acids, <italic>A</italic> and <italic>B</italic>, are found. Let us imagine that that <italic>A</italic> is the more advantageous amino acid, that is, organisms with <italic>A</italic> at this site have a higher fitness, while organisms with <italic>B</italic> at this site has relative fitness <inline-formula><inline-graphic xlink:href="pcbi.1000564.e008.jpg" mimetype="image"></inline-graphic></inline-formula>. Let us also imagine that the mutation rate from <italic>A</italic> to <italic>B</italic> μ<italic><sub>AB</sub></italic> is equal to the reverse mutation rate μ<italic><sub>BA</sub></italic> = μ. We imagine a number of different lineages that have diverged, each with effective population size <italic>N</italic><sub>eff</sub>. Assuming that the mutation rate relative to the population is reasonably small, <italic>A</italic> or <italic>B</italic> will become fixed in each lineage. For haploid organisms, the probability that <italic>A</italic> would become fixed in a given lineage is given by <xref rid="pcbi.1000564-Crow1" ref-type="bibr">[43]</xref><disp-formula><graphic xlink:href="pcbi.1000564.e009.jpg" mimetype="image" position="float"></graphic><label>(1)</label></disp-formula>where we have recognised that this probability is simply the equilibrium frequency of <italic>A</italic> in the ensemble of diverged organisms, with <inline-formula><inline-graphic xlink:href="pcbi.1000564.e010.jpg" mimetype="image"></inline-graphic></inline-formula>.</p><p>The substitution rate of <italic>A</italic> by <italic>B</italic> is just the mutation rate <italic>μ</italic> times the fixation probability, given by Kimura's formula for small <italic>s</italic><xref rid="pcbi.1000564-Crow1" ref-type="bibr">[43]</xref>.<disp-formula><graphic xlink:href="pcbi.1000564.e011.jpg" mimetype="image" position="float"></graphic><label>(2)</label></disp-formula>We can compare these expressions with <inline-formula><inline-graphic xlink:href="pcbi.1000564.e012.jpg" mimetype="image"></inline-graphic></inline-formula> as used in phylogenetic analyses. As we are only dealing with two different residues, <inline-formula><inline-graphic xlink:href="pcbi.1000564.e013.jpg" mimetype="image"></inline-graphic></inline-formula> is a simple multiplicative constant and can be set equal to one, resulting in <inline-formula><inline-graphic xlink:href="pcbi.1000564.e014.jpg" mimetype="image"></inline-graphic></inline-formula>. Equating these two expressions for <italic>Q<sub>BA</sub></italic> and solving for ν yields<disp-formula><graphic xlink:href="pcbi.1000564.e015.jpg" mimetype="image" position="float"></graphic><label>(3)</label></disp-formula><disp-formula><graphic xlink:href="pcbi.1000564.e016.jpg" mimetype="image" position="float"></graphic><label>(4)</label></disp-formula></p><p>Similar results are obtained, as would be expected, when we express <inline-formula><inline-graphic xlink:href="pcbi.1000564.e017.jpg" mimetype="image"></inline-graphic></inline-formula>.</p><p>We can now consider the cases of neutral, adaptive (positive), and purifying (negative) selection. Neutral selection is simply the case when <italic>N</italic><sub>eff</sub><italic>s</italic> is small and <inline-formula><inline-graphic xlink:href="pcbi.1000564.e018.jpg" mimetype="image"></inline-graphic></inline-formula>. For both neutral and negative selection, we can consider the overall rate at which substitutions occur, given by <inline-formula><inline-graphic xlink:href="pcbi.1000564.e019.jpg" mimetype="image"></inline-graphic></inline-formula>, which is equal to <italic>N</italic><sub>eff</sub> μ in the case of neutral selection. Positive selection involves the situation where we are not at equilibrium, but rather, at least in this case, we have the less-fit residue occupying the given position. In this case, assuming again that <italic>A</italic> is the favoured residue, <inline-formula><inline-graphic xlink:href="pcbi.1000564.e020.jpg" mimetype="image"></inline-graphic></inline-formula>.</p></sec><sec id="s4f"><title>Characterising the magnitude of selective constraints</title><p>We characterise the selection constraints by how far the equilibrium amino acid frequencies <bold>π</bold> differ from what would be expected under no selection <bold>π<sup>0</sup></bold> through the relative entropy, defined as<disp-formula><graphic xlink:href="pcbi.1000564.e021.jpg" mimetype="image" position="float"></graphic><label>(5)</label></disp-formula>which is, as is desired, zero when <bold>π</bold> equals <bold>π<sup>0</sup></bold>. Unfortunately, it is difficult to estimate <bold>π<sup>0</sup></bold>, as there is little of the virus genome that is not under some degree of selective constraints. We estimate <bold>π<sup>0</sup></bold> by averaging the amino acid frequencies over our entire database, with the expectation that specific selection constraints will, at least approximately, average out.</p></sec></sec><sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title><supplementary-material content-type="local-data" id="pcbi.1000564.s001"><label>Table S1</label><caption><p>Sites identified as undergoing changes in selective pressure with a false discovery rate of 0.20. Position: Sites in H1, H2, and H3 are identified with respect to their H3 positions. Sites in N1 and N2 are identified with respect to their own indices. P: <italic>P</italic> value using LRT for non-homogeneous model (Model 2) compared with homogeneous model (Model 3). delta: Minimum FDR value required for this site to be identified. This should not be equated with the probability that this identification is a false positive. Residues Identities: Residue identity is shown unadorned if its equilibrium frequency is greater than 0.5, in single parentheses if its equilibrium frequency is between 0.1 and 0.5, and in double parentheses if its equilibrium frequency is between 0.01 and 0.1. Residues with equilibrium frequencies below 0.01 are not listed. Selective constraints represent the strength of the selective pressure in the two different hosts as defined in Equation 5 in the text. Highlighted locations are significant with a false discovery rate of 0.05.</p><p>(0.11 MB XLS)</p></caption><media xlink:href="pcbi.1000564.s001.xls"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pcbi.1000564.s002"><label>Figure S1</label><caption><p>Phylogenetic tree of HA genomic segment for subtype H1. Avian section of the tree is in black, human in red.</p><p>(2.81 MB PDF)</p></caption><media xlink:href="pcbi.1000564.s002.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pcbi.1000564.s003"><label>Figure S2</label><caption><p>Phylogenetic tree of HA genomic segment for subtype H2. Avian section of the tree is in black, human in red.</p><p>(0.71 MB PDF)</p></caption><media xlink:href="pcbi.1000564.s003.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pcbi.1000564.s004"><label>Figure S3</label><caption><p>Phylogenetic tree of HA genomic segment for subtype H3. Avian section of the tree is in black, human in red.</p><p>(2.94 MB PDF)</p></caption><media xlink:href="pcbi.1000564.s004.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pcbi.1000564.s005"><label>Figure S4</label><caption><p>Phylogenetic tree of NA genomic segment for subtype N1. Avian section of the tree is in black, human in red.</p><p>(3.72 MB PDF)</p></caption><media xlink:href="pcbi.1000564.s005.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pcbi.1000564.s006"><label>Figure S5</label><caption><p>Phylogenetic tree of NA genomic segment for subtype N2. Avian section of the tree is in black, human in red.</p><p>(3.29 MB PDF)</p></caption><media xlink:href="pcbi.1000564.s006.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pcbi.1000564.s007"><label>Figure S6</label><caption><p>Phylogenetic tree of MP genomic segment. Avian section of the tree is in black, human in red.</p><p>(3.02 MB PDF)</p></caption><media xlink:href="pcbi.1000564.s007.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pcbi.1000564.s008"><label>Figure S7</label><caption><p>Phylogenetic tree of NS genomic segment. Avian section of the tree is in black, human in red.</p><p>(2.94 MB PDF)</p></caption><media xlink:href="pcbi.1000564.s008.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pcbi.1000564.s009"><label>Figure S8</label><caption><p>Phylogenetic tree of NP genomic segment. Avian section of the tree is in black, human in red.</p><p>(3.28 MB PDF)</p></caption><media xlink:href="pcbi.1000564.s009.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pcbi.1000564.s010"><label>Figure S9</label><caption><p>Phylogenetic tree of PA genomic segment. Avian section of the tree is in black, human in red.</p><p>(3.25 MB PDF)</p></caption><media xlink:href="pcbi.1000564.s010.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pcbi.1000564.s011"><label>Figure S10</label><caption><p>Phylogenetic tree of PB1 genomic segment. Avian section of the tree is in black, human in red.</p><p>(3.03 MB PDF)</p></caption><media xlink:href="pcbi.1000564.s011.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pcbi.1000564.s012"><label>Figure S11</label><caption><p>Phylogenetic tree of PB2 genomic segment. Avian section of the tree is in black, human in red.</p><p>(3.21 MB PDF)</p></caption><media xlink:href="pcbi.1000564.s012.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><ack><p>We would like to thank Seena Shah, John McCauley, Ziheng Yang, and Ian Holmes for advice and helpful comments. Michael T. Flanagan provided the Java optimisation library used in this work.</p></ack><fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>Funding was provided by the NIMR (MRC), the Wellcome Trust, and the European FP6 FLUPOL contract (SPSB-CT-2007-044263). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group><ref-list><title>References</title><ref id="pcbi.1000564-Webster1"><label>1</label><element-citation publication-type="journal"><person-group person-group-type="author"><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></person-group><year>1992</year><article-title>Evolution and ecology of influenza A viruses.</article-title><source>Microbiol Rev</source><volume>56</volume><fpage>152</fpage><lpage>179</lpage><pub-id pub-id-type="pmid">1579108</pub-id></element-citation></ref><ref id="pcbi.1000564-Schafer1"><label>2</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schafer</surname><given-names>JR</given-names></name><name><surname>Kawaoka</surname><given-names>Y</given-names></name><name><surname>Bean</surname><given-names>WJ</given-names></name><name><surname>Suss</surname><given-names>J</given-names></name><name><surname>Senne</surname><given-names>D</given-names></name><etal></etal></person-group><year>1993</year><article-title>Origin of the pandemic 1957 H2 influenza A virus and the persistence of its possible progenitors in the avian reservoir.</article-title><source>Virology</source><volume>194</volume><fpage>781</fpage><lpage>788</lpage><pub-id pub-id-type="pmid">7684877</pub-id></element-citation></ref><ref id="pcbi.1000564-Reid1"><label>3</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reid</surname><given-names>AH</given-names></name><name><surname>Taubenberger</surname><given-names>JK</given-names></name><name><surname>Fanning</surname><given-names>TG</given-names></name></person-group><year>2004</year><article-title>Evidence of an absence: the genetic origins of the 1918 pandemic influenza virus.</article-title><source>Nat Rev Microbiol</source><volume>2</volume><fpage>909</fpage><lpage>914</lpage><pub-id pub-id-type="pmid">15494747</pub-id></element-citation></ref><ref id="pcbi.1000564-Taubenberger1"><label>4</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Taubenberger</surname><given-names>JK</given-names></name></person-group><year>2006</year><article-title>The origin and virulence of the 1918 “Spanish” influenza virus.</article-title><source>Proceedings of the American Philosophical Society</source><volume>150</volume><fpage>86</fpage><lpage>112</lpage><pub-id pub-id-type="pmid">17526158</pub-id></element-citation></ref><ref id="pcbi.1000564-Taubenberger2"><label>5</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Taubenberger</surname><given-names>JK</given-names></name><name><surname>Reid</surname><given-names>AH</given-names></name><name><surname>Lourens</surname><given-names>RM</given-names></name><name><surname>Wang</surname><given-names>R</given-names></name><name><surname>Jin</surname><given-names>G</given-names></name><etal></etal></person-group><year>2005</year><article-title>Characterization of the 1918 influenza virus polymerase genes.</article-title><source>Nature</source><volume>437</volume><fpage>889</fpage><lpage>893</lpage><pub-id pub-id-type="pmid">16208372</pub-id></element-citation></ref><ref id="pcbi.1000564-Antonovics1"><label>6</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Antonovics</surname><given-names>J</given-names></name><name><surname>Hood</surname><given-names>ME</given-names></name><name><surname>Baker</surname><given-names>CH</given-names></name></person-group><year>2006</year><article-title>Molecular virology: was the 1918 flu avian in origin?</article-title><source>Nature</source><volume>440</volume><fpage>E9</fpage><comment>discussion E9–10</comment><pub-id pub-id-type="pmid">16641950</pub-id></element-citation></ref><ref id="pcbi.1000564-dosReis1"><label>7</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>dos Reis</surname><given-names>M</given-names></name><name><surname>Hay</surname><given-names>AJ</given-names></name><name><surname>Goldstein</surname><given-names>RA</given-names></name></person-group><year>2009</year><article-title>Using non-homogeneous models of nucleotide substitution to identify host shift events: application to the origin of the 1918 ‘Spanish’ influenza pandemic virus.</article-title><source>J Mol Evol</source><comment>DOI 10.1007/s00239-009-9282-x</comment></element-citation></ref><ref id="pcbi.1000564-Gibbs1"><label>8</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gibbs</surname><given-names>M</given-names></name><name><surname>Gibbs</surname><given-names>A</given-names></name></person-group><year>2006</year><article-title>Molecular virology: Was the 1918 pandemic caused by a bird flu?</article-title><source>Nature</source><volume>440</volume><fpage>E8</fpage><lpage>E8</lpage><pub-id pub-id-type="pmid">16641948</pub-id></element-citation></ref><ref id="pcbi.1000564-Smith1"><label>9</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname><given-names>GJD</given-names></name><name><surname>Bahl</surname><given-names>J</given-names></name><name><surname>Vijaykrishna</surname><given-names>D</given-names></name><name><surname>Zhang</surname><given-names>JZ</given-names></name><name><surname>Poon</surname><given-names>LLM</given-names></name><etal></etal></person-group><year>2009</year><article-title>Dating the emergence of pandemic influenza viruses.</article-title><source>Proc Natl Acad Sci U S A</source><volume>106</volume><fpage>11709</fpage><lpage>11712</lpage><pub-id pub-id-type="pmid">19597152</pub-id></element-citation></ref><ref id="pcbi.1000564-Lin1"><label>10</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lin</surname><given-names>YP</given-names></name><name><surname>Shaw</surname><given-names>M</given-names></name><name><surname>Gregory</surname><given-names>V</given-names></name><name><surname>Cameron</surname><given-names>K</given-names></name><name><surname>Lim</surname><given-names>W</given-names></name><etal></etal></person-group><year>2000</year><article-title>Avian-to-human transmission of H9N2 subtype influenza A viruses: relationship between H9N2 and H5N1 human isolates.</article-title><source>Proc Natl Acad Sci U S A</source><volume>97</volume><fpage>9654</fpage><lpage>9658</lpage><pub-id pub-id-type="pmid">10920197</pub-id></element-citation></ref><ref id="pcbi.1000564-Gambaryan1"><label>11</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gambaryan</surname><given-names>AS</given-names></name><name><surname>Matrosovich</surname><given-names>MN</given-names></name><name><surname>Yamnikova</surname><given-names>SS</given-names></name><name><surname>Lvov</surname><given-names>DK</given-names></name><name><surname>Robertson</surname><given-names>JS</given-names></name><etal></etal></person-group><year>2002</year><article-title>Differences in receptor specificity between the influenza viruses of duck, chicken, and human.</article-title><source>Molecular Biology</source><volume>36</volume><fpage>429</fpage><lpage>435</lpage></element-citation></ref><ref id="pcbi.1000564-Connor1"><label>12</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Connor</surname><given-names>RJ</given-names></name><name><surname>Kawaoka</surname><given-names>Y</given-names></name><name><surname>Webster</surname><given-names>RG</given-names></name><name><surname>Paulson</surname><given-names>JC</given-names></name></person-group><year>1994</year><article-title>Receptor Specificity in Human, Avian, and Equine H2 and H3 Influenza-Virus Isolates.</article-title><source>Virology</source><volume>205</volume><fpage>17</fpage><lpage>23</lpage><pub-id pub-id-type="pmid">7975212</pub-id></element-citation></ref><ref id="pcbi.1000564-Matrosovich1"><label>13</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Matrosovich</surname><given-names>M</given-names></name><name><surname>Tuzikov</surname><given-names>A</given-names></name><name><surname>Bovin</surname><given-names>N</given-names></name><name><surname>Gambaryan</surname><given-names>A</given-names></name><name><surname>Klimov</surname><given-names>A</given-names></name><etal></etal></person-group><year>2000</year><article-title>Early alterations of the receptor-binding properties of H1, H2, and H3 avian influenza virus hemagglutinins after their introduction into mammals.</article-title><source>Journal of Virology</source><volume>74</volume><fpage>8502</fpage><lpage>8512</lpage><pub-id pub-id-type="pmid">10954551</pub-id></element-citation></ref><ref id="pcbi.1000564-Nobusawa1"><label>14</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nobusawa</surname><given-names>E</given-names></name><name><surname>Aoyama</surname><given-names>T</given-names></name><name><surname>Kato</surname><given-names>H</given-names></name><name><surname>Suzuki</surname><given-names>Y</given-names></name><name><surname>Tateno</surname><given-names>Y</given-names></name><etal></etal></person-group><year>1991</year><article-title>Comparison of Complete Amino-Acid-Sequences and Receptor-Binding Properties among 13 Serotypes of Hemagglutinins of Influenza a-Viruses.</article-title><source>Virology</source><volume>182</volume><fpage>475</fpage><lpage>485</lpage><pub-id pub-id-type="pmid">2024485</pub-id></element-citation></ref><ref id="pcbi.1000564-Rogers1"><label>15</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rogers</surname><given-names>GN</given-names></name><name><surname>Paulson</surname><given-names>JC</given-names></name><name><surname>Daniels</surname><given-names>RS</given-names></name><name><surname>Skehel</surname><given-names>JJ</given-names></name><name><surname>Wilson</surname><given-names>IA</given-names></name><etal></etal></person-group><year>1983</year><article-title>Single Amino-Acid Substitutions in Influenza Hemagglutinin Change Receptor-Binding Specificity.</article-title><source>Nature</source><volume>304</volume><fpage>76</fpage><lpage>78</lpage><pub-id pub-id-type="pmid">6191220</pub-id></element-citation></ref><ref id="pcbi.1000564-Vines1"><label>16</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vines</surname><given-names>A</given-names></name><name><surname>Wells</surname><given-names>K</given-names></name><name><surname>Matrosovich</surname><given-names>M</given-names></name><name><surname>Castrucci</surname><given-names>MR</given-names></name><name><surname>Ito</surname><given-names>T</given-names></name><etal></etal></person-group><year>1998</year><article-title>The role of influenza a virus hemagglutinin residues 226 and 228 in receptor specificity and host range restriction.</article-title><source>Journal of Virology</source><volume>72</volume><fpage>7626</fpage><lpage>7631</lpage><pub-id pub-id-type="pmid">9696865</pub-id></element-citation></ref><ref id="pcbi.1000564-Baigent1"><label>17</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Baigent</surname><given-names>SJ</given-names></name><name><surname>McCauley</surname><given-names>JW</given-names></name></person-group><year>2003</year><article-title>Influenza type A in humans, mammals and birds: determinants of virus virulence, host-range and interspecies transmission.</article-title><source>Bioessays</source><volume>25</volume><fpage>657</fpage><lpage>671</lpage><pub-id pub-id-type="pmid">12815721</pub-id></element-citation></ref><ref id="pcbi.1000564-Naffakh1"><label>18</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Naffakh</surname><given-names>N</given-names></name><name><surname>Tomoiu</surname><given-names>A</given-names></name><name><surname>Rameix-Welti</surname><given-names>MA</given-names></name><name><surname>van der Werf</surname><given-names>S</given-names></name></person-group><year>2008</year><article-title>Host Restriction of Avian Influenza Viruses at the Level of the Ribonucleoproteins.</article-title><source>Annual Review of Microbiology</source><volume>62</volume><fpage>403</fpage><lpage>424</lpage></element-citation></ref><ref id="pcbi.1000564-Hatta1"><label>19</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hatta</surname><given-names>M</given-names></name><name><surname>Gao</surname><given-names>P</given-names></name><name><surname>Halfmann</surname><given-names>P</given-names></name><name><surname>Kawaoka</surname><given-names>Y</given-names></name></person-group><year>2001</year><article-title>Molecular basis for high virulence of Hong Kong H5N1 influenza A viruses.</article-title><source>Science</source><volume>293</volume><fpage>1840</fpage><lpage>1842</lpage><pub-id pub-id-type="pmid">11546875</pub-id></element-citation></ref><ref id="pcbi.1000564-Steel1"><label>20</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Steel</surname><given-names>J</given-names></name><name><surname>Lowen</surname><given-names>A</given-names></name><name><surname>Mubareka</surname><given-names>S</given-names></name><name><surname>Palese</surname><given-names>P</given-names></name><name><surname>Baric</surname><given-names>R</given-names></name></person-group><year>2009</year><article-title>Transmission of Influenza Virus in a Mammalian Host Is Increased by PB2 Amino Acids 627K or 627E/701N.</article-title><source>PLoS Pathog</source><volume>5</volume><fpage>e1000252</fpage><pub-id pub-id-type="pmid">19119420</pub-id></element-citation></ref><ref id="pcbi.1000564-Subbarao1"><label>21</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Subbarao</surname><given-names>EK</given-names></name><name><surname>London</surname><given-names>W</given-names></name><name><surname>Murphy</surname><given-names>BR</given-names></name></person-group><year>1993</year><article-title>A Single Amino-Acid in the Pb2-Gene of Influenza-a Virus Is a Determinant of Host Range.</article-title><source>Journal of Virology</source><volume>67</volume><fpage>1761</fpage><lpage>1764</lpage><pub-id pub-id-type="pmid">8445709</pub-id></element-citation></ref><ref id="pcbi.1000564-Chen1"><label>22</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname><given-names>GW</given-names></name><name><surname>Chang</surname><given-names>SC</given-names></name><name><surname>Mok</surname><given-names>CK</given-names></name><name><surname>Lo</surname><given-names>YL</given-names></name><name><surname>Kung</surname><given-names>YN</given-names></name></person-group><year>2006</year><article-title>Genomic signatures of human versus avian influenza A viruses.</article-title><source>Emerg Infect Dis</source></element-citation></ref><ref id="pcbi.1000564-Finkelstein1"><label>23</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Finkelstein</surname><given-names>D</given-names></name><name><surname>Mukatira</surname><given-names>S</given-names></name><name><surname>Mehta</surname><given-names>P</given-names></name><name><surname>Obenauer</surname><given-names>J</given-names></name><name><surname>Su</surname><given-names>X</given-names></name><etal></etal></person-group><year>2007</year><article-title>Persistent Host Markers in Pandemic and H5N1 Influenza Viruses.</article-title><source>J Virol</source><volume>81</volume><fpage>10292</fpage><lpage>10299</lpage><pub-id pub-id-type="pmid">17652405</pub-id></element-citation></ref><ref id="pcbi.1000564-Miotto1"><label>24</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Miotto</surname><given-names>O</given-names></name><name><surname>Heiny</surname><given-names>A</given-names></name><name><surname>Tan</surname><given-names>T</given-names></name><name><surname>August</surname><given-names>J</given-names></name><name><surname>Brusic</surname><given-names>V</given-names></name></person-group><year>2008</year><article-title>Identification of human-to-human transmissibility factors in PB2 proteins of influenza A by large-scale mutual information analysis.</article-title><source>BMC Bioinformatics</source><volume>9</volume><fpage>S18</fpage></element-citation></ref><ref id="pcbi.1000564-Mayr1"><label>25</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Mayr</surname><given-names>E</given-names></name></person-group><year>1942</year><source>Systematics and the Origin of Species</source><publisher-loc>New York</publisher-loc><publisher-name>Columbia University Press</publisher-name></element-citation></ref><ref id="pcbi.1000564-Bruno1"><label>26</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bruno</surname><given-names>W</given-names></name></person-group><year>1996</year><article-title>Modeling residue usage in aligned protein sequences via maximum likelihood.</article-title><source>Molecular Biology and Evolution</source><volume>13</volume><fpage>1368</fpage><lpage>1374</lpage><pub-id pub-id-type="pmid">8952081</pub-id></element-citation></ref><ref id="pcbi.1000564-FreireMaia1"><label>27</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Freire-Maia</surname><given-names>N</given-names></name></person-group><year>1979</year><article-title>Neutral hypothesis is Darwinian.</article-title><source>Annals of Human Genetics</source><volume>42</volume><fpage>531</fpage><pub-id pub-id-type="pmid">382974</pub-id></element-citation></ref><ref id="pcbi.1000564-Blouin1"><label>28</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Blouin</surname><given-names>C</given-names></name><name><surname>Boucher</surname><given-names>Y</given-names></name><name><surname>Roger</surname><given-names>AJ</given-names></name></person-group><year>2003</year><article-title>Inferring functional constraints and divergence in protein families using 3D mapping of phylogenetic information.</article-title><source>Nucleic Acids Research</source><volume>31</volume><fpage>790</fpage><lpage>797</lpage><pub-id pub-id-type="pmid">12527789</pub-id></element-citation></ref><ref id="pcbi.1000564-Dorman1"><label>29</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dorman</surname><given-names>K</given-names></name></person-group><year>2007</year><article-title>Identifying dramatic selection shifts in phylogenetic trees.</article-title><source>BMC Evol Biol</source><volume>7</volume><fpage>S10</fpage></element-citation></ref><ref id="pcbi.1000564-Gu1"><label>30</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gu</surname><given-names>X</given-names></name></person-group><year>1999</year><article-title>Statistical methods for testing functional divergence after gene duplication.</article-title><source>Molecular Biology and Evolution</source><volume>16</volume><fpage>1664</fpage><lpage>1674</lpage><pub-id pub-id-type="pmid">10605109</pub-id></element-citation></ref><ref id="pcbi.1000564-Gu2"><label>31</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gu</surname><given-names>X</given-names></name></person-group><year>2001</year><article-title>Maximum-likelihood approach for gene family evolution under functional divergence.</article-title><source>Molecular Biology and Evolution</source><volume>18</volume><fpage>453</fpage><lpage>464</lpage><pub-id pub-id-type="pmid">11264396</pub-id></element-citation></ref><ref id="pcbi.1000564-Gu3"><label>32</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gu</surname><given-names>X</given-names></name><name><surname>Fu</surname><given-names>Y-X</given-names></name><name><surname>Li</surname><given-names>W-H</given-names></name></person-group><year>1995</year><article-title>Maximum Likelihood Estimation of the Heterogeneity of Substitution Rate among Nucleotide Sites.</article-title><source>Mol Biol Evol</source><volume>12</volume><fpage>546</fpage><lpage>557</lpage><pub-id pub-id-type="pmid">7659011</pub-id></element-citation></ref><ref id="pcbi.1000564-Knudsen1"><label>33</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Knudsen</surname><given-names>B</given-names></name><name><surname>Miyamoto</surname><given-names>MM</given-names></name></person-group><year>2001</year><article-title>A likelihood ratio test for evolutionary rate shifts and functional divergence among proteins.</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>98</volume><fpage>14512</fpage><lpage>14517</lpage><pub-id pub-id-type="pmid">11734650</pub-id></element-citation></ref><ref id="pcbi.1000564-Penn1"><label>34</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Penn</surname><given-names>O</given-names></name><name><surname>Stern</surname><given-names>A</given-names></name><name><surname>Rubinstein</surname><given-names>N</given-names></name><name><surname>Dutheil</surname><given-names>J</given-names></name><name><surname>Bacharach</surname><given-names>E</given-names></name><etal></etal></person-group><year>2008</year><article-title>Evolutionary Modeling of Rate Shifts Reveals Specificity Determinants in HIV-1 Subtypes.</article-title><source>PLoS Comput Biol</source><volume>4</volume><fpage>e1000214</fpage><pub-id pub-id-type="pmid">18989394</pub-id></element-citation></ref><ref id="pcbi.1000564-Pupko1"><label>35</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pupko</surname><given-names>T</given-names></name><name><surname>Galtier</surname><given-names>N</given-names></name></person-group><year>2002</year><article-title>A covarion-based method for detecting molecular adaptation: application to the evolution of primate mitochondrial genomes.</article-title><source>Proceedings of the Royal Society of London Series B-Biological Sciences</source><volume>269</volume><fpage>1313</fpage><lpage>1316</lpage></element-citation></ref><ref id="pcbi.1000564-Guindon1"><label>36</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Guindon</surname><given-names>S</given-names></name><name><surname>Rodrigo</surname><given-names>AG</given-names></name><name><surname>Dyer</surname><given-names>KA</given-names></name><name><surname>Huelsenbeck</surname><given-names>JP</given-names></name></person-group><year>2004</year><article-title>Modeling the site-specific variation of selection patterns along lineages.</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>101</volume><fpage>12957</fpage><lpage>12962</lpage><pub-id pub-id-type="pmid">15326304</pub-id></element-citation></ref><ref id="pcbi.1000564-Yang1"><label>37</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname><given-names>ZH</given-names></name><name><surname>Nielsen</surname><given-names>R</given-names></name></person-group><year>2002</year><article-title>Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages.</article-title><source>Molecular Biology and Evolution</source><volume>19</volume><fpage>908</fpage><lpage>917</lpage><pub-id pub-id-type="pmid">12032247</pub-id></element-citation></ref><ref id="pcbi.1000564-Zhang1"><label>38</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname><given-names>JZ</given-names></name><name><surname>Nielsen</surname><given-names>R</given-names></name><name><surname>Yang</surname><given-names>ZH</given-names></name></person-group><year>2005</year><article-title>Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level.</article-title><source>Molecular Biology and Evolution</source><volume>22</volume><fpage>2472</fpage><lpage>2479</lpage><pub-id pub-id-type="pmid">16107592</pub-id></element-citation></ref><ref id="pcbi.1000564-Forsberg1"><label>39</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Forsberg</surname><given-names>R</given-names></name><name><surname>Christiansen</surname><given-names>FB</given-names></name></person-group><year>2003</year><article-title>A codon-based model of host-specific selection in parasites, with an application to the influenza A virus.</article-title><source>Mol Biol Evol</source><volume>20</volume><fpage>1252</fpage><lpage>1259</lpage><pub-id pub-id-type="pmid">12777510</pub-id></element-citation></ref><ref id="pcbi.1000564-Blackburne1"><label>40</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Blackburne</surname><given-names>BP</given-names></name><name><surname>Hay</surname><given-names>AJ</given-names></name><name><surname>Goldstein</surname><given-names>RA</given-names></name></person-group><year>2008</year><article-title>Changing Selective Pressure during Antigenic Changes in Human Influenza H3.</article-title><source>PLoS Pathogens</source><volume>4</volume><fpage>e1000058</fpage><pub-id pub-id-type="pmid">18451985</pub-id></element-citation></ref><ref id="pcbi.1000564-Whelan1"><label>41</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Whelan</surname><given-names>S</given-names></name><name><surname>Goldman</surname><given-names>N</given-names></name></person-group><year>2001</year><article-title>A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach Molecular Biology and Evolution</article-title><volume>18</volume><fpage>691</fpage><lpage>699</lpage></element-citation></ref><ref id="pcbi.1000564-Benjamini1"><label>42</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Benjamini</surname><given-names>Y</given-names></name><name><surname>Hochberg</surname><given-names>T</given-names></name></person-group><year>1995</year><article-title>Controlling the false discovery rate: A practical and powerful approach to multiple testing.</article-title><source>J Royal Stat Soc B</source><volume>85</volume><fpage>289</fpage><lpage>300</lpage></element-citation></ref><ref id="pcbi.1000564-Crow1"><label>43</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Crow</surname><given-names>JF</given-names></name><name><surname>Kimura</surname><given-names>M</given-names></name></person-group><year>1970</year><source>An introduction to population genetics theory</source><publisher-loc>New York</publisher-loc><publisher-name>Harper & Row</publisher-name></element-citation></ref><ref id="pcbi.1000564-Dawood1"><label>44</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dawood</surname><given-names>FS</given-names></name><name><surname>Jain</surname><given-names>S</given-names></name><name><surname>Finelli</surname><given-names>L</given-names></name><name><surname>Shaw</surname><given-names>MW</given-names></name><name><surname>Lindstrom</surname><given-names>SE</given-names></name><etal></etal></person-group><year>2009</year><article-title>Emergence of a novel swine-origin influenza A (H1N1) virus in humans.</article-title><source>New England Journal of Medicine</source><volume>361</volume><comment>10.1056/NEJMoa0903810</comment></element-citation></ref><ref id="pcbi.1000564-Sheerar1"><label>45</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sheerar</surname><given-names>M</given-names></name><name><surname>Easterday</surname><given-names>B</given-names></name><name><surname>Hinshaw</surname><given-names>V</given-names></name></person-group><year>1989</year><article-title>Antigenic conservation of H1N1 swine influenza viruses.</article-title><source>Journal of General Virology</source><volume>70</volume><fpage>3297</fpage><lpage>3303</lpage><pub-id pub-id-type="pmid">2558159</pub-id></element-citation></ref><ref id="pcbi.1000564-Vincent1"><label>46</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vincent</surname><given-names>A</given-names></name><name><surname>Lager</surname><given-names>K</given-names></name><name><surname>Ma</surname><given-names>W</given-names></name><name><surname>Lekcharoensuk</surname><given-names>P</given-names></name><name><surname>Gramer</surname><given-names>M</given-names></name><etal></etal></person-group><year>2006</year><article-title>Evaluation of hemagglutinin subtype 1 swine influenza viruses from the United States.</article-title><source>Veterinary Microbiology</source><volume>118</volume><fpage>212</fpage><lpage>222</lpage><pub-id pub-id-type="pmid">16962262</pub-id></element-citation></ref><ref id="pcbi.1000564-Boni1"><label>47</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Boni</surname><given-names>MF</given-names></name><name><surname>Zhou</surname><given-names>Y</given-names></name><name><surname>Taubenberger</surname><given-names>JK</given-names></name><name><surname>Holmes</surname><given-names>EC</given-names></name></person-group><year>2008</year><article-title>Homologous recombination is very rare or absent in human influenza A virus.</article-title><source>Journal of Virology</source><volume>82</volume><fpage>4807</fpage><lpage>4811</lpage><pub-id pub-id-type="pmid">18353939</pub-id></element-citation></ref><ref id="pcbi.1000564-Chare1"><label>48</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chare</surname><given-names>ER</given-names></name><name><surname>Gould</surname><given-names>EA</given-names></name><name><surname>Holmes</surname><given-names>EC</given-names></name></person-group><year>2003</year><article-title>Phylogenetic analysis reveals a low rate of homologous recombination in negative-sense RNA viruses.</article-title><source>J Gen Virol</source><volume>84</volume><fpage>2691</fpage><lpage>2703</lpage><pub-id pub-id-type="pmid">13679603</pub-id></element-citation></ref><ref id="pcbi.1000564-Krasnitz1"><label>49</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Krasnitz</surname><given-names>M</given-names></name><name><surname>Levine</surname><given-names>AJ</given-names></name><name><surname>Rabadan</surname><given-names>R</given-names></name></person-group><year>2008</year><article-title>Anomalies in the influenza virus genome database: New biology or laboratory errors?</article-title><source>Journal of Virology</source><volume>82</volume><fpage>8947</fpage><lpage>8950</lpage><pub-id pub-id-type="pmid">18579605</pub-id></element-citation></ref><ref id="pcbi.1000564-Yang2"><label>50</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Yang</surname><given-names>Z</given-names></name></person-group><year>2006</year><source>Computational molecular evolution</source><publisher-loc>Oxford</publisher-loc><publisher-name>Oxford University Press</publisher-name></element-citation></ref><ref id="pcbi.1000564-Felsenstein1"><label>51</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Felsenstein</surname><given-names>J</given-names></name></person-group><year>1973</year><article-title>Maximum likelihood and minimum steps methods for estimating evolutionary trees from data on discrete characters.</article-title><source>Syst Zool</source><volume>22</volume><fpage>240</fpage><lpage>249</lpage></element-citation></ref><ref id="pcbi.1000564-Felsenstein2"><label>52</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Felsenstein</surname><given-names>J</given-names></name></person-group><year>1981</year><article-title>Evolutionary trees from DNA sequences: a maximum likelihood approach.</article-title><source>J Mol Evol</source><volume>17</volume><fpage>368</fpage><lpage>376</lpage><pub-id pub-id-type="pmid">7288891</pub-id></element-citation></ref><ref id="pcbi.1000564-Bao1"><label>53</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bao</surname><given-names>YM</given-names></name><name><surname>Bolotov</surname><given-names>P</given-names></name><name><surname>Dernovoy</surname><given-names>D</given-names></name><name><surname>Kiryutin</surname><given-names>B</given-names></name><name><surname>Zaslavsky</surname><given-names>L</given-names></name><etal></etal></person-group><year>2008</year><article-title>The influenza virus resource at the national center for biotechnology information.</article-title><source>Journal of Virology</source><volume>82</volume><fpage>596</fpage><lpage>601</lpage><pub-id pub-id-type="pmid">17942553</pub-id></element-citation></ref><ref id="pcbi.1000564-Edgar1"><label>54</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Edgar</surname><given-names>RC</given-names></name></person-group><year>2004</year><article-title>MUSCLE: multiple sequence alignment with high accuracy and high throughput.</article-title><source>Nucl Acids Res</source><volume>32</volume><fpage>1792</fpage><lpage>1797</lpage><pub-id pub-id-type="pmid">15034147</pub-id></element-citation></ref><ref id="pcbi.1000564-Suyama1"><label>55</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Suyama</surname><given-names>M</given-names></name><name><surname>Torrents</surname><given-names>D</given-names></name><name><surname>Bork</surname><given-names>P</given-names></name></person-group><year>2006</year><article-title>PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments.</article-title><source>Nucleic Acids Research</source><volume>34</volume><fpage>W609</fpage><lpage>W612</lpage><pub-id pub-id-type="pmid">16845082</pub-id></element-citation></ref><ref id="pcbi.1000564-Guindon2"><label>56</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Guindon</surname><given-names>S</given-names></name><name><surname>Gascuel</surname><given-names>O</given-names></name></person-group><year>2003</year><article-title>A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood.</article-title><source>Systematic Biology</source><volume>52</volume><fpage>696</fpage><lpage>704</lpage><pub-id pub-id-type="pmid">14530136</pub-id></element-citation></ref><ref id="pcbi.1000564-Hasegawa1"><label>57</label><element-citation publication-type="journal"><person-group person-group-type="author"><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>TA</given-names></name></person-group><year>1985</year><article-title>Dating of the Human Ape Splitting by a Molecular Clock of Mitochondrial-DNA.</article-title><source>Journal of Molecular Evolution</source><volume>22</volume><fpage>160</fpage><lpage>174</lpage><pub-id pub-id-type="pmid">3934395</pub-id></element-citation></ref><ref id="pcbi.1000564-Yang3"><label>58</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname><given-names>Z</given-names></name></person-group><year>1997</year><article-title>PAML: a program package for phylogenetic analysis by maximum likelihood.</article-title><source>Comput Appl Biosci</source><volume>13</volume><fpage>555</fpage><lpage>556</lpage><pub-id pub-id-type="pmid">9367129</pub-id></element-citation></ref><ref id="pcbi.1000564-Yang4"><label>59</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname><given-names>Z</given-names></name></person-group><year>2007</year><article-title>PAML 4: phylogenetic analysis by maximum likelihood.</article-title><source>Molecular Biology and Evolution</source><volume>24</volume><fpage>1586</fpage><lpage>1591</lpage><pub-id pub-id-type="pmid">17483113</pub-id></element-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 000A41 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000A41 | 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. |  |