# Phylogeography and Antigenic Diversity of Low-Pathogenic Avian Influenza H13 and H16 Viruses

Josanne H Verhagen, Marjolein Poen, David E Stallknecht, Stefan van der Vliet, Pascal Lexmond, Srinand Sreevatsan, Rebecca L Poulson, Ron A. M. Fouchier, Camille Lebarbenchon

► To cite this version:

Josanne H Verhagen, Marjolein Poen, David E Stallknecht, Stefan van der Vliet, Pascal Lexmond, et al.. Phylogeography and Antigenic Diversity of Low-Pathogenic Avian Influenza H13 and H16 Viruses. Journal of Virology, 2020, 94, 10.1128/jvi.00537-20. hal-04109937

HAL Id: hal-04109937

https://hal.univ-reunion.fr/hal-04109937

Submitted on 30 May 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Phylogeography and Antigenic Diversity of Low-Pathogenic Avian Influenza H13 and H16 Viruses

Josanne H. Verhagen,ᵃ,ᵇ Marjolein Poen,ª David E. Stallknecht,ᶜ Stefan van der Vliet,ª Pascal Lexmond,ª Srinand Sreevatsan,ᵈ Rebecca L. Poulson,ᶜ Ron A. M. Fouchier,ª Camille Lebarbenchon,ᶜ,ᵉ

**ABSTRACT** Low-pathogenic avian influenza viruses (LPAIVs) are genetically highly variable and have diversified into multiple evolutionary lineages that are primarily associated with wild-bird reservoirs. Antigenic variation has been described for mammalian influenza viruses and for highly pathogenic avian influenza viruses that circulate in poultry, but much less is known about antigenic variation of LPAIVs. In this study, we focused on H13 and H16 LPAIVs that circulate globally in gulls. We investigated the evolutionary history and intercontinental gene flow based on the hemagglutinin (HA) gene and used representative viruses from genetically distinct lineages to determine their antigenic properties by hemagglutination inhibition assays. For H13, at least three distinct genetic clades were evident, while for H16, at least two distinct genetic clades were evident. Twenty and ten events of intercontinental gene flow were identified for H13 and H16 viruses, respectively. At least two antigenic variants of H13 and at least one antigenic variant of H16 were identified. Amino acid positions in the HA protein that may be involved in the antigenic variation were inferred, and some of the positions were located near the receptor binding site of the HA protein, as they are in the HA protein of mammalian influenza A viruses. These findings suggest independent circulation of H13 and H16 subtypes in gull populations, as antigenic patterns do not overlap, and they contribute to the understanding of the genetic and antigenic variation of LPAIVs naturally circulating in wild birds.

**IMPORTANCE** Wild birds play a major role in the epidemiology of low-pathogenic avian influenza viruses (LPAIVs), which are occasionally transmitted—directly or indirectly—from them to other species, including domestic animals, wild mammals, and humans, where they can cause subclinical to fatal disease. Despite a multitude of genetic studies, the antigenic variation of LPAIVs in wild birds is poorly understood. Here, we investigated the evolutionary history, intercontinental gene flow, and antigenic variation among H13 and H16 LPAIVs. The circulation of subtypes H13 and H16 seems to be maintained by a narrower host range, in particular gulls, than the majority of LPAIV subtypes and may therefore serve as a model for evolution and epidemiology of H1 to H12 LPAIVs in wild birds. The findings suggest that H13 and H16 LPAIVs circulate independently of each other and emphasize the need to investigate within-clade antigenic variation of LPAIVs in wild birds.

**KEYWORDS** avian viruses, influenza, evolution, epidemiology, ecology, antigenic variation, seabird, gulls, wild birds

**Citation** Verhagen JH, Poen M, Stallknecht DE, van der Vliet S, Lexmond P, Sreevatsan S, Poulson RL, Fouchier RAM, Lebarbenchon C. 2020. Phylogeography and antigenic diversity of low-pathogenic avian influenza H13 and H16 viruses. J Virol 94:e00537-20. https://doi.org/10.1128/JVI.00537-20.

**Editor** Colin R. Parrish, Cornell University

**Copyright** © 2020 Verhagen et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Address correspondence to Josanne H. Verhagen, josanne.verhagen@lnu.se.

**Received** 26 March 2020

**Accepted** 13 April 2020

**Accepted manuscript posted online** 22 April 2020

**Published** 16 June 2020

<footnote>ªErasmus Medical Center, Department of Viroscience, Rotterdam, The Netherlands</footnote>

<footnote>ᵇLinnaeus University, Department of Biology and Environmental Science, Kalmar, Sweden</footnote>

<footnote>ᶜSoutheastern Cooperative Wildlife Disease Study, College of Veterinary Medicine, Department of Population Health, University of Georgia, Athens, Georgia, USA</footnote>

<footnote>ᵈMichigan State University, College of Veterinary Medicine, Department of Pathobiology and Diagnostic Investigation, East Lansing, Michigan, USA</footnote>

<footnote>ᵉUniversité de La Réunion, UMR Processus Infectieux en Milieu Insulaire Tropical, INSERM 1187, CNRS 9192, IRD 249, Sainte-Clotilde, La Réunion, France</footnote>
Wild birds of the orders Anseriformes (mainly ducks, geese, and swans) and Charadriiformes (mainly gulls, terns, and waders) play a major role in the epidemiology of low-pathogenic avian influenza viruses (LPAIVs). LPAIVs are characterized into subtypes based on the surface proteins hemagglutinin (HA; H1 to H16) and neuraminidase (NA; N1 to N9), e.g., H13N6. Ducks play an important role in the epidemiology of most LPAIV subtypes. However, birds of the order Charadriiformes—in particular gulls—are the major reservoir for subtypes H13 and H16 (see Table S1 in the supplemental material) (1–4). High prevalence of H13 and/or H16 LPAIVs has been observed in juvenile gulls at breeding colony sites (5–7) and in adults during spring and/or fall migration (8, 9). H13 and H16 viruses have a global distribution. Since their first detection in 1977, H13 viruses have been detected in North America, South America, Europe, Asia, Africa, and Oceania. Since their first detection in 1975, H16 viruses have been detected in North America, South America, Europe, and Asia. The spatial isolation of host populations has shaped LPAIV evolution and led to the independent circulation of different virus gene pools between the Western and Eastern hemispheres (10). However, some pelagic gull populations connect multiple continents through seasonal migration and overlapping distributions and could facilitate rapid and long-distance dispersal of LPAIV genomes (2, 9, 11–14). For instance, great black-backed gulls (*Larus marinus*) migrate between Europe and the east coast of North America, and LPAIVs consisting of both North American and Eurasian genes have been isolated from this species (12). Upon intercontinental gene flow, i.e., the movement of genes between the different continents, some LPAIV genes seem to have become established in the population, e.g., H6 (15).

Influenza A viruses (IAVs) belong to the family *Orthomyxoviridae* and are negative-sense single-stranded RNA viruses with a segmented genome. The genome consists of eight segments encoding 12 proteins or more, including the surface proteins HA and NA. The HA protein of IAV is a major determinant for virus binding to cells and subsequent cell entry and for generation of IAV-specific antibodies, and it is thus subjected to strong selective pressure (16). Indeed, in wild birds—in particular mallards (*Anas platyrhynchos*)—LPAIV infection dynamics between LPAIV subtypes seem to be shaped partially by preexisting homo- or heterologous antibodies (17). Furthermore, within other host systems, evasion of IAV-specific antibodies by IAVs—called antigenic variation—has been described for seasonal human IAVs (18, 19), swine IAVs (20–22), and equine IAVs (23) and for highly pathogenic avian influenza viruses (HPAIVs) that circulate in poultry (24, 25). Despite numerous studies on the genetic variation of LPAIVs in wild birds, the antigenic variation within LPAIV subtypes that circulate in wild birds has barely been investigated (26, 27).

To better understand LPAIV epidemiology in gulls, we investigated the global distribution of H13 and H16 LPAIVs and the antigenic variation of a representative subset of H13 and H16 LPAIVs. Based on the sequencing of HA genes of 84 viruses and hemagglutination inhibition assays, we showed that intercontinental H13 and H16 gene flow occurred frequently and that H16 genetic lineages did not form antigenic clusters, suggesting that clade-defining mutations were not in critical epitopes (i.e., part of the antigen that binds to specific antibodies). In contrast, the H13 genetic clades partially corresponded with the antigenic variation of H13 LPAIVs, suggesting that some of the clade-defining mutations were in critical epitopes.

## RESULTS

**Phylogeographic structure and intercontinental gene flow.** Phylogenetic analyses supported the idea that the H13 HA was structured in three major genetic lineages (lineages A to C) (Fig. 1; also, see Fig. S1 and S2). The time to the most recent common ancestor (tMRCA) of the H13 HA gene was dated to 1927 (±95% highest posterior density [HPD], 1920–1934). The tMRCA of viruses of clade A (1963 [1958–1966]) was longer than those of clade B (1975 [1974–1976]) and C (1977 [1976–1978]). Our analyses support the idea that the geographic origins of H13 viruses of clades B and C could be North America and Europe, respectively (posterior probabilities for the
<figure><img src="image_3.png" /><figcaption>**FIG 1** Maximum-clade-credibility (MCC) tree for influenza A virus H13 hemagglutinin subtype (n = 338). Branches were colored according to most probable geographic origin (red, North America; orange, South America; dark blue, Europe; light blue, Asia; green, Oceania; gray, not identified). Black node bars represent the 95% highest posterior densities for times to the common ancestors. Numbers highlight intercontinental gene flow events as detailed in Table 1 and Fig. S3. Virus strain names and posterior probabilities are detailed in Fig. S2.</figcaption></figure>

geographic origin of the MRCA, 1 for clade B and 1 for clade C). For clade A, limited historical data on viruses from different locations as well as low posterior probability (0.62) precludes a conclusion on the geographic origin of the MRCA.

