# Genetic Evolution of Avian Influenza A (H9N2) Viruses Isolated from Domestic Poultry in Uganda Reveals Evidence of Mammalian Host Adaptation, Increased Virulence and Reduced Sensitivity to Baloxavir

Gladys Atim, Titus Tugume, Qouilazoni A Ukuli, Bernard Erima, Andrew Mubiru, Hannah Kibuuka, Edison Mworozi, Pamela Mckenzie, Jasmine C M Turner, David Walker, et al.

► To cite this version:

Gladys Atim, Titus Tugume, Qouilazoni A Ukuli, Bernard Erima, Andrew Mubiru, et al. Genetic Evolution of Avian Influenza A (H9N2) Viruses Isolated from Domestic Poultry in Uganda Reveals Evidence of Mammalian Host Adaptation, Increased Virulence and Reduced Sensitivity to Baloxavir. Viruses, 2022, 14 (9), pp.2074. 10.3390/v14092074. hal-04163674

HAL Id: hal-04163674

https://hal.inrae.fr/hal-04163674

Submitted on 17 Jul 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
Article

# Genetic Evolution of Avian Influenza A (H9N2) Viruses Isolated from Domestic Poultry in Uganda Reveals Evidence of Mammalian Host Adaptation, Increased Virulence and Reduced Sensitivity to Baloxavir

Gladys Atim<fnref n="1" /> ID, Titus Tugume<fnref n="1" />, Qouilazoni A. Ukuli<fnref n="1" />, Bernard Erima<fnref n="1" />, Andrew Mubiru<fnref n="1" />, Hannah Kibuuka<fnref n="1" />, Edison Mworozi<fnref n="1" />, Pamela McKenzie<fnref n="2" />, Jasmine C. M. Turner<fnref n="2" /> ID, David Walker<fnref n="2" />, Trushar Jeevan<fnref n="2" />, Robert G. Webster<fnref n="2" />, Jeremy Jones<fnref n="2" />, Richard J. Webby<fnref n="2" /> ID, Mariette F. Ducatez<fnref n="3" /> ID, Fred Wabwire-Mangen<fnref n="1" />,<fnref n="4" /> and Denis K. Byarugaba<fnref n="1" />,<fnref n="5" />,* ID

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

check for
updates

**Citation:** Atim, G.; Tugume, T.; Ukuli, Q.A.; Erima, B.; Mubiru, A.; Kibuuka, H.; Mworozi, E.; McKenzie, P.; Turner, J.C.M.; Walker, D.; et al. Genetic Evolution of Avian Influenza A (H9N2) Viruses Isolated from Domestic Poultry in Uganda Reveals Evidence of Mammalian Host Adaptation, Increased Virulence and Reduced Sensitivity to Baloxavir. *Viruses* **2022**, *14*, 2074. https://doi.org/10.3390/v14092074

Academic Editor: Hualan Chen

Received: 25 August 2022

Accepted: 14 September 2022

Published: 18 September 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** A (H9N2) avian influenza A viruses were first detected in Uganda in 2017 and have since established themselves in live bird markets. The aim of this study was to establish the subsequent genetic evolution of H9N2 viruses in Uganda. Cloacal samples collected from live bird market stalls in Kampala from 2017 to 2019 were screened by RT-PCR for influenza A virus and H9N2 viruses were isolated in embryonated eggs. One hundred and fifty H9N2 isolates were subjected to whole genome sequencing on the Illumina MiSeq platform. The sequence data analysis and comparison with contemporary isolates revealed that the virus was first introduced into Uganda in 2014 from ancestors in the Middle East. There has since been an increase in nucleotide substitutions and reassortments among the viruses within and between live bird markets, leading to variations in phylogeny of the different segments, although overall diversity remained low. The isolates had several mutations such as HA-Q226L and NS-I106M that enable mammalian host adaptation, NP-M105V, PB1-D3V, and M1-T215A known for increased virulence/pathogenicity and replication, and PA-E199D, NS-P42S, and M2-S31N that promote drug resistance. The PA-E199D substitution in particular confers resistance to the endonuclease inhibitor Baloxavir acid, which is one of the new anti-influenza drugs. Higher EC50 was observed in isolates with a double F105L+E199D substitution that may suggest a possible synergistic effect. These H9N2 viruses have established an endemic situation in live bird markets in Uganda because of poor biosecurity practices and therefore pose a zoonotic threat. Regular surveillance is necessary to further generate the needed evidence for effective control strategies and to minimize the threats.

**Keywords:** influenza; phylogenetics; molecular markers; reassortment

## 1. Introduction

Avian influenza viruses (AIVs) are negative-sense RNA viruses of the *Orthomyxoviridae* family whose natural hosts are birds [1]. They constantly undergo evolution mediated through mutation and reassortment resulting in viruses with varying pathogenicity, virulence, host specificity, antigenicity, and drug susceptibility [2]. The incorporation of mutations is driven by the influenza RNA polymerase which lacks proofreading functions [3]. Reassortment is a result of the segmented nature of the genome and mixing of genes following coinfections [3–5].

<footnote n="1">¹ Makerere University Walter Reed Project, Kampala P.O. Box 16524, Uganda</footnote>

<footnote n="2">² St Jude Children's Research Hospital, MS330, Memphis, TN 38105, USA</footnote>

<footnote n="3">³ Interactions Hôtes-Agents Pathogène, Ecole Nationale Vétérinaire de Toulouse, 31076 Toulouse, France</footnote>

<footnote n="4">⁴ School of Public Health, College of Health Sciences, Makerere University, Kampala P.O. Box 7062, Uganda</footnote>

<footnote n="5">⁵ College of Veterinary Medicine, Makerere University, Kampala P.O. Box 7062, Uganda</footnote>

<footnote>* Correspondence: denis.byarugaba@mak.ac ug</footnote>
Of the diversity of viruses in avian hosts, A(H9N2) viruses are of particular interest. While being of low pathogenicity for birds, they have donated gene segments to novel influenza strains, such as H5N1, that appear to support infection of mammalian hosts including humans [6–11].

First isolated in the United States in 1966 [12], H9N2 viruses have been detected in many regions of the world and those belonging to the G1 lineage have been endemic in the Middle East (ME) since the late 1990s and in North Africa since the late 2000s. They were not, however, documented in sub-Saharan Africa until 2017, despite surveillance activities [13,14].

In Uganda, H9N2 viruses were first detected in poultry in 2017 and have since become established in live bird markets (LBMs) in part due to poor biosecurity [15]. A recent study found that H9N2 viruses from Uganda had reassorted more than those from West Africa [16], warranting further investigation.

Although they only cause mild to moderate illness in birds, H9N2 viruses impart a substantial economic burden by reducing productivity and performance of affected flocks. In addition, as of June 2022, there have been 102 laboratory-confirmed cases of human infection with H9N2 and 2 recorded deaths [17] in various countries including China [10,18], India [19,20], and Egypt [13], and the same in Senegal as of 2021; the human virus from Senegal was genetically similar to Ugandan H9N2 viruses [21]. In almost all reported human cases, previous contact with poultry was identified as a risk factor so H9N2 viruses remain a constant zoonotic threat [1,22]. In this study, we investigated H9N2 virus evolution in Uganda, identified further reassortment events and detected several mutations of interest.

