key: cord-1043292-g8ky8fck
authors: Lian, Jiamin; Wang, Zhanxin; Xu, Zhouyi; Chen, Tong; Shao, Guanming; Zhang, Xinheng; Qin, Jianping; Xie, Qingmei; Lin, Wencheng
title: Distribution and molecular characterization of avian infectious bronchitis virus in southern China
date: 2021-03-27
journal: Poult Sci
DOI: 10.1016/j.psj.2021.101169
sha: c5217e0459a250a60b51706491331f445826a712
doc_id: 1043292
cord_uid: g8ky8fck

Avian infectious bronchitis virus (IBV) is causing considerable economic losses in the world poultry industry. The main difficulty of prevention and control of IB disease is the numerous genotypes and serotypes. The genetic analysis of IBV was mainly based on the S1 gene which played an important role in infectivity. In the study, One hundred and thirty-nine strains of avian infectious bronchitis virus were isolated from chickens showing signs of disease in southern China during the period from April 2019 to March 2020. The nucleotide and amino acid sequences from the isolated field strains were compared to 22 published references. Nucleotide homologies ranged from 64.5% to 100% and amino acid homologies ranging from 70% to 99.8%. Six genotype IBV strains were co-circulating in southern China. QX-type was still the most dominant genotype. Alignment of nucleotide and amino acid sequences of S1 gene revealed that the substitutions, insertions and deletions are widely among isolated strains. Recombination analysis showed that there is a large number of recombinant strains amongst these isolates, forming new sub branches, subtypes and variants. Therefore, long-term continuing surveillance is necessary for IBV prevention and control.

Avian infectious bronchitis (IB) is an acute and highly contagious disease characterized by upper respiratory symptoms, air sacculitis, proventriculitis, nephritis, enteritis and egg production drop, resulting in considerable economic losses in the poultry industry, globally (Yu et al., 2001; Cook et al., 2012; Sjaak de Wit et al., 2011) . The causative agent of the diseases avian infectious bronchitis virus (IBV), which belongs to the genus gamma Coronavirus of the family Coronaviridae of the order Nidovirales (Cavanagh, 2005) . Similar to other coronaviruses, avian infectious bronchitis virus genome consists of a linear, non-segmented, positivesense, and single-stranded RNA approximately 27 kilobases (kb) in length and encodes four structural proteins: nucleocapsid (N) protein, envelope (E) protein, membrane (M) glycoprotein, and spike (S) glycoprotein (Sutou et al., 1988; Liu et al., 2009) .

The glycoprotein S is a precursor protein, and subsequently cleaved into two non-covalently bound polypeptides (S1 and S2) by proteases (Stern and Sefton, 1982) . S1 subunit forming the tip of the viral spike, and S2 subunit anchoring the S1 into viral membrane. The S1 glycoprotein is the main infectious and pathogenic protein of IBV, and plays a critical role in viral attachment and neutralizing antibodies induction (Cavanagh, 1983; Cavanagh et al., 1986; Shang et al., 2018) . Interestingly, the differences in the S1 subunit contribute to significant divergences in the viral serotype as well as the poor cross-protection (Cavanagh et al., 1986; Koch et al., 1990) , leading to serotype-specific host immunity against IBV. Currently, the rapid variation of S1 gene and the poor cross-protection between the different serotypes pose significant challenges to control this disease.

Considering the antigenic variability of IBV and the poor cross-protection of commercial vaccine against different IBV strains, the characterization of IBV strains based 1 on genotype or serotype is beneficial to prevent this disease (Jones et al., 2004) . However, the lack of standard positive serum makes the identification of serotype impossible. The variability in the S1 protein is considered as the basis for genotype and serotype classification, definition of antigen characteristics, induction of poor crossprotection of commercial vaccines, and emergence of IBV variants. Thus, genotype classification based on the S1 gene becomes the primary method for classifying IBV strains because of the rapid evolution of S1 gene and the importance of S1 in virus-host interactions and immune responses (Mo et al., 2013; Valastro et al., 2016) . Till now, at least 32 genotypes have been identified based on the S1 gene worldwide, including QX-type, 4/91-type, LDT3-type, HN08-type, TW-type, TC07-2-type and Mass-type (Yu et al., 2001; Ji et al., 2011) .