Since the first isolation of an H13 IAV from a gull in 1977, 20 potential events of intercontinental gene flow were identified (indicated by the numerals 1 to 20 in Fig. 1 and Table 1 and also in Fig. S3). Clade A supports the maintenance of H13 in European gulls, with evidence of multiple introductions to North America and Asia (events 3, 5, 6, 9, and 10), and a reverse introduction from North America to Asia (event 8). Clade C was also composed mainly of viruses circulating in Europe, with evidence of multiple introductions to North America (events 12, 15, and 19) and Asia (events 13, 16, and 17).
<table><caption>TABLE 1 Intercontinental gene flow events for influenza A virus H13 hemagglutinin<sup>a</sup></caption><thead><tr><th>H13 clade</th><th>Event</th><th>Time of MRCA ± 95% HPD</th><th>Geographic origin of the MRCA (posterior density)</th><th>Location of introduction</th></tr></thead><tbody><tr><td rowspan="10">A</td><td>1</td><td>1963 (1958–1966)</td><td>North America (0.62)</td><td>Oceania</td></tr><tr><td>2</td><td>1974 (1972–1975)</td><td>North America (0.73)</td><td>Europe</td></tr><tr><td>3</td><td>1988 (1987–1989)</td><td>Europe (1)</td><td>North America</td></tr><tr><td>4</td><td>1990 (1988–1991)</td><td>Europe (0.82)</td><td>South America</td></tr><tr><td>5</td><td>1996 (1995–1997)</td><td>Europe (0.75)</td><td>Asia</td></tr><tr><td>6</td><td>2003 (2003–2004)</td><td>Europe (1)</td><td>Asia</td></tr><tr><td>7</td><td>2005 (2004–2005)</td><td>Asia (0.48)</td><td>North America</td></tr><tr><td>8</td><td>2009 (2009–2010)</td><td>North America (0.9)</td><td>Asia</td></tr><tr><td>9</td><td>2006 (2006–2007)</td><td>Europe (0.96)</td><td>Asia</td></tr><tr><td>10</td><td>2011 (2010–2011)</td><td>Europe (1)</td><td>Asia</td></tr><tr><td>B</td><td>11</td><td>2013 (2012–2014)</td><td>North America (0.96)</td><td>South America</td></tr><tr><td rowspan="10">C</td><td>12</td><td>1987 (1985–1988)</td><td>Europe (0.99)</td><td>North America</td></tr><tr><td>13</td><td>2002 (2002–2003)</td><td>Europe (1)</td><td>Asia</td></tr><tr><td>14</td><td>2005 (2004–2005)</td><td>Asia (0.55)</td><td>North America</td></tr><tr><td>15</td><td>2010 (2009–2010)</td><td>Europe (1)</td><td>North America</td></tr><tr><td>16</td><td>2004 (2003–2005)</td><td>Europe (0.97)</td><td>Asia</td></tr><tr><td>17</td><td>2013 (2013–2014)</td><td>Europe (0.99)</td><td>Asia</td></tr><tr><td>18</td><td>2014 (2013–2014)</td><td>North America (0.39)</td><td>Asia</td></tr><tr><td>19</td><td>2011 (2010–2011)</td><td>Europe (0.99)</td><td>North America</td></tr><tr><td>20</td><td>2012 (2011–2012)</td><td>North America (0.94)</td><td>South America</td></tr></tbody></table>

The introduction of clade C H13 HA in North America (event 19) was followed by an introduction to South America (event 20). Evidence for intercontinental gene flow among North American H13 IAVs occurred among eastern and western North American isolates (events 3, 12, 15, and 19). Clade B was composed almost exclusively of viruses circulating in North America, although one gene flow event to South America occurred recently (event 11).

The H16 HA was structured in at least two major genetic lineages (Fig. 2; also, see Fig. S4 and S5). The maximum-clade-credibility (MCC) tree was structured in three main clades (A to C) (see Fig. S5), while the maximum-likelihood (ML) tree provided support for only two main genetic clades (A and B/C merged) (see Fig. S4). The tMRCA of the H16 HA gene was dated to 1924 (1914–1932). Clade A included only viruses from Europe and was dated to 1977 (1975–1980); clade B included only viruses from North America, with a tMRCA estimated at 1969 (1967–1971). Our analyses supported the idea that the geographic origins of clades A and B were Europe and North America, respectively (posterior probabilities for the geographic origin of the MRCA, 0.99 for clade A and 1 for clade B). The tMRCA of clade C was estimated to be 1965 (1962–1968). Clade C may have arisen in Europe (posterior probability for the geographic origin of the MRCA, 0.87) and consisted of viruses of mixed origin, i.e., Europe, Asia, and North America.

Since the first isolation of an H16 IAV from a black-legged Kittiwake (*Rissa tridactyla*) in 1975, ten intercontinental gene flow events were identified for viruses of clade C (indicated by the numerals 1 to 10 in Fig. 2 and Table 2 and in Fig. S6). As for the H13 subtype, strong support for gene flow between Europe and North America was found, in particular from northwestern European countries (Denmark) to northeastern America (Delaware, New Hampshire, and Quebec) and from Iceland to Newfoundland (events 6 and 10). Evidence for intercontinental gene flow among North American H16 IAVs occurred among eastern and western North American isolates (events 3, 6, 8, and 10). In particular, intercontinental gene flow 8 seems to have been maintained in North America after its initial introduction in 2006 (2005–2006) for at least 10 years and may have replaced clade B of H16 HA (Fig. 2).

High rates of nucleotide substitution obtained for the H13 HA genetic lineages were consistent with those previously reported for H4, H6, and H7 subtypes circulating in wild ducks (Table 3) (15, 38). However, the nucleotide substitution rate of clade

<footnote><sup>a</sup>MRCA, most recent common ancestor; HPD, highest posterior density. Event numbers correspond to the numbers in Fig. 1 and Fig. S3.</footnote>
<figure><img src="image_4.png" /><figcaption>**FIG 2** Maximum-clade-credibility (MCC) tree for influenza A virus H16 hemagglutinin subtype (n = 192). Branches were colored according to most probable geographic origin (red, North America; orange, South America; dark blue, Europe; light blue, Asia; green, Oceania; gray, not identified). Black node bars represent the 95% highest posterior densities for times to the common ancestors. Numbers highlight intercontinental gene flow events as detailed in Table 2 and Fig. S6. Virus strain names and posterior probabilities are presented in Fig. S5.</figcaption></figure>

B—which consists exclusively of North American IAV—was lower than mean rates and HPD obtained for the other two H13 clades. The mean *dN/dS* (nonsynonymous substitutions/synonymous substitutions) ratios obtained for the three H13 genetic clades were comparable to those previously reported for other subtypes, which suggests that the HA was under strong purifying selection (Table 3). Nonetheless, a slightly higher *dN/dS* ratio obtained for clade B and C compared to other lineages suggests that they may be subjected to a more neutral selection. The mean nucleotide substitution and

<table><caption>**TABLE 2** Intercontinental gene flow events for influenza A virus H16 hemagglutinin<sup>a</sup></caption><thead><tr><th>H16 clade</th><th>Event</th><th>Time of MRCA ± 95% HPD</th><th>Geographic origin of the MRCA (posterior density)</th><th>Location of introduction</th></tr></thead><tbody><tr><td rowspan="10">C</td><td>1</td><td>1971 (1968–1972)</td><td>Europe (0.97)</td><td>Asia</td></tr><tr><td>2</td><td>1976 (1976–1976)</td><td>Asia (0.71)</td><td>Europe</td></tr><tr><td>3</td><td>1976 (1972–1980)</td><td>Europe (0.86)</td><td>North America</td></tr><tr><td>4</td><td>1999 (1999–1999)</td><td>Europe (1)</td><td>Asia</td></tr><tr><td>5</td><td>2003 (2002–2004)</td><td>Europe (1)</td><td>Asia</td></tr><tr><td>6</td><td>1999 (1998–2000)</td><td>Europe (0.99)</td><td>North America</td></tr><tr><td>7</td><td>2008 (2007–2009)</td><td>Europe (0.99)</td><td>Asia</td></tr><tr><td>8</td><td>2006 (2005–2006)</td><td>Europe (0.97)</td><td>North America</td></tr><tr><td>9</td><td>2006 (2006–2007)</td><td>North America (0.55)</td><td>South America</td></tr><tr><td>10</td><td>2008 (2007–2009)</td><td>Europe (0.63)</td><td>North America</td></tr></tbody></table>