## 2. Materials and Methods

### 2.1. Sample Collection

Cloacal-oropharyngeal samples were collected from birds in LBMs around Kampala city (Kalerwe, Kasubi, Nakasero and Nakawa markets which receive birds from all over the country) from 2017 to 2019. Samples were collected once a month on the same day, from up to 20 birds that were randomly selected from multiple stalls in each of the four markets. A total of 80 samples were collected each month. Samples were collected using Dacron swabs and stored in antibiotic enriched virus transport media.

RNA was extracted from swabs by using a QIAamp Viral RNA extraction kit (Qiagen, Germany) according to the manufacturer's instructions, and they were tested for Influenza A by real time RT-PCR as previously described [16]. Positive samples were further tested for various subtypes including H9 subtype by RT-PCR.

Samples that were H9 positive were inoculated into 9-day old embryonated chicken eggs. These were incubated at 37 °C for at least 72 h. The eggs were then refrigerated for at least 4 h before allantoic fluid was harvested.

### 2.2. Whole Genome Sequencing and Sequence Identification

Following a confirmatory PCR for H9, isolates with a CT value below 30 were selected for whole genome sequencing; in total, there were 150. RNA was extracted and all eight gene segments were amplified in an RT-PCR with SuperScript III One-Step RT-PCR System with Platinum® Taq High Fidelity DNA Polymerase (Invitrogen, Waltham, MA, USA). Duplicate multiplexed paired-end sequencing libraries were generated using an Illumina Nextera XT library prep kit (FC-131-1096). The libraries were then fragmented and tagged with Illumina MiSeq-compatible barcode adapters using the NEBNext Ultra DNA Library Preparation kit (New England Biolabs, Ipswich, MA, USA) as previously described [16]. Paired-end sequencing was conducted on the Illumina MiSeq platform. Sequences were deposited into GenBank under accession numbers MW205009 to MW205198, MZ542846 to MZ543305, OP388450 to OP388679 and OP388789 to OP388827.
## 2.3. Sequence Analysis

Consensus sequences for each segment were generated. A BLAST search was conducted to find the most related sequences by E-score. In order to determine whether the isolates carried any specific Sequence Features (SF) of importance, SF variant type analysis was performed in the Influenza Research Database (IRD) [23].

Groups of sequences were aligned using the MUSCLE alignment tool offered in the MEGA11 software along with calculations of maximum composite likelihood, pairwise and overall mean distances [24]. Phylogenetic analysis was performed using BEAST v1.10.4. This was carried out with a relaxed (uncorrelated lognormal) molecular clock for a constant population size coalescent tree prior. Other pre-set conditions were a Hasegawa–Kishino–Yano (HYK) + Gamma (4 categories) nucleotide substitution model, a consideration of 3 codon positions and an unlinked base frequency across all codon positions. The length of Markov chain Monte Carlo (MCMC) was set at 100 million states.

The output was assessed in Tracer v1.7.1 to determine various evolution parameters such as the annual mean rate of nucleotide substitution per site in the segments, their time to most common recent ancestor (tMRCA) as well as the number of lineages through time. Maximum Clade Credibility (MCC) trees with a 10% burn-in were generated using the TreeAnnotator v1.10.4 and visualized using FigTree v1.4.4 [25].

## 2.4. Assessment of H9N2 Susceptibility to Baloxivir (BXA)

Full plaque reduction EC50 analysis was performed on randomly selected PA-E199D H9N2 viruses to determine susceptibility to BXA. MDCK cells ($1 \times 10^6$) were plated in six-well plates, inoculated with influenza viruses (multiplicity of infection [MOI] of 0.001), cultured in MEM containing 1% bovine serum albumin (BSA), 1 µg/mL TPCK trypsin, and BXA (0, 1, 10, or 100 nM). At 1 h post-virus inoculation, cells were washed, overlaid with MEM containing 0.45% immunodiffusion-grade agarose (MP Biomedicals, Santa Ana, SC, USA), 1% BSA, 1.8 µg/mL TPCK trypsin, and BXA (0.001 nM to 500 nM). At 72 hpi to 120 hpi, the overlays were removed, and the cell monolayers were stained with 1% crystal violet–10% formaldehyde. The number of PFUs per well was calculated, and EC50s were determined by using the log (inhibitor) versus response logistic nonlinear regression equation in GraphPad Prism 8.0 software.

## 2.5. Hemagglutination Inhibition (HI) Assay

Antigenic analysis was also carried out on the selected viruses with additional antigens and serum samples that were provided by members of the World Health Organization Global Influenza Surveillance and Response System. A total of 25 µL of PBS was dispensed into each well of a 96-well microtiter plate. A total of 25 µL of reference anti-serum was added to the first well the plate and two-fold dilutions made. Four hemagglutinating units of virus was added to each well and left for a minimum of 30 min at room temperature, then 0.5% ($v/v$) turkey RBCs were added to each well and allowed to settle for 40 min at room temperature. The HI titer recorded was the highest dilution of serum causing complete inhibition of hemagglutination [26].

# 3. Results

## 3.1. H9N2 Detection

A total of 2670 samples were collected between 2017 and 2019. The overall prevalence of H9N2 virus was 49%, with a steady increase from 47% in 2017 to 61% in 2019. Positivity rates in individual markets ranged from 33% to 68% [16].

## 3.2. Phylogenetic Analyses

In total, 150 isolates were selected for whole genome sequencing. Of these, 48 were from Nakawa market, 42 from Nakasero market, 25 from Kalerwe market, and 35 from Kasubi market. All were from chicken and belonged to the G1 lineage. Phylogenetic analysis of the HA segment showed that the Uganda H9N2 viruses clustered in the same
genetic clade as viruses from Kenya, (Figure 1) and were similarly related to viruses from
Saudi Arabia, the United Arab Emirates and those from Western and Northern Africa:
Algeria, Benin, Burkina Faso, Ghana, Libya, Morocco, Senegal, Togo and Tunisia.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Phylogenetic tree for HA. Ugandan viruses are highlighted in blue, Kenyan in red, North and West African countries in pink, and the Middle East in the green clade. Posterior probabilities are shown on the branches. The top scale represents the estimated number of substitutions and the bottom scale shows time of divergence. HKY + G (4) nucleotide substitution model and an MCMC of 100 million states was applied in BEAST to generate this tree.</figcaption></figure>

The other gene segments had similar patterns of relatedness as the HA as shown
in Figures A1–A7 (Appendix A) of other segments, respectively. The segments had very
similar trees with no major differences in the sequences they were most closely related to.
There were, however, minor differences in clade arrangement and general tree topology.
### 3.3. Molecular Clock Analysis

The tMRCA of the Uganda H9N2 HA gene as calculated by BEAST was March 2014 (interval from 2013 to 2015 for 95% of the samples). The tMRCA range for all the segments also fell within the same time range of 2013 to 2014. It was, however, difficult to pinpoint a true ancestor as the most closely related sequences (from A/Pheasant/United_Arab_Emirates/D1307/2011) were quite distant from the Uganda clade in terms of branch length and isolation date. This suggests a gap in surveillance over this time period.

