Analysis of a QX-like avian infectious bronchitis virus genome identified recombination in the region containing the ORF 5a, ORF 5b, and nucleocapsid protein gene sequences

Genome Sequence Variations of Infectious Bronchitis Virus Serotypes From Commercial Chickens in Mexico

New variants of infectious bronchitis viruses (IBVs; Coronaviridae) continuously emerge despite routine vaccinations. Here, we report genome sequence variations of IBVs identified by random non-targeted next generation sequencing (NGS) of vaccine and field samples collected on FTA cards from commercial flocks in Mexico in 2019–2021. Paired-ended sequencing libraries prepared from rRNA-depleted RNAs were sequenced using Illumina MiSeq. IBV RNA was detected in 60.07% (n = 167) of the analyzed samples, from which 33 complete genome sequences were de novo assembled. The genomes are organized as 5'UTR-[Rep1a-Rep1b-S-3a-3b-E-M-4b-4c-5a-5b-N-6b]-3'UTR, except in eight sequences lacking non-structural protein genes (accessory genes) 4b, 4c, and 6b. Seventeen sequences have auxiliary S2' cleavage site located 153 residues downstream the canonically conserved primary furin-specific S1/S2 cleavage site. The sequences distinctly cluster into lineages GI-1 (Mass-type; n = 8), GI-3 (Holte/Iowa-97; n = 2), GI-9 (Arkansas-like; n = 8), GI-13 (793B; n = 14), and GI-17 (California variant; CAV; n = 1), with regional distribution in Mexico; this is the first report of the presence of 793B- and CAV-like strains in the country. Various point mutations, substitutions, insertions and deletions are present in the S1 hypervariable regions (HVRs I-III) across all 5 lineages, including in residues 38, 43, 56, 63, 66, and 69 that are critical in viral attachment to respiratory tract tissues. Nine intra-/inter-lineage recombination events are present in the S proteins of three Mass-type sequences, two each of Holte/Iowa-97 and Ark-like sequence, and one each of 793B-like and CAV-like sequences. This study demonstrates the feasibility of FTA cards as an attractive, adoptable low-cost sampling option for untargeted discovery of avian viral agents in field-collected clinical samples. Collectively, our data points to co-circulation of multiple distinct IBVs in Mexican commercial flocks, underscoring the need for active surveillance and a review of IBV vaccines currently used in Mexico and the larger Latin America region.

Molecular investigation of S2-3a/3b-E-M-4b/4c-5a/5b-N gene of QX-like and variant genotype infectious bronchitis virus isolated in Thailand reveals a distinct E gene

The QX-like infectious bronchitis virus (IBV) and variant genotype have been discovered worldwide including Thailand. In order to know the origin of QX-like and variant genotype IBV in Thailand, the genetic analysis on multiple genes was investigated. Seven IBVs including four QX-like and three variant genotype were randomly selected from IBVs isolated in Thailand during 2008 and 2010. Phylogenetic analysis of the S2-3a/3b-E-M-4b/4c-5a/5b-N gene showed that Thai QX-like and variant genotype IBV were grouped together in a separate branch from other IBV strains. The isolated IBVs shared nucleotide identities of 96-99.9% with each other. They exhibited a high level of similarity (93.8%) with KM91 strain in South Korea. Phylogenetic analysis of the S2 and 3a/3b gene showed a relationship to KM91 strain. The E gene was distinct from other IBV strains. The M, 4a/4b and 5a/5b gene were closely related to Massachusetts type. The N gene was classified into two groups which were a group of unique to Thailand (variant genotype) and a relationship with Massachusetts type (QX-like). Recombination analysis identified the occurrence of recombination events in the genome of viruses. These findings demonstrated that the QX-like IBV and variant genotype isolates in Thailand were the recombinant viruses. Thai QX-like IBV had a genetic relationship with KM91 strain, Massachusetts type and unknown IBV whereas variant genotype had a genetic relationship with Thai QX-like IBV and Connecticut strain.

