1. Introduction
The H9N2 AIV was initially isolated from turkeys in Wisconsin in the United States of America in 1966 (A/turkey/Wisconsin/1/1966(H9N2)) [1]. In Asia, the first cases of H9N2 infection in poultry were recorded in Hong Kong in the 1970s [2]. In the 1990s, the H9N2 spread became more extensive with several outbreaks affecting industrial poultry in China, Germany, South Korea, Italy, Iran, Pakistan, the USA, South Africa, and Ireland [3,4,5]. This lead to a persistent viral circulation in many countries in Asia, the Middle East, and North Africa [6], where the virus has attained endemic status in poultry populations [7,8,9]. Today, H9N2 is considered the most prevalent, and damaging, low-pathogenic avian influenza (LPAI) subtype in poultry worldwide [7], presenting a threat to both animal and human health [9,10]. H9N2 infections usually cause minor clinical disease in poultry, which makes its identification in the field rather challenging [11]. However, despite their low pathogenicity, H9N2 viruses may cause variable pathogenicity and mortality in domestic poultry [5,6,12], especially when associated with immunosuppression or coinfection with other viruses or bacteria [3,11,13,14], resulting in significant economic losses [5,6,12].
The H9N2 virus has a zoonotic potential with interspecies transmission between poultry and mammals, including swine and humans [15,16,17,18]. At least 59 human cases have been reported to date [9]. Furthermore, H9N2 viruses have proved able to provide partial or even complete sets of internal gene segments to zoonotic avian influenza viruses (AIVs), leading to the emergence of human-lethal reassortants, as in the case of China in 2013 with the appearance of H7N9 and H10N8 viruses, and, a year later of the H5N6 subtype [19,20,21,22,23,24,25,26] posing serious risks to public health. Moreover, there are multiple reports of H9N2 natural infection in other mammals, such as ferrets [27], mice [28], minks, foxes, raccoon dogs [29], cats, and dogs [30].
Based on their genetic characteristics, the H9N2 AIVs are divided into three major lineages: the G1-like lineage (A/Quail/Hong Kong/G1/97), the Korean-like lineage (A/Chicken/Korea/38349-p96323/96 and A/Duck/Hong Kong/Y439/97) and the Y280-like lineage (A/Duck/Hong Kong/Y280/97, A/Chicken/Beijing/1/94, and A/Chicken/Hong Kong /G9/97) [3,4,5,31]. Until now, the G1-like lineage has been the only one found in Africa [32].
Since the early 2000s, the emergence of H9N2 viruses has been reported in several North African countries, including Libya in 2006 [33], Egypt in 2006 [34], and Tunisia in 2009 [35]. Additionally, following the first isolation of the LPAI H9N2 virus in Morocco in January 2016 [36], H9N2 viruses have been extensively circulating in North Africa, spreading to West Africa, as well as to East Africa. Since 2017, H9N2 viruses have been detected in Algeria [37] and in several Sub-Saharan African countries, including Burkina Faso [38], Ghana [32], Uganda, Kenya, Senegal [39], Benin Republic, and Togo [40]. All H9N2 viruses detected in poultry in Africa belong to the G1 lineage and are related to those isolated in the United Arab Emirates.
In Morocco, following its first identification in 2016, the H9N2 virus has spread rapidly across the country, causing significant economic losses in poultry [36]. A previous study revealed that the Moroccan isolate might be related to the Middle East isolate A/chicken/Dubai/D2506.A/2015 [41]. The virus showed high pathogenicity under field conditions with aggravation of the disease in case of coinfection with other respiratory diseases [41,42]. Despite the swift implementation of a mass poultry vaccination program following the onset of the disease, the virus became endemic in all poultry sectors of the country [43]. Knowledge of the evolutionary history and the circulation of the Moroccan H9N2 viruses is crucial for public and animal health in Morocco. A better understanding of the epidemiological and geographic origins of the Moroccan H9N2 viruses might help to improve the control and management of future outbreaks. In this study, we characterized the whole genome of H9N2 AIVs collected in Morocco from 2016 to 2021, exploring the spatiotemporal dynamics and the acquisition of mutations of particular concern for animal and human health.
2. Materials and Methods
2.1. Sample Collection
From 2016 to 2021, thanks to the establishment of a program to monitor the presence of influenza viruses in poultry in Morocco, samples were collected from 81 poultry farms located in seven different regions in the country, including layers, breeders, turkeys, and broiler flocks of varying ages. Most flocks (72.84%, 59/81) were vaccinated against AIV H9N2. Oropharyngeal, cloacal swabs and organ samples were collected from chickens and turkeys, stored at 2–8 °C, and sent to the laboratory within 24 h for further analysis.
2.2. Real Time RT-PCR for H9 Gene Detection
Swabs, as well as tissue samples, were re-suspended in 1 mL of phosphate-buffered saline (pH 7.4 +/− 0.2, supplemented with 10,000 IU/mL penicillin, 10 mg/mL streptomycin, 0.25 mg/mL gentamycin, and 5000 IU/mL nystatin), thoroughly homogenized by using Bead Blaster 24 Homogenizer (Benchmark Scientific, 2600, Sayreville, NJ, USA). After centrifugation at 1500× g for 10 min at 4 °C, total RNA was purified with the NucleoSpin® RNA Virus Kit (Macherey-Nagel, Düren, Germany) following the manufacturer’s instructions. RNA was then eluted in 50 µL RNase-free water. All samples were tested by RT PCR for the influenza A matrix gene as described by Spackman, E. et al. (2002) [44]. The positive samples were then tested for the H9 gene detection by RT-PCR according to Monne, I. et al. (2008) [45] by using the AgPath-IDTM One-Step RT-PCR Kit (Ambion, Applied Biosystems, Grand Island, NY, USA). The rRT-PCR was performed with Aria Mx (Agilent Technologies, Santa Clara, CA, USA) and ABI 7500 Fast (Applied Biosystems, Foster City, CA, USA) 96-well plate real-time PCR devices.
2.3. Whole Genome Sequencing
The complete genomes were amplified using SuperScript™ III One-Step RT-PCR System with Platinum™ Taq High Fidelity DNA Polymerase (Invitrogen, Carlsbard, CA, USA), as previously described by Zhou B. et al. (2009) [46]. Amplicons were purified by using Agencourt AMPure XP (Beckman Coulter Inc., Brea, CA, USA) and quantified with Qubit™ DNA HS Assay (Thermo Fisher Scientific, Waltham, MA, USA). Sequencing libraries were prepared with the Illumina Nextera XT DNA Sample Preparation Kit (Illumina, San Diego, CA, USA) and sequenced on the Illumina MiSeq platform (2 × 250 bp Paired-End; Illumina, San Diego, CA, USA).
Read quality was assessed using FastQC v0.11.2 (
2.4. Phylogenetic and Evolutionary Analyses
After sequence alignment using MAFFT v7 [52], IQTREE v1.6 was used to perform ultrafast bootstrap resampling analyses (1000 replications) to generate the maximum likelihood phylogenetic trees [53,54]. Trees were visualized with FigTree v1.4.4 (
Phylogenetic network analyses using the Median Joining method implemented in the program NETWORK 10 (
The evolutionary rate and the time to the most common recent ancestor (tMRCA) were estimated for all the gene segments using BEAST v1.10.4 [56]. Markov chain Monte Carlo (MCMC) sampling analyses employing an uncorrelated lognormal relaxed molecular clock, that allows for rate variation across lineages, and the SRD06 substitution model (HKY85 + Γ4 with two partitions—1st + 2nd positions vs. 3rd position -, base frequencies, and Γ-rate heterogeneity unlinked across all codon positions) were implemented [57]. Tree Annotator v1.10.4 (
2.5. Analysis of Selection Pressures
Site-specific selection pressures for the HA gene segment of the Moroccan H9N2 viruses were measured as the ratio of nonsynonymous (dN) to synonymous (dS) nucleotide substitutions per site. Specifically, the Mixed Effects Model of Evolution (MEME) and the fixed-effects likelihood (FEL) methods were used to infer sites subjected to episodic or pervasive positive selection [58,59]. Analyses were performed using the Datamonkey online version of the Hy-Phy package [60]. All analyses used the GTR model of nucleotide substitution and ML phylogenetic trees.
3. Results
3.1. Phylogenetic and Evolutionary Analyses
In the framework of the surveillance activity implemented in Morocco between 2016 and 2021, field samples collected from 81 poultry farms were submitted to the laboratory to be tested for the presence of AIV. Overall, 76 out of 81 (76/81) tested positive for AIV and specifically for the H9 subtype (Figure 1 shows the location of the sampled and the positive farms). In total, 28 representative positive samples were selected based on the year and area of collection for whole genome sequencing. Specifically, we obtained the complete genome of 22 viruses and the partial genome of 6 viruses (Table 1).
Phylogenetic analysis of the hemagglutinin (HA) gene segment revealed that all the H9N2 viruses collected in the country since 2016 belong to the G1 lineage and cluster with H9N2 viruses identified in the Middle East (United Arab Emirates and Oman), North Africa (Algeria), and West Africa (Burkina Faso, Senegal, Togo, Ghana, Benin) between 2015 and 2020 (Figure 2). The phylogenetic trees of the other gene segments confirm the same genetic clustering (see Supplementary Materials Figures S1–S7). The tree topologies indicate that the H9N2 viruses from Morocco had evolved into multiple genetic groups within the country. However, the lack of a consistency of the groups among the different phylogenies is suggestive of the occurrence of multiple inter-subtype reassortment events.
Virus evolution into multiple genetic clusters is evident also from the HA median-joining network analysis. Although all the 2016–2017 sequences (yellow in Figure 3) show a high genetic relationship (from 1 to 16 nucleotide differences) and cluster together on the centre of the network, the viruses collected in 2018–2019 (green in Figure 3) and 2020–2021 (blue in Figure 3) separate from the 2016–2017 progenitor viruses by long branches (from 15 to 54 nucleotide differences) that spread from the centre in multiple directions, forming separate clusters. The high number of median vectors, representing the lost ancestral sequences, clearly indicate a relevant data gap. No particular grouping was observed according to the geographical origin, which suggests a wide virus dispersal within the country.
To investigate the role of Morocco in the geographic spread of the H9N2 viruses of the G1 lineage, we performed a Bayesian phylogeographic analysis. Our results showed a virus spread from the Middle East to Morocco, and subsequently from Morocco to Algeria and from Morocco to West Africa. All of these transitions were well supported by Bayes Factors higher than 20 (BF > 20, Figure 4). This finding suggests that Morocco might have played a central role in the spread of the H9N2 virus in the African continent, although this interpretation may have been biased by the higher number of sequences from Morocco compared to the ones from other African countries.
3.2. Molecular Analysis and Positively Selected Sites
The molecular analysis revealed that the H9N2 viruses from Morocco present the same cleavage site (RSSR/GLF) as the previously identified North African H9N2 viruses of the G1 lineage. In addition, most of the Moroccan viruses possess the amino acid Leucine at position 226 and Threonine at position 155 (H3 numbering) of the HA receptor binding site; these mutations imply preferential binding to human-like α2-6-linked sialic acid (SA α2–6) receptors [61,62]. Of note, site 226 has been identified as under positive selection by both MEME and FEL methods (p-value < 0.05) (Table 2). The H9N2 Moroccan viruses show three different amino acids in this position: Leucine (L), Glutamine (Q), and Arginine (R). These amino acid substitutions may affect not only the virus receptor binding properties but also influence its antigenic specificity [63]. Although amino acids L and Q have frequently been identified in viruses belonging to the G1 lineage, R has rarely been observed in the H9N2 viruses (one virus from Algeria and one from Morocco). Previous studies demonstrated that this mutation is associated with change of the pandemic 2009 (H1N1) viral HA binding preference, from the human-type receptor, α2,6-linked sialic acid, to the avian-type receptor, α2,3-linked sialic acid [64,65].
In addition to 226, in the HA gene we identified three further sites (sites 190, 160, and 186) under positive selection: one (site 190) using both methods, while the other two were recognized only by FEL. Of note, all these sites fall within the receptor-binding domain and/or antigenic sites [63,66,67].
4. Discussion
In Morocco, the first outbreak of LPAI H9N2 was described in broiler flocks in January 2016 [36]. The virus was later isolated from layers and breeders in different regions of the country. As part of a surveillance program, the national veterinary services implemented mass vaccination against H9N2 using homologous inactivated vaccines in all sectors of poultry production, first for breeders and layers, and a few days later for turkey and broiler flocks. The outbreak was controlled in 3–4 months, after which the vaccination was authorized for all types of poultry production using autogenous and heterologous inactivated vaccines. Over the past 5 years, 13 H9N2 vaccines originating from different H9 strains, i.e., from Morocco, UAE, Korea, Turkey, Italy, Egypt, were available and commercialized in Morocco (see Supplementary Materials Table S2). Despite this vaccine availability, a recent study on the effectiveness of the four most used vaccines in Moroccan poultry farms has confirmed the need to continuously adapt the vaccine strains to the virus circulating in the field [69]. In the neighboring North African countries, vaccination against AIV H9N2 is also implemented in poultry production sectors using different vaccine types, such as inactivated homologous vaccines and vectored vaccines in Algeria. In Tunisia, vaccination using inactivated homologous vaccines for chicken and turkey breeders and layers has been performed since 2015, but to the present day, no vaccination has been used in broiler flocks. In Egypt, vaccination is currently implemented using local and non-local strains, with variable efficacy. In contrast, vaccination in the field has never been reported in Libya [70]. Despite such vaccination approaches, LPAI H9N2 has become enzootic in North Africa and continues to cause economic losses even in the vaccinated flocks. Thus, circulating H9N2 viruses in North Africa need to be further investigated, and at the same time, revision of employed vaccines and of regional surveillance activities should be implemented.
Our results indicate that all Moroccan H9N2 viruses from 2016 to 2021 belong to the G1 like-lineage and are related to the viruses from the Middle East and North and West Africa, as described in previous reports [43,71]. Results obtained from the spatial analysis demonstrated that AIV H9N2 had been introduced into Morocco from the Middle East, probably through the import of falcons for hunting, and has been circulating in the country for the last six years. The study also showed that the H9N2 G1 like-lineage is the predominant H9 genotype in Moroccan poultry farms. From Morocco, the virus spread to Algeria and West Africa, suggesting the key role of Morocco in the H9N2 virus dispersal on the continent, as reported in previous studies [32,37,38,39,40]. Such virus-spread might have occurred through the movement of poultry professionals [38] and/or the trade of poultry products, since the export market of day-old-chicks and hatching eggs from Morocco to West Africa is well documented [72]. However, the gap in epidemiological data, as well as the limited number of sequences available from the other African countries, might have affected our results and biased our interpretation. Evolutionary and network analyses of the H9N2 viruses from Morocco indicate a high similarity among all the 2016–2017 HA sequences, which is in accordance with previous studies [43,71]. Long branches separate the 2018–2019 and 2020–2021 viruses from their ancestors (2016–2017). This might be explained by the implementation of vaccination from 2016 onwards. The high vaccination coverage and the use of multiple types of vaccines may have driven the virus evolution in multiple directions, with the consequent generation of multiple separate genetic groups in the country. This was further confirmed by the identification of several sites in the HA protein under positive selection located in the RBS or antigenic sites, which may facilitate vaccine escapes, considering the fact that most (71.05%, 54/76) of the H9N2 positive flocks were vaccinated.
As previously observed in the H9N2 African viruses [32,38,40,73], most of Moroccan H9N2 HA protein harbored mammalian-adaptive mutations, including the hemagglutinin human-like marker 216L (226 in H3 numbering), which promotes preferential binding to the α2–6-linked sialic acid receptors of human epithelial cells [60] and was described as an immune escape mechanism [62]. The evidence for positive selection of the 226L marker during serial passage of the virus in chickens was provided by Jegede et al., 2018 [74], suggesting that the H9N2 viruses acquired the mammalian adaptation while circulating in poultry. However, specific receptor-binding studies will be instrumental to assess the zoonotic potential of the H9N2 Moroccan viruses, as other mutations in the HA and internal proteins may be necessary for a successful avian to human adaptation. The present work is the first comprehensive study that provides an overview of the historical origins and evolutionary features of Moroccan H9N2 AIV, based on the Bayesian phylodynamic analysis and whole genome characterization. Over the past five years, the H9N2 AIV in Morocco has shown a rapid genetic diversification, likely directed by selective pressure exerted by vaccine-induced antibodies. The observed mutations and reassortments might extend the range of host species and could lead to the emergence of new strains [75,76]. In addition, our study demonstrates that viruses with a marker well-described in literature associated with zoonotic potential (226L, H3 numbering) may represent a serious risk for public health in Morocco. Hence, the circulation of H9N2 AIV should be investigated in other species, such as wild birds and mammals, including humans, to better understand virus ecology in Morocco. All the more so as a recent molecular study on the presence of the AIV in wild birds staying in the corridor of wetlands during their winter crossing in Morocco, showed a high positivity-rate of the influenza virus in different birds species [77]. Continuous systematic surveillance of molecular evolution of H9N2 AIV is strongly needed to anticipate the emergence of novel H9N2 strains escaping vaccine protection or with an increased human pandemic potential.
Conceptualization, F.E.M., A.F. and A.B.; Data curation, F.E.M. and B.Z.; Formal analysis, F.E.M., B.Z. and H.Z.; Funding acquisition, A.F. and I.M.; Investigation, M.M., H.Z., A.E.K. and H.R.I.; Methodology, F.E.M., A.F., H.Z. and E.P.; Project administration, H.R.I. and A.B.; Resources, F.E.M., M.M. and A.F.; Software, B.Z. and E.G.; Supervision, A.F. and A.B.; Validation, F.E.M. and A.F.; Visualization, F.E.M.; Writing—original draft, F.E.M., A.F. and B.Z.; Writing—review and editing, M.M., A.F., B.Z. and I.M. All authors have read and agreed to the published version of the manuscript.
This research received no external funding.
Not applicable.
Not applicable.
The primary data used to support the findings of this study are available from the corresponding author upon request.
We acknowledge the International Atomic Energy Agency (IAEA) for providing the funding for Fatiha El Mellouli’s training on Sequencing and Phylogeny of Influenza viruses at the IZSVE, Padova, Italy. We are grateful to Mounir Khayli from the Epidemiology and Health Surveillance Unit (SEVS/ONSSA), Rabat and Mohamed Rida Salam from the Immunology and Biodiversity Laboratory, Faculty of Science Ain chock, Hassan II University of Casablanca for their technical support and we also thank Francesca Ellero from the Istituto Zooprofilattico Sperimentale delle Venezie, Padova, Italy for editing the manuscript.
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure 2. Maximum Likelihood phylogenetic tree of the HA gene (IQ-TREE v.1.6.8). The H9N2 viruses from Morocco analyzed in this study are marked in blue.Ultra-fast bootstrap supports equal to or higher than 80% are indicated next to the nodes.
Figure 3. The HA median-joining network showing high genetic relationship between all 2016 and 2017 H9N2 Moroccan viruses while the viruses of 2018–2019 and 2020–2021 are separated from their 2016–2017 ancestors by long branches.
Figure 4. Spatiotemporal dynamics describing the geographic origins of Moroccan H9N2 viruses and the virus spread from Morocco to other African countries. Supported transitions are indicated by continuous lines (BF > 20).
H9N2 isolates related information and GISAID accession numbers.
Year | Virus | Subtype | Flock | Age (d/w 1) | Collection Date | Region | Location | GISAID Accession Number |
---|---|---|---|---|---|---|---|---|
2016 | A/chicken/Morocco/19RS1944-1/2016 | H9N2 | Breeders | na | 25 January 2016 | Fès-Meknès | Fès | EPI1950238–EPI1950245 |
A/chicken/Morocco/19RS1944-2/2016 | H9N2 | Broilers | na | 26 January 2016 | Rabat-Salé-Kénitra | Kenitra | EPI1950325–EPI1950332 | |
2019 | A/pheasant/Morocco/19RS1944-13/2019 | H9N2 | na | na | 19 March 2019 | Casablanca-Settat | Ain Borja | EPI1950209–212, 214-216 and EPI1950384 |
A/pheasant/Morocco/19RS1944-14/2019 | H9N2 | na | na | 19 March 2019 | Casablanca-Settat | Ain Borja | EPI1950217–EPI1950224 | |
A/chicken/Morocco/1_456_21RS1333-1/2019 | H9N2 | Layers | 26 w | 16 February 2019 | Rabat-Salé-Kénitra | Skhirat Temara | EPI1950349–EPI1950356 | |
A/chicken/Morocco/2_924_21RS1333-2/2019 | H9N2 | Breeders | 33 w | 3 April 2019 | Casablanca-Settat | El Jadida | EPI1950230–233, 235–237 and EPI1950385 | |
A/chicken/Morocco/3_1004_21RS1333-3/2019 | H9N2 | Broilers | 15 d | 11 April 2019 | Fès-Meknès | Fès | EPI1950246–EPI1950253 | |
A/chicken/Morocco/4_2573_21RS1333-4/2019 | H9N2 | Broilers | 34 d | 5 September 2019 | Fès-Meknès | Sefrou | EPI1950277–EPI1950284 | |
A/chicken/Morocco/5_2781_21RS1333-5/2019 | H9N2 | Breeders | 29 w | 24 September 2019 | Fès-Meknès | Fès | EPI1950254–EPI1950261 | |
A/chicken/Morocco/7_3046_21RS1333-7/2019 | H9N2 | Broilers | 24 d | 19 October 2019 | Fès-Meknès | Ifrane | EPI1950262–265, 267–269 and EPI1950386 | |
A/chicken/Morocco/8_3465_21RS1333-8/2019 | H9N2 | Layers | 22 w | 29 November 2019 | Rabat-Salé-Kénitra | Bouznika | EPI1950309–EPI1950316 | |
2020 | A/chicken/Morocco/21_4374_21RS1333-23/2020 | H9N2 | Breeders | 42 w | 25 February 2020 | Fès-Meknès | Meknès | EPI1950270–EPI1950276 |
A/chicken/Morocco/9_4410_21RS1333-9/2020 | H9N2 | Breeders | 43 w | 26 February 2020 | Fès-Meknès | EL Hajeb | EPI1950317–320, 322–324 and EPI1950387 | |
A/chicken/Morocco/10_4437_21RS1333-10/2020 | H9N2 | Broilers | 34 d | 29 February 2020 | Fès-Meknès | Sefrou | EPI1950285–EPI1950292 | |
A/chicken/Morocco/11_4438_21RS1333-11/2020 | H9N2 | Broilers | 34 d | 29 February 2020 | Fès-Meknès | Sefrou | EPI1950293–EPI1950300 | |
A/chicken/Morocco/12_4476_21RS1333-12/2020 | H9N2 | Layers | 33 w | 3 March 2020 | Rabat-Salé-Kénitra | Bouznika | EPI1950301–EPI1950308 | |
A/chicken/Morocco/13_4562_21RS1333-13/2020 | H9N2 | Broilers | 34 d | 12 March 2020 | Casablanca-Settat | El Jadida | EPI1950225–EPI1950229 | |
A/chicken/Morocco/14-1_4636_21RS1333-14/2020 | H9N2 | Broilers | 27 d | 20 March 2020 | Casablanca-Settat | Casablanca | EPI1950194–EPI1950200 | |
A/chicken/Morocco/16-1_6595_21RS1333-17/2020 | H9N2 | Broilers | 39 d | 23 November 2020 | Rabat-Salé-Kénitra | Rommani | EPI1950365–EPI1950369 | |
A/chicken/Morocco/16-2_6595_21RS1333-18/2020 | H9N2 | Broilers | 39 d | 23 November 2020 | Rabat-Salé-Kénitra | Rommani | EPI1950370–EPI1950375 | |
A/chicken/Morocco/17_6650_21RS1333-19/2020 | H9N2 | Broilers | 41 d | 28 November 2020 | Rabat-Salé-Kénitra | Bouznika | EPI1950186–EPI1950193 | |
A/chicken/Morocco/18_6657_21RS1333-20/2020 | H9N2 | Broilers | 37 d | 30 November 2020 | Rabat-Salé-Kénitra | Skhirat Temara | EPI1950357–EPI1950364 | |
A/chicken/Morocco/19_6700_21RS1333-21/2020 | H9N2 | Layers | 44 w | 4 December 2020 | Casablanca-Settat | Casablanca | EPI1950870–EPI1950875 | |
2021 | A/chicken/Morocco/22-1_72_21RS1333-24/2021 | H9N2 | Broilers | 21 d | 11 February 2021 | Tanger-Tétouan- |
Tangeir | EPI1950376–EPI1950383 |
A/chicken/Morocco/24_520_21RS1333-28/2021 | H9N2 | Broilers | 40 d | 18 March 2021 | Casablanca-Settat | Casablanca | EPI1950201–EPI1950208 | |
A/chicken/Morocco/25-1_3155_21RS1333-29/2021 | H9N2 | Broilers | 40 d | 9 April 2021 | Rabat-Salé-Kénitra | Sidi Boubker El Haj | EPI1950333–EPI1950340 | |
A/chicken/Morocco/25-2_3155_21RS1333-30/2021 | H9N2 | Broilers | 40 d | 9 April 2021 | Rabat-Salé-Kénitra | Sidi Boubker El Haj | EPI1950341–EPI1950348 | |
A/chicken/Morocco/26-2_917_21RS1333-32/2021 | H9N2 | Layers | 19 w | 28 April 2021 | Béni Mellal-Khénifra | Béni Mellal | EPI1950178–EPI1950185 |
1 d: days; w: weeks; na: not available.
Sites under positive selection (p-value < 0.05) in the Moroccan H9N2 viruses HA protein.
Method | Site |
Site |
p-Value | Aminoacids | References | |
---|---|---|---|---|---|---|
MEME | 180 | 190 | 0.03 | A, T, V | RBS, antigenic site B | [ |
216 | 226 | 0.03 | L, Q, R | RBS, antigenic site D | [ |
|
FEL | 180 | 190 | 0.02 | A, T, V | RBS, antigenic site B | [ |
216 | 226 | 0.02 | L, Q, R | RBS, antigenic site D | [ |
|
150 | 160 | 0.03 | G, N, V | antigenic site B | [ |
|
176 | 186 | 0.04 | P, A | RBS | [ |
* H9 numbering.
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Homme, P.; Easterday, B. Avian influenza virus infections. I. Characteristics of influenza A/Turkey/Wisconsin/1966 virus. Avian Dis.; 1970; pp. 66-74. [DOI: https://dx.doi.org/10.2307/1588557]
2. Shortridge, K.F. Pandemic influenza: A zoonosis?. Semin. Respir. Infect.; 1992; 7, pp. 11-25. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/1609163]
3. Sun, Y.; Liu, J. H9N2 influenza virus in China: A cause of concern. Protein Cell.; 2015; 6, pp. 18-25. [DOI: https://dx.doi.org/10.1007/s13238-014-0111-7] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25384439]
4. Butt, A.M.; Siddique, S.; Idrees, M.; Tong, Y. Avian influenza A (H9N2): Computational molecular analysis and phylogenetic characterization of viral surface proteins isolated between 1997 and 2009 from the human population. Virol. J.; 2010; 7, 319. [DOI: https://dx.doi.org/10.1186/1743-422X-7-319] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21078137]
5. Lee, C.-W.; Song, C.-S.; Lee, Y.-J.; Mo, I.-P.; Garcia, M.; Suarez, D.L.; Kim, S.-J. Sequence analysis of the hemagglutinin gene of H9N2 Korean avian influenza viruses and assessment of the pathogenic potential of isolate MS96. Avian Dis.; 2000; 44, pp. 527-535. [DOI: https://dx.doi.org/10.2307/1593091]
6. Alexander, D.J. An overview of the epidemiology of avian influenza. Vaccine; 2007; 25, pp. 5637-5644. [DOI: https://dx.doi.org/10.1016/j.vaccine.2006.10.051]
7. Nagy, A.; Mettenleiter, T.; Abdelwhab, E. 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. J. Infect. Chemother.; 2017; 145, pp. 3320-3333. [DOI: https://dx.doi.org/10.1017/S0950268817002576]
8. Yang, J.; Müller, N.F.; Bouckaert, R.; Xu, B.; Drummond, A.J. Bayesian phylodynamics of avian influenza A virus H9N2 in Asia with time-dependent predictors of migration. PLoS Comput. Biol.; 2019; 15, e1007189. [DOI: https://dx.doi.org/10.1371/journal.pcbi.1007189]
9. Peacock, T.T.P.; James, J.; Sealy, J.E.; Iqbal, M. A global perspective on H9N2 avian influenza virus. Viruses; 2019; 11, 620. [DOI: https://dx.doi.org/10.3390/v11070620]
10. Peiris, M.; Yuen, K.; Leung, C.; Chan, K.; Ip, P.; Lai, R.; Orr, W.; Shortridge, K. Human infection with influenza H9N2. Lancet; 1999; 354, pp. 916-917. [DOI: https://dx.doi.org/10.1016/S0140-6736(99)03311-5]
11. Samy, A.; Naguib, M.M. Avian respiratory coinfection and impact on avian influenza pathogenicity in domestic poultry: Field and experimental findings. Vet. Sci.; 2018; 5, 23. [DOI: https://dx.doi.org/10.3390/vetsci5010023] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29495276]
12. Naeem, K.; Ullah, A.; Manvell, R.; Alexander, D. Avian influenza A subtype H9N2 in poultry in Pakistan. Vet. Rec.; 1999; 145, 560. [DOI: https://dx.doi.org/10.1136/vr.145.19.560]
13. Umar, S.; Guerin, J.L.; Ducatez, M.F. Low pathogenic avian influenza and coinfecting pathogens: A review of experimental infections in avian models. Avian Dis.; 2017; 61, pp. 3-15. [DOI: https://dx.doi.org/10.1637/11514-101316-Review] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28301244]
14. Zhang, P.; Tang, Y.; Liu, X.; Peng, D.; Liu, W.; Liu, H.; Lu, S.; Liu, X. Characterization of H9N2 influenza viruses isolated from vaccinated flocks in an integrated broiler chicken operation in eastern China during a 5 year period (1998–2002). J. Gen. Virol.; 2008; 89, pp. 3102-3112. [DOI: https://dx.doi.org/10.1099/vir.0.2008/005652-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19008399]
15. Butt, K.; Smith, G.J.; Chen, H.; Zhang, L.; Leung, Y.C.; Xu, K.; Lim, W.; Webster, R.G.; Yuen, K.; Peiris, J.M. Human infection with an avian H9N2 influenza A virus in Hong Kong in 2003. J. Clin. Microbiol.; 2005; 43, pp. 5760-5767. [DOI: https://dx.doi.org/10.1128/JCM.43.11.5760-5767.2005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16272514]
16. Matrosovich, M.N.; Krauss, S.; Webster, R.G. H9N2 influenza A viruses from poultry in Asia have human virus-like receptor specificity. Virology; 2001; 281, pp. 156-162. [DOI: https://dx.doi.org/10.1006/viro.2000.0799] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11277689]
17. Peiris, J.; Guan, Y.; Markwell, D.; Ghose, P.; Webster, R.; Shortridge, K. Cocirculation of avian H9N2 and contemporary “human” H3N2 influenza A viruses in pigs in southeastern China: Potential for genetic reassortment?. J. Virol.; 2001; 75, pp. 9679-9686. [DOI: https://dx.doi.org/10.1128/JVI.75.20.9679-9686.2001]
18. Gu, M.; Xu, L.; Wang, X.; Liu, X. Current situation of H9N2 subtype avian influenza in China. Vet. Res.; 2017; 48, pp. 1-10. [DOI: https://dx.doi.org/10.1186/s13567-017-0453-2]
19. Gu, M.; Chen, H.; Li, Q.; Huang, J.; Zhao, M.; Gu, X.; Jiang, K.; Wang, X.; Peng, D.; Liu, X. Enzootic genotype S of H9N2 avian influenza viruses donates internal genes to emerging zoonotic influenza viruses in China. Vet. Microbiol.; 2014; 174, pp. 309-315. [DOI: https://dx.doi.org/10.1016/j.vetmic.2014.09.029]
20. Guan, Y.; Shortridge, K.; Krauss, S.; Chin, P.; Dyrting, K.; Ellis, T.; Webster, R.; Peiris, M. H9N2 influenza viruses possessing H5N1-like internal genomes continue to circulate in poultry in southeastern China. J. Virol.; 2000; 74, pp. 9372-9380. [DOI: https://dx.doi.org/10.1128/JVI.74.20.9372-9380.2000]
21. Shen, Y.-Y.; Ke, C.-W.; Li, Q.; Yuan, R.-Y.; Xiang, D.; Jia, W.-X.; Yu, Y.-D.; Liu, L.; Huang, C.; Qi, W.-B. Novel reassortant avian influenza A (H5N6) viruses in humans, Guangdong, China, 2015. Emerg. Infect. Dis.; 2016; 22, 1507. [DOI: https://dx.doi.org/10.3201/eid2208.160146] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27331418]
22. Zhang, Z.; Li, R.; Jiang, L.; Xiong, C.; Chen, Y.; Zhao, G.; Jiang, Q. The complexity of human infected AIV H5N6 isolated from China. BMC Infect. Dis.; 2016; 16, 600. [DOI: https://dx.doi.org/10.1186/s12879-016-1932-1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27782815]
23. RahimiRad, S.; Alizadeh, A.; Alizadeh, E.; Hosseini, S.M. The avian influenza H9N2 at avian-human interface: A possible risk for the future pandemics. J. Res. Med. Sci.; 2016; 21, 51. [DOI: https://dx.doi.org/10.4103/1735-1995.187253] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28083072]
24. Cui, L.; Liu, D.; Shi, W.; Pan, J.; Qi, X.; Li, X.; Guo, X.; Zhou, M.; Li, W.; Li, J. Dynamic reassortments and genetic heterogeneity of the human-infecting influenza A (H7N9) virus. Nat. Commun.; 2014; 5, 3142. [DOI: https://dx.doi.org/10.1038/ncomms4142]
25. Di Liu, W.S.; Gao, G.F. Poultry carrying H9N2 act as incubators for novel human avian influenza viruses. Lancet; 2014; 384, 869. [DOI: https://dx.doi.org/10.1016/S0140-6736(14)60386-X]
26. Shi, W.; Li, W.; Li, X.; Haywood, J.; Ma, J.; Gao, G.F.; Liu, D. Phylogenetics of varied subtypes of avian influenza viruses in China: Potential threat to humans. Protein Cell.; 2014; 5, pp. 253-257. [DOI: https://dx.doi.org/10.1007/s13238-014-0036-1]
27. Wan, H.; Sorrell, E.M.; Song, H.; Hossain, M.J.; Ramirez-Nieto, G.; Monne, I.; Stevens, J.; Cattoli, G.; Capua, I.; Chen, L.-M. Replication and transmission of H9N2 influenza viruses in ferrets: Evaluation of pandemic potential. PLoS ONE; 2008; 3, e2923. [DOI: https://dx.doi.org/10.1371/journal.pone.0002923]
28. Choi, Y.; Ozaki, H.; Webby, R.; Webster, R.; Peiris, J.; Poon, L.; Butt, C.; Leung, Y.; Guan, Y. Continuing evolution of H9N2 influenza viruses in Southeastern China. J. Virol.; 2004; 78, pp. 8609-8614. [DOI: https://dx.doi.org/10.1128/JVI.78.16.8609-8614.2004]
29. Yong-Feng, Z.; Fei-Fei, D.; Jia-Yu, Y.; Feng-Xia, Z.; Chang-Qing, J.; Jian-Li, W.; Shou-Yu, G.; Kai, C.; Chuan-Yi, L.; Xue-Hua, W. Intraspecies and interspecies transmission of mink H9N2 influenza virus. Sci. Rep.; 2017; 7, 7429. [DOI: https://dx.doi.org/10.1038/s41598-017-07879-1]
30. Zhang, K.; Zhang, Z.; Yu, Z.; Li, L.; Cheng, K.; Wang, T.; Huang, G.; Yang, S.; Zhao, Y.; Feng, N. Domestic cats and dogs are susceptible to H9N2 avian influenza virus. Virus Res.; 2013; 175, pp. 52-57. [DOI: https://dx.doi.org/10.1016/j.virusres.2013.04.004]
31. Li, C.; Wang, S.; Bing, G.; Carter, R.A.; Wang, Z.; Wang, J.; Wang, C.; Wang, L.; Wu, G.; Webster, R.G. Genetic evolution of influenza H9N2 viruses isolated from various hosts in China from 1994 to 2013. Emerg. Microbes Infect.; 2017; 6, e106. [DOI: https://dx.doi.org/10.1038/emi.2017.94] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29184157]
32. Awuni, J.A.; Bianco, A.; Dogbey, O.J.; Fusaro, A.; Yingar, D.T.; Salviato, A.; Ababio, P.T.; Milani, A.; Bonfante, F.; Monne, I. Avian influenza H9N2 subtype in Ghana: Virus characterization and evidence of co-infection. Avian Pathol.; 2019; 48, pp. 470-476. [DOI: https://dx.doi.org/10.1080/03079457.2019.1624687]
33. 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, pp. 422-430. [DOI: https://dx.doi.org/10.1637/11068-032215-ResNote.1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26478162]
34. El-Zoghby, E.F.; Arafa, A.-S.; Hassan, M.K.; Aly, M.M.; Selim, A.; Kilany, W.H.; Selim, U.; Nasef, S.; Aggor, M.G.; Abdelwhab, E. Isolation of H9N2 avian influenza virus from bobwhite quail (Colinus virginianus) in Egypt. Arch. Virol.; 2012; 157, pp. 1167-1172. [DOI: https://dx.doi.org/10.1007/s00705-012-1269-z] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22426861]
35. Tombari, W.; Nsiri, J.; Larbi, I.; Guerin, J.L.; Ghram, A. Genetic evolution of low pathogenecity H9N2 Avian influenza viruses in Tunisia: Acquisition of new mutations. Virol. J.; 2011; 8, 467. [DOI: https://dx.doi.org/10.1186/1743-422X-8-467]
36. 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. [DOI: https://dx.doi.org/10.1186/s12985-016-0596-1]
37. Jeevan, T.; Darnell, D.; Gradi, E.A.; Benali, Y.; Kara, R.; Guetarni, D.; Rubrum, A.; Seiler, P.J.; Crumpton, J.C.; Webby, R.J. A (H9N2) influenza viruses associated with chicken mortality in outbreaks in Algeria 2017. Influenza Other Respir. Viruses; 2019; 13, pp. 622-626. [DOI: https://dx.doi.org/10.1111/irv.12675]
38. Zecchin, B.; Minoungou, G.; Fusaro, A.; Moctar, S.; Ouedraogo-Kaboré, A.; Schivo, A.; Salviato, A.; Marciano, S.; Monne, I. Influenza A (H9N2) Virus, Burkina Faso. Emerg. Infect. Dis.; 2017; 23, 2118. [DOI: https://dx.doi.org/10.3201/eid2312.171294]
39. 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. [DOI: https://dx.doi.org/10.1016/j.meegid.2019.104074]
40. 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. 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, pp. 753-761. [DOI: https://dx.doi.org/10.1080/22221751.2021.1908097]
41. Boumart, Z.; Bamouh, Z.; Jazouli, M.; Zecchin, B.; Fusaro, A.; Salviato, A.; Monne, I.; Tadlaoui, K.O.; Harrak, M.E. Pathogenicity and full genome sequencing of the avian influenza H9N2 Moroccan isolate 2016. Avian Dis.; 2019; 63, pp. 24-30. [DOI: https://dx.doi.org/10.1637/11941-080418-Reg.1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31251516]
42. Belkasmi, S.F.; Fellahi, S.; Touzani, C.D.; Faraji, F.Z.; Maaroufi, I.; Delverdier, M.; Guérin, J.-L.; Fihri, O.F.; El Houadfi, M.; Ducatez, M.F. Co-infections of chickens with avian influenza virus H9N2 and Moroccan Italy 02 infectious bronchitis virus: Effect on pathogenesis and protection conferred by different vaccination programmes. Avian Pathol.; 2020; 49, pp. 21-28. [DOI: https://dx.doi.org/10.1080/03079457.2019.1656328] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31412705]
43. El Khantour, A.; Soulaymani, A.; Salek, M.; Maltouf, A.F.; El Mellouli, F.; Ducatez, M.F.; Fellahi, S. Molecular characterization of the hemagglutinin gene of H9N2 avian influenza viruses isolated from broiler flocks in Morocco from 2016 to 2018. Vet. Arh.; 2020; 90, pp. 477-484. [DOI: https://dx.doi.org/10.24099/vet.arhiv.0724]
44. Spackman, E.; Senne, D.A.; Myers, T.J.; Bulaga, L.L.; Garber, L.P.; Perdue, M.L.; Lohman, K.; Daum, L.T.; Suarez, D.L. Development of a Real-Time Reverse Transcriptase PCR Assay for Type A Influenza Virus and the Avian H5 and H7 Hemagglutinin Subtypes. J. Clin. Microbiol.; 2002; 40, pp. 3256-3260. [DOI: https://dx.doi.org/10.1128/JCM.40.9.3256-3260.2002]
45. 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, pp. 1769-1773. [DOI: https://dx.doi.org/10.1128/JCM.02204-07]
46. Zhou, B.; Donnelly, M.E.; Scholes, D.T.; George, K.S.; Hatta, M.; Kawaoka, Y.; Wentworth, D.E. Single-Reaction Genomic Amplification Accelerates Sequencing and Vaccine Production for Classical and Swine Origin Human Influenza A Viruses. J. Virol.; 2009; 83, pp. 10309-10313. [DOI: https://dx.doi.org/10.1128/JVI.01109-09]
47. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics; 2014; 30, pp. 2114-2120. [DOI: https://dx.doi.org/10.1093/bioinformatics/btu170]
48. Li, H.; Durbin, R. Fast and accurate long-read alignment with Burrows–Wheeler transform. Bioinformatics; 2010; 26, pp. 589-595. [DOI: https://dx.doi.org/10.1093/bioinformatics/btp698]
49. DePristo, M.A.; Banks, E.; Poplin, R.; Garimella, K.V.; Maguire, J.R.; Hartl, C.; Philippakis, A.A.; Del Angel, G.; Rivas, M.A.; Hanna, M. et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet.; 2011; 43, pp. 491-498. [DOI: https://dx.doi.org/10.1038/ng.806]
50. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res.; 2010; 20, pp. 1297-1303. [DOI: https://dx.doi.org/10.1101/gr.107524.110]
51. Wilm, A.; Aw, P.P.K.; Bertrand, D.; Yeo, G.H.T.; Ong, S.H.; Wong, C.H.; Khor, C.C.; Petric, R.; Hibberd, M.L.; Nagarajan, N. LoFreq: A sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets. Nucleic Acids Res.; 2012; 40, pp. 11189-11201. [DOI: https://dx.doi.org/10.1093/nar/gks918] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23066108]
52. Katoh, K.; Standley, D.M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol.; 2013; 30, pp. 772-780. [DOI: https://dx.doi.org/10.1093/molbev/mst010] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23329690]
53. Nguyen, L.-T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol.; 2014; 32, pp. 268-274. [DOI: https://dx.doi.org/10.1093/molbev/msu300] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25371430]
54. Hoang, D.T.; Chernomor, O.; Von Haeseler, A.; Minh, B.Q.; Vinh, L.S. UFBoot2: Improving the Ultrafast Bootstrap Approximation. Mol. Biol. Evol.; 2017; 35, pp. 518-522. [DOI: https://dx.doi.org/10.1093/molbev/msx281] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29077904]
55. Bandelt, H.J.; Forster, P.; Röhl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol.; 1999; 16, pp. 37-48. [DOI: https://dx.doi.org/10.1093/oxfordjournals.molbev.a026036]
56. 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. [DOI: https://dx.doi.org/10.1093/ve/vey016]
57. Shapiro, B.; Rambaut, A.; Drummond, A.J. Choosing Appropriate Substitution Models for the Phylogenetic Analysis of Protein-Coding Sequences. Mol. Biol. Evol.; 2005; 23, pp. 7-9. [DOI: https://dx.doi.org/10.1093/molbev/msj021]
58. Kosakovsky Pond, S.L.; Frost, S.D.W. Not So Different After All: A Comparison of Methods for Detecting Amino Acid Sites Under Selection. Mol. Biol. Evol.; 2005; 22, pp. 1208-1222. [DOI: https://dx.doi.org/10.1093/molbev/msi105]
59. Murrell, B.; Wertheim, J.O.; Moola, S.; Weighill, T.; Scheffler, K.; Kosakovsky Pond, S.L. Detecting Individual Sites Subject to Episodic Diversifying Selection. PLoS Genet.; 2012; 8, e1002764. [DOI: https://dx.doi.org/10.1371/journal.pgen.1002764]
60. Delport, W.; Poon, A.F.Y.; Frost, S.D.W.; Kosakovsky Pond, S.L. Datamonkey 2010: A suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics; 2010; 26, pp. 2455-2457. [DOI: https://dx.doi.org/10.1093/bioinformatics/btq429]
61. 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, pp. 5181-5191. [DOI: https://dx.doi.org/10.1128/JVI.02827-06] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17344280]
62. 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, e11. [DOI: https://dx.doi.org/10.1038/emi.2016.139] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28325922]
63. Peacock, T.P.; Harvey, W.T.; Sadeyen, J.-R.; Reeve, R.; Iqbal, M. The molecular basis of antigenic variation among A(H9N2) avian influenza viruses. Emerg. Microbes Infect.; 2018; 7, 176. [DOI: https://dx.doi.org/10.1038/s41426-018-0178-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30401826]
64. Suptawiwat, O.; Jeamtua, W.; Boonarkart, C.; Kongchanagul, A.; Puthawathana, P.; Auewarakul, P. Effects of the Q223R mutation in the hemagglutinin (HA) of egg-adapted pandemic 2009 (H1N1) influenza A virus on virus growth and binding of HA to human- and avian-type cell receptors. Acta Virol.; 2013; 57, pp. 333-338. Available online: http://europepmc.org/abstract/MED/24020758 (accessed on 4 January 2022).
65. He, L.; Jiang, K.; Wu, Q.; Duan, Z.; Xu, H.; Liu, J.; Cui, Z.; Gu, M.; Wang, X.; Liu, X. Two amino acid substitutions in the haemagglutinin of the 2009 pandemic H1N1 virus decrease direct-contact transmission in guinea pigs. J. Gen. Virol.; 2014; 95, pp. 2612-2617. [DOI: https://dx.doi.org/10.1099/vir.0.067694-0]
66. Peacock, T.; Reddy, K.; James, J.; Adamiak, B.; Barclay, W.; Shelton, H.; Iqbal, M. Antigenic mapping of an H9N2 avian influenza virus reveals two discrete antigenic sites and a novel mechanism of immune escape. Sci. Rep.; 2016; 6, 18745. [DOI: https://dx.doi.org/10.1038/srep18745]
67. Kaverin, N.V.; Rudneva, I.A.; Ilyushina, N.A.; Lipatov, A.S.; Krauss, S.; Webster, R.G. Structural Differences among Hemagglutinins of Influenza A Virus Subtypes Are Reflected in Their Antigenic Architecture: Analysis of H9 Escape Mutants. J. Virol.; 2004; 78, pp. 240-249. [DOI: https://dx.doi.org/10.1128/JVI.78.1.240-249.2004]
68. Wan, Z.; Ye, J.; Xu, L.; Shao, H.; Jin, W.; Qian, K.; Wan, H.; Qin, A.; Lyles, D.S. Antigenic Mapping of the Hemagglutinin of an H9N2 Avian Influenza Virus Reveals Novel Critical Amino Acid Positions in Antigenic Sites. J. Virol.; 2014; 88, pp. 3898-3901. [DOI: https://dx.doi.org/10.1128/JVI.03440-13]
69. El Khantour, A.; El Houadfi, M.; Nassik, S.; Tligui, N.S.; El Mellouli, F.; Sikht, F.-Z.; Ducatez, M.F.; 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, pp. 351-357. [DOI: https://dx.doi.org/10.1637/aviandiseases-D-21-00015]
70. Kaboudi, K. Low pathogenic avian influenza virus subtype H9N2 in poultry in north Africa: Current status. Vet. Sci. Res. Rev.; 2019; 5, pp. 73-79. [DOI: https://dx.doi.org/10.17582/journal.vsrr/2019/5.2.73.79]
71. Amal, E.-B.; Saâdi, N.; Asma, F.; Moncef, B.; Ouafae, F.F. Characterization and Phylogenetic Analysis of the Hemagglutinin Gene in H9 Influenza Viruses from Chickens in Morocco from 2017 to 2019. Avian Dis.; 2020; 64, pp. 310-314. [DOI: https://dx.doi.org/10.1637/aviandiseases-D-20-00009] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33205168]
72. Hassan, N.; Kaoutar, M.; Sara, G.; Mustapha, M. Focus on Poultry in Morocco. Adv. Environ. Biol.; 2017; 11, 209678240.Available online: https://link.gale.com/apps/doc/A518588590/AONE?u=anon~de5f3275&sid=googleScholar&xid=6178e5a7 (accessed on 6 January 2022).
73. 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, pp. 1527-1540. [DOI: https://dx.doi.org/10.1007/s00705-020-04624-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32335769]
74. Jegede, A.; Fu, Q.; Berhane, Y.; Lin, M.; Kumar, A.; Guan, J. H9N2 avian influenza virus retained low pathogenicity after serial passage in chickens. Can. J. Vet. Res.; 2018; 82, pp. 131-138. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29755193]
75. Jin, H.; Wang, W.; Yang, X.; Su, H.; Fan, J.; Zhu, R.; Wang, S.; Shi, H.; Liu, X. Evolution of H9N2 avian influenza virus in embryonated chicken eggs with or without homologous vaccine antibodies. BMC Vet. Res.; 2018; 14, 71. [DOI: https://dx.doi.org/10.1186/s12917-018-1391-6]
76. Park, K.J.; Kwon, H.-I.; 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. Rapid evolution of low-pathogenic H9N2 avian influenza viruses following poultry vaccination programmes. J. Gen. Virol.; 2011; 92, pp. 36-50. [DOI: https://dx.doi.org/10.1099/vir.0.024992-0]
77. El Mellouli, F.; Abouchoaib, N.; Zekhnini, H.; Khayli, M.; Fusaro, A.; Idrissi, H.R.; Benhoussa, A. Molecular Detection of Avian Influenza Virus in Wild Birds in Morocco, 2016–2019. Avian Dis.; 2021; 66, pp. 1-10. [DOI: https://dx.doi.org/10.1637/aviandiseases-D-21-00070]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 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/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
The H9N2 virus continues to spread in wild birds and poultry worldwide. At the beginning of 2016, the H9N2 Avian influenza virus (AIV) was detected in Morocco for the first time; despite the implementation of vaccination strategies to control the disease, the virus has become endemic in poultry in the country. The present study was carried out to investigate the origins, zoonotic potential, as well as the impact of vaccination on the molecular evolution of Moroccan H9N2 viruses. Twenty-eight (28) H9N2 viruses collected from 2016 to 2021 in Moroccan poultry flocks were isolated and their whole genomes sequenced. Phylogenetic and evolutionary analyses showed that Moroccan H9N2 viruses belong to the G1-like lineage and are closely related to viruses isolated in Africa and the Middle East. A high similarity among all the 2016–2017 hemagglutinin sequences was observed, while the viruses identified in 2018–2019 and 2020–2021 were separated from their 2016–2017 ancestors by long branches. Mutations in the HA protein associated with antigenic drift and increased zoonotic potential were also found. The Bayesian phylogeographic analyses revealed the Middle East as being the region where the Moroccan H9N2 virus may have originated, before spreading to the other African countries. Our study is the first comprehensive analysis of the evolutionary history of the H9N2 viruses in the country, highlighting their zoonotic potential and pointing out the importance of implementing effective monitoring systems.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details





1 Biodiversity, Ecology and Genome Laboratory, Faculty of Sciences, Mohammed V University in Rabat, 4 Avenue Ibn Battouta, Rabat 10106, Morocco;
2 Mouahid’s Veterinary Clinic, Temara 12000, Morocco;
3 Istituto Zooprofilattico Sperimentale delle Venezie, 35020 Legnaro, Italy;
4 Immunology and Biodiversity Laboratory, Faculty of Science Ain chock, Hassan II University of Casablanca, Casablanca 20100, Morocco;
5 Division of Pharmacy and Veterinary Inputs, Food Safety National Office (ONSSA), Rue Ikhlas Cym, PB 4509 Akkari, Rabat 10050, Morocco;
6 Biodiversity, Ecology and Genome Laboratory, Faculty of Sciences, Mohammed V University in Rabat, 4 Avenue Ibn Battouta, Rabat 10106, Morocco;