Two genotypes of infectious bronchitis virus are responsible for serological variation in KwaZulu-Natal poultry flocks prior to 2012

Copyright: © 2014. The Authors. Licensee: AOSIS OpenJournals. This work is licensed under the Creative Commons Attribution License. This study describes the isolation, serotyping and genotyping of 54 infectious bronchitis virus (IBV) cases predominantly in KwaZulu-Natal and compared to several isolates from other South African provinces between 2011 and 2012 and several historic isolates. The results indicate the division of isolates into two different genotypes of IBV within the province, Massachusetts (Mass)-like and QX-like. The IBV Mass-like genotype was the most prevalent and was detected in 79% of the full spike protein S1 gene sequences. Variation up to 22.3% was detected within local Mass-type strains, supporting the hypothesis that multiple IBV serotypes may co-circulate in the same region simultaneously. Additionally, more conservation was observed amongst Mass serotypes versus QX-like serotypes, implying that vaccine use can influence the variability within the IBV population; this is deduced from the fact that the only live vaccine registered for use in South Africa at the time of the study was of Mass origin and no QX-like vaccines were available for use. This study offers the first published consolidation of IBV isolates from an area of South Africa and identifies variation within the IBV population of the broiler flock within the study area over a 2-year period.


Introduction
Infectious bronchitis virus (IBV) is a highly contagious coronavirus of poultry that causes disease characterised by upper respiratory or urogenital lesions.IBV is a listed disease according to the World Organisation for Animal Health (OIE) and can result in large losses due to mortality or loss of production.The disease has a worldwide distribution.The virus is shed both through the upper respiratory system and in faeces and may be detected within the bird's gastrointestinal system for several weeks or months (Cavanaugh & Naqi 1997).Clinical signs and lesions include respiratory symptoms, effects on egg production and egg shell quality and kidney pathology.Mortality rates vary between 14% and 82% depending on the IBV serotype and coinfections with other pathogens, particularly Escherichia coli and Mycoplasma spp.(Cavanaugh & Naqi 1997;Cook 1983;Cook, Huggins & Ellis 1991;Dhinakar Raj & Jones 1997).
Coronaviruses are enveloped single-stranded positive-sense RNA viruses.The IBV genome is 27.569 kbp in length and contains four coronavirus structural protein genes, spike (S), envelope (E), membrane (M) and nucleocapsid (N), plus non-structural protein genes 3 and 5 and its viral replicase gene (Rep).The order of these genes is 5′ Rep-S-3-E-M-5-N 3′ (Boursnell et al. 1987;Spaan, Cavanaugh & Horzinek 1988).IBV S protein undergoes post-translational cleavage to S1 and S2 subproteins.The S1 protein is responsible for attachment to the host cell membrane and also has epitopes to which neutralising and haemagglutination-inhibiting antibodies bind (Kant et al. 1992;Koch et al. 1990).Two hypervariable regions (HVRs) have been demonstrated within the 50-159 residue regions of the S1 gene.HVR1 lies within region 38-51 whilst HVR 2 is located between amino acids 99 and 115 (Cavanaugh, Davis & Mockett 1988;Cavanaugh et al. 1992;Wang, Xu & Collisson 1997).A third HVR between amino acids 274 and 387 has also been described (Lee & Jackwood 2001).The S1, S2 and M proteins have all been shown to induce cell-mediated immunity with the N gene proposed as carrying the predominant T cell epitope (Ignjatovic & Galli 1994).The high degree of variability within the S1 gene has been attributed to the lack of proofreading by the viral RNA-dependent RNA polymerase and to recombination events during replication (Toro, Van Santen & Jackwood 2012).
The emergence of variation within the IBV population has been widely described since 1951, when Jungherr, Chomiak and Luginbuhl (1956) demonstrated that the Connecticut isolate neither cross-neutralised nor cross-protected against the traditional H120 Massachusetts (Mass) isolate of the 1940s.The first published reports of IBV from South Africa were by Morley and Thompson in 1984.This was confirmed as an unknown variant by Cook et al. in 1999. These viruses Page 1 of 10 Scan this QR code with your smart phone or mobile device to read online.

Read online:
Two genotypes of infectious bronchitis virus are responsible for serological variation in KwaZulu-Natal poultry flocks prior to 2012 were associated with swollen-head syndrome observed in infected flocks throughout southern Africa.Novel IBV has been described in Nigeria (De Wit, Cook & Van der Heijden 2011) as well as Zimbabwe, where a QX-like IBV strain has been reported (Jackwood 2012).
The South African poultry industry is the single largest agricultural subsector, accounting for an annual turnover of R31.7 billion per annum (South African Poultry Association 2012) with flock sizes ranging from several hundred to large complexes of 600 000 birds.Perdue and Seal (2000) listed the estimated losses per flock from an IBV infection as 10% to 20% of market value.This can devastate a poultry farm, where margins are relatively small.The economic impact of this disease is thus of great importance.Currently available research findings do not describe variation over a large collection of South African IBV isolates and thus do not adequately describe the degree of variation within the South African broiler flock.The core question revolves around the extent of variation present in the South African IBV population.Objectives of this study included identifying variation within IBV isolates, describing the correlation between serotyping and genotyping, as a method of variation identification and describing variation within the KwaZulu-Natal IBV isolates over time.