<footnote><sup>a</sup>MRCA, most recent common ancestor; HPD, highest posterior density. Event numbers correspond to the numbers in Fig. 2 and Fig. S6.</footnote>
<table><caption>TABLE 3 Molecular evolution of the HA gene of influenza A virus subtypes H13 and H16</caption><thead><tr><th rowspan="2">Genetic lineage</th><th rowspan="2">N<sup>a</sup></th><th rowspan="2">Time period (yrs)</th><th colspan="3">Substitution rate<sup>b</sup></th></tr><tr><th>Mean</th><th>95% HPD</th><th>Mean dN/dS</th></tr></thead><tbody><tr><td>H13</td><td>338</td><td>40</td><td>3.8</td><td>3.6-4.1</td><td>0.13</td></tr><tr><td>H13 A</td><td>54</td><td>39</td><td>3.8</td><td>2.3-4.9</td><td>0.09</td></tr><tr><td>H13 B</td><td>76</td><td>39</td><td>0.8</td><td>0.6-1.0</td><td>0.18</td></tr><tr><td>H13 C</td><td>208</td><td>37</td><td>5.5</td><td>5.0-6.0</td><td>0.16</td></tr><tr><td>H16</td><td>192</td><td>41</td><td>3.1</td><td>2.8-3.4</td><td>0.09</td></tr><tr><td>H16 A</td><td>56</td><td>33</td><td>4.5</td><td>3.9-5.2</td><td>0.10</td></tr><tr><td>H16 B</td><td>19</td><td>35</td><td>4.6</td><td>3.9-5.2</td><td>0.06</td></tr><tr><td>H16 C</td><td>117</td><td>40</td><td>1.5</td><td>1.2-1.8</td><td>0.11</td></tr></tbody></table>

dN/dS ratios for the H16 gene were also consistent with H13 HA as well as with H4, H6, and H7 subtypes from wild ducks. However, H16 clade C (European mixed), which consisted of viruses of a geographically more mixed origin, had slightly lower nucleotide substitution rates and higher dN/dS ratios than clade A (European) and clade B (North American) (Table 3).

**Antigenic diversity between H13 and H16 LPAIVs.** As expected from two different HA subtypes, the H13 and H16 viruses formed two separate antigenic variants. The H13 and H16 viruses were generally well separated, forming groups on opposite sides of the antigenic map (Fig. 3; Table 4). A total of nine amino acid positions within or near the receptor binding site of the HA were identified that differed consistently between H13 and H16 viruses (based on alignments of 338 H13 and 192 H16 HA indicated in Table 5); of those, amino acid position 145 was located in the 130 loop, 200 and 208 were in the 190 helix, and 231 and 233 were in the 220 loop of the receptor binding site of the HA (HA numbering based on references 28 and 29). Of those, amino acid position 233 was listed previously as being involved in differences in receptor-binding site between HAs originating from Laridae and Anatidae (30). Additionally, the amino acid at position 196 differed between H13 (valine [V]) and H16 (aspartic acid [D]) viruses; this position may contribute to receptor binding specificity as identified previously based on crystal structures of H5 and H13 LPAIV (31). Due to nonspecific cross-reactivity, two H13 viruses (i.e., HEGU/AK/458/85 and HEGU/AK/479/85) had unexpected high titers against H16 antisera (Table 4); they were therefore positioned in the center of the map and served to pull H13 and H16 together.

**Antigenic diversity among H13 LPAIVs.** The representative H13 viruses formed at least two different antigenic variants (Fig. 3; Table 4). The viruses of H13 clades A and B were genetically distinct (Fig. 1) but were antigenically similar (Fig. 3), based on the H13 clade A antisera cross-reacting with H13 clade B viruses and vice versa. In contrast, H13 clade C viruses reacted poorly—if at all—with antisera that were raised against clade A and B viruses, and, conversely, antisera against clade C viruses rarely reacted with substantial titers with viruses of clade A and B. Thus, H13 clade A/B and H13 clade C viruses formed two different antigenic variants. The antigenic diversity of H13 clades A and B combined is about the same as the antigenic diversity of the H13 clade C. One H13 clade B virus, i.e., LAGU/DB/1370/86, could not be placed well in the map due to hemagglutination inhibition (HI) titers of 40 or lower (Table 4).

To gain insight into the molecular basis of the antigenic variation between H13 clades A/B and C, amino acids that differed consistently among the different clades of H13 viruses were identified (based on the alignment of 338 H13 [Table 5]). A total of four amino acid positions within or near the receptor binding site of the HA were identified that differed consistently for clades A, B, and/or C. Of those, amino acids at positions 149 and 254 differed consistently between clades A/B and C. Viruses belonging to clade C—except a single virus from South America that had an arginine (R) at position 149—had a deletion at position 149 (previously identified using a smaller data set as position 154 [12]), in contrast to viruses of clade A or B, which had an aspartic acid

<footnote><sup>a</sup>Number of nucleotide sequences included in the analysis.</footnote>

<footnote><sup>b</sup>Values are substitutions (10<sup>-3</sup>) per site per year. HPD, highest posterior density.</footnote>
<figure><img src="image_5.png" /><figcaption>**FIG 3** Antigenic map of H13 and H16 influenza A viruses (n = 44). Different subtypes and genetic clades are indicated with colors (yellow, H13 clade A; orange, H13 clade B; red, H13 clade C; blue, H16 clade A; purple, H16 clade B; green, H16 clade C). White circles indicate the antisera (n = 10). Virus strain names are abbreviated; the full names can be found in Table 4. Asterisks indicate antigens BHGU/NL/20/09, BHGU/SE/1/06, BHGU/SE/1/03, GBBG/AK/1421/79, BHGU/NL/1/07, HEGU/NY/AI00-532/00, and LAGU/NJ/AI08-0714/08, which had only two numerical HI titers to the tested sera, and hence, their placement in the map is not robust. In this map, the distance between the points represents antigenic distance as measured by the hemagglutination inhibition (HI) assay, in which the distances between antigens and antisera are inversely related to the log₂ HI titer. Each square in the grid of the antigenic map equals a 2-fold difference in the HI assay.</figcaption></figure>

(D), glutamic acid (E), asparagine (N), or serine (S) at this position. The correlation between the antigenic distance of H13 representative viruses from A/gull/MD/704/1977 (H13N6) (clade A)—the first detected H13 virus—and the number of HA1 amino acid substitutions from A/gull/MD/704/1977 was 0.87 and was statistically significant (*P* < 0.0001; Pearson correlation).

**Antigenic diversity among H16 LPAIVs.** The representative H16 viruses formed at least one antigenic variant (Fig. 3 and Table 4). The genetically distinct H16 clades A, B, and C did not form separate antigenic clusters in the map, which reflects the raw HI data, as there are no patterns for any of the four H16 antisera tested that correspond to the genetic lineages. The antigenic diversity of the H16 viruses is within eight antigenic units, with BHGU/NL/1/07 being on the edge of this antigenic space (i.e., low titers to all sera). The antigenic diversity of H16 clade A/B/C is about the same as the antigenic diversity of the H13 clades A and B combined and similar to the antigenic diversity of the H13 clade C.

