key: cord-0936501-78m8ylz5
authors: Mondal, Shankar P; Cardona, Carol J
title: Comparison of four regions in the replicase gene of heterologous infectious bronchitis virus strains
date: 2004-06-20
journal: Virology
DOI: 10.1016/j.virol.2004.03.032
sha: e78f1c1d9cd697f8b9141342fa2bb93209942e83
doc_id: 936501
cord_uid: 78m8ylz5

Infectious bronchitis virus (IBV) produces six subgenomic (sg) mRNAs, each containing a 64 nucleotide (nt) leader sequence, derived from the 5′ end of the genome by a discontinuous process. Several putative functional domains such as a papain-like proteinase (PL(pro)), main protease (M(pro)), RNA-dependent RNA polymerase (RdRp), and RNA helicase encoded by the replicase gene are important for virus replication. We have sequenced four regions of the replicase genes corresponding to the 5′-terminal sequence, PL(pro), M(pro), and RdRp domains from 20 heterologous IBV strains, and compared them with previously published coronavirus sequences. All the coronavirus 5′-termini and PL(pro) domains were divergent, unlike the M(pro) and the RdRp domains that were highly conserved with 28% and 48% conserved residues, respectively. Among IBV strains, the 5′ untranslated region including the leader sequence was highly conserved (>94% identical); whereas, the N-terminal coding region and the PL(pro) domains were highly variable ranging from 84.6% to 100%, and 77.6% to 100% identity, respectively. The IBV M(pro) and RdRp domains were highly conserved with 82.7% and 92.7% conserved residues, respectively. The BJ strain was the most different from other IBVs in all four regions of the replicase. Phylogeny-based clustering based on replicase genes was identical to the antigen-based classification of coronaviruses into three groups. However, the IBV strain classification based on replicase gene domains did not correlate with that of the type-specific antigenic groups. The replicase gene sequences of many IBVs recovered from infected chickens were identical to those of vaccine viruses irrespective of serotype, suggesting that either there has been an exchange of genetic material among vaccine and field isolates or that there is a convergent evolution to a specific replicase genotype. There was no correlation between the genotype of any region of the replicase gene and pathotype, suggesting that the replicase is not the sole determinant of IBV pathogenicity.

Avian infectious bronchitis virus (IBV), the prototype species of the family Coronaviridae, causes an acute and highly contagious disease of chickens. Numerous IBV strains exist, and clinical signs of the disease they cause vary from mild to severe. More than 20 IBV serotypes have been recognized worldwide, which are defined by crossreactivity among the structural proteins (Cowen and Hitchner, 1975; Hopkins, 1974) . Live vaccines containing strains of IBV from multiple serotypes are routinely used in commercial chicken flocks. Disease outbreaks in commer-cial flocks occur when new viruses emerge from existing viruses by mutations, insertions, deletions, or recombination (Cavanagh et al., 1992; Gelb et al., 1991; Kusters et al., 1990) . However, because these conclusions are based on available sequence from structural genes, viral evolution has only been described for this small part of the genome. Because only two IBVs have been fully sequenced, it has been impossible to fully understand the viral evolutionary process including the extent to which vaccine viruses have contributed genetic material to pathogenic field strains.

