key: cord-0004278-rqsn0g8q
authors: Li, Shuyun; Du, Lijing; Xia, Jing; Du, Jiteng; You, Guojin; Wen, Yiping; Huang, Xiaobo; Zhao, Qing; Han, Xinfeng; Yan, Qigui; Wu, Rui; Cui, Min; Cao, Sanjie; Huang, Yong
title: Antigenic and Pathogenic Characteristics of QX-Type Avian Infectious Bronchitis Virus Strains Isolated in Southwestern China
date: 2019-12-13
journal: Viruses
DOI: 10.3390/v11121154
sha: a4f95bc90ea589d614f0e2a70daf60fbd7d23c31
doc_id: 4278
cord_uid: rqsn0g8q

The QX-type avian infectious bronchitis virus (IBV) is still a prevalent genotype in Southwestern China. To analyze the antigenicity and pathogenicity characteristics of the dominant genotype strains (QX-type), S1 gene sequence analysis, virus cross-neutralization tests, and pathogenicity test of eight QX-type IBV isolates were conducted. Sequence analysis showed that the nucleotide homology between the eight strains was high, but distantly related to H120 and 4/91 vaccine strains. Cross-neutralization tests showed that all eight strains isolated from 2015 and 2017 belonged to the same serotype, but exhibited antigenic variations over time. The pathogenicity test of the five QX-type IBV isolates showed that only three strains, CK/CH/SC/DYW/16, CK/CH/SC/MS/17, and CK/CH/SC/GH/15, had a high mortality rate with strong respiratory and renal pathogenicity, whereas CK/CH/SC/PZ/17 and CK/CH/SC/DYYJ/17 caused only mild clinical symptoms and tissue lesions. Our results indicate that the prevalent QX-type IBVs displayed antigenic variations and pathogenicity difference. These findings may provide reference for research on the evolution of IBV and vaccine preparation of infectious bronchitis (IB).

Infectious bronchitis (IB) is a highly contagious disease that causes significant economic losses to the poultry industry worldwide [1, 2] . Furthermore, chickens infected with infectious bronchitis virus (IBV) also become susceptible to secondary infections with bacteria, or other pathogens, because of the tracheal cilia damage, and display a higher mortality rate [3] . IBV, an etiological agent of IB, is an enveloped virus with a single-stranded positive-sense non-segmented RNA genome of approximately 27.6 kb in length and belongs to the Gammacoronavirus genus [4] . Due to the incomplete proofreading mechanism of coronavirus RNA polymerase and gene recombination during viral replication, new genotypes and serotypes of IBV variant strains appear continuously [5, 6] . An array of serotypes and strains of IBV infect chickens exist throughout the world. Therefore, continuous testing of pathogenicity and serotype determination of new isolates in regions and/or countries remains crucial for better epidemiological understanding and control of IB. Of these isolates, the Sczy3 strain was identified as QX-type and SCZY3-serotype by a previous study (GenBank number: JF732903.1) [12] . The other seven field strains were isolated from clinical samples (trachea, lung, and kidney) of dead/diseased chickens from Southwestern China that displayed respiratory symptoms and/or nephritis during 2015-2017. Virus isolation was carried out by being inoculated into the allantoic cavity of 9-to-10-day-old specific-pathogen-free (SPF) chicken embryos with 0.2 mL of 10% tissue homogenates. The embryos were incubated at 37 • C and examined twice daily for their viability. The allantoic fluids wereharvestedafter36 h of incubation, and three blind passages were conducted. Detection methods of presence of IBV were determined based on the previous literature [12] . The 50% embryo infectious doses (EID 50 ) were determined as previously described (Reed and Muench, 1938) [22] . SPF chicken embryos were obtained from the Beijing Merial Vital Laboratory Animal Technology Co., Ltd. (Beijing, China) and were hatched at our laboratory (Lab of Chicken Infectious Disease, College of Veterinary Medicine, Sichuan Agricultural University, Chengdu, China). The details of the seven field isolates are shown in Table 1 . 

