# Avian Influenza a H9N2 Viruses in Morocco, 2018–2019

Fatima-Zohra Sikht, Mariette Ducatez, Charifa Drissi Touzani, Adam Rubrum, Richard Webby, Mohammed El Houadfi, Nour-Said Tligui, Christelle Camus, Siham Fellahi

► To cite this version:

Fatima-Zohra Sikht, Mariette Ducatez, Charifa Drissi Touzani, Adam Rubrum, Richard Webby, et al.. Avian Influenza a H9N2 Viruses in Morocco, 2018–2019. Viruses, 2022, 14 (3), pp.529. 10.3390/v14030529 . hal-03601568

HAL Id: hal-03601568

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

Submitted on 8 Mar 2022

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

# Avian Influenza a H9N2 Viruses in Morocco, 2018–2019

Fatima-Zohra Sikht <fnref n="1" />,<fnref n="2" />, Mariette Ducatez <fnref n="2" />, Charifa Drissi Touzani <fnref n="1" />, Adam Rubrum <fnref n="3" />, Richard Webby <fnref n="3" /> ID,
Mohammed El Houadfi <fnref n="1" />, Nour-Said Tligui <fnref n="4" />, Christelle Camus <fnref n="2" />,*,† ID and Siham Fellahi <fnref n="1" />,†

**Abstract:** Low pathogenic H9N2 avian influenza (LPAI H9N2) is considered one of the most important diseases found in poultry (broiler, laying hens, breeding chickens, and turkeys). This infection causes considerable economic losses. The objective of this work was to monitor and assess the presence of avian influenza virus (AIV) H9N2 in eight different regions of Morocco using real-time RT-PCR, and to assess the phylogenetic and molecular evolution of the H9N2 viruses between 2016 and 2019. Field samples were collected from 108 farms suspected of being infected with LPAI H9N2 virus. Samples were analyzed using H9N2-specific real-time RT-PCR. Highly positive samples were subjected to virus isolation and seven isolates were fully sequenced. Low pathogenic H9N2 avian influenza virus was introduced in Morocco in 2016. We show that in 2018–2019, the virus was still present irrespective of vaccination status. Phylogenetic and molecular analyses showed mutations related to virulence, although our viruses were related to 2016 Moroccan viruses and grouped in the G1 lineage. Specific amino acid substitutions were identified in Moroccan H9N2 viruses that are believed to lead to increased resistance to antiviral drugs.

**Keywords:** low pathogenic avian influenza virus; H9N2; Morocco; sequencing; full genome

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

check for
updates

**Citation:** Sikht, F.-Z.; Ducatez, M.; Touzani, C.D.; Rubrum, A.; Webby, R.; El Houadfi, M.; Tligui, N.-S.; Camus, C.; Fellahi, S. Avian Influenza a H9N2 Viruses in Morocco, 2018–2019. *Viruses* **2022**, *14*, 529. https://doi.org/10.3390/v14030529

Academic Editor: Feng Li

Received: 26 January 2022

Accepted: 2 March 2022

Published: 4 March 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/).

## 1. Introduction

Low pathogenic H9N2 avian influenza is an emerging disease that causes important economic losses in the poultry sector and is considered a threat to both poultry farms and public health.

Being a member of the genus *Alphainfluenzavirus*, and of the family *Orthomyxoviridae*, avian influenza viruses are enveloped RNA virus, with a genome composed of eight negative-sense RNA segments [1]. They are classified as low or highly pathogenic, on the basis of their virulence and hemagglutinin (HA) and neuraminidase (NA) sequences [2]. There are currently 18 HA and 11 NA identified, including the bat-specific H17–H18 and N10–N11 [3]. Avian influenza virus subtype H9N2 is pathotyped as a low pathogenic virus (LPAI). However, co-infections with other pathogens can lead to severe outbreaks with high mortality rates and severe economic losses [4,5].

H9N2 LPAIV had first been described in 1966 in a turkey farm in the USA [6]. Since then, it has been reported in numerous countries around the world. Between 1992 and 1994, an H9N2 outbreak occurred in Guangdong Province, China, and affected chicken farms with a mortality rate of 10% to 40%, with a reduction in egg-laying rate of 14% to 75% [7].

In 1996, H9N2 LPAIV was reported in South Korea [8]. In 1998, it was isolated from most provinces in China and, as a result, it was considered to be one of the most widespread

<footnote n="1">¹ Avian Pathology Unit, Department of Veterinary Pathology and Public Health, Agronomy and Veterinary Institute Hassan II, Rabat B.P. 6202, Morocco; fz.sikht@iav.ac.ma (F.-Z.S.); charifadrissi@gmail.com (C.D.T.); elhouadfimohammed@yahoo.fr (M.E.H.); fellahisiham2015@gmail.com (S.F.)</footnote>

<footnote n="2">² IHAP, Toulouse University, INRAE, ENVT, 31300 Toulouse, France; mariette.ducatez@envt.fr</footnote>

<footnote n="3">³ Department of Infectious Diseases, St. Jude Children's Research Hospital, Memphis, TN 38105, USA; adam.rubrum@stjude.org (A.R.); richard.webby@stjude.org (R.W.)</footnote>

<footnote n="4">⁴ Anatomo-Pathology Unit, Department of Veterinary Pathology and Public Health, Agronomy and Veterinary Institute Hassan II, Rabat B.P. 6202, Morocco; n.tligui@iav.ac.ma</footnote>

<footnote>* Correspondence: christelle.camus@envt.fr; Tel.: +33-5-61-19-38-80</footnote>

<footnote>† These authors contributed equally to this article.</footnote>
avian influenza virus in China [9]. Other countries in the Middle East and North Africa have been affected by this virus including Iran [10], Saudi Arabia [11], Jordan [12], the United Arab Emirates [13], Tunisia [14], Egypt [15], Sultanate of Oman [16], and Libya [17].

Phylogenetic analysis of the genome of LPAI H9N2 strains allowed to classify this virus into two distinct lineages: Eurasian and American. Though the Eurasian lineage contains several clades, most of the strains detected were classified in two clades (G1 and Y280) [18]. The G1 clade is represented by the A/Quail/Hong Kong/G1/1997 prototype virus, which mainly circulates in South China, Central Asia and the Middle East, while the Y280 clade viruses circulate throughout China and are represented by the A/Duck/Hong Kong/Y280/1997 prototype.

The main sources of LPAI H9N2 infections are the domestic and wild avian species. Wild birds are considered one of the reservoirs of the virus, and can transmit it over long distances. Transmission of the H9N2 virus can occur through direct contact with infected animals, and the infection can spread between farms through the movement of infected birds, vehicles, contaminated equipment or people with contaminated shoes or clothing [19]. In Pakistan, sparrows were shown to play a very important role in the transmission of the virus between farms [20]. In general, the sensitivity and receptivity of H9N2 is strongly dependent on the avian species (chicken and turkey). However, other species raised for consumption and/or hunting, such as guinea fowl, quail, pheasant, partridge, duck, goose, and ostrich are also considered sensitive. LPAI H9N2 virus has also been reported to be transmissible to mammals including dogs and cats [21] as well as humans [22–25].

The LPAIV H9N2 strain (SF1, GenBank accession number SCA48100) introduced in Morocco in January 2016, belongs to the G1 lineage, and is closely related to viruses circulating in the Middle East. As a response, the competent authorities authorized vaccination of any type of chicken as the best choice to limit the rapid spread of this disease [26]. However, in order to guide veterinarians to a rational choice of vaccines, it is important to determine and to phylogenetically analyze the circulating strains.

The aim of this study was to monitor the presence of LPAI H9N2 viruses in farms where animals with respiratory signs are reported using real-time reverse transcription PCR. Sequencing of isolates was performed in order to detect potential mutations that might affect the efficacy of commercial vaccines.

## 2. Materials and Methods

### 2.1. *Specimen Collection*

In collaboration with private veterinarians, a total of 151 samples, which included organs (trachea, lungs) and tracheal swabs, were collected from 108 commercial broilers farms (vaccinated and non-vaccinated) in eight regions of Morocco. Our sampling was based on chickens suspected of being infected with LPAI H9N2 virus, and presenting respiratory signs (rales, sneezing), associated with a decrease in food consumption and a drop in production. The specimens were collected in a period of 11 months, from 28/06/2018 to 31/05/2019.

### 2.2. *Samples Processing*

#### 2.2.1. RNA Extraction and Real Time RT-PCR

RNA extraction was performed using the NucleoSpin® RNA Virus Kit (Macherey-Nagel, Düren, Germany, No. 740956.250), following the manufacturer's instructions. In order to detect the H9N2 virus, the extracted RNA was amplified on the 7500 Fast Real-Time PCR System thermal cycler (Applied Biosystems, Foster City, CA, USA), using the primers and probe for generic detection of H9 subtypes described by [27], which target a conserved region in the HA2 subunit of the HA gene sequence.