The genetic diversity of the Ugandan viruses calculated by maximum composite likelihood ranged from 0.02% to 0.1%. The annual mean rate of substitutions per site was higher overall in the Ugandan virus subtree than in the rest of the H9N2 viruses with the HA segment having the most substitutions, while the MP had the lowest under the uncorrelated lognormal clock (Figure 2). When plotted against time, there was evidence of an increasing substitution rate as shown for the HA segment (Figure 3).

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Comparing the mean rate of substitutions of the segments in Ugandan viruses and the rest. Ugandan viruses had a higher-than-average substitution rate.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Mean rate of nucleotide substitution per site in the HA segment of Uganda H9N2 viruses over time. Each dot represents a branch in the tree and the age(root) corresponds to time of divergence.</figcaption></figure>

3.4. Reassortment

Pairwise distances between the Ugandan H9N2 viruses gave rise to multiple clusters in different gene segments; however, the specific clustering differed from one segment to another and did not segregate by LBM of origin. This was also evident in the tree topology with multiple sub-clades within the Ugandan clusters that differed from segment to segment. Together, these findings are consistent with reassortment between H9N2 viruses within Uganda. To investigate this possibility, all clades containing branches of lengths < 1 were collapsed, and a representative taxon (the oldest) was picked from each clade. The clades in HA were used as points of reference and color-coded in creating Figure 4, to visualize reassortment events. All twelve representative isolates had unique constellations (different patterns). Only two isolates, A/Chicken/Uganda/MUWRP/713/2017 and A/Chicken/Uganda/MUWRP/756/2017, showed no reassortment at all. Besides HA (the reference segment), the PB2 segment had the most reassortment, with the greatest number of unique clades, while the NP segment had the least reassortment, with only four clades, all from the same market. This complexity is consistent with introduction of genetic diversity from different markets. These findings suggest that there is a significant level of reassortment within and between the different LBMs in our study.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Reassortment within Ugandan H9N2 viruses. Clades were created in each segment by collapsing all branches with length < 1. From these, the oldest virus isolate was selected and those from the HA segment were given a color code. Representative isolates were cross-referenced to see if they fell within the same clades across all segments. Symbols for markets of origin were also assigned and overlaid.</figcaption></figure>

### 3.5. *Molecular Marker Analysis*

A detailed examination of sequence features within the Ugandan H9N2 virus sequences revealed multiple substitutions in the H9N2 viruses that are associated with increased virulence in both mammals and avian species, increased virus replication, and reduced sensitivity to antiviral drugs. Details of the relevant substitutions in the different segments found in each segment are presented in Table 1.

Among the sequence feature variants identified were those that have been associated with infection in mammals such as Q226L in HA, which increases binding to sialic acid receptors preferred by mammalian-adapted viruses [27]. Others that increase virulence in mammals were N30D and T215A in M1 [28] and NS-P42S [29].

Multiple mutations that are known to increase replication via enhanced polymerase activity were also found and included NP-E210D, S37A and N383D in PA, D3V and D622G in PB1, D9N and K226R in PB2 [30–35].

Sequence markers associated with increased resistance to amantadine were also present in M2 (S31N) [36]. Eight H9N2 viruses were found to have a PA-E199D mutation, which is a potential marker associated with reduced sensitivity to Baloxivir (BXA) [37]. Wild-type (WT) E199 and E199D viruses were tested for their susceptibility to BXA. An almost three-fold reduction in the median EC50 value was observed for H9N2 viruses with the E199D mutation (0.7–1.2 nM for E119D compared to 0.1–0.5 nM for WT) (Figure 5).
<table><caption>**Table 1.** Sequence Feature Variants in H9N2 virus proteins and their phenotypic function.</caption><thead><tr><th>Protein</th><th>Sequence Features</th><th>Phenotype</th><th>Frequency</th></tr></thead><tbody><tr><td rowspan="2">HA</td><td>T197, G520</td><td>Acquisition of respiratory droplet transmission and human-like clinical symptoms</td><td rowspan="2">All</td></tr><tr><td>Q226L</td><td>Increased virus binding to α2-6, enhanced replication in mammalian cells and ferrets, enhanced contact transmission in ferrets</td></tr><tr><td>M1</td><td>N30D, I43M, T215A</td><td>Increased virulence in mice</td><td rowspan="2">S31N absent in 1</td></tr><tr><td>M2</td><td>S31N</td><td>Increased Amantadine resistance</td></tr><tr><td rowspan="2">NP</td><td>M105V, A184K</td><td>Increased virulence in chicken [38,39]</td><td rowspan="2">All</td></tr><tr><td>E210D</td><td>Increased polymerase activity in mammalian cell line</td></tr><tr><td rowspan="2">NS</td><td>P42S</td><td>Increased virulence and decreased antiviral response in mice</td><td rowspan="2">All</td></tr><tr><td>I106M</td><td>Increased viral replication in mammalian cells and virulence in mice [40]</td></tr><tr><td rowspan="2">PA</td><td>S37A, N383D</td><td>Increased polymerase activity in mammalian and avian cell line</td><td rowspan="2">All<br>E199D in 8</td></tr><tr><td>E199D</td><td>Decreased sensitivity to Baloxivir</td></tr><tr><td rowspan="2">PB1</td><td>D3V</td><td>Increased polymerase activity and replication in mammalian and avian cell lines</td><td rowspan="2">All</td></tr><tr><td>D622G</td><td>Increased polymerase activity in mammalian cell lines</td></tr><tr><td rowspan="3">PB2</td><td>L89V, G309D, T339K, R477G, I495V, A676V</td><td>Increased polymerase activity in mammalian cell line and increased virulence in mice</td><td rowspan="3">All</td></tr><tr><td>D9N</td><td>Increased virulence in mice</td></tr><tr><td>K256R</td><td>Increased polymerase activity in mammalian cell lines</td></tr></tbody></table>

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

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

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

<figure><img src="image_11.png" /><figcaption>**Figure 5.** *Cont.*</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>F. A/Chicken/MUWRP-Uganda/1064/2019 (PA E199D)</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>G. A/Chicken/MUWRP-Uganda/1082/2019 (PA E199D)</figcaption></figure>

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

**Figure 5.** BXA Susceptibility of H9N2 viruses from Uganda containing the PA E199D substitution. The wild type (WT) viruses in (A–C) have EC50s ranging from 0.2–0.5 nM while those with the ED substitution (E–G) have increased EC50s of 0.74–0.83 nM. (D) shows a virus which had both PA F105L+E199D had the highest EC50 of 1.28nM, which suggests a synergistic effect. Data is representative of an average of three independent dose-response curves.

### 3.6. Antigenic Analysis

