Cocirculation of Four Infectious Bronchitis Virus Lineages in Broiler Chickens in the Eastern Region of Saudi Arabia from 2012 to 2014

Al-Mubarak, Abdullah I. A.; Al-Kubati, Anwar A. G.

doi:https://doi.org/10.1155/2020/6037893

Veterinary Medicine International

On this page

Abstract Introduction Materials and Methods Results Discussion Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 6037893 | https://doi.org/10.1155/2020/6037893

Cocirculation of Four Infectious Bronchitis Virus Lineages in Broiler Chickens in the Eastern Region of Saudi Arabia from 2012 to 2014

Abdullah I. A. Al-Mubarak¹and Anwar A. G. Al-Kubati²

Academic Editor: Douglas Morck

Received31 Jul 2019

Revised17 Nov 2019

Accepted11 Dec 2019

Published31 Mar 2020

Abstract

Avian infectious bronchitis virus (IBV) is an evolving and dynamic virus that causes major economic losses for the poultry industry worldwide. Continuous evolution and emergence of new variants of this virus are the major challenges for controlling the disease with routine vaccination. Successful vaccination usually requires the use of a homologous vaccine, which in turn necessitates continuous investigation of the circulating strains. Herein, we performed a reverse transcriptase-polymerase chain reaction- (RT-PCR-) based investigation in broiler chicken flocks of the Eastern Region of Saudi Arabia. IBV was detected in 36.5% of the tested flocks (42 out of 115) from January 2012 to March 2014. Direct sequencing of hypervariable region-3 (HVR-3) of the Spike (S)-1 gene was performed, followed by phylogenetic analysis to determine the circulating IBV genotypes. Four lineages appear to coexist in this region, including the GI-13 or 4/91 IBV (31%), GI-16 or CK/CH/LDL/97I IBV (28.6%), GI-1 or Mass IBV (19%), and GI-23 or Middle East IBV (21.4%). The latter lineage include two subgroups: IS/720/99 IBV (16.7%) and IS/Variant2/98 IBV (4.7%). Some of the detections made in the 4/91 and Mass lineages are expected to belong to the vaccine strains. Lineages without a homologous vaccine in use (CK/CH/LDL/97I and Middle East) represent 50% of the isolates recovered in this study. Based on identity with the vaccine sequences, field observations, and frequent detection, these two lineages appear to be out of coverage of the IBV vaccines used in Saudi Arabia. This is the first time to identify Middle East lineage (IS/720/99 IBV and IS/Variant2/98 IBV) in the Eastern Region of Saudi Arabia.

1. Introduction

Infectious bronchitis (IB) is an acute highly contagious respiratory disease of chickens (Gallus gallus domesticus) that is present in virtually all regions with an intensive poultry industry. In addition to respiratory involvement, IB may also affect the urogenital and/or alimentary tracts. IB causes significant economic losses by adversely affecting weight gain and the quality and quantity of produced eggs as well as predisposing birds to secondary infections [1].

The causative agent, infectious bronchitis virus (IBV), belongs to group III Coronaviruses of family Coronaviridae. IBV is an enveloped virus with a single-stranded positive-sense RNA genome of approximately 27.6 kb. Four structural proteins (spike (S), envelope (E), membrane (M), and nucleocapsid (N)) were reported to be encoded by the IBV genome [2]. The glycosylated S protein is posttranslationally cleaved by cellular proteases into S1, which forms the spike bulb, and S2, which anchors S1 to the virion membrane. The S1 and S2 subunits have been shown to mediate viral attachment and membrane fusion, respectively. Virus neutralizing epitopes are located in the first and third quarters of the S1 glycoprotein [2, 3]. Three hypervariable regions (HVR)-1, 2, and 3 were reported to be located in the S1 subunit between amino acids 38 to 67, 91 to 141, and 274 to 387, respectively [4]. Genetic variability of these HVRs due to mutations and/or recombination was incriminated in the continuous emergence of new IBV serotypes/genotypes [5]. In this regard, substitution of 10 to 15 amino acids (2-3%) in the S1 subunit has been shown to result in generation of a new serotype [6].

Large numbers of IBV serotypes have been reported worldwide, some of which remain restricted to particular geographic areas, whereas others tend to spread over the majority of the world [5, 7]. Vaccination is the primary measure used to control IBV infection, but caution must be taken when choosing a vaccine strain(s) that should be homologous with the circulating IBV serotype(s). Vaccination-challenge studies have shown that protection by heterologous vaccines ranges between very poor to moderate depending on the criteria used. Poor cross protection may be induced by as little as a 5% difference in the S1 amino acid sequence, while most IBV serotypes differ from each other by 20–25% and up to 50% [6]. Accordingly, S1 nucleotide sequence similarity and genotyping approaches have been used successfully to predict cross-protection and circulating serotypes, with some exceptions [8, 9]. Recently, IBV was genotyped based on S1 sequences into 6 major genotypes (Genotype I (GI) to GVI] containing 32 lineages [10].

Few reports have evaluated IBV infection in Saudi Arabia. IBV infection in this region extends back to 1984 when Zwaagstra and others reported an IBV isolate without specifying its serotype [11]. The 4/91 IBV serotype (GI-13 lineage) was detected in Saudi Arabia in 2000 [12]. Vaccine strains belonging to this serotype and the Mass serotype (GI-1 lineage), such as H120 and Ma5, are in use to control IB occurrence; however, IB still occurs due to the CH/CK/LDL/97I (GI-16) lineage, the Middle East (GI-23) lineage (IS/720/99 and IS/Variant2/98), and D274 (GI-12) lineage [13, 14]. Strains similar to those circulating in Egypt, India, China, and Italy have also been reported in the Eastern Region of Saudi Arabia [15]. In the present study, we report the IBV lineages circulating in 115 broiler chicken flocks in the Eastern Region of Saudi Arabia.

2. Materials and Methods

2.1. Ethical Approval

All experimental procedures and management conditions used in this study were approved by the Ethics Committee at King Faisal University, Saudi Arabia.

2.2. Sampling and Nucleic Acid Extraction