Sample collection
Over 100 IBV isolates were collected during the 18-month period; 46 of these were selected for the present study based on clinical history and pathology caused by the strain, geographic distribution and serum neutralisation test results.In addition, eight prior isolates were also included in the study.These isolates included six Rainbow Veterinary Laboratory IBV reference strains and two historical isolates (Table 1).Seventy-five percent of isolates selected for the study originated from KwaZulu-Natal (KZN), 11% from the North West province, 3.5% each from the Eastern Cape and Gauteng and 6% from the Western Cape.The weighting towards KZN isolates demonstrated the higher isolation rate from this province in addition to the higher suspected infection rate of this subpopulation within the study population.

Virus isolation
Samples for isolation were collected from postmortem samples on suspicion of IBV infection.Organs collected included either tracheal samples or pooled tracheal and kidney samples.Samples were collected from various postmortem facilities around South Africa but all virus isolations were performed at the KZN laboratory.
The standard virus isolation procedures are based on the published laboratory manual for the isolation, identification and characterisation of avian pathogens in 9 to 10-day-old embryonated chicken eggs via the allantoic sac inoculation route (Gelb 1989).Dead embryos were opened after candling; diagnosis of the cause of death was conducted and in the case of evidence that the embryo was infected with IBV, the allantoic fluid was collected.Typical signs of IBV infection in the embryo, namely stunting, urate deposits in the mesonephros and clubbing of the down feathers, were noted.In the cases where no pathological signs were seen in the embryo, the allantoic fluid was still collected and reinoculated for an additional passage.A total of three passages were performed until the sample was confirmed negative.Inoculations with supporting IBV clinical signs were confirmed for the presence of IBV through a chloroform inactivation test as well as diagnostic polymerase chain reaction (PCR) assay.The chloroform inactivation test differentiates between enveloped and non-enveloped viruses, and in combination with embryo pathology and haemagglutination tests a diagnosis of IBV can be made.PCR-confirmed IBV-positive allantoic fluid isolates were then stored at -70 °C until further processing.

Virus neutralisation tests
Viral antigen-containing allantoic fluid was treated in serial dilution with antisera collected against each reference isolate (Gelb 1989).These antisera were produced through the hyperimmunisation of specific-pathogen-free leghorn chickens.

RNA extraction and reverse transcription polymerase chain reaction
RNA was extracted using Trisure reagent (Bioline) according to the procedure recommended by the manufacturer.A twostep PCR procedure was followed: viral cDNA was generated by combining 5 µg of extracted PCR products were electrophoretically separated on a 1% agarose gel stained with ethidium bromide.A KAPA™ universal ladder was included to identify the 1700 bp S1 gene amplicon.Bands of the correct size were excised and DNA purified using QiaQuick PCR purification kit (Qiagen) following the recommended procedure.
Purified 1.7 kb PCR amplicons were submitted to Inqaba Biotech (Pretoria) for Sanger sequencing.Results were edited using Chromas Lite and then aligned using BioEdit version 7.1.11(Hall 1999).The Basic Local Alignment Search Tool (BLAST) server was used to retrieve similar sequences selected from Genbank for comparative analysis, together with important reference viruses.Pairwise similarity was calculated using BioEdit (Hall 1999) and phylogenetic trees were constructed in MEGA 5.2 (Tamura et al. 2011).Full S1 gene sequences were deposited in Genbank under the accession numbers KJ200273-KJ200291.Partial sequences are available from the corresponding author on request.
Phylogenetic relationships were inferred using the neighbour-joining tree inference method with a Jukes-Cantor model of sequence evolution using MEGA 5.2 (Tamura et al. 2011).One thousand bootstrap samplings were performed to assign confidence values to branching orders.

Ethical considerations
This study was approved by the University of Pretoria's Animal Ethics Committee under reference number V006/12.