Selected viruses were evaluated for their antigenic properties using a panel of reference antigens and antisera. The selected viruses were antigenically similar to each other and reacted within four-fold to post-infection antiserum raised against A/Oman/2747/2019 IDCDC-RG66A, a human H9N2 virus (Table 2).
<table><caption>**Table 2.** Hemagglutination inhibition assay of recent H9N2 viruses from Uganda (0.5% turkey RBC).</caption><thead><tr><th></th><th></th><th></th><th>Clade</th><th>Oman/2747</th><th>Bd/0994</th><th>HK/1073</th><th>HK/33982</th><th>qa/Bd/19462</th><th>ck/Benin</th><th>ck/Bd/46240</th><th>ck/BD/46129</th><th>Collection Date</th><th>Passage History</th></tr></thead><tbody><tr><td colspan="14"><strong>Reference antigen</strong></td></tr><tr><td>1</td><td>A/Oman/2747/2019 IDCDC-RG66A</td><td>G1</td><td>2560</td><td>320</td><td>&lt;</td><td>40</td><td>640</td><td>320</td><td>640</td><td>640</td><td></td><td>V1E2/E1</td></tr><tr><td>2</td><td>A/Bangladesh/0994/2011-IDCDC-RG31</td><td>G1</td><td>1280</td><td>1280</td><td>&lt;</td><td>80</td><td>640</td><td>640</td><td>1280</td><td>2560</td><td></td><td>V1E2/E1</td></tr><tr><td>3</td><td>A/Hong Kong/1073/99</td><td>G1</td><td>&lt;</td><td>&lt;</td><td>160</td><td>640</td><td>&lt;</td><td>&lt;</td><td>&lt;</td><td>&lt;</td><td></td><td>C4</td></tr><tr><td>4</td><td>A/Hong Kong/33982/2009-PR8-IDCDC-RG-26</td><td>G1</td><td>&lt;</td><td>&lt;</td><td>160</td><td>1280</td><td>&lt;</td><td>&lt;</td><td>&lt;</td><td>&lt;</td><td></td><td>V1E3/E2/E1</td></tr><tr><td>5</td><td>A/quail/Bangladesh/19462/2013</td><td>G1</td><td>320</td><td>40</td><td>10</td><td>80</td><td>320</td><td>80</td><td>160</td><td>160</td><td></td><td>E2</td></tr><tr><td>6</td><td>A/chicken/Benin/19-A-01-145-E/2019</td><td>G1</td><td>1280</td><td>320</td><td>&lt;</td><td>80</td><td>640</td><td>640</td><td>640</td><td>1280</td><td></td><td>E1/E2</td></tr><tr><td>7</td><td>A/chicken/Bangladesh/46240/2020</td><td>G1</td><td>80</td><td>80</td><td>&lt;</td><td>40</td><td>80</td><td>160</td><td>640</td><td>1280</td><td></td><td>E1</td></tr><tr><td>8</td><td>A/chicken/Bangladesh/46129/2020</td><td>G1</td><td>640</td><td>640</td><td>&lt;</td><td>80</td><td>640</td><td>320</td><td>1280</td><td>2560</td><td></td><td>E1</td></tr><tr><td colspan="14"><strong>Test antigens</strong></td></tr><tr><td>9</td><td>A/chicken/MUWRP-UGANDA/944/2019</td><td>G1</td><td>1280</td><td>2560</td><td>&lt;</td><td>20</td><td>640</td><td>320</td><td>1280</td><td>2560</td><td>03/13/19</td><td>X/E1</td></tr><tr><td>10</td><td>A/chicken/MUWRP-UGANDA/951/2019</td><td>G1</td><td>2560</td><td>2560</td><td>&lt;</td><td>40</td><td>1280</td><td>640</td><td>1280</td><td>2560</td><td>03/13/19</td><td>X/E1</td></tr><tr><td>11</td><td>A/chicken/MUWRP-UGANDA/957/2019</td><td>G1</td><td>1280</td><td>1280</td><td>&lt;</td><td>20</td><td>320</td><td>320</td><td>1280</td><td>2560</td><td>04/24/19</td><td>X/E1</td></tr><tr><td>12</td><td>A/chicken/MUWRP-UGANDA/977/2019</td><td>G1</td><td>1280</td><td>1280</td><td>&lt;</td><td>20</td><td>640</td><td>320</td><td>1280</td><td>2560</td><td>06/19/19</td><td>X/E1</td></tr><tr><td>13</td><td>A/chicken/MUWRP-UGANDA/1038/2019</td><td>G1</td><td>640</td><td>1280</td><td>&lt;</td><td>20</td><td>320</td><td>320</td><td>640</td><td>1280</td><td>09/25/19</td><td>X/E1</td></tr><tr><td>14</td><td>A/chicken/MUWRP-UGANDA/1064/2019</td><td>G1</td><td>1280</td><td>1280</td><td>10</td><td>20</td><td>640</td><td>320</td><td>1280</td><td>2560</td><td>11/20/19</td><td>X/E1</td></tr><tr><td>15</td><td>A/chicken/MUWRP-UGANDA/1082/2019</td><td>G1</td><td>1280</td><td>1280</td><td>10</td><td>10</td><td>640</td><td>320</td><td>2560</td><td>2560</td><td>12/19/19</td><td>X/E1</td></tr><tr><td colspan="14">Serum production: P = prime</td></tr><tr><td></td><td></td><td></td><td>CDC</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td></td><td></td></tr></tbody></table>

<figure><img src="image_15.png" /></figure>
## 4. Discussion

Non-human Influenza A viruses were first reported in Uganda during a country-wide survey carried out in 2014 [41] with a prevalence of 1.1% in livestock and poultry. Regular surveillance did not reveal any remarkable findings until 2017 when H9N2 viruses were first documented in LBMs. As a consequence of this finding, we conducted monthly surveillance in four LBMs from 2017 to 2019 as previously reported [16]. It is possible that infections are amplified in LBMs due to poor biosecurity and management practices that foster virus maintenance [15]. Ugandan LBMs receive birds from different sources, have insufficient cleaning practices, and have no all-in-all-out policy with new birds introduced into cages where unsold birds are kept. These practices create an environment suitable for AIV maintenance, amplification, and evolution that consequently becomes a source of potential zoonotic viruses [42,43].

Genetic analysis of the Uganda H9N2 viruses isolated during our studies showed that they belong to the G1 lineage, more related to Kenyan viruses than to viruses from Northern and Western Africa. The Ugandan viruses also clustered with viruses from the Middle East and Asia (Saudi Arabia, Pakistan, Iran and United Arab Emirates). Unfortunately, few countries in Eastern, Central and Southern Africa conduct regular surveillance for AIVs, making it difficult to accurately determine the spread of viruses across the continent and identify ancestral viruses. We can, however, infer from the available data that H9N2 viruses in Kenya and Uganda most likely originated from the Middle East via poultry movement or human activities. Viruses were introduced into Uganda around 2014 and became firmly established in commercial markets by 2017.

The overall genetic diversity in each segment of the Ugandan H9N2 viruses was less than 0.5%. This is likely influenced by the limited geographic distribution of the markets we sampled, where spread might be localized and not indicative of overall diversity across the country. Nevertheless, the nucleotide substitution rate within the H9N2 viruses has been steadily increasing, and pairwise distance analysis revealed clusters within the different gene segments that showed a random pattern of distribution by market. The topologies of the Ugandan H9N2 sequences within the phylogenetic trees differed by segment as did the rate of substitution. This demonstrates a differing pattern of evolution in the gene segments in addition to segment-specific differences in rates of nucleotide substitution.

