Identification, characterization, and natural selection of mutations driving airborne transmission of A/H5N1 virus.

Recently, A/H5N1 influenza viruses were shown to acquire airborne transmissibility between ferrets upon targeted mutagenesis and virus passage. The critical genetic changes in airborne A/Indonesia/5/05 were not yet identified. Here, five substitutions proved to be sufficient to determine this airborne transmission phenotype. Substitutions in PB1 and PB2 collectively caused enhanced transcription and virus replication. One substitution increased HA thermostability and lowered the pH of membrane fusion. Two substitutions independently changed HA binding preference from α2,3-linked to α2,6-linked sialic acid receptors. The loss of a glycosylation site in HA enhanced overall binding to receptors. The acquired substitutions emerged early during ferret passage as minor variants and became dominant rapidly. Identification of substitutions that are essential for airborne transmission of avian influenza viruses between ferrets and their associated phenotypes advances our fundamental understanding of virus transmission and will increase the value of future surveillance programs and public health risk assessments.

Defining the risk of SARS-CoV-2 variants on immune protection.

The global emergence of many severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants jeopardizes the protective antiviral immunity induced after infection or vaccination. To address the public health threat caused by the increasing SARS-CoV-2 genomic diversity, the National Institute of Allergy and Infectious Diseases within the National Institutes of Health established the SARS-CoV-2 Assessment of Viral Evolution (SAVE) programme. This effort was designed to provide a real-time risk assessment of SARS-CoV-2 variants that could potentially affect the transmission, virulence, and resistance to infection- and vaccine-induced immunity. The SAVE programme is a critical data-generating component of the US Government SARS-CoV-2 Interagency Group to assess implications of SARS-CoV-2 variants on diagnostics, vaccines and therapeutics, and for communicating public health risk. Here we describe the coordinated approach used to identify and curate data about emerging variants, their impact on immunity and effects on vaccine protection using animal models. We report the development of reagents, methodologies, models and notable findings facilitated by this collaborative approach and identify future challenges. This programme is a template for the response to rapidly evolving pathogens with pandemic potential by monitoring viral evolution in the human population to identify variants that could reduce the effectiveness of countermeasures.

Species-specific emergence of H7 highly pathogenic avian influenza virus is driven by intrahost selection differences between chickens and ducks.

Highly pathogenic avian influenza viruses (HPAIVs) cause severe hemorrhagic disease in terrestrial poultry and are a threat to the poultry industry, wild life, and human health. HPAIVs arise from low pathogenic avian influenza viruses (LPAIVs), which circulate in wild aquatic birds. HPAIV emergence is thought to occur in poultry and not wild aquatic birds, but the reason for this species-restriction is not known. We hypothesized that, due to species-specific tropism and replication, intrahost HPAIV selection is favored in poultry and disfavored in wild aquatic birds. We tested this hypothesis by co-inoculating chickens, representative of poultry, and ducks, representative of wild aquatic birds, with a mixture of H7N7 HPAIV and LPAIV, mimicking HPAIV emergence in an experimental setting. Virus selection was monitored in swabs and tissues by RT-qPCR and immunostaining of differential N-terminal epitope tags that were added to the hemagglutinin protein. HPAIV was selected in four of six co-inoculated chickens, whereas LPAIV remained the major population in co-inoculated ducks on the long-term, despite detection of infectious HPAIV in tissues at early time points. Collectively, our data support the hypothesis that HPAIVs are more likely to be selected at the intrahost level in poultry than in wild aquatic birds and point towards species-specific differences in HPAIV and LPAIV tropism and replication levels as possible explanations.

Substitutions near the HA receptor binding site explain the origin and major antigenic change of the B/Victoria and B/Yamagata lineages.