#### 2.2.2. Virus Isolation

In order to obtain a maximum viral load detectable by conventional RT-PCR and for full genome sequencing purposes, 17 samples among those with the highest Ct in
RT-qPCR, from 17 different farms, were grown on 10-day-old, specific pathogen free (SPF) embryonated eggs. Briefly, the eggs were mirrored and the air chamber was delimited. The viral inoculums were prepared by mixing 0.2 mL of the viral suspension, 0.6 mL of sterile PBS and 0.2 mL of antibiotic OXY-Kel 20 L.A (oxytetracycline) and injected via allantoic cavity route using a sterile needle into the air chamber of the embryonated eggs. After viral inoculation, the eggs were incubated at 37 °C and examined daily for five days to assess the viability of the embryos. After the death of the embryo, eggs were refrigerated at 4 °C for 4 h. Then, the lesions on the embryos were observed and the allantoic fluids were collected, clarified, and stored at -80 °C until use.

### 2.2.3. Full Genome Amplification and Sequencing of H9N2 Moroccan Isolates

Viral RNA was extracted from allantoic fluids harvested from the 7 SPF embryonated eggs with the highest viral load, using the Macherey Nagel kit (Duren, Germany, No. 740956.250). Whole genome sequencing of 2018–2019 Moroccan isolates was performed with an Illumina MiSeq System (Illumina, San Diego, CA, USA) as previously described [28]. The preparation of libraries was performed using an Illumina Nextera XT library prep kit (FC-131-1096) (Illumina, San Diego, CA, USA) following the manufacturer's instructions. A tape station was used to verify the library quantity and quality. CLC genomic workbench was used for genomes assembly. The nucleotides sequences of all characterized strains in this study are submitted in the GenBank database under accession numbers summarized in Table 1.

<table><caption>**Table 1.** Accession numbers of segments sequences of studied Moroccan H9N2 viruses.</caption><thead><tr><th rowspan="2">Segment</th><th colspan="7">Strain</th></tr><tr><th>AS13</th><th>AS14</th><th>AS29</th><th>AS32</th><th>AS71</th><th>AS76</th><th>AS77</th></tr></thead><tbody><tr><td>PB2</td><td>MW165151</td><td>MW165079</td><td>MW165121</td><td>MW165089</td><td>MW165136</td><td>MW165110</td><td>MW165106</td></tr><tr><td>PB1</td><td>MW165154</td><td>MW165125</td><td>MW165122</td><td>MW165088</td><td>MW165142</td><td>MW165113</td><td>MW165101</td></tr><tr><td>PA</td><td>MW165158</td><td>MW165082</td><td>MW165117</td><td>MW165085</td><td>MW165139</td><td>MW165116</td><td>MW165103</td></tr><tr><td>HA</td><td>MW165152</td><td>MW165084</td><td>MW165120</td><td>MW165090</td><td>MW165137</td><td>MW165111</td><td>MW165105</td></tr><tr><td>NP</td><td>MW165157</td><td>MW165083</td><td>MW165124</td><td>MW165086</td><td>MW165140</td><td>MW165109</td><td>MW165108</td></tr><tr><td>NA</td><td>MW165155</td><td>MW165078</td><td>MW165119</td><td>MW165092</td><td>MW165135</td><td>MW165115</td><td>MW165104</td></tr><tr><td>NS</td><td>MW165156</td><td>MW165080</td><td>MW165123</td><td>MW165091</td><td>MW165141</td><td>MW165114</td><td>MW165102</td></tr><tr><td>M</td><td>MW165153</td><td>MW165081</td><td>MW165118</td><td>MW165087</td><td>MW165138</td><td>MW165112</td><td>MW165107</td></tr></tbody></table>

AS13, A/chicken/Morocco/AS13/2018; AS14, A/chicken/Morocco/AS14/2018; AS29, A/chicken/Morocco/AS29/2018; AS32, A/chicken/Morocco/AS32/2019; AS71, A/chicken/Morocco/AS71/2019; AS76, A/chicken/Morocco/AS76/2019; AS77, A/chicken/Morocco/AS77/2019.

### 2.3. Sequences and Phylogenetic Analyses

Bioedit 7.2.5 software [29] and ClustalW (version 1.83) [30] were used to compare and align nucleotide sequences of the complete genomes of seven Moroccan H9N2 isolates.

The phylogenetic tree was constructed by the maximum likelihood method, using the Mega 6.06 software [31]. The Blast [32] and Bioedit programs [30] were used to determine the sequence identity and compare the Moroccan strains with those selected from Genbank.

### 2.4. Statistical Analysis

Statistics describing the correlation between H9N2 positivity and the different factors: regions and vaccination status were calculated for each variable, including the mean and percentage distribution of frequencies. A non-parametric test (chi-squared test) was used to calculate the correlation between the H9N2 frequency in farms and their vaccination status.
# 3. Results

## 3.1. Case History and H9N2 Detection

One hundred and fifty-one samples from respiratory tissues and tracheal swabs were collected between June 2018 and May 2019 from different areas of Morocco: Fes-Meknes, Rabat-Sale-Kenitra, Casablanca-Settat, Draa-Tafilalet, Benimellal-Khenifra, Souss-Massa, Marrakech-Safi, and the eastern region. The samples were tested by real time RT-PCR to detect the presence of influenza virus. A total of 83 were positive for AIV with cycle threshold (Ct) values varying from 12 to 39 (Table A1), of which 40%, 56%, and 4% of the samples had a Ct below or equal to 25, between 25 and 35, and above 35, respectively. The epidemiological survey resulted in a positivity rate of the disease of 58% (63 positive farms out of 108 sampled farms) (Table 2).

<table><caption>**Table 2.** Farms positivity rate.</caption><thead><tr><th></th><th>Number of Farms</th><th>Positive Farms</th><th>Positivity Rate</th></tr></thead><tbody><tr><td>Fes-Meknes</td><td>34</td><td>20</td><td>59%</td></tr><tr><td>Rabat-Sale-Kenitra</td><td>18</td><td>8</td><td>44%</td></tr><tr><td>Casablanca-Settat</td><td>16</td><td>13</td><td>81%</td></tr><tr><td>Draa-Tafilalet</td><td>1</td><td>1</td><td>100%</td></tr><tr><td>BeniMellal-Khenifra</td><td>4</td><td>3</td><td>75%</td></tr><tr><td>Souss-Massa</td><td>26</td><td>17</td><td>65%</td></tr><tr><td>Oriental</td><td>4</td><td>1</td><td>25%</td></tr><tr><td>Marrakech-Safi</td><td>5</td><td>0</td><td>0%</td></tr><tr><td>Morocco (Total)</td><td>108</td><td>63</td><td>58%</td></tr><tr><td>Vaccinated</td><td>44 (41%)</td><td>31</td><td>70%</td></tr><tr><td>Unvaccinated</td><td>64 (59%)</td><td>32</td><td>50%</td></tr></tbody></table>

## 3.2. Vaccination Status

The positivity rate (relative prevalence) of LPAI H9N2 positive farms was estimated to be 50% in unvaccinated farms (32 positive farms out of 64 chicken unvaccinated farms tested), while it was 70% in vaccinated farms (31 positive farms out of 44 chicken vaccinated farms tested). The overall vaccination rate reached 41% (44 out of 108 farms tested against H9N2) (Table 2).

The presence of LPAIV H9N2 was detected differently between the groups of vaccinated and unvaccinated farms, but the difference was not statistically significant (95% CI, *p* value: 0.9).

## 3.3. Molecular Characterization and Phylogenetic Analysis of the Eight Viral Segments

Viruses from highly positive samples were isolated. The genome of 7 of them was fully sequenced with an IlluminaMiSeq System [28]. Phylogenetic analysis showed that our Moroccan H9N2 viruses isolated from chickens were in the same cluster as the other Moroccan viruses detected in 2016, and grouped into G1 lineage. They were compared with relevant virus sequences available in GenBank.

Based on HA and NA phylogenetic trees, our Moroccan viruses were closely related to viruses previously isolated in the Emirates (2015), Morocco (2016), Burkina Faso (2017), and Algeria (2017), with bootstrap values of 100 and 60 for HA and NA, respectively (Figure 1). Regarding the internal genes, they grouped with the Moroccan viruses of 2016–2017, Algerian viruses of 2017, and Ghana viruses of 2017–2018 (Figure A1).
(a)

(b)

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Phylogenetic trees of Moroccan HA (a) and NA (b) gene segments. The nucleotide sequences of Moroccan H9N2 viruses (black squares) characterized in this study were compared with relevant virus sequences available in GenBank and GISAID databases, reference viruses, and relevant sequences from neighboring areas.</figcaption></figure>