A major driver for evolution and risk of influenza viruses lies in reassortment. Over the years, multiple novel influenza A reassortants have been discovered that are more virulent and can infect multiple species [44–46]. In order for reassortment to happen, there must be co-infection within the host cells [4]. This can be between viruses of different strains or subtypes, or viruses from different locations or host species [46]. In this case, we investigated reassortment within the isolates and reassortment by market, and found significant levels in both instances. Cocirculation of highly pathogenic AIV in the same population has been demonstrated to provide an opportunity for these viruses to reassert as has been observed with H5N1 in Cambodia [47] and in Egypt [44,45,48]. With H6N2 viruses having recently been discovered in Uganda [49], monitoring for reassortment must be intensified. Consistent with the limited genetic diversity in HA, selected H9N2 viruses were antigenically similar to each other. In addition, they were antigenically similar to H9N2 | A/Oman/2747/2019, a virus responsible for a recent human infection [19].

Molecular markers that have been reported to increase mammalian host susceptibility were also found in the Ugandan H9N2 viruses. These included the HA Q226L (H3 numbering) substitution which changes the flexibility of the 220-loop allowing it to interact more favorably with receptors on mammalian cells [50,51] and also promoting immune escape [52,53]. Also present in the HA were T197/G520 (H9 numbering), which have been associated with increased transmission of H9N2 viruses via respiratory droplets [54,55]. The canonical PB2 E627V/K mutations, responsible for increased virulence and adaptation in mice [56,57] were absent; however, this is compensated for by the presence of VDKGVV substitutions at positions 89, 309, 339, 477, 495 and 676 respectively [58,59].
Mutations that increase drug resistance against antivirals such as Amantadine and Baloxavir, which are the drugs of choice for treating severe Influenza infections, were also found in the MP and PA segments. This was in addition to mutations that increased pathogenicity in avian and mammalian cell lines, and in experimental animals such as ferrets and mice.

Our H9N2 viruses with PA-E199D demonstrated marginally lower susceptibility to the endonuclease inhibitor Baloxavir acid as has been reported for other viruses [37]. The virus with the combined F105L+E199D substitution had the highest EC50 which suggests a possible synergistic effect warranting further investigation. Many of the markers identified in the Ugandan viruses seem to be widespread in Africa, having been found in viruses from Kenya, Tunisia, and Algeria [60–62].

**5. Conclusions**

Our findings indicate that G1-lineage H9N2 viruses were introduced into Uganda from the Middle East around 2014. The viruses subsequently became established, quickly reaching a prevalence of 49%. Through comparative whole genome analysis, the viruses were found to have evolved extensively through reassortment events within and between the LBMs in segment-specific manners. There were also multiple molecular markers identified that have been associated with cross species infection; these include HA Q226L and those that enhance pathogenicity and reduce antiviral drug sensitivity. The documented zoonotic capacity of G1-lineage H9N2 viruses highlights the need for regular and exhaustive surveillance of all influenza viruses in circulation within the country and subsequent risk-assessment to provide a complete picture of prevalence, distribution, and risk.

**Author Contributions:** Conceptualization, G.A. and D.K.B.; Data curation, G.A.; Formal analysis, G.A.; Funding acquisition, D.K.B.; Investigation, G.A., T.T., Q.A.U., B.E., J.J. and R.J.W.; Methodology, G.A., T.T., Q.A.U., A.M., P.M., J.C.M.T., D.W., T.J. and R.G.W.; Project administration, B.E., H.K., E.M. and F.W.-M.; Supervision, M.F.D. and D.K.B.; Writing—original draft, G.A.; Writing—review and editing, R.J.W., M.F.D. and D.K.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** The study was funded by the United States Department of Defense Global Emerging Infections Surveillance Program: grant number P0136_19_KY_08 to DKB. Sequencing was supported by the National Institute of Allergy and Infectious Diseases, the National Institutes of Health, under contract number HHSN272201400006C.

**Institutional Review Board Statement:** All components of this study were approved by the Institution Review Board (IRB) of the College of Veterinary Medicine, Animal resources and Biosecurity, Makerere University under IRB reference; VET-0804.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Sequences analyzed in this work can be found in the GeneBank under accession numbers MW205009 to MW205198, MZ542846 to MZ543305, OP388450 to OP388679 and OP388789 to OP388827.

**Acknowledgments:** Thank you to Ashley Webb for support in whole genome sequencing and submission of sequences. Special thanks to Frederick Mulyazawo and Gift Ahimbisibwe.

**Conflicts of Interest:** The authors declare no conflict of interest.
Appendix A

<figure><img src="image_16.png" /><figcaption>**Figure A1.** Phylogenetic tree for PB2. Ugandan viruses are highlighted in blue, Kenyan in red, North and West African countries in pink, and the Middle East in the green clade. Posterior probabilities are shown on the branches. tMRCA for Uganda is November 2014. Scales represent substitutions and time.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>**Figure A2.** Phylogenetic tree for PB1. Ugandan viruses are highlighted in blue. Kenyan in red, North and West African countries in pink, and the Middle East, Asia and a few European viruses in the green clade. tMRCA for Uganda is August 2013. Posterior probabilities are shown on the branches. Scales represent substitutions and time.</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>**Figure A3.** Phylogenetic tree for PA. Ugandan viruses are highlighted in blue, Kenyan in red, North and West African countries in pink, and the Middle East in the green clade Including various non H9N2 viruses from Asia and Europe tMRCA for Uganda is September 2014. Posterior probabilities are shown on the branches. Scales represent substitutions and time.</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>**Figure A4.** Phylogenetic tree for NP. Ugandan viruses are highlighted in blue, Kenyan in red, North and West African countries in pink, and the Middle East and Asia in the green clade. tMRCA for Uganda is December 2013. Posterior probabilities are shown on the branches. Scales represent substitutions and time.</figcaption></figure>
<figure><img src="image_20.png" /><figcaption>**Figure A5.** Phylogenetic tree for NA. Ugandan viruses are highlighted in blue, Kenyan in red, North and West African countries in pink. tMRCA for Uganda is May 2014. Posterior probabilities are shown on the branches. Scales represent substitutions and time.</figcaption></figure>
<figure><img src="image_21.png" /><figcaption>**Figure A6.** Phylogenetic tree for MP. Ugandan viruses are highlighted in blue, Kenyan in red, North and West African countries in pink, and the Middle East in the green clade. tMRCA for Uganda is December 2014. Posterior probabilities are shown on the branches. Scales represent substitutions and time.</figcaption></figure>
<figure><img src="image_22.png" /></figure>

**Figure A7.** Phylogenetic tree for NS. Ugandan viruses are highlighted in blue, Kenyan in red, North and West African countries in pink, and the Middle East and Asia in the green clade. tMRCA for Uganda is May 2014. Posterior probabilities are shown on the branches. Scales represent substitutions and time.

## References

1. Peacock, T.P.; James, J.; Sealy, J.E.; Iqbal, M. A Global Perspective on H9N2 Avian Influenza Virus. *Viruses* **2019**, *11*, 620. [CrossRef] [PubMed]

2. Mertens, E.; Dugan, V.G.; Stockwell, T.B.; Lindsay, L.L.; Plancarte, M.; Boyce, W.M. Evaluation of Phenotypic Markers in Full Genome Sequences of Avian Influenza Isolates from California. *Comp. Immunol. Microbiol. Infect. Dis.* **2013**, *36*, 521–536. [CrossRef] [PubMed]

3. Shao, W.; Li, X.; Goraya, M.; Wang, S.; Chen, J.-L. Evolution of Influenza A Virus by Mutation and Re-Assortment. *Int. J. Mol. Sci.* **2017**, *18*, 1650. [CrossRef] [PubMed]