Viral RNA was extracted from fresh allantoic-fluid-infected IBVs, using a TRIzol agent, as per the manufacture's recommendation (Invitrogen, Carlsbad, CA, USA). Reverse transcription was performed, using a reverse-transcriptase (RT) kit (Takara Bio-Inc, Dalian, China), according to the manufacturer's protocol. PCR amplification of IBV S1 gene sequencing was referred to in the previous study [12, 15] . The nucleotide sequences of the S1 gene, including isolates and reference strains, were aligned by using the ClustalW method in MEGA version 7.0.14 [23] . The information of 18 references is shown in Table S1 . The phylogenetic tree of the S1 gene created by a neighbor-joining tree method was drawn in MEGA version 7.0.14 [24] . Bootstrap values were determined from 1000 replicates. Recombination events of the potential recombinant strain were estimated, using the Recombination Detection Program (RDP 4.72) and SimPlot version 3.5.1 software, as previously reported [12] .

Chicken antisera were prepared against the eight strains. Briefly, fresh allantoic of 10 4 EID 50 IBVs were inactivated, using 0.1% formal, at 37 • C, for 18 h. The inactivated mixture was emulsified with an equal volume of complete Freund's adjuvant (Sigma, Saint Louis, MO, USA). Following this, 0.1 mL of the mixture was used for subcutaneous injection in four-week-old SPF chickens (n = 3) as first immunization. Subsequently, the same amount of antigen was emulsified in Freund's incomplete adjuvant (Sigma, Saint Louis, MO, USA) for the following two booster injections (at two-week intervals). A serum sample of each immunized chicken was collected 14 days after immunization and inactivated at 56 • C for 30 min.

CEK cells were prepared and cultured as previously described [15] . Primary CEK cells were prepared from the kidneys of SPF chicken embryos at 18-20 days old and maintained in Dulbecco's Modified Eagle Medium (DMEM; Gibco, Grand Island, NY, USA) supplemented with 10% fetal bovine serum (FBS; Zhejiang Tian-hang Biological Technology Stock Co., Ltd., Zhejiang, China) [25] . Those IBVs were propagated in primary CEK cells cultured in DMEM supplemented with 2% FBS at 37 • C under 5% CO 2 . The supernatant was harvested at 40 h post-inoculation and passaged blindly in CEK cells until a characteristic cytopathic effect (CPE), such as syncytia, was observed. The presence of IBV in the supernatant of the CEK culture was verified by RT-PCR, as previously reported [15] . Determination of the TCID 50 of CEK-adapted IBVs in CEK cells was conducted as per the method reported by Reed and Muench (1938) 

The antigenic relatedness (R) of those eight strains was determined via double-direction viral cross-neutralization (VN) tests with constant viral titers and diluted serum in CEK cells [12] . Briefly, sera were serially diluted two-fold with sterile phosphate-buffered saline (PBS), and then equal volumes of 100 TCID 50 of the CEK-adapted strains were added. The mixtures were then incubated at 37 • C for 1 h. Virus-serum mixtures were used to infect CEK cell cultures in 24-well plates (6 wells for each dilution) and incubated for 72 h. Following this, 50% end-point neutralizing titers were calculated, using the method described by Reed and Muench (1938) . A negative serum sample was always included in the analysis.

