Skip to main content
  • Research article
  • Open access
  • Published:

Predominance of low pathogenic avian influenza virus H9N2 in the respiratory co-infections in broilers in Tunisia: a longitudinal field study, 2018–2020


Respiratory diseases are a health and economic concern for poultry production worldwide. Given global economic exchanges and migratory bird flyways, respiratory viruses are likely to emerge continuously in new territories. The primary aim of this study was to investigate the major pathogens involved in respiratory disease in Tunisian broiler poultry and their epidemiology. Between 2018 and 2020, broilers farms in northeastern Tunisia were monitored, and 39 clinically diseased flocks were sampled. Samples were screened for five viral and three bacterial respiratory pathogens using a panel of real-time PCR assays. The reemergence of H9N2 low pathogenic avian influenza virus (LPAIV) in commercial poultry was reported, and the Northern and Western African GI lineage strain was typed. The infectious bronchitis virus (IBV) GI-23 lineage and the avian metapneumovirus (aMPV) subtype B also were detected for the first time in broilers in Tunisia. H9N2 LPAIV was the most detected pathogen in the flocks tested, but rarely alone, as 15 of the 16 H9N2 positive flocks were co-infected. Except for infectious laryngotracheitis virus (ILTV), all of the targeted pathogens were detected, and in 61% of the respiratory disease cases, a combination of pathogens was identified. The major combinations were H9N2 + aMPV (8/39) and H9N2 + IBV (6/39), showing the high contribution of H9N2 LPAIV to the multifactorial respiratory diseases. This field survey provided evidence of the emergence of new respiratory viruses and the complexity of respiratory disease in Tunisia. A comprehensive and continuous surveillance strategy therefore is needed to better control respiratory pathogens in Tunisia.


Over the past 20 years, the commercial poultry sector has grown rapidly in Tunisia. Production levels have doubled, and poultry have become the main source of protein in the country (FAO). Respiratory diseases pose a major health threat to commercial poultry worldwide, causing tremendous economic losses. The H9N2 low pathogenic avian influenza virus (LPAIV) is the most prevalent avian influenza virus in the world [1]. In the late 1990s, H9N2 LPAIV already had been detected in Southeast Asia and the Middle East in domestic poultry. Since the 2000s, H9N2 LPAIV has become enzootic in Asia, the Middle East, and North and West Africa [2].

H9N2 LPAIV was detected in Tunisia for the first time in 2009. The strain, assumed to have originated from Pakistan, was replaced in 2012 by a genetically related virus that originated from the United Arab Emirates (UAE) and spread to Tunisia through Libya [3,4,5]. In early 2018, the surveillance of H9N2 LPAIV in migratory birds in Tunisia enabled the detection of two different strains of H9N2 LPAIV from the same lagoon in the northeast region of the country. The first was closely related to the Tunisian H9N2 strain of 2012 and was isolated from wild birds, whereas the second fell in the Northern and Western African H9N2 cluster and was isolated from lagoon water [6]. The earliest virus in the latter cluster was detected in Morocco in 2016, and has since become endemic in several North and West African countries like Algeria, Burkina-Faso, Ghana, Togo and Benin [5, 7,8,9,10]. An H9N2 virus from the same cluster was detected in an infant in Senegal in 2019, exhibiting its potential zoonotic risk [11]. Over the past decade, other respiratory viruses affecting poultry also have emerged in Tunisia, including Newcastle disease virus (NDV) genotype VII.2 in 2013, a novel variant of infectious bronchitis virus (IBV) in 2016, and infectious laryngotracheitis virus (ILTV) between 2013 and 2016 [12, 13].

Previous studies have provided evidence of the multifactorial origin of respiratory diseases. A primary infection can be complicated by environmental factors (ammonia, dust, moisture etc.) and/or by co-infection with other pathogens [14,15,16]. Co-infection is likely to enhance the clinical signs and mortality of each of the pathogens involved through synergistic mechanisms operating between them. This seems to be true for H9N2 LPAIV, which under experimental conditions barely induces clinical signs in chickens, while in the field is correlated with severe clinical signs and high mortality [10, 17]. Several experimental infections succeeded in reproducing H9N2 LPAIV clinical signs in the field by co-infecting H9N2 alongside another pathogen such as IBV, NDV or E. coli [18,19,20,21].

Given the repeated emergence and circulation of respiratory pathogens in the Tunisian commercial poultry sector, a comprehensive and longitudinal surveillance study is needed. In this study, we monitored broiler flocks in northeast Tunisia over 3 years to investigate the pathogens involved in the respiratory disease outbreaks observed and assess their longitudinal circulation. In the end we brought to light the emergence of three respiratory viruses on broiler farms in Tunisia between 2018 and 2020. We also conveyed that the respiratory diseases have a multifactorial etiology as 2/3rd of the flocks tested were co-infected.

Materials and methods

Context and sampling protocol

From January 2018 to June 2020, six farms belonging to the same broiler rearing company in Tunisia were monitored for respiratory outbreaks. These farms were located either in Ben Arous or in Nabeul, two northeastern Tunisian governorates (Figure 1). The farms were named using letters from A to F. Flocks showing acute respiratory signs were sampled within the first 5 days of an outbreak. The number of flocks sampled from the same farm during an outbreak depended on the health status of each one of the flocks. For each flock, 12 tracheal swabs were collected from diseased or freshly dead birds and immediately smeared onto the four circles of an FTA® card, with three swabs per circle. FTA® cards were dried for at least five minutes, labelled, placed in a Ziploc bag and transported on ice to be stored at −20 °C until later use [22]. A written form including information about the GPS location, age, date of the first respiratory signs, daily mortality and vaccination history was provided with each sampled flock.

Figure 1
figure 1

Geographic distribution of tested flocks on monitored and external broiler farms in Tunisia. Circles on the map represents the surveyed farms and the numerical value represents the cumulative number of sampled flocks in each farm. The monitored farms were distributed in two governorates, Nabeul in the east of the country, and Ben Arous in the west. The map was created using QGIS3 software

