key: cord-296167-np0b9a7o authors: Mardani, Karim; Noormohammadi, Amir H.; Ignjatovic, Jagoda; Browning, Glenn F. title: Naturally occurring recombination between distant strains of infectious bronchitis virus date: 2010-06-24 journal: Arch Virol DOI: 10.1007/s00705-010-0731-z sha: doc_id: 296167 cord_uid: np0b9a7o New variants of infectious bronchitis virus (IBV) have emerged in Australia despite its geographical isolation and intensive vaccination programs. In the present study, the 3′ terminal 7.2 kb of the genome of a recently isolated variant of IBV (N1/03) was sequenced and compared with the sequences of classical and novel strains of IBV, the two main groups of these viruses in Australia. The comparison revealed that recombination between classical and novel IBVs was responsible for the emergence of the new variant. It was concluded that novel IBVs, which have not been detected since 1993, and which are phylogenically more distant from classical IBVs than turkey coronaviruses, might still be circulating and contributing to the evolution of IBV in Australia. Infectious bronchitis virus (IBV), the prototype of the family Coronaviridae, is an important pathogen in chickens and infects the respiratory tract, kidneys and oviduct, causing reduced performance, reduced egg quality and quantity, increased susceptibility to infection with other pathogens and increased mortality [4] . The genome of IBV is single-stranded, positive-sense RNA of about 27.6 kb, encoding four structural proteins, including the spike glycoprotein (S), the membrane glycoprotein (M), the phosphorylated nucleocapsid protein (N) and the small membrane protein (E). New serotypes and genotypes of IBV emerge frequently in different parts of the world [1, 6, 9, 10, 15, 18, 30, 31] . A number of factors, including mutation and recombination, and the widespread use of live attenuated vaccines, play an important role in increasing the number of new genetic variants [3, 8, 15, 16, 28] . Recombination and mutation are two major forces in coronavirus evolution, and polymerase jumping during coinfection may promote recombination events between coronaviruses [22] . Molecular characterization using sequencing and phylogenetic analysis are effective methods for detection and characterization of recombination events among RNA viruses [29] . N1/62, a nephropathogenic strain of IBV, was the first isolate of IBV obtained in Australia, in 1962 [7] . Since 1962, a number of distinct strains of IBV have been isolated, including a group of strains that were found to have no epitopes in common with other IBV strains on either their N or M proteins [12] . In a study on the pathogenicity of Australian strains of IBV, a change was detected in the prevalence of IBV strains, from highly nephropathogenic strains in the 1960s and 1970s to strains predominantly pathogenic for the respiratory tract in the 1980s and early 1990s, indicating that IBV strains in Australia have undergone a significant change since 1962 [13] . IBV strains in Australia have been classified into two groups, classical and novel [21] , based on their genotypes. The novel IBVs are phylogenetically distant from the classical IBVs, and in fact, classical IBVs and turkey coronaviruses are more similar to each other than they are to the novel IBVs. Recently, the isolation of a new variant of IBV in Australia, from broiler farms located in New South Wales, that had low S1 gene identity to those of classical and novel IBVs but N gene and 3 0 untranslated region sequences similar to those of classical IBVs, suggested the emergence of a new subgroup of IBVs in Australia [14] . In this study, the new variant of IBV that circulated in NSW during 2002 and 2003 was further analysed by comparing the 3 0 terminal 7.2 kb of its genome (all the structural protein genes) with those of the other two IBV groups in Australia. The IBV strain N1/03 used in this study has been described previously [14] . Virus propagated in allantoic fluid was used for RNA extraction. Extraction of RNA was performed using an RNeasy kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. Approximately 50 ll of allantoic fluid was used for each extraction, and RNA was eluted in 30 ll of elution buffer. The extracted RNA was used as template in a reverse transcription reaction. For synthesis of cDNA, 5 ll of extracted RNA was denatured at 100°C for 1 min, cooled by placing on ice for 5 min, and then mixed with 20 ll of reaction buffer containing 10 ll of diethylpyrocarbonate-treated water, 0.5 lM oligo (dT) (Promega, Madison, WI, USA), 0.65 U of RNAguard (Amersham Pharmacia Biotech, Sydney, Australia), 50 lM each of dATP, dTTP, dGTP and dCTP, 5 ll of 59 reaction buffer (Promega), and 50 U of Moloney murine leukaemia virus reverse transcriptase (Promega). The reaction mixture was incubated at 42°C for 1 h and subsequently incubated at 100°C for 5 min to inactivate the reverse transcriptase. The resultant cDNA was immediately used or stored at -70°C for later use. For the amplification reaction, two primers, POLY-F1 (GATTGTGCATGGTGGACAATG) and UTR-R1 (CTGT ACCCTCGATCGTACTC), binding to the 3 0 end of the polymerase gene (nucleotides 20,070-20,090, Beaudette strain, GenBank accession number NC_001451) and the 3 0 untranslated region (nucleotides 27,489-27,508, Beaudette strain) of the IBV genome, respectively, were used to amplify a 7.2-kb fragment of the IBV genome that contains all of the genes for the structural proteins. The PCR reaction was carried out in 50-l mixtures containing 50 lM each of dATP, dTTP, dGTP and dCTP, 0.5 lM of each primer, 5 ll of 10 9 High Fidelity PCR buffer (Invitrogen, Carlsbad, CA, USA), 2 mM magnesium sulfate, 1.5 U Platinum Taq High Fidelity DNA polymerase (Invitrogen) and 5 ll of cDNA as template. Amplification was performed using 35 cycles of incubation at 94°C for 30 s, 57°C for 30 s and 68°C for 8 min, with a final extension at 68°C for 10 min. The resultant PCR products were separated in a 0.8% agarose gel, and the gel was photographed. PCR product was purified using a PCR purification kit (Ultra Clean PCR Clean-Up, Mo Bio Laboratories, Solana Beach, CA, USA) according to the manufacturer's instructions. The 3 0 7.2-kb PCR product from IBV strain N1/03 was cloned in the pGEM-T plasmid using a pGEM-T Vector System Cloning Kit (Promega). The cloned PCR product was sequenced using a Big Dye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems), and the reaction products were sent to the Australian Genomic Research Facility, Walter and Eliza Hall Institute of Medical Research, Melbourne, for analysis on an ABI Prism 3100 Genetic Analyzer. Sequences of different genes of N1/03 were aligned with the same genes from other Australian IBVs and the Beaudette strain (Table 1) using ClustalW [11] . Evolutionary relationships between Australian IBV genes were inferred using the neighbour-joining method [24] . The evolutionary distances were computed using the Maximum Composite Likelihood method [26] and measured as the number of base substitutions per site. Phylogenetic analyses were conducted in MEGA 4 [27] . The 3 0 -terminal 7.2 kb of the genomes of the representative classical and novel Australian IBV strains were aligned with the sequences from the same region of the N1/03 isolate, and the multiple alignment results were introduced into SimPlot version 3.5.1 to identify likely recombination sites [19] . Bootscanning analysis was performed in Simplot to identify and map the putative recombinant IBV genome [25] . This program was also used to find phylogenetically informative sites, as described previously by Robertson et al. [23] . Statistical analysis of the distribution of phylogenetically informative sites was performed using Minitab version 15 (Minitab Inc, US). A two-proportion comparison was performed, and a P value of less than 0.05 was considered significant. Amplification, sequencing and analysis of the 3 0 7.2-kb region of N1/03 The 3 0 7.2 kb of the N1/03 genome was amplified, cloned and sequenced successfully and compared with the same region of a number of classical and novel strains. Phylogenetic trees were constructed for each gene. Analysis of the S gene sequences of Australian IBV strains showed that, among classical strains, they differed by 3-12.5% and, among novel strains, by 3.2-16.4%. The classical and novel IBVs differed by 36.9-39.9%, indicating a major difference between the S genes of these two groups of IBVs. Interestingly, the S1 gene sequence of strain N1/03 differed from those of both classical and novel IBVs by 29.2-32%. Phylogenetic analysis of the S gene grouped N1/03 in a separate cluster from classical and novel IBVs (Fig. 1) . E gene sequences of classical strains differed by 0-9.8%, and those of novel strains by 24-27%. The difference between classical and novel IBVs ranged from 66.5 to 78.8%. More interestingly, the E gene sequence of N1/03 was very similar to those of the classical strains (differences of 3.8-8.9%) and very different from those of novel IBVs (differences of 72.2 to 73.4%). In the E gene phylogenetic tree, N1/03 clustered with the classical IBVs (Fig. 2) . M gene sequences of classical strains differed by 0-10.4%, and those of novel IBVs differed by 3-3.6%. The M gene sequences of classical and novel IBVs differed by 54.2-60.9%. The N1/03 M gene clustered with the classical strains (differences of 2.4-9.8%) (Fig. 3) . The N gene of classical IBVs differed by 0.2-8.3%, and those of novel IBVs differed by 9.9-16.4%. The N gene sequences of classical and novel IBVs differed by 50.7-54.4%. The N gene of N1/03 was more similar to the N genes of the classical IBVs (differences of 2.4-7.7%) than of novel IBVs (differences of 54.3-55.6%) and clustered with the classical strains in phylogenetic analyses (Fig. 4) . To identify the region likely to have been involved in recombination in N1/03, similarity plots and bootscan analysis were performed using strains VicS and N1/88 as representatives of the two main groups of IBV in Australia and the Armidale strain as a query. In both the similarity plots and the bootscan graph, N1/03 had greater similarity to N1/88 (a novel IBV strain) in its S gene region (almost 3 kb), while the reminder of the 3 0 -terminal 7.2-kb region had greater similarity to VicS (a classical IBV strain) (Fig. 5, 6 ). Data from phylogenetic analyses also suggested that N1/03 had sequence similarity to both main groups of Australian IBVs. These findings suggested that there had been a recombination event with a crossover point at the end of the S gene (base number 3000) of N1/03 (Fig. 5, 6) . To obtain a precise picture of this possible crossover point, the distribution of phylogenetically informative sites along the 3 0 -terminal 7.2-kb sequences of N1/88 and VicS were examined. Using Simplot, 456 informative sites were found that could support one of three possible trees. For almost 90% of the S gene (3 kb), 365 informative sites were found, and 84% of these informative sites supported a tree clustering N1/03 with the novel IBV strain N1/88. In contrast, only 8% of the informative sites supported each of the other two possible trees in this region (P \ 0.0005). For the rest of sequences (3,000-7,200 bp), 95 informative sites were found, with 60% supporting the tree clustering N1/03 with VicS and only 14% supporting the tree clustering N1/03 with N1/88 (P \ 0.05). This indicated that N1/03 clustered with the classical IBV strain VicS throughout the 3 0 -terminal 4.2 kb of its genome and confirmed that a recombination event had occurred at the end of the S gene region (at around position 3000) and that N1/03 is a recombinant virus that emerged from recombination between these two groups of IBVs in Australia. New IBV isolates have been reported in recent years in Australia and have been classified as subgroup 3 viruses. The emergence of these new isolates was considered to be unexplained, as was the earlier emergence of subgroup 2, or novel viruses. N1/03 is one of these new viruses, which have been isolated from broiler farms using one of the four commercial IBV vaccines [14] . This study was undertaken to explore the mechanism behind the emergence of the subgroup 3 IBVs in Australia. Based on phylogenetic analysis of the complete 3 0 7.2-kb region of the genome, N1/03 was distinct from both classical and novel IBVs. However, the S gene sequence of Fig. 4 Evolutionary relationships between N genes of Australian IBV strains inferred using the neighbour-joining method. The evolutionary distances were computed using the Maximum Composite Likelihood method and are given in the number of base substitutions per site N1/03 was equally different from both the classical and novel IBVs, while the rest of the 3 0 -terminal 7.2-kb region of the genome clustered closely with those of the classical IBVs and was very different from those of the novel IBVs. The similarity of the N1/03 E, M and N genes to those of classical IBVs suggested that this isolate was closely related to the most commonly used Australian IBV vaccine strains. As these strains have been used for a long time in Australia, it is possible that they have contributed to the emergence of variant IBVs in the field. The S gene of N1/03 clustered more closely with the novel IBVs, and similarity plotting and bootscan analysis confirmed the close relationship between the S genes of N1/03 and the novel IBVs and suggested that the crossover event occurred near the 3 0 end of this gene. These analyses provide convincing evidence that recombination has occurred between classical and novel IBVs, leading to emergence of variant IBVs in the field. Novel IBVs have not been isolated since 1992, but this finding suggests that novel IBVs are still circulating in the field in Australia. The failure to detect them since 1992 may be a result of the relatively slow growth of these viruses [20] . These results confirm the hypothesis of Ignjatovic et al. [14] that recombination may have played a role in the emergence of this new variant of IBV in Australia. Recombination between IBVs and its role in emergence of new IBV variants has been reported previously [2, 15, 17, 21] and may occur at multiple sites [15, 17] , although crossover events occur more frequently at the 3 0 end of the S gene [5, 15, 21] , as seen in this study. However, it is particularly notable here because of the very significant phylogenic distance between the novel and classical Australian strains of IBV. Indeed, the novel strains of IBV are more distinct from the classical strains than is turkey coronavirus, and they lack several of the smaller nonstructural genes found in classical strains. In conclusion, this study showed that the recombination was involved in the emergence of the new variant of IBV in Australia and that a rapid and reliable technique is needed to determine whether the novel IBVs are still circulating in poultry farms. It would be appropriate to sequence and analyse the polymerase genes of classical, novel and new variant strains of IBV to obtain further information about the relationships between the different Australian IBVs. Isolation and characterization of new infectious bronchitis virus variants in Hungary Comparisons of envelope through 5B sequences of infectious bronchitis coronaviruses indicates recombination occurs in the envelope and membrane genes Recent advances in avian virology Severe acute respiratory syndrome vaccine development: experiences of vaccination against avian infectious bronchitis coronavirus Identification of Taiwan and China-like recombinant avian infectious bronchitis viruses in Taiwan Isolation of a new serotype of infectious bronchitis-like virus from chickens in England Infectious avian nephrosis (uraemia) in Australia A recombination event, induced in ovo, between a low passage infectious bronchitis virus field isolate and a highly embryo adaptedvaccine strain Variant serotypes of infectious bronchitis virus isolated from commercial layer and broiler chickens A 'new' strain of infectious bronchitis virus infecting domestic fowl in Great Britain BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT A long-term study of Australian infectious bronchitis viruses indicates a major antigenic change in recently isolated strains Pathogenicity of Australian strains of avian infectious bronchitis virus Isolation of a variant infectious bronchitis virus in Australia that further illustrates diversity among emerging strains A novel variant of avian infectious bronchitis virus resulting from recombination among three different strains Sequence evidence for RNA recombination in field isolates of avian coronavirus infectious bronchitis virus Evidence of genetic diversity generated by recombination among avian coronavirus IBV A new genotype of nephropathogenic infectious bronchitis virus circulating in vaccinated and nonvaccinated flocks in China Fulllength 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 Coronavirus replication and pathogenesis: Implications for the recent outbreak of severe acute respiratory syndrome (SARS), and the challenge for vaccine development Recombination in AIDS viruses The neighbor-joining method: a new method for reconstructing phylogenetic trees Identification of breakpoints in intergenotypic recombinants of HIV type 1 by bootscanning Prospects for inferring very large phylogenies by using the neighbor-joining method MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0 Evidence of natural recombination within the S1 gene of infectious bronchitis virus Evolutionary aspects of recombination in RNA viruses Characterization of three infectious bronchitis virus isolates from China associated with proventriculus in vaccinated chickens Avian infectious bronchitis: characterization of new isolates from Italy Acknowledgments The senior author would like to thank the Department of Education, Employment and Workplace Relations of Australia for providing him with an Endeavour Research Fellowship Award, which enabled this work to be completed. He also thanks Urmia University, Iran, for granting him leave to go to Australia to undertake this work.