Long-term continuing surveillance is significantly important to control this disease. With the aim of understanding the epidemiology and the characteristics of IBV strains circulating in southern China, a long-term surveillance program was performed in southern China. This study covered the monitoring data from April 2019 to March 2020. Here, we described in detail the results of this study, including the distribution, characterization, and recombination events of IBV isolates in southern China. The present work provides information on the molecular epidemiology of IBV, which will be beneficial for the control of IBV in Chinese poultry industry.

All animal experiments were approved by the Committee of the Ethics on Animal Care and Experiments at South China Agricultural University (approval ID: SYXK-2019-0136). Permission was issued for the field studies, and all of the owners of the farms gave informed consent to conduct the study on this site.

From April 2019 to March 2020, an epidemiological survey was performed on yellow-feathered broiler chicken farms in southern China, including Guangdong, Guangxi, Fujian, Yunnan, Hunan, Hubei, Sichuan, Jiangsu, Anhui and Zhejiang provinces. All the farms were visited monthly to collect clinical samples (including kidney and trachea) from sick birds displaying respiratory signs, nephritis and egg production drop.

A total of 420 swab samples were collected for IBV detection. Virus isolation was performed using SPF chicken embryonated eggs according to the previous description (Han et al., 2011) . Briefly, clinical samples were frozen and thawed three times, homogenized in PBS containing penicillin (100 U/mL) and streptomycin (100 mg/mL), clarified by centrifugation, and sterilized by filtration using a 0.22-mm filter. The supernatant was inoculated into the 9-day-old SPF chicken embryonated eggs by the allantoic sac route, and blind passage till the occurrence of embryonic lesions. Virus purification was performed using the limiting dilution method (Hebberecht et al., 2019) . Briefly, the allantoic fluid was filtered by the filter, and added 0.1 mL into 0.9 mL PBS for 10 times to obtain 10 À1~1 0 À8 viral dilutions. The 9day-old SPF chicken embryonated eggs was inoculated with each dilution, and harvested after 72 hours. According to the pathological characteristics and the reverse transcription-polymerase chain reaction (RT-PCR) assay results, the highest dilution of allantoic fluid confirmed IBV positive was used to inoculating chicken embryonated eggs three times.

The presence of IBV was identified and verified by RT-PCR as previously described (Liu et al., 2006; Feng et al., 2018) . Briefly, specific primer pairs (sense primer: 5'-CCTAAGAACGGTTGGAAT-3' and antisense primer: 5'-TACTCTCTACACACACAC-3') were designated based on the M gene sequence, RNA was extracted from the allantoic fluids using a Viral DNA/ RNA Miniprep Kit (Axygen, USA) following the manufacturer's instructions, and subjected to RT-PCR for IBV detection.

Cloning and Sequencing of S1 Gene Specific primer sets for amplifying the entire S1 gene (sense primer: 5'-AACTGAACAAAAGACCGACT-3' and anti-sense primer: 5'-CAAAACCTGCCATAAC-TAACA-3') were designated based on the nucleotide sequence of LX4 strain (Accession No. AY338732). RT-PCR was carried out to amply the entire S1 gene using a high-fidelity pfu DNA polymerase (Q5 TM High-Fidelity DNA polymerase, NEB) under the following conditions: 98°C for 10 s, 45 cycles of 98°C for 10 s, 56°C for 30 s, and 72°C for 30 s, and 72°C for 10 min.The amplicons were cloned into a pEASY-Blunt plasmid, and sequenced commercially (Sangon Biotech, Shanghai). Sequence analyses were carried out using the DNAStar software (Madison, WI, USA).