The sequence analyses of the seven Moroccan isolates showed several substitutions in both HA and NA sequences when compared to 2016 strain SF1 (Tables A2 and A3).
All seven Moroccan isolates had the RSSR*GLF motif at the HA cleavage site, which is a characteristic and signature of the low pathogenic H9N2 viruses.

Potential HA glycosylation sites were identical to 2016 Moroccan viruses (29 NSTE, 82 NPSC, 105 NGTC, 141 NVTY, 298 NSTM, 492 NGTY, H3 numbering throughout), except for position site (297 NISK→NVSK) for four out of seven samples.

Our viruses did not present HA Receptor Binding Site (RBS) sequence associated with greater affinity for 6'-sialylacetyllactosamine (6SLN) [33] (Table 3), nor, when compared to Moroccan 2016 viruses, new critical amino acids defined as supporting mammalian replication [34–36].

<table><caption>**Table 3.** Residues associated to 6'-sialylacetyllactosamine-affinity and to drug resistance.</caption><thead><tr><th>HA * RBS</th><th>190</th><th>Q226L</th><th>227</th><th>NA</th><th>M2</th></tr></thead><tbody><tr><td>SF1 **</td><td>V</td><td>L</td><td>I</td><td>Q</td><td>N</td></tr><tr><td>AS13</td><td>T</td><td>L</td><td>I</td><td>Q</td><td>N</td></tr><tr><td>AS14</td><td>A</td><td>L</td><td>I</td><td>Q</td><td>N</td></tr><tr><td>AS29</td><td>V</td><td>L</td><td>I</td><td>Q</td><td>N</td></tr><tr><td>AS32</td><td>A</td><td>L</td><td>I</td><td>Q</td><td>N</td></tr><tr><td>AS71</td><td>A</td><td>L</td><td>I</td><td>Q</td><td>N</td></tr><tr><td>AS76</td><td>A</td><td>Q</td><td>I</td><td>Q</td><td>N</td></tr><tr><td>AS77</td><td>A</td><td>Q</td><td>I</td><td>Q</td><td>N</td></tr></tbody></table>

HA, hemagglutinin; RBS, receptor-binding site. * H3 numbering; ** GenBank accession number: SCA48100.

Among mutations associated to resistance to antiviral molecules, no changes from SF1 strain were identified for NA H274Y substitution [37] or M2 S31N mutation [38] (Table 3).

## 4. Discussion

Our analysis included 151 field samples from 108 poultry farms suspected of being infected by the LPAI H9N2 virus and presenting mainly respiratory signs as well as poor zootechnical performances (decrease in production, decrease in feed consumption and mortalities), as reported previously [39].

The results revealed that 58% of the samples were positive for LPAIV H9N2. However, we cannot extrapolate these results on the epidemiological profile of the LPAI H9N2 virus in Morocco since our sampling was not representative enough of the national territory and only 8 regions out of 12 (Fes-Meknes, Rabat-Sale-Kenitra, Casablanca-Settat, Draa-Tafilalet, Béni Mellal-Khenifra, Souss-Massa, Marrakech-Safi and the Oriental) were sampled.

The spread of LPAI H9N2 virus in Morocco can be explained mainly by the movement of farmers, workers, and feed suppliers without compliance with biosecurity rules, to which can be added the transport of live chickens [13,15,26]. It should be noted that the application of sanitary biosecurity measures in broiler farms has been shown to be insufficient to prevent the entry of the virus [26]. In addition, mutations associated with resistance to antiviral molecules are still present in our LPAI H9N2 strains. The M2 S31N mutation is known to increase resistance to antiviral molecules, especially amantadine and rimantadine [38]. Likewise, other studies have shown that the absence of the H274Y substitution in the NA protein can confer to the virus a sensitivity to neuraminidase inhibitors such as oseltamivir [37].

The positivity rate of AI H9N2 positive farms was estimated to be 50% in unvaccinated chicken farms (32 positive farms out of 64 chicken-unvaccinated farms tested), while it reached 70% in vaccinated chicken farms (31 positive farms out of 44 chicken-vaccinated farms tested). A recent study evaluating the efficacy of four different commercial vaccines on H9N2 LPAIV SF1 strain has shown that they conferred a very limited protection against
the infection [40]. Park and collaborators indicated that vaccination against H9N2 virus coupled with continuous infection of vaccinated flocks is an advantage for mutant viruses selection [41], whereas other studies report that vaccination decreases viral pressure in the field by reducing the level and duration of viral shedding [42]. Other explanations include the quality of the vaccine (either that it was not of the same strain as the virus currently circulating, or that it had a low antigen concentration [43,44]). We were not able to gather sufficient information relative to the vaccines used and their composition to be able to conclude on this point.

In addition, the vaccination rate against H9N2 was low (41%), which can be explained by the high cost of vaccination or by the fact that some farmers consider the vaccine ineffective.

For unvaccinated specimens that tested negative, despite respiratory signs, other respiratory diseases, including BI or NDV, might be the cause of the observed clinical signs.

In this study, we demonstrated a relationship between our viruses, isolated in 2018–2019, and those isolated in Morocco in 2016, which all belong to the G1 lineage [26]. In order to evaluate the evolution of the Moroccan H9N2 virus over time (i.e., after its first introduction into Morocco), phylogenetic and genetic analyses were carried out.

On the HA and NA phylogeny, the 2018–2019 Moroccan viruses were close to those from Emirates (2015), Morocco (2016), Burkina Faso (2017), and Algeria (2017). As for internal genes, they were grouped in the same cluster as the Moroccan viruses of 2016–2017, Algerian viruses of 2017, and Ghana viruses of 2017–2018. This similarity can be explained by the common border between Morocco and Algeria, and by the history of commercial exchanges within western Africa countries. The evolution of the influenza virus directly depends on its genomic properties, which leads us to follow and verify the presence of possible mutations over time (especially on the HA and NA genes, which are the main proteins targeted by antibodies). Moroccan viruses harbor several mutations in HA and NA. Some have already been characterized, such as HA Q226L, which is known to enhance binding to mammalian-like receptors [45]. For other mutations, further studies are necessary in order to determine whether they could affect the virulence of the virus in poultry, or increase transmissibility to human. Potential glycosylation sites were identified in our Moroccan strains. As compared to Moroccan strains isolated in 2016, there was one amino acid change within a glycosylation site (297 NISK→NVSK, H3 numbering throughout) in four of the seven sequenced isolates. Changes in glycosylation sites may affect the host range and virulence of influenza viruses [46], though we do not know if it is the case here.

**5. Conclusions**

The low pathogenic avian influenza virus H9N2 is endemic within the country despite vaccination. Biosecurity issues in farm management, combined with high mutation potential are likely to cause dynamic changes in LPAI H9N2 strains. This prompts us to propose appropriate surveillance and adaptation of vaccines to circulating strains in order to better understand and fight public health risks.

**Author Contributions:** F.-Z.S. participated in the design of the study, drafted the manuscript, isolated the viruses and realized the genetic and the sequence alignment phylogenetic studies. S.F. participated in the design of the study, helped in isolation of the viruses and corrected the manuscript. A.R. and R.W. carried out the sequencing. C.C. helped analyzing the data and corrected the manuscript. M.D. corrected the manuscript. C.D.T., N.-S.T., and M.E.H. helped in the draft of the manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by PHC TOUBKAL 19/76, grant number 41439WM and by the National Institute of Allergy and Infectious Diseases, National Institutes of Health (CEIRS HHSN272201400006C).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.
**Data Availability Statement:** Not applicable.

**Acknowledgments:** We thank partners and veterinarians for their help in this project. This project was financially supported by Ministry of Europe and Foreign Affairs, Ministry of Higher Education, Research and Innovation and the French Institute of Rabat (PHC TOUBKAL 2019 (French-Morocco bilateral program) Grant Number: 41439WM. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

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

## Appendix A

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

<figure><img src="image_6.png" /><figcaption>**Figure A1.** *Cont.*</figcaption></figure>
(C)

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

(D)

<figure><img src="image_8.png" /><figcaption>**Figure A1.** *Cont.*</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>(F)</figcaption></figure>

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

**Figure A1.** Phylogenetic trees of Moroccan PB2 (A), PB1 (B), PA (C), NP (D), M (E), and NS (F) segments. The nucleotide sequences of Moroccan H9N2 viruses characterized in this study (black squares) were compared with relevant virus sequences available in GenBank and GISAID databases, reference viruses, and relevant sequences from neighboring areas.
<table><caption>Table A1. List and status of samples used in this study.</caption><thead><tr><th>Sample</th><th>Farm</th><th>Sampling Date</th><th>Location of Farm</th><th>Age of Birds</th><th>H9N2 Vaccination Status</th><th>H9N2 RT-qPCR Ct Value</th><th>Farms Status</th></tr></thead><tbody><tr><td>AS1</td><td>F1</td><td>28/06/2018</td><td>Rabat</td><td>27D</td><td>V</td><td>25</td><td>P</td></tr><tr><td>AS2</td><td rowspan="3">F2</td><td rowspan="3">16/07/2018</td><td rowspan="3">Rabat</td><td>46D</td><td rowspan="3">V</td><td>-</td><td rowspan="3">N</td></tr><tr><td>AS3</td><td>46D</td><td>-</td></tr><tr><td>AS4</td><td>46D</td><td>-</td></tr><tr><td>AS5</td><td rowspan="2">F3</td><td rowspan="2">08/08/2018</td><td rowspan="2">Oriental</td><td>45D</td><td rowspan="2">NV</td><td>-</td><td rowspan="2">N</td></tr><tr><td>AS6</td><td>45D</td><td>-</td></tr><tr><td>AS7</td><td>F4</td><td>08/08/2018</td><td>Oriental</td><td>36D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS8</td><td>F5</td><td>17/07/2018</td><td>Casablanca</td><td>28D</td><td>V</td><td>-</td><td>N</td></tr><tr><td>AS9</td><td>F6</td><td>17/09/2018</td><td>Temara</td><td>43D</td><td>V</td><td>-</td><td>N</td></tr><tr><td>AS11</td><td>F7</td><td>19/10/2018</td><td>Meknes</td><td>32D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS12B1 *</td><td rowspan="2">F8</td><td>23/10/2018</td><td rowspan="2">Kenitra</td><td>23D</td><td rowspan="2">NV</td><td>20</td><td rowspan="2">P</td></tr><tr><td>AS12B2</td><td>23/10/2018</td><td>23D</td><td>33</td></tr><tr><td>AS13 *+</td><td>F9</td><td>18/10/2018</td><td>El hajeb</td><td>32D</td><td>NV</td><td>17</td><td>P</td></tr><tr><td>AS14 *+</td><td>F10</td><td>24/10/2018</td><td>Meknes</td><td>30D</td><td>NV</td><td>18</td><td>P</td></tr><tr><td>AS15</td><td>F11</td><td>21/10/2018</td><td>Meknes</td><td>37D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS16 *</td><td>F12</td><td>17/10/2018</td><td>Midelt</td><td>40D</td><td>NV</td><td>18</td><td>P</td></tr><tr><td>AS17</td><td>F13</td><td>11/12/2018</td><td>Fes</td><td>38D</td><td>NV</td><td>26</td><td>P</td></tr><tr><td>AS18</td><td>F14</td><td>11/12/2018</td><td>Fes</td><td>35D</td><td>NV</td><td>26</td><td>P</td></tr><tr><td>AS19</td><td>F15</td><td>11/12/2018</td><td>Fes</td><td>41D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS20</td><td>F16</td><td>23/11/2018</td><td>Fes</td><td>32D</td><td>V</td><td>25</td><td>P</td></tr><tr><td>AS21</td><td>F17</td><td>12/11/2018</td><td>Salé</td><td>36D</td><td>V</td><td>-</td><td>N</td></tr><tr><td>AS22</td><td rowspan="3">F18</td><td>07/12/2019</td><td rowspan="3">Meknes</td><td>34D</td><td rowspan="3">NV</td><td>22</td><td rowspan="3">P</td></tr><tr><td>AS23</td><td>07/12/2019</td><td>34D</td><td>26</td></tr><tr><td>AS24</td><td>07/12/2019</td><td>34D</td><td>22</td></tr><tr><td>AS26</td><td>F19</td><td>19/02/2019</td><td>Benslimane</td><td>30D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS27</td><td>F20</td><td>20/02/2019</td><td>Rabat</td><td>24D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AAS28</td><td>F21</td><td>04/11/2018</td><td>Ait brahim (Fes)</td><td>32D</td><td>V</td><td>24</td><td>P</td></tr><tr><td>AS29 *+</td><td>F22</td><td>16/11/2018</td><td>Fes</td><td>36D</td><td>V</td><td>16</td><td>P</td></tr><tr><td>AS30</td><td>F23</td><td>26/01/2019</td><td>Hajeb</td><td>20D</td><td>V</td><td>24</td><td>P</td></tr><tr><td>AS31</td><td>F24</td><td>11/02/2019</td><td>Sefrou</td><td>30-36D</td><td>V</td><td>17</td><td>P</td></tr><tr><td>AS32 *+</td><td>F25</td><td>12/02/2019</td><td>Ain chegag</td><td>37D</td><td>V</td><td>19</td><td>P</td></tr><tr><td>AS33 *</td><td>F26</td><td>13/02/2019</td><td>Zerarda tahla</td><td>40D</td><td>NV</td><td>16</td><td>P</td></tr><tr><td>AS34</td><td>F27</td><td>25/01/2019</td><td>Khemisset</td><td>33D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS35</td><td>F28</td><td>04/03/2019</td><td>Meknes</td><td>32D</td><td>-</td><td>-</td><td>N</td></tr><tr><td>AS36</td><td>F29</td><td>13/02/2019</td><td>Fes</td><td>34D</td><td>-</td><td>14</td><td>P</td></tr><tr><td>AS37</td><td>F30</td><td>13/02/2019</td><td>Fes</td><td>37D</td><td>V</td><td>-</td><td>N</td></tr><tr><td>AS38</td><td>F31</td><td>01/02/2019</td><td>Meknes</td><td>32D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS39</td><td>F32</td><td>12/02/2019</td><td>Hadeb</td><td>34D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS40</td><td>F33</td><td>12/02/2019</td><td>Salé</td><td>42D</td><td>NV</td><td>12</td><td>P</td></tr><tr><td>AS41</td><td>F34</td><td>25/02/2019</td><td>Khemisset</td><td>44D</td><td>V</td><td>-</td><td>N</td></tr><tr><td>AS42</td><td>F35</td><td>25/02/2019</td><td>Meknes</td><td>35D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS43 *</td><td>F36</td><td>07/02/2019</td><td>Casablanca</td><td>29D</td><td>V</td><td>16</td><td>P</td></tr><tr><td>AS44</td><td>F37</td><td>07/02/2019</td><td>Casablanca</td><td>32D</td><td>V</td><td>16</td><td>P</td></tr><tr><td>AS45 *</td><td rowspan="2">F38</td><td>10/02/2019</td><td rowspan="2">Rabat</td><td>45D</td><td rowspan="2">V</td><td>17</td><td rowspan="2">P</td></tr><tr><td>AS46</td><td>10/02/2019</td><td>45D</td><td>14</td></tr></tbody></table>
<table><caption>Table A1. *Cont.*</caption><thead><tr><th>Sample</th><th>Farm</th><th>Sampling Date</th><th>Location of Farm</th><th>Age of Birds</th><th>H9N2 Vaccination Status</th><th>H9N2 RT-qPCR Ct Value</th><th>Farms Status</th></tr></thead><tbody><tr><td>AS47</td><td>F39</td><td>11/02/2019</td><td>Tiflet</td><td>28D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS54</td><td>F40</td><td>20/02/2019</td><td>Rabat</td><td>32D</td><td>V</td><td>24</td><td>P</td></tr><tr><td>AS55</td><td>F41</td><td>20/02/2019</td><td>Oriental</td><td>37D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS56 *</td><td>F42</td><td>20/02/2019</td><td>Oriental</td><td>38D</td><td>NV</td><td>26</td><td>P</td></tr><tr><td>AS57</td><td>F43</td><td>21/02/2019</td><td>Taza</td><td>29D</td><td>NV</td><td>26</td><td>P</td></tr><tr><td>AS58</td><td>F44</td><td>21/02/2019</td><td>Tahla</td><td>36D</td><td>NV</td><td>-</td><td>P</td></tr><tr><td>AS59</td><td>F44</td><td>21/02/2019</td><td>Tahla</td><td>36D</td><td>NV</td><td>26</td><td>P</td></tr><tr><td>AS60</td><td>F45</td><td>21/02/2019</td><td>Tahla</td><td>39D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS61 *</td><td>F46</td><td>24/02/2019</td><td>Meknes</td><td>44D</td><td>NV</td><td>21</td><td>P</td></tr><tr><td>AS62</td><td>F47</td><td>24/02/2019</td><td>Elhajeb</td><td>37D</td><td>NV</td><td>23</td><td>P</td></tr><tr><td>AS63</td><td>F47</td><td>24/02/2019</td><td>Elhajeb</td><td>37D</td><td>NV</td><td>12</td><td>P</td></tr><tr><td>AS64</td><td>F48</td><td>24/02/2019</td><td>Elhajeb</td><td>35D</td><td>NV</td><td>26</td><td>P</td></tr><tr><td>AS65</td><td>F49</td><td>25/02/2019</td><td>Meknes</td><td>41D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS66</td><td>F50</td><td>25/02/2019</td><td>Meknes</td><td>35D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS67</td><td>F51</td><td>25/02/2019</td><td>Meknes</td><td>32D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS68 *</td><td>F52</td><td>26/02/2019</td><td>Khenifra</td><td>38D</td><td>NV</td><td>14</td><td>P</td></tr><tr><td>AS69</td><td>F53</td><td>25/02/2019</td><td>Khenifra</td><td>28D</td><td>NV</td><td>25</td><td>P</td></tr><tr><td>AS70</td><td>F54</td><td>28/02/2019</td><td>Salé</td><td>23D</td><td>NV</td><td>22</td><td>P</td></tr><tr><td>AS71 **</td><td>F55</td><td>01/04/2019</td><td>Sidi slimane</td><td>14D</td><td>NV</td><td>22</td><td>P</td></tr><tr><td>AS72</td><td>F56</td><td>11/12/2017</td><td>Khenifra</td><td>33D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS73</td><td>F57</td><td>20/02/2019</td><td>Khenifra</td><td>36D</td><td>NV</td><td>14</td><td>P</td></tr><tr><td>AS74</td><td>F58</td><td>11/04/2019</td><td>Sidi slimane</td><td>37D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS75</td><td>F59</td><td>06/03/2019</td><td>Salé</td><td>29D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS76 **</td><td>F60</td><td>23/02/2019</td><td>Meknes</td><td>38D</td><td>NV</td><td>20</td><td>P</td></tr><tr><td>AS77 **</td><td>F61</td><td>17/01/2019</td><td>Meknes</td><td>29D</td><td>NV</td><td>16</td><td>P</td></tr><tr><td>AS78</td><td>F62</td><td>01/04/2019</td><td>Tiznit</td><td>25D</td><td>NV</td><td>-</td><td>P</td></tr><tr><td>AS79</td><td>F62</td><td>01/04/2019</td><td>Tiznit</td><td>25D</td><td>NV</td><td>30</td><td>P</td></tr><tr><td>AS80</td><td>F63</td><td>01/04/2019</td><td>Tiznit</td><td>32D</td><td>NV</td><td>30</td><td>P</td></tr><tr><td>AS81</td><td>F63</td><td>01/04/2019</td><td>Tiznit</td><td>32D</td><td>NV</td><td>28</td><td>P</td></tr><tr><td>AS82</td><td>F64</td><td>01/04/2019</td><td>Tiznit</td><td>32D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS83</td><td>F64</td><td>01/04/2019</td><td>Tiznit</td><td>37D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS84</td><td>F65</td><td>01/04/2019</td><td>Tiznit</td><td>37D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS85 *</td><td>F66</td><td>01/04/2019</td><td>Tiznit</td><td>28D</td><td>NV</td><td>16</td><td>P</td></tr><tr><td>AS86</td><td>F67</td><td>01/04/2019</td><td>Tiznit</td><td>37D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS87</td><td>F67</td><td>01/04/2019</td><td>Tiznit</td><td>34D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS88</td><td>F68</td><td>01/04/2019</td><td>Tiznit</td><td>34D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS89</td><td>F69</td><td>01/04/2019</td><td>Tiznit</td><td>15D</td><td>NV</td><td>26</td><td>P</td></tr><tr><td>AS90</td><td></td><td>31/05/2019</td><td></td><td>32D</td><td></td><td>-</td><td></td></tr><tr><td>AS91</td><td></td><td>31/05/2019</td><td></td><td>32D</td><td></td><td>-</td><td></td></tr><tr><td>AS92</td><td></td><td>31/05/2019</td><td></td><td>32D</td><td></td><td>-</td><td></td></tr><tr><td>AS93</td><td></td><td>31/05/2019</td><td></td><td>32D</td><td></td><td>-</td><td></td></tr><tr><td>AS94</td><td>F70</td><td>31/05/2019</td><td>Rabat</td><td>32D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>AS95</td><td></td><td>31/05/2019</td><td></td><td>32D</td><td></td><td>-</td><td></td></tr><tr><td>AS96</td><td></td><td>31/05/2019</td><td></td><td>32D</td><td></td><td>-</td><td></td></tr><tr><td>AS97</td><td></td><td>31/05/2019</td><td></td><td>32D</td><td></td><td>-</td><td></td></tr><tr><td>AS98</td><td></td><td>31/05/2019</td><td></td><td>32D</td><td></td><td>-</td><td></td></tr><tr><td>BS1</td><td>F71</td><td>09/04/2019</td><td>Fes</td><td>28D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>BS2</td><td>F72</td><td>09/04/2019</td><td>Fes</td><td>36D</td><td>NV</td><td>-</td><td>N</td></tr></tbody></table>
<table><caption>Table A1. *Cont.*</caption><thead><tr><th>Sample</th><th>Farm</th><th>Sampling Date</th><th>Location of Farm</th><th>Age of Birds</th><th>H9N2 Vaccination Status</th><th>H9N2 RT-qPCR Ct Value</th><th>Farms Status</th></tr></thead><tbody><tr><td>BS3</td><td rowspan="2">F73</td><td>16/09/2019</td><td rowspan="2">Ait moussa</td><td>35D</td><td rowspan="2">NV</td><td>-</td><td rowspan="2">N</td></tr><tr><td>BS4</td><td>16/09/2019</td><td>35D</td><td>-</td></tr><tr><td>BS5</td><td rowspan="2">F74</td><td>08/10/2019</td><td rowspan="2">Ait moussa</td><td>36D</td><td rowspan="2">V</td><td>-</td><td rowspan="2">N</td></tr><tr><td>BS6</td><td>09/10/2019</td><td>36D</td><td>-</td></tr><tr><td>BS7</td><td>F75</td><td>23/10/2019</td><td>Marrakech</td><td>18D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>BS8</td><td>F76</td><td>24/10/2019</td><td>Haouz</td><td>29D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>BS9</td><td>F77</td><td>30/10/2019</td><td>Marrakech</td><td>34D</td><td>V</td><td>-</td><td>N</td></tr><tr><td>BS10</td><td>F78</td><td>31/10/2019</td><td>Rhamna</td><td>13D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>BS11</td><td>F79</td><td>05/11/2019</td><td>Marrakech</td><td>12D</td><td>NV</td><td>-</td><td>N</td></tr><tr><td>BS12</td><td>F80</td><td>23/11/2019</td><td>Casablanca</td><td>30D</td><td>NV</td><td>29</td><td>P</td></tr><tr><td>BS13</td><td>F81</td><td>23/11/2019</td><td>Casablanca</td><td>21D</td><td>NV</td><td>28</td><td>P</td></tr><tr><td>BS14</td><td>F82</td><td>23/11/2019</td><td>Casablanca</td><td>24D</td><td>NV</td><td>33</td><td>P</td></tr><tr><td>BS15</td><td>F83</td><td>23/11/2019</td><td>Casablanca</td><td>24D</td><td>NV</td><td>31</td><td>P</td></tr><tr><td>BS18</td><td rowspan="3">F84</td><td>25/11/2019</td><td rowspan="3">Tiznit</td><td>34D</td><td rowspan="3">V</td><td>-</td><td rowspan="3">N</td></tr><tr><td>BS19</td><td>25/11/2019</td><td>34D</td><td>-</td></tr><tr><td>BS20</td><td>25/11/2019</td><td>34D</td><td>-</td></tr><tr><td>BS21</td><td rowspan="3">F85</td><td>05/12/2019</td><td rowspan="3">Tiznit</td><td>29D</td><td rowspan="3">V</td><td>-</td><td rowspan="3">N</td></tr><tr><td>BS22</td><td>05/12/2019</td><td>29D</td><td>-</td></tr><tr><td>BS23</td><td>05/12/2019</td><td>29D</td><td>-</td></tr><tr><td>BS47</td><td>F86</td><td>14/11/2019</td><td>Casablanca</td><td>28D</td><td>V</td><td>26</td><td>P</td></tr><tr><td>BS48</td><td>F87</td><td>14/11/2019</td><td>Casablanca</td><td>28D</td><td>V</td><td>26</td><td>P</td></tr><tr><td>BS49</td><td>F88</td><td>14/11/2019</td><td>Casablanca</td><td>28D</td><td>V</td><td>27</td><td>P</td></tr><tr><td>BS50</td><td>F89</td><td>24/10/2019</td><td>Mohammedia</td><td>28D</td><td>V</td><td>-</td><td>N</td></tr><tr><td>BS57</td><td rowspan="2">F90</td><td>05/03/2020</td><td rowspan="2">Casablanca</td><td>34D</td><td rowspan="2">V</td><td>29</td><td rowspan="2">P</td></tr><tr><td>BS58</td><td>05/03/2020</td><td>34D</td><td>-</td></tr><tr><td>BS59</td><td rowspan="2">F91</td><td>05/03/2020</td><td rowspan="2">Casablanca</td><td>30D</td><td rowspan="2">V</td><td>-</td><td rowspan="2">P</td></tr><tr><td>BS60</td><td>05/03/2020</td><td>30D</td><td>33</td></tr><tr><td>BS61</td><td>F92</td><td>02/07/2020</td><td>Tiznit</td><td>42D</td><td>V</td><td>-</td><td>N</td></tr><tr><td>BS62</td><td>F93</td><td>02/07/2020</td><td>Tiznit</td><td>34D</td><td>V</td><td>30</td><td>P</td></tr><tr><td>BS63</td><td>F94</td><td>02/07/2020</td><td>Tiznit</td><td>30D</td><td>V</td><td>30</td><td>P</td></tr><tr><td>BS64</td><td>F95</td><td>02/07/2020</td><td>Tiznit</td><td>22D</td><td>V</td><td>30</td><td>P</td></tr><tr><td>BS65</td><td>F96</td><td>02/07/2020</td><td>Tiznit</td><td>30D</td><td>V</td><td>29</td><td>P</td></tr><tr><td>BS66</td><td>F97</td><td>02/07/2020</td><td>Tiznit</td><td>35D</td><td>V</td><td>32</td><td>P</td></tr><tr><td>BS67</td><td>F98</td><td>02/07/2020</td><td>Tiznit</td><td>29D</td><td>V</td><td>27</td><td>P</td></tr><tr><td>BS68</td><td>F99</td><td>02/07/2020</td><td>Tiznit</td><td>21D</td><td>V</td><td>29</td><td>P</td></tr><tr><td>BS69</td><td>F100</td><td>02/07/2020</td><td>Tiznit</td><td>13D</td><td>V</td><td>27</td><td>P</td></tr><tr><td>BS70</td><td>F101</td><td>02/07/2020</td><td>Tiznit</td><td>31D</td><td>V</td><td>28</td><td>P</td></tr><tr><td>BS71</td><td>F102</td><td>02/07/2020</td><td>Tiznit</td><td>44D</td><td>V</td><td>30</td><td>P</td></tr><tr><td>BS72</td><td>F103</td><td>02/07/2020</td><td>Tiznit</td><td>36D</td><td>V</td><td>28</td><td>P</td></tr><tr><td>BS73</td><td>F104</td><td>02/07/2020</td><td>Tiznit</td><td>24D</td><td>V</td><td>29</td><td>P</td></tr><tr><td>BS76</td><td rowspan="6">F105</td><td>02/07/2020</td><td rowspan="6">Tiznit</td><td>42D</td><td rowspan="6">V</td><td>30</td><td rowspan="6">P</td></tr><tr><td>BS77</td><td>02/07/2020</td><td>42D</td><td>28</td></tr><tr><td>BS78</td><td>02/07/2020</td><td>42D</td><td>29</td></tr><tr><td>BS79</td><td>02/07/2020</td><td>42D</td><td>27</td></tr><tr><td>BS80</td><td>02/07/2020</td><td>42D</td><td>30</td></tr><tr><td>BS81</td><td>02/07/2020</td><td>42D</td><td>32</td></tr><tr><td>BS82</td><td rowspan="3">F105</td><td>02/07/2020</td><td rowspan="3">Tiznit</td><td>42D</td><td rowspan="3">V</td><td>39</td><td rowspan="3">P</td></tr><tr><td>BS83</td><td>02/07/2020</td><td>42D</td><td>39</td></tr><tr><td>BS84</td><td>02/07/2020</td><td>42D</td><td>36</td></tr></tbody></table>
<table><caption>**Table A1. Cont.**</caption><thead><tr><th>Sample</th><th>Farm</th><th>Sampling Date</th><th>Location of Farm</th><th>Age of Birds</th><th>H9N2 Vaccination Status</th><th>H9N2 RT-qPCR Ct Value</th><th>Farms Status</th></tr></thead><tbody><tr><td>BS85</td><td rowspan="2">F106</td><td>09/09/2020</td><td rowspan="2">Casablanca</td><td>33D</td><td rowspan="2">NV</td><td>31</td><td rowspan="2">P</td></tr><tr><td>BS86</td><td>09/09/2020</td><td>33D</td><td>30</td></tr><tr><td>BS87</td><td rowspan="4">F107</td><td>09/09/2020</td><td rowspan="4">Casablanca</td><td>27D</td><td rowspan="4">NV</td><td>30</td><td rowspan="4">P</td></tr><tr><td>BS88</td><td>09/09/2020</td><td>27D</td><td>30</td></tr><tr><td>BS89</td><td>09/09/2020</td><td>27D</td><td>29</td></tr><tr><td>BS90</td><td>09/09/2020</td><td>27D</td><td>30</td></tr><tr><td>BS91 *</td><td rowspan="3">F108</td><td>10/09/2020</td><td rowspan="3">Rabat</td><td>29D</td><td rowspan="3">V</td><td>24</td><td rowspan="3">P</td></tr><tr><td>BS92</td><td>10/09/2020</td><td>29D</td><td>29</td></tr><tr><td>BS93</td><td>10/09/2020</td><td>29D</td><td>34</td></tr></tbody></table>

F: Farm; D: day; V: vaccinated; NV: Non-Vaccinated; P: positive; N: negative; Ct: Cycle Threshold; *: isolated sample; +: sequenced sample.

<table><caption>**Table A2.** HA mutations as compared to 2016 Moroccan strain SF1 (H3 numbering).</caption><thead><tr><th></th><th>137</th><th>188</th><th>190</th><th>222</th><th>226</th><th>298</th><th>325</th><th>364</th><th>375</th><th>397</th><th>402</th><th>408</th><th>493</th></tr></thead><tbody><tr><td>SF1</td><td>T</td><td>D</td><td>A</td><td>L</td><td>L</td><td>I</td><td>H</td><td>M</td><td>V</td><td>D</td><td>E</td><td>N</td><td>T</td></tr><tr><td>AS13</td><td>T</td><td>D</td><td>T</td><td>L</td><td>L</td><td>I</td><td>H</td><td>I</td><td>V</td><td>D</td><td>E</td><td>N</td><td>T</td></tr><tr><td>AS14</td><td>T</td><td>D</td><td>A</td><td>L</td><td>L</td><td>I</td><td>Q</td><td>M</td><td>V</td><td>D</td><td>E</td><td>N</td><td>I</td></tr><tr><td>AS29</td><td>T</td><td>D</td><td>V</td><td>L</td><td>L</td><td>I</td><td>H</td><td>M</td><td>V</td><td>N</td><td>E</td><td>N</td><td>T</td></tr><tr><td>AS32</td><td>T</td><td>D</td><td>A</td><td>L</td><td>L</td><td>V</td><td>H</td><td>M</td><td>V</td><td>D</td><td>D</td><td>N</td><td>T</td></tr><tr><td>AS71</td><td>T</td><td>D</td><td>A</td><td>L</td><td>L</td><td>V</td><td>H</td><td>M</td><td>I</td><td>D</td><td>E</td><td>S</td><td>T</td></tr><tr><td>AS76</td><td>T</td><td>N</td><td>A</td><td>L</td><td>Q</td><td>V</td><td>H</td><td>M</td><td>V</td><td>D</td><td>E</td><td>N</td><td>T</td></tr><tr><td>AS77</td><td>T</td><td>D</td><td>A</td><td>L</td><td>Q</td><td>V</td><td>H</td><td>M</td><td>V</td><td>D</td><td>E</td><td>N</td><td>T</td></tr></tbody></table>

<table><caption>**Table A3.** NA mutations as compared to 2016 Moroccan strain SF1 (N2 numbering).</caption><thead><tr><th></th><th>8</th><th>16</th><th>31</th><th>42</th><th>46</th><th>50</th><th>56</th><th>57</th><th>58</th><th>60</th><th>65</th><th>80</th><th>88</th><th>101</th><th>116</th><th>127</th><th>254</th><th>261</th><th>290</th><th>329</th><th>332</th><th>368</th><th>385</th><th>390</th><th>400</th><th>416</th></tr></thead><tbody><tr><td>SF1</td><td>I</td><td>T</td><td>T</td><td>Y</td><td>S</td><td>A</td><td>I</td><td>I</td><td>I</td><td>R</td><td>I</td><td>T</td><td>S</td><td>S</td><td>V</td><td>S</td><td>I</td><td>K</td><td>V</td><td>N</td><td>S</td><td>K</td><td>T</td><td>S</td><td>N</td><td>I</td></tr><tr><td>AS13</td><td>I</td><td>T</td><td>T</td><td>Y</td><td>S</td><td>A</td><td>I</td><td>I</td><td>I</td><td>R</td><td>I</td><td>T</td><td>S</td><td>S</td><td>V</td><td>S</td><td>I</td><td>R</td><td>V</td><td>D</td><td>S</td><td>K</td><td>N</td><td>A</td><td>S</td><td>I</td></tr><tr><td>AS14</td><td>I</td><td>T</td><td>T</td><td>Y</td><td>S</td><td>T</td><td>T</td><td>T</td><td>T</td><td>R</td><td>I</td><td>P</td><td>L</td><td>S</td><td>V</td><td>N</td><td>V</td><td>K</td><td>V</td><td>N</td><td>S</td><td>T</td><td>T</td><td>S</td><td>N</td><td>M</td></tr><tr><td>AS29</td><td>I</td><td>T</td><td>T</td><td>H</td><td>P</td><td>A</td><td>I</td><td>I</td><td>I</td><td>K</td><td>I</td><td>T</td><td>S</td><td>S</td><td>V</td><td>S</td><td>I</td><td>R</td><td>V</td><td>N</td><td>S</td><td>K</td><td>T</td><td>A</td><td>S</td><td>I</td></tr><tr><td>AS32</td><td>I</td><td>I</td><td>T</td><td>Y</td><td>S</td><td>A</td><td>I</td><td>I</td><td>I</td><td>R</td><td>T</td><td>T</td><td>S</td><td>S</td><td>I</td><td>S</td><td>I</td><td>K</td><td>V</td><td>N</td><td>S</td><td>K</td><td>T</td><td>S</td><td>S</td><td>I</td></tr><tr><td>AS71</td><td>M</td><td>I</td><td>T</td><td>Y</td><td>S</td><td>A</td><td>I</td><td>I</td><td>I</td><td>R</td><td>T</td><td>T</td><td>S</td><td>A</td><td>V</td><td>S</td><td>I</td><td>K</td><td>A</td><td>N</td><td>F</td><td>K</td><td>T</td><td>S</td><td>S</td><td>I</td></tr><tr><td>AS76</td><td>M</td><td>I</td><td>T</td><td>Y</td><td>S</td><td>A</td><td>I</td><td>I</td><td>I</td><td>R</td><td>T</td><td>T</td><td>S</td><td>A</td><td>V</td><td>S</td><td>I</td><td>K</td><td>V</td><td>N</td><td>F</td><td>K</td><td>T</td><td>S</td><td>S</td><td>I</td></tr><tr><td>AS77</td><td>I</td><td>I</td><td>M</td><td>Y</td><td>S</td><td>A</td><td>I</td><td>I</td><td>I</td><td>R</td><td>T</td><td>T</td><td>S</td><td>S</td><td>I</td><td>S</td><td>I</td><td>K</td><td>V</td><td>N</td><td>S</td><td>K</td><td>T</td><td>S</td><td>S</td><td>I</td></tr></tbody></table>

## References

1. Taxonomy. Available online: https://talk.ictvonline.org/taxonomy/ (accessed on 3 August 2021).

2. Capua, I.; Alexander, D.J. Avian Influenza: Recent Developments. *Avian Pathol.* **2004**, *33*, 393–404. [CrossRef] [PubMed]

3. Tong, S.; Zhu, X.; Li, Y.; Shi, M.; Zhang, J.; Bourgeois, M.; Yang, H.; Chen, X.; Recuenco, S.; Gomez, J.; et al. New World Bats Harbor Diverse Influenza A Viruses. *PLoS Pathog.* **2013**, *9*, e1003657. [CrossRef] [PubMed]

4. Azizpour, A.; Goudarzi, H.; Charkhkar, S.; Momayez, R.; Hablolvarid, M.H. Experimental Study on Tissue Tropism and Dissemination of H9N2 Avian Influenza Virus and Ornithobacterium Rhinotracheale Co-Infection in SPF Chickens. *J. Anim. Plant Sci.* **2014**, *24*, 1655–1662.

5. Seifi, S.; Asasi, K.; Mohammadi, A. Short Paper: An Experimental Study on Broiler Chicken Co-Infected with the Specimens Containing Avian Influenza (H9 Subtype) and Infectious Bronchitis (4/91 Strain) Viruses. *Iran. J. Vet. Res.* **2012**, *13*, 5.

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

7. Chen, B.L.; Zhang, Z.J.; Chen, W.B. Isolation and Preliminary Serological Characterization of Type A Influenza Viruses from Chickens. *Chin. J. Vet. Med.* **1994**, *22*, 3–5.

8. Lee, D.-H.; Swayne, D.E.; Sharma, P.; Rehmani, S.F.; Wajid, A.; Suarez, D.L.; Afonso, C.L. H9N2 Low Pathogenic Avian Influenza in Pakistan (2012–2015). *Vet. Rec. Open* **2016**, *3*, e000171. [CrossRef]

9. Pan, Q.; Liu, A.; Zhang, F.; Ling, Y.; Ou, C.; Hou, N.; He, C. Co-Infection of Broilers with Ornithobacterium Rhinotracheale and H9N2 Avian Influenza Virus. *BMC Vet. Res.* **2012**, *8*, 104. [CrossRef]

10. Nili, H.; Asasi, K. Avian Influenza (H9N2) Outbreak in Iran. *Avian Dis.* **2003**, *47*, 828–831. [CrossRef]

11. Monne, I.; Fusaro, A.; Al-Blowi, M.H.; Ismail, M.M.; Khan, O.A.; Dauphin, G.; Tripodi, A.; Salviato, A.; Marangon, S.; Capua, I.; et al. Co-Circulation of Two Sublineages of HPAI H5N1 Virus in the Kingdom of Saudi Arabia with Unique Molecular Signatures Suggesting Separate Introductions into the Commercial Poultry and Falconry Sectors. *J. Gen. Virol.* **2008**, *89*, 2691–2697. [CrossRef]
12. Roussan, D.A.; Khawaldeh, G.Y.; Al Rifai, R.H.; Totanji, W.S.; Shaheen, I.A. Avian Influenza Virus H9 Subtype in Poultry Flocks in Jordan. *Prev. Vet. Med.* **2009**, *88*, 77–81. [CrossRef] [PubMed]

13. Wernery, U.; Shanmuganatham, K.K.; Krylov, P.S.; Joseph, S.; Friedman, K.; Krauss, S.; Webster, R.G. H9N2 Influenza Viruses from Birds Used in Falconry. *Influenza Other Respir. Viruses* **2013**, *7*, 1241–1245. [CrossRef] [PubMed]

14. Tombari, W.; Paul, M.; Bettaieb, J.; Larbi, I.; Nsiri, J.; Elbehi, I.; Gribaa, L.; Ghram, A. Risk Factors and Characteristics of Low Pathogenic Avian Influenza Virus Isolated from Commercial Poultry in Tunisia. *PLoS ONE* **2013**, *8*, e53524. [CrossRef] [PubMed]

15. Kandeil, A.; El-Shesheny, R.; Maatouq, A.M.; Moatasim, Y.; Shehata, M.M.; Bagato, O.; Rubrum, A.; Shanmuganatham, K.; Webby, R.J.; Ali, M.A.; et al. Genetic and Antigenic Evolution of H9N2 Avian Influenza Viruses Circulating in Egypt between 2011 and 2013. *Arch. Virol.* **2014**, *159*, 2861–2876. [CrossRef]

16. Body, M.H.; Alrarawahi, A.H.; Alhubsy, S.S.; Saravanan, N.; Rajmony, S.; Mansoor, M.K. Characterization of Low Pathogenic Avian Influenza Virus Subtype H9N2 Isolated from Free-Living Mynah Birds (*Acridotheres tristis*) in the Sultanate of Oman. *Avian Dis.* **2015**, *59*, 329–334. [CrossRef]

17. Kammon, A.; Heidari, A.; Dayhum, A.; Eldaghayes, I.; Sharif, M.; Monne, I.; Cattoli, G.; Asheg, A.; Farhat, M.; Kraim, E. Characterization of Avian Influenza and Newcastle Disease Viruses from Poultry in Libya. *Avian Dis.* **2015**, *59*, 422–430. [CrossRef]

18. Group, T.S.H.W.; Schultz-Cherry, S.; Thomas, P. Assessing the Fitness of Distinct Clades of Influenza A (H9N2) Viruses. *Emerg. Microbes Infect.* **2013**, *2*, e75. [CrossRef]

19. Jestin, V.; Manuguerra, J.C.; Eterradossi, N. Risque de Transmission à l'homme Des Virus Influenza Aviaires. *Bull. Epidemiol. AFSSA* **2003**, *11*, 1–3.

20. Iqbal, M.; Yaqub, T.; Mukhtar, N.; Shabbir, M.Z.; McCauley, J.W. Infectivity and Transmissibility of H9N2 Avian Influenza Virus in Chickens and Wild Terrestrial Birds. *Vet. Res.* **2013**, *44*, 100. [CrossRef]

21. Zhang, K.; Zhang, Z.; Yu, Z.; Li, L.; Cheng, K.; Wang, T.; Huang, G.; Yang, S.; Zhao, Y.; Feng, N.; et al. Domestic Cats and Dogs Are Susceptible to H9N2 Avian Influenza Virus. *Virus Res.* **2013**, *175*, 52–57. [CrossRef]

22. Peiris, M.; Yam, W.C.; Chan, K.H.; Ghose, P.; Shortridge, K.F. Influenza A H9N2: Aspects of Laboratory Diagnosis. *J. Clin. Microbiol.* **1999**, *37*, 3426–3427. [CrossRef] [PubMed]

23. Peiris, M.; Yuen, K.; Leung, C.; Chan, K.; Ip, P.; Lai, R.; Orr, W.; Shortridge, K. Human Infection with Influenza H9N2. *Lancet* **1999**, *354*, 916–917. [CrossRef]

24. Guo, Y.; Li, J.; Cheng, X. Discovery of Men Infected by Avian Influenza A (H9N2) Virus. *Zhonghua Shi Yan He Lin Chuang Bing Du Xue Za Zhi* **1999**, *13*, 105–108. [PubMed]

25. Lin, Y.P.; Shaw, M.; Gregory, V.; Cameron, K.; Lim, W.; Klimov, A.; Subbarao, K.; Guan, Y.; Krauss, S.; Shortridge, K.; et al. Avian-to-Human Transmission of H9N2 Subtype Influenza A Viruses: Relationship between H9N2 and H5N1 Human Isolates. *Proc. Natl. Acad. Sci. USA* **2000**, *97*, 9654–9658. [CrossRef] [PubMed]

26. EL Houadfi, M.; Fellahi, S.; Nassik, S.; Guérin, J.-L.; Ducatez, M.F. First Outbreaks and Phylogenetic Analyses of Avian Influenza H9N2 Viruses Isolated from Poultry Flocks in Morocco. *Virol. J.* **2016**, *13*, 140. [CrossRef]

27. Monne, I.; Ormelli, S.; Salviato, A.; Battisti, C.D.; Bettini, F.; Salomoni, A.; Drago, A.; Zecchin, B.; Capua, I.; Cattoli, G. Development and Validation of a One-Step Real-Time PCR Assay for Simultaneous Detection of Subtype H5, H7, and H9 Avian Influenza Viruses. *J. Clin. Microbiol.* **2008**, *46*, 5. [CrossRef] [PubMed]

28. Barman, S.; Turner, J.C.M.; Hasan, M.K.; Akhtar, S.; El-Shesheny, R.; Franks, J.; Walker, D.; Seiler, P.; Friedman, K.; Kercher, L.; et al. Continuing Evolution of Highly Pathogenic H5N1 Viruses in Bangladeshi Live Poultry Markets. *Emerg. Microbes Infect.* **2019**, *8*, 650–661. [CrossRef]

29. Hall, T.A. BioEdit: A User-Friendly Biological Sequence Alignment Editor and Analysis Program for Windows 95/98/NT. *Nucleic Acids Symp. Ser.* **1999**, *41*, 95–98.

30. Thompson, J.D.; Higgins, D.G.; Gibson, T.J. CLUSTAL W: Improving the Sensitivity of Progressive Multiple Sequence Alignment through Sequence Weighting, Position-Specific Gap Penalties and Weight Matrix Choice. *Nucleic Acids Res.* **1994**, *22*, 4673–4680. [CrossRef]

31. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. *Mol. Biol. Evol.* **2013**, *30*, 2725–2729. [CrossRef]

32. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic Local Alignment Search Tool. *J. Mol. Biol.* **1990**, *215*, 403–410. [CrossRef]

33. Peacock, T.P.; Benton, D.J.; Sadeyen, J.-R.; Chang, P.; Sealy, J.E.; Bryant, J.E.; Martin, S.R.; Shelton, H.; McCauley, J.W.; Barclay, W.S.; et al. Variability in H9N2 Haemagglutinin Receptor-Binding Preference and the PH of Fusion. *Emerg. Microbes Infect.* **2017**, *6*, 1–7. [CrossRef] [PubMed]

34. 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] [PubMed]

35. Li, X.; Shi, J.; Guo, J.; Deng, G.; Zhang, Q.; Wang, J.; He, X.; Wang, K.; Chen, J.; Li, Y.; et al. Genetics, Receptor Binding Property, and Transmissibility in Mammals of Naturally Isolated H9N2 Avian Influenza Viruses. *PLoS Pathog.* **2014**, *10*, e1004508. [CrossRef]
36. Teng, Q.; Xu, D.; Shen, W.; Liu, Q.; Rong, G.; Li, X.; Yan, L.; Yang, J.; Chen, H.; Yu, H.; et al. A Single Mutation at Position 190 in Hemagglutinin Enhances Binding Affinity for Human Type Sialic Acid Receptor and Replication of H9N2 Avian Influenza Virus in Mice. *J. Virol.* **2016**, *90*, 9806–9825. [CrossRef]

37. Hurt, A.C.; Holien, J.K.; Parker, M.W.; Barr, I.G. Oseltamivir Resistance and the H274Y Neuraminidase Mutation in Seasonal, Pandemic and Highly Pathogenic Influenza Viruses. *Drugs* **2009**, *69*, 2523–2531. [CrossRef]

38. Suttie, A.; Deng, Y.-M.; Greenhill, A.R.; Dussart, P.; Horwood, P.F.; Karlsson, E.A. Inventory of Molecular Markers Affecting Biological Characteristics of Avian Influenza A Viruses. *Virus Genes* **2019**, *55*, 739–768. [CrossRef]

