key: cord-0887587-w8kpsi2s
authors: Modiri Hamadan, Amir; Ghalyanchilangeroudi, Arash; Hashemzadeh, Masoud; Hosseini, Hossein; Karimi, Vahid; Yahyaraeyat, Ramak; Najafi, Hamideh
title: Genotyping of Avian infectious bronchitis viruses in Iran (2015–2017) reveals domination of IS-1494 like virus
date: 2017-08-15
journal: Virus Res
DOI: 10.1016/j.virusres.2017.08.002
sha: abb2a92fe176ca007e19c2d2ba2dffdbb4cc0908
doc_id: 887587
cord_uid: w8kpsi2s

Avian infectious bronchitis virus (IBV) is causing major economic losses to the poultry industry. The analysis of the S1 gene has been used to determine IBV genotype. The aim of this study was genotyping of IBVs circulating among the Iranian broiler flocks in the period between 2015 to 2017. Trachea samples from 278 broiler flocks were collected from broiler farms in eight provinces of Iran. After Real-time RT-PCR, IBV-positive samples were further characterized based on S1 gene. The results of the Real-time RT-PCR showed that 52.16% of flocks were IBV positive. Four genotypes were detected and the frequency of occurrence rates of IS-1494-like, 793/B, QX and Massachusetts IBV genotypes were 70.34%, 19.31%, 7.58% and 2.75%, respectively. Sequence analysis revealed that nucleotide identities within IS-1494-like group ranged between 98.86–100%, while each of the QX, Massachusetts and 793/B groups were 98.05–100%, 98.20–100% and 93.29–100% respectively. These results show that the IS-1494-like IBV is the dominant IBV genotype in Iran. Proper control strategies are essential to overcoming the high frequency of occurrence of IS-1494-like IBV. The phylogenetic relationship of the strains with respect to different sequences and geographical regions displayed complexity and diversity. Further studies are needed and should include the isolation and full-length molecular characterization of IBV in Iran.

The Avian infectious bronchitis virus (IBV) continues to be one of the most economically-impactful diseases in poultry despite intensive control (Jackwood, 2012) . The virus affects mainly the respiratory tract, causing tracheal rales, sneezing, coughing, reduced weight gain and mortality, particularly in broiler chickens (Cook et al., 2012) . The etiologic agent is a gamma coronavirus. IBV is an enveloped virus approximately 120 nm in diameter with a single-stranded, non-segmented, positive sense RNA genome. The RNA genome of 27.6 kb has at least 10 open reading frames (ORF), from 5′ to 3′ are as follows: 5′-1a-1b-S (S1,S2)-3a,b,c(E)-M-5a,b-N-Poly(A)-3′. The genome encodes four main structural proteins (spike glycoprotein (S), small envelope protein (E), membrane glycoprotein (M), and nucleocapsid protein (N)) and numerous non-structural proteins (Liu et al., 2009) . The S1 spike glycoprotein carries virus-neutralizing and serotype-specific antigenic determinants, which can evolve rapidly especially within the three hypervariable regions (HVR) (Abro et al., 2012) . Variations in the S1 protein have led to the emergence of numerous variants over time, some of which persist and spread, while others disappear after a short time. (Jackwood, 2012) . Improving vaccination efficiency and the understanding of IBV evolution requires knowledge of the IBV variants in the region (Najafi et al., 2016) . In Iran, several IBV genotypes have been reported; the first isolation of IBV in Iranian chicken flocks was reported in 1994 (Aghakhan et al., 1994) . Later, several Iranian researchers identified the 793/B genotype (Seyfi Abad Shapouri et al., 2004) . Genotypes of IBV strains isolated in Iran during 2014-2015 were classified into seven distinct phylogenetic groups (Massachusetts (Mass), 793/B, IS-1494, IS-720, QX, IR-1, and IR-2) (Najafi et al., 2016) . Currently, Ma5, H120, and attenuated 793/B like IBV-based vaccination strategies are applied for IB control on poultry farms in Iran. Despite the use of vaccines, IBV still occurs in Iranian poultry farms. Vaccine failures are often associated with the emergence of antigenic variants that differ from the vaccine viruses (Huang and Wang, 2006) , thus, continuous surveillance programs are necessary to control IB. The present study was conducted to evaluate the frequency of occurrence of different IBV genotypes circulating in Iranian poultry farms in the period 2015-2017.