Influenza B virus primarily infects humans, causing seasonal epidemics globally. Two antigenic variants-Victoria-like and Yamagata-like-were detected in the 1980s, of which the molecular basis of emergence is still incompletely understood. Here, the antigenic properties of a unique collection of historical virus isolates, sampled from 1962 to 2000 and passaged exclusively in mammalian cells to preserve antigenic properties, were determined with the hemagglutination inhibition assay and an antigenic map was built to quantify and visualize the divergence of the lineages. The antigenic map revealed only three distinct antigenic clusters-Early, Victoria, and Yamagata-with relatively little antigenic diversity in each cluster until 2000. Viruses with Victoria-like antigenic properties emerged around 1972 and diversified subsequently into two genetic lineages. Viruses with Yamagata-like antigenic properties evolved from one lineage and became clearly antigenically distinct from the Victoria-like viruses around 1988. Recombinant mutant viruses were tested to show that insertions and deletions (indels), as observed frequently in influenza B virus hemagglutinin, had little effect on antigenic properties. In contrast, amino-acid substitutions at positions 148, 149, 150, and 203, adjacent to the hemagglutinin receptor binding site, determined the main antigenic differences between the Early, Victoria-like, and Yamagata-like viruses. Surprisingly, substitutions at two of the four positions reverted in recent viruses of the Victoria lineage, resulting in antigenic properties similar to viruses circulating ∼50 y earlier. These data shed light on the antigenic diversification of influenza viruses and suggest there may be limits to the antigenic evolution of influenza B virus.

Effect of 2020-21 and 2021-22 Highly Pathogenic Avian Influenza H5 Epidemics on Wild Birds, the Netherlands.

The number of highly pathogenic avian influenza (HPAI) H5-related infections and deaths of wild birds in Europe was high during October 1, 2020-September 30, 2022. To quantify deaths among wild species groups with known susceptibility for HPAI H5 during those epidemics, we collected and recorded mortality data of wild birds in the Netherlands. HPAI virus infection was reported in 51 bird species. The species with the highest numbers of reported dead and infected birds varied per epidemic year; in 2020-21, they were within the Anatidae family, in particular barnacle geese (Branta leucopsis) and in 2021-22, they were within the sea bird group, particularly Sandwich terns (Thalasseus sandvicensis) and northern gannet (Morus bassanus). Because of the difficulty of anticipating and modeling the future trends of HPAI among wild birds, we recommend monitoring live and dead wild birds as a tool for surveillance of the changing dynamics of HPAI.

4-Octyl itaconate reduces influenza A replication by targeting the nuclear export protein CRM1.

IMPORTANCE: Itaconate derivates, as well as the naturally produced metabolite, have been proposed as antivirals against influenza virus. Here, the mechanism behind the antiviral effects of exogenous 4-octyl itaconate (4-OI), a derivative of itaconate, against the influenza A virus replication is demonstrated. The data indicate that 4-OI targets the cysteine at position 528 of the CRM1 protein, resulting in inhibition of the nuclear export of viral ribonucleoprotein complexes in a similar manner as previously described for other selective inhibitors of nuclear export. These results postulate a mechanism not observed before for this immuno-metabolite derivative. This knowledge is helpful for the development of derivatives of 4-OI as potential antiviral and anti-inflammatory therapeutics.

Field Research Is Essential to Counter Virological Threats.

The interface between humans and wildlife is changing and, with it, the potential for pathogen introduction into humans has increased. Avian influenza is a prominent example, with an ongoing outbreak showing the unprecedented expansion of both geographic and host ranges. Research in the field is essential to understand this and other zoonotic threats. Only by monitoring dynamic viral populations and defining their biology in situ can we gather the information needed to ensure effective pandemic preparation.

Author Correction: Ancient hepatitis B viruses from the Bronze Age to the Medieval period.

In Fig. 2 of this Letter, the 'E' and 'G' clade labels were inadvertently reversed, and in Table 2 the genotype of DA27 was 'D1' instead of 'D5'. These have been corrected online.

Ancient hepatitis B viruses from the Bronze Age to the Medieval period.