The study was conducted between January 2012 and March 2014. Samples were collected from 115 broiler chicken flocks in the Eastern Region (Al-Hassa and Dammam) of Saudi Arabia by visiting broilers farms, broilers slaughterhouse, and the poultry clinic at the Veterinary Teaching Hospital of King Faisal University. Tissue samples, such as trachea, lung, spleen, kidney, and bursa of Fabricius, were collected from birds showing signs of respiratory illness. Although not always successful, attempts were made to collect various organs and a complete history from the sampled flocks. For each flock, each organ was pooled together from all sampled birds. Tissue samples were collected in 10 volumes of phosphate-buffered saline (PBS) containing gentamicin and nystatin (50 μg of each/ml). The samples were stored at −80°C until being homogenized. Homogenization was performed using the Omni International Ceramic bead kit and the BioSpecMini-Beadbeater-16 instrument. For each organ, viral nucleic acids were extracted from homogenized tissues with the IQEasy Plus Viral DNA/RNA Extraction Kit (Cat #17153, iNtRON Biotechnology, South Korea) according to the manufacturer’s instructions. The extracted nucleic acids were stored at −80°C until being tested by RT-PCR. The isolates were named according to the Cavanagh nomenclature system [16] that includes the country/sample name/sampling year.

2.3. Synthesis of cDNA

First-strand cDNA synthesis was performed using a Reverse Transcription System (Cat #A3500, Promega, USA) according to the manufacturer’s instructions. The reaction was performed in a final volume of 20 μl containing the following ingredients: 5 mM MgCl₂, 1x RT buffer, 1 mM of each dNTP, 1 U/μl of Recombinant RNasin® Ribonuclease Inhibitor, 0.75 U/μl of AMV reverse transcriptase, 25 ng/μl of Random Primers, and 50 ng/μl of total RNA. The reaction mixture was incubated at 42°C for 45 minutes, followed by heat inactivation of the enzyme at 95°C for 5 minutes and 4°C for 5 minutes. The cDNA was diluted with nuclease-free water to a final volume of 50 μl and stored at −80°C until being tested by PCR.

2.4. N Gene PCR

For diagnostic purposes, a nested PCR targeting the nucleocapsid (N) gene of IBV was carried out. The reaction was performed using the previously published primers N784, N1145, N791, and N1129 [17], which were obtained from Integrated DNA Technologies (IDT, Coralville, IA, USA) and presented in Table 1. The nested PCR reactions were performed using the Go Taq® Green Master Mix (Cat #M7122, Promega) in a final volume of 25 μl containing the primers at a final concentration of 0.8 μM with 2 μl of the template cDNA. The Bio-Rad T-100 Thermal Cycler was used to perform an amplification reaction that consisted of 94°C for 5 minutes and 35 cycles of 94°C for 45 seconds, 60°C for 1 minute, and 72°C for 2 minutes, followed by a final extension step of 72°C for 7 minutes. A similar thermal profile was used for the second round of the nested PCR except that the annealing step was performed at 53°C for 1 minute and the extension step lasted for 1 minute and 10 minutes for the cyclic and final extension steps, respectively. The second-round PCR products were separated on a 1.5% agarose gel in Tris-Borate-EDTA (TBE) buffer. Red Safe (Cat #21141, iNtRON Biotechnology) was used as a DNA stain. Positive samples were subjected to the S1 gene PCR.

2.5. S1 Gene PCR

For sequencing purposes, another nested PCR targeting HVR-3 of the IBV-S1 gene was performed. We used previously published primers, including SX1, SX2, SX3, and SX4 [18], which were obtained from Integrated DNA Technologies (IDT, Coralville, IA) and presented in Table 1. A PCR was performed using i-StarMAX II (Cat #25174, iNtRON Biotechnology) in a final volume of 50 μl. One microliter (final concentration 0.2 μM) of each primer, 2 μl of template cDNA, and 21 μl of nuclease-free water were added to 25 μl of the master mix. Thereafter, the reactions were treated in a similar manner as the second round of the N gene PCR except that the annealing temperature was 55°C for 1 minute for both rounds of the nested PCR. The second-round PCR product was separated on a 1.5% low melting point agarose gel in TBE buffer. The targeted bands were excised, and the amplicons were purified from the agarose gel using the Wizard SV Gel and PCR Cleanup System (Cat #A1460, Promega) according to the manufacturer’s instructions. Purified amplicons were sequenced by Macrogen Sequencing Service (South Korea).

2.6. Isolation

The isolation protocol described by Gelb and Jackwood [19] was performed to propagate the virus in embryonated Specific Pathogen Free (SPF) eggs (Nile SPF, Egypt). Isolation was attempted from samples (Table 2) that gave positive results in the N gene PCR but turned either negative or weakly positive in the S1 gene PCR. Before inoculation, homogenized tissue samples were centrifuged at 1000 ×g for 10 minutes. The centrifugation-derived supernatant was passed through a 0.2 μm sterile nylon syringe filter (Thermo Scientific, Nalgene®, Cat #195-2520, USA). A 200 μl volume of filtrated liquid was inoculated into the allantoic cavity in 10-day-old embryonated eggs. Forty-eight hours postinoculation, inoculated eggs were placed at 4°C overnight. Allantoic fluids (AFs) were collected and used as inocula, at a rate of 100 μl/egg, for additional two passages. AFs collected from the third passage were retested by RT-PCR.

2.7. Sequence Analysis

The sequence analysis was performed using the Molecular Evolutionary Genetic Analysis (MEGA) version X software [20]. Sequences were aligned using the Clustal W method. The aligned sequences were trimmed to a length of 327 nt. (from nt. 718 to 1044 based on the 793/B S1 gene sequence (GB#Z83979)). A phylogenetic tree was constructed using the maximum likelihood method with bootstrap values of 1000 replicates. Twenty-nine IBV reference sequences were included (Figure 1). The overall mean distance was used as an indicator of sequence divergence for the Saudi isolates and was shown as the number of different bases per sequence from averaged overall sequence pairs. Sequence identity was calculated using the BioEdit version 7.1.7 software [21]. A BLAST search was performed for all sequences to determine the most related sequence in GenBank.

2.8. Vaccination History

Vaccination histories for the flocks in which IBV was detected are presented in Table 2. Regarding the 13 flocks in which 4/91 IBV was detected, four flocks received only the 4/91 vaccine, one flock received both the Mass and 4/91 vaccines, and four flocks received only the Mass vaccine twice. Histories were not available for the other four flocks. For the 12 flocks in which the CK/CH/LDL/97I IBV was detected, eight flocks received the Mass vaccine either once (n = 5) or twice (n = 3), one flock received the 4/91 vaccine once, and histories were not available for the additional three flocks. In the 8 flocks that were positive for the Mass IBV, five flocks received the Mass vaccine, and histories were not available for the other three flocks. Concerning the seven flocks positive for the IS/720/99 IBV, 2 flocks received the H120 vaccine and one received the 4/91 vaccine; histories were not available for the other 4 flocks. Finally, the two flocks in which IS/Variant2/98 was detected had received the Mass vaccine.

2.9. Differentiation between Field and Vaccine Viruses

The methods used in this study could not definitely determine the fraction contributed by the vaccine viruses; however, removal of sequences with 100% nucleotide identity to the reference vaccine sequences was used only as a guide, as previously reported by [14, 18, 22]. Comparisons were performed with the two Mass vaccines (H120 and Ma5) and with the 4/91 vaccine.

2.10. GenBank Accession Numbers

The sequences of the 42 IBV isolates from the present study have been deposited in the GenBank database with accession numbers MH648687 to MH648727 and MH449644, as shown in Figure 1 and Table 2.

3. Results

3.1. Detection Rate and Genotyping

During the study period, tissue samples were collected from 115 broiler flocks in the Eastern Region of Saudi Arabia. Samples from forty-two flocks (36.52%) were positive for IBV in the N and S1 genes RT-PCRs (Table 2) either directly or after isolation in SPF eggs. Isolation were performed for 7 samples including SA/IH1/12, SA/IH5/12, SA/IH11/13, SA/IH19/13, SA/IH21/13, SA/IC8/13, and SA/IC15/13. Phylogenetic analysis based on partial S1 gene sequences of these forty-two Saudi isolates and twenty-nine reference sequences showed segregation of these isolates into four lineages (Figure 1 and Supplementary Table 1) including GI-13 (4/91), GI-16 (CK/CH/LDL/97I), GI-1 (Mass), and GI-23 (Middle East) that include two subdivisions: IS/720/99 and IS/Variant2/98 IBVs.

3.2. GI-13 or 4/91 Lineage

The first cluster included thirteen sequences that were grouped with 4/91, IS/Variant1, IS/1366, and other reference sequences. A BLAST search for 4/91-like sequences showed that the 4/91 strain (GB#AF093794) was the most similar sequence, with identities ranging between 96.9% and 100% over the sequenced region of the S1 gene. Sequences of this lineage showed an overall mean distance of 3.000 ± 0.720 nucleotides and an average nucleotide identity of 99.5% with the 4/91 vaccine (Table 3).

3.3. GI-16 or CK/CH/LDL/97I Lineage

The second lineage included twelve sequences that clustered around the CK/CH/LDL/97I, CK/CH/SCYA/101, Q1, and J2 reference sequences. A BLAST search for CK/CH/LDL/97I-like sequences showed that Chinese Q1 was the most similar sequence, with nucleotide identities ranging between 99.0% and 100% over the sequenced part of the S1 gene. The overall mean distance was 0.667 ± 0.316 nucleotides. The average nucleotide identity of these sequences with the IBV vaccines was 79.5% with the H120 and Ma5 vaccines and 81.7% with the 4/91 vaccine (Table 3).

3.4. GI-1 or Mass Lineage

Eight Saudi sequences were included in this lineage together with the Mass serotype strains, including the H120, Ma5, M41, Beaudette, and Connecticut reference sequences. A BLAST search for Mass-like sequences revealed H120 as the most similar sequence, with identities ranging between 97.8% and 100% over the sequenced part of the S1 gene. The overall mean distance of this group was 2.000 ± 0.674 nucleotides. The presented sequences showed a 99.7% average nucleotide identity with the H120 and Ma5 vaccines (Table 3).

3.5. GI-23 or Middle East Lineage

The fourth lineage included nine Saudi isolates and further subdivided into two subgroups: IS/720/00 IBV and IS/Variant2/98 IBV. IS/720/99 IBV included seven Saudi isolates that were assembled around the Middle East Sul/01/09, IR/59/2010, IS/885/00, and IS/720/99 reference sequences. A BLAST search for IS/720/99-like sequences revealed that IR/IS720/H140/11 was the most related sequence, with identities ranging between 97.8% and 99.3% over the sequenced region of the S1 gene. The overall mean distance for these sequences was 4.952 ± 1.463 nucleotides, and their nucleotide identities were on average 81.9% with both the Mass vaccines and the 4/91 vaccine (Table 3). In the other subgroup, IS/Variant2/98 IBV, only two Saudi sequences were located with the IS/Variant2/98 and IS/1494/06 reference sequences. A BLAST search for the IS/Variant2/98-like sequences showed that IS/1494/06 was the most similar sequence, with identities ranging between 97.8% and 100% over the sequenced part of the S1 gene. The overall mean distance was 7.000 ± 2.552 nucleotides. The average nucleotide identity was 82.2% with the Mass vaccines and 83.0% with the 4/91 vaccine (Table 3).

3.6. Field and Vaccine Viruses

At least two types of IBV vaccines (Mass and 4/91) are in use in Saudi Arabia. Hence, isolates in these two lineages could belong to vaccine viruses. We compared the partial S1-gene sequences of isolates in these lineages with those of the relevant vaccines. Nine out of the 13 sequences in 4/91 lineage and 6 out of the 8 sequences in Mass lineage showed complete nucleotide identities with corresponding vaccine sequences over the sequenced part of S1 gene (Table 2).

4. Discussion

The present investigation documents the IBV lineages circulating in broiler chicken flocks of the Eastern Region of Saudi Arabia during the period from January 2012 to March 2014. For IBV detection, we used a validated and highly sensitive nested PCR that targeted the conserved part of the N gene and was used frequently in other studies [13, 17, 23, 24].

In the current work, IBV was detected in 36.5% of the tested flocks. A slightly higher detection rate (42.7%) was reported in Riyadh between October 2009 and May 2010 [25]. Caution should be taken in interpreting these results because samples were collected only from flocks showing signs of respiratory illness. Similar investigations revealed detection rates as high as 74% in south Iraq [26], 64% in Egypt [27], 60% in Jordan [28], 59% in western Europe [18], 52% in Iran [29], 37.5% in eastern Iran [30], 34% in Russia [22], and 32% in middle-south Iraq [31]. Testing flocks regardless of their clinical picture are expected to reveal lower detection rates as previously reported by Chen and others [32]. Furthermore, vaccine isolates represent a fraction of this percentage.

Comparison of sequences of the present Saudi isolates with those of the IBV vaccines revealed that six and nine isolates from the Mass and 4/91 lineages, respectively, had vaccine-identical sequences (Table 2). Comparable vaccine detection rates were reported in other RT-PCR-based investigations, where vaccines represented approximately half or more of the total detections made in the vaccines lineages [14, 18, 22]. If the vaccine-identical sequences were removed, the CK/CH/LDL/97I lineage would become the dominant lineage in the Eastern Region of Saudi Arabia.

In Saudi Arabia, the CK/CH/LDL/97I (GI-16) lineage was first reported in Riyadh in 2009–2010 [25]. Then, this strain was reported in 2011 in Saudi Arabia, Jordan, and Iraq by Ababneh and colleagues [13]. In the present study, the CK/CH/LDL/97I lineage was detected at a rate of 28.6% (12/42) between January 2012 and March 2014 in the Eastern Region of Saudi Arabia. A detection rate of 17% for this lineage was reported from Saudi Arabia during an overlapping period from 2009 to 2014 [14]. Conversely, this strain was either absent or had a low prevalence rate in surrounding countries because it was not detected in studies in Egypt, Lebanon, Jordan, Kuwait, the United Arab Emirate (UAE) and Oman [14, 33], Iraq [26, 31], and Iran [29, 30].

Our data showed that 8 of the 12 flocks in which the CK/CH/LDL/97I lineage was detected received only the Mass vaccine and that one received only the 4/91 vaccine; histories were not available for the remaining three flocks. This finding was not surprising because the original CK/CH/LDL/97I strain was isolated from an H120-vaccinated flock. Likewise, a vaccination-challenge experiment showed that heterologous vaccines did not provide sufficient protection against CK/CH/LDL/97I IBV. In contrast, complete protection was provided by a homologous vaccine derived from the same type and attenuated by passaging in embryonated hen eggs [34]. In Saudi Arabia, the high prevalence of this lineage in broiler chickens appears to be facilitated by an immunity gap created by use of only heterologous vaccines. A similar scenario was suggested for spread of the IT-02 IBV strain in Spanish 4/91-vaccinated flocks [9].

The 4/91 (G1-13) lineage was first reported in Saudi Arabia in 2000 [12]. In the present investigation, the 4/91 lineage was detected at a rate of 31% in the Eastern Region of Saudi Arabia between 2012 and 2014. Similarly, it was the most prevalent lineage with a detection rate of 43% in Saudi Arabia during the period from 2009 to 2014 [14] and had reported detection rates of 40% and 50% in Iraq [26, 31], 64% and 81% in Oman, and 85% in the UAE [14, 33]. The 4/91 lineage was also the second most prevalent lineage in Iran (19% and 27%), Egypt (17%), and Lebanon (13%) but was undetectable in Jordan and Kuwait [14, 29, 30].

In the current study, approximately two-thirds of the 4/91 lineage isolates (9 out of 13) showed 100% sequence identity with the 4/91 vaccine virus. These isolates most likely belong to the vaccine virus. This speculation is supported by the history of recent use of the 4/91 vaccine in 4 out of these 9 flocks. Conversely, four of these nine flocks received only the Mass vaccine, and the vaccination history was not available for the ninth flock. Keeping in mind that approximately 70% of the 4/91-like sequences were vaccine-identical, and this lineage contains some isolates that displayed relatively high sequence divergence (Table 3). Heavy vaccine application may contribute to such divergence, as previously suggested by Lee and Jackwood [35]. In contrast, the CK/CH/LDL/97I lineage showed the lowest divergence, which could further support our justification of the low immune pressure exerted by the used vaccines on CK/CH/LDL/97I lineage viruses.

The Mass (GI-1) lineage was detected at a rate of 19% in the present study. A similar rate (17%) was reported by Ganapathy and colleagues (2015) in Saudi Arabia during the 2009–2014 period. A relatively lower prevalence was reported in surrounding countries, including Egypt (8%), Lebanon (4%), the UAE (4%), Iraq (6%), Oman (6% and 3%), and Iran (3%), and it was undetectable in Jordan and Kuwait [14, 26, 29, 33].

Out of the eight isolates recovered for the Mass lineage, six showed vaccine-identical sequences and were likely to belong to the Mass vaccine strain. In favor of this opinion, four out of these six vaccine-identical isolates were obtained from flocks that received the Mass vaccines; histories were not available for the remaining two flocks. On the contrary, two out of the eight Mass isolates showed incomplete identity with Mass vaccines (SA/IH20/13 and SA/IC79/13). SA/IH20/13 showed 99.6% nucleotide identity with the Mass vaccine sequence and was recovered from a Mass-vaccinated flock (Table 2). Similarly, SA/IH3/12 (4/91-like isolate) showed 99.6% nucleotide identity with the 4/91 vaccine sequence and was recovered from a 4/91-vaccinated flock. Whether these isolates belong to field viruses or to mutant vaccine viruses is not clear. Nucleotide mutations and subsequent amino acid substitutions may affect vaccine protectivity or even lead to vaccine reversion. An escape mutant of the Mass M41 strain was reported to occur due to a single nucleotide mutation at position 134, leading to an amino acid change at position 45 of HVR-1 of the S1 gene [36].

In the current study, the IS/720/99 IBV of the Middle East (GI-23) lineage was detected in 7 flocks (16.7%) in the Eastern Region of Saudi Arabia during the period from 2012 to 2014. This subdivision was first reported in Saudi Arabia in 2010 by Ganapathy and others and up to 2014 was detected in 11% of IBV-positive samples [14]. In the present study, all seven isolates were recovered in 2013 (Table 2). Similarly, most of the detections by Ganapathy and colleagues were made in 2013, whereas this virus was almost undetectable in 2011 and 2012 [14]. This difference in the detection rate may reflect the time of introduction and ongoing spread of this genotype inside the country. The IS/720/99 IBV was shown to have a low prevalence in Oman (3%) and was not detected in Jordan, Lebanon, Kuwait, the UAE, Iraq, and Iran. In contrast, IS/720/99 IBV was the most prevalent (71%) in Egypt [14, 26, 29–31, 33].

Vaccination histories were available for 3 of the 7 flocks in which the IS/720/99 IBV was detected. Two flocks received the Mass vaccine, and one flock received the 4/91 vaccine. A vaccination-challenge experiment showed that the H120 vaccine provided little protection against IS/720/99 [37]. Over the sequenced region in the present study, genetic relatedness between the IS/720/99-like isolates and the Ma5 or 4/91 vaccines was not better than that with H120. Therefore, we can rationally expect poor cross-protection with these heterologous vaccines.

The other subdivision of the Middle East (GI-23) lineage detected in this study was the IS/Variant2/98 IBV. Our data showed that the IS/Variant2/98 IBV had the lowest prevalence (4.7%) in Eastern Region of Saudi Arabia. From 2009 to 2014, a prevalence of 11% was reported in Saudi Arabia and 4% in Egypt, and the strain was not detected in the UAE [14]. Similarly, variable prevalence rates of 21% and 3% were reported in Oman during overlapping periods from 2009 to 2014 and during 2012, respectively [14, 33]. In contrast, this strain was reported to be the most prevalent in the neighboring north countries, including prevalence rates of 100% in Jordan and Kuwait; 82% in Lebanon; 67%, 70%, and 75% in Iran; and 47% in Iraq [14, 29–31, 38].

The two isolates within the IS/Variant2/98 subdivision reported in this study were recovered from two Mass-vaccinated broiler flocks. Based on reported field observations of H120-vaccinated flocks [39] and sequence similarity of the detected isolates with available vaccines (Mass and 4/91), which showed ≤83% sequence identity (Table 3), these vaccines seem unlikely to be able to provide protection against IS/variant2-like strains.

As stated above, considerable variation existed between the detection rates reported from different countries and even from the same countries, for example, the difference in the detection rates of 4/91 lineage in Saudi Arabia reported by this study and Ganapathy and others (2015) and the above mentioned difference in the detection rates of IS/Variant2/98 IBV in Oman. Likewise, other lineages were reported in surrounding countries but not in Saudi Arabia or vice versa; for example, the QX (or GI-19) lineage was detected at rates of 7% and 8% in Iran and 9% in Iraq [29–31]. Similarly, Ganapathy and colleagues (2015) detected the D274 (or GI-12) lineage at a frequency of once out of 236 samples during the period from 2009 to 2014, whereas we did not detect this lineage at the Eastern Region of Saudi Arabia during the period from 2012 to 2014. Similarly, the CK/CH/LDL/97I lineage was prevalent in Saudi Arabia but not in surrounding countries. Several factors may contribute to this inconsistency, including but not limited to the following: (1) differences in the sampled population, such as those that arise from sampling different birds (broilers, layers, breeders etc.); (2) serotype(s) of the used vaccines, frequency of their application, and the protection they provide against the circulating genotypes; (3) the effect of migratory birds and localization of the sampled flocks on the migration route [40, 41]; (4) applied biosecurity measures; and (5) exchange of poultry products between countries and regions [13]. Additionally, testing different tissues would also contribute to variation in the detection rates. IBV vaccine and field viruses, or their genome/antigen, were detected in variety of tissues [42–44] and reported to persist in some viscera, particularly cecal tonsils [18, 45–48]. Consequently, testing the latter tissue would increase the detection rate on the one hand and increase the fraction contributed by vaccine viruses, on the other hand, as previously speculated by Callison and others [49].

In the present study, IBV genotyping was performed by direct sequencing of amplicon produced by nested RT-PCR that targets the HVR-3 of the S1 gene [18]. Since being introduced, this method has been frequently used for IBV genotyping [14, 29, 50–52]. Advantages and disadvantages of full and partial sequencing of IBV-S1 gene was previously reported [53]. According to the phylogenetic analysis of the partial (HVR-3) S1 gene sequences published by Valastro et al. [10], all of the lineages involved in the present study (G1-1, -13, -16, and -23) formed distinct and well-supported clusters (Shimodaira–Hasegawa-like test support values of ≥0.93). Indeed, Valastro et al. specify 8 lineages (GI-5, -7, -10, -18, -22, -24, -25, and -27) that loss their unity in phylogeny when using the partial 342 nt. encoding the HVR3 [10]. However, the main limitation of the current work is the partial sequence of the S1 gene that was used for phylogenetic analysis. Full S1 gene sequence would provide more accurate analysis especially for differentiation between field and vaccine viruses.

In conclusion, our data show that IB is a prevalent disease in broiler chicken flocks of the Eastern Region of Saudi Arabia. Four IBV lineages were found to cocirculate in this region. Half of the detected isolates belong to two of these lineages (CK/CH/LDL/97I and Middle East) against which no homologues vaccines are available, and these lineages are expected to keep causing IBV outbreaks and economic losses in the absence of such vaccines. Further studies to elucidate the serotypic, pathotypic, and protectotypic properties of the endemic types are required. Continuous monitoring of the IBV types circulating in Saudi Arabia as well as continuous update of the vaccination strategy seem to be necessary to control IBV infection.

Data Availability

The epidemiological data used to support findings of this article are included within the article. The sequence data used to support findings of this article were deposited in GenBank under accession numbers shown in the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to thank the Deanship of Scientific Research, King Faisal University, Saudi Arabia, for the moral and financial support of this work (Grant number: 130026). The authors also thank Dr. Mustafa Ababneh and Mr. Mohamed Saleem.

Supplementary Materials

Supplementary Table 1: nucleotide identity matrix based on partial S1 gene sequences from seventeen IBV isolates and the most related reference sequences. The reference sequences are shown in bold font. Sequences with 100% identity are not shown, including the following: SA/IH9/13, SA/IH11/13, SA/IH12/13, SA/IH15/13, SA/IH16/13, SA/IC10/13, SA/IC71/13, SA/IC78/13, and SA/IC80/13 resembling SA/IH7/12; SA/IH19/13, SA/IH25/13, SA/IH26/13, SA/IH27/13, SA/IH28/13, SA/IH29/13, SA/IC7/13, and SA/IC16/13 resembling SA/IH8/12; SA/H21/13, SA/IH22/13, SA/IC3/13, SA/IC24/13, and SA/IC72/13 resembling SA/IH6/12; SA/H14/13 resembling SA/H10/13; and SA/C8/13 and SA/IC9/13 resembling SA/IC6/13. (Supplementary Materials)

References

D. Cavanagh and J. Gelb, “Infectious bronchitis,” in Diseases of Poultry, Y. M. SAIF, A. M. Fadly, J. R. Glisson, L. R. Mcdougald, L. K. Nolan, and D. E. Swayne, Eds., pp. 117–135, Blackwell Publishing, Ames, Iowa, USA, 12th edition, 2008.
View at: Google Scholar
D. Cavanagh, “Coronavirus avian infectious bronchitis virus,” Veterinary Research, vol. 38, no. 2, pp. 281–297, 2007.
View at: Publisher Site | Google Scholar
M. W. Jackwood, D. A. Hilt, S. A. Callison, C.-W. Lee, H. Plaza, and E. Wade, “Spike glycoprotein cleavage recognition site analysis of infectious bronchitis virus,” Avian Diseases, vol. 45, no. 2, pp. 366–372, 2001.
View at: Publisher Site | Google Scholar
A. Rimondi, M. I. Craig, A. Vagnozzi, G. König, M. Delamer, and A. Pereda, “Molecular characterization of avian infectious bronchitis virus strains from outbreaks in Argentina (2001–2008),” Avian Pathology, vol. 38, no. 2, pp. 149–153, 2009.
View at: Publisher Site | Google Scholar
J. J. de Wit, J. K. A. Cook, and H. M. J. F. van der Heijden, “Infectious bronchitis virus variants: a review of the history, current situation and control measures,” Avian Pathology, vol. 40, no. 3, pp. 223–235, 2011.
View at: Publisher Site | Google Scholar
D. Cavanagh, “Severe acute respiratory syndrome vaccine development: experiences of vaccination against avian infectious bronchitis coronavirus,” Avian Pathology, vol. 32, no. 6, pp. 567–582, 2003.
View at: Publisher Site | Google Scholar
J. Ignjatovic and S. Sapats, “Avian infectious bronchitis virus,” Revue Scientifique et Technique de l’OIE, vol. 19, no. 2, pp. 493–508, 2000.
View at: Publisher Site | Google Scholar
B. S. Ladman, A. B. Loupos, and J. Gelb JR., “Infectious bronchitis virus S1 gene sequence comparison is a better predictor of challenge of immunity in chickens than serotyping by virus neutralization,” Avian Pathology, vol. 35, no. 2, pp. 127–133, 2006.
View at: Publisher Site | Google Scholar
R. Dolz, J. Pujols, G. Ordóñez, R. Porta, and N. Majó, “Molecular epidemiology and evolution of avian infectious bronchitis virus in Spain over a fourteen-year period,” Virology, vol. 374, no. 1, pp. 50–59, 2008.
View at: Publisher Site | Google Scholar
V. Valastro, E. C. Holmes, P. Britton et al., “S1 gene-based phylogeny of infectious bronchitis virus: an attempt to harmonize virus classification,” Infection, Genetics and Evolution, vol. 39, pp. 349–364, 2016.
View at: Publisher Site | Google Scholar
K. A. Zwaagstra, B. A. van der Zeijst, and J. G. Kusters, “Rapid detection and identification of avian infectious bronchitis virus,” Journal of Clinical Microbiology, vol. 30, no. 1, pp. 79–84, 1992.
View at: Publisher Site | Google Scholar
D. Cavanagh, J.-P. Picault, R. E. Gough, M. Hess, K. Mawditt, and P. Britton, “Variation in the spike protein of the 793/B type of infectious bronchitis virus, in the field and during alternate passage in chickens and embryonated eggs,” Avian Pathology, vol. 34, no. 1, pp. 20–25, 2005.
View at: Publisher Site | Google Scholar
M. Ababneh, A. E. Dalab, S. Alsaad, and M. Al-Zghoul, “Presence of infectious bronchitis virus strain CK/CH/LDL/97I in the Middle East,” ISRN Veterinary Science, vol. 2012, Article ID 201721, 6 pages, 2012.
View at: Publisher Site | Google Scholar
K. Ganapathy, C. Ball, and A. Forrester, “Genotypes of infectious bronchitis viruses circulating in the Middle East between 2009 and 2014,” Virus Research, vol. 210, pp. 198–204, 2015.
View at: Publisher Site | Google Scholar
M. G. Hemida, M. A. Al-Hammadi, A. H. S. Daleb, and C. R. Gonsalves, “Molecular characterization and phylogenetic analyses of virulent infectious bronchitis viruses isolated from chickens in Eastern Saudi Arabia,” VirusDisease, vol. 28, no. 2, pp. 189–199, 2017.
View at: Publisher Site | Google Scholar
D. Cavanagh, “A nomenclature for avian coronavirus isolates and the question of species status,” Avian Pathology, vol. 30, no. 2, pp. 109–115, 2001.
View at: Publisher Site | Google Scholar
A. Farsang, C. Ros, L. H. M. Renström, C. Baule, T. Soós, and S. Belák, “Molecular epizootiology of infectious bronchitis virus in Sweden indicating the involvement of a vaccine strain,” Avian Pathology, vol. 31, no. 3, pp. 229–236, 2002.
View at: Publisher Site | Google Scholar
K. J. Worthington, R. J. W. Currie, and R. C. Jones, “A reverse transcriptase-polymerase chain reaction survey of infectious bronchitis virus genotypes in Western Europe from 2002 to 2006,” Avian Pathology, vol. 37, no. 3, pp. 247–257, 2008.
View at: Publisher Site | Google Scholar
J. J. Gelb and M. W. Jackwood, “Infectious broncitis,” in A Laboratory Manual for the Isolation and Identification of Avian Pathogens, D. E. Swayne, J. R. Glisson, M. W. Jackwood, J. E. Pearson, and W. R. Reed, Eds., pp. 169–174, American Association of Avian Pathologists, University of Pennsylvania, Kennett Square, PA, USA, 4th edition, 1998.
View at: Google Scholar
S. Kumar, G. Stecher, M. Li, C. Knyaz, and K. Tamura, “Mega X: molecular evolutionary genetics analysis across computing platforms,” Molecular Biology and Evolution, vol. 35, no. 6, pp. 1547–1549, 2018.
View at: Publisher Site | Google Scholar
T. A. Hall, “BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT,” Nucleic Acids Symposium Series, vol. 41, pp. 95–98, 1999.
View at: Google Scholar
E. V. Ovchinnikova, Y. A. Bochkov, L. O. Shcherbakova et al., “Molecular characterization of infectious bronchitis virus isolates from Russia and neighbouring countries: identification of intertypic recombination in the S1 gene,” Avian Pathology, vol. 40, no. 5, pp. 507–514, 2011.
View at: Publisher Site | Google Scholar
Z. H. Mahmood, R. R. Sleman, and A. U. Uthman, “Isolation and molecular characterization of Sul/01/09 avian infectious bronchitis virus, indicates the emergence of a new genotype in the Middle East,” Veterinary Microbiology, vol. 150, no. 1-2, pp. 21–27, 2011.
View at: Publisher Site | Google Scholar
V. Sumi, S. D. Singh, K. Dhama, V. Gowthaman, R. Barathidasan, and K. Sukumar, “Isolation and molecular characterization of infectious bronchitis virus from recent outbreaks in broiler flocks reveals emergence of novel strain in India,” Tropical Animal Health and Production, vol. 44, no. 7, pp. 1791–1795, 2012.
View at: Publisher Site | Google Scholar
Y. M. AL-Hammad, A. I. Al-Afaleq, and M. H. A. Mohamed, “Molecular survey and phylogenic analysis of infectious bronchitis virus (IBV) circulating among chicken flocks in Riyadh province, Saudi Arabia,” Journal of Animal and Veterinary Advances, vol. 13, pp. 1002–1008, 2014.
View at: Google Scholar
W. Seger, A. Ghalyanchi Langeroudi, V. Karimi, O. Madadgar, M. Vasfi Marandi, and M. Hashemzadeh, “Prevalence of avian infectious bronchitis virus in broiler chicken farms in south of Iraq, 2014 – 2015,” Vet Res Forum, vol. 7, no. 4, pp. 317–321, 2016a.
View at: Google Scholar
K. Selim, A. S. Arafa, H. A. Hussein, and A. A. El-Sanousi, “Molecular characterization of infectious bronchitis viruses isolated from broiler and layer chicken farms in Egypt during 2012,” International Journal of Veterinary Science and Medicine, vol. 1, no. 2, pp. 102–108, 2013.
View at: Publisher Site | Google Scholar
D. A. Roussan, W. S. Totanji, and G. Y. Khawaldeh, “Molecular subtype of infectious bronchitis virus in broiler flocks in Jordan,” Poultry Science, vol. 87, no. 4, pp. 661–664, 2008.
View at: Publisher Site | Google Scholar
A. Hamadana, A. Ghalyanchilangerodui, M. Hashemzadeh et al., “Genotyping of avian infectious bronchitis viruses in Iran (2015–2017) reveals domination of IS-1494 like virus,” Virus Research, vol. 240, pp. 101–106, 2017.
View at: Publisher Site | Google Scholar
A. Ghalyanchi-Langeroudi, V. Karimi, A. Jannat et al., “Genotyping of infectious bronchitis viruses in the East of Iran, 2015,” Iranian Journal of Virology, vol. 9, no. 2, pp. 31–35, 2015.
View at: Publisher Site | Google Scholar
W. Seger, A. Ghalyanchilangeroudi, V. Karimi, O. Madadgar, M. V. Marandi, and M. Hashemzadeh, “Genotyping of infectious bronchitis viruses from broiler farms in Iraq during 2014-2015,” Archives of Virology, vol. 161, no. 5, pp. 1229–1237, 2016b.
View at: Publisher Site | Google Scholar
H. W. Chen, Y. P. Huang, and C. H. Wang, “Identification of intertypic recombinant infectious bronchitis viruses from slaughtered chickens,” Poultry Science, vol. 89, no. 3, pp. 439–446, 2010.
View at: Publisher Site | Google Scholar
T. Al-Shekaili, M. Baylis, and K. Ganapathy, “Molecular detection of infectious bronchitis and avian metapneumoviruses in Oman backyard poultry,” Research in Veterinary Science, vol. 99, pp. 46–52, 2015.
View at: Publisher Site | Google Scholar
S. Liu, X. Zhang, Y. Wang et al., “Evaluation of the protection conferred by commercial vaccines and attenuated heterologous isolates in China against the CK/CH/LDL/97I strain of infectious bronchitis coronavirus,” The Veterinary Journal, vol. 179, no. 1, pp. 130–136, 2009.
View at: Publisher Site | Google Scholar
C.-W. Lee and M. W. Jackwood, “Origin and evolution of Georgia 98 (GA98), a new serotype of avian infectious bronchitis virus,” Virus Research, vol. 80, no. 1-2, pp. 33–39, 2001.
View at: Publisher Site | Google Scholar
D. Cavanagh, P. J. Davis, and A. P. A. MOCKETT, “Amino acids within hypervariable region 1 of avian coronavirus IBV (Massachusetts serotype) spike glycoprotein are associated with neutralization epitopes,” Virus Research, vol. 11, no. 2, pp. 141–150, 1988.
View at: Publisher Site | Google Scholar
J. Gelb Jr, Y. Weisman, B. S. Ladman, and R. Meir, “S1 gene characteristics and efficacy of vaccination against infectious bronchitis virus field isolates from the United States and Israel (1996 to 2000),” Avian Pathology, vol. 34, no. 3, pp. 194–203, 2005.
View at: Publisher Site | Google Scholar
Y. Saadat, M. H. Bozorgmehri Fard, S. Charkhkar, H. Hosseini, N. Shaikhi, and B. Akbarpour, “Molecular characterization of infectious bronchitis viruses isolated from broiler flocks in Bushehr province, Iran: 2014 - 2015,” Veterinary Research Forum, vol. 8, pp. 195–201, 2017.
View at: Google Scholar
S. Kahya, F. Coven, S. Temelli, A. Eyigor, and K. T. Carli, “Presence of IS/1494/06 genotype-related infectious bronchitis virus in breeder and broiler flocks in Turkey,” Ankara Üniversitesi Veteriner Fakültesi Dergisi, vol. 60, no. 1, pp. 27–31, 2013.
View at: Publisher Site | Google Scholar
L. Sun, G.-H. Zhang, J.-W. Jiang et al., “A Massachusetts prototype like coronavirus isolated from wild peafowls is pathogenic to chickens,” Virus Research, vol. 130, no. 1-2, pp. 121–128, 2007.
View at: Publisher Site | Google Scholar
L. A. Hughes, C. Savage, C. Naylor, M. Bennett, J. Chantrey, and R. Jones, “Genetically diverse coronaviruses in wild bird populations of northern England,” Emerging Infectious Diseases, vol. 15, no. 7, pp. 1091–1094, 2009.
View at: Publisher Site | Google Scholar
A. G. Ambali and R. C. Jones, “Early pathogenesis in chicks of infection with an enterotropic strain of infectious bronchitis virus,” Avian Diseases, vol. 34, no. 4, pp. 809–817, 1990.
View at: Publisher Site | Google Scholar
A. Abdel-Moneim, H. M. Madbouly, and M. F. El-Kady, “In vitro characterization and pathogenesis of Egypt/Beni-Suef/01; a novel genotype of infectious bronchitis virus,” Beni Suef Veterinary Medical Journal, vol. 15, no. 2, 2005.
View at: Google Scholar
W. Fan, H. Wang, Y. Zhang, and E. Al, “Comparative dynamic distribution of avian infectious bronchitis virus M41, H120, and SAIBK strains by quantitative real-time RT-PCR in SPF chickens,” Bioscience, Biotechnology, and Biochemistry, vol. 76, no. 12, pp. 2255–2260, 2012.
View at: Google Scholar
J. J. De Wit, “Detection of infectious bronchitis virus,” Avian Pathology, vol. 29, no. 2, pp. 71–93, 2000.
View at: Publisher Site | Google Scholar
S. Naqi, K. Gay, P. Patalla, S. Mondal, and R. Liu, “Establishment of persistent avian infectious bronchitis virus infection in antibody-free and antibody-positive chickens,” Avian Diseases, vol. 47, no. 3, pp. 594–601, 2003.
View at: Publisher Site | Google Scholar
I. R. Alvarado, P. Villegas, J. El-Attrache, and M. W. Jackwood, “Detection of Massachusetts and Arkansas serotypes of infectious bronchitis virus in broilers,” Avian Diseases, vol. 50, no. 2, pp. 292–297, 2006.
View at: Publisher Site | Google Scholar
F. Amjad, N. Hassan, and H. Arsalan, “Detection the 4/91 strain of infectious bronchitis virus in testicular tissue from experimentally infected rooster by reverse transcription-polymerase chain reaction,” African Journal of Agricultural Research, vol. 4, pp. 1093–1096, 2009.
View at: Google Scholar
S. A. Callison, D. A. Hilt, T. O. Boynton et al., “Development and evaluation of a real-time Taqman RT-PCR assay for the detection of infectious bronchitis virus from infected chickens,” Journal of Virological Methods, vol. 138, no. 1-2, pp. 60–65, 2006.
View at: Publisher Site | Google Scholar
B. Vidović, M. Šekler, D. Rogan et al., “Molecular characterization of infectious bronchitis virus strains isolated from vaccinated flocks in Serbia and their comparison with the isolated strains from neighboring countries,” Kafkas Universitesi Veteriner Fakultesi Dergisi, vol. 24, 2018.
View at: Google Scholar
N. Sadri, A. Ghalyanchilangeroudi, M. H. Fallah Mehrabadi et al., “Genotyping of Avian Infectious Bronchitis Virus in Afghanistan (2016-2017): the first report,” Iranian Journal of Veterinary Research, vol. 20, no. 1, pp. 60–63, 2019.
View at: Google Scholar
M. H. Wibowo, T. E. Ginting, and W. ASMARA, “Molecular characterization of pathogenic 4/91-like and QX-like infectious bronchitis virus infecting commercial poultry farms in Indonesia,” Veterinary World, vol. 12, no. 2, pp. 277–287, 2019.
View at: Publisher Site | Google Scholar
C. Ball, B. Manswr, A. Forrester, J. Chantrey, and K. Ganapathy, “Full versus partial S1 gene sequencing of classical and variant avian infectious bronchitis viruses,” in Proceedings of the 9th International Symposium on Avian Corona- and Pneumoviruses, pp. 156–161, Utrecht, Netherlands, June 2016.
View at: Google Scholar

Copyright

Copyright © 2020 Abdullah I. A. Al-Mubarak and Anwar A. G. Al-Kubati. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1955

Downloads

1315

Citations