Trachea samples (10 tracheas from each flock) from 278 broiler flocks were collected from eight provinces of Iran from March 2015 to February 2017. The investigated broiler chicks showed signs of respiratory problems. Samples from 35 flocks were obtained from each of the 7 provinces of Golestan, Kerman, Ardabil, Isfahan, Ghazvin, Kurdestan, Razavi Khorasan, and 33 flocks from Khuzestan province. All flocks received the IBV vaccines (Massachusetts + 793/B type vaccine). Geographical locations of these provinces are shown in Fig. 1 .

Ten tracheas from each flock (n = 278) were pooled together and homogenized. RNA was extracted from the pooled samples using RNA easy mini kit (Qiagen) as recommended by the supplier. The cDNA was synthesized using a RevertAid First Strand cDNA synthesis Kit (Thermo Scientific) (Seger et al., 2016) .

Real-time PCR assay was done in order to amplify a conserved sequence within the 5′-untranslated region (UTR) of the IBV genome described by Callison et al. (Callison et al., 2006) . The amplification reaction was performed using a nested PCR assay targeting the hypervariable region 3 of spike gene (393 bp) sequences on cDNA of all positive samples in Real-time PCR step. SX1 (5′-CACCTAGAGGTTT-GYTWGCATG-3′) and SX2 (5′-TCCACCTCTATAAACACCYTTAC-3′) primer set were used in first reaction while SX3 (5′-TAATACTGGYAA-TTTTTCAGATGG-3′) and SX4 (5′ AATACAGATTGCTTACAACCACC-3′) primers were employed in the second round (Worthington et al., 2008) .

The cycling conditions for both first and second rounds of PCRs were an initial denaturation step at 94°C for 2 min and 35 cycles of denaturation at 94°C for 15 s, annealing at 58°C for 30 s, and polymerization at 72°C for 30 s. We sent all positive samples for sequencing. The Accu-Prep ® PCR purification Kit (Bioneer Co., Korea) was used for the purification of the PCR products. Sequencing was performed with the primers (both directions) (Bioneer Co., Korea). Chromatograms were evaluated with CromasPro (CromasPro Version 1.5). After sequences editing, we did NCBI BLAST on the results and did primary identifications. In the following, we draw the phylogenetic tree. Because of the high number of sequences in tree (Data has not been shown), the identical sequences have been removed in each group and finally selected 51 sequences for drawing the final tree and submission in gene bank. Phylogenetic tree for the S1 glycoprotein was generated using the neighbor-joining method with 1000 bootstrap replicates to assign confidence levels to branches using MEGA 7 software package by Kimura-2-parameter model (Kumar et al., 2016) . Sequences reported in this paper were added to the GenBank database under accession number (MF322802-MF322852).

From 2015 to 2017, 278 flocks were examined; of these, 145 (52.16%) were found to be IBV-positive. Based on phylogenetic analysis, the strains were clustered in four distinct groups (Fig. 2) . The distributions of different IBV genotypes expressed as a percentage of total IBV detected is shown in Table 1 . The most frequently occurring type was IS-1494-like IBV with 70.34% frequency of occurrence rate; 793/B was the second most common type with 19.31% frequency of occurrence. The third and fourth most frequently detected IBV types were Massachusetts and QX with rates of 7.58% and 2.75% respectively (Fig. 3) . QX IBV was not detected in Isfahan, Kurdistan, Golestan, and Khuzestan provinces. The highest frequently occurring of IBV was detected in Isfahan province (68.57%) and the lowest in Khuzestan province (39.39%). Sequence similarity among IS-1494-like IBV strains of this study varied between 98.86% to 100% and they had more than Table 2 ).