The cross-reactivity R values were calculated according to a previously reported method, representing the relationship between antigens: R ≥ 0.70 proves antigenic identity, 0.70 > R ≥ 0.33 proves antigenic relatedness (minor subtype difference), and 0.33 > R ≥ 0.11 proves loose relatedness (major subtype difference), whereas R < 0.11 indicates no relatedness at all (serotype difference) [26, 27] . Antigenic cartography was also performed by using the program AntigenMap (http://sysbio.cvm.msstate.edu/AntigenMap), which uses matrix completion multidimensional scaling to map microneutralization titers in two dimensions [28] . Viruses with high antigenic relevance clustered closely on the map, while viruses with low antigenic relevance stayed far away from each other. Detailed settings were set as follows: Low Reactor Threshold: 20; Normalization Method: N1; Temporal Model: Normal; Rank: 2; Number of Iterations: 2000. 

Viral RNA was extracted from fresh allantoic fluid, using TRIzol (Invitrogen, Carlsbad, CA, USA) according to the manual instructions. Viral cDNA was transcribed by using the reverse-transcriptase (RT) kit (Takara Bio-Inc, Dalian, China), following the manufacturer's instructions. The fragments of IBV strains were amplified by using SF and SR by PCR reaction. The amplicons were then inserted into pMD-19T cloning vector (Takara Bio-Inc, Dalian, China) and transformed into Escherichia coli TOP10 competent cells (Takara Bio-Inc, Dalian, China). The recombinant plasmid DNA was extracted, and its copy numbers were calculated according to a previous report [29] .

The RT-qPCR was performed by using the SYBR Premix Ex Taq TM II (Takara Bio-Inc, Dalian, China). The 20 µL reaction mix consisted of 10 µL of SYBR ® Premix Ex Taq TM II (TliRNaseH Plus, Takara), 0.8 µL of each specific primer (QX-F and QX-R: 10 pmol/µL), 1 µL of standard plasmid, and 7.4 µL of RNase-free ddH 2 O. The optimum thermal cycling parameters were as follows: 95 • C for 30 s, and 40 cycles of 95 • C for 10 s, 54 • C for 15 s, and 72 • C for 20 s. The dissolution curve was generated by the default parameters at the end of the cycle. Under optimum thermal cycling parameters, the standard curves of RT-qPCR were established by linear regression analysis of the threshold cycle (Ct) value (y axis) versus the log of the initial copy number presented in each sample dilution (x axis) [30] .

Five QX-type IBVs, GH/15, DYW/16, MS/17, DYYJ/17, and PZ/17, were chosen to infect SPF chickens. A total of 122 ten-day-old SPF chickens were randomly divided into six groups (A-E, n = 22 chickens each; F, n = 12 chickens). The birds of Groups A-E were inoculated with 100 µL of five QX-type strains (GH/15, DYW/16, MS/17, DYYJ/17, and PZ/17) containing 10 5 EID 50 , via intraocular and intranasal routes. The birds of Group F were inoculated with sterilized PBS as control. Chickens were held in separate biosafety level 2 (BSL2) isolators in the Laboratory Animal Center of Sichuan Agricultural University (Ya'an, Sichuan, China) with ad libitum access to feed and water; and they were maintained under uniform standard management conditions. The challenge schedule is shown in Table S2 .

Ten chickens in Groups A-E and six chickens in Group F were used for the observation of clinical signs and tissue lesion. In detail, IB clinical symptoms, including sneezing, tracheal rales, chills, and increased drinking, were observed and recorded daily, until 8 days post-challenge (DPC). The morbidity and mortality of birds were calculated at the end of the study period.

These chickens that died during the experimental period in each group were necropsied immediately, and the survival birds in each group were euthanized and necropsied at 8 DPC. Samples of the trachea and kidney were collected and saved in 10% neutral formalin for 24 h, embedded in paraffin, and stained with hematoxylin and eosin (HE) before being observed under a standard light microscope. Gross lesions of trachea were scored as follows: 0 points, normal; 1 point, abundant mucus and minor bleeders in the trachea; and 2 points, large hemorrhage in the trachea (Figure 1 ). Gross lesions of the kidney were scored as follows: 0 points, normal; 1 point, kidney swelling; and 2 points, mottled kidney ( Figure 1 ). The mean lesion scores (MLS) in the trachea and kidney were calculated for each group. 

The twelve chickens in Groups A-E and six chickens in Group F were used for the quantification of virus shedding; four chickens were randomly selected for sampling the trachea and kidney at 4, 6, and 8 DPC. The samples of trachea and kidney were collected and processed as follows: after flushing the tracheal mucosa with 3 mL of PBS and grinding the kidney with 2 mL PBS, the fluids were collected and stored at −70 • C for further analysis. After freezing, thawing, and centrifugation at 1000× g for 30 min at 4 • C, three times, 200 µL of the supernatant was used for RNA extraction. Virus copies were determined by RT-qPCR, as described above.

Data analysis of survival, lesion scores, and RT-qPCR were performed by using IBM SPSS Statistics 22 software. Statistical significance was considered as follows: probability (p) value < 0.05 was considered statistically significant.

All animal experiments were conducted in compliance with protocols approved by the Sichuan provincial Laboratory Animal Management Committee (Permit Number: XYXK (Sichuan) 2014-187). The protocols for this experiment were performed according to the guidelines of the Ethics and Animal Welfare Committee (EAWC) of the Sichuan Agricultural University.

Phylogenetic trees were constructed based on ORF nucleotide sequences of the S1 gene ( Figure 2 ). All eight IBV field strains were QX genotype, sharing 94.3%-99.2% nucleotide identity of the S1 gene with each other and 90.6%-98.6% with other QX-type reference strains. The S1 gene identities of these two IBV isolates, DYW/16 and PZ/17, were highest 99.2%, while the lowest similarity was 94.3% between Sczy3 and XSH/17 isolates. Additionally, the nucleotide identity for the S1 gene was only 76.5%-78.7% between the eight IBV isolates and commercial vaccine strains (H120 and 4/91). Recombination analysis of the S1 gene sequence indicated no recombination event in the S1 genes of the seven field strains.

On analyzing the amino acid sequence, the S1 domain showed 92%-98.3% identity, and the hypervariable regions (HVR) were located at residues 23-65 and 120-393 among the eight isolates. Amino acid substitution between S (Ser) residues and other amino acid residues was frequently observed, with the highest substitution rate occurring in S and N (Asn). Interestingly, the S1 protein of DYW/16 had six amino acid differences at the sites 150 D (Aspartic acid), 151N (Asparagine), 163W (Tryptophan), 193I (Isoleucine), 240K (Lysine), and 241C (Cystine), compared with that of the other six field strains and Sczy3, H120, and 4/91.

Phylogenetic tree was constructed based on the ORF nucleotide sequences of the S1 gene from seven field strains (filled triangles) and 18 reference strains, using MEGA version 7.0.14, with the neighbor-joining method and 1000 bootstrap replicates. The scale bar corresponds to 0.05 estimated amino acid substitutions per site.

In terms of the neutralization titer of the eight QX-type IBVs, the immune sera of all strains could neutralize each other, but the neutralization ability varied with different IBV strains. Immune sera of DYYJ/17 showed high neutralizing ability against MS/17, XSH/17, DYYJ/17, GH/15, PZ/17, Sczy3, and DYW/16 strains. Sczy3 and DYW/16 antisera showed high neutralizing ability against at least four strains. However, antisera of PZ/17 and GH/15 showed low neutralizing ability against other isolates; their mean neutralization titer was significantly lower than that of other antisera ( Table 2) .

The seven QX-type isolates were classified into the same serotype as the Sczy3 strain, a representative strain of the SCZY3-serotype identified as the predominant serotype in Southwestern China during 2012-2016 [12, 25] . Therefore, seven QX-type isolates could be clustered into the SCZY3-serotype, but the XSH/17 strain showed a loose relatedness with Sczy3 (0.11 ≤ R < 0.32) ( Table 3) . Of these seven IBV isolates, MS/15 and GH/15 showed a loose relatedness with DYW/16 (0.11 ≤ R < 0.32) ( Table 3 ). AntigenMap analysis showed that eight QX-type strains (including Sczy3) could be classified into three antigenic subgroups: six IBVs including the five IBV field isolates MS/17, XSH/17, DYYJ/17, GH/15, and PZ/17 and Sczy3 strain were grouped into subgroup 2, while MS/15 while MS/15 strainwas placed in subgroup 1 and DYW/16 strain belonged to subgroup 3 ( Figure 3 ). 

Standard curves were plotted with the log of target gene copy numbers as the x axis and cycle threshold (Ct value) as the y axis. The efficiency of the standard curves was 94%, the correlation coefficient (R2) of the standard curve was 1.0000, and the concluding equation of the standard curve was Y = − 3.659lgX + 4.16 (y axis is the threshold cycle (Ct); x axis is the initial copy number present in each sample). To verify the specificity of a particular primer, a dissolution curve was generated at the end of the cycle, and the result showed only the target product as a single melt peak.

At 2-3 DPC, infected chickens in all groups began to display clinical symptoms such as sneezing, listlessness, and huddling. Different morbidity rates were observed in different groups: DYW/16, GH/15, MS/17, and DYYJ/17 showed 100% morbidity, and PZ/17 showed 70% morbidity. Dead chickens were observed in DYW/16, GH/15, and MS/17 groups at 3-6 DPC, whereas no dead chickens were observed in DYYJ/17 and PZ/17 groups ( Figure 4A ).The percentage survival of DYW/16isolates was 50% during the eight-day observation period. Significant differences were analyzed between DYW/16 group and group of DYYJ/17 and PZ/17 (p < 0.05). Survival rates of MS/17 and GH/15 isolates were 70% and 80%, respectively. There were no clinical symptoms or deaths observed in the control group birds.

For routine gross anatomy in challenge groups, obvious pathogenic differences in the trachea and kidney were observed in five QX-type fields groups, while no gross lesion was observed in the control group. However, different degrees of gross lesions were clearly observed in each group by calculating MLS ( Figure 4B The microscopic lesions in the trachea and kidney collected from each of the challenged chickens were further characterized. Except for the control group, microscopic lesions were both observed the groups DYW/16, MS/17, GH/15, PZ/17, and DYYJ/17, which were typically characterized as degeneration and necrosis of the ciliated epithelial cells (Figure 4C ), as well as degeneration and necrosis in the renal tubules of the kidney ( Figure 4D) . The difference in trends of microscopic lesions between the five groups was consistent with the result of gross lesions observation in each group. 

The viral load in the trachea and kidneys from the infected groups are shown in Figure 5 . For the DYW/16 strain, the tracheal viral load peak was observed at six DPC, and the kidney viral load peak was observed at four DPC. The viral loads of the DYW/16 strain in the trachea were significantly higher than other strains at 6-8 DPC (p < 0.05), and the viral loads in the kidney was significantly higher than other strains at four DCP. For the MS/17 strain, the peak of tracheal viral loads was six DCP, but the most virus loads of kidneys at eight DCP. The viral loads of the MS/17 strain were higher than other strains in infected chicken kidneys at 6-8 DCP (p< 0.05). For the GH/15, PZ/17, and DYYJ/17 strains, the viral loads in the trachea and kidney peaked at six DCP. In terms of total viral loads in the trachea and kidney, the DYW/16 strain showed higher values than other strains in the trachea, and the MS/17 strain showed higher values than other strains in the kidney. The viral loads of the DYYJ/17 in the trachea strain were lower than that of other strains, and the viral loads of the PZ/17 strain in the kidney were the lowest during the experiment period. Figure 5 . Virus loads of trachea and kidney in 12 chickens per group after challenge IBV were measured by RT-qPCR. Statistical differences between different groups were determined by using the analysis of one-way ANOVA method. * p < 0.05 indicates a significant difference.

IB is an economically important, highly contagious, acute, upper-respiratory tract disease of chickens and other fowl, caused by the avian coronavirus IBV [31, 32] . Although multiple genotypes of IBV have continuously been detected, e.g., QX-type, TW-type, and TC07-2, QX-type IBV is still the prevalent genotype in China [12, 15, 33] . Therefore, it is important to understand the antigenicity and pathogenicity of QX-type IBV field isolates for the prevention and control of IB. It has been reported that the antigenicity and virulence vary between different genotypes [12, 20] , and no systematic comparison among QX-type IBVs has been conducted. In our study, eight QX-type IBVs were isolated from H120 and 4/91 vaccinated chickens in China, and their antigenicity and pathogenicity were analyzed for understanding the epidemiology and evolution of IBVs.

The S1 phylogenetic analyses revealed that eight QX-type isolates shared 94.3%-99.2% nucleotide identity with each other and were distantly related to H120 and 4/91 strains with 76.5%-78.7% nucleotide identity. The cross-reactivity R values indicated that the seven QX-type isolates were classified into the same SCZY3-serotype, which was different from H120 and 4/91 (R < 0.11) [12, 25] . It has been reported that some QX-type strains showed distant relationships with other representative strains of the classical IBV serotypes [20] . The results indicated that different QX-type isolates were important agents for IBV outbreaks in China, and a new vaccine strain should be used for the prevention and control of IB. In the present study, higher neutralization ability against all analyzed strains came from immune serum of field isolates DYYJ/17 and DYW/16, implying that these strains have the potential to be used for the development of vaccine candidates.

The program AntigenMaphas been shown to effectively detect antigenic variationevents, and it is useful in antigenic characterization, which was applied for influenza, human enteroviruses, noroviruses, and porcine epidemic diarrhea virus [28, 34] . In our study, eight QX-type IBVs with 94.3%-99.2% nucleotide identity were classified into the same serotype, but showed a minor difference in antigenicity among those strains. AntigenMap analysis showed that subgroup 2, which includes the Sczy3 strain and field isolates MS/17, XSH/17, DYYJ/17, GH/15, and PZ/17 had tight antigenicity, which has been majorly isolated since 2017. Meanwhile, subgroup 1 (MS/15 isolate) and 3 (DYW/16 isolate) were separated in 2015 and 2016. These results indicated that there areantigenic variations in different IBV strains over time, andminor changes in the S1 gene maycause antigenic variation. For example, the two strains of DYW/16 and PZ/17 shared 99.2% nucleotide identity, but were classified into different subgroups. However, the data of QX-IBVs used in this study are limited, and this map will be expanded further when more data becomes available in the future. It has been reported that antigenic relatedness values were influenced by hypervariable regions (HVR) in the S1 subunit of IBV S protein [8, 35] . Therefore, antigenic variations mapped by AntigenMap may be related to the nucleotide mutations in HVR, which were located at residues 23aa-65aa and 120aa-393aa in the S1 protein of the eight QX-type isolates. Furthermore, amino acid substitution between S (Ser) residues and other amino acids' residues, such as N (Asn) and P (Pro) in HVR, were frequently observed among the seven field isolates, but whether this has an effect on antigenic variations is unclear.

As we all know, pathogenicity of strains of IBV is variable, depending on the strain, age at exposure, and route of inoculation. In the present study, five QX-type isolates showed different pathogenicity to ten-old-day SPF chickens. The DYW/16 strain caused a mortality rate of 50%, and the MS/17 and GH/15 strains were respectively associated with a mortality rate of 30% and 20%, whereas the DYYJ/17 and PZ/17 strains only caused clinical symptoms and tissue lesion. At necropsy, the DYW/16, MS/17, and GH/15 strains could cause more server lesions in tracheas and kidneys than the DYYJ/17 and PZ/17 strains. However, the S1 gene between five QX-type IBVs with different pathogenicity shared 95.2%-99.2% identical. Such as the S1 gene, similarity is 99.2% between the DYW/16 and PZ/17 strains, indicating some key nucleotide mutation in the IBV S1 gene can lead to difference in pathogenicity [21] . Previous studies have shown that the binding properties of S proteins have been hypothesized to be associated with differences in pathogenicity [36] . The receptor-binding domain (RBD) in the spike protein of the M41 strain overlapped with HVR in the S1 gene was critical for attachment in the spike protein of IBV, and Y43H change in RBD of recombinant ArkDPC S1 proteins allows for enhanced binding to the trachea [10, 35] . HVRs of five QX-type field isolates were located at residues 23aa-65aa and 120aa-393aa, and the six unique amino acids (D150, N151, W163, I193, K240, and C241) in HVRs of the virulent DYW/16 strain were observed, which may perform an important role in the pathogenicity.

The pathogenicity of IBV is determined by several factors, which may involve replication efficacy, tissue tropism, and the ability to deal with the host's immune system. In the present study, RT-qPCR examination and tissue lesions both showed that the elevated levels of viral replication efficacy and tissue lesions in the trachea and kidney are associated with increased pathogenicity of IBV strains.The mechanisms of tissue tropism characteristics in different IBVs were elucidated, nephropathogenicB1648 strain disseminated to internal organs via a cell-free and -associated viremia with KUL01 + monocyticcells as important carrier cells, resulting in serious kidney lesions, whereas respiratory M41 does not do this [37] . The mechanism of replication in monocytic cells might be recruited to cause kidney lesion and even high mortality in chickens infected with five QX-type IBVs. Besides, it has been reported that QX-type IBV caused greater innate immune responses and induced greater apoptosis in CEK than the M41 strain, which is also a major contributing factor to the pathological outcome [38] . Therefore, QX-type IBVs have a strong ability to combat the host immune system, causing tissue lesion and even death. However, the mechanism responsible for the differences in the combat abilityremains unclear for different QX-type IBVs. Although, the S gene, replicase gene, and 5a accessory gene have also been reported to be responsible for IBV virulence [11, 39] , this still cannot completely explain the mechanism behind the virulence differences of IBVs. The key gene and amino acid residues associated with IBV virulence and antigenicity should be verified via new methods, such as the reverse-genetics technique.

In summary, QX-type isolates were important agents for IBV outbreaks in China; their antigenic and pathogenic characteristics were analyzed in this study. The results showed that the eight QX-type IBV isolates belonged to the same serotype, but the antigenic variation over time was observed. Different pathogenicity, replication efficacy, and tissue lesions in 10-day-old SPF chickens were also observed among QX-type isolates. In addition, the eight strains were distantly associated with H120 and 4/91 vaccine strains with 76.5%-78.7% nucleotide identity, suggesting the need for a new vaccine strain to prevent and control IB.The immune serum of two IBVs, DYYJ/17 and DYW/16, displayed good neutralization activity; therefore they have the potential to be developed as candidate vaccine strains.

Supplementary Materials: The following are available online at http://www.mdpi.com/1999-4915/11/12/1154/s1: Table S1 : The information of 18 IBV reference; Table S2 : The challenge schedule of 5 QX-genotype IBVs. 

The authors declare no conflicts of interest.

Isolation, propagation, identification and comparative pathogenicity of five Egyptian field strains of Eimeriatenella from broiler chickens in five different provinces in Egypt

Assessing the economic burden of avian infectious bronchitis on poultry farms in Brazil

Identification of a novel recombinant virulent avian infectious bronchitis virus

Safety and efficacy of an attenuated Chinese QX-like infectious bronchitis virus strain as a candidate vaccine

Nidovirales: A new order comprising Coronaviridae and Arteriviridae

Serotype and genotype diversity of infectious bronchitis viruses isolated during 1985-2008 in Guangxi

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

Glycosylation of the viral attachment protein of avian coronavirus is essential for host cell and receptor binding

Mapping of the receptor-binding domain and amino acids critical for attachment in the spike protein of avian coronavirus infectious bronchitis virus

Effects of Hypervariable Regions in Spike Protein on Pathogenicity, Tropism, and Serotypes of Infectious Bronchitis Virus

S gene and 5a accessory gene are responsible for the attenuation of virulent infectious bronchitis coronavirus

Phylogenetic analysis of the S1 glycoprotein gene of infectious bronchitis viruses isolated in China during

Phylogenetic and antigenic analysis of avian infectious bronchitis virus in southwestern China

Evolution of infectious bronchitis virus in China over the past two decades

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

Characterization and analysis of the full-length genome of a strain of the European QX-like genotype of infectious bronchitis virus

The pathogenesis of a new variant genotype and QX-like infectious bronchitis virus isolated from chickens in Thailand

Monne, I. S1 gene-based phylogeny of infectious bronchitis virus: An attempt to harmonize virus classification

Progress and Challenges toward the Development of Vaccines against Avian Infectious Bronchitis

Analysis of antigenicity and pathogenicity reveals major differences among QX-like infectious bronchitis viruses and other serotypes

Comparative features of infections of two Massachusetts (Mass) infectious bronchitis virus (IBV) variants isolated from Western Canadian layer flocks

A Simple Method of Estimating Fifty Per Cent Endpoints

Molecular Evolutionary Genetics Analysis version 6.0

The neighbor-joining method: A new method for reconstructing phylogenetic trees

Preparation and protective efficacy of a chicken embryo kidney cell-attenuation GI-19/QX-like avian infectious bronchitis virus vaccine

Antigenic diversity of psittacine herpesviruses: Cluster analysis of antigenic differences obtained from crossâ-neutralization tests

Antigenic variation of LaSota and genotype VII Newcastle disease virus (NDV) and their efficacy against challenge with velogenic NDV

AntigenMap 3D: An online antigenic cartography resource

Dynamic changes of virus load in supernatant of primary CEK cell culture infected with different generations of avian infectious bronchitis virus strains Sczy3 as revealed by real-time reverse transcription-polymerase chain reaction

Development of a SYBR Green real-time RT-PCR assay for the detection of avian encephalomyelitis virus

Effect of pH on the growth and cytopathogenicity of avian infectious bronchitis virus in chick kidney cells

Sequence analysis of nephropathogenic infectious bronchitis virus strains of the Massachusetts genotype in Beijing

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

Evolution, antigenicity and pathogenicity of global porcine epidemic diarrhea virus strains

Polymorphisms in the S1 spike glycoprotein of Arkansas-type infectious bronchitis virus (IBV) show differential binding to host tissues and altered antigenicity

Binding of avian coronavirus spike proteins to host factors reflects virus tropism and pathogenicity

Productive replication of nephropathogenic infectious bronchitis virus in peripheral blood monocytic cells, a strategy for viral dissemination and kidney infection in chickens

Pathogenicity and tissue tropism of infectious bronchitis virus is associated with elevated apoptosis and innate immune responses

The replicase gene of avian coronavirus infectious bronchitis virus is a determinant of pathogenicity