Results
A total of 54 isolates confirmed as IBV positive by serum neutralisation and diagnostic PCR were analysed by reverse transcription polymerase chain reaction (RT-PCR).The isolates were named in accordance with the system proposed by Cavanaugh (2001) (Table 1).
Forty-one samples produced a band of approximately 1700 base pairs (data not shown), of which full S1 gene sequences could be assembled for 19 isolates, with the remainder generating either partial forward or reverse sequences.Despite numerous attempts on the service provider's part to troubleshoot sequencing reactions, they were unable to generate the complementary forward or reverse reactions for 30 of the amplicons submitted.The reasons for this are unknown, since the exact same forward and reverse primers were used for sequencing as for RT-PCR, and numerous re-submissions of fresh amplified RT-PCR product in sufficient molar quantities were made.Phylogenetic analysis of the full S1 gene sequences are presented in Figure 1.Partial sequence trees are not shown, but were used to genotypically classify the isolates (results are presented in Table 1).Isolates were classified as either QX-like (n = 4) or Mass-type (n = 16).
In this QX-like group, residues 77-85 contained the motif KDVYNQSVA compared to the Mass grouping with a motif of HGGRVVNA.A large insertion was also detected between residues 122 and 129 with a motif of SSGSGSCP, which is not present in any of the Mass-type isolates.Two point mutations within CK/ZA/3665/11 resulted in the replacement of asparagine with histidine at position 24 and then aspartic acid with tyrosine at position 248.Variation was observed between the two 2011 and the 2012 isolates within the QX-like group; the protein sequence at position 18-22 of the 2012 isolates was FALWK and CALCS as compared to the 2011 isolates at the same position.The most variable portion of the amino acid sequence for this group occurred at residue 400, where two of four isolates had consensus and variation between asparagine, histidine and tyrosine was seen.
The remaining isolates were closely related to reference sequence H120 (Mass serotype).Single-point mutations were observed at various points across the isolates.No variation in       Pairwise similarities were calculated between the study isolates (Table 2) and selected serotype sequences obtained from Genbank (Table 3).Pairwise comparisons against reference sequences support the description of two distinct groups of isolates.QX-like isolates shared 97.2% -98.4% nucleotide sequence identities, corresponding to 97.6% -98.8% identity at the amino acid level.These QX-like isolates only shared 77.7% -79% nucleotide sequence and 85.9% -88.6% amino acid sequence identities with the Masslike isolates.Sequence variation within the Mass-like group was less, at 98.9% -100% and 99% -100% nucleotide and amino acid similarities respectively (Table 2).Four isolates (CK/ZA/0890/81, CK/ZA/6689/12, CK/ZA/5315/11 and CK/ZA/2287/11) were very closely related.These sequences differed by 0.3% or four point mutations at positions 294, 302, 374 and 483.
The remainder of the Mass-type South African isolates shared high degrees of similarity with two point mutations common to all isolates when compared to the reference strain.The replacement of thymine with cysteine at residue 11 and thymine with adenine at residue 138 was found amongst all Mass isolates.The highest frequency of variation in this case was found between base pairs 706 and 818, with 10-point mutations found in this region.

Discussion
Phylogenetic analysis of the S1 genes revealed two distinct groupings of IBV within the study isolates.The larger group contained 15 of the initial 19 full sequences and shared 99.0% -99.7% amino acid sequence similarity with the H120 serotype.Within this group there was 98.8% -100% similarity at the nucleotide sequence level.Reported S1 gene mutation rates under vaccine selective pressure are 2.5% (Lee & Jackwood 2001) and reports suggest that in the absence of a selective vaccine, mutation rates are twice as high than when under selective pressure (McKinley et al. 2011).
When compared against the year or area of isolation it was found that the highly similar Mass grouping originated from KZN between April 2011 and September 2011 and then again in October 2012.One isolate from the North West province from July 2011 and two from the Western Cape collected in October 2011 also fell within this group.Strong correlations were observed between the isolates that originated outside KZN.Despite the geographic areas being far apart there was movement of eggs and embryos between provinces during the study period.The variability of 0.4% observed between CK/ZA/6743a/11 and CK/ZA/6743b/11 is interesting, since these samples were collected from the same farm on the same day but originated from separate chicken houses or units.The existence of an insertion at 673 and point mutations at 4, 6, 296 and 879 demonstrate the existence of variability and selection of subpopulations on a single farm.The selective criteria for these subpopulations are unknown, but variables such as maternally derived antibody levels and timing of vaccination could apply greater selective pressure for the emergence of variants of the epitopes found on vaccinederived viruses (De Herdt et al. 2001;Meulemans et al. 2001;Mondal & Naqi 2001).An alternative hypothesis is that the individual house conditions with regard to air quality and the degree of damage to the chicken's tracheal epithelium and cilia affect the selection of S1 protein subpopulations.