Hepatitis B virus (HBV) is a major cause of human hepatitis. There is considerable uncertainty about the timescale of its evolution and its association with humans. Here we present 12 full or partial ancient HBV genomes that are between approximately 0.8 and 4.5 thousand years old. The ancient sequences group either within or in a sister relationship with extant human or other ape HBV clades. Generally, the genome properties follow those of modern HBV. The root of the HBV tree is projected to between 8.6 and 20.9 thousand years ago, and we estimate a substitution rate of 8.04 × 10-6-1.51 × 10-5 nucleotide substitutions per site per year. In several cases, the geographical locations of the ancient genotypes do not match present-day distributions. Genotypes that today are typical of Africa and Asia, and a subgenotype from India, are shown to have an early Eurasian presence. The geographical and temporal patterns that we observe in ancient and modern HBV genotypes are compatible with well-documented human migrations during the Bronze and Iron Ages1,2. We provide evidence for the creation of HBV genotype A via recombination, and for a long-term association of modern HBV genotypes with humans, including the discovery of a human genotype that is now extinct. These data expose a complexity of HBV evolution that is not evident when considering modern sequences alone.

Contemporary human H3N2 influenza A viruses require a low threshold of suitable glycan receptors for efficient infection.

Recent human H3N2 influenza A viruses have evolved to employ elongated glycans terminating in α2,6-linked sialic acid as their receptors. These glycans are displayed in low abundancies by (humanized) Madin-Darby Canine Kidney cells, which are commonly employed to propagate influenza A virus, resulting in low or no viral propagation. Here, we examined whether the overexpression of the glycosyltransferases ß-1,3-N-acetylglucosaminyltransferase and ß-1,4-galactosyltransferase 1, which are responsible for the elongation of poly-N-acetyllactosamines (LacNAcs), would result in improved A/H3N2 propagation. Stable overexpression of ß-1,3-N-acetylglucosaminyltransferase and ß-1,4-galactosyltransferase 1 in Madin-Darby Canine Kidney and "humanized" Madin-Darby Canine Kidney cells was achieved by lentiviral integration and subsequent antibiotic selection and confirmed by qPCR and protein mass spectrometry experiments. Flow cytometry and glycan mass spectrometry experiments using the ß-1,3-N-acetylglucosaminyltransferase and/or ß-1,4-galactosyltransferase 1 knock-in cells demonstrated increased binding of viral hemagglutinins and the presence of a larger number of LacNAc repeating units, especially on "humanized" Madin-Darby Canine Kidney-ß-1,3-N-acetylglucosaminyltransferase cells. An increase in the number of glycan receptors did, however, not result in a greater infection efficiency of recent human H3N2 viruses. Based on these results, we propose that H3N2 influenza A viruses require a low number of suitable glycan receptors to infect cells and that an increase in the glycan receptor display above this threshold does not result in improved infection efficiency.

US National Institutes of Health Prioritization of SARS-CoV-2 Variants.

Since late 2020, SARS-CoV-2 variants have regularly emerged with competitive and phenotypic differences from previously circulating strains, sometimes with the potential to escape from immunity produced by prior exposure and infection. The Early Detection group is one of the constituent groups of the US National Institutes of Health National Institute of Allergy and Infectious Diseases SARS-CoV-2 Assessment of Viral Evolution program. The group uses bioinformatic methods to monitor the emergence, spread, and potential phenotypic properties of emerging and circulating strains to identify the most relevant variants for experimental groups within the program to phenotypically characterize. Since April 2021, the group has prioritized variants monthly. Prioritization successes include rapidly identifying most major variants of SARS-CoV-2 and providing experimental groups within the National Institutes of Health program easy access to regularly updated information on the recent evolution and epidemiology of SARS-CoV-2 that can be used to guide phenotypic investigations.

Continued adaptation of A/H2N2 viruses during pandemic circulation in humans.