Molecular Characterization of QX-Like and Variant Infectious Bronchitis Virus Strains in Malaysia Based on Partial Genomic Sequences Comprising the S-3a/3b-E-M-Intergenic Region-5a/5b-N Gene Order

Infectious bronchitis virus (IBV) is one of the major poultry pathogens of global importance. However, the prevalence of IBV strains in Malaysia is poorly characterized. The partial genomic sequences (6.8 kb) comprising the S-3a/3b-E-M-intergenic region-5a/5b-N gene order of 11 Malaysian IBVs isolated in 2014 and 2015 were sequenced using next-generation sequencing technology. Phylogenetic and pairwise sequence comparison analysis showed that the isolated IBVs are divided into two groups. Group 1 (IBS124/2015, IBS125/2015, IBS126/2015, IBS130/2015, IBS131/2015, IBS138/2015, and IBS142/2015) shared 90%-95% nucleotide and deduced amino acid similarities to the QX-like strain. Among these isolates, IBS142/2015 is the first IBV detected in Sarawak state located in East Malaysia (Borneo Island). Meanwhile, IBV isolates in Group 2 (IBS037A/2015, IBS037B/2015, IBS051/2015, and IBS180/2015) were 91.62% and 89.09% identical to Malaysian variant strain MH5365/95 (EU086600) at nucleotide and amino acid levels, respectively. In addition, all studied IBVs were distinctly separate from Massachusetts (70%-72% amino acid similarity) and European strains including 793/B, Italy-02, and D274 (68%-73% amino acid similarity). Viruses in Group 1 have the insertion of three amino acids at positions 23, 121, and 122 of the S1 protein and recombinant events detected at nucleotide position 4354-5864, with major parental sequence derived from QX-like (CK-CH-IBYZ-2011) and a minor parental sequence derived from Massachusetts vaccine strain (H120). This study demonstrated coexistence of the IBV Malaysian variant strain along with the QX-like strain in Malaysia.

Molecular characteristic and pathogenicity analysis of a virulent recombinant avain infectious bronchitis virus isolated in China

A virulent infectious bronchitis virus (IBV), designated as CK/CH/GD/QY16 (referred as QY16), was isolated from a diseased chicken farm in Guangdong province, China, in 2016. The complete genome of the strain was sequenced and analyzed. The results show that the genome of QY16 consists of 27,670 nucleotides, excluding poly (A) tail, and that its genome organization is 5’ UTR-1a-1b-S-3a-3b-E-M-4b-4c-5a-5b-N-6b-3’ UTR-poly (A) tail. Sequence comparison among QY16 and other IBV strains was conducted and its results demonstrate that the S1 gene of QY16 has the highest nucleotide sequence identity with that of 4/91, and the other part of its genome is highly similar to that of YX10. The results of the phylogenic analysis show that the entire genome of QY16 and most of the QY16 genes are located in the same cluster as those of YX10, except for the S1 gene which is located in the same cluster with that of 4/91. It has been further confirmed by the RDP and SimPlot analysis that QY16 is a recombinant strain deriving from YX10 (as the major parental sequence) and 4/91 (as the minor parental sequence), and that the recombination occurs in a region which includes the 3’-terminal 1b sequence (85 nt) and the 5’-terminal S1 protein gene sequence (1,466 nt). The results of the vaccination-challenge test suggest that QY16 is a nephropathogenic strain of IBV and that the vaccine strains–H120 and 4/91—cannot provide effective protection against it. These results indicate that the continuing evolution of IBV strains by genetic drift and genetic recombination may lead to IBV outbreaks even among the vaccinated chickens in China.

Genetic diversity of avian infectious bronchitis virus in China in recent years