Though clades A, B, and C did not form separate antigenic clusters in our analysis, amino acids that differed consistently among the different clades of H16 viruses were identified (based on the alignment of 192 H16 HA [Table 5]). A total of three amino acid positions within or near the receptor binding site of the HA were identified that differed
<table><caption>TABLE 4 Hemagglutinin inhibition data of H13 and H16 influenza A viruses (n = 44)</caption><thead><tr><th rowspan="3">Subtype<br>and clade</th><th rowspan="3">Virus isolate</th><th rowspan="3">Subtype</th><th rowspan="3">Virus abbreviation</th><th colspan="10">Antisera<sup>a</sup>:</th><th colspan="5">H16</th></tr><tr><th colspan="3">H13</th><th colspan="3">H14</th><th colspan="3">H15</th><th colspan="3">H16</th></tr><tr><th>BHG/U/NL/<br>2/07</th><th>GULL/MD/<br>704/77</th><th>B</th><th>GULL/MN/<br>1352/81</th><th>BRG/U/GE/<br>A100-2658/00</th><th>C</th><th>GBBG/AK/<br>1420/79</th><th>C</th><th>LAGU/NJ/<br>A108-714/08</th><th>A</th><th>BHG/U/SE/<br>2/99</th><th>BHG/U/SE/<br>A10-532/00</th><th>C</th><th>BHG/U/SE/<br>5/99</th><th>C</th><th>BHG/U/TM/<br>13/76</th></tr></thead><tbody><tr><td rowspan="10">H13/A</td><td>A/Black-headed gull/Netherlands/2/07</td><td>H13N6</td><td>BHG/U/NL/2/07</td><td>320</td><td>280</td><td>80</td><td>&lt;10</td><td>&lt;10</td><td>20</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>25</td></tr><tr><td>A/Black-headed gull/Netherlands/4/07</td><td>H13N6</td><td>BHG/U/NL/4/07</td><td>1,280</td><td>400</td><td>320</td><td>&lt;10</td><td>35</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>10</td><td>10</td><td>40</td><td>40</td><td>15</td></tr><tr><td>A/Black-headed gull/Netherlands/7/09</td><td>H13N2</td><td>BHG/U/NL/7/09</td><td>10</td><td>160</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>15</td></tr><tr><td>A/Black-headed gull/Sweden/10/05</td><td>H13N2</td><td>BHG/U/SE/10/05</td><td>240</td><td>320</td><td>40</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>15</td></tr><tr><td>A/Great-black headed gull/Sweden/1/03</td><td>H13N6</td><td>GBBG/SE/1/03</td><td>80</td><td>240</td><td>20</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/gull/MD/704/77</td><td>H13N6</td><td>GULL/MD/704/77</td><td>40</td><td>240</td><td>20</td><td>&lt;10</td><td>&lt;10</td><td>&lt;20</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td rowspan="10">H13/B</td><td>A/gull/MN/1352/81</td><td>H13N6</td><td>GULL/MN/1352/81</td><td>120</td><td>160</td><td>320</td><td>&lt;10</td><td>20</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/gull/NJ/34/92</td><td>H13N6</td><td>GULL/NJ/34/92</td><td>80</td><td>240</td><td>80</td><td>&lt;10</td><td>240</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/Herring gull/DB/13/90</td><td>H13N2</td><td>HEGU/DB/13/90</td><td>40</td><td>140</td><td>140</td><td>10</td><td>25</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/Laughing gull/DB/1370/86</td><td>H13N2</td><td>LAGU/DB/1370/86</td><td>10</td><td>40</td><td>&lt;10</td><td>10</td><td>40</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/ling-billed gull/GE/A100-2658/00</td><td>H13N6</td><td>RBGU/GE/A100-2658/00</td><td>60</td><td>40</td><td>640</td><td>15</td><td>15</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/ling-billed gull/MN/A10-1708/10</td><td>H13N6</td><td>RBGU/MN/A10-1708/10</td><td>80</td><td>200</td><td>120</td><td>10</td><td>10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td rowspan="10">H13/C</td><td>A/Black-headed gull/Netherlands/1/00</td><td>H13N8</td><td>BHG/U/NL/1/00</td><td>35</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>1,280</td><td>120</td><td>&lt;10</td><td>30</td><td>&lt;10</td><td>&lt;10</td><td>30</td><td>30</td><td>35</td></tr><tr><td>A/Black-headed gull/Netherlands/20/09</td><td>H13N2</td><td>BHG/U/NL/20/09</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>280</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>35</td><td>35</td></tr><tr><td>A/Black-headed gull/Netherlands/4/08</td><td>H13N8</td><td>BHG/U/NL/4/08</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>140</td><td>80</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>25</td></tr><tr><td>A/Black-headed gull/Sweden/1/03</td><td>H13N8</td><td>BHG/U/SE/1/03</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>560</td><td>40</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/Black-headed gull/Sweden/1/06</td><td>H13N8</td><td>BHG/U/SE/1/06</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>120</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/Black-headed gull/Sweden/2/03</td><td>H13N6</td><td>BHG/U/SE/1/03</td><td>10</td><td>&lt;10</td><td>10</td><td>30</td><td>160</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/Great-black headed gull/AK/1420/79</td><td>H13N2</td><td>BHG/U/SE/2/03</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>200</td><td>50</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/Great-black headed gull/AK/1421/79</td><td>H13N2</td><td>GBBG/AK/1420/79</td><td>10</td><td>35</td><td>10</td><td>&lt;10</td><td>160</td><td>10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/Great-black headed gull/AK/591/82</td><td>H13N2</td><td>GBBG/AK/591/82</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>140</td><td>80</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/Great-black headed gull/GJ/76/83</td><td>H13N2</td><td>GBBG/GJ/76/83</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>480</td><td>100</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/Herring gull/AK/458/85</td><td>H13N6</td><td>HEGU/AK/458/85</td><td>30</td><td>20</td><td>&lt;10</td><td>&lt;10</td><td>1,920</td><td>480</td><td>70</td><td>&lt;10</td><td>80</td><td>80</td><td>80</td><td>80</td><td>30</td></tr><tr><td>A/Herring gull/AK/479/85</td><td>H13N6</td><td>HEGU/AK/479/85</td><td>140</td><td>35</td><td>10</td><td>&lt;10</td><td>1,920</td><td>640</td><td>280</td><td>120</td><td>280</td><td>120</td><td>120</td><td>120</td><td>&lt;10</td></tr><tr><td>A/Laughing gull/NJ/A108-714/08</td><td>H13N9</td><td>LAGU/NJ/A108-714/08</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>320</td><td>560</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td rowspan="10">H16/A</td><td>A/Black-headed gull/Netherlands/5/07</td><td>H16N3</td><td>BHG/U/NL/5/07</td><td>35</td><td>25</td><td>&lt;10</td><td>&lt;10</td><td>140</td><td>&lt;10</td><td>960</td><td>160</td><td>320</td><td>640</td><td>40</td><td>40</td><td>40</td></tr><tr><td>A/Black-headed gull/Netherlands/1/07</td><td>H16N3</td><td>BHG/U/NL/1/07</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>80</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td></tr><tr><td>A/Black-headed gull/Netherlands/10/09</td><td>H16N3</td><td>BHG/U/NL/10/09</td><td>20</td><td>80</td><td>&lt;10</td><td>&lt;10</td><td>280</td><td>15</td><td>1,280</td><td>160</td><td>640</td><td>640</td><td>640</td><td>280</td><td>280</td></tr><tr><td>A/Black-headed gull/Netherlands/21/09</td><td>H16N3</td><td>BHG/U/NL/21/09</td><td>70</td><td>200</td><td>20</td><td>&lt;10</td><td>240</td><td>&lt;10</td><td>480</td><td>&lt;10</td><td>240</td><td>240</td><td>240</td><td>120</td><td>120</td></tr><tr><td>A/Black-headed gull/Networks/3/07</td><td>H16N3</td><td>BHG/U/NL/3/07</td><td>100</td><td>90</td><td>20</td><td>&lt;10</td><td>100</td><td>&lt;10</td><td>120</td><td>140</td><td>60</td><td>120</td><td>120</td><td>120</td><td>120</td></tr><tr><td>A/Black-headed gull/Sweden/2/99</td><td>H16N3</td><td>BHG/U/SE/2/99</td><td>10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>10</td><td>&lt;10</td><td>960</td><td>80</td><td>35</td><td>380</td><td>140</td><td>140</td><td>140</td></tr><tr><td>A/Black-headed gull/Sweden/8/05</td><td>H16N3</td><td>BHG/U/SE/8/05</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>10</td><td>&lt;10</td><td>1,280</td><td>&lt;10</td><td>30</td><td>30</td><td>120</td><td>20</td><td>1,600</td></tr><tr><td rowspan="10">H16/C</td><td>A/Herring gull/DB/2617/87</td><td>H16N3</td><td>HEGU/DB/2617/87</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>320</td><td>&lt;10</td><td>320</td><td>320</td><td>1920</td></tr><tr><td>A/Herring gull/INY/A10-532/00</td><td>H16N3</td><td>HEGU/INY/A10-532/00</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>80</td><td>20</td><td>20</td><td>1920</td><td>1920</td></tr><tr><td>A/Laughing gull/DB/2839/87</td><td>H16N3</td><td>LAGU/DB/2839/87</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>160</td><td>80</td><td>20</td><td>20</td><td>1920</td></tr><tr><td>A/Black-headed gull/Netherlands/26/09</td><td>H16N3</td><td>BHG/U/NL/26/09</td><td>10</td><td>25</td><td>&lt;10</td><td>&lt;10</td><td>20</td><td>&lt;10</td><td>30</td><td>80</td><td>20</td><td>1,600</td><td>400</td><td>1,600</td><td>1,600</td></tr><tr><td>A/Black-headed gull/Sweden/5/99</td><td>H16N3</td><td>BHG/U/SE/2/99</td><td>10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>70</td><td>&lt;10</td><td>560</td><td>30</td><td>320</td><td>100</td><td>400</td><td>1,600</td><td>400</td></tr><tr><td>A/Black-headed gull/TM/13/76</td><td>H16N3</td><td>BHG/U/TM/13/76</td><td>25</td><td>30</td><td>&lt;10</td><td>&lt;10</td><td>27.5</td><td>&lt;10</td><td>50</td><td>320</td><td>100</td><td>1,280</td><td>640</td><td>1,280</td><td>1,280</td></tr><tr><td>A/environment/Sweden/2/05</td><td>H16N3</td><td>ENV/SE/2/05</td><td>20</td><td>30</td><td>10</td><td>&lt;10</td><td>140</td><td>30</td><td>960</td><td>320</td><td>1,280</td><td>640</td><td>640</td><td>200</td><td>1,280</td></tr><tr><td>A/Little tern/Sweden/1/05</td><td>H16N3</td><td>LITE/SE/1/05</td><td>&lt;10</td><td>15</td><td>&lt;10</td><td>&lt;10</td><td>15</td><td>&lt;10</td><td>10</td><td>30</td><td>20</td><td>20</td><td>200</td><td>200</td><td>1,280</td></tr><tr><td>A/shorebird/DB/172/05</td><td>H16N3</td><td>SB/DB/172/05</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>30</td><td>&lt;10</td><td>240</td><td>60</td><td>60</td><td>200</td><td>200</td><td>200</td><td>1,280</td></tr><tr><td>A/shorebird/DB/195/06</td><td>H16N3</td><td>SB/DB/195/06</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>&lt;10</td><td>30</td><td>20</td><td>20</td><td>20</td><td>560</td><td>560</td></tr><tr><td>A/Slander-billed gull/AK/28/76</td><td>H16N3</td><td>SBGU/AK/28/76</td><td>20</td><td>140</td><td>10</td><td>&lt;10</td><td>50</td><td>&lt;10</td><td>80</td><td>160</td><td>100</td><td>100</td><td>100</td><td>1,280</td><td>1,280</td></tr></tbody></table>

