key: cord-279495-zxerb7de authors: Liu, Xiaoli; Shao, Yuhao; Ma, Huijie; Sun, Chuyang; Zhang, Xiaonan; Li, Chengren; Han, Zongxi; Yan, Baolong; Kong, Xiangang; Liu, Shengwang title: Comparative analysis of four Massachusetts type infectious bronchitis coronavirus genomes reveals a novel Massachusetts type strain and evidence of natural recombination in the genome date: 2012-11-21 journal: Infect Genet Evol DOI: 10.1016/j.meegid.2012.09.016 sha: doc_id: 279495 cord_uid: zxerb7de Four Massachusetts-type (Mass-type) strains of infectious bronchitis coronavirus (IBV) were compared genetically with the pathogenic M41 and H120 vaccine strains using the complete genomic sequences. The results revealed that strains ck/CH/LNM/091017 and ck/CH/LDL/101212 were closely related to the H120 vaccine, which suggests that they might represent re-isolations of vaccine strains or variants of vaccine strains that have resulted from the accumulated point mutations after several passages in chickens. In contrast, strains ck/CH/LHLJ/07VII and ck/CH/LHLJ/100902 had a close genetic relationship with the pathogenic M41 strain. In addition, molecular markers have been identified that distinguish between field and vaccine (or vaccine-like) Mass-type viruses, which may be able to differentiate between field and vaccine strains for diagnostic purposes. Phylogenetic analysis, and pairwise comparison of full-length genomes and the nine genes, identified the occurrence of recombination events in the genome of strain CK/VH/LHLJ/07VII, which suggests that this virus originated from recombination events between M41- and H120-like strains at the switch site located at the 3′ end of the nucleocapsid (N) genes. To our knowledge, this is the first time that evidence for the evolution and natural recombination under field conditions between Mass-type pathogenic and vaccinal IBV strains has been documented. These findings provide insights into the emergence and evolution of the Mass-type IB coronaviruses and may help to explain the emergence of Mass-type IBV in chicken flocks all over the world. In 1931, Schalk and Hawn described ''an apparently new respiratory disease of chicks'' in North Dakota in the United States, which was considered to be infectious bronchitis (IB) by later researchers of avian respiratory diseases (Schalk and Hawn, 1931) . Currently, IB still occurs in nearly all poultry-producing countries; it is a highly contagious, acute, and economically important viral disease of chickens. The etiology of IB, which was first demonstrated by Beach and Schalm (1936) , is infectious bronchitis virus (IBV). IBV is grouped in the genus Gammacoronavirus of the family Coronaviridae in the order Nidovirales (de Groot et al., 2011) . The coronavirus genomes are the largest among the known RNA viruses and are polycistronic, generating a nested set of subgenomic RNAs with common 5 0 and 3 0 sequences (Masters, 2006) . Like those of all other coronaviruses, the 5 0 two-thirds of the IBV genome consists of two large replicase open reading frames (ORFs), ORF1a and ORF1b. The ORF1a polyprotein (pp1a) can be extended with ORF1b-encoded sequences via a À1 ribosomal frameshift at a conserved slippery site (Brierley et al., 1987) , which generates the polyprotein pp1ab, comprising more than 7000 amino acids, which includes the putative RNA-dependent RNA polymerase (RdRp) and RNA helicase (HEL) activity (Ziebuhr et al., 2001) . The pp1a and pp1ab of IBV are processed autocatalytically by two different viral proteases, encoded by a papain-like protease (PLP) and a 3C-loke protease (3CL pro ) (Lee et al., 1991; Ziebuhr et al., 2000 Ziebuhr et al., , 2001 . Other putative domains, presumably associated with a 3 0 -to-5 0 exonuclease (ExoN) activity, a poly(U)-specific endo-RNase (XendoU) activity, and a 2 0 -O-methyltransferase (2 0 -O-MT) activity, have been predicted in pp1ab (Ivanov et al., 2004; Snijder et al., 2003) . The 3 0 end of a coronavirus genome includes the viral structural and accessory protein genes: a spike (S) glycoprotein gene; an envelope (E) protein gene; a membrane (M) glycoprotein gene; a nucleocapsid (N) phosphoprotein gene; and several ORFs that encode putative non-structural accessory proteins (Masters, 2006) . Of the virus-encoded proteins, the S1 subunit of the S protein carries virus-neutralizing activity, determines the serotype of IBV and is responsible for viral attachment to cells. It is also a major determinant of cell tropism in culture (Casais et al., 2003) . The accumulation of point mutations, deletions, insertions and recombination events that have been observed in multiple structural genes, especially the S1 gene, of IBV recovered from naturally occurring infections have been considered to contribute to the genetic diversity and evolution of IBV, and consequently, to a number of IBV serotypes (Cavanagh, 2007) . The occurrence and emergence of multiple serotypes of the virus have complicated control by vaccination because many serotypes and variants do not confer complete cross-protection against each other (Cavanagh and Gelb, 2008) . The originally discovered Massachusetts (Mass) type of IBV was identified in the United States, beginning in the 1950s (Fabricant, 2000; Johnson and Marquardt, 1975; Mondal et al., 2001) . Mass-type strains have been isolated in Europe and Asia since the 1950s and up to the present day (Cavanagh and Gelb, 2008) , together with dozens of other serotypes that have been isolated in Africa, Asia, India, Australia, Europe, and South America (Cavanagh, 2001 (Cavanagh, , 2003 (Cavanagh, , 2005 . The first Mass-type ''H'' vaccines were developed in about 1960. They include H120 and H52 (Bijlenga et al., 2004) , and are used very commonly and widely around the world. However, virus of this type is occasionally isolated from Massachusetts-vaccinated and -unvaccinated flocks with respiratory clinical signs. Some of the viruses have shown close genetic relationships with pathogenic Mass-type, rather than vaccine, strains by S1 gene analysis. However, conclusions based on the genetic analysis of a single gene sequence, and sometimes even a partial gene sequence, require caution because the true phylogeny can only be demonstrated by analyzing complete genomic sequences. Herein, we sequenced the complete genome of four IBV Mass-type strains that showed S1 gene diversity (Liu et al., 2009; Ma et al., 2012; Sun et al., 2011) , and we present evidence for in-field recombination between pathogenic and vaccinal strains. Furthermore, we characterized the molecular variability of the four Mass-type strains to gain insight into the emergence and evolution of these viruses. Four Mass-type IBV strains were used for complete genomic sequence comparison and analysis in this study. Strain ck/CH/LHLJ/ 07VII was isolated in 2007 from the kidney of a layer hen vaccinated with H120 and 4/91 in Heilongjiang province, China (Liu et al., 2009) . Strain ck/CH/LNM/091017 was isolated in 2009 from the swollen proventricular tissues of a broiler vaccinated with H120 in Neimenggu province, China . Strains ck/ CH/LDL/101212 and ck/CH/LHLJ/100902, both of which were isolated in 2010, were obtained from laying hens in Dalian and Heilongjiang provinces, respectively, in China; the birds were suffering from nephropathogenic lesions and respiratory signs, respectively. In addition, the diseased birds in both flocks were suffering from proventriculitis (Ma et al., 2012) . All of the IBV strains have been associated with various IB outbreaks in recent years in China and were assigned to the Mass-type strains by S1 sequence analysis. To avoid the possible mutation in the viral genome after serial passages in specific-pathogen-free (SPF) embryonated chicken eggs, the first passage of each original virus stock was used and purified once by propagating in 9-to 11-day-old SPF chicken eggs with a dose of l0 4 -fold dilutions per egg, and the presence of viral particles in the allantoic fluids of inoculated eggs was confirmed with a negative contrast electronic microscope (JEM-1200, EX) and reverse transcriptase-polymerase chain reaction (RT-PCR) as described previously . In addition, since these viruses were isolated from chickens vaccinated with H120, it is possible that mixed IBV infections are present in one chicken flock. To exclude this, nine clones of S1 gene of each virus obtained from three independent PCR reactions were sequenced and compared. Sequences of each virus identical to the previously results were obtained (Liu et al., 2009; Ma et al., 2012; Sun et al., 2011) . Fertile White Leghorn SPF chicken eggs were obtained from the Laboratory Animal Center, Harbin Veterinary Research Institute, the Chinese Academy of Agricultural Sciences, China. To determine the full-length genomic sequences of the four viruses, 15 pairs of overlapping primers encompassing the entire genome were used. The primers were designed in regions that are conserved among most of the IBV strains available in the Gen-Bank database. The sequences and locations of the primers used in RT and PCR in this study are presented in Table 1 . Viral RNA was extracted from 200 ll of infectious allantoic fluid using TRIzol reagents (Invitrogen, Grand Island, USA), following the manufacturer's protocol. Complementary DNA (cDNA) was synthesized using 80 ll of the first strand mixture (Invitrogen) containing 20 lM of primers N (À), 0.5 mM each of dNTP (TaKaRa, Dalian, China) and 40 ll of total RNA. The mixture was incubated at 70°C for 5 min and then quick-chilled on ice for 2 min. The RT master mix was composed of 16 ll 5Â RT buffer (Invitrogen), 8 ll 10 mM DTT, 200 U of M-MLV Reverse Transcriptase (Invitrogen), and 20 U RNAse inhibitors (Invitrogen). This RT master mix was incubated at 37°C for 3 h. The reaction was terminated by heating at 70°C for 10 min then chilling on ice for 5 min. The PCR was performed in a 25 ll reaction containing 2 ll first strand cDNA; 15 nmol each of downstream and upstream primers; 5 ll of 10Â PCR buffer (Mg 2+ Plus, TaKaRa); 4 ll of 2.5 mmol dNTPs; 2 U Taq polymerase (TaKaRa); and 18 ll of water. The reaction was conducted at 95°C for 5 min, and 30 cycles of 94°C for 1 min; 50°C for 1 min; 72°C for 2 min, and a final extension step of 72°C for 10 min. A product, detectable by ethidium bromide staining, of the expected size was generated. 2.4. The 5 0 -and 3 0 -ends of the genome A cDNA clone representing the 5 0 and 3 0 ends of the genome of the four IBV strains was synthesized according to the 5 0 RACE and 3 0 RACE System for rapid amplification of cDNA ends (TaKaRa). PCR was performed according to the instructions accompanying the kits. The sense and antisense primers used to amplifying the 5 0and 3 0 -ends of the genome had been designed on the basis of the sequences obtained above that were constant in the four IBV strains, respectively. The outer and inner primers used to amplify the 5 0 -end of the four IBV strains were 5 0 -CAGCTATGGCAATGCG CAG-3 0 and 5 0 -CATCTTTGGTGTCTCA/TCC-3 0 , respectively. The primer used to amplify the 3 0 -end was 5 0 -GAGGAGAGGAACAATGC ACA-3 0 . The DNA generated by PCR amplification was cloned using a Ttailed vector, pMD18-T (TaKaRa), and transformed using JM109 competent cells (TaKaRa) according to the manufacturer's instructions. At least five clones of each fragment in each strain were sequenced and the consensus sequence was determined. The sequences were analyzed using the Sequencher 4.5 sequence analysis program, and a single contiguous sequence comprising the entire IBV genome of each of the four IBV strains was constructed. The nucleotide and amino acid sequences of the entire genome of the four IBV strains were assembled, aligned, and compared with those of other reference IBV and Turkey coronavirus (TCoV) strains using the MEGALIGN program in DNAStar (version 7, Lasergene Corp, Madison, WI). The ORFs were determined using the Gene Runner program version 3.00 (http://www.generunner.com) by comparison with those of other reference IBV and TCoV strains. A total of 39 IBV and 7 TCoV reference strains, for which entire genomic sequences were available in GenBank database, were selected for phylogenetic analysis of full-length genomes. The selected avian coronavirus reference strains and their accession numbers are provided in Table 2 . Phylogenetic analysis, accurate estimation and comparison of the 5 0 -UTR, Gene 1, S1, S2, Gene 3, M, Gene 5, N and 3 0 -UTR of the four IBV strains was conducted with those of the Mass-type strains selected in this study using the Clustal V method of DNAStar software and MEGA4 (Liu et al., 2009) , and the alignments were edited manually and adjusted to remove mistakes. Deletion, insertion and gene recombination were determined according to the results of the phylogenetic analysis and pairwise comparisons. The full genomic sequences of the four Mass-type IBV strains described in this report have been deposited in the GenBank database with accession numbers ck/CH/LNM/091017 JF330899, CK/ CH/LHLJ/07VII JF274479, ck/CH/LDL/101212 JF828981 and CK/ CH/LHLJ/100902 JF828980. Four Mass-type IBV strains were subjected to genome sequencing and phylogenetic analysis in this study. The sequences of each the four strains were assembled into one contiguous sequence to represent the entire viral genomes. Sequences of 27630, 27678, 27630 and 27473 nucleotides were obtained from strains ck/CH/ LNM/091017, ck/CH/LHLJ/07VII, ck/CH/LDL/101212 and ck/CH/ LHLJ/100902, respectively, excluding the polyadenylation tail at the 3 0 end. The genomes of the viruses were similar overall in their coding capacity and genomic organization to those of other IBVs. The genome of each of the viruses contained two large slightly overlapping ORFs in the 5 0 two-thirds of the genome and multiple additional ORFs in the 3 0 one-third of the genome. Both termini were flanked with untranslated regions (UTRs). Ten ORFs were identified within the genome. Gene 1 contained motifs common to all coronaviruses, including ribosomal frameshifting and slippery sequences, because ORF1b is translated in the À1 frame. The typical coronavirus structural genes encoding the S, E, M and N proteins were identified following Gene 1 (Fig. 1) . The genome organization was determined to be as follows: 5 0 -UTR-Gene 1 (ORF1a, 1b)-S-Gene 3 (ORFs 3a, 3b, E)-M-Gene 5 (ORFs 5a, 5b)-N-UTR-3 0 . The analysis of the complete genome showed that strains ck/CH/ LNM/091017 and ck/CH/LDL/101212 possessed 99.9% and 99.8% nucleotide identity with H120, respectively. However, they shared 91.3% and 91.4% identity with M41, respectively. Phylogenetic analysis using the full-length genome and the 5 0 -UTR, Gene 1, S1, S2, Gene 3, M, Gene 5, N and 3 0 -UTR showed that the IBV strains ck/ CH/LNM/091017 and ck/CH/LDL/101212 consistently formed the same clade with vaccine-related strains of Mass-type (Figs 2 and 3 ). The analysis of the S1 gene showed that strains ck/CH/LNM/ 091017 and ck/CH/LDL/101212 had high nucleotide identities (99.8% and 99.8%, respectively) with H120, while they had 97.7% and 97.7% identity with M41. Multiple alignments revealed that there were 2 and 1 nucleotide mutations within the S1 gene between strains ck/CH/LNM/091017 and ck/CH/LDL/101212 and H120; however, there were 36 and 39 mutations between strains ck/CH/LNM/091017 and ck/CH/LDL/101212 and M41. All these results suggest that strains ck/CH/LNM/091017 and ck/CH/LDL/ 101212 are closely related to the H120 vaccine strain. The percent nucleotide similarity between strain ck/CH/LHLJ/ 100902 and H120 for the full-length genomes was 91.3%; however, the percent similarity was up to 99.7% between strain ck/CH/LHLJ/ 100902 and M41. In addition, in all of the nine trees constructed for the 5 0 -UTR, Gene 1, S1, S2, Gene 3, M, ORF5, N and 3 0 -UTR, the strain ck/CH/LHLJ/100902 constantly fell into the same clusters as the pathogenic M41 strain, and both belonged to the Mass-type. Pairwise comparison of the S1 protein gene revealed that strain ck/ CH/LHLJ/100902 had 3 and 41 nucleotide mutations with respect to M41 and H120, respectively. Taken together, these results demonstrate that strain ck/CH/LHLJ/100902 exhibits a close genetic relationship to pathogenic M41. Pathogenic and non-pathogenic Mass-type strains were clustered into different clades by phylogenetic analysis of full-length genomic sequences and the 5 0 -UTR, Gene 1, S1, S2, Gene 3, M, Gene 5, N and 3 0 -UTR. In addition, insertions and deletions were also observed that distinguished between the genomes of pathogenic and non-pathogenic Mass-type strains, as illustrated in Fig. 4 and Supplementary material 1. In non-pathogenic strains, five deletions: of 31 nucleotides, 5 nucleotides, 9 nucleotides, 10 nucleotides, and 9 nucleotides, respectively, were observed to be located in the nsp3 of Gene 1. They were found to occur between genomic positions 2897-2929, 3219-3225, 3241-3251, 3257-3268, and 3398-3407, respectively , by comparing the sequences with the homologous regions of pathogenic strains. In contrast, a 3-nucleotide and a 9nucleotide insertion were found in nsp6 and between the M gene and gene 5, respectively. Additionally, a cluster of insertions was found at the 3 0 -UTR region in the non-pathogenic strains. These changes might not only account, at least partly, for viral fitness when the pathogenic virus has become adapted to egg embryos Hewson et al., 2012) , but may act also as molecular markers, able to differentiate between vaccine and field strains, for diagnostic purposes. Comparative sequence analysis based on full-length genomic sequences and the sequences of the 5 0 -UTR, Gene 1, S1, S2, Gene 3, M, Gene 5 and N showed that strain CK/CH/LHLJ/07VII clustered with pathogenic Mass-type strains. The exceptions were the trees constructed using the 3 0 -UTR and S1 gene, in which CK/CH/LHLJ/07VII was grouped with non-pathogenic strains; this suggests that a possible recombination event may have occurred. Thus the N and 3 0 -UTR of CK/CH/LHLJ/07VII were carefully compared pairwise with those of strains ck/CH/LNM/091017, ck/CH/LDL/101212, ck/CH/ LNM/091017, H120 and M41. Parallel to the result of the phylogenetic analysis, ck/CH/LHLJ/07VII showed high similarity with M41 at the 5 0 -end of the N gene; however, it showed high similarity with vaccine strain H120 at the 3 0 -end of the N gene (Supplementary material 2). The data strongly suggest that CK/CH/LHLJ/07VII arose from a homologous RNA recombinant event that involved a template switch between Massachusetts pathogenic M41-like and non-pathogenic H120-like strains. We located the switch site at the 3 0 -end of the N gene (Supplementary material 2), which implies that the template switch occurred within the N gene. The percent nucleotide similarity between strain CK/CH/LHLJ/ 07VII and H120, and CK/CH/LHLJ/07VII and M41, for the full-length genomes was 92.5% and 99.4%, respectively. Alignment revealed that a 3-nucleotide insertion was located in nsp6 of the CK/CH/ LHLJ/07VII strain between genomic positions 24009 and 24013 (Supplementary material 3). In addition, the S1 gene of strain CK/ CH/LHLJ/07VII showed extensive mutations by pairwise comparison (Supplementary material 4) though it was grouped with H120 by S1 gene phylogenic analysis. These and our previous results (Liu et al., 2009) showed that, with the exception of the occurrence of recombination events, CK/CH/LHLJ/07VII has experienced multiple mutations and deletions in the genome over time. Understanding the evolution of Mass-type IBV is important because not only is this virus circulating worldwide but information on virus genomics will aid our understanding of the evolution and emergence of IBV with infectious potential in vaccinated chicken flocks. In this study, we focused on the full-length genomic sequences of four IBV isolates which had been shown to be of the Mass-type by S1 gene analysis (Liu et al., 2009; Ma et al., 2012; Sun et al., 2011) . Based on the high degree of similarity in the full genomic sequence, it could be concluded that two IBV strains, ck/ CH/LNM/091017 and ck/CH/LDL/101212, were very similar to the vaccine strain H120. They might therefore represent re-isolations of vaccine strains, although they were isolated from vaccinated chickens with respiratory disease. Similarly, IBV strains that showed a very close relationship to the H120 vaccine strain have been isolated from unvaccinated broiler flocks in Slovenia with respiratory problems (Krapež et al., 2010) . Alternatively, these strains might be variants of vaccine strains that have resulted from accumulated point mutations after several passages in chickens. A few key mutations in the S1 subunit of the spike protein might result in a change to a new serotype, which is defined as a lack of cross-neutralization with specific sera against different IBV serotypes (Cavanagh et al., 1992) . The point mutations found in the genome that distinguish between the two isolates and vaccine strain H120 might be the result of adaptive evolution driven by the host immune response when the vaccine strain was transmitted among chickens. Adaptive evolution is the process by which genetic changes in the viral genome leading to a more fit virus population become fixed over time, and it has been reported to occur in many coronaviruses (Hasoksuz et al., 2007; Lee and Jackwood, 2001; Shi et al., 2006; Tang et al., 2009; Zhang et al., 2006) . As shown in this study, and as also occurs in other countries (Dolz et al., 2008; Rimondi et al., 2009; Roussan et al., 2009) , the isolation of Mass-type IBV is expected, because attenuated vaccine strains are used extensively in chicken flocks in China. However, vaccination is not likely to be the only explanation for the circulation of Mass-type virus, because ck/CH/LHLJ/100902 and CK/CH/ LHLJ/07VII were most closely related to a Massachusetts pathogenic type strain, M41. The isolation of a Massachusetts pathogenic strain from H120-vaccinated chicken flocks may be due to vaccination failure in these flocks . Alternatively, molecular studies have shown that only a few changes in the amino acid composition of the S1 spike protein can result in immune failure, even when the majority of the virus genome remains unchanged (Cavanagh et al., 1992) . Our findings showed that mutations had occurred in the genomes of both ck/CH/LHLJ/100902 and CK/CH/ LHLJ/07VII, especially in the S1 genes of CK/CH/LHLJ/07VII though it was in the same group with H120 strain in the phylogenetic analysis, implicating that strain CK/CH/LHLJ/07VII has experienced evolution over time. It has been reported that amino acid changes may result from immunological pressure caused by the widespread use of vaccines (Cavanagh et al., 1988 . The occurrence of recombination events is another process that allows new strains to emerge, and this has been well documented in IBV (Hughes, 2011; Jia et al., 1995; Kottier et al., 1995; Kusters et al., 1990; Wang et al., 1993) and other coronaviruses Baric, 1992, 1994; Makino et al., 1986) . It is believed that the conditions for recombination amongst IBV strains in the field are as follows: an extremely large number of chickens, most maintained at high density; the ease of spread of the virus; and serotype cocirculation, including proof of co-infection with more than one serotype in a given flock (Cavanagh, 2007) . In China, intensive chicken farms are concentrated in many provinces, including Heilongjiang, where ck/CH/LHLJ/07VII was isolated (Liu et al., 2009) . Almost all the chickens in China receive Mass-type vaccines at a very young age and subsequently receive this vaccine a couple more times during the rest of their life span. Therefore the vaccine virus exists constantly in chickens; it may persist in various internal organs for 163 days or longer (Cavanagh and Gelb, 2008) . Generally, vaccination using the H120 vaccine provides full protection against pathogenic Mass-type pathogenic IBVs and prevents the same type of pathogenic strain from being replicated and spreading in the flocks. However, a single amino acid substitution at position 63 of the S1 subunit of the spike has resulted in escape mutants of Mass 41 (Cavanagh et al., 1988) . This may have occurred in the case of the ck/CH/LHLJ/07VII S1 gene (Liu et al., 2009) , which might have made it possible for both pathogenic and vaccine strains to co-exist in a given flock, leading to the occurrence of recombination. Similarly, an escape mutant could be a result of adaptive evolution driven by the host immune response. Consequently, it is likely that genetic changes due to adaptive evolution and recombination both contributed to the origin and evolution of strain ck/CH/LHLJ/07VII: it is possible that adaptive evolution created a mutant, followed by recombination between Mass 41-and H120-like strains to create a novel virus. The recombination events from which the CK/CH/LHLJ/07VII virus resulted can be explained by a scenario in which the recombination may have involved two parental viral strains, with initiation of RNA replication in a M41-like template of either negative or positive polarity (Liao and Lai, 1992) . This would be followed by switching of the polymerase-nascent cRNA complex to an H120like virus template. The switch may have occurred at the 3 0 -end of the N gene. In general, for a recombinant virus to emerge and establish itself in the field, it must be viable and have selective advantages. It has been reported that uptake of canine coronavirus (CCV) sequences by type II feline coronavirus (FCoV) may have led to increased viral fitness when compared with type I FCoV (Herrewegh et al., 1998) . Recombination can also result in the emergence of new strains with distinct characteristics, such as pathogenicity and tissue tropism (Worobey and Holmes, 1999) . In addition, in the CoV genome, as with most RNA viruses, the 3 0 -UTRs usually harbor important structural elements that are involved in replication and/or translation (Chang et al., 1994; Raman et al., 2003; Raman and Brian, 2005; Goebel et al., 2006; Züst et al., 2008) . In IBV, the 3 0 -UTR-binds to the N protein, which is essential for the synthesis of negative-strain viral RNA. Perhaps the acquisition of the 3 0 end of the N gene and the 3 0 -UTR from H120-like virus by an M41-like virus (e.g. CK/CH/LHLJ/07VII) can alter the efficiency of viral replication. This alteration may in turn affect pathogenicity. However, it remains unknown whether this is the true origin of CK/ CH/LHLJ/07VII, and therefore this strain is of particular importance to the surveillance of IBV in China. It will be of equal importance to examine future outbreaks of IBV in chickens by full-length genomic sequence analysis in the context of novel recombination events among IBV strains. Furthermore, investigations using reverse genetic systems might provide further insight into these issues and increase our understanding of IBV pathogenesis. A filterable virus, distinct from that of laryngotracheitis, the cause of a respiratory disease of chicks Development and use of the H strain of avian infectious bronchitis virus from the Netherlands as a vaccine: a review An efficient ribosomal frame-shifting signal in the polymeraseencoding region of the coronavirus IBV Recombinant avian infectious bronchitis virus expressing a heterologous spike gene demonstrates that the spike protein is a determinant of cell tropism Commentary: a nomenclature for avian coronavirus isolates and the question of species status Severe acute respiratory syndrome vaccine development: experiences of vaccination against avian infectious bronchitis coronavirus Coronaviruses in poultry and other birds Coronavirus avian infectious bronchitis virus Location of the amino acid differences in the S1 spike glycoprotein subunit of closely related serotypes of infectious bronchitis virus Amino acids within hypervariable region 1 of avian coronavirus IBV (Massachusetts serotype) spike glycoprotein are associated with neutralization epitopes Infectious bronchitis Variations 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 A cis-acting function for the coronavirus leader in defective interfering RNA replication Virus Taxonomy, Classification and Nomenclature of Viruses, Ninth Report of the International Committee on Taxonomy of Viruses, International Union of Microbiological Societies, Virology Division Molecular epidemiology and evolution of avian infectious bronchitis virus in Spain over a fourteen-year period The early history of infectious bronchitis Evidence for variable rates of recombination in the MHV genome Map locations of mouse hepatitis virus temperature sensitive mutants: confirmation of variable rates of recombination A hypervariable region within the 3 0 cis-acting element of the murine coronavirus genome is nonessential for RNA synthesis but affects pathogenesis A 15-year analysis of molecular epidemiology of avian infectious bronchitis coronavirus in China Biologic, antigenic, and full-length genomic characterization of a bovine-like coronavirus isolated from a giraffe Feline coronavirus type II strains 79-1683 and 79-1146 originate from a double recombination between feline coronavirus type I and canine coronavirus The present of viral subpopulations in an infectious bronchitis virus vaccine with differing pathogenicity -a preliminary study Recombinational histories of avian infectious bronchitis virus and turkey coronavirus Multiple enzymatic activities associated with severe acute respiratory syndrome coronavirus helicase A novel variant of infectious bronchitis virus resulting from recombination among three different strains The neutralizing characteristics of strains of infectious bronchitis virus as measured by the constant virus variable serum methods in chicken tracheal cultures Experimental evidence of recombination in coronavirus infectious bronchitis virus Molecular analysis of infectious bronchitis viruses isolated in Slovenia between 1990 and 2005: a retrospective study Sequence evidence for RNA recombination in field isolates of avian coronavirus infectious bronchitis virus Origin and evolution of Georgia 98 (GA98), a new serotype of avian infectious bronchitis virus The complete sequence (22 kilobases) of murine coronavirus gene 1 encoding the putative proteases and RNA polymerase RNA recombination in a coronavirus: recombination between viral genomic RNA and transfected RNA fragments Molecular characterization and pathogenicity of infectious bronchitis coronaviruses: complicated evolution and epidemiology in China caused by cocirculation of multiple types of infectious bronchitis coronaviruses Genetic diversity of avian infectious bronchitis coronavirus in recent years in China High-frequency RNA recombination of murine coronaviruses The molecular biology of coronaviruses Isolation and characterization of a novel antigenic subtype of infectious bronchitis virus serotype DE072 Stem-loop III in the 5 0 untranslated region is a cis-acting element in bovine coronavirus defective interfering RNA replication Stem-loop IV in the 5 0 untranslated region is a cisacting element in bovine coronavirus defective interfering RNA replication Molecular characterization of avian infectious bronchitis virus strains from outbreaks in Argentina Infectious bronchitis virus in Jordanian chickens: seroprevalence and detection An apparently new respiratory disease of chicks Evolutionary implications of avian infectious bronchitis virus (AIBV) analysis Unique and conserved features of genome and proteonome of SARS-coronavirus, an early split-off from the coronavirus group 2 lineage Phylogenetic analysis of infectious bronchitis coronaviruses newly isolated in China, and pathogenicity and evaluation of protection induced by Mass serotype H120 vaccine against strains of the LX4-type (QX) Differential stepwise evolution of SARS coronavirus functional proteins in different host species Evidence of natural recombination within the S1 gene of infectious bronchitis virus Evolutionary aspects of recombination in RNA viruses Adaptive evolution of the spike gene of SARS coronavirus: changes in positively selected sites in different epidemic groups Virus-encoded proteinases and proteolytic processing in the Nidovirales The autocatalytic release of a putative RNA virus transcription factor from its polyprotein precursor involves two paralogous papain-like proteases that cleave the same peptide bond Genetic interactions between an essential 3 0 cis-acting RNA pseudoknot, replicase gene products, and the extreme 3 0 end of the mouse coronavirus genome Supplementary data associated with this article can be found, in the online version, at http://dx.doi.org/10.1016/j.meegid.2012. 09.016.