Sequence alignments were performed with nucleotide sequences of IBV isolates and other IBV strains retrieved from GenBank using DNASTAR software (Madison, WI, USA). Phylogenetic analysis was performed based on the nucleotide sequences of S1 gene by neighbor-joining method by using MEGA version 7.0. The bootstrap values were determined from 1,000 replicates of the original data. All the reference IBV strains used in this study are provided in Table 1 .

Putative recombination analysis was performed based on the S1 gene of IBV isolates using the Recombination Detection Program 4 (RDP v.4.97) software. A total of 9 methods (including RDP, GENECONV, BootScan, MaxChi, Chimaera, SiScan, PhyLPro, LARD and 3Seq) were employed to analyze the putative recombination events. Recombination events only supported by five of these methods with a P-value adjusted to 0.05 which were regarded as positive. The potential recombination events were further verified by SimPlot software (version 3.5.1).

From April 2019 to March 2020, a total of 420 clinical samples were collected from sick chickens in southern China. The onset age of sick chickens mainly ranged between 5 and 60 days. All the IBV strains were isolated from chickens showing typical respiratory clinical signs such as difficulty in breathing, nasal discharge, sneezing, tracheal rales, watery eyes, lethargy, and mucosal thickening with serous or catarrhal exudates in the nasal passages, sinuses, and trachea, or showing typical nephritis symptoms such as enlarged and pale kidneys, frequently with urate deposits in the tubules, severe dehydration and weight loss. All the 420 field strains were confirmed IBV positive by observation of curled and dwarfed embryos, RT-PCR and sequencing, and 139 IBV strains were identified.

The information of IBV isolation and identification was statistically analyzed: 139 field IBV strains circulating in poultry industry were isolated from dead or diseased chickens in 10 provinces in southern China, including Yunnan (n = 26), Hubei (n = 29), Hunan (n = 12), Guangxi (n = 15), Jiangsu (n = 8), Fujian (n = 16 ), Guangdong (n = 12), Zhejiang (n = 5), Anhui (n = 8) and Sichuan (n = 8). Yunnan, Hubei, Fujian, and Guangxi were the four regions with the largest number of IBV isolates ( Figure 1B ). Moreover, 109 strains (78.4%) were collected in winter and spring seasons ( Figure 1A ).

Due to the importance of S1 protein in serotyping and cross-protection, S1 gene of IBV isolates were sequenced for diversity analysis. As a result, a total of 420 nucleotide sequences of S1 gene were obtained, but only 139 sequences were identified because many IBV isolates from the same chicken flock had identical sequences. Within the S1 genes obtained in this study, the nucleotide sequences were 1614~1640 nt in length, and the corresponding deduced amino acid (aa) numbers of S1 protein were 538 aa (2.2%, 3/139), 539 aa (5.8%, 8/ 139), 540 aa (39.6%, 55/139), 541 aa (5%, 7/139), 542 aa (7.9%, 11/139), 543 aa (18%, 25/139), 544 aa (7.2%, 10/139), 545 aa (9.4%, 13/139) and 546 aa (5%, 7/139) (Supplementary material Table 2 ). Interestingly, the diverse S1 genes contained mutations, insertions and/or deletions, resulting in the differences of nucleotides and different numbers of nucleotides and deduced amino acids. All the sequences of S1 gene of IBV isolates have been submitted to GenBank, and the accession numbers and information of IBV strains were listed (Supplementary material Table 1) .

Because the protease cleavage site motif is essential for cleaving during viral maturation period , we analyzed the motif sequence of the cleavage site in this study. As a result, a total of 8 evolution of IBV in south China. The cleavage site motifs of IBV have become complex and diverse in many regions with the spread of IBV in worldwide. Multiple sequence alignment based on the nucleotide sequence of S1 gene of the 139 IBV studied strains from this study and other strains retrieved from GenBank was performed using the Clustal X program. The nucleotide homology of S1 gene between isolates was 64.5% to 100%, while that between isolates and reference strains was 70% to 99.8%. These results indicated that the S1 gene of the field isolates was highly variable.

Furthermore, several insertions, deletions and mutations were observed in S1 gene of 139 isolated field strains (Supplementary material Table 2 ). Consistent with previous reports, three hypervariable regions (HVRs) were found according to the sequence analysis of S1 gene. Most variation in the deduced amino acid sequences of 139 field strains occurred in the regions 3-79, 61-136 and 154-432 (comparison with the S1 sequence of H120 strain). Generally, most of the isolates had deletions of N at site 25, G/S at site 130, G at site 96, and S at site 123.

To determine the genetic relationships among the IBV strains, phylogenetic analysis was performed using the nucleotide sequences of the S1 gene of the 139 isolates and 22 reference strains. As a result, 132 IBV isolates and most reference strains were clustered into six groups: QX-type (66.9%, 93/139), 4/91-type (13.7%, 19/139), TW-type (6.5%, 9/139), LDT3-type (5.8%, 8/ 139), HN08-type (0.7%, 1/139), TC07-2-type (1.4%, 2/ 139). Additionally, seven field isolates CK/CH/ YNML04-3, CK/CH/GXNN12-1, CK/CH/GDCG12-2, CK//CH/JSHA12-1, CK/CH/AHWH11-1, CK/CH/ YNXY04-1 and CK/CH/FJFS12-8 were clustered into individual genotypes without reference strains based on their S1 genes and termed variants (Figure 2) .

QX-type strains were the most prevalent strain among the surveyed provinces ( Figure 1B) . The proportion and geographical distribution indicated that QXtype strain was the predominant genotype of IBV circulating in southern China. All the QX-type isolates shared a high degree of nucleotide (92%) and amino acid (92.6%) homology with QX-type vaccine strains, but shared low homology at nucleotide (86%) and amino acid (86.2%) levels with other vaccine strains (H120, H52, Ma5, 4/91, and LDT3-A).

The 4/91-type isolates were the second prevalent strain in southern China. A total of 19 strains (13.7%) were isolated from Guangxi, Hunan, Hubei, Fujian, Guangdong and Sichuan. The homology of nucleotide and amino acid among the 4/91-type isolates were 93.5% and 93%, respectively. The homology of nucleotide and amino acid sequences between the 4/91-type isolates and vaccine strains were 95.5% and 93.7%, respectively. In contrast, the 4/91-type isolates shared a lower degree of nucleotide (80%) and amino acid (78.8%) similarity with the vaccine strains (H120, H52, Ma5, QX and LDT3-A).

The TW-type isolates were the third most prevalent strain in southern China. A total of 9 strains were isolated from Guangxi, Hubei, Yunnan and Zhejiang. The nucleotides and amino acid homologies between TWtype isolates and TW2575/98 reference strain were 92.8%~99.4% and 94.6%~97.8%, but their nucleotide and amino acid homology with other Chinese mainland strains and vaccine strains (H120, H52, 4/91, QX and LDT3) were lower than 84.9% and 83.3% respectively. These observations may suggest different serotypes and Figure 2 . Phylogenetic tree of 139 isolates and 22 reference strains for S1 genes of IBV (starting at the AUG translation initiation codon and ending at the cleavage recognition motifs). A) Phylogenetic tree of S1 gene of isolates from April 2019 to March 2020 (Traditional). B) Phylogenetic tree of S1 gene of isolates in from April 2019 to March 2020 (Circle). The phylogenetic tree was constructed using the MEGA version 5.0 by the neighbor joining method with No. of differences model and setting bootstrap 1,000 replicates. The reference strains are marked with~. Seven variant strains that were marked with . Unmarked strains are field strains isolated in this study. low protection between vaccine strains and TW-type isolated strains.

The other 11 isolates belonged to LDT3-type (n = 8), TC07-2-type (n = 2), and HN08-type (n = 1). The LDT3-type isolates shared high identity with reference strains LDT3 at nucleotide and amino acid levels, and TC072-type isolates shared highly identity with TC07-2 strains. In contrast, the homology of nucleotide and amino acid sequences between the HN08-type isolates and vaccine strain HN08 was 88.7% and 87.1%, respectively.