<footnote>aBoldface indicates homologous titers.</footnote>
<table><caption>TABLE 5 Amino acid differences within or near the receptor binding site of the HA protein among H13 and H16 subtypes and clades<sup>a</sup></caption><thead><tr><th rowspan="2">Clade</th><th colspan="12">Amino acid(s) at position:</th></tr><tr><th>139</th><th>142</th><th>145</th><th>149</th><th>166</th><th>176</th><th>177</th><th>196</th><th>198</th><th>200</th><th>208</th><th>217</th><th>218</th><th>224</th><th>231</th><th>233</th></tr></thead><tbody><tr><td>H13 A</td><td>D</td><td>A, T, S</td><td>A</td><td>D, E, N, S</td><td>K, Q</td><td>K</td><td>T</td><td>V, L</td><td>V</td><td>E</td><td>S, G</td><td>K</td><td>S, L</td><td>K</td><td>P</td><td>Y</td></tr><tr><td>H13 B</td><td>D</td><td>A, T, S</td><td>A</td><td>D, N, S</td><td>K, R</td><td>G, R</td><td>T</td><td>V, I</td><td>T, A</td><td>E</td><td>S, G</td><td>S, R, N, H</td><td>S, L</td><td>K, N</td><td>P, L</td><td>Y, Q</td></tr><tr><td>H13 C</td><td>D</td><td>V, A</td><td>A</td><td>DEL, R</td><td>K, R, S</td><td>G, R</td><td>T, A, V</td><td>V, I</td><td>T, A, E</td><td>E</td><td>D, N, S</td><td>S, R, G</td><td>S, T</td><td>N, T, K</td><td>P</td><td>Y</td></tr><tr><td>H16 A</td><td>E</td><td>T</td><td>S</td><td>DEL</td><td>L</td><td>G</td><td>E</td><td>D</td><td>E</td><td>T</td><td>K</td><td>K</td><td>E</td><td>E</td><td>I</td><td>D</td></tr><tr><td>H16 B</td><td>D</td><td>V</td><td>S</td><td>DEL</td><td>DEL</td><td>G</td><td>D</td><td>D</td><td>E, ?</td><td>T, V</td><td>K</td><td>K, E</td><td>E</td><td>E</td><td>I</td><td>D, E, N</td></tr><tr><td>H16 C</td><td>D</td><td>V, A</td><td>S</td><td>DEL</td><td>K, DEL</td><td>G</td><td>E, D</td><td>D</td><td>E</td><td>T</td><td>K</td><td>K</td><td>E</td><td>E</td><td>I, V</td><td>D, N</td></tr></tbody></table>

consistently among the three H16 clades and were not associated with antigenic variation. The correlation between the antigenic distance of the representative viruses from A/Black-headed gull/TM/13/76 (H16N3) (clade C)—one of the first detected H16 viruses—and the number of HA1 amino acid substitutions from A/Black-headed gull/TM/13/76 was 0.67 and was statistically significant (P = 0.003; Pearson correlation).

DISCUSSION

We investigated the evolutionary history and intercontinental gene flow based on the hemagglutinin (HA) gene of H13 and H16 LPAIV and selected representative viruses from genetically distinct lineages to determine their antigenic properties by HI assays. H13 formed at least three distinct genetic clades as suggested previously based on smaller data sets (9, 32–35), while H16 formed at least two distinct genetic clades. Twenty and ten events of intercontinental gene flow were identified for H13 and for H16 viruses, respectively. At least two antigenic variants of H13 and at least one antigenic variant of H16 were identified. The presence of different antigenic variants among viruses of a single LPAIV subtype is in contrast to previous findings based on antigenic characterization of LPAIV H3 (26) and implies that antigenic variation within LPAIV subtypes occurs.

The frequency of intercontinental gene flow of the HA gene of H13 and H16 viruses was similar to that of the HA gene of H6 viruses but lower than that of internal genes (2, 27, 36, 37). Previously, intercontinental gene flow was described extensively for the H6 HA genes, while no intercontinental gene flow was detected for the H4 and H7 subtypes (15, 38). For the H6 subtype, gene flow has been described to have occurred ten times with four established genes during a period of 31 years (1975 to 2006) (15). Also, evidence for intercontinental gene flow among North American H13 and H16 genes occurred among eastern and western North American LPAIVs, in contrast to eastern North American LPAIVs only, as reported previously (39). Given the relatively high number of intercontinental flow events of IAV internal genes associated with shorebirds and gulls (2, 27, 36, 37), one might have expected a higher gene flow of gull-associated H13 and H16 HA genes, compared to, e.g., H6. However, a higher intercontinental gene flow was apparent only with H13 (i.e., 20 events during a period of 35 years). This may suggest one or more of the following: (i) H13 has a broader host range, host population size, and/or host distribution than H16; (ii) local H13-specific herd immunity is lower than H16-specific herd immunity and therefore less limiting of establishment opportunities in host populations of H13; (iii) H13 has higher environmental survival than H16; and (iv) introduced H13 HA genes may be less affected by strong subtype-dependent competition with endemic HA genes (e.g., with respect to linkage to NS1 and NP, as these contain most gull-specific features [33]) than introduced H16 genes. Interestingly, no H13 or H16 gene flow from Asia to Europe was described, which is in contrast to, e.g., H5 HPAIVs that have been introduced from Asia to Europe several times (40, 41). The relatively low frequency of detection of intercontinental gene flow of H13 or H16 genes out of North America and in particular Asia, relative to Europe, may be due to a bias in IAV surveillance and sequencing (i.e., the number of available IAV sequences from gulls isolated in Europe is higher than those from North America and in particular Asia).

<footnote><sup>a</sup>Based on the HA gene of H13 (n = 338) and H16 (n = 192) LPAIVs, including the 130 loop (positions 136 to 147 according to reference 29), the 190 helix (200 to 208), and the 220 loop (230 to 240). DEL, deletion of the amino acid.</footnote>
Antigenic diversity of LPAIVs depends partially on the host population size and structure. In this study, H13 and H16 LPAIVs formed at least three and two distinct genetic clades, respectively, that did not correspond or only partially corresponded with antigenic clusters. The H16 genetic clades did not form antigenic clusters, suggesting that clade-defining mutations were not in critical epitopes. In contrast, the H13 genetic clades partially corresponded with the antigenic variation of H13 LPAIV, suggesting that some of the clade-defining mutations were in critical epitopes. Also, given that the H13 antigenic space is larger than the antigenic space covered by H16 viruses, the host population of H13 may be larger and more widely distributed than the host population of H16 LPAIV, facilitating the circulation of more than one antigenic variant of a single LPAIV subtype. Strong genetic and antigenic divergence between two cocirculating lineages could be the product of a very large host meta-population size and relatively low cross-species transmission rate (42). Globally, viruses of the H13 subtype seem to be more common than viruses of the H16 subtype (2, 4), which is consistent with the finding that H13 LPAIV consists of multiple antigenic variants. Besides increased host population size and host distribution, prolonged virus survival may shape LPAIV epidemiology and evolution. Antigenic diversity within H13 LPAIV may be shaped by amino acid substitutions near the receptor binding site of the HA protein. In this study, we found evidence that amino acids or deletions at positions 149 and 254 of the HA protein may be involved in antigenic diversity among H13 strains. In addition, position 149 could be involved in H16 LPAIV antigenic diversity, as all H16 viruses had a deletion at this position and H16 clades A, B, and C were antigenically similar.

