Detection and molecular characterization of infectious bronchitis virus isolated from recent outbreaks in broiler flocks in Thailand

Tawatchai Pohuang; Niwat Chansiripornchai; Achara Tawatsin; Jiroj Sasipreeyajan

doi:10.4142/jvs.2009.10.3.219

Abstract

Thirteen field isolates of infectious bronchitis virus (IBV) were isolated from broiler flocks in Thailand between January and June 2008. The 878-bp of the S1 gene covering a hypervariable region was amplified and sequenced. Phylogenetic analysis based on that region revealed that these viruses were separated into two groups (I and II). IBV isolates in group I were not related to other IBV strains published in the GenBank database. Group 1 nucleotide sequence identities were less than 85% and amino acid sequence identities less than 84% in common with IBVs published in the GenBank database. This group likely represents the strains indigenous to Thailand. The isolates in group II showed a close relationship with Chinese IBVs. They had nucleotide sequence identities of 97-98% and amino acid sequence identities 96-98% in common with Chinese IBVs (strain A2, SH and QXIBV). This finding indicated that the recent Thai IBVs evolved separately and at least two groups of viruses are circulating in Thailand.

Introduction

Infectious bronchitis (IB), caused by the IB virus (IBV), is a highly contagious disease and results in a significant economic loss to the commercial chicken industry. The disease frequently causes respiratory signs including gasping, coughing, sneezing, tracheal rales, and nasal discharge. In layer fowls, respiratory distress and a decrease in egg production have been reported [10]. In addition, some strains have been associated with kidney lesions [18,19,22].

IBV, the causative agent of IB, belongs to the genus Coronavirus in the family Coronaviridae [6]. It is an enveloped virus and has a positive-sense, single-stranded, RNA genome, approximately 27 kb in length. The virion has three major structural proteins namely the nucleocapsid (N) protein, the membrane (M) protein, and the spike (S) glycoprotein. The S glycoprotein is post-translationally cleaved into the S1 and S2 subunits [6]. The S1 subunit, located on the outside of the virion, is responsible for the fusion between the virus envelope and the host cell membrane. Moreover, it is responsible for neutralizing serotype-specific antibodies in chickens [4]. The S1 subunit demonstrates more sequence variability than S2 [14]. Neutralizing and serotype-specific epitopes are associated with the defined hypervariable region (HVR) in the S1 subunit; therefore, the molecular characterization of IBV is based on analysis of the S1 gene [12].

During 1953-1954, the first IBV outbreak was reported in Thailand [7]. Although many IB vaccine strains such as Connecticut, H120, Ma5, M41 and Armidale A3, have been used in Thailand for many years, IBV outbreaks have been ongoing [1]. Moreover, the relationships between recent Thai IBV isolates and foreign IBV isolates are not known.

The objective of this study was to characterize IBV field isolates in the recent outbreak of the disease in Thailand by analyzing the S1 gene and compared them with those that have been published previously.

Materials and Methods

Viruses

Between January and June, 2008, thirteen poultry farms in the eastern part of Thailand had an outbreak of a mild-to-moderate respiratory disease (Table 1). All flocks had been vaccinated against IB with commercial live attenuated H120. Chickens showed respiratory symptoms including gasping, coughing, sneezing, and tracheal rales. Sick chickens were selected and sent to the Department of Veterinary Medicine, Faculty of Veterinary Science, Chulalongkorn University. Necropsy was performed and gross lesions were evaluated. Gross lesions showed mild- to moderated tracheitis and non-purulent airsacculitis. No gross lesion were found in the kidneys. The trachea and lung samples were taken as pools of chickens from the same farm. The samples were prepared as 10% w/v suspensions in phosphate-buffered saline (pH 7.4) and clarified at 1,800 × g for 10 min; the supernatants were then collected for analysis.

RNA extraction

Viral RNA was extracted by using Viral Nucleic Acid Extraction Kit (Real Biotech, Taiwan) following the manufacturer's instructions directly from the supernatants of 10% w/v sample suspensions and from the allantoic fluid of embryonated chicken eggs used for virus isolation.

Virus screening with nested RT-PCR