In this study, 213 infectious bronchitis viruses (IBVs) were isolated from samples collected from 801 flocks suspected to be infected with IBV from January 2016 to December 2017 in China. By using complete nucleotide sequences of S1 gene we determined the phylogeny of these IBV isolates, which in turn allowed us to define six lineages/genotypes, a number of recombinants and a novel variant. The GI-19 lineage was the most frequently isolated type in China in recent years. Although scattered mutations in the S1 gene among the GI-19 lineage viruses were observed, we also noted different sublineages in the GI-19 lineage with unique mutations, suggesting a high degree of S1 gene variation since they were first isolated in the mid-1990s. We also isolated a number of vaccine-like viruses from vaccinated diseased chickens, although more work is needed to differentiate the reisolation of vaccine strains from field strains of the same serotype. One of the important findings in this study is that the prevalence of the TW I type viruses in GI-7 lineage has been increasing in recent years in China. Another important finding is that recombination events occurred between the predominant GI-19 lineage and the commonly used 4/91 vaccine, which gave rise to distinct IBV isolates. In addition, a novel IBV isolate, together with a reference strain in GenBank database, were found to form a novel lineage/genotype that was remarkably distinct from established lineages. The characteristics of the antigenicity, tissue tropism, pathogenicity and complete genome were required for further investigation for the recombinants and the viruses in different sublineages and novel lineages. Meanwhile, permanent monitoring of circulating strains was needed to monitor the emerging viruses and rationally modify vaccination strategies in the field situation.

•

The GI-19 lineage was the most frequently isolated type in China in recent years.

•

Different sub-lineages were found in the GI-19 lineage with unique mutations.

•

A number of vaccine-like viruses were isolated from vaccinated diseased chickens.

•

The TW I type viruses in GI-13 lineage has been increasing in recent years in China.

•

A novel IBV isolate was found to form a novel clade that was distinct from established lineages.

Complete genome analysis of Iranian IS-1494 like avian infectious bronchitis virus

The nephropathogenic infectious bronchitis virus (IBV) strain IS-1494 like (variant-2; GI-23) was first isolated in the Middle East (1998). Despite intensive vaccinations, IS-1494 like IBVs are still circulating in Iran (the dominant genotype) and spread to other countries. Here, the full-length genome of this Iranian IS-1494 like IBV was (Mahed) determined to understand its evolutionary relationships. The genome consists of 27,652 nucleotides, with mutations in most of the structural genes. Thirteen open reading frames (ORFs) were predicted in the Mahed isolate (5' UTR-1a-1b-S-3a-3b-E-M-4b-4c-5a-5b-N-6b-3' UTR). ORFs 4b, 4c, and 6b, which has rarely been reported, were present in the Mahed genome. According to phylogenetic analysis of the full-length genome, 1a, S2, M, E, N protein, Mahed isolate clustered with the QX type strain. Based on the partial 1b, S1, Mahed clustered with the Q1 strain. The full-length genome of Mahed isolate shared the highest sequence homology with Gray and JMK (90.06-90.07%) and was least related to the Vic-s (86.21%). These data show that evolutionary variation because of recombination in IBV plays a major role in the adaptation and origin of IBV leading to new genetic and types of the virus strain.

Characterization of the complete genome, antigenicity, pathogenicity, tissue tropism, and shedding of a recombinant avian infectious bronchitis virus with a ck/CH/LJL/140901-like backbone and an S2 fragment from a 4/91-like virus

