key: cord-330260-xuw31zfn authors: Chen, Hui-Wen; Huang, Yuan-Pin; Wang, Ching-Ho title: Identification of Taiwan and China-like recombinant avian infectious bronchitis viruses in Taiwan date: 2009-01-20 journal: Virus Res DOI: 10.1016/j.virusres.2008.11.012 sha: doc_id: 330260 cord_uid: xuw31zfn Infectious bronchitis virus (IBV) infections in poultry cause great economic losses to the poultry industry worldwide. The emergence of viral variants complicates disease control. The IBV strains in Taiwan were clustered into two groups, Taiwan group I and Taiwan group II, based on the S1 gene. A variant was previously identified and showed a distinct S1 gene homology with other local strains. This study investigated the 3′ 7.3 kb genome of eight Taiwan strains isolated from 1992 to 2007. The genes of interest were directly sequenced. Sequence analyses were performed to detect any recombination event among IBVs. The results demonstrated that all of the examined viruses maintained the typical IBV genome organization as 5′-S-3a-3b-E-M-5a-5b-N-UTR-3′. In the phylogenetic analyses, various genes from one strain were clustered into separate groups. Moreover, frequent recombination events were identified in the Simplot analyses among the Taiwan and China CK/CH/LDL/97I-type strains. Putative crossover sites were located in the S1, S2, 3b, M genes and the intergenic region between the M and 5a genes. All of the recombinants showed chimeric IBV genome arrangements originated from Taiwan and China-like parental strains. Field IBVs in Taiwan undergo genetic recombination and evolution. Avian infectious bronchitis virus (IBV), with a large infectious RNA genome (27.6 kb), belongs to the family Coronaviridae, the group 3 Coronavirus. The RNA molecule of the IBV is linear, singlestranded, positive sense and possesses a 5 cap and a 3 poly-A tail. IBVs can replicate in respiratory, alimentary and urogenital tracts in chickens of all ages, resulting in great economic losses to the poultry industry worldwide (Cavanagh and Naqi, 2003) . Clinical signs associated with respiratory and enteric tracts and renal damage might be found in infected chickens (Cavanagh, 2007) . During infection, IBV possesses a unique discontinuous transcription system, i.e. a nested set of 3 co-terminal subgenomic mRNAs sharing a common leader sequence in the 5 end are transcribed in the presence of negative-stranded RNA intermediates. The viral polymerase "jumping" possibly contributes to the high RNA recombination frequency in coronaviruses (Lai, 1992; Lai and Holmes, 2001) . The IBV genome encodes four structural proteins: spike glycoprotein (S), envelope protein (E; also known as sM), membrane glycoprotein (M), and nucleocapsid protein (N) (Cavanagh and Naqi, 2003) . The spike glycoprotein is post-translationally cleaved into S1 * Corresponding author. Tel.: +886 2 23690628; fax: +886 2 23631542. E-mail address: chingho@ntu.edu.tw (C.-H. Wang) . and S2 subunits. The S1 subunit anchors onto the viral outer membrane using the S2 subunit to form a club-shaped projection on the mature virion (Cavanagh, 1983) . The S1 glycoprotein is involved in cell attachment and carries epitopes for serotype-specific hemagglutination-inhibition and virus-neutralization antibodies (Cavanagh and Davis, 1986; Hodgson et al., 2004; Ignjatovic and Sapats, 2005; Koch et al., 1990) . The S2 glycoprotein, in which two antigenic determinants were identified (Koch et al., 1990) , may possess different secondary structures that affect the S1 specific antibody binding (Callison et al., 1999) . The membrane glycoprotein is associated with virus assembling and budding (Lai and Holmes, 2001) . The nucleocapsid protein interacts with genomic RNA to form the viral nucleocapsid, playing a role in viral RNA synthesis and cell immunity (Lai and Holmes, 2001) . In addition, four non-structural proteins of unknown function are expressed by the polycistronic genes, gene 3 and 5 . Virus mutants carrying truncated 3b genes demonstrate increased virulence and growth advantages in vitro and in ovo (Shen et al., 2003) . The open reading frames (ORFs) of 3a, 3b, 5a, and 5b encode accessory proteins not essential for IBV replication (Casais et al., 2005; Hodgson et al., 2006; Youn et al., 2005) . As a signature of avian coronaviruses, IBV strains are continuously evolving through point mutations and recombination of their genomes. Those variants have better adaptation or increased virulence advantageous to IB outbreaks. To date, a large number of IBV sero-or genotypes have been identified worldwide (Cavanagh, 2007 1171/92 1992 Broiler 3 Taoyuan Nephropathogenic TW-I 2296/95 1995 Broiler 2 Taoyuan Nephropathogenic TW-II 2575/98 1998 Broiler 4 Changhua Nephropathogenic TW-I 2992/ spike glycoprotein gene. (Bochkov et al., 2006; Dolz et al., 2006; Dolz et al., 2008; Jackwood et al., 2007) . It was reported that slight sequence differences in the S1 gene probably lead to poor crossprotection (Cavanagh et al., 1997) . Viruses of different types can co-circulate within a region (Capua et al., 1999; Liu et al., 2006) , raising the inter-strain RNA recombination frequency (Bochkov et al., 2007; Jia et al., 1995; Lee and Jackwood, 2000) . The wide use of live virus vaccine may also critically contribute to the genetic evolution of IBVs by acting as a heterologous RNA donor template (Kusters et al., 1990; Wang et al., 1993) . The emergence of viral variants has complicated disease control requiring persistent IBV molecular surveys. IBV strains in Taiwan were previously clustered into two groups, Taiwan group I (TW-I) and Taiwan group II (TW-II), on the basis of the S1 gene (Wang and Tsai, 1996) . However, a variant isolated in 2002 showed an unusually high S1 gene homology with China strains, but not in the N gene, suggesting an inter-strain recombination event (Huang et al., 2004) . In this study, the 3 7.3 kb genomes from Taiwan IBV strains were investigated to elucidate the genetic diversity of viruses. Eight IBVs isolated in Taiwan from 1992 to 2007 (Huang et al., 2004) and the vaccine strain H120 (ABIC Biological Laboratories Teva Ltd., Israel) were recovered for this study. The case histories of local strains are listed in Table 1 . Viruses were propagated in the allantoic cavity of 9-11-day-old specific pathogen free embryonated eggs (Animal Health Research Institute, Council of Agriculture, Tamsui, Taiwan). Each egg received 0.1-0.2 ml inoculation. After 48-72 h incubation, allantoic fluid was collected and frozen at −80 • C until use. Viral RNA was extracted from 200 L of virus-infected allantoic fluid using a viral nucleic acid extraction kit (Geneaid Biotech Ltd., Taipei, Taiwan) following the manufacturer's protocol. Previously published primers (Huang and Wang, 2007) were employed in this study to amplify the gene fragments. For sequencing the 3 untranslated region (UTR) of H120, one additional forward primer was designed as 5 -GGAACAATGCACAGCTGGAA-3 from strain H120-GD (GenBank accession no. AY028296). Reverse transcriptase-polymerase chain reaction (RT-PCR) was performed with one step in a reaction volume of 50 L containing 0.5 L of 5 U/L RealTaq DNA polymerase (Real Biotech, Taipei, Taiwan DNA sequencing in both strands from separate RT-PCR products was conducted by a commercial service (Mission Biotech). Each nucleotide was determined from at least four identical results. The obtained nucleotide sequences were compiled and the amino acid sequences were deduced using DNAStar software (DNAStar, Madison, WI). The IBV reference strains were retrieved from the GenBank database with the accession number listed in Table 2 . The AlignX program of the Vector NTI Suite 8 software Table 2 IBV reference strains included in this study. Country Accession numbers Armidale Australia DQ490205 (S1-N) a Vic Australia DQ490221 (S1-N) BJ China AY319651 (S1-N) Taiwan DQ646406 (S1-N) 2296/95 Taiwan DQ646404 (S1-N) 2575/98 Taiwan DQ646405 (S1-N) Beaudette USA NC 001451 (S1-N) Cal99 USA AY514485 (S1-N) CU-T2 USA U49858 (S1-N) DE072 USA AF274435 (S1), AY024337 (S2) (InforMax, North Bethesda, ML) was used to generate multiple sequence alignments and determine the nucleotide identity. Phylogenetic analyses were performed with the neighbor-joining method using MEGA version 4 (Tamura et al., 2007) . The bootstrap values were determined from 1000 replicates of the original data. The consecutive IBV nucleotide sequences from the S to N genes (6.8 kb) based on the multiple alignment results were introduced into similarity plots with SimPlot version 3.5.1 (Lole et al., 1999) . The nucleotide identity was calculated using the Kimura 2-parameter method with a transition-transversion ratio of 2 in each window of 500 bp. The window was successively furthered along the alignment using a 20-bp increment. At least four sequences were required to initiate an analysis. The IBV sequences resulting from this study were submitted to the GenBank database. The accession numbers are EU822336 (strain 3468/07), EU822337 (strain 3374/05), EU822338 (strain 3263/04), EU822339 (strain 3071/03), EU822340 (strain 2992/02), and EU822341 (strain H120). Fig. 1 . Phylogenetic analyses of the Taiwan strains (᭹), H120 ( ) and reference strains for structural and non-structural protein genes of IBVs. The phylogenetic trees were constructed using the MEGA version 4 by the neighbor-joining method (bootstrapping for 1000 replicates with its value >70%). Sequences that covered the 3 7.3 kb genome were determined from the Taiwan IBVs and strain H120. All of the examined viruses were found to maintain the typical IBV genome organization as 5 -S-3a-3b-E-M-5a-5b-N-3 (data not shown). The S gene size ranged from a minimum of 3471 nucleotides (strain 3071/03) to a maximum of 3501 nucleotides (strains 2992/02 and 3374/05). Compared with other local strains, single base mutations in the 3 end of the S and 3b genes from the strain 3071/03 changed the genetic code from GAA to TAA (Glutamine → stop codon), resulting in 27-and 3-base truncated ORFs, respectively. Similarly, in the 3374/05, a single base T-insertion in the 3b gene created an early stop codon, leading to a 48-base truncation. However, the 3a, 5b, and N gene ORF sizes were conserved among IBVs. All of the virus genomes carried an intergenic (IG) region located between the M and 5a genes with a size of 351-362 nucleotides. The IG region of the strain H120 was 55 bases longer than that of the strain Beaudette. The 3 UTR, downstream of the stop codon in the N gene, a region of 475 and 500 nucleotides was sequenced from local strains and H120, respectively. All of the local strains shared a high sequence identity (96-100%) in this region. However, only 52% and 73% identity were observed between H120 and its closely related M41 and Beaudette strains, respectively. Phylogenetic analyses were performed based on the nucleotide sequence alignment using each ORF from the S to N genes among eight Taiwan and reference strains (Fig. 1) . In the S1 gene, the Taiwan strains showed the highest identity (84%) with the Beaudette reference strain and the lowest (60%) with the DE072. The Taiwan strains (except for 2992/02 and 3374/05) could be clustered into two groups, TW-I and TW-II, based on the S1 gene. Strains 2992/02 and 3374/05 were closely related (>95% identity) to the China genotype VII strains (CK/CH/LDL/97I, CK/CH/LDL/98I, Q1, J2, and T3) (Liu et al., 2006) . In the S2 gene, all of the local strains were grouped with CK/CH/LDL/97I and CK/CH/LDL/98I except for the two viruses isolated before 1995 (strains 1171/92 and 2296/95). Strains 2575/98 and 2992/02 were classified with CK/CH/LDL/97I and CK/CH/LDL/98I in the analyses of the 3a and 3b genes. The E gene analysis segregated the 2992/02 into the China group. Strains 2992/02 and 3468/07 were distributed with the China strains in the M gene analysis. The 5a, 5b and N genes analyses revealed that all of the local strains were in the same group. Nucleotide sequences from the three ORFs were conserved among the Taiwan strains with >92% homology. The 3 6.8 kb genome sequence (S-N gene) of six Taiwan IBVs were queried in the Simplot analysis. Strains 1171/92, H120 and the China strain CK/CH/LDL/97I were used as putative parental strains when strains 2575/98, 2992/02, 3071/03, 3374/05, 3468/07 (TW-I) were queried. The parental strain 1171/92 was replaced with 2296/95 when strain 3263/04 (TW-II) was queried. The similarity plot displays the consecutive nucleotide identity (%) among the queried strain and parental strains. Strains were considered as recombinants if any crossover event took place between two putative parental strains. The breakpoint in which the parental strains have equal identity to the query strain is the predicted recom-bination site. As Fig. 2 demonstrates, crossover events between parental 1171/92 or 2296/95 and CK/CH/LDL/97I were detected in each plot. The recombination sites were located in the S1, S2, 3b, M genes, and the IG region between the M and 5a genes. Each putative recombinant was schematically assembled using Taiwan and China-like sequence fragments. The genomic positions of those crossover sites were indicated within the plot (numbers in red). In Taiwan, field outbreaks are frequently reported despite routine vaccine use. There were originally two IBV genotypes (TW-I and TW-II) circulating in the field Wang and Tsai, 1996) until genetic variants 2992/02 (Huang et al., 2004) and 3374/05 emerged. Both variants showed a high S1 gene homology with the proventriculus pathogenic strains J1, Q2, T3, CK/CH/LDL/97I and CK/CH/LDL/98I, belonging to the China genotype VII (Liu et al., 2006; Yu et al., 2001) . The 3 6.8 kb gene of CK/CH/LDL/97I was used as a putative parental strain. Mass type viruses have been widely used as vaccine strains in Taiwan. Our primary concern was the possible viral recombination resulting from the co-circulation of heterologous vaccine strains in flocks. The consecutive sequence of 3 7.3 kb genome in strain H120 was obtained first in this study and served as another putative parental strain. To our surprise, six local strains experienced crossover events with the strain CK/CH/LDL/97I instead of H120. Thus, inter-strain recombination events had occurred between the IBVs from Taiwan and China. Taiwan and China are geographically separate. Neither live poultry nor processed poultry products from China have been allowed to import into Taiwan for years. How those recombinants arose remains unknown. A recently identified IBV isolate in Korea (Kr/D64/05) also revealed a close relationship to the China CK/CH/LDL/97I-type strains (Lee et al., 2008) . Furthermore, the appearance of the Chinalike strains in Taiwan is reminiscent of the spread of the China QXIBV strain in European countries (Beato et al., 2005; Domanska-Blicharz et al., 2006) . In this case, it could be speculated that migrating birds provide the genetic sources of IBV variants in Taiwan. A chicken-nephropathogenic IBV strain was identified from a non-diseased teal (Anas sp.) in China (Liu et al., 2005) . Thus, the transport of IBVs over long distances by other avian species appears to be possible (Cavanagh, 2005) . However, we still cannot overlook illegal trafficking or unapproved vaccine use in the field. All six recombinants defined in this study were inter-genotypic recombination and the "China-like" sequence substitutions took place in multiple genes. In particular, partial S gene replacement was observed in every recombinant. Since coronaviruses possess different host range or cell tropism through the variance in the S gene (Casais et al., 2003; Kuo et al., 2000) , alterations in the antigenic characteristic in those variants could be expected. In this study, however, nearly all of the strains showed nephropathogenicity (lesions of kidney) in chickens, rather than the proventricular lesions caused by the CK/CH/LDL/97I-type strains. How a variant with a chimeric genome arrangement from heterologous strains demonstrates its tissue tropism or pathogenicity in a host is not clear. In addition, it was found that chicks challenged with virulent CK/CH/LDL/97I were incompletely protected by commercial vaccines and other heterologous strains (Liu et al., 2007) . To effectively control the IB disease in Taiwan, the protective effect of vaccines against challenges from those recombinants awaits to be investigated. To our knowledge, this is the first use of Simplot for genetic analyses of IBV strains. The similarity plot can depict the genetic distance among the aligned sequences in a graphical window. In this study IBV recombinants could be defined directly from the Simplot analyses, and the crossover events and corresponding genome positions were readily observed. In phylogenetic analyses, strains were deduced as recombinants if different genes from the genomes were clustered into separate phylogenetic groups. Parallel results were obtained from both analyses. The emergence of IBV variants through RNA recombination was previously described (Brooks et al., 2004; Jia et al., 1995; Lee and Jackwood, 2000; Mondal and Cardona, 2007; Wang et al., 1993) . Recombination events occurred in multiple genes. The consensus IG sequences CT(T/G)AACAA or the conserved regions around were assumed as the recombination "hot spots" in IBVs (Lee and Jackwood, 2000) . In addition, the CTTTTG sequence was observed adjacent to the putative recombinant junction (Brooks et al., 2004; Wang et al., 1993) . In this study, T-rich motifs (A/T/G)TTTTG, consensus among a recombinant and its parental strains, were located upstream of several crossover sites. For instance, the motifs were found in the S2 gene of the 3263/04, the 3b gene of the 2575/98, and the IG region of the 2992/02 and 3468/07. All of the recombinants experienced crossover events in the 3 end of the S gene, indicating a possible region for template switching in viral RNA synthesis. Similar results were obtained from the IBV strain CU-T2 and Cal99, whose recombination sites were deduced in the 3 751 and 700 nucleotide of the S gene, respectively (Jia et al., 1995; Mondal and Cardona, 2007) . In addition to viral recombination, point mutation is a way of virus evolution. The single base nonsense mutations resulted in truncations of the S and 3b genes in the Taiwan strains 3071/03 and 3374/05. The naturally occurring IBV strains with mutations or complete deletions in non-structural ORFs 3a, 3b, 5a, and 5b revealed lower virus growth rate or titer in ovo (Liu et al., 2008) and in vivo (Mardani et al., 2008) . This study describes the dual viral genomes incorporation in the separated areas having no known interactions between poultry. The molecular investigation in the 3 7.3 kb of the IBV genome demonstrates that inter-strain recombination events contribute to the genetic diversity of the Taiwan strains. Field IBVs in Taiwan undergo genetic recombination and evolution, which might lead to disease control difficulty. The biological and antigenic characteristics of the IBV variants await further studies. Evidence of circulation of a Chinese strain of infectious bronchitis virus (QXIBV) in Italy Molecular epizootiology of avian infectious bronchitis in Russia Phylogenetic analysis of partial S1 and N gene sequences of infectious bronchitis virus isolates from Italy revealed genetic diversity and recombination Genes 3 and 5 of infectious bronchitis virus are accessory protein genes Comparisons of envelope through 5B sequences of infectious bronchitis coronaviruses indicates recombination occurs in the envelope and membrane genes Infectious bronchitis virus S2 gene sequence variability may affect S1 subunit specific antibody binding Co-circulation of four types of infectious bronchitis virus (793/B, 624/I, B1648 and Massachusetts) Gene 5 of the avian coronavirus infectious bronchitis virus is not essential for replication Recombinant avian infectious bronchitis virus expressing a heterologous spike gene demonstrates that the spike protein is a determinant of cell tropism Coronavirus IBV: structural characterization of the spike protein Coronaviruses in poultry and other birds Coronavirus avian infectious bronchitis virus Coronavirus IBV: removal of spike glycopolypeptide S1 by urea abolishes infectivity and haemagglutination but not attachment to cells Relationship between sequence variation in the S1 spike protein of infectious bronchitis virus and the extent of crossprotection in vivo Infectious bronchitis Antigenic and molecular characterization of isolates of the Italy 02 infectious bronchitis virus genotype Molecular epidemiology and evolution of avian infectious bronchitis virus in Spain over a fourteen-year period New variant of IBV in Poland Neither the RNA nor the proteins of open reading frames 3a and 3b of the coronavirus infectious bronchitis virus are essential for replication Recombinant infectious bronchitis coronavirus Beaudette with the spike protein gene of the pathogenic M41 strain remains attenuated but induces protective immunity S1 and N gene analysis of avian infectious bronchitis viruses in Taiwan Sequence changes of infectious bronchitis virus isolates in the 3 7.3 kb of the genome after attenuating passage in embryonated eggs Identification of previously unknown antigenic epitopes on the S and N proteins of avian infectious bronchitis virus Molecular and serologic characterization, pathogenicity, and protection studies with infectious bronchitis virus field isolates from California A novel variant of avian infectious bronchitis virus resulting from recombination among three different strains Antigenic domains on the peplomer protein of avian infectious bronchitis virus: correlation with biological functions Retargeting of coronavirus by substitution of the spike glycoprotein ectodoma: crossing the host cell species barrier Sequence evidence for RNA recombination in field isolates of avian coronavirus infectious bronchitis virus RNA recombination in animal and plant viruses Coronaviridae: the virus and their replication Evidence of genetic diversity generated by recombination among avian coronavirus IBV Genetic diversity of avian infectious bronchitis virus isolates in Korea between Detection of infectious bronchitis virus by multiplex polymerase chain reaction and sequence analysis Isolation of avian infectious bronchitis coronavirus from domestic peafowl (Pavo cristatus) and teal (Anas) Identification of the avian infectious bronchitis coronaviruses with mutations in gene 3 Evaluation of the protection conferred by commercial vaccines and attenuated heterologous isolates in China against the CK/CH/LDL/97I strain of infectious bronchitis coronavirus Genetic diversity of avian infectious bronchitis coronavirus strains isolated in China between Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination Infectious bronchitis viruses with a novel genomic organization Genotypic and phenotypic characterization of the California 99 (Cal99) variant of infectious bronchitis virus Emergence of a coronavirus infectious bronchitis virus mutant with a truncated 3b gene: functional characterization of the 3b protein in pathogenesis and replication MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0 Genetic grouping for the isolates of avian infectious bronchitis virus in Taiwan Evidence of natural recombination within the S1 gene of infectious bronchitis virus In vitro assembled, recombinant infectious bronchitis viruses demonstrate that the 5a open reading frame is not essential for replication Characterization of three infectious bronchitis virus isolates from China associated with proventriculus in vaccinated chickens The financial support of the National Science Council and Council of Agriculture, Executive Yuan, Taiwan, is greatly appreciated.