Long-term adaptation of the influenza A virus by escaping cytotoxic T-cell recognition
Identifieur interne : 000105 ( Pmc/Corpus ); précédent : 000104; suivant : 000106Long-term adaptation of the influenza A virus by escaping cytotoxic T-cell recognition
Auteurs : Rutger G. Woolthuis ; Christiaan H. Van Dorp ; Can Ke Mir ; Rob J. De Boer ; Michiel Van BovenSource :
- Scientific Reports [ 2045-2322 ] ; 2016.
Abstract
The evolutionary adaptation of the influenza A virus (IAV) to human antibodies is well characterised. Much less is known about the long-term evolution of cytotoxic T lymphocyte (CTL) epitopes, which are important antigens for clearance of infection. We construct an antigenic map of IAVs of all human subtypes using a compendium of 142 confirmed CTL epitopes, and show that IAV evolved gradually in the period 1932–2015, with infrequent antigenic jumps in the H3N2 subtype. Intriguingly, the number of CTL epitopes per virus decreases with more than one epitope per three years in the H3N2 subtype (from 84 epitopes per virus in 1968 to 64 in 2015), mostly attributed to the loss of HLA-B epitopes. We confirm these observations with epitope predictions. Our findings indicate that selection pressures imposed by CTL immunity shape the long-term evolution of IAV.
Url:
DOI: 10.1038/srep33334
PubMed: 27629812
PubMed Central: 5024124
Links to Exploration step
PMC:5024124Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Long-term adaptation of the influenza A virus by escaping cytotoxic T-cell recognition</title><author><name sortKey="Woolthuis, Rutger G" sort="Woolthuis, Rutger G" uniqKey="Woolthuis R" first="Rutger G." last="Woolthuis">Rutger G. Woolthuis</name><affiliation><nlm:aff id="a1"><institution>Theoretical Biology, Utrecht University</institution>, Utrecht,<country>The Netherlands</country></nlm:aff></affiliation><affiliation><nlm:aff id="a2"><institution>National Institute for Public Health and the Environment</institution>, Bilthoven,<country>The Netherlands</country></nlm:aff></affiliation></author><author><name sortKey="Van Dorp, Christiaan H" sort="Van Dorp, Christiaan H" uniqKey="Van Dorp C" first="Christiaan H." last="Van Dorp">Christiaan H. Van Dorp</name><affiliation><nlm:aff id="a1"><institution>Theoretical Biology, Utrecht University</institution>, Utrecht,<country>The Netherlands</country></nlm:aff></affiliation><affiliation><nlm:aff id="a2"><institution>National Institute for Public Health and the Environment</institution>, Bilthoven,<country>The Netherlands</country></nlm:aff></affiliation></author><author><name sortKey="Ke Mir, Can" sort="Ke Mir, Can" uniqKey="Ke Mir C" first="Can" last="Ke Mir">Can Ke Mir</name><affiliation><nlm:aff id="a1"><institution>Theoretical Biology, Utrecht University</institution>, Utrecht,<country>The Netherlands</country></nlm:aff></affiliation></author><author><name sortKey="De Boer, Rob J" sort="De Boer, Rob J" uniqKey="De Boer R" first="Rob J." last="De Boer">Rob J. De Boer</name><affiliation><nlm:aff id="a1"><institution>Theoretical Biology, Utrecht University</institution>, Utrecht,<country>The Netherlands</country></nlm:aff></affiliation></author><author><name sortKey="Van Boven, Michiel" sort="Van Boven, Michiel" uniqKey="Van Boven M" first="Michiel" last="Van Boven">Michiel Van Boven</name><affiliation><nlm:aff id="a2"><institution>National Institute for Public Health and the Environment</institution>, Bilthoven,<country>The Netherlands</country></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">27629812</idno><idno type="pmc">5024124</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC5024124</idno><idno type="RBID">PMC:5024124</idno><idno type="doi">10.1038/srep33334</idno><date when="2016">2016</date><idno type="wicri:Area/Pmc/Corpus">000105</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000105</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Long-term adaptation of the influenza A virus by escaping cytotoxic T-cell recognition</title><author><name sortKey="Woolthuis, Rutger G" sort="Woolthuis, Rutger G" uniqKey="Woolthuis R" first="Rutger G." last="Woolthuis">Rutger G. Woolthuis</name><affiliation><nlm:aff id="a1"><institution>Theoretical Biology, Utrecht University</institution>, Utrecht,<country>The Netherlands</country></nlm:aff></affiliation><affiliation><nlm:aff id="a2"><institution>National Institute for Public Health and the Environment</institution>, Bilthoven,<country>The Netherlands</country></nlm:aff></affiliation></author><author><name sortKey="Van Dorp, Christiaan H" sort="Van Dorp, Christiaan H" uniqKey="Van Dorp C" first="Christiaan H." last="Van Dorp">Christiaan H. Van Dorp</name><affiliation><nlm:aff id="a1"><institution>Theoretical Biology, Utrecht University</institution>, Utrecht,<country>The Netherlands</country></nlm:aff></affiliation><affiliation><nlm:aff id="a2"><institution>National Institute for Public Health and the Environment</institution>, Bilthoven,<country>The Netherlands</country></nlm:aff></affiliation></author><author><name sortKey="Ke Mir, Can" sort="Ke Mir, Can" uniqKey="Ke Mir C" first="Can" last="Ke Mir">Can Ke Mir</name><affiliation><nlm:aff id="a1"><institution>Theoretical Biology, Utrecht University</institution>, Utrecht,<country>The Netherlands</country></nlm:aff></affiliation></author><author><name sortKey="De Boer, Rob J" sort="De Boer, Rob J" uniqKey="De Boer R" first="Rob J." last="De Boer">Rob J. De Boer</name><affiliation><nlm:aff id="a1"><institution>Theoretical Biology, Utrecht University</institution>, Utrecht,<country>The Netherlands</country></nlm:aff></affiliation></author><author><name sortKey="Van Boven, Michiel" sort="Van Boven, Michiel" uniqKey="Van Boven M" first="Michiel" last="Van Boven">Michiel Van Boven</name><affiliation><nlm:aff id="a2"><institution>National Institute for Public Health and the Environment</institution>, Bilthoven,<country>The Netherlands</country></nlm:aff></affiliation></author></analytic><series><title level="j">Scientific Reports</title><idno type="eISSN">2045-2322</idno><imprint><date when="2016">2016</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>The evolutionary adaptation of the influenza A virus (IAV) to human antibodies is well characterised. Much less is known about the long-term evolution of cytotoxic T lymphocyte (CTL) epitopes, which are important antigens for clearance of infection. We construct an antigenic map of IAVs of all human subtypes using a compendium of 142 confirmed CTL epitopes, and show that IAV evolved gradually in the period 1932–2015, with infrequent antigenic jumps in the H3N2 subtype. Intriguingly, the number of CTL epitopes per virus decreases with more than one epitope per three years in the H3N2 subtype (from 84 epitopes per virus in 1968 to 64 in 2015), mostly attributed to the loss of HLA-B epitopes. We confirm these observations with epitope predictions. Our findings indicate that selection pressures imposed by CTL immunity shape the long-term evolution of IAV.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Van De Sandt, C E" uniqKey="Van De Sandt C">C. E. van de Sandt</name></author><author><name sortKey="Kreijtz, J H C M" uniqKey="Kreijtz J">J. H. C. M. Kreijtz</name></author><author><name sortKey="Rimmelzwaan, G F" uniqKey="Rimmelzwaan G">G. F. Rimmelzwaan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Westgeest, K B" uniqKey="Westgeest K">K. B. Westgeest</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Luksza, M" uniqKey="Luksza M">M. Łuksza</name></author><author><name sortKey="L Ssig, M" uniqKey="L Ssig M">M. Lässig</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Smith, D J" uniqKey="Smith D">D. J. Smith</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bedford, T" uniqKey="Bedford T">T. Bedford</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fonville, J M" uniqKey="Fonville J">J. M. Fonville</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Thomas, P G" uniqKey="Thomas P">P. G. Thomas</name></author><author><name sortKey="Keating, R" uniqKey="Keating R">R. Keating</name></author><author><name sortKey="Hulse Post, D J" uniqKey="Hulse Post D">D. J. Hulse-Post</name></author><author><name sortKey="Doherty, P C" uniqKey="Doherty P">P. C. Doherty</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="La Gruta, N L" uniqKey="La Gruta N">N. L. La Gruta</name></author><author><name sortKey="Turner, S J" uniqKey="Turner S">S. J. Turner</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lalvani, A" uniqKey="Lalvani A">A. Lalvani</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sallusto, F" uniqKey="Sallusto F">F. Sallusto</name></author><author><name sortKey="Geginat, J" uniqKey="Geginat J">J. Geginat</name></author><author><name sortKey="Lanzavecchia, A" uniqKey="Lanzavecchia A">A. Lanzavecchia</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mcmichael, A J" uniqKey="Mcmichael A">A. J. McMichael</name></author><author><name sortKey="Gotch, F M" uniqKey="Gotch F">F. M. Gotch</name></author><author><name sortKey="Noble, G R" uniqKey="Noble G">G. R. Noble</name></author><author><name sortKey="Beare, P A S" uniqKey="Beare P">P. A. S. Beare</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sridhar, S" uniqKey="Sridhar S">S. Sridhar</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Weinfurter, J T" uniqKey="Weinfurter J">J. T. Weinfurter</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Quiones Parra, S" uniqKey="Quiones Parra S">S. Quiones-Parra</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="La Gruta, N" uniqKey="La Gruta N">N. La Gruta</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Price, G E" uniqKey="Price G">G. E. Price</name></author><author><name sortKey="Ou, R" uniqKey="Ou R">R. Ou</name></author><author><name sortKey="Jiang, H" uniqKey="Jiang H">H. Jiang</name></author><author><name sortKey="Huang, L" uniqKey="Huang L">L. Huang</name></author><author><name sortKey="Moskophidis, D" uniqKey="Moskophidis D">D. Moskophidis</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Berkhoff, E G M" uniqKey="Berkhoff E">E. G. M. Berkhoff</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Berkhoff, E" uniqKey="Berkhoff E">E. Berkhoff</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gong, L I" uniqKey="Gong L">L. I. Gong</name></author><author><name sortKey="Bloom, J D" uniqKey="Bloom J">J. D. Bloom</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Valkenburg, S A" uniqKey="Valkenburg S">S. A. Valkenburg</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kashiwagi, T" uniqKey="Kashiwagi T">T. Kashiwagi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gog, J R" uniqKey="Gog J">J. R. Gog</name></author><author><name sortKey="Rimmelzwaan, G F" uniqKey="Rimmelzwaan G">G. F. Rimmelzwaan</name></author><author><name sortKey="Osterhaus, A D M E" uniqKey="Osterhaus A">A. D. M. E. Osterhaus</name></author><author><name sortKey="Grenfell, B T" uniqKey="Grenfell B">B. T. Grenfell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Machkovech, H M" uniqKey="Machkovech H">H. M. Machkovech</name></author><author><name sortKey="Bedford, T" uniqKey="Bedford T">T. Bedford</name></author><author><name sortKey="Suchard, M A" uniqKey="Suchard M">M. A. Suchard</name></author><author><name sortKey="Bloom, J D" uniqKey="Bloom J">J. D. Bloom</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bui, H H" uniqKey="Bui H">H. H. Bui</name></author><author><name sortKey="Peters, B" uniqKey="Peters B">B. Peters</name></author><author><name sortKey="Assarsson, E" uniqKey="Assarsson E">E. Assarsson</name></author><author><name sortKey="Mbawuike, I" uniqKey="Mbawuike I">I. Mbawuike</name></author><author><name sortKey="Sette, A" uniqKey="Sette A">A. Sette</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Assarsson, E" uniqKey="Assarsson E">E. Assarsson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bogner, P" uniqKey="Bogner P">P. Bogner</name></author><author><name sortKey="Capua, I" uniqKey="Capua I">I. Capua</name></author><author><name sortKey="Lipman, D J" uniqKey="Lipman D">D. J. Lipman</name></author><author><name sortKey="Cox, N J" uniqKey="Cox N">N. J. Cox</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vita, R" uniqKey="Vita R">R. Vita</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rimmelzwaan, G F" uniqKey="Rimmelzwaan G">G. F. Rimmelzwaan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kotturi, M F" uniqKey="Kotturi M">M. F. Kotturi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Webster, R G" uniqKey="Webster R">R. G. Webster</name></author><author><name sortKey="Bean, W J" uniqKey="Bean W">W. J. Bean</name></author><author><name sortKey="Gorman, O T" uniqKey="Gorman O">O. T. Gorman</name></author><author><name sortKey="Chambers, T M" uniqKey="Chambers T">T. M. Chambers</name></author><author><name sortKey="Kawaoka, Y" uniqKey="Kawaoka Y">Y. Kawaoka</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Herfst, S" uniqKey="Herfst S">S. Herfst</name></author><author><name sortKey="Imai, M" uniqKey="Imai M">M. Imai</name></author><author><name sortKey="Kawaoka, Y" uniqKey="Kawaoka Y">Y. Kawaoka</name></author><author><name sortKey="Fouchier, R A M" uniqKey="Fouchier R">R. a. M. Fouchier</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hoof, I" uniqKey="Hoof I">I. Hoof</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wilkinson, T M" uniqKey="Wilkinson T">T. M. Wilkinson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Worobey, M" uniqKey="Worobey M">M. Worobey</name></author><author><name sortKey="Han, G Z" uniqKey="Han G">G.-Z. Han</name></author><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A. Rambaut</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Goldfield, M" uniqKey="Goldfield M">M. Goldfield</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sidney, J" uniqKey="Sidney J">J. Sidney</name></author><author><name sortKey="Peters, B" uniqKey="Peters B">B. Peters</name></author><author><name sortKey="Frahm, N" uniqKey="Frahm N">N. Frahm</name></author><author><name sortKey="Brander, C" uniqKey="Brander C">C. Brander</name></author><author><name sortKey="Sette, A" uniqKey="Sette A">A. Sette</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mukherjee, S" uniqKey="Mukherjee S">S. Mukherjee</name></author><author><name sortKey="Chandra, N" uniqKey="Chandra N">N. Chandra</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sievers, F" uniqKey="Sievers F">F. Sievers</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">Sci Rep</journal-id><journal-id journal-id-type="iso-abbrev">Sci Rep</journal-id><journal-title-group><journal-title>Scientific Reports</journal-title></journal-title-group><issn pub-type="epub">2045-2322</issn><publisher><publisher-name>Nature Publishing Group</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">27629812</article-id><article-id pub-id-type="pmc">5024124</article-id><article-id pub-id-type="pii">srep33334</article-id><article-id pub-id-type="doi">10.1038/srep33334</article-id><article-categories><subj-group subj-group-type="heading"><subject>Article</subject></subj-group></article-categories><title-group><article-title>Long-term adaptation of the influenza A virus by escaping cytotoxic T-cell recognition</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Woolthuis</surname><given-names>Rutger G.</given-names></name><xref ref-type="corresp" rid="c1">a</xref><xref ref-type="aff" rid="a1">1</xref><xref ref-type="aff" rid="a2">2</xref></contrib><contrib contrib-type="author"><name><surname>van Dorp</surname><given-names>Christiaan H.</given-names></name><xref ref-type="aff" rid="a1">1</xref><xref ref-type="aff" rid="a2">2</xref></contrib><contrib contrib-type="author"><name><surname>Keşmir</surname><given-names>Can</given-names></name><xref ref-type="aff" rid="a1">1</xref></contrib><contrib contrib-type="author"><name><surname>de Boer</surname><given-names>Rob J.</given-names></name><xref ref-type="aff" rid="a1">1</xref><xref ref-type="author-notes" rid="n1">*</xref></contrib><contrib contrib-type="author"><name><surname>van Boven</surname><given-names>Michiel</given-names></name><xref ref-type="aff" rid="a2">2</xref><xref ref-type="author-notes" rid="n1">*</xref></contrib><aff id="a1"><label>1</label><institution>Theoretical Biology, Utrecht University</institution>, Utrecht,<country>The Netherlands</country></aff><aff id="a2"><label>2</label><institution>National Institute for Public Health and the Environment</institution>, Bilthoven,<country>The Netherlands</country></aff></contrib-group><author-notes><corresp id="c1"><label>a</label><email>rutger@woolthuis.nl</email></corresp><fn id="n1"><label>*</label><p>These authors contributed equally to this work.</p></fn></author-notes><pub-date pub-type="epub"><day>15</day><month>09</month><year>2016</year></pub-date><pub-date pub-type="collection"><year>2016</year></pub-date><volume>6</volume><elocation-id>33334</elocation-id><history><date date-type="received"><day>01</day><month>06</month><year>2016</year></date><date date-type="accepted"><day>23</day><month>08</month><year>2016</year></date></history><permissions><copyright-statement>Copyright © 2016, The Author(s)</copyright-statement><copyright-year>2016</copyright-year><copyright-holder>The Author(s)</copyright-holder><license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by/4.0/"><pmc-comment>author-paid</pmc-comment>
<license-p>This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link></license-p></license></permissions><abstract><p>The evolutionary adaptation of the influenza A virus (IAV) to human antibodies is well characterised. Much less is known about the long-term evolution of cytotoxic T lymphocyte (CTL) epitopes, which are important antigens for clearance of infection. We construct an antigenic map of IAVs of all human subtypes using a compendium of 142 confirmed CTL epitopes, and show that IAV evolved gradually in the period 1932–2015, with infrequent antigenic jumps in the H3N2 subtype. Intriguingly, the number of CTL epitopes per virus decreases with more than one epitope per three years in the H3N2 subtype (from 84 epitopes per virus in 1968 to 64 in 2015), mostly attributed to the loss of HLA-B epitopes. We confirm these observations with epitope predictions. Our findings indicate that selection pressures imposed by CTL immunity shape the long-term evolution of IAV.</p></abstract></article-meta></front><body><p>IAV persists by continuously escaping pre-existing immunity in the population. Most attention has been on the evolution of surface proteins hemagglutinin (HA) and neuraminidase (NA) that form the main targets of neutralising antibodies<xref ref-type="bibr" rid="b1">1</xref><xref ref-type="bibr" rid="b2">2</xref><xref ref-type="bibr" rid="b3">3</xref>. Antibody mediated immunity is subtype specific and lasts for 2–7 years due to rapid evolution of the antigenic sites on the HA and NA proteins<xref ref-type="bibr" rid="b4">4</xref><xref ref-type="bibr" rid="b5">5</xref><xref ref-type="bibr" rid="b6">6</xref>. IAV also elicits CTL immune responses<xref ref-type="bibr" rid="b7">7</xref><xref ref-type="bibr" rid="b8">8</xref>, which reduce viral spread within the host by killing infected cells. As with memory B cells, memory CTLs mount a fast immune response upon recognition of epitopes years after the primary infection<xref ref-type="bibr" rid="b9">9</xref><xref ref-type="bibr" rid="b10">10</xref>, such that individuals with pre-existing CTLs develop less severe disease<xref ref-type="bibr" rid="b11">11</xref><xref ref-type="bibr" rid="b12">12</xref>. CTLs also provide heterosubtypic immunity<xref ref-type="bibr" rid="b13">13</xref><xref ref-type="bibr" rid="b14">14</xref>, which could be an attractive feature for universal vaccines<xref ref-type="bibr" rid="b15">15</xref>.</p><p>Viruses escape CTL recognition by mutating amino acid residues within CTL epitopes. Such immune escape mutations play an important role in the within-host dynamics of chronic pathogens (e.g. HIV) and are also observed during acute IAV infection<xref ref-type="bibr" rid="b16">16</xref><xref ref-type="bibr" rid="b17">17</xref>. While immune escape mutations in IAV cripple the virus<xref ref-type="bibr" rid="b18">18</xref><xref ref-type="bibr" rid="b19">19</xref>, these mutations can persist in a prolonged infection<xref ref-type="bibr" rid="b20">20</xref>, and at the population level despite the high polymorphism of human leukocyte antigen (HLA)<xref ref-type="bibr" rid="b21">21</xref><xref ref-type="bibr" rid="b22">22</xref>. Recently, positive selection in CTL epitopes has been shown in the nucleoprotein (NP) by comparing human and swine viruses in a phylogenetic analysis<xref ref-type="bibr" rid="b23">23</xref>. Many CTL epitopes have been identified in IAV<xref ref-type="bibr" rid="b24">24</xref><xref ref-type="bibr" rid="b25">25</xref>, but a framework capturing the dynamics of CTL epitopes in all proteins over long evolutionary time is lacking. Here we analyse historical and contemporary IAV sequence data spanning the period 1932–2015, using 142 empirically confirmed CTL epitopes known to date<xref ref-type="bibr" rid="b26">26</xref><xref ref-type="bibr" rid="b27">27</xref> (<xref ref-type="supplementary-material" rid="S1">Supplementary Tables S1, S2 and S3</xref>; Methods).</p><sec disp-level="1"><title>Results</title><sec disp-level="2"><title>Antigenic cartography based on CTL epitopes</title><p>We combine 295 representative human IAVs and the compendium of CTL epitopes into an antigenic map that tracks the long-term evolution of CTL epitopes in IAV across the H1N1, H2N2 and H3N2 subtypes (<xref ref-type="fig" rid="f1">Fig. 1</xref>). Each virus contains a subset of the CTL epitopes (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S1</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S2</xref>), with on average 74 epitopes per virus (summed over all class I HLAs). In total, we find 134 out of the 142 epitopes in these viruses, of which 24 are conserved in the study period (marked in <xref ref-type="supplementary-material" rid="S1">Supplementary Table S2</xref>). At seven loci (positions in the proteome) we find more than one confirmed epitope, i.e. at these loci epitope variants have mutated at some point in time to another epitope variant.</p><p>As a measure of immune similarity we use the Jaccard index, defined as the number of epitopes shared by a pair of viruses divided by their number of unique epitopes (Methods). This measure is attractive biologically as it is based on overlaps of epitopes between viruses, and counts any mutation in an epitope as a CTL escape. Multidimensional scaling (MDS) based on Jaccard distances yields a map in which the distance between any pair of viruses represents the expected number of different epitopes (<xref ref-type="fig" rid="f1">Fig. 1</xref>). The map accurately visualises the expected cross-immunity between viruses, even across subtypes (R<sup>2</sup> = 0.93, <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S3</xref>). We find similar results using Manhattan distances (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S3</xref>). Similar maps based on antibody-mediated immunity can only be constructed for each subtype separately, because hemagglutination inhibition assays are subtype specific<xref ref-type="bibr" rid="b4">4</xref><xref ref-type="bibr" rid="b5">5</xref>.</p><p>In the CTL antigenic map, a time directional evolutionary path runs from the early A/1932 (H1N1) to the recent A/2015 (H3N2), with the H2N2 subtype located between the H1N1 and H3N2 subtypes (<xref ref-type="fig" rid="f1">Fig. 1</xref>). The highly reassorted 2009 pandemic H1N1 (pH1N1) virus is distinct from all other viruses, and is genetically (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S4</xref>) and antigenically (<xref ref-type="fig" rid="f1">Fig. 1</xref>) most closely related to the early twentieth century H1N1 viruses.</p><p>The early H3N2 viruses (1968–1970) are remarkably close to the previously circulating H2N2 viruses (1957–1967; <xref ref-type="fig" rid="f1">Fig. 1</xref>). Conversely, the H3N2 subtype has evolved by more than 35 CTL epitopes over the past 46 years (<xref ref-type="fig" rid="f1">Fig. 1</xref>), with noticeable antigenic jumps in 1993 and 2003–2004. The near-continuous H2N2-H3N2 transition results from the fact that only HA and polymerase basic 1 protein (PB1) were reassorted, while HA carries few known CTL epitopes and PB1 is a highly conserved protein. The antigenic jump in 1993 results from a R384G mutation in NP, affecting four epitopes that are restricted by the HLA alleles B*44, B*08:01, B*27:02, and B*27:05<xref ref-type="bibr" rid="b28">28</xref>. In 2003–2004, five epitopes mutated in the NP, M1, and PA proteins. The CTL antigenic map differs completely from a map based on amino acid sequences of the viruses (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S4</xref>). Here subtypes are clearly separated, mostly due to the variable HA and NA proteins. Excluding HA and NA from the amino acid map decreases the difference with the epitope map, in particular with respect to the subtype clustering (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S4</xref>). Nevertheless, there are still important differences between the antigenic and genetic maps, like the absence of the antigenic jumps in the amino acid map. Overall, IAV gradually drifts away from ancestral viruses by escaping CTL epitopes.</p><p>Not all epitopes are expected to impose equal immune pressures on the IAV. First, epitopes vary in immunodominance due to heterogeneity in HLA-allele frequencies in the population. However, when weighting the epitopes with the corresponding HLA supertype frequencies of five main ethnicities, we find that the distances between viruses are hardly affected (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S5</xref>). Similarly, some epitopes are intrinsically immunodominant by evoking stronger immune responses than others, e.g., due to differences in precursor frequencies<xref ref-type="bibr" rid="b29">29</xref>. Since the relative immunodominance of epitopes is not known, weighting for their exact contribution is impossible. However, when we assign random weights to the relative contribution of the different epitopes, we find that even with substantial immunodominance the results are unaffected (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S6</xref>). Thus, while individual immune-responses might be affected by dominant CTL epitope escapes, the analyses show that at the population level the pattern of viral CTL epitope evolution in <xref ref-type="fig" rid="f1">Fig. 1</xref> is robust.</p><p>The antigenic cartography allows us to investigate avian IAVs that cause human infections and may pose a pandemic threat<xref ref-type="bibr" rid="b14">14</xref><xref ref-type="bibr" rid="b30">30</xref><xref ref-type="bibr" rid="b31">31</xref>. To characterise the antigenic relatedness of avian viruses to human viruses in the context of human CTL immunity, we superimpose representative avian viruses of the subtypes H5N1, H7N9 and H9N2 onto the antigenic map (<xref ref-type="fig" rid="f1">Fig. 1</xref>; Methods; <xref ref-type="supplementary-material" rid="S1">Supplementary Table S4</xref>). The avian viruses are more similar to pH1N1 viruses than to recent H3N2 viruses, suggesting that prior infection with a pH1N1-like virus is likely to yield more CTL immunity to the avian viruses than infection with recent H3N2 viruses. This illustrates that positioning of a virus in the antigenic map is possible on the basis of its sequence, enabling assessment of the extent of pre-existing CTL immunity in the human population.</p></sec><sec disp-level="2"><title>Escape of CTL mediated immunity</title><p>We hypothesise that IAV evolution is in part driven by adaptation to pre-existing CTL immunity. To test this hypothesis, we compile a data set of 62 HLA binding peptides of IAV origin that do not elicit a CTL response (Methods; <xref ref-type="supplementary-material" rid="S1">Supplementary Table S3</xref>). With the exception of a small number of early H1N1 viruses in the period 1932–1935, the confirmed epitopes are significantly more variable than the non-immunogenic HLA-binding peptides (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S7</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S8</xref>). Additionally, the viruses are located considerably closer together in a map obtained using distances based on the non-immunogenic peptides than in the antigenic map (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S9</xref>). The difference in the evolutionary rates between CTL epitopes and non-immunogenic peptides suggests that there is selection for escape of CTL epitopes.</p><p>To substantiate this further, we use 7,347 whole-proteome sequenced viruses (Methods), and find that the total number of epitopes per virus decreases over time (<xref ref-type="fig" rid="f2">Fig. 2a</xref>). In H3N2 approximately one epitope is lost every two-and-a-half years, and in H1N1 approximately one epitope is lost every eight years. Since these time series are ancestrally dependent, we also consider the annual (relative) change in the number of epitopes (Methods), and reject neutral evolution for H3N2 (P = 0.043, Wilcoxon signed-rank test). We cannot distinguish neutral from directed evolution for H1N1 and H2N2, perhaps because the decrease is too small for H1N1 and the time series is too short for H2N2. The number of epitopes decreases most markedly in NP, while remaining almost constant in PB1 (<xref ref-type="fig" rid="f2">Fig. 2b</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S10</xref>).</p><p>The decrease in the number of CTL epitopes could possibly be due to an underrepresentation of recently emerged epitopes which may be less likely to be identified than CTL epitopes in old viruses. However, the number of non-immunogenic peptides is constant over time (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S11</xref>), suggesting that the decrease is unlikely the result of observational bias.</p><p>In addition, we obtain an independent test by computationally predicting potential CTL epitopes. These predictions are free from any observational bias that may be present in the set of empirical epitopes. In the analysis we use prediction tools based on HLA-peptide binding affinities<xref ref-type="bibr" rid="b32">32</xref>, using a suite of binding thresholds and focussing on the top 5% binders. With this threshold 88% of the empirical epitopes is present in the predicted set (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S12</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S14</xref>). For NP in H3N2, which harbours most of the empirical epitopes (<xref ref-type="fig" rid="f2">Fig. 2b</xref>), we find a decrease in the number of strongly binding peptides for the large majority of HLA-B alleles (<xref ref-type="fig" rid="f3">Fig. 3</xref>). We find a similar though non-significant trend when combining the NP, M1, and PB1 proteins in H3N2 (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S14</xref>). These observations are in agreement with a significant decrease in the number of empirical epitopes restricted to HLA-B alleles (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S15</xref>). The decrease is two-fold larger in HLA-B restricted epitopes compared to HLA-A restricted epitopes. The predictions show that the decrease in the number of empirical epitopes per virus (<xref ref-type="fig" rid="f2">Fig. 2a</xref>) is not restricted to the empirical set of epitopes.</p><p>The loss of CTL epitopes is the long-term net result of a continuous and rapid turnover of epitopes (<xref ref-type="fig" rid="f4">Fig. 4</xref>). In the short-term, genetic drift of the virus results in similar (and correlated) losses and gains of epitopes (<xref ref-type="fig" rid="f4">Fig. 4</xref>). Here we calculated the losses and gains of the epitopes weighted by their frequency in the 7,347 viruses. Of the two proteins containing most epitopes, NP has a higher turnover of epitopes than the more conserved PB1, in H1N1 as well as in H3N2. As in the antigenic map we find no unusual losses or gains during the H2N2-H3N2 transition. The turnover of epitopes shows that only a fraction of the escapes results in a lower number of epitopes per virus.</p></sec></sec><sec disp-level="1"><title>Discussion</title><p>By combining CTL epitope data with IAV sequence data, we provide evidence that IAV evolves to escape preexisting cellular immunity in the population. Namely, mutations accumulate more rapidly in epitopes than in non-immunogenic peptides, and the total number of epitopes decreases over time (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S7</xref> and <xref ref-type="fig" rid="f2">Fig. 2</xref>).</p><p>Mutations in epitope sites can result in two different types of escape. First, T-cell receptor (TCR) escapes prevent memory T cells to recognise the mutated epitope. In this case the virus circumvents preexisting memory responses against the epitope in the population, but can trigger the expansion of other immune responses upon infection. The second type are escapes that prevent loading of the epitope. This happens when the mutation either prevents HLA binding or the processing of the epitope. In this case the host cannot elicit any T-cell response against the epitope. Combinations of these escapes are also possible when the mutation lowers the TCR as well as the HLA binding affinity.</p><p>In our analysis we treat both escapes equally, since the exact impact of the mutations is unknown. This might be an oversimplification, because some of the mutations may hardly affect the binding with either TCR or HLA. However, experiments have shown that CTL escape mutations usually affect the immune response<xref ref-type="bibr" rid="b16">16</xref><xref ref-type="bibr" rid="b17">17</xref><xref ref-type="bibr" rid="b18">18</xref><xref ref-type="bibr" rid="b28">28</xref>. In addition, the decrease over time in the number of strongly binding predicted peptides indicates that the overall HLA binding affinity decreases as well, confirming that mutations in epitopes lower the HLA binding. Therefore, we believe that most mutations do affect the CTL response and argue that the simplification of treating any mutation in epitopes as escape is justified.</p><p>The adaptation to preexisting CTL immunity that we observe is in agreement with a recent study showing positive selection at CTL epitope sites in NP by a comparison of evolutionary rates in swine and human IAVs<xref ref-type="bibr" rid="b23">23</xref>. Using a phylogenetic analysis the authors find an increased rate of epitope mutations on the trunk of the tree in human compared to swine IAVs. In line with this, we find that the change in epitopes is larger than the change in non-immunogenic peptides (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S7</xref>). In addition, we find that CTL epitope loci disappear over time, which is a direct indication of adaptive evolution (<xref ref-type="fig" rid="f2">Figs 2</xref> and <xref ref-type="fig" rid="f3">3</xref>). In agreement with Machkovech <italic>et al</italic>.<xref ref-type="bibr" rid="b23">23</xref>, the decrease in the number of CTL epitopes is most pronounced in NP, although we see similar trends in other proteins.</p><p>At this stage the impact of the lower number of CTL epitopes in IAV is unknown, as there are no studies done showing changes in the virulence of influenza. Correlating incidence data with epitope escapes is complex, since other factors, like vaccination and the rapid antigenic drift due to neutralising antibodies<xref ref-type="bibr" rid="b5">5</xref>, interfere with the CTL mediated drift. While additional experiments have to unravel the phenotypic effect of the adaptation, circumventing preexisting CTL immunity could in principle affect the severity of disease, since the lack of preexisting cellular immunity results in more severe infection<xref ref-type="bibr" rid="b11">11</xref><xref ref-type="bibr" rid="b12">12</xref><xref ref-type="bibr" rid="b33">33</xref>.</p><p>We speculate that the decrease in the number of CTL epitopes has been ongoing since the global selective sweep of the internal proteins of IAV in the late 1800s<xref ref-type="bibr" rid="b34">34</xref>. Possibly, the number of CTL epitopes could return to a high value when zoonotic rearrangements introduce proteins that have not undergone evolutionary adaptation to human cellular immunity<xref ref-type="bibr" rid="b23">23</xref>. Examples are the increase of the number of epitopes in PB1 of H3N2 viruses in 1968, and the increase of the number of epitopes found in NP in the pandemic H1N1 viruses in 2009 (<xref ref-type="fig" rid="f2">Fig. 2b</xref>). Underpinning this further, we find a high number of epitopes in the avian viruses, especially in H5N1 (<xref ref-type="fig" rid="f2">Fig. 2a</xref>).</p><p>Escaping humoral immunity releases IAV from antibody-mediated selection pressures generated by viruses circulating in the recent past<xref ref-type="bibr" rid="b4">4</xref><xref ref-type="bibr" rid="b5">5</xref>. This is not the case for cellular immune responses that have the potential to impact the long-term evolution of the virus. Since human IAVs seem to be adapted to CTL immune responses, novel viruses with unadapted reassorted internal genes would be at a selective disadvantage in the human population. This may explain why emerging IAVs with reassorted internal genes historically have only been successful when introduced in viruses that have also reassorted their HA and/or NA genes.</p></sec><sec disp-level="1"><title>Methods</title><sec disp-level="2"><title>Complete virus set</title><p>We collected human IAV sequences (H1N1, H2N2, and H3N2) from the GISAID EpiFlu database from the period 1932–2015 (<ext-link ext-link-type="uri" xlink:href="http://www.gisaid.org">www.gisaid.org</ext-link>; July 2015)<xref ref-type="bibr" rid="b26">26</xref>. We selected viruses for which all protein sequences (except PB1-F2) are available and of which the length of each protein was not shorter than 20 amino acids of the consensus protein length. Moreover, we excluded any virus with sequences containing ambiguous amino acids (X, J, B). We excluded all H1N1 viruses collected in the year 1976, because these viruses originated from an exceptional swine flu outbreak in a military basis in New Jersey<xref ref-type="bibr" rid="b35">35</xref>. Moreover, A/Victoria/36/1988 was excluded, because this virus is a rare reassorted virus, and A/Canada/720/2005 (H2N2), because the subtype H2N2 did not circulate in 2005. This procedure resulted in 3,050 H1N1 (including pH1N1), 65 H2N2, and 4,213 H3N2 viruses.</p></sec><sec disp-level="2"><title>Characteristic virus set</title><p>For computational reasons we compiled a subset with characteristic viruses spanning the period 1932–2015, and used this set for constructing the antigenic map. For the H1N1 and H3N2 subtypes, we included all viruses used by Westgeest <italic>et al</italic>.<xref ref-type="bibr" rid="b2">2</xref> and Bedford <italic>et al</italic>.<xref ref-type="bibr" rid="b5">5</xref>, together with A/Michigan/14/2014 (H3N2), A/Pennsylvania/44/2014 (H3N2), A/Sweden/1/2015 (H3N2) and A/Michigan/02/2015 (H3N2) to cover recent years, and included H2N2 viruses. In total, the characteristic set contains 295 viruses (<xref ref-type="supplementary-material" rid="S1">Supplementary Table S1</xref>), of which 48 are H1N1, 52 are H2N2, 185 are H3N2, and 10 are pH1N1.</p></sec><sec disp-level="2"><title>Epitope compendium</title><p>We downloaded the epitope data from IEDB by using the filters ‘human’ as host organism and ‘Influenza A virus’ as source organism and selecting ‘Linear epitopes’, ‘MHC Class I’ and ‘T Cell Assays’ (<ext-link ext-link-type="uri" xlink:href="http://www.iedb.org">www.iedb.org</ext-link>; September 2015)<xref ref-type="bibr" rid="b27">27</xref>. Next, we gather unique epitopes with all corresponding immune assays. Some epitopes are embedded in longer epitopes belonging to the same HLA supertype<xref ref-type="bibr" rid="b36">36</xref>. In these cases, we kept only the embedded (shorter) epitopes. Thereafter, we removed epitopes with length larger than 12 amino acids, because the exact (shorter) epitope within these longer peptides are not reported. We also excluded epitopes that were not present in any of the IAVs in the complete set described above. Finally, we characterised the candidate epitopes either as immunological epitopes or as non-immunogenic peptides. The candidate epitopes with at least one reported non-negative T-cell assay are added to our epitope compendium, while candidates having negative T-cell responses in all assays are classified as non-immunogenic peptides. The list of 142 epitopes and 62 non-immunogenic peptides is provided in <xref ref-type="supplementary-material" rid="S1">Supplementary Tables S2 and S3</xref> respectively.</p></sec><sec disp-level="2"><title>Jaccard distance</title><p>As a measure for the antigenic distance between two viruses we use a scaled adjusted Jaccard distance, reflecting the expected number of different epitopes between viruses <italic>i</italic> and <italic>j</italic>,</p><p><disp-formula id="eq1"><inline-graphic id="d33e470" xlink:href="srep33334-m1.jpg"></inline-graphic></disp-formula></p><p>Here <italic>S</italic><sub><italic>i</italic></sub> and <italic>S</italic><sub><italic>j</italic></sub> represent the epitopes in viruses <italic>i</italic> and <italic>j</italic>, with <inline-formula id="d33e492"><inline-graphic id="d33e493" xlink:href="srep33334-m2.jpg"></inline-graphic></inline-formula> the number of epitopes that two viruses <italic>i</italic> and <italic>j</italic> have in common. <italic>T</italic><sub><italic>i</italic></sub> and <italic>T</italic><sub><italic>j</italic></sub> represent the loci of viruses <italic>i</italic> and <italic>j</italic> carrying a functional epitope, with <inline-formula id="d33e520"><inline-graphic id="d33e521" xlink:href="srep33334-m3.jpg"></inline-graphic></inline-formula> the number of unique loci in viruses <italic>i</italic> and <italic>j</italic>. Furthermore, <inline-formula id="d33e529"><inline-graphic id="d33e530" xlink:href="srep33334-m4.jpg"></inline-graphic></inline-formula> is the expected number of unique epitope loci, with <italic>N</italic> the total number of viruses. To account for the fact that two different alleles at the same locus is a unit difference in epitopes, we deviate from the classical Jaccard distance by dividing by |<italic>T</italic><sub><italic>i</italic></sub> ∪ <italic>T</italic><sub><italic>j</italic></sub>| instead of |<italic>S</italic><sub><italic>i</italic></sub> ∪ <italic>S</italic><sub><italic>j</italic></sub>|. Notice that |<italic>S</italic><sub><italic>i</italic></sub>| = |<italic>T</italic><sub><italic>i</italic></sub>|, because a virus can at most have one epitope at a given locus, and that the difference between the classical Jaccard distance, e.g., as used elsewhere<xref ref-type="bibr" rid="b37">37</xref>, and our distance is small, since only a few epitopes in our compendium contain multiple alleles.</p></sec><sec disp-level="2"><title>Manhattan distance</title><p>As an alternative to the Jaccard distance, we computed the Manhattan distance by counting the fraction of all possible epitope loci differing between pairs of viruses. In contrast to the Jaccard distance, this measure also counts epitope loci that are absent in both viruses, and therefore depends on the size of the epitope compendium. Additionally, the Jaccard index better counts the number of epitopes that would be recognised when encountering a second virus after infection with the first virus of a pair. Therefore, we believe that the Manhattan distance is less appropriate immunologically. Nevertheless, we found that there is an almost one to one relation between the two distance measures (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S3</xref>), such that the results do not change qualitatively, and to a large extent quantitatively, when using the Manhattan distance.</p></sec><sec disp-level="2"><title>Pairwise proteome distance</title><p>We concatenated the protein sequences for each virus and perform a global alignment using the Clustal Omega package<xref ref-type="bibr" rid="b38">38</xref>. The whole-proteome pair-wise distances of the aligned viruses are calculated with Clustal Omega, which uses the <italic>k</italic>-tuple distance metric (with values between 0 to 1).</p></sec><sec disp-level="2"><title>Antigenic map construction using multi-dimensional scaling</title><p>We follow Smith <italic>et al</italic>.<xref ref-type="bibr" rid="b4">4</xref> by using geometric multi-dimensional scaling (MDS) to represent <italic>N</italic> viruses as elements of a low dimensional Euclidean space (<inline-formula id="d33e606"><inline-graphic id="d33e607" xlink:href="srep33334-m5.jpg"></inline-graphic></inline-formula>), while preserving the distances between viruses as much as possible. Positions <inline-formula id="d33e609"><inline-graphic id="d33e610" xlink:href="srep33334-m6.jpg"></inline-graphic></inline-formula> for the <italic>N</italic> viruses are calculated by minimising the error function <italic>E</italic>.</p><p><disp-formula id="eq7"><inline-graphic id="d33e621" xlink:href="srep33334-m7.jpg"></inline-graphic></disp-formula></p><p>Here <inline-formula id="d33e624"><inline-graphic id="d33e625" xlink:href="srep33334-m8.jpg"></inline-graphic></inline-formula> is the Euclidean norm on the low dimensional space <inline-formula id="d33e627"><inline-graphic id="d33e628" xlink:href="srep33334-m9.jpg"></inline-graphic></inline-formula>. We minimised the error for different initial configurations to ensure that the global minimum is found. With increasing target dimension <italic>n</italic>, the error decreases strongly and the variance in the distances captured (R<sup>2</sup>) increases (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S3</xref>).</p></sec><sec disp-level="2"><title>Superimposing avian viruses</title><p>From GISAID, we collected 21 avian viruses of the H5N1, H7N9, and H9N2 subtypes. To assess the extent of human cellular immunity against avian viruses, we determine their position in the antigenic map given the positions of the human viruses, i.e. we minimise <italic>E</italic> in <xref ref-type="disp-formula" rid="eq7">Eq. 2</xref> for each of the 21 avian viruses independently. See <xref ref-type="supplementary-material" rid="S1">Supplementary Table S4</xref> for an overview of CTL epitopes in the three avian subtypes causing human infections.</p></sec><sec disp-level="2"><title>Predictions</title><p>In the HLA-peptide binding prediction analyses, we use all unique HLA alleles available in the netMHCpan 2.8 prediction tool<xref ref-type="bibr" rid="b32">32</xref>, and measure the relative change in the number of top-binders over a period of time. We focus on the top 5% and the top 10% strongest binders, and use the >50% (weakest binders) as comparative control (% is relative to a large set of peptides used in netMHCpan). Most empirical epitopes belong to the 5% top binders (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S12</xref>).</p></sec><sec disp-level="2"><title>Statistics</title><p>The slopes of the decrease in number of epitopes are calculated by linear regression (<xref ref-type="fig" rid="f2">Fig. 2a</xref>). Since the number of epitopes depends on the number of epitopes in previous years, we also calculated the difference in the number of epitopes from year-to-year (excluding missing years) and used the conservative Wilcoxon signed-rank test (one tailed) to estimate the significance of the decrease. The analysis with non-immunogenic peptides was performed similarly (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S11</xref>).</p><p>Since we expected that the change in the number of binders (<bold>x</bold>) to be alike for HLA alleles with similar binding motifs, we used a statistical model that corrects for these dependencies when determining the significance of the mean decrease per HLA allele. We assume that <inline-formula id="d33e679"><inline-graphic id="d33e680" xlink:href="srep33334-m10.jpg"></inline-graphic></inline-formula>, where <inline-formula id="d33e682"><inline-graphic id="d33e683" xlink:href="srep33334-m11.jpg"></inline-graphic></inline-formula>, and <bold>∑</bold> is the covariance matrix describing the above mentioned dependencies. We calculate <bold>∑</bold> using the similarities between the binding motifs as follows: a HLA allele <italic>i</italic> can be represented as a binary string <italic>b</italic><sup><italic>i</italic></sup> indexed by all 9-mers <italic>p</italic> so that <inline-formula id="d33e704"><inline-graphic id="d33e705" xlink:href="srep33334-m12.jpg"></inline-graphic></inline-formula> if <italic>p</italic> is predicted to bind <italic>i</italic>, and 0 otherwise. We then choose <inline-formula id="d33e714"><inline-graphic id="d33e715" xlink:href="srep33334-m13.jpg"></inline-graphic></inline-formula> for some constant <italic>s</italic> > 0, and obtain maximum-likelihood estimates for the mean decrease <italic>μ</italic> and scaling parameter <italic>s</italic>. This model is selected in favour of a model that assumes independence of the HLA alleles (<italic>P</italic> < 10<sup>−12</sup>, likelihood-ratio test). The maximum likelihood estimates of the mean decrease and corresponding 95% confidence intervals are given in <xref ref-type="fig" rid="f3">Fig. 3</xref>.</p></sec></sec><sec disp-level="1"><title>Additional Information</title><p><bold>How to cite this article</bold>: Woolthuis, R. G. <italic>et al</italic>. Long-term adaptation of the influenza A virus by escaping cytotoxic T-cell recognition. <italic>Sci. Rep.</italic><bold>6</bold>, 33334; doi: 10.1038/srep33334 (2016).</p></sec><sec sec-type="supplementary-material" id="S1"><title>Supplementary Material</title><supplementary-material id="d33e46" content-type="local-data"><caption><title>Supplementary Information</title></caption><media xlink:href="srep33334-s1.pdf"></media></supplementary-material><supplementary-material id="d33e49" content-type="local-data"><caption><title>Supplementary Dataset 1</title></caption><media xlink:href="srep33334-s2.csv"></media></supplementary-material><supplementary-material id="d33e52" content-type="local-data"><caption><title>Supplementary Dataset 2</title></caption><media xlink:href="srep33334-s3.csv"></media></supplementary-material><supplementary-material id="d33e55" content-type="local-data"><caption><title>Supplementary Dataset 3</title></caption><media xlink:href="srep33334-s4.csv"></media></supplementary-material><supplementary-material id="d33e58" content-type="local-data"><caption><title>Supplementary Dataset 4</title></caption><media xlink:href="srep33334-s5.csv"></media></supplementary-material></sec></body><back><ack><p>We thank Josine van Beek, José Ferreira, Marion Hendriks, Scott McDonald, Adam Meijer, and Johannes Textor for discussion and/or technical support. We acknowledge the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO grant 645.000.002 and 823.02.014) and the National Institute for Public Health and the Environment for funding support.</p></ack><ref-list><ref id="b1"><mixed-citation publication-type="journal"><name><surname>van de Sandt</surname><given-names>C. E.</given-names></name>, <name><surname>Kreijtz</surname><given-names>J. H. C. M.</given-names></name> & <name><surname>Rimmelzwaan</surname><given-names>G. F.</given-names></name><article-title>Evasion of influenza A viruses from innate and adaptive immune responses</article-title>. <source>Viruses</source><volume>4</volume>, <fpage>1438</fpage>–<lpage>1476</lpage> (<year>2012</year>).<pub-id pub-id-type="pmid">23170167</pub-id></mixed-citation></ref><ref id="b2"><mixed-citation publication-type="journal"><name><surname>Westgeest</surname><given-names>K. B.</given-names></name><etal></etal>. <article-title>Genetic evolution of the neuraminidase of influenza A (H3N2) viruses from 1968 to 2009 and its correspondence to haemagglutinin evolution</article-title>. <source>Journal of General Virology</source><volume>93</volume>, <fpage>1996</fpage>–<lpage>2007</lpage> (<year>2012</year>).<pub-id pub-id-type="pmid">22718569</pub-id></mixed-citation></ref><ref id="b3"><mixed-citation publication-type="journal"><name><surname>Łuksza</surname><given-names>M.</given-names></name> & <name><surname>Lässig</surname><given-names>M.</given-names></name><article-title>A predictive fitness model for influenza</article-title>. <source>Nature</source><volume>507</volume>, <fpage>57</fpage>–<lpage>61</lpage> (<year>2014</year>).<pub-id pub-id-type="pmid">24572367</pub-id></mixed-citation></ref><ref id="b4"><mixed-citation publication-type="journal"><name><surname>Smith</surname><given-names>D. J.</given-names></name><etal></etal>. <article-title>Mapping the antigenic and genetic evolution of influenza virus</article-title>. <source>Science</source><volume>305</volume>, <fpage>371</fpage>–<lpage>376</lpage> (<year>2004</year>).<pub-id pub-id-type="pmid">15218094</pub-id></mixed-citation></ref><ref id="b5"><mixed-citation publication-type="journal"><name><surname>Bedford</surname><given-names>T.</given-names></name><etal></etal>. <article-title>Integrating influenza antigenic dynamics with molecular evolution</article-title>. <source>eLife</source><volume>3</volume>, <fpage>e01914</fpage> (<year>2014</year>).<pub-id pub-id-type="pmid">24497547</pub-id></mixed-citation></ref><ref id="b6"><mixed-citation publication-type="journal"><name><surname>Fonville</surname><given-names>J. M.</given-names></name><etal></etal>. <article-title>Antibody landscapes after influenza virus infection or vaccination</article-title>. <source>Science</source><volume>346</volume>, <fpage>996</fpage>–<lpage>1000</lpage> (<year>2014</year>).<pub-id pub-id-type="pmid">25414313</pub-id></mixed-citation></ref><ref id="b7"><mixed-citation publication-type="journal"><name><surname>Thomas</surname><given-names>P. G.</given-names></name>, <name><surname>Keating</surname><given-names>R.</given-names></name>, <name><surname>Hulse-Post</surname><given-names>D. J.</given-names></name> & <name><surname>Doherty</surname><given-names>P. C.</given-names></name><article-title>Cell-mediated protection in influenza infection</article-title>. <source>Emerging Infectious Diseases</source><volume>12</volume>, <fpage>48</fpage>–<lpage>54</lpage> (<year>2006</year>).<pub-id pub-id-type="pmid">16494717</pub-id></mixed-citation></ref><ref id="b8"><mixed-citation publication-type="journal"><name><surname>La Gruta</surname><given-names>N. L.</given-names></name> & <name><surname>Turner</surname><given-names>S. J.</given-names></name><article-title>T cell mediated immunity to influenza: mechanisms of viral control</article-title>. <source>Trends in Immunology</source><volume>35</volume>, <fpage>396</fpage>–<lpage>402</lpage> (<year>2014</year>).<pub-id pub-id-type="pmid">25043801</pub-id></mixed-citation></ref><ref id="b9"><mixed-citation publication-type="journal"><name><surname>Lalvani</surname><given-names>A.</given-names></name><etal></etal>. <article-title>Rapid effector function in CD8+ memory T cells</article-title>. <source>The Journal of Experimental Medicine</source><volume>186</volume>, <fpage>859</fpage>–<lpage>865</lpage> (<year>1997</year>).<pub-id pub-id-type="pmid">9294140</pub-id></mixed-citation></ref><ref id="b10"><mixed-citation publication-type="journal"><name><surname>Sallusto</surname><given-names>F.</given-names></name>, <name><surname>Geginat</surname><given-names>J.</given-names></name> & <name><surname>Lanzavecchia</surname><given-names>A.</given-names></name><article-title>Central memory and effector memory T cell subsets: function, generation, and maintenance</article-title>. <source>Annual Review of Immunology</source><volume>22</volume>, <fpage>745</fpage>–<lpage>763</lpage> (<year>2004</year>).</mixed-citation></ref><ref id="b11"><mixed-citation publication-type="journal"><name><surname>McMichael</surname><given-names>A. J.</given-names></name>, <name><surname>Gotch</surname><given-names>F. M.</given-names></name>, <name><surname>Noble</surname><given-names>G. R.</given-names></name> & <name><surname>Beare</surname><given-names>P. A. S.</given-names></name><article-title>Cytotoxic T-cell immunity to influenza</article-title>. <source>New England Journal of Medicine</source><volume>309</volume>, <fpage>13</fpage>–<lpage>17</lpage> (<year>1983</year>).<pub-id pub-id-type="pmid">6602294</pub-id></mixed-citation></ref><ref id="b12"><mixed-citation publication-type="journal"><name><surname>Sridhar</surname><given-names>S.</given-names></name><etal></etal>. <article-title>Cellular immune correlates of protection against symptomatic pandemic influenza</article-title>. <source>Nature Medicine</source><volume>19</volume>, <fpage>1305</fpage>–<lpage>1312</lpage> (<year>2013</year>).</mixed-citation></ref><ref id="b13"><mixed-citation publication-type="journal"><name><surname>Weinfurter</surname><given-names>J. T.</given-names></name><etal></etal>. <article-title>Cross-reactive T cells are involved in rapid clearance of 2009 pandemic H1N1 influenza virus in nonhuman primates</article-title>. <source>PLoS Pathogens</source><volume>7</volume>, <fpage>e1002381</fpage> (<year>2011</year>).<pub-id pub-id-type="pmid">22102819</pub-id></mixed-citation></ref><ref id="b14"><mixed-citation publication-type="other"><name><surname>Quiones-Parra</surname><given-names>S.</given-names></name><etal></etal>. Preexisting CD8+ T-cell immunity to the H7N9 influenza A virus varies across ethnicities. <italic>Proceedings of the National Academy of Sciences</italic> 201322229 (2014).</mixed-citation></ref><ref id="b15"><mixed-citation publication-type="journal"><name><surname>La Gruta</surname><given-names>N.</given-names></name><etal></etal>. <article-title>Role of CD8(+) T-cell immunity in influenza infection: potential use in future vaccine development</article-title>. <source>Expert Review of Respiratory Medicine</source><volume>3</volume>, <fpage>523</fpage>–<lpage>537</lpage> (<year>2009</year>).<pub-id pub-id-type="pmid">20477341</pub-id></mixed-citation></ref><ref id="b16"><mixed-citation publication-type="journal"><name><surname>Price</surname><given-names>G. E.</given-names></name>, <name><surname>Ou</surname><given-names>R.</given-names></name>, <name><surname>Jiang</surname><given-names>H.</given-names></name>, <name><surname>Huang</surname><given-names>L.</given-names></name> & <name><surname>Moskophidis</surname><given-names>D.</given-names></name><article-title>Viral escape by selection of cytotoxic T cell-resistant variants in influenza A virus pneumonia</article-title>. <source>The Journal of Experimental Medicine</source><volume>191</volume>, <fpage>1853</fpage>–<lpage>1867</lpage> (<year>2000</year>).<pub-id pub-id-type="pmid">10839802</pub-id></mixed-citation></ref><ref id="b17"><mixed-citation publication-type="journal"><name><surname>Berkhoff</surname><given-names>E. G. M.</given-names></name><etal></etal>. <article-title>A mutation in the HLA-B*2705-restricted NP383-391 epitope affects the human influenza A virus-specific cytotoxic t-lymphocyte response <italic>in vitro</italic></article-title>. <source>Journal of Virology</source><volume>78</volume>, <fpage>5216</fpage>–<lpage>5222</lpage> (<year>2004</year>).<pub-id pub-id-type="pmid">15113903</pub-id></mixed-citation></ref><ref id="b18"><mixed-citation publication-type="journal"><name><surname>Berkhoff</surname><given-names>E.</given-names></name><etal></etal>. <article-title>Fitness costs limit escape from cytotoxic T lymphocytes by influenza A viruses</article-title>. <source>Vaccine</source><volume>24</volume>, <fpage>6594</fpage>–<lpage>6596</lpage> (<year>2006</year>).<pub-id pub-id-type="pmid">16837112</pub-id></mixed-citation></ref><ref id="b19"><mixed-citation publication-type="journal"><name><surname>Gong</surname><given-names>L. I.</given-names></name> & <name><surname>Bloom</surname><given-names>J. D.</given-names></name><article-title>Epistatically interacting substitutions are enriched during adaptive protein evolution</article-title>. <source>PLoS Genet</source><volume>10</volume>, <fpage>e1004328</fpage> (<year>2014</year>).<pub-id pub-id-type="pmid">24811236</pub-id></mixed-citation></ref><ref id="b20"><mixed-citation publication-type="journal"><name><surname>Valkenburg</surname><given-names>S. A.</given-names></name><etal></etal>. <article-title>Acute emergence and reversion of influenza A virus quasispecies within CD8+ T cell antigenic peptides</article-title>. <source>Nature Communications</source><volume>4</volume>, <fpage>2663</fpage> (<year>2013</year>).</mixed-citation></ref><ref id="b21"><mixed-citation publication-type="journal"><name><surname>Kashiwagi</surname><given-names>T.</given-names></name><etal></etal>. <article-title>Emergence of new influenza A viruses which carry an escape mutation of the HLA-B27-restricted CTL epitope of NP in japan</article-title>. <source>Microbiology and Immunology</source><volume>44</volume>, <fpage>867</fpage>–<lpage>870</lpage> (<year>2000</year>).<pub-id pub-id-type="pmid">11128072</pub-id></mixed-citation></ref><ref id="b22"><mixed-citation publication-type="journal"><name><surname>Gog</surname><given-names>J. R.</given-names></name>, <name><surname>Rimmelzwaan</surname><given-names>G. F.</given-names></name>, <name><surname>Osterhaus</surname><given-names>A. D. M. E.</given-names></name> & <name><surname>Grenfell</surname><given-names>B. T.</given-names></name><article-title>Population dynamics of rapid fixation in cytotoxic T lymphocyte escape mutants of influenza A</article-title>. <source>Proceedings of the National Academy of Sciences</source><volume>100</volume>, <fpage>11143</fpage>–<lpage>11147</lpage> (<year>2003</year>).</mixed-citation></ref><ref id="b23"><mixed-citation publication-type="journal"><name><surname>Machkovech</surname><given-names>H. M.</given-names></name>, <name><surname>Bedford</surname><given-names>T.</given-names></name>, <name><surname>Suchard</surname><given-names>M. A.</given-names></name> & <name><surname>Bloom</surname><given-names>J. D.</given-names></name><article-title>Positive selection in CD8+ T-cell epitopes of influenza nucleoprotein revealed by a comparative analysis of human and swine viral lineages</article-title>. <source>Journal of Virology</source><volume>89</volume>, <fpage>01571</fpage>–<lpage>15</lpage> (<year>2015</year>).</mixed-citation></ref><ref id="b24"><mixed-citation publication-type="journal"><name><surname>Bui</surname><given-names>H. H.</given-names></name>, <name><surname>Peters</surname><given-names>B.</given-names></name>, <name><surname>Assarsson</surname><given-names>E.</given-names></name>, <name><surname>Mbawuike</surname><given-names>I.</given-names></name> & <name><surname>Sette</surname><given-names>A.</given-names></name><article-title>Ab and T cell epitopes of influenza A virus, knowledge and opportunities</article-title>. <source>Proceedings of the National Academy of Sciences</source><volume>104</volume>, <fpage>246</fpage>–<lpage>251</lpage> (<year>2007</year>).</mixed-citation></ref><ref id="b25"><mixed-citation publication-type="journal"><name><surname>Assarsson</surname><given-names>E.</given-names></name><etal></etal>. <article-title>Immunomic analysis of the repertoire of T-cell specificities for influenza A virus in humans</article-title>. <source>Journal of Virology</source><volume>82</volume>, <fpage>12241</fpage>–<lpage>12251</lpage> (<year>2008</year>).<pub-id pub-id-type="pmid">18842709</pub-id></mixed-citation></ref><ref id="b26"><mixed-citation publication-type="journal"><name><surname>Bogner</surname><given-names>P.</given-names></name>, <name><surname>Capua</surname><given-names>I.</given-names></name>, <name><surname>Lipman</surname><given-names>D. J.</given-names></name>, <name><surname>Cox</surname><given-names>N. J.</given-names></name><etal></etal>. <article-title>A global initiative on sharing avian flu data</article-title>. <source>Nature</source><volume>442</volume>, <fpage>981</fpage>–<lpage>981</lpage> (<year>2006</year>).</mixed-citation></ref><ref id="b27"><mixed-citation publication-type="journal"><name><surname>Vita</surname><given-names>R.</given-names></name><etal></etal>. <article-title>The immune epitope database 2.0</article-title>. <source>Nucleic Acids Research</source><volume>38</volume>, <fpage>D854</fpage>–<lpage>862</lpage> (<year>2010</year>).<pub-id pub-id-type="pmid">19906713</pub-id></mixed-citation></ref><ref id="b28"><mixed-citation publication-type="journal"><name><surname>Rimmelzwaan</surname><given-names>G. F.</given-names></name><etal></etal>. <article-title>Sequence variation in the influenza A virus nucleoprotein associated with escape from cytotoxic T lymphocytes</article-title>. <source>Virus Research</source><volume>103</volume>, <fpage>97</fpage>–<lpage>100</lpage> (<year>2004</year>).<pub-id pub-id-type="pmid">15163496</pub-id></mixed-citation></ref><ref id="b29"><mixed-citation publication-type="journal"><name><surname>Kotturi</surname><given-names>M. F.</given-names></name><etal></etal>. <article-title>Naive precursor frequencies and MHC binding rather than the degree of epitope diversity shape CD8+ T cell immunodominance</article-title>. <source>Journal of Immunology (Baltimore, Md.: 1950)</source><volume>181</volume>, <fpage>2124</fpage>–<lpage>2133</lpage> (<year>2008</year>).</mixed-citation></ref><ref id="b30"><mixed-citation publication-type="journal"><name><surname>Webster</surname><given-names>R. G.</given-names></name>, <name><surname>Bean</surname><given-names>W. J.</given-names></name>, <name><surname>Gorman</surname><given-names>O. T.</given-names></name>, <name><surname>Chambers</surname><given-names>T. M.</given-names></name> & <name><surname>Kawaoka</surname><given-names>Y.</given-names></name><article-title>Evolution and ecology of influenza A viruses</article-title>. <source>Microbiological Reviews</source><volume>56</volume>, <fpage>152</fpage>–<lpage>179</lpage> (<year>1992</year>).<pub-id pub-id-type="pmid">1579108</pub-id></mixed-citation></ref><ref id="b31"><mixed-citation publication-type="journal"><name><surname>Herfst</surname><given-names>S.</given-names></name>, <name><surname>Imai</surname><given-names>M.</given-names></name>, <name><surname>Kawaoka</surname><given-names>Y.</given-names></name> & <name><surname>Fouchier</surname><given-names>R. a. M.</given-names></name><article-title>Avian influenza virus transmission to mammals</article-title>. <source>Current Topics in Microbiology and Immunology</source><volume>385</volume>, <fpage>137</fpage>–<lpage>155</lpage> (<year>2014</year>).<pub-id pub-id-type="pmid">25048542</pub-id></mixed-citation></ref><ref id="b32"><mixed-citation publication-type="journal"><name><surname>Hoof</surname><given-names>I.</given-names></name><etal></etal>. <article-title>NetMHCpan, a method for MHC class I binding prediction beyond humans</article-title>. <source>Immunogenetics</source><volume>61</volume>, <fpage>1</fpage>–<lpage>13</lpage> (<year>2009</year>).<pub-id pub-id-type="pmid">19002680</pub-id></mixed-citation></ref><ref id="b33"><mixed-citation publication-type="journal"><name><surname>Wilkinson</surname><given-names>T. M.</given-names></name><etal></etal>. <article-title>Preexisting influenza-specific CD4+ T cells correlate with disease protection against influenza challenge in humans</article-title>. <source>Nature Medicine</source><volume>18</volume>, <fpage>274</fpage>–<lpage>280</lpage> (<year>2012</year>).</mixed-citation></ref><ref id="b34"><mixed-citation publication-type="journal"><name><surname>Worobey</surname><given-names>M.</given-names></name>, <name><surname>Han</surname><given-names>G.-Z.</given-names></name> & <name><surname>Rambaut</surname><given-names>A.</given-names></name><article-title>A synchronized global sweep of the internal genes of modern avian influenza virus</article-title>. <source>Nature</source><volume>508</volume>, <fpage>254</fpage>–<lpage>257</lpage> (<year>2014</year>).<pub-id pub-id-type="pmid">24531761</pub-id></mixed-citation></ref><ref id="b35"><mixed-citation publication-type="journal"><name><surname>Goldfield</surname><given-names>M.</given-names></name><etal></etal>. <article-title>Influenza in New Jersey in 1976: isolations of influenza A/New Jersey/76 virus at fort dix</article-title>. <source>The Journal of Infectious Diseases</source><volume>136</volume> Suppl, <fpage>S347</fpage>–<lpage>355</lpage> (<year>1977</year>).<pub-id pub-id-type="pmid">606758</pub-id></mixed-citation></ref><ref id="b36"><mixed-citation publication-type="journal"><name><surname>Sidney</surname><given-names>J.</given-names></name>, <name><surname>Peters</surname><given-names>B.</given-names></name>, <name><surname>Frahm</surname><given-names>N.</given-names></name>, <name><surname>Brander</surname><given-names>C.</given-names></name> & <name><surname>Sette</surname><given-names>A.</given-names></name><article-title>HLA class I supertypes: a revised and updated classification</article-title>. <source>BMC Immunology</source><volume>9</volume>, <fpage>1</fpage> (<year>2008</year>).<pub-id pub-id-type="pmid">18211710</pub-id></mixed-citation></ref><ref id="b37"><mixed-citation publication-type="journal"><name><surname>Mukherjee</surname><given-names>S.</given-names></name> & <name><surname>Chandra</surname><given-names>N.</given-names></name><article-title>Grouping of large populations into few CTL immune ‘response-types’ from influenza H1N1 genome analysis</article-title>. <source>Clinical & Translational Immunology</source><volume>3</volume>, <fpage>e24</fpage> (<year>2014</year>).<pub-id pub-id-type="pmid">25505972</pub-id></mixed-citation></ref><ref id="b38"><mixed-citation publication-type="journal"><name><surname>Sievers</surname><given-names>F.</given-names></name><etal></etal>. <article-title>Fast, scalable generation of high-quality protein multiple sequence alignments using clustal omega</article-title>. <source>Molecular Systems Biology</source><volume>7</volume>, <fpage>539</fpage> (<year>2011</year>).<pub-id pub-id-type="pmid">21988835</pub-id></mixed-citation></ref></ref-list><fn-group><fn><p><bold>Author Contributions</bold> R.G.W., C.H.v.D., C.K., R.J.d.B. and M.v.B. designed the study. R.G.W. and C.H.v.D. collected data and performed the analyses. R.G.W., C.H.v.D., C.K., R.J.d.B. and M.v.B. discussed the results and helped writing the manuscript.</p></fn></fn-group></back><floats-group><fig id="f1"><label>Figure 1</label><caption><title>CTL epitope evolution in the influenza A virus.</title><p>Antigenic map of 295 representative influenza A viruses spanning the period 1932–2015 (□H1N1, △H2N2, ○H3N2, <inline-formula id="d33e765"><inline-graphic id="d33e766" xlink:href="srep33334-m15.jpg"></inline-graphic></inline-formula> pH1N1) based on 134 CTL epitopes. The H3N2 subtype has evolved extensively over the period 1968–2015, while the H2N2 and H3N2 viruses circulating in the late 1960s are antigenically close. Recent avian viruses are superimposed independently onto the antigenic map (<inline-formula id="d33e768"><inline-graphic id="d33e769" xlink:href="srep33334-m16.jpg"></inline-graphic></inline-formula> H5N1, <inline-formula id="d33e771"><inline-graphic id="d33e772" xlink:href="srep33334-m17.jpg"></inline-graphic></inline-formula> H7N9 and <inline-formula id="d33e774"><inline-graphic id="d33e775" xlink:href="srep33334-m18.jpg"></inline-graphic></inline-formula> H9N2), using the 134 CTL epitopes of human IAV origin. The map is constructed using multi-dimensional scaling (MDS) based on Jaccard distances, explaining 93% of the antigenic distances (Methods). Scale bar denotes expected differences in the number of epitopes; colours indicate collection year of the virus.</p></caption><graphic xlink:href="srep33334-f1"></graphic></fig><fig id="f2"><label>Figure 2</label><caption><title>Decrease in the number of CTL epitopes over time.</title><p>(<bold>a</bold>) The number of CTL epitopes per virus decreases with 0.14 ± 0.01, 0.20 ± 0.06, and 0.40 ± 0.02, epitopes per year (SEM) in H1N1, H2N2, and H3N2, respectively (Methods). The number of CTL epitopes in avian viruses is larger than in recent human viruses (<inline-formula id="d33e786"><inline-graphic id="d33e787" xlink:href="srep33334-m19.jpg"></inline-graphic></inline-formula> H5N1, <inline-formula id="d33e789"><inline-graphic id="d33e790" xlink:href="srep33334-m20.jpg"></inline-graphic></inline-formula> H7N9 and <inline-formula id="d33e792"><inline-graphic id="d33e793" xlink:href="srep33334-m21.jpg"></inline-graphic></inline-formula> H9N2). (<bold>b</bold>) Protein specific CTL epitope dynamics showing a decrease in the number of epitopes in NP, matrix 1 protein (M1) and non-structural protein (NS1). The other IAV proteins are laid down in <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S10</xref>. Reassortment events resulting in the appearance of novel proteins are marked by triangles. Subtype replacements are marked by vertical lines. The analysis is based on 7,347 whole-proteome sequenced viruses (Methods).</p></caption><graphic xlink:href="srep33334-f2"></graphic></fig><fig id="f3"><label>Figure 3</label><caption><title>Decrease in the number of predicted epitopes.</title><p>The change in number of top 5% strongest binding peptides in NP of H3N2 viruses for HLA-A (<bold>a</bold>) and HLA-B (<bold>b</bold>) alleles is shown over the period 1968–2015 (histograms; Methods). These peptides decrease in number for the large majority of HLA-B alleles (P = 0.006; mean change <inline-formula id="d33e814"><inline-graphic id="d33e815" xlink:href="srep33334-m14.jpg"></inline-graphic></inline-formula> (95% CI: −15.6%, −2.6%) over the study period, Methods), but not for HLA-A alleles (P = 0.46). Despite the high level of heterogeneity in binding motifs, we find a decrease in the number of top-binders for a broad range of HLA-B alleles (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S12</xref>). The position of the twelve HLA supertype representatives are indicated with arrows. Fitted Gaussian curves indicate the decrease calculated for the thresholds ≤5%, ≤10% (top binders), and ≥50% (weak binders).</p></caption><graphic xlink:href="srep33334-f3"></graphic></fig><fig id="f4"><label>Figure 4</label><caption><title>CTL epitope dynamics.</title><p>Epitopes appear and disappear continuously in H1N1 (<bold>a</bold>) and H2N2-H3N2 (<bold>b</bold>), with the highest turnover in NP. Within years epitope gains and losses are correlated, due to annual variation in genetic drift. We calculate for each epitope the difference in the fraction of viruses harbouring the epitope in a given year and the year before. Summing over all positive (negative) differences of the fractions gives the total gain (loss).</p></caption><graphic xlink:href="srep33334-f4"></graphic></fig></floats-group></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 000105 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000105 | 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é= PMC:5024124
|texte= Long-term adaptation of the influenza A virus by escaping cytotoxic T-cell recognition
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:27629812" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a H2N2V1
| This area was generated with Dilib version V0.6.33. |  |