The avian infectious bronchitis virus (IBV) first isolated in 1937 can be destructive and cause severe economic losses in the poultry industry worldwide (De Wit et al., 2011) . Because of mutations or recombination events or both occurring during virus replication, dozens of serotypes and genotypes of IBV have been detected around the world (Cook et al., 2012) . The first isolation of IBV in Iran was reported by Aghakhan et al. (1994) . Genotyping of IBV strains isolated in Iran were classified into seven distinct phylogenetic groups (Massachusetts, 793/B, IS-1494, IS-720, QX, IR-1, and IR-2) (Najafi et al., 2016) .

In the present study, IS-1494-like IBV was the most dominant type, followed by the 793/B type. Massachusetts and QX type IBVs were the third and fourth most frequently detected types, respectively. In comparison to our last survey (2014-2015) and other reports from Iran the IR1, IS-720, and IR-2 IB genotypes have disappeared (Fig. 6) (Najafi et al., 2016) .

Although IS-1494-like IBV was first reported in Israel in 2006, Egypt/Beni-Suef/01 was isolated in Egypt, with 99% similarity to IS/ 1494/06 of the Israeli variant 2 isolates, in 2001 (Gelb et al., 2005) . Since then IS-1494-like IBV type was detected in Jordan, Turkey and other Middle East countries . In Iran, IS-1494like IBV was detected in 2010 for the first time and was reported as the second most frequently occurring IBV from 2010 to 2014 (Hosseini et al., 2015) , thereafter it turned out to be the most dominant type in the molecular study done from 2014 to 2015 (Najafi et al., 2016) . Based on the data, not only is IS-1494-like IBV type still the predominant IBV type, its appearance has remarkably increased up to two times more than that reported in the previous study (Najafi et al., 2016) . Since the IS-1494-like IBV type vaccines are not used in Iran and because current vaccines and vaccination strategy do not protect completely against IS-1494-like, it is highly likely to become more widespread (Habibi et al., 2017) . 793/B IBV type was first found in 1998 in Iran (Vasfi-Marandi and Bozorgmehrifard, 2001) . According to the last two epidemiological data, the frequency of occurrence was 8.4% from 2010 to 2014, rising up to 21% from 2014 to 2015 (Najafi et al., 2016) . According to our results, its presence is still high (19.33%) despite vaccination. As most 793/B IBVs detected in this study have the same nucleotide sequences with those of vaccine, the high frequency of occurrence probably reflects re-isolation of 793/B like IB vaccine strains. Like previous data, some 793/B type viruses were detected which differ from the vaccine Fig. 2 . Phylogenetic tree based on a partial sequence of the S1 gene, showing the relationship between the Iranian strains and other IBV strains. The neighbor-joining method was used with the Kimura2-parameter substitution model and 1000 bootstrap replicates to assign confidence level to the branches of the phylogenetic tree. Strains detected in the current study are indicated by black circles while Iranian previously identified strains marked by white squares. The vertical lines are for spacing branches and labels. The scale bar represents the distance unit between sequence pairs. 

Percent identity of partial nucleotide sequences of the S1 glycoprotein genes of some Iranian IBVs to those of IBV reference strains.

A. Modiri Hamadan et al. Virus Research 240 (2017) 101-106 viruses. QX type IBV was first isolated by Bozorgmehri-Fard et al. in Iran, the incidence was 9.6% from 2010 to 2014 and 10% from 2014 to 2015 (Najafi et al., 2016) .The present study shows a reduction in the frequency of occurrence rate of QX IBV, which may reflect the efficacy of combined Massachusetts and 793/B type IB vaccines. Epidemiological analysis of IBV in Italy and Spain indicated a causal relationship between Massachusetts plus 793/B vaccination and a decreased prevalence of QX (Franzo et al., 2016) . Kahya et al. detected IS-1494 -06 like IBV strain in turkey, 2011 (Kahya et al., 2013 . S1 analysis of Turkish IBVs (2014) (2015) showed 17.2% of the sequences were closely related to MA5 vaccine strain, 49.5% to H120, and 48.9% of those were similar to the IS-1494-like IBV (Yilmaz et al., 2016 (Mahmood et al., 2011) . In early 2011, Iraqi QX (IBV/CH/Kurdistan-Iraq/7266-10/ 2011) was isolated (Amin et al., 2012) .The percentage of IS-1494-like, 793/B, QX, and DY12-2 genotypes in five Iraqi governorates was 46.87%, 40.62%, 9.37% and 3.12% in 2014, respectively (Seger et al., 2016) . The presence of antibody against Massachusetts (M41), Arkansas, Connecticut, D-274, D-1466, and 4/91 IBV strains have been detected (Ahmed et al., 2007; Munir et al., 2012) . Of 39 isolates obtained from 243 Oman backyard poultry in 2012, 66.67% of samples showed the greatest homology to genotype 793/B, 5.12% related closely to M41, 5.12% clustered with D274 and 2.56% belonged to IS-1494-like (Al-Shekaili et al., 2015) . In another study from 2009 to 2014, the prevalence of 793/B, IS-1494-like, and Massachusetts genotypes were reported to be 63.7%, 21.2%, and 3%, respectively in Oman Fig. 4 . Percent identity of partial S1 gene sequences of some IBVs from the current study to that of the H120 strain. Fig. 5 . Percent identity of partial S1 gene sequences of some IBVs from the current study to that of the 4/ 91 vaccine strain. . In summary, IS-1494-like and 793/B are the most prevalent IBV genotypes circulating in our region such as Oman, Iraq, and Turkey. Due to intensive trading and uncontrolled movement of inhabitants and animals across borders, approximately the same distribution of IBV genotypes exists. The updated data from continuing molecular surveillance of IBV in Iran completes the puzzle of circulating IBV genotypes in the region. The study helps us to monitor the IBV strains to see whether new variants have emerged and we can evaluate the efficacy of our national IB vaccination program. Further, it is necessary to develop a vaccine in order to control the rising spread of the IS-1494-like IBV type.

All authors declare that they have no conflict of interest.

Bioinformatics and evolutionary insight on the spike glycoprotein gene of QX-like and Massachusetts strains of infectious bronchitis virus

Studies on avian viral infections in Iran

Detection and seroprevalence of infectious bronchitis virus strains in commercial poultry in Pakistan

Molecular detection of infectious bronchitis and avian metapneumoviruses in Oman backyard poultry

Circulation of QX-like infectious bronchitis virus in the Middle East

Development and evaluation of a real-time Taqman RT-PCR assay for the detection of infectious bronchitis virus from infected chickens

The long view: 40 years of infectious bronchitis research

Infectious bronchitis virus variants: a review of the history, current situation and control measures

Effect of different vaccination strategies on IBV QX population dynamics and clinical outbreaks

Genotypes of infectious bronchitis viruses circulating in the Middle East between

S1 gene characteristics and efficacy of vaccination against infectious bronchitis virus field isolates from the United States and Israel

Combination of H120 and 1/96 avian infectious bronchitis virus vaccine strains protect chickens against challenge with IS/1494/06 (variant 2)-like infectious bronchitis virus

Epidemiology of avian infectious bronchitis virus genotypes in Iran

Development of attenuated vaccines from Taiwanese infectious bronchitis virus strains

Review of infectious bronchitis virus around the world

Presence of IS/1494/06 genotype-related infectious bronchitis virus in breeder and broiler flocks in Turkey

MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets

Complete genome sequence analysis of a predominant infectious bronchitis virus (IBV) strain in China

Isolation and molecular characterization of Sul/01/09 avian infectious bronchitis virus, indicates the emergence of a new genotype in the Middle East

Quantification of antibodies against poultry haemagglutinating viruses by haemagglutination inhibition test in Lahore

Molecular characterization of infectious bronchitis viruses isolated from broiler chicken farms in Iran

Genotyping of infectious bronchitis viruses from broiler farms in Iraq during

A survey of the prevalence of infectious bronchitis virus type 4/91 in Iran

Isolation and identification of infectious bronchitis viruses in chickens between 1997 and 2000 in Iran

A reverse transcriptase-polymerase chain reaction survey of infectious bronchitis virus genotypes in Western Europe from

Phylogeny and S1 gene variation of infectious bronchitis virus detected in broilers and layers in Turkey

We would like to thank Ghalyanchi lab experts (