Strains were considered as recombinants if any crossover event appeared between two putative parental strains. It has been reported that the recombination events can occur in S1 gene of IBV (Wang et al., 1993) , so we have identified putative recombination events in this study. As a result, evidence of recombination was found in 7 IBV strains using RDP package ( Figure 3 ). As shown in Table 2 , the variant CK/CH/GDCG12-2 was derived from recombination between SAIBK (major parent) and CK/CH/AHHF12-2 (minor parent), with recombination breakpoints mapping to positions 106 (beginning breakpoint) and 629 (ending breakpoint); The variant CK/CH/YNML04-3 was derived from recombination between CK/CH/GXNN11-3 (major parent) and CK/CH/HBHC09-3 (minor parent), the recombination region located in the nucleotide site 646 −1468, with 99.1% similar to CK/CH/HBHC09-3 (4/ 91-type); The isolate CK/CH/GXNN12-1 appears to be a recombinant between strain CK/CH/HBXN08-3 as the major parent and strain CK/CH/HBXN08-1 as the The S1 gene of CK/CH/YNML04-3 was from a recombination of the reference strain CK/CH/GXNN11-3 (TW-type) and reference strain CK/CH/HBHC09-3 (4/91-type). (B) The S1 gene of CK/CH/GXNN12-1 was from a recombination of the reference strain CK/CH/HBXN08-3 (LDT3-type) and reference strain CK/CH/HBXN08-1 (4/91-type). (C) The S1 gene of CK/CH/ GDCG12-2 was from a recombination of the reference strain SAIBK (HN08-type) and reference strain CK/CH/AHHF12-2 (4/91-type). (D) The S1 gene of CK//CH/JSHA12-1 was from a recombination of the reference strain CK/CH/JSYC1912-6 (QX-type) and reference strain CK/CH/ HBXN08-3 (LDT3-type). (E) The S1 gene of CK/CH/AHWH11-1 was from a recombination of the reference strain A2 (QX-type) and reference strain HN08 (HN08-type). (F) The S1 gene of CK/CH/HBYNXY04-1 was from a recombination of the reference strain 793 (491-type) and reference strain CK/CH//YNXY04-3 TW-type). (G) The S1 gene of CK/CH/FJFS12-8 was from a recombination of the reference strain TW2575-98 (TWtype) and reference strain CK/CH/YNSL04 (QX-type). The y-axis gives the percent identity within a sliding window 200 bp wide centered on the position plotted, with a step size between plots of 20 bp. minor parent; the recombination occurring between nucleotides 723−1091, which was 99.7% similar to CK/ CH/HBXN08-1 (4/91-type), while the rest of the S1 gene was 97.8% similar to CK/CH/HBXN08-3 (LDT3type. The variant CK/CH/JSHA12-1 was derived from recombination between CK/CH/JSYC1912-6 (major parent) and CK/CH/HBXN08-3 (minor parent), with recombination breakpoints mapping to positions 1395 (beginning breakpoint) and 1615 (ending breakpoint). The variant CK/CH/AHWH11-1 was a recombinant between A2 (major parent) and HN08 (minor parent), the recombination occurred between nucleotides 762-845; The isolate CK/CH/HBYNXY04-1 appears to be a recombinant with strain CK/CH/YNXY04-3(TWtype) as the major parent and 793 (4/91-type) as the minor parent, with recombination breakpoints mapping to positions 646 (beginning breakpoint) and 1621 (ending breakpoint). The isolate CK/CH/FJFS12-8 appears to be a recombinant with strain TW2575-98 as the major parent and CK/CH/YNSL04 (QX-type) as the minor parent, with recombination breakpoints mapping to positions 1101 (beginning breakpoint) and 1270 (ending breakpoint) ( Table 2) .

There are two major factors accounting for the genetic diversity of IBV. Firstly, mutations, insertions and deletions frequently occurred in the IBV genomes due to the inaccuracy of the coronavirus RNA-dependent RNA polymerase (RDRP) (Mahmood et al., 2011; Feng et al., 2017) . Secondly, unique random template switching during RNA replication may result in a high frequency of homologous RNA recombination . Thereby, the variant strains and new genotypes of IBV are continuously emerging, worldwide.

The commercial attenuated vaccines against IBV (including QX, H120, H52, Ma5, and 4/91), have been widely used in Chinese poultry industry. However, IB disease still frequently break out in commercial poultry farms in China due to the poor cross-protection Lim et al., 2011a) . The genotypes or serotypes of IBV strains circulating in Chinese poultry industry are different from the currently used vaccine strains (Yu et al., 2001; Xu et al., 2007) . Currently, the commercial vaccines against IBV used in southern China are attenuated strains, and always used with multi-genotypes combination, such as VG / GA + H120 + 491 + LDT3, LaSota + H120 + 491 and Clone30 + H120 + D90 + 491 (Xu et al., 2021) . However, the attenuated strains in commercial vaccines still need to be updated due to the evolutionary adaptation of the virus strain. Therefore, the long-term molecular epidemiological surveillance is of great significance for the prevention and control of IB disease.

This study is a part of our ongoing surveillance program, which focuses on investigation of IBV molecular epidemiology and identification of circulating IBV Table 2 . S1 gene recombination in variant IBV strains. DISTRIBUTION AND MOLECULAR CHARACTERIZATION OF IBV strains in commercial poultry farms in south China. In this study, a total of 139 IBV field strains were isolated from diseased or dead chickens in southern China during April 2019 to March 2020. Most diseased chickens were found at the age of 10-25 days, and few chickens were infected with IBV before 10 days. We speculated that the maternal antibody might provide a degree of protection against the circulating strains. Most of IBV infected broilers showed typical respiratory and nephritis symptoms, and few IBV infected broiler showed air sacculitis and peritonitis symptoms, indicating that the respiratory and nephron pathogenic IBV strains remained to be prevalent in south China (Li et al., 2010; Jackwood, 2012) .

S protein was cleaved into S1 and S2 subunits when it was in viral maturation period, and the protease cleavage site motif was essential for cleaving . In the study, a total of 8 cleavage site motifs were identified in 139 IBV isolates, including HRRRR, RRFRR, HRHRR, HRRKR, RRLRR, RRFSR, and RRCRR. Previous studies showed that HRRRR was the main existing cleavage site motif of isolated IBV strains in China (Bing et al., 2007; Li et al., 2010) . Consistent with these reports, HRRRR (69/139), RRF(S) RR (38/139), HRHRR (21/ 139) and HRRKR (7/139) were the predominant cleavage site motifs. Interestingly, three cleavage site motifs (RRLRR, RRFSR and RRCRR) were first identified in the present study, indicating the continuously evolution of IBV in south China. S1 glycoprotein containing the neutralizing epitopes could induce the neutralizing antibodies and involve the serotype evolution, phenotype change and genetic diversity. In this study, we found that S1 genes of the IBV isolates possess different nucleotide in length, with series of mutations, insertions and deletions. It has been reported that the serotype differences of IBVs generally correlated with variations in HVR Hypervariable region of S1 protein, and the division of genotype was partly according to the HVR (Cavanagh and Davis, 1988; Moore et al., 1997; Cavanagh, 2007) . Considering the variations in these IBV isolates, the genetic mutation of IBV HVRs may provide the explanation about the evolution of IBV in south China.