FTA card preparation and nucleic acid extraction

Three disks of 2 mm diameter were punched out of each FTA® card’s circle using a Whatman Harris Uni-Core puncher and pooled together in a 2 mL microcentrifuge tube with 300 µL of TE Buffer (10 mM Tris, 1 mM EDTA, pH8.0) to soak the disks. The puncher was cleaned using pure ethanol between circles and cards. For nucleic acid elution, microcentrifuge tubes were vortexed for two hours with Vortex-Genie® 2 then briefly centrifuged. Then 150 µL of eluent were used for RNA and DNA extraction using the NucleoSpin RNA Virus extraction kit (Macherey–Nagel, Düren, Germany) following the manufacturer’s instructions.

Screening of respiratory pathogens

Samples were screened for a panel of eight respiratory pathogens including H9N2 LPAIV, IBV, NDV, aMPV and ILTV, Mycoplasma gallisepticum (MG), Mycoplasma synoviae (MS), and Ornithobacterium rhinotracheale (ORT). Primers used for this study are listed in Additional file 1 [23,24,25,26,27,28,29,30]. Depending on the different pathogens’ nucleic acid class, real-time PCR (qPCR) or RT-real-time PCR (RT-qPCR) assays were performed on an Agilent Bravo Automated Liquid Handling Platform (Santa Clara, USA) and an Applied Biosystems ViiA7 thermocycler (Foster City, USA).

RT-qPCR for RNA amplification

The RNA viruses were screened using the iTaq universal SYBR® green one-step RT-qPCR kit (Bio-Rad, Hercules, USA). The reaction mix consisted of 5 µL of 2 × iTaq mix, 0.125 µL of iScript reverse transcriptase, 0.3 µL of each primer (10 µM), 2.275 µL of nuclease free water and 2 µL of RNA with the following program: reverse transcription at 50 °C for 10 min and polymerase activation at 95 °C for 1 min followed by 35 cycles of denaturation at 95 °C for 15 s and annealing/extension at 60 °C for 60 s. AIV M gene positive samples were further subtyped using [31] H9 primers following the same reaction setup and thermal cycling conditions of one-step RT-qPCR on a LightCycler 96® system (Roche Applied Science, Penzberg, Germany).

qPCR for DNA amplification

The DNA viruses and bacteria were screened using LightCycler® 480 SYBR Green I Master qPCR kit (Roche Applied Science, Mannheim, Germany). The reaction mix consisted of 5 µL of Master mix, 0.2 µL of each primer (10 µM), 2.6 µL of nuclease free water and 2 µL of DNA with the following program: pre-incubation at 95 °C for 10 min followed by 45 cycles of denaturation at 95 °C for 10 s, annealing at 60 °C for 15 s and extension at 72 °C for 15 s.

Sanger sequencing

Sequencing was performed only on material from FTA cards. The HA gene of H9N2 LPAIV was sequenced using universal primers described by Hoffmann et al. [32]. The complete fusion gene of NDV was sequenced using Esmaelizad et al. [33] primers to differentiate vaccine type strains from field challenge strains. A partial sequence (~950 base pairs) overlapping the small hydrophobic (SH) and attachment (G) gene of aMPV was generated using SH1 + and G6- primers described by Cecchinato et al. [34]. The partial S1 gene sequence of IBV (~300 base pairs) spanning the hypervariable region 3 (HVR3) was generated using Worthington et al. [35] nested PCR. For all viruses, PCR products were loaded on agarose gel and DNA fragments with the expected length were excised and purified using NucleoSpin® Gel and PCR Clean-up kit (Macherey–Nagel, Düren, Germany) according to the supplier’s instructions. Purified DNA was prepared for Sanger sequencing and shipped to the sequencing external service of Eurofins genomics (Cologne, Germany).

Sequence treatment and phylogenetic analysis

Fasta format sequences were manually treated using the BioEdit v7.2.5 software package and consensus sequences were generated. The BLAST program [36] was used for comparison with sequences in the NCBI database. Sequences with high identity percentages were used to construct phylogenetic trees. MAFFT multiple sequence alignment program online version 7 was used to align sequences [37]. Maximum likelihood phylogenetic trees with 1000 bootstraps were created using MEGA 10 software: Molecular Evolutionary Genetics Analysis version 10.


Farms and flocks

During the study period, a total of 39 flocks were investigated and sampled between January 2018 and June 2020. Due to logistical reasons, monitoring and sampling activities were interrupted between April and October 2019. Of these 39 flocks, 34 originated from the six monitored farms in the Ben Arous and Nabeul governorates (Figure 1). The number of flocks sampled from farms A, B, C, D, E and F were 1, 15, 2, 10, 4 and 2, respectively (Additional file 2). The five additional flocks came from farms that were not affiliated with the rearing company (referred to as “external farms”) (Additional file 2) and were located in the Ben Arous, Nabeul, Bizerte and Manouba governorates (Figure 1). The median age of all 39 flocks studied was 31 days, with a minimum age of 15 days. The daily mortality of the flocks studied ranged from 0.05 to 14.36% during the outbreak (Additional file 3). The respiratory signs observed included sneezing, coughing, foamy eye, and dyspnea. In terms of lesions, tracheitis, fibrinous cast in the trachea, pneumonia, petechial proventriculitis and inflammation of Harderian gland, spleen and caeca tonsils were recorded. All flocks were vaccinated against IBV and NDV. The IBV vaccination schedule comprised priming at the hatchery with a Mass-like vaccine and boosting on the farm with a 793B-like vaccine. For NDV, birds were primed at the hatchery with either a combination of a recombinant and a live virus vaccine or with a live virus vaccine only, and were boosted every 10 days using a live virus vaccine. None of the flocks had been vaccinated for H9N2 LPAIV or aMPV.

Molecular epidemiology