Influenza A viruses of the H2N2 subtype sparked a pandemic in 1957 and circulated in humans until 1968. Because A/H2N2 viruses still circulate in wild birds worldwide and human population immunity is low, the transmissibility of six avian A/H2N2 viruses was investigated in the ferret model. None of the avian A/H2N2 viruses was transmitted between ferrets, suggesting that their pandemic risk may be low. The transmissibility, receptor binding preference and haemagglutinin (HA) stability of human A/H2N2 viruses were also investigated. Human A/H2N2 viruses from 1957 and 1958 bound to human-type α2,6-linked sialic acid receptors, but the 1958 virus had a more stable HA, indicating adaptation to replication and spread in the new host. This increased stability was caused by a previously unknown stability substitution G205S in the 1958 H2N2 HA, which became fixed in A/H2N2 viruses after 1958. Although individual substitutions were identified that affected the HA receptor binding and stability properties, they were not found to have a substantial effect on transmissibility of A/H2N2 viruses via the air in the ferret model. Our data demonstrate that A/H2N2 viruses continued to adapt during the first year of pandemic circulation in humans, similar to what was previously shown for the A/H1N1pdm09 virus.

Strong host phylogenetic and ecological effects on host competency for avian influenza in Australian wild birds.

Host susceptibility to parasites is mediated by intrinsic and external factors such as genetics, ecology, age and season. While waterfowl are considered central to the reservoir community for low pathogenic avian influenza A viruses (LPAIV), the role of host phylogeny has received limited formal attention. Herein, we analysed 12 339 oropharyngeal and cloacal swabs and 10 826 serum samples collected over 11 years from wild birds in Australia. As well as describing age and species-level differences in prevalence and seroprevalence, we reveal that host phylogeny is a key driver in host range. Seasonality effects appear less pronounced than in the Northern Hemisphere, while annual variations are potentially linked to El Niño-Southern Oscillation. Our study provides a uniquely detailed insight into the evolutionary ecology of LPAIV in its avian reservoir community, defining distinctive processes on the continent of Australia and expanding our understanding of LPAIV globally.

Insertions of codons encoding basic amino acids in H7 hemagglutinins of influenza A viruses occur by recombination with RNA at hotspots near snoRNA binding sites.

The presence of multiple basic amino acids in the protease cleavage site of the hemagglutinin (HA) protein is the main molecular determinant of virulence of highly pathogenic avian influenza (HPAI) viruses. Recombination of HA RNA with other RNA molecules of host or virus origin is a dominant mechanism of multibasic cleavage site (MBCS) acquisition for H7 subtype HA. Using alignments of HA RNA sequences from documented cases of MBCS insertion due to recombination, we show that such recombination with host RNAs is most likely to occur at particular hotspots in ribosomal RNAs (rRNAs), transfer RNAs (tRNAs), and viral RNAs. The locations of these hotspots in highly abundant RNAs indicate that RNA recombination is facilitated by the binding of small nucleolar RNA (snoRNA) near the recombination points.

Long-Term Protective Effect of Serial Infections with H5N8 Highly Pathogenic Avian Influenza Virus in Wild Ducks.

Highly pathogenic avian influenza viruses (HPAIVs) of the Goose/Guangdong (Gs/Gd) lineage are an emerging threat to wild birds. In the 2016-2017 H5N8 outbreak, unexplained variability was observed in susceptible species, with some reports of infected birds dying in high numbers and other reports of apparently subclinical infections. This experimental study was devised to test the hypothesis that previous infection with a less-virulent HPAIV (i.e., 2014 H5N8) provides long-term immunity against subsequent infection with a more-virulent HPAIV (i.e., 2016 H5N8). Therefore, two species of wild ducks-the more-susceptible tufted duck (Aythya fuligula) and the more-resistant mallard (Anas platyrhynchos)-were serially inoculated, first with 2014 H5N8 and after 9 months with 2016 H5N8. For both species, a control group of birds was first sham inoculated and after 9 months inoculated with 2016 H5N8. Subsequent infection with the more-virulent 2016 H5N8 caused no clinical signs in tufted ducks that had previously been infected with 2014 H5N8 (n = 6) but caused one death in tufted ducks that had been sham inoculated (n = 7). In mallards, 2016 H5N8 infection caused significant body weight loss in previously sham-inoculated birds (n = 8) but not in previously infected birds (n = 7). IMPORTANCE This study showed that ducks infected with a less-virulent HPAIV developed immunity that was protective against a subsequent infection with a more-virulent HPAIV 9 months later. Following 2014 H5N8 infection, the proportion of birds with detectable influenza nucleoprotein antibody declined from 100% (8/8) in tufted ducks and 78% (7/9) in mallards after 1 month to 33% (2/6) in tufted ducks and 29% (2/7) in mallards after 9 months. This finding helps predict the expected impact that an HPAIV outbreak may have on wild bird populations, depending on whether they are immunologically naive or have survived previous infection with HPAIV.

