key: cord-328032-unwehrr5 authors: Zhou, Haisheng; Zhang, Meihong; Tian, Xue; Shao, Hongxia; Qian, Kun; Ye, Jianqiang; Qin, Aijian title: Identification of a novel recombinant virulent avian infectious bronchitis virus date: 2017-01-03 journal: Vet Microbiol DOI: 10.1016/j.vetmic.2016.12.038 sha: doc_id: 328032 cord_uid: unwehrr5 The emergence of new infectious bronchitis virus (IBV) variants is often disastrous in the poultry industry. In this study, an IBV, CK/CH/2010/JT-1, was isolated from an H120- and 4/91-IBV-vaccinated flock in China. Antisera against vaccine strains H120 and 4/91 could not provide effective protection against CK/CH/2010/JT-1 in virus neutralization assays. CK/CH/2010/JT-1 could cause 43.75% mortality with respiratory and severe renal lesions in inoculated chickens. Phylogenetic analysis of the S1 gene showed that CK/CH/2010/JT-1 and 31 other isolates could be grouped as a new genotypic cluster. Recombination analysis revealed that three recombination events could be found in the genome of CK/CH/2010/JT-1 at positions 24709-365, 17160-19811 and 21136-21770. Whole-genome sequence analysis showed that CK/CH/2010/JT-1 originated from multiple template switches among QX-like, CK/CH/LSC/99I-, tl/CH/LDT3/03- and 4/91-type IBVs. All of these data demonstrated that CK/CH/2010/JT-1 is a new recombinant genotype IBV with high virulence. Our findings suggest that the surveillance of new genotype strains of IBV is very important for developing more effective anti-IBV strategies. Infectious bronchitis (IB) is a highly contagious disease caused by infectious bronchitis virus (IBV), which belongs to the family Coronaviridae and the genus Coronavirus. IBV infection can cause clinical pathological signs mainly in the respiratory tract, kidney and reproductive tract of chickens, resulting in huge economic losses in affected flocks (Cook et al., 2012) . Little or no crossprotection occurs between different serotypes of IBV, and the increasing number of new serotypes of IBV is a major challenge for the prevention and control of IB (Cavanagh, 2007) . Although the IBV vaccine plays a vital role in controlling IB, it is still an epidemic across the world. IBV variants have been continuously emerging in China Xue et al., 2012; Zhao et al., 2014) . The nature of the large, single-RNA genome of IBV makes the virus variable through recombination or mutation. Recombination can cause the emergence and evolution of different IBV genotypes and different species of coronaviruses (Jackwood et al., 2010) . More and more recombination events have been reported in IBV, and they are distributed throughout the entire genome (Brooks et al., 2004; Kuo et al., 2013; Thor et al., 2011) . Recently, a new recombinant cluster of nephropathogenic IBVs able to cause death, respiratory signs and nephritis in infected chickens emerged in Korea (Lim et al., 2011 (Lim et al., , 2015 . Therefore, studies on the function of recombination in the antigenicity and pathogenicity of IBV are very important, potentially allowing IBV evolution to be predicted and better strategies for IBV control to be developed. In this study, an IBV, CK/CH/2010/JT-1, was isolated from a vaccinated flock in China, and its complete genomic sequence was determined. Genomic sequence analysis and pathogenicity studies revealed that CK/CH/2010/JT-1 is a virulent recombinant strain that belongs to a novel IB genotype. In this study, we identified and isolated an IBV (CK/CH/2014/JT-1) from a broiler flock vaccinated with IBV H120 and 4/91, in which the birds showed severe respiratory symptoms and some of them died. Gross examination showed lesions in the different organs, especially in trachea and kidney. The morbidity and mortality in this flock were as high as approximate 50% and 10% respectively. Trachea, lung and kidney specimens of five 25-day-old chickens were collected from the broiler flock. Sample suspensions (10% w/ v) of tissues prepared in sterile phosphate-buffered saline (PBS) were used for virus isolation . Ten passes were performed, and characteristic embryo changes, such as dwarfing, stunting, curling or death of the embryos, were observed between 2 and 7 days post-inoculation (dpi). Cross virus neutralization tests as described in the OIE manual (http://www.oie.int/international-standard-setting/terrestrialmanual/access-online/) were performed using anti-sera against Mass serotype vaccine strain H120 and 4/91 serotype vaccine strain 4/91 to determine their antigenic relationship. In brief, the viral titres of the Mass-type M41 and 4/91-type CK/CH/2014/TM98 strains were determined by inoculation of 10-fold dilutions into groups of five 10-day-old embryonated specific-pathogen-free (SPF) chicken eggs. The b method of cross-virus neutralization testing was performed using constant (10 2 50% embryo infectious doses [EID 50 s]) viral titres and 2-fold diluted serum against H120 and 4/91 in SPF chickens embryos. The EID 50 and the end-point of each serum sample were calculated using the methods of Reed and Muench. Viral RNA was extracted from 200 ml of allantoic fluid from the inoculated eggs using TRIzol reagent (Introvigen, Grand Island, USA). First-strand cDNA was synthesized using random hexamers (Introvigen, Grand Island, USA). The complete genomic sequence of strain CK/CH/2010/JT-1 was amplified with 26 pairs of primers, including the 5 0 -and 3 0 -terminal segments. Primer sequences and locations were designed corresponding to the genomic sequence of strains ck/CH/IBTZ/2012 and ck/CH/LJL/110302 (GenBank accession numbers KF663559 and KC136209, respectively). The polymerase chain reaction (PCR) conditions for amplification included 95 C for 5 min; 30 cycles of 94 C for 1 min, 55 C for 1 min for P1-P16 or 52 C for 1 min for P17-P23, 72 C for 2 min; followed by 72 C for 10 min. The 5 0 -and 3 0 -ends of the viral genome were amplified using 5 0 -and 3 0 -random amplification of cDNA ends (RACE) Kits (TaKaRa, Dalian, China) according to the operating instructions. The amplified products were analysed by 1.0% agarose gel electrophoresis, then purified with a Gel Extraction Kit (Qiagen, Hilden, Germany) and ligated into the pGEM-T Easy vector (Promega, Madison, USA). Each fragment of the viral genome was sequenced at least three times, and the consensus sequence was determined. The complete genome and gene sequences of isolate CK/CH/ 2010/JT-1 and of IBV reference strains obtained from GenBank were aligned and analysed using the ClustalW multiple alignment method in the MegAlign program of DNASTAR software (version 7.1; DNAstar, Madison, USA). A phylogenetic tree was constructed from the nucleotide sequences using MEGA4.0 software (www. megasoftware.net) with UPGMA or neighbour-joining statistical methods. The phylogeny test method chosen was bootstrap with 1000 bootstrap replicates. Evolutionary distances were computed by the pairwise distance method using the maximum composite likelihood model. BLASTN analysis was performed by services available on http://blast.ncbi.nlm.nih.gov/Blast.cgi The aligned nucleotide sequences of the complete genomes were analysed with the Recombination Detection Program (RDP4, Version 4.36) to detect potential within-gene recombination events. Seven detection methods in RDP v.4.36, including RDP, GENECONV, BootScan, MaxChi, Chimaera, SiScan and 3Seq, were used to confirm the recombination events. Only transferred gene fragments where at least 5 detection methods p-value 1 Â10 À14 were accepted. Recombination events and recombination breakpoints were further confirmed by simplot and bootscan analysis using the SimPlot program (version 3.5.1.). Nucleotide identity was performed by the Kimura (2-parameter) method with a transitiontransversion ratio of 2, and the window width and step size were 200 and 20 bp, respectively. The genomic sequence of IBV strain CK/CH/2010/JT-1 was submitted to the GenBank database and assigned the accession number KU361187. The GenBank accession numbers of strains M41 and CK/CH/2014/TM98 are DQ834384 and KU361199, respectively. Twenty-one 3-day-old SPF white leghorn chickens were randomly divided into two groups. Sixteen birds in group 1 were inoculated intranasally with 10 5 EID 50 s of IBV strain CK/CH/2010/ JT-1 in 0.2 ml PBS. Five birds in group 2 were inoculated with 0.2 ml PBS as non-infected controls. The birds were housed in isolators. All birds were observed daily for signs of disease (e.g., dishevelled feathers, depression, respiratory signs or diarrhoea) and death for 28 dpi. Gross pathologic changes in tissues such as trachea, lung and kidney were observed. Serum samples collected from chickens that survived were checked for IBV antibodies using a commercial enzyme-linked immunosorbent assay (ELISA) kit (IDEXX Laboratories, Westbrook, ME, USA). The endpoint titres were calculated according to the manufacturer's instructions, and titres of more than 396 were considered positive for IBV antibody. All experiments complied with the institutional animal care guidelines and were approved by the University of Yangzhou Animal Care Committee. A histopathological assay and haematoxylin-eosin (HE) staining of the tissues were performed as previously described (Zhang et al., 2013) . Briefly, the tracheae, lungs, and kidneys of the sick chickens were fixed in 10% neutral formalin, dehydrated in alcohol, and embedded in paraffin. They were then stained with HE and screened via light microscopy. 3.1. Sera against mass-and 4/91-type IBV could not completely neutralize isolate CK/CH/2010/JT-1 The last dilution of each serum against H120 and 4/91, which protected 50% of the embryos against 10 2 EID 50 s of strains M41 and CK/CH/2014/TM98, was tested at 1:9.85 and 1:25.46, respectively. However, serum against H120 could not neutralize 10 2 EID 50 s of CK/CH/2010/JT-1, while the end-point of serum against 4/91 that could neutralize 10 2 EID 50 s of CK/CH/2010/JT-1 was only 1:2.14. These results indicate that the CK/CH/2010/JT-1 isolate was antigenically distinct from the IBV Mass and 4/91 serotypes. By the third passage, typical symptoms, such as curling, dwarfing and stunting, were observed in the eggs. Inoculated egg death appeared at the 4th passage, and the inoculated eggs had all died by the 6th-10th passages. The virus was identified as IBV by RT-PCR and named CK/CH/2010/JT-1. The complete genome sequence of CK/CH/2010/JT-1 was obtained by assembling 25 overlapping sequences, with a 27670-nucleotide genome. The complete genome encodes several different genes, and components of the related genes are ordered 5 0 UTR -1a-1b-S-3a,b,c(E)-M-4b-5a,b-N-UTR3 0 . Sequence analysis showed that CK/CH/2010/ JT-1 shares the highest complete genome sequence identity with GX-YL9 (96.3%) and GX-YL5 (96.1%) but only 44.3%-89.2% S1 nucleic acid identity with 59 reference strains (Fig. 1) . BLASTN analysis showed that the S1 gene of CK/CH/2010/JT-1 shares 90%-99% similarity with 78 isolates from China in recent years. Phylogenetic analysis of the S1 gene of IBVs revealed that CK/CH/2010/JT-1 and 31 other isolates described in BLASTN analysis are grouped as a new cluster (Fig. 1) , which is separated from the previously identified serotypes and genotypes. To examine the sequence characteristics of isolate CK/CH/2010/ JT-1, potential recombination events were analysed. Three recombination events were identified in the genome of CK/CH/ 2010/JT-1. The assessment using RDP4 software revealed four recombination events in the CK/CH/2010/JT-1 genome ( Table 1) . The parent strains CQ04-1, SAIBK and GX-YL5 were grouped into genotype CK/CH/LSC/99I, while LX4 and ck/CH/LDL/091022 were QX-like genotypes, and 4/91 and partridge/GD/S14/2003 were genotypes 4/91 and tl/CH/LDT3/03 (Fig. 1) . Three recombinant events were further confirmed by Simplot ( Fig. 2A) and Bootscan analysis (Fig. 2B) . Recombinant gene fragments were considered recombinants if any crossover event appeared between two putative parental strains. Obvious recombination signals were found in simplot and bootscan analyses of recombinant events 1, 3 and 4. The similarities and breakpoints also agreed with the RDP software. However, the recombinant gene fragments shared low homology with the parental sequences, similarity analysis could not further confirm recombination in event 2. Furthermore, whole non-structural protein (nsp) 15; genes 5a, 5b and N; the 3 0 UTR and parts of the 5 0 UTR; nsp 14 and nsp 16; and genes S and M were involved in the recombination events. To further track the origins of the other parts of isolate CK/CH/ 2010/JT-1 0 s genome, we cut the complete genome sequence into different fragments according to recombination analysis and conducted a pairwise comparison of the fragment sequences of isolate CK/CH/2010/JT-1 with four QX-like, three CK/CH/LSC/99Itype, two tl/CH/LDT3/03-type and two 4/91-type strains. Genomic positions at 366-17159 of CK/CH/2010/JT-1 appeared very similar to CK/CH/LSC/99I-type strains GX-YL5 and GX-YL9, while genomic positions at 19812-21135 and 21171-23867 fell into the same group with tl/CH/LDT3/03-type IBVs. However, the CK/CH/2010/JT-1 isolate fell into the group of QX-like strains from location 23868-24708 in the phylogenetic tree (Fig. 3A) . These results strongly suggested that CK/CH/2010/JT-1 originated from homologous RNA recombination events from multiple template switches among QX-like, CK/CH/LSC/99I-, tl/CH/LDT3/ 03-and 4/91-type IBVs (Fig. 3B) , and the genome approximate positions 1-17159 and 24709-27670, 17160-19811 and 23868-24708, 19812-21135 and 21171-23867, and 21136-21770 possibly arose from CK/CH/LSC/99I-type, QX-like, tl/CH/LDT3/03-type and 4/91-type IBVs, respectively. Animal experiments for CK/CH/2010/JT-1 in chicken were performed to evaluate the pathogenicity of isolate CK/CH/2010/ Fig. 1 . Phylogenetic tree of the S1 gene from infectious bronchitis viruses (IBVs). A phylogenetic tree was constructed with the UPGMA method using MEGA version 5.05. Bootstrap values were determined from 1000 replicates of the original data. The CK/CH/2010/JT-1 strain is marked with a black solid triangle, and strains described in BLASTN analysis are marked with black hollow triangles. RDP, GENECONV, MaxChi, Chimaera, 3Seq (1.056 Â 10 À38 , 2.420 Â 10 À39 , 2.134 Â 10 À19 , 3.918 Â 10 À21 , and 1.623 Â 10 À54 ) a Only transferred gene fragments where at least 5 detection methods p-value 1 Â10 À14 are included in the table. b "Genes" indicates the coding sequences contained within the fragment introduced by recombination. c The "minor parent" is the sequence closely related to that from which sequences in the proposed recombinant region may have been derived. d The "major parent" is the sequence closely related to that from which the greater part of the recombinant's sequence may have been derived. JT-1. At 3 dpi, birds in group 1 inoculated with CK/CH/2010/JT-1 appeared depressed with ruffled feathers and huddled together. At 4 dpi, the birds in group 1 developed respiratory symptoms, most of the birds breathing with head shaking and a few breathing with open mouth and lifted up neck. Notably, infected birds began to die at 4 dpi, which continued until 13 dpi, with the death rate reaching 43.75% (Fig. 4A) . No deaths or clinical signs were observed in group 2. Among the seven dead birds in group 1, slight haemorrhage with serous catarrhal exudates could be seen in the trachea (Fig. 5B) . Typical kidney lesions were found in all dead chickens, the affected kidneys from the dead chickens showing predominant gross lesions that were pale, mottled, and swollen, with renal tubules and ureters that were distended with excess urate (Fig. 5D) . Clinical signs in the surviving birds tended to disappear gradually, and the rest of the pathological changes were present by 28 dpi. No gross lesions were observed in any bird in group 2. Antibody responses in birds that survived at 28 dpi were measured. All of the chickens showed positive reactions (Fig. 4B) , and the mean titre induced by the CK/CH/2010/JT-1 strain was 702.38 at 28 dpi. These data clearly demonstrated that strain CK/CH/2010/JT-1 isolated here was a virulent IBV. Microscopic examination of the tracheal tissues revealed extensive degeneration and necrosis of the ciliated epithelial cells, sometimes with pseudoacinar structures resulting from dropout of dead cells (Fig. 6B) . Lung lesions included haemorrhage, congestion, and lymphocytic infiltration in the alveolar lumen (Fig. 6D) . Severe renal lesions were characterized by degeneration and necrosis of renal tubular epithelial cells, lymphocytic infiltration in the interstitium, and exfoliated renal tubular epithelial cells, with erythrocytes being frequently observed (Fig. 6F) . Recently, various live-attenuated and inactivated vaccines have been widely and extensively used to prevent IB disease in chicken farms in China. However, the IB vaccine program has not been very effective, as new IBV serotypes, genotypes and variants are sometimes isolated from vaccinated chickens (Feng et al., 2015; Han et al., 2011) . This study identified an IBV, CK/CH/2010/JT-1, and 31 other isolates grouped into a novel genotypic cluster different from previously identified genotypes. Although 23 of the 31 isolates were reported in other five papers, we found that among the 23 isolates, 13 isolates were clustered into a genotype with reference strains A2 and QXIBV (Luo et al., 2012) , 9 isolates were grouped with reference strains 4/91 and UK 7/93 (Feng et al., 2014; Ji et al., 2011; Li et al., 2010) , and one isolate showed far evolutionary distances with reference tl/CH/LDT3/03 strains, but grouped into tl/CH/LDT3/03-type . In our study, the CK/CH/2010/JT-1 and 31 strains showed a low S1 nucleic acid identity with reference strains. Phylogenetic analysis of the S1 gene showed these strains could be grouped as a new genotypic cluster which was different from the previously identified genotypes IBV in China. Some vaccines, such as Mass-type H120 and 4/91-type 4/91 vaccine strains, could not provide effective protection against CK/CH/2010/JT-1 infection. In addition, the CK/ CH/2010/JT-1-like isolates form a new genotypic cluster, which is similar to a new IBV cluster in Korea (Lim et al., 2011) . Coronaviruses are identified in a wide variety of animals and divided into three distinct groups based on genotypic and serological characterization (Brian and Baric, 2005; Lai and Cavanagh, 1997) . Coronaviruses have a high frequency of recombination due to the unique mechanism of viral replication. IBV is a member of the genus Gamma Coronavirus. Novel IBV variants could emerge mainly via two approaches: recombination and mutation Toro et al., 2012) . In this paper, we found the CK/CH/2010/JT-1-like strains underwent three recombination events. The recombination fragments were located in approximate positions 17160-24708, 19812-23876 and 21136-21770 . The genome approximate positions 1-17159 and 24709-27670, 17160-19811 and 23868-24708, 19812-21135 and 21171-23867, and 21136-21770 of the CK/CH/2010/JT-1-like IBV possibly arose from CK/CH/LSC/99I-type, QX-like, tl/CH/LDT3/03-type and 4/91-type IBVs respectively (Fig. 3B ). However the chronological order of the recombination events needs further investigation. Recently, nine genotypes of IBV have been circulating in China. The predominant genotype is QX-like IBV, which has spread widely in Asia, Europe and Africa since it was first reported in 1996 in China (Cook et al., 2012; Knoetze et al., 2014; Toffan et al., 2011) . The second dominant genotype detected in China is the CK/CH/LSC/ 99I-type . The co-circulation of QX-like, CK/CH/ LSC/99I, vaccine, tl/CH/LDT3/03 and 4/91 genotype strains in the same chicken farm provided the possibility of new recombinants like CK/CH/2010/JT-1. Previous studies found that the recombination hot spots tended to be located in the 1a (nsp2 and nsp3), 1b (nsp16), S1, envelope, membrane genes and the 3 0 UTR . Most of the recombination breakpoints identified in isolate CK/CH/2010/JT-1 are located in these genes. To our knowledge, this recombinant virus is different from those in other reports Wu et al., 2016; Xu et al., 2016; Zhao et al., 2014) . All these data indicate that the new isolate is recombined with homologous RNA recombination events from multiple template switches among QX-like, CK/CH/LSC/99I-, tl/CH/LDT3/ 03-and 4/91-type IBV strains. A new genotype virus could acquire new biological characteristics, including virulence from gene recombination. The IBV spike and replicase genes are known to be determinants of pathogenicity (Armesto et al., 2009; Wickramasinghe et al., 2011) . The spike gene also plays a key role in IBV tissue tropism (Wickramasinghe et al., 2011) . The Mass-type M41 and 4/91-type 793/B could cause the respiratory signs without death in infected chickens (Benyeda et al., 2009) . Recently, the Chinese QX-like IBV virulent strains were found to cause 10%-50% mortality in infected chickens (Feng et al., 2015; Liu et al., 2009; Shi et al., 2011; Sun et al., 2011) , while the virulent CK/CH/LSC/99I-type strains YN and ck/CH/LDL/07I and the tl/CH/LDT3/03-type strain tl/CH/LDT3/03 caused 65%, 30%, and 80% mortality with severe kidney lesions in infected chickens, respectively (Feng et al., 2012 (Feng et al., , 2015 Liu et al., 2005) . The isolate CK/CH/2010/JT-1 could cause respiratory and severe renal disease and resulted in deaths of 43.75% experimental infections of 3-day- old SPF chickens. These results indicate the CK/CH/2010/JT-1 isolate is highly virulent. The replicase gene of isolate CK/CH/2010/ JT-1 was derived from QX-like, CK/CH/LSC/99I-and tl/CH/LDT3/03type strains, and the spike gene possibly originated from 4/91-and tl/CH/LDT3/03-type strains. All of these features yield 43.75% mortality with respiratory and severe renal disease in CK/CH/2010/ JT-1 infected chickens. The altered 4/91 fragment in the spike gene may not affect the tissue tropism of isolate CK/CH/2010/JT-1; the tissue tropism of the virus may have more than one pathogenicity factor, an issue to be studied in the future. In this study, an IBV, named CK/CH/2010/JT-1, was identified from an IBV-vaccinated flock in China. CK/CH/2010/JT-1-like viruses are grouped as a new genotypic cluster different from previously identified genotypes. Genome sequence analysis suggested that isolate CK/CH/2010/JT-1 originated from QX-like, CK/CH/LSC/99I, tl/CH/LDT3/03 and 4/91 genotypes. CK/CH/2010/ JT-1 of IBV is very virulent in chickens. Such novel recombinants should receive more attention during viral surveillance for better development of anti-IBV strategies. Comparison of the pathogenicity of QX-like, M41 and 793/B infectious bronchitis strains from different pathological conditions Coronavirus genome structure and replication Comparisons of envelope through 5B sequences of infectious bronchitis coronaviruses indicates recombination occurs in the envelope and membrane genes Coronavirus avian infectious bronchitisvirus The long view: 40 years of infectious bronchitis research Virulent avian infectious bronchitis virus, People's Republic of China Analysis of S1 gene of avian infectious bronchitis virus isolated in southern China during Development and efficacy of a novel live-attenuated QX-like nephropathogenic infectious bronchitis virus vaccine in China A 15-year analysis of molecular epidemiology of avian infectious bronchitis coronavirus in China Emergence of a group 3 coronavirus through recombination Molecular evolution and emergence of avian gammacoronaviruses Phylogenetic distribution and predominant genotype of the avian infectious bronchitis virus in China during Two genotypes of infectious bronchitis virus are responsible for serological variation in KwaZulu-Natal poultry flocks prior to 2012 Evolution of infectious bronchitis virus in Taiwan: positively selected sites in the nucleocapsid protein and their effects on RNA-binding activity The molecular biology of coronaviruses 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 Successful cross-protective efficacy induced by heatadapted live attenuated nephropathogenic infectious bronchitis virus derived from a natural recombinant strain Isolation of avian infectious bronchitis coronavirus from domestic peafowl (Pavo cristatus) and teal (Anas) Histopathologic analysis (haematoxylin-eosin stain) of chickens infected with IBV strain CK/CH/2010/JT-1 A, normal trachea tissue control B, trachea tissue infected with CK/CH/2010/JT-1 displays extensive dropout, degeneration, and necrosis of ciliated epithelial cells and lymphocytic infiltration (black arrow). Scale bar = 50 mm D, lung tissue infected with CK/CH/2010/JT-1 displays haemorrhage and congestion (black arrow) JT-1 displays severe renal lesions, including degeneration, necrosis of renal tubular epithelial cells, and exfoliated renal tubular epithelial cells (black arrow) Molecular characterization and pathogenicity of infectious bronchitis coronaviruses: complicated evolution and epidemiology in china caused by cocirculation of multiple types of infectious bronchitis coronaviruses Origin and characteristics of the recombinant novel avian infectious bronchitis coronavirus isolate ck/CH/LJL/111054 Phylogenetic analysis of the S1 glycoprotein gene of infectious bronchitis viruses isolated in China during Evaluation of recombinant fowlpox virus expressing infectious bronchitis virus S1 gene and chicken interferon-gamma gene for immune protection against heterologous strains Phylogenetic analysis of infectious bronchitis coronaviruses newly isolated in China, and pathogenicity and evaluation of protection induced by Massachusetts serotype H120 vaccine against QX-like strains Recombination in avian gamma-coronavirus infectious bronchitis virus QX-like infectious bronchitis virus in Africa Genetic diversity and selection regulates evolution of infectious bronchitis virus Binding of avian coronavirus spike proteins to host factors reflects virus tropism and pathogenicity Genome sequence and origin analyses of the recombinant novel IBV virulent isolate SAIBK2 Characterization and analysis of an infectious bronchitis virus strain isolated from southern China in 2013 Complete genome sequence of a recombinant nephropathogenic infectious bronchitis virus strain in China Comparison of hematoxylin-eosin staining and methyl violet staining for displaying ghost cells Molecular characterization of an infectious bronchitis virus strain isolated from northern China in 2012 This study was made possible by funding from the Priority Academic Program Development of Jiangsu Province Higher Education Institutions (PAPD). The authors declared that they have no conflict of interest.