Cocirculating and newly introduced H13 or H16 LPAIV can be either antigenically similar or antigenically different. In the Northern hemisphere, H13 and H16 IAV subtypes circulate most extensively in breeding colonies in hatch-year birds at the end of summer and early fall (5-7). In black-headed gulls (which in Europe are one of the main hosts for H13 and H16 LPAIVs), infection with H13 or H16 results in strong protection against reinfection with the same virus; however, susceptibility to infection with the other subtype or with another strain of the same subtype is unknown (43, 44). Our findings support the independent long-term maintenance and cocirculation of at least two genetically distinct lineages of H13 and of H16 in Eurasia. This pattern is similar to the one that has been described for the H3 IAV subtype in ducks in North America (42). Our analysis showed that these genetically distinct cocirculating lineages may belong to the same antigenic variant. Here, we found evidence that genetically distinct cocirculating H13 or H16 LPAIV in a black-headed gull breeding colony site in the Netherlands may be either antigenically different (e.g., H13 clade A virus A/BHGU/NL/7/2009 [H13N2] and H13 clade C virus A/BHGU/NL/20/2009 [H13N2]) or antigenically similar (e.g., H16 clade A viruses A/BHGU/NL/10/2009 [H16N3] and A/BHGU/NL/21/2009 [H16N3] and H16 clade C virus A/BHGU/NL/26/2009 [H16N3]). Similarly, intercontinental gene flow occurred with HA genes that were antigenically similar to local circulating viruses (i.e., H16 clade C viruses that were genetically closely related to SB/DB/172/06 and SB/DB/195/06 versus local circulating H16 clade B viruses) and HA genes that were antigenically different from local circulating viruses (i.e., H13 clade C viruses genetically closely related to LAGU/NJ/AI08-0714/08 versus local circulating H13 clade B viruses.

Antigenic variation within an LPAIV subtype at the clade level (i.e., H13 clades A and B combined versus H13 clade C) was described here, yet less is known about antigenic variation within genetic clades of H13, H16, or other LPAIV subtypes. For H13, genetic diversity within clades seemed stable—i.e., viruses of clade A, B, or C collected over 3 decades were antigenically closely related—suggesting no major genetic differences; this is in contrast to the few mutations needed for antigenic change in seasonal human IAV. Similarly, a study on antigenic variation of H3 LPAIV isolated in North America suggested that genetically diverse viruses were antigenically stable (26). Major antigenic changes in seasonal human IAV were due to amino acid substitutions immediately adjacent to the receptor binding site (18); this could potentially also explain antigenic variation between antigenically different viruses of H13 clade A/B and clade
<table><caption>**TABLE 6** Representative viruses selected to generate ferret antisera used to map the antigenic diversity of H13 and H16 influenza A viruses</caption><thead><tr><th>Subtype</th><th>Clade</th><th>Virus strain name</th></tr></thead><tbody><tr><td rowspan="6">H13</td><td>A</td><td>A/gull/Maryland/704/1977 (H13N6)</td></tr><tr><td>A</td><td>A/Black-headed gull/Netherlands/2/2007 (H13N6)</td></tr><tr><td>B</td><td>A/Ring-billed gull/Georgia/AI00-2658/2000 (H13N6)</td></tr><tr><td>B</td><td>A/gull/Minnesota/1352/1981 (H13N6)</td></tr><tr><td>C</td><td>A/Laughing gull/New Jersey/AI08-0714/2008 (H13N9)</td></tr><tr><td>C</td><td>A/Great black-headed gull/Astrakhan/1420/1979 (H13N2)</td></tr><tr><td rowspan="4">H16</td><td>A</td><td>A/Black-headed gull/Sweden/2/1999 (H16N3)</td></tr><tr><td>B</td><td>A/Herring gull/New York/AI00-532/2000 (H16N3)</td></tr><tr><td>C</td><td>A/Black-headed gull/Turkmenistan/13/1976 (H16N3)</td></tr><tr><td>C</td><td>A/Black-headed gull/Sweden/5/1999 (H16N3)</td></tr></tbody></table>

C (i.e., amino acid position 149 of the HA). Future work on antigenic variation of LPAIV should include within-clade genetic and antigenic variation.

## MATERIALS AND METHODS

**Viruses.** The HA sequences of H13 (n = 64) and H16 (n = 20) viruses isolated from wild birds in North America (n = 39 and n = 5, respectively) and Europe (n = 25 and n = 15, respectively) between 1976 and 2010 were determined at the University of Minnesota (Saint Paul, MN) and at the Department of Viroscience of the Erasmus Medical Center (Rotterdam, the Netherlands). Details on virus isolates, including GenBank accession numbers, are summarized in Tables S2 and S3; details related to the Sanger sequencing method are available upon request. The HA sequences were supplemented with full-length nucleotide sequences of the HA gene of H13 and H16 viruses isolated from wild birds between 1975 and 2017 and downloaded from GenBank (https://www.ncbi.nlm.nih.gov). The full data set included sequences of H13 (n = 519) and H16 (n = 276) HA genes and was biased toward virus strains collected since 2000 due to increased surveillance and sequencing since 2000.

Of this full data set, viruses representing the genetically distinct clades were selected (n = 44; H13 clades A, B, and C and H16 clades A, B, and C; see Results for clade definition) to investigate the antigenic diversity of H13 and H16 viruses. Of those viruses, viruses that were genetically most divergent were selected (n = 10) to generate ferret antisera (Table 6). The antigenic properties of all representative viruses (n = 44) were analyzed in hemagglutination inhibition (HI) assays using a panel of ten ferret antisera.

**Genetic analyses.** The nucleotide sequences of the coding region of H13 and H16 HA were aligned with the program CLC 8.0 (CLC bio, Aarhus, Denmark). Neighbor-joining trees were then generated, with 1,000 bootstraps, in order to assess the overall genetic structure of the H13 (n = 519) and H16 (n = 276) HA sequences. To lower the bias in species and geography (e.g., black-headed gulls [Chroicocephalus ridibundus] from the Netherlands and glaucous-winged gulls [Larus glaucescens] from Alaska), duplicate sequences (i.e., identical sequences of the same host species, location, and date) were identified with mothur 1.39.5 (45) and removed, resulting in final alignments of H13 (n = 338) and H16 (n = 192) HA. To identify the genetic structure of H13 and H16 virus subtypes, maximum-likelihood trees with 1,000 bootstraps were generated with the software PhyML 3.1 (46). The general time-reversible (GTR) evolutionary model, an estimation of the proportion of invariable sites (I) and of the nucleotide heterogeneity of substitution rate (α), was used as selected by Model Generator 0.85 (47). To investigate the evolutionary history of H13 and H16 virus subtypes, Bayesian Markov chain Monte Carlo coalescent analyses were performed. The temporal structure of the data set was assessed with the program TempEst 1.5.3 (48). Both data sets showed a positive correlation between genetic divergence and sampling time and appear to be suitable for phylogenetic molecular clock analyses. Time to the most recent common ancestors (MRCA) as well as geographic ancestral states (i.e., continent) and their associated posterior probabilities were obtained based on the method described by Lemey et al. with the program BEAST 1.10.1 (49, 50). A strict molecular clock model was selected, as relaxed clock models (uncorrelated exponential and uncorrelated log-normal) resulted in low effective sample sizes (< 200) in spite of high chain length (> 200 million states). In all simulations, a Bayesian skyline coalescent tree prior (51) was selected. The Shapiro-Rambaut-Drummond-2006 nucleotide substitution model was selected (52), and it has been used in population dynamic studies of other IAV subtypes (15, 38, 42, 53). Overall, a method similar to that in previous studies on IAV evolutionary dynamics of subtypes H4, H6, and H7 (15, 38, 54) was used. Analyses were performed with two independent chain lengths of 100 million generations sampled every 1,000 iterations; the first 10% of trees were discarded as burn-in. Substitutions rates based on independent analyses of the major H13 and H16 clades were obtained using the program BEAST 1.10.1. Rates of nonsynonymous substitutions (dN) and synonymous substitutions (dS) were obtained using the single likelihood ancestor counting method implemented in HyPhy (55). Computations were performed with the Datamonkey Web server (56, 57).

**Antisera.** Postinfection antisera were prepared upon nasal inoculation of ferrets (>1 year of age; male; two ferrets per virus) with virus (cultured on embryonated chicken eggs; each ferret received 106 to 107 median egg infectious doses [EID$_{50}$]/100 µl) and blood collection by exsanguination 14 days later.
An overview of antisera used in this study is provided in Table 6. Antisera were pretreated overnight at 37°C with receptor-destroying enzyme (Vibrio cholerae neuraminidase), followed by inactivation for 1 h at 56°C before use in HI assays.

**Antigenic analyses.** HI assays were performed according to standard procedures (58). The HI titer is expressed as the reciprocal value of the highest serum dilution that completely inhibited hemagglutination. To investigate antigenic variation among and within H13 and H16 viruses, antigenic cartography methods were used as described previously (19). Briefly, antigenic cartography is a method to analyze and visualize HI assay data. The titers in an HI table can be thought of as specifying target distances between antigens and antisera. In an antigenic map, the distance between antigen point A and antiserum point S corresponds to the difference between the log₂ value of the maximum observed titer to antiserum S from any antigen and the titer of antigen A to antiserum S. Modified multidimensional scaling methods are used to arrange the antiserum and antigen points in an antigenic map to best satisfy the target distances specified by the HI data (18). Because antigens are tested against multiple antisera and antisera are tested against multiple antigens, many measurements can be used to determine the position of the antigens and antisera in an antigenic map, thus improving the resolution of the HI data.

**Ethics statement.** This study was approved by the independent animal experimentation ethical review committee Stichting DEC consult (Erasmus MC permits 122-98-01, 122-08-04, and 15-340-03) and was performed under animal biosafety level 2 (ABSL2) conditions. Animal welfare was monitored daily, and all animal handling was performed under light anesthesia (ketamine) to minimize animal discomfort.

**Data availability.** Sequences are available in GenBank under accession numbers KF612922 to KF612965, KR087564, KR087572, KR087577 to KR087595, KR087597 to KR087601, KR087604 to KR087615, MK027211, and MK027212.

## SUPPLEMENTAL MATERIAL

Supplemental material is available online only.

**SUPPLEMENTAL FILE 1, PDF file, 6.1 MB.**

## ACKNOWLEDGMENTS

This work was funded by the Swedish Research Council Vetenskapsrådet (2015-03877) and the National Institute of Allergy and Infectious Diseases (NIAID), National Institutes of Health (NIH), Department of Health and Human Services, under contract no. HHSN266200700007C and HHSN272201400008C. C.L. is supported by Chaire mixte: Université de La Réunion-INSERM. The funding agencies did not have any involvement in the study design, implementation, or publishing of this study, and the research presented herein represents the opinions of the authors but not necessarily the opinions of the funding agencies.

We gratefully acknowledge the following researchers for sharing and preparing virus isolates and sequences, among others: Scott Krauss, Janice C. Pedersen, Shinichiro Enomoto, Justin D. Brown, Jonathan Runstadler, Nichola Hill, Nicola Lewis, Alexander Shestopalov, Neus Latorre-Margalef, and Jonas Waldenström.

## REFERENCES

1. Olsen B, Munster VJ, Wallensten A, Waldenstrom J, Osterhaus AD, Fouchier RA. 2006. Global patterns of influenza A virus in wild birds. Science 312:384–388. https://doi.org/10.1126/science.1122438.

2. Wille M, Robertson GJ, Whitney H, Bishop MA, Runstadler JA, Lang AS. 2011. Extensive geographic mosaicism in avian influenza viruses from gulls in the northern hemisphere. PLoS One 6:e20664. https://doi.org/10.1371/journal.pone.0020664.

3. Lang AS, Lebarbenchon C, Ramey AM, Robertson GJ, Waldenstrom J, Wille M. 2016. Assessing the role of seabirds in the ecology of influenza A viruses. Avian Dis 60:378–386. https://doi.org/10.1637/11135-050815 -RegR.

4. Zhang Y, Aevermann BD, Anderson TK, Burke DF, Dauphin G, Gu Z, He S, Kumar S, Larsen CN, Lee AJ, Li X, Macken C, Mahaffey C, Pickett BE, Reardon B, Smith T, Stewart L, Suloway C, Sun G, Tong L, Vincent AL, Walters B, Zaremba S, Zhao H, Zhou L, Zmasek C, Klem EB, Scheuermann RH. 2017. Influenza Research Database: an integrated bioinformatics resource for influenza virus research. Nucleic Acids Res 45:D466–D474. https://doi.org/10.1093/nar/gkw857.

5. Verhagen JH, Majoor F, Lexmond P, Vuong O, Kasemir G, Lutterop D, Osterhaus AD, Fouchier RA, Kuiken T. 2014. Epidemiology of influenza A virus among black-headed gulls, the Netherlands, 2006–2010. Emerg Infect Dis 20:138–141. https://doi.org/10.3201/eid2001.130984.

6. Velarde R, Calvin SE, Ojkic D, Barker IK, Nagy E. 2010. Avian influenza virus H13 circulating in ring-billed gulls (Larus delawarensis) in southern Ontario, Canada. Avian Dis 54:411–419. https://doi.org/10.1637/8808 -040109-Reg.1.

7. Graves IL. 1992. Influenza viruses in birds of the Atlantic flyway. Avian Dis 36:1–10. https://doi.org/10.2307/1591706.

8. Lewis NS, Javakhishvili Z, Russell CA, Machabishvili A, Lexmond P, Verhagen JH, Vuong O, Onashvili T, Donduashvili M, Smith DJ, Fouchier RA. 2013. Avian influenza virus surveillance in wild birds in Georgia: 2009–2011. PLoS One 8:e58534. https://doi.org/10.1371/journal.pone .0058534.

9. Huang Y, Wille M, Benkaroun J, Munro H, Bond AL, Fifield DA, Robertson GJ, Ojkic D, Whitney H, Lang AS. 2014. Perpetuation and reassortment of gull influenza A viruses in Atlantic North America. Virology 456–457: 353–363. https://doi.org/10.1016/j.virol.2014.04.009.

10. Obenauer JC, Denson J, Mehta PK, Su X, Mukatira S, Finkelstein DB, Xu X, Wang J, Ma J, Fan Y, Rakestraw KM, Webster RG, Hoffmann E, Krauss S, Zheng J, Zhang Z, Naeve CW. 2006. Large-scale sequence analysis of avian influenza isolates. Science 311:1576–1580. https:// doi.org/10.1126/science.1121586.

11. Van Borm S, Rosseel T, Vangeluwe D, Vandenbussche F, van den Berg T, Lambrecht B. 2012. Phylogeographic analysis of avian influenza viruses
isolated from Charadriiformes in Belgium confirms intercontinental re- assortment in gulls. Arch Virol 157:1509-1522. https://doi.org/10.1007/ s00705-012-1323-x.

12. Wille M, Robertson GJ, Whitney H, Ojkic D, Lang AS. 2011. Reassortment of American and Eurasian genes in an influenza A virus isolated from a great black-backed gull (Larus marinus), a species demonstrated to move between these regions. Arch Virol 156:107-115. https://doi.org/ 10.1007/s00705-010-0839-1.

13. Ratanakorn P, Wiratsudakul A, Wiriyaratte W, Eiamampai K, Farmer AH, Webster RG, Chaichoune K, Suwanpakdee S, Pothieng D, Puthathavathana P. 2012. Satellite tracking on the flyways of brown-headed gulls and their potential role in the spread of highly pathogenic avian influenza H5N1 virus. PLoS One 7:e49939. https://doi.org/10.1371/journal.pone .0049939.

14. Hall JS, Teslaa JL, Nashold SW, Halpin RA, Stockwell T, Wentworth DE, Dugan V, Ip HS. 2013. Evolution of a reassortant North American gull influenza virus lineage: drift, shift and stability. Virol J 10:179. https:// doi.org/10.1186/1743-422X-10-179.

15. Bahl J, Vijaykrishna D, Holmes EC, Smith GJ, Guan Y. 2009. Gene flow and competitive exclusion of avian influenza A virus in natural reservoir hosts. Virology 390:289-297. https://doi.org/10.1016/j.virol.2009.05.002.

16. Webster RG, Bean WJ, Gorman OT, Chambers TM, Kawaoka Y. 1992. Evolution and ecology of influenza A viruses. Microbiol Rev 56:152-179.

17. Latorre-Margalef N, Grosbois V, Wahlgren J, Munster VJ, Tolf C, Fouchier RAM, Osterhaus ADME, Olsen B, Waldenström J. 2013. Heterosubtypic immunity to influenza A virus infections in mallards may explain existence of multiple virus subtypes. PLoS Pathog 9:e1003443. https://doi .org/10.1371/journal.ppat.1003443.

18. Koel BF, Burke DF, Bestebroer TM, van der Vliet S, Zondag GC, Vervaet G, Skepner E, Lewis NS, Spronken MI, Russell CA, Eropkin MY, Hurt AC, Barr IG, de Jong JC, Rimmelzwaan GF, Osterhaus AD, Fouchier RA, Smith DJ. 2013. Substitutions near the receptor binding site determine major antigenic change during influenza virus evolution. Science 342:976-979. https://doi.org/10.1126/science.1244730.

19. Smith DJ, Lapedes AS, de Jong JC, Bestebroer TM, Rimmelzwaan GF, Osterhaus AD, Fouchier RA. 2004. Mapping the antigenic and genetic evolution of influenza virus. Science 305:371-376. https://doi.org/10. 1126/science.1097211.

20. Lorusso A, Vincent AL, Harland ML, Alt D, Bayles DO, Swenson SL, Gramer MR, Russell CA, Smith DJ, Lager KM, Lewis NS. 2011. Genetic and antigenic characterization of H1 influenza viruses from United States swine from 2008. J Gen Virol 92:919-930. https://doi.org/10.1099/vir.0 .027557-0.

21. Marozin S, Gregory V, Cameron K, Bennett M, Valette M, Aymard M, Foni E, Barigazzi G, Lin Y, Hay A. 2002. Antigenic and genetic diversity among swine influenza A H1N1 and H1N2 viruses in Europe. J Gen Virol 83: 735-745. https://doi.org/10.1099/0022-1317-83-4-735.

22. Cong YL, Pu J, Liu QF, Wang S, Zhang GZ, Zhang XL, Fan WX, Brown EG, Liu JH. 2007. Antigenic and genetic characterization of H9N2 swine influenza viruses in China. J Gen Virol 88:2035-2041. https://doi.org/10. 1099/vir.0.82783-0.

23. Lewis NS, Daly JM, Russell CA, Horton DL, Skepner E, Bryant NA, Burke DF, Rash AS, Wood JL, Chambers TM, Fouchier RA, Mumford JA, Elton DM, Smith DJ. 2011. Antigenic and genetic evolution of equine influenza A (H3N8) virus from 1968 to 2007. J Virol 85:12742-12749. https://doi .org/10.1128/JVI.05319-11.

24. Cattoli G, Milani A, Temperton N, Zecchin B, Buratin A, Molesti E, Aly MM, Arafa A, Capua I. 2011. Antigenic drift in H5N1 avian influenza virus in poultry is driven by mutations in major antigenic sites of the hemagglutinin molecule analogous to those for human influenza virus. J Virol 85:8718-8724. https://doi.org/10.1128/JVI.02403-10.

25. Koel BF, van der Vliet S, Burke DF, Bestebroer TM, Bharoto EE, Yasa IW, Herliana I, Laksono BM, Xu K, Skepner E, Russell CA, Rimmelzwaan GF, Perez DR, Osterhaus AD, Smith DJ, Prajitno TY, Fouchier RA. 2014. Antigenic variation of clade 2.1 H5N1 virus is determined by a few amino acid substitutions immediately adjacent to the receptor binding site. mBio 5:e01070-14. https://doi.org/10.1128/mBio.01070-14.

26. Bailey E, Long LP, Zhao N, Hall JS, Baroch JA, Nolting J, Senter L, Cunningham FL, Pharr GT, Hanson L, Slemons R, DeLiberto TJ, Wan XF. 2016. Antigenic characterization of H3 subtypes of avian influenza A viruses from North America. Avian Dis 60:346-353. https://doi.org/10.1637/11086-041015-RegR.

27. Krauss S, Obert CA, Franks J, Walker D, Jones K, Seiler P, Niles L, Pryor SP, Obenauer JC, Naeve CW, Widjaja L, Webby RJ, Webster RG. 2007.

Influenza in migratory birds and evidence of limited intercontinental virus exchange. PLoS Pathog 3:e167. https://doi.org/10.1371/journal .ppat.0030167.

28. de Graaf M, Fouchier RA. 2014. Role of receptor binding specificity in influenza A virus transmission and pathogenesis. EMBO J 33:823-841. https://doi.org/10.1002/embj.201387442.

29. Burke DF, Smith DJ. 2014. A recommended numbering scheme for influenza A HA subtypes. PLoS One 9:e112302. https://doi.org/10.1371/ journal.pone.0112302.

30. Yamnikova SS, Gambaryan AS, Tuzikov AB, Bovin NV, Matrosovich MN, Fedyakina IT, Grinev AA, Blinov VM, Lvov DK, Suarez DL, Swayne DE. 2003. Differences between HA receptor-binding sites of avian influenza viruses isolated from Laridae and Anatidae. Avian Dis 47:1164-1168. https://doi.org/10.1637/0005-2086-47.3.1164.

31. Lu X, Qi J, Shi Y, Wang M, Smith DF, Heimburg-Molinaro J, Zhang Y, Paulson JC, Xiao H, Gao GF. 2013. Structure and receptor binding specificity of hemagglutinin H13 from avian influenza A virus H13N6. J Virol 87:9077-9085. https://doi.org/10.1128/JVI.00235-13.

32. Benkarino J, Shoham D, Kroyer ANK, Whitney H, Lang AS. 2016. Analysis of influenza A viruses from gulls: an evaluation of inter-regional movements and interactions with other avian and mammalian influenza A viruses. Cogent Biol 2:1234957. https://doi.org/10.1080/23312025.2016.1234957.

33. Tonnessen R, Hauge AG, Hansen EF, Rimstad E, Jonassen CM. 2013. Host restrictions of avian influenza viruses: in silico analysis of H13 and H16 specific signatures in the internal proteins. PLoS One 8:e63270. https:// doi.org/10.1371/journal.pone.0063270.

34. lamnikova SS, Gambaryan AS, Aristova VA, L'Vov DK, Lomakina NF, Munster V, Lexmond P, Foucher RA. 2009. A/H13 and A/H16 influenza viruses: different lines of one precursors. Vopr Virusol 54:10-18. (In Russian.)

35. Wang ZJ, Kikutani Y, Nguyen LT, Hiono T, Matsuno K, Okamatsu M, Krauss S, Webby R, Lee YJ, Kida H, Sakoda Y. 2018. H13 influenza viruses in wild birds have undergone genetic and antigenic diversification in nature. Virus Genes 54:543-549. https://doi.org/10.1007/ s11262-018-1573-0.

36. Dugan VG, Chen R, Spiro DJ, Sengamalay N, Zaborsky J, Ghedin E, Nolting J, Swayne DE, Runstadler JA, Happ GM, Senne DA, Wang R, Slemons RD, Holmes EC, Taubenberger JK. 2008. The evolutionary genetics and emergence of avian influenza viruses in wild birds. PLoS Pathog 4:e1000076. https://doi.org/10.1371/journal.ppat.1000076.

37. Fries AC, Nolting JM, Danner A, Webster RG, Bowman AS, Krauss S, Slemons RD. 2013. Evidence for the circulation and inter-hemispheric movement of the H14 subtype influenza A virus. PLoS One 8:e59216. https://doi.org/10.1371/journal.pone.0059216.

38. Lebarbenchon C, Stallknecht DE. 2011. Host shifts and molecular evolution of H7 avian influenza virus hemagglutinin. Virol J 8:328. https://doi.org/10.1186/1743-422X-8-328.

39. Pearce JM, Ramey AM, Ip HS, Gill RE, Jr. 2010. Limited evidence of trans-hemispheric movement of avian influenza viruses among contemporary North American shorebird isolates. Virus Res 148:44-50. https:// doi.org/10.1016/j.viruses.2009.12.002.

40. The Global Consortium for H5N8 and Related Influenza Viruses. 2016. Role for migratory wild birds in the global spread of avian influenza H5N8. Science 354:213-217. https://doi.org/10.1126/science.aaf8852.

41. Adlhoch C, Gossner C, Koch G, Brown I, Bouwstra R, Verdonck F, Penttinen P, Harder T. 2014. Comparing introduction to Europe of highly pathogenic avian influenza viruses A(H5N8) in 2014 and A(H5N1) in 2005. Eurosurveillance 19:20996. https://doi.org/10.2807/ 1560-7917.ES2014.19.50.20996.

42. Bahl J, Krauss S, Kuhnert D, Fourment M, Raven G, Pryor SP, Niles LJ, Danner A, Walker D, Mendenhall IH, Su YC, Dugan VG, Halpin RA, Stockwell TB, Webby RJ, Wentworth DE, Drummond AJ, Smith GJ, Webster RG. 2013. Influenza A virus migration and persistence in North American wild birds. PLoS Pathog 9:e1003570. https://doi.org/10.1371/ journal.ppat.1003570.

43. Verhagen JH, Hofle U, van Amerongen G, van de Bildt M, Majoor F, Fouchier RA, Kuiken T. 2015. Long-term effect of serial infections with H13 and H16 low-pathogenic avian influenza viruses in black-headed gulls. J Virol 89:11507-11522. https://doi.org/10.1128/JVI.01765-15.

44. Brown J, Poulson R, Carter D, Lebarbenchon C, Pantin-Jackwood M, Spackman E, Shepherd E, Killian M, Stallknecht D. 2012. Susceptibility of avian species to North American H13 low pathogenic avian influ-
enza viruses. Avian Dis 56:969-975. https://doi.org/10.1637/10158-040912-Reg.1.

45. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thallinger GG, Van Horn DJ, Weber CF. 2009. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol 75:7537-7541. https://doi.org/10.1128/AEM.01541-09.

46. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. 2010. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol 59: 307-321. https://doi.org/10.1093/sysbio/syq010.

47. Keane TM, Creevey CJ, Pentony MM, Naughton TJ, McLnerney JO. 2006. Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evol Biol 6:29. https://doi.org/10.1186/1471-2148-6-29.

48. Rambaut A, Lam TT, Carvalho LM, Pybus OG. 2016. Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evol 2:vew007. https://doi.org/10.1093/ve/vew007.

49. Drummond AJ, Rambaut A. 2007. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 7:214. https://doi.org/10.1186/1471-2148-7-214.

50. Lemey P, Rambaut A, Drummond AJ, Suchard MA. 2009. Bayesian phylogeography finds its roots. PLoS Comput Biol 5:e1000520. https://doi.org/10.1371/journal.pcbi.1000520.

51. Drummond AJ, Rambaut A, Shapiro B, Pybus OG. 2005. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol 22:1185-1192. https://doi.org/10.1093/molbev/msi103.

52. Shapiro B, Rambaut A, Drummond AJ. 2006. Choosing appropriate substitution models for the phylogenetic analysis of protein-coding sequences. Mol Biol Evol 23:7-9. https://doi.org/10.1093/molbev/msj021.

53. Vijaykrishna D, Bahl J, Riley S, Duan L, Zhang JX, Chen H, Peiris JS, Smith GJ, Guan Y. 2008. Evolutionary dynamics and emergence of panzootic H5N1 influenza viruses. PLoS Pathog 4:e1000161. https://doi.org/10.1371/journal.ppat.1000161.

54. Lebarbenchon C, Brown JD, Stallknecht DE. 2013. Evolution of influenza A virus H7 and N9 subtypes, Eastern Asia. Emerg Infect Dis 19: 1635-1638. https://doi.org/10.3201/eid1910.130609.

55. Kosakovsky Pond SL, Frost SD, Muse SV. 2005. HyPhy: hypothesis testing using phylogenies. Bioinformatics 21:676-679. https://doi.org/10.1093/bioinformatics/bti079.

56. Delport W, Poon AF, Frost SD, Kosakovsky Pond SL. 2010. Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics 26:2455-2457. https://doi.org/10.1093/bioinformatics/btq429.

57. Kosakovsky Pond SL, Frost SD. 2005. Not so different after all: a comparison of methods for detecting amino acid sites under selection. Mol Biol Evol 22:1208-1222. https://doi.org/10.1093/molbev/msi105.

58. Hirst GK. 1943. Studies of antigenic differences among strains of influenza A by means of red cell agglutination. J Exp Med 78:407-423. https://doi.org/10.1084/jem.78.5.407.