4. Lowen, A.C. It's in the Mix: Reassortment of Segmented Viral Genomes. *PLoS Pathog.* **2018**, *14*, e1007200. [CrossRef]

5. Rossman, J.S.; Lamb, R.A. Influenza Virus Assembly and Budding. *Virology* **2011**, *411*, 229–236. [CrossRef] [PubMed]
6. Arai, Y.; Ibrahim, M.S.; Elgendy, E.M.; Daidoji, T.; Ono, T.; Suzuki, Y.; Nakaya, T.; Matsumoto, K.; Watanabe, Y. Genetic Compatibility of Reassortants between Avian H5N1 and H9N2 Influenza Viruses with Higher Pathogenicity in Mammals. *J. Virol.* **2019**, *93*, e01969-18. [CrossRef] [PubMed]

7. Guo, Y.J.; Krauss, S.; Senne, D.A.; Mo, I.P.; Lo, K.S.; Xiong, X.P.; Norwood, M.; Shortridge, K.F.; Webster, R.G.; Guan, Y. Characterization of the Pathogenicity of Members of the Newly Established H9N2 Influenza Virus Lineages in Asia. *Virology* **2000**, *267*, 279-288. [CrossRef]

8. Naguib, M.M.; Ulrich, R.; Kasbohm, E.; Eng, C.L.P.; Hoffmann, D.; Grund, C.; Beer, M.; Harder, T.C. Natural Reassortants of Potentially Zoonotic Avian Influenza Viruses H5N1 and H9N2 from Egypt Display Distinct Pathogenic Phenotypes in Experimentally Infected Chickens and Ferrets. *J. Virol.* **2017**, *91*, e01300-17. [CrossRef]

9. Pu, J.; Wang, S.; Yin, Y.; Zhang, G.; Carter, R.A.; Wang, J.; Xu, G.; Sun, H.; Wang, M.; Wen, C.; et al. Evolution of the H9N2 Influenza Genotype That Facilitated the Genesis of the Novel H7N9 Virus. *Proc. Natl. Acad. Sci. USA* **2015**, *112*, 548-553. [CrossRef]

10. WHO. *Influenza at the Human-Animal Interface Summary and Assessment 27 June 2022*; World Health Organization: Geneva, Switzerland, 2022.

11. Zhang, M.; Zhao, C.; Chen, H.; Teng, Q.; Jiang, L.; Feng, D.; Li, X.; Yuan, S.; Xu, J.; Zhang, X.; et al. Internal Gene Cassette From a Human-Origin H7N9 Influenza Virus Promotes the Pathogenicity of H9N2 Avian Influenza Virus in Mice. *Front. Microbiol.* **2020**, *11*, 1441. [CrossRef]

12. Homme, P.J.; Easterday, B.C. Avian Influenza Virus Infection. I. Characteristics of Influenza A/Turkey/Wisconsin/1966 Virus. *Avian Dis.* **1970**, *14*, 66-74. [CrossRef] [PubMed]

13. Kaboudi, K. Low Pathogenic Avian Influenza Virus Subtype H9N2 in Poultry in North Africa: Current Status. *Vet. Sci. Res. Rev.* **2019**, *5*, 73-79. [CrossRef]

14. Nagy, A.; Mettenleiter, T.C.; Abdelwhab, E.M. A Brief Summary of the Epidemiology and Genetic Relatedness of Avian Influenza H9N2 Virus in Birds and Mammals in the Middle East and North Africa. *Epidemiol. Infect.* **2017**, *145*, 3320-3333. [CrossRef]

15. Kirunda, H.; Kibuuka, H.; Byaruhanga, A.; Mworozi, E.; Luswa, L.; Millard, M.; Wabwire-Mangen, F. Poor Biosecurity in Live Bird Markets in Uganda: A Potential Risk for Highly Pathogenic Avian Influenza Disease Outbreak in Poultry and Spread to Humans. *Int. J. Public Health Epidemiol.* **2014**, *3*, 67-74.

16. Fusade-Boyer, M.; Djegui, F.; Batawui, K.; Byuragaba, D.K.; Jones, J.C.; Wabwire-Mangeni, F.; Erima, B.; Atim, G.; Ukuli, Q.A.; Tugume, T.; et al. Antigenic and Molecular Characterization of Low Pathogenic Avian Influenza A(H9N2) Viruses in Sub-Saharan Africa from 2017 through 2019. *Emerg. Microbes Infect.* **2021**, *10*, 753-761. [CrossRef] [PubMed]

17. Public Health Agency of Canada. *Human Emerging Respiratory Pathogens Bulletin*; Public Health Agency of Canada: Ottawa, ON, Canada, 2022.

18. Qi, Y.; Guo, W.; Liu, C.; Li, W.; Gu, Y.; Li, S.; Chen, X. Seroprevalence of Influenza A (H9N2) Virus Infection among Humans in China: A Meta-Analysis. *Microb. Pathog.* **2021**, *155*, 104881. [CrossRef]

19. Almayahi, Z.K.; Al Kindi, H.; Davies, C.T.; Al-Rawahi, B.; Al-Jardani, A.; Al-Yaqoubi, F.; Jang, Y.; Jones, J.; Barnes, J.R.; Davis, W.; et al. First Report of Human Infection with Avian Influenza A(H9N2) Virus in Oman: The Need for a One Health Approach. *Int. J. Infect. Dis.* **2020**, *91*, 169-173. [CrossRef]

20. Potdar, V.; Hinge, D.; Satav, A.; Simões, E.A.F.; Yadav, P.D.; Chadha, M.S. Laboratory-Confirmed Avian Influenza A(H9N2) Virus Infection, India, 2019. *Emerg. Infect. Dis.* **2019**, *25*, 2328-2330. [CrossRef]

21. Jallow, M.M.; Fall, A.; Barry, M.A.; Diop, B.; Sy, S.; Goudiaby, D.; Fall, M.; Enouf, V.; Niang, M.N.; Dia, N. Genetic Characterization of the First Detected Human Case of Low Pathogenic Avian Influenza A/H9N2 in Sub-Saharan Africa, Senegal. *Emerg. Microbes Infect.* **2020**, *9*, 1092-1095. [CrossRef]

22. McCrone, J.T.; Woods, R.J.; Martin, E.T.; Malosh, R.E.; Monto, A.S.; Lauring, A.S. Stochastic Processes Constrain the within and between Host Evolution of Influenza Virus. *eLife* **2018**, *7*, e35962. [CrossRef]

23. Zhang, Y.; Aevermann, B.D.; Anderson, T.K.; Burke, D.F.; Dauphin, G.; Gu, Z.; He, S.; Kumar, S.; Larsen, C.N.; Lee, A.J.; et al. Influenza Research Database: An Integrated Bioinformatics Resource for Influenza Virus Research. *Nucleic Acids Res.* **2017**, *45*, D466-D474. [CrossRef] [PubMed]

24. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. *Mol. Biol. Evol.* **2021**, *38*, 3022-3027. [CrossRef] [PubMed]

25. Suchard, M.A.; Lemey, P.; Baele, G.; Ayres, D.L.; Drummond, A.J.; Rambaut, A. Bayesian Phylogenetic and Phylodynamic Data Integration Using BEAST 1.10. *Virus Evol.* **2018**, *4*, vey016. [CrossRef] [PubMed]