Contribution of Neuraminidase to the Efficacy of Seasonal Split Influenza Vaccines in the Ferret Model.

Seasonal influenza vaccination takes into account primarily hemagglutinin (HA)-specific neutralizing antibody responses. However, the accumulation of substitutions in the antigenic regions of HA (i.e., antigenic drift) occasionally results in a mismatch between the vaccine and circulating strains. To prevent poor vaccine performance, we investigated whether an antigenically matched neuraminidase (NA) may compensate for reduced vaccine efficacy due to a mismatched HA. Ferrets were vaccinated twice with adjuvanted split inactivated influenza vaccines containing homologous HA and NA (vacH3N2), only homologous HA (vacH3N1), only homologous NA (vacH1N2), heterologous HA and NA (vacH1N1), or phosphate-buffered saline (vacPBS), followed by challenge with H3N2 virus (A/Netherlands/16190/1968). Ferrets vaccinated with homologous HA (vacH3N2 and vacH3N1) displayed minimum fever and weight loss compared to vacH1N1 and vacPBS ferrets, while ferrets vaccinated with NA-matched vacH1N2 displayed intermediate fever and weight loss. Vaccination with vacH1N2 further led to a reduction in virus shedding from the nose and undetectable virus titers in the lower respiratory tract, similarly to when the homologous vacH3N2 was used. Some protection was observed upon vacH1N1 vaccination, but this was not comparable to that observed for vacH1N2, again highlighting the important role of NA in vaccine-induced protection. These results illustrate that NA antibodies can prevent severe disease caused by influenza virus infection and that an antigenically matched NA in seasonal vaccines might prevent lower respiratory tract complications. This underlines the importance of considering NA during the yearly vaccine strain selection process, which may be particularly beneficial in seasons when the HA component of the vaccine is mismatched. IMPORTANCE Despite the availability of vaccines, influenza virus infections continue to cause substantial morbidity and mortality in humans. Currently available influenza vaccines take primarily the hemagglutinin (HA) into account, but the highly variable nature of this protein as a result of antigenic drift has led to a recurrent decline in vaccine effectiveness. While the protective effect of neuraminidase (NA) antibodies has been highlighted by several studies, there are no requirements with regard to quantity or quality of NA in licensed vaccines, and NA immunity remains largely unexploited. Since antigenic changes in HA and NA are thought to occur asynchronously, NA immunity could compensate for reduced vaccine efficacy when drift in HA occurs. By matching and mismatching the HA and NA components of monovalent split inactivated vaccines, we demonstrated the potential of NA immunity to protect against disease, virus replication in the lower respiratory tract, and virus shedding in the ferret model.

Mass Mortality Caused by Highly Pathogenic Influenza A(H5N1) Virus in Sandwich Terns, the Netherlands, 2022.

We collected data on mass mortality in Sandwich terns (Thalasseus sandvicensis) during the 2022 breeding season in the Netherlands. Mortality was associated with at least 2 variants of highly pathogenic avian influenza A(H5N1) virus clade 2.3.4.4b. We report on carcass removal efforts relative to survival in colonies. Mitigation strategies urgently require structured research.

HIGHLY PATHOGENIC AVIAN INFLUENZA VIRUS (H5N8) OUTBREAK IN A WILD BIRD RESCUE CENTER, THE NETHERLANDS: CONSEQUENCES AND RECOMMENDATIONS.