In this study, we isolated an infectious bronchitis virus, designated I1101/16, from broiler breeders in China. Analysis of the S1 gene showed that isolate I1101/16 was genetically close to strain ck/CH/LJL/140901, which belongs to the TW I genotype (also known as lineage GI-7 based on the recent IBV classification), however the S2 gene showed genetic diversity comparing to that of S1 gene. Comparison of the genomic sequences showed that the genome of isolate I1101/16 was similar to that of strain ck/CH/LJL/140901 from the 5' end of the genome to the 5' end of the S2 gene and from the 5' end of the 3a gene to the end of the genome, whereas the remaining parts of the genome sequences were more closely related to those of strain 4/91 than those of ck/CH/LJL/140901, thereby suggesting that recombination might have occurred during the origin of the virus. SimPlot and Bootscan analysis of the complete genomic sequence confirmed this hypothesis, where it showed that isolate I1101/16 arose through recombination events between ck/CH/LJL/140901- and 4/91-like viruses. Isolate I1101/16 and strain ck/CH/LJL/140901 shared identical amino acids in almost all five of their B cell epitopes, but the two viruses had a serotype relatedness value of 65, which is well below 80, i.e., the lower cutoff value for viruses of the same serotype. In addition, pathogenicity tests demonstrated that isolate I1101/16 was more pathogenic to 1-day-old specific-pathogen-free chickens than strain ck/CH/LJL/140901, according to analysis of the clinical signs, whereas strain ck/CH/LJL/140901 exhibited prolonged replication and shedding after challenge compared with isolate I1101/16. This study did not provide evidence that recombination can directly alter the antigenicity, virulence, replication, shedding, and tissue tropism of a virus, but it did show that recombination events are likely to be major determinants of viral evolution.

Infectious Bronchitis Virus Variants: Molecular Analysis and Pathogenicity Investigation

Infectious bronchitis virus (IBV) variants constantly emerge and pose economic threats to poultry farms worldwide. Numerous studies on the molecular and pathogenic characterization of IBV variants have been performed between 2007 and 2017, which we have reviewed herein. We noted that viral genetic mutations and recombination events commonly gave rise to distinct IBV genotypes, serotypes and pathotypes. In addition to characterizing the S1 genes, full viral genomic sequencing, comprehensive antigenicity, and pathogenicity studies on emerging variants have advanced our understanding of IBV infections, which is valuable for developing countermeasures against IBV field outbreaks. This review of IBV variants provides practical value for understanding their phylogenetic relationships and epidemiology from both regional and worldwide viewpoints.

Emergence of novel nephropathogenic infectious bronchitis viruses currently circulating in Chinese chicken flocks

The emergence of novel infectious bronchitis viruses (IBVs) has been reported worldwide. Between 2011 and 2014, eight IBV isolates were identified from disease outbreaks in northeast China. In the current study we analysed the S1 gene of these eight IBV isolates in addition to the complete genome of five of them. We confirmed that these isolates emerged through the recombination of LX4 and Taiwan group 1 (TW1) viruses at two switch sites, one was in the Nsp 16 region and the other in the spike protein gene. The S1 gene in these viruses exhibited high nucleotide similarity with TW1-like viruses; the TW1 genotype was found to be present in southern China from 2009. Pathogenicity experiments in chickens using three of the eight virus isolates revealed that they were nephropathogenic and had similar pathogenicity to the parental viruses. The results of our study demonstrate that recombination, coupled with mutations, is responsible for the emergence of novel IBVs.

A duplex RT-PCR assay for detection of H9 subtype avian influenza viruses and infectious bronchitis viruses

H9 subtype avian influenza virus (AIV) and infectious bronchitis virus (IBV) are major pathogens circulating in poultry and have resulted in great economic losses due to respiratory disease and reduced egg production. As similar symptoms are elicited by the two pathogens, it is difficult for their differential diagnosis. So far, no reverse transcription-polymerase chain reaction (RT-PCR) assay has been found to differentiate between H9 AIV and IBV in one reaction. Therefore, developing a sensitive and specific method is of importance to simultaneously detect and differentiate H9 AIV and IBV. In this study, a duplex RT-PCR (dRT-PCR) was established. Two primer sets target the hemagglutinin (HA) gene of H9 AIV and the nucleocapsid (N) gene of IBV, respectively. Specific PCR products were obtained from all tested H9 AIVs and IBVs belonging to the major clades circulating in China, but not from AIVs of other subtypes or other infectious avian viruses. The sensitivity of the dRT-PCR assay corresponding to H9 AIV, IBV and mixture of H9 AIV and IBV were at a concentration of 1×10¹, 1.5×10¹ and 1.5×10¹ 50% egg infective doses (EID₅₀) mL^-1, respectively. The concordance rates between the dRT-PCR and virus isolation were 99.1 and 98.2%, respectively, for detection of samples from H9N2 AIV or IBV infected chickens, while the concordance rate was 99.1% for detection of samples from H9N2 AIV and IBV co-infected chickens. Thus, the dRT-PCR assay reported herein is specific and sensitive, and suitable for the differential diagnosis of clinical infections and surveillance of H9 AIVs and IBVs.