Viral RNA, extracted directly from the supernatants of 10% w/v sample suspensions, was screened for the presence of IBV by using a nested RT-PCR. The reaction was performed with AccessQuick RT-PCR System (Promega, USA). The first amplification reaction was carried out with one-step RT-PCR using the primer sets of FOR1 (5'-CTT TTG TTT GCA CTA TGT AG-3') and RE3 (5'-TAA TAA CCA CTC TGA GCT GT-3'). The second amplification reaction was carried out using the primer sets of FOR2 (5'-CAG TGT TTG TCA CAC ATT GT-3') and RE2 (5'-CCA TCT GAA AAA TTG CCA GT-3'). Amplification products were analyzed in 1.5% agarose gel. The predicted size of nested RT-PCR product was about 400-bp.

Virus isolation and propagation

For virus isolation, the supernatants of IBV-positive samples determined by RT-PCR were inoculated into 10-day-old embryonated chicken eggs. For each sample to be examined, five embryonated chicken eggs were used. The eggs were inoculated with 0.2 mL of the sample into the allantoic cavity. The inoculated eggs were incubated at 37℃ and candled daily. Allantoic fluids were harvested at 96 h postinoculation. A further blind serial passage was performed in a similar way. All of the allantoic fluids were harvested and stored at -70℃.

RT-PCR amplification for sequencing

The allantoic fluids from the second passage of each sample positive for virus screening was submitted to another RT-PCR for amplification of a segment of 878-bp of the S1 gene coding region using a primer combination of FOR1 and RE3. The amplification reaction was carried out with one-step RT-PCR and the reaction conditions were the same as the first amplification of nested RT-PCR described above.

Product purification and sequencing

The RT-PCR products were cut from the gel and purified using the Wizard SV Gel and PCR Clean-Up system (Promega, USA) according to the manufacturer's protocol. Purified RT-PCR products were sequenced in a forward direction using primer FOR1 and in a reverse direction using primer RE3. Sequencing reactions were performed with the ABI Prism BigDye Terminator Cycle Sequencing Ready Reaction Kit (Applied Biosystems, USA) as described by the manufacturer. Sequencing reactions were run on an ABI Prism 310 Genetic Analyzer.

Sequences and phylogenetic analysis

To identify the Thai IBV isolates, sequences of the S1 gene of the Thai IBV isolates were compared with published IBV sequences deposited in the GenBank database using a BLAST search via the National Center of Biotechnology Information (USA). Sequence identities by BLAST analysis were included in alignment and phylogenetic construction. A phylogenetic tree of the nucleotide sequences was constructed using MEGA version 3.1 [13]. The S1 gene sequences of the thirteen IBV isolates were submitted to the GenBank database (Table 1). The other S1 gene sequences from the GenBank database were used for comparison or phylogenetic analysis in this study included M41 (AY561711), Ma5 (AY561713), H120 (M21970), IBN (AAW83034), W93 (AY427818), Connecticut 46 (L18990), Florida 18288 (AF027512), JMK (L14070), Spain/99/319 (DQ064810), Spain/00/337 (DQ064813), J2 (AF286303), BJQ (DQ070839), QXIBV (AF193423), LC2 (DQ480154), A2 (AY043312), SH (DQ480156), K069-01 (AY257061), 4/91 (AF093794), UK2/91 (Z83976), Ark DPI (AF006624), Australian T (AY775779), N1/62 (AIU29522), Armidale (DQ490205), GA/7994/99 (AF338717), GA/8077/99 (AF338718), DLD (EU589323), THA001 (DQ449628).

Results

Virus screening and isolation

For virus screening, pooled trachea and lung samples from each flocks suspected of IBV infection were determined to be positive for IBV by screening with nested RT-PCR. A 400-bp fragment of the S1 gene was amplified in all 13 samples tested (data not shown). The allantoic fluid from the second passage of each sample screened positive for the virus was also determined to be positive with RT-PCR amplification and a segment of 878-bp of the S1 gene was obtained (data not shown).

Phylogenetic analysis

To assess the genetic relationship among the IBV isolates, a phylogenetic tree was constructed from the nucleotide sequences of S1 genes. The results are shown in Fig.1. The thirteen IBV isolates were separated into two distinct groups. Group I consisted of five isolates including THA20151, THA40151, THA50151, THA60151, and THA90151. The isolates in group I showed evolutionary distances from each other. Group II consisted of eight isolates including THA30151, THA 70151, THA 80151, THA100151, THA110351, THA120351, THA130551, and THA140551, which had a close relationship with Chinese IBV isolates (strain A2, SH and QXIBV).