Since the emergence of the Goose/Guangdong H5 lineage in 1996 and spillover of highly pathogenic avian influenza (HPAI) from poultry to wild birds, outbreaks have become increasingly frequent in wild birds. The latest outbreak in the Netherlands occurred in the fall-winter of 2020-2021 and was linked to incursions of HPAI H5N8 virus. During the larger national outbreak, wild birds in rehabilitation center "Vogelklas Karel Schot (VKS)" in Rotterdam presented with clinical signs compatible with HPAI, including head shaking, torticollis, and abnormal gait. During an epidemiologic investigation at VKS, water samples from the pools in the enclosures and oropharyngeal and cloacal swabs from 128 birds of different species were analyzed for the presence of H5N8 virus. Forty-five birds and the pool water tested positive for the virus. The outbreak at VKS was likely introduced by one or more infected geese (Anser anser, Anser anser domesticus, Branta leucopsis), after which the virus spread via pool water and with the relocation of infected birds within the center. In principle, such outbreaks are preventable. Recent updates about HPAI to provide guidance to help avoid future incursions of HPAI into wildlife rescue centers are reported.

Outbreak Severity of Highly Pathogenic Avian Influenza A(H5N8) Viruses Is Inversely Correlated to Polymerase Complex Activity and Interferon Induction.

Highly pathogenic avian influenza A(H5N8) viruses first emerged in China in 2010 and in 2014 spread throughout Asia and to Europe and the United States via migrating birds. Influenza A(H5N8) viruses were first detected in the Netherlands in 2014 and caused five outbreaks in poultry farms but were infrequently detected in wild birds. In 2016, influenza A(H5N8) viruses were reintroduced into the Netherlands, resulting in eight poultry farm outbreaks. This outbreak resulted in numerous dead wild birds with severe pathology. Phylogenetic analysis showed that the polymerase genes of these viruses had undergone extensive reassortment between outbreaks. Here, we investigated the differences in virulence between the 2014-15 and the 2016-17 outbreaks by characterizing the polymerase complex of influenza A(H5N8) viruses from both outbreaks. We found that viruses from the 2014-15 outbreak had significantly higher polymerase complex activity in both human and avian cell lines than did those from the 2016-17 outbreak. No apparent differences in the balance between transcription and replication of the viral genome were observed. Interestingly, the 2014-15 polymerase complexes induced significantly higher levels of interferon beta (IFN-ß) than the polymerase complexes of the 2016-17 outbreak viruses, mediated via retinoic acid-inducible gene I (RIG-I). Inoculation of primary duck cells with recombinant influenza A(H5N8) viruses, including viruses with reassorted polymerase complexes, showed that the polymerase complexes from the 2014-15 outbreak induced higher levels of IFN-ß despite relatively minor differences in replication capacity. Together, these data suggest that despite the lower levels of polymerase activity, the higher 2016-17 influenza A(H5N8) virus virulence may be attributed to the lower level of activation of the innate immune system.IMPORTANCE Compared to the 2014-15 outbreak, the 2016-17 outbreak of influenza A(H5N8) viruses in the Netherlands and Europe was more virulent; the number of dead or diseased wild birds found and the severity of pathological changes were higher during the 2016-17 outbreak. The polymerase complex plays an important role in influenza virus virulence, and the gene segments of influenza A(H5N8) viruses reassorted extensively between the outbreaks. In this study, the 2014-15 polymerase complexes were found to be more active, which is counterintuitive with the observed higher virulence of the 2016-17 outbreak viruses. Interestingly, the 2014-15 polymerase complexes also induced higher levels of IFN-ß. These findings suggest that the higher virulence of influenza A(H5N8) viruses from the 2016-17 outbreak may be related to the lower induction of IFN-ß. An attenuated interferon response could lead to increased dissemination, pathology, and mortality, as observed in (wild) birds infected during the 2016-2017 outbreak.

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

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.