The qPCR assays showed a high detection of viral pathogens with at least one of the five screened viruses per flock. For the 39 tested flocks, 56 positive results for viruses were recorded. H9N2 LPAIV, IBV, aMPV and NDV were detected in 21 (54%), 20 (53%), 11 (28%), and 4 (10%) flocks, respectively. ILTV was not detected in any of the flocks. MS, ORT and MG were detected in 13 (33%), 12 (31%) and one (3%) flock, respectively (Table 1).

Table 1 qPCR positiveresults and number of pathogens detected per flock

Although monitoring began in January 2018, it was not until April 2018 that H9N2 LPAIV was first detected. This virus continued to be detected until the end of the study on all of the monitored farms, exhibiting two waves of outbreaks. IBV was detected throughout the 3 year study on all of the monitored farms except farms A and C. NDV was recorded in two flocks of farm D and in one external flock during the same period in February 2019. aMPV was first detected in February 2020, and continued to be detected until the end of the study on monitored and external farms (Table 1).

Respiratory co-infections

Altogether, the PCR tests performed on the 39 study flocks showed a total of 82 positive results, covering 4 viruses and 3 bacteria included in the panel, which confirms frequent co-infections and the diversity of respiratory pathogens in these Tunisian poultry holdings. The co-infections were detected in 24 (61%) of the 39 flocks and involved viruses and/or bacteria in combinations of two, three, four or five pathogens in 14, 2, 7 and 1 flock, respectively (Table 1). Co-infections involving only viruses occurred in 16 flocks. These viral-viral co-infection combinations were mainly H9N2 + aMPV, H9N2 + IBV, IBV + aMPV and H9N2 + IBV + aMPV, profiled in 8, 6, 1 and 1 co-infections, respectively. H9N2 LPAIV was the most involved pathogen, found in 15 of the 16 viral co-infected flocks. Co-infections involving viruses and bacteria occurred in 18 of the 24 co-infected flocks, and were divided in two groups. The first consisted of a bacterial superinfection of a viral-viral co-infection and occurred in 10 flocks; the second consisted of a bacterial superinfection of a single virus infection and occurred in 8 flocks. Again, H9N2 LPAIV was the most involved pathogen in the co-infections, with 18/21 H9N2 LPAIV positive flocks. H9N2 LPAIV co-infections were detected with ORT, MS, aMPV and IBV in 11, 10, 9, 7 flocks, respectively.

Phylogenetic analysis


Complete (1.7 kb) or partial (1 kb) HA sequences of H9N2 LPAIVs from seven flocks were generated. HA gene sequences and related metadata were deposited in the GenBank database under accession numbers OQ179924-OQ179930 (Additional file 4). The phylogenetic tree showed that Tunisian H9N2 LPAIV viruses belonged to the h9.4.2 clade of the G1-like lineage and have a common ancestor with the United Arab Emirates H9N2 LPAIV virus (JX273562.1) (Figure 2). Our Tunisian H9N2 viruses clustered with those detected in Algeria and Morocco. The H9N2 LPAIV viruses detected during this study exhibit a nucleotide pairwise distance of 10% from the Tunisian H9N2 LPAIV of 2010.

Figure 2
figure 2

Phylogenetic tree of H9N2 LPAIV. The tree was constructed using the HA gene nucleotide sequences. Sequences from this study were labelled with red triangle, Tunisian sequences available in GenBank were labelled with black circle and G1 lineage sequences from Africa and the Middle East and reference sequences of other H9N2 lineages were unlabelled. The evolutionary history was inferred using the Maximum Likelihood method with General Time Reversible model and 1000 bootstrap replications in MEGA X. All nucleotide positions containing gaps and missing data were eliminated. A total of 924 positions were included in the final dataset. Only bootstrap values higher than 60 were conserved

Other viruses


Complete F gene sequences of NDVs from one monitored flock and one external flock were generated. F gene sequences and related metadata were deposited in the GenBank database under accession numbers OQ199524 and OQ199525 (Additional file 4). Referring to the revised nomenclature recommended by [38], the phylogenetic analysis showed that NDV sequences belonged to the genotype VII.2 (formerly VII-i) and shared 97% of nucleotide identity with the 2013 and 2015 Tunisian genotype VII.2 strains (Figure 3).

Figure 3
figure 3

Phylogenetic tree of NDV. The tree was constructed using F gene nucleotide sequences of Class II NDV. Sequences from this study were labelled with filled red triangle, previous NDV Tunisian sequences available in GenBank were labelled with black circle, 2019 NDV sequences from Tunisian backyard chickens were labelled with empty red triangle and reference sequences were unlabelled. The evolutionary history was inferred using the Maximum Likelihood method with Tamura 3-parameter model and 1000 bootstrap replications in MEGA X. All nucleotide positions containing gaps and missing data were eliminated. A total of 1675 positions were included in the final dataset. Only bootstrap values higher than 60 were conserved


Partial SH and G gene sequences of aMPVs were generated from eight flocks. G gene sequences and related metadata were deposited in the GenBank database under accession numbers OQ199508 and OQ199515 (Additional file 4). The phylogenetic analysis sorted G gene sequences into two different clusters within the subtype B, and were distinct from aMPV vaccine strains (Figure 4). The alignment of the SH and G gene sequences of the detected aMPVs and the reference sequence of aMPV revealed some genetic polymorphic features. These features consist of two or four nucleotide insertions in the noncoding region between the SH and G genes. Based on these two features, we were able to retrieve the two clusters obtained by the phylogenetic analysis of the G gene sequences. Two motifs of nucleotide insertions were found; the first involved the insertion of two nucleotides and the second the insertion of four nucleotides (Figure 5).

Figure 4
figure 4

Phylogenetic tree of aMPV. The tree was constructed using G gene nucleotide sequences. Sequences from the present study were labelled with red triangle, previous Tunisian sequence were labelled with black circle, vaccine sequences were labelled with green circle and other aMPV-B sequences available in GenBank were unlabelled. The evolutionary history was inferred using the Maximum Likelihood method with Tamura 3-parameter model and 1000 bootstrap replications in MEGA X. All nucleotide positions containing gaps and missing data were eliminated. A total of 292 positions were included in the final dataset. Only bootstrap values higher than 65 were conserved