Nucleotide and amino acid sequence comparison

The S1 gene of the thirteen IBV isolates was sequenced to characterize the isolates. The nucleotide and deduced amino acid sequences were determined and compared among each other and with other IBV strains published in the GenBank database. Group I Thai IBV isolates had nucleotide and amino acid sequence identities between 99~100% with each other. They had nucleotide sequence identities less than 85% and amino acid sequence identities less than 84% with IBVs published in the GenBank database. Group II Thai IBV isolates had nucleotide and amino acid sequence identities between 99~100% with each other. They had nucleotide sequence identities of 97~98% and amino acid sequence identities of 96~98% with Chinese IBVs (strain A2, SH and QXIBV).

Discussion

One of the major problems of IBV is the frequent emergence of new variants [21]. Different serotypes have been reported world wide and new variant serotypes continue to be recognized [9,18]. Thus, it is necessary and important to be able to diagnose these new serotypes. Furthermore, determining the type as well as field isolates is important to select an appropriate vaccine against IBV infection in the next flock. Therefore, several tests have been employed to identify the isolates into serotypes or genotypes [15]. Typing with RT-PCR and sequencing of the S1 gene is easier and faster than the more traditional virus neutralization methods, but it is difficult to design PCR primers that can be used for detect all of IBV isolates [16]. In this study, IBV isolates in group I could not be detected with the primer set and PCR method described by Gelb et al. [8]. To overcome this problem, the new primer sets were designed and the nested RT-PCR was developed in order to increase the sensitivity and specificity of the test used to detect IBV from the infected samples. Moreover, the region amplified by the outer primer covering the HVR of the S1 gene (nucleotide position 31~908) could be used for typing of IBV.

Genotyping of IBV on the basis of the S1 gene sequence, particularly the HVR region of the S1 gene, is the most common way to classify IBV isolates. It has been shown that the genetic typing based on HVR I (nucleotide position 114~325) of the S1 gene could represent the grouping method based on the whole S1 gene [20]. Furthermore, in a recent study using an approximately 450-bp region covering HVR I and HVR II for IBV typing, it was found that typing by this region correlates with virus neutralization results [17]. In the present study, an approximately 878-bp region of S1 gene covering HVR I and HVR II were amplified and used for typing the field isolates in Thailand.

Although many different strains of live attenuated and inactivated vaccines have been widely used to control IB, the outbreaks of the disease have continued to be a problem in Thailand [1]. In this study, thirteen IBV isolates from different commercial poultry farms in Thailand were analyzed by sequencing of the HVR in S1 gene. The molecular data indicated that the recent IBV isolated in Thailand evolved separately into two groups. The IBV isolates in group I were not related to other IBV strains published in the GenBank database. This genotype likely represents strains indigenous to Thailand. IBV isolates in group II were genetically related to Chinese strains.

The results from both S1 gene comparison and phylogenetic analysis showed that IBV isolates in group I had a distant relation to vaccine strains used in Thailand including Ma5, H120, M41, and Connecticut. New serotypes or variant strains can emerge as a result of only a few changes in the amino acid composition of S1 gene [5]. The rationale for these changes could be due to immunological pressure caused by the wide spread use of vaccines, recombination as a consequence of mixed infections, or the decrease of dominant serotypes as a result of vaccination, allowing other field strains to emerge [19].

Group II Thai isolates and Chinese strains shared more than 97% nucleotide identity and 96% amino acid identity; therefore, they were grouped into the same genotype. Interestingly, group II Thai isolates were closely related to Chinese QXIBV which appeared to become widespread in several countries in the world including the UK [11], Russia [3], and Italy [2]. Although no information is available currently for the introduction of QXIBV from China to other countries, it had been hypothesized that wild birds were the source of introduction based on the evidence that IBV may replicate in Anseriformes [3].

In summary, this finding indicated that the recent Thai IBVs evolved separately and at least two groups of viruses are cocirculating in Thailand. The distinctive dissimilarity between group I Thai isolates and the widely used IBV vaccine indicated that the antigenic drift is likely to occur among some Thai IBVs under the long-term immune pressure. Group II Thai isolates were closely related to Chinese strains, suggesting that the transfers of IBV infection occur among countries.