39. Jonas, M.; Sahesti, A.; Murwijati, T.; Lestariningsih, C.L.; Irine, I.; Ayesda, C.S.; Prihartini, W.; Mahardika, G.N. Identification of Avian Influenza Virus Subtype H9N2 in Chicken Farms in Indonesia. *Prev. Vet. Med.* **2018**, *159*, 99–105. [CrossRef]

40. El Khantour, A.; El Houadfi, M.; Saâdia, N.; Tligui, N.; el Mellouli, F.; Sikht, F.-Z.; Ducatez, M.; Soulaymani, A.; Fellahi, S. Protective Efficacy Evaluation of Four Inactivated Commercial Vaccines Against Low Pathogenic Avian Influenza H9N2 Virus Under Experimental Conditions in Broiler Chickens. *Avian Dis.* **2021**, *65*, 351–357. [CrossRef]

41. Park, K.J.; Kwon, H.; Song, M.-S.; Pascua, P.N.Q.; Baek, Y.H.; Lee, J.H.; Jang, H.-L.; Lim, J.-Y.; Mo, I.-P.; Moon, H.-J.; et al. Rapid Evolution of Low-Pathogenic H9N2 Avian Influenza Viruses Following Poultry Vaccination Programmes. *J. Gen. Virol.* **2011**, *92*, 36–50. [CrossRef]

42. Alexander, D.J. An Overview of the Epidemiology of Avian Influenza. *Vaccine* **2007**, *25*, 5637–5644. [CrossRef] [PubMed]

43. Van den Berg, T.; Lambrecht, B.; Marché, S.; Steensels, M.; Van Borm, S.; Bublot, M. Influenza Vaccines and Vaccination Strategies in Birds. *Comp. Immunol. Microbiol. Infect. Dis.* **2008**, *31*, 121–165. [CrossRef] [PubMed]

44. Choi, J.G.; Lee, Y.J.; Kim, Y.J.; Lee, E.K.; Jeong, O.M.; Sung, H.W.; Kim, J.H.; Kwon, J.H. An Inactivated Vaccine to Control the Current H9N2 Low Pathogenic Avian Influenza in Korea. *J. Vet. Sci.* **2008**, *9*, 67. [CrossRef]

45. 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. *JVI* **2007**, *81*, 5181–5191. [CrossRef]

46. Kawaoka, Y.; Webster, R.G. Sequence Requirements for Cleavage Activation of Influenza Virus Hemagglutinin Expressed in Mammalian Cells. *Proc. Natl. Acad. Sci. USA* **1988**, *85*, 324–328. [CrossRef] [PubMed]