Figure 5
figure 5

The genetic polymorphic features of Tunisian aMPV sequences. The alignment of the noncoding region between SH and G genes of Tunisian aMPV sequences from this study with VCO3 reference sequence showed the presence of two motifs of nucleotide insertions highlighted in the red frames


The partial S1 gene sequence of IBV was generated from 13 flocks. These sequences and related metadata were deposited in the GenBank database under accession numbers OQ199494 and OQ199507 (Additional file 4). To distinguish between vaccine and field type viruses, nucleotide identity between the detected IBV sequences and IBV vaccine sequences was used. The detected IBVs were considered a vaccine type when their nucleotide identity with a vaccine strain was 99% or more, and were considered a field challenge when their nucleotide identity was less than 99% [35]. Adopting [39] nomenclature, the phylogenetic analysis showed the presence of three GI lineages. Eleven, two and one sequence respectively fell in the GI-13, GI-23 and GI-1 lineages. GI-13 were 793B vaccine-like sequences and homologous to either 4/91 (eight sequences) or CR88 vaccine strains (three sequences). The GI-23 sequences shared 95% of nucleotide identity with both the Israeli variant 2 strain (AF093796) and the Libyan GI-23 strains detected in 2012 (Figure 6).

Figure 6
figure 6

Phylogenetic tree of IBV. The tree was constructed using IBV HVR3 nucleotide sequences. The sequences from the present study were labelled with red triangle and reference sequences of the different IBV lineages available in GenBank were unlabelled. The evolutionary history was inferred using the Maximum Likelihood method with General Time Reversible model and 1000 bootstrap replications in MEGA X. A total of 1875 positions were included in the final dataset. Only bootstrap values higher than 60 were conserved


In this study, we monitored over a 3 year-period spontaneous acute respiratory cases on broiler farms operated by the same rearing company in Tunisia. Our aim was to investigate the major respiratory pathogens involved in respiratory diseases in the northeastern region of Tunisia and their longitudinal circulation.

The design of our study, which involved targeting a large panel of respiratory pathogens and sampling clinically diseased birds during the early stage of a respiratory outbreak, enabled us to detect at least one viral pathogen per flock with high viral loads. We even detected aMPV (Ct 15), which is known to be difficult to detect due to its short shedding window (Additional file 3) [40].

Importantly, we found that H9N2 LPAIV was the most detected respiratory pathogen (54%) on the farms monitored. H9N2 LPAIV was enzootic on the monitored farms in 2018 and 2020, with outbreaks occurring during spring and summer. However, we cannot confirm the seasonality of H9N2 LPAIV due to the interruption of field monitoring that occurred during the study period. Further investigation is therefore required. The phylogenetic analysis revealed that the H9N2 LPAIV viruses detected in 2018 and 2020 were monophyletic, belonged to the G1 lineage and clustered with Northern and Western African viruses introduced from the Middle East in 2016 [5]. In addition, the H9N2 LPAIV strain reported here was distant from the old Tunisian G1 lineage that emerged in 2009 and that was detected until 2016 on poultry farms, and until 2018 in migratory birds [3, 4]. The phylogenetics of H9N2 LPAIV suggest a new introduction of H9N2 LPAIV in Tunisia. This study is the first to report the circulation of the Northern and Western African clade of H9N2 LPAIV in poultry in Tunisia. Interestingly, during the same period in 2018, the first detection of H9N2 LPAIV coincided with the detection of a very closely related virus (99%) from the water of a lagoon in the same region of Tunisia [6]. Our findings support the hypothesis proposed by Larbi et al. [6] that the transmission of H9N2 LPAIV occurs from poultry farms to migratory birds via their contamination of the environment as strains isolated from the water and the migratory birds on the same site were different. Considering the close phylogenetic relatedness between the H9N2 LPAIV introduced in Tunisia and those circulating in neighboring countries during the same period, we assume that the new introduction occurred through the trade of poultry or poultry products.

The second most detected virus on the farms monitored was IBV (53%), with an enzootic circulation over the study period. The phylogenetic analysis showed the dominance of GI-13 vaccine-like strains (793B). This finding is common as GI-13 live virus vaccines are used widely on poultry. The rolling infection of birds with vaccine viruses within a flock, when the initial coverage is incomplete, can lead to vaccination reactions, thus complicating diagnostic efforts [41]. Given the epidemiological context in Tunisia, an infection with a rolling IBV vaccine strain could be co-infected by H9N2 LPAIV which leads to the exacerbation of the clinical signs [42, 43]. In this study, we report the detection of the GI-23 lineage (variant 2) for the first time in Tunisia. Over the 3 years of the study period, variant 2 was detected twice, 5 months apart, on the same farm. The lack of spread of the GI-23 virus to the other farms monitored and its occasional detection suggest a cross protection of Massachusetts and 793B combination against variant 2, as also described in previous studies [44, 45]. However, we cannot draw conclusions about the prevalence of this lineage in Tunisia; active and continuous surveillance of this variant therefore is required.

The detection of aMPV started in February 2020 and lasted until June 2020. During this period, aMPV was detected on both monitored and external farms (28%) in the different governorates, indicating the emergence and circulation of a new aMPV in broilers. Lachheb et al. [46] also reported the detection of an aMPV field strain in 2019 in broilers. Unfortunately, we were not able to determine the phylogenetic relatedness of the previously described virus with the detected aMPVs during this study as the sequence was not published. The phylogenetic analysis grouped the detected aMPVs in two different clusters within the subtype B and outside of the vaccine strains cluster, suggesting the circulation of two different field strains. We noticed that these aMPVs held polymorphic features in their SH-G non-coding region which do not exist in other aMPVs. Noncoding regions have different biological roles, including the regulation of viral replication, viral persistence, host immune evasion, and cellular transformation [47]. These insertions consequently should not be neglected, and pathogenicity studies are required.