Full-length genome sequencing analysis of avian infectious bronchitis virus isolate associated with nephropathogenic infection

Infectious bronchitis virus (IBV) produces infectious bronchitis (IB) disease in poultry worldwide. In spite of proper vaccinations against the IBV, new IBV strains are continually emerging worldwide. In this study, a new highly virulent nephropathogenic IBV strain named CK/CH/XDC-2/2013 was identified from a vaccinated flock with clinical signs of IB in the Jiangsu province of China. The full-length genome sequence of the isolate was 27,714 nucleotides long, and the genome was organized similarly to classical IBV strains. Minimum divergence, phylogenetic analysis, and distance matrix of the genome showed that the CK/CH/XDC-2/2013 isolate had the highest similarity to the IBV BJ strain. The spike glycoprotein (S) gene had the greatest similarity to the nephropathogenic BJ strain and showed an 8 amino acid insertion (YSNGNSDV) at 73 to 80 sites and 3 amino acid deletion at sites 126 to 128 compared to the IBV vaccine strains. A recombination analysis of the S gene showed that the new isolate evolved from the IBV BJ strain and the KM91 vaccine strain. An animal challenge experiment showed a mortality of 60 to 80% in early-age chickens by different inoculation routes. Pathological examinations of the kidneys revealed inflammation, distention with uric acid deposits, and tubular degeneration. It indicated that the CK/CH/XDC-2/2013 isolate has robust kidney tissue tropism, and new nephropathogenic IBV strains are continuously evolving in China.

Whole genome analysis of Japanese bovine toroviruses reveals natural recombination between porcine and bovine toroviruses

Bovine toroviruses (BToVs), belong to the subfamily Toroviridae within the family Coronaviridae, and are pathogens, causing enteric disease in cattle. In Japan, BToVs are distributed throughout the country and cause gastrointestinal infection of calves and cows. In the present study, complete genome sequences of two Japanese BToVs and partial genome sequences of two Japanese BToVs and one porcine torovirus (PToV) from distant regions in Japan were determined and genetic analyses were performed. Pairwise nucleotide comparison and phylogenetic analyses revealed that Japanese BToVs shared high identity with each other and showed high similarities with BToV Breda1 strain in S, M, and HE coding regions. Japanese BToVs showed high similarities with porcine toroviruses in ORF1a, ORF1b, and N coding regions and the 5′ and 3′ untranslated regions, suggestive of a natural recombination event. Recombination analyses mapped the putative recombinant breakpoints to the 3′ ends of the ORF1b and HE regions. These findings suggest that the interspecies recombinant nature of Japanese BToVs resulted in a closer relationship between BToV Breda1 and PToVs.

•

Recombination events between porcine and bovine torovirus were identified.

•

Recombinant breakpoints were mapped at ORF1b and HE coding regions.

•

These recombinant viruses are prevalent throughout Japan.

Altered pathogenicity of a tl/CH/LDT3/03 genotype infectious bronchitis coronavirus due to natural recombination in the 5'- 17kb region of the genome