There are 6 different groups of IBV strains circulating in southern China, including QX-type, HN08-type, TWtype and LDT3-type, 4/91-type and Mass-type strains. Besides China, QX IBV has spread all over the world in the past 20 years, causing serious losses to laying hens and broilers on commercial farms (Franzo et al., 2017) . In 1996, QX-type strain was first isolated in Shandong province, China (Kiss et al., 2015) . Our long-term surveillance indicated QX-type strain was mainly associated with typical nephritis and pseudostratified clinical symptoms during 2011 to 2015 (Feng et al., 2017) . In this study, we found that QX is the most predominant genotype in southern China, accounting for 66.9% of the isolates from April 2019 to March 2020. All these data provide evidence that QX-type strain is the dominant genotype in China. Currently, vaccines against QX-type strain are widely used to prevent and control IBV (Ji et al., 2020; Ren et al., 2020) . However, the recombination and mutation frequently occurred under immunization pressure and natural selection, resulting in the continuous emergence of IBV variants; The 4/91-type and LDT3-type strains accounted for 13.7% and 5.8% of the total isolates, respectively. Both genotypes have their corresponding commercial vaccine, which are approved for use in Chinese poultry industry, since vaccination reduces the incidence of IB disease, but also increases the risk of recombination with field strains. The common circulation of multiple IBV genotypes and the increase in the proportion of IBV variants accelerate the possibility of virus recombination (Ali et al., 2018; Yan et al., 2019; Zhang et al., 2020) , which poses a great challenge to IBV vaccination control. Therefore, it is necessary to continue monitoring the prevalence and evolution of the 4/91-type and LDT3-type IBV strains and other co-circulating viral strains in the field. TWtype IBV strains accounted for 6.5% of the field isolates in this study. Consistent with previous studies, the proportion of TW-type strains increased year by year (Yan et al., 2016) . Therefore, it is urgent to develop new TW type IB vaccine to prevent the outbreak of TW-IB.

Recombination is an important mechanism for IBV evolution Su et al., 2016) . The unique discontinuous transcription system and the viral polymerase "jumping" contribute to the high RNA recombination frequency in the genome of IBV (Zou et al., 2010) . The recombination always occurs when the viral polymerase switches from one template to another during genomic synthesis in a host . In this study, seven variants with different parental genotypes were detected by recombination analysis. The recombination events took place between the 4/91 and TW, LDT3 or HN08 genotype, and between the QX and TW, LDT3 or HN08 genotype, producing new variant strains and subgroups. These results suggested that recombination between different genotypes may contribute to the emergence of novel IBV variants. Additionally, our data indicated that 4/91-type strains involved the recombination events, providing clues that 4/91 has become an important gene donor in the genetic evolution of IBV. More effort will be required to monitor the recombination of IBV in the future.

Currently, the attenuated vaccine strains against IBV are widely used to control IB disease globally, controlling the occurrence of IB disease to a certain extent. However, the immune selective pressure is considered as an important reason for the occurrence of novel genotypes or serotypes of IBV . IBV possesses highly genomic mutation and recombination characteristics under immune selective pressure, causing the development of new genotypic groups under high-intensity vaccination (Pohuang et al., 2011) . Actually, the probability of gene recombination between vaccine strains and epidemic strains is increasing in recent years (Lim et al, 2011b) . More effort will be required to focus on the occurrence of IBV novel strains in the future.

This work was supported by the Guangdong Provincial Promotion Project on Preservation and Utilization of Local Breed of Livestock and Poultry (4000-F18260), National Natural Science Foundation of China (31802206, 32072842) , and the Construction Project of Modern Agricultural Science and Technology Innovation Alliance in Guangdong Province (2020KJ128).

Ethical Approval and Consent to Participate: Our animal research was conducted under the guidance of the SCAU's Institutional Animal Care and Use Committee. The chicken sampling procedures were approved by the Animal Care and Use Committee of Guangdong Province, China.

Availability of Data and Materials: The nucleotide sequences of 139 IBV isolated strains have been deposited in GenBank under the accession no. MT766911 to MT769003 and MW042820-MW042867. The authors declare that all data supporting the findings of this study are available within the article.

We declare that we do not have any commercial or associative interest that represents a conflict of interest in connection with the work submitted.

Supplementary material associated with this article can be found in the online version at doi:10.1016/j. psj.2021.101169.

Safety and efficacy of attenuated classic and variant 2 infectious bronchitis virus candidate vaccines

Different genotypes of nephropathogenic infectious bronchitis viruses co-circulating in chicken population in China

Coronavirus IBV: structural characterization of the spike protein

Coronaviruses in poultry and other birds

Coronavirus avian infectious bronchitis virus

Evolution of avian coronavirus IBV: sequence of the matrix glycoprotein gene and intergenic region of several serotypes

Coronavirus IBV: virus retaining spike glycopolypeptide S2 but not S1 is unable to induce virus-neutralizing or haemagglutination-inhibiting antibody, or induce chicken tracheal protection

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

The long view: 40 years of infectious bronchitis research

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

Epidemiology and characterization of avian infectious bronchitis virus strains circulating in southern China during the period from 2013-2015

Think globally, act locally: phylodynamic reconstruction of infectious bronchitis virus (IBV) QX genotype (GI-19 lineage) reveals different population dynamics and spreading patterns when evaluated on different epidemiological scales

A 15-year analysis of molecular epidemiology of avian infectious bronchitis coronavirus in China

Single genome sequencing of near full-length HIV-1 RNA using a limiting dilution approach

Review of infectious bronchitis virus around the world

Spike glycoprotein cleavage recognition site analysis of infectious bronchitis virus

Epidemiological investigation of avian infectious bronchitis and locally determined genotype diversity in central China

Phylogenetic distribution and predominant genotype of the avian infectious bronchitis virus in China during

Possible North African progenitor of the major European infectious bronchitis variant (793B, 4/91, CR88)

Survey indicates circulation of 4/91 and QXtype infectious bronchitis viruses in Hungary in 2014 -Short communication

Antigenic domains on the peplomer protein of avian infectious bronchitis virus: correlation with biological functions

Origin and evolution of Georgia 98 (GA98), a new serotype of avian infectious bronchitis virus

Isolation and genetic analysis revealed no predominant new strains of avian infectious bronchitis virus circulating in South China during

An emerging recombinant cluster of nephropathogenic strains of avian infectious bronchitis virus in Korea

Live attenuated nephropathogenic infectious bronchitis virus vaccine provides broad cross protection against New Variant Strains

Genetic diversity of avian infectious bronchitis coronavirus strains isolated in China between

Molecular characterization and pathogenicity of infectious bronchitis coronaviruses: Complicated evolution and epidemiology in China caused by cocirculation of multiple types of infectious bronchitis coronaviruses

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

Molecular characterization of major structural protein genes of avian coronavirus infectious bronchitis virus isolates in Southern China

Identification of amino acids involved in a serotype and neutralization specific epitope within the s1 subunit of avian infectious bronchitis virus

Sequence analysis of S1 genes of infectious bronchitis virus Isolated in Thailand during

Pathogenicity of a QX-like avian infectious bronchitis virus isolated in China

Cryo-EM structure of infectious bronchitis coronavirus spike protein reveals structural and functional evolution of coronavirus spike proteins

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

Coronavirus proteins: structure and function of the oligosaccharides of the avian infectious bronchitis virus glycoproteins

Epidemiology, genetic recombination, and pathogenesis of coronaviruses

Cloning and sequencing of genes encoding structural proteins of avian infectious bronchitis virus

S1 gene-based phylogeny of infectious bronchitis virus: An attempt to harmonize virus classification

Evidence of natural recombination within the S1 gens of infectious bronchitis virus

Isolation and identification of four infectious bronchitis virus strains in China and analyses of their S1 glycoprotein gene

An attenuated TW-like infectious bronchitis virus strain has potential to become a candidate vaccine and S gene is responsible for its attenuation

Pathogenicity of a TW-like strain of infectious bronchitis virus and evaluation of the protection induced against it by a QX-like strain

Molecular characteristics and pathogenicity analysis of QX-like avian infectious bronchitis virus

Characterization of three infectious bronchitis virus isolates from China associated with proventriculus in vaccinated chickens

Molecular characterization of variant infectious bronchitis virus in China, 2019: Implications for control programmes

Genetic analysis revealed LX4 genotype strains of avian infectious bronchitis virus became predominant in recent years in Sichuan area