NDV genotype VII.2 was detected during the same period in two flocks on the same monitored farm and on one external farm, showing an epizootic circulation of this virus. This sporadic detection of NDV could be explained by the reinforcement of NDV control measures. After genotype VII.2 emerged in 2013 in Tunisia, it was recommended to add NDV live vaccine boosters every 10 days in broilers and to use NDV vector vaccines. The phylogenetic analysis of NDV showed that the detected viruses are identical and share the most recent common ancestor with previously reported viruses in Tunisian (KU175357) and Libya [48].

Bacterial pathogens also were highly detected, even though samples were collected in the acute respiratory phase showing an early bacterial superinfection. To the best of our knowledge, mycoplasma prevalence in poultry in Tunisia has not yet been studied. Based on the global prevalence of MG and MS, the low detection of MG (one flock) in contrast with the high detection of MS (33%) is in line with their reported occurrence in all poultry populations. However, our results are not in line with the global occurrence of MG and MS in broiler poultry, which are reported to be approximately the same (around 25%) [49]. The lower rate of MG compared to MS could be explained by control measures that are oriented more toward MG than MS in breeding birds, thereby reducing MG vertical transmission in broilers. Considering these results, we are tempted to say that control measures in Tunisia have reduced the occurrence of MG on broiler farms. Nevertheless, more studies on mycoplasma in different poultry populations are needed to monitor the effectiveness of MG and MS control measures and adapt them accordingly. Surprisingly, the detection rate of ORT was very high (100%) in the 12 flocks tested. The screening of ORT was included in our respiratory panel in 2020 and only flocks sampled that year were tested. Even so, these results draw our attention to the likely contribution of this pathogen to respiratory disease in broilers. Our results are consistent with the few existing ORT studies in other countries suggesting the underdiagnosis of this pathogen [21, 50,51,52].

Interestingly, co-infections between respiratory pathogens were very frequent as 61% of the flocks studied were positive for two or more pathogens. Our results are consistent with the results of similar studies [15, 16] confirming the complexity of respiratory disease and its diagnosis in the field. Viral co-infections occurred in 16 of the 39 flocks studied; of these, H9N2 LPAIV was the major contributor (15/16). The dominance of co-infections by H9N2 LPAIV could be explained by the high detection rate of H9N2 LPAIV in the flocks studied, and could explain the clinical signs associated with H9N2 LPAIV positive flocks. In these co-infections, the major viral combinations were H9N2 + IBV (6/16) and H9N2 + aMPV (8/16). Several experimental studies have shown that co-infections with H9N2 LPAIV and IBV consistently result in more severe clinical signs than single infections with either H9N2 LPAIV or IBV, regardless of the virulence of the IBV strain [17, 42, 43]. This synergy between H9N2 LPAIV and IBV could be explained either by the severe inflammatory response induced by IBV [20], or by the trypsin-like serine protease encoded by IBV which probably facilitates the cleavage activation of the hemagglutinin of the H9N2 LPAIV virus, and thus increases its replication [53, 54]. In contrast, little is known about co-infection mechanisms between H9N2 LPAIV and aMPV as to the best of our knowledge this co-infection has not been previously reported. Clinical manifestations of aMPV in broilers are absent or mild, so co-infections with H9N2 LPAIV could explain the clinical manifestation of flocks positive for aMPV [55]. All of the flocks co-infected with H9N2 and aMPV also were positive for ORT, making it even more difficult to understand the sequential order of co-infections and the synergistic effect between the three pathogens. Only dual co-infections between ORT and H9N2 or aMPV and ORT were previously studied by experimental infections. Pan et al. [21] showed that in broilers, a co-infection involving H9N2 and ORT, either simultaneously or primed by ORT, increased mortality by 60–70% respectively comparing to a single infection with H9N2 LPAIV. Similarly, Marien et al. [56] showed that a prior infection of ORT with aMPV enabled the reproduction of respiratory disease in specific pathogen free (SPF) turkeys, which was not possible by the inoculation of ORT alone via the natural infection route.