An infectious bronchitis coronavirus, designated as ck/CH/LGX/130530, was isolated from an IBV strain H120-vaccinated chicken in this study. Analysis of the S1 gene showed that isolate ck/CH/LGX/130530 was a tl/CH/LDT3/03-like virus, with a nucleotide sequence similarity of 99%. However, a complete genomic sequence analysis showed that ck/CH/LGX/130530 was more closely related to a Massachusetts type strain (95% similarity to strain H120) than to the tl/CH/LDT3/03 strain (86%), suggesting that recombination might have occurred during the origin of the virus. A SimPlot analysis of the complete genomic sequence confirmed this hypothesis, and it showed that isolate ck/CH/LGX/130530 emerged from a recombination event between parental IBV H120 strain and pathogenic tl/CH/LDT3/03-like virus. The results obtained from the pairwise comparison and nucleotide similarity showed that the recombination breakpoint was located in the nsp14 gene at nucleotides 17055-17083. In line with the high S1 gene sequence similarity, the ck/CH/LGX/130530 isolate was serotypically close to that of the tl/CH/LDT3/03 strain (73% antigenic relatedness). Furthermore, vaccination with the LDT3-A vaccine, which was derived from the tl/CH/LDT3/03 strain by serial passaging in chicken eggs, provided good protection against challenge with the tl/CH/LDT3/03 strain, in contrast to the poor protection offered with the H120 vaccine. Interestingly, isolate ck/CH/LGX/130530 exhibited low pathogenicity toward specific-pathogen-free chickens compared with the nephropathogenic tl/CH/LDT3/03 strain, which was likely due to natural recombination in the 5' 17-kb region of the genome. Our results also indicate that the replicase gene of IBV isolate ck/CH/LGX/130530 is associated with viral pathogenicity.

Development and characterization of neutralizing monoclonal antibodies against the S1 subunit protein of QX-like avian infectious bronchitis virus strain Sczy3

Infectious bronchitis (IB) is a highly contagious disease in chickens caused by infectious bronchitis virus (IBV). The present study was carried out with the aim to develop anti-spike 1 (S1) subunit monoclonal antibodies (MAbs) that could react with IBV strains of different genotypes. The high antigenicity region of S1 gene of an QX-like IBV strain Sczy3 was amplified and ligated into the prokaryotic expression vector pET-32a(+), and the recombinant His-S1 fusion proteins were expressed and purified. The purified whole viral antigen of Sczy3 strain was used to immunize BALB/c mice to produce hybridoma-secreting anti-IBV MAbs. Eleven anti-IBV MAbs were generated, and two MAbs 1C8 and 2C10 were positive in indirect ELISA against both His-S1 protein and the purified whole viral antigen. These two MAbs showed positive reaction with IBV in Western blot, and the isotype were both IgM. These two MAbs react specifically with IBV but not with Newcastle disease virus (NDV) or avian influenza virus (AIV) subtype H9 or H5, and could cross-react with other 10 IBV strains in five different genotypes. End-point neutralizing assay performed in chicken embro kidney (CEK) cells revealed that the neutralization titer of 1C8 and 2C10 against Sczy3 reached 1:2.82 and 1:4.70, respectively. The anti-S1 MAbs produced in the present work may be valuable in developing an antigen-capture ELISA test for antigen detection or a competitive ELISA test for antibody detection or therapeutic medicine for IB in poultry.

Two novel neutralizing antigenic epitopes of the s1 subunit protein of a QX-like avian infectious bronchitis virus strain Sczy3 as revealed using a phage display peptide library

The spike (S) protein of the infectious bronchitis virus (IBV) plays a central role in the pathogenicity, the immune antibody production, serotype and the tissue tropism. In this study, we generate 11 monoclonal antibodies (mAbs) against S1 subunit of IBV Sczy3 strain, and two mAbs 1D5 and 6A12 were positive in indirect ELISA against both His-S1 protein and the purified whole viral antigen. MAb 6A12 and 1D5 could recognized by other 10 IBV strains (IBVs) from five different genotypes, except that 1D5 had a relatively low reaction with two of the 10 tested IBVs. End-point neutralizing assay performed in chicken embro kidney (CEK) cells revealed that the neutralization titer of 6A12 and 1D5 against Sczy3 reached 1:44.7 and 1:40.6, respectively. After screening a phage display peptide library and peptide scanning, we identified two linear B-cell epitopes that were recognized by the mAbs 1D5 and 6A12, which corresponded to the amino acid sequences (87)PPQGMAW(93) and (412)IQTRTEP(418), respectively, in the IBV S1 subunit. Sequences comparison revealed that epitope (412)IQTRTEP(418) was conserved among IBVs, while the epitope (87)PPQGMAW(93) was relatively variable among IBVs. The novel mAbs and the epitopes identified will be useful for developing diagnostic assays for IBV infections.

Genomic and single nucleotide polymorphism analysis of infectious bronchitis coronavirus

Infectious bronchitis virus (IBV) is a Gammacoronavirus that causes a highly contagious respiratory disease in chickens. A QX-like strain was analysed by high-throughput Illumina sequencing and genetic variation across the entire viral genome was explored at the sub-consensus level by single nucleotide polymorphism (SNP) analysis. Thirteen open reading frames (ORFs) in the order 5'-UTR-1a-1ab-S-3a-3b-E-M-4b-4c-5a-5b-N-6b-3'UTR were predicted. The relative frequencies of missense: silent SNPs were calculated to obtain a comparative measure of variability in specific genes. The most variable ORFs in descending order were E, 3b, 5'UTR, N, 1a, S, 1ab, M, 4c, 5a, 6b. The E and 3b protein products play key roles in coronavirus virulence, and RNA folding demonstrated that the mutations in the 5'UTR did not alter the predicted secondary structure. The frequency of SNPs in the Spike (S) protein ORF of 0.67% was below the genomic average of 0.76%. Only three SNPS were identified in the S1 subunit, none of which were located in hypervariable region (HVR) 1 or HVR2. The S2 subunit was considerably more variable containing 87% of the polymorphisms detected across the entire S protein. The S2 subunit also contained a previously unreported multi-A insertion site and a stretch of four consecutive mutated amino acids, which mapped to the stalk region of the spike protein. Template-based protein structure modelling produced the first theoretical model of the IBV spike monomer. Given the lack of diversity observed at the sub-consensus level, the tenet that the HVRs in the S1 subunit are very tolerant of amino acid changes produced by genetic drift is questioned.

Serotype shift of a 793/B genotype infectious bronchitis coronavirus by natural recombination

•

Infectious bronchitis virus causes a respiratory disease in domestic chickens worldwide.

•

Recombination is thought to contribute to the emergence of IBV variants.

•

Strain ck/CH/LHLJ/140906 is originated from recombination events between 4/91- and H120-like strains.

•

Recombination of the S1 domain resulted in the emergence of a novel serotype of IBV.

An infectious bronchitis coronavirus, designated as ck/CH/LHLJ/140906, was isolated from an infectious bronchitis virus (IBV) strain H120-vaccinated chicken flock, which presented with a suspected infectious bronchitis virus (IBV) infection. A phylogenetic analysis based on the S1 gene clustered ck/CH/LHLJ/140906 with the 793/B group; however, a pairwise comparison showed that the 5′ terminal of the S1 gene (containing hypervariable regions I and II) had high sequence identity with the H120 strain, while the 3′ terminal sequence was very similar to that of IBV 4/91 strain. A SimPlot analysis of the complete genomic sequence, which was confirmed by a phylogenetic analysis and nucleotide similarities using the corresponding gene fragments, suggested that isolate ck/CH/LHLJ/140906 emerged from multiple recombination events between parental IBV strains 4/91 and H120. Although the isolate ck/CH/LHLJ/140906 had slightly higher S1 amino acid sequence identity to strain 4/91 (88.2%) than to strain H120 (86%), the serotype of the virus was more closely related to that of the H120 strain (32% antigenic relatedness) than to the 4/91 strain (15% antigenic relatedness). Whereas, vaccination of specific pathogen-free chickens with the 4/91 vaccine provided better protection against challenge with ck/CH/LHLJ/140906 than did vaccination with the H120 strain according to the result of virus re-isolation. As the spike protein, especially in the hypervariable regions of the S1 domain, of IBVs contains viral neutralizing epitopes, the results of this study showed that recombination of the S1 domain resulted in the emergence of a new serotype.