Six subgenomic (sg) mRNA species are transcribed from the IBV genome in virus-infected cells. This includes a genome-length mRNA (mRNA 1) of 27.6 kb and five sg species (mRNA 2 -6) with sizes ranging from 2 to 7 kb. These mRNAs form a 3V-co-terminal nested structure (Stern and Sefton, 1984) , each containing a short nontranslated leader sequence (64 nucleotides (nt)), derived from the 5V end of the genome by a discontinuous process. The phenomenon of leader switching has been described for several coronaviruses including IBV (Chang et al., 1994; Makino et al., 1988; Stirrups et al., 2000) . The mRNA 1 contains two large overlapping open reading frames (ORF 1a and 1b) encoding 441-and 300-kDa polyproteins, respectively . By a unique frameshifting mechanism, the downstream ORF 1b is produced as a fusion protein of 741 kDa with 1a (Brierley et al., 1987 (Brierley et al., , 1989 . These polyproteins are proteolytically processed into smaller products required for RNA synthesis and other aspects of viral replication. Four main structural proteins are encoded by sgmRNAs: the spike (S) glycoprotein (mRNA 2), the small envelope (E) protein (mRNA 3), the membrane (M) glycoprotein (mRNA 4), and the nucleocapsid (N) protein (mRNA 6) (Lai, 1990; Spaan et al., 1988; Sutou et al., 1988) . The S1 fragment of the spike protein carries antigenic epitopes that induce virus-neutralizing (VN) antibody, and also determines virus serotype .

Two overlapping papain-like proteinase (PL pro ) domains encoded between nucleotides (nt) 4243 and 5553 (amino acids (aa) 1239-1675), and a picornavirus 3C-like proteinase (3C-LP) domain encoded between nt 8937 and 9357 (aa 2804 -2943) , are involved in the production of mature viral proteins required for replication Liu et al., 1994 Liu et al., , 1997 . Proteolytic processing of the polyproteins by equivalent protease activities has been demonstrated in murine and human coronaviruses (reviewed in Lim et al., 2000; Ng and Liu, 2000) . Despite the conservation of function, there are some organizational differences in PL pro domains between IBV and some other coronaviruses (Herold et al., 1993; Lee et al., 1991) . Because of its key role in replicase gene expression, 3C-LP has been termed the coronavirus main protease (M pro ; Ziebuhr et al., 2000) . The IBV PL pro domain-1 (PL1 pro ), encoded between nt 4243 and 5019 (aa 1239-1497), is responsible for cleavage of the N-terminal regions of 1a and 1a/1b polyproteins in IBV-infected cells, whereas the PL pro domain-2 (PL2 pro ) encoded between nt 4681 -5553 (aa 1385 -1675) is reported to be functionally inactive in IBV Lim et al., 2000; . A more important role in the processing of polyproteins is played by the M pro , which is responsible for the cleavage of other regions of the ORF 1a and 1a/1b polyproteins (Liu et al., 1994 (Liu et al., , 1997 Liu, 1998, 2002) . The catalytic center of both proteinases consists of His and Cys residues in all coronaviruses (reviewed in Lim et al., 2000; Ng and Liu, 2000) . The processed products constitute a group of replicative enzymes, including the RNA-dependent RNA polymerase (RdRp) and RNA helicase. Therefore, the significance of the coronavirus leader sequence and several putative functional domains essential for viral replication made them logical choices for sequence analysis and phylogenetic comparisons to correlate with the serological or pathological classification of coronaviruses.

Phylogeny-based clustering based on replicase genes is identical to the original antigen-based classification of coronaviruses (Chouljenko et al., 2001; Hegyi and Ziebuhr, 2002; Stephensen et al., 1999) . However, within the group 3 coronaviruses, no comparisons of replicase-based and antigen-based phylogenies have been made. One report has described the alignment of the 5V-proximal region of replicase genes from a limited number of IBV strains (Stirrups et al., 2000) , and a few reports describe IBV M pro and RdRp domains from the Beaudette strain only (Hegyi and Ziebuhr, 2002; Stephensen et al., 1999) . Here, we have sequenced four regions of the replicase gene corresponding to the 5Vterminal end, PL pro , M pro , and RdRp domains from 20 heterologous IBV strains, and compared then with previously published IBV and other coronavirus sequences. We present evidence that (i) the M pro and RdRp regions are highly conserved among all coronaviruses unlike the 5Vterminal sequence and PL pro domains, (ii) the clustering of heterologous IBV strains based on the replicase gene sequences does not correlate with the antigen-based S1 phylogeny, (iii) several widely used IBV vaccines and field strains isolated over the past decade have closely related replicase genes suggestive of a possible common ancestry or a converging evolution, and (iv) the replicase gene is unlikely to be the sole determinant of IBV pathogenicity.

The 5V -terminal sequences (nt 1-1263)

The nt sequence alignment data obtained from the IBV 5V-termini are presented in Table 1 , Figs. 1 and 2a. The leader sequences (nt 1 -64) of H52, CU705, CU510, Ark99, ArkDPI, CU805, CU570, GA98, Cal99, CU994, K1699, Gray, and BJ strains are identical to each other and each have 3 nt substitutions compared to the Beaudette strain (Fig. 1 ). The DE072 virus has four, the vaccine and challenge M41 strains, and CU-T2 have three, Ma5, K0751, and Florida have two, and Conn has a single nt substitution compared to the Beaudette strain (Fig. 1) . The leader junction sequences termed as leader transcriptionregulating sequences (TRSs) are absolutely conserved in all of the IBV strains compared. The sequence distance data (Table 1 ) obtained from the whole 5V UTR (nt 1 -528) indicate that nt sequence identities were 94 -100% among all IBV strains in this region. The N-terminal coding region (nt 529-1263) of the 5V terminal sequence is more variable with the lowest identity represented by the BJ strain at 84.6-86% identity with the other strains. Pairwise comparisons of the N-terminal sequences from other strains revealed 89.7-100% identity irrespective of serotype, geographic source, and virulence (Table 1 ). The challenge strain NV M41 is identical to the vaccine strain M41 in the noncoding region but in the coding region has a single nt change that results in an amino acid change (data not shown). These relationships are graphically shown in the phylogenetic tree created from the nt sequences of the entire 5V terminal region (Fig. 2a) . No detectable sequence similarity was found when IBV 5V terminal sequences were compared to those of coronaviruses from other groups (data not shown).

The PL pro domains (nt 4243 -5553: aa 1239-1675)

Sequence distance data derived from multiple alignments of the two overlapping PL pro domains of IBV strains revealed 77.6 -100% identity in PL1 pro and 79 -100% in Fig. 1 . Alignment of the first 80 nucleotides from cDNA sequences derived from the 5V end of 22 heterologous IBV strains. The gray shading represents differences between sequences. The leader sequence is denoted by a solid black line above the sequence and the leader junction sequence is boxed. Sequences with >95% identity are indicated in bold letters. a Because of the overall sequence identity of the vaccine and challenge strains of M41, for brevity, only the vaccine strain is included in this table. PL2 pro , with the former being more variable than the latter ( Table 2 ). The BJ strain differed the most from other strains by 18.1 -22.4% in PL1 pro and 16.5-21% in PL2 pro . The catalytic Cys 1274 and His 1437 residues in the PL pro domains are conserved in all IBV strains (Fig. 3 ). In the overlapping area of the PL pro domains, a region of 24 residues with more heterogeneity has been identified, where only two residues (Phe 1470 and Lys 1485 ) are conserved among all IBV strains (Fig. 3) . In this heterogeneous region, the vaccine and challenge M41 and BJ strains each differed by 12 aa (50%), Conn, CU510, Florida, CU805, and Gray each differed by 11 aa (45%), CU705, Ark99, and Cal99 strains each differed by 10 aa (40%), and all other strains are identical to each other and differed by 9 aa (37.5%) when compared to the Beaudette strain sequence. In the complete PL pro domain, the NV M41 challenge strain has 1 nt and 1 aa difference from the M41 vaccine strain. The relationships between the different sequences are shown in the phylogenetic tree derived from analysis of the aa sequences of the entire PL pro region (aa 1239 -1675) (Fig. 2b) . No detectable sequence similarity was found when IBV PL pro sequences were compared to the PL pro sequences of coronaviruses from other groups (data not shown).

The M pro (aa ORF1a 2804 -2943) and the RdRp (aa ORF1b 548 -780) domains

The alignments of the deduced aa sequences of IBV M pro and RdRp domains are summarized in Table 3 . The overall sequences are highly conserved in which there are 82.7% (116 aa) absolutely conserved residues in M pro and 92.7% (216 aa) absolutely conserved residues in RdRp (data not shown). The remaining 24 aa differences found in M pro and 17 aa differences in RdRp are distributed throughout the domains (Table 3 ). The M pro and RdRp domains of ArkDPI, GA98, and Gray strains are identical to each other and those of the challenge NV M41 and vaccine M41 strains are also identical to each other. These relationships are graphically shown in the phylogenetic tree created from the aa sequences of the M pro and RdRp domains (Figs. 2c and d).

By aligning with the corresponding regions of prototype members of other coronavirus groups, all M pro had 28% (40 aa) absolutely conserved residues including catalytic Cys-His dyad and three substrate-binding sites (Tyr-Met-His), and the RdRp had 47.6% (111 aa) absolutely conserved residues (data not shown). The four RdRp motifs (DXXXD, SGXXXTXXXN, SDD, and K) appear to be fully conserved Sequences with >95% identity are indicated in bold letters. a Because of the overall sequence identity of the vaccine and challenge strains of M41, for brevity, only the vaccine strain is included in this table. 2829 2831 2835 2840 2845 2847 2848 2849 2851 2855 2860 2874 2881 2894 2898 2900 2906 2925 2928 2929 2935 2937 569 572 592 596 606 608 611 674 678 681 698 708 744 752 754 756 765 Beau

. 

The findings of the present study can be summarized as: (i) the M pro and RdRp regions are highly conserved among all coronaviruses unlike the 5V-terminal sequence and PL pro domains, (ii) the clustering of heterologous IBV strains based on the replicase gene sequences does not correlate with the antigen-based S1 phylogeny, (iii) several widely used IBV vaccines and field strains isolated over the past decade have closely related replicase genes suggestive of a possible common ancestry or a converging evolution, and (iv) the replicase gene is unlikely to be the sole determinant of IBV pathogenicity.

Coronaviruses fall into three antigenic groups (Dea et al., 1990) , which are defined by cross-reactivity among the structural proteins. Unlike the structural proteins, the coronavirus replicase gene is not subjected to immune selective pressure. Therefore, one would expect that the replicase gene would be highly conserved among coronavirus groups. In fact, this is the case when the M pro and RdRp sequences are compared (Table 3 ) and as others have reported for other coronavirus groups (Chouljenko et al., 2001; Gonzalez et al., 2003; Hegyi and Ziebuhr, 2002; Stephensen et al., 1999) . We found a high percentage of sequence similarities in coronavirus M pro (with 28% absolutely conserved residues), including the substrate recognition site, and an even higher similarity in the RdRp (with 48% absolutely conserved residues), including the polymerase motifs. However, two-way sequence comparisons of the coronavirus 5V-termini and the PL pro domains of ORF1a had very low sequence identity between coronaviruses of different groups as has been reported by others (Herold et al., 1999) . Despite the fact that these regions have varying degrees of genetic diversity, the relationships between the coronaviruses of different groups are the same, no matter what region of the genome is compared. When recombination occurs (Cavanagh and Davis, 1988) , it would seem to occur within the major coronaviruses antigenic groups, but not between them. This is expected because in order for recombination to occur, two viruses must be present in the same cell, which is unlikely between coronaviruses adapted to different host species. The highly conserved M pro -mediated processing pathways as well as the substrate binding sites of this enzyme in all coronaviruses (Anand et al., 2003; Hegyi and Ziebuhr, 2002) make this proteinase an attractive target for the development of drugs directed against coronaviruses. The active domains of RdRp also maintain a high degree of homology among all coronaviruses, as has been found in other RNA viruses (Koonin and Dolja, 1993; Otsuka et al., 1999) , making this proteinase another potential target for drugs controlling coronaviruses.

With IBV, the phylogenies of the four regions of the replicase gene did not correlate with those of the S1 gene, which correspond to the type-specific antigenic groups (Fig. 2) . For example, unlike the S1 genes, the replicase genes of the Ma5 and H52 Mass serotype viruses consis-tently clustered with those of non-Mass strains (Figs. 2ad) . The common clustering of replicase genes of these heterologous IBV strains could be due to an introduction of a replicase gene from a common source through introduction. Others have speculated that IBV strains substitute large genomic fragments in multiple genes (Lee and Jackwood, 2000) because the viability of the progeny viruses depends on the specific interactions of more than one gene working in concert to maintain structural and replicative integrity. Alternatively, it is probably more likely that the smaller S1 genes are exchanged between viruses through recombination and the stable replicase backbone is maintained. It is entirely possible that recombination events have occurred to introduce new replicase or new S1 genes into these viruses.

Leader acquisition in mRNA synthesis appears to be TRS dependent, and the leader-to-body joining is guided by a base-pairing interaction involving leader and body TRSs (Hiscox et al., 1995; Pasternak et al., 2001) . The phenomenon of leader switching involves a recombination event over a region comprising the leader sequence and the adjacent 5V UTR (Chang et al., 1996) . Previous reports have suggested that the leader sequences of IBV strains may contain higher (17.2%) substitution rates when compared to the complete 5V UTR (4.3%; Stirrups et al., 2000) . In our study, the 5V UTR of all 22 IBV strains are highly conserved (94 -100% identity), indicating the exchange of leader or adjacent 5V UTR fragments has occurred in mixed infections. This observation is supported by experimental evidence of the exchange of leader sequences among heterologous IBV strains (Stirrups et al., 2000) and in mixed MHV infections (Makino et al., 1986) . The sequence comparison scores are more variable in the N-terminal coding region (Table 1) , and the most in PL pro domains (Table 2) including a more variable region in the latter (Fig.  3 ) than the IBV M pro and RdRp domains ( Table 3 ). The reason for sequence variation in the N-terminal half of ORF1a among IBV strains is not clear.

The genetic variations in the 5V terminal region and PL pro gene suggest that there is some selective pressure on these regions. We speculate that there are some replicase genotypes that may confer a selective advantage to the virus. We examined both the 5V terminal region and the PL pro gene by year and location of isolation, but excluded the M pro and RdRp genes because they are so highly conserved that no clear differences between isolates can be discerned. In the 5V terminal region and the PL pro gene, the IBV strains isolated from commercial chickens over the past decade cluster with three specific vaccine virus strains, ArkDPI, H52, and Ma5, no matter where in the United States they were isolated. There are three exceptions to this clustering, the 5V terminal region of the Gray virus isolated in 1960 (Winterfield and Hitchner, 1962) clusters with modern field strains (Fig. 2a) , although its PL pro gene is more like the Conn and Florida reference viruses isolated in 1956 and 1971, respectively (Jungherr et al., 1956; Winterfield et al., 1971) (Fig. 2b ). In addition, there are two modern field isolates, one from China (BJ) and one from California (K1699) isolated in 2003 and 2001, respectively, that are outliers and do not appear to have 5V termini or PL pro genes that are related to other viruses (Figs. 2a and b) . This common clustering of three widely used vaccines and IBV field isolates from the past decade suggests to us that there may be a common evolutionary progression in the 5V terminal region and in the PL pro gene. Although the 5V terminal region genotype found in modern isolates was present in 1960, when the Gray virus was isolated, it appears to be more prevalent today. The question arises as to why this type of progression to a specific genotype would have occurred and how it could have happened. As to the question of why, we speculate that the more modern genotype must offer some fitness advantage to the virus. As to how it could have occurred, we consider the widespread use of multiple types of modified live vaccines in both the broiler and egg layer commercial poultry industries as the possible mechanism by which this modern replicase gene genotype was introduced. The Ark, Conn, and Mass (Ma5, H52) strains have been used simultaneously as live vaccines among commercial poultry (Cavanagh and Naqi, 2003) . Although we cannot prove or disprove that recombination has occurred, we can point to the evolutionary trends in the S1 gene as compared to the evolutionary trends in the replicase gene. The trends in S1 evolution have been toward ever increasing diversity (Jia et al., 2002) , while it appears that the replicase gene is evolving to less variety between isolates. These juxtaposed evolutionary trends suggest that (1) these two regions are under very different selective pressures and thus undergo much different rates and modes of change, and (2) the exchange of genetic material may have occurred between vaccine strains and field isolates.

The genetic basis of IBV pathogenicity is not known. The role of the replicase gene in the virulence of IBV has not been established, and it is likely that some changes in the M pro or viral polymerase would influence viral replication rate, and thus the pathogenic potential of a virus. In some RNA viruses, replicase genes are known to contain the determinants of species tropism and virulence factors (Brandt et al., 2001; Hatta et al., 2001; Yao et al., 2001) . In our study, there were no differences in any essential residues in M pro and RdRp domains. In addition, we found that there was no clear relationship between viral pathotype and replicase genotype. As an example, the H52 vaccine strain is very closely related to the virulent DE072 isolate in all regions of the replicase (99.8% identity in the 5V terminus, 99.5% in PL pro , 99.3% in M pro , 99.6% in RdRp) despite their differences in pathotype. In addition, a direct comparison of the replicase genes of the M41 challenge and vaccine strains where there were very few aa differences between the two strains with clearly different pathotypes. None of the aa changes between the two M41 strains were consistently observed in either all virulent or all avirulent strains. These findings would support our contention that the genotype of the replicase is not the sole determinant of pathogenicity among strains of IBV.

Twenty IBV strains from nine different serotypes, Mass (M41, NV M41, Ma5, H52, K0751, CU705), Conn (Conn, CU510), Florida, CU-T2, Arkansas (Ark-99, ArkDPI), Delaware 072 (DE072, CU805, CU570), Georgia 98 (GA98), California 99 (Cal99, CU994, K1699), and Gray, were used in this study (Table 4 ). These virus strains were selected because they are very common vaccine and field strains that are frequently isolated from commercial chickens and are associated with clinical disease. All the reference IBV strains were obtained from Dr. Syed Naqi's laboratory (Cornell University, Ithaca), where they had been propagated in specific pathogen-free (SPF) embryonated chicken eggs to maintain the stock. In our laboratory, all of the viruses were also propagated in 9-day-old embryonated SPF eggs (Charles River SPAFAS, Wilmington, MA). Four regions of the replicase gene from 20 heterologous IBV strains, which correspond to the 5V-termini, and the PL pro , the M pro , and the RdRp domains (Fig. 4) , were sequenced and compared with those of previously published IBV (Beaudette, BJ) and other coronavirus sequences.

Genomic RNA was extracted from virus-inoculated allantoic fluid with TRIzol reagent (Invitrogen Co., Carlsbad, CA) following the manufacturer's instructions. The first strand cDNA synthesis and subsequent PCR were performed using GeneAmp RNA PCR kit (Perkin-Elmer Cetus, Norwalk, CT) following the manufacturer's instructions. The selected regions of the replicase genes were amplified with primers synthesized from conserved regions designed by aligning IBV Beaudette sequence and sequences from other coronaviruses. The PCR profiles involved an initial denaturation for 4 min at 94 jC followed by 35 cycles of annealing at 55 jC for 30 s, extension at 72 jC for 1 min, and melting at 94 jC for 1 min. The amplified products were directly sequenced at a commercial sequencing facility (Davis Sequencing, Davis, CA). The genome sense Pol32 (5V-CACTAGCCTTGCGC-TAGA-3V) and anti-sense Pol1501 (5V-GACCAACC-TTCTGGTTCAAC-3V) primers were used to amplify the 5V-termini that correspond to nt 32 -1501 of the Beaudette sequence . The 5V ends of all IBV strains were amplified with FirstChoice RLM-RACE cDNA amplification kit (Ambion Inc., Austin, TX) using genomic antisense Pol1263 (5V-GCTTGCAAGACAAGTTCCTGC-3V) primer according to the manufacturer's protocol. The Pol1263 paired with another inner site primer Pol133 (5V-GGCACCTGGCCACCTGTCAGG-3V) selected from conserved areas (by aligning 11 IBV sequences including the Beaudette virus) could amplify nt 133 -1263 (according to Beaudette sequence) of all IBV strains tested in our laboratory (data not shown). Three more consensus primer pairs were used to amplify other regions as follows: primers for the overlapping PL pro amplification are PL55 (5V-AGGATAAA-GAAATYCTCTTC-3V) and PL58 (5V-GGACCACAYAAA-GAACCCTC-3V); primers for M pro amplification are Mp51 (5V-CGCCACGTTACTCTATTGGT-3V) and Mp33 (5V-GCCGCATAGAGCCATGCTAC-3V); and primers for RdRp amplification are Rp52 (5V-CTACTATGACTAATAGG-CAG-3V) and Rp31 (5V-CTGAGAAAGCTCTTGATAGAG-3V).

Assembly of contiguous sequences, translation of nt sequence into protein sequence, and initial multiple sequence alignments were performed with Vector NTI Suite 9 software (InforMax, North Bethesda, MD). Selected sequences from GenBank were included in the alignment. After the alignment, the AlignX program constructs a phylogenetic tree employing the neighbor-joining method. Comparisons with published sequences were made by performing BlastN search (GenBank). Pairwise blast searches were also performed when there were no significant hits with BlastN search.

The accession numbers for IBV replicase gene sequences are as follows: (a) the 5V-termini: M41, AY392047; NV M41, AY561719; Ma5, AY561720; H52, AY392048; K0751, AY561721; CU705, AY561722; Conn, AY392049; CU510, AY561723; Florida, AY392050; CU-T2, AY561724; Ark99, AY392051; ArkDPI, AY392052; DE072, AY392054; CU805, AY561725; CU570, AY561726; GA98, AY392053; Cal99, AY392055; CU994, AY561727; K1699, Fig. 4 . Diagram of the coronavirus replicase gene (encoding ORF 1a and 1b) illustrating the 5V untranslated region (UTR), papain-like proteinase (PL pro ), main protease (M pro ), RNA-dependent RNA polymerase (RdRp), metal-binding (MB), and RNA helicase (RH) domains. The arrows originating from PL pro and M pro indicate the cleavage sites. The four regions amplified in this study are indicated as shading areas with arrow at each end (bottom). AY319651 a Except for the Gray strain that is nephrotropic and the BJ strain whose tissue tropism is unknown, all other strains are pneumotropic. b Based on S1 gene sequence. c Extensively propagated in vitro since the first isolation (Beaudette and Hudson, 1937 ) and known to be attenuated. d Introduced to the US market at about 1990 (Ma5) and 1980 (H52). e IBV recovered from infected chicken flocks and associated with clinical disease and their serotype determined by S1 sequence.

BoCV (NC _ 003045), and SARS (NC _ 004718) were obtained from GenBank. The S1 sequences of IBV strains were

Coronavirus main proteinase (3CLpro) structure: basis for design of anti-SARS drugs

Cultivation of the virus of infectious bronchitis

Completion of the sequence of the genome of the coronavirus avian infectious bronchitis virus

Molecular determinants of virulence, cell tropism, and pathogenic phenotype of infectious bursal disease virus

An efficient ribosomal frame-shifting signal in the polymerase-encoding region of the coronavirus IBV

Characterization of an efficient coronavirus ribosomal frameshifting signal: requirement for an RNA pseudoknot

Coronavirus IBV: removal of spike glycopolypeptide S1 by urea abolishes infectivity and haemagglutination but not attachment to cells

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

Infectious bronchitis

Coronavirus IBV glycopolypeptides: locational studies using proteases and saponin, a membrane permeabilizer

Location of the amino acid differences in the S1 spike glycoprotein subunit of closely related serotypes of infectious bronchitis virus

A cis-acting function for the coronavirus leader in defective interfering RNA replication

The UCUAAAC promoter motif is not required for high-frequency leader recombination in bovine coronavirus defective interfering RNA

Comparison of genomic and predicted amino acid sequences of respiratory and enteric bovine coronaviruses isolated from the same animal with fatal shipping pneumonia

Serotyping of avian infectious bronchitis viruses by the virus-neutralization test

Antigenic and genomic relationships among turkey and bovine enteric coronaviruses

Serologic and cross-protection studies with several infectious bronchitis virus isolates from Delmarava-reared broiler chickens

Variant serotypes of infectious bronchitis virus isolated from commercial layer and broiler chickens

Antigenic and S-1 genomic characterization of the Delaware variant serotype of infectious bronchitis virus

A comparative sequence analysis to revise the current taxonomy of the family Coronaviridae

Molecular basis for high virulence of Hong Kong H5N1 influenza A viruses

Conservation of substrate specificities among coronavirus main proteases

Nucleotide sequence of the human coronavirus 229E RNA polymerase locus

A human RNA viral cysteine proteinase that depends upon a unique Zn2+-binding finger connecting the two domains of a papain-like fold

Investigation of the control of coronavirus subgenomic mRNA transcription by using T7-generated negative-sense RNA transcripts

Serological comparisons of strains of infectious bronchitis virus using plaque purified isolants

A novel variant of avian infectious bronchitis virus resulting from recombination among three different strains

Genetic and antigenic diversity in avian infectious bronchitis virus isolates of the 1940s

A new serotype of infectious bronchitis virus responsible for respiratory disease in Arkansas broiler flocks

Immunological differences in strains of infectious bronchitis virus

Evolution and taxonomy of positive-strand RNA viruses: implications of comparative analysis of amino acid sequences

Sequence evidence for RNA recombination in field isolates of avian coronavirus infectious bronchitis virus

Background paper. Transcription and replication of coronavirus RNA: a 1989 update

Evidence of genetic diversity generated by recombination among avian coronavirus IBV

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

The complete sequence (22 kilobases) of murine coronavirus gene 1 encoding the putative proteases and RNA polymerase

Characterization of the two overlapping papain-like proteinase domains encoded in gene 1 of the coronavirus infectious bronchitis virus and determination of the C-terminal cleavage site of an 87-kDa protein

Identification of a novel cleavage activity of the first papain-like proteinase domain encoded by open reading frame 1a of the coronavirus Avian infectious bronchitis virus and characterization of the cleavage products

Characterisation and mutational analysis of an ORF 1a-encoding proteinase domain responsible for proteolytic processing of the infectious bronchitis virus 1a/1b polyprotein

A 100-kilodalton polypeptide encoded by open reading frame (ORF) 1b of the coronavirus infectious bronchitis virus is processed by ORF 1a products

Identification, expression, and processing of an 87-kDa polypeptide encoded by ORF 1a of the coronavirus infectious bronchitis virus

Proteolytic processing of the coronavirus infectious bronchitis virus 1a polyprotein: identification of a 10-kilodalton polypeptide and determination of its cleavage sites

Regulation of mRNA 1 expression by the 5V-untranslated region (5V-UTR) of the coronavirus infectious bronchitis virus (IBV)

Leader sequences of murine coronavirus mRNAs can be freely reassorted: evidence for the role of free leader RNA in transcription

Defective-interfering particles of murine coronavirus: mechanism of synthesis of defective viral RNAs

Evaluation of commercially produced infectious bronchitis virus vaccines against an IBV field isolate obtained from broilers in California

Isolation and characterization of a novel antigenic subtype of infectious bronchitis virus serotype DE072

Identification of a 24-kDa polypeptide processed from the coronavirus infectious bronchitis virus 1a polyprotein by the 3C-like proteinase and determination of its cleavage sites

Further characterization of the coronavirus infectious bronchitis virus 3C-like proteinase and determination of a new cleavage site

Membrane association and dimerization of a cysteine-rich, 16-kilodalton polypeptide released from the C-terminal region of the coronavirus infectious bronchitis virus 1a polyprotein

Similarity relations of DNA and RNA polymerases investigated by the principal component analysis of amino acid sequences

Sequence requirements for RNA strand transfer during nidovirus discontinuous subgenomic RNA synthesis

Coronaviruses: structure and genome expression

Phylogenetic analysis of a highly conserved region of the polymerase gene from 11 coronaviruses and development of a consensus polymerase chain reaction assay

Coronavirus multiplication: locations of genes for virion proteins on the avian infectious bronchitis virus genome

Leader switching occurs during the rescue of defective RNAs by heterologous strains of the coronavirus infectious bronchitis virus

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

Etiology of an infectious nephritis -nephrosis syndrome of chickens

Characteristics of an isolate of infectious bronchitis virus from chickens in Florida

Sequences in influenza A virus PB2 protein that determine productive infection for an avian influenza virus in mouse and human cell lines

Virus-encoded proteinases and proteolytic processing in the Nidovirales

We thank Drs. Syed A. Naqi and Udeni Balasuriya for their important suggestions in preparing the manuscript. We also thank Ms. Kelly Weaver for her excellent technical assistance.