TABLE 2:
Percentage similarity of study isolates with selected reference strains.The use of especially the live day-old H120 vaccine and the close similarity of Mass-type isolates to H120 must not be overlooked as the most likely influence on the prevalence of the H120 group.This implies that the most likely contributors to IBV serotypes within these flocks are Mass-type live vaccines.However, there is the possibility that either through genetic drift or quasi-species selection, field-adapted virulent Mass-type IBV was able to proliferate in suitable hosts.The basis of this statement is that all the samples collected were from clinical cases of respiratory disease and the presence of suggestive pathology on postmortem examination.Recombination and reversion to virulence within Mass-type IBV vaccines have been described and field experience over the years of use has found that poor vaccination technique, concurrent viral or bacterial infections as well as on-farm management conditions can lead to IBV pathology in light of vaccine usage (Cavanaugh, Davis & Cook 1992).

Isolate
Variability within the QX-like genotype was greater at 97.7% -98.5% than that seen in the Mass grouping at 99.0% -99.7%.A temporal trend in genetic drift was observed.This observation is in contrast to evolutionary patterns within the Mass group and supports the existence of genetic drift within this IBV subpopulation.The greatest similarity was noted when comparing it to the QX-like strains of European origin (97%) as opposed to the Asian QX-like strains (94%), which provides circumstantial evidence of the spread of QX-like IBV from Asia into Europe and then downward into Africa.
The results of the KZN isolates would support the high host adaptability of the QX-like variant and identify it as the leading 'variant' within KZN poultry flocks.Shaw, Britton and Cavanaugh (1996) investigated Belgian nephrotropic IBV strains and reported that the S1 gene differed by 21% -45% from vaccine strains whilst the S2 gene only had 10% -12% differentiation.The 793/B strain was shown to display between 21% and 25% variation amongst field isolates, whilst still retaining the same serotype (Adzhar & Gough 1997).The CK/CH/LDL/971 strain of IBV was shown to have < 80% homology between the vaccine and pathogenic strain as well as between itself and other IBV types (Liu et al. 2009).Host-driven selection has been postulated and can occur as demonstrated during egg passage of several 793/B isolates; although this selection did not alter the pathogenicity, it does illustrate the quasispecies theorem regarding IBV populations (Cavanaugh et al. 2005).Jackwood, Hilt and Callison (2003) utilised real-time PCR techniques and melting curves to visualise the existence of quasi-species within the IBV Beaudette strain.Re-isolated virus has been demonstrated to contain subpopulations, again supporting the quasi-species theorem (McKinley, Hilt & Jackwood 2008).
The correlation between QX-like serotypes and the actual sequence was inconsistent, as only 4 of the 11 serotype QXlike isolates were classified within the QX-like genotype.
All seven remaining Mass serotypes genotyped into the same group.By comparison, when examining partial S1 gene sequence data, none of the QX-like partial sequences were serotyped as QX-like.Seven serotype QX-like isolates genotyped as Mass, whilst the two isolates that did not fall into a known serotype at the time were split between the Mass and QX-like genotype.The above results cast doubt on the correlation between serotype and genotype results.It has been described that small sequence variations may alter the virus neutralising epitope and thereby the serotype of the virus (Cavanaugh et al. 1992).An alternative hypothesis could be that coinfection with multiple strains of IBV is responsible for clinical disease, and between the culture or serum neutralisation test and the subsequent successive passages and RT-PCR selection for host-specific strains there is selection for an alternate dominant strain (Liu et al. 2009;Xu et al. 2007).Comparing these protein sequences to the isolates that both genotyped and serotyped as Massachusetts, there is no obvious genetic marker that serotyped as QX-like in contrast to the group that serotyped as Mass-like.

Conclusion
This study describes the isolation, identification and serotyping of a representative selection of IBV, predominantly in KwaZulu-Natal, between 2011 and 2012.The Massachusetts (H120) serotype was the most prevalent serotype within the KZN poultry flock, according to the samples analysed in this study.The evidence suggested that the most commonly isolated Mass-types are primarily derived from variation or mutation of vaccine strains administered to flocks.The QXlike serotype is the only non-vaccine serotype identified in the sample pool and had a greater degree of variation within the group than the Mass-type isolates.The presumption therefore is that greater genotypic variation is seen amongst IBV serotypes that are not exposed to vaccine selective pressure.

FIGURE 2 :
FIGURE 2: Multiple amino acid sequence alignment for full S1 proteins.Hypervariable regions are indicated in grey.The first sequence is used as a reference and identities are plotted as '.'.

Figure 2
Figure 2 continues on next page 

CK
similarity; AAS, Amino Acid similarity.†, comparing the same two isolates and thus similarity is 100%.

Table 2
continues on next column 

TABLE 3 :
Percentage similarity of study isolates with selected reference strains.Denotes highest comparitive similarity in study, ‡, Nucleotide sequence identity of the 1460 bp sequence; §, 427 amino acid translation similarity.