In conclusion, this study brought to light the emergence of three new respiratory viruses in Tunisia, as well as several respiratory co-infections. The introduction of a new G1 lineage of H9N2 LPAIV, GI-23 lineage of IBV and subtype B of aMPV was shown, which highlights the importance of active and continuous surveillance and points to an urgent need to adjust control measures to the current situation. The high frequency of the co-infection between H9N2 LPAIV and aMPV suggests that it would be relevant to study the mechanisms of a likely synergism between these two viruses. The complexity of respiratory disease demonstrated in the study emphasizes the need to adopt a comprehensive diagnostic approach for multifactorial respiratory diseases.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author upon reasonable request. The dataset supporting the conclusions of this article is included within the article.


  1. Iqbal M, Yaqub T, Reddy K, McCauley JW (2009) Novel genotypes of H9N2 influenza A viruses isolated from poultry in Pakistan containing NS genes similar to highly pathogenic H7N3 and H5N1 viruses. PLoS One 4:e5788.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Pusch E, Suarez D (2018) The multifaceted zoonotic risk of H9N2 avian influenza. Vet Sci 5:82.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Arbi M, Souiai O, Rego N, Larbi I, Naya H, Ghram A, Houimel M (2020) Historical origins and zoonotic potential of avian influenza virus H9N2 in Tunisia revealed by Bayesian analysis and molecular characterization. Arch Virol 165:1527–1540.

    Article  CAS  PubMed  Google Scholar 

  4. Tombari W, Nsiri J, Larbi I, Guerin JL, Ghram A (2011) Genetic evolution of low pathogenecity H9N2 Avian influenza viruses in Tunisia: acquisition of new mutations. Virol J 8:467.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Zecchin B, Minoungou G, Fusaro A, Moctar S, Ouedraogo-Kaboré A, Schivo A, Salviato A, Marciano S, Monne I (2017) Influenza A(H9N2) virus, Burkina Faso. Emerg Infect Dis 23:2118–2119.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Larbi I, Ghedira K, Arbi M, Butcher GD, Rego N, Naya H, Tougorti H, Lachhab J, Behi IE, Nsiri J, Ghram A (2022) Phylogenetic analysis and assessment of the pathogenic potential of the first H9N2 avian influenza viruses isolated from wild birds and Lagoon water in Tunisia. Virus Res 322:198929.

    Article  CAS  PubMed  Google Scholar 

  7. Awuni JA, Bianco A, Dogbey OJ, Fusaro A, Yingar DT, Salviato A, Ababio PT, Milani A, Bonfante F, Monne I (2019) Avian influenza H9N2 subtype in Ghana: virus characterization and evidence of co-infection. Avian Pathol 48:470–476.

    Article  CAS  PubMed  Google Scholar 

  8. Fusade-Boyer M, Djegui F, Batawui K, Byuragaba DK, Jones JC, Wabwire-Mangeni F, Erima B, Atim G, Ukuli QA, Tugume T, Dogno K, Adjabli K, Nzuzi M, Adjin R, Jeevan T, Rubrum A, Go-Maro W, Kayali G, McKenzie P, Webby RJ, Ducatez MF (2021) Antigenic and molecular characterization of low pathogenic avian influenza A(H9N2) viruses in Sub-Saharan Africa from 2017 through 2019. Emerg Microbes Infect 10:753–761.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Jeevan T, Darnell D, Gradi EA, Benali Y, Kara R, Guetarni D, Rubrum A, Seiler PJ, Crumpton JC, Webby RJ, Derrar F (2019) A(H9N2) influenza viruses associated with chicken mortality in outbreaks in Algeria 2017. Influenza Other Respir Viruses 13:622–626.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. EL Houadfi M, Fellahi S, Nassik S, Guérin J-L, Ducatez MF (2016) First outbreaks and phylogenetic analyses of avian influenza H9N2 viruses isolated from poultry flocks in Morocco. Virol J 13:140.

    Article  CAS  PubMed  Google Scholar 

  11. Jallow MM, Fall A, Barry MA, Diop B, Sy S, Goudiaby D, Fall M, Enouf V, Niang MN, Dia N (2020) Genetic characterization of the first detected human case of low pathogenic avian influenza A/H9N2 in sub-Saharan Africa, Senegal. Emerg Microbes Infect 9:1092–1095.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Kaboudi K, Nciri J, Amara A, Laarbi I, Moalla N, Bouzouaia M, Ghram A (2016) Histopathological and molecular diagnosis of infectious laryngotracheitis in Tunisia: first report. Int J Livest Res 6:34–45.

    Article  Google Scholar 

  13. Lachheb J, Turki A, Nsiri J, Fathallah I, El behi I, Larbi I, Ghram A, (2019) Molecular characterization of a unique variant of avian infectious bronchitis virus in Tunisia. Poult Sci 98:4338–4345.

    Article  CAS  PubMed  Google Scholar 

  14. Gowthaman V, Singh SD, Dhama K, Srinivasan P, Saravanan S, Murthy TRGK, Ramakrishnan MA (2017) Molecular survey of respiratory and immunosuppressive pathogens associated with low pathogenic avian influenza H9N2 subtype and virulent Newcastle disease viruses in commercial chicken flocks. J Poult Sci 54:179–184.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Hassan KE, Eady MF, Eawah AAA, Luttermann C, Parvin R, Shany S, Beer M, Harder T (2019) Respiratory disease due to mixed viral infections in poultry flocks in Egypt between 2017 and 2018: upsurge of highly pathogenic avian influenza virus subtype H5N8 since 2018. Transbound Emerg Dis 68:21–36.

    Article  CAS  PubMed  Google Scholar 

  16. Roussan DA, Haddad R, Khawaldeh G (2008) Molecular survey of avian respiratory pathogens in commercial broiler chicken flocks with respiratory diseases in Jordan. Poult Sci 87:444–448.

    Article  CAS  PubMed  Google Scholar 

  17. Belkasmi SFZ, Fellahi S, Touzani CD, Faraji FZ, Maaroufi I, Delverdier M, Guérin J-L, Fihri OF, El Houadfi M, Ducatez MF (2020) 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 49:21–28.

    Article  CAS  PubMed  Google Scholar 

  18. Barbour EK, Mastori FA, Nour AMA, Shaib HA, Jaber LS, Yaghi RH, Sabra A, Sleiman FT, Sawaya RK, Niedzwieck A, Tayeb IT, Kassaify ZG, Rath M, Harakeh S, Barbour KE (2009) Standardisation of a new model of H9N2/Escherichia coli challenge in broilers in the Lebanon. Vet Ital 45:317–322

    PubMed  Google Scholar 

  19. Bonfante F, Cattoli G, Leardini S, Salomoni A, Mazzetto E, Davidson I, Haddas R, Terregino C (2017) Synergy or interference of a H9N2 avian influenza virus with a velogenic Newcastle disease virus in chickens is dose dependent. Avian Pathol 46:488–496.

    Article  CAS  PubMed  Google Scholar 

  20. Kong L, You R, Zhang D, Yuan Q, Xiang B, Liang J, Lin Q, Ding C, Liao M, Chen L, Ren T (2021) Infectious bronchitis virus infection increases pathogenicity of H9N2 avian influenza virus by inducing severe inflammatory response. Front Vet Sci 8:824179.

    Article  PubMed  Google Scholar 

  21. Pan Q, Liu A, Zhang F, Ling Y, Ou C, Hou N, He C (2012) Co-infection of broilers with Ornithobacterium rhinotracheale and H9N2 avian influenza virus. BMC Vet Res 8:104.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Sakai T, Ishii A, Segawa T, Takagi Y, Kobayashi Y, Itou T (2015) Establishing conditions for the storage and elution of rabies virus RNA using FTA® cards. J Vet Med Sci 77:461–465.

    Article  CAS  PubMed  Google Scholar 

  23. Croville G, Foret C, Heuillard P, Senet A, Delpont M, Mouahid M, Ducatez MF, Kichou F, Guerin J-L (2018) Disclosing respiratory co-infections: a broad-range panel assay for avian respiratory pathogens on a nanofluidic PCR platform. Avian Pathol 47:253–260.

    Article  CAS  PubMed  Google Scholar 

  24. Fouchier RAM, Bestebroer TM, Herfst S, Van Der Kemp L, Rimmelzwaan GF, Osterhaus ADME (2000) Detection of influenza A viruses from different species by PCR amplification of conserved sequences in the matrix gene. J Clin Microbiol 38:4096–4101

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Kho CL, Mohd-Azmi ML, Arshad SS, Yusoff K (2000) Performance of an RT-nested PCR ELISA for detection of newcastle disease virus. J Virol Methods 86:71–83.

    Article  CAS  PubMed  Google Scholar 

  26. Callison SA, Hilt DA, Boynton TO, Sample BF, Robison R, Swayne DE, Jackwood MW (2006) Development and evaluation of a real-time Taqman RT-PCR assay for the detection of infectious bronchitis virus from infected chickens. J Virol Methods 138:60–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Cecchinato M, Lupini C, Pogoreltseva OSM, Listorti V, Mondin A, Drigo M, Catelli E (2013) Development of a real-time RT-PCR assay for the simultaneous identification, quantitation and differentiation of avian metapneumovirus subtypes A and B. Avian Pathol 42:283–289.

    Article  CAS  PubMed  Google Scholar 

  28. Mahmoudian A, Kirkpatrick NC, Coppo M, Lee S-W, Devlin JM, Markham PF, Browning GF, Noormohammadi AH (2011) Development of a SYBR green quantitative polymerase chain reaction assay for rapid detection and quantification of infectious laryngotracheitis virus. Avian Pathol 40:237–242.

    Article  CAS  PubMed  Google Scholar 

  29. Wetzel AN, Lefevre KM, Raviv Z (2010) Revised Mycoplasma synoviae vlhA PCRs. Avian Dis 54:1292–1297.

    Article  PubMed  Google Scholar 

  30. García M, Ikuta N, Levisohn S, Kleven SH (2005) Evaluation and comparison of various PCR methods for detection of Mycoplasma gallisepticum infection in chickens. Avian Dis 49:125–132.

    Article  PubMed  Google Scholar 

  31. Shabat M, Meir R, Haddas R, Lapin E, Shkoda I, Raibstein I, Perk S, Davidson I (2010) Development of a real-time TaqMan RT-PCR assay for the detection of H9N2 avian influenza viruses. J Virol Methods 168:72–77.

    Article  CAS  PubMed  Google Scholar 

  32. Hoffmann E, Stech J, Guan Y, Webster RG, Perez DR (2001) Universal primer set for the full-length amplification of all influenza A viruses. Arch Virol 146:2275–2289.

    Article  CAS  PubMed  Google Scholar 

  33. Esmaelizad M, Mayahi V, Pashaei M, Goudarzi H (2017) Identification of novel newcastle disease virus sub-genotype VII-(j) based on the fusion protein. Arch Virol 162:971–978.

    Article  CAS  PubMed  Google Scholar 

  34. Cecchinato M, Catelli E, Lupini C, Ricchizzi E, Clubbe J, Battilani M, Naylor CJ (2010) Avian metapneumovirus (AMPV) attachment protein involvement in probable virus evolution concurrent with mass live vaccine introduction. Vet Microbiol 146:24–34.

    Article  CAS  PubMed  Google Scholar 

  35. Worthington KJ, Currie RJW, Jones RC (2008) A reverse transcriptase-polymerase chain reaction survey of infectious bronchitis virus genotypes in Western Europe from 2002 to 2006. Avian Pathol 37:247–257.

    Article  CAS  PubMed  Google Scholar 

  36. BLAST: Basic Local Alignment Search Tool.

  37. Katoh K, Rozewicki J, Yamada KD (2019) MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform 20:1160–1166.

    Article  CAS  PubMed  Google Scholar 

  38. Dimitrov KM, Abolnik C, Afonso CL, Albina E, Bahl J, Berg M, Briand F-X, Brown IH, Choi K-S, Chvala I, Diel DG, Durr PA, Ferreira HL, Fusaro A, Gil P, Goujgoulova GV, Grund C, Hicks JT, Joannis TM, Torchetti MK, Kolosov S, Lambrecht B, Lewis NS, Liu H, Liu H, McCullough S, Miller PJ, Monne I, Muller CP, Munir M et al (2019) Updated unified phylogenetic classification system and revised nomenclature for newcastle disease virus. Infect Genet Evol 74:103917.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Valastro V, Holmes EC, Britton P, Fusaro A, Jackwood MW, Cattoli G, Monne I (2016) S1 gene-based phylogeny of infectious bronchitis virus: an attempt to harmonize virus classification. Infect Genet Evol 39:349–364.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Aldous EW, Mynn JK, Banks J, Alexander DJ (2003) A molecular epidemiological study of avian paramyxovirus type 1 (Newcastle disease virus) isolates by phylogenetic analysis of a partial nucleotide sequence of the fusion protein gene. Avian Pathol 32:239–257.

    Article  CAS  PubMed  Google Scholar 

  41. Jackwood MW, Lee D-H (2017) Different evolutionary trajectories of vaccine-controlled and non-controlled avian infectious bronchitis viruses in commercial poultry. PLoS One 12:e0176709.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Haghighat-Jahromi M, Asasi K, Nili H, Dadras H, Shooshtari AH (2008) Coinfection of avian influenza virus (H9N2 subtype) with infectious bronchitis live vaccine. Arch Virol 153:651–655.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Ismail ZM, EL-Deeb AH, EL-Safty MM, Hussein HA (2018) Enhanced pathogenicity of low-pathogenic H9N2 avian influenza virus after vaccination with infectious bronchitis live attenuated vaccine. Vet World 11:977–985.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Cook JK, Orbell SJ, Woods MA, Huggins MB (1999) Breadth of protection of the respiratory tract provided by different live-attenuated infectious bronchitis vaccines against challenge with infectious bronchitis viruses of heterologous serotypes. Avian Pathol 28:477–485.

    Article  CAS  PubMed  Google Scholar 

  45. de Wit JJ, Boelm GJ, van Gerwe TJWM, Swart WAJM (2013) The required sample size in vaccination-challenge experiments with infectious bronchitis virus, a meta-analysis. Avian Pathol 42:9–16.

    Article  CAS  PubMed  Google Scholar 

  46. Lachheb J, Bouslama Z, Nsiri J, Badr C, Al Gallas N, Souissi N, Khazri I, Larbi I, Kaboudi K, Ghram A (2022) Phylogenetic and phylodynamic analyses of subtype-B metapneumovirus from chickens in Tunisia. Poult Sci 102:102253.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Tycowski KT, Guo YE, Lee N, Moss WN, Vallery TK, Xie M, Steitz JA (2015) Viral noncoding RNAs: more surprises. Genes Dev 29:567–584.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Kammon A, Heidari A, Dayhum A, Eldaghayes I, Sharif M, Monne I, Cattoli G, Asheg A, Farhat M, Kraim E (2015) Characterization of avian influenza and newcastle disease viruses from poultry in Libya. Avian Dis 59:422–430.

    Article  PubMed  Google Scholar 

  49. Chaidez-Ibarra MA, Velazquez DZ, Enriquez-Verdugo I, Castro del Campo N, Rodriguez-Gaxiola MA, Montero-Pardo A, Diaz D, Gaxiola SM (2022) Pooled molecular occurrence of Mycoplasma gallisepticum and Mycoplasma synoviae in poultry: a systematic review and meta-analysis. Transbound Emerg Dis 69:2499–2511.

    Article  PubMed  Google Scholar 

  50. Al-Hasan BA, Alhatami AO, Abdulwahab HM, Bustani GS, Alkuwaity EAW (2021) The first isolation and detection of Ornithobacterium rhinotracheale from swollen head syndrome-infected broiler flocks in Iraq. Vet World 14:2346.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Bhuiyan ZA, Ali MZ, Moula MM, Bary MA, Arefin N, Giasuddin M, Khan ZUM (2019) Seroprevalence of major avian respiratory diseases in broiler and sonali chicken in selected areas of Bangladesh. J Adv Vet Anim Res 6:561–566.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Chansiripornchai N, Wanasawaeng W, Sasipreeyajan J (2007) Seroprevalence and identification of Ornithobacterium rhinotracheale from broiler and broiler breeder flocks in Thailand. Avian Dis 51:777–780.[777:SAIOOR]2.0.CO;2

    Article  PubMed  Google Scholar 

  53. Liu DX, Brierley I, Brown TDK (1995) Identification of a trypsin-like serine proteinase domain encoded by ORF la of the coronavirus IBV. In: Talbot PJ, Levy GA (eds) Corona- and related viruses: current concepts in molecular biology and pathogenesis. Springer, Boston, pp 405–411

    Chapter  Google Scholar 

  54. Ng LFP, Liu DX (2000) Further characterization of the coronavirus infectious bronchitis virus 3C-like proteinase and determination of a new cleavage site. Virology 272:27–39.

    Article  CAS  PubMed  Google Scholar 

  55. Hosseini H, Ziafati Kafi Z, Malekan M, Ghafouri SA, Fallah Mehrabadi MH, Sadri N, Hojabr Rajeoni A, Ghalyanchilangeroudi A (2021) Molecular characterization of circulating avian metapneumovirus, subgroup B, in broiler chickens, Iran, 2016–2018. Iran J Vet Res 22:217–221.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Marien M, Decostere A, Martel A, Chiers K, Froyman R, Nauwynck H (2005) Synergy between avian pneumovirus and Ornithobacterium rhinotracheale in turkeys. Avian Pathol 34:204–211.

    Article  PubMed  Google Scholar 