26. WOAH Chapter 3.3.4 Avian Influenza (Including Infection with High Pathogenicity Avian Influenza Viruses). In *Manual of Diagnostic Tests and Vaccines for Terrestrial Animals* 2022; World Organisation for Animal Health: Paris, France, 2022.

27. Wan, H.; Perez, D.R. Amino Acid 226 in the Hemagglutinin of H9N2 Influenza Viruses Determines Cell Tropism and Replication in Human Airway Epithelial Cells. *J. Virol.* **2007**, *81*, 5181-5191. [CrossRef] [PubMed]

28. Fan, S.; Deng, G.; Song, J.; Tian, G.; Suo, Y.; Jiang, Y.; Guan, Y.; Bu, Z.; Kawaoka, Y.; Chen, H. Two Amino Acid Residues in the Matrix Protein M1 Contribute to the Virulence Difference of H5N1 Avian Influenza Viruses in Mice. *Virology* **2009**, *384*, 28-32. [CrossRef]

29. Jiao, P.; Tian, G.; Li, Y.; Deng, G.; Jiang, Y.; Liu, C.; Liu, W.; Bu, Z.; Kawaoka, Y.; Chen, H. A Single-Amino-Acid Substitution in the NS1 Protein Changes the Pathogenicity of H5N1 Avian Influenza Viruses in Mice. *J. Virol.* **2008**, *82*, 1146-1154. [CrossRef]
30. Zhu, W.; Zou, X.; Zhou, J.; Tang, J.; Shu, Y. Residues 41V and/or 210D in the NP Protein Enhance Polymerase Activities and Potential Replication of Novel Influenza (H7N9) Viruses at Low Temperature. *Virol. J.* **2015**, *12*, 71. [CrossRef]

31. Yamayoshi, S.; Yamada, S.; Fukuyama, S.; Murakami, S.; Zhao, D.; Uraki, R.; Watanabe, T.; Tomita, Y.; Macken, C.; Neumann, G.; et al. Virulence-Affecting Amino Acid Changes in the PA Protein of H7N9 Influenza A Viruses. *J. Virol.* **2014**, *88*, 3127–3134. [CrossRef]

32. Elgendy, E.M.; Arai, Y.; Kawashita, N.; Daidoji, T.; Takagi, T.; Ibrahim, M.S.; Nakaya, T.; Watanabe, Y. Identification of Polymerase Gene Mutations That Affect Viral Replication in H5N1 Influenza Viruses Isolated from Pigeons. *J. Gen. Virol.* **2017**, *98*, 6–17. [CrossRef]

33. Feng, X.; Wang, Z.; Shi, J.; Deng, G.; Kong, H.; Tao, S.; Li, C.; Liu, L.; Guan, Y.; Chen, H. Glycine at Position 622 in PB1 Contributes to the Virulence of H5N1 Avian Influenza Virus in Mice. *J. Virol.* **2016**, *90*, 1872–1879. [CrossRef]

34. Graef, K.M.; Vreede, F.T.; Lau, Y.-F.; McCall, A.W.; Carr, S.M.; Subbarao, K.; Fodor, E. The PB2 Subunit of the Influenza Virus RNA Polymerase Affects Virulence by Interacting with the Mitochondrial Antiviral Signaling Protein and Inhibiting Expression of Beta Interferon. *J. Virol.* **2010**, *84*, 8433–8445. [CrossRef] [PubMed]

35. Kim, J.H.; Hatta, M.; Watanabe, S.; Neumann, G.; Watanabe, T.; Kawaoka, Y. Role of Host-Specific Amino Acids in the Pathogenicity of Avian H5N1 Influenza Viruses in Mice. *J. Gen. Virol.* **2010**, *91*, 1284–1289. [CrossRef] [PubMed]

36. Lan, Y.; Zhang, Y.; Dong, L.; Wang, D.; Huang, W.; Xin, L.; Yang, L.; Zhao, X.; Li, Z.; Wang, W.; et al. A Comprehensive Surveillance of Adamantane Resistance among Human Influenza A Virus Isolated from Mainland China between 1956 and 2009. *Antivir. Ther.* **2010**, *15*, 853–860. [CrossRef] [PubMed]

37. Hashimoto, T.; Baba, K.; Inoue, K.; Okane, M.; Hata, S.; Shishido, T.; Naito, A.; Wildum, S.; Omoto, S. Comprehensive Assessment of Amino Acid Substitutions in the Trimeric RNA Polymerase Complex of Influenza A Virus Detected in Clinical Trials of Baloxavir Marboxil. *Influenza Other Respir. Viruses* **2021**, *15*, 389–395. [CrossRef]

38. Tada, T.; Suzuki, K.; Sakurai, Y.; Kubo, M.; Okada, H.; Itoh, T.; Tsukamoto, K. NP Body Domain and PB2 Contribute to Increased Virulence of H5N1 Highly Pathogenic Avian Influenza Viruses in Chickens. *J. Virol.* **2011**, *85*, 1834–1846. [CrossRef] [PubMed]

39. Wasilenko, J.L.; Sarmento, L.; Pantin-Jackwood, M.J. A Single Substitution in Amino Acid 184 of the NP Protein Alters the Replication and Pathogenicity of H5N1 Avian Influenza Viruses in Chickens. *Arch. Virol.* **2009**, *154*, 969–979. [CrossRef]

40. Ayllon, J.; Domingues, P.; Rajsbaum, R.; Miorin, L.; Schmolke, M.; Hale, B.G.; Garcia-Sastre, A. A Single Amino Acid Substitution in the Novel H7N9 Influenza A Virus NS1 Protein Increases CPSF30 Binding and Virulence. *J. Virol.* **2014**, *88*, 12146–12151. [CrossRef] [PubMed]

41. Kirunda, H.; Erima, B.; Tumushabe, A.; Kiconco, J.; Tugume, T.; Mulei, S.; Mimbe, D.; Mworozi, E.; Bwogi, J.; Luswa, L.; et al. Prevalence of Influenza A Viruses in Livestock and Free-Living Waterfowl in Uganda. *BMC Vet. Res.* **2014**, *10*, 50. [CrossRef] [PubMed]

42. Sulaiman, L.; Shittu, I.; Fusaro, A.; Inuwa, B.; Zecchin, B.; Gado, D.; Schivo, A.; Bianco, A.; Laleye, A.; Gobbo, F.; et al. Live Bird Markets in Nigeria: A Potential Reservoir for H9N2 Avian Influenza Viruses. *Viruses* **2021**, *13*, 1445. [CrossRef]

43. Youk, S.; Lee, D.-H.; Jeong, J.-H.; Pantin-Jackwood, M.J.; Song, C.; Swayne, D.E. Live Bird Markets as Evolutionary Epicentres of H9N2 Low Pathogenicity Avian Influenza Viruses in Korea. *Emerg. Microbes Infect.* **2020**, *9*, 616–627. [CrossRef]

44. Elsayed, M.; Arafa, A.; Abdelwahab, S.; Hashish, A.; Youssef, A. Novel Reassortant of H9N2 Avian Influenza Viruses Isolated from Chickens and Quails in Egypt. *Vet. World* **2021**, *14*, 2142–2149. [CrossRef] [PubMed]

45. Kandeil, A.; El-Shesheny, R.; Maatouq, A.; Moatasim, Y.; Cai, Z.; McKenzie, P.; Webby, R.; Kayali, G.; Ali, M.A. Novel Reassortant H9N2 Viruses in Pigeons and Evidence for Antigenic Diversity of H9N2 Viruses Isolated from Quails in Egypt. *J. Gen. Virol.* **2017**, *98*, 548–562. [CrossRef]

46. Li, X.; Sun, J.; Lv, X.; Wang, Y.; Li, Y.; Li, M.; Liu, W.; Zhi, M.; Yang, X.; Fu, T.; et al. Novel Reassortant Avian Influenza A(H9N2) Virus Isolate in Migratory Waterfowl in Hubei Province, China. *Front. Microbiol.* **2020**, *11*, 220. [CrossRef] [PubMed]

47. Suttie, A.; Karlsson, E.A.; Deng, Y.-M.; Horm, S.V.; Yann, S.; Tok, S.; Sorn, S.; Holl, D.; Tum, S.; Hurt, A.C.; et al. Influenza A(H5N1) Viruses with A(H9N2) Single Gene (Matrix or PB1) Reassortment Isolated from Cambodian Live Bird Markets. *Virology* **2018**, *523*, 22–26. [CrossRef] [PubMed]

48. Mostafa, A.; Mahmoud, S.H.; Shehata, M.; Müller, C.; Kandeil, A.; El-Shesheny, R.; Nooh, H.Z.; Kayali, G.; Ali, M.A.; Pleschka, S. PA from a Recent H9N2 (G1-Like) Avian Influenza A Virus (AIV) Strain Carrying Lysine 367 Confers Altered Replication Efficiency and Pathogenicity to Contemporaneous H5N1 in Mammalian Systems. *Viruses* **2020**, *12*, 1046. [CrossRef]

49. Mercan, Y.; Atim, G.; Kayed, A.E.; Azbazdar, M.E.; Kandeil, A.; Ali, M.A.; Rubrum, A.; McKenzie, P.; Webby, R.J.; Erima, B.; et al. Molecular Characterization of Closely Related H6N2 Avian Influenza Viruses Isolated from Turkey, Egypt, and Uganda. *Viruses* **2021**, *13*, 607. [CrossRef]

50. Sun, Y.; Cong, Y.; Yu, H.; Ding, Z.; Cong, Y. Assessing the Effects of a Two-Amino Acid Flexibility in the Hemagglutinin 220-Loop Receptor-Binding Domain on the Fitness of Influenza A(H9N2) Viruses. *Emerg. Microbes Infect.* **2021**, *10*, 822–832. [CrossRef]

51. Zhang, J.; Ma, K.; Li, B.; Chen, Y.; Qiu, Z.; Xing, J.; Huang, J.; Hu, C.; Huang, Y.; Li, H.; et al. A Risk Marker of Tribasic Hemagglutinin Cleavage Site in Influenza A (H9N2) Virus. *Commun. Biol.* **2021**, *4*, 71. [CrossRef]

52. Li, Y.; Liu, M.; Sun, Q.; Zhang, H.; Zhang, H.; Jiang, S.; Liu, S.; Huang, Y. Genotypic Evolution and Epidemiological Characteristics of H9N2 Influenza Virus in Shandong Province, China. *Poult. Sci.* **2019**, *98*, 3488–3495. [CrossRef]
53. Sealy, J.E.; Yaqub, T.; Peacock, T.P.; Chang, P.; Ermetal, B.; Clements, A.; Sadeyen, J.-R.; Mehboob, A.; Shelton, H.; Bryant, J.E.; et al. Association of Increased Receptor-Binding Avidity of Influenza A(H9N2) Viruses with Escape from Antibody-Based Immunity and Enhanced Zoonotic Potential. *Emerg. Infect. Dis.* **2018**, *25*, 63–72. [CrossRef]

54. Sorrell, E.M.; Wan, H.; Araya, Y.; Song, H.; Perez, D.R. Minimal Molecular Constraints for Respiratory Droplet Transmission of an Avian-Human H9N2 Influenza A Virus. *Proc. Natl. Acad. Sci. USA* **2009**, *106*, 7565–7570. [CrossRef] [PubMed]

55. Xu, H.; Qian, J.; Song, Y.; Ming, D. The Adaptability of H9N2 Avian Influenza A Virus to Humans: A Comparative Docking Simulation Study. *Biochem. Biophys. Res. Commun.* **2020**, *529*, 963–969. [CrossRef] [PubMed]

56. Arai, Y.; Kawashita, N.; Ibrahim, M.S.; Elgendy, E.M.; Daidoji, T.; Ono, T.; Takagi, T.; Nakaya, T.; Matsumoto, K.; Watanabe, Y. PB2 Mutations Arising during H9N2 Influenza Evolution in the Middle East Confer Enhanced Replication and Growth in Mammals. *PLoS Pathog.* **2019**, *15*, e1007919. [CrossRef]

57. Sun, X.; Belser, J.A.; Maines, T.R. Adaptation of H9N2 Influenza Viruses to Mammalian Hosts: A Review of Molecular Markers. *Viruses* **2020**, *12*, 541. [CrossRef]

58. Chen, Z.; Huang, Q.; Yang, S.; Su, S.; Li, B.; Cui, N.; Xu, C. A Well-Defined H9N2 Avian Influenza Virus Genotype with High Adaption in Mammals Was Prevalent in Chinese Poultry Between 2016 to 2019. *Viruses* **2020**, *12*, 432. [CrossRef]

59. Li, J.; Ishaq, M.; Prudence, M.; Xi, X.; Hu, T.; Liu, Q.; Guo, D. Single Mutation at the Amino Acid Position 627 of PB2 That Leads to Increased Virulence of an H5N1 Avian Influenza Virus during Adaptation in Mice Can Be Compensated by Multiple Mutations at Other Sites of PB2. *Virus Res.* **2009**, *144*, 123–129. [CrossRef] [PubMed]

60. Kariithi, H.M.; Welch, C.N.; Ferreira, H.L.; Pusch, E.A.; Ateya, L.O.; Binepal, Y.S.; Apopo, A.A.; Dulu, T.D.; Afonso, C.L.; Suarez, D.L. Genetic Characterization and Pathogenesis of the First H9N2 Low Pathogenic Avian Influenza Viruses Isolated from Chickens in Kenyan Live Bird Markets. *Infect. Genet. Evol.* **2020**, *78*, 104074. [CrossRef]

61. Arbi, M.; Souiai, O.; Rego, N.; Larbi, I.; Naya, H.; Ghram, A.; Houimel, M. Historical Origins and Zoonotic Potential of Avian Influenza Virus H9N2 in Tunisia Revealed by Bayesian Analysis and Molecular Characterization. *Arch. Virol.* **2020**, *165*, 1527–1540. [CrossRef]

62. Barberis, A.; Boudaoud, A.; Gorrill, A.; Loupias, J.; Ghram, A.; Lachheb, J.; Alloui, N.; Ducatez, M.F. Full-Length Genome Sequences of the First H9N2 Avian Influenza Viruses Isolated in the Northeast of Algeria. *Virol. J.* **2020**, *17*, 108. [CrossRef]