Download references


The authors would like to thank the Tunisian rearing company for offering us the possibility to collect samples and data from its farms. We also thank Professor Moncef Bouzouaia, DVM, from the Ecole Nationale Vétérinaire de Sidi Thabet, Tunisia for his help and effort dedicated to the project.


This study was performed in the framework of the PhD program of Adam Jbenyeni, hosted by the Ecole Nationale Vétérinaire de Toulouse (ENVT) and funded by Ceva Santé Animale S.A., with the cofunding by the Association Nationale de la Recherche et de la Technologie (Grant: CIFRE—2019/1605).

Author information

Authors and Affiliations



JLG, GC and CC: contributed to the conception and the design of the study. Samples collected by AJ. PCR analyses were conducted by AJ and were interpreted by AJ and GC. Data was curated by AJ. AJ drafted the first manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jean-Luc Guérin.

Ethics declarations

Ethics approval and consent to participate

All farms enrolled in the study followed their own management practices. No protocol approval from the ethics committee of Ecole Nationale Vétérinaire de Toulouse (ENVT) was required.

Competing interests

The authors declare that they have no competing interests.

Additional information

Handling editor: Stéphane Biacchesi.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Oligonucleotides set used for the PCR screening of respiratory pathogens. The targeted gene, the primers’ nucleotide sequences and the amplicon size for each one of the eight screened pathogens were listed in this file.

Additional file 2:

Investigated flocks’ information. The date of sampling, the farm and governorate of origin, the size and the age of the investigated flocks were showed in this file.

Additional file 3:

The qPCR Ct values of detected pathogens. The Age, the mortality, and the qPCR Ct values of the detected pathogens per flocks were listed in this file.

Additional file 4:

GenBank accession numbers of nucleotide sequences generated during this study. Virus sequences generated in this study were submitted to GenBank database under the accession numbers listed in this additional file.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jbenyeni, A., Croville, G., Cazaban, C. et al. Predominance of low pathogenic avian influenza virus H9N2 in the respiratory co-infections in broilers in Tunisia: a longitudinal field study, 2018–